Next Article in Journal
Pterosaur Tracks from the Upper Cretaceous Anacleto Formation (Neuquén Basin), Northern Patagonia, Argentina: Insights into Campanian Pterosaur Diversity in Gondwana
Previous Article in Journal
A New Genus and Two New Species of Fireflies from South America (Lampyridae: Lampyrinae: Photinini)
Previous Article in Special Issue
OTU Delimitation with Earthworm DNA Barcodes: A Comparison of Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phylogenetic Relationships in Earthworm Megascolex Species (Oligochaeta: Megascolecidae) with Addition of Two New Species

by
Azhar Rashid Lone
1,
Samrendra Singh Thakur
2,
Pooja Tiwari
1,
Samuel Wooster James
3 and
Shweta Yadav
1,*
1
Department of Zoology, School of Biological Sciences, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar 470003, India
2
Department of Biotechnology, School of Biological Sciences, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar 470003, India
3
Regenerative Agriculture Department, Maharishi International University, Fairfield, IA 52557, USA
*
Author to whom correspondence should be addressed.
Diversity 2022, 14(11), 1006; https://doi.org/10.3390/d14111006
Submission received: 26 October 2022 / Revised: 16 November 2022 / Accepted: 17 November 2022 / Published: 20 November 2022
(This article belongs to the Special Issue Molecular Phylogeny and Evolution of Earthworms)

Abstract

:
Megascolex (Oligochaeta: Megascolecidae) are endemic species to India and Sri Lanka, however, to date their molecular taxonomy and phylogenetic relationships have not been reported. We applied the first integrative approach using morpho-anatomical features and a COI dataset to unveil species delimitation (SD), molecular taxonomy, and phylogenetic relationships in Megascolex species. Our morpho-anatomical results revealed nine Megascolex species, namely, M. auriculata, M. cochinensis cochinensis, M. filiciseta, M. ratus, M. travancorensis travancorensis, M. triangularis, M. konkanensis konkanensis, M. polytheca polytheca, and M. polytheca zonatus. We also reported the occurrence of two new species, namely, M. papparensis sp. nov, and M. vazhichlensis sp. nov. Such findings were also supported by the analysed COI dataset, in which these new species appeared distinct on the phylogenetic trees with strong support. The studied Megascolex species appeared paraphyletic and formed three subclades on Bayesian inference (BI) and Maximum Likelihood (ML) phylogenetic trees. The first clade consisted of six species: M. cochinensis cochinensis, M. polytheca polytheca, M. polytheca zonatus, M. konkanensis konkanensis, M. filiciseta, and M. auriculata with strong posterior probability support. The second clade consisted of M. travancorensis travancorensis, M. papparensis sp. nov, and M. vazhichlensis sp. nov with strong support. The third clade consisted of M. ratus and M. triangularis with good support. In addition, the validation of species was confirmed by SD methods, in which the congruence among OTUs was observed with the clear barcode gap of 12–14% suggested by ABGD analysis. However, the species M. ratus and M. travancorensis travancorensis show deep intraspecific divergence and, therefore, require more sampling data. Such findings are essential to study the phylogenetics and evolution of the genus and, nonetheless, demand larger COI datasets to make concrete conclusions.

1. Introduction

Earthworms are the group of Oligochaeta worms that are known for their soil-dwelling properties and ecosystem functions that include soil turnover, soil fertility, and biomass [1,2]. Megascolecidae is the most diverse family of earthworms [3] comprising over 1000 species [3,4], most of them being native to Asia and Australia [4]. Out of many genera of the family Megascolecidae, the genus Megascolex was first described by Templeton [5] under the name M. caeruleus from the alpine regions of Ceylon (Sri Lanka), characterised by the presence of a row of small numerous spines or setae on each segment. Since then, 37 species and sub-species of Megascolex have been reported from Sri Lanka [6]. In India, the genus is comprised of 33 species [7], out of which 24 species were presented in a checklist by Narayanan et al. [8]. However, with the addition of M. lawsone from the state [9], the list was updated to 25 species with the presence of M. insignis as the only near-endemic species which occurs both in India and Sri Lanka [6]. The distribution of Megascolex species in India has been reported from peninsular parts, particularly in the southern part of India [9,10,11,12]. Kerala is an important constituent part of the Western Ghats and houses 99 species of earthworms [9]; however, just a few molecular studies have been conducted on Drawida and Moniligaster species [13] and there are no molecular records for Megascolex species.
Since earthworm taxonomy can be quite challenging particularly due to conservative morphological features [14], the taxonomy and diversity assessments are unstable and in need of revision. With the availability of online databases such as the Barcode of Life Database (BOLD) and the National Centre for Biotechnology Information (NCBI) that hold thousands of DNA barcodes, there are opportunities to use these information resources to improve the study of the taxonomy, evolution, and ecology of many organisms [15,16,17]. Nonetheless, in earthworms, these data have only been used to clarify a few taxonomic groups and the application of DNA barcoding to most of the Indian earthworm species is still deficient. This limits inferences about their actual diversity and correct classification [18]. Progress in molecular taxonomy of some Indian genera/species complexes such as Amynthas [19], Aporrectodea caliginosa species complex [20], Drawida and Moniligaster [13], Eutyphoeus [21,22], Kanchuria [23], as well as the molecular diversity and genetic variability of some other earthworm species [24,25,26,27] has significantly added to work in the field over the past few years. Nevertheless, most of the species remained untouched by the molecular approach and have only been described by morpho-anatomical observations that require strong diagnostic features, extensive labour, and expertise to segregate complex species and cryptic species.
The present investigation aimed to provide the DNA barcodes as well as the morpho-anatomical descriptions using an integrative approach in Megascolex species collected from different sites of the Western Ghats in Kerala to unveil their phylogenetic relationships on the tree of life.

2. Material and Methods

2.1. Study Area and Collection of Earthworms

Kerala is a small state located on the western tip of peninsular India (Figure 1), occupying 38,863 sq. km and comprising approximately 1.18% of the landmass of India [28]. The state is blessed with diversified micro- and macro-environmental conditions contributing to diverse flora and fauna. This mega biodiversity hot-spot area of the Western Ghats hosts about 30%, 24%, and 35% of animal, plant, and fresh water fish species of India [28] in spite of its small size. Earthworms were collected from different sites using the digging and hand-sorting techniques following the method described by Lewis and Taylor [29].

2.2. Morpho-Anatomical Observation

The preserved earthworm specimens were screened for morphological observations following dissections under a stereomicroscope (Leica Model No. M60). The specimens were stored in the museum, at the Department of Zoology, Dr. Harisingh Gour Vishwavidyalaya (A Central University) Sagar, MP. The species’ names with sample ID, process ID, collection sites, GPS coordinates, date of collection, name of collectors, and BOLD accession numbers are provided with the description of each species. The comparative account of morphological characters of Megascolex species is provided in Table 1.
Abbreviations used: sp., spermathecal pores; fp., female pore; mp., male pore; amp, ampulla; div., diverticulum; sp.d., spermathecal duct.

2.3. DNA Extraction, Amplification, and Sequencing

For sequencing, 28 samples of collected Megascolex species were sent to Barcode of Life Data Systems (BOLD), Biodiversity Institute of Ontario, University of Guelph, Canada to obtain DNA barcodes and accession numbers. The COI dataset of various earthworm species used in sequence alignments and data analysis for the current study is given in Supplementary Materials Data S1.

2.4. Sequence Data Curation, Alignment, and Genetic Distances

The alignment of the final COI dataset was performed in MEGA X [30] using MUSCLE [31]. Additionally, sequence composition was calculated and the sequences with insertions, deletions, and stop codons were removed, and all sequences were checked for ambiguous nucleotide sites and saturation. The uncorrected pairwise genetic distances were estimated within and between species [32] in MEGA X. In addition, the cluster sequencing analysis was performed on BOLD to infer intraspecific distances amongst Megascolex species compared to the distances of their nearest neighbours (NN).

2.5. Phylogenetic Analysis

Maximum likelihood (ML) and Bayesian inference (BI) methods were applied to reconstruct phylogenetic trees. The final dataset for nucleotide substitution model fit for phylogenetic tree reconstruction was tested by using jModelTest 2 [33] with the GTR+I+G4 model, at seed 419195. The maximum likelihood analysis was undertaken using raxmlGUI [34]. The bootstrap resampling with 1000 replicates was performed to support the individual branches of the ML tree. The final tree was visualised and edited on iTol v5 (Interactive Tree of Life (https://itol.embl.de; assessed on 4 October 2022) [35]. For BI analysis, the program BEAST [36] was employed and the parameters were set to 1,000,000 MCMC (Markov chain Monte Carlo) steps with the tree sampled at every 10,000 generations while discarding the first 10% as burn-in. Using Tracer v1.6 [37], the output was examined for all parameters and the maximum clade credibility tree was accessed using Tree Annotator v1.5.3 [38]. The final tree was visualised and edited using Fig Tree v1.4.2 [39]. The posterior probabilities with >0.90 were set for each branch in BI analysis. The topological differences between ML and BI phylogenetic trees were compared at the level of resolution obtained by each taxon and its bootstrap support.

2.6. Species Delimitation (SD) Analysis

Species delimitation was performed using three standardised methods; Automated Barcode Gap Discovery (ABGD) [40], Assemble Species by Automatic Partitioning (ASAP) [41], and the Refined Single Linkage (RESL) algorithm, which forms the basis of the Barcode Index Number system (BIN) [42]. ABGD was implemented with the MUSCLE aligned matrix as the input file in the ABGD package of iTaxoTools [43], adopting Kimura (k80) model = 2.0, X (relative gap width) = 1.0, and keeping other parameters as default values (Pmin = 0.001, Pmax = 0.1; Steps 10; Nb bins = 20). For ASAP, the output file was used in the ASAP program with default settings of iTaxoTools. BIN analysis was conducted in data analysis option of the BOLD database.

3. Results

3.1. Megascolex auriculata Aiyer, 1929

Megascolex auriculata Aiyer, 1929: 64.
Material examined: Sample ID: KERL0272A4; Process ID: IEW414-17; towards south of Periyar National Park, submerged area, Kottayam (9°27′22.6′′ N 77°13′53.8′′ E), Kerala, India; 28 October 2015; Coll. Samuel W James and Shweta Yadav; BOLD accession No: ADH2012.
Description:
Length, 180 mm; diameter, 2 mm; segments, 318. Prostomium epilobous. Clitellum in xiv–1/2 xvii, dark-pigmented. First dorsal pore at 9/10. Setae lumbricine loses orientation posterior. Setae ab closely and cd widely paired; the distance between ab same throughout the body, while cd varied after lxi segment and d moved outwardly. At the extreme lower end, setae 9 or 10 in number and setae e very close to d without maintaining specific orientation. On xviii, two opposite-oriented separate trenches present. Male pores bounded with pit-like depression, encircled by circular lobe-like depression. Female pores paired in xiv between aa encircled by whitish oval circle. Spermathecal pores two pairs in 7/8 and 8/9 in lines b (Figure 2).
Septa 6/7–9/10 very thick and 10/11–11/12 less thickened. Large barrel-shaped gizzard in vi. Oesophageal swelling at vii-xiii with numerous villi-like structures on its inner surface. Testes and funnels single pair in segment xi. Seminal vesicles large follicular in xii. Intestine begins in xvi. Last heart in xiii. The ectal half of the duct is wider than ental. Penial setae absent. The spermathecae are in viii and ix, ampullae long club-shaped with cylindrical diverticulum. Nephridia tufts are found, with up to six segments with one-two pair tufts, 8–9 tufts with long coiled tubes in the clitellar area, and 5–6 small-sized tufts in the post-clitellar region. Prostates racemose, loosely lobed in xviii-xxii segments, duct arises from the anterior part of the gland after 3–4 coils and opens outside. The duct emerges from the anterior portion of the gland after 3–4 coils and opens outside.
Distribution: The species is endemic to Kerala. Dist. Idukki, Kumili; Vandiperiyar Kottayam; Athirampuzha [8,10].

3.2. Megascolex cochinensis cochinensis Stephenson, 1915

Megascolex cochinensis Stephenson, 1915: 96–97.
Megascolex cochinensis cochinensis Stephenson. Blakemore 2007: 33.
Material examined: Sample ID: KERL0265A4; Process ID: IEW378-17; near waterfall, Adimali (10°02′03.8′′N 76°56′24.0′′E), Idukki, Kerala, India; 15 September 2015; Coll. Samuel W James and Shweta Yadav; BOLD Accession No: ADH2818. Sample ID: KEL17-773-27A25; Process ID: IEW727-17; Range Mannamangalam Peechi-Vazhani Wildlife Sanctuary (10°31′40.8′′N 76°22′26.0′′E), Thrissur, Kerala, India; 1 September 2017; Coll. Shweta Yadav; BOLD Accession No: ADL2102. Sample ID: KEL17-881-33A17; Process ID: IEW835-17; Vavlla sector, Chimmini Wildlife Sanctuary (10°26′55.3′′ N 76°27′32.8′′ E), Thrissur, Kerala, India; 1 September 2017; Coll. Shweta Yadav; BOLD Accession No: ADL1460. Sample ID: KEL17-882-33A18; Process ID: IEW836-17; Vavlla sector, Chimmini Wildlife Sanctuary (10°26′55.3′′ N 76°27′32.8′′ E), Thrissur, Kerala, India; 1 September 2017; Coll. Shweta Yadav; BOLD Accession No ADL1460. Sample ID: KEL17-886-33A22; Process ID: IEW840 17; Vavlla sector, Chimmini Wildlife Sanctuary (10°26′55.3′ N 76°27′32.8′′ E) Thrissur, Kerala, India; 1 September 2017; Coll. Shweta Yadav; BOLD Accession No ADL1460.
Description:
Large-sized worm; length, 190–208 mm; diameter 3–4 mm. Segments 125–180. Colour darkly pigmented in three specimens collected from same location and one non-pigmented collected from different location. Clitellum pale yellow in dark-coloured specimens extending over ½ xiii–1/3 xvii, setae present in clitellar region. Prostomium closed epilobic, first dorsal pore 5/6. Setae perichaetine, ring broken, intersetal intervals irregular, pre-clitellar aa = 2 or 3 ab, clitellar aa = ab, post-clitellar aa = 3 ab, setal counts were 46/v, 50/ix 54/xii, 42/xx, 36/xl. Male pores wavy slit-like on 18 segment, wavy slits approaching each other posteriorly within a light-coloured oval swollen area. Female pores single (in one specimen seems paired) mid-ventrally in line of setal ring in an oval whitish patch, on 14 segment extends laterally beyond aa. Spermathecal pores two pairs in segments 7/8 and 8/9, at a and sometimes on a/b or b lines (Figure 3).
Septa 4/5 and 5/6 are thin; 6/7–11/12 moderately thickened; 12/13–13/14 thickened. Gizzard large barrel-shaped in v. Calciferous glands absent. Intestine begins in xiv. Last hearts in xiii. Excretory system is micronephridal; large thick tufts from oesophagus to gizzard; 5–6 pairs of bushy tufts in clitellar region. Testes and male funnels free in x and xi. Seminal vesicles large in xi and xii attached to anterior septum. Prostates’ tightly packed lobules appear as a solid structure in xviii; duct thick, wide, and straight. Ovaries in xiii. Spermathecae two pairs in viii and ix, ampulla ovoid, duct long, ental part shiny. No penial setae.
Distribution: The species is near endemic and has been located in Thrissur: Forest Tramway (nr. Vazhachal); Dist. Kottayam, Athirampuzha [8].

3.3. Megascolex filiciseta Stephenson, 1915

Megascolex filiciseta Stephenson, 1915: 94–96.
Material examined: Sample ID: KERL0272A2; Process ID: IEW412-17; towards south of Periyar National Park, submerged area, Kottayam (9°27′22.6′′ N 77°13′53.76′′ E), Kerala, India; 28 October 2015; Coll. Samuel W James and Shweta Yadav; BOLD accession No: ADH2018.
Description:
Length, 86 mm; diameter, 3 mm; segments, 140. Colour dorsally dark-bluish grey; extremely dark at posterior ends; 3–4 fine dark stripes present dorsally. Thick and stiff body wall. Clitellum not very clear; dorsally slight difference in colour observed in xiii-xvi. Prostomium epibolic marked by median groove, and major part hidden in transverse groove. First dorsal pore 5/6. In pre-clitellar region aa = 2 ab, while distance between all setae was much reduced from clitellar region, it was difficult to observe the accurate distance between the setae; the longitudinal lines of setae broke dorsally while being ventrally distinct, numbers 28/v, 30/ix, and 22/xxv. Male pores in xviii on small porophores between a and b lines. Female pore was difficult to observe; appears to be in a whitish presetal ovoid patch present on xiv between aa. Spermathecal pores minute; two pairs in 7/8 and 8/9 at a (Figure 4).
Septa 5/6–7/8 delicate; 8/9–13/14 relatively thickened. Gizzard large barrel-shaped in vi. Calciferous glands absent. Intestine begins in xvi. Last hearts in xiii. Excretory system micronephridial; bushy tufts present behind iv with numerous tubules emerging from the main stem; loops prominent behind xiv, not attached to septs. Testes and male funnels are free in x and xi; seminal vesicles small; lobed in ix and xii. Prostates small, flattened, and confined to xviii; duct not visible. Spermatheca two pairs in viii and ix close to ventral nerve cord; ampulla ovoid; duct short; cylindrical diverticulum almost half length of ampullae. Penial setae present; shaft bow-shaped; tapers towards the distal ends; tip slightly curved with stout teeth.
Distribution: Only known from the type locality, i.e., Palakkad: Parambikulam, Kerala [8].

3.4. Megascolex konkanensis konkanensis Fedarb, 1898

Megascolex konkanensis Fedarb, 1898: 434–436. Michaelsen 1910: 75; Stephenson 1916: 328.
Megascolex konkanensis konkanensis Fedarb. Blakemore, 2007: 34.
Material examined: Sample ID: KERL0267A4; Process ID: IEW381-17; Grassland, Prambikulam Road, Muthalamada (10°23′34.4′′ N 76°46′32.2′′ E), Kerala, India; 17 September 2015; Coll. Samuel W James and Shweta Yadav; BOLD accession No: ACU5977. Sample ID: KERL0267A6; Process ID: IEW383-17; Grassland, Prambikulam Road, Muthalamada (10°23′34.4′′ N 76°46′32.2′′ E), Kerala, India; 17 September 2015; Coll. Samuel W James and Shweta Yadav; BOLD accession No: ACU5977. Sample ID: KERL0270A5; Process ID: IEW397-17; Thenmala Reservoir, Shendurney Wildlife Sanctuary (8°53′03.7′′ N 77°10′03.4′′ E), Kollam, Kerala, India; 26 October 2015; Coll. Shweta Yadav; BOLD accession No: ACS2283. Sample ID: KERL0270A7; Process ID: IEW399-17; Thenmala Reservoir, Shendurney Wildlife Sanctuary (8°53′03.7′′ N 77°10′03.4′′ E), Kollam, Kerala, India; 26 October 2015; Coll. Shweta Yadav; BOLD accession No: ACS2283. Sample ID: KERL0271A1; Process ID: IEW411-17; Mlappara, Periyar National Park (9°27′22.6′′ N 77°13′53.8′′ E), Kollam, Kerala, India; 27 October 2015; Coll. Shweta Yadav; BOLD accession No: ACU6070.
Description:
Large-sized worm; length, 170–316 mm; diameter, 4–5 mm; segments, 260–316. Colour grey with bluish irregular marks. Prostomium epilobic with two short longitudinal grooves on the dorsal surface of the first segment. First dorsal pore in 5/6. Dorsal setal gaps; before clitellum zz equal to two times larger than yz, while behind clitellum less than a gap zz = 1/2 yz. Ventrally in front of clitellum aa = 2 ab and behind aa = 3 ab/4 ab. Setae visible from ii. The distance interval between ab and bc is also irregular. Setal numbers were 36/vi, 32/ix, 38/xii, and 34/xxi. Clitellum 14–17. Male pores on xviii on transverse oval papillae placed towards lateral side of the segment (the actual pore unrecognizable); papillae enlarged dumbbell-shaped; surrounded by darker area of corresponding shape; occupy xvii–ix segments; anteriorly a little close towards the middle line. Female pores minute in oval white space between aa in xiv, while exact pore is not visible. Spermathecal pores tiny in 7/8 and 8/9 in line with setae e (Figure 5).
Septa 5/6–11/12 moderately thickened. Gizzard large in vi. Pharynx large with one pair mucous glands in v. Last heart in xiii. Intestine begins in xviii. Spermatheca two pairs in viii and ix; pear-shaped elongated ampullae with little small/equal-sized duct. In two specimens, the anterior part of duct as wide as ampullae. The club-shaped large diverticulum; length equal to or little larger than duct. Well-developed globular seminal vesicles in xi-xii approaching towards each other. Prostate large, mop-like, and bushy with numerous finger-shaped lobules. The duct was thick, glistening, and straight-inwardly directed, without muscular sac. Without calciferous glands and penial setae.
Distribution: The species is native peregrine to Kerala and has been reported from different regions of the country: Dist. Kozhikode: Tiruvallur (Thiruvallur), Calicut (Kozhikode); Dist. Palakkad: Chitoor (Chittur); Dist. Malappuram: Tirur, Kanjikode, Palghat (Palakkad); Dist. Ernakulam: Kalady; Dist. Thrissur: Kavalai; Dist. Kottayam: Athirampuzha, Kottayam; Dist. Kollam: Kulattupuzha (Kulathupuzha), Maddathoray (Madathara), Pathanapuram, Quilon (Kollam), Shasthancottah (Sasthamkotta); Dist. Alappuzha: Kerumaadi; Dist. Thiruvananthapuram: Trivandrum (Thiruvananthapuram); Travancore [8].

3.5. Megascolex polytheca polytheca Stephenson, 1915

Megascolex polytheca Stephenson, 1915: 89–90.
Megascolex polytheca polytheca Stephenson. Blakemore 2007: 36.
Material examined: Sample ID: KERL0276A5; Process ID: IEW441-17; Reserve Forest, Kanan Devan Hills (10°06′39.2′′ N 77°05′28′′ E), Munnar Kerala, India; 27 October 2015, Coll. Shweta Yadav; BOLD accession No: ADH2014. Sample ID: KERL0264A5; Process ID: IEW373-17; (10°6′56.16′′ N 77°5′13.56′′ E), Eravikulam National Park, Munnar, Kerala, India; 14 September 2015, Coll. Shweta Yadav, BOLD accession No: ADH2014.
Description:
Length, 105–140 mm; diameter, 3 mm; segments, 159–194. Colour grey, slightly darker at anterior end. The anterior part of the body before clitellum relatively thick and stout with sharply demarcated segments, while post-clitellar region smoother and cylindrical. Setae on segment iii-xii are arranged on raised segmental equators. Clitellum not profoundly distinguishable; in 13–17 with visible setal rings. Prostomium epilobous with hinder end of open tongue, and a transverse groove at the front end of the tongue present. First dorsal pore at 4/5. Setae ring closed dorsally zz = 1/2 yz and ventrally aa = 3 ab, behind the male aperture aa = 4 ab or 5 ab, ab > bc. Setae a and b are relatively larger than other setae. Setae a to e are in regular longitudinal lines. Setal numbers counted were 51/v, 52/ix, 53/xiii, 48/xx, and 46/xxxi. Male pores on large circular raised papillae, which are enclosed within biconcave depressions in xviii. The papillae are confined to xviii in one specimen, while in the other partly in xvii. A dark-coloured oval spot present at aa in xiv, while female aperture not distinguished. Spermathecal pores 4–7 in numbers 7/8 and 8/9. On separating the lips of the groove, a row of 4–7 white dots visible on each side. These points are surrounded by a dark area and begin internally between the lines b or c(Figure 6).
Septa 4/5 very delicate; 5/6 and 6/7 slightly thickened; 7/8–11/12 highly thickened. The large barrel-shaped gizzard in v. Calciferous glands absen;, oesophagus dilated deep yellow and vascular in xii. Intestine begins in xix. The last heart in xiii. Micronephridial excretory system. Spermathecae are small, circular, and 4–7 in number. Ampullae of spermathecae club-shaped organ with a long stalk with backwardly directed dilated end and more or less parallel to each other in a close set row. Male funnels in x and xi and seminal vesicles attached to anterior walls of xi and xii and made up of ovoid lobules. The ovaries are composed of finger-like lobes and break up at their free ends into strings of ova. Prostates large in xviii composed of small lobes and closely compacted together; duct short, muscular, stout, and widened near its termination.
Distribution: This species is endemic to Kerala. Thrissur: Kavalai in forest tramway [8].

3.6. Megascolex polytheca zonatus Stephenson, 1915

Megascolex polytheca var. zonatus Stephenson, 1915: 90–91.
Megascolex polytheca zonatus Stephenson. Blakemore 2007: 36.
Material examined: Sample ID: KERL0273A3; Process ID: IEW418-17; (10°6′56.88′′N 77°5′20.4′′E), Eravikulam National Park, Kannan Devan Hills, Kerala, India; 29 October 2015, Coll. Shweta Yadav; BOLD accession No: ADH2015. Sample ID: KERL0273A4; Process ID: IEW419-17; near Agraharam resort (10°6′56.88′′ N 77°5′20.4′′ E), Eravikulam National Park, Kannan Devan Hills, Kerala, India; 29 October 2015, Coll. Shweta Yadav; BOLD accession No: ADH2015.
Description:
Length, 112–118 mm; diameter, 2.5–2.75 mm; segments, 115–122. Colour grey. The anterior part of the body before clitellum relatively thick and stout with sharply demarcated segments, while post-clitellar region relatively constant in diameter. Setae perichaetine, clitellum browner than body in 13–17 with visible setal ring. Prostomium similar to M.p. polytheca epilobous with hinder end of the tongue open, and a transverse groove at the front end of the tongue present. First dorsal pore at 5/6. Setae ring closed dorsally zz = 1/2 yz and ventrally aa = 2 ab, behind the male aperture aa = 3 ab, ab > bc. Setae a and b are relatively larger than other setae. Setae a to e are in regular longitudinal lines. Setal numbers counted were 48/v, 50/ix, 52/xiii, 48/xx, and 48/xxxi. Male pores on large circular in lines with setae b; raised papillae which are enclosed within biconcave depressions in xviii segment. The male region quite similar to M. polytheca polytheca occupies xviii. A dark-coloured oval spot present at aa in xiv, while female aperture not distinguished. Spermathecal pores numerous in 7/8 and 8/9. On separating the lips of the groove, a row of 4–7 white dots visible on each side. These points are surrounded by a dark area and begin internally between the lines b or c (Figure 7).
Spermathecae small, five in each side of spermathecal groove. Each is a club-shaped organ with a long stalk with backwardly directed dilated end and more or less parallel to each other in a close set row. Male funnels in x and xi and seminal vesicles attached to anterior walls of xi and xii and made up of ovoid lobules. The ovaries are composed of finger-like lobes and break up at their free ends into strings of ova. Prostates large in xviii composed of small lobes and closely compacted together; duct short, muscular, stout, and widened near its termination.
The internal anatomy agrees closely with M. polytheca polytheca, however, the spermathecae are relatively smaller than M. polytheca polytheca, with four spermathecae on each side of the spermathecal groove. The ampulla of the spermathecae is distinguishable from the duct, is ovoid, dark purple in colour, and can be easily differentiated from M. p. polytheca; the diverticulum is club-shaped and simple; duct is cylindrical.
Distribution: The species is endemic to Kerala; Dist. Palakkad: Parambikula [8].

3.7. Megascolex ratus Cognetti de Martiis, 1911

Megascolex ratus Cognetti, 1911: 500–502. Michaelsen 1913: 87; Stephenson 1916: 327; Aiyer 1929: 68.
Material examined: Sample ID: KERL0268A4; Process ID: IEW387-17; Neyyar Dam and Wildlife Sanctuary (8°32′00.4′′N 77°08′53.5′′E), Thiruvanathapuram, Kerala, India; 26 October 2015; Coll. Shweta Yadav; BOLD accession No: ADH2820. Sample ID: KERL0270A1; Process ID: IEW395-17; Rosemala, Shendurney Wildlife Sanctuary (8°53′03.7′′ N 77°10′03.4′′ E), Kollam, Kerala, India; 24 October 2015; Coll. Shweta Yadav; BOLD accession No: ADH2016. Sample ID: KERL0270A9; Process ID: IEW401-17; Rosemala, Shendurney Wildlife Sanctuary (8°53′03.7′′ N 77°10′03.4′′ E), Kollam, Kerala, India; 26 October 2015; Coll. Shweta Yadav; BOLD accession No: ADH2016. Sample ID: KERL0274A1; Process ID: IEW422-17; Mannoorkara, close to Peppara Dam Wildlife Sanctuary, (8°38′35.2′′ N 77°10′50.5′′ E), Thiruvananthapuram, Kerala, India; 26 November 2015; Coll. Shweta Yadav; BOLD accession No: ADH2819. Sample ID: KERL0274A9; Process ID: IEW427-17; Mannoorkara, close to Peppara Dam Wildlife Sanctuary (8°38′35.2′′N 77°10′50.5′′E), Thiruvananthapuram, Kerala, India; 26 November 2015; Coll. Shweta Yadav; BOLD accession No: ADH2019. Sample ID: KERL0275A1; Process ID: IEW430-17; Mannoorkara, close to Neyyar Wild Life Sanctuary (8°33′00.2′′ N 77°14′33.0′′ E), Vazhichal, Kerala, India; 27 November 2015; Coll. Shweta Yadav; BOLD accession No: ADH2820.
Description:
Length, 140–350 mm; diameter, 7–8 mm; segments, 160–180. Colour dorsally violet-brown, ventrally grey. Prostomium short, broad, tanylobous, dorsally with longitudinal furrows that do not reach the posterior margin. Segments x-xiv biannulate. Setae more crowded in the ventral line than in the dorsal region, circles of setae are interrupted in the mid-ventral line of genital region where aa = 1/2 ab. There is no dorsal break, while in pre-clitellar aa = 2 ab and behind the clitellum aa = 3 ab. At x segment 108 setae and at xviii 130 setae. Clitellum prominent in one specimen out of seven collected specimens while setae present in entire clitellar region. Clitellum saddle-shaped occupying xiv- xviii and without intersegmental furrows. On xiv it is little extended towards xiii. Male pores on xviii, in the lines of setae h; between the male pores the setae absent. Male pores are placed on whitish tubercles supported by swollen papillae. Paired genital markings in 17, 19–22, sometimes also in 16 and 23. Female pores in xiv presetal on a, transversely extended. Spermathecal pores two pairs in intersegmental furrows 7/8 and 8/9 in the lines of the setae f (Figure 8).
Septa 7/8 and 8/9 thickened. Gizzard large in vi. Intestine begins in xiv. Last hearts in xiii. Two pair of testes in x and xi encapsulated in a large lobulated structure, which is compressed between strong septa. Prostates with strong, muscular cylindrical duct, which terminates into little folded lobular structure. Two pairs of Spermathecae in viii and ix with sac-shaped transversely striped ampullae, broader distally and rounded proximally. The duct opens into a tiny diverticulum, which is enclosed in the duct wall.
Distribution: The species is near endemic to Kerala and has been reported in Dist. Thiruvananthapuram: Bonaccord (Bonacaud, Bonakkad), Chimungi (Chemmunji), Coorloon, Mukkunni Reserve Forest, Trivandrum (Thiruvananthapuram) [8].

3.8. Megascolex travancorensis travancorensis Michaelsen, 1910

Megascolex travancorensis Michaelsen, 1910: 72–73.
Megascolex travancorensis travancorensis Michaelsen. Blakemore 2007: 37.
Material examined: Sample ID: KERL0275A2; Process ID: IEW431-17; close to Neyyar Wild Life Sanctuary (8°33′0.216′′ N 77°14′33′′ E), Vazhichal, Kerala, India; 27 October 2017; Coll. Shweta Yadav; BOLD accession No: ADH2010. Sample ID: KERL0269A3; Process ID: IEW389-17; Peppara Wildlife Sanctuary (8°37′14.7′′ N 77°10′3.36′′ E), Thiruvananthapuram, Kerala, India; 25 October 2017; Coll. Shweta Yadav; BOLD accession No: ADH2017.
Description:
Length, 118–122 mm; diameter, 3.0–3.2 mm; segments, 182–198. Colour grey. Prostomium epilobic. Setae ii-viii enlarged and closely paired anteriorly, while onwards loses orientation, more or less in longitudinal lines. Setal counts were 12/v, 12/xiii, 16/xxii, 20 xl. Clitellum not distinct. Dorsal pores start at 4/5. Male pores in xviii on raised egg-shaped cushions in lines b, within longitudinal curved slits in each cushion, run posterior to anterior from the base of circular rings. Female pore tiny in xiv, very hard to recognise. Spermathecal pores hard to recognise, apparently two pairs in 7/8 and 8/9 in lines b (Figure 9).
Septa 6/7–12/13 thickened. Gizzard in vi. Last hearts in xiii. Micronephridial system. Two pairs of testis funnels free in x and xi. Seminal vesicles small; densely packed racemose in xi and xii. Prostates fairly large, elongated, rectangular, and deeply incised with a cracked uneven surface occupying xvi-xxii; duct fairly long, shiny, thick distal part goes straight forward, and proximal part goes back to the irregular curve. Proximal part of the duct longer and covered by the glandular part. Spermathecae large, ampullae pear-shaped, dark in colour, distally narrowed, and severely bent at the opening. The duct is narrow than ampullae. A slender club-shaped distally somewhat bent diverticulum present. No penial setae.
Distribution: The species is endemic to Kerala and reported so far from Dist. Kollam: Kottarakkara, Kulathupuzha, Dist.Thiruvananthapuram: Killipalam, Pallode (Palode) [8].

3.9. Megascolex papparensis sp. nov.

LSIDurn:lsid:zoobank.org:act:BC2E4221-2984-4697-BC73-EEF09894E051
Holotype: Clitellate specimen (Sample ID: KERL0274A14; Process ID: IEW366-17); registration number: DHSGV-ZDM-H004; BOLD accession No: ADH2299, Mannoorkara, close to Peppara Dam and Wildlife Sanctuary (8°38′35.2′′ N 77°10′50.5′′ E), Thiruvananthapuram, Kerala, India; 26 November 2017, Coll. Shweta Yadav.
Paratypes: Clitellate specimen (Sample ID: KERL0274A14-1; Process ID: IEW366-17); registration number: DHSGV-ZDM-272015-013; BOLD accession number: ADH2299, Mannoorkara, close to Peppara Dam and Wildlife Sanctuary (8°38′35.2′′ N 77°10′50.5′′ E), Thiruvananthapuram, Kerala, India; 26 November 2017, Coll. Shweta Yadav.
Clitellate specimen (Sample ID: KERL0274A14-2; Process ID: IEW366-17); registration number: DHSGV-ZDM-272015-014; BOLD accession No: ADH2299; collection site similar to other paratypes.
Description:
Length, 85 mm; diameter, 3 mm; segments, 196. Prostomium epiobic with raised growth. Colour light grey. Setae paired in longitudinal lines aa = 2 ab = 3 ab in pre-clitellar region, aa = 1.5 ab in post-clitellar region, while setae absent in clitellar region. Setal counts were 12/v, 17/xii, 20/xx, and 21/xxx. Clitellum dark-brown coloured in 13–17. Dorsal pores start at 5/6. Male pores ca ½ circumference apart in a groove in xviii on oval glandular papillae, both papillae united by v-shaped groove at the distal end. The male region is slightly extended, sometimes appears in xix, while on deep observation it remains in xviii and touches the boundary of xix. Female pores paired in light-coloured oval patch in aa. Spermathecal pores two pairs in 7/8 and 8/9 in a lines (Figure 10).
Septa 6/7–12/13 thickened. Gizzard large in vi. Last hearts in xiii. Funnels free in xi and xii. Seminal vesicles compact in xi and xii. Prostate large, occupies xviii-xxiv segments; duct proximally very thin and curved. Spermathecae ampulla oval and slightly bent at anterior end; duct long and slightly bent; diverticulum shorter than ampulla.
Etymology: “papparensis” is derived from its type of habitation, Peppara Dam, Kerala.
Remarks: Particularly in terms of the prostatic duct, this species is similar to Megascolex travancorensis var. quilonensis [44]. However, the structure of spermathecae varies, especially in diverticulum length.
Variations: The main variation observed between holotype and paratype was the size of the prostate. In the holotype, the prostate covered six segments (xviii–xxiv) and in the paratypes it occupied eight segments, xvii-xxiv and xviii-xxv, respectively. Further, the colour of the clitellum in the holotype was darker as compared to paratypes, and in the holotype the clitellum occupied xiii-xvii segments and in paratypes it occupied xiii ½, xvi ½, or xvii segments.

3.10. Megascolex triangularis Stephenson, 1925

Megascolex triangularis Stephenson, 1925: 56–57.
Material examined: Sample ID: KERL0264A6; Process ID: IEW374-17; Reserve Forest (10°06′56.2′′ N 77°05′13.6′′ E), Kannan Devan Hills, Kerala, India; 29 October 2017; Coll. Shweta Yadav; BOLD accession No: ADH2301; Sample ID: KERL0273A1; Process ID: IEW416-17; Reserve Forest (10°06′56.9′′ N 77°05′20.4′′ E), Kannan Devan Hills, Kerala, India; 29 October 2017; Coll. Shweta Yadav; BOLD accession No: ADH2301.
Description:
Length, 138–150 mm; diameter, 5–6 mm; segments, 190–210. Dorsal pores start in 6/7. Clitellum 14–17. Prostomium? Colour dark brownish. Setae more than 50 per segment. Male pores paired, slit-like structure, situated at a in centre of the curved grooves, which are concave laterally and span 17/18 to 18/19. Female pores in 14th in whitish oval patch. Spermathecal pores two pairs in 7/8 and 8/9 at ab lines (Figure 11).
Septa 5/6–11/12 thickened. Gizzard large in v. Intestine begins in iv. No calciferous glands present. Male funnels in 10 and 11. Seminal vesicles in 11 and 13. Last heart in 13. Prostate racemose in 18; no prostatic duct seen. Spermathecae ampulla large ovate; duct short and wide without any demarcation; diverticula multiloculate. Penial setae absent.
Distribution: The species is known from the original description; Dist. Thrissur: Kavalai, Kerala [8].

3.11. Megascolex vazhichlensis sp. nov.

LSID. urn:lsid:zoobank.org:act:7C257761-1B02-4100-A738-9C39A2D5B857
Material examined: Holotype Clitellate specimen (Sample ID: KERL0275A3; Process ID: IEW449-17); registration number: DHSGV-ZDM-H005; BOLD accession No: ADH2821. Forest land (8°33′00.2′′ N 77°14′33.0′′ E), Vazhichal Kerala, India; 27 November 2017, Coll. Shweta Yadav.
Description:
Length, 175 mm; diameter, 3 mm; colour mild grey; segments, 256; secondary annulation in vii, viii, and ix. The anterior end is truncated (not tapering); prostomium is small and triangular, the pointed posterior angle directed upwards. First dorsal pore in 4/5. In pre-clitellar region aa > 2 ab; in clitellar region aa = 1/2 ab = 2 ab. The ventral setae are in definite longitudinal lines; setae of viii and ix are relatively small. The dorsal setae are in irregular intervals zz = 4–5 yz. The numbers were 16/v,18/xi, 21/xx, and 30/xxx. Setae at irregular intervals in post-clitellar region with smaller dorsal break in posterior region ca zz = 3 yz. Clitellum ½ xiii–xvi. On segment xviii, two projected U-shaped grooves connecting two male pore papillae. The two grooves are united horizontally at the base. It is difficult to observe exact position of male aperture, apparently in lines of b and c. The female pores (?) paired represented by two separate white spots between aa. The spermathecal apertures are small in 7/8 and 8/9 in lines of setae b, while slits laterally extended (Figure 12).
Septa 5/6 very thin; 6/7–10/11 relatively thickened; from 13/14 considerably thickened. Gizzard large barrel-shaped in v. Intestine begins in xvi. The last hearts in segment xiii. The micronephridia are present as large tufts from v to ix; thick covering of micro nephridia present in clitellar region and behind the prostates. Testes and funnels are free in x and xi. Seminal vesicles composed of lobules, attached to the anterior walls of segments xi and xii. Those in xi are small and those in xii are moderate in size. Prostates are flat, occupy xvi-ix; duct long, thick, curved, shining, and backwardly directed. Ovaries are not visible. Spermathecae are sausage-shaped, bent towards their free end, and dilated at the extremity. The duct is short, moderately stout, half as thick as the ampulla, single posteriorly curved elongated and club-shaped; shining diverticulum present.
Etymology: “vazhichlensis” is derived from its type of habitation, the Vazhichal Kerala.
Remarks: The species resembles Megascolex pepparensis sp. nov., particularly in the male region, where the prostate of M. vazhichlensis sp. nov occupies four segments (xvi–xix) with a consistently thickened duct. A large-sized prostate (xviii-xxiv) with a thin curved ental end are characteristic of Megascolex pepparensis sp. Further, spermathecae of Megascolex pepparensis sp. nov. have oval-shaped and slightly bent ampullae and M. vazhichlensis sp. nov. have club-shaped ampullae.

3.12. Intraspecific and Interspecific Genetic Distances

The mean nucleotide sequence composition was A = 28.6; T = 30.4; G = 18.6; and C = 22.4 with GC% = 40.04, GC% at codon 1 = 53.25, GC% at codon 2 = 42.26, and GC% at codon 3 = 24.80, which has been reported in many earthworm species [13,45]. The values of uncorrected pairwise genetic distances are provided in Figure 13. In Megascolex species, the mean intraspecific distance recorded was 1.232% and the minimum interspecific distance recorded was 15.18%. The minimum interspecific distance recorded was between M. polytheca polytheca and M. polytheca zonatus (14.3%) and between M. polytheca polytheca and M. cochinensis cochinensis (14.5%). Contrary to this, the maximum interspecific distance was found between M. vazhichlensis sp. nov and M. auriculta (20.3%) and between M. vazhichlensis sp. nov and M. triangularis (19.8%), respectively. Additionally, we recorded highest intraspecific genetic distance in M. travancorensis travancorensis (9.0%), followed by M. ratus (6.0%). The result of the cluster sequencing analysis indicated no overlapping interactions between genetic distances among Megascolex species (Table 2).

3.13. Phylogenetic Studies and Species Delimitation

The phylogenetic studies based on the COI dataset fell into distinct clusters in BI (Figure 14) and ML trees (Figure 15). Except for the species represented by single COI sequences, the phylogenetic analysis shows full support for each taxon as they were fully recovered on both BI and ML trees with strong supports. Moreover, there were minor differences in topologies of BI and ML trees. Furthermore, the three species delimitation methods produced different OTUs (Figure 16). The ABGD method found 15 OTUs, while ASAP analysis at 2.0 and 3.0 ASAP scores gave 18 and 13 OTUs, respectively. The BIN analysis also revealed 18 OTUs. Additionally, the ABGD depicted a clear barcode gap of 12–14% within Megascolex species (Figure 17).

4. Discussion

In the present investigation, the collected species from the Kerala part of the Western Ghats were discriminated into nine species, namely, Megascolex auriculata, M. cochinensis cochinensis, M. filiciseta, M. ratus, M. travancorensis travancorensis, M. triangularis, M. konkanensis konkanensis, M. polytheca polytheca, and M. polytheca zonatus, based on their morpho-anatomical observations. In addition, two new species were added to the genus, i.e., M. papparensis sp. nov. and M. vazhichlensis sp. nov. This segregation was further confirmed based on the DNA barcoding approach. In DNA barcoding, the reliability depends on a clear discontinuity between values of intraspecific and interspecific genetic divergences. In our report, the ABGD reflects the clear and no overlap barcode gap of 12–14%, which supports the validation and accuracy of DNA barcoding in the delimitation of Megascolex species. Moreover, in DNA barcoding, the species were considered distinct if their maximum intraspecific distances were less than the distances to their nearest neighbours (NN). In our cluster sequencing report, the maximum intraspecific genetic distance in each species was less than the distance to their nearest neighbour. Furthermore, based on such a limited COI dataset, certain conclusions were inferred. Firstly, looking at the family levels in both BI and ML phylogenetic trees, the species of the family Moniligastridae formed a distinct monophyletic clade in both phylogenetic trees (BI and ML), within a phylogenetic Megascolecidae (a monophyletic family of Crassiclitellata), to which the studied Megascolex species belongs. These results are not in accordance with the findings of James et al [46], which showed that Moniligastridae is the sister taxon to the Crassiclitellata. Secondly, the genus Megascolex appeared paraphyletic as revealed by BI and ML phylogenetic trees. On the BI tree, three clades were found with good posterior probability supports. Clade one with a high posterior probability of 0.98 was composed of M. polytheca polytheca, M. polytheca zonatus, M. konkanensis konkanensis, M. cochinensis cochinensis, M. filiciseta, and M. auriculata, respectively. Subsequently, the second clade with a good posterior probability value of 0.75 was formed by M. travancorensis travancorensis, M. ratus, and M. triangularis species. The last clade consisted of M. travancorensis travancorensis, M. papparensis sp. nov., and M. vazhichlensis sp. nov with a higher posterior probability value of 0.99. Although similar clades were observed on the ML tree, only the clade formed by M. travancorensis travancorensis, M. papparensis sp. nov., and M. vazhichlensis sp. nov show good clade support of 75, whereas the other two clades show weak supports. Additionally, in a phylogenetic tree study, a model of short internal and long external branches was detected in species, namely, M. ratus and M. travancorensis travancorensis, which is due to the high genetic divergence commonly found in some earthworm species. We followed three SD approaches to delimit species, and the results show incongruence in OTUs (BIN (18), ABGD (15), ASAP 2.0 (18), and ASAP 3.0 (13)). All the approaches show congruence in terms of the number of OTUs for species M. triangularis, M. polytheca polytheca, and M. polytheca zonatus, respectively. Contrary to this, M. travancorensis travancorensis was split by BIN, ASAP 2.0, and ASAP 3.0 into putative species, although it was merged by ABGD analysis. Similarly, in M. cochinensis cochinensis, M. konkanensis konkanensis, and M. ratus, the BIN, ASAP 2.0, and ASAP 3.0 split the individuals of these species into putative species. Generally, the incongruence in OTUs was seen in species of M.cochinensis cochinensis, M. konkanensis konkanensis, and M. ratus with high intraspecific genetic divergence. Moreover, the high intraspecific divergence in M. travancorensis travancorensis, M. cochinensis cochinensis, M. konkanensis konkanensis, and M. ratus was retained in each of them as single species due to their lesser intraspecific distance compared to their nearest neighbours (NN).
The high divergence in M. cochinensis cochinensis could be explained in that some of the individuals of this species were collected from Chimmini Wild Life Sanctuary, Kerala, and a few others were taken from Peechi-Vazhani Wildlife Sanctuary and Adimali Idukki Waterfall, respectively, with an average distance of 100 km and varied habitats. Such variations in habitats may perhaps induce selection [47], and, due to poor dispersal capability [48], they show large intraspecific divergence. A similar case is with M. konkanensis konkanensis and M. ratus, in which the individuals of these species differed in the habitats where they had been collected. Moreover, owing to a long evolutionary history [49], a direct effect of selection forces the soil-dwelling invertebrates to evolve to survive in specific habitats and stabilizing selection may perhaps select forms that diverge from the morphological optimum [47]. Conversely, the high genetic divergence found in earthworms can be explained by poor dispersal capability as they travel only limited distances/year [48], excluding the occurrence of passive dispersal by waterways or vertebrate predators [50]. This low dispersal ability is generally reflected as isolation by distance flow, where the genetic segregation is extremely linked to the geographical distance of the species.
The taxa which occur naturally, occupy limited geographical ranges, and are restricted to a specific geographical region are termed endemic species [51]. India is known to host diverse endemic fauna and flora because of its two major hotspots, namely the Himalayas and the Western Ghats. The Western Ghats mountain range located in the southwestern portion of India along with Sri Lanka is known for its biodiversity hotspots [52,53] and contains several endemic species. The level of endemism in these regions is quite common among various taxa such as fishes, land snails, trees, amphibians, odonates, and reptiles [54,55]. This sort of endemism is also true for earthworms that exceptionally show a high level of endemism (71.6% in Sri Lanka and 77% in the Western Ghats) as reported by Narayanan et al. [6,56]. The Megascolex species with their origins in Gondwanaland are ancient lineages with 36 and 25 species [9] in Sri Lanka and Kerala, a part of Western Ghats, respectively. It is noteworthy that the earthworm fauna of Sri Lanka possesses a close relation to the Western Ghats mountains of India, especially in the Kerala state [8]. Given the Indian mainland and the present Sri Lanka island were together before the breakup of the Gondwana land [57], comparative studies may highlight an understanding of the evolutionary history of Megascolex species, which could be achieved with the supplementations of their molecular data. Since the high rate of endemism is of immense bio-geographical significance nonetheless, the endemism of species is predominantly influenced by poor regional survey and taxonomic impediments and, therefore, the status of given species may alter with its distribution range and expansion. In such cases, the present study serves as a barcode reference library to infer the study of their molecular taxonomy and phylogenetic relationships. Moreover, the integrative methods are essential to work on the species that are considered difficult to discriminate exclusively on taxonomical features alone. Therefore, involving integrative methods not only provides species delimitations but also deciphers their molecular phylogeny, evolution, and population ecology. Such studies have already begun in Southeast Asian countries [58,59], America [60,61], and Europe [62,63], which has opened new avenues in the fields of their ecology, conservation, and sustainable development. Therefore, to achieve such objectives, there is an urgent demand for molecular data on Indian earthworms that are not only diverse in the country but also endemic to the sub-continent.

5. Conclusions

In biological research, species play a central role in all its branches and are the fundamental unit for measuring biodiversity, particularly soil invertebrates, which are of immense importance. Nonetheless, precise estimation of species and their boundaries is challenging. With the emergence of integrative methods to delimit species both at the classical taxonomical platform, using morpho-anatomical screening as well as DNA barcoding, the true level of biodiversity at few interfamilial/generic levels and species complexes has been estimated. The current report provided the first barcode library of endemic Megascolex species of India that have been collected from the Kerala part of Western Ghats. Moreover, we reported two new Megascolex species, i.e., papparensis sp. nov. and M. vazhichlensis sp. nov., along with the occurrence of nine species, namely, M. auriculata, M. cochinensis cochinensis, M. filiciseta, M. ratus, M. travancorensis travancorensis, M. triangularis, M. konkanensis konkanensis, M. polytheca polytheca, and M. polytheca zonatus from the Kerala state of the Western Ghats. Additionally, the phylogenetic relationship of Megascolex species fell into three distinct clusters and appeared paraphyletic on the trees with strong branch supports. We also reported a clear barcode gap of 12–14% in studied Megascolex species, reflecting the potential of DNA barcoding in these species. The study may aid in unveiling a better understanding of the phylogenetic relationships and evolution of Megascolex species that are important endemic resources of India.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/d14111006/s1, The COI dataset used in phylogenetic analysis is provided in Supplementary Materials Data S1.

Author Contributions

A.R.L. identified worms using integrated approach of taxonomy, performed phylogenetic analysis and drafted the manuscript; S.S.T. performed molecular analysis; P.T. drawn illustration; S.W.J. and S.Y. collected the samples, edited and approved the final draft of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

Authors are highly thankful to the Department of Biotechnology (Grant No.BT/PR7998/AGR/21/366/2013) Ministry of Science and Technology, Government of India, New Delhi for the financial support to carry out the study.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The molecular data used in the study are available on DS-IEWMG on Diversity studies in earthworms of India using DNA Barcode approach on Barcode of Life Data System (BOLD).

Acknowledgments

The authors are thankful to Principal Chief Conservator Forest (PCCF), Kerala, and Forest Officers of Periyar National Park, Peechi-Vazhani Wildlife Sanctuary, Chimmini Wildlife Sanctuary, Prambkulam Tiger Reserve, Shendurney Wildlife Sanctuary, Neyyar Wild Life Sanctuary and Peppara Wild Life Sanctuary, Kerala, India for granting permission to access reserve forests.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lavelle, P.; Spain, A.V. Soil Ecology; Kluwer Scientific Publications: Amsterdam, The Netherlands, 2001. [Google Scholar]
  2. Lavelle, P.; Decaëns, T.; Aubert, M.; Barot, S.; Blouin, M.; Bureau, F.; Margerie, P.; Mora, P.; Rossi, J.P. Soil invertebrates and ecosystem services. Eur. J. Soil Biol. 2006, 42, 3–15. [Google Scholar] [CrossRef]
  3. Blakemore, R.J. Cosmopolitan Earthworms—An Eco-Taxonomic Guide to the Peregrine Species of the World; CD-ROM: Canberra, Australia, 2002; p. 500. [Google Scholar]
  4. Sims, R.W.; Gerard, B.M. Earthworms. Keys and Notes to the Identification and Study of the Species; Synopsis of the British Fauna (New Series); Linnean Society of London Burlington House: London, UK, 1985; Volume 31. [Google Scholar]
  5. Templeton, R. Description of Megascolex caeruleus. Proc. Zool. Soc. Lond. 1844, 12, 89–91. [Google Scholar]
  6. Narayanan, S.P.; Kumari, S.; Kurien, V.T.; Thomas, A.P.; Paliwal, R.; Julka, J.M. A comprehensive checklist of the earthworms (Annelida: Clitellata: Megadrili) of Sri Lanka; a component of the Western Ghats–Sri Lanka biodiversity hotspot. Trav. Du Muséum Natl. D’histoire Nat. Grigore Antipa 2021, 64, 7–36. [Google Scholar] [CrossRef]
  7. Thakur, S.S.; Yadav, S. Exploration of Earthworms of India through Online Digital Library; Earthworms—The Ecological Engineers of Soil; IntechOpen: London, UK, 2018. [Google Scholar] [CrossRef] [Green Version]
  8. Narayanan, S.P.; Sathrumithra, S.; Christopher, G.; Thomas, A.P.; Julka, J.M. Checklist of the earthworms (Oligochaeta) of Kerala; a constituent of Western Ghats biodiversity hotspot, India. Zootaxa 2016, 4193, 117. [Google Scholar] [CrossRef]
  9. Narayanan, S.P.; Thomas, B.; Sreeraj, P.R.; Joseph, R.; Sathrumithra, S.; Kurien, V.T.; Anuja, R.; Moorkoth, S.; Kunnath, J.J.; Thomas, A.P.; et al. The first record of Megascolex lawsoni (Bourne; 1886) (Clitellata: Megascolecidae) from the state of Kerala, India. Megadrilogica 2019, 24, 68–72. [Google Scholar]
  10. Narayanan, S.P.; Sathrumithra, S.; Kuriakose, D.; Christopher, G.; Thomas, A.P.; Julka, J.M. Earthworms (Oligochaeta: Megadrile) from the Mahatma Gandhi University campus, Kottayam, Kerala. Malabar Trogon 2014, 12, 1–3. [Google Scholar]
  11. Kumar, H.T.S.; Sreepada, K.S.; Narayanan, S.P.; Reynolds, J.W. Megadrile earthworms (Annelida: Oligochaeta) around Udupi power corporation limited (upcl), Udupi district; Karnataka, south-west coast of India. Megadrilogica 2018, 23, 79–91. [Google Scholar]
  12. Sathrumithra, S.; Narayanan, S.P.; Anuja, R.; Kurian, P.; Thomas, A.P.; Julka, J.M.; Reynolds, J.W. Diversity of earthworms (Annelida: Oligochaeta) in Konni, a part of the Western Ghats of Kerala India. Megadrilogica 2018, 23, 57–68. [Google Scholar]
  13. Thakur, S.S.; Lone, A.R.; Tiwari, N.; Jain, S.K.; James, S.W.; Yadav, S. A contribution to the earthworm diversity (Clitellata; Moniligastridae) of Kerala, a component of the Western Ghats biodiversity hotspot, India, using integrated taxonomy. Anim. Biodivers. Conserv. 2021, 44, 117–137. [Google Scholar] [CrossRef]
  14. Pérez-Losada, M.; Breinholt, J.W.; Aira, M.; Domínguez, J. An updated multilocus phylogeny of the Lumbricidae (Annelida: Clitellata: Oligochaeta) earthworms. J. Phylogenetics Evol. Biol. 2015, 3, 140. [Google Scholar] [CrossRef] [Green Version]
  15. Hubert, N.; Hanner, R.; Holm, E.; Mandrak, N.E.; Taylor, E.; Burridge, M.; Watkinson, D.; Dumont, P.; Curry, A.; Bentzen, P.; et al. Identifying Canadian freshwater fishes through DNA barcodes. PLoS ONE 2008, 3, e2490. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Laopichienpong, N.; Muangmai, N.; Supikamolseni, A.; Twilprawat, P.; Chanhome, L.; Suntrarachun, S.; Peyachoknagul, S.; Srikulnath, K. Assessment of snake DNA barcodes based on mitochondrial COI and Cytb genes revealed multiple putative cryptic species in Thailand. Gene 2016, 594, 238–247. [Google Scholar] [CrossRef] [PubMed]
  17. Li, Q.; Kong, L.; Yu, H.; Zheng, X.; Yu, R.; Dai, L.; Sun, Y.; Chen, J.; Liu, J.; Ni, L.; et al. DNA barcoding reveal patterns of species diversity among north-western Pacific molluscs. Sci. Rep. 2016, 6, 1–17. [Google Scholar]
  18. Lalthanzara, H.; Lalfelpuii, R.; Zothansanga, C.; Vabeiryureilai, M.; Senthil, K.N.; Gurusubramanium, G. Oligochaete taxonomy—The rise of earthworm DNA barcode in India. Sci. Vis. 2018, 18, 1–10. [Google Scholar] [CrossRef] [Green Version]
  19. Vabeiryureilai, M.; Zothansanga, C.; Lalchhanhima, M.; Kumar, N.S.; Lalthanzara, H. Study on the Amynthas (Kinberg; 1867) earthworm (Megascolecidae: Oligochaeta) diversity through DNA barcoding from Northeast India. J. Environ. Biol. 2020, 41, 867–873. [Google Scholar] [CrossRef]
  20. Kumari, S.K.; Narayanan, S.P.; Julka, J.M. DNA barcoding of Aporrectodea caliginosa species-complex earthworms revealed the presence of Aporrectodea tuberculata (eisen; 1874) and Aporrectodea trapezoides (Dugés, 1828) from the Western Himalaya, India. Ann. Agri Bio Res. 2021, 26, 20–24. [Google Scholar]
  21. Thakur, S.S.; Lone, A.R.; Tiwari, N.; Yadav, S. Exploring new records of Eutyphoeus sp. (haplotaxida: Octochaetidae) from Garo hills; Meghalaya, North Eastern state of India with use of DNA barcodes. Mitochondrial DNA Part A 2020, 31, 265–272. [Google Scholar] [CrossRef]
  22. Tiwari, N.; Lone, A.R.; Thakur, S.S.; James, S.W.; Yadav, S. Three uncharted endemic earthworm species of the genus Eutyphoeus (Oligochaeta: Octochaetidae) from Mizoram, India. Zootaxa 2021, 5005, 41–61. [Google Scholar] [CrossRef]
  23. Lone, A.R.; Tiwari, N.; Thakur, S.S.; Pearlson, O.; Pavlíček, T.; Yadav, S. Exploration of four new Kanchuria sp. of earthworms (Oligochaeta: Megascolecidae) from the North Eastern Region of India using DNA bar-coding approach. J. Asia Pac. Biodivers. 2020, 13, 268–281. [Google Scholar] [CrossRef]
  24. Lalthanzara, H.; Zothansanga, C.; Lalchhanhima, M.; Kumar, N.S.; Ngukir, J.; Kimsing, A.; Vabeiryureilai, M. Diversity and New Records of Earthworms in Arunachal Pradesh, Northeast India. J. Environ. Biol. 2020, 41, 874–883. [Google Scholar] [CrossRef]
  25. Kumari, S.; Kumari, C.; Narayanan, S.P.; Dev, K.; Julka, J. Phylogenetic Relatedness of Drawida japonica (Michaelsen, 1892) and Drawida nepalensis Michaelsen, 1907 (Clitellata: Moniligastridae) from the Western Himalaya, India. Ann. Biol. 2021, 37, 36–41. [Google Scholar]
  26. Tiwari, N.; Lone, A.R.; Thakur, S.S.; Yadav, S. Interrogation of earthworm (Clitellata: Haplotaxida) taxonomy and the DNA sequence database. J. Asia Pac. Biodivers. 2021, 14, 40–52. [Google Scholar] [CrossRef]
  27. Lone, A.R.; Thakur, S.S.; Tiwari, N.; Sokefun, O.B.; Yadav, S. DNA barcoding and genetic variability of earthworms (Clitellata: Oligochaeta) with new records from Mizoram, India. Org. Divers. Evol. 2021, 4, 1–15. [Google Scholar] [CrossRef]
  28. Sreedharan, T.P. Biological Diversity of Kerala: A survey of Kalliasseri Panchayat; Kannur district. Kerala Research Programme on Local Level Development; Centre for Development Studies: Thiruvananthapuram, India, 2004. [Google Scholar]
  29. Lewis, T.; Taylor, L.R. Introduction to Experimental Ecology; Academic Press: Cambridge, MA, USA, 1968; p. 401. [Google Scholar]
  30. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. Mega X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  31. Edgar, R.C. Muscle: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [Green Version]
  32. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: Molecular evolutionary genetics analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [Green Version]
  33. Darriba, D.; Taboada, G.L.; Doallo, R.; Posada, D. jModelTest 2: More models, new heuristics and parallel computing. Nat. Methods 2012, 9, 772. [Google Scholar] [CrossRef]
  34. Edler, D.; Klein, J.; Antonelli, A.; Silvestro, D. raxmlGUI 2.0: A graphical interface and toolkit for phylogenetic analyses using RAxML. Methods Ecol. Evol. 2021, 12, 373–377. [Google Scholar] [CrossRef]
  35. Letunic, I.; Bork, P. Interactive Tree of Life (iTOL) v5: An online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021, 49, 293–296. [Google Scholar] [CrossRef]
  36. Suchard, M.A.; Lemey, P.; Baele, G.; Ayres, D.L.; Drummond, A.J.; Rambaut, A. Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol. 2018, 4, vey016. [Google Scholar] [CrossRef] [Green Version]
  37. Rambaut, A.; Drummond, A.J.; Xie, D.; Baele, G.; Suchard, M.A. Posterior summarisation in Bayesian phylogenetics using Tracer 1.7. Syst. Biol. 2018, 7, 901–904. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Rambaut, A.; Drummond, A.J. Tree Annotator v1. 5.3: MCMC Output Analysis. 2009. Available online: https://beast.community/treeannotator (accessed on 4 October 2022).
  39. Rambaut, A. Fig Tree v1.4.2: Tree Figure Drawing Tool. 2014. Available online: https://tree.bio.ed.ac.uk (accessed on 4 October 2022).
  40. Puillandre, N.; Lambert, A.; Brouillet, S.; Achaz, G.J.M.E. ABGD, Automatic Barcode Gap Discovery for primary species delimitation. Mol. Ecol. 2012, 21, 1864–1877. [Google Scholar] [CrossRef] [PubMed]
  41. Puillandre, N.; Sophie, B.; Guillaume, A. ASAP: Assemble species by automatic partitioning. Mol. Ecol. Resour. 2021, 21, 609–620. [Google Scholar] [CrossRef] [PubMed]
  42. Ratnasingham, S.; Hebert, P.D. A DNA-based registry for all animal species: The Barcode Index Number (BIN) system. PLoS ONE 2013, 8, e66213. [Google Scholar] [CrossRef] [Green Version]
  43. Vences, M.; Miralles, A.; Brouillet, S.; Ducasse, J.; Fedosov, A.; Kharchev, V.; Kostadinov, I.; Kumari, S.; Patmanidis, S.; Scherz, M.D.; et al. iTaxoTools 0.1: Kickstarting a specimen-based software toolkit for taxonomists. Megataxa 2021, 6, 77–92. [Google Scholar] [CrossRef]
  44. Michaelsen, W. The oligochaete fauna of the Indian-Ceylonian region. Repr. Verh. Nat. Sci. Assoc. Hambg. 1910, 19, 1–108. [Google Scholar]
  45. Huang, J.; Qin, X.; Zhen, J.S.; Gui, L.T.; Zi, Y.S. Identifying earthworms through DNA barcodes. Pedobiologia 2007, 51, 301–309. [Google Scholar] [CrossRef]
  46. James, S.W.; Davidson, S.K. Molecular phylogeny of earthworms (Annelida: Crassiclitellata) based on 28S, 18S and 16S gene sequences. Invertebr. Syst. 2012, 26, 213–229. [Google Scholar] [CrossRef]
  47. Bickford, D.; Lohman, D.J.; Sodhi, N.S.; Ng, P.K.; Meier, R.; Winker, K.; Ingram, K.K.; Das, I. Cryptic species as a window on diversity and conservation. Trends Ecol. Evol. 2007, 22, 148–155. [Google Scholar] [CrossRef]
  48. Novo, M.; Almodóvar, A.; Fernández, R.; Trigo, D.; Cosín, D.J. Cryptic speciation of hormogastrid earthworms revealed by mitochondrial and nuclear data. Mol. Phylogenetics Evol. 2010, 56, 507–512. [Google Scholar] [CrossRef]
  49. Chang, C.H.; James, S. A critique of earthworm molecular phylogenetics. Pedobiologia 2011, 54, S3–S9. [Google Scholar] [CrossRef]
  50. Cameron, E.K.; Bayne, E.M.; Clapperton, M.J. Human-facilitated invasion of exotic earthworms into northern boreal forests. Ecoscience 2007, 14, 482–490. [Google Scholar] [CrossRef]
  51. Gaston, K.J. How large is a species’ geographic range. Oikos 1991, 61, 434–438. [Google Scholar] [CrossRef]
  52. Myers, N.; Mittermeier, R.A.; Mittermeier, C.G.; da Fonseca, G.A.B.; Kent, J. Biodiversity hotspots for conservation priorities. Nature 2000, 403, 853–858. [Google Scholar] [CrossRef] [PubMed]
  53. Mittermeier, R.A.; Turner, W.R.; Larsen, F.W.; Broos, T.M.; Gscon, C. Global Biodiversity Conservation: The Critical Role of Hotspots; ZaChos, F.E., Habel, J.C., Eds.; Biodiversity Hotspots; Springer Publishers: London, UK, 2011; pp. 3–22. [Google Scholar]
  54. Gunawardene, N.R.; Daniels, A.E.D.; Gunatilleke, I.A.U.N.; Gunatilleke, C.V.S.; Karunakaran, P.V.; Nayak, K.G.; Prasad, S.; Puyravaud, P.; Ramesh, B.R.; Subramanian, K.A.; et al. A brief overview of the Western Ghats—Sri Lanka biodiversity hotpsot. Curr. Sci. 2007, 93, 1567–1572. [Google Scholar]
  55. Sudasinghe, H.; Pethiyagoda, R.; Meegaskumbura, M. A molecular phylogeny of the genus Laubuka (Teleostei: Cyprinidae) in Sri Lanka reveals multiple origins and a cryptic species. Syst. Biodivers. 2020, 18, 592–613. [Google Scholar] [CrossRef]
  56. Narayanan, S.P.; Paliwal, R.; Kumari, S.; Ahmed, S.; Thomas, A.P.; Julka, J.M. Annelida: Oligochaeta. In Faunal Diversity of Biogeographic Zones of India: Western Ghats; Chandra, K., Raghunathan, C., Suresh, P.M., Subramanian, K.A., Rizvi, A.N., Eds.; Zoological Survey of India: Kolkata, India, 2020; pp. 87–102. [Google Scholar]
  57. Biswas, S. Did biotic impoverishment facilitate phenomenal diversification in Sri Lanka? Curr. Sci. 2008, 95, 1021–1025. [Google Scholar]
  58. Seesamut, T.; Jirapatrasilp, P.; Chanabun, R.; Oba, Y.; Panha, S. Size variation and geographical distribution of the luminous earthworm Pontodrilus litoralis (Grube, 1855) (Clitellata, Megascolecidae) in Southeast Asia and Japan. ZooKeys 2019, 862, 23. [Google Scholar] [CrossRef]
  59. Seesamut, T.; Jirapatrasilp, P.; Sutcharit, C.; Tongkerd, P.; Panha, S. Mitochondrial genetic population structure and variation of the littoral earthworm Pontodrilus longissimus Seesamut and Panha, 2018 along the coast of Thailand. Eur. J. Soil Biol. 2019, 93, 103091. [Google Scholar] [CrossRef]
  60. Somers, C.M.; Neudorf, K.; Jones, K.L.; Lance, S.L. Novel microsatellite loci for the compost earthworm Eisenia fetida: A genetic comparison of three North American vermiculture stocks. Pedobiologia 2011, 54, 111–117. [Google Scholar] [CrossRef]
  61. Ortíz-Gamino, D.; Cunha, L.; Martínez-Romero, E.; Flores-Estévez, N.; Ortíz-Ceballos, Á.I. Genetic diversity and population structure in an invasive pantropical earthworm along an altitudinal gradient. Peer J. Prepr. 2017, 5, e3024v1. [Google Scholar]
  62. King, R.A.; Tibble, A.L.; Symondson, W.O. Opening a can of worms: Unprecedented sympatric cryptic diversity within British lumbricid earthworms. Mol. Ecol. 2008, 17, 4684–4698. [Google Scholar] [CrossRef] [PubMed]
  63. Marchán, D.F.; Csuzdi, C.; Decaëns, T.; Szederjesi, T.; Pizl, V.; Domínguez, J. The disjunct distribution of relict earthworm genera clarifies the early historical biogeography of the Lumbricidae (Crassiclitellata, Annelida). J. Zool. Syst. Evol. Res. 2021, 59, 1703–1717. [Google Scholar] [CrossRef]
Figure 1. Collection of Megascolex species from various geographical spots of Kerala, Western Ghats, India.
Figure 1. Collection of Megascolex species from various geographical spots of Kerala, Western Ghats, India.
Diversity 14 01006 g001
Figure 2. Camera lucida illustration of Megascolex auriculata (a) genital region; (b) spermathecal region; (c) spermathecae.
Figure 2. Camera lucida illustration of Megascolex auriculata (a) genital region; (b) spermathecal region; (c) spermathecae.
Diversity 14 01006 g002
Figure 3. Camera lucida illustration of Megascolex cochinensis cochinensis (a) genital region; (b) spermathecal region; (c) spermathecae.
Figure 3. Camera lucida illustration of Megascolex cochinensis cochinensis (a) genital region; (b) spermathecal region; (c) spermathecae.
Diversity 14 01006 g003
Figure 4. Camera lucida illustration of Megascolex filiciseta (a) genital region; (b) spermathecal pore; (c) spermathecae.
Figure 4. Camera lucida illustration of Megascolex filiciseta (a) genital region; (b) spermathecal pore; (c) spermathecae.
Diversity 14 01006 g004
Figure 5. Camera lucida illustration of Megascolex konkanensis konkanensis (a) genital region; (b) spermathecal pores; (c) spermathecae.
Figure 5. Camera lucida illustration of Megascolex konkanensis konkanensis (a) genital region; (b) spermathecal pores; (c) spermathecae.
Diversity 14 01006 g005
Figure 6. Camera lucida illustration of Megascolex polytheca polytheca (a) genital region; (b) spermathecal pore; (c) spermathecae.
Figure 6. Camera lucida illustration of Megascolex polytheca polytheca (a) genital region; (b) spermathecal pore; (c) spermathecae.
Diversity 14 01006 g006
Figure 7. Camera lucida illustration of Megascolex polytheca zonatus (a) genital region; (b) spermathecal pore; (c) spermathecae.
Figure 7. Camera lucida illustration of Megascolex polytheca zonatus (a) genital region; (b) spermathecal pore; (c) spermathecae.
Diversity 14 01006 g007
Figure 8. Camera lucida illustration of Megascolex ratus (a) genital region; (b) spermathecal pore; (c) spermathecae.
Figure 8. Camera lucida illustration of Megascolex ratus (a) genital region; (b) spermathecal pore; (c) spermathecae.
Diversity 14 01006 g008
Figure 9. Camera lucida illustration of Megascolex travancorensis travancorensis (a) genital region; (b) spermathecae.
Figure 9. Camera lucida illustration of Megascolex travancorensis travancorensis (a) genital region; (b) spermathecae.
Diversity 14 01006 g009
Figure 10. Camera lucida illustration of Megascolex papparensis sp. nov. (a) genital region; (b) spermathecal region; (c) spermathecae.
Figure 10. Camera lucida illustration of Megascolex papparensis sp. nov. (a) genital region; (b) spermathecal region; (c) spermathecae.
Diversity 14 01006 g010
Figure 11. Camera lucida illustration of Megascolex triangularis (a) genital region; (b) spermathecal region; (c) spermathecae.
Figure 11. Camera lucida illustration of Megascolex triangularis (a) genital region; (b) spermathecal region; (c) spermathecae.
Diversity 14 01006 g011
Figure 12. Camera lucida illustration of Megascolex vazhichlensis sp. nov. (a) genital region (b) spermathecal region; (c) spermathecae.
Figure 12. Camera lucida illustration of Megascolex vazhichlensis sp. nov. (a) genital region (b) spermathecal region; (c) spermathecae.
Diversity 14 01006 g012
Figure 13. Intraspecific (bold diagonally) and interspecific uncorrected pairwise genetic distances of Megascolex species. Species with single COI sequence are represented by N/A.
Figure 13. Intraspecific (bold diagonally) and interspecific uncorrected pairwise genetic distances of Megascolex species. Species with single COI sequence are represented by N/A.
Diversity 14 01006 g013
Figure 14. BI phylogenetic tree, depicting paraphyletic clades of Megascolex species with posterior probabilities node values at each node. The families Moniligastridae and Megascolecidae are represented in blue and pink colours; grey colour shows outgroups.
Figure 14. BI phylogenetic tree, depicting paraphyletic clades of Megascolex species with posterior probabilities node values at each node. The families Moniligastridae and Megascolecidae are represented in blue and pink colours; grey colour shows outgroups.
Diversity 14 01006 g014
Figure 15. Maximum Likelihood (ML) phylogenetic tree with highest log likelihood (−4616.68), containing COI barcodes of earthworm species.
Figure 15. Maximum Likelihood (ML) phylogenetic tree with highest log likelihood (−4616.68), containing COI barcodes of earthworm species.
Diversity 14 01006 g015
Figure 16. Maximum likelihood dendrogram showing species delimitation analyses of the COI dataset of Megascolex species of Kerala (black rectangles). MORPHO = Morphology, BIN = Barcode Gap Index Number, ABGD = Automatic Barcode Gap Discovery, and ASAP = Assemble Species by Automatic Partitioning at 2.0 and 3.0 ASAP scores.
Figure 16. Maximum likelihood dendrogram showing species delimitation analyses of the COI dataset of Megascolex species of Kerala (black rectangles). MORPHO = Morphology, BIN = Barcode Gap Index Number, ABGD = Automatic Barcode Gap Discovery, and ASAP = Assemble Species by Automatic Partitioning at 2.0 and 3.0 ASAP scores.
Diversity 14 01006 g016
Figure 17. ABGD analysis showing a clear barcode gap in Megascolex species (X axis = genetic distances Y axis = frequency).
Figure 17. ABGD analysis showing a clear barcode gap in Megascolex species (X axis = genetic distances Y axis = frequency).
Diversity 14 01006 g017
Table 1. Comparative account of morphological characters of selected Megascolex species which belong to different clades in phylogenetic trees.
Table 1. Comparative account of morphological characters of selected Megascolex species which belong to different clades in phylogenetic trees.
Species
/Characters
M. papparensis sp. nov.M. vazhichlensis sp. nov.M. cochinensis cochinensisM. ratusM. triangularis
Size mm85175190–208140–350138–150
PigmentationLight greyMild greyDarkly colouredDorsally violet brown, ventrally greyDark brownish
Diameter mm333–47–85–6
No. of segments196256125–180160–180190–210
ClitellumDark brown, xiii–xvii½ xiii–xvi½ xiii–1/3 xviiSaddle shaped, xiv–xviiixiv–xvii
SetaeAbsent in clitellar regionIrrugular post-clitellur with small dorsal breakSetae ring brokenMore crowded ventrally, no break ?
Dorsal pore start5/64/55/6?6/7
Spermathecal poresa linesb linesa lines, sometimes on a/b or bf linesab lines
Male pores½ xviii xviii, in lines b and cxviiixviii, on hAt a
Female poresPaired, aaPaired, aaSingle, aapaired on aPaired, aa
Gizzardvivvviv
Testes and funnels Free, xi and xiiFree,
x and xi
Free,
x and xi
Free,
x and xi
Free,
x and xi
Seminal vesiclesxi and xiixi and xiix and xii?xi and xiii
Calciferous glands??absent?absent
SpermathecaAmpula oval-shaped, short diverticulum Suasage-shaped without shining diverticulumAmpula ovoid with shining diverticulumLarge sac-shaped ampula, tiny diverticulumAmpula large, multiloculate diverticulum
ProstatesLarge in xviii–xxivxvi–ixTightly packed in xviiiIn xviiixviii, no prostatic ducts seen
Penial setae??absent??
Table 2. Cluster Sequencing of COI barcodes of Megascolex species of Kerala depicting distances to nearest neighbour (NN).
Table 2. Cluster Sequencing of COI barcodes of Megascolex species of Kerala depicting distances to nearest neighbour (NN).
OTUsProcess IDAverage DistanceMax DistanceTaxonCountDistance to NN
OTU-1IEW366-1700M. papparensis sp. nov114.9206
OTU-2IEW373-17
IEW441-17
00M.polytheca polytheca2
2
13.7614
OTU-3IEW374-17
IEW416-17
00M.triangularis217.4603
OTU-4IEW378-1700M. cochinensis cochinensis13.2051
OTU-5IEW381-17
IEW383-17
0.31440.3144M.konkanensis konkanensis22.5396
OTU-6IEW387-17
IEW430-17
0.80120.8012M.ratus22.5641
OTU-7IEW389-1700M. travancorensis travancorensis19.0909
OTU-8IEW395-17
IEW401-17
0.32050.3205M.ratus26.0897
OTU-9IEW397-17
IEW399-17
00M.konkanensis konkanensis22.2222
OTU-10IEW411-1700M.konkanensis konkanensis12.2222
OTU-11IEW412-1700M.filiciseta115.0793
OUT-12IEW414-1700M. auriculata115.5555
OTU-13
IEW418-17
IEW419-17
00M. polytheca zonatus214.3540
OTU-14IEW422-1700M.ratus16.0897
OTU-15IEW427-1700M.ratus12.5641
OTU-16IEW449-1700M. vazhichlensis sp. nov.114.9371
OTU-17IEW431-1700M.travancorensis travancorensis19.0909
OTU-18
IEW835-17
IEW836-17
IEW840-17
00M. cochinensis cochinensis33.2051
OTU-19IEW727-1700M. cochinensis cochinensis14.6474
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lone, A.R.; Thakur, S.S.; Tiwari, P.; James, S.W.; Yadav, S. Phylogenetic Relationships in Earthworm Megascolex Species (Oligochaeta: Megascolecidae) with Addition of Two New Species. Diversity 2022, 14, 1006. https://doi.org/10.3390/d14111006

AMA Style

Lone AR, Thakur SS, Tiwari P, James SW, Yadav S. Phylogenetic Relationships in Earthworm Megascolex Species (Oligochaeta: Megascolecidae) with Addition of Two New Species. Diversity. 2022; 14(11):1006. https://doi.org/10.3390/d14111006

Chicago/Turabian Style

Lone, Azhar Rashid, Samrendra Singh Thakur, Pooja Tiwari, Samuel Wooster James, and Shweta Yadav. 2022. "Phylogenetic Relationships in Earthworm Megascolex Species (Oligochaeta: Megascolecidae) with Addition of Two New Species" Diversity 14, no. 11: 1006. https://doi.org/10.3390/d14111006

APA Style

Lone, A. R., Thakur, S. S., Tiwari, P., James, S. W., & Yadav, S. (2022). Phylogenetic Relationships in Earthworm Megascolex Species (Oligochaeta: Megascolecidae) with Addition of Two New Species. Diversity, 14(11), 1006. https://doi.org/10.3390/d14111006

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop