Rhizosphere Fungal Dynamics in Sugarcane during Different Growth Stages

Understanding the normal variation of the sugarcane rhizosphere fungal community throughout its life cycle is essential for the development of agricultural practices for fungal and ecological health associated with the microbiota. Therefore, we performed high-throughput sequencing of 18S rDNA of soil samples using the Illumina sequencing platform for correlation analysis of rhizosphere fungal community time series, covering information from 84 samples in four growth periods. The results revealed that the sugarcane rhizosphere fungi possessed the maximum fungal richness in Tillering. Rhizosphere fungi were closely associated with sugarcane growth, including Ascomycota, Basidiomycota, and Chytridiomycota, which showed high abundance in a stage-specific manner. Through the Manhattan plots, 10 fungal genera showed a decreasing trend throughout the sugarcane growth, and two fungal genera were significantly enriched at three stages of sugarcane growth (p < 0.05) including Pseudallescheria (Microascales, Microascaceae) and Nectriaceae (Hypocreales, Nectriaceae). In addition, soil pH, soil temperature, total nitrogen, and total potassium were critical drivers of fungal community structure at different stages of sugarcane growth. We also found that sugarcane disease status showed a significant and strong negative effect on selected soil properties by using structural equation modeling (SEM), suggesting that poor soil increases the likelihood of sugarcane disease. In addition, the assembly of sugarcane rhizosphere fungal community structure was mainly influenced by stochastic factors, but after the sugarcane root system became stable (Maturity), the stochastic contribution rate decreased to the lowest value. Our work provides a more extensive and solid basis for the biological control of sugarcane potential fungal diseases.


Introduction
The microbial community in the rhizosphere environment is critical for the health of terrestrial plants and sustainable soil development [1]. As more attention is paid to the relationship between plants and soil microbiomes, comprising the interactions between microbial and plant root systems and genetic components [2][3][4], it is apparent that more exploration is needed of what constitutes normal temporal variations in rhizosphere microbial community structures and compositions. In particular, research related to rhizosphere fungal communities that may produce soil-borne diseases of plants is essential. Variations in the rhizosphere fungal group and between different plants and life-cycle stages have a lot of space for exploration. As we all know, after germinating in the soil, plant roots begin to recruit microbial communities closely related to their growth through the accumulation of root secretions, rhizosphere deposition, and nutrient uptake. Thus, the growth activities of plants during different periods cause changes in rhizosphere fungal communities and functions. For instance, model plants (Arabidopsis thaliana and Medicago truncatula) can retain certain fungal populations in the soil through certain chemicals secreted by their root systems [5]. In addition, advances in sequencing technology provide a foundation for research on the dynamic changes of microbes in the rhizosphere of plants in space and time [6]. Similarly, these types of studies are needed to understand the immigration and emigration patterns of microorganisms between different parts of the plant, between different points of time and between different plants and soil environments [7][8][9], especially the related research on time series. For example, Zhang et al. (2018) identified biomarker taxa and established a model to correlate root microbiota with rice resident time in the field using a machine learning approach [10], and Li et al. (2020) analyzed the variation of the fungal community structure in the rhizosphere of Gardenia, which provided a new theoretical framework for further study on the rhizosphere mechanism [11].
Furthermore, as an important economic crop of food and biofuel, sugarcane is constantly expanding its planting area under the condition of global warming [12,13]. The study of sugarcane rhizosphere fungi temporal variation has great significance in better understanding the interaction between sugarcane and rhizosphere microbiota [14]. For example, sugarcane growth is facing a huge threat due to Yellow Canopy Syndrome (YCS) in Australia. YCS is a kind of largely undiagnosed plant disease that is impacting sugarcane growth across Queensland, Australia, causing huge yield losses [15]. Studying the temporal variation of sugarcane rhizosphere microorganisms provides a new way of thinking about the solution of YCS. Additionally, given the growing importance of sugarcane [16,17], identifying ways to promote the healthy growth of sugarcane and sustainably increase productivity is critical, and hence, the potential to harness microorganisms from the sugarcane rhizosphere has recently gained more attention [18][19][20]. For instance, rhizosphere microorganisms are involved in the solubilization of phosphorus and potassium containing minerals [21]. Sugarcane smut disease, caused by a fungus called Sporisorium scitamineum (Basidiomycota, Ustilaginales), is a limiting factor to cane production and a potential threat to the sugar industry. Juma and Musyimi's research showed that the selected isolates from sugarcane rhizosphere microflora had an evidently antagonistic activity against the Sporisorium scitamineum, and is recommended as a potential biocontrol agent for this pathogen [22]. In a related study, She et al. (2021) found maize rhizosphere microorganisms are critical for facilitating microbiome bioremediation for soil affected by neutral-alkaline mining [23]. Moreover, understanding how the sugarcane rhizosphere fungal population change with stages would provide valuable data and could in the long-term result in the identification of potential management strategies for healthy cultivation of sugarcane, including an improved breeding process taking the microbe into account or better targeted biological prevention and treatments [24]. Here, we used sugarcane as an experimental plant to fill critical knowledge gaps that will accelerate our exploration process to successfully harness rhizosphere fungal groups to sustainably enhance sugarcane productivity. This also has reference significance for other crops in nature. Our study had three main objectives: (i) to determine the drivers of the sugarcane fungal community under field conditions; (ii) to understand the seasonal trends of fungi in the rhizosphere of sugarcane; and (iii) to explore the assembly and interaction relationships of rhizosphere fungal communities during sugarcane development. To achieve this, we examined fungal community assemblages of sugarcane rhizosphere soil through Illumina amplicon sequencing of the samples collected from young to mature sugarcane (Seeding, Tillering, Elongation, and Maturity). Our results provide more foundational insight into the process of sugarcane root microbiota exploration.

Sugarcane Rhizosphere Nutrients Produce Differential Changes at Different Stages
The results showed that the height and stem diameter of sugarcane plants increased continuously over time, with the rate of diseased plants reaching their highest (19.8%) in Elongation (Table 1). We also investigated the temporal variation of soil biochemical properties collected in Seeding, Tillering, Elongation, and Maturity during the different growth stages of sugarcane. Our analysis showed that soil biochemical properties varied with time. Soil properties such as soil pH and available nitrogen (AN) decreased considerably (p < 0.05), with Elongation and Maturity recording the lowest soil pH, while Seeding had recorded the highest, followed by Tillering. In addition, the soil organic matter (OM) and total phosphorus (TP) were significantly higher in Elongation compared to Seeding, Tillering, and Maturity (p < 0.05). We also observed that Seeding recorded a significant amount of soil total nitrogen (TN) relative to the other stages. Compared to Elongation, a considerable amount of soil available phosphorus (AP) was recorded in Seeding and Tillering, and the available phosphorus nutrient decreased by 20.07% and 25.37%, respectively. In comparison to Tillering and Maturity, soil total potassium (TK) revealed significant improvement (p < 0.05) in Seeding, followed by Elongation. However, soil available potassium (AK) also showed significant differences among the different time points, with the lowest AK values in Maturity (p < 0.05), which is significantly reduced by 24.68% compared with the highest value. Soil temperature peaked significantly in Tillering and Elongation than in Seeding and Maturity (Table 2).

Rhizosphere Fungal Diversity Differs between Different Stages and Correlates with Soil Properties
The total number of sequences in the raw sequencing data is 4,592,374 with an average length of 401.21. The shortest sequence length is 221 and the longest sequence length is 453. For details, see Annexes-1. Fungal diversity (Shannon and Simpson), richness (ACE), and Sobs (number of observed species) were assessed during the different time points. It was observed that fungal diversity increased substantially (p < 0.05) in Tillering compared with the other stages (Table 3 and Figure S1A). Moreover, compared to Maturity, fungal richness and Sobs were significantly enriched (p < 0.05) in Tillering, followed by Seeding and Elongation (Table 3, Figure S1C,D). In addition, the dilution curve showed a sufficient amount of sequencing data ( Figure S2). We carried out Principal Coordinates Analysis (PCoA) to explore and visualise similarities or dissimilarities in fungal community composition in all the samples collected during the different growth stages of sugarcane. Based on the growth changes of FN 41 sugarcane throughout its life cycle, simulation map with distinctive features for different stages were constructed ( Figure 1A). The analysis revealed that fungal community composition was clustered together in the different samples, and the soil samples shifted with different stages in the first axis ( Figure 1B).  The relationship between soil environmental variables and fungal alpha diversity was investigated using Pearson's correlation coefficients (Figure 2A). The analysis demonstrated that soil temperature (Tem) exhibited a strong and positive association with fungal Sobs and ACE, Shannon and PC1. However, soil Tem was negatively associated with PC2. Furthermore, soil pH revealed a significant positive correlation with PC1, Sobs, and ACE. Soil TN and AN were significantly and positively associated with PC1. Soil TN also exhibited a positive relationship with PC2, while TK had a strong positive association with PC1 and PC2. We also noticed that soil TP demonstrated a positive correlation with PC2. On the other hand, soil TN and TK were negatively related to Shannon. Meanwhile, regression analysis was conducted to further confirm the association among soil environmental variables, fungal ACE (richness), Shannon (diversity), PC1 and PC2 during the different time points ( Figure S3). The analysis demonstrated that soil pH and AN showed a positive relationship with fungal richness. Moreover, soil Tem exhibited a positive relationship with fungal richness, particularly in Tillering and Elongation ( Figure S3A). Soil pH and C/N showed a positive correlation with fungal richness in Tillering, whereas AN showed a strong positive correlation with fungal richness in Seeding. In Tillering and Elongation, soil Tem revealed a strong positive relationship with fungal richness ( Figure S3B). Furthermore, soil pH, TN, AP, AN, and AK had positive relationships with PC1 in Seeding, while in Tillering and Elongation, soil Tem was positively related to PC1 ( Figure S3C), whereas in Seeding, TN and TK had a positive association with PC2, while in Elongation and Maturity soil TP and AK demonstrated a positive correlation with PC2, respectively ( Figure S3D).

Variability in Fungal Community Composition at Different Growth Stages
Fungal relative abundance during the stages was assessed at the phylum level. The analysis showed that Ascomycota accounted for more than 50% of the absolute dominance in the entire sample present in the rhizosphere soil of sugarcane, while Basidiomycota, Chytridiomycota, and Glomeromycota were present in smaller dominant groups in the

Variability in Fungal Community Composition at Different Growth Stages
Fungal relative abundance during the stages was assessed at the phylum level. The analysis showed that Ascomycota accounted for more than 50% of the absolute dominance in the entire sample present in the rhizosphere soil of sugarcane, while Basidiomycota, Chytridiomycota, and Glomeromycota were present in smaller dominant groups in the hot zone of sugarcane roots. However, Chytridiomycota and Glomeromycota were significantly higher (p < 0.05) in Tillering than the other stages ( Figure S4A,B). Basidiomycota increased significantly in Tillering and Seeding compared to other phases and the mean sample abundance peaked in Tillering ( Figure S4C). Venn diagram analysis demonstrated that 15, 23, 11, and 11 unique fungal OTUs were detected in Seeding, Tillering, Elongation, and Maturity, respectively ( Figure 3B). Moreover, among the soil samples collected in the various stages, 451 OTUs were shared among the different time points. In Tillering, the highest OTUs were recorded compared with the other sampling times.
Basidiomycota increased significantly in Tillering and Seeding compared to other phases and the mean sample abundance peaked in Tillering ( Figure S4C). Venn diagram analysis demonstrated that 15, 23, 11, and 11 unique fungal OTUs were detected in Seeding, Tillering, Elongation, and Maturity, respectively ( Figure 3B). Moreover, among the soil samples collected in the various stages, 451 OTUs were shared among the different time points. In Tillering, the highest OTUs were recorded compared with the other sampling times.

Fungal Community Composition, FUNGuild, and Soil Physicochemical Properties Are Interrelated
The different environmental variables were significantly correlated with the variation of selected soil fungal genera ( Figure 4A). RDA demonstrated the relationship between soil properties and fungal communities, which had an eigenvalue of 0.2846 for the first axis and 0.0802 for the second axis, respectively. The vectors indicated the following: total N, total P, total K, pH, available N and soil temperature played a greater role than

Fungal Community Composition, FUNGuild, and Soil Physicochemical Properties Are Interrelated
The different environmental variables were significantly correlated with the variation of selected soil fungal genera ( Figure 4A). RDA demonstrated the relationship between soil properties and fungal communities, which had an eigenvalue of 0.2846 for the first axis and 0.0802 for the second axis, respectively. The vectors indicated the following: total N, total P, total K, pH, available N and soil temperature played a greater role than available K, available P and organic matter for sugarcane rhizosphere fungi genera. The soil fungus community was dominated by Basidiomycota and Ascomycota, which showed stronger associations with higher pH, available N and soil temperature, while Chytridiomycota slightly revealed the opposite trend. We then used the Source Model of Plant Microbiome (SMPM) to estimate the proportion of sugarcane rhizosphere fungal communities from "adjacent period" and "unknown" sources and took the adjacent periods as the source and library in turn ( Figure 4B). The results showed that the related fungal communities mainly came from the transmission of adjacent periods, and the transfer proportion decreased gradually with the migration of time. In addition, Tremellales and Saccharomycetales, as the dominating orders in fungi, were positively associated with pH and AN. However, Sordariomycetes and Pseudallescheria showed a stronger negative association with pH and AN ( Figure 4D). To further explore the environmental driving factors of the sugarcane fungus community, we conducted a Mantel-test analysis and analyzed the correlation between the microbial matrix and the soil property matrix, and then we correlated distancecorrected dissimilarities of taxonomic and functional community composition with those of environmental factors. Overall, pH, TN and TK exhibited the strongest correlations with the taxonomic composition in the sugarcane rhizosphere soil ( Figure 4C), while no significant correlation was found for other soil factors (p < 0.05). In addition, soil temperature and pH were only weakly correlated with taxonomic and functional community composition. Besides, almost all the environmental variables exhibited no significant correlations, except for AN.

Fungal Communities Differ in Co-Occurrence Networks with Changes in Growth Periods
We employed Manhattan plots to examine fungal community composition OTUs differences in sugarcane rhizosphere soil during the growth stages ( Figure 5 and Table S1). The analysis revealed that Basidiomycota, Ascomycota, and Chytridiomycota were significantly more abundant in Tillering than in Seeding ( Figure 5A). Additionally, Ascomycota increased considerably in Tillering relative to Seeding; however, Basidiomycota and Chytridiomycota decreased profoundly in Seeding ( Figure 5B). In Elongation, Basidiomycota and Ascomycota increased significantly, while Chytridiomycota diminished considerably in Maturity ( Figure 5C). Meanwhile, Venn diagram analysis was adopted to gain a deeper understanding of the genera of fungal OTUs from one time point to another. The analysis demonstrated that between Seeding and Tillering, the highest number of enriched fungal OTUs was detected, followed by Elongation to Maturity ( Figure 5E). However, between Tillering and Elongation the highest amount of depleted fungal OTUs were identified, followed by Elongation to Maturity ( Figure 5D). Meanwhile, the co-occurrence network showed the interactions of rhizosphere fungi in sugarcane during four critical periods, where the fungi belonging to mainly Talaromycetes, Talaromyces, Fusarium, Sordariales, and Pseudallescheria were dominated ( Figure 5F). The role of each fungus in the network changed over time during different reproductive periods, with the degree and importance changing. Among the four critical fertility networks, Seeding had the largest average density and mean degree, followed by Maturity (Table S9).

Construction of Structural Equation Models Related to Pokkah Boeng Disease of Sugarcane
We constructed an SEM to assess the direct and indirect effects of soil properties (pH, AN and TN), microbial genera (Fusarium and Talaromyces), and microbial diversity (ACE and Shannon) on the incidence of sugarcane ( Figure 6 and S5 and Table S2). We used a multigroup modelling method to assess which relationships exist between soil properties, microbial communities, and sugarcane disease during the growing process. The results of model-1 showed that soil pH was significantly and directly affected the abundance of Fusarium, sugarcane disease and fungal diversity. However, soil pH negatively regulated the disease rate of sugarcane (the lower the pH, the more disease occurs) ( Figure 6A and Table S3). Meanwhile, AN, TN and disease had similar significant regulatory relationships as pH and disease ( Figure S5B,C and Tables S6 and S7). In addition, soil temperature affected the diversity of rhizosphere fungi, and had both direct and indirect effects on the abundance of Fusarium and Talaromyces.

Construction of Structural Equation Models Related to Pokkah Boeng Disease of Sugarcane
We constructed an SEM to assess the direct and indirect effects of soil properties (pH, AN and TN), microbial genera (Fusarium and Talaromyces), and microbial diversity (ACE and Shannon) on the incidence of sugarcane (Figures 6 and S5 and Table S2). We used a multigroup modelling method to assess which relationships exist between soil properties, microbial communities, and sugarcane disease during the growing process. The results of model-1 showed that soil pH was significantly and directly affected the abundance of Fusarium, sugarcane disease and fungal diversity. However, soil pH negatively regulated the disease rate of sugarcane (the lower the pH, the more disease occurs) ( Figure 6A and Table S3). Meanwhile, AN, TN and disease had similar significant regulatory relationships as pH and disease ( Figure S5B,C, Tables S6 and S7). In addition, soil temperature affected the diversity of rhizosphere fungi, and had both direct and indirect effects on the abundance of Fusarium and Talaromyces.

The Assembly Process of Fungal Communities Is Influenced by Changes in the Growth Period
The results showed that the mean nearest-taxon index (NTI) was greater than 0, and the mean nearest taxon distance between samples (β-NTI) values among the samples from the critical fertility period of sugarcane were mainly concentrated in the interval of −2 to 2, indicating that the changes in microbial community structure during the critical fertility period were mainly influenced by stochastic factors. The existence of a certain number of samples with β-NTI values greater than 2 among the soils in different periods indicated that the changes in the rhizosphere microbial community structure of some soil samples were mainly influenced by deterministic factors ( Figure 7B). In addition, the β-NTI values of samples between stages (Seeding-Tillering, Tillering-Elongation, and Elongation-Maturity) were between (−2, 2), indicating that the changes in the rhizosphere microbial community structure of sugarcane with time were mainly influenced by stochastic factors, but there were differences in the contribution of stochastic factors with time (Table S8).

The Assembly Process of Fungal Communities Is Influenced by Changes in the Growth Period
The results showed that the mean nearest-taxon index (NTI) was greater than 0, and the mean nearest taxon distance between samples (β-NTI) values among the samples from the critical fertility period of sugarcane were mainly concentrated in the interval of −2 to 2, indicating that the changes in microbial community structure during the critical fertility period were mainly influenced by stochastic factors. The existence of a certain number of samples with β-NTI values greater than 2 among the soils in different periods indicated that the changes in the rhizosphere microbial community structure of some soil samples were mainly influenced by deterministic factors ( Figure 7B). In addition, the β-NTI values of samples between stages (Seeding-Tillering, Tillering-Elongation, and Elongation-Maturity) were between (−2, 2), indicating that the changes in the rhizosphere microbial community structure of sugarcane with time were mainly influenced by stochastic factors, but there were differences in the contribution of stochastic factors with time (Table S8).

The Assembly Process of Fungal Communities Is Influenced by Changes in the Growth Period
The results showed that the mean nearest-taxon index (NTI) was greater than 0, and the mean nearest taxon distance between samples (β-NTI) values among the samples from the critical fertility period of sugarcane were mainly concentrated in the interval of −2 to 2, indicating that the changes in microbial community structure during the critical fertility period were mainly influenced by stochastic factors. The existence of a certain number of samples with β-NTI values greater than 2 among the soils in different periods indicated that the changes in the rhizosphere microbial community structure of some soil samples were mainly influenced by deterministic factors ( Figure 7B). In addition, the β-NTI values of samples between stages (Seeding-Tillering, Tillering-Elongation, and Elongation-Maturity) were between (−2, 2), indicating that the changes in the rhizosphere microbial community structure of sugarcane with time were mainly influenced by stochastic factors, but there were differences in the contribution of stochastic factors with time (Table S8).

Functional and Evolutionary Relationships among the Major Fungal Genera of Rhizosphere Soil
As for functional classification, the fungal communities in rhizosphere soil were classified by using the trophic mode (Figure 8). A maximum likelihood phylogenetic tree was conducted to further visualized the relationship among crucial fungal communities. Following the procedure, the top 28 genera, were classified into five guilds. They were divided into Undefined Saprotroph, Dung Saprotroph, Animal Pathogen, Plant Pathogen, and Endophyte, respectively. Based on the classification, the results revealed that 33.3% of these genera belonged to Undefined Saprotroph, 3.3 % belonged to Plant Pathogen, 6.6 % belonged to Animal Pathogen, 6% belonged to Dung Saprotroph, 3.3% belonged to Endophyte while 46.6% of these genera were unclassified. We also observed that these top 28 genera, were detected in 8 different phyla, of which five were identified as Ascomycota, Chytridiomycota, Glomeromycota, and Basidiomycota, respectively. Moreover, the plant pathogenic fungi observed in this study was Cochliobolus, and it is worth noting that the abundance of Fusarium genera was higher in Tillering and Maturity. In addition, according to the results of the FUNGuild classification difference ( Figure S6), in the comparison of the Seeding and Tillering of sugarcane growth ( Figure S6A), there were more classifications with significant differences (p < 0.05). They were Dung Saprotroph-Soil Saprotroph-Wood Saprotrop, Animal Pathogen-Endophyte-Lichen Parasite-Plant, Orchid Mycorrhizal-Plant Pathogen-Wood Saprotroph, Endophyte-Plant Pathogen, Animal Pathogen-Soil Saprotroph, respectively. However, with the shift in sugarcane stages, the number of significantly different FUNGuild classifications gradually decreased, accompanied by a decrease in the significance of differences between some FUNGuild classifications. (Figure S6B,C). Sci. 2023, 24, x FOR PEER REVIEW 13 of 22 Figure 8. Phylogram with fungal guilds of the top 28 fungal genera, color the branches according to the phylum level to which the species belongs; maximum likelihood tree for the sequences obtained through high-throughput sequencing. The relative abundance data were normalized using the Zscore method, and then the average of the relative abundance of each genus for each group of samples was plotted as a heat map, the bar chart showed the percentage of Reads for the species in different stages. Guild annotation used the FUNGuild database.

Changes in Sugarcane Growth Periods Alter the Composition and Diversity of Soil Properties and Fungal Communities
Soil microbial communities associated with plants can have strong influences on plant growth as well as contribute to soil health and sustainable production [25][26][27]. Therefore, understanding the temporal progression of rhizosphere microbiota is a prerequisite for plant and soil environmental improvement [28,29]. Previous research showed that plants' root microbiota composition varied with plant developmental stage [30,31], but these studies were carried out either with other model crops or under Figure 8. Phylogram with fungal guilds of the top 28 fungal genera, color the branches according to the phylum level to which the species belongs; maximum likelihood tree for the sequences obtained through high-throughput sequencing. The relative abundance data were normalized using the Z-score method, and then the average of the relative abundance of each genus for each group of samples was plotted as a heat map, the bar chart showed the percentage of Reads for the species in different stages. Guild annotation used the FUNGuild database.

Changes in Sugarcane Growth Periods Alter the Composition and Diversity of Soil Properties and Fungal Communities
Soil microbial communities associated with plants can have strong influences on plant growth as well as contribute to soil health and sustainable production [25][26][27]. Therefore, understanding the temporal progression of rhizosphere microbiota is a prerequisite for plant and soil environmental improvement [28,29]. Previous research showed that plants' root microbiota composition varied with plant developmental stage [30,31], but these studies were carried out either with other model crops or under greenhouse conditions. Our findings provide a detailed description of the rhizosphere fungal population during the entire sugarcane growth period in the field, as well as insights into how sugarcane growth and the soil environment influence the development of the rhizosphere fungal population. We showed temporal shifts in fungal community composition during the life of sugarcane, and these results are missing in sugarcane-related research in the field ( Figure 1A,B). Our findings showed that the nutrients in sugarcane rhizosphere soil exhibited significant changes during different growth periods (p < 0.05), which we assumed was precipitated by the difference in sugarcane requirements for different nutrients at different growth stages [32]. Moreover, although the alpha diversity of the sugarcane rhizosphere fungal community varies significantly between different growth periods, it has a gradually stable trend ( Figure S1A-D). It showed that sugarcane has the ability to regulate its own rhizosphere fungal environment to maintain a stable state under the natural growth conditions. Whether this is applicable to other crop systems requires more extensive and in-depth studies. Additionally, the species composition showed that rhizosphere fungi were mainly Ascomycota and Basidiomycota during the growth period of sugarcane ( Figure 3A), which was similar to the research data of Zeng et al. (2020) [33]. Stursova et al.
(2012) also found that compared with Basidiomycota, Ascomycota are more involved in cellulose decomposition [34]. Meanwhile, Chytridiomycota and Glomeromycota, which had a lower relative abundance, also showed regular and significant changes over time (p < 0.05). Determining the role, if any, of these low-abundance microbes in responding to changes in plant diseases, soil health, and so on will be a fascinating challenge for future studies. The Venn diagram showed that 451 OTUs were shared by sugarcane during the growth periods ( Figure 3B). The number of OTU unique to each of these periods varies considerably. To further explore, we used the source model of microorganisms to analyse the fungal transmission ratio between stages. Since more than 90% (91-95%, Figure 4B) of the fungal community were passed to the next period in each stage, it indicated that after the formation of the sugarcane rhizosphere fungal flora, although different microorganisms would be recruited or consumed at different periods, the overall structure would remain stable. Such a result also verifies the idea that plants can maintain resident soil fungal populations but not non-resident soil fungal populations, as has been previously verified in model plants [5].

Fungal Taxonomy, Function, and Soil Traits Were Interrelated, and the Transfer of Fungal Communities between Changes in Sugarcane Growth Periods Showed Regularity
Additionally, there was a strong correlation between soil environmental factors and sugarcane rhizosphere fungi populations in different growth stages. RDA and network map showed that TN, TK, pH, and AN were the main factors driving sugarcane fungal communities ( Figure 4A,D). These soil nutrients directly or indirectly affected the survival and growth of sugarcane fungal communities, which is in agreement with the results reported by Zhang et al. (2016). They indicated that organic matter, total N and total P significantly affected soil fungal community composition in the southeastern Tengger Desert [35]. To further explore the influence of these nutrient factors on the fungal community during the growth period of sugarcane, we used the mantel test to calculate the correlation between three matrices and further validated that TN, TK, and pH were still the main environmental factors driving the OTU composition of the fungal community ( Figure 4C). There were more genera that were decreasing in relative abundance in the rhizosphere over the life cycle of sugarcane, while fewer genera were increasing in relative abundance. It may imply the degree of sugarcane's control over the rhizosphere fungal population under natural conditions. Additionally, the stabilization of the fungal population makes it difficult to re-enrich once it is lost from the rhizosphere environment, and the difficulty of plant recruitment to soil fungi also may cause this phenomenon ( Figure 5D,E). This is similar to the results of tracking changes in the rice root microbial flora [28]. These data reinforce the separation role that the growth period plays in distinguishing plant rhizosphere microbial flora as observed in other studies [36][37][38], but it needs more research to illustrate. During the whole growth period of sugarcane, there are 10 fungal genera that have been decreasing, mainly including Ascomycota, Chytridiomycota, Eurotiales, Hypocreales, Sordariomycetes, and Tremellales. We speculate that the decline of these fungal genera is closely related to the genotype and developmental stage of sugarcane [39,40]. This requires us to conduct further verification. Such results provide a list of sugarcane rhizosphere fungi, which should now be targeted to elucidate their potential functions in the roots of sugarcane (symptomless colonizers vs. plant growth-promoters vs. pathogens) by targeted separation and sequencing technology. Furthermore, according to our results, Pseudallescheria had been enriched during the first three growth periods of sugarcane growth (Seeding-Elongation). However, Nectriaceae were significantly (p < 0.05) increased in the Tillering, Elongation and maturation stages ( Figure 5A-C). Studies have reported that Pseudallescheria is a filamentous pathogenic fungus. Moreover, it is not only a potential human and animal pathogen, but also exists in the soil environment [41,42]. Pseudallescheria is of special importance for biological health and whether its significant variation in the rhizosphere could lead to local and disseminated infection of the host is a matter of alarm [43]. Studying the time trend of this pathogenic fungus in the sugarcane rhizosphere has great significance for the defense against certain new diseases that occur in the early growth of sugarcane in the future. Similarly, the ascomycete family Nectriaceae also includes numerous plant and human pathogens [44]. The enrichment of Nectriaceae is presumed to be potentially related to certain diseases that occurred in the late stage of sugarcane growth. However, evidence of pathogenicity does not necessarily exist, which needs to be established through more in-depth and mechanistic studies. Nonetheless, our findings revealed the most basic information and provided the possibility of sugarcane disease research. Additionally, it was found in the co-occurrence network that some fungal interactions between genera disappeared and then reappeared with time, presumably because the "autonomous consciousness" of sugarcane and environmental factors adjusted the balance of fungal interactions. The specific regulatory pathways and the types and sources of metabolites involved in them are not known. Whether this phenomenon is associated with the development of pokkah boeng disease remains to be proven.

Structural Equation Modeling Demonstrates a Correlation between Key Fungi, Soil Properties, and Sugarcane Disease Rates
Structural equation models (SEMs) indicated that soil temperature is a critical factor affecting the α-diversity of rhizosphere fungal communities of sugarcane, which showed significance in all models that we constructed (Figure 6 and S5). Soil temperature changes triggered by seasonal changes had a certain perturbing effect on the rhizosphere community [45]. In addition, the strong direct negative effects of soil pH, AN and TN on the disease status of sugarcane suggest that soil infertile would increase the likelihood of sugarcane disease [46,47]. At the same time, Fusarium and Talaromyces also responded to changes in the soil nutrient environment, thus interacting adversely with the plant. The high-quality genome sequence of Fusarium andiyazi in China, published by Bao et al. (2021), provides more clarity on the role of Fusarium in sugarcane pokkah boeng disease [48].

Rhizosphere Fungal Community Assembly Is Mainly Influenced by Stochastic Factors
Furthermore, the results of the rhizosphere fungal community assembly process showed that the majority of the βNTI values among the samples were less than |2| ( Figure 7B), indicating that the changes observed in the sugarcane rhizosphere fungal community over time were primarily influenced by stochastic factors [49], but the stochastic contribution increased and then decreased with sugarcane growth and development, eventually reaching the lowest value at the maturity stage (Table S8). This is due to stochastic changes in the probability distribution and the relative abundance of species of rhizosphere fungi in the early stages of sugarcane growth (ecological drift) [50], and the formation of an adaptive system of root and soil environment in the later stages of sugarcane growth, resulting in abiotic and biotic factors playing an increasingly pivotal role in the presence or absence and relative abundance of rhizosphere fungi. Furthermore, mutual inhibition between rhizosphere fungi has an effect on community assembly, but the specific mechanism of action and whether it leads to a deterministic process of rhizosphere fungal communities needs to be supported by more research evidence. Phylogenetic trees and FUNGuild functional prediction analysis were applied to the study of sugarcane rhizosphere fungal communities to simplify and visualize the complex community functions and evolutionary relationships (Figure 8 and Figure S6). FUNGuild analysis revealed that the abundance of fungi showed differential responses across stages, for example, Fusarium was abundant in Seeding and Elongation and Talaromyces in Maturity. These results correlate with the work of Li et al. (2021), who reported that some potential biocontrol genera change with plant growth and tillage [51]. While showing dominant species richness, it also provided new ideas for seasonal prediction of rhizosphere pathogens of sugarcane and disease control during the different growth periods, such as annotated phytopathogens (Cochliobolus) that directly or indirectly contribute to the disease of sugarcane in a given stage.

Sugarcane Planting
In the spring of 2017, sugarcane cultivar, FN 41 was grown in a separate sugarcane experimental field in China to track the rhizosphere fungal population change procedure during the entire sugarcane growth cycle. Before planting, sugarcane seed-canes were soaked in 0.1% carbendazim (w/v) for 10 min to avoid surface-associated microbes followed by soaking in water for 24 h [52]. After soaking seed-canes, they were transferred to the fields at the Baisha experimental station (119 • 06 E, 26 • 23 N). Before the start of the experiment, the same variety of sugarcane crop was grown in the test plot for a period of three years. Abundant heat resources in the test area to supply sugarcane growth. The field growth of sugarcane was managed uniformly according to the farming conditions of local growers. A randomized group design was used, with an area of 1008 m 2 , 8 rows per zone, 5 m row length, 1.2 m row spacing, 48 m 2 per plot, and 3 replications, distributed in randomized groups in each plot. The planting height of sugarcane was 1/2 of the row height. Twenty-one sugarcane stalks were randomly selected in each field and measured with a measuring tape and vernier caliper to determine stalk height and stem diameter in each period. Meanwhile, the number of effective plants, diseased plants and total plants were also counted.

Sample Collection
The experiment started on the 3 January 2017. The soil samples were collected at the Baisha Experimental Station of the National Sugarcane Research Center of Fujian Agriculture and Forestry University (subtropical monsoon climate, the annual average temperature is 19.5 • C, the annual average precipitation is 1673.9 mm). At different fertility stages, we used the "S-shaped sampling method" to select 21 sampling points in the sugarcane field [53]. According to the shaking-off method of Riley and Barber [54], the soil adhered to the sugarcane roots was brushed with a small sterile brush, and soil samples were finally collected, sieved at 2.0 mm, and stored in a refrigerator at −20 • C [55]. The samples were collected during the sugarcane seeding stage (Seeding 4, 2017), tillering period (Tillering 2, 2017), elongation period (Elongation 17, 2017), and maturity period (Maturity 29, 2017), with a total of 84 samples (4 periods, 21 repetitions).

Determination of Pokkah Boeng Disease of Sugarcane and Soil Physio-Chemical Properties
The physicochemical indicators in the rhizosphere soil that had a strong influence on the sugarcane and varied greatly during the succession of growth periods were selected for measurement. Soil suspension with water (1:2.5 WV −1 ) was prepared to estimate soil pH using a pH meter (PHS-3C, INESA Scientific Instrument Co., Ltd., Shanghai, China) [56]. The soil temperature (Tem) was measured with the Soil Temperature Detector (Model: JC-TW, Shandong, China). The available nitrogen was measured using the alkaline hydrolyzable diffusion method [57], Soil total nitrogen and organic matter were determined by Kjeldahl digestion and determined by the oil bath-K 2 CrO 7 titration method [58,59]. C:N is the ratio of soil total N to organic matter. Soil total potassium and total phosphorus were determined by digestion with HF-HClO 4 , followed by flame photometry and molybdenumblue colorimetry, respectively [60,61]. Available potassium was extracted by ammonium acetate and determined by flame photometry [62]. Available phosphorus was extracted by sodium bicarbonate and determined using the molybdenum blue method [63]. Field judgment of pokkah boeng disease of sugarcane was mainly based on the symptoms described in previous studies [64] and divided into three types: the chlorotic phase, the acute phase or top-rot phase and the knife-cut phase (associate with top rot phase).
The effective number of stems/ha = Effective number of stems/m 2 × 666.67/average row spacing.

Data Quality Control and Filtering
The sequencing extracted data were saved in fastq format, while the paired reads were spliced (merged) into one sequence based on the overlap between PE reads, and quality control was performed on the read quality and merging effect. Valid sequences were obtained from the barcodes and primers at the first and last ends of the sequences, and then the double-ended sequences were spliced (Flash, 1.2.11) to generate an abundance table for each taxonomy (QIIME, 1.9.1). The complete data sets generated in our study have been deposited in the NCBI Sequence Read Archive database under BioProject ID PRJNA721464.

Sequences and Statistical Analysis
The UPARSE standard pipeline was used to analyze the sequence data [66]. Briefly, sequences with short reads (<250 bp) were filtered out before for downstream analysis. Sequences with ≥97% similarity were clustered into OTUs. All sequences were assigned using the RDP classifier to identify taxa with a confidence threshold of 0.8 [67]. We selected these OTU with 97% similarity, and then calculated the Alpha diversity index under different random sampling using Mothur [68]. Later, we used R to draw the rarefaction curves. The DPS software was used to analyze the variance of the soil physical and chemical properties, and the significance was calculated based on the Bonferroni test (FDR adjusted p < 0.05) [69]. RDA was used to visualize the relationship between fungal communities and soil environmental factors. Network analysis was performed using R to calculate the correlation between the factors (Spearman correlation), and Cytoscape (version 3.6.1) was used to adjust and visualize the results [70]. Analysis of differential OTU abundance and taxa was performed using a DESeq2 of the R package, and then we used a Manhattan plot to visualize the results (R 3.6.0). Structural equation modeling (SEM) was performed using IBM SPSS Amos 26. Maximum likelihood estimation with standard errors was also used [71]. The nearest-taxon index (NTI) and βNTI (999 random) were used to quantify changes in rhizosphere fungal phylogeny over time, and the two indices were calculated using the package "picante" [72,73]. The heatmap and functional annotation and maximum likelihood trees were created using the majorbio platform (http://cloud.majorbio.com, accessed on 12 September 2022). FUNGuild classification map was completed using STAMP (version 2.1.1), comparing the two stages, using Welch's t-test, 95% confidence intervals, p < 0.05. Network diagramming and parameter were calculated using R, Cytoscape and UCINET 6 together [74].

Conclusions
When focusing on soil productivity, rhizosphere microbiomes are a critical factor and closely related to plant health, especially the fungal community. Under the background conditions of this study, both the growth stages and soil characteristics of sugarcane had significant effects on the composition and function of the rhizosphere fungal community. The pH, TN, TK and AN decreased significantly during the growth of sugarcane. TN, TK, and AN were significantly decreased by 12.41%, 26.11%, and 31.32%, respectively, compared to the early growth stage. The dominant fungal phyla were Ascomycota, Basidiomycota, and Chytridiomycota, and the dominance of dominant fungal groups can be explained by both the temporal dynamics of sugarcane roots and the awareness of plant-autonomous regulation. We also discovered that diseases had significant and strong negative effects on selected soil traits (pH, Tem, and AN) using SEM, while soil temperature had a strong and direct positive effect on fungal α-diversity, and Fusarium and Talaromyces had some direct and indirect associations with sugarcane diseases. It indicates soil nutrients affected plant health and rhizosphere fungal interactions. The variation in rhizosphere fungal structure was mainly influenced by stochastic factors, and the decrease in stochastic contribution rate was hypothesized to be due to the coupling of sensitive fungal genera or modules in the rhizosphere fungal network with sugarcane growth dynamics to coordinate the growth process. These results contribute to the understanding of plant-rhizosphere fungal interactions and provide additional opportunities for the development of green agriculture and preventive tools for crop pests and diseases.