Search

Microbes are key mediators in maintaining environmental stability and the migration of geological elements. Although the effect of environmental factors on the geochemical distribution of elements and microbes in Xikuangshan have been studied, the connection between the above two remains to be clarified. In this study, the distribution of elemental forms, especially heavy metals, various geochemical conditions, and the microbial community compositions of the soil, sediment and water samples collected from different sites of Xikuangshan were investigated to study the correlation between representative elemental geochemistry and microbial diversity. The results show that the pH values of the soil samples were significantly lower than those of the sediment and water samples, and the soil had the highest microbial biomass. The contents of As and Sb in the samples are higher than those of other heavy metal elements, and the most abundant occurrence form of heavy metals is reducible. The most abundant type of the microorganisms is bacteria, followed by fungi and archaea. The correlation analyses between environmental factors and microbial diversities indicated that the microbial community structure composition of Xikuangshan was strictly correlated with the content and speciation of elements. The main factors impacting microbial diversity are pH and redox potential, as well as the content and occurrence forms of Sb, As and Fe. This study is of value for revealing the geochemical characteristics of heavy metal occurrence forms associated with microbial diversity in the Xikuangshan mine and for pollution control of Sb mines.

Elton's biotic resistance hypothesis, which posits that diverse communities should be more resistant to biological invasions, has received considerable experimental support. However, it remains unclear whether such a negative diversity-invasibility relationship would persist under anthropogenic environmental change. By using the common ragweed (Ambrosia artemisiifolia) as a model invader, our 4-year grassland experiment demonstrated consistently negative relationships between resident species diversity and community invasibility, irrespective of nitrogen addition, a result further supported by a meta-analysis. Importantly, our experiment showed that plant diversity consistently resisted invasion simultaneously through increased resident biomass, increased trait dissimilarity among residents, and increased community-weighted means of resource-conservative traits that strongly resist invasion, pointing to the importance of both trait complementarity and sampling effects for invasion resistance even under resource enrichment. Our study provides unique evidence that considering species' functional traits can help further our understanding of biotic resistance to biological invasions in a changing environment.

A wide array of microorganisms, including many novel, phylogenetically deeply rooted taxa, survive and thrive in extreme environments. These unique and reduced-complexity ecosystems offer a tremendous opportunity for studying the structure, function and evolution of natural microbial communities. Marker gene surveys have resolved patterns and ecological drivers of these extremophile assemblages, revealing a vast uncultured microbial diversity and the often predominance of archaea in the most extreme conditions. New omics studies have uncovered linkages between community function and environmental variables, and have enabled discovery and genomic characterization of major new lineages that substantially expand microbial diversity and change the structure of the tree of life. These efforts have significantly advanced our understanding of the diversity, ecology and evolution of microorganisms populating Earth’s extreme environments, and have facilitated the exploration of microbiota and processes in more complex ecosystems. Microbial life can thrive in extreme environments such as terrestrial hot springs and deep sea hydrothermal vents, glaciers and permafrost, hypersaline habitats, acid mine drainage and the subsurface. In this Review, Shu and Huang explore the diversity, functions and evolution of bacteria and archaea inhabiting Earth’s major extreme environments.

CandidatusParvarchaeales, representing a DPANN archaeal group with limited metabolic potentials 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, 28Parvarchaeales-associated metagenome-assembled genomes (MAGs) representing two orders and five genera were recovered. Among them, we identified three new genera and proposed the namesCandidatusJingweiarchaeum,CandidatusHaiyanarchaeum, andCandidatusRehaiarchaeum with the former two belonging to a new orderCandidatusJingweiarchaeales. Further analyses of metabolic potentials revealed substantial niche differentiation between Jingweiarchaeales and Parvarchaeales. Jingweiarchaeales may rely on fermentation, salvage pathways, partial glycolysis, and pentose phosphate pathway (PPP) for energy reservation, while the metabolic potentials of Parvarchaeales might be more versatile. Comparative genomic analyses suggested that Jingweiarchaeales are more favorable to habitats with higher temperatures andParvarchaealesare better adapted to acidic environments. We further revealed that the thermal adaptation of these lineages especially for Haiyanarchaeum might rely on innate genomic features such as the usage of specific amino acids, genome streamlining, and hyperthermal featured genes such asrgy. Notably, the acidic adaptation of Parvarchaeales was possibly driven by horizontal gene transfer (HGT). Reconstruction of ancestral states demonstrated that both may originate from thermal and neutral environments and later spread to mesothermal and acidic environments. These evolutionary processes may also be accompanied by adaptation toward oxygen-rich environments via HGT.ImportanceCandidatusParvarchaeales 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 - Jingweiarchaeales and in-depth comparative genomic analysis, we unveiled the functional differentiation of these lineages. Further, we show that the adaptation to high-temperature and acidic environments of these lineages was driven by different strategies, with the prior relying more on innate genomic characteristics and the latter more on the acquisition of genes associated with acid tolerance. Finally, by reconstruction of ancestral states of OGT andpI, we showed the potential evolutionary process of Parvarchaeales-related lineages with regard to the shift from a high-temperature environment of their common ancestors to low-temperature (potentially acidic) environments.

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.

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.

Understanding the influences of global climate change on soil microbial communities is essential in evaluating the terrestrial biosphere's feedback to this alarming anthropogenic disturbance. However, little is known about how intra-site historical climate variability can mediate the influences of current climate differences on community dissimilarity and assembly. To fill this gap, we examined and disentangled the interactive effects of historical climate variability and current climate differences on the soil bacterial community dissimilarity and stochasticity of community assembly among 143 sites from 28 forests across eastern China. We hypothesize that the relative importance of stochasticity and community dissimilarity are related to historical climate variability and that an increasing sum of intra-site historical variability enhances stochasticity while reduces dissimilarity between two communities. To test our hypothesis, we statistically controlled for covariates between sites including differences in soil chemistry, plant diversity, spatial distance, and seasonal climate variations at annual timescales. We observed that an increase in inter-site current climate differences led to a reduced impact of stochasticity in community assembly and a pronounced divergence between communities. In stark contrast, when communities were subjected to a high level of intra-site historical climate fluctuation, the observed impact incurred from current climate differences was substantially weakened. Moreover, the influence of increased historical variability was consistent along the gradient of current temperature differences between sites. However, effects induced by historical fluctuation in precipitation were disproportional and only evident when small inter-site differences were observed. Consequently, if the prior climate variability is ignored, especially regarding environmental factors like temperature, we assert that the influence current climate differentiation has on regulating community dissimilarity and assembly stochasticity will be underestimated. Together, our findings highlight the importance and need of explicitly controlling the mean and the historical variability of climate factors for the next "generation" of climate change experiments to come.

The haloalkalitolerant bacterium Egicoccus halophilus EGI 80432T exhibits high adaptability to saline–alkaline environment. The salinity adaptation mechanism of E. halophilus EGI 80432T was fully understood based on transcriptome analyses and physiological responses; however, the alkaline response mechanism has not yet been investigated. Here, we investigated the alkaline response mechanism of E. halophilus EGI 80432T by a transcriptomic comparison. In this study, the genes involved in the glycolysis, TCA cycle, starch, and trehalose metabolism for energy production and storage, were up-regulated under highly alkaline condition. Furthermore, genes responsible for the production of acidic and neutral metabolites, i.e., acetate, pyruvate, formate, glutamate, threonine, and ectoine, showed increased expression under highly alkaline condition, compared with the control pH condition. In contrast, the opposite results were observed in proton capture or retention gene expression profiles, i.e., cation/proton antiporters and ATP synthases. The above results revealed that E. halophilus EGI 80432T likely tended to adopt an “acidic metabolites production” strategy in response to a highly alkaline condition. These findings would pave the way for further studies in the saline–alkaline adaptation mechanisms of E. halophilus EGI 80432T, and hopefully provide a new insight into the foundational theory and application in ecological restoration with saline–alkaline strains.

Biological nitrogen fixation (BNF) by cyanobacteria is of significant importance for the Earth’s biogeochemical nitrogen cycle but is restricted to a few genera that do not form monophyletic group. To explore the evolutionary trajectory of BNF and investigate the driving forces of its evolution, we analyze 650 cyanobacterial genomes and compile the database of diazotrophic cyanobacteria based on the presence of nitrogen fixation gene clusters (NFGCs). We report that 266 of 650 examined genomes are NFGC-carrying members, and these potentially diazotrophic cyanobacteria are unevenly distributed across the phylogeny of Cyanobacteria, that multiple independent losses shaped the scattered distribution. Among the diazotrophic cyanobacteria, two types of NFGC exist, with one being ancestral and abundant, which have descended from diazotrophic ancestors, and the other being anaerobe-like and sparse, possibly being acquired from anaerobic microbes through horizontal gene transfer. Interestingly, we illustrate that the origin of BNF in Cyanobacteria coincide with two major evolutionary events. One is the origin of multicellularity of cyanobacteria, and the other is concurrent genetic innovations with massive gene gains and expansions, implicating their key roles in triggering the evolutionary transition from nondiazotrophic to diazotrophic cyanobacteria. Additionally, we reveal that genes involved in accelerating respiratory electron transport (coxABC), anoxygenic photosynthetic electron transport (sqr), as well as anaerobic metabolisms (pfor, hemN, nrdG, adhE) are enriched in diazotrophic cyanobacteria, representing adaptive genetic signatures that underpin the diazotrophic lifestyle. Collectively, our study suggests that multicellularity, together with concurrent genetic adaptations contribute to the evolution of diazotrophic cyanobacteria.

Egicoccus halophilus EGI 80432T, a halotolerant bacterium isolated from a saline-alkaline soil, belongs to a member of the class Nitriliruptoria, which exhibits high adaptability to salt environments. At present, the detailed knowledge of the salinity adaptation strategies of Nitriliruptoria was limited except for one research by using comparative genomics analysis. Here, we investigated the salinity adaptation mechanism of E. halophilus EGI 80432T by comparative physiological and transcriptomic analyses. The results of physiological analyses showed that trehalose and glutamate were accumulated by salt stress and showed the maximum at moderate salinity condition. Furthermore, the contents of histidine, threonine, proline, and ectoine were increased with increasing salt concentration. We found that both 0% and 9% NaCl conditions resulted in increased expressions of genes involved in carbohydrate and energy metabolisms, but negatively affected the Na+ efflux, iron, and molybdate transport. Moreover, the high salt condition led to enhancement of transcription of genes required for the synthesis of compatible solutes, e.g., glutamate, histidine, threonine, proline, and ectoine, which agree with the results of physiological analyses. The above results revealed that E. halophilus EGI 80432T increased inorganic ions uptake and accumulated trehalose and glutamate in response to moderate salinity condition, while the salinity adaptation strategy was changed from a "salt-in-cytoplasm" strategy to a "compatible solute" strategy under high salinity condition. The findings in this study would promote further studies in salt tolerance molecular mechanism of Nitriliruptoria and provide a theoretical support for E. halophilus EGI 80432T's application in ecological restoration.Key Points• Salt stress affected gene expressions responsible for carbohydrate and energy metabolisms of E. halophilus EGI 8042T.• E. halophilus EGI 80432T significantly accumulated compatible solutes under salt stress.• E. halophilus EGI 80432T adopted a "compatible solute" strategy to withstand high salt stress.

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 SO

Mining is among the human activities with widest environmental impacts, and mining-impacted environments are characterized by high levels of metals that can co-select for antibiotic resistance genes (ARGs) in microorganisms. However, ARGs in mining-impacted environments are still poorly understood. Here, we conducted a comprehensive study of ARGs in such environments worldwide, taking advantage of 272 metagenomes generated from a global-scale data collection and two national sampling efforts in China. The average total abundance of the ARGs in globally distributed studied mine sites was 1572 times per gigabase, being rivaling that of urban sewage but much higher than that of freshwater sediments. Multidrug resistance genes accounted for 40% of the total ARG abundance, tended to co-occur with multimetal resistance genes, and were highly mobile (e.g. on average 16% occurring on plasmids). Among the 1848 high-quality metagenome-assembled genomes (MAGs), 85% carried at least one multidrug resistance gene plus one multimetal resistance gene. These high-quality ARG-carrying MAGs considerably expanded the phylogenetic diversity of ARG hosts, providing the first representatives of ARG-carrying MAGs for the Archaea domain and three bacterial phyla. Moreover, 54 high-quality ARG-carrying MAGs were identified as potential pathogens. Our findings suggest that mining-impacted environments worldwide are underexplored hotspots of multidrug resistance genes.

Cyanobacteria are photosynthetic prokaryotes that inhabit diverse aquatic and terrestrial environments. However, the evolutionary mechanisms involved in the cyanobacterial habitat adaptation remain poorly understood. Here, based on phylogenetic and comparative genomic analyses of 650 cyanobacterial genomes, we investigated the genetic basis of cyanobacterial habitat adaptation (marine, freshwater, and terrestrial). We show: (1) the expansion of gene families is a common strategy whereby terrestrial cyanobacteria cope with fluctuating environments, whereas the genomes of many marine strains have undergone contraction to adapt to nutrient-poor conditions. (2) Hundreds of genes are strongly associated with specific habitats. Genes that are differentially abundant in genomes of marine, freshwater, and terrestrial cyanobacteria were found to be involved in light sensing and absorption, chemotaxis, nutrient transporters, responses to osmotic stress, etc., indicating the importance of these genes in the survival and adaptation of organisms in specific habitats. (3) A substantial fraction of genes that facilitate the adaptation of Cyanobacteria to specific habitats are contributed by horizontal gene transfer, and such genetic exchanges are more frequent in terrestrial cyanobacteria. Collectively, our results further our understandings of the adaptations of Cyanobacteria to different environments, highlighting the importance of ecological constraints imposed by the environment in shaping the evolution of Cyanobacteria.

The interaction between airway microbiome and host in chronic obstructive pulmonary disease (COPD) is poorly understood. Here we used a multi-omic meta-analysis approach to characterize the functional signature of airway microbiome in COPD. We retrieved all public COPD sputum microbiome datasets, totaling 1640 samples from 16S rRNA gene datasets and 26 samples from metagenomic datasets from across the world. We identified microbial taxonomic shifts using random effect meta-analysis and established a global classifier for COPD using 12 microbial genera. We inferred the metabolic potentials for the airway microbiome, established their molecular links to host targets, and explored their effects in a separate meta-analysis on 1340 public human airway transcriptome samples for COPD. 29.6% of differentially expressed human pathways were predicted to be targeted by microbiome metabolism. For inferred metabolite-host interactions, the flux of disease-modifying metabolites as predicted from host transcriptome was generally concordant with their predicted metabolic turnover in microbiome, suggesting a synergistic response between microbiome and host in COPD. The meta-analysis results were further validated by a pilot multi-omic study on 18 COPD patients and 10 controls, in which airway metagenome, metabolome, and host transcriptome were simultaneously characterized. 69.9% of the proposed "microbiome-metabolite-host" interaction links were validated in the independent multi-omic data. Butyrate, homocysteine, and palmitate were the microbial metabolites showing strongest interactions with COPD-associated host genes. Our meta-analysis uncovered functional properties of airway microbiome that interacted with COPD host gene signatures, and demonstrated the possibility of leveraging public multi-omic data to interrogate disease biology.

Microbes, similar to plants and animals, exhibit biogeographic patterns. However, in contrast with the considerable knowledge on the island biogeography of higher organisms, we know little about the distribution of microorganisms within and among islands. Here, we explored insular soil bacterial and fungal biogeography and underlying mechanisms, using soil microbiota from a group of land-bridge islands as a model system. Similar to island species-area relationships observed for many macroorganisms, both island-scale bacterial and fungal diversity increased with island area; neither diversity, however, was affected by island isolation. By contrast, bacterial and fungal communities exhibited strikingly different assembly patterns within islands. The loss of bacterial diversity on smaller islands was driven primarily by the systematic decline of diversity within samples, whereas the loss of fungal diversity on smaller islands was driven primarily by the homogenization of community composition among samples. Lower soil moisture limited within-sample bacterial diversity, whereas smaller spatial distances among samples restricted among-sample fungal diversity, on smaller islands. These results indicate that among-island differences in habitat quality generate the bacterial island species-area relationship, whereas within-island dispersal limitation generates the fungal island species-area relationship. Together, our study suggests that different mechanisms underlie similar island biogeography patterns of soil bacteria and fungi.

The group Nitriliruptoria, recently classified as a separate class of phylum Actinobacteria, has five members at present, which belong to halophilic or halotolerant Actinobacteria. Here, we sequenced the genomes of Egicoccus halophilus EGI 80432T and Egibacter rhizosphaerae EGI 80759T, and performed a comparative genomics approach to analyze the genomic differences and salt adaptation mechanisms in Nitriliruptoria. Phylogenetic analysis suggested that Euzebya tangerina F10T has a closer phylogenetic relationship to Euzebya rosea DSW09T, while genomic analysis revealed highest genomic similarity with Nitriliruptor alkaliphilus ANL-iso2T and E. halophilus EGI 80432T. Genomic differences of Nitriliruptoria were mainly observed in genome size, gene contents, and the amounts of gene in per functional categories. Furthermore, our analysis also revealed that Nitriliruptoria possess similar synthesis systems of solutes, such as trehalose, glutamine, glutamate, and proline. On the other hand, each member of Nitriliruptoria species possesses specific mechanisms, K+ influx and efflux, betaine and ectoine synthesis, and compatible solutes transport to survive in various high-salt environments.

It is well known that speciation transformations of As(III) vs. As(V) in acid mine drainage (AMD) are mainly driven by microbially mediated redox reactions of Fe and S. However, these processes are rarely investigated. In this study, columns containing mine water were inoculated with two typical acidophilic Fe/S-oxidizing/reducing bacteria [the chemoautotrophic

Plant roots harbor and interact with diverse fungal species. By changing these belowground fungal communities, focal plants can affect the performance of surrounding individuals and the outcome of coexistence. Although highly host related, the roles of these root-associated fungal communities per se in host plant spatial co-occurrence is largely unknown. Here, we evaluated the host dependency of root-associated communities for 39-plant species spatially mapped throughout a 50-ha subtropical forest plot with relevant environmental properties. In addition, we explored whether the differentiation in root fungal associations among plant species can reflect their observed co-occurrence patterns. We demonstrated a strong host-dependency by discriminating the differentiation of root-associated fungal communities regardless of background soil heterogeneity. Furthermore, Random Forest modeling indicated that these nonrandom root fungal associations significantly increased our ability to explain spatial co-occurrence patterns, and to a greater degree than the relative abundance, phylogenetic relatedness, and functional traits of the host plants. Our results further suggested that plants harbor more abundant shared, “generalist” pathogens are likely segregated, while hosting more abundant unique, “specialist” ectomycorrhizal fungi might be an important strategy for promoting spatial aggregation, particularly between early established trees and the heterospecific adults. Together, we provide a conceptual and testable approach to integrate this host-dependent root fungal “fingerprinting” into the plant diversity patterns. We highlight that this approach is complementary to the classic cultivation-based scheme and can deepen our understanding of the community-level effect from overall fungi and its contribution to the pairwise plant dynamics in local species-rich communities.

Investigations of microbial biogeography in extreme environments provide unique opportunities to disentangle the roles of environment and space in microbial community assembly. Here, we reported a comprehensive microbial biogeographic survey of 90 acid mine drainage (AMD) sediment samples from 18 mining sites of various mineral types across southern China. We found that environmental selection was strong in determining the AMD habitat species pool. However, microbial alpha diversity was primarily explained by mining sites rather than environmental factors, and microbial beta diversity correlated more strongly with geographic than environmental distance at both large and small spatial scales. Particularly, the presence/absence of widespread AMD habitat generalists was only correlated with geographic distance and independent of environmental variation. These distance–decay patterns suggested that spatial processes played a more important role in determining microbial compositional variation across space; which could be explained by the reinforced impacts of dispersal limitation in less fluid, spatially structured sediment habitat with diverse pre-existing communities. In summary, our findings suggested that the deterministic assembling and spatial constraints interact to shape microbial biogeography in AMD sediments; and provided implications that spatial processes should be considered when predicting microbial dynamics in response to severe environmental change across large spatial scales.

“Candidatus Bathyarchaeia” is a phylogenetically diverse and widely distributed lineage often in high abundance in anoxic submarine sediments; however, their evolution and ecological roles in terrestrial geothermal habitats are poorly understood. In the present study, 35 Ca. Bathyarchaeia metagenome-assembled genomes (MAGs) were recovered from hot spring sediments in Tibet and Yunnan, China. Phylogenetic analysis revealed all MAGs of Ca. Bathyarchaeia can be classified into 7 orders and 15 families. Among them, 4 families have been first discovered in the present study, significantly expanding the known diversity of Ca. Bathyarchaeia. Comparative genomics demonstrated Ca. Bathyarchaeia MAGs from thermal habitats to encode a large variety of genes related to carbohydrate degradation, which are likely a metabolic adaptation of these organisms to a lifestyle at high temperatures. At least two families are potential methanogens/alkanotrophs, indicating a potential for the catalysis of short-chain hydrocarbons. Three MAGs from Family-7.3 are identified as alkanotrophs due to the detection of an Mcr complex. Family-2 contains the largest number of genes relevant to alkyl-CoM transformation, indicating the potential for methylotrophic methanogenesis, although their evolutionary history suggests the ancestor of Ca. Bathyarchaeia was unable to metabolize alkanes. Subsequent lineages have acquired the ability via horizontal gene transfer. Overall, our study significantly expands our knowledge and understanding of the metabolic capabilities, habitat adaptations, and evolution of Ca. Bathyarchaeia in thermal environments. IMPORTANCECa. Bathyarchaeia MAGs from terrestrial hot spring habitats are poorly revealed, though they have been studied extensively in marine ecosystems. In this study, we uncovered the metabolic capabilities and ecological role of Ca. Bathyarchaeia in hot springs and give a comprehensive comparative analysis between thermal and nonthermal habitats to reveal the thermal adaptability of Ca. Bathyarchaeia. Also, we attempt to determine the evolutionary history of methane/alkane metabolism in Ca. Bathyarchaeia, since it appears to be the first archaea beyond Euryarchaeota which contains the mcrABG genes. The reclassification of Ca. Bathyarchaeia and significant genomic differences among different lineages largely expand our knowledge on these cosmopolitan archaea, which will be beneficial in guiding the future studies.

"Candidatus Aenigmarchaeota" ("Ca. Aenigmarchaeota") represents one of the earliest proposed evolutionary branches within the Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota (DPANN) superphylum. However, their ecological roles and potential host-symbiont interactions are still poorly understood. Here, eight metagenome-assembled genomes (MAGs) were reconstructed from hot spring ecosystems, and further in-depth comparative and evolutionary genomic analyses were conducted on these MAGs and other genomes downloaded from public databases. Although with limited metabolic capacities, we reported that "Ca. Aenigmarchaeota" in thermal environments harbor more genes related to carbohydrate metabolism than "Ca. Aenigmarchaeota" in nonthermal environments. Evolutionary analyses suggested that members from the Thaumarchaeota, Aigarchaeota, Crenarchaeota, and Korarchaeota (TACK) superphylum and Euryarchaeota contribute substantially to the niche expansion of "Ca. Aenigmarchaeota" via horizontal gene transfer (HGT), especially genes related to virus defense and stress responses. Based on co-occurrence network results and recent genetic exchanges among community members, we conjectured that "Ca. Aenigmarchaeota" may be symbionts associated with one MAG affiliated with the genus Pyrobaculum, though host specificity might be wide and variable across different "Ca. Aenigmarchaeota" organisms. This study provides significant insight into possible DPANN-host interactions and ecological roles of "Ca. Aenigmarchaeota." IMPORTANCE Recent advances in sequencing technology promoted the blowout discovery of super tiny microbes in the Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota (DPANN) superphylum. However, the unculturable properties of the majority of microbes impeded our investigation of their behavior and symbiotic lifestyle in the corresponding community. By integrating horizontal gene transfer (HGT) detection and co-occurrence network analysis on "Candidatus Aenigmarchaeota" ("Ca. Aenigmarchaeota"), we made one of the first attempts to infer their putative interaction partners and further decipher the potential functional and genetic interactions between the symbionts. We revealed that HGTs contributed by members from the Thaumarchaeota, Aigarchaeota, Crenarchaeota, and Korarchaeota (TACK) superphylum and Euryarchaeota conferred "Ca. Aenigmarchaeota" with the ability to survive under different environmental stresses, such as virus infection, high temperature, and oxidative stress. This study demonstrates that the interaction partners might be inferable by applying informatics analyses on metagenomic sequencing data.

Low soil phosphorus (P) bioavailability causes the widespread occurrence of P-limited terrestrial ecosystems around the globe. Exploring the factors influencing soil P bioavailability at large spatial scales is critical for managing these ecosystems. However, previous studies have mostly focused on abiotic factors. In this study, we explored the effects of microbial factors on soil P bioavailability of terrestrial ecosystems using a country-scale sampling effort. Our results showed that soil microbial biomass carbon (MBC) and acid phosphatase were important predictors of soil P bioavailability of agro- and natural ecosystems across China although they appeared less important than total soil P. The two microbial factors had a positive effect on soil P bioavailability of both ecosystem types and were able to mediate the effects of several abiotic factors (e.g., mean annual temperature). Meanwhile, we revealed that soil phytase could affect soil P bioavailability at the country scale via ways similar to those of soil MBC and acid phosphatase, a pattern being more pronounced in agroecosystems than in natural ecosystems. Moreover, we obtained evidence for the positive effects of microbial genes encoding these enzymes on soil P bioavailability at the country scale although their effect sizes varied between the two ecosystem types. Taken together, this study demonstrated the remarkable effects of microbial factors on soil P bioavailability at a large spatial scale, highlighting the importance to consider microbial factors in managing the widespread P-limited terrestrial ecosystems.

Fungi in acid mine drainage (AMD) environments are of great concern due to their potentials of decomposing organic carbon, absorbing heavy metals and reducing AMD acidity. Based on morphological analysis and ITS/18S high-throughput sequencing technology, previous studies have provided deep insights into the diversity and community composition of fungi in AMD environments. However, knowledge about physiology, metabolic potential and transcriptome profiles of fungi inhabiting AMD environments is still scarce. Here, we reported the physiological, genomic, and transcriptomic characterization of Acidiella bohemica SYSU C17045 to improve our understanding of the physiological, genomic, and transcriptomic mechanisms underlying fungal adaptation to AMD environments. A. bohemica was isolated from an AMD environment, which has been proved to be an acidophilic fungus in this study. The surface of A. bohemica cultured in AMD solutions was covered with a large number of minerals such as jarosite. We thus inferred that the A. bohemica might have the potential of biologically induced mineralization. Taking advantage of PacBio single-molecule real-time sequencing, we obtained the high-quality genome sequences of A. bohemica (50 Mbp). To our knowledge, this was the first attempt to employ a third-generation sequencing technology to explore the genomic traits of fungi isolated from AMD environments. Moreover, our transcriptomic analysis revealed that a series of genes in the A. bohemica genome were related to its metabolic pathways of C, N, S, and Fe as well as its adaptation mechanisms, including the response to acid stress and the resistance to heavy metals. Overall, our physiological, genomic, and transcriptomic data provide a foundation for understanding the metabolic potential and adaptation mechanisms of fungi in AMD environments.

Phosphate-solubilizing microbes (PSMs) drive the biogeochemical cycling of phosphorus (P) and hold promise for sustainable agriculture. However, their global distribution, overall diversity and application potential remain unknown. Here, we present the first synthesis of their biogeography, diversity and utility, employing data from 399 papers published between 1981 and 2017, the results of a nationwide field survey in China consisting of 367 soil samples, and a genetic analysis of 12986 genome-sequenced prokaryotic strains. We show that at continental to global scales, the population density of PSMs in environmental samples is correlated with total P rather than pH. Remarkably, positive relationships exist between the population density of soil PSMs and available P, nitrate-nitrogen and dissolved organic carbon in soil, reflecting functional couplings between PSMs and microbes driving biogeochemical cycles of nitrogen and carbon. More than 2704 strains affiliated with at least nine archaeal, 88 fungal and 336 bacterial species were reported as PSMs. Only 2.59% of these strains have been tested for their efficiencies in improving crop growth or yield under field conditions, providing evidence that PSMs are more likely to exert positive effects on wheat growing in alkaline P-deficient soils. Our systematic genetic analysis reveals five promising PSM genera deserving much more attention.

The haloalkalitolerant bacterium Egicoccus halophilus EGI 80432