Tundra soil carbon is vulnerable to rapid microbial decomposition under climate warming

Microbial decomposition of soil carbon in high-latitude tundra underlain with permafrost is one of the most important, but poorly understood, potential positive feedbacks of greenhouse gas emissions from terrestrial ecosystems into the atmosphere in a warmer world. Using integrated metagenomic technologies, we showed that the microbial functional community structure in the active layer of tundra soil was significantly altered after only 1.5 years of warming, a rapid response demonstrating the high sensitivity of this ecosystem to warming. The abundances of microbial functional genes involved in both aerobic and anaerobic carbon decomposition were also markedly increased by this short-term warming. Consistent with this, ecosystem respiration (R eco) increased up to 38%. In addition, warming enhanced genes involved in nutrient cycling, which very likely contributed to an observed increase (30%) in gross primary productivity (GPP). However, the GPP increase did not offset the extra R eco, resulting in significantly more net carbon loss in warmed plots compared with control plots. Altogether, our results demonstrate the vulnerability of active-layer soil carbon in this permafrost-based tundra ecosystem to climate warming and the importance of microbial communities in mediating such vulnerability

Genomic fluidity: an integrative view of gene diversity within microbial populations

<p>Abstract</p> <p>Background</p> <p>The dual concepts of pan and core genomes have been widely adopted as means to assess the distribution of gene families within microbial species and genera. The core genome is the set of genes shared by a group of organisms; the pan genome is the set of all genes seen in any of these organisms. A variety of methods have provided drastically different estimates of the sizes of pan and core genomes from sequenced representatives of the same groups of bacteria.</p> <p>Results</p> <p>We use a combination of mathematical, statistical and computational methods to show that current predictions of pan and core genome sizes may have no correspondence to true values. Pan and core genome size estimates are problematic because they depend on the estimation of the occurrence of rare genes and genomes, respectively, which are difficult to estimate precisely because they are rare. Instead, we introduce and evaluate a robust metric - genomic fluidity - to categorize the gene-level similarity among groups of sequenced isolates. Genomic fluidity is a measure of the dissimilarity of genomes evaluated at the gene level.</p> <p>Conclusions</p> <p>The genomic fluidity of a population can be estimated accurately given a small number of sequenced genomes. Further, the genomic fluidity of groups of organisms can be compared robustly despite variation in algorithms used to identify genes and their homologs. As such, we recommend that genomic fluidity be used in place of pan and core genome size estimates when assessing gene diversity within genomes of a species or a group of closely related organisms.</p

Methanocella conradii sp. nov., a Thermophilic, Obligate Hydrogenotrophic Methanogen, Isolated from Chinese Rice Field Soil

BACKGROUND: Methanocellales contributes significantly to anthropogenic methane emissions that cause global warming, but few pure cultures for Methanocellales are available to permit subsequent laboratory studies (physiology, biochemistry, etc.). METHODOLOGY/PRINCIPAL FINDINGS: By combining anaerobic culture and molecular techniques, a novel thermophilic methanogen, strain HZ254(T) was isolated from a Chinese rice field soil located in Hangzhou, China. The phylogenetic analyses of both the 16S rRNA gene and mcrA gene (encoding the α subunit of methyl-coenzyme M reductase) confirmed its affiliation with Methanocellales, and Methanocella paludicola SANAE(T) was the most closely related species. Cells were non-motile rods, albeit with a flagellum, 1.4-2.8 µm long and by 0.2-0.3 µm in width. They grew at 37-60 °C (optimally at 55 °C) and salinity of 0-5 g NaCl l(-1) (optimally at 0-1 g NaCl l(-1)). The pH range for growth was 6.4-7.2 (optimum 6.8). Under the optimum growth condition, the doubling time was 6.5-7.8 h, which is the shortest ever observed in Methanocellales. Strain HZ254(T) utilized H(2)/CO(2) but not formate for growth and methane production. The DNA G+C content of this organism was 52.7 mol%. The sequence identities of 16S rRNA gene and mcrA gene between strain HZ254(T) and SANAE(T) were 95.0 and 87.5% respectively, and the genome based Average Nucleotide Identity value between them was 74.8%. These two strains differed in phenotypic features with regard to substrate utilization, possession of a flagellum, doubling time (under optimal conditions), NaCl and temperature ranges. Taking account of the phenotypic and phylogenetic characteristics, we propose strain HZ254(T) as a representative of a novel species, Methanocella conradii sp. nov. The type strain is HZ254(T) ( = CGMCC 1.5162(T) = JCM 17849(T) = DSM 24694(T)). CONCLUSIONS/SIGNIFICANCE: Strain HZ254(T) could potentially serve as an excellent laboratory model for studying Methanocellales due to its fast growth and consistent cultivability

Genomic Expansion of Magnetotactic Bacteria Reveals an Early Common Origin of Magnetotaxis with Lineage-specific Evolution

The origin and evolution of magnetoreception, which in diverse prokaryotes and protozoa is known as magnetotaxis and enables these microorganisms to detect Earth’s magnetic field for orientation and navigation, is not well understood in evolutionary biology. The only known prokaryotes capable of sensing the geomagnetic field are magnetotactic bacteria (MTB), motile microorganisms that biomineralize intracellular, membrane-bounded magnetic single-domain crystals of either magnetite (Fe3O4) or greigite (Fe3S4) called magnetosomes. Magnetosomes are responsible for magnetotaxis in MTB. Here we report the first large-scale metagenomic survey of MTB from both northern and southern hemispheres combined with 28 genomes from uncultivated MTB. These genomes expand greatly the coverage of MTB in the Proteobacteria, Nitrospirae, and Omnitrophica phyla, and provide the first genomic evidence of MTB belonging to the Zetaproteobacteria and “Candidatus Lambdaproteobacteria” classes. The gene content and organization of magnetosome gene clusters, which are physically grouped genes that encode proteins for magnetosome biosynthesis and organization, are more conserved within phylogenetically similar groups than between different taxonomic lineages. Moreover, the phylogenies of core magnetosome proteins form monophyletic clades. Together, these results suggest a common ancient origin of iron-based (Fe3O4 and Fe3S4) magnetotaxis in the domain Bacteria that underwent lineage-specific evolution, shedding new light on the origin and evolution of biomineralization and magnetotaxis, and expanding significantly the phylogenomic representation of MTB

Bacterial microevolution and the Pangenome

The comparison of multiple genome sequences sampled from a bacterial population reveals considerable diversity in both the core and the accessory parts of the pangenome. This diversity can be analysed in terms of microevolutionary events that took place since the genomes shared a common ancestor, especially deletion, duplication, and recombination. We review the basic modelling ingredients used implicitly or explicitly when performing such a pangenome analysis. In particular, we describe a basic neutral phylogenetic framework of bacterial pangenome microevolution, which is not incompatible with evaluating the role of natural selection. We survey the different ways in which pangenome data is summarised in order to be included in microevolutionary models, as well as the main methodological approaches that have been proposed to reconstruct pangenome microevolutionary history

From learning taxonomies to phylogenetic learning: Integration of 16S rRNA gene data into FAME-based bacterial classification

<p>Abstract</p> <p>Background</p> <p>Machine learning techniques have shown to improve bacterial species classification based on fatty acid methyl ester (FAME) data. Nonetheless, FAME analysis has a limited resolution for discrimination of bacteria at the species level. In this paper, we approach the species classification problem from a taxonomic point of view. Such a taxonomy or tree is typically obtained by applying clustering algorithms on FAME data or on 16S rRNA gene data. The knowledge gained from the tree can then be used to evaluate FAME-based classifiers, resulting in a novel framework for bacterial species classification.</p> <p>Results</p> <p>In view of learning in a taxonomic framework, we consider two types of trees. First, a FAME tree is constructed with a supervised divisive clustering algorithm. Subsequently, based on 16S rRNA gene sequence analysis, phylogenetic trees are inferred by the NJ and UPGMA methods. In this second approach, the species classification problem is based on the combination of two different types of data. Herein, 16S rRNA gene sequence data is used for phylogenetic tree inference and the corresponding binary tree splits are learned based on FAME data. We call this learning approach 'phylogenetic learning'. Supervised Random Forest models are developed to train the classification tasks in a stratified cross-validation setting. In this way, better classification results are obtained for species that are typically hard to distinguish by a single or flat multi-class classification model.</p> <p>Conclusions</p> <p>FAME-based bacterial species classification is successfully evaluated in a taxonomic framework. Although the proposed approach does not improve the overall accuracy compared to flat multi-class classification, it has some distinct advantages. First, it has better capabilities for distinguishing species on which flat multi-class classification fails. Secondly, the hierarchical classification structure allows to easily evaluate and visualize the resolution of FAME data for the discrimination of bacterial species. Summarized, by phylogenetic learning we are able to situate and evaluate FAME-based bacterial species classification in a more informative context.</p

A new approach to in silico SNP detection and some new SNPs in the Bacillus anthracis genome

<p>Abstract</p> <p>Background</p> <p><it>Bacillus anthracis </it>is one of the most monomorphic pathogens known. Identification of polymorphisms in its genome is essential for taxonomic classification, for determination of recent evolutionary changes, and for evaluation of pathogenic potency.</p> <p>Findings</p> <p>In this work three strains of the <it>Bacillus anthracis </it>genome are compared and previously unpublished single nucleotide polymorphisms (SNPs) are revealed. Moreover, it is shown that, despite the highly monomorphic nature of <it>Bacillus anthracis</it>, the SNPs are (1) abundant in the genome and (2) distributed relatively uniformly across the sequence.</p> <p>Conclusions</p> <p>The findings support the proposition that SNPs, together with indels and variable number tandem repeats (VNTRs), can be used effectively not only for the differentiation of perfect strain data, but also for the comparison of moderately incomplete, noisy and, in some cases, unknown <it>Bacillus anthracis </it>strains. In the case when the data is of still lower quality, a new DNA sequence fingerprinting approach based on recently introduced markers, based on combinatorial-analytic concepts and called cyclic difference sets, can be used.</p

Going Deeper: Metagenome of a Hadopelagic Microbial Community

The paucity of sequence data from pelagic deep-ocean microbial assemblages has severely restricted molecular exploration of the largest biome on Earth. In this study, an analysis is presented of a large-scale 454-pyrosequencing metagenomic dataset from a hadopelagic environment from 6,000 m depth within the Puerto Rico Trench (PRT). A total of 145 Mbp of assembled sequence data was generated and compared to two pelagic deep ocean metagenomes and two representative surface seawater datasets from the Sargasso Sea. In a number of instances, all three deep metagenomes displayed similar trends, but were most magnified in the PRT, including enrichment in functions for two-component signal transduction mechanisms and transcriptional regulation. Overrepresented transporters in the PRT metagenome included outer membrane porins, diverse cation transporters, and di- and tri-carboxylate transporters that matched well with the prevailing catabolic processes such as butanoate, glyoxylate and dicarboxylate metabolism. A surprisingly high abundance of sulfatases for the degradation of sulfated polysaccharides were also present in the PRT. The most dramatic adaptational feature of the PRT microbes appears to be heavy metal resistance, as reflected in the large numbers of transporters present for their removal. As a complement to the metagenome approach, single-cell genomic techniques were utilized to generate partial whole-genome sequence data from four uncultivated cells from members of the dominant phyla within the PRT, Alphaproteobacteria, Gammaproteobacteria, Bacteroidetes and Planctomycetes. The single-cell sequence data provided genomic context for many of the highly abundant functional attributes identified from the PRT metagenome, as well as recruiting heavily the PRT metagenomic sequence data compared to 172 available reference marine genomes. Through these multifaceted sequence approaches, new insights have been provided into the unique functional attributes present in microbes residing in a deeper layer of the ocean far removed from the more productive sun-drenched zones above

Mining Virulence Genes Using Metagenomics

When a bacterial genome is compared to the metagenome of an environment it inhabits, most genes recruit at high sequence identity. In free-living bacteria (for instance marine bacteria compared against the ocean metagenome) certain genomic regions are totally absent in recruitment plots, representing therefore genes unique to individual bacterial isolates. We show that these Metagenomic Islands (MIs) are also visible in bacteria living in human hosts when their genomes are compared to sequences from the human microbiome, despite the compartmentalized structure of human-related environments such as the gut. From an applied point of view, MIs of human pathogens (e.g. those identified in enterohaemorragic Escherichia coli against the gut metagenome or in pathogenic Neisseria meningitidis against the oral metagenome) include virulence genes that appear to be absent in related strains or species present in the microbiome of healthy individuals. We propose that this strategy (i.e. recruitment analysis of pathogenic bacteria against the metagenome of healthy subjects) can be used to detect pathogenicity regions in species where the genes involved in virulence are poorly characterized. Using this approach, we detect well-known pathogenicity islands and identify new potential virulence genes in several human pathogens

Horizontal DNA transfer mechanisms of bacteria as weapons of intragenomic conflict

Horizontal DNA transfer (HDT) is a pervasive mechanism of diversification in many microbial species, but its primary evolutionary role remains controversial. Much recent research has emphasised the adaptive benefit of acquiring novel DNA, but here we argue instead that intragenomic conflict provides a coherent framework for understanding the evolutionary origins of HDT. To test this hypothesis, we developed a mathematical model of a clonally descended bacterial population undergoing HDT through transmission of mobile genetic elements (MGEs) and genetic transformation. Including the known bias of transformation toward the acquisition of shorter alleles into the model suggested it could be an effective means of counteracting the spread of MGEs. Both constitutive and transient competence for transformation were found to provide an effective defence against parasitic MGEs; transient competence could also be effective at permitting the selective spread of MGEs conferring a benefit on their host bacterium. The coordination of transient competence with cell-cell killing, observed in multiple species, was found to result in synergistic blocking of MGE transmission through releasing genomic DNA for homologous recombination while simultaneously reducing horizontal MGE spread by lowering the local cell density. To evaluate the feasibility of the functions suggested by the modelling analysis, we analysed genomic data from longitudinal sampling of individuals carrying Streptococcus pneumoniae. This revealed the frequent within-host coexistence of clonally descended cells that differed in their MGE infection status, a necessary condition for the proposed mechanism to operate. Additionally, we found multiple examples of MGEs inhibiting transformation through integrative disruption of genes encoding the competence machinery across many species, providing evidence of an ongoing "arms race." Reduced rates of transformation have also been observed in cells infected by MGEs that reduce the concentration of extracellular DNA through secretion of DNases. Simulations predicted that either mechanism of limiting transformation would benefit individual MGEs, but also that this tactic's effectiveness was limited by competition with other MGEs coinfecting the same cell. A further observed behaviour we hypothesised to reduce elimination by transformation was MGE activation when cells become competent. Our model predicted that this response was effective at counteracting transformation independently of competing MGEs. Therefore, this framework is able to explain both common properties of MGEs, and the seemingly paradoxical bacterial behaviours of transformation and cell-cell killing within clonally related populations, as the consequences of intragenomic conflict between self-replicating chromosomes and parasitic MGEs. The antagonistic nature of the different mechanisms of HDT over short timescales means their contribution to bacterial evolution is likely to be substantially greater than previously appreciated