De novo transcriptome sequencing and analysis of freshwater snail (Radix balthica) to discover genes and pathways affected by exposure to oxazepam

Environmental Risks of Pharmaceutical Mixtures in Aquatic Ecosystems: Reflections on a Decade of Research

Pharmaceuticals and personal care products (PPCPs) occur as variable mixtures in surface waters receiving discharges of human and animal wastes. A key question identified a decade ago is how to assess the effects of long-term exposures of these PPCP mixtures on nontarget organisms. We review the recent progress made on assessing the aquatic ecotoxicity of PPCP mixtures-with a focus on active pharmaceutical ingredients-and the challenges and research needs that remain. New knowledge has arisen from the use of whole-mixture testing combined with component-based approaches, and these studies show that mixtures often result in responses that meet the concentration addition model. However, such studies have mainly been done on individual species over shorter time periods, and longer-term, multispecies assessments remain limited. The recent use of targeted and nontargeted gene analyses has improved our understanding of the diverse pathways that are impacted, and there are promising new "read-across" methods that use mammalian data to predict toxicity in wildlife. Risk assessments remain challenging given the paucity of ecotoxicological and exposure data on PPCP mixtures. As such, the assessment of PPCP mixtures in aquatic environments should remain a priority given the potential for additive-as well as nontarget-effects in nontarget organisms. In addition, we need to improve our understanding of which species, life stages, and relevant endpoints are most sensitive to which types of PPCP mixtures and to expand our knowledge of environmental PPCP levels in regions of the globe that have been poorly studied to date. We recommend an increased use of new approach methodologies, in particular "omics," to advance our understanding of the molecular mechanics of mixture effects. Finally, we call for systematic research on the role of PPCP mixtures in the development of antimicrobial resistance. Environ Toxicol Chem 2024;43:549-558. © 2023 SETAC.

Biochemical and respiratory parameters in a gastropod Radix balthica exposed to diclofenac

Diclofenac, an anti-inflammatory drug, is often detected in natural waters in the ng/L to μg/L range, posing a threat to aquatic organisms. The study focused on the effects of diclofenac in a gastropod mollusk Radix balthica. A 72-h exposure to environmentally relevant concentrations of diclofenac caused deviations from the baseline activities of the studied enzymes in the digestive gland of snails. Acetylcholinesterase activity was induced by the end of exposure, with the most pronounced increase at 3 μg/L. Results on glutathione-S-transferase activity were nonuniform, and no significant variations were observed in thiobarbituric acid reactive substances concentrations, indicating that diclofenac did not cause oxidative stress in the digestive gland of R. balthica at 0.04-4 μg/L range. Diclofenac lowered the oxygen consumption rate in snails in a concentration-dependent manner. At concentrations ≥0.9 μg/L, animals attempted to switch aquatic respiration to breathing air to regulate their metabolic needs. The study showed that diclofenac at environmentally relevant concentrations affected the fitness of R. balthica.

Effects of the anxiolytic benzodiazepine oxazepam on freshwater gastropod reproduction: a prospective study

Psychoactive drugs have emerged as contaminants over the last few decades. These drugs are frequently prescribed and poorly eliminated by wastewater treatment plants, and many are present at non-negligible concentrations in surface waters. Several studies have investigated the non-target organism toxicity of one such drug, oxazepam, a benzodiazepine anxiolytic frequently detected in rivers. However, very little is known about the impact of this drug on reproduction. We investigated the effects of environmentally relevant concentrations of oxazepam on Radix balthica, a freshwater gastropod widespread in Europe. We identified the reproductive organs of Radix balthica. We then exposed this gastropod to oxazepam for two months and assessed several reproductive parameters, from reproductive organ status to behavioral parameters. We found that adults exposed to 10 µg/L oxazepam display an increase in the density of spermatozoa, and that adults exposed to 0.8 µg/L oxazepam displayed a decrease in the number of eggs per egg mass over time. By contrast, oxazepam had no effect on shell length, the size of male reproductive organs or social interactions. Finally, a locomotor activity analysis showed the distance covered over time decreased in all conditions of exposure to oxazepam, potentially reflecting a disturbance of exploratory activity. These results shed light on the effects of oxazepam on the reproduction of a non-target freshwater mollusk.

Low concentrations of oxazepam induce feeding and molecular changes in Radix balthica juveniles

Psychotropics, especially benzodiazepines, are commonly prescribed worldwide. Poorly eliminated at wastewater treatment plants, they belong to a group of emerging contaminants. Due to their interaction with the GABA_A receptor, they may affect the function of the nervous system of non-target organisms, such as aquatic organisms. The toxicity of oxazepam, a very frequently detected benzodiazepine in continental freshwater, has been largely studied in aquatic vertebrates over the last decade. However, its effects on freshwater non-vertebrates have received much less attention. We aimed to evaluate the long-term effects of oxazepam on the juvenile stage of a freshwater gastropod widespread in Europe, Radix balthica. Juveniles were exposed for a month to environmentally-relevant concentrations of oxazepam found in rivers (0.8 μg/L) and effluents (10 μg/L). Three main physiological functions were studied: feeding, growth, and locomotion. Additionally, gene expression analysis was performed to provide insights into toxicity mechanisms. There was a strong short-term activation of the feeding rate at low concentration, whereas the high dose resulted in long-term inhibition of food intake. A significant decrease in mortality rate was observed in juveniles exposed to the lowest dose. Shell growth and locomotor activity did not appear to be affected by oxazepam. Transcriptomic analysis revealed global over-expression of genes involved in the nervous regulation of the feeding, digestive, and locomotion systems after oxazepam exposure. The molecular analysis also revealed a possible interference of animal manipulation with the molecular effects induced by oxazepam exposure. Overall, these results improve our understanding of the effects of the psychoactive drug oxazepam on an aquatic mollusc gastropod.

Gene expression profiles at different stages for formation of pearl sac and pearl in the pearl oyster Pinctada fucata

Background

The most critical step in the pearl formation during aquaculture is issued to the proliferation and differentiation of outer epithelial cells of mantle graft into pearl sac. This pearl sac secretes various matrix proteins to produce pearls by a complex physiological process which has not been well-understood yet. Here, we aimed to unravel the genes involved in the development of pearl sac and pearl, and the sequential expression patterns of different shell matrix proteins secreted from the pearl sac during pearl formation by pearl oyster Pinctada fucata using high-throughput transcriptome profiling.

Results

Principal component analysis (PCA) showed clearly different gene expression profiles between earlier (before 1 week) and later stages (1 week to 3 months) of grafting. Immune-related genes were highly expressed between 0 h – 24 h (donor dependent) and 48 h – 1 w (host dependent), and in the course of wound healing process pearl sac was developed by two weeks of graft transplantation. Moreover, for the first time, we identified some stem cell marker genes including ABCG2, SOX2, MEF2A, HES1, MET, NRP1, ESR1, STAT6, PAX2, FZD1 and PROM1 that were expressed differentially during the formation of pearl sac. The expression profiling of 192 biomineralization-related genes demonstrated that most of the shell matrix proteins (SMPs) involved in prismatic layer formation were first up-regulated and then gradually down-regulated indicating their involvement in the development of pearl sac and the onset of pearl mineralization. Most of the nacreous layer forming SMPs were up-regulated at 2 weeks after the maturation of pearl sac. Nacrein, MSI7 and shematrin involved in both layer formation were highly expressed during 0 h – 24 h, down-regulated up to 1 week and then up-regulated again after accomplishment of pearl sac formation.

Conclusions

Using an RNA-seq approach we unraveled the expression pattern of the key genes involved in the development of pearl sac and pearl as a result of host immune response after grafting. These findings provide valuable information in understanding the molecular mechanism of pearl formation and immune response in P. fucata.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5579-3) contains supplementary material, which is available to authorized users.

Combining the assessment of apical endpoints and gene expression in the freshwater snail Physa acuta after exposure to reclaimed water

Post-treatment wastewater reuses are diverse. Recreational and environmental restoration uses of reclaimed water (RW) can be potentially harmful to aquatic organisms. In this work the freshwater snail Physa acuta was exposed to RW (100%) and its dilution (RW 50%). A simple laboratory mixture of three emerging pollutants was used to address the complex problem of mixture toxicity of RW. Hence fortified reclaimed water (FRW), obtained by adding fluoxetine (400 μg FLX/L), perfluorooctane sulphonic acid (90 μg PFOS/L) and methylparaben (9 μg MP/L), was tested at two dilution percentages: 100% and 50%. The effects of the laboratory mixture of FLX, PFOS and MP on the test medium were also studied. Long-lasting effects, together with early molecular responses, were assessed. Fecundity (cumulative egg production) over 21 days and the hatching of produced eggs (F1) after another 21-day embryonic exposure were monitored. The gene expression of three genes was analysed after 24 h of exposure: two endocrine-related nuclear receptors (ERR and RXR) and one stress protein gene (Hsp70). This reproduction test, with additional assessments of the F1 recovered eggs' hatching success, showed that both RW and FRW significantly reduced fecundity. F1 hatching was affected only by FRW. The gene expression results showed that the RXR response was strikingly similar to the fecundity response, which suggests that this nuclear receptor is involved in the reproductive pathways of gastropods. ERR remained virtually unaltered. Hsp70 was overexpressed by the laboratory mixture in the test medium, but no effect was observed in the fortification of RW. This opposite effect and lack of response for F1 hatching produced by the laboratory mixture in the test medium highlighted the difficulty of predicting mixture effects. The experimental approach allowed us to test the effects caused by RW on P. acuta at different biological organisation levels. Thus, the combination of molecular biomarkers and ecological relevant endpoints is a good strategy to test complex mixtures like RW as it provides a framework to link mechanisms of action and whole organism effects when it is almost impossible to detect the pollutant(s) that cause toxic effects.

Evaluation of psychiatric hospital wastewater toxicity: what is its impact on aquatic organisms?

The primary source of pharmaceuticals to the aquatic environment is the discharge of wastewater effluents. Pharmaceuticals are a large and diverse group of compounds. Among them, psychotropic substances are particularly interesting to study due to their specific known mode of action. The present study was performed to investigate the effects of wastewater effluents from a psychiatric hospital wastewater treatment plant (WWTP) on several aquatic organisms. All the analyzed pharmaceuticals (10 compounds) were detected in WWTP effluents as well as in the receiving river. Although the environmental concentrations were generally at trace levels (ng L^-1 to μg L^-1), induce toxic effects were observed. This study showed the effects of the WWTP effluents on the oogenesis and/or embryogenesis of amphipod crustacean Gammarus fossarum, Japanese fish medaka Oryzias latipes, mollusk Radix peregra, and planarian Schmidtea polychroa. A decrease of the number of oocytes and produced embryos was observed for G. fossarum and S. polychroa. Similarly, the hatching rate of R. peregra was affected by effluents. In the receiving river, the macroinvertebrate community was affected by the wastewater effluents discharge.

The root transcriptome of Achyranthes bidentata and the identification of the genes involved in the replanting benefit

The transcriptome profiling in replanting roots revealed that expression pattern changes of key genes promoted important metabolism pathways, antioxidant and pathogen defense systems, adjusted phytohormone signaling and inhibited lignin biosynthesis. The yield of the medicinal plant Achyranthes bidentata could be significantly increased when replanted into a field cultivated previously for the same crop, but the biological basis of this so-called "replanting benefit" is unknown. Here, the RNA-seq technique was used to identify candidate genes responsible for the benefit. The analysis of RNA-seq libraries prepared from mRNA extracted from the roots of first year planting (normal growth, NG) and second year replanting (consecutive monoculture, CM) yielded about 40.22 GB sequencing data. After de novo assembly, 87,256 unigenes were generated with an average length of 1060 bp. Among these unigenes, 55,604 were annotated with public databases, and 52,346 encoding sequences and 2881 transcription factors were identified. A contrast between the NG and CM libraries resulted in a set of 3899 differentially transcribed genes (DTGs). The DTGs related to the replanting benefit and their expression profiles were further analyzed by bioinformatics and qRT-PCR approaches. The major differences between the NG and CM transcriptomes included genes encoding products involved in glycolysis/gluconeogenesis, glutathione metabolism and antioxidant defense, in aspects of the plant/pathogen interaction, phytohormone signaling and phenylpropanoid biosynthesis. The indication was that replanting material enjoyed a stronger level of defense systems, a balance regulation of hormone signals and a suppression of lignin formation, thereby promoting root growth and development. The study provides considerable significant insights for a better understanding of the molecular mechanism of the replanting benefit and suggests their possible application in developing methods to reinforce the effects in medicinal plants.

Deep transcriptome profiling sheds light on key players in nucleus implantation induced immune response in the pearl oyster Pinctada martensii

Immunological rejection of the pearl oysters following nucleus implantation is a major issue limiting the successful rate of cultured pearls. To date, the molecular mechanism of immune tolerance during pearl formation in the pearl oysters is still largely unknown. Through the RNA sequencing platform and comparative transcriptomic analysis, we investigated the chronic gene expression changes at seven time points (0, 5, 10, 15, 20, 30, 60 days post implantation or dpi) over a period of 60 days following nucleus implantation in the pearl oyster Pinctada martensii. A total of 81,390 unique transcripts (or unigenes) with a combined length of 96.8 million bp and a N50 value of 2227 bp were obtained. When compared with sequences in the nr, nt, Swiss-Prot, KEGG, COG and GO databases, 36,380 unigenes can find homologous genes. Pairwise comparison of gene expression among all the samples showed that the largest number (or 6846) of differentially expressed genes was observed at 10 dpi. The number then decreased to below 5000 at 15, 20 and 30 dpi and increased again to 6679 at 60 dpi. PCA analysis further showed that the seven time points can be roughly divided into four groups. Comparative transcriptomic analysis between the four groups identified a variety of genes showing differential expression at different time points, including many immune-related genes such as those encoding for toll-like receptor, lectin, scavenger receptor, and peroxidase. In addition, GO and KEGG enrichment analysis revealed that these differentially expressed genes were mainly associated with metabolism, ribosome function, immune response, signaling transduction, and cytoskeleton organization. Notably, two KEGG pathways, namely "cell adhesion molecules" and "primary immunodeficiency" were significantly enriched during the whole process. This finding indicates that genes in these pathways are likely to play critical roles in the immune tolerance of the pearl oysters. To conclude, the data obtained contribute to a better understanding of the molecular mechanisms of nucleus implantation induced immune response in the pearl oysters, and will facilitate the development of effective measures to improve the performance of pearl culture.