Distribution of rhizosphere fungi of Kobresia humilis on the Qinghai-Tibet Plateau

Kobresia humilis is a major species in the alpine meadow communities of the Qinghai-Tibet Plateau (QTP); it plays a crucial role in maintaining the ecological balance of these meadows. Nevertheless, little is known about the rhizosphere fungi associated with K. humilis on the Qinghai Tibet Plateau. In this study, we used Illumina Miseq to investigate the fungal diversity, community structure, and ecological types in the root and rhizosphere soil of K. humilis across eight areas on the QTP and analyzed the correlation between rhizosphere fungi of K. humilis and environmental factors. A total of 19,423 and 25,101 operational taxonomic units (OTUs) were obtained from the roots and rhizosphere soil of K. humilis. These were classified into seven phyla, 25 classes, 68 orders, 138 families, and 316 genera in the roots, and nine phyla, 31 classes, 76 orders, 152 families, and 407 genera in the rhizosphere soil. There were 435 and 415 core OTUs identified in root and rhizosphere soil, respectively, which were categorized into 68 and 59 genera, respectively, with 25 shared genera. Among them, the genera with a relative abundance >1% included Mortierella, Microscypha, Floccularia, Cistella, Gibberella, and Pilidium. Compared with the rhizosphere soil, the roots showed five differing fungal community characteristics, as well as differences in ecological type, and in the main influencing environmental factors. First, the diversity, abundance, and total number of OTUs in the rhizosphere soil of K. humilis were higher than for the endophytic fungi in the roots by 11.85%, 9.85%, and 22.62%, respectively. The composition and diversity of fungal communities also differed between the eight areas. Second, although saprotroph-symbiotrophs were the main ecological types in both roots and rhizosphere soil; there were 62.62% fewer pathotrophs in roots compared to the rhizosphere soil. Thirdly, at the higher altitude sites (3,900–4,410 m), the proportion of pathotroph fungi in K. humilis was found to be lower than at the lower altitude sites (3,200–3,690 m). Fourthly, metacommunity-scale network analysis showed that during the long-term evolutionary process, ZK (EICZK = 1) and HY (EICHY = 1) were critical sites for development of the fungal community structure in the roots and rhizosphere soil of K. humilis, respectively. Fifthly, canonical correspondence analysis (CCA) showed that key driving factors in relation to the fungal community were longitude (R2 = 0.5410) for the root community and pH (R2 = 0.5226) for the rhizosphere soil community. In summary, these results show that K. humilis fungal communities are significantly different in the root and rhizosphere soil and at the eight areas investigated, indicating that roots select for specific microorganisms in the soil. This is the first time that the fungal distribution of K. humilis on the QTP in relation to long-term evolutionary processes has been investigated. These findings are critical for determining the effects of environmental variables on K. humilis fungal communities and could be valuable when developing guidance for ecological restoration and sustainable utilization of the biological resources of the QTP.