Community transcriptomics reveals universal patterns of protein sequence conservation in natural microbial communities

Metatranscriptomes reveal the diverse responses of Thaumarchaeota ecotypes to environmental variations in the northern slope of the South China Sea

Thaumarchaeota are among the most abundant prokaryotes in the ocean, playing important roles in carbon and nitrogen cycling. Marine Thaumarchaeota ecotypes exhibit depth-related diversification and seasonal changes. However, transcriptomic activities concerning niche partitioning among thaumarchaeal ecotypes remain unclear. Here, we examined the variations in the distribution and transcriptomic activity of marine Thaumarchaeota ecotypes. Three primary ecotypes were identified: a Nitrosopumilus-like clade; a Nitrosopelagicus-like water column A (WCA) clade, thriving in epipelagic water; and a water column B (WCB) clade, dominant in deep water. Depth-related partitioning of the three ecotypes and the seasonal variability of the WCA and WCB ecotypes were observed. Nutrient concentrations, chlorophyll α and salinity were the primary environmental factors. The relative abundance of the WCA ecotype and its transcript abundance of amoA gene were positively correlated with chlorophyll α and salinity, while the WCB ecotype was positively correlated with nitrate and phosphate. Based on high-quality metagenome-assembled genomes, transcriptomic analysis revealed that the three ecotypes exhibited various co-occurring expression patterns of the elemental cycling genes in the nitrogen, carbon, phosphorus, and sulfur cycles. Our results provide transcriptomic evidence of the niche differentiation of marine Thaumarchaeota ecotypes, highlighting the diverse roles of ecotypes and WCA subclades in biogeochemical cycles.

Metabolic Pathway Prediction Using Non-Negative Matrix Factorization with Improved Precision

Machine learning provides a probabilistic framework for metabolic pathway inference from genomic sequence information at different levels of complexity and completion. However, several challenges, including pathway features engineering, multiple mapping of enzymatic reactions, and emergent or distributed metabolism within populations or communities of cells, can limit prediction performance. In this article, we present triUMPF (triple non-negative matrix factorization [NMF] with community detection for metabolic pathway inference), which combines three stages of NMF to capture myriad relationships between enzymes and pathways within a graph network. This is followed by community detection to extract a higher-order structure based on the clustering of vertices that share similar statistical properties. We evaluated triUMPF performance by using experimental datasets manifesting diverse multi-label properties, including Tier 1 genomes from the BioCyc collection of organismal Pathway/Genome Databases and low complexity microbial communities. Resulting performance metrics equaled or exceeded other prediction methods on organismal genomes with improved precision on multi-organismal datasets.

Unlocking the potential of insect and ruminant host symbionts for recycling of lignocellulosic carbon with a biorefinery approach: a review

Uprising fossil fuel depletion and deterioration of ecological reserves supply have led to the search for alternative renewable and sustainable energy sources and chemicals. Although first generation biorefinery is quite successful commercially in generating bulk of biofuels globally, the food versus fuel debate has necessitated the use of non-edible feedstocks, majorly waste biomass, for second generation production of biofuels and chemicals. A diverse class of microbes and enzymes are being exploited for biofuels production for a series of treatment process, however, the conversion efficiency of wide range of lignocellulosic biomass (LCB) and consolidated way of processing remains challenging. There were lot of research efforts in the past decade to scour for potential microbial candidate. In this context, evolution has developed the gut microbiota of several insects and ruminants that are potential LCB degraders host eco-system to overcome its host nutritional constraints, where LCB processed by microbiomes pretends to be a promising candidate. Synergistic microbial symbionts could make a significant contribution towards recycling the renewable carbon from distinctly abundant recalcitrant LCB. Several studies have assessed the bioprospection of innumerable gut symbionts and their lignocellulolytic enzymes for LCB degradation. Though, some reviews exist on molecular characterization of gut microbes, but none of them has enlightened the microbial community design coupled with various LCB valorization which intensifies the microbial diversity in biofuels application. This review provides a deep insight into the significant breakthroughs attained in enrichment strategy of gut microbial community and its molecular characterization techniques which aids in understanding the holistic microbial community dynamics. Special emphasis is placed on gut microbial role in LCB depolymerization strategies to lignocellulolytic enzymes production and its functional metagenomic data mining eventually generating the sugar platform for biofuels and renewable chemicals production.

Metabolic pathway inference using multi-label classification with rich pathway features

Metabolic inference from genomic sequence information is a necessary step in determining the capacity of cells to make a living in the world at different levels of biological organization. A common method for determining the metabolic potential encoded in genomes is to map conceptually translated open reading frames onto a database containing known product descriptions. Such gene-centric methods are limited in their capacity to predict pathway presence or absence and do not support standardized rule sets for automated and reproducible research. Pathway-centric methods based on defined rule sets or machine learning algorithms provide an adjunct or alternative inference method that supports hypothesis generation and testing of metabolic relationships within and between cells. Here, we present mlLGPR, multi-label based on logistic regression for pathway prediction, a software package that uses supervised multi-label classification and rich pathway features to infer metabolic networks in organismal and multi-organismal datasets. We evaluated mlLGPR performance using a corpora of 12 experimental datasets manifesting diverse multi-label properties, including manually curated organismal genomes, synthetic microbial communities and low complexity microbial communities. Resulting performance metrics equaled or exceeded previous reports for organismal genomes and identify specific challenges associated with features engineering and training data for community-level metabolic inference.

Species specific gene expression dynamics during harmful algal blooms

Harmful algal blooms are caused by specific members of microbial communities. Understanding the dynamics of these events requires comparing the strategies developed by the problematic species to cope with environmental fluctuations to the ones developed by the other members of the community. During three consecutive years, the meta-transcriptome of micro-eukaryote communities was sequenced during blooms of the toxic dinoflagellate Alexandrium minutum. The dataset was analyzed to investigate species specific gene expression dynamics. Major shifts in gene expression were explained by the succession of different species within the community. Although expression patterns were strongly correlated with fluctuation of the abiotic environment, and more specifically with nutrient concentration, transcripts specifically involved in nutrient uptake and metabolism did not display extensive changes in gene expression. Compared to the other members of the community, A. minutum displayed a very specific expression pattern, with lower expression of photosynthesis transcripts and central metabolism genes (TCA cycle, glucose metabolism, glycolysis…) and contrasting expression pattern of ion transporters across environmental conditions. These results suggest the importance of mixotrophy, cell motility and cell-to-cell interactions during A. minutum blooms.

Diel Oscillation of Microbial Gene Transcripts Declines With Depth in Oligotrophic Ocean Waters

Diel oscillations in primary and secondary production, growth, metabolic activity, and gene expression commonly occur in marine microbial communities in ocean surface waters. Diel periodicity of gene transcription has been demonstrated in photoautotrophic and heterotrophic microbes in both coastal and open ocean environments. To better define the spatiotemporal distribution and patterns of these daily oscillations, we investigated how diel periodicity in gene transcripts changed with depth from the surface waters to the upper mesopelagic. We postulated that diel oscillation of transcript abundances would diminish at greater depths across the collective microbial community due to decreasing light availability. The results showed that the number and total proportion of gene transcripts and taxa exhibiting diel periodicity were greatest in the shallow sunlit mixed layer, diminished rapidly with increasing depth to the base of the euphotic zone, and could not be detected in the mesopelagic. The results confirmed an overall decrease in microbial diel transcript oscillation with depth through the euphotic zone and suggested a relationship between abundance of diel oscillating transcripts and the daily integrated light exposure experienced by planktonic microbes in the water column. Local dissolved macronutrient concentration also appeared to influence the diel transcriptional patterns of specific microbial genes. The diminishing diel transcript oscillations found at increasing depths suggest that diel patterns of other microbial processes and interactions may likewise be attenuated at depth.

Taxonomic and Functional Compositions Impacted by the Quality of Metatranscriptomic Assemblies

Metatranscriptomics has recently been applied to investigate the active biogeochemical processes and elemental cycles, and in situ responses of microbiomes to environmental stimuli and stress factors. De novo assembly of RNA-Sequencing (RNA-Seq) data can reveal a more detailed description of the metabolic interactions amongst the active microbial communities. However, the quality of the assemblies and the depiction of the metabolic network provided by various de novo assemblers have not yet been thoroughly assessed. In this study, we compared 15 de novo metatranscriptomic assemblies for a fracture fluid sample collected from a borehole located at 1.34 km below land surface in a South African gold mine. These assemblies were constructed from total, non-coding, and coding reads using five de novo transcriptomic assemblers (Trans-ABySS, Trinity, Oases, IDBA-tran, and Rockhopper). They were evaluated based on the number of transcripts, transcript length, range of transcript coverage, continuity, percentage of transcripts with confident annotation assignments, as well as taxonomic and functional diversity patterns. The results showed that these parameters varied considerably among the assemblies, with Trans-ABySS and Trinity generating the best assemblies for non-coding and coding RNA reads, respectively, because the high number of transcripts assembled covered a wide expression range, and captured extensively the taxonomic and metabolic gene diversity, respectively. We concluded that the choice of de novo transcriptomic assemblers impacts substantially the taxonomic and functional compositions. Care should be taken to obtain high-quality assemblies for informing the in situ metabolic landscape.

Comparative Metatranscriptomics of Wheat Rhizosphere Microbiomes in Disease Suppressive and Non-suppressive Soils for Rhizoctonia solani AG8

The soilborne fungus Rhizoctonia solani anastomosis group (AG) 8 is a major pathogen of grain crops resulting in substantial production losses. In the absence of resistant cultivars of wheat or barley, a sustainable and enduring method for disease control may lie in the enhancement of biological disease suppression. Evidence of effective biological control of R. solani AG8 through disease suppression has been well documented at our study site in Avon, South Australia. A comparative metatranscriptomic approach was applied to assess the taxonomic and functional characteristics of the rhizosphere microbiome of wheat plants grown in adjacent fields which are suppressive and non-suppressive to the plant pathogen R. solani AG8. Analysis of 12 rhizosphere metatranscriptomes (six per field) was undertaken using two bioinformatic approaches involving unassembled and assembled reads. Differential expression analysis showed the dominant taxa in the rhizosphere based on mRNA annotation were Arthrobacter spp. and Pseudomonas spp. for non-suppressive samples and Stenotrophomonas spp. and Buttiauxella spp. for the suppressive samples. The assembled metatranscriptome analysis identified more differentially expressed genes than the unassembled analysis in the comparison of suppressive and non-suppressive samples. Suppressive samples showed greater expression of a polyketide cyclase, a terpenoid biosynthesis backbone gene (dxs) and many cold shock proteins (csp). Non-suppressive samples were characterised by greater expression of antibiotic genes such as non-heme chloroperoxidase (cpo) which is involved in pyrrolnitrin synthesis, and phenazine biosynthesis family protein F (phzF) and its transcriptional activator protein (phzR). A large number of genes involved in detoxifying reactive oxygen species (ROS) and superoxide radicals (sod, cat, ahp, bcp, gpx1, trx) were also expressed in the non-suppressive rhizosphere samples most likely in response to the infection of wheat roots by R. solani AG8. Together these results provide new insight into microbial gene expression in the rhizosphere of wheat in soils suppressive and non-suppressive to R. solani AG8. The approach taken and the genes involved in these functions provide direction for future studies to determine more precisely the molecular interplay of plant-microbe-pathogen interactions with the ultimate goal of the development of management options that promote beneficial rhizosphere microflora to reduce R. solani AG8 infection of crops.

Metagenomic insights into the ecology and physiology of microbes in bioelectrochemical systems

In bioelectrochemical systems (BESs), electrons are transferred between electrochemically active microbes (EAMs) and conductive materials, such as electrodes, via extracellular electron transfer (EET) pathways, and electrons thus transferred stimulate intracellular catabolic reactions. Catabolic and EET pathways have extensively been studied for several model EAMs, such as Shewanella oneidensis MR-1 and Geobacter sulfurreducens PCA, whereas it is also important to understand the ecophysiology of EAMs in naturally occurring microbiomes, such as those in anode biofilms in microbial fuel cells treating wastewater. Recent studies have exploited metagenomics and metatranscriptomics (meta-omics) approaches to characterize EAMs in BES-associated microbiomes. Here we review recent BES studies that used meta-omics approaches and show that these studies have discovered unexpected features of EAMs and deepened our understanding of functions and behaviors of microbes in BESs. It is desired that more studies will employ meta-omics approaches for advancing our knowledge on microbes in BESs.

SAR11 bacteria linked to ocean anoxia and nitrogen loss

Bacteria of the SAR11 clade constitute up to one half of all microbial cells in the oxygen-rich surface ocean. SAR11 bacteria are also abundant in oxygen minimum zones (OMZs), where oxygen falls below detection and anaerobic microbes have vital roles in converting bioavailable nitrogen to N2 gas. Anaerobic metabolism has not yet been observed in SAR11, and it remains unknown how these bacteria contribute to OMZ biogeochemical cycling. Here, genomic analysis of single cells from the world's largest OMZ revealed previously uncharacterized SAR11 lineages with adaptations for life without oxygen, including genes for respiratory nitrate reductases (Nar). SAR11 nar genes were experimentally verified to encode proteins catalysing the nitrite-producing first step of denitrification and constituted ~40% of OMZ nar transcripts, with transcription peaking in the anoxic zone of maximum nitrate reduction activity. These results link SAR11 to pathways of ocean nitrogen loss, redefining the ecological niche of Earth's most abundant organismal group.

Using "Omics" and Integrated Multi-Omics Approaches to Guide Probiotic Selection to Mitigate Chytridiomycosis and Other Emerging Infectious Diseases

Emerging infectious diseases in wildlife are responsible for massive population declines. In amphibians, chytridiomycosis caused by Batrachochytrium dendrobatidis, Bd, has severely affected many amphibian populations and species around the world. One promising management strategy is probiotic bioaugmentation of antifungal bacteria on amphibian skin. In vivo experimental trials using bioaugmentation strategies have had mixed results, and therefore a more informed strategy is needed to select successful probiotic candidates. Metagenomic, transcriptomic, and metabolomic methods, colloquially called "omics," are approaches that can better inform probiotic selection and optimize selection protocols. The integration of multiple omic data using bioinformatic and statistical tools and in silico models that link bacterial community structure with bacterial defensive function can allow the identification of species involved in pathogen inhibition. We recommend using 16S rRNA gene amplicon sequencing and methods such as indicator species analysis, the Kolmogorov-Smirnov Measure, and co-occurrence networks to identify bacteria that are associated with pathogen resistance in field surveys and experimental trials. In addition to 16S amplicon sequencing, we recommend approaches that give insight into symbiont function such as shotgun metagenomics, metatranscriptomics, or metabolomics to maximize the probability of finding effective probiotic candidates, which can then be isolated in culture and tested in persistence and clinical trials. An effective mitigation strategy to ameliorate chytridiomycosis and other emerging infectious diseases is necessary; the advancement of omic methods and the integration of multiple omic data provide a promising avenue toward conservation of imperiled species.

Metatranscriptomic analyses of plant cell wall polysaccharide degradation by microorganisms in the cow rumen

The bovine rumen represents a highly specialized bioreactor where plant cell wall polysaccharides (PCWPs) are efficiently deconstructed via numerous enzymes produced by resident microorganisms. Although a large number of fibrolytic genes from rumen microorganisms have been identified, it remains unclear how they are expressed in a coordinated manner to efficiently degrade PCWPs. In this study, we performed a metatranscriptomic analysis of the rumen microbiomes of adult Holstein cows fed a fiber diet and obtained a total of 1,107,083 high-quality non-rRNA reads with an average length of 483 nucleotides. Transcripts encoding glycoside hydrolases (GHs) and carbohydrate binding modules (CBMs) accounted for 1% and 0.1% of the total non-rRNAs, respectively. The majority (98%) of the putative cellulases belonged to four GH families (i.e., GH5, GH9, GH45, and GH48) and were primarily synthesized by Ruminococcus and Fibrobacter. Notably, transcripts for GH48 cellobiohydrolases were relatively abundant compared to the abundance of transcripts for other cellulases. Two-thirds of the putative hemicellulases were of the GH10, GH11, and GH26 types and were produced by members of the genera Ruminococcus, Prevotella, and Fibrobacter. Most (82%) predicted oligosaccharide-degrading enzymes were GH1, GH2, GH3, and GH43 proteins and were from a diverse group of microorganisms. Transcripts for CBM10 and dockerin, key components of the cellulosome, were also relatively abundant. Our results provide metatranscriptomic evidence in support of the notion that members of the genera Ruminococcus, Fibrobacter, and Prevotella are predominant PCWP degraders and point to the significant contribution of GH48 cellobiohydrolases and cellulosome-like structures to efficient PCWP degradation in the cow rumen.

Diversification and niche adaptations of Nitrospina-like bacteria in the polyextreme interfaces of Red Sea brines

Nitrite-oxidizing bacteria (NOB) of the genus Nitrospina have exclusively been found in marine environments. In the brine–seawater interface layer of Atlantis II Deep (Red Sea), Nitrospina-like bacteria constitute up to one-third of the bacterial 16S ribosomal RNA (rRNA) gene sequences. This is much higher compared with that reported in other marine habitats (~10% of all bacteria), and was unexpected because no NOB culture has been observed to grow above 4.0% salinity, presumably due to the low net energy gained from their metabolism that is insufficient for both growth and osmoregulation. Using phylogenetics, single-cell genomics and metagenomic fragment recruitment approaches, we document here that these Nitrospina-like bacteria, designated as Candidatus Nitromaritima RS, are not only highly diverged from the type species Nitrospina gracilis (pairwise genome identity of 69%) but are also ubiquitous in the deeper, highly saline interface layers (up to 11.2% salinity) with temperatures of up to 52 °C. Comparative pan-genome analyses revealed that less than half of the predicted proteome of Ca. Nitromaritima RS is shared with N. gracilis. Interestingly, the capacity for nitrite oxidation is also conserved in both genomes. Although both lack acidic proteomes synonymous with extreme halophiles, the pangenome of Ca. Nitromaritima RS specifically encodes enzymes with osmoregulatory and thermoprotective roles (i.e., ectoine/hydroxyectoine biosynthesis) and of thermodynamic importance (i.e., nitrate and nitrite reductases). Ca. Nitromaritima RS also possesses many hallmark traits of microaerophiles and high-affinity NOB. The abundance of the uncultured Ca. Nitromaritima lineage in marine oxyclines suggests their unrecognized ecological significance in deoxygenated areas of the global ocean.

Microbial metabolic networks in a complex electrogenic biofilm recovered from a stimulus-induced metatranscriptomics approach

Microorganisms almost always exist as mixed communities in nature. While the significance of microbial community activities is well appreciated, a thorough understanding about how microbial communities respond to environmental perturbations has not yet been achieved. Here we have used a combination of metagenomic, genome binning, and stimulus-induced metatranscriptomic approaches to estimate the metabolic network and stimuli-induced metabolic switches existing in a complex microbial biofilm that was producing electrical current via extracellular electron transfer (EET) to a solid electrode surface. Two stimuli were employed: to increase EET and to stop EET. An analysis of cell activity marker genes after stimuli exposure revealed that only two strains within eleven binned genomes had strong transcriptional responses to increased EET rates, with one responding positively and the other responding negatively. Potential metabolic switches between eleven dominant members were mainly observed for acetate, hydrogen, and ethanol metabolisms. These results have enabled the estimation of a multi-species metabolic network and the associated short-term responses to EET stimuli that induce changes to metabolic flow and cooperative or competitive microbial interactions. This systematic meta-omics approach represents a next step towards understanding complex microbial roles within a community and how community members respond to specific environmental stimuli.

Mining microbial metatranscriptomes for expression of antibiotic resistance genes under natural conditions

Antibiotic resistance genes are found in a broad range of ecological niches associated with complex microbiota. Here we investigated if resistance genes are not only present, but also transcribed under natural conditions. Furthermore, we examined the potential for antibiotic production by assessing the expression of associated secondary metabolite biosynthesis gene clusters. Metatranscriptome datasets from intestinal microbiota of four human adults, one human infant, 15 mice and six pigs, of which only the latter have received antibiotics prior to the study, as well as from sea bacterioplankton, a marine sponge, forest soil and sub-seafloor sediment, were investigated. We found that resistance genes are expressed in all studied ecological niches, albeit with niche-specific differences in relative expression levels and diversity of transcripts. For example, in mice and human infant microbiota predominantly tetracycline resistance genes were expressed while in human adult microbiota the spectrum of expressed genes was more diverse, and also included β-lactam, aminoglycoside and macrolide resistance genes. Resistance gene expression could result from the presence of natural antibiotics in the environment, although we could not link it to expression of corresponding secondary metabolites biosynthesis clusters. Alternatively, resistance gene expression could be constitutive, or these genes serve alternative roles besides antibiotic resistance.

Ecological roles of dominant and rare prokaryotes in acid mine drainage revealed by metagenomics and metatranscriptomics

High-throughput sequencing is expanding our knowledge of microbial diversity in the environment. Still, understanding the metabolic potentials and ecological roles of rare and uncultured microbes in natural communities remains a major challenge. To this end, we applied a 'divide and conquer' strategy that partitioned a massive metagenomic data set (>100 Gbp) into subsets based on K-mer frequency in sequence assembly to a low-diversity acid mine drainage (AMD) microbial community and, by integrating with an additional metatranscriptomic assembly, successfully obtained 11 draft genomes most of which represent yet uncultured and/or rare taxa (relative abundance <1%). We report the first genome of a naturally occurring Ferrovum population (relative abundance >90%) and its metabolic potentials and gene expression profile, providing initial molecular insights into the ecological role of these lesser known, but potentially important, microorganisms in the AMD environment. Gene transcriptional analysis of the active taxa revealed major metabolic capabilities executed in situ, including carbon- and nitrogen-related metabolisms associated with syntrophic interactions, iron and sulfur oxidation, which are key in energy conservation and AMD generation, and the mechanisms of adaptation and response to the environmental stresses (heavy metals, low pH and oxidative stress). Remarkably, nitrogen fixation and sulfur oxidation were performed by the rare taxa, indicating their critical roles in the overall functioning and assembly of the AMD community. Our study demonstrates the potential of the 'divide and conquer' strategy in high-throughput sequencing data assembly for genome reconstruction and functional partitioning analysis of both dominant and rare species in natural microbial assemblages.

Differential assemblage of functional units in paddy soil microbiomes

Flooded rice fields are not only a global food source but also a major biogenic source of atmospheric methane. Using metatranscriptomics, we comparatively explored structural and functional succession of paddy soil microbiomes in the oxic surface layer and anoxic bulk soil. Cyanobacteria, Fungi, Xanthomonadales, Myxococcales, and Methylococcales were the most abundant and metabolically active groups in the oxic zone, while Clostridia, Actinobacteria, Geobacter, Anaeromyxobacter, Anaerolineae, and methanogenic archaea dominated the anoxic zone. The protein synthesis potential of these groups was about 75% and 50% of the entire community capacity, respectively. Their structure-function relationships in microbiome succession were revealed by classifying the protein-coding transcripts into core, non-core, and taxon-specific transcripts based on homologous gene distribution. The differential expression of core transcripts between the two microbiomes indicated that structural succession is primarily governed by the cellular ability to adapt to the given oxygen condition, involving oxidative stress, nitrogen/phosphorus metabolism, and fermentation. By contrast, the non-core transcripts were expressed from genes involved in the metabolism of various carbon sources. Among those, taxon-specific transcripts revealed highly specialized roles of the dominant groups in community-wide functioning. For instance, taxon-specific transcripts involved in photosynthesis and methane oxidation were a characteristic of the oxic zone, while those related to methane production and aromatic compound degradation were specific to the anoxic zone. Degradation of organic matters, antibiotics resistance, and secondary metabolite production were detected to be expressed in both the oxic and anoxic zones, but by different taxonomic groups. Cross-feeding of methanol between members of the Methylococcales and Xanthomonadales was suggested by the observation that in the oxic zone, they both exclusively expressed homologous genes encoding methanol dehydrogenase. Our metatranscriptomic analysis suggests that paddy soil microbiomes act as complex, functionally coordinated assemblages whose taxonomic composition is governed by the prevailing habitat factors and their hierarchical importance for community succession.

Microbial community transcriptional networks are conserved in three domains at ocean basin scales

Planktonic microbial communities in the ocean are typically dominated by several cosmopolitan clades of Bacteria, Archaea, and Eukarya characterized by their ribosomal RNA gene phylogenies and genomic features. Although the environments these communities inhabit range from coastal to open ocean waters, how the biological dynamics vary between such disparate habitats is not well known. To gain insight into the differential activities of microbial populations inhabiting different oceanic provinces we compared the daily metatranscriptome profiles of related microbial populations inhabiting surface waters of both a coastal California upwelling region (CC) as well as the oligotrophic North Pacific Subtropical Gyre (NPSG). Transcriptional networks revealed that the dominant photoautotrophic microbes in each environment (Ostreococcus in CC, Prochlorococcus in NPSG) were central determinants of overall community transcriptome dynamics. Furthermore, heterotrophic bacterial clades common to both ecosystems (SAR11, SAR116, SAR86, SAR406, and Roseobacter) displayed conserved, genome-wide inter- and intrataxon transcriptional patterns and diel cycles. Populations of SAR11 and SAR86 clades in particular exhibited tightly coordinated transcriptional patterns in both coastal and pelagic ecosystems, suggesting that specific biological interactions between these groups are widespread in nature. Our results identify common diurnally oscillating behaviors among diverse planktonic microbial species regardless of habitat, suggesting that highly conserved temporally phased biotic interactions are ubiquitous among planktonic microbial communities worldwide.

Evaluation of metatranscriptomic protocols and application to the study of freshwater microbial communities

Metatranscriptomics of environmental samples enables the identification of community activities without a priori knowledge of taxonomic or functional composition. However, several technical challenges associated with the RNA preparation protocols can affect the relative representation of transcripts and data interpretation. Here, seven replicate metatranscriptomes from planktonic freshwater samples (Lake Lanier, USA) were sequenced to evaluate technical and biological reproducibility of different RNA extraction protocols. Organic versus bead-beating extraction showed significant enrichment for low versus high G + C% mRNA populations respectively. The sequencing data were best modelled by a negative binomial distribution to account for the large technical and biological variation observed. Despite the variation, the transcriptional activities of populations that persisted in year-round metagenomes from the same site consistently showed distinct expression patterns, reflecting different ecologic strategies and allowing us to test prevailing models on the contribution of both rare biosphere and abundant members to community activity. For instance, abundant members of the Verrucomicrobia phylum systematically showed low transcriptional activity compared with other abundant taxa. Our results provide a practical guide to the analysis of metatranscriptomes and advance understanding of the activity and ecology of abundant and rare members of temperate freshwater microbial communities.

Towards an integrated understanding of gut microbiota using insects as model systems

Metazoans form symbioses with microorganisms that synthesize essential nutritional compounds and increase their efficiency to digest and absorb nutrients. Despite the growing awareness that microbes within the gut play key roles in metabolism, health and development of metazoans, symbiotic relationships within the gut are far from fully understood. Insects, which generally harbor a lower microbial diversity than vertebrates, have recently emerged as potential model systems to study these interactions. In this review, we give a brief overview of the characteristics of the gut microbiota in insects in terms of low diversity but high variability at intra- and interspecific levels and we investigate some of the ecological and methodological factors that might explain such variability. We then emphasize how studies integrating an array of techniques and disciplines have the potential to provide new understanding of the biology of this micro eco-system.

Metabolic pathways for the whole community

Background

A convergence of high-throughput sequencing and computational power is transforming biology into information science. Despite these technological advances, converting bits and bytes of sequence information into meaningful insights remains a challenging enterprise. Biological systems operate on multiple hierarchical levels from genomes to biomes. Holistic understanding of biological systems requires agile software tools that permit comparative analyses across multiple information levels (DNA, RNA, protein, and metabolites) to identify emergent properties, diagnose system states, or predict responses to environmental change.

Results

Here we adopt the MetaPathways annotation and analysis pipeline and Pathway Tools to construct environmental pathway/genome databases (ePGDBs) that describe microbial community metabolism using MetaCyc, a highly curated database of metabolic pathways and components covering all domains of life. We evaluate Pathway Tools’ performance on three datasets with different complexity and coding potential, including simulated metagenomes, a symbiotic system, and the Hawaii Ocean Time-series. We define accuracy and sensitivity relationships between read length, coverage and pathway recovery and evaluate the impact of taxonomic pruning on ePGDB construction and interpretation. Resulting ePGDBs provide interactive metabolic maps, predict emergent metabolic pathways associated with biosynthesis and energy production and differentiate between genomic potential and phenotypic expression across defined environmental gradients.

Conclusions

This multi-tiered analysis provides the user community with specific operating guidelines, performance metrics and prediction hazards for more reliable ePGDB construction and interpretation. Moreover, it demonstrates the power of Pathway Tools in predicting metabolic interactions in natural and engineered ecosystems.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-619) contains supplementary material, which is available to authorized users.

Land use type significantly affects microbial gene transcription in soil

Soil microorganisms play an essential role in sustaining biogeochemical processes and cycling of nutrients across different land use types. To gain insights into microbial gene transcription in forest and grassland soil, we isolated mRNA from 32 sampling sites. After sequencing of generated complementary DNA (cDNA), a total of 5,824,229 sequences could be further analyzed. We were able to assign nonribosomal cDNA sequences to all three domains of life. A dominance of bacterial sequences, which were affiliated to 25 different phyla, was found. Bacterial groups capable of aromatic compound degradation such as Phenylobacterium and Burkholderia were detected in significantly higher relative abundance in forest soil than in grassland soil. Accordingly, KEGG pathway categories related to degradation of aromatic ring-containing molecules (e.g., benzoate degradation) were identified in high abundance within forest soil-derived metatranscriptomic datasets. The impact of land use type forest on community composition and activity is evidently to a high degree caused by the presence of wood breakdown products. Correspondingly, bacterial groups known to be involved in lignin degradation and containing ligninolytic genes such as Burkholderia, Bradyrhizobium, and Azospirillum exhibited increased transcriptional activity in forest soil. Higher solar radiation in grassland presumably induced increased transcription of photosynthesis-related genes within this land use type. This is in accordance with high abundance of photosynthetic organisms and plant-infecting viruses in grassland.

Bioinformatic approaches reveal metagenomic characterization of soil microbial community

As is well known, soil is a complex ecosystem harboring the most prokaryotic biodiversity on the Earth. In recent years, the advent of high-throughput sequencing techniques has greatly facilitated the progress of soil ecological studies. However, how to effectively understand the underlying biological features of large-scale sequencing data is a new challenge. In the present study, we used 33 publicly available metagenomes from diverse soil sites (i.e. grassland, forest soil, desert, Arctic soil, and mangrove sediment) and integrated some state-of-the-art computational tools to explore the phylogenetic and functional characterizations of the microbial communities in soil. Microbial composition and metabolic potential in soils were comprehensively illustrated at the metagenomic level. A spectrum of metagenomic biomarkers containing 46 taxa and 33 metabolic modules were detected to be significantly differential that could be used as indicators to distinguish at least one of five soil communities. The co-occurrence associations between complex microbial compositions and functions were inferred by network-based approaches. Our results together with the established bioinformatic pipelines should provide a foundation for future research into the relation between soil biodiversity and ecosystem function.

The transcriptome landscape of Prochlorococcus MED4 and the factors for stabilizing the core genome

BACKGROUND

Gene gain and loss frequently occurs in the cyanobacterium Prochlorococcus, a phototroph that numerically dominates tropical and subtropical open oceans. However, little is known about the stabilization of its core genome, which contains approximately 1250 genes, in the context of genome streamlining. Using Prochlorococcus MED4 as a model organism, we investigated the constraints on core genome stabilization using transcriptome profiling.

RESULTS

RNA-Seq technique was used to obtain the transcriptome map of Prochlorococcus MED4, including operons, untranslated regions, non-coding RNAs, and novel genes. Genome-wide expression profiles revealed that three factors contribute to core genome stabilization. First, a negative correlation between gene expression levels and protein evolutionary rates was observed. Highly expressed genes were overrepresented in the core genome but not in the flexible genome. Gene necessity was determined as a second powerful constraint on genome evolution through functional enrichment analysis. Third, quick mRNA turnover may increase corresponding proteins' fidelity among genes that were abundantly expressed. Together, these factors influence core genome stabilization during MED4 genome evolution.

CONCLUSIONS

Gene expression, gene necessity, and mRNA turnover contribute to core genome maintenance during cyanobacterium Prochlorococcus genus evolution.

In situ gene expression profiling of the thermoacidophilic alga Cyanidioschyzon in relation to visible and ultraviolet irradiance

Ultraviolet and high-intensity visible radiation generate reactive intermediates that damage phototrophic microorganisms. In Yellowstone National Park, the thermoacidophilic alga Cyanidioschyzon exhibits an annual seasonal biomass fluctuation referred to as 'mat decline', where algal viability decreases as ultraviolet and visible irradiances increase during summer. We examined the role irradiance might play in mat decline using irradiance filters that uncouple ultraviolet and visible effects along with custom microarrays to study gene expression in situ. Of the 6507 genes, 88% showed no response to ultraviolet or visible, implying that at the biomolecular level, these algae inhabit a chemostat-like environment and is consistent with the near constant aqueous chemistry measured. The remaining genes exhibited expression changes linked to ultraviolet exposure, to increased visible radiation, or to the apparent combined effects of ultraviolet and visible. Expression of DNA repetitive elements was synchronized, being repressed by visible but also influenced by ultraviolet. At highest irradiance levels, these algae reduced transcription of genes encoding functions involved with DNA replication, photosynthesis and cell cycle progression but exhibited an uptick in activities related to repairing DNA damage. This corroborates known physiological responses to ultraviolet and visible radiation, and leads us to provisionally conclude that mat decline is linked to photoinhibition.

Metatranscriptomic analyses of plankton communities inhabiting surface and subpycnocline waters of the Chesapeake Bay during oxic-anoxic-oxic transitions

We used metatranscriptomics to study the gene transcription patterns of microbial plankton (0.2 to 64 μm) at a mesohaline station in the Chesapeake Bay under transitions from oxic to anoxic waters in spring and from anoxic to oxic waters in autumn. Samples were collected from surface (i.e., above pycnocline) waters (3 m) and from waters beneath the pycnocline (16 to 22 m) in both 2010 and 2011. Metatranscriptome profiles based on function and potential phylogeny were different between 2010 and 2011 and strongly variable in 2011. This difference in variability corresponded with a highly variable ratio of eukaryotic to bacterial sequences (0.3 to 5.5), reflecting transient algal blooms in 2011 that were absent in 2010. The similarity between metatranscriptomes changed at a lower rate during the transition from oxic to anoxic waters than after the return to oxic conditions. Transcripts related to photosynthesis and low-affinity cytochrome oxidases were significantly higher in shallow than in deep waters, while in deep water genes involved in anaerobic metabolism, particularly sulfate reduction, succinyl coenzyme A (succinyl-CoA)-to-propionyl-CoA conversion, and menaquinone synthesis, were enriched relative to in shallow waters. Expected transitions in metabolism between oxic and anoxic deep waters were reflected in elevated levels of anaerobic respiratory reductases and utilization of propenediol and acetoin. The percentage of archaeal transcripts increased in both years in late summer (from 0.1 to 4.4% of all transcripts in 2010 and from 0.1 to 6.2% in 2011). Denitrification-related genes were expressed in a predicted pattern during the oxic-anoxic transition. Overall, our data suggest that Chesapeake Bay microbial assemblages express gene suites differently in shallow and deep waters and that differences in deep waters reflect variable redox states.

A novel metatranscriptomic approach to identify gene expression dynamics during extracellular electron transfer

Microbial respiration via extracellular electron transfer (EET) is a ubiquitous reaction that occurs throughout anoxic environments and is a driving force behind global biogeochemical cycling of metals. Here we identify specific EET-active microbes and genes in a diverse biofilm using an innovative approach to analyse the dynamic community-wide response to changing EET rates. We find that the most significant gene expression responses to applied EET stimuli occur in only two microbial groups, Desulfobulbaceae and Desulfuromonadales. Metagenomic analyses reveal high coverage draft genomes of these abundant and active microbes. Our metatranscriptomic results show known and unknown genes that are highly responsive to EET stimuli and associated with our identified draft genomes. This new approach yields a comprehensive image of functional microbes and genes related to EET activity in a diverse community, representing the next step towards unravelling complex microbial roles within a community and how microbes adapt to specific environmental stimuli.

Combined analyses of the ITS loci and the corresponding 16S rRNA genes reveal high micro- and macrodiversity of SAR11 populations in the Red Sea

Bacteria belonging to the SAR11 clade are among the most abundant prokaryotes in the pelagic zone of the ocean. 16S rRNA gene-based analyses indicate that they constitute up to 60% of the bacterioplankton community in the surface waters of the Red Sea. This extremely oligotrophic water body is further characterized by an epipelagic zone, which has a temperature above 24 °C throughout the year, and a remarkable uniform temperature (~22 °C) and salinity (~41 psu) from the mixed layer (~200 m) to the bottom at over 2000 m depth. Despite these conditions that set it apart from other marine environments, the microbiology of this ecosystem is still vastly understudied. Prompted by the limited phylogenetic resolution of the 16S rRNA gene, we extended our previous study by sequencing the internal transcribed spacer (ITS) region of SAR11 in different depths of the Red Sea's water column together with the respective 16S fragment. The overall diversity captured by the ITS loci was ten times higher than that of the corresponding 16S rRNA genes. Moreover, species estimates based on the ITS showed a highly diverse population of SAR11 in the mixed layer that became diminished in deep isothermal waters, which was in contrast to results of the related 16S rRNA genes. While the 16S rRNA gene-based sequences clustered into three phylogenetic subgroups, the related ITS fragments fell into several phylotypes that showed clear depth-dependent shifts in relative abundances. Blast-based analyses not only documented the observed vertical partitioning and universal co-occurrence of specific phylotypes in five other distinct oceanic provinces, but also highlighted the influence of ecosystem-specific traits (e.g., temperature, nutrient availability, and concentration of dissolved oxygen) on the population dynamics of this ubiquitous marine bacterium.

Unifying concepts in anaerobic respiration: insights from dissimilatory sulfur metabolism

Behind the versatile nature of prokaryotic energy metabolism is a set of redox proteins having a highly modular character. It has become increasingly recognized that a limited number of redox modules or building blocks appear grouped in different arrangements, giving rise to different proteins and functionalities. This modularity most likely reveals a common and ancient origin for these redox modules, and is obviously reflected in similar energy conservation mechanisms. The dissimilation of sulfur compounds was probably one of the earliest biological strategies used by primitive organisms to obtain energy. Here, we review some of the redox proteins involved in dissimilatory sulfur metabolism, focusing on sulfate reducing organisms, and highlight links between these proteins and others involved in different processes of anaerobic respiration. Noteworthy are links to the complex iron-sulfur molybdoenzyme family, and heterodisulfide reductases of methanogenic archaea. We discuss how chemiosmotic and electron bifurcation/confurcation may be involved in energy conservation during sulfate reduction, and how introduction of an additional module, multiheme cytochromes c, opens an alternative bioenergetic strategy that seems to increase metabolic versatility. Finally, we highlight new families of heterodisulfide reductase-related proteins from non-methanogenic organisms, which indicate a widespread distribution for these protein modules and may indicate a more general involvement of thiol/disulfide conversions in energy metabolism. This article is part of a Special Issue entitled: The evolutionary aspects of bioenergetic systems.

Expression patterns reveal niche diversification in a marine microbial assemblage

Resolving the ecological niches of coexisting marine microbial taxa is challenging due to the high species richness of microbial communities and the apparent functional redundancy in bacterial genomes and metagenomes. Here, we generated over 11 million Illumina reads of protein-encoding transcripts collected from well-mixed southeastern US coastal waters to characterize gene expression patterns distinguishing the ecological roles of hundreds of microbial taxa sharing the same environment. The taxa with highest in situ growth rates (based on relative abundance of ribosomal protein transcripts) were typically not the greatest contributors to community transcription, suggesting strong top-down ecological control, and their diverse transcriptomes indicated roles as metabolic generalists. The taxa with low in situ growth rates typically had low diversity transcriptomes dominated by specialized metabolisms. By identifying protein-encoding genes with atypically high expression for their level of conservation, unique functional roles of community members emerged related to substrate use (such as complex carbohydrates, fatty acids, methanesulfonate, taurine, tartrate, ectoine), alternative energy-conservation strategies (proteorhodopsin, AAnP, V-type pyrophosphatases, sulfur oxidation, hydrogen oxidation) and mechanisms for negotiating a heterogeneous environment (flagellar motility, gliding motility, adhesion strategies). On average, the heterotrophic bacterioplankton dedicated 7% of their transcriptomes to obtaining energy by non-heterotrophic means. This deep sequencing of a coastal bacterioplankton transcriptome provides the most highly resolved view of bacterioplankton niche dimensions yet available, uncovering a spectrum of unrecognized ecological strategies.

Cytoplasmic sulfurtransferases in the purple sulfur bacterium Allochromatium vinosum: evidence for sulfur transfer from DsrEFH to DsrC

While the importance of sulfur transfer reactions is well established for a number of biosynthetic pathways, evidence has only started to emerge that sulfurtransferases may also be major players in sulfur-based microbial energy metabolism. Among the first organisms studied in this regard is the phototrophic purple sulfur bacterium Allochromatium vinosum. During the oxidation of reduced sulfur species to sulfate this Gammaproteobacterium accumulates sulfur globules. Low molecular weight organic persulfides have been proposed as carrier molecules transferring sulfur from the periplasmic sulfur globules into the cytoplasm where it is further oxidized via the “Dsr” (dissimilatory sulfite reductase) proteins. We have suggested earlier that the heterohexameric protein DsrEFH is the direct or indirect acceptor for persulfidic sulfur imported into the cytoplasm. This proposal originated from the structural similarity of DsrEFH with the established sulfurtransferase TusBCD from E. coli. As part of a system for tRNA modification TusBCD transfers sulfur to TusE, a homolog of another crucial component of the A. vinosum Dsr system, namely DsrC. Here we show that neither DsrEFH nor DsrC have the ability to mobilize sulfane sulfur directly from low molecular weight thiols like thiosulfate or glutathione persulfide. However, we demonstrate that DsrEFH binds sulfur specifically to the conserved cysteine residue DsrE-Cys78 in vitro. Sulfur atoms bound to cysteines in DsrH and DsrF were not detected. DsrC was exclusively persulfurated at DsrC-Cys111 in the penultimate position of the protein. Most importantly, we show that persulfurated DsrEFH indeed serves as an effective sulfur donor for DsrC in vitro. The active site cysteines Cys78 of DsrE and Cys20 of DsrH furthermore proved to be essential for sulfur oxidation in vivo supporting the notion that DsrEFH and DsrC are part of a sulfur relay system that transfers sulfur from a persulfurated carrier molecule to the dissimilatory sulfite reductase DsrAB.

Comparative metatranscriptomics reveals widespread community responses during phenanthrene degradation in soil

Soil microbial community response to phenanthrene was evaluated by metatranscriptomics. A marked increase in transcripts involved in aromatic compound metabolism, respiration and stress responses, and concurrent decreases in virulence, carbohydrate, DNA metabolism and phosphorus metabolism transcripts was revealed. Phenanthrene addition led to a 1.8-fold to 33-fold increase in the abundance of dioxygenase, stress response and detoxification transcripts, whereas those of general metabolism were little affected. Heavy metal P-type ATPases and thioredoxin transcripts were more abundant in the phenanthrene-amended soil, and this is the first time these proteins have been associated with polycyclic aromatic hydrocarbon (PAH) stress in microorganisms. Annotation with custom databases constructed with bacterial or fungal PAH metabolism protein sequences showed that increases in PAH-degradatory gene expression occurred for all gene groups investigated. Taxonomic determination of mRNA transcripts showed widespread changes in the bacteria, archaea and fungi, and the actinobacteria were responsible for most of the de novo expression of transcripts associated with dioxygenases, stress response and detoxification genes. This is the first report of an experimental metatranscriptomic study detailing microbial community responses to a pollutant in soil, and offers information on novel in situ effects of PAHs on soil microbes that can be explored further.

Generation and analysis of a mouse intestinal metatranscriptome through Illumina based RNA-sequencing

With the advent of high through-put sequencing (HTS), the emerging science of metagenomics is transforming our understanding of the relationships of microbial communities with their environments. While metagenomics aims to catalogue the genes present in a sample through assessing which genes are actively expressed, metatranscriptomics can provide a mechanistic understanding of community inter-relationships. To achieve these goals, several challenges need to be addressed from sample preparation to sequence processing, statistical analysis and functional annotation. Here we use an inbred non-obese diabetic (NOD) mouse model in which germ-free animals were colonized with a defined mixture of eight commensal bacteria, to explore methods of RNA extraction and to develop a pipeline for the generation and analysis of metatranscriptomic data. Applying the Illumina HTS platform, we sequenced 12 NOD cecal samples prepared using multiple RNA-extraction protocols. The absence of a complete set of reference genomes necessitated a peptide-based search strategy. Up to 16% of sequence reads could be matched to a known bacterial gene. Phylogenetic analysis of the mapped ORFs revealed a distribution consistent with ribosomal RNA, the majority from Bacteroides or Clostridium species. To place these HTS data within a systems context, we mapped the relative abundance of corresponding Escherichia coli homologs onto metabolic and protein-protein interaction networks. These maps identified bacterial processes with components that were well-represented in the datasets. In summary this study highlights the potential of exploiting the economy of HTS platforms for metatranscriptomics.

Dissimilatory sulfur cycling in oxygen minimum zones: an emerging metagenomics perspective

Biological diversity in marine OMZs (oxygen minimum zones) is dominated by a complex community of bacteria and archaea whose anaerobic metabolisms mediate key steps in global nitrogen and carbon cycles. Molecular and physiological studies now confirm that OMZs also support diverse micro-organisms capable of utilizing inorganic sulfur compounds for energy metabolism. The present review focuses specifically on recent metagenomic data that have helped to identify the molecular basis for autotrophic sulfur oxidation with nitrate in the OMZ water column, as well as a cryptic role for heterotrophic sulfate reduction. Interpreted alongside marker gene surveys and process rate measurements, these data suggest an active sulfur cycle with potentially substantial roles in organic carbon input and mineralization and critical links to the OMZ nitrogen cycle. Furthermore, these studies have created a framework for comparing the genomic diversity and ecology of pelagic sulfur-metabolizing communities from diverse low-oxygen regions.

Combined phylogenetic and genomic approaches for the high-throughput study of microbial habitat adaptation

High-throughput sequencing technologies provide new opportunities to address longstanding questions about habitat adaptation in microbial organisms. How have microbes managed to adapt to such a wide range of environments, and what genomic features allow for such adaptation? We review recent large-scale studies of habitat adaptation, with emphasis on those that utilize phylogenetic techniques. On the basis of current trends, we summarize methodological challenges faced by investigators, and the tools, techniques and analytical approaches available to overcome them. Phylogenetic approaches and detailed information about each environmental sample will be crucial as the ability to collect genome sequences continues to expand.

Metatranscriptomic analysis of sulfur oxidation genes in the endosymbiont of solemya velum

Thioautotrophic endosymbionts in the Domain Bacteria mediate key sulfur transformations in marine reducing environments. However, the molecular pathways underlying symbiont metabolism and the extent to which these pathways are expressed in situ are poorly characterized for almost all symbioses. This is largely due to the difficulty of culturing symbionts apart from their hosts. Here, we use pyrosequencing of community RNA transcripts (i.e., the metatranscriptome) to characterize enzymes of dissimilatory sulfur metabolism in the model symbiosis between the coastal bivalve Solemya velum and its intracellular thioautotrophic symbionts. High-throughput sequencing of total RNA from the symbiont-containing gill of a single host individual generated 1.6 million sequence reads (500 Mbp). Of these, 43,735 matched Bacteria protein-coding genes in BLASTX searches of the NCBI database. The taxonomic identities of the matched genes indicated relatedness to diverse species of sulfur-oxidizing Gammaproteobacteria, including other thioautotrophic symbionts and the purple sulfur bacterium Allochromatium vinosum. Manual querying of these data identified 28 genes from diverse pathways of sulfur energy metabolism, including the dissimilatory sulfite reductase (Dsr) pathway for sulfur oxidation to sulfite, the APS pathway for sulfite oxidation, and the Sox pathway for thiosulfate oxidation. In total, reads matching sulfur energy metabolism genes represented 7% of the Bacteria mRNA pool. Together, these data highlight the dominance of thioautotrophy in the context of symbiont community metabolism, identify the likely pathways mediating sulfur oxidation, and illustrate the utility of metatranscriptome sequencing for characterizing community gene transcription of uncultured symbionts.