Continental-scale patterns of hyper-cryptic diversity within the freshwater model taxon Gammarus fossarum (Crustacea, Amphipoda)

Wattier, Remi; Mamos, Tomasz; Copilaş-Ciocianu, Denis; Jelić, Mišel; Ollivier, Anthony; Chaumot, Arnaud; Danger, Michael; Felten, Vincent; Piscart, Christophe; Žganec, Krešimir; Rewicz, Tomasz; Wysocka, Anna; Rigaud, Thierry; Grabowski, Michał

doi:10.1038/s41598-020-73739-0

Download PDF

Article
Open access
Published: 06 October 2020

Continental-scale patterns of hyper-cryptic diversity within the freshwater model taxon Gammarus fossarum (Crustacea, Amphipoda)

Remi Wattier¹,
Tomasz Mamos^2,3,
Denis Copilaş-Ciocianu⁴,
Mišel Jelić⁵,
Anthony Ollivier¹,
Arnaud Chaumot⁶,
Michael Danger⁷,
Vincent Felten⁷,
Christophe Piscart⁸,
Krešimir Žganec⁹,
Tomasz Rewicz^2,10,
Anna Wysocka¹¹,
Thierry Rigaud¹ &
…
Michał Grabowski²

Scientific Reports volume 10, Article number: 16536 (2020) Cite this article

3373 Accesses
44 Citations
5 Altmetric
Metrics details

Subjects

Abstract

Traditional morphological diagnoses of taxonomic status remain widely used while an increasing number of studies show that one morphospecies might hide cryptic diversity, i.e. lineages with unexpectedly high molecular divergence. This hidden diversity can reach even tens of lineages, i.e. hyper cryptic diversity. Even well-studied model-organisms may exhibit overlooked cryptic diversity. Such is the case of the freshwater crustacean amphipod model taxon Gammarus fossarum. It is extensively used in both applied and basic types of research, including biodiversity assessments, ecotoxicology and evolutionary ecology. Based on COI barcodes of 4926 individuals from 498 sampling sites in 19 European countries, the present paper shows (1) hyper cryptic diversity, ranging from 84 to 152 Molecular Operational Taxonomic Units, (2) ancient diversification starting already 26 Mya in the Oligocene, and (3) high level of lineage syntopy. Even if hyper cryptic diversity was already documented in G. fossarum, the present study increases its extent fourfold, providing a first continental-scale insight into its geographical distribution and establishes several diversification hotspots, notably south-eastern and central Europe. The challenges of recording hyper cryptic diversity in the future are also discussed.

Disentangling the taxonomy of the subfamily Rasborinae (Cypriniformes, Danionidae) in Sundaland using DNA barcodes

Article Open access 18 February 2020

Evolving in islands of mud: old and structured hidden diversity in an endemic freshwater crayfish from the Chilean hotspot

Article Open access 21 April 2021

Geographic destiny trumps taxonomy in the Roundtail Chub, Gila robusta species complex (Teleostei, Leuciscidae)

Article Open access 22 September 2023

Introduction

In many areas of biology, including biodiversity assessments, eco-toxicology, environmental monitoring, and behavioural ecology, the species status of the studied organisms relies only on a traditional morphological definition¹. However, a continually increasing number of studies of plants and animals show that one morphospecies might hide lineages with molecular divergences far exceeding what is expected at the intra-specific level in a given taxon^2,3. This hidden diversity, nicknamed cryptic diversity (CD), can reach even tens of lineages for one morphospecies, i.e. showing hyper cryptic diversity (HCD)¹. Even well-studied model-organisms may exhibit an overlooked CD⁴. This hidden diversity occurs at various geographic scales, and it has often been revealed to have pre-Pleistocene origins⁵.

Cryptic diversity has only recently been taken into account in both basic and applied ecology, as well as in evolutionary biology^6,7,8. Even basic evolutionary studies were shown to be impacted by CD if two or more cryptic species are merged into one analysis. Such is the case for sexual selection analyses⁹, studies upon host-parasite co-evolution^6,10,11, functional ecology^12,13 or population genetics¹⁴, and is likely to impact genome assembly, even in the case of model organisms. Similarly, results of biodiversity assessments may be severely underestimated if ignoring the presence of CD¹⁵. At the macroecological scale, considering CD often generates a more nuanced view on determinants of biodiversity patterns⁸, and it should be taken into account while defining conservation units, especially in case of endemic cryptic lineages^16,17. Simulations of distributional patterns under future climatic global change underestimate the real scale of biodiversity loss when CD is not taken into account¹⁸. Finally, ecotoxicological and biomonitoring studies are at risk of being flawed, as some authors reported that different cryptic species might differ in their tolerance to contaminants and ecosystem deterioration^{19,20,21,22,23}.

While CD was detected across a wide range of taxa³, some may be more prone than others². As a rule of thumb, morphospecies with broad geographic distributions, living in patchy habitats and having limited dispersal abilities tend to contain more divergent lineages^2,24,25,26. For example, many freshwater macroinvertebrates²⁷ and especially amphipod crustaceans are known to be associated with a high level of CD^{5,19,28,29,30,31,32,33,34,35,36,37,38,39}.

Here, we focus on the freshwater amphipod morphospecies Gammarus fossarum Koch, 1836, widespread in continental Europe and recently identified in the United Kingdom⁴⁰ (Fig. 1). The species is mostly epigean, being present in a wide range of habitats, but sometimes is found even in shallow groundwaters. It is a keystone species in streams and rivers, where it occurs in high densities and plays a major role in the food web as the primary scavenger/shredder of organic matter as well as a food source for fish and macroinvertebrates (e.g.^41,42). Thus, it is commonly used as a sentinel species in freshwater risk assessment of ecosystem quality, especially in the frame of the European Union's Water Framework Directive (EU WFD)⁴³. It is also used in eco-physiology, e.g. in association with global warming studies^44,45, as well as, extensively, in eco-toxicology^46,47,48,49. This species has also been emerging as a model for molecular reproduction physiology⁵⁰. Finally, G. fossarum is often used in evolutionary ecology studies such as behavioural ecology⁹, host-parasite relationships^10,51 or at the crossroads between parasitology and eco-toxicology^45,52,53.

Already two decades ago, three cryptic lineages (A, B and C) were identified in German populations of G. fossarum⁵⁴, and two of them were shown to have different habitat requirements¹³ and, likely, sensitivity to pollutants^20,21. In recent years, the level of CD was further explored in several parts of Europe and has dramatically increased, establishing G. fossarum as a complex of highly divergent lineages that may represent cryptic species^{55,56,57,58,59}. These studies revealed a substantial amount of CD occurring even at small geographical scales. Intriguingly, the number of lineages seems to be the highest in south-eastern Europe, but the highest phylogenetic diversity recorded so far is in central Europe⁵⁹. However, the extent of diversity in G. fossarum is currently inadequately known and most likely underestimated at the continental level, because of numerous and wide sampling gaps. Indeed, most of the previous studies were restricted geographically compared to the known distribution of G. fossarum in Europe.

Thus, the present study is based on a large number of individuals (almost 5000) and a geographically extensive sampling scheme (ca. 500 sites), including both original and literature data, covering the morphospecies' geographic range in Europe. The mitochondrial DNA cytochrome C oxidase subunit I (COI) gene barcode region sequence⁶⁰ has been chosen to assess the extent of CD. The reasons for that were: (1) sequences of this marker were widely available for G. fossarum in the literature and (2) this marker was successfully used at identifying CD both in G. fossarum^{55,56,57,58,59} and in other gammarids (e.g.²⁸). Following Kekkonen and Hebert⁶¹, we used COI as a quantifier of diversity and an "efficient start for taxonomic workflow", only targeting here the delineation of Molecular Operational Taxonomic Units (MOTUs) as a way to propose testable "species hypotheses" (e.g.¹⁷).

Specifically, we aimed to:

(1)
estimate the level of hidden diversity within G. fossarum at the continental scale using MOTUs delimited with distance and tree-based methods,
(2)
reveal and interpret the critical temporal and geographic patterns of the observed diversification allowed by the pan-European sampling,
(3)
discuss HCD in relation to evolutionary, ecotoxicology, biodiversity surveys and biomonitoring as well as to the future of HCD assessment.

Results

A global overview of MOTU diversity

COI sequences for the DNA-barcode region were either newly generated or derived from the literature for 4926 Gammarus fossarum individuals, from 498 sampling sites distributed throughout 19 European countries (Fig. 1), totalling 691 haplotypes (Table S1, Fig S1). An extraordinarily high level of cryptic diversity was evidenced by a total of 146, 152 and 84 MOTUs obtained with BIN, bPTP and ABGD methods, respectively (Fig. 2). TheK2P distance between haplotypes ranged from almost zero up to 0.35 (Fig. 3a). Although the distribution of genetic distance values was generally bimodal, only a small proportion of values centred around 0.03, while most concentrated around 0.22, illustrating that most haplotypes are highly divergent (Fig. 3a). Only 35 of the 146 BINs have already been known from literature or derived from the sequences mined from GenBank by BOLD, resulting in 111 new BINs in the present study, a fourfold increase. In addition, eight short sequences (397–468 bp) derived from the literature, divergent enough from the existing BINs not to be ascribed to any of them, but too short to be attributed their own BINs (following BOLD standard), should be considered as eight extra BIN-level MOTUs, extending the total to 154 BIN-equivalent MOTUs. Both BIN- and ABGD-MOTUs, shared a bell-shaped distribution ofK2P distances centred on 0.22, illustrating that most MOTUs are highly divergent (Fig. 3b,c). However, while the lowestK2P value among ABGD-MOTUs was 0.05, this value for BIN-MOTUs was occasionally down to < 0.01.

In 62 instances, MOTUs based on BIN, bPTP and ABGD methods were concordant. Single ABGD-MOTUs comprised from two to up to fifteen BINs, clearly showing that BIN was the least conservative among the applied MOTU delimitation methods (Fig. 2).

We observed 32 lineages being both sister clades in the phylogeny (Fig. 2) and diverging by > 0.2 K2P distance. These 32 Lagrue-MOTUs are a conservative estimate of the number of lineages for which we could expect pre-zygotic reproductive isolation as suggested by Lagrue et al.⁵⁷.

Spatial and temporal pattern of diversification

The Bayesian maximum clade credibility chronogram (Fig. 2) based on BINs, each represented by one randomly chosen haplotype, revealed the phylogenetic relationships and divergence times between MOTUs of the G. fossarum morphospecies (Fig. 2). Three significant points are emphasized. First, the estimated age of the entire lineage complex is ca. 26 Mya, i.e. middle-Oligocene. Second, all the major-clades appeared already in the early-middle Miocene, ca. 21–13 Mya. These two features illustrate the ancient age of the G. fossarum complex and its early diversification. Third, the Lineage Through Time analysis (LTT, inset Fig. 2) indicated a rather continuous evolutionary diversification, although at a somewhat slower pace from ca. 12.5 Mya.

Six highly divergent major clades (Fig. 2) SEE, FRA-CE A, EE Q, CWE, CE B and CEE (see below for explanation of acronyms) were recovered according to the rationale presented in Material and Methods , matching the patterns observed in previous multilocus studies. All major clades diverged anciently, around 20 Mya (Figs. 2, 4b). In the following interpretation, a clade was considered as (1) "narrowly distributed" when it was present at only one sampling site or several geographically close sites (< 100 km apart), (2) "broadly distributed" when sites harbouring this clade were > 100 km and < 1000 km apart, and (3) "widespread" when locations were spanning > 1000 km. The same rule was applied to BIN- and ABGD-MOTUs.

Among these six clades, three clades, i.e. CWE (Central-Western Europe), CEE (Central-Eastern Europe) and SEE (South-Eastern Europe), are widespread and rich in BIN- and ABGD-MOTUs. However, the MOTUs harboured by these clades can drastically differ in geographic range size, ranging from narrowly endemic to widespread (Fig. 4a). The CWE clade included 44 BIN- or BIN-equivalent MOTUs and eight ABGD-MOTUs. The clade itself is widespread, ranging from the UK and France in Western Europe, to south-western Poland and western Slovakia in central Europe (Fig. 4). Four ABGD-MOTUs within this clade are broadly distributed. Among them, ABGD-MOTU 3 is restricted to France but split between north-western Alps and south-eastern Massif Central, two areas > 500 km away from each other, each harbouring a specific BIN (Table S1). The ABGD-MOTU 9 is present in Slovakia, Hungary and Poland. The three other ABGD-MOTUs 1, 2 and 6 are more widespread. They would (see “Discussion”) correspond to the G. fossarum types A, B and C, respectively, described by Müller (2000)⁵⁴ (Fig. 2) and is presented in detail below. Only three BINs in this CWE clade were presumably narrowly distributed (Fig. 4a,b, Table S1).

The CEE clade encompassed the highest number of MOTUs: 72 BIN-, or BIN-equivalent MOTUs and 48 ABGD-MOTUs (Fig. 2). Like CWE, the CEE clade is widespread although with a peculiar, discontinuous distribution. While it is present almost exclusively in Central and Eastern Europe, two BINs are nevertheless found in the southern Massif Central in France. The CEE clade includes seven ABGD-MOTUs that are broadly distributed to a various extent (Fig. 4a,b). The ABGD-MOTUs 48 and 61 were present in a few sampling sites ca 250 km away from each other. The ABGD-MOTUs 36, 40, 44 and 45 were present on each side of the Pannonian Basin, (in the south-eastern Alps and the Southern Carpathians), at least ca 500 km apart. The last ABGD-MOTU, 77 is the most widely distributed, being present in Central Europe, namely in Hungary, Slovakia, Poland and eastern Germany. The remaining 48 BINs are presumably endemic to eastern and south-eastern Europe (Fig. 4, Table S1).

The SEE clade harbours 35 BIN- or BIN-equivalent MOTUs and 28 ABGD-MOTUs (Fig. 2). This clade is present in seven countries, from northern Italy in the west, to Albania and Bulgaria in the southeast of Europe. Contrary to CWE and CEE clades, none of these MOTUs was either widespread or even broadly distributed to any extent (Fig. 4a,b).

The major clade FRA-CE A (France-Central-Europe-A) was represented by two ABGD-MOTUs, each being also a single BIN-MOTU (Fig. 2). These two MOTUs are located in two distinct geographical regions, distant by > 1000 km: one in the northern French Alps, the other in western Slovakia and southern Poland (Western Carpathians) (Fig. 4a,b).

The two remaining major clades EE Q (Eastern-Europe-Q) and CE B (Central-Europe-B) were represented only by a single ABGD-MOTU, 78 and 80, respectively, each being also a single BIN-MOTU (Fig. 2). The CE B clade is present only in the Western Carpathians, while the EE Q clade is slightly more broadly distributed in eastern Czech-Republic and westernmost Slovakia (Fig. 4a,b).

ABGD-MOTUs latitudinal and altitudinal distribution

The diversity of ABGD-MOTUs also varied according to altitude (Fig. 4c). Although numerous MOTUs inhabit the zone between 100 and 500 m asl, some ABGD-MOTUs were found below 100 m and many above 600 m, even up to 1500 m. Among the 15 ABGD-MOTUs found to be broadly distributed (Fig. 4d), most had a rather wide altitudinal span. For example, ABGD-MOTU 2 occurs from the sea level to as high as 1500 m asl, while MOTU 78 was found in sub-mountain areas and some ABGD-MOTUs (e.g. 36, 40, 44 and 77) are mainly restricted to lowlands (Fig. 4d).

The diversity of ABGD-MOTUs also varied according to latitude (Fig. 4e), with two peaks – the most prominent from 45.5°N to 46.5°N and another, much lower, from 48.5°N to 50°N. The first one is associated predominantly with the diversity hotspots in the southern part of G. fossarum distribution range, namely the western and south-eastern outskirts of the Alps, the south-western part of the Pannonian Basin and the southern part of the Western Carpathians. The other peak illustrates the lineage diversity present in the Western Carpathians predominantly. The lowest lineage diversity is observed between 50.5°N and 54.0°N, corresponding mainly to the western and central parts of the Great European Plain.

Historical types A, B and C

As shown on Figs. 2 and 4, the ABGD-MOTUs 1, 2 and 6 would correspond, respectively, to the G. fossarum types A, B and C initially defined by Müller⁵⁴ (see “Discussion” for details). Since these MOTUs were the most commonly used for various ecological and ecotoxicological studies (e.g.^20,21), it is important to highlight their biogeographical history. They include large numbers of BINs, 12, 15 and 10, respectively (Figs. 2, 4, 5b–d for details) and their geographic ranges overlap either fully (B and C) or only partially (Fig. 5a). Numerous cases of sympatry between these types were also observed in France, especially in case of ABGD-MOTU 2 (type B) and 6 (type C), and ABGD-MOTU 2 (type B) and 6 (type C) (Fig. 5b–d). Altitudinal ranges for each BIN-MOTUs present in at least three sampled sites are shown in box-plot insets in Fig. 5b–d. The ranges observed at the ABGD-MOTU levels are overall as wide at the BIN-MOTU level for ABGD-MOTUs with high sampling effort. It shows that the ranges of ABGD-MOTUs were not artefacts resulting from combining BIN-MOTU ranges which would have non-overlapping distributions.

Geographic overlap and co-occurrence of BIN-, ABGD-MOTUs and major clades

We found a regional co-occurrence of BIN-MOTUs, ABGD-MOTUs and major clades (Figs. 4, 5). To further quantify within-site co-occurrence, sites with two or more sampled individuals (321 out of 498 sites) were used. Within-site co-occurrence of at least two BIN-MOTUs, ABGD-MOTUs and highly divergent clades was observed in 32.4%, 17.8% and 7.2% of the sites, respectively. Also, at least three BIN-MOTUs co-occurred in 19 locations and at least three ABGD-MOTUs in four sites (Table S1). Sympatry between these different genotypes is therefore quite frequent. Sampling size is likely to impact the ability to detect co-occurrence, particularly if MOTUs are unevenly distributed, which seems a regular feature for G. fossarum (e.g.⁵⁵). For example, in our dataset, rare BIN-MOTU had a frequency of no more than 0.2 in 73.6% of the sites with more than ten individuals sampled, and no more than 0.1 in 49.1% of such sites.

Geographic hotspots of divergence

The landscape genetic map (Fig. 6a) of molecular divergence showed that the Carpathian, Dinaric and Pannonian regions exhibit the highest levels of lineage divergence. This pattern holds regardless if one or more of the six major clades are present in a given region.

If highly divergent endemic Carpathian, Dinaric and Pannonian lineages are excluded (i.e. explaining most of the previously described hotspots), the Western Carpathians still appear as a hotspot of high divergence with five major clades being sympatric in this area: CEE (ABGD-MOTU 77), CWE, EE Q, FRA-CE A and CE B (Fig. 6b). Many regions in France also appeared as local hotspots either, associated with the sympatric occurrence of two major clades, e.g. north-western Alps (CWE and FRA-CE A) and southern Massif Central (CWE and CEE) or divergent ABGD-MOTUs from the same major clade (e.g. 1, 2 and 6, i.e. type A, B and C in central-eastern France) (Fig. 6b). On the contrary, the Great European Plain (northern France, Germany, Poland) of the current distribution of G. fossarum appeared as an area of low divergence (Fig. 6a,b).

Discussion

Our results, based on extensive continental scale sampling of G. fossarum, revealed a remarkable level of hyper cryptic diversity (HCD) with 146, 152 and 84 MOTUs obtained with BIN, bPTP and ABGD methods, respectively. The age of the species complex was estimated at ca. 26 Mya and the divergence was predominantly high, most lineages exhibiting a K2P distance of 0.22. Finally, a relatively high level of sympatry was observed even at the site level. Although CD was already documented in G. fossarum^{54,55,56,57,58,59}, including HCD⁵⁸, the present study considerably increases its extent, offering for the first time the opportunity to discuss its origin on a continental scale.

CD is a common phenomenon in freshwater morphospecies of the genus Gammarus, being reported in all Eurasian wide-ranged taxa studied thus far: e.g. G. balcanicus^5,33,62, G. komareki³², G. ochridensis²⁹, G. pulex^57,63,64, G. roeselii²⁸, G. lacustris³², G. kischineffensis⁶⁴, G. leopoliensis⁶⁵, but also in North American (e.g. G. pecos³⁷) or East Asian (e.g. G. nekkensis⁶⁶) taxa. HCD was detected in three studies^5,28,32. Two studies identified 15 and 35 ABGD MOTUs in G. roeselii²⁶, and in G. komareki³⁰, respectively. A survey covering the entire European range of G. balcanicus⁵ presented 49 highly divergent lineages (likely to represent ABGD-MOTUs given their degree of divergence). In comparison, a previous study⁶² on the same taxon reported up to 22 lineages in the Carpathian Mountains alone. These studies shared a high sampling effort, both in the number of sampling sites and individuals, pointing out the need for appropriate sampling strategies to detect HCD. The age of these species complexes was usually shown to be old, ranging from 4–6 Mya in G. kischineffensis⁶⁴ or G. leopoliensis⁶⁵ to 15–17 Mya for G. balcanicus⁵ and G. roeselii²⁶. Intriguingly, the high frequency of sympatry observed between G. fossarum MOTUs was not detected in the other species complexes^5,32.

Cryptic diversity is common in crustaceans, even overrepresented when compared to other metazoans². The CD in amphipods is far from being restricted to the genus Gammarus. It has been reported for many other epigean freshwater amphipod families, e.g. North American Hyalellidae^30,38,39,67, New Zealand Paracalliopiidae³¹ or Australian Chiltonidae³⁴. It is also extensively known from hypogean freshwater taxa, e.g. Niphargidae^68,69, Paramelitidae⁷⁰, or Crangonyctidae⁷¹, as well as in marine amphipods, including abyssal habitats, e.g. Epimeriidae⁷², Eusiridae⁷³, Eurytheneidae^74,75. Among all these families, examples of both high/old divergence^69,70,75 and HCD^34,38,70 are not rare. Therefore amphipods, whatever the habitat, are highly prone to CD. The present, often decades old, morphological definitions of amphipod species seem insufficiently detailed in many cases. Either morphological stasis, convergent or parallel evolution could be reasons explaining this pattern, but discussing evolutionary processes in detail is beyond the scope of the present paper. It might also be a matter of taxonomic impediment, associated with an increasing shortage of taxonomists⁷⁶, amphipods ranking amongst the least taxonomically studied taxa⁷⁷.

Our data reveal an old and highly complex pattern of molecular diversity within a morphospecies at the local and the continental scale. For example, the ancestors of three clades likely had a pan-European distribution. First, while the CEE clade geographic distribution is centred in central and eastern Europe, three ABGD-MOTUs 68, 71 and 72 were located in southern France, suggesting that the ancestor of this clade had a pan-European distribution. Additional supporting evidence comes from AGBD-MOTU 71, which is in a well-supported sister relationship with the ABGD-MOTUs 65–67, all the latter currently restricted to Romania. The estimated age of this CEE group of lineages is ca 11 Mya, i.e. middle Miocene. Second and similarly, the two MOTUs of the clade FRA-CE A, one restricted to France, and the other to Czechia and Slovakia, are well-supported sister MOTUs, which diverged in Mid-Miocene, ca 12.5 Mya. Third, MOTUs of the CWE clade were present in Western Europe (ABGD-MOTUs 2, 3, 6) and Central Europe (ABGD-MOTUs 1, 9), and also diverged during the Miocene.

Local presence of high CD has already been known for several, predominantly submontane regions of southern and central Europe⁵⁸. We show here a general trend that confirms this pattern and provides further evidence for other local diversity hotspots in the south-western Pannonian Basin and the south-western outskirt of the Alps. However, we point out that the Balkan region is still poorly covered and likely harbours a significant amount of diversity. The probable middle Oligocene (ca. 26 Mya) to late Pliocene (ca. 3 Mya) timeframe of diversification for the main clades and MOTUs corresponds well with the temporal scheme obtained for other widespread European freshwater amphipods, such as G. balcanicus⁵ or G. roeselii²⁸. Our results advocate Alpide orogeny, i.e. simultaneous uplift of the Alps, Carpathians, and Dinarides, as the main factor shaping the phylogeography of G. fossarum. This period was also associated with the presence of extensive lacustrine systems and marshlands covering vast areas of the continent, subsequently fragmented as a result of the regional uplifts⁷⁸. This could enhance the spread of some lineages over vast areas, promoting the subsequent diversification process. That is the case of the sister ABGD-MOTUs 7 (FRA) and 10 (CE A), the former being endemic to the western outskirts of Alps and the latter to the Western Carpathians, while both are in a sister relationship to another highly divergent ABGD-MOTU from the Western Carpathians. We also reveal the presence of wide areas of reduced genetic diversity, predominantly on the Great European Plain, that was partly covered by glaciers during the Last Glacial Maximum (LGM) and partly by a permafrost tundra with few scattered patches of forest⁷⁹. Thus they provided no favourable conditions for most of the Gammarus species, which are known to play a major role in decomposing broadleaved tree debris^41,80. In consequence, the area was probably colonised only after the LGM, as is evidenced already for two largely sympatric BINs of the ABGD-MOTU 77 that spread over the central part of the Great European Plain from the glacial refugium/diversity hotspot in the Western Carpathians⁸¹. Similarly, the other widespread ABGD-MOTU 2 is represented by three BINs in the western part of the Great European Plain. Overall, our findings stress the fundamental role played by a combination of regional historical and environmental factors shaping the local genetic diversity of the G. fossarum complex in different parts of Europe and point out to the need for precise identification of the sampled populations with DNA barcodes.

The G. fossarum types A, B and C were defined in by Müller⁵⁴ in 2000 based on 16S rRNA sequences and allozyme data from samples collected mainly in Germany. More recent publications referring to these three types were based on COI, but the sampling remained restricted to or centred around Germany^20,55,56,63. Our study shows that samples from these studies^20,55,56,63 are associated exclusively with one or two BINs per type (type A: ADH0573²⁰; type B: ACG7784 and ACG8343^20,52,53; type C: ACG8246⁵³). We conclude that our ABGD-MOTUs 1, 2 and 6 correspond to types A, B and C, respectively, each type comprising many BIN level MOTUs (11, 15 and 10, respectively). Our study confirms that the geographic distribution of type A does not extend westward more than already reported in the literature^54,58,82. However, our results substantially extend the distribution of type B (but see⁵⁵) and even more notably type C. Type B was known to be widely present in Germany, Switzerland and north-eastern France^54,57,63 and was recently recorded in the UK⁴⁰. Our study extends its distribution to north-western and southern France. Initial studies restricted the distribution of type C to the French-German border marked by the Rhine River^20,54,63,82. The present study shows it is widely present also in the whole eastern part of France. Some studies pointed out that geography could alone be predictive of the distribution of ABGD level MOTUs, e.g. type A and B within Germany⁵⁴. However, the co-occurrence of individuals associated with different types was also pointed out in the literature (e.g.¹³). Based on our substantial sampling effort and extended geographic cover, we conclude that co-occurrence of MOTUs (ABGDs, BINs) may be the rule for the distributional pattern of G. fossarum. Initially proposed to be of pre-Pleistocene origin^54,82, the evolutionary splits between types A, B and C are in our study suggested at early Miocene, matching the early diversification of the CWE clade itself. Phylogeographic history of type A was already extensively studied⁵⁴. Detailed analysis of other types is beyond the scope of the present paper. However, intricate diversification patterns including multiple refugia and more or less extensive expansion during each interglacial in the Pleistocene and part of the Pliocene are likely, as was already proposed for two stygobiotic amphipods (Niphargus virei⁸³, N. rhenorhodanensis⁸⁴) in the same area.

It is essential that HCD is taken into account in both basic and applied studies while ignoring it might have multiple potential impacts. Recently, accumulating studies upon G. fossarum show that CD may impact various types of biological processes, including (1) sexual selection with impact on the outcome of mate choice and male-male competition⁹, (2) intensity of infection by acanthocephalan parasites^10,11 and parasite-induced mortality¹⁰, (3) functional ecology associated with habitat partitioning¹³ and (4) ecotoxicology, with possibly different response to a contaminant^20,21. It might also be largely involved in the within-species biodiversity effects, increasingly acknowledged in the biodiversity-functioning relationships⁸⁵. Generally, how many different Evolutionary Significant Units (ESU) should be considered within G. fossarum is a relevant question and a challenge as exemplified by other amphipods^35,86. On the one hand, our results suggest that individuals of different types will barely interbreed as they belong to different Lagrue-MOTUs. Indeed, it was revealed that individuals from BINs diverging by ca 16% K2P distance (ACY7276, type C vs either ACY7784 or ACH6569, both type B) predominantly discriminate between each other^9,57. On the other hand, the same authors also found out that individuals from the two studied BINs within type B (ACY7784 or ACH6569) diverging by ca. 4% K2P distance, mate randomly, exhibiting no apparent pre-zygotic reproductive barriers. Our study illustrates that many BINs within each type are sympatric, so the question of how frequent between-MOTU mating remains open but testable. Applying nuclear markers, such as microsatellites (which are available for G. fossarum), might help to answer this question in the future¹⁴.

Our study points out a tremendous haplotypic (691 haplotypes) and MOTU-wise diversity as well as to the co-occurrence of MOTUs at the same sampling site. Also, such syntopy often includes rare MOTUs. Both phenomena will impact the future exploration of CD in G. fossarum, either for basic (e.g. phylogeography) or applied (e.g. monitoring or ecotoxicology) topics. Undoubtedly, the sampling strategy should include enough individuals per site to detect rare MOTUs. The Sanger sequencing of single PCR amplicons⁸⁷, targeting the full COI barcode (as in the present study), might convey highly informative content but, simultaneously, it is a costly low throughput method⁸⁷. Pyrosequencing of mitochondrial 16S targeting diagnostic Single Nucleotide Polymorphisms (SNPs) has been applied to identify type A, B and C⁸². Being high-throughput and fully valid in the restricted geographic area where it was initially implemented (mainly Switzerland), the method would seem inappropriate if applied on a continental scale. More recently, e-DNA-metabarcoding was used to detect the presence of G. fossarum in UK rivers. Environmental DNA (eDNA) was amplified for the entire macroinvertebrate community⁴⁰, and a 333 bp COI mini-barcode was sequenced using high throughput sequencing (HTS⁸⁸). The method proved (1) to be sensitive at detecting G. fossarum even at very low abundance (2.6% of total biomass), and (2) discriminant, as the mini-barcode proved informative enough compared to the full barcode to allow unambiguous phylogenetic assignment at the ABGD-MOTUs level (i.e. type B). Would this mini barcode be sufficiently informative to discriminate all BIN-MOTUs and ABGD-MOTUs diversity at the European scale? Although being outside the scope of the present paper, the large amount of data associated with the current study would easily allow testing it in silico.

Given the HCD described for G. fossarum in the present paper, we recommend any research, either related to applied or basic topics, to obtain both COI barcodes and BINs for the individuals under study. Our research also points out the necessity of a comprehensive DNA barcode reference library⁸⁹ for G. fossarum at the European scale. Even if our research is based on an extensive sampling effort, such sampling is heterogeneous, and many areas of G. fossarum distribution are unsampled or undersampled. Even if the task will be tremendous, particularly given the shortage of morpho-taxonomists, defining and naming the MOTUs would be desirable⁹⁰ as recently occurred for ABGD-MOTU 77, which was ascribed a formal taxonomic name: Gammarus jazdzewskii⁸¹. In the same line, as G. fossarum offers the opportunity to relatively easily set up crossing experiments⁵⁷, a feature rarely explored in amphipods in the context of CD, it should be further explored³⁷.

Methods

Data overview

Sequences for COI were derived from a total of 4926 individuals collected from 498 sites (Fig. 1, Table S1). The sampling included 19 countries, covering the whole range of G. fossarum in continental Europe as well as two sites from United Kingdom, an area considered as recently invaded⁴⁰. The collected data included: (1) original data (206 sites, 2071 sequences), (2) data derived from the co-authors' previously published studies^57,58,59,91 (198 sites, 682 sequences, although 463 sequences from⁵⁴ were initially unreleased in GenBank) and (3) data derived from other published papers^{20,40,55,56,63,92,93} (94 sites and 2173 sequences, although 2086 of these sequences⁵⁶ originated from a very restricted area in Germany). All sequences retrieved from GenBank were checked to be phylogenetically part the Gammarus fossarum species complex using both BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi⁹⁴) and assessing monophyly relative to sister clades⁹³. Sampling size per site was on average 10 individuals but was highly variable with a standard error of 14. A total of 180 sites contained only one individual, 95% originating from already published studies (Table S1). Sampling effort was variable (Fig. 1). Some countries e.g. Slovakia and The Czech Republic were evenly and highly sampled but with low number of samples per site while other countries e.g. Germany accounted for few sites but heavily sampled.

Processing of newly gathered material (sampling, taxonomic identification, molecular methods)

The newly gathered material (Fig. 1) was collected using kick-sampling with a benthic hand-net and fixed in 96% ethanol on site. In the laboratory, the animals were identified to the morphospecies level using characters described in available keys^{95,96,97,98,99}. Genomic DNA was extracted from pieces of muscle tissue using either standard phenol-chloroform¹⁰⁰ or Chelex¹⁰¹ protocols. For all individuals the barcoding region of the mtDNA cytochrome C oxidase subunit I (COI) was amplified for 659–710 bp long fragment (including primers) using any of the following primer pairs: LCO1490 and HCO2198¹⁰², UCOIF and UCOIR¹⁰³, COIGrF and COIGrR2²⁸, LCO1490-JJ and HCO2198-JJ¹⁰⁴. Details about primer sequence and location are given in Table S2. All details of the molecular procedures followed those described by Mamos et al.⁵. The amplicons were sequenced using the BigDye sequencing protocol (Applied Biosystems 3730xl) by Macrogen Inc., Korea. Sequences were edited and aligned with CLUSTALW 2.0 (https://www.clustal.org/)¹⁰⁵ using either MEGA 7 (https://www.megasoftware.net/)¹⁰⁶ or Geneious 6.0.5 (https://www.geneious.com/). Based on a 530 nucleotides long alignment, haplotypes were retrieved using DnaSP v5 (https://www.ub.edu/dnasp/)¹⁰⁷. All newly produced sequences were checked to be phylogenetically part the Gammarus fossarum species complex using both BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi)⁹⁴) and assessing monophyly in relation to sister clades⁹³.

Cryptic diversity—MOTU delimitation

To explore the number of MOTUs that could represent putative cryptic species within G. fossarum, we applied four different approaches. Two methods were purely genetic distance-based methods, i.e. the Automatic Barcode Gap Discovery (ABGD)¹⁰⁸ and the Barcode Index Number (BIN)¹⁰⁹. A third one was a tree based phylogenetic approach, using the Bayesian implementation of the Poisson Tree Processes (bPTP)¹¹⁰. The fourth method takes advantage of the study of Lagrue et al.⁵⁷ which combined genetic distance and reproductive isolation experiments. Both BIN and ABGD share the principle of clustering sequences into MOTUs according to their molecular divergence using algorithms aiming to find discontinuities, i.e. a so-called barcode gap separating intra- and interspecific genetic distances. For the ABGD (https://wwwabi.snv.jussieu.fr/public/abgd/) method, we used primary partitions as a basis for group definition. Primary partitions are typically stable over a wide range of prior values, minimise the number of false-positive (over-split species) and, in the absence of cryptic diversity, are usually found to reflect the number of taxa described by morpho-taxonomists¹¹¹. The Kimura two-parameter (K2P) substitution model was applied⁶⁰, and the default value of 0.001 was used as the minimum intraspecific distance. Given that neither morphology^32,55,103 nor reproductive isolation⁵⁷ provides any consensus about which maximum intraspecific distance allows for reliable species delimitation in amphipods, we explored a set of values of up to 0.1. The BIN algorithm is implemented as part of The Barcode of Life Data Systems (BOLD; https://www.boldsystems.org/)¹¹². In this algorithm, the newly submitted COI sequences are aligned and compared pairwise and also to each sequence already deposited in BOLD. Then, the initial single-linkage clustering is performed on the aligned sequence data in such a way that each cluster is allowed for a maximum intra-cluster distance of 0.022, and its distance to any other cluster is more than twice the threshold (> 0.044). In the subsequent step, the clusters are refined in the sense that clusters whose members show high sequence variation, but lack discontinuity remain as a single MOTU, while those in which sequence variation shows clear internal partitions are assigned to two or more MOTUs, even if their separation is less than 0.022. Then, finally, each refined cluster is assigned a globally unique and specific identifier. For details on the BIN system see Ratnasingham and Hebert (2003)¹⁰⁹. Like ABGD, in the absence of cryptic diversity, BINs are claimed to be congruent with the number of taxa described by morpho-taxonomists¹⁰⁹. In addition, BINs may be registered and publicly available through BOLD, allowing tractability and have similar properties as classic taxonomy descriptors (e.g. can be split or synonymised). Only sequences over 500 bp can be included in the BIN clustering, while shorter sequences which are over 300 bp may be only ascribed to an existing BIN, but will not create a BIN or split an existing one. For example, our data set included eight short sequences derived from the literature, divergent enough from the existing BINs not to be ascribed to any of them but too short to be ascribed their own BIN (following the BOLD standard). In the result, we considered them as separate BIN-equivalent MOTUs.

The tree-based bPTP method uses non-ultrametric phylogenies. It incorporates the number of substitutions in the model of speciation and assumes that the probability that a substitution gives rise to a speciation event follows a Poisson distribution. The branch lengths of the input tree are assumed to be generated by two independent classes of Poisson events, one corresponding to speciation and the other to coalescence¹¹⁰. For the input tree, we used a Maximum likelihood (ML) tree obtained with RAxML HPC 8.2.9 (https://cme.h-its.org/exelixis/web/software/raxml/)¹¹³ from the haplotype data. The ML analysis was run under a thorough tree search with the GTR + G substitution model applied to each codon partition. The bPTP (https://species.h-its.org/ptp/) analysis was performed with 500,000 iterations of MCMC and 10% burn-in. Three runs were performed. The convergence of each run was verified through stationary pattern of the MCMC iterations trace plot. All runs provided congruent results.

The fourth method is newly introduced in the present paper and will be designated as the Lagrue-MOTUs. It takes advantage of the study of Lagrue et al.⁵⁷, which combined COI K2Pdistance and reproductive isolation in Gammarus fossarum. Lagrue et al.⁵⁷ tested two distance classes ca. 4% (for which reproductive isolation was absent) and ca 16% (for which reproductive isolation was almost exclusively present). Conservatively, sister clades that diverged by > 20% K2P distance in the present study were considered as different Lagrue-MOTUs and a surrogate of the minimum number of biological species.

Time-calibrated phylogenetic tree

The time-calibrated phylogenetic tree was reconstructed using Bayesian inference (BI) in BEAST 2.4.8 (https://www.beast2.org/)¹¹⁴. A Gammarus roeselii sequence (KP789695) was used as an outgroup for the analysis. Two sets of data were used, one at the haplotype level the other at the BIN level. The outgroup was removed after the analysis. Priors for the substitution models were selected using bModelTest (https://github.com/BEAST2-Dev/bModelTest/)¹¹⁵. The best-fitting model of substitution in both cases was the transversion model (TVM) with gamma-distributed rate heterogeneity (G) and a given proportion of invariable sites (I). The log-normal relaxed clock with the Birth–Death speciation model was set as priors following the results of path sampling selection¹¹⁶. The recent cross-validated analyses of molecular clock, based on geological events for G. fossarum^58,59 and other gammarids (e. g.^5,28,29), indicated that the widely used arthropod COI rate of 0.0115 substitutions per site per My¹¹⁷ is useful in linking present-day patterns with historical processes. Thus, we have used it to calibrate the molecular clock in our phylogeny reconstruction. Four MCMC chains were run for 50 M iterations and sampled every 5000 iterations. The Effective Sampling Size (ESS) of each parameter was verified to be above 200 in Tracer 1.7 (https://beast.community/tracer)¹¹⁸. The runs were combined in LogCombiner 2.4.8, with 25% burn-in, and the maximum clade credibility chronogram was annotated using TreeAnnotator 2.4.8, both programs being part of BEAST 2.4.8 (https://www.beast2.org/)¹¹⁴.

In order to provide additional support for the BI topology, we have also reconstructed a phylogeny using the Maximum Likelihood approach (ML) in RAxML HPC 8.2.9 (https://cme.h-its.org/exelixis/web/software/raxml/)¹¹³. The best-scoring ML tree was searched under the GTR + G substitution model which was applied to each codon partition. Statistical support was estimated with the GTRCAT model of rate heterogeneity and 1,000 rapid bootstrap (RBS) replicates¹¹⁹.

The history of diversification was explored with the help of a Lineage Through Time (LTT) plot, generated with Tracer 1.7¹¹⁸, using 1,500 post burn-in trees that resulted from the BEAST analysis on the BIN dataset.

Definition of major clades

We defined the major clades following studies from Copilaş-Ciocianu and Petrusek⁵⁸ and Copilaş-Ciocianu et al.⁵⁹ which were based on mitochondrial (COI and 16S) and nuclear markers (18S, 28S and EF1α). These studies have robustly recovered these deep clades in multilocus concatenation and species tree analyses. We acknowledge that, in the present study, node supports are not high in all cases as only the highly variable COI marker was used, which has reduced ability to resolve deep nodes (see “Results”). In addition, the names and acronyms were kept in order to maintain continuity with clade nomenclature from these previous studies^58,59. Such clades refer to geography (i.e. SEE: south-eastern Europe), illustrating rather in which part of Europe a given clade has a centre of its distribution than willing to match any predefined strict geographic units. In addition, it is to be noted that extending the present study analysis at a continental scale might also have been associated with extending the original distribution of some clades.

Geographic patterns of genetic diversity

The spatial diversity pattern of G. fossarum in Europe was illustrated with genetic landscape shapes generated with Alleles in Space (AIS; https://www.marksgeneticsoftware.net) software¹²⁰. The genetic landscape visualizes the abrupt transitions between populations and groups of populations characterized by divergent haplotypes. The analysis followed the approach already applied to amphipods (e.g.⁵), i.e. converting the genetic distances between sites into a hypsometric map. First, using the AIS software, the genetic distances between sampling sites were calculated based on the COI sequences longer than 496 bp (in order to obtain a gap-free alignment) and connected into a network based on the Delaunay Triangulation. The genetic distance values were set in the midpoints of each connection in the network. The raw genetic distances acquired from the program were interpolated afterwards and the matrix of the 'elevation' values, with their respective latitude and longitude coordinates, was then imported into QGIS (https://qgis.osgeo.org) software to produce a genetic divergence surface image using the Triangulated Irregular Network (TIN) algorithm. The resulting image was plotted onto a relief map of Europe to create a final map in which the hypsometric tints reflect the genetic distance between population pairs. Two maps were constructed, one that included all the sampling sites (498) and all haplotypes (691) and another one which excluded the highly divergent lineages endemic to the Carpathian, Dinaric and Pannonian regions (the hotspots of the most ancient divergences, see “Results”), to better illustrate the level of genetic divergence in the remaining area. A total of 336 sites and 469 haplotypes were included in this latter analysis.

Ethical approval

All the applicable international, national, and/or institutional guidelines for the care and use of animals were followed and were in agreement with recommendations from the "Comité d'Ethique et de l'Expérimentation Animale (C2EADijon Grand Campus) from Dijon University, France. No experiments on alive animals were performed for this study.

Data availability

All newly produced COI sequences are released in GenBank under accession numbers: MT978656–MT980726. The 463 sequences from Lagrue et al.⁵⁷ which were initially unreleased in Genbank are now available under accession numbers: MT411018-MT411480. Sequences and metadata are available in BOLD dataset DS-GFOSCDEU (dx.doi.org/10.5883/DS-GFOSCDEU). All the new material used in this study has been stored in the permanent collection of either Biogeosciences Laboratory, University of Burgundy-Franche-Comté or the Department of Invertebrate Zoology and Hydrobiology, University of Lodz, Poland.

References

Adams, M., Raadik, T. A., Burridge, C. P. & Georges, A. Global biodiversity assessment and hyper-cryptic species complexes: more than one species of elephant in the room?. Syst. Biol. 63, 518–533 (2014).
Article PubMed Google Scholar
Pérez-Ponce de León, G. & Poulin, R. Taxonomic distribution of cryptic diversity among metazoans: not so homogeneous after all. Biol. Lett. 12, 20160371. https://doi.org/10.1098/rsbl.2016.0371 (2016).
Article PubMed PubMed Central Google Scholar
Pfenninger, M. & Schwenk, K. Cryptic animal species are homogeneously distributed among taxa and biogeographical regions. BMC Evol. Biol. 7, 121. https://doi.org/10.1186/1471-2148-7-121 (2007).
Article PubMed PubMed Central Google Scholar
Keshavmurthy, S. et al. DNA barcoding reveals the coral “laboratory-rat”, Stylophora pistillata encompasses multiple identities. Sci. Rep. 3, 1520. https://doi.org/10.1038/srep01520 (2013).
Article CAS PubMed PubMed Central Google Scholar
Mamos, T., Wattier, R., Burzyński, A. & Grabowski, M. The legacy of a vanished sea: a high level of diversification within a European freshwater amphipod species complex driven by 15 My of Paratethys regression. Mol. Ecol. 25, 785–810 (2016).
Article Google Scholar
Bickford, D. et al. Cryptic species as a window on diversity and conservation. Trends Ecol. Evol. 22, 148–155 (2007).
Article PubMed Google Scholar
Brodersen, J. & Seehausen, O. Why evolutionary biologists should get seriously involved in ecological monitoring and applied biodiversity assessment programs. Evol. Appl. 7, 968–983 (2014).
Article PubMed PubMed Central Google Scholar
Eme, D. et al. Do cryptic species matter in macroecology? Sequencing European groundwater crustaceans yields smaller ranges but does not challenge biodiversity determinants. Ecography 41, 424–436 (2018).
Article Google Scholar
Galipaud, M., Gauthey, Z., Turlin, J., Bollache, L. & Lagrue, C. Mate choice and male–male competition among morphologically cryptic but genetically divergent amphipod lineages. Behav. Ecol. Sociobiol. 69, 1907–1916 (2015).
Article Google Scholar
Galipaud, M., Bollache, L. & Lagrue, C. Variations in infection levels and parasite-induced mortality among sympatric cryptic lineages of native amphipods and a congeneric invasive species: Are native hosts always losing?. Int. J. Parasitol. Parasites Wildl. 6, 439–447 (2017).
Article PubMed PubMed Central Google Scholar
Westram, A. M., Baumgartner, C., Keller, I. & Jokela, J. Are cryptic host species also cryptic to parasites? Host specificity and geographical distribution of acanthocephalan parasites infecting freshwater Gammarus. Infect. Genet. Evol. 11, 1083–1090 (2011).
Article CAS PubMed Google Scholar
Fišer, Ž, Altermatt, F., Zakšek, V., Knapič, T. & Fišer, C. Morphologically cryptic amphipod species are “ecological clones” at regional but not at local scale: a case study of four Niphargus species. PLoS ONE 10, e0134384. https://doi.org/10.1371/journal.pone.0134384 (2015).
Article CAS PubMed PubMed Central Google Scholar
Eisenring, M., Altermatt, F., Westram, A. M. & Jokela, J. Habitat requirements and ecological niche of two cryptic amphipod species at landscape and local scales. Ecosphere 7, e01319. https://doi.org/10.1002/ecs2.1319 (2016).
Article Google Scholar
Westram, A. M., Jokela, J. & Keller, I. Hidden biodiversity in an ecologically important freshwater amphipod: Differences in genetic structure between two cryptic species. PLoS ONE 8, e69576. https://doi.org/10.1371/journal.pone.0069576 (2013).
Article ADS CAS PubMed PubMed Central Google Scholar
Mutanen, M., Kaila, L. & Tabell, J. Wide-ranging barcoding aids discovery of one-third increase of species richness in presumably well-investigated moths. Sci. Rep. 3, 2901. https://doi.org/10.1038/srep02901 (2013).
Article ADS PubMed PubMed Central Google Scholar
Cook, B. D., Page, T. J. & Hughes, J. M. Importance of cryptic species for identifying ‘representative’ units of biodiversity for freshwater conservation. Biol. Conserv. 141, 2821–2831 (2008).
Article Google Scholar
Fišer, C., Robinson, C. T. & Malard, F. Cryptic species as a window into the paradigm shift of the species concept. Mol. Ecol. 27, 613–635 (2018).
Article PubMed Google Scholar
Bálint, M., Barnard, P. C., Schmitt, T., Ujvárosi, L. & Popescu, O. Differentiation and speciation in mountain streams: a case study in the caddisfly Rhyacophila aquitanica (Trichoptera). J. Zool. Syst. Evol. Res. 46, 340–345 (2008).
Article Google Scholar
Major, K., Soucek, D. J., Giordano, R., Wetzel, M. J. & Soto-Adames, F. The common ecotoxicology laboratory strain of Hyalella azteca is genetically distinct from most wild strains sampled in eastern North America: common lab strain of H. azteca is distinct from wild strains. Environ. Toxicol. Chem. https://doi.org/10.1002/etc.2355 (2013).
Article PubMed Google Scholar
Feckler, A., Thielsch, A., Schwenk, K., Schulz, R. & Bundschuh, M. Differences in the sensitivity among cryptic lineages of the Gammarus fossarum complex. Sci. Total Environ. 439, 158–164 (2012).
Article ADS CAS PubMed Google Scholar
Feckler, A. et al. Cryptic species diversity: an overlooked factor in environmental management?. J. Appl. Ecol. 51, 958–967 (2014).
Article Google Scholar
Caputo, D. R., Robson, S. C., Werner, I. & Ford, A. T. Complete transcriptome assembly and annotation of a critically important amphipod species in freshwater ecotoxicological risk assessment: Gammarus fossarum. Environ. Int. 137, 105319 (2020).
Article CAS PubMed Google Scholar
Cogne, Y. et al. De novo transcriptomes of 14 gammarid individuals for proteogenomic analysis of seven taxonomic groups. Sci. Data 6, 184. https://doi.org/10.1038/s41597-019-0192-5 (2019).
Article CAS PubMed PubMed Central Google Scholar
Costa-Silva, G. J., Rodriguez, M. S., Roxo, F. F., Foresti, F. & Oliveira, C. Using different methods to access the difficult task of delimiting species in a complex Neotropical hyperdiverse group. PLoS ONE 10, e0135075 (2015).
Article PubMed PubMed Central CAS Google Scholar
Dincă, V. et al. DNA barcode reference library for Iberian butterflies enables a continental-scale preview of potential cryptic diversity. Sci. Rep. 5, 12395. https://doi.org/10.1371/journal.pone.0135075 (2015).
Article ADS CAS PubMed PubMed Central Google Scholar
Trontelj, P. et al. A molecular test for cryptic diversity in ground water: how large are the ranges of macro-stygobionts?. Freshw. Biol. 54, 727–744 (2009).
Article CAS Google Scholar
Jackson, J. K. et al. Cryptic biodiversity in streams: a comparison of macroinvertebrate communities based on morphological and DNA barcode identifications. Freshw. Sci. 33, 312–324 (2014).
Article Google Scholar
Grabowski, M., Mamos, T., Bacela-Spychalska, K., Rewicz, T. & Wattier, R. A. Neogene paleogeography provides context for understanding the origin and spatial distribution of cryptic diversity in a widespread balkan freshwater amphipod. PeerJ https://doi.org/10.7717/peerj.3016 (2017).
Article PubMed PubMed Central Google Scholar
Grabowski, M., Wysocka, A. & Mamos, T. Molecular species delimitation methods provide new insight into taxonomy of the endemic gammarid species flock from the ancient Lake Ohrid. Zool. J. Linn. Soc. https://doi.org/10.1093/zoolinnean/zlw025 (2017).
Article Google Scholar
Hogg, I. D., Larose, C., de Lafontaine, Y. & Doe, K. G. Genetic evidence for a Hyalella species complex within the Great Lakes—St. Lawrence River drainage basin: implications for ecotoxicology and conservation biology. Can. J. Zool. 76, 1134–1140 (1998).
Article Google Scholar
Hogg, I. D., Stevens, M. I., Schnabel, K. E. & Ann Chapman, M. Deeply divergent lineages of the widespread New Zealand amphipod Paracalliope fluviatilis revealed using allozyme and mitochondrial DNA analyses. Freshw. Biol. 51, 236–248 (2006).
Article CAS Google Scholar
Katouzian, A.-R. et al. Drastic underestimation of amphipod biodiversity in the endangered Irano-Anatolian and Caucasus biodiversity hotspots. Sci. Rep. 6, 22507. https://doi.org/10.1038/srep22507 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Mamos, T., Wattier, R., Majda, A., Sket, B. & Grabowski, M. Morphological vs. molecular delineation of taxa across montane regions in Europe: the case study of Gammarusbalcanicus Schäferna, (Crustacea: Amphipoda). J. Zool. Syst. Evol. Res. 52, 237–248 (2014).
Article Google Scholar
Murphy, N. P., Adams, M., Guzik, M. T. & Austin, A. D. Extraordinary micro-endemism in Australian desert spring amphipods. Mol. Phylogenet. Evol. 66, 645–653 (2013).
Article CAS PubMed Google Scholar
Murphy, N. P., King, R. A. & Delean, S. Species, ESUs or populations? Delimiting and describing morphologically cryptic diversity in Australian desert spring amphipods. Invertebr. Syst. 29, 457 (2015).
Article Google Scholar
Seidel, R. A., Lang, B. K. & Berg, D. J. Phylogeographic analysis reveals multiple cryptic species of amphipods (Crustacea: Amphipoda) in Chihuahuan Desert springs. Biol. Conserv. 142, 2303–2313 (2009).
Article Google Scholar
Sutherland, D. L., Hogg, I. D. & Waas, J. R. Phylogeography and species discrimination in the Paracalliopefluviatilis species complex (Crustacea: Amphipoda): can morphologically similar heterospecifics identify compatible mates?: Mate discrimination in P. fluviatilis. Biol. J. Linn. Soc. 99, 196–205 (2009).
Article Google Scholar
Witt, J. D. S., Threloff, D. L. & Hebert, P. D. N. DNA barcoding reveals extraordinary cryptic diversity in an amphipod genus: implications for desert spring conservation. Mol. Ecol. 15, 3073–3082 (2006).
Article CAS PubMed Google Scholar
Witt, J. D. S. & Hebert, P. D. N. Cryptic species diversity and evolution in the amphipod genus Hyalella within central glaciated North America: a molecular phylogenetic approach. Can. J. Fish Aquat. Sci. 57, 12 (2000).
Article Google Scholar
Blackman, R. et al. Detection of a new non-native freshwater species by DNA metabarcoding of environmental samples—first record of Gammarus fossarum in the UK. Aquat. Invasions 12, 177–189 (2017).
Article Google Scholar
Dangles, O., Gessner, M. O., Guerold, F. & Chauvet, E. Impacts of stream acidification on litter breakdown: implications for assessing ecosystem functioning. J. Appl. Ecol. 41, 365–378 (2004).
Article CAS Google Scholar
Felten, V., Tixier, G., Guérold, F., De Crespin De Billy, V. & Dangles, O. Quantification of diet variability in a stream amphipod: implications for ecosystem functioning. Fundam. Appl. Limnol. Arch. Für Hydrobiol. 170, 303–313 (2008).
Article Google Scholar
Besse, J.-P., Geffard, O. & Coquery, M. Relevance and applicability of active biomonitoring in continental waters under the Water Framework Directive. TrAC Trends Anal. Chem. 36, 113–127 (2012).
Article CAS Google Scholar
Schmidlin, L., von Fumetti, S. & Nagel, P. Temperature effects on the feeding and electron transport system (ETS) activity of Gammarus fossarum. Aquat. Ecol. 49, 71–80 (2015).
Article Google Scholar
Labaude, S., Rigaud, T. & Cézilly, F. Additive effects of temperature and infection with an acanthocephalan parasite on the shredding activity of Gammarus fossarum (Crustacea: Amphipoda): the importance of aggregative behavior. Glob. Change Biol. 23, 1415–1424 (2017).
Article ADS Google Scholar
Kunz, P. Y., Kienle, C. & Gerhardt, A. Gammarus spp. in aquatic ecotoxicology and water quality assessment: Toward integrated multilevel tests. In Reviews of Environmental Contamination and Toxicology Vol. 205 (ed. Whitacre, D. M.) 1–76 (Springer, New York, 2010).
Google Scholar
Mehennaoui, K. et al. Gammarus fossarum (Crustacea, Amphipoda) as a model organism to study the effects of silver nanoparticles. Sci. Total Environ. 566–567, 1649–1659 (2016).
Article ADS PubMed CAS Google Scholar
Besse, J.-P. et al. Caged Gammarus fossarum (Crustacea) as a robust tool for the characterization of bioavailable contamination levels in continental waters: towards the determination of threshold values. Water Res. 47, 650–660 (2013).
Article CAS PubMed Google Scholar
Gouveia, D. et al. Ecotoxico-Proteomics for aquatic environmental monitoring: First in situ application of a new proteomics-based multibiomarker assay using caged amphipods. Environ. Sci. Technol. 51, 13417–13426 (2017).
Article ADS CAS PubMed Google Scholar
Trapp, J. et al. Proteogenomics of Gammarus fossarum to document the reproductive system of amphipods. Mol. Cell. Proteomics 13, 3612–3625 (2014).
Article CAS PubMed PubMed Central Google Scholar
Sanchez-Thirion, K. et al. High food quality increases infection of Gammarus pulex (Crustacea: Amphipoda) by the acanthocephalan parasite Pomphorhynchus laevis. Int. J. Parasitol. 49, 805–817 (2019).
Article CAS PubMed Google Scholar
Bigot-Clivot, A. et al. Bioaccumulation of toxoplasma and cryptosporidium by the freshwater crustacean Gammarus fossarum: involvement in biomonitoring surveys and trophic transfer. Ecotoxicol. Environ. Saf. 133, 188–194 (2016).
Article CAS PubMed Google Scholar
Chen, H.-Y., Grabner, D. S., Nachev, M., Shih, H.-H. & Sures, B. Effects of the acanthocephalan Polymorphus minutus and the microsporidian Dictyocoela duebenum on energy reserves and stress response of cadmium exposed Gammarus fossarum. PeerJ 3, e1353. https://doi.org/10.7717/peerj.1353 (2015).
Article CAS PubMed PubMed Central Google Scholar
Müller, J. Mitochondrial DNA variation and the evolutionary history of cryptic Gammarus fossarum types. Mol. Phylogenet. Evol. 15, 260–268 (2000).
Article PubMed CAS Google Scholar
Weiss, M., Macher, J. N., Seefeldt, M. A. & Leese, F. Molecular evidence for further overlooked species within the Gammarus fossarum complex (Crustacea: Amphipoda). Hydrobiologia 721, 165–184 (2014).
Article CAS Google Scholar
Weiss, M. & Leese, F. Widely distributed and regionally isolated! Drivers of genetic structure in Gammarus fossarum in a human-impacted landscape. BMC Evol. Biol. 16, 153 (2016).
Article PubMed PubMed Central CAS Google Scholar
Lagrue, C. et al. Confrontation of cryptic diversity and mate discrimination within Gammarus pulex and Gammarus fossarum species complexes. Freshw. Biol. 59, 2555–2570 (2014).
Article Google Scholar
Copilaş-Ciocianu, D. & Petrusek, A. The southwestern Carpathians as an ancient centre of diversity of freshwater gammarid amphipods: insights from the Gammarus fossarum species complex. Mol. Ecol. 24, 3980–3992 (2015).
Article PubMed Google Scholar
Copilaş-Ciocianu, D., Rutová, T., Pařil, P. & Petrusek, A. Epigean gammarids survived millions of years of severe climatic fluctuations in high latitude refugia throughout the Western Carpathians. Mol. Phylogenet. Evol. 112, 218–229 (2017).
Article PubMed Google Scholar
Hebert, P. D. N., Ratnasingham, S. & de Waard, J. R. Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc. R. Soc. Lond. B Biol. Sci. 270, S96–S99 (2003).
CAS Google Scholar
Kekkonen, M. & Hebert, P. D. N. DNA barcode-based delineation of putative species: efficient start for taxonomic workflows. Mol. Ecol. Resour. 14, 706–715 (2014).
Article PubMed PubMed Central Google Scholar
Copilaş-Ciocianu, D. & Petrusek, A. Phylogeography of a freshwater crustacean species complex reflects a long-gone archipelago. J. Biogeogr. 44, 421–432 (2017).
Article Google Scholar
Grabner, D. S. et al. Invaders, natives and their enemies: distribution patterns of amphipods and their microsporidian parasites in the Ruhr Metropolis, Germany. Parasit. Vectors 8, 419 (2015).
Article PubMed PubMed Central Google Scholar
Copilaș-Ciocianu, D., Zimța, A. & Petrusek, A. Integrative taxonomy reveals a new Gammarus species (Crustacea, Amphipoda) surviving in a previously unknown southeast European glacial refugium. J. Zool. Syst. Evol. Res. 57, 272–297 (2019).
Article Google Scholar
Copilaş-Ciocianu, D., Zimţa, A.-A., Grabowski, M. & Petrusek, A. Survival in northern microrefugia in an endemic Carpathian gammarid (Crustacea: Amphipoda). Zool. Scr. 47, 357–372 (2018).
Article Google Scholar
Hou, Z., Li, J. & Li, S. Diversification of low dispersal crustaceans through mountain uplift: a case study of Gammarus (Amphipoda: Gammaridae) with descriptions of four novel species: Diversification of Gammarus species. Zool. J. Linn. Soc. 170, 591–633 (2014).
Article Google Scholar
Wellborn, G. A. & Cothran, R. D. Niche diversity in crustacean cryptic species: complementarity in spatial distribution and predation risk. Oecologia 154, 175–183 (2007).
Article ADS PubMed Google Scholar
Delić, T., Švara, V., Coleman, C. O., Trontelj, P. & Fišer, C. The giant cryptic amphipod species of the subterranean genus Niphargus (Crustacea, Amphipoda). Zool. Scr. 46, 740–752 (2017).
Article Google Scholar
McInerney, C. E. et al. The ancient Britons: groundwater fauna survived extreme climate change over tens of millions of years across NW Europe. Mol. Ecol. 23, 1153–1166 (2014).
Article PubMed Google Scholar
Cooper, S. J. B. et al. Subterranean archipelago in the Australian arid zone: mitochondrial DNA phylogeography of amphipods from central Western Australia: Phylogeography of subterranean amphipods. Mol. Ecol. 16, 1533–1544 (2007).
Article CAS PubMed Google Scholar
Ethridge, J. Z., Gibson, J. R. & Nice, C. C. Cryptic diversity within and amongst spring-associated Stygobromus amphipods (Amphipoda: Crangonyctidae): Stygobromus amphipod cryptic diversity. Zool. J. Linn. Soc. 167, 227–242 (2013).
Article Google Scholar
Beermann, J. et al. Cryptic species in a well-known habitat: applying taxonomics to the amphipod genus Epimeria (Crustacea, Peracarida). Sci. Rep. 8, 6893. https://doi.org/10.1038/s41598-018-25225-x (2018).
Article ADS CAS PubMed PubMed Central Google Scholar
Baird, H. P., Miller, K. J. & Stark, J. S. Evidence of hidden biodiversity, ongoing speciation and diverse patterns of genetic structure in giant Antarctic amphipods: Genetic diversity in antarticamphipods. Mol. Ecol. 20, 3439–3454 (2011).
Article PubMed Google Scholar
Havermans, C. et al. Genetic and morphological divergences in the cosmopolitan deep-sea amphipod Eurythenes gryllus reveal a diverse abyss and a bipolar species. PLoS ONE 8, e74218. https://doi.org/10.1371/journal.pone.0074218 (2013).
Article ADS CAS PubMed PubMed Central Google Scholar
Havermans, C. Have we so far only seen the tip of the iceberg? Exploring species diversity and distribution of the giant amphipod Eurythenes. Biodiversity 17, 12–25 (2016).
Article Google Scholar
Coleman, C. O. Taxonomy in times of the taxonomic impediment—examples from the community of experts on amphipod crustaceans. J. Crustac. Biol. 35, 729–740 (2015).
Article Google Scholar
Arfianti, T., Wilson, S. & Costello, M. J. Progress in the discovery of amphipod crustaceans. PeerJ 6, e5187 (2018).
Article PubMed PubMed Central Google Scholar
Popov, S. V. et al. Lithological-Paleogeographic maps of Partethys—10 maps late Eocene to Plioccene. Courier Foschung-Intitut Senckberg 1–46 (2004).
Frenzel, B., Pécsi, M. & Velichko, A. A. Atlas of Paleoclimates and Paleoenvironments of the Northern Hemisphere: Late Pleistocene, Holocene (Geographical Research Institute Hungarian Academy of Science; G. Fischer, Budapest, 1992).
Google Scholar
Macneil, C., Dick, J. T. A. & Elwood, R. W. The trophic ecology of freshwater Gammarus spp (Crustacea; Amphipoda): Problems and perspective concerning the functional feeding group concept. Biol. Rev. 72, 349–364 (2007).
Article Google Scholar
Rudolph, K., Coleman, C. O., Mamos, T. & Grabowski, M. Description and post-glacial demography of Gammarusjazdzewskii sp. nov. (Crustacea: Amphipoda) from Central Europe. Syst. Biodivers. 16, 587–603 (2018).
Article Google Scholar
Westram, A. M., Jokela, J., Baumgartner, C. & Keller, I. Spatial distribution of cryptic species diversity in European freshwater amphipods (Gammarus fossarum) as revealed by pyrosequencing. PLoS ONE 6, e23879. https://doi.org/10.1371/journal.pone.0023879 (2011).
Article ADS CAS PubMed PubMed Central Google Scholar
Lefébure, T. et al. Phylogeography of a subterranean amphipod reveals cryptic diversity and dynamic evolution in extreme environments: cryptic and dynamic evolution in subsurface. Mol. Ecol. 15, 1797–1806 (2006).
Article PubMed CAS Google Scholar
Lefébure, T., Douady, C. J., Malard, F. & Gibert, J. Testing dispersal and cryptic diversity in a widely distributed groundwater amphipod (Niphargus rhenorhodanensis). Mol. Phylogenet. Evol. 42, 676–686 (2007).
Article PubMed CAS Google Scholar
Raffard, A., Santoul, F., Cucherousset, J. & Blanchet, S. The community and ecosystem consequences of intraspecific diversity: a meta-analysis: the ecological effects of intraspecific diversity. Biol. Rev. 94, 648–661 (2019).
Article PubMed Google Scholar
Jabłońska, A., Wrzesińska, W., Zawal, A., Pešić, V. & Grabowski, M. Long-term within-basin isolation patterns, different conservation units, and interspecific mitochondrial DNA introgression in an amphipod endemic to the ancient Lake Skadar system, Balkan Peninsula. Freshw. Biol. 65, 209–225 (2020).
Article CAS Google Scholar
Stein, E. D., Martinez, M. C., Stiles, S., Miller, P. E. & Zakharov, E. V. Is DNA barcoding actually cheaper and faster than traditional morphological methods: results from a survey of freshwater bioassessment efforts in the United States?. PLoS ONE 9, e95525. https://doi.org/10.1371/journal.pone.0095525 (2014).
Article ADS CAS PubMed PubMed Central Google Scholar
Valentini, A., Pompanon, F. & Taberlet, P. DNA barcoding for ecologists. Trends Ecol. Evol. 24, 110–117 (2009).
Article PubMed Google Scholar
Weigand, H. et al. DNA barcode reference libraries for the monitoring of aquatic biota in Europe: Gap-analysis and recommendations for future work. Sci. Total Environ. 678, 499–524 (2019).
Article ADS CAS PubMed Google Scholar
Delić, T., Trontelj, P., Rendoš, M. & Fišer, C. The importance of naming cryptic species and the conservation of endemic subterranean amphipods. Sci. Rep. 7, 3391. https://doi.org/10.1038/s41598-017-02938-z (2017).
Article ADS CAS PubMed PubMed Central Google Scholar
Wysocka, A. et al. Origin of the Lake Ohrid gammarid species flock: ancient local phylogenetic lineage diversification. J. Biogeogr. 41, 1758–1768 (2014).
Article Google Scholar
Esmaeili-Rineh, S., Sari, A., Delić, T., Moškrič, A. & Fišer, C. Molecular phylogeny of the subterranean genus Niphargus (Crustacea: Amphipoda) in the Middle East: a comparison with European Niphargids. Zool. J. Linn. Soc. 175, 812–826 (2015).
Article Google Scholar
Hou, Z., Sket, B., Fiser, C. & Li, S. Eocene habitat shift from saline to freshwater promoted Tethyan amphipod diversification. Proc. Natl. Acad. Sci. 108, 14533–14538 (2011).
Article ADS CAS PubMed PubMed Central Google Scholar
AltschuP, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic Local Alignment Search Tool. 8.
Karaman, G. & Pinkster, S. Freshwater Gammarus species from Europe, North Africa and adjacent regions of Asia (Crustacea-Amphipoda). Part I. Gammarus pulex-group and related species. Bijdragen tot de Dierkunde 1–97 (1977).
Piscart, C. & Bollache, L. Crustacés amphipodes de surface : gammares d’eau douce (Association Française de Limnologie, Thonon-les-Bains, 2012).
Google Scholar
Eggers, T. & Martens, A. A key of the freshwater Amphipods of Germany. Lauterbornia 1–68 (2001).
Jazdzewski, K. Morfologia, taksonomia i wystepowanie w Polsce kielzy z rodzajów Gammarus Fabr. i Chaetogammarus Mart. (Crustacea, Amphipoda). Acta Universitatis Lodziensis 1–187 (1975).
Karaman, G. S. Crustacea (Amphipoda di acqua dolce. Calderini, Bologna, 1993).
Google Scholar
Hillis, D. M., Moritz, C. & Mable, B. K. Molecular Systematics (Sinauer Associates Inc, Sunderland, 1996).
Google Scholar
Casquet, J., Thebaud, C. & Gillespie, R. G. Chelex without boiling, a rapid and easy technique to obtain stable amplifiable DNA from small amounts of ethanol-stored spiders: Quick and easy dsDNA extraction for barcoding. Mol. Ecol. Resour. 12, 136–141 (2012).
Article CAS PubMed Google Scholar
Folmer, O., Black, M., Hoeh, W., Lutz, R. & Vrijenhoek, R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol. Mar. Biol. Biotechnol 6, 294–299 (1994).
Google Scholar
Costa, F. O., Henzler, C. M., Lunt, D. H., Whiteley, N. M. & Rock, J. Probing marine Gammarus (Amphipoda) taxonomy with DNA barcodes. Syst. Biodivers. 7, 365–379 (2009).
Article Google Scholar
Astrin, J. J. & Stüben, P. E. Phylogeny in cryptic weevils: molecules, morphology and new genera of western Palaearctic Cryptorhynchinae (Coleoptera:Curculionidae). Invertebr. Syst. 22, 503–522 (2008).
Article Google Scholar
Larkin, M. A. et al. Clustal W and Clustal X version 2.0. Bioinformatics 23, 2947–2948 (2007).
Article CAS PubMed Google Scholar
Kumar, S., Stecher, G. & Tamura, K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874 (2016).
Article CAS PubMed PubMed Central Google Scholar
Librado, P. & Rozas, J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25, 1451–1452 (2009).
Article CAS PubMed Google Scholar
Puillandre, N., Lambert, A., Brouillet, S. & Achaz, G. ABGD, Automatic Barcode Gap Discovery for primary species delimitation. Mol. Ecol. 21, 1864–1877 (2012).
Article CAS PubMed Google Scholar
Ratnasingham, S. & Hebert, P. D. N. A DNA-based registry for all animal species: the barcode index number (BIN) System. PLoS ONE 8, e66213. https://doi.org/10.1371/journal.pone.0066213 (2013).
Article ADS CAS PubMed PubMed Central Google Scholar
Zhang, J., Kapli, P., Pavlidis, P. & Stamatakis, A. A general species delimitation method with applications to phylogenetic placements. Bioinformatics 29, 2869–2876 (2013).
Article CAS PubMed PubMed Central Google Scholar
Puillandre, N. et al. Large-scale species delimitation method for hyperdiverse groups: large-scale species delimitation. Mol. Ecol. 21, 2671–2691 (2012).
Article CAS PubMed Google Scholar
Ratnasingham, S. & Hebert, P. D. N. BOLD: The Barcode of life data system (https://www.barcodinglife.org). Mol. Ecol. Notes7, 355–364 (2007).
Stamatakis, A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313 (2014).
Article CAS PubMed PubMed Central Google Scholar
Bouckaert, R. et al. BEAST 2: a software platform for Bayesian evolutionary analysis. PLoS Comput. Biol. 10, e1003537. https://doi.org/10.1371/journal.pcbi.1003537 (2014).
Article CAS PubMed PubMed Central Google Scholar
Bouckaert, R. R. & Drummond, A. J. bModelTest: Bayesian phylogenetic site model averaging and model comparison. BMC Evol. Biol. 17, 42 (2017).
Article PubMed PubMed Central Google Scholar
Baele, G., Lemey, P. & Vansteelandt, S. Make the most of your samples: Bayes factor estimators for high-dimensional models of sequence evolution. BMC Bioinform. 14, 85 (2013).
Article Google Scholar
Brower, A. V. Rapid morphological radiation and convergence among races of the butterfly Heliconius erato inferred from patterns of mitochondrial DNA evolution. Proc. Natl. Acad. Sci. 91, 6491–6495 (1994).
Article ADS CAS PubMed PubMed Central Google Scholar
Rambaut, A., Drummond, A. J., Xie, D., Baele, G. & Suchard, M. A. Posterior summarization in Bayesian phylogenetics using Tracer 1.7. Syst. Biol. 67, 901–904 (2018).
Article CAS PubMed PubMed Central Google Scholar
Felsenstein, J. Confidence limits on phylogenies: an approach using the boostrap. Evolution 39, 783–791 (1985).
Article PubMed Google Scholar
Miller, M. P. Alleles In Space (AIS): computer software for the joint analysis of interindividual spatial and genetic information. J. Hered. 96, 722–724 (2005).
Article CAS PubMed Google Scholar

Download references

Acknowledgements

The authors thank the French National Agency for Research for funding (ANR-13-BSV7-0004 ANR MULTISTRESS; ANR-14-CE21-0006 PROTEOGAM), the Adour-Garonne water agency (program Dylaq) and the French Office for Biodiversity (OFB) for funding. Study was financially supported by the Polish Ministry of Science and Education, grant number N N303 579439 and partially from the internal funds of the University of Lodz. Additional financing for MJ’s stay in Poland was provided by the Short-Term Scientific Mission (STSM) of COST Action CA15219-Developing new genetic tools for bioassessment of aquatic ecosystems in Europe—DNAqua-Net (Grant no. 39779)".

Author information

Authors and Affiliations

UMR CNRS 6282 Biogéosciences, Université Bourgogne Franche Comté, Dijon, France
Remi Wattier, Anthony Ollivier & Thierry Rigaud
Department of Invertebrate Zoology and Hydrobiology, University of Lodz, Lodz, Poland
Tomasz Mamos, Tomasz Rewicz & Michał Grabowski
Zoological Institute, University of Basel, Basel, Switzerland
Tomasz Mamos
Institute of Ecology, Nature Research Centre, Vilnius Nature Research Centre, Institute of Ecology, Vilnius, Lithuania
Denis Copilaş-Ciocianu
Department of Natural Sciences, Varaždin City Museum, Varaždin, Croatia
Mišel Jelić
Laboratoire d’écotoxicologie, INRAE, UR RiverLy, Villeurbanne, France
Arnaud Chaumot
UMR CNRS 73602 LIEC, Université de Lorraine, Metz, France
Michael Danger & Vincent Felten
Univ Rennes, CNRS, ECOBIO-UMR 6553, 35000, Rennes, France
Christophe Piscart
Department of Teacher Education Studies in Gospić, University of Zadar, Gospić, Croatia
Krešimir Žganec
University of Guelph, Centre for Biodiversity Genomics, Guelph, ON, Canada
Tomasz Rewicz
Department of Genetics and Biosystematics, University of Gdansk, Gdansk, Poland
Anna Wysocka

Authors

Remi Wattier
View author publications
You can also search for this author in PubMed Google Scholar
Tomasz Mamos
View author publications
You can also search for this author in PubMed Google Scholar
Denis Copilaş-Ciocianu
View author publications
You can also search for this author in PubMed Google Scholar
Mišel Jelić
View author publications
You can also search for this author in PubMed Google Scholar
Anthony Ollivier
View author publications
You can also search for this author in PubMed Google Scholar
Arnaud Chaumot
View author publications
You can also search for this author in PubMed Google Scholar
Michael Danger
View author publications
You can also search for this author in PubMed Google Scholar
Vincent Felten
View author publications
You can also search for this author in PubMed Google Scholar
Christophe Piscart
View author publications
You can also search for this author in PubMed Google Scholar
Krešimir Žganec
View author publications
You can also search for this author in PubMed Google Scholar
Tomasz Rewicz
View author publications
You can also search for this author in PubMed Google Scholar
Anna Wysocka
View author publications
You can also search for this author in PubMed Google Scholar
Thierry Rigaud
View author publications
You can also search for this author in PubMed Google Scholar
Michał Grabowski
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

Designed the study: R.W. and M.G. Sampling: R.W., M.G., V.F., C.P., A.C., M.D., M.J., K.Z., T.M. and D.C.-C. Performed the laboratory work: R.W., A.O., M.J., D.C.-C., T.M., T.Re. and A.W. Contributed reagents: R.W., T.Ri and M.G. Data analysis: R.W., T.M., D.C.-C., M.G., and M.J. Discussion and interpretation of the results: R.W., T.M., M.G., D.C.-C., T. Ri and M.J. Wrote the manuscript: R.W. and M.G. All authors have read and approved the final version of the manuscript.

Corresponding authors

Correspondence to Remi Wattier or Michał Grabowski.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher's note

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

Supplementary information

Supplementary Figure S1.

Supplementary Captions.

Supplementary Table S1.

Supplementary Table S2.

Rights and permissions

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

Reprints and permissions

About this article

Cite this article

Wattier, R., Mamos, T., Copilaş-Ciocianu, D. et al. Continental-scale patterns of hyper-cryptic diversity within the freshwater model taxon Gammarus fossarum (Crustacea, Amphipoda). Sci Rep 10, 16536 (2020). https://doi.org/10.1038/s41598-020-73739-0

Download citation

Received: 15 May 2020
Accepted: 16 September 2020
Published: 06 October 2020
DOI: https://doi.org/10.1038/s41598-020-73739-0

This article is cited by

Multiple genetic species in a halophilic non-marine ostracod (Crustacea)
- Adrianna Kilikowska
- Isa Schön
- Tadeusz Namiotko
Hydrobiologia (2024)
Evidence for a single population expansion event across 24,000 km: the case of the deep-sea scavenging amphipod Abyssorchomene distinctus
- Deborah W. E. Dupont
- Tasnim Patel
- Isa Schön
Hydrobiologia (2024)
Organ-oriented proteogenomics functional atlas of three aquatic invertebrate sentinel species
- Maxime Leprêtre
- Davide Degli Esposti
- Arnaud Chaumot
Scientific Data (2023)
Conservation of freshwater biodiversity in Tunisia in a climate change context: combining amphipod distribution data and molecular analyses to improve priorities
- Khaoula Ayati
- Kamil Hupało
- Christophe Piscart
Biodiversity and Conservation (2023)
DNA barcode-based survey documents underestimated diversity and intricate phylogeographic patterns of aquatic Heteroptera in an endangered Balkan biodiversity hotspot: ancient Lake Skadar basin
- Tomasz Rewicz
- Grzegorz Tończyk
- Michał Grabowski
Biodiversity and Conservation (2023)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.