Biosynthetic capacity, metabolic variety and unusual biology in the CPR and DPANN radiations

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.

Comprehensive analysis of insertion sequences within rRNA genes of CPR bacteria and biochemical characterization of a homing endonuclease encoded by these sequences

The Candidate Phyla Radiation (CPR) represents an extensive bacterial clade comprising primarily uncultured lineages and is distinguished from other bacteria by a significant prevalence of insertion sequences (ISs) within their rRNA genes. However, our understanding of the taxonomic distribution and characteristics of these ISs remains limited. In this study, we used a comprehensive approach to systematically determine the nature of the rRNA ISs in CPR bacteria. The analysis of hundreds of rRNA gene sequences across 65 CPR phyla revealed that ISs are present in 48% of 16S rRNA genes and 82% of 23S rRNA genes, indicating a broad distribution across the CPR clade, with exceptions in the 16S and 23S rRNA genes of Candidatus (Ca.) Saccharibacteria and the 16S rRNA genes of Ca. Peregrinibacteria. Over half the ISs display a group-I-intron-like structure, whereas specific 16S rRNA gene ISs display features reminiscent of group II introns. The ISs frequently encode proteins with homing endonuclease (HE) domains, centered around the LAGLIDADG motif. The LAGLIDADG HE (LHE) proteins encoded by the rRNA ISs of CPR bacteria predominantly have a single-domain structure, deviating from the usual single- or double-domain configuration observed in typical prokaryotic LHEs. Experimental analysis of one LHE protein, I-ShaI from Ca. Shapirobacteria, confirmed that its endonuclease activity targets the DNA sequence of its insertion site, and chemical cross-linking experiments demonstrated its capacity to form homodimers. These results provide robust evidence supporting the hypothesis that the explosive proliferation of rRNA ISs in CPR bacteria was facilitated by mechanisms involving LHEs.

IMPORTANCE

Insertion sequences (ISs) in rRNA genes are relatively limited and infrequent in most bacterial phyla. With a comprehensive bioinformatic analysis, we show that in CPR bacteria, these ISs occur in 48% of 16S rRNA genes and 82% of 23S rRNA genes. We also report the systematic and biochemical characterization of the LAGLIDADG homing endonucleases (LHEs) encoded by these ISs in the first such analysis of the CPR bacteria. This study significantly extends our understanding of the phylogenetic positions of rRNA ISs within CPR bacteria and the biochemical features of their LHEs.

Swimming towards each other: the role of chemotaxis in bacterial interactions

Chemotaxis allows microorganisms to direct movement in response to chemical stimuli. Bacteria use this behaviour to develop spatial associations with animals and plants, and even larger microbes. However, current theory suggests that constraints imposed by the limits of chemotactic sensory systems will prevent sensing of chemical gradients emanating from cells smaller than a few micrometres, precluding the utility of chemotaxis in interactions between individual bacteria. Yet, recent evidence has revealed surprising levels of bacterial chemotactic precision, as well as a role for chemotaxis in metabolite exchange between bacterial cells. If indeed widespread, chemotactic sensing between bacteria could represent an important, but largely overlooked, phenotype within interbacterial interactions, and play a significant role in shaping cooperative and competitive relationships.

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.

Vertical distribution of arsenic and bacterial communities in calcareous farmland amending by organic fertilizer and iron-oxidizing bacteria: Field experiment on concomitant remediation

The migration and transformation mechanisms of arsenic (As) in soil environments necessitate an understanding of its influencing processes. Here, we investigate the subsurface biogeochemical transformation of As and iron (Fe) through amended in the top 20 cm with iron oxidizing bacteria (FeOB) and organic fertilizer (OF). Our comprehensive 400-day field study, conducted in a calcareous soil profile sectioned into 20 cm increments, involved analysis by sequential extraction and assessment of microbial properties. The results reveal that the introduction of additional OF increased the release ratio of As/Fe from the non-specific adsorption fraction (136.47 %) at the subsoil depth (40-60 cm), underscoring the importance of sampling at various depths and time points to accurately elucidate the form, instability, and migration of As within the profile. Examination of bacterial interaction networks indicated a disrupted initial niche in the bottom layer, resulting in a novel cooperative symbiosis. While the addition of FeOB did not lead to the dominance of specific bacterial species, it did enhance the relative abundance of As-tolerant Acidobacteria and Gemmatimonadetes in both surface (39.2 % and 38.76 %) and deeper soils (44.29 % and 23.73 %) compared to the control. Consequently, the amendment of FeOB in conjunction with OF facilitated the formation of poorly amorphous Fe (hydr)oxides in the soil, achieved through abiotic and biotic sequestration processes. Throughout the long-term remediation process, the migration coefficient of bioavailable As within the soil profile decreased, indicating that these practices did not exacerbate As mobilization. This study carries significant implications for enhancing biogeochemical cycling in As-contaminated Sierozem soils and exploring potential bioremediation strategies. ENVIRONMENTAL IMPLICATION: The long-term exposure of sewage irrigation has potential adverse effects on the local ecosystem, causing serious environmental problems. Microorganisms play a vital role in the migration and transformation of arsenic in calcareous soil in arid areas, which highlights the necessity of understanding its dynamics. The vertical distribution, microbial community and fate of arsenic in calcareous farmland soil profile in northwest China were studied through field experiments. The results of this work have certain significance for the remediation of arsenic-contaminated soil in arid areas, and provide new insights for the migration, transformation and remediation of arsenic in this kind of soil.

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.

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.

Effects of diets for three growing stages by rumen inocula donors on in vitro rumen fermentation and microbiome

Hanwoo and Jeju Black cattle (Jeju Black) are native breeds of Korean cattle. Jeju Black cattle are recognized as natural monuments and are known to exhibit slower growth rates compared to Hanwoo. While several studies have analyzed the genetic characteristics of these cattle, there has been limited research on the differences in their microbiome. In this study, rumen fluid was obtained from three Hanwoo steers and three Jeju Black steers, and three different diets (total mixed rations [TMRs] for growing, early fattening, and late fattening periods) were used as substrates for in vitro fermentation. The in vitro incubation was conducted for 3 h and 24 h following a 2 × 3 factorial arrangement. After both incubation periods, fermentation characteristics were analyzed, and ruminal microbiome analysis was performed using 16S rRNA gene sequencing, employing both QIIME2 and PICRUSt2. The results revealed significant differences in the ruminal microbiota due to the inoculum effect. At the phylum level, Patescibacteria and Synergistota were found to be enriched in the Jeju Black inoculum-treated group. Additionally, using different inocula also affected the relative abundance of major taxa, including Ruminococcus, Pseudoramibacter, Ruminococcaceae CAG-352, and the [Eubacterium] ruminantium group. These microbial differences induced by the inoculum may have originated from varying levels of domestication between the two subspecies of donor animals, which mainly influenced the fermentation and microbiome features in the early incubation stages, although this was only partially offset afterward. Furthermore, predicted commission numbers of microbial enzymes, some of which are involved in the biosynthesis of secondary metabolites, fatty acids, and alpha amylase, differed based on the inoculum effect. However, these differences may account for only a small proportion of the overall metabolic pathway. Conversely, diets were found to affect protein biosynthesis and its related metabolism, which showed differential abundance in the growing diet and were potentially linked to the growth-promoting effects in beef cattle during the growing period. In conclusion, this study demonstrated that using different inocula significantly affected in vitro fermentation characteristics and microbiome features, mainly in the early stages of incubation, with some effects persisting up to 24 h of incubation.

Selective lipid recruitment by an archaeal DPANN symbiont from its host

The symbiont Ca. Nanohaloarchaeum antarcticus is obligately dependent on its host Halorubrum lacusprofundi for lipids and other metabolites due to its lack of certain biosynthetic genes. However, it remains unclear which specific lipids or metabolites are acquired from its host, and how the host responds to infection. Here, we explored the lipidome dynamics of the Ca. Nha. antarcticus - Hrr. lacusprofundi symbiotic relationship during co-cultivation. By using a comprehensive untargeted lipidomic methodology, our study reveals that Ca. Nha. antarcticus selectively recruits 110 lipid species from its host, i.e., nearly two-thirds of the total number of host lipids. Lipid profiles of co-cultures displayed shifts in abundances of bacterioruberins and menaquinones and changes in degree of bilayer-forming glycerolipid unsaturation. This likely results in increased membrane fluidity and improved resistance to membrane disruptions, consistent with compensation for higher metabolic load and mechanical stress on host membranes when in contact with Ca. Nha. antarcticus cells. Notably, our findings differ from previous observations of other DPANN symbiont-host systems, where no differences in lipidome composition were reported. Altogether, our work emphasizes the strength of employing untargeted lipidomics approaches to provide details into the dynamics underlying a DPANN symbiont-host system.

Exploring the prokaryote-eukaryote interplay in microbial mats from an Andean athalassohaline wetland

Microbial community assembly results from the interaction between biotic and abiotic factors. However, environmental selection is thought to predominantly shape communities in extreme ecosystems. Salar de Huasco, situated in the high-altitude Andean Altiplano, represents a poly-extreme ecosystem displaying spatial gradients of physicochemical conditions. To disentangle the influence of abiotic and biotic factors, we studied prokaryotic and eukaryotic communities from microbial mats and underlying sediments across contrasting areas of this athalassohaline ecosystem. The prokaryotic communities were primarily composed of bacteria, notably including a significant proportion of photosynthetic organisms like Cyanobacteria and anoxygenic photosynthetic members of Alpha- and Gammaproteobacteria and Chloroflexi. Additionally, Bacteroidetes, Verrucomicrobia, and Deltaproteobacteria were abundantly represented. Among eukaryotes, photosynthetic organisms (Ochrophyta and Archaeplastida) were predominant, alongside relatively abundant ciliates, cercozoans, and flagellated fungi. Salinity emerged as a key driver for the assembly of prokaryotic communities. Collectively, abiotic factors influenced both prokaryotic and eukaryotic communities, particularly those of algae. However, prokaryotic communities strongly correlated with photosynthetic eukaryotes, suggesting a pivotal role of biotic interactions in shaping these communities. Co-occurrence networks suggested potential interactions between different organisms, such as diatoms with specific photosynthetic and heterotrophic bacteria or with protist predators, indicating influences beyond environmental selection. While some associations may be explained by environmental preferences, the robust biotic correlations, alongside insights from other ecosystems and experimental studies, suggest that symbiotic and trophic interactions significantly shape microbial mat and sediment microbial communities in this athalassohaline ecosystem.IMPORTANCEHow biotic and abiotic factors influence microbial community assembly is still poorly defined. Here, we explore their influence on prokaryotic and eukaryotic community assembly within microbial mats and sediments of an Andean high-altitude polyextreme wetland system. We show that, in addition to abiotic elements, mutual interactions exist between prokaryotic and eukaryotic communities. Notably, photosynthetic eukaryotes exhibit a strong correlation with prokaryotic communities, specifically diatoms with certain bacteria and other protists. Our findings underscore the significance of biotic interactions in community assembly and emphasize the necessity of considering the complete microbial community.

Microbiological insights into anaerobic phenol degradation mechanisms and bulking phenomenon in a mesophilic upflow anaerobic sludge blanket reactor in long-term operation

In this study, a long-term operation of 2,747 days was conducted to evaluate the performance of the upflow anaerobic sludge blanket (UASB) reactor and investigated the degradation mechanisms of high-organic loading phenol wastewater. During the reactor operation, the maximum chemical oxygen demand (COD) removal rate of 6.1 ± 0.6 kg/m3/day under 1,680 mg/L phenol concentration was achieved in the mesophilic UASB reactor. After a significant change in the operating temperature from 24.0 ± 4.1 °C to 35.9 ± 0.6 °C, frequent observations of floating and washout of the bloated granular sludge (novel types of the bulking phenomenon) were made in the UASB reactor, suggesting that the change in operating temperature could be a trigger for the bulking phenomenon. Through the metagenomic analysis, phenol degradation mechanisms were predicted that phenol was converted to 4-hydroxybenzoate via two possible routes by Syntrophorhabdaceae and Pelotomaculaceae bacteria. Furthermore, the degradation of 4-hydroxybenzoate to benzoyl-CoA was carried out by members of Syntrophorhabdaceae and Smithellaceae. In the bulking sludge, a predominant presence of Nanobdellota, belonging to DPANN archaea, was detected. The metagenome-assembled genome of the Nanobdellota lacks many biosynthetic pathways and has several genes for the symbiotic lifestyle such as trimeric autotransporter adhesin-related protein. Furthermore, the Nanobdellota have significant correlations with several methanogenic archaea that are predominantly present in the UASB reactor. Considering the results of this study, the predominant Nanobdellota may negatively affect the growth of the methanogens through the parasitic lifestyle and change the balance of microbial interactions in the granular sludge ecosystem.

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.

How do bacterial endosymbionts work with so few genes?

The move from a free-living environment to a long-term residence inside a host eukaryotic cell has profound effects on bacterial function. While endosymbioses are found in many eukaryotes, from protists to plants to animals, the bacteria that form these host-beneficial relationships are even more diverse. Endosymbiont genomes can become radically smaller than their free-living relatives, and their few remaining genes show extreme compositional biases. The details of how these reduced and divergent gene sets work, and how they interact with their host cell, remain mysterious. This Unsolved Mystery reviews how genome reduction alters endosymbiont biology and highlights a "tipping point" where the loss of the ability to build a cell envelope coincides with a marked erosion of translation-related genes.

Expanded phylogeny of extremely halophilic archaea shows multiple independent adaptations to hypersaline environments

Extremely halophilic archaea (Haloarchaea, Nanohaloarchaeota, Methanonatronarchaeia and Halarchaeoplasmatales) thrive in saturating salt concentrations where they must maintain osmotic equilibrium with their environment. The evolutionary history of adaptations enabling salt tolerance remains poorly understood, in particular because the phylogeny of several lineages is conflicting. Here we present a resolved phylogeny of extremely halophilic archaea obtained using improved taxon sampling and state-of-the-art phylogenetic approaches designed to cope with the strong compositional biases of their proteomes. We describe two uncultured lineages, Afararchaeaceae and Asbonarchaeaceae, which break the long branches at the base of Haloarchaea and Nanohaloarchaeota, respectively. We obtained 13 metagenome-assembled genomes (MAGs) of these archaea from metagenomes of hypersaline aquatic systems of the Danakil Depression (Ethiopia). Our phylogenomic analyses including these taxa show that at least four independent adaptations to extreme halophily occurred during archaeal evolution. Gene-tree/species-tree reconciliation suggests that gene duplication and horizontal gene transfer played an important role in this process, for example, by spreading key genes (such as those encoding potassium transporters) across extremely halophilic lineages.

Global abundance patterns, diversity, and ecology of Patescibacteria in wastewater treatment plants

BACKGROUND

Microorganisms are responsible for nutrient removal and resource recovery in wastewater treatment plants (WWTPs), and their diversity is often studied by 16S rRNA gene amplicon sequencing. However, this approach underestimates the abundance and diversity of Patescibacteria due to the low coverage of commonly used PCR primers for this highly divergent bacterial phylum. Therefore, our current understanding of the global diversity, distribution, and ecological role of Patescibacteria in WWTPs is very incomplete. This is particularly relevant as Patescibacteria are considered to be associated with microbial host cells and can therefore influence the abundance and temporal variability of other microbial groups that are important for WWTP functioning.

RESULTS

Here, we evaluated the in silico coverage of widely used 16S rRNA gene-targeted primer pairs and redesigned a primer pair targeting the V4 region of bacterial and archaeal 16S rRNA genes to expand its coverage for Patescibacteria. We then experimentally evaluated and compared the performance of the original and modified V4-targeted primers on 565 WWTP samples from the MiDAS global sample collection. Using the modified primer pair, the percentage of ASVs classified as Patescibacteria increased from 5.9 to 23.8%, and the number of detected patescibacterial genera increased from 560 to 1576, while the detected diversity of the remaining microbial community remained similar. Due to this significantly improved coverage of Patescibacteria, we identified 23 core genera of Patescibacteria in WWTPs and described the global distribution pattern of these unusual microbes in these systems. Finally, correlation network analysis revealed potential host organisms that might be associated with Patescibacteria in WWTPs. Interestingly, strong indications were found for an association between Patescibacteria of the Saccharimonadia and globally abundant polyphosphate-accumulating organisms of the genus Ca. Phosphoribacter.

CONCLUSIONS

Our study (i) provides an improved 16S rRNA gene V4 region-targeted amplicon primer pair inclusive of Patescibacteria with little impact on the detection of other taxa, (ii) reveals the diversity and distribution patterns of Patescibacteria in WWTPs on a global scale, and (iii) provides new insights into the ecological role and potential hosts of Patescibacteria in WWTPs. Video Abstract.

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.

Microscopic and metatranscriptomic analyses revealed unique cross-domain parasitism between phylum Candidatus Patescibacteria/candidate phyla radiation and methanogenic archaea in anaerobic ecosystems

To verify whether members of the phylum Candidatus Patescibacteria parasitize archaea, we applied cultivation, microscopy, metatranscriptomic, and protein structure prediction analyses on the Patescibacteria-enriched cultures derived from a methanogenic bioreactor. Amendment of cultures with exogenous methanogenic archaea, acetate, amino acids, and nucleoside monophosphates increased the relative abundance of Ca. Patescibacteria. The predominant Ca. Patescibacteria were families Ca. Yanofskyibacteriaceae and Ca. Minisyncoccaceae, and the former showed positive linear relationships (r2 ≥ 0.70) Methanothrix in their relative abundances, suggesting related growth patterns. Methanothrix and Methanospirillum cells with attached Ca. Yanofskyibacteriaceae and Ca. Minisyncoccaceae, respectively, had significantly lower cellular activity than those of the methanogens without Ca. Patescibacteria, as extrapolated from fluorescence in situ hybridization-based fluorescence. We also observed that parasitized methanogens often had cell surface deformations. Some Methanothrix-like filamentous cells were dented where the submicron cells were attached. Ca. Yanofskyibacteriaceae and Ca. Minisyncoccaceae highly expressed extracellular enzymes, and based on structural predictions, some contained peptidoglycan-binding domains with potential involvement in host cell attachment. Collectively, we propose that the interactions of Ca. Yanofskyibacteriaceae and Ca. Minisyncoccaceae with methanogenic archaea are parasitisms.IMPORTANCECulture-independent DNA sequencing approaches have explored diverse yet-to-be-cultured microorganisms and have significantly expanded the tree of life in recent years. One major lineage of the domain Bacteria, Ca. Patescibacteria (also known as candidate phyla radiation), is widely distributed in natural and engineered ecosystems and has been thought to be dependent on host bacteria due to the lack of several biosynthetic pathways and small cell/genome size. Although bacteria-parasitizing or bacteria-preying Ca. Patescibacteria have been described, our recent studies revealed that some lineages can specifically interact with archaea. In this study, we provide strong evidence that the relationship is parasitic, shedding light on overlooked roles of Ca. Patescibacteria in anaerobic habitats.

CPR bacteria and DPANN archaea play pivotal roles in response of microbial community to antibiotic stress in groundwater

The accumulation of antibiotics in the natural environment can disrupt microbial population dynamics. However, our understanding of how microbial communities adapt to the antibiotic stress in groundwater ecosystems remains limited. By recovering 2675 metagenome-assembled genomes (MAGs) from 66 groundwater samples, we explored the effect of antibiotics on bacterial, archaeal, and fungal communities, and revealed the pivotal microbes and their mechanisms in coping with antibiotic stress. The results indicated that antibiotics had the most significant influence on bacterial and archaeal communities, while the impact on the fungal community was minimal. Analysis of co-occurrence networks between antibiotics and microbes revealed the critical roles of Candidate Phyla Radiation (CPR) bacteria and DPANN archaea, two representative microbial groups in groundwater ecosystem, in coping with antibiotic resistance and enhancing network connectivity and complexity. Further genomic analysis demonstrated that CPR bacteria carried approximately 6 % of the identified antibiotic resistance genes (ARGs), indicating their potential to withstand antibiotics on their own. Meanwhile, the genomes of CPR bacteria and DPANN archaea were found to encode diverse biosynthetic gene clusters (BGCs) responsible for producing antimicrobial metabolites, which could not only assist CPR and DPANN organisms but also benefit the surrounding microbes in combating antibiotic stress. These findings underscore the significant impact of antibiotics on prokaryotic microbial communities in groundwater, and highlight the importance of CPR bacteria and DPANN archaea in enhancing the overall resilience and functionality of the microbial community in the face of antibiotic stress.

Cell surface architecture of the cultivated DPANN archaeon Nanobdella aerobiophila

The DPANN archaeal clade includes obligately ectosymbiotic species. Their cell surfaces potentially play an important role in the symbiotic interaction between the ectosymbionts and their hosts. However, little is known about the mechanism of ectosymbiosis. Here, we show cell surface structures of the cultivated DPANN archaeon Nanobdella aerobiophila strain MJ1T and its host Metallosphaera sedula strain MJ1HA, using a variety of electron microscopy techniques, i.e., negative-staining transmission electron microscopy, quick-freeze deep-etch TEM, and 3D electron tomography. The thickness, unit size, and lattice symmetry of the S-layer of strain MJ1T were different from those of the host archaeon strain MJ1HA. Genomic and transcriptomic analyses highlighted the most highly expressed MJ1T gene for a putative S-layer protein with multiple glycosylation sites and immunoglobulin-like folds, which has no sequence homology to known S-layer proteins. In addition, genes for putative pectin lyase- or lectin-like extracellular proteins, which are potentially involved in symbiotic interaction, were found in the MJ1T genome based on in silico 3D protein structure prediction. Live cell imaging at the optimum growth temperature of 65°C indicated that cell complexes of strains MJ1T and MJ1HA were motile, but sole MJ1T cells were not. Taken together, we propose a model of the symbiotic interaction and cell cycle of Nanobdella aerobiophila.IMPORTANCEDPANN archaea are widely distributed in a variety of natural and artificial environments and may play a considerable role in the microbial ecosystem. All of the cultivated DPANN archaea so far need host organisms for their growth, i.e., obligately ectosymbiotic. However, the mechanism of the ectosymbiosis by DPANN archaea is largely unknown. To this end, we performed a comprehensive analysis of the cultivated DPANN archaeon, Nanobdella aerobiophila, using electron microscopy, live cell imaging, transcriptomics, and genomics, including 3D protein structure prediction. Based on the results, we propose a reasonable model of the symbiotic interaction and cell cycle of Nanobdella aerobiophila, which will enhance our understanding of the enigmatic physiology and ecological significance of DPANN archaea.

Comamonas-dominant microbial community in carbon poor aquitard sediments revealed by metagenomic-based growth rate investigation

The microbiological ecology of a low-nutrient shallow aquifer with high arsenic content in the Yinchuan Plain was investigated in this study. Amplicon sequencing data from five samples (depths: 1.5 m, 3.5 m, 11.2 m, 19.3 m, and 25.5 m) revealed diverse and adaptable microbial community. Among the microbial community, Comamonas was the most prominent, accounting for 10.52 % of the total. This genus displayed high growth rates, with a maximum growth rate of 12.06 d-1 and a corresponding doubling time of 1.38 days, as determined through an analysis of codon usage bias. Functional annotation of Metagenome-Assembled Genomes (MAGs) for samples at 1.5 m and 11.2 m depths revealed Comamonas' metabolic versatility, including various carbon pathways, assimilative sulfate reduction (ASR), and dissimilatory reduction to ammonium (DNRA). The TPM (Transcripts Per Kilobase of exon model per Million mapped reads) of MAGs at 11.2 m sample was 15.7 and 12.3. The presence of arsenic resistance genes in Comamonas aligns with sediment arsenic levels (65.8 mg/kg for 1.5 m depth, 32.8 mg/kg for 11.2 m depth). This study highlights the role of Comamonas as a 'generalist' bacteria in challenging oligotrophic sediments, emphasizing the significance of such organisms in community stability and ecological functions. ENVIRONMENTAL IMPLICATION: Low-biomass limits the microbial activity and biogeochemical study in oligotrophic environments, which is the typical condition for underground aquatic ecosystems. Facilitated by growth rate estimation, our research focuses on active functional microorganisms and their biogeochemical metabolic in oligotrophic aquifer sediments, revealing their impact on the environment and response to arsenic threats. Findings illuminate the metabolic advantage of a 'generalist life-style' in carbon-scarce environments and contribute to a broader understanding of bacterial ecosystems and environmental impacts in oligotrophic aquifer sediments worldwide.

Bacterial ecology and evolution converge on seasonal and decadal scales

Ecology and evolution are distinct theories, but the short lifespans and large population sizes of microbes allow evolution to unfold along contemporary ecological time scales. To document this in a natural system, we collected a two-decade, 471-metagenome time series from a single site in a freshwater lake, which we refer to as the TYMEFLIES dataset. This massive sampling and sequencing effort resulted in the reconstruction of 30,389 metagenomic-assembled genomes (MAGs) over 50% complete, which dereplicated into 2,855 distinct genomes (>96% nucleotide sequence identity). We found both ecological and evolutionary processes occurred at seasonal time scales. There were recurring annual patterns at the species level in abundances, nucleotide diversities (π), and single nucleotide variant (SNV) profiles for the majority of all taxa. During annual blooms, we observed both higher and lower nucleotide diversity, indicating that both ecological differentiation and competition drove evolutionary dynamics. Overlayed upon seasonal patterns, we observed long-term change in 20% of the species' SNV profiles including gradual changes, step changes, and disturbances followed by resilience. Most abrupt changes occurred in a single species, suggesting evolutionary drivers are highly specific. Nevertheless, seven members of the abundant Nanopelagicaceae family experienced abrupt change in 2012, an unusually hot and dry year. This shift coincided with increased numbers of genes under selection involved in amino acid and nucleic acid metabolism, suggesting fundamental organic nitrogen compounds drive strain differentiation in the most globally abundant freshwater family. Overall, we observed seasonal and decadal trends in both interspecific ecological and intraspecific evolutionary processes. The convergence of microbial ecology and evolution on the same time scales demonstrates that understanding microbiomes requires a new unified approach that views ecology and evolution as a single continuum.

Functional and evolutionary significance of unknown genes from uncultivated taxa

Many of the Earth's microbes remain uncultured and understudied, limiting our understanding of the functional and evolutionary aspects of their genetic material, which remain largely overlooked in most metagenomic studies1. Here we analysed 149,842 environmental genomes from multiple habitats2-6 and compiled a curated catalogue of 404,085 functionally and evolutionarily significant novel (FESNov) gene families exclusive to uncultivated prokaryotic taxa. All FESNov families span multiple species, exhibit strong signals of purifying selection and qualify as new orthologous groups, thus nearly tripling the number of bacterial and archaeal gene families described to date. The FESNov catalogue is enriched in clade-specific traits, including 1,034 novel families that can distinguish entire uncultivated phyla, classes and orders, probably representing synapomorphies that facilitated their evolutionary divergence. Using genomic context analysis and structural alignments we predicted functional associations for 32.4% of FESNov families, including 4,349 high-confidence associations with important biological processes. These predictions provide a valuable hypothesis-driven framework that we used for experimental validatation of a new gene family involved in cell motility and a novel set of antimicrobial peptides. We also demonstrate that the relative abundance profiles of novel families can discriminate between environments and clinical conditions, leading to the discovery of potentially new biomarkers associated with colorectal cancer. We expect this work to enhance future metagenomics studies and expand our knowledge of the genetic repertory of uncultivated organisms.

Viral AMGs-driven pentose phosphate pathway in natural wetland

Viruses exist anywhere on earth where there is life, and among them, virus-encoded auxiliary metabolic genes (AMGs) can maintain ecosystem balance and play a major role in the global ecosystem. Although the function of AMGs has been widely reported, the genetic diversity of AMGs in natural ecosystems is still poorly understood. Exploring the genetic diversity of viral community-wide AMGs is essential to gain insight into the complex interactions between viruses and hosts. In this article, we studied the phylogenetic tree, principal co-ordinates analysis (PCoA), α diversity, and metabolic pathways of viral auxiliary metabolism genes involved in the pentose phosphate pathway (PPP) through metagenomics, and the changes of metabolites and genes of host bacteria were further studied by using Pseudomonas mandelii SW-3 and its lytic phage based on metabolic flow and AMGs expression. We found that the viral AMGs in the Napahai plateau wetland were created by a combination of various external forces, which contributed to the rich genetic diversity, uniqueness, and differences of the virus, which promoted the reproduction of offspring and better adaptation to the environment. Overall, this study systematically describes the genetic diversity of AMGs associated with the PPP in plateau wetland ecosystems and further expands the understanding of phage-host unique interactions.

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.

Microbial sensor variation across biogeochemical conditions in the terrestrial deep subsurface

Microbes can be found in abundance many kilometers underground. While microbial metabolic capabilities have been examined across different geochemical settings, it remains unclear how changes in subsurface niches affect microbial needs to sense and respond to their environment. To address this question, we examined how microbial extracellular sensor systems vary with environmental conditions across metagenomes at different Deep Mine Microbial Observatory (DeMMO) subsurface sites. Because two-component systems (TCSs) directly sense extracellular conditions and convert this information into intracellular biochemical responses, we expected that this sensor family would vary across isolated oligotrophic subterranean environments that differ in abiotic and biotic conditions. TCSs were found at all six subsurface sites, the service water control, and the surface site, with an average of 0.88 sensor histidine kinases (HKs) per 100 genes across all sites. Abundance was greater in subsurface fracture fluids compared with surface-derived fluids, and candidate phyla radiation (CPR) bacteria presented the lowest HK frequencies. Measures of microbial diversity, such as the Shannon diversity index, revealed that HK abundance is inversely correlated with microbial diversity (r2 = 0.81). Among the geochemical parameters measured, HK frequency correlated most strongly with variance in dissolved organic carbon (r2 = 0.82). Taken together, these results implicate the abiotic and biotic properties of an ecological niche as drivers of sensor needs, and they suggest that microbes in environments with large fluctuations in organic nutrients (e.g., lacustrine, terrestrial, and coastal ecosystems) may require greater TCS diversity than ecosystems with low nutrients (e.g., open ocean).IMPORTANCEThe ability to detect extracellular environmental conditions is a fundamental property of all life forms. Because microbial two-component sensor systems convert information about extracellular conditions into biochemical information that controls their behaviors, we evaluated how two-component sensor systems evolved within the deep Earth across multiple sites where abiotic and biotic properties vary. We show that these sensor systems remain abundant in microbial consortia at all subterranean sampling sites and observe correlations between sensor system abundances and abiotic (dissolved organic carbon variation) and biotic (consortia diversity) properties. These results suggest that multiple environmental properties may drive sensor protein evolution and highlight the need for further studies of metagenomic and geochemical data in parallel to understand the drivers of microbial sensor evolution.

Distinct life cycle stages of an ectosymbiotic DPANN archaeon

DPANN archaea are a diverse group of microorganisms that are thought to rely on an ectosymbiotic lifestyle; however, the cell biology of these cell-cell interactions remains largely unknown. We applied live-cell imaging and cryo-electron tomography to the DPANN archaeon Nanobdella aerobiophila and its host, revealing two distinct life cycle stages. Free cells possess archaella and are motile. Ectobiotic cells are intimately linked with the host through an elaborate attachment organelle. Our data suggest that free cells may actively seek a new host, while the ectobiotic state is adapted to mediate intricate interaction with the host.

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.

Metabolic profiling of petroleum-degrading microbial communities incubated under high-pressure conditions

While pressure is a significant characteristic of petroleum reservoirs, it is often overlooked in laboratory studies. To clarify the composition and metabolic properties of microbial communities under high-pressure conditions, we established methanogenic and sulfate-reducing enrichment cultures under high-pressure conditions using production water from the Jilin Oilfield in China. We utilized a metagenomics approach to analyze the microbial community after a 90-day incubation period. Under methanogenic conditions, Firmicutes, Deferribacteres, Ignavibacteriae, Thermotogae, and Nitrospirae, in association with the hydrogenotrophic methanogen Archaeoglobaceae and acetoclastic Methanosaeta, were highly represented. Genomes for Ca. Odinarchaeota and the hydrogen-dependent methylotrophic Ca. Methanosuratus were also recovered from the methanogenic culture. The sulfate-reducing community was dominated by Firmicutes, Thermotogae, Nitrospirae, Archaeoglobus, and several candidate taxa including Ca. Bipolaricaulota, Ca. Aminicenantes, and Candidate division WOR-3. These candidate taxa were key pantothenate producers for other community members. The study expands present knowledge of the metabolic roles of petroleum-degrading microbial communities under high-pressure conditions. Our results also indicate that microbial community interactions were shaped by syntrophic metabolism and the exchange of amino acids and cofactors among members. Furthermore, incubation under in situ pressure conditions has the potential to reveal the roles of microbial dark matter.

Host translation machinery is not a barrier to phages that interact with both CPR and non-CPR bacteria

IMPORTANCE

Here, we profiled putative phages of Saccharibacteria, which are of particular importance as Saccharibacteria influence some human oral diseases. We additionally profiled putative phages of Gracilibacteria and Absconditabacteria, two Candidate Phyla Radiation (CPR) lineages of interest given their use of an alternative genetic code. Among the phages identified in this study, some are targeted by spacers from both CPR and non-CPR bacteria and others by both bacteria that use the standard genetic code as well as bacteria that use an alternative genetic code. These findings represent new insights into possible phage replication strategies and have relevance for phage therapies that seek to manipulate microbiomes containing CPR bacteria.

Scaling of Protein Function across the Tree of Life

Scaling laws are a powerful way to compare genomes because they put all organisms onto a single curve and reveal nontrivial generalities as genomes change in size. The abundance of functional categories across genomes has previously been found to show power law scaling with respect to the total number of functional categories, suggesting that universal constraints shape genomic category abundance. Here, we look across the tree of life to understand how genome evolution may be related to functional scaling. We revisit previous observations of functional genome scaling with an expanded taxonomy by analyzing 3,726 bacterial, 220 archaeal, and 79 unicellular eukaryotic genomes. We find that for some functional classes, scaling is best described by multiple exponents, revealing previously unobserved shifts in scaling as genome-encoded protein annotations increase or decrease. Furthermore, we find that scaling varies between phyletic groups at both the domain and phyla levels and is less universal than previously thought. This variability in functional scaling is not related to taxonomic phylogeny resolved at the phyla level, suggesting that differences in cell plan or physiology outweigh broad patterns of taxonomic evolution. Since genomes are maintained and replicated by the functional proteins encoded by them, these results point to functional degeneracy between taxonomic groups and unique evolutionary trajectories toward these. We also find that individual phyla frequently span scaling exponents of functional classes, revealing that individual clades can move across scaling exponents. Together, our results reveal unique shifts in functions across the tree of life and highlight that as genomes grow or shrink, proteins of various functions may be added or lost.

First shotgun metagenomics study of Juan de Fuca deep-sea sediments reveals distinct microbial communities above, within, between, and below sulfate methane transition zones

The marine deep subsurface is home to a vast microbial ecosystem, affecting biogeochemical cycles on a global scale. One of the better-studied deep biospheres is the Juan de Fuca (JdF) Ridge, where hydrothermal fluid introduces oxidants into the sediment from below, resulting in two sulfate methane transition zones (SMTZs). In this study, we present the first shotgun metagenomics study of unamplified DNA from sediment samples from different depths in this stratified environment. Bioinformatic analyses showed a shift from a heterotrophic, Chloroflexota-dominated community above the upper SMTZ to a chemolithoautotrophic Proteobacteria-dominated community below the secondary SMTZ. The reintroduction of sulfate likely enables respiration and boosts active cells that oxidize acetate, iron, and complex carbohydrates to degrade dead biomass in this low-abundance, low-diversity environment. In addition, analyses showed many proteins of unknown function as well as novel metagenome-assembled genomes (MAGs). The study provides new insights into microbial communities in this habitat, enabled by an improved DNA extraction protocol that allows a less biased view of taxonomic composition and metabolic activities, as well as uncovering novel taxa. Our approach presents the first successful attempt at unamplified shotgun sequencing samples from beyond 50 meters below the seafloor and opens new ways for capturing the true diversity and functional potential of deep-sea sediments.

ATP synthase evolution on a cross-braced dated tree of life

The timing of early cellular evolution, from the divergence of Archaea and Bacteria to the origin of eukaryotes, is poorly constrained. The ATP synthase complex is thought to have originated prior to the Last Universal Common Ancestor (LUCA) and analyses of ATP synthase genes, together with ribosomes, have played a key role in inferring and rooting the tree of life. We reconstruct the evolutionary history of ATP synthases using an expanded taxon sampling set and develop a phylogenetic cross-bracing approach, constraining equivalent speciation nodes to be contemporaneous, based on the phylogenetic imprint of endosymbioses and ancient gene duplications. This approach results in a highly resolved, dated species tree and establishes an absolute timeline for ATP synthase evolution. Our analyses show that the divergence of ATP synthase into F- and A/V-type lineages was a very early event in cellular evolution dating back to more than 4 Ga, potentially predating the diversification of Archaea and Bacteria. Our cross-braced, dated tree of life also provides insight into more recent evolutionary transitions including eukaryogenesis, showing that the eukaryotic nuclear and mitochondrial lineages diverged from their closest archaeal (2.67-2.19 Ga) and bacterial (2.58-2.12 Ga) relatives at approximately the same time, with a slightly longer nuclear stem-lineage.

Current knowledge and clinical perspectives for a unique new phylum: Nanaorchaeota

Nanoarchaea measuring less than 500 nm and encasing an average 600-kb compact genome have been studied for twenty years, after an estimated 4193-million-year evolution. Comprising only four co-cultured representatives, these symbiotic organisms initially detected in deep-sea hydrothermal vents and geothermal springs, have been further distributed in various environmental ecosystems worldwide. Recent isolation by co-culture of Nanopusillus massiliensis from the unique ecosystem of the human oral cavity, prompted us to review the evolutionary diversity of nanaorchaea resulting in a rapidly evolving taxonomiy. Regardless of their ecological niche, all nanoarchaea share limited metabolic capacities correlating with an obligate ectosymbiotic or parasitic lifestyle; focusing on the dynamics of nanoarchaea-bacteria nanoarchaea-archaea interactions at the morphological and metabolic levels; highlighting proteins involved in nanoarchaea attachment to the hosts, as well metabolic exchanges between both organisms; and highlighting clinical nanoarchaeology, an emerging field of research in the frame of the recent discovery of Candidate Phyla radiation (CPR) in human microbiota. Future studies in clinical nanobiology will expand knowledge of the nanaorchaea repertoire associated with human microbiota and diseases, to improve our understanding of the diversity of these nanoorganims and their intreactions with microbiota and host tissues.

Resource partitioning and amino acid assimilation in a terrestrial geothermal spring

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

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

Archaea belonging to the DPANN superphylum have been found within 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 15-26% 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 belonging to DPANN phyla Nanoarchaeota, Pacearchaeota, Woesarchaeota, Undinarchaeota, Iainarchaeota, and SpSt-1190 from pelagic ODZs in the Eastern Tropical North Pacific and Arabian Sea. We find these archaea to be permanent, stable residents of all 3 major ODZs only within anoxic depths, comprising up to 1% of the total microbial community and up to 25-50% of archaea. ODZ DPANN appear 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.

Genetic manipulation of Patescibacteria provides mechanistic insights into microbial dark matter and the epibiotic lifestyle

Patescibacteria, also known as the candidate phyla radiation (CPR), are a diverse group of bacteria that constitute a disproportionately large fraction of microbial dark matter. Its few cultivated members, belonging mostly to Saccharibacteria, grow as epibionts on host Actinobacteria. Due to a lack of suitable tools, the genetic basis of this lifestyle and other unique features of Patescibacteira remain unexplored. Here, we show that Saccharibacteria exhibit natural competence, and we exploit this property for their genetic manipulation. Imaging of fluorescent protein-labeled Saccharibacteria provides high spatiotemporal resolution of phenomena accompanying epibiotic growth, and a transposon-insertion sequencing (Tn-seq) genome-wide screen reveals the contribution of enigmatic Saccharibacterial genes to growth on their 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.

Metagenome-derived virus-microbe ratios across ecosystems

It is generally assumed that viruses outnumber cells on Earth by at least tenfold. Virus-to-microbe ratios (VMR) are largely based on counts of fluorescently labelled virus-like particles. However, these exclude intracellular viruses and potentially include false positives (DNA-containing vesicles, gene-transfer agents, unspecifically stained inert particles). Here, we develop a metagenome-based VMR estimate (mVRM) that accounts for DNA viruses across all stages of their replication cycles (virion, intracellular lytic and lysogenic) by using normalised RPKM (reads per kilobase of gene sequence per million of mapped metagenome reads) counts of the major capsid protein (MCP) genes and cellular universal single-copy genes (USCGs) as proxies for virus and cell counts, respectively. After benchmarking this strategy using mock metagenomes with increasing VMR, we inferred mVMR across different biomes. To properly estimate mVMR in aquatic ecosystems, we generated metagenomes from co-occurring cellular and viral fractions (>50 kDa-200 µm size-range) in freshwater, seawater and solar saltern ponds (10 metagenomes, 2 control metaviromes). Viruses outnumbered cells in freshwater by ~13 fold and in plankton from marine and saline waters by ~2-4 fold. However, across an additional set of 121 diverse non-aquatic metagenomes including microbial mats, microbialites, soils, freshwater and marine sediments and metazoan-associated microbiomes, viruses, on average, outnumbered cells by barely two-fold. Although viruses likely are the most diverse biological entities on Earth, their global numbers might be closer to those of cells than previously estimated.

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.

A predicted CRISPR-mediated symbiosis between uncultivated archaea

CRISPR-Cas systems defend prokaryotic cells from invasive DNA of viruses, plasmids and other mobile genetic elements. Here, we show using metagenomics, metatranscriptomics and single-cell genomics that CRISPR systems of widespread, uncultivated archaea can also target chromosomal DNA of archaeal episymbionts of the DPANN superphylum. Using meta-omics datasets from Crystal Geyser and Horonobe Underground Research Laboratory, we find that CRISPR spacers of the hosts Candidatus Altiarchaeum crystalense and Ca. A. horonobense, respectively, match putative essential genes in their episymbionts' genomes of the genus Ca. Huberiarchaeum and that some of these spacers are expressed in situ. Metabolic interaction modelling also reveals complementation between host-episymbiont systems, on the basis of which we propose that episymbionts are either parasitic or mutualistic depending on the genotype of the host. By expanding our analysis to 7,012 archaeal genomes, we suggest that CRISPR-Cas targeting of genomes associated with symbiotic archaea evolved independently in various archaeal lineages.

Metabolic implications for predatory and parasitic bacterial lineages in activated sludge wastewater treatment systems

Deciphering unclear microbial interactions is key to improving biological wastewater treatment processes. Microbial predation and parasitism in wastewater treatment ecosystems are unexplored survival strategies that have long been known and have recently attracted attention because these interspecies interactions may contribute to the reduction of excess sludge. Here, microbial community profiling of 600 activated sludge samples taken from six industrial and one municipal wastewater treatment processes (WWTPs) was conducted. To identify the shared lineages in the WWTPs, the shared microbial constituents were defined as the family level taxa that had ≥ 0.1% average relative abundance and detected in all processes. The microbial community analysis assigned 106 families as the shared microbial constituents in the WWTPs. Correlation analysis showed that 98 of the 106 shared families were significantly correlated with total carbon (TC) and/or total nitrogen (TN) concentrations, suggesting that they may contribute to wastewater remediation. Most possible predatory or parasitic bacteria belonging to the phyla Bdellovibrionota, Myxococcota, and Candidatus Patescibacteria were found to be the shared families and negatively correlated with TC/TN; thus, they were frequently present in the WWTPs and could be involved in the removal of carbon/nitrogen derived from cell components. Shotgun metagenome-resolved metabolic reconstructions indicated that gene homologs associated with predation or parasitism are conserved in the Bdellovibrionota, Myxococcota, and Ca. Patescibacteria genomes (e.g., host interaction (hit) locus, Tad-like secretion complexes, and type IV pilus assembly proteins). This study provides insights into the complex microbial interactions potentially linked to the reduction of excess sludge biomass in these processes.

CheckM2: a rapid, scalable and accurate tool for assessing microbial genome quality using machine learning

Advances in sequencing technologies and bioinformatics tools have dramatically increased the recovery rate of microbial genomes from metagenomic data. Assessing the quality of metagenome-assembled genomes (MAGs) is a critical step before downstream analysis. Here, we present CheckM2, an improved method of predicting genome quality of MAGs using machine learning. Using synthetic and experimental data, we demonstrate that CheckM2 outperforms existing tools in both accuracy and computational speed. In addition, CheckM2's database can be rapidly updated with new high-quality reference genomes, including taxa represented only by a single genome. We also show that CheckM2 accurately predicts genome quality for MAGs from novel lineages, even for those with reduced genome size (for example, Patescibacteria and the DPANN superphylum). CheckM2 provides accurate genome quality predictions across bacterial and archaeal lineages, giving increased confidence when inferring biological conclusions from MAGs.

Metallo-Beta-Lactamase-like Encoding Genes in Candidate Phyla Radiation: Widespread and Highly Divergent Proteins with Potential Multifunctionality

The Candidate Phyla Radiation (CPR) was found to harbor a vast repertoire of genes encoding for enzymes with potential antibiotic resistance activity. Among these, as many as 3349 genes were predicted in silico to contain a metallo-beta-lactamase-like (MBL-like) fold. These proteins were subject to an in silico functional characterization by comparing their protein profiles (presence/absence of conserved protein domains) to other MBLs, including 24 already expressed in vitro, along with those of the beta-lactamase database (BLDB) (n = 761). The sequence similarity network (SSN) was then used to predict the functional clusters of CPR MBL-like sequences. Our findings showed that CPR MBL-like sequences were longer and more diverse than bacterial MBL sequences, with a high content of functional domains. Most CPR MBL-like sequences did not show any SSN connectivity with expressed MBLs, indicating the presence of many potential, yet unidentified, functions in CPR. In conclusion, CPR was shown to have many protein functions and a large sequence variability of MBL-like folds, exceeding all known MBLs. Further experimental and evolutionary studies of this superfamily of hydrolyzing enzymes are necessary to illustrate their functional annotation, origin, and expansion for adaptation or specialization within a given niche or compared to a specific substrate.

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.

Preliminary landscape of Candidatus Saccharibacteria in the human microbiome

Introduction

Candidate Phyla Radiation (CPR) and more specifically Candidatus Saccharibacteria (TM7) have now been established as ubiquitous members of the human oral microbiota. Additionally, CPR have been reported in the gastrointestinal and urogenital tracts. However, the exploration of new human niches has been limited to date.

Methods

In this study, we performed a prospective and retrospective screening of TM7 in human samples using standard PCR, real-time PCR, scanning electron microscopy (SEM) and shotgun metagenomics.

Results

Using Real-time PCR and standard PCR, oral samples presented the highest TM7 prevalence followed by fecal samples, breast milk samples, vaginal samples and urine samples. Surprisingly, TM7 were also detected in infectious samples, namely cardiac valves and blood cultures at a low prevalence (under 3%). Moreover, we observed CPR-like structures using SEM in all sample types except cardiac valves. The reconstruction of TM7 genomes in oral and fecal samples from shotgun metagenomics reads further confirmed their high prevalence in some samples.

Conclusion

This study confirmed, through their detection in multiple human samples, that TM7 are human commensals that can also be found in clinical settings. Their detection in clinical samples warrants further studies to explore their role in a pathological setting.

Arsenic (As) oxidation by core endosphere microbiome mediates As speciation in Pteris vittata roots

Pteris vittata is an arsenic(As)-hyperaccumulator that may be employed in phytoremediation of As-contaminated soils. P. vittata-associated microbiome are adapted to elevated As and may be important for host survival under stresses. Although P. vittata root endophytes could be critical for As biotransformation in planta, their compositions and metabolisms remain elusive. The current study aims to characterize the root endophytic community composition and As-metabolizing potentials in P. vittata. High As(III) oxidase gene abundances and rapid As(III) oxidation activity indicated that As(III) oxidation was the dominant microbial As-biotransformation processes compared to As reduction and methylization in P. vittata roots. Members of Rhizobiales were the core microbiome and the dominant As(III) oxidizers in P. vittata roots. Acquasition of As-metabolising genes, including both As(III) oxidase and As(V) detoxification reductase genes, through horizontal gene transfer was identified in a Saccharimonadaceae genomic assembly, which was another abundant population residing in P. vittata roots. Acquisition of these genes might improve the fitness of Saccharimonadaceae population to elevated As concentrations in P. vittata. Diverse plant growth promoting traits were encoded by the core root microbiome populations Rhizobiales. We propose that microbial As(III) oxidation and plant growth promotion are critical traits for P. vittata survival in hostile As-contaiminated sites.

A Reduction of Transcriptional Regulation in Aquatic Oligotrophic Microorganisms Enhances Fitness in Nutrient-Poor Environments

In this review, we consider the regulatory strategies of aquatic oligotrophs, microbial cells that are adapted to thrive under low-nutrient concentrations in oceans, lakes, and other aquatic ecosystems. Many reports have concluded that oligotrophs use less transcriptional regulation than copiotrophic cells, which are adapted to high nutrient concentrations and are far more common subjects for laboratory investigations of regulation. It is theorized that oligotrophs have retained alternate mechanisms of regulation, such as riboswitches, that provide shorter response times and smaller amplitude responses and require fewer cellular resources. We examine the accumulated evidence for distinctive regulatory strategies in oligotrophs. We explore differences in the selective pressures copiotrophs and oligotrophs encounter and ask why, although evolutionary history gives copiotrophs and oligotrophs access to the same regulatory mechanisms, they might exhibit distinctly different patterns in how these mechanisms are used. We discuss the implications of these findings for understanding broad patterns in the evolution of microbial regulatory networks and their relationships to environmental niche and life history strategy. We ask whether these observations, which have emerged from a decade of increased investigation of the cell biology of oligotrophs, might be relevant to recent discoveries of many microbial cell lineages in nature that share with oligotrophs the property of reduced genome size.

Effects of macrophytes on micro - And nanoplastic retention and cycling in constructed wetlands

Macrophytes play the important roles in purifying pollutants of constructed wetlands (CWs), while their effects on CWs exposed to micro/nano plastics are not clear. Therefore, planted and unplanted CWs were established to reveal the impacts of macrophytes (Iris pseudacorus) on the overall performance of CWs under polystyrene micro/nano plastics (PS MPs/NPs) exposure. Results showed that macrophytes effectively enhanced the interception capacities of CWs to PS NPs, and significantly promoted the removal of nitrogen and phosphorus after exposed to PS MPs/NPs. Meanwhile, macrophytes improved the activities of dehydrogenase, urease, and phosphatase. Sequencing analysis showed that macrophytes optimized the composition of microbial communities in CWs and stimulated the growth of functional bacteria involved in nitrogen and phosphorus transformation. Moreover, macrophytes further altered the absolute abundance of nitrogen transformation functional genes (amoA, nxrA, narG and nirS). Functional annotation analysis revealed that macrophytes promoted metabolic functions such as Xenobiotics, Amino acids, Lipids metabolism and Signal transduction, ensuring the metabolic balance and homeostasis of microbes under PS MPs/NPs stress. These results exhibited profound implications for the comprehensive evaluation on the roles of macrophytes in CWs for treating wastewater containing PS MPs/NPs.