Transcriptomic Changes in Zebrafish Embryos and Larvae Following Benzo[a]pyrene Exposure

Intervention mechanism of marine-based chito-oligosaccharide on acute liver injury induced by AFB1 in rats

Aflatoxin B1 (AFB1) is extremely hepatotoxic, a causative agent of liver cancer, and can cause symptoms of acute or chronic liver damage. Chito-oligosaccharides (COS), obtained from the degradation of chitosan derived from shrimp and crab shells, is a natural antioxidant substance and its antitumor properties have been widely studied, but less research has been done on the prevention of AFB1-induced acute liver injury. In this study, rats were acutely exposed to 1 mg/kg BW AFB1 and simultaneously gavaged with different doses of COS for 8 days. The results showed that COS attenuated the hepatic histopathological changes and reduced serum biochemical indices (ALT, AST, ALP, and TBIL) in rats. It significantly inhibited MDA content and promoted SOD and GSH-Px activity production. Moreover, it also improved hepatocyte apoptosis. Furthermore, AFB1-vs-HCOS differential genes were enriched with 622 GO entries, and 380 were Biological Processes, 170 were Molecular Functions, 72 were Cellular Components. Differentially expressed genes (DEGs) analyzed by KEGG enrichment were more enriched in pathways, such as metabolism, PPAR signaling pathway, and peroxisome. Q-PCR technique verified that Lama5, Egr1, Cyp2b1, and Gadd45g in DEGs were associated with oxidative stress damage and apoptosis. In conclusion, COS intervention reduces the effect of AFB1 on hepatic genes and thus reduces the changes in hepatic gene function.

A complex mixture of polycyclic aromatic compounds causes embryotoxic, behavioral, and molecular effects in zebrafish larvae (Danio rerio), and in vitro bioassays

Polycyclic aromatic compounds (PACs) are prevalent in the environment, typically found in complex mixtures and high concentrations. Our understanding of the effects of PACs, excluding the 16 priority polycyclic aromatic hydrocarbons (16 PAHs), remains limited. Zebrafish embryos and in vitro bioassays were utilized to investigate the embryotoxic, behavioral, and molecular effects of a soil sample from a former gasworks site in Sweden. Additionally, targeted chemical analysis was conducted to analyze 87 PACs in the soil, fish, water, and plate material. CALUX® assays were used to assess the activation of aryl hydrocarbon and estrogen receptors, as well as the inhibition of the androgen receptor. Larval behavior was measured by analyzing activity during light and darkness and in response to mechanical stimulation. Furthermore, qPCR analyses were performed on a subset of 36 genes associated with specific adverse outcomes, and the total lipid content in the larvae was measured. Exposure to the sample resulted in embryotoxic effects (LC50 = 0.480 mg dry matter soil/mL water). The mixture also induced hyperactivity in darkness and hypoactivity in light and in response to the mechanical stimulus. qPCR analysis revealed differential regulation of 15 genes, including downregulation of opn1sw1 (eye pigmentation) and upregulation of fpgs (heart failure). The sample caused significant responses in three bioassays (ERα-, DR-, and PAH-CALUX), and the exposed larvae exhibited elevated lipid levels. Chemical analysis identified benzo[a]pyrene as the predominant compound in the soil and approximately half of the total PAC concentration was attributed to the 16 PAHs. This study highlights the value of combining in vitro and in vivo methods with chemical analysis to assess toxic mechanisms at specific targets and to elucidate the possible interactions between various pathways in an organism. It also enhances our understanding of the risks associated with environmental mixtures of PACs and their distribution during toxicity testing.

Alkylation of Benz[a]anthracene Affects Toxicity to Early-Life Stage Zebrafish and In Vitro Aryl Hydrocarbon Receptor 2 Transactivation in a Position-Dependent Manner

Polycyclic aromatic hydrocarbons (PAHs) are structurally diverse organic chemicals that can have adverse effects on the health of fishes through activation of aryl hydrocarbon receptor 2 (AhR2). They are ubiquitous in the environment, but alkyl PAHs are more abundant in some environmental matrices. However, relatively little is known regarding the effects of alkylation on the toxicity of PAHs to fishes in vivo and how this relates to potency for activation of AhR2 in vitro. Therefore, the objectives of the present study were to determine the toxicity of benz[a]anthracene and three alkylated homologs representing various alkylation positions to early life stages of zebrafish (Danio rerio) and to assess the potency of each for activation of the zebrafish AhR2 in a standardized in vitro AhR transactivation assay. Exposure of embryos to each of the PAHs caused a dose-dependent increase in mortality and malformations characteristic of AhR2 activation. Each alkyl homolog had in vivo toxicities and in vitro AhR2 activation potencies different from those of the parent PAH in a position-dependent manner. However, there was no statistically significant linear relationship between responses measured in these assays. The results suggest a need for further investigation into the effect of alkylation on the toxicity of PAHs to fishes and greater consideration of the contribution of alkylated homologs in ecological risk assessments. Environ Toxicol Chem 2022;41:1993-2002. © 2022 SETAC.

Transcriptome aberration associated with altered locomotor behavior of zebrafish (Danio rerio) caused by Waterborne Benzo[a]pyrene

Waterborne Benzo[a]pyrene (B[a]P) pollution is a global threat to aquatic organisms. The exposure to waterborne B[a]P can disrupt the normal locomotor behavior of zebrafish (Danio rerio), however, how it affect the locomotor behavior of adult zebrafish remains unclear. Herein, B[a]P at two concentrations (0.8 μg/L and 2.0 μg/L) were selected to investigate the molecular mechanisms of the affected locomotor behavior of zebrafish by B[a]P based on transcriptome profiling. Adverse effects of B[a]P exposure affecting locomotor behavior in zebrafish were studied by RNA sequencing, and the locomotion phenotype was acquired. The gene enrichment results showed that the differentially highly expressed genes (atp2a1, cdh2, aurka, fxyd1, clstn1, apoc1, mt-co1, tnnt3b, and fads2) of zebrafish are mainly enriched in adrenergic signaling in cardiomyocytes (dre04261) and locomotory behavior (GO:0007626). The movement trajectory plots showed an increase in the locomotor distance and velocity of zebrafish in the 0.8 μg/L group and the opposite in the 2.0 μg/L group. The results showed that B[a]P affects the variety of genes in zebrafish, including motor nerves, muscles, and energy supply, and ultimately leads to altered locomotor behavior.

The effect of benzo[a]pyrene on the gut microbiota of Nile tilapia (Oreochromis niloticus)

Benzo[a]pyrene (BaP) is a highly toxic and carcinogenic polycyclic aromatic hydrocarbon (PAH) whose toxicological effects in the gut microbiota of aquatic organisms have not yet been fully revealed. Therefore, in this study, we used high-throughput 16S rRNA gene sequencing to evaluate the effects of BaP in the gut microbiome of Oreochromis niloticus, including its possible participation in the process of detoxification and its ability to recover. The fish were injected with a single intraperitoneal dose of 20 mg kg^-1 of BaP, and the effects in the microbiome were evaluated at 24, 72, and 120 h post-injection. The results indicate a clear dysbiosis (in composition, relative abundance, diversity, and interaction networks) of the gut microbiota during 24 h post-injection, dominated by Fusobacteria and Bacteroidetes and a decrease in Proteobacteria and Spirochaetae. Interestingly, a slight recovery of the microbiome begins at 72 h and stabilises at 120 h post-injection. Pathway analysis revealed the participation of the gut microbiome in PAH degradation mainly at 24 h post-injection. This study provides new insights in the toxicology of BaP in O. niloticus and the first evidence of the ability of the gut microbiome to recovery after a chemical disturbance. KEY POINTS: • Benzo[a]pyrene caused a dysbiosis in the gut microbiota of Oreochromis niloticus. • We observed an enrichment of bacteria involved in the metabolism of xenobiotics. • The gut microbiota was recovered after exposure to benzo[a]pyrene.

Knockdown of DNA methyltransferase 1 reduces DNA methylation and alters expression patterns of cardiac genes in embryonic cardiomyocytes

We previously found that DNA methyltransferase 3a (DNMT3a) plays an important role in regulating embryonic cardiomyocyte gene expression, morphology, and function. In this study, we investigated the role of the most abundant DNMT in mammalian cells, DNMT1, in these processes. It is known that DNMT1 is essential for embryonic development, during which it is involved in regulating cardiomyocyte DNA methylation and gene expression. We used siRNA to knock down DNMT1 expression in primary cultures of mouse embryonic cardiomyocytes. Immunofluorescence staining and multielectrode array were, respectively, utilized to evaluate cardiomyocyte growth and electrophysiology. RNA sequencing (RNA‐Seq) and multiplex bisulfite sequencing were, respectively, performed to examine gene expression and promoter methylation. At 72 h post‐transfection, reduction of DNMT1 expression decreased the number and increased the size of embryonic cardiomyocytes. Beat frequency and the amplitude of field action potentials were decreased by DNMT1 siRNA. RNA‐Seq analysis identified 801 up‐regulated genes and 494 down‐regulated genes in the DNMT1 knockdown cells when compared to controls. Pathway analysis of the differentially expressed genes revealed pathways that were associated with cell death and survival, cell morphology, cardiac function, and cardiac disease. Alternative splicing analysis identified 929 differentially expressed exons, including 583 up‐regulated exons and 308 down‐regulated exons. Moreover, decreased methylation levels were found in the promoters of cardiac genes Myh6, Myh7, Myh7b, Tnnc1, Tnni3, Tnnt2, Nppa, Nppb, mef2c, mef2d, Camta2, Cdkn1A, and Cdkn1C. Of these 13 genes, 6 (Myh6, Tnnc1, Tnni3, Tnnt2, Nppa, Nppb) and 1 (Cdkn1C) had increased or decreased gene expression, respectively. Altogether, these data show that DNMT1 is important in embryonic cardiomyocytes by regulating DNA methylation, gene expression, gene splicing, and cell function.

In vivo assessment of respiratory burst inhibition by xenobiotic exposure using larval zebrafish

Currently, assessment of the potential immunotoxicity of a given agent involves a tiered approach for hazard identification and mechanistic studies, including observational studies, evaluation of immune function, and measurement of susceptibility to infectious and neoplastic diseases. These studies generally use costly low-throughput mammalian models. Zebrafish, however, offer an excellent alternative due to their rapid development, ease of maintenance, and homology to mammalian immune system function and development. Larval zebrafish also are a convenient model to study the innate immune system with no interference from the adaptive immune system. In this study, a respiratory burst assay (RBA) was utilized to measure reactive oxygen species (ROS) production after developmental xenobiotic exposure. Embryos were exposed to non-teratogenic doses of chemicals and at 96 h post-fertilization, the ability to produce ROS was measured. Using the RBA, 12 compounds with varying immune-suppressive properties were screened. Seven compounds neither suppressed nor enhanced the respiratory burst; five reproducibly suppressed global ROS production, but with varying potencies: benzo[a]pyrene, 17β-estradiol, lead acetate, methoxychlor, and phenanthrene. These five compounds have all previously been reported as immunosuppressive in mammalian innate immunity assays. To evaluate whether the suppression of ROS by these compounds was a result of decreased immune cell numbers, flow cytometry with transgenic zebrafish larvae was used to count the numbers of neutrophils and macrophages after chemical exposure. With this assay, benzo[a]pyrene was found to be the only chemical that induced a change in the number of immune cells by increasing macrophage but not neutrophil numbers. Taken together, this work demonstrates the utility of zebrafish larvae as a vertebrate model for identifying compounds that impact innate immune function at non-teratogenic levels and validates measuring ROS production and phagocyte numbers as metrics for monitoring how xenobiotic exposure alters the innate immune system.

The chemical defensome of five model teleost fish

How an organism copes with chemicals is largely determined by the genes and proteins that collectively function to defend against, detoxify and eliminate chemical stressors. This integrative network includes receptors and transcription factors, biotransformation enzymes, transporters, antioxidants, and metal- and heat-responsive genes, and is collectively known as the chemical defensome. Teleost fish is the largest group of vertebrate species and can provide valuable insights into the evolution and functional diversity of defensome genes. We have previously shown that the xenosensing pregnane x receptor (pxr, nr1i2) is lost in many teleost species, including Atlantic cod (Gadus morhua) and three-spined stickleback (Gasterosteus aculeatus), but it is not known if compensatory mechanisms or signaling pathways have evolved in its absence. In this study, we compared the genes comprising the chemical defensome of five fish species that span the teleosteii evolutionary branch often used as model species in toxicological studies and environmental monitoring programs: zebrafish (Danio rerio), medaka (Oryzias latipes), Atlantic killifish (Fundulus heteroclitus), Atlantic cod, and three-spined stickleback. Genome mining revealed evolved differences in the number and composition of defensome genes that can have implication for how these species sense and respond to environmental pollutants, but we did not observe any candidates of compensatory mechanisms or pathways in cod and stickleback in the absence of pxr. The results indicate that knowledge regarding the diversity and function of the defensome will be important for toxicological testing and risk assessment studies.

Gene coexpression networks reveal molecular interactions underlying cichlid jaw modularity

Background

The oral and pharyngeal jaw of cichlid fishes are a classic example of evolutionary modularity as their functional decoupling boosted trophic diversification and contributed to the success of cichlid adaptive radiations. Most studies until now have focused on the functional, morphological, or genetic aspects of cichlid jaw modularity. Here we extend this concept to include transcriptional modularity by sequencing whole transcriptomes of the two jaws and comparing their gene coexpression networks.

Results

We show that transcriptional decoupling of gene expression underlies the functional decoupling of cichlid oral and pharyngeal jaw apparatus and the two units are evolving independently in recently diverged cichlid species from Lake Tanganyika. Oral and pharyngeal jaw coexpression networks reflect the common origin of the jaw regulatory program as there is high preservation of gene coexpression modules between the two sets of jaws. However, there is substantial rewiring of genetic architecture within those modules. We define a global jaw coexpression network and highlight jaw-specific and species-specific modules within it. Furthermore, we annotate a comprehensive in silico gene regulatory network linking the Wnt and AHR signalling pathways to jaw morphogenesis and response to environmental cues, respectively. Components of these pathways are significantly differentially expressed between the oral and pharyngeal jaw apparatus.

Conclusion

This study describes the concerted expression of many genes in cichlid oral and pharyngeal jaw apparatus at the onset of the independent life of cichlid fishes. Our findings suggest that – on the basis of an ancestral gill arch network—transcriptional rewiring may have driven the modular evolution of the oral and pharyngeal jaws, highlighting the evolutionary significance of gene network reuse. The gene coexpression and in silico regulatory networks presented here are intended as resource for future studies on the genetics of vertebrate jaw morphogenesis and trophic adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01787-9.

Associations between weekly air pollution exposure and congenital heart disease

BACKGROUND

The topic of congenital heart diseases (CHDs) has attracted more and more attentions. Accumulating evidence suggests that exposure to air pollutants during pregnancy is associated with CHDs, yet the results are inconsistent and study about weekly exposure is few. Our study evaluated the association between weekly air pollution and CHDs in Hefei, China.

MATERIALS AND METHODS

Daily CHDs admission data were obtained from eight hospitals in Hefei from October 2015 to September 2017. Meteorological data and air quality were collected from China Meteorological Data Network. Distributed lag nonlinear model (DLNM) considering both the lag effect of exposure factors and the nonlinear relationship of exposure-reaction was used to assess the effect of weekly air pollutants exposure on CHDs admission.

RESULTS

During the study period, totally 47,046 cases of perinatal infants were recruited, and the incidence of CHDs was 9.71 per thousand. The findings showed PM_2.5, PM₁₀, SO₂ and NO₂ significantly increased the risk of CHDs. Each 10 μg/m³ increase in PM_2.5 during gestational weeks 20-26 increased the risk of CHDs. The susceptibility windows of PM₁₀ (weeks 0-2 and weeks 25-29 of pregnancy), SO₂ (weeks 8-16 and weeks 29-38) and NO₂ (week 40), while the strongest effects of these 4 pollutants on CHDs were observed in week 22 (RR = 1.034, 95% CI: 1.007-1.062), week 0 (RR = 1.081, 95% CI: 1.02-1.146), week 37 (RR = 1.528, 95% CI: 1.085-2.153) and week 40 (RR = 1.171, 95% CI: 1.006-1.364), respectively.

CONCLUSIONS

Air pollutants (SO₂, NO₂, PM₁₀, and PM_2.5) exposure could increase the risk of CHDs, while the most crucial susceptibility windows for the exposure were mainly in the second and third trimesters. Boys seemed to be more sensitive to air pollution. Our study contributes to the knowledge of the association between maternal exposure to air pollution and CHDs, but the associations need to be verified by further studies.

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.

Acute exposure to oil induces age and species-specific transcriptional responses in embryo-larval estuarine fish

Because oil spills frequently occur in coastal regions that serve as spawning habitat, characterizing the effects of oil in estuarine fish carries both economic and environmental importance. There is a breadth of research investigating the effects of crude oil on fish, however few studies have addressed how transcriptional responses to oil change throughout development or how these responses might be conserved across taxa. To investigate these effects, we performed RNA-seq and pathway analysis following oil exposure 1) in a single estuarine species (Cyprinodon variegatus) at three developmental time points (embryos, yolk-sack larvae, free-feeding larvae), and 2) in two ecologically similar species (C. variegatus and Fundulus grandis), immediately post-hatch (yolk-sack stage). Our results indicate that C. variegatus embryos mount a diminished transcriptional response to oil compared to later stages, and that few transcriptional responses are conserved throughout development. Pathway analysis of larval C. variegatus revealed dysregulation of similar biological processes at later larval stages, including alteration of cholesterol biosynthesis pathways, cardiac development processes, and immune functions. Our cross-species comparison showed that F. grandis exhibited a reduced transcriptional response compared to C. variegatus. Pathway analysis revealed that the two species shared similar immune and cardiac responses, however pathways related to cholesterol biosynthesis exhibited a divergent response as they were activated in C. variegatus but inhibited in F. grandis. Our results suggest that examination of larval stages may provide a more sensitive estimate of oil-impacts than examination of embryos, and challenge assumptions that ecologically comparable species respond to oil similarly.

Cytochrome P450 1A mRNA in the Gambusia affinis and Response to Several PAHs

Cytochrome P4501A (CYP1A) has been used as a specific biomarker for monitoring water contamination such as PAHs, PCBs and dioxins. In the present study, the cyp1a gene of Gambusia affinis was cloned and sequenced and their expressions under PAHs exposure were characterized. The newly identified cyp1a encodes a protein with 521 amino acids that shared 96-80% identity with other Cyprinodontiformes fishes. RT-PCR analysis revealed that the basal mRNA level of cyp1a was highly expressed in liver and intestine. The expression level of cyp1a was significantly induced by exposure to 100 μg/L 3, 4-Benzopyrene (BaP) for 5 days in the muscle, testis, brain, liver and intestine of adult male fish. Except in the testis, the induced mRNA level of cyp1a ultimately decreased after prolonging the exposure time to 25 days. As for testis, the induced mRNA level of cyp1a was maintained at a high level during the entire exposure time under 100 μg/L BaP exposure. Furthermore, the expression of cyp1a increased with exposure time under a relatively low exposure concentrations 1 μg/L. Regarding the effects of other PAHs, D(a,h)A, BbF, and BaA showed a statistically significant effect of induction on mRNA level of cyp1a in the muscle, testis, brain, liver and intestine.

Transcriptomic changes in western mosquitofish (Gambusia affinis) liver following benzo[a]pyrene exposure

Widely distributed western mosquitofish (Gambusia affinis) has been used as a new model species for hazard assessment of environmental stressors such as polycyclic aromatic hydrocarbons (PAHs). However, most of the PAH studies using G. affinis rely on targeted biomarker-based analysis, and thus may not adequately address the complexity of the toxic mechanisms of the stressors. In the present study, the whole transcriptional sequencing of G. affinis liver after exposure to a PAH model, benzo[a]pyrene (BaP) (100 μg/L), for 20 days was performed by using the HiSeq XTen sequencers. In total, 58,156,233 and 51,825,467 clean nucleotide reads were obtained in the control and BaP-exposed libraries, respectively, with average N50 lengths of 1419 bp. In addition, after G. affinis was exposed for 20 days, 169 genes were upregulated, and 176 genes were downregulated in liver. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were applied to all the genes to determine the genes' biological functions and processes. The results clearly showed that the differentially expressed genes were mainly related to immune pathways and metabolic correlation pathways. Interestingly, almost all the pathways related with the immunity were upregulated, while the metabolism pathways were downregulated. Lastly, quantitative real-time PCR (qRT-PCR) was performed to measure expressional levels of twelve genes confirmed through the DGE analysis. These results demonstrate that BaP damages immunity and enhances the consumption of all available energy storage to activate mechanisms of the detoxification in G. affinis. Up until now, the present study is the first time that a whole transcriptome sequencing analysis in the liver of G. affinis exposed to BaP has been reported.

Parental exposure to bisphenol A affects pharyngeal cartilage development and causes global transcriptomic changes in zebrafish (Danio rerio) offspring

BACKGROUND

As one of the most common endocrine-disrupting chemicals (EDCs), bisphenol A (BPA) is a threat to aquatic ecosystems. Despite a rich literature addressing the adverse effects of BPA on various systems in fish models, the potential impact of parental BPA exposure on offspring pharyngeal cartilage development is poorly understood.

METHODS

Adult zebrafish (F0) were exposed to BPA (1.0 μM) or control for 7 days. Eggs (F1) were collected and exposed to BPA (control, 0.05, 0.1, 1, 10 μM) until 120 h post-fertilization. Histomorphometrical essay was used to quantitatively and qualitatively assess the effects of BPA on pharyngeal cartilage development. RNA sequencing (RNA-seq) was used to discover differentially expressed genes (DEGs), and KEGG pathway and GO enrichment analysis were performed to interpret functional ontology.

RESULTS

Parental BPA exposure affected hatchability and heart rates of F1 progeny. By pathology analysis, parental BPA exposure caused craniofacial deformity, characterized by wider angles of cartilage elements, disrupted pharyngeal chondrocytes and promoted apoptosis and elongation of head length. RNA-seq suggested that many DEGs were involved in multiple biological processes and signaling pathways; defense responses, reactive oxygen species metabolic process, apoptosis, p53 signaling pathway and MAPK signaling pathway were closely associated with the toxicity of parental BPA exposure.

CONCLUSIONS

Parental BPA exposure affected chondrogenesis in the viscerocranium of zebrafish offspring and led to global transcriptomic changes involved in apoptosis, hyperplasia and oxidative stress. These newly identified gene expression patterns, pathways and gene networks of zebrafish eleutheroembryos after early-life waterborne BPA exposure, may lead to severe and permanent morphological and functional consequences.

Benzo[a]pyrene and 2,3-benzofuran induce divergent temporal patterns of AhR-regulated responses in zebrafish embryos (Danio rerio)

Biotests like the fish embryo toxicity test have become increasingly popular in risk assessment and evaluation of chemicals found in the environment. The large range of possible endpoints is a big advantage when researching on the mode of action of a certain substance. Here, we utilized the frequently used model organism zebrafish (Danio rerio) to examine regulative mechanisms in the pathway of the aryl-hydrocarbon receptor (AHR) in early development. We exposed embryos to representatives of two chemical classes known to elicit dioxin-like activity: benzo[a]pyrene for polycyclic aromatic hydrocarbons (PAHs) and 2,3-benzofuran for polar O-substituted heterocycles as a member of heterocyclic compounds in general (N-, S-, O-heterocycles; NSO-hets). We measured gene transcription of the induced P450 cytochromes (cyp1), their formation of protein and biotransformation activity throughout the whole embryonic development until 5 days after fertilization. The results show a very specific time course of transcription depending on the chemical properties (e.g. halogenation, planarity, K_ow), the physical decay and the biodegradability of the tested compound. However, although this temporal pattern was not precisely transferable onto the protein level, significant regulation in enzymatic activity over time could be detected. We conclude, that a careful choice of time and end point as well as consideration of the chemical properties of a substance are fairly important when planning, conducting and especially evaluating biotests.

Development of a Zebrafish S1500+ Sentinel Gene Set for High-Throughput Transcriptomics

Sentinel gene sets have been developed with the purpose of maximizing the information from targeted transcriptomic platforms. We recently described the development of an S1500+ sentinel gene set, which was built for the human transcriptome, utilizing a data- and knowledge-driven hybrid approach to select a small subset of genes that optimally capture transcriptional diversity, correlation with other genes based on large-scale expression profiling, and known pathway annotation within the human genome. While this detailed bioinformatics approach for gene selection can in principle be applied to other species, the reliability of the resulting gene set depends on availability of a large body of transcriptomics data. For the model organism zebrafish, we aimed to create a similar sentinel gene set (Zf S1500+ gene set); however, there is insufficient standardized expression data in the public domain to train the gene correlation model. Therefore, our strategy was to use human-zebrafish ortholog mapping of the human S1500+ genes and nominations from experts in the zebrafish scientific community. In this study, we present the bioinformatics curation and refinement process to produce the final Zf S1500+ gene set, explore whole transcriptome extrapolation using this gene set, and assess pathway-level inference. This gene set will add value to targeted high-throughput transcriptomics in zebrafish for toxicogenomic screening and other research domains.

Transcriptomic analysis reveals dose-dependent modes of action of benzo(a)pyrene in polar cod (Boreogadus saida)

Polar cod (Boreogadus saida) has been used as a model Arctic species for hazard assessment of environmental stressors such as polycyclic aromatic hydrocarbons (PAHs). However, most of the PAH studies using polar cod rely on targeted biomarker-based analysis thus may not adequately address the complexity of the toxic mechanisms of the stressors. The present study was performed to develop a broad-content transcriptomic platform for polar cod and apply it for understanding the toxic mechanisms of a model PAH, benzo(a)pyrene (BaP). Hepatic transcriptional analysis using a combination of high-density polar cod oligonucleotide microarray and quantitative real-time RT-PCR was conducted to characterize the stress responses in polar cod after 14d repeated dietary exposure to 0.4 (Low) and 20.3 μg/g fish/feeding (High) BaP doses. Bile metabolic analysis was performed to identify the storage of a key BaP hepatic biotransformation product, 3-hydroxybenzo(a)pyrene (3-OH-BaP). The results clearly showed that 3-OH-BaP was detected in the bile of polar cod after both Low and High BaP exposure. Dose-dependent hepatic stress responses were identified, with Low BaP suppressing genes involved in the defense mechanisms and High BaP inducing genes associated with these pathways. The results suggested that activation of the aryl hydrocarbon receptor signaling, induction of oxidative stress, DNA damage and apoptosis were the common modes of action (MoA) of BaP between polar cod or other vertebrates, whereas induction of protein degradation and disturbance of mitochondrial functions were proposed as novel MoAs. Furthermore, conceptual toxicity pathways were proposed for BaP-mediated effects in Arctic fish. The present study has for the first time reported a transcriptome-wide analysis using a polar cod-specific microarray and suggested novel MoAs of BaP. The analytical tools, bioinformatics solutions and mechanistic knowledge generated by this study may facilitate mechanistically-based hazard assessment of environmental stressors in the Arctic using this important fish as a model species.

RNA-Seq analysis of transcriptome responses in Atlantic cod (Gadus morhua) precision-cut liver slices exposed to benzo[a]pyrene and 17α-ethynylestradiol

Polycyclic aromatic hydrocarbons such as benzo[a]pyrene (BaP) that activate the aryl hydrocarbon receptor (Ahr) pathway, and endocrine disruptors acting through the estrogen receptor pathway are among environmental pollutants of major concern. In this work, we exposed Atlantic cod (Gadus morhua) precision-cut liver slices (PCLS) to BaP (10 nM and 1000 nM), ethynylestradiol (EE2) (10 nM and 1000 nM), and equimolar mixtures of BaP and EE2 (10 nM and 1000 nM) for 48 h, and performed RNA-Seq based transcriptome mapping followed by systematic bioinformatics analyses. Our gene expression analysis showed that several genes were differentially expressed in response to BaP and EE2 treatments in PCLS. Strong up-regulation of genes coding for the cytochrome P450 1a (Cyp1a) enzyme and the Ahr repressor (Ahrrb) was observed in BaP treated PCLS. EE2 treatment of liver slices strongly up-regulated genes coding for precursors of vitellogenin (Vtg) and eggshell zona pellucida (Zp) proteins. As expected, pathway enrichment and network analysis showed that the Ahr and estrogen receptor pathways are among the top affected by BaP and EE2 treatments, respectively. Interestingly, two genes coding for fibroblast growth factor 3 (Fgf3) and fibroblast growth factor 4 (Fgf4) were up-regulated by EE2 in this study. To our knowledge, the fgf3 and fgf4 genes have not previously been described in relation to estrogen signaling in fish liver, and these results suggest the modulation of the FGF signaling pathway by estrogens in fish. The signature expression profiles of top differentially expressed genes in response to the single compound (BaP or EE2) treatment were generally maintained in the expression responses to the equimolar binary mixtures. However, in the mixture-treated groups, BaP appeared to have anti-estrogenic effects as observed by lower number of differentially expressed putative EE2 responsive genes. Our in-depth quantitative analysis of changes in liver transcriptome in response to BaP and EE2, using PCLS tissue culture provides further mechanistic insights into effects of the compounds. Moreover, the analyses demonstrate the usefulness of PCLS in cod for omics experiments.

Transcriptome analysis reveals novel insights into the response of low-dose benzo(a)pyrene exposure in male tilapia

Despite a wide number of toxicological studies that describe benzo[a]pyrene (BaP) effects, the metabolic mechanisms that underlie these effects in fish are largely unknown. Of great concern is the presence of BaP in aquatic systems, especially those in close proximity to human activity leading to consumption of potentially contaminated foods. BaP is a known carcinogen and it has been reported to have adverse effects on the survival, development and reproduction of fish. The purpose of this study was to investigate if a low dose of BaP can alter genes and key metabolic pathways in the liver and testis in male adult tilapia, and whether these could be associated with biological endpoints disruption. We used both high-throughput RNA-Sequencing to assess whole genome gene expression following repeated intraperitoneal injections of 3 mg/kg of BaP (every 6 days for 26 days) and morphometric endpoints as indicators of general health. Condition factor (K) along with hepatosomatic and gonadosomatic indices (morphometric parameters) were significantly lower in BaP-treated fish than in controls. BaP exposure induced important changes in the gene expression pattern in liver and testis as revealed by both Pathway and Gene Ontology (GO) analyses. Alterations that were shared by both tissues included arachidonic acid metabolism, androgen receptor to prostate-specific antigen signaling, and insulin-associated effects on lipogenesis. The most salient liver-specific effects included: biological processes involved in detoxification, IL6-associated insulin resistance, mTOR hyperactivation, mitotic cytokinesis, spindle pole and microtubule binding. BaP effects that were confined to the testis included: immune system functions, inflammatory response, estrogen and androgen metabolic pathways. Taken together, gene expression and morphometric end point data indicate that the reproductive success of adult male tilapia could be compromised as a result of BaP exposure. These results constitute new insights on the mechanism of action of low dose BaP in a non-model organism (tilapia).

Omics in Zebrafish Teratogenesis

The genome revolution represents a complete change on our view of biological systems. The quantitative determination of changes in all major molecular components of the living cells, the "omics" approach, opened whole new fields for all health sciences. Genomics, transcriptomics, proteomics, metabolomics, and others, together with appropriate prediction and modeling tools, will mark the future of developmental toxicity assessment both for wildlife and humans. This is especially true for disciplines, like teratology, which rely on studies in model organisms, as studies at lower levels of organization are difficult to implement. Rodents and frogs have been the favorite models for studying human reproductive and developmental disorders for decades. Recently, the study of the development of zebrafish embryos (ZE) is becoming a major alternative tool to adult animal testing. ZE intrinsic characteristics makes this model a unique system to analyze in vivo developmental alterations that only can be studied applying in toto approaches. Moreover, under actual legislations, ZE is considered as a replacement model (and therefore, excluded from animal welfare regulations) during the first 5 days after fertilization. Here we review the most important components of the zebrafish toolbox available for analyzing early stages of embryotoxic events that could eventually lead to teratogenesis.

Functional roles of eriocalyxin B in zebrafish revealed by transcriptome analysis

Eriocalyxin B (EriB) is a naturalent-kaurane diterpenoid obtained fromIsodon eriocalyxvar.laxiflora(family Lamiaceae), which exerted multiple biological activities (e.g.anti-tumor and anti-inflammatory)viathe alteration of gene expression and signaling transduction.

Benzo[a]pyrene Induces Autophagic and Pyroptotic Death Simultaneously in HL-7702 Human Normal Liver Cells

As a common polycyclic aromatic hydrocarbon compound, benzo[a]pyrene (BaP) is readily produced in processing of oil and fatty foods. It is not only a strong carcinogen but also a substance with strong immunotoxicity and reproduction toxicity. Autophagy and pyroptosis are two types of programmed cell death. Whether or not BaP damages body tissues via autophagy or pyroptosis remains unknown. The present study investigated the effects of BaP on autophagy and pyroptosis in HL-7702 cells. The results showed that BaP induced cell death in HL-7702 cells enhanced the intracellular levels of ROS and arrested the cell cycle at the S phase. Additionally, BaP resulted in cell death through autophagy and pyroptosis. Compared with the BaP group, the autophagy inhibitor 3-MA significantly (p < 0.01) inhibited the release of LDH by 70.53% ± 0.46 and NO by 50.36% ± 0.80, the increase of electrical conductivity by 12.08% ± 0.55, and the expressions of pyroptotic marker proteins (Caspase-1, Cox-2, IL-1β, IL-18). The pyroptosis inhibitor Ac-YVAD-CM also notably (p < 0.01) blocked BaP-induced autophagic cell death characterized by the increase of autophagic vacuoles and overexpression of Beclin-1 and LC3-II. In conclusion, BaP led to injury by inducing autophagy and pyroptosis simultaneously, the two of which coexisted and promoted each other in HL-7702 cells.

The single and joint toxicity effects of chlorpyrifos and beta-cypermethrin in zebrafish (Danio rerio) early life stages

Water environment pollution caused by the widespread application of beta-cypermethrin (BCP) and chlorpyrifos (CPF) in agriculture has attracted extensive concern of the world. In this study, zebrafish was used as a model to investigate the individual and joint toxicity of BCP and CPF. In the acute toxicity test, 3 hpf embryos were exposed to various concentrations of CPF, BCP and their binary mixtures (MIX) for 96h. The results indicated that these two pesticides and mixtures induced malformation and death in larvae, and affected hatchability. These two pesticides in mixtures were verified to act together in a synergistic manner under experimental conditions. Oxidative stress assaying manifested that CPF, BCP and MIX altered CAT, SOD and GST activities and MDA content, resulting in oxidative damage in larvae. By pathology analysis, CPF (236μg/L), BCP (5.9μg/L) and MIX (236μg/L CPF+5.9μg/L BCP) were found to trigger liver lesions and promote apoptosis in tissues. The transcriptome sequencing suggested that ECM- receptor interaction, focal adhesion, cell cycle, DNA replication, phototransduction and adherens junction pathways were closely associated with the toxicity of these two pesticides.

Molecular Mechanisms of Crude Oil Developmental Toxicity in Fish

With major oil spills in Korea, the United States, and China in the past decade, there has been a dramatic increase in the number of studies characterizing the developmental toxicity of crude oil and its associated polycyclic aromatic compounds (PACs). The use of model fish species with associated tools for genetic manipulation, combined with high throughput genomics techniques in nonmodel fish species, has led to significant advances in understanding the cellular and molecular bases of functional and morphological defects arising from embryonic exposure to crude oil. Following from the identification of the developing heart as the primary target of crude oil developmental toxicity, studies on individual PACs have revealed a diversity of cardiotoxic mechanisms. For some PACs that are strong agonists of the aryl hydrocarbon receptor (AHR), defects in heart development arise in an AHR-dependent manner, which has been shown for potent organochlorine agonists, such as dioxins. However, crude oil contains a much larger fraction of compounds that have been found to interfere directly with cardiomyocyte physiology in an AHR-independent manner. By comparing the cellular and molecular responses to AHR-independent and AHR-dependent toxicity, this review focuses on new insights into heart-specific pathways underlying both acute and secondary adverse outcomes to crude oil exposure during fish development.

Transgenerational inheritance of neurobehavioral and physiological deficits from developmental exposure to benzo[a]pyrene in zebrafish

Benzo[a]pyrene (B[a]P) is a well-known genotoxic polycylic aromatic compound whose toxicity is dependent on signaling via the aryl hydrocarbon receptor (AHR). It is unclear to what extent detrimental effects of B[a]P exposures might impact future generations and whether transgenerational effects might be AHR-dependent. This study examined the effects of developmental B[a]P exposure on 3 generations of zebrafish. Zebrafish embryos were exposed from 6 to 120h post fertilization (hpf) to 5 and 10μM B[a]P and raised in chemical-free water until adulthood (F0). Two generations were raised from F0 fish to evaluate transgenerational inheritance. Morphological, physiological and neurobehavioral parameters were measured at two life stages. Juveniles of the F0 and F2 exhibited hyper locomotor activity, decreased heartbeat and mitochondrial function. B[a]P exposure during development resulted in decreased global DNA methylation levels and generally reduced expression of DNA methyltransferases in wild type zebrafish, with the latter effect largely reversed in an AHR2-null background. Adults from the F0 B[a]P exposed lineage displayed social anxiety-like behavior. Adults in the F2 transgeneration manifested gender-specific increased body mass index (BMI), increased oxygen consumption and hyper-avoidance behavior. Exposure to benzo[a]pyrene during development resulted in transgenerational inheritance of neurobehavioral and physiological deficiencies. Indirect evidence suggested the potential for an AHR2-dependent epigenetic route.

In Utero Caffeine Exposure Induces Transgenerational Effects on the Adult Heart

Each year millions of pregnant woman are exposed to caffeine, which acts to antagonize adenosine action. The long-term consequences of this exposure on the developing fetus are largely unknown, although in animal models we have found adverse effects on cardiac function. To assess if these effects are transmitted transgenerationally, we exposed pregnant mice to caffeine equivalent to 2–4 cups of coffee at two embryonic stages. Embryos (F1 generation) exposed to caffeine early from embryonic (E) day 6.5–9.5 developed a phenotype similar to dilated cardiomyopathy by 1 year of age. Embryos exposed to caffeine later (E10.5–13.5) were not affected. We next examined the F2 generation and F3 generation of mice exposed to caffeine from E10.5–13.5, as this coincides with germ cell development. These F2 generation adult mice developed a cardiac phenotype similar to hypertrophic cardiomyopathy. The F3 generation exhibited morphological changes in adult hearts, including increased mass. This report shows that in utero caffeine exposure has long-term effects into adulthood and that prenatal caffeine exposure can exert adverse transgenerational effects on adult cardiac function.

Mechanistic Evaluation of Benzo[a]pyrene's Developmental Toxicities Mediated by Reduced Cyp19a1b Activity

Benzo[a]pyrene (BaP) is a ubiquitous environmental contaminant that is both an endocrine disruptor and a carcinogen. Aromatase (CYP19) is a key enzyme in steroidogenesis that is responsible for conversion of androgens to estrogens and thus plays a key role in steroid homeostasis. We hypothesized that some of the adverse outcomes of early developmental exposure to BaP are the result of reduced Cyp19a1b activity. Our goal was to investigate the role of estrogen homeostasis during early development and determine the role of aromatase inhibition as a relevant mechanism in BaP's developmental toxicities. One-cell zebrafish embryos were injected with a Cyp19a1b-morpholino (MO) or control-MO. Other non-injected embryos were exposed to waterborne BaP, fadrozole (a Cyp19 inhibitor), estradiol (E2), BaP + E2, Cyp19a1b MO + E2, or fadrozole + E2 for 96 hours post-fertilization (hpf). Adverse outcomes were compared between treatments, and the ability of E2 co-exposure to rescue each observed dysmorphology was assessed. BaP significantly decreased cyp19a1b gene expression in 96 hpf zebrafish larvae homogenates. Concentrations of E2 in 48 hpf larvae were significantly decreased by BaP, fadrozole and Cyp19a1b-MO. Cumulative mortality of zebrafish larvae was significantly increased following BaP or fadrozole exposure or Cyp19a1b knockdown compared to controls. E2 co-treatment rescued mortality caused by 10 μg/L BaP, 10 μg/L fadrozole, or Cyp19a1b-MO. In a treatment-blinded morphological assessment of larvae at 96 hpf, several phenotypes were negatively impacted by BaP, fadrozole, or Cyp19a1b knockdown and rescued by exogenous E2 co-treatment; these included body length, optic vesicle size, swim bladder inflation, pericardial and abdominal edema, and incidence of normal larval tail shape. Abnormal pectoral fins were caused by BaP exposure only. Uninflated swim bladders were caused by all treatments including E2 alone. Our results indicate that certain BaP-mediated adverse developmental outcomes were mechanistically in accordance with BaP-mediated Cyp19a1b inhibition.

Knockdown of DNA methyltransferase 3a alters gene expression and inhibits function of embryonic cardiomyocytes

We previously found that in utero caffeine exposure causes down-regulation of DNA methyltransferases (DNMTs) in embryonic heart and results in impaired cardiac function in adulthood. To assess the role of DNMTs in these events, we investigated the effects of reduced DNMT expression on embryonic cardiomyocytes. siRNAs were used to knock down individual DNMT expression in primary cultures of mouse embryonic cardiomyocytes. Immunofluorescence staining was conducted to evaluate cell morphology. A video-based imaging assay and multielectrode array were used to assess cardiomyocyte contractility and electrophysiology, respectively. RNA-Seq and multiplex bisulfite sequencing were performed to examine gene expression and promoter methylation, respectively. At 72 h after transfection, reduced DNMT3a expression, but not DNMT1 or -3b, disrupted sarcomere assembly and decreased beating frequency, contractile movement, amplitude of field action potential, and cytosolic calcium signaling of cardiomyocytes. RNA-Seq analysis revealed that the DNMT3a-deficient cells had deactivated gene networks involved in calcium, endothelin-1, renin-angiotensin, and cardiac β-adrenergic receptor signaling, which were not inhibited by DNMT3b siRNA. Moreover, decreased methylation levels were found in the promoters of Myh7, Myh7b, Tnni3, and Tnnt2, consistent with the up-regulation of these genes by DNMT3a siRNA. These data show that DNMT3a plays an important role in regulating embryonic cardiomyocyte gene expression, morphology and function.-Fang, X., Poulsen, R. R., Wang-Hu, J., Shi, O., Calvo, N. S., Simmons, C. S., Rivkees, S. A., Wendler, C. C. Knockdown of DNA methyltransferase 3a alters gene expression and inhibits function of embryonic cardiomyocytes.

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.

Impact of prenatal polycyclic aromatic hydrocarbon exposure on behavior, cortical gene expression and DNA methylation of the Bdnf gene

Prenatal exposure to polycyclic aromatic hydrocarbons (PAH) has been associated with sustained effects on the brain and behavior in offspring. However, the mechanisms have yet to be determined. We hypothesized that prenatal exposure to ambient PAH in mice would be associated with impaired neurocognition, increased anxiety, altered cortical expression of Bdnf and Grin2b, and greater DNA methylation of Bdnf. Our results indicated that during open-field testing, prenatal PAH exposed offspring spent more time immobile and less time exploring. Females produced more fecal boli. Offspring prenatally exposed to PAH displayed modest reductions in overall exploration of objects. Further, prenatal PAH exposure was associated with lower cortical expression of Grin2b and Bdnf in males, and greater Bdnf IV promoter methylation. Epigenetic differences within the Bdnf IV promoter correlated with Bdnf gene expression, but not with the observed behavioral outcomes, suggesting that additional targets may account for these PAH-associated effects.