Cultivation of novel Atribacterota from oil well provides new insight into their diversity, ecology, and evolution in anoxic, carbon-rich environments.

BACKGROUND: The Atribacterota are widely distributed in the subsurface biosphere. Recently, the first Atribacterota isolate was described and the number of Atribacterota genome sequences retrieved from environmental samples has increased significantly; however, their diversity, physiology, ecology, and evolution remain poorly understood. RESULTS: We report the isolation of the second member of Atribacterota, Thermatribacter velox gen. nov., sp. nov., within a new family Thermatribacteraceae fam. nov., and the short-term laboratory cultivation of a member of the JS1 lineage, Phoenicimicrobium oleiphilum HX-OS.bin.34TS, both from a terrestrial oil reservoir. Physiological and metatranscriptomics analyses showed that Thermatribacter velox B11T and Phoenicimicrobium oleiphilum HX-OS.bin.34TS ferment sugars and n-alkanes, respectively, producing H2, CO2, and acetate as common products. Comparative genomics showed that all members of the Atribacterota lack a complete Wood-Ljungdahl Pathway (WLP), but that the Reductive Glycine Pathway (RGP) is widespread, indicating that the RGP, rather than WLP, is a central hub in Atribacterota metabolism. Ancestral character state reconstructions and phylogenetic analyses showed that key genes encoding the RGP (fdhA, fhs, folD, glyA, gcvT, gcvPAB, pdhD) and other central functions were gained independently in the two classes, Atribacteria (OP9) and Phoenicimicrobiia (JS1), after which they were inherited vertically; these genes included fumarate-adding enzymes (faeA; Phoenicimicrobiia only), the CODH/ACS complex (acsABCDE), and diverse hydrogenases (NiFe group 3b, 4b and FeFe group A3, C). Finally, we present genome-resolved community metabolic models showing the central roles of Atribacteria (OP9) and Phoenicimicrobiia (JS1) in acetate- and hydrocarbon-rich environments. CONCLUSION: Our findings expand the knowledge of the diversity, physiology, ecology, and evolution of the phylum Atribacterota. This study is a starting point for promoting more incisive studies of their syntrophic biology and may guide the rational design of strategies to cultivate them in the laboratory. Video Abstract.

Resource partitioning and amino acid assimilation in a terrestrial geothermal spring.

High-temperature geothermal springs host simplified microbial communities; however, the activities of individual microorganisms and their roles in the carbon cycle in nature are not well understood. Here, quantitative stable isotope probing (qSIP) was used to track the assimilation of 13C-acetate and 13C-aspartate into DNA in 74 °C sediments in Gongxiaoshe Hot Spring, Tengchong, China. This revealed a community-wide preference for aspartate and a tight coupling between aspartate incorporation into DNA and the proliferation of aspartate utilizers during labeling. Both 13C incorporation into DNA and changes in the abundance of taxa during incubations indicated strong resource partitioning and a significant phylogenetic signal for aspartate incorporation. Of the active amplicon sequence variants (ASVs) identified by qSIP, most could be matched with genomes from Gongxiaoshe Hot Spring or nearby springs with an average nucleotide similarity of 99.4%. Genomes corresponding to aspartate primary utilizers were smaller, near-universally encoded polar amino acid ABC transporters, and had codon preferences indicative of faster growth rates. The most active ASVs assimilating both substrates were not abundant, suggesting an important role for the rare biosphere in the community response to organic carbon addition. The broad incorporation of aspartate into DNA over acetate by the hot spring community may reflect dynamic cycling of cell lysis products in situ or substrates delivered during monsoon rains and may reflect N limitation.

Mcr-dependent methanogenesis in Archaeoglobaceae enriched from a terrestrial hot spring.

The preeminent source of biological methane on Earth is methyl coenzyme M reductase (Mcr)-dependent archaeal methanogenesis. A growing body of evidence suggests a diversity of archaea possess Mcr, although experimental validation of hypothesized methane metabolisms has been missing. Here, we provide evidence of a functional Mcr-based methanogenesis pathway in a novel member of the family Archaeoglobaceae, designated Methanoglobus nevadensis, which we enriched from a terrestrial hot spring on the polysaccharide xyloglucan. Our incubation assays demonstrate methane production that is highly sensitive to the Mcr inhibitor bromoethanesulfonate, stimulated by xyloglucan and xyloglucan-derived sugars, concomitant with the consumption of molecular hydrogen, and causing a deuterium fractionation in methane characteristic of hydrogenotrophic and methylotrophic methanogens. Combined with the recovery and analysis of a high-quality M. nevadensis metagenome-assembled genome encoding a divergent Mcr and diverse potential electron and carbon transfer pathways, our observations suggest methanogenesis in M. nevadensis occurs via Mcr and is fueled by the consumption of cross-fed byproducts of xyloglucan fermentation mediated by other community members. Phylogenetic analysis shows close affiliation of the M. nevadensis Mcr with those from Korarchaeota, Nezhaarchaeota, Verstraetearchaeota, and other Archaeoglobales that are divergent from well-characterized Mcr. We propose these archaea likely also use functional Mcr complexes to generate methane on the basis of our experimental validation in M. nevadensis. Thus, divergent Mcr-encoding archaea may be underestimated sources of biological methane in terrestrial and marine hydrothermal environments.

Proposed minimal standards for description of methanogenic archaea.

Methanogenic archaea are a diverse, polyphyletic group of strictly anaerobic prokaryotes capable of producing methane as their primary metabolic product. It has been over three decades since minimal standards for their taxonomic description have been proposed. In light of advancements in technology and amendments in systematic microbiology, revision of the older criteria for taxonomic description is essential. Most of the previously recommended minimum standards regarding phenotypic characterization of pure cultures are maintained. Electron microscopy and chemotaxonomic methods like whole-cell protein and lipid analysis are desirable but not required. Because of advancements in DNA sequencing technologies, obtaining a complete or draft whole genome sequence for type strains and its deposition in a public database are now mandatory. Genomic data should be used for rigorous comparison to close relatives using overall genome related indices such as average nucleotide identity and digital DNA-DNA hybridization. Phylogenetic analysis of the 16S rRNA gene is also required and can be supplemented by phylogenies of the mcrA gene and phylogenomic analysis using multiple conserved, single-copy marker genes. Additionally, it is now established that culture purity is not essential for studying prokaryotes, and description of Candidatus methanogenic taxa using single-cell or metagenomics along with other appropriate criteria is a viable alternative. The revisions to the minimal criteria proposed here by the members of the Subcommittee on the Taxonomy of Methanogenic Archaea of the International Committee on Systematics of Prokaryotes should allow for rigorous yet practical taxonomic description of these important and diverse microbes.

An essential role for tungsten in the ecology and evolution of a previously uncultivated lineage of anaerobic, thermophilic Archaea.

Trace metals have been an important ingredient for life throughout Earth's history. Here, we describe the genome-guided cultivation of a member of the elusive archaeal lineage Caldarchaeales (syn. Aigarchaeota), Wolframiiraptor gerlachensis, and its growth dependence on tungsten. A metagenome-assembled genome (MAG) of W. gerlachensis encodes putative tungsten membrane transport systems, as well as pathways for anaerobic oxidation of sugars probably mediated by tungsten-dependent ferredoxin oxidoreductases that are expressed during growth. Catalyzed reporter deposition-fluorescence in-situ hybridization (CARD-FISH) and nanoscale secondary ion mass spectrometry (nanoSIMS) show that W. gerlachensis preferentially assimilates xylose. Phylogenetic analyses of 78 high-quality Wolframiiraptoraceae MAGs from terrestrial and marine hydrothermal systems suggest that tungsten-associated enzymes were present in the last common ancestor of extant Wolframiiraptoraceae. Our observations imply a crucial role for tungsten-dependent metabolism in the origin and evolution of this lineage, and hint at a relic metabolic dependence on this trace metal in early anaerobic thermophiles.

Genomics, Exometabolomics, and Metabolic Probing Reveal Conserved Proteolytic Metabolism of Thermoflexus hugenholtzii and Three Candidate Species From China and Japan.

Thermoflexus hugenholtzii JAD2T, the only cultured representative of the Chloroflexota order Thermoflexales, is abundant in Great Boiling Spring (GBS), NV, United States, and close relatives inhabit geothermal systems globally. However, no defined medium exists for T. hugenholtzii JAD2T and no single carbon source is known to support its growth, leaving key knowledge gaps in its metabolism and nutritional needs. Here, we report comparative genomic analysis of the draft genome of T. hugenholtzii JAD2T and eight closely related metagenome-assembled genomes (MAGs) from geothermal sites in China, Japan, and the United States, representing "Candidatus Thermoflexus japonica," "Candidatus Thermoflexus tengchongensis," and "Candidatus Thermoflexus sinensis." Genomics was integrated with targeted exometabolomics and 13C metabolic probing of T. hugenholtzii. The Thermoflexus genomes each code for complete central carbon metabolic pathways and an unusually high abundance and diversity of peptidases, particularly Metallo- and Serine peptidase families, along with ABC transporters for peptides and some amino acids. The T. hugenholtzii JAD2T exometabolome provided evidence of extracellular proteolytic activity based on the accumulation of free amino acids. However, several neutral and polar amino acids appear not to be utilized, based on their accumulation in the medium and the lack of annotated transporters. Adenine and adenosine were scavenged, and thymine and nicotinic acid were released, suggesting interdependency with other organisms in situ. Metabolic probing of T. hugenholtzii JAD2T using 13C-labeled compounds provided evidence of oxidation of glucose, pyruvate, cysteine, and citrate, and functioning glycolytic, tricarboxylic acid (TCA), and oxidative pentose-phosphate pathways (PPPs). However, differential use of position-specific 13C-labeled compounds showed that glycolysis and the TCA cycle were uncoupled. Thus, despite the high abundance of Thermoflexus in sediments of some geothermal systems, they appear to be highly focused on chemoorganotrophy, particularly protein degradation, and may interact extensively with other microorganisms in situ.

Diversity and Distribution of a Novel Genus of Hyperthermophilic Aquificae Viruses Encoding a Proof-Reading Family-A DNA Polymerase.

Despite the high abundance of Aquificae in many geothermal systems, these bacteria are difficult to culture and no viruses infecting members of this phylum have been isolated. Here, we describe the complete, circular dsDNA Uncultivated Virus Genome (UViG) of Thermocrinis Octopus Spring virus (TOSV), derived from metagenomic data, along with eight related UViGs representing three additional viral species. Despite low overall similarity among viruses from different hot springs, the genomes shared a high degree of synteny, and encoded numerous genes for nucleotide metabolism, including a PolA-type DNA polymerase polyprotein with likely accessory functions, a DNA Pol III sliding clamp, a thymidylate kinase, a DNA gyrase, a helicase, and a DNA methylase. Also present were conserved genes predicted to code for phage capsid, large and small subunits of terminase, portal protein, holin, and lytic transglycosylase, all consistent with a distant relatedness to cultivated Caudovirales. These viruses are predicted to infect Aquificae, as multiple CRISPR spacers matching the viral genomes were identified within the genomes and metagenomic contigs from these bacteria. Based on the predicted atypical bi-directional replication strategy, low sequence similarity to known viral genomes, and unique position in gene-sharing networks, we propose a new putative genus, "Pyrovirus," in the order Caudovirales.

Diverse respiratory capacity among Thermus strains from US Great Basin hot springs.

Thermus species are thermophilic heterotrophs, with most capable of using a variety of organic and inorganic electron donors for respiration. Here, a combined cultivation-independent and -dependent approach was used to explore the diversity of Thermus in Great Boiling Spring (GBS) and Little Hot Creek (LHC) in the US Great Basin. A cultivation-independent 16S rRNA gene survey of ten LHC sites showed that Thermus made up 0-3.5% of sequences and were predominately Thermus thermophilus. 189 Thermus isolates from GBS and LHC were affiliated with T. aquaticus (73.0%), T. oshimai (25.4%), T. sediminis (1.1%), and T. thermophilus (0.5%), with T. aquaticus and T. oshimai forming biogeographic clusters. 22 strains were selected for characterization, including chemolithotrophic oxidation of thiosulfate and arsenite, and reduction of ferric iron, polysulfide, and nitrate, revealing phenotypic diversity and broad respiratory capability within each species. PCR demonstrated the wide distribution of aerobic arsenite oxidase genes. A GBS sediment metaproteome contained sulfite oxidase and Fe3+ ABC transporter permease peptides, suggesting sulfur and iron transformations in situ. This study expands our knowledge of the physiological diversity of Thermus, suggesting widespread chemolithotrophic and anaerobic respiration phenotypes, and providing a foundation for better understanding the ecology of this genus in thermal ecosystems.

Position-Specific Metabolic Probing and Metagenomics of Microbial Communities Reveal Conserved Central Carbon Metabolic Network Activities at High Temperatures.

Temperature is a primary driver of microbial community composition and taxonomic diversity; however, it is unclear to what extent temperature affects characteristics of central carbon metabolic pathways (CCMPs) at the community level. In this study, 16S rRNA gene amplicon and metagenome sequencing were combined with 13C-labeled metabolite probing of the CCMPs to assess community carbon metabolism along a temperature gradient (60-95°C) in Great Boiling Spring, NV. 16S rRNA gene amplicon diversity was inversely proportional to temperature, and Archaea were dominant at higher temperatures. KO richness and diversity were also inversely proportional to temperature, yet CCMP genes were similarly represented across the temperature gradient and many individual metagenome-assembled genomes had complete pathways. In contrast, genes encoding cellulosomes and many genes involved in plant matter degradation and photosynthesis were absent at higher temperatures. In situ 13C-CO2 production from labeled isotopomer pairs of glucose, pyruvate, and acetate suggested lower relative oxidative pentose phosphate pathway activity and/or fermentation at 60°C, and a stable or decreased maintenance energy demand at higher temperatures. Catabolism of 13C-labeled citrate, succinate, L-alanine, L-serine, and L-cysteine was observed at 85°C, demonstrating broad heterotrophic activity and confirming functioning of the TCA cycle. Together, these results suggest that temperature-driven losses in biodiversity and gene content in geothermal systems may not alter CCMP function or maintenance energy demands at a community level.

Thermus sediminis sp. nov., a thiosulfate-oxidizing and arsenate-reducing organism isolated from Little Hot Creek in the Long Valley Caldera, California.

Thermus species are widespread in natural and artificial thermal environments. Two new yellow-pigmented strains, L198T and L423, isolated from Little Hot Creek, a geothermal spring in eastern California, were identified as novel organisms belonging to the genus Thermus. Cells are Gram-negative, rod-shaped, and non-motile. Growth was observed at temperatures from 45 to 75 °C and at salinities of 0-2.0% added NaCl. Both strains grow heterotrophically or chemolithotrophically by oxidation of thiosulfate to sulfate. L198T and L423 grow by aerobic respiration or anaerobic respiration with arsenate as the terminal electron acceptor. Values for 16S rRNA gene identity (≤ 97.01%), digital DNA-DNA hybridization (≤ 32.7%), OrthoANI (≤ 87.5%), and genome-to-genome distance (0.13) values to all Thermus genomes were less than established criteria for microbial species. The predominant respiratory quinone was menaquinone-8 and the major cellular fatty acids were iso-C15:0, iso-C17:0 and anteiso-C15:0. One unidentified phospholipid (PL1) and one unidentified glycolipid (GL1) dominated the polar lipid pattern. The new strains could be differentiated from related taxa by ß-galactosidase and ß-glucosidase activity and the presence of hydroxy fatty acids. Based on phylogenetic, genomic, phenotypic, and chemotaxonomic evidence, the novel species Thermus sediminis sp. nov. is proposed, with the type strain L198T (= CGMCC 1.13590T = KCTC XXX).

Localized high abundance of Marine Group II archaea in the subtropical Pearl River Estuary: implications for their niche adaptation.

Marine Group II archaea are widely distributed in global oceans and dominate the total archaeal community within the upper euphotic zone of temperate waters. However, factors controlling the distribution of MGII are poorly delineated and the physiology and ecological functions of these still-uncultured organisms remain elusive. In this study, we investigated the planktonic MGII associated with particles and in free-living forms in the Pearl River Estuary (PRE) over a 10-month period. We detected high abundance of particle-associated MGII in PRE (up to ∼108 16S rRNA gene copies/l), which was around 10-fold higher than the free-living MGII in the same region, and an order of magnitude higher than previously reported in other marine environments. 10‰ salinity appeared to be a threshold value for these MGII because MGII abundance decreased sharply below it. Above 10‰ salinity, the abundance of MGII on the particles was positively correlated with phototrophs and MGII in the surface water was negatively correlated with irradiance. However, the abundances of those free-living MGII showed positive correlations with salinity and temperature, suggesting the different physiological characteristics between particle-attached and free-living MGIIs. A nearly completely assembled metagenome, MGIIa_P, was recovered using metagenome binning methods. Compared with the other two MGII genomes from surface ocean, MGIIa_P contained higher proportions of glycoside hydrolases, indicating the ability of MGIIa_P to hydrolyse glycosidic bonds in complex sugars in PRE. MGIIa_P is the first assembled MGII metagenome containing a catalase gene, which might be involved in scavenging reactive oxygen species generated by the abundant phototrophs in the eutrophic PRE. Our study presented the widespread and high abundance of MGII in the water columns of PRE, and characterized the determinant abiotic factors affecting their distribution. Their association with heterotrophs, preference for particles and resourceful metabolic traits indicate MGII might play a significant role in metabolising organic matters in the PRE and other temperate estuarine systems.

Genomic Comparison of Two Family-Level Groups of the Uncultivated NAG1 Archaeal Lineage from Chemically and Geographically Disparate Hot Springs.

Recent progress based on single-cell genomics and metagenomic investigations of archaea in a variety of extreme environments has led to significant advances in our understanding of the diversity, evolution, and metabolic potential of archaea, yet the vast majority of archaeal diversity remains undersampled. In this work, we coordinated single-cell genomics with metagenomics in order to construct a near-complete genome from a deeply branching uncultivated archaeal lineage sampled from Great Boiling Spring (GBS) in the U.S. Great Basin, Nevada. This taxon is distantly related (distinct families) to an archaeal genome, designated "Novel Archaeal Group 1" (NAG1), which was extracted from a metagenome recovered from an acidic iron spring in Yellowstone National Park (YNP). We compared the metabolic predictions of the NAG1 lineage to better understand how these archaea could inhabit such chemically distinct environments. Similar to the NAG1 population previously studied in YNP, the NAG1 population from GBS is predicted to utilize proteins as a primary carbon source, ferment simple carbon sources, and use oxygen as a terminal electron acceptor under oxic conditions. However, GBS NAG1 populations contained distinct genes involved in central carbon metabolism and electron transfer, including nitrite reductase, which could confer the ability to reduce nitrite under anaerobic conditions. Despite inhabiting chemically distinct environments with large variations in pH, GBS NAG1 populations shared many core genomic and metabolic features with the archaeon identified from YNP, yet were able to carve out a distinct niche at GBS.

Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea.

We present two standards developed by the Genomic Standards Consortium (GSC) for reporting bacterial and archaeal genome sequences. Both are extensions of the Minimum Information about Any (x) Sequence (MIxS). The standards are the Minimum Information about a Single Amplified Genome (MISAG) and the Minimum Information about a Metagenome-Assembled Genome (MIMAG), including, but not limited to, assembly quality, and estimates of genome completeness and contamination. These standards can be used in combination with other GSC checklists, including the Minimum Information about a Genome Sequence (MIGS), Minimum Information about a Metagenomic Sequence (MIMS), and Minimum Information about a Marker Gene Sequence (MIMARKS). Community-wide adoption of MISAG and MIMAG will facilitate more robust comparative genomic analyses of bacterial and archaeal diversity.

High-Quality Draft Genome Sequence of Thermocrinis jamiesonii GBS1T Isolated from Great Boiling Spring, Nevada.

The draft genome of Thermocrinis jamiesonii GBS1T is 1,315,625 bp in 10 contigs and encodes 1,463 predicted genes. The presence of sox genes and various glycoside hydrolases and the absence of uptake NiFe hydrogenases (hyaB) are consistent with a requirement for thiosulfate and suggest the ability to use carbohydrate polymers.

Single-Cell-Genomics-Facilitated Read Binning of Candidate Phylum EM19 Genomes from Geothermal Spring Metagenomes.

The vast majority of microbial life remains uncatalogued due to the inability to cultivate these organisms in the laboratory. This "microbial dark matter" represents a substantial portion of the tree of life and of the populations that contribute to chemical cycling in many ecosystems. In this work, we leveraged an existing single-cell genomic data set representing the candidate bacterial phylum "Calescamantes" (EM19) to calibrate machine learning algorithms and define metagenomic bins directly from pyrosequencing reads derived from Great Boiling Spring in the U.S. Great Basin. Compared to other assembly-based methods, taxonomic binning with a read-based machine learning approach yielded final assemblies with the highest predicted genome completeness of any method tested. Read-first binning subsequently was used to extract Calescamantes bins from all metagenomes with abundant Calescamantes populations, including metagenomes from Octopus Spring and Bison Pool in Yellowstone National Park and Gongxiaoshe Spring in Yunnan Province, China. Metabolic reconstruction suggests that Calescamantes are heterotrophic, facultative anaerobes, which can utilize oxidized nitrogen sources as terminal electron acceptors for respiration in the absence of oxygen and use proteins as their primary carbon source. Despite their phylogenetic divergence, the geographically separate Calescamantes populations were highly similar in their predicted metabolic capabilities and core gene content, respiring O2, or oxidized nitrogen species for energy conservation in distant but chemically similar hot springs.

Phylogeny and physiology of candidate phylum 'Atribacteria' (OP9/JS1) inferred from cultivation-independent genomics.

The 'Atribacteria' is a candidate phylum in the Bacteria recently proposed to include members of the OP9 and JS1 lineages. OP9 and JS1 are globally distributed, and in some cases abundant, in anaerobic marine sediments, geothermal environments, anaerobic digesters and reactors and petroleum reservoirs. However, the monophyly of OP9 and JS1 has been questioned and their physiology and ecology remain largely enigmatic due to a lack of cultivated representatives. Here cultivation-independent genomic approaches were used to provide a first comprehensive view of the phylogeny, conserved genomic features and metabolic potential of members of this ubiquitous candidate phylum. Previously available and heretofore unpublished OP9 and JS1 single-cell genomic data sets were used as recruitment platforms for the reconstruction of atribacterial metagenome bins from a terephthalate-degrading reactor biofilm and from the monimolimnion of meromictic Sakinaw Lake. The single-cell genomes and metagenome bins together comprise six species- to genus-level groups that represent most major lineages within OP9 and JS1. Phylogenomic analyses of these combined data sets confirmed the monophyly of the 'Atribacteria' inclusive of OP9 and JS1. Additional conserved features within the 'Atribacteria' were identified, including a gene cluster encoding putative bacterial microcompartments that may be involved in aldehyde and sugar metabolism, energy conservation and carbon storage. Comparative analysis of the metabolic potential inferred from these data sets revealed that members of the 'Atribacteria' are likely to be heterotrophic anaerobes that lack respiratory capacity, with some lineages predicted to specialize in either primary fermentation of carbohydrates or secondary fermentation of organic acids, such as propionate.