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.
Keywords: Heshang Cave; ITS1 region amplicons; Illumina HiSeq sequencing; geomicrobiology; mycobiomes; subterranean karst ecosystem.