A Statistical Framework for Improving Genomic Annotations of Prokaryotic Essential Genes

Large-scale systematic analysis of gene essentiality is an important step closer toward unraveling the complex relationship between genotypes and phenotypes. Such analysis cannot be accomplished without unbiased and accurate annotations of essential genes. In current genomic databases, most of the essential gene annotations are derived from whole-genome transposon mutagenesis (TM), the most frequently used experimental approach for determining essential genes in microorganisms under defined conditions. However, there are substantial systematic biases associated with TM experiments. In this study, we developed a novel Poisson model–based statistical framework to simulate the TM insertion process and subsequently correct the experimental biases. We first quantitatively assessed the effects of major factors that potentially influence the accuracy of TM and subsequently incorporated relevant factors into the framework. Through iteratively optimizing parameters, we inferred the actual insertion events occurred and described each gene’s essentiality on probability measure. Evaluated by the definite mapping of essential gene profile in Escherichia coli, our model significantly improved the accuracy of original TM datasets, resulting in more accurate annotations of essential genes. Our method also showed encouraging results in improving subsaturation level TM datasets. To test our model’s broad applicability to other bacteria, we applied it to Pseudomonas aeruginosa PAO1 and Francisella tularensis novicida TM datasets. We validated our predictions by literature as well as allelic exchange experiments in PAO1. Our model was correct on six of the seven tested genes. Remarkably, among all three cases that our predictions contradicted the TM assignments, experimental validations supported our predictions. In summary, our method will be a promising tool in improving genomic annotations of essential genes and enabling large-scale explorations of gene essentiality. Our contribution is timely considering the rapidly increasing essential gene sets. A Webserver has been set up to provide convenient access to this tool. All results and source codes are available for download upon publication at http://research.cchmc.org/essentialgene/.     

Methionine Synthase A2756G Polymorphism and Risk of Colorectal Adenoma and Cancer: Evidence Based on 27 Studies

Methionine synthase (MTR), which plays a central role in maintaining adequate intracellular folate, methionine and normal homocysteine concentrations, was thought to be involved in the development of colorectal cancer (CRC) and colorectal adenoma (CRA) by affecting DNA methylation. However, studies on the association between MTR A2756G polymorphism and CRC/CRA remain conflicting. We conducted a meta-analysis of 27 studies, including 13465 cases and 20430 controls for CRC, and 4844 cases and 11743 controls for CRA. Potential sources of heterogeneity and publication bias were also systematically explored. Overall, the summary odds ratio of G variant for CRC was 1.03 (95% CI: 0.96–1.09) and 1.05 (95% CI: 0.99–1.12) for CRA. No significant results were observed in heterozygous and homozygous when compared with wild genotype for these polymorphisms. In the stratified analyses according to ethnicity, source of controls, sample size, sex, and tumor site, no evidence of any gene-disease association was obtained. Results from the meta-analysis of four studies on MTR stratified according to smoking and alcohol drinking status showed an increased CRC risk in heavy smokers (OR = 2.06, 95% CI: 1.32–3.20) and heavy drinkers (OR = 2.00, 95% CI: 1.28–3.09) for G allele carriers. This meta-analysis suggests that the MTR A2756G polymorphism is not associated with CRC/CRA susceptibility and that gene-environment interaction may exist.     

New Insights into the Phylogeny and Worldwide Dispersion of Two Closely Related Nematode Species,Bursaphelenchus xylophilus and Bursaphelenchus mucronatus

The pinewood nematode, Bursaphelenchus xylophilus, is one of the greatest threats to coniferous forests worldwide, causing severe ecological damage and economic loss. The biology of B. xylophilus is similar to that of its closest relative, B. mucronatus, as both species share food resources and insect vectors, and have very similar morphological characteristics, although little pathogenicity to conifers has been associated with B. mucronatus. Using both nuclear and mitochondrial DNA markers, we show that B. xylophilus and B. mucronatus form distinct phylogenetic groups with contrasting phylogeographic patterns. B. xylophilus presents lower levels of intraspecific diversity than B. mucronatus, as expected for a species that evolved relatively recently through geographical or reproductive isolation. Genetic diversity was particularly low in recently colonised areas, such as in southwestern Europe. By contrast, B. mucronatus displays high levels of genetic diversity and two well-differentiated clades in both mitochondrial and nuclear DNA phylogenies. The lack of correlation between genetic and geographic distances in B. mucronatus suggests intense gene flow among distant regions, a phenomenon that may have remained unnoticed due to the reduced pathogenicity of the species. Overall, our findings suggest that B. xylophilus and B. mucronatus have different demographic histories despite their morphological resemblance and ecological overlap. These results suggest that Bursaphelenchus species are a valuable model for understanding the dispersion of invasive species and the risks posed to native biodiversity and ecosystems.     

Evaluating Purifying Selection in the Mitochondrial DNA of Various Mammalian Species

Mitochondrial DNA (mtDNA), the circular DNA molecule inside the mitochondria of all eukaryotic cells, has been shown to be under the effect of purifying selection in several species. Traditional testing of purifying selection has been based simply on ratios of nonsynonymous to synonymous mutations, without considering the relative age of each mutation, which can be determined by phylogenetic analysis of this non-recombining molecule. The incorporation of a mutation time-ordering from phylogeny and of predicted pathogenicity scores for nonsynonymous mutations allow a quantitative evaluation of the effects of purifying selection in human mtDNA. Here, by using this additional information, we show that purifying selection undoubtedly acts upon the mtDNA of other mammalian species/genera, namely Bos sp., Canis lupus, Mus musculus, Orcinus orca, Pan sp. and Sus scrofa. The effects of purifying selection were comparable in all species, leading to a significant major proportion of nonsynonymous variants with higher pathogenicity scores in the younger branches of the tree. We also derive recalibrated mutation rates for age estimates of ancestors of these various species and proposed a correction curve in order to take into account the effects of selection. Understanding this selection is fundamental to evolutionary studies and to the identification of deleterious mutations.     

BmPAH Catalyzes the Initial Melanin Biosynthetic Step in Bombyx mori

Pigmentation during insect development is a primal adaptive requirement. In the silkworm, melanin is the primary component of larval pigments. The rate limiting substrate in melanin synthesis is tyrosine, which is converted from phenylalanine by the rate-limiting enzyme phenylalanine hydroxylase (PAH). While the role of tyrosine, derived from phenylalanine, in the synthesis of fiber proteins has long been known, the role of PAH in melanin synthesis is still unknown in silkworm. To define the importance of PAH, we cloned the cDNA sequence of BmPAH and expressed its complete coding sequence using the Bac-to-Bac baculovirus expression system. Purified recombinant protein had high PAH activity, some tryptophan hydroxylase activity, but no tyrosine hydroxylase activity, which are typical properties of PAH in invertebrates. Because melanin synthesis is most robust during the embryonic stage and larval integument recoloring stage, we injected BmPAH dsRNA into silkworm eggs and observed that decreasing BmPAH mRNA reduced neonatal larval tyrosine and caused insect coloration to fail. In vitro cultures and injection of 4^th instar larval integuments with PAH inhibitor revealed that PAH activity was essential for larval marking coloration. These data show that BmPAH is necessary for melanin synthesis and we propose that conversion of phenylalanine to tyrosine by PAH is the first step in the melanin biosynthetic pathway in the silkworm.     

Genome-Wide Association Studies for Growth and Meat Production Traits in Sheep

Background

Growth and meat production traits are significant economic traits in sheep. The aim of the study is to identify candidate genes affecting growth and meat production traits at genome level with high throughput single nucleotide polymorphisms (SNP) genotyping technologies.

Methodology and Results

Using Illumina OvineSNP50 BeadChip, we performed a GWA study in 329 purebred sheep for 11 growth and meat production traits (birth weight, weaning weight, 6-month weight, eye muscle area, fat thickness, pre-weaning gain, post-weaning gain, daily weight gain, height at withers, chest girth, and shin circumference). After quality control, 319 sheep and 48,198 SNPs were analyzed by TASSEL program in a mixed linear model (MLM). 36 significant SNPs were identified for 7 traits, and 10 of them reached genome-wise significance level for post-weaning gain. Gene annotation was implemented with the latest sheep genome Ovis_aries_v3.1 (released October 2012). More than one-third SNPs (14 out of 36) were located within ovine genes, others were located close to ovine genes (878bp-398,165bp apart). The strongest new finding is 5 genes were thought to be the most crucial candidate genes associated with post-weaning gain: s58995.1 was located within the ovine genes MEF2B and RFXANK, OAR3_84073899.1, OAR3_115712045.1 and OAR9_91721507.1 were located within CAMKMT, TRHDE, and RIPK2 respectively. GRM1, POL, MBD5, UBR2, RPL7 and SMC2 were thought to be the important candidate genes affecting post-weaning gain too. Additionally, 25 genes at chromosome-wise significance level were also forecasted to be the promising genes that influencing sheep growth and meat production traits.

Conclusions

The results will contribute to the similar studies and facilitate the potential utilization of genes involved in growth and meat production traits in sheep in future.     

Altered Pattern of Spontaneous Brain Activity in the Patients with End-Stage Renal Disease: A Resting-State Functional MRI Study with Regional Homogeneity Analysis

Purpose

To investigate the pattern of spontaneous neural activity in patients with end-stage renal disease (ESRD) with and without neurocognitive dysfunction using resting-state functional magnetic resonance imaging (rs-fMRI) with a regional homogeneity (ReHo) algorithm.

Materials and Methods

rs-fMRI data were acquired in 36 ESRD patients (minimal nephro-encephalopathy [MNE], n = 19, 13 male, 37±12.07 years; non-nephro-encephalopathy [non-NE], n = 17, 11 male, 38±12.13 years) and 20 healthy controls (13 male, 7 female, 36±10.27 years). Neuropsychological (number connection test type A [NCT-A], digit symbol test [DST]) and laboratory tests were performed in all patients. The Kendall''s coefficient of concordance (KCC) was used to measure the regional homogeneity for each subject. The regional homogeneity maps were compared using ANOVA tests among MNE, non-NE, and healthy control groups and post hoc t -tests between each pair in a voxel-wise way. A multiple regression analysis was performed to evaluate the relationships between ReHo index and NCT-A, DST scores, serum creatinine and urea levels, disease and dialysis duration.

Results

Compared with healthy controls, both MNE and non-NE patients showed decreased ReHo in the multiple areas of bilateral frontal, parietal and temporal lobes. Compared with the non-NE, MNE patients showed decreased ReHo in the right inferior parietal lobe (IPL), medial frontal cortex (MFC) and left precuneus (PCu). The NCT-A scores and serum urea levels of ESRD patients negatively correlated with ReHo values in the frontal and parietal lobes, while DST scores positively correlated with ReHo values in the bilateral PCC/precuneus, MFC and inferior parietal lobe (IPL) (all P<0.05, AlphaSim corrected). No significant correlations were found between any regional ReHo values and disease duration, dialysis duration and serum creatinine values in ESRD patients (all P>0.05, AlphaSim corrected).

Conclusion

Diffused decreased ReHo values were found in both MNE and non-NE patients. The progressively decreased ReHo in the default mode network (DMN), frontal and parietal lobes might be trait-related in MNE. The ReHo analysis may be potentially valuable for elucidating neurocognitive abnormalities of ESRD patients and detecting the development from non-NE to MNE.     

高沥水性钝顶螺旋藻新品系的选育及超微结构与RAPD分析

[目的]选育高沥水性的螺旋藻新品系,显著降低藻粉生产中干燥的能耗。[方法]以用于工厂化培植的钝顶螺旋藻ZJU0115为出发品系,用组织匀浆-离心沉降法制得其原生质球,并先以0.6%EMS处理30 min再用2.4 kGy的~(60)Coγ射线辐照,经含0.02%黄原胶(xanthan gum)的Zarrouk's培养液筛选、藻丝单体分离培养、藻泥持水率和胞外多糖(EPS)等检测及生产培植试验。[结果]获得了一株产量、蛋白质和多糖含量与ZJU0115相当,而藻泥持水率和EPS含量分别下降5.9%和29.7%的突变体,命名为ZJU0115(HD)。超微结构与随机扩增多态性DNA标记(random amplified polymorphic DNA,RAPD)分析结果显示,与其亲本ZJU0115相比,ZJU0115(HD)藻丝表面更光滑,可能为酸性多糖的乳白状粘附物也更少;ZJU0115(HD)细胞内的多磷酸盐颗粒显著变小且呈弥散状;基因组DNA在随机引物S90的扩增产物中显示出多态性差异。[结论]ZJU0115(HD)在工厂化培植中生产性状好、高沥水性能稳定,藻泥的干燥能耗降低了近50%,它的育成与应用,有助于推进当前螺旋藻产业迈向高效、节能、绿色、环保发展的新阶段。     

Untangling taxonomy: a DNA barcode reference library for Canadian spiders

Approximately 1460 species of spiders have been reported from Canada, 3% of the global fauna. This study provides a DNA barcode reference library for 1018 of these species based upon the analysis of more than 30 000 specimens. The sequence results show a clear barcode gap in most cases with a mean intraspecific divergence of 0.78% vs. a minimum nearest‐neighbour (NN) distance averaging 7.85%. The sequences were assigned to 1359 Barcode index numbers (BINs) with 1344 of these BINs composed of specimens belonging to a single currently recognized species. There was a perfect correspondence between BIN membership and a known species in 795 cases, while another 197 species were assigned to two or more BINs (556 in total). A few other species (26) were involved in BIN merges or in a combination of merges and splits. There was only a weak relationship between the number of specimens analysed for a species and its BIN count. However, three species were clear outliers with their specimens being placed in 11–22 BINs. Although all BIN splits need further study to clarify the taxonomic status of the entities involved, DNA barcodes discriminated 98% of the 1018 species. The present survey conservatively revealed 16 species new to science, 52 species new to Canada and major range extensions for 426 species. However, if most BIN splits detected in this study reflect cryptic taxa, the true species count for Canadian spiders could be 30–50% higher than currently recognized.