Lake microbiome composition determines community adaptability to warming perturbations

Background Microbial communities are crucial for ecosystems. A central goal of microbial ecological research is to simplify the biodiversity in natural environments and quantify the functional roles to explore how the ecosystems respond under different environmental conditions. However, the research on the stability and dynamics of lake microbes in response to repeated warming stress is limited. Methods To exclude confounding environmental factors, we conducted a 20-day repeated warming simulation experiment to examine the composition and function dynamics of lake microbial communities. Results Experimental warming significantly altered the community structure of bacteria instead of fungi. Microbial community structure, together with microbial biomass, jointly regulated the function of microbial communities. The plummeting of aerobic denitrifiers Pseudomonadaceae decreased by 99% ( P < 0.001) after high temperature, leading to reduced microbial nitrogen metabolism on nitrogen respiration and nitrate respiration. Under warming conditions, the microbial community with higher adaptability showed more positive correlations and less competitive relationships in co-occurrence networks to acclimate to warming. Conclusion Microbiome composition controlled carbon and nitrogen metabolism, thus determining lake microbial communities’ adaptability to heat stress. This study extended our insights on the lake microbial community response and adaptability under warming drivers.


Introduction
As the core of the ecosystem, microorganisms are not only numerous but also have greater diversity than other organisms, which have made significant contributions to the biogeochemical cycles (e.g., nitrogen fixation, oxygen generation, and methane generation) (Shilky et al. 2023) and stability of ecosystem services (Fuhrman 2009;Graham et al. 2016).Understanding how climate changes determine the spatial and temporal diversity of microbial species is a central goal of microbial ecology.
Temperature, a main driving factor of global environmental changes, can affect the structures and ecological functions of microbiota (Bai et al. 2013;Gillooly et al. 2001).The response of microorganisms to warming varies greatly in both taxonomy and geography.For example, archaea methanogens are more abundant in lakes in warmer regions, and future methane emissions will be more controlled by archaea methanogens under global warming (Winder et al. 2023).Over the globe, the microbial carbon use efficiency is lower in lower-latitude regions, reflecting that microbes reduce the proportion of carbon allocated to organic synthesis to adapt to the high-temperature environment to maintain more energy required for metabolism (Tao et al. 2023).In the soil vertical profile, warming does not change the carbon use efficiency of microorganisms in surface and deep soil layers (Zhang et al. 2023a).From the secular variation of more than half a century, warming accelerates the growth, respiration, and uptake of microbial communities (Walker et al. 2018), while the short-term and local-area effects of temperature on ecosystems are still complex.On one hand, extreme warming could result in biodiversity loss and halting of biodiversity recovery (Haase et al. 2023), and weaken the ecosystem functions and services related to biodiversity (Guo et al. 2018(Guo et al. , 2019)).On the other hand, warming can alter microbial enzyme activity, accelerate microbial turnover, and promote carbon accumulation in the environment (Hagerty et al. 2014).The microbial network response could also be enhanced by warming, which results in higher community stability (Thébault and Fontaine 2010;Yuan et al. 2021).Even under stress conditions, well-adapted microbial communities have special mechanisms to maintain important biogeochemical processes (Lay et al. 2013;Zorz et al. 2019).At present, studies on the impact of warming on microbial communities are mostly focused on soil ecosystems, while only very few studies have a full understanding of the response of microorganisms to warming in aquatic environments.As one of the most sensitive ecosystems on Earth to environmental change (Dudgeon et al. 2006), lake ecosystems have species that are more susceptible to the negative consequences of environmental change than terrestrial and marine species (O'Reilly et al. 2015;Woolway and Maberly 2020).Therefore, understanding how microbial species respond to sustained climate change is crucial for assessing freshwater ecosystem vulnerability and avoiding potential biodiversity losses.
Previous studies have conducted many long-term in situ warming experiments to explore the warminginduced changes and mechanisms in microbial communities.Warming and concomitant decrease in soil moisture exacerbate the reduction of soil microbial biodiversity (Wu et al. 2022).In aquatic environments, warming-induced acceleration of organic carbon decomposition and the promotion of algae growth affect the response of microorganisms to environmental changes (Hu et al. 2024;Sarmento et al. 2016).In addition, the interaction between eutrophic water bodies and temperature can amplify the effects of warming and may also lead to high variability in the composition of planktonic bacterial communities (Cross et al. 2022;Ren et al. 2017).What's more, the decrease in microbial diversity caused by early warming will change the microbial community composition and lead to the enrichment of thermophilic taxa (Nottingham et al. 2022).However, microbial communities are highly complex under inhabiting intricate natural environments.In situ experiments would be influenced by various factors so the work on establishing causal relationships between microbial communities, environmental conditions, and ecosystem functions is difficult (Pascual-Garcia and Bell 2020).Other big databased research evaluated the microbial data and environmental variables at the global scale (Delgado-Baquerizo et al. 2016;Wu et al. 2019;Yang et al. 2019), while such inference of the correlations between species and climate conditions ignored potential biological mechanisms, such as biomass changes, species interactions, and so on.Laboratory microcosm experiments have a greater potential for understanding ecological communities than generally thought (Cadotte et al. 2005;Gilpin et al. 1986).The environment for laboratory community studies can be precisely controlled and easily replicated and can answer many fundamental ecological questions through the construction of rapidly growing species and communities (Jessup et al. 2004;Lawton 1998).
To investigate the adaptability of microbes and their metabolisms and functional potential in differential lake environments under repeated warming, we conducted a 20-day microcosm experiment to simulate the responses of river microbes to warming.According to the microbial activity monitoring and amplicon sequencing analyses, we identified the differences between microbial populations and communities in the ability to respond to warming.Combined with Biolog EcoPlate experiments and metabolomics, we identified which taxonomic groups may play a key role, and described the potential ecological important nutrient cycles that may exist.We investigated the adaptability of the microbial community level and the responsive changes of the microbial community function to repeated heat stress.The mechanistic insight identified by this study corroborates the stability and recovery dynamics of freshwater microbes under heat stress.

Samples collection and microcosm experiments
Due to the minimum temperature difference between the north and south of China, samples were collected during the summer (August 2020) at a depth of 0.5 m of 8 lakes in China (Additional file 1: Table S1) to present different microbial communities and physicochemical properties.All samples were transported on ice to Tianjin Key Laboratory of Environmental Remediation and Pollution Control in China and stored at 4 °C.
In this experiment, there was significant heterogeneity in the physicochemical property data of 8 lakes, making it difficult to elucidate the impact of environmental changes on microbial adaptability using natural lake samples.Microcosm experiments allow for the study of individual environmental change factors through simplified microbial ecosystems (Lü et al. 2022), providing strong support for studying the potential relationship between warming and microbial community changes.Microcosm experiments were guided by an original batch incubation (Valenzuela et al. 2018) to explore the adaptability of the microbial population.Based on meteorological data from the website of Chinese Weather (https:// lishi.tianqi.com/), we calculated the monthly environmental temperature at sampling sites in August from 2011 to 2020.Considering that the highest monthly temperature at sampling sites was 35 °C, we set the ambient temperature at 25 °C and the warming temperature at 35 °C for the experiment (Renes et al. 2020;Yang et al. 2023).Lake samples were cultured in 10 mL round bottom centrifuge tubes with sterile artificial freshwater medium (Additional file 1: Table S2; Reasoner and Geldreich 1985) with a mean light intensity of 3000 lx under a 12:12 h lightdark cycle.The culturing tubes were covered with sterile sealing films to avoid microbial contamination but allow gas exchange.Microbes were first cultured for 5 days initial stage (Stage 0) to allow the microbial population to reach the stationary phase (Qu et al. 2023).Then, experiments were performed in triplicates for two control and warmed groups for 4 stages (every 5 days was one stage) and were kept under sterile and nutrient-replete growth conditions.During incubation, the microbial inoculum was homogenized by shaking centrifuge tubes every day.At the end of each stage, count living microbial cells with the flow cytometer (Accuri ™ C6 Plus, BD, USA) to ensure the initial cell density in about 1 × 10 3 cells mL −1 .Then transfer a small amount of cultures (approximately 5-30 μL) into a fresh medium to restart the growth in a new stage.The control groups were cultured under constant temperature at 25 °C, while warmed groups were under temperature stress (the temperature was 35 °C during stage 1 and stage 3) (Fig. 1A).In warmed groups, communities were at ambient temperature after hightemperature condition to assess whether microbes from the previous stage could recover from disturbance to an appropriate population growth state (Fig. 1A).
After the 20-day microcosm experiment, the natural logarithm of response ratio (RR) was used to measure the impact of experimental warming on microbes (Zhang et al. 2023b).RR was defined as the microbial living cells in the warmed groups divided by those in the control groups.If the lnRR value is greater than zero, it indicates that the experimental warming has a promoting effect on microbial growth.Previous researchers distinguished freshwater microbial communities based on net growth Fig. 1 Experimental design.A The microcosm experiment began with the pre-incubation (stage 0) of lake samples.Then the culture was separated as the initial communities into control and repeated warmed groups for 20 days of incubation (stages 1-4).B According to the changing trend of microbial biomass and response ratio (see Additional file 1: Figure S1), two types of microbial communities (vulnerable and adaptive microbial communities) were selected in the microcosm experiment.In four-line charts, the gray column represents temperature variation; the solid lines and dotted lines represent microbial biomass trends rate under temperature pulse conditions (Qu et al. 2023), and we defined the adaptability of microbial communities under experimental warming using the positive or negative response rate (Additional file 1: Fig. S1).In the four stages, if the lnRR values of two or more stages are greater than zero, it indicates that the microbial community has recovery and adaptability to experimental warming.On the contrary, it indicates that the microbial community is vulnerable and difficult to recover to the original state.Thus, the adaptive classification of microbial communities in our experiment was divided into two types: adaptive and vulnerable (Adger 2006;Brooks and Adger 2005).After incubation, two representative groups with significant differences were chosen to represent vulnerable and adaptive communities (Additional file 1: Fig. S1).3.75 mL of cultures in each centrifuge tube (n = 54) was collected and centrifuged, and collected sediment was stored at -80 °C for DNA extraction.Another 2.0 mL of cultures were collected, centrifuged, washed once by phosphate-buffered saline (PBS), and centrifuged, and the collected sediment was stored at − 80 °C for metabolites extraction.A total of 54 samples were constructed in this experiment, including 3 replicated samples for each of the two types of lake microbial communities at the initial stage, and 24 samples were collected in 4 stages of microcosmic experiment (3 × 4 controls and 3 × 4 repeated high-temperature treatments).

DNA extraction and amplicon sequencing
Microbial DNA was extracted from samples using the E.Z.N.A. ® Soil DNA Kit (Omega Biotek, Norcross, GA, USA).Using a NanoDrop 2000 UV-Vis spectrophotometer (Thermo Scientific, Wilmington, USA) to concentrate and purify DNA, with the DNA quality checking by 1% agarose gel electrophoresis.
The bacterial 16S rRNA and the fungal ITS genes were amplified with the primers 338F/806R and ITS3F/ ITS4R, respectively.Thermal cycling conditions were as follows: initial denaturation for 3 min at 95 °C, followed by 27 cycles for 30 s at 95 °C, annealing for 30 s at 55 °C and elongation for 45 s at 72 °C, with final extension for 10 min at 72 °C.The triplicate PCRs were performed, and PCR amplification was carried out in a 20 μL mixture according to the manufacturer's protocol.Amplicons were extracted from 2% agarose gels and further purified.The purified amplicons were pooled in equimolar amounts and paired-end sequenced (2 × 300) by Majorbio Bio-Pharm Technology Co., Ltd.(Shanghai, China) on a platform (Illumina, San Diego, USA).QIIME (version 1.17) was used to process raw FASTQ files.
After demultiplexed, quality-filtered by Trimmomatic and merged, the raw 16S rRNA and ITS gene sequencing reads generated 2,843,696 high-qualities sequences with 434 bp in length on average for 54 bacterial samples and 4,147,841 sequences with 319 bp in length on average for 54 fungal samples.Operational taxonomic units (OTUs) with a 97% similarity level were clustered by UPARSE (version 7.0.1090).RDP Classifier (version 2.11, with a confidence threshold of 0.7) analyzed the taxonomy of each representative OTU sequence against the Silva database (Release 138) and fungal OTU sequences against Unite database (Release 8.0).Finally, the microbial amplicon reads were normalized based on the minimum data of samples.

Microbial metabolic capacity experiments and analysis
The carbon metabolic fingerprints of the microbial community were indicated by EcoPlates (Biolog Inc., USA).The EcoPlates consist of 96 wells with 31 different carbon sources.Microbial suspensions were dispensed into 96 wells (150 μL per well with the initial microbial concentration of 1 × 10 4 cells mL −1 ) with another 3 blank samples.Then the incubation was at 25 °C in the dark for 5 days, and the optical denticity at 590 nm (OD 590 ) of each well was read with a microplate reader (Synergy H4, Bio-Tec, USA) in 0, 24, 48, 72, 96, and 120 h.Average well color development (AWCD) change rate of each hole on the microplates was used to evaluate the ability of microbial communities to utilize a single carbon source, i.e. the overall metabolic activity of microorganisms.Simpson's diversity index (D) was calculated to evaluate the dominance of the most common species and to further analyze the functional diversity of communities.The equations were as follows: where R i is the OD 590 of each well on the EcoPlates, R 0 is the OD 590 value of the control well (if R i -R 0 ≤ 0, R 0 is regarded as 0), and N is the number of carbon sources in 32 holes on the EcoPlates, i.e.N = 31.
In the above formula, P i is the ratio of the difference of OD 590 between the hole with carbon source and the control hole to the total difference of the whole plate, i.e.
. R s i is the standardized absorbance value.Based on the AWCD value, the utilization of various carbon substrates (at 120 h) by microorganisms was analyzed.

Metabolite extraction and LC/MS analysis
Refer to previous research (Meyer-Cifuentes et al. 2020), the metabolites in metabolomics analysis were extracted from the microbial pellets.Briefly, the cell pellets for metabolic extraction gained at the end of Sect.Samples collection and microcom experiments were resuspended in methanol and removing residual plastic film/ particulate material.The supernatant was concentrated by using a vacuum centrifugal concentrator (SCIENTZ-1LS, Scientz, China).After two-step derivatization of adding methoxyamine hydrochloride dissolved in pyridine (20 mg/mL, 50 μL) and N-methyl-N-(trimethylsilyl) trifluoroacetamide (80 μL), the injected samples were analyzed by gas chromatography mass spectrometry (GC-MS, 6890A/5977A, Agilent, USA).Metabolic data was performed by a platform, called MetaboAnalyst (version 5.0, https:// www.metab oanal yst.ca).Metabolites in significant differences were screened by volcano plots (FC > 2.0 or < 0.5) and T test with P < 0.05.

Statistical analyses
The microbial community structure composition at different taxonomic levels (e.g., phylum, family, and OTU) was performed by R (version 3.6.3).Qualified OTU data were used to calculate Chao1 index (microbial richness) and Shannon index (microbial diversity) by Mothur (version 1.30.2).Based on the normalized intensities of OTUs, the distribution of top OTUs in different groups was shown by enrichment charts made on an online tool ImageGP (Chen et al. 2022), and used hierarchical clustering analysis of pairwise Bray-Curtis distances to explore the similarity and difference among different groups of samples.Beta diversity (β-diversity) analysis was assessed and mapped by principal coordinate analysis (PCoA) and partial least squares discriminant analysis (PLS-DA).Using linear discriminant analysis (LDA) coupled with effect size measurements (LEfSe) analysis to describe the significant differences from phylum to OTU level.The OTUs data of bacteria and fungi were annotated with function information by matching with the functional annotation of prokaryotic taxa (FAPROTAX) database (Louca et al. 2016) and the FUNGuild database (http:// www.fungu ild.org/).Co-occurrence analyses were constructed using the R package "Hmisc".To reduce complexity, OTUs with a total abundance more than 50 were retained, and the OTU-based pairwise Spearman's correlations were calculated (with the absolute value of a correlation coefficient > 0.6 and P < 0.05), then networks were visualized by Gephi (version 0.9.7).Mantel-test analysis was conducted using the R package "ggcor" and Spearman correlation analysis was used to present the relationship between microbial taxa and differential metabolites, which was performed by clustering correlation heatmaps using the OmicStudio (https:// www.omics tudio.cn).
All experiments in our study were performed in triplicate, with the results as the means ± standard deviation (SD).The data was tested for normality and variance homogeneity, and then, data analysis was processed by IBM SPSS Statistics software (version 22.0) and Graphpad prism (version 8.0.2).The significant differences (P < 0.05) between the control and warmed groups were calculated by T-test, and one-way ANOVA test was used to seek significant differences between multiple groups.

Changes in microbial community structure under warming
In the microcosm experiment, warming significantly influenced the relative abundance of the vulnerable community (P < 0.0001) with plummeting in microbial numbers, while that of adaptive communities increased relatively after disturbance (Fig. 1B).Compared with the control, microbial biomass of both communities was in great hit under the first high-temperature disturbance (Fig. 1B, stage 1).The microbial biomass of the vulnerable community was near to collapse and hard to recover during cushioning periods (hollow circles and dash lines in red, in stage 2 and stage 4), while the biomass of the adaptive community was in steady-state growth after warming (hollow circles and dash lines in blue).
The 16S and ITS rRNA gene amplicons identified 142 and 212 taxa in the 54 samples of bacterial and fungal communities, respectively.In general, the abundance and diversity of bacteria were lower than that of fungi (Fig. 2A, B).For OTU-level profiling and hierarchical clustering analysis (Additional file 1: Fig. S2), the differences caused by high temperature to bacteria were more differentiable than those to fungi.This result was also reflected in partial least squares discriminant analysis (PLS-DA) and β-diversity (Additional file 1: Fig. S3).
Both types of initial microbial communities were mainly represented by phyla Proteobacteria (Fig. 2C).Warming had no significant impact on the community structure of adaptive communities, but greatly altered the vulnerable community, whose dominant phylum turned to Firmicutes (with abundance from 17.00% to 99.38%).At the family level of bacteria, the majority of phyla Firmicutes were dominated by the family Paenibacillaceae in control groups, while dominated by the family Bacillaceae in warmed groups (Fig. 2C, Additional file 1: Fig. S4A, C; P < 0.00001).It can be inferred from the changes in the community composition of the adaptation group at the family level (Fig. 2C) that the community biomass fluctuations were primarily driven by the abundance of the bacterial family Pseudomonadaceae, whose relative biomass plummeted in vulnerable groups (Additional file 1: Fig. S4A; 82.22% vs. 0.51% in controls and treatments, respectively; P < 0.001) but increased in adaptive groups (Additional file 1: Fig. S4C; 29.25% vs. 57.53% in controls and treatments, respectively; P < 0.05) after warming.The linear discriminant analysis (LDA) effect size (LEfSe) cladograms showed more detailed changes and differences in community structure after repeated warming from phylum-level to OTU-level (Additional file 1: Fig. S5, S6).

Changes in microbial network and carbon metabolic capability
Warming had different effects on the complexity of the associated networks in different microbial communities.After warming disturbance, the microbial network complexity of adaptive communities decreased with lower richness and connectivity but the cohesiveness increased (0.55 vs. 0.82 cohesiveness in control and warmed samples, respectively), while the vulnerable community networks exhibited opposite characteristics (Fig. 3A and Additional file 1: Table S3).What's more, all the relative links in the adaptive community after warming were positive links (Additional file 1: Table S3).
Warming disturbance had a greater impact on bacterial function than fungi.The Faprotax function prediction results for bacteria showed significant variation among the treatments for over half of the bacterial functions (Fig. 3B, P < 0.05), while the FUNGuild function prediction for fungi showed significant differences only in two functional groups between control and warmed groups (Additional file 1: Fig. S7).After the high-temperature disturbance, the bacterial functional groups like nitrogen respiration and nitrate respiration significantly decreased in vulnerable groups (Fig. 3C, P < 0.05) but showed no significant changes in adaptive groups (P > 0.05).There were some functional groups, such as dark hydrogen oxidation, nitrite respiration, and so on, with significant small increases in both two microbial communities after warming (Fig. 3C, P < 0.01).
According to the Biolog results, the utilization of carbon sources fluctuated in different microbial communities under high-temperature disturbance.The results of the corresponding ratio of heatwaves to microbial carbon metabolism activity showed that warming decreased the microbial diversity and carbon metabolic capacity of vulnerable groups (Fig. 4A, B).Warming had slightly different effects on carbon metabolic capacity in the adaptive groups, with statistically significant at stage 4 (Fig. 4A).High temperature impaired the metabolic capability of labile carbon in vulnerable groups (Fig. 4C), while enhancing such capability in adaptive communities (Fig. 4D).

The intracellular metabolite changes of microbial communities
In revealing microbial metabolism associated with microbial communities, Mantel test analysis (Fig. 5) suggested that the significant correlations between bacteria and metabolites were more than that between fungi and that the majority of correlations were positive correlations between metabolites (80.95% of all 231 correlations).To further clarify these associations, we constructed links between representative differentially abundant metabolites and microbial species described earlier.In the correlation data set of differential metabolites and microbes, only 11.57% and 5.71% of all associations with bacteria and fungi, respectively, were significant (Additional file 1: Fig. S7, P < 0.05).Of these statistically significant associations, the connection between differential metabolites and bacterial communities was stronger than that of fungal communities.The correlation heatmap and microbial relative abundance bars indicated that the relative abundance of microbial species had little effect on such correlations and that the significant relationships between microbes and metabolites were mostly concentrated in metabolites related to amino acid metabolism pathways (Additional file 1: Fig. S8).
Under high-temperature disturbance, the types of differential metabolites in vulnerable groups were far more than those in adaptive groups (Fig. 6A, stage 1 and stage 3).What's more, the significant effects of high-temperature on the metabolic pathway were mainly concentrated on the metabolism of amino acids.High-temperature generally downregulated the amino acid metabolism in vulnerable communities, but upregulated the amino acid metabolism in adaptive communities (Fig. 6B).The differential metabolites related to carbohydrate and energy metabolism only appear during the recovery phase after high-temperature disturbance (Fig. 6A, stage 2 and stage 4).

Bacterial community is less stable than fungal community under warming
The warming-induced changes in the microbial communities are different, most intuitively reflected in their impact on the structure of microbial communities.Our study verified that the fungal community structure was much more stable than the bacterial community in the face of high-temperature disturbance.In our experiment, the overall richness of fungi was higher compared to bacteria, and species-rich communities were more able to maintain their ecosystem functions and had higher resistance to fluctuating environmental conditions (Bestion et al. 2021).The changes in community structure also indicate that bacteria are more sensitive to short-term environmental temperature changes and have a faster evolutionary response (Huang et al. 2021).Although short-term warming had no significant impact on bacterial and fungal OTU richness, long-term warming can decrease the richness of bacteria and fungi (Wu et al. 2022).
The dynamic changes in community structure were partly due to the direct impact of temperature on the growth status of microorganisms, confirmed research that Proteobacteria is the main core species in colder lakes, and the population gradually decreases and eventually disappears as the temperature increases (Guo et al. 2023).For phyla Firmicutes, as the optimum growth Fig. 6 The differential metabolites and metabolomics analysis of microbes under warmed disturbance.A Relative abundance of differential metabolites between groups.These metabolites were mainly engaged in carbohydrate, energy, amino acid, lipid, and nucleotide metabolism.The abundance of differential metabolites (**P < 0.01, *P < 0.05) was standardized by Z-score.The size of the black circles represented the relative abundance of differential metabolites, which was the proportion of the substance to the total abundance of all differential metabolites in the current stage.Stages from 1 to 4 represent temperature conditions in treatment at 35, 25, 35, 25 °C, respectively.B Metabolic pathway analysis.Metabolites marked green were differential metabolites, and those indicated by arrows were those significantly different between the control and warmed groups.The number below the arrow represents the stage that the metabolite changed temperature of the family Bacillaceae was higher than that of Paenibacillaceae and they can also survive in harsh environmental conditions, so family Bacillaceae dominated the phyla Firmicutes after warming (Kakagianni 2022; Patowary and Deka 2020).Other studies have shown that under competitive conditions, Bacillaceae has a more adaptive advantage than Paenibacillaceae (Sun et al. 2022).Another reason might be the indirect influence of temperature on growth conditions.Because the bacterial family Pseudomonadaceae's optimal growth temperature is above 20 °C (Chakravarty and Anderson 2015) and they are not acid-tolerant and unable to grow below pH 4.5 (Dodd 2014), we inferred that the temperature indirectly affected the acidity of the environment, leading to changes in its relative abundance.Therefore, community structure changes to drive better adaptation of microorganisms to warming.For the adaptive microbial community, the community structure can adapt to the changes in environmental conditions in a short time, while the vulnerable community needs a longer time to recover.Although the biomass of vulnerable communities did not recover in this short-term experiment, as species with higher heat tolerance limits became more abundant, the temperature range of biological activity may be amplified or transferred (Jabiol et al. 2020;Manning et al. 2018).

Highly stable microbial communities adapt to warming by positive species correlations and adaptive functional changes
In the microbial network, the "clustering phenomenon" (low connectivity and high cohesiveness) implied that microbes altered the clustering mode to maintain microbial function under high-temperature stress (Wang et al. 2023a).What's more, taxonomic groups tend to construct a dominant position of positive correlations to face high-pressure environments (Hernandez et al. 2021;Wang et al. 2023b).Such positive correlations between nodes in microbial co-occurrence networks could be the result of functional interdependencies among microbial communities under extreme environmental conditions (Mandakovic et al. 2018;Neilson et al. 2017).Our study also showed that microbes in adaptive communities tend to establish positive correlations to adapt to experimental warming.
The results of the bacterial functional prediction showed that warming significantly impaired nitrogen respiration and nitrate respiration in vulnerable groups, evidenced by the inhibition of the growth and activity of denitrifying communities.Our results are consistent with previous studies which found that the higher temperature ascribed to climate change might decrease the reduction of nitrate in the lakes (Jiang et al. 2023b).It has been reported that denitrifiers display a notable sensitivity to changes in temperature (Kuypers et al. 2018;Palacin-Lizarbe et al. 2018).In this study, we also observed that a significant decrease in aerobic denitrifiers Pseudomonadaceae after high temperature led to such responsive changes in nitrogen metabolism.Because of the diverse forms of nitrogen in the process of nitrogen metabolism (Gruber and Galloway 2008;Ollivier et al. 2011), nitrogen cycle-based metabolic functions are particularly vulnerable to environmental fluctuations, such as temperature, nutrient, pH and dissolved oxygen (Jiang et al. 2023a;Palacin-Lizarbe et al. 2018).The stability of microbial communities has an impact on their maintenance of ecosystem functions because highly stable microbial communities (strong resistance to warming) play a more important role in preserving the initial "functional optimal" communities and ecosystem functions (Hammill et al. 2018;Pasari et al. 2013).

Moderate warming enhances the metabolic capacity of microbial communities by stimulating fungal metabolism
Warming can affect the carbon use efficiency of microbes through changes in microbial growth and the microbial community structure (Dijkstra et al. 2015).Warming significantly weakened the microbial biomass of the vulnerable communities, resulting in the decrease of degrading enzymes, which leads to a decrease in carbon metabolism.Microbial adaptation or changes in microbial communities may lead to an increase in carbon use efficiency, offsetting the decline in microbial biomass (Allison et al. 2010).Generally, the decomposition of carbon increases in warming conditions (Hu et al. 2024).From the view of the carbon sources, molecules with higher energy availability are easier decomposed by microbes as the temperature increases (Williams and Plante 2018).The energy source of most microorganisms is labile carbon (Huang et al. 2021;Nottingham et al. 2019).Warming at suitable temperatures can increase the rapid utilization of labile carbon (Sullivan et al. 2020).
The positive correlations between bacteria and intracellular metabolites indicated a possibility that the metabolite promoted the growth of the species, or that the metabolite was generated from the species (Franzosa et al. 2019).However, most differential metabolites tended not to be linked mechanistically to most microbes, and vice versa (Franzosa et al. 2019).High temperatures reduced the metabolic capacity of vulnerable communities because the significant decrease in microbial biomass caused by excessively high temperatures has reduced microbial respiration and metabolic capacity.Nevertheless, warming increased the metabolic capacity of the adaptive community, because a 5-10 °C warming stimulated fungi metabolism (Ferreira and Chauvet 2011) and improved the nutrient utilization efficiency and organic matter turnover compacity of microbes (Fernandes et al. 2014).The structure of bacterial communities (Wan et al. 2023) and labile carbon were major factors in the changes in bacterial metabolism, while recalcitrant carbon was mainly influencing the fungi community (Wang et al. 2022).Another reason that the warming-induced effect on different communities was different in microbial metabolism and function is that the environmental temperatures of different microbial communities are inherently different, and the optimal growth temperature dominates the growth and growing states of species (García et al. 2018).Therefore, there is a certain temperature preference range for specific species within the community.The adaptive community with high species richness (Albrecht et al. 2021) and high stability would contribute more to resisting environmental changes and maintaining ecosystem functions than vulnerable communities (Pasari et al. 2013).

Conclusions
In our study, the microbial community originated from the freshwater environment in nature, and the extent to which the two types of community classifications determined in the experiment can be transferred to more complex natural communities remains a challenge.Nevertheless, our research showed that lake bacteria were more likely to be impacted by warming than fungi based on the short-term response of microbial communities under fluctuating environmental temperature.Bacterial communities were characterized by low stability under disturbance, such as biomass reduction, community composition dynamics, species diversity decreasing, and microbial function changes.Although the exact nature of the predicted interactions between bacteria and fungi has not been fully elucidated, our data revealed that microbial communities with higher adaptive ability have closer bacterial-fungal connections and more positive correlations to be interdependent to adapt to disturbances.We also found that the degree of influence of warming varies in the microbial community.For example, in China, the resilience and resistance of lake microbial communities in colder areas are far inferior to those in warmer areas when facing the same degree of heatwaves.As temperature's indirect effects varied between microbial communities, these responses reflect the differences in the functionality and metabolic activity of microbial communities to some extent.This study provides useful information about the effect of short-term experimental warming on lake microbial communities and provides more insights to fully understand microbial adaptive

Fig. 2
Fig. 2 The taxonomic composition of bacteria and fungi.A The mean richness of bacterial and fungal OTUs, evaluated by the Chao1 index, was mainly shown by the gathering of control (n = 12) and warmed groups (n = 12).Statistical significance between the two groups was calculated by the Wilcoxon rank-sum test, showing no significant difference.B The differences of microbial diversity between control and warmed were analyzed by Student's t-tests (***P < 0.001, **P < 0.01, *P < 0.05).Bacterial (C) and fungal (D) communities of vulnerable and adaptive communities at phylum and family levels under initial (Init.),control, and repeated warmed groups.Stages from 0 to 4 represent temperature conditions in treatment at 25, 35, 25, 35, 25 °C, respectively

Fig. 3
Fig. 3 The cooccurrence networks and function prediction for microbial communities.A The nodes in different colors indicated individual bacterial or fungal OTU, colored according to phylum.Larger nodes indicate the OTU was relatively more abundant within that group.The lines represented significant Spearman correlations (|R|< 0.6 and P < 0.05) between OTUs.The gray and black lines indicated positive and negative correlations, respectively.B Functional characteristics of bacterial samples (n = 54) were predicted by Functional Annotation of Prokaryotic Taxonomy (FAPROTAX).C Pairwise comparisons and significant testing (***P < 0.001, **P < 0.01, *P < 0.05) for vulnerable (red bar) and adaptive (blue bar) communities

Fig. 4
Fig. 4 The warming effect on carbon metabolic capacity.A, B The microbial carbon metabolic activity (evaluated by AWCD) and diversity (evaluated by Simpson's diversity index) at 120 h in the Biolog experiment.Response ratio (RR) was calculated as warmed/control.C, D Carbon metabolic capacity of different types of carbon sources.The differences between groups were analyzed by one-way ANOVA (***P < 0.001, **P < 0.01, *P < 0.05), marked by black stars.E The carbon metabolic capacity of three labile carbon sources.The change in carbon metabolism was calculated as (warmed − control)/control × 100%.Error bars represent standard error (n = 3).The differences between control and warmed were analyzed by Student's t-tests (***P < 0.01, **P < 0.01, *P < 0.05), and marked by red stars.Stages from 1 to 4 represent temperature conditions in treatment at 35, 25, 35, 25 °C, respectively

Fig. 5
Fig. 5 Variations of microbes and metabolites were analyzed by Mantel-test analysis.Edge width corresponds to Mantel's r statistics, edge color represents the statistical significance, and color gradient indicates Spearman's correlation coefficient mechanisms under climate warming and to illustrate how these ecosystems are at risk from climate change.This work was financially supported by the Ministry of Science and Technology of People's Republic of China as a key technology research and development program project (Grant number 2023YFC3709001), the Ministry of Education of People's Republic of China as a discipline innovation and intelligence introduction project (Grant number B17025), and the Tianjin Science and Technology Bureau as a key science and technology supporting project (Grant number S19ZC60133).