Skip to main content

Phylogeography and genetic diversity of the copepod family Cyclopidae (Crustacea: Cyclopoida) from freshwater ecosystems of Southeast Nigeria



Copepods are key components of aquatic ecosystems and can help regulate the global carbon cycle. Much attention has been paid to the species diversity of copepods worldwide, but the phylogeography and genetic diversity of copepods in Nigeria is unexplored.


Using a mitochondrial cytochrome c oxidase subunit I marker, we preformed phylogenetic and phylogeographic analyses for Cyclopidae copepods in Southeast Nigeria. A high species diversity of Cyclopidae in Nigeria: 5 species of Tropocyclops, 5 species of Mesocyclops and 2 species of Thermocyclops from Cyclopidae were identified in 15 populations. Moreover, we detected 18 unique haplotypes, which fell into two distinct clades. Pairwise genetic distances (uncorrected p-distances) among the species of Cyclopidae ranged from 0.05 to 0.257. Several species co-existed in the same lake, and some haplotypes were shared among different geographic populations, suggesting a dispersal of Cyclopidae in our sampling region. Finally, we found that the population genetic diversity for each species of Cyclopidae was low in Nigeria.


Our findings explored the species diversity and distribution of copepods within the family Cyclopidae for 15 Nigerian freshwater ecosystems: a high species diversity of Cyclopidae copepods was detected over a small geographic sampling range. Results from this study contribute to a better understanding of copepod diversity of Nigerian freshwater ecosystems.


Copepods are one of the most taxonomically diverse groups of crustaceans, containing approximately 14,000 described species globally [1]. Copepods can be found in most kinds of aquatic habitats because of their remarkable evolutionary adaptability [1, 2]. They are key components in aquatic ecosystems, playing an important role in food webs [3, 4] and living as endo- or ectoparasites associated with aquatic animals [2, 5, 6]. Many previous studies have shown that copepods are sensitive to climate change [7, 8], because the range of copepods could track the rate of climate change [7]. Copepods can also help regulate the global carbon cycle [9, 10], and they can be used as indicators to natural and anthropogenic environmental stressors by tracing their responses to the elevation of atmospheric CO2 levels [11]. Thus, much attention has been paid to the bio-diversity of copepods in aquatic ecosystems [12, 13].

Copepods are the intermediate hosts of the parasitic nematode Dracunculus medinensis, which causes a serious Guinea-worm disease in Nigeria and elsewhere [14]. Humans could become infected by drinking unfiltered water containing copepods which are infected with larvae of D. medinensis. Therefore, most studies on copepods from Nigeria have focused on their role in the dispersal of the pathogen [15, 16]. Only a few regional biogeographic studies have been conducted on copepods based on morphological species identification [17]. For example, based on the morphology, a previous study showed the occurrence of the genera Mesocyclops Sars, 1914 and Thermocyclops Kiefer, 1927 in Nigerian freshwater ecosystems: six Mesocyclops species and three Thermocyclops species were identified [18]. Moreover, it was believed that M. aspericornis was one of the most abundant species of Mesocyclops in Nigerian waterbodies, and T. decipiens was the most abundant species of Thermocyclops from Nigeria [18]. However, the identification of different species of copepods solely based on morphology has technical limitations [19], as cryptic species are often detected. Therefore, more discerning methods such as DNA barcoding are needed to investigate copepod taxonomy, especially to recognize morphologically cryptic genetic lineages [20].

DNA barcoding has already been successfully applied to estimate the species/genetic diversity in many zooplankton taxa [21], as it can be used for rapid, accurate, reliable and remote identification of specimens of all metazoan [22]. A fragment of the mitochondrial cytochrome c oxidase subunit I (COI) gene has proved to be a useful marker for many biodiversity studies [23, 24], as COI has advantages of being effective for species identification from a wider range of metazoan phyla and possessing a phylogenetic signal which can be used over a wider range of taxonomic levels [22]. The COI marker has been successfully applied to species identification for cladocerans [25, 26]. For example, DNA barcoding was used to identify sibling cryptic species of the Ceriodaphnia cornuta species complex from Australia [27] and to examine 61 species of Cladocera, such as Daphnia, Diaphanosoma, Ceriodaphnia, Moina and Alona, from Mexico and Guatemala [28]. This approach has also been applied to copepods [29]. For instance, a study reported 800 new sequences of 63 marine copepod species by using a COI marker [30]. Elías-Gutierrez et al. [28] examined 21 species of Copepoda from Mexico and Guatemala by applying a COI marker. Using COI is highly advantageous because it can also detect cryptic species, a phenomenon that is very common in copepod assemblages [31]. For example, three genetically divergent but morphologically similar forms of Hemidiaptomus (Occidodiaptomus) ingens were detected throughout the distribution range of this species complex [32]. Moreover, Oithona similis s.l. was found to be a complex of nine cryptic species instead of a single cosmopolitan species, according to a COI and a nuclear ribosomal 28S genetic marker [20]. Similarly, the nominal species “Eudiaptomus hadzici” in the Western Balkans consists of four cryptic species according to a mitochondrial (COI) and a nuclear (nH3) marker [33]. DNA barcoding often reveals differences between allopatric populations. In that situation, it is difficult to decide whether this indicates different genetic lineages or simply geographical intraspecific variation. For instance, several different genetic lineages of Moina which were allopatric in a phylogeny were assigned to a single species, because they had similar morphology [34].

The phylogeny of copepods had been widely studied using molecular data. Recently, a comprehensive study from Asia showed a high species diversity of copepods in South Korea [29]. In that study, 133 sequenced individuals represented 94 species belonging to six different orders [29]. Another study has shown that Sinocalanus tenellus consists of two very distinct clades in China, suggesting they are parts of a complex of cryptic species [35]. Moreover, Karanovic [36] detected a new species of Schizopera from Japan, which was the first member of its genus reported in Japanese freshwater ecosystems, and it had no close relatives from elsewhere in the world. Another study has revised the higher systematics of copepods and proposed the new taxa Canuelloida ordo. nov., Smirnovipinidae fam. nov. and Cyclopicinidae fam. nov. [37]. Use of molecular data has not been restricted to species-level taxonomy [20], for example, the phylogeography of copepods has been also frequently investigated. They focused on the genealogical lineages of closely related species of copepods and their geographic distributions, by combining the information from phylogenetics, molecular genetics, population genetics, geology, paleontology, demography, ethology and historical biogeography [38]. For instance, two species of copepods (i.e. Neodiaptomus schmackeri and Mongolodiaptomus birulai) occur in Chinese Taiwan: there was little gene flow among populations for both species [39]. Additionally, four populations of Leptodiaptomus cf. sicilis  in Mexico were found to diverge into 3 distinct phenotypes, and their specialization was further supported by molecular data which showed persistence of a founder effect, limited gene flow, and a pattern of allopatric speciation [40].

There have been no studies on phylogeography and genetic diversity of copepods from Nigeria. In this study, we analyzed 15 copepod populations (out of 32 pools/lakes sampled) from Nigeria. By analyzing sequence variation in the COI gene, we aimed to explore the species diversity and distribution of copepods among these populations. Our expectation was to detect several members of the Cyclopidae, as it is commonly observed worldwide [29, 37]. We also investigated the phylogeography of Cyclopidae in Nigeria.


Species and COI genetic diversity

One to 9 specimens of Cyclopidae were sequenced per location, and a total of 88 Cyclopidae COI sequences were successfully obtained from 15 freshwater lakes around Southeast Nigeria, of which 18 unique haplotypes were detected (Tables 1 and 2). None of the COI sequences exhibited characteristics of nuclear pseudogenes (frame shifts or premature stop codons). Two independent species-delimitation methods (i.e. GMYC and bPTP) based on the COI Bayesian tree consistently identified 12 Cyclopidae species from Nigeria: 5 species of Tropocyclops (i.e. T. cf. confinis, T. cf. onabamiroi, T. cf. prasinus, T. cf. prasinus shagamiensis and T. cf. mellanbyi), 5 species of Mesocyclops (i.e. M. cf. aspericornis, M. cf. dussarti, M. cf. ogunnus, M. cf. aequatorialis similis and M. cf. salinus) and 2 species of Thermocyclops (i.e. T. decipiens and T. cf. crassus; Figs. 1 and 2). Species identified through molecular analyses fell into 2 distinct clades (i.e. clade I: T. cf. confinis, T. cf. onabamiroi, T. cf. prasinus, T. cf. prasinus shagamiensis and T. cf. mellanbyi; clade II: M. cf. aspericornis, M. cf. dussarti, M. cf. ogunnus, M. cf. aequatorialis similis, M. cf. salinus, T. decipiens and T. cf. crassus). Two Thermocyclops species were in the same clade as the Mesocyclops species. Pairwise genetic distances (uncorrected p-distances) based on COI sequence analysis ranged from 0.05 to 0.257 between species (Table 3). For each species, the population haplotype diversity (Hd) of COI ranged from 0 to 0.533, and the population nucleotide diversity (π) ranged from 0 to 6.86 × 10− 3 (Table 2).

Table 1 List of localities inhabited by Cyclopidae (name, abbreviation, geographical position), sampling time and water surface temperature
Table 2 Genetic characterization of sequenced individuals of Cyclopidae in each population
Fig. 1

a Geographic locations of sampling sites for Cyclopidae in Nigeria. b Haplotype network of Cyclopidae from Nigeria, based on the COI gene (544 bp). Each circle represents a unique haplotype and its size reflects the number of individuals expressing that haplotype. Color codes denote geographic location of populations. Portion of circles indicate distribution of haplotypes among different populations. The number of marks on connecting lines indicates the number of mutations between haplotypes. For lake abbreviations see Table 1. The map was obtained from ArcGIS and edited in Adobe Illustrator

Fig. 2

Bayesian phylogenetic tree and species delimitation results for Cyclopidae from Nigeria, based on the COI gene (547 bp). The IDs for shared haplotypes are provided in Table 2; for origin of reference sequence IDs see Table S1. Only posterior probabilities > 0.70 are shown. Species delimitation according to the GMYC and bPTP methods are indicated. For the bPTP method, the statistical support (PP) for species membership is also shown. Paracalanus parvus was used as an outgroup

Table 3 Uncorrected pairwise genetic distances (p-distances) among species of Cyclopidae copepods based on mitochondrial COI sequence analysis. The number in bracket denote the sample size of each species

Geographic distribution of species

Based on the haplotype network, seven out of 12 species detected through analysis of the COI gene occurred at more than one locality in Nigeria (Fig. 1b). The most frequently occurring species in this study was T. cf. prasinus, which had 4 haplotypes and was distributed in 4 lakes, including A5G, AOR, O3M and U1H, and one of the 4 haplotypes was shared by 3 lakes (A5G, AOR, O3M). Such a pattern was also observed in species M. cf. dussarti, which had 2 haplotypes and one of them was shared by 3 lakes (N1O, N2O and UBS) (Fig. 1b). Different Cyclopidae species co-existed in the same lake. For example, three species (i.e. T. cf. confinis, T. cf. prasinus and T. cf. crassus) co-existed in Lake A5G (Fig. 1b). Similarly, T. cf. onabamiroi, M. cf. aequatorialis similis and M. cf. salinus co-existed in Lake UII (Fig. 1b). Moreover, five out of 18 haplotypes were shared by different lakes (Fig. 1b). The most abundant haplotype was CTH1, including 21 specimens shared by A1G, A2G and IHE. This was followed by CMS1, shared by N1O, N2O and UBS, and CTR1, shared by A5G, AOR and O3M (Fig. 1b). Four species (i.e. T. cf. onabamiroi, T. cf. mellanbyi, M. cf. ogunnus and M. cf. aequatorialis similis) expressed only one haplotype with a single individual (Fig. 1b).


Through analysis of COI sequence variation, we explored the species diversity and distribution of copepods within the family Cyclopidae for 15 Nigerian freshwater ecosystems, the first such study for West Africa. Our results suggested a high species diversity of Cyclopidae copepods over a small geographic sampling range.

High species diversity has already been reported in the copepods from Nigeria [18, 19]. Forty species of Cyclopidae copepods from Nigeria were described based on morphological characteristics in the 1990s [19]. Here, we did not detect any new species based on molecular data; all the 12 species identified in the present study were described in [19]. In agreement with a previous study based on morphology [18], we found that T. decipiens was the most abundant species of Thermocyclops from Nigeria. However, M. aspericornis was recorded as the most abundant species of Mesocyclops in Nigerian waterbodies [18], whereas we found that M. cf. dussarti is the most abundant species of the genus Mesocyclops. This inconsistency might be explained by the relatively small sampling region in our study in Southeast Nigeria.

Globally, high levels of species diversity of copepods have also been reported [29, 41]. For example, 53 Caligus species were present in Chinese Taiwan, and many more species remained to be discovered from this region [5]. Similarly, thirteen species of Copepoda, including three members of Calanoida (Diaptomidae) and ten members of Cyclopoida (Eucyclopinae and Cyclopinae), were recorded in Chiapas, Mexico [42]. Indeed, high species diversity, even in a relatively small area, has often been observed in copepods [41, 43, 44]. For example, a study identified 43 species that belonged to 11 genera of copepods in Sagami Bay [43]. Another study identified 48 species of copepods in Tolo Harbour, Hong Kong, and Oithona rigida, O. simplex and Paracalanus crassirostris were found to be the most abundant species [44]. Here, we detected 12 species with several species and species complexes across three genera in Nigeria; suggesting a high species diversity of Cyclopidae in Southeast Nigeria.

In agreement with a previous study of Cyclopoida in Nigeria [18], our results showed that the same species could be found in geographically separate populations, which also suggests that there are not extensive and common cryptic species in these sampled lakes. Thermocyclops decipiens has also been detected in Antilles, Central America, Columbia, Venezuela, east of the Andes, Brazil [45] and Congo [46], indicating that this species has a wide distribution across continents. In contrast to several copepod species with high genetic divergence over short distances, e.g. Tigriopus californicus [47], our data showed genetic similarity of the T. decipiens populations from different continents. A similar phenomenon has been detected in some open-ocean copepods which have more obvious dispersal routes. For example, it was found that several mtCOI haplotypes of Calanus pacificus were distributed across multiple sampling location from the North Pacific Ocean [48]. By using restriction site-associated DNA sequencing, no significant genetic differentiation was found among Centropages typicus samples collected from different NW Atlantic regions with clear connectivity [49]. Zooplankton species often have vast ranges [34, 50]. For example, Daphnia galeata has been detected in both China and Europe with some haplotypes shared across large distances [51]. Birds are often regarded as the key vectors for the dispersal of resting eggs of aquatic zooplankton [52], across geographical barriers.

We found that different sibling species of Cyclopidae co-existed in the same Nigerian lake, a common finding in copepods [29]. For example, a study of the genus Mesocyclops conducted in Africa reported that M. major and M. ogunnus often co-existed in the same waterbody [53]. Similarly, another study from Nigeria reported that it was common for up to 3 Mesocyclops species co-existed in a single locality [18]. Sympatry provides a possibility for interspecific hybridization, which is believed to be a common phenomenon in zooplankton [54]. Hybridization has often been observed in copepods. For example, hybrids between Calanus glacialis and C. finmarchicus were detected along the Atlantic and Arctic Canadian coast [55]. Another study also found that hybridization occurred between a female Neocalanus cristatus and a male N. plumchrus, and was then followed by backcrossing to a N. plumchrus individual [56]. Our molecular data indicated paraphyly between Thermocyclops and Mesocyclops, which might reflect introgression resulting from hybridization [57]. Paraphyly has also been observed in other zooplankton, for example, in the Daphnia pulex species complex [58] and Moina [34]. Future studies and nuclear markers are needed to investigate gene introgression among the copepod species from Nigeria.

Here, Cyclopidae in Nigeria showed a high species diversity, but for each species, the haplotype diversity and nucleotide diversity were rather low. Consistently, low haplotype and nucleotide diversities were observed in the copepods Calanus finmarchicus [59], C. agulhensis and C. sinicus [60]. Patterns of population genetic diversity could be caused by different evolutionary forces, such as mutation, migration, genetic drift and natural selection. How these evolutionary forces affect the population genetic diversity depends on many factors, including species’ response to environmental changes, and the past and present sizes of the population [59]. Another explanation might be that the structuring of a metapopulation together with founder effects resulted in low population genetic diversity [61, 62]. However, the limited sampling of each species among populations cannot be ruled out as a cause for their low genetic diversity in this study.


In conclusion, our data revealed a high species diversity of Cyclopidae in Southeast Nigeria: twelve species were detected. Our geographical sampling scale in this study was quite small, and therefore, further studies are called for a comprehensive understanding of species distribution and genetic diversity of Cyclopidae in West Africa.



Copepod specimens were collected from 15 freshwater lakes around Southeast Nigeria (Fig. 1a and Table 1). Samples were collected using a 125-μm plankton net hauled vertically through the water column at three different sites per location. Samples collected from different sites in the same location were pooled together and preserved in 95% ethanol. All specimens were identified morphologically according to the morphological description of copepods in Nigeria [17,18,19, 63], which also worked as taxonomic keys in this paper.

DNA extraction, PCR amplification and sequencing

Copepods were processed for molecular analyses (Table 2). The cephalosome portion of the prosome was obtained from each individual copepod to avoid DNA contamination from prey items in the gut, using a microscopic tweezer and a sharp blade under the stereomicroscope. Total genomic DNA was extracted from the head using H3 buffer with proteinase K (30 μL), containing 10 mM Tris-HCl, 0.05 M KCl, 0.005% Tween 20, 0.005% NP-40 and 10 mg/ml proteinase K (MERCK, Germany). Samples were incubated overnight at 55 °C in a water-bath with mild shaking. The proteinase K was irreversibly denatured after a 12-min incubation at 95 °C. The homogenate was centrifuged briefly and stored at 4 °C before use. A 680 base-pair fragment of the COI gene was amplified using the consensus primer pair (forward: 5′- GGT CAA CAA ATC ATA AAG ATA TTG G − 3′; reverse: 5′- TAA ACT TCA GGG TGA CCA AAA AAT CA -3′ [64];). Polymerase chain reaction (PCR) was carried out in a total volume of 20 μL, consisting of 2 μL 10 X PCR buffer (10 mM Tris-HCl, pH 8.3, 5 mM MgCl2, 50 mM KCl), 2 μL 2.5 mM of each dNTP, 1 μL 0.5 μM of each primer, 11.6 μL water, 2 units of Taq DNA polymerase (SuperTherm DNA polymerase, Taq HS from TAKARA BIO INC., California, USA) and 2 μL of genomic DNA. The PCR thermocycle protocol was as follows: denaturation at 94 °C for 1 min, then 40 cycles of 1 min at 94 °C, 1.5 min at 40 °C and 1.5 min at 72 °C; followed by a final elongation at 72 °C for 6 min. The success of amplification was checked using agarose gel electrophoresis. Afterwards, the PCR products were purified (High Pure PCR Product Purification Kit, Roche Diagnostics) and sequenced in the forward direction on an ABI PRISM 3730 DNA capillary sequencer by Invitrogen Trading Co., Ltd. (China). Chromatograms were checked for ambiguous base calling and errors in base calling were corrected using MEGA 6 [65], and the Phred quality scores of the sequences were examined with Chromas Lite Version 2.1 (Technelysium Pty. Ltd., South Brisbane, Australia). Sequences with double peaks or noise were re-sequenced in the reverse direction, and only chromatograms of high quality were applied to the following genetic analyses. All newly obtained sequences were submitted to GenBank under accession numbers MN990471-MN990488.

Sequence alignment and genetic diversity

We identified unique haplotypes in DnaSP 5.10 [66]. MUSCLE [67] implemented in MEGA 6 was used to align the sequences that were subsequently translated into amino acids to examine the presence of stop codons. Afterwards, all haplotypes were aligned, together with 8 reference sequences obtained from GenBank (Table S1), using the Clustal W algorithm [68] in MEGA 6. Then, all the sequences were timed to a uniform length of 544 bp in MEGA 6. For each species, the number of haplotypes (N2), haplotype diversity (Hd) and nucleotide diversity (π) per population (populations with sample size less than 3 were excluded) were calculated in DnaSP 5.10.

Phylogenetic analyses

The test of Xia et al. [69] implemented in DAMBE 5 [70] was used to inspect potential loss of phylogenetic signal resulting from substitution saturation at the COI locus. A phylogenetic tree was then constructed using the Bayesian method in BEAST 2 [71], with a tree sampled every 1000 generations among 10,000,000, a burn-in of 25%, and the final 10,000 trees summarized using TreeAnnotator. The best-fitting substitution model was GTR + G + I according to the corrected Akaike Information Criterion in jModeltest v. 2.1.7 [72]. We applied a strict molecular clock and set other tree priors to their default values. Tracer v1.6 [73] was applied to ensure that enough generations were computed. Paracalanus parvus, a member of the Calanoida phylogenetically close to Cyclopoida, was used as an outgroup.

Species identification and phylogeographic analyses

To test the hypothesis that the Cyclopidae in Nigeria contains high biodiversity, two independent species delimitation methods were applied: the general mixed Yule coalescent model (GMYC [74]) and Poisson tree processes methods (PTP [75]). The GMYC model is a likelihood-based method using an ultrametric tree to delimit species by fitting within- and between-species branching models to reconstruct gene trees. We performed the GMYC modeling using the “splits” package [76] in R 2.15 [77] and conducted the PTP calculations on the bPTP webserver (, with 100,000 Markov Chain Monte Carlo (MCMC) generations, thinning set to 100, burnin at 25% and a Bayesian search performed. The input phylogenetic tree was generated using BEAST 2 (see above). A network of COI haplotypes was then constructed to visualize genetic relationships among populations using Haploviewer [78]. The maximum likelihood tree inferred with MEGA 6 using the best model GTR + G + I (by jModeltest v. 2.1.7) was applied as input. Uncorrected pairwise genetic distances between species were calculated in MEGA 6 based on COI.

Availability of data and materials

All the sequencing data are available via NCBI (under accession numbers MN990471-MN990488).



Mitochondrial cytochrome c oxidase subunit I gene


Polymerase chain reaction


Base pairs


The general mixed Yule coalescent model


Poisson tree processes methods


  1. 1.

    Boxshall GA, Halsey SH. An introduction to copepod diversity. London: Ray Society; 2004.

  2. 2.

    Huys R, Boxshall GA. Copepod Evolution. London: The Ray Society; 1991.

  3. 3.

    Roemmich D, Mcgowan J. Climatic warming and the decline of zooplankton in the California current. Science. 1995;267(5202):1324-6.

  4. 4.

    Schminke H. Entomology for the copepodologist. J Plankton Res. 2007;29:i149–62.

    Article  Google Scholar 

  5. 5.

    Ho J, Lin C, Liu W. High diversity of Caligus species (Copepoda: Siphonostomatoida: Caligidae) in Taiwanese waters. Zootaxa. 2016;4174:114–21.

    PubMed  Article  Google Scholar 

  6. 6.

    Humes AG. How many copepods. Hydrobiologia. 1994;293:1–7.

    Article  Google Scholar 

  7. 7.

    Chivers WJ, Walne AW, Hays GC. Mismatch between marine plankton range movements and the velocity of climate change. Nat Commun. 2017;8:14434.

  8. 8.

    Beaugrand G, Reid PC, Ibanez F, Lindley A, Edwards M. Reorganization of North Atlantic marine copepod biodiversity and climate. Science. 2002;296:1692–4.

  9. 9.

    Costanza R, d'Arge R, de Groot R, Farber S, Grasso M, Hannon B, Limburg K, Naeem S, O'Neill R, Paruelo J, et al. The value of the world's ecosystem services and natural capital. Nature. 1997;387:253–60.

    Article  CAS  Google Scholar 

  10. 10.

    Pond D, Tarling G. Phase transitions of wax esters adjust buoyancy in diapausing Calanoides acutus. Limnol Oceanogr. 2011;56:1310–8.

    Article  CAS  Google Scholar 

  11. 11.

    Hays GC, Richardson AJ, Robinson C. Climate change and marine plankton. Trends Ecol Evol. 2005;20(6):337–44.

    PubMed  Article  Google Scholar 

  12. 12.

    Blanco-Bercial L, Bradford-Grieve J, Bucklin A. Molecular phylogeny of the Calanoida (Crustacea: Copepoda). Mol Phylogenet Evol. 2011;59(1):103–13.

    PubMed  Article  Google Scholar 

  13. 13.

    Bradford-Grieve JM, Boxshall GA, Ahyong ST, Ohtsuka S. Cladistic analysis of the calanoid Copepoda. Invertebr Syst. 2010;24(3):291–321.

    Article  Google Scholar 

  14. 14.

    Okoye S, Onwuliri C, Anosike J. A survey of predilection sites and degree of disability associated with guineaworm (Dracunculus medinensis). Int J Parasitol. 1995;25:1127–9.

    PubMed  Article  CAS  Google Scholar 

  15. 15.

    Anosike J, Azoro V, Nwoke B, Keke R, Okere A, Oku E, Ogbulie J, Tony-Njoku R, Okoro O, Nwosu D. Dracunculiasis in the north eastern border of Ebonyi state, south eastern Nigeria. Int J Hyg Envir Heal. 2003;206:45–51.

    Article  Google Scholar 

  16. 16.

    Morenikeji O, Asiatu A. Progress in dracunculiasis eradication in Oyo state, south-West Nigeria: a case study. Afr Health Sci. 2010;10:297–301.

    PubMed  PubMed Central  CAS  Google Scholar 

  17. 17.

    Onabamiro S. Four new species of Cyclops sensu lat. (Crustacea: Copepoda) from Nigeria. Proc Zool Soc Lond. 1952;122:253–66.

    Article  Google Scholar 

  18. 18.

    Jeje C. A revision of the Nigerian species of the genera Mesocyclops Sars, 1914 and Thermocyclops Kiefer, 1927 (Copepoda: Cyclopoida). Hydrobiologia. 1988;164:171–84.

    Article  Google Scholar 

  19. 19.

    Boxshall GA, Braide E. The freshwater cyclopoid copepods of Nigeria, with an illustrated key to all species. Bull Br Mus Nat Hist Zool. 1991;57:185–212.

    Google Scholar 

  20. 20.

    Cornils A, Wend-Heckmann B, Held C. Global phylogeography of Oithona similis s.l. (Crustacea, Copepoda, Oithonidae) - a cosmopolitan plankton species or a complex of cryptic lineages? Mol Phylogenet Evol. 2017;107:473–85.

    PubMed  Article  Google Scholar 

  21. 21.

    Hebert PDN, Ratnasingham S, deWaard JR. Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. P Roy Soc B-Biol Sci. 2003;270:S96–9.

  22. 22.

    Bucklin A, Steinke D, Blanco-Bercial L. DNA barcoding of marine Metazoa. Annu Rev Mar Sci. 2011;3:471–508.

  23. 23.

    Costa F, deWaard J, Boutillier J, Ratnasingham S, Dooh R, Hajibabaei M, Hebert PDN. Biological identifications through DNA barcodes: the case of the Crustacea. Can J Fish Aquat Sci. 2007;64:272–95.

  24. 24.

    Hebert PDN, Cywinska A, Ball S, deWaard JR. Biological identifications through DNA barcode. P Roy Soc B-Biol Sci. 2003;270:313–21.

    Article  CAS  Google Scholar 

  25. 25.

    deWaard JR, Sacherova V, Cristescu ME, Remigio EA, Crease T, Hebert PDN. Probing the relationships of the branchiopod crustaceans. Mol Phylogenet Evol. 2006;39:491–502.

    PubMed  Article  CAS  Google Scholar 

  26. 26.

    Jeffery NW, Elias-Gutierrez M, Adamowicz SJ. Species diversity and phylogeographical affinities of the Branchiopoda (Crustacea) of Churchill, Manitoba, Canada. PLOS One. 2011;6(5):e18364.

  27. 27.

    Sharma P, Kotov AA. Molecular approach to identify sibling species of the Ceriodaphnia cornuta complex (Cladocera: Daphniidae) from Australia with notes on the continental endemism of this group. Zootaxa. 2013;3702:79–89.

    PubMed  Article  Google Scholar 

  28. 28.

    Elias-Gutierrez M, Jeronimo FM, Ivanova NV, Valdez-Moreno M, Hebert PDN. DNA barcodes for Cladocera and Copepoda from Mexico and Guatemala, highlights and new discoveries. Zootaxa. 2008;1839:1–42.

  29. 29.

    Baek SY, Jang KH, Choi EH, Ryu SH, Kim SK, Lee JH, Lim YJ, Lee J, Jun J, Kwak M, et al. DNA barcoding of metazoan zooplankton Copepods from South Korea. PLOS One. 2016;11(7):e0157307.

  30. 30.

    Blanco-Bercial L, Cornils A, Copley N, Bucklin A. DNA barcoding of marine copepods: assessment of analytical approaches to species identification. PLOS Curr. 2014.

  31. 31.

    Dippenaar S, Mathibela R, Bloomer P. Cytochrome oxidase I sequences reveal possible cryptic diversity in the cosmopolitan symbiotic copepod Nesippus orientalis Heller, 1868 (Pandaridae: Siphonostomatoida) on elasmobranch hosts from the KwaZulu-Natal coast of South Africa. Exp Parasitol. 2010;125:42–50.

    PubMed  Article  CAS  Google Scholar 

  32. 32.

    Marrone F, Lo Brutto S, Hundsdoerfer AK, Arculeo M. Overlooked cryptic endemism in copepods: systematics and natural history of the calanoid subgenus Occidodiaptomus Borutzky 1991 (Copepoda, Calanoida, Diaptomidae). Mol Phylogenet Evol. 2013;66(1):190–202.

    PubMed  Article  Google Scholar 

  33. 33.

    Previsic A, Gelemanovic A, Urbanic G, Ternjej I. Cryptic diversity in the Western Balkan endemic copepod: four species in one? Mol Phylogenet Evol. 2016;100:124–34.

    PubMed  Article  Google Scholar 

  34. 34.

    Ni YJ, Ma XL, Hu W, Blair D, Yin MB. New lineages and old species: lineage diversity and regional distribution of Moina (Crustacea: Cladocera) in China. Mol Phylogenet Evol. 2019;134:87–98.

    PubMed  Article  Google Scholar 

  35. 35.

    Wang XY, Ma XL, Hu W, Yin MB. Genetic diversity and population differentiation of the freshwater copepod Sinocalanus tenellus (Calanoida, Centropagidae) in China. J Limnol. 2018;77(2):300–7.

    Article  Google Scholar 

  36. 36.

    Karanovic T, Kim K, Grygier M. A new species of Schizopera (Copepoda: Harpacticoida) from Japan, its phylogeny based on the mtCOI gene and comments on the genus Schizoperopsis. J Nat Hist. 2015;49:2493–526.

    Article  Google Scholar 

  37. 37.

    Khodami S, McArthur J, Blanco-Bercial L, Arbizu P. Molecular phylogeny and revision of copepod orders (Crustacea: Copepoda). Sci Rep. 2017;7:9164.

  38. 38.

    Avise J. Phylogeography: The History and Formation of Species. Cambridge: Harvard University Press; 2000.

  39. 39.

    Young S, Lin S, Liu M. Genetic diversity and population structure of two freshwater copepods (Copepoda: Diaptomidae), Neodiaptomus schmackeri (Poppe and Richard, 1892) and Mongolodiaptomus birulai (Rylov, 1922) from Taiwan. Diversity. 2013;5:796–810.

    Article  Google Scholar 

  40. 40.

    Barrera-Moreno O, Ciros-Perez J, Ortega-Mayagoitia E, Alcantara-Rodriguez J, Piedra-Ibarra E. From local adaptation to ecological speciation in copepod populations from neighboring lakes. PLoS One. 2015;10:e0125524.

  41. 41.

    Brancelj A. Microdistribution and high diversity of Copepoda (Crustacea) in a small cave in Central Slovenia. Hydrobiologia. 2002;477:59–72.

    Article  Google Scholar 

  42. 42.

    Gutierrez-Aguirre M, Cervantes-Mrtinez A. Diversity of freshwater copepods (Maxillopoda: Copepoda: Calanoida, Cyclopoida) from Chiapas, Mexico with a description of Mastigodiaptomus suarezmoralesi sp. nov. J Nat Hist. 2013;47:479–98.

  43. 43.

    Nishida S, Kuriyama M. Species diversity and niche-partitioning in the pelagic copepods of the family Scolecitrichidae (Calanoida). Crustaceana. 2006;79:293–317.

  44. 44.

    Wong C, Chan A, Chen Q. Planktonic copepods of Tolo harbour, Hong Kong. Crustaceana. 1993;64:76–84.

    Article  Google Scholar 

  45. 45.

    Landa G, Barbosa F, Rietzler A, Barbosa P. Thermocyclops decipiens (Kiefer, 1929) (Copepoda, Cyclopoida) as indicator of water quality in the state of Minas Gerais, Brazil. Braz Arch Biol Techn. 2007;50:695–705.

  46. 46.

    Mirabdullayev I, Reid J, Ueda H. Copepoda: Cyclopoida genera Mesocyclops and Thermocyclops. Genus Thermocyclops Kiefer, 1927. In: Guides to the identification of the microinvertebrates of the continental waters of the world. Leiden: Backhuys Publishers; 2003. p. 214–302.

  47. 47.

    Willett C, Ladner JT. Investigations of fine-scale phylogeography in Tigriopus californicus reveal historical patterns of population divergence. BMC Evol Biol. 2009;9:139.

  48. 48.

    Nuwer M, Frost B, Armbrust V. Population structure of the planktonic copepod Calanus pacificus in the North Pacific Ocean. Mar Biol. 2008;156:107–15.

    Article  Google Scholar 

  49. 49.

    Blanco-Bercial L, Bucklin A. New view of population genetics of zooplankton: RADseq analysis reveals population structure of the North Atlantic planktonic copepod Centropages typicus. Mol Ecol. 2016;25:1566–80.

    PubMed  Article  CAS  Google Scholar 

  50. 50.

    Ma X, Petrusek A, Wolinska J, Hu W, Yin M. Lineage diversity and reproductive modes of the Daphnia pulex group in Chinese lakes and reservoirs. Mol Phylogenet Evol. 2019;130:424–33.

    PubMed  Article  Google Scholar 

  51. 51.

    Yin M, Wang X, Ma X, Giessler S, Petrusek A, Griebel J, Hu W, Wolinska J. Cytonuclear diversity and shared mitochondrial haplotypes among Daphnia galeata populations separated by seven thousand kilometres. BMC Evol Biol. 2018;18:130.

    PubMed  PubMed Central  Article  Google Scholar 

  52. 52.

    Green AJ, Figuerola J. Recent advances in the study of long-distance dispersal of aquatic invertebrates via birds. Divers Distrib. 2005;11(2):149–56.

    Article  Google Scholar 

  53. 53.

    Van de Velde I. Revision of the African species of the genus Mesocyclops Sars, 1914 (Copepoda: Cyclopidae). Hydrobiologia. 1984;109:3–66.

    Article  Google Scholar 

  54. 54.

    Hebert PDN. Interspecific hybridization between cyclic parthenogens. Evolution. 1985;39:216–20.

    PubMed  Article  Google Scholar 

  55. 55.

    Parent GJ, Plourde S, Joly P, Turgeon J. Phenology and fitness of Calanus glacialis, C. finmarchicus (Copepoda), and their hybrids in the St. Lawrence estuary. Mar Ecol Prog Ser. 2015;524:1–9.

    Article  Google Scholar 

  56. 56.

    Machida RJ, Lin Y. Occurrence of mitochondrial CO1 pseudogenes in Neocalanus plumchrus (Crustacea: Copepoda): hybridization indicated by recombined nuclear mitochondrial pseudogenes. PLoS One. 2017;12:e0172710.

  57. 57.

    Funk DJ, Omland KE. Species-level paraphyly and polyphyly: frequency, causes, and consequences, with insights from animal mitochondrial DNA. Annu Rev Ecol Evol S. 2003;34:397–423.

    Article  Google Scholar 

  58. 58.

    Colbourne JK, Crease TJ, Weider LJ, Hebert PDN, Dufresne F, Hobaek A. Phylogenetics and evolution of a circumarctic species complex (Cladocera : Daphnia pulex). Biol J Linn Soc. 1998;65(3):347–65.

    Google Scholar 

  59. 59.

    Bucklin A, Wiebe PH. Low mitochondrial diversity and small effective population sizes of the copepods Calanus finmarchicus and Nannocalanus minor: possible impact of climatic variation during recent glaciation. J Hered. 1998;89(5):383–92.

    PubMed  Article  CAS  Google Scholar 

  60. 60.

    Kozol R, Blanco-Bercial L, Bucklin A. Multi-gene analysis reveals a lack of genetic divergence between Calanus agulhensis and C. sinicus (Copepoda; Calanoida). PLOS One. 2012;7:e45710.

  61. 61.

    Hanski I, Gaggiotti O. Ecology, genetics, and evolution of metapopulations. Amsterdam: Elsevier Academic Press; 2004.

  62. 62.

    Pannell J, Charlesworth B. Neutral genetic diversity in a metapopulation with recurrent local extinction and recolonization. Evolution. 1999;53:664–76.

    PubMed  Article  Google Scholar 

  63. 63.

    Onabamiro S. Some new species of Cyclops sensu lat (Crustacea: Copepoda) from Nigeria. Zool J Linn Soc-Lond. 1957;43:123–33.

  64. 64.

    Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3(5):294–9.

  65. 65.

    Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  66. 66.

    Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25(11):1451–2.

    PubMed  Article  CAS  Google Scholar 

  67. 67.

    Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  68. 68.

    Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22(22):4673–80.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  69. 69.

    Xia XH, Xie Z, Salemi M, Chen L, Wang Y. An index of substitution saturation and its application. Mol Phylogenet Evol. 2003;26(1):1–7.

    PubMed  Article  CAS  Google Scholar 

  70. 70.

    Xia XH. DAMBE5: a comprehensive software package for data analysis in molecular biology and evolution. Mol Biol Evol. 2013;30(7):1720–8.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  71. 71.

    Bouckaert R, Heled J, Kuhnert D, Vaughan T, Wu CH, Xie D, Suchard MA, Rambaut A, Drummond AJ. BEAST 2: a software platform for Bayesian evolutionary analysis. PLOS Comput Biol. 2014;10(4):e1003537.

  72. 72.

    Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9(8):772.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  73. 73.

    Rambaut A, Drummond A, Xie D, Baele G, Suchard M. Posterior summarization in Bayesian phylogenetics using Tracer 1.7. Syst Biol. 2018;67(5):901-4.

  74. 74.

    Pons J, Barraclough TG, Gomez-Zurita J, Cardoso A, Duran DP, Hazell S, Kamoun S, Sumlin WD, Vogler AP. Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Syst Biol. 2006;55(4):595–609.

    PubMed  Article  Google Scholar 

  75. 75.

    Zhang JJ, Kapli P, Pavlidis P, Stamatakis A. A general species delimitation method with applications to phylogenetic placements. Bioinformatics. 2013;29(22):2869–76.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  76. 76.

    Ezard T, Fujisawa T, Barraclough T. Splits: SPecies’ LImits by threshold statistics R package; 2009.

    Google Scholar 

  77. 77.

    R Development Core Team: R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. 2009.

    Google Scholar 

  78. 78.

    Walter S, Ewing G, von Haeseler A. The performance of phylogenetic algorithms in estimating haplotype genealogies with migration. Mol Ecol. 2011;20:1952–63.

Download references


We thank David Blair for critical discussion and the linguistic help, and the anonymous reviewer for useful comments on the earlier version of this manuscript.


This research was funded by the National Natural Science Foundation of China (31670380) and Natural Science Foundation of Shanghai (16ZR1402900) to MY. The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information




MY designed the study, CE, EO, CO, FA, CN, JO and PU collected samples, YN and JW carried out the molecular work. YN, WH and MY contributed to data analyses. MY wrote the manuscript with the help of YN. All authors read and approved the final version.

Corresponding author

Correspondence to Mingbo Yin.

Ethics declarations

Ethics approval and consent to participate

Collection of zooplankton (copepods in this study) did not require specific permissions because these samples were obtained from unprotected lakes that are open for public activities. Our study did not involve the use or collection of endangered or protected species.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1:

Table S1. List of reference COI sequences of Cyclopidae (from South Korea, Brazil and China) and the outgroup used in this study.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ni, Y., Ebido, C.C., Odii, E.C. et al. Phylogeography and genetic diversity of the copepod family Cyclopidae (Crustacea: Cyclopoida) from freshwater ecosystems of Southeast Nigeria. BMC Evol Biol 20, 45 (2020).

Download citation


  • Cyclopidae
  • Species diversity
  • COI
  • Nigeria