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

Aerobic nitrous oxide production through N-nitrosating hybrid formation in ammonia-oxidizing archaea

Soil emissions are largely responsible for the increase of the potent greenhouse gas nitrous oxide (N2O) in the atmosphere and are generally attributed to the activity of nitrifying and denitrifying bacteria. However, the contribution of the recently discovered ammonia-oxidizing archaea (AOA) to N2O production from soil is unclear as is the mechanism by which they produce it. Here we investigate the potential of Nitrososphaera viennensis, the first pure culture of AOA from soil, to produce N2O and compare its activity with that of a marine AOA and an ammonia-oxidizing bacterium (AOB) from soil. N. viennensis produced N2O at a maximum yield of 0.09% N2O per molecule of nitrite under oxic growth conditions. N2O production rates of 4.6±0.6 amol N2O cell(-1) h(-1) and nitrification rates of 2.6±0.5 fmol NO2(-) cell(-1) h(-1) were in the same range as those of the AOB Nitrosospira multiformis and the marine AOA Nitrosopumilus maritimus grown under comparable conditions. In contrast to AOB, however, N2O production of the two archaeal strains did not increase when the oxygen concentration was reduced, suggesting that they are not capable of denitrification. In (15)N-labeling experiments we provide evidence that both ammonium and nitrite contribute equally via hybrid N2O formation to the N2O produced by N. viennensis under all conditions tested. Our results suggest that archaea may contribute to N2O production in terrestrial ecosystems, however, they are not capable of nitrifier-denitrification and thus do not produce increasing amounts of the greenhouse gas when oxygen becomes limiting.

Linking activity and function to ecosystem dynamics in a coastal bacterioplankton community

For bacterial communities containing hundreds to thousands of distinct populations, connecting functional processes and environmental dynamics at high taxonomic resolution has remained challenging. Here we use the expression of ribosomal proteins (%RP) as a proxy for in situ activity of 200 taxa within 20 metatranscriptomic samples in a coastal ocean time series encompassing both seasonal variability and diel dynamics. %RP patterns grouped the taxa into seven activity clusters with distinct profiles in functional gene expression and correlations with environmental gradients. Clusters 1-3 had their highest potential activity in the winter and fall, and included some of the most active taxa, while Clusters 4-7 had their highest potential activity in the spring and summer. Cluster 1 taxa were characterized by gene expression for motility and complex carbohydrate degradation (dominated by Gammaproteobacteria and Bacteroidetes), and Cluster 2 taxa by transcription of genes for amino acid and aromatic compound metabolism and aerobic anoxygenic phototrophy (Roseobacter). Other activity clusters were enriched in transcripts for proteorhodopsin and methylotrophy (Cluster 4; SAR11 and methylotrophs), photosynthesis and attachment (Clusters 5 and 7; Synechococcus, picoeukaryotes, Verucomicrobia, and Planctomycetes), and sulfur oxidation (Cluster 7; Gammaproteobacteria). The seasonal patterns in activity were overlain, and sometimes obscured, by large differences in %RP over shorter day-night timescales. Seventy-eight taxa, many of them heterotrophs, had a higher %RP activity index during the day than night, indicating a strong diel activity rhythm at this coastal site. Emerging from these taxonomically- and time-resolved estimates of in situ microbial activity are predictions of specific ecological groupings of microbial taxa in a dynamic coastal environment.

Differential contributions of archaeal ammonia oxidizer ecotypes to nitrification in coastal surface waters

The occurrence of nitrification in the oceanic water column has implications extending from local effects on the structure and activity of phytoplankton communities to broader impacts on the speciation of nitrogenous nutrients and production of nitrous oxide. The ammonia-oxidizing archaea, responsible for carrying out the majority of nitrification in the sea, are present in the marine water column as two taxonomically distinct groups. Water column group A (WCA) organisms are detected at all depths, whereas Water column group B (WCB) are present primarily below the photic zone. An open question in marine biogeochemistry is whether the taxonomic definition of WCA and WCB organisms and their observed distributions correspond to distinct ecological and biogeochemical niches. We used the natural gradients in physicochemical and biological properties that upwelling establishes in surface waters to study their roles in nitrification, and how their activity--ascertained from quantification of ecotype-specific ammonia monooxygenase (amoA) genes and transcripts--varies in response to environmental fluctuations. Our results indicate a role for both ecotypes in nitrification in Monterey Bay surface waters. However, their respective contributions vary, due to their different sensitivities to surface water conditions. WCA organisms exhibited a remarkably consistent level of activity and their contribution to nitrification appears to be related to community size. WCB activity was less consistent and primarily constrained to colder, high nutrient and low chlorophyll waters. Overall, the results of our characterization yielded a strong, potentially predictive, relationship between archaeal amoA gene abundance and the rate of nitrification.

Comparison of assembly algorithms for improving rate of metatranscriptomic functional annotation

BACKGROUND

Microbiome-wide gene expression profiling through high-throughput RNA sequencing ('metatranscriptomics') offers a powerful means to functionally interrogate complex microbial communities. Key to successful exploitation of these datasets is the ability to confidently match relatively short sequence reads to known bacterial transcripts. In the absence of reference genomes, such annotation efforts may be enhanced by assembling reads into longer contiguous sequences ('contigs'), prior to database search strategies. Since reads from homologous transcripts may derive from several species, represented at different abundance levels, it is not clear how well current assembly pipelines perform for metatranscriptomic datasets. Here we evaluate the performance of four currently employed assemblers including de novo transcriptome assemblers - Trinity and Oases; the metagenomic assembler - Metavelvet; and the recently developed metatranscriptomic assembler IDBA-MT.

RESULTS

We evaluated the performance of the assemblers on a previously published dataset of single-end RNA sequence reads derived from the large intestine of an inbred non-obese diabetic mouse model of type 1 diabetes. We found that Trinity performed best as judged by contigs assembled, reads assigned to contigs, and number of reads that could be annotated to a known bacterial transcript. Only 15.5% of RNA sequence reads could be annotated to a known transcript in contrast to 50.3% with Trinity assembly. Paired-end reads generated from the same mouse samples resulted in modest performance gains. A database search estimated that the assemblies are unlikely to erroneously merge multiple unrelated genes sharing a region of similarity (<2% of contigs). A simulated dataset based on ten species confirmed these findings. A more complex simulated dataset based on 72 species found that greater assembly errors were introduced than is expected by sequencing quality. Through the detailed evaluation of assembly performance, the insights provided by this study will help drive the design of future metatranscriptomic analyses.

CONCLUSION

Assembly of metatranscriptome datasets greatly improved read annotation. Of the four assemblers evaluated, Trinity provided the best performance. For more complex datasets, reads generated from transcripts sharing considerable sequence similarity can be a source of significant assembly error, suggesting a need to collate reads on the basis of common taxonomic origin prior to assembly.

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.

Seasonal variation in the metatranscriptomes of a Thaumarchaeota population from SE USA coastal waters

We used a combination of metatranscriptomic analyses and quantitative PCR (qPCR) to study seasonal changes in Thaumarchaeota populations from a salt marsh-dominated estuary. Surface waters (0.5 m depth) were sampled quarterly at Marsh Landing, Sapelo Island, GA, USA over a 3-year period. We found a mid-summer peak in Thaumarchaeota abundance measured by qPCR of either 16S rRNA or amoA genes in each of the 3 years. Thaumarchaeota were 100-1000-fold more abundant during the peak than at other times of the year, whereas the abundance of ammonia- and nitrite-oxidizing Bacteria varied <10-fold over the same period. Analysis of the microdiversity of several highly transcribed genes in 20 metatranscriptomes from a 1-year subset of these samples showed that the transcriptionally active population consisted of 2 or 3 dominant phylotypes that differed between successive summers. This shift appeared to have begun during the preceding winter and spring. Transcripts from the same genes dominated the Thaumarchaeota mRNA pool throughout the year, with genes encoding proteins believed to be involved in nitrogen uptake and oxidation, and two hypothetical proteins being the most abundant transcripts in all libraries. Analysis of individual genes over the seasonal cycle suggested that transcription was tied more closely to variation in growth rates than to seasonal changes in environmental conditions. Day-night differences in the relative abundance of transcripts for ribosomal proteins suggested diurnal variation in Thaumarchaeota growth.

Transcriptional response of the archaeal ammonia oxidizer Nitrosopumilus maritimus to low and environmentally relevant ammonia concentrations

The ability of chemoautotrophic ammonia-oxidizing archaea to compete for ammonia among marine microorganisms at low ambient concentrations has been in part attributed to their extremely high affinity for ammonia, but as yet there is no mechanistic understanding of supporting metabolism. We examined transcription of selected genes for anabolic functions (CO2 fixation, ammonia transport, and cell wall synthesis) and a central catabolic function (ammonia oxidation) in the thaumarchaeon Nitrosopumilus maritimus SCM1 growing at two ammonia concentrations, as measured by combined ammonia and ammonium, one well above the Km for ammonia oxidation (∼500 μM) and the other well below the Km (<10 nM). Transcript levels were generally immediately and differentially repressed when cells transitioned from ammonia-replete to ammonia-limiting conditions. Transcript levels for ammonia oxidation, CO2 fixation, and one of the ammonia transport genes were approximately the same at high and low ammonia availability. Transcripts for all analyzed genes decreased with time in the complete absence of ammonia, but with various rates of decay. The new steady-state mRNA levels established are presumably more reflective of the natural physiological state of ammonia-oxidizing archaea and offer a reference for interpreting message abundance patterns in the natural environment.

An analysis of thaumarchaeota populations from the northern gulf of Mexico

We sampled Thaumarchaeota populations in the northern Gulf of Mexico, including shelf waters under the Mississippi River outflow plume that are subject to recurrent hypoxia. Data from this study allowed us to: (1) test the hypothesis that Thaumarchaeota would be abundant in this region; (2) assess phylogenetic composition of these populations for comparison with other regions; (3) compare the efficacy of quantitative PCR (qPCR) based on primers for 16S rRNA genes (rrs) with primers for genes in the ammonia oxidation (amoA) and carbon fixation (accA, hcd) pathways; (4) compare distributions obtained by qPCR with the relative abundance of Thaumarchaeota rrs in pyrosequenced libraries; (5) compare Thaumarchaeota distributions with environmental variables to help us elucidate the factors responsible for the distributions; (6) compare the distribution of Thaumarchaeota with Nitrite-Oxidizing Bacteria (NOB) to gain insight into the coupling between ammonia and nitrite oxidation. We found up to 10⁸ copies L⁻¹ of Thaumarchaeota rrs in our samples (up to 40% of prokaryotes) by qPCR, with maximum abundance in slope waters at 200–800 m. Thaumarchaeota rrs were also abundant in pyrosequenced libraries and their relative abundance correlated well with values determined by qPCR (r² = 0.82). Thaumarchaeota populations were strongly stratified by depth. Canonical correspondence analysis using a suite of environmental variables explained 92% of the variance in qPCR-estimated gene abundances. Thaumarchaeota rrs abundance was correlated with salinity and depth, while accA abundance correlated with fluorescence and pH. Correlations of Archaeal amoA abundance with environmental variables were primer-dependent, suggesting differential responses of sub-populations to environmental variables. Bacterial amoA was at the limit of qPCR detection in most samples. NOB and Euryarchaeota rrs were found in the pyrosequenced libraries; NOB distribution was correlated with that of Thaumarchaeota (r² = 0.49).

Differential responses of ammonia-oxidizing archaea and bacteria to long-term fertilization in a New England salt marsh

Since the discovery of ammonia-oxidizing archaea (AOA), new questions have arisen about population and community dynamics and potential interactions between AOA and ammonia-oxidizing bacteria (AOB). We investigated the effects of long-term fertilization on AOA and AOB in the Great Sippewissett Marsh, Falmouth, MA, USA to address some of these questions. Sediment samples were collected from low and high marsh habitats in July 2009 from replicate plots that received low (LF), high (HF), and extra high (XF) levels of a mixed NPK fertilizer biweekly during the growing season since 1974. Additional untreated plots were included as controls (C). Terminal restriction fragment length polymorphism analysis of the amoA genes revealed distinct shifts in AOB communities related to fertilization treatment, but the response patterns of AOA were less consistent. Four AOB operational taxonomic units (OTUs) predictably and significantly responded to fertilization, but only one AOA OTU showed a significant pattern. Betaproteobacterial amoA gene sequences within the Nitrosospira-like cluster dominated at C and LF sites, while sequences related to Nitrosomonas spp. dominated at HF and XF sites. We identified some clusters of AOA sequences recovered primarily from high fertilization regimes, but other clusters consisted of sequences recovered from all fertilization treatments, suggesting greater physiological diversity. Surprisingly, fertilization appeared to have little impact on abundance of AOA or AOB. In summary, our data reveal striking patterns for AOA and AOB in response to long-term fertilization, and also suggest a missing link between community composition and abundance and nitrogen processing in the marsh.

Anoxic carbon flux in photosynthetic microbial mats as revealed by metatranscriptomics

Photosynthetic microbial mats possess extraordinary phylogenetic and functional diversity that makes linking specific pathways with individual microbial populations a daunting task. Close metabolic and spatial relationships between Cyanobacteria and Chloroflexi have previously been observed in diverse microbial mats. Here, we report that an expressed metabolic pathway for the anoxic catabolism of photosynthate involving Cyanobacteria and Chloroflexi in microbial mats can be reconstructed through metatranscriptomic sequencing of mats collected at Elkhorn Slough, Monterey Bay, CA, USA. In this reconstruction, Microcoleus spp., the most abundant cyanobacterial group in the mats, ferment photosynthate to organic acids, CO2 and H2 through multiple pathways, and an uncultivated lineage of the Chloroflexi take up these organic acids to store carbon as polyhydroxyalkanoates. The metabolic reconstruction is consistent with metabolite measurements and single cell microbial imaging with fluorescence in situ hybridization and NanoSIMS.

Ecophysiology of an ammonia-oxidizing archaeon adapted to low-salinity habitats

Ammonia oxidation in marine and terrestrial ecosystems plays a pivotal role in the cycling of nitrogen and carbon. Recent discoveries have shown that ammonia-oxidizing archaea (AOA) are both abundant and diverse in these systems, yet very little is known about their physiology. Here we report a physiological analysis of a novel low-salinity-type AOA enriched from the San Francisco Bay estuary, Candidatus Nitrosoarchaeum limnia strain SFB1. N. limnia has a slower growth rate than Nitrosopumilus maritimus and Nitrososphaera viennensis EN76, the only pure AOA isolates described to date, but the growth rate is comparable to the growth of marine AOA enrichment cultures. The growth rate only slightly decreased when N. limnia was grown under lower-oxygen conditions (5.5 % oxygen in the headspace). Although N. limnia was capable of growth at 75 % of seawater salinity, there was a longer lag time, incomplete oxidation of ammonia to nitrite, and slower overall growth rate. Allylthiourea (ATU) only partially inhibited growth and ammonia oxidation by N. limnia at concentrations known to completely inhibit bacterial ammonia oxidation. Using electron microscopy, we confirmed the presence of flagella as suggested by various flagellar biosynthesis genes in the N. limnia genome. We demonstrate that N. limnia is representative of a low-salinity estuarine AOA ecotype and that more than 85 % of its proteins have highest identity to other coastal and estuarine metagenomic sequences. Our findings further highlight the physiology of N. limnia and help explain its ecological adaptation to low-salinity niches.

The genome of the ammonia-oxidizing Candidatus Nitrososphaera gargensis: insights into metabolic versatility and environmental adaptations

The cohort of the ammonia-oxidizing archaea (AOA) of the phylum Thaumarchaeota is a diverse, widespread and functionally important group of microorganisms in many ecosystems. However, our understanding of their biology is still very rudimentary in part because all available genome sequences of this phylum are from members of the Nitrosopumilus cluster. Here we report on the complete genome sequence of Candidatus Nitrososphaera gargensis obtained from an enrichment culture, representing a different evolutionary lineage of AOA frequently found in high numbers in many terrestrial environments. With its 2.83 Mb the genome is much larger than that of other AOA. The presence of a high number of (active) IS elements/transposases, genomic islands, gene duplications and a complete CRISPR/Cas defence system testifies to its dynamic evolution consistent with low degree of synteny with other thaumarchaeal genomes. As expected, the repertoire of conserved enzymes proposed to be required for archaeal ammonia oxidation is encoded by N. gargensis, but it can also use urea and possibly cyanate as alternative ammonia sources. Furthermore, its carbon metabolism is more flexible at the central pyruvate switch point, encompasses the ability to take up small organic compounds and might even include an oxidative pentose phosphate pathway. Furthermore, we show that thaumarchaeota produce cofactor F420 as well as polyhydroxyalkanoates. Lateral gene transfer from bacteria and euryarchaeota has contributed to the metabolic versatility of N. gargensis. This organisms is well adapted to its niche in a heavy metal-containing thermal spring by encoding a multitude of heavy metal resistance genes, chaperones and mannosylglycerate as compatible solute and has the genetic ability to respond to environmental changes by signal transduction via a large number of two-component systems, by chemotaxis and flagella-mediated motility and possibly even by gas vacuole formation. These findings extend our understanding of thaumarchaeal evolution and physiology and offer many testable hypotheses for future experimental research on these nitrifiers.

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.

Diversity, physiology, and niche differentiation of ammonia-oxidizing archaea

Nitrification, the aerobic oxidation of ammonia to nitrate via nitrite, has been suggested to have been a central part of the global biogeochemical nitrogen cycle since the oxygenation of Earth. The cultivation of several ammonia-oxidizing archaea (AOA) as well as the discovery that archaeal ammonia monooxygenase (amo)-like gene sequences are nearly ubiquitously distributed in the environment and outnumber their bacterial counterparts in many habitats fundamentally revised our understanding of nitrification. Surprising insights into the physiological distinctiveness of AOA are mirrored by the recognition of the phylogenetic uniqueness of these microbes, which fall within a novel archaeal phylum now known as Thaumarchaeota. The relative importance of AOA in nitrification, compared to ammonia-oxidizing bacteria (AOB), is still under debate. This minireview provides a synopsis of our current knowledge of the diversity and physiology of AOA, the factors controlling their ecology, and their role in carbon cycling as well as their potential involvement in the production of the greenhouse gas nitrous oxide. It emphasizes the importance of activity-based analyses in AOA studies and formulates priorities for future research.

Diversity and evolution of bioenergetic systems involved in microbial nitrogen compound transformations

Nitrogen is an essential element of life that needs to be assimilated in its most reduced form, ammonium. On the other hand, nitrogen exists in a multitude of oxidation states and, consequently, nitrogen compounds (NCs) serve as electron donor and/or acceptors in many catabolic pathways including various forms of microbial respiration that contribute to the global biogeochemical nitrogen cycle. Some of these NCs are also known as reactive nitrogen species able to cause nitrosative stress because of their high redox reactivity. The best understood processes of the nitrogen cycle are denitrification and ammonification (both beginning with nitrate reduction to nitrite), nitrification (aerobic oxidation of ammonium and nitrite) and anaerobic ammonium oxidation (anammox). This review presents examples of the diverse architecture, either elucidated or anticipated, and the high degree of modularity of the corresponding respiratory electron transport processes found in Bacteria and Archaea, and relates these to their respective bioenergetic mechanisms of proton motive force generation. In contrast to the multiplicity of enzymes that catalyze NC transformations, the number of proteins or protein modules involved in connecting electron transport to and from these enzymes with the quinone/quinol pool is comparatively small. These quinone/quinol-reactive protein modules consist of cytochromes b and c and iron-sulfur proteins. Conclusions are drawn towards the evolutionary relationships of bioenergetic systems involved in NC transformation and deduced aspects of the evolution of the biogeochemical nitrogen cycle are presented. This article is part of a Special Issue entitled: The evolutionary aspects of bioenergetic systems.

Genome-enabled transcriptomics reveals archaeal populations that drive nitrification in a deep-sea hydrothermal plume

Ammonia-oxidizing Archaea (AOA) are among the most abundant microorganisms in the oceans and have crucial roles in biogeochemical cycling of nitrogen and carbon. To better understand AOA inhabiting the deep sea, we obtained community genomic and transcriptomic data from ammonium-rich hydrothermal plumes in the Guaymas Basin (GB) and from surrounding deep waters of the Gulf of California. Among the most abundant and active lineages in the sequence data were marine group I (MGI) Archaea related to the cultured autotrophic ammonia-oxidizer, Nitrosopumilus maritimus. Assembly of MGI genomic fragments yielded 2.9 Mb of sequence containing seven 16S rRNA genes (95.4-98.4% similar to N. maritimus), including two near-complete genomes and several lower-abundance variants. Equal copy numbers of MGI 16S rRNA genes and ammonia monooxygenase genes and transcription of ammonia oxidation genes indicates that all of these genotypes actively oxidize ammonia. De novo genomic assembly revealed the functional potential of MGI populations and enhanced interpretation of metatranscriptomic data. Physiological distinction from N. maritimus is evident in the transcription of novel genes, including genes for urea utilization, suggesting an alternative source of ammonia. We were also able to determine which genotypes are most active in the plume. Transcripts involved in nitrification were more prominent in the plume and were among the most abundant transcripts in the community. These unique data sets reveal populations of deep-sea AOA thriving in the ammonium-rich GB that are related to surface types, but with key genomic and physiological differences.

Diversity, abundance and expression of nitrite reductase (nirK)-like genes in marine thaumarchaea

Ammonia-oxidizing archaea (AOA) are widespread and abundant in aquatic and terrestrial habitats and appear to have a significant impact on the global nitrogen cycle. Like the ammonia-oxidizing bacteria, AOA encode a gene homologous to copper-containing nitrite reductases (nirK), which has been studied very little to date. In this study, the diversity, abundance and expression of thaumarchaeal nirK genes from coastal and marine environments were investigated using two mutually excluding primer pairs, which amplify the nirK variants designated as AnirKa and AnirKb. Only the AnirKa variant could be detected in sediment samples from San Francisco Bay and these sequences grouped with the nirK from Candidatus Nitrosopumilus maritimus and Candidatus Nitrosoarchaeum limnia. The two nirK variants had contrasting distributions in the water column in Monterey Bay and the California Current. AnirKa was more abundant in the epi- to mesopelagic Monterey Bay water column, whereas AnirKb was more abundant in the meso- to bathypelagic California Current water. The abundance and community composition of AnirKb, but not AnirKa, followed that of thaumarchaeal amoA, suggesting that either AnirKa is not exclusively associated with AOA or that commonly used amoA primers may be missing a significant fraction of AOA diversity in the epipelagic. Interestingly, thaumarchaeal nirK was expressed 10-100-fold more than amoA in Monterey Bay. Overall, this study provides valuable new insights into the distribution, diversity, abundance and expression of this alternative molecular marker for AOA in the ocean.

Metatranscriptomics of the marine sponge Geodia barretti: tackling phylogeny and function of its microbial community

Geodia barretti is a marine cold-water sponge harbouring high numbers of microorganisms. Significant rates of nitrification have been observed in this sponge, indicating a substantial contribution to nitrogen turnover in marine environments with high sponge cover. In order to get closer insights into the phylogeny and function of the active microbial community and the interaction with its host G. barretti, a metatranscriptomic approach was employed, using the simultaneous analysis of rRNA and mRNA. Of the 262 298 RNA-tags obtained by pyrosequencing, 92% were assigned to ribosomal RNA (ribo-tags). A total of 109 325 SSU rRNA ribo-tags revealed a detailed picture of the community, dominated by group SAR202 of Chloroflexi, candidate phylum Poribacteria and Acidobacteria, which was different in its composition from that obtained in clone libraries prepared form the same samples. Optimized assembly strategies allowed the reconstruction of full-length rRNA sequences from the short ribo-tags for more detailed phylogenetic studies of the dominant taxa. Cells of several phyla were visualized by FISH analyses for confirmation. Of the remaining 21 325 RNA-tags, 10 023 were assigned to mRNA-tags, based on similarities to genes in the databases. A wide range of putative functional gene transcripts from over 10 different phyla were identified among the bacterial mRNA-tags. The most abundant mRNAs were those encoding key metabolic enzymes of nitrification from ammonia-oxidizing archaea as well as candidate genes involved in related processes. Our analysis demonstrates the potential and limits of using a combined rRNA and mRNA approach to explore the microbial community profile, phylogenetic assignments and metabolic activities of a complex, but little explored microbial community.

Effective gene collection from the metatranscriptome of marine microorganisms

Background

Metagenomic studies, accelerated by the evolution of sequencing technologies and the rapid development of genomic analysis methods, can reveal genetic diversity and biodiversity in various samples including those of uncultured or unknown species. This approach, however, cannot be used to identify active functional genes under actual environmental conditions. Metatranscriptomics, which is similar in approach to metagenomics except that it utilizes RNA samples, is a powerful tool for the transcriptomic study of environmental samples. Unlike metagenomic studies, metatranscriptomic studies have not been popular to date due to problems with reliability, repeatability, redundancy and cost performance. Here, we propose a normalized metatranscriptomic method that is suitable for the collection of genes from samples as a platform for comparative transcriptomics.

Results

We constructed two libraries, one non-normalized and the other normalized library, from samples of marine microorganisms taken during daylight hours from Hiroshima bay in Japan. We sequenced 0.6M reads for each sample on a Roche GS FLX, and obtained 0.2M genes after quality control and assembly. A comparison of the two libraries showed that the number of unique genes was larger in the normalized library than in the non-normalized library. Functional analysis of genes revealed that a small number of gene groups, ribosomal RNA genes and chloroplast genes, were dominant in both libraries. Taxonomic distribution analysis of the libraries suggests that Stramenopiles form a major taxon that includes diatoms. The normalization technique thus increases unique genes, functional categories of genes, and taxonomic richness.

Conclusions

Normalization of the marine metatranscriptome could be useful in increasing the number of genes collected, and in reducing redundancies among highly expressed genes. Gene collection through the normalization method was effective in providing a foundation for comparative transcriptomic analysis.

Environmental factors determining ammonia-oxidizing organism distribution and diversity in marine environments

Ammonia-oxidizing bacteria (AOB) and archaea (AOA) play a vital role in bridging the input of fixed nitrogen, through N-fixation and remineralization, to its loss by denitrification and anammox. Yet the major environmental factors determining AOB and AOA population dynamics are little understood, despite both groups having a wide environmental distribution. This study examined the relative abundance of both groups of ammonia-oxidizing organisms (AOO) and the diversity of AOA across large-scale gradients in temperature, salinity and substrate concentration and dissolved oxygen. The relative abundance of AOB and AOA varied across environments, with AOB dominating in the freshwater region of the Chesapeake Bay and AOA more abundant in the water column of the coastal and open ocean. The highest abundance of the AOA amoA gene was recorded in the oxygen minimum zones (OMZs) of the Eastern Tropical South Pacific (ETSP) and the Arabian Sea (AS). The ratio of AOA : AOB varied from 0.7 in the Chesapeake Bay to 1600 in the Sargasso Sea. Relative abundance of both groups strongly correlated with ammonium concentrations. AOA diversity, as determined by phylogenetic analysis of clone library sequences and archetype analysis from a functional gene DNA microarray, detected broad phylogenetic differences across the study sites. However, phylogenetic diversity within physicochemically congruent stations was more similar than would be expected by chance. This suggests that the prevailing geochemistry, rather than localized dispersal, is the major driving factor determining OTU distribution.

Measuring unbiased metatranscriptomics in suboxic waters of the central Baltic Sea using a new in situ fixation system

An analysis of the microbial metabolism is fundamental to understanding globally important element transformations. One culture-independent approach to deduce those prokaryotic metabolic functions is to analyze metatranscriptomes. Unfortunately, since mRNA is extremely labile, it is unclear whether the abundance patterns detected in nature are vulnerable to considerable modification in situ simply due to sampling procedures. Exemplified on comparisons of metatranscriptomes retrieved from pelagic suboxic zones of the central Baltic Sea (70-120 m depth), earlier identified as areas of high aerobic ammonium oxidation activity, and quantification of specific transcripts in them, we show that different sampling techniques significantly influence the relative abundance of transcripts presumably diagnostic of the habitat. In situ fixation using our newly developed automatic flow injection sampler resulted in an abundance of thaumarchaeal ammonia monooxygenase transcripts that was up to 30-fold higher than that detected in samples obtained using standard oceanographic sampling systems. By contrast, the abundance of transcripts indicative of cellular stress was significantly greater in non-fixed samples. Thus, the importance of in situ fixation in the reliable evaluation of distinct microbial activities in the ecosystem based on metatranscriptomics is obvious. In consequence, our data indicate that the significance of thaumarchaeota to aerobic ammonium oxidation could yet have been considerably underestimated. Taken these results, this could in general also be the case in attempts aimed at an unbiased gene expression analysis of areas below the epipelagic zone, which cover 90% of the world's oceans.

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.