key: cord-010149-zupx8qn1
authors: Ni, Yijun; Ebido, Chike Chukwuenyem; Odii, Elijah Chibueze; Wang, Jinhui; Orakwelu, Chinemerem Hodges; Abonyi, Francis Chukwuemeka; Ngene, Chinedu Innocent; Okoro, Joseph Onyekwere; Ubachukwu, Patience Obiageli; Hu, Wei; Yin, Mingbo
title: Phylogeography and genetic diversity of the copepod family Cyclopidae (Crustacea: Cyclopoida) from freshwater ecosystems of Southeast Nigeria
date: 2020-04-21
journal: BMC Evol Biol
DOI: 10.1186/s12862-020-01608-5
sha: 
doc_id: 10149
cord_uid: zupx8qn1

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

CO 2 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.

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 iden- (Table 3) . For each species, the population haplotype diversity (H d ) of COI ranged from 0 to 0.533, and the population nucleotide diversity (π) ranged from 0 to 6.86 × 10 − 3 ( Table 2) .

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 (See figure on previous page.) 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 it was common for up to 3 Mesocyclops species coexisted 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.

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 MgCl 2 , 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.

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 (N 2 ), haplotype diversity (H d ) and nucleotide diversity (π) per population (populations with sample size less than 3 were excluded) were calculated in DnaSP 5.10.

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.

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 likelihoodbased 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 (http://species.h-its.org/ptp/), 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.

Supplementary information accompanies this paper at https://doi.org/10. 1186/s12862-020-01608-5.

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.

Abbreviations COI: Mitochondrial cytochrome c oxidase subunit I gene; PCR: Polymerase chain reaction; bp: Base pairs; GMYC: The general mixed Yule coalescent model; PTP: Poisson tree processes methods

An introduction to copepod diversity. London: Ray Society

Copepod Evolution. London: The Ray Society

Climatic warming and the decline of zooplankton in the California current

Entomology for the copepodologist

High diversity of Caligus species (Copepoda: Siphonostomatoida: Caligidae) in Taiwanese waters

How many copepods

Mismatch between marine plankton range movements and the velocity of climate change

Reorganization of North Atlantic marine copepod biodiversity and climate

The value of the world's ecosystem services and natural capital

Phase transitions of wax esters adjust buoyancy in diapausing Calanoides acutus

Climate change and marine plankton

Molecular phylogeny of the Calanoida (Crustacea: Copepoda)

Cladistic analysis of the calanoid Copepoda

A survey of predilection sites and degree of disability associated with guineaworm (Dracunculus medinensis)

Dracunculiasis in the north eastern border of Ebonyi state, south eastern Nigeria

Progress in dracunculiasis eradication in Oyo state, south-West Nigeria: a case study

Four new species of Cyclops sensu lat. (Crustacea: Copepoda) from Nigeria

A revision of the Nigerian species of the genera Mesocyclops Sars

The freshwater cyclopoid copepods of Nigeria, with an illustrated key to all species

Global phylogeography of Oithona similis s.l. (Crustacea, Copepoda, Oithonidae) -a cosmopolitan plankton species or a complex of cryptic lineages?

Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species

Biological identifications through DNA barcodes: the case of the Crustacea

Biological identifications through DNA barcode

Probing the relationships of the branchiopod crustaceans

Species diversity and phylogeographical affinities of the Branchiopoda (Crustacea) of Churchill

Molecular approach to identify sibling species of the Ceriodaphnia cornuta complex (Cladocera: Daphniidae) from Australia with notes on the continental endemism of this group

DNA barcodes for Cladocera and Copepoda from Mexico and Guatemala, highlights and new discoveries

DNA barcoding of metazoan zooplankton Copepods from South Korea

DNA barcoding of marine copepods: assessment of analytical approaches to species identification

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

Overlooked cryptic endemism in copepods: systematics and natural history of the calanoid subgenus Occidodiaptomus Borutzky 1991 (Copepoda, Calanoida, Diaptomidae)

Cryptic diversity in the Western Balkan endemic copepod: four species in one?

New lineages and old species: lineage diversity and regional distribution of Moina (Crustacea: Cladocera) in China

Genetic diversity and population differentiation of the freshwater copepod Sinocalanus tenellus (Calanoida, Centropagidae) in China

A new species of Schizopera (Copepoda: Harpacticoida) from Japan, its phylogeny based on the mtCOI gene and comments on the genus Schizoperopsis

Molecular phylogeny and revision of copepod orders

The History and Formation of Species

Genetic diversity and population structure of two freshwater copepods (Copepoda: Diaptomidae), Neodiaptomus schmackeri (Poppe and Richard, 1892) and Mongolodiaptomus birulai (Rylov, 1922) from Taiwan

From local adaptation to ecological speciation in copepod populations from neighboring lakes

Microdistribution and high diversity of Copepoda (Crustacea) in a small cave in Central Slovenia

Diversity of freshwater copepods (Maxillopoda: Copepoda: Calanoida, Cyclopoida) from Chiapas, Mexico with a description of Mastigodiaptomus suarezmoralesi sp. nov

Species diversity and niche-partitioning in the pelagic copepods of the family Scolecitrichidae (Calanoida)

Planktonic copepods of Tolo harbour, Hong Kong

Thermocyclops decipiens (Kiefer, 1929) (Copepoda, Cyclopoida) as indicator of water quality in the state of Minas Gerais

Cyclopoida genera Mesocyclops and Thermocyclops

Investigations of fine-scale phylogeography in Tigriopus californicus reveal historical patterns of population divergence

Population structure of the planktonic copepod Calanus pacificus in the North Pacific Ocean

New view of population genetics of zooplankton: RADseq analysis reveals population structure of the North Atlantic planktonic copepod Centropages typicus

Lineage diversity and reproductive modes of the Daphnia pulex group in Chinese lakes and reservoirs

Cytonuclear diversity and shared mitochondrial haplotypes among Daphnia galeata populations separated by seven thousand kilometres

Recent advances in the study of long-distance dispersal of aquatic invertebrates via birds

Revision of the African species of the genus Mesocyclops Sars, 1914 (Copepoda: Cyclopidae)

Interspecific hybridization between cyclic parthenogens

Phenology and fitness of Calanus glacialis, C. finmarchicus (Copepoda), and their hybrids in the St. Lawrence estuary

Occurrence of mitochondrial CO1 pseudogenes in Neocalanus plumchrus (Crustacea: Copepoda): hybridization indicated by recombined nuclear mitochondrial pseudogenes

Species-level paraphyly and polyphyly: frequency, causes, and consequences, with insights from animal mitochondrial DNA

Phylogenetics and evolution of a circumarctic species complex (Cladocera : Daphnia pulex)

Low mitochondrial diversity and small effective population sizes of the copepods Calanus finmarchicus and Nannocalanus minor: possible impact of climatic variation during recent glaciation

Multi-gene analysis reveals a lack of genetic divergence between Calanus agulhensis and C. sinicus (Copepoda

Ecology, genetics, and evolution of metapopulations

Neutral genetic diversity in a metapopulation with recurrent local extinction and recolonization

Some new species of Cyclops sensu lat (Crustacea: Copepoda) from Nigeria

DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates

MEGA6: molecular evolutionary genetics analysis version 6.0

DnaSP v5: a software for comprehensive analysis of DNA polymorphism data

MUSCLE: multiple sequence alignment with high accuracy and high throughput

CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice

An index of substitution saturation and its application

DAMBE5: a comprehensive software package for data analysis in molecular biology and evolution

BEAST 2: a software platform for Bayesian evolutionary analysis

jModelTest 2: more models, new heuristics and parallel computing

Posterior summarization in Bayesian phylogenetics using Tracer 1.7

Sequence-based species delimitation for the DNA taxonomy of undescribed insects

A general species delimitation method with applications to phylogenetic placements

Splits: SPecies' LImits by threshold statistics R package

R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing

The performance of phylogenetic algorithms in estimating haplotype genealogies with migration

Publisher's Note

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

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.

Authors' contributions 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.

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.

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

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.

Not applicable.

The authors declare that they have no competing interests.