Characterisation of the transcriptome of male and female wild-type guppy brains with RNA-Seq and consequences of exposure to the pharmaceutical pollutant, 17α-ethinyl estradiol

Hepatic transcriptome, transcriptional effects and antioxidant responses in Poecilia vivipara exposed to sanitary sewage

Aquatic environments are subject to threats from multiple human activities, particularly through the release of untreated sanitary sewage into the coastal environments. These effluents contain a large group of natural or synthetic compounds referred to as emerging contaminants. Monitoring the types and quantities of toxic substances in the environment, especially complex mixtures, is an exhausting and challenging task. Integrative effect-based tools, such as biomarkers, are recommended for environmental quality monitoring programs. In this study, fish Poecilia vivipara were exposed for 24 and 96 h to raw untreated sewage diluted 33 % (v/v) in order to identify hepatic genes to be used as molecular biomarkers. Through a de novo hepatic transcriptome assembly, using Illumina MiSeq, 54,285 sequences were assembled creating a reference transcriptome for this guppy species. Transcripts involved in biotransformation systems, antioxidant defenses, ABC transporters, nuclear and xenobiotic receptors were identified and evaluated by qPCR. Sanitary sewage induced transcriptional changes in AhR, PXR, CYP2K1, CYP3A30, NQO1, UGT1A1, GSTa3, GSTmu, ST1C1, SOD, ABCC1 and SOX9 genes from liver of fish, particularly after 96 h of exposure. Changes in hepatic enzyme activities were also observed. The enzymes showed differences in fish exposed to both periods, while in the gills there was a prevalence of significant results after 96 h. The observed differences were associated to gender and/or to sewage exposure. The obtained results support the use of P. vivipara as sentinel and model organism for ecotoxicological studies and evidence the importance of understanding the differential responses associated to gender.

Sub-acute exposure of male guppies (Poecilia reticulata) to environmentally relevant concentrations of PFOA and GenX induces significant changes in the testis transcriptome and reproductive traits

Poly- and perfluoroalkyl substances (PFAS) are frequently detected in the environment and are linked to adverse reproductive health outcomes in humans. Although legacy PFAS have been phased out due to their toxicity, alternative PFAS are increasingly used despite the fact that information on their toxic effects on reproductive traits is particularly scarce. Here, we exposed male guppies (Poecilia reticulata) for a short period (21 days) to an environmentally realistic concentration (1 ppb) of PFOA, a legacy PFAS, and its replacement compound, GenX, to assess their impact on reproductive traits and gene expression. Exposure to PFAS did not impair survival but instead caused sublethal effects. Overall, PFAS exposure caused changes in male sexual behaviour and had detrimental effects on sperm motility. Sublethal variations were also seen at the transcriptional level, with the modulation of genes involved in immune regulation, spermatogenesis, and oxidative stress. We also observed bioaccumulation of PFAS, which was higher for PFOA than for GenX. Our results offer a comprehensive comparison of these two PFAS and shed light on the toxicity of a newly emerging alternative to legacy PFAS. It is therefore evident that even at low concentrations and with short exposure, PFAS can have subtle yet significant effects on behaviour, fertility, and immunity. These findings underscore the potential ramifications of pollution under natural conditions and their impact on fish populations.

Bioinformatics approaches and big data analytics opportunities in improving fisheries and aquaculture

Aquaculture has witnessed an excellent growth rate during the last two decades and offers huge potential to provide nutritional as well as livelihood security. Genomic research has contributed significantly toward the development of beneficial technologies for aquaculture. The existing high throughput technologies like next-generation technologies generate oceanic data which requires extensive analysis using appropriate tools. Bioinformatics is a rapidly evolving science that involves integrating gene based information and computational technology to produce new knowledge for the benefit of aquaculture. Bioinformatics provides new opportunities as well as challenges for information and data processing in new generation aquaculture. Rapid technical advancements have opened up a world of possibilities for using current genomics to improve aquaculture performance. Understanding the genes that govern economically relevant characteristics, necessitates a significant amount of additional research. The various dimensions of data sources includes next-generation DNA sequencing, protein sequencing, RNA sequencing gene expression profiles, metabolic pathways, molecular markers, and so on. Appropriate bioinformatics tools are developed to mine the biologically relevant and commercially useful results. The purpose of this scoping review is to present various arms of diverse bioinformatics tools with special emphasis on practical translation to the aquaculture industry.

Gene expression and biochemical patterns in the digestive gland of the mussel Mytilus galloprovincialis (Lamarck, 1819) exposed to 17α-ethinylestradiol

Contaminants of emerging concern (CECs) are a class of chemicals that can spread throughout the environment and may cause adverse biological and ecological effects. While there are many different classes of CECs, one of the most well documented in the aquatic environment are pharmaceutical drugs, such as natural and synthetic estrogens. In particular, the widespread presence of the synthetic estrogen 17 α-Ethinylestradiol (EE2) in water may lead to bioaccumulation in sediment and biota. EE2 is the primary component in contraceptive pills, and is a derivative of the natural hormone estradiol (E2). In this study, the mussel Mytilus galloprovincialis was exposed to EE2 in a semi-static and time-dependent experiment, for a total exposure period of 28 days. Biochemical and transcriptomics analyses were performed on mussel digestive glands after exposure for 14 (T14) and 28 (T28) days. Metabolic and DNA impairments, as well as activation of antioxidant and biotransformation enzymes activation, were detected in T28 exposed mussels. RNA-Seq analysis showed significant differential expression of 160 (T14 compared to controls), 33 (T28 compared to controls) and 79 (T14 compared to T28) genes. Signs of stress after EE2 treatment included up-regulation of gene/proteins involved with immune function, lipid transport, and metabolic and antibacterial properties. This study elucidates the underlying mechanisms of EE2 in a filter feeding organisms to elucidate the effects of this human pharmaceutical on aquatic biota.

Integrated Genomic and Bioinformatics Approaches to Identify Molecular Links between Endocrine Disruptors and Adverse Outcomes

Exposure to Endocrine Disrupting Chemicals (EDC) has been linked with several adverse outcomes. In this review, we examine EDCs that are pervasive in the environment and are of concern in the context of human, animal, and environmental health. We explore the consequences of EDC exposure on aquatic life, terrestrial animals, and humans. We focus on the exploitation of genomics technologies and in particular whole transcriptome sequencing. Genome-wide analyses using RNAseq provides snap shots of cellular, tissue and whole organism transcriptomes under normal physiological and EDC perturbed conditions. A global view of gene expression provides highly valuable information as it uncovers gene families or more specifically, pathways that are affected by EDC exposures, but also reveals those that are unaffected. Hypotheses about genes with unknown functions can also be formed by comparison of their expression levels with genes of known function. Risk assessment strategies leveraging genomic technologies and the development of toxicology databases are explored. Finally, we review how the Adverse Outcome Pathway (AOP) has exploited this high throughput data to provide a framework for toxicology studies.

Clustering of Sex-Biased Genes and Transposable Elements in the Genome of the Medaka Fish Oryzias latipes

Although genes with similar expression patterns are sometimes found in the same genomic regions, almost nothing is known about the relative organization in genomes of genes and transposable elements (TEs), which might influence each other at the regulatory level. In this study, we used transcriptomic data from male and female gonads of the Japanese medaka Oryzias latipes to define sexually biased genes and TEs and analyze their relative genomic localization. We identified 20,588 genes expressed in the adult gonads of O. latipes. Around 39% of these genes are differentially expressed between male and female gonads. We further analyzed the expression of TEs using the program SQuIRE and showed that more TE copies are overexpressed in testis than in ovaries (36% vs. 10%, respectively). We then developed a method to detect genomic regions enriched in testis- or ovary-biased genes. This revealed that sex-biased genes and TEs are not randomly distributed in the genome and a part of them form clusters with the same expression bias. We also found a correlation of expression between TE copies and their closest genes, which increases with decreasing intervening distance. Such a genomic organization suggests either that TEs hijack the regulatory sequences of neighboring sexual genes, allowing their expression in germ line cells and consequently new insertions to be transmitted to the next generation, or that TEs are involved in the regulation of sexual genes, and might therefore through their mobility participate in the rewiring of sex regulatory networks.

Transcriptome-wide changes associated with the reproductive behaviour of male guppies exposed to 17α-ethinyl estradiol

Although many pharmaceutical compounds (and their metabolites) can induce harmful impacts at the molecular, physiological and behavioural levels, their underlying mechanistic associations have remained largely unexplored. Here, we utilized RNA-Seq to build a whole brain transcriptome profile to examine the impact of a common endocrine disrupting pharmaceutical (17α-ethinyl estradiol, EE2) on reproductive behaviour in wild guppies (Poecilia reticulata). Specifically, we annotated 16,791 coding transcripts in whole brain tissue in relation to the courtship behaviour (i.e. sigmoid display) of EE2 exposed (at environmentally relevant concentration of 8 ng/L for 28-days) and unexposed guppies. Further, we obtained 10,960 assembled transcripts matching in the non-coding orthologous genomes. Behavioural responses were assessed using a standard mate choice experiment, which allowed us to disentangle chemical cues from visual cues. We found that a high proportion of the RNAseq reads aligned back to our de novo assembled transcriptome with 80.59% mapping rate. Behavioural experiments showed that when males were presented only with female visual cues, there was a significant interaction between male treatment and female treatment in the time spent in the preference zone. This is one of the first studies to show that transcriptome-wide changes are associated with the reproductive behaviour of fish: EE2 exposed male guppies that performed high levels of courtship had a gene profile that deviated the most from the other treatment groups, while both non-courting EE2 and control males had similar gene signatures. Using Gene Ontology pathway analysis, our study shows that EE2-exposed males had gene transcripts enriched for pathways associated with altered immunity, starvation, altered metabolism and spermatogenesis. Our study demonstrates that multiple gene networks orchestrate courting behaviour, emphasizing the importance of investigating impacts of pharmaceuticals on gene networks instead of single genes.

De Novo Transcriptome Assembly and Annotation of Liver and Brain Tissues of Common Brushtail Possums (Trichosurus vulpecula) in New Zealand: Transcriptome Diversity after Decades of Population Control

The common brushtail possum (Trichosurus vulpecula), introduced from Australia in the mid-nineteenth century, is an invasive species in New Zealand where it is widespread and forms the largest self-sustained reservoir of bovine tuberculosis (Mycobacterium bovis) among wild populations. Conservation and agricultural authorities regularly apply a series of population control measures to suppress brushtail possum populations. The evolutionary consequence of more than half a century of intensive population control operations on the species’ genomic diversity and population structure is hindered by a paucity of available genomic resources. This study is the first to characterise the functional content and diversity of brushtail possum liver and brain cerebral cortex transcriptomes. Raw sequences from hepatic cells and cerebral cortex were assembled into 58,001 and 64,735 transcripts respectively. Functional annotation and polymorphism assignment of the assembled transcripts demonstrated a considerable level of variation in the core metabolic pathways that represent potential targets for selection pressure exerted by chemical toxicants. This study suggests that the brushtail possum population in New Zealand harbours considerable variation in metabolic pathways that could potentially promote the development of tolerance against chemical toxicants.

Twenty years of transcriptomics, 17alpha-ethinylestradiol, and fish

In aquatic toxicology, perhaps no pharmaceutical has been investigated more intensely than 17alpha-ethinylestradiol (EE2), the active ingredient of the birth control pill. At the turn of the century, the fields of comparative endocrinology and endocrine disruption research witnessed the emergence of omics technologies, which were rapidly adapted to characterize potential hazards associated with exposures to environmental estrogens, such as EE2. Since then, significant advances have been made by the scientific community, and as a result, much has been learned about estrogen receptor signaling in fish from environmental xenoestrogens. Vitellogenin, the egg yolk precursor protein, was identified as a major estrogen-responsive gene, establishing itself as the premier biomarker for estrogenic exposures. Omics studies have identified a plethora of estrogen responsive genes, contributing to a wealth of knowledge on estrogen-mediated regulatory networks in teleosts. There have been ~40 studies that report on transcriptome responses to EE2 in a variety of fish species (e.g., zebrafish, fathead minnows, rainbow trout, pipefish, mummichog, stickleback, cod, and others). Data on the liver and testis transcriptomes dominate in the literature and have been the subject of many EE2 studies, yet there remain knowledge gaps for other tissues, such as the spleen, kidney, and pituitary. Inter-laboratory genomics studies have revealed transcriptional networks altered by EE2 treatment in the liver; networks related to amino acid activation and protein folding are increased by EE2 while those related to xenobiotic metabolism, immune system, circulation, and triglyceride storage are suppressed. EE2-responsive networks in other tissues are not as comprehensively defined which is a knowledge gap as regulated networks are expected to be tissue-specific. On the horizon, omics studies for estrogen-mediated effects in fish include: (1) Establishing conceptual frameworks for incorporating estrogen-responsive networks into environmental monitoring programs; (2) Leveraging in vitro and computational toxicology approaches to identify chemicals associated with estrogen receptor-mediated effects in fish (e.g., male vitellogenin production); (3) Discovering new tissue-specific estrogen receptor signaling pathways in fish; and (4) Developing quantitative adverse outcome pathway predictive models for estrogen signaling. As we look ahead, research into EE2 over the past several decades can serve as a template for the array of hormones and endocrine active substances yet to be fully characterized or discovered.

Environmentally relevant concentrations of bifenthrin affect the expression of estrogen and glucocorticoid receptors in brains of female western mosquitofish

In recent decades, pyrethroid pesticides have been deemed a safer alternative to previously used pesticides. While some evidence supports this assumption in mammals and birds, exposure to certain pyrethroids can affect concentrations of hormones vital to reproduction in fish. Thus, we hypothesized that pyrethroid exposure impacts fish reproductive behavior and the expression of genes associated with reproduction. We tested our hypothesis by examining effects of the widely used pyrethroid pesticide, bifenthrin, on the reproductive behaviors of the broadly distributed livebearing western mosquitofish, Gambusia affinis. We exposed sexually mature female fish to one of five environmentally relevant concentrations of bifenthrin and conducted behavioral assays to assess reproductive, social, and space use behaviors before and after exposure. We did not detect changes in behaviors measured in response to bifenthrin. However, exposure was associated with increased expression of an estrogen receptor gene (ER-α) and glucocorticoid receptor (GR) in brain tissue at bifenthrin concentrations at concentrations of 5.90 and 24.82 ng/L, and 5.90 and 12.21 ng/L, respectively. Our study supports the perspective that the use of multiple endpoints through integrative approaches is essential for understanding the cumulative impact of pollutants. Integrating physiological, morphological, and behavioral investigations of nonlethal concentrations of pollutants like bifenthrin may heighten our potential to predict their impact on individuals, populations, and communities.

Early-life exposure to the endocrine disruptor 17-α-ethinylestradiol induces delayed effects in adult brain, liver and ovotestis proteomes of a self-fertilizing fish

Early-life represents a critically sensitive window to endocrine disrupting chemicals, potentially leading to long-term repercussions on the phenotype later in life. The mechanisms underlying this phenomenon, referred to as the Developmental Origins of Health and Disease (DOHaD), are still poorly understood. To gain molecular understanding of these effects, we exposed mangrove rivulus (Kryptolebias marmoratus) for 28 days post hatching (dph) to 4 and 120 ng/L 17-α-ethinylestradiol, a model xenoestrogen. After 28 days, fish were raised for 140 days in clean water and we performed quantitative label-free proteomics on brain, liver and ovotestis of 168 dph adults. A total of 820, 888 and 420 proteins were robustly identified in the brain, liver and ovotestis, respectively. Effects of 17-α-ethinylestradiol were tissue and dose-dependent: a total of 31, 51 and 18 proteins were differentially abundant at 4 ng/L in the brain, liver and ovotestis, respectively, compared to 20, 25 and 39 proteins at 120 ng/L. Our results suggest that estrogen-responsive pathways, such as lipid metabolism, inflammation, and the innate immune system were affected months after the exposure. In addition, the potential perturbation of S-adenosylmethionine metabolism encourages future studies to investigate the role of DNA methylation in mediating the long-term effects of early-life exposures. SIGNIFICANCE: The Developmental Origins of Health and Disease (DOHaD) states that early life stages of humans and animals are sensitive to environmental stressors and can develop health issues later in life, even if the stress has ceased. Molecular mechanisms supporting DOHaD are still unclear. The mangrove rivulus is a new fish model species naturally reproducing by self-fertilization, making it possible to use isogenic lineages in which all individuals are highly homozygous. This species therefore permits to strongly reduce the confounding factor of genetic variability in order to investigate the effects of environmental stress on the phenotype. After characterizing the molecular phenotype of brain, liver and ovotestis, we obtained true proteomic reaction norms of these three organs in adults after early life stages have been exposed to the common endocrine disruptor 17-α-ethinylestradiol (EE2). Our study demonstrates long-term effects of early-life endocrine disruption at the proteomic level in diverse estrogen-responsive pathways 5 months after the exposure. The lowest tested and environmentally relevant concentration of 4 ng/L had the highest impact on the proteome in brain and liver, highlighting the potency of endocrine disruptors at low concentrations.

Testis transcriptome alterations in zebrafish (Danio rerio) with reduced fertility due to developmental exposure to 17α-ethinyl estradiol

17α-Ethinylestradiol (EE₂) is a ubiquitous aquatic contaminant shown to decrease fish fertility at low concentrations, especially in fish exposed during development. The mechanisms of the decreased fertility are not fully understood. In this study, we perform transcriptome analysis by RNA sequencing of testes from zebrafish with previously reported lowered fertility due to exposure to low concentrations of EE₂ during development. Fish were exposed to 1.2 and 1.6 ng/L (measured concentration; nominal concentrations 3 and 10 ng/L) of EE₂ from fertilization to 80 days of age, followed by 82 days of remediation in clean water. RNA sequencing analysis revealed 249 and 16 genes to be differentially expressed after exposure to 1.2 and 1.6 ng/L, respectively; a larger inter-sample variation was noted in the latter. Expression of 11 genes were altered by both exposures and in the same direction. The coding sequences most affected could be categorized to the putative functions cell signalling, proteolysis, protein metabolic transport and lipid metabolic process. Several homeobox transcription factors involved in development and differentiation showed increased expression in response to EE₂ and differential expression of genes related to cell death, differentiation and proliferation was observed. In addition, several genes related to steroid synthesis, testis development and function were differentially expressed. A number of genes associated with spermatogenesis in zebrafish and/or mouse were also found to be differentially expressed. Further, differences in non-coding sequences were observed, among them several differentially expressed miRNA that might contribute to testis gene regulation at post-transcriptional level. This study has generated insights of changes in gene expression that accompany fertility alterations in zebrafish males that persist after developmental exposure to environmental relevant concentrations of EE₂ that persist followed by clean water to adulthood. Hopefully, this will generate hypotheses to test in search for mechanistic explanations.