Next Article in Journal
Copper Sulfate Combined with Photodynamic Therapy Enhances Antifungal Effect by Downregulating AIF1
Previous Article in Journal
Chromosome-Level Assembly and Comparative Genomic Analysis of Suillus bovinus Provides Insights into the Mechanism of Mycorrhizal Symbiosis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diversity, Community Structure, and Antagonism of Endophytic Fungi from Asymptomatic and Symptomatic Mongolian Pine Trees

1
Beijing Key Laboratory of Forest Pest Control, Beijing Forestry University, Beijing 100083, China
2
Honghuaerji Forestry Field, Honghuaerji Forestry Bureau, Hulunbeier 165456, China
*
Author to whom correspondence should be addressed.
J. Fungi 2024, 10(3), 212; https://doi.org/10.3390/jof10030212
Submission received: 24 January 2024 / Revised: 4 March 2024 / Accepted: 5 March 2024 / Published: 13 March 2024
(This article belongs to the Section Fungal Pathogenesis and Disease Control)

Abstract

:
Diplodia tip blight, caused by Diplodia sapinea (=Sphaeropsis sapinea), are widely distributed in Honghuaerji, Inner Mongolia, China, causing severe damage on natural Mongolian pine (Pinus sylvestris var. mongolica). D. sapinea is an endophyte that becomes pathogenic under conditions of drought, hail damage, or temperature-associated stress. The role of the endophytic community inhabiting different pine tissues in the expression of disease is still unknown. In this study, the diversity and community structure of endophytic fungi among asymptomatic and symptomatic Mongolian pine were detected using culture-based isolation and high-throughput sequencing (HTS), and the potential antagonistic endophytes against D. sapinea were also screened. The results indicated that 198 and 235 strains of endophytic fungi were isolated from different tissues of symptomatic and asymptomatic Mongolian pine, respectively. D. sapinea was the most common endophyte isolated from the current-year needles and shoots of symptomatic trees, and Diplodia was also the most common in the HTS data. There were no significant differences in the endophytic fungal species richness among asymptomatic and symptomatic trees, but there were differences observed within specific sampled tissues. The ANOSIM analysis confirmed that the endophytic fungi community structure significantly differed between sampling tissues among symptomatic and asymptomatic Mongolian pine. Furthermore, the antagonism study revealed Penicillium fructuariae-cellae with the ability to inhibit the growth of D. sapinea in vitro, and the potential performance of this fungus, acting as biological control agent, was evaluated under greenhouse. Our findings can pave the way to a better understanding of the interactions between D. sapinea, other endophytic fungi and their hosts, and provide helpful information for more efficient disease management strategies.

1. Introduction

Mongolian pine (Pinus sylvestris var. mongolica), a significant pine species indigenous to Asia, has been extensively cultivated in Liaoning, Inner Mongolia, Heilongjiang, and Hebei Provinces of China, due primarily to its rapid growth and exceptional adaptation to cold climates [1]. Since 2018, Diplodia tip blight, caused by Diplodia sapinea (Fr.) Fuckel, has been widely distributed in Honghuaerji, Inner Mongolia, China, causing severe damage on natural Mongolian pine. The first evidence of Diplodia tip blight on Mongolian pine was reported in Heilongjiang Province, China in 1980s [2]. This disease is also observed on various Pinus species, especially Scots pine (Pinus sylvestris), across parts of Europe such as Germany, Sweden, and Finland, and its incidence has significantly increased over the past decades [3,4,5,6,7,8,9]. Climate change is expected to exacerbate the damage caused by this fungus, as it flourishes in higher temperatures. The reported increase in temperature and drought associated with climate change would also make pine more susceptible to D. sapinea, potentially leading to more frequent disease outbreaks. Currently, there is no established and effective method for the disease control.
The pathogen D. sapinea, is also known as Sphaeropsis sapinea (Fr.: Fr.) Dyko and Sutton. The correct name of this fungus is still in discussion, and the current name in Index Fungorum is S. sapinea. The preferred name after the EPPO Global Database, however, is D. sapinea [6,8,10] D. sapinea is indeed an endophytic fungus that can transition to a pathogenic state under specific environmental conditions, particularly during periods of drought or temperature-associated stress, leading to the development of symptoms, such as tip blight, dieback of current year shoots, stem cankers, blue staining of the sapwood, and ultimately, the decline or death of the pine tree [4,8].
Endophytic fungi are a group of fascinating host-associated fungal communities that live asymptotically inside host tissues for the entire or at least a significant part of their life cycle, without causing apparent negative harm to the host [6,11,12]. It is widely accepted that these fungi have co-evolved with their host trees, and often establish mutualistic relationships that can influence the host fitness [12,13]. The endophytic stage represents a balanced interaction between the fungus and its host. However, endophytic fungi can become pathogen when this balance is disturbed [3,6]. Climate change may influence the lifestyle switch of D. sapinea. When climate shifts towards conditions less favorable for host trees (e.g., drought/hailstorm damage), and more conducive for pathogenic fungi (e.g., warmer temperatures), the expression of disease was triggered [4,14,15,16]. The mechanism underlying D. sapinea’s transition into a pathogenic state, inducing disease expression, remains incompletely understood.Sherwood et al. [17] proposed that the disease expression may be linked to D. sapinea’s ability to exploit metabolites produced by the host in response to stress, such as proline, using these metabolites as a nitrogen source for growth. Furthermore, in order to investigate the role of competitive interactions among endophytes and latent pathogens in determining disease in stressed pine, Oliva et al. [8] concluded that hail promoted D. sapinea and other endophytes with a rapid colonization strategy of N-rich substrates. Rapid niche occupation is critical for D. sapinea to cause disease after a stress-inducing event. However, competition with other endophytes for key metabolites can suppress the pathogen and prevent trees from developing symptoms.
Despite extensive investigation into the endophytic fungal community of Scots pine, the intricate interactions among different endophytes and their host, as well as with other microbiomes, have yet to be fully elucidated [3,4,6,18]. Bußkamp et al. [4] isolated a total of 103 outgrowing endophytic fungal species from twigs of healthy and diseased Scots pine. Interestingly, the comparison of fungal endophyte communities between twigs from diseased and non-diseased Scots pine trees revealed few differences, except D. sapinea, which was the most common endophyte. Blumenstein et al. [6] conducted a comparison of the mycobiome in twigs of healthy and diseased Scots pines using high-throughput sequencing (HTS) of the ITS2 region. They indicated that the mycobiome community composition significantly differed between growth years and sampling time, but not between healthy and diseased trees. The composition of endophytic fungal communities exhibited considerable variation among different host tissues due to differences in microenvironments, nutrient availability, and physiological conditions [19]. Kovalchuk et al. [20] analyzed the mycobiome of different anatomic parts of asymptomatic and symptomatic Norway spruce trees naturally infected by Heterobasidion, and they demonstrated that not only does each individual tree tissue (wood, bark, needles, and roots) harbor a unique fungal community, but also that the structure of fungal communities residing in the wood differed significantly among symptomatic and asymptomatic trees. In that sense, differences in the mycobiota among distinct tissues may influence the processes leading to disease outbreak and various disease symptoms. Furthermore, endophytic fungi residing in different tissues of hosts may demonstrate diverse functional traits. For example, root endophytes might directly influence nutrient uptake and root architecture, while leaf endophytes could primarily impact foliar physiology and defense responses [21,22]. The endophytic community, and competition between endophytes can play a role in the Diplodia tip blight disease expression of different tissue (e.g., needles, shoots). Therefore, comprehending the variances in endophytic fungi across diverse host tissues is crucial for deciphering the intricate dynamics of plant-microbe relationships and maximizing their potential advantages in disease control.
The recognition of fungal endophytes’ capacity to enhance host fitness and stress tolerance has sparked the notion of employing these organisms as biocontrol agents [23,24,25]. Studies on conifer trees have demonstrated that inoculations with fungal endophytes can protect the host from natural infections by other pathogens, and these endophytes can be used as antagonists against potential pathogens. Various metabolites with antifungal properties have also been extracted from foliar endophytes of Picea [26,27,28]. The antagonism study revealed 13 possible endophytic fungi with the ability to inhibit the growth of D. sapinea in vitro, for example Sydowia polyspora (Bref. & Tavel) E. Müll. [18]. Oliva et al. [8] also deliberated on additional endophytic fungi with potential antagonistic properties, such as Alternaria sp. and Epicoccum nigrum Link.
The objectives of this study were to (i) examine differences in the composition of the endophytic fungal communities associated with various tissues of asymptomatic and symptomatic Mongolian pine, using culture-based isolation and high-throughput sequencing (HTS); (ii) identify the presence of antagonist endophytes in asymptomatic pine; and (iii) screen antagonist endophytes with the ability to inhibit the growth of D. sapinea in vitro and under greenhouse inoculation experiment.

2. Materials and Methods

2.1. Study Sites and Sample Collection

Two natural forest sites of Mongolian pine (Baogentu Forest Farm and Toudaoqiao Forest Farm) in Honghuaerji, Hulunbuir, Inner Mongolia Autonomous Region, China, were chosen for sampling. The two sites are 70 km apart and growing under similar conditions (with a cold temperate continental climate: −1.5 °C–3.7 °C annual mean temperature and 260–490 mm annual mean precipitation, an elevation of 95 m above sea level, and predominantly sandy soils) (https://data.cma.cn/) (accessed on 15 January 2024). The pine stands at both sites were naturally regenerated and approximately 40 years old at the time of sampling.
In July 2022, 30 sampled trees in Baogentu Forest Farm (119°55′57″ E, 48°24′12″ N) (Figure 1A) were classified as ‘asymptomatic’ (<10% crown displaying blight symptoms), while the 30 sampled trees from Toudaoqiao Forest Farm (120°13′38″ E, 48°16′44″ N) (Figure 1B) were considered symptomatic (with 30% or higher percentage of Diplodia tip blight symptoms) [8]. A total of 180 samples were collected, including current-year needles, current-year shoots, and phloem (at breast height) (Figure 1C) from 30 asymptomatic and 30 symptomatic trees. These samples obtained from the current-year needles, current-year shoots, and phloem of asymptomatic trees were labeled as HA, HB, HC, while those from symptomatic trees were labeled as DA, DB, DC. The diameter at breast height of the selected trees ranged from 20 to 25 cm. The samples for culture-based isolation were kept 4 °C and analyzed within 2 days, and the samples for high-throughput sequencing was immediately placed in liquid nitrogen at the site and stored at −80 °C.

2.2. Culture-Based Isolation and Molecular Identification

Small sections (5 × 5 mm) were cut from the collected tissue samples, surface-sterilized in 70% ethanol for 1 min, sterilized in 3% (v/v) sodium hypochlorite solution for 5 min and in 70% ethanol for 30 s, followed by three rinses in sterile distilled water and finally dried on sterilized filter paper [8]. In total, 4–5 cutting segments were then placed on the 0.3% malt extract agar (MEA) media amended with chloramphenicol. All cultures were incubated at 25 °C in the dark for 3 days. Distinct colonies were picked for further purification on 2% MEA media and incubated for seven days [8]. The colony characteristics of each isolate were observed and recorded. One representative per morphotype was used for molecular identification.
Mycelium was scraped from the surface of the colony and DNA was extracted using the cetyltrime thylammonium bromide (CTAB) method [29]. ITS region was amplified using universal ITS primers ITS1/ITS4 [30]. The PCR mixture consisted of 10 μL of 2 × Hieff Master Mix, 7 μL of nucleic acid-free H2O, 1 μL of each primer (10 μM), and 1 μL of DNA samples were made up to the final volume of 20 μL. The PCR procedure consisted of an initial denaturation step at 94 °C for 5 min, followed by 35 cycles comprising 30 s at 94 °C, 50 s at 48 °C, and 1 min at 72 °C, with a final elongation step of 7 min at 72 °C. After the PCR amplification, the PCR products were delivered to Beijing Nuosai Genome Research Center Co., Ltd. for sequencing. The DNA sequences were cleaned and blasted in GenBank. The OTU identities were assigned according to match thresholds of >98.5% for species, >97% for genus, >95% for family, >92% for order, >90% for class, and >80% for phylum [8]. The isolated endophytes were accurately identified based on a combination of their morphological characteristics. All isolated strains were stored at 4 °C.

2.3. High-Throughput Sequencing

2.3.1. DNA Extraction, Amplification of ITS1 Region, and Sequencing

A total of 60 tissue samples from 10 asymptomatic and 10 symptomatic trees were surface-sterilized, defrosted, and re-sterilized for 1 min in a 3% NaOCl solution. DNA was extracted from 30 mg of the homogenized tissue samples using the TGuide S96 Magnetic Soil /Stool DNA Kit (Tiangen Biotech (Beijing, China) Co., Ltd.) according to the manufacturer’s instructions. The quantification of DNA samples was performed using the Qubit dsDNA HS Assay Kit and Qubit 4.0 Fluorometer (Invitrogen, Thermo Fisher Scientific, Eugene, OR, USA). The primer pair of ITS1F/2 (ITS1:5′-CTTGGTCATTTAGAGGAAGTAA-3′ and ITS2:5′-GCTGCGTTCTTCATCGATGC-3′) was used to amplify the ITS1 region. The PCR was performed in a total reaction volume of 10 μL that included 50 ng of DNA template, 0.3 μL of each primer (10μM), 5 μL of KOD FX Neo buffer, 1 μL of each dNTP2 (2 mM), 0.2 μL of KOD FX Neo, and ddH2O to make up the volume to 10 μL. The PCR conditions were as follows: an initial denaturation at 95 °C for 5 min, followed by 25 cycles of denaturation at 95 °C for 30 s, annealing at 50 °C for 30 s, and extension at 72 °C for 40 s, and a final step at 72 °C for 7 min. The PCR products were purified with Agencourt AMPure XP Beads (Beckman Coulter, Indianapolis, IN, USA) and sequenced with Illumina Novaseq 6000 (Illumina, Santiago CA, USA).

2.3.2. Preprocessing and Analysis of ITS1 Sequences

Raw sequences were preprocessed using the BMK Cloud (Biomarker Technologies Co., Ltd., Beijing, China). The raw reads obtained were filtered using Trimmomatic v0.33, and the primer sequences were removed using Cutadapt v1.9.1 to obtain high-quality clean reads. To acquire final valid data (effective reads), sequences from each sample were spliced and length-filtered using Usearch v10, and chimeric sequences were detected and removed using UCHIME v4.2. Clean reads were then denoised using the DADA2 method [31] in QIIME2 [32], which generated amplicon sequence variants (ASVs). Taxonomic analysis of the ASV sequences was performed using a simple Bayesian classifier with UNITE as the reference database. This produced taxonomic data about the species associated with each feature, which subsequently made it possible to analyze the composition of the microbial community at several levels, including phylum, class, order, family, genus, and species.
A rarefaction curve based on sample depth was generated to assess alpha diversity, with ACE, Chao1, Simpson, and Shannon indices employed to gauge species richness and community diversity of endophytic fungi across different tissues of symptomatic and asymptomatic pines. Chao1 and ACE indices were used to estimate species richness, with higher values indicating greater richness. Shannon and Simpson indices were utilized to measure community diversity, where higher values indicate a more diverse community [33]. Beta diversity was assessed to determine the similarity of endophytic fungal communities across samples using QIIME2. Principal coordinate analysis (PCoA) was employed to visualize the endophytic fungal community structure based on Bray–Curtis similarity. PERMANOVA test and analysis of similarities (Anosim) were conducted to determine the significant differences in community structure between different tissues in the pine trees. Additionally, LEfSe analysis was utilized to detect statistically significant differences between various tissues among symptomatic and asymptomatic pines. The linear discriminant analysis (LDA) > 4 and p < 0.05 was were considered significantly enriched in that group compared to other groups. FUNGuild was employed to predict the functional profiles of endophytic fungal communities. Endophytes were categorized into three groups based on their nutritional mode: pathotroph (which acquires nutrients by damaging host cells), symbiotroph (which obtains nutrients through resource exchange with host cells), and saprotroph (which obtains nutrients by decomposing deceased host cells).

2.4. Antagonism Assay

To identify endophytic fungal isolates with the ability to antagonize the pathogen D. sapinea, all endophytic fungal isolates were cultured on MEA at 25 °C for seven days prior to the antagonism experiment in vitro. Mycelial plugs (5 mm in diameter) of the endophytic fungi and the pathogenic fungi D. sapinea were placed 4 cm apart on the surface of MEA media [34]. For the control groups (CK), mycelial plugs of D. sapinea were placed in the center of MEA media. The Petri dishes were then incubated in the dark at 25 °C. Observations were recorded every 24 h. All treatment experiments were carried out in triplicate. The ability of an endophyte to antagonize the pathogen was determined based on the inhibition level (defined as pathogen growth with and without the endophyte) over a given period of time [18], and it was categorized into four groups: (1) inhibition of D. sapinea growth with the presence of a reaction zone (Figure 2A), (2) endophyte superior = the endophyte has overgrown and inhibited mycelial growth (Figure 2B), (3) neutral or ”mutually intermingling growth” = equal mycelial growth capacity between the endophyte and D. sapinea, or no apparent mutual inhibition between them (Figure 2C), and (4) D. sapinea superior = overgrowth of D. sapinea inhibited endophyte mycelial growth (Figure 2D). Growth was assessed 3 and 7 days after inoculation. The inhibition rate was calculated using the formula:
  Inhibition   rate   = ( R 1 R 2 ) R 1 × 100 %  
R1 refers to the radius of the control pathogen colony; R2 refers to the radius of the pathogen colony pointing towards the antagonist fungus.
One-way ANOVA analysis in SPSS 24.0 (IBM Corp., New York, NY, USA) was used to compare the inhibition rates and screen the antagonistic fungi with the highest inhibition rate.

2.5. Co-Culture Experiments with Preferential Placement of Endophytic Fungi

After screening for endophytic fungi with the highest inhibition rate, D. sapinea and the selected endophytic fungi were cultivated on MEA for three days in the dark at 25 °C. A mycelial plug of the pathogenic fungi D. sapinea was placed in the center of the Petri dish, while two mycelial plugs of the selected endophytic fungi were positioned on either side of the Petri dish, 2 cm away from the center (Figure 3). The treatments were as follows: (1) mycelial plugs of the selected endophytic fungi were placed on the MEA medium for 12 h, 24 h, 36 h, and 48 h, followed by the positioning of mycelial plug of D. sapinea; (2) simultaneously placed mycelial plugs of the selected endophytic fungi and D sapinea onto the MEA medium. For the control group (CK), only mycelial plugs of D. sapinea were placed onto the MEA medium. All cultures were then incubated in the dark at 25 °C for seven days. Colony growth was measured on the seventh day, and the inhibition rate of the selected endophyte was calculated using the same formula as mentioned.

2.6. Preparation of Antagonistic Endophytic Fungal Fermentation Broth

The endophytic fungus exhibiting the highest inhibition rate against D. sapinea was chosen and incubated in the dark at 25 °C for three days. Four mycelial plugs (5 mm in diameter) were then placed in a 250 mL conical flask containing 100 mL of autoclaved MEB liquid media. Fermentation was carried out at 25 °C and 180 rpm on a shaker for 7 days. Following filtration of the fermentation broth through gauze, it was transferred to 2 mL centrifuge tubes and centrifuged at 10,000 rpm for 10 min. 1.5 mL supernatant from each centrifuge tube was removed and filtered through a 0.22 um filter membrane to eliminate mycelium and spores, after which the filtered fermentation broth was stored at 4 °C [35].

2.7. Inoculation

Fifteen healthy three-year-old Mongolian pine seedlings were used for the greenhouse inoculation experiment. D. sapinea were cultured on 2% MEA and grown at 25 °C for 7 days prior to the experimental inoculations. Wounding inoculations were performed on the current year’s shoots, and the wounds were created using a sterilized needle. Lesion length was recorded daily after the inoculation.
In Group I, the wounded shoots of five pine seedlings were inoculated with mycelial plug (5 mm in diameter) of D. sapinea. The mycelium was positioned facing the wound, and the area was sealed with Parafilm to maintain moisture. In Group II, two days after the inoculation of D. sapinea, the fermentation broth (5 mL) of the selected antagonistic endophytic fungus, as mentioned in 2.6, was sprayed onto the current-year shoots of five pine seedlings. In Group III (CK), mock-control, the wounded shoots of five pine seedlings were inoculated with sterile MEA plugs (5 mm in diameter).

3. Results

3.1. Isolation and Identification of Endophytic Fungi from Mongolian Pine

A total of 433 strains were successfully isolated from 2700 tissue segments of 60 pine trees in the 2 sampling plots (Table 1). Among the 235 isolates isolated from the asymptomatic trees, 70 were from needles (HA), 108 from current-year shoots (HB), and 57 from phloem (HC). Meanwhile, the 198 isolates from 30 symptomatic trees included 64 from needles (DA), 86 from shoots (DB), and 48 from phloem (DC). These 433 isolates were initially classified into 57 representative morphotypes according to their cultural characteristics. Based on the ITS sequences similarity, 57 isolates were categorized at the genus level (Table 1; Figure 4). The majority of endophytic fungi were identified as Ascomycota (98.2%), which represented 34 genera. In addition, one genus belonged to Basidiomycota (Coprinopsis atramentaria). The genera of Sy. polyspora (22.01%), D. sapinea (15.22%), Cladosporium sp. (11.48%), and Penicillium sp. (10.54%) were most frequently isolated from the asymptomatic and symptomatic Mongolian Pine.
The most abundant species was Sy. polyspora (28.51%) in asymptomatic pines, followed by Cladosporium sp. (17.45%), and D. sapinea (11.49%). Sy. polyspora was also most commonly isolated from the current-year needles and shoots of asymptomatic pines, with 45.71% and 23.15%, respectively. Cladosporium sp., on the other hand, was most frequently found in the phloem (24.56%) of asymptomatic pines. In symptomatic pines, the most abundant species was D. sapinea (19.19%), followed by Sy. polyspora (13.64%), and Penicillium sp. (11.62%). Within the symptomatic shoots, D. sapinea was the most abundant, with a frequency of 25.00%. D. sapinea and Penicillium sp. were both isolated at a frequency of 15.63%, from the needles of symptomatic trees. In the phloem of symptomatic trees, Talaromyces pseudofuniculosus was the most frequently encountered at 18.75%, with D. sapinea following at a frequency of 16.67%.

3.2. HTS Analysis of Endophytic Fungal Community of Mongolian Pine

3.2.1. Analysis of Endophytic Fungal Community Composition of Mongolian Pine

A total of 28,022,878 raw reads were generated from 60 samples, and 27,872,009 clean reads were obtained after data cleaning. Each sample yielded a minimum of 253,437 clean reads and an average of 464,533 clean reads. The rarefaction curves reached stability with increasing sequencing volume, suggesting that the generated sequencing data adequately represent the community composition and structure of endophytic fungi in various tissues of both asymptomatic and symptomatic pines (Figure 5).
The reads were assigned to 30,130 amplicon sequence variants (ASVs). The sampled tissues of the pine trees shared 636 (2.11%) of the total ASVs (Figure 6). The proportion of the ASVs unique to a certain tissue ranged from 0.06% (1937 ASVs; DB) to 16.20% (4884 ASVs; HB), with the shoots of asymptomatic pine trees exhibiting the highest number of specific ASVs. According to the taxonomy annotation, these ASVs were classified into 17 phyla, 60 classes, 158 orders, 383 families, 954 genera, and 1707 species.
At the phylum level, Ascomycota and Basidiomycota were the predominant constituents of the endophytic fungal communities in various tissues of both asymptomatic and symptomatic pines (Figure 7A). The relative abundances of Ascomycota in DA, DB, DC, HA, HB, and HC were 83%, 98%, 79%, 64%, 82%, and 80%, respectively. Following Ascomycota, Basidiomycota had relative abundances of 12%, 22%, and 11% in DC, HA, and HC, respectively, while it accounted for less than 10% in DA, DB, and HB. Relative abundances of other phyla were below 10%.
Symptomatic and asymptomatic pines exhibited different dominant endophytic fungal taxa at the genus level. Diplodia abundance was statistically highest (48.88%) in symptomatic pines, followed by Others (27.28%). While in asymptomatic pines, the “Others” fungi had the highest relative abundance (47.30%), followed by Hormonema (14.63%) (Figure 7C). The endophytic fungal communities in various tissues of symptomatic and asymptomatic pines displayed different dominant taxa (Figure 7B). In symptomatic pines, Diplodia had higher relative abundance of 46.67% and 74.24% in current-year needles (DA) and shoots (DB), respectively. In the phloem (DC), the “Others” fungi had the highest relative abundance (49.80%), followed by Diplodia (25.39%). However, the “Others” was the most abundant in all three tissues of asymptomatic pines (HA: 66.52%, HB: 39.17%, and HC: 39.17%). But the abundance of Diplodia was statistically higher in the phloem (10.38%) compared to the other two tissues (HA: 5.33%, HB: 4.47%).

3.2.2. Alpha Diversity of Endophytic Fungal Community in Mongolian Pine

Alpha diversity was assessed using the ACE and Chao1 indices (Figure 8) to measure species richness and the Simpson and Shannon indices (Figure 9) to evaluate community diversity. The highest richness of endophytic fungal communities in asymptomatic and symptomatic trees were observed in current-year needles (HA) and in phloem (DC), respectively, with the current-year shoots in symptomatic (DB) and the phloem in asymptomatic trees (HC) having the lowest numbers of ASVs, respectively (Figure 8). There were significant differences observed in the species richness of endophytic fungi between shoots of asymptomatic and symptomatic pine trees (HB vs. DB). However, no significant differences were found in the species richness of endophytic fungi among asymptomatic and symptomatic trees in needles and phloem (HA vs. DA and HC vs. DC).
Needles had the highest community diversity in both symptomatic and asymptomatic trees. The lowest community diversity was observed in shoots of symptomatic trees (DB) and in phloem of asymptomatic trees (HC), respectively (Figure 9). There was statistically significant difference in endophytic fungal diversity among asymptomatic and symptomatic trees found in current-year needles and shoots, but no significant differences in phloem of asymptomatic and symptomatic trees (HC vs. DC).

3.2.3. Beta Diversity of Endophytic Fungal Community in Mongolian Pine

Beta diversity was analyzed using Principal Coordinate Analysis (PCoA), with statistical significance determined using ANOSIM. The PCoA, based on the relative abundance of ASVs, accounted for 41.38% of the observed variation (Figure 10). Significant differences were found in the endophytic fungal community structure between asymptomatic and symptomatic trees (H vs. D) (Figure S1). Additionally, ANOSIM analysis revealed significant differences in the endophytic fungal community structure among the sampled pine tissues (A vs. B vs. C) (p-value = 0.001) (Figure S2). Notably, there were significant differences related to needles and shoots among asymptomatic and symptomatic trees (HB vs. DB; HA vs. DA), but no distinction was observed between HC and DC (R2 = 0.08, p-value = 0.115) (Table 2).

3.2.4. LEfSe Analysis of Endophytic Fungal Community in the Shoots of Asymptomatic and Symptomatic Mongolian Pine

To identify the taxonomic fungi with significantly different abundances among shoots of asymptomatic and symptomatic Mongolian Pine, LEfSe analysis was used for biomarker analysis (LDA score > 4, p value < 0.05). As shown in Figure 11, Diplodia was representative genus (higher relative abundance) in the shoots of symptomatic pine, while Perusta, Cladosporium, and Alternaria were the significant genus in the shoots of asymptomatic pine.

3.2.5. Functional Prediction of Endophytic Fungal Community in Mongolian Pine

FUNGuild was used to assess the ecological and functional status of 30,130 ASVs obtained from the sampled trees, of which 12,870 ASVs (42.7%) were successfully defined. Across various tissues of both asymptomatic and symptomatic pines, saprotrophs constituted the largest proportion, followed by pathotrophs and symbiotrophs (Figure 12). The proportion of ASVs classified as saprotrophs within a certain tissue ranged from 63.15% (DA) to 92.49% (DB) (Figure 12). Additionally, Undefined saprotrophs exhibited the highest prevalence in all tissue samples (Figure 13).

3.3. In Vitro Study: Antagonism Assay

Fifty-seven endophytes isolated from the sampled pine trees were employed in the antagonism assay (Table 3). Three isolated endophytes (Penicillium fructuariae-cellae, Taifanglania sp., and Paracamarosporium hawaiiense) inhibited the growth of D. sapinea (Figure 14A and Figure 15). Two isolated endophytes (Aspergillus flavus and Phoma sp.) displayed similar growth capability (Figure 14B), while four isolated endophytes (Aspergillus sp., Tricharina sp., Trichoderma afroharzianum, and Trichoderma harzianum) exhibited superior growth over D. sapinea (Figure 14C). The remaining tested endophytes showed inferior growth against D. sapinea (Figure 14D).
Three selected antagonistic fungi, Taifanglania sp., Pa. hawaiiense, and Pe. fructuariae-cellae, were tested further. Taifanglania sp. was isolated from DA, while Pa. hawaiiense originated from HB, and Pe. fructuariae-cellae was obtained from HB, HC, DB, and DC. All three strains exhibited inhibition rates of over 60% on the mycelial growth of D. sapinea (Table 4). Pe. fructuariae-cellae (EMD32) and Pa. Hawaiiense (EMD29) demonstrated notably higher inhibition rates against D. sapinea compared to Taifanglania sp. (EMD48). However, the inhibition rates of Pe. fructuariae-cellae and Pa. Hawaiiense against D. sapinea did not show significant differences. Overall, Pe. fructuariae-cellae emerged as the most effective endophytic fungus in inhibiting the pathogen D. sapinea, and was selected for further analysis.

3.4. Co-Culture Experiment with Preferential Placement of Pe. fructuariae-cellae

Pe. fructuariae-cellae was initially cultured on MEA medium, after which D. sapinea was added to the culture after 12, 24, 36, and 48 h. The findings demonstrated that Pe. fructuariae-cellae consistently inhibited D. sapinea by over 70% at each tested time point. There was no notable variance in the inhibition rates at 24 and 36 h against D. sapinea. The 48-h experiment exhibited the most substantial inhibitory effect, with an inhibition rate of 92.86% (Figure 16).

3.5. Greenhouse Inoculation Experiment

Pe. fructuariae-cellae (EMD32) was chosen for the greenhouse inoculation experiment. The mock-inoculated control pine seedlings (Group III, CK) stayed healthy during the experiment, with no dead shoots were observed (Figure 17). In contrast, inoculation with D. sapinea alone (Group I) led to the highest shoot mortality (60%) among the inoculated pine twigs after 21 days of inoculation. Necrosis was significantly greater (5–6 mm after 7 days inoculation) in this treatment, compared to co-infection with the potential antagonistic fungus and D. sapinea. The 5 pine seedlings in Group II, co-infected with Pe. fructuariae-cellae and D. sapinea, remained healthy after 7 days of spraying, with an average lesion length of 2 mm, and no dead shoots and yellow needles were observed after 21 days of inoculation.

4. Discussion

4.1. Endophytic Fungi Detected by Culture-Based Isolation and HTS Method

Endophytic fungal diversity and community composition of different tissues among asymptomatic and symptomatic Mongolian pine were investigated using both culture-based isolation and HTS methods. Our study revealed a broader range of endophytic fungi (1707 species) detected in the HTS data compared to only 57 species isolated through cultivation methods. In contrast to the HTS data, only one species of Basidiomycota fungi was isolated by the culture-based method in this study. This may be partially explained by the fast growth of D. sapinea in culture [16,36], which possibly resulted in an under representation of slow-growing fungal species.
Previous studies have shown that HTS can detect a wide range of fungi, including rare and unculturable species, providing insights into fungal community structure and dynamics. Culture-based methods, on the other hand, allow for the isolation of individual strains, which grow rapidly on the provided nutrient medium, facilitating further studies like biocontrol experiments and pathogenicity tests [6,37,38]. Blumenstein et al. [6] indicated that 65% of the species isolated from Scots pine twigs using cultivation methods were detected in the HTS data, while some isolated species were not detected by HTS. In our study, 45% isolated endophytic fungi were obtained from the shoots of the asymptomatic trees (HB), and 43% from the shoots of the symptomatic trees (DB), while in the HTS data, HB exhibiting the highest number of ASVs, but DB with the lowest numbers of ASVs. The selection of primers and databases remains a limiting factor in fully uncovering all species [6]. Therefore, each approach exhibited benefits and limitations. By combining culture-based methods with HTS, a more comprehensive understanding of endophytic fungal communities can be obtained.
Similar to previous studies [4,6,14,39,40], D.sapinea was isolated from both asymptomatic and symptomatic pine trees, indicating its endophytic mode. It was the most common endophyte isolated from the symptomatic pine trees and the most common in the HTS data of the symptomatic pine trees. However, its relative abundance varied among various tissues, with D. sapinea being most abundant in the current-year shoots of symptomatic pine trees (DB) and less abundant in the shoots and phloem of asymptomatic pine trees (HB & HC). The observations might indicated that D. sapinea accumulated mainly in the shoot/twigs and played a role in determine disease development. Bußkamp et al. [4] indicated that during periods of stress such as drought, common endophytes within the Scots pine tree may faced disadvantages, allowing D. sapinea to thrive as a main pathogen. It grew into the peridermis and cortex, eventually reaching the vascular tissues in diseased plants, and can quickly occupy more tissues in the host. This could explain why D. sapinea is found in higher abundance in the twigs of diseased pines.
In our study, Sy. polyspora was isolated from all tissues of both asymptomatic and symptomatic pine trees, with the highest frequency observed in the asymptomatic pine trees. It was also the second most common endophyte isolated from the symptomatic pine trees. Hormonem, which includes the species Sy. polyspora (formerly known as Hormonema dematioides Lagerb. and Melin), was the second most common genus in asymptomatic pine trees according to the HTS data. Additionally, Sy. polyspora exhibited the highest isolation frequency in the current-year needles of asymptomatic pine trees (HB), and the lowest frequency in the phloem of symptomatic pine trees (DC). Several studies indicated that Sy. polyspora, as a common foliar endephyte of Scots pine, was the most abundant fungus in all disease classes identified with HTS [3,41], and was also the second most common fungus identified by the culture-based method in Blumenstein et al.’s (2021) [18] study. However, Sy. polyspora is also known to infect Pinus pinea L. in Portugal [42] and Pinus yunnanensis Franch. in Southwestern China [43], resulting in symptoms like tip dieback, needles with tan- to yellow-colored lesions, and eventual shoot death.

4.2. Endophytic Fungi Community in Different Tissues among Asymptomatic and Symptomatic Mongolian Pine

Alpha diversity analysis unveiled higher species richness and community diversity of endophytic fungi in asymptomatic Mongolian pine compared to symptomatic pine (H vs. D) (Figures S3 and S4). Despite no statistically significant differences in richness indices between asymptomatic and symptomatic trees (Figure S3), significant differences were observed in diversity indices, with the Simpson index (p = 0.0012) and Shannon index (p = 0.0027) (Figure S4). This implies that although the number of different endophytic fungi may not vary greatly between the two groups, there are notable differences in the composition and evenness of species within the asymptomatic and symptomatic pine. Furthermore, the structure of endophytic fungal communities differed significantly among symptomatic and asymptomatic pine tree. It could be attributed to differences in the presence and abundance of core endophytes. It is possible that the ability of the mycobiome to cause disease symptoms depends on the dominance of certain endophytic fungal species. In other words, the relative abundance and interactions among specific endophytic fungal species may play a crucial role in determining the health status of the pine trees. Differing findings in Scots pine (P. sylvestris) were reported by Blumenstein et al. [6], suggesting that diversity indices did not exhibit statistical differences between disease classes, indicating a similar mycobiome in Scots pine of different health classes at the already diseased (Diplodia tip blight) forest site. It appears that host species and tree genotype may influence the composition and community structure of the tree mycobiome [44,45].
The mycobiome of different tissues of asymptomatic and symptomatic pine trees was further analyzed. The lowest richness and community diversity of endophytic fungal communities was observed in the current-year shoots of symptomatic trees (DB). Statistically significant differences were exclusively observed in richness indices among asymptomatic and symptomatic trees in current-year shoots, and there are significant differences in the fungal communities among asymptomatic and symptomatic trees in the needles and shoots. Furthermore, there were no statistically significant differences in richness and community diversity among the three tissues in asymptomatic trees (Figures S5 and S6). Similarly, there were no significant differences in richness among the three tissues were observed in symptomatic tree (Figure S7). However, statistically significant differences in community diversity were observed only between shoots and phloem in symptomatic trees (Figure S8). Furthermore, there are significant differences in community structure related to shoots among asymptomatic and symptomatic trees. These results suggest that the structure and composition of the endophytic fungal community in the specific tissues where typical disease symptoms occur, such as current-year shoots, undergo significant changes. Similar to the study by Kovalchuk et al. [20], they indicated that Heterobasidion infection affects the fungal communities in the parts of the tree adjacent to the tissues colonized by the pathogen, with no significant effect on more distant regions.

4.3. Antagonism Assay In Vitro and In Vivo

A previous study demonstrated that certain endophytic fungi (A.alternata, Preussia sp., and Sy.polyspora) inhibited the growth of D. sapinea in antagonism tests [18]. However, in both our study and in research conducted by Bußkamp [3], Sy. polyspora did not exhibit inhibition when cultured alongside D. sapinea. This suggests that the observed antagonistic effects of fungal endophytes might be strain specific [8].
In our study, three strains (Pe. fructuariae-cellae, Taifanglania sp., and Pa. hawaiiense) demonstrated inhibitory effects on D. sapinea. These strains inhibited D. sapinea by over 60%, with Pe. fructuariae-cellae exhibiting the highest inhibition rate (73%). Lorenzini et al. [46] discovered that a novel Penicillium species from Italy, Pe. fructuariae-cellae, which infects grapes but with lower infectivity compared to Botrytis cinerea, the primary pathogen. In our study, Pe. fructuariae-cellae was assessed as an antagonistic endophytic fungus. Previous research has identified fungi belonging to the genus Penicillium as antagonistic endophytes in studies involving other plants. For example, Penicillium ehrlichii was isolated from Camellia sinensis roots as an antagonistic endophyte against two pathogens, Colletotrichum gloeosporioides and Pseudopestalotiopsis camelliae [47]. Our findings in vitro and in vivo suggest that Pe. fructuariae-cellae holds significant potential for suppressing fungal pathogens of D. sapinea. Therefore, future research should focus on extracting and identifying active compounds, as discovering and applying antibacterial agents could aid in the biocontrol of Diplodia tip blight.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jof10030212/s1. Figure S1: ANOSIM analysis of endophytic fungal communities from symptomatic and asymptomatic trees. Figure S2: ANOSIM analysis of endophytic fungal communities in A (needles), B (shoots), and C (phloem) samples. Figure S3: Box plots of ACE and Chao1 indices for samples of symptomatic and asymptomatic pine. Figure S4: Boxplots of Simpson and Shannon indices for samples of symptomatic and asymptomatic pine. Figure S5: Box plots of ACE and Chao1 indices for needles, shoots and phloem samples of asymptomatic pine. Figure S6: Box plots of Simpson and Shannon indices for needles, shoots and phloem samples of asymptomatic pine. Figure S7: Box plots of ACE and Chao1 indices for needles, shoots and phloem samples of symptomatic pine. Figure S8: Box plots of Simpson and Shannon indices for needles, shoots and phloem samples of symptomatic pine.

Author Contributions

C.Y. and L.W. conceived and designed the experiments; N.R. performed the experiments and wrote the manuscript; C.Y. revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Natural National Key Research and Development Program of China (2021YFC2600400).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

We would like to thank Yuxiang Liu from Beijing Forestry University for joining us to collect samples in Honghuaerji, Hulunbuir City, Inner Mongolia Autonomous Region.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kang, H.Z.; Zhu, J.J.; Xu, M.L. Study on water physiological properties of the artificially-planned saplings of Pinus sylvestris var. mongolica in the Horqin sand land. Arid Zone Res. 2007, 24, 15–22. [Google Scholar] [CrossRef]
  2. Xiang, C.T.; Yuan, S.Z.; Meng, F.R.; Jiang, J.Q. Study on diplodia tip blight of Pinus sylvestris var. mongolica. J. Northeast. For. Coll. 1981, 1–10, 129. [Google Scholar]
  3. Bußkamp, J. Schadenserhebung, Kartierung und Charakterisierung des “Diplodia-Triebsterbens“ der Kiefer, Insbesondere des Endophytischen Vorkommens in den Klimasensiblen Räumen und Identifikation von den in Kiefer (Pinus sylvestris) Vorkommenden Endophyten. [Damage Assessment, Mapping and Characterization of “Diplodia Tip Blight” of Pines, Especially the Endophytic Presence in Climate-Sensitive Areas and Identification of Pine Endophytes]. Ph.D. Thesis, Universität Kassel, Kassel, Germany, 2018. [Google Scholar]
  4. Bußkamp, J.; Langer, G.J.; Langer, E.J. Diplodia sapinea and fungal endophyte diversity in twigs of Scots pine (Pinus sylvestris) in Germany. Mycol. Progress. 2020, 19, 985–999. [Google Scholar] [CrossRef]
  5. Brodde, L.; Adamson, K.; Julio Camarero, J.; Castaño, C.; Drenkhan, R.; Lehtijärvi, A.; Luchi, N.; Migliorini, D.; Sánchez-Miranda, Á.; Stenlid, J. Diplodia tip blight on its way to the north: Drivers of disease emergence in northern Europe. Front. Plant Sci. 2019, 9, 414144. [Google Scholar] [CrossRef] [PubMed]
  6. Blumenstein, K.; Bußkamp, J.; Langer, G.J.; Langer, E.J.; Terhonen, E. The diplodia tip blight pathogen Sphaeropsis sapinea is the most common fungus in Scots pines’ mycobiome, irrespective of health status—A case study from Germany. J. Fungi 2021, 7, 607. [Google Scholar] [CrossRef]
  7. Oliva, J.; Boberg, J.; Stenlid, J. First Report of Sphaeropsis sapinea on Scots Pine (Pinus sylvestris) and Austrian Pine (P. nigra) in Sweden. New Dis. Rep. 2013, 27, 23. [Google Scholar] [CrossRef]
  8. Oliva, J.; Ridley, M.; Redondo, M.A.; Caballol, M. Competitive exclusion amongst endophytes determines shoot blight severity on pine. Funct. Ecol. 2021, 35, 239–254. [Google Scholar] [CrossRef]
  9. Terhonen, E.; Babalola, J.; Kasanen, R.; Jalkanen, R.; Blumenstein, K. Diplodia sapinea found as symptomless endophyte in Finland. Silva Fenn. 2021, 55, 10420. [Google Scholar] [CrossRef]
  10. Phillips, A.J.; Alves, A.; Abdollahzadeh, J.; Slippers, B.; Wingfield, M.J.; Groenewald, J.Z.; Crous, P.W. The Botryosphaeriaceae: Genera and species known from culture. Stud. Mycol. 2013, 76, 51–167. [Google Scholar] [CrossRef] [PubMed]
  11. Saikkonen, K.; Faeth, S.H.; Helander, M.; Sullivan, T.J. Fungal endophytes: A continuum of interactions with host plants. Ann. Rev. Ecol. Syst. 1998, 29, 319–343. [Google Scholar] [CrossRef]
  12. Alam, B.; Lǐ, J.; Gě, Q.; Khan, M.A.; Gōng, J.; Mehmood, S.; Yuan, Y.; Gǒng, W. Endophytic fungi: From symbiosis to secondary metabolitecommunications or vice versa? Front. Plant Sci. 2021, 12, 791033. [Google Scholar] [CrossRef] [PubMed]
  13. Vandenkoornhuyse, P.; Quaiser, A.; Duhamel, M.; Le Van, A.; Dufresne, A. The importance of the microbiome of the plant holobiont. New Phytologist. 2015, 206, 1196–1206. [Google Scholar] [CrossRef] [PubMed]
  14. Fabre, B.; Piou, D.; Desprez-Loustau, M.L.; Marçais, B. Can the emergence of pine Diplodia shoot blight in France be explained by changes in pathogen pressure linked to climate change? Glob. Chang. Biol. 2011, 17, 3218–3227. [Google Scholar] [CrossRef]
  15. Bosso, L.; Luchi, N.; Maresi, G.; Cristinzio, G.; Smeraldo, S.; Russo, D. Predicting current and future disease outbreaks of Diplodia sapinea shoot blight in Italy: Species distribution models as a tool for forest management planning. For. Ecol. Manag. 2017, 400, 655–664. [Google Scholar] [CrossRef]
  16. Slippers, B.; Wingfield, M.J. Botryosphaeriaceae as endophytes and latent pathogens of woody plants: Diversity, ecology and impact. Fungal Biol. Rev. 2007, 21, 90–106. [Google Scholar] [CrossRef]
  17. Sherwood, P.; Villari, C.; Capretti, P.; Bonello, P. Mechanisms of induced susceptibility to Diplodia tip blight in drought-stressed Austrian pine. Tree Physiol. 2015, 35, 549–562. [Google Scholar] [CrossRef] [PubMed]
  18. Blumenstein, K.; Bußkamp, J.; Langer, G.J.; Schlößer, R.; Parra Rojas, N.M.; Terhonen, E. Diplodia sapinea and associated endophytes in scots pine: Interactions and effect on the host under variable water content. Front. For. Glob. Chang. 2021, 4, 655769. [Google Scholar] [CrossRef]
  19. Peršoh, D. Factors shaping community structure of endophytic fungi—Evidence from the Pinus–Viscum-system. Fungal Divers. 2013, 60, 55–69. [Google Scholar] [CrossRef]
  20. Kovalchuk, A.; Mukrimin, M.; Zeng, Z.; Raffaello, T.; Liu, M.X.; Kasanen, R. Mycobiome analysis of asymptomatic and symptomatic Norway spruce trees naturally infected by the conifer pathogens Heterobasidion spp. Environ. Microbiol. Rep. 2018, 10, 532–541. [Google Scholar] [CrossRef]
  21. Chaudhary, P.; Chaudhary, A.; Bhatt, P.; Kumar, G.; Khatoon, H.; Rani, A.; Kumar, S.; Sharma, A. Assessment of soil health indicators under the influence of nanocompounds and Bacillus spp. in Field Condition. Front. Environ. Sci. 2022, 9, 769871. [Google Scholar] [CrossRef]
  22. Verma, S.K.; Sahu, P.K.; Kumar, K.; Pal, G.; Gond, S.K.; Kharwar, R.N.; White, J.F. Endophyte roles in nutrient acquisition, root system architecture development and oxidative stress tolerance. J. Appl. Microbiol. 2021, 131, 2161–2177. [Google Scholar] [CrossRef]
  23. Mejía, L.C.; Rojas, E.I.; Maynard, Z.; Van Bael, S.; Arnold, A.E.; Hebbar, P.; Samuels, G.J.; Robbins, N.; Herre, E.A. Endophytic fungi as biocontrol agents of Theobroma cacao pathogens. Biol. Control 2018, 46, 4–14. [Google Scholar] [CrossRef]
  24. Blumenstein, K.; Albrectsen, B.R.; Martín, J.A.; Hultberg, M.; Sieber, T.N.; Helander, M.; Witzell, J. Nutritional niche overlap potentiates the use of endophytes in biocontrol of a tree disease. BioControl 2015, 60, 655–667. [Google Scholar] [CrossRef]
  25. Pautasso, M.; Schlegel, M.; Holdenrieder, O. Forest health in a changing world. Microb. Ecol. 2015, 69, 826–842. [Google Scholar] [CrossRef] [PubMed]
  26. Witzell, J.; Martín, J.A.; Blumenstein, K. Ecological aspects of endophyte-based biocontrol of forest diseases. In Advances in Endophytic Research; Verma, V.C., Gange, A.C., Eds.; Springer: New Delhi, India, 2014; pp. 321–333. [Google Scholar] [CrossRef]
  27. Terhonen, E.; Kovalchuk, A.; Zarsav, A.; Asiegbu, F.O. Biocontrol potential of forest tree endophytes. In Endophytes of Forest Trees: Biology and Applications, Forestry Sciences; Pirttilä, A.M., Frank, A.C., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 283–318. [Google Scholar] [CrossRef]
  28. McMullin, D.R.; Nguyen, H.D.T.; Daly, G.J.; Menard, B.S.; Miller, J.D. Detection of foliar endophytes and their metabolites in Picea and Pinus seedling needles. Fungal Ecol. 2018, 31, 1–8. [Google Scholar] [CrossRef]
  29. Doyle, J.J. Isolation of plant DNA from fresh tissue. Focus 1990, 12, 13–15. [Google Scholar]
  30. White, T.J.; Bruns, T.; Lee, S.; Taylor, J. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In PCR Protocols: A Guide to Methods and Applications; Innis, M.A., Gelfand, D.H., Sninsky, J.J., White, T.J., Eds.; Academic Press: San Diego, CA, USA, 1990; Volume 18, pp. 315–322. [Google Scholar]
  31. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.; Holmes, S.P. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef]
  32. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F.; et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. [Google Scholar] [CrossRef]
  33. Grice, E.A.; Kong, H.H.; Conlan, S.; Deming, C.B.; Davis, J.; Young, A.C.; Bouffard, G.G.; Blakesley, R.W.; Murray, P.R.; Green, E.D.; et al. Topographical and temporal diversity of the human skin microbiome. Science 2009, 324, 1190–1192. [Google Scholar] [CrossRef]
  34. Blumenstein, K. Endophytic Fungi in Elms: Implications for the Integrated Management of Dutch Elm Disease; Bangor University: Bangor, UK, 2015. [Google Scholar]
  35. Wu, H.; Liu, W.; Liu, L.Y.; Ren, Y.L. Isolation and identification of endophytic fungi in junggar hawthorn fruit and study of their antibacterial metabolites. Food Res. Dev. 2021, 42, 143–148. [Google Scholar] [CrossRef]
  36. Decourcelle, T.; Piou, D.; Desprez-Loustau, M.-L. Detection of Diplodia sapinea in Corsican pine seeds. Plant Pathol. 2015, 64, 442–449. [Google Scholar] [CrossRef]
  37. Lindahl, B.D.; Nilsson, R.H.; Tedersoo, L.; Abarenkov, K.; Carlsen, T.; Kjøller, R.; Kõljalg, U.; Pennanen, T.; Rosendahl, S.; Stenlid, J.; et al. Fungal community analysis by high-throughput sequencing of amplified markers—A user’s guide. New Phytologist. 2013, 199, 288–299. [Google Scholar] [CrossRef] [PubMed]
  38. Basili, M.; Perini, L.; Zaggia, L.; Luna, G.M.; Quero, G.M. Integrating culture-based and molecular methods provides an improved assessment of microbial quality in a coastal lagoon. Environ. Pollution. 2023, 334, 122140. [Google Scholar] [CrossRef]
  39. Langer, G.; Bressem, U.; Habermann, M. Diplodia-Triebsterben der Kiefer und endophytischer Nachweis des Erregers Sphaeropsis sapinea. [Diplodia Tip Blight of pines and the endophytic proof of the pathogen Sphaeropsis sapinea]. AFZ-Der Wald 2011, 11, 28–31. [Google Scholar]
  40. Luchi, N.; Oliveira Longa, C.M.; Danti, R.; Capretti, P.; Maresi, G. Diplodia sapinea: The main fungal species involved in the colonization of pine shoots in Italy. For. Pathol. 2014, 44, 372–381. [Google Scholar] [CrossRef]
  41. Terhonen, E.; Marco, T.; Sun, H.; Jalkanen, R.; Kasanen, R.; Vuorinen, M.; Asiegbu, F. The effect of latitude, season and needle-age on the mycota of Scots pine (Pinus sylvestris) in Finland. Silva Fenn. 2011, 45, 301–317. [Google Scholar] [CrossRef]
  42. Silva, A.C.; Henriques, J.; Diogo, E.; Ramos, A.P.; Bragança, H. First report of Sydowia polyspora causing disease on Pinus pinea shoots. For. Pathol. 2020, 50, e12570. [Google Scholar] [CrossRef]
  43. Pan, Y.; Ye, H.; Lu, J.; Chen, P.; Zhou, X.D.; Qiao, M.; Yu, Z.F. Isolation and identification of Sydowia polyspora and its pathogenicity on Pinus yunnanensis in Southwestern China. J. Phytopathol. 2018, 166, 386–395. [Google Scholar] [CrossRef]
  44. Leopold, D.R.; Busby, P.E. Host genotype and colonist arrival order jointly govern plant microbiome composition and function. Curr. Biol. 2020, 30, 3260–3266.e3265. [Google Scholar] [CrossRef] [PubMed]
  45. Molina, L.; Rajchenberg, M.R.; Errasti, A.; Voge, B.; Coetzee, M.P.A.; Aime, M.C.; Pildain, M.B. Sapwood mycobiome vary across host, plant compartment and environments in Nothofagus forests from Northern Patagonia. Mol. Ecol. 2023, 32, 6599–6618. [Google Scholar] [CrossRef] [PubMed]
  46. Lorenzini, M.; Cappello, M.S.; Perrone, G.; Logrieco, A.; Zapparoli, G. New records of Penicillium and Aspergillus from withered grapes in Italy, and description of Penicillium fructuariae-cellae sp. nov. Phytopathol. Mediterr. 2019, 58, 323–340. [Google Scholar]
  47. Li, J.X. Screening and Bioactivity Evaluation of Antagonistic Endophytic Fungi from Camellia sinensis. Master’s Thesis, Anhui Agricultural University, Hefei, China, 2022. [Google Scholar]
Figure 1. Sampled Mongolian pine trees from the two sampling sites. (A) Asymptomatic Mongolian pine from Baogentu Forest Farm; (B) symptomatic Mongolian pine from Toudaoqiao Forest Farm; (C) schematic diagram illustrating the Mongolian pine tree tissues sampled for analyzing the associated endophytic fungal community.
Figure 1. Sampled Mongolian pine trees from the two sampling sites. (A) Asymptomatic Mongolian pine from Baogentu Forest Farm; (B) symptomatic Mongolian pine from Toudaoqiao Forest Farm; (C) schematic diagram illustrating the Mongolian pine tree tissues sampled for analyzing the associated endophytic fungal community.
Jof 10 00212 g001
Figure 2. Depiction of the four different reactions during the antagonism assays. (A) Inhibition of D. sapinea growth; (B) endophyte superior; (C) neutral or “mutually intermixing growth”; (D) D. sapinea superior.
Figure 2. Depiction of the four different reactions during the antagonism assays. (A) Inhibition of D. sapinea growth; (B) endophyte superior; (C) neutral or “mutually intermixing growth”; (D) D. sapinea superior.
Jof 10 00212 g002
Figure 3. Schematic diagram of a co-culture experiment with preferential placement of endophytic fungi.
Figure 3. Schematic diagram of a co-culture experiment with preferential placement of endophytic fungi.
Jof 10 00212 g003
Figure 4. Morphological characteristics of endophytic fungi isolated from different tissues of Mongolian pine (25 °C, seven days in the dark). (A) Alternaria alternata; (B) Alternaria infectoria; (C) Alternaria solani; (D) Alternaria tenuissima; (E) Aspergillus flavus; (F) Aspergillus sp.; (G) Botryotrichum murorum; (H) Chaetomium globosum; (I) Cladosporium sp.; (J) Collariella bostrychodes; (K) Coprinopsis atramentaria; (L) Diplodia sapinea; (M) Dothiorella gregaria; (N) Fusarium fujikuroi; (O) Fusarium sp.; (P) Paracamarosporium hawaiiense; (Q) Penicillium fructuariae-cella; (R) Phoma sp.; (S) Sarocladium zeae; (T)Sydowia polyspora.
Figure 4. Morphological characteristics of endophytic fungi isolated from different tissues of Mongolian pine (25 °C, seven days in the dark). (A) Alternaria alternata; (B) Alternaria infectoria; (C) Alternaria solani; (D) Alternaria tenuissima; (E) Aspergillus flavus; (F) Aspergillus sp.; (G) Botryotrichum murorum; (H) Chaetomium globosum; (I) Cladosporium sp.; (J) Collariella bostrychodes; (K) Coprinopsis atramentaria; (L) Diplodia sapinea; (M) Dothiorella gregaria; (N) Fusarium fujikuroi; (O) Fusarium sp.; (P) Paracamarosporium hawaiiense; (Q) Penicillium fructuariae-cella; (R) Phoma sp.; (S) Sarocladium zeae; (T)Sydowia polyspora.
Jof 10 00212 g004
Figure 5. Rarefaction curves for all analyzed samples.
Figure 5. Rarefaction curves for all analyzed samples.
Jof 10 00212 g005
Figure 6. Venn diagram displaying the ASVs shared by different tissues of asymptomatic and symptomatic pines.
Figure 6. Venn diagram displaying the ASVs shared by different tissues of asymptomatic and symptomatic pines.
Jof 10 00212 g006
Figure 7. Relative abundances of the endophytic fungal community at phylum (A) and genus (B) levels in needles, shoots, and phloem of Mongolian pine; relative abundances of the endophytic fungal community at genus (C) levels in symptomatic and asymptomatic trees. (H: samples from asymptomatic pine trees, and D: samples from symptomatic trees).
Figure 7. Relative abundances of the endophytic fungal community at phylum (A) and genus (B) levels in needles, shoots, and phloem of Mongolian pine; relative abundances of the endophytic fungal community at genus (C) levels in symptomatic and asymptomatic trees. (H: samples from asymptomatic pine trees, and D: samples from symptomatic trees).
Jof 10 00212 g007
Figure 8. Boxplots of ACE and Chao1 indices (** p < 0.01; ‘ns’ represents no significant difference).
Figure 8. Boxplots of ACE and Chao1 indices (** p < 0.01; ‘ns’ represents no significant difference).
Jof 10 00212 g008
Figure 9. Boxplots of Simpson and Shannon indices (** p < 0.01; *** p < 0.001; ‘ns’ represents no significant difference).
Figure 9. Boxplots of Simpson and Shannon indices (** p < 0.01; *** p < 0.001; ‘ns’ represents no significant difference).
Jof 10 00212 g009
Figure 10. PCoA based on the relative abundance of fungal ASVs, showing the differences in fungal community structure in sampled tissues.
Figure 10. PCoA based on the relative abundance of fungal ASVs, showing the differences in fungal community structure in sampled tissues.
Jof 10 00212 g010
Figure 11. LEfSe rank plots of shoots among asymptomatic and symptomatic Mongolian Pine.
Figure 11. LEfSe rank plots of shoots among asymptomatic and symptomatic Mongolian Pine.
Jof 10 00212 g011
Figure 12. FUNGuild analysis.
Figure 12. FUNGuild analysis.
Jof 10 00212 g012
Figure 13. FUNGuild analysis of ASVs from the three tissues of Mongolian pine. The top ten fungi with the highest percentage have been shown.
Figure 13. FUNGuild analysis of ASVs from the three tissues of Mongolian pine. The top ten fungi with the highest percentage have been shown.
Jof 10 00212 g013
Figure 14. Four different kinds of interaction between the endophytic fungi (on the right) and D. sapinea (on the left). (A) D. sapinea vs. Pe. fructuariae-cellae; (B) D. sapinea vs. Alternaria alternata; (C) D. sapinea vs. Trichoderma sp.; (D) D. sapinea vs. Phoma sp.
Figure 14. Four different kinds of interaction between the endophytic fungi (on the right) and D. sapinea (on the left). (A) D. sapinea vs. Pe. fructuariae-cellae; (B) D. sapinea vs. Alternaria alternata; (C) D. sapinea vs. Trichoderma sp.; (D) D. sapinea vs. Phoma sp.
Jof 10 00212 g014
Figure 15. Dual culture of the selected endophytic fungi (on the right) with D. sapinea (on the left), showing variations in the diameter of the inhibition zone. (A) D. sapinea vs. Taifanglania sp.; (B) D. sapinea vs. Pa. hawaiiense; (C) D. sapinea vs. Pe. fructuariae-cellae.
Figure 15. Dual culture of the selected endophytic fungi (on the right) with D. sapinea (on the left), showing variations in the diameter of the inhibition zone. (A) D. sapinea vs. Taifanglania sp.; (B) D. sapinea vs. Pa. hawaiiense; (C) D. sapinea vs. Pe. fructuariae-cellae.
Jof 10 00212 g015
Figure 16. The inhibition rate of Pe. fructuariae-cellae against D. sapinea at different time periods. (A) Concurrently culturing of Pe. fructuariae-cellae and D. sapinea; (BE) preferential culturing of Pe. fructuariae-cellae for 12 h (B), 24 h (C), 36 h (D), and 48 h (E), followed by co-culturing with D. sapinea; (F) D. sapinea alone. Different lower case letter indicate significant differences (p < 0.05).
Figure 16. The inhibition rate of Pe. fructuariae-cellae against D. sapinea at different time periods. (A) Concurrently culturing of Pe. fructuariae-cellae and D. sapinea; (BE) preferential culturing of Pe. fructuariae-cellae for 12 h (B), 24 h (C), 36 h (D), and 48 h (E), followed by co-culturing with D. sapinea; (F) D. sapinea alone. Different lower case letter indicate significant differences (p < 0.05).
Jof 10 00212 g016
Figure 17. The infected Mongolian pine seedlings on the 7th day after inoculation (The red circle in group I shows the necrosis on the inoculated shoot infected with D. sapinea).
Figure 17. The infected Mongolian pine seedlings on the 7th day after inoculation (The red circle in group I shows the necrosis on the inoculated shoot infected with D. sapinea).
Jof 10 00212 g017
Table 1. Taxa isolated from asymptomatic and symptomatic Mongolian Pine with the culture-based method.
Table 1. Taxa isolated from asymptomatic and symptomatic Mongolian Pine with the culture-based method.
Isolation Strain No.StrainsGenbank
Accession Number
HAHBHCDADBDCTotal
NIFNIFNIFNIFNIFNIFNIF
EMD1Alternaria alternataOR002004.122.86 32.78 23.51 --22.33 510.42143.23
EMD2Alternaria infectoriaOR067375.1--21.85 11.75 ------30.69
EMD3Alternaria seleniiphilaMK140693.1--------44.65 --40.92
EMD4Alternaria solaniMN871615.1--10.93 11.75 11.56 ----30.69
EMD5Alternaria sp. MN096578.1--------11.16 36.25 40.92
EMD6Alternaria tenuissimaKM921667.1--10.93 --11.56 --12.08 30.69
EMD7Aspergillus flavusMH931826.1--10.93 --57.81 44.65 24.17 122.77
EMD8Aspergillus sp. KR154911.1--10.93 11.75 ------20.46
EMD9Botryotrichum murorumMG770259.1--10.93 --------10.23
EMD10Chaetomium globosumMG885806.111.43 ----------10.23
EMD11Chromolaenicola clematidisMT310601.1--10.93 --------10.23
EMD12Cladosporium cladosporioidesOQ555096.1------34.69 --12.08 40.92
EMD13Cladosporium perangustumMK722299.1----11.75 ------10.23
EMD14Cladosporium sp.MG975642.11014.291715.74 1424.56 57.81 66.98 36.25 5512.70
EMD15Collariella bostrychodesMH931826.1--10.93 --------10.23
EMD16Coprinopsis atramentariaKJ817302.1--10.93 11.75 --11.16 --30.69
EMD17Curvularia senegalensisMT476857.1--------22.33 --20.46
EMD18Deniquelata sp.ON705536.1--------11.16 --10.23
EMD19Diplodia sapineaMT763348.1710.00 1513.89 58.77 1015.63 2023.26 816.67 6515.01
EMD20Dothiorella gregariaMH791151.1--10.93 --------10.23
EMD21Fusarium fujikuroiMH084746.1--21.85 --------20.46
EMD22Fusarium sp.OQ422002.122.86 ----------20.46
EMD23Fusarium verticillioidesOM956059.1--43.70 --------40.92
EMD24Gregarithecium sp. DQD-2016aKX364281.1--21.85 11.75 --33.49 --61.39
EMD25Iodophanus carneusMF161095.1--10.93 --------10.23
EMD26Microsphaeropsis sp. MN153956.1--------11.16 --10.23
EMD27Neocamarosporium salicorniicolaMK809918.1------11.56 ----10.23
EMD28Nothophoma quercinaON527450.111.43 ----------10.23
EMD29Paracamarosporium hawaiienseHM751092.1--10.93 --------10.23
EMD30Paraphaeosphaeria sp. GU985234.1------23.13 ----20.46
EMD31Penicillium sp. MK817609.145.71 54.63 1322.81 1015.63 1112.79 24.17 4510.39
EMD32Penicillium fructuariae-cellaeOP703391.1--10.93 23.51 --22.33 24.17 71.62
EMD33Penicillium glabrumOP681429.111.43 ----------10.23
EMD34Penicillium oxalicumMT597864.1--------11.16 --10.23
EMD35Periconia macrospinosaMH345963.111.43 ----------10.23
EMD36Phialophora cyclaminisMW447030.1--------11.16 --10.23
EMD37Phoma sp.MT366774.1------23.13 ----20.46
EMD38Podospora australisKX015765.1------34.69 ----30.69
EMD39Podospora sp.MW349896.1----------12.08 10.23
EMD40Preussia africanaOM743863.1--------11.16 --10.23
EMD41Preussia persicaMW081300.1--10.93 --------10.23
EMD42Preussia sp.MT862280.1--76.48 --------71.62
EMD43Rhizopus arrhizusON920727.1----------36.25 30.69
EMD44Sarocladium zeaeMZ969671.1--10.93 --------10.23
EMD45Subramaniula cuniculorumMT072073.1----11.75 --22.33 --30.69
EMD46Sydowia polysporaMN900630.13245.71 2523.15 1017.54 812.50 1416.28 510.42 9421.71
EMD47Taifanglania parvisporaKF719170.1--------44.65 --40.92
EMD48Taifanglania sp.KP143100.1------34.69 ----30.69
EMD49Talaromyces amestolkiaeKT445914.134.29 --11.75 34.69 --12.08 81.85
EMD50Talaromyces cecidicolaMN889416.1--10.93 --------10.23
EMD51Talaromyces funiculosusMT367866.111.43 ----------10.23
EMD52Talaromyces pseudofuniculosusOP482383.134.29 98.33 35.26 46.25 33.49 918.75 317.16
EMD53Torula ficiOW988181.1----------24.17 20.46
EMD54Tricharina sp.MF055701.111.43 ----------10.23
EMD55Trichoderma afroharzianumMT102402.1--10.93 --------10.23
EMD56Trichoderma harzianumKY495199.1--10.93 --34.69 ----40.92
EMD57Trichoderma sp.MK871263.111.43 ------22.33 --30.69
Total 70100.00 108100.00 57100.00 64100.00 86100.00 48100.00 433100.00
“-” means the corresponding endophytic fungal species was not isolated from the specific tissue; “No.” represents the occurrence of certain endophytic fungi in different tissues; the isolation frequency: IF (%) = No. of Isolates/Total No. of Isolates (%).
Table 2. Pair_adonis analysis of endophytic fungal community structure at different tissues.
Table 2. Pair_adonis analysis of endophytic fungal community structure at different tissues.
PairsR2p-Value
HA vs. DA0.160197190.005
HB vs. DB0.406956180.001
HC vs. DC0.078956050.115
Table 3. Observations of All Isolated endophytic fungi against D. sapinea.
Table 3. Observations of All Isolated endophytic fungi against D. sapinea.
Isolation Strain No.StrainsVisual Observation
EMD1Alternaria alternataPathogen superior
EMD2Alternaria infectoriaPathogen superior
EMD3Alternaria seleniiphilaPathogen superior
EMD4Alternaria solaniPathogen superior
EMD5Alternaria sp.Pathogen superior
EMD6Alternaria tenuissimaPathogen superior
EMD7Aspergillus flavusEqual growth capability
EMD8Aspergillus sp.Endophyte superior
EMD9Botryotrichum murorumPathogen superior
EMD10Chaetomium globosumPathogen superior
EMD11Chromolaenicola clematidisPathogen superior
EMD12Cladosporium cladosporioidesPathogen superior
EMD13Cladosporium perangustumPathogen superior
EMD14Cladosporium sp.Pathogen superior
EMD15Collariella bostrychodesPathogen superior
EMD16Coprinopsis atramentariaPathogen superior
EMD17Curvularia senegalensisPathogen superior
EMD18Deniquelata sp.Pathogen superior
EMD20Dothiorella gregariaPathogen superior
EMD21Fusarium fujikuroiPathogen superior
EMD22Fusarium sp.Pathogen superior
EMD23Fusarium verticillioidesPathogen superior
EMD24Gregarithecium sp. DQD-2016aPathogen superior
EMD25Iodophanus carneusPathogen superior
EMD26Microsphaeropsis sp.Pathogen superior
EMD27Neocamarosporium salicorniicolaPathogen superior
EMD28Nothophoma quercinaPathogen superior
EMD29Paracamarosporium hawaiienseInhibition of pathogen growth
EMD30Paraphaeosphaeria sp.Pathogen superior
EMD31Penicillium sp.Pathogen superior
EMD32Penicillium fructuariae-cellaeInhibition of pathogen growth
EMD33Penicillium glabrumPathogen superior
EMD34Penicillium oxalicumPathogen superior
EMD35Periconia macrospinosaPathogen superior
EMD36Phialophora cyclaminisPathogen superior
EMD37Phoma sp.Equal growth capability
EMD38Podospora australisPathogen superior
EMD39Podospora sp.Pathogen superior
EMD40Preussia africanaPathogen superior
EMD41Preussia persicaPathogen superior
EMD42Preussia sp.Pathogen superior
EMD43Rhizopus arrhizusPathogen superior
EMD44Sarocladium zeaePathogen superior
EMD45Subramaniula cuniculorumPathogen superior
EMD46Sydowia polysporaPathogen superior
EMD47Taifanglania parvisporaPathogen superior
EMD48Taifanglania sp.Inhibition of pathogen growth
EMD49Talaromyces amestolkiaePathogen superior
EMD50Talaromyces cecidicolaPathogen superior
EMD51Talaromyces funiculosusPathogen superior
EMD52Talaromyces pseudofuniculosusPathogen superior
EMD53Torula ficiPathogen superior
EMD54Tricharina sp.Endophyte superior
EMD55Trichoderma afroharzianumEndophyte superior
EMD56Trichoderma harzianumEndophyte superior
EMD57Trichoderma sp.Pathogen superior
Table 4. The Inhibition Rate of the Selected Endophytic Fungi against D. sapinea after 7 days.
Table 4. The Inhibition Rate of the Selected Endophytic Fungi against D. sapinea after 7 days.
Isolation Strain No.StrainsInhibition Rate (%)Isolated Tissue
EMD48Taifanglania sp.(61.00 ± 1.73) bDA
EMD29Paracamarosporium hawaiiense(72.00. ± 4.00) aHB
EMD32Penicillium fructuariae-cellae(73.00. ± 4.62) aHB and HC
Different letters indicate significant differences (p < 0.05).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ren, N.; Wang, L.; You, C. Diversity, Community Structure, and Antagonism of Endophytic Fungi from Asymptomatic and Symptomatic Mongolian Pine Trees. J. Fungi 2024, 10, 212. https://doi.org/10.3390/jof10030212

AMA Style

Ren N, Wang L, You C. Diversity, Community Structure, and Antagonism of Endophytic Fungi from Asymptomatic and Symptomatic Mongolian Pine Trees. Journal of Fungi. 2024; 10(3):212. https://doi.org/10.3390/jof10030212

Chicago/Turabian Style

Ren, Ninghong, Lei Wang, and Chongjuan You. 2024. "Diversity, Community Structure, and Antagonism of Endophytic Fungi from Asymptomatic and Symptomatic Mongolian Pine Trees" Journal of Fungi 10, no. 3: 212. https://doi.org/10.3390/jof10030212

APA Style

Ren, N., Wang, L., & You, C. (2024). Diversity, Community Structure, and Antagonism of Endophytic Fungi from Asymptomatic and Symptomatic Mongolian Pine Trees. Journal of Fungi, 10(3), 212. https://doi.org/10.3390/jof10030212

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop