Oral spirochetes implicated in dental diseases are widespread in normal human subjects and carry extremely diverse integron gene cassettes

The NIH Human Microbiome Project (HMP) has produced several hundred metagenomic data sets, allowing studies of the many functional elements in human-associated microbial communities. Here, we survey the distribution of oral spirochetes implicated in dental diseases in normal human individuals, using recombination sites associated with the chromosomal integron in Treponema genomes, taking advantage of the multiple copies of the integron recombination sites (repeats) in the genomes, and using a targeted assembly approach that we have developed. We find that integron-containing Treponema species are present in ～80% of the normal human subjects included in the HMP. Further, we are able to de novo assemble the integron gene cassettes using our constrained assembly approach, which employs a unique application of the de Bruijn graph assembly information; most of these cassette genes were not assembled in whole-metagenome assemblies and could not be identified by mapping sequencing reads onto the known reference Treponema genomes due to the dynamic nature of integron gene cassettes. Our study significantly enriches the gene pool known to be carried by Treponema chromosomal integrons, totaling 826 (598 97% nonredundant) genes. We characterize the functions of these gene cassettes: many of these genes have unknown functions. The integron gene cassette arrays found in the human microbiome are extraordinarily dynamic, with different microbial communities sharing only a small number of common genes.     

Compositional and genetic alterations in Graves’ disease gut microbiome reveal specific diagnostic biomarkers

Graves’ Disease is the most common organ-specific autoimmune disease and has been linked in small pilot studies to taxonomic markers within the gut microbiome. Important limitations of this work include small sample sizes and low-resolution taxonomic markers. Accordingly, we studied 162 gut microbiomes of mild and severe Graves’ disease (GD) patients and healthy controls. Taxonomic and functional analyses based on metagenome-assembled genomes (MAGs) and MAG-annotated genes, together with predicted metabolic functions and metabolite profiles, revealed a well-defined network of MAGs, genes and clinical indexes separating healthy from GD subjects. A supervised classification model identified a combination of biomarkers including microbial species, MAGs, genes and SNPs, with predictive power superior to models from any single biomarker type (AUC = 0.98). Global, cross-disease multi-cohort analysis of gut microbiomes revealed high specificity of these GD biomarkers, notably discriminating against Parkinson’s Disease, and suggesting that non-invasive stool-based diagnostics will be useful for these diseases.Subject terms: Microbiome, Biomarkers, Population genetics     

<Emphasis Type="Italic">In silico</Emphasis> Identification of Metagenomic Signature Describing Neurometabolic Potential of Normal Human Gut Microbiota

A great amount of attention has been paid to the study of the microbiota–gut–brain axis in recent years. Gut microbiota can affect development and functioning of the brain through synthesis of various neuroactive metabolites, such as neurotransmitters, hormones, and other compounds. In the present study, the presence and distribution are analyzed for the genes controlling the synthesis of enzymes involved in production of neuroactive compounds in 147 gut metagenomes of healthy people from Human Microbiome Project database and synthetic metagenome artificially assembled from 508 bacterial genomes. The analysis is conducted using the collected catalog of orthologs for 17 key enzymes and an algorithm developed for their search. As a result of analyses of genomic and metagenomic data of healthy people, seven bacterial genera containing the greatest number of enzyme genes and 8 enzymes out of 17 that are observed the most frequently are chosen. It is assumed that the selected “core” genera and enzymes form a metagenomic signature reflecting the neurometabolic potential of the human intestinal microbiota in the norm.     

Archaeology and evolution of transfer RNA genes in the Escherichia coli genome

Transfer RNA genes tend to be presented in multiple copies in the genomes of most organisms, from bacteria to eukaryotes. The evolution and genomic structure of tRNA genes has been a somewhat neglected area of molecular evolution. Escherichia coli, the first phylogenetic species for which more than two different strains have been sequenced, provides an invaluable framework to study the evolution of tRNA genes. In this work, a detailed analysis of the tRNA structure of the genomes of Escherichia coli strains K12, CFT073, and O157:H7, Shigella flexneri 2a 301, and Salmonella typhimurium LT2 was carried out. A phylogenetic analysis of these organisms was completed, and an archaeological map depicting the main events in the evolution of tRNA genes was drawn. It is shown that duplications, deletions, and horizontal gene transfers are the main factors driving tRNA evolution in these genomes. On average, 0.64 tRNA insertions/duplications occur every million years (Myr) per genome per lineage, while deletions occur at the slower rate of 0.30 per million years per genome per lineage. This work provides a first genomic glance at the problem of tRNA evolution as a repetitive process, and the relationship of this mechanism to genome evolution and codon usage is discussed.     

High proportions of deleterious polymorphisms in constrained human genes

Previous studies on human mitochondrial genomes showed that the ratio of intra-specific diversities at nonsynonymous-to-synonymous positions was two to ten times higher than the ratio of interspecific divergences at these positions, suggesting an excess of slightly deleterious nonsynonymous polymorphisms. However, such an overabundance of nonsynonymous single nucleotide polymorphisms (SNPs) was not found in human nuclear genomes. Here, genome-wide estimates using >14,000 human-chimp nuclear genes and 1 million SNPs from four human genomes showed a significant proportion of deleterious nonsynonymous SNPs (～ 15%). Importantly, this study reveals a negative correlation between the magnitude of selection pressure and the proportion of deleterious SNPs on human genes. The proportion of deleterious amino acid replacement polymorphisms is 3.5 times higher in genes under high purifying selection compared with that in less constrained genes (28% vs. 8%). These results are explained by differences in the extent of contribution of mildly deleterious mutations to diversity and substitution.     

The complete chloroplast genome of Gracilariopsis lemaneiformis (Rhodophyta) gives new insight into the evolution of family Gracilariaceae

The complete chloroplast genome of Gracilariopsis lemaneiformis was recovered from a Next Generation Sequencing data set. Without quadripartite structure, this chloroplast genome (183,013 bp, 27.40% GC content) contains 202 protein‐coding genes, 34 tRNA genes, 3 rRNA genes, and 1 tmRNA gene. Synteny analysis showed plasmid incorporation regions in chloroplast genomes of three species of family Gracilariaceae and in Grateloupia taiwanensis of family Halymeniaceae. Combined with reported red algal plasmid sequences in nuclear and mitochondrial genomes, we postulated that red algal plasmids may have played an important role in ancient horizontal gene transfer among nuclear, chloroplast, and mitochondrial genomes. Substitution rate analysis showed that purifying selective forces maintaining stability of protein‐coding genes of nine red algal chloroplast genomes over long periods must be strong and that the forces acting on gene groups and single genes of nine red algal chloroplast genomes were similar and consistent. The divergence of Gp. lemaneiformis occurred ~447.98 million years ago (Mya), close to the divergence time of genus Pyropia and Porphyra (443.62 Mya).     

Transposases are the most abundant,most ubiquitous genes in nature

Genes, like organisms, struggle for existence, and the most successful genes persist and widely disseminate in nature. The unbiased determination of the most successful genes requires access to sequence data from a wide range of phylogenetic taxa and ecosystems, which has finally become achievable thanks to the deluge of genomic and metagenomic sequences. Here, we analyzed 10 million protein-encoding genes and gene tags in sequenced bacterial, archaeal, eukaryotic and viral genomes and metagenomes, and our analysis demonstrates that genes encoding transposases are the most prevalent genes in nature. The finding that these genes, classically considered as selfish genes, outnumber essential or housekeeping genes suggests that they offer selective advantage to the genomes and ecosystems they inhabit, a hypothesis in agreement with an emerging body of literature. Their mobile nature not only promotes dissemination of transposable elements within and between genomes but also leads to mutations and rearrangements that can accelerate biological diversification and—consequently—evolution. By securing their own replication and dissemination, transposases guarantee to thrive so long as nucleic acid-based life forms exist.     

Long- and Short-Term Selective Forces on Malaria Parasite Genomes

Plasmodium parasites, the causal agents of malaria, result in more than 1 million deaths annually. Plasmodium are unicellular eukaryotes with small ∼23 Mb genomes encoding ∼5200 protein-coding genes. The protein-coding genes comprise about half of these genomes. Although evolutionary processes have a significant impact on malaria control, the selective pressures within Plasmodium genomes are poorly understood, particularly in the non-protein-coding portion of the genome. We use evolutionary methods to describe selective processes in both the coding and non-coding regions of these genomes. Based on genome alignments of seven Plasmodium species, we show that protein-coding, intergenic and intronic regions are all subject to purifying selection and we identify 670 conserved non-genic elements. We then use genome-wide polymorphism data from P. falciparum to describe short-term selective processes in this species and identify some candidate genes for balancing (diversifying) selection. Our analyses suggest that there are many functional elements in the non-genic regions of these genomes and that adaptive evolution has occurred more frequently in the protein-coding regions of the genome.     

33种鱼类微型反向重复转座元件鉴定

微型反向重复转座元件(MITEs)是一类短的非自主DNA转座子, 其分布的位置会对宿主产生影响。文章使用生物信息学的方法对无颌类、软骨鱼纲、肉鳍鱼纲和辐鳍鱼纲鱼类基因组进行了MITEs预测, 最终在33种鱼类基因组中鉴定出2433个MITEs家族。文章发现鱼类基因组中MITEs含量存在较大差异(0.11%—21.18%), 并且MITEs含量与鱼类基因组大小呈正相关关系。根据末端重复序列(TIRs)和靶位点重复序列(TSDs)的特征将MITEs分为10个超家族, 其中TC1-Mariner超家族的含量最高。MITEs在鱼类基因组中的插入事件主要发生在4百万年前至今, 大多数物种的MITEs在2百万—0.5百万年前发生了爆发式扩增。鱼类基因组中的MITEs多数插入到基因内部或附近, 这些转座子可能在基因的表达调控方面存在重要作用。     

Leishmania: conserved evolution--diverse diseases

The landmark completion of the Leishmania major genome sequence and the recent publication of the L. infantum and L. braziliensis genomes revealed the surprising result that, although separated by 15-50 million years of evolution, the Leishmania genomes are highly conserved and have less than 1% species-specific genes. Yet, these three species of Leishmania cause distinctive and diverse diseases in humans. Here, we discuss these findings together with recent microarray and proteomics studies and highlight their importance in understanding Leishmania disease phenotypes.     

Genomics of a phylum distant from flowering plants: conifers

Conifers are evolutionarily distant from angiosperms, separated by 300 million years of evolution. The genomes of coniferous species are very large, among the largest of any nonpolyploid plant species. Their genomes are characterized by reduced evolutionary rate for coding genes, accumulation of noncoding DNA, and evolutionarily distance from angiosperms. I highlight both the advantages and disadvantages for conifers as model organism for genomics. With advances of new high-throughput sequencing technologies, we are at a watershed in conifer genomics.     

Landscape of mobile genetic elements and their antibiotic resistance cargo in prokaryotic genomes

Prokaryotic Mobile Genetic Elements (MGEs) such as transposons, integrons, phages and plasmids, play important roles in prokaryotic evolution and in the dispersal of cargo functions like antibiotic resistance. However, each of these MGE types is usually annotated and analysed individually, hampering a global understanding of phylogenetic and environmental patterns of MGE dispersal. We thus developed a computational framework that captures diverse MGE types, their cargos and MGE-mediated horizontal transfer events, using recombinases as ubiquitous MGE marker genes and pangenome information for MGE boundary estimation. Applied to ∼84k genomes with habitat annotation, we mapped 2.8 million MGE-specific recombinases to six operational MGE types, which together contain on average 13% of all the genes in a genome. Transposable elements (TEs) dominated across all taxa (∼1.7 million occurrences), outnumbering phages and phage-like elements (<0.4 million). We recorded numerous MGE-mediated horizontal transfer events across diverse phyla and habitats involving all MGE types, disentangled and quantified the extent of hitchhiking of TEs (17%) and integrons (63%) with other MGE categories, and established TEs as dominant carriers of antibiotic resistance genes. We integrated all these findings into a resource (proMGE.embl.de), which should facilitate future studies on the large mobile part of genomes and its horizontal dispersal.     

Estimate of the degradation potentials of cellulose,xylan, and chitin across global prokaryotic communities

Complex polysaccharides (e.g. cellulose, xylan, and chitin), the most abundant renewable biomass resources available on Earth, are mainly degraded by microorganisms in nature. However, little is known about the global distribution of the enzymes and microorganisms responsible for the degradation of cellulose, xylan, and chitin in natural environments. Through large-scale alignments between the sequences released by the Earth Microbiome Project and sequenced prokaryotic genomes, we determined that almost all prokaryotic communities have the functional potentials to degrade cellulose, xylan, and chitin. The median abundances of genes encoding putative cellulases, xylanases, and chitinases in global prokaryotic communities are 0.51 (0.17–1.01), 0.24 (0.05–0.57), and 0.33 (0.11–0.71) genes/cell, respectively, and the composition and abundance of these enzyme systems are environmentally varied. The taxonomic sources of the three enzymes are highly diverse within prokaryotic communities, and the main factor influencing the diversity is the community's alpha diversity index rather than gene abundance. Moreover, there are obvious differences in taxonomic sources among different communities, and most genera with degradation potentials are narrowly distributed. In conclusion, our analysis preliminarily depicts a panorama of cellulose-, xylan-, and chitin-degrading enzymatic systems across global prokaryotic communities.     

Recently mobilized transposons in the human and chimpanzee genomes

Transposable genetic elements are abundant in the genomes of most organisms, including humans. These endogenous mutagens can alter genes, promote genomic rearrangements, and may help to drive the speciation of organisms. In this study, we identified almost 11,000 transposon copies that are differentially present in the human and chimpanzee genomes. Most of these transposon copies were mobilized after the existence of a common ancestor of humans and chimpanzees, approximately 6 million years ago. Alu, L1, and SVA insertions accounted for >95% of the insertions in both species. Our data indicate that humans have supported higher levels of transposition than have chimpanzees during the past several million years and have amplified different transposon subfamilies. In both species, approximately 34% of the insertions were located within known genes. These insertions represent a form of species-specific genetic variation that may have contributed to the differential evolution of humans and chimpanzees. In addition to providing an initial overview of recently mobilized elements, our collections will be useful for assessing the impact of these insertions on their hosts and for studying the transposition mechanisms of these elements.     

Divergences in gene repertoire among the reference Prevotella genomes derived from distinct body sites of human

Background

The community composition of the human microbiome is known to vary at distinct anatomical niches. But little is known about the nature of variations, if any, at the genome/sub-genome levels of a specific microbial community across different niches. The present report aims to explore, as a case study, the variations in gene repertoire of 28 Prevotella reference genomes derived from different body-sites of human, as reported earlier by the Human Microbiome Consortium.

Results

The pan-genome for Prevotella remains “open”. On an average, 17% of predicted protein-coding genes of any particular Prevotella genome represent the conserved core genes, while the remaining 83% contribute to the flexible and singletons. The study reveals exclusive presence of 11798, 3673, 3348 and 934 gene families and exclusive absence of 17, 221, 115 and 645 gene families in Prevotella genomes derived from human oral cavity, gastro-intestinal tracts (GIT), urogenital tract (UGT) and skin, respectively. Distribution of various functional COG categories differs significantly among the habitat-specific genes. No niche-specific variations could be observed in distribution of KEGG pathways.

Conclusions

Prevotella genomes derived from different body sites differ appreciably in gene repertoire, suggesting that these microbiome components might have developed distinct genetic strategies for niche adaptation within the host. Each individual microbe might also have a component of its own genetic machinery for host adaptation, as appeared from the huge number of singletons.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1350-6) contains supplementary material, which is available to authorized users.     

人体微生物组计划

人体微生物组计划又称第二人类基因组计划。由美国国立卫生研究院立项资助,2007年正式启动。计划用5年时间耗资1.5亿美元完成900个人体微生物基因组测序。其目标是探索研究人类微生物组的可行性;研究人体微生物组变化与疾病健康的关系;同时为其它科学研究提供信息和技术支持。     

Continuity of the dsRNA genome of bluetongue virus

The duplex RNA genome of bluetongue virus, extracted under acidic conditions with sodium dodecyl sulfate and phenol, is an unfragmented continuous structure. Most genomes appeared as a rosette with 10 loops (genes) emanating from a central crescent- or doughnutshaped area. The genome has a mean length of ca. 10 μm which corresponds to a molecular weight of ca. 23 million.     

Echinochloa Chloroplast Genomes: Insights into the Evolution and Taxonomic Identification of Two Weedy Species

The genus Echinochloa (Poaceae) includes numerous problematic weeds that cause the reduction of crop yield worldwide. To date, DNA sequence information is still limited in the genus Echinochloa. In this study, we completed the entire chloroplast genomes of two Echinochloa species (Echinochloa oryzicola and Echinochloa crus-galli) based on high-throughput sequencing data from their fresh green leaves. The two Echinochloa chloroplast genomes are 139,891 and 139,800 base pairs in length, respectively, and contain 131 protein-coding genes, 79 indels and 466 substitutions helpful for discrimination of the two species. The divergence between the genus Echinochloa and Panicum occurred about 21.6 million years ago, whereas the divergence between E. oryzicola and E. crus-galli chloroplast genes occurred about 3.3 million years ago. The two reported Echinochloa chloroplast genome sequences contribute to better understanding of the diversification of this genus.