Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Divergent clades or cryptic species? Mito-nuclear discordance in a Daphnia species complex

BMC Evolutionary Biology volume 17, Article number: 227 (2017) Cite this article

Abstract

Background

Genetically divergent cryptic species are frequently detected by molecular methods. These discoveries are often a byproduct of molecular barcoding studies in which fragments of a selected marker are used for species identification. Highly divergent mitochondrial lineages and putative cryptic species are even detected in intensively studied animal taxa, such as the crustacean genus Daphnia. Recently, eleven such lineages, exhibiting genetic distances comparable to levels observed among well-defined species, were recorded in the D. longispina species complex, a group that contains several key taxa of freshwater ecosystems. We tested if three of those lineages represent indeed distinct species, by analyzing patterns of variation of ten nuclear microsatellite markers in six populations.

Results

We observed a discordant pattern between mitochondrial and nuclear DNA, as all individuals carrying one of the divergent mitochondrial lineages grouped at the nuclear level with widespread, well-recognized species coexisting at the same localities (Daphnia galeata, D. longispina, and D. cucullata).

Conclusions

A likely explanation for this pattern is the introgression of the mitochondrial genome of undescribed taxa into the common species, either in the distant past or after long-distance dispersal. The occurrence of highly divergent but rare mtDNA lineages in the gene pool of widespread species would suggest that hybridization and introgression in the D. longispina species complex is frequent even across large phylogenetic distances, and that discoveries of such distinct clades must be interpreted with caution. However, maintenance of ancient polymorphisms through selection is another plausible alternative that may cause the observed discordance and cannot be entirely excluded.

Background

Since the plea for DNA taxonomy [1] and the introduction of the DNA barcoding concept [2] the use of mitochondrial DNA (mtDNA) as genetic marker in animals has become an inherent part not only of taxonomic research but is also widely used in systematics, phylogeography, population genetics, ecology and evolutionary biology. One initial goal of integrating genetic markers into traditional taxonomy was to uncover hidden diversity [2,3,4]. Especially genetically divergent cryptic species, i.e., those pooled under a single name due to their morphological similarity [5, 6], could be easily detected by the use of molecular markers. Since the introduction of the polymerase chain reaction (PCR) in the 1980s, the number of reports of cryptic species increased exponentially [5].

Interestingly, divergent mitochondrial lineages are still detected even in well-studied groups. This is also true for the freshwater crustacean genus Daphnia (Cladocera: Anomopoda), in which previously undescribed lineages were recovered in most biogeographical regions [7]. Within the Daphnia longispina species complex, a dominant taxonomic group of Daphnia in lake plankton of northern temperate zone, an extensive amount of mitochondrial lineages was reported recently [8,9,10,11], with eleven divergent clades without applicable names. The levels of genetic differentiation of seven lineages reported from Europe [8, 11] are as high or even higher than the genetic differentiation among well-recognized species of that complex. Specifically, the sequence divergence from their closest sister clade, expressed as Kimura 2-parameter distance calculated for the mitochondrial gene for the 12S ribosomal ribonucleic acid (12S rRNA), ranged from 10.8 to 17.6% [8]. It was thus suggested that these lineages may represent distinct biological species [8].

However, not every divergent mitochondrial DNA lineage should be considered as species [e.g., [12,13,14]. Several studies show that the divergence detected at the mitochondrial level is not always supported by nuclear DNA data [e.g., 15, 16]. The reasons for such a discrepancy in mitochondrial and nuclear DNA patterns are diverse [17]. On the one hand, nuclear mitochondrial pseudogenes (numts) or pseudogenes from the mitochondrial genome that originate from duplication events might be amplified rather than the functional target gene [18]. On the other hand, cyto-nuclear discordance may be caused by various biological processes, like hybridization and subsequent introgression of the mitochondrial genome [15, 16, 19], incomplete lineage sorting of ancestral polymorphisms [19,20,21], or selection acting on mitochondrial genes, either direct due to, e.g., environmental factors [22, 23], or through indirect selection induced by maternally inherited symbionts [reviewed by 24]. Thus, firm conclusions about the delimitations of species based on genetic data should be drawn from a combination of different approaches, e.g., combining mtDNA and nuclear DNA, and assessing evidence on the level of gene flow between putative species.

For most of the unnamed European lineages of the D. longispina complex, the available information is scarce, mostly limited to sequences of one or more mitochondrial genes [8, 11]. For one of the lineages from Norway (Fig. 1: lineage N), restriction fragment length polymorphism (RFLP) of the nuclear ribosomal internal transcribed spacer (ITS) suggested its distinctness from other taxa [25]. In contrast, allozyme- and ITS-RFLP-based analyses of Daphnia population samples from a Czech reservoir where one of these divergent lineages was later found (Fig. 1: lineage III), did not provide any evidence for the presence of a previously undiscovered species [11, 26]. However, as the samples from that reservoir were analyzed for a different purpose, it is not known whether any individual of this divergent mitochondrial lineage was actually processed in the previous studies [11, 26].

Fig. 1
figure 1

Comparison of nuclear and mitochondrial DNA patterns in the Daphnia longispina species complex. Factorial correspondence analysis (a) demonstrating the position of the 49 individuals belonging to mtDNA clades I, II, or III in relation to the reference dataset (grey circles) consisting of 312 individuals belonging to D. galeata, D. longispina, and D. cucullata. Genetic relationships are depicted using the first two factors of an FCA based on multilocus genotypes from ten microsatellite loci (see Additional file 1 Fig. S1 for illustrations of factor 1 and 3 as well as 2 and 3). All individuals that were characterized as pure taxa according to their assignment probabilities calculated in Structure 2.3.4 [34] are encased by an oval. Patterns of variation of the 12S rDNA [derived and modified from 8] is shown in the inset Maximum Likelihood tree (b). The six divergent mitochondrial lineages firstly described by [8] are labelled I-VI and the remaining five divergent lineages are labelled according to their origin (N: Norway; J1, J2: Japan; R1, R2: Russia)

Full size image

The main aim of our study was to analyze at least some of the distinct mitochondrial lineages of the Daphnia longispina species complex reported recently from Europe [8], to gather additional evidence in favor or against the hypothesis that they represent cryptic species. Therefore, we studied individuals belonging to three distinct mitochondrial lineages from six geographically distant localities using nuclear microsatellite markers [27] in addition to the mitochondrial gene for the 12S rRNA.

Methods

We analyzed 49 individuals of the Daphnia longispina complex from six different lakes (Table 1) that belonged to one of the previously undescribed lineages [reported in 8]. While some DNA samples used in the previous study [8] were reanalyzed at the nuclear level (lineage I: lake Druksiai, Lithuania; III: Želivka reservoir, Czechia), further samples were collected from Želivka reservoir and from lake Norrviken, Sweden (lineage II) for analysis of mitochondrial and nuclear DNA. In addition, individuals belonging to some of the previously undescribed lineages were discovered in further localities (Seragah, Iran; Neuhofener Altrhein, Germany; and Esch-sur-Sûre, Luxembourg), and these were included in our analysis as well.

Table 1 Data on populations from which analyzed individuals belonging to mtDNA clades I, II and III were obtained.
Full size table

As no distinguishing morphological characters were reported for the respective mitochondrial lineages, we screened for these lineages by sequencing mitochondrial 12S rDNA in randomly selected individuals from the respective populations (in total, 508 individuals were studied at the mitochondrial level; for more information on sample sizes per population see Table 1). Total DNA was extracted either by using HotSHOT [28] or H3 extraction method [according to 29].

We aimed for the amplification of an approximately 600 bp fragment of the mitochondrial gene for 12S rRNA [following 11]. If DNA samples had low quality and the amplification of the target marker failed, shorter fragments that still allowed discrimination within the Daphnia longispina complex were amplified, either according to a previously published protocol [30; fragment length approximately 150 bp] or using newly designed primer pairs (segment 1 with an approximate length of 190 bp: forward 5′ CAGGGTATCTAATCCTGG 3′, reverse 5′ GCGACGGCTGGCACGATTT 3′; segment 2 with an approximate length of 230 bp: forward 5′ CTGCACCTTGACCTGAAGT 3′, reverse 5′ CAGGGTATCTAATCCTGG 3′). All amplification products were sequenced using forward and reverse primer on a capillary sequencer (either 3130xl or 3730 Genetic Analyzer; Applied Biosystems). The sequences were manually checked using Geneious 4 (Biomatters Ltd) with a particular focus on possible presence of double peaks or other ambiguities, and aligned using MUSCLE algorithm [31] implemented in MEGA 6 [32].

Individuals which belonged to one of the previously undescribed mtDNA lineages (clade I = 3 individuals, clade II = 11 individuals, and clade III = 35 individuals) were subsequently analyzed at the nuclear level by characterizing the variation at ten microsatellite loci (DaB17/17, Dgm105, Dgm109, Dgm112, Dp196NB, Dp281NB, Dp519, SwiD6, SwiD14, SwiD18) that were originally developed for D. galeata, D. longispina, D. cucullata and other taxa of the D. longispina complex [27, 33]. Each individual was assigned to a multilocus genotype (MLG), based on the allelic variation at all ten loci. We compared the allelic variation with a reference dataset consisting of 312 individuals of D. galeata, D. longispina, and D. cucullata from 21 European locations [see Additional file 1: Table S1 for information on reference populations; partly published in 33].

To test if individuals with previously undescribed mtDNA lineages were also divergent at the nuclear level, we used a model-based assignment test implemented in Structure 2.3.4 [34] to identify the number of genetic clusters in the total dataset (i.e., 312 reference individuals and 49 individuals belonging to the three divergent mitochondrial lineages). As we were only interested in the interspecific patterns, we limited the number of expected groups to a range of K = 1–10 (the upper bound already exceeding the number of expected taxa within the dataset).

For each value of K we conducted ten independent replicates using 100,000 burn-in iterations and 900,000 Markov Chain Monte Carlo (MCMC) steps and we used the admixture modus (ancestry model) as well as independent allele frequencies (allele frequency model). For the assessment of the most likely value of K, we employed the ΔK method [35] implemented in Structure Harvester [36] and evaluated Ln P(D), convergence between replicates as well as individual assignment probabilities. Individuals with a probability to belong to a certain cluster exceeding 90% were regarded as pure-taxon individuals while the others were regarded as hybrids or individuals showing some level of backcrossing or introgression [according to 37]. In addition to the model-based assignment test, the pattern of similarity among individuals from the whole dataset, based on their MLGs, was characterized using an factorial correspondence analysis (FCA) implemented in Genetix 4.05 [38].

Furthermore, we conducted a hierarchical Structure analysis using all individuals assigned to D. longispina, including 103 reference individuals, eight individuals belonging to clade II and nine individuals belonging to clade III. The number of K was set to 1–15 with ten replicates for each K using 50,000 burn-in iterations and 450,000 MCMC steps for each replicate. As in the preceding Structure analysis description, we used the admixture modus as well as independent allele frequencies among populations and employed ΔK, Ln P(D) and convergence between replicates to assess the most likely number of K in the D. longispina data set. With this additional analysis, we wanted to evaluate if within one taxon a substructure is apparent that would support the distinct status of the divergent lineages.

Results

We sequenced 508 individuals from the six studied localities and detected 49 individuals belonging to one of the previously undescribed mtDNA clades (Table 1). Three were assigned to clade I (Druksiai, Lithuania), eleven to clade II (three from Norrviken, Sweden; and eight from Seragah, Iran), and 35 to clade III (seven from Želivka, Czechia; 27 from Neuhofener Altrhein, Germany; and one from Esch-sur-Sûre, Luxembourg). Clades I and III were only represented by a single already reported haplotype [8]. Clade II, however, was represented by two haplotypes, one found in the Swedish population [8] and the other one in the Iranian population (GenBank accession number: KY652589). In all locations, other coexisting Daphnia taxa were detected when sequencing the mitochondrial 12S rDNA (Table 1); in three sites (Druksiai, Želivka, and Neuhofener Altrhein), three different hybridizing species coexisted (D. longispina, D. galeata, and D. cucullata).

All microsatellite loci amplified well in all three divergent mtDNA clades. The clustering approach used by Structure 2.3.4 [34] suggested that K = 3 was appropriate for our dataset based on Evanno’s ΔK [35] as well as the evaluation of Ln P(D) and the convergence between replicates (Fig. 2). The three clusters were associated with the three widespread species of the complex, D. galeata, D. longispina, and D. cucullata, and the 49 individuals belonging to mtDNA clade I, II, or III assorted to either one of the above-mentioned species or exhibited an admixed genotype indicative of individuals with hybrid status (Fig. 3). The factorial correspondence analysis (Fig. 1, Additional file 1: Fig. S1) supported this pattern as well. All 49 individuals grouped within the clusters representing D. galeata, D. longispina, and D. cucullata or their hybrids.

Fig. 2
figure 2

Detection of the uppermost hierarchical level of genetic structure in the complete microsatellite dataset. Ln P(D) and convergence between replicates (open circles with error bars) as well as Delta K (filled circles) were used for the detection of the most likely number of K. The graph is modified from the output derived from Structure Harvester [36]. According to these results, K = 3 is adequate to describe the structure in the dataset

Full size image
Fig. 3
figure 3

Results of the admixture analysis of all 49 individuals belonging to mtDNA clades I, II, and III. Analysis was performed using Structure 2.3.4 [34] with K = 3 where each cluster associates to one of the well-recognized species D. galeata (green), D. longispina (blue), and D. cucullata (purple); reference data are not shown. The population identification code as well as the corresponding mtDNA clade are given on the x axis and the posterior probabilities on the y axis. For information on population codes see Table 1

Full size image

Analysis of the structure within D. longispina (including 103 reference individuals and individuals of clade II and III) resulted in a most likely number of seven or eight clusters, respectively (Additional file 1: Figure S2). These clusters corresponded mainly to geographically distinct populations, and the resulting pattern did not support a separation of D. longispina, clade II and clade III (Additional file 1: Figure S3).

Discussion

The distinctness of the 49 focal individuals, each belonging to one of the three studied mtDNA clades according to 12S rDNA, was not confirmed at the nuclear level. On the contrary, all these individuals were assorted to one of the coexisting species of the complex (Fig. 1): D. galeata (clade II), D. longispina (clade II and III), and D. cucullata (clade I). This finding explains the lack of earlier nuclear evidence (based on allozyme and ITS-RFLP studies) of an unknown taxon in the reservoir Želivka [see 8]. An additional study analyzing in detail Daphnia population structure from the same locality by ten microsatellite markers also failed to find any evidence of an additional distinct taxon [39]. Our data also show that the divergent mtDNA lineages are not a short-term phenomenon but rather persist over longer time periods: Lineage III was detected at different times over one growing season in Neuhofener Altrhein (Table 1) and lineage II was still detected in lake Norrviken after almost two decades.

The reason for these discordant patterns of mitochondrial and nuclear DNA could be, in principle, a methodological artifact in case we failed to amplify the functional mtDNA gene under study and sequenced pseudogenes instead. Especially, the so-called numts (nuclear copies of mitochondrial DNA) are reported frequently [18]. These numts are, in general, observed for protein-coding genes [but see 40 for an example of ribosomal DNA] and some of them still look functional (i.e., exhibit no signs of stop codons or frameshifts within the coding region). In addition, cases are known where those numts are diverged enough that they were inappropriately considered as evidence for the presence of different species [see 18 for more information]. Besides numts, mitochondrial pseudogenes due to gene duplication have also been observed [41] but they are assumed to be very rare.

However, we consider the pseudogene scenarios an unlikely explanation for the observed patterns. First, the original study reporting the divergent lineages [8] already evaluated whether the secondary structures of obtained sequences correspond to that expected from functional ribosomal RNAs. Second, for some individuals from the reservoir Želivka [8] and the lake Seragah (unpublished data), additional partial sequences of mitochondrial genes for the cytochrome c oxidase subunit 1 (COI) and for the reduced nicotine adenine dinucleotide dehydrogenase subunit 2 (ND2) were obtained, confirming the distinctness of the mitochondrial genomes of those lineages.

There are, however, further explanations for the discordant patterns of mitochondrial and nuclear DNA that should be considered when evaluating the hypothesis that divergent mitochondrial lineages represent cryptic species within the D. longispina complex.

For example, cases of mito-nuclear discordances are observed after hybridization and subsequent introgression of the mitochondrial genome [15, 16, reviewed by 17]. The authors of a recent meta-analysis on freshwater and marine fish emphasized that discordance is expected to be higher in freshwater fish due to greater historical hybridization and introgression within this group [42]. Individuals of the genus Daphnia and especially of the D. longispina complex are known to hybridize regularly [43]. In the D. longispina complex, D. galeata is reported as the most “aggressively” hybridizing taxon [44] and its interspecific hybrids with at least five other taxa (D. mendotae, D. lacustris, D. dentifera, D. longispina, and D. cucullata) have been detected [e.g., 9, 33, 45, 46]. Most hybrids of the D. longispina complex are categorized as F1 hybrids, but also later-generation hybrids and backcrosses have been frequently detected [e.g., 33, [47,48,50], which confirms that gene flow between the parental species, although reduced [51], is not entirely restricted. This is also supported by studies that detected mitochondrial introgression within the D. longispina complex [50, 52,53,54] as well as mitochondrial capture in other Daphnia species [55]. Thus, the discordant patterns found in our study may indeed be explained by introgressive hybridization. However, the source taxa of the mitochondrial genomes are not known, either because individuals of the original taxon have not yet been discovered, or because they are extinct. This would not be an exception; the authors of other studies concluded in comparable cases that mitochondrial introgression did cause cyto-nuclear discrepancies, although they failed to identify the source species [13: leaf beetle Chrysomela, 56: Drosophila quinaria].

Incomplete lineage sorting (ILS) of ancestral polymorphisms is another possible scenario [17, 19] as assumed for example for the mito-nuclear discordance observed in the Drosophila serido haplogroup [21]. Distinguishing between ILS and introgression is difficult and often the only indications might be drawn from genetic data. It was suggested that haplotypes which hold a phylogenetically basal position are hints for ancestral polymorphisms and that no predictable geographical pattern is expected from ILS [19]. However, these criteria are inadequate to distinguish ancient mitochondrial introgression from ILS. Most of the previously reported divergent mtDNA lineages in the D. longispina complex [8] are phylogenetically basal compared to the mtDNA lineages of well-recognized species and so far no geographic pattern is visible (but the number of observations is very low at the moment). This would support the idea that the mtDNA variation we observed is due to ILS according to the above-mentioned criteria. However, the divergence among mtDNA lineages in this complex [8] is so high (10.8 to 17.6% sequence divergence from their closest sister clade expressed as Kimura 2-parameter distance) that stochastically driven ILS is highly unlikely as lineage sorting in mitochondrial loci progresses fast [19].

Direct balancing selection may preserve ancestral polymorphisms within a population indefinitely [19]. Often, it is assumed that mitochondrial DNA is selectively neutral; however, various recent studies show deviations from these assumptions [e.g. 22, 23, 57] and propose that mitochondrial DNA variation might be determined also by natural selection. Even in the Daphnia longispina complex, experimental evidence was collected indicating a higher fitness of certain haplotypes under various temperature regimes [58]. Additionally, simulations have shown that uniparentally inherited loci may display strong phylogeographic structures generated by weak selection [59]. Therefore, it is possible that we detected rare ancestral mitochondrial lineages in the D. longispina species complex that are favored under specific environmental conditions, like temperature, predation, food quality or quantity. However, we have no experimental evidence yet to verify the potential adaptive value of the divergent mtDNA lineages.

Besides direct selection, indirect selection induced by maternally inherited symbionts (including parasites), might also explain mito-nuclear discordance [reviewed by 24]. In these cases, certain mitochondrial genomes might be linked to the transmission of those symbionts. The maintenance of different symbionts via natural selection could, therefore, result in mitochondrial differentiation between populations despite ongoing gene flow. However, also selective sweeps, occurring in one population but not in others, might induce this pattern [24]. This is reported for the intracellular bacteria of the genus Wolbachia (which is not known to infect Daphnia [60]), but a few studies of other vertically transmitted symbionts indicate that this is not genus-specific [24]. Most symbionts so far studied in Daphnia (mostly microparasites) are nevertheless transmitted horizontally and not vertically [61], even though many Daphnia species, including those of the D. longispina complex [e.g., 62], carry them. Intracellular parasites that are known to transmit vertically in Daphnia are species of the microsporidian genus Hamiltosporidium [63]. Yet, to our knowledge, there is no published study presenting a case of mito-nuclear discordance caused by a microsporidian. Irrespective thereof, the mitochondrial divergence caused by such a process is expected to be lower than the deep divergence we observed in Daphnia.

Conclusions

In the studied populations, we observed high mitochondrial divergence in the absence of nuclear divergence in the Daphnia longispina species complex. For the reasons discussed above, we would exclude methodological artifacts (like the amplification and sequencing of numts or pseudogenes instead of the functional gene) as explanation. In addition, the level of divergence among lineages is so high that we find stochastic ILS as well as indirect selection that arose from linkage disequilibrium with maternally transmitted symbionts highly unlikely. Likely explanations are that we observe polymorphisms resulting from past introgression from unknown source taxa, or maintained through selection. Both explanations are plausible, as hybridization and introgression are common phenomena in this species complex [43] and experimental evidence hints that selection of mitochondrial genomes may occur in Daphnia [58].

Altogether, we consider past hybridization and introgression as most probable explanation for the observed polymorphisms in the lineages I, II and III, assuming that the original parental lineages did not evolve efficient reproductive barriers. This is commonly observed also for other taxa in the D. longispina species complex and long-lasting introgressed lineages have been recorded before; a well-described example is D. mendotae, which arose from hybridization between D. galeata and D. dentifera but subsequently maintained its distinct character [44]. The observed patterns in our study might be the legacy of long-range dispersal events that ended up in mitochondrial DNA introgression rather than establishment of separated lineages. If the presumed parental species of lineages I-III still exist somewhere as entities with distinct nuclear genomes, their distribution is either geographically restricted or scattered [as hypothesized in 8], possibly in understudied regions like Siberia or parts of eastern and north-eastern Europe.

The results of our study illustrate that the discovery of highly divergent mitochondrial DNA lineages does not necessarily translate into the discovery of previously undescribed species. Additional evidence, like the analyses of nuclear markers, morphology, biogeography, or behavior, is mandatory. Detailed studies of divergent mitochondrial lineages frequently reveal relevant, but small, phenotypic differences, supporting the conclusion that they are indeed distinct taxa [e.g., 64, 65]. Also in the genus Daphnia two such species have been recently described from Europe: Daphnia hrbaceki [66] and Daphnia inopinata [67]. Considering the substantial yet undescribed and often unexplored variation within the genus [7], we assume that other new species descriptions will follow. However, our study indicates that high mtDNA divergence represents a precondition, but not a sufficient criterion for species delimitation.

References

  1. Tautz D, Arctander P, Minelli A, Thomas RH, Vogler AP. A plea for DNA taxonomy. Trends Ecol Evol. 2003;18(2):70–4.

  2. Hebert PDN, Cywinska A, Ball SL, DeWaard JR. Biological identifications through DNA barcodes. Proc R Soc Lond Ser B Biol Sci. 2003;270(1512):313–21.

    Article  CAS  Google Scholar 

  3. Hebert PDN, Penton EH, Burns JM, Janzen DH, Hallwachs W. Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. P Natl Acad Sci USA. 2004;101(41):14812–7.

    Article  CAS  Google Scholar 

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

  5. Bickford D, Lohman DJ, Sodhi NS, Ng PKL, Meier R, Winker K, Ingram KK, Das I. Cryptic species as a window on diversity and conservation. Trends Ecol Evol. 2007;22(3):148–55.

    Article  PubMed  Google Scholar 

  6. Pfenninger M, Schwenk K. Cryptic animal species are homogeneously distributed among taxa and biogeographical regions. BMC Evol Biol. 2007;7:121.

  7. Adamowicz SJ, Petrusek A, Colbourne JK, Hebert PDN, Witt JDS. The scale of divergence: a phylogenetic appraisal of intercontinental allopatric speciation in a passively dispersed freshwater zooplankton genus. Mol Phylogenet Evol. 2009;50(3):423–36.

  8. Petrusek A, Thielsch A, Schwenk K. Mitochondrial sequence variation suggests extensive cryptic diversity within Western Palaearctic Daphnia longispina complex. Limnol Oceanogr. 2012;57(4):1838–45.

  9. Ishida S, Takahashi A, Matsushima N, Yokoyama J, Makino W, Urabe J, Kawata M. The long-term consequences of hybridization between the two Daphnia species, D. galeata and D. dentifera, in mature habitats. BMC Evol Biol. 2011;11:209.

  10. Zuykova EI, Bochkarev NA, Katokhin AV. Molecular genetic identification and phylogeny of Daphnia species (Crustacea, Cladocera) from water bodies of the Lake Chany basin. Russ J Genet. 2013;49(2):206–13.

    Article  CAS  Google Scholar 

  11. Petrusek A, Hobæk A, Nilssen JP, Skage M, Černý M, Brede N, Schwenk K. Taxonomic reappraisal of the European Daphnia longispina complex (Crustacea, Cladocera, Anomopoda). Zool Scr. 2008;37(5):507–19.

  12. Hernández-López A, Rougerie R, Augustin S, Lees DC, Tomov R, Kenis M, Çota E, Kullaj E, Hansson C, Grabenweger G, et al. Host tracking or cryptic adaptation? Phylogeography of Pediobius saulius (Hymenoptera, Eulophidae), a parasitoid of the highly invasive horse-chestnut leafminer. Evol Appl. 2012;5(3):256–69.

  13. Mardulyn P, Othmezouri N, Mikhailov YE, Pasteels JM. Conflicting mitochondrial and nuclear phylogeographic signals and evolution of host-plant shifts in the boreo-montane leaf beetle Chrysomela lapponica. Mol Phylogenet Evol. 2011;61(3):686–96.

  14. Dasmahapatra KK, Elias M, Hill RI, Hoffman JI, Mallet J. Mitochondrial DNA barcoding detects some species that are real, and some that are not. Mol Ecol Resour. 2010;10(2):264–73.

    Article  CAS  PubMed  Google Scholar 

  15. Gompert Z, Forister ML, Fordyce JA, Nice CC. Widespread mito-nuclear discordance with evidence for introgressive hybridization and selective sweeps in Lycaeides. Mol Ecol. 2008;17(24):5231–44.

  16. Linnen CR, Farrell BD. Mitonuclear discordance is caused by rampant mitochondrial introgression in Neodiprion (Hymenoptera: Diprionidae) sawflies. Evolution. 2007;61(6):1417–38.

  17. Toews DPL, Brelsford A. The biogeography of mitochondrial and nuclear discordance in animals. Mol Ecol. 2012;21(16):3907–30.

    Article  CAS  PubMed  Google Scholar 

  18. Song H, Buhay JE, Whiting MF, Crandall KA. Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified. P Natl Acad Sci USA. 2008;105(36):13486–91.

    Article  CAS  Google Scholar 

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

  20. Mckay BD, Zink RM. The causes of mitochondrial DNA gene tree paraphyly in birds. Mol Phylogenet Evol. 2010;54(2):647–50.

  21. Franco FF, Lavagnini TC, Sene FM, Manfrin MH. Mito-nuclear discordance with evidence of shared ancestral polymorphism and selection in cactophilic species of Drosophila. Biol J Linn Soc. 2015;116(1):197–210.

    Article  Google Scholar 

  22. Pavlova A, Amos JN, Joseph L, Loynes K, Austin JJ, Keogh JS, Stone GN, Nicholls JA, Sunnucks P. Perched at the mito-nuclear crossroads: divergent mitochondrial lineages correlate with environment in the face of ongoing nuclear gene flow in an Australian bird. Evolution. 2013;67(12):3412–28.

  23. Cheviron ZA, Brumfield RT. Migration-selection balance and local adaptation of mitochondrial haplotypes in rufous-collared sparrows (Zonotrichia capensis) along an elevational gradient. Evolution. 2009;63(6):1593–605.

    Article  PubMed  Google Scholar 

  24. Hurst GDD, Jiggins FM. Problems with mitochondrial DNA as a marker in population, phylogeographic and phylogenetic studies: the effects of inherited symbionts. Proc R Soc Lond Ser B Biol Sci. 2005;272(1572):1525–34.

  25. Skage M, Hobæk A, Ruthová S, Keller B, Petrusek A, Seďa J, Spaak P. Intra-specific rDNA-ITS restriction site variation and an improved protocol to distinguish species and hybrids in the Daphnia longispina complex. Hydrobiologia. 2007;594:19–32.

  26. Dlouhá S, Thielsch A, Kraus RHS, Seda J, Schwenk K, Petrusek A. Identifying hybridizing taxa within the Daphnia longispina species complex: a comparison of genetic methods and phenotypic approaches. Hydrobiologia. 2010;643:107–22.

  27. Brede N, Thielsch A, Sandrock C, Spaak P, Keller B, Streit B, Schwenk K. Microsatellite markers for European Daphnia. Mol Ecol Notes. 2006;6(2):536–9.

    Article  CAS  Google Scholar 

  28. Montero-Pau J, Gómez A, Muñoz J. Application of an inexpensive and high-throughput genomic DNA extraction method for the molecular ecology of zooplanktonic diapausing eggs. Limnol Oceanogr Methods. 2008;6:218–22.

    Article  CAS  Google Scholar 

  29. Schwenk K, Sand A, Boersma M, Brehm M, Mader E, Offerhaus D, Spaak P. Genetic markers, genealogies and biogeographic patterns in the Cladocera. Aquat Ecol. 1998;32:37–51.

  30. Hamrová E, Goliáš V, Petrusek A. Identifying century-old long-spined Daphnia: species replacement in a mountain lake characterised by paleogenetic methods. Hydrobiologia. 2010;643:97–106.

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

  32. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Thielsch A, Völker E, Kraus RHS, Schwenk K. Discrimination of hybrid classes using cross-species amplification of microsatellite loci: methodological challenges and solutions in Daphnia. Mol Ecol Resour. 2012;12(4):697–705.

    Article  CAS  PubMed  Google Scholar 

  34. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59.

    CAS  PubMed  PubMed Central  Google Scholar 

  35. Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14(8):2611–20.

    Article  CAS  PubMed  Google Scholar 

  36. Earl DA, vonHoldt BM. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4(2):359–61.

    Article  Google Scholar 

  37. Vähä JP, Primmer CR. Efficiency of model-based Bayesian methods for detecting hybrid individuals under different hybridization scenarios and with different numbers of loci. Mol Ecol. 2006;15(1):63–72.

    Article  PubMed  Google Scholar 

  38. Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F. GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Laboratoire Génome, Populations, Interactions, CNRS UMR 5171, Université de Montpellier II, Montpellier (France) 1996–2004.

  39. Stodola J. Detailed taxonomic and clonal structure of the Daphnia longispina species complex on the longitudinal gradient of the Želivka reservoir. MSc. thesis. Prague: Charles University in Prague; 2013.

  40. Sword GA, Senior LB, Gaskin JF, Joern A. Double trouble for grasshopper molecular systematics: intra-individual heterogeneity of both mitochondrial 12S-valine-16S and nuclear internal transcribed spacer ribosomal DNA sequences in Hesperotettix viridis (Orthoptera: Acrididae). Syst Entomol. 2007;32(3):420–8.

  41. Campbell NJH, Barker SC. The novel mitochondrial gene arrangement of the cattle tick, Boophilus microplus: fivefold tandem repetition of a coding region. Mol Biol Evol. 1999;16(6):732–40.

    Article  CAS  PubMed  Google Scholar 

  42. Wallis GP, Cameron-Christie SR, Kennedy HL, Palmer G, Sanders TR, Winter DJ. Interspecific hybridization causes long-term phylogenetic discordance between nuclear and mitochondrial genomes in freshwater fishes. Mol Ecol. 2017;26:3116–27.

    Article  CAS  PubMed  Google Scholar 

  43. Schwenk K, Spaak P. Evolutionary and ecological consequences of interspecific hybridization in cladocerans. Experientia. 1995;51(5):465–81.

    Article  CAS  Google Scholar 

  44. Taylor DJ, Sprenger HL, Ishida S. Geographic and phylogenetic evidence for dispersed nuclear introgression in a daphniid with sexual propagules. Mol Ecol. 2005;14(2):525–37.

    Article  CAS  PubMed  Google Scholar 

  45. Hobæk A, Skage M, Schwenk K. Daphnia galeata × D. longispina hybrids in western Norway. Hydrobiologia. 2004;526(1):55–62.

  46. Taylor DJ, Hebert PDN. Cryptic intercontinental hybridization in Daphnia (Crustacea): the ghost of introductions past. Proc R Soc Lond Ser B Biol Sci. 1993;254(1340):163–8.

    Article  Google Scholar 

  47. Wolf HG, Mort MA. Interspecific hybridization underlies phenotypic variability in Daphnia populations. Oecologia. 1986;68:507–11.

    Article  PubMed  Google Scholar 

  48. Yin M, Petrusek A, Seda J, Wolinska J. Fine-scale genetic analysis of Daphnia host populations infected by two virulent parasites - strong fluctuations in clonal structure at small temporal and spatial scales. Int J Parasitol. 2012;42(1):115–21.

    Article  PubMed  Google Scholar 

  49. Yin M, Wolinska J, Gießler S. Clonal diversity, clonal persistence and rapid taxon replacement in natural populations of species and hybrids of the Daphnia longispina complex. Mol Ecol. 2010;19(19):4168–78.

  50. Alric B, Möst M, Domaizon I, Pignol C, Spaak P, Perga ME. Local human pressures influence gene flow in a hybridizing Daphnia species complex. J Evol Biol. 2016;29(4):720–35.

    Article  CAS  PubMed  Google Scholar 

  51. Keller B, Wolinska J, Tellenbach C, Spaak P. Reproductive isolation keeps hybridizing Daphnia species distinct. Limnol Oceanogr. 2007;52(3):984–91.

    Article  Google Scholar 

  52. Jankowski T, Straile D. Comparison of egg-bank and long-term plankton dynamics of two Daphnia species, D. hyalina and D. galeata: potentials and limits of reconstruction. Limnol Oceanogr. 2003;48(5):1948–55.

  53. Brede N, Sandrock C, Straile D, Spaak P, Jankowski T, Streit B, Schwenk K. The impact of human-made ecological changes on the genetic architecture of Daphnia species. Proc Natl Acad Sci USA. 2009;106(12):4758–63.

  54. Rellstab C, Keller B, Girardclos S, Anselmetti FS, Spaak P. Anthropogenic eutrophication shapes the past and present taxonomic composition of hybridizing Daphnia in unproductive lakes. Limnol Oceanogr. 2011;56(1):292–302.

    Article  CAS  Google Scholar 

  55. Marková S, Dufresne F, Manca M, Kotlik P. Mitochondrial capture misleads about ecological speciation in the Daphnia pulex complex. PLoS One. 2013;8(7):e69497.

  56. Dyer KA, Burke C, Jaenike J. Wolbachia-mediated persistence of mtDNA from a potentially extinct species. Mol Ecol. 2011;20(13):2805–17.

    Article  PubMed  Google Scholar 

  57. Irwin DE. Local adaptation along smooth ecological gradients causes phylogeographic breaks and phenotypic clustering. Am Nat. 2012;180(1):35–49.

  58. Brede N. The reconstruction of evolutionary patterns from Daphnia resting egg banks. Ph.D. thesis. Frankfurtam Main: Goethe University; 2008.

  59. Irwin DE. Local adaptation along smooth ecological gradients causes phylogeographic breaks and phenotypicclustering. Am Nat. 2012;180(1):35–49.

  60. Fitzsimmons JM, Innes DJ. No evidence of Wolbachia among Great Lakes area populations of Daphnia pulex (Crustacea : Cladocera). J Plankton Res. 2005;27(1):121–4.

  61. Ebert D. Host-parasite coevolution: insights from the Daphnia-parasite model system. Curr Opin Microbiol. 2008;11(3):290–301.

    Article  CAS  PubMed  Google Scholar 

  62. Wolinska J, Gießler S, Koerner H. Molecular identification and hidden diversity of novel Daphnia parasites from European lakes. Appl Environ Microbiol. 2009;75(22):7051–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Navia D, Mendonca RS, Ferragut F, Miranda LC, Trincado RC, Michaux J, Navajas M. Cryptic diversity in Brevipalpus mites (Tenuipalpidae). Zool Scr. 2013;42(4):406–26.

  64. Lajus D, Sukhikh N, Alekseev V. Cryptic or pseudocryptic: can morphological methods inform copepodtaxonomy? An analysis of publications and a case study of the Eurytemora affinis species complex. Ecol Evol.2015;5(12):2374–85.

  65. Navia D, Mendonca RS, Ferragut F, Miranda LC, Trincado RC, Michaux J, Navajas M. Cryptic diversity inBrevipalpus mites (Tenuipalpidae). Zool Scr. 2013;42(4):406–26.

  66. Popova EV, Petrusek A, Kořínek V, Mergeay J, Bekker EI, Karabanov DP, Galimov YR, Neretina TV, Taylor DJ, Kotov AA. Revision of the Old World Daphnia (Ctenodaphnia) similis group (Cladocera: Daphniidae). Zootaxa. 2016;4161(1):1–40.

  67. Popova EV, Petrusek A, Kořínek V, Mergeay J, Bekker EI, Karabanov DP, Galimov YR, Neretina TV,Taylor DJ, Kotov AA. Revision of the Old World Daphnia (Ctenodaphnia) similis group (Cladocera: Daphniidae).Zootaxa. 2016;4161(1):1–40.

Download references

Acknowledgements

We thank Kęstutis Arbačiauskas for the zooplankton sample from Lake Druksiai, Isabelle Thys for the sample from reservoir Esch-sur-Sûre, and Melanie Sinn, Irene Kregul, Linda Peters, and Jasna Vukić for laboratory support. In addition, we thank four anonymous reviewers for comments improving the manuscript.

Funding

The authors acknowledge the support by the Czech Science Foundations (project P506/10/P167) and by the Biodiversity and Climate Research Centre Frankfurt am Main (BiKF; “LOEWE–Landes-Offensive zur Entwicklung Wissenschaftlich-ökonomischer Exzellenz” of Hesse’s Ministry of Higher Education, Research and the Arts). In addition, AK received a support from the Charles University Mobility Fund.

Availability of data and materials

The microsatellite dataset of the 49 individuals belonging to the divergent mitochondrial lineages I, II and III as well as the reference data are published in Dryad Digital Repository (doi:https://doi.org/10.5061/dryad.6jd34) as input files for the programs Structure and Genetix. The Structure parameter files as well as the results of both analyses (complete dataset and D. longispina subset) are also stored in Dryad Digital Repository. The sequence data are stored in GenBank (accession number KY652589).

Author information

Authors and Affiliations

  1. Institute for Environmental Sciences, Molecular Ecology, University of Koblenz-Landau, Landau in der Pfalz, Germany

    Anne Thielsch, Alexis Knell & Klaus Schwenk

  2. Faculty of Science, Department of Biology, Ferdowsi University of Mashhad, Mashhad, Iran

    Ali Mohammadyari

  3. Faculty of Science, Department of Ecology, Charles University, Prague, Czechia

    Adam Petrusek

Authors
  1. Anne Thielsch
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Alexis Knell
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Ali Mohammadyari
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Adam Petrusek
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Klaus Schwenk
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

AP, AT and KS had the original idea for the study. AM, AP and AT obtained the samples. AK, AM and AT conducted the genetic analyses. AK and AT carried out the statistical analyses, interpreted the results with input from AP and KS. AT wrote the manuscript with input from all authors. All authors read and approved the final version.

Corresponding author

Correspondence to Anne Thielsch.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional file

Additional file 1:

Mito-nuclear discordance in a Daphnia species complex. The additional file contains information on reference populations (Table S1), displays the comparison of nuclear and mitochondrial DNA patterns in the Daphnia longispina species complex using the first and the third as well as the second and the third factor of a factorial correspondence analysis based on multilocus genotype data from ten microsatellite loci (Figure S1), and presents the results of the hierarchical Structure analysis of individuals belonging to D. longispina, clade II and III that were assigned to D. longispina by the Structure analysis using the whole dataset (Figures S2 and S3). (DOCX 420 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Thielsch, A., Knell, A., Mohammadyari, A. et al. Divergent clades or cryptic species? Mito-nuclear discordance in a Daphnia species complex. BMC Evol Biol 17, 227 (2017). https://doi.org/10.1186/s12862-017-1070-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12862-017-1070-4

Keywords

BMC Ecology and Evolution

ISSN: 2730-7182

Contact us