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

Proteogenomic reconstruction of organ-specific metabolic networks in an environmental sentinel species, the amphipod Gammarus fossarum

Metabolic pathways are affected by the impacts of environmental contaminants underlying a large variability of toxic effects across different species. However, the systematic reconstruction of metabolic pathways remains limited in environmental sentinel species due to the lack of available genomic data in many taxa of animal diversity. In this study we used a multi-omics approach to reconstruct the most comprehensive map of metabolic pathways for a crustacean model in biomonitoring, the amphipod Gammarus fossarum in order to improve the knowledge of the metabolism of this sentinel species. We revisited the assembly of RNA-seq data by de novo approaches to reduce RNA contaminants and transcript redundancy. We also acquired extensive mass spectrometry shotgun proteomic data on several organs from a reference population of G. fossarum males and females to identify organ-specific metabolic profiles. The G. fossarum metabolic pathway reconstruction (available through the metabolic database GamfoCyc) was performed by adapting the genomic tool CycADS and we identified 377 pathways representing 7630 annotated enzymes, 2610 enzymatic reactions and the expression of 858 enzymes was experimentally validated by proteomics. To our knowledge, our analysis provides for the first time a systematic metabolic pathway reconstruction and the proteome profiles of these pathways at the organ level in this sentinel species. As an example, we show an elevated abundance in enzymes involved in ATP biosynthesis and fatty acid beta-oxidation indicative of the high-energy requirement of the gills, or the key anabolic and detoxification role of the hepatopancreatic caeca, as exemplified by the specific expression of the retinoid biosynthetic pathways and glutathione synthesis. In conclusion, the multi-omics data integration performed in this study provides new resources to investigate metabolic processes in crustacean amphipods and their role in mediating the effects of environmental contaminant exposures in sentinel species. SYNOPSIS: This study provide the first evidence that it is possible to combine multiple omics data to exhaustively describe the metabolic network of a model species in ecotoxicology, Gammarus fossarum, for which a reference genome is not yet available.

A broad-taxa approach as an important concept in ecotoxicological studies and pollution monitoring

Aquatic invertebrates play a pivotal role in (eco)toxicological assessments because they offer ethical, cost-effective and repeatable testing options. Additionally, their significance in the food chain and their ability to represent diverse aquatic ecosystems make them valuable subjects for (eco)toxicological studies. To ensure consistency and comparability across studies, international (eco)toxicology guidelines have been used to establish standardised methods and protocols for data collection, analysis and interpretation. However, the current standardised protocols primarily focus on a limited number of aquatic invertebrate species, mainly from Arthropoda, Mollusca and Annelida. These protocols are suitable for basic toxicity screening, effectively assessing the immediate and severe effects of toxic substances on organisms. For more comprehensive and ecologically relevant assessments, particularly those addressing long-term effects and ecosystem-wide impacts, we recommended the use of a broader diversity of species, since the present choice of taxa exacerbates the limited scope of basic ecotoxicological studies. This review provides a comprehensive overview of (eco)toxicological studies, focusing on major aquatic invertebrate taxa and how they are used to assess the impact of chemicals in diverse aquatic environments. The present work supports the use of a broad-taxa approach in basic environmental assessments, as it better represents the natural populations inhabiting various ecosystems. Advances in omics and other biochemical and computational techniques make the broad-taxa approach more feasible, enabling mechanistic studies on non-model organisms. By combining these approaches with in vitro techniques together with the broad-taxa approach, researchers can gain insights into less-explored impacts of pollution, such as changes in population diversity, the development of tolerance and transgenerational inheritance of pollution responses, the impact on organism phenotypic plasticity, biological invasion outcomes, social behaviour changes, metabolome changes, regeneration phenomena, disease susceptibility and tissue pathologies. This review also emphasises the need for harmonised data-reporting standards and minimum annotation checklists to ensure that research results are findable, accessible, interoperable and reusable (FAIR), maximising the use and reusability of data. The ultimate goal is to encourage integrated and holistic problem-focused collaboration between diverse scientific disciplines, international standardisation organisations and decision-making bodies, with a focus on transdisciplinary knowledge co-production for the One-Health approach.

Microplastics in river sediment: Chronic exposure of the amphipod Gammarus fossarum to polyethylene terephthalate in a microcosm

Microplastics constitute a form of particulate matter in aquatic environments, where they are a widespread pollutant. The broad range of particle sizes facilitates interactions with diverse species assemblages. Exposure to microplastics can negatively impact organisms, but similar effects also arise from exposure to naturally occurring particles, such as increased oxidative stress. It therefore remains uncertain, what effects are specific to microplastic particles, and how these effects manifest as a consequence of chronic exposure. Here we show in microcosm experiments that long-term exposure (111 days) to irregularly shaped polyethylene terephthalate (PET) fragments (10-400 µm) added to riverine sediments did not negatively impact the amphipod Gammarus fossarum's group size, and oxygen consumption, and minimally affected proteome composition. We found that these results were consistent for male and female specimens when exposed to an environmentally relevant concentration (0.004% of sediment dry weight; dw) and an environmentally less realistic one (4% dw). In female specimens' whole proteomes, we identified two highly differentially abundant proteins, which have been associated with an organism's response to xenobiotics. We conclude that in this sentinel species exposure to PET microplastic fragments mixed into the sediment does not elicit significant stress, even at concentrations exceeding current exposure levels in the environment.

Examination of gammarid transcriptomes reveals a widespread occurrence of key metabolic genes from epibiont bdelloid rotifers in freshwater species

Previous data revealed the unexpected presence of genes encoding for long-chain polyunsaturated fatty acid (LC-PUFA) biosynthetic enzymes in transcriptomes from freshwater gammarids but not in marine species, even though closely related species were compared. This study aimed to clarify the origin and occurrence of selected LC-PUFA biosynthesis gene markers across all published gammarid transcriptomes. Through systematic searches, we confirmed the widespread occurrence of sequences from seven elongases and desaturases involved in LC-PUFA biosynthesis, in transcriptomes from freshwater gammarids but not marine species, and clarified that such occurrence is independent from the gammarid species and geographical origin. The phylogenetic analysis established that the retrieved elongase and desaturase sequences were closely related to bdelloid rotifers, confirming that multiple transcriptomes from freshwater gammarids contain contaminating rotifers' genetic material. Using the Adineta steineri genome, we investigated the genomic location and exon-intron organization of the elongase and desaturase genes, establishing they are all genome-anchored and, importantly, identifying instances of horizontal gene transfer. Finally, we provide compelling evidence demonstrating Bdelloidea desaturases and elongases enable these organisms to perform all the reactions for de novo biosynthesis of PUFA and, from them, LC-PUFA, an advantageous trait when considering the low abundance of these essential nutrients in freshwater environments.

Effects of Diclofenac on the Reproductive Health, Respiratory Rate, Cardiac Activity, and Heat Tolerance of Aquatic Animals

Diclofenac is an important pharmaceutical present in the water cycle of wastewater treatment and one of the most distributed drugs in aquatic ecosystems. Despite the great interest in the fate of diclofenac in freshwaters, the effects of environmentally relevant concentrations on invertebrates are still unclear. Two species of freshwater invertebrates, the amphipod Gmelinoides fasciatus and the bivalve mollusk Unio pictorum, were exposed to diclofenac concentrations of 0.001-2 μg/L (environmentally relevant levels) for 96 h. A set of biological endpoints (survival, fecundity, embryo abnormalities, respiration and heart rates, heat tolerance, and cardiac stress tolerance) were estimated in exposed invertebrates. Effects of diclofenac on amphipod metabolic rate and reproduction (number and state of embryos) and adaptive capacity (cardiac stress tolerance) in both species were evident. The oxygen consumption of amphipods exposed to diclofenac of 0.1-2 μg/L was 1.5-2 times higher than in the control, indicating increased energy requirements for standard metabolism in the presence of diclofenac (>0.1 μg/L). The heart rate recovery time in mollusks after heating to critical temperature (30 °C) was 1.7 and 9 times greater in mollusks exposed to 0.1 and 0.9 μg/L, respectively, than in the control (24 min). A level of diclofenac >0.9 μg/L adversely affected amphipod embryos, leading to an increase in the number of embryos with impaired development, which subsequently died. Thus, the lowest effective concentration of diclofenac (0.1 μg/L) led to increased energy demands of animals while reducing cardiac stress tolerance, and at a level close to 1 μg/L reproductive disorders (elevated mortality of the embryos) occurred. Environ Toxicol Chem 2022;41:677-686. © 2021 SETAC.

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

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

Trends in the Application of "Omics" to Ecotoxicology and Stress Ecology

Our ability to predict and assess how environmental changes such as pollution and climate change affect components of the Earth's biome is of paramount importance. This need positioned the fields of ecotoxicology and stress ecology at the center of environmental monitoring efforts. Advances in these interdisciplinary fields depend not only on conceptual leaps but also on technological advances and data integration. High-throughput "omics" technologies enabled the measurement of molecular changes at virtually all levels of an organism's biological organization and thus continue to influence how the impacts of stressors are understood. This bibliometric review describes literature trends (2000-2020) that indicate that more different stressors than species are studied each year but that only a few stressors have been studied in more than two phyla. At the same time, the molecular responses of a diverse set of non-model species have been investigated, but cross-species comparisons are still rare. While transcriptomics studies dominated until 2016, a shift towards proteomics and multiomics studies is apparent. There is now a wealth of data at functional omics levels from many phylogenetically diverse species. This review, therefore, addresses the question of how to integrate omics information across species.

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

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

Transcriptomic data on the transgenerational exposure of the keystone amphipod Gammarus locusta to simvastatin

The use of transcriptomics data brings new insights and works as a powerful tool to explore the molecular mode of action (MoA) of transgenerational inheritance effects of contaminants of emerging concern. Therefore, in this dataset, we present the transcriptomic data of the transgenerational effects of environmentally relevant simvastatin levels, one of the most prescribed human pharmaceuticals, in the keystone amphipod species Gammarus locusta. In summary, G. locusta juveniles were maintained under simvastatin exposure up to adulthood (exposed group - F0E) and the offspring of F0E were transferred to control water for the three subsequent generations (transgenerational group - F1T, F2T and F3T). To gain insights into the biological functions and canonical pathways transgenerationally disrupted by simvastatin, a G. locusta de novo transcriptome assembly was produced and the transcriptomic profiles of three individual G. locusta females, per group, over the four generations (F0 to F3) - solvent control groups (F0.C, F1.C, F2.C and F3.C), F0 320 ng/L simvastatin exposed group (F0.320E) and F1 to F3 320 transgenerational group (F1.320T; F2.320T and F3.320T) - were analyzed. Briefly, Illumina HiSeq™ 2500 platform was used to perform RNA sequencing, and due to the unavailability of G. locusta genome, the RNA-seq datasets were assembled de novo using Trinity and annotated with Trinotate software. After assembly and post-processing steps, 106093 transcripts with N50 of 2371 bp and mean sequence length of 1343.98 bp was produced. BUSCO analyses showed a transcriptome with gene completeness of 97.5 % Arthropoda library profile. The Bowtie2, RSEM and edgeR tools were used for the differential gene expression (DEGs) analyses that allowed the identification of a high quantity of genes differentially expressed in all generations. Finally, to identify the main metabolic pathways affected by the transgenerational effects of SIM across all generations, the DGEs genes were blasted onto KEGG pathways database using the KAAS webserver. The data furnished in this article allows a better molecular understanding of the transgenerational effects produced by simvastatin in the keystone amphipod G. locusta and has major implications for hazard and risk assessment of pharmaceuticals and other emerging contaminants. This article is related to the research article entitled "Transgenerational inheritance of chemical-induced signature: a case study with simvastatin [1].