Thermokarst lakes are a ubiquitous and important landscape feature in cold regions and are changing tremendously due to the accelerated climate change. In thermokarst lakes, sediment and water are two distinct but highly interconnected habitats, harboring different bacterial communities in terms of taxonomic composition. However, the co-occurrence networks of these bacterial communities remain unclear. Here, we investigate the co-occurrence ecological networks of sediment and water bacterial communities for thermokarst lakes in the Yellow River Source Area on the Qinghai-Tibet Plateau. The results show that the bacterial communities construct distinct co-occurrence networks in sediment and water. The metacommunity network was parsed into four major modules formed by the operational taxonomic units (OTUs) enriched in sediment or water independently, and water-enriched OTUs exhibited much closer interconnections than sediment-enriched OTUs. When considering the sediment and water bacterial networks separately, different topological properties and modular patterns present: the sediment bacterial network was more clustered while the modules less responded to the environmental variables. On the contrary, the water bacterial network was more complex with the OTUs more interconnected and its modules more responded to the environmental variables. Moreover, the results of the structural equation model suggest that, by the influence of environmental variations on individual modules, the water bacterial communities would be more vulnerable under the fact of accelerating climate change. This study provides insights beyond a conventional taxonomic perspective, adding our knowledge of the potential mechanisms structuring bacterial community assembly and improving our prediction of the responses of this fast-changing ecosystem to future climate change.

Thermokarst lakes are forming from permafrost thaw and are severely affected by accelerating climate change. Sediment and water in these lakes are distinct habitats but closely connected. However, our understanding of the differences and linkages between sediment and water in thermokarst lakes remains largely unknown, especially from the perspective of community assembly mechanisms. Here, we examined bacterial communities in sediment and water in thermokarst lakes in the Yellow River Source area, Qinghai-Tibet Plateau. Bacterial taxa were divided into abundant and rare according to their relative abundance, and the Sorensen dissimilarity (β_sor) was partitioned into turnover (β_turn) and nestedness (β_nest). The whole bacterial communities and the abundant and rare subcommunities differed substantially between sediment and water in taxonomical composition, α-diversity, and β-diversity. Sediment had significantly lower α-diversity indexes but higher β-diversity than water. In general, bacterial communities are predominantly governed by strong turnover processes (β_turn/β_sor ratio of 0.925). Bacterial communities in sediment had a significantly higher β_turn/β_sor ratio than in water. Abundant subcommunities were significantly lower in the β_turn/β_sor ratio compared with rare subcommunities. The results suggest that the bacterial communities of thermokarst lakes, especially rare subcommunities or particularly in sediment, might be strongly structured by heterogeneity in the source material, environmental filtering, and geographical isolation, leading to compositionally distinct communities. This integral study increased our current knowledge of thermokarst lakes, enhancing our understanding of the community assembly rules and ecosystem structures and processes of these rapidly changing and vulnerable ecosystems.

Thermokarst lakes are widely distributed in cold regions as a result of ice-rich permafrost thaw. Disentangling the biogeography of abundant and rare microbes is essential to understanding the environmental influences, assembly mechanisms, and responses to climate change of bacterial communities in thermokarst lakes. In light of this, we assessed the abundant and rare bacterial subcommunities in sediments from thermokarst lakes across the Qinghai-Tibet Plateau (QTP). The operational taxonomic unit (OTU) richness was more strongly associated with location and climate factors for abundant subcommunities, while more strongly associated with physicochemical variables for rare subcommunities. The relative abundance of abundant and rare taxa showed opposite patterns with abundant taxa having greater relative abundance at higher latitude and pH, but at lower mean annual precipitation and nutrients. Both the abundant and rare subcommunities had a clear distribution pattern along the gradient of latitude and mean annual precipitation. Abundant subcommunities were dominantly shaped by dispersal limitation processes (80.9%), while rare subcommunities were shaped almost equally by deterministic (47.3%) and stochastic (52.7%) processes. The balance between stochastic and deterministic processes was strongly environmentally adjusted for rare subcommunities, while not associated with environmental changes for abundant subcommunities. The results shed light on biogeography patterns and structuring mechanisms of bacterial communities in thermokarst lakes, improving our ability to predict the influences of future climate change on these lakes.