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

Cryptic species complex shows population-dependent, rather than lineage-dependent tolerance to a neonicotinoid

Cryptic species are rarely considered in ecotoxicology, resulting in misleading outcomes when using a single morphospecies that encompasses multiple cryptic species. This oversight contributes to the lack of reproducibility in ecotoxicological experiments and promotes unreliable extrapolations. The important question of ecological differentiation and the sensitivity of cryptic species is rarely tackled, leaving a substantial knowledge gap regarding the vulnerability of individual cryptic species within species complexes. In times of agricultural intensification and the frequent use of pesticides, there is an urgent need for a better understanding of the vulnerability of species complexes and possible differences in adaptive processes. We used the cryptic species complex of the aquatic amphipod Gammarus roeselii, which comprises at least 13 genetic mtDNA lineages and spans from small-scale endemic lineages in Greece to a large-scale widely distributed lineage in central Europe. We exposed eleven populations belonging to four lineages to the neonicotinoid thiacloprid in an acute toxicity assay. We recorded various environmental variables in each habitat to assess the potential pre-exposure of the populations to contaminants. Our results showed that the populations differed up to 4-fold in their tolerances. The lineage identity had a rather minor influence, suggesting that the cryptic species complex G. roeselii does not differ significantly in tolerance to the neonicotinoid thiacloprid. However, the observed population differentiation implies that recent pre-exposure to thiacloprid (or similar substances) or general habitat contamination has triggered adaptive processes. Though, the extent to which these mechanisms are equally triggered in all lineages needs to be addressed in the future. Our study provides two key findings: Firstly, it shows that observed phylogenetic differences within the G. roeselii species complex did not reveal differences in thiacloprid tolerance. Second, it confirms that differentiation occurs at the population level, highlighting that susceptibility to toxicants is population-dependent. The population-specific differences were within the range of accepted intraspecific variability from a regulatory standpoint. From an evolutionary-ecological perspective, it remains intriguing to observe how persistent stresses will continue to influence tolerance and whether different populations are on distinct pathways of adaptation. Given that the potential selection process has only lasted a relatively short number of generations, it is crucial to monitor these populations in the future, as even brief exposure periods significantly impact evolutionary responses.

Chronic metal contamination shapes the size structure of Gammarus fossarum populations in French headwater rivers

Assessing the effects of multigenerational exposure of aquatic animal populations to chemical contamination is essential for ecological risk assessment. However, beyond rare examples reporting the sporadic emergence of a toxicological tolerance within populations that persist in contaminated environments, conclusive results are even more limited from field studies when it comes to the alteration of life-history traits. Here, we investigated whether long-term exposure to cadmium (Cd) influences size-related life-history traits (i.e., size at puberty, median adult size, maximum size) in Gammarus fossarum, a keystone species of European stream ecosystems. We studied 13 field populations of G. fossarum (cryptic lineage B) living in headwater rivers located in natural areas scattered at a large geographical scale and exposed to contrasted bioavailable Cd contamination levels due to different local geochemical backgrounds. We achieved a detailed description of the physical and physicochemical conditions of the river reaches investigated. Land-use parameters, hydrological characteristics (flow, slope, river width, flow structure, mosaic of substrates), and physicochemical conditions (temperature, conductivity, dissolved oxygen) were measured. Metallic bioavailable contamination was assessed using a standardized active biomonitoring procedure (Gammarus caging). Based on the field demographic census of the 13 populations, our results demonstrated that chronic Cd contamination significantly influences life-history in the G. fossarum species, with a significant reduction in all size traits of populations (size at puberty, median adult size, maximum size). In addition, we confirmed Cd-tolerance in contaminated populations during exposure tests in the laboratory. Various hypotheses can be then put forward to explain the modification of size-related life-history traits: a direct toxic effect of Cd, a cost of Cd-tolerance, or an adaptive evolution of life-history exposed to toxic pressure.

Pollination biology of Impatiens capensis Meerb. in non-native range

Pollination biology in the widespread species Impatiens capensis Meerb. has only been studied in America, specifically in zones of the U.S.A. and Canada. In this study, we investigated the pollination biology of I. capensis using an integrative identification approach using morphological and molecular tools in four populations of Northwest Poland. We also determined and compared the functional characteristics of the pollinators of the introduced species from the study sites and the native ones reported, for the latter collecting information from bibliographic sources. Visitors were identified using standard morphological keys, including identifying and classifying insect mouthparts. Molecular identification was carried out using mitochondrial DNA's cytochrome oxidase subunit I (COI). We morphologically identified 20 species of visitors constituted by 17 pollinators and three nectar robbers. DNA barcoding of 59 individuals proved the identification of 18 species (also 18 BINs). The frequency of pollinator species was primarily made up of representatives of both Hymenoptera (75%) and Diptera (21%). The morphological traits, such as the chewing and sucking mouthparts, small and big body height, and robber and pollinator behavior explained mainly the native and introduced visitors' arrangements that allow pollination success. However, to understand the process comprehensively, further investigation of other causalities in pollination success and understanding the diversity of pollinators in outer native ranges are necessary.

One like all? Behavioral response range of native and invasive amphipods to neonicotinoid exposure

Native and invasive species often occupy similar ecological niches and environments where they face comparable risks from chemical exposure. Sometimes, invasive species are phylogenetically related to native species, e.g. they may come from the same family and have potentially similar sensitivities to environmental stressors due to phylogenetic conservatism and ecological similarity. However, empirical studies that aim to understand the nuanced impacts of chemicals on the full range of closely related species are rare, yet they would help to comprehend patterns of current biodiversity loss and species turnover. Behavioral sublethal endpoints are of increasing ecotoxicological interest. Therefore, we investigated behavioral responses (i.e., change in movement behavior) of the four dominant amphipod species in the Rhine-Main area (central Germany) when exposed to the neonicotinoid thiacloprid. Moreover, beyond species-specific behavioral responses, ecological interactions (e.g. parasitation with Acanthocephala) play a crucial role in shaping behavior, and we have considered these infections in our analysis. Our findings revealed distinct baseline behaviors and species-specific responses to thiacloprid exposure. Notably, Gammarus fossarum exhibited biphasic behavioral changes with hyperactivity at low concentrations that decreased at higher concentrations. Whereas Gammarus pulex, Gammarus roeselii and the invasive species Dikerogammarus villosus, showed no or weaker behavioral responses. This may partly explain why G. fossarum disappears in chemically polluted regions while the other species persist there to a certain degree. But it also shows that potential pre-exposure in the habitat may influence behavioral responses of the other amphipod species, because habituation occurs, and potential hyperactivity would be harmful to individuals in the habitat. The observed responses were further influenced by acanthocephalan parasites, which altered baseline behavior in G. roeselii and enhanced the behavioral response to thiacloprid exposure. Our results underscore the intricate and diverse nature of responses among closely related amphipod species, highlighting their unique vulnerabilities in anthropogenically impacted freshwater ecosystems.

A DNA barcode library of Austrian geometridae (Lepidoptera) reveals high potential for DNA-based species identification

Situated in the Eastern section of the European Alps, Austria encompasses a great diversity of different habitat types, ranging from alpine to lowland Pannonian ecosystems, and a correspondingly high level of species diversity, some of which has been addressed in various DNA barcoding projects. Here, we report a DNA barcode library of all the 476 species of Geometridae (Lepidoptera) that have been recorded in Austria. As far as possible, species were sampled from different Austrian regions in order to capture intraspecific genetic variation. In total, 2500 DNA barcode sequences, representing 438 species, were generated in this study. For complete coverage of Austrian geometrid species in the subsequent analyses, the dataset was supplemented with DNA barcodes from specimens of non-Austrian origin. Species delimitations by ASAP, BIN and bPTP methods yielded 465, 510 and 948 molecular operational taxonomic units, respectively. Congruency of BIN and ASAP partitions with morphospecies assignments was reasonably high (85% of morphospecies in unique partitions), whereas bPTP appeared to overestimate the number of taxonomic units. The study furthermore identified taxonomically relevant cases of morphospecies splitting and sharing in the molecular partitions. We conclude that DNA barcoding and sequence analysis revealed a high potential for accurate DNA-based identification of the Austrian Geometridae species. Additionally, the study provides an updated checklist of the geometrid moths of Austria.

Acclimation and transgenerational plasticity support increased cadmium tolerance in Gammarus populations exposed to natural metal contamination in headwater streams

Considering long-term population effects of chronic exposure to contaminants remains limited in ecological risk assessment. Field evidence that multigenerational exposure influences organisms' sensitivity is still scarce, and mechanisms have yet to be elucidated in the environmental context. This study focuses on the crustacean Gammarus fossarum, for which an increased tolerance to cadmium (Cd) has previously been reported in a naturally low-contaminated headwater stream. Our objectives were to investigate whether Cd tolerance is a common phenomenon in headwater populations, and to elucidate the nature of the tolerance and its intergenerational transmission. For this, we carried out an in-depth in situ characterization of Cd exposure (gammarids' caging) and levels of tolerance in nine populations on a regional scale, as well as laboratory maintenance and cross-breeding of contaminated and uncontaminated populations. Acute tolerance levels correlate positively with bioavailable Cd contamination levels among streams. The contaminated and non-contaminated populations differ about two-fold in sensitivity to Cd. Tolerance was found in all age classes of contaminated populations, it can be transiently lost during the year, and it was transmissible to offspring. In addition, tolerance levels dropped significantly when organisms were transferred to a Cd-free environment for two months. These organisms also ceased producing tolerant offspring, confirming a non-genetic transmission of Cd tolerance between generations. These findings support that Cd tolerance corresponds to non-genetic acclimation combined with transgenerational plasticity. Moreover, cross-breeding revealed that tolerance transmission to offspring is not limited to maternal effect. We suggest epigenetics as a plausible mechanism for the plasticity of Cd sensitivity observed in the field. Our results therefore highlight the neglected role of plasticity and non-genetic transmission of modified sensitivities during the long-term exposure of natural populations to environmental contamination.

Is the relationship between body length and body mass consistent across habitats? A case study on Niphargus (Crustacea: Amphipoda)

Body size is one of the main characteristics of any organism and influences various aspects of individual's biology. In animal ecology, it represents a key functional trait that can be quantified using different measures and is often used as a proxy for different organismal functions. The way we quantify body size is critical in any study using this measure alone or to scale other organismal traits. It is especially important in groups that act as model systems across different fields of biological research. One of such groups are amphipods, which are at focus in many ecological studies where appropriate quantification of body size is needed. Here, we explored the relationship between body length and body mass in the largest freshwater amphipod genus Niphargus, and evaluated whether the two measures lead to different conclusions in a putative ecological study of species coexistence. We selected 16 species inhabiting two different subterranean habitats, cave lakes and cave streams. The relationship between log-transformed body mass and body length was linear in all species, but body mass increased steeper among species from cave lakes than from cave streams, reflecting the stouter body shape of the former. In the simulated ecological study, the comparisons of the two measures showed that they may yield different results: in 10 % of cases, body length detected differences between species when body mass did not and vice versa (13 %). Usage of body length or body mass can thus lead to different conclusions. We recommend avoiding direct transformations between body length and body mass in ecological studies. Whenever needed, such transformations should be done with caution using habitat-specific body mass - body length ratios.

Organ-oriented proteogenomics functional atlas of three aquatic invertebrate sentinel species

Proteogenomic methodologies have enabled the identification of protein sequences in wild species without annotated genomes, shedding light on molecular mechanisms affected by pollution. However, proteomic resources for sentinel species are limited, and organ-level investigations are necessary to expand our understanding of their molecular biology. This study presents proteomic resources obtained from proteogenomic analyses of key organs (hepatopancreas, gills, hemolymph) from three established aquatic sentinel invertebrate species of interest in ecotoxicological/ecological research and environmental monitoring: Gammarus fossarum, Dreissena polymorpha, and Palaemon serratus. Proteogenomic analyses identified thousands of proteins for each species, with over 90% of them being annotated to putative function. Functional analysis validated the relevance of the proteomic atlases by revealing similarities in functional annotation of catalogues of proteins across analogous organs in the three species, while deep contrasts between functional profiles are delimited across different organs in the same organism. These organ-level proteomic atlases are crucial for future research on these sentinel animals, aiding in the evaluation of aquatic environmental risks and providing a valuable resource for ecotoxicological studies.

Cryptic Species in Ecotoxicology

The advent of genetic methods has led to the discovery of an increasing number of species that previously could not be distinguished from each other on the basis of morphological characteristics. Even though there has been an exponential growth of publications on cryptic species, such species are rarely considered in ecotoxicology. Thus, the particular question of ecological differentiation and the sensitivity of closely related cryptic species is rarely addressed. Tackling this question, however, is of key importance for evolutionary ecology, conservation biology, and, in particular, regulatory ecotoxicology. At the same time, the use of species with (known or unknown) cryptic diversity might be a reason for the lack of reproducibility of ecotoxicological experiments and implies a false extrapolation of the findings. Our critical review includes a database and literature search through which we investigated how many of the species most frequently used in ecotoxicological assessments show evidence of cryptic diversity. We found a high proportion of reports indicating overlooked species diversity, especially in invertebrates. In terrestrial and aquatic realms, at least 67% and 54% of commonly used species, respectively, were identified as cryptic species complexes. The issue is less prominent in vertebrates, in which we found evidence for cryptic species complexes in 27% of aquatic and 6.7% of terrestrial vertebrates. We further exemplified why different evolutionary histories may significantly determine cryptic species' ecology and sensitivity to pollutants. This in turn may have a major impact on the results of ecotoxicological tests and, consequently, the outcome of environmental risk assessments. Finally, we provide a brief guideline on how to deal practically with cryptic diversity in ecotoxicological studies in general and its implementation in risk assessment procedures in particular. Environ Toxicol Chem 2023;42:1889-1914. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Avoidance behaviour of aquatic macroinvertebrates for real-time detection of micropollutant surge in wastewater effluents

Micropollutants are regularly detected at the outlets of wastewater treatment plants (WWTPs). Across urban and industrial WWTPs, monitoring directives only require assessment for a handful of chemicals via sampling methods that fail to capture the temporal variability in micropollutant discharge. In this study, we develop a biotest for real-time on-line monitoring of micropollutant discharge dynamics in WWTPs effluents. The selected biomonitoring device ToxMate uses videotracking of invertebrate movement, which was used to deduce avoidance behaviour of the amphipod Gammarus fossarum. Organism conditioning was set up to induce a state of minimal locomotor activity in basal conditions to maximise avoidance signal sensitivity to micropollutant spikes. We showed that with a standardised protocol, it was possible to minimise both overall movement and sensitivity to physio-chemical variations typical to WWTP effluents, as well as capture the spikes of two micropollutants upon exposure (copper and methomyl). Spikes in avoidance behaviour were consistently seen for the two chemicals, as well as a strong correlation between avoidance intensity and spiked concentration. A two-year effluent monitoring case study also illustrates how this biomonitoring method is suitable for real-time on-site monitoring, and shows a promising non-targeted approach for characterising complex micropollutant discharge variability at WWTP effluents, which today remains poorly understood.

Uncovering the Grinnellian niche space of the cryptic species complex Gammarus roeselii

BACKGROUND

The discovery of cryptic species complexes within morphologically established species comes with challenges in the classification and handling of these species. We hardly know to what extent species within a species complex differ ecologically. Such knowledge is essential to assess the vulnerability of individual genetic lineages in the face of global change. The abiotic conditions, i.e., the Grinnellian niche that a genetic lineage colonizes, provides insights into how diverse the ecological requirements of each evolutionary lineage are within a species complex.

MATERIAL AND METHODS

We sampled the cryptic species complex of the amphipod Gammarus roeselii from Central Germany to Greece and identified genetic lineages based on cytochrome c oxidase subunit I (COI) barcoding. At the same time, we recorded various abiotic parameters and local pollution parameters using a series of in vitro assays to then characterize the Grinnellian niches of the morphospecies (i.e., Gammarus roeselii sensu lato) as well as each genetic lineage. Local pollution can be a significant factor explaining current and future distributions in times of increasing production and release of chemicals into surface waters.

RESULTS

We identified five spatially structured genetic lineages in our dataset that differed to varying degrees in their Grinnellian niche. In some cases, the niches were very similar despite the geographical separation of lineages, supporting the hypothesis of niche conservatism while being allopatrically separated. In other cases, we found a small niche that was clearly different from those of other genetic lineages.

CONCLUSION

The variable niches and overlaps of different dimensions make the G. roeselii species complex a promising model system to further study ecological, phenotypic and functional differentiation within this species complex. In general, our results show that the Grinnellian niches of genetically distinct molecular operational taxonomic units (MOTUs) within a cryptic species complex can differ significantly between each other, calling for closer inspection of cryptic species in a conservational and biodiversity context.

Widespread infection, diversification and old host associations of Nosema Microsporidia in European freshwater gammarids (Amphipoda)

The microsporidian genus Nosema is primarily known to infect insects of economic importance stimulating high research interest, while other hosts remain understudied. Nosema granulosis is one of the formally described Nosema species infecting amphipod crustaceans, being known to infect only two host species. Our first aim was to characterize Nosema spp. infections in different amphipod species from various European localities using the small subunit ribosomal DNA (SSU) marker. Second, we aimed to assess the phylogenetic diversity, host specificity and to explore the evolutionary history that may explain the diversity of gammarid-infecting Nosema lineages by performing a phylogenetic reconstruction based on RNA polymerase II subunit B1 (RPB1) gene sequences. For the host species Gammarus balcanicus, we also analyzed whether parasites were in excess in females to test for sex ratio distortion in relation with Nosema infection. We identified Nosema spp. in 316 individuals from nine amphipod species being widespread in Europe. The RPB1-based phylogenetic reconstruction using newly reported sequences and available data from other invertebrates identified 39 haplogroups being associated with amphipods. These haplogroups clustered into five clades (A-E) that did not form a single amphipod-infecting monophyletic group. Closely related sister clades C and D correspond to Nosema granulosis. Clades A, B and E might represent unknown Nosema species infecting amphipods. Host specificity seemed to be variable with some clades being restricted to single hosts, and some that could be found in several host species. We show that Nosema parasite richness in gammarid hosts is much higher than expected, illustrating the advantage of the use of RPB1 marker over SSU. Finally, we found no hint of sex ratio distortion in Nosema clade A infecting G. balcanicus. This study shows that Nosema spp. are abundant, widespread and diverse in European gammarids. Thus, Nosema is as diverse in aquatic as in terrestrial hosts.

Phylogenetically-conserved candidate genes unify biodiversity-ecosystem function relationships and eco-evolutionary dynamics across biological scales

The intra- and interspecific facets of biodiversity have traditionally been analysed separately, limiting our understanding of how evolution has shaped biodiversity, how biodiversity (as a whole) alters ecological dynamics and hence eco-evolutionary feedbacks at the community scale. Here, we propose using candidate genes phylogenetically-conserved across species and sustaining functional traits as an inclusive biodiversity unit transcending the intra- and interspecific boundaries. This framework merges knowledge from functional genomics and functional ecology, and we first provide guidelines and a concrete example for identifying phylogenetically-conserved candidate genes (PCCGs) within communities and for measuring biodiversity from PCCGs. We then explain how biodiversity measured at PCCGs can be linked to ecosystem functions, which unifies recent observations that both intra- and interspecific biodiversity are important for ecosystem functions. We then highlight the eco-evolutionary processes shaping PCCG diversity patterns and argue that their respective role can be inferred from concepts derived from population genetics. Finally, we explain how PCCGs may shift the field of eco-evolutionary dynamics from a focal-species approach to a more realistic focal-community approach. This framework provides a novel perspective to investigate the global ecosystem consequences of diversity loss across biological scales, and how these ecological changes further alter biodiversity evolution.

Morphology, nuclear SNPs and mate selection reveal that COI barcoding overestimates species diversity in a Mediterranean freshwater amphipod by an order of magnitude

DNA sequence information has revealed many morphologically cryptic species worldwide. For animals, DNA-based assessments of species diversity usually rely on the mitochondrial cytochrome c oxidase subunit I (COI) gene. However, a growing amount of evidence indicate that mitochondrial markers alone can lead to misleading species diversity estimates due to mito-nuclear discordance. Therefore, reports of putative species based solely on mitochondrial DNA should be verified by other methods, especially in cases where COI sequences are identical for different morphospecies or where divergence within the same morphospecies is high. Freshwater amphipods are particularly interesting in this context because numerous putative cryptic species have been reported. Here, we investigated the species status of the numerous mitochondrial molecular operational taxonomic units (MOTUs) found within Echinogammarus sicilianus. We used an integrative approach combining DNA barcoding with mate selection observations, detailed morphometrics and genome-wide double digest restriction site-associated DNA sequencing (ddRAD-seq). Within a relatively small sampling area, we detected twelve COI MOTUs (divergence = 1.8-20.3%), co-occurring in syntopy at two-thirds of the investigated sites. We found that pair formation was random and there was extensive nuclear gene flow among the ten MOTUs co-occurring within the same river stretch. The four most common MOTUs were also indistinguishable with respect to functional morphology. Therefore, the evidence best fits the hypothesis of a single, yet genetically diverse, species within the main river system. The only two MOTUs sampled outside the focal area were genetically distinct at the nuclear level and may represent distinct species. Our study reveals that COI-based species delimitation can significantly overestimate species diversity, highlighting the importance of integrative taxonomy for species validation, especially in hyperdiverse complexes with syntopically occurring mitochondrial MOTUs.

A Blind-Identification Test on Criconema annuliferum (de Man, 1921) Micoletzky, 1925 Species Complex Corroborate the Hyper-Cryptic Species Diversity Using Integrative Taxonomy

Ring nematodes are obligate ectoparasites on crops and natural herbaceous and woody plants, and some species are of economic importance and cause damage to roots of several crops. Recent integrative taxonomical analyses recognized the existence of two cryptic species within the Criconema annuliferum morphotype in Spain. In this study, we corroborated that morphometric, morphological and a multi-locus analysis (including the ribosomal markers D2-D3 expansion segments of 28S rRNA, ITS rRNA, 18S RNA, and the mitochondrial DNA cytochrome oxidase I gene) identified a new lineage clearly separated from C. annuliferum, C. paraannuliferum and C. plesioannuliferum. The new lineage was described herein as Criconema pseudoannuliferum sp. nov., confirming that C. annuliferum species complex species complex comprises a hyper-cryptic species complex. This research analysed soil samples from the rhizosphere of maritime pine (Pinus pinaster Ait.) forests in Bermeja-Crestellina Mountain, located at the western part of Málaga province, southern Spain. The integrative taxonomical analyses revealed the occurrence of a new cryptic species identified using females, males and juveniles with detailed morphology, morphometry and molecular markers, described herein as Criconema pseudoannuliferum sp. nov. All molecular markers (D2-D3, ITS, 18S and COI) were obtained from the same individual that was also used for morphological and morphometric analyses. This research demonstrated the hidden diversity within the C. annuliferum species complex species complex can reach to four lineages under ribosomal and mitochondrial gene markers for one morphospecies group, which includes four species, viz. C. annuliferum, C. paraannuliferum, C. plesioannuliferum, and C. pseudoannuliferum sp. nov. Criconema pseudoannuliferum sp. nov. was detected in moderate soil density in two maritime pine forests (5 and 25 nematodes/500 cm³ of soil) suggesting that does not cause damage to maritime pine.

Optimizing laboratory cultures of Gammarus fossarum (Crustacea: Amphipoda) as a study organism in environmental sciences and ecotoxicology

Amphipods are among the most abundant macroinvertebrates in freshwater ecosystems of the Palaearctic and crucial for ecosystem functioning. Furthermore, their high sensitivity to environmental change and pollutants makes them widely used model organisms in environmental sciences and ecotoxicology. In field studies and surveys across Eurasia, species of the genus Gammarus are commonly used, yet laboratory-based studies and ecotoxicological tests are often restricted to the in most parts of the world non-native Hyalella azteca, as Gammarus is much harder to breed and maintain under laboratory conditions. However, for direct comparisons and extrapolations of results of field- vs. laboratory-based studies, the use of the same species would be desirable. Here, we investigated different settings with respect to feeding, shelter and day length to successfully increase survival, juvenile production and their respective growth and survival, and ultimately multi-generation breeding of the amphipod Gammarus fossarum. Amphipod populations persisted and reproduced successfully under optimized husbandry conditions for 12 months and were partially maintained for another year in populations up to a few hundred individuals. Specifically, supplementing diet with protein-rich food sources as well as the provisioning of shelters improved survival rate of G. fossarum significantly. However, we found no significant effect of different day length treatments on the overall relative reproductive activity or on the total amphipod abundance maintained. We conclude that G. fossarum can be kept and reared under standardized conditions. Despite the longer generation times of G. fossarum and higher effort required for maintenance compared to H. azteca, direct ecological relevance and comparability of results to natural systems may justify its future use and development as a study organism for environmental sciences and ecotoxicology.

Multiple reaction monitoring mass spectrometry for the discovery of environmentally modulated proteins in an aquatic invertebrate sentinel species, Gammarus fossarum

Multiple reaction monitoring (MRM) mass spectrometry is emerging as a relevant tool for measuring customized molecular markers in freshwater sentinel species. While this technique is typically used for the validation of protein molecular markers preselected from shotgun experiments, recent gains of MRM multiplexing capacity offer new possibilities to conduct large-scale screening of animal proteomes. By combining the strength of active biomonitoring strategies and MRM technologies, this study aims to propose a new strategy for the discovery of candidate proteins that respond to environmental variability. For this purpose, 249 peptides derived from 147 proteins were monitored by MRM in 273 male gammarids caged in 56 environmental sites, representative of the diversity of French water bodies. A methodology is here proposed to identify a set of customized housekeeping peptides (HKPs) used to correct analytical batch effects and allow proper comparison of peptide levels in gammarids. A comparative analysis performed on HKPs-normalized data resulted in the identification of peptides highly modulated in the environment and derived from proteins likely involved in the environmental stress response. Overall, this study proposes a breakthrough approach to screen and identify potential proteins responding to relevant environmental conditions in sentinel species.

Experimental Crossing Confirms Reproductive Isolation between Cryptic Species within Eulimnogammarus verrucosus (Crustacea: Amphipoda) from Lake Baikal

Ancient lakes are known speciation hotspots. One of the most speciose groups in the ancient Lake Baikal are gammaroid amphipods (Crustacea: Amphipoda: Gammaroidea). There are over 350 morphological species and subspecies of amphipods in Baikal, but the extent of cryptic variation is still unclear. One of the most common species in the littoral zone of the lake, Eulimnogammarus verrucosus (Gerstfeldt, 1858), was recently found to comprise at least three (pseudo)cryptic species based on molecular data. Here, we further explored these species by analyzing their mitogenome-based phylogeny, genome sizes with flow cytometry, and their reproductive compatibility. We found divergent times of millions of years and different genome sizes in the three species (6.1, 6.9 and 8 pg), further confirming their genetic separation. Experimental crossing of the western and southern species, which are morphologically indistinguishable and have adjacent ranges, showed their separation with a post-zygotic reproductive barrier, as hybrid embryos stopped developing roughly at the onset of gastrulation. Thus, the previously applied barcoding approach effectively indicated the separate biological species within E. verrucosus. These results provide new data for investigating genome evolution and highlight the need for precise tracking of the sample origin in any studies in this morphospecies.

A DNA barcode reference library for endemic Ponto-Caspian amphipods

The Ponto-Caspian region is an endemicity hotspot that harbours several crustacean radiations, among which amphipods are the most diverse. These poorly known species are severely threatened in their native range, while at the same time they are invading European inland waters with significant ecological consequences. A proper taxonomic knowledge of this fauna is paramount for its conservation within the native region and monitoring outside of it. Here, we assemble a DNA barcode reference library for nearly 60% of all known Ponto-Caspian amphipod species. We use several methods to define molecular operational taxonomic units (MOTUs), based on two mitochondrial markers (COI and 16S), and assess their congruence with current species-level taxonomy based on morphology. Depending on the method, we find that 54–69% of species had congruent morpho-molecular boundaries. The cases of incongruence resulted from lumping distinct morphospecies into a single MOTU (7–27%), splitting a morphospecies into several MOTUs (4–28%), or both (4–11%). MOTUs defined by distance-based methods without a priori divergence thresholds showed the highest congruence with morphological taxonomy. These results indicate that DNA barcoding is valuable for clarifying the diversity of Ponto-Caspian amphipods, but reveals that extensive work is needed to resolve taxonomic uncertainties. Our study advances the DNA barcode reference library for the European aquatic biota, paving the way towards improved taxonomic knowledge needed to enhance monitoring and conservation efforts.

DNA barcodes evidence the contact zone of eastern and western caddisfly lineages in the Western Carpathians

The region of the Western Carpathians is, among other aspects, very important for survival and diversity of European freshwater fauna due to the presence of a large number of (sub)mountain springs and streams. However, these ecologically and faunistically diversified habitats are still understudied in the context of genetic diversity and population structure of their inhabitants. This study focuses on genetic diversity and distribution patterns of the caddisfly Rhyacophila tristis, common and widespread representative of mountain freshwater fauna. Analysis of the COI mitochondrial marker revealed presence of the western and eastern lineages, with samples from both lineages being grouped in BOLD (Barcode of Life Data System) into separate BINs (Barcode Index Numbers). Our data indicates that eastern lineage (BIN_E) is more closely related to the Balkan populations than to co-occurring western lineage (BIN_W), and that the contact zone of the lineages passes through the W Carpathians. The study revealed phylogeographic and demographic differences between lineages, supporting hypothesis of their evolutionary independence and specific ecological preferences. The obtained genetic data of the R. tristis population from W Carpathians improved our knowledge about population genetics of this aquatic species and can contribute to understanding the state and evolution of biodiversity of freshwater ecosystems in Europe.

First Glimpse at the Diverse Aquaporins of Amphipod Crustaceans

The importance of aquaporins (AQPs) in the transport of water and solutes through cell membranes is well recognized despite being relatively new. To date, despite their abundance, diversity, and presence in disparate environments, amphipods have only been mentioned in studies about the AQPs of other animals and have never been further investigated. In this work, we aimed to recover from public data available AQPs of these crustaceans and reconstruct phylogenetic affinities. We first performed BLAST searches with several queries of diverse taxa against different NCBI databases. Then, we selected the clades of AQPs retrieving the amphipod superfamily Gammaroidea as monophyletic and ran phylogenetic analyses to assess their performances. Our results show how most of the AQPs of amphipods are similar to those of other crustaceans, despite the Prip-like displayed different paralogs, and report for the first time a putative Aqp8-like for arthropods. We also found that the candidate genes of Prip-like, Bib-like, Aqp12-like, and Glp-like help solve deeper relationships in phylogenies of amphipods while leaving uncertainties in shallower parts. With our findings, we hope to increase attention to the study of amphipods as models for AQP functioning and evolution.

DNA metabarcoding reveals differences in distribution patterns and ecological preferences among genetic variants within some key freshwater diatom species

Our study evaluates differences in the distribution and ecology of genetic variants within several ecologically important diatom species that are also key for Water Framework Directive monitoring of European rivers: Fistulifera saprophila (FSAP), Achnanthidium minutissimum (ADMI), Nitzschia inconspicua (NINC) and Nitzschia soratensis (NSTS). We used DADA2 to infer amplicon sequence variants (ASVs) of a short rbcL barcode in 531 environmental samples from biomonitoring campaigns in Catalonia and France. ASVs within each species showed different distribution patterns. Threshold Indicator Taxa ANalysis revealed three ecological groupings of ASVs in both ADMI and FSAP. Two of these in each species were separated by opposite responses to calcium and conductivity. Boosted regression trees additionally showed that both variables greatly influenced the occurrence of these groupings. A third grouping in FSAP was characterized by a negative response to total organic carbon and hence was better represented in waters with higher ecological status than the other FSAP ASVs, contrasting with what is generally assumed for the species. In the two Nitzschia species, our analyses confirmed earlier studies: NINC preferred higher levels of calcium and conductivity. Our findings suggest that the broad ecological tolerance of some diatom species results from overlapping preferences among genetic variants, which individually show much more restricted preferences and distributions. This work shows the importance of studying the ecological preferences of genetic variants within species complexes, now possible with DNA metabarcoding. The results will help reveal and understand biogeographical distributions and facilitate the development of more accurate biological indexes for biomonitoring programmes.

Global gap-analysis of amphipod barcode library

In the age of global climate change and biodiversity loss there is an urgent need to provide effective and robust tools for diversity monitoring. One of the promising techniques for species identification is the use of DNA barcoding, that in Metazoa utilizes the so called 'gold-standard' gene of cytochrome c oxidase (COI). However, the success of this method relies on the existence of trustworthy barcode libraries of the species. The Barcode of Life Data System (BOLD) aims to provide barcodes for all existing organisms, and is complemented by the Barcode Index Number (BIN) system serving as a tool for potential species recognition. Here we provide an analysis of all public COI sequences available in BOLD of the diverse and ubiquitous crustacean order Amphipoda, to identify the barcode library gaps and provide recommendations for future barcoding studies. Our gap analysis of 25,702 records has shown that although 3,835 BINs (indicating putative species) were recognised by BOLD, only 10% of known amphipod species are represented by barcodes. We have identified almost equal contribution of both records (sequences) and BINs associated with freshwater and with marine realms. Three quarters of records have a complete species-level identification provided, while BINs have just 50%. Large disproportions between identification levels of BINs coming from freshwaters and the marine environment were observed, with three quarters of the former possessing a species name, and less than 40% for the latter. Moreover, the majority of BINs are represented by a very low number of sequences rendering them unreliable according to the quality control system. The geographical coverage is poor with vast areas of Africa, South America and the open ocean acting as "white gaps". Several, of the most species rich and highly abundant families of Amphipoda (e.g., Phoxocephalidae, Ampeliscidae, Caprellidae), have very poor representation in the BOLD barcode library. As a result of our study we recommend stronger effort in identification of already recognised BINs, prioritising the studies of families that are known to be important and abundant components of particular communities, and targeted sampling programs for taxa coming from geographical regions with the least knowledge.

Fuzzy species borders of glacial survivalists in the Carpathian biodiversity hotspot revealed using a multimarker approach

The Carpathians are one of the key biodiversity hotspots in Europe. The mountain chain uplifted during Alpine orogenesis and is characterised by a complex geological history. Its current biodiversity was highly influenced by Pleistocene glaciations. The goal of the current study was to examine the phylogenetic and demographic history of Gammarus balcanicus species complex in the Carpathians using multiple markers as well as to delimit, using an integrative approach, and describe new species hidden so far under the name G. balcanicus. Results showed that divergence of the studied lineages reaches back to the Miocene, which supports the hypothesis of their survival in multiple micro refugia. Moreover, the increase of their diversification rate in the Pleistocene suggests that glaciation was the driving force of their speciation. The climatic changes during and after the Pleistocene also played a major role in the demography of the local Carpathian lineages. Comparison of diversity patterns and phylogenetic relationships of both, the mitochondrial and nuclear markers, provide evidence of putative hybridisation and retention of ancient polymorphism (i.e., incomplete lineage sorting). The morphological examination supported the existence of two morphological types; one we describe as a G. stasiuki sp. nov. and another we redescribe as a G. tatrensis (S. Karaman, 1931).

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.

The dark side of pseudoscorpion diversity: The German Barcode of Life campaign reveals high levels of undocumented diversity in European false scorpions

DNA barcoding is particularly useful for identification and species delimitation in taxa with conserved morphology. Pseudoscorpions are arachnids with high prevalence of morphological crypsis. Here, we present the first comprehensive DNA barcode library for Central European Pseudoscorpiones, covering 70% of the German pseudoscorpion fauna (35 out of 50 species). For 21 species, we provide the first publicly available COI barcodes, including the rare Anthrenochernes stellae Lohmander, a species protected by the FFH Habitats Directive. The pattern of intraspecific COI variation and interspecific COI variation (i.e., presence of a barcode gap) generally allows application of the DNA barcoding approach, but revision of current taxonomic designations is indicated in several taxa. Sequences of 36 morphospecies were assigned to 74 BINs (barcode index numbers). This unusually high number of intraspecific BINs can be explained by the presence of overlooked cryptic species and by the accelerated substitution rate in the mitochondrial genome of pseudoscorpions, as known from previous studies. Therefore, BINs may not be an appropriate proxy for species numbers in pseudoscorpions, while partitions built with the ASAP algorithm (Assemble Species by Automatic Partitioning) correspond well with putative species. ASAP delineated 51 taxonomic units from our data, an increase of 42% compared with the present taxonomy. The Neobisium carcionoides complex, currently considered a polymorphic species, represents an outstanding example of cryptic diversity: 154 sequences from our dataset were allocated to 23 BINs and 12 ASAP units.

Mitochondrial Genomes, Phylogenetic Associations, and SNP Recovery for the Key Invasive Ponto-Caspian Amphipods in Europe

The Ponto-Caspian region is the main donor of invasive amphipods to freshwater ecosystems, with at least 13 species successfully established in European inland waters. Dikerogammarus spp. and Pontogammarus robustoides are among the most successful, due to their strong invasive impact on local biota. However, genomic knowledge about these invaders is scarce, while phylogeography and population genetics have been based on short fragments of mitochondrial markers or nuclear microsatellites. In this study, we provide: (i) a reconstruction of six mitogenomes for four invasive gammarids (D. villosus, D. haemobaphes, D. bispinosus, and P. robustoides); (ii) a comparison between the structure of the newly obtained mitogenomes and those from the literature; (iii) SNP calling rates for individual D. villosus and D. haemobaphes from different invasion sites across Europe; and (iv) the first time-calibrated full mitogenome phylogeny reconstruction of several Ponto-Caspian taxa. We found that, in comparison to other gammarids, the mitogenomes of Ponto-Caspian species show a translocation between the tRNA-E and tRNA-R positions. Phylogenetic reconstruction using the mitogenomes identified that Ponto-Caspian gammarids form a well-supported group that originated in the Miocene. Our study supports paraphyly in the family Gammaridae. These provided mitogenomes will serve as vital genetic resources for the development of new markers for PCR-based identification methods and demographic studies.

Aquatic Organisms Research with DNA Barcodes

Since their inception, DNA barcodes have become a powerful tool for understanding the biodiversity and biology of aquatic species, with multiple applications in diverse fields such as food security, fisheries, environmental DNA, conservation, and exotic species detection. Nevertheless, most aquatic ecosystems, from marine to freshwater, are understudied, with many species disappearing due to environmental stress, mostly caused by human activities. Here we highlight the progress that has been made in studying aquatic organisms with DNA barcodes, and encourage its further development in assisting sustainable use of aquatic resources and conservation.

Wide geographic distribution of overlooked parasites: Rare Microsporidia in Gammarus balcanicus, a species complex with a high rate of endemism

Parasites and other symbionts deeply influence host organisms, and no living organism can be considered to have evolved independent of its symbionts. The first step towards understanding symbiotic influences upon host organisms is a strong supporting knowledge of parasite/symbiont diversity. Parasites of freshwater amphipods are diverse, with Microsporidia being a major group. These intracellular parasites impact gammarid fitness in different ways, ranging from reduced fitness to increased fecundity. Many Microsporidia have been recorded using molecular data, with multiple taxa pending formal taxonomic description. While some parasites are common, others are known only through sporadic records of single infections. In this study, we focus on rare/sporadic microsporidian infections within Gammarus balcanicus, a host species complex with a high level of endemism. In addition to enriching our knowledge on Microsporidia parasite diversity in amphipod hosts, we test whether these symbionts are specific to G. balcanicus or if they are the same taxa infecting other gammarid species. Of 2231 hosts from 87 sites, we catalogued 29 sequences of "rare" Microsporidia clustering into 19 haplogroups. These haplogroups cluster into 11 lineages: four pre-described taxa (Cucumispora roeselum, C. ornata, C. dikerogammari and Enterocytospora artemiae) and seven 'Molecular Operational Taxonomic Units', which are known from previously published studies to infect other European amphipod species. Our study significantly widens the geographic range of these Microsporidia and expands the known spectrum of hosts infected. Our results suggest that these parasites are ancient infections of European gammarids. For some host-parasite systems, we hypothesize that the common parasite ancestors that infected the hosts' common ancestors, diversified alongside host diversification. For others, we observe Microsporidia taxa with wide host ranges that do not follow host phylogeny.

Shotgun lipidomics and mass spectrometry imaging unveil diversity and dynamics in Gammarus fossarum lipid composition

Sentinel species are playing an indispensable role in monitoring environmental pollution in aquatic ecosystems. Many pollutants found in water prove to be endocrine disrupting chemicals that could cause disruptions in lipid homeostasis in aquatic species. A comprehensive profiling of the lipidome of these species is thus an essential step toward understanding the mechanism of toxicity induced by pollutants. Both the composition and spatial distribution of lipids in freshwater crustacean Gammarus fossarum were extensively examined herein. The baseline lipidome of gammarids of different sex and reproductive stages was established by high throughput shotgun lipidomics. Spatial lipid mapping by high resolution mass spectrometry imaging led to the discovery of sulfate-based lipids in hepatopancreas and their accumulation in mature oocytes. A diverse and dynamic lipid composition in G. fossarum was uncovered, which deepens our understanding of the biochemical changes during development and which could serve as a reference for future ecotoxicological studies.

•

Baseline lipidome profiling of G. fossarum of different sex and reproductive stages

•

Spatial localization of lipids in gammarid tissue by mass spectrometry imaging

•

SIMS imaging guided discovery of sulfate-based lipids in hepatopancreas epithelium

•

Disclosure of a dynamic lipid composition in maturing female oocytes