Genome-resolved metagenomics reveals site-specific diversity of episymbiotic CPR bacteria and DPANN archaea in groundwater ecosystems

Gene horizontal transfers and functional diversity negatively correlated with bacterial taxonomic diversity along a nitrogen gradient

Horizontal gene transfer (HGT) mediated diversification is a critical force driving evolutionary and ecological processes. However, how HGT might relate to anthropogenic activity such as nitrogen addition, and its subsequent effect on functional diversity and cooccurrence networks remain unknown. Here we approach this knowledge gap by blending bacterial 16S rRNA gene amplicon and shotgun metagenomes from a platform of cessation of nitrogen additions and continuous nitrogen additions. We found that bacterial HGT events, functional genes, and virus diversities increased whereas bacterial taxonomic diversity decreased by nitrogen additions, resulting in a counterintuitive strong negative association between bacterial taxonomic and functional diversities. Nitrogen additions, especially the ceased one, complexified the cooccurrence network by increasing the contribution of vitamin B12 auxotrophic Acidobacteria, indicating cross-feeding. These findings advance our perceptions of the causes and consequences of the diversification process in community ecology.

Genetic potential for aerobic respiration and denitrification in globally distributed respiratory endosymbionts

The endosymbiont Candidatus Azoamicus ciliaticola was proposed to generate ATP for its eukaryotic host, an anaerobic ciliate of the Plagiopylea class, fulfilling a function analogous to mitochondria in other eukaryotic cells. The discovery of this respiratory endosymbiosis has major implications for both evolutionary history and ecology of microbial eukaryotes. However, with only a single species described, knowledge of its environmental distribution and diversity is limited. Here we report four complete, circular metagenome assembled genomes (cMAGs) representing respiratory endosymbionts inhabiting groundwater in California, Ohio, and Germany. These cMAGs form two lineages comprising a monophyletic clade within the uncharacterized gammaproteobacterial order UBA6186, enabling evolutionary analysis of their key protein complexes. Strikingly, all four cMAGs encode a cytochrome cbb3 oxidase, which indicates that these endosymbionts have the capacity for aerobic respiration. Accordingly, we detect these respiratory endosymbionts in diverse habitats worldwide, thus further expanding the ecological scope of this respiratory symbiosis.

Metagenomic characterization of viruses and mobile genetic elements associated with the DPANN archaeal superphylum

The archaeal superphylum DPANN (an acronym formed from the initials of the first five phyla discovered: Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanohaloarchaeota and Nanoarchaeota) is a group of ultrasmall symbionts able to survive in extreme ecosystems. The diversity and dynamics between DPANN archaea and their virome remain largely unknown. Here we use a metagenomic clustered regularly interspaced short palindromic repeats (CRISPR) screening approach to identify 97 globally distributed, non-redundant viruses and unclassified mobile genetic elements predicted to infect hosts across 8 DPANN phyla, including 7 viral groups not previously characterized. Genomic analysis suggests a diversity of viral morphologies including head-tailed, tailless icosahedral and spindle-shaped viruses with the potential to establish lytic, chronic or lysogenic infections. We also find evidence of a virally encoded Cas12f1 protein (probably originating from uncultured DPANN archaea) and a mini-CRISPR array, which could play a role in modulating host metabolism. Many metagenomes have virus-to-host ratios >10, indicating that DPANN viruses play an important role in controlling host populations. Overall, our study illuminates the underexplored diversity, functional repertoires and host interactions of the DPANN virome.

Eutrophication impacts the distribution and functional traits of viral communities in lakes

Viruses play a crucial role in aquatic ecosystems by regulating microbial composition and impacting biogeochemical cycling. While the response of viral diversity to the trophic status has been preliminarily explored in lake ecosystems, there is limited integrated exploration of the biogeography of viruses, host associations, and the auxiliary metabolic genes (AMGs), particularly for plateau lakes. Therefore, this research investigated the viral biogeography, virus-host association, and AMGs in the surface waters of 11 lakes varying in trophic levels (eutrophic and oligo-mesotrophic) in the Yunnan-Guizhou plateau region of China. A total of 73,105 viral operational taxonomic units were obtained from 11 samples, with 84.8 % remaining unannotated at the family level, indicating a predominance of novel viruses within these lakes. The most abundant viral family was Kyanoviridae (24.4 %), recognized as a common cyanophage. The vast majority of cyanobacteria and several eukaryotic algae were predicted as hosts for the viruses, with a lytic lifestyle predominating the life strategy of these cyanophages, implying the potential influence of the virus on algae. The viral community structure significantly correlated with both trophic status and the bacterial community. The structure equation model analysis revealed chlorophyll a was the primary factor affecting viral communities. Moreover, numerous AMGs linked to carbon metabolism, phosphorus metabolism, sulfur metabolism, and photosynthesis were found in these lakes, some of which showed virus preference for the trophic statuses, suggesting a vital role of the virus in driving biogeochemical cycling in the lake crossing different nutrient levels. In addition, a restricted presence of viruses was found to infect humans or harbor antibiotic resistance genes in the lakes, suggesting a subtle yet potential link to human health. Overall, these findings offer insights into the response of viral communities to eutrophication and their potential role in biogeochemical cycling and controlling algal propagation.

Metaproteogenomics resolution of a high-CO2 aquifer community reveals a complex cellular adaptation of groundwater Gracilibacteria to a host-dependent lifestyle

BACKGROUND

Bacteria of the candidate phyla radiation (CPR), constituting about 25% of the bacterial biodiversity, are characterized by small cell size and patchy genomes without complete key metabolic pathways, suggesting a symbiotic lifestyle. Gracilibacteria (BD1-5), which are part of the CPR branch, possess alternate coded genomes and have not yet been cultivated. The lifestyle of Gracilibacteria, their temporal dynamics, and activity in natural ecosystems, particularly in groundwater, has remained largely unexplored. Here, we aimed to investigate Gracilibacteria activity in situ and to discern their lifestyle based on expressed genes, using the metaproteogenome of Gracilibacteria as a function of time in the cold-water geyser Wallender Born in the Volcanic Eifel region in Germany.

RESULTS

We coupled genome-resolved metagenomics and metaproteomics to investigate a cold-water geyser microbial community enriched in Gracilibacteria across a 12-day time-series. Groundwater was collected and sequentially filtered to fraction CPR and other bacteria. Based on 725 Gbps of metagenomic data, 1129 different ribosomal protein S3 marker genes, and 751 high-quality genomes (123 population genomes after dereplication), we identified dominant bacteria belonging to Gallionellales and Gracilibacteria along with keystone microbes, which were low in genomic abundance but substantially contributing to proteomic abundance. Seven high-quality Gracilibacteria genomes showed typical limitations, such as limited amino acid or nucleotide synthesis, in their central metabolism but no co-occurrence with potential hosts. The genomes of these Gracilibacteria were encoded for a high number of proteins involved in cell to cell interaction, supporting the previously surmised host-dependent lifestyle, e.g., type IV and type II secretion system subunits, transporters, and features related to cell motility, which were also detected on protein level.

CONCLUSIONS

We here identified microbial keystone taxa in a high-CO2 aquifer, and revealed microbial dynamics of Gracilibacteria. Although Gracilibacteria in this ecosystem did not appear to target specific organisms in this ecosystem due to lack of co-occurrence despite enrichment on 0.2-µm filter fraction, we provide proteomic evidence for the complex machinery behind the host-dependent lifestyle of groundwater Gracilibacteria. Video Abstract.

N-linked protein glycosylation in Nanobdellati (formerly DPANN) archaea and their hosts

Members of the kingdom Nanobdellati, previously known as DPANN archaea, are characterized by ultrasmall cell sizes and reduced genomes. They primarily thrive through ectosymbiotic interactions with specific hosts in diverse environments. Recent successful cultivations have emphasized the importance of adhesion to host cells for understanding the ecophysiology of Nanobdellati. Cell adhesion is often mediated by cell surface carbohydrates, and in archaea, this may be facilitated by the glycosylated S-layer protein that typically coats their cell surface. In this study, we conducted glycoproteomic analyses on two co-cultures of Nanobdellati with their host archaea, as well as on pure cultures of both host and non-host archaea. Nanobdellati exhibited various glycoproteins, including archaellins and hypothetical proteins, with glycans that were structurally distinct from those of their hosts. This indicated that Nanobdellati autonomously synthesize their glycans for protein modifications probably using host-derived substrates, despite the high energy cost. Glycan modifications on Nanobdellati proteins consistently occurred on asparagine residues within the N-X-S/T sequon, consistent with patterns observed across archaea, bacteria, and eukaryotes. In both host and non-host archaea, S-layer proteins were commonly modified with hexose, N-acetylhexosamine, and sulfonated deoxyhexose. However, the N-glycan structures of host archaea, characterized by distinct sugars such as deoxyhexose, nonulosonate sugar, and pentose at the nonreducing ends, were implicated in enabling Nanobdellati to differentiate between host and non-host cells. Interestingly, the specific sugar, xylose, was eliminated from the N-glycan in a host archaeon when co-cultured with Nanobdella. These findings enhance our understanding of the role of protein glycosylation in archaeal interactions.IMPORTANCENanobdellati archaea, formerly known as DPANN, are phylogenetically diverse, widely distributed, and obligately ectosymbiotic. The molecular mechanisms by which Nanobdellati recognize and adhere to their specific hosts remain largely unexplored. Protein glycosylation, a fundamental biological mechanism observed across all domains of life, is often crucial for various cell-cell interactions. This study provides the first insights into the glycoproteome of Nanobdellati and their host and non-host archaea. We discovered that Nanobdellati autonomously synthesize glycans for protein modifications, probably utilizing substrates derived from their hosts. Additionally, we identified distinctive glycosylation patterns that suggest mechanisms through which Nanobdellati differentiate between host and non-host cells. This research significantly advances our understanding of the molecular basis of microbial interactions in extreme environments.

The fate of antibiotics and antibiotic resistance genes in Large-Scale chicken farm Environments: Preliminary view of the performance of National veterinary Antimicrobial use reduction Action in Guangdong, China

In 2018, China implemented the Veterinary Antimicrobial Use Reduction Action to curb the rapid development of antibiotic resistance (AR). However, the AR-related pollutions in animal farms after the reduction policy has been poorly investigated. Here, we performed a comprehensive investigation combining UPLC-MS/MS, metagenomic, and bacterial genomic analyses in eight representative large-scale chicken farms in Guangdong, China. Our results showed that antibiotics and ARGs contaminations were more severe in broiler farms than in layer farms. Notably, diverse tet(X) variants were prevalent in the chicken farms. These tet(X)s was carried by diverse E. coli lineages and obviously correlated with ISCR2 and IS1B transposases. The resistomes in chicken farms was significantly correlated with microbial community, and multiple factor analyses indicated that the joint effect of antibiotics-microbial community-MGEs was the most dominant driver of ARGs. Host tracking identified a variety of ARG bacterial hosts and the co-occurrence of ARGs-MRGs-MGEs. Source tracking indicated that the inherent component represented the main feature of resistomes in different hosts, while ARG transfer between the chicken gut and farm environments were frequent. A multiperspective evaluation of AR risk revealed that the early effect of antibiotic reduction was exhibited by the mitigation of maximum level of risky ARGs, prevalence of environmental AR pathogens, and HGT potential of ARGs mediated by phage structures. Overall, our findings provide insights into the antibiotic and ARG profiles in large-scale chicken farms with different rearing strategies and demonstrate a preliminary view of the performance of antibiotic reduction actions in China.

The dynamic history of prokaryotic phyla: discovery, diversity and division

Here, I review the dynamic history of prokaryotic phyla. Following leads set by Darwin, Haeckel and Woese, the concept of phylum has evolved from a group sharing common phenotypes to a set of organisms sharing a common ancestry, with modern taxonomy based on phylogenetic classifications drawn from macromolecular sequences. Phyla came as surprising latecomers to the formalities of prokaryotic nomenclature in 2021. Since then names have been validly published for 46 prokaryotic phyla, replacing some established names with neologisms, prompting criticism and debate within the scientific community. Molecular barcoding enabled phylogenetic analysis of microbial ecosystems without cultivation, leading to the identification of candidate divisions (or phyla) from diverse environments. The introduction of metagenome-assembled genomes marked a significant advance in identifying and classifying uncultured microbial phyla. The lumper-splitter dichotomy has led to disagreements, with experts cautioning against the pressure to create a profusion of new phyla and prominent databases adopting a conservative stance. The Candidatus designation has been widely used to provide provisional status to uncultured prokaryotic taxa, with phyla named under this convention now clearly surpassing those with validly published names. The Genome Taxonomy Database (GTDB) has offered a stable, standardized prokaryotic taxonomy with normalized taxonomic ranks, which has led to both lumping and splitting of pre-existing phyla. The GTDB framework introduced unwieldy alphanumeric placeholder labels, prompting recent publication of over 100 user-friendly Latinate names for unnamed prokaryotic phyla. Most candidate phyla remain 'known unknowns', with limited knowledge of their genomic diversity, ecological roles, or environments. Whether phyla still reflect significant evolutionary and ecological partitions across prokaryotic life remains an area of active debate. However, phyla remain of practical importance for microbiome analyses, particularly in clinical research. Despite potential diminishing returns in discovery of biodiversity, prokaryotic phyla offer extensive research opportunities for microbiologists for the foreseeable future.

The Saint-Leonard Urban Glaciotectonic Cave Harbors Rich and Diverse Planktonic and Sedimentary Microbial Communities

The terrestrial subsurface harbors unique microbial communities that play important biogeochemical roles and allow for studying a yet unknown fraction of the Earth's biodiversity. The Saint-Leonard cave in Montreal City (Canada) is of glaciotectonic origin. Its speleogenesis traces back to the withdrawal of the Laurentide Ice Sheet 13,000 years ago, during which the moving glacier dislocated the sedimentary rock layers. Our study is the first to investigate the microbial communities of the Saint-Leonard cave. By using amplicon sequencing, we analyzed the taxonomic diversity and composition of bacterial, archaeal and eukaryote communities living in the groundwater (0.1 µm- and 0.2 µm-filtered water), in the sediments and in surface soils. We identified a microbial biodiversity typical of cave ecosystems. Communities were mainly shaped by habitat type and harbored taxa associated with a wide variety of lifestyles and metabolic capacities. Although we found evidence of a geochemical connection between the above soils and the cave's galleries, our results suggest that the community assembly dynamics are driven by habitat selection rather than dispersal. Furthermore, we found that the cave's groundwater, in addition to being generally richer in microbial taxa than sediments, contained a considerable diversity of ultra-small bacteria and archaea.

Globally distributed marine Gemmatimonadota have unique genomic potentials

BACKGROUND

Gemmatimonadota bacteria are widely distributed in nature, but their metabolic potential and ecological roles in marine environments are poorly understood.

RESULTS

Here, we obtained 495 metagenome-assembled genomes (MAGs), and associated viruses, from coastal to deep-sea sediments around the world. We used this expanded genomic catalog to compare the protein composition and update the phylogeny of these bacteria. The marine Gemmatimonadota are phylogenetically different from those previously reported from terrestrial environments. Functional analyses of these genomes revealed these marine genotypes are capable of degradation of complex organic carbon, denitrification, sulfate reduction, and oxidizing sulfide and sulfite. Interestingly, there is widespread genetic potential for secondary metabolite biosynthesis across Gemmatimonadota, which may represent an unexplored source of novel natural products. Furthermore, viruses associated with Gemmatimonadota have the potential to "hijack" and manipulate host metabolism, including the assembly of the lipopolysaccharide in their hosts.

CONCLUSIONS

This expanded genomic diversity advances our understanding of these globally distributed bacteria across a variety of ecosystems and reveals genetic distinctions between those in terrestrial and marine communities. Video Abstract.

Unveiling the unknown viral world in groundwater

Viruses as the prevailing biological entities are poorly understood in underground realms. Here, we establish the first metagenomic Groundwater Virome Catalogue (GWVC) comprising 280,420 viral species ( ≥ 5 kb) detected from 607 monitored wells in seven geo-environmental zones throughout China. In expanding ~10-fold the global portfolio of known groundwater viruses, we uncover over 99% novel viruses and about 95% novel viral clusters. By linking viruses to hosts from 119 prokaryotic phyla, we double the number of microbial phyla known to be virus-infected in groundwater. As keystone ultrasmall symbionts in aquifers, CPR bacteria and DPANN archaea are susceptible to virulent viruses. Certain complete CPR viruses even likely infect non-CPR bacteria, while partial CPR/DPANN viruses harbor cell-surface modification genes that assist symbiont cell adhesion to free-living microbes. This study reveals the unknown viral world and auxiliary metabolism associated with methane, nitrogen, sulfur, and phosphorus cycling in groundwater, and highlights the importance of subsurface virosphere in viral ecology.

Genome streamlining in Parcubacteria transitioning from soil to groundwater

BACKGROUND

To better understand the influence of habitat on the genetic content of bacteria, with a focus on members of Candidate Phyla Radiation (CPR) bacteria, we studied the effects of transitioning from soil via seepage waters to groundwater on genomic composition of ultra-small Parcubacteria, the dominating CPR class in seepage waters, using genome resolved metagenomics.

RESULTS

Bacterial metagenome-assembled genomes (MAGs), (318 total, 32 of Parcubacteria) were generated from seepage waters and compared directly to groundwater counterparts. The estimated average genome sizes of members of major phyla Proteobacteria, Bacteroidota and Cand. Patescibacteria (Candidate Phyla Radiation - CPR bacteria) were significantly higher in soil-seepage water as compared to their groundwater counterparts. Seepage water Parcubacteria (Paceibacteria) exhibited 1.18-fold greater mean genome size and 2-fold lower mean proportion of pseudogenes than those in groundwater. Bacteroidota and Proteobacteria also showed a similar trend of reduced genomes in groundwater compared to seepage. While exploring gene loss and adaptive gains in closely related CPR lineages in groundwater, we identified a membrane protein, and a lipoglycopeptide resistance gene unique to a seepage Parcubacterium genome. A nitrite reductase gene was also identified and was unique to the groundwater Parcubacteria genomes, likely acquired from other planktonic microbes via horizontal gene transfer.

CONCLUSIONS

Overall, our data suggest that bacteria in seepage waters, including ultra-small Parcubacteria, have significantly larger genomes and higher metabolic enrichment than their groundwater counterparts, highlighting possible genome streamlining of the latter in response to habitat selection in an oligotrophic environment.

Minimal and hybrid hydrogenases are active from archaea

Microbial hydrogen (H2) cycling underpins the diversity and functionality of diverse anoxic ecosystems. Among the three evolutionarily distinct hydrogenase superfamilies responsible, [FeFe] hydrogenases were thought to be restricted to bacteria and eukaryotes. Here, we show that anaerobic archaea encode diverse, active, and ancient lineages of [FeFe] hydrogenases through combining analysis of existing and new genomes with extensive biochemical experiments. [FeFe] hydrogenases are encoded by genomes of nine archaeal phyla and expressed by H2-producing Asgard archaeon cultures. We report an ultraminimal hydrogenase in DPANN archaea that binds the catalytic H-cluster and produces H2. Moreover, we identify and characterize remarkable hybrid complexes formed through the fusion of [FeFe] and [NiFe] hydrogenases in ten other archaeal orders. Phylogenetic analysis and structural modeling suggest a deep evolutionary history of hybrid hydrogenases. These findings reveal new metabolic adaptations of archaea, streamlined H2 catalysts for biotechnological development, and a surprisingly intertwined evolutionary history between the two major H2-metabolizing enzymes.

Contamination Status of Novel Organophosphate Esters Derived from Organophosphite Antioxidants in Soil and the Effects on Soil Bacterial Communities

The contamination status of novel organophosphate esters (NOPEs) and their precursors organophosphite antioxidants (OPAs) and hydroxylated/diester transformation products (OH-OPEs/di-OPEs) in soils across a large-scale area in China were investigated. The total concentrations of the three test NOPEs in soil were 82.4-716 ng g-1, which were considerably higher than those of traditional OPEs (4.50-430 ng g-1), OPAs (n.d.-30.8 ng g-1), OH-OPEs (n.d.-0.49 ng g-1), and di-OPEs (0.57-21.1 ng g-1). One NOPE compound, i.e., tris(2,4-di-tert-butylphenyl) phosphate (AO168 = O) contributed over 65% of the concentrations of the studied OPE-associated contaminants. A 30-day soil incubation experiment was performed to confirm the influence of AO168 = O on soil bacterial communities. Specific genera belonging to Proteobacteria, such as Lysobacter and Ensifer, were enriched in AO168 = O-contaminated soils. Moreover, the ecological function of methylotrophy was observed to be significantly enhanced (t-test, p < 0.01) in soil treated with AO168 = O, while nitrogen fixation was significantly inhibited (t-test, p < 0.01). These findings comprehensively revealed the contamination status of OPE-associated contaminants in the soil environment and provided the first evidence of the effects of NOPEs on soil microbial communities.

Comparative proteomics of a versatile, marine, iron-oxidizing chemolithoautotroph

This study conducted a comparative proteomic analysis to identify potential genetic markers for the biological function of chemolithoautotrophic iron oxidation in the marine bacterium Ghiorsea bivora. To date, this is the only characterized species in the class Zetaproteobacteria that is not an obligate iron-oxidizer, providing a unique opportunity to investigate differential protein expression to identify key genes involved in iron-oxidation at circumneutral pH. Over 1000 proteins were identified under both iron- and hydrogen-oxidizing conditions, with differentially expressed proteins found in both treatments. Notably, a gene cluster upregulated during iron oxidation was identified. This cluster contains genes encoding for cytochromes that share sequence similarity with the known iron-oxidase, Cyc2. Interestingly, these cytochromes, conserved in both Bacteria and Archaea, do not exhibit the typical β-barrel structure of Cyc2. This cluster potentially encodes a biological nanowire-like transmembrane complex containing multiple redox proteins spanning the inner membrane, periplasm, outer membrane, and extracellular space. The upregulation of key genes associated with this complex during iron-oxidizing conditions was confirmed by quantitative reverse transcription-PCR. These findings were further supported by electromicrobiological methods, which demonstrated negative current production by G. bivora in a three-electrode system poised at a cathodic potential. This research provides significant insights into the biological function of chemolithoautotrophic iron oxidation.

Genome-resolved metagenomics reveals diverse taxa and metabolic complexity in Antarctic lake microbial structures

Lake Untersee, a lake in Antarctica that is perennially covered with ice, is home to unique microbial structures that are not lithified. We have evaluated the structure of the community and its metabolic potential across the pigmented upper layers and the sediment-enriched deeper layers in these pinnacle and cone-shaped microbial structures using metagenomics. These microbial structures are inhabited by distinct communities. The upper layers of the cone-shaped structures have a higher abundance of the cyanobacterial MAG Microcoleus, while the pinnacle-shaped structures have a higher abundance of Elainellacea MAG. This suggests that cyanobacteria influence the morphologies of the mats. We identified stark contrasts in the composition of the community and its metabolic potential between the upper and lower layers of the mat. The upper layers of the mat, which receive light, have an increased abundance of photosynthetic pathways. In contrast, the lower layer has an increased abundance of heterotrophic pathways. Our results also showed that Lake Untersee is the first Antarctic lake with a substantial presence of ammonia-oxidizing Nitrospiracea and amoA genes. The genomic capacity for recycling biological molecules was prevalent across metagenome-assembled genomes (MAGs) that cover 19 phyla. This highlights the importance of nutrient scavenging in ultra-oligotrophic environments. Overall, our study provides new insights into the formation of microbial structures and the potential metabolic complexity of Antarctic laminated microbial mats. These mats are important environments for biodiversity that drives biogeochemical cycling in polar deserts.

Analysis of nearly 3000 archaeal genomes from terrestrial geothermal springs sheds light on interconnected biogeochemical processes

Terrestrial geothermal springs are physicochemically diverse and host abundant populations of Archaea. However, the diversity, functionality, and geological influences of these Archaea are not well understood. Here we explore the genomic diversity of Archaea in 152 metagenomes from 48 geothermal springs in Tengchong, China, collected from 2016 to 2021. Our dataset is comprised of 2949 archaeal metagenome-assembled genomes spanning 12 phyla and 392 newly identified species, which increases the known species diversity of Archaea by ~48.6%. The structures and potential functions of the archaeal communities are strongly influenced by temperature and pH, with high-temperature acidic and alkaline springs favoring archaeal abundance over Bacteria. Genome-resolved metagenomics and metatranscriptomics provide insights into the potential ecological niches of these Archaea and their potential roles in carbon, sulfur, nitrogen, and hydrogen metabolism. Furthermore, our findings illustrate the interplay of competition and cooperation among Archaea in biogeochemical cycles, possibly arising from overlapping functional niches and metabolic handoffs. Taken together, our study expands the genomic diversity of Archaea inhabiting geothermal springs and provides a foundation for more incisive study of biogeochemical processes mediated by Archaea in geothermal ecosystems.

Metagenomics of two aquifers with thermal anomalies in Mallorca Island, and proposal of new uncultivated taxa named following the rules of SeqCode

Groundwater offers an intriguing blend of distinctive physical and chemical conditions, constituting a challenge for microbial life. In Mallorca, the largest island of Balearic archipelago, harbours a variety of thermal anomalies (i.e., geothermal manifestation where surface aquifers exhibiting temperatures exceeding the regional average). The metagenomes of two aquifers in the centre and southern of the island showed Pseudomonadota to be the most represented phylum when using extracted 16S rRNA gene sequences. However, the microbial structures within and between aquifers were remarkably diverse but similar in their metabolic profiles as revealed by the metagenome-assembled genomes (MAGs) pointing to a prevalence of aerobic chemolithoautotrophic and heterotrophic metabolisms, especially in the Llucmajor aquifer. Also, some evidences of anaerobic lifestyles were detected, which would indicate that these environments either could suffer episodes of oxygen depletion or the anaerobes had been transported from deeper waters. We believe that the local environmental factors (temperature, external inputs or chemistry) seem to be more relevant than the connection and, eventually, transport of microbial cells within the aquifer in determining the highly divergent structures. Notably, most of the reconstructed genomes belonged to undescribed bacterial lineages and from them two high-quality MAGs could be classified as novel taxa named following the rules of the Code for Nomenclature of Prokaryotes Described from Sequence Data (SeqCode). Accordingly, we propose the new species and genus Costitxia debesea gen. nov., sp. nov., affiliated with the novel family Costitxiaceae fam. nov., order Costitxiales ord. nov. and class Costitxiia class. nov.; and the new new species and genus Lloretia debesea gen. nov. sp. nov. affiliated with the novel family Lloretiaceae fam. nov.

Experimental evidence of d-glutamate racemase activity in the uncultivated bacterium Candidatus Saccharimonas aalborgensis

The Candidate Phyla Radiation (CPR) encompasses widespread uncultivated bacteria with reduced genomes and limited metabolic capacities. Most CPR bacteria lack the minimal set of enzymes required for peptidoglycan (PG) synthesis, leaving it unclear how these bacteria produce this essential envelope component. In this study, we analysed the distribution of d-amino acid racemases that produce the universal PG components d-glutamate (d-Glu) or d-alanine (d-Ala). We also examined moonlighting enzymes that synthesize d-Glu or d-Ala. Unlike other phyla in the domain Bacteria, CPR bacteria do not exhibit these moonlighting activities and have, at most, one gene encoding either a Glu or Ala racemase. One of these 'orphan' racemases is a predicted Glu racemase (MurICPR) from the CPR bacterium Candidatus Saccharimonas aalborgenesis. The expression of MurICPR restores the growth of a Salmonella d-Glu auxotroph lacking its endogenous racemase and results in the substitution of l-Ala by serine as the first residue in a fraction of the PG stem peptides. In vitro, MurICPR exclusively racemizes Glu as a substrate. Therefore, Ca. Saccharimonas aalborgensis may couple Glu racemization to serine and d-Glu incorporation into the stem peptide. Our findings provide the first insights into the synthesis of PG by an uncultivated environmental bacterium and illustrate how to experimentally test enzymatic activities from CPR bacteria related to PG metabolism.

CyanoStrainChip: A Novel DNA Microarray Tool for High-Throughput Detection of Environmental Cyanobacteria at the Strain Level

Detecting cyanobacteria in environments is an important concern due to their crucial roles in ecosystems, and they can form blooms with the potential to harm humans and nonhuman entities. However, the most widely used methods for high-throughput detection of environmental cyanobacteria, such as 16S rRNA sequencing, typically provide above-species-level resolution, thereby disregarding intraspecific variation. To address this, we developed a novel DNA microarray tool, termed the CyanoStrainChip, that enables strain-level comprehensive profiling of environmental cyanobacteria. The CyanoStrainChip was designed to target 1277 strains; nearly all major groups of cyanobacteria are included by implementing 43,666 genome-wide, strain-specific probes. It demonstrated strong specificity by in vitro mock community experiments. The high correlation (Pearson's R > 0.97) between probe fluorescence intensities and the corresponding DNA amounts (ranging from 1-100 ng) indicated excellent quantitative capability. Consistent cyanobacterial profiles of field samples were observed by both the CyanoStrainChip and next-generation sequencing methods. Furthermore, CyanoStrainChip analysis of surface water samples in Lake Chaohu uncovered a high intraspecific variation of abundance change within the genus Microcystis between different severity levels of cyanobacterial blooms, highlighting two toxic Microcystis strains that are of critical concern for Lake Chaohu harmful blooms suppression. Overall, these results suggest a potential for CyanoStrainChip as a valuable tool for cyanobacterial ecological research and harmful bloom monitoring to supplement existing techniques.

Uncultivated DPANN archaea are ubiquitous inhabitants of global oxygen-deficient zones with diverse metabolic potential

Archaea belonging to the DPANN (Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota) superphylum have been found in an expanding number of environments and perform a variety of biogeochemical roles, including contributing to carbon, sulfur, and nitrogen cycling. Generally characterized by ultrasmall cell sizes and reduced genomes, DPANN archaea may form mutualistic, commensal, or parasitic interactions with various archaeal and bacterial hosts, influencing the ecology and functioning of microbial communities. While DPANN archaea reportedly comprise a sizeable fraction of the archaeal community within marine oxygen-deficient zone (ODZ) water columns, little is known about their metabolic capabilities in these ecosystems. We report 33 novel metagenome-assembled genomes (MAGs) belonging to the DPANN phyla Nanoarchaeota, Pacearchaeota, Woesearchaeota, Undinarchaeota, Iainarchaeota, and SpSt-1190 from pelagic ODZs in the Eastern Tropical North Pacific and the Arabian Sea. We find these archaea to be permanent, stable residents of all three major ODZs only within anoxic depths, comprising up to 1% of the total microbial community and up to 25%-50% of archaea as estimated from read mapping to MAGs. ODZ DPANN appear to be capable of diverse metabolic functions, including fermentation, organic carbon scavenging, and the cycling of sulfur, hydrogen, and methane. Within a majority of ODZ DPANN, we identify a gene homologous to nitrous oxide reductase. Modeling analyses indicate the feasibility of a nitrous oxide reduction metabolism for host-attached symbionts, and the small genome sizes and reduced metabolic capabilities of most DPANN MAGs suggest host-associated lifestyles within ODZs.

IMPORTANCE

Archaea from the DPANN (Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota) superphylum have diverse metabolic capabilities and participate in multiple biogeochemical cycles. While metagenomics and enrichments have revealed that many DPANN are characterized by ultrasmall genomes, few biosynthetic genes, and episymbiotic lifestyles, much remains unknown about their biology. We report 33 new DPANN metagenome-assembled genomes originating from the three global marine oxygen-deficient zones (ODZs), the first from these regions. We survey DPANN abundance and distribution within the ODZ water column, investigate their biosynthetic capabilities, and report potential roles in the cycling of organic carbon, methane, and nitrogen. We test the hypothesis that nitrous oxide reductases found within several ODZ DPANN genomes may enable ultrasmall episymbionts to serve as nitrous oxide consumers when attached to a host nitrous oxide producer. Our results indicate DPANN archaea as ubiquitous residents within the anoxic core of ODZs with the potential to produce or consume key compounds.

BASALT refines binning from metagenomic data and increases resolution of genome-resolved metagenomic analysis

Metagenomic binning is an essential technique for genome-resolved characterization of uncultured microorganisms in various ecosystems but hampered by the low efficiency of binning tools in adequately recovering metagenome-assembled genomes (MAGs). Here, we introduce BASALT (Binning Across a Series of Assemblies Toolkit) for binning and refinement of short- and long-read sequencing data. BASALT employs multiple binners with multiple thresholds to produce initial bins, then utilizes neural networks to identify core sequences to remove redundant bins and refine non-redundant bins. Using the same assemblies generated from Critical Assessment of Metagenome Interpretation (CAMI) datasets, BASALT produces up to twice as many MAGs as VAMB, DASTool, or metaWRAP. Processing assemblies from a lake sediment dataset, BASALT produces ~30% more MAGs than metaWRAP, including 21 unique class-level prokaryotic lineages. Functional annotations reveal that BASALT can retrieve 47.6% more non-redundant opening-reading frames than metaWRAP. These results highlight the robust handling of metagenomic sequencing data of BASALT.

Peeling off the layers from microbial dark matter (MDM): recent advances, future challenges, and opportunities

Microbes represent the most common organisms on Earth; however, less than 2% of microbial species in the environment can undergo cultivation for study under laboratory conditions, and the rest of the enigmatic, microbial world remains mysterious, constituting a kind of "microbial dark matter" (MDM). In the last two decades, remarkable progress has been made in culture-dependent and culture-independent techniques. More recently, studies of MDM have relied on culture-independent techniques to recover genetic material through either unicellular genomics or shotgun metagenomics to construct single-amplified genomes (SAGs) and metagenome-assembled genomes (MAGs), respectively, which provide information about evolution and metabolism. Despite the remarkable progress made in the past decades, the functional diversity of MDM still remains uncharacterized. This review comprehensively summarizes the recently developed culture-dependent and culture-independent techniques for characterizing MDM, discussing major challenges, opportunities, and potential applications. These activities contribute to expanding our knowledge of the microbial world and have implications for various fields including Biotechnology, Bioprospecting, Functional genomics, Medicine, Evolutionary and Planetary biology. Overall, this review aims to peel off the layers from MDM, shed light on recent advancements, identify future challenges, and illuminate the exciting opportunities that lie ahead in unraveling the secrets of this intriguing microbial realm.

Cutting edge tools in the field of soil microbiology

The study of the whole of the genetic material contained within the microbial populations found in a certain environment is made possible by metagenomics. This technique enables a thorough knowledge of the variety, function, and interactions of microbial communities that are notoriously difficult to research. Due to the limitations of conventional techniques such as culturing and PCR-based methodologies, soil microbiology is a particularly challenging field. Metagenomics has emerged as an effective technique for overcoming these obstacles and shedding light on the dynamic nature of the microbial communities in soil. This review focuses on the principle of metagenomics techniques, their potential applications and limitations in soil microbial diversity analysis. The effectiveness of target-based metagenomics in determining the function of individual genes and microorganisms in soil ecosystems is also highlighted. Targeted metagenomics, including high-throughput sequencing and stable-isotope probing, is essential for studying microbial taxa and genes in complex ecosystems. Shotgun metagenomics may reveal the diversity of soil bacteria, composition, and function impacted by land use and soil management. Sanger, Next Generation Sequencing, Illumina, and Ion Torrent sequencing revolutionise soil microbiome research. Oxford Nanopore Technology (ONT) and Pacific Biosciences (PacBio)'s third and fourth generation sequencing systems revolutionise long-read technology. GeoChip, clone libraries, metagenomics, and metabarcoding help comprehend soil microbial communities. The article indicates that metagenomics may improve environmental management and agriculture despite existing limitations.Metagenomics has revolutionised soil microbiology research by revealing the complete diversity, function, and interactions of microorganisms in soil. Metagenomics is anticipated to continue defining the future of soil microbiology research despite some limitations, such as the difficulty of locating the appropriate sequencing method for specific genes.

Autotrophic biofilms sustained by deeply sourced groundwater host diverse bacteria implicated in sulfur and hydrogen metabolism

BACKGROUND

Biofilms in sulfide-rich springs present intricate microbial communities that play pivotal roles in biogeochemical cycling. We studied chemoautotrophically based biofilms that host diverse CPR bacteria and grow in sulfide-rich springs to investigate microbial controls on biogeochemical cycling.

RESULTS

Sulfide springs biofilms were investigated using bulk geochemical analysis, genome-resolved metagenomics, and scanning transmission X-ray microscopy (STXM) at room temperature and 87 K. Chemolithotrophic sulfur-oxidizing bacteria, including Thiothrix and Beggiatoa, dominate the biofilms, which also contain CPR Gracilibacteria, Absconditabacteria, Saccharibacteria, Peregrinibacteria, Berkelbacteria, Microgenomates, and Parcubacteria. STXM imaging revealed ultra-small cells near the surfaces of filamentous bacteria that may be CPR bacterial episymbionts. STXM and NEXAFS spectroscopy at carbon K and sulfur L2,3 edges show that filamentous bacteria contain protein-encapsulated spherical elemental sulfur granules, indicating that they are sulfur oxidizers, likely Thiothrix. Berkelbacteria and Moranbacteria in the same biofilm sample are predicted to have a novel electron bifurcating group 3b [NiFe]-hydrogenase, putatively a sulfhydrogenase, potentially linked to sulfur metabolism via redox cofactors. This complex could potentially contribute to symbioses, for example, with sulfur-oxidizing bacteria such as Thiothrix that is based on cryptic sulfur cycling. One Doudnabacteria genome encodes adjacent sulfur dioxygenase and rhodanese genes that may convert thiosulfate to sulfite. We find similar conserved genomic architecture associated with CPR bacteria from other sulfur-rich subsurface ecosystems.

CONCLUSIONS

Our combined metagenomic, geochemical, spectromicroscopic, and structural bioinformatics analyses of biofilms growing in sulfide-rich springs revealed consortia that contain CPR bacteria and sulfur-oxidizing Proteobacteria, including Thiothrix, and bacteria from a new family within Beggiatoales. We infer roles for CPR bacteria in sulfur and hydrogen cycling. Video Abstract.

Genesis of ectosymbiotic features based on commensalistic syntrophy

The symbiogenetic origin of eukaryotes with mitochondria is considered a major evolutionary transition. The initial interactions and conditions of symbiosis, along with the phylogenetic affinity of the host, are widely debated. Here, we focus on a possible evolutionary path toward an association of individuals of two species based on unidirectional syntrophy. With the backing of a theoretical model, we hypothesize that the first step in the evolution of such symbiosis could be the appearance of a linking structure on the symbiont's membrane, using which it forms an ectocommensalism with its host. We consider a commensalistic model based on the syntrophy hypothesis in the framework of coevolutionary dynamics and mutant invasion into a monomorphic resident system (evolutionary substitution). We investigate the ecological and evolutionary stability of the consortium (or symbiotic merger), with vertical transmissions playing a crucial role. The impact of the 'effectiveness of vertical transmission' on the dynamics is also analyzed. We find that the transmission of symbionts and the additional costs incurred by the mutant determine the conditions of fixation of the consortia. Additionally, we observe that small and highly metabolically active symbionts are likely to form the consortia.

High niche specificity and host genetic diversity of groundwater viruses

Viruses are key members of microbial communities that exert control over host abundance and metabolism, thereby influencing ecosystem processes and biogeochemical cycles. Aquifers are known to host taxonomically diverse microbial life, yet little is known about viruses infecting groundwater microbial communities. Here, we analysed 16 metagenomes from a broad range of groundwater physicochemistries. We recovered 1571 viral genomes that clustered into 468 high-quality viral operational taxonomic units. At least 15% were observed to be transcriptionally active, although lysis was likely constrained by the resource-limited groundwater environment. Most were unclassified (95%), and the remaining 5% were Caudoviricetes. Comparisons with viruses inhabiting other aquifers revealed no shared species, indicating substantial unexplored viral diversity. In silico predictions linked 22.4% of the viruses to microbial host populations, including to ultra-small prokaryotes, such as Patescibacteria and Nanoarchaeota. Many predicted hosts were associated with the biogeochemical cycling of carbon, nitrogen, and sulfur. Metabolic predictions revealed the presence of 205 putative auxiliary metabolic genes, involved in diverse processes associated with the utilization of the host's intracellular resources for biosynthesis and transformation reactions, including those involved in nucleotide sugar, glycan, cofactor, and vitamin metabolism. Viruses, prokaryotes overall, and predicted prokaryotic hosts exhibited narrow spatial distributions, and relative abundance correlations with the same groundwater parameters (e.g. dissolved oxygen, nitrate, and iron), consistent with host control over viral distributions. Results provide insights into underexplored groundwater viruses, and indicate the large extent to which viruses may manipulate microbial communities and biogeochemistry in the terrestrial subsurface.

Large attachment organelle mediates interaction between Nanobdellota archaeon YN1 and its host

DPANN archaea are an enigmatic superphylum that are difficult to isolate and culture in the laboratory due to their specific culture conditions and apparent ectosymbiotic lifestyle. Here, we successfully isolated and cultivated a coculture system of a novel Nanobdellota archaeon YN1 and its host Sulfurisphaera ohwakuensis YN1HA. We characterized the coculture system by complementary methods, including metagenomics and metabolic pathway analysis, fluorescence microscopy, and high-resolution electron cryo-tomography (cryoET). We show that YN1 is deficient in essential metabolic processes and requires host resources to proliferate. CryoET imaging revealed an enormous attachment organelle present in the YN1 envelope that forms a direct interaction with the host cytoplasm, bridging the two cells. Together, our results unravel the molecular and structural basis of ectosymbiotic relationship between YN1 and YN1HA. This research broadens our understanding of DPANN biology and the versatile nature of their ectosymbiotic relationships.

Genome-resolved metaproteogenomic and nanosolid characterization of an inactive vent chimney densely colonized by enigmatic DPANN archaea

Recent successes in the cultivation of DPANN archaea with their hosts have demonstrated an episymbiotic lifestyle, whereas the lifestyle of DPANN archaea in natural habitats is largely unknown. A free-living lifestyle is speculated in oxygen-deprived fluids circulated through rock media, where apparent hosts of DPANN archaea are lacking. Alternatively, DPANN archaea may be detached from their hosts and/or rock surfaces. To understand the ecology of rock-hosted DPANN archaea, rocks rather than fluids should be directly characterized. Here, we investigated a deep-sea hydrothermal vent chimney without fluid venting where our previous study revealed the high proportion of Pacearchaeota, one of the widespread and enigmatic lineages of DPANN archaea. Using spectroscopic methods with submicron soft X-ray and infrared beams, the microbial habitat was specified to be silica-filled pores in the inner chimney wall comprising chalcopyrite. Metagenomic analysis of the inner wall revealed the lack of biosynthetic genes for nucleotides, amino acids, cofactors, and lipids in the Pacearchaeota genomes. Genome-resolved metaproteomic analysis clarified the co-occurrence of a novel thermophilic lineage actively fixing carbon and nitrogen and thermophilic archaea in the inner chimney wall. We infer that the shift in metabolically active microbial populations from the thermophiles to the mesophilic DPANN archaea occurs after the termination of fluid venting. The infilling of mineral pores by hydrothermal silica deposition might be a preferred environmental factor for the colonization of free-living Pacearchaeota with ultrasmall cells depending on metabolites synthesized by the co-occurring thermophiles during fluid venting.

Microbial ecology of the deep terrestrial subsurface

The terrestrial subsurface hosts microbial communities that, collectively, are predicted to comprise as many microbial cells as global surface soils. Although initially thought to be associated with deposited organic matter, deep subsurface microbial communities are supported by chemolithoautotrophic primary production, with hydrogen serving as an important source of electrons. Despite recent progress, relatively little is known about the deep terrestrial subsurface compared to more commonly studied environments. Understanding the composition of deep terrestrial subsurface microbial communities and the factors that influence them is of importance because of human-associated activities including long-term storage of used nuclear fuel, carbon capture, and storage of hydrogen for use as an energy vector. In addition to identifying deep subsurface microorganisms, recent research focuses on identifying the roles of microorganisms in subsurface communities, as well as elucidating myriad interactions-syntrophic, episymbiotic, and viral-that occur among community members. In recent years, entirely new groups of microorganisms (i.e. candidate phyla radiation bacteria and Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoloarchaeota, Nanoarchaeota archaea) have been discovered in deep terrestrial subsurface environments, suggesting that much remains unknown about this biosphere. This review explores the historical context for deep terrestrial subsurface microbial ecology and highlights recent discoveries that shape current ecological understanding of this poorly explored microbial habitat. Additionally, we highlight the need for multifaceted experimental approaches to observe phenomena such as cryptic cycles, complex interactions, and episymbiosis, which may not be apparent when using single approaches in isolation, but are nonetheless critical to advancing our understanding of this deep biosphere.

Implications for nitrogen and sulphur cycles: phylogeny and niche-range of Nitrospirota in terrestrial aquifers

Increasing evidence suggests Nitrospirota are important contributors to aquatic and subsurface nitrogen and sulphur cycles. We determined the phylogenetic and ecological niche associations of Nitrospirota colonizing terrestrial aquifers. Nitrospirota compositions were determined across 59 groundwater wells. Distributions were strongly influenced by oxygen availability in groundwater, marked by a trade-off between aerobic (Nitrospira, Leptospirillum) and anaerobic (Thermodesulfovibrionia, unclassified) lineages. Seven Nitrospirota metagenome-assembled genomes (MAGs), or populations, were recovered from a subset of wells, including three from the recently designated class 9FT-COMBO-42-15. Most were relatively more abundant and transcriptionally active in dysoxic groundwater. These MAGs were analysed with 743 other Nitrospirota genomes. Results illustrate the predominance of certain lineages in aquifers (e.g. non-nitrifying Nitrospiria, classes 9FT-COMBO-42-15 and UBA9217, and Thermodesulfovibrionales family UBA1546). These lineages are characterized by mechanisms for nitrate reduction and sulphur cycling, and, excluding Nitrospiria, the Wood-Ljungdahl pathway, consistent with carbon-limited, low-oxygen, and sulphur-rich aquifer conditions. Class 9FT-COMBO-42-15 is a sister clade of Nitrospiria and comprises two families spanning a transition in carbon fixation approaches: f_HDB-SIOIB13 encodes rTCA (like Nitrospiria) and f_9FT-COMBO-42-15 encodes Wood-Ljungdahl CO dehydrogenase (like Thermodesulfovibrionia and UBA9217). The 9FT-COMBO-42-15 family is further differentiated by its capacity for sulphur oxidation (via DsrABEFH and SoxXAYZB) and dissimilatory nitrate reduction to ammonium, and gene transcription indicated active coupling of nitrogen and sulphur cycles by f_9FT-COMBO-42-15 in dysoxic groundwater. Overall, results indicate that Nitrospirota are widely distributed in groundwater and that oxygen availability drives the spatial differentiation of lineages with ecologically distinct roles related to nitrogen and sulphur metabolism.

Enrichment of different taxa of the enigmatic candidate phyla radiation bacteria using a novel picolitre droplet technique

The candidate phyla radiation (CPR) represents a distinct monophyletic clade and constitutes a major portion of the tree of life. Extensive efforts have focused on deciphering the functional diversity of its members, primarily using sequencing-based techniques. However, cultivation success remains scarce, presenting a significant challenge, particularly in CPR-dominated groundwater microbiomes characterized by low biomass. Here, we employ an advanced high-throughput droplet microfluidics technique to enrich CPR taxa from groundwater. Utilizing a low-volume filtration approach, we successfully harvested a microbiome resembling the original groundwater microbial community. We assessed CPR enrichment in droplet and aqueous bulk cultivation for 30 days using a novel CPR-specific primer to rapidly track the CPR fraction through the cultivation attempts. The combination of soil extract and microbial-derived necromass provided the most supportive conditions for CPR enrichment. Employing these supplemented conditions, droplet cultivation proved superior to bulk cultivation, resulting in up to a 13-fold CPR enrichment compared to a 1- to 2-fold increase in bulk cultivation. Amplicon sequencing revealed 10 significantly enriched CPR orders. The highest enrichment in CPRs was observed for some unknown members of the Parcubacteria order, Cand. Jorgensenbacteria, and unclassified UBA9983. Furthermore, we identified co-enriched putative host taxa, which may guide more targeted CPR isolation approaches in subsequent investigations.

Viromes vs. mixed community metagenomes: choice of method dictates interpretation of viral community ecology

Background

Viruses, the majority of which are uncultivated, are among the most abundant biological entities on Earth. From altering microbial physiology to driving community dynamics, viruses are fundamental members of microbiomes. While the number of studies leveraging viral metagenomics (viromics) for studying uncultivated viruses is growing, standards for viromics research are lacking. Viromics can utilize computational discovery of viruses from total metagenomes of all community members (hereafter metagenomes) or use physical separation of virus-specific fractions (hereafter viromes). However, differences in the recovery and interpretation of viruses from metagenomes and viromes obtained from the same samples remain understudied.

Results

Here, we compare viral communities from paired viromes and metagenomes obtained from 60 diverse samples across human gut, soil, freshwater, and marine ecosystems. Overall, viral communities obtained from viromes were more abundant and species rich than those obtained from metagenomes, although there were some exceptions. Despite this, metagenomes still contained many viral genomes not detected in viromes. We also found notable differences in the predicted lytic state of viruses detected in viromes vs metagenomes at the time of sequencing. Other forms of variation observed include genome presence/absence, genome quality, and encoded protein content between viromes and metagenomes, but the magnitude of these differences varied by environment.

Conclusions

Overall, our results show that the choice of method can lead to differing interpretations of viral community ecology. We suggest that the choice of whether to target a metagenome or virome to study viral communities should be dependent on the environmental context and ecological questions being asked. However, our overall recommendation to researchers investigating viral ecology and evolution is to pair both approaches to maximize their respective benefits.

Phylogenetic and ecophysiological novelty of subsurface mercury methylators in mangrove sediments

Mangrove sediment is a crucial component in the global mercury (Hg) cycling and acts as a hotspot for methylmercury (MeHg) production. Early evidence has documented the ubiquity of well-studied Hg methylators in mangrove superficial sediments; however, their diversity and metabolic adaptation in the more anoxic and highly reduced subsurface sediments are lacking. Through MeHg biogeochemical assay and metagenomic sequencing, we found that mangrove subsurface sediments (20-100 cm) showed a less hgcA gene abundance but higher diversity of Hg methylators than superficial sediments (0-20 cm). Regional-scale investigation of mangrove subsurface sediments spanning over 1500 km demonstrated a prevalence and family-level novelty of Hg-methylating microbial lineages (i.e., those affiliated to Anaerolineae, Phycisphaerae, and Desulfobacterales). We proposed the candidate phylum Zixibacteria lineage with sulfate-reducing capacity as a currently understudied Hg methylator across anoxic environments. Unlike other Hg methylators, the Zixibacteria lineage does not use the Wood-Ljungdahl pathway but has unique capabilities of performing methionine synthesis to donate methyl groups. The absence of cobalamin biosynthesis pathway suggests that this Hg-methylating lineage may depend on its syntrophic partners (i.e., Syntrophobacterales members) for energy in subsurface sediments. Our results expand the diversity of subsurface Hg methylators and uncover their unique ecophysiological adaptations in mangrove sediments.

VBCG: 20 validated bacterial core genes for phylogenomic analysis with high fidelity and resolution

BACKGROUND

Phylogenomic analysis has become an inseparable part of studies of bacterial diversity and evolution, and many different bacterial core genes have been collated and used for phylogenomic tree reconstruction. However, these genes have been selected based on their presence and single-copy ratio in all bacterial genomes, leaving out the gene's 'phylogenetic fidelity' unexamined.

RESULTS

From 30,522 complete genomes covering 11,262 species, we examined 148 bacterial core genes that have been previously used for phylogenomic analysis. In addition to the gene presence and single-copy rations, we evaluated the gene's phylogenetic fidelity by comparing each gene's phylogeny with its corresponding 16S rRNA gene tree. Out of the 148 bacterial genes, 20 validated bacterial core genes (VBCG) were selected as the core gene set with the highest bacterial phylogenetic fidelity. Compared to the larger gene set, the 20-gene core set resulted in more species having all genes present and fewer species with missing data, thereby enhancing the accuracy of phylogenomic analysis. Using Escherichia coli strains as examples of prominent bacterial foodborne pathogens, we demonstrated that the 20 VBCG produced phylogenies with higher fidelity and resolution at species and strain levels while 16S rRNA gene tree alone could not.

CONCLUSION

The 20 validated core gene set improves the fidelity and speed of phylogenomic analysis. Among other uses, this tool improves our ability to explore the evolution, typing and tracking of bacterial strains, such as human pathogens. We have developed a Python pipeline and a desktop graphic app (available on GitHub) for users to perform phylogenomic analysis with high fidelity and resolution. Video Abstract.

Hidden from plain sight: Novel Simkaniaceae and Rhabdochlamydiaceae diversity emerging from screening genomic and metagenomic data

Chlamydiota are an ancient and hyperdiverse phylum of obligate intracellular bacteria. The best characterized representatives are pathogens or parasites of mammals, but it is thought that their most common hosts are microeukaryotes like Amoebozoa. The diversity in taxonomy, evolution, and function of non-pathogenic Chlamydiota are slowly being described. Here we use data mining techniques and genomic analysis to extend our current knowledge of Chlamydiota diversity and its hosts, in particular the Order Parachlamydiales. We extract one Rhabdochlamydiaceae and three Simkaniaceae Metagenome-Assembled Genomes (MAGs) from NCBI Short Read Archive deposits of ciliate and algal genome sequencing projects. We then use these to identify a further 14 and 8 MAGs respectively amongst existing, unidentified environmental assemblies. From these data we identify two novel clades with host associated data, for which we propose the names "Sacchlamyda saccharinae" (Family Rhabdochlamydiaceae) and "Amphrikana amoebophyrae" (Family Simkaniaceae), as well as a third new clade of environmental MAGs "Acheromyda pituitae" (Family Rhabdochlamydiaceae). The extent of uncharacterized diversity within the Rhabdochlamydiaceae and Simkaniaceae is indicated by 16 of the 22 MAGs being evolutionarily distant from currently characterised genera. Within our limited data, there was great predicted diversity in Parachlamydiales metabolism and evolution, including the potential for metabolic and defensive symbioses as well as pathogenicity. These data provide an imperative to link genomic diversity in metagenomics data to their associated eukaryotic host, and to develop onward understanding of the functional significance of symbiosis with this hyperdiverse clade.

A genome catalogue of lake bacterial diversity and its drivers at continental scale

Lakes are heterogeneous ecosystems inhabited by a rich microbiome whose genomic diversity is poorly defined. We present a continental-scale study of metagenomes representing 6.5 million km2 of the most lake-rich landscape on Earth. Analysis of 308 Canadian lakes resulted in a metagenome-assembled genome (MAG) catalogue of 1,008 mostly novel bacterial genomospecies. Lake trophic state was a leading driver of taxonomic and functional diversity among MAG assemblages, reflecting the responses of communities profiled by 16S rRNA amplicons and gene-centric metagenomics. Coupling the MAG catalogue with watershed geomatics revealed terrestrial influences of soils and land use on assemblages. Agriculture and human population density were drivers of turnover, indicating detectable anthropogenic imprints on lake bacteria at the continental scale. The sensitivity of bacterial assemblages to human impact reinforces lakes as sentinels of environmental change. Overall, the LakePulse MAG catalogue greatly expands the freshwater genomic landscape, advancing an integrative view of diversity across Earth's microbiomes.

Habitat Transition in the Evolution of Bacteria and Archaea

Related groups of microbes are widely distributed across Earth's habitats, implying numerous dispersal and adaptation events over evolutionary time. However, relatively little is known about the characteristics and mechanisms of these habitat transitions, particularly for populations that reside in animal microbiomes. Here, we review the literature concerning habitat transitions among a variety of bacterial and archaeal lineages, considering the frequency of migration events, potential environmental barriers, and mechanisms of adaptation to new physicochemical conditions, including the modification of protein inventories and other genomic characteristics. Cells dependent on microbial hosts, particularly bacteria from the Candidate Phyla Radiation, have undergone repeated habitat transitions from environmental sources into animal microbiomes. We compare their trajectories to those of both free-living cells-including the Melainabacteria, Elusimicrobia, and methanogenic archaea-and cellular endosymbionts and bacteriophages, which have made similar transitions. We conclude by highlighting major related topics that may be worthy of future study.

"Candidatus Subterrananammoxibiaceae," a New Anammox Bacterial Family in Globally Distributed Marine and Terrestrial Subsurfaces

Bacteria specialized in anaerobic ammonium oxidation (anammox) are widespread in many anoxic habitats and form an important functional guild in the global nitrogen cycle by consuming bio-available nitrogen for energy rather than biomass production. Due to their slow growth rates, cultivation-independent approaches have been used to decipher their diversity across environments. However, their full diversity has not been well recognized. Here, we report a new family of putative anammox bacteria, "Candidatus Subterrananammoxibiaceae," existing in the globally distributed terrestrial and marine subsurface (groundwater and sediments of estuary, deep-sea, and hadal trenches). We recovered a high-quality metagenome-assembled genome of this family, tentatively named "Candidatus Subterrananammoxibius californiae," from a California groundwater site. The "Ca. Subterrananammoxibius californiae" genome not only contains genes for all essential components of anammox metabolism (e.g., hydrazine synthase, hydrazine oxidoreductase, nitrite reductase, and nitrite oxidoreductase) but also has the capacity for urea hydrolysis. In an Arctic ridge sediment core where redox zonation is well resolved, "Ca. Subterrananammoxibiaceae" is confined within the nitrate-ammonium transition zone where the anammox rate maximum occurs, providing environmental proof of the anammox activity of this new family. Phylogenetic analysis of nitrite oxidoreductase suggests that a horizontal transfer facilitated the spreading of the nitrite oxidation capacity between anammox bacteria (in the Planctomycetota phylum) and nitrite-oxidizing bacteria from Nitrospirota and Nitrospinota. By recognizing this new anammox family, we propose that all lineages within the "Ca. Brocadiales" order have anammox capacity. IMPORTANCE Microorganisms called anammox bacteria are efficient in removing bioavailable nitrogen from many natural and human-made environments. They exist in almost every anoxic habitat where both ammonium and nitrate/nitrite are present. However, only a few anammox bacteria have been cultured in laboratory settings, and their full phylogenetic diversity has not been recognized. Here, we present a new bacterial family whose members are present across both the terrestrial and marine subsurface. By reconstructing a high-quality genome from the groundwater environment, we demonstrate that this family has all critical enzymes of anammox metabolism and, notably, also urea utilization. This bacterium family in marine sediments is also preferably present in the niche where the anammox process occurs. These findings suggest that this novel family, named "Candidatus Subterrananammoxibiaceae," is an overlooked group of anammox bacteria, which should have impacts on nitrogen cycling in a range of environments.

Identifying and tracking mobile elements in evolving compost communities yields insights into the nanobiome

Microbial evolution is driven by rapid changes in gene content mediated by horizontal gene transfer (HGT). While mobile genetic elements (MGEs) are important drivers of gene flux, the nanobiome-the zoo of Darwinian replicators that depend on microbial hosts-remains poorly characterised. New approaches are necessary to increase our understanding beyond MGEs shaping individual populations, towards their impacts on complex microbial communities. A bioinformatic pipeline (xenoseq) was developed to cross-compare metagenomic samples from microbial consortia evolving in parallel, aimed at identifying MGE dissemination, which was applied to compost communities which underwent periodic mixing of MGEs. We show that xenoseq can distinguish movement of MGEs from demographic changes in community composition that otherwise confounds identification, and furthermore demonstrate the discovery of various unexpected entities. Of particular interest was a nanobacterium of the candidate phylum radiation (CPR) which is closely related to a species identified in groundwater ecosystems (Candidatus Saccharibacterium), and appears to have a parasitic lifestyle. We also highlight another prolific mobile element, a 313 kb plasmid hosted by a Cellvibrio lineage. The host was predicted to be capable of nitrogen fixation, and acquisition of the plasmid coincides with increased ammonia production. Taken together, our data show that new experimental strategies combined with bioinformatic analyses of metagenomic data stand to provide insight into the nanobiome as a driver of microbial community evolution.

Global transmission of broad-host-range plasmids derived from the human gut microbiome

Broad-host-range (BHR) plasmids in human gut bacteria are of considerable interest for their ability to mediate horizontal gene transfer (HGT) across large phylogenetic distance. However, the human gut plasmids, especially the BHR plasmids, remain largely unknown. Here, we identified the plasmids in the draft genomes of gut bacterial isolates from Chinese and American donors, resulting in 5372 plasmid-like clusters (PLCs), of which, 820 PLCs (comPLCs) were estimated with > 60% completeness genomes and only 155 (18.9%) were classified to known replicon types (n = 37). We observed that 175 comPLCs had a broad host range across distinct bacterial genera, of which, 71 were detected in at least two human populations of Chinese, American, Spanish, and Danish, and 13 were highly prevalent (>10%) in at least one human population. Haplotype analyses of two widespread PLCs demonstrated their spreading and evolutionary trajectory, suggesting frequent and recent exchanges of the BHR plasmids in environments. In conclusion, we obtained a large collection of plasmid sequences in human gut bacteria and demonstrated that a subset of the BHR plasmids can be transmitted globally, thus facilitating extensive HGT (e.g. antibiotic resistance genes) events. This study highlights the potential implications of the plasmids for global human health.

Influence of soil nutrients on the presence and distribution of CPR bacteria in a long-term crop rotation experiment

Bacteria affiliated with the Candidate Phyla Radiation (CPR) are a hyper-diverse group of ultra-small bacteria with versatile yet sparse metabolisms. However, most insights into this group come from a surprisingly small number of environments, and recovery of CPR bacteria from soils has been hindered due to their extremely low abundance within complex microbial assemblages. In this study we enriched soil samples from 14 different soil fertility treatments for ultra-small (<0.45 μm) bacteria in order to study rare soil CPR. 42 samples were sequenced, enabling the reconstruction of 27 quality CPR metagenome-assembled genomes (MAGs) further classified as Parcubacteria/Paceibacteria, Saccharibacteria/Saccharimonadia and ABY1, in addition to representative genomes from Gemmatimonadetes, Dependentiae and Chlamydae phyla. These genomes were fully annotated and used to reconstruct the CPR community across all 14 plots. Additionally, for five of these plots, the entire microbiota was reconstructed using 16S amplification, showing that specific soil CPR may form symbiotic relationships with a varied and circumstantial range of hosts. Cullars CPR had a prevalence of enzymes predicted to degrade plant-derived carbohydrates, which suggests they have a role in plant biomass degradation. Parcubacteria appear to be more apt at microfauna necromass degradation. Cullars Saccharibacteria and a Parcubacteria group were shown to carry a possible aerotolerance mechanism coupled with potential for aerobic respiration, which appear to be a unique adaptation to the oxic soil environment. Reconstruction of CPR communities across treatment plots showed that they were not impacted by changes in nutrient levels or microbiota composition, being only impacted by extreme conditions, causing some CPR to dominate the community. These findings corroborate the understanding that soil-dwelling CPR bacteria have a very broad symbiont range and have metabolic capabilities associated to soil environments which allows them to scavenge resources and form resilient communities. The contributions of these microbial dark matter species to soil ecology and plant interactions will be of significant interest in future studies.

May microbial ecological baseline exist in continental groundwater?

BACKGROUND

Microbes constitute almost the entire biological community in subsurface groundwater and play an important role in ecological evolution and global biogeochemical cycles. Ecological baseline as a fundamental reference with less human interference has been investigated in surface ecosystems such as soils, rivers, and ocean, but the existence of groundwater microbial ecological baseline (GMEB) is still an open question so far.

RESULTS

Based on high-throughput sequencing information derived from national monitoring of 733 newly constructed wells, we find that bacterial communities in pristine groundwater exhibit a significant lateral diversity gradient and gradually approach the topsoil microbial latitudinal diversity gradient with decreasing burial depth of phreatic water. Among 74 phyla dominated by Proteobacteria in groundwater, Patescibacteria act as keystone taxa that harmonize microbes in shallower aquifers and accelerate decline in bacterial diversity with increasing well-depth. Decreasing habitat niche breadth with increasing well-depth suggests a general change in the relationship among key microbes from closer cooperation in shallow to stronger competition in deep groundwater. Unlike surface-water microbes, microbial communities in pristine groundwater are predominantly shaped by deterministic processes, potentially associated with nutrient sequestration under dark and anoxic environments in aquifers.

CONCLUSIONS

By unveiling the biogeographic patterns and mechanisms controlling the community assembly of microbes in pristine groundwater throughout China, we firstly confirm the existence of GMEB in shallower aquifers and propose Groundwater Microbial Community Index (GMCI) to evaluate anthropogenic impact, which highlights the importance of GMEB in groundwater water security and health diagnosis. Video Abstract.

The vertically-stratified resistomes in mangrove sediments was driven by the bacterial diversity

Early evidence has elucidated that the spread of antibiotic (ARGs) and metal resistance genes (MRGs) are mainly attributed to the selection pressure in human-influenced environments. However, whether and how biotic and abiotic factors mediate the distribution of ARGs and MRGs in mangrove sediments under natural sedimentation is largely unclear. Here, we profiled the abundance and diversity of ARGs and MRGs and their relationships with sedimental microbiomes in 0-100 cm mangrove sediments. Our results identified multidrug-resistance and multimetal-resistance as the most abundant ARG and MRG classes, and their abundances generally decreased with the sediment depth. Instead of abiotic factors such as nutrients and antibiotics, the bacterial diversity was significantly negatively correlated with the abundance and diversity of resistomes. Also, the majority of resistance classes (e.g., multidrug and arsenic) were carried by more diverse bacterial hosts in deep layers with low abundances of resistance genes. Together, our results indicated that bacterial diversity was the most important biotic factor driving the vertical profile of ARGs and MRGs in the mangrove sediment. Given that there is a foreseeable increasing human impact on natural environments, this study emphasizes the important role of biodiversity in driving the abundance and diversity of ARGs and MRGs.

"Candidatus Nealsonbacteria" Are Likely Biomass Recycling Ectosymbionts of Methanogenic Archaea in a Stable Benzene-Degrading Enrichment Culture

The Candidate Phyla Radiation (CPR), also referred to as superphylum Patescibacteria, is a very large group of bacteria with no pure culture representatives discovered by 16S rRNA sequencing or genome-resolved metagenomic analyses of environmental samples. Within the CPR, candidate phylum Parcubacteria, previously referred to as OD1, is prevalent in anoxic sediments and groundwater. Previously, we had identified a specific member of the Parcubacteria (referred to as DGGOD1a) as an important member of a methanogenic benzene-degrading consortium. Phylogenetic analyses herein place DGGOD1a within the clade "Candidatus Nealsonbacteria." Because of its persistence over many years, we hypothesized that "Ca. Nealsonbacteria" DGGOD1a must play an important role in sustaining anaerobic benzene metabolism in the consortium. To try to identify its growth substrate, we amended the culture with a variety of defined compounds (pyruvate, acetate, hydrogen, DNA, and phospholipid), as well as crude culture lysate and three subfractions thereof. We observed the greatest (10-fold) increase in the absolute abundance of "Ca. Nealsonbacteria" DGGOD1a only when the consortium was amended with crude cell lysate. These results implicate "Ca. Nealsonbacteria" in biomass recycling. Fluorescence in situ hybridization and cryogenic transmission electron microscope images revealed that "Ca. Nealsonbacteria" DGGOD1a cells were attached to larger archaeal Methanothrix cells. This apparent epibiont lifestyle was supported by metabolic predictions from a manually curated complete genome. This is one of the first examples of bacterial-archaeal episymbiosis and may be a feature of other "Ca. Nealsonbacteria" found in anoxic environments. IMPORTANCE An anaerobic microbial enrichment culture was used to study members of candidate phyla that are difficult to grow in the lab. We were able to visualize tiny "Candidatus Nealsonbacteria" cells attached to a large Methanothrix cell, revealing a novel episymbiosis.

N/S element transformation modulating lithospheric microbial communities by single-species manipulation

BACKGROUND

The lithospheric microbiome plays a vital role in global biogeochemical cycling, yet their mutual modulation mechanisms remain largely uncharted. Petroleum reservoirs are important lithosphere ecosystems that provide desirable resources for understanding microbial roles in element cycling. However, the strategy and mechanism of modulating indigenous microbial communities for the optimization of community structures and functions are underexplored, despite its significance in energy recovery and environmental remediation.

RESULTS

Here we proposed a novel selective stimulation of indigenous functional microbes by driving nitrogen and sulfur cycling in petroleum reservoirs using injections of an exogenous heterocycle-degrading strain of Pseudomonas. We defined such bacteria capable of removing and releasing organically bound sulfur and nitrogen from heterocycles as "bioredox triggers". High-throughput 16S rRNA amplicon sequencing, metagenomic, and gene transcription-level analyses of extensive production water and sandstone core samples spanning the whole oil production process clarified the microbiome dynamics following the intervention. These efforts demonstrated the feasibility of in situ N/S element release and electron acceptor generation during heterocycle degradation, shifting microbiome structures and functions and increasing phylogenetic diversity and genera engaged in sulfur and nitrogen cycling, such as Desulfovibrio, Shewanella, and Sulfurospirillum. The metabolic potentials of sulfur- and nitrogen-cycling processes, particularly dissimilatory sulfate reduction and dissimilatory nitrate reduction, were elevated in reservoir microbiomes. The relative expression of genes involved in sulfate reduction (dsrA, dsrB) and nitrate reduction (napA) was upregulated by 85, 28, and 22 folds, respectively. Field trials showed significant improvements in oil properties, with a decline in asphaltenes and aromatics, hetero-element contents, and viscosity, hence facilitating the effective exploitation of heavy oil.

CONCLUSIONS

The interactions between microbiomes and element cycling elucidated in this study will contribute to a better understanding of microbial metabolic involvement in, and response to, biogeochemical processes in the lithosphere. The presented findings demonstrated the immense potential of our microbial modulation strategy for green and enhanced heavy oil recovery. Video Abstract.

A new strategy for the enrichment of ammonia-oxidizing archaea in wastewater treatment systems: The positive role of quorum-sensing signaling molecules

Ammonia-oxidizing archaea (AOA) play an important role in natural nitrogen cycle, but are difficult to be enriched in wastewater treatment systems. In this experiment, under ambient temperature and high dissolved oxygen, different types of acyl-homoserine lactones (C₆-HSL, C₈-HSL, C₁₀-HSL, C₁₄-HSL and 3-oxo-C₁₄-HSL) were added to five wastewater nitrification systems to achieve AOA enrichment. Results showed that AOA couldn't be detected in the blank group without the addition of signaling molecules, while the AOA could be detected in all the reactors with the addition. The enrichment effect of AOA was not obvious with added 100 or 200 nmol/L signaling molecules, while the enrichment effect was both obvious with added C₈-HSL of 400 nmol/L and C₁₀-HSL of 800 nmol/L. And relative abundance of AOA increased from undetected in the control group to 1.10 % and 0.96 %, respectively. The exogenous signaling molecules may provide new view for AOA enrichment in wastewater treatment systems.

Genetic manipulation of candidate phyla radiation bacteria provides functional insights into microbial dark matter

The study of bacteria has yielded fundamental insights into cellular biology and physiology, biotechnological advances and many therapeutics. Yet due to a lack of suitable tools, the significant portion of bacterial diversity held within the candidate phyla radiation (CPR) remains inaccessible to such pursuits. Here we show that CPR bacteria belonging to the phylum Saccharibacteria exhibit natural competence. We exploit this property to develop methods for their genetic manipulation, including the insertion of heterologous sequences and the construction of targeted gene deletions. Imaging of fluorescent protein-labeled Saccharibacteria provides high spatiotemporal resolution of phenomena accompanying epibiotic growth and a transposon insertion sequencing genome-wide screen reveals the contribution of enigmatic Saccharibacterial genes to growth on their Actinobacteria hosts. Finally, we leverage metagenomic data to provide cutting-edge protein structure-based bioinformatic resources that support the strain Southlakia epibionticum and its corresponding host, Actinomyces israelii , as a model system for unlocking the molecular underpinnings of the epibiotic lifestyle.

Using Genomics to Decipher the Enigmatic Properties and Survival Adaptation of Candidate Phyla Radiation

Microbial ecology is a critical field for understanding the composition, diversity, and functions of microorganisms in various environmental and health-related processes. The discovery of Candidate Phyla Radiation (CPR) through culture-independent methods has introduced a new division of microbes characterized by a symbiotic/parasitic lifestyle, small cell size, and small genome. Despite being poorly understood, CPRs have garnered significant attention in recent years due to their widespread detection in a variety of environmental and clinical samples. These microorganisms have been found to exhibit a high degree of genetic diversity compared to other microbes. Several studies have shed light on their potential importance in global biogeochemical cycles and their impact on various human activities. In this review, we provide a systematic overview of the discovery of CPRs. We then focus on describing how the genomic characteristics of CPRs have helped them interact with and adapt to other microbes in different ecological niches. Future works should focus on discovering the metabolic capacities of CPRs and, if possible, isolating them to obtain a better understanding of these microorganisms.