A potential central role of Thaumarchaeota in N-Cycling in a semi-arid environment, Fort Stanton Cave, Snowy River passage, New Mexico, USA

The geomicrobiology of limestone, sulfuric acid speleogenetic, and volcanic caves: basic concepts and future perspectives

Caves are ubiquitous subterranean voids, accounting for a still largely unexplored surface of the Earth underground. Due to the absence of sunlight and physical segregation, caves are naturally colonized by microorganisms that have developed distinctive capabilities to thrive under extreme conditions of darkness and oligotrophy. Here, the microbiomes colonizing three frequently studied cave types, i.e., limestone, sulfuric acid speleogenetic (SAS), and lava tubes among volcanic caves, have comparatively been reviewed. Geological configurations, nutrient availability, and energy flows in caves are key ecological drivers shaping cave microbiomes through photic, twilight, transient, and deep cave zones. Chemoheterotrophic microbial communities, whose sustenance depends on nutrients supplied from outside, are prevalent in limestone and volcanic caves, while elevated inorganic chemical energy is available in SAS caves, enabling primary production through chemolithoautotrophy. The 16S rRNA-based metataxonomic profiles of cave microbiomes were retrieved from previous studies employing the Illumina platform for sequencing the prokaryotic V3-V4 hypervariable region to compare the microbial community structures from different cave systems and environmental samples. Limestone caves and lava tubes are colonized by largely overlapping bacterial phyla, with the prevalence of Pseudomonadota and Actinomycetota, whereas the co-dominance of Pseudomonadota and Campylobacterota members characterizes SAS caves. Most of the metataxonomic profiling data have so far been collected from the twilight and transient zones, while deep cave zones remain elusive, deserving further exploration. Integrative approaches for future geomicrobiology studies are suggested to gain comprehensive insights into the different cave types and zones. This review also poses novel research questions for unveiling the metabolic and genomic capabilities of cave microorganisms, paving the way for their potential biotechnological applications.

Comparative genomics reveals the adaptation of ammonia-oxidising Thaumarchaeota to arid soils

Thaumarchaeota are predominant in oligotrophic habitats such as deserts and arid soils, but their adaptations to these arid conditions are not well understood. In this study, we assembled 23 Thaumarchaeota genomes from arid and semi-arid soils collected from the Inner Mongolia Steppe and the Qinghai-Tibet Plateau. Using a comparative genomics approach, integrated with 614 Thaumarchaeota genomes from public databases, we identified the traits and evolutionary forces that contribute to their adaptations to aridity. Our results showed that the newly assembled genomes represent an early diverging group within the lineage of ammonia-oxidising Thaumarchaeota. While the genomic functions previously identified in arid soil lineages were conserved across terrestrial, shallow-ocean and deep-ocean lineages, several traits likely contribute to Thaumarchaeota's adaptation to aridity. These include chlorite dismutase, arsenate reductase, V-type ATPase and genes dealing with oxidative stresses. The acquisition and loss of traits at the last common ancestor of arid soil lineages may have facilitated the specialisation of Thaumarchaeota in arid soils. Additionally, the acquisition of unique adaptive traits, such as a urea transporter, Ca2+ :H+ antiporter, mannosyl-3-phosphoglycerate synthase and phosphatase, DNA end-binding protein Ku and phage shock protein A, further distinguishes arid soil Thaumarchaeota. This study provides evidence for the adaptations of Thaumarchaeota to arid soil, enhancing our understanding of the nitrogen and carbon cycling driven by Thaumarchaeota in drylands.

Metagenome Analysis of Speleothem Microbiome from Subterranean Cave Reveals Insight into Community Structure, Metabolic Potential, and BGCs Diversity

The Borra caves, the second largest subterranean karst cave ecosystem in the Indian sub-continent, are located at the Ananthagiri hills of Araku Valley in the Alluri district of Andhra Pradesh, India. The present investigation applied a shotgun metagenomic approach to gain insights into the microbial community structure, metabolic potential, and biosynthetic gene cluster (BGC) diversity of the microbes colonizing the surface of the speleothems from the aphotic zone of Borra caves. The taxonomic analysis of the metagenome data illustrated that the speleothem-colonizing core microbial community was dominated mainly by Alpha-, Beta-, and Gamma-Proteobacteria, Actinobacteria, Firmicutes, and Bacteroidetes. The key energy metabolic pathways analysis provides strong evidence of chemolithoautotrophic and chemoheterotrophic modes of nutrition in the speleothem-colonizing microbial community. Metagenome data suggests that sulfur reducers and sulfur-disproportionating microbes might play a vital role in energy generation in this ecosystem. Our metagenome data also suggest that the dissimilatory nitrifiers and nitrifying denitrifiers might play an essential role in conserving nitrogen pools in the ecosystem. Furthermore, metagenome-wide BGCs mining retrieved 451 putative BGCs; NRPS was the most abundant (24%). Phylogenetic analysis of the C domain of NRPS showed that sequences were distributed across all six function categories of the known C domain, including several novel subclades. For example, a novel subclade had been recovered within the LCL domain clade as a sister subclade of immunosuppressant cyclosporin encoding C domain sequences. Our result suggested that subterranean cave microbiomes might be a potential reservoir of novel microbial metabolites.

Nitrate determines the bacterial habitat specialization and impacts microbial functions in a subsurface karst cave

Karst caves are usually considered as natural laboratories to study pristine microbiomes in subsurface biosphere. However, effects of the increasingly detected nitrate in underground karst ecosystem due to the acid rain impact on microbiota and their functions in subsurface karst caves have remained largely unknown. In this study, samples of weathered rocks and sediments were collected from the Chang Cave, Hubei province and subjected to high-throughput sequencing of 16S rRNA genes. The results showed that nitrate significantly impacted bacterial compositions, interactions, and functions in different habitats. Bacterial communities clustered according to their habitats with distinguished indicator groups identified for each individual habitat. Nitrate shaped the overall bacterial communities across two habitats with a contribution of 27.2%, whereas the pH and TOC, respectively, structured bacterial communities in weathered rocks and sediments. Alpha and beta diversities of bacterial communities increased with nitrate concentration in both habitats, with nitrate directly affecting alpha diversity in sediments, but indirectly on weathered rocks by lowering pH. Nitrate impacted more on bacterial communities in weathered rocks at the genus level than in sediments because more genera significantly correlated with nitrate concentration in weathered rocks. Diverse keystone taxa involved in nitrogen cycling were identified in the co-occurrence networks such as nitrate reducers, ammonium-oxidizers, and N₂-fixers. Tax4Fun2 analysis further confirmed the dominance of genes involved in nitrogen cycling. Genes of methane metabolism and carbon fixation were also dominant. The dominance of dissimilatory and assimilatory nitrate reduction in nitrogen cycling substantiated nitrate impact on bacterial functions. Our results for the first time revealed the impact of nitrate on subsurface karst ecosystem in terms of bacterial compositions, interactions, and functions, providing an important reference for further deciphering the disturbance of human activities on the subsurface biosphere.

Archaea and their interactions with bacteria in a karst ecosystem

Karst ecosystems are widely distributed around the world, accounting for 15-20% of the global land area. However, knowledge on microbial ecology of these systems does not match with their global importance. To close this knowledge gap, we sampled three niches including weathered rock, sediment, and drip water inside the Heshang Cave and three types of soils overlying the cave (forest soil, farmland soil, and pristine karst soil). All these samples were subjected to high-throughput sequencing of V4-V5 region of 16S rRNA gene and analyzed with multivariate statistical analysis. Overall, archaeal communities were dominated by Thaumarchaeota, whereas Actinobacteria dominated bacterial communities. Thermoplasmata, Nitrosopumilaceae, Aenigmarchaeales, Crossiella, Acidothermus, and Solirubrobacter were the important predictor groups inside the Heshang Cave, which were correlated to NH₄ ⁺ availability. In contrast, Candidatus Nitrososphaera, Candidatus Nitrocosmicus, Thaumarchaeota Group 1.1c, and Pseudonocardiaceae were the predictors outside the cave, whose distribution was correlated with pH, Ca²⁺, and NO₂ ^-. Tighter network structures were found in archaeal communities than those of bacteria, whereas the topological properties of bacterial networks were more similar to those of total prokaryotic networks. Both chemolithoautotrophic archaea (Candidatus Methanoperedens and Nitrosopumilaceae) and bacteria (subgroup 7 of Acidobacteria and Rokubacteriales) were the dominant keystone taxa within the co-occurrence networks, potentially playing fundamental roles in obtaining energy under oligotrophic conditions and thus maintaining the stability of the cave ecosystem. To be noted, all the keystone taxa of karst ecosystems were related to nitrogen cycling, which needs further investigation, particularly the role of archaea. The predicted ecological functions in karst soils mainly related to carbohydrate metabolism, biotin metabolism, and synthesis of fatty acid. Our results offer new insights into archaeal ecology, their potential functions, and archaeal interactions with bacteria, which enhance our understanding about the microbial dark matter in the subsurface karst ecosystems.

Convergent Community Assembly among Globally Separated Acidic Cave Biofilms

Acidophilic bacteria and archaea inhabit extreme geochemical "islands" that can tell us when and how geographic barriers affect the biogeography of microorganisms. Here, we describe microbial communities from extremely acidic (pH 0 to 1) biofilms, known as snottites, from hydrogen sulfide-rich caves. Given the extreme acidity and subsurface location of these biofilms, and in light of earlier work showing strong geographic patterns among snottite Acidithiobacillus populations, we investigated their structure and diversity in order to understand how geography might impact community assembly. We used 16S rRNA gene cloning and fluorescence in situ hybridization (FISH) to investigate 26 snottite samples from four sulfidic caves in Italy and Mexico. All samples had very low biodiversity and were dominated by sulfur-oxidizing bacteria in the genus Acidithiobacillus. Ferroplasma and other archaea in the Thermoplasmatales ranged from 0 to 50% of total cells, and relatives of the bacterial genera Acidimicrobium and Ferrimicrobium were up to 15% of total cells. Rare phylotypes included Sulfobacillus spp. and members of the phyla "Candidatus Dependentiae" and "Candidatus Saccharibacteria" (formerly TM6 and TM7). Although the same genera of acidophiles occurred in snottites on separate continents, most members of those genera represent substantially divergent populations, with 16S rRNA genes that are only 95 to 98% similar. Our findings are consistent with a model of community assembly where sulfidic caves are stochastically colonized by microorganisms from local sources, which are strongly filtered through environmental selection for extreme acid tolerance, and these different colonization histories are maintained by dispersal restrictions within and among caves. IMPORTANCE Microorganisms that are adapted to extremely acidic conditions, known as extreme acidophiles, are catalysts for rock weathering, metal cycling, and mineral formation in naturally acidic environments. They are also important drivers of large-scale industrial processes such as biomining and contaminant remediation. Understanding the factors that govern their ecology and distribution can help us better predict and utilize their activities in natural and engineered systems. However, extremely acidic habitats are unusual in that they are almost always isolated within circumneutral landscapes. So where did their acid-adapted inhabitants come from, and how do new colonists arrive and become established? In this study, we took advantage of a unique natural experiment in Earth's subsurface to show how isolation may have played a role in the colonization history, community assembly, and diversity of highly acidic microbial biofilms.

Impact of environmental factors on the ammonia-oxidizing and denitrifying microbial community and functional genes along soil profiles from different ecologically degraded areas in the Siding mine

Ammonia oxidizers (ammonia-oxidizing bacteria (AOB amoA) and ammonia-oxidizing archaea (AOA amoA)) and denitrifiers (encoded by nirS, nirK and nosZ) in the soil nitrogen cycle exist in a variety of natural ecosystems. However, little is known about the contribution of these five N-related functional genes to nitrification and denitrification in the soil profile in severely ecologically degraded areas. Therefore, in the present study, the abundance, diversity and community composition of AOA, AOB, nirS, nirK and nosZ were investigated in the soil profiles of different ecologically degraded areas in the Siding mine. The results indicated that, at the phylum level, the dominant archaea were Crenarchaeota and Thaumarchaeota and the dominant bacteria were Proteobacteria. Heavy metal contents had a great impact on AOA amoA, nirS and nirK gene abundances. AOA amoA contributed more during the ammonia oxidation process and was better adapted for survival in heavy metal-contaminated environments. In addition to heavy metals, the soil organic matter (SOM) content and C/N ratio had strong effects on the AOA and AOB community diversity and structure. In addition, variations in the net ammonification and nitrification rates were proportional to AOA amoA abundance along the soil profile. The soil C/N ratio, soil available phosphorus content and soil moisture influenced the denitrification process. Both soil available phosphorus and moisture were more strongly related to nosZ than to nirS and nirK. In addition, nosZ presented a higher correlation with the nosZ/(nirS + nirK) ratio. Moreover, nosZ/(nirS + nirK) was the key functional gene group that drove the major processes for NH₄⁺-N and NO₃^--N transformation. This study demonstrated the role and importance of soil property impacts on N-related microbes in the soil profile and provided a better understanding of the role and importance of N-related functional genes and their contribution to soil nitrification and denitrification processes in highly degraded areas in the Siding mine.

Major impacts of caving activities on cave microbial diversity: case study of Morca Cave, Turkey

Some of microorganisms identified in cave ecosystems have been reported to play a permanent and strategic role for maintaining life of these environments. Human entrance into caves can induce some changes of cave physicochemical parameters which ultimately impacts the living organisms. In these facts, for the first time, Morca Cave was explored in a purpose to evaluate the impacts that can be caused by the human activities on microbial diversity in a limited period of time. Morca is a karts cave located in the Taşeli Plateau in the Middle Taurus mountains in Turkey. The first entrance into this cave was in 2018, and expeditors reached to - 500 m. During the second expedition in 2019, a camp was established at the - 1040-m depth during 4 days. In order to evaluate the human impacts in a new explored cave, this camping depth is chosen to be our studied area because it was its first entrance. Before the installation and at the end of the camp, sediments and surface samples were taken from different points of the camp area and around. Sequencing of 16 s rRNA of each sample to isolate DNA by using the next-generation sequencing (NGS) method was performed. The profile of the microbial diversity before the camping revealed that the class Thermoplasmata was dominated the archaea group and Gammaproteobacteria and Alphaproteobacteria were the most dominant bacterial classes. After the camp, most studied sites were noted with a decrease of microbial diversity especially the previous cited classes strains. Bacteria belonging to Bacilli class have increased after the camp. Increase of bacteria that are belonging to Bacteroidia has also observed in the most active areas. This present study highlight how cave microbial diversity can respond to the human activities within a short period inside a closed cave. Furthermore, it may constitute a solid basis and support on the improvement of techniques for cave management and expedition planning for the conservation of cave nature.

Microbial roles in cave biogeochemical cycling

Among fundamental research questions in subterranean biology, the role of subterranean microbiomes playing in key elements cycling is a top-priority one. Karst caves are widely distributed subsurface ecosystems, and cave microbes get more and more attention as they could drive cave evolution and biogeochemical cycling. Research have demonstrated the existence of diverse microbes and their participance in biogeochemical cycling of elements in cave environments. However, there are still gaps in how these microbes sustain in caves with limited nutrients and interact with cave environment. Cultivation of novel cave bacteria with certain functions is still a challenging assignment. This review summarized the role of microbes in cave evolution and mineral deposition, and intended to inspire further exploration of microbial performances on C/N/S biogeocycles.

Islands Within Islands: Bacterial Phylogenetic Structure and Consortia in Hawaiian Lava Caves and Fumaroles

Lava caves, tubes, and fumaroles in Hawai‘i present a range of volcanic, oligotrophic environments from different lava flows and host unexpectedly high levels of bacterial diversity. These features provide an opportunity to study the ecological drivers that structure bacterial community diversity and assemblies in volcanic ecosystems and compare the older, more stable environments of lava tubes, to the more variable and extreme conditions of younger, geothermally active caves and fumaroles. Using 16S rRNA amplicon-based sequencing methods, we investigated the phylogenetic distinctness and diversity and identified microbial interactions and consortia through co-occurrence networks in 70 samples from lava tubes, geothermal lava caves, and fumaroles on the island of Hawai‘i. Our data illustrate that lava caves and geothermal sites harbor unique microbial communities, with very little overlap between caves or sites. We also found that older lava tubes (500–800 yrs old) hosted greater phylogenetic diversity (Faith's PD) than sites that were either geothermally active or younger (<400 yrs old). Geothermally active sites had a greater number of interactions and complexity than lava tubes. Average phylogenetic distinctness, a measure of the phylogenetic relatedness of a community, was higher than would be expected if communities were structured at random. This suggests that bacterial communities of Hawaiian volcanic environments are phylogenetically over-dispersed and that competitive exclusion is the main driver in structuring these communities. This was supported by network analyses that found that taxa (Class level) co-occurred with more distantly related organisms than close relatives, particularly in geothermal sites. Network “hubs” (taxa of potentially higher ecological importance) were not the most abundant taxa in either geothermal sites or lava tubes and were identified as unknown families or genera of the phyla, Chloroflexi and Acidobacteria. These results highlight the need for further study on the ecological role of microbes in caves through targeted culturing methods, metagenomics, and long-read sequence technologies.

Comparative microbial antibiotic resistome between urban and deep forest environments

A paradoxical result of using antibiotics to eradicate microbial pathogens is the emergence of a vast number of resistant microbes in various environments. The concern that environmental microbes will inevitably become resistant to virtually every clinically usable antibiotics has been exacerbated by the spread of these resistance genes across different environments and the emergence of multidrug resistant phenotypes. Here, we provide metagenomic insights into the microbiomes and resistomes of 16 soil samples collected from hospitals, residential areas, and forest parks in the megacity of Beijing and deep forests in the Yunnan province. Using Illumina HiSeq sequencing, we investigated the microbial diversity within the metagenomic shotgun reads and identified 486 antibiotic-resistant genes (ARGs) classified into 30 types from these samples, among which multidrug resistance genes were the most abundant. Our results present an important reference and direct comparison of microbial antibiotic resistomes of soil samples from a megacity and deep forests and extend our understanding of the spread of ARGs in modern urban and natural environments.

Bacteria and Metabolic Potential in Karst Caves Revealed by Intensive Bacterial Cultivation and Genome Assembly

Karst caves are widely distributed subsurface systems, and the microbiomes therein are proposed to be the driving force for cave evolution and biogeochemical cycling. In past years, culture-independent studies on the microbiomes of cave systems have been conducted, yet intensive microbial cultivation is still needed to validate the sequence-derived hypothesis and to disclose the microbial functions in cave ecosystems. In this study, the microbiomes of two karst caves in Guizhou Province in southwest China were examined. A total of 3,562 bacterial strains were cultivated from rock, water, and sediment samples, and 329 species (including 14 newly described species) of 102 genera were found. We created a cave bacterial genome collection of 218 bacterial genomes from a karst cave microbiome through the extraction of 204 database-derived genomes and de novo sequencing of 14 new bacterial genomes. The cultivated genome collection obtained in this study and the metagenome data from previous studies were used to investigate the bacterial metabolism and potential involvement in the carbon, nitrogen, and sulfur biogeochemical cycles in the cave ecosystem. New N₂-fixing Azospirillum and alkane-oxidizing Oleomonas species were documented in the karst cave microbiome. Two pcaIJ clusters of the β-ketoadipate pathway that were abundant in both the cultivated microbiomes and the metagenomic data were identified, and their representatives from the cultivated bacterial genomes were functionally demonstrated. This large-scale cultivation of a cave microbiome represents the most intensive collection of cave bacterial resources to date and provides valuable information and diverse microbial resources for future cave biogeochemical research.IMPORTANCE Karst caves are oligotrophic environments that are dark and humid and have a relatively stable annual temperature. The diversity of bacteria and their metabolisms are crucial for understanding the biogeochemical cycling in cave ecosystems. We integrated large-scale bacterial cultivation with metagenomic data mining to explore the compositions and metabolisms of the microbiomes in two karst cave systems. Our results reveal the presence of a highly diversified cave bacterial community, and 14 new bacterial species were described and their genomes sequenced. In this study, we obtained the most intensive collection of cultivated microbial resources from karst caves to date and predicted the various important routes for the biogeochemical cycling of elements in cave ecosystems.

Microbiomes of different ages in Rendzic Leptosols in the Crimean Peninsula

Rendzic Leptosols are intrazonal soils formed on limestone bedrock. The specialty of these soils is that parent rock material is more influential in shaping soil characteristics than zonal factors such as climate, especially during soil formation. Unlike fast evolving Podzols due to their leaching regime, Leptosols do not undergo rapid development due to the nature of the limestone. Little is known how microbiome reflects this process, so we assessed microbiome composition of Rendzic Leptosols of different ages, arising from disruption and subsequent reclamation. The mountains and foothills that cover much of the Crimean Peninsula are ideal for this type of study, as the soils were formed on limestone and have been subjected to anthropogenic impacts through much of human history. Microbiomes of four soil sites forming a chronosequence, including different soil horizons, were studied using sequencing of 16S rRNA gene libraries and quantitative PCR. Dominant phyla for all soil sites were Actinobacteria, Proteobacteria, Acidobacteria, Bacteroidetes, Thaumarchaeota, Planctomycetes, Verrucomicrobia and Firmicutes. Alpha diversity was similar across sites and tended to be higher in topsoil. Beta diversity showed that microbiomes diverged according to the soil site and the soil horizon. The oldest and the youngest soils had the most similar microbiomes, which could have been caused by their geographic proximity. Oligotrophic bacteria from Chitinophagaceae, Blastocatellaceae and Rubrobacteriaceae dominated the microbiome of these soils. The microbiome of 700-year old soil was the most diverse. This soil was from the only study location with topsoil formed by plant litter, which provided additional nutrients and could have been the driving force of this differentiation. Consistent with this assumption, this soil was abundant in copiotrophic bacteria from Proteobacteria and Actinobacteria phyla. The microbiome of 50-year old Leptosol was more similar to the microbiome of benchmark soil than the microbiome of 700-year old soil, especially by weighted metrics. CCA analysis, in combination with PERMANOVA, linked differences in microbiomes to the joint change of all soil chemical parameters between soil horizons. Local factors, such as parent material and plant litter, more strongly influenced the microbiome composition in Rendzic Leptosols than soil age.

Mapping metabolic activity at single cell resolution in intact volcanic fumarole sediment

Interactions among microorganisms and their mineralogical substrates govern the structure, function and emergent properties of microbial communities. These interactions are predicated on spatial relationships, which dictate metabolite exchange and access to key substrates. To quantitatively assess links between spatial relationships and metabolic activity, this study presents a novel approach to map all organisms, the metabolically active subset and associated mineral grains, all while maintaining spatial integrity of an environmental microbiome. We applied this method at an outgassing fumarole of Vanuatu's Marum Crater, one of the largest point sources of several environmentally relevant gaseous compounds, including H2O, CO2 and SO2. With increasing distance from the sediment-air surface and from mineral grain outer boundaries, organism abundance decreased but the proportion of metabolically active organisms often increased. These protected niches may provide more stable conditions that promote consistent metabolic activity of a streamlined community. Conversely, exterior surfaces accumulate more organisms that may cover a wider range of preferred conditions, implying that only a subset of the community will be active under any particular environmental regime. More broadly, the approach presented here allows investigators to see microbial communities ‘as they really are’ and explore determinants of metabolic activity across a range of microbiomes.