Taxonomy group
Fig 4 illustrates the main fungal groups in healthy and Diseased sites across all sites. In Koto Rami Healthy (H) soil, the dominant fungi were Ascomycota and Basidiomycota and smaller groups like Mucoromycota and Zoopagomycota. The Ascomycota phylum includes Eurotiomycetes. When drawn into the family Aspergillaceae (Fig 4C), it was found in classes such as Tremellomycetes and notable genera like Aspergillus and Penicillium.
The quantitative distribution of Aspergillus and Penicillium is illustrated as follows (Fig 5). In Kumpeh H, out of 45.200 observed microbes, 22.400 were from the Didymosphaeriaceae family. The Trichocomaceae family had 4.110 members and the Aspergillaceae family included Aspergillus (2.620), Penicillium (3.530) and Talaromyces (4.080). In Koto Rami H, Clavariaceae led with 6.520 entities, followed by Aspergillus (3.690), Penicillium (2.510) and Ganoderma (741). Key species included
Penicillium flaviroseum (414) and
Aspergillus fumigatus (523). In Rano H, Aspergillus comprised 12.400 entities, while Penicillium had 31.500 entities, including
Aspergillus nomiae (2.950) and
Penicillium citrinum (5.900). In the Diseased site, Aspergillus counted 3.310 and Ganoderma 3.470, with various species identified.
Fungal relative abundance
Table 2 presents fungal relative abundance (%) at the phylum level under Diseased site conditions and three healthy locations: Koto Rami, Kumpeh and Rano. Ascomycota and Basidiomycota were the dominant phyla across all conditions, with Ascomycota abundance ranging from 0,4482 to 0,5972% and Basidiomycota from 0,3359 to 0,4849%. Kumpeh (H) had the highest Ascomycota abundance (0,5972%), while Koto Rami (H) had the highest Basidiomycota (0,4849%). Other phyla, such as Blastocla-diomycota, Chytridiomycota, Cryptomycota, Olpidiomycota and Zoopagomycota, were less abundant, each below 0.05%. Blastocladiomycota was absent in Rano (H), while Chytridiomycota was highest at Kumpeh (H) (0.0122%). Mucoromycota increased at Rano (H) (0,1616%) and Zoopagomycota was most abundant at Kumpeh (H) (0,0161%). The diseased site soil, Ascomycota (0.455%) and Basidiomycota (0.4776%), displays a more balanced ratio than the Healthy sites.
Fungal Alpha diversity
The soils sampled from Diseased sites demonstrated the highest observed species richness, quantified at 3,636. This significant level of diversity is further supported by high estimates of potential richness, with Chao1 calculated at 5,370.40 and ACE at 5,361.10. The Shannon Index for these samples was recorded at 39.60, marking it the highest among all sites. Simpson’s Index was measured at 5.01, indicating moderate species-level dominance within these communities. Additionally, the InvSimpson’s Index, at 0.96, suggests a less equitable species distribution compared to Kumpeh_H. Despite the Diseased site showcasing the highest observed species richness, it remains evident that its potential richness is lower relative to other studied areas.
Koto Rami documented 3105 species, surpassing Rano_H (2900) but still behind Diseased site (3636) and Kumpeh_H (3297), indicating a diverse ecosystem. Chao1 and ACE estimates suggested potential species richness of around 4693,01 and 4586,20, respectively. The Shannon index was 36,79, slightly lower than Kumpeh_H (38,35), while both the Simpson’s Index (5,38) and InvSimpson Index (0,97) indicated an even distribution of species. Koto Rami’s Fisher’s Alpha of 34,81 suggested a high potential for species richness. In contrast, Rano_H had the lowest species richness, 2900 species and potential richness estimated at 4482,63. Its Shannon Index of 35,74 and Simpson Dominance Index of 4,65 revealed an uneven distribution, dominated by certain species. Rano_H’s Fisher’s Alpha of 30,40 indicated limited diversity. Rano_H was the least diverse, with an imbalanced species distribution (Fig 6).
Site distinctive of fungal diversity
Each site had unique microbial communities, as shown in Fig 8. The Diseased site soil contains the highest number of exclusive fungi (805) (Fig 7 B) and has a similarity with Kumpeh H (296). Overall, 1064 microbes were identified. Kumpeh H had the highest variety with 1036 unique species, while Rano H and Koto Rami H had 741 and 980, respectively. Shared species included 409 between Rano H and Koto Rami H, 545 between Rano H and Kumpeh H and 511 between Koto Rami H and Kumpeh H, totaling 1205 common microbial phyla across all sites.
Beta diversity
Beta diversity describes the variation or differences between healthy and Diseased site of duku plantations. Fig 8 presents a principal coordinate analysis (PCoA) based on the Bray-Curtis distance, which evaluates the soil microbial community structure at Kumpeh H, Koto Rami H, Rano H and the Diseased site. The analysis showed distinct microbial distributions among the sites, with PCoA1 (40,16%) and PCoA2 (24,54%) on the left and PCoA1 and PCoA3 (25,21%) on the right.
This study compares soil fungal profiles in healthy and Diseased duku plantations at various altitudes in Jambi. It highlights how soil properties affect fungal abundance. The results show that differences in soil characteristics-such as pH, organic matter, nutrient availability and altitude-influence fungal communities across the study sites.
Rano H, 5.80 m above sea level, showed the highest abundance of Aspergillus and Penicillium. It also has elevated Phosphorus levels and a relatively neutral pH. This finding is in line with previous studies suggesting that Phosphorus availability boosts fungal growth, especially among Ascomycota species such as Aspergillus and Penicillium
(Wang et al., 2018). In contrast, Kumpeh H, located 10.32 m above sea level, had a moderately acidic pH and low organic matter content, leading to lower counts of Aspergillus and Penicillium, while Didymosphaeriaceae thrived instead. The thriving of Didymosphaeriaceae indicates that certain fungi can flourish in nutrient-poor soils in lowland areas, which may contribute to soil degradation and increased susceptibility to diseases (
Treseder and Lennon, 2015).
At an elevation of 623.46 m above sea level, Koto Rami H exhibited distinct soil characteristics, including the highest levels of organic matter and total Nitrogen. Despite low nutrient levels of Phosphorus and Kalium, Koto Rami H supported a diverse fungal community, with higher counts of Aspergillus and Penicillium compared to Rano H. Typically, as altitude increases, temperatures drop, which slows the decomposition of organic matter. This slowdown leads to a greater buildup of organic materials and, in turn, fosters a richer diversity of fungi
(Hartmann et al., 2015; Zahra et al., 2021). Additionally, Koto Rami H showed the highest concentrations of Iron and Zinc, which can enhance the enzymatic activities of fungi. However, excessive iron can negatively affect fungal communities
(Luo et al., 2024).
The presence of Aspergillus and Ganoderma at the Diseased site hints at a fungal-driven disease progression, as elevated microbial respiration and CEC are frequently linked to microbial competition and shifts in fungal community structure in degraded soils
(Zhou et al., 2023). The disease site’s intermediate altitude, due to its specific soil conditions, may have created an environment that favors fungal communities associated with plant disease.
At the Diseased site, high microbial activity but lower diversity was observed. It reflects stress conditions that can favor pathogenic fungi over beneficial mutualists
(Gupta et al., 2022). This finding resonates with studies indicating that Diseased plant soil often hosts more complex fungal networks than healthy soil
(Jia et al., 2022). The elevated CEC in Diseased soil implies improved nutrient cycling; however, this does not necessarily translate to enhanced plant health, as high microbial respiration may signify nutrient competition or stress (
Meimaroglou and Mouzakis, 2019;
Lalkhumliana et al., 2024).
Across all sites, Ascomycota and Basidiomycota were the dominant fungal groups, though variations were observed depending on environmental conditions
(Feng et al., 2024). The ratio of Ascomycota to Basidiomycota can serve as a soil health indicator, as a dominance of Ascomycota in organic-rich soils may limit the presence of beneficial fungal groups such as Basidiomycota and arbuscular mycorrhizal fungi
(Manici et al., 2024; Martínez-García et al., 2018). The lower Nitrogen content in Kumpeh H could help explain its reduced fungal abundance (
Di Lonardo et al., 2020). Higher Phosphorus levels in Rano H might support the growth of Mucoromycota, which is important for hyphal development
(Bhalla et al., 2022). Furthermore, the greater moisture levels in Koto Rami H and Kumpeh H may be linked to increasing Zoopagomycota, which thrive in damp conditions
(Chen et al., 2023).
Fungi play vital roles in soil health as decomposers, mutualists and pathogens
(Li et al., 2019; Yiallouris et al., 2024). The prevalence of Penicillium and Aspergillus across various sites contributes to organic matter decomposition, Phosphorus solubilization and microbial antagonism (
Díaz-Urbano et al., 2023). While Basidiomycota populations can stabilize microbial interactions, disruptions such as land use changes and agricultural practices may diminish their abundance, ultimately impacting soil health
(Egidi et al., 2019; Xiang et al., 2024). Moreover, fungal activity significantly influences soil structure by producing microbial exudates that improve soil aggregation
(Arias et al., 2023; Sudheer et al., 2024).
The metagenomic analysis conducted in this study successfully identified a diverse fungal community across all sampling sites, including healthy and diseased duku plantations. Dominant genera such as
Aspergillus,
Penicillium,
Ganoderma and members of the
Didymosp-haeriaceae and
Clavariaceae families were detected in varying abundance. However, it is notable that previously reported pathogens associated with duku diseases, such as
Phytophthora palmivora (
Handoko, 2014) and
Ceratocystis sp (
Firmanto, 2023), were not detected in the current dataset.
Several factors may account for this discrepancy. First, metagenomic sequencing is highly dependent on the DNA extraction method and sequencing depth. It is possible that the DNA of
Phytophthora or
Ceratocystis was present at very low concentrations, falling below the detection limit of the sequencing pipeline (
Sharpton, 2014). Additionally, the reference database used for taxonomic classification (ITS RefSeq) may have limitations in resolving certain oomycete taxa like
Phytophthora, which often require other markers such as
cox1 or
cox2 for accurate identification
(Robideau et al., 2011).
Second, the symptoms observed in the diseased site may not be caused solely by fungal or oomycete pathogens, but could be the result of complex biotic interactions, including latent infections, bacterial co-infections, or abiotic stressors that mimic pathogen-induced symptoms (
Scholthof, 2007). Without pathogen isolation and pathogenicity testing (fulfilling Koch’s postulates), it remains speculative to attribute symptoms to specific taxa based solely on sequencing data.
Finally, spatial and temporal variation in pathogen presence could have influenced the results. Pathogens like
P. palmivora may not be uniformly distributed in soil or plant tissues and may only be detectable during specific stages of infection (
Hardham and Blackman, 2018). The findings highlight the need to combine metagenomic approaches with traditional plant pathology methods like culturing, microscopy and pathogenicity tests for better disease diagnostics.