Research Article |
Corresponding author: Ekachai Chukeatirote ( ekachai@mfu.ac.th ) Corresponding author: Saowaluck Tibpromma ( saowaluckfai@gmail.com ) Academic editor: Ajay Kumar Gautam
© 2025 Mei-Yan Han, Samantha C. Karunarathna, Li Lu, De-Ge Zheng, Nakarin Suwannarach, Abdallah M. Elgorban, Dong-Qin Dai, Li-Juan Zhang, Wan-Tong Zhao, Ekachai Chukeatirote, Saowaluck Tibpromma.
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:
Han M-Y, Karunarathna SC, Lu L, Zheng D-G, Suwannarach N, Elgorban AM, Dai D-Q, Zhang L-J, Zhao W-T, Chukeatirote E, Tibpromma S (2025) Three new species of Camporesiomyces (Tubeufiaceae, Tubeufiales, Dothideomycetes) associated with coffee in Yunnan, China. MycoKeys 117: 315-336. https://doi.org/10.3897/mycokeys.117.154573
|
During our surveys of microfungi associated with coffee plants in Yunnan Province, China, three saprobic fungi were isolated from dead coffee branches. Multigene phylogenetic analyses (ITS, LSU, tef1-α, and rpb2) and morphological characteristics resulted in the identification of three novel species in Camporesiomyces, namely C. bhatii, C. coffeae, and C. puerensis. Detailed morphological descriptions, illustrations, and phylogenetic analyses of these three new species are provided, along with morphological comparisons to closely related taxa. These findings have global implications for understanding the diversity of microfungi associated with coffee trees.
Coffee-associated microfungi, morphology, multigene phylogeny, taxonomy, three novel species
Coffee (roasted beans of Coffea) is one of the most widely consumed beverages globally and ranks as the second most valuable traded commodity after crude oil (
The history of coffee cultivation in China can be traced back to 1884, when British merchants first introduced coffee to Taiwan Province of China. The earliest coffee cultivation on mainland China began in the early 20th century when French missionaries brought coffee seedlings to Binchuan County, Yunnan Province, for cultivation (
Microfungi play a crucial role in agricultural ecosystems and are key nutrient regulating organisms in coffee ecosystems. They are globally distributed and exhibit high diversity (Bahram & Netherway, 2022). In recent years, scientists have conducted in-depth research on coffee-related fungi, identifying approximately 600 species (
Tubeufiales was established by
Tubeufiaceae was established by
Based on integrated morphological and phylogenetic analyses of combined ITS, LSU, tef1-α, and rpb2 sequences, Camporesiomyces was proposed by
In this study, we aim to introduce three new species of Camporesiomyces found on dead coffee branches collected from Pu’er and Xishuangbanna, Yunnan Province, China. Morphological characteristics and multi-locus phylogeny analyses, based on combined ITS, LSU, tef1-α, and rpb2 sequence data, were used to confirm the taxonomic placement of the new species.
Dead branches of Arabica Coffee and Liberica Coffee with fungal fruiting bodies were collected from Pu’er City and Xishuangbanna Autonomous Prefecture, Yunnan Province, China. The specimens were placed in self-sealing bags, with sampling information recorded on the cover (
Fresh mycelia grown on PDA (15–30 days) were scraped and transferred to 1.5 ml microcentrifuge tubes for DNA extraction. The Biospin Fungus Genomic DNA Extraction Kit BSC14S1 (BioFlux, China) was used to obtain DNA products from the above mycelia, according to the manufacturer’s protocol. The extracted DNA was then preserved at -20 °C for future use. Polymerase chain reaction (PCR) was performed for four loci; the details of the different loci, primers, and PCR thermal cycle programs used for amplification are presented in Table
Loci | Primers | PCR conditions | References |
---|---|---|---|
ITS | ITS5/ITS4 | 94 °C: 3 mins, (94 °C: 45 s, 55 °C: 50 s, 72 °C: 1 min) × 35 cycles, 72 °C: 10 mins |
|
LSU | LR5/LR0R |
|
|
tef1-α | 983F/2218R |
|
|
rpb2 | 5F/7cR | 95 °C: 5 mins, (95 °C: 1 min, 55 °C: 2 min, 72 °C: 90 s) × 40 cycles, 72 °C: 10 mins |
|
Names, strain numbers, and corresponding GenBank accession numbers of the taxa used in the phylogenetic analyses. The newly generated sequences are in bold black. After the strain number, “T” indicates the type strains. “NA” indicates sequence unavailability.
Fungal species | Strain numbers | ITS | LSU | tef1-α | rpb2 |
---|---|---|---|---|---|
Acanthohelicospora aurea | GZCC 16-0059 | KY321322 | KY321325 | KY792599 | MF589910 |
Acanthohelicospora aurea | GZCC 16-0060 | KY321323 | KY321326 | KY792600 | MF589911 |
Acanthostigma chiangmaiense | MFLUCC 10-0125T | JN865209 | JN865197 | KF301560 | NA |
Acanthostigma multiseptatum | ANM 475T | GQ856145 | GQ850492 | NA | NA |
Acanthostigma multiseptatum | ANM 665 | GQ856144 | GQ850493 | NA | NA |
Acanthostigma perpusillum | UAMH 7237T | AY916492 | AY856892 | NA | NA |
Acanthostigma scopulum | ANM 95 | GQ856142 | GQ850490 | NA | NA |
Acanthostigma scopulum | ANM 386 | GQ856141 | GQ850489 | NA | NA |
Berkleasmium fusiforme | MFLUCC 17-1978T | MH558693 | MH558820 | MH550884 | MH551007 |
Berkleasmium fusiforme | MFLUCC 17-1979 | MH558694 | MH558821 | MH550885 | MH551008 |
Bezerromyces pernambucoensis | URM7414 | KX470393 | KX518626 | KX518634 | NA |
Bezerromyces pernambucoensis | URM7412 | KX470391 | KX518624 | KX518632 | NA |
Camporesiomyces bhatii | GMBCC 1120T | PQ763360 | PQ842543 | PV388888 | PV388894 |
Camporesiomyces bhatii | GMBCC 1125 | PQ763361 | PQ842544 | PV388889 | PV388895 |
Camporesiomyces coffeae | GMBCC 1130T | PQ763358 | PQ842545 | PV388890 | PV388896 |
Camporesiomyces coffeae | GMBCC 1131 | PQ763359 | PQ842546 | PV388891 | PV388897 |
Camporesiomyces mali | KUMCC 19-0216T | NR_169709 | NG_075312 | MN794018 | NA |
Camporesiomyces patagoniensis | BBB MVB 573 | JN127358 | JN127359 | NA | NA |
Camporesiomyces puerensis | GMBCC 1113T | PQ763356 | PQ842541 | PV388886 | PV388892 |
Camporesiomyces puerensis | GMBCC 1114 | PQ763357 | PQ842542 | PV388887 | PV388893 |
Camporesiomyces vaccinii | CBS 216.90 | MH862204 | MH873889 | NA | NA |
Chlamydotubeufia cylindrica | MFLUCC 16-1130T | MH558702 | MH558830 | MH550893 | MH551018 |
Chlamydotubeufia huaikangplaensis | MFLUCC 16-0023 | KY678766 | KY678758 | KY792597 | MF535259 |
Dematiohelicosporum guttulatum | MFLUCC 17-2011T | MH558705 | MH558833 | MH550896 | MH551021 |
Dictyospora thailandica | MFLUCC 16-0001T | KY873627 | KY873622 | KY873286 | MH551023 |
Dictyospora thailandica | MFLUCC 18-0641 | MH558706 | MH558834 | MH550897 | MH551022 |
Helicoarctatus aquaticus | MFLUCC 17-1996T | MH558707 | MH558835 | MH550898 | MH551024 |
Helicoma muelleri | CBS 964.69T | AY916453 | MH871278 | NA | NA |
Helicoma nematosporum | MFLUCC 16-0011 | MH558722 | MH558848 | MH550913 | MH551039 |
Helicoma rubriappendiculatum | MFLUCC 18-0491T | MH558723 | NG_069593 | MH550914 | MH551040 |
Helicomyces hyalosporus | GZCC 16 0070 | MH558728 | MH558854 | MH550919 | MH551044 |
Helicomyces torauatus | MELUCC16-0217 | MH558732 | MH558858 | MH550923 | MH551048 |
Helicosporium nanningense | GZCC 22-2175T | OR066418 | OR066425 | OR058864 | OR058857 |
Helicosporium nanningense | GZCC 23-0588 | OR066419 | OR066426 | OR058865 | OR058858 |
Helicosporium setiferum | MFLUCC 17-1994T | MH558735 | MH558861 | MH550926 | MH551051 |
Helicosporium setiferum | MFLUCC 17-2007 | MH558737 | MH558863 | MH550928 | MH551053 |
Helicosporium sexuale | GZCC 22-2007 | OP508731 | OP508771 | OP698082 | OP698071 |
Helicosporium sexuale | MFLUCC 16-1244 | MZ538503 | MZ538537 | MZ567082 | MZ567111 |
Helicosporium vegetum | BCC 8125 | AY916491 | MH871277 | NA | NA |
Helicosporium vegetum | CBS 254.75 | NA | DQ470982 | DQ471105 | NA |
Helicosporium vegetum | CBS 941.72T | AY916488 | AY856883 | NA | NA |
Helicosporium vesicarium | MFLUCC 17-1795T | MH558739 | MH558864 | MH550930 | MH551055 |
Helicosporium viridiflavum | MFLUCC 17-2336T | MH558738 | NA | MH550929 | MH551054 |
Helicotubeufia guangxiensis | MFLUCC 17-0040T | MH290018 | MH290024 | MH290029 | MH290034 |
Helicotubeufia hydei | MFLUCC 17-1980T | MH290021 | MH290026 | MH290031 | MH290036 |
Kamalomyces thailandicus | MFLUCC 13-0233T | MF506884 | MF506882 | MF506886 | NA |
Kamalomyces thailandicus | MFLUCC 11-0158 | MF506883 | MF506881 | MF506885 | MF506887 |
Neoacanthostigma fusiforme | MFLUCC 11-0510T | KF301529 | KF301537 | NA | NA |
Neohelicomyces aquaticus | KUMCC 15-0463 | KY320529 | KY320546 | KY320562 | MH551065 |
Neotubeufia krabiensis | MFUCC 16-1125T | MG012031 | MG012024 | MG012010 | MG012017 |
Thaxteriella inthanonensis | MFLUCC 11-0003T | JN865211 | JN865199 | NA | NA |
Thaxteriellopsis lignicola | MFLUCC 16-0026 | MH558768 | MH558893 | MH550960 | MH551094 |
Thaxteriellopsis lignicola | MELUCC 16-0024 | MH558767 | MH558892 | MH550959 | MH551093 |
Tubeufia bambusicola | MFLUCC 17-1803T | MH558771 | MH558896 | MH550963 | MH551097 |
Tubeufia javanica | MFLUCC 12-0545T | KJ880034 | KJ880036 | KJ880037 | NA |
The assembly of the forward and reverse primers for the recently obtained sequence was accomplished using BioEdit version 7.0.5.3 and SeqMan version 7.0.0 software packages (DNASTAR, Madison, WI) (
Three new species (Camporesiomyces bhatii, C. coffeae, and C. puerensis) formed a distinct clade within Camporesiomyces, with strong statistical support (ML = 100%, PP = 1.00). Camporesiomyces clade was sister to Helicosporium within the family Tubeufiaceae, and all genera in Tubeufiaceae are distinctly separated (Fig.
The ML analysis and Bayesian inference (BI) analyses yielded nearly identical tree topologies, with bootstrap support values for ML equal to or greater than 50%, and BI analysis values equal to or greater than 0.90 PP are provided at each node. Newly generated sequences are in red, while the ex-type strains are in bold.
The best-scoring RAxML tree, with a final ML optimization likelihood value of -26167.982465 is presented. The matrix contained 1,266 distinct alignment patterns, with 21.89% of the characters being undetermined or gaps. Estimated base frequencies were as follows: A = 0.243800, C = 0.255403, G = 0.265360, T = 0.235437; substitution rates: AC = 0.936178, AG = 5.168669, AT = 2.168333, CG = 0.811291, CT = 10.165688, GT = 1.000000; gamma distribution shape parameter α = 0.816594. The phylogenetic tree resulting from RAxML analysis is shown in Fig.
According to the phylogenetic tree, our three new species (Camporesiomyces bhatii, C. coffeae, and C. puerensis) stand out as distinct entities within the genus Camporesiomyces. Among these, C. bhatii (GMBCC 1120 (ex-type) and GMBCC 1125) formed a sister lineage to C. vaccinii (CBS 216.90) with 96% ML/0.99 PP support. The other two new species, C. puerensis (GMBCC 1113 (ex-type) and GMBCC 1114) and C. coffeae (GMBCC 1130 (ex-type) and GMBCC 1131), clustered together but formed well-separated branches (98% ML/1.00 PP support), further emphasizing their distinctiveness. These two new species create a distinct branch that diverges from C. patagoniensis (BBB MVB 573) with 70% ML/0.99 PP support.
The species epithet “bhatii” honors Prof. Jayarama Darbhe Bhat for his immense contributions to mycology.
Camporesiomyces bhatii (GMB-W1176, holotype) a, b colonies on the natural substrate c, d conidiophore with conidia and conidiogenous cells e, f conidiogenous cells g–m conidia n a germinated conidium o, p Culture on PDA from above and reverse. Scale bars: 50 μm (c); 20 μm (d); 10 μm (e, g–n); 5 μm (f).
GMB-W1176
Saprobic on dead branches of Coffea liberica. Sexual morph: Undetermined. Asexual morph: Hyphomycetous. Colonies on natural substrate superficial, solitary or in clusters, hairy, yellow at apex, velvety. Mycelium exposed on the surface of substrate, except on the roots. Conidiophores 23–87.4 × 2.4–5.7 μm (x̄ = 63 × 3.7 μm, n = 30) macronematous, mononematous, either solitary or forming small clusters, smooth or occasionally verruculose, cylindrical, brown, slightly flexuous, simple, unbranched, 2–8-septate, sometimes slightly constricted at the septa. Conidiogenous cells 8.8–21.7 × 3–4.7 μm (x̄ = 17 × 3.5 μm, n = 20), polyblastic, terminal, cylindrical, slightly tapering, conspicuously denticulate on conidial secession, pale brown. Conidia 16–30 × 3.3–6.3 μm (x̄ = 24 × 6.5 μm, n = 30), solitary, acrogenous, cylindrical, obclavate or fusiform, 3–8-septate, subhyaline to pale brown.
Conidia germinate on PDA within 24 h at 28 °C, colony on PDA reaching 3 cm diam. after two weeks, circular or irregular, umbonate, with the entire margin dark brown to black. The reverse side displays predominantly black colonies with brown peripheral edges.
China • Yunnan Province, Xishuangbanna Autonomous Prefecture (21°41'N, 101°25'E 570 m), on dead branches of Coffea liberica, 27 August 2024, M.Y. Han & Tibpromma, (DL4 = GMB-W1176, holotype), ex-type GMBCC 1120, other living culture GMBCC 1125.
In the phylogenetic tree, our collection Camporesiomyces bhatii [GMBCC 1120 (ex-type) and GMBCC 1125] formed a well-separated lineage, sister to C. vaccinii [CBS 216.90, ex-type] with 96% ML/0.99 PP support (Fig.
The species epithet “coffeae” refers to the host genus Coffea.
GMB-W1181
Saprobic on dead branches of Coffea arabica. Sexual morph: Undetermined. Asexual morph: Hyphomycetous. Colonies on natural substrate superficial, solitary or clusters, hairy, subhyaline at apex, bunch of flowers-like, shiny. Mycelium exposed on surface of substrate, except on the roots. Conidiophores 43–97 × 2.8–4.5 μm (x̄ = 68 × 3.7 μm, n = 30) macronematous, mononematous, either solitary or forming small clusters, smooth or occasionally verruculose, cylindrical, dark brown, slightly flexuous, simple, unbranched, with longitudinal striations in the upper part, 2–9-septate, sometimes slightly constricted at the septa. Conidiogenous cells 8–21 × 2.3–4.6 μm (x̄ = 15.4 × 3.4 μm, n = 30), polyblastic, terminal, cylindrical, with longitudinal striations, with several conspicuous denticles at apex, brown. Conidia 20–50 × 3.3–6.5 μm (x̄ = 26.7 × 4.5 μm, n = 30), solitary, acrogenous, guttules, cylindrical, fusiform, 3–7-septate, subhyaline or hyaline, obtuse or conical at both ends.
Conidia germinating on PDA within 24 h at 28 °C, colony on PDA reaching 2 cm diam. after two weeks, circular or irregular, umbonate, with entire margin dark grey to black, irregular, the reverse side displays predominantly black colonies with brown peripheral edges.
China • Yunnan Province, Pu’er City, Simao District (22°35'53"N, 100°59'17"E, 1186.4 m), on dead branches of Coffea arabica, 8 August 2024, M.Y. Han & Tibpromma, (YYT15 = GMB-W1181, holotype), ex-type GMBCC 1130; other living culture GMBCC 1131.
Based on multi-gene phylogenetic analysis, Camporesiomyces coffeae [GMBCC 1130 (ex-type) and GMBCC 1131] is a distinct species, sister to C. puerensis [GMBCC 1113 (ex-type) and GMBCC 1114] with 98% ML/1.00 PP support (Fig.
Nucleotide comparisons of C. bhatii (GMBCC 1120, ex-type), C. coffeae (GMBCC 1130, ex-type), and C. puerensis (GMBCC 1113, ex-type) based on ITS, LSU, tef1-α, and rpb2; all of them were compared, excluding gaps.
Species | ITS (%) | LSU (%) | tef1-α (%) | rpb2 (%) |
---|---|---|---|---|
C. puerensis vs. C. coffeae | 44/628(7.0%) | 15/900(1.67%) | 61/914 (6.65%) | 118/1076 (10.97%) |
C. puerensis vs. C. bhatii | 42/484 (8.68%) | 17/895 (1.9%) | 70/924 (7.58%) | 165/1076 (15.76%) |
C. coffeae vs. C. bhatii | 37/515 (7.18%) | 12/901 (1.33%) | 76/969 (7.84%) | 172/1070 (16.07%) |
The name reflects the type location, “Pu’er” City, China.
Camporesiomyces puerensis (GMB-W1121, holotype) a, b colonies on the natural substrate c, d conidiophores and conidiogenous cells e, f conidiogenous cells g–l conidia m a germinated conidium n, o culture on PDA from above and reverse. Scale bars: 50 μm (c); 30 μm (d); 10 μm (e, f); 20 μm (g–m).
GMB-W1121
Saprobic on dead branches of Coffea arabica. Sexual morph: Undetermined. Asexual morph: Hyphomycetous. Colonies on the natural substrate are superficial, solitary or clusters, hairy, yellow at apex velvety. Mycelium exposed on the surface of the substrate, with the exception of the roots. Conidiophores 52–176.5 × 2.8–5.6 μm (x̄ = 88.4 × 4 μm, n = 30), macronematous, mononematous, solitary or forming small clusters, smooth or occasionally verruculose, cylindrical, brown, flexuous, simple, unbranched, 3–13-septate, sometimes slightly constricted at the septa. Conidiogenous cells 6.8–26 × 2.3–4.3 μm (x̄ = 15 × 3.6 μm, n = 20), polyblastic, terminal, cylindrical to slightly tapering, denticulate, smooth, slightly curved, pale brown. Conidia 21.7–83 × 4–9.4 μm (x̄ = 43.3 × 6.3 μm, n = 30), solitary, acrogenous, cylindrical, or fusiform, sometimes slightly curved, 4–9-septate, subhyaline to yellow, hyaline at both rostrate ends.
Conidia germinating on PDA within 24 h at 28 °C. Colony on PDA reaching 2 cm diam. after two weeks, circular or irregular, umbonate, with entire margin dark brown to black, reverse side displays predominantly black colonies with brown peripheral.
China • Yunnan Province, Pu’er City, Simao District (22°36'36"N, 101°0'14"E, 1189 m), on dead branches of Coffea arabica, 8 August 2024, M.Y. Han & Tibpromma, (BG25 = GMB-W1121, holotype), ex-type GMBCC 1113, other living culture GMBCC 1114.
In the concatenated phylogenetic analysis, our strain Camporesiomyces puerensis [GMBCC 1113 (ex-type), and GMBCC 1114] formed a sister branch with C. coffeae [GMBCC 1130 (ex-type), and GMBCC 1131] with 98% ML/1.00 PP bootstrap support (Fig.
In recent years, research on the diversity of microfungi on coffee plants has made significant progress, revealing their complexity and importance within the coffee ecosystem (
In this study, we introduce three new asexual species of Camporesiomyces, viz., C. bhatii, C. coffeae, and C. puerensis, which were isolated from decaying coffee branches based on morphological and phylogenetic analyses. Camporesiomyces currently comprises three recognized species: C. mali, C. patagonicus (sexual morph), and C. vaccinii (asexual morph) (
Species | Conidiophores | Conidiogenous cells | Conidia | Host/location | References |
---|---|---|---|---|---|
C. bhatii (GMB-W1176) | 23–87.4 × 2.4–5.7 μm, brown, 2–8-septate, sometimes slightly constricted at the septa. | 8.8–21.7 × 3–4.7 μm, cylindrical, slightly tapering, conspicuously denticulate on conidial secession, pale brown. | 16–30 × 3.3–6.3 μm, obclavate or fusiform, 3–8-septate, subhyaline to pale brown. | Dead branches of Coffea liberica/China | This study |
C. coffeae (GMB-W1181) | 43–97 × 2.8–4.5 μm, dark brown, with longitudinal striations in the upper part, 2–9-septate, sometimes slightly constricted at the septa. | 8–21 × 2.3–4.6 μm, cylindrical, with longitudinal striations, with several conspicuous denticles at apex, brown. | 20–50 × 3.3–6.5 μm, guttules, cylindrical, fusiform, 3–7-septate, subhyaline or hyaline, obtuse or conical at both ends. | Dead branches of Coffea arabica/China | This study |
C. puerensis (GMB-W1121) | 52–176.5 × 2.8–5.6 μm, brown, 3–13-septate, sometimes slightly constricted at the septa. | 6.8–26 × 2.3–4.3 μm, cylindrical to slightly tapering, denticulate, smooth, slightly curved, pale brown. | 21.7–83 × 4–9.4 μm, or fusiform, sometimes slightly curved, 4–9-septate, subhyaline to yellow, hyaline at both rostrate ends. | Dead branches of Coffea arabica/China | This study |
C. vaccinii (LMC 0043) | 64–145 μm long, 4.2–5.0 μm wide at base, 2.5–3.3 μm wide at apex, simple or rarely branched, dark brown with paleapical cells, occasionally roughened, 4–10 septate, usually with l–2 percurrent proliferations. | 1.3–2.7 × 1.0–1.3 μm, mono-or polyblastic, denticulate, denticles. | 8.0–13.0 μm, conidial filament coiled l.5–l.75 times, 4–8-septate, rounded at apex and tapering to a truncate base, 2.0–4.0 μm, hyaline to pale brown, smooth-walled. | On stems of Vaccinium elliotii/America |
|
We present descriptions, illustrations, and phylogenetic analysis results to validate and confirm the novelty of the three species: C. bhatii, C. coffeae, and C. puerensis. This is the first report of Camporesiomyces fungi associated with coffee, contributing three new asexual species to the genus Camporesiomyces. Our findings significantly enrich the fungal diversity of Camporesiomyces and hold substantial scientific value for enhancing our understanding of the fungal community inhabiting coffee hosts.
Samantha C. Karunarathan and Saowaluck Tibpromma thank the High-Level Talent Recruitment Plan of Yunnan Province (High-End Foreign Experts program) and the Key Laboratory of Yunnan Provincial Department of Education of the Deep-Time Evolution on Biodiversity from the Origin of the Pearl River for their support. Nakarin Suwannarach thanks Chiang Mai University for the support. The authors extend their appreciation to the Researchers Supporting Project Number (RSP2025R56), King Saud University, Riyadh, Saudi Arabia. The authors thank Dr. Shaun Pennycook for his advice on new fungal names.
The authors have declared that no competing interests exist.
No ethical statement was reported.
This research was supported by the Yunnan Provincial Department of Science and Technology “Zhihui Yunnan” plan (202403AM140023), the National Natural Science Foundation of China (Number 32260004) and the Researchers Supporting Project Number (RSP2025R56), King Saud University, Riyadh, Saudi Arabia.
Conceptualization: SCK, ST, and LL. Data curation: MYH, LL, WTZ, and LJZ. Formal analysis: SCK, ST, and NS. Funding acquisition: SCK, ST, and AME. Investigation: SCK and MYH. Methodology: MYH, ST, EC. Project administration: ST, SCK. Resources: MYH. Software: MYH, LL, DGZ, WTZ and LJZ. Supervision: ST and SCK. Validation: ST and SCK. Visualization: ST and MYH. Writing - original draft: MYH. Review, and editing: SCK, ST, EC, LL, DGZ, NS, AME, DQD, LJZ, and WTZ.
Mei-Yan Han https://orcid.org/0009-0004-3705-5408
Samantha C. Karunarathna https://orcid.org/0000-0001-7080-0781
Li Lu https://orcid.org/0000-0003-0977-6414
De-Ge Zheng https://orcid.org/0009-0001-9199-6784
Nakarin Suwannarach https://orcid.org/0000-0002-2653-1913
Abdallah M. Elgorban https://orcid.org/0000-0003-3664-7853
Dong-Qin Dai https://orcid.org/0000-0001-8935-8807
Li-Juan Zhang https://orcid.org/0000-0002-3234-6757
Wan-Tong Zhao https://orcid.org/0009-0006-7298-0150
Ekachai Chukeatirote https://orcid.org/0000-0002-9968-5841
Saowaluck Tibpromma https://orcid.org/0000-0002-4706-6547
All of the data that support the findings of this study are available in the main text.