Spatial distribution of cryptic species diversity in european freshwater amphipods (Gammarus fossarum) as revealed by pyrosequencing

Diversity and distribution of air-breathing sea slug genus Peronia Fleming, 1822 (Gastropoda: Onchidiidae) in southern Japanese waters

Species of the genus Peronia Fleming, 1822, are air-breathing onchidiid sea slugs that inhabit intertidal reef flats of temperate to tropical zones. In the Ryukyu Islands of southern subtropical Japan, Peronia species are a traditional food source for local people. To date, there have been three species recorded around Okinawajima Island; P. verruculata and P. peronii, along with recently described P. okinawensis, which was described as possibly endemic to Okinawajima Island. This study aimed to map the distribution ranges of these three Peronia species within the Ryukyu Islands using molecular analyses in order to understand the specific distribution of each species. Since Peronia species are generally indistinguishable by gross external morphology, a DNA barcoding approach was employed to identify specimens. The molecular data showed that there are four species present in the Ryukyu Islands. P. verruculata (unit #1 sensu Dayrat et al., 2020) was dominant at almost all locations, while P. peronii was present in much lower numbers than P. verruculata, but found across a relatively wide range in the Ryukyu Islands. We newly record P. okinawensis and P. setoensis from Amami Oshima Island and from several places around Okinawajima Island, and also identified high levels of genetic variation within P. setoensis. Peronia okinawensis and P. setoensis have been thought to be endemic to Okinawajima Island and to Honshu, mainland Japan, respectively. However, as both species were observed around Okinawajima and Amami Oshima islands, other islands of the Ryukyus are also likely to harbor these species, and their distribution ranges are wider than previously thought. Based on the results from molecular analyses, we provide general descriptions of each species. Sizes of specimens were consistently smaller for P. setoensis and relatively larger for P. peronii specimens. On the other hand, P. verruculata and P. okinawensis had similar size ranges, but P. okinawensis had comparatively much more distinct papillae. This study revealed that the Ryukyu Islands are the only region currently known with four sympatric Peronia species, and this work provides a basis for future research on these Peronia species throughout the northwest Pacific Ocean, representing the first step in more effective management of the local Peronia fisheries in the Ryukyu Islands.

Temporal Patterns of Bacterial and Viral Communities during Algae Blooms of a Reservoir in Macau

Compositions of microbial communities associated with blooms of algae in a storage reservoir in Macau, China were investigated between 2013 and 2016. Algae were enumerated by visible light microscopy. Profiles of organisms in water were examined by 16S rRNA sequences and viral metagenomics, based on next generation sequencing. Results of 16S rRNA sequencing indicated that majority of the identified organisms were bacteria closely related to Proteobacteria, Cyanobacteria, Verrucomicrobia, Bacteroidetes, and Actinobacteria. Metagenomics sequences demonstrated that the dominant virus was Phycodnavirus, accounting for 70% of the total population. Patterns of relative numbers of bacteria in the microbial community and their temporal changes were determined through alpha diversity indices, principal coordinates analysis (PCoA), relative abundance, and visualized by Venn diagrams. Ways in which the bacterial and viral communities are influenced by various water-related variables were elucidated based on redundancy analysis (RDA). Relationships of the relative numbers of bacteria with trophic status in a reservoir used for drinking water in Macau, provided insight into associations of Phycodnavirus and Proteobacteria with changes in blooms of algae.

High diversity and local endemism in Aotearoa New Zealand's groundwater crustacean fauna

We used DNA barcoding to assess the diversity and distribution of New Zealand's groundwater amphipods and isopods (Crustacea) and to determine whether biodiversity and endemism within tectonically active New Zealand are similar to those of more tectonically stable continents. Sixty-five wells were sampled in seven aquifers across four regions within the North and South islands of New Zealand, and resident invertebrates were morphologically identified and then assessed using sequencing of the mitochondrial DNA cytochrome c oxidase subunit one (COI) gene. Invertebrates were found in 54 wells. Of the 228 individual amphipods and isopods found in 36 of the wells, 154 individuals were successfully sequenced for COI (68% success rate) from 25 wells, with at least one well in each aquifer containing sequenced individuals. Of the 45 putative species identified using Barcode Index Numbers (BINs), 30 BINs (78% of all taxa and 83% of amphipods) were previously unrecorded. Substantial morphologically cryptic, species-level diversity was revealed, particularly within the amphipod Family Paraleptamphopidae. Similarly, one isopod taxon morphologically identified as Cruregens fontanus was assigned to five well-separated BINs based on COI sequences. Endemism appeared high, with all taxa regionally endemic; 87% of species were restricted to one aquifer and more than 50% restricted to one well. Non-saturated species accumulation curves indicated that, while additional sampling may increase the range of some currently identified taxa, additional range-restricted taxa are also likely to be discovered. Patterns of diversity and short-range endemism were similar to those found elsewhere, including locations which are more tectonically stable. The predominance of local endemism within New Zealand's groundwater fauna suggests that land-use activities and groundwater extraction require careful evaluation to minimize threats to groundwater biodiversity.

Dispersal behaviour and riverine network connectivity shape the genetic diversity of freshwater amphipod metapopulations

Theory predicts that the distribution of genetic diversity in a landscape is strongly dependent on the connectivity of the metapopulation and the dispersal of individuals between patches. However, the influence of explicit spatial configurations such as dendritic landscapes on the genetic diversity of metapopulations is still understudied, and theoretical corroborations of empirical patterns are largely lacking. Here, we used microsatellite data and stochastic simulations of two metapopulations of freshwater amphipods in a 28,000 km² riverine network to study the influence of spatial connectivity and dispersal strategies on the spatial distribution of their genetic diversity. We found a significant imprint of the effects of riverine network connectivity on the local and global genetic diversity of both amphipod species. Data from 95 sites showed that allelic richness significantly increased towards more central nodes of the network. This was also seen for observed heterozygosity, yet not for expected heterozygosity. Genetic differentiation increased with instream distance. In simulation models, depending on the mutational model assumed, upstream movement probability and dispersal rate, respectively, emerged as key factors explaining the empirically observed distribution of local genetic diversity and genetic differentiation. Surprisingly, the role of site‐specific carrying capacities, for example by assuming a direct dependency of population size on local river size, was less clear cut: while our best fitting model scenario included this feature, over all simulations, scaling of carrying capacities did not increase data‐model fit. This highlights the importance of dispersal behaviour along spatial networks in shaping population genetic diversity.

Dictyocoela microsporidia diversity and co-diversification with their host, a gammarid species complex (Crustacea, Amphipoda) with an old history of divergence and high endemic diversity

BACKGROUND

Although the processes of co-evolution between parasites and their hosts are well known, evidence of co-speciation remains scarce. Microsporidian intracellular parasites, due to intimate relationships with their hosts and mixed mode of transmission (horizontal but also vertical, from mother to offspring), may represent an interesting biological model for investigating co-speciation. Amphipod crustaceans, especially gammarids, are regular hosts of microsporidian parasites, in particular the Dictyocoela spp., which have so far been found limited to these amphipods and are known to use a vertical mode of transmission. The amphipod genus Gammarus has a diversification history spanning the last 50-60 Mya and an extensive cryptic diversity in most of the nominal species. Here, we investigated the degree of co-diversification between Dictyocoela and Gammarus balcanicus, an amphipod with high degrees of ancient cryptic diversification and lineage endemism, by examining the genetic diversity of these parasites over the entire geographic range of the host. We hypothesised that the strong host diversification and vertical transmission of Dictyocoela would promote co-diversification.

RESULTS

Using the parasite SSU rDNA as a molecular marker, analyzing 2225 host specimens from 88 sites covering whole host range, we found 31 haplogroups of Dictyocoela, 30 of which were novel, belonging to four Dictyocoela species already known to infect other Gammarus spp. The relationships between Dictyocoela and gammarids is therefore ancient, with the speciation in parasites preceding those of the hosts. Each novel haplogroup was nevertheless specific to G. balcanicus, leaving the possibility for subsequent co-diversification process during host diversification. A Procrustean Approach to Co-phylogeny (PACo) analysis revealed that diversification of Dictyocoela was not random with respect to that of the host. We found high degrees of congruence between the diversification of G. balcanicus and that of Dictyocoela roeselum and D. muelleri. However, we also found some incongruences between host and Dictyocoela phylogenies, e.g. in D. duebenum, probably due to host shifts between different G. balcanicus cryptic lineages.

CONCLUSION

The evolutionary history of Dictyocoela and Gammarus balcanicus represents an example of an overall host-parasite co-diversification, including cases of host shifts.

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

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.

How to quantify the links between bioavailable contamination in watercourses and pressures of anthropogenic land cover, contamination sources and hydromorphology at multiple scales?

Active biomonitoring permits the quantification of biological exposure to chemicals through measurements of bioavailable concentrations in biota and biological markers of toxicity in organisms. It enables respective comparison of the levels of contamination between sites and sampling campaigns. Caged gammarids are recently proposed as relevant probes for measuring bioavailable contamination in freshwater systems. The purpose of the present study was to develop a multi-pressure and multiscale approach, considering metallic contamination levels (from data based on active biomonitoring) as a response to pressures (combination of individual stressors). These pressures were anthropogenic land cover, industry density, wastewater treatment plant density, pressures on stream hydromorphological functioning, riverside vegetation and bioavailability factors. A dataset combining active biomonitoring and potentially related pressures was established at the French national scale, with 196 samplings from 2009 to 2016. The links between pressures and metallic contamination were defined and modelled via structural equation modeling (more specifically partial least squares - path modeling). The model enabled the understanding of the respective influences of pressures on metallic bioconcentration in caged sentinel organisms. Beyond validating the local influence of industries and wastewater treatment plants on metallic contamination, this model showed a complementary effect of driving forces of anthropogenic land cover (leading to human activities). It also quantified a significant influence of pressures on stream hydromorphological functioning, presence of vegetation and physico-chemical parameters on metal bioconcentration. This hierarchical multi-pressure approach could serve as a concept on how pressures and contamination (assessed by active biomonitoring) can be connected. Its future application will enable better understanding of environmental pressures leading to contamination in freshwater ecosystems.

Origin of the natural variation in the storage of dietary carotenoids in freshwater amphipod crustaceans

Carotenoids are diverse lipophilic natural pigments which are stored in variable amounts by animals. Given the multiple biological functions of carotenoids, such variation may have strong implications in evolutionary biology. Crustaceans such as Gammarus amphipods store large amounts of these pigments and inter-population variation occurs. While differences in parasite selective pressure have been proposed to explain this variation, the contribution of other factors such as genetic differences in the gammarid ability to assimilate and/or store pigments, and the environmental availability of carotenoids cannot be dismissed. This study investigates the relative contributions of the gammarid genotype and of the environmental availability of carotenoids in the natural variability in carotenoid storage. It further explores the link of this natural variability in carotenoid storage with major crustacean immune parameters. We addressed these aspects using the cryptic diversity in the amphipod crustacean Gammarus fossarum and a diet supplementation protocol in the laboratory. Our results suggest that natural variation in G. fossarum storage of dietary carotenoids results from both the availability of the pigments in the environment and the genetically-based ability of the gammarids to assimilate and/or store them, which is associated to levels of stimulation of cellular immune defences. While our results may support the hypothesis that carotenoids storage in this crustacean may evolve in response to parasitic pressure, a better understanding of the specific roles of this large pigment storage in the crustacean physiology is needed.

Complete transcriptome assembly and annotation of a critically important amphipod species in freshwater ecotoxicological risk assessment: Gammarus fossarum

Because of their crucial role in ecotoxicological risk assessment, amphipods (Crustacea) are commonly employed as model species in a wide range of studies. However, despite their ecological importance, their genome has not yet been completely annotated and molecular mechanisms underlying key pathways, such as the serotonin pathway, in development of ecotoxicological biomarkers of exposure to neuroactive pharmaceuticals are still poorly understood. Furthermore, genetic similarities and discrepancies with other model arthropods (e.g., Drosophila melanogaster) have not been completely clarified. In this report, we present a new transcriptome assembly of Gammarus fossarum, an important amphipod species, widespread in Central Europe. RNA-Seq with Illumina HiSeq technology was used to analyse samples extracted from total internal tissues. We used the Trinity and Trinotate software suites for transcriptome assembly and annotation, respectively. The quality of this assembly and the affiliated targeted homology searches greatly enrich the molecular knowledge on this species. Because of the lack of publicly available molecular information on the serotonin pathway, we also highlighted sequence homologies and divergences of the genes encoding the serotonin pathway components of the well-annotated arthropod D. melanogaster, and Crustacea with the corresponding genes of our assembly. An inferior number of hits was found when running a BLAST analysis of both D. melanogaster and Crustacea mRNA sequences encoding serotonin receptors available in GenBank against the total assembly, compared to other serotonin pathway components. A lack of information on important components for serotonin biosynthesis and vesicle endocytosis (i.e., tryptophan hydroxylase and vesicular monoamine transporter) in Crustacea was also brought to light. Our results will provide an extensive transcriptional resource for this important species in ecotoxicological risk assessment and highlight the need for a more detailed categorization of neuronal pathways components in invertebrates.

Acanthocephalan infection patterns in amphipods: a reappraisal in the light of recently discovered host cryptic diversity

Amphipods are model species in studies of pervasive biological patterns such as sexual selection, size assortative pairing and parasite infection patterns. Cryptic diversity (i.e. morphologically identical but genetically divergent lineages) has recently been detected in several species. Potential effects of such hidden diversity on biological patterns remain unclear, but potentially significant, and beg the question of whether we have missed part of the picture by involuntarily overlooking the occurrence and effects of cryptic diversity on biological patterns documented by previous studies. Here we tested for potential effects of cryptic diversity on parasite infection patterns in amphipod populations and discuss the implications of our results in the context of previously documented host-parasite infection patterns, especially amphipod-acanthocephalan associations. We assessed infection levels (prevalence and abundance) of 3 acanthocephalan species (Pomphorhynchus laevis, P. tereticollis and Polymorphus minutus) among cryptic lineages of the Gammarus pulex/G. fossarum species complex and G. roeseli from sampling sites where they occur in sympatry. We also evaluated potential differences in parasite-induced mortality among host molecular operational taxonomic units (MOTUs)-parasite species combinations. Acanthocephalan prevalence, abundance and parasite-induced mortality varied widely among cryptic MOTUs and parasite species; infection patterns were more variable among MOTUs than sampling sites. Overall, cryptic diversity in amphipods strongly influenced apparent infection levels and parasite-induced mortality. Future research on species with cryptic diversity should account for potential effects on documented biological patterns. Results from previous studies may also need to be reassessed in light of cryptic diversity and its pervasive effects.

De novo transcriptomes of 14 gammarid individuals for proteogenomic analysis of seven taxonomic groups

Gammarids are amphipods found worldwide distributed in fresh and marine waters. They play an important role in aquatic ecosystems and are well established sentinel species in ecotoxicology. In this study, we sequenced the transcriptomes of a male individual and a female individual for seven different taxonomic groups belonging to the two genera Gammarus and Echinogammarus: Gammarus fossarum A, G. fossarum B, G. fossarum C, Gammarus wautieri, Gammarus pulex, Echinogammarus berilloni, and Echinogammarus marinus. These taxa were chosen to explore the molecular diversity of transcribed genes of genotyped individuals from these groups. Transcriptomes were de novo assembled and annotated. High-quality assembly was confirmed by BUSCO comparison against the Arthropod dataset. The 14 RNA-Seq-derived protein sequence databases proposed here will be a significant resource for proteogenomics studies of these ecotoxicologically relevant non-model organisms. These transcriptomes represent reliable reference sequences for whole-transcriptome and proteome studies on other gammarids, for primer design to clone specific genes or monitor their specific expression, and for analyses of molecular differences between gammarid species.

Measurement(s)	transcription profiling assay
Technology Type(s)	RNA sequencing
Factor Type(s)	sex • species
Sample Characteristic - Organism	Gammarus • Echinogammarus
Sample Characteristic - Environment	habitat

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.9777905

Multisubstance Indicators Based on Caged Gammarus Bioaccumulation Reveal the Influence of Chemical Contamination on Stream Macroinvertebrate Abundances across France

Most anthropogenic stressors affecting freshwater systems are qualitatively known. However, the quantitative assessment of contaminant exposure and effects to aquatic communities is still difficult, limiting the understanding of consequences on aquatic ecosystem functioning and the implementation of effective management plans. Here, multisubstance indicators based on caged gammarid bioaccumulated contamination data are proposed (for metals and persistent organic pollutants, POPs) to map the bioavailable contamination level of freshwater ecosystems at a large spatial scale. We assessed the ability of these indicators to highlight the relationships between chemical exposure gradients and alteration in the abundance of macroinvertebrate populations on a data set of 218 watercourses distributed throughout France. We identified spatial regional heterogeneities in the levels of bioavailable contamination of metals (18 compounds) and POPs (43 compounds). Besides this, a degradation of Gammaridae, Ephemeridae, and Hydrobiidae densities with increasing levels of metal contamination are identified relative to Baetidae, Chironomidae, and Hydropsychidae. We show here that active biomonitoring allows the establishment of multisubstance indicators of bioavailable contamination, which reliably quantify chemical exposure gradients in freshwater ecosystems. Our ability to identify species-specific responses to chemical exposure gradients demonstrates the promising possibility to further decipher the effects of chemical contamination on macroinvertebrate assemblages through this type of indicator.

Effects of the herbicide Roundup® on the metabolic activity of Gammarus fossarum Koch, 1836 (Crustacea; Amphipoda)

Pesticides can easily reach surface waters via runoff and their potential to have detrimental impacts on freshwater organisms is high. Not much is known about how macroinvertebrates react to glyphosate contamination. In this study we investigated lethal and sublethal effects of the exposure of Gammarus fossarum to Roundup®, a glyphosate-based herbicide. The LC10 and LC50 values after 96 h were determined to be 0.65 ml/L Roundup® (230 mg/L glyphosate) and 0.96 ml/L Roundup® (340 mg/L glyphosate), respectively. As a sublethal measure of toxicity we conducted eight experiments with the feeding activity and the respiratory electron transport system (ETS) activity as endpoints. All experiments lasted seven days. Although the LC10 concentration of Roundup® was used for the feeding activity tests, 49% of the gammarids died before the end of the experiments, which is inconsistent with the calculated LC10-values. The feeding activity was significantly higher in Roundup®-enriched water (mean = 0.18 mg/mg x d) in comparison to pure spring water (mean = 0.079 mg/mg x d). No significant difference was observed between the ETS activity, which was determined after 24, 48 or 96 h after the start of the experiment, of the gammarids in Roundup® solution and in the control. The LC-values determined here are rather high, and exceed background glyphosate concentrations in most anthropogenically influenced surface waters. The increased feeding activity when exposed to Roundup® in combination with an unchanged ETS activity suggests effects on the metabolic efficiency of G. fossarum. We argue that Roundup® enhances the anabolic activity (feeding activity) in order to maintain the catabolic activity (ETS activity).

Living apart-together: Microhabitat differentiation of cryptic nematode species in a saltmarsh habitat

Coexistence of highly similar species is at odds with ecological theory of competition; coexistence, then, requires stabilizing mechanisms such as differences in ecological niche. In the bacterivore nematode Litoditis marina species complex, which occurs associated with macro-algae, four cryptic lineages (Pm I-IV) co-occur in the field along the south-western coast and estuaries of The Netherlands. Here we investigate the temporal and/or spatial niche differentiation in their natural environment using a qPCR-based detection and relative quantification method. We collected different algal species (i.e. two Fucus species and Ulva sp.) and separated algal structures (i.e. receptacula, thalli, non-fertile tips and bladders) at different sampling months and times (i.e. twice per sampling month), to examine differences in microhabitat use between coexisting L. marina species. Results demonstrate that the cryptic species composition varied among different algal species and algal structures, which was also subject to temporal shifts. Pm I dominated on Fucus spp., Pm II showed dominance on Ulva sp., while Pm III overall had the lowest frequencies. Microhabitat partitioning was most pronounced between the two cryptic species which had similar microbiomes (Pm I and Pm II), and less so between the two species which had significantly different microbiomes (Pm I and Pm III), suggesting that species which share the same microhabitats may avoid competition through resource partitioning. The interplay of microhabitat differentiation and temporal dynamics among the cryptic species of L. marina implies that there is a complex interaction between biotic components and abiotic factors which contributes to their coexistence in the field.

UV-irradiation and leaching in water reduce the toxicity of imidacloprid-contaminated leaves to the aquatic leaf-shredding amphipod Gammarus fossarum

Systemic neonicotinoid insecticides such as imidacloprid are increasingly applied against insect pest infestations on forest trees. However, leaves falling from treated trees may reach nearby surface waters and potentially represent a neonicotinoid exposure source for aquatic invertebrates. Given imidacloprid's susceptibility towards photolysis and high water solubility, it was hypothesized that the leaves' toxicity might be modulated by UV-irradiation during decay on the forest floor, or by leaching and re-mobilization of the insecticide from leaves within the aquatic ecosystem. To test these hypotheses, the amphipod shredder Gammarus fossarum was fed (over 7 d; n = 30) with imidacloprid-contaminated black alder (Alnus glutinosa) leaves that had either been pre-treated (i.e., leached) in water for up to 7 d or UV-irradiated for 1 d (at intensities relevant during autumn in Central Europe) followed by a leaching duration of 1 d. Gammarids' feeding rate, serving as sublethal response variable, was reduced by up to 80% when consuming non-pretreated imidacloprid-contaminated leaves compared to imidacloprid-free leaves. Moreover, both leaching of imidacloprid from leaves (for 7 d) as well as UV-irradiation reduced the leaves' imidacloprid load (by 46 and 90%) thereby mitigating the effects on gammarids' feeding rate to levels comparable to the respective imidacloprid-free controls. Therefore, natural processes, such as UV-irradiation and re-mobilization of foliar insecticide residues in water, might be considered when evaluating the risks systemic insecticide applications in forests might pose for aquatic organisms in nearby streams.

Epigean gammarids survived millions of years of severe climatic fluctuations in high latitude refugia throughout the Western Carpathians

Isolated glacial refugia have been documented in Central Europe for a number of taxa, but conclusive evidence for epigean aquatic species has remained elusive. Using molecular data (mitochondrial and nuclear markers), we compared the spatial patterns of lineage diversity of the widely distributed Gammarus fossarum species complex between two adjacent biogeographically and geomorphologically distinct Central European regions: the Bohemian Massif and the Western Carpathians. We investigated if the observed patterns of spatial diversity are more likely to stem from historical or present-day factors. Phylogenetic and phylogeographic analyses revealed eight phylogenetically diverse lineages: two exhibiting local signatures of recent demographic expansion inhabit both regions, while the other six display a relict distributional pattern and are found only in the Western Carpathians. Molecular dating indicates that these lineages are old and probably diverged throughout the Miocene (7-18Ma). Furthermore, their distribution does not seem to be constrained by the present boundaries of river catchments or topography. The contrasting spatial patterns of diversity observed between the two regions thus more likely result from historical rather than contemporaneous or recent factors. Our results indicate that despite the high latitude and proximity to the Pleistocene ice sheets, the Western Carpathians functioned as long-term glacial refugia for permanent freshwater fauna, allowing the uninterrupted survival of ancient lineages through millions of years of drastic climatic fluctuations.

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?

Shared parasites can strongly influence the outcome of competition between congeneric, sympatric hosts, and thus host population dynamics. Parasite-mediated competition is commonly hypothesized as an important factor in biological invasion success; invasive species often experience lower infection levels and/or parasite-induced mortality than native congeneric hosts. However, variation in infection levels among sympatric hosts can be due to contrasting abilities to avoid infection or different parasite-induced mortality rates following infection. Low parasite infection levels in a specific host can be due to either factor but have drastically different implications in interaction outcomes between sympatric hosts.

We assessed acanthocephalan infection levels (prevalence and abundance) among cryptic molecular taxonomic units (MOTU) of the native G. pulex/G. fossarum species complex from multiple populations where they occur in sympatry. We concomitantly estimated the same parameters in the invasive Gammarus roeseli commonly found in sympatry with G. pulex/G. fossarum MOTUs. We then tested for potential differences in parasite-induced mortality among these alternative hosts. As expected, the invasive G. roeseli showed relatively low infection level and was not subject to parasite-induced mortality. We also found that both acanthocephalan infection levels and parasite-induced mortality varied greatly among cryptic MOTUs of the native amphipods. Contrary to expectations, some native MOTUs displayed levels of resistance to their local parasites similar to those observed in the invasive G. roeseli. Overall, cryptic diversity in native amphipods coupled with high levels of variability in infection levels and parasite-induced mortality documented here may strongly influence inter-MOTU interactions and native population dynamics as well as invasion success and population dynamics of the congeneric invasive G. roeseli.

•

Parasite-mediated competition is an important factor in interspecific interactions.

•

Acanthocephalan infection levels in native and invasive amphipods were assessed.

•

Native amphipods also comprised sympatric, cryptic lineages.

•

Infection levels and host mortality varied greatly among native cryptic host lineages.

•

Some native amphipod lineages were also as resistant to parasites as invasive hosts.

Additive effects of temperature and infection with an acanthocephalan parasite on the shredding activity of Gammarus fossarum (Crustacea: Amphipoda): the importance of aggregative behavior

Climate change can have critical impacts on the ecological role of keystone species, leading to subsequent alterations within ecosystems. The consequences of climate change may be best predicted by understanding its interaction with the cumulative effects of other stressors, although this approach is rarely adopted. However, whether this interaction is additive or interactive can hardly be predicted from studies examining a single factor at a time. In particular, biotic interactions are known to induce modifications in the functional role of many species. Here, we explored the effect of temperature on leaf consumption by a keystone freshwater shredder, the amphipod Gammarus fossarum. This species is found at high densities in the wild and relies on aggregation as an antipredator behavior. In addition, gammarids regularly harbor acanthocephalan parasites that are known to induce multiple effects on their hosts, including modifications on their functional role. We thus assessed the cumulative effect of both intraspecific interactions and parasitism. Consumption tests were conducted on gammarids, either naturally infected with Pomphorhynchus tereticollis or uninfected, feeding alone or in groups. Our results show that increased temperatures induced a significant increase in consumption, but only to a certain extent. Interestingly, consumption at the highest temperature depended on amphipod density: Whereas a decrease was observed for single individuals, no such effect on feeding was observed for individuals in groups. In addition, infection by acanthocephalan parasites per se significantly negatively impacted the shredding role of gammarids. Overall, the combined effects of parasitism and temperature appeared to be additive. Thus, future studies focusing on the impact of climate change on the functional role of keystone species may benefit from a multimodal approach under realistic conditions to derive accurate predictions.

Neogene paleogeography provides context for understanding the origin and spatial distribution of cryptic diversity in a widespread Balkan freshwater amphipod

BACKGROUND

The Balkans are a major worldwide biodiversity and endemism hotspot. Among the freshwater biota, amphipods are known for their high cryptic diversity. However, little is known about the temporal and paleogeographic aspects of their evolutionary history. We used paleogeography as a framework for understanding the onset of diversification in Gammarus roeselii: (1) we hypothesised that, given the high number of isolated waterbodies in the Balkans, the species is characterised by high level of cryptic diversity, even on a local scale; (2) the long geological history of the region might promote pre-Pleistocene divergence between lineages; (3) given that G. roeselii thrives both in lakes and rivers, its evolutionary history could be linked to the Balkan Neogene paleolake system; (4) we inspected whether the Pleistocene decline of hydrological networks could have any impact on the diversification of G. roeselii.

MATERIAL AND METHODS

DNA was extracted from 177 individuals collected from 26 sites all over Balkans. All individuals were amplified for ca. 650 bp long fragment of the mtDNA cytochrome oxidase subunit I (COI). After defining molecular operational taxonomic units (MOTU) based on COI, 50 individuals were amplified for ca. 900 bp long fragment of the nuclear 28S rDNA. Molecular diversity, divergence, differentiation and historical demography based on COI sequences were estimated for each MOTU. The relative frequency, geographic distribution and molecular divergence between COI haplotypes were presented as a median-joining network. COI was used also to reconstruct time-calibrated phylogeny with Bayesian inference. Probabilities of ancestors' occurrence in riverine or lacustrine habitats, as well their possible geographic locations, were estimated with the Bayesian method. A Neighbour Joining tree was constructed to illustrate the phylogenetic relationships between 28S rDNA haplotypes.

RESULTS

We revealed that G. roeselii includes at least 13 cryptic species or molecular operational taxonomic units (MOTUs), mostly of Miocene origin. A substantial Pleistocene diversification within-MOTUs was observed in several cases. We evidenced secondary contacts between very divergent MOTUs and introgression of nDNA. The Miocene ancestors could live in either lacustrine or riverine habitats yet their presumed geographic localisations overlapped with those of the Neogene lakes. Several extant riverine populations had Pleistocene lacustrine ancestors.

DISCUSSION

Neogene divergence of lineages resulting in substantial cryptic diversity may be a common phenomenon in extant freshwater benthic crustaceans occupying areas that were not glaciated during the Pleistocene. Evolution of G. roeselii could be associated with gradual deterioration of the paleolakes. The within-MOTU diversification might be driven by fragmentation of river systems during the Pleistocene. Extant ancient lakes could serve as local microrefugia during that time.

Widely distributed and regionally isolated! Drivers of genetic structure in Gammarus fossarum in a human-impacted landscape

Background

The actual connectivity between populations of freshwater organisms is largely determined by species biology, but is also influenced by many area- and site-specific factors, such as water pollution and habitat fragmentation. Therefore, the prediction of effective gene flow, even for well-studied organisms, is difficult. The amphipod crustacean Gammarus fossarum is a key invertebrate in freshwater ecosystems and contains many cryptic species. One of these species is the broadly distributed G. fossarum clade 11 (type B). In this study, we tested for factors driving the genetic structure of G. fossarum clade 11 in a human-impacted landscape at local and regional scales. To determine population structure, we analyzed the mitochondrial cytochrome c oxidase 1 (CO1) gene of 2,086 specimens from 54 sampling sites and microsatellite loci of 420 of these specimens from ten sites.

Results

We detected strong overall genetic differentiation between populations at regional and local scales with both independent marker systems, often even within few kilometers. Interestingly, we observed only a weak correlation of genetic distances with geographic distances or catchment boundaries. Testing for factors explaining the observed population structure revealed, that it was mostly the colonization history, which has influenced the structure rather than any of the chosen environmental factors. Whereas the number of in-stream barriers did not explain population differentiation, the few large water reservoirs in the catchment likely act as dispersal barriers.

Conclusions

We showed that populations of Gammarus fossarum clade 11 are strongly isolated even at local scales in the human-impacted region. The observed genetic structure was best explained by the effects of random genetic drift acting independently on isolated populations after historical colonization events. Genetic drift in isolated populations was probably further enhanced by anthropogenic impacts, as G. fossarum is sensitive to many anthropogenic stressors. These findings highlight the importance of small-scale genetic studies to determine barriers restricting gene flow to prevent further loss of genetic diversity and maintain intact freshwater ecosystems.

Electronic supplementary material

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

Spatial patterns of genetic diversity, community composition and occurrence of native and non-native amphipods in naturally replicated tributary streams

Background

Worldwide, natural communities are invaded by non-native species, with potentially devastating effects on the native communities. A large part of past research aimed at finding traits and characteristics of the invading species or the invaded community explaining observed invasions. Only recently, the focus shifted on the spatial patterns during invasions per se. Empirical data, however, are limited, as invasions are often unique incidences of a complex spatio-temporal process. In order to identify generalities of invasion patterns, we studied 13 naturally replicated tributary streams draining into Lake Constance, and studied the occurrence of native and non-native amphipods along linear transects from the stream outlets to the upstream headwater reaches.

Results

We found repeated spatial patterns of community composition and the occurrence of native and non-native amphipod species across two different years. Specifically, occurrence as well as abundance of two non-native amphipod species decreased from the stream outlets at the lake site towards upstream headwater reaches. Populations of the most common native amphipod species were largest at the uppermost headwater reaches. All populations of this native species, however, showed significant signals of recent genetic bottlenecks, irrespective of the stream position and occurrence of non-native species. Contrary to our expectations, this native species also showed no longitudinal genetic differentiation within individual tributaries as postulated for headwater versus outlet populations.

Conclusions

Our results indicate that invasions of river-systems may overall follow predictable patterns on the level of spatial distributions and community composition. However, effects of invading organisms on the genetic diversity and genetic structure of native populations observed at larger scales may not necessarily be directly reflected at the scale of smaller tributaries.

Electronic supplementary material

The online version of this article (doi:10.1186/s12898-016-0079-7) contains supplementary material, which is available to authorized users.

Identifying a key host in an acanthocephalan-amphipod system

Trophically transmitted parasites may use multiple intermediate hosts, some of which may be 'key-hosts', i.e. contributing significantly more to the completion of the parasite life cycle, while others may be 'sink hosts' with a poor contribution to parasite transmission. Gammarus fossarum and Gammarus roeseli are sympatric crustaceans used as intermediate hosts by the acanthocephalan Pomphorhynchus laevis. Gammarus roeseli suffers higher field prevalence and is less sensitive to parasite behavioural manipulation and to predation by definitive hosts. However, no data are available on between-host differences in susceptibility to P. laevis infection, making it difficult to untangle the relative contributions of these hosts to parasite transmission. Based on results from estimates of prevalence in gammarids exposed or protected from predation and laboratory infections, G. fossarum specimens were found to be more susceptible to P. laevis infection. As it is more susceptible to both parasite infection and manipulation, G. fossarum is therefore a key host for P. laevis transmission.

Invaders, natives and their enemies: distribution patterns of amphipods and their microsporidian parasites in the Ruhr Metropolis, Germany

Background

The amphipod and microsporidian diversity in freshwaters of a heterogeneous urban region in Germany was assessed. Indigenous and non-indigenous host species provide an ideal framework to test general hypotheses on potentially new host-parasite interactions, parasite spillback and spillover in recently invaded urban freshwater communities.

Methods

Amphipods were sampled in 17 smaller and larger streams belonging to catchments of the four major rivers in the Ruhr Metropolis (Emscher, Lippe, Ruhr, Rhine), including sites invaded and not invaded by non-indigenous amphipods. Species were identified morphologically (hosts only) and via DNA barcoding (hosts and parasites). Prevalence was obtained by newly designed parasite-specific PCR assays.

Results

Three indigenous and five non-indigenous amphipod species were detected. Gammarus pulex was further distinguished into three clades (C, D and E) and G. fossarum more precisely identified as type B. Ten microsporidian lineages were detected, including two new isolates (designated as Microsporidium sp. nov. RR1 and RR2). All microsporidians occurred in at least two different host clades or species. Seven genetically distinct microsporidians were present in non-invaded populations, six of those were also found in invaded assemblages. Only Cucumispora dikerogammari and Dictyocoela berillonum can be unambiguously considered as non-indigenous co-introduced parasites. Both were rare and were not observed in indigenous hosts. Overall, microsporidian prevalence ranged from 50 % (in G. roeselii and G. pulex C) to 73 % (G. fossarum) in indigenous and from 10 % (Dikerogammarus villosus) to 100 % (Echinogammarus trichiatus) in non-indigenous amphipods. The most common microsporidians belonged to the Dictyocoela duebenum- /D. muelleri- complex, found in both indigenous and non-indigenous hosts. Some haplotype clades were inclusive for a certain host lineage.

Conclusions

The Ruhr Metropolis harbours a high diversity of indigenous and non-indigenous amphipod and microsporidian species, and we found indications for an exchange of parasites between indigenous and non-indigenous hosts. No introduced microsporidians were found in indigenous hosts and prevalence of indigenous parasites in non-indigenous hosts was generally low. Therefore, no indication for parasite spillover or spillback was found. We conclude that non-indigenous microsporidians constitute only a minimal threat to the native amphipod fauna. However, this might change e.g. if C. dikerogammari adapts to indigenous amphipod species or if other hosts and parasites invade.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-1036-6) contains supplementary material, which is available to authorized users.

Copper sulphate reduces the metabolic activity of Gammarus fossarum in laboratory and field experiments

The specialised fauna of freshwater springs is affected by contamination of the water with xenobiotics from human activities in the surrounding landscape. We assessed the effects of exposure to toxins in laboratory and field experiments by using copper sulphate as a model substance and Gammarus fossarum Koch, 1836, as the model organism. This amphipod is a common representative of the European spring fauna and copper is a widespread contaminant, mainly from agricultural practice. The experiments were conducted in test chambers placed in flow channels and directly in a spring. The gammarids were fed with conditioned beech leaf discs, which had been exposed to a 0.8 mg Cu/L solution for 96 h. The feeding activity of the amphipods was quantified on the level of the organism; and the respiratory electron transport system (ETS) assay was conducted in order to determine changes on the cellular level in the test organisms. The results show that the feeding activity, when the leaf discs were contaminated with copper, was not significantly different from the control. The ETS activity of the gammarids, which had been feeding on the copper contaminated leaf discs was however significantly reduced. The results followed the same pattern for gammarids from both the laboratory and the spring. By conducting the experiments not only in a laboratory but also directly in a spring in the field, we took a crucial step towards a more realistic approach when examining environmental pollutants on an organism. Our findings demonstrate the importance of conducting experiments out in the field, in natural conditions, as well as in the laboratory.

Diversity and distribution of freshwater amphipod species in Switzerland (Crustacea: Amphipoda)

Amphipods are key organisms in many freshwater systems and contribute substantially to the diversity and functioning of macroinvertebrate communities. Furthermore, they are commonly used as bioindicators and for ecotoxicological tests. For many areas, however, diversity and distribution of amphipods is inadequately known, which limits their use in ecological and ecotoxicological studies and handicaps conservation initiatives. We studied the diversity and distribution of amphipods in Switzerland (Central Europe), covering four major drainage basins, an altitudinal gradient of>2,500 m, and various habitats (rivers, streams, lakes and groundwater). We provide the first provisional checklist and detailed information on the distribution and diversity of all amphipod species from Switzerland. In total, we found 29 amphipod species. This includes 16 native and 13 non-native species, one of the latter (Orchestia cavimana) reported here for the first time for Switzerland. The diversity is compared to neighboring countries. We specifically discuss species of the genus Niphargus, which are often receiving less attention. We also found evidence of an even higher level of hidden diversity, and the potential occurrence of further cryptic species. This diversity reflects the biogeographic past of Switzerland, and suggests that amphipods are ideally suited to address questions on endemism and adaptive radiations, post-glaciation re-colonization and invasion dynamics as well as biodiversity-ecosystem functioning relationships in aquatic systems.

Importance of toxicokinetics for interspecies variation in sensitivity to chemicals

Interspecies variation in sensitivity to synthetic chemicals can be orders of magnitude large. Species traits causing the variation can be related to toxicokinetics (uptake, distribution, biotransformation, elimination) or toxicodynamics (interaction with biological target sites). We present an approach to systematically measure and model the contribution of uptake, biotransformation, internal distribution, and elimination kinetics toward species sensitivity differences. The aim is to express sensitivity as target tissue specific, internal lethal concentrations. A case study with the pesticides diazinon, imidacloprid, and propiconazole and the aquatic invertebrates Gammarus pulex, Gammarus fossarum, and Lymnaea stagnalis illustrates the approach. L. stagnalis accumulates more pesticides than Gammaridae when measured in whole organisms but less in target tissues such as the nervous system. Toxicokinetics, i.e. biotransformation and distribution, explain the higher tolerance of L. stagnalis to the insecticide diazinon when compared to Gammaridae. L. stagnalis was again more tolerant to the other neurotoxicant imidacloprid; however, the difference in sensitivity could not be explained by toxicokinetics alone, indicating the importance of toxicodynamic differences. Sensitivity to propiconazole was comparable among all species and, when expressed as internal lethal concentrations, falls in the range of baseline toxicity.

Genetic heterogeneity reveals on-going speciation and cryptic taxonomic diversity of stream-dwelling gudgeons (Teleostei, Cyprinidae) in the middle Danubian hydrosystem (Hungary)

Although stream-dwelling gudgeons (Cyprinidae, genus: Gobio) are widespread in Central Europe, the taxonomy of this group and the distribution of its species are still unexplored in detail. The aims of our study are to ascertain taxonomic composition and distribution of the former Gobio gobio superspecies in the inner area of the Carpathian Basin. Since the presence of cryptic species is suspected in this area, we examined the taxonomic and phylogenetic relationships of Central European Gobio taxa by sequencing the mitochondrial DNA control region (mtCR). Additionally, we characterized the genetic structure of 27 stream-dwelling gudgeon populations of this area by Amplified Fragment Length Polymorphism (AFLP). Results of mtCR analysis proved the presence of three species already known as G. obtusirostris (dominant in NW-Hungary), G. gobio (sporadic) and G. carpathicus (sporadic). Additionally, the analysis revealed the existence of one doubtful taxon, G. sp1 (dominant in NE-Hungary), and a new isolated haplogroup (dominant in SW-Hungary). Although Network analysis showed significant detachment among haplogroups, their genetic distances were quite small. Therefore Bayesian phylogenetic analysis showed weak nodal support for the branching pattern both for newly described haplotypes, and for the already accepted species. AFLP data showed distinct population structure and a clear pattern of isolation was revealed by distance of stocks. At the same time, level of separation was not affected by the altitudinal position of sites. Moreover we found three major clusters of populations which were separated according to hydrographic regions, and corresponded to the findings of mtCR analysis. Our results suggest the on-going speciation of gudgeons in the Carpathian Basin, however the separation of haplogroups seems to only be an intermediate phase. The discovered natural pattern seems to be only slightly influenced by anthropogenic impacts. Additionally our results put into question the suitability of the recently accepted within Gobio genus taxonomy.

Pyrosequencing: powerful and quantitative sequencing technology

Pyrosequencing is a sequencing-by-synthesis method for DNA analysis that has emerged as a platform not only for de novo sequencing applications, but also for quantitative analysis of genomic methylation, single-nucleotide polymorphisms, and allele quantification. In this unit, we describe a complete workflow from sample to result that is suitable for each of these applications. As cytosine conversion is a key element of successful methylation analysis using pyrosequencing, a support protocol for bisulfite treatment is also included.

Hidden biodiversity in an ecologically important freshwater amphipod: differences in genetic structure between two cryptic species

Cryptic species, i.e. species that are morphologically hard to distinguish, have been detected repeatedly in various taxa and ecosystems. In order to evaluate the importance of this finding, we have to know in how far cryptic species differ in various aspects of their biology. The amphipod Gammarus fossarum is a key invertebrate in freshwater streams and contains several cryptic species. We examined the population genetic structure, genetic diversity and demographic history of two of them (type A and type B) using microsatellite markers and asked whether they show significant differences. We present results of population genetic analyses based on a total of 37 populations from the headwaters of two major European drainages, Rhine and Rhone. We found that, in both species, genetic diversity was geographically structured among and within drainages. For type A in the Rhine and type B in the Rhone, we detected significant patterns of isolation by distance. The increase of genetic differentiation with geographical distance, however, was much higher in type A than in type B. This result indicates substantial interspecific differences in population history and/or the extent of current gene flow between populations. In the Rhine, type B does not show evidence of isolation by distance, and population differentiation is relatively low across hundreds of kilometres. The majority of these populations also show signatures of recent bottlenecks. These patterns are consistent with a recent expansion of type B into the Rhine drainage. In summary, our results suggest considerable and previously unrecognized interspecific differences in the genetic structure of these cryptic keystone species.

Multidimensionality in parasite-induced phenotypic alterations: ultimate versus proximate aspects

In most cases, parasites alter more than one dimension in their host phenotype. Although multidimensionality in parasite-induced phenotypic alterations (PIPAs) seems to be the rule, it has started to be addressed only recently. Here, we critically review some of the problems associated with the definition, quantification and interpretation of multidimensionality in PIPAs. In particular, we confront ultimate and proximate accounts, and evaluate their own limitations. We end up by introducing several suggestions for the development of future research, including some practical guidelines for the quantitative analysis of multidimensionality in PIPAs.

Differences in the sensitivity among cryptic lineages of the Gammarus fossarum complex

Lineages that are at least superficially morphologically identical but genetically distinct are usually misclassified as a single nominal species and, hence, belong to a cryptic species complex, as for example observed for Gammarus fossarum. Since genetic differentiation between cryptic lineages of the G. fossarum complex is comparable to an interspecific level, deviations in physiological and behavioral characteristics are conceivable. The present study investigates for the first time deviations in the sensitivity between two cryptic G. fossarum lineages - namely lineage A and lineage B. Two size classes of both cryptic lineages were exposed in six independent experiments to different concentrations of the fungicide tebuconazole or the insecticide thiacloprid (each n=20) for seven days. Feeding rate on leaf discs was used as a measure of sensitivity. By combining the outcome of all bioassays in a meta-analysis, cryptic lineage A showed - with an approximate 50% more pronounced decline in feeding - a significantly higher overall sensitivity compared to lineage B (n=17). This was also obvious, when each chemical stressor was considered individually; however, data were not significant for thiacloprid. Although other reasons, like land-use of habitats around sampling sites and physiological characteristics of organisms, could not be excluded as driving factors, the results indicate that affiliation to a cryptic lineage is a potential explanatory variable for the observed deviations in sensitivity. As the present study considered only one population for each cryptic lineage, a final conclusion cannot yet be drawn. Hence, it is recommended to investigate behavioral and physiological deviations between cryptic lineages of the G. fossarum complex in further studies.

Vitellogenin-like proteins among invertebrate species diversity: potential of proteomic mass spectrometry for biomarker development

Cost-effective methodologies along with cross-species applicability constitute key points for biomarker development in ecotoxicology. With the advent of cheaper affordable genomic techniques and high throughput sequencing, omics tools could facilitate the assessment of effects of environmental contaminants for all taxa biodiversity. We assessed the potential of absolute quantification of proteins using mass spectrometry to develop vitellogenin (Vg)-like protein assays for invertebrates. We used available sequences in public databases to rapidly identify Vg-proteotypic peptides in seven species from different main taxa of protostome invertebrates (mollusk bivalves, crustacean amphipods, branchiopods, copepods and isopods, and insect diptera). Functional validation was performed by comparing proteomic signals from reproductive female tissue samples and negative controls (male or juvenile tissues). In a second part, we demonstrate in gammarids, daphnids, drosophilids, and gastropods that the assay validated in Vg-sequenced species can be applied to Vg-unsequenced species thanks to the evolutionary conservation of Vg-proteotypic peptide motifs. Finally, we discuss the relevance of mass spectrometry for biomarker development (specific measurement, rapid development, transferability across species). Our study supplies an illustration of the promising strategy to address the challenge of biodiversity in ecotoxicology, which consists in employing omics tools from comparative and evolutionary perspectives.