Prokaryotic Super Program Advisory Committee DOE Joint Genome Institute, Walnut Creek, CA, March 27, 2013

The Prokaryotic Super Program Advisory Committee met on March 27, 2013 for their annual review the Prokaryotic Super Program at the DOE Joint Genome Institute. As is the case with any site visit or program review, the objective is to evaluate progress in meeting organizational objectives, provide feedback to from the user-community and to assist the JGI in formulating plans for the coming year. The advisors want to commend the JGI for its central role in developing new technologies and capabilities, and for catalyzing the formation of new collaborative user communities. Highlights of the post-meeting exchanges among the advisors focused on the importance of programmatic initiatives including:

• GEBA, which serves as a phylogenetic “base-map” on which our knowledge of functional diversity can be layered.

• FEBA, which promises to provide new insights into the physiological capabilities of prokaryotes under highly standardized conditions.

• Single-cell genomics technology, which is seen to significantly enhance our ability to interpret genomic and metagenomic data and broaden the scope of the GEBA program to encompass at least a part of the microbial “dark-matter”.

• IMG, which is seen to play a central role in JGI programs and is viewed as a strategically important asset in the JGI portfolio.

On this latter point, the committee encourages the formation of a strategic relationship between IMG and the Kbase to ensure that the intelligence, deep knowledge and experience captured in the former is not lost. The committee strongly urges the DOE to continue its support for maintaining this critical resource.

Bacterial and archaeal community structure in the surface microlayer of high mountain lakes examined under two atmospheric aerosol loading scenarios

Bacteria and Archaea of the air-water surface microlayer (neuston) and plankton from three high mountain lakes (Limnological Observatory of the Pyrenees, Spain) were analysed by 16S rRNA gene 454 pyrosequencing (V6 region) in two dates with different atmospheric aerosol loading conditions: (1) under a Saharan dust plume driven by southern winds; and (2) under northern winds with oceanic influence. In general, bacterial communities were richer than archaea, with estimated total richness of c. 2500 OTUs for Bacteria and c. 900 OTUs for Archaea equivalent to a sequencing effort of c. 250,000 and c. 20,000 sequences, respectively. The dominant bacterial OTU was affiliated to Actinobacteria. Archaea were one to two orders of magnitude less abundant than bacteria but were more evenly distributed. Apparently, Bacteroidetes and Thaumarchaeota sequences were preferentially found at the neuston, but no consistent pattern in either total microbial abundance or richness was found in any sample. However, we observed more marked changes in microbial relative abundances between neuston and plankton in the dust-influenced scenario. Higher community dissimilarities between neuston and plankton were also found during the Saharan dust episode, and such differences were higher for Bacteria than for Archaea. Nonetheless, relatively few (< 0.05%) of the neuston sequences matched previously identified airborne microorganisms, and none became important in the dates analysed.

Transcriptional changes underlying elemental stoichiometry shifts in a marine heterotrophic bacterium

Marine bacteria drive the biogeochemical processing of oceanic dissolved organic carbon (DOC), a 750-Tg C reservoir that is a critical component of the global C cycle. Catabolism of DOC is thought to be regulated by the biomass composition of heterotrophic bacteria, as cells maintain a C:N:P ratio of ∼50:10:1 during DOC processing. Yet a complicating factor in stoichiometry-based analyses is that bacteria can change the C:N:P ratio of their biomass in response to resource composition. We investigated the physiological mechanisms of resource-driven shifts in biomass stoichiometry in continuous cultures of the marine heterotrophic bacterium Ruegeria pomeroyi (a member of the Roseobacter clade) under four element limitation regimes (C, N, P, and S). Microarray analysis indicated that the bacterium scavenged for alternate sources of the scarce element when cells were C-, N-, or P-limited; reworked the ratios of biomolecules when C- and P- limited; and exerted tighter control over import/export and cytoplasmic pools when N-limited. Under S limitation, a scenario not existing naturally for surface ocean microbes, stress responses dominated transcriptional changes. Resource-driven changes in C:N ratios of up to 2.5-fold and in C:P ratios of up to sixfold were measured in R. pomeroyi biomass. These changes were best explained if the C and P content of the cells was flexible in the face of shifting resources but N content was not, achieved through the net balance of different transcriptional strategies. The cellular-level metabolic trade-offs that govern biomass stoichiometry in R. pomeroyi may have implications for global carbon cycling if extendable to other heterotrophic bacteria. Strong homeostatic responses to N limitation by marine bacteria would intensify competition with autotrophs. Modification of cellular inventories in C- and P-limited heterotrophs would vary the elemental ratio of particulate organic matter sequestered in the deep ocean.

Community analysis of high- and low-nucleic acid-containing bacteria in NW Mediterranean coastal waters using 16S rDNA pyrosequencing

The ecological significance of the marine bacterial populations distinguishable by flow cytometry on the basis of the fluorescence (FL) of their nucleic acid (NA) content and proxies of cell size (such as side scatter, SSC) remains largely unknown. Some studies have suggested that cells with high NA (HNA) content and high SSC (HS) represent the active members of the community, while the low NA (LNA) cells are inactive members of the same phylogenetic groups. But group-specific activity measurements and phylogenetic assignment after cell sorting have suggested this is not be the case, particularly in open-ocean communities. To test the extent to which the different NA subgroups are similar, and consequently the extent to which they likely have similar ecological and biogeochemical roles in the environment, we analysed the phylogenetic composition of three populations after cell sorting [high NA-high SC (HNA-HS), high NA-low SC (HNA-LS), low NA (LNA)] by 454 pyrosequencing in two contrasting periods of the year in NW Mediterranean coastal waters (BBMO, Blanes Bay Microbial Observatory) where these three populations have recurrent seasonal patterns. Statistical analyses showed that summer and winter samples were significantly different and, importantly, the sorted populations within a sample were composed of different taxa. The majority of taxa were associated with one NA fraction only, and the degree of overlap (i.e. OTUs present simultaneously in 2 fractions) between HNA and LNA and between summer and winter communities was very small. Rhodobacterales, SAR116 and Bacteroidetes contributed primarily to the HNA fraction, whereas other groups such as SAR11 and SAR86 contributed largely to the LNA fractions. Gammaproteobacteria other than SAR86 showed less preference for one particular NA fraction. An increase in diversity was observed from the LNA to the HNA-HS fraction for both sample dates. Our results suggest that, in Blanes Bay, flow cytometric signatures of natural communities track their phylogenetic composition.

Environmental, biochemical and genetic drivers of DMSP degradation and DMS production in the Sargasso Sea

Dimethylsulfide (DMS) is a climatically relevant trace gas produced and cycled by the surface ocean food web. Mechanisms driving intraannual variability in DMS production and dimethylsulfoniopropionate (DMSP) degradation in open-ocean, oligotrophic regions were investigated during a 10-month time-series at the Bermuda Atlantic Time-series Study site in the Sargasso Sea. Abundance and transcription of bacterial DMSP degradation genes, DMSP lyase enzyme activity, and DMS and DMSP concentrations, consumption rates and production rates were quantified over time and depth. This interdisciplinary data set was used to test current hypotheses of the role of light and carbon supply in regulating upper-ocean sulfur cycling. Findings supported UV-A-dependent phytoplankton DMS production. Bacterial DMSP degraders may also contribute significantly to DMS production when temperatures are elevated and UV-A dose is moderate, but may favour DMSP demethylation under low UV-A doses. Three groups of bacterial DMSP degraders with distinct intraannual variability were identified and niche differentiation was indicated. The combination of genetic and biochemical data suggest a modified 'bacterial switch' hypothesis where the prevalence of different bacterial DMSP degradation pathways is regulated by a complex set of factors including carbon supply, temperature and UV-A dose.

Structures of dimethylsulfoniopropionate-dependent demethylase from the marine organism Pelagabacter ubique

Dimethylsulfoniopropionate (DMSP) is a ubiquitous algal metabolite and common carbon and sulfur source for marine bacteria. DMSP is a precursor for the climatically active gas dimethylsulfide that is readily oxidized to sulfate, sulfur dioxide, methanesulfonic acid, and other products that act as cloud condensation nuclei. Although the environmental importance of DMSP metabolism has been known for some time, the enzyme responsible for DMSP demethylation by marine bacterioplankton, dimethylsufoniopropionate-dependent demethylase A (DmdA, EC 2.1.1.B5), has only recently been identified and biochemically characterized. In this work, we report the structure for the apoenzyme DmdA from Pelagibacter ubique (2.1 Å), as well as for DmdA co-crystals soaked with substrate DMSP (1.6 Å) or the cofactor tetrahydrofolate (THF) (1.6 Å). Surprisingly, the overall fold of the DmdA is not similar to other enzymes that typically utilize the reduced form of THF and in fact is a triple domain structure similar to what has been observed for the glycine cleavage T protein or sarcosine oxidase. Specifically, while the THF binding fold appears conserved, previous biochemical studies have shown that all enzymes with a similar fold produce 5,10-methylene-THF, while DmdA catalyzes a redox-neutral methyl transfer reaction to produce 5-methyl-THF. On the basis of the findings presented herein and the available biochemical data, we outline a mechanism for a redox-neutral methyl transfer reaction that is novel to this conserved THF binding domain.

Genomic insights into bacterial DMSP transformations

Genomic and functional genomic methods applied to both model organisms and natural communities have rapidly advanced understanding of bacterial dimethylsulfoniopropionate (DMSP) degradation in the ocean. The genes for the two main pathways in bacterial degradation, routing DMSP to distinctly different biogeochemical fates, have recently been identified. The genes dmdA, -B, -C, and -D mediate the demethylation of DMSP and facilitate retention of carbon and sulfur in the marine microbial food web. The genes dddD, -L, -P, -Q, -W, and -Y mediate the cleavage of DMSP to dimethylsulfide (DMS), with important consequences for ocean-atmosphere sulfur flux. In ocean metagenomes, sufficient copies of these genes are present for approximately 60% of surface ocean bacterial cells to directly participate in DMSP degradation. The factors that regulate these two competing pathways remain elusive, but gene transcription analyses of natural bacterioplankton communities are making headway in unraveling the intricacies of bacterial DMSP processing in the ocean.

Metatranscriptomic signature of exogenous polyamine utilization by coastal bacterioplankton

The polyamines putrescine (PUT) and spermidine (SPD) are ubiquitous in seawater, but mechanisms that drive the degradation of these important nitrogen sources by marine bacteria remain unclear. We employed a comparative metatranscriptomics approach to compare gene transcription patterns between coastal bacterioplankton communities with and without amendments of PUT or SPD, in an effort to understand how bacterial communities and their genes shape polyamine biogeochemistry in the ocean. Statistically different transcript categories in the PUT (25 COG groups) and SPD (23 COG groups) samples, relative to controls that received no amendment (CTRL), indicated that genes encoding the cellular translation machinery and the metabolism of organic nitrogen and carbon became enriched in the community transcriptome when polyamine availability increased. Of the three known pathways for bacterial polyamine degradation, only genes in the transamination pathway were enriched in the PUT and SPD libraries, suggesting that this route dominated polyamine degradation. Taxonomic affiliation of significantly enriched diagnostic genes in the PUT and SPD libraries pointed to roseobacter- and SAR11-affiliated bacteria as the predominant taxa driving transformation in this coastal ocean, although other diverse marine bacterioplankton groups (Gammaproteobacteria, Betaproteobacteria, Actinobacteria and Bacteroidetes) also contributed to polyamine-related gene transcription.

Bacterial community transcription patterns during a marine phytoplankton bloom

Bacterioplankton consume a large proportion of photosynthetically fixed carbon in the ocean and control its biogeochemical fate. We used an experimental metatranscriptomics approach to compare bacterial activities that route energy and nutrients during a phytoplankton bloom compared with non-bloom conditions. mRNAs were sequenced from duplicate bloom and control microcosms 1 day after a phytoplankton biomass peak, and transcript copies per litre of seawater were calculated using an internal mRNA standard. Transcriptome analysis revealed a potential novel mechanism for enhanced efficiency during carbon-limited growth, mediated through membrane-bound pyrophosphatases [V-type H(+)-translocating; hppA]; bloom bacterioplankton participated less in this metabolic energy scavenging than non-bloom bacterioplankton, with possible implications for differences in growth yields on organic substrates. Bloom bacterioplankton transcribed more copies of genes predicted to increase cell surface adhesiveness, mediated by changes in bacterial signalling molecules related to biofilm formation and motility; these may be important in microbial aggregate formation. Bloom bacterioplankton also transcribed more copies of genes for organic acid utilization, suggesting an increased importance of this compound class in the bioreactive organic matter released during phytoplankton blooms. Transcription patterns were surprisingly faithful within a taxon regardless of treatment, suggesting that phylogeny broadly predicts the ecological roles of bacterial groups across 'boom' and 'bust' environmental backgrounds.

Bromodeoxyuridine (BrdU) labeling and subsequent fluorescence activated cell sorting for culture-independent identification of dissolved organic carbon-degrading bacterioplankton

Microbes are major agents mediating the degradation of numerous dissolved organic carbon (DOC) substrates in aquatic environments. However, identification of bacterial taxa that transform specific pools of DOC in nature poses a technical challenge. Here we describe an approach that couples bromodeoxyuridine (BrdU) incorporation, fluorescence activated cell sorting (FACS), and 16S rRNA gene-based molecular analysis that allows culture-independent identification of bacterioplankton capable of degrading a specific DOC compound in aquatic environments. Triplicate bacterioplankton microcosms are set up to receive both BrdU and a model DOC compound (DOC amendments), or only BrdU (no-addition control). BrdU substitutes the positions of thymidine in newly synthesized bacterial DNA and BrdU-labeled DNA can be readily immunodetected. Through a 24-hr incubation, bacterioplankton that are able to use the added DOC compound are expected to be selectively activated, and therefore have higher levels of BrdU incorporation (HI cells) than non-responsive cells in the DOC amendments and cells in no-addition controls (low BrdU incorporation cells, LI cells). After fluorescence immunodetection, HI cells are distinguished and physically separated from the LI cells by fluorescence activated cell sorting (FACS). Sorted DOC-responsive cells (HI cells) are extracted for DNA and taxonomically identified through subsequent 16S rRNA gene-based analyses including PCR, clone library construction and sequencing.

Bacterial Catabolism of Dimethylsulfoniopropionate (DMSP)

Dimethylsulfoniopropionate (DMSP) is a metabolite produced primarily by marine phytoplankton and is the main precursor to the climatically important gas dimethylsulfide (DMS). DMS is released upon bacterial catabolism of DMSP, but it is not the only possible fate of DMSP sulfur. An alternative demethylation/demethiolation pathway results in the eventual release of methanethiol, a highly reactive volatile sulfur compound that contributes little to the atmospheric sulfur flux. The activity of these pathways control the natural flux of sulfur released to the atmosphere. Although these biochemical pathways and the factors that regulate them are of great interest, they are poorly understood. Only recently have some of the genes and pathways responsible for DMSP catabolism been elucidated. Thus far, six different enzymes have been identified that catalyze the cleavage of DMSP, resulting in the release of DMS. In addition, five of these enzymes appear to produce acrylate, while one produces 3-hydroxypropionate. In contrast, only one enzyme, designated DmdA, has been identified that catalyzes the demethylation reaction producing methylmercaptopropionate (MMPA). The metabolism of MMPA is performed by a series of three coenzyme-A mediated reactions catalyzed by DmdB, DmdC, and DmdD. Interestingly, CandidatusPelagibacter ubique, a member of the SAR11 clade of Alphaproteobacteria that is highly abundant in marine surface waters, possessed functional DmdA, DmdB, and DmdC enzymes. Microbially mediated transformations of both DMS and methanethiol are also possible, although many of the biochemical and molecular genetic details are still unknown. This review will focus on the recent discoveries in the biochemical pathways that mineralize and assimilate DMSP carbon and sulfur, as well as the areas for which a comprehensive understanding is still lacking.

Metatranscriptomic analysis of ammonia-oxidizing organisms in an estuarine bacterioplankton assemblage

Quantitative PCR (qPCR) analysis revealed elevated relative abundance (1.8% of prokaryotes) of marine group 1 Crenarchaeota (MG1C) in two samples of southeastern US coastal bacterioplankton, collected in August 2008, compared with samples collected from the same site at different times (mean 0.026%). We analyzed the MG1C sequences in metatranscriptomes from these samples to gain an insight into the metabolism of MG1C population growing in the environment, and for comparison with ammonia-oxidizing bacteria (AOB) in the same samples. Assemblies revealed low diversity within sequences assigned to most individual MG1C open reading frames (ORFs) and high homology with 'Candidatus Nitrosopumilus maritimus' strain SCM1 genome sequences. Reads assigned to ORFs for ammonia uptake and oxidation accounted for 37% of all MG1C transcripts. We did not recover any reads for Nmar_1354-Nmar_1357, proposed to encode components of an alternative, nitroxyl-based ammonia oxidation pathway; however, reads from Nmar_1259 and Nmar_1667, annotated as encoding a multicopper oxidase with homology to nirK, were abundant. Reads assigned to two homologous ORFs (Nmar_1201 and Nmar_1547), annotated as hypothetical proteins were also abundant, suggesting that their unknown function is important to MG1C. Superoxide dismutase and peroxiredoxin-like transcripts were more abundant in the MG1C transcript pool than in the complete metatranscriptome, suggesting an enhanced response to oxidative stress by the MG1C population. qPCR indicated low AOB abundance (0.0010% of prokaryotes), and we found no transcripts related to ammonia oxidation and only one RuBisCO transcript among the transcripts assigned to AOB, suggesting they were not responding to the same environmental cues as the MG1C population.

Quantitative analysis of a deeply sequenced marine microbial metatranscriptome

The potential of metatranscriptomic sequencing to provide insights into the environmental factors that regulate microbial activities depends on how fully the sequence libraries capture community expression (that is, sample-sequencing depth and coverage depth), and the sensitivity with which expression differences between communities can be detected (that is, statistical power for hypothesis testing). In this study, we use an internal standard approach to make absolute (per liter) estimates of transcript numbers, a significant advantage over proportional estimates that can be biased by expression changes in unrelated genes. Coastal waters of the southeastern United States contain 1 × 10(12) bacterioplankton mRNA molecules per liter of seawater (~200 mRNA molecules per bacterial cell). Even for the large bacterioplankton libraries obtained in this study (~500,000 possible protein-encoding sequences in each of two libraries after discarding rRNAs and small RNAs from >1 million 454 FLX pyrosequencing reads), sample-sequencing depth was only 0.00001%. Expression levels of 82 genes diagnostic for transformations in the marine nitrogen, phosphorus and sulfur cycles ranged from below detection (<1 × 10(6) transcripts per liter) for 36 genes (for example, phosphonate metabolism gene phnH, dissimilatory nitrate reductase subunit napA) to >2.7 × 10(9) transcripts per liter (ammonia transporter amt and ammonia monooxygenase subunit amoC). Half of the categories for which expression was detected, however, had too few copy numbers for robust statistical resolution, as would be required for comparative (experimental or time-series) expression studies. By representing whole community gene abundance and expression in absolute units (per volume or mass of environment), 'omics' data can be better leveraged to improve understanding of microbially mediated processes in the ocean.

In situ reverse transcription, an approach to characterize genetic diversity and activities of prokaryotes

Reverse transcription of RNA molecules inside intact bacterial cells was carried out by using reverse transcriptase with a single oligonucleotide complementary to specific 16S rRNA or mRNA sequences. Fluorescently labeled nucleotides were incorporated into each transcribed cDNA inside cells. This protocol is termed in situ reverse transcription (ISRT). In this study, by using species-specific primers targeting unique regions of the 16S rRNA sequences, ISRT was used successfully to detect and enumerate the two lignin-degrading bacteria Microbulbifer hydrolyticus IRE-31 and Sagittula stellata E-37 in culture mixtures and complex enrichment communities selected for lignin degradation. Image analysis revealed that M. hydrolyticus IRE-31 and S. stellata E-37 accounted for approximately 30 and 2%, respectively, of the total bacterial cells in lignin enrichment communities. Populations estimated by ISRT were comparable to those estimated by in situ hybridization (ISH) techniques and to those estimated by hybridization against extracted community DNA. ISRT was also successfully used to detect Pseudomonas putida F1 expressing the todC1 gene in seawater exposed to toluene vapor. ISRT provided a higher signal intensity than ISH, especially when targeting mRNA. The calculated pixel intensities resulting from ISRT were up to 4.2 times greater than those from ISH. This suggests that multiple incorporation of fluorescently labeled nucleotides into cDNA provides a high sensitivity for phylogenetic identification of bacterial populations as well as detection of cells expressing a specific functional gene within complex bacterial communities.

Transporter genes expressed by coastal bacterioplankton in response to dissolved organic carbon

Coastal ocean bacterioplankton control the flow of dissolved organic carbon (DOC) from terrestrial and oceanic sources into the marine food web, and regulate the release of inorganic carbon to atmospheric and offshore reservoirs. While the fate of the chemically complex coastal DOC reservoir has long been recognized as a critical feature of the global carbon budget, it has been problematic to identify both the compounds that serve as major conduits for carbon flux and the roles of individual bacterioplankton taxa in mediating that flux. Here we analyse random libraries of expressed genes from a coastal bacterial community to identify sequences representing DOC-transporting proteins. Predicted substrates of expressed transporter genes indicated that carboxylic acids, compatible solutes, polyamines and lipids may be key components of the biologically labile DOC pool in coastal waters, in addition to canonical bacterial substrates such as amino acids, oligopeptides and carbohydrates. Half of the expressed DOC transporter sequences in this coastal ocean appeared to originate from just eight taxa: Roseobacter, SAR11, Flavobacteriales and five orders of γ-Proteobacteria. While all major taxa expressed transporter genes for some DOC components (e.g. amino acids), there were indications of specialization within the bacterioplankton community for others (e.g. carbohydrates, carboxylic acids and polyamines). Experimental manipulations of the natural DOC pool that increased the concentration of phytoplankton- or vascular plant-derived compounds invoked a readily measured response in bacterial transporter gene expression. This highly resolved view of the potential for carbon flux into heterotrophic bacterioplankton cells identifies possible bioreactive components of the coastal DOC pool and highlights differing ecological roles in carbon turnover for the resident bacterial taxa.

Microbial community response to seawater amendment in low-salinity tidal sediments

Rising sea levels and excessive water withdrawals upstream are making previously freshwater coastal ecosystems saline. Plant and animal responses to variation in the freshwater-saline interface have been well studied in the coastal zone; however, microbial community structure and functional response to seawater intrusion remains relatively unexplored. Here, we used molecular approaches to evaluate the response of the prokaryotic community to controlled changes in porewater salinity levels in freshwater sediments from the Altamaha River, Georgia, USA. This work is a companion to a previously published study describing results from an experiment using laboratory flow-through sediment core bioreactors to document biogeochemical changes as porewater salinity was increased from 0 to 10 over 35 days. As reported in Weston et al. (Biogeochemistry, 77:375-408, 62), porewater chemistry was monitored, and cores were sacrificed at 0, 9, 15, and 35 days, at which time we completed terminal restriction fragment length polymorphism and 16S rRNA clone library analyses of sediment microbial communities. The biogeochemical study documented changes in mineralization pathways in response to artificial seawater additions, with a decline in methanogenesis, a transient increase in iron reduction, and finally a dominance of sulfate reduction. Here, we report that, despite these dramatic and significant changes in microbial activity at the biogeochemical level, no significant differences were found between microbial community composition of control vs. seawater-amended treatments for either Bacterial or Archaeal members. Further, taxa in the seawater-amended treatment community did not become more "marine-like" through time. Our experiment suggests that, as seawater intrudes into freshwater sediments, observed changes in metabolic activity and carbon mineralization on the time scale of weeks are driven more by shifts in gene expression and regulation than by changes in the composition of the microbial community.

The diverse bacterial community in intertidal, anaerobic sediments at Sapelo Island, Georgia

The phylogenetic diversity and composition of the bacterial community in anaerobic sediments from Sapelo Island, GA, USA were examined using 16S rRNA gene libraries. The diversity of this community was comparable to that of soil, and 1,186 clones formed 817 OTUs at 99% sequence similarity. Chao1 estimators for the total richness were also high, at 3,290 OTUs at 99% sequence similarity. The program RDPquery was developed to assign clones to taxonomic groups based upon comparisons to the RDP database. While most clones could be assigned to describe phyla, fewer than 30% of the clones could be assigned to a described order. Similarly, nearly 25% of the clones were only distantly related (<90% sequence similarity) to other environmental clones, illustrating the unique composition of this community. One quarter of the clones were related to one or more undescribed orders within the gamma-Proteobacteria. Other abundant groups included the delta-Proteobacteria, Bacteroidetes, and Cyanobacteria. While these phyla were abundant in other estuarine sediments, the specific members at Sapelo Island appeared to be different from those previously described in other locations, suggesting that great diversity exists between as well as within estuarine intertidal sediments. In spite of the large differences in pore water chemistry with season and depth, differences in the bacterial community were modest over the temporal and spatial scales examined and generally restricted to only certain taxa.

Analyzing gene expression from marine microbial communities using environmental transcriptomics

Analogous to metagenomics, environmental transcriptomics (metatranscriptomics) retrieves and sequences environmental mRNAs from a microbial assemblage without prior knowledge of what genes the community might be expressing. Thus it provides the most unbiased perspective on community gene expression in situ. Environmental transcriptomics protocols are technically difficult since prokaryotic mRNAs generally lack the poly(A) tails that make isolation of eukaryotic messages relatively straightforward ¹ and because of the relatively short half lives of mRNAs ². In addition, mRNAs are much less abundant than rRNAs in total RNA extracts, thus an rRNA background often overwhelms mRNA signals. However, techniques for overcoming some of these difficulties have recently been developed. A procedure for analyzing environmental transcriptomes by creating clone libraries using random primers to reverse-transcribe and amplify environmental mRNAs was recently described was successful in two different natural environments, but results were biased by selection of the random primers used to initiate cDNA synthesis ³. Advances in linear amplification of mRNA obviate the need for random primers in the amplification step and make it possible to use less starting material decreasing the collection and processing time of samples and thereby minimizing RNA degradation ⁴. In vitro transcription methods for amplifying mRNA involve polyadenylating the mRNA and incorporating a T7 promoter onto the 3 end of the transcript. Amplified RNA (aRNA) can then be converted to double stranded cDNA using random hexamers and directly sequenced by pyrosequencing ⁵. A first use of this method at Station ALOHA demonstrated its utility for characterizing microbial community gene expression ⁶.

Comparative day/night metatranscriptomic analysis of microbial communities in the North Pacific subtropical gyre

Metatranscriptomic analyses of microbial assemblages (< 5 microm) from surface water at the Hawaiian Ocean Time-Series (HOT) revealed community-wide metabolic activities and day/night patterns of differential gene expression. Pyrosequencing produced 75 558 putative mRNA reads from a day transcriptome and 75 946 from a night transcriptome. Taxonomic binning of annotated mRNAs indicated that Cyanobacteria contributed a greater percentage of the transcripts (54% of annotated sequences) than expected based on abundance (35% of cell counts and 21% 16S rRNA of libraries), and may represent the most actively transcribing cells in this surface ocean community in both the day and night. Major heterotrophic taxa contributing to the community transcriptome included alpha-Proteobacteria (19% of annotated sequences, most of which were SAR11-related) and gamma-Proteobacteria (4%). The composition of transcript pools was consistent with models of prokaryotic gene expression, including operon-based transcription patterns and an abundance of genes predicted to be highly expressed. Metabolic activities that are shared by many microbial taxa (e.g. glycolysis, citric acid cycle, amino acid biosynthesis and transcription and translation machinery) were well represented among the community transcripts. There was an overabundance of transcripts for photosynthesis, C1 metabolism and oxidative phosphorylation in the day compared with night, and evidence that energy acquisition is coordinated with solar radiation levels for both autotrophic and heterotrophic microbes. In contrast, housekeeping activities such as amino acid biosynthesis, membrane synthesis and repair, and vitamin biosynthesis were overrepresented in the night transcriptome. Direct sequencing of these environmental transcripts has provided detailed information on metabolic and biogeochemical responses of a microbial community to solar forcing.