Transcriptomic responses of the endangered freshwater mussel Margaritifera margaritifera to trace metal contamination in the Dronne River, France

Transcriptome sequencing of the endangered land snail Karaftohelix adamsi from the Island Ulleung: De novo assembly, annotation, valuation of fitness genes and SSR markers

BACKGROUND

The Bradybaenidae snail Karaftohelix adamsi is endemic to Korea, with the species tracked from Island Ulleung in North Gyeongsang Province of South Korea. K. adamsi has been classified under the Endangered Wildlife Class II species of Korea and poses a severe risk of extinction following habitat disturbances. With no available information at the DNA (genome) or mRNA (transcriptome) level for the species, conservation by utilizing informed molecular resources seems difficult.

OBJECTIVE

In this study, we used the Illumina short-read sequencing and Trinity de novo assembly to draft the reference transcriptome of K. adamsi.

RESULTS

After assembly, 13,753 unigenes were obtained of which 10,511 were annotated to public databases (a maximum of 10,165 unigenes found homologs in PANM DB). A total of 6,351, 3,535, 358, and 3,407 unigenes were ascribed to the functional categories under KOG, GO, KEGG, and IPS, respectively. The transcripts such as the HSP 70, aquaporin, TLR, and MAPK, among others, were screened as putative functional resources for adaptation. DNA transposons were found to be thickly populated in comparison to retrotransposons in the assembled unigenes. Further, 2,164 SSRs were screened with the promiscuous presence of dinucleotide repeats such as AC/GT and AG/CT.

CONCLUSION

The transcriptome-guided discovery of molecular resources in K. adamsi will not only serve as a basis for functional genomics studies but also provide sustainable tools to be utilized for the protection of the species in the wild. Moreover, the development of polymorphic SSRs is valuable for the identification of species from newer habitats and cross-species genotyping.

The transcriptome assembly of the European freshwater mussel Unio elongatulus C. Pfeiffer, 1825

Freshwater mussels of the order Unionida are a global conservation concern. Species of this group are strictly freshwater, sessile, slow-growing animals and, extremely sensitive to environmental changes. Human-mediated changes in freshwater habitats are imposing enormous pressure on the survival of freshwater mussels. Although a few flagship species are protected in Europe, other highly imperilled species receive much less attention. Moreover, knowledge about biology, ecology, and evolution and proper conservation assessments of many European species are still sparse. This knowledge gap is further aggravated by the lack of genomic resources available, which are key tools for conservation. Here we present the transcriptome assembly of Unio elongatulus C. Pfeiffer, 1825, one of the least studied European freshwater mussels. Using the individual sequencing outputs from eight physiologically representative mussel tissues, we provide an annotated panel of tissue-specific Relative Gene Expression profiles. These resources are pivotal to studying the species' biological and ecological features, as well as helping to understand its vulnerability to current and future threats.

Pathology and infectious agents of unionid mussels: A primer for pathologists in disease surveillance and investigation of mortality events

Freshwater mussels are one of the most imperiled groups of organisms in the world, and more than 30 species have gone extinct in the last century. While habitat alteration and destruction have contributed to the declines, the role of disease in mortality events is unclear. In an effort to involve veterinary pathologists in disease surveillance and the investigation of freshwater mussel mortality events, we provide information on the conservation status of unionids, sample collection and processing techniques, and unique and confounding anatomical and physiological differences. We review the published accounts of pathology and infectious agents described in freshwater mussels including neoplasms, viruses, bacteria, fungi, fungal-like agents, ciliated protists, Aspidogastrea, Digenea, Nematoda, Acari, Diptera, and Odonata. Of the identified infectious agents, a single viral disease, Hyriopsis cumingii plague disease, that occurs only in cultured mussels is known to cause high mortality. Parasites including ciliates, trematodes, nematodes, mites, and insects may decrease host fitness, but are not known to cause mortality. Many of the published reports identify infectious agents at the light or ultrastructural microscopy level with no lesion or molecular characterization. Although metagenomic analyses provide sequence information for infectious agents, studies often fail to link the agents to tissue changes at the light or ultrastructural level or confirm their role in disease. Pathologists can bridge this gap between identification of infectious agents and confirmation of disease, participate in disease surveillance to ensure successful propagation programs necessary to restore decimated populations, and investigate mussel mortality events to document pathology and identify causality.

A review of the endangered mollusks transcriptome under the threatened species initiative of Korea

Transcriptome studies for conservation of endangered mollusks is a proactive approach towards managing threats and uncertainties facing these species in natural environments. The population of these species is declining due to habitat destruction, illicit wildlife trade, and global climate change. These activities risk the free movement of species across the wild landscape, loss of breeding grounds, and restrictions in displaying the physiological attributes so crucial for faunal welfare. Gastropods face the most negative ecological effects and have been enlisted under Korea's protective species consortium based on their population dynamics in the last few years. Moreover, with the genetic resources restricted for such species, conservation by informed planning is not possible. This review provides insights into the activities under the threatened species initiative of Korea with special reference to the transcriptome assemblies of endangered mollusks. The gastropods such as Ellobium chinense, Aegista chejuensis, Aegista quelpartensis, Incilaria fruhstorferi, Koreanohadra kurodana, Satsuma myomphala, and Clithon retropictus have been represented. Moreover, the transcriptome summary of bivalve Cristaria plicata and Caenogastropoda Charonia lampas sauliae is also discussed. Sequencing, de novo assembly, and annotation identified transcripts or homologs for the species and, based on an understanding of the biochemical and molecular pathways, were ascribed to predictive gene function. Mining for simple sequence repeats from the transcriptome have successfully assisted genetic polymorphism studies. A comparison of the transcriptome scheme of Korean endangered mollusks with the genomic resources of other endangered mollusks have been discussed with homologies and analogies for dictating future research.

PacBio Hi-Fi genome assembly of the Iberian dolphin freshwater mussel Unio delphinus Spengler, 1793

Mussels of order Unionida are a group of strictly freshwater bivalves with nearly 1,000 described species widely dispersed across world freshwater ecosystems. They are highly threatened showing the highest record of extinction events within faunal taxa. Conservation is particularly concerning in species occurring in the Mediterranean biodiversity hotspot that are exposed to multiple anthropogenic threats, possibly acting in synergy. That is the case of the dolphin freshwater mussel Unio delphinus Spengler, 1793, endemic to the western Iberian Peninsula with recently strong population declines. To date, only four genome assemblies are available for the order Unionida and only one European species. We present the first genome assembly of Unio delphinus. We used the PacBio HiFi to generate a highly contiguous genome assembly. The assembly is 2.5 Gb long, possessing 1254 contigs with a contig N50 length of 10 Mbp. This is the most contiguous freshwater mussel genome assembly to date and is an essential resource for investigating the species' biology and evolutionary history that ultimately will help to support conservation strategies.

Harnessing the omics revolution to address the global biodiversity crisis

Human disturbances are altering global biodiversity in unprecedented ways. We identify three fundamental challenges underpinning our understanding of global biodiversity (namely discovery, loss, and preservation), and discuss how the omics revolution (e.g. genomics, transcriptomics, proteomics, metabolomics, and meta-omics) can help address these challenges. We also discuss how omics tools can illuminate the major drivers of biodiversity loss, including invasive species, pollution, urbanization, overexploitation, and climate change, with a special focus on highly diverse tropical environments. Although omics tools are transforming the traditional toolkit of biodiversity research, their application to addressing the current biodiversity crisis remains limited and may not suffice to offset current rates of biodiversity loss. Despite technical and logistical challenges, omics tools need to be fully integrated into global biodiversity research, and better strategies are needed to improve their translation into biodiversity policy and practice. It is also important to recognize that although the omics revolution can be considered the biologist's dream, socioeconomic disparity limits their application in biodiversity research.

The gill transcriptome of threatened European freshwater mussels

Genomic tools applied to non-model organisms are critical to design successful conservation strategies of particularly threatened groups. Freshwater mussels of the Unionida order are among the most vulnerable taxa and yet almost no genetic resources are available. Here, we present the gill transcriptomes of five European freshwater mussels with high conservation concern: Margaritifera margaritifera, Unio crassus, Unio pictorum, Unio mancus and Unio delphinus. The final assemblies, with N50 values ranging from 1069–1895 bp and total BUSCO scores above 90% (Eukaryote and Metazoan databases), were structurally and functionally annotated, and made available. The transcriptomes here produced represent a valuable resource for future studies on these species’ biology and ultimately guide their conservation.

Measurement(s)	transcriptomics
Technology Type(s)	Illumina sequencing
Sample Characteristic - Organism	Margaritifera margaritifera • Unio crassus • Unio delphinus • Unio mancus • Unio pictorum
Sample Characteristic - Location	Europe

Evaluation of transcriptional biomarkers using a high-resolution regression approach: Concentration-dependence of selected transcripts in copper-exposed freshwater mussels (Anodonta anatina)

We tested concentration-dependence of selected gene transcripts (cat, gst, hsp70, hsp90, mt and sod) for evaluation as biomarkers of chemical stress. Contrary to the common approach of factorial designs and few exposure concentrations, we used regression across a high-resolution concentration series. Specifically, freshwater mussels (Anodonta anatina) were acutely (96 h) exposed to Cu (13 nominal concentrations, measuring 0.13-1 600 µg/L), and transcripts were measured by RT-qPCR. In digestive glands, cat, hsp90 and mt decreased with water Cu (p < 0.05), but response magnitudes saturated at < 2-fold decreases. In gills, gst, hsp70, hsp90 and mt increased with water Cu (p < 0.05). While hsp70, hsp90 and mt exceeded 2-fold increases within the exposure range, high Cu concentrations were required (38-160 µg/L). Although gill responses were generally more robust compared to digestive glands, overall small response magnitudes and moderate sensitivity may set limit for potential application as general biomarkers of chemical stress.

The Crown Pearl: a draft genome assembly of the European freshwater pearl mussel Margaritifera margaritifera (Linnaeus, 1758)

Since historical times, the inherent human fascination with pearls turned the freshwater pearl mussel Margaritifera margaritifera (Linnaeus, 1758) into a highly valuable cultural and economic resource. Although pearl harvesting in M. margaritifera is nowadays residual, other human threats have aggravated the species conservation status, especially in Europe. This mussel presents a myriad of rare biological features, e.g. high longevity coupled with low senescence and Doubly Uniparental Inheritance of mitochondrial DNA, for which the underlying molecular mechanisms are poorly known. Here, the first draft genome assembly of M. margaritifera was produced using a combination of Illumina Paired-end and Mate-pair approaches. The genome assembly was 2.4 Gb long, possessing 105,185 scaffolds and a scaffold N50 length of 288,726 bp. The ab initio gene prediction allowed the identification of 35,119 protein-coding genes. This genome represents an essential resource for studying this species’ unique biological and evolutionary features and ultimately will help to develop new tools to promote its conservation.

A Commentary on the Use of Bivalve Mollusks in Monitoring Metal Pollution Levels

The objective of this commentary is to promote the use of bivalves as biomonitors, which is a part of the continual efforts of the International Mussel Watch. This commentary is an additional discussion on "Bivalve mollusks in metal pollution studies: From bioaccumulation to biomonitoring" by Zuykov et al., published in Chemosphere 93, 201-208. The present discussion can serve as a platform for further insights to provide new thoughts and novel ideas on how to make better use of bivalves in biomonitoring studies. The certainty of better and more extensive applications of mollusks in environmental monitoring in the future is almost confirmed but more studies are urgently needed. With all the reported studies using bivalves as biomonitors of heavy metal pollution, the effectiveness of using Mussel Watch is beyond any reasonable doubts. The challenge is the development of more accurate methodologies for of heavy metal data interpretation, and the precision of the biomonitoring studies using bivalves as biomonitors, whether in coastal or freshwater ecosystems. Lastly, inclusion of human health risk assessment of heavy metals in commercial bivalves would make the research papers of high public interest.

A High-Quality Reference Genome for a Parasitic Bivalve with Doubly Uniparental Inheritance (Bivalvia: Unionida)

From a genomics perspective, bivalves (Mollusca: Bivalvia) have been poorly explored with the exception for those of high economic value. The bivalve order Unionida, or freshwater mussels, has been of interest in recent genomic studies due to their unique mitochondrial biology and peculiar life cycle. However, genomic studies have been hindered by the lack of a high-quality reference genome. Here, I present a genome assembly of Potamilus streckersoni using Pacific Bioscience single-molecule real-time long reads and 10X Genomics-linked read sequencing. Further, I use RNA sequencing from multiple tissue types and life stages to annotate the reference genome. The final assembly was far superior to any previously published freshwater mussel genome and was represented by 2,368 scaffolds (2,472 contigs) and 1,776,755,624 bp, with a scaffold N50 of 2,051,244 bp. A high proportion of the assembly was comprised of repetitive elements (51.03%), aligning with genomic characteristics of other bivalves. The functional annotation returned 52,407 gene models (41,065 protein, 11,342 tRNAs), which was concordant with the estimated number of genes in other freshwater mussel species. This genetic resource, along with future studies developing high-quality genome assemblies and annotations, will be integral toward unraveling the genomic bases of ecologically and evolutionarily important traits in this hyper-diverse group.

Immunological Responses of Marine Bivalves to Contaminant Exposure: Contribution of the -Omics Approach

The increasing number of data studies on the biological impact of anthropogenic chemicals in the marine environment, together with the great development of invertebrate immunology, has identified marine bivalves as a key invertebrate group for studies on immunological responses to pollutant exposure. Available data on the effects of contaminants on bivalve immunity, evaluated with different functional and molecular endpoints, underline that individual functional parameters (cellular or humoral) and the expression of selected immune-related genes can distinctly react to different chemicals depending on the conditions of exposure. Therefore, the measurement of a suite of immune biomarkers in hemocytes and hemolymph is needed for the correct evaluation of the overall impact of contaminant exposure on the organism's immunocompetence. Recent advances in -omics technologies are revealing the complexity of the molecular players in the immune response of different bivalve species. Although different -omics represent extremely powerful tools in understanding the impact of pollutants on a key physiological function such as immune defense, the -omics approach has only been utilized in this area of investigation in the last few years. In this work, available information obtained from the application of -omics to evaluate the effects of pollutants on bivalve immunity is summarized. The data shows that the overall knowledge on this subject is still quite limited and that to understand the environmental relevance of any change in immune homeostasis induced by exposure to contaminants, a combination of both functional assays and cutting-edge technology (transcriptomics, proteomics, and metabolomics) is required. In addition, the utilization of metagenomics may explain how the complex interplay between the immune system of bivalves and its associated bacterial communities can be modulated by pollutants, and how this may in turn affect homeostatic processes of the host, host–pathogen interactions, and the increased susceptibility to disease. Integrating different approaches will contribute to knowledge on the mechanism responsible for immune dysfunction induced by pollutants in ecologically and economically relevant bivalve species and further explain their sensitivity to multiple stressors, thus resulting in health or disease.

Transfer and Transcriptomic Profiling in Liver and Brain of European Eels (Anguilla anguilla) After Diet-borne Exposure to Gold Nanoparticles

A nanometric revolution is underway, promising technical innovations in a wide range of applications and leading to a potential boost in environmental discharges. The propensity of nanoparticles (NPs) to be transferred throughout trophic chains and to generate toxicity was mainly assessed in primary consumers, whereas a lack of knowledge for higher trophic levels persists. The present study focused on a predatory fish, the European eel (Anguilla anguilla) exposed to gold NPs (AuNPs; 10 nm, polyethylene glycol-coated) for 21 d at 3 concentration levels in food: 0 (NP0), 1 (NP1), and 10 (NP10) mg Au kg^-1 . Transfer was assessed by Au quantification in eel tissues, and transcriptomic responses in the liver and brain were revealed by a high-throughput RNA-sequencing approach. Eels fed at NP10 presented an erratic feeding behavior, whereas Au quantification only indicated transfer to intestine and kidney of NP1-exposed eels. Sequencing of RNA was performed in NP0 and NP1 eels. A total of 258 genes and 156 genes were significantly differentially transcribed in response to AuNP trophic exposure in the liver and brain, respectively. Enrichment analysis highlighted modifications in the immune system-related processes in the liver. In addition, results pointed out a shared response of both organs regarding 13 genes, most of them being involved in immune functions. This finding may shed light on the mode of action and toxicity of AuNPs in fish. Environ Toxicol Chem 2020;39:2450-2461. © 2020 SETAC.

Omics approaches for conservation biology research on the bivalve Chamelea gallina

The striped venus (Chamelea gallina) is an important economic resource in the Mediterranean Basin; this species has exhibited a strong quantitative decline in the Adriatic Sea. The aim of this work was to provide a comprehensive view of the biological status of C. gallina to elucidate the bioecological characteristics and genetic diversity of wild populations. To the best of our knowledge, this investigation is the first to perform a multidisciplinary study on C. gallina based on two omics approaches integrated with histological, ecotoxicological, and chemical analyses and with the assessment of environmental parameters. The results obtained through RNA sequencing indicated that the striped venus has a notable ability to adapt to different environmental conditions. Moreover, the stock reduction exhibited by this species in the last 2 decades seems not to have negatively affected its genetic diversity. Indeed, the high level of genetic diversity that emerged from our ddRAD dataset analyses is ascribable to the high larval dispersal rate, which might have played a “compensatory role” on local fluctuations, conferring to this species a good adaptive potential to face the environmental perturbations. These findings may facilitate the efforts of conservation biologists to adopt ad hoc management plans for this fishery resource.

Disentangling of the ecotoxicological signal using "omics" analyses, a lesson from the survey of the impact of cyanobacterial proliferations on fishes

Omics technologies offer unprecedented perspectives for the rational investigation of complex biological systems. Indeed, omics present the ability of offering an extensive perception of the biochemistry and physiology of the cell and of any perturbing consequences of contaminants through the joint investigation of thousands of molecular responses simultaneously; then it has recently conducted to a fervent attention by research ecotoxicologists. Beyond the presentation of latest advances, exemplified here by omics investigation of cyanobacterial deleterious effects on various fishes (at various experimental and biological scales and with various analytical tools and pipeline), the present review paper re-explores the promising perspectives and also the pitfalls of such holistic investigations of the ecotoxicological response of organisms for environmental assessment.

Acute Toxicity of Sodium Chloride, Nitrates, Ortho-Phosphates, Cadmium, Arsenic and Aluminum for Juveniles of the Freshwater Pearl Mussel: Margaritifera Margaritifera (L.1758)

The freshwater pearl mussel (FWPM) Margaritifera margaritifera (L.1758) is critically endangered in Europe and ecotoxicological studies on the species are scares. Here, 96 h acute toxicity tests were conducted at 16 °C with sodium chloride (NaCl), nitrates (NO3−), ortho-phosphates (PO43−), cadmium (Cd), aluminum (Al) and arsenic (As) on 13- to 28-month-old post-parasitic juveniles. The experimental protocol was developed according to conditions described in a standard guide and was modified in order to assess toxicity thresholds for the Dronne River. Results showed that juveniles were tolerant to concentrations of Al, Cd, As, PO43−, NO3− and NaCl, largely higher than those found in the Dronne river, since 96 h EC50s (effective concentrations) were >954 µg/L for Al; >110 µg/L for Cd; >127 µg/L for As; >5.01 mg/L for PO43−; between 1000 and 1500 mg/L for NO3− and were 1.19 and 1.33 g/L for NaCl. Moreover, the use of a substrate in experiments was found not to affect juvenile responses and younger juveniles seemed more sensitive than older individuals. This study thus provides new data about the tolerance of FWPM to environmental pollution and suggests that pollutant concentrations in the river are significantly lower than levels affecting them.

Macroimmunology: The drivers and consequences of spatial patterns in wildlife immune defence

The prevalence and intensity of parasites in wild hosts varies across space and is a key determinant of infection risk in humans, domestic animals and threatened wildlife. Because the immune system serves as the primary barrier to infection, replication and transmission following exposure, we here consider the environmental drivers of immunity. Spatial variation in parasite pressure, abiotic and biotic conditions, and anthropogenic factors can all shape immunity across spatial scales. Identifying the most important spatial drivers of immunity could help pre-empt infectious disease risks, especially in the context of how large-scale factors such as urbanization affect defence by changing environmental conditions. We provide a synthesis of how to apply macroecological approaches to the study of ecoimmunology (i.e. macroimmunology). We first review spatial factors that could generate spatial variation in defence, highlighting the need for large-scale studies that can differentiate competing environmental predictors of immunity and detailing contexts where this approach might be favoured over small-scale experimental studies. We next conduct a systematic review of the literature to assess the frequency of spatial studies and to classify them according to taxa, immune measures, spatial replication and extent, and statistical methods. We review 210 ecoimmunology studies sampling multiple host populations. We show that whereas spatial approaches are relatively common, spatial replication is generally low and unlikely to provide sufficient environmental variation or power to differentiate competing spatial hypotheses. We also highlight statistical biases in macroimmunology, in that few studies characterize and account for spatial dependence statistically, potentially affecting inferences for the relationships between environmental conditions and immune defence. We use these findings to describe tools from geostatistics and spatial modelling that can improve inference about the associations between environmental and immunological variation. In particular, we emphasize exploratory tools that can guide spatial sampling and highlight the need for greater use of mixed-effects models that account for spatial variability while also allowing researchers to account for both individual- and habitat-level covariates. We finally discuss future research priorities for macroimmunology, including focusing on latitudinal gradients, range expansions and urbanization as being especially amenable to large-scale spatial approaches. Methodologically, we highlight critical opportunities posed by assessing spatial variation in host tolerance, using metagenomics to quantify spatial variation in parasite pressure, coupling large-scale field studies with small-scale field experiments and longitudinal approaches, and applying statistical tools from macroecology and meta-analysis to identify generalizable spatial patterns. Such work will facilitate scaling ecoimmunology from individual- to habitat-level insights about the drivers of immune defence and help predict where environmental change may most alter infectious disease risk.

Selection of suitable reference genes for gene expression analysis in gills and liver of fish under field pollution conditions

To understand the role of gene expression in adaptive variation, it is necessary to examine expression variation in an ecological context. Quantitative real-time PCR (qPCR) is considered the most accurate and reliable technique to measure gene expression and to validate the data obtained by RNA-seq; however, accurate normalization is crucial. In Chile, the freshwater silverside fish Basilichthys microlepidotus inhabits both polluted and nonpolluted areas, showing differential gene expression related to pollution. In this study, we infer the stability of six potential reference genes (tubulin alpha, hypoxanthine-guanine phosphoribosyltransferase, glyceraldehyde-3-phosphate dehydrogenase, beta-actin, 60S ribosomal protein L13, and 60S ribosomal protein L8) in the gills and liver of silverside individuals inhabiting polluted and nonpolluted areas. To validate the reference genes selected, the most and least stable reference genes were used to normalize two target transcripts, one for each organ. The RefFinder tool was used to analyze and identify the most stably expressed genes. The 60S ribosomal protein L8 gene was ranked as the most stable gene for both organs. Our results show that reference gene selection influences the detection of differences in the expression levels of target genes in different organs and, also highlighting candidate reference genes that could be used in field studies.

Transcriptomic analysis of zebrafish (Danio rerio) embryos to assess integrated biotoxicity of Xitiaoxi River waters

Assessing the toxicity posed by mixtures of unknown chemicals to aquatic organisms is challenging. In this study, water samples from six cross-sections along the Xitiaoxi River Basin (XRB) were monthly or bimonthly collected in 2014. The year-period physiochemical parameters as well as one-month-water sample based acute biotoxicity tests showed that the river water quality of the year was generally in a good status. High performance liquid chromatography (HPLC) screening based on one-month-water samples suggested that the organic pollutants might be non-to-moderately-polar chemicals in very low concentrations. One-month-water sample based RNA-seq was performed to measure the mRNA differential expression profile of zebrafish larvae to furtherly explore the potential bioeffect and the spatial water quality change of the river. Result indicated that the number of deferentially expressed genes (DEGs) tended to increase along the downstream direction of the river. Gene ontology (GO) enrichment analysis implied that the key pollutants might mainly be the function disruptors of biological processes. Principle components analysis (PCA) combining with transcripts and one-month-water sample based physiochemical parameters indicated that the pollution might be similar at TP, DP and CTB sites while pollution homology existed on some extent between YBQ and JW sites. Although the water quality of the river had a complex time-space alternation during the year, and the one-month-data based RNA-seq could not reflex the whole year-water quality of a watershed, the gene expression profile via RNA-seq provided an alternative way for assessing integrated biotoxicity of surface water, and it was relatively fit for early-warning of water quality of a watershed with unobservable acute toxicity. However, the identification of detail toxicants and the links between DEGs and pollution level as well as physiological-biochemical toxicity needed further investigation.