Assessing the effect of prevalence on the predictive performance of species distribution models using simulated data

Aim The proportion of sampled sites where a species is present is known as prevalence. Empirical studies have shown that prevalence can affect the predictive performance of species distribution models. This paper uses simulated species data to examine how prevalence and the form of species environmental dependence affect the assessment of the predictive performance of models. Methods Simulated species data were based on various functions of simulated environmental data with differing degrees of spatial correlation. Seven model performance measures – sensitivity, specificity, class‐average (CA), overall prediction success, kappa (κ), normalized mutual information (NMI) and area under the receiver operating characteristic curve (AUC) – were applied to species models fitted by three regression methods. The response of the performance measures to prevalence was then assessed. Three probability threshold selection methods used to convert fitted logistic model values to presence or absence were also assessed. Results The study shows that the extent to which prevalence affects model performance depends on the modelling technique and its degree of success in capturing dominant environmental determinants. It also depends on the statistic used to measure model performance and the probability threshold method. The response based on κ generally preferred models with medium prevalence. All performance measures were least affected by prevalence when the probability threshold was chosen to maximize predictive performance or was based directly on prevalence. In these cases, the responses based on AUC, CA and NMI generally preferred models with small or large prevalence. Main conclusions The effect of prevalence on the predictive performance of species distribution models has a methodological basis. Relevant factors include the success of the fitted distribution model in capturing the dominant environmental determinant, the model performance measure and the probability threshold selection method. The fixed probability threshold method yields a marked response of model performance to prevalence and is therefore not recommended. The study explains previous empirical results obtained with real data.     

Codon optimization for high level expression of human bone morphogenetic protein-2 in Escherichia coli

Codons in the open reading frame (ORF) encoding for human bone morphogenetic protein-2 (hBMP-2) were optimized to reach high level expression in Escherichia coli. The optimization was done by the computer programs DNA works and DNA Star according to Thermodynamically Balanced Inside Out (TBIO) approach. The ORF consisting of 342 base pairs (bp) was assembled using two-steps Polymerase Chain Reaction, cloned into a pGEM-T vector with a mutation rate of 6.38 bp per kb and transformed into E. coli JM109. After a DNA sequence confirmation, mutation-free ORF was subcloned into pET32b and transformed into E. coli BL21(DE3). The rhBMP-2 was produced as a thioredoxin-his-tag fusion protein at relatively high level, approximately 60% of total intracellular proteins as inclusion bodies (IB), with a yield of 1.39 g per liter culture. Solubilization of IB gave soluble monomer rhBMP-2 with a recovery of 13.6% and refolding of soluble rhBMP-2 produced dimeric forms with a yield of 8.7%. The size and identity of the purified rhBMP-2 was confirmed by nano-LC-MS/MS2 analysis. Our work demonstrates for the first time that by using TBIO approach, a codon-optimized ORF encoding for rhBMP-2 protein can be expressed at high level in E. coli expression system.     

In silico prioritisation of candidate genes for prokaryotic gene function discovery: an application of phylogenetic profiles

Background  

In silico candidate gene prioritisation (CGP) aids the discovery of gene functions by ranking genes according to an objective relevance score. While several CGP methods have been described for identifying human disease genes, corresponding methods for prokaryotic gene function discovery are lacking. Here we present two prokaryotic CGP methods, based on phylogenetic profiles, to assist with this task.     

Improved reduced representation bisulfite sequencing for epigenomic profiling of clinical samples

Background

DNA methylation plays crucial roles in epigenetic gene regulation in normal development and disease pathogenesis. Efficient and accurate quantification of DNA methylation at single base resolution can greatly advance the knowledge of disease mechanisms and be used to identify potential biomarkers. We developed an improved pipeline based on reduced representation bisulfite sequencing (RRBS) for cost-effective genome-wide quantification of DNA methylation at single base resolution. A selection of two restriction enzymes (Taq^αI and MspI) enables a more unbiased coverage of genomic regions of different CpG densities. We further developed a highly automated software package to analyze bisulfite sequencing results from the Solexa GAIIx system.

Results

With two sequencing lanes, we were able to quantify ~1.8 million individual CpG sites at a minimum sequencing depth of 10. Overall, about 76.7% of CpG islands, 54.9% of CpG island shores and 52.2% of core promoters in the human genome were covered with at least 3 CpG sites per region.

Conclusions

With this new pipeline, it is now possible to perform whole-genome DNA methylation analysis at single base resolution for a large number of samples for understanding how DNA methylation and its changes are involved in development, differentiation, and disease pathogenesis.     

Oncogenic comparison of human papillomavirus type 58 E7 variants

Human papillomavirus 58 (HPV58) ranks the second or third in East Asian cervical cancers. Current studies on HPV58 are scarce and focus on the prototype. Previously, we identified the three most common circulating HPV58 E7 strains contained amino acid alterations: G41R/G63D (51%), T20I/G63S (22%) and T74A/D76E (14%) respectively. Among them, the T20I/G63S variant (V1) had a stronger epidemiological association with cervical cancer. We therefore suggested that V1 possessed stronger oncogenicity than the other two variants. Here, we performed phenotypic assays to characterize and compare their oncogenicities with HPV58 E7 prototype. Our results showed that overexpression of V1 conferred a higher colony‐forming ability to primary murine epithelial cells than prototype (P < 0.05) and other variants, implicating its higher immortalising potential. Further experiments showed that both V1 and prototype enhanced the anchorage‐independent growth of NIH/3T3 cells (P < 0.001), implicating their stronger transforming power than the two other variants. Moreover, they possessed an increased ability to degrade pRb (P < 0.001), which is a major effector pathway of E7‐driven oncogenesis. Our work represents the first study to compare the oncogenicities of HPV58 E7 prototype and variants. These findings deepened our understanding of HPV58 and might inform clinical screening and follow‐up strategy.     

Structural basis for the specificity of PPM1H phosphatase for Rab GTPases

LRRK2 serine/threonine kinase is associated with inherited Parkinson’s disease. LRRK2 phosphorylates a subset of Rab GTPases within their switch 2 motif to control their interactions with effectors. Recent work has shown that the metal‐dependent protein phosphatase PPM1H counteracts LRRK2 by dephosphorylating Rabs. PPM1H is highly selective for LRRK2 phosphorylated Rabs, and closely related PPM1J exhibits no activity towards substrates such as Rab8a phosphorylated at Thr72 (pThr72). Here, we have identified the molecular determinant of PPM1H specificity for Rabs. The crystal structure of PPM1H reveals a structurally conserved phosphatase fold that strikingly has evolved a 110‐residue flap domain adjacent to the active site. The flap domain distantly resembles tudor domains that interact with histones in the context of epigenetics. Cellular assays, crosslinking and 3‐D modelling suggest that the flap domain encodes the docking motif for phosphorylated Rabs. Consistent with this hypothesis, a PPM1J chimaera with the PPM1H flap domain dephosphorylates pThr72 of Rab8a both in vitro and in cellular assays. Therefore, PPM1H has acquired a Rab‐specific interaction domain within a conserved phosphatase fold.