Spatial patterns and light-driven variation of microbial population gene expression in surface waters of the oligotrophic open ocean

A Reduction of Transcriptional Regulation in Aquatic Oligotrophic Microorganisms Enhances Fitness in Nutrient-Poor Environments

In this review, we consider the regulatory strategies of aquatic oligotrophs, microbial cells that are adapted to thrive under low-nutrient concentrations in oceans, lakes, and other aquatic ecosystems. Many reports have concluded that oligotrophs use less transcriptional regulation than copiotrophic cells, which are adapted to high nutrient concentrations and are far more common subjects for laboratory investigations of regulation. It is theorized that oligotrophs have retained alternate mechanisms of regulation, such as riboswitches, that provide shorter response times and smaller amplitude responses and require fewer cellular resources. We examine the accumulated evidence for distinctive regulatory strategies in oligotrophs. We explore differences in the selective pressures copiotrophs and oligotrophs encounter and ask why, although evolutionary history gives copiotrophs and oligotrophs access to the same regulatory mechanisms, they might exhibit distinctly different patterns in how these mechanisms are used. We discuss the implications of these findings for understanding broad patterns in the evolution of microbial regulatory networks and their relationships to environmental niche and life history strategy. We ask whether these observations, which have emerged from a decade of increased investigation of the cell biology of oligotrophs, might be relevant to recent discoveries of many microbial cell lineages in nature that share with oligotrophs the property of reduced genome size.

The microbial contribution to reactive oxygen species dynamics in marine ecosystems

This review surveys the current state of knowledge of the concentrations, sources and sinks of reactive oxygen species (ROS) in the ocean. Both abiotic and biotic factors contribute to ROS dynamics in seawater, and ROS can feature prominently in marine microbe-microbe interactions. The sun plays a key role in the production of ROS in the ocean, and consequently ROS concentrations are typically maximal in the sun-exposed surface. However, microbes can also contribute significantly to extracellular ROS. Production of superoxide is widespread within the microbial community, and may benefit the producers as antimicrobial agents or perhaps more generally, as a means of nutrient scavenging. Decomposition of hydrogen peroxide is a community-wide activity, though some members may play less significant roles in this process. The more reactive forms of ROS, singlet oxygen and the hydroxyl radical, may be less important as microbial stressors, as they tend to react with the chemicals in seawater before they can contact the cells. However, exceptions may exist for microbes attached to singlet oxygen-generating sinking particulate matter. Extracellular ROS thus plays an important role in the ecology of marine microbes, the full extent to which we are only beginning to appreciate.

Cross-protection from hydrogen peroxide by helper microbes: the impacts on the cyanobacterium Prochlorococcus and other beneficiaries in marine communities

Hydrogen peroxide (HOOH) is a reactive oxygen species, derived from molecular oxygen, that is capable of damaging microbial cells. Surprisingly, the HOOH defence systems of some aerobes in the oxygenated marine environments are critically depleted, relative to model aerobes. For instance, the gene encoding catalase is absent in the numerically dominant photosynthetic cyanobacterium, Prochlorococcus. Accordingly, Prochlorococcus is highly susceptible to HOOH when exposed as pure cultures. Pure cultures do not exist in the marine environment, however. Catalase-positive community members can remove HOOH from the seawater medium, thus lowering the threat to Prochlorococcus and any other member that likewise lacks their own catalase. This cross-protection may constitute a loosely defined symbiosis, whereby the catalase-positive helper cells may benefit through the acquisition of nutrients released by the beneficiaries such as Prochlorococcus. Other members of the community that may be helped by the catalase-positive cells may include some lineages of Synechococcus - the sister genus of Prochlorococcus - as well as some lineages of SAR11 and ammonia oxidizing archaea and bacteria. The co-occurrence of catalase-positive and -negative members suggests that cross-protection from HOOH-mediated oxidative stress may play an important role in the construction of the marine microbial community.

Bacterial Diversity and Nitrogen Utilization Strategies in the Upper Layer of the Northwestern Pacific Ocean

Nitrogen (N) is a primary limiting nutrient for bacterial growth and productivity in the ocean. To better understand bacterial community and their N utilization strategy in different N regimes of the ocean, we examined bacterial diversity, diazotrophic diversity, and N utilization gene expressions in the northwestern Pacific Ocean (NWPO) using a combination of high-throughput sequencing and real-time qPCR methods. 521 and 204 different operational taxonomic units (OTUs) were identified in the 16s rRNA and nifH libraries from nine surface samples. Of the 16s rRNA gene OTUs, 11.9% were observed in all samples while 3.5 and 15.9% were detected only in N-sufficient and N-deficient samples. Proteobacteria, Cyanobacteria and Bacteroidetes dominated the bacterial community. Prochlorococcus and Pseudoalteromonas were the most abundant at the genus level in N-deficient regimes, while SAR86, Synechococcus and SAR92 were predominant in the Kuroshio-Oyashio confluence region. The distribution of the nifH gene presented great divergence among sampling stations: Cyanobacterium_UCYN-A dominated the N-deficient stations, while clusters related to the Alpha-, Beta-, and Gamma-Proteobacteria were abundant in other stations. Temperature was the main factor that determined bacterial community structure and diversity while concentration of NO_X-N was significantly correlated with structure and distribution of N₂-fixing microorganisms. Expression of the ammonium transporter was much higher than that of urea transporter subunit A (urtA) and ferredoxin-nitrate reductase, while urtA had an increased expression in N-deficient surface water. The predicted ammonium transporter and ammonium assimilation enzymes were most abundant in surface samples while urease and nitrogenase were more abundant in the N-deficient regions. These findings underscore the fact that marine bacteria have evolved diverse N utilization strategies to adapt to different N habitats, and that urea metabolism is of vital ecological importance in N-deficient regimes.

mdRNA-Seq analysis of marine microbial communities from the northern Red Sea

Metatranscriptomic differential RNA-Seq (mdRNA-Seq) identifies the suite of active transcriptional start sites at single-nucleotide resolution through enrichment of primary transcript 5′ ends. Here we analyzed the microbial community at 45 m depth at Station A in the northern Gulf of Aqaba, Red Sea, during 500 m deep mixing in February 2012 using mdRNA-Seq and a parallel classical RNA-Seq approach. We identified promoters active in situ for five different pico-planktonic genera (the SAR11 clade of Alphaproteobacteria, Synechococcus of Cyanobacteria, Euryarchaeota, Thaumarchaeota, and Micromonas as an example for picoeukaryotic algae), showing the applicability of this approach to highly diverse microbial communities. 16S rDNA quantification revealed that 24% of the analyzed community were group II marine Euryarchaeota in which we identified a highly abundant non-coding RNA, Tan1, and detected very high expression of genes encoding intrinsically disordered proteins, as well as enzymes for the synthesis of specific B vitamins, extracellular peptidases, carbohydrate-active enzymes, and transport systems. These results highlight previously unknown functions of Euryarchaeota with community-wide relevance. The complementation of metatranscriptomic studies with mdRNA-Seq provides substantial additional information regarding transcriptional start sites, promoter activities, and the identification of non-coding RNAs.

Genetic Diversity Affects the Daily Transcriptional Oscillations of Marine Microbial Populations

Marine microbial communities are genetically diverse but have robust synchronized daily transcriptional patterns at the genus level that are similar across a wide variety of oceanic regions. We developed a microarray-inspired gene-centric approach to resolve transcription of closely-related but distinct strains/ecotypes in high-throughput sequence data. Applying this approach to the existing metatranscriptomics datasets collected from two different oceanic regions, we found unique and variable patterns of transcription by individual taxa within the abundant picocyanobacteria Prochlorococcus and Synechococcus, the alpha Proteobacterium Pelagibacter and the eukaryotic picophytoplankton Ostreococcus. The results demonstrate that marine microbial taxa respond differentially to variability in space and time in the ocean. These intra-genus individual transcriptional patterns underlie whole microbial community responses, and the approach developed here facilitates deeper insights into microbial population dynamics.

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.

Marine microbial metagenomics: from individual to the environment

Microbes are the most abundant biological entities on earth, therefore, studying them is important for understanding their roles in global ecology. The science of metagenomics is a relatively young field of research that has enjoyed significant effort since its inception in 1998. Studies using next-generation sequencing techniques on single genomes and collections of genomes have not only led to novel insights into microbial genomics, but also revealed a close association between environmental niches and genome evolution. Herein, we review studies investigating microbial genomics (largely in the marine ecosystem) at the individual and community levels to summarize our current understanding of microbial ecology in the environment.

Metaproteomic survey of six aquatic habitats: discovering the identities of microbial populations active in biogeochemical cycling

Our goal is to strengthen the foundations of metaproteomics as a microbial community analysis tool that links the functional identity of actively expressed gene products with host phylogeny. We used shotgun metaproteomics to survey waters in six disparate aquatic habitats (Cayuga Lake, NY; Oneida Lake, NY; Gulf of Maine; Chesapeake Bay, MD; Gulf of Mexico; and the South Pacific). Peptide pools prepared from filter-gathered microbial biomass, analyzed by nano-liquid chromatography-mass spectrometry (MS/MS) generating 9,693 ± 1,073 mass spectra identified 326 ± 107 bacterial proteins per sample. Distribution of proteobacterial (Alpha and Beta) and cyanobacterial (Prochlorococcus and Synechococcus spp.) protein hosts across all six samples was consistent with the previously published biogeography for these microorganisms. Marine samples were enriched in transport proteins (TRAP-type for dicarboxylates and ATP binding cassette (ABC)-type for amino acids and carbohydrates) compared with the freshwater samples. We were able to match in situ expression of many key proteins catalyzing C-, N-, and S-cycle processes with their bacterial hosts across all six habitats. Pelagibacter was identified as the host of ABC-type sugar-, organic polyanion-, and glycine betaine-transport proteins; this extends previously published studies of Pelagibacter's in situ biogeochemical role in marine C- and N-metabolism. Proteins matched to Ruegeria confirmed these organism's role in marine waters oxidizing both carbon monoxide and sulfide. By documenting both processes expressed in situ and the identity of host cells, metaproteomics tested several existing hypotheses about ecophysiological processes and provided fodder for new ones.

A microarray for assessing transcription from pelagic marine microbial taxa

Metagenomic approaches have revealed unprecedented genetic diversity within microbial communities across vast expanses of the world's oceans. Linking this genetic diversity with key metabolic and cellular activities of microbial assemblages is a fundamental challenge. Here we report on a collaborative effort to design MicroTOOLs (Microbiological Targets for Ocean Observing Laboratories), a high-density oligonucleotide microarray that targets functional genes of diverse taxa in pelagic and coastal marine microbial communities. MicroTOOLs integrates nucleotide sequence information from disparate data types: genomes, PCR-amplicons, metagenomes, and metatranscriptomes. It targets 19 400 unique sequences over 145 different genes that are relevant to stress responses and microbial metabolism across the three domains of life and viruses. MicroTOOLs was used in a proof-of-concept experiment that compared the functional responses of microbial communities following Fe and P enrichments of surface water samples from the North Pacific Subtropical Gyre. We detected transcription of 68% of the gene targets across major taxonomic groups, and the pattern of transcription indicated relief from Fe limitation and transition to N limitation in some taxa. Prochlorococcus (eHLI), Synechococcus (sub-cluster 5.3) and Alphaproteobacteria SAR11 clade (HIMB59) showed the strongest responses to the Fe enrichment. In addition, members of uncharacterized lineages also responded. The MicroTOOLs microarray provides a robust tool for comprehensive characterization of major functional groups of microbes in the open ocean, and the design can be easily amended for specific environments and research questions.

Integrated metagenomic and metaproteomic analyses of marine biofilm communities

Metagenomic and metaproteomic analyses were utilized to determine the composition and function of complex air-water interface biofilms sampled from the hulls of two US Navy destroyers. Prokaryotic community analyses using PhyloChip-based 16S rDNA profiling revealed two significantly different and taxonomically rich biofilm communities (6,942 taxa) in which the majority of unique taxa were ascribed to members of the Gammaproteobacteria, Alphaproteobacteria and Clostridia. Although metagenomic sequencing indicated that both biofilms were dominated by prokaryotic sequence reads (> 91%) with the majority of the bacterial reads belonging to the Alphaproteobacteria, the Ship-1 metagenome harbored greater organismal and functional diversity and was comparatively enriched for sequences from Cyanobacteria, Bacteroidetes and macroscopic eukaryotes, whereas the Ship-2 metagenome was enriched for sequences from Proteobacteria and microscopic photosynthetic eukaryotes. Qualitative liquid chromatography-tandem mass spectrometry metaproteome analyses identified 678 unique proteins, revealed little overlap in species and protein composition between the ships and contrasted with the metagenomic data in that ~80% of classified and annotated proteins were of eukaryotic origin and dominated by members of the Bacillariophyta, Cnidaria, Chordata and Arthropoda (data deposited to the ProteomeXchange, identifier PXD000961). Within the shared metaproteome, quantitative (18)O and iTRAQ analyses demonstrated a significantly greater abundance of structural proteins from macroscopic eukaryotes on Ship-1 and diatom photosynthesis proteins on Ship-2. Photosynthetic pigment composition and elemental analyses confirmed that both biofilms were dominated by phototrophic processes. These data begin to provide a better understanding of the complex organismal and biomolecular composition of marine biofilms while highlighting caveats in the interpretation of stand-alone environmental '-omics' datasets.

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.

Microbial and viral metagenomes of a subtropical freshwater reservoir subject to climatic disturbances

Extreme climatic activities, such as typhoons, are widely known to disrupt our natural environment. In particular, studies have revealed that typhoon-induced perturbations can result in several long-term effects on various ecosystems. In this study, we have conducted a 2-year metagenomic survey to investigate the microbial and viral community dynamics associated with environmental changes and seasonal variations in an enclosed freshwater reservoir subject to episodic typhoons. We found that the microbial community structure and the associated metagenomes continuously changed, where microbial richness increased after typhoon events and decreased during winter. Among the environmental factors that influenced changes in the microbial community, precipitation was considered to be the most significant. Similarly, the viral community regularly showed higher relative abundances and diversity during summer in comparison to winter, with major variations happening in several viral families including Siphoviridae, Myoviridae, Podoviridae and Microviridae. Interestingly, we also found that the precipitation level was associated with the terrestrial viral abundance in the reservoir. In contrast to the dynamic microbial community (L-divergence 0.73 ± 0.25), we found that microbial metabolic profiles were relatively less divergent (L-divergence 0.24 ± 0.04) at the finest metabolic resolution. This study provides for the first time a glimpse at the microbial and viral community dynamics of a subtropical freshwater ecosystem, adding a comprehensive set of new knowledge to aquatic environments.

Pattern and synchrony of gene expression among sympatric marine microbial populations

Planktonic marine microbes live in dynamic habitats that demand rapid sensing and response to periodic as well as stochastic environmental change. The kinetics, regularity, and specificity of microbial responses in situ, however, are not well-described. We report here simultaneous multitaxon genome-wide transcriptome profiling in a naturally occurring picoplankton community. An in situ robotic sampler using a Lagrangian sampling strategy enabled continuous tracking and repeated sampling of coherent microbial populations over 2 d. Subsequent RNA sequencing analyses yielded genome-wide transcriptome profiles of eukaryotic (Ostreococcus) and bacterial (Synechococcus) photosynthetic picoplankton as well as proteorhodopsin-containing heterotrophs, including Pelagibacter, SAR86-cluster Gammaproteobacteria, and marine Euryarchaea. The photosynthetic picoplankton exhibited strong diel rhythms over thousands of gene transcripts that were remarkably consistent with diel cycling observed in laboratory pure cultures. In contrast, the heterotrophs did not cycle diurnally. Instead, heterotrophic picoplankton populations exhibited cross-species synchronous, tightly regulated, temporally variable patterns of gene expression for many genes, particularly those genes associated with growth and nutrient acquisition. This multitaxon, population-wide gene regulation seemed to reflect sporadic, short-term, reversible responses to high-frequency environmental variability. Although the timing of the environmental responses among different heterotrophic species seemed synchronous, the specific metabolic genes that were expressed varied from taxon to taxon. In aggregate, these results provide insights into the kinetics, diversity, and functional patterns of microbial community response to environmental change. Our results also suggest a means by which complex multispecies metabolic processes could be coordinated, facilitating the regulation of matter and energy processing in a dynamically changing environment.

Metagenomics and metatranscriptomics: windows on CF-associated viral and microbial communities

BACKGROUND

Samples collected from CF patient airways often contain large amounts of host-derived nucleic acids that interfere with recovery and purification of microbial and viral nucleic acids. This study describes metagenomic and metatranscriptomic methods that address these issues.

METHODS

Microbial and viral metagenomes, and microbial metatranscriptomes, were successfully prepared from sputum samples from five adult CF patients.

RESULTS

Contaminating host DNA was dramatically reduced in the metagenomes. Each CF patient presented a unique microbiome; in some Pseudomonas aeruginosa was replaced by other opportunistic bacteria. Even though the taxonomic composition of the microbiomes is very different, the metabolic potentials encoded by the community are very similar. The viral communities were dominated by phages that infect major CF pathogens. The metatranscriptomes reveal differential expression of encoded metabolic potential with changing health status.

CONCLUSIONS

Microbial and viral metagenomics combined with microbial transcriptomics characterize the dynamic polymicrobial communities found in CF airways, revealing both the taxa present and their current metabolic activities. These approaches can facilitate the development of individualized treatment plans and novel therapeutic approaches.

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.

Metagenomic and metabolic profiling of nonlithifying and lithifying stromatolitic mats of Highborne Cay, The Bahamas

BACKGROUND

Stromatolites are laminated carbonate build-ups formed by the metabolic activity of microbial mats and represent one of the oldest known ecosystems on Earth. In this study, we examined a living stromatolite located within the Exuma Sound, The Bahamas and profiled the metagenome and metabolic potential underlying these complex microbial communities.

METHODOLOGY/PRINCIPAL FINDINGS

The metagenomes of the two dominant stromatolitic mat types, a nonlithifying (Type 1) and lithifying (Type 3) microbial mat, were partially sequenced and compared. This deep-sequencing approach was complemented by profiling the substrate utilization patterns of the mats using metabolic microarrays. Taxonomic assessment of the protein-encoding genes confirmed previous SSU rRNA analyses that bacteria dominate the metagenome of both mat types. Eukaryotes comprised less than 13% of the metagenomes and were rich in sequences associated with nematodes and heterotrophic protists. Comparative genomic analyses of the functional genes revealed extensive similarities in most of the subsystems between the nonlithifying and lithifying mat types. The one exception was an increase in the relative abundance of certain genes associated with carbohydrate metabolism in the lithifying Type 3 mats. Specifically, genes associated with the degradation of carbohydrates commonly found in exopolymeric substances, such as hexoses, deoxy- and acidic sugars were found. The genetic differences in carbohydrate metabolisms between the two mat types were confirmed using metabolic microarrays. Lithifying mats had a significant increase in diversity and utilization of carbon, nitrogen, phosphorus and sulfur substrates.

CONCLUSION/SIGNIFICANCE

The two stromatolitic mat types retained similar microbial communities, functional diversity and many genetic components within their metagenomes. However, there were major differences detected in the activity and genetic pathways of organic carbon utilization. These differences provide a strong link between the metagenome and the physiology of the mats, as well as new insights into the biological processes associated with carbonate precipitation in modern marine stromatolites.

Description of viral assemblages associated with the Gorgonia ventalina holobiont

The diversity and function of viruses in coral holobionts has only recently received attention. The non-reef building gorgonian octocoral, Gorgonia ventalina, is a major constituent of Caribbean reefs. We investigated viral communities associated with G. ventalina tissues to understand their role in gorgonian ecology. Pyrosequencing was used to prepare a total of 514,632 sequence reads of DNA- and RNA-based mixed-community viral genomes (metaviromes). RNA viral assemblages were comprised of primarily unidentifiable reads, with most matching host transcripts and other RNA metaviromes. DNA metaviromes were similar between healthy and diseased tissues and comprised of contiguous sequences (contigs) that matched primarily metazoan and bacterial proteins. Only ~5% of contigs matched viral proteins that were primarily cyanophage and viruses of Chlorella and Ostreococcus. Our results confirm that DNA and RNA viruses comprise a component of the gorgonian holobiont, suggesting that they may play a role in the ecology of G. ventalina.

Transcriptome fingerprinting analysis: an approach to explore gene expression patterns in marine microbial communities

Microbial transcriptomics are providing new insights into the functional processes of microbial communities. However, analysis of each sample is still expensive and time consuming. A rapid and low cost method that would allow the identification of the most interesting samples for posterior in-depth metatranscriptomics analysis would be extremely useful. Here we present Transcriptome Fingerprinting Analysis (TFA) as an approach to fulfill this objective in microbial ecology studies. We have adapted the differential display technique for mRNA fingerprinting based on the PCR amplification of expressed transcripts to interrogate natural microbial eukaryotic communities. Unlike other techniques, TFA does not require prior knowledge of the mRNA sequences to be detected. We have used a set of arbitrary primers coupled with a fluorescence labeled primer targeting the poly(A) tail of the eukaryotic mRNA, with further detection of the resulting labeled cDNA products in an automated genetic analyzer. The output represented by electropherogram peak patterns allowed the comparison of a set of genes expressed at the time of sampling. TFA has been optimized by testing the sensitivity of the method for different initial RNA amounts, and the repeatability of the gene expression patterns with increasing time after sampling both with cultures and environmental samples. Results show that TFA is a promising approach to explore the dynamics of gene expression patterns in microbial communities.

Metatranscriptomic analysis of autonomously collected and preserved marine bacterioplankton

Planktonic microbial activity and community structure is dynamic, and can change dramatically on time scales of hours to days. Yet for logistical reasons, this temporal scale is typically under-sampled in the marine environment. In order to facilitate higher-resolution, long-term observation of microbial diversity and activity, we developed a protocol for automated collection and fixation of marine microbes using the Environmental Sample Processor (ESP) platform. The protocol applies a preservative (RNALater) to cells collected on filters, for long-term storage and preservation of total cellular RNA. Microbial samples preserved using this protocol yielded high-quality RNA after 30 days of storage at room temperature, or onboard the ESP at in situ temperatures. Pyrosequencing of complementary DNA libraries generated from ESP-collected and preserved samples yielded transcript abundance profiles nearly indistinguishable from those derived from conventionally treated replicate samples. To demonstrate the utility of the method, we used a moored ESP to remotely and autonomously collect Monterey Bay seawater for metatranscriptomic analysis. Community RNA was extracted and pyrosequenced from samples collected at four time points over the course of a single day. In all four samples, the oxygenic photoautotrophs were predominantly eukaryotic, while the bacterial community was dominated by Polaribacter-like Flavobacteria and a Rhodobacterales bacterium sharing high similarity with Rhodobacterales sp. HTCC2255. However, each time point was associated with distinct species abundance and gene transcript profiles. These laboratory and field tests confirmed that autonomous collection and preservation is a feasible and useful approach for characterizing the expressed genes and environmental responses of marine microbial communities.

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.

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

Background

Combined metagenomic and metatranscriptomic datasets make it possible to study the molecular evolution of diverse microbial species recovered from their native habitats. The link between gene expression level and sequence conservation was examined using shotgun pyrosequencing of microbial community DNA and RNA from diverse marine environments, and from forest soil.

Results

Across all samples, expressed genes with transcripts in the RNA sample were significantly more conserved than non-expressed gene sets relative to best matches in reference databases. This discrepancy, observed for many diverse individual genomes and across entire communities, coincided with a shift in amino acid usage between these gene fractions. Expressed genes trended toward GC-enriched amino acids, consistent with a hypothesis of higher levels of functional constraint in this gene pool. Highly expressed genes were significantly more likely to fall within an orthologous gene set shared between closely related taxa (core genes). However, non-core genes, when expressed above the level of detection, were, on average, significantly more highly expressed than core genes based on transcript abundance normalized to gene abundance. Finally, expressed genes showed broad similarities in function across samples, being relatively enriched in genes of energy metabolism and underrepresented by genes of cell growth.

Conclusions

These patterns support the hypothesis, predicated on studies of model organisms, that gene expression level is a primary correlate of evolutionary rate across diverse microbial taxa from natural environments. Despite their complexity, meta-omic datasets can reveal broad evolutionary patterns across taxonomically, functionally, and environmentally diverse communities.

Microbial metatranscriptomics in a permanent marine oxygen minimum zone

Simultaneous characterization of taxonomic composition, metabolic gene content and gene expression in marine oxygen minimum zones (OMZs) has potential to broaden perspectives on the microbial and biogeochemical dynamics in these environments. Here, we present a metatranscriptomic survey of microbial community metabolism in the Eastern Tropical South Pacific OMZ off northern Chile. Community RNA was sampled in late austral autumn from four depths (50, 85, 110, 200 m) extending across the oxycline and into the upper OMZ. Shotgun pyrosequencing of cDNA yielded 180,000 to 550,000 transcript sequences per depth. Based on functional gene representation, transcriptome samples clustered apart from corresponding metagenome samples from the same depth, highlighting the discrepancies between metabolic potential and actual transcription. BLAST-based characterizations of non-ribosomal RNA sequences revealed a dominance of genes involved with both oxidative (nitrification) and reductive (anammox, denitrification) components of the marine nitrogen cycle. Using annotations of protein-coding genes as proxies for taxonomic affiliation, we observed depth-specific changes in gene expression by key functional taxonomic groups. Notably, transcripts most closely matching the genome of the ammonia-oxidizing archaeon Nitrosopumilus maritimus dominated the transcriptome in the upper three depths, representing one in five protein-coding transcripts at 85 m. In contrast, transcripts matching the anammox bacterium Kuenenia stuttgartiensis dominated at the core of the OMZ (200 m; 1 in 12 protein-coding transcripts). The distribution of N. maritimus-like transcripts paralleled that of transcripts matching ammonia monooxygenase genes, which, despite being represented by both bacterial and archaeal sequences in the community DNA, were dominated (> 99%) by archaeal sequences in the RNA, suggesting a substantial role for archaeal nitrification in the upper OMZ. These data, as well as those describing other key OMZ metabolic processes (e.g. sulfur oxidation), highlight gene-specific expression patterns in the context of the entire community transcriptome, as well as identify key functional groups for taxon-specific genomic profiling.

Integrated metatranscriptomic and metagenomic analyses of stratified microbial assemblages in the open ocean

As part of an ongoing survey of microbial community gene expression in the ocean, we sequenced and compared ∼38 Mbp of community transcriptomes and ∼157 Mbp of community genomes from four bacterioplankton samples, along a defined depth profile at Station ALOHA in North Pacific subtropical gyre (NPSG). Taxonomic analysis suggested that the samples were dominated by three taxa: Prochlorales, Consistiales and Cenarchaeales, which comprised 36-69% and 29-63% of the annotated sequences in the four DNA and four cDNA libraries, respectively. The relative abundance of these taxonomic groups was sometimes very different in the DNA and cDNA libraries, suggesting differential relative transcriptional activities per cell. For example, the 125 m sample genomic library was dominated by Pelagibacter (∼36% of sequence reads), which contributed fewer sequences to the community transcriptome (∼11%). Functional characterization of highly expressed genes suggested taxon-specific contributions to specific biogeochemical processes. Examples included Roseobacter relatives involved in aerobic anoxygenic phototrophy at 75 m, and an unexpected contribution of low abundance Crenarchaea to ammonia oxidation at 125 m. Read recruitment using reference microbial genomes indicated depth-specific partitioning of coexisting microbial populations, highlighted by a transcriptionally active high-light-like Prochlorococcus population in the bottom of the photic zone. Additionally, nutrient-uptake genes dominated Pelagibacter transcripts, with apparent enrichment for certain transporter types (for example, the C4-dicarboxylate transport system) over others (for example, phosphate transporters). In total, the data support the utility of coupled DNA and cDNA analyses for describing taxonomic and functional attributes of microbial communities in their natural habitats.