Research Article |
Corresponding author: Ruttapon Srisonchai ( ruttasr@kku.ac.th ) Academic editor: Luiz Felipe Iniesta
© 2025 Theemaporn Benchapong, Henrik Enghoff, Natdanai Likhitrakarn, Ratmanee Chanabun, Anuwat Aoonkum, Ruttapon Srisonchai.
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:
Benchapong T, Enghoff H, Likhitrakarn N, Chanabun R, Aoonkum A, Srisonchai R (2025) Three new species of the millipede genus Enghoffosoma Golovatch, 1993 from Thailand (Diplopoda, Polydesmida, Paradoxosomatidae). Zoosystematics and Evolution 101(2): 509-532. https://doi.org/10.3897/zse.101.144277
|
An examination of newly collected material from northeastern Thailand revealed three new millipede species in the genus Enghoffosoma Golovatch, 1993: E. furca sp. nov., E. parvispina sp. nov., and E. rubellum sp. nov. We utilized morphological evidence together with genetic analysis of the COI gene to validate species identification. Each of the new species has unique characteristics, particularly in the structure of the gonopodal solenophore, that easily distinguish them from other species in the genus. The interspecific p-distances calculated from 658-bp barcoding sequences between the new species and their congeners show ranges of 18.56% to 22.78% for E. furca sp. nov., 15.09% to 21.25% for E. parvispina sp. nov., and 15.09% to 22.37% for E. rubellum sp. nov., supporting the morphological distinctions. Preliminary phylogenetic trees using Bayesian inference (BI) and maximum likelihood (ML) indicate that the genus Enghoffosoma is monophyletic, with two of the newly described species forming distinct clades. Our findings increase the total number of Enghoffosoma species to 16 (including six in Thailand) and underscore a greater diversity within the genus than previously known. Morphological illustrations of the new species, derived from SEM micrographs, and a distribution map of all known species are included.
Biodiversity, identification, phylogeny, South-East Asia, taxonomy
Currently, Thailand stands out as a significant hotspot for millipede diversity with a total of 282 species identified across various habitats. Notably, approximately 84% of these species are endemic, highlighting the unique biodiversity in the country’s ecosystems. Among the families of millipede in the class Diplopoda, the family Paradoxosomatidae Daday, 1889 has long been recognized as the most species-rich group not only globally (1,600 species) but also in Thailand, comprising 106 out of the total 267 described species (
Research on the systematics of several groups of millipedes has progressed rapidly in Thailand in recent years (
The millipede genus Enghoffosoma Golovatch, 1993 is probably one of the most taxonomically diverse genera of Paradoxosomatidae in Asia, especially in South-East Asia (
In the present study, we examined recently collected specimens of Enghoffosoma from the northeastern region of Thailand. These specimens appear to be geographically isolated from other species. Morphological evidence obtained from scanning electron micrographs and molecular data support the identification of three distinct morphotypes, which are described here as new species.
The specimens were collected by hand from the northeastern region of Thailand. Live animals were photographed in the field using a Canon EOS 90D digital camera equipped with a Canon EF-S 60 mm f/2.8 Macro USM lens. All specimens were placed in 10% ethanol for proper euthanization based on the procedure of the American Veterinary Medical Association (
The collecting sites (latitude, longitude, and elevation) were recorded by a Garmin GPSMAP 60CSx, and all exact coordinates were double-checked with Google Earth. The background image for the distribution map was taken from the website Elastic Terrain Map (http://elasticterrain.xyz/) and altered with Adobe Photoshop CC2024.
This work was conducted under the approval of the Khon Kaen University’s Animal Care and Use Committee (Protocol Reviews No. IACUC-KKU-70/67) and with permission to conduct research in protected areas (Department of National Parks, Wildlife and Plant Conservation; Protocols No. 23103 and 23104).
All specimens were studied under a Nikon SMZ445 stereo microscope and were measured with a TIGA 6x1/128 in a vernier caliper. For scanning electron micrographs, dissected body parts were mounted on aluminum stubs, air-dried, coated with gold at 30 mAmp for 240 seconds, and then studied under a JEOL, JSM-5410 LV microscope. After imaging, all specimens were removed from the stubs and preserved in ethanol in Eppendorf tubes for long-term preservation. The images were processed and edited with Adobe Photoshop CC2024 and line drawings were sketched in Adobe Illustrator 2018.
The species identification and morphological descriptions followed
cl claw of leg,
cx coxa of leg,
ep epiproct,
fe femur of leg,
hy hypoproct,
oz ozopore,
pfe prefemur of leg,
pk pleurosternal keel,
pof postfemur of leg,
pp paraproct,
sc sensory cones,
ta tarsus of leg,
ti tibia of leg,
vb ventral brushes.
a process a on solenophore,
ca cannula of gonopod,
cox coxa of gonopod,
fm femur,
ls lateral sulcus,
p process p on femur,
prf prefemur of gonopod,
sg seminal groove of gonopod,
sl solenomere of gonopod,
sp spatula-like process,
sph solenophore of gonopod,
sph-l small process located on base of solenophore, clearly seen in lateral view,
sph-m small process located on base of solenophore, clearly seen in mesal view.
NHMD Natural History Museum of Denmark, University of Copenhagen, Denmark
A total of 24 sequences were used for phylogenetic study, comprising 17 sequences from this study and 7 sequences available from GenBank (Suppl. material
Genomic DNA was extracted from body rings and leg tissues using the NucleoSpin Tissue Kit. The COI gene was amplified via polymerase chain reactions (PCR) using the universal primers LCO1490 and HCO2198 — LCO-1490 (5ʹ-GGT CAA CAA ATC ATA AAG ATA TTG G-3ʹ) and HCO-2198 (5ʹ-TAA ACT TCA GGG TGA CCA AAA AAT CA-3ʹ) (
The DNA sequences were checked for errors by putting them into the same reading frame, translating into proteins, and all sequences were aligned using MEGA 11 (
The phylogenetic trees were reconstructed using, Bayesian inference (BI) and maximum likelihood (ML), through the online CIPRES Science Gateway (
Uncorrected pairwise genetic distances (p-distances) were also calculated in MEGA 11.
The COI sequence alignments consisted of 658 base pairs with 243 variable sites, 415 conserved sites, and 187 parsimony-informative sites. The range of percentual distances between outgroups (A. festiva, D. planata, O. setosa, O. gracilis, and T. roseiparaterga) and Enghoffosoma species was 17.82–30.85%. The range of interspecific distances within the genus Enghoffosoma was 15.09–22.78%. The interspecific distances between other Enghoffosoma species (vs the three new species) were 18.56–22.78% for E. furca sp. nov., 15.09–21.25% for E. parvispina sp. nov. and 15.09–22.37% for E. rubellum sp. nov. The largest intraspecific distance was found in E. parvispina with 1.88%. Additional information for genetic distances is provided in Suppl. material
The Bayesian inference and maximum likelihood analyses based on the COI gene fragment resulted in the same topology for the ingroup (Fig.
Phylogenetic tree of COI gene fragment based on Bayesian Inference (BI) and maximum likelihood (ML). Nodal numbers refer to Bayesian posterior probabilities (PP) from BI and bootstrap support (BS) from ML. Different colors represent three new species (Enghoffosoma furca sp. nov., Enghoffosoma parvispina sp. nov. and Enghoffosoma rubellum sp. nov.), corresponding to the live photographs. Scale bar indicates substitutions per site.
Class Diplopoda de Blainville in Gervais, 1844
Order Polydesmida Pocock, 1887
Family Paradoxosomatidae Daday, 1889
Subfamily Paradoxosomatinae Daday, 1889
Tribe Paradoxosomatini Daday, 1889
Enghoffosoma Golovatch, 1993: 8.
Enghoffosoma
— Golovatch, 2011: 262;
Enghoffosoma can be differentiated from the other seven genera within the tribe Paradoxosomatini by the following combination of characteristics: paraterga conspicuous as short keels; pleurosternal keel present on rings 2–17; gonopod telopodite suberect, typically branched, with a very long and slender solenophore; solenomere also very long; femur often bearing a distal process; a distinct demarcation between the femur and the postfemoral portion; solenophore frequently featuring process(es) or lobe(s) near its base. Enghoffosoma is distinguished from Ciliciosoma Verhoeff, 1940, Strongylosoma Brandt, 1883, and Substrongylosoma Golovatch, 1984 by the presence of a conspicuous demarcation between the femur and the postfemoral part of the gonopod. The gonopod telopodite in Enghoffosoma is branched, which contrasts with Haplogonomorpha Mršić, 1996, whose species lack an apical branch. Enghoffosoma is separated from Lohmanderodesmus Schubart, 1934 and Stosatea Gray, 1843 by the typically long and slender solenophore in Enghoffosoma compared to the stout lobe in Lohmanderodesmus and the broad lamella in Stosatea. Enghoffosoma differs from Tetrarthrosoma Verhoeff, 1898 in having a distal process on the femur and a solenophore that branches into 2 or 3 processes.
Sundanina spinipleura Carl, 1941.
1) E. bispinum Likhitrakarn, Golovatch & Panha, 2014. 2) E. funda Likhitrakarn, Golovatch & Panha, 2014. 3) E. furca sp. nov. 4) E. parvispina sp. nov. 5) E. rubellum sp. nov. 6) E. zebra Likhitrakarn, Golovatch & Panha, 2014.
Enghoffosoma bispinum
Likhitrakarn, Golovatch & Panha, 2014 in
Enghoffosoma bispinum
— Nguyen & Golovatch, 2016: 151;
Thailand • 3 ♂ 2 ♀ (
E. bispinum is currently known from only two locations in Rayong Province, situated just 10 km apart. All newly collected specimens were observed walking on the ground, with some found on thin layers of leaf litter.
Enghoffosoma funda
Likhitrakarn, Golovatch & Panha, 2014 in
Enghoffosoma funda
— Nguyen & Golovatch, 2016: 151;
Thailand • 15 ♂ 8 ♀ (
This species is restricted to a small area within the dry deciduous forests of Sisaket Province. The type locality and its surrounding area are characterized by sandstone habitats. A recent visit to the type locality confirmed the presence of several topotypic specimens, observed on humid ground covered with a thin layer of leaf litter.
Holotype
: Thailand • ♂ (
Paratypes
: Thailand • 7 ♂♂, 5 ♀♀ (
Thailand • ♂ broken specimen (
The specific name is a Latin noun in apposition and refers to the two-pronged fork solenomere.
Distal region of femoral part of gonopod without process, and tip of solenomere forklike (Fig.
Holotype: length 30.4 mm, height of midbody 3.5 mm, width of midbody prozonum and metazonum 3.0 and 3.5 mm, respectively. Males: length 27.5–29.4 mm, height of midbody 3.3–3.5 mm, width of midbody prozona and metazona 3.1–3.2 and 3.6–3.7 mm, respectively. Females: length 28.8–32.9 mm, height of midbody 3.5 mm, width of midbody prozona and metazona 3.5–3.6 and 4.0–4.1 mm, respectively. Width of head < collum = ring 2 = 3= 4 < 5–13 (male); head < ring 2 = 3 < 4 < collum < 5–13 (female), thereafter gradually tapering towards telson.
Coloration
(Fig.
Head and clypeolabral region
(Figs
Scanning electron micrographs of Enghoffosoma furca sp. nov., male paratype (
Scanning electron micrographs of Enghoffosoma furca sp. nov., male paratype (
Antennae
(Figs
Collum
(Fig.
Tegument
(Figs
Metaterga
(Figs
Limbus
(Fig.
Paraterga
(Figs
Pleurosternal keels
(Figs
Telson
(Fig.
Sterna
(Fig.
Legs
(Fig.
Gonopods
(Fig.
Scanning electron micrographs of Enghoffosoma furca sp. nov., right gonopod, male paratype (
(Fig.
In general appearance, E. furca sp. nov. is morphologically similar to E. anchoriforme
Holotype
: Thailand • ♂ (
Paratypes
: Thailand • 15 ♂♂, 15 ♀♀ (
Thailand • 56 ♂♂, 4 ♀♀; same data as for holotype • ♂ (
The specific name is a noun in apposition referring to a small process on the lateral base of the solenophore.
Femoral part of gonopod with a long process (p) distally (Fig.
Holotype: length 28.5 mm, height of midbody 3.1 mm, width of midbody prozonum and metazonum 3.0 and 3.5 mm, respectively. Males: length 26.5–28.9 mm, height of midbody 2.7–3.2 mm, width of midbody prozona and metazona 2.7–3.1 and 3.3–3.6 mm, respectively. Females: length 28.1–30.6 mm, height of midbody 3.4–3.7 mm, width of midbody prozona and metazona 3.3–3.6 and 3.6–3.7 mm, respectively. Width of head < collum < ring 4 < 2 = 3 < 5–13 (male); head < ring 2 = 3 = 4 < collum < 5–15 (female), thereafter gradually tapering towards telson.
Coloration
(Fig.
Head and clypeolabral region
(Figs
Scanning electron micrographs of Enghoffosoma parvispina sp. nov., male paratype (
Scanning electron micrographs of Enghoffosoma parvispina sp. nov., male paratype (
Antennae
(Fig.
Collum
(Fig.
Tegument
(Figs
Metaterga
(Figs
Scanning electron micrographs of Enghoffosoma parvispina sp. nov., holotype, right gonopod male paratype (
Limbus
(Fig.
Paraterga
(Figs
Pleurosternal keels
(Figs
Telson
(Fig.
Sterna
(Fig.
Legs
(Fig.
Gonopods
(Fig.
(Fig.
Based on gonopod solely, this species is very similar to E. rubellum sp. nov. It can be distinguished by the lateral small spine at base of solenophore (sph-l) (Fig.
Holotype
: Thailand • ♂ (
Paratypes
: Thailand • ♂ (
Thailand • 2 ♂♂, 3 ♀♀ (
The specific name is an adjective and refers to the remarkable reddish brown body color of living specimens.
Femoral part of gonopod with a long process (p) distally (Fig.
Holotype: length 25.7 mm, height of midbody 3.2 mm, width of midbody prozonum and metazonum 2.9 and 3.7 mm, respectively. Males: length 25.7–27.1 mm, height of midbody 3.2–3.6 mm, width of midbody prozona and metazona 2.7–3.1 and 3.3–3.6 mm, respectively. Females: length 28–35.4 mm, height of midbody 3.6–4.3 mm, width of midbody prozona and metazona 3.6–4.1 and 4.0–4.7 mm, respectively. Width of head < collum < ring 3 = 4 < 2 < 5–15 (male), head < collum < ring 2 = 3 = 4 < 5–15 (female), thereafter gently tapering towards telson.
Coloration
(Fig.
Head and clypeolabral region
(Figs
Scanning electron micrographs of Enghoffosoma rubellum sp. nov., male paratype (
Scanning electron micrographs of Enghoffosoma rubellum sp. nov., male paratype (
Antennae
(Fig.
Collum
(Fig.
Tegument
(Figs
Metaterga
(Figs
Limbus
(Fig.
Paraterga
(Figs
Pleurosternal keels
(Figs
Telson
(Fig.
Sterna
(Fig.
Legs
(Fig.
Gonopods
(Fig.
Scanning electron micrographs of Enghoffosoma rubellum sp. nov., right gonopod, male paratype (
(Fig.
Among all recorded Enghoffosoma species, E. rubellum sp. nov. shows such vivid brownish red/red body color that it is easily distinguished from all congeners.
Enghoffosoma zebra
Likhitrakarn, Golovatch & Panha, 2014 in
Enghoffosoma zebra
—Nguyen & Golovatch, 2016: 151;
Thailand (Nakhon Ratchasima and Ubon Ratchathani provinces) and Laos (Champasak Province).
Outline of gonopod telopodite of some Enghoffosoma species that show similar conformation in the distal part, all in mesal view. A. E. bispinum
Outline of gonopod telopodite of some Enghoffosoma species that show similar conformation in the distal part, all in lateral view A. E. bispinum
This is the first treatment of the genus Enghoffosoma that integrates morphological characteristics and molecular data. It is important to note that the details of the gonopod telopodite clearly differentiate all species from each other. E. furca sp. nov. differs from the other species by the presence of three short processes, while E. parvispina sp. nov. and E. rubellum sp. nov. can be distinguished by the presence of small spinelike processes sph-l and sph-m, at the base of the solenophore, respectively. This aligns with the results of previous studies on paradoxosomatids and xystodesmids, where the integration of morphological and molecular approaches provided sufficient data for accurate species identification (
The interspecific genetic distances of the barcoding COI gene seem to agree with the morphological evidence: the three new species differ significantly from the others, showing high genetic distances of over 15.09% (18.56–22.78% for E. furca sp. nov., 15.09–21.25% for E. parvispina sp. nov. and 15.09–22.37% for E. rubellum sp. nov.). These values closely match the interspecific distances observed in other Paradoxosomatidae, which typically range from 11.4% to above 25% (
The number of described Enghoffosoma species found exclusively in Thailand has now reached six. Of these, five species inhabit sandstone environments, while only one is known to occur in limestone habitats. Field surveys conducted between 2020 and 2024 have revealed a greater diversity of Enghoffosoma in Thailand than previously thought. These findings indicate that northeastern Thailand is an important area for millipede biodiversity (Fig.
Distribution of all Enghoffosoma species. Color of symbol: black = type locality, white = other localities. Inset photo indicates locations of three new species. Type localities of E. digitatum Nguyen & Golovatch, 2016, E. retrosum Nguyen & Golovatch, 2016 and E. triangulare Nguyen & Golovatch, 2016 are in the same location and indicated by the same symbol.
The funding for this work was obtained from the National Science, Research and Innovation Fund through the Fundamental Fund of Khon Kaen University (2568), and partly supported by the Office of the Permanent Secretary, Ministry of Higher Education, Science, Research and Innovation (OPS MHESI), Thailand Science Research and Innovation (TSRI) (Grant No. RGNS65-056). The research capability enhancement program through graduate student scholarship (TB), Faculty of Science, Khon Kaen University, is gratefully acknowledged. Thanks are due to staff members of the Department of National Parks, Wildlife and Plant Conservation in Thailand who supported all available surveys in protected areas. We give special gratitude to ASRU and
List of species names used in this study and their COI accession numbers
Data type: docx
Uncorrected p-distances of cytochrome c oxidase subunit I (COI) sequence divergences in the genus Enghoffosoma
Data type: xlsx