Research Article |
Corresponding author: Chirasak Sutcharit ( jirasak4@yahoo.com ) Academic editor: Pavel Stoev
© 2024 Natdanai Likhitrakarn, Sergei I. Golovatch, Ruttapon Srisonchai, Parin Jirapatrasilp, Pichsinee Sapparojpattana, Ekgachai Jeratthitikul, Somsak Panha, Chirasak Sutcharit.
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:
Likhitrakarn N, Golovatch SI, Srisonchai R, Jirapatrasilp P, Sapparojpattana P, Jeratthitikul E, Panha S, Sutcharit C (2024) A new species of the pill millipede genus Rhopalomeris Verhoeff, 1906 (Diplopoda, Glomerida, Glomeridae) from Myanmar, and notes on Rhopalomeris carnifex (Pocock, 1889). ZooKeys 1215: 235-257. https://doi.org/10.3897/zookeys.1215.130919
|
The taxonomy of the pill millipede genus Rhopalomeris Verhoeff, 1906, which is restricted to Indochina and currently comprises six described species, is refined and updated. An integrative taxonomic approach was employed that combines morphological examination with DNA barcoding using the cytochrome c oxidase subunit I (COI) gene for species identification and delineation. The first objective was to confirm the identity of Rhopalomeris carnifex (Pocock, 1889), a charismatic species known as the “candy pill millipede” due to its vivid coloration, based on specimens collected near the type locality in Myanmar. The second objective was to describe a new species, Rhopalomeris nigroflava Likhitrakarn, sp. nov., discovered in Linno Gu, Kayin State, Myanmar. This new species is distinguished by its small body size (5.1–9.7 mm long) and yellow body with contrasting brown to blackish markings on certain terga. In addition, the position of the telopod syncoxital lobe relative to the lateral syncoxite horns separates it from other Rhopalomeris species. The interspecific divergence between R. nigroflava Likhitrakarn, sp. nov. and other congeners ranges from 10.85% to 16.13%, based on uncorrected COI p-distances, while the intraspecific divergence was 0%–7.44%. A distribution map of and a revised identification key to all known species of Rhopalomeris are also provided.
Biodiversity, candy pill millipede, key, systematics, taxonomy
The Oriental genus Rhopalomeris Verhoeff, 1906 consists of only six species, all of which are restricted to Indochina (
Distributions of all seven currently known Rhopalomeris species. Open triangles Rhopalomeris sauda Nguyen, Sierwald & Marek, 2019; filled circle Rhopalomeris nagao Nguyen, Nguyen & Eguchi, 2021; inverted filled triangle Rhopalomeris tonkinensis Silvestri, 1917; crossed circle Rhopalomeris nigroflava sp. nov.; filled squares Rhopalomeris carnifex (Pocock, 1889); red square Elphinstone Island; green square Kala Island; open circle Rhopalomeris variegata Golovatch, 2016; open diamond Rhopalomeris monacha Silvestri, 1917.
The genus Rhopalomeris belongs to the family Glomeridae. A total of 43 species of Glomeridae have so far been identified in Indochina and classified into six genera: Annameris Verhoeff, 1915 (two species), Hyleoglomeris Verhoeff, 1910 (23 species), Hyperglomeris Silvestri, 1917 (eight species), Rhopalomeris Verhoeff, 1906 (six species), Peplomeris Silvestri, 1917 and Tonkinomeris Nguyen, Sierwald & Marek, 2019 (one species each) (
Two unique morphological characters could be used to distinguish Rhopalomeris from the other glomerid genera: (1) antennomere 6 conspicuously enlarged, axe-shaped, exceeding the size of antennomeres 3–5 combined; (2) antennomere 7 also wide, topped by a disc-shaped antennomere 8 with numerous sensory cones, vs usually only four apical cones in other genera (except Peplomeris).
A well-known Rhopalomeris species is R. carnifex, commonly referred to as the “candy pill millipede” or “rainbow candy pill millipede” because of its striking and vibrant patterns (Fig.
The present study employs an integrative taxonomic approach, combining both traditional morphological examinations and DNA barcodes derived from a fragment of the COI gene. The aims of this study are to re-evaluate the taxonomy of R. carnifex by examining specimens collected from Koh Kala, Tanintharyi Division, Myanmar, and to describe a new species discovered at Linno Gu, Kayin State, Myanmar. We also provide a comprehensive distribution map and a revised identification key to all species currently recognized in this genus.
The new material was collected in Myanmar in 2015 and 2016 by SP and members of the Animal Systematics Research Unit, Chulalongkorn University, as well as by a French collecting team led by Louis Deharveng, of the
Muséum national d’Histoire naturelle (
The holotype and most paratypes are housed in the
Museum of Zoology, Chulalongkorn University (
The terminology used to describe the morphological structures is consistent with that applied in the most recent publications (
In the catalogue sections, D stands for the original description and subsequent descriptive notes; K for the appearance in a key; L for the appearance in a species list; M for a mere mention; MI for molecular information; and R for new subsequent records.
Total genomic DNA was extracted from the legs and part of the thoracic tissues using a DNA extraction kit for animal tissue (NucleoSpin Tissue Extraction Kit, Macherey-Nagel, Germany) following the standard procedure. The mitochondrial cytochrome c oxidase subunit I gene (COI: 660 bp) fragments were amplified using the primers LCO1490 and HCOoutout (
All nucleotide sequences obtained in this study have been deposited in the GenBank Nucleotide sequence database under the accession numbers PQ219547–PQ219550. The collecting localities and GenBank accession numbers for each nominal species are listed in Table
List of the species used for molecular phylogenetic analyses and their relevant information. * = paratype.
Species | Voucher number | Locality | GenBank accession number COI | Reference |
---|---|---|---|---|
Rhopalomeris carnifex | GLO016-1; GLO016-2 | Myanmar | PQ219547; PQ219548 | This study |
Rhopalomeris nigroflava sp. nov. | GLO093-1*; GLO093-2* | Myanmar | PQ219549; PQ219550 | This study |
Rhopalomeris sauda | IEBR-801; IEBR-706; IEBR-654; IEBR-533 | Vietnam | MT749398; MT749400; MT749401; MT749404 |
|
Rhopalomeris nagao | IEBR-852; IEBR-854 | Vietnam | MT749411; MT749392 |
|
Hyleoglomeris tongkerdae | MUMNH-GLO071-1* | Thailand | P493218 |
|
Hyleoglomeris bomba | MUMNH-GLO096-1* | Thailand | PP493219 |
|
Hyleoglomeris suwannakhuhensis | MUMNH-GLO039* | Thailand | PP493220 |
|
Hyleoglomeris nigromaculata | MUMNH-GLO019-1*; MUMNH-GLO019-2*; MUMNH-GLO019-3* | Thailand | PP493221; PP493222; PP493223 |
|
Hyleoglomeris dracosphaera | MUMNH-GLO001-1*; MUMNH-GLO001-2*; MUMNH-GLO001-3*; MUMNH-GLO035-M2*; MUMNH-GLO035-F4*; MUMNH-GLO035-M1*; MUMNH-GLO035-F3* | Thailand | PP493224; PP493225; PP493226; PP493227; PP493228; PP493229; PP493230 |
|
Hyleoglomeris krasoon | MUMNH-GLO059* | Thailand | PP493231 |
|
Hyleoglomeris hongkhraiensis | MUMNH-GLO029-2*; MUMNH-GLO029-3*; MUMNH-GLO031-3* | Thailand | PP493232; PP493233; PP493234 |
|
Hyleoglomeris awaumi | EG20210711-227-01; EG20210711-227-03; KS20210513-04; KS20210513-07 | Japan | LC713407; LC713409; LC713416; LC713419 |
|
Hyleoglomeris insularum | EG20201213-09 | Japan | LC713421 |
|
Hyleoglomeris japonica | MS20210617-01; MS20210617-02; MS20210617-03 | Japan | LC713422; LC713423; LC713424 |
|
Hyleoglomeris lucida | EG20210718-240-01; MS20210426-11 | Japan | LC713425; LC713426 |
|
Hyleoglomeris sulcata | MS20210521B-05 | Japan | LC713428 |
|
Hyleoglomeris uenoi | ST20211028 | Japan | LC713429 |
|
Hyleoglomeris halang | IEBR-Myr898P; IEBR-Myr926 | Vietnam | ON704753; ON704754 |
|
Hyleoglomeris lobus | SVE-204; IEBR-653; IEBR-678 | Vietnam | MT749391; MT749402; MT749406 |
|
Hyperglomeris bicaudata | CUMZ-GLO004*; CUMZ-GLO007* | Laos | OQ661871; OQ661872 |
|
Hyperglomeris inkhavilayi | CUMZ-GLO095/1*; CUMZ-GLO095/2* | Laos | OQ661873; OQ661874 |
|
Hyperglomeris simplex | IEBR-605; SVE-102 | Vietnam | MT749403; MT749410 |
|
Peplomeris magna | IEBR-677; IEBR-656 | Vietnam | MT749405; MT749408 |
|
Tonkinomeris napoensis | IEBR-804b; IEBR-804a | Vietnam | MT749396; MT749397 |
|
Trachysphaera costata | Tcost8-MK | Slovakia | KX467622 |
|
Glomeris marginata | ZFMK-TIS-18977; ZFMK-TIS-2517216 | France | MG892125; MG892167 |
|
Trachysphaera lobata | ZFMK:MYR TW01 | United; Kingdom | KJ408484 |
|
Trachysphaera schmidtii | ZFMK:MYR BGIMyr16 | Croatia | KJ408481 |
|
Eupeyerimhoffia archimedis | ZFMK:MYR1876 | Italy | KP205574 |
|
Sphaerobelum truncatum | CUMZ:2010.18 | Thailand | JN885184 |
|
Zephronia laotica | ZFMK:MYR3502 | Laos | MK330977 |
|
The sequences were aligned using MEGA7 (
The best-fit substitution model was determined using PartitionFinder2 v. 2.3.4 (
Intraspecific genetic distances within taxa that contained more than one individual and interspecific genetic distances based on the COI sequences were also calculated using uncorrected p-distances, as implemented in MEGA7 (
Family Glomeridae Leach, 1815
Rhopalomeris Verhoeff, 1906: 188 (D).
Rhopalomeris
–
The genus Rhopalomeris could be recognized through numerous apical sensory cones on the antennal tip, and antennomere 6 being particularly enlarged and rather strongly curved. In addition, the posterior telopods are rather strongly enlarged and stout, supplied with both prefemoral and femoral trichosteles. The femur has a distinctive and particular distocaudal process. The body is relatively large, with adults ranging from 11 to 20 mm in length. The body coloration is variable, but often useful for species identification.
Glomeris carnifex Pocock, 1889, fixed under Art. 70.3 (ICZN 1999) in
Rhopalomeris monacha Silvestri, 1917; R. tonkinensis Silvestri, 1917; R. variegata Golovatch & Semenyuk, 2016; R. sauda Nguyen, Sierwald & Marek, 2019; and R. nagao Nguyen, Nguyen & Eguchi, 2021.
The genus Rhopalomeris was originally typified by
Furthermore,
Subsequently,
Finally,
Peplomeris was originally described as a subgenus of Rhopalomeris (Silvestri, 1917). However, it was later raised to a genus level by
Glomeris carnifex Pocock, 1889: 290 (D). Type locality: Tenasserim.
Glomeris carnifex var. pallida Pocock, 1889: 290 (D). Type locality: Elphinstone Island.
Glomeris carnifex
–
Glomeris carnifex var. pallida
–
“Glomeris” bicolor [non Wood, 1865] –
Rhopalomeris carnifex
–
Rhopalomeris carnifex var. pallida
–
Myanmar, south Tenasserim (
Myanmar – Tanintharyi Division • 2 ♂; Myeik, Kala Island; elev. 5 m a.s.l.; 12°29'38″N, 98°30'53″E; 5 Apr. 2016; C. Sutcharit, W. Siriwut, R. Srisonchai leg.; CUMZ-GLO016-1, 16-2.
Body length of unrolled specimens, 17.5–17.9 mm (♂), width 8.9–9.1 (♂).
Color faded after 15 years of preservation in alcohol (Fig.
Labrum sparsely setose (Fig.
Rhopalomeris carnifex (Pocock, 1889), ♂ specimen from Kala Island A antenna, anterior view B leg 17, anterior view C leg 18 anterior view D, E telopod, anterior and posterior views F tip of syncoxital lobes (without scaling). Scale bars: 0.5 mm. Abbreviations: cx coxa, cxl coxal lobe, fe femur, fp femoral process, ft femoral trichostele, pf prefemur, pft prefemoral trichostele of telopod, sh syncoxital horn of telopod, sl syncoxital lobe of telopod, sn syncoxite notch, sx syncoxite, ta tarsus, ti tibia, tp tibial process.
Collum as usual, with two transverse striae. Thoracic shield with a small hyposchism field not projecting past rear tergal margin (Fig.
♂ legs 17 (Fig.
♂ legs 18 (Fig.
Telopods (= ♂ legs 19) (Fig.
The taxonomic status of R. carnifex presents a challenge.
In contrast, the description of G. carnifex var. pallida from Elphinstone Island contained far more detail.
Subsequently,
However, our recently obtained specimens from Kala Island, Myanmar (Fig.
Although there were previous reports of R. carnifex from several localities in southern Thailand, preliminary analyses of a number of R. carnifex specimens from this area reveal notable intraspecific variation in coloration, morphology, and molecular genetics, suggesting an occurrence of cryptic species (unpublished data). Therefore, a comprehensive redescription of newly retrieved male specimens from Kala Island, Myanmar in this study, comparing them with the original description of G. carnifex var. pallida from the nearby Elphinstone Island, is essential before any taxonomic revisions of other Thai specimens can be proposed. Furthermore, the morphological redescription of R. carnifex above is thus based only on these Myanmarese specimens.
Rhopalomeris (s.s.) monacha Silvestri, 1917: 143 (D).
Rhopalomeris monacha
–
This species was described from Perak State, western Malaysia (
Rhopalomeris (s.s.) tonkinensis Silvestri, 1917: 144 (D).
Rhopalomeris tonkinensis
–
This species was described from Tonkin, Montes Mauson, 2,000–3,000 ft. a.s.l., Lang Son Province, northern Vietnam (
Rhopalomeris variegata Golovatch & Semenyuk, 2016: 411 (D, K).
Rhopalomeris variegata
–
This species was described from Vietnam, Gia Lai Province, Kon Chu Rang Nature Reserve, 14°30'54″N, 108°32'47″E, ca 1,000 m a.s.l. (
Rhopalomeris sauda
Rhopalomeris sauda
–
This species was described from Vietnam, Bac Kan Province, Ba Be National Park, 400–500 m a.s.l.; Vinh Phuc Province, Phuc Yen Town, Ngoc Thanh Commune, Me Linh Station for Biodiversity, 21.385°N, 105.7119°E; Tam Dao district, Tam Dao National Park, 21.460945°N, 105.647021°E (
Rhopalomeris nagao
This species was described from Vietnam, Cao Bang Province, Pia Oac – Pia Den National Park, 22.5540°N, 105.8622°E, 850–1,600 m a.s.l. (
Holotype : Myanmar – Kayin State • ♂; Linno Gu (Lateral small cave); 16°50'52.9"N, 097°36'37.7"E; 25 Nov. 2015; F. Bréhier leg.; MY15-13/01-CUMZ-GLO093. Paratypes: Myanmar – Kayin State • 7 ♂♂ 5 ♀♀; same locality as holotype; MY15-13/01-CUMZ-GLO093) • 2 ♂♂ 2 ♀♀; same locality as holotype; MNHN-MY15-13/01.
Differs from other species of Rhopalomeris by the yellowish body with contrasting brown to blackish markings on terga 4–9 (Fig.
Body length of unrolled holotype,7.3 mm, width 4.1 mm. Body length of unrolled paratypes, 5.6–9.3 mm (♂), 5.1–9.7 mm (♀), width 3.1–4.8 (♂), 2.6–4.7 mm (♀).
Color faded after nine years of preservation in alcohol (Fig.
Labrum sparsely setose (Figs
Rhopalomeris nigroflava sp. nov., ♂ holotype A head and anterior part of body, ventral view B thoracic shield, lateral view C left antenna, frontal view D leg 18, anterior view. Scale bars: 1 mm (A, B), 0.2 mm (C), 0.4 mm (D). Abbreviations: cxl coxal lobe, e eggs, hy hyposchism field, K caudomedial tubercle, sc schism, sn syncoxite notch, sx syncoxite.
Collum as usual, with two transverse striae (Fig.
♂ legs 17 (Fig.
Rhopalomeris nigroflava sp. nov., A, B ♂ paratypes C–F ♂ holotype A–C leg 17, anterior view D, E telopod, posterior and anterior views, respectively F, G tip of syncoxital lobes (without scaling). Scale bars: 0.2 mm. Abbreviations: cx coxa, cxl coxal lobe, fe femur, fp femoral process, ft femoral trichostele, pf prefemur, pft prefemoral trichostele of telopod, sh syncoxital horn of telopod, sl syncoxital lobe of telopod, ta tarsus, ti tibia, tp tibial process.
♂ legs 18 (Fig.
Telopods (= ♂ legs 19) (Fig.
It seems noteworthy that a female and two male (Fig.
The specific epithet nigroflava is derived from the Latin niger meaning black and flavus meaning yellow, in reference to the dark bands on a yellowish dorsum, adjective in feminine gender.
1 | Thoracic shield yellowish or yellowish brown, contrasting to dark body background | 2 |
– | Thoracic shield dark, variegated | 3 |
2 | Body larger, ~ 12 mm in length and 6 mm in width. Head, collum and thoracic shield light yellowish, body mostly blackish with a light, rather broad, axial stripe. Perak State, western Malaysia | R. monacha |
– | Body the smallest, 5.1–9.7 mm in length and 2.6–4.7 mm in width. Body yellowish to yellowish brown-, contrasting to brown to blackish terga 4–9 (Fig. |
R. nigroflava sp. nov. |
3 | Body large-sized, ~ 20 mm in length and 11 mm in width (♀). Dorsum: mostly blackish, not variegated. Mount Mau Son, Lang Son Province, northern Vietnam | R. tonkinensis |
– | Body < 20 mm in length and 11 mm in width. Dorsum: with a varied color pattern or contrasting colors | 4 |
4 | Dorsum mostly dark, lateral edges of terga contrasting reddish or carmine. Telopod syncoxital lobe clearly higher than lateral syncoxite horns (Fig. |
R. carnifex |
– | Dorsum dark or light, sometimes variegated, lateral edges of terga neither reddish nor carmine. Telopod syncoxital lobe clearly shorter than lateral syncoxite horns, each latter without filament on top. Vietnam | 5 |
5 | Body almost entirely dark with contrasting four yellow lateral spots on each of terga 3–11. Telopods with a medially slightly concave syncoxite lobe. Prefemoral trichostele short, not extending to about half the distal boundary of femur. Femoral process (fp) long, narrow, erect, acute at tip. Tibial process short and lobuliform | R. nagao |
– | Body color variegated. Telopods with a medially slightly convex syncoxite lobe. Prefemoral trichostele long, extending to about half or distal boundary of femur. Femoral process large, subtrapeziform, rounded at tip. Tibial process long and sigmoid mesally | 6 |
6 | Body light brown to blackish, with variegated, marbled, brown-yellow to yellowish markings. Body larger, 15–18 mm in length and 7–8 mm in width. Antennomere 6 slightly shorter, ~ 2.0× as long as high. Thoracic shield with 8–10 striae. Central Vietnam | R. variegata |
– | Body blackish, with contrasting yellowish, lateral, oval bands on each of terga 2–11. Body smaller, 11 mm in length and 6 mm in width. Antennomere 6 longer, 2.5–3× as long as high. Thoracic shield with 2–3 striae. Northern Vietnam | R. sauda |
The COI alignment (Table
Bayesian inference tree (BI) of pill millipedes in the family Glomeridae based on 660 bp of COI gene. Clades of Rhopalomeris species in this study are highlighted in green. Numbers above branches indicate bipartition posterior probability (bpp) from Bayesian inference analysis (BI) and numbers below branches are Bootstrap Support (BS) values from the ML analysis.
The COI tree revealed a sister relationship between R. carnifex and R. nigroflava sp. nov., forming a well-supported clade 98% BS for ML and 0.99 bpp for BI. However, all Rhopalomeris species, including R. sauda and R. nagao, were not retrieved together as monophyletic (Fig.
The interspecific divergences based on COI uncorrected p-distance among the glomerid species in this study ranged from 9.74 to 19.87%, with an average of 14.87% (data not show). The interspecific divergences among Rhopalomeris species ranged from 10.85 to 16.13%, with an average of 13.32% (Table
Matrix of the average interspecific genetic divergence (uncorrected p-distance) for the 660 bp barcoding region of the COI gene between Rhopalomeris species.
Taxa | Rhopalomeris sauda | Rhopalomeris nagao | Rhopalomeris carnifex | Rhopalomeris nigroflava sp. nov. |
---|---|---|---|---|
Rhopalomeris sauda | 0.0744 | |||
Rhopalomeris nagao | 0.1289 | 0.0521 | ||
Rhopalomeris carnifex | 0.1613 | 0.1397 | 0.0175 | |
Rhopalomeris nigroflava sp. nov. | 0.1376 | 0.1235 | 0.1085 | 0 |
Currently, the genus Rhopalomeris comprises seven species distributed across Vietnam (four species), Myanmar (two species), Thailand and Malaysia (one species each), with a notable absence of documented sympatry of species. The distribution patterns (Fig.
Rhopalomeris species occur at elevations ranging between 5 and 1,600 meters above sea-level. Apparently, most are narrow endemics restricted to their type localities and are only rare to be encountered. The exceptions are R. sauda which boasts a wider distribution of roughly 180 kilometers, and R. carnifex that demonstrates a remarkably extensive range exceeding 1,200 kilometers and stretching from southern Myanmar through Thailand to northern Malaysia (Fig.
It seems noteworthy that our preliminary surveys of millipede diversity in southern Thailand have yielded a high level of variation in the shape and coloration of R. carnifex, and high genetic diversity, suggesting a high-level intraspecific variation or the occurrence of cryptic species (unpublished data). This observation underscores the need for further research to comprehensively understand the extent of this variation, ultimately paving the way for future studies to definitively identify the R. carnifex complex.
The interspecific divergence based on COI uncorrected p-distance among the Rhopalomeris species in this study ranged between 10.85–16.13%, aligning with previous findings for European Glomeris species (11.5–17.1%;
This consistency suggests that COI proves effective in identifying species-level differentiation within Glomeridae. Even such species as R. carnifex and R. nigroflava sp. nov. that are very closely related and form a well-supported clade show a significant p-distance of 10.85%. This indicates that there can be considerable variability in the COI gene even among closely related glomerid species.
This study investigates the intraspecific COI divergence within Glomeridae millipedes. The low genetic intraspecific differences observed in Rhopalomeris carnifex (1.75%) and the newly described R. nigroflava sp. nov. (0%) are consistent with previous reports on Peplomeris magna (0.2%;
Analyzing the COI gene sequence is highly valuable in determining species boundaries and enabling precise classifications of glomerid species. Unsurprisingly, most recent taxonomic studies on millipedes frequently employ this technique to distinguish between taxa. Unfortunately, the phylogenetic relationships in this study appear insufficient to resolve genus-level relationships within the family, as shown in this study and others (
We thank the members of the Animal Systematics Research Unit for their invaluable assistance in the field. We owe an obligation of gratitude to Henrik Enghoff, Duc-Anh Nguyen, and Thomas Wesener, the reviewers, whose constructive comments, grammatic corrections, and assistance have significantly improved the paper.
The authors have declared that no competing interests exist.
No ethical statement was reported.
This study was funded by the Thailand Science Research and Innovation Fund of Chulalongkorn University (grant number DIS66230007) and the TRF (Thailand Research Fund) Senior Research Scholar.
NL, RS, PJ, PS, EJ, SP, and CS collected and prepared specimens in the field. SP and CS provided financial and intellectual support. NL, RS, and EJ prepared specimens and wrote the manuscript. NL, RS, PJ, EJ, and CS conceived, designed, supervised the study, prepared figures, and approved and edited the final manuscript. SIG and SP reviewed, advised, and approved the final manuscript.
Natdanai Likhitrakarn https://orcid.org/0000-0002-1306-317X
Sergei I. Golovatch https://orcid.org/0000-0001-7159-5484
Ruttapon Srisonchai https://orcid.org/0000-0002-7142-0999
Parin Jirapatrasilp https://orcid.org/0000-0002-5591-6724
Pichsinee Sapparojpattana https://orcid.org/0009-0005-9134-4910
Ekgachai Jeratthitikul https://orcid.org/0000-0002-3477-9548
Somsak Panha https://orcid.org/0000-0002-4431-2458
Chirasak Sutcharit https://orcid.org/0000-0001-7670-9540
All of the data that support the findings of this study are available in the main text.