Multi-locus phylogenetic analyses uncover species boundaries and reveal the occurrence of two new entomopathogenic nematode species, Heterorhabditis ruandica n. sp. and Heterorhabditis zacatecana n. sp.
Ribosomal RNA (rRNA) gene sequences such as ITS sequences and the sequences of the D2–D3 expansion segments of the 28S rRNA are traditionally used for identification purposes and for novel taxonomic status descriptions of the species of the genus Heterorhabditis (Adams et al., 1998; Campos‐Herrera et al., 2011; Li et al., 2012; Malan et al., 2008; Nguyen et al., 2008; Rana et al., 2020; Spiridonov and Subbotin, 2016). As a recently evolved group, marginal variations in the rRNA gene sequences are expected in this genus, which limits the use of these genetic markers for taxonomic purposes, especially of closely related species (Blaxter et al., 1998; Blouin, 2002; Haag et al., 2018). In addition, the use of sequences containing several ambiguous nucleotides, potentially arisen from sequencing errors and/or poor quality-control, leads to erroneous taxonomic affiliations, as it is exemplified by the relatively high number of synonym species in the genus Heterorhabditis (Dhakal et al., 2020; Hunt and Nguyen, 2016). The use of mitochondrial DNA such as COI sequences, the gold standard taxonomic marker for species delimitation in the Kingdom Animalia, may help to overcome the taxonomic limitations of rRNA gene sequences. However, this taxonomic marker has been used only sporadically for identification purposes, barely used for taxonomy, and never used to describe new Heterorhabditis species (Chaubey et al., 2016; Hebert et al., 2003; Joyce et al., 1994a; Kuwata et al., 2007). As a consequence, the availability of COI sequences for this genus remained very limited for several years, limiting our understanding of the phylogenetic relationships of this genus (Chaubey et al., 2016; Dhakal et al., 2020; Kuwata et al., 2007).
To improve our understanding on the phylogenetic relationships of Heterorhabditis nematodes, to determine the most suitable genetic markers for the rapid and reliable identification of the species of this genus, specially of closely related species, and to determine species boundaries in this genus, we generated nucleotide sequences of several phylogenetically relevant gene markers and used them to reconstruct phylogenetic trees, to calculate sequence similarity scores, and to determine signatures of species- and population-specific genetic polymorphism. To improve our understanding on the taxonomic relationships of Heterorhabditis nematodes, we revisited the current literature related to the description, synonymisation, and declaration as species inquirendae of Heterorhabditis species to compile taxonomically relevant morphological and morphometric characters, characterized new nematode isolates at the morphological and morphometrical level, and conducted self-crossing and cross-hybridization experiments. Our study illustrates the importance of multi-locus sequencing for the characterization of new species within the genus Heterorhabditis, serves to clarify the phylogenetic relationships of these important biological control agents, and can inform future species descriptions to advance our efforts towards developing more tools for sustainable and environmentally friendly agriculture.
Self-crossing and cross-hybridization experiments were carried out on lipid agar plates as described by Dix et al. (1992). Heterorhabditis ruandica n. sp. Rw14_N-C4a and H. zacatecana n. sp. MEX-39 were self-crossed, hybridized with each other and with H. bacteriophora CH21 (Rana et al., 2020). For this, one second–generation male and one second–generation virgin female were placed on lipid agar plates (35 mm diam.) and incubated at 27°C. Ten independent plates per crossing type were set. Progeny production was observed daily for a period of five consecutive days. Experiments were repeated three times under the same conditions.
Nematode molecular characterization and phylogenetic relationships
Genomic DNA from about 10 to 20 thousand nematodes was extracted using the genomic DNA isolation kit following manufacturer’s instructions (Norgen Biotek Corp., Thorold, Ontario, Canada). The following genes/genomic regions were amplified by polymerase chain reaction (PCR): the D2–D3 expansion segments of the 28S rRNA, the internal transcribed spacer (ITS) region of the rRNA, the cytochrome c oxidase I (COI), the thin filament (F-actin)-associated protein (unc-87), and the calmodulin 1 (cmd-1). Primers used were selected based on previous publications (Dhakal et al., 2020; Joyce et al., 1994b; Regeai et al., 2009; Subbotin et al., 2006) (Table S2). PCR reactions consisted of 1 µL of genomic DNA, 12.5 µL of EmeraldAmp GT PCR Master Mix (Takara Bio, Shiga, Japan), 0.5 µL of both forward and reverse primers at 10 mM and 10.5 µL of dH2O. The PCR reaction was performed using a thermocycler (Mastercycler nexus gradient, Eppendorf, Germany) with the following settings: (i) for ITS and D2–D3, 1 cycle of 1 min at 98°C followed by 35 cycles of 10 sec at 98°C, 30 sec at 50°C, 1 min 30 sec at 72°C, and by a single final elongation step at 72°C for 10 min; (ii) for cmd-1 and unc-87, 1 cycle of 1 min at 98°C followed by 40 cycles of 10 sec at 98°C, 30 sec at 50°C, 30 sec at 72°C, and by a single final elongation step at 72°C for 10 min; (iii) for COI, 1 cycle of 1 min at 98°C followed by 40 cycles of 10 sec at 98°C, 30 sec at 40°C, 30 sec at 72°C, and by a single final elongation step at 72°C for 10 min. PCR was followed by electrophoresis (45 min, 100 V) of 5 μl of PCR products in a 1% TBA (Tris–boric acid–EDTA) buffered agarose gel stained with SYBR Safe DNA Gel Stain (Invitrogen, Carlsbad, California, USA). PCR products were purified using the FastGene Gel/PCR extraction kit (Nippon Genetics Co., Japan) and sequenced using reverse and forward primers by Sanger sequencing (Microsynth AG, Balgach, Switzerland). Obtained sequences were manually curated and trimmed and deposited in the NCBI under the accession numbers given in Table S3. Sequences of the following nematode strains were obtained in this study: Heterorhabditis ruandica n. sp. (isolates Rw18_M-Hr1a, Rw18_M-Hr1b, and Rw14_N-C4a), H. zacatecana n. sp. (isolates MEX-39, MEX-40, and MEX-41), H. bacteriophora (isolates DE2, DE6, PT1, IT6, EN01, and TT01); H. georgiana Hbb, H. beicherriana H06, H. indica CH7, and H. atacamensis MEX-20. To complete this data set and to obtain genomic sequences of nematodes that belong to all the validly described species of the genus Heterorhabditis, we searched the database of the National Center for Biotechnology Information (NCBI) by the Basic Local Alignment Search Tool (BLAST) using the accession numbers of the sequences obtained previously (Dhakal et al., 2020) (Table S3). Resulting sequences were used to reconstruct phylogenetic relationships by the Maximum Likelihood method based on the following nucleotide substitution models: Hasegawa–Kishino–Yano (HKY + I) (cmd-1), Tamura–Nei (TN93 + G + I) (COI), Kimura 2-parameter (K2 + G) (D2–D3 and ITS), and Tamura 3-parameter (T92) (unc-87). To select the best substitution model, best-fit nucleotide substitution model analyses were carried out in MEGA 7 (Hasegawa et al., 1985; Kimura, 1980; Kumar et al., 2016; Nei and Kumar, 2000). Sequences were aligned with MUSCLE (v3.8.31) (Edgar, 2004). The trees with the highest log likelihood are shown. The percentage of trees in which the associated taxa clustered together is shown next to the branches. Initial tree(s) for the heuristic search were obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using the Maximum Composite Likelihood (MCL) approach, and then selecting the topology with superior log likelihood value. In some cases, a discrete Gamma distribution (+G) was used to model evolutionary rate differences among sites and the rate variation model allowed for some sites to be evolutionarily (+I). The trees are drawn to scale, with branch lengths measured in the number of substitutions per site. Graphical representation and edition of the phylogenetic tree were performed with Interactive Tree of Life (v3.5.1) (Chevenet et al., 2006; Letunic and Bork, 2016).
Symbiotic relationships
The Photorhabdus entomopathogenic bacteria associated with H. ruandica n. sp. Rw14_N-C4a and H. zacatecana n. sp. MEX-39 nematodes were isolated as described by Machado et al. (2019), (2021b). Briefly, Galleria mellonella larvae (Lepidoptera: Pyralidae) were exposed to 100 nematode infective juveniles. Three to four days later, insect cadavers were surface–sterilized and cut open with a blade. Bacteria-digested internal organs were spread onto LB agar plates and incubated at 28°C for 24 to 48 h. Photorhabdus-like colonies were sub-cultured until monocultures were obtained. A single primary form colony was then selected and used for further experiments. Bacteria primary forms were determined by examining colony characteristics and by examining pigments uptake on NBTA plates (LB agar plates supplemented with 25 mg l–1 bromothymol blue and 4 mg l–1 triphenyl-2,3,5-tetrazolium chloride). The strains were further sub-cultured and maintained on LB agar plates at 28°C. To establish their taxonomic identities, we reconstructed phylogenetic relationships based on whole genome sequences of the isolated bacteria and all the different species/subspecies of the genus Photorhabdus (Machado et al., 2021a, b). To obtain genomic sequences, genomic DNA was extracted and purified using the GenElute Bacterial Genomic DNA Kit (Sigma-Aldrich, Switzerland) following manufacturer’s instructions. The resulting DNA was used for library preparation using the TruSeq DNA PCR-Free LT Library Prep (FC-121-3003) kit. Indexed libraries were then pooled at equimolar concentrations and sequenced (2 × 150 bp) on an Illumina HiSeq 3000 instrument. Genomes were assembled using the Bactopia pipeline (Petit and Read, 2020). Briefly, the raw Illumina reads were quality trimmed using Trimmomatic 0.39 (Bolger et al., 2014). The resulting reads were assembled with SPAdes 3.14.1 (k-mer sizes of 31, 51, 71, 91, and 111 bp) (Bankevich et al., 2012). Scaffolds with a mean read-depth smaller than 20% of the median read-depth of the longer scaffolds (≥5,000 bp) as well as scaffolds that were shorter than 200 bp were removed. The final assemblies were polished using Pilon 1.22 (Walker et al., 2014). Genome sequences were deposited in the National Centre for Biotechnology Information. Accession numbers are listed in Table S4. Phylogenetic relationships were reconstructed based on the assembled genomes and the genome sequences of all validly published species of the genus (Machado et al., 2021a, b). For this, core genome alignments were created using Roary 3.6.2 (Page et al., 2015). Using this alignment, a maximum likelihood tree was constructed using Fasttree 2.1.10 based on the Jukes-Cantor + CAT nucleotide evolution model (Price et al., 2009, 2010). These analyses were carried out in Galaxy (Afgan et al., 2018). Whole genome sequence similarities were calculated by the digital DNA-DNA hybridization (dDDH) method using the recommended formula 2 of the genome-to-genome distance calculator (GGDC) web service of the Deutsche Sammlung von Mikroorganismen und Zellkulturen GmbH (DSMZ) (Auch et al., 2010a, 2010b; Meier-Kolthoff et al., 2013, 2014).
Morphometrics of infective juveniles and adult generations of Heterorhabditis ruandica n. sp.
Male
Characters
Holotype
Paratypes
Hermaphrodite (1st Gen) paratypes
Female (2nd Gen) paratypes
Infective juvenile paratypes
n
1
20
20
20
25
Body length (L)
760
769 ± 60 (652-863)
3295 ± 286 (2907-4123)
1366 ± 123 (1131-1608)
544 ± 29 (496-591)
a (L/BD)
20.3
17 ± 1.5 (15-21)
14.1 ± 1.1 (11.7-16.1)
18 ± 1.4 (15-20)
24 ± 3.0 (20-29)
b (L/NL)
7.8
8.1 ± 1.0 (5.8-9.7)
23 ± 1.8 (21-27)
11.4 ± 1.2 (9.0-13.6)
4.7 ± 0.4 (4.1-5.4)
c (L/T)
26.2
31 ± 3.6 (23-36)
42 ± 5.7 (34-51)
20 ± 2.2 (16-24)
8.2 ± 1.0 (7.6-8.6)
c’ (T/ABW)
1.1
1.4 ± 0.2 (0.6-1.7)
2.2 ± 0.3 (1.7-2.6)
2.8 ± 0.5 (1.9-3.6)
4.6 ± 0.8 (3.4-5.8)
V (VA/L × 100)
–
–
48 ± 2.5 (45-55)
48 ± 2.3 (41-51)
–
Max. Body Width (MBD)
37.5
44 ± 3.0 (40-51)
233 ± 17 (209-274)
77 ± 4.0 (68-83)
23 ± 2.7 (18-27)
Lip region width
6.5
7.2 ± 0.8 (5.7-8.4)
12.4 ± 0.8 (11.0-14.0)
10.3 ± 0.9 (8.8-12.2)
–
Stoma length
9.5
11.1 ± 1.6 (8.7-13.9)
14.9 ± 1.4 (13-18)
13.6 ± 1.8 (10.4-16.0)
13.8 ± 1.2 (12.1-16.0)
Bulb length (BL)
18.5
20 ± 1.8 (18-25)
35 ± 3.6 (29-42)
27 ± 2 (23-30)
13.8 ± 1.8 (11.0-19.0)
Pharynx length (PL)
95.2
84 ± 7.1 (74-107)
128 ± 6.3 (118-142)
107 ± 6.9 (91-120)
102 ± 7.0 (91-115)
Nerve ring – ant. end (NR)
68
63 ± 5.2 (56-74)
93 ± 7.5 (78-108)
81 ± 6.4 (69-97)
55 ± 3.6 (52-64)
Excretory pore– ant. end (EP)
84.3
81 ± 10.1 (61-109)
121 ± 11 (106-153)
111 ± 10.8 (92-129)
78 ± 3.4 (70-89)
Neck length (Stoma+Pharynx, NL)
98
96 ± 7.3 (84-117)
143 ± 6.3 (134-159)
120 ± 6.0 (107-132)
115 ± 7.3 (103-131)
Body width at neck base
36
34 ± 1.9 (30-37)
119 ± 8.9 (101-138)
58 ± 4.3 (50-66)
18 ± 3.0 (15-24)
Vagina length
–
–
28 ± 4.0 (20-38)
19.2 ± 2.9 (15-26)
–
Body width at vulva
–
–
240 ± 21 (199-278)
78 ± 3.8 (72-85)
–
Vulva – ant. end (VA)
–
–
1581 ± 151 (1369-1882)
655 ± 47 (572-706)
–
Vulva – post. End (PV)
–
–
1713 ± 178 (1453-2241)
710 ± 89 (559-949)
–
Rectum length
–
–
36 ± 4.6 (29-49)
30 ± 3.8 (24-35)
8.5 ± 1.9 (6.1-13.7)
Anal body diam. (ABD)
26.1
18 ± 2.4 (15-25)
37 ± 5.5 (29-51)
25 ± 4.5 (18-34)
12.4 ± 1.8 (9.2-16.0)
Tail with sheath length (T)
–
–
–
–
56 ± 4.9 (49-64)
Tail without sheath length
29
25 ± 3.2 (21-29)
80 ± 7.9 (63-98)
68 ± 6.5 (62-88)
30.4 ± 4.5 (22-39)
Spicule length (SL)
49
43 ± 4.1 (34-50)
–
–
–
Gubernaculum length (GL)
20.2
18 ± 1.5 (15-21)
–
–
–
Stoma length/lip region width
1.5
–
1.2 ± 0.2 (1.0-1.6)
1.3 ± 0.1 (1.1-1.6)
–
Nerve ring % (NR/NL × 100)
69.4
67 ± 4.4 (58-75)
65 ± 5.0 (56-78)
67 ± 3.9 (61-75)
–
Excretory pore % (EP/NL × 100)
86
85 ± 8.5 (61-97)
85 ± 8.3 (67-103)
92 ± 7.6 (74-104)
–
Rectum % (R/ABD × 100)
–
–
90 ± 10 (80-130)
128 ± 29 (90-181)
–
D % (EP /PL × 100)
88.5
96 ± 9.2 (69-111)
95 ± 9.3 (74-114)
104 ± 9.3 (82-118)
78 ± 7.6 (66-98)
E % (EP/T × 100)
290
325 ± 49 (232-413)
153 ± 24 (120-205)
164 ± 23 (111-203)
139 ± 13.4 (112-168)
SW % (SL/ABD × 100)
242
243 ± 47 (150-306)
–
–
–
GS % (GL/SL × 100)
41.2
42 ± 5.2 (35-57)
–
–
–
H % (H/T × 100
–
–
–
–
46 ± 4.0 (37–55)
Morphometrics of infective juveniles and adult generations of Heterorhabditis zacatecana n. sp.
Male
Characters
Holotype
Paratypes
Hermaphrodite (1st Gen) paratypes
Female (2nd Gen) paratypes
Infective juvenile paratypes
n
1
20
22
22
25
Body length (L)
808.1
861 ± 29 (811-914)
5127 ± 494 (4408-6179)
2244 ± 203 (1954-2798)
539 ± 21 (493-578)
a (L/BD)
19.0
18 ± 1.6 (15-22)
16 ± 2.0 (13-20)
12.3 ± 1.2 (10.5-15.0)
22 ± 1.2 (19-24)
b (L/NL)
8.1
9.1 ± 1.1 (7.6-12)
26 ± 4.3 (20-34)
18 ± 1.8 (16-21)
5.0 ± 0.4 (4.4-5.9)
c (L/T)
28.9
34 ± 4.2 (26-43)
70 ± 10.4 (52-90)
39 ± 7.4 (31-63)
9.4 ± 0.6 (8.2-10.5)
c’ (T/ABW)
1.4
1.6 ± 0.3 (1.2-2.5)
1.6 ± 0.3 (1.2-2.4)
1.7 ± 0.2 (1.3-2.0)
5.3 ± 0.6 (4.3-6.7)
V (VA/L × 100)
–
–
48 ± 4.3 (36-57)
53 ± 4.2 (43-61)
–
Max. Body Width (MBD)
42.5
48 ± 3.6 (41-56)
319 ± 41 (235-385)
183 ± 23 (160-228)
24 ± 0.9 (23-27)
Lip region width
6.2
7.4 ± 0.7 (6.2-8.8)
11.7 ± 2.4 (9.2-19.2)
10.1 ± 1.0 (7.7-11.4)
4.0 ± 0.5 (3.2-5.2)
Stoma length
10
9.3 ± 1.0 (6.3-11)
19 ± 2.0 (14-23)
11.5 ± 1.7 (8.0-15.2)
13.5 ± 1.0 (12.0-15.3)
Bulb length (BL)
20.2
22 ± 2.4 (19-28)
40 ± 4.6 (28-49)
30 ± 2.6 (28-38)
20 ± 1.4 (17.1-23.0)
Pharynx length (PL)
95.2
86 ± 9.8 (57-100)
182 ± 23 (155-211)
113 ± 9.5 (101-133)
95 ± 7.2 (82-111)
Nerve ring – ant. end (NR)
65.4
66 ± 5.3 (60-78)
131 ± 22 (96-169)
83 ± 7.3 (71-96)
81 ± 6.3 (69-72)
Excretory pore– ant. end (EP)
96.2
93 ± 9.6 (77-109)
150 ± 24 (108-190)
113 ± 11 (100-133)
89 ± 6.8 (72-99)
Neck length (Stoma+Pharynx, NL)
99.3
96 ± 9.6 (71-108)
201 ± 21 (174-231)
124 ± 10 (112-148)
109 ± 6.9 (96-124)
Body width at neck base
34.5
36 ± 2.3 (31-40)
167 ± 13 (133-188)
95 ± 13.9 (74-121)
23 ± 1.3 (19-26)
Vagina length
–
–
31 ± 4.0 (24-36)
25 ± 6.4 (17-42)
–
Body width at vulva
–
–
331 ± 33 (257-379)
185 ± 27 (153-230)
–
Vulva – ant. end (VA)
–
–
2470 ± 279 (1959-3038)
1182 ± 129 (910-1397)
–
Vulva – post. end (PV)
–
–
2657 ± 279 (1990-3938)
1062 ± 147 (860-1455)
–
Rectum length
–
–
36 ± 4.6 (30-41)
27 ± 4.1 (19-39)
–
Anal body diam. (ABD)
19.6
17 ± 2.3 (13-22)
47 ± 8.1 (34-58)
35 ± 3.2 (31-41)
11.1 ± 1.3 (8.6-14.1)
Tail with sheath length (T)
–
–
–
–
58 ± 3.1 (52-63)
Tail without sheath length
28
26 ± 3.3 (21-33)
74 ± 8.3 (63-87)
58 ± 8.2 (45-75)
29.4 ± 2.5 (25-34)
Spicule length (SL)
54.1
45 ± 3.7 (38-55)
–
–
–
Gubernaculum length (GL)
18.7
20 ± 2.1 (15-25)
–
–
–
Stoma length/lip region width
1.6
–
1.6 ± 0.3 (1.1-2.1)
1.2 ± 0.2 (0.8-1.7)
–
Nerve ring % (NR/NL × 100)
65.9
69 ± 9.9 (61-96)
65 ± 9.4 (49-86)
67 ± 4.8 (60-82)
–
Excretory pore % (EP/NL × 100)
96.9
98 ± 17 (78-134)
75 ± 11 (51-95)
67 ± 4.8 (60-82)
–
Rectum % (R/ABD × 100)
–
–
79 ± 17 (54-112)
76 ± 13 (52-106)
–
D % (EP /NL × 100)
101.05
109 ± 21 (83-156)
75 ± 11 (55-95)
92 ± 7.9 (80-111)
94 ± 12 (68-120)
E % (EP/T × 100)
343.6
365 ± 68 (236-503)
206 ± 46 (145-303)
197 ± 27 (145-246)
154 ± 14 (128-184)
SW % (SL/ABD × 100)
276
270 ± 50 (170-320)
–
–
–
GS % (GL/SL × 100)
34.56
40 ± 10 (40-60)
–
–
–
H % (H/T × 100
–
–
–
–
47 ± 5.6 (35-56)
Comparative morphometrics of adult males of Heterorhabditis ruandica n. sp., H. zacatecana n. sp., and of different closely related Heterorhabditis species. All measurements are in µm (except ratios and percentages)
Comparative morphometrics of hermaphrodite females of Heterorhabditis ruandica n. sp., H. zacatecana n. sp., and of different closely related Heterorhabditis species.
Comparative morphometrics of adult females of Heterorhabditis ruandica n. sp., H. zacatecana n. sp., and of different closely related Heterorhabditis species.
Comparative morphometrics of infective juveniles of Heterorhabditis ruandica n. sp., H. zacatecana n. sp., and of different closely related Heterorhabditis species.
Body 0.65 to 0.86 mm long, C-shaped after fixation. Cuticle almost smooth, with transversal striae poorly developed. Lateral field not visible. Lip region with six lips developed but not fused bearing six acute labial papillae at oral margin and four rounded cephalic papillae at the base of lips. Oral opening almost rounded with thick margins. Amphidial apertures pore-like, ovoid and located posterior to lateral labial papillae. Stoma rhabditoid type, 1.2 to 2.3 times the lip region width, with short cheilostom with poorly refringent rounded cheilorhabdia, short gymnostom with refringent bar-like rhabdia, and long stegostom surrounded by the pharyngeal collar and bearing bar-like pro-mesorhabdia and small poorly refringent meta-telorhabdia. Pharynx poorly developed with robust corpus without differentiated metacorpus, short and slightly narrow isthmus and pyriform bulb with poorly visible valvular apparatus. Nerve ring encircling the isthmus at 58 to 75% of neck length, just anterior to basal bulb. Excretory pore located at basal bulb level, located at 61 to 97% of neck length. Cardia poorly developed, surrounded by intestinal tissue. Intestine without differentiations. Cardiac anterior end with thin walls. Genital system monorchic, laterally reflexed. Spicules well-developed, separate, with small angular manubrium, calamus poorly developed, and robust lamina with acute tip, scarcely prominent dorsal hump and poorly developed ventral velum. Gubernaculum with manubrium straight and slightly ventrally curved corpus, 40 to 50% of spicule length. Tail conoid with acute tip, ventrally curved posteriorly, flanked by the bursa. Bursa peloderan, with nine pairs of genital papillae (1 + 2/3 + 3), one of them probably the phasmid: three pairs pre-cloacal (GP1–GP3) and six pairs post-cloacal being three pairs at mid tail length (GP4–GP6) and three pairs (GP7–GP9) terminal; GP1 and GP2 more spaced, GP2 and GP3 closely spaced (Figs. 1–4).
Hermaphroditic females
Body 2.91 to 4.12 mm long, arcuate with general morphology similar to male, having labial papillae very acute and prominent. Nerve ring encircling the isthmus at 56 to 78% of neck length. Excretory pore located at or posterior to basal bulb, located at 67 to 103% of neck length. Genital reproductive system didelphic-amphidelphic with ovaries well developed, reflexed, oviducts and uteri not well visible, vagina very short and vulva small having transverse slit opening. Rectum slender, 0.8 to 1.3 times longer than the anal body diameter. Anus with prominent lips. Tail conoid with acute tip lacking mucro, having cellular part simple at its junction with the hyaline part. Phasmids inconspicuous (Figs. 1–4).
Amphimictic females
Body similar to, but usually smaller than hermaphroditic females, 1.13–1.61 mm long. Rectum very long, almost twice longer than the anal body diameter. Anus with posterior lip very prominent. Tail conoid with acute tip lacking mucro, having cellular part bifurcated at its junction with the hyaline part (Figs. 1–4).
Infective sheathed juveniles (J3 stage envolved by the J2 stage cuticle)
Body 0.5 to 0.6 mm long, with habitus slightly ventral curved after fixation. Cuticle with transversal striae at anterior end, with both transversal and longitudinal striae at neck region and only with longitudinal striae at rest of body. Lip region lacking differentiate lips, bearing six labial papillae and cephalic papillae not visible. Amphidial apertures very reduced. Oral opening closed, having triradial symmetry. Stoma tubular, about twice the lip region wide. Pharynx slender, with long and narrow corpus, very narrow isthmus and pyriform basal bulb. Nerve ring surrounding the isthmus. Excretory pore at or just posterior to basal bulb. Cardia reduced, surrounded by intestinal tissue. Reproductive system absent. Rectum poorly visible. Anus closed. Tail conoid elongate with acute tip without mucro. Terminal hyaline part 37 to 54% of tail length (Figs. 1–4).
Infective non-sheathed juveniles (J3 stage)
Body 0.47 to 0.56 mm long, with habitus almost straight after fixation. Cuticle with only transversal striae. Lip region lacking differentiate lips, and labial and cephalic papillae not visible. Oral opening rounded, closed, bearing a large, very refringent dorsal tooth. Amphidial apertures very prominent. Stoma tubular, slightly longer than the lip region wide. Pharynx, nerve ring and excretory pore location similar to the sheathed stage. Cardia reduced, surrounded by intestinal tissue. Rectum poorly visible. Anus closed. Tail conoid with very acute tip without mucro. Terminal hyaline part absent (Figs. 1–4).
Diagnosis of Heterorhabditis ruandica n. sp. and morphological relationships with other species
Heterorhabditis ruandica n. sp. is characterized by having hermaphrodite females 2.91 to 4.12 mm long, amphimictic females 1.13 to 1.61 mm long, males 0.65 to 0.86 mm long, and IJs 0.50 to 0.59 mm long. Cuticle with poorly visible annuli in adults, with longitudinal crests in IJ2 and with well-developed annuli in IJ3. Lip region with six low lips having thin and acute lipplets in adults. Lips are poorly developed in IJ2 and bearing a large refringent dorsal tooth in IJ3. Stoma reduced in adults and tubular in IJ. Pharynx robust and short in adults and narrow and slender in IJ. Female reproductive system didelphic–amphidelphic. Anal body diameter in hermaphrodites 29 to 51 µm long, in amphimictic females 18 to 34 µm long, and in males 15 to 25 µm long. Tail short and conoid with acute terminus at cellular part in hermaphrodite females (63-98 µm long, c = 34-51, c′ = 1.7-2.6) and slightly bifurcated in amphimictic females (62-88 µm long, c = 16-24, c′ = 1.9-3.6). Tail conoid-elongate in IJ2 (49-64 µm long, c = 8.0-12, c′ = 3.1-6.2) and IJ3 (22-39 µm long, c = 7.6-8.6, c′ = 3.4-5.8). Male reproductive system monorchic, with spicules 34 to 50 µm long having reduced manubrium 15 to 21 µm long, bursa peloderan bearing nine pairs of genital papillae (1 + 2/3 + 3).
Heterorhabditis ruandica n. sp. is morphologically similar to H. egyptii, H. bacteriophora, H. georgiana, and H. beicherriana, and can be distinguished from these species mainly by adult and infective juvenile characters (Tables 3–6). Heterorhabditis ruandica n. sp. can be distinguished from H. egyptii by the distance from the anterior end to the nerve ring in IJs (52-64 vs. 78-100 µm), the presence of a cephalic tooth in IJs (large vs. apparently small or absent). Additionally, hermaphroditic females of these two species differ in size (2.91-4.12 vs. 2.10-3.10), body diameter (209-274 vs. 107-164 µm), and in the distance from the anterior end to the excretory pore (106-153 vs. 154-205). Amphimictic female of H. ruandica n. sp. and H. egyptii differ in the size of their tails (62-88 vs. 56-78 µm).
Heterorhabditis ruandica n. sp. IJs can be distinguished from the IJs of H. bacteriophora by the distance between the anterior end and the excretory pore (67-90 vs. 87-110 µm) and the distance between the anterior end and the nerve ring (52-64 vs. 72-93 µm), and by the tail length (49-65 vs. 83-112 µm). The males of Heterorhabditis ruandica n. sp. can be distinguished from the males of H. bacteriophora by the distance from the anterior end to the excretory pore (61-109 vs. 114-130 µm) and by the lower D% value (61-97 vs. 117). Hermaphroditic and amphimictic females also show various morphometric differences (Tables 3–6).
Heterorhabditis ruandica n. sp. can be distinguished from H. beicherriana by the size of IJs (496-591 vs. 566-687 µm), the distance between the anterior end and the excretory pore (67-90 vs. 100-122 µm) and between the anterior end and the nerve ring (52-64 vs. 85-106 µm), and by neck (103-131 vs. 118-146 µm) and tail lengths of IJs (49-65 vs. 86-111 µm). The body length of H. ruandica n. sp. males is shorter than the body length of H. beicherriana males (652-863 vs 889-1192 µm). Males can also be distinguished by body diameter (40-51 vs. 51-73 µm), and by the distance between the anterior end and the excretory pore (61-109 vs. 130-157 µm) and between the anterior end and the nerve ring (56-74 vs. 81-108 µm), by neck (84-117 vs. 116-143 µm), tail (21-29 vs. 32-45 µm) and gubernaculum (15-21 vs. 22-27 µm) lengths, and by the D% value (61-97 vs. 102-120). Several other morphometric differences were also observed in hermaphroditic and amphimictic females (Tables 3–6).
Heterorhabditis ruandica n. sp. can be distinguished from H. georgiana by the anterior end to the excretory pore (67-90 vs. 97-113 µm) and the anterior end to the nerve ring (52-64 vs. 74-94 µm) distances, and by the tail length (49-65 vs. 86-108 µm) of IJs. The males can be distinguished by the anterior end to the excretory pore (61-109 vs. 101-145 µm) and the anterior end to the nerve ring (56-74 vs. 72-93 µm) distances, and by tail (21-29 vs. 29-41 µm) and gubernaculum (15-21 vs. 20-28 µm) length, and by D% values (61-97 vs. 100-122). Several other morphometric characters of hermaphroditic and amphimictic females differ between these two species (Tables 3–6).
Type host and locality
The type hosts are unknown as the nematodes of this genus can be hosted by different insect species and were isolated from soil samples by the Galleria baiting technique (Bedding and Akhurst, 1975; White, 1927). Nematode strains H. ruandica n. sp. Rw18_M-Hr1a and Rw18_M-Hr1b were collected in the district of Karongi, Western province of the Republic of Rwanda (Decimal degrees coordinates: -2.131500, 29.325467) in a moist habitat along a river bench covered with sweet potato plants. Heterorhabditis ruandica n. sp. Rw14_N-C4a nematodes were collected in a ploughed cropland on terraces in a hilly area near Kanyirandori village, Tare sector, Nyamagabe district, Southern province of the Republic of Rwanda (Decimal degrees coordinates: -2.500000, 29.483333).
Type material
Rw14_N-C4a nematodes are the type material for Heterorhabditis ruandica n. sp. Holotype male, and 15 paratype hermaphrodites, males and amphimictic females and 15 third stage juveniles were deposited in the National Nematode Collection of India, IARI, New Delhi, India. Additional specimens were deposited at the nematode collection of the Department of Animal Biology, Plant Biology and Ecology of the University of Jaén, Spain, under the following slide numbers: Rwa001-01 to -12 (25 hermaphrodite females and 6 juveniles), Rwa002-01 to -05 (8 amphimictic females and 9 males), and Rwa003-01 to -02 (8 juveniles). Nematode cultures are maintained in the Institute of Biology, University of Neuchatel, Switzerland and in the Rwanda Agriculture and Animal Resource Development Board, Rubona, Rwanda.
Etymology
The specific name refers to the country, the Republic of Rwanda (Africa), where the type material, Heterorhabditis ruandica n. sp. Rw14_N-C4a nematodes, used to phenotypically characterize the species, were collected.
Body 0.81 to 0.91 mm long, J-shaped after heat killing and body arcuate posteriorly. Cuticle almost smooth, with transversal striae poorly developed. Lateral field not visible. Lip region with six lips poorly developed bearing six acute labial papillae at oral margin and four rounded cephalic papillae at base of lips. Oral opening almost rounded with thick margin. Amphidial apertures pore-like, ovoid and located posterior to lateral labial papillae. Stoma rhabditoid type, 0.9 to 1.6 times the lip region width, with short cheilostom with poorly refringent rounded cheilorhabdia, short gymnostom with refringent bar-like rhabdia, and long stegostom surrounded by the pharyngeal collar and bearing bar-like pro-mesorhabdia and small poorly refringent meta-telorhabdia. Pharynx poorly developed with robust corpus without differentiated metacorpus, short and slightly narrow isthmus and robust pyriform bulb with poorly visible valvular apparatus. Nerve ring encircling the isthmus at 61% to 96% of neck length, just anterior to basal bulb. Excretory pore located at or posterior to the basal bulb, located at 78% to 134% of neck length. Cardia poorly developed, surrounding by intestinal tissue. Intestine without differentiations. Genital system monorchic, laterally reflexed. Spicules well developed, separate, with more or less rounded manubrium, calamus poorly developed, and thinner and slender lamina with acute tip, scarcely prominent dorsal hump and poorly developed ventral velum. Gubernaculum with manubrium slightly ventral curved and straight corpus, 40% to 60% of spicule length. Tail conoid with acute tip, ventrally curved posteriorly, flanked by the bursa. Bursa peloderan, with nine pairs of genital papillae (1 + 2/3 + 3), one of them probably the phasmid: three pairs pre-cloacal (GP1–GP3) and six pairs post-cloacal being three pairs at mid tail length (GP4–GP6) and three pairs (GP7–GP9) terminal; GP1 and GP2 more spaced, GP2 and GP3 closely spaced (Figs. 5–8).
Hermaphroditic females
Body 4.41 to 6.18 mm long, arcuate with general morphology similar to male, having labial papillae more acute and prominent. Genital reproductive system didelphic–amphidelphic with ovaries well developed, reflexed, oviducts and uteri not well visible, vagina very short and vulva small having transverse slit opening. Rectum slender, about 1.5 times longer than the anal body diameter. Anus with prominent lips. Tail conoid with acute tip lacking mucro, having cellular part simple at its junction with the hyaline part. Phasmids inconspicuous (Figs. 5–8).
Amphimictic females
Body similar to, but usually smaller than hermaphrodites, 1.95 to 2.80 mm long. Rectum very long, about twice longer than the anal body diameter. Anus with posterior lip more prominent. Tail conoid with acute tip lacking mucro, having cellular part simple at its junction with the hyaline part (Figs. 5–8).
Infective sheathed juveniles (J3 stage envolved by the J2 stage cuticle)
Body 0.49–0.58 mm long, with habitus slightly ventral curved after fixation. Cuticle with transversal striae at anterior end, with both transversal and longitudinal striae at neck region and only with longitudinal striae at rest of body. Lip region lacking differentiate lips, bearing six labial papillae and cephalic papillae not visible. Amphidial apertures very reduced. Oral opening closed, having triradial symmetry. Stoma tubular, about twice the lip region wide. Pharynx slender, with long and narrow corpus, very narrow isthmus and pyriform basal bulb. Nerve ring surrounding the isthmus. Excretory pore at or just posterior to basal bulb. Cardia reduced, surrounded by intestinal tissue. Reproductive system absent. Rectum poorly visible. Anus closed. Tail conoid elongate with acute tip without mucro. Terminal hyaline part 31% to 56% of tail length (Figs. 5–8).
Infective non-sheathed juveniles (J3 stage)
Body 0.47 to 0.55 mm long, with habitus slightly ventral curved after fixation. Cuticle with only transversal striae. Lip region lacking differentiate lips, and labial and cephalic papillae not visible. Oral opening rounded, closed, bearing a small dorsal tooth. Amphidial apertures very prominent. Stoma tubular, slightly longer than the lip region wide. Pharynx, nerve ring and excretory pore location similar to the sheathed stage. Cardia reduced, surrounded by intestinal tissue. Rectum poorly visible. Anus closed. Tail conoid with acute tip without mucro. Terminal hyaline part absent (Figs. 5–8).
Diagnosis of Heterorhabditis zacatecana n. sp. and relationships with other species
Heterorhabditis zacatecana n. sp. is characterized by having hermaphrodite females 4.41 to 6.18 mm long, amphimictic females 1.9 to 2.7 mm long, males 0.81 to 0.91 mm long, and IJs 0.49 to 0.57 mm long. Cuticle with poorly visible annuli in adults, with longitudinal crests in IJ2 and with well-developed annuli in IJ3. Lip region with six low lips having thick and acute lipplets in adults. Lips are poorly developed in IJ2 and bearing a small refringent dorsal tooth in IJ3. Stoma reduced in adults and tubular in IJs. Pharynx robust and short in adults, and narrow and slender in IJs. Female reproductive system didelphic–amphidelphic. Anal body diameter in hermaphrodites 34 to 58 µm long, in amphimictic females 31 to 41 µm long, and in males 13 to 22 µm long. Tail short and conoid with acute terminus at cellular part in hermaphrodite females (63-87 µm long, c = 52-90, c′ = 1.2-2.4). and in amphimictic females (45-75 µm long, c = 31-63, c′ = 1.3-2.0). Tail conoid-elongate in IJ2 (52-63 µm long, c = 7.9–9.8, c′ = 4.0-6.5) and in IJ3 (25-34 µm long, c = 8.2-10.5, c′ = 4.3-6.7). Male reproductive system monorchic, with spicules 38 to 55 µm long having conoid manubrium 15 to 25 µm long, bursa peloderan bearing nine pairs of genital papillae (1 + 2/3 + 3).
Heterorhabditis zacatecana n. sp. is morphologically similar to H. ruandica n. sp., H. amazonensis, H. bacteriophora, H. georgiana and H. beicherriana, and can be distinguished from these species mainly by adults and infective juvenile characters (Tables 3–6). Heterorhabditis zacatecana n. sp. can be distinguished from H. ruandica n. sp., one of the morphologically most similar species, by the shape of the male spicule (slender vs. robust) and the manubrium size (large vs. small), the size of hermaphrodites (4.41-6.18 vs. 2.91-4.12 mm), the hermaphrodite neck length (174-231 vs. 134-159 µm), and the hermaphrodite c ratio (52-90 vs. 34-51). The size of amphimictic females (1.95-2.80 vs. 1.13-1.61 µm), the shape of the tail tip (acute and longer vs. with mucro), the type of cellular–hyaline junction part (simple vs. bifurcated), the body diameter (160-228 vs. 68-83 µm), the a (11-15 vs. 15-20), b (16-21 vs. 9-14), and c ratios (31-63 vs. 16-24) and the anal body diameter (31-41 vs. 18-34 µm) differ also between H. zacatecana n. sp. and H. ruandica n. sp. IJs anterior ends also differ between these two species (small vs. large), and the presence of a cephalic tooth (small or absent vs. refringent and large).
Morphologically, the IJs of H. zacatecana n. sp. can be distinguished from the IJs of H. amazonensis by their size (493-578 vs. 567-612 µm), the distance from the anterior end to the nerve ring (59-72 vs. 76-93 µm), the neck length (96-124 vs. 107-132 µm), the tail length (52-63 vs. 98-115 µm), the a (19-24 vs. 24-29), c (8.2-10.5 vs. 5.1-6.1), and c′ (4.3-6.7 vs. ca. 7.3 µm) ratios and the E% (128-184 vs. 89-109). Moreover, hermaphroditic females differ in body size (4.41-6.12 vs. 3.52-5.59 mm), tail length (62-87 vs. 104-154 µm) and anal body diameter (34-58 vs. 59-85 µm). Amphimictic females of these two species differ also in body size (1.95-2.80 vs. 1.28-2.07 µm), tail length (45-75 vs. 25-38 µm), and body diameter (160-228 vs. 70-122). Male sizes differ between H. zacatecana n. sp. and H. amazonensis (0.81-0.91 vs. 0.69 vs. 0.83 mm) and body diameter (41-56 vs. 36-43 µm).
Heterorhabditis zacatecana n. sp. IJ can be distinguished from H. bacteriophora by the distance from the anterior end to the nerve ring (59-72 vs. 72-93 µm), and the tail length (52-63 vs. 83-112 µm). In the case of males, they differ in the distance from the excretory pore to the anterior end (77-109 vs. 114-130 µm) and in body diameter (41-56 vs. 38-46 µm). Several morphometric differences were also observed in hermaphrodites and amphimictic females (Tables 2–6).
Heterorhabditis zacatecana n. sp. IJs can be distinguished from H. beicherriana IJs by the distance from anterior end to the excretory pore (72-99 vs. 100-122 µm) and the distance from the anterior end to the nerve ring (59-72 vs. 85-106 µm), the tail length (52-63 vs. 86-111), values of a (19-24 vs. 24-29), c´ (4.3-6.7 vs. 6.0-7.4), and c (8.2-10.5 vs. 5.9-6.8) ratios, and the E% value (128-184 vs. 103-121). Males can be differentiated by differences in neck (71-108 vs. 116-143 µm) and tail (21-33 vs. 32-35 µm) lengths, the distance from the anterior end to the excretory pore (77-109 vs, 130-157 µm) and from the anterior end to the nerve ring (60-78 vs. 81-100 µm). Several morphometric differences were also observed in hermaphrodites and amphimictic females of these two species (Tables 2–6).
Heterorhabditis zacatecana n. sp. IJs can be distinguished from H. georgiana IJs by differences in body diameter (23-27 vs. 17-26 µm), tail length (52-63 vs. 86-108 µm), and anterior end to excretory pore (72-99 vs. 97-113 µm) and anterior end to nerve ring distances (59-72 vs. 74-94 µm). The a, b and c ratios, E% and D% of IJs differ also in these two species. The males of these two species differ in anterior end to excretory pore (77-109 vs. 101-145 µm) and anterior end to nerve ring distances (60-78 vs. 72-93 µm), and neck (71-108 vs. 100-122 µm) and tail (21-33 vs. 29-41 µm) lengths. Several morphometric characters of hermaphroditic and amphimictic females differ between these two species (Tables 2–6).
Type host and locality
The type hosts are unknown as the nematodes of this genus can be hosted by different insect species and were isolated from soil samples by the Galleria baiting technique (Bedding and Akhurst, 1975; White, 1927). Heterorhabditis zacatecana n. sp. MEX-39 and MEX-40 nematodes were collected in maize fields in Villanueva (Zacatecas, Mexico; decimal degrees coordinates: 22.161371, -102.887940), and Heterorhabditis zacatecana n. sp. MEX-41 nematodes were collected in maize fields in Apaseo el Alto (Guanajuato, Mexico; decimal degrees coordinates: 20.470774, -100.59571).
Type material
MEX-39 nematodes are the type material for Heterorhabditis zacatecana n. sp. Holotype male, 15 paratype and 15 third stage juveniles were deposited in the National Nematode Collection of India, IARI, New Delhi. Additional specimens were deposited in the nematode collection of the Department of Animal Biology, Plant Biology and Ecology of the University of Jaén, Spain, under the following slide numbers: Mex001-01 to -03 (6 hermaphrodite females), Mex002-01 to -04 (8 amphimictic females and 3 males), and Mex003-01 to -04 (14 juveniles). Nematode cultures are maintained in the Institute of Biology, University of Neuchatel, Switzerland.
Etymology
The specific name refers to the Mexican state, Zacatecas, where the type material, Heterorhabditis zacatecana n. sp. MEX-39 nematodes, used to phenotypically characterize the species were collected.
Cross-hybridization experiments
No progeny was observed when males and females of H. ruandica n. sp. Rw14_N-C4a and of H. zacatecana n. sp. MEX-39 were left to interact. No progeny was observed when males and females of H. ruandica n. sp. Rw14_N-C4a and of H. bacteriophora CH21 were left to interact. No progeny was observed when males and females of H. zacatecana n. sp. MEX-39 and of H. bacteriophora CH21 were left to interact. When males and females of H. ruandica n. sp. Rw14_N-C4a were crossed, fertile progeny was observed. When males and females of H. zacatecana n. sp. MEX-39 were crossed, fertile progeny was observed. When males and females of H. bacteriophora CH21 were crossed, fertile progeny was observed. Similarly, H. zacatecana n. sp. MEX-39 and H. zacatecana n. sp. MEX-40 nematodes produced fertile progeny, and H. ruandica n. sp. Rw18_M-Hr1a and H. ruandica n. sp. Rw14_N-C4a nematodes produced fertile progeny. These results provide further support for the heterospecific status of the Rwandan and the Mexican nematode populations.
Nematode molecular characterization and phylogenetic relationships
Phylogenetic reconstructions based on nuclear and mitochondrial genes (ITS, D2–D3, COI, umc-87, and cmd–1), either individually or concatenated, confirm that the nematodes of the genus Heterorhabditis are grouped into three major clades: the “Megidis-group”, the “Indica-group” and the “Bacteriophora-group”, which is consistent with previous studies (Dhakal et al., 2020) (Fig. 9, Fig. S1). The clade of the “Bacteriophora-group” is, in turn, separated into five subclades. Three of them are composed of already described species: H. beicherriana, H. georgiana, and H. bacteriophora, and two of them are composed of two new, undescribed species, which we named here H. zacatecana n. sp., and H. ruandica n. sp. (Fig. 9, Fig. S1). Clearer phylogenetic separations within the species of the clade of the “Bacteriophora–group” were observed when phylogenies were reconstructed based on COI, ITS, or on concatenated sequences of COI, ITS, and D2–D3 (Fig. 9, Fig. S1). Closer inspection at the ITS, D2–D3 and COI sequences reveals unambiguous genetic differences between the nematodes of the “Bacteriophora–group” (Fig. 10). Sequence similarity scores and nucleotide difference counts show a closer relationship between H. bacteriophora, H. ruandica n. sp., and H. zacatecana n. sp. nematodes (Fig. 11 and Figs. S2-S6). Heterorhabditis ruandica n. sp. and H. bacteriophora share 99.1% and differ in 6 nucleotide positions in the ITS sequences flanked by primers TW81 and AB28, share 99.8% and differ in 1 nucleotide position in the D2–D3 sequences flanked by primers D2A and D3B, and share 94.1 to 94.7% and differ in 18 to 19 nucleotide positions in the COI sequences flanked by primers HCF and HCR (Fig. 11 and Figs. S2-S6). Heterorhabditis zacatecana n. sp. and H. bacteriophora share 99.4% and differ in 4 nucleotide positions in the ITS sequences flanked by primers TW81 and AB28, share 99.8% and differ in 1 nucleotide position in the D2–D3 sequences flanked by primers D2A and D3B, and share 94.1 to 94.4% and differ in 19 to 20 nucleotide positions in the COI sequences flanked by primers HCF and HCR (Fig. 11 and Figs. S2-S6). Heterorhabditis ruandica n. sp. and H. zacatecana share 99.7% and differ in 2 nucleotide positions in the ITS sequences flanked by primers TW81 and AB28, share 100% and differ in no nucleotide position in the D2–D3 sequences flanked by primers D2A and D3B, and share 97.6% to 98.2% and differ in 6–8 nucleotide positions in the COI sequences flanked by primers HCF and HCR (Fig. 11 and Figs. S2-S6). Noteworthy, we observed almost no intraspecific variations within the nematodes of the “Bacteriophora-group” at different genetic loci (Figs. 10, 11, and Figs. S2–S6). However, the sequences of the COI gene show very interesting signatures of population–specific polymorphism (Figs. 10D-F, 11). Specifically, Heterorhabditis ruandica n. sp. Rw18_M-Hr1a and Rw18_M-Hr1b nematodes that were collected in the same western Rwandan region differ from the Heterorhabditis ruandica n. sp. Rw14_N-C4a nematodes collected in a southern Rwandan region in a transitional nucleotide change (g.1212A > G) (Fig. 10D). Moreover, H. zacatecana n. sp. MEX-39 and MEX-40 nematodes collected in north-central Mexico and H. zacatecana n. sp. MEX-41 nematodes collected in central Mexico differ in three transitional nucleotide changes (g.1257T > C, g.1324T > C, and g.1464A > G) (Fig. 10D-F). Hence, due to its highly conserved species–specific polymorphism, and the consistent population-specific polymorphic patterns, the COI gene emerges as an important phylogenetic marker also for the genus Heterorhabditis, in a similar manner as it is for many other taxonomic groups (Hebert et al., 2003; Pentinsaari et al., 2016).
Interspecific genetic variability within the H. bacteriophora clade
In a recent study, Dhakal et al. (2020) studied several hundreds of ITS sequences of Heterorhabitis nematodes and recognized that nematodes identified as H. bacteriophora are represented by at least three haplotypes, some of which, the authors hypothesized, could represent new species. We contrasted their results and the sequences they used with the sequences we generated and found out that H. bacteriophora DE2, DE6, EN01, HP88, IT6, and PT1 nematodes represent Dhakal’s haplotype 1 (Figs. S7 and S8). Dhakal’s haplotype 2 is actually a mixture of two haplotypes: one represented by H. zacatecana n. sp. MEX-39, MEX-40, MEX-41, and by strains N2 and MK with identical ITS sequences; and one represented by strains UP2A2, 267, 269, 270, 271, 275, and 276 with identical ITS sequences. Strains MEX-39, MEX-40, MEX-41, N2 and MK differ in a transitional nucleotide change (g.2049A > G) with strains UP2A2, 267, 269, 270, 271, 275, and 276. Hence, strains N2, MK are likely H. zacatecana, and strains UP2A2, 267, 269, 270, 271, 275, and 276 might represent a new species. However, full characterization is needed to prove this hypothesis. Dhakal’s haplotype 3, represented by NGPS20, among others isolates, might also represent a new species, but again full characterization is needed to prove this hypothesis. In addition, our analyses reveal what we call a fourth haplotype, to follow Dhakal’s system, which is represented by H. ruandica Rw18_M-Hr1a, Rw18_M-Hr1b, and Rw14_N-C4a nematodes (Figs. S7 and S8). Phylogenetic reconstructions show a clear phylogenetic separation between all these haplotypes (Fig. S8). Hence, some of the haplotypes described by Dhakal et al. (2020) represent new species, closely related to H. bacteriophora, and some others likely represent new species, which highlights the power of statistical parsimony network analyses to uncover undescribed species of the genus Heterorhabditis, and supporting previous hypothesis regarding the taxonomic status of these nematode isolates (Bruno et al., 2020; Dhakal et al., 2020; Fallet et al., 2020).
Symbiotic relationships
Up to now, the bacterial genus Photorhabdus Boemare, Akhurst and Mourtant 1993 contains 27 taxa, including species and subspecies (Machado et al., 2021b). Phylogenetic relationship reconstructions based on whole genome sequences show that the bacterial symbionts isolated from H. zacatecana n. sp. MEX-39 and H. ruandica n. sp. Rw14_N-C4a nematodes, named here as MEX-39 and RW14-46, respectively, show high similarity with two of the already described Photorhabdus species: Photorhabdus kleinii and P. laumondii subsp. laumondii, respectively (Fig. 12). Photorhabdus kleinii MEX-39 shares 87–88% digital DNA–DNA hybridization (dDDH) with other members of the same species, while P. laumondii subsp. laumondii RW14-46 shares 89% digital DNA–DNA hybridization (dDDH) with other members of the same species, (Fig. S9).
On the synonymization and declaration of species inquirendae of some species
We revised the original publications of all synonymized species and based on their morphology and molecular data (when available), we reinforce the synonymized status of most of them (Khan et al., 1976; Wouts, 1979; Stock, 1993; Gardner et al., 1994; Liu, 1994; Stock et al., 1996; Plichta et al., 2009; Stock et al., 2009; Maneesakorn et al., 2015; Hunt and Nguyen, 2016; Shahina et al., 2017; Dhakal et al., 2020). However, the original description of H. bacteriophora provided by Poinar (1976) shows males with very anterior GP1 while in its synonymized species H. heliothidis (Khan, Brooks & Hirschmann, 1976) Poinar, Thomas & Hess, 1977 (=Chromonema heliothidisKhan et al., 1976) the GP1 appears more posterior (Khan et al., 1976; Poinar, 1976). Hence, it is likely that both species are not conspecific. Therefore, we declare H. heliothidis (Khan, Brooks & Hirschmann, 1976) Poinar, Thomas & Hess, 1977 as species inquirenda. Heterorhabditis hoptha and H. poinari were poorly described (Turco, 1970; Kakulia and Mikaia, 1997). Original descriptions lack differentiated description of all diagnostic characters of adult and larval stages. According to this, both species should remain in the list of species inquirendae. Heterorhabditis egyptii and H. hambletoni were described showing all diagnostic characters of adults and larvae stages. According to this, both species are considered valid herein (Pereira, 1937; Abd-Elgawad and Ameen, 2005). The lack of molecular data, however, impairs their inclusion in future phylogenetic studies. Nevertheless, new species description should contrast morphological characters with these species. An updated dichotomous key to identify the species of the genus Heterorhabditis is provided (Fig. 13, Tables 3-6).
On the species of the genus Heterorhabditis
Considering the results of this study and the analyses of all the literature that describes new species of the genus Heterorhabditis, the updated list of the species of the genus, including their status, is as follows.
H. beicherriana Li, Liu, Nermut, Půža & Mráček, 2012
H. egyptiiAbd-Elgawad & Ameen, 2005. This species was declared species inquirenda by Nguyen and Hunt (2007) but considered valid by Sudhaus (2011). As this species was described showing all diagnostic characters of adults and larvae stages, and it is morphologically distinct from all the other valid species, this species is also considered valid herein. The lack of molecular data, however, impairs its inclusion in future phylogenetic studies. Nevertheless, new species description should contrast morphological characters with this species.
= Rhabditis hambletoniPereira, 1937. This species was described showing all diagnostic characters of adults and larvae stages. It was transferred to the genus Heterorhabditis by Poinar (1976). As this species was described showing all diagnostic characters of adults and larvae stages, and it is morphologically distinct from all the other valid species, this species is considered valid herein. The lack of molecular data, however, impairs its inclusion in future phylogenetic studies. Nevertheless, new species description should contrast morphological characters with this species.
= Heterorhabditis brevicaudisLiu, 1994. Several important diagnostic characters are missing and no molecular data are provided in the description of this species, although, it appears to be morphologically different from H. downesi, H. baujardi, and H. mexicana (Stock et al., 2002; Phan et al., 2003; Nguyen et al., 2004). Perhaps due to this reason, it was declared species inquirenda by Nguyen and Hunt (2007). A nematode population that shares several morphological characters with the original population used to describe the species was characterized more recently (Hsieh et al., 2009). ITS sequences are almost identical to the sequences of H. indica, justifying its synonymization (Hunt and Nguyen, 2016; Dhakal et al., 2020).
= Heterorhabditis pakistanensis Shahina, Tabassum, Salma, Mehreen & Knoetze, 2016. Synonymized by Hunt and Nguyen (2016) based on the minor ITS sequence divergencies between Heterorhabditis pakistanensis and H. indica. Synonymisation status is further supported by molecular data analyses carried out by Dhakal et al. (2020).
This species was poorly described. The original description lacks differentiated description of all diagnostic characters of adult and larval stages. According to this, this species should remain on the list of species inquirendae.
H. poinari Kakuliya and Mikaia, 1997. This species was poorly described. The original description lacks differentiated description of all diagnostic characters of adult and larval stages. According to this, this species should remain on the list of species inquirendae.
This species was synonimized by Akhurst 1987 based on differential electrophoretic patterns of nematode lysates. However, the original description of H. bacteriophora carried out by Poinar (1976) shows males with very anterior GP1, while in its synonymized species H. heliothidis (Khan, Brooks & Hirschmann, 1976) Poinar, Thomas & Hess, 1977 the GP1 appears more posterior (Khan et al., 1976; Poinar, 1976). Probably both species are not conspecific. We therefore declare H. heliothidis (Khan, Brooks & Hirschmann, 1976) Poinar, Thomas & Hess, 1977 species inquirenda.
The results of our study uncover the low levels of interspecific variation in some regions of the rRNA genes, especially in the D2–D3 expansion segments of the 28S rRNA, and also uncover the almost absent intraspecific variation of these sequences in the nematodes of the “Bacteriophora-group”. Mitochondrial genes such as COI provide better phylogenetic resolutive power, even at the population level, highlighting their great potential for the taxonomic characterization of closely related species of the genus Heterorhabditis. The threshold for species delimitation using COI sequences has been proposed to be around 94% (Pentinsaari et al., 2016). Using this threshold, we can clearly assign the Mexican and the Rwandan nematodes to new taxa within the “Bacteriophora group”. However, the sequence similarity scores of the Mexican and the Rwandan nematodes is between 97.6% and 98.2%. These scores are higher than the proposed 94% threshold but are consistent across nematode isolates and significantly lower than the intraspecific variations, prompting the question of whether the Rwandan and the Mexican nematodes should be classified into two different species, or into the same. Based on the results of the self-crossing and cross-hybridization experiments, and on the evident morphological and morphometric differences between these two groups of nematodes, we conclude that they indeed represent two distinct biological species. Thus, the boundary that delimits species in the genus Heterorhabditis is around 97% to 98% sequence similarity in the COI genomic sequence, and the Rwandan and the Mexican nematodes represent two new species, Heterorhabditis ruandica n. sp and H. zacatecana n. sp.
Supplementary Material
Supplementary figures and tables can be retrieved from: https://doi.org/10.5281/zenodo.5614704