Comparative mitogenomic analysis of subterranean and surface amphipods (Crustacea, Amphipoda) with special reference to the family Crangonyctidae

Mitochondrial genomes play important roles in studying genome evolution, phylogenetic analyses, and species identification. Amphipods (Class Malacostraca, Order Amphipoda) are one of the most ecologically diverse crustacean groups occurring in a diverse array of aquatic and terrestrial environments globally, from freshwater streams and lakes to groundwater aquifers and the deep sea, but we have a limited understanding of how habitat influences the molecular evolution of mitochondrial energy metabolism. Subterranean amphipods likely experience different evolutionary pressures on energy management compared to surface-dwelling taxa that generally encounter higher levels of predation and energy resources and live in more variable environments. In this study, we compared the mitogenomes, including the 13 protein-coding genes involved in the oxidative phosphorylation (OXPHOS) pathway, of surface and subterranean amphipods to uncover potentially different molecular signals of energy metabolism between surface and subterranean environments in this diverse crustacean group. We compared base composition, codon usage, gene order rearrangement, conducted comparative mitogenomic and phylogenomic analyses, and examined evolutionary signals of 35 amphipod mitogenomes representing 13 families, with an emphasis on Crangonyctidae. Mitogenome size, AT content, GC-skew, gene order, uncommon start codons, location of putative control region (CR), length of rrnL and intergenic spacers differed between surface and subterranean amphipods. Among crangonyctid amphipods, the spring-dwelling Crangonyx forbesi exhibited a unique gene order, a long nad5 locus, longer rrnL and rrnS loci, and unconventional start codons. Evidence of directional selection was detected in several protein-encoding genes of the OXPHOS pathway in the mitogenomes of surface amphipods, while a signal of purifying selection was more prominent in subterranean species, which is consistent with the hypothesis that the mitogenome of surface-adapted species has evolved in response to a more energy demanding environment compared to subterranean amphipods. Overall, gene order, locations of non-coding regions, and base-substitution rates points to habitat as an important factor influencing the evolution of amphipod mitogenomes.

Functional trait dataset of European groundwater Amphipoda: Niphargidae and Typhlogammaridae

Groundwater represents a vast, but mostly hidden and inaccessible ecosystem. Although often overlooked in freshwater research, groundwater organisms form a significant part of freshwater biodiversity, whereas their functions are crucial in different ecosystem processes. Knowledge on functional traits is generally lacking for most groundwater species worldwide, yet European groundwater amphipods, particularly the family Niphargidae, are an exception. They are well-researched and used as a model system in ecological and evolutionary studies. We focused on this group to assemble a first functional trait dataset dedicated to groundwater species. We gathered data for eight morphological functional traits quantified through 27 measurements for 1123 individuals which represent 180 species and 314 MOTUs. Besides functional trait data, every entry is accompanied with locality information, including habitat type, and DNA sequences if available. The structure of the dataset and data processing information provided along enable wide applicability and extension to other amphipod taxa. When coupled with phylogeny, the dataset may further enhance different aspects of groundwater research, including biodiversity patterns, community assembly processes, and trait evolution.

Continent-island boundary and environment-shaped evolution in the marine amphipod Ampithoe marcuzzii complex (Crustacea: Eumalacostraca: Ampithoidae)

Marine amphipods are crustaceans that lack a larval phase and consequently have low dispersion rates. Despite that, these crustaceans present a remarkable ability to be transported by rafting on natural floating substrata, especially macroalgae, where they find shelter, food and a mating ground. The species Ampithoe marcuzzii is widely distributed throughout the western Atlantic Ocean. Here, it was used as a model to study seascape genomics and phylogeography in invertebrates with low dispersion capacities. We anticipated that the lineages would present isolation-by-distance patterns. However, surface currents and other abiotic variables could facilitate connectivity among distant sites. Based on mitochondrial and nuclear genes, SNPs, and environmental associations, we observed the presence of a species complex within A. marcuzzii, separating mainland and insular populations. Each species showed an independent evolutionary history, with a strong latitudinal population structure and evidence of isolation-by-distance and isolation-by-environment, characterizing the 'continent' species. Historical expansion and environmental variables were observed associated with the southeastern population, and ecological niche modeling corroborated the region as a paleorefuge. Conversely, populations from 'islands' presented complicated evolutionary histories, with closer localities genetically isolated and distant localities connected. These findings indicate that insular populations with low dispersion capacity might be more susceptible to spatial connectivity by floating substrata and to changes in surface currents. In contrast, mainland populations might be more vulnerable to local climate changes due to lack of gene flow.

Identification and organ-specific patterns of expression of two metallothioneins in the sentinel species Gammarus fossarum

Metal pollution is a major concern for aquatic environments. Widespread contamination by various trace metal ions has been described in freshwater streams as well as their subsequent bioaccumulation, potentially leading to toxicity and trophic transfer. Metallothioneins constitute an evolutionary conserved family of low molecular weight, cysteine-rich, metal-chelating proteins, whose known physiological functions are the maintenance of the homeostasis of essential metals, the detoxification of non-essential metals, and the protection against oxidative stress and free radicals. In this study, we identified two metallothionein-coding transcripts, mt1 and mt2, in the transcriptome of the amphipod Gammarus fossarum, a sentinel species widely used to assess the quality of watersheds. For the first time, we investigated the organ-specific patterns of expression of these two mt transcripts at the individual level in the gills and the caeca of this small crustacean. In silico analysis and experimental exposures to environmentally relevant concentrations of cadmium, zinc and silver showed that G. fossarum mt1 induction is stronger after Cd exposure compared to the other tested metals. G. fossarum mt1 was more significantly induced in the caeca than in the gills of exposed organisms for any metal exposure, while G. fossarum mt2 was, at least at the individual level, more inducible in the gills than in the caeca of G. fossarum exposed to Cd and Zn. Our results provide new genetic resources that will help to improve the understanding of metal homeostasis in this sentinel species.

Integrating citizen science and environmental DNA metabarcoding to study biodiversity of groundwater amphipods in Switzerland

Groundwater is the physically largest freshwater ecosystem, yet one of the least explored habitats on earth, both because of accessing difficulties and the scarcity of the organisms inhabiting it. Here, we demonstrate how a two-fold approach provides complementary information on the occurrence and diversity of groundwater amphipods. Firstly, we used a citizen science approach in collaboration with municipal water providers who sampled groundwater organisms in their spring catchment boxes over multiple weeks, followed by DNA barcoding. Secondly, we collected four 10 L water samples at each site, in one sampling event, for environmental DNA (eDNA) metabarcoding. We found that citizen science was very effective in describing the distribution and abundance of groundwater amphipods. Although the single time-point of eDNA sampling did not detect as many amphipods, it allowed the assessment of the entire groundwater community, including microorganisms. By combining both methods, we found different amphipod species co-occurring with distinct sequences from the eDNA-metabarcoding dataset, representing mainly micro-eukaryotic species. We also found a distinct correlation between the diversity of amphipods and the overall biodiversity of groundwater organisms detected by eDNA at each site. We thus suggest that these approaches can be used to get a better understanding of subterranean biodiversity.

Parasite manipulation of host phenotypes inferred from transcriptional analyses in a trematode-amphipod system

Manipulation of host phenotypes by parasites is hypothesized to be an adaptive strategy enhancing parasite transmission across hosts and generations. Characterizing the molecular mechanisms of manipulation is important to advance our understanding of host-parasite coevolution. The trematode (Levinseniella byrdi) is known to alter the colour and behaviour of its amphipod host (Orchestia grillus) presumably increasing predation of amphipods which enhances trematode transmission through its life cycle. We sampled 24 infected and 24 uninfected amphipods from a salt marsh in Massachusetts to perform differential gene expression analysis. In addition, we constructed novel genomic tools for O. grillus including a de novo genome and transcriptome. We discovered that trematode infection results in upregulation of amphipod transcripts associated with pigmentation and detection of external stimuli, and downregulation of multiple amphipod transcripts implicated in invertebrate immune responses, such as vacuolar ATPase genes. We hypothesize that suppression of immune genes and the altered expression of genes associated with coloration and behaviour may allow the trematode to persist in the amphipod and engage in further biochemical manipulation that promotes transmission. The genomic tools and transcriptomic analyses reported provide new opportunities to discover how parasites alter diverse pathways underlying host phenotypic changes in natural populations.

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.

Toxicity and mutagenicity studies of 6PPD-quinone in a marine invertebrate species and bacteria

N-(1,3-dimethylbutyl)-N'-phenyl-p-phenylenediamine-quinone (6PPD-quinone), an oxidation product of the tire additive, 6PPD, has been associated with high mortality of salmonids (0.1 μg/L). The objective of this study was to determine the acute toxicity using neonates and mutagenicity (micronuclei in hemolymph of exposed adults) of 6PPD-quinone in the marine amphipod Parhyale hawaiensis. Also, we studied its mutagenicity in the Salmonella/microsome assay using five strains of Salmonella with and without metabolic system (rat liver S9, 5%). 6PPD-quinone did not present acute toxicity to P. hawaiensis from 31.25 to 500 μg/L. Micronuclei frequency increased after 96 h-exposure to 6PPD-quinone (250 and 500 μg/L) when compared to the negative control. 6PPD-quinone also showed a weak mutagenic effect for TA100 only in the presence of S9. We conclude that 6PPD-quinone is mutagenic to P. hawaiensis and weakly mutagenic to bacteria. Our work provides information for future risk assessment of the presence of 6PPD-quinone in the aquatic environment.

Using DNA Metabarcoding to Characterize the Prey Spectrum of Two Co-Occurring Themisto Amphipods in the Rapidly Changing Atlantic-Arctic Gateway Fram Strait

The two congeneric hyperiids Themisto libellula and T. abyssorum provide an important trophic link between lower and higher trophic levels in the rapidly changing Arctic marine ecosystem. These amphipods are characterized by distinct hydrographic affinities and are hence anticipated to be impacted differently by environmental changes, with major consequences for the Arctic food web. In this study, we applied DNA metabarcoding to the stomach contents of these Themisto species, to comprehensively reveal their prey spectra at an unprecedented-high-taxonomic-resolution and assess the regional variation in their diet across the Fram Strait. Both species feed on a wide variety of prey but their diet strongly differed in the investigated summer season, showing overlap for only a few prey taxa, such as calanoid copepods. The spatially structured prey field of T. libellula clearly differentiated it from T. abyssorum, of which the diet was mainly dominated by chaetognaths. Our approach also allowed the detection of previously overlooked prey in the diet of T. libellula, such as fish species and gelatinous zooplankton. We discuss the reasons for the differences in prey spectra and which consequences these may have in the light of ongoing environmental changes.

Predicting Resistance: Quantifying the Relationship between Urban Development, Agricultural Pesticide Use, and Pesticide Resistance in a Nontarget Amphipod

Resistance alleles within the voltage-gated sodium channel (vgsc) have been correlated with pyrethroid resistance in wild populations of the nontarget amphipod, Hyalella azteca from California (CA), U.S.A. In the present study, we expand upon the relationship between land use and the evolution of pesticide resistance in H. azteca to develop a quantitative methodology to target and screen novel populations for resistance allele genotypes in a previously uninvestigated region of the U.S. (New England: NE). By incorporating urban land development and toxicity-normalized agricultural pesticide use indices into our site selection, we successfully identified three amino acid substitutions associated with pyrethroid resistance. One of the resistance mutations has been described in H. azteca from CA (L925I). We present the remaining two (vgsc I936F and I936V) as novel pyrethroid-resistance alleles in H. azteca based on previous work in insects and elevated cyfluthrin resistance in one NE population. Our results suggest that urban pesticide use is a strong driver in the evolution of resistance alleles in H. azteca. Furthermore, our method for resistance allele screening provides an applied framework for detecting ecosystem impairment on a nationwide scale that can be incorporated into ecological risk assessment decisions.

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.

Chromosomal painting in Charadrius collaris Vieillot, 1818 and Vanellus chilensis Molina, 1782 and an analysis of chromosomal signatures in Charadriiformes

Charadriiformes represent one of the largest orders of birds; members of this order are diverse in morphology, behavior and reproduction, making them an excellent model for studying evolution. It is accepted that the avian putative ancestral karyotype, with 2n = 80, remains conserved for about 100 million years. So far, only a few species of Charadriiformes have been studied using molecular cytogenetics. Here, we performed chromosome painting on metphase chromosomes of two species of Charadriidae, Charadrius collaris and Vanellus chilensis, with whole chromosome paint probes from Burhinus oedicnemus. Charadrius collaris has a diploid number of 76, with both sex chromosomes being submetacentric. In V. chilensi a diploid number of 78 was identified, and the Z chromosome is submetacentric. Chromosome painting suggests that chromosome conservation is a characteristic common to the family Charadriidae. The results allowed a comparative analysis between the three suborders of Charadriiformes and the order Gruiformes using chromosome rearrangements to understand phylogenetic relationships between species and karyotypic evolution. However, the comparative analysis between the Charadriiformes suborders so far has not revealed any shared rearrangements, indicating that each suborder follows an independent evolutionary path, as previously proposed. Likewise, although the orders Charadriiformes and Gruiformes are placed on sister branches, they do not share any signature chromosomal rearrangements.

Gut Microbial Divergence Between Three Hadal Amphipod Species from the Isolated Hadal Trenches

Amphipods are the dominant scavenging metazoan species in the hadal trenches at water depths below 6,000 m. The gut microbiota have been considered to be contribution to the adaptation of deep-sea organisms; however, few comparative analyses of animal gut microbiota between different isolated hadal environments have been done so far. Here, we employed high-throughput 16S rRNA sequencing to compare the gut microbial taxonomic composition and functional potential diversity of three hadal amphipod species, Hirondellea gigas, Bathycallisoma schellenbergi, and Alicella gigantea, collected from the Mariana Trench, Marceau Trench, and New Britain Trench in the Pacific Ocean, respectively. Results showed that different community compositions were detected across all the amphipod specimens based on the analyses of alpha-diversity, hierarchical cluster tree, and PCoA (principal coordinate analysis). Moreover, almost no correlation was observed between genera overrepresented in different amphipods by microbe-microbe correlations analysis, which suggested that the colonization of symbionts were host-specific. At genus level, Psychromonas was dominant in H. gigas, and Candidatus Hepatoplasma was overall dominant in A. gigantea and B. schellenbergi. Comparison of the functional potential showed that, though three hadal amphipod species shared the same predominant functional pathways, the abundances of those most shared pathways showed distinct differences across all the specimens. These findings pointed to the enrichment of particular functional pathways in the gut microbiota of the different isolated trench amphipods. Moreover, in terms of species relative abundance, alpha-diversity and beta-diversity, there was high similarity of gut microbiota between the two A. gigantea populations, which dwelled in two different localities of the same hadal trench. Altogether, this study provides an initial investigation into the gut-microbial interactions and evolution at the hadal depths within amphipod. Each of these three amphipod species would be a model taxa for future studies investigating the influence habitat difference and geography on gut-microbial communities.

Regulation of gene expression during ontogeny of physiological function in the brackishwater amphipod Gammarus chevreuxi

Embryonic development is a complex process involving the co-ordinated onset and integration of multiple morphological features and physiological functions. While the molecular basis of morphological development in embryos is relatively well known for traditional model species, the molecular underpinning of the development of physiological functions is not. Here, we used global gene expression profiling to investigate the transcriptional changes associated with the development of morphological and physiological function in the amphipod crustacean Gammarus chevreuxi. We compared the transcriptomes at three timepoints during the latter half of development, characterised by different stages of the development of heart form and function: 10 days post fertilisation (dpf, Early: no heart structure visible), 15 dpf (Middle: heart present but not fully functional), and 18 dpf (Late: regular heartbeat). Gene expression profiles differed markedly between developmental stages, likely representing a change in the activity of different processes throughout the latter period of G. chevreuxi embryonic development. Differentially expressed genes belonged to one of three distinct clusters based on their expression patterns across development. One of these clusters, which included key genes relating to cardiac contractile machinery and calcium handling, displayed a pattern of sequential up-regulation throughout the developmental period studied. Further analyses of these transcripts could reveal genes that may influence the onset of a regular heartbeat. We also identified morphological and physiological processes that may occur alongside heart development, such as development of digestive caeca and the cuticle. Elucidating the mechanisms underpinning morphological and physiological development of non-model organisms will support improved understanding of conserved mechanisms, addressing the current phylogenetic gap between relatively well known model species.

Micro-endemic species of snails and amphipods show population genetic structure across very small geographic ranges

Understanding variation in population genetic structure, even across small distances and for species with extremely limited ranges, is critical for conservation planning and the development of effective management strategies for imperiled species. Organisms that occupy the same geographic extent can maintain different population structures, ranging from highly diverged to panmictic. Such differences can result from differences in biological characteristics such as dispersal ability or demographic history. We used microsatellite loci to evaluate population genetic structure and variation of four desert spring invertebrates having high to low dispersal ability: the lung snail Physa acuta, two species of gilled snails (Juturnia kosteri and Pyrgulopsis roswellensis; family Hydrobiidae) and the amphipod Gammarus desperatus. The study location represents entire species ranges for the micro-endemic hydrobiids and G. desperatus, while P. acuta is ubiquitous throughout much of North America. We found little evidence of significant population genetic structure for P. acuta and J. kosteri, but much more for P. roswellensis and G. desperatus. Our results demonstrate differences in habitat preference and/or dispersal ability between the species. This information provides insight into how gene flow shapes varying population genetic structure between species across small spatial scales (<100 km2). Most importantly, our results suggest that conservation agencies should not consider these micro-endemic species to be composed of single populations, but rather, that management plans for such species should account for population genetic variation across the species' ranges.

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

The contribution of detoxification pathways to pyrethroid resistance in Hyalella azteca

Chronic exposure to pyrethroid insecticides can result in strong selective pressures on non-target species in aquatic systems and drive the evolution of resistance and population-level changes. Characterizing the underlying mechanisms of resistance is essential to better understanding the potential consequences of contaminant-driven microevolution. The current study found that multiple mechanisms enhance the overall tolerance of Hyalella azteca to the pyrethroid permethrin. In H. azteca containing mutations in the voltage-gated sodium channel (VGSC), both adaptation and acclimation played a role in mitigating the adverse effects of pyrethroid exposures. Pyrethroid resistance is primarily attributed to the heritable mutation at a single locus of the VGSC, resulting in reduced target-site sensitivity. However, additional pyrethroid tolerance was conferred through enhanced enzyme-mediated detoxification. Cytochrome P450 monooxygenases (CYP450) and general esterases (GE) significantly contributed to the detoxification of permethrin in H. azteca. Over time, VGSC mutated H. azteca retained most of their pyrethroid resistance, though there was some increased sensitivity from parent to offspring when reared in the absence of pyrethroid exposure. Permethrin median lethal concentrations (LC50s) declined from 1809 ng/L in parent (P₀) individuals to 1123 ng/L in the first filial (F₁) generation, and this reduction in tolerance was likely related to alterations in acclimation mechanisms, rather than changes to target-site sensitivity. Enzyme bioassays indicated decreased CYP450 and GE activity from P₀ to F₁, whereas the VGSC mutation was retained. The permethrin LC50s in resistant H. azteca were still two orders-of-magnitude higher than non-resistant populations indicating that the largest proportion of resistance was maintained through the inherited VGSC mutation. Thus, the noted variation in tolerance in H. azteca is likely associated with inducible traits controlling enzyme pathways. A better understanding of the mechanistic and genomic basis of acclimation is necessary to more accurately predict the ecological and evolutionary consequences of contaminant-driven change in H. azteca.

The diversity of opsins in Lake Baikal amphipods (Amphipoda: Gammaridae)

BACKGROUND

Vision is a crucial sense for the evolutionary success of many animal groups. Here we explore the diversity of visual pigments (opsins) in the transcriptomes of amphipods (Crustacea: Amphipoda) and conclude that it is restricted to middle (MWS) and long wavelength-sensitive (LWS) opsins in the overwhelming majority of examined species.

RESULTS

We evidenced (i) parallel loss of MWS opsin expression in multiple species (including two independently evolved lineages from the deep and ancient Lake Baikal) and (ii) LWS opsin amplification (up to five transcripts) in both Baikal lineages. The number of LWS opsins negatively correlated with habitat depth in Baikal amphipods. Some LWS opsins in Baikal amphipods contained MWS-like substitutions, suggesting that they might have undergone spectral tuning.

CONCLUSIONS

This repeating two-step evolutionary scenario suggests common triggers, possibly the lack of light during the periods when Baikal was permanently covered with thick ice and its subsequent melting. Overall, this observation demonstrates the possibility of revealing climate history by following the evolutionary changes in protein families.

Chemical Pollution Levels in a River Explain Site-Specific Sensitivities to Micropollutants within a Genetically Homogeneous Population of Freshwater Amphipods

Anthropogenic micropollutants alter chemical and ecological conditions of freshwater ecosystems and impact aquatic species that live along the pollution gradient of a river. Species sensitivity to micropollutants depends on the site-specific exposure; however, it remains unclear to what degree this sensitivity relates to the species' genetic structure. Here, we explored the relationship between the toxic sensitivity and genetic structure of the amphipod species Gammarus pulex (Linnaeus, 1758) along an organic micropollutant gradient in the Holtemme River in central Germany. We determined the river's site-specific micropollutant patterns and analyzed the genetic structure of G. pulex using nuclear and mitochondrial genetic markers. Furthermore, we examined the exposure sensitivities and bioaccumulation of the commonly detected insecticide imidacloprid in G. pulex from different sites. Our results show that throughout the Holtemme River, G. pulex forms a well-connected and homogeneous population with no observable pollution-related differences in the genetic structure. However, G. pulex from polluted sites responded more sensitively to imidacloprid; survival times for half of the amphipods were up to 54% shorter, the percentage of immobile individuals increased up to 65%, and the modeled imidacloprid depuration rate was lower in comparison to amphipods from non-polluted sites. Altogether, these results suggest that the level of sensitivity of G. pulex amphipods to micropollutants in the river depends on the degree of pollution: amphipods may thrive in food-rich but polluted habitats; yet, their sensitivity is increased when chronically exposed to organic micropollutants.

Biosynthesis of Long-Chain Polyunsaturated Fatty Acids in Marine Gammarids: Molecular Cloning and Functional Characterisation of Three Fatty Acyl Elongases

Long-chain (C20-24) polyunsaturated fatty acids (LC-PUFAs) are essential nutrients that are mostly produced in marine ecosystems. Previous studies suggested that gammarids have some capacity to endogenously produce LC-PUFAs. This study aimed to investigate the repertoire and functions of elongation of very long-chain fatty acid (Elovl) proteins in gammarids. Our results show that gammarids have, at least, three distinct elovl genes with putative roles in LC-PUFA biosynthesis. Phylogenetics allowed us to classify two elongases as Elovl4 and Elovl6, as they were bona fide orthologues of vertebrate Elovl4 and Elovl6. Moreover, a third elongase was named as "Elovl1/7-like" since it grouped closely to the Elovl1 and Elovl7 found in vertebrates. Molecular analysis of the deduced protein sequences indicated that the gammarid Elovl4 and Elovl1/7-like were indeed polyunsaturated fatty acid (PUFA) elongases, whereas Elovl6 had molecular features typically found in non-PUFA elongases. This was partly confirmed in the functional assays performed on the marine gammarid Echinogammarus marinus Elovl, which showed that both Elovl4 and Elovl1/7-like elongated PUFA substrates ranging from C18 to C22. E. marinus Elovl6 was only able to elongate C18 PUFA substrates, suggesting that this enzyme does not play major roles in the LC-PUFA biosynthesis of gammarids.

Comparative Mitogenomics in Hyalella (Amphipoda: Crustacea)

We present the sequencing and comparative analysis of 17 mitochondrial genomes of Nearctic and Neotropical amphipods of the genus Hyalella, most from the Andean Altiplano. The mitogenomes obtained comprised the usual 37 gene-set of the metazoan mitochondrial genome showing a gene rearrangement (a reverse transposition and a reversal) between the North and South American Hyalella mitogenomes. Hyalella mitochondrial genomes show the typical AT-richness and strong nucleotide bias among codon sites and strands of pancrustaceans. Protein-coding sequences are biased towards AT-rich codons, with a preference for leucine and serine amino acids. Numerous base changes (539) were found in tRNA stems, with 103 classified as fully compensatory, 253 hemi-compensatory and the remaining base mismatches and indels. Most compensatory Watson–Crick switches were AU -> GC linked in the same haplotype, whereas most hemi-compensatory changes resulted in wobble GU and a few AC pairs. These results suggest a pairing fitness increase in tRNAs after crossing low fitness valleys. Branch-site level models detected positive selection for several amino acid positions in up to eight mitochondrial genes, with atp6 and nad5 as the genes displaying more sites under selection.

Phylogenomics of the Hyalella amphipod species-flock of the Andean Altiplano

Species diversification in ancient lakes has enabled essential insights into evolutionary theory as they embody an evolutionary microcosm compared to continental terrestrial habitats. We have studied the high-altitude amphipods of the Andes Altiplano using mitogenomic, nuclear ribosomal and single-copy nuclear gene sequences obtained from 36 Hyalella genomic libraries, focusing on species of the Lake Titicaca and other water bodies of the Altiplano northern plateau. Results show that early Miocene South American lineages have recently (late Pliocene or early Pleistocene) diversified in the Andes with a striking morphological convergence among lineages. This pattern is consistent with the ecological opportunities (access to unoccupied resources, initial relaxed selection on ecologically-significant traits and low competition) offered by the lacustrine habitats established after the Andean uplift.

Molecular data suggest multiple origins and diversification times of freshwater gammarids on the Aegean archipelago

Our main aim was to investigate the diversity, origin and biogeographical affiliations of freshwater gammarids inhabiting the Aegean Islands by analysing their mtDNA and nDNA polymorphism, thereby providing the first insight into the phylogeography of the Aegean freshwater gammarid fauna. The study material was collected from Samothraki, Lesbos, Skyros, Evia, Andros, Tinos and Serifos islands as well as from mainland Greece. The DNA extracted was used for amplification of two mitochondrial (COI and 16S) and two nuclear markers (28S and EF1-alpha). The multimarker time-calibrated phylogeny supports multiple origins and different diversification times for the studied taxa. Three of the sampled insular populations most probably represent new, distinct species as supported by all the delimitation methods used in our study. Our results show that the evolution of freshwater taxa is associated with the geological history of the Aegean Basin. The biogeographic affiliations of the studied insular taxa indicate its continental origin, as well as the importance of the land fragmentation and the historical land connections of the islands. Based on the findings, we highlight the importance of studying insular freshwater biota to better understand diversification mechanisms in fresh waters as well as the origin of studied Aegean freshwater taxa.

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.

Widespread Torix Rickettsia in New Zealand amphipods and the use of blocking primers to rescue host COI sequences

Endosymbionts and intracellular parasites are common in arthropod hosts. As a consequence, (co)amplification of untargeted bacterial sequences has been occasionally reported as a common problem in DNA barcoding. While identifying amphipod species with universal COI primers, we unexpectedly detected rickettsial endosymbionts belonging to the Torix group. To map the distribution and diversity of Rickettsia species among amphipod hosts, we conducted a nationwide molecular screening of seven families of New Zealand freshwater amphipods. In addition to uncovering a diversity of Torix Rickettsia species across multiple amphipod populations from three different families, our research indicates that: (1) detecting Torix Rickettsia with universal primers is not uncommon, (2) obtaining 'Rickettsia COI sequences' from many host individuals is highly likely when a population is infected, and (3) obtaining 'host COI' may not be possible with a conventional PCR if an individual is infected. Because Rickettsia COI is highly conserved across diverse host taxa, we were able to design blocking primers that can be used in a wide range of host species infected with Torix Rickettsia. We propose the use of blocking primers to circumvent problems caused by unwanted amplification of Rickettsia and to obtain targeted host COI sequences for DNA barcoding, population genetics, and phylogeographic studies.

Excessive Parallelism in Protein Evolution of Lake Baikal Amphipod Species Flock

Repeated emergence of similar adaptations is often explained by parallel evolution of underlying genes. However, evidence of parallel evolution at amino acid level is limited. When the analyzed species are highly divergent, this can be due to epistatic interactions underlying the dynamic nature of the amino acid preferences: The same amino acid substitution may have different phenotypic effects on different genetic backgrounds. Distantly related species also often inhabit radically different environments, which makes the emergence of parallel adaptations less likely. Here, we hypothesize that parallel molecular adaptations are more prevalent between closely related species. We analyze the rate of parallel evolution in genome-size sets of orthologous genes in three groups of species with widely ranging levels of divergence: 46 species of the relatively recent lake Baikal amphipod radiation, a species flock of very closely related cichlids, and a set of significantly more divergent vertebrates. Strikingly, in genes of amphipods, the rate of parallel substitutions at nonsynonymous sites exceeded that at synonymous sites, suggesting rampant selection driving parallel adaptation. At sites of parallel substitutions, the intraspecies polymorphism is low, suggesting that parallelism has been driven by positive selection and is therefore adaptive. By contrast, in cichlids, the rate of nonsynonymous parallel evolution was similar to that at synonymous sites, whereas in vertebrates, this rate was lower than that at synonymous sites, indicating that in these groups of species, parallel substitutions are mainly fixed by drift.

Transcriptome analysis in Parhyale hawaiensis reveal sex-specific responses to AgNP and AgCl exposure

Analysis of the transcriptome of organisms exposed to toxicants offers new insights for ecotoxicology, but further research is needed to enhance interpretation of results and effectively incorporate them into useful environmental risk assessments. Factors that must be clarified to improve use of transcriptomics include assessment of the effect of organism sex within the context of toxicant exposure. Amphipods are well recognized as model organisms for toxicity evaluation because of their sensitivity and amenability to laboratory conditions. To investigate whether response to metals in crustaceans differs according to sex we analyzed the amphipod Parhyale hawaiensis after exposure to AgCl and Ag nanoparticles (AgNP) via contaminated food. Gene specific analysis and whole genome transcriptional profile of male and female organisms were performed by both RT-qPCR and RNA-seq. We observed that expression of transcripts of genes glutathione transferase (GST) did not differ among AgCl and AgNP treatments. Significant differences between males and females were observed after exposure to AgCl and AgNP. Males presented twice the number of differentially expressed genes in comparison to females, and more differentially expressed were observed after exposure to AgNP than AgCl treatments in both sexes. The genes that had the greatest change in expression relative to control were those genes related to peptidase and catalytic activity and chitin and carbohydrate metabolic processes. Our study is the first to demonstrate sex specific differences in the transcriptomes of amphipods upon exposure to toxicants and emphasizes the importance of considering gender in ecotoxicology.

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.

First draft genome for the sand-hopper Trinorchestia longiramus

Crustacean amphipods are important trophic links between primary producers and higher consumers. Although most amphipods occur in or around aquatic environments, the family Talitridae is the only family found in terrestrial and semi-terrestrial habitats. The sand-hopper Trinorchestia longiramus is a talitrid species often found in the sandy beaches of South Korea. In this study, we present the first draft genome assembly and annotation of this species. We generated ~380.3 Gb of sequencing data assembled in a 0.89 Gb draft genome. Annotation analysis estimated 26,080 protein-coding genes, with 89.9% genome completeness. Comparison with other amphipods showed that T. longiramus has 327 unique orthologous gene clusters, many of which are expanded gene families responsible for cellular transport of toxic substances, homeostatic processes, and ionic and osmotic stress tolerance. This first talitrid genome will be useful for further understanding the mechanisms of adaptation in terrestrial environments, the effects of heavy metal toxicity, as well as for studies of comparative genomic variation across amphipods.

Deep-sea amphipods around cobalt-rich ferromanganese crusts: Taxonomic diversity and selection of candidate species for connectivity analysis

The aim of this study was to select a candidate deep-sea amphipod species suitable for connectivity analyses in areas around cobalt-rich ferromanganese crusts (CRCs). We applied DNA barcoding based on the mitochondrial protein-coding gene, cytochrome c oxidase subunit I (COI), to specimens collected from the Xufu Guyot (the JA06 Seamount) off southeastern Minami-Torishima Island in the North Pacific, where CRCs are distributed. We used baited traps to collect 37 specimens. Comparison of COI sequences with public reference databases (GenBank, BOLD) showed that almost all of the specimens belonged to the superfamily Lysianassoidea, which is known to be ubiquitous in deep-sea areas. In a molecular taxonomic analysis of these sequences, we detected 11 clades. One of these clades (group 9) composed of 18 sequences and was identified by DNA barcoding as a putative species belonging to Abyssorchomene, which has been reported from the New Hebrides Trench in the South Pacific. We considered this species to be a candidate for connectivity analysis and analyzed its genome by restriction site–associated DNA sequencing. The results showed that the genetic variation in this species is adequate for analyzing connectivity patterns in CRC areas in the future.

DNA epigenetic marks are linked to embryo aberrations in amphipods

Linking exposure to environmental stress factors with diseases is crucial for proposing preventive and regulatory actions. Upon exposure to anthropogenic chemicals, covalent modifications on the genome can drive developmental and reproductive disorders in wild populations, with subsequent effects on the population persistence. Hence, screening of chemical modifications on DNA can be used to provide information on the probability of such disorders in populations of concern. Using a high-resolution mass spectrometry methodology, we identified DNA nucleoside adducts in gravid females of the Baltic amphipods Monoporeia affinis, and linked the adduct profiles to the frequency of embryo malformations in the broods. Twenty-three putative nucleoside adducts were detected in the females and their embryos, and eight modifications were structurally identified using high-resolution accurate mass data. To identify which adducts were significantly associated with embryo malformations, partial least squares regression (PLSR) modelling was applied. The PLSR model yielded three adducts as the key predictors: methylation at two different positions of the DNA (5-methyl-2'-deoxycytidine and N6-methyl-2'-deoxyadenosine) representing epigenetic marks, and a structurally unidentified nucleoside adduct. These adducts predicted the elevated frequency of the malformations with a high classification accuracy (84%). To the best of our knowledge, this is the first application of DNA adductomics for identification of contaminant-induced malformations in field-collected animals. The method can be adapted for a broad range of species and evolve as a new omics tool in environmental health assessment.

Impacts of ionising radiation on sperm quality, DNA integrity and post-fertilisation development in marine and freshwater crustaceans

Crustaceans have been designated as internationally important model organisms in the development of environmental radioprotection measures. Despite the known sensitivity of sperm to ionizing radiation, the impacts of chronic radiation exposure on male fertility in crustaceans have not been studied. For the first time, the present study aimed to assess the impacts of chronic radiation exposure on male fertility, sperm DNA damage and concomitant impacts on breeding in two amphipod crustaceans. Echinogammarus marinus and Gammarus pulex (male fertility only) were exposed to phosphorus-32 at dose rates of 0, 0.1, 1 and 10 mGy/d and sperm parameters, DNA damage and knock-on impacts on breeding were assessed. Sperm quality parameters and DNA damage were assessed using a fluorescent staining method and single cell gel electrophoresis respectively. Concomitant effects of male exposure to radiation on fecundity were determined by pairing phosphorus-32 exposed males to unexposed sexually mature females. In E. marinus, a statistically significant reduction of 9 and 11% in the quality of sperm was recorded at dose rates of 1 and 10 mGy/d respectively, with no significant effects recorded on sperm counts. Conversely in the freshwater G. pulex, no significant impact of radiation on sperm quantity or quality was recorded. For E. marinus, a statistically significant increase in DNA damage was recorded at doses of 10 mGy/d. Reduced fecundity and an increase in the frequency of abnormal embryos was recorded in female E. marinus breeding with males exposed to radiation. These findings suggest sperm quality may be a sensitive indicator of radiation exposure in invertebrates with potential impacts on the unexposed embryo, though unclear dose-response and differences between two closely related species necessitate further study before robust conclusions can be drawn.

Isolation and characterization of eleven novel microsatellite markers for fine-scale population genetic analyses of Gammarus pulex (Crustacea: Amphipoda)

The freshwater amphipod species Gammarus pulex (Linnaeus, 1758) is widespread across Europe and Asia and is able to live in a broad range of environmental conditions. Yet, it is of great interest to which degree it is able to tolerate and adapt to the current rapid anthropogenic environmental changes affecting its habitat, such as pollution, changes in river morphology, and invasions of alien gammarid species. Microsatellite primers for genetic population studies with G. pulex have been developed but due to the existence of several genetically different lineages within the species, the application of these primers is not always successful. In order to investigate the impacts of anthropogenic environmental changes on the spatio-genetic patterns of G. pulex lineage E in streams in the Saale river catchment in Germany, we designed eleven novel polymorphic microsatellites for this lineage using a high-throughput sequencing approach. These microsatellites enabled highly specific characterization of three closely related populations. The results show genetically distinct populations reflected by both a principal coordinates analysis and an analysis of molecular variance. Several of the newly designed microsatellite primers also enabled successful cross-amplification of the respective microsatellites in specimens of G. pulex lineage C, while only two microsatellites were amplified successfully and showed polymorphisms for all of the analyzed specimens of G. fossarum Koch, 1836. The microsatellites identified here are suitable for future assessments of micro-evolutionary dynamics of G. pulex from central Germany.

Genetic Diversity and Population Structure Analysis of Three Deep-Sea Amphipod Species from Geographically Isolated Hadal Trenches in the Pacific Ocean

Amphipods of the superfamily Lysianassoidea that inhabit the hadal zone ( > 6000 m) have large bathymetric ranges and play a key role in deep ocean ecosystems. The endemism of these amphipod species makes them a good model for investigating potent natural selection and restricted dispersal in deep ocean trenches. Here, we describe genetic diversity and intraspecific population differentiation among three amphipod species from four Pacific trenches based on a mtDNA concatenated dataset (CO Ι and 16S rRNA genes) from 150 amphipod individuals. All amphipod populations had low genetic diversity, as indicated by haplotype and nucleotide diversity values. Population geographic relationship analysis of two Alicella gigantea populations revealed no genetic differentiation between these two localities (pairwise genetic differentiation coefficient = 0.00032, gene flow = 784.58), and the major variation (99.97%) was derived from variation within the populations. Historical demographic events were investigated using Tajima's D and Fu's F neutrality tests and analysis of mismatch distribution. Consistent results provided strong evidence to support the premise that demographic expansion occurred only for the Mariana population of Hirondellea gigas, possibly within the last 2.1-3.4 million years. These findings suggest that the formation of amphipod population structure might be the result of multiple factors including high hydrostatic pressure, food distribution, trench topographic forcing and potential ecological interactions.

CO2 leakage simulation: Effects of the decreasing pH to the survival and reproduction of two crustacean species

The effects of CO₂-related acidification on two crustacean populations, the isopod Cyathura carinata and the amphipod Elasmopus rapax, were studied. Three pH levels were tested: artificial seawater without CO₂ injection and two levels of reduced pH. Even though RNA:DNA ratio was reduced for both species, no statistical significant differences were found between the control and the treatments. Both species experienced a reduction in survivorship, longevity and the body length of surviving animals; although the impairment observed in E. rapax was more severe than in C. carinata. The long life span isopod and the short life span amphipod experienced a high degree of impairment in the reproduction, likely due to the reallocation of resources from reproduction to body maintenance and increasing survival by postponing the brood production. Regardless of the underlying processes and the energetic pathways, both experienced failure to reproduce, which could lead to the local extinction of these species.

Characterization of the mitochondrial genome of an ancient amphipod Halice sp. MT-2017 (Pardaliscidae) from 10,908 m in the Mariana Trench

Small amphipods (Halice sp. MT-2017) with body length <1 cm were collected from the Challenger Deep (~10,920 m below sea level). The divergence time of their lineage was approximately 109 Mya, making this group ancient compared to others under study. The mitochondrial genome of Halice sp. shared the usual gene components of metazoans, comprising 13 protein coding genes (PCGs), 22 transfer RNAs (tRNAs), and 2 ribosomal RNAs (rRNAs). The arrangement of these genes, however, differed greatly from that of other amphipods. Of the 15 genes that were rearranged with respect to the pancrustacean gene pattern, 12 genes (2 PCGs, 2 rRNAs, and 8 tRNAs) were both translocated and strand-reversed. In contrast, the mitochondrial genomes in other amphipods never show so many reordered genes, and in most instances, only tRNAs were involved in strand-reversion-coupled translocation. Other characteristics, including reversed strand nucleotide composition bias, relatively higher composition of non-polar amino acids, and lower evolutionary rate, were also identified. Interestingly, the latter two features were shared with another hadal amphipod, Hirondellea gigas, suggesting their possible associations with the adaptation to deep-sea extreme habitats. Overall, our data provided a useful resource for future studies on the evolutionary and adaptive mechanisms of hadal faunas.

DNA metabarcoding unveils multiscale trophic variation in a widespread coastal opportunist

A thorough understanding of ecological networks relies on comprehensive information on trophic relationships among species. Since unpicking the diet of many organisms is unattainable using traditional morphology-based approaches, the application of high-throughput sequencing methods represents a rapid and powerful way forward. Here, we assessed the application of DNA metabarcoding with nearly universal primers for the mitochondrial marker cytochrome c oxidase I in defining the trophic ecology of adult brown shrimp, Crangon crangon, in six European estuaries. The exact trophic role of this abundant and widespread coastal benthic species is somewhat controversial, while information on geographical variation remains scant. Results revealed a highly opportunistic behaviour. Shrimp stomach contents contained hundreds of taxa (>1,000 molecular operational taxonomic units), of which 291 were identified as distinct species, belonging to 35 phyla. Only twenty ascertained species had a mean relative abundance of more than 0.5%. Predominant species included other abundant coastal and estuarine taxa, including the shore crab Carcinus maenas and the amphipod Corophium volutator. Jacobs' selectivity index estimates based on DNA extracted from both shrimp stomachs and sediment samples were used to assess the shrimp's trophic niche indicating a generalist diet, dominated by crustaceans, polychaetes and fish. Spatial variation in diet composition, at regional and local scales, confirmed the highly flexible nature of this trophic opportunist. Furthermore, the detection of a prevalent, possibly endoparasitic fungus (Purpureocillium lilacinum) in the shrimp's stomach demonstrates the wide range of questions that can be addressed using metabarcoding, towards a more robust reconstruction of ecological networks.

Morphology and life history divergence in cave and surface populations of Gammarus lacustris (L.)

Cave animals provide a unique opportunity to study contrasts in phenotype and life history in strikingly different environments when compared to surface populations, potentially related to natural selection. As such, we compared a permanent cave-living Gammarus lacustris (L.) population with two lake-resident surface populations analyzing morphology (eye- and antennal characters) and life-history (size at maturity, fecundity and egg-size). A part of the cytochrome c oxidase subunit I gene in the mitochondrion (COI) was analyzed to contrast genetic relationship of populations and was compared to sequences in GenBank to assess phylogeography and colonization scenarios. In the cave, a longer life cycle was implied, while surface populations seemed to have a shorter life cycle. Egg size, and size at maturity for both sexes, were larger in the cave than in surface populations, while fecundity was lower in the cave than in surface populations. The cave population had longer first- and second antennae with more articles, longer first- and second peduncles, and fewer ommatidia than surface populations. The cold low-productive cave environment may facilitate different phenotypic and life-history traits than in the warmer and more productive surface lake environments. The trait divergences among cave and surface populations resembles other cave-surface organism comparisons and may support a hypothesis of selection on sensory traits. The cave and Lake Ulvenvann populations grouped together with a sequence from Slovenia (comprising one genetic cluster), while Lake Lille Lauarvann grouped with a sequence from Ukraine (comprising another cluster), which are already recognized phylogenetic clusters. One evolutionary scenario is that the cave and surface populations were colonized postglacially around 9 000–10 000 years ago. We evaluate that an alternative scenario is that the cave was colonized during an interstadial during the last glaciation or earlier during the warm period before onset of the last glaciation.

Species identification and connectivity of marine amphipods in Canada's three oceans

Monitoring the distribution of marine biodiversity is a crucial step to better assess the impacts of global changes. Arctic marine fauna is dominated by amphipods in terms of abundance and biomass. These peracarids are an important marine order of crustaceans but the number of species found in the different Canadian oceans is currently unknown. Furthermore, most species are difficult to identify due to poor taxonomic descriptions and morphological convergence. The aim of this study was to assess the species diversity of marine amphipods in the three Canadian oceans using DNA barcoding. To do so, we produced a database of DNA barcodes of amphipods from the three Canadian Oceans publicly available from the BOLD website to which we added 310 new sequences from the Canadian Arctic Archipelago. We first delimited amphipod species based on barcode gap detection techniques and tree based method (bPTP) and then compared the composition of amphipods among the three oceans in order to assess the influence of past transarctic exchanges on Arctic diversity. Our analysis of 2309 sequences which represent more than 250 provisional species revealed a high connectivity between the Atlantic and Arctic Oceans. Our results also suggest that a single threshold to delimitate species is not suitable for amphipods. This study highlights the challenges involved in species delimitation and the need to obtain complete barcoding inventories in marine invertebrates.

Cryptic species in a well-known habitat: applying taxonomics to the amphipod genus Epimeria (Crustacea, Peracarida)

Taxonomy plays a central role in biological sciences. It provides a communication system for scientists as it aims to enable correct identification of the studied organisms. As a consequence, species descriptions should seek to include as much available information as possible at species level to follow an integrative concept of 'taxonomics'. Here, we describe the cryptic species Epimeria frankei sp. nov. from the North Sea, and also redescribe its sister species, Epimeria cornigera. The morphological information obtained is substantiated by DNA barcodes and complete nuclear 18S rRNA gene sequences. In addition, we provide, for the first time, full mitochondrial genome data as part of a metazoan species description for a holotype, as well as the neotype. This study represents the first successful implementation of the recently proposed concept of taxonomics, using data from high-throughput technologies for integrative taxonomic studies, allowing the highest level of confidence for both biodiversity and ecological research.

Identification, expression, and endocrine-disruption of three ecdysone-responsive genes in the sentinel species Gammarus fossarum

Taking advantage of a large transcriptomic dataset recently obtained in the sentinel crustacean amphipod Gammarus fossarum, we developed an approach based on sequence similarity and phylogenetic reconstruction to identify key players involved in the endocrine regulation of G. fossarum. Our work identified three genes of interest: the nuclear receptors RXR and E75, and the regulator broad-complex (BR). Their involvement in the regulation of molting and reproduction, along with their sensitivity to chemical contamination were experimentally assessed by studying gene expression during the female reproductive cycle, and after laboratory exposure to model endocrine disrupting compounds (EDCs): pyriproxyfen, tebufenozide and piperonyl butoxide. RXR expression suggested a role of this gene in ecdysis and post-molting processes. E75 presented two expression peaks that suggested a role in vitellogenesis, and molting. BR expression showed no variation during molting/reproductive cycle. After exposure to the three EDCs, a strong inhibition of the inter-molt E75 peak was observed with tebufenozide, and an induction of RXR after exposure to pyriproxyfen and piperonyl butoxide. These results confirm the implication of RXR and E75 in hormonal regulation of female reproductive cycles in G. fossarum and their sensitivity towards EDCs opens the possibility of using them as specific endocrine disruption biomarkers.

Application of cDNA-AFLP to biomarker exploration in a non-model species Grandidierella japonica

Biomarkers of exposure can be used to identify specific contaminants that are adversely affecting aquatic organisms. However, it remains prohibitively costly to investigate multiple novel biomarkers of exposure in a non-model species, despite the development of next-generation sequencing technology. In this study, we focused on the use of cDNA-amplified fragment length polymorphism (AFLP) as a cost-effective biomarker discovery tool to test whether it could identify biomarkers of exposure in the non-model amphipod species Grandidierella japonica. Loci were identified that were differentially expressed in amphipods exposed to reference chemicals (Cu, Zn, and nicotine) and to an environmental sample (road dust) at sublethal concentrations. Eight loci were shown to respond consistently to nicotine at different concentrations, but not to Cu or Zn. Some of the loci also responded to an environmental road dust sample containing nicotine. These findings suggest that loci identified using cDNA-AFLP could be used as biomarkers of nicotine exposure in environmental samples with complex matrices. Further studies with other organisms and toxicants are needed, but we have demonstrated that the use of cDNA-AFLP to identify biomarkers for ecotoxicological studies of non-model species is at least feasible.

Testing Local Adaptation in Five Populations of Hyalella azteca in Northern Alberta's Oil Sands Region

Canada's oil sands hold the third largest petroleum reserves worldwide and have experienced rapid economic growth. The oil sands region provides an ideal location for studying local adaptations through reciprocal transplant (RT) because populations within the region have been historically exposed to naturally occurring bitumen. Our objectives were to (1) determine if Hyalella azteca from habitats within the oil sands region exhibited increased tolerance to constituents associated with industrial bitumen extraction compared with H. azteca from habitats outside the region; and (2) determine if any observed tolerance was attributable to local adaptation. Five populations of H. azteca were reciprocally transplanted from reclaimed and reference wetlands: four from local wetlands plus one naïve laboratory population. Survival, toxicity, and behaviour were measured before and after the RT period. Survival varied by population and site. These results show that the differences in responses among populations are likely not attributable to local adaptation and that laboratory populations of H. azteca provide ecologically relevant results when tested in the field.

Evolution of mitochondrial genomes in Baikalian amphipods

BACKGROUND

Amphipods (Crustacea) of Lake Baikal are a very numerous and diverse group of invertebrates generally believed to have originated by adaptive radiation. The evolutionary history and phylogenetic relationships in Baikalian amphipods still remain poorly understood. Sequencing of mitochondrial genomes is a relatively feasible way for obtaining a set of gene sequences suitable for robust phylogenetic inferences. The architecture of mitochondrial genomes also may provide additional information on the mechanisms of evolution of amphipods in Lake Baikal.

RESULTS

Three complete and four nearly complete mitochondrial genomes of Baikalian amphipods were obtained by high-throughput sequencing using the Illumina platform. A phylogenetic inference based on the nucleotide sequences of all mitochondrial protein coding genes revealed the Baikalian species to be a monophyletic group relative to the nearest non-Baikalian species with a completely sequenced mitochondrial genome - Gammarus duebeni. The phylogeny of Baikalian amphipods also suggests that the shallow-water species Eulimnogammarus has likely evolved from a deep-water ancestor, however many other species have to be added to the analysis to test this hypothesis. The gene order in all mitochondrial genomes of studied Baikalian amphipods differs from the pancrustacean ground pattern. Mitochondrial genomes of four species possess 23 tRNA genes, and in three genomes the extra tRNA gene copies have likely undergone remolding. Widely varying lengths of putative control regions and other intergenic spacers are typical for the mitochondrial genomes of Baikalian amphipods.

CONCLUSIONS

The mitochondrial genomes of Baikalian amphipods display varying organization suggesting an intense rearrangement process during their evolution. Comparison of complete mitochondrial genomes is a potent approach for studying the amphipod evolution in Lake Baikal.

Implications of Cu and Ni toxicity in two members of the Hyalella azteca cryptic species complex: Mortality, growth, and bioaccumulation parameters

Hyalella azteca, an amphipod crustacean, is frequently used in freshwater toxicity tests. Since the mid-1980s, numerous organizations have collected and established cultures of H. azteca originating from localities across North America. However, H. azteca is actually a large cryptic species complex whose members satisfy both the biological and the phylogenetic species concepts. Genetic analysis at the mitochondrial COI gene has revealed that only 2 clades are cultured in 17 North American laboratories; however, there are 85 genetically divergent lineages within this complex in the wild. In the present study, 2 members (clades 1 and 8) of the H. azteca species complex were identified using the mitochondrial COI gene. These 2 clades were exposed to Cu or Ni for 14 d. A saturation-based mortality model and the general growth model were used to determine mortality (lethal concentration, 25% and 50% [LC25 and LC50], lethal body concentration, 25% and 50% [LBC25 and LBC50]) and growth (inhibitory concentration, 25% [IC25, IBC25]) endpoints, respectively. A modified saturation-based model was used to estimate metal bioaccumulation parameters. Clade 8 was significantly more tolerant than clade 1, with differences in LC50s. However, the effects of the metals on growth were not significantly different between clades, even though clade 1 was significantly larger than then clade 8. Differences in Cu or Ni bioaccumulation were not observed between clades 1 and 8. The differences in Cu and Ni LC50s may have implications for risk assessments, and it is recommended that toxicity experiments should only be performed with properly identified members of the H. azteca complex to maintain consistency among laboratories. Environ Toxicol Chem 2016;35:2817-2826. © 2016 SETAC.

Characterization, localization and temporal expression of crustacean hyperglycemic hormone (CHH) in the behaviorally rhythmic peracarid crustaceans, Eurydice pulchra (Leach) and Talitrus saltator (Montagu)

Crustacean hyperglycemic hormone (CHH) has been extensively studied in decapod crustaceans where it is known to exert pleiotropic effects, including regulation of blood glucose levels. Hyperglycemia in decapods seems to be temporally gated to coincide with periods of activity, under circadian clock control. Here, we used gene cloning, in situ hybridization and immunohistochemistry to describe the characterization and localization of CHH in two peracarid crustaceans, Eurydice pulchra and Talitrus saltator. We also exploited the robust behavioral rhythmicity of these species to test the hypothesis that CHH mRNA expression would resonate with their circatidal (12.4h) and circadian (24h) behavioral phenotypes. We show that both species express a single CHH transcript in the cerebral ganglia, encoding peptides featuring all expected, conserved characteristics of other CHHs. E. pulchra preproCHH is an amidated 73 amino acid peptide N-terminally flanked by a short, 18 amino acid precursor related peptide (CPRP) whilst the T. saltator prohormone is also amidated but 72 amino acids in length and has a 56 residue CPRP. The localization of both was mapped by immunohistochemistry to the protocerebrum with axon tracts leading to the sinus gland and into the tritocerebrum, with striking similarities to terrestrial isopod species. We substantiated the cellular position of CHH immunoreactive cells by in situ hybridization. Although both species showed robust activity rhythms, neither exhibited rhythmic transcriptional activity indicating that CHH transcription is not likely to be under clock control. These data make a contribution to the inventory of CHHs that is currently lacking for non-decapod species.

Ecotoxicogenomics: Microarray interlaboratory comparability

Transcriptomic analysis can complement traditional ecotoxicology data by providing mechanistic insight, and by identifying sub-lethal organismal responses and contaminant classes underlying observed toxicity. Before transcriptomic information can be used in monitoring and risk assessment, it is necessary to determine its reproducibility and detect key steps impacting the reliable identification of differentially expressed genes. A custom 15K-probe microarray was used to conduct transcriptomics analyses across six laboratories with estuarine amphipods exposed to cyfluthrin-spiked or control sediments (10 days). Two sample types were generated, one consisted of total RNA extracts (Ex) from exposed and control samples (extracted by one laboratory) and the other consisted of exposed and control whole body amphipods (WB) from which each laboratory extracted RNA. Our findings indicate that gene expression microarray results are repeatable. Differentially expressed data had a higher degree of repeatability across all laboratories in samples with similar RNA quality (Ex) when compared to WB samples with more variable RNA quality. Despite such variability a subset of genes were consistently identified as differentially expressed across all laboratories and sample types. We found that the differences among the individual laboratory results can be attributed to several factors including RNA quality and technical expertise, but the overall results can be improved by following consistent protocols and with appropriate training.

Challenges for using quantitative PCR test batteries as a TIE-type approach to identify metal exposure in benthic invertebrates

The epibenthic amphipod Melita plumulosa shows unique gene expression profiles when exposed to different contaminants. We hypothesized that specific changes in transcript abundance could be used in a battery of quantitative polymerase chain reaction (qPCR) assays as a toxicity identification evaluation (TIE)-like approach to identify the most relevant stressor in field-contaminated sediments. To test this hypothesis, seven candidate transcriptomic markers were selected, and their specificity following metal exposure was confirmed. The performance of these markers across different levels of added metals was verified. The ability of these transcripts to act as markers was tested by exposing amphipods to metal-contaminated field-collected sediments and measuring changes in transcript abundance via qPCR. For two of the three sediments tested, at least some of the transcriptomic patterns matched our predictions, suggesting that they would be effective in helping to identify metal exposure in field sediments. However, following exposure to the third sediment, transcriptomic patterns were unlike our predictions. These results suggest that the seven transcripts may be insufficient to discern individual contaminants from complex mixtures and that microarray or RNA-Seq global gene expression profiles may be more effective for TIE. Changes in transcriptomics based on laboratory exposures to single compounds should be carefully validated before the results are used to analyze mixtures.

Genomic resources notes accepted 1 October 2014-30 November 2014

This article documents the public availability of (i) RAD sequencing data and validated SNPs for the American mink Neovison vison and (ii) Transcriptome resources for two nonmodel freshwater crustacean species, the copepod Eucyclops serrulatus and the amphipod Echinogammarus veneris.