EMIRGE: reconstruction of full-length ribosomal genes from microbial community short read sequencing data

Inoculum source determines the stress resistance of electroactive functional taxa in biofilms: A metagenomic perspective

The inoculum has a crucial impact on bioreactor initialization and performance. However, there is currently a lack of guidance on selecting appropriate inocula for applications in environmental biotechnology. In this study, we applied microbial electrolysis cells (MECs) as models to investigate the differences in the functional potential of electroactive microorganisms (EAMs) within anodic biofilms developed from four different inocula (natural or artificial), using shotgun metagenomic techniques. We specifically focused on extracellular electron transfer (EET) function and stress resistance, which affect the performance and stability of MECs. Community profiling revealed that the family Geobacteraceae was the key EAM taxon in all biofilms, with Geobacter as the dominant genus. The c-type cytochrome gene imcH showed universal importance for Geobacteraceae EET and was utilized as a marker gene to evaluate the EET potential of EAMs. Additionally, stress response functional genes were used to assess the stress resistance potential of Geobacter species. Comparative analysis of imcH gene abundance revealed that EAMs with comparable overall EET potential could be enriched from artificial and natural inocula (P > 0.05). However, quantification of stress response gene copy numbers in the genomes demonstrated that EAMs originating from natural inocula possessed superior stress resistance potential (196 vs. 163). Overall, this study provides novel perspectives on the inoculum effect in bioreactors and offers theoretical guidance for selecting inoculum in environmental engineering applications.

Unique episymbiotic relationship between Candidatus Patescibacteria and Zoogloea in activated sludge flocs at a municipal wastewater treatment plant

Candidatus Patescibacteria, also known as candidate phyla radiation (CPR), including the class-level uncultured clade JAEDAM01 (formerly a subclass of Gracilibacteria/GN02/BD1-5), are ubiquitous in activated sludge. However, their characteristics and relationships with other organisms are largely unknown. They are believed to be episymbiotic, endosymbiotic or predatory. Despite our understanding of their limited metabolic capacity, their precise roles remain elusive due to the difficulty in cultivating and identifying them. In previous research, we successfully recovered high-quality metagenome-assembled genomes (MAGs), including a member of JAEDAM01 from activated sludge flocs. In this study, we designed new probes to visualize the targeted JAEDAM01-associated MAG HHAS10 and identified its host using fluorescence in situ hybridization (FISH). The FISH observations revealed that JAEDAM01 HHAS10-like cells were located within dense clusters of Zoogloea, and the fluorescence brightness of zoogloeal cells decreased in the vicinity of the CPR cells. The Zoogloea MAGs possessed genes related to extracellular polymeric substance biosynthesis, floc formation and nutrient removal, including a polyhydroxyalkanoate (PHA) accumulation pathway. The JAEDAM01 MAG HHAS10 possessed genes associated with type IV pili, competence protein EC and PHA degradation, suggesting a Zoogloea-dependent lifestyle in activated sludge flocs. These findings indicate a new symbiotic relationship between JAEDAM01 and Zoogloea.

METASEED: a novel approach to full-length 16S rRNA gene reconstruction from short read data

BACKGROUND

With the emergence of Oxford Nanopore technology, now the on-site sequencing of 16S rRNA from environments is available. Due to the error level and structure, the analysis of such data demands some database of reference sequences. However, many taxa from complex and diverse environments, have poor representation in publicly available databases. In this paper, we propose the METASEED pipeline for the reconstruction of full-length 16S sequences from such environments, in order to improve the reference for the subsequent use of on-site sequencing.

RESULTS

We show that combining high-precision short-read sequencing of both 16S and full metagenome from the same samples allow us to reconstruct high-quality 16S sequences from the more abundant taxa. A significant novelty is the carefully designed collection of metagenome reads that matches the 16S amplicons, based on a combination of uniqueness and abundance. Compared to alternative approaches this produces superior results.

CONCLUSION

Our pipeline will facilitate numerous studies associated with various unknown microorganisms, thus allowing the comprehension of the diverse environments. The pipeline is a potential tool in generating a full length 16S rRNA gene database for any environment.

The impact of antibiotics on the gut microbiota of children recovering from watery diarrhoea

Infectious diarrhoeal diseases remain a substantial health burden in young children in low- and middle-income countries. The disease and its variable treatment options significantly alter the gut microbiome, which may affect clinical outcomes and overall gut health. Antibiotics are often prescribed, but their impact on the gut microbiome during recovery is unclear. Here, we used 16S rRNA sequencing to investigate changes in the gut microbiota in Vietnamese children with acute watery diarrhoea, and highlight the impact of antibiotic treatment on these changes. Our analyses identified that, regardless of treatment, recovery was characterised by reductions in Streptococcus and Rothia species and expansion of Bacteroides/Phocaeicola, Lachnospiraceae and Ruminococcacae taxa. Antibiotic treatment significantly delayed the temporal increases in alpha- and beta-diversity within patients, resulting in distinctive patterns of taxonomic change. These changes included a pronounced, transient overabundance of Enterococcus species and depletion of Bifidobacterium pseudocatenulatum. Our findings demonstrate that antibiotic treatment slows gut microbiota recovery in children following watery diarrhoea.

Time-shifted expression of acetoclastic and methylotrophic methanogenesis by a single Methanosarcina genomospecies predominates the methanogen dynamics in Philippine rice field soil

BACKGROUND

The final step in the anaerobic decomposition of biopolymers is methanogenesis. Rice field soils are a major anthropogenic source of methane, with straw commonly used as a fertilizer in rice farming. Here, we aimed to decipher the structural and functional responses of the methanogenic community to rice straw addition during an extended anoxic incubation (120 days) of Philippine paddy soil. The research combined process measurements, quantitative real-time PCR and RT-PCR of particular biomarkers (16S rRNA, mcrA), and meta-omics (environmental genomics and transcriptomics).

RESULTS

The analysis methods collectively revealed two major bacterial and methanogenic activity phases: early (days 7 to 21) and late (days 28 to 60) community responses, separated by a significant transient decline in microbial gene and transcript abundances and CH4 production rate. The two methanogenic activity phases corresponded to the greatest rRNA and mRNA abundances of the Methanosarcinaceae but differed in the methanogenic pathways expressed. While three genetically distinct Methanosarcina populations contributed to acetoclastic methanogenesis during the early activity phase, the late activity phase was defined by methylotrophic methanogenesis performed by a single Methanosarcina genomospecies. Closely related to Methanosarcina sp. MSH10X1, mapping of environmental transcripts onto metagenome-assembled genomes (MAGs) and population-specific reference genomes revealed this genomospecies as the key player in acetoclastic and methylotrophic methanogenesis. The anaerobic food web was driven by a complex bacterial community, with Geobacteraceae and Peptococcaceae being putative candidates for a functional interplay with Methanosarcina. Members of the Methanocellaceae were the key players in hydrogenotrophic methanogenesis, while the acetoclastic activity of Methanotrichaceae members was detectable only during the very late community response.

CONCLUSIONS

The predominant but time-shifted expression of acetoclastic and methylotrophic methanogenesis by a single Methanosarcina genomospecies represents a novel finding that expands our hitherto knowledge of the methanogenic pathways being highly expressed in paddy soils. Video Abstract.

Polyploidy is widespread in Microsporidia

Microsporidia are obligate intracellular eukaryotic parasites with an extremely broad host range. They have both economic and public health importance. Ploidy in microsporidia is variable, with a few species formally identified as diploid and one as polyploid. Given the increase in the number of studies sequencing microsporidian genomes, it is now possible to assess ploidy levels across all currently explored microsporidian diversity. We estimate ploidy for all microsporidian data sets available on the Sequence Read Archive using k-mer-based analyses, indicating that polyploidy is widespread in Microsporidia and that ploidy change is dynamic in the group. Using genome-wide heterozygosity estimates, we also show that polyploid microsporidian genomes are relatively homozygous, and we discuss the implications of these findings on the timing of polyploidization events and their origin.IMPORTANCEMicrosporidia are single-celled intracellular parasites, distantly related to fungi, that can infect a broad range of hosts, from humans all the way to protozoans. Exploiting the wealth of microsporidian genomic data available, we use k-mer-based analyses to assess ploidy status across the group. Understanding a genome's ploidy is crucial in order to assemble it effectively and may also be relevant for better understanding a parasite's behavior and life cycle. We show that tetraploidy is present in at least six species in Microsporidia and that these polyploidization events are likely to have occurred independently. We discuss why these findings may be paradoxical, given that Microsporidia, like other intracellular parasites, have extremely small, reduced genomes.

Comparative genomic analysis of Microcystis strain diversity using conserved marker genes

Microcystis-dominated cyanobacterial harmful algal blooms (cyanoHABs) have a global impact on freshwater environments, affecting both wildlife and human health. Microcystis diversity and function in field samples and laboratory cultures can be determined by sequencing whole genomes of cultured isolates or natural populations, but these methods remain computationally and financially expensive. Amplicon sequencing of marker genes is a lower cost and higher throughput alternative to characterize strain composition and diversity in mixed samples. However, the selection of appropriate marker gene region(s) and primers requires prior understanding of the relationship between single gene genotype, whole genome content, and phenotype. To identify phylogenetic markers of Microcystis strain diversity, we compared phylogenetic trees built from each of 2,351 individual core genes to an established phylogeny and assessed the ability of these core genes to predict whole genome content and bioactive compound genotypes. We identified single-copy core genes better able to resolve Microcystis phylogenies than previously identified marker genes. We developed primers suitable for current Illumina-based amplicon sequencing with near-complete coverage of available Microcystis genomes and demonstrate that they outperform existing options for assessing Microcystis strain composition. Results showed that genetic markers can be used to infer Microcystis gene content and phenotypes such as potential production of bioactive compounds , although marker performance varies by bioactive compound gene and sequence similarity. Finally, we demonstrate that these markers can be used to characterize the Microcystis strain composition of laboratory or field samples like those collected for surveillance and modeling of Microcystis-dominated cyanobacterial harmful algal blooms.

MRGM: an enhanced catalog of mouse gut microbial genomes substantially broadening taxonomic and functional landscapes

Mouse gut microbiome research is pivotal for understanding the human gut microbiome, providing insights into disease modeling, host-microbe interactions, and the dietary influence on the gut microbiome. To enhance the translational value of mouse gut microbiome studies, we need detailed and high-quality catalogs of mouse gut microbial genomes. We introduce the Mouse Reference Gut Microbiome (MRGM), a comprehensive catalog with 42,245 non-redundant mouse gut bacterial genomes across 1,524 species. MRGM marks a 40% increase in the known taxonomic diversity of mouse gut microbes, capturing previously underrepresented lineages through refined genome quality assessment techniques. MRGM not only broadens the taxonomic landscape but also enriches the functional landscape of the mouse gut microbiome. Using deep learning, we have elevated the Gene Ontology annotation rate for mouse gut microbial proteins from 3.2% with orthology to 60%, marking an over 18-fold increase. MRGM supports both DNA- and marker-based taxonomic profiling by providing custom databases, surpassing previous catalogs in performance. Finally, taxonomic and functional comparisons between human and mouse gut microbiota reveal diet-driven divergences in their taxonomic composition and functional enrichment. Overall, our study highlights the value of high-quality microbial genome catalogs in advancing our understanding of the co-evolution between gut microbes and their host.

Phylogenetic diversity of functional genes in deep-sea cold seeps: a novel perspective on metagenomics

BACKGROUND

Leakages of cold, methane-rich fluids from subsurface reservoirs to the sea floor are termed cold seeps. Recent exploration of the deep sea has shed new light on the microbial communities in cold seeps. However, conventional metagenomic methods largely rely on reference databases and neglect the phylogeny of functional genes.

RESULTS

In this study, we developed the REMIRGE program to retrieve the full-length functional genes from shotgun metagenomic reads and fully explored the phylogenetic diversity in cold seep sediments. The abundance and diversity of functional genes involved in the methane, sulfur, and nitrogen cycles differed in the non-seep site and five cold seep sites. In one Haima cold seep site, the divergence of functional groups was observed at the centimeter scale of sediment depths, with the surface layer potentially acting as a reservoir of microbial species and functions. Additionally, positive correlations were found between specific gene sequence clusters of relevant genes, indicating coupling occurred within specific functional groups.

CONCLUSION

REMIRGE revealed divergent phylogenetic diversity of functional groups and functional pathway preferences in a deep-sea cold seep at finer scales, which could not be detected by conventional methods. Our work highlights that phylogenetic information is conducive to more comprehensive functional profiles, and REMIRGE has the potential to uncover more new insights from shotgun metagenomic data. Video Abstract.

Metagenomic analysis reveals indole signaling effect on microbial community in sequencing batch reactors: Quorum sensing inhibition and antibiotic resistance enrichment

Indole is an essential signal molecule in microbial studies. However, its ecological role in biological wastewater treatments remains enigmatic. This study explores the links between indole and complex microbial communities using sequencing batch reactors exposed to 0, 15, and 150 mg/L indole concentrations. A concentration of 150 mg/L indole enriched indole degrader Burkholderiales, while pathogens, such as Giardia, Plasmodium, and Besnoitia were inhibited at 15 mg/L indole concentration. At the same time, indole reduced the abundance of predicted genes in the "signaling transduction mechanisms" pathway via the Non-supervised Orthologous Groups distributions analysis. Indole significantly decreased the concentration of homoserine lactones, especially C₁₄-HSL. Furthermore, the quorum-sensing signaling acceptors containing LuxR, the dCACHE domain, and RpfC showed negative distributions with indole and indole oxygenase genes. Signaling acceptors' potential origins were mainly Burkholderiales, Actinobacteria, and Xanthomonadales. Meanwhile, concentrated indole (150 mg/L) increased the total abundance of antibiotic resistance genes by 3.52 folds, especially on aminoglycoside, multidrug, tetracycline, and sulfonamide. Based on Spearman's correlation analysis, the homoserine lactone degradation genes which were significantly impacted by indole negatively correlated with the antibiotic resistance gene abundance. This study brings new insights into the effect of indole signaling on in biological wastewater treatment plants.

Improved Assembly of Metagenome-Assembled Genomes and Viruses in Tibetan Saline Lake Sediment by HiFi Metagenomic Sequencing

With the development and reduced costs of high-throughput sequencing technology, environmental dark matter, such as novel metagenome-assembled genomes (MAGs) and viruses, is now being discovered easily. However, due to read length limitations, MAGs and viromes often suffer from genome discontinuity and deficiencies in key functional elements. Here, by applying long-read sequencing technology to sediment samples from a Tibetan saline lake, we comprehensively analyzed the performance of high-fidelity (HiFi) reads and the possibility of integration with short-read next-generation sequencing (NGS) data. In total, 207 full-length nonredundant 16S rRNA gene sequences and 19 full-length nonredundant 18S rRNA genes were directly obtained from HiFi reads, which greatly surpassed the retrieval performance of NGS technology. We carried out a cross-sectional comparison among multiple assembly strategies, referred to as 'NGS', 'Hybrid (NGS+HiFi)', and 'HiFi'. Two MAGs and 29 viruses with circular genomes were reconstructed using HiFi reads alone, indicating the great power of the 'HiFi' approach to assemble high-quality microbial genomes. Among the 3 strategies, the 'Hybrid' approach produced the highest number of medium/high-quality MAGs and viral genomes, while the ratio of MAGs containing 16S rRNA genes was significantly improved in the 'HiFi' assembly results. Overall, our study provides a practical metagenomic resolution for analyzing complex environmental samples by taking advantage of both the short-read and HiFi long-read sequencing methods to extract the maximum amount of information, including data on prokaryotes, eukaryotes, and viruses, via the 'Hybrid' approach. IMPORTANCE To expand the understanding of microbial dark matter in the environment, we did the first comparative evaluation of multiple assembly strategies based on high-throughput short-read and HiFi data from lake sediments metagenomic sequencing. The results demonstrated great improvement of the 'Hybrid' assembly method (short-read next-generation sequencing data plus HiFi data) in the recovery of medium/high-quality MAGs and viral genomes. Further analysis showed that HiFi data is important to retrieve the complete circular prokaryotic and viral genomes. Meanwhile, hundreds of full-length 16S/18S rRNA genes were assembled directly from HiFi data, which facilitated the species composition studies of complex environmental samples, especially for understanding micro-eukaryotes. Therefore, the application of the latest HiFi long-read sequencing could greatly improve the metagenomic assembly integrity and promote environmental microbiome research.

Profiling Microbial Communities in Idiopathic Granulomatous Mastitis

Idiopathic granulomatous mastitis (IGM) is a rare and benign inflammatory breast disease with ambiguous aetiology. Contrastingly, lactational mastitis (LM) is commonly diagnosed in breastfeeding women. To investigate IGM aetiology, we profiled the microbial flora of pus and skin in patients with IGM and LM. A total of 26 patients with IGM and 6 patients with LM were included in the study. The 16S rRNA sequencing libraries were constructed from 16S rRNA gene amplified from total DNA extracted from pus and skin swabs in patients with IGM and LM controls. Constructed libraries were multiplexed and paired-end sequenced on HiSeq4000. Metagenomic analysis was conducted using modified microbiome abundance analysis suite customised R-resource for paired pus and skin samples. Microbiome multivariable association analyses were performed using linear models. A total of 21 IGM and 3 LM paired pus and skin samples underwent metagenomic analysis. Bray−Curtis ecological dissimilarity distance showed dissimilarity across four sample types (IGM pus, IGM skin, LM pus, and LM skin; PERMANOVA, p < 0.001). No characteristic dominant genus was observed across the IGM samples. The IGM pus samples were more diverse than corresponding IGM skin samples (Shannon and Simpson index; Wilcoxon paired signed-rank tests, p = 0.022 and p = 0.07). Corynebacterium kroppenstedtii, reportedly associated with IGM in the literature, was higher in IGM pus samples than paired skin samples (Wilcoxon, p = 0.022). Three other species and nineteen genera were statistically significant in paired IGM pus−skin comparison after antibiotic treatment adjustment and multiple comparisons correction. Microbial profiles are unique between patients with IGM and LM. Inter-patient variability and polymicrobial IGM pus samples cannot implicate specific genus or species as an infectious cause for IGM.

Metagenome-assembled genome extraction and analysis from microbiomes using KBase

Uncultivated Bacteria and Archaea account for the vast majority of species on Earth, but obtaining their genomes directly from the environment, using shotgun sequencing, has only become possible recently. To realize the hope of capturing Earth's microbial genetic complement and to facilitate the investigation of the functional roles of specific lineages in a given ecosystem, technologies that accelerate the recovery of high-quality genomes are necessary. We present a series of analysis steps and data products for the extraction of high-quality metagenome-assembled genomes (MAGs) from microbiomes using the U.S. Department of Energy Systems Biology Knowledgebase (KBase) platform ( http://www.kbase.us/ ). Overall, these steps take about a day to obtain extracted genomes when starting from smaller environmental shotgun read libraries, or up to about a week from larger libraries. In KBase, the process is end-to-end, allowing a user to go from the initial sequencing reads all the way through to MAGs, which can then be analyzed with other KBase capabilities such as phylogenetic placement, functional assignment, metabolic modeling, pangenome functional profiling, RNA-Seq and others. While portions of such capabilities are available individually from other resources, the combination of the intuitive usability, data interoperability and integration of tools in a freely available computational resource makes KBase a powerful platform for obtaining MAGs from microbiomes. While this workflow offers tools for each of the key steps in the genome extraction process, it also provides a scaffold that can be easily extended with additional MAG recovery and analysis tools, via the KBase software development kit (SDK).

Predominance of Methanomicrobiales and diverse hydrocarbon-degrading taxa in the Appalachian coalbed biosphere revealed through metagenomics and genome-resolved metabolisms

Coalbed deposits are a unique subsurface environment and represent an underutilized resource for methane generation. Microbial communities extant in coalbed deposits are responsible for key subsurface biogeochemical cycling and could be utilized to enhance methane production in areas where existing gas wells have depleted methane stores, or in coalbeds that are unmined, or conversely be utilized for mitigation of methane release. Here we utilize metagenomics and metagenome-assembled genomes (MAGs) to identify extant microbial lineages and genome-resolved microbial metabolisms of coalbed produced water, which has not yet been explored in the Appalachian Basin (AppB). Our analyses resulted in the recovery of over 40 MAGs from 8 coalbed methane wells. The most commonly identified taxa among samples were hydrogenotrophic methanogens from the order Methanomicrobiales and these dominant MAGs were highly similar to one another. Conversely, low-abundance coalbed bacterial populations were taxonomically and functionally diverse, mostly belonging to a variety of Proteobacteria classes, and encoding various hydrocarbon solubilization and degradation pathways. The data presented herein provides novel insights into AppB coalbed microbial ecology, and our findings provide new perspectives on underrepresented Methanocalculus species and low-relative abundance bacterial assemblages in coalbed environments, and their potential roles in stimulation or mitigation of methane release.

Microbial carbon, sulfur, iron, and nitrogen cycling linked to the potential remediation of a meromictic acidic pit lake

Cueva de la Mora is a permanently stratified acidic pit lake and a model system for extreme acid mine drainage (AMD) studies. Using a combination of amplicon sequencing, metagenomics and metatranscriptomics we performed a taxonomically resolved analysis of microbial contributions to carbon, sulfur, iron, and nitrogen cycling. We found that active green alga Coccomyxa onubensis dominated the upper layer and chemocline. The chemocline had activity for iron(II) oxidation carried out by populations of Ca. Acidulodesulfobacterium, Ferrovum, Leptospirillium, and Armatimonadetes. Predicted activity for iron(III) reduction was only detected in the deep layer affiliated with Proteobacteria. Activity for dissimilatory nitrogen cycling including nitrogen fixation and nitrate reduction was primarily predicted in the chemocline. Heterotrophic archaeal populations with predicted activity for sulfide oxidation related to uncultured Thermoplasmatales dominated in the deep layer. Abundant sulfate-reducing Desulfomonile and Ca. Acidulodesulfobacterium populations were active in the chemocline. In the deep layer, uncultured populations from the bacterial phyla Actinobacteria, Chloroflexi, and Nitrospirae contributed to both sulfate reduction and sulfide oxidation. Based on this information we evaluated the potential for sulfide mineral precipitation in the deep layer as a tool for remediation. We argue that sulfide precipitation is not limited by microbial genetic potential but rather by the quantity and quality of organic carbon reaching the deep layer as well as by oxygen additions to the groundwater enabling sulfur oxidation. Addition of organic carbon and elemental sulfur should stimulate sulfate reduction and limit reoxidation of sulfide minerals.

Differences in gene expression patterns between cultured and natural Haloquadratum walsbyi ecotypes

Solar crystallizer ponds are characterized by high population density with a relatively simple community structure in terms of species composition. The microbial community in the solar saltern of Santa Pola (Alicante, Spain), is largely dominated by the hyperhalophilic square archaeon Haloquadratum walsbyi. Here we studied metatranscriptomes retrieved from a crystallizer pond during the winter of 2012 and summer of 2014 and compared Hqr. walsbyi's transcription patterns with that of the cultured strain Hqr. walsbyi HBSQ001. Significant differences were found between natural and the cultured grown strain in the distribution of transcript levels per gene. This likely reflects the adaptation of the cultured strain to the relative homogeneous growth conditions while the natural species, which is represented by multiple ecotypes, is adapted to heterogeneous environmental conditions and challenges of nutrient competition, viral attack, and other stressors. An important consequence of this study is that expression patterns obtained under artificial cultivation conditions cannot be directly extrapolated to gene expression under natural conditions. Moreover, we found 195 significantly differential expressed genes between the seasons, with 140 genes being higher expressed in winter and mainly encode proteins involved in energy and carbon source acquiring processes, and in stress responses.

Insights into the microbial life in silica-rich subterranean environments: microbial communities and ecological interactions in an orthoquartzite cave (Imawarì Yeuta, Auyan Tepui, Venezuela)

Microbial communities inhabiting caves in quartz-rich rocks are still underexplored, despite their possible role in the silica cycle. The world’s longest orthoquartzite cave, Imawarì Yeuta, represents a perfect arena for the investigation of the interactions between microorganisms and silica in non-thermal environments due to the presence of extraordinary amounts of amorphous silica speleothems of different kinds. In this work, the microbial diversity of Imawarì Yeuta was dissected by analyzing nineteen samples collected from different locations representative of different silica amorphization phases and types of samples. Specifically, we investigated the major ecological patterns in cave biodiversity, specific taxa enrichment, and the main ecological clusters through co-occurrence network analysis. Water content greatly contributed to the microbial communities’ composition and structures in the cave leading to the sample clustering into three groups DRY, WET, and WATER. Each of these groups was enriched in members of Actinobacteriota, Acidobacteriota, and Gammaproteobacteria, respectively. Alpha diversity analysis showed the highest value of diversity and richness for the WET samples, while the DRY group had the lowest. This was accompanied by the presence of correlation patterns including either orders belonging to various phyla from WET samples or orders belonging to the Actinobacteriota and Firmicutes phyla from DRY group samples. The phylogenetic analysis of the dominant species in WET and DRY samples showed that Acidobacteriota and Actinobacteriota strains were affiliated with uncultured bacteria retrieved from various oligotrophic and silica/quartz-rich environments, not only associated with subterranean sites. Our results suggest that the water content greatly contributes to shaping the microbial diversity within a subterranean quartzite environment. Further, the phylogenetic affiliation between Imawarì Yeuta dominant microbes and reference strains retrieved from both surface and subsurface silica- and/or CO₂/CO-rich environments, underlines the selective pressure applied by quartz as rock substrate. Oligotrophy probably in association with the geochemistry of silica/quartz low pH buffering activity and alternative energy sources led to the colonization of specific silica-associated microorganisms. This study provides clues for a better comprehension of the poorly known microbial life in subsurface and surface quartz-dominated environments.

Heterotrophic Bacteria Dominate Catalase Expression during Microcystis Blooms

In the oligotrophic oceans, key autotrophs depend on "helper" bacteria to reduce oxidative stress from hydrogen peroxide (H₂O₂) in the extracellular environment. H₂O₂ is also a ubiquitous stressor in freshwaters, but the effects of H₂O₂ on autotrophs and their interactions with bacteria are less well understood in freshwaters. Naturally occurring H₂O₂ in freshwater systems is proposed to impact the proportion of microcystin-producing (toxic) and non-microcystin-producing (nontoxic) Microcystis in blooms, which influences toxin concentrations and human health impacts. However, how different strains of Microcystis respond to naturally occurring H₂O₂ concentrations and the microbes responsible for H₂O₂ decomposition in freshwater cyanobacterial blooms are unknown. To address these knowledge gaps, we used metagenomics and metatranscriptomics to track the presence and expression of genes for H₂O₂ decomposition by microbes during a cyanobacterial bloom in western Lake Erie in the summer of 2014. katG encodes the key enzyme for decomposing extracellular H₂O₂ but was absent in most Microcystis cells. katG transcript relative abundance was dominated by heterotrophic bacteria. In axenic Microcystis cultures, an H₂O₂ scavenger (pyruvate) significantly improved growth rates of one toxic strain while other toxic and nontoxic strains were unaffected. These results indicate that heterotrophic bacteria play a key role in H₂O₂ decomposition in Microcystis blooms and suggest that their activity may affect the fitness of some Microcystis strains and thus the strain composition of Microcystis blooms but not along a toxic versus nontoxic dichotomy. IMPORTANCE Cyanobacterial harmful algal blooms (CHABs) threaten freshwater ecosystems globally through the production of toxins. Toxin production by cyanobacterial species and strains during CHABs varies widely over time and space, but the ecological drivers of the succession of toxin-producing species remain unclear. Hydrogen peroxide (H₂O₂) is ubiquitous in natural waters, inhibits microbial growth, and may determine the relative proportions of Microcystis strains during blooms. However, the mechanisms and organismal interactions involved in H₂O₂ decomposition are unexplored in CHABs. This study shows that some strains of bloom-forming freshwater cyanobacteria benefit from detoxification of H₂O₂ by associated heterotrophic bacteria, which may impact bloom development.

Uptake of Phytoplankton-Derived Carbon and Cobalamins by Novel Acidobacteria Genera in Microcystis Blooms Inferred from Metagenomic and Metatranscriptomic Evidence

Interactions between bacteria and phytoplankton can influence primary production, community composition, and algal bloom development. However, these interactions are poorly described for many consortia, particularly for freshwater bloom-forming cyanobacteria. Here, we assessed the gene content and expression of two uncultivated Acidobacteria from Lake Erie Microcystis blooms. These organisms were targeted because they were previously identified as important catalase producers in Microcystis blooms, suggesting that they protect Microcystis from H₂O₂. Metatranscriptomics revealed that both Acidobacteria transcribed genes for uptake of organic compounds that are known cyanobacterial products and exudates, including lactate, glycolate, amino acids, peptides, and cobalamins. Expressed genes for amino acid metabolism and peptide transport and degradation suggest that use of amino acids and peptides by Acidobacteria may regenerate nitrogen for cyanobacteria and other organisms. The Acidobacteria genomes lacked genes for biosynthesis of cobalamins but expressed genes for its transport and remodeling. This indicates that the Acidobacteria obtained cobalamins externally, potentially from Microcystis, which has a complete gene repertoire for pseudocobalamin biosynthesis; expressed them in field samples; and produced pseudocobalamin in axenic culture. Both Acidobacteria were detected in Microcystis blooms worldwide. Together, the data support the hypotheses that uncultured and previously unidentified Acidobacteria taxa exchange metabolites with phytoplankton during harmful cyanobacterial blooms and influence nitrogen available to phytoplankton. Thus, novel Acidobacteria may play a role in cyanobacterial physiology and bloom development. IMPORTANCE Interactions between heterotrophic bacteria and phytoplankton influence competition and successions between phytoplankton taxa, thereby influencing ecosystem-wide processes such as carbon cycling and algal bloom development. The cyanobacterium Microcystis forms harmful blooms in freshwaters worldwide and grows in buoyant colonies that harbor other bacteria in their phycospheres. Bacteria in the phycosphere and in the surrounding community likely influence Microcystis physiology and ecology and thus the development of freshwater harmful cyanobacterial blooms. However, the impacts and mechanisms of interaction between bacteria and Microcystis are not fully understood. This study explores the mechanisms of interaction between Microcystis and uncultured members of its phycosphere in situ with population genome resolution to investigate the cooccurrence of Microcystis and freshwater Acidobacteria in blooms worldwide.

Novel Microorganisms Contribute to Biosulfidogenesis in the Deep Layer of an Acidic Pit Lake

Cueva de la Mora is a permanently stratified acidic pit lake with extremely high concentrations of heavy metals at depth. In order to evaluate the potential for in situ sulfide production, we characterized the microbial community in the deep layer using metagenomics and metatranscriptomics. We retrieved 18 high quality metagenome-assembled genomes (MAGs) representing the most abundant populations. None of the MAGs were closely related to either cultured or non-cultured organisms from the Genome Taxonomy or NCBI databases (none with average nucleotide identity >95%). Despite oxygen concentrations that are consistently below detection in the deep layer, some archaeal and bacterial MAGs mapped transcripts of genes for sulfide oxidation coupled with oxygen reduction. Among these microaerophilic sulfide oxidizers, mixotrophic Thermoplasmatales archaea were the most numerous and represented 24% of the total community. Populations associated with the highest predicted in situ activity for sulfate reduction were affiliated with Actinobacteria, Chloroflexi, and Nitrospirae phyla, and together represented about 9% of the total community. These MAGs, in addition to a less abundant Proteobacteria MAG in the genus Desulfomonile, contained transcripts of genes in the Wood-Ljungdahl pathway. All MAGs had significant genetic potential for organic carbon oxidation. Our results indicate that novel acidophiles are contributing to biosulfidogenesis in the deep layer of Cueva de la Mora, and that in situ sulfide production is limited by organic carbon availability and sulfur oxidation.

RiboTaxa: combined approaches for rRNA genes taxonomic resolution down to the species level from metagenomics data revealing novelties

Metagenomic classifiers are widely used for the taxonomic profiling of metagenomics data and estimation of taxa relative abundance. Small subunit rRNA genes are a gold standard for phylogenetic resolution of microbiota, although the power of this marker comes down to its use as full-length. We aimed at identifying the tools that can efficiently lead to taxonomic resolution down to the species level. To reach this goal, we benchmarked the performance and accuracy of rRNA-specialized versus general-purpose read mappers, reference-targeted assemblers and taxonomic classifiers. We then compiled the best tools (BBTools, FastQC, SortMeRNA, MetaRib, EMIRGE, VSEARCH, BBMap and QIIME 2’s Sklearn classifier) to build a pipeline called RiboTaxa. Using metagenomics datasets, RiboTaxa gave the best results compared to other tools (i.e. Kraken2, Centrifuge, METAXA2, phyloFlash, SPINGO, BLCA, MEGAN) with precise taxonomic identification and relative abundance description without false positive detection (F-measure of 100% and 83.7% at genus level and species level, respectively). Using real datasets from various environments (i.e. ocean, soil, human gut) and from different approaches (e.g. metagenomics and gene capture by hybridization), RiboTaxa revealed microbial novelties not discerned by current bioinformatics analysis opening new biological perspectives in human and environmental health.

MarkerMAG: linking metagenome-assembled genomes (MAGs) with 16S rRNA marker genes using paired-end short reads

MOTIVATION

Metagenome-assembled genomes (MAGs) have substantially extended our understanding of microbial functionality. However, 16S rRNA genes, which are commonly used in phylogenetic analysis and environmental surveys, are often missing from MAGs. Here, we developed MarkerMAG, a pipeline that links 16S rRNA genes to MAGs using paired-end sequencing reads.

RESULTS

Assessment of MarkerMAG on three benchmarking metagenomic datasets with various degrees of complexity shows substantial increases in the number of MAGs with 16S rRNA genes and a 100% assignment accuracy. MarkerMAG also estimates the copy number of 16S rRNA genes in MAGs with high accuracy. Assessments on three real metagenomic datasets demonstrates 1.1- to 14.2-fold increases in the number of MAGs with 16S rRNA genes. We also show that MarkerMAG-improved MAGs increase the accuracy of functional prediction from 16S rRNA gene amplicon data. MarkerMAG is helpful in connecting information in MAG database with those in 16S rRNA databases and surveys and hence contributes to our increasing understanding of microbial diversity, function, and phylogeny.

AVAILABILITY

MarkerMAG is implemented in Python3 and freely available at https://github.com/songweizhi/MarkerMAG.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification of nosZ-expressing microorganisms consuming trace N2O in microaerobic chemostat consortia dominated by an uncultured Burkholderiales

Microorganisms possessing N₂O reductases (NosZ) are the only known environmental sink of N₂O. While oxygen inhibition of NosZ activity is widely known, environments where N₂O reduction occurs are often not devoid of O₂. However, little is known regarding N₂O reduction in microoxic systems. Here, 1.6-L chemostat cultures inoculated with activated sludge samples were sustained for ca. 100 days with low concentration (<2 ppmv) and feed rate (<1.44 µmoles h⁻¹) of N₂O, and the resulting microbial consortia were analyzed via quantitative PCR (qPCR) and metagenomic/metatranscriptomic analyses. Unintended but quantified intrusion of O₂ sustained dissolved oxygen concentration above 4 µM; however, complete N₂O reduction of influent N₂O persisted throughout incubation. Metagenomic investigations indicated that the microbiomes were dominated by an uncultured taxon affiliated to Burkholderiales, and, along with the qPCR results, suggested coexistence of clade I and II N₂O reducers. Contrastingly, metatranscriptomic nosZ pools were dominated by the Dechloromonas-like nosZ subclade, suggesting the importance of the microorganisms possessing this nosZ subclade in reduction of trace N₂O. Further, co-expression of nosZ and ccoNO/cydAB genes found in the metagenome-assembled genomes representing these putative N₂O-reducers implies a survival strategy to maximize utilization of scarcely available electron acceptors in microoxic environmental niches.

Metabolic Diversity and Aero-Tolerance in Anammox Bacteria from Geochemically Distinct Aquifers

Anaerobic ammonium oxidation (anammox) is important for converting bioavailable nitrogen into dinitrogen gas, particularly in carbon-poor environments. However, the diversity and prevalence of anammox bacteria in the terrestrial subsurface-a typically oligotrophic environment-are little understood. To determine the distribution and activity of anammox bacteria across a range of aquifer lithologies and physicochemistries, we analyzed 16S rRNA genes and quantified hydrazine synthase genes and transcripts sampled from 59 groundwater wells and metagenomes and metatranscriptomes from an oxic-to-dysoxic subset. Data indicate that anammox and anammox-associated bacteria (class "Candidatus Brocadiae") are prevalent in the aquifers studied, and that anammox community composition is strongly differentiated by dissolved oxygen (DO), but not ammonia/nitrite. While "Candidatus Brocadiae" diversity decreased with increasing DO, "Candidatus Brocadiae" 16S rRNA genes and hydrazine synthase (hzsB) genes and transcripts were detected across a wide range of bulk groundwater DO concentrations (0 to 10 mg/L). Anammox genes and transcripts correlated significantly with those involved in aerobic ammonia oxidation (amoA), potentially representing a major source of nitrite for anammox. Eight "Candidatus Brocadiae" genomes (63 to 95% complete), representing 2 uncharacterized families and 6 novel species, were reconstructed. Six genomes have genes characteristic of anammox, all for chemolithoautotrophy. Anammox and aerotolerance genes of up to four "Candidatus Brocadiae" genomes were transcriptionally active under oxic and dysoxic conditions, although activity was highest in dysoxic groundwater. The coexpression of nrfAH nitrite reductase genes by "Candidatus Brocadiae" suggests active regeneration of ammonia for anammox. Our findings indicate that anammox bacteria contribute to loss of fixed N across diverse anoxic-to-oxic aquifer conditions, which is likely supported by nitrite from aerobic ammonia oxidation. IMPORTANCE Anammox is increasingly shown to play a major role in the aquatic nitrogen cycle and can outcompete heterotrophic denitrification in environments low in organic carbon. Given that aquifers are characteristically oligotrophic, anammox may represent a major route for the removal of fixed nitrogen in these environments, including agricultural nitrogen, a common groundwater contaminant. Our research confirms that anammox bacteria and the anammox process are prevalent in aquifers and occur across diverse lithologies (e.g., sandy gravel, sand-silt, and volcanic) and groundwater physicochemistries (e.g., various oxygen, carbon, nitrate, and ammonium concentrations). Results reveal niche differentiation among anammox bacteria largely driven by groundwater oxygen contents and provide evidence that anammox is supported by proximity to oxic niches and handoffs from aerobic ammonia oxidizers. We further show that this process, while anaerobic, is active in groundwater characterized as oxic, likely due to the availability of anoxic niches.

Incorporating genome-based phylogeny and functional similarity into diversity assessments helps to resolve a global collection of human gut metagenomes

Tree-based diversity measures incorporate phylogenetic or functional relatedness into comparisons of microbial communities. This can improve the identification of explanatory factors compared to tree-agnostic diversity measures. However, applying tree-based diversity measures to metagenome data is more challenging than for single-locus sequencing (e.g., 16S rRNA gene). Utilizing the Genome Taxonomy Database (GTDB) for species-level metagenome profiling allows for functional diversity measures based on genomic content or traits inferred from it. Still, it is unclear how metagenome-based assessments of microbiome diversity benefit from incorporating phylogeny or function into measures of diversity. We assessed this by measuring phylogeny-based, function-based, and tree-agnostic diversity measures from a large, global collection of human gut metagenomes composed of 30 studies and 2943 samples. We found tree-based measures to explain phenotypic variation (e.g., westernization, disease status, and gender) better or equivalent to tree-agnostic measures. Ecophylogenetic and functional diversity measures provided unique insight into how microbiome diversity was partitioned by phenotype. Tree-based measures greatly improved machine learning model performance for predicting westernization, disease status, and gender, relative to models trained solely on tree-agnostic measures. Our findings illustrate the usefulness of tree- and function-based measures for metagenomic assessments of microbial diversity, which is a fundamental component of microbiome science. This article is protected by copyright. All rights reserved.

RNA Virosphere in a Marine Zooplankton Community in the Subtropical Western North Pacific

Zooplankton and viruses play a key role in marine ecosystems; however, their interactions have not been examined in detail. In the present study, the diversity of viruses associated with zooplankton collected using a plankton net (mesh size: 100‍ ‍μm) in the subtropical western North Pacific was investigated by fragmented and primer ligated dsRNA sequencing. We obtained 21 and 168 operational taxonomic units (OTUs) of ssRNA and dsRNA viruses, respectively, containing RNA-dependent RNA polymerase (RdRp). These OTUs presented average amino acid similarities of 43.5 and 44.0% to the RdRp genes of known viruses in ssRNA viruses and dsRNA viruses, respectively. Dominant OTUs mainly belonged to narna-like and picorna-like ssRNA viruses and chryso-like, partiti-like, picobirna-like, reo-like, and toti-like dsRNA viruses. Phylogenetic ana-lyses of the RdRp gene revealed that OTUs were phylogenetically diverse and clustered into distinct clades from known viral groups. The community structure of the same zooplankton sample was investigated using small subunit (SSU) rRNA sequences assembled from the metatranscriptome of single-stranded RNA. More than 90% of the sequence reads were derived from metazoan zooplankton; copepods comprised approximately 70% of the sequence reads. Although this ana-lysis provided no direct evidence of the host species of RNA viruses, these dominant zooplankton are expected to be associated with the RNA viruses detected in the present study. The present results indicate that zooplankton function as a reservoir of diverse RNA viruses and suggest that investigations of zooplankton viruses will provide a more detailed understanding of the role of viruses in marine ecosystems.

Discrimination of Bacterial Community Structures among Healthy, Gingivitis, and Periodontitis Statuses through Integrated Metatranscriptomic and Network Analyses

Periodontal disease is an inflammatory condition caused by polymicrobial infection. The inflammation is initiated at the gingiva (gingivitis) and then extends to the alveolar bone, leading to tooth loss (periodontitis). Previous studies have shown differences in bacterial composition between periodontal healthy and diseased sites. However, bacterial metabolic activities during the health-to-periodontitis microbiome shift are still inadequately understood. This study was performed to investigate the bacterial characteristics of healthy, gingivitis, and periodontitis statuses through metatranscriptomic analysis. Subgingival plaque samples of healthy, gingivitis, and periodontitis sites in the same oral cavity were collected from 21 patients. Bacterial compositions were then determined based on 16S rRNA reads; taxonomic and functional profiles derived from genes based on mRNA reads were estimated. The results showed clear differences in bacterial compositions and functional profiles between healthy and periodontitis sites. Co-occurrence networks were constructed for each group by connecting two bacterial species if their mRNA abundances were positively correlated. The clustering coefficient values were 0.536 for healthy, 0.600 for gingivitis, and 0.371 for periodontitis sites; thus, network complexity increased during gingivitis development, whereas it decreased during progression to periodontitis. Taxa, including Eubacterium nodatum, Eubacterium saphenum, Filifactor alocis, and Fretibacterium fastidiosum, showed greater transcriptional activities than those of red complex bacteria, in conjunction with disease progression. These taxa were associated with periodontal disease progression, and the health-to-periodontitis microbiome shift was accompanied by alterations in bacterial network structure and complexity.

IMPORTANCE The characteristics of the periodontal microbiome influence clinical periodontal status. Gingivitis involves reversible gingival inflammation without alveolar bone resorption. In contrast, periodontitis is an irreversible disease characterized by inflammatory destruction in both soft and hard tissues. An imbalance of the microbiome is present in both gingivitis and periodontitis. However, differences in microbiomes and their functional activities in the healthy, gingivitis, and periodontitis statuses are still inadequately understood. Furthermore, some inflamed gingival statuses do not consistently cause attachment loss. In this study, metatranscriptomic analyses were used to investigate the specific bacterial composition and gene expression patterns of the microbiomes of the healthy, gingivitis, and periodontitis statuses. In addition, co-occurrence network analysis revealed that the gingivitis site included features of networks observed in both the healthy and periodontitis sites. These results provide transcriptomic evidence to support gingivitis as an intermediate state between the healthy and periodontitis statuses.

Revealing microbial species diversity using sequence capture by hybridization

Targeting small parts of the 16S rDNA phylogenetic marker by metabarcoding reveals microorganisms of interest but cannot achieve a taxonomic resolution at the species level, precluding further precise characterizations. To identify species behind operational taxonomic units (OTUs) of interest, even in the rare biosphere, we developed an innovative strategy using gene capture by hybridization. From three OTU sequences detected upon polyphenol supplementation and belonging to the rare biosphere of the human gut microbiota, we revealed 59 nearly full-length 16S rRNA genes, highlighting high bacterial diversity hidden behind OTUs while evidencing novel taxa. Inside each OTU, revealed 16S rDNA sequences could be highly distant from each other with similarities down to 85 %. We identified one new family belonging to the order Clostridiales, 39 new genera and 52 novel species. Related bacteria potentially involved in polyphenol degradation have also been identified through genome mining and our results suggest that the human gut microbiota could be much more diverse than previously thought.

Community RNA-Seq: multi-kingdom responses to living versus decaying roots in soil

Roots are a primary source of organic carbon input in most soils. The consumption of living and detrital root inputs involves multi-trophic processes and multiple kingdoms of microbial life, but typical microbial ecology studies focus on only one or two major lineages. We used Illumina shotgun RNA sequencing to conduct PCR-independent SSU rRNA community analysis ("community RNA-Seq") and simultaneously assess the bacteria, archaea, fungi, and microfauna surrounding both living and decomposing roots of the annual grass, Avena fatua. Plants were grown in 13CO2-labeled microcosms amended with 15N-root litter to identify the preferences of rhizosphere organisms for root exudates (13C) versus decaying root biomass (15N) using NanoSIMS microarray imaging (Chip-SIP). When litter was available, rhizosphere and bulk soil had significantly more Amoebozoa, which are potentially important yet often overlooked top-down drivers of detritusphere community dynamics and nutrient cycling. Bulk soil containing litter was depleted in Actinobacteria but had significantly more Bacteroidetes and Proteobacteria. While Actinobacteria were abundant in the rhizosphere, Chip-SIP showed Actinobacteria preferentially incorporated litter relative to root exudates, indicating this group's more prominent role in detritus elemental cycling in the rhizosphere. Our results emphasize that decomposition is a multi-trophic process involving complex interactions, and our methodology can be used to track the trajectory of carbon through multi-kingdom soil food webs.

Multispecies populations of methanotrophic Methyloprofundus and cultivation of a likely dominant species from the Iheya North deep-sea hydrothermal field

The Methyloprofundus clade is represented by uncultivated methanotrophic bacterial endosymbionts of deep-sea bathymodiolin mussels, but only a single free-living species has been cultivated to date. This study reveals the existence of free-living Methyloprofundus variants in the Iheya North deep-sea hydrothermal field in the mid-Okinawa Trough. A clade-targeted amplicon analysis of the particulate methane monooxygenase gene (pmoA) detected 647 amplicon sequence variants (ASVs) of the Methyloprofundus clade in microbial communities newly formed in in situ colonization systems. Such systems were deployed at colonies of bathymodiolin mussels and a galatheoid crab in diffuse-flow areas. These ASVs were classified into 161 species-like groups. The proportion of the species-like groups representing endosymbionts of mussels was unexpectedly low. A methanotrophic bacterium designated INp10, a likely dominant species in the Methyloprofundus population in this field, was enriched in a biofilm formed in a methane-fed cultivation system operated at 10°C. Genomic characterization with the gene transcription data set of INp10 from the biofilm suggested traits advantageous to niche competition in environments, such as mobility, chemotaxis, biofilm formation, offensive and defensive systems, and hypoxia tolerance. The notable metabolic traits that INp10 shares with some Methyloprofundus members are the use of lanthanide-dependent XoxF as the sole methanol dehydrogenase due to the absence of the canonical MxaFI, the glycolytic pathway using fructose-6-phosphate aldolase instead of fructose-1,6-bisphosphate aldolase, and the potential to perform partial denitrification from nitrate under oxygen-limited conditions. These findings help us better understand the ecological strategies of this possibly widespread marine-specific methanotrophic clade.

IMPORTANCE The Iheya North deep-sea hydrothermal field in the mid-Okinawa Trough is characterized by abundant methane derived from organic-rich sediments and diverse chemosynthetic animal species, including those harboring methanotrophic bacterial symbionts, such as bathymodiolin mussels Bathymodiolus japonicus and “Bathymodiolus” platifrons and a galatheoid crab, Shinkaia crosnieri. Symbiotic methanotrophs have attracted significant attention, and yet free-living methanotrophs in this environment have not been studied in detail. We focused on the free-living Methyloprofundus spp. that thrive in this hydrothermal field and identified an unexpectedly large number of species-like groups in this clade. Moreover, we enriched and characterized a methanotroph whose genome sequence indicated that it corresponds to a new species in the genus Methyloprofundus. This species might be a dominant member of the indigenous Methyloprofundus population. New information on free-living Methyloprofundus populations suggests that the hydrothermal field is a promising locale at which to investigate the adaptive capacity and associated genetic diversity of Methyloprofundus spp.

Oligotyping and metagenomics reveal distinct Candidatus Accumulibacter communities in side-stream versus conventional full-scale enhanced biological phosphorus removal (EBPR) systems

Candidatus Accumulibacter phosphatis (CAP) and its clade-level micro-diversity has been associated with and implicated in functional differences in phosphorus removal performance in enhanced biological phosphorus removal (EBPR) systems. Side-stream EBPR (S2EBPR) is an emerging process that has been shown to present a suite of advantages over the conventional EBPR design, however, large knowledge gaps remain in terms of its underlying ecological mechanisms. Here, we compared and revealed the higher-resolution differences in microbial ecology of CAP between a full-scale side-stream EBPR configuration and a conventional A2O EBPR process that were operated in parallel and with the same influent feed. Even though the relative abundance of CAP, revealed by 16S rRNA gene amplicon sequencing, was similar in both treatment trains, a clade-level analysis, using combined 16S rRNA-gene based amplicon sequencing and oligotyping analysis and metagenomics analysis, revealed the distinct CAP microdiversity between the S2EBPR and A2O configurations that likely attributed to the improved performance in S2EBPR in comparison to conventional EBPR. Furthermore, genome-resolved metagenomics enabled extraction of three metagenome-assembled genomes (MAGs) belonging to CAP clades IIB (RCAB4-2), IIC (RC14) and II (RC18), from full-scale EBPR sludge for the first time, including a distinct Ca. Accumulibacter clade that is dominant and associated only with the S2EBPR configuration. The results also revealed the temporally increasing predominance of RC14, which belonged to Clade IIC, during the implementation of the S2EBPR configuration. Finally, we also show the existence of previously uncharacterized diversity of clades of CAP, namely the clades IIB and as yet unidentified clade of type II, in full-scale EBPR communities, highlighting the unknown diversity of CAP communities in full-scale EBPR systems.

Genome Streamlining, Plasticity, and Metabolic Versatility Distinguish Co-occurring Toxic and Nontoxic Cyanobacterial Strains of Microcoleus

Harmful cyanobacterial bloom occurrences have increased worldwide due to climate change and eutrophication, causing nuisance and animal deaths. Species from the benthic cyanobacterial genus Microcoleus are ubiquitous and form thick mats in freshwater systems, such as rivers, that are sometimes toxic due to the production of potent neurotoxins (anatoxins). Anatoxin-producing (toxic) strains typically coexist with non-anatoxin-producing (nontoxic) strains in mats, although the reason for this is unclear. To determine the genetic mechanisms differentiating toxic and nontoxic Microcoleus, we sequenced and assembled genomes from 11 cultures and compared these to another 31 Microcoleus genomes. Average nucleotide identities (ANI) indicate that toxic and nontoxic strains are distinct species (ANI, <95%), and only 6% of genes are shared across all 42 genomes, suggesting a high level of genetic divergence among Microcoleus strains. Comparative genomics showed substantial genome streamlining in toxic strains and a potential dependency on external sources for thiamine and sucrose. Toxic and nontoxic strains are further differentiated by an additional set of putative nitrate transporter (nitrogen uptake) and cyanophycin (carbon and nitrogen storage) genes, respectively. These genes likely confer distinct competitive advantages based on nutrient availability and suggest nontoxic strains are more robust to nutrient fluctuations. Nontoxic strains also possess twice as many transposable elements, potentially facilitating greater genetic adaptation to environmental changes. Our results offer insights into the divergent evolution of Microcoleus strains and the potential for cooperative and competitive interactions that contribute to the co-occurrence of toxic and nontoxic species within mats. IMPORTANCE Microcoleus autumnalis, and closely related Microcoleus species, compose a geographically widespread group of freshwater benthic cyanobacteria. Canine deaths due to anatoxin-a poisoning, following exposure to toxic proliferations, have been reported globally. While Microcoleus proliferations are on the rise, the mechanisms underpinning competition between, or coexistence of, toxic and nontoxic strains are unknown. This study identifies substantial genetic differences between anatoxin-producing and non-anatoxin-producing strains, pointing to reduced metabolic flexibility in toxic strains, and potential dependence on cohabiting nontoxic strains. Results provide insights into the metabolic and evolutionary differences between toxic and nontoxic Microcoleus, which may assist in predicting and managing aquatic proliferations.

Music of metagenomics-a review of its applications, analysis pipeline, and associated tools

This humble effort highlights the intricate details of metagenomics in a simple, poetic, and rhythmic way. The paper enforces the significance of the research area, provides details about major analytical methods, examines the taxonomy and assembly of genomes, emphasizes some tools, and concludes by celebrating the richness of the ecosystem populated by the "metagenome."

Metagenomic insights into co-proliferation of Vibrio spp. and dinoflagellates Prorocentrum during a spring algal bloom in the coastal East China Sea

Coastal harmful algal blooms (HABs), commonly termed 'red tides', have severe undesirable consequences to the marine ecosystems and local fishery and tourism industries. Increase in nitrogen and/or phosphorus loading is often regarded as the major culprits of increasing frequency and intensity of the coastal HAB; however, fundamental understanding is lacking as to the causes and mechanism of bloom formation despite decades of intensive investigation. In this study, we interrogated the prokaryotic microbiomes of surface water samples collected at two neighboring segments of East China Sea that contrast greatly in terms of the intensity and frequency of Prorocentrum-dominated HAB. Mantel tests identified significant correlations between the structural and functional composition of the microbiomes and the physicochemical state and the algal biomass density of the surface seawater, implying the possibility that prokaryotic microbiota may play key roles in the coastal HAB. A conspicuous feature of the microbiomes at the sites characterized with high trophic state index and eukaryotic algal cell counts was disproportionate proliferation of Vibrio spp., and their complete domination of the functional genes attributable to the dissimilatory nitrate reduction to ammonia (DNRA) pathway substantially enriched at these sites. The genes attributed to phosphorus uptake function were significantly enriched at these sites, presumably due to the P_i-deficiency induced by algal growth; however, the profiles of the phosphorus mineralization genes lacked consistency, barring any conclusive evidence with regard to contribution of prokaryotic microbiota to phosphorus bioavailability. The results of the co-occurrence network analysis performed with the core prokaryotic microbiome supported that the observed proliferation of Vibrio and HAB may be causally associated. The findings of this study suggest a previously unidentified association between Vibrio proliferation and the Prorocentrum-dominated HAB in the subtropical East China Sea, and opens a discussion regarding a theoretically unlikely, but still possible, involvement of Vibrio-mediated DNRA in Vibrio-Prorocentrum symbiosis. Further experimental substantiation of this supposed symbiotic mechanism may prove crucial in understanding the dynamics of explosive local algal growth in the region during spring algal blooms.

Microbial river-to-sea continuum: gradients in benthic and planktonic diversity, osmoregulation and nutrient cycling

BACKGROUND

Coastal aquatic ecosystems include chemically distinct, but highly interconnected environments. Across a freshwater-to-marine transect, aquatic communities are exposed to large variations in salinity and nutrient availability as tidal cycles create periodic fluctuations in local conditions. These factors are predicted to strongly influence the resident microbial community structure and functioning, and alter the structure of aquatic food webs and biogeochemical cycles. Nevertheless, little is known about the spatial distribution of metabolic properties across salinity gradients, and no study has simultaneously surveyed the sediment and water environments. Here, we determined patterns and drivers of benthic and planktonic prokaryotic and microeukaryotic community assembly across a river and tidal lagoon system by collecting sediments and planktonic biomass at nine shallow subtidal sites in the summer. Genomic and transcriptomic analyses, alongside a suite of complementary geochemical data, were used to determine patterns in the distribution of taxa, mechanisms of salt tolerance, and nutrient cycling.

RESULTS

Taxonomic and metabolic profiles related to salt tolerance and nutrient cycling of the aquatic microbiome were found to decrease in similarity with increasing salinity, and distinct trends in diversity were observed between the water column and sediment. Non-saline and saline communities adopted divergent strategies for osmoregulation, with an increase in osmoregulation-related transcript expression as salinity increased in the water column due to lineage-specific adaptations to salt tolerance. Results indicated a transition from phosphate limitation in freshwater habitats to nutrient-rich conditions in the brackish zone, where distinct carbon, nitrogen and sulfur cycling processes dominated. Phosphorus acquisition-related activity was highest in the freshwater zone, along with dissimilatory nitrate reduction to ammonium in freshwater sediment. Activity associated with denitrification, sulfur metabolism and photosynthesis were instead highest in the brackish zone, where photosynthesis was dominated by distinct microeukaryotes in water (Cryptophyta) and sediment (diatoms). Despite microeukaryotes and archaea being rare relative to bacteria, results indicate that they contributed more to photosynthesis and ammonia oxidation, respectively.

CONCLUSIONS

Our study demonstrates clear freshwater-saline and sediment-water ecosystem boundaries in an interconnected coastal aquatic system and provides a framework for understanding the relative importance of salinity, planktonic-versus-benthic habitats and nutrient availability in shaping aquatic microbial metabolic processes, particularly in tidal lagoon systems. Video abstract.

Capture by hybridization for full-length barcode-based eukaryotic and prokaryotic biodiversity inventories of deep sea ecosystems

Biodiversity inventory of marine systems remains limited due to unbalanced access to the three ocean dimensions. The use of environmental DNA (eDNA) for metabarcoding allows fast and effective biodiversity inventory and is forecast as a future biodiversity research and biomonitoring tool. However, in poorly understood ecosystems, eDNA results remain difficult to interpret due to large gaps in reference databases and PCR bias limiting the detection of some major phyla. Here, we aimed to circumvent these limitations by avoiding PCR and recollecting larger DNA fragments to improve assignment of detected taxa through phylogenetic reconstruction. We applied capture by hybridization (CBH) to enrich DNA from deep-sea sediment samples and compared the results with those obtained through an up-to-date metabarcoding PCR-based approach (MTB). Originally developed for bacterial communities and targeting 16S rDNA, the CBH approach was applied to 18S rDNA to improve the detection of species forming benthic communities of eukaryotes, with a particular focus on metazoans. The results confirmed the possibility of extending CBH to metazoans with two major advantages: (i) CBH revealed a broader spectrum of prokaryotic, eukaryotic, and particularly metazoan diversity, and (ii) CBH allowed much more robust phylogenetic reconstructions of full-length barcodes with up to 1900 base pairs. This is particularly important for taxa whose assignment is hampered by gaps in reference databases. This study provides a database and probes to apply 18S CBH to diverse marine systems, confirming this promising new tool to improve biodiversity assessments in data-poor ecosystems such as those in the deep sea.

Seasonality affects function and complexity but not diversity of the rhizosphere microbiome in European temperate grassland

Knowledge on how grassland microbiota responds on gene expression level to winter-summer change of seasons is poor. Here, we used a combination of quantitative PCR-based assays and metatranscriptomics to assess the impact of seasonality on the rhizospheric microbiota in temperate European grassland. Bacteria dominated, being at least one order of magnitude more abundant than fungi. Despite a fivefold summer increase in bacterial community size, season had nearly no effect on microbiome diversity. It, however, had a marked impact on taxon-specific gene expression, with 668 genes significantly differing in relative transcript abundance between winter and summer samples. Acidobacteria, Bacteroidetes, Planctomycetes, and Proteobacteria showed a greater relative gene expression activity in winter, while mRNA of Actinobacteria and Fungi was, relative to other taxa, significantly enriched in summer. On functional level, mRNA involved in protein turnover (e.g., transcription and translation) and cell maintenance (e.g., chaperones that protect against cell freezing damage such as GroEL and Hsp20) were highly enriched in winter. By contrast, mRNA involved in central carbon and amino acid metabolisms had a greater abundance in summer. Among carbohydrate-active enzymes, transcripts of GH36 family (hemicellulases) were highly enriched in winter, while those encoding GH3 family (cellulases) showed increased abundance in summer. The seasonal differences in plant polymer breakdown were linked to a significantly greater microbial network complexity in winter than in summer. Conceptually, the winter-summer change in microbiome functioning can be well explained by a shift from stress-tolerator to high-yield life history strategy.

Multiomic Approach to Analyze Infant Gut Microbiota: Experimental and Analytical Method Optimization

Background: The human intestinal microbiome plays a central role in overall health status, especially in early life stages. 16S rRNA amplicon sequencing is used to profile its taxonomic composition; however, multiomic approaches have been proposed as the most accurate methods for study of the complexity of the gut microbiota. In this study, we propose an optimized method for bacterial diversity analysis that we validated and complemented with metabolomics by analyzing fecal samples. Methods: Forty-eight different analytical combinations regarding (1) 16S rRNA variable region sequencing, (2) a feature selection approach, and (3) taxonomy assignment methods were tested. A total of 18 infant fecal samples grouped depending on the type of feeding were analyzed by the proposed 16S rRNA workflow and by metabolomic analysis. Results: The results showed that the sole use of V4 region sequencing with ASV identification and VSEARCH for taxonomy assignment produced the most accurate results. The application of this workflow showed clear differences between fecal samples according to the type of feeding, which correlated with changes in the fecal metabolic profile. Conclusion: A multiomic approach using real fecal samples from 18 infants with different types of feeding demonstrated the effectiveness of the proposed 16S rRNA-amplicon sequencing workflow.

A synbiotic intervention modulates meta-omics signatures of gut redox potential and acidity in elective caesarean born infants

BACKGROUND

The compromised gut microbiome that results from C-section birth has been hypothesized as a risk factor for the development of non-communicable diseases (NCD). In a double-blind randomized controlled study, 153 infants born by elective C-section received an infant formula supplemented with either synbiotic, prebiotics, or unsupplemented from birth until 4 months old. Vaginally born infants were included as a reference group. Stool samples were collected from day 3 till week 22. Multi-omics were deployed to investigate the impact of mode of delivery and nutrition on the development of the infant gut microbiome, and uncover putative biological mechanisms underlying the role of a compromised microbiome as a risk factor for NCD.

RESULTS

As early as day 3, infants born vaginally presented a hypoxic and acidic gut environment characterized by an enrichment of strict anaerobes (Bifidobacteriaceae). Infants born by C-section presented the hallmark of a compromised microbiome driven by an enrichment of Enterobacteriaceae. This was associated with meta-omics signatures characteristic of a microbiome adapted to a more oxygen-rich gut environment, enriched with genes associated with reactive oxygen species metabolism and lipopolysaccharide biosynthesis, and depleted in genes involved in the metabolism of milk carbohydrates. The synbiotic formula modulated expression of microbial genes involved in (oligo)saccharide metabolism, which emulates the eco-physiological gut environment observed in vaginally born infants. The resulting hypoxic and acidic milieu prevented the establishment of a compromised microbiome.

CONCLUSIONS

This study deciphers the putative functional hallmarks of a compromised microbiome acquired during C-section birth, and the impact of nutrition that may counteract disturbed microbiome development.

TRIAL REGISTRATION

The study was registered in the Dutch Trial Register (Number: 2838 ) on 4th April 2011.

Strenuous Physical Training, Physical Fitness, Body Composition and Bacteroides to Prevotella Ratio in the Gut of Elderly Athletes

Regular physical activity seems to have a positive effect on the microbiota composition of the elderly, but little is known about the added possible benefits of strenuous endurance training. To gain insight into the physiology of the elderly and to identify biomarkers associated with endurance training, we combined different omics approaches. We aimed to investigate the gut microbiome, plasma composition, body composition, cardiorespiratory fitness, and muscle strength of lifetime elderly endurance athletes (LA) age 63.5 (95% CI 61.4, 65.7), height 177.2 (95% CI 174.4, 180.1) cm, weight 77.8 (95% CI 75.1, 80.5) kg, VO2max 42.4 (95% CI 39.8, 45.0) ml.kg^–1.min^–1 (n = 13) and healthy controls age 64.9 (95% CI 62.1, 67.7), height 174.9 (95% CI 171.2, 178.6) cm, weight 83.4 (95% CI 77.1, 89.7) kg, VO2max 28.9 (95% CI 23.9, 33.9), ml.kg^–1.min^–1 (n = 9). Microbiome analysis was performed on collected stool samples further subjected to 16S rRNA gene analysis. NMR-spectroscopic analysis was applied to determine and compare selected blood plasma metabolites mostly linked to energy metabolism. The machine learning (ML) analysis discriminated subjects from the LA and CTRL groups using the joint predictors Bacteroides 1.8E + 00 (95% CI 1.1, 2.5)%, 3.8E + 00 (95% CI 2.7, 4.8)% (p = 0.002); Prevotella 1.3 (95% CI 0.28, 2.4)%, 0.1 (95% CI 0.07, 0.3)% (p = 0.02); Intestinimonas 1.3E-02 (95% CI 9.3E-03, 1.7E-02)%, 5.9E-03 (95% CI 3.9E-03, 7.9E-03)% (p = 0.002), Subdoligranulum 7.9E-02 (95% CI 2.5E-02, 1.3E-02)%, 3.2E-02 (95% CI 1.8E-02, 4.6E-02)% (p = 0.02); and the ratio of Bacteroides to Prevotella 133 (95% CI -86.2, 352), 732 (95% CI 385, 1079.3) (p = 0.03), leading to an ROC curve with AUC of 0.94. Further, random forest ML analysis identified VO2max, BMI, and the Bacteroides to Prevotella ratio as appropriate, joint predictors for discriminating between subjects from the LA and CTRL groups. Although lifelong endurance training does not bring any significant benefit regarding overall gut microbiota diversity, strenuous athletic training is associated with higher cardiorespiratory fitness, lower body fat, and some favorable gut microbiota composition, all factors associated with slowing the rate of biological aging.

Genomic adaptations enabling Acidithiobacillus distribution across wide-ranging hot spring temperatures and pHs

BACKGROUND

Terrestrial hot spring settings span a broad spectrum of physicochemistries. Physicochemical parameters, such as pH and temperature, are key factors influencing differences in microbial composition across diverse geothermal areas. Nonetheless, analysis of hot spring pools from the Taupo Volcanic Zone (TVZ), New Zealand, revealed that some members of the bacterial genus, Acidithiobacillus, are prevalent across wide ranges of hot spring pHs and temperatures. To determine the genomic attributes of Acidithiobacillus that inhabit such diverse conditions, we assembled the genomes of 19 uncultivated hot spring Acidithiobacillus strains from six geothermal areas and compared these to 37 publicly available Acidithiobacillus genomes from various habitats.

RESULTS

Analysis of 16S rRNA gene amplicons from 138 samples revealed that Acidithiobacillus comprised on average 11.4 ± 16.8% of hot spring prokaryotic communities, with three Acidithiobacillus amplicon sequence variants (ASVs) (TVZ_G1, TVZ_G2, TVZ_G3) accounting for > 90% of Acidithiobacillus in terms of relative abundance, and occurring in 126 out of 138 samples across wide ranges of temperature (17.5-92.9 °C) and pH (1.0-7.5). We recovered 19 environmental genomes belonging to each of these three ASVs, as well as a fourth related group (TVZ_G4). Based on genome average nucleotide identities, the four groups (TVZ_G1-TVZ_G4) constitute distinct species (ANI < 96.5%) of which three are novel Acidithiobacillus species (TVZ_G2-TVZ_G4) and one belongs to Acidithiobacillus caldus (TVZ_G1). All four TVZ Acidithiobacillus groups were found in hot springs with temperatures above the previously known limit for the genus (up to 40 °C higher), likely due to significantly higher proline and GC contents than other Acidithiobacillus species, which are known to increase thermostability. Results also indicate hot spring-associated Acidithiobacillus have undergone genome streamlining, likely due to thermal adaptation. Moreover, our data suggest that Acidithiobacillus prevalence across varied hot spring pHs is supported by distinct strategies, whereby TVZ_G2-TVZ_G4 regulate pH homeostasis mostly through Na⁺/H⁺ antiporters and proton-efflux ATPases, whereas TVZ_G1 mainly relies on amino acid decarboxylases.

CONCLUSIONS

This study provides insights into the distribution of Acidithiobacillus species across diverse hot spring physichochemistries and determines genomic features and adaptations that potentially enable Acidithiobacillus species to colonize a broad range of temperatures and pHs in geothermal environments. Video Abstract.

Evaluation of saliva as a source of accurate whole-genome and microbiome sequencing data

This study sets out to establish the suitability of saliva-based whole-genome sequencing (WGS) through a comparison against blood-based WGS. To fully appraise the observed differences, we developed a novel technique of pseudo-replication. We also investigated the potential of characterizing individual salivary microbiomes from non-human DNA fragments found in saliva. We observed that the majority of discordant genotype calls between blood and saliva fell into known regions of the human genome that are typically sequenced with low confidence; and could be identified by quality control measures. Pseudo-replication demonstrated that the levels of discordance between blood- and saliva-derived WGS data were entirely similar to what one would expect between technical replicates if an individual's blood or saliva had been sequenced twice. Finally, we successfully sequenced salivary microbiomes in parallel to human genomes as demonstrated by a comparison against the Human Microbiome Project.

Comparative genomics of the ADA clade within the Nostocales

The ADA clade of Nostocales cyanobacteria, a group that is prominent in current harmful algal bloom events, now includes over 40 genome sequences with the recent addition of sixteen novel sequenced genomes (Dreher et al., Harmful Algae, 2021). Fourteen genomes are complete (closed), enabling highly detailed assessments of gene content and genome architecture. ADA genomes contain 5 rRNA operons, genes expected to support a photoautotrophic and diazotrophic lifestyle, and a varied array of genes for the synthesis of bioactive secondary metabolites. Genes for the production of the taste-and-odor compound geosmin and the four major classes of cyanotoxins - anatoxin-a, cylindrospermopsin, microcystin and saxitoxin - are represented in members of the ADA clade. Notably, the gene array for the synthesis of cylindrospermopsin by Dolichospermum sp. DET69 was located on a plasmid, raising the possibility of facile horizontal transmission. However, genes supporting independent conjugative transfer of this plasmid are lacking. Further, analysis of genomic loci containing this and other cyanotoxin gene arrays shows evidence that these arrays have long-term stability and do not appear to be genomic islands easily capable of horizontal transmission to other cells. There is considerable diversity in the gene complements of individual ADA genomes, including the variable presence of physiologically important genes: genomes in three species-level subclades lack the gas vesicle genes that facilitate a planktonic lifestyle, and, surprisingly, the genome of Cuspidothrix issatschenkoi CHARLIE-1, a reported diazotroph, lacks the genes for nitrogen fixation. Notably, phylogenetically related genomes possess limited synteny, indicating a prominent role for chromosome rearrangements during ADA strain evolution. The genomes contain abundant insertion sequences and repetitive transposase genes, which could be the main drivers of genome rearrangement through active transposition and homologous recombination. No prophages were found, and no evidence of viral infection was observed in the bloom population samples from which the genomes discussed here were derived. Phages thus seem to have a limited influence on ADA evolution.

Opening up the Single-Cell Toolbox for Microbial Natural Products Research

The diverse microbes that produce natural products represent an important source of novel therapeutics, drug leads, and scientific tools. However, the vast majority have not been grown in axenic culture and are members of complex communities. While meta-'omic methods such as metagenomics, -transcriptomics, and -proteomics reveal collective molecular features of this "microbial dark matter", the study of individual microbiome members can be challenging. To address these limits, a number of techniques with single-bacterial resolution have been developed in the last decade and a half. While several of these are embraced by microbial ecologists, there has been less use by researchers interested in mining microbes for natural products. In this review, we discuss the available and emerging techniques for targeted single-cell analysis with a particular focus on applications to the discovery and study of natural products.

Viral RNA Genomes Identified from Marine Macroalgae and a Diatom

Protists provide insights into the diversity and function of RNA viruses in marine systems. Among them, marine macroalgae are good targets for RNA virome analyses because they have a sufficient biomass in nature. However, RNA viruses in macroalgae have not yet been examined in detail, and only partial genome sequences have been reported for the majority of RNA viruses. Therefore, to obtain further insights into the distribution and diversity of RNA viruses associated with marine protists, we herein examined RNA viruses in macroalgae and a diatom. We report the putative complete genome sequences of six novel RNA viruses from two marine macroalgae and one diatom holobiont. Four viruses were not classified into established viral genera or families. Furthermore, a virus classified into Totiviridae showed a genome structure that has not yet been reported in this family. These results suggest that a number of distinct RNA viruses are widespread in a broad range of protists.

The diversity and evolution of microbial dissimilatory phosphite oxidation

Phosphite is the most energetically favorable chemotrophic electron donor known, with a half-cell potential (Eo') of -650 mV for the PO43-/PO33- couple. Since the discovery of microbial dissimilatory phosphite oxidation (DPO) in 2000, the environmental distribution, evolution, and diversity of DPO microorganisms (DPOMs) have remained enigmatic, as only two species have been identified. Here, metagenomic sequencing of phosphite-enriched microbial communities enabled the genome reconstruction and metabolic characterization of 21 additional DPOMs. These DPOMs spanned six classes of bacteria, including the Negativicutes, Desulfotomaculia, Synergistia, Syntrophia, Desulfobacteria, and Desulfomonilia_A Comparing the DPO genes from the genomes of enriched organisms with over 17,000 publicly available metagenomes revealed the global existence of this metabolism in diverse anoxic environments, including wastewaters, sediments, and subsurface aquifers. Despite their newfound environmental and taxonomic diversity, metagenomic analyses suggested that the typical DPOM is a chemolithoautotroph that occupies low-oxygen environments and specializes in phosphite oxidation coupled to CO2 reduction. Phylogenetic analyses indicated that the DPO genes form a highly conserved cluster that likely has ancient origins predating the split of monoderm and diderm bacteria. By coupling microbial cultivation strategies with metagenomics, these studies highlighted the unsampled metabolic versatility latent in microbial communities. We have uncovered the unexpected prevalence, diversity, biochemical specialization, and ancient origins of a unique metabolism central to the redox cycling of phosphorus, a primary nutrient on Earth.

Transition from unclassified Ktedonobacterales to Actinobacteria during amorphous silica precipitation in a quartzite cave environment

The orthoquartzite Imawarì Yeuta cave hosts exceptional silica speleothems and represents a unique model system to study the geomicrobiology associated to silica amorphization processes under aphotic and stable physical-chemical conditions. In this study, three consecutive evolution steps in the formation of a peculiar blackish coralloid silica speleothem were studied using a combination of morphological, mineralogical/elemental and microbiological analyses. Microbial communities were characterized using Illumina sequencing of 16S rRNA gene and clone library analysis of carbon monoxide dehydrogenase (coxL) and hydrogenase (hypD) genes involved in atmospheric trace gases utilization. The first stage of the silica amorphization process was dominated by members of a still undescribed microbial lineage belonging to the Ktedonobacterales order, probably involved in the pioneering colonization of quartzitic environments. Actinobacteria of the Pseudonocardiaceae and Acidothermaceae families dominated the intermediate amorphous silica speleothem and the final coralloid silica speleothem, respectively. The atmospheric trace gases oxidizers mostly corresponded to the main bacterial taxa present in each speleothem stage. These results provide novel understanding of the microbial community structure accompanying amorphization processes and of coxL and hypD gene expression possibly driving atmospheric trace gases metabolism in dark oligotrophic caves.

Genomic Insights Into the Lifestyles of Thaumarchaeota Inside Sponges

Sponges are among the oldest metazoans and their success is partly due to their abundant and diverse microbial symbionts. They are one of the few animals that have Thaumarchaeota symbionts. Here we compare genomes of 11 Thaumarchaeota sponge symbionts, including three new genomes, to free-living ones. Like their free-living counterparts, sponge-associated Thaumarchaeota can oxidize ammonia, fix carbon, and produce several vitamins. Adaptions to life inside the sponge host include enrichment in transposases, toxin-antitoxin systems and restriction modifications systems, enrichments previously reported also from bacterial sponge symbionts. Most thaumarchaeal sponge symbionts lost the ability to synthesize rhamnose, which likely alters their cell surface and allows them to evade digestion by the host. All but one archaeal sponge symbiont encoded a high-affinity, branched-chain amino acid transporter system that was absent from the analyzed free-living thaumarchaeota suggesting a mixotrophic lifestyle for the sponge symbionts. Most of the other unique features found in sponge-associated Thaumarchaeota, were limited to only a few specific symbionts. These features included the presence of exopolyphosphatases and a glycine cleavage system found in the novel genomes. Thaumarchaeota have thus likely highly specific interactions with their sponge host, which is supported by the limited number of host sponge species to which each of these symbionts is restricted.

Usage of Metatranscriptomics to Understand Oral Disease

Metatranscriptomics is a method used to comprehensively capture bacterial activity within microbiota at the transcription level. It has become an alternative to the 16S rDNA sequencing, which uses only the 16S rRNA gene for predicting bacterial composition. By conducting metatranscriptomics, investigators can obtain substantial information about what types of genes are transcribed at the time of sampling and which bacterial taxa are responsible for their transcription. Here, I describe a protocol for metatranscriptomics for oral microbiota by using high-throughput sequencing technology. A remarkable feature of this protocol is that it uses the level of rRNA expression as the internal control for measuring transcriptional activity of each bacterial taxon. The normalized mRNA level is given by the mRNA/rRNA ratio, which indicates the extent of transcriptional activity.