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

Studying microbial functionality within the gut ecosystem by systems biology

Humans are not autonomous entities. We are all living in a complex environment, interacting not only with our peers, but as true holobionts; we are also very much in interaction with our coexisting microbial ecosystems living on and especially within us, in the intestine. Intestinal microorganisms, often collectively referred to as intestinal microbiota, contribute significantly to our daily energy uptake by breaking down complex carbohydrates into simple sugars, which are fermented to short-chain fatty acids and subsequently absorbed by human cells. They also have an impact on our immune system, by suppressing or enhancing the growth of malevolent and beneficial microbes. Our lifestyle can have a large influence on this ecosystem. What and how much we consume can tip the ecological balance in the intestine. A "western diet" containing mainly processed food will have a different effect on our health than a balanced diet fortified with pre- and probiotics. In recent years, new technologies have emerged, which made a more detailed understanding of microbial communities and ecosystems feasible. This includes progress in the sequencing of PCR-amplified phylogenetic marker genes as well as the collective microbial metagenome and metatranscriptome, allowing us to determine with an increasing level of detail, which microbial species are in the microbiota, understand what these microorganisms do and how they respond to changes in lifestyle and diet. These new technologies also include the use of synthetic and in vitro systems, which allow us to study the impact of substrates and addition of specific microbes to microbial communities at a high level of detail, and enable us to gather quantitative data for modelling purposes. Here, we will review the current state of microbiome research, summarizing the computational methodologies in this area and highlighting possible outcomes for personalized nutrition and medicine.

Proteorhodopsin variability and distribution in the North Pacific Subtropical Gyre

Proteorhodopsin is a light-activated retinal-containing proton pump found in many marine bacteria. These photoproteins are globally distributed in the ocean’s photic zone and are capable of generating a proton motive force across the cell membrane. We investigated the phylogenetic diversity, distribution, and abundance of proteorhodopsin encoding genes in free-living bacterioplankton in the North Pacific Subtropical Gyre, leveraging a gene catalog derived from metagenomic samples from the ocean’s surface to 1000 m depth. Proteorhodopsin genes were identified at all depths sampled, but were most abundant at depths shallower than 200 m. The majority of proteorhodopsin gene sequences (60.9%) belonged to members of the SAR11 lineage, with remaining sequences distributed among other diverse taxa. We observed variations in the conserved residues involved in ion pumping and spectral tuning, and biochemically confirmed four different proton pumping proteorhodopsin motifs, including one unique to deep-water SAR11. We also identified a new group of putative proteorhodopsins having unknown function. Our results reveal a broad organismal and unexpected depth distribution for different proteorhodopsin types, as well as substantial within-taxon variability. These data provide a framework for exploring the ecological relevance of proteorhodopsins and their spatiotemporal variation and function in heterotrophic bacteria in the open ocean.

Candidatus Nitrosocaldus cavascurensis, an Ammonia Oxidizing, Extremely Thermophilic Archaeon with a Highly Mobile Genome

Ammonia oxidizing archaea (AOA) of the phylum Thaumarchaeota are widespread in moderate environments but their occurrence and activity has also been demonstrated in hot springs. Here we present the first enrichment of a thermophilic representative with a sequenced genome, which facilitates the search for adaptive strategies and for traits that shape the evolution of Thaumarchaeota. Candidatus Nitrosocaldus cavascurensis has been enriched from a hot spring in Ischia, Italy. It grows optimally at 68°C under chemolithoautotrophic conditions on ammonia or urea converting ammonia stoichiometrically into nitrite with a generation time of approximately 23 h. Phylogenetic analyses based on ribosomal proteins place the organism as a sister group to all known mesophilic AOA. The 1.58 Mb genome of Ca. N. cavascurensis harbors an amoAXCB gene cluster encoding ammonia monooxygenase and genes for a 3-hydroxypropionate/4-hydroxybutyrate pathway for autotrophic carbon fixation, but also genes that indicate potential alternative energy metabolisms. Although a bona fide gene for nitrite reductase is missing, the organism is sensitive to NO-scavenging, underlining the potential importance of this compound for AOA metabolism. Ca. N. cavascurensis is distinct from all other AOA in its gene repertoire for replication, cell division and repair. Its genome has an impressive array of mobile genetic elements and other recently acquired gene sets, including conjugative systems, a provirus, transposons and cell appendages. Some of these elements indicate recent exchange with the environment, whereas others seem to have been domesticated and might convey crucial metabolic traits.

Stress response of a marine ammonia-oxidizing archaeon informs physiological status of environmental populations

High representation by ammonia-oxidizing archaea (AOA) in marine systems is consistent with their high affinity for ammonia, efficient carbon fixation, and copper (Cu)-centric respiratory system. However, little is known about their response to nutrient stress. We therefore used global transcriptional and proteomic analyses to characterize the response of a model AOA, Nitrosopumilus maritimus SCM1, to ammonia starvation, Cu limitation and Cu excess. Most predicted protein-coding genes were transcribed in exponentially growing cells, and of ~74% detected in the proteome, ~6% were modified by N-terminal acetylation. The general response to ammonia starvation and Cu stress was downregulation of genes for energy generation and biosynthesis. Cells rapidly depleted transcripts for the A and B subunits of ammonia monooxygenase (AMO) in response to ammonia starvation, yet retained relatively high levels of transcripts for the C subunit. Thus, similar to ammonia-oxidizing bacteria, selective retention of amoC transcripts during starvation appears important for subsequent recovery, and also suggests that AMO subunit transcript ratios could be used to assess the physiological status of marine populations. Unexpectedly, cobalamin biosynthesis was upregulated in response to both ammonia starvation and Cu stress, indicating the importance of this cofactor in retaining functional integrity during times of stress.

Metatranscriptomic and metagenomic description of the bacterial nitrogen metabolism in waste water wet oxidation effluents

Anaerobic digestion is a common method for reducing the amount of sludge solids in used waters and enabling biogas production. The wet oxidation process (WOX) improves anaerobic digestion by converting carbon into methane through oxidation of organic compounds. WOX produces effluents rich in ammonia, which must be removed to maintain the activity of methanogens. Ammonia removal from WOX could be biologically operated by aerobic granules. To this end, granulation experiments were conducted in 2 bioreactors containing an activated sludge (AS). For the first time, the dynamics of the microbial community structure and the expression levels of 7 enzymes of the nitrogen metabolism in such active microbial communities were followed in regard to time by metagenomics and metatranscriptomics. It was shown that bacterial communities adapt to the wet oxidation effluent by increasing the expression level of the nitrogen metabolism, suggesting that these biological activities could be a less costly alternative for the elimination of ammonia, resulting in a reduction of the use of chemicals and energy consumption in sewage plants. This study reached a strong sequencing depth (from 4.4 to 7.6 Gb) and enlightened a yet unknown diversity of the microorganisms involved in the nitrogen pathway. Moreover, this approach revealed the abundance and expression levels of specialised enzymes involved in nitrification, denitrification, ammonification, dissimilatory nitrate reduction to ammonium (DNRA) and nitrogen fixation processes in AS.

Integrated structural biology and molecular ecology of N-cycling enzymes from ammonia-oxidizing archaea

Knowledge of the molecular ecology and environmental determinants of ammonia-oxidizing organisms is critical to understanding and predicting the global nitrogen (N) and carbon cycles, but an incomplete biochemical picture hinders in vitro studies of N-cycling enzymes. Although an integrative structural and dynamic characterization at the atomic scale would advance our understanding of function tremendously, structural knowledge of key N-cycling enzymes from ecologically relevant ammonia oxidizers is unfortunately extremely limited. Here, we discuss the challenges and opportunities for examining the ecology of ammonia-oxidizing organisms, particularly uncultivated Thaumarchaeota, through (meta)genome-driven structural biology of the enzymes ammonia monooxygenase (AMO) and nitrite reductase (NirK).

Excess labile carbon promotes the expression of virulence factors in coral reef bacterioplankton

Coastal pollution and algal cover are increasing on many coral reefs, resulting in higher dissolved organic carbon (DOC) concentrations. High DOC concentrations strongly affect microbial activity in reef waters and select for copiotrophic, often potentially virulent microbial populations. High DOC concentrations on coral reefs are also hypothesized to be a determinant for switching microbial lifestyles from commensal to pathogenic, thereby contributing to coral reef degradation, but evidence is missing. In this study, we conducted ex situ incubations to assess gene expression of planktonic microbial populations under elevated concentrations of naturally abundant monosaccharides (glucose, galactose, mannose, and xylose) in algal exudates and sewage inflows. We assembled 27 near-complete (>70%) microbial genomes through metagenomic sequencing and determined associated expression patterns through metatranscriptomic sequencing. Differential gene expression analysis revealed a shift in the central carbohydrate metabolism and the induction of metalloproteases, siderophores, and toxins in Alteromonas, Erythrobacter, Oceanicola, and Alcanivorax populations. Sugar-specific induction of virulence factors suggests a mechanistic link for the switch from a commensal to a pathogenic lifestyle, particularly relevant during increased algal cover and human-derived pollution on coral reefs. Although an explicit test remains to be performed, our data support the hypothesis that increased availability of specific sugars changes net microbial community activity in ways that increase the emergence and abundance of opportunistic pathogens, potentially contributing to coral reef degradation.

Minimal tool set for a prokaryotic circadian clock

BACKGROUND

Circadian clocks are found in organisms of almost all domains including photosynthetic Cyanobacteria, whereby large diversity exists within the protein components involved. In the model cyanobacterium Synechococcus elongatus PCC 7942 circadian rhythms are driven by a unique KaiABC protein clock, which is embedded in a network of input and output factors. Homologous proteins to the KaiABC clock have been observed in Bacteria and Archaea, where evidence for circadian behavior in these domains is accumulating. However, interaction and function of non-cyanobacterial Kai-proteins as well as homologous input and output components remain mainly unclear.

RESULTS

Using a universal BLAST analyses, we identified putative KaiC-based timing systems in organisms outside as well as variations within Cyanobacteria. A systematic analyses of publicly available microarray data elucidated interesting variations in circadian gene expression between different cyanobacterial strains, which might be correlated to the diversity of genome encoded clock components. Based on statistical analyses of co-occurrences of the clock components homologous to Synechococcus elongatus PCC 7942, we propose putative networks of reduced and fully functional clock systems. Further, we studied KaiC sequence conservation to determine functionally important regions of diverged KaiC homologs. Biochemical characterization of exemplary cyanobacterial KaiC proteins as well as homologs from two thermophilic Archaea demonstrated that kinase activity is always present. However, a KaiA-mediated phosphorylation is only detectable in KaiC1 orthologs.

CONCLUSION

Our analysis of 11,264 genomes clearly demonstrates that components of the Synechococcus elongatus PCC 7942 circadian clock are present in Bacteria and Archaea. However, all components are less abundant in other organisms than Cyanobacteria and KaiA, Pex, LdpA, and CdpA are only present in the latter. Thus, only reduced KaiBC-based or even simpler, solely KaiC-based timing systems might exist outside of the cyanobacterial phylum, which might be capable of driving diurnal oscillations.

Imaging mass spectrometry and MS/MS molecular networking reveals chemical interactions among cuticular bacteria and pathogenic fungi associated with fungus-growing ants

The fungus-growing ant-microbe symbiosis is an ideal system to study chemistry-based microbial interactions due to the wealth of microbial interactions described, and the lack of information on the molecules involved therein. In this study, we employed a combination of MALDI imaging mass spectrometry (MALDI-IMS) and MS/MS molecular networking to study chemistry-based microbial interactions in this system. MALDI IMS was used to visualize the distribution of antimicrobials at the inhibition zone between bacteria associated to the ant Acromyrmex echinatior and the fungal pathogen Escovopsis sp. MS/MS molecular networking was used for the dereplication of compounds found at the inhibition zones. We identified the antibiotics actinomycins D, X2 and X_0β, produced by the bacterium Streptomyces CBR38; and the macrolides elaiophylin, efomycin A and efomycin G, produced by the bacterium Streptomyces CBR53.These metabolites were found at the inhibition zones using MALDI IMS and were identified using MS/MS molecular networking. Additionally, three shearinines D, F, and J produced by the fungal pathogen Escovopsis TZ49 were detected. This is the first report of elaiophylins, actinomycin X_0β and shearinines in the fungus-growing ant symbiotic system. These results suggest a secondary prophylactic use of these antibiotics by A. echinatior because of their permanent production by the bacteria.

Marine Group II Dominates Planktonic Archaea in Water Column of the Northeastern South China Sea

Temperature, nutrients, and salinity are among the important factors constraining the distribution and abundance of microorganisms in the ocean. Marine Group II (MGII) belonging to Euryarchaeota commonly dominates the planktonic archaeal community in shallow water and Marine Group I (MGI, now is called Thaumarchaeota) in deeper water in global oceans. Results of quantitative PCR (qPCR) and 454 sequencing in our study, however, showed the dominance of MGII in planktonic archaea throughout the water column of the northeastern South China Sea (SCS) that is characterized by strong water mixing. The abundance of ammonia-oxidizing archaea (AOA) representing the main group of Thaumarchaeota in deeper water in the northeastern SCS was significantly lower than in other oceanic regions. Phylogenetic analysis showed that the top operational taxonomic units (OTUs) of the MGII occurring predominantly below 200 m depth may be unique in the northeastern SCS based on the observation that they are distantly related to known sequences (identity ranging from 90–94%). The abundance of MGII was also significantly correlated with total bacteria in the whole column, which may indicate that MGII and bacteria may have similar physiological or biochemical properties or responses to environmental variation. This study provides valuable information about the dominance of MGII over AOA in both shallow and deep water in the northeastern SCS and highlights the need for comprehensive studies integrating physical, chemical, and microbial oceanography.

Interpreting Microbial Biosynthesis in the Genomic Age: Biological and Practical Considerations

Genome mining has become an increasingly powerful, scalable, and economically accessible tool for the study of natural product biosynthesis and drug discovery. However, there remain important biological and practical problems that can complicate or obscure biosynthetic analysis in genomic and metagenomic sequencing projects. Here, we focus on limitations of available technology as well as computational and experimental strategies to overcome them. We review the unique challenges and approaches in the study of symbiotic and uncultured systems, as well as those associated with biosynthetic gene cluster (BGC) assembly and product prediction. Finally, to explore sequencing parameters that affect the recovery and contiguity of large and repetitive BGCs assembled de novo, we simulate Illumina and PacBio sequencing of the Salinispora tropica genome focusing on assembly of the salinilactam (slm) BGC.

Warming enhances old organic carbon decomposition through altering functional microbial communities

Soil organic matter (SOM) stocks contain nearly three times as much carbon (C) as the atmosphere and changes in soil C stocks may have a major impact on future atmospheric carbon dioxide concentrations and climate. Over the past two decades, much research has been devoted to examining the influence of warming on SOM decomposition in topsoil. Most SOM, however, is old and stored in subsoil. The fate of subsoil SOM under future warming remains highly uncertain. Here, by combining a long-term field warming experiment and a meta-analysis study, we showed that warming significantly increased SOM decomposition in subsoil. We also showed that a decade of warming promoted decomposition of subsoil SOM with turnover times of decades to millennia in a tall grass prairie and this effect was largely associated with shifts in the functional gene structure of microbial communities. By coupling stable isotope probing with metagenomics, we found that microbial communities in warmed soils possessed a higher relative abundance of key functional genes involved in the degradation of organic materials with varying recalcitrance than those in control soils. These findings suggest warming may considerably alter the stability of the vast pool of old SOM in subsoil, contributing to the long-term positive feedback between the C cycle and climate.

Arsenic metabolism in high altitude modern stromatolites revealed by metagenomic analysis

Modern stromatolites thrive only in selected locations in the world. Socompa Lake, located in the Andean plateau at 3570 masl, is one of the numerous extreme Andean microbial ecosystems described over recent years. Extreme environmental conditions include hypersalinity, high UV incidence, and high arsenic content, among others. After Socompa's stromatolite microbial communities were analysed by metagenomic DNA sequencing, taxonomic classification showed dominance of Proteobacteria, Bacteroidetes and Firmicutes, and a remarkably high number of unclassified sequences. A functional analysis indicated that carbon fixation might occur not only by the Calvin-Benson cycle, but also through alternative pathways such as the reverse TCA cycle, and the reductive acetyl-CoA pathway. Deltaproteobacteria were involved both in sulfate reduction and nitrogen fixation. Significant differences were found when comparing the Socompa stromatolite metagenome to the Shark Bay (Australia) smooth mat metagenome: namely, those involving stress related processes, particularly, arsenic resistance. An in-depth analysis revealed a surprisingly diverse metabolism comprising all known types of As resistance and energy generating pathways. While the ars operon was the main mechanism, an important abundance of arsM genes was observed in selected phyla. The data resulting from this work will prove a cornerstone for further studies on this rare microbial community.

Proteomics and comparative genomics of Nitrososphaera viennensis reveal the core genome and adaptations of archaeal ammonia oxidizers

Ammonia-oxidizing archaea (AOA) are among the most abundant microorganisms and key players in the global nitrogen and carbon cycles. They share a common energy metabolism but represent a heterogeneous group with respect to their environmental distribution and adaptions, growth requirements, and genome contents. We report here the genome and proteome of Nitrososphaera viennensis EN76, the type species of the archaeal class Nitrososphaeria of the phylum Thaumarchaeota encompassing all known AOA. N. viennensis is a soil organism with a 2.52-Mb genome and 3,123 predicted protein-coding genes. Proteomic analysis revealed that nearly 50% of the predicted genes were translated under standard laboratory growth conditions. Comparison with genomes of closely related species of the predominantly terrestrial Nitrososphaerales as well as the more streamlined marine Nitrosopumilales [Candidatus (Ca.) order] and the acidophile "Ca. Nitrosotalea devanaterra" revealed a core genome of AOA comprising 860 genes, which allowed for the reconstruction of central metabolic pathways common to all known AOA and expressed in the N. viennensis and "Ca Nitrosopelagicus brevis" proteomes. Concomitantly, we were able to identify candidate proteins for as yet unidentified crucial steps in central metabolisms. In addition to unraveling aspects of core AOA metabolism, we identified specific metabolic innovations associated with the Nitrososphaerales mediating growth and survival in the soil milieu, including the capacity for biofilm formation, cell surface modifications and cell adhesion, and carbohydrate conversions as well as detoxification of aromatic compounds and drugs.

High Throughput Sequencing: An Overview of Sequencing Chemistry

In the present century sequencing is to the DNA science, what gel electrophoresis was to it in the last century. From 1977 to 2016 three generation of the sequencing technologies of various types have been developed. Second and third generation sequencing technologies referred commonly to as next generation sequencing technology, has evolved significantly with increase in sequencing speed, decrease in sequencing cost, since its inception in 2004. GS FLX by 454 Life Sciences/Roche diagnostics, Genome Analyzer, HiSeq, MiSeq and NextSeq by Illumina, Inc., SOLiD by ABI, Ion Torrent by Life Technologies are various type of the sequencing platforms available for second generation sequencing. The platforms available for the third generation sequencing are Helicos™ Genetic Analysis System by SeqLL, LLC, SMRT Sequencing by Pacific Biosciences, Nanopore sequencing by Oxford Nanopore's, Complete Genomics by Beijing Genomics Institute and GnuBIO by BioRad, to name few. The present article is an overview of the principle and the sequencing chemistry of these high throughput sequencing technologies along with brief comparison of various types of sequencing platforms available.

A Sample-to-Sequence Protocol for Genus Targeted Transcriptomic Profiling: Application to Marine Synechococcus

Recent studies using whole community metagenomic and metatranscriptomic approaches are revealing important new insights into the functional potential and activity of natural marine microbial communities. Here, we complement these approaches by describing a complete ocean sample-to-sequence protocol, specifically designed to target a single bacterial genus for purposes of both DNA and RNA profiling using fluorescence activated cell sorting (FACS). The importance of defining and understanding the effects of a sampling protocol are critical if we are to gain meaningful data from environmental surveys. Rigorous pipeline trials with a cultured isolate, Synechococcus sp. BL107 demonstrate that water filtration has a well-defined but limited impact on the transcriptomic profile of this organism, whilst sample storage and multiple rounds of cell sorting have almost no effect on the resulting RNA sequence profile. Attractively, the means to replicate the sampling strategy is within the budget and expertise of most researchers.

Minimal Tool Set for a Prokaryotic Circadian Clock.. [DOI: 10.1101/075291]

Background

Circadian clocks are found in organisms of almost all domains including photosynthetic Cyanobacteria, whereby large diversity exists within the protein components involved. In the model cyanobacteriumSynechococcus elongatusPCC 7942 circadian rhythms are driven by a unique KaiABC protein clock, which is embedded in a network of input and output factors. Homologous proteins to the KaiABC clock have been observed in Bacteria and Archaea, where evidence for circadian behavior in these domains is accumulating. However, interaction and function of non-cyanobacterial Kai-proteins as well as homologous input and output components remain mainly unclear.

Results

Using a universal BLAST analyses, we identified putative KaiC-based timing systems in organisms outside as well as variations within Cyanobacteria. A systematic analyses of publicly available microarray data elucidated interesting variations in circadian gene expression between different cyanobacterial strains, which might be correlated to the diversity of genome encoded clock components. Based on statistical analyses of co-occurrences of the clock components homologous toSynechococcus elongatusPCC 7942, we propose putative networks of reduced and fully functional clock systems. Further, we studied KaiC sequence conservation to determine functionally important regions of diverged KaiC homologs. Biochemical characterization of exemplary cyanobacterial KaiC proteins as well as homologs from two thermophilic Archaea demonstrated that kinase activity is always present. However, a KaiA-mediated phosphorylation is only detectable in KaiC1 orthologs.

Conclusion

Our analysis of 11,264 genomes clearly demonstrates that components of theSynechococcus elongatusPCC 7942 circadian clock are present in Bacteria and Archaea. However, all components are less abundant in other organisms than Cyanobacteria and KaiA, Pex, LdpA, and CdpA are only present in the latter. Thus, only reduced KaiBC-based or even simpler, solely KaiC-based timing systems might exist outside of the cyanobacterial phylum, which might be capable of driving diurnal oscillations.

Marine Bacterial and Archaeal Ion-Pumping Rhodopsins: Genetic Diversity, Physiology, and Ecology

The recognition of a new family of rhodopsins in marine planktonic bacteria, proton-pumping proteorhodopsin, expanded the known phylogenetic range, environmental distribution, and sequence diversity of retinylidene photoproteins. At the time of this discovery, microbial ion-pumping rhodopsins were known solely in haloarchaea inhabiting extreme hypersaline environments. Shortly thereafter, proteorhodopsins and other light-activated energy-generating rhodopsins were recognized to be widespread among marine bacteria. The ubiquity of marine rhodopsin photosystems now challenges prior understanding of the nature and contributions of "heterotrophic" bacteria to biogeochemical carbon cycling and energy fluxes. Subsequent investigations have focused on the biophysics and biochemistry of these novel microbial rhodopsins, their distribution across the tree of life, evolutionary trajectories, and functional expression in nature. Later discoveries included the identification of proteorhodopsin genes in all three domains of life, the spectral tuning of rhodopsin variants to wavelengths prevailing in the sea, variable light-activated ion-pumping specificities among bacterial rhodopsin variants, and the widespread lateral gene transfer of biosynthetic genes for bacterial rhodopsins and their associated photopigments. Heterologous expression experiments with marine rhodopsin genes (and associated retinal chromophore genes) provided early evidence that light energy harvested by rhodopsins could be harnessed to provide biochemical energy. Importantly, some studies with native marine bacteria show that rhodopsin-containing bacteria use light to enhance growth or promote survival during starvation. We infer from the distribution of rhodopsin genes in diverse genomic contexts that different marine bacteria probably use rhodopsins to support light-dependent fitness strategies somewhere between these two extremes.

Species Divergence vs. Functional Convergence Characterizes Crude Oil Microbial Community Assembly

Oil reservoirs exhibit extreme environmental conditions such as high salinity and high temperature. Insights into microbial community assemblages in oil reservoirs and their functional potentials are important for understanding biogeochemical cycles in the subterranean biosphere. In this study, we performed shotgun metagenomic sequencing of crude oil samples from two geographically distant oil reservoirs in China, and compared them with all the 948 available environmental metagenomes deposited in IMG database (until October 2013). Although the dominant bacteria and the proportion of hydrogenotrophic and acetoclastic methanogens were different among oil metagenomes, compared with the metagenomes from other environments, all the oil metagenomes contained higher abundances of genes involved in methanogenic hydrocarbon degradation and stress response systems. In addition, a "shape-sorting" manner was proposed for the assembly of microbial communities in oil reservoirs, with the oil reservoir acting as a function sorter to select microbes with special functions from its endemic pool of microorganisms. At the functional level, we found that environmental metagenomes were generally clustered according to their isolation environments but not their geographical locations, suggesting selective processes to be involved in the assembly of microbial communities based on functional gene composition.

Simultaneous DNA-RNA Extraction from Coastal Sediments and Quantification of 16S rRNA Genes and Transcripts by Real-time PCR

Real Time Polymerase Chain Reaction also known as quantitative PCR (q-PCR) is a widely used tool in microbial ecology to quantify gene abundances of taxonomic and functional groups in environmental samples. Used in combination with a reverse transcriptase reaction (RT-q-PCR), it can also be employed to quantify gene transcripts. q-PCR makes use of highly sensitive fluorescent detection chemistries that allow quantification of PCR amplicons during the exponential phase of the reaction. Therefore, the biases associated with 'end-point' PCR detected in the plateau phase of the PCR reaction are avoided. A protocol to quantify bacterial 16S rRNA genes and transcripts from coastal sediments via real-time PCR is provided. First, a method for the co-extraction of DNA and RNA from coastal sediments, including the additional steps required for the preparation of DNA-free RNA, is outlined. Second, a step-by-step guide for the quantification of 16S rRNA genes and transcripts from the extracted nucleic acids via q-PCR and RT-q-PCR is outlined. This includes details for the construction of DNA and RNA standard curves. Key considerations for the use of RT-q-PCR assays in microbial ecology are included.

Meta-omics in Inflammatory Bowel Disease Research: Applications, Challenges, and Guidelines

Meta-omics [metagenomics, metatranscriptomics, and metaproteomics] are rapidly expanding our knowledge of the gut microbiota in health and disease. These technologies are increasingly used in inflammatory bowel disease [IBD] research. Yet, meta-omics data analysis, interpretation, and among-study comparison remain challenging. In this review we discuss the role these techniques are playing in IBD research, highlighting their strengths and limitations. We give guidelines on proper sample collection and preparation methods, and on performing the analyses and interpreting the results, reporting available user-friendly tools and pipelines.

Ecophysiological Distinctions of Haloarchaea from a Hypersaline Antarctic Lake as Determined by Metaproteomics

Deep Lake in the Vestfold Hills is hypersaline and the coldest system in Antarctica known to support microbial growth (temperatures as low as -20°C). It represents a strong experimental model because the lake supports a low-complexity community of haloarchaea, with the three most abundant species totaling ∼72%. Moreover, the dominant haloarchaea are cultivatable, and their genomes are sequenced. Here we use metaproteomics linked to metagenome data and the genome sequences of the isolates to characterize the main pathways, trophic strategies, and interactions associated with resource utilization. The dominance of the most abundant member, Halohasta litchfieldiae, appears to be predicated on competitive utilization of substrates (e.g., starch, glycerol, and dihydroxyacetone) produced by Dunaliella, the lake's primary producer, while also possessing diverse mechanisms for acquiring nitrogen and phosphorus. The second most abundant member, strain DL31, is proficient in degrading complex proteinaceous matter. Hht. litchfieldiae and DL31 are inferred to release labile substrates that are utilized by Halorubrum lacusprofundi, the third most abundant haloarchaeon in Deep Lake. The study also linked genome variation to specific protein variants or distinct genetic capacities, thereby identifying strain-level variation indicative of specialization. Overall, metaproteomics revealed that rather than functional differences occurring at different lake depths or through size partitioning, the main lake genera possess major trophic distinctions, and phylotypes (e.g., strains of Hht. litchfieldiae) exhibit a more subtle level of specialization. This study highlights the extent to which the lake supports a relatively uniform distribution of taxa that collectively possess the genetic capacity to effectively exploit available nutrients throughout the lake.

IMPORTANCE

Life on Earth has evolved to colonize a broad range of temperatures, but most of the biosphere (∼85%) exists at low temperatures (≤5°C). By performing unique roles in biogeochemical cycles, environmental microorganisms perform functions that are critical for the rest of life on Earth to survive. Cold environments therefore make a particularly important contribution to maintaining healthy, stable ecosystems. Here we describe the main physiological traits of the dominant microorganisms that inhabit Deep Lake in Antarctica, the coldest aquatic environment known to support life. The hypersaline system enables the growth of halophilic members of the Archaea: haloarchaea. By analyzing proteins of samples collected from the water column, we determined the functions that the haloarchaea were likely to perform. This study showed that the dominant haloarchaea possessed distinct lifestyles yet formed a uniform community throughout the lake that was collectively adept at using available light energy and diverse organic substrates for growth.

Metagenomics, Metatranscriptomics, and Metabolomics Approaches for Microbiome Analysis

Microbiomes are ubiquitous and are found in the ocean, the soil, and in/on other living organisms. Changes in the microbiome can impact the health of the environmental niche in which they reside. In order to learn more about these communities, different approaches based on data from multiple omics have been pursued. Metagenomics produces a taxonomical profile of the sample, metatranscriptomics helps us to obtain a functional profile, and metabolomics completes the picture by determining which byproducts are being released into the environment. Although each approach provides valuable information separately, we show that, when combined, they paint a more comprehensive picture. We conclude with a review of network-based approaches as applied to integrative studies, which we believe holds the key to in-depth understanding of microbiomes.

Connectivity between surface and deep waters determines prokaryotic diversity in the North Atlantic Deep Water

To decipher the influence of depth stratification and surface provincialism on the dark ocean prokaryotic community composition, we sampled the major deep‐water masses in the eastern North Atlantic covering three biogeographic provinces. Their diversity was evaluated using ordination and canonical analysis of 454 pyrotag sequences. Variance partitioning suggested that 16% of the variation in the bacterial community composition was based on depth stratification while 9% of the variation was due to geographic location. General linear mixed effect models showed that the community of the subsurface waters was connected to the dark ocean prokaryotic communities in different biogeographic provinces. Cluster analysis indicated that some prokaryotic taxa are specific to distinct regions in bathypelagic water masses. Taken together, our data suggest that the dark ocean prokaryotic community composition of the eastern North Atlantic is primed by the formation and the horizontal transport of water masses.

Distinct Spatial Patterns of SAR11, SAR86, and Actinobacteria Diversity along a Transect in the Ultra-oligotrophic South Pacific Ocean

Distinct distribution patterns of members of the major bacterial clades SAR11, SAR86, and Actinobacteria were observed across a transect from the Marquesas islands through the ultra-oligotrophic South Pacific Gyre into the Chilean upwelling using 16S rRNA gene sequencing and RNA–DNA fingerprinting. Three different Actinobacteria sequence clusters belonging to “Candidatus Actinomarinidae” were localized in the western half of the transect, one was limited to the gyre deep chlorophyll maximum (DCM) and sequences affiliated to the OCS155 clade were unique to the upwelling. The structure of the surface bacterial community was highly correlated with water mass and remained similar across the whole central gyre (1300 nautical miles). The surface hyperoligotrophic gyre was dominated (>70% of all sequences) by highly diverse SAR11 and SAR86 operational taxonomic units and these communities were significantly different from those in the DCM. Analysis of 16S rRNA fingerprints generated from RNA allowed insights into the potential activity of assigned bacterial groups. SAR11 and Prochlorococcus showed the highest potential activity in all water masses except for the upwelling, accounting together for 65% of the total bacterial 16S rRNA in the gyre surface waters in equal proportions whereas the contribution of SAR11 decreased significantly at the DCM.

Deciphering ocean carbon in a changing world

Dissolved organic matter (DOM) in the oceans is one of the largest pools of reduced carbon on Earth, comparable in size to the atmospheric CO2 reservoir. A vast number of compounds are present in DOM, and they play important roles in all major element cycles, contribute to the storage of atmospheric CO2 in the ocean, support marine ecosystems, and facilitate interactions between organisms. At the heart of the DOM cycle lie molecular-level relationships between the individual compounds in DOM and the members of the ocean microbiome that produce and consume them. In the past, these connections have eluded clear definition because of the sheer numerical complexity of both DOM molecules and microorganisms. Emerging tools in analytical chemistry, microbiology, and informatics are breaking down the barriers to a fuller appreciation of these connections. Here we highlight questions being addressed using recent methodological and technological developments in those fields and consider how these advances are transforming our understanding of some of the most important reactions of the marine carbon cycle.

Molecular Ecology of Hypersaline Microbial Mats: Current Insights and New Directions

Microbial mats are unique geobiological ecosystems that form as a result of complex communities of microorganisms interacting with each other and their physical environment. Both the microorganisms present and the network of metabolic interactions govern ecosystem function therein. These systems are often found in a range of extreme environments, and those found in elevated salinity have been particularly well studied. The purpose of this review is to briefly describe the molecular ecology of select model hypersaline mat systems (Guerrero Negro, Shark Bay, S’Avall, and Kiritimati Atoll), and any potentially modulating effects caused by salinity to community structure. In addition, we discuss several emerging issues in the field (linking function to newly discovered phyla and microbial dark matter), which illustrate the changing paradigm that is seen as technology has rapidly advanced in the study of these extreme and evolutionally significant ecosystems.

Experimental Design and Bioinformatics Analysis for the Application of Metagenomics in Environmental Sciences and Biotechnology

Recent advances in DNA sequencing technologies have prompted the widespread application of metagenomics for the investigation of novel bioresources (e.g., industrial enzymes and bioactive molecules) and unknown biohazards (e.g., pathogens and antibiotic resistance genes) in natural and engineered microbial systems across multiple disciplines. This review discusses the rigorous experimental design and sample preparation in the context of applying metagenomics in environmental sciences and biotechnology. Moreover, this review summarizes the principles, methodologies, and state-of-the-art bioinformatics procedures, tools and database resources for metagenomics applications and discusses two popular strategies (analysis of unassembled reads versus assembled contigs/draft genomes) for quantitative or qualitative insights of microbial community structure and functions. Overall, this review aims to facilitate more extensive application of metagenomics in the investigation of uncultured microorganisms, novel enzymes, microbe-environment interactions, and biohazards in biotechnological applications where microbial communities are engineered for bioenergy production, wastewater treatment, and bioremediation.

Functional group-specific traits drive phytoplankton dynamics in the oligotrophic ocean

A diverse microbial assemblage in the ocean is responsible for nearly half of global primary production. It has been hypothesized and experimentally demonstrated that nutrient loading can stimulate blooms of large eukaryotic phytoplankton in oligotrophic systems. Although central to balancing biogeochemical models, knowledge of the metabolic traits that govern the dynamics of these bloom-forming phytoplankton is limited. We used eukaryotic metatranscriptomic techniques to identify the metabolic basis of functional group-specific traits that may drive the shift between net heterotrophy and autotrophy in the oligotrophic ocean. Replicated blooms were simulated by deep seawater (DSW) addition to mimic nutrient loading in the North Pacific Subtropical Gyre, and the transcriptional responses of phytoplankton functional groups were assayed. Responses of the diatom, haptophyte, and dinoflagellate functional groups in simulated blooms were unique, with diatoms and haptophytes significantly (95% confidence) shifting their quantitative metabolic fingerprint from the in situ condition, whereas dinoflagellates showed little response. Significantly differentially abundant genes identified the importance of colimitation by nutrients, metals, and vitamins in eukaryotic phytoplankton metabolism and bloom formation in this system. The variable transcript allocation ratio, used to quantify transcript reallocation following DSW amendment, differed for diatoms and haptophytes, reflecting the long-standing paradigm of phytoplankton r- and K-type growth strategies. Although the underlying metabolic potential of the large eukaryotic phytoplankton was consistently present, the lack of a bloom during the study period suggests a crucial dependence on physical and biogeochemical forcing, which are susceptible to alteration with changing climate.

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

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

Challenges of metagenomics and single-cell genomics approaches for exploring cyanobacterial diversity

Cyanobacteria have played a crucial role in the history of early earth and continue to be instrumental in shaping our planet, yet applications of cutting edge technology have not yet been widely used to explore cyanobacterial diversity. To provide adequate background, we briefly review current sequencing technologies and their innovative uses in genomics and metagenomics. Next, we focus on current cell capture technologies and the challenges of using them with cyanobacteria. We illustrate the utility in coupling breakthroughs in DNA amplification with cell capture platforms, with an example of microfluidic isolation and subsequent targeted amplicon sequencing from individual terrestrial thermophilic cyanobacteria. Single cells of thermophilic, unicellular Synechococcus sp. JA-2-3-B'a(2-13) (Syn OS-B') were sorted in a microfluidic device, lysed, and subjected to whole genome amplification by multiple displacement amplification. We amplified regions from specific CRISPR spacer arrays, which are known to be highly diverse, contain semi-palindromic repeats which form secondary structure, and can be difficult to amplify. Cell capture, lysis, and genome amplification on a microfluidic device have been optimized, setting a stage for further investigations of individual cyanobacterial cells isolated directly from natural populations.

Physiology and evolution of nitrate acquisition in Prochlorococcus

Prochlorococcus is the numerically dominant phototroph in the oligotrophic subtropical ocean and carries out a significant fraction of marine primary productivity. Although field studies have provided evidence for nitrate uptake by Prochlorococcus, little is known about this trait because axenic cultures capable of growth on nitrate have not been available. Additionally, all previously sequenced genomes lacked the genes necessary for nitrate assimilation. Here we introduce three Prochlorococcus strains capable of growth on nitrate and analyze their physiology and genome architecture. We show that the growth of high-light (HL) adapted strains on nitrate is ∼17% slower than their growth on ammonium. By analyzing 41 Prochlorococcus genomes, we find that genes for nitrate assimilation have been gained multiple times during the evolution of this group, and can be found in at least three lineages. In low-light adapted strains, nitrate assimilation genes are located in the same genomic context as in marine Synechococcus. These genes are located elsewhere in HL adapted strains and may often exist as a stable genetic acquisition as suggested by the striking degree of similarity in the order, phylogeny and location of these genes in one HL adapted strain and a consensus assembly of environmental Prochlorococcus metagenome sequences. In another HL adapted strain, nitrate utilization genes may have been independently acquired as indicated by adjacent phage mobility elements; these genes are also duplicated with each copy detected in separate genomic islands. These results provide direct evidence for nitrate utilization by Prochlorococcus and illuminate the complex evolutionary history of this trait.

Differential assemblage of functional units in paddy soil microbiomes

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

The perinatal microbiome and pregnancy: moving beyond the vaginal microbiome

The human microbiome, the collective genome of the microbial community that is on and within us, has recently been mapped. The initial characterization of healthy subjects has provided investigators with a reference population for interrogating the microbiome in metabolic, intestinal, and reproductive health and disease states. Although it is known that bacteria can colonize the vagina, recent metagenomic studies have shown that the vaginal microbiome varies among reproductive age women. Similarly, the richness and diversity of intestinal microbiota also naturally fluctuate among gravidae in both human and nonhuman primates, as well as mice. Moreover, recent evidence suggests that microbiome niches in pregnancy are not limited to maternal body sites, as the placenta appears to harbor a low biomass microbiome that is presumptively established in early pregnancy and varies in association with a remote history of maternal antenatal infection as well as preterm birth. In this article, we will provide a brief overview on metagenomics science as a means to investigate the microbiome, observations pertaining to both variation and the presumptive potential role of a varied microbiome during pregnancy, and how future studies of the microbiome in pregnancy may lend to a better understanding of human biology, reproductive health, and parturition.

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

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

Multi -omics and metabolic modelling pipelines: challenges and tools for systems microbiology

Integrated -omics approaches are quickly spreading across microbiology research labs, leading to (i) the possibility of detecting previously hidden features of microbial cells like multi-scale spatial organization and (ii) tracing molecular components across multiple cellular functional states. This promises to reduce the knowledge gap between genotype and phenotype and poses new challenges for computational microbiologists. We underline how the capability to unravel the complexity of microbial life will strongly depend on the integration of the huge and diverse amount of information that can be derived today from -omics experiments. In this work, we present opportunities and challenges of multi -omics data integration in current systems biology pipelines. We here discuss which layers of biological information are important for biotechnological and clinical purposes, with a special focus on bacterial metabolism and modelling procedures. A general review of the most recent computational tools for performing large-scale datasets integration is also presented, together with a possible framework to guide the design of systems biology experiments by microbiologists.

Comparative metagenomic and metatranscriptomic analyses of microbial communities in acid mine drainage

The microbial communities in acid mine drainage have been extensively studied to reveal their roles in acid generation and adaption to this environment. Lacking, however, are integrated community- and organism-wide comparative gene transcriptional analyses that could reveal the response and adaptation mechanisms of these extraordinary microorganisms to different environmental conditions. In this study, comparative metagenomics and metatranscriptomics were performed on microbial assemblages collected from four geochemically distinct acid mine drainage (AMD) sites. Taxonomic analysis uncovered unexpectedly high microbial biodiversity of these extremely acidophilic communities, and the abundant taxa of Acidithiobacillus, Leptospirillum and Acidiphilium exhibited high transcriptional activities. Community-wide comparative analyses clearly showed that the AMD microorganisms adapted to the different environmental conditions via regulating the expression of genes involved in multiple in situ functional activities, including low-pH adaptation, carbon, nitrogen and phosphate assimilation, energy generation, environmental stress resistance, and other functions. Organism-wide comparative analyses of the active taxa revealed environment-dependent gene transcriptional profiles, especially the distinct strategies used by Acidithiobacillus ferrivorans and Leptospirillum ferrodiazotrophum in nutrients assimilation and energy generation for survival under different conditions. Overall, these findings demonstrate that the gene transcriptional profiles of AMD microorganisms are closely related to the site physiochemical characteristics, providing clues into the microbial response and adaptation mechanisms in the oligotrophic, extremely acidic environments.

Metatranscriptomics of N2-fixing cyanobacteria in the Amazon River plume

Biological N2 fixation is an important nitrogen source for surface ocean microbial communities. However, nearly all information on the diversity and gene expression of organisms responsible for oceanic N2 fixation in the environment has come from targeted approaches that assay only a small number of genes and organisms. Using genomes of diazotrophic cyanobacteria to extract reads from extensive meta-genomic and -transcriptomic libraries, we examined diazotroph diversity and gene expression from the Amazon River plume, an area characterized by salinity and nutrient gradients. Diazotroph genome and transcript sequences were most abundant in the transitional waters compared with lower salinity or oceanic water masses. We were able to distinguish two genetically divergent phylotypes within the Hemiaulus-associated Richelia sequences, which were the most abundant diazotroph sequences in the data set. Photosystem (PS)-II transcripts in Richelia populations were much less abundant than those in Trichodesmium, and transcripts from several Richelia PS-II genes were absent, indicating a prominent role for cyclic electron transport in Richelia. In addition, there were several abundant regulatory transcripts, including one that targets a gene involved in PS-I cyclic electron transport in Richelia. High sequence coverage of the Richelia transcripts, as well as those from Trichodesmium populations, allowed us to identify expressed regions of the genomes that had been overlooked by genome annotations. High-coverage genomic and transcription analysis enabled the characterization of distinct phylotypes within diazotrophic populations, revealed a distinction in a core process between dominant populations and provided evidence for a prominent role for noncoding RNAs in microbial communities.

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.

GeoChip-based insights into the microbial functional gene repertoire of marine sponges (high microbial abundance, low microbial abundance) and seawater

The GeoChip 4.2 gene array was employed to interrogate the microbial functional gene repertoire of sponges and seawater collected from the Red Sea and the Mediterranean. Complementary amplicon sequencing confirmed the microbial community composition characteristic of high microbial abundance (HMA) and low microbial abundance (LMA) sponges. By use of GeoChip, altogether 20,273 probes encoding for 627 functional genes and representing 16 gene categories were identified. Minimum curvilinear embedding analyses revealed a clear separation between the samples. The HMA/LMA dichotomy was stronger than any possible geographic pattern, which is shown here for the first time on the level of functional genes. However, upon inspection of individual genes, very few specific differences were discernible. Differences were related to microbial ammonia oxidation, ammonification, and archaeal autotrophic carbon fixation (higher gene abundance in sponges over seawater) as well as denitrification and radiation-stress-related genes (lower gene abundance in sponges over seawater). Except for few documented specific differences the functional gene repertoire between the different sources appeared largely similar. This study expands previous reports in that functional gene convergence is not only reported between HMA and LMA sponges but also between sponges and seawater.

Dynamics of metabolic activities and gene expression in the Roseobacter clade bacterium Phaeobacter sp. strain MED193 during growth with thiosulfate

Metagenomic analyses of surface seawater reveal that genes for sulfur oxidation are widespread in bacterioplankton communities. However, little is known about the metabolic processes used to exploit the energy potentially gained from inorganic sulfur oxidation in oxic seawater. We therefore studied the sox gene system containing Roseobacter clade isolate Phaeobacter sp. strain MED193 in acetate minimal medium with and without thiosulfate. The addition of thiosulfate enhanced the bacterial growth yields up to 40% in this strain. Concomitantly, soxB and soxY gene expression increased about 8-fold with thiosulfate and remained 11-fold higher than that in controls through stationary phase. At stationary phase, thiosulfate stimulated protein synthesis and anaplerotic CO2 fixation rates up to 5- and 35-fold, respectively. Several genes involved in anaplerotic CO2 fixation (i.e., pyruvate carboxylase, propionyl coenzyme A [CoA], and crotonyl-CoA carboxylase) were highly expressed during active growth, coinciding with high CO2 fixation rates. The high expression of key genes in the ethylmalonyl-CoA pathway suggests that this is an important pathway for the utilization of two-carbon compounds in Phaeobacter sp. MED193. Overall, our findings imply that Roseobacter clade bacteria carrying sox genes can use their lithotrophic potential to gain additional energy from sulfur oxidation for both increasing their growth capacity and improving their long-term survival.

The functional potential of microbial communities in hydraulic fracturing source water and produced water from natural gas extraction characterized by metagenomic sequencing

Microbial activity in produced water from hydraulic fracturing operations can lead to undesired environmental impacts and increase gas production costs. However, the metabolic profile of these microbial communities is not well understood. Here, for the first time, we present results from a shotgun metagenome of microbial communities in both hydraulic fracturing source water and wastewater produced by hydraulic fracturing. Taxonomic analyses showed an increase in anaerobic/facultative anaerobic classes related to Clostridia, Gammaproteobacteria, Bacteroidia and Epsilonproteobacteria in produced water as compared to predominantly aerobic Alphaproteobacteria in the fracturing source water. The metabolic profile revealed a relative increase in genes responsible for carbohydrate metabolism, respiration, sporulation and dormancy, iron acquisition and metabolism, stress response and sulfur metabolism in the produced water samples. These results suggest that microbial communities in produced water have an increased genetic ability to handle stress, which has significant implications for produced water management, such as disinfection.