Cardiovascular gene expression profiles of dioxin exposure in zebrafish embryos

Toxicity Evaluation and Transcriptome Analysis of Yellowstripe Goby (Mugilogobius chulae) in Response to 2,7-Dibromocarbazole Exposure during Early Development

Polyhalogenated carbazoles (PHCZs) are a class of nitrogen-containing heterocyclic compounds that are widely distributed throughout the marine environment and sediment. These compounds share structural and toxicity similarities with dioxins. However, our understanding of the toxicological effects of PHCZs on marine organisms and their underlying molecular mechanisms remains limited. In this study, we employed the marine model organism Mugilogobius chulae as the experimental subject and selected 2,7-dibromocarbazole (2,7-DBCZ), a compound known for its high toxicity and detection frequency, to conduct both an acute toxicity test and transcriptome analysis on M. chulae embryos. Our findings revealed that the 96 h median lethal concentration (LC50) of 2,7-DBCZ for M. chulae embryos was 174 μg/L, with a median effective concentration (EC50) resulting in pericardial edema deformity of 88.82 μg/L. Transcriptome analysis revealed significant impacts on various systems in M. chulae embryos following exposure to 2,7-DBCZ, including the sensory, cardiovascular, immune, and endocrine systems. Furthermore, this compound perturbed signaling pathways such as phototransduction, protein folding and processing, amino acid metabolism, lipid transport, and exogenous compound metabolism. Notably, transcript abundance of the CYP1A gene associated with the activation of the AhR signaling pathway, similar to dioxin-like compounds, was 18.18 times higher than that in the control group. This observation suggests that M. chulae embryos mount a stress response when exposed to PHCZs. In summary, this study contributes to our understanding of the toxicological implications of PHCZ in marine fish and offers a theoretical foundation for risk assessment and regulatory frameworks for PHCZs in the marine environment.

Comprehensive Characterization of Multitissue Expression Landscape, Co-Expression Networks and Positive Selection in Pikeperch

Promising efforts are ongoing to extend genomics resources for pikeperch (Sander lucioperca), a species of high interest for the sustainable European aquaculture sector. Although previous work, including reference genome assembly, transcriptome sequence, and single-nucleotide polymorphism genotyping, added a great wealth of genomic tools, a comprehensive characterization of gene expression across major tissues in pikeperch still remains an unmet research need. Here, we used deep RNA-Sequencing of ten vital tissues collected in eight animals to build a high-confident and annotated trancriptome atlas, to detect the tissue-specificity of gene expression and co-expression network modules, and to investigate genome-wide selective signatures in the Percidae fish family. Pathway enrichment and protein–protein interaction network analyses were performed to characterize the unique biological functions of tissue-specific genes and co-expression modules. We detected strong functional correlations and similarities of tissues with respect to their expression patterns—but also significant differences in the complexity and composition of their transcriptomes. Moreover, functional analyses revealed that tissue-specific genes essentially play key roles in the specific physiological functions of the respective tissues. Identified network modules were also functionally coherent with tissues’ main physiological functions. Although tissue specificity was not associated with positive selection, several genes under selection were found to be involved in hypoxia, immunity, and gene regulation processes, that are crucial for fish adaption and welfare. Overall, these new resources and insights will not only enhance the understanding of mechanisms of organ biology in pikeperch, but also complement the amount of genomic resources for this commercial species.

PCB and TCDD derived embryonic cardiac defects result from a novel AhR pathway

Polychlorinated biphenyls (PCBs) and 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) are environmental contaminants known to impact cardiac development, a key step in the embryonic development of most animals. To date, little is understood of the molecular mechanism driving the observed cardiac defects in exposed fishes. The literature shows PCB & TCDD derived cardiac defects are concurrent with, but not caused by, expression of cyp1A, due to activation of the aryl hydrocarbon receptor (AhR) gene activation pathway. However, in this study, detailed visualization of fish hearts exposed to PCBs and TCDD show that, in addition to a failure of cardiac looping in early heart development, the inner endocardial lining of the heart fails to maintain proper cell adhesion and tissue integrity. The resulting gap between the endocardium and myocardium in both zebrafish and Atlantic sturgeon suggested functional faults in endothelial adherens junction formation. Thus, we explored the molecular mechanism triggering cardiac defects using immunohistochemistry to identify the location and phosphorylation state of key regulatory and adhesion molecules. We hypothesized that PCB and TCDD activates AhR, phosphorylating Src, which then phosphorylates the endothelial adherens junction protein, VEcadherin. When phosphorylated, VEcadherin dimers, found in the endocardium and vasculature, separate, reducing tissue integrity. In zebrafish, treatment with PCB and TCDD contaminants leads to higher phosphorylation of VEcadherin in cardiac tissue suggesting that these cells have reduced connectivity. Small molecule inhibition of Src phosphorylation prevents contaminant stimulated phosphorylation of VEcadherin and rescues both cardiac function and gross morphology. Atlantic sturgeon hearts show parallels to contaminant exposed zebrafish cardiac phenotype at the tissue level. These data suggest that the mechanism for PCB and TCDD action in the heart is, in part, distinct from the canonical mechanism described in the literature and that cardiac defects are impacted by this nongenomic mechanism.

Identification of molecular toxicity pathways across early life-stages of zebrafish exposed to PCB126 using a whole transcriptomics approach

Although withdrawn from the market in the 1980s, polychlorinated biphenyls (PCBs) are still found ubiquitously in the aquatic environment and pose a serious risk to biota due to their teratogenic potential. In fish, early life-stages are often considered most sensitive with regard to their exposure to PCBs and other dioxin-like compounds. However, little is known about the molecular drivers of the frequently observed teratogenic effects. Therefore, the aims of our study were to: (1) characterize the baseline transcriptome profiles at different embryonic life-stages in zebrafish (Danio rerio); and (2) to identify the molecular response to PCB exposure and life-stage specific-effects of the chemical on associated processes. For both objectives, embryos were sampled at 12, 48, and 96 h post-fertilization (hpf) and subjected to Illumina sequence-by-synthesis and RNAseq analysis. Results revealed that with increasing age more genes and related pathways were upregulated both in terms of number and magnitude. Yet, other transcripts followed an opposite pattern with greater transcript abundance at the earlier time points. Additionally, embryos were exposed to PCB126, a potent agonist of the aryl hydrocarbon receptor (AHR). ClueGO network analysis revealed significant enrichment of genes associated with basic cell metabolism, communication, and homeostasis as well as eye development, muscle formation, and skeletal formation. We selected eight genes involved in the affected pathways for an in-depth characterization of their regulation throughout normal embryogenesis and after exposure to PCB126 by quantification of transcript abundances every 12 h until 118 hpf. Among these, fgf7 and c9 stood out because of their strong upregulation by PCB126 exposure at 48 and 96 hpf, respectively. Cyp2aa12 was upregulated from 84 hpf on. Fabp10ab, myhz1.1, col8a1a, sulf1, and opn1sw1 displayed specific regulation depending on the developmental stage. Overall, we demonstrate that (1) the developmental transcriptome of zebrafish is highly dynamic, and (2) dysregulation of gene expression by exposure to PCB126 was significant and in several cases not directly connected to AHR-signaling. Hence, this study improves the understanding of linkages between molecular events and apical outcomes that are of regulatory relevance.

A Review of the Functional Roles of the Zebrafish Aryl Hydrocarbon Receptors

Over the last 2 decades, the zebrafish (Danio rerio) has emerged as a stellar model for unraveling molecular signaling events mediated by the aryl hydrocarbon receptor (AHR), an important ligand-activated receptor found in all eumetazoan animals. Zebrafish have 3 AHRs-AHR1a, AHR1b, and AHR2, and studies have demonstrated the diversity of both the endogenous and toxicological functions of the zebrafish AHRs. In this contemporary review, we first highlight the evolution of the zebrafish ahr genes, and the characteristics of the receptors including developmental and adult expression, their endogenous and inducible roles, and the predicted ligands from homology modeling studies. We then review the toxicity of a broad spectrum of AHR ligands across multiple life stages (early stage, and adult), discuss their transcriptomic and epigenetic mechanisms of action, and report on any known interactions between the AHRs and other signaling pathways. Through this article, we summarize the promising research that furthers our understanding of the complex AHR pathway through the extensive use of zebrafish as a model, coupled with a large array of molecular techniques. As much of the research has focused on the functions of AHR2 during development and the mechanism of TCDD (2,3,7,8-tetrachlorodibenzo-p-dioxin) toxicity, we illustrate the need to address the considerable knowledge gap in our understanding of both the mechanistic roles of AHR1a and AHR1b, and the diverse modes of toxicity of the various AHR ligands.

Retene, pyrene and phenanthrene cause distinct molecular-level changes in the cardiac tissue of rainbow trout (Oncorhynchus mykiss) larvae, part 1 - Transcriptomics

Polycyclic aromatic hydrocarbons (PAHs) are contaminants of concern that impact every sphere of the environment. Despite several decades of research, their mechanisms of toxicity are still poorly understood. This study explores the mechanisms of cardiotoxicity of the three widespread model PAHs retene, pyrene and phenanthrene in the rainbow trout (Oncorhynchus mykiss) early life stages. Newly hatched larvae were exposed to each individual compound at sublethal doses causing no significant increase in the prevalence of deformities. Changes in the cardiac transcriptome were assessed after 1, 3, 7 and 14 days of exposure using custom Salmo salar microarrays. The highest number of differentially expressed genes was observed after 1 or 3 days of exposure, and retene was the most potent compound in that regard. Over-representation analyses suggested that genes related to cardiac ion channels, calcium homeostasis and muscle contraction (actin binding, troponin and myosin complexes) were especially targeted by retene. Pyrene was also able to alter similar myosin-related genes, but at a different timing and in an opposite direction, suggesting compound-specific mechanisms of toxicity. Pyrene and to a lesser extent phenanthrene were altering key genes linked to the respiratory electron transport chain and to oxygen and iron metabolism. Overall, phenanthrene was not very potent in inducing changes in the cardiac transcriptome despite being apparently metabolized at a slower rate than retene and pyrene. The present study shows that exposure to different PAHs during the first few days of the swim-up stage can alter the expression of key genes involved into the cardiac development and function, which could potentially affect negatively the fitness of the larvae in the long term.

Phase 0 of the Xenobiotic Response: Nuclear Receptors and Other Transcription Factors as a First Step in Protection from Xenobiotics

This mini-review examines the crucial importance of transcription factors as a first line of defense in the detoxication of xenobiotics. Key transcription factors that recognize xenobiotics or xenobiotic-induced stress such as reactive oxygen species (ROS), include AhR, PXR, CAR, MTF, Nrf2, NF-κB, and AP-1. These transcription factors constitute a significant portion of the pathways induced by toxicants as they regulate phase I-III detoxication enzymes and transporters as well as other protective proteins such as heat shock proteins, chaperones, and anti-oxidants. Because they are often the first line of defense and induce phase I-III metabolism, could these transcription factors be considered the phase 0 of xenobiotic response?

Evaluation of the developmental toxicity of 2,7-dibromocarbazole to zebrafish based on transcriptomics assay

Polyhalogenated carbazoles (PHCZs), which have the similar structure of dioxin, have been reported ubiquitous in the environments and drawn wide concerns. However, their potential ecological and health risks are still poorly understood. Here, wildtype zebrafish embryos were used to evaluate the environmental risks of 2,7-dibromocarbazole (2,7-DBCZ), 3,6-dibromocarbazole (3,6-DBCZ), and 3,6-dichlorocarbazole (3,6-DCCZ). 2,7-DBCZ was the most toxic compound with the 96-h LC₅₀ value of 581.8 ± 29.3 μg·L^-1 and the EC₅₀ value of 201.5 ± 6.5 μg·L^-1 for pericardial edema. The teratogenic effects of 2,7-DBCZ were tested using transgenic zebrafish larvae. The transcriptomic analysis revealed that 90 genes in zebrafish expressed differently after exposure to 2,7-DBCZ, and many pathways were related to aryl hydrocarbon receptor (AhR) activation. The qRT-PCR also showed that expression levels of AhR1 and CYP1 A in zebrafish were significantly up-regulated after exposure to 2,7-DBCZ. In conclusion, 2,7-DBCZ exhibited more potent toxicity and cardiac teratogenic effects, and presented developmental toxicity partially consistent with AhR activation. Our results will be of great help to the risk assessment and regulation-making of PHCZs. Meanwhile, further studies should be promoted to illustrate the potential mechanism between PHCZs and AhR in the near future.

Fishing for contaminants: identification of three mechanism specific transcriptome signatures using Danio rerio embryos

In ecotoxicology, transcriptomics is an effective way to detect gene expression changes in response to environmental pollutants. Such changes can be used to identify contaminants or contaminant classes and can be applied as early warning signals for pollution. To do so, it is important to distinguish contaminant-specific transcriptomic changes from genetic alterations due to general stress. Here we present a first step in the identification of contaminant class-specific transcriptome signatures. Embryos of zebrafish (Danio rerio) were exposed to three substances (methylmercury, chlorpyrifos and Aroclor 1254, each from 24 to 48 hpf exposed) representing sediment typical contaminant classes. We analyzed the altered transcriptome to detect discriminative genes significantly regulated in reaction to the three applied contaminants. By comparison of the results of the three contaminants, we identified transcriptome signatures and biologically important pathways (using Cytoscape/ClueGO software) that react significantly to the contaminant classes. This approach increases the chance of finding genes that play an important role in contaminant class-specific pathways rather than more general processes.

Histopathological changes in zebrafish embryos exposed to DLPCBs extract from Zhanjiang coastal sediment

Dioxin-like polychlorinated biphenyls (DLPCBs) are ubiquitous persistent pollutants that cause adverse effects in many environmental organisms. DLPCBs in marine sediments can be absorbed by benthic organisms, bioaccumulate, and biomagnify through the food chain and threaten animal and human health. There are no reports of DLPCBs concentrations in the Zhanjiang Gulf seabed. This study was designed to investigate the concentration of DLPCBs in the Zhanjiang coastal sediment and histopathological changes in zebrafish (Diano rerio) embryos exposed to environmentally relevant concentrations of DLPCBs. Of the five sites selected, two sites TS and JSW contained DLPCBs at concentrations of 0.08 and 22.54 ng/g dry sediment, respectively. Two groups of zebrafish embryos were used. One group was exposed to 3.75, 7.5, 15, 30, and 60 mg/ml of DLPCBs extracted from the sediments sampled from the TS site and the second group to 4.375, 8.75, 17.5, 35, and 70 mg/ml of DLPCBs from JSW site from 0.75 h post-fertilization (hpf) to 96 hpf. The zebrafish exposed to 60 and 70 mg/ml of DLPCBs at 96 hpf displayed gross histopathological changes with cardiac lesions including pericardial edema being the most deleterious. Other changes observed were hydropic degeneration of gill filaments and hepatocytes, loss of intestinal folds, and uninflated swim bladder. It appears that only a few sites of the Zhanjiang gulf are contaminated with DLPCBs. This is the first report of histopathological changes in the gills, hepatocytes, intestines, heart, and the swim bladder in zebrafish embryos exposed to DLPCBs from a coastal sediment. Further studies with sampling at different stages of development are required to identify which organ/tissue is most sensitive to DLPCBs.

Combined Transcriptomic and Proteomic Approach to Identify Toxicity Pathways in Early Life Stages of Japanese Medaka (Oryzias latipes) Exposed to 1,2,5,6-Tetrabromocyclooctane (TBCO)

Currently, the novel brominated flame retardant 1,2,5,6-tetrabromocyclooctane (TBCO) is considered a potential replacement for hexabromocyclododecane (HBCD). Therefore, use of TBCO could increase in the near future. To assess potential toxicological risks to aquatic organisms, embryos of Japanese medaka (Oryzias latipes) were exposed to 10, 100, or 1000 μg/L TBCO from 2 h postfertilization until 1 day post-hatch. TBCO accumulated in embryos in the order of 0.43-1.3 × 10(4)-fold, and the rate constant of accumulation was 1.7-1.8 per day. The number of days to hatch and the hatching success of embryos exposed to the medium and the greatest concentrations of TBCO were impaired. Responses of the transcriptome (RNA-seq) and proteome were characterized in embryos exposed to 100 μg/L TBCO because this was the least concentration of TBCO that caused an effect on hatching. Consistent with effects on hatching, proteins whose abundances were reduced by exposure to TBCO were enriched in embryo development and hatching pathways. Also, on the basis of the responses of transcriptome and proteome, it was predicted that TBCO might impair vision and contraction of cardiac muscle, respectively, and these effects were confirmed by targeted bioassays. This study provided a comprehensive understanding of effects of TBCO on medaka at early life stages and illustrated the power of "omics" to explain and predict phenotypic responses to chemicals.

Time- and Oil-Dependent Transcriptomic and Physiological Responses to Deepwater Horizon Oil in Mahi-Mahi (Coryphaena hippurus) Embryos and Larvae

The Deepwater Horizon (DWH) oil spill contaminated the spawning habitats for numerous commercially and ecologically important fishes. Exposure to the water accommodated fraction (WAF) of oil from the spill has been shown to cause cardiac toxicity during early developmental stages across fishes. To better understand the molecular events and explore new pathways responsible for toxicity, RNA sequencing was performed in conjunction with physiological and morphological assessments to analyze the time-course (24, 48, and 96 h post fertilization (hpf)) of transcriptional and developmental responses in embryos/larvae of mahi-mahi exposed to WAF of weathered (slick) and source DWH oils. Slick oil exposure induced more pronounced changes in gene expression over time than source oil exposure. Predominant transcriptomic responses included alteration of EIF2 signaling, steroid biosynthesis, ribosome biogenesis and activation of the cytochrome P450 pathway. At 96 hpf, slick oil exposure resulted in significant perturbations in eye development and peripheral nervous system, suggesting novel targets in addition to the heart may be involved in the developmental toxicity of DHW oil. Comparisons of changes of cardiac genes with phenotypic responses were consistent with reduced heart rate and increased pericardial edema in larvae exposed to slick oil but not source oil.

Developmental toxicity, EROD, and CYP1A mRNA expression in zebrafish embryos exposed to dioxin-like PCB126

Dioxin-like PCB126 is a persistent organic pollutant that causes a range of syndromes including developmental toxicity. Dioxins have a high affinity for aryl hydrocarbon receptor (AhR) and induce cytochrome P4501A (CYP1A). However, the role of CYP1A activity in developmental toxicity is less clear. To better understand dioxin induced developmental toxicity, we exposed zebrafish (Danio rerio) embryos to PCB126 at concentrations of 0, 16, 32, 64, and 128 μg L(-1) from 3-h post-fertilization (hpf) to 168 hpf. The embryonic survival rate decreased at 144 and 168 hpf. The fry at 96 hpf displayed gross developmental malformations, including pericardial and yolk sac edema, spinal curvature, abnormal lower jaw growth, and non-inflated swim bladder. The pericardial and yolk sac edema rate significantly increased and the heart rate declined from 96 hpf compared with the controls. PCB126 did not alter the hatching rate. To elucidate the mechanism of PCB126-induced developmental toxicity, we conducted ethoxyresorufin-O-deethylase (EROD) in vivo assay to determine CYP1A enzyme activity, and real-time PCR to study the induction of CYP1A mRNA gene expression in embryo/larval zebrafish at 24, 72, 96, and 132 hpf. In vivo EROD activity was induced by PCB126 at 16 μg L(-1) concentration as early as 72 hpf but significant increases were observed only in zebrafish exposed to 64 and 128 μg L(-1) doses (p < 0.005) at 72, 96, and 132 hpf. Induction of CYP1A mRNA expression was significantly upregulated in zebrafish exposed to 32 and 64 μg L(-1) at 24, 72, 96, and 132 hpf. Overall, the severe pericardial and yolk sac edema and reduced heart rate suggest that heart defects are a sensitive endpoint, and the general trend of dose-dependent increase in EROD activity and induction of CYP1A mRNA gene expression provide evidence that the developmental toxicity of PCB126 to zebrafish embryos is mediated by activation of AhR.

Retene causes multifunctional transcriptomic changes in the heart of rainbow trout (Oncorhynchus mykiss) embryos

Fish are particularly sensitive to aryl hydrocarbon receptor (AhR)-mediated developmental toxicity. The molecular mechanisms behind these adverse effects have remained largely unresolved in salmonids, and for AhR-agonistic polycyclic aromatic hydrocarbons (PAHs). This study explored the cardiac transcriptome of rainbow trout (Oncorhynchus mykiss) eleuteroembryos exposed to retene, an AhR-agonistic PAH. The embryos were exposed to retene (nominal concentration 32 μg/L) and control, their hearts were collected before, at and after the onset of the visible signs of developmental toxicity, and transcriptomic changes were studied by microarray analysis. Retene up- or down-regulated 122 genes. The largest Gene Ontology groups were signal transduction, transcription, apoptosis, cell growth, cytoskeleton, cell adhesion/mobility, cardiovascular development, xenobiotic metabolism, protein metabolism, lipid metabolism and transport, and amino acid metabolism. Together these findings suggest that retene affects multiple signaling cascades in the heart of rainbow trout embryos, and potentially disturbs processes related to cardiovascular development and function.

Ah Receptor Activation by Dioxin Disrupts Activin, BMP, and WNT Signals During the Early Differentiation of Mouse Embryonic Stem Cells and Inhibits Cardiomyocyte Functions

The AHR is a ligand-activated transcription factor that mediates gene-environment interactions. Genome-wide expression profiling during differentiation of mouse ES cells into cardiomyocytes showed that AHR activation by 2,3,7,8-tetrachlorodibenzo-p-dioxin; Dioxin (TCDD), its prototypical ligand, disrupted the expression of multiple homeobox transcription factors and inhibited cardiomyocyte contractility. Here we treated ES cells with TCDD at daily differentiation intervals to investigate whether TCDD-induced loss of contractility had a developmental window of sensitivity. Surprisingly, contractility was an AHR-dependent TCDD target solely between differentiation days 0 and 3 during the period of panmesoderm development, when TCDD also disrupted expression of genes in the TGFβ/BMP2/4 and wingless-type MMTV integration site (WNT)signaling pathways, suppressed the secretion of bone morphogenetic protein (BMP4), WNT3a, and WNT5a and elevated the secretion of Activin A, as determined by ELISA of the secreted proteins in the culture medium. Supplementing the culture medium with BMP4, WNT3a, or WNT5a during the first 3 days of differentiation successfully countered TCDD-induced impairment of contractility, while anti-WNT3a, or anti-WNT5a antibodies or continuous Noggin (a BMP4 antagonist) or Activin A treatment inhibited the contractile phenotype. In Ahr(+/+), but not in Ahr(-) (/) (-) ES cells, TCDD treatment significantly increased mitochondrial copy number, suggestive of mitochondrial stress and remodeling. Sustained AHR activation during ES cell differentiation appears to disrupt the expression of signals critical to the ontogeny of cardiac mesoderm and cause the loss of contractility in the resulting cardiomyocyte lineage.

Adaptation of a Bioinformatics Microarray Analysis Workflow for a Toxicogenomic Study in Rainbow Trout

Sex steroids play a key role in triggering sex differentiation in fish, the use of exogenous hormone treatment leading to partial or complete sex reversal. This phenomenon has attracted attention since the discovery that even low environmental doses of exogenous steroids can adversely affect gonad morphology (ovotestis development) and induce reproductive failure. Modern genomic-based technologies have enhanced opportunities to find out mechanisms of actions (MOA) and identify biomarkers related to the toxic action of a compound. However, high throughput data interpretation relies on statistical analysis, species genomic resources, and bioinformatics tools. The goals of this study are to improve the knowledge of feminisation in fish, by the analysis of molecular responses in the gonads of rainbow trout fry after chronic exposure to several doses (0.01, 0.1, 1 and 10 μg/L) of ethynylestradiol (EE2) and to offer target genes as potential biomarkers of ovotestis development. We successfully adapted a bioinformatics microarray analysis workflow elaborated on human data to a toxicogenomic study using rainbow trout, a fish species lacking accurate functional annotation and genomic resources. The workflow allowed to obtain lists of genes supposed to be enriched in true positive differentially expressed genes (DEGs), which were subjected to over-representation analysis methods (ORA). Several pathways and ontologies, mostly related to cell division and metabolism, sexual reproduction and steroid production, were found significantly enriched in our analyses. Moreover, two sets of potential ovotestis biomarkers were selected using several criteria. The first group displayed specific potential biomarkers belonging to pathways/ontologies highlighted in the experiment. Among them, the early ovarian differentiation gene foxl2a was overexpressed. The second group, which was highly sensitive but not specific, included the DEGs presenting the highest fold change and lowest p-value of the statistical workflow output. The methodology can be generalized to other (non-model) species and various types of microarray platforms.

A rapid in vivo zebrafish model to elucidate oxidative stress-mediated PCB126-induced apoptosis and developmental toxicity

Dioxin-like 3,3',4,4',5-pentachlorobiphenyl (PCB126) is one of the most potent and widespread environmental pollutants. Although PCB126-induced toxicity is related to the aryl hydrocarbon receptor pathway, there is still no study that has constructed an in vivo visual model to clarify the role of the Nrf2/ARE signaling pathway in the oxidative stress mechanism of PCB126-induced toxicity. In the present study, an in vivo zebrafish model of nrf2a fused to enhanced green fluorescent protein (nrf2a-eGFP) was constructed. The zebrafish embryos microinjected with nrf2a-eGFP (72h postfertilization) were exposed to various concentrations of PCB126 (0, 25, 50, 100, 200μg/L) or 30mMN-acetylcysteine (NAC)+200μg/L PCB126. After 72h exposure, PCB126 significantly increased the malformation rates and induced eGFP expression in a dose-dependent manner in several zebrafish tissue types. The distribution of eGFP fluorescence coincided with developmental deformity sites. NAC pretreatment effectively counteracted PCB126-induced developmental toxicity including heart rate, pericardial edema, and body length. The highest PCB126 dose, 200μg/L, produced marked apoptosis in the eye, gill, and trunk detected by the terminal deoxynucleotidyl transferase-mediated dUTP nick-end labeling assay. At 48 and 72h exposure, 200μg/L PCB126 affected glutathione metabolism as evidenced by decreased glutathione and increased glutathione disulfide concentrations, indicative of oxidative stress. These effects were also counteracted by NAC pretreatment. Furthermore, the Nrf2-regulated genes gclc, gpx, gstp1, and hmox1 were significantly induced at 24, 48, and 72h at the highest PCB126 exposures but not in the NAC-pretreated group. In addition, a significant increase in ROS generation was detected in zebrafish larvae at 72h PCB126 exposure, which might offer a link for future mechanistic studies. Collectively, these data suggest that PCB126-induced developmental toxicity and apoptosis in the nrf2a-eGFP-injected zebrafish model are due to oxidative stress mediated by disruption to glutathione metabolism and changes in Nrf2-regulated gene expression.

A genome scan for candidate genes involved in the adaptation of turbot (Scophthalmus maximus)

Partitioning phenotypic variance in genotypic and environmental variance may benefit from the population genomic assignment of genes putatively involved in adaptation. We analyzed a total of 256 markers (120 microsatellites and 136 Single Nucleotide Polymorphisms - SNPs), several of them associated to Quantitative Trait Loci (QTL) for growth and resistance to pathologies, with the aim to identify potential adaptive variation in turbot Scophthalmus maximus L. The study area in the Northeastern Atlantic Ocean, from Iberian Peninsula to the Baltic Sea, involves a gradual change in temperature and an abrupt change in salinity conditions. We detected 27 candidate loci putatively under selection. At least four of the five SNPs identified as outliers are located within genes coding for ribosomal proteins or directly related with the production of cellular proteins. One of the detected outliers, previously identified as part of a QTL for growth, is a microsatellite linked to a gene coding for a growth factor receptor. A similar set of outliers was detected when natural populations were compared with a sample subjected to strong artificial selection for growth along four generations. The observed association between FST outliers and growth-related QTL supports the hypothesis of changes in growth as an adaptation to differences in temperature and salinity conditions. However, further work is needed to confirm this hypothesis.

Understanding different facets of cardiovascular diseases based on model systems to human studies: a proteomic and metabolomic perspective

Cardiovascular disease has remained as the largest cause of morbidity and mortality worldwide. From dissecting the disease aetiology to identifying prognostic markers for better management of the disease is still a challenge for researchers. In the post human genome sequencing era much of the thrust has been focussed towards application of advanced genomic tools along with evaluation of traditional risk factors. With the advancement of next generation proteomics and metabolomics approaches it has now become possible to understand the protein interaction network & metabolic rewiring which lead to the perturbations of the disease phenotype. Further, elucidating different post translational modifications using advanced mass spectrometry based methods have provided an impetus towards in depth understanding of the proteome. The past decade has observed a plethora of studies where proteomics has been applied successfully to identify potential prognostic and diagnostic markers as well as to understand the disease mechanisms for various types of cardiovascular diseases. In this review, we attempted to document relevant proteomics based studies that have been undertaken either to identify potential biomarkers or have elucidated newer mechanistic insights into understanding the patho-physiology of cardiovascular disease, primarily coronary artery disease, cardiomyopathy, and myocardial ischemia. We have also provided a perspective on the potential of proteomics in combating this deadly disease.

BIOLOGICAL SIGNIFICANCE

This review has catalogued recent studies on proteomics and metabolomics involved in understanding several cardiovascular diseases (CVDs). A holistic systems biology based approach, of which proteomics and metabolomics are two very important components, would help in delineating various pathways associated with complex disorders like CVD. This would ultimately provide better mechanistic understanding of the disease biology leading to development of prognostic biomarkers. This article is part of a Special Issue entitled: Proteomics in India.

Transcriptional and physiological response of fathead minnows (Pimephales promelas) exposed to urban waters entering into wildlife protected areas

The mission of protected areas is to conserve biodiversity and improve human welfare. To assess the effect of urban waters entering into protected areas, we performed 48-h whole-effluent exposures with fathead minnows, analyzing changes in steady state levels of mRNAs in the livers of exposed fish. Raw wastewater, treated city wastewater, and treated wastewater from a university were collected for exposures. All exposed fish showed altered mRNA levels of DNA damage-repair genes. Fish exposed to raw and treated wastewaters showed down-regulation of transcripts for key intermediates of cholesterol biosynthesis and elevated plasma cholesterol. The type of wastewater treatment influenced the response of gene transcription. Because of the relevance of some of the altered cellular pathways, we suggest that these effluents may cause deleterious effects on fish inside protected areas that receive these waters. Inclusion of research and mitigation efforts for this type of threat in protected areas management is advised.

Aryl hydrocarbon receptor activation and developmental toxicity in zebrafish in response to soil extracts containing unsubstituted and oxygenated PAHs

Many industrial sites are polluted by complex mixtures of polycyclic aromatic compounds (PACs). Besides polycyclic aromatic hydrocarbons (PAHs), these mixtures often contain significant amounts of more polar PACs including oxygenated PAHs (oxy-PAHs). The effects of oxy-PAHs are, however, poorly known. Here we used zebrafish embryos to examine toxicities and transcriptional changes induced by PAC containing soil extracts from three different industrial sites: a gasworks (GAS), a former wood preservation site (WOOD), and a coke oven (COKE), and to PAH and oxy-PAH containing fractions of these. All extracts induced aryl hydrocarbon receptor (Ahr)-regulated mRNAs, malformations, and mortality. The WOOD extract was most toxic and the GAS extract least toxic. The extracts induced glutathione transferases and heat shock protein 70, suggesting that the toxicity also involved oxidative stress. With all extracts, Ahr2-knock-down reduced the toxicity, indicating a significant Ahr2-dependence on the effects. Ahr2-knock-down was most effective with the PAH fraction of the WOOD extract and with the oxy-PAH fraction of the COKE extract. Our results indicate that oxy-PAH containing mixtures can be as potent Ahr activators and developmental toxicants as PAHs. In addition to Ahr activating potency, the profile of cytochrome P4501 inhibitors may also determine the toxic potency of the extracts.

Role of pregnane X receptor and aryl hydrocarbon receptor in transcriptional regulation of pxr, CYP2, and CYP3 genes in developing zebrafish

Ligand-activated receptors regulate numerous genes, and mediate effects of a broad set of endogenous and exogenous chemicals in vertebrates. Understanding the roles of these transcription factors in zebrafish (Danio rerio) is important to the use of this non-mammalian model in toxicological, pharmacological, and carcinogenesis research. Response to a potential agonist for the pregnane X receptor (Pxr) [pregnenolone (PN)] was examined in developing zebrafish, to assess involvement of Pxr in regulation of selected genes, including genes in cytochrome P450 subfamilies CYP2 and CYP3. We also examined interaction of Pxr and the aryl hydrocarbon receptor (Ahr) signaling pathways. Pregnenolone caused a dose-dependent increase in mRNA levels of pxr, ahr2, CYP1A, CYP2AA1, CYP2AA12, CYP3A65, and CYP3C1, most of which peaked at 3 µM PN. The well-known Ahr agonist 3,3',4,4',5-pentachlorobiphenyl (PCB126) also upregulated expression of pxr, ahr2, CYP1A, CYP2AA12, CYP3A65, and CYP3C1 in a dose-dependent manner. Inhibition of pxr translation by morpholino antisense oligonucleotides (MO) suppressed PN-induced expression of pxr, ahr2, CYP3A65, and CYP3C1 genes. Levels of CYP2AA1 and CYP2AA12 mRNA were increased in the control-MO group exposed to PN; this was prevented by knocking down Pxr. Similarly, Ahr2-MO treatment blocked PCB126-induced mRNA expression of pxr, CYP1A, CYP2AA12, CYP3A65, and CYP3C1. The present study shows self-regulation of pxr by PN in developing zebrafish. Selected zebrafish CYP1, CYP2 (including several CYP2AAs) and CYP3 genes appear to be under the regulation of both Pxr and Ahr2.

AHR2-Mediated transcriptomic responses underlying the synergistic cardiac developmental toxicity of PAHs

Polycyclic aromatic hydrocarbons (PAHs) induce developmental defects including cardiac deformities in fish. The aryl hydrocarbon receptor (AHR) mediates the toxicity of some PAHs. Exposure to a simple PAH mixture during embryo development consisting of an AHR agonist (benzo(a)pyrene-BaP) with fluoranthene (FL), an inhibitor of cytochrome p450 1(CYP1)--a gene induced by AHR activation--results in cardiac deformities. Exposure to BaP or FL alone at similar concentrations alters heart rates, but does not induce morphological deformities. Furthermore, AHR2 knockdown prevents the toxicity of BaP + FL mixture. Here, we used a zebrafish microarray analysis to identify heart-specific transcriptomic changes during early development that might underlie cardiotoxicity of BaP + FL. We used AHR2 morphant embryos to determine the role of this receptor in mediating toxicity. Control and knockdown embryos at 36 h post-fertilization were exposed to DMSO, 100 μg/l BaP, 500 μg/l FL, or 100 μg/l BaP + 500 μg/l FL, and heart tissues for RNA were extracted at 2, 6, 12, and 18 h-post-exposure (hpe), prior to the appearance of cardiac deformities. Data show AHR2-dependent BaP + FL effects on expression of genes involved in protein biosynthesis and neuronal development in addition to signaling molecules and their associated molecular pathways. Ca(2+)-cycling and muscle contraction genes were the most significantly differentially expressed category of transcripts when comparing BaP + FL-treated AHR2 morphant and control embryos. These differences were most prominent at 2 and 6 hpe. Therefore, we postulate that BaP + FL may affect cellular Ca(2+) levels and subsequently cardiac muscle function, potentially underlying BaP + FL cardiotoxicity.

The transcriptional response to oxidative stress during vertebrate development: effects of tert-butylhydroquinone and 2,3,7,8-tetrachlorodibenzo-p-dioxin

Oxidative stress is an important mechanism of chemical toxicity, contributing to teratogenesis and to cardiovascular and neurodegenerative diseases. Developing animals may be especially sensitive to chemicals causing oxidative stress. The developmental expression and inducibility of anti-oxidant defenses through activation of NF-E2-related factor 2 (NRF2) affect susceptibility to oxidants, but the embryonic response to oxidants is not well understood. To assess the response to chemically mediated oxidative stress and how it may vary during development, zebrafish embryos, eleutheroembryos, or larvae at 1, 2, 3, 4, 5, and 6 days post fertilization (dpf) were exposed to DMSO (0.1%), tert-butylhydroquinone (tBHQ; 10 µM) or 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD; 2 nM) for 6 hr. Transcript abundance was assessed by real-time qRT-PCR and microarray. qRT-PCR showed strong (4- to 5-fold) induction of gstp1 by tBHQ as early as 1 dpf. tBHQ also induced gclc (2 dpf), but not sod1, nqo1, or cyp1a. TCDD induced cyp1a but none of the other genes. Microarray analysis showed that 1477 probes were significantly different among the DMSO-, tBHQ-, and TCDD-treated eleutheroembryos at 4 dpf. There was substantial overlap between genes induced in developing zebrafish and a set of marker genes induced by oxidative stress in mammals. Genes induced by tBHQ in 4-dpf zebrafish included those involved in glutathione synthesis and utilization, signal transduction, and DNA damage/stress response. The strong induction of hsp70 determined by microarray was confirmed by qRT-PCR and by use of transgenic zebrafish expressing enhanced green fluorescent protein (EGFP) under control of the hsp70 promoter. Genes strongly down-regulated by tBHQ included mitfa, providing a molecular explanation for the loss of pigmentation in tBHQ-exposed embryos. These data show that zebrafish embryos are responsive to oxidative stress as early as 1 dpf, that responsiveness varies with development in a gene-specific manner, and that the oxidative stress response is substantially conserved in vertebrate animals.

Dioxin-induced thrombocyte aggregation in zebrafish

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a canonical member of a group of dioxins which are byproducts of industrial combustion and are dangerous environmental pollutants. TCDD has been shown to cause several abnormalities in humans and wildlife, and recently, some dioxins have been found to activate platelets. However, TCDD-mediated platelet activation pathways are elusive and virtually nothing is known about TCDD activation of fish thrombocytes. To investigate TCDD effect on thrombocyte function, we tested zebrafish blood in presence of TCDD using a thrombocyte functional assay. We found that TCDD activated thrombocytes. Further experiments showed that thrombocytes of fish treated with TCDD formed both aggregates and filopodia. To investigate the mechanism of TCDD-mediated activation of thrombocytes we used inhibitors for Gq, cyclooxygenase-1, aryl hydrocarbon receptor (AHR), c-src, Akt, and ERK1/2. We found that TCDD induces AHR which activates c-src and signals the activation of Akt and ERK1/2 which are ultimately involved in generation of thromboxane A2. Furthermore, we found that ADP potentiates TCDD action, which led to the discovery that ADP itself activates AHR in the absence of TCDD. Taken together, these results resolved the pathway of TCDD activation of thrombocytes and led to the finding that ADP is an activator of AHR.

Identification, modeling and ligand affinity of early deuterostome CYP51s, and functional characterization of recombinant zebrafish sterol 14α-demethylase

BACKGROUND

Sterol 14α-demethylase (cytochrome P450 51, CYP51, P45014DM) is a microsomal enzyme that in eukaryotes catalyzes formation of sterols essential for cell membrane function and as precursors in biosynthesis of steroid hormones. Functional properties of CYP51s are unknown in non-mammalian deuterostomes.

METHODS

PCR-cloning and sequencing and computational analyses (homology modeling and docking) addressed CYP51 in zebrafish Danio rerio, the reef fish sergeant major Abudefduf saxatilis, and the sea urchin Strongylocentrotus purpuratus. Following N-terminal amino acid modification, zebrafish CYP51 was expressed in Escherichia coli, and lanosterol 14α-demethylase activity and azole inhibition of CYP51 activity were characterized using GC-MS.

RESULTS

Molecular phylogeny positioned S. purpuratus CYP51 at the base of the deuterostome clade. In zebrafish, CYP51 is expressed in all organs examined, most strongly in intestine. The recombinant protein bound lanosterol and catalyzed 14α-demethylase activity, at 3.2nmol/min/nmol CYP51. The binding of azoles to zebrafish CYP51 gave KS (dissociation constant) values of 0.26μM for ketoconazole and 0.64μM for propiconazole. Displacement of carbon monoxide also indicated zebrafish CYP51 has greater affinity for ketoconazole. Docking to homology models showed that lanosterol docks in fish and sea urchin CYP51s with an orientation essentially the same as in mammalian CYP51s. Docking of ketoconazole indicates it would inhibit fish and sea urchin CYP51s.

CONCLUSIONS

Biochemical and computational analyses are consistent with lanosterol being a substrate for early deuterostome CYP51s.

GENERAL SIGNIFICANCE

The results expand the phylogenetic view of animal CYP51, with evolutionary, environmental and therapeutic implications.

RNA-sequencing analysis of TCDD-induced responses in zebrafish liver reveals high relatedness to in vivo mammalian models and conserved biological pathways

TCDD is one of the most persistent environmental toxicants in biological systems and its effect through aryl hydrocarbon receptor (AhR) has been well characterized. However, the information on TCDD-induced toxicity in other molecular pathways is rather limited. To fully understand molecular toxicity of TCDD in an in vivo animal model, adult zebrafish were exposed to TCDD at 10 nM for 96 h and the livers were sampled for RNA-sequencing based transcriptomic profiling. A total of 1,058 differently expressed genes were identified based on fold-change>2 and TPM (transcripts per million) >10. Among the top 20 up-regulated genes, 10 novel responsive genes were identified and verified by RT-qPCR analysis on independent samples. Transcriptomic analysis indicated several deregulated pathways associated with cell cycle, endocrine disruptors, signal transduction and immune systems. Comparative analyses of TCDD-induced transcriptomic changes between fish and mammalian models revealed that proteomic pathway is consistently up-regulated while calcium signaling pathway and several immune-related pathways are generally down-regulated. Finally, our study also suggested that zebrafish model showed greater similarity to in vivo mammalian models than in vitro models. Our study indicated that the zebrafish is a valuable in vivo model in toxicogenomic analyses for understanding molecular toxicity of environmental toxicants relevant to human health. The expression profiles associated with TCDD could be useful for monitoring environmental dioxin and dioxin-like contamination.

Statistically enhanced spectral counting approach to TCDD cardiac toxicity in the adult zebrafish heart

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a persistent environmental pollutant and teratogen that produces cardiac toxicity in the developing zebrafish. Here we adopted a label free quantitative proteomic approach based on normalized spectral abundance factor (NSAF) to investigate the disturbance of the cardiac proteome induced by TCDD in the adult zebrafish heart. The protein expression level changes between heart samples from TCDD-treated and control zebrafish were systematically evaluated by a large scale MudPIT analysis, which incorporated triplicate analyses for both control and TCDD-exposed heart proteomic samples to overcome the data-dependent variation in shotgun proteomic experiments and obtain a statistically significant protein data set with improved quantification confidence. A total of 519 and 443 proteins were identified in hearts collected from control and TCDD-treated zebrafish, respectively, among which 106 proteins showed statistically significant expression changes. After correcting for the experimental variation between replicate analyses by statistical evaluation, 55 proteins exhibited NSAF ratios above 2 and 43 proteins displayed NSAF ratios smaller than 0.5, with statistical significance by t test (p < 0.05). The proteins identified as altered by TCDD encompass a wide range of biological functions including calcium handling, myocardium cell architecture, energy production and metabolism, mitochondrial homeostasis, and stress response. Collectively, our results indicate that TCDD exposure alters the adult zebrafish heart in a way that could result in cardiac hypertrophy and heart failure and suggests a potential mechanism for the diastolic dysfunction observed in TCDD-exposed embryos.

Developmental effects of aerosols and coal burning particles in zebrafish embryos

Embryo toxicity of particles generated by combustion processes is of special concern for human health. A significant part of these toxic effects is linked to the binding of some pollutants (like polycyclic aromatic hydrocarbons or PAHs) to the Aryl hydrocarbon Receptor (AhR) and the activation of target genes, like the cytochrome P4501A. This activity was analyzed for ambient air and coal-combustion particle extracts in zebrafish embryos (the cyp1aDarT assay) and in two single-cell bioassays: the yeast-based YCM-RYA and the DR-luc (rat cells) assay. Observed AhR ligand activity of samples generally correlated to the predicted toxic effect according to their PAH composition, except for one of the coal combustion samples with an anomalously high activity in the cyp1aDarT assay. This sample induced deformities in zebrafish embryos. We concluded that the combination of morphological and molecular assays may detect embryonic toxic effects that cannot be predicted from chemical analyses or single-cell bioassays.

GFP transgenic medaka (Oryzias latipes) under the inducible cyp1a promoter provide a sensitive and convenient biological indicator for the presence of TCDD and other persistent organic chemicals

Persistent organic pollutants (POPs) are resistant to environmental degradation and can cause multitude of health problems. Cytochrome P450 1A (Cyp1a) is often up-regulated by POPs through the activation of aryl hydrocarbon receptor (AhR) pathway and is thus usually used as a biomarker for xenobiotics exposure. To develop a convenient in vivo tool to monitor xenobiotic contamination in the water, we have established GFP transgenic medaka using the inducible cyp1a promoter, Tg(cyp1a:gfp). Here we tested Tg(cyp1a:gfp) medaka at three different stages, prehatching embryos, newly hatched fry and adult with 2,3,7,8-tetrachlorodiebnzo-p-dioxin (TCDD), a dioxin. While GFP induction was observed in all three stages, newly hatched fry were the most sensitive with the lowest observed effective concentration of 0.005 nM or 16.1 ng/L. The highly sensitive organs included the kidney, liver and intestine. With high concentrations of TCDD, several other organs such as the olfactory pit, tail fin, gills, lateral line neuromast cells and blood vessels also showed GFP expression. In addition, Tg(cyp1a:gfp) medaka fry also responded to two other AhR agonists, 3-methylcholanthrene and benzo[a]pyrene, for GFP induction, but no significant GFP induction was observed towards several other chemicals tested, indicating the specificity of this transgenic line. The GFP inducibility of Tg(cyp1a:gfp) medaka at both fry and adult stages may be useful for development of high-throughput assays as well as online water monitoring system to detect xenobiotic toxicity.

Relative potency of PCB126 to TCDD for sublethal embryotoxicity in the mummichog (Fundulus heteroclitus)

The relative potency (ReP) of 3,3',4,4',5-pentachlorobiphenyl (PCB126) to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) for sublethal responses was assessed in Fundulus heteroclitus embryos. Eggs were treated with intravitelline injections of graded sublethal doses of PCB126 (312-5000 pg g(-1) wet weight, ww) or TCDD (5-1280 pg g(-1) ww). At 16 days post-fertilization (DPF), craniofacial deformities were observed in larvae hatched from eggs treated with the two highest doses of PCB126 (2500-5000 pg g(-1) ww). Both compounds caused a dose-responsive reduction of larval growth and prey capture ability (at ≥1250 pg g(-1) ww), and induction of ethoxyresorufin-O-deethylase (EROD) activity (at ≥80 pg g(-1) ww). The dose-response relationships for EROD activity for PCB126 and TCDD had similar slopes and the ReP of PCB126 to TCDD for EROD activity was estimated at 0.71. This is 140-fold higher than the World Health Organization (WHO) TCDD equivalency factor (TEF) of PCB126 for fish (0.005), which is based on rainbow trout (Oncorhynchus mykiss) embryolethality data. The slope of the dose-response relationship for prey capture ability for PCB126 was steeper than for TCDD, suggesting different mechanisms of action. Expression levels of several genes were also studied by quantitative real-time polymerase chain reaction (qPCR) following exposure to single doses of TCDD or PCB126 (1280 and 1250 pg g(-1) ww, respectively) causing similar EROD induction. A different pattern of responses was observed between PCB126 and TCDD: PCB126 appeared to induce antioxidant responses by inducing sod2 expression, while TCDD did not. These results suggest that relative potencies are species-specific and that the current ReP for PCB126 underestimates its toxicity for some fish species. It is recommended to develop species-specific RePs for a variety of sublethal endpoints and at environmentally relevant doses.

Basal and 3-methylcholanthrene-induced expression of cytochrome P450 1A, 1B and 1C genes in the Brazilian guppy, Poecilia vivipara

In fish there are four cytochrome P450 (CYP1) subfamilies: CYP1A, CYP1B, CYP1C, and CYP1D. Here we cloned Poecilia vivipara CYP1A, with an inferred amino acid sequence 91% identical to CYP1A from the killifish Fundulus heteroclitus, another member of the Cypriniformes, and an important model in ecotoxicology. In addition, we examined the expression of CYP1A, CYP1B1, and CYP1C1 by qPCR in liver, gill, and intestine of adult P. vivipara injected with 3-methylcholanthrene (3-MC) or held in clean water (control group) for 24h. All three tissues examined showed basal expression of the three CYP1 genes. CYP1A was most strongly expressed in the liver, while CYP1B1, and CYP1C1 were most strongly expressed in the gill and intestine respectively. 3-MC induced CYP1A, CYP1B1, and CYP1C1 significantly (20-120-fold) in the three organs, consistent with the regulation of CYP1A, CYP1B1 and CYP1C1 via the aryl hydrocarbon receptor. Validation of CYP1 gene biomarkers in fish collected from a contaminated urban mangrove environment was confirmed with significant induction of CYP1A and CYP1C1 in gills (10-15-fold) and CYP1B1 in liver (23-fold), relative to fish from a control site. The responsiveness of these CYP1 genes indicates P. vivipara is suitable as a model for environmental toxicology studies and environmental assessment in Brazil.

Ahr2-dependence of PCB126 effects on the swim bladder in relation to expression of CYP1 and cox-2 genes in developing zebrafish

The teleost swim bladder is assumed a homolog of the tetrapod lung. Both swim bladder and lung are developmental targets of persistent aryl hydrocarbon receptor (AHR(2)) agonists; in zebrafish (Danio rerio) the swim bladder fails to inflate with exposure to 3,3',4,4',5-pentachlorobiphenyl (PCB126). The mechanism for this effect is unknown, but studies have suggested roles of cytochrome P450 1 (CYP1) and cyclooxygenase 2 (Cox-2) in some Ahr-mediated developmental effects in zebrafish. We determined relationships between swim bladder inflation and CYP1 and Cox-2 mRNA expression in PCB126-exposed zebrafish embryos. We also examined effects on β-catenin dependent transcription, histological effects, and Ahr2 dependence of the effect of PCB126 on swim bladder using morpholinos targeting ahr2. One-day-old embryos were exposed to waterborne PCB126 or carrier (DMSO) for 24h and then held in clean water until day 4, a normal time for swim bladder inflation. The effects of PCB126 were concentration-dependent with EC(50) values of 1.4 to 2.0 nM for induction of the CYP1s, 3.7 and 5.1 nM (or higher) for cox-2a and cox-2b induction, and 2.5 nM for inhibition of swim bladder inflation. Histological defects included a compaction of the developing bladder. Ahr2-morpholino treatment rescued the effect of PCB126 (5 nM) on swim bladder inflation and blocked induction of CYP1A, cox-2a, and cox-2b. With 2nM PCB126 approximately 30% of eleutheroembryos(3) failed to inflate the swim bladder, but there was no difference in CYP1 or cox-2 mRNA expression between those embryos and embryos showing inflated swim bladder. Our results indicate that PCB126 blocks swim bladder inflation via an Ahr2-mediated mechanism. This mechanism seems independent of CYP1 or cox-2 mRNA induction but may involve abnormal development of swim bladder cells.

Cardiosulfa induces heart deformation in zebrafish through the AhR-mediated, CYP1A-independent pathway

Heart development is a complicated and elaborate biological process. To study this and similar complicated process and diseases, the discovery and use of small molecules for probing biological events is invaluable. As part of such an investigation, we have identified cardiosulfa, a small molecule that induces severely impaired heart morphology and function in zebrafish. The results of the present study show that cardiosulfa-promoted heart deformation is protected by negative regulators of the aryl hydrocarbon receptor (AhR) signaling pathway, such as the AhR antagonist CH-223191 and an AhR2-morpholino antisense oligonucleotide, zfahr2-MO. However, the toxic effect of cardiosulfa is not alleviated by zfcyp1a-MO, a morpholino antisense oligo for cytochrome P450 1A (CYP1A), which is the most well-characterized gene of the AhR pathway. Similar results were obtained for the known AhR agonist PCB126. These observations suggest that cardiosulfa causes heart deformation in zebrafish through the AhR-mediated, CYP1A-independent pathway. Our results indicate that cardiosulfa has potential as a novel type of a biological probe to investigate the AhR pathway.

Reproductive and developmental toxicity of dioxin in fish

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD or dioxin) is a global environmental contaminant and the prototypical ligand for investigating aryl hydrocarbon receptor (AHR)-mediated toxicity. Environmental exposure to TCDD results in developmental and reproductive toxicity in fish, birds and mammals. To resolve the ecotoxicological relevance and human health risks posed by exposure to dioxin-like AHR agonists, a vertebrate model is needed that allows for toxicity studies at various levels of biological organization, assesses adverse reproductive and developmental effects and establishes appropriate integrative correlations between different levels of effects. Here we describe the reproductive and developmental toxicity of TCDD in feral fish species and summarize how using the zebrafish model to investigate TCDD toxicity has enabled us to characterize the AHR signaling in fish and to better understand how dioxin-like chemicals induce toxicity. We propose that such studies can be used to predict the risks that AHR ligands pose to feral fish populations and provide a platform for integrating risk assessments for both ecologically relevant organisms and humans.

EZR1: a novel family of highly expressed retroelements induced by TCDD and regulated by a NF-κB-like factor in embryos of zebrafish (Danio rerio)

Transcript profiling using a zebrafish heart cDNA library previously revealed abundant expressed sequence tags (ESTs) upregulated in zebrafish embryos treated with the aryl hydrocarbon receptor (AHR) agonist 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Here, we identify those ESTs as LTR-containing retroelements termed EZR1 (Expressed-Zebrafish-Retroelement group 1). EZR1 is highly redundant in the genome and includes canonical long terminal repeats (LTRs) flanking an integrase-like open reading frame and a region similar to retroviral envelope protein genes. EZR1 sequences lack reverse transcriptase, RNase H, or protease, indicating retrotransposition would be nonautonomous. No AHR binding motifs were found in the EZR1 promoter region. A putative NF-κB-binding site was found, and TCDD-treated zebrafish embryos had significantly increased levels of nuclear protein(s) binding to this sequence. Protein-EZR1 DNA complex formation was partially competed by a mammalian consensus κB sequence, consistent with NF-κB-like activation contributing to increased protein binding to this site. Mobility of the TCDD-induced protein-EZR1 complex differed from that of authentic NF-κB protein bound to the consensus κB site. The results suggest that EZR1 is regulated by interaction with NF-κB or NF-κB-like protein(s) different from the NF-κB protein binding to the consensus κB site. The nature of the NF-κB-like protein and the relationship between EZR1 induction and cardiovascular toxicity caused by TCDD warrant further investigation.

A combined DNA-microarray and mechanism-specific toxicity approach with zebrafish embryos to investigate the pollution of river sediments

The zebrafish embryo has repeatedly proved to be a useful model for the analysis of effects by environmental toxicants. This proof-of-concept study was performed to investigate if an approach combining mechanism-specific bioassays with microarray techniques can obtain more in-depth insights into the ecotoxicity of complex pollutant mixtures as present, e.g., in sediment extracts. For this end, altered gene expression was compared to data from established bioassays as well as to results from chemical analysis. Mechanism-specific biotests indicated a defined hazard potential of the sediment extracts, and microarray analysis revealed several classes of significantly regulated genes which could be related to the hazard potential. Results indicate that potential classes of contaminants can be assigned to sediment extracts by both classical biomarker genes and corresponding expression profile analyses of known substances. However, it is difficult to distinguish between specific responses and more universal detoxification of the organism.

Transcriptional response of hepatic largemouth bass (Micropterus salmoides) mRNA upon exposure to environmental contaminants

Microarrays enable gene transcript expression changes in near-whole genomes to be assessed in response to environmental stimuli. We utilized oligonucleotide microarrays and subsequent gene set enrichment analysis (GSEA) to assess patterns of gene expression changes in male largemouth bass (Micropterus salmoides) hepatic tissues after a 96 h exposure to common environmental contaminants. Fish were exposed to atrazine, cadmium chloride, PCB 126, phenanthrene and toxaphene via intraperitoneal injection with target body burdens of 3.0, 0.00067, 2.5, 50 and 100 µg g(-1), respectively. This was conducted in an effort to identify potential biomarkers of exposure. The expressions of 4, 126, 118, 137 and 58 mRNA transcripts were significantly (P ≤ 0.001, fold change ≥2×) affected by exposure to atrazine, cadmium chloride, PCB 126, phenanthrene and toxaphene exposures, respectively. GSEA revealed that none, four, five, five and three biological function gene ontology categories were significantly influenced by exposure to these chemicals, respectively. We observed that cadmium chloride elicited ethanol metabolism responses, and along with PCB 126 and phenanthrene affected transcripts associated with protein biosynthesis. PCB 126, phenanthrene and toxaphene also influenced one-carbon compound metabolism while PCB 126 and phenanthrene affected mRNA transcription and mRNA export from the nucleus and may have induced an antiestrogenic response. Atrazine was found to alter the expression of few hepatic transcripts. This work has highlighted several biological processes of interest that may be helpful in the development of gene transcript biomarkers of chemical exposure in fish.

Genetic architecture of susceptibility to PCB126-induced developmental cardiotoxicity in zebrafish

Variability in risk of developmental defects caused by dioxin-like compounds (DLCs) has been demonstrated within and among several vertebrate species. Beyond our knowledge of the aryl hydrocarbon receptor (AHR) and its role in mediating toxicity for this class of compounds, little else is known concerning precise downstream targets influencing this vulnerability. In the present study, zebrafish with divergent genetic backgrounds were screened for susceptibility to developmental cardiotoxicity caused by the prototypical DLC, 3,3',4,4',5-pentachlorobiphenyl (PCB126); a range up to ∼40-fold differences was observed. Differentially sensitive zebrafish were chosen for a genetic cross, and the recombinant generation was used for genome-wide quantitative trait loci (QTL) mapping. Multiple QTLs were identified--several acting alone, one additively, and two others via epistatic interaction. Together, these QTLs account for 24% of the phenotypic variance observed in cardioteratogenicity resulting from PCB126 exposure (logarithm of the odds = 13.55, p = 1.89 × 10⁻¹⁰). Candidate genes in these QTL regions include the following: ahr2, bcor, and capn1 (Chr 22); e2f1 and pdyn (Chr 23); ctnnt2, plcg1, eno3, tgm1, and tgm2 (interacting on Chr 23); and vezf1 (Chr 15). These data demonstrate that DLC-induced cardiac teratogenicity is a multifactorial complex trait influenced by gene × gene and gene × environment interactions. The identified QTLs harbor many DLC-responsive genes critical to cardiovascular development and provide insight into the genetic basis of susceptibility to AHR-mediated developmental toxicity.

Role of zebrafish cytochrome P450 CYP1C genes in the reduced mesencephalic vein blood flow caused by activation of AHR2

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) causes various signs of toxicity in early life stages of vertebrates through activation of the aryl hydrocarbon receptor (AHR). We previously reported a sensitive and useful endpoint of TCDD developmental toxicity in zebrafish, namely a decrease in blood flow in the dorsal midbrain, but downstream genes involved in the effect are not known. The present study addressed the role of zebrafish cytochrome P450 1C (CYP1C) genes in association with a decrease in mesencephalic vein (MsV) blood flow. The CYP1C subfamily was recently discovered in fish and includes the paralogues CYP1C1 and CYP1C2, both of which are induced via AHR2 in zebrafish embryos. We used morpholino antisense oligonucleotides (MO or morpholino) to block initiation of translation of the target genes. TCDD-induced mRNA expression of CYP1Cs and a decrease in MsV blood flow were both blocked by gene knockdown of AHR2. Gene knockdown of CYP1C1 by two different morpholinos and CYP1C2 by two different morpholinos, but not by their 5 nucleotide-mismatch controls, was effective in blocking reduced MsV blood flow caused by TCDD. The same CYP1C-MOs prevented reduction of blood flow in the MsV caused by β-naphthoflavone (BNF), representing another class of AHR agonists. Whole-mount in situ hybridization revealed that mRNA expression of CYP1C1 and CYP1C2 was induced by TCDD most strongly in branchiogenic primordia and pectoral fin buds. In situ hybridization using head transverse sections showed that TCDD increased the expression of both CYP1Cs in endothelial cells of blood vessels, including the MsV. These results indicate a potential role of CYP1C1 and CYP1C2 in the local circulation failure induced by AHR2 activation in the dorsal midbrain of the zebrafish embryo.

Comparative transcriptomics implicates mechanisms of evolved pollution tolerance in a killifish population

Wild populations of the killifish Fundulus heteroclitus resident in heavily contaminated North American Atlantic coast estuaries have recently and independently evolved dramatic, heritable, and adaptive pollution tolerance. We compared physiological and transcriptome responses to embryonic polychlorinated biphenyl (PCB) exposures between one tolerant population and a nearby sensitive population to gain insight into genomic, physiological and biochemical mechanisms of evolved tolerance in killifish, which are currently unknown. The PCB exposure concentrations at which developmental toxicity emerged, the range of developmental abnormalities exhibited, and global as well as specific gene expression patterns were profoundly different between populations. In the sensitive population, PCB exposures produced dramatic, dose-dependent toxic effects, concurrent with the alterations in the expression of many genes. For example, PCB-mediated cardiovascular system failure was associated with the altered expression of cardiomyocyte genes, consistent with sarcomere mis-assembly. In contrast, genome-wide expression was comparatively refractory to PCB induction in the tolerant population. Tolerance was associated with the global blockade of the aryl hydrocarbon receptor (AHR) signalling pathway, the key mediator of PCB toxicity, in contrast to the strong dose-dependent up-regulation of AHR pathway elements observed in the sensitive population. Altered regulation of signalling pathways that cross-talk with AHR was implicated as one candidate mechanism for the adaptive AHR signalling repression and the pollution tolerance that it affords. In addition to revealing mechanisms of PCB toxicity and tolerance, this study demonstrates the value of comparative transcriptomics to explore molecular mechanisms of stress response and evolved adaptive differences among wild populations.

Dynamic zebrafish interactome reveals transcriptional mechanisms of dioxin toxicity

Background

In order to generate hypotheses regarding the mechanisms by which 2,3,7,8-tetrachlorodibenzo-p-dioxin (dioxin) causes toxicity, we analyzed global gene expression changes in developing zebrafish embryos exposed to this potent toxicant in the context of a dynamic gene network. For this purpose, we also computationally inferred a zebrafish (Danio rerio) interactome based on orthologs and interaction data from other eukaryotes.

Methodology/Principal Findings

Using novel computational tools to analyze this interactome, we distinguished between dioxin-dependent and dioxin-independent interactions between proteins, and tracked the temporal propagation of dioxin-dependent transcriptional changes from a few genes that were altered initially, to large groups of biologically coherent genes at later times. The most notable processes altered at later developmental stages were calcium and iron metabolism, embryonic morphogenesis including neuronal and retinal development, a variety of mitochondria-related functions, and generalized stress response (not including induction of antioxidant genes). Within the interactome, many of these responses were connected to cytochrome P4501A (cyp1a) as well as other genes that were dioxin-regulated one day after exposure. This suggests that cyp1a may play a key role initiating the toxic dysregulation of those processes, rather than serving simply as a passive marker of dioxin exposure, as suggested by earlier research.

Conclusions/Significance

Thus, a powerful microarray experiment coupled with a flexible interactome and multi-pronged interactome tools (which are now made publicly available for microarray analysis and related work) suggest the hypothesis that dioxin, best known in fish as a potent cardioteratogen, has many other targets. Many of these types of toxicity have been observed in mammalian species and are potentially caused by alterations to cyp1a.

Dioxin exposure disrupts the differentiation of mouse embryonic stem cells into cardiomyocytes

Experimental exposure of fish, birds, and rodents to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD; dioxin) causes multiple Ah receptor-mediated developmental abnormalities, an observation consistent with compelling evidence in human populations that TCDD exposure is responsible for a significant incidence of birth defects. To characterize molecular mechanisms that might explain the developmental effects of dioxin, we have studied the consequences of TCDD exposure on the differentiation of mouse embryonic stem (ES) cells in culture and on the expression of genes, including those coding for homeodomain containing transcription factors, with a role in progression of tissue differentiation and embryonic identity during development. We find that TCDD treatment causes expression changes in a number of homeobox genes concomitant with Ah receptor recruitment to the promoters of many of these genes, whether under naïve or dioxin-activated conditions. TCDD exposure also derails temporal expression trajectories of developmentally regulated genes in a wide diversity of differentiation pathways, including genes with functions in neural and cardiovascular development, self-renewal, hematopoiesis and mesenchymal lineage specification, and Notch and Wnt pathways. Among these, we find that TCDD represses the expression of the cardiac development-specific Nkx2.5 homeobox transcription factor, of cardiac troponin-T and of alpha- and beta-myosin heavy chains, inhibiting the formation of beating cardiomyocytes, a characteristic phenotype of differentiating mouse ES cells in culture. These data identify potential pathways for dioxin to act as a developmental teratogen, possibly critical to cardiovascular development and disease, and provide molecular targets that may help us understand the molecular basis of Ah receptor-mediated developmental toxicity.

Receptor- and reactive intermediate-mediated mechanisms of teratogenesis

Drugs and environmental chemicals can adversely alter the development of the fetus at critical periods during pregnancy, resulting in death, or in structural and functional birth defects in the surviving offspring. This process of teratogenesis may not be evident until a decade or more after birth. Postnatal functional abnormalities include deficits in brain function, a variety of metabolic diseases, and cancer. Due to the high degree of fetal cellular division and differentiation, and to differences from the adult in many biochemical pathways, the fetus is highly susceptible to teratogens, typically at low exposure levels that do not harm the mother. Insights into the mechanisms of teratogenesis come primarily from animal models and in vitro systems, and involve either receptor-mediated or reactive intermediate-mediated processes. Receptor-mediated mechanisms involving the reversible binding of xenobiotic substrates to a specific receptor are exemplified herein by the interaction of the environmental chemical 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD or "dioxin") with the cytosolic aryl hydrocarbon receptor (AHR), which translocates to the nucleus and, in association with other proteins, binds to AH-responsive elements (AHREs) in numerous genes, initiating changes in gene transcription that can perturb development. Alternatively, many xenobiotics are bioactivated by fetal enzymes like the cytochromes P450 (CYPs) and prostaglandin H synthases (PHSs) to highly unstable electrophilic or free radical reactive intermediates. Electrophilic reactive intermediates can covalently (irreversibly) bind to and alter the function of essential cellular macromolecules (proteins, DNA), causing developmental anomalies. Free radical reactive intermediates can enhance the formation of reactive oxygen species (ROS), resulting in oxidative damage to cellular macromolecules and/or altered signal transduction. The teratogenicity of reactive intermediates is determined to a large extent by the balance among embryonic and fetal pathways of xenobiotic bioactivation, detoxification of the xenobiotic reactive intermediate, detoxification of ROS, and repair of oxidative macromolecular damage.

The zebrafish gill model: induction of CYP1A, EROD and PAH adduct formation

Zebrafish (Danio rerio) is a common model species in fish toxicology, and the zebrafish gill is potentially useful in screening waterborne pollutants. We have previously developed a gill-based ethoxyresorufin-O-deethylase (EROD) assay, and proposed gill EROD activity as a biomarker for exposure to waterborne aryl hydrocarbon receptor (AHR) agonists. In this study we modified the gill EROD assay for use in zebrafish. We used immunohistochemistry to localize CYP1A induction, and microautoradiography to localize irreversible binding of the prototypic polycyclic aromatic hydrocarbon, 7,12-dimethylbenz[a]anthracene (DMBA) in zebrafish gills. Gill filament and liver microsomal EROD activities were measured after waterborne exposure of zebrafish and rainbow trout to benzo[a]pyrene (BaP) or beta-naphthoflavone (betaNF). The results showed considerably lower relative EROD induction by betaNF (1microM) in zebrafish than in rainbow trout, both in gills (13-fold versus 230-fold compared to control) and in liver (5-fold versus 320-fold compared to control). The induced hepatic EROD activity was similar in the two species, whereas the basal activity was considerably higher in zebrafish than in rainbow trout. In zebrafish gills, betaNF enhanced DMBA adduct formation and CYP1A immunostaining. Ellipticine blocked DMBA adduct formation and EROD activity following betaNF exposure but had no effect on CYP1A immunostaining. A notable finding was that the localization of DMBA adducts differed from that of CYP1A protein in betaNF-induced fish; CYP1A immunoreactivity was evenly distributed in the gills whereas DMBA adduction was confined to the leading edges of the filaments and the gill rakers, i.e., structures being highly exposed to DMBA-containing inhaled water. The results show that the modified method is suitable for determination of gill EROD activity in zebrafish, although rainbow trout seems more sensitive. They also imply that the sites of DMBA adduct formation in zebrafish gills are markedly influenced by kinetic factors.

Cytochrome P450 1D1: a novel CYP1A-related gene that is not transcriptionally activated by PCB126 or TCDD

Enzymes in the cytochrome P450 1 family oxidize many common environmental toxicants. We identified a new CYP1, termed CYP1D1, in zebrafish. Phylogenetically, CYP1D1 is paralogous to CYP1A and the two share 45% amino acid identity and similar gene structure. In adult zebrafish, CYP1D1 is most highly expressed in liver and is relatively highly expressed in brain. CYP1D1 transcript levels were higher at 9h post-fertilization than at later developmental times. Treatment of zebrafish with potent aryl hydrocarbon receptor (AHR) agonists (3,3',4,4',5-pentachlorobiphenyl or 2,3,7,8-tetrachlorodibenzo-p-dioxin) did not induce CYP1D1 transcript expression. Morpholino oligonucleotide knockdown of AHR2, which mediates induction of other CYP1s, did not affect CYP1D1 expression. Zebrafish CYP1D1 heterologously expressed in yeast exhibited ethoxyresorufin- and methoxyresorufin-O-dealkylase activities. Antibodies against a CYP1D1 peptide specifically detected a single electrophoretically-resolved protein band in zebrafish liver microsomes, distinct from CYP1A. CYP1D1 in zebrafish is a CYP1A-like gene that could have metabolic functions targeting endogenous compounds.

Protein expression profiling in the African clawed frog Xenopus laevis tadpoles exposed to the polychlorinated biphenyl mixture aroclor 1254

Exposure to environmental pollutants such as polychlorinated biphenyls (PCBs) is now taken into account to partly explain the worldwide decline of amphibians. PCBs induce deleterious effects on developing amphibians including deformities and delays in metamorphosis. However, the molecular mechanisms by which they express their toxicity during the development of tadpoles are still largely unknown. A proteomics analysis was performed on developing Xenopus laevis tadpoles exposed from 2 to 5 days postfertilization to either 0.1 or 1 ppm Aroclor 1254, a PCB mixture. Two-dimensional DIGE with a minimal labeling method coupled to nanoflow liquid chromatography-tandem mass spectrometry was used to detect and identify proteins differentially expressed under PCBs conditions. Results showed that 59 spots from the 0.1 ppm Aroclor 1254 condition and 57 spots from the 1 ppm Aroclor 1254 condition displayed a significant increase or decrease of abundance compared with the control. In total, 28 proteins were identified. The results suggest that PCBs induce mechanisms against oxidative stress (peroxiredoxins 1 and 2), adaptative changes in the energetic metabolism (enolase 1, glycerol-3-phosphate dehydrogenase, and creatine kinase muscle and brain types), and the implication of the unfolded protein response system (glucose-regulated protein, 58 kDa). They also affect, at least at the highest concentration tested, the synthesis of proteins involved in normal cytogenesis (alpha-tropomyosin, myosin heavy chain, and alpha-actin). For the first time, proteins such as aldehyde dehydrogenase 7A1, CArG binding factor-A, prolyl 4-hydroxylase beta, and nuclear matrix protein 200 were also shown to be up-regulated by PCBs in developing amphibians. These data argue that protein expression reorganization should be taken into account while estimating the toxicological hazard of wild amphibian populations exposed to PCBs.