Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mBio
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Ecological and Evolutionary Science

Atribacteria Reproducing over Millions of Years in the Atlantic Abyssal Subseafloor

Aurèle Vuillemin, Sergio Vargas, Ömer K. Coskun, Robert Pockalny, Richard W. Murray, David C. Smith, Steven D’Hondt, William D. Orsi
Nicole Dubilier, Editor
Aurèle Vuillemin
aDepartment of Earth and Environmental Sciences, Paleontology and Geobiology, Ludwig-Maximilians-Universität München, Munich, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sergio Vargas
aDepartment of Earth and Environmental Sciences, Paleontology and Geobiology, Ludwig-Maximilians-Universität München, Munich, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ömer K. Coskun
aDepartment of Earth and Environmental Sciences, Paleontology and Geobiology, Ludwig-Maximilians-Universität München, Munich, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Robert Pockalny
cGraduate School of Oceanography, University of Rhode Island, Narragansett, Rhode Island, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Richard W. Murray
dWoods Hole Oceanographic Institution, Woods Hole, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
David C. Smith
cGraduate School of Oceanography, University of Rhode Island, Narragansett, Rhode Island, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Steven D’Hondt
cGraduate School of Oceanography, University of Rhode Island, Narragansett, Rhode Island, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
William D. Orsi
aDepartment of Earth and Environmental Sciences, Paleontology and Geobiology, Ludwig-Maximilians-Universität München, Munich, Germany
bLMU
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for William D. Orsi
Nicole Dubilier
Max Planck Institute for Marine Microbiology
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mBio.01937-20
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

How microbial metabolism is translated into cellular reproduction under energy-limited settings below the seafloor over long timescales is poorly understood. Here, we show that microbial abundance increases an order of magnitude over a 5 million-year-long sequence in anoxic subseafloor clay of the abyssal North Atlantic Ocean. This increase in biomass correlated with an increased number of transcribed protein-encoding genes that included those involved in cytokinesis, demonstrating that active microbial reproduction outpaces cell death in these ancient sediments. Metagenomes, metatranscriptomes, and 16S rRNA gene sequencing all show that the actively reproducing community was dominated by the candidate phylum “Candidatus Atribacteria,” which exhibited patterns of gene expression consistent with fermentative, and potentially acetogenic, metabolism. “Ca. Atribacteria” dominated throughout the 8 million-year-old cored sequence, despite the detection limit for gene expression being reached in 5 million-year-old sediments. The subseafloor reproducing “Ca. Atribacteria” also expressed genes encoding a bacterial microcompartment that has potential to assist in secondary fermentation by recycling aldehydes and, thereby, harness additional power to reduce ferredoxin and NAD+. Expression of genes encoding the Rnf complex for generation of chemiosmotic ATP synthesis were also detected from the subseafloor “Ca. Atribacteria,” as well as the Wood-Ljungdahl pathway that could potentially have an anabolic or catabolic function. The correlation of this metabolism with cytokinesis gene expression and a net increase in biomass over the million-year-old sampled interval indicates that the “Ca. Atribacteria” can perform the necessary catabolic and anabolic functions necessary for cellular reproduction, even under energy limitation in millions-of-years-old anoxic sediments.

IMPORTANCE The deep subseafloor sedimentary biosphere is one of the largest ecosystems on Earth, where microbes subsist under energy-limited conditions over long timescales. It remains poorly understood how mechanisms of microbial metabolism promote increased fitness in these settings. We discovered that the candidate bacterial phylum “Candidatus Atribacteria” dominated a deep-sea subseafloor ecosystem, where it exhibited increased transcription of genes associated with acetogenic fermentation and reproduction in million-year-old sediment. We attribute its improved fitness after burial in the seabed to its capabilities to derive energy from increasingly oxidized metabolites via a bacterial microcompartment and utilize a potentially reversible Wood-Ljungdahl pathway to help meet anabolic and catabolic requirements for growth. Our findings show that “Ca. Atribacteria” can perform all the necessary catabolic and anabolic functions necessary for cellular reproduction, even under energy limitation in anoxic sediments that are millions of years old.

INTRODUCTION

Marine sediments contain a ubiquitous “deep biosphere” (1) extending at least as far as 2,500 m below the seafloor (mbsf) (2), which consists of active and dormant cells (3–5) with measurable impacts on subseafloor biogeochemical processes (6). At abyssal water depths in the deep sea, the subseafloor communities generally are less sampled (7) than those in continental shelf sediments that have higher activities and rates of microbial sulfate reduction (8–10). At abyssal depths under the oligotrophic ocean gyres, the sedimentation rates are low, ranging from 1 to 5 m of sediment deposited per million years (11, 12). These abyssal subseafloor communities have extremely low metabolic activity (13) and live near the energy limit to life (6, 14). As a result of this, the deep biosphere of marine sediment is generally characterized by net death, and it remains poorly understood to what extent microbial activity is translated into cellular reproduction in the deep subseafloor (15, 16).

Many subseafloor microbes exhibit viability, since they actively take up carbon and nitrogen in incubation experiments (4, 17), indicating potential for microbial growth in energy-limited anoxic subseafloor sediments. Microbial activities can also be stimulated at redox interfaces deep below the seafloor over geological timescales (3, 18). However, the capacity of microbes to reproduce in abyssal subseafloor ecosystems close to the energy limit to life (6, 14) is particularly unconstrained, given the extreme scarcity of organic substrates in these settings (19). There is reason to suspect that cellular reproduction in the abyssal subseafloor is minimal, since microbial biomass tends to decrease an order of magnitude over the top 10 m of sediment in all abyssal locations yet sampled, reaching the detection limit for life at relatively shallow subseafloor depths of ca. 15 mbsf (12). This follows the global trend whereby subseafloor microbes tend to die faster than they grow (1), particularly in the top 10 m of marine sediment.

Here, we report an exception to this global trend in anoxic deep-sea clay recovered from an abyssal water depth of >5,500 m in the North Atlantic, characterized by an ultraslow sedimentation rate of ca. 3 m per million years. In contrast to oxic abyssal red clay where microbial abundance decreases several orders of magnitude over the top 10 mbsf (12, 20), we show here that microbial abundance in the anoxic abyssal clay increases an order of magnitude from the seafloor down to 15 mbsf (spanning ca. 5 million years). We then proceeded to use metatranscriptomics to further investigate the anaerobic metabolic mechanisms that explain this net growth in the size of the subseafloor microbial ecosystem over multimillion-year timescales.

RESULTS AND DISCUSSION

Sediment biogeochemistry.We obtained deep-sea clay sediment from a 5,515-m water depth in the ultraoligotrophic open ocean of the North Atlantic. This coring site (KN223-15) is characterized by a mean sedimentation rate of ca. 3 m per million years (11). Samples ranged from 0.1 to 30 m below seafloor (mbsf). Given the mean sedimentation rate, the deepest sample has an approximate age of 9 to 10 million years. Oxygen and nitrate penetration into the sediment is restricted to the top millimeter of sediment, as they were detectable in the bottom water but below detection in the uppermost portion of the core at 0.02 and 0.03 mbsf, respectively (see Fig. S1 in the supplemental material). The abyssal sediments of the North Atlantic are typically oxic red clay that tend to have O2 penetrating many meters into the seafloor (11, 21), but the subseafloor microbial ecosystem sampled here is unique in the sense that the sediments are anoxic despite having ultraslow sedimentation rates. Moreover, while the sediment of our abyssal subseafloor core displays sulfate (SO42−) concentrations of approximately 29 mM at 0.02 mbsf and is fully anoxic downward, the rates of anaerobic microbial SO42− reduction over the top 10 m of sampled sediment are low (−3.8 × 10−4 mol SO42− m3 yr−1) and below detection underneath, resulting in pore water SO42− remaining >20 mM throughout the core (Fig. 1A). This profile is very similar to profiles observed previously in anoxic sediments from other oligotrophic regions, such as the Eastern Equatorial Pacific, where sediments are anoxic but community metabolic activity is too slow to consume all of the available SO42− (6, 22).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Biogeochemistry, microbial diversity, and abundance in the anoxic sediment of North Atlantic site KN223-15. (A) Profiles of mean net SO42− reduction rates and dissolved SO42− concentrations. (B) Quantitative PCR (qPCR) of 16S rRNA (black dots) and dsrB (orange stars) genes, and the summed qPCR-normalized abundance of 16S rRNA gene OTUs affiliated with “Ca. Atribacteria” (blue), Chloroflexi (green), and Deltaproteobacteria (light orange). Error bars correspond to standard deviations (2 σ) based on four biological replicates. (C) Diversity of 16S rRNA genes based on three to four biological replicates. (D) Number of ORFs attributed to “Ca. Atribacteria” in the metatranscriptomes as a function of depth. Light blue circles are the number of ORFs detected in individual metatranscriptome replicates, and the larger dark blue circles are the total number of unique ORFs detected when summing across all replicates. The top two samples in panels B and C (0.1 to 0.3 mbsf) were recovered by gravity coring, the deeper samples (0.5 to 29 mbsf) were recovered via long-piston coring. Note in panels B and C that an increase in abundance of “Ca. Atribacteria” in the upper 1 mbsf coincides with their higher level of gene expression in panel D.

FIG S1

Map of sample site 15 with bathymetry, sedimentation rates, and pore water O2 profiles. (Left) Global map of subseafloor sedimentation rates and corresponding O2 penetration is modified from D’Hondt et al. (11). (Right) Note that at site 11, there is O2 penetrating to 30 mbsf, whereas at site 15 (the core sampled in this study), O2 is consumed immediately at the seafloor surface and the entire sediment sequence is anoxic. Download FIG S1, JPG file, 0.8 MB.
Copyright © 2020 Vuillemin et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Abundance and diversity of the microbial communities.The density of 16S rRNA gene copies per gram of wet sediment exhibits a subsurface peak in abundance, increasing 1 order of magnitude (2.5 × 105 to 2.3 × 106 copies) from 0.2 to 3.5 mbsf and remaining >106 between 3 and 10 mbsf. Thereafter, microbial 16S rRNA gene abundances decrease gradually to a minimum of 9.9 × 104 copies at 17 mbsf (Fig. 1B). In the uppermost sediment samples (0.1 to 0.4 mbsf), Actinobacteria, Planctomycetes, Deltaproteobacteria, and Gammaproteobacteria dominate the community (Fig. 1C). Below this depth, starting at 0.5 mbsf, the relative abundance of “Candidatus Atribacteria” rapidly increased from <5% to >40% (Fig. 1B and C) in all four biological replicates sampled (see Fig. S2 and S3). The “Ca. Atribacteria” is only represented by 3 operational taxonomic units (OTUs), with one single OTU accounting for up to 40% of the whole community throughout the record (Fig. 2A and B), which was consistent across three to four biological replicates (Fig. S2). This dominant OTU (Fig. 2A) was most closely affiliated with unpublished 16S rRNA gene sequences from deep subseafloor sediments from the Nankai Trough (23) and a potential gas hydrate region of southwest Taiwan subseafloor sediments (Fig. 2B).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Relative abundance and phylogenetic analyses of the most abundant OTU (“Ca. Atribacteria”) and its phylogenomic protein markers. (A) The histogram shows the relative abundance of the most abundant OTU (“Ca. Atribacteria”), error bars are standard deviations from 3 to 4 biological replicates (see Fig. S2 in the supplemental material). (B) Phylogenetic analysis of 16S rRNA genes (V4 hypervariable region) showing the affiliation of all three “Ca. Atribacteria” OTUs to the JS1 clade. (C) Phylogenetic tree based on 31 concatenated protein markers identified with AMPHORA2 (33). The monophyletic clade including our samples in the tree has >95% similarity to the sister clade with taxa from the Atlantic Ocean, Ross Sea, and North Pacific methane hydrate subseafloor sediment.

FIG S2

Microbial diversity in the North Atlantic abyssal anoxic subseafloor. Diversity of 16S rRNA genes obtained from the sequencing of four biological replicates. Download FIG S2, JPG file, 0.8 MB.
Copyright © 2020 Vuillemin et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S3

Phylogenetic analysis of 16S rRNA genes (V4 hypervariable region). The tree shows the 30 most abundant operational taxonomic units (OTUs) detected in this study. The pie chart shows the 16S rRNA gene richness of the whole assemblage across all depths. Download FIG S3, JPG file, 1.0 MB.
Copyright © 2020 Vuillemin et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

According to a protocol for quantitative normalization of barcoded 16S gene sequence data (24), we normalized group-specific 16S rRNA gene abundances by dividing the total number of 16S rRNA gene copies determined via quantitative PCR (qPCR) by their affiliated 16S rRNA gene relative abundances (fractional percentage of total sequences per sample). This showed that “Ca. Atribacteria” dominates the community (Fig. 1B), which is attributed to a single 16S OTU (Fig. 2). We considered the potential influence of multiple copies of 16S rRNA genes (25) on the result and searched for all “Ca. Atribacteria” genomes in the JGI database sequenced to date, which revealed that all sequenced genomes from this group have only one 16S rRNA gene copy, whereas Deltaproteobacteria and Chloroflexi have a median of two 16S rRNA gene copies (up to four) per genome (25). Thus, despite having an average lower 16S rRNA gene copy number than Chloroflexi and Deltaproteobacteria, “Ca. Atribacteria” relative abundance was high in our 16S rRNA gene sequence data set, even suggesting that our data actually underestimate the abundance of “Ca. Atribacteria.” These profiles show exponential net increases of “Ca. Atribacteria” and Chloroflexi within the upper 10 m of sediment, whereas those of Deltaproteobacteria increase at 0.4 mbsf but recede rapidly over time from 104 to 103 copies or less, reaching our limit of detection (102 gene copies g−1 [wet weight] sediment) at the bottom of the core. The vertical profile of the beta subunit of the dissimilatory sulfate reductase (dsrB) gene runs parallel to that of Deltaproteobacteria 16S rRNA genes, which points to Deltaproteobacteria as the main sulfate-reducing bacteria (SRB) in this anoxic sediment. This inference is supported by the abrupt 10-fold drop in dsrB and Deltaproteobacteria 16S rRNA gene sequences at 10 mbsf (Fig. 1B), the same depth at which net SO42− reduction decreases to below detection levels (Fig. 1A). The abundance of SRB is lower than in organic-rich shelf sediments (26, 27), which may be related to the relatively low sedimentation rate at our sampling location, since SO42− reduction rates are correlated with the square of the sedimentation rate (9). The relative slowness of SO42− reduction at our sampling location is readily evident in the profile of dissolved SO42− concentrations, which remained >20 mM throughout the core (Fig. 1A).

Metagenomic analysis.Metagenomes from five depths were sequenced at an average depth of 8.4 million reads (± 2.5 million) (see Table S3). Because the DNA was amplified by multiple displacement amplification (MDA), we only considered presence/absence of genes and did not interpret the relative abundance of genes in the MDA-amplified metagenomes to be indicative of taxon abundance. De novo “binning”-based methods for creating metagenome-assembled genomes (MAGs) are useful for discovering new taxa when distantly related genomes in databases preclude similarity-based searches (28–31). Because of the ultralow DNA concentrations extractable from these abyssal clay sediments, we were only able to obtain metagenomic data after amplifying the extracted DNA using multiple displacement amplification (MDA). Our attempts at de novo binning of the MDA products revealed a selective amplification of short fragments that precluded binning and completion of high-quality MAGs (see Table S2). Specifically, manual curation of MaxBin results using Anvi’o (32) produced 12 bins at relatively low levels of genome completeness (17% to 44%) that were able to be assigned to “Ca. Atribacteria” (see Fig. S4). Thus, in order to increase the annotation of putative functions of open reading frames (ORFs) in metagenomes from the “Ca. Atribacteria” that could not be assembled into bins, we also applied a bioinformatics pipeline, whereby ORFs harbored by de novo assembled contigs were searched for similarity against a large aggregated genome database of predicted proteins from all atribacterial MAGs and single-cell genomes (SAGs) sequenced to date, including all published data from subsurface metagenomics studies (see Materials and Methods). We then extracted all the ORFs having a predicted protein from a “Ca. Atribacteria” genome as a top BLASTp hit and ran a phylogenomic analysis based on 31 phylogenomic markers using AMPHORA (33). This phylogenomic analysis demonstrates that our previously published method (34, 35) recovers ORFs with high similarity (>95% amino acid similarity) to predicted proteins in previously sequenced “Ca. Atribacteria” genomes (Fig. 2C). The phylogenomic analysis shows that the atribacterial ORFs from all sampled depths form a monophyletic clade sister to those from deep-sea sediments of the Atlantic Ocean, Ross Sea, and Pacific Ocean methane hydrates within the JS1 clade (Fig. 2C). This high level of similarity to existing “Ca. Atribacteria” genomes enabled similarity-based assignment of “Ca. Atribacteria” ORFs in our samples (see Fig. S5). As further evidence of this, the median similarity of ORFs to their top hit in the database in all metatranscriptome and metagenome data sets of >60% indicates that the ORFs in our samples had relatively high similarity to existing predicted proteins in the database (see Fig. S6).

FIG S4

Circular representation of refined binning results obtained with DAS Tool for all bins (top) and MaxBin for “Ca. Atribacteria” only (bottom). For detailed results, see Table S2. Download FIG S4, JPG file, 0.5 MB.
Copyright © 2020 Vuillemin et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S5

Affiliation of ORFs encoded in metagenomes (A) and metatranscriptomes (B) based on similarity. The pie charts display the relative taxonomic abundances assigned to metagenomes and transcriptomes with sediment depth, which show consistency with the 16S rRNA gene assemblages. Download FIG S5, JPG file, 0.6 MB.
Copyright © 2020 Vuillemin et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S6

Box and whisker plot of the similarity of ORFs encoded in the metatranscriptomes and metagenomes compared to their closest top hit in the database after BLASTp searches with DIAMOND. In the box plots, the thick horizontal lines represent the median values, the ends of the boxes represent the upper and lower quartiles (25% data greater or less than this value), and the vertical lines extending from the boxes indicate the maximum and minimum values. Individual points represent outliers. The individual replicates are plotted separately for the metatranscriptomes. The dashed horizontal line indicates the threshold that we used to assign an ORF to a particular group (e.g., “Ca. Atribacteria”). Download FIG S6, JPG file, 0.6 MB.
Copyright © 2020 Vuillemin et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

We confirmed the accuracy of this previously published approach (34, 35) by an in silico test for true- and false-positive annotations based on 151 randomly selected peptide fragments extracted from “Ca. Atribacteria” predicted proteomes as well as other bacterial and archaeal genomes (see Fig. S7). In this analysis, 50 randomly selected predicted proteins were randomly cut into peptide fragments ranging from 20 to 140 amino acid residues in length in order to replicate partial ORFs typically recovered in metagenomes. The random peptide fragments were then searched against the large aggregate database (34, 35) for their top hits with BLASTp. This in silico test showed that 100% of all randomly cut peptide fragments from atribacteria were true positives; they had the same atribacterial genome as a top BLASTp hit. This shows that our use of a similarity-based approach for ORF annotations (34–36) is adequate for assigning ORFs harbored by de novo assembled contigs to groups at high taxonomic levels, including those derived from the candidate phylum Atribacteria.

FIG S7

An in silico experiment testing the rate of true-positive and false-negative annotations of 151 simulated ORFs (randomly selected peptide fragments) from the predicted proteomes of a cultivated bacterium and three metagenomic-assembled genomes. (A) Most (89%) of the simulated ORFs (n = 151 randomized peptide fragments) were correctly annotated to the species level (“Thalassospira xiamenensis”), and 93% of simulated ORFs were correctly annotated to the genus level when comparing against MetaProt. The same experiment applied to the predicted proteomes of the archaeon “Ca. Prometheoarchaeum syntrophicum” (B), Chloroflexi (C), and “Ca. Atribacteria” (D) MAGs shows that the false-negative rate is much lower than that for the bacterial genome: only 1 of 150 peptide fragments was incorrectly assigned. The red and blue colors in the histogram plots represent those random peptide fragments that were assigned as false negatives and true positives, respectively: e.g., those that had highest similarity (best BLASTp hit) to an organism other than “T. xiamenensis”, “Ca. Prometheoarchaeum syntrophicum,” Chloroflexi, and “Ca. Atribacteria” in the MetaProt database (e.g., the individual thin slices in the pie chart figure above). Download FIG S7, JPG file, 0.9 MB.
Copyright © 2020 Vuillemin et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Gene expression analysis.Metatranscriptomes from eight depths were produced in biological replicates (two to seven replicates per depth) (Fig. 1D) and sequenced at an average depth of 4.0 million reads (±1.5 million) with 91,199 contigs across all samples sequenced (Table S3). Similar to that for the 16S rRNA gene abundances, there was a subsurface peak in the number of unique expressed ORFs assigned to “Ca. Atribacteria” that increases exponentially between 3.5 and 10 mbsf, which was consistent across multiple replicate metatranscriptomes from each depth, and at 16 m, our RNA limit of detection was reached (Fig. 1D). We defined this depth as our RNA detection limit, because the number of unique ORFs assigned to “Ca. Atribacteria” was no longer detectable and, below this depth, the only ORFs annotated were those assigned to groups of known contaminants from molecular kits, including those from human skin and soil (37). Many of these same groups are common laboratory contaminants found in dust samples from our lab in 16S rRNA gene surveys (38) and include Pseudomonas, Rhizobium, Acinetobacter, and Staphylococcus. Presumably, the detection limit was reached because a smaller amount of extracted RNA from the in situ active community becomes overprinted by background “noise” from contaminating DNA (or RNA) derived from the kits, aerosols, or other laboratory contaminants.

There was a statistically significant correlation (r values = 0.59 and 0.61, P values = 0.016 and 0.014) between the abundance of atribacterial 16S rRNA genes and expressed ORFs (see Table S1). The number of expressed ORFs correlated with 16S rRNA gene quantities, which was consistent for the entire data set (total bacteria) and when comparing the number of ORFs expressed per group to the qPCR normalized abundance of the same taxonomic group for the four groups with the highest numbers of detected ORFs in the metatranscriptomes (atribacteria, Deltaproteobacteria, Chloroflexi, and Archaea) (Fig. 3A).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Group-specific transcriptional activity correlates with abundance. (A) Abundance of 16S rRNA genes (x axis) normalized to total bacteria (white), “Ca. Atribacteria” (blue), Chloroflexi (green), Deltaproteobacteria (light orange), and Archaea (red) 16S rRNA sequences plotted against their respective number of ORFs assigned to the same groups in the metatranscriptomes (y axis). Error bars represent standard deviations across biological replicates (in the case of three or more) or ranges (in the case of two replicates). (B) Abundance of 16S rRNA genes normalized to “Ca. Atribacteria” (blue) and Chloroflexi (green) plotted against their respective number of ORFs in the transcripts assigned to cell division (FtsAEKQWZ, MreBC, and RodA), showing that cell division is only detectable at 16S rRNA gene densities of >105 copies.

TABLE S1

Statistical correlation between the abundance of 16S rRNA genes and the number of ORFs expressed in the metatranscriptomes. Download Table S1, JPG file, 0.2 MB.
Copyright © 2020 Vuillemin et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S2

Binning results using different algorithms. (A) Summary from DAS Tool workflow with which we produced 6 taxonomically unassigned bins. (B) We then obtained 14 bins from manual curation of DAS Tool results. (C) Summary of MaxBin results after manual curation, which allowed us to assign 16 bins to “Ca. Atribacteria,” 4 of which have 0.0% genome completeness and 12 have 17% to 44% genome completeness. Download Table S2, JPG file, 2.2 MB.
Copyright © 2020 Vuillemin et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S3

Sequencing and assembly statistics for the anoxic subseafloor metagenomes and transcriptomes. Download Table S3, XLSX file, 0.1 MB.
Copyright © 2020 Vuillemin et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Evidence for subseafloor reproduction.The steadily increasing abundance of “Ca. Atribacteria” over time (since sediment deposition) is apparently due to cells undergoing cytokinesis, as the 16S rRNA gene abundances correlated with the transcription of genes encoding proteins involved in cell division and cell shape determination such as FtsAEKQWZ, MreBC, and RodA (Fig. 3B) (39, 40). The expressed genes encoding the Fts proteins form the “divisome” (41), which includes the FtsZ ring, a cytokinetic protein ring that localizes at the cell division site prior to cytokinesis in dividing bacterial cells and pumps the replicated chromosome into the daughter cell (42). These genes are expressed during cellular division in “Ca. Atribacteria,” indicating the exponential growth phase (43). Thus, while we do not report measurements of biomass turnover such as amino acid racemization (44, 45) or stable isotope probing (4, 17), the steadily increasing microbial biomass of “Ca. Atribacteria” with increasing sediment depth together with the detected correlation between their abundance and transcriptional activity (Fig. 3A) and the transcription of ORFs encoding FtsZ ring proteins (Fig. 3B) strongly suggest that the higher abundances of “Ca. Atribacteria” between 0.5 and 10 mbsf are due to actively dividing cells. In this zone of apparent increased reproduction, the number of unique expressed ORFs correlated significantly with the number of 16S rRNA gene copies from the dominant groups (Fig. 1 and 3). The most parsimonious explanation for this is that more expressed ORFs from atribacteria are detected at these depths because there are higher numbers of metabolically active atribacteria producing mRNA transcripts.

The sediment ages in this interval span several million years; thus, this higher number of metabolically active microbes originated from “Ca. Atribacteria” that slowly reproduced through binary fission and cytokinesis over multimillion-year timescales. Because genes involved in the formation of the divisome, FstZ ring, and cytokinesis were only expressed at depths where the highest 16S rRNA gene copies were found (Fig. 4B), the higher number of 16S rRNA gene copies at those depths is likely due to higher numbers of actively reproducing atribacterial cells. Active cells have higher numbers of ribosomes and thus higher copies of 16S rRNA per cell (46), but our qPCR assay targeted DNA (the bacterial chromosome), not RNA (expressed genes), thus discarding active rRNA synthesis as a confounding factor in our bacterial abundance estimations. Since “Ca. Atribacteria” has only one 16S rRNA gene copy, it is possible to conclude that the observed 10-fold increase in atribacterial 16S rRNA gene copy numbers over the top 10 mbsf requires chromosomal (genome) replication. In bacteria, genome replication occurs in actively dividing cells prior to cytokinesis and binary fission (46), assuming steady-state microbial input in sediment over millions of years. The only alternative explanation to our observations would require an increase in 16S rRNA copy number in the chromosome of “Ca. Atribacteria” over the top 10 mbsf followed by its subsequent decrease in the same chromosome below 10 mbsf. Such an incredibly high rate of genome evolution affecting the highly conserved 16S gene and then acting in a reversible way after 10 mbsf is inconceivable, especially considering this gene evolves at a rate of 1% every 100 million years (47) and our deepest sampled depth is roughly 10 million years old. Thus, the most likely explanation is that the detection of more 16S rRNA gene copies from the dominant OTU of “Ca. Atribacteria” (Fig. 2A) is due to more cells, each containing one chromosome with one copy of the 16S rRNA gene. The correlation of these 16S rRNA genes with higher numbers of expressed genes from “Ca. Atribacteria” over the top 15 mbsf (Fig. 1D and 3) can be attributed to a single clade (Fig. 2) that has been slowly reproducing and increasing in abundance over millions of years.

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Metabolic potential and transcriptional activity of “Ca. Atribacteria” and its microcompartment. (A) Presence (colored) or absence (white) of ORFs assigned to “Ca. Atribacteria” encoding predicted proteins in metagenomes (circles) and metatranscriptomes (squares). (B) Bubble plots showing the coverage values (percent total reads) of expressed genes identified for “Ca. Atribacteria” at eight different depths. The numbering of the groups in panel A corresponds to the same numbering in panel B. DERA, 2-deoxy-d-ribose 5-phosphate aldolase. (C) Sequences of gene synteny related to the BMC present in the cell of “Ca. Atribacteria.” These three syntenies provide evidence for metabolic use of aldehydes and alcohols via dehydrogenases with biosynthesis of acetyl-coenzyme A (acetyl-CoA) and regeneration of NAD+.

Predicted metabolism for “Ca. Atribacteria.”Our finding that “Ca. Atribacteria” expressed the highest number of protein-encoding genes among all other groups (Fig. 3) is furthermore consistent with their high levels of gene transcription in deep subseafloor sediments of the Baltic Sea (48). Below 15 mbsf, the microbial ecosystem transitions to net death, since microbial abundances decrease by 2 orders of magnitude, RNA levels decline to below detection, and the net SO42− reduction rate is below our detection limit (Fig. 1). However, “Ca. Atribacteria” remains the dominant group throughout the core, even below 10 mbsf. Its continued dominance is consistent with its abundance in other deep subseafloor settings (49–52) and the findings that dominant taxa in the subseafloor community do not necessarily require cellular reproduction to outcompete other taxa but can reach higher relative abundances from lower mortality than their competitors over long timescales (53, 54). This is supported by the lack of expression of ORFs by atribacteria harboring cellular division proteins below 10 mbsf, where atribacterial 16S rRNA gene abundances are <105 copies per g (Fig. 3B). Similar to what has been predicted from genomic studies (55–58), the gene transcription data strongly indicate that the “Ca. Atribacteria” dominating throughout the core actively utilizes a sugar-based acetogenic metabolism (Fig. 4A). This is also consistent with gene transcription from “Ca. Atribacteria” in deep subseafloor Baltic Sea sediments that showed active utilization of trehalose (48).

The Wood-Ljungdahl pathway (WLP) has been found to be a widespread and potentially important metabolic pathway in subsurface microbial life (57, 59, 60). The WLP can be used catabolically to achieve redox balance and regenerate NAD+ and oxidized ferredoxin to thereby increase anaerobic metabolic efficiency (61). In acetogenic bacteria, this is coupled to the Rnf complex at the membrane that utilizes either ferredoxin or NAD+ as the terminal electron acceptor, driving a Na+ ion gradient and ATP synthesis (Fig. 5) (62). Several lines of evidence in our gene transcription data indicate that “Ca. Atribacteria” utilizes a similar metabolism.

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

Potential “Ca. Atribacteria” physiology based on metagenomes and transcriptomes. The figure is based on the known interlinking of energy metabolism pathways for (homo)acetogenic bacteria (61) and BMC metabolic functions (57), as well as the function of the FstZ ring and divisome (41). Colored dots correspond to the successive enzymatic steps (as listed in Fig. 4), e.g., fermentation corresponds to glycolysis (11 dots) and lactate fermentation (3 dots). Actively expressed genes encoding the divisome complex (i.e., FtsAEKQWZ) are indicated at the top right. nPPP, nonoxidative pentose phosphate pathway; SLP, substrate-level phosphorylation; DERA, 2-deoxy-d-ribose 5-phosphate aldolase; ALDH, aldehyde dehydrogenase; ADH, alcohol dehydrogenase; FBEB, flavin-based electron bifurcation; FHL, formate hydrogen lyase; WLP, Wood-Ljungdahl pathway.

Specifically, “Ca. Atribacteria” expressed transcripts encoding ORFs with similarity to proteins involved in sugar fermentation and the WLP (Fig. 4B). Genes encoding enzymes of the WLP were all present (Fig. 4A) and differently expressed by “Ca. Atribacteria” as were those encoding proteins involved in glycolysis, fermentation, and electron bifurcation (Fig. 4B). Transcription of ORFs by “Ca. Atribacteria” with similarity to enzymes known to be involved in electron bifurcation that results in production of molecular hydrogen included the beta subunit of the heterodisulfide reductase, formate hydrogen lyase (63), gamma subunit of the NiFe hydrogenase (64), and subunit alpha of the reversible formate dehydrogenase (65). The formate hydrogen lyase complex combines NiFe hydrogenase and soluble formate dehydrogenase to couple formate and/or CO-dependent hydrogen production to the generation of the Na+ motive force to generate ATP. This enzymatic complex allows growth on formate disproportionation (66), during which two electrons are passed to hydrogenase and two protons are neutralized (67). The gene expression data indicate that the “Ca. Atribacteria” produces the chemiosmotic Na+ gradient using the Rnf complex for ATP synthesis (61). Atribacterial ORFs encoding the Rnf complex were expressed at the majority of depths (Fig. 4B), indicating this is indeed an important mechanism of anaerobic ATP synthesis.

We also detected active transcription of genes encoding pyruvate-formate lyase, acetate kinase, and formate dehydrogenase, which is reversible (65); therefore, the WLP could be used for either catabolic or anabolic purposes (Fig. 5). However, the expression of an H4F:B12 methyltransferase from “Ca. Atribacteria” is a direct indication that the WLP is functioning as an acetogenic pathway as opposed to one for methanogenesis, because acetogens use 5-methyltetrahydrofolate:corrinoid (H4F:B12) methyltransferase, whereas methanogens use tetrahydromethanopterin methyltransferase in the last step of methyl synthesis prior to acetyl synthesis (68). Expression of the H4F:B12 methyltransferase from “Ca. Atribacteria” was relatively high in the deeper samples at 8.2 and 9.7 mbsf (Fig. 4B), where the steadily increasing abundance of 16S rRNA genes with depth reached their peak values (Fig. 1B).

“Ca. Atribacteria” expression of transporter-encoding genes for multiple types of sugar (i.e., hexose, hexulose, pentose, pentulose, and ribose) and glycosidases (Fig. 4B) indicates a potential preference for sugar-based substrates. Sugars are in anoxic sediments as an energy substrate from bacterial necromass (e.g., ribose-containing nucleic acids DNA and RNA) and could be used as a fermentation substrate (69). The only detectable gene expression from “Ca. Atribacteria” observed at a 15.9-m depth indicates sugar transport, peptidases, fructose biphosphate, and phosphogluconate aldolase activities (Fig. 4B), again pointing that utilization of necromass (sugars and peptides from dead cells) is an important activity for long-term survival.

Fermentation products and other oxidized organic substrates (i.e., alcohols, ketones, aldehydes, and carboxylic acids) become toxic if they accumulate in the cell (70, 71). The gene expression from “Ca. Atribacteria” shows a metabolic potential for aldehydes to be reoxidized in secondary fermentations via a bacterial microcompartment (BMC), regenerating additional reduced ferredoxin and NADH in the process (72, 73). We obtained three contigs from “Ca. Atribacteria” with ORFs encoding BMC shell proteins syntenous with Rnf, aldehyde dehydrogenase, and 2-deoxy-d-ribose 5-phosphate aldolase (DERA) (Fig. 4C), similar to the already reported genome assembled from “Ca. Atribacteria” (57). These genes were coexpressed at multiple depths (Fig. 4B), suggesting that the BMC and Rnf-based establishment of the Na+ ion gradient for ATP production are connected (Fig. 5). Collectively, these findings indicate that “Ca. Atribacteria” may use a BMC to recycle toxic intermediates produced during secondary fermentations for additional energy (Fig. 4 and 5). The ability to use a BMC to sustain secondary fermentations, in addition to primary fermentation of sugars (Fig. 4 and 5), may optimize energy conservation (e.g., increase ATP yields per mole substrate oxidized) under the extreme energy-limiting subseafloor conditions. Similar to our results, a dominance of “Ca. Atribacteria” in organic lean sediments has been observed previously during the IODP expedition 313 to the New Jersey shelf (27). However, we acknowledge that biases inherent to all methods for determining the abundance of Bacteria and Archaea make quantitative comparisons difficult (5) and thus interpret the qPCR data showing a relatively lower abundance of Archaea (Fig. 1B) with caution.

Although fermenters and acetogens are often outcompeted by SRB in organic-rich subseafloor sediments, the poor fitness of SRB in the sampled environment (Fig. 1B and C), compared to that of “Ca. Atribacteria,” may be related to the relatively low SO42− reduction rates (Fig. 1A) and relatively low rates of organic matter burial at our sampled site (11, 19). In contrast, Chloroflexi related to Dehalococcoidia, which are predicted (homo)acetogens with metabolic potential to degrade complex organic compounds (59, 60), are relatively abundant in these deep-sea anoxic clays (Fig. 1B and C) and transcribe genes involved in cell division (Fig. 3B). Our results lead us to speculate that organic-lean, anoxic abyssal subseafloor sediment thus represents a niche that favors the reproduction of Chloroflexi and “Ca. Atribacteria” over that of other microbes, including SRB. Moreover, in culture experiments from the first cultivated representative of the atribacteria, the presence of an H2-consuming methanogenic partner increased the growth rate of atribacteria by more than 100-fold (43). Based on this experimental result, we speculate that the relatively high abundance of the atribacteria in our sampled ancient anoxic sediment could be attributed to an as-of-yet-unidentified syntrophic or semisyntrophic partner organism, from the Chloroflexi, for example, Dehalococcoidia (74, 75).

Our data show that reproduction can occur in the subseafloor over multimillion-year timescales. Although RNA from the “Ca. Atribacteria” was no longer detectable below 16 mbsf, metagenomes and 16S sequencing show that they subsist and remain the dominant group down to at least 29 mbsf. The extractable DNA with our protocol does not target extracellular DNA (20); thus, the DNA in these deeper sediments likely does not derive from “dead” DNA preserved from once-living cells. Therefore, we assume that “Ca. Atribacteria” remains viable, and potentially active, in the deeper sediments but at cell abundances and transcriptional activities that are too low for our current RNA-based methods. Thus, the survival of the “Ca. Atribacteria,” associated with gene expression at relatively low levels, supports the hypothesis that reduced metabolic activity is a fitness advantage (16) in the energy-starved subseafloor.

MATERIALS AND METHODS

Sampling.All samples were taken by cruise KN223 of the R/V Knorr in the North Atlantic, from 26 October to 3 December 2014 (Woods Hole, MA, to Woods Hole, MA). At site KN223-15 (33°29.0′N, 54°10.0′W; water depth, 5,515 m), successively longer sediment cores were retrieved using a multicorer (∼0.4 m), gravity corer (∼3 m), and the Woods Hole Oceanographic Institution (WHOI) piston-coring device (∼28 m). Additional details of sampling are published (11). Dissolved oxygen concentrations in the core sections (see Fig. S1 in the supplemental material) were measured with optical O2 sensors as described previously (76). Concentrations of dissolved SO42− were measured as published previously (6), and SO42− reduction rates were calculated as previously described (14). Sediment subcores were retrieved on the ship aseptically using end-cut sterile syringes and kept frozen at −80°C until extraction in the home laboratory.

DNA extraction, quantitative PCR, and 16S rRNA gene libraries.For each sampled depth, we performed four biological replicates of DNA extraction. Total DNA was extracted from 0.7 g of sediment as previously described (20). DNA templates were diluted 10 times in ultrapure PCR water (Roche) and used in qPCR amplifications with updated 16S rRNA gene primer pair 515F (5′-GTG YCA GCM GCC GCG GTA A-3′) and 806R (5′-GGA CTA CNV GGG TWT CTA AT-3′) to increase our coverage of Archaea and marine clades (77) and run as previously described (38). All qPCRs were set up in 20-μl volumes with 4 μl of DNA template, 20 μl SsoAdvanced SYBR green Supermix (Bio-Rad), 4.8 μl nuclease-free H2O (Roche), 0.4 μl primers (10 μM; biomers.net), and 0.4 μl MgCl2 and carried out on a CFX-Connect qPCR machine for gene quantification. For 16S rRNA genes, we ran 40 PCR cycles of two steps corresponding to denaturation at 95°C for 15 s and annealing and extension at 55°C for 30 s. To measure the abundance of dsrB genes, we used a previously described assay (78) with the primer pair dsrB4-R (5′-GTG TAG CAG TTA CCG CA-3′) and dsrB2060F (5′-CAA CAT CGT YCA YAC CCA GGG-3′). All qPCRs were set up in 20-μl volumes with 4 μl of DNA template and performed as previously described (79). Gel purified amplicons of the 16S rRNA, dsrB, and mcrA genes were quantified in triplicates using QuantiT dsDNA reagent (Life Technologies) and used as a standard. An EpMotion 5070 automated liquid handler (Eppendorf) was used to set up all qPCRs and prepare the standard curve dilution series spanning from 107 to 101 gene copies. Reaction efficiency values in all qPCR assays were between 90% and 110% with R2 values of >0.95% for the standards.

For 16S rRNA gene library preparation, qPCR runs were performed with barcoded primer pair 515F and 806R. All 16S rRNA gene amplicons were purified from 1.5% agarose gels using the QIAquick gel extraction kit (Qiagen), quantified with the Qubit dsDNA HS assay kit (Thermo Fisher Scientific), normalized to 1 nM solutions, and pooled (38). Library preparation was carried out according to the MiniSeq System Denature and Dilute Libraries guide (Illumina). Sequencing was performed on all four biological replicates (Fig. S2) on the Illumina MiniSeq platform at the Geo-Bio LMU Center. We used USEARCH version 10.0.240 for MiniSeq read trimming and assembly, OTU picking, and 97% sequence identity clustering (80), which we showed previously captures an accurate diversity represented within mock communities sequenced on the same platform (38). OTU representative sequences were identified by BLASTn searches against SILVA database version 132 (81). To identify contaminants, 16S rRNA genes from extraction blanks and dust samples from the lab were also sequenced in triplicates (38). These 16S rRNA gene sequences were used to identify any contaminating bacteria (e.g., Acinetobacter, Bacillus, and Staphylococcus) and selectively curate the OTU table of our anoxic abyssal clay samples prior to downstream analysis. The 30 most abundant OTUs were aligned with SINA online v.1.2.11 (82) and plotted in a maximum likelihood RAxML phylogenetic tree (83) (Fig. S3) against the SILVA 16S rRNA SSU NR99 reference database version 132 (81) using ARB (84). Closest environmental sequences with nearly full-length sequences (>1,400 bp) were selected as taxonomic references and used to calculate trees implemented with the bacterial and archaeal filter and advanced bootstrap refinement selecting the best tree among 300 replicates (84). Partial OTU sequences were added to the tree using the maximum parsimony algorithm without allowing changes of tree topology (Fig. 2 and S3).

Metagenome libraries.Whole-genome amplifications were performed on DNA extracts at 10× dilutions through a multiple displacement amplification (MDA) step of 6 to 7 h, using the REPLI-g Midi kit (Qiagen) and according to the manufacturer’s instructions. Dilution of the extracted DNA was necessary, because the MDA reaction was inhibited in undiluted DNA extracts (presumably from coextracted inhibitory chemicals). We added SYBR green I (Invitrogen) at 1,000× concentration to visualize the amplification in real time on the CFX-Connect qPCR machine with a fluorescence reading taken every 10 min. Amplification was stopped after reaching the exponential increase by heating to 65°C for 3 min. Metagenomic libraries were prepared using the Nextera XT DNA Library Prep kit (Illumina), quantified on an Agilent 2100 bioanalyzer system (Agilent Genomics), and normalized with the Select-a-Size DNA Clean and Concentrator MagBead kit (Zymo Research) as previously described (20).

RNA extractions and metatranscriptome libraries.Total RNA extractions were obtained according to a previously published protocol (85). In brief, RNA was extracted from 3 g of sediments using the FastRNA Pro Soil-Direct kit (MP Biomedicals) according to the manufacturer’s instructions, with the addition of 4 μl glycogen (0.1 g ml−1) to increase yield during precipitation of the RNA pellet and final elution in 40 μl PCR-grade water (Roche). Extraction blanks (no RNA added) were processed alongside to assess laboratory contamination, and sequencing of these contamination controls revealed common laboratory contaminants found in dust samples from our lab in 16S rRNA gene surveys (38), including Pseudomonas, Rhizobium, Acinetobacter, and Staphylococcus. RNA extracts were quantified using the QuBit RNA HS assay kit (Thermo Fisher Scientific). DNase treatment, synthesis of cDNA, and library construction were performed on the same day from 10 μl of RNA templates using the Trio RNA-Seq kit protocol (NuGEN Technologies). Libraries were quantified as described above. All libraries were diluted to 1 nM and pooled for further sequencing on the MiniSeq platform (Illumina).

Gene identification and normalization in metagenomes.The SqueezeMeta (86) metagenomic analysis pipeline was used for downstream analysis of metagenomic reads in coassembly mode. For adapter removing, trimming, and quality filtering, we used Trimmomatic set to the following: leading, 8; trailing, 8; sliding window, 10:15; and minimum length, 30 (87). Contigs were assembled to minimum length of 200 bp using Megahit assembler (88). ORFs for genes and rRNAs were called using Prodigal (89), and rRNA genes were determined by barrnap (90). RDP classifier were used for the classification of 16S rRNA genes (91). Diamond software (92) was deployed for taxonomic assignment of retrieved gene homologies against the GenBank, eggNOG v. 4.5 (93), and KEGG (94) databases. Cutoff values for assigning hits to specific taxa were performed at an E value of 1 × e−3, minimum amino acid similarity of 40 for taxa and 30 for functional assignment, using SqueezeMeta with default settings. Reads were mapped onto contigs and genes by using Bowtie2 (95). Coverage and transcripts per million (TPM) values were calculated using SqueezeMeta. For binning, we used MaxBin 2.0 (96) and MetaBAT (97), and bins generated by the two different algorithms were merged into one single data set using DAS Tool (98). Bin completeness and contamination were checked using CheckM (99). Further analysis of metagenome-assembled genomes (MAGs) based on results from SqueezeMeta was achieved using Anvi’o v. 6.2 (32) and MAGs selected by DAS tool further refined manually based on hierarchical clustering of contigs. The complete workflow using DAS Tool separates 14 bins whose taxonomic affiliations are uncertain. Manual curation of these results points to “Candidatus Aerophobetes,” Chloroflexi, and “Ca. Atribacteria” as potential bin affiliations (Table S2). In comparison, from the manual curation of MaxBin results, we produced 31 bins (17% to 59%), among which, 12 bins are clearly assigned to “Ca. Atribacteria” with different levels of genome completeness (Fig. S4).

Because the MDA step produced many short fragments that did not allow high-quality binning and full genome completion, assigning a taxonomic affiliation to metagenomic and metatranscriptomic data is challenging (100). For annotating putative functions of ORFs in metagenomes and metatranscriptomes from particular “higher-level” taxonomic groups of microorganisms (34, 35), we also applied a bioinformatics pipeline whereby protein-encoding ORFs were extracted from de novo assembled contigs using FragGeneScan v. 1.30 (101) and functionally annotated against a large aggregated genome database (20, 34, 35) containing predicted proteins from all protist, fungal, bacterial, and archaeal genomes and MAGs in the JGI and NCBI databases using DIAMOND version 0.9.24 (92). This database, which we refer to as MetaProt, also contained all ORFs from all of the transcriptomes of microbial eukaryotes from the MMETS project (102), and we removed any hits to photosynthetic eukaryotic algae as contaminants. This custom MetaProt database that we used for this study is available as a single 32 GB amino acid fasta file on the LMU Open Data website (https://data.ub.uni-muenchen.de/183/). Cutoff values for assigning hits to specific taxa were performed at a minimum bit score of 50, minimum amino acid similarity of 60, and an alignment length of 50 residues. All scripts and code used to produce the analysis have been posted on GitHub (https://github.com/williamorsi/MetaProt-database), and we provide a link to the MetaProt on the GitHub page as well as instructions within the scripts regarding how to conduct the workflows that we used. This approach assigns ORFs to higher-level taxonomic groups (35). As is the case in all metagenomic studies, the incomplete nature of genomes in databases, together with the lower representation of sequenced genomes from candidate clades than from cultured ones, makes it likely that our pipeline misses annotation of ORFs that are derived from as-yet-unsequenced atribacterial genomes. We acknowledge that some genes in databases annotated as being present in “Ca. Atribacteria” might have been assigned to bins according to criteria that differ from study to study (Fig. S7).

Gene identification and normalization in metatranscriptomes.Paired-end reads were trimmed and assembled into contigs using CLC Genomics Workbench 9.5.4 (Qiagen, Hilden, Germany), using a word size of 20, bubble size of 50, and a minimum contig length of 300 nucleotides. Reads were then mapped to the contigs using the following parameters (mismatch penalty, 3; insertion penalty, 3; deletion penalty, 3; minimum alignment length, 50% of read length; minimum percent identity, 95%). Coverage values were obtained from the number of reads mapped to a contig divided by its length (i.e., average coverage). Only contigs with an average coverage of >5 were selected for ORF searches and downstream analysis (20, 34, 85). This protocol does not assemble rRNA; thus, results are only discussed in terms of mRNA. We then performed even further stringency controls by removing any contig that had less than 5× coverage, e.g., reads per kilobase mapped (RPKM). The final resulting data set of contigs was then used for ORF searches and annotation against the MetaProt aggregate database, as described above.

For the metatranscriptomes, normalization of the relative abundance of ORFs was based on the number of unique ORFs assigned to a group (e.g., “Ca. Atribacteria”), a fractional percentage of total ORFs detected. We normalized expression in this manner as opposed to more conventional procedures such as RPKM because we found that the transcriptome sequencing (RNA-seq) kit we used has an amplification step (SPIA amplification, Trio RNAseq Ovation kit; NuGen) that biases the relative abundance of reads mapping to contigs when normalized using RPKM. For example, the RPKM value for the same ORF across technical replicates was found to have very large (orders of magnitude) variability in RPKM. In contrast, the total number of unique ORFs (e.g., presence/absence of an expressed ORF) assigned to specific groups (e.g., “Ca. Atribacteria”) was highly consistent between technical replicates. We assume that this technical variation in the RNA-seq data is associated with randomized SPIA amplification of different transcripts and/or fluctuations in the number of mRNA molecules in technical replicate tubes due to the highly labile nature of RNA during the extraction and library prep procedure. For this reason, we normalized the relative abundance of ORFs assigned to a specific group based on presence/absence of expressed ORFs, which was highly consistent between technical replicates despite the SPIA amplification. If significantly higher numbers of unique ORFs are detected from a particular group than from other groups, it can be attributed to relatively higher transcriptional activity.

COG categories were assigned by searching the ORFs against the COG database (103) using BLASTp. Metagenomic raw reads and metatranscriptomic atribacterial contigs were mapped with high stringency against a previously sequenced subseafloor MAG of “Ca. Atribacteria” that has relatively high completeness (88%) as a reference (104), using Geneious 8.1.9. The average coverage after mapping against GenBank reference no. NCRO00000000.1 to 0.700 (104) was 224 (±131), and 93% of ORFs in the reference MAG were detected in our subseafloor metagenomes and metatranscriptomes. Consensus sequences were exported and annotated using the online tool RAST v. 2.0 (105). Taxonomic assignment of protein-encoding genes to “Ca. Atribacteria” clade JS1 was further confirmed in our metagenomes (Fig. S6) by selecting and aligning 31 phylogenomic markers from the corresponding metagenomic ORFs and 36 atribacterial reference genomes obtained from the NCBI database by using AMPHORA2 (33). Statistical analyses of beta-diversity were performed using RStudio v. 3.3.3 with the Bioconductor package (106).

Data availability.Data are publicly available through NCBI BioProject PRJNA590088. Metagenomes and metatranscriptomes have accession numbers SAMN13317858 to SAMN13317880. The 16S data are available in SRA BioSample accessions SAMN10929403 to SAMN10929517 and SAMN13324854 to SAMN13324920. Additional data related to this paper may be requested from the authors.

ACKNOWLEDGMENTS

This work was supported primarily by the Deutsche Forschungsgemeinschaft (DFG) project OR 417/1-1 granted to W.D.O. Preliminary work was supported by the Center for Dark Energy Biosphere Investigations project OCE-0939564 also granted to W.D.O. The expedition was funded by the US National Science Foundation through grant NSF-OCE-1433150 to S.D. and R.P. R.W.M. led the expedition. Shipboard microbiology efforts were supported by the Center for Dark Energy Biosphere Investigations (C-DEBI grant NSF-OCE-0939564). This is C-DEBI publication 545. This is a contribution of the Deep Carbon Observatory (DCO).

W.D.O. conceived the work and experimental approach. A.V., W.D.O., Ö.K.C., and S.V. contributed to the laboratory/bioinformatics analyses and experimental work. R.W.M., D.C.S., and R.P. obtained the samples during the KN223 R/V Knorr oceanographic expedition. W.D.O., R.W.M., D.C.S., A.V., S.V., and S.D. discussed and wrote the manuscript and commented on the paper.

We declare no competing interests.

FOOTNOTES

    • Received 15 July 2020
    • Accepted 8 September 2020
    • Published 6 October 2020
  • Copyright © 2020 Vuillemin et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Parkes RJ,
    2. John Parkes R,
    3. Cragg BA,
    4. Wellsbury P
    . 2000. Recent studies on bacterial populations and processes in subseafloor sediments: a review. Hydrogeol J 8:11–28. doi:10.1007/PL00010971.
    OpenUrlCrossRefWeb of Science
  2. 2.↵
    1. Inagaki F,
    2. Hinrichs K-U,
    3. Kubo Y,
    4. Bowles MW,
    5. Heuer VB,
    6. Hong W-L,
    7. Hoshino T,
    8. Ijiri A,
    9. Imachi H,
    10. Ito M,
    11. Kaneko M,
    12. Lever MA,
    13. Lin Y-S,
    14. Methe BA,
    15. Morita S,
    16. Morono Y,
    17. Tanikawa W,
    18. Bihan M,
    19. Bowden SA,
    20. Elvert M,
    21. Glombitza C,
    22. Gross D,
    23. Harrington GJ,
    24. Hori T,
    25. Li K,
    26. Limmer D,
    27. Liu C-H,
    28. Murayama M,
    29. Ohkouchi N,
    30. Ono S,
    31. Park Y-S,
    32. Phillips SC,
    33. Prieto-Mollar X,
    34. Purkey M,
    35. Riedinger N,
    36. Sanada Y,
    37. Sauvage J,
    38. Snyder G,
    39. Susilawati R,
    40. Takano Y,
    41. Tasumi E,
    42. Terada T,
    43. Tomaru H,
    44. Trembath-Reichert E,
    45. Wang DT,
    46. Yamada Y
    . 2015. Exploring deep microbial life in coal-bearing sediment down to ∼2.5 km below the ocean floor. Science 349:420–424. doi:10.1126/science.aaa6882.
    OpenUrlAbstract/FREE Full Text
  3. 3.↵
    1. Schippers A,
    2. Neretin LN,
    3. Kallmeyer J,
    4. Ferdelman TG,
    5. Cragg BA,
    6. John Parkes R,
    7. Jørgensen BB
    . 2005. Prokaryotic cells of the deep sub-seafloor biosphere identified as living bacteria. Nature 433:861–864. doi:10.1038/nature03302.
    OpenUrlCrossRefPubMedWeb of Science
  4. 4.↵
    1. Morono Y,
    2. Terada T,
    3. Nishizawa M,
    4. Ito M,
    5. Hillion F,
    6. Takahata N,
    7. Sano Y,
    8. Inagaki F
    . 2011. Carbon and nitrogen assimilation in deep subseafloor microbial cells. Proc Natl Acad Sci U S A 108:18295–18300. doi:10.1073/pnas.1107763108.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Lloyd KG,
    2. May MK,
    3. Kevorkian RT,
    4. Steen AD
    . 2013. Meta-analysis of quantification methods shows that archaea and bacteria have similar abundances in the subseafloor. Appl Environ Microbiol 79:7790–7799. doi:10.1128/AEM.02090-13.
    OpenUrlAbstract/FREE Full Text
  6. 6.↵
    1. D’Hondt S,
    2. Jørgensen BB,
    3. Miller DJ,
    4. Batzke A,
    5. Blake R,
    6. Cragg BA,
    7. Cypionka H,
    8. Dickens GR,
    9. Ferdelman T,
    10. Hinrichs K-U,
    11. Holm NG,
    12. Mitterer R,
    13. Spivack A,
    14. Wang G,
    15. Bekins B,
    16. Engelen B,
    17. Ford K,
    18. Gettemy G,
    19. Rutherford SD,
    20. Sass H,
    21. Skilbeck CG,
    22. Aiello IW,
    23. Guèrin G,
    24. House CH,
    25. Inagaki F,
    26. Meister P,
    27. Naehr T,
    28. Niitsuma S,
    29. Parkes RJ,
    30. Schippers A,
    31. Smith DC,
    32. Teske A,
    33. Wiegel J,
    34. Padilla CN,
    35. Acosta JLS
    . 2004. Distributions of microbial activities in deep subseafloor sediments. Science 306:2216–2221. doi:10.1126/science.1101155.
    OpenUrlAbstract/FREE Full Text
  7. 7.↵
    1. Orsi WD
    . 2018. Ecology and evolution of seafloor and subseafloor microbial communities. Nat Rev Microbiol 16:671–683. doi:10.1038/s41579-018-0046-8.
    OpenUrlCrossRef
  8. 8.↵
    1. Bowles MW,
    2. Mogollon JM,
    3. Kasten S,
    4. Zabel M,
    5. Hinrichs K-U
    . 2014. Global rates of marine sulfate reduction and implications for sub-sea-floor metabolic activities. Science 344:889–891. doi:10.1126/science.1249213.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Jørgensen BB
    . 1982. Mineralization of organic matter in the sea bed—the role of sulphate reduction. Nature 296:643–645. doi:10.1038/296643a0.
    OpenUrlCrossRefWeb of Science
  10. 10.↵
    1. Jørgensen BB,
    2. Findlay AJ,
    3. Pellerin A
    . 2019. The biogeochemical sulfur cycle of marine sediments. Front Microbiol 10:849. doi:10.3389/fmicb.2019.00849.
    OpenUrlCrossRef
  11. 11.↵
    1. D’Hondt S,
    2. Inagaki F,
    3. Zarikian CA,
    4. Abrams LJ,
    5. Dubois N,
    6. Engelhardt T,
    7. Evans H,
    8. Ferdelman T,
    9. Gribsholt B,
    10. Harris RN,
    11. Hoppie BW,
    12. Hyun J-H,
    13. Kallmeyer J,
    14. Kim J,
    15. Lynch JE,
    16. McKinley CC,
    17. Mitsunobu S,
    18. Morono Y,
    19. Murray RW,
    20. Pockalny R,
    21. Sauvage J,
    22. Shimono T,
    23. Shiraishi F,
    24. Smith DC,
    25. Smith-Duque CE,
    26. Spivack AJ,
    27. Steinsbu BO,
    28. Suzuki Y,
    29. Szpak M,
    30. Toffin L,
    31. Uramoto G,
    32. Yamaguchi YT,
    33. Zhang G-L,
    34. Zhang X-H,
    35. Ziebis W
    . 2015. Presence of oxygen and aerobic communities from sea floor to basement in deep-sea sediments. Nature Geosci 8:299–304. doi:10.1038/ngeo2387.
    OpenUrlCrossRef
  12. 12.↵
    1. Kallmeyer J,
    2. Pockalny R,
    3. Adhikari RR,
    4. Smith DC,
    5. D’Hondt S
    . 2012. Global distribution of microbial abundance and biomass in subseafloor sediment. Proc Natl Acad Sci U S A 109:16213–16216. doi:10.1073/pnas.1203849109.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. D’Hondt S,
    2. Spivack AJ,
    3. Pockalny R,
    4. Ferdelman TG,
    5. Fischer JP,
    6. Kallmeyer J,
    7. Abrams LJ,
    8. Smith DC,
    9. Graham D,
    10. Hasiuk F,
    11. Schrum H,
    12. Stancin AM
    . 2009. Subseafloor sedimentary life in the South Pacific gyre. Proc Natl Acad Sci U S A 106:11651–11656. doi:10.1073/pnas.0811793106.
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    1. D’Hondt S,
    2. Pockalny R,
    3. Fulfer VM,
    4. Spivack AJ
    . 2019. Subseafloor life and its biogeochemical impacts. Nat Commun 10:3519. doi:10.1038/s41467-019-11450-z.
    OpenUrlCrossRef
  15. 15.↵
    1. Hoehler TM,
    2. Jørgensen BB
    . 2013. Microbial life under extreme energy limitation. Nat Rev Microbiol 11:83–94. doi:10.1038/nrmicro2939.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Lever MA,
    2. Rogers KL,
    3. Lloyd KG,
    4. Overmann J,
    5. Schink B,
    6. Thauer RK,
    7. Hoehler TM,
    8. Jørgensen BB
    . 2015. Life under extreme energy limitation: a synthesis of laboratory- and field-based investigations. FEMS Microbiol Rev 39:688–728. doi:10.1093/femsre/fuv020.
    OpenUrlCrossRefPubMed
  17. 17.↵
    1. Trembath-Reichert E,
    2. Morono Y,
    3. Ijiri A,
    4. Hoshino T,
    5. Dawson KS,
    6. Inagaki F,
    7. Orphan VJ
    . 2017. Methyl-compound use and slow growth characterize microbial life in 2-km-deep subseafloor coal and shale beds. Proc Natl Acad Sci U S A 114:E9206–E9215. doi:10.1073/pnas.1707525114.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    1. Zhao R,
    2. Hannisdal B,
    3. Mogollon JM,
    4. Jørgensen SL
    . 2019. Nitrifier abundance and diversity peak at deep redox transition zones. Sci Rep 9:8633. doi:10.1038/s41598-019-44585-6.
    OpenUrlCrossRef
  19. 19.↵
    1. Estes ER,
    2. Pockalny R,
    3. D’Hondt S,
    4. Inagaki F,
    5. Morono Y,
    6. Murray RW,
    7. Nordlund D,
    8. Spivack AJ,
    9. Wankel SD,
    10. Xiao N,
    11. Hansel CM
    . 2019. Persistent organic matter in oxic subseafloor sediment. Nature Geosci 12:126–131. doi:10.1038/s41561-018-0291-5.
    OpenUrlCrossRef
  20. 20.↵
    1. Vuillemin A,
    2. Wankel SD,
    3. Coskun ÖK,
    4. Magritsch T,
    5. Vargas S,
    6. Estes ER,
    7. Spivack AJ,
    8. Smith DC,
    9. Pockalny R,
    10. Murray RW,
    11. D’Hondt S,
    12. Orsi WD
    . 2019. Archaea dominate oxic subseafloor communities over multimillion-year time scales. Sci Adv 5:eaaw4108. doi:10.1126/sciadv.aaw4108.
    OpenUrlFREE Full Text
  21. 21.↵
    1. Ziebis W,
    2. McManus J,
    3. Ferdelman T,
    4. Schmidt-Schierhorn F,
    5. Bach W,
    6. Muratli J,
    7. Edwards KJ,
    8. Villinger H
    . 2012. Interstitial fluid chemistry of sediments underlying the North Atlantic gyre and the influence of subsurface fluid flow. Earth Planet Sci Lett 323-324:79–91. doi:10.1016/j.epsl.2012.01.018.
    OpenUrlCrossRef
  22. 22.↵
    1. Parkes RJ,
    2. John Parkes R,
    3. Webster G,
    4. Cragg BA,
    5. Weightman AJ,
    6. Newberry CJ,
    7. Ferdelman TG,
    8. Kallmeyer J,
    9. Jørgensen BB,
    10. Aiello IW,
    11. Fry JC
    . 2005. Deep sub-seafloor prokaryotes stimulated at interfaces over geological time. Nature 436:390–394. doi:10.1038/nature03796.
    OpenUrlCrossRefPubMedWeb of Science
  23. 23.↵
    1. Newberry CJ,
    2. Webster G,
    3. Cragg BA,
    4. John Parkes R,
    5. Weightman AJ,
    6. Fry JC
    . 2004. Diversity of prokaryotes and methanogenesis in deep subsurface sediments from the Nankai Trough, ocean drilling program leg 190. Environ Microbiol 6:274–287. doi:10.1111/j.1462-2920.2004.00568.x.
    OpenUrlCrossRefPubMedWeb of Science
  24. 24.↵
    1. Lloyd KG,
    2. Bird JT,
    3. Buongiorno J,
    4. Deas E,
    5. Kevorkian R,
    6. Noordhoek T,
    7. Rosalsky J,
    8. Roy T
    . 25 March 2020. Evidence for a growth zone for deep subsurface microbial clades in near-surface anoxic sediments. bioRxiv doi:10.1101/2020.03.24.005512.
    OpenUrlCrossRef
  25. 25.↵
    1. Větrovský T,
    2. Baldrian P
    . 2013. The variability of the 16S rRNA gene in bacterial genomes and its consequences for bacterial community analyses. PLoS One 8:e57923. doi:10.1371/journal.pone.0057923.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Schippers A,
    2. Kock D,
    3. Höft C,
    4. Köweker G,
    5. Siegert M
    . 2012. Quantification of microbial communities in subsurface marine sediments of the Black Sea and off Namibia. Front Microbiol 3:e16. doi:10.3389/fmicb.2012.00016.
    OpenUrlCrossRef
  27. 27.↵
    1. Breuker A,
    2. Stadler S,
    3. Schippers A
    . 2013. Microbial community analysis of deeply buried marine sediments of the New Jersey shallow shelf (IODP expedition 313). FEMS Microbiol Ecol 85:578–592. doi:10.1111/1574-6941.12146.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Brown CT,
    2. Hug LA,
    3. Thomas BC,
    4. Sharon I,
    5. Castelle CJ,
    6. Singh A,
    7. Wilkins MJ,
    8. Wrighton KC,
    9. Williams KH,
    10. Banfield JF
    . 2015. Unusual biology across a group comprising more than 15% of domain Bacteria. Nature 523:208–211. doi:10.1038/nature14486.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Castelle CJ,
    2. Wrighton KC,
    3. Thomas BC,
    4. Hug LA,
    5. Brown CT,
    6. Wilkins MJ,
    7. Frischkorn KR,
    8. Tringe SG,
    9. Singh A,
    10. Markillie LM,
    11. Taylor RC,
    12. Williams KH,
    13. Banfield JF
    . 2015. Genomic expansion of domain Archaea highlights roles for organisms from new phyla in anaerobic carbon cycling. Curr Biol 25:690–701. doi:10.1016/j.cub.2015.01.014.
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Spang A,
    2. Saw JH,
    3. Jørgensen SL,
    4. Zaremba-Niedzwiedzka K,
    5. Martijn J,
    6. Lind AE,
    7. van Eijk R,
    8. Schleper C,
    9. Guy L,
    10. Ettema TJG
    . 2015. Complex archaea that bridge the gap between prokaryotes and eukaryotes. Nature 521:173–179. doi:10.1038/nature14447.
    OpenUrlCrossRefPubMed
  31. 31.↵
    1. Anantharaman K,
    2. Brown CT,
    3. Hug LA,
    4. Sharon I,
    5. Castelle CJ,
    6. Probst AJ,
    7. Thomas BC,
    8. Singh A,
    9. Wilkins MJ,
    10. Karaoz U,
    11. Brodie EL,
    12. Williams KH,
    13. Hubbard SS,
    14. Banfield JF
    . 2016. Thousands of microbial genomes shed light on interconnected biogeochemical processes in an aquifer system. Nat Commun 7:13219. doi:10.1038/ncomms13219.
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. Eren AM,
    2. Murat Eren A,
    3. Esen ÖC,
    4. Quince C,
    5. Vineis JH,
    6. Morrison HG,
    7. Sogin ML,
    8. Delmont TO
    . 2015. Anvi’o: an advanced analysis and visualization platform for ‘omics data. PeerJ 3:e1319. doi:10.7717/peerj.1319.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Wu M,
    2. Scott AJ
    . 2012. Phylogenomic analysis of bacterial and archaeal sequences with AMPHORA2. Bioinformatics 28:1033–1034. doi:10.1093/bioinformatics/bts079.
    OpenUrlCrossRefPubMedWeb of Science
  34. 34.↵
    1. Orsi WD,
    2. Richards TA,
    3. Francis WR
    . 2018. Predicted microbial secretomes and their target substrates in marine sediment. Nat Microbiol 3:32–37. doi:10.1038/s41564-017-0047-9.
    OpenUrlCrossRef
  35. 35.↵
    1. Orsi WD,
    2. Vuillemin A,
    3. Rodriguez P,
    4. Coskun ÖK,
    5. Gomez-Saez GV,
    6. Lavik G,
    7. Mohrholz V,
    8. Ferdelman TG
    . 2020. Metabolic activity analyses demonstrate that Lokiarchaeon exhibits homoacetogenesis in sulfidic marine sediments. Nat Microbiol 5:248–255. doi:10.1038/s41564-019-0630-3.
    OpenUrlCrossRef
  36. 36.↵
    1. Orsi WD,
    2. Morard R,
    3. Vuillemin A,
    4. Eitel M,
    5. Wörheide G,
    6. Milucka J,
    7. Kucera M
    . 2020. Anaerobic metabolism of Foraminifera thriving below the seafloor. ISME J 14:2580–2594. doi:10.1038/s41396-020-0708-1.
    OpenUrlCrossRef
  37. 37.↵
    1. Salter SJ,
    2. Cox MJ,
    3. Turek EM,
    4. Calus ST,
    5. Cookson WO,
    6. Moffatt MF,
    7. Turner P,
    8. Parkhill J,
    9. Loman NJ,
    10. Walker AW
    . 2014. Reagent and laboratory contamination can critically impact sequence-based microbiome analyses. BMC Biol 12:87. doi:10.1186/s12915-014-0087-z.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Pichler M,
    2. Coskun ÖK,
    3. Ortega-Arbulú A-S,
    4. Conci N,
    5. Wörheide G,
    6. Vargas S,
    7. Orsi WD
    . 2018. A 16S rRNA gene sequencing and analysis protocol for the Illumina MiniSeq platform. Microbiologyopen 7:e00611. doi:10.1002/mbo3.611.
    OpenUrlCrossRef
  39. 39.↵
    1. Pichoff S,
    2. Lutkenhaus J
    . 2007. Identification of a region of FtsA required for interaction with FtsZ. Mol Microbiol 64:1129–1138. doi:10.1111/j.1365-2958.2007.05735.x.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. Loose M,
    2. Mitchison TJ
    . 2014. The bacterial cell division proteins FtsA and FtsZ self-organize into dynamic cytoskeletal patterns. Nat Cell Biol 16:38–46. doi:10.1038/ncb2885.
    OpenUrlCrossRefPubMedWeb of Science
  41. 41.↵
    1. den Blaauwen T,
    2. den Blaauwen T,
    3. Hamoen LW,
    4. Levin PA
    . 2017. The divisome at 25: the road ahead. Curr Opin Microbiol 36:85–94. doi:10.1016/j.mib.2017.01.007.
    OpenUrlCrossRef
  42. 42.↵
    1. Erickson HP,
    2. Anderson DE,
    3. Osawa M
    . 2010. FtsZ in bacterial cytokinesis: cytoskeleton and force generator all in one. Microbiol Mol Biol Rev 74:504–528. doi:10.1128/MMBR.00021-10.
    OpenUrlAbstract/FREE Full Text
  43. 43.↵
    1. Katayama T,
    2. Nobu MK,
    3. Kusada H,
    4. Meng X-Y,
    5. Yoshioka H,
    6. Kamagata Y,
    7. Tamaki H
    . 8 August 2019. Membrane-bounded nucleoid discovered in a cultivated bacterium of the candidate phylum “Atribacteria.” bioRxiv doi:10.1101/728279.
    OpenUrlAbstract/FREE Full Text
  44. 44.↵
    1. Braun S,
    2. Mhatre SS,
    3. Jaussi M,
    4. Røy H,
    5. Kjeldsen KU,
    6. Pearce C,
    7. Seidenkrantz M-S,
    8. Jørgensen BB,
    9. Lomstein BA
    . 2017. Microbial turnover times in the deep seabed studied by amino acid racemization modelling. Sci Rep 7:e5680. doi:10.1038/s41598-017-05972-z.
    OpenUrlCrossRef
  45. 45.↵
    1. Møller MH,
    2. Glombitza C,
    3. Lever MA,
    4. Deng L,
    5. Morono Y,
    6. Inagaki F,
    7. Doll M,
    8. Su C-C,
    9. Lomstein BA
    . 2018. d:l-Amino acid modeling reveals fast microbial turnover of days to months in the subsurface hydrothermal sediment of Guaymas Basin. Front Microbiol 9:e967. doi:10.3389/fmicb.2018.00967.
    OpenUrlCrossRef
  46. 46.↵
    1. Madigan MT,
    2. Martinko JM,
    3. Stahl DA,
    4. Clark DP
    . 2012. Chapter 5: Microbial growth, p 118–119. In Espinoza D (ed), Brock biology of microorganisms. Pearson Education, Inc., San Francisco, CA.
  47. 47.↵
    1. Ochman H,
    2. Elwyn S,
    3. Moran NA
    . 1999. Calibrating bacterial evolution. Proc Natl Acad Sci U S A 96:12638–12643. doi:10.1073/pnas.96.22.12638.
    OpenUrlAbstract/FREE Full Text
  48. 48.↵
    1. Bird JT,
    2. Tague ED,
    3. Zinke L,
    4. Schmidt JM,
    5. Steen AD,
    6. Reese B,
    7. Marshall IPG,
    8. Webster G,
    9. Weightman A,
    10. Castro HF,
    11. Campagna SR,
    12. Lloyd KG
    . 2019. Uncultured microbial phyla suggest mechanisms for multi-thousand-year subsistence in Baltic Sea sediments. mBio 10:e02376-18. doi:10.1128/mBio.02376-18.
    OpenUrlAbstract/FREE Full Text
  49. 49.↵
    1. Inagaki F,
    2. Suzuki M,
    3. Nealson KH,
    4. Horikoshi K,
    5. D’Hondt SL,
    6. Jørgensen BB
    . 2003. Subseafloor microbial diversity in the Peru Margin (ODP Leg. 201). Geochim Cosmochim Acta 67:171.
    OpenUrlCrossRefWeb of Science
  50. 50.↵
    1. Inagaki F,
    2. Suzuki M,
    3. Takai K,
    4. Oida H,
    5. Sakamoto T,
    6. Aoki K,
    7. Nealson KH,
    8. Horikoshi K
    . 2003. Microbial communities associated with geological horizons in coastal subseafloor sediments from the Sea of Okhotsk. Appl Environ Microbiol 69:7224–7235. doi:10.1128/aem.69.12.7224-7235.2003.
    OpenUrlAbstract/FREE Full Text
  51. 51.↵
    1. Inagaki F,
    2. Nunoura T,
    3. Nakagawa S,
    4. Teske A,
    5. Lever M,
    6. Lauer A,
    7. Suzuki M,
    8. Takai K,
    9. Delwiche M,
    10. Colwell FS,
    11. Nealson KH,
    12. Horikoshi K,
    13. D’Hondt S,
    14. Jørgensen BB
    . 2006. Biogeographical distribution and diversity of microbes in methane hydrate-bearing deep marine sediments on the Pacific Ocean margin. Proc Natl Acad Sci U S A 103:2815–2820. doi:10.1073/pnas.0511033103.
    OpenUrlAbstract/FREE Full Text
  52. 52.↵
    1. Fry JC,
    2. John Parkes R,
    3. Cragg BA,
    4. Weightman AJ,
    5. Webster G
    . 2008. Prokaryotic biodiversity and activity in the deep subseafloor biosphere. FEMS Microbiol Ecol 66:181–196. doi:10.1111/j.1574-6941.2008.00566.x.
    OpenUrlCrossRefPubMedWeb of Science
  53. 53.↵
    1. Starnawski P,
    2. Bataillon T,
    3. Ettema TJG,
    4. Jochum LM,
    5. Schreiber L,
    6. Chen X,
    7. Lever MA,
    8. Polz MF,
    9. Jørgensen BB,
    10. Schramm A,
    11. Kjeldsen KU
    . 2017. Microbial community assembly and evolution in subseafloor sediment. Proc Natl Acad Sci U S A 114:2940–2945. doi:10.1073/pnas.1614190114.
    OpenUrlAbstract/FREE Full Text
  54. 54.↵
    1. Kirkpatrick JB,
    2. Walsh EA,
    3. D’Hondt S
    . 2019. Microbial selection and survival in subseafloor sediment. Front Microbiol 10:956. doi:10.3389/fmicb.2019.00956.
    OpenUrlCrossRef
  55. 55.↵
    1. Dodsworth JA,
    2. Blainey PC,
    3. Murugapiran SK,
    4. Swingley WD,
    5. Ross CA,
    6. Tringe SG,
    7. Chain PSG,
    8. Scholz MB,
    9. Lo C-C,
    10. Raymond J,
    11. Quake SR,
    12. Hedlund BP
    . 2013. Single-cell and metagenomic analyses indicate a fermentative and saccharolytic lifestyle for members of the OP9 lineage. Nat Commun 4:e1854. doi:10.1038/ncomms2884.
    OpenUrlCrossRef
  56. 56.↵
    1. Rinke C,
    2. Schwientek P,
    3. Sczyrba A,
    4. Ivanova NN,
    5. Anderson IJ,
    6. Cheng J-F,
    7. Darling A,
    8. Malfatti S,
    9. Swan BK,
    10. Gies EA,
    11. Dodsworth JA,
    12. Hedlund BP,
    13. Tsiamis G,
    14. Sievert SM,
    15. Liu W-T,
    16. Eisen JA,
    17. Hallam SJ,
    18. Kyrpides NC,
    19. Stepanauskas R,
    20. Rubin EM,
    21. Hugenholtz P,
    22. Woyke T
    . 2013. Insights into the phylogeny and coding potential of microbial dark matter. Nature 499:431–437. doi:10.1038/nature12352.
    OpenUrlCrossRefPubMedWeb of Science
  57. 57.↵
    1. Nobu MK,
    2. Dodsworth JA,
    3. Murugapiran SK,
    4. Rinke C,
    5. Gies EA,
    6. Webster G,
    7. Schwientek P,
    8. Kille P,
    9. John Parkes R,
    10. Sass H,
    11. Jørgensen BB,
    12. Weightman AJ,
    13. Liu W-T,
    14. Hallam SJ,
    15. Tsiamis G,
    16. Woyke T,
    17. Hedlund BP
    . 2016. Phylogeny and physiology of candidate phylum “Atribacteria” (OP9/JS1) inferred from cultivation-independent genomics. ISME J 10:273–286. doi:10.1038/ismej.2015.97.
    OpenUrlCrossRefPubMed
  58. 58.↵
    1. Liu Y-F,
    2. Qi Z-Z,
    3. Shou L-B,
    4. Liu J-F,
    5. Yang S-Z,
    6. Gu J-D,
    7. Mu B-Z
    . 2019. Anaerobic hydrocarbon degradation in candidate phylum “Atribacteria” (JS1) inferred from genomics. ISME J 13:2377–2390. doi:10.1038/s41396-019-0448-2.
    OpenUrlCrossRef
  59. 59.↵
    1. Sewell HL,
    2. Kaster A-K,
    3. Spormann AM
    . 2017. Homoacetogenesis in deep-sea Chloroflexi, as inferred by single-cell genomics, provides a link to reductive dehalogenation in terrestrial Dehalococcoidetes. mBio 8:e02022-17. doi:10.1128/mBio.02022-17.
    OpenUrlAbstract/FREE Full Text
  60. 60.↵
    1. Fincker M,
    2. Huber JA,
    3. Orphan VJ,
    4. Rappé MS,
    5. Teske A,
    6. Spormann AM
    . 2020. Metabolic strategies of marine subseafloor Chloroflexi inferred from genome reconstructions. Environ Microbiol 22:3188–3204. doi:10.1111/1462-2920.15061.
    OpenUrlCrossRef
  61. 61.↵
    1. Schuchmann K,
    2. Müller V
    . 2016. Energetics and application of heterotrophy in acetogenic bacteria. Appl Environ Microbiol 82:4056–4069. doi:10.1128/AEM.00882-16.
    OpenUrlAbstract/FREE Full Text
  62. 62.↵
    1. Hess V,
    2. Schuchmann K,
    3. Müller V
    . 2013. The ferredoxin: NAD+ oxidoreductase (Rnf) from the acetogen Acetobacterium woodii requires Na+ and is reversibly coupled to the membrane potential. J Biol Chem 288:31496–31502. doi:10.1074/jbc.M113.510255.
    OpenUrlAbstract/FREE Full Text
  63. 63.↵
    1. McDowall JS,
    2. Murphy BJ,
    3. Haumann M,
    4. Palmer T,
    5. Armstrong FA,
    6. Sargent F
    . 2014. Bacterial formate hydrogenlyase complex. Proc Natl Acad Sci U S A 111:E3948–E3956. doi:10.1073/pnas.1407927111.
    OpenUrlAbstract/FREE Full Text
  64. 64.↵
    1. Vignais PM,
    2. Billoud B,
    3. Meyer J
    . 2001. Classification and phylogeny of hydrogenases. FEMS Microbiol Rev 25:455–501.
    OpenUrlCrossRefPubMedWeb of Science
  65. 65.↵
    1. Finney AJ,
    2. Sargent F
    . 2019. Formate hydrogenlyase: a group 4 [NiFe]-hydrogenase in tandem with a formate dehydrogenase. Adv Microb Physiol 74:465–486. doi:10.1016/bs.ampbs.2019.02.004.
    OpenUrlCrossRef
  66. 66.↵
    1. Hakobyan M,
    2. Sargsyan H,
    3. Bagramyan K
    . 2005. Proton translocation coupled to formate oxidation in anaerobically grown fermenting Escherichia coli. Biophys Chem 115:55–61. doi:10.1016/j.bpc.2005.01.002.
    OpenUrlCrossRefPubMed
  67. 67.↵
    1. Schauer NL,
    2. Ferry JG
    . 1982. Properties of formate dehydrogenase in Methanobacterium formicicum. J Bacteriol 150:1–7. doi:10.1128/JB.150.1.1-7.1982.
    OpenUrlAbstract/FREE Full Text
  68. 68.↵
    1. Sousa FL,
    2. Martin WF
    . 2014. Biochemical fossils of the ancient transition from geoenergetics to bioenergetics in prokaryotic one carbon compound metabolism. Biochim Biophys Acta 1837:964–981. doi:10.1016/j.bbabio.2014.02.001.
    OpenUrlCrossRef
  69. 69.↵
    1. Orsi WD,
    2. Schink B,
    3. Buckel W,
    4. Martin WF
    . 2020. Physiological limits to life in anoxic subseafloor sediment. FEMS Microbiol Rev 44:219–231. doi:10.1093/femsre/fuaa004.
    OpenUrlCrossRef
  70. 70.↵
    1. Wilbanks B,
    2. Trinh CT
    . 2017. Comprehensive characterization of toxicity of fermentative metabolites on microbial growth. Biotechnol Biofuels 10:262. doi:10.1186/s13068-017-0952-4.
    OpenUrlCrossRef
  71. 71.↵
    1. LaRowe DE,
    2. Van Cappellen P
    . 2011. Degradation of natural organic matter: a thermodynamic analysis. Geochim Cosmochim Acta 75:2030–2042. doi:10.1016/j.gca.2011.01.020.
    OpenUrlCrossRefWeb of Science
  72. 72.↵
    1. Axen SD,
    2. Erbilgin O,
    3. Kerfeld CA
    . 2014. A taxonomy of bacterial microcompartment loci constructed by a novel scoring method. PLoS Comput Biol 10:e1003898. doi:10.1371/journal.pcbi.1003898.
    OpenUrlCrossRefPubMed
  73. 73.↵
    1. Kerfeld CA,
    2. Erbilgin O
    . 2015. Bacterial microcompartments and the modular construction of microbial metabolism. Trends Microbiol 23:22–34. doi:10.1016/j.tim.2014.10.003.
    OpenUrlCrossRefPubMed
  74. 74.↵
    1. Vuillemin A,
    2. Kerrigan Z,
    3. D’Hondt S,
    4. Orsi WD
    . 27 May 2020. Chloroflexi persisting for millions of years in oxic and anoxic deep-sea clay. bioRxiv doi:10.1101/2020.05.26.116590.
    OpenUrlAbstract/FREE Full Text
  75. 75.↵
    1. Zhuang W-Q,
    2. Yi S,
    3. Bill M,
    4. Brisson VL,
    5. Feng X,
    6. Men Y,
    7. Conrad ME,
    8. Tang YJ,
    9. Alvarez-Cohen L
    . 2014. Incomplete Wood-Ljungdahl pathway facilitates one-carbon metabolism in organohalide-respiring Dehalococcoides mccartyi. Proc Natl Acad Sci U S A 111:6419–6424. doi:10.1073/pnas.1321542111.
    OpenUrlAbstract/FREE Full Text
  76. 76.↵
    1. Fischer JP,
    2. Ferdelman TG,
    3. D’Hondt S,
    4. Røy H,
    5. Wenzhöfer F
    . 2009. Oxygen penetration deep into the sediment of the South Pacific gyre. Biogeosciences 6:1467–1478. doi:10.5194/bg-6-1467-2009.
    OpenUrlCrossRef
  77. 77.↵
    1. Parada AE,
    2. Needham DM,
    3. Fuhrman JA
    . 2016. Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ Microbiol 18:1403–1414. doi:10.1111/1462-2920.13023.
    OpenUrlCrossRefPubMed
  78. 78.↵
    1. Jochum LM,
    2. Chen X,
    3. Lever MA,
    4. Loy A,
    5. Jørgensen BB,
    6. Schramm A,
    7. Kjeldsen KU
    . 2017. Depth distribution and assembly of sulfate-reducing microbial communities in marine sediments of Aarhus Bay. Appl Environ Microbiol 83:e01547-17. doi:10.1128/AEM.01547-17.
    OpenUrlAbstract/FREE Full Text
  79. 79.↵
    1. Coskun ÖK,
    2. Özen V,
    3. Wankel SD,
    4. Orsi WD
    . 2019. Quantifying population-specific growth in benthic bacterial communities under low oxygen using H218O. ISME J 13:1546–1559. doi:10.1038/s41396-019-0373-4.
    OpenUrlCrossRef
  80. 80.↵
    1. Edgar RC
    . 2013. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods 10:996–998. doi:10.1038/nmeth.2604.
    OpenUrlCrossRefPubMedWeb of Science
  81. 81.↵
    1. Quast C,
    2. Pruesse E,
    3. Yilmaz P,
    4. Gerken J,
    5. Schweer T,
    6. Yarza P,
    7. Peplies J,
    8. Glöckner FO
    . 2012. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res 41:D590–D596. doi:10.1093/nar/gks1219.
    OpenUrlCrossRefPubMedWeb of Science
  82. 82.↵
    1. Pruesse E,
    2. Quast C,
    3. Knittel K,
    4. Fuchs BM,
    5. Ludwig W,
    6. Peplies J,
    7. Glockner FO
    . 2007. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res 35:7188–7196. doi:10.1093/nar/gkm864.
    OpenUrlCrossRefPubMedWeb of Science
  83. 83.↵
    1. Stamatakis A
    . 2006. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22:2688–2690. doi:10.1093/bioinformatics/btl446.
    OpenUrlCrossRefPubMedWeb of Science
  84. 84.↵
    1. Ludwig W,
    2. Strunk O,
    3. Westram R,
    4. Richter L,
    5. Meier H,
    6. Yadhukumar Buchner A,
    7. Lai T,
    8. Steppi S,
    9. Jobb G,
    10. Förster W,
    11. Brettske I,
    12. Gerber S,
    13. Ginhart AW,
    14. Gross O,
    15. Grumann S,
    16. Hermann S,
    17. Jost R,
    18. König A,
    19. Liss T,
    20. Lüssmann R,
    21. May M,
    22. Nonhoff B,
    23. Reichel B,
    24. Strehlow R,
    25. Stamatakis A,
    26. Stuckmann N,
    27. Vilbig A,
    28. Lenke M,
    29. Ludwig T,
    30. Bode A,
    31. Schleifer K-H
    . 2004. ARB: a software environment for sequence data. Nucleic Acids Res 32:1363–1371. doi:10.1093/nar/gkh293.
    OpenUrlCrossRefPubMedWeb of Science
  85. 85.↵
    1. Ortega-Arbulú A,
    2. Pichler M,
    3. Vuillemin A,
    4. Orsi WD
    . 2019. Effects of organic matter and low oxygen on the mycobenthos in a coastal lagoon. Environ Microbiol 21:374–388. doi:10.1111/1462-2920.14469.
    OpenUrlCrossRef
  86. 86.↵
    1. Tamames J,
    2. Puente-Sánchez F
    . 2019. SqueezeMeta, a highly portable, fully automatic metagenomic analysis pipeline. Front Microbiol 9:e3349. doi:10.3389/fmicb.2018.03349.
    OpenUrlCrossRef
  87. 87.↵
    1. Bolger AM,
    2. Lohse M,
    3. Usadel B
    . 2014. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30:2114–2120. doi:10.1093/bioinformatics/btu170.
    OpenUrlCrossRefPubMedWeb of Science
  88. 88.↵
    1. Li D,
    2. Liu C-M,
    3. Luo R,
    4. Sadakane K,
    5. Lam T-W
    . 2015. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 31:1674–1676. doi:10.1093/bioinformatics/btv033.
    OpenUrlCrossRefPubMed
  89. 89.↵
    1. Hyatt D,
    2. Chen G-L,
    3. LoCascio PF,
    4. Land ML,
    5. Larimer FW,
    6. Hauser LJ
    . 2010. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics 11:119. doi:10.1186/1471-2105-11-119.
    OpenUrlCrossRefPubMed
  90. 90.↵
    1. Seemann T
    . 2014. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30:2068–2069. doi:10.1093/bioinformatics/btu153.
    OpenUrlCrossRefPubMedWeb of Science
  91. 91.↵
    1. Wang Q,
    2. Garrity GM,
    3. Tiedje JM,
    4. Cole JR
    . 2007. Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol 73:5261–5267. doi:10.1128/AEM.00062-07.
    OpenUrlAbstract/FREE Full Text
  92. 92.↵
    1. Buchfink B,
    2. Xie C,
    3. Huson DH
    . 2015. Fast and sensitive protein alignment using DIAMOND. Nat Methods 12:59–60. doi:10.1038/nmeth.3176.
    OpenUrlCrossRefPubMed
  93. 93.↵
    1. Huerta-Cepas J,
    2. Szklarczyk D,
    3. Heller D,
    4. Hernández-Plaza A,
    5. Forslund SK,
    6. Cook H,
    7. Mende DR,
    8. Letunic I,
    9. Rattei T,
    10. Jensen LJ,
    11. von Mering C,
    12. Bork P
    . 2019. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res 47:D309–D314. doi:10.1093/nar/gky1085.
    OpenUrlCrossRefPubMed
  94. 94.↵
    1. Kanehisa M,
    2. Goto S
    . 2000. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res 28:27–30. doi:10.1093/nar/28.1.27.
    OpenUrlCrossRefPubMedWeb of Science
  95. 95.↵
    1. Langmead B,
    2. Salzberg SL
    . 2012. Fast gapped-read alignment with Bowtie 2. Nat Methods 9:357–359. doi:10.1038/nmeth.1923.
    OpenUrlCrossRefPubMedWeb of Science
  96. 96.↵
    1. Wu Y-W,
    2. Simmons BA,
    3. Singer SW
    . 2016. MaxBin 2.0: an automated binning algorithm to recover genomes from multiple metagenomic datasets. Bioinformatics 32:605–607. doi:10.1093/bioinformatics/btv638.
    OpenUrlCrossRefPubMed
  97. 97.↵
    1. Kang DD,
    2. Froula J,
    3. Egan R,
    4. Wang Z
    . 2015. MetaBAT, an efficient tool for accurately reconstructing single genomes from complex microbial communities. PeerJ 3:e1165. doi:10.7717/peerj.1165.
    OpenUrlCrossRefPubMed
  98. 98.↵
    1. Sieber CMK,
    2. Probst AJ,
    3. Sharrar A,
    4. Thomas BC,
    5. Hess M,
    6. Tringe SG,
    7. Banfield JF
    . 2018. Recovery of genomes from metagenomes via a dereplication, aggregation and scoring strategy. Nat Microbiol 3:836–843. doi:10.1038/s41564-018-0171-1.
    OpenUrlCrossRef
  99. 99.↵
    1. Parks DH,
    2. Imelfort M,
    3. Skennerton CT,
    4. Hugenholtz P,
    5. Tyson GW
    . 2015. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res 25:1043–1055. doi:10.1101/gr.186072.114.
    OpenUrlAbstract/FREE Full Text
  100. 100.↵
    1. Breitwieser FP,
    2. Lu J,
    3. Salzberg SL
    . 2019. A review of methods and databases for metagenomic classification and assembly. Brief Bioinform 20:1125–1136. doi:10.1093/bib/bbx120.
    OpenUrlCrossRef
  101. 101.↵
    1. Rho M,
    2. Tang H,
    3. Ye Y
    . 2010. FragGeneScan: predicting genes in short and error-prone reads. Nucleic Acids Res 38:e191. doi:10.1093/nar/gkq747.
    OpenUrlCrossRefPubMed
  102. 102.↵
    1. Keeling PJ,
    2. Burki F,
    3. Wilcox HM,
    4. Allam B,
    5. Allen EE,
    6. Amaral-Zettler LA,
    7. Armbrust EV,
    8. Archibald JM,
    9. Bharti AK,
    10. Bell CJ,
    11. Beszteri B,
    12. Bidle KD,
    13. Cameron CT,
    14. Campbell L,
    15. Caron DA,
    16. Cattolico RA,
    17. Collier JL,
    18. Coyne K,
    19. Davy SK,
    20. Deschamps P,
    21. Dyhrman ST,
    22. Edvardsen B,
    23. Gates RD,
    24. Gobler CJ,
    25. Greenwood SJ,
    26. Guida SM,
    27. Jacobi JL,
    28. Jakobsen KS,
    29. James ER,
    30. Jenkins B,
    31. John U,
    32. Johnson MD,
    33. Juhl AR,
    34. Kamp A,
    35. Katz LA,
    36. Kiene R,
    37. Kudryavtsev A,
    38. Leander BS,
    39. Lin S,
    40. Lovejoy C,
    41. Lynn D,
    42. Marchetti A,
    43. McManus G,
    44. Nedelcu AM,
    45. Menden-Deuer S,
    46. Miceli C,
    47. Mock T,
    48. Montresor M,
    49. Moran MA,
    50. Murray S, et al
    . 2014. The Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP): illuminating the functional diversity of eukaryotic life in the oceans through transcriptome sequencing. PLoS Biol 12:e1001889. doi:10.1371/journal.pbio.1001889.
    OpenUrlCrossRefPubMed
  103. 103.↵
    1. Galperin MY,
    2. Makarova KS,
    3. Wolf YI,
    4. Koonin EV
    .2015. Expanded microbial genome coverage and improved protein family annotation in the COG database. Nucleic Acids Res 43:D261–D269. doi:10.1093/nar/gku1223.
    OpenUrlCrossRefPubMed
  104. 104.↵
    1. Lee YM,
    2. Hwang K,
    3. Lee JI,
    4. Kim M,
    5. Hwang CY,
    6. Noh H-J,
    7. Choi H,
    8. Lee HK,
    9. Chun J,
    10. Hong SG,
    11. Shin SC
    . 2018. Genomic insight into the predominance of candidate phylum Atribacteria JS1 lineage in marine sediments. Front Microbiol 9:2909. doi:10.3389/fmicb.2018.02909.
    OpenUrlCrossRef
  105. 105.↵
    1. Aziz RK,
    2. Bartels D,
    3. Best AA,
    4. DeJongh M,
    5. Disz T,
    6. Edwards RA,
    7. Formsma K,
    8. Gerdes S,
    9. Glass EM,
    10. Kubal M,
    11. Meyer F,
    12. Olsen GJ,
    13. Olson R,
    14. Osterman AL,
    15. Overbeek RA,
    16. McNeil LK,
    17. Paarmann D,
    18. Paczian T,
    19. Parrello B,
    20. Pusch GD,
    21. Reich C,
    22. Stevens R,
    23. Vassieva O,
    24. Vonstein V,
    25. Wilke A,
    26. Zagnitko O
    . 2008. The RAST Server: rapid annotations using subsystems technology. BMC Genomics 9:75. doi:10.1186/1471-2164-9-75.
    OpenUrlCrossRefPubMed
  106. 106.↵
    1. Huber W,
    2. Carey VJ,
    3. Gentleman R,
    4. Anders S,
    5. Carlson M,
    6. Carvalho BS,
    7. Bravo HC,
    8. Davis S,
    9. Gatto L,
    10. Girke T,
    11. Gottardo R,
    12. Hahne F,
    13. Hansen KD,
    14. Irizarry RA,
    15. Lawrence M,
    16. Love MI,
    17. MacDonald J,
    18. Obenchain V,
    19. Oleś AK,
    20. Pagès H,
    21. Reyes A,
    22. Shannon P,
    23. Smyth GK,
    24. Tenenbaum D,
    25. Waldron L,
    26. Morgan M
    . 2015. Orchestrating high-throughput genomic analysis with Bioconductor. Nat Methods 12:115–121. doi:10.1038/nmeth.3252.
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Download PDF
Citation Tools
Atribacteria Reproducing over Millions of Years in the Atlantic Abyssal Subseafloor
Aurèle Vuillemin, Sergio Vargas, Ömer K. Coskun, Robert Pockalny, Richard W. Murray, David C. Smith, Steven D’Hondt, William D. Orsi
mBio Oct 2020, 11 (5) e01937-20; DOI: 10.1128/mBio.01937-20

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mBio article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Atribacteria Reproducing over Millions of Years in the Atlantic Abyssal Subseafloor
(Your Name) has forwarded a page to you from mBio
(Your Name) thought you would be interested in this article in mBio.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Atribacteria Reproducing over Millions of Years in the Atlantic Abyssal Subseafloor
Aurèle Vuillemin, Sergio Vargas, Ömer K. Coskun, Robert Pockalny, Richard W. Murray, David C. Smith, Steven D’Hondt, William D. Orsi
mBio Oct 2020, 11 (5) e01937-20; DOI: 10.1128/mBio.01937-20
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS AND DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

deep biosphere
energy limit to life
atribacteria
acetogenesis
metagenomics
transcriptomics
fermentation
bacterial microcompartment
clade JS1
metatranscriptomics
subseafloor life

Related Articles

Cited By...

About

  • About mBio
  • Editor in Chief
  • Board of Editors
  • AAM Fellows
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Article Types
  • Ethics
  • Contact Us

Follow #mBio

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2150-7511