Defining, cataloguing, mapping and preserving biodiversity is generally accepted to be one of the key challenges for the 21st Century. The Kingdom of Morocco possesses the richest and most varied herpetofauna of the Maghreb and the western Mediterranean, is characterized by high richness of reptiles, endemism and European relict species (del Mármol et al., 2019). Phylogenetic analyses performed over recent decades have identified notable evolutionary lineages in the Moroccan herpetofauna, several of which represent new (cryptic) species or species complexes (e.g. Barata et al., 2012; Salvi et al., 2018), demonstrating the value of incorporating molecular tools into diversity assessments. Still, despite the well-known biodiversity of reptiles, parasite diversity associated with these hosts remains poorly known. Part of the problem in estimating parasite diversity has historically been that, especially for endoparasites, identification based on morphological characters and life-cycle traits is often difficult. Parasites often present a simplified morphology (Jorge et al., 2011), and while DNA sequencing approaches should help overcome these problems, molecular studies have lagged behind those of free-ranging organisms (Criscione et al., 2005). However, application of molecular tools is starting to show that some parasites currently considered as a single species actually consist of genetically different lineages or cryptic species (e.g. Jorge et al., 2012, 2013a).
Reptiles are parasitized by various helminth species, typically occurring in depauperate communities (Aho, 1990), possibly due to characteristics of the hosts such as the simplicity of the alimentary canal, their low vagility, a nonspecialized diet, and characteristics of the parasites such as direct life cycles (Roca & Hornero, 1994). Various studies have suggested that the typical helminth infection pattern in reptiles is that few species occur frequently, while many species are rare (Birlik et al., 2015). In studies of reptiles from the Iberian Peninsula and the North of Africa, helminth fauna of assessed lizards was poor, and mainly composed by members of the family Pharyngodonidae that are often detected in insectivorous reptiles (Chabaud & Golvan, 1957; Ibrahim et al., 2005; Carretero et al., 2011; Roca et al., 2020). These are usually identified on the basis of male morphology, since females are generally similar among species (Jorge et al., 2014). Although this family includes 21 genera, only a few of these are typically found in insectivorous lizards, with
Recent nematode surveys of reptiles in the Mediterranean region resulted in several new host records and descriptions of new species (Jorge et al., 2011, 2013a, 2014). Despite the high diversity of reptiles in Morocco, parasites from these hosts remain poorly studied, with a few incidental assessments as part of surveys at the Canary Islands (e.g. Jorge et al., 2011, 2018). In this sense, no survey including molecular data has been performed for nematodes of reptiles from Morocco. Indeed, with one or two exceptions (e.g. Carretero et al., 2011), there have been few molecular studies in nematodes from reptiles across the whole of the Maghreb. The aim of this study was to fill this gap, by assessing diversity of Pharyngodonidae nematodes from six lizard species from Morocco. By sequencing exemplars for two molecular markers (fragments of the 18S rRNA and Cytochrome Oxidase 1 gene), we aim to place these in a phylogenetic framework, to potentially identify cryptic forms, and to examine patterns of host-specificity.
Helminths were collected from pellets, which were obtained through spontaneous defecation or by gentle abdominal massage of six lizard species:
Fig. 1
Map of Northern Morocco with the geographic locations of sampled nematodes. Their respective hosts species are in Table 1.

Fecal samples from which parasitic nematodes were recovered and included in the genetic analysis.
Code | Nematode species | Host species | Locality | 18S | CoI |
---|---|---|---|---|---|
26 | Ounagha, Essaouira | ||||
32 | Dar Bouzza, Casa | OP548559 | OP558784 | ||
111 | Oukaïmeden, Marrakech | OP548558 | |||
112 | Oukaïmeden, Marrakech | OP548557 | |||
131 | Oukaïmeden, Marrakech | OP548556 | |||
164 | Oukaïmeden, Marrakech | OP548555 | |||
205 | Fom Jrana, Chichaoua | OP548554 | |||
219 | Ait El Qaq, Marrakech | OP548553 | |||
234 | Aït Aammour ou Ali, Azrou | OP548552 | OP558785 | ||
283 | Ait El Qaq, Marrakech | OP548551 | |||
377 | Imlil, Marrakech | OP548550 | |||
1149 | Admin Forest, Agadir | OP548549 | |||
1213 | El Ghazoua, Essaouira | OP548548 | |||
1447 | Bigoudine, Argana | OP548547 | OP558786 | ||
1476 | Oued Tensift, Marrakech | OP548546 | OP558787 | ||
1480 | Oued Tensift, Marrakech | OP548545 | |||
1550 | Oued Tensift, Marrakech | OP548544 |
Extractions of genomic DNA were performed using individual nematodes, according to the saline method (Maia et al., 2014). Two partial genes were amplified: the nuclear 18S rRNA (18S) gene and the mitochondrial Cytochrome c Oxidase subunit I (COI). The COI fragment was amplified using the primers LCO and HCO from Folmer et al. (1994), while the 18S was amplified using Nem 18SF and Nem 18SR from Floyd et al. (2005). Polymerase chain reactions were performed in a total volume of 15 μl, consisting of PCR buffer at 1 × concentration; MgCl2 at 1.5 mM; dNTPs at a concentration of 0.2 mM for each nucleotide; primers at 0.5 μM each; BSA at 0.4 μg/μl (bovine serum albumin) (Roche Applied Science); and Taq DNA polymerase (Invitrogen Corporation) 0.025 units/μl and 1 μl of DNA template. PCR conditions were 35 cycles of: 30 sec at 94°C, 30 sec at 50 – 54°C and 1 min at 72°C, with an additional denaturation step of 3 min at 94°C and ending with a final extension at 72°C for 10 min. Successful amplified products, confirmed through electrophoresis, were cleaned and sequenced by a commercial facility (Gene Wiz, Germany).
Sequences obtained were compared with those from GenBank using BLAST to confirm the taxonomic identity of the amplified products. Related sequences of
The research related to animal handling complied with all the relevant national regulations and institutional policies for the care and use of animals. The authorization for sampling of wild animals was granted by Cadi Ayyad University, Marrakech, Morocco. A field-work permit was issued by “Haut-Commissariat aux Eaux et Forêts et à la Lutte Contre la Désertification (HCEFLCD)”.
Three different genera were identified using microscopy. Identification of
Following previous descriptions (Lucker, 1952; Skrjabin et al., 1967; Ashour et al., 1992; Amer & Bursey, 2008; Mašová et al., 2009; Pereira et al., 2017; Er-Rguibi et al., 2022) samples from hosts of the genus
In principle, the geographic distribution of parasite diversity is predicted to match that of host diversity (Jorge & Poulin, 2018). Under this expectation, Morocco with its high diversity of reptile species should also harbor extensive nematode diversity within these hosts. However, this diversity remains essentially unassessed. In this study, we made a preliminary assessment of nematodes in pellets from six reptile host species. While it is known that there is a significantly lower detectability of nematodes from pellets compared to studies of intestines (Jorge et al., 2013b), this noninvasive approach can be used to give some baseline data on the parasites occurring in reptiles in this region.
As expected, the nematodes identified all belonged to three genera,
Fig. 2
Estimate ot relationships derived from the 18S rRNA gene sequences using a Bayesian approach. Values above branches represent Bayesian posterior probabilities and those below represent ML bootstrap support values (both given as percentages). Specimen code descriptions are given in Table 1. * These samples are considered as belonging to Patapharyngodon following de Sousa et al. (2018), but are listed on GenBank as

Fig. 3
Estimate of relationships derived from the CO1 mitochondrial gene sequences using a Bayesian approach. Values above branches represent Bayesian posterior probabilities and those below represent ML bootstrap support values (both given as percentages). Specimen code descriptions are given in Table 1.

Using universal 18S rRNA primers, we were able to amplify and sequence a region of this gene for all but one of our samples. The specimens of
While the distinction of
Although we only obtained 4 sequences from
To conclude, our new data clarifies several aspects regarding the presence of nematodes in reptiles from Morocco. Each lineage of parasite is quite specific to a related group of hosts, with
Fig. 1

Fig. 2

Fig. 3

Fecal samples from which parasitic nematodes were recovered and included in the genetic analysis.
Code | Nematode species | Host species | Locality | 18S | CoI |
---|---|---|---|---|---|
26 | Ounagha, Essaouira | ||||
32 | Dar Bouzza, Casa | OP548559 | OP558784 | ||
111 | Oukaïmeden, Marrakech | OP548558 | |||
112 | Oukaïmeden, Marrakech | OP548557 | |||
131 | Oukaïmeden, Marrakech | OP548556 | |||
164 | Oukaïmeden, Marrakech | OP548555 | |||
205 | Fom Jrana, Chichaoua | OP548554 | |||
219 | Ait El Qaq, Marrakech | OP548553 | |||
234 | Aït Aammour ou Ali, Azrou | OP548552 | OP558785 | ||
283 | Ait El Qaq, Marrakech | OP548551 | |||
377 | Imlil, Marrakech | OP548550 | |||
1149 | Admin Forest, Agadir | OP548549 | |||
1213 | El Ghazoua, Essaouira | OP548548 | |||
1447 | Bigoudine, Argana | OP548547 | OP558786 | ||
1476 | Oued Tensift, Marrakech | OP548546 | OP558787 | ||
1480 | Oued Tensift, Marrakech | OP548545 | |||
1550 | Oued Tensift, Marrakech | OP548544 |
Prevalence of ruminant paramphistomosis and comparative histopathology of the infected rumens in Narowal district, Punjab, Pakistan Inducing systemic acquired resistance (SAR) against root-knot nematode Meloidogyne javanica and evaluation of biochemical changes in cucumber rootEnterobius vermicularis infection: a cross-sectional study in preschool and school children in the North-Western part of SloveniaAdvantages and limitations of microscopy and molecular detections for diagnosis of soil-transmitted helminths: An overview First report of Aphelenchoides bicaudatus (Imamura, 1931 ) Filipjev and Schuurmans Stekhoven, 1941 associated with grass in South AfricaEnterobius vermicularis infections in IraqDifferential activity of human leukocyte extract on systemic immune response and cyst growth in mice with Echinococcus multilocularis infection after oral, subcutaneous and intraperitoneal routes of administrationNew Parasitic records of Aves: Phasianidae ( Alectoris chukar ) in Malakand division northern, PakistanEpidemiological survey on gastrointestinal and pulmonary parasites in cats around Toulouse (France) A large adnexal tumor caused by Enterobius vermicularis mimicking malignancy