Research Article |
Corresponding author: Nakarin Suwannarach ( suwan.462@gmail.com ) Corresponding author: Hao-Han Wang ( wanghh@eastern-himalaya.cn ) Academic editor: Danushka Sandaruwan Tennakoon
© 2025 Xian Zhang, Saowaluck Tibpromma, Samantha C. Karunarathna, Tian-Ye Du, Li-Su Han, Abdallah M. Elgorban, Jaturong Kumla, Chanokned Senwanna, Dong-Qin Dai, Nakarin Suwannarach, Hao-Han Wang.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Zhang X, Tibpromma S, Karunarathna SC, Du T-Y, Han L-S, Elgorban AM, Kumla J, Senwanna C, Dai D-Q, Suwannarach N, Wang H-H (2025) Additions to the saprobic fungi (Ascomycota) associated with macadamia trees from the Greater Mekong Subregion. MycoKeys 113: 1-29. https://doi.org/10.3897/mycokeys.113.140031
|
Macadamia trees, the most economically important Proteaceae perennial crop, are globally renowned for their edible nuts. During our surveys of microfungi associated with macadamia in China and Thailand, we isolated three saprobic fungi from dead macadamia branches. Our multigene phylogenetic analyses (ITS, LSU, SSU, tef1-α, TUB2, and ACT loci), genealogical concordance phylogenetic species recognition (GCPSR) with a pairwise homoplasy index (PHI) test, and morphological characteristics led to the discovery of two new species, Dothiorella macadamiae and Phaeoacremonium chiangmaiense, and one new record, Melomastia puerensis. We provide morphological descriptions, photo plates, phylogenetic analysis results, and PHI test results of the two new species, along with comparisons with closely related taxa. These findings have global implications for understanding the diversity of microfungi associated with macadamia trees.
Dothideomycetes, morphology, multigene phylogeny, new taxa, saprobes, Sordariomycetes, taxonomy
Macadamia F. Muell. is one of the perennial crops that produces an edible nut, which is native to Australia and later widely distributed in frost-free tropical and subtropical regions (
Macadamia was initially introduced from Australia and trial-planted in China in the 1970s (
Microfungi are one of the key organisms in forest ecosystems, and they are distributed worldwide with a very high diversity (
In this study, we aim to introduce two new species and one new host record of saprobic fungi found on dead twigs of Macadamia spp. based on morphological characteristics, multi-locus phylogeny analyses of the combined ITS, LSU, SSU, tef1-α, TUB2, and ACT sequence data, and genealogical concordance phylogenetic species recognition (GCPSR) with a pairwise homoplasy index (PHI).
Dead twigs of Macadamia spp. with fungal fruiting bodies were collected from the Yunnan Province of China and the Chiang Mai Province of Thailand. Each sample was placed in a separate plastic bag, marked with information such as collection site, collection date, collection altitude, and global positioning system (GPS) (
DNA was extracted from the three-week-old pure cultures growing on PDA, or fungal fruiting bodies were used when the pure cultures could not be obtained. Genomic DNA extraction was obtained using the Biospin Fungus Genomic DNA Extraction Kit-BSC14M1 (BioFlux®, P.R. China), following the manufacturer’s guidelines and preserved at -20 °C for long-term use. The different gene regions, primers, and PCR thermal cycle programs for amplification are detailed in Table
Genes | Primers | PCR conditions | References |
---|---|---|---|
ITS | ITS5/ITS4 | 95 °C: 3 mins, (94 °C: 30 s, 55 °C: 50 s, 72 °C: 90 s) × 35 cycles, 72 °C: 10 mins (Phaeoacremonium) |
|
95 °C: 2 mins, (95 °C: 30 s, 52 °C: 30 s, 72 °C: 60 s) × 35 cycles, 72 °C: 10 mins (Dothiorella) | |||
LSU | LR0R/LR5 | 95 °C: 3 mins, (94 °C: 30 s, 55 °C: 50 s, 72 °C: 90 s) × 35 cycles, 72 °C: 10 mins (Phaeoacremonium, Melomastia) |
|
SSU | NS1/NS4 | 95 °C: 3 mins, (94 °C: 30 s, 55 °C: 50 s, 72 °C: 90 s) × 35 cycles, 72 °C: 10 mins (Melomastia) |
|
tef1-α | 983F/2218R | 95 °C: 3 mins, (95 °C: 30 s, 55 °C: 50 s, 72 °C: 90 s) × 35 cycles, 72 °C: 10 mins (Phaeoacremonium, Melomastia) |
|
728F/986R | 95 °C: 3 mins, (95 °C: 30 s, 55 °C: 30 s, 72 °C: 60 s) × 40 cycles, 72 °C: 10 mins (Dothiorella) |
|
|
TUB2 | T1/Bt2b | 94 °C: 4 mins, (94 °C: 40 s, 52 °C: 30 s, 72 °C: 60 s) × 35 cycles, 72 °C: 10 mins (Phaeoacremonium, Dothiorella) |
|
ACT | ACT-512F/ACT-783R | 94 °C: 4 mins, (94 °C: 40 s, 52 °C: 30 s, 72 °C: 60 s) × 35 cycles, 72 °C: 10 mins (Phaeoacremonium) |
|
Sequences were assembled by the Geneious program (9.1.2) (https://www.geneious.com/), and the newly generated assembled sequences were copied to BLASTn for searches. The related sequences used in phylogenetic analyses were retrieved from GenBank following BLASTn search results and the latest publication (
The phylogenetic analyses of combined genes (ITS, tef1-α, and TUB2 for the Dothiorella dataset; ITS+TUB2+ACT+tef1-α+LSU for the Phaeoacremonium dataset; and LSU+SSU+tef1-α for the Melomastia dataset) were based on ML and BI analyses. The ML trees were performed using RAxML-HPC2 on XSEDE (8.2.12) (
The pairwise homoplasy index (PHI) test using SplitsTree V4 (
Dothiorella (Botryosphaeriaceae Theiss. & Syd.) was introduced by
“macadamiae” refers to the host plant genus Macadamia.
CMUB 40066.
Saprobic on dead twigs of Macadamia sp. Sexual morph: Ascomata 80–120 × 160–220 µm (x̄ = 100 × 193 µm, n = 10), immersed, visible as dark dots on the host surface, under to clypeus, solitary, uni-loculate, ampulliform, papillate, without ostiole. Peridium 50–190 µm wide (x̄ = 95 µm, n = 25), comprising three section layers, the inner section layer composed of hyaline cells of textura angularis, the outer section layer with brown to dark brown cells of textura angularis, and the outermost layer of cells surrounding the ascomata is composed of brown cells of textura prismatica. Hamathecium 4.5–8.5 µm wide (x̄ = 6.5 μm, n = 30), comprising cylindrical, hyaline, septate, cellular pseudoparaphyses. Asci 110–235 × 23–38 µm (x̄ = 173 × 32 µm, n = 30), 6–8-spored, bitunicate, clavate to broadly fusoid-ellipsoid, with furcated pedicel, apically rounded, with a well-developed ocular chamber. Ascospores (27–)30–37(–40) × 14–19 µm (x̄ = 33 × 17 µm, n = 55), overlapping, uniseriate, oval to ellipsoid, hyaline to yellowish brown, aseptate when young, becoming brown to dark brown, 1- or 2- septate at maturity, slightly constricted at the septum, smooth-walled, granular, with mucilaginous polar appendages at one or both ends. Asexual morph: Conidiomata pycnidial produced on PDA within seven weeks, solitary or aggregated, superficial, brown, hairy, globose to subglobose, covered with hyphal hairs, unilocular. Conidiophores reduce to conidiogenous cells. Conidiogenous cells holoblastic, discrete, cylindrical, hyaline, smooth-walled, proliferating percurrently. Conidia 19–26.5 × 10–13.5 µm (x̄ = 22.6 × 11.2 µm, n = 50), hyaline and aseptate when immature, brown to dark brown and one-septate when mature, oblong to ovoid, granular, one end obtuse to slightly rounded ends, one cell slightly wider or same width. Chlamydospores hyaline to brown, branched, with thickened, septate, brown to dark brown at the septa.
Ascospores germinating on PDA within 24 h at 28 °C, colony on PDA reaching 9 cm diam. after two weeks at 28 °C, rough surface, hairy, cottony, and pale olivaceous grey from above, and grey to black in reverse.
Thailand • Chiang Mai Province, Doi Saket District, 18°52'43"N, 99°13'15"E, 384 m elevation, on a dead branch of Macadamia sp., 24 November 2023, Xian Zhang, TCMM25, CMUB 40066, holotype; ex-type living culture, SDBR-CMU512, other living culture SDBR-CMU513.
In the phylogenetic analyses, our isolate D. macadamiae forms an independent branch sister to D. albiziae and D. thailandica with 57% ML and 1.00 PP support (Fig.
RAxML tree based on a combined dataset of ITS, tef1-α, and TUB2 gene sequences data, which comprised 1265 base pairs (ITS = 1–527 bp, tef1-α = 528–818 bp, TUB2 = 819–1265 bp). The best scoring RAxML tree with a final ML optimization likelihood value of -8702.022443 is presented. The matrix had 696 distinct alignment patterns, with 16.83% of undetermined characters or gaps. Estimated base frequencies were as follows: A = 0.209011, C = 0.306672, G = 0.250984, T = 0.233333; substitution rates: AC = 1.037189, AG = 2.254040, AT = 1.059283, CG = 1.045839, CT = 4.497195, GT = 1.000000; proportion of invariable sites I = 0.499857; gamma distribution shape parameter α = 0.610523. The ML analysis and Bayesian inference (BI) analyses showed nearly identical tree topologies, bootstrap support values for ML equal to or greater than 60%, and BI analysis values equal to or greater than 0.90 PP are given at each node. The tree is rooted with Neofusicoccum luteum (CBS 562.92) and N. luteum (CMW 41365). Newly generated species are shown in red, while the ex-type species are shown in bold. Remarks: The ML bootstrap value less than 60% is presented on the node of the new taxon.
Results of the PHI test of Dothiorella macadamiae and closely related species using both LogDet transformation and splits decomposition. The PHI test results (Φw) < 0.05 indicate significant recombination within the dataset. The new taxa are in red font, and bold indicates holotype or ex-type strains.
Dothiorella macadamiae (CMUB 40066, holotype) a–c appearance of ascomata on the host surface d vertical section of an ascoma e peridium f hamathecium g–i asci j–n ascospores (arrows indicate mucilaginous polar appendages) o a germinating ascospore p, q colony on PDA (p-front and q-reverse views). Scale bars: 200 µm (d); 100 µm (i); 50 µm (e, g, h); 10 µm (f, j–o).
Nucleotide comparisons of Dothiorella macadamiae (SDBR-CMU512) with D. albiziae (MFLUCC 22-0057) and D. thailandica (MFLUCC 11-0438) based on ITS, tef1-α, and TUB2.
Taxa | ITS | tef1-α | TUB2 |
---|---|---|---|
D. albiziae (MFLUCC 22-0057) | 5/471 bp (1.0%, without gaps) | 25/244 bp (10.2%, 8 gaps) | 3/340 bp (0.8%, without gaps) |
D. thailandica (CBS 133991) | 8/519 bp (1.5%, without gaps) | 30/285 bp (10.5%, 7 gaps) | 6/371 bp (1.6%, without gaps) |
Phaeoacremonium (Togniniaceae Réblová, L. Mostert, W. Gams & Crous) was introduced by
“chiangmaiense” refers to the location “Chiang Mai,” from where the holotype was collected.
CMUB 40065.
Saprobic on dead twigs of Macadamia sp. Sexual morph: Ascomata 80–165 × 115–170 µm (x̄ = 136 × 146 µm, n = 20), immersed, solitary, globose to subglobose, dark brown to black, glabrous, ostiole with a long neck, neck straight or flexuous. Perithecial necks 75–160 µm high × 15–35 µm diam. (x̄ = 119 × 25 µm, n = 20), cylindrical, periphysate, ostiolar canals sulcate. Peridium 18–40 µm wide (x̄ = 29.4 µm, n = 25), comprising two section layers, the inner section layer composed of hyaline cells of textura prismatica, the outer section layer comprising brown to dark brown cells of textura prismatica. Hamathecium 3.5–5.5 µm wide (x̄ = 4.6 μm, n = 50), comprising cylindrical, hyaline, septate paraphyses, slightly inflated between the septa near their base and slightly contracted at the septa, longer than the asci. Asci 17–27 × 4–6 µm (x̄ = 20 × 5 µm, n = 50), 8-spored, arising in acropetal succession, unitunicate, apedicellate, cylindrical to clavate, apically rounded to truncate. Ascogenous hyphae hyaline, smooth-walled, septate, simple, 4–5 µm (x̄ = 4.7 μm, n = 10) at the base. Ascospores 4–7 × 1.2–2.5 µm (x̄ = 5.5 × 2 µm, n = 35), overlapping, hyaline, oblong to allantoid, aseptate, smooth-walled, rounded, and small guttules at both ends. Asexual morph: Undetermined.
Ascospores germinating on PDA within 24 h at 28 °C, colony on PDA reaching 3 cm diam. after two weeks, culture from above flat, smooth surface, entire edges, white-yellow, low convex at the middle, forming tufts on the surface, wrinkled, reverse white to light reddish-brown from the edge to the center, wrinkled.
Thailand • Chiang Mai Province, 18°52'43"N, 99°13'15"E, 384 m elevation, on a dead branch of Macadamia sp., 24 November 2023, Xian Zhang, TCMM19, CMUB 40065 holotype; ex-type living culture, SDBR-CMU510, other living culture SDBR-CMU511.
The phylogenetic analyses showed that our isolates of Phaeoacremonium chiangmaiense formed an independent lineage that is basal to three species of Phaeoacremonium (P. iranianum (CBS 101357, CBS 117114), P. minimum (CBS 246.91, CBS 100397), and P. tuscanum (CBS 123033)) with 87% ML and 1.00 PP support (Fig.
RAxML tree based on a combined dataset of ITS+TUB2+ACT+tef1-α+LSU gene sequences data, which comprised 2709 base pairs (ITS = 1–600 bp, TUB2 = 601–1237 bp, ACT = 1238–1499 bp, tef1-α = 1500–1816 bp, LSU = 1817–2709 bp). The best-scoring RAxML tree with a final ML optimization likelihood value of -29286.426085 is presented. The matrix had 1400 distinct alignment patterns, with 45.73% of undetermined characters or gaps. Estimated base frequencies were as follows: A = 0.224680, C = 0.290688, G = 0.252826, T = 0.231806; substitution rates: AC = 1.466922, AG = 3.673841, AT = 1.453901, CG = 1.188614, CT = 5.325402, GT = 1.000000; proportion of invariable sites I = 0.385193; gamma distribution shape parameter α = 0.691947. The bootstrap support values for ML are equal to or greater than 60%, and BI analysis values are equal to or greater than 0.90 PP at each node. The tree is rooted with Flabellascus tenuirostris (CBS 138680) and Jattaea algeriensis (STE-U 6201). Newly generated species are shown in red, while the ex-type species are shown in bold.
Phaeoacremonium chiangmaiense (CMUB 40065, holotype) a, b appearance of ascomata on the host surface c vertical section of ascomata d ostiolar canal e section of peridium f, g hamathecium h–j ascogenous hyphae with asci attached k–n ascospores o germinated ascospores p, q colonies on PDA (p-front and q-reverse views). Scale bars: 100 µm (c); 40 µm (d); 20 µm (e, f, h–j); 10 µm (g, o); 5 µm (k–n).
Nucleotide comparisons of P. chiangmaiense (SDBR-CMU510) with P. iranianum (CBS 101357), P. minimum (CBS 246.91), and P. tuscanum (CBS 123033), based on ITS, TUB2, ACT, and tef1-α; all of them were compared excluding gaps.
Taxa | ITS | TUB2 | ACT | tef1-α |
---|---|---|---|---|
P. iranianum (CBS 101357) | 20/544 bp (3.7%) | 60/528 bp (11.3%) | 37/260 bp (14.2%) | 43/264 bp (16.3%) |
P. minimum (CBS 246.91) | 22/522 bp (4.2%) | 78/498 bp (15.7%) | 42/260 bp (16.2%) | 42/267 bp (15.7%) |
P. tuscanum (CBS 123033) | 22/585 bp (3.8%) | 70/596 bp (11.7%) | 40/237 bp (16.9%) | 38/269 bp (14.1%) |
Melomastia (Pleurotremataceae Walt. Watson) was established by
ZHKU 23-0106.
Saprobic on dead twigs of Macadamia integrifolia. Sexual morph: Ascomata 300–650 × 430–605 µm (x̄ = 521 × 516 µm, n = 20), visible as black raised dots on the host surface, solitary, semi-immersed, dark brown to black, subglobose or irregular pyriform, carbonaceous, papillate. Ostiolar central, carbonaceous, brown to dark brown to black. Peridium 30–80 µm wide (x̄ = 52 µm, n = 25), comprising two section layers: inner section layers hyaline to brown cells of textura prismatica, outer section layer, brown to black cells of textura prismatica. Hamathecium 2–4 µm wide (x̄ = 3 μm, n = 60), hyaline, filiform, septate, branched, pseudoparaphyses, longer than asci. Asci 166–235 × 5.6–9.9 µm (x̄ = 196 × 8.1 µm, n = 30), 8-spored, bitunicate, cylindrical, flexuous, smooth-walled, apically obtuse, with an ocular chamber, short-pedicellate. Ascospores 19–30 × 5–8 µm (x̄ = 26 × 6.5 µm, n = 50), hyaline, fusiform, uniseriate, 3-septate, narrow towards the apex and obtuse or conical ends, constricted at the septa, smooth-walled, without mucilaginous sheath or appendages, with guttules in each cell. Asexual morph: Undetermined.
China • Yunnan Province, Baoshan City, 24°48'18"N, 99°22'36"E, 1199.6 m elevation, on a dead branch of Macadamia integrifolia, 30 July 2022, Xian Zhang, MBC85, GMB1173, GMB1174, new host record.
On a dead branch of Hevea brasiliensis in China (
Macadamia puerensis was reported by
RAxML tree based on a combined dataset of LSU+SSU+tef1-α gene sequences. The topology of the trees generated by both maximum likelihood (ML) and Bayesian inference (BI) analyses exhibited high similarity. The RAxML tree with a final ML optimization likelihood value of -12726.117171. The aligned matrix had 933 distinct alignment patterns, with 22.57% of undetermined characters or gaps. Parameters for the GTR+I+G model of the combined LSU, tef1-α, and the SYM +I+G model of the combined SSU were as follows: estimated base frequencies A = 0.239385, C = 0.263871, G = 0.289634, T = 0.207110; substitution rates AC = 0.851550, AG = 2.071762, AT = 1.124966, CG = 0.971448, CT = 7.978925, GT = 1.000000; the proportion of invariable sites I = 0.478880; and gamma distribution shape parameter α = 0.661129. Bootstrap support values for ML equal to or greater than 60% and PP equal to or greater than 0.90 are given above the nodes. New records are in blue, while the ex-type species are in bold.
Melomastia puerensis (GMB1173, new host record) a, b appearance of ascomata on the host surface c section of peridium d vertical section of an ascoma e ocular chamber in lactophenol cotton blue reagent f, k hamathecium g–j asci l–p ascospores. Scale bars: 200 µm (d); 100 µm (g–j); 50 µm (k); 20 µm (c, f); 10 µm (l–p).
Melomastia species | Host records | Location | References |
---|---|---|---|
M. antarctica | Pernettya mucronata | Argentina |
|
M. aquatica | Unknown | China |
|
M. aquilegiae | Aquilegia karelini | Kirghiz SSR | Index Fungorum (2024) |
M. calami | Calamus sp. | Philippine |
|
M. beihaiensis | Chromolaena odorata | China |
|
M. calligoni | Calligonum sp. | Central Asia |
|
M. carinata | Ephedra | Iran |
|
M. chilensis | Sophora macrocarpa | Chile |
|
M. chromolaenae | Chromolaena odorata | Thailand |
|
M. clematidis | Clematis sikkimensis | Thailand |
|
M. clypeata | Salix martiana | Brazil |
|
M. coffeae | Coffea robusta | Central African Republic |
|
M. constricta | Malus turkmenorum | Turkmen SSR |
|
Cydonia oblonga | Central Asia |
|
|
M. corylina | Corylus sp. | Luxemburg |
|
M. distoseptata | Undetermined dead branch | India |
|
M. fulvicomae | Clematis fulvicoma | Thailand |
|
M. fusispora | Olea europaea | China |
|
M. graminicola | Sorghum vulgare | French |
|
M. haloxyli | Haloxylon aphyllum | Kazakh SSR | Index Fungorum (2024) |
M. heteroderma | Unknown | Cuba |
|
M. heveae | Hevea brasiliensis | Africa |
|
M. hyalostoma | Cola vera | Ivory Coast |
|
M. italica | Vitis vinifera | Italy |
|
M. jaapiana | Betulaceae | Germany |
|
M. kazachstanica | Ammodendron conollyi Haloxylon persicum | Central Asia |
|
M. lignicola | Betula pendula | Germany |
|
M. loropetalicola | Loropetalum chinense | China |
|
M. mastoidea | Chaenomeles speciosa | Ukraine |
|
Cornus sanguinea | Denmark |
|
|
Deutzia corymbosa | India |
|
|
Fraxinus excelsior | Poland |
|
|
Fraxinus sp. | Denmark |
|
|
Lantana involucrata | Bermuda |
|
|
Lonicera periclymenum | Denmark |
|
|
Lonicera quinquelocularis | India |
|
|
Lonicera xylosteum | Denmark |
|
|
Poland |
|
||
Russia |
|
||
Osmanthus fragrans | Ukraine |
|
|
Populus tremula | Denmark |
|
|
Rubia peregrina | Portugal |
|
|
Sambucus nigra | Denmark |
|
|
Symphoricarpos sp. | Denmark |
|
|
Syringa sp. | Denmark |
|
|
England |
|
||
Viburnum opulus | Denmark |
|
|
Poland |
|
||
Metasphaeria macounii | Canada, British, Columbia |
|
|
M. metasequoiae | Metasequoia glyptostroboides | Ukraine |
|
M. mangrovei | Rhizophora sp. | Thailand |
|
M. maolanensis | Undetermined dead branch | China |
|
M. marinospora | Kandelia candel | Brunei |
|
M. neothailandica | Rhizophora sp. | Thailand |
|
M. nigrificans | Salicis | Luxemburg |
|
M. oleae | Olea europaea | China |
|
M. pallidispora | Trematosphaeria pallidispora | Italy |
|
M. phetchaburiensis | Rhizophora apiculata | Thailand |
|
M. popuschoji | Amygdalus turcomanica | Turkmen SSR |
|
M. prorumpens = Trematosphaeria prorumpens = Zignoëlla prorumpens | Pine | Germany |
|
M. puerensis | Hevea brasiliensis | China |
|
M. pyriformis | Undetermined dead branch | China |
|
M. rhizophorae | Rhizophora apiculata | Thailand |
|
M. salicicola = Zignoëlla salicicola | Salix alba | Vaucluse Galliae |
|
M. saxauli | Haloxylon persicum, Salsola arbuscula, Salsola rigida |
|
|
M. sedi | Sedum acre | Crimean SSR |
|
M. septata | Undetermined dead branch | Thailand |
|
M. sichuanensis | Olea europaea | China |
|
M. sinensis | Camellia sinensis | Thailand |
|
M. shastensis | Abies magnifica var. shastensis | California |
|
M. thailandica | Marina cvicennia | Thailand |
|
M. thamplaensis | Undetermined dead branch | Thailand |
|
M. tiomanensis | Rhizophora sp. | Malaysia |
|
M. winteri | Olea europaea | China |
|
M. yezoensis | Sasa kurilensis | Japan |
|
This study introduces two new species and one new record of microfungi isolated from macadamia based on morphological and phylogenetic analyses. Botryosphaeriaceae contains 22 genera (
Togniniaceae contains one genus (Phaeoacremonium), and Pleurotremataceae contains three genera (Dyfrolomyces, Melomastia, and Pleurotrema) (
Xian Zhang thanks Faculty of Science, Chiang Mai University for granting a tuition-fee scholarship (Active Recruitment 2023) for her master’s study. The authors extend their appreciation to the researchers supporting Project Number (RSP2025R56), King Saud University, Riyadh, Saudi Arabia. Nakarin Suwannarach thanks Chiang Mai University for the partial support. The authors would like to thank the Program of Doctoral Innovation Research Team from Qujing Normal University for support. We thank Dr. Shaun Pennycook for suggestions for Latin names for the new taxa.
The authors have declared that no competing interests exist.
No ethical statement was reported.
This research was supported by the Yunnan Revitalization Talents Support Plan (Young Talents Program and High-End Foreign Experts Program), Yunnan Provincial Department of Science and Technology “Zhihui Yunnan” plan (202403AM140023), the National Natural Science Foundation of China (No. 32260004 and 3246002), Yunnan Fundamental Research projects [Grant No. 202201AU070017], the Special Basic Cooperative Research Programs of Yunnan Provincial Undergraduate Universities (Grant No. 202101BA070001-209, 202101BA070001-279), Mee-mann Chang Academician Workstation in Yunnan Province (Grant No. 202205AF150002), the researchers supporting Project Number (RSP2025R56), King Saud University, Riyadh, Saudi Arabia, and the Key Laboratory of Yunnan Provincial Department of Education of the Deep-Time Evolution on Biodiversity from the Origin of the Pearl River.
Conceptualization: SCK, NS, and ST. Data curation: XZ, LSH. Formal analysis: SCK, JK. Funding acquisition: ST, SCK, DQD, NS, HHW. Investigation: ST, XZ, NS. Methodology: CS, TYD, AME, XZ, ST, JK. Project administration: NS, SCK. Resources: XZ. Software: TYD, AME, LSH, XZ. Supervision: HHW, SCK. Validation: NS. Visualization: ST, XZ. Writing - original draft: XZ. Writing - review and editing: DQD, CS, NS, JK, AME, LSH, TYD, HHW, XZ, SCK, ST.
Xian Zhang https://orcid.org/0000-0001-6097-8922
Saowaluck Tibpromma https://orcid.org/0000-0002-4706-6547
Samantha C. Karunarathna https://orcid.org/0000-0001-7080-0781
Tian-Ye Du https://orcid.org/0000-0003-2105-1803
Li-Su Han https://orcid.org/0000-0001-5380-9928
Abdallah M. Elgorban https://orcid.org/0000-0003-3664-7853
Jaturong Kumla https://orcid.org/0000-0002-3673-6541
Chanokned Senwanna https://orcid.org/0000-0002-1008-4514
Dong-Qin Dai https://orcid.org/0000-0001-8935-8807
Nakarin Suwannarach https://orcid.org/0000-0002-2653-1913
Hao-Han Wang https://orcid.org/0000-0002-2128-7894
All of the data that support the findings of this study are available in the main text.