Local and regional founder effects in lake zooplankton persist after thousands of years despite high dispersal potential

DNA metabarcoding reveals impact of local recruitment, dispersal, and hydroperiod on assembly of a zooplankton metacommunity

Understanding the environmental impact on the assembly of local communities in relation to their spatial and temporal connectivity is still a challenge in metacommunity ecology. This study aims to unravel underlying metacommunity processes and environmental factors that result in observed zooplankton communities. Unlike most metacommunity studies, we jointly examine active and dormant zooplankton communities using a DNA metabarcoding approach to overcome limitations of morphological species identification. We applied two-fragment (COI and 18S) metabarcoding to monitor communities of 24 kettle holes over a two-year period to unravel (i) spatial and temporal connectivity of the communities, (ii) environmental factors influencing local communities, and (iii) dominant underlying metacommunity processes in this system. We found a strong separation of zooplankton communities from kettle holes of different hydroperiods (degree of permanency) throughout the season, while the community composition within single kettle holes did not differ between years. Species richness was primarily dependent on pH and permanency, while species diversity (Shannon Index) was influenced by kettle hole location. Community composition was impacted by kettle hole size and surrounding field crops. Environmental processes dominated temporal and spatial processes. Sediment communities showed a different composition compared to water samples but did not differ between ephemeral and permanent kettle holes. Our results suggest that communities are mainly structured by environmental filtering based on pH, kettle hole size, surrounding field crops, and permanency. Environmental filtering based on specific conditions in individual kettle holes seems to be the dominant process in community assembly in the studied zooplankton metacommunity.

Survivors and colonizers: Contrasting biogeographic histories reconciled in the Antarctic freshwater copepod Boeckella poppei

Two main hypotheses have been proposed to explain the contemporary distribution of Antarctic terrestrial biota. We assess whether the current distribution of maritime Antarctic populations of the freshwater copepod Boeckella poppei is the result of (1) a post-Last Glacial Maximum (LGM) colonization, or whether (2) the species survived in regional glacial refugia throughout the LGM and earlier glaciations. Using 438 specimens from 34 different sampling sites across Southern South America, South Georgia, South Orkney Islands, South Shetland Islands, and the Antarctic Peninsula, we analyzed mitochondrial and nuclear sequences to uncover patterns of genetic diversity and population structure. We also performed median-joining haplotype network, phylogenetic reconstruction, and divergence time analyses. Finally, we evaluated past demographic changes and historical scenarios using the Approximate Bayesian Computation (ABC) method. Our data support the existence of two clades with different and contrasting biogeographic histories. The first clade has been present in maritime Antarctica since at least the mid-Pleistocene, with the South Orkney Islands the most likely refugial area. The second clade has a broader distribution including southern South America, South Georgia, South Shetland Islands, and the Antarctic Peninsula. The ABC method identified long-distance dispersal (LDD) colonization event(s) from southern South America to South Georgia and the maritime Antarctic after the LGM deglaciation, supporting more recent colonization of Antarctic locations. The current Antarctic and sub-Antarctic distribution of B. poppei is likely derived from two independent biogeographic events. The combination of both (1) post-LGM colonization from southern South America and (2) longer-term persistence in in situ regional refugia throughout glacial periods challenges current understanding of the biogeographic history of Antarctic freshwater biota. Re-colonization of ice-impacted Antarctic areas would have occurred following a LDD and Establishment model, pointing to the existence of possible post-dispersal barriers, despite widely assumed high passive dispersal capacity in freshwater invertebrates.

Long-read sequencing reveals the evolutionary drivers of intra-host diversity across natural RNA mycovirus infections

Intra-host dynamics are a core component of virus evolution but most intra-host data come from a narrow range of hosts or experimental infections. Gaining broader information on the intra-host diversity and dynamics of naturally occurring virus infections is essential to our understanding of evolution across the virosphere. Here we used PacBio long-read HiFi sequencing to characterize the intra-host populations of natural infections of the RNA mycovirus Cryphonectria hypovirus 1 (CHV1). CHV1 is a biocontrol agent for the chestnut blight fungus (Cryphonectria parasitica), which co-invaded Europe alongside the fungus. We characterized the mutational and haplotypic intra-host virus diversity of thirty-eight natural CHV1 infections spread across four locations in Croatia and Switzerland. Intra-host CHV1 diversity values were shaped by purifying selection and accumulation of mutations over time as well as epistatic interactions within the host genome at defense loci. Geographical landscape features impacted CHV1 inter-host relationships through restricting dispersal and causing founder effects. Interestingly, a small number of intra-host viral haplotypes showed high sequence similarity across large geographical distances unlikely to be linked by dispersal.

Freshwater zooplankton metapopulations and metacommunities respond differently to environmental and spatial variation

Theory predicts that population genetic structure and metacommunity structure are linked by the common processes of drift and migration, but how population genetic structure and metacommunity structure are related in nature is still unknown. Deeper understanding of the processes influencing both genetic and community diversity is vital for better predicting how environmental change will impact biodiversity patterns. We examined how crustacean zooplankton and rotifer species' metapopulation genetic structure and metacommunities respond to environmental and spatial variation both within and across four regions of boreal Canada. Metapopulation and metacommunity variation partitioning results were compared within and across the four regions. Metapopulations and metacommunities responded differently to environmental variation and spatial structure both within and across regions, as metapopulations were influenced by different environmental variables compared to metacommunities. At larger spatial scales both metapopulations and metacommunities exhibited greater spatial and environmental structuring, again responding to a different subset of environmental variables. Our findings suggest that even though both genetic and species diversity are linked by the same processes, regional variation in environmental characteristics and spatial structure influence resulting biodiversity patterns differently. To date, no other empirical research has explored relationships between entire metapopulation and metacommunity assemblages at large regional spatial scales.

Breaking ecological barriers: Anthropogenic disturbance leads to habitat transitions, hybridization, and high genetic diversity

Genetic diversity is expected to erode in disturbed habitats through strong selection, local extinctions, and recolonization associated with genetic bottlenecks and restricted gene flow. Despite this general prediction and over three decades of population genetics studies, our understanding of the long-term effect of environmental disturbance on local and regional genetic diversity remains limited. We conducted a population genetic survey of the microcrustacean Daphnia across a landscape subject to anthropogenic stressors from a century of industrial mining. At the local scale we found moderate genetic diversity (i.e., low clonal diversity), characteristic of habitat-specific selective sweeps and local extinctions, but high diversity and strong genetic structure at the regional scale despite the shared watershed of many lakes and exceptional dispersal ability of daphniids. Many habitats experienced changes in species assemblages, with the obligate asexual Daphnia pulex lineages-known only to inhabit ponds-dominating disrupted urban lakes. This habitat transition (pond to lake) was likely facilitated by the disruption of ecological barriers maintaining the genomic separation of these young species. Thus, disrupted habitats can exhibit complex and unexpected genetic patterns of local extinctions and recolonizations, followed by habitat transitions, hybridization and potential speciation events that are difficult to predict and should not be underestimated.

Multiple glacial refugia and contemporary dispersal shape the genetic structure of an endemic amphibian from the Pyrenees

Historical factors (colonization scenarios, demographic oscillations) and contemporary processes (population connectivity, current population size) largely contribute to shaping species' present-day genetic diversity and structure. In this study, we use a combination of mitochondrial and nuclear DNA markers to understand the role of Quaternary climatic oscillations and present-day gene flow dynamics in determining the genetic diversity and structure of the newt Calotriton asper (Al. Dugès, 1852), endemic to the Pyrenees. Mitochondrial DNA did not show a clear phylogeographic pattern and presented low levels of variation. In contrast, microsatellites revealed five major genetic lineages with admixture patterns at their boundaries. Approximate Bayesian computation analyses and linear models indicated that the five lineages likely underwent separate evolutionary histories and can be tracked back to distinct glacial refugia. Lineage differentiation started around the Last Glacial Maximum at three focal areas (western, central and eastern Pyrenees) and extended through the end of the Last Glacial Period in the central Pyrenees, where it led to the formation of two more lineages. Our data revealed no evidence of recent dispersal between lineages, whereas borders likely represent zones of secondary contact following expansion from multiple refugia. Finally, we did not find genetic evidence of sex-biased dispersal. This work highlights the importance of integrating past evolutionary processes and present-day gene flow and dispersal dynamics, together with multilocus approaches, to gain insights into what shaped the current genetic attributes of amphibians living in montane habitats.

Limited long-distance dispersal success in a Western European fairy shrimp evidenced by nuclear and mitochondrial lineage structuring

Anostraca are known by their ability for long-distance dispersal, but the existence in several species of deep, geographically structured mtDNA lineages suggests their populations are subjected to allopatric differentiation, isolation, and prevalence of local scale dispersion. Tanymastix stagnalis is one of the most widespread species of Anostraca and previous studies revealed an unclear geographical pattern of mtDNA genetic diversity. Here, we analyze populations from the Iberian and Italian Peninsulas, Central Europe, and Scandinavia, with the aim to characterize the patterns of genetic diversity in a spatio-temporal framework using mtDNA and nuclear markers to test gene flow among close populations. For these aims we built a time-calibrated phylogeny and carried out Bayesian phylogeographic analyses using a continuous diffusion model. Our results indicated that T. stagnalis presents a deeply structured genetic diversity, including 7 ancient lineages, some of them even predating the Pleistocene. The Iberian Peninsula harbors high diversity of lineages, with strong isolation and recent absence of gene flow between populations. Dispersal at local scale seems to be the prevailing dispersal mode of T. stagnalis, which exhibits a pattern of isolation-by-distance in the Iberian Peninsula. We remark the vulnerability of most of these lineages, given the limited known geographic distribution of some of them, and the high risk of losing important evolutionary potential for the species.

Facing Adversity: Dormant Embryos in Rotifers

An in-depth look at the basic aspects of dormancy in cyclic parthenogenetic organisms is now possible thanks to research efforts conducted over the past two decades with rotifer dormant embryos. In this review, we assemble and compose the current knowledge on four central themes: (1) distribution of dormancy in animals, with an overview on the phylogenetic distribution of embryo dormancy in metazoans, and (2) physiological and cellular processes involved in dormancy, with a strong emphasis on the dormant embryos of cyclically parthenogenetic monogonont rotifers; and discussions of (3) the selective pressures and (4) the evolutionary and population implications of dormancy in these animals. Dormancy in metazoans is a widespread phenomenon with taxon-specific features, and rotifers are among the animals in which dormancy is an intrinsic feature of their life cycle. Our review shows that embryo dormancy in rotifers shares common functional pathways with other taxa at the molecular and cellular level, despite the independent evolution of dormancy across phyla. These pathways include the arrest of similar metabolic routes and the usage of common metabolites for the stabilization of cellular structures and to confer stress resistance. We conclude that specific features of recurrent harsh environmental conditions are a powerful selective pressure for the fine-tuning of dormancy patterns in rotifers. We hypothesize that similar mechanisms at the organism level will lead to similar adaptive consequences at the population level across taxa, among which the formation of egg banks, the coexistence of species, and the possibility of differentiation among populations and local adaptation stand out. Our review shows how studies of rotifers have contributed to improved knowledge of all of these aspects.

Unexpected endemism in the Daphnia longispina complex (Crustacea: Cladocera) in Southern Siberia

The biological significance of regional cladoceran morphotypes in the montane regions of the central Palearctic remains poorly understood. In the Holarctic Daphnia longispina complex (Cladocera: Daphniidae), several variants, lineages and species have been proposed as endemic for Southern Siberia. Daphnia turbinata Sars, for example, named after its unusual head shape, is known only from Southern Siberia. Here we sequence DNA of Daphnia from three mitochondrial genes (12S rRNA, 16S rRNA, and NADH dehydrogenase subunit 2, ND2) from 57 localities in Russia and Mongolia (the majority being from Southern Siberia) and place them in evolutionary context with existing data. Our aim was to examine regional endemism of the Daphnia longispina complex in Southern Siberian; to improve the phylogenetic understanding with improved taxonomic and regional sampling, and to better understand the influence of Pleistocene glaciation on the biogeography of these lineages. At least three lineages showed genetic evidence for endemism in Southern Siberia. There was strong support for D. turbinata as a sister lineage to to D. longispina/D. dentifera. Another endemic, Siberian D. cf. longispina, is a sister group to the longispina group in general. Within D. longispina s. str. there was an endemic Siberian clade with a western range boundary near the Yenisei River Basin. Gene flow estimates among populations (based on F_ST values) were very low for clades of D. longispina on a regional (the original 12S dataset), and on a pan-Eurasian (the extended 12S dataset) scale. Negative values of Fu’s F_S and Tajima’s D tests prevailed for the species examined with significant values found for two D. longispina clades, D. dentifera, D. galeata and D. cristata. Our results support the notion that Southern Siberia is an important biogeographic region for cladocerans as it contained unexpected diversity of endemics (such as D. turbinata, D. cf. longispina and lineages of D. umbra and D. longsipina s.str.) and from being the geographic meeting place of expanding postglacial lineages from eastern and western refugia.

Contrasting patterns of divergence at the regulatory and sequence level in European Daphnia galeata natural populations

Understanding the genetic basis of local adaptation has long been a focus of evolutionary biology. Recently, there has been increased interest in deciphering the evolutionary role of Daphnia's plasticity and the molecular mechanisms of local adaptation. Using transcriptome data, we assessed the differences in gene expression profiles and sequences in four European Daphnia galeata populations. In total, ~33% of 32,903 transcripts were differentially expressed between populations. Among 10,280 differentially expressed transcripts, 5,209 transcripts deviated from neutral expectations and their population-specific expression pattern is likely the result of local adaptation processes. Furthermore, a SNP analysis allowed inferring population structure and distribution of genetic variation. The population divergence at the sequence level was comparatively higher than the gene expression level by several orders of magnitude consistent with strong founder effects and lack of gene flow between populations. Using sequence homology, the candidate transcripts were annotated using a comparative genomics approach. Additionally, we also performed a weighted gene co-expression analysis to identify population-specific regulatory patterns of transcripts in D. galeata. Thus, we identified candidate transcriptomic regions for local adaptation in this key species of aquatic ecosystems in the absence of any laboratory-induced stressor.

Contrasting phylogeographic patterns and demographic history in closely related species of Daphnia longispina group (Crustacea: Cladocera) with focus on North-Eastern Eurasia

Species with large geographic distributions present a challenge for phylogeographic studies due to the logistic difficulties of obtaining adequate samples. Daphnia O.F. Müller (Anomopoda: Daphniidae) is a model genus for evolutionary biology and ecology, but many regions such as the remote areas of Siberia, remain poorly studied. Here we examined genetic polymorphism in the ribosomal 12S and the protein-coding ND2 mitochondrial genes of three closely related taxa of the Daphnia (Daphnia) longispina complex, namely D. galeata Sars, D. longispina O.F. Müller and D. dentifera Forbes. We estimated the phylogenetic relationships among these taxa based on a concatenated alignment of these two genes. Using sequences from the present study and those available in GenBank, we investigated the geographic distributions of the mitochondrial haplotypes of these species and proposed an evolutionary scenario for each taxon. Network structures, haplotype distribution patterns, and FST values indicated significant differences in the evolutionary history of the examined species. Our analysis of D. galeata populations confirmed its recent and fast expansion, without a previous phase of a strong population disconnection. In contrast, the high haplotype diversity in D. dentifera and D. longispina could be explained by the survival of different phylogroups in several glacial refugia located in different geographic regions. For all studied species, maximum haplotype diversity was recorded in the remote regions of Siberia-lakes of the Yenisei River and Transbaikalia. Our study is an important step in our understanding of the evolutionary history of the Daphnia longispina group and provides further evidence of the biogeographic significance of Siberia for freshwater taxa.

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

BACKGROUND

The zooplanktonic cladocerans Daphnia, present in a wide range of water bodies, are an important component of freshwater ecosystems. In contrast to their high dispersal capacity through diapausing eggs carried by waterfowl, Daphnia often exhibit strong population genetic differentiation. Here, to test for common patterns in the population genetic structure of a widespread Holarctic species, D. galeata, we genotyped two sets of populations collected from geographically distant areas: across 13 lakes in Eastern China and 14 lakes in Central Europe. The majority of these populations were genotyped at two types of markers: a mitochondrial gene (for 12S rRNA) and 15 nuclear microsatellite loci.

RESULTS

Mitochondrial DNA demonstrated relatively shallow divergence within D. galeata, with distinct haplotype compositions in the two study regions but one widely distributed haplotype shared between several of the Chinese as well as European populations. At microsatellite markers, clear separation was observed at both large (between China and Europe) and small (within Europe) geographical scales, as demonstrated by Factorial Correspondence Analyses, Bayesian assignment and a clustering method based on genetic distances. Genetic diversity was comparable between the sets of Chinese and European D. galeata populations for both types of markers. Interestingly, we observed a significant association between genetic distance and geographical distance for D. galeata populations in China but not in Europe.

CONCLUSIONS

Our results indicate relatively recent spread of D. galeata across wide expanses of the Palaearctic, with one mtDNA lineage of D. galeata successfully establishing over large distances. Despite a clear differentiation of Chinese and European D. galeata at a nuclear level, the pattern of genetic variation is nevertheless similar between both regions. Overall, our findings provide insights into the genetic population structure of a cladoceran species with extremely wide geographical range.

Daphnia diversity on the Tibetan Plateau measured by DNA taxonomy

Daphnia on the Tibetan Plateau has been little studied, and information on species diversity and biogeography is lacking. Here, we conducted a 4-year survey using the barcoding fragment of the mitochondrial COI gene to determine the distribution and diversity of Daphnia species found across the Plateau. Our results show that species richness is higher than previously thought, with total described and provisional species number doubling from 5 to 10. Six of the taxonomic units recovered by DNA taxonomy agreed well with morphology, but DNA barcoding distinguished three clades each for the D. longispina (D. galeata, D. dentifera, and D. longispina) and D. pulex (D. pulex, D. cf. tenebrosa, and D. pulicaria) complexes. The sequence divergence between congeneric species varied within a large range, from 9.25% to 30.71%. The endemic D. tibetana was the most common and widespread species, occurring in 12 hyposaline to mesosaline lakes. The lineage of D. longispina is the first confirmed occurrence in west Tibet.

Diversification of mitochondrial genome of Daphnia galeata (Cladocera, Crustacea): Comparison with phylogenetic consideration of the complete sequences of clones isolated from five lakes in Japan

To characterize genetic diversity and gene flow among Daphnia galeata populations, the complete nucleotide (nt) sequences of the mitochondrial (mt) DNAs of D. galeata clones isolated from five lakes in Japan (Lakes Shirakaba, Suwa, Kizaki, Kasumigaura, and Biwa) were determined. Comparison of non-synonymous (amino acid altering) substitution rates with synonymous substitution rates of D. galeata mt protein-coding genes demonstrated that ATPase8 and COI genes were the most and least susceptible, respectively, to the evolutional forces selecting the aa substitutions. Several non-synonymous substitutions were found in ATPase8 and ATPase6 even in the comparison that no synonymous substitution was found. Comparison of the total number of nt variations among the mt DNAs suggested the phylogenetic relationship ((((Shirakaba/Suwa, Kizaki), Kasumigaura), Biwa), D. pulex). Maximum-likelihood analysis using the total nt sequences of mt protein-coding genes confirmed this relationship with bootstrap values higher than 98%. All the mtDNAs of the analyzed Japanese D. galeata clones contained a control region of essentially the same structure that is distinct from those of the previously reported European Daphnia species of the D. longispina complex. The two control regions of different structures spread among mtDNAs of the Japanese and European Daphnia species, respectively, probably after the divergence of the Japanese D. galeata under different selection pressures associated with their habitats.

Characterization of genome-wide SNPs for the water flea Daphnia pulicaria generated by genotyping-by-sequencing (GBS)

The keystone aquatic herbivore Daphnia has been studied for more than 150 years in the context of evolution, ecology and ecotoxicology. Although it is rapidly becoming an emergent model for environmental and population genomics, there have been limited genome-wide level studies in natural populations. We report a unique resource of novel Single Nucleotide Polymorphic (SNP) markers for Daphnia pulicaria using the reduction in genomic complexity with the restriction enzymes approach, genotyping-by-sequencing. Using the genome of D. pulex as a reference, SNPs were scored for 53 clones from five natural populations that varied in lake trophic status. Our analyses resulted in 32,313 highly confident and bi-allelic SNP markers. 1,364 outlier SNPs were mapped on the annotated D. pulex genome, which identified 2,335 genes, including 565 within functional genes. Out of 885 EuKaryotic Orthologous Groups that we found from outlier SNPs, 294 were involved in three metabolic and four regulatory pathways. Bayesian-clustering analyses showed two distinct population clusters representing the possible combined effects of geography and lake trophic status. Our results provide an invaluable tool for future population genomics surveys in Daphnia targeting informative regions related to physiological processes that can be linked to the ecology of this emerging eco-responsive taxon.

Colonization history and clonal richness of asexual Daphnia in periglacial habitats of contrasting age in West Greenland

Due to climate change, Arctic ice sheets are retreating. This leads to the formation of numerous new periglacial ponds and lakes, which are being colonized by planktonic organisms such as the water flea Daphnia. This system provides unique opportunities to test genotype colonization dynamics and the genetic assemblage of populations. Here, we studied clonal richness of the Daphnia pulex species complex in novel periglacial habitats created by glacial retreat in the Jakobshavn Isbrae area of western Greenland. Along a 10 km transect, we surveyed 73 periglacial habitats out of which 61 were colonized by Daphnia pulex. Hence, for our analysis, we used 21 ponds and 40 lakes in two clusters of habitats differing in age (estimated <50 years vs. >150 years). We tested the expectation that genetic diversity would be low in recently formed (i.e. young), small habitats, but would increase with increasing age and size. We identified a total of 42 genetically distinct clones belonging to two obligately asexual species of the D. pulex species complex: D. middendorffiana and the much more abundant D. pulicaria. While regional clonal richness was high, most clones were rare: 16 clones were restricted to a single habitat and the five most widespread clones accounted for 68% of all individuals sampled. On average, 3·2 clones (range: 1-12) coexisted in a given pond or lake. There was no relationship between clonal richness and habitat size when we controlled for habitat age. Whereas clonal richness was statistically higher in the cluster of older habitats when compared with the cluster of younger ponds and lakes, most young habitats were colonized by multiple genotypes. Our data suggest that newly formed (periglacial) ponds and lakes are colonized within decades by multiple genotypes via multiple colonization events, even in the smallest of our study systems (4 m(2) ).

Geographical genetic structure of Schistosoma japonicum revealed by analysis of mitochondrial DNA and microsatellite markers

BACKGROUND

Schistosoma japonicum is a significant public health risk in parts of China and elsewhere in Southeast Asia. To gain an insight into the epidemiology of schistosomiasis japonica, a detailed investigation of S. japonicum genetic population structure is needed.

METHODS

Using three mitochondrial DNA fragments and ten microsatellite loci, we investigated the genetic diversity within and structure among twelve populations of S. japonicum sampled on a geographical scale covering most major endemic areas.

RESULTS

Schistosoma japonicum lineages from Indonesia, the Philippines and Chinese Taiwan were clearly distinct from each other and from those in mainland China. Within mainland China, there was some evidence for genetic divergence between populations from the mountain and lake regions. However, the analysis inferred no clear sub-population structure in the lake region of mainland China. High genetic diversity was found among S. japonicum populations of mainland China and this was significantly higher than those from island regions.

CONCLUSIONS

High genetic diversity within and substantial differentiation among populations were demonstrated in S. japonicum.

Population structure of a microparasite infecting Daphnia: spatio-temporal dynamics

Background

Detailed knowledge of spatial and temporal variation in the genetic population structure of hosts and parasites is required for understanding of host − parasite coevolution. As hot-spots of contemporary coevolution in natural systems are difficult to detect and long-term studies are restricted to few systems, additional population genetic data from various host − parasite systems may provide important insights into the topic. This is particularly true for parasites, as these players have been under-investigated so far due to the lower availability of suitable molecular markers. Here, we traced genetic variation (based on sequence variants in the internal transcribed spacer region, ITS) among seven geographically isolated populations of the ichthyosporean Caullerya mesnili, a common microparasite of the cladoceran Daphnia (here, the D. longispina hybrid complex). At three sites, we also studied parasite genetic variation over time (three to four sampling points) and tested for associations between parasite genotypes and host species.

Results

Parasite (and host) populations were significantly structured across space, indicating limited dispersal. Moreover, the frequency of parasite genotypes varied significantly over time, suggesting rapid evolutionary change in Caullerya. However, the distribution of parasite genotypes was similar across different host species, which might in turn have important consequences for parasite epidemiology.

Conclusions

The approach proposed here can be applied to track spatial and temporal changes in the population structure of other microparasite species for which sequence variation in the ITS or other highly variable genome regions has been documented but other types of polymorphic markers are lacking. Screening of parasite sequence variants allows for reliable detection of cross-species infections and, using advanced sequencing techniques in the near future, for detailed studies of parasite evolution in natural host − parasite systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0247-3) contains supplementary material, which is available to authorized users.