Metagenomic sequencing of a Patescibacteria-containing enrichment from Zodletone spring in Oklahoma, USA

An enrichment of sulfidic sediments from Zodletone spring was sequenced as a metagenome. Draft genomes representing Cloacimonadota, Deltabacterota, Firmicutes, and Patescibacteria were binned and annotated and will aid functional genomics and cultivation efforts.

From wolves to humans: oral microbiome resistance to transfer across mammalian hosts

The mammalian mouth is colonized by complex microbial communities, adapted to specific niches, and in homeostasis with the host. Individual microbes interact metabolically and rely primarily on nutrients provided by the host, with which they have potentially co-evolved along the mammalian lineages. The oral environment is similar across mammals, but the diversity, specificity, and evolution of community structure in related or interacting mammals are little understood. Here, we compared the oral microbiomes of dogs with those of wild wolves and humans. In dogs, we found an increased microbial diversity relative to wolves, possibly related to the transition to omnivorous nutrition following domestication. This includes a larger diversity of Patescibacteria than previously reported in any other oral microbiota. The oral microbes are most distinct at bacterial species or strain levels, with few if any shared between humans and canids, while the close evolutionary relationship between wolves and dogs is reflected by numerous shared taxa. More taxa are shared at higher taxonomic levels including with humans, supporting their more ancestral common mammalian colonization followed by diversification. Phylogenies of selected oral bacterial lineages do not support stable human-dog microbial transfers but suggest diversification along mammalian lineages (apes and canids). Therefore, despite millennia of cohabitation and close interaction, the host and its native community controls and limits the assimilation of new microbes, even if closely related. Higher resolution metagenomic and microbial physiological studies, covering a larger mammalian diversity, should help understand how oral communities assemble, adapt, and interact with their hosts.IMPORTANCENumerous types of microbes colonize the mouth after birth and play important roles in maintaining oral health. When the microbiota-host homeostasis is perturbed, proliferation of some bacteria leads to diseases such as caries and periodontitis. Unlike the gut microbiome, the diversity of oral microbes across the mammalian evolutionary space is not understood. Our study compared the oral microbiomes of wild wolves, dogs, and apes (humans, chimpanzees, and bonobos), with the aim of identifying if microbes have been potentially exchanged between humans and dogs as a result of domestication and cohabitation. We found little if any evidence for such exchanges. The significance of our research is in finding that the oral microbiota and/or the host limit the acquisition of exogenous microbes, which is important in the context of natural exclusion of potential novel pathogens. We provide a framework for expanded higher-resolution studies across domestic and wild animals to understand resistance/resilience.

The gut microbiome mediates adaptation to scarce food in Coleoptera

Beetles are ubiquitous cave invertebrates worldwide that adapted to scarce subterranean resources when they colonized caves. Here, we investigated the potential role of gut microbiota in the adaptation of beetles to caves from different climatic regions of the Carpathians. The beetles' microbiota was host-specific, reflecting phylogenetic and nutritional adaptation. The microbial community structure further resolved conspecific beetles by caves suggesting microbiota-host coevolution and influences by local environmental factors. The detritivore species hosted a variety of bacteria known to decompose and ferment organic matter, suggesting turnover and host cooperative digestion of the sedimentary microbiota and allochthonous-derived nutrients. The cave Carabidae, with strong mandibula, adapted to predation and scavenging of animal and plant remains, had distinct microbiota dominated by symbiotic lineages Spiroplasma or Wolbachia. All beetles had relatively high levels of fermentative Carnobacterium and Vagococcus involved in lipid accumulation and a reduction of metabolic activity, and both features characterize adaptation to caves.

Unraveling Vulnerabilities in Endocrine Therapy-Resistant HER2+/ER+ Breast Cancer

Breast tumors overexpressing human epidermal growth factor receptor (HER2) confer intrinsic resistance to endocrine therapy (ET), and patients with HER2/estrogen receptor-positive (HER2+/ER+) breast cancer (BCa) are less responsive to ET than HER2-/ER+. However, real-world evidence reveals that a large subset of patients with HER2+/ER+ receive ET as monotherapy, positioning this treatment pattern as a clinical challenge. In the present study, we developed and characterized 2 in vitro models of ET-resistant (ETR) HER2+/ER+ BCa to identify possible therapeutic vulnerabilities. To mimic ETR to aromatase inhibitors (AIs), we developed 2 long-term estrogen deprivation (LTED) cell lines from BT-474 (BT474) and MDA-MB-361 (MM361). Growth assays, PAM50 subtyping, and genomic and transcriptomic analyses, followed by validation and functional studies, were used to identify targetable differences between ET-responsive parental and ETR-LTED HER2+/ER+ cells. Compared to their parental cells, MM361 LTEDs grew faster, lost ER, and increased HER2 expression, whereas BT474 LTEDs grew slower and maintained ER and HER2 expression. Both LTED variants had reduced responsiveness to fulvestrant. Whole-genome sequencing of aggressive MM361 LTEDs identified mutations in genes encoding transcription factors and chromatin modifiers. Single-cell RNA sequencing demonstrated a shift towards non-luminal phenotypes, and revealed metabolic remodeling of MM361 LTEDs, with upregulated lipid metabolism and ferroptosis-associated antioxidant genes, including GPX4. Combining a GPX4 inhibitor with anti-HER2 agents induced significant cell death in both MM361 and BT474 LTEDs. The BT474 and MM361 AI-resistant models capture distinct phenotypes of HER2+/ER+ BCa and identify altered lipid metabolism and ferroptosis remodeling as vulnerabilities of this type of ETR BCa.

Tectonic and geological setting influence hot spring microbiology

Hydrothermal systems form at divergent and convergent boundaries of lithospheric plates and within plates due to weakened crust and mantle plumes, playing host to diverse microbial ecosystems. Little is known of how differences in tectonic setting influence the geochemical and microbial compositions of these hydrothermal ecosystems. Here, coordinated geochemical and microbial community analyses were conducted on 87 high-temperature (>65°C) water and sediment samples from hot springs in Yellowstone National Park, Wyoming, USA (n = 41; mantle plume setting), Iceland (n = 41, divergent boundary), and Japan (n = 5; convergent boundary). Region-specific variation in geochemistry and sediment-associated 16S rRNA gene amplicon sequence variant (ASV) composition was observed, with 16S rRNA gene assemblages being nearly completely distinguished by region and pH being the most explanatory parameter within regions. Several low abundance ASVs exhibited cosmopolitan distributions across regions, while most high-abundance ASVs were only identified in specific regions. The presence of some cosmopolitan ASVs across regions argues against dispersal limitation primarily shaping the distribution of taxa among regions. Rather, the results point to local tectonic and geologic characteristics shaping the geochemistry of continental hydrothermal systems that then select for distinct microbial assemblages. These results provide new insights into the co-evolution of hydrothermal systems and their microbial communities.

Unraveling Vulnerabilities in Endocrine Therapy-Resistant HER2+/ER+ Breast Cancer

Background

Breast tumors overexpressing human epidermal growth factor receptor (HER2) confer intrinsic resistance to endocrine therapy (ET), and patients with HER2/ estrogen receptor-positive (HER2+/HR+) breast cancer (BCa) are less responsive to ET than HER2-/ER+. However, real-world evidence reveals that a large subset of HER2+/ER+ patients receive ET as monotherapy, positioning this treatment pattern as a clinical challenge. In the present study, we developed and characterized two distinct in vitro models of ET-resistant (ETR) HER2+/ER+ BCa to identify possible therapeutic vulnerabilities.

Methods

To mimic ETR to aromatase inhibitors (AI), we developed two long-term estrogen-deprived (LTED) cell lines from BT-474 (BT474) and MDA-MB-361 (MM361). Growth assays, PAM50 molecular subtyping, genomic and transcriptomic analyses, followed by validation and functional studies, were used to identify targetable differences between ET-responsive parental and ETR-LTED HER2+/ER+ cells.

Results

Compared to their parental cells, MM361 LTEDs grew faster, lost ER, and increased HER2 expression, whereas BT474 LTEDs grew slower and maintained ER and HER2 expression. Both LTED variants had reduced responsiveness to fulvestrant. Whole-genome sequencing of the more aggressive MM361 LTED model system identified exonic mutations in genes encoding transcription factors and chromatin modifiers. Single-cell RNA sequencing demonstrated a shift towards non-luminal phenotypes, and revealed metabolic remodeling of MM361 LTEDs, with upregulated lipid metabolism and antioxidant genes associated with ferroptosis, including GPX4. Combining the GPX4 inhibitor RSL3 with anti-HER2 agents induced significant cell death in both the MM361 and BT474 LTEDs.

Conclusions

The BT474 and MM361 AI-resistant models capture distinct phenotypes of HER2+/ER+ BCa and identify altered lipid metabolism and ferroptosis remodeling as vulnerabilities of this type of ETR BCa.

The gut microbiome mediates adaptation to scarce food in Coleoptera

Beetles are ubiquitous cave invertebrates worldwide that adapted to scarce subterranean resources when they colonized caves. Here, we investigated the potential role of gut microbiota in the adaptation of beetles to caves from different climatic regions of the Carpathians. The beetles' microbiota was host-specific, reflecting phylogenetic and nutritional adaptation. The microbial community structure further resolved conspecific beetles by caves suggesting microbiota-host coevolution and influences by local environmental factors. The detritivore species hosted a variety of bacteria known to decompose and ferment organic matter, suggesting turnover and host cooperative digestion of the sedimentary microbiota and allochthonous-derived nutrients. The cave Carabidae, with strong mandibulae adapted to predation and scavenging of animal and plant remains, had distinct microbiota dominated by symbiotic lineages Spiroplasma or Wolbachia . All beetles had relatively high levels of fermentative Carnobacterium and Vagococcus involved in lipid accumulation and a reduction of metabolic activity, and both features characterize adaptation to caves.

Complete Genome Sequence of Desulfomicrobium sp. Strain ZS1 from Zodletone Spring in Oklahoma, USA

Desulfomicrobium sp. strain ZS1 is an obligate anaerobic, sulfate-reducing member of the Desulfobacterota from Zodletone Spring, an anoxic sulfide-rich spring in southwestern Oklahoma. Its complete genome was sequenced using a combination of Illumina and Oxford Nanopore platforms and encodes 3,364 proteins and 81 RNAs on a single chromosome.

Transcriptional Control of hgcAB by an ArsR-Like Regulator in Pseudodesulfovibrio mercurii ND132

The hgcAB gene pair encodes mercury (Hg) methylation capability in a diverse group of microorganisms, but its evolution and transcriptional regulation remain unknown. Working from the possibility that the evolutionary function of HgcAB may not be Hg methylation, we test a possible link to arsenic resistance. Using model Hg methylator Pseudodesulfovibrio mercurii ND132, we evaluated transcriptional control of hgcAB by a putative ArsR encoded upstream and cotranscribed with hgcAB. This regulator shares homology with ArsR repressors of arsenic resistance and S-adenosylhomocysteine (SAH)-responsive regulators of methionine biosynthesis but is distinct from other ArsR/SahR proteins in P. mercurii. Using quantitative PCR (qPCR) and RNA sequencing (RNA-seq) transcriptome analyses, we confirmed this ArsR regulates hgcAB transcription and is responsive to arsenic and SAH. Additionally, RNA-seq indicated a possible link between hgcAB activity and arsenic transformations, with significant upregulation of other ArsR-regulated arsenic resistance operons alongside hgcAB. Interestingly, wild-type ND132 was less sensitive to As(V) (but not As(III)) than an hgcAB knockout strain, supporting the idea that hgcAB may be linked to arsenic resistance. Arsenic significantly impacted rates of Hg methylation by ND132; however, responses varied with culture conditions. Differences in growth and metabolic activity did not account for arsenic impacts on methylation. While arsenic significantly increased hgcAB expression, hgcAB gene and transcript abundance was not a good predictor of Hg methylation rates. Taken together, these results support the idea that Hg and As cycling are linked in P. mercurii ND132. Our results may hold clues to the evolution of hgcAB and the controls on Hg methylation in nature. IMPORTANCE This work reveals a link between microbial mercury methylation and arsenic resistance and may hold clues to the evolution of mercury methylation genes (hgcAB). Microbes with hgcAB produce methylmercury, a strong neurotoxin that readily accumulates in the food web. This study addresses a critical gap in our understanding about the environmental factors that control hgcAB expression. We show that hgcAB expression is controlled by an ArsR-like regulator responsive to both arsenic and S-adenosylhomocysteine in our model organism, Pseudodesulfovibrio mercurii ND132. Exposure to arsenic also significantly impacted Pseudodesulfovibrio mercurii ND132 mercury methylation rates. However, expression of hgcAB was not always a good predictor of Hg methylation rates, highlighting the roles of Hg bioavailability and other biochemical mechanisms in methylmercury production. This study improves our understanding of the controls on hgcAB expression, which is needed to better predict environmental methylmercury production.

Benchmarking second and third-generation sequencing platforms for microbial metagenomics

Shotgun metagenomic sequencing is a common approach for studying the taxonomic diversity and metabolic potential of complex microbial communities. Current methods primarily use second generation short read sequencing, yet advances in third generation long read technologies provide opportunities to overcome some of the limitations of short read sequencing. Here, we compared seven platforms, encompassing second generation sequencers (Illumina HiSeq 300, MGI DNBSEQ-G400 and DNBSEQ-T7, ThermoFisher Ion GeneStudio S5 and Ion Proton P1) and third generation sequencers (Oxford Nanopore Technologies MinION R9 and Pacific Biosciences Sequel II). We constructed three uneven synthetic microbial communities composed of up to 87 genomic microbial strains DNAs per mock, spanning 29 bacterial and archaeal phyla, and representing the most complex and diverse synthetic communities used for sequencing technology comparisons. Our results demonstrate that third generation sequencing have advantages over second generation platforms in analyzing complex microbial communities, but require careful sequencing library preparation for optimal quantitative metagenomic analysis. Our sequencing data also provides a valuable resource for testing and benchmarking bioinformatics software for metagenomics.

Incorporating concentration-dependent sediment microbial activity into methylmercury production kinetics modeling

In anoxic environments, anaerobic microorganisms carrying the hgcAB gene cluster can mediate the transformation of inorganic mercury (Hg(II)) to monomethylmercury (MMHg). The kinetics of Hg(II) transformation to MMHg in periphyton from East Fork Poplar Creek (EFPC) in Oak Ridge, TN have previously been modeled using a transient availability model (TAM). The TAM for Hg(II) methylation combines methylation/demethylation kinetics with kinetic expressions for processes that decrease Hg(II) and MMHg availability for methylation and demethylation (multisite sorption of Hg(II) and MMHg, Hg(II) reduction/Hg(0) oxidation). In this study, the TAM is used for the first time to describe MMHg production in sediment. We assessed MMHg production in sediment microcosms using two different sediment types from EFPC: a relatively anoxic, carbon-rich sediment with higher microbial activity (higher CO₂ production from sediment) and a relatively oxic, sandy, carbon-poor sediment with lower microbial activity (lower CO₂ production from sediment). Based on 16s rRNA sequencing, the overall microbial community structure in the two sediments was retained during the incubations. However, the hgcA containing methanogenic Euryarchaeota communities differed between sediment types and their growth followed different trajectories over the course of incubations, potentially contributing to the distinct patterns of MMHg production observed. The general TAM paradigm performed well in describing MMHg production in the sediments. However, the MMHg production and ancillary data suggested the need to revise the model structure to incorporate terms for concentration-dependent microbial activity over the course of the incubations. We modified the TAM to include Monod-type kinetics for methylation and demethylation and observed an improved fit for the carbon-rich, microbially active sediment. Overall our work shows that the TAM can be applied to describe Hg(II) methylation in sediments and that including expressions accounting for concentration-dependent microbial activity can improve the accuracy of the model description of the data in some cases.

Morphine and high-fat diet differentially alter the gut microbiota composition and metabolic function in lean versus obese mice

There are known associations between opioids, obesity, and the gut microbiome, but the molecular connection/mediation of these relationships is not understood. To better clarify the interplay of physiological, genetic, and microbial factors, this study investigated the microbiome and host inflammatory responses to chronic opioid administration in genetically obese, diet-induced obese, and lean mice. Samples of feces, urine, colon tissue, and plasma were analyzed using targeted LC-MS/MS quantification of metabolites, immunoassays of inflammatory cytokine levels, genome-resolved metagenomics, and metaproteomics. Genetic obesity, diet-induced obesity, and morphine treatment in lean mice each showed increases in distinct inflammatory cytokines. Metagenomic assembly and binning uncovered over 400 novel gut bacterial genomes and species. Morphine administration impacted the microbiome's composition and function, with the strongest effect observed in lean mice. This microbiome effect was less pronounced than either diet or genetically driven obesity. Based on inferred microbial physiology from the metaproteome datasets, a high-fat diet transitioned constituent microbes away from harvesting diet-derived nutrients and towards nutrients present in the host mucosal layer. Considered together, these results identified novel host-dependent phenotypes, differentiated the effects of genetic obesity versus diet induced obesity on gut microbiome composition and function, and showed that chronic morphine administration altered the gut microbiome.

A consensus protocol for the recovery of mercury methylation genes from metagenomes

Mercury (Hg) methylation genes (hgcAB) mediate the formation of the toxic methylmercury and have been identified from diverse environments, including freshwater and marine ecosystems, Arctic permafrost, forest and paddy soils, coal-ash amended sediments, chlor-alkali plants discharges and geothermal springs. Here we present the first attempt at a standardized protocol for the detection, identification and quantification of hgc genes from metagenomes. Our Hg-MATE (Hg-cycling Microorganisms in Aquatic and Terrestrial Ecosystems) database, a catalogue of hgc genes, provides the most accurate information to date on the taxonomic identity and functional/metabolic attributes of microorganisms responsible for Hg methylation in the environment. Furthermore, we introduce "marky-coco", a ready-to-use bioinformatic pipeline based on de novo single-metagenome assembly, for easy and accurate characterization of hgc genes from environmental samples. We compared the recovery of hgc genes from environmental metagenomes using the marky-coco pipeline with an approach based on co-assembly of multiple metagenomes. Our data show similar efficiency in both approaches for most environments except those with high diversity (i.e., paddy soils) for which a co-assembly approach was preferred. Finally, we discuss the definition of true hgc genes and methods to normalize hgc gene counts from metagenomes.

Abstract 1776: Loss of estrogen receptor alters drug responsiveness and supports a basal-like phenotype in endocrine-resistant HER2+/ER+ breast cancer

At least 50% of human epidermal growth factor 2 (HER2) enriched breast cancer (BC) is estrogen receptor positive (ER+). In clinical trials, patients with HER2+ BC expressing ER respond poorly to neoadjuvant anti-HER2 therapy versus those lacking ER expression. Additionally, combining anti-HER2 with endocrine therapy (ET), like aromatase inhibitor (AI), does not significantly improve pathological complete response in HER2+/ER+ patients. In advanced or de novo HER2+/ER+ metastatic BC (MBC), treatment with ET plus anti-HER2 yields the highest 5-year overall survival. Unfortunately, in the real-world setting outside of clinical trials, just 23% of patients receive this combination, and nearly 40% receive ET alone. The clinical challenge associated with treating advanced HER2+/ER+ BC demands a better understanding of ER signaling and downstream pathways in HER2-driven BC. Our goal is to identify alternative targeted therapies and/or tailor therapeutic combinatorial strategies to treat and/or delay progression of HER2+/ER+ BC. In this study, we established two long-term estrogen (E2) deprived (LTED) HER2+/ER+ cell line models from BT-474 (BT) and MDA-MB-361 (MM) to mimic endocrine therapy resistance (ETR) to AI, and characterized the response of these cell lines to anti-HER2 and other ETs. We also performed whole-genome sequencing (WGS) and single-cell RNA sequencing (scRNAseq) to identify differences between parental and derived ETR-LTEDs. The ER expression was retained in BT-LTEDs but lost in MM-LTEDs as compared to parental cells. Additionally, ER transcriptional activity was confirmed by increased expression of ER-target genes upon E2 stimulation only in BT-LTEDs. However, both LTED variants were less responsive to fulvestrant and showed upregulation of pro-survival AKT signaling. Focusing on the now ER- MM-LTED model, WGS and targeted duplex sequencing identified a significant increase in exonic missense mutations, notably C>T and C>A. Several of the mutated genes encode for transcription and chromatin regulatory factors. scRNAseq analysis showed that MM-LTEDs shift from a luminal- to more basal-like phenotype, with enrichment of genes that regulate immune response and cell motility. Thus, loss of ER expression in MM-LTEDs on the mRNA and protein levels may explain observed intrinsic resistance to fulvestrant and gain of the basal-like phenotype. In summary, we report that our ETR BT- and MM-LTED models capture distinct ER-associated phenotypes of HER2+/ER+ BC. Additional studies are necessary to understand the functional significance of these missense mutations on the development of drug resistance in HER2+/ER+ BC. Ongoing studies are testing pharmacological inhibition of key upregulated genes associated with the basal phenotype to potentially delay disease progression and provide better treatments for ETR HER2+/ER+ BC.

Citation Format: Shaymaa Bahnassy, Hillary Stires, Lu Jin, Stanley Tam, Yasmin J. Dunn, Brendan F. Kohrn, Lawrence A. Loeb, Manasi Balachandran, Mircea Podar, Matthew D. McCoy, Robert A. Beckman, Rebecca B. Riggins. Loss of estrogen receptor alters drug responsiveness and supports a basal-like phenotype in endocrine-resistant HER2+/ER+ breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1776.

Pseudodesulfovibrio mercurii sp. nov., a mercury-methylating bacterium isolated from sediment

The sulfate-reducing, mercury-methylating strain ND132^T was isolated from the brackish anaerobic bottom sediments of Chesapeake Bay, USA. Capable of high levels of mercury (Hg) methylation, ND132^T has been widely used as a model strain to study the process and to determine the genetic basis of Hg methylation. Originally called Desulfovibrio desulfuricans ND132^T on the basis of an early partial 16S rRNA sequence, the strain has never been formally described. Phylogenetic and physiological traits place this strain within the genus Pseudodesulfovibrio, in the recently reclassified phylum Desulfobacterota (formerly Deltaproteobacteria). ND132^T is most closely related to Pseudodesulfovibrio hydrargyri BerOc1^T and Pseudodesulfovibrio indicus J2^T. Analysis of average nucleotide identity (ANI) of whole-genome sequences showed roughly 88 % ANI between P. hydrargyri BerOc1^T and ND132^T, and 84 % similarity between ND132^T and P. indicus J2^T. These cut-off scores <95 %, along with a multi-gene phylogenetic analysis of members of the family Desulfovibrionacea, and differences in physiology indicate that all three strains represent separate species. The Gram-stain-negative cells are vibrio-shaped, motile and not sporulated. ND132^T is a salt-tolerant mesophile with optimal growth in the laboratory at 32 °C, 2 % salinity, and pH 7.8. The DNA G+C content of the genomic DNA is 65.2 %. It is an incomplete oxidizer of short chain fatty acids, using lactate, pyruvate and fumarate with sulfate or sulfite as the terminal electron acceptors. ND132^T can respire fumarate using pyruvate as an electron donor. The major fatty acids are iso-C_15 : 0, anteiso-C_15 : 0, iso-C_17 : 0, iso-C_17 : 1ω9c and anteiso-C_17 : 0. We propose the classification of strain ND132^T (DSM 110689, ATCC TSD-224) as the type strain Pseudodesulfovibrio mercurii sp. nov.

An Improved hgcAB Primer Set and Direct High-Throughput Sequencing Expand Hg-Methylator Diversity in Nature

The gene pair hgcAB is essential for microbial mercury methylation. Our understanding of its abundance and diversity in nature is rapidly evolving. In this study we developed a new broad-range primer set for hgcAB, plus an expanded hgcAB reference library, and used these to characterize Hg-methylating communities from diverse environments. We applied this new Hg-methylator database to assign taxonomy to hgcA sequences from clone, amplicon, and metagenomic datasets. We evaluated potential biases introduced in primer design, sequence length, and classification, and suggest best practices for studying Hg-methylator diversity. Our study confirms the emerging picture of an expanded diversity of HgcAB-encoding microbes in many types of ecosystems, with abundant putative mercury methylators Nitrospirae and Chloroflexi in several new environments including salt marsh and peat soils. Other common microbes encoding HgcAB included Phycisphaerae, Aminicenantes, Spirochaetes, and Elusimicrobia. Combined with high-throughput amplicon specific sequencing, the new primer set also indentified novel hgcAB sequences similar to Lentisphaerae, Bacteroidetes, Atribacteria, and candidate phyla WOR-3 and KSB1 bacteria. Gene abundance data also corroborate the important role of two "classic" groups of methylators (Deltaproteobacteria and Methanomicrobia) in many environments, but generally show a scarcity of hgcAB+ Firmicutes. The new primer set was developed to specifically target hgcAB sequences found in nature, reducing degeneracy and providing increased sensitivity while maintaining broad diversity capture. We evaluated mock communities to confirm primer improvements, including culture spikes to environmental samples with variable DNA extraction and PCR amplification efficiencies. For select sites, this new workflow was combined with direct high-throughput hgcAB sequencing. The hgcAB diversity generated by direct amplicon sequencing confirmed the potential for novel Hg-methylators previously identified using metagenomic screens. A new phylogenetic analysis using sequences from freshwater, saline, and terrestrial environments showed Deltaproteobacteria HgcA sequences generally clustered among themselves, while metagenome-resolved HgcA sequences in other phyla tended to cluster by environment, suggesting horizontal gene transfer into many clades. HgcA from marine metagenomes often formed distinct subtrees from those sequenced from freshwater ecosystems. Overall the majority of HgcA sequences branch from a cluster of HgcAB fused proteins related to Thermococci, Atribacteria (candidate division OP9), Aminicenantes (OP8), and Chloroflexi. The improved primer set and library, combined with direct amplicon sequencing, provide a significantly improved assessment of the abundance and diversity of hgcAB+ microbes in nature.

Characterization of subsurface media from locations up- and down-gradient of a uranium-contaminated aquifer

The processing of sediment to accurately characterize the spatially-resolved depth profiles of geophysical and geochemical properties along with signatures of microbial density and activity remains a challenge especially in complex contaminated areas. This study processed cores from two sediment boreholes from background and contaminated core sediments and surrounding groundwater. Fresh core sediments were compared by depth to capture the changes in sediment structure, sediment minerals, biomass, and pore water geochemistry in terms of major and trace elements including pollutants, cations, anions, and organic acids. Soil porewater samples were matched to groundwater level, flow rate, and preferential flows and compared to homogenized groundwater-only samples from neighboring monitoring wells. Groundwater analysis of nearby wells only revealed high sulfate and nitrate concentrations while the same analysis using sediment pore water samples with depth was able to suggest areas high in sulfate- and nitrate-reducing bacteria based on their decreased concentration and production of reduced by-products that could not be seen in the groundwater samples. Positive correlations among porewater content, total organic carbon, trace metals and clay minerals revealed a more complicated relationship among contaminant, sediment texture, groundwater table, and biomass. The fluctuating capillary interface had high concentrations of Fe and Mn-oxides combined with trace elements including U, Th, Sr, Ba, Cu, and Co. This suggests the mobility of potentially hazardous elements, sediment structure, and biogeochemical factors are all linked together to impact microbial communities, emphasizing that solid interfaces play an important role in determining the abundance of bacteria in the sediments.

Comparative Analysis of Microbial Diversity Across Temperature Gradients in Hot Springs From Yellowstone and Iceland

Geothermal hot springs are a natural setting to study microbial adaptation to a wide range of temperatures reaching up to boiling. Temperature gradients lead to distinct microbial communities that inhabit their optimum niches. We sampled three alkaline, high temperature (80-100°C) hot springs in Yellowstone and Iceland that had cooling outflows and whose microbial communities had not been studied previously. The microbial composition in sediments and mats was determined by DNA sequencing of rRNA gene amplicons. Over three dozen phyla of Archaea and Bacteria were identified, representing over 1700 distinct organisms. We observed a significant non-linear reduction in the number of microbial taxa as the temperature increased from warm (38°C) to boiling. At high taxonomic levels, the community structure was similar between the Yellowstone and Iceland hot springs. We identified potential endemism at the genus level, especially in thermophilic phototrophs, which may have been potentially driven by distinct environmental conditions and dispersal limitations.

Structure determination of the HgcAB complex using metagenome sequence data: insights into microbial mercury methylation

Bacteria and archaea possessing the hgcAB gene pair methylate inorganic mercury (Hg) to form highly toxic methylmercury. HgcA consists of a corrinoid binding domain and a transmembrane domain, and HgcB is a dicluster ferredoxin. However, their detailed structure and function have not been thoroughly characterized. We modeled the HgcAB complex by combining metagenome sequence data mining, coevolution analysis, and Rosetta structure calculations. In addition, we overexpressed HgcA and HgcB in Escherichia coli, confirmed spectroscopically that they bind cobalamin and [4Fe-4S] clusters, respectively, and incorporated these cofactors into the structural model. Surprisingly, the two domains of HgcA do not interact with each other, but HgcB forms extensive contacts with both domains. The model suggests that conserved cysteines in HgcB are involved in shuttling Hg^II, methylmercury, or both. These findings refine our understanding of the mechanism of Hg methylation and expand the known repertoire of corrinoid methyltransferases in nature.

Connor J. Cooper et al. expressed HgcA and HgcB in Escherichia coli and modeled the structure of the HgcAB complex by combining metagenome sequence data, coevolution analysis, and ab initio structure calculations. This study provides insights into the biochemical mechanism of mercury (Hg) methylation.

A Microbe Associated with Sleep Revealed by a Novel Systems Genetic Analysis of the Microbiome in Collaborative Cross Mice

The microbiome influences health and disease through complex networks of host genetics, genomics, microbes, and environment. Identifying the mechanisms of these interactions has remained challenging. Systems genetics in laboratory mice (Mus musculus) enables data-driven discovery of biological network components and mechanisms of host-microbial interactions underlying disease phenotypes. To examine the interplay among the whole host genome, transcriptome, and microbiome, we mapped QTL and correlated the abundance of cecal messenger RNA, luminal microflora, physiology, and behavior in a highly diverse Collaborative Cross breeding population. One such relationship, regulated by a variant on chromosome 7, was the association of Odoribacter (Bacteroidales) abundance and sleep phenotypes. In a test of this association in the BKS.Cg-Dock7m +/+ Leprdb/J mouse model of obesity and diabetes, known to have abnormal sleep and colonization by Odoribacter, treatment with antibiotics altered sleep in a genotype-dependent fashion. The many other relationships extracted from this study can be used to interrogate other diseases, microbes, and mechanisms.

Single-Cell Genomics and the Oral Microbiome

The human oral cavity is one of the first environments where microbes have been discovered and studied since the dawn of microbiology. Nevertheless, approximately 200 types of bacteria from the oral microbiota have remained uncultured in the laboratory. Some are associated with a healthy oral microbial community, while others are linked to oral diseases, from dental caries to gum disease. Single-cell genomics has enabled inferences on the physiology, virulence, and evolution of such uncultured microorganisms and has further enabled isolation and cultivation of several novel oral bacteria, including the discovery of novel interspecies interactions. This review summarizes some of the more recent advances in this field, which is rapidly moving toward physiologic characterization of single cells and ultimately cultivation of the yet uncultured. A combination of traditional microbiological approaches with genomic-based physiologic predictions and isolation strategies may lead to the oral microbiome being the first complex microbial community to have all its members cultivable in the laboratory. Studying the biology of the individual microbes when in association with other members of the community, in controlled laboratory conditions and in vivo, should lead to a better understanding of oral dysbiosis and its prevention and reversion.

Targeted isolation and cultivation of uncultivated bacteria by reverse genomics

Most microorganisms from all taxonomic levels are uncultured. Single-cell genomes and metagenomes continue to increase the known diversity of Bacteria and Archaea, but while ‘omics can be used to infer physiological or ecological roles for species in a community, most of those hypothetical roles remain unvalidated. Here we report an approach to capture specific microorganisms from complex communities into pure cultures using genome-informed antibody engineering. We apply our reverse genomics approach to isolate and sequence single cells and to cultivate three different species-level lineages of human oral Saccharibacteria/TM7. Using our pure cultures we show that all three saccharibacteria species are epibionts of diverse Actinobacteria. We also isolate and cultivate human oral SR1 bacteria, which are members of a lineage of previously uncultured bacteria. Reverse-genomics-enabled cultivation of microorganisms can be applied to any species from any environment and has the potential to unlock the isolation, cultivation and characterization of species from as-yet-uncultured branches of the microbial tree of life.

Determining the Reliability of Measuring Mercury Cycling Gene Abundance with Correlations with Mercury and Methylmercury Concentrations

Methylmercury (MeHg) is a bioaccumulative toxic contaminant in many ecosystems, but factors governing its production are poorly understood. Recent work has shown that the anaerobic microbial conversion of mercury (Hg) to MeHg requires the Hg-methylation genes hgcAB and that these genes can be used as biomarkers in PCR-based estimators of Hg-methylator abundance. In an effort to determine reliable methods for assessing hgcA abundance and diversity and linking them to MeHg concentrations, multiple approaches were compared including metagenomic shotgun sequencing, 16S rRNA gene pyrosequencing and cloning/sequencing hgcAB gene products. Hg-methylator abundance was also determined by quantitative hgcA qPCR amplification and metaproteomics for comparison to the above measurements. Samples from eight sites were examined covering a range of total Hg (HgT; 0.03-14 mg kg^-1 dry wt. soil) and MeHg (0.05-27 μg kg^-1 dry wt. soil) concentrations. In the metagenome and amplicon sequencing of hgcAB diversity, the Deltaproteobacteria were the dominant Hg-methylators while Firmicutes and methanogenic Archaea were typically ∼50% less abundant. This was consistent with metaproteomics estimates where the Deltaproteobacteria were steadily higher. The 16S rRNA gene pyrosequencing did not have sufficient resolution to identify hgcAB⁺ species. Metagenomic and hgcAB results were similar for Hg-methylator diversity and clade-specific qPCR-based approaches for hgcA are only appropriate when comparing the abundance of a particular clade across various samples. Weak correlations between Hg-methylating bacteria and soil Hg concentrations were observed for similar environmental samples, but overall total Hg and MeHg concentrations poorly correlated with Hg-cycling genes.

Toward unrestricted use of public genomic data

Publication interests should not limit access to public data

Single-cell genomics of co-sorted Nanoarchaeota suggests novel putative host associations and diversification of proteins involved in symbiosis

BACKGROUND

Nanoarchaeota are obligate symbionts of other Archaea first discovered 16 years ago, yet little is known about this largely uncultivated taxon. While Nanoarchaeota diversity has been detected in a variety of habitats using 16S rRNA gene surveys, genome sequences have been available for only three Nanoarchaeota and their hosts. The host range and adaptation of Nanoarchaeota to a wide range of environmental conditions has thus largely remained elusive. Single-cell genomics is an ideal approach to address these questions as Nanoarchaeota can be isolated while still attached to putative hosts, enabling the exploration of cell-cell interactions and fine-scale genomic diversity.

RESULTS

From 22 single amplified genomes (SAGs) from three hot springs in Yellowstone National Park, we derived a genome-based phylogeny of the phylum Nanoarchaeota, linking it to global 16S rRNA gene diversity. By exploiting sequencing of co-sorted tightly attached cells, we associated Nanoarchaeota with 6 novel putative hosts, 2 of which were found in multiple SAGs, and showed that the same host species may associate with multiple species of Nanoarchaeota. Comparison of single nucleotide polymorphisms (SNPs) within a population of Nanoarchaeota SAGs indicated that Nanoarchaeota attached to a single host cell in situ are likely clonal. In addition to an overall pattern of purifying selection, we found significantly higher densities of non-synonymous SNPs in hypothetical cell surface proteins, as compared to other functional categories. Genes implicated in interactions in other obligate microbe-microbe symbioses, including those encoding a cytochrome bd-I ubiquinol oxidase and a FlaJ/TadC homologue possibly involved in type IV pili production, also had relatively high densities of non-synonymous SNPs.

CONCLUSIONS

This population genetics study of Nanoarchaeota greatly expands the known potential host range of the phylum and hints at what genes may be involved in adaptation to diverse environments or different hosts. We provide the first evidence that Nanoarchaeota cells attached to the same host cell are clonal and propose a hypothesis for how clonality may occur despite diverse symbiont populations.

Development and characterization of stable anaerobic thermophilic methanogenic microbiomes fermenting switchgrass at decreasing residence times

BACKGROUND

Anaerobic fermentation of lignocellulose occurs in both natural and managed environments, and is an essential part of the carbon cycle as well as a promising route to sustainable production of fuels and chemicals. Lignocellulose solubilization by mixed microbiomes is important in these contexts.

RESULTS

Here, we report the development of stable switchgrass-fermenting enrichment cultures maintained at various residence times and moderately high (55 °C) temperatures. Anaerobic microbiomes derived from a digester inoculum were incubated at 55 °C and fed semi-continuously with medium containing 30 g/L mid-season harvested switchgrass to achieve residence times (RT) of 20, 10, 5, and 3.3 days. Stable, time-invariant cellulolytic methanogenic cultures with minimal accumulation of organic acids were achieved for all RTs. Fractional carbohydrate solubilization was 0.711, 0.654, 0.581 and 0.538 at RT = 20, 10, 5 and 3.3 days, respectively, and glucan solubilization was proportional to xylan solubilization at all RTs. The rate of solubilization was described well by the equation r = k(C - C₀f_r), where C represents the concentration of unutilized carbohydrate, C₀ is the concentration of carbohydrate (cellulose and hemicellulose) entering the bioreactor and f_r is the extrapolated fraction of entering carbohydrate that is recalcitrant at infinite residence time. The 3.3 day RT is among the shortest RT reported for stable thermophilic, methanogenic digestion of a lignocellulosic feedstock. 16S rDNA phylotyping and metagenomic analyses were conducted to characterize the effect of RT on community dynamics and to infer functional roles in the switchgrass to biogas conversion to the various microbial taxa. Firmicutes were the dominant phylum, increasing in relative abundance from 54 to 96% as RT decreased. A Clostridium clariflavum strain with genetic markers for xylose metabolism was the most abundant lignocellulose-solubilizing bacterium. A Thermotogae (Defluviitoga tunisiensis) was the most abundant bacterium in switchgrass digesters at RT = 20 days but decreased in abundance at lower RTs as did multiple Chloroflexi. Synergistetes and Euryarchaeota were present at roughly constant levels over the range of RTs examined.

CONCLUSIONS

A system was developed in which stable methanogenic steady-states were readily obtained with a particulate biomass feedstock, mid-season switchgrass, at laboratory (1 L) scale. Characterization of the extent and rate of carbohydrate solubilization in combination with 16S rDNA and metagenomic sequencing provides a multi-dimensional view of performance, species composition, glycoside hydrolases, and metabolic function with varying residence time. These results provide a point of reference and guidance for future studies and organism development efforts involving defined cultures.

A new symbiotic nanoarchaeote (Candidatus Nanoclepta minutus) and its host (Zestosphaera tikiterensis gen. nov., sp. nov.) from a New Zealand hot spring

Three thermophilic Nanoarchaeota-Crenarchaeota symbiotic systems have been described. We obtained another stable anaerobic enrichment culture at 80°C, pH 6.0 from a New Zealand hot spring. The nanoarchaeote (Ncl-1) and its host (NZ3^T) were isolated in co-culture and their genomes assembled. The small (∼200nm) flagellated cocci were often attached to larger cocci. Based on 16S rRNA gene similarity (88.4%) and average amino acid identity (52%), Ncl-1 is closely related to Candidatus Nanopusillus acidilobi. Their genomes both encode for archaeal flagella and partial glycolysis and gluconeogenesis pathways, but lack ATP synthase genes. Like Nanoarchaeum equitans, Ncl-1 has a CRISPR-Cas system. Ncl-1 also relies on its crenarchaeotal host for most of its biosynthetic needs. The host NZ3^T was isolated and grows on proteinaceous substrates but not on sugars, alcohols, or fatty acids. NZ3^T requires thiosulfate and grows best at 82°C, pH 6.0. NZ3^T is most closely related to the Desulfurococcaceae, Ignisphaera aggregans (∼92% 16S rRNA gene sequence similarity, 45% AAI). Based on phylogenetic, physiological and genomic data, Ncl-1 and NZ3^T represent novel genera in the Nanoarchaeota and the Desulfurococcaceae, respectively, with the proposed names Candidatus Nanoclepta minutus and Zestosphaera tikiterensis gen. nov., sp. nov., type strain NZ3^T (=DSMZ 107634^T=OCM 1213^T).

Genomics of the Uncultivated, Periodontitis-Associated Bacterium Tannerella sp. BU045 (Oral Taxon 808)

Despite decades of research into the human oral microbiome, many species remain uncultivated. The technique of single-cell whole-genome amplification and sequencing provides a means of deriving genome sequences for species that can be informative on biological function and suggest pathways to cultivation. Tannerella forsythia has long been known to be highly associated with chronic periodontitis and to cause periodontitis-like symptoms in experimental animals, and Tannerella sp. BU045 (human oral taxon 808) is an uncultivated relative of this organism. In this work, we extend our previous sequencing of the Tannerella sp. BU063 (human oral taxon 286) genome by sequencing amplified genomes from 11 cells of Tannerella sp. BU045, including 3 genomes that are at least 90% complete. Tannerella sp. BU045 is more closely related to Tannerella sp. BU063 than to T. forsythia by gene content and average nucleotide identity. However, two independent data sets of association with periodontitis, one based on 16S rRNA gene abundance and the other based on gene expression in a metatranscriptomic data set, show that Tannerella sp. BU045 is more highly associated with disease than Tannerella sp. BU063. Comparative genomics shows genes and functions that are shared or unique to the different species, which may direct further research of the pathogenesis of chronic periodontitis. IMPORTANCE Periodontitis (gum disease) affects 47% of adults over 30 in the United States (P. I. Eke, B. A. Dye, L. Wei, G. O. Thornton-Evans, R. J. Genco, et al., J Dent Res 91:914-920, 2012), and it cost between $39 and $396 billion worldwide in 2015 (A. J. Righolt, M. Jevdjevic, W. Marcenes, and S. Listl, J Dent Res, 17 January 2018, https://doi.org/10.1177/0022034517750572). Many bacteria associated with the disease are known only by the DNA sequence of their 16S rRNA gene. In this publication, amplification and sequencing of DNA from single bacterial cells are used to obtain nearly complete genomes of Tannerella sp. BU045, a species of bacteria that is more prevalent in patients with periodontitis than in healthy patients. Comparing the complete genome of this bacterium to genomes of related bacterial species will help to better understand periodontitis and may help to grow this organism in pure culture, which would allow a better understanding of its role in the mouth.

Genus-Wide Assessment of Lignocellulose Utilization in the Extremely Thermophilic Genus Caldicellulosiruptor by Genomic, Pangenomic, and Metagenomic Analyses

Metagenomic data from Obsidian Pool (Yellowstone National Park, USA) and 13 genome sequences were used to reassess genus-wide biodiversity for the extremely thermophilic Caldicellulosiruptor The updated core genome contains 1,401 ortholog groups (average genome size for 13 species = 2,516 genes). The pangenome, which remains open with a revised total of 3,493 ortholog groups, encodes a variety of multidomain glycoside hydrolases (GHs). These include three cellulases with GH48 domains that are colocated in the glucan degradation locus (GDL) and are specific determinants for microcrystalline cellulose utilization. Three recently sequenced species, Caldicellulosiruptor sp. strain Rt8.B8 (renamed here Caldicellulosiruptor morganii), Thermoanaerobacter cellulolyticus strain NA10 (renamed here Caldicellulosiruptor naganoensis), and Caldicellulosiruptor sp. strain Wai35.B1 (renamed here Caldicellulosiruptor danielii), degraded Avicel and lignocellulose (switchgrass). C. morganii was more efficient than Caldicellulosiruptor bescii in this regard and differed from the other 12 species examined, both based on genome content and organization and in the specific domain features of conserved GHs. Metagenomic analysis of lignocellulose-enriched samples from Obsidian Pool revealed limited new information on genus biodiversity. Enrichments yielded genomic signatures closely related to that of Caldicellulosiruptor obsidiansis, but there was also evidence for other thermophilic fermentative anaerobes (Caldanaerobacter, Fervidobacterium, Caloramator, and Clostridium). One enrichment, containing 89.8% Caldicellulosiruptor and 9.7% Caloramator, had a capacity for switchgrass solubilization comparable to that of C. bescii These results refine the known biodiversity of Caldicellulosiruptor and indicate that microcrystalline cellulose degradation at temperatures above 70°C, based on current information, is limited to certain members of this genus that produce GH48 domain-containing enzymes.IMPORTANCE The genus Caldicellulosiruptor contains the most thermophilic bacteria capable of lignocellulose deconstruction, which are promising candidates for consolidated bioprocessing for the production of biofuels and bio-based chemicals. The focus here is on the extant capability of this genus for plant biomass degradation and the extent to which this can be inferred from the core and pangenomes, based on analysis of 13 species and metagenomic sequence information from environmental samples. Key to microcrystalline hydrolysis is the content of the glucan degradation locus (GDL), a set of genes encoding glycoside hydrolases (GHs), several of which have GH48 and family 3 carbohydrate binding module domains, that function as primary cellulases. Resolving the relationship between the GDL and lignocellulose degradation will inform efforts to identify more prolific members of the genus and to develop metabolic engineering strategies to improve this characteristic.

Microbial community structure with trends in methylation gene diversity and abundance in mercury-contaminated rice paddy soils in Guizhou, China

Paddy soils from mercury (Hg)-contaminated rice fields in Guizhou, China were studied with respect to total mercury (THg) and methylmercury (MeHg) concentrations as well as Bacterial and Archaeal community composition. Total Hg (0.25-990 μg g-1) and MeHg (1.3-30.5 ng g-1) varied between samples. Pyrosequencing (454 FLX) of the hypervariable v1-v3 regions of the 16S rRNA genes showed that Proteobacteria, Actinobacteria, Chloroflexi, Acidobacteria, Euryarchaeota, and Crenarchaeota were dominant in all samples. The Bacterial α-diversity was higher in samples with relatively Low THg and MeHg and decreased with increasing THg and MeHg concentrations. In contrast, Archaeal α-diversity increased with increasing of MeHg concentrations but did not correlate with changes in THg concentrations. Overall, the methylation gene hgcAB copy number increased with both increasing THg and MeHg concentrations. The microbial communities at High THg and High MeHg appear to be adapted by species that are both Hg resistant and carry hgcAB genes for MeHg production. The relatively high abundance of both sulfate-reducing δ-Proteobacteria and methanogenic Archaea, as well as their positive correlations with increasing THg and MeHg concentrations, suggests that these microorganisms are the primary Hg-methylators in the rice paddy soils in Guizhou, China.

Robust Mercury Methylation across Diverse Methanogenic Archaea

Methylmercury (MeHg) production was compared among nine cultured methanogenic archaea that contain hgcAB, a gene pair that codes for mercury (Hg) methylation. The methanogens tested produced MeHg at inherently different rates, even when normalized to growth rate and Hg availability. Eight of the nine tested were capable of MeHg production greater than that of spent- and uninoculated-medium controls during batch culture growth. Methanococcoides methylutens, an hgcAB+ strain with a fused gene pair, was unable to produce more MeHg than controls. Maximal conversion of Hg to MeHg through a full batch culture growth cycle for each species (except M. methylutens) ranged from 2 to >50% of the added Hg(II) or between 0.2 and 17 pmol of MeHg/mg of protein. Three of the species produced >10% MeHg. The ability to produce MeHg was confirmed in several hgcAB+ methanogens that had not previously been tested (Methanocella paludicola SANAE, Methanocorpusculum bavaricum, Methanofollis liminatans GKZPZ, and Methanosphaerula palustris E1-9c). Maximal methylation was observed at low sulfide concentrations (<100 μM) and in the presence of 0.5 to 5 mM cysteine. For M. hollandica, the addition of up to 5 mM cysteine enhanced MeHg production and cell growth in a concentration-dependent manner. As observed for bacterial Hg methylators, sulfide inhibited MeHg production. An initial evaluation of sulfide and thiol impacts on bioavailability showed methanogens responding to Hg complexation in the same way as do Deltaproteobacteria The mercury methylation rates of several methanogens rival those of the better-studied Hg-methylating sulfate- and iron-reducing DeltaproteobacteriaIMPORTANCEArchaea, specifically methanogenic organisms, play a role in mercury methylation in nature, but their global importance to MeHg production and the subsequent risk to ecosystems are not known. Methanogenesis has been linked to Hg methylation in several natural habitats where methylmercury production incurs risk to people and ecosystems, including rice paddies and permafrost. In this study, we confirm that most methanogens carrying the hgcAB gene pair are capable of Hg methylation. We found that methylation rates vary inherently among hgcAB+ methanogens but that several species are capable of MeHg production at rates that rival those of the better-know Hg-methylating sulfate- and iron-reducing bacteria. Methanogens may need to be considered equally with sulfate and iron reducers in evaluations of MeHg production in nature.

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

Standards for sequencing the microbial 'uncultivated majority', namely bacterial and archaeal single-cell genome sequences, and genome sequences from metagenomic datasets, are proposed.

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.

Hypersaline sapropels act as hotspots for microbial dark matter

Present-day terrestrial analogue sites are crucial ground truth proxies for studying life in geochemical conditions close to those assumed to be present on early Earth or inferred to exist on other celestial bodies (e.g. Mars, Europa). Although hypersaline sapropels are border-of-life habitats with moderate occurrence, their microbiological and physicochemical characterization lags behind. Here, we study the diversity of life under low water activity by describing the prokaryotic communities from two disparate hypersaline sapropels (Transylvanian Basin, Romania) in relation to geochemical milieu and pore water chemistry, while inferring their role in carbon cycling by matching taxa to known taxon-specific biogeochemical functions. The polyphasic approach combined deep coverage SSU rRNA gene amplicon sequencing and bioinformatics with RT-qPCR and physicochemical investigations. We found that sapropels developed an analogous elemental milieu and harbored prokaryotes affiliated with fifty-nine phyla, among which the most abundant were Proteobacteria, Bacteroidetes and Chloroflexi. Containing thirty-two candidate divisions and possibly undocumented prokaryotic lineages, the hypersaline sapropels were found to accommodate one of the most diverse and novel ecosystems reported to date and may contribute to completing the phylogenetic branching of the tree of life.

A Complex Endomembrane System in the Archaeon Ignicoccus hospitalis Tapped by Nanoarchaeum equitans

Based on serial sectioning, focused ion beam scanning electron microscopy (FIB/SEM), and electron tomography, we depict in detail the highly unusual anatomy of the marine hyperthermophilic crenarchaeon, Ignicoccus hospitalis. Our data support a complex and dynamic endomembrane system consisting of cytoplasmic protrusions, and with secretory function. Moreover, we reveal that the cytoplasm of the putative archaeal ectoparasite Nanoarchaeum equitans can get in direct contact with this endomembrane system, complementing and explaining recent proteomic, transcriptomic and metabolomic data on this inter-archaeal relationship. In addition, we identified a matrix of filamentous structures and/or tethers in the voluminous inter-membrane compartment (IMC) of I. hospitalis, which might be responsible for membrane dynamics. Overall, this unusual cellular compartmentalization, ultrastructure and dynamics in an archaeon that belongs to the recently proposed TACK superphylum prompts speculation that the eukaryotic endomembrane system might originate from Archaea.

Multi-omics analysis provides insight to the Ignicoccus hospitalis-Nanoarchaeum equitans association

BACKGROUND

Studies of interspecies interactions are inherently difficult due to the complex mechanisms which enable these relationships. A model system for studying interspecies interactions is the marine hyperthermophiles Ignicoccus hospitalis and Nanoarchaeum equitans. Recent independently-conducted 'omics' analyses have generated insights into the molecular factors modulating this association. However, significant questions remain about the nature of the interactions between these archaea.

METHODS

We jointly analyzed multiple levels of omics datasets obtained from published, independent transcriptomics, proteomics, and metabolomics analyses. DAVID identified functionally-related groups enriched when I. hospitalis is grown alone or in co-culture with N. equitans. Enriched molecular pathways were subsequently visualized using interaction maps generated using STRING.

RESULTS

Key findings of our multi-level omics analysis indicated that I. hospitalis provides precursors to N. equitans for energy metabolism. Analysis indicated an overall reduction in diversity of metabolic precursors in the I. hospitalis-N. equitans co-culture, which has been connected to the differential use of ribosomal subunits and was previously unnoticed. We also identified differences in precursors linked to amino acid metabolism, NADH metabolism, and carbon fixation, providing new insights into the metabolic adaptions of I. hospitalis enabling the growth of N. equitans.

CONCLUSIONS

This multi-omics analysis builds upon previously identified cellular patterns while offering new insights into mechanisms that enable the I. hospitalis-N. equitans association.

GENERAL SIGNIFICANCE

Our study applies statistical and visualization techniques to a mixed-source omics dataset to yield a more global insight into a complex system, that was not readily discernable from separate omics studies.

Temporal Dynamics of In-Field Bioreactor Populations Reflect the Groundwater System and Respond Predictably to Perturbation

Temporal variability complicates testing the influences of environmental variability on microbial community structure and thus function. An in-field bioreactor system was developed to assess oxic versus anoxic manipulations on in situ groundwater communities. Each sample was sequenced (16S SSU rRNA genes, average 10,000 reads), and biogeochemical parameters are monitored by quantifying 53 metals, 12 organic acids, 14 anions, and 3 sugars. Changes in dissolved oxygen (DO), pH, and other variables were similar across bioreactors. Sequencing revealed a complex community that fluctuated in-step with the groundwater community and responded to DO. This also directly influenced the pH, and so the biotic impacts of DO and pH shifts are correlated. A null model demonstrated that bioreactor communities were driven in part not only by experimental conditions but also by stochastic variability and did not accurately capture alterations in diversity during perturbations. We identified two groups of abundant OTUs important to this system; one was abundant in high DO and pH and contained heterotrophs and oxidizers of iron, nitrite, and ammonium, whereas the other was abundant in low DO with the capability to reduce nitrate. In-field bioreactors are a powerful tool for capturing natural microbial community responses to alterations in geochemical factors beyond the bulk phase.

Enrichment of Root Endophytic Bacteria from Populus deltoides and Single-Cell-Genomics Analysis

Bacterial endophytes that colonize Populus trees contribute to nutrient acquisition, prime immunity responses, and directly or indirectly increase both above- and below-ground biomasses. Endophytes are embedded within plant material, so physical separation and isolation are difficult tasks. Application of culture-independent methods, such as metagenome or bacterial transcriptome sequencing, has been limited due to the predominance of DNA from the plant biomass. Here, we describe a modified differential and density gradient centrifugation-based protocol for the separation of endophytic bacteria from Populus roots. This protocol achieved substantial reduction in contaminating plant DNA, allowed enrichment of endophytic bacteria away from the plant material, and enabled single-cell genomics analysis. Four single-cell genomes were selected for whole-genome amplification based on their rarity in the microbiome (potentially uncultured taxa) as well as their inferred abilities to form associations with plants. Bioinformatics analyses, including assembly, contamination removal, and completeness estimation, were performed to obtain single-amplified genomes (SAGs) of organisms from the phyla Armatimonadetes, Verrucomicrobia, and Planctomycetes, which were unrepresented in our previous cultivation efforts. Comparative genomic analysis revealed unique characteristics of each SAG that could facilitate future cultivation efforts for these bacteria.

IMPORTANCE

Plant roots harbor a diverse collection of microbes that live within host tissues. To gain a comprehensive understanding of microbial adaptations to this endophytic lifestyle from strains that cannot be cultivated, it is necessary to separate bacterial cells from the predominance of plant tissue. This study provides a valuable approach for the separation and isolation of endophytic bacteria from plant root tissue. Isolated live bacteria provide material for microbiome sequencing, single-cell genomics, and analyses of genomes of uncultured bacteria to provide genomics information that will facilitate future cultivation attempts.

Culture Independent Genomic Comparisons Reveal Environmental Adaptations for Altiarchaeales

The recently proposed candidatus order Altiarchaeales remains an uncultured archaeal lineage composed of genetically diverse, globally widespread organisms frequently observed in anoxic subsurface environments. In spite of 15 years of studies on the psychrophilic biofilm-producing Candidatus Altiarchaeum hamiconexum and its close relatives, very little is known about the phylogenetic and functional diversity of the widespread free-living marine members of this taxon. From methanogenic sediments in the White Oak River Estuary, NC, USA, we sequenced a single cell amplified genome (SAG), WOR_SM1_SCG, and used it to identify and refine two high-quality genomes from metagenomes, WOR_SM1_79 and WOR_SM1_86-2, from the same site. These three genomic reconstructions form a monophyletic group, which also includes three previously published genomes from metagenomes from terrestrial springs and a SAG from Sakinaw Lake in a group previously designated as pMC2A384. A synapomorphic mutation in the Altiarchaeales tRNA synthetase β subunit, pheT, caused the protein to be encoded as two subunits at non-adjacent loci. Consistent with the terrestrial spring clades, our estuarine genomes contained a near-complete autotrophic metabolism, H2 or CO as potential electron donors, a reductive acetyl-CoA pathway for carbon fixation, and methylotroph-like NADP(H)-dependent dehydrogenase. Phylogenies based on 16S rRNA genes and concatenated conserved proteins identified two distinct sub-clades of Altiarchaeales, Alti-1 populated by organisms from actively flowing springs, and Alti-2 which was more widespread, diverse, and not associated with visible mats. The core Alti-1 genome suggested Alti-1 is adapted for the stream environment with lipopolysaccharide production capacity and extracellular hami structures. The core Alti-2 genome suggested members of this clade are free-living with distinct mechanisms for energy maintenance, motility, osmoregulation, and sulfur redox reactions. These data suggested that the hamus structures found in Candidatus Altiarchaeum hamiconexum are not present outside of stream-adapted Altiarchaeales. Homologs to a Na(+) transporter and membrane bound coenzyme A disulfide reductase that were unique to the brackish sediment Alti-2 genomes, could indicate adaptations to the estuarine, sulfur-rich environment.

Genomics-informed isolation and characterization of a symbiotic Nanoarchaeota system from a terrestrial geothermal environment

Biological features can be inferred, based on genomic data, for many microbial lineages that remain uncultured. However, cultivation is important for characterizing an organism's physiology and testing its genome-encoded potential. Here we use single-cell genomics to infer cultivation conditions for the isolation of an ectosymbiotic Nanoarchaeota (‘Nanopusillus acidilobi') and its host (Acidilobus, a crenarchaeote) from a terrestrial geothermal environment. The cells of ‘Nanopusillus' are among the smallest known cellular organisms (100–300 nm). They appear to have a complete genetic information processing machinery, but lack almost all primary biosynthetic functions as well as respiration and ATP synthesis. Genomic and proteomic comparison with its distant relative, the marine Nanoarchaeum equitans illustrate an ancient, common evolutionary history of adaptation of the Nanoarchaeota to ectosymbiosis, so far unique among the Archaea.

Many microbial lineages have not yet been cultured, which hampers our understanding of their physiology. Here, Wurch et al. use single-cell genomics to infer cultivation conditions for the isolation of a tiny ectosymbiotic nanoarchaeon and its crenarchaeota host from a geothermal spring.

A comprehensive benchmarking study of protocols and sequencing platforms for 16S rRNA community profiling

Background

In the last 5 years, the rapid pace of innovations and improvements in sequencing technologies has completely changed the landscape of metagenomic and metagenetic experiments. Therefore, it is critical to benchmark the various methodologies for interrogating the composition of microbial communities, so that we can assess their strengths and limitations. The most common phylogenetic marker for microbial community diversity studies is the 16S ribosomal RNA gene and in the last 10 years the field has moved from sequencing a small number of amplicons and samples to more complex studies where thousands of samples and multiple different gene regions are interrogated.

Results

We assembled 2 synthetic communities with an even (EM) and uneven (UM) distribution of archaeal and bacterial strains and species, as metagenomic control material, to assess performance of different experimental strategies. The 2 synthetic communities were used in this study, to highlight the limitations and the advantages of the leading sequencing platforms: MiSeq (Illumina), The Pacific Biosciences RSII, 454 GS-FLX/+ (Roche), and IonTorrent (Life Technologies). We describe an extensive survey based on synthetic communities using 3 experimental designs (fusion primers, universal tailed tag, ligated adaptors) across the 9 hypervariable 16S rDNA regions. We demonstrate that library preparation methodology can affect data interpretation due to different error and chimera rates generated during the procedure. The observed community composition was always biased, to a degree that depended on the platform, sequenced region and primer choice. However, crucially, our analysis suggests that 16S rRNA sequencing is still quantitative, in that relative changes in abundance of taxa between samples can be recovered, despite these biases.

Conclusion

We have assessed a range of experimental conditions across several next generation sequencing platforms using the most up-to-date configurations. We propose that the choice of sequencing platform and experimental design needs to be taken into consideration in the early stage of a project by running a small trial consisting of several hypervariable regions to quantify the discriminatory power of each region. We also suggest that the use of a synthetic community as a positive control would be beneficial to identify the potential biases and procedural drawbacks that may lead to data misinterpretation. The results of this study will serve as a guideline for making decisions on which experimental condition and sequencing platform to consider to achieve the best microbial profiling.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2194-9) contains supplementary material, which is available to authorized users.

Global prevalence and distribution of genes and microorganisms involved in mercury methylation

Mercury (Hg) methylation produces the neurotoxic, highly bioaccumulative methylmercury (MeHg). The highly conserved nature of the recently identified Hg methylation genes hgcAB provides a foundation for broadly evaluating spatial and niche-specific patterns of microbial Hg methylation potential in nature. We queried hgcAB diversity and distribution in >3500 publicly available microbial metagenomes, encompassing a broad range of environments and generating a new global view of Hg methylation potential. The hgcAB genes were found in nearly all anaerobic (but not aerobic) environments, including oxygenated layers of the open ocean. Critically, hgcAB was effectively absent in ~1500 human and mammalian microbiomes, suggesting a low risk of endogenous MeHg production. New potential methylation habitats were identified, including invertebrate digestive tracts, thawing permafrost soils, coastal "dead zones," soils, sediments, and extreme environments, suggesting multiple routes for MeHg entry into food webs. Several new taxonomic groups capable of methylating Hg emerged, including lineages having no cultured representatives. Phylogenetic analysis points to an evolutionary relationship between hgcA and genes encoding corrinoid iron-sulfur proteins functioning in the ancient Wood-Ljungdahl carbon fixation pathway, suggesting that methanogenic Archaea may have been the first to perform these biotransformations.

Untargeted metabolomics studies employing NMR and LC-MS reveal metabolic coupling between Nanoarcheum equitans and its archaeal host Ignicoccus hospitalis

Interspecies interactions are the basis of microbial community formation and infectious diseases. Systems biology enables the construction of complex models describing such interactions, leading to a better understanding of disease states and communities. However, before interactions between complex organisms can be understood, metabolic and energetic implications of simpler real-world host-microbe systems must be worked out. To this effect, untargeted metabolomics experiments were conducted and integrated with proteomics data to characterize key molecular-level interactions between two hyperthermophilic microbial species, both of which have reduced genomes. Metabolic changes and transfer of metabolites between the archaea Ignicoccus hospitalis and Nanoarcheum equitans were investigated using integrated LC-MS and NMR metabolomics. The study of such a system is challenging, as no genetic tools are available, growth in the laboratory is challenging, and mechanisms by which they interact are unknown. Together with information about relative enzyme levels obtained from shotgun proteomics, the metabolomics data provided useful insights into metabolic pathways and cellular networks of I. hospitalis that are impacted by the presence of N. equitans, including arginine, isoleucine, and CTP biosynthesis. On the organismal level, the data indicate that N. equitans exploits metabolites generated by I. hospitalis to satisfy its own metabolic needs. This finding is based on N. equitans's consumption of a significant fraction of the metabolite pool in I. hospitalis that cannot solely be attributed to increased biomass production for N. equitans. Combining LC-MS and NMR metabolomics datasets improved coverage of the metabolome and enhanced the identification and quantitation of cellular metabolites.

Contrasting taxonomic stratification of microbial communities in two hypersaline meromictic lakes

Hypersaline meromictic lakes are extreme environments in which water stratification is associated with powerful physicochemical gradients and high salt concentrations. Furthermore, their physical stability coupled with vertical water column partitioning makes them important research model systems in microbial niche differentiation and biogeochemical cycling. Here, we compare the prokaryotic assemblages from Ursu and Fara Fund hypersaline meromictic lakes (Transylvanian Basin, Romania) in relation to their limnological factors and infer their role in elemental cycling by matching taxa to known taxon-specific biogeochemical functions. To assess the composition and structure of prokaryotic communities and the environmental factors that structure them, deep-coverage small subunit (SSU) ribosomal RNA (rDNA) amplicon sequencing, community domain-specific quantitative PCR and physicochemical analyses were performed on samples collected along depth profiles. The analyses showed that the lakes harbored multiple and diverse prokaryotic communities whose distribution mirrored the water stratification patterns. Ursu Lake was found to be dominated by Bacteria and to have a greater prokaryotic diversity than Fara Fund Lake that harbored an increased cell density and was populated mostly by Archaea within oxic strata. In spite of their contrasting diversity, the microbial populations indigenous to each lake pointed to similar physiological functions within carbon degradation and sulfate reduction. Furthermore, the taxonomy results coupled with methane detection and its stable C isotope composition indicated the presence of a yet-undescribed methanogenic group in the lakes' hypersaline monimolimnion. In addition, ultrasmall uncultivated archaeal lineages were detected in the chemocline of Fara Fund Lake, where the recently proposed Nanohaloarchaeota phylum was found to thrive.

Structure, mineralogy, and microbial diversity of geothermal spring microbialites associated with a deep oil drilling in Romania

Modern mineral deposits play an important role in evolutionary studies by providing clues to the formation of ancient lithified microbial communities. Here we report the presence of microbialite-forming microbial mats in different microenvironments at 32°C, 49°C, and 65°C around the geothermal spring from an abandoned oil drill in Ciocaia, Romania. The mineralogy and the macro- and microstructure of the microbialites were investigated, together with their microbial diversity based on a 16S rRNA gene amplicon sequencing approach. The calcium carbonate is deposited mainly in the form of calcite. At 32°C and 49°C, the microbialites show a laminated structure with visible microbial mat-carbonate crystal interactions. At 65°C, the mineral deposit is clotted, without obvious organic residues. Partial 16S rRNA gene amplicon sequencing showed that the relative abundance of the phylum Archaea was low at 32°C (<0.5%) but increased significantly at 65°C (36%). The bacterial diversity was either similar to other microbialites described in literature (the 32°C sample) or displayed a specific combination of phyla and classes (the 49°C and 65°C samples). Bacterial taxa were distributed among 39 phyla, out of which 14 had inferred abundances >1%. The dominant bacterial groups at 32°C were Cyanobacteria, Gammaproteobacteria, Firmicutes, Bacteroidetes, Chloroflexi, Thermi, Actinobacteria, Planctomycetes, and Defferibacteres. At 49°C, there was a striking dominance of the Gammaproteobacteria, followed by Firmicutes, Bacteroidetes, and Armantimonadetes. The 65°C sample was dominated by Betaproteobacteria, Firmicutes, [OP1], Defferibacteres, Thermi, Thermotogae, [EM3], and Nitrospirae. Several groups from Proteobacteria and Firmicutes, together with Halobacteria and Melainabacteria were described for the first time in calcium carbonate deposits. Overall, the spring from Ciocaia emerges as a valuable site to probe microbes-minerals interrelationships along thermal and geochemical gradients.

Community analysis of plant biomass-degrading microorganisms from Obsidian Pool, Yellowstone National Park

The conversion of lignocellulosic biomass into biofuels can potentially be improved by employing robust microorganisms and enzymes that efficiently deconstruct plant polysaccharides at elevated temperatures. Many of the geothermal features of Yellowstone National Park (YNP) are surrounded by vegetation providing a source of allochthonic material to support heterotrophic microbial communities adapted to utilize plant biomass as a primary carbon and energy source. In this study, a well-known hot spring environment, Obsidian Pool (OBP), was examined for potential biomass-active microorganisms using cultivation-independent and enrichment techniques. Analysis of 33,684 archaeal and 43,784 bacterial quality-filtered 16S rRNA gene pyrosequences revealed that archaeal diversity in the main pool was higher than bacterial; however, in the vegetated area, overall bacterial diversity was significantly higher. Of notable interest was a flooded depression adjacent to OBP supporting a stand of Juncus tweedyi, a heat-tolerant rush commonly found growing near geothermal features in YNP. The microbial community from heated sediments surrounding the plants was enriched in members of the Firmicutes including potentially (hemi)cellulolytic bacteria from the genera Clostridium, Anaerobacter, Caloramator, Caldicellulosiruptor, and Thermoanaerobacter. Enrichment cultures containing model and real biomass substrates were established at a wide range of temperatures (55-85 °C). Microbial activity was observed up to 80 °C on all substrates including Avicel, xylan, switchgrass, and Populus sp. Independent of substrate, Caloramator was enriched at lower (<65 °C) temperatures while highly active cellulolytic bacteria Caldicellulosiruptor were dominant at high (>65 °C) temperatures.