Research Article |
|
Corresponding author: Nikolay A. Poyarkov ( n.poyarkov@gmail.com ) Corresponding author: Pavel P. Skutschas ( p.skutschas@spbu.ru ) Academic editor: Deepak Veerappan
© 2025 Nikolay A. Poyarkov, Dana D. Skorinova, Andrey M. Bragin, Veniamin V. Kolchanov, Vladislav A. Gorin, Alexey V. Trofimets, Alexander P. Yuzefovich, Dac Xuan Le, Tan Van Nguyen, Pavel P. Skutschas.
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:
Poyarkov NA, Skorinova DD, Bragin AM, Kolchanov VV, Gorin VA, Trofimets AV, Yuzefovich AP, Le DX, Nguyen TV, Skutschas PP (2025) Integrative taxonomy reveals a new unstriped Ichthyophis Fitzinger, 1826 from Vietnam and provides new data on diagnostic osteological traits for Asian tailed caecilians (Gymnophiona: Ichthyophiidae). Vertebrate Zoology 75: 405-440. https://doi.org/10.3897/vz.75.e149399
|
Abstract
Herein we examined the cranial osteology of 15 species of Ichthyophis (I. asplenius, I. beddomei, I. glutinosus, I. kohtaoensis, I. larutensis, I. mindanaoensis, I. multicolor, I. nguyenorum, I. nigroflavus, I. sikkimensis, I. singaporensis, I. supachaii, I. tricolor, I. weberi, and Ichthyophis sp. from northern Vietnam) with a special emphasis on the temporal region. We presented the first detailed description of the cranium and the atlas of an Ichthyophis species based on micro-CT scanning data. We discuss the implications of temporal region composition for the systematics of this group and the evolution of the cranium in Gymnophiona as a whole. We further provided comments on a jaw-closing mechanism and reported on the presence of phylogenetically basal cranial features in ichthyophiids that are also found in stem caecilians. Our detailed morphological description was based on a specimen from a previously unknown population of unstriped Ichthyophis from northern Vietnam. We consequently described this population as a new species based on morphological and molecular (3967 bp from cyt b, 12S rRNA, and 16S rRNA mitochondrial DNA genes) lines of evidence. We provide comparisons of external morphological traits of the new species with its congeners and further compare its cranial osteological features with other Ichthyophis for which skull descriptions exist. The new species differs from the morphologically similar species I. yangi and I. chaloensis by a significant divergence in cyt b and 16S rRNA mitochondrial DNA gene sequences (p = 6.5%–6.9% and p = 4.5%, respectively). The new species is currently known only from evergreen forests of Xuan Lien National Park (Thanh Hoa Province) and Pu Hoat (Nghe An Province) Nature Reserve, northern Vietnam, and was recorded at elevations of 700–800 m asl. We suggest the new species be considered Data Deficient (DD), following the IUCN’s Red List categories.
Ichthyophis griseivermis sp. nov., micro-CT scanning, molecular analyses, morphology, Nghe An, osteology, Pu Hoat, taxonomy, Thanh Hoa, Xuan Lien
Caecilians (Gymnophiona) have a worm-like body plan lacking limbs and limb girdles; most species (at least as adults) are adapted to burrowing underground, though some members of this group have aquatic or semi-aquatic lifestyles or occasionally can be observed actively moving on the surface (e.g.,
The Asian tailed caecilians of the family Ichthyophiidae Taylor, 1968, comprise terrestrial limbless amphibians with aquatic larvae that are widely distributed across South and Southeast Asia from Sri Lanka and India through southern China, mainland Indochina, and Southeast Asian islands west of Wallace’s Line (
Ichthyophiids represent an ancient lineage of extant caecilians, forming a sister group to all other families except the Rhinatrematidae Nussbaum, 1977 (
The most speciose genus of Ichthyophiidae is Ichthyophis, which currently comprises 50 nominal species, 12 of which occur in Indochina (including Vietnam, Cambodia, Laos, and Thailand) and southern China (
Presently, four named Ichthyophis species are recorded from Vietnam: I. catlocensis, I. chaloensis, I. kohtaoensis Taylor, 1960 (including I. bannanicus as its junior synonym following
Fieldwork in Pu Hoat Nature Reserve in Nghe An Province, Vietnam, was conducted by N. A. Poyarkov and D. X. Le in May 2019, and in Xuan Lien National Park in Thanh Hoa Province, Vietnam, by A. M. Bragin, A. P. Yuzefovich, and D. X. Le in October 2023 (see Fig.
Distribution of the unstriped species of the genus Ichthyophis from Indochina and southern China. Symbols indicating Ichthyophis species are identical to those used in Figure
To study the skeletal morphology of Ichthyophis sp., we examined X-ray projections for two collected specimens (
Terminology and identification of the osteological features. We followed
Measurements and counts generally followed
In addition to measurements, counts were made of the following meristic characters: total number of annuli in dorsal count, from posterior margin of second collar (third nuchal groove) to tail cap (except the latter) (TAD); total number of annuli (except the tail cap) in ventral count (TAV); total number of annuli interrupted (at least one annular groove) by the cloacal disc (AV); total number of annuli posterior to the cloacal disc (TAT); the number of transverse grooves on dorsal surface of collar (TG); the number of premaxillary / maxillary teeth (PMM); the number of vomero-palatine teeth (VP); the number of outer mandibular teeth (DE, also sometimes referred to as dentary teeth); and the number of inner mandibular teeth (IM, also sometimes referred to as splenial teeth). All tooth counts included ankylosed teeth and empty sockets (counted from micro-CT scans). Comparative data for other species of Ichthyophis were obtained from the literature (e.g.,
Other abbreviations. Dist. – District; Is. – Island; NP: – National Park; NR – Nature Reserve; Prov.: – Province; WS – Wildlife Sanctuary.
For the molecular phylogenetic analyses, total genomic DNA was isolated using the standard phenol chloroform-proteinase K extraction procedures with consequent isopropanol precipitation for a final concentration of about 1 mg/mL (protocols followed
For both cyt b and 12S–16S rRNA fragments we used the modified PCR protocol of
To estimate the phylogenetic relationships of the genus Ichthyophis, we used the newly obtained cyt b, 12S rRNA, and 16S rRNA sequences of two unstriped specimens from Xuan Lien NP and Pu Hoat NR, together with previously published sequences of these genes for the family Ichthyophiidae, including 25 nominal species and six unnamed species of the genus Ichthyophis (of them 23 striped and eight unstriped species) and four species of the genus Uraeotyphlus Peters, 1880. Altogether, the concatenated alignment of cyt b, 12S rRNA, and 16S rRNA fragments included sequences from 46 representatives of Ichthyophiidae; the sequences of Epicrionops marmoratus Taylor, 1968 and Rhinatrema bivittatum (Guérin-Méneville, 1838) (family Rhinatrematidae) were used to root the tree (data summarized in Table S2).
We aligned the nucleotide sequences using the default parameters in MAFFT online (
We estimated phylogenetic trees for the concatenated mtDNA fragments data set. We inferred the phylogenetic relationships of Ichthyophis using Bayesian inference (BI) and maximum likelihood (ML) approaches. We used the IQ-TREE online server (
The following description is based on micro-computed tomographic reconstruction of the skull and the atlas of the female specimen
General skull morphology of the holotype of Ichthyophis griseivermis sp. nov. (
Braincase morphology of the holotype of Ichthyophis griseivermis sp. nov. (
3D reconstruction of the atlas of the holotype of Ichthyophis griseivermis sp. nov. (
General features of the skull. The cranium is heavily ossified, dorsoventrally flattened, and more like V-shaped in the dorsal view (with the maximum width at the level where the squamosal overlies the contact between the pterygoid and the quadrate) (Fig.
The braincase consists of the two well ossified compound bones: the sphenethmoid anteriorly and the os basale posteriorly (Fig.
The cranium bears the following ten major external openings: the foramen magnum and the paired external nostrils, orbital apertures, upper and lower temporal fossae, and choanae. Additionally, the palatal aspect of the cranium contains the openings of the large mediopalatinal cavities (= interpterygoid vacuities) and the lower temporal fossae (= ventral openings of the adductor chambers) (both openings are much larger than the orbital apertures). The nostrils are relatively large, about the same size as the orbital apertures, rounded with terminal and anterior orientation bordered by the premaxillae (medially and ventrally), nasals (dorsally), and septomaxillae (laterally) (Fig.
The dorsally oriented upper temporal fossae are elongated and slit-like, widening posteriorly (Fig.
The large foramen magnum (notably larger than the orbital aperture) faces posteriorly, is rounded and is entirely bordered by the os basale (Fig.
As in all caecilians (
Dentition. The dentigerous elements of the cranium (upper jaw) are the premaxillae, the maxillopalatine, and the vomers. The only dentigerous element of the lower jaw is the pseudodentary. The dentition is organized as two (labial and lingual) continuous sub-parallel rows of teeth on both the cranium and the lower jaw. The upper jaw dentition includes a labial ‘premaxillary-maxillary’ tooth row, which is present on the premaxillae and the maxillary portion of the maxillopalatine, and a lingual ‘vomero-palatine’ row (sensu
The lower jaw dentition includes a labial ‘dentary’ (or outer mandibular) tooth row that is located on the dentary portion of the pseudodentary and a lingual ‘splenial’ (or inner mandibular) tooth row (Fig.
All functional teeth are bicuspid and pedicellate. The crowns of functional teeth are moderately sized, with gently posteriorly curved tips. The teeth of labial rows (premaxillary-maxillary row on the cranium and outer mandibular tooth row on the lower jaw) are approximately the same size as the teeth of lingual rows (vomero-palatine row on the cranium and the inner mandibular row on the lower jaw). The replacement teeth are present lingually to the base of the corresponding functional teeth.
Premaxillae. The paired dentigerous premaxillae articulate with each other medially to form the anterior margin of the snout (Fig. S1). Each premaxilla contacts with the corresponding nasal dorsally, septomaxilla posterodorsally, vomer posteromedially, and maxillopalatine posteriorly, and consists of three main parts: the pars dorsalis (= dorsal or alary process), the pars dentalis, and the pars palatina. Pars dorsalis is relatively long and narrow with a dorsal apex that is sharply pointed and set within a deep notch of the nasal (Fig. S1A). There are two large neurovascular canals presumably for the ophthalmic division of the trigeminal nerve and the associated blood vessels, which nurture the skin, passing through the base and middle parts of the dorsal process. In the posterior view, these canals are visible as large foramina and in the medial view as short deep grooves (Fig. S1B). These canals branch within the bone and emerge on the external surface through several openings. The pars dentalis of each premaxilla is short and slightly curved along the contour of the snout (Fig. S1A). Its lingual surface bears closely spaced pedicellate teeth that form the anterior part of the labial tooth row.
The posterior end of each pars dentalis is penetrated by one oblique neurovascular canal that enters the bone posteriorly and almost immediately exits onto the dorsal surface, opening into a deep short groove. This canal is presumably for the passage of the medial branch of the maxillary division of the trigeminal nerve and the associated blood vessels. The pars palatina is a wide ledge along the lingual aspect of the pars dentalis (Fig. S1C, E). In ventral view, the inner margin of this bony ledge runs parallel to the curvature of the outer edge of the premaxilla, except the median part, where it has a pronounced projection (Fig. S1C). The pars palatina articulates along its entire length with the vomer to contribute to the anterior part of the palate.
Maxillopalatines. The paired maxillopalatines are compound dentigerous bones with a complex shape (Fig. S2). Each bone consists of a lateral maxillary part and a medial palatine part. Each maxillopalatine articulates with the premaxilla anteriorly, the septomaxilla anterodorsally, the prefrontal dorsally, the squamosal posteriorly, and the vomer and the pterygoid medially to contribute to the lateral part of the palate.
The maxillary part of each maxillopalatine has a short anterior process with a blunt anterior end (for the contact with the premaxilla anteriorly and the septomaxilla dorsally), the pars facialis (= facial, dorsal process) in the anterior half of the maxilla (for the contact with the prefrontal) (Fig. S2A), and the posterior process for the articulation with the squamosal posteriorly (Fig. S2D). The central part of each pars facialis is penetrated by a large foramen (Fig. S2A, B). The medial surface of the pars facialis around this large foramen bears a deep oval depression that is presumably holding the vomeronasal organ (Fig. S2A).
The posterior part of each pars facialis has a narrow, sharp, and posterodorsally oriented process that borders an anterodorsal edge of the orbital aperture (Fig. S2A). The posterior process of the maxillary part forms the posteroventral border of the orbital aperture. There is a wide ventral gutter-like anterolaterally oriented groove between the sharp process of pars facialis and the posterior process that is confluent with the orbital aperture and corresponds to the tentacular canal (therefore, the tentacular canal is an open groove, but not a closed canal separated from the orbital aperture) (Fig. S2B).
The pars palatina of each maxillary part is a wide ledge that contacts the vomer (at the level of the anterior process and the pars facialis of the maxillary part, anterior to the palatine part of the bone) and fuses to the palatine part of the bone posteriorly. The maxillary part bears a labial tooth row (with 13 pedicellate teeth on the left and 15 teeth on the right bone).
The palatine part of each maxillopalatine has the anterior and medial processes that form the anterior, lateral, and posterior borders of the choana and the posterior process for the articulation with the pterygoid posteromedially (Fig. S2C, D). The end of the anterior process of the palatine part is dorsally curved and contacts with a similarly curved crest on the dorsal surface of the vomer; these two structures border a large concavity dorsally. According to
Septomaxillae. The paired septomaxillae are small bony elements that form the lateral borders of the nostrils (Fig. S3). Each septomaxilla contacts the nasal dorsally and posterodorsally, the premaxilla anteroventrally, the maxillopalatine posteroventrally, and has a narrow contact with the prefrontal posteriorly. The anterior portion of each bone is rounded in the anterior view and has shaped like an open tube, with dorsal and ventral processes (Fig. S3A). The anterior dorsal process has a ledge-like medial portion. The posterior part of each bone is triangular in lateral view. The ventral edge of each bone is slightly inflected medially. The lateral surface bears numerous neurovascular foramina that are connected with a longitudinal canal on the medial surface of the posterior part of the bone; this canal is open medially and looks like a deep longitudinal groove (Fig. S3B). The neurovascular foramina and the corresponding canal/groove are presumably for the passing of the branch of the ophthalmic division of the trigeminal nerve and the associated blood vessels.
Nasals. Each nasal is a flat and anteroposteriorly elongated bone (the ratio of its maximum width to the midline length is about 0.44) (Fig. S4). Each nasal is a large bone nearly the same size as the parietal. Each bone is trapezoidal, with a tapering posterior process. Each nasal contacts the premaxilla anteriorly, the septomaxilla anterolaterally, the prefrontal laterally, and the frontal posterolaterally. The paired nasals articulate with each other medially along their entire length. Anteriorly, the nasals contribute to the dorsal margin of the nostril. The dorsal surface of each bone bears several neurovascular foramina in the anterior and central portions, presumed to be for branches of the ophthalmic division of the trigeminal nerve and the associated blood vessels (Fig. S4B).
Frontals. Each frontal is a flat bone, longer than wide (the ratio of its maximum width to the maximum length is about 0.46), with its maximum width approximately at the level of mid-length (Fig. S5). Each frontal is somewhat smaller than the nasal and the parietal. The anterior portion is triangular in dorsal and ventral views, with the top oriented anterolaterally. The posterior portion is blunt and rounded in dorsal and ventral views. On the dorsal surface, the anteromedial edge of the anterior portion bears a wide facet for contact with the nasal (Fig. S5B, nf), and the anterolaterally edge of the anterior portion bears a facet for contact with the prefrontal (Fig. S5B, C). On the ventral surface, the anteromedial edge bears a deep groove-like indentation, forming the dorsal part of the canal presumably for the dorsal olfactory nerve (Fig. S5A). Each frontal contacts the nasal anteromedially, prefrontal anterolaterally, and the parietal posteriorly. The paired frontals contact each other along about one-third of their lengths; they are widely separated by the posterior process of the nasal anteriorly and by the anteromedial portions of the parietals posteriorly. The lateral edge of each bone forms the medial border of the temporal fossa. The ridge (Fig. S5C), formed by the lateral edge of each bone, presumably acts as an attachment point for the long adductor of the lower jaw (m. adductor mandibulae longus) (see
Parietals. Each parietal is a large trapezoidal and flat bone (Fig. S6); its lateral edge is slightly curved ventrally; parietals are the largest bones of the cranial roof. Each parietal is anteroposteriorly elongated (the ratio of its maximum width to the midline length is about 0.48), with the maximum width approximately at the level of the posterior third of its length. Each parietal is slightly larger than the nasal. Each parietal contacts the frontal anteriorly, the os basale laterally, and overlaps the tip of the dorsomedial process of the sphenethmoid anteromedially and the os basale posteriorly. The paired parietals contact with each other medially along their entire lengths.
The anterior part of the dorsal surface of each bone has a wide, semicircular facet for the overlying frontal (Fig. S6B). Posteriorly to the frontal facet, the dorsal surface divides into the flat medial area, tapering posteriorly, and the slightly curved lateral area. The anterolateral edge of the flat medial area bears a rounded ridge. This ridge is a direct continuation of the similar ridge on the frontal, extending the presumed attachment point for the long adductor of the lower jaw. The lateral curved area is presumed to be for the attachment of the m. depressor mandibulae (
Squamosals. Each squamosal is a large (about the same size as the parietal), trapezoidal, and slightly convex bone that forms most of the lateral surface of the temporal region of the cranial roof (Fig. S7). Each squamosal contacts the circumorbital anterodorsally, the maxillopalatine anteroventrally, and overlaps the quadrate posteroventrally. The dorsal edge of each squamosal forms the lateral border of the temporal fossa. The anteroventral part of the lateral surface of each bone has a triangular socket-like facet for contact with the posterior process of the maxillopalatine. The posterior part of the lateral surface has a wide depression, presumably for the attachment of the m. depressor mandibulae (
Prefrontals. Each prefrontal is a flat and relatively small bone (slightly larger than the septomaxilla), elongated anteroposteriorly (the ratio of its maximum width to the midline length is about 0.3) (Fig. S8). Each prefrontal contacts the septomaxilla anteriorly, the nasal anterodorsally, the frontal dorsally and posterodorsally, and the maxillopalatine anteroventrally, ventrally, and posteroventrally. The ventral edge of each bone is thickened and has a facet for contact with the pars facialis (namely, an elongated sharp process) of the maxillopalatine. On the left side of the cranium, each prefrontal appears to be fused with this process.
Circumorbitals. Each circumorbital is a small, crescent-shaped, and curved bone, widely open ventrally, that forms the dorsal and posterodorsal border of the orbital aperture (Fig. S9); it was referred to as postfrontal in some earlier works (
Vomers. Each vomer is a large (about the same size as the nasal) and flat dentigerous bone that forms most of the anterior part of the palate (Fig. S10). Each vomer is subtriangular in the ventral view, with the maximum width at the level just behind the vomerine foramen. Each vomer contacts the premaxilla anteriorly, the maxillopalatine anterolaterally and laterally, the sphenethmoid dorsally, and overlaps the os basale posteriorly. The paired vomers articulate with each other medially along their entire length. Each vomer forms the medial border of the choana. The posterior edge of each vomer is situated slightly anteriorly to the center of the palate (about 44–45% of the cranium length).
The ventral surface of each vomer bears an anterior part of the lingual (= vomero-palatine) tooth row (with 10 functional teeth on the left and 11 on the right bone) that is curved along the anterolateral edge of the bone and is parallel to the labial tooth row. There is a large vomerine foramen in the central part of the ventral surface (Fig. S10A).
The dorsal surface of each vomer is complex. The anterior part bears an anterodorsally oriented premaxillary process that is in articulation with the pars dorsalis of the premaxilla (Fig. S10B). The base of the process is perforated by a longitudinal canal presumably for the ophthalmic division of the trigeminal nerve and the associated blood vessels that continue into the corresponding canal in the premaxilla. Posteriorly, there is a longitudinal groove that runs parallel to the medial line and almost reaches the level of the middle of each bone. Laterally and posterolaterally to this groove, the dorsal surface has a bulge. This bulge is limited posterolaterally by the curved crest. The bulge contains the inner vomerine cavity that opens ventrally through the vomerine foramen and posteriorly and posterolaterally through the neurovascular canals. The curved crest separates a large concavity from the bulge (Fig. S10B). As mentioned above, this concavity is dorsally bordered by the anterior process of the palatine part of the maxillopalatine and the curved crest on the dorsal surface of the vomer and presumably contains the ‘Choanenschleimbeutal’ sensu
Pterygoids. Each pterygoid is a small and anteroposteriorly elongated bone (Fig. S11). Each pterygoid contacts with the maxillopalatine anterolaterally and with the quadrate posterodorsally. Each pterygoid forms the lateral border of the mediopalatinal cavities (= interpterygoid vacuities) and the medial border of the lower temporal fossae. Each pterygoid consists of the anteromedially oriented anterior process, the dorsomedially oriented basal process, and the posteroventrally oriented posterior process (Fig. S11A). In lateral view, the ventral edge of the posterior process is situated barely below the level of the tips of the teeth of the labial (premaxillary-maxillary) row. The basal process and the posterior process comprise the expanded and subvertical portion of each bone. The anterior process is tapering anteriorly, and it contacts the posterior portion of the palatal part of the maxillopalatine anterolaterally. The basal process is low, and its posterodorsal part contacts the quadrate. There is no pronounced contact between the basal process of each pterygoid and the os basale.
Quadrates. Each quadrate is a relatively small bone that forms a joint with the lower jaw (= the articular facet for the pseudoangular) (Fig. S12). Each quadrate contacts with the pterygoid anteroventrally and is largely overlapped by the squamosal dorsolaterally. Each quadrate consists of two distinct parts: the condylar portion and the moderately elongated processus ascendens (sensu
Sphenethmoid. The sphenethmoid consists of the laterally expanded main body, the elongated anteriorly projecting nasal septum, and the posteriorly projecting dorsomedial process (Fig.
Os basale. The os basale consists of the main body, formed by the fusion of the occipital and otic bones, the ossified antotic walls, and the ventral part formed by the parasphenoid (Fig.
The antotic wall has a weak incisure presumably for the ophthalmic nerve on its anterior edge and a large subcircular foramen in its posterior portion, through which the branches of the trigeminal nerve presumably exit the braincase (Fig.
The medial wall (= inner surface) of the otic capsule region bears a number of foramina. The anterior-most foramen is presumed to be for the trunk of the facial nerve (Fig.
Stapes. Each stapes is a small bone that overlaps the fenestra ovalis (Fig. S13). Each stapes consists of a footplate and a columellar process (= style) (Fig. S13A). The footplate is oval and dorsoventrally compressed. There is a small foramen perforating the footplate dorsally to the base of the columellar process. The columellar process is anterodorsally oriented and extends in the direction of the processus oticus of the quadrate (but does not contact it). The columellar process ends in a depression indicating that in life its distal-most end was finished in cartilage. The base of the columellar process is dorsoventrally expanded, forming a robust circular ridge (Fig. S13C, E). The stapedial foramen is absent; however, there is a deep groove formed by the ventral part of the footplate and the ridge of the base of the columellar process (Fig. S13E); the stapedial artery presumably passes through this groove.
Pseudodentaries. The paired elongated dentigerous pseudodentaries articulate with each other medially in the symphyseal area to form the anterior tip of the mandible (Fig. S14). The posterior half of each bone is overlapped by the pseudoangular medially (lingually) (Fig. S14C). Each bone is the widest in the anterior presymphyseal part (in dorsal and ventral views). Each bone has a relatively straight occlusal margin and a nearly similar depth along its entire length except for the tapering posterior end. Each pseudodentary has a pronounced ‘splenial’ ridge (sensu
The labial outer mandibular tooth row is longer than the lingual inner mandibular tooth row and terminates at the level of the posterior end of the ‘splenial’ ridge. There are numerous neurovascular foramina (Fig. S14A) that presumably represent the openings of the inner canals for the branches of the mandibular division of the trigeminal nerve (external and alveolar branches) and the alveolar branch of the facial nerve (
The external branch of the mandibular division of the trigeminal nerve (r. mandibularis externus) presumably enters the pseudodentary through two large foramina on the lingual side posterior to the level of the posterior ends of the outer and inner mandibular tooth rows. After passing through the bone, it exits through numerous foramina at the base of the teeth on the lingual side and along the entire length of the pseudodentary on the labial side (Fig. S14A, C). The combined alveolar branches of the mandibular division of the trigeminal nerve and of the facial nerve presumably pass in a groove below the facet for contact with pseudoangular.
Pseudoangulars. Each pseudoangular is an elongate bone that overlaps the pseudodentary medially (lingually) (Fig. S14C). Each pseudoangular consists of the anterior process, the processus condyloideus, the processus internus, and the retroarticular process. The anterior process is tapering anteriorly. The processus condyloideus (= condyle) is situated on the dorsal margin of the bone, anterior to the retroarticular process, and forms a joint with the cranium (= the articular facet for the quadrate) (Fig. S14A). The processus internus is a small lingually oriented projection on the lingual surface of the bone, ventrally to the processus condyloideus (Fig. S14B). Anteriorly to the processus condyloideus, the dorsal margin of each bone has a high curved ridge and a high dorsal projection. The posterior end of the curved ridge terminates between the processus condyloideus and the processus internus, and the anterior end of this ridge is fused with the dorsal projection. The anterior end of the curved ridge and the dorsal projection border a large foramen dorsolaterally (Fig. S14B). This is the foramen of the canalis primordialis presumably for the passing of the mandibular division of the trigeminal nerve and the mandibular artery. Ventrally to the foramen of the canalis primordialis, the labial surface of the bone bears a large foramen presumably for the passage of the intermandibular branch of the mandibular division of the trigeminal nerve (Fig. S14C).
The retroarticular process is elongated, with its rounded posterior end oriented nearly dorsally (Fig. S14F). At the base of the retroarticular process, posterior to the processus internus, there is a large oblique foramen presumably for the passage of the alveolar branch of the facial nerve (Fig. S14F).
Atlas. The atlantal centrum is very short (Fig.
The atlantal neural arch is high. In the lateral view, the anterior edge of the neural arch is oblique; its dorsal part projects anteriorly beyond the level of the anterior edge of the anterior cotyles. The neural spine is not well pronounced, but there is a very low and swollen median ridge. The anterior part of the dorsal surface of the neural arch is smooth. The posterior part of the dorsal surface of the neural arch is lower than the anterior part and bears a low median depression. The postzygapophyseal processes are oriented ventrolaterally (at approximately 20° from the midline) and are elongate and narrow in ventral outline. The neural canal is broad, high, and nearly round in its posterior outline. The inner surface of the base of the neural arch bears a pair of spinal cord supports that extend medially into the neural canal (Fig.
Sequence variation. A total of 3967 aligned base pairs (1143 bp from cyt b, and 2824 bp from the fragment containing 12S rRNA and 16S rRNA genes) were obtained for our study. Sequence characteristics, including the estimated transition/transversion bias, nucleotide frequencies, and suggested best-fit models of DNA evolution for each genetic marker and partition, are summarized in Table S3. We deposited the newly obtained sequences in GenBank (see Table S2 for the accession numbers and voucher details; the final alignment is available upon request).
Phylogenetic relationships. ML and BI analyses of our mtDNA alignment recovered trees with identical topologies, except for a few relationships that do not affect our interpretations of the results (Fig.
Maximum Likelihood (ML) phylogenetic tree of the genus Ichthyophis derived from the analysis of 3967 bp of cyt b, 12S rRNA, and 16S rRNA mtDNA gene sequences. For voucher specimen information and GenBank accession numbers, see Table S2. Numbers at the nodes correspond to ML UFBS/BI PP support values, respectively. Species symbols and locality numbers correspond to those provided on the map in Figure
Intragenerically, our analyses revealed four deeply-diverged clades with unclear phylogenetic relationships (see Fig.
Clade D (99/1.0) was well-supported and regrouped Ichthyophis taxa from East and Southeast Asia (including northeastern India). It was further subdivided into three major subclades (Fig.
Pairwise distances. The uncorrected p distances for the cyt b gene and the 16S rRNA mtDNA fragment among the members of the genus Ichthyophis are summarized in Table S4. The interspecific distances among Ichthyophis species in the cyt b sequences varied from p = 5.3% (between I. supachaii and I. hypocyaneus) to p = 19.8% (between I. elongatus and I. cf. asplenius, and I. elongatus and I. pauli), while in the 16S rRNA sequences, p distances varied from p = 1.5% (between I. supachaii and I. hypocyaneus) to p = 9.9% (between I. orthoplicatus and Ichthyophis sp. A). The genetic divergence of the newly discovered population of Ichthyophis sp. from northern Vietnam in cyt b gene sequences varied from p = 6.5% (with I. yangi) to p = 14.6% (with Ichthyophis sp. E), while in the 16S rRNA gene sequences p distances varied from p = 4.5% (with I. chaloensis) to p = 7.6% (with I. khumhzi, Ichthyophis sp. D, Ichthyophis sp. E). Inter-group genetic differentiation between the two samples of Ichthyophis sp. from northern Vietnam was small and comprised p = 0.9% (for the 16S rRNA gene fragment); the two samples of Ichthyophis sp. shared the same cyt b haplotype.
Our molecular phylogenetic analyses strongly suggested that the newly discovered caecilian populations from Xuan Lien NP and Pu Hoat NR in northern Vietnam belong to the genus Ichthyophis, within which they formed, among the sampled species, a distinct species-level lineage, which was closely related to two recently discovered unstriped species: I. yangi from Yunnan Province, China, and I. chaloensis from central Vietnam (Fig.
Overall, based on the external morphology and geographic proximity of their distributions, the unstriped species I. chaloensis, I. yangi, and Ichthyophis sp. from northern Vietnam appeared to be most closely related to the enigmatic species I. laosensis, described from central Laos and to date reliably known only from the holotype
Morphological examination of the two specimens of Ichthyophis sp. from northern Vietnam indicated the presence of several diagnostically important traits related both to body proportions and to meristic characters, including the number of annuli, teeth, and vertebrae, which allowed distinguishing this population from all other unstriped species of Ichthyophis in mainland Southeast Asia, as well as from all other congeners (summarized in Table 1). As a basis for further taxonomic research on Southeast Asian Ichthyophiidae, we also provided a detailed description of osteological traits of this population, including a detailed morphological description of skull bones and the atlas. Overall, our results supported the hypothesis that the recently discovered population of Ichthyophis sp. from northern Vietnam represents an undescribed species, which we formally describe below.
Comparison of morphological characters of Ichthyophis griseivermis sp. nov. with other unstriped species of the genus Ichthyophis from Indochina and China. Diagnostic differences from the new species are marked in bold. Morphological data from larval specimens were not used in comparisons with the new species type series. For character abbreviations see Materials and methods. Type status abbreviations: H – holotype; P – paratype. Morphological data are taken from
| Species | Ichthyophis griseivermis sp. nov. | I. yangi | I. catlocensis | I. chaloensis | I. laosensis | |||
| Museum IDs |
A-8208 |
VRTC NAP-08953 |
|
2023048 |
|
88976 |
IEBR A.2011.16 |
|
| Type status | H | P | H | P | P | H | P | H |
| Sex | F | M | F | F | F | F | F | F |
| Morphometry (in mm) | ||||||||
| TL | 206.0 | 242.0 | 321.4 | 307.4 | 329.5 | 183.5 | 215.7 | 318.0 |
| TAL | 2.4 | 3.6 | 2.9 | 3.0 | 3.9 | 2.1 | 3.7 | 3.7 |
| BW1 | 6.6 | 6.9 | 10.6 | 10.5 | 8.0 | 5.3 | 6.9 | 11.4 |
| BW2 | 9.8 | 11.2 | 13.9 | 13.8 | 10.5 | 7.1 | 7.6 | 16.1 |
| BW3 | 4.9 | 5.0 | 5.8 | 5.7 | 4.2 | 3.7 | 3.2 | 6.3 |
| HL | 8.8 | 9.8 | 14.6 | 15.4 | 17.0 | 7.3 | 9.4 | 10.9 |
| HW | 6.1 | 6.9 | 8.7 | 8.1 | 9.5 | 4.8 | 6.3 | 9.7 |
| UJL | 6.5 | 8.0 | 10.8 | 9.6 | 10.5 | 5.4 | 7.2 | 10.8 |
| LJL | 6.2 | 7.5 | 10.3 | 9.2 | 10.1 | 5.2 | 6.3 | 10.3 |
| SP | 0.6 | 0.8 | 0.8 | 0.6 | 0.7 | 0.5 | 1.3 | 0.6 |
| EL | 0.7 | 0.8 | 0.8 | 0.7 | 1.0 | 0.6 | 0.8 | 1.7 |
| ES | 4.0 | 4.7 | 5.7 | 5.5 | 5.6 | 3.6 | 4.8 | 6.3 |
| EE | 4.6 | 5.8 | 5.9 | 5.1 | 6.6 | 4.3 | 4.8 | 7.3 |
| NN | 1.8 | 2.2 | 2.0 | 1.7 | 2.6 | 2.0 | 1.8 | 3.2 |
| TT | 4.2 | 4.9 | 6.1 | 5.5 | 5.4 | 4.2 | 4.7 | 7.5 |
| EN | 2.8 | 3.7 | 4.2 | 4.3 | 5.4 | 2.9 | 3.6 | 5.0 |
| ET | 1.1 | 1.5 | 1.2 | 1.1 | 1.2 | 0.6 | 1.3 | 1.6 |
| TN | 2.4 | 3.1 | 1.3 | 1.3 | 1.3 | 2.7 | 2.8 | 4.3 |
| ED | 0.6 | 0.7 | 1.0 | 1.1 | 1.2 | 0.5 | 0.3 | 1.0 |
| ETH | 0.7 | 0.9 | – | – | – | – | – | – |
| C1 | 2.2 | 2.8 | – | – | – | – | – | – |
| C2 | 2.4 | 3.0 | – | – | – | – | – | – |
| STTA | 3.6 | 4.5 | – | – | – | – | – | – |
| LTA | 0.6 | 0.7 | – | – | – | – | – | – |
| CD | 1.5 | 3.0 | – | – | – | – | – | – |
| Measurement ratios | ||||||||
| TN/ET | 2.2 | 2.1 | 1.1 | 1.1 | 1.1 | 4.5 | 2.2 | 2.7 |
| HL/ED | 15.5 | 14.5 | 15.2 | 13.8 | 14.2 | 14.6 | 31.3 | 10.9 |
| SP/HL | 0.06 | 0.08 | 0.05 | 0.04 | 0.04 | 0.07 | 0.14 | 0.06 |
| TL/TAL | 85.5 | 66.7 | 110.8 | 103.9 | 83.6 | 87.4 | 58.3 | 85.9 |
| TL/BW2 | 21.1 | 21.6 | 23.2 | 22.4 | 31.4 | 25.8 | 28.4 | 19.8 |
| ES/HL | 0.46 | 0.48 | 0.39 | 0.36 | 0.33 | 0.49 | 0.51 | 0.58 |
| EL/HL | 0.08 | 0.08 | 0.05 | 0.04 | 0.06 | 0.08 | 0.09 | 0.16 |
| Meristic characters | ||||||||
| TAD | 306 | 301 | 369 | 372 | 372 | 342 | 344 | 346 |
| TAV | 298 | 292 | 367 | 369 | 368 | 340 | 342 | 345 |
| AV | 4 | 4 | 6 | 6 | 6 | 5 | 3 | 3 |
| TAT | 1 | 1 | 2 | 2 | 2 | 5 | 5 | 2 |
| PMM | 44 | 48 | – | 51 | 53 | 44 | 37 | 33 |
| VP | 43 | 48 | – | 44 | 44 | 51 | 54 | 36 |
| DE | 38 | 38 | – | 48 | 49 | 27 | 26 | 35 |
| IM | 25 | 33 | – | 36 | 35 | 16 | 11 | 30 |
| VERT | 112 | 111 | 114 | 114 | – | 110 | 110 | 112 |
| TG | 0 | 0 | – | – | – | – | – | – |
| Species | I. youngorum | I. acuminatus | I. cardamomensis | |||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Museum IDs |
|
|
|
|
|
|
|
|
|
|
|
|
| Type status | H | P | P | P | H | P | P | P | P | H | P | P |
| Sex | M | M | – (larva) | F (larva) | M | F | F | F | – (larva) | F | F | F |
| Morphometry (in mm) | ||||||||||||
| TL | 208.0 | 217.0 | 82.0 | 189.0 | 294.0 | 199.0 | 203.0 | 215.0 | 172.0 | 183.0 | 321.7 | 289.1 |
| TAL | 2.9 | 3.1 | 1.8 | 3.5 | 3.9 | 2.8 | 2.4 | 3.6 | 2.0 | 3.1 | 3.8 | 6.0 |
| BW1 | 8.1 | 8.0 | – | 7.8 | 10.6 | 7.8 | 6.4 | 8.1 | 6.7 | 7.1 | 9.2 | 8.9 |
| BW2 | 11.2 | 10.5 | 4.8 | 10.6 | 13.7 | 10.9 | 8.1 | 10.6 | 8.0 | 7.9 | 8.6 | 8.0 |
| BW3 | 5.5 | 5.0 | 2.5 | 5.3 | 5.4 | 4.3 | 3.6 | 5.9 | 3.9 | 4.1 | 5.4 | 4.9 |
| HL | 10.4 | 10.4 | – | 9.3 | 14.6 | 8.8 | 9.1 | 10.1 | 8.8 | 8.1 | 11.7 | 9.2 |
| HW | 7.6 | 7.3 | 3.9 | 7.0 | 9.8 | 6.7 | 6.2 | 7.0 | 6.3 | 5.8 | 7.9 | 7.1 |
| UJL | 7.6 | 7.3 | 2.8 | 4.5 | 10.6 | 6.4 | 6.2 | 6.2 | 5.9 | 6.9 | 10.2 | 9.5 |
| LJL | 6.7 | 6.9 | 2.2 | 4.5 | 9.3 | 6.0 | 5.6 | 5.6 | 5.9 | 6.0 | 9.7 | 9.0 |
| SP | 0.1 | 0.1 | 0.3 | 0.7 | 0.8 | 0.4 | 0.1 | 0.3 | 0.4 | 0.7 | 0.5 | 0.6 |
| EL | – | – | – | – | – | – | – | – | – | 1.0 | 1.0 | – |
| ES | – | – | – | – | 5.8* | – | – | – | – | 4.1 | 6.0 | 5.5 |
| EE | 4.5 | 4.8 | 2.4 | 4.2 | 6 | 4.3 | 4 | 4.6 | 4.3 | 4.3 | 5.6 | 5.2 |
| NN | 2.4 | 2.2 | 1.4 | 1.7 | 2.6 | 1.3 | 1.1 | 1.4 | 1.3 | 2.2 | 2.8 | 2.6 |
| TT | 5.9 | 4.9 | – | – | 7.2 | 5.0 | 4.6 | 5.5 | 5.0 | 4.7 | 6.2 | 5.8 |
| EN | 3.8 | 3.6 | 1.7 | 2.8 | 4.8 | 3.1 | 3.2 | 3.2 | 3.1 | 3.5 | 4.7 | 4.2 |
| ET | 1.1 | 1.0 | – | – | 1.3 | 0.7 | 0.8 | 0.8 | 0.8 | 0.9 | 1.2 | 1.1 |
| TN | 2.7 | 2.5 | – | – | 3.1 | 2.0 | 2.0 | 2.1 | 2.5 | 2.9 | 3.5 | 3.1 |
| ED | – | – | – | – | 0.7* | – | – | – | – | 0.7 | 1.0 | 0.8 |
| ETH | – | – | – | – | – | – | – | – | – | – | – | – |
| C1 | – | – | – | – | – | – | – | – | – | – | – | – |
| C2 | – | – | – | – | – | – | – | – | – | – | – | – |
| STTA | – | – | – | – | – | – | – | – | – | – | – | – |
| LTA | – | – | – | – | – | – | – | – | – | – | – | – |
| CD | – | – | – | – | – | – | – | – | – | – | – | – |
| Measurement ratios | ||||||||||||
| TN/ET | 2.5 | 2.5 | – | – | 2.4 | 2.9 | 2.5 | 2.6 | 3.1 | 3.2 | 2.9 | 2.8 |
| HL/ED | – | – | – | – | 20.9 | – | – | – | – | 11.6 | 13.2 | 11.5 |
| SP/HL | 0.01 | 0.01 | – | 0.08 | 0.05 | 0.05 | 0.01 | 0.03 | 0.05 | 0.09 | 0.04 | 0.06 |
| TL/TAL | 71.7 | 70 | 45.6 | 54.0 | 43.0 | 71.1 | 84.6 | 59.7 | 86 | 59 | 84.7 | 48.1 |
| TL/BW2 | 18.6 | 20.7 | 17.1 | 17.8 | 21.5 | 18.3 | 25.1 | 25.1 | 21.5 | 23.2 | 37.4 | 36.1 |
| ES/HL | – | – | – | – | 0.40 | – | – | – | – | 0.51 | 0.51 | 0.60 |
| EL/HL | – | – | – | – | – | – | – | – | – | 0.12 | 0.09 | – |
| Meristic characters | ||||||||||||
| TAD | 318 | 328 | – | 317 | 302 | 296 | – | 301 | 327 | 322 | 364 | 340 |
| TAV | – | – | – | – | – | – | – | – | – | 320 | 359 | 338 |
| AV | 5 | 6 | – | 6 | 6 | 7 | 5 | 6 | 3 | 5 | 4 | 5 |
| TAT | 4 | 4 | – | 4 | 3 | 2 | 3 | 2 | 3 | 6 | 4 | 6 |
| PMM | 28 | 22** | – | 23** | 40 | 43 | 37 | 41 | 35 | 23 | 38 | 38 |
| VP | 33** | 40 | – | 33** | 46 | 39 | 38 | 43 | 37 | 29 | 28 | 28 |
| DE | 28** | 29** | – | 28** | 42 | 43 | 38 | 42 | 31 | 20 | 34 | 34 |
| IM | 19** | 18** | – | 15** | 32 | 26 | 26 | 26 | 19 | 19 | 21 | 22 |
| VERT | 108 | 108 | 104 | 108 | 109 | 108 | – | 111 | 111 | 120 | – | 120 |
| TG | – | – | – | – | – | – | – | – | – | – | – | – |
| * Based on |
||||||||||||
VRTC NAP08953 (field tag NAP-08953), adult male, from Pu Hoat Nature Reserve, Tien Phong Commune, Nghe An Province, Vietnam (elevation 815 m a.s.l.; geographic coordinates: 19.755°N, 104.796°E), collected by N. A. Poyarkov from the bank of a river on May 15, 2019.
The specific name “griseivermis” is a Latin noun in the nominative singular, given in apposition, derived from the Latin adjective “griseus” for “grey” and the Latin noun “vermis” for “worm.” The new species is named in reference to its characteristic uniform grey body coloration. The specific epithet also alludes to Grey Worm, the commander of the Unsullied, the warrior-eunuchs of Astapor with an unparalleled reputation for combat in George R. R. Martin’s fictional work “A Song of Ice and Fire” (also known as “Game of Thrones”). We suggest the following common names for the new species: “Grey Worm Caecilian” (in English), “?ch giun xám B?c Trung B?” (in Vietnamese), and “Seryi rybozmey” (“????? ????????,” in Russian).
The new species Ichthyophis griseivermis sp. nov. differs from other members of the genus Ichthyophis by the following combination of the morphological characters: unstriped body lacking lateral yellow stripe; adult total length 206–242 mm (based on two available specimens); snout blunt and rounded (snout length/head length ratio 0.06–0.08); tentacle aperture located closer to eye than to naris (tentacle aperture-naris distance/tentacle aperture-eye distance ratio 2.1–2.2); premaxillary and maxillary teeth 44–48, vomero-palatine teeth 43–48, dentary teeth 38, inner mandibular teeth 25–33; tail very short, acuminate, ending in a nipple-like cap; annuli angulate, total 301–306 (dorsal count), four interupted by cloacal disc, one posterior to cloacal disc, the degree of annuli angulation decreasing from head to cloaca with grooves appearing almost orthoplicate at mid-body and posteriorly; vertebrae 111–112; scales in one series per annulus (dosolaterally), present only in the posterior half of body.
Adult female (Figs
The holotype of Ichthyophis griseivermis sp. nov. in life (
Details of external morphology of the holotype of Ichthyophis griseivermis sp. nov. in life (
In life (Figs
Variation in measurements and meristic characters of the type series is presented in Table 1. The paratype (VRTC NAP08953, adult male) was found as a desiccated specimen and is in a moderate condition of preservation (Fig. S16). Paratype body significantly dorsoventrally flattened; a large transverse incision present in the posterior one-third of the specimen length on the ventral side. The copulatory organ was extruded from the cloaca and damaged prior to specimen collection; the remains of the everted organ are visible on the ventral and lateral aspects of the specimen (Fig. S16A, B), but do not allow us to provide a description of its morphology. Overall, the male paratype VRTC NAP08953 is generally very similar in most morphological traits to the holotype; however, it is c. 17% longer than the holotypeis (TL 242 mm), and has a slightly lower number of annuli (TAD 301, TAV 292), one less vertebra (VERT 111), and a greater number of teeth in all tooth series (PMM 48; VP 48; DE 38; IM 33) all within the range of intraspecific variation expected for an Ichthyophis. It is identical to the holotype in in the number of annuli interrupted by the cloacal disc (AV 4) and the number of annuli posterior to the cloacal disc (TAT 1); and similar in body proportions. Due to desiccation and storage in ethanol for over five years, the original dark coloration of the paratype has significantly faded; the specimen is uniformly brown; the head and ventral surfaces are light-brown (Fig. S16).
The new species lacks light lateral stripes, so it can be easily distinguished from all striped members of the genus Ichthyophis, and its comparisons with the unstriped congeners are the most pertinent. We will first compare Ichthyophis griseivermis sp. nov. with seven currently recognized unstriped species of the genus Ichthyophis from the Indochinese region (including Vietnam, Cambodia, Laos, and Thailand) and China; the main diagnostic characters separating the new species from these species are summarized in Table 1.
Ichthyophis griseivermis sp. nov. is a comparatively small-sized species (TL 206.0–242.0 mm): though only two specimens of the new species are known to date and it is impossible to be confident about its maximal body size, the presence of mature ova in the holotype at least indicates that this specimen is adult at the TL of 206.0 mm. This can arguably distinguish the new species from its sister species I. yangi (endemic to Yunnan Province, China) (TL 307.4–329.5 mm) and, with less confidence, from I. laosensis (known only from northern Laos) (TL of the only known holotype 318.0 mm).
In body proportions, by the relatively longer eye-snout distance (ES/HL ratio 0.46–0.48), the new species can be further distinguished from I. yangi, in which the snout is much shorter (ES/HL ratio 0.33–0.39), and, though with less confidence, from I. acuminatus (distributed in northwestern Thailand and northwestern Laos) (ES/HL ratio 0.40; however, this comparison should be taken with caution as its calculation is based on the ES measurement by
Differences between the new species and its congeners observed in meristic characters should be taken with caution due to a small sample size available for our examination. Nevertheless, I. griseivermis sp. nov. has notably fewer total annuli both in dorsal count (TAD 301–306) and in ventral count (TAV 292–298) than I. catlocensis (TAD 342; TAV 340), I. chaloensis (TAD 344; TAV 342), I. laosensis (TAD 346; TAV 345), I. cardamomensis (TAD 322–364; TAV 320–359), and I. yangi (TAD 369–372; TAV 367–369). The new species has generally fewer annuli interrupted by the cloacal disc (AV 4) than in I. acuminatus (AV 5–7), I. youngorum (AV 5–7), and I. yangi (AV 6), slightly fewer annuli interrupted by the cloacal disc than I. catlocensis (AV 5), and slightly more than I. chaloensis (AV 3) and I. laosensis (AV 3). By having only a single annulus posterior to the cloacal disc (TAT 1), the new species is distinguished from I. acuminatus (TAT 2–3), I. cardamomensis (TAT 2–6), I. catlocensis (TAT 5), I. chaloensis (TAT 5), and I. youngorum (TAT 4). The new species has a slightly higher number of labial premaxillary-maxillary teeth (PMM 44–48) than I. acuminatus (PMM 37–43), I. cardamomensis (PMM 23–38), I. chaloensis (PMM 37), I. laosensis (PMM 33), and I. youngorum (PMM 22–28), but fewer premaxillary-maxillary teeth than in I. yangi (PMM 51–53). Ichthyophis griseivermis sp. nov. has slightly fewer vomero-palatine teeth (VP 43–48) than I. catlocensis (VP 51) and I. chaloensis (VP 54), but more than in I. laosensis (VP 36), I. cardamomensis (VP 28–29), and I. youngorum (VP 33–40; note that the lower value was obtained from a partially damaged specimen and may be erroneous). The new species has more dentary (labial) teeth (DE 38) than I. catlocensis (DE 27), I. chaloensis (DE 26), and I. youngorum (DE 28–29; note that these values were obtained from partially damaged specimens and may be erroneous), but fewer than I. yangi (DE 48–49). At the same time, the new species has more inner mandibular teeth (IM 25–33) than I. catlocensis (16), I. chaloensis (IM 11), I. youngorum (IM 18–19; note that these values were obtained from partially damaged specimens and may be erroneous), and I. cardamomensis (IM 19–22). Ichthyophis griseivermis sp. nov. further differs from I. cardamomensis by fewer vertebrae (VERT 111–112 vs. 120).
Furthermore, I. griseivermis sp. nov. can be readily diagnosed from the following unstriped species of Ichthyophis which occur outside the Indochinese Region and southern China. In particular, the new species differs from I. lakimi (Sabah, Borneo) by having more premaxillary-maxillary (labial) teeth (PMM 44–48 vs. 16–25) and by having more inner mandibular teeth (IM 25–33 vs. 14) (though the tooth counts for I. lakimi reported by Nishikawa et al. 2012 appear to be much lower than are typically known for Ichthyophis and require a re-examination); from I. billitonensis Taylor, 1965 (Belitung Is., Indonesia) by having more inner mandibular teeth (IM 25–33 vs. 2); and by having more annuli in dorsal count (TAD 301–306 vs. 251–254); from I. dulitensis Taylor, 1960 (Sarawak, Borneo) by the absence of scales in the anterior half of body (vs. present); by having more inner mandibular teeth (IM 25–33 vs. 8); and by the absence of a light marking on the throat (vs. present); from I. glandulosus Taylor, 1923 (Basilan Is., Philippines) by having more annuli in dorsal count (TAD 301–306 vs. 273–286); and by a higher number of vertebrae (VERT 111–112 vs. 102); from I. javanicus Taylor, 1960 (Java Is., Indonesia) by having less annuli both in dorsal count (TAD 301–306 vs. 351), and in ventral count (TAV 292–298 vs. 348); and fewer annuli posterior to the cloacal disc (TAT 1 vs. 10); from I. larutensis (Peninsular Malaysia) by the presence of inner mandibular teeth (vs. absent); and by a slightly higher number of vertebrae (VERT 111–112 vs. 107); from I. mindanaoensis (Mindanao Is., Philippines) by having more inner mandibular teeth (IM 25–33 vs. 16–22); from I. monochrous (Bleeker, 1858) (Borneo Is., Indonesia) by having more annuli in dorsal count (TAD 301–306 vs. 247); by having more inner mandibular teeth (IM 25–33 vs. 8); and by a greater number of vertebrae (VERT 111–112 vs. 103); from I. orthoplicatus (Sri Lanka) by having fewer annuli posterior to the cloacal disc (TAT 1 vs. 7); and by having slightly more inner mandibular teeth (IM 25–33 vs. 18–20); from I. sikkimensis Taylor, 1960 (northeastern India) by having more annuli in dorsal count (TAD 301–306 vs. 276–292); by a greater number of vertebrae (VERT 111–112 vs. 106–108); and by having slightly more inner mandibular teeth (IM 25–33 vs. 18–20); from I. singaporensis Taylor, 1960 (Singapore) by the absence of scales on the anterior half of body (vs. present); by having more annuli in dorsal count (TAD 301–306 vs. 260–273); by having fewer annuli posterior to the cloacal disc (TAT 1 vs. 7); and by having more inner mandibular teeth (IM 25–33 vs. 6–10); from I. sumatranus Taylor, 1960 (Sumatra Is., Indonesia) by the absence of scales on the anterior half of body (vs. present); and by having fewer annuli posterior to the cloacal disc (TAT 1 vs. 7); from I. weberi Taylor, 1920 (Palawan Is., Philippines) by the presence of inner mandibular teeth (vs. absent); and by having fewer both in dorsal count (TAD 301–306 vs. 313–329), and in ventral count (TAV 292–298 vs. 304–322).
Overall, the holotype specimen of I. griseivermis sp. nov. (
Below, we compare the cranial features of I. griseivermis sp. nov. with the few congeners for which skull morphology reconstructions are available via MorphoSource or from previous publications: I. asplenius, I. kohtaoensis, I. tricolor, I. multicolor, I. nguyenorum, and I. supachaii (based on
The new species differs from all other Ichthyophis species examined in having a small, crescent-shaped circumorbital, lacking the ventral portion, and bordering only the upper posterior corner of the orbital aperture (vs. a larger circumorbital with a posterior ventral process). The smooth, gradually tapering retroarticular process with its posterior end oriented dorsally is also a characteristic feature of the new species. Additionally, the new species differs from I. asplenius, I. beddomei, I. glutinosus, I. kohtaoensis, I. larutensis, I. mindanaoensis, I. multicolor, I. nguyenorum, I. nigroflavus, I. singaporensis, I. supachaii, I. tricolor, and I. weberi, but resembles I. sikkimensis in having a weakly zygokrotaphic skull (vs. the stegokrotaphic or weakly stegokrotaphic condition of the skull). Furthermore, the new species differs from I. asplenius, I. kohtaoensis, I. tricolor, I. multicolor, I. nguyenorum, and I. supachaii in lacking the stapedial foramen (vs. present in all these species; for other species the structure of the stapes has not been described or illustrated).
Additionally, I. griseivermis sp. nov. further differs from I. kohtaoensis in having a more like V-shaped cranium in the dorsal view (vs. more like U-shaped), in having the posterior edge of the vomer situated slightly anteriorly to the center of the palate, about 44–45% of the cranium length (vs. the vomer is situated around the center of the palate, about 50% of the cranium length), and by the laterally expanded foramen magnum with widely spaced occipital condyles (vs. dorsoventrally expanded foramen magnum with closely positioned condyles). The new species further differs from I. tricolor in having a laterally open tentacular groove confluent with the orbital aperture (vs. a laterally closed tentacular groove not confluent with the orbital aperture), by having moderately sized teeth with gently posteriorly curved tips (vs. enlarged teeth with strongly curved tips), by the ventral edge of the posterior process of the pterygoid situated barely below the level of premaxillary-maxillary teeth (vs. far below the level of premaxillary-maxillary teeth), and by the absence of the contact between circumorbital and frontal (vs. present). The new species further differs from I. multicolor in having the posterior edge of the vomer situated slightly anteriorly to the center of the palate, about 44–45% of the cranium length (vs. vomer situated around the center of palate, about 50% of the cranium length), and in the absence of the contact between circumorbital and frontal (vs. present). The new species further differs from I. asplenius in having a more rounded tip of the snout (vs. more blunt), in more ventrally located nostrils in lateral view (vs. more dorsally), in having a laterally open tentacular canal confluent with the orbital aperture (vs. partly laterally closed tentacular canal not confluent with the orbital aperture), in rounded choanae (vs. subtriangular), and in the absence of the lateroventral process of the squamosal (vs. present). The new species further differs from I. nguyenorum in having widely spaced occipital condyles in ventral view (vs. almost confluent) and in a comparatively wider anterior part of the parasphenoid portion of os basale (vs. narrow). The new species further differs from I. supachaii in having widely spaced occipital condyles in ventral view (vs. almost confluent) and in having a laterally open tentacular groove confluent with the orbital aperture (vs. partly laterally closed tentacular groove not confluent with the orbital aperture). The new species differs from I. sikkimensis in having a longer contact between the frontal and prefrontal (vs. shorter, L-shaped contact). The new species further differs from I. beddomei in having a more like V-shaped cranium in the dorsal view (vs. more like U-shaped), in lacking the contact between frontal and circumrobital (vs. wide contact), in having oblique and narrow posterior edge of prefrontal in dorsal view (vs. transverse), and in having long anterior process of pterygoid (vs. short).
Ichthyophis griseivermis sp. nov. further differs from I. glutinosus in having a more like V-shaped cranium in the dorsal view (vs. more like U-shaped), in having rounded snout tip in dorsal aspect (vs. blunt), in having a laterally open tentacular groove confluent with the orbital aperture (vs. a laterally closed tentacular canal), in having palatine tooth series terminating at the level of the anterior border of the adductor chamber (vs. extends behind the level of the anterior border of the adductor chamber), in having widely spaced occipital condyles in ventral view (vs. closely positioned), in having shorter midline contact between frontals (vs. longer), and in having oblique and narrow posterior edge of prefrontal in dorsal view (vs. rounded). The new species further differs from I. larutensis in having rounded snout tip in dorsal aspect (vs. blunt), in having a laterally open tentacular groove confluent with the orbital aperture (vs. a laterally closed tentacular canal), in having palatine tooth series terminating at the level of the anterior border of the adductor chamber (vs. terminates anterior to the adductor chamber), in having widely spaced occipital condyles in ventral view (vs. closely positioned), in having shorter midline contact between frontals (vs. longer), and in having long anterior process of pterygoid (vs. short), and in a comparatively wider anterior part of the parasphenoid portion of os basale (vs. narrow tapering). The new species further differs from I. mindanaoensis in having a more like V-shaped cranium in the dorsal view (vs. more like U-shaped), in having rounded choanae in ventral view (vs. distinctly triangular), in having a laterally open tentacular groove confluent with the orbital aperture (vs. a laterally closed tentacular canal), in having shorter midline contact between frontals (vs. longer), and in a comparatively wider anterior part of the parasphenoid portion of os basale (vs. narrow tapering). The new species further differs from I. nigroflavus in having a more like V-shaped cranium in the dorsal view (vs. more like U-shaped), in having rounded choanae in ventral view (vs. anteroposteriorly elongated, oval), in having widely spaced occipital condyles in ventral view (vs. in contact), in having long anterior process of pterygoid (vs. short), and in a comparatively wider anterior part of the parasphenoid portion of os basale (vs. narrow tapering).
Ichthyophis griseivermis sp. nov. further differs from I. singaporensis in having a more like V-shaped cranium in the dorsal view (vs. more like U-shaped), in having rounded snout tip in dorsal aspect (vs. blunt), in having rounded choanae in ventral view (vs. subtriangular), in having a laterally open tentacular groove confluent with the orbital aperture (vs. a laterally closed tentacular canal), in having palatine tooth series terminating at the level of the anterior border of the adductor chamber (vs. extends behind the level of the anterior border of the adductor chamber), in having widely spaced occipital condyles in ventral view (vs. in contact), and in a comparatively wider anterior part of the parasphenoid portion of os basale (vs. narrow tapering). The new species further differs from I. weberi in having a more like V-shaped cranium in the dorsal view (vs. more like U-shaped), and in having long anterior process of pterygoid (vs. short).
Finally, a scan of the skull of the I. laosensis specimen from Laos (
Currently, I. griseivermis sp. nov. is known from two protected areas in northern Vietnam: from Xuan Lien NP in Thanh Hoa Province and Pu Hoat NR in Nghe An Province (Fig.
The type locality is at an elevation of 800 m asl. The holotype was found on the bank of a small forest stream with a stony bottom and steep clay bank (Fig. S17). The surrounding secondary montane evergreen forest is severely damaged by regular logging, but some old trees remain, including Cunninghamia konishii Hayata (Cupressaceae), forming mixed polydominant forests with Symingtonia populnea (R. Br. ex Griff.) (Hamamelidaceae), Carallia suffruticosa Ridl. (Rhizophoraceae), Engelhardtia roxburghiana Wall (Juglandaceae), Guarea excelsa Kunth (Meliaceae), Castanopsis ferox (Rosb.) (Fagaceae), Michelia mediocris Dandy (Magnoliaceae), Pellionia radicans var. grande (Gagnep.) H. Schroter (Urticaceae), Ardisia quinquegona Blume (Primulaceae), Litsea acutivena Hayata and Litsea yunnanensis Y. C. Yang & P. H. Huang (Lauraceae), and Alniphyllum fortunei (Hemsley) Makino (Styracaceae). The undergrowth is well developed; soil is densely covered with grasses, sedges, and fern thickets. The holotype was found in a small ravine at the foothill of a low mountain range composed of clays and shales with a limestone base. The holotype was found at night (ca. 23:00 h) while slowly crawling among the stones on the banks of a small mountain stream during a rain with an ambient air temperature around 19°C. The holotype was observed crawling along the unflooded part of the bank covered with silt and large pebbles about 30 cm from the water’s edge. The paratype male was found in the daytime (10:00 h) as a dead, desiccated specimen on a sandy bank of a river among large stones at an elevation of ca. 815 m asl. The two localities are separated by a direct distance of 31 km.
Though we don’t have direct observations on the new species reproductive biology, we can assume that, like the other members of the genus Ichthyophis, I. griseivermis sp. nov. is oviparous with aquatic larvae, which is also supported by the large size of the eggs in the new species. Like all caecilians, the new species is carnivorous; however, we do not have specific information about its diet. Other species of amphibians recorded at the same habitat in the Xuan Lien NP and in Pu Hoat NR included: Boulenophrys palpebralespinosa (Bourret, 1937), Leptobrachella eos (Ohler et al., 2011), Ophryophryne pachyproctus Kou, 1985, Xenophrys lancangica Lyu, Wang & Wang, 2023 (Megophryidae); Polypedates megacephalus Hallowell, 1861 (Rhacophoridae); Amolops tanfuilianae Sheridan et al., 2023, Hylarana cubitalis (Smith, 1917), Odorrana cf. chloronota (Günther, 1876) (Ranidae); Quasipaa ohlerae Pham et al., 2025, Limnonectes bannaensis Ye et al., 2007, Fejervarya limnocharis (Gravenhorst, 1829) (Dicroglossidae); and Tylototriton thaiorum Poyarkov, Nguyen & Arkhipov, 2021 (Salamandridae).
At the present moment, the habitat of the new species at the type locality in Xuan Lien NP is subject to serious anthropogenic pressure. The lowland areas at the foot of the ridges and hills are subjected to overgrazing of cattle, buffalos, and pigs, which destroy the leaf litter and upper layers of soil; pigs also dig up the banks of rivers and streams. In the forested areas along the rivers and in the valleys between the ridges, small mammals, birds, reptiles, and amphibians are actively hunted by locals. Forests at the foot of the mountain are being cut down for firewood; the open territories are subsequently adapted for planting vegetables, establishing paddy fields, or pastures. Xuan Lien NP was established as a protected area in part to protect the last remaining groves of Fokienia hodginsii (Dunn) A. Henry & H. H. Thomas (Cupressaceae), also known as “Po Mu,” a precious species of timber because of its characteristic aroma and its exceptional density. Despite the protective measures to preserve F. hodginsii and tourist routes to the most impressive old trees, during our surveys we have repeatedly observed illegal logging and transportation of F. hodginsii wood by the local population. These observations suggest that nature conservation measures in the Xuan Lien NP should be strengthened.
Given the information provided above, we suggest that I. griseivermis sp. nov., as with the majority of nominal caecilian species (
We presented the most complete matrilineal genealogy for the family Ichthyophiidae published to date, which includes four nominal species of the genus Uraeotyphlus and 26 nominal species of the genus Ichthyophis along with six potentially unnamed new species (see Fig.
Our study confirms that Southeast Asia represents an important secondary center of the Ichthyophiidae diversification, and all Southeast Asian Ichthyophis are monophyletic and belong to a single Clade D (see Fig.
As in the previous phylogenetic studies of the Ichthyophiidae (
Our matrilineal genealogy of Ichthyophiidae highlighted some issues within the genus Ichthyophis, which may result from either incomplete taxonomy of the group or from misidentification(s) (Fig.
Our study further highlights that the taxonomy of Southeast Asian Ichthyophis is far from complete. In addition to the three putative unnamed species of Ichthyophis from Malaysia reported by
In our matrilineal genealogy, three unstriped taxa of Ichthyophis from Indochina and southern China formed a group of closely related species: Ichthyophis griseivermis sp. nov., I. chaloensis, and I. yangi (Subclade D1, Fig.
There is a certain confusion regarding the origin of the I. laosensis holotype.
According to the terrestrial ecoregion classification of
Vietnam is recognized as a major center of amphibian diversity in Indochina (see
In recent years, the growing number of field surveys has led to advancements in herpetofauna assessments of northern and central Vietnam, which resulted in a rapid increase in the number of recorded amphibian species (e.g.,
The scarcity of materials in the herpetological collections and the small number of informative external morphological features represent some of the most important challenges in caecilian taxonomy. The novel non-destructive micro-CT scanning facilities represent an invaluable tool for obtaining data on osteology without damaging the precious specimens, which is especially important for caecilians where the samples are often minimal (
Despite the above-mentioned gaps in our knowledge, the current study revealed new features of the skull structure of the ichthyophiids, which suggest some evolutionary transformations in the cranial morphology of caecilians.
The traditional distinction between the two variants of skull structure—the stegokrotaphy (the complete coverage of the temporal region and the jaw-closing muscles by the squamosal) and the zygokrotaphy (the incomplete coverage of the temporal region by the squamosal and the presence of the temporal fossa)—appears to be simplified. It is possible to distinguish an intermediate variant of the structure of the temporal region, which we herein refer to as the weak zygokrotaphy, when the narrow temporal fossae are present, but only some of the posterior fibers of the long adductor of the jaw (m. adductor mandibulae longus) are attached to the dorsal surface of the cranial roof. At the same time, the weak zygokrotaphy found in ichthyophiids corresponds functionally (by the arrangement and involvement of the jaw-closing muscles) to that found in advanced caecilians (
It is widely recognized (
All of the Gymnophiona species have a peculiar dual jaw-closing mechanism, mainly provided by two muscles: the m. adductor mandibulae longus and the m. interhyoideus posterior (
In Rhinatrematidae, which represent the sister group of all other living caecilians, the m. adductor mandibulae longus is the main adductor muscle that is dominant over the m. interhyoideus posterior. It passes through the large temporal fossae and attaches to the sagittal crest of the parietals, leaving it exposed to external mechanical pressure (
Regardless of the condition of the temporal region (stegokrotaphy/weak zygokrotaphy), the ichthyophiids show a decrease in the importance of the m. adductor mandibulae longus in the jaw-closing mechanism compared to rhinatrematids. The m. interhyoideus posterior likely takes over the role of the main adductor of the jaw, and the retroarticular process of the pseudoangular is upturned to facilitate this function by moving the point of attachment to a more mechanically advantageous position. This condition of the retroarticular process is also present in all members of other, more advanced, groups of Gymnophiona. However, in the ichtyophiids, the upturned position of the retroarticular process is well pronounced, probably because the m. interhyoideus posterior itself is less developed and strong (
According to the classical works by
In addition, a comparison of the skulls of Ichthyophiidae and Rhinatrematidae revealed the presence of primitive cranial features in ichthyophiids that are found in stem caecilians (
We deposited the scans for this study into MorphoSource using the following links:
Ichthyophis griseivermis sp. nov. (
Ichthyophis nguyenorum (
Ichthyophis supachaii (
The fieldwork was completed within the framework and with partial financial support from the research project “Conservation, restoration, and sustainable use of tropical forest ecosystems based on the study of their structural and functional organization” of the Joint Vietnam-Russia Tropical Science and Technology Research Center for 2024 and 2025. Permission to conduct fieldwork in Vietnam was granted by the Department of Forestry, the Ministry of Agriculture and Rural Development of Vietnam, and local administrations of Thanh Hoa (permit numbers #562/GP of 01.06.2022 and #179/SNN&PTNT-CCKL of 05.10.2023) and Nghe An provinces (permit number #2089/UBND.VX of 03.04.2019). The authors are grateful to Andrey N. Kuznetsov (VRTC, Hanoi), Hoi Dang Nguyen (VRTC, Hanoi), Leonid P. Korzoun (MSU, Moscow), Bao Nguyen Le (DTU, Da Nang), as well as the Xuan Lien NP and Pu Hoat NR ranger staff (Thanh Hoa and Nghe An, Vietnam) for their support, organization, and assistance with the fieldwork. We express our gratitude to Alexey A. Polilov (MSU, Moscow) and Evgeny Scherbakov (MSU, Moscow) for the opportunity to conduct the osteological study using the Bruker Skyscan 1272 tomograph of the Biological Faculty of Moscow State University; and to the MSU HerpLab members, including Tang Van Duong (VNMN, Hanoi) for assistance in the lab. We thank Mark Wilkinson, Alexander Kupfer and an anonymous reviewer for numerous helpful comments and suggestions that allowed us to improve the previous version of the manuscript. This work was supported by the Russian Science Foundation to N. A. Poyarkov (Grant No. 22-14-00037-P, molecular and morphological analyses, data analysis) and in part by the Rufford Foundation to T. V. Nguyen (Grant No. 45888-2, data analysis).
List of abbreviations of the osteological features.
acot – anterior cotyles
addr – ridge for the attachment of the long adductor of the lower jaw
alp – anterolateral process of the quadrate
amp – anteromedial process of the quadrate
ap – anterior process
apm – anterior process of the maxillary part of the maxillopalatine
app – anterior process of the palatal part of the maxillopalatine
bas – os basale
bca – basicranial articulation
bp – basal process of the pterygoid
canp – canalis primordialis
car – foramen for the carotid artery
car+VIIpal – foramen for the carotid artery and the palatal branch of the facial nerve
ccr – curved crest of the vomer
ch – choana
con – condyle
corb – circumorbital
cp – columellar process of the stapes
cpr – ridge of the columellar process of the stapes
dm – depression for the attachment of the m. depressor mandibulae
dmp – dorsomedial process of the sphenethmoid
dp – dorsal process of the septomaxilla
dv? – incisure for the dorsal vein
fe – endolymphatic foramen
ff – facet for contact of the nasal or the parietal with the frontal
fm – foramen magnum
fp – facial process of the maxillopalatine
fper – perilymphatic foramen
fpl – footplate of the stapes
front – frontal
fv – fenestra vestibuli
jf – jugular foramen
lg – longitudinal grooves of the vomer
ltf – lower temporal fossa
mo – medial outgrowth of the atlantal centrum
mp – medial portion of the dorsal process of the septomaxilla
mpc - mediopalatinal cavity
mpp – medial process of the palatal part of the maxillopalatine
mxf – facet for contact of the prefrontal or the squamosal with the maxillopalatine
mxpal – maxillopalatine
na – neural arch
nas – nasal
nf – facet for contact of the frontal with the nasal
nos – nostrils
ns – nasal septum
oc – occipital condyles
orb+tc – orbit and tentacular canal
p – posterior process of the maxillary part of the maxillopalatine
pa – processus ascendens of the quadrate
par – parietal
pc – processus condyloideus of the pseudoarticular
pcot – posterior condyle
pdent – pars dentalis of the premaxilla
pdor – pars dorsalis of the premaxilla
pff – facet for contact of the frontal with the prefrontal
pi – processus internus of the pseudoarticular
pmp – premaxillary process of the vomer
pmx – premaxilla
po – processus oticus of the quadrate
pp – posterior process of the palatal part of the maxillopalatine
ppal – pars palatina
ppn – posterior process of the nasal
prfr – prefrontal
psa – pseudoangular
psaf – facet for contact of the pseudodentary with the pseudoangular
psd – pseudodentary
pt – pterygoid
pzp – postzygopophyseal processes
q – quadrate
qf – facet for contact of the squamosal with the quadrate
rap – retroarticular process
scs – spinal cord support
smx – septomaxilla
sn – sola nasi (processus conchoides)
snf – spinal nerve foramen
spf – facet for contact of the parietal with the sphenethmoid
sphen – sphenethmoid
sq – squamosal
sqf – facet for contact of the quadrate with the squamosal
st – stapes
sta – groove for the stapedial artery
tc – tentacular canal
utf – upper temporal fossa
vcf – foramina leading to the inner vomerine cavity
vf – vomeral foramen
vlr – ventrolateral ridge of the frontal
vo – vomeronasal organ cavity
vom – vomer
vp – ventral process of the septomaxilla
Id – incisure for the dorsal branch of the olfactory nerve
Iv – foramen for the ventral branch of the olfactory nerve
II – incisure for the optic nerve
Vim – foramen for the intermandibular branch of the trigeminal nerve
Vmd – foramina for the mandibular branch of the trigeminal nerve
Vmde – foramina for the external branch of the mandibular division of the trigeminal nerve
Vmx – foramen for the maxillary branch of the trigeminal nerve
Vop – foramen for the deep ophthalmic branch of the trigeminal nerve
Vop+mx – foramina for the maxillary and the deep ophthalmic branches of the trigeminal nerve
V+VII – groove for combined alveolar branches of the trigeminal and facial nerves
VII – foramen for the trunk of the facial nerve
VIIalv – foramina for the alveolar branch of the facial nerve
VIIIa – foramen for the anterior branch of the auditory nerve
VIIIp – foramen for the posterior branch of the auditory nerve
VIIImd – foramina for the medial branches of the auditory nerve
Figures S1–S17
Data type: .pdf
Explanation notes: Figure S1–S14. 3D reconstructions of the holotype of Ichthyophis griseivermis sp. nov. — Figure S15. The holotype of Ichthyophis griseivermis sp. nov. in preservative. — Figure S16. The paratype of Ichthyophis griseivermis sp. nov. in preservative. — Figure S17. Natural habitat of Ichthyophis griseivermis sp. nov. at the type locality in Xuan Lien NR, Thanh Hoa Province, northern Vietnam.
Tables S1–S5
Data type: .pdf
Explanation notes: Table S1. Primers used in this study. — Table S2. Sequences and voucher specimens of the family Ichthyophiidae and outgroup taxa used in this study. — Table S3. Characteristics of analyzed DNA sequences and the proposed optimal evolutionary models for gene and codon partitions as estimated in PartitionFinder 2.1.1. — Table S4. Uncorrected p–distances (percentage) between the sequences of cyt b mtDNA gene (above the diagonal) and between the sequences of 16S rRNA mtDNA gene (below the diagonal) of species of the genus Ichthyophis included in the phylogenetic analyses. — Table S5. Potentially diagnostic characters of skull morphology that differ among 15 species of the genus Ichthyophis.