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.