Enhanced monsoon duration and soil acidification from acid rain are expected to impact the distribution of microbial communities in surface and subsurface environments, although these impacts are poorly understood for most systems. In central China, soluble carbonate bedrock forms extensive karst landscapes. Current predictions are that the amount of monsoonal precipitation and acid rainfall in central China will increase, which is expected to lead to changes in the pH balance of karst ecosystems. To evaluate the role of pH, total organic carbon, and other geochemical parameters (e.g., Ca2+, Mg2+, NH4+, NOx, SO42-) in shaping bacterial communities within a single karst system in central China, samples were collected from the thin surface soils overlying Heshang Cave, cave sediments, and weathered cave passage rocks from the entrance, twilight, and dark zones, as well as from epikarstic drip waters inside the cave. Illumina sequencing of 16S rRNA genes and multivariate statistical analyses revealed that each tested community was distinct and the community variability was significantly correlated with pH, total organic carbon, and potassium concentrations. Specifically, surface soils were dominated by Acidobacteria, Verrucomicrobia and Planctomycetes, and diversity significantly decreased with acidic pH values. Nitrospirae, Gemmatimonadetes, Firmicutes, and Chloroflexi were unique to cave sediments, while Actinobacteria and Proteobacteria dominated weathered rocks and drip waters, respectively. The results reveal important implications regarding the effects of acidification on bacterial communities in karst areas, and on the control of pH in shaping bacterial communities throughout a karst system. Increased water flux into and through karst habitats due to monsoonal precipitation may result in deeper penetration of acidic solutions into karst and shift the bacterial communities inside the cave in the future.

High abundance and widespread distribution of the archaeal phylum Bathyarchaeota in marine environment have been recognized recently, but knowledge about Bathyarchaeota in terrestrial settings and their correlation with environmental parameters is fairly limited. Here we reported the abundance of Bathyarchaeota members across different ecosystems and their correlation with environmental factors by constructing 16S rRNA clone libraries of peat from the Dajiuhu Peatland, coupling with bioinformatics analysis of 16S rRNA data available to date in NCBI database. In total, 1456 Bathyarchaeota sequences from 28 sites were subjected to UniFrac analysis based on phylogenetic distance and multivariate regression tree analysis of taxonomy. Both phylogenetic and taxon-based approaches showed that salinity, total organic carbon and temperature significantly influenced the distribution of Bathyarchaeota across different terrestrial habitats. By applying the ecological concept of ‘indicator species’, we identify 9 indicator groups among the 6 habitats with the most in the estuary sediments. Network analysis showed that members of Bathyarchaeota formed the “backbone” of archaeal community and often co-occurred with Methanomicrobia. These results suggest that Bathyarchaeota may play an important ecological role within archaeal communities via a potential symbiotic association with Methanomicrobia. Our results shed light on understanding of the biogeography, potential functions of Bathyarchaeota and environment conditions that influence Bathyarchaea distribution in terrestrial settings.

Caves are nutrient-limited and dark subterranean ecosystems. To date, attention has been focused on geological research of caves in China, whilst indigenous microbial diversity has been insufficiently characterized. Here, we report the fungal diversity in the pristine, oligotrophic, karst Heshang Cave, central China, using a culture-dependent method coupled with the analysis of the fungal rRNA-ITS gene sequences. A total of 194 isolates were obtained with six different media from 14 sampling sites of sediments, weathered rocks, and bat guanos. Phylogenetic analysis clustered the 194 sequenced isolates into 33 genera within 15 orders of three phyla, Ascomycota, Basidiomycota, and Zygomycota, indicating a high degree of fungal diversity in the Heshang Cave. Notably, 16 out of the 36 fungal genera were also frequently observed in solution caves around the world and 23 genera were previously found in carbonate cave, indicating potential similarities among fungal communities in cave ecosystems. However, 10 genera in this study were not reported previously in any solution caves, thus expanding our knowledge about fungal diversity in cave ecosystems. Moreover, culturable fungal diversity varied from one habitat to another within the cave, being the highest in sediments, followed by weathered rocks and bat guanos as indicated by α-diversity indexes. At the genus level, Penicillium accounted for 40, 54, and 52% in three habitats of sediments, weathered rocks, and bat guanos, respectively. Trichoderma, Paecilomyces, and Aspergillus accounted for 9, 22, and 37% in the above habitats, correspondingly. Despite of the dominance of Penicillium in all samples, β-diversity index indicated significant differences between each two fungal communities in the three habitats in view of both the composition and abundance. Our study is the first report on fungal communities in a natural pristine solution cave system in central China and sheds light on fungal diversity and functions in cave ecosystems.

Deciphering of the mycobiome in pristine karst caves has been impeded by constraints of remote locations, inaccessibility to specimens and technical limitations, which greatly restricted in-depth understanding of mycobiomes in subterranean ecosystem. Here, mycobiomes of Heshang Cave in south-western karst region of China were investigated by Illumina HiSeq sequencing of fungal rRNA-ITS1 gene across different habitats. In total 793,502 ITS1 reads and 2,179 OTUs from 778 Mb reads after stringent quality control (Q30) and 453 genera, 72 orders and 19 classes within 6 phyla were detected. Ascomycota (42% OTUs) dominated across the five habitats. Shannon-Wiener index varied from 1.25 to 7.62 and community richness was highest in drip waters, followed by weathered rocks, bat guanos, sediments, and air samples. Mycobiomes displayed specificity to five habitats and more distinct OTUs were found in weathered rocks (12%) and drip waters (9%). In contrast, only 6.60% core OTUs were shared by five habitats. Notably, weathered rocks possessed more indicator groups and were revealed for the first time to be dominated by Sordariomycetes (43%). The community richness of air mycobiomes increased from cave entrance to the innermost part and dominated by the indicator groups of Penicillium mallochii (>30%) and P. herquei (>9%). Our work represents the largest attempt to date to a systematical investigation of oligotrophic solution-cave-associated mycobiomes in China. Our discovery of high diversity of mycobiomes in Heshang Cave also suggests that eukaryotic microorganisms may play a crucial role in subsurface environments.