*3.1. Diversity Analysis of Fungal Microbiomes in MOR and AOF Samples*

A total of 1,413,703 valid ITS1 sequences with an average length of 235 bp were obtained from 18 samples. The sequences were clustered into 579 OTUs (≥97% similarity, Table S1). Venn analysis exhibited that 91 and 76 OTUs were, respectively, unique for AOF and MOR groups, and the remaining 412 OTUs were common in the two groups (Figure 1a). Moreover, 282 OTUs were shared by 4 groups, and 23, 34, 43, and 28 OTUs were, respectively, unique for NM, MM, NA, and MA groups (Figure 1b). With the number of sequences sampled increasing, the rarefaction curves of all samples were parallel to the *x*-axis, which indicates the reliability of the sequencing depth employed (Figure 1c). The 5 alpha-diversity indices were employed to analyze the richness, diversity, and coverage of fungal microbiomes in 18 samples (Table 2). The Good's coverage was over 99.8%, which shows that the sampling depth satisfies the analysis requirements. In the MOR group, the ACE and Chao 1 indices of GDM were the highest, and the community richness was the highest. In contrast, GXM had the lowest community richness. The community diversity of GDM samples was the highest with the highest Shannon and Simpson indices. GXM had the lowest community diversity. Similarly, HNA and ACK in the AOF group had the highest community diversity and community richness, respectively. PCoA analysis showed that AOF and MOR groups were distinguishable (Figure 1d).

**Figure 1.** Analysis of fungal diversity in MOR and AOF samples. (**a**) Venn diagram of OTUs in the MOR and AOF groups; (**b**) Venn diagram of OTUs in the moldy and non-moldy samples; (**c**) Rarefaction curves for OTUs in all samples; and (**d**) PCoA diagram. of fungal compositions in samples.


**Table 2.** Alpha-diversity indices of samples.
