Transcriptomic and Network Analyses Reveal Mechanistic-Based Biomarkers of Endocrine Disruption in the Marine Mussel, Mytilus edulis.

Transcriptomics, high-throughput assays, and adverse outcome pathways (AOP) are promising approaches applied to toxicity monitoring in the 21st century, but development of these methods is challenging for nonmodel organisms and emerging contaminants. For example, Endocrine Disrupting Compounds (EDCs) may cause reproductive impairments and feminization of male bivalves; however, the mechanism linked to this adverse outcome is unknown. To develop mechanism-based biomarkers that may be linked through an AOP, we exposed Mytilus edulis to 17-alpha-ethinylestradiol (5 and 50 ng/L) and 4-nonylphenol (1 and 100 µg/L) for 32 and 39 days. When mussels were exposed to these EDCs, we found elevated female specific transcripts and significant female-skewed sex ratios using a RT-qPCR assay. We performed gene expression analysis on digestive gland tissue using an M. edulis microarray and through network and targeted analyses identified the nongenomic estrogen signaling pathway and steroidogenesis pathway as the likely mechanisms of action for a putative AOP. We also identified several homologues to genes within the vertebrate steroidogenesis pathway including the cholesterol side chain cleavage complex. From this AOP, we designed the Coastal Biosensor for Endocrine Disruption (C-BED) assay which was confirmed in the laboratory and tested in the field.

The Toxicogenome of Hyalella azteca: A Model for Sediment Ecotoxicology and Evolutionary Toxicology.

Hyalella azteca is a cryptic species complex of epibenthic amphipods of interest to ecotoxicology and evolutionary biology. It is the primary crustacean used in North America for sediment toxicity testing and an emerging model for molecular ecotoxicology. To provide molecular resources for sediment quality assessments and evolutionary studies, we sequenced, assembled, and annotated the genome of the H. azteca U.S. Lab Strain. The genome quality and completeness is comparable with other ecotoxicological model species. Through targeted investigation and use of gene expression data sets of H. azteca exposed to pesticides, metals, and other emerging contaminants, we annotated and characterized the major gene families involved in sequestration, detoxification, oxidative stress, and toxicant response. Our results revealed gene loss related to light sensing, but a large expansion in chemoreceptors, likely underlying sensory shifts necessary in their low light habitats. Gene family expansions were also noted for cytochrome P450 genes, cuticle proteins, ion transporters, and include recent gene duplications in the metal sequestration protein, metallothionein. Mapping of differentially expressed transcripts to the genome significantly increased the ability to functionally annotate toxicant responsive genes. The H. azteca genome will greatly facilitate development of genomic tools for environmental assessments and promote an understanding of how evolution shapes toxicological pathways with implications for environmental and human health.

Multiple origins of pyrethroid insecticide resistance across the species complex of a nontarget aquatic crustacean, Hyalella azteca.

Use of pesticides can have substantial nonlethal impacts on nontarget species, including driving evolutionary change, often with unknown consequences for species, ecosystems, and society. Hyalella azteca, a species complex of North American freshwater amphipods, is widely used for toxicity testing of water and sediment and has frequently shown toxicity due to pyrethroid pesticides. We demonstrate that 10 populations, 3 from laboratory cultures and 7 from California water bodies, differed by at least 550-fold in sensitivity to pyrethroids. The populations sorted into four phylogenetic groups consistent with species-level divergence. By sequencing the primary pyrethroid target site, the voltage-gated sodium channel, we show that point mutations and their spread in natural populations were responsible for differences in pyrethroid sensitivity. At least one population had both mutant and WT alleles, suggesting ongoing evolution of resistance. Although nonresistant H. azteca were susceptible to the typical neurotoxic effects of pyrethroids, gene expression analysis suggests the mode of action in resistant H. azteca was not neurotoxicity but was oxidative stress sustained only at considerably higher pyrethroid concentrations. The finding that a nontarget aquatic species has acquired resistance to pesticides used only on terrestrial pests is troubling evidence of the impact of chronic pesticide transport from land-based applications into aquatic systems. Our findings have far-reaching implications for continued uncritical use of H. azteca as a principal species for monitoring and environmental policy decisions.

Toxicogenomic responses of nanotoxicity in Daphnia magna exposed to silver nitrate and coated silver nanoparticles.

Applications for silver nanomaterials in consumer products are rapidly expanding, creating an urgent need for toxicological examination of the exposure potential and ecological effects of silver nanoparticles (AgNPs). The integration of genomic techniques into environmental toxicology has presented new avenues to develop exposure biomarkers and investigate the mode of toxicity of novel chemicals. In the present study we used a 15k oligonucleotide microarray for Daphnia magna, a freshwater crustacean and common indicator species for toxicity, to differentiate between particle specific and ionic silver toxicity and to develop exposure biomarkers for citrate-coated and PVP-coated AgNPs. Gene expression profiles revealed that AgNO(3) and AgNPs have distinct expression profiles suggesting different modes of toxicity. Major biological processes disrupted by the AgNPs include protein metabolism and signal transduction. In contrast, AgNO(3) caused a downregulation of developmental processes, particularly in sensory development. Metal responsive and DNA damage repair genes were induced by the PVP AgNPs, but not the other treatments. In addition, two specific biomarkers were developed for the environmental detection of PVP AgNPs; although further verification under different environmental conditions is needed.

Assessing legacy and endocrine disrupting pollutants in Boston Harbor with transcriptomic biomarkers.

Within monitoring frameworks, biomarkers provide several benefits because they serve as intermediates between pollutant exposure and effects, and integrate the responses of contaminants that operate through the same mechanism of action. This study was designed to verify the use of transcriptomic biomarkers developed in our prior work (i.e., Coastal Biosensor of Endocrine Disruption; C-BED assay) on Mytilus edulis and identify additional biomarkers for legacy pollutants. M. edulis were collected from a reference site in Pemaquid, ME, USA and deployed by the Massachusetts Water Resources Authority (MWRA) at locations in and outside Boston Harbor, MA, USA: including (1) Boston Inner Harbor (IH), (2) the current outfall (OS), (3) 1 km away from the current outfall (LNB), and (4) Deer Island (DI), the site where untreated wastewater was formerly discharged into the bay. Differential gene expression was quantified with a high density microarray. Seven genes significantly correlated with whole tissue concentration of PAHs, and six genes significantly correlated with whole body concentrations of PCBs, two groups of legacy contaminants that were elevated at stations IH, OS, and DI. Enrichment analysis indicated that IH mussels had the highest induction of stress response genes, which correlated with the higher levels of contaminants measured at this site. Based on the C-BED assay gene analysis, stations IH and OS exhibited signs of endocrine disruption, which were further confirmed by incorporating the results for the C-BED assay within the Integrated Biomarker Response (IBR) approach. This study successfully demonstrated the potential use of transcriptomic biomarkers within a monitoring program to identify the presence and organismal responses to endocrine disrupting and legacy contaminant classes.

Toward Sustainable Environmental Quality: Priority Research Questions for North America.

Anticipating, identifying, and prioritizing strategic needs represent essential activities by research organizations. Decided benefits emerge when these pursuits engage globally important environment and health goals, including the United Nations Sustainable Development Goals. To this end, horizon scanning efforts can facilitate identification of specific research needs to address grand challenges. We report and discuss 40 priority research questions following engagement of scientists and engineers in North America. These timely questions identify the importance of stimulating innovation and developing new methods, tools, and concepts in environmental chemistry and toxicology to improve assessment and management of chemical contaminants and other diverse environmental stressors. Grand challenges to achieving sustainable management of the environment are becoming increasingly complex and structured by global megatrends, which collectively challenge existing sustainable environmental quality efforts. Transdisciplinary, systems-based approaches will be required to define and avoid adverse biological effects across temporal and spatial gradients. Similarly, coordinated research activities among organizations within and among countries are necessary to address the priority research needs reported here. Acquiring answers to these 40 research questions will not be trivial, but doing so promises to advance sustainable environmental quality in the 21st century. Environ Toxicol Chem 2019;38:1606-1624. © 2019 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals, Inc. on behalf of SETAC.

Bioconcentration and depuration of (14)C-labeled 17α-ethinyl estradiol and 4-nonylphenol in individual organs of the marine bivalve Mytilus edulis L. .

Endocrine-disrupting compounds (EDCs), including 17α-ethinyl estradiol (EE2) and 4-nonylphenol (4-NP), enter coastal environments primarily in effluents of wastewater treatment facilities and have become ubiquitous in marine surface waters, sediments, and biota. Although EE2 and 4-NP have been detected in marine shellfish, the kinetics of bioconcentration and their tissue distribution have not been thoroughly investigated. The authors performed bioconcentration and depuration experiments in the blue mussel, Mytilus edulis, with 3.37 nM EE2 (0.999 µg/L) and 454 nM 4-NP (100.138 µg/L). Mussels and seawater were sampled throughout a 38-d exposure and a 35-d depuration period, and 6 tissues were individually assayed. Uptake of EE2 and 4-NP was curvilinear throughout exposure and followed a similar uptake pattern: digestive gland > gill ≥ remaining viscera > gonad > adductor > plasma. Depuration varied, however, with half-lives ranging from 2.7 d (plasma) to 92 d (gill) for EE2 and 15 d (plasma) to 57 d (gill) for 4-NP. An innovative modeling approach, with 3 coupled mathematical models, was developed to differentiate the unique roles of the gill and plasma in distributing the EDCs to internal tissues. Plasma appears pivotal in regulating EDC uptake and depuration within the whole mussel.

Correlation of transcriptomic responses and metal bioaccumulation in Mytilus edulis L. reveals early indicators of stress.

Marine biomonitoring programs in the U.S. and Europe have historically relied on monitoring tissue concentrations of bivalves to monitor contaminant levels and ecosystem health. By integrating 'omic methods with these tissue residue approaches we can uncover mechanistic insight to link tissue concentrations to potential toxic effects. In an effort to identify novel biomarkers and better understand the molecular toxicology of metal bioaccumulation in bivalves, we exposed the blue mussel, Mytilus edulis L., to sub-lethal concentrations (0.54 µM) of cadmium, lead, and a Cd+Pb mixture. Metal concentrations were measured in gill tissues at 1, 2, and 4 weeks, and increased linearly over the 4 week duration. In addition, there was evidence that Pb interfered with Cd uptake in the mixture treatment. Using a 3025 sequence microarray for M. edulis, we performed transcriptomic analysis, identifying 57 differentially expressed sequences. Hierarchical clustering of these sequences successfully distinguished the different treatment groups demonstrating that the expression profiles were reproducible among the treatments. Enrichment analysis of gene ontology terms identified several biological processes that were perturbed by the treatments, including nucleoside phosphate biosynthetic processes, mRNA metabolic processes, and response to stress. To identify transcripts whose expression level correlated with metal bioaccumulation, we performed Pearson correlation analysis. Several transcripts correlated with gill metal concentrations including mt10, mt20, and contig 48, an unknown transcript containing a wsc domain. In addition, three transcripts directly involved in the unfolded protein response (UPR) were induced in the metal treatments at 2 weeks and were further up-regulated at 4 weeks. Overall, correlation of tissue concentrations and gene expression responses indicates that as mussels accumulate higher concentrations of metals, initial stress responses are mobilized to protect tissues. However, given the role of UPR in apoptosis, it serves as an early indicator of stress, which once overwhelmed will result in adverse physiological effects.