Abstract
Introduction:
Endophytic fungi exhibit diverse interactions with plants, from pathogenic to mutualistic symbiosis, and the community composition is regulated by phytohormones. Yet, the composition and dynamics of endophytic fungi in Paris polyphylla var. chinensis (Franch.) Hara (PPC) during fresh seed (FD), sand-stored seed (SSD), and seedling (SS) stages remain unclear. Similarly, the overall impact of phytohormones on the management of endophytic fungal communities is yet to be elucidated.
Methods:
We carried out a pot experiment to examine the effects of various stages of PPC seeds and the external addition of three phytohormones, namely, melatonin (MT), strigolactone (SL), and 24-epibrassinolide (BR) on the endophytic fungi of PPC seedlings. This was done through internal transcribed spacer (ITS) amplicon sequencing.
Results:
The study of the endophytic fungal microbiome in FD, SSD, and SS stages of PPC revealed an increased richness and diversity of fungi during the SS stage, with significant changes in community composition observed. We found that Sordariomycetes played a crucial role in this process, potentially contributing to the establishment and growth of PPC seedlings. Additionally, this study investigated the influence of phytohormones on the phenotypic and physiological characteristics of PPC and its endophytic fungal community. Our results demonstrated that MT and SL significantly increased PPC biomass by 69.32 and 15.23%, respectively, while 2 mg/L of BR hindered the growth of PPC roots. MT, SL, and BR not only induced significant changes in the composition and diversity of the endophytic fungal community in PPC but also affected biomass potentially through specific regulation of potential biomarkers. Furthermore, phytohormones were shown to indirectly modify the endophytic fungal community by altering antioxidant system in plants.
Conclusion:
This study provides novel insights into the dynamic changes of microbial communities in the FD, SSD, and SS stages. Furthermore, the differences among various phytohormones ultimately enhance our predictive understanding of how to directly or indirectly manipulate the plant microbiome to improve plant health.
1 Introduction
Paris polyphylla var. chinensis (Franch.) Hara (PPC) is a perennial flowering herb that belongs to the Melanthiaceae family. Formerly, it was classified under the Liliaceae family and is endemic to China. This species is listed in the International Union for Conservation of Nature’s Red List (IUCN) (). Additionally, the 2020 edition of Chinese Pharmacopeia includes PPC as one of the two original medicinal plants utilized in Paridis Rhizoma (). Modern pharmacological studies have demonstrated that the primary active component, steroid saponins, exhibits significant biological activities such as anti-tumor, anti-inflammatory, and hemostatic properties. Its efficacy and scope of application continue to expand (). Research on the Chinese pharmaceutical market indicates that PPC is a key ingredient in over 80 traditional Chinese medicine formulations, with annual sales ranging from 800 to 1,050 tons. The annual consumption far exceeds its natural growth capacity, posing a threat of depletion to its wild medicinal resources ().
Endophytic microorganisms in plants primarily originate from two sources: the external environment on the plant surface, and seeds (). Endophytic microorganisms, including bacteria, fungi, and actinomycetes, are transmitted vertically between successive plant generations (). The relationship between endophytic fungi and plants is extensive, ranging from being potential pathogens or saprophytes to mutualistic symbiosis (). Consequently, this intergenerational transmission significantly impacts plant health, quality, productivity, and microecology (; ). In agricultural production, due to the lack of cost-effective tissue culture techniques, PPC can only be propagated via seeds. Fresh PPC seeds are immature, both morphologically and physiologically. Following 80 days of warm wet sand stratification before sowing, the PPC seeds complete their entire process of embryo morphological maturation, from nearly spherical to elliptical to cylindrical embryos. After an additional 60 days of low-temperature stratification, the seeds reach physiological maturation. During the sand stratification process, the seeds may be invaded and colonized by external fungi (). After reaching morphological and physiological maturity, changes in the internal nutritional and hormonal physiological levels of seeds alter the habitat and promote the growth of certain microorganisms within the seeds (). Following seed germination, the seed coat gradually sheds, triggering the reassembly of microbial communities on the plant surface and within its tissues. These communities progressively converge toward a similar composition, with the presence of endophytic fungi serving as a key indicator of the completion of this assembly process (). However, compared to the seed stage, this phase typically undergoes significant compositional changes. Studies have shown that the endophytic bacterial communities in seeds vary not only due to the seed development process and environmental conditions (), but also exhibit significant differences among different plant species (). Due to their positional advantage, endophytic fungi within the seeds may influence plant growth and adaptability, from seeds to seedling emergence, and continue to affect plant development over time. Therefore, understanding the dynamics of fungi during the processes of fresh seeds, seed storage, and early seedling germination is crucial for the selection and maintenance of endogenous microorganisms beneficial to plant growth and health, as well as monitoring changes in pathogenic bacteria.
The radicle begins to develop after seed germination, leading to seedling growth. PPC typically grows in valleys, beneath forests, or amid shrubbery, where plants compete for sunlight, air, and water. Relative to the fragility of seedlings during the early stages of growth, their swift and robust development influences subsequent growth and yield (). Therefore, the early growth of seedlings is crucial because it enhances the plant’s ability to compete with other species in the field. Phytohormones have long been used to address sluggish early growth of plant seedlings. Research shows that phytohormones play a broad role in plant physiological processes, regulating plant growth and development, and they are critical in shaping the plant microbiome (). These hormones can act as nutrients and signaling molecules to regulate the growth and metabolism of endophytic microorganisms (), or selectively restrict the growth of certain members within the plant microbiome by modulating the reactive oxygen species (ROS) system (). Brassinosteroids (BR) are a type of efficient, broad-spectrum steroidal plant hormone involved in processes like cell elongation, tissue differentiation, and responses to abiotic stress (). Research confirms that BRs promote plant growth by enhancing antioxidant capacity and photosynthetic efficiency (). Moreover, BRs reportedly increase rhizosphere fungal diversity (), and positively impact arbuscular mycorrhizal (AM) symbiosis (). Strigolactone (SL), a novel carotenoid-derived phytohormones is extensively involved in plant growth, including seed germination, seedling development, bud branching, stem elongation, and lateral root formation (). In terms of interactions with microorganisms, SL, similar to BR, has also been reported to participate in AM symbiosis. It reportedly influences significant differences in fungal community composition, but without significant differences in bacterial community composition (). SL can enhance plant resistance to pathogenic fungi (). Melatonin (N-acetyl-5-methoxytryptamine) is a natural indoleamine compound that aids in ROS scavenging, photosynthetic system regulation, leaf senescence delay, and tolerance promotion to some biotic stresses by reducing the abundance of several harmful fungi in the rhizosphere soil (; ). While previous studies successfully manipulated specific endophytic fungi in plants through the exogenous application of synthetic plant hormones, there have been no reports on manipulating endophytic fungi at the overall community level using phytohormones.
In this study, we delineated the endophytic fungal communities in PPC throughout three stages: fresh seeds (FD), sand-stored seed (SSD), and seedling (SS), and in response to the treatment of three phytohormones (MT, SL, and BR). The study objectives were as follows: (1) to examine the dynamic changes in endophytic fungi within PPC seeds from the fresh to the stored and seedling stages; (2) to scrutinize the specific impacts of the exogenous addition of three phytohormones on these endophytic fungi; and (3) to further analyze the correlation between the PPC endophytic fungal microbiota and certain physiological parameters, to understand the influence of these phytohormones on the assembly of PPC endophytic fungal microbiota.
2 Materials and methods
2.1 Site description
The study was carried out at Biological Research Institute of Gansu Academy of Sciences and Heping testing base, Longnan (36.008 N, 103.970 E), Gansu, China. The experimental site is in the eastern monsoon region, and has a cold, semi-humid, and rainy climate, with an altitude of 1,700 m. The growing season of PPC seedlings is from December 2018 to October 2019. The soil is classified as Calcaric Cambisol according to the Food and Agricultural Organization (FAO) classification system. The initial physicochemical characteristics of the soil were: pH = 7.84, organic matter = 1.19 g/kg, TN = 0.63 g/kg and available K = 3.42 mg/kg. During the process of collecting the soil, large particles of impurities, including plant litter, plant roots, and gravel, were removed. The soil was thoroughly mixed before being filled into the flower pots.
2.2 Plant materials and experimental design
In this study, PPC samples were collected at three stages: fresh seeds (FD), sand-stored seeds (SSD), and the seedling (SS) stage. The fresh fruits (Figure 1A) were collected from Guanshan, Huating City, Gansu Province, China, at an elevation of 2,000 m (35°200 N, 106°396 E). To obtain fresh seeds, the fruits were carefully processed by removing the pulp and peel (Figure 1B). Before sowing, fresh seeds require a two-step temperature-alternating sand storage protocol to break dormancy. Initially, seeds undergo a warm stratification phase of approximately 40 days at 18°C, during which PPC seeds complete the full process of embryonic morphological maturation, transitioning from nearly spherical to elliptical, and finally to a cylindrical embryo. After this, seeds are subjected to a cold stratification period of 60 days at 4°C, at which point the seeds attain physiological maturity. Seeds collected at this stage are used as the post-stratification sample (Figure 1C). The sand-stored seeds are then transplanted into plastic pots (17.6 cm × 26.5 cm × 23 cm). The pre-germination treatment of seeds is conducted under moist and cool conditions (with a 70% shade net) to ensure rapid and uniform germination following sowing. Post-germination, the seedlings are grown under natural temperature conditions, and seedling roots at the thirteenth week are collected as samples for the seedling stage (Figure 1D).
FIGURE 1
The present study implemented four treatments, which comprised individual inoculations with MT (MT), SL (SL), BR (BR), and an equivalent volume of sterile distilled water as the control (CK). The treatments were arranged in a randomized block design with three replicates per treatment, containing three seedlings per pot. The concentrations of MT, SL, and BR used in the experiment were in line with previous research by , all at 2 mg/mL and 15 mL per pot. Starting from the 8th week of the seedling stage, applications were made once a week for five consecutive weeks. Seedling roots were collected in the 13th week to be used as samples for the seedling stage evaluation.
All collected samples underwent surface sterilization. Initially, the samples were rinsed with sterile distilled water at least three times, or until no turbidity was observed. Subsequently, the samples were immersed in a solution of 5% NaClO for 10 min, followed by three additional rinses with sterile distilled water. After this, 150 mL of the final rinse water was plated onto potato dextrose agar plates. The plates were then incubated at 28°C for 72 h to confirm the effectiveness of the seed surface sterilization. The samples were promptly stored at −80°C for subsequent utilization. The experiment included six treatments, named FD, SSD, SS (CK), MT, GR, and BR. Each treatment had three replicates, resulting in a total of 18 samples.
2.3 Measurement of phenotypic, physiological, and hormone content of PPC
After phytohormone treatment, we randomly selected three pots per treatment group to measure root length, shoot length, and root fresh weight. Plant samples were placed in an oven at 105°C for 30 min to inactivate the enzymes. Next, they were dried at 70°C until they reached a constant weight; we then recorded the biomass. We calculated all indicators based on an average of three plants per pot. We transferred 1.0 g of mixed plant samples to a plastic centrifuge tube containing 18 mL of precooled 10 mM phosphate-buffered saline (130 mM NaCl, 7 mM Na2HPO4, and 3 mM NaH2PO4; pH 7.4). The plants were then homogenized using a homogenizer. Afterward, we centrifuged the homogenized sample at 3,000 rpm/min for 20 min at 4°C. We then collected and stored the supernatant at 4°C for physiological index measurement. In this study, we determined four physiological indices: malondialdehyde (MDA) content, hydrogen peroxide (H2O2), superoxide dismutase (SOD) activity, and catalase (CAT) activity. We measured all indices using plant Enzyme-Linked Immunosorbent Assay (ELISA) kits (Mibio, Shanghai, China) following the kit instructions. We purchased test kits for MDA (Lot. No.20200521), H2O2 (Lot. No.20200526), SOD (Lot. No.20200619), and CAT (Lot. No.20200603) from Suzhou Keming Biotechnology Co., Ltd. Net photosynthetic rate (Pn), stomatal conductance (Gs), inter-cellular CO2 (Ci), and transpiration rate (Tr) were determined by a portable infra-red gas analyzer—photosynthesis system, LI-6400XT (Lincoln, NE, United States) between 10 a.m. and 12 noon. The contents of Salicylic acid (SA), MT, SL, and BR in plant were measured using plant ELISA kits (Mibio, Shanghai, China) according to the provided instructions.
2.4 Soil microbial DNA extraction and sequencing
The total DNA of each sample was extracted using a DNeasy PowerSoil DNA Extraction Kit, following the manufacturer’s instructions (QIAGEN). The identification region for the analysis of fungal diversity, ITS I variable regions were amplified with universal primers (primers: ITS1F 5′-CTTGGTCATTTAGAGGAAGTAA-3′; ITS2 5′-GCTGCGTTCTTCATCGATGC-3′). Sequencing was conducted on an Illumina NovaSeq 6,000 instrument post PCR amplification and purification. Sequencing reads were processed to obtain effective sequences using barcodes for all samples. Then, alignment, filtration, and removal of chimeras from reads were accomplished using software applications such as Pandaseq, Prinseq, and Usearch. Sequences were clustered using Usearch at a similarity threshold of 0.97, chimeras were filtered from the resultant cluster sequences, and OTUs (operational taxonomic units) were collected for species classification. Then, sequences with a similarity threshold of 0.97 were aligned against the representative sequences of the OTU for additional analysis. Taxonomic data annotation was done for representative sequences of fungi using UNITE (version 8.0)1 databases, respectively. The microbial community structure was investigated using the Illumina NovaSeq 6,000 platform (Nanjing GENEPIONEER Biotechnology Co., Ltd., China) for high-throughput sequencing of fungal ITS rRNA.
2.5 Statistical analysis
Alpha diversity indices, including community richness (Chao1) and community diversity (Shannon), were utilized to indicate dissimilarity among samples. The beta diversities across the different stages were determined by principal coordinates analysis (PCoA), and the significant differences were determined by permutational multivariate analysis of variance (PERMANOVA). Non-metric multidimensional scaling (NMDS) and multivariate PERMANOVA were employed to reveal the succession of endophytic fungi, as well as to measure differences in community structure among different treatments. Co-occurrence network analysis and visualization were completed using Gephi software (v.0.9.2), and only significant relationships were included (R > 0.85, p < 0.05). Four node-level topological features were calculated per node: degree, betweenness, closeness, and eigenvector centrality. To investigate relationships between phytohormone contents, physiological properties, and endophytic fungal communities on PPC biomass, partial least squares path modeling (PLS-PM) was performed using the “plspm” package in R.3.1.3. Random forest analysis was performed using the “randomForest” package in R.3.1.3. Significant pearson correlation (p < 0.05) between phytohormone and microbial taxa were analyzed and visualized using the interactive networks in R.3.1.3, Cytoscape 3.3.0. Only significant relationships were included (R > 0.5, p < 0.05) in co-occurrence network analysis.
The data was analyzed using Analysis of Variance (ANOVA) in SPSS 17.0 for Windows (IBM SPSS Inc., Chicago, United States), and the results are expressed as the mean ± standard deviation. A p-value less than 0.05 suggested significant differences, while a p-value less than 0.01 indicated extremely significant differences.
3 Results
3.1 Diversity and endophytic fungi communities composition of different growth stages
In this study, the ITS gene was sequenced to identify endophytic fungi at three key growth stages of PPC: FS, SSD, and SS. A total of 831,954 high-quality sequences were generated from nine samples. These were classified into 209 fungal OTUs at a 97% clustering level, with 131, 141, and 185 OTUs in FS, SSD, and SS, respectively. Alpha diversity showed that the fungal richness (Chao1) and diversity (Shannon) during the SS stage were significantly higher than those of FS and SSD, with no significant difference between FS and SSD (Figure 2A). Furthermore, the PCoA, based on the Bray-Curtis metric (Figure 2B), revealed substantial differences in fungal composition between SS and the other two stages. Specifically, FS and SSD samples partially overlapped after clustering, while SS samples were separate from the other two stages on the first coordinate axis. The overall network of endophytic fungal communities at different stages portrayed a clear symbiotic pattern. The fungal communities of FS, SSD, and SS had meta-network node numbers of 130, 142, and 185, respectively, with significant edge numbers of 949, 951, and 1,700, respectively (p < 0.05) (Figure 2C). Generally, almost all links between OTUs exhibited positive correlations. Further comparison of the specific node topological characteristics showed that the betweenness centrality and closeness centrality values of SS were notably greater than those of FS and SSD (Figure 2D). Overall, FS and SSD endophytic fungi were significantly different from the SS stage. The dominant phyla of endophytic fungal communities under different treatments were Sordariomycetes (79.048%), Saccharomycetes (5.948%), Leotiomycetes (0.118%), Eurotiomycetes (0.116%), Mortierellomycetes (0.091%), Dothideomycetes (0.024%), Agaricomycetes (0.017%), Chytridiomycetes (0.002%), Glomeromycetes (0.002%), Spizellomycetes, and Ustilaginomycetes (Figure 2E).
FIGURE 2
3.2 Discrimination of endophytic fungi at different stages of PPC seeds
We used random forest analysis to identify the endophytic fungal taxa distinguishing seeds at different stages (Figure 3A). Ranked by their importance value, the top 20 fungal OTUs mainly belonged to Ascomycota, specifically Stachybotryaceae, Claussenomyces, Nectriaceae, Sordariales, Fusarium, and Trichoderma. The relative abundance results showed that the four pathogenic fungi: Stachybotryaceae, Claussenomyces, Sordariales, and Fusarium, could be viewed as key taxa in the SS stage with a corresponding significant increase in the biocontrol fungus Trichoderma. The pathogen Nectriaceae could be considered a key taxon in the FS stage, showing a significantly higher abundance compared to the other two stages (Figure 3B).
FIGURE 3
3.3 Phenotypic, physiological characteristics and hormone content of PPC
We observed significant changes in the phenotype of PPC treated with sterile water, MT, SL, and BR (Table 1). In general, compared to CK, plants treated with MT depicted significantly enhanced robustness, with respective increases of 39.45, 10.34, 67.11, and 42.48% in root length, shoot length, fresh weight, and biomass (p < 0.05). SL did not affect the primary root growth significantly, but it led to a significant increase in shoot length, fresh weight, and biomass. BR notably improved shoot growth and exerted a distinct inhibitory effect on primary root growth, and did not significantly influence biomass accumulation. The activities of SOD and CAT were 32.81–108.48% and 20.63–42.86% higher in the three hormone-treated groups than the control, respectively (Table 1). Plants treated with MT exhibited low levels of MDA and H2O2. In contrast, SL and BR treatments significantly raised H2O2 levels, suggesting that the plants under these treatments might be experiencing oxidative stress. MT, SL and BR improved the Pn by10.74–40.00% than CK. Gs and Tr were 14.08–49.32% and 12.88–15.02% higher in MT and SL treatments, than CK, respectively, while Ci was significantly reduced. Treatment with BR significantly increased the gs and Ci transpiration rates compared to CK, while TR significantly decreased. With respect to the control, the phytohormone contents in the MT, SL, and BR treatments were 0.51, 0.69, and 0.34 mg/kg, respectively, marking an increase of 66.66, 56.25, and 66.52%, respectively. The content of SA significantly increased under the influence of all phytohormone treatments.
TABLE 1
| Plant index | Treatments | |||
| CK | MT | SL | BR | |
| Root length (cm) | 3.63 ± 0.15b | 5.07 ± 0.15a | 3.47 ± 0.06bc | 3.3 ± 0.1c |
| Shoot length (cm) | 0.97 ± 0.06c | 1.07 ± 0.06b | 1.21 ± 0.02a | 1.07 ± 0.02b |
| Fresh weight (g) | 7.50 ± 0.20c | 12.53 ± 0.40a | 9.50 ± 0.10b | 7.23 ± 0.25c |
| Biomass (g) | 4.40 ± 0.10c | 7.45 ± 0.15a | 5.07 ± 0.40b | 4.87 ± 0.25bc |
| SOD activity (U/g) | 20.05 ± 1.51d | 31.96 ± 0.79b | 41.80 ± 2.54a | 26.63 ± 1.19c |
| CAT activity (U/g) | 854.28 ± 48.89c | 1211.36 ± 10.36a | 1030.56 ± 52.95b | 1220.4 ± 17.94a |
| H2O2 content (μmol/g) | 473.86 ± 4.25c | 421.75 ± 19.78c | 612.28 ± 1.32a | 539.12 ± 58.68b |
| MDA content (nmol/g) | 66.56 ± 2.68a | 32.25 ± 0.93d | 55.38 ± 1.42b | 38.36 ± 4.73c |
| Pn (μmol⋅CO2⋅m–2⋅s–1) | 3.92 ± 0.32c | 5.49 ± 0.24a | 4.62 ± 0.37b | 4.35 ± 0.31bc |
| Gs (mmol⋅CO2⋅m–2⋅s–1) | 0.031 ± 0.00c | 0.03 ± 0.00b | 0.05 ± 0.00a | 0.05 ± 0.00a |
| Ci (μmol CO2 mol–1) | 48.96 ± 2.85b | 36.86 ± 0.95d | 44.75 ± 3.44c | 53.16 ± 2.06a |
| Tr (mmol⋅H2O⋅m–2⋅s–1) | 1.96 ± 0.10b | 2.21 ± 0.02a | 2.26 ± 0.04a | 1.78 ± 0.08c |
| Salicylic acid content (mg/kg) | 15.39 ± 0.31c | 26.16 ± 0.89a | 19.30 ± 3.52b | 22.83 ± 0.37a |
| Melatonin content (mg/kg) | 0.31 ± 0.05b | 0.51 ± 0.06a | 0.21 ± 0.03b | 0.21 ± 0.06b |
| Strigolactone content (mg/kg) | 0.44 ± 0.04b | 0.49 ± 0.03b | 0.69 ± 0.14a | 0.42 ± 0.03b |
| 24-epibrassinolide content (mg/kg) | 0.21 ± 0.02c | 0.24 ± 0.01bc | 0.29 ± 0.04ab | 0.34 ± 0.03a |
Phenotype, physiological properties and phytohormone content of Paris polyphylla var. chinensis (Franch.) Hara (PPC) treated with different phytohormones.
Data are the means of three replicates with standard deviations (±SD). Different letters (a, b, c, and d) indicate significant differences among treatments according to Dunn’s test (p < 0.05), while the same letter signifies no significant difference. CK, sterile water; MT, Melatonin; SL, strigolactone; BR, 24-epibrassinolide.
3.4 Variations in the diversity of rhizosphere soil microbial communities of PPC following phytohormone treatments
We found that the Chao1 and Shannon indices of endophytic fungi in PPC under MT and SL treatments were significantly lower in comparison to CK (Figure 4A). We further analyzed the variations in soil microbial community structure under phytohormone treatments using NMDS analysis (Figure 4B). The NMDS plot revealed a statistically significant separation (p < 0.05) between the microbial communities of the CK treatment and those associated with MT, SL, and BR treatments, indicating that phytohormone application significantly altered microbial community composition.
FIGURE 4
At the class level, Sordariomycetes, Eurotiomycetes, Mortierellomycetes, and Leotiomycetes were dominant in all samples (Figure 4C). Based on pertinent findings, we derived the co-occurrence network of the PPC endophytic fungal community across CK and various hormone treatments (Figure 4C). CK had 290 edges and 76 nodes, MT had 128 edges and 65 nodes, SL had 210 edges and 68 nodes, and BR had 240 edges and 67 nodes. MT demonstrated the fewest nodes and edges. We studied unique community-level topological traits at the node level, chiefly focusing on degree, betweenness centrality, closeness centrality, and eigenvector centrality (Figure 4D). The degree, betweenness centrality, and eigenvector centrality of BR were notably higher than those of CK, suggesting a more pronounced level of symbiotic connection.
3.5 Potential effects of phytohormone, antioxidant system, photosynthetic rate, and fungal communities on PPC biomass
To examine the potential direct and indirect effects of phytohormone, antioxidant systems, photosynthetic rate, and fungal communities on PPC biomass, we conducted a PLS-PM analysis based on the effects and relationships of known predictors (Figure 5). The final model fit our dataset well and suggested that the MT content could be the most significant factor directly influencing the PPC biomass. MT, SL, and BR appeared to directly influence the photosynthetic rate. Additionally, MT and BR might modulate the diversity of endophytic fungi indirectly through their influence on antioxidant enzyme activities. Furthermore, BR exerts a direct effect on the diversity of endophytic fungi. Interestingly, neither the α-diversity of endophytic fungi (including Chao1 and Shannon) nor the structure of their community (NMDS1) yielded a direct and significant influence on biomass.
FIGURE 5
3.6 Discrimination of PPC endophytic fungi at different phytohormone treatments
We applied the random forest model to regress the relative abundance of root endophytic fungi at the genus level in PPC against different phytohormone treatments, determining which endophytic fungal taxa set apart from the different hormone treatments (Figure 6A). The model depicted the top 20 root microbiota with a variation index associated with different hormones. Figure 6B illustrates the specific shifts in significant fungal taxa under different hormone treatments, with MT augmenting the abundance of Funneliformis. SL notably increased the abundance of Funneliformis and significantly intensified the abundance of the pathogenic fungus Xylaria. BR markedly increased the abundance of Cryptosphaeria. A linear regression analysis revealed that among the key genera, Cryptosphaeria and Xylaria significantly negatively influenced biomass.
FIGURE 6
Pearson correlation analysis showed that indicated that endophytic fungi had a stronger association with phytohormone and antioxidant systems than with photosynthetic systems (Figure 6C). The abundances of fungal genera Arnium, Claussenomyces, and Mortierella were found to negatively correlate with PPC BR content. The abundance of genera Spizellomyces, Claussenomyces, and Gibellulopsis showed a positive correlation with plant MDA content. SL content was found to have a significant negative correlation with Xylaria. Furthermore, there was a significant positive correlation between Pn and Aspergillus, as well as a significant negative correlation between Tr and Cryptosphaeria.
Owing to the notable changes within the endophytic fungal community in PPC spurred by phytohormones, we delved deeper into their potentially interlinked consequences on the increase in PPC biomass (Figure 6D). Our findings point out that among the key taxa significantly influencing alterations in the endophytic community under specified treatments with various phytohormones, species notably impacting biomass include Cryptosphaeria (R = −0.63, p = 0.028), Xylaria (R = −0.78, p = 0.002), Simplicillium (R = −0.59, p = 0.043), Humicola (R = −0.59, p = 0.045), Gibellulopsis (R = −0.7, p = 0.011), and Fusarium (R = −0.61, p = 0.035). Each of these revealed significant negative correlations with PPC biomass.
4 Discussion
4.1 Endophytic fungal dynamics in fresh seeds, sand-stored seeds, and seedlings
The study results indicate that there is no significant change in endophytic fungal diversity after PPC seeds are stored in sand. This can be attributed to the protective function of the seed coat, the outermost layer of PPC seeds, which shields the seeds from mechanical damage and prevents the invasion of diseases and pests (). Notably, Fusarium and Plectosphaerella are significantly enriched as potential biomarkers during this stage, with reports suggesting that these fungi are reputed for causing crop root rot and wilt through widespread colonization (; ).
Paris polyphylla var. chinensis (Franch.) Hara seedlings significantly alter the diversity and community structure of their endophytic fungi. It is posited that this shift is associated with the shedding of the seed coat, which initiates endophytic fungal interaction with the surrounding environmental microbiota. In this study, the increased relative abundance of symbiotic microorganisms such as Aspergillaceae, Coniochaetales, and Trichoderma during seed germination is ascribed to the efficient nutrient supply from the endosperm formed during germination, fostering the growth of the microbial communities within the seed. These microbial communities, possibly in a dormant state inside the seed, utilize the nutrients released during germination to break dormancy (; ).
The phylum Ascomycota was the most abundantly annotated group of endophytic fungi at various stages of PPC seed germination. In our study, the early seedling growth stage led to the enrichment of multiple Ascomycota biomarkers, including Gibellulopsis and Cryptosphaeria, two potential plant pathogens that may detriment growth (; ), as well as Sordariomycetes, a prevalent category with anti-pathogenic activities in plants (). Furthermore, we found that the topological values of the Sordariomycetes taxa were significantly higher than those of other groups, indicating that Sordariomycetes often occupy central positions in the network. Therefore, we presuppose that Sordariomycetes plays a significant role during the PPC seedling growth stage.
4.2 Responses of plant phenotype and physiology to phytohormones addition
The pot experiment demonstrated that both MT and SL significantly promoted PPC biomass accumulation. On the one hand, it has been shown that three phytohormones enhanced the Pn rate to varying degrees, providing plants with more carbon sources for synthesis and positively affecting plant growth (). On the other hand, all three phytohormones can regulate the production and elimination of reactive oxygen, which in turn affects various aspects of plant growth and development (). The three hormones notably increased the levels of SOD and CAT. MT exhibited lower levels of MDA and H2O2, implying that the oxidation-reduction state of the plant is relatively stable. Numerous reports have suggested that MT primarily functions to protect plants from various abiotic and biotic stresses by scavenging ROS and regulating antioxidant levels (). Moreover, studies have found that MT does not significantly impact the membrane lipid peroxidation system (MDA and H2O2) of the plant under normal growth conditions but can decrease their MDA and H2O2 content under oxidative stress conditions (; ). Based on this evidence, we posit that the natural slow growth of PPC seedlings may be due to oxidative stress in their environment. Therefore, additional supplementation of phytohormones for scavenging reactive oxygen species is deemed necessary. The inhibition of root growth in the BR treatment may be due to the relatively high concentration chosen for this study, which may not be optimal for PPC plants (). Under SL treatment, a significant increase in H2O2 content and root length was observed, while shoot growth declined. The study revealed that the SL-induced accumulation of H2O2 initiates various reactions such as transcription and oxidative modification, ultimately stimulating root development (). The restriction of shoot growth aligns with the conclusions drawn by . Encouragingly, MT exerted a stronger effect than SL and BR on nearly all growth parameters. To some extent, this supports the role of melatonin as a primary regulator in PPC, in agreement with previous research findings ().
4.3 Responses of plant endophytic fungal microbial communities and key taxa to phytohormones addition
Manipulating the endophytic microbial communities in natural environments is critical for leveraging ecological interactions to enhance plant productivity. Therefore, it is essential to conduct pot experiments in uncontrolled, natural conditions, to examine the impact of phytohormones on the endophytic fungal community of PPC. In this study, all three phytohormones significantly altered the structure of the endophytic fungal community. We used PLS-PM to assess the relationships between hormone concentrations, physiological parameters (such as antioxidant enzyme activities and photosynthetic rates), endophytic fungal community composition (diversity and structure), and plant biomass. In the current study, MT and BR indirectly affected endophytic fungal diversity by enhancing the activity of antioxidant enzymes. Fungi like Gibellulopsis, Spizellomyces, and Claussenomyces showed significant positive correlations with MDA content, possibly due to the hypersensitive response triggered by plant infection and fungal recognition, which leads to increased lipid peroxidation of the cell membrane (). SOD and CAT are vital protective enzymes involved in scavenging ROS in plants. The coordination of ROS levels with the abundance and diversity of microbial populations is well-documented (). Xylaria and Humicola demonstrated significant positive correlations with SOD activity, which may be due to the tendency of commensal fungal genomes to enrich genes for ROS-scavenging enzymes (). MT content had a direct positive effect on biomass, but the diversity and overall structure of the endophytic fungal community did not vastly influence biomass. We conducted a subsequent random forest analysis to identify potential biomarkers of the endophytic fungal community responding to various phytohormone treatments. These key taxa are generally considered critical factors in differentiating microbial community structure and they exert varied effects on plant growth. In this study, MT exerted a significant positive impact on the biomass of PPC by markedly inhibiting the abundance of the pathogenic fungus Fusarium. The study showed the significant enrichment in two key taxa, Funneliformis and Xylaria, in that SL. Xylaria is identified as a pathogen that negatively influences crop root development (). The results of linear analysis also indicate a negative correlation between Xylaria and biomass. In the BR treatment, Cryptosphaeria was significantly enriched, and reports have demonstrated that this pathogen can cause cankers in crops and forest plants (). Consequently, Cryptosphaeria could be a key taxon responsible for the adverse effects of BR on PPC growth. In conclusion, the effects of various plant hormones on the endophytic fungi within PPC are specific and complex. MT, SL, and BR can regulate the antioxidant system, thereby influencing diversity and structural alterations within the endophytic fungal community. The findings of this study underline the view that phytohormones indirectly impact microbial communities by modifying plant physiological properties (; ). The effects of MT, SL, and BR on the biomass of PPC are independent of changes in microbial community structure. Nevertheless, it is connected to the succession of specific key biomarker species affected by each hormone.
5 Conclusion
In summary, the study of the endophytic fungal microbiome in FS, SSD, and SS stages of PPC indicates that the seedling stage exhibits the greatest microbial richness and diversity, along with a significant shift in composition. We also identified Sordariomycetes as essential microorganisms in this process, potentially contributing to the establishment and growth of PPC seedlings. Moreover, this study explored the effects of phytohormones on the phenotypic, physiological characteristics, and endophytic fungal community of PPC. Our findings showed that MT and SL significantly boosted PPC biomass by 69.32 and 15.23%, respectively. However, 2 mg/L of BR hindered PPC root development. MT, SL, and BR notably induce changes in the endophytic fungal community composition and diversity within PPC, potentially affecting biomass through the specific regulation of potential biomarkers. Moreover, phytohormones can indirectly impact the endophytic fungal community by updating the plant antioxidant system. This research provides novel insights into the dynamics of microbial communities within FS, SSD, and SS. The variations among distinct phytohormones may ultimately enhance our predictive understanding regarding strategies to control plant microbiota, either directly or indirectly, for improving plant health.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA1193989.
Author contributions
TP: Investigation, Methodology, Writing – original draft. TY: Funding acquisition, Resources, Supervision, Writing – review and editing. JS: Data curation, Project administration, Writing – review and editing. JZ: Formal Analysis, Investigation, Writing – review and editing. JS: Investigation, Writing – original draft.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was funded by the Key research and development project of Gansu Academy of Sciences (2023ZDYF-01), Industrialization application Project of Gansu Academy of Sciences (2024YLCY-01).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The authors declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
References
1
Al-BabiliS.BouwmeesterH. J. (2015). Strigolactones, a novel carotenoid-derived plant hormone.Annu. Rev. Plant Biol.66161–186. 10.1146/annurev-arplant-043014-114759
2
BackK. (2021). Melatonin metabolism, signaling and possible roles in plants.Plant J.105376–391. 10.1111/tpj.14915
3
BarretM.BriandM.BonneauS.PréveauxA.ValièreS.BouchezO.et al (2015). Emergence shapes the structure of the seed microbiota.Appl. Environ. Microbiol.811257–1266. 10.1128/aem.03722-14
4
BerriosL.ElyB. (2020). Plant growth enhancement is not a conserved feature in the Caulobacter genus.Plant Soil44981–95. 10.1007/s11104-020-04472-w
5
BerriosL.RentschJ. D. (2022). Linking reactive oxygen species (ROS) to abiotic and biotic feedbacks in plant microbiomes: The dose makes the poison.Int. J. Mol. Sci.231422–1467. 10.3390/ijms23084402
6
Carrera-CastañoG.Calleja-CabreraJ.PernasM.GómezL.Oñate-SánchezL. (2020). An updated overview on the regulation of seed germination.Plants (Basel)92223–7747. 10.3390/plants9060703
7
CarvalhaisL. C.Rincon-FlorezV. A.BrewerP. B.BeveridgeC. A.DennisP. G.SchenkP. M. (2019). The ability of plants to produce strigolactones affects rhizosphere community composition of fungi but not bacteria.Rhizosphere918–26. 10.1016/j.rhisph.2018.10.002
8
ChenC.ZouJ.ZhangS.ZaitlinD.ZhuL. (2009). Strigolactones are a new-defined class of plant hormones which inhibit shoot branching and mediate the interaction of plant-AM fungi and plant-parasitic weeds.Sci. China C Life Sci.52693–700. 10.1007/s11427-009-0104-6
9
CunninghamA. B.BrinckmannJ. A.BiY. F.PeiS. J.SchippmannU.LuoP. (2018). Paris in the spring: A review of the trade, conservation and opportunities in the shift from wild harvest to cultivation of Paris polyphylla (Trilliaceae).J. Ethnopharmacol.222208–216. 10.1016/j.jep.2018.04.048
10
DarréM.ZaroM. J.Guijarro-FuertesM.CareriL.ConcellónA. (2024). Melatonin combined with wax treatment enhances tolerance to chilling injury in red bell pepper.Metabolites142218–1989. 10.3390/metabo14060330
11
DevireddyA. R.ZandalinasS. I.FichmanY.MittlerR. (2021). Integration of reactive oxygen species and hormone signaling during abiotic stress.Plant J.105459–476. 10.1111/tpj.15010
12
EichmannR.RichardsL.SchäferP. (2021). Hormones as go-betweens in plant microbiome assembly.Plant J.105518–541. 10.1111/tpj.15135
13
Franke-WhittleI. H.ManiciL. M.InsamH.StresB. (2015). Rhizosphere bacteria and fungi associated with plant growth in soils of three replanted apple orchards.Plant Soil395317–333. 10.1007/s11104-015-2562-x
14
FullerA. W.YoungP.PierceB. D.Kitson-FinuffJ.JainP.SchneiderK.et al (2017). Redox-mediated quorum sensing in plants.PLoS One12:e0182655. 10.1371/journal.pone.0182655
15
Garcia-ArocaT.PriceP. P.Tomaso-PetersonM.AllenT. W.WilkersonT. H.SpurlockT. N.et al (2021). Xylaria necrophora, sp. nov., is an emerging root-associated pathogen responsible for taproot decline of soybean in the southern United States.Mycologia113326–347. 10.1080/00275514.2020.1846965
16
GuanL.ZhengZ.GuoZ.XiaoS.LiuT.ChenL.et al (2024). Steroidal saponins from rhizome of Paris polyphylla var. chinensis and their anti-inflammatory, cytotoxic effects.Phytochemistry219:113994. 10.1016/j.phytochem.2024.113994
17
HeX.YinB.ZhangJ.ZhouS.LiZ.ZhangX.et al (2023). Exogenous melatonin alleviates apple replant disease by regulating rhizosphere soil microbial community structure and nitrogen metabolism.Sci. Total Environ.884:163830. 10.1016/j.scitotenv.2023.163830
18
LiaoK.ChuT.ShiY.XieX.LiL.FanT.et al (2024). First report of gibellulopsis nigrescens causing yellow wilt on chinese cabbage in China.Plant Dis.108:3197. 10.1094/pdis-05-24-0989-pdn
19
LinksM. G.DemekeT.GräfenhanT.HillJ. E.HemmingsenS. M.DumonceauxT. J. (2014). Simultaneous profiling of seed-associated bacteria and fungi reveals antagonistic interactions between microorganisms within a shared epiphytic microbiome on riticum and rassica seeds.New Phytol.202542–553. 10.1111/nph.12693
20
LiuY.YanH.ZhangX.ZhangR.LiM.XuT.et al (2020). Investigating the endophytic bacterial diversity and community structures in seeds of genetically related maize (Zea mays L.) genotypes.3 Biotech10:27. 10.1007/s13205-019-2034-8
21
LvB.TianH.ZhangF.LiuJ.LuS.BaiM.et al (2018). Brassinosteroids regulate root growth by controlling reactive oxygen species homeostasis and dual effect on ethylene synthesis in Arabidopsis.PLoS Genet.14:e1007144. 10.1371/journal.pgen.1007144
22
MarzecM. (2016). Strigolactones as part of the plant defence system.Trends Plant Sci.21900–903. 10.1016/j.tplants.2016.08.010
23
NelsonE. B. (2018). The seed microbiome: Origins, interactions, and impacts.Plant Soil4227–34. 10.1007/s11104-017-3289-7
24
OkunishiS.SakoK.ManoH.ImamuraA.MorisakiH. (2005). Bacterial Flora of Endophytes in the Maturing Seed of Cultivated Rice (Oryza sativa).Microbes Environ.20168–177. 10.1264/jsme2.20.168
25
RenY.CheX.LiangJ.WangS.HanL.LiuZ.et al (2021). Brassinosteroids benefit plants performance by augmenting arbuscular mycorrhizal symbiosis.Microbiol. Spectrum.9:e0164521. 10.1128/spectrum.01645-21
26
ShadeA.JacquesM.-A.BarretM. (2017). Ecological patterns of seed microbiome diversity, transmission, and assembly.Curr. Opin. Microbiol.3715–22. 10.1016/j.mib.2017.03.010
27
ShahzadR.KhanA. L.BilalS.AsafS.LeeI. J. (2018). What is there in seeds? vertically transmitted endophytic resources for sustainable improvement in plant growth.Front. Plant Sci.9:24. 10.3389/fpls.2018.00024
28
ShaoQ.RanQ.LiX.DongC.ZhangY.HanY. (2024). Differential responses of the phyllosphere abundant and rare microbes of Eucommia ulmoides to phytohormones.Microbiol. Res.286:127798. 10.1016/j.micres.2024.127798
29
SongX. E.CaoJ.GuoS.WangH.DongQ.GuoP.et al (2023). Effect of brassinolide on soil microorganisms in millet field polluted by tribenuron-methyl.Microorganisms112076–2607. 10.3390/microorganisms11071829
30
State Pharmacopoeia Commission of the Prc. (2020). Pharmacopoeia of People’s Republic of China.China Med. Sci. Technol. Press1274–275.
31
StringlisI. A.PieterseC. M. J. (2021). Evolutionary “hide and seek” between bacterial flagellin and the plant immune system.Cell Host Microbe29548–550. 10.1016/j.chom.2021.03.010
32
SytarO.KumariP.YadavS.BresticM.RastogiA. (2019). Phytohormone priming: Regulator for heavy metal stress in plants.J. Plant Growth Regul.38739–752. 10.1007/s00344-018-9886-8
33
TalaatN. B.ShawkyB. T. (2012). 24-Epibrassinolide ameliorates the saline stress and improves the productivity of wheat (Triticum aestivum L.).Environ. Exp. Bot.8280–88. 10.1016/j.envexpbot.2012.03.009
34
TanD. X.HardelandR.ManchesterL. C.KorkmazA.MaS.Rosales-CorralS.et al (2012). Functional roles of melatonin in plants, and perspectives in nutritional and agricultural science.J. Exp. Bot.63577–597. 10.1093/jxb/err256
35
TianY.FanM.QinZ.LvH.WangM.ZhangZ.et al (2018). Hydrogen peroxide positively regulates brassinosteroid signaling through oxidation of the Brassinazole-Resistant1 transcription factor.Nat. Commun.9:1063. 10.1038/s41467-018-03463-x
36
Torres-CortésG.BonneauS.BouchezO.GenthonC.BriandM.JacquesM. A.et al (2018). Functional microbial features driving community assembly during seed germination and emergence.Front. Plant Sci.9:902. 10.3389/fpls.2018.00902
37
TrouillasF. P.GublerW. D. (2016). Cryptosphaeria dieback of fremont cottonwood caused by Cryptosphaeria pullmanensis and C. multicontinentalis in California.Plant Dis.100777–783. 10.1094/pdis-09-15-0972-re
38
TruyensS.WeyensN.CuypersA.VangronsveldJ. (2015). Bacterial seed endophytes: Genera, vertical transmission and interaction with plants.Environ. Microbiol. Rep.740–50. 10.1111/1758-2229.12181
39
WangX.PengC.LiangJ.LiangQ.XuC.GuoW. (2019). The complete chloroplast genome of Paris polyphylla var. chinensis, an endemic medicinal herb in China.Mitochondrial DNA B Resour.43888–3889. 10.1080/23802359.2019.1687351
40
WangX.WanasingheD. N.ZhangJ.MaJ.ZhouP.ZhangL.et al (2023). Insights from the Endophytic Fungi in Amphisphaeria (Sordariomycetes): A. orixae sp. nov. from Orixa japonica and its secondary metabolites.Microorganisms112076–2607. 10.3390/microorganisms11051268
41
XuG.YangS.MengL.WangB.-G. (2018). The plant hormone abscisic acid regulates the growth and metabolism of endophytic fungus Aspergillus nidulans.Sci. Rep.8:6504. 10.1038/s41598-018-24770-9
42
YanC.HaoH.WangZ.ShaS.ZhangY.WangQ.et al (2023). Prediction of suitable habitat distribution of Cryptosphaeria pullmanensis in the world and China under climate change.J. Fungi (Basel)9739. 10.3390/jof9070739
43
YanD.WangJ.LuZ.LiuR.HongY.SuB.et al (2023). Melatonin-mediated enhancement of photosynthetic capacity and photoprotection improves salt tolerance in wheat.Plants (Basel)122223–7747. 10.3390/plants12233984
44
YangF.JiangH.ChangG.LiangS.MaK.CaiY.et al (2023). Effects of rhizosphere microbial communities on cucumber fusarium wilt disease suppression.Microorganisms11:1576. 10.3390/microorganisms11061576
45
YeF.JiangM.ZhangP.LiuL.LiuS.ZhaoC.et al (2022). Exogenous melatonin reprograms the Rhizosphere microbial community to modulate the responses of barley to drought stress.Int. J. Mol. Sci.23:9665. 10.3390/ijms23179665
46
ZhangJ.HuH.LiS.ShangW.JiangJ.XuX.et al (2023). Diversity of fungal endophytes in american ginseng seeds.Plant Dis.1072784–2791. 10.1094/PDIS-10-22-2312-RE
47
ZhangX.HuangT.LiangY.HussainS.PengR.WangT.et al (2024). Melatonin and 14-hydroxyed brassinosteroid combined promote kiwifruit seedling growth by improving soil microbial distribution, enzyme activity and nutrients uptake.Front. Plant Sci.15:1336116. 10.3389/fpls.2024.1336116
Summary
Keywords
Paris polyphylla var. chinensis (Franch.) Hara, endophytic fungi, phytohormones, phenotypic, antioxidant system
Citation
Peng T, Yang T, Sha J, Zhao J and Shi J (2025) Dynamics of endophytic fungi composition in paris polyphylla var. chinensis (franch.) hara seeds during storage and growth, and responses of seedlings to phytohormones. Front. Microbiol. 16:1540651. doi: 10.3389/fmicb.2025.1540651
Received
06 December 2024
Accepted
14 February 2025
Published
26 February 2025
Volume
16 - 2025
Edited by
Anukool Vaishnav, GLA University, India
Reviewed by
Zeba I. Seraj, University of Dhaka, Bangladesh
Prachi Singh, Banaras Hindu University, India
Updates
Copyright
© 2025 Peng, Yang, Sha, Zhao and Shi.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tao Yang, yangtao8383@163.comJianwu Shi, jwshi@ntu.edu.cn
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.