A compendium of multi-omic sequence information from the Saanich Inlet water column

Microbial community respiration kinetics and their dynamics in coastal seawater

Oxygen (O2) concentrations in coastal seawater have been declining for decades and models predict continued deoxygenation into the future. As O2 declines, metabolic energy use is progressively channelled from higher trophic levels into microbial community respiration, which in turn influences coastal ecology and biogeochemistry. Despite its critical role in deoxygenation and ecosystem functioning, the kinetics of microbial respiration at low O2 concentrations in coastal seawater remain uncertain and are mostly modeled based on parameters derived from laboratory cultures and a limited number of environmental observations. To explore microbial responses to declining O2, we measured respiration kinetics in coastal microbial communities in Hong Kong over the course of an entire year. We found the mean maximum respiration rate (Vmax) ranged between 560 ± 280 and 5930 ± 800 nmol O2 L-1 h-1, with apparent half-saturation constants (Km) for O2 uptake of between 50 ± 40 and 310 ± 260 nmol O2 L-1. These kinetic parameters vary seasonally in association with shifts in microbial community composition that were linked to nutrient availability, temperature, and biological productivity. In general, coastal communities in Hong Kong exhibited low affinities for O2, yet communities in the dry season had higher affinities, which may play a key role in shaping the relationship between community size, biomass, and O2 consumption rates through respiration. Overall, parameters derived from these experiments can be employed in models to predict the expansion of deoxygenated waters and associated effects on coastal ecology and biogeochemistry.

Development of Nested Polymerase Chain Reaction with Novel Specific Primers for Detection of Tritrichomonas muris Infection in Laboratory Mice

A variety of rodent ceca are parasitized by Tritrichomonas muris (T. muris), a flagellated protozoan. To date, there are no ideal methods for the detection of T. muris infections in laboratory mice; thus, new molecular methodologies for its specific detection need to be developed. In this study, using staining and SEM, it was observed that T. muris has a pear-shaped body and contains three anterior flagella. A nested PCR system with novel specific primers was designed based on the conserved regions of the SSU rRNA gene of T. muris. The nested PCR system for T. muris showed good specificity and high sensitivity for at least 100 T. muris trophozoites/mL and 0.1 ng/μL of fecal genomic DNA, which means that 176 trophozoites per gram of mouse feces could be detected. When using this nested PCR system, the detection rate was 18.96% (58/306), which was higher than the detection rate of 14.05% (43/306) detected via smear microscopy in fecal samples from five mouse strains. The sensitivity and specificity of nested PCR in detecting T. muris was found to be 100%, and it demonstrated a 26% increase in diagnostic sensitivity compared to the smear microscopy method in the present study. In conclusion, the nested PCR developed with novel primers based on the SSU rRNA gene of T. muris has good accuracy, specificity, and sensitivity for the detection of T. muris infections in laboratory mice.

Community- and genome-based evidence for a shaping influence of redox potential on bacterial protein evolution

Despite deep interest in how environments shape microbial communities, whether redox conditions influence the sequence composition of genomes is not well known. We predicted that the carbon oxidation state (ZC) of protein sequences would be positively correlated with redox potential (Eh). To test this prediction, we used taxonomic classifications for 68 publicly available 16S rRNA gene sequence data sets to estimate the abundances of archaeal and bacterial genomes in river & seawater, lake & pond, geothermal, hyperalkaline, groundwater, sediment, and soil environments. Locally, ZC of community reference proteomes (i.e., all the protein sequences in each genome, weighted by taxonomic abundances but not by protein abundances) is positively correlated with Eh corrected to pH 7 (Eh7) for the majority of data sets for bacterial communities in each type of environment, and global-scale correlations are positive for bacterial communities in all environments. In contrast, archaeal communities show approximately equal frequencies of positive and negative correlations in individual data sets, and a positive pan-environmental correlation for archaea only emerges after limiting the analysis to samples with reported oxygen concentrations. These results provide empirical evidence that geochemistry modulates genome evolution and may have distinct effects on bacteria and archaea. IMPORTANCE The identification of environmental factors that influence the elemental composition of proteins has implications for understanding microbial evolution and biogeography. Millions of years of genome evolution may provide a route for protein sequences to attain incomplete equilibrium with their chemical environment. We developed new tests of this chemical adaptation hypothesis by analyzing trends of the carbon oxidation state of community reference proteomes for microbial communities in local- and global-scale redox gradients. The results provide evidence for widespread environmental shaping of the elemental composition of protein sequences at the community level and establish a rationale for using thermodynamic models as a window into geochemical effects on microbial community assembly and evolution.

A compendium of bacterial and archaeal single-cell amplified genomes from oxygen deficient marine waters

Oxygen-deficient marine waters referred to as oxygen minimum zones (OMZs) or anoxic marine zones (AMZs) are common oceanographic features. They host both cosmopolitan and endemic microorganisms adapted to low oxygen conditions. Microbial metabolic interactions within OMZs and AMZs drive coupled biogeochemical cycles resulting in nitrogen loss and climate active trace gas production and consumption. Global warming is causing oxygen-deficient waters to expand and intensify. Therefore, studies focused on microbial communities inhabiting oxygen-deficient regions are necessary to both monitor and model the impacts of climate change on marine ecosystem functions and services. Here we present a compendium of 5,129 single-cell amplified genomes (SAGs) from marine environments encompassing representative OMZ and AMZ geochemical profiles. Of these, 3,570 SAGs have been sequenced to different levels of completion, providing a strain-resolved perspective on the genomic content and potential metabolic interactions within OMZ and AMZ microbiomes. Hierarchical clustering confirmed that samples from similar oxygen concentrations and geographic regions also had analogous taxonomic compositions, providing a coherent framework for comparative community analysis.

The OceanDNA MAG catalog contains over 50,000 prokaryotic genomes originated from various marine environments

Marine microorganisms are immensely diverse and play fundamental roles in global geochemical cycling. Recent metagenome-assembled genome studies, with particular attention to large-scale projects such as Tara Oceans, have expanded the genomic repertoire of marine microorganisms. However, published marine metagenome data is still underexplored. We collected 2,057 marine metagenomes covering various marine environments and developed a new genome reconstruction pipeline. We reconstructed 52,325 qualified genomes composed of 8,466 prokaryotic species-level clusters spanning 59 phyla, including genomes from the deep-sea characterized as deeper than 1,000 m (n = 3,337), low-oxygen zones of <90 μmol O2 per kg water (n = 7,884), and polar regions (n = 7,752). Novelty evaluation using a genome taxonomy database shows that 6,256 species (73.9%) are novel and include genomes of high taxonomic novelty, such as new class candidates. These genomes collectively expanded the known phylogenetic diversity of marine prokaryotes by 34.2%, and the species representatives cover 26.5-42.0% of prokaryote-enriched metagenomes. Thoroughly leveraging accumulated metagenomic data, this genome resource, named the OceanDNA MAG catalog, illuminates uncharacterized marine microbial 'dark matter' lineages.

Variation in genomic traits of microbial communities among ecosystems

Free-living bacteria in nutrient limited environments often exhibit traits which may reduce the cost of reproduction, such as smaller genome size, low GC content and fewer sigma (σ) factor and 16S rRNA gene copies. Despite the potential utility of these traits to detect relationships between microbial communities and ecosystem-scale properties, few studies have assessed these traits on a community-scale. Here, we analysed these traits from publicly available metagenomes derived from marine, soil, host-associated and thermophilic communities. In marine and thermophilic communities, genome size and GC content declined in parallel, consistent with genomic streamlining, with GC content in thermophilic communities generally higher than in marine systems. In contrast, soil communities averaging smaller genomes featured higher GC content and were often from low-carbon environments, suggesting unique selection pressures in soil bacteria. The abundance of specific σ-factors varied with average genome size and ecosystem type. In oceans, abundance of fliA, a σ-factor controlling flagella biosynthesis, was positively correlated with community average genome size—reflecting known trade-offs between nutrient conservation and chemotaxis. In soils, a high abundance of the stress response σ-factor gene rpoS was associated with smaller average genome size and often located in harsh and/or carbon-limited environments—a result which tracks features observed in culture and indicates an increased capacity for stress response in nutrient-poor soils. This work shows how ecosystem-specific constraints are associated with trade-offs which are embedded in the genomic features of bacteria in microbial communities, and which can be detected at the community level, highlighting the importance of genomic features in microbial community analysis.

Pangenomics reveals alternative environmental lifestyles among chlamydiae

Chlamydiae are highly successful strictly intracellular bacteria associated with diverse eukaryotic hosts. Here we analyzed metagenome-assembled genomes of the "Genomes from Earth's Microbiomes" initiative from diverse environmental samples, which almost double the known phylogenetic diversity of the phylum and facilitate a highly resolved view at the chlamydial pangenome. Chlamydiae are defined by a relatively large core genome indicative of an intracellular lifestyle, and a highly dynamic accessory genome of environmental lineages. We observe chlamydial lineages that encode enzymes of the reductive tricarboxylic acid cycle and for light-driven ATP synthesis. We show a widespread potential for anaerobic energy generation through pyruvate fermentation or the arginine deiminase pathway, and we add lineages capable of molecular hydrogen production. Genome-informed analysis of environmental distribution revealed lineage-specific niches and a high abundance of chlamydiae in some habitats. Together, our data provide an extended perspective of the variability of chlamydial biology and the ecology of this phylum of intracellular microbes.

Mercury methylation by metabolically versatile and cosmopolitan marine bacteria

Microbes transform aqueous mercury (Hg) into methylmercury (MeHg), a potent neurotoxin that accumulates in terrestrial and marine food webs, with potential impacts on human health. This process requires the gene pair hgcAB, which encodes for proteins that actuate Hg methylation, and has been well described for anoxic environments. However, recent studies report potential MeHg formation in suboxic seawater, although the microorganisms involved remain poorly understood. In this study, we conducted large-scale multi-omic analyses to search for putative microbial Hg methylators along defined redox gradients in Saanich Inlet, British Columbia, a model natural ecosystem with previously measured Hg and MeHg concentration profiles. Analysis of gene expression profiles along the redoxcline identified several putative Hg methylating microbial groups, including Calditrichaeota, SAR324 and Marinimicrobia, with the last the most active based on hgc transcription levels. Marinimicrobia hgc genes were identified from multiple publicly available marine metagenomes, consistent with a potential key role in marine Hg methylation. Computational homology modelling predicts that Marinimicrobia HgcAB proteins contain the highly conserved amino acid sites and folding structures required for functional Hg methylation. Furthermore, a number of terminal oxidases from aerobic respiratory chains were associated with several putative novel Hg methylators. Our findings thus reveal potential novel marine Hg-methylating microorganisms with a greater oxygen tolerance and broader habitat range than previously recognized.

Mercury methylation by metabolically versatile and cosmopolitan marine bacteria

Microbes transform aqueous mercury (Hg) into methylmercury (MeHg), a potent neurotoxin that accumulates in terrestrial and marine food webs, with potential impacts on human health. This process requires the gene pair hgcAB, which encodes for proteins that actuate Hg methylation, and has been well described for anoxic environments. However, recent studies report potential MeHg formation in suboxic seawater, although the microorganisms involved remain poorly understood. In this study, we conducted large-scale multi-omic analyses to search for putative microbial Hg methylators along defined redox gradients in Saanich Inlet, British Columbia, a model natural ecosystem with previously measured Hg and MeHg concentration profiles. Analysis of gene expression profiles along the redoxcline identified several putative Hg methylating microbial groups, including Calditrichaeota, SAR324 and Marinimicrobia, with the last the most active based on hgc transcription levels. Marinimicrobia hgc genes were identified from multiple publicly available marine metagenomes, consistent with a potential key role in marine Hg methylation. Computational homology modelling predicts that Marinimicrobia HgcAB proteins contain the highly conserved amino acid sites and folding structures required for functional Hg methylation. Furthermore, a number of terminal oxidases from aerobic respiratory chains were associated with several putative novel Hg methylators. Our findings thus reveal potential novel marine Hg-methylating microorganisms with a greater oxygen tolerance and broader habitat range than previously recognized.

Evaluating and Improving Small Subunit rRNA PCR Primer Coverage for Bacteria, Archaea, and Eukaryotes Using Metagenomes from Global Ocean Surveys

Small subunit rRNA (SSU rRNA) amplicon sequencing can quantitatively and comprehensively profile natural microbiomes, representing a critically important tool for studying diverse global ecosystems. However, results will only be accurate if PCR primers perfectly match the rRNA of all organisms present. To evaluate how well marine microorganisms across all 3 domains are detected by this method, we compared commonly used primers with >300 million rRNA gene sequences retrieved from globally distributed marine metagenomes. The best-performing primers compared to 16S rRNA of bacteria and archaea were 515Y/926R and 515Y/806RB, which perfectly matched over 96% of all sequences. Considering cyanobacterial and chloroplast 16S rRNA, 515Y/926R had the highest coverage (99%), making this set ideal for quantifying marine primary producers. For eukaryotic 18S rRNA sequences, 515Y/926R also performed best (88%), followed by V4R/V4RB (18S rRNA specific; 82%)-demonstrating that the 515Y/926R combination performs best overall for all 3 domains. Using Atlantic and Pacific Ocean samples, we demonstrate high correspondence between 515Y/926R amplicon abundances (generated for this study) and metagenomic 16S rRNA (median R² = 0.98, n = 272), indicating amplicons can produce equally accurate community composition data compared with shotgun metagenomics. Our analysis also revealed that expected performance of all primer sets could be improved with minor modifications, pointing toward a nearly completely universal primer set that could accurately quantify biogeochemically important taxa in ecosystems ranging from the deep sea to the surface. In addition, our reproducible bioinformatic workflow can guide microbiome researchers studying different ecosystems or human health to similarly improve existing primers and generate more accurate quantitative amplicon data. IMPORTANCE PCR amplification and sequencing of marker genes is a low-cost technique for monitoring prokaryotic and eukaryotic microbial communities across space and time but will work optimally only if environmental organisms match PCR primer sequences exactly. In this study, we evaluated how well primers match globally distributed short-read oceanic metagenomes. Our results demonstrate that primer sets vary widely in performance, and that at least for marine systems, rRNA amplicon data from some primers lack significant biases compared to metagenomes. We also show that it is theoretically possible to create a nearly universal primer set for diverse saline environments by defining a specific mixture of a few dozen oligonucleotides, and present a software pipeline that can guide rational design of primers for any environment with available meta'omic data.

High-throughput molecular analyses of microbiomes as a tool to monitor the wellbeing of aquatic environments

Aquatic environments are the recipients of many sources of environmental stress that trigger both local and global changes. To evaluate the associated risks to organisms and ecosystems more sensitive and accurate strategies are required. The analysis of the microbiome is one of the most promising candidates for environmental diagnosis of aquatic systems. Culture-independent interconnected meta-omic approaches are being increasing used to fill the gaps that classical microbial approaches cannot resolve. Here, we provide a prospective view of the increasing application of these high-throughput molecular technologies to evaluate the structure and functional activity of microbial communities in response to changes and disturbances in the environment, mostly of anthropogenic origin. Some relevant topics are reviewed, such as: (i) the use of microorganisms for water quality assessment, highlighting the incidence of antimicrobial resistance as an increasingly serious threat to global public health; (ii) the crucial role of microorganisms and their complex relationships with the ongoing climate change, and other stress threats; (iii) the responses of the environmental microbiome to extreme pollution conditions, such as acid mine drainage or oil spills. Moreover, protists and viruses, due to their huge impacts on the structure of microbial communities, are emerging candidates for the assessment of aquatic environmental health.

An integrated, modular approach to data science education in microbiology

We live in an increasingly data-driven world, where high-throughput sequencing and mass spectrometry platforms are transforming biology into an information science. This has shifted major challenges in biological research from data generation and processing to interpretation and knowledge translation. However, postsecondary training in bioinformatics, or more generally data science for life scientists, lags behind current demand. In particular, development of accessible, undergraduate data science curricula has the potential to improve research and learning outcomes as well as better prepare students in the life sciences to thrive in public and private sector careers. Here, we describe the Experiential Data science for Undergraduate Cross-Disciplinary Education (EDUCE) initiative, which aims to progressively build data science competency across several years of integrated practice. Through EDUCE, students complete data science modules integrated into required and elective courses augmented with coordinated cocurricular activities. The EDUCE initiative draws on a community of practice consisting of teaching assistants (TAs), postdocs, instructors, and research faculty from multiple disciplines to overcome several reported barriers to data science for life scientists, including instructor capacity, student prior knowledge, and relevance to discipline-specific problems. Preliminary survey results indicate that even a single module improves student self-reported interest and/or experience in bioinformatics and computer science. Thus, EDUCE provides a flexible and extensible active learning framework for integration of data science curriculum into undergraduate courses and programs across the life sciences.

Non-denitrifier nitrous oxide reductases dominate marine biomes

Microbial enzymes often occur as distinct variants that share the same substrate but differ in substrate affinity, sensitivity to environmental conditions, or phylogenetic ancestry. Determining where variants occur in the environment helps identify thresholds that constrain microbial cycling of key chemicals, including the greenhouse gas nitrous oxide (N₂O). To understand the enzymatic basis of N₂O cycling in the ocean, we mined metagenomes to characterize genes encoding bacterial nitrous oxide reductase (NosZ) catalyzing N₂O reduction to N₂. We examined data sets from diverse biomes but focused primarily on those from oxygen minimum zones where N₂O levels are often elevated. With few exceptions, marine nosZ data sets were dominated by 'atypical' clade II gene variants. Atypical nosZ has been associated with low oxygen, enhanced N₂O affinity, and organisms lacking enzymes for complete denitrification, i.e., non-denitrifiers. Atypical nosZ often occurred in metagenome-assembled genomes (MAGs) with nitrate or nitrite respiration genes, although MAGs with genes for complete denitrification were rare. We identified atypical nosZ in several taxa not previously associated with N₂O consumption, in addition to known N₂O-associated groups. The data suggest that marine environments generally select for high N₂O-scavenging ability across diverse taxa and have implications for how N₂O concentration may affect N₂O removal rates.

Heterotrophic Thaumarchaea with Small Genomes Are Widespread in the Dark Ocean

The Thaumarchaeota is a diverse archaeal phylum comprising numerous lineages that play key roles in global biogeochemical cycling, particularly in the ocean. To date, all genomically characterized marine thaumarchaea are reported to be chemolithoautotrophic ammonia oxidizers. In this study, we report a group of putatively heterotrophic marine thaumarchaea (HMT) with small genome sizes that is globally abundant in the mesopelagic, apparently lacking the ability to oxidize ammonia. We assembled five HMT genomes from metagenomic data and show that they form a deeply branching sister lineage to the ammonia-oxidizing archaea (AOA). We identify this group in metagenomes from mesopelagic waters in all major ocean basins, with abundances reaching up to 6% of that of AOA. Surprisingly, we predict the HMT have small genomes of ∼1 Mbp, and our ancestral state reconstruction indicates this lineage has undergone substantial genome reduction compared to other related archaea. The genomic repertoire of HMT indicates a versatile metabolism for aerobic chemoorganoheterotrophy that includes a divergent form III-a RuBisCO, a 2M respiratory complex I that has been hypothesized to increase energetic efficiency, and a three-subunit heme-copper oxidase complex IV that is absent from AOA. We also identify 21 pyrroloquinoline quinone (PQQ)-dependent dehydrogenases that are predicted to supply reducing equivalents to the electron transport chain and are among the most highly expressed HMT genes, suggesting these enzymes play an important role in the physiology of this group. Our results suggest that heterotrophic members of the Thaumarchaeota are widespread in the ocean and potentially play key roles in global chemical transformations.IMPORTANCE It has been known for many years that marine Thaumarchaeota are abundant constituents of dark ocean microbial communities, where their ability to couple ammonia oxidation and carbon fixation plays a critical role in nutrient dynamics. In this study, we describe an abundant group of putatively heterotrophic marine Thaumarchaeota (HMT) in the ocean with physiology distinct from those of their ammonia-oxidizing relatives. HMT lack the ability to oxidize ammonia and fix carbon via the 3-hydroxypropionate/4-hydroxybutyrate pathway but instead encode a form III-a RuBisCO and diverse PQQ-dependent dehydrogenases that are likely used to conserve energy in the dark ocean. Our work expands the scope of known diversity of Thaumarchaeota in the ocean and provides important insight into a widespread marine lineage.

Examining the predictive accuracy of metabolomics for small-for-gestational-age babies: a systematic review

INTRODUCTION

To date, there is no robust enough test to predict small-for-gestational-age (SGA) infants, who are at increased lifelong risk of morbidity and mortality.

OBJECTIVE

To determine the accuracy of metabolomics in predicting SGA babies and elucidate which metabolites are predictive of this condition.

DATA SOURCES

Two independent researchers explored 11 electronic databases and grey literature in February 2018 and November 2018, covering publications from 1998 to 2018. Both researchers performed data extraction and quality assessment independently. A third researcher resolved discrepancies.

STUDY ELIGIBILITY CRITERIA

Cohort or nested case-control studies were included which investigated pregnant women and performed metabolomics analysis to evaluate SGA infants. The primary outcome was birth weight <10th centile-as a surrogate for fetal growth restriction-by population-based or customised charts.

STUDY APPRAISAL AND SYNTHESIS METHODS

Two independent researchers extracted data on study design, obstetric variables and sampling, metabolomics technique, chemical class of metabolites, and prediction accuracy measures. Authors were contacted to provide additional data when necessary.

RESULTS

A total of 9181 references were retrieved. Of these, 273 were duplicate, 8760 were removed by title or abstract, and 133 were excluded by full-text content. Thus, 15 studies were included. Only two studies used the fifth centile as a cut-off, and most reports sampled second-trimester pregnant women. Liquid chromatography coupled to mass spectrometry was the most common metabolomics approach. Untargeted studies in the second trimester provided the largest number of predictive metabolites, using maternal blood or hair. Fatty acids, phosphosphingolipids and amino acids were the most prevalent predictive chemical subclasses.

CONCLUSIONS AND IMPLICATIONS

Significant heterogeneity of participant characteristics and methods employed among studies precluded a meta-analysis. Compounds related to lipid metabolism should be validated up to the second trimester in different settings.

PROSPERO REGISTRATION NUMBER

CRD42018089985.

Diverse Marinimicrobia bacteria may mediate coupled biogeochemical cycles along eco-thermodynamic gradients

Microbial communities drive biogeochemical cycles through networks of metabolite exchange that are structured along energetic gradients. As energy yields become limiting, these networks favor co-metabolic interactions to maximize energy disequilibria. Here we apply single-cell genomics, metagenomics, and metatranscriptomics to study bacterial populations of the abundant "microbial dark matter" phylum Marinimicrobia along defined energy gradients. We show that evolutionary diversification of major Marinimicrobia clades appears to be closely related to energy yields, with increased co-metabolic interactions in more deeply branching clades. Several of these clades appear to participate in the biogeochemical cycling of sulfur and nitrogen, filling previously unassigned niches in the ocean. Notably, two Marinimicrobia clades, occupying different energetic niches, express nitrous oxide reductase, potentially acting as a global sink for the greenhouse gas nitrous oxide.

Monitoring microbial responses to ocean deoxygenation in a model oxygen minimum zone

Today in Scientific Data, two compendia of geochemical and multi-omic sequence information (DNA, RNA, protein) generated over almost a decade of time series monitoring in a seasonally anoxic coastal marine setting are presented to the scientific community. These data descriptors introduce a model ecosystem for the study of microbial responses to ocean deoxygenation, a phenotype that is currently expanding due to climate change. Public access to this time series information is intended to promote scientific collaborations and the generation of new hypotheses relevant to microbial ecology, biogeochemistry and global change issues.

A compendium of geochemical information from the Saanich Inlet water column

Extensive and expanding oxygen minimum zones (OMZs) exist at variable depths in coastal and open ocean waters. As oxygen levels decline, nutrients and energy are increasingly diverted away from higher trophic levels into microbial community metabolism, resulting in fixed nitrogen loss and production of climate active trace gases including nitrous oxide and methane. While ocean deoxygenation has been reported on a global scale, our understanding of OMZ biology and geochemistry is limited by a lack of time-resolved data sets. Here, we present a historical dataset of oxygen concentrations spanning fifty years and nine years of monthly geochemical time series observations in Saanich Inlet, a seasonally anoxic fjord on the coast of Vancouver Island, British Columbia, Canada that undergoes recurring changes in water column oxygenation status. This compendium provides a unique geochemical framework for evaluating long-term trends in biogeochemical cycling in OMZ waters.