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.

Aerobic Anoxygenic Phototrophic Bacteria in the Marine Environments Revealed by Raman/Fluorescence-Guided Single-Cell Sorting and Targeted Metagenomics.

Aerobic anoxygenic phototrophic bacteria (AAPB) contribute profoundly to the global carbon cycle. However, most AAPB in marine environments are uncultured and at low abundance, hampering the recognition of their functions and molecular mechanisms. In this study, we developed a new culture-independent method to identify and sort AAPB using single-cell Raman/fluorescence spectroscopy. Characteristic Raman and fluorescent bands specific to bacteriochlorophyll a (Bchl a) in AAPB were determined by comparing multiple known AAPB with non-AAPB isolates. Using these spectroscopic biomarkers, AAPB in coastal seawater, pelagic seawater, and hydrothermal sediment samples were screened, sorted, and sequenced. 16S rRNA gene analysis and functional gene annotations of sorted cells revealed novel AAPB members and functional genes, including one species belonging to the genus Sphingomonas, two genera affiliated to classes Betaproteobacteria and Gammaproteobacteria, and function genes bchCDIX, pucC2, and pufL related to Bchl a biosynthesis and photosynthetic reaction center assembly. Metagenome-assembled genomes (MAGs) of sorted cells from pelagic seawater and deep-sea hydrothermal sediment belonged to Erythrobacter sanguineus that was considered as an AAPB and genus Sphingomonas, respectively. Moreover, multiple photosynthesis-related genes were annotated in both MAGs, and comparative genomic analysis revealed several exclusive genes involved in amino acid and inorganic ion metabolism and transport. This study employed a new single-cell spectroscopy method to detect AAPB, not only broadening the taxonomic and genetic contents of AAPB in marine environments but also revealing their genetic mechanisms at the single-genomic level.

Human genetic associations of the airway microbiome in chronic obstructive pulmonary disease.

Little is known about the relationships between human genetics and the airway microbiome. Deeply sequenced airway metagenomics, by simultaneously characterizing the microbiome and host genetics, provide a unique opportunity to assess the microbiome-host genetic associations. Here we performed a co-profiling of microbiome and host genetics with the identification of over 5 million single nucleotide polymorphisms (SNPs) through deep metagenomic sequencing in sputum of 99 chronic obstructive pulmonary disease (COPD) and 36 healthy individuals. Host genetic variation was the most significant factor associated with the microbiome except for geography and disease status, with its top 5 principal components accounting for 12.11% of the microbiome variability. Within COPD individuals, 113 SNPs mapped to candidate genes reported as genetically associated with COPD exhibited associations with 29 microbial species and 48 functional modules (P < 1 × 10-5), where Streptococcus salivarius exhibits the strongest association to SNP rs6917641 in TBC1D32 (P = 9.54 × 10-8). Integration of concurrent host transcriptomic data identified correlations between the expression of host genes and their genetically-linked microbiome features, including NUDT1, MAD1L1 and Veillonella parvula, TTLL9 and Stenotrophomonas maltophilia, and LTA4H and Haemophilus influenzae. Mendelian randomization analyses revealed a potential causal link between PARK7 expression and microbial type III secretion system, and a genetically-mediated association between COPD and increased relative abundance of airway Streptococcus intermedius. These results suggest a previously underappreciated role of host genetics in shaping the airway microbiome and provide fresh hypotheses for genetic-based host-microbiome interactions in COPD.

Hidden diversity and potential ecological function of phosphorus acquisition genes in widespread terrestrial bacteriophages.

Phosphorus (P) limitation of ecosystem processes is widespread in terrestrial habitats. While a few auxiliary metabolic genes (AMGs) in bacteriophages from aquatic habitats are reported to have the potential to enhance P-acquisition ability of their hosts, little is known about the diversity and potential ecological function of P-acquisition genes encoded by terrestrial bacteriophages. Here, we analyze 333 soil metagenomes from five terrestrial habitat types across China and identify 75 viral operational taxonomic units (vOTUs) that encode 105 P-acquisition AMGs. These AMGs span 17 distinct functional genes involved in four primary processes of microbial P-acquisition. Among them, over 60% (11/17) have not been reported previously. We experimentally verify in-vitro enzymatic activities of two pyrophosphatases and one alkaline phosphatase encoded by P-acquisition vOTUs. Thirty-six percent of the 75 P-acquisition vOTUs are detectable in a published global topsoil metagenome dataset. Further analyses reveal that, under certain circumstances, the identified P-acquisition AMGs have a greater influence on soil P availability and are more dominant in soil metatranscriptomes than their corresponding bacterial genes. Overall, our results reinforce the necessity of incorporating viral contributions into biogeochemical P cycling.

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.

Temperature, pH, and oxygen availability contributed to the functional differentiation of ancient Nitrososphaeria.

Ammonia-oxidizing Nitrososphaeria are among the most abundant archaea on Earth and have profound impacts on the biogeochemical cycles of carbon and nitrogen. In contrast to these well-studied ammonia-oxidizing archaea (AOA), deep-branching non-AOA within this class remain poorly characterized because of a low number of genome representatives. Here, we reconstructed 128 Nitrososphaeria metagenome-assembled genomes from acid mine drainage and hot spring sediment metagenomes. Comparative genomics revealed that extant non-AOA are functionally diverse, with capacity for carbon fixation, carbon monoxide oxidation, methanogenesis, and respiratory pathways including oxygen, nitrate, sulfur, or sulfate, as potential terminal electron acceptors. Despite their diverse anaerobic pathways, evolutionary history inference suggested that the common ancestor of Nitrososphaeria was likely an aerobic thermophile. We further surmise that the functional differentiation of Nitrososphaeria was primarily shaped by oxygen, pH, and temperature, with the acquisition of pathways for carbon, nitrogen, and sulfur metabolism. Our study provides a more holistic and less biased understanding of the diversity, ecology, and deep evolution of the globally abundant Nitrososphaeria.

Strong succession in prokaryotic association networks and community assembly mechanisms in an acid mine drainage-impacted riverine ecosystem.

Acid mine drainage (AMD) serves as an ideal model system for investigating microbial ecology, interaction, and assembly mechanism in natural environments. While previous studies have explored the structure and function of microbial communities in AMD, the succession patterns of microbial association networks and underlying assembly mechanisms during natural attenuation processes remain elusive. Here, we investigated prokaryotic microbial diversity and community assembly along an AMD-impacted river, from the extremely acidic, heavily polluted headwaters to the nearly neutral downstream sites. Microbial diversity was increased along the river, and microbial community composition shifted from acidophile-dominated to freshwater taxa-dominated communities. The complexity and relative modularity of the microbial networks were also increased, indicating greater network stability during succession. Deterministic processes, including abiotic selection of pH and high contents of sulfur and iron, governed community assembly in the headwaters. Although the stochasticity ratio was increased downstream, manganese content, microbial negative cohesion, and relative modularity played important roles in shaping microbial community structure. Overall, this study provides valuable insights into the ecological processes that govern microbial community succession in AMD-impacted riverine ecosystems. These findings have important implications for in-situ remediation of AMD contamination.

In situ enrichment of sulphate-reducing microbial communities with different carbon sources stimulating the acid mine drainage sediments.

The applications of sulphate-reducing microorganisms (SRMs) in acid mine drainage (AMD) treatment systems have received extensive attention due to their ability to reduce sulphate and stabilize metal(loid)s. Despite great phylogenetic diversity of SRMs, only a few have been used in AMD treatment bioreactors. In situ enrichment could be an efficient approach to select new effective SRMs for AMD treatment. Here, we performed in situ enrichment of SRMs in highly stratified AMD sediment cores using different kinds of carbon source mixture. The dsrAB (dissimilatory sulfite reductase) genes affiliated with nine phyla (two archaeal and seven bacterial phyla) and 26 genera were enriched. Remarkably, those genes affiliated with Aciduliprofundum and Vulcanisaeta were enriched in situ in AMD-related environments for the first time, and their relative abundances were negatively correlated with pH. Furthermore, 107 dsrAB-containing metagenome-assembled genomes (MAGs) were recovered from metagenomic datasets, with 14 phyla (two archaeal and 12 bacterial phyla) and 15 genera. The relative abundances of MAGs were positively correlated with total carbon and sulphate contents. Our findings expanded the diversity of SRMs that can be enriched in AMD sediment, and revealed the physiochemical properties that might affect the growth of SRMs, which provided guidance for AMD treatment bioreators.

Long-term elevated precipitation induces grassland soil carbon loss via microbe-plant-soil interplay.

Global climate models predict that the frequency and intensity of precipitation events will increase in many regions across the world. However, the biosphere-climate feedback to elevated precipitation (eP) remains elusive. Here, we report a study on one of the longest field experiments assessing the effects of eP, alone or in combination with other climate change drivers such as elevated CO2 (eCO2 ), warming and nitrogen deposition. Soil total carbon (C) decreased after a decade of eP treatment, while plant root production decreased after 2 years. To explain this asynchrony, we found that the relative abundances of fungal genes associated with chitin and protein degradation increased and were positively correlated with bacteriophage genes, suggesting a potential viral shunt in C degradation. In addition, eP increased the relative abundances of microbial stress tolerance genes, which are essential for coping with environmental stressors. Microbial responses to eP were phylogenetically conserved. The effects of eP on soil total C, root production, and microbes were interactively affected by eCO2 . Collectively, we demonstrate that long-term eP induces soil C loss, owing to changes in microbial community composition, functional traits, root production, and soil moisture. Our study unveils an important, previously unknown biosphere-climate feedback in Mediterranean-type water-limited ecosystems, namely how eP induces soil C loss via microbe-plant-soil interplay.

The facilitating role of phycospheric heterotrophic bacteria in cyanobacterial phosphonate availability and Microcystis bloom maintenance.

BACKGROUND: Phosphonates are the main components in the global phosphorus redox cycle. Little is known about phosphonate metabolism in freshwater ecosystems, although rapid consumption of phosphonates has been observed frequently. Cyanobacteria are often the dominant primary producers in freshwaters; yet, only a few strains of cyanobacteria encode phosphonate-degrading (C-P lyase) gene clusters. The phycosphere is defined as the microenvironment in which extensive phytoplankton and heterotrophic bacteria interactions occur. It has been demonstrated that phytoplankton may recruit phycospheric bacteria based on their own needs. Therefore, the establishment of a phycospheric community rich in phosphonate-degrading-bacteria likely facilitates cyanobacterial proliferation, especially in waters with scarce phosphorus. We characterized the distribution of heterotrophic phosphonate-degrading bacteria in field Microcystis bloom samples and in laboratory cyanobacteria "phycospheres" by qPCR and metagenomic analyses. The role of phosphonate-degrading phycospheric bacteria in cyanobacterial proliferation was determined through coculturing of heterotrophic bacteria with an axenic Microcystis aeruginosa strain and by metatranscriptomic analysis using field Microcystis aggregate samples. RESULTS: Abundant bacteria that carry C-P lyase clusters were identified in plankton samples from freshwater Lakes Dianchi and Taihu during Microcystis bloom periods. Metagenomic analysis of 162 non-axenic laboratory strains of cyanobacteria (consortia cultures containing heterotrophic bacteria) showed that 20% (128/647) of high-quality bins from eighty of these consortia encode intact C-P lyase clusters, with an abundance ranging up to nearly 13%. Phycospheric bacterial phosphonate catabolism genes were expressed continually across bloom seasons, as demonstrated through metatranscriptomic analysis using sixteen field Microcystis aggregate samples. Coculturing experiments revealed that although Microcystis cultures did not catabolize methylphosphonate when axenic, they demonstrated sustained growth when cocultured with phosphonate-utilizing phycospheric bacteria in medium containing methylphosphonate as the sole source of phosphorus. CONCLUSIONS: The recruitment of heterotrophic phosphonate-degrading phycospheric bacteria by cyanobacteria is a hedge against phosphorus scarcity by facilitating phosphonate availability. Cyanobacterial consortia are likely primary contributors to aquatic phosphonate mineralization, thereby facilitating sustained cyanobacterial growth, and even bloom maintenance, in phosphate-deficient waters. Video Abstract.

The airway microbiome mediates the interaction between environmental exposure and respiratory health in humans.

Exposure to environmental pollution influences respiratory health. The role of the airway microbial ecosystem underlying the interaction of exposure and respiratory health remains unclear. Here, through a province-wide chronic obstructive pulmonary disease surveillance program, we conducted a population-based survey of bacterial (n = 1,651) and fungal (n = 719) taxa and metagenomes (n = 1,128) from induced sputum of 1,651 household members in Guangdong, China. We found that cigarette smoking and higher PM2.5 concentration were associated with lung function impairment through the mediation of bacterial and fungal communities, respectively, and that exposure was associated with an enhanced inter-kingdom microbial interaction resembling the pattern seen in chronic obstructive pulmonary disease. Enrichment of Neisseria was associated with a 2.25-fold increased risk of high respiratory symptom burden, coupled with an elevation in Aspergillus, in association with occupational pollution. We developed an individualized microbiome-based health index, which covaried with exposure, respiratory symptoms and diseases, with potential generalizability to global datasets. Our results may inform environmental risk prevention and guide interventions that harness airway microbiome.

Generation of zero-valent sulfur from dissimilatory sulfate reduction in sulfate-reducing microorganisms.

Dissimilatory sulfate reduction (DSR) mediated by sulfate-reducing microorganisms (SRMs) plays a pivotal role in global sulfur, carbon, oxygen, and iron cycles since at least 3.5 billion y ago. The canonical DSR pathway is believed to be sulfate reduction to sulfide. Herein, we report a DSR pathway in phylogenetically diverse SRMs through which zero-valent sulfur (ZVS) is directly generated. We identified that approximately 9% of sulfate reduction was directed toward ZVS with S8 as a predominant product, and the ratio of sulfate-to-ZVS could be changed with SRMs' growth conditions, particularly the medium salinity. Further coculturing experiments and metadata analyses revealed that DSR-derived ZVS supported the growth of various ZVS-metabolizing microorganisms, highlighting this pathway as an essential component of the sulfur biogeochemical cycle.

Airway dysbiosis accelerates lung function decline in chronic obstructive pulmonary disease.

Progressive lung function decline is a hallmark of chronic obstructive pulmonary disease (COPD). Airway dysbiosis occurs in COPD, but whether it contributes to disease progression remains unknown. Here, we show, through a longitudinal analysis of two cohorts involving four UK centers, that baseline airway dysbiosis in COPD patients, characterized by the enrichment of opportunistic pathogenic taxa, associates with a rapid forced expiratory volume in 1 s (FEV1) decline over 2 years. Dysbiosis associates with exacerbation-related FEV1 fall and sudden FEV1 fall at stability, contributing to long-term FEV1 decline. A third cohort in China further validates the microbiota-FEV1-decline association. Human multi-omics and murine studies show that airway Staphylococcus aureus colonization promotes lung function decline through homocysteine, which elicits a neutrophil apoptosis-to-NETosis shift via the AKT1-S100A8/A9 axis. S. aureus depletion via bacteriophages restores lung function in emphysema mice, providing a fresh approach to slow COPD progression by targeting the airway microbiome.

The contrasting effects of fluctuating temperature on bacterial diversity and performances in temperate and subtropical soils.

The extremely high species diversity of soil bacterial community has fascinated and puzzled community ecologists. Although theory predicts that fluctuations in environments can facilitate diversity maintenance, the effects of fluctuating temperature on species diversity have rarely been investigated in species-rich microbial communities. Here, we examined whether fluctuating temperature had positive effects on species diversity relative to constant temperatures in soil bacterial communities, and investigated the effects of fluctuating temperature on bacterial performances (changes in relative abundance). We performed a temperature manipulation experiment with soils collected from temperate and subtropical zones, where the soils were subjected to constant high, low or fluctuating temperatures. We found that fluctuating temperatures showed significant positive effects on species diversity. The time-averaged effect of fluctuating temperatures (i.e., averaging out the differences between species in their environment-dependent performances) appeared to delay species loss in both the temperate and the subtropical communities. In addition, we found that the performances of temperature-responsive species at fluctuating temperatures significantly deviated from their time-weighted average performances at constant high and low temperatures, which was defined as fluctuation-dependent effects in our study. Intriguingly, fluctuation-dependent effects beyond time-averaged effect led to an opposite trend: differences in temperature-responsive species' performances decreased in the temperate communities, but increased in the subtropical communities. Our findings provide new insights into diversity maintenance in soil bacterial communities, and imply that the effects of fluctuating temperature on species diversity in soil bacterial community might vary across latitude.

Metagenomic Discovery of "Candidatus Parvarchaeales"-Related Lineages Sheds Light on Adaptation and Diversification from Neutral-Thermal to Acidic-Mesothermal Environments.

"Candidatus Parvarchaeales" microbes, representing a DPANN archaeal group with limited metabolic potential and reliance on hosts for their growth, were initially found in acid mine drainage (AMD). Due to the lack of representatives, however, their ecological roles and adaptation to extreme habitats such as AMD as well as how they diverge across the lineage remain largely unexplored. By applying genome-resolved metagenomics, 28 Parvarchaeales-associated metagenome-assembled genomes (MAGs) representing two orders and five genera were recovered. Among them, we identified three new genera and proposed the names "Candidatus Jingweiarchaeum," "Candidatus Haiyanarchaeum," and "Candidatus Rehaiarchaeum," with the former two belonging to a new order, "Candidatus Jingweiarchaeales." Further analyses of the metabolic potentials revealed substantial niche differentiation between Jingweiarchaeales and Parvarchaeales. Jingweiarchaeales may rely on fermentation, salvage pathways, partial glycolysis, and the pentose phosphate pathway (PPP) for energy conservation reservation, while the metabolic potentials of Parvarchaeales might be more versatile. Comparative genomic analyses suggested that Jingweiarchaeales favor habitats with higher temperatures and that Parvarchaeales are better adapted to acidic environments. We further revealed that the thermal adaptation of these lineages, especially Haiyanarchaeum, might rely on genomic features such as the usage of specific amino acids, genome streamlining, and hyperthermophile featured genes such as rgy. Notably, the adaptation of Parvarchaeales to acidic environments was possibly driven by horizontal gene transfer (HGT). The reconstruction of ancestral states demonstrated that both may have originated from thermal and neutral environments and later spread to mesothermal and acidic environments. These evolutionary processes may also be accompanied by adaptation to oxygen-rich environments via HGT. IMPORTANCE "Candidatus Parvarchaeales" microbes may represent a lineage uniquely distributed in extreme environments such as AMD and hot springs. However, little is known about the strategies and processes of how they adapted to these extreme environments. By the discovery of potential new order-level lineages, "Ca. Jingweiarchaeales," and in-depth comparative genomic analysis, we unveiled the functional differentiation of these lineages. Furthermore, we show that the adaptation of these lineages to high-temperature and acidic environments was driven by different strategies, with the former relying more on genomic characteristics such as genome streamlining and amino acid compositions and the latter relying more on the acquisition of genes associated with acid tolerance. Finally, by the reconstruction of the ancestral states of the optimal growth temperature (OGT) and isoelectric point (pI), we showed the potential evolutionary process of Parvarchaeales-related lineages with regard to the shift from the high-temperature environment of their common ancestors to low-temperature (potentially acidic) environments.

Panguiarchaeum symbiosum, a potential hyperthermophilic symbiont in the TACK superphylum.

The biology of Korarchaeia remains elusive due to the lack of genome representatives. Here, we reconstruct 10 closely related metagenome-assembled genomes from hot spring habitats and place them into a single species, proposed herein as Panguiarchaeum symbiosum. Functional investigation suggests that Panguiarchaeum symbiosum is strictly anaerobic and grows exclusively in thermal habitats by fermenting peptides coupled with sulfide and hydrogen production to dispose of electrons. Due to its inability to biosynthesize archaeal membranes, amino acids, and purines, this species likely exists in a symbiotic lifestyle similar to DPANN archaea. Population metagenomics and metatranscriptomic analyses demonstrated that genes associated with amino acid/peptide uptake and cell attachment exhibited positive selection and were highly expressed, supporting the proposed proteolytic catabolism and symbiotic lifestyle. Our study sheds light on the metabolism, evolution, and potential symbiotic lifestyle of Panguiarchaeum symbiosum, which may be a unique host-dependent archaeon within the TACK superphylum.

Diverse Methylmercury (MeHg) Producers and Degraders Inhabit Acid Mine Drainage Sediments, but Few Taxa Correlate with MeHg Accumulation.

Methylmercury (MeHg) is a notorious neurotoxin, and its production and degradation in the environment are mainly driven by microorganisms. A variety of microbial MeHg producers carrying the gene pair hgcAB and degraders carrying the merB gene have been separately reported in recent studies. However, surprisingly little attention has been paid to the simultaneous investigation of the diversities of microbial MeHg producers and degraders in a given habitat, and no studies have been performed to explore to what extent these two contrasting microbial groups correlate with MeHg accumulation in the habitat of interest. Here, we collected 86 acid mine drainage (AMD) sediments from an area spanning approximately 500,000 km2 in southern China and profiled the sediment-borne putative MeHg producers and degraders using genome-resolved metagenomics. 46 metagenome-assembled genomes (MAGs) containing hgcAB and 93 MAGs containing merB were obtained, including those from various taxa without previously known MeHg-metabolizing microorganisms. These diverse MeHg-metabolizing MAGs were formed largely via multiple independent horizontal gene transfer (HGT) events. The putative MeHg producers from Deltaproteobacteria and Firmicutes as well as MeHg degraders from Acidithiobacillia were closely correlated with MeHg accumulation in the sediments. Furthermore, these three taxa, in combination with two abiotic factors, explained over 60% of the variance in MeHg accumulation. Most of the members of these taxa were characterized by their metabolic potential for nitrogen fixation and copper tolerance. Overall, these findings improve our understanding of the ecology of MeHg-metabolizing microorganisms and likely have implications for the development of management strategies for the reduction of MeHg accumulation in the AMD sediments. IMPORTANCE Microorganisms are the main drivers of MeHg production and degradation in the environment. However, little attention has been paid to the simultaneous investigation of the diversities of microbial MeHg producers and degraders in a given habitat. We used genome-resolved metagenomics to reveal the vast phylogenetic and metabolic diversities of putative MeHg producers and degraders in AMD sediments. Our results show that the diversity of MeHg-metabolizing microorganisms (particularly MeHg degraders) in AMD sediments is much higher than was previously recognized. Via multiple linear regression analysis, we identified both microbial and abiotic factors affecting MeHg accumulation in AMD sediments. Despite their great diversity, only a few taxa of MeHg-metabolizing microorganisms were closely correlated with MeHg accumulation. This work underscores the importance of using genome-resolved metagenomics to survey MeHg-metabolizing microorganisms and provides a framework for the illumination of the microbial basis of MeHg accumulation via the characterization of physicochemical properties, MeHg-metabolizing microorganisms, and the correlations between them.

Genomic Legacies of Ancient Adaptation Illuminate GC-Content Evolution in Bacteria.

Bacterial evolution is characterized by strong purifying selection as well as rapid adaptive evolution in changing environments. In this context, the genomic GC content (genomic GC) varies greatly but presents some level of phylogenetic stability, making it challenging to explain based on current hypotheses. To illuminate the evolutionary mechanisms of the genomic GC, we analyzed the base composition and functional inventory of 11,083 representative genomes. A phylogenetically constrained bimodal distribution of the genomic GC, which mainly originated from parallel divergences in the early evolution, was demonstrated. Such variation of the genomic GC can be well explained by DNA replication and repair (DRR), in which multiple pathways correlate with the genomic GC. Furthermore, the biased conservation of various stress-related genes, especially the DRR-related ones, implies distinct adaptive processes in the ancestral lineages of high- or low-GC clades which are likely induced by major environmental changes. Our findings support that the mutational biases resulting from these legacies of ancient adaptation have changed the course of adaptive evolution and generated great variation in the genomic GC. This highlights the importance of indirect effects of natural selection, which indicates a new model for bacterial evolution. IMPORTANCE GC content has been shown to be an important factor in microbial ecology and evolution, and the genomic GC of bacteria can be characterized by great intergenomic heterogeneity, high intragenomic homogeneity, and strong phylogenetic inertia, as well as being associated with the environment. Current hypotheses concerning direct selection or mutational biases cannot well explain these features simultaneously. Our findings of the genomic GC showing that ancient adaptations have transformed the DRR system and that the resulting mutational biases further contributed to a bimodal distribution of it offer a more reasonable scenario for the mechanism. This would imply that, when thinking about the evolution of life, diverse processes of adaptation exist, and combined effects of natural selection should be considered.

Phytostabilization mitigates antibiotic resistance gene enrichment in a copper mine tailings pond.

Mining-impacted environments are distributed globally and have become increasingly recognized as hotspots of antibiotic resistance genes (ARGs). However, there are currently no reports on treatment technologies to deal with such an important environmental problem. To narrow this knowledge gap, we implemented a phytostabilization project in an acidic copper mine tailings pond and employed metagenomics to explore ARG characteristics in the soil samples. Our results showed that phytostabilization decreased the total ARG abundance in 0-10 cm soil layer by 75 %, which was companied by a significant decrease in ARG mobility, and a significant increase in ARG diversity and microbial diversity. Phytostabilization was also found to drastically alter the ARG host composition and to significantly reduce the total abundance of virulence factor genes of ARG hosts. Soil nutrient status, heavy metal toxicity and SO42- concentration were important physicochemical factors to affect the total ARG abundance, while causal mediation analysis showed that their effects were largely mediated by the changes in ARG mobility and microbial diversity. The increase in ARG diversity associated with phytostabilization was mainly mediated by a small subgroup of ARG hosts, most of which could not be classified at the genus level and deserve further research in the future.

Functional Guilds, Community Assembly, and Co-occurrence Patterns of Fungi in Metalliferous Mine Tailings Ponds in Mainland China.

Metalliferous mine tailings ponds are generally characterized by low levels of nutrient elements, sustained acidic conditions, and high contents of toxic metals. They represent one kind of extreme environments that are believed to resemble the Earth's early environmental conditions. There is increasing evidence that the diversity of fungi inhabiting mine tailings ponds is much higher than previously thought. However, little is known about functional guilds, community assembly, and co-occurrence patterns of fungi in such habitats. As a first attempt to address this critical knowledge gap, we employed high-throughput sequencing to characterize fungal communities in 33 mine tailings ponds distributed across 18 provinces of mainland China. A total of 5842 fungal phylotypes were identified, with saprotrophic fungi being the major functional guild. The predictors of fungal diversity in whole community and sub-communities differed considerably. Community assembly of the whole fungal community and individual functional guilds were primarily governed by stochastic processes. Total soil nitrogen and total phosphorus mediated the balance between stochastic and deterministic processes of the fungal community assembly. Co-occurrence network analysis uncovered a high modularity of the whole fungal community. The observed main modules largely consisted of saprotrophic fungi as well as various phylotypes that could not be assigned to known functional guilds. The richness of core fungal phylotypes, occupying vital positions in co-occurrence network, was positively correlated with edaphic properties such as soil enzyme activity. This indicates the important roles of core fungal phylotypes in soil organic matter decomposition and nutrient cycling. These findings improve our understanding of fungal ecology of extreme environments.