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

Emerging disinfection byproducts 3-bromine carbazole induces cardiac developmental toxicity via aryl hydrocarbon receptor activation in zebrafish larvae

3-bromine carbazole (3-BCZ) represents a group of emerging aromatic disinfection byproducts (DBP) detected in drinking water; however, limited information is available regarding its potential cardiotoxicity. To assess its impacts, zebrafish embryos were exposed to 0, 0.06, 0.14, 0.29, 0.58, 1.44 or 2.88 mg/L of 3-BCZ for 120 h post fertilization (hpf). Our results revealed that ≥1.44 mg/L 3-BCZ exposure induced a higher incidence of heart malformation and an elevated pericardial area in zebrafish larvae; it also decreased the number of cardiac muscle cells and thins the walls of the ventricle and atrium while increasing cardiac output and impeding cardiac looping. Furthermore, 3-BCZ exposure also exhibited significant effects on the transcriptional levels of genes related to both cardiac development (nkx2.5, vmhc, gata4, tbx5, tbx2b, bmp4, bmp10, and bmp2b) and cardiac function (cacna1ab, cacna1da, atp2a1l, atp1b2b, atp1a3b, and tnnc1a). Notably, N-acetyl-L-cysteine, a reactive oxygen species scavenger, may alleviate the failure of cardiac looping induced by 3-BCZ but not the associated cardiac dysfunction or malformation; conversely, the aryl hydrocarbon receptor agonist CH131229 can completely eliminate the cardiotoxicity caused by 3-BCZ. This study provides new evidence for potential risks associated with ingesting 3-BCZ as well as revealing underlying mechanisms responsible for its cardiotoxic effects on zebrafish embryos.

Low-dose PCB126 exposure disrupts cardiac metabolism and causes hypertrophy and fibrosis in mice

The residue of polychlorinated biphenyls (PCBs) exists throughout the environment and humans are subject to long-term exposure. As such, the potential environmental and health risk caused by low-dose exposure to PCBs has attracted much attention. 3, 3', 4, 4', 5-pentachlorobiphenyl (PCB126), the highest toxicity compound among dioxin-like-PCBs, has been widely used and mass-produced. Cardiotoxicity is PCB126's crucial adverse effect. Maintaining proper metabolism underlies heart health, whereas the impact of PCB126 exposure on cardiac metabolic patterns has yet to be elucidated. In this study, we administered 0.5 and 50 μg/kg bw of PCB126 to adult male mice weekly by gavage for eight weeks. Pathological results showed that low-dose PCB126 exposure induced heart injury. Metabolomic analysis of the heart tissue exposed to low-dose PCB126 identified 59 differential metabolites that were involved in lipid metabolism, amino acid metabolism, and the tricarboxylic acid (TCA) cycle. Typical metabolomic characteristic of cardiac hypertrophy was reflected by accumulation of fatty acids (e.g. palmitic, palmitoleic, and linoleic acid), and disturbance of carbohydrates including D-glucose and intermediates in TCA cycle (fumaric, succinic, and citric acid). Low-dose PCB126 exposure increased glycine and threonine, the amino acids necessary for the productions of collagen and elastin. Besides, PCB126-exposed mice exhibited upregulation of collagen synthesis enzymes and extracellular matrix proteins, indicative of cardiac fibrosis. Moreover, the expression of genes related to TGFβ/PPARγ/MMP-2 signaling pathway was perturbed in the PCB126-treated hearts. Together, our results reveal that low-dose PCB126 exposure disrupts cardiac metabolism correlated with hypertrophy and fibrosis. This study sheds light on the underlying mechanism of PCBs' cardiotoxicity and identifies potential sensitive biomarkers for environmental monitoring.

Toxic Effects of Polychlorinated Biphenyl Congeners and Aroclors on Embryonic Growth and Development

Polychlorinated biphenyls (PCBs) cause significant health and reproductive problems in many vertebrates. Exposure during embryogenesis likely leads to defects in organ development, compromising survival and growth through adulthood. The present study identifies the impact of PCBs on the embryonic development of key organs and resulting consequences on survival and growth. Zebrafish embryos were treated with individual PCB congeners (126 or 104) or one of 4 Aroclor mixtures (1016, 1242, 1254, or 1260) and analyzed for changes in gross embryonic morphology. Specific organs were assessed for defects during embryonic development, using a variety of transgenic zebrafish to improve organ visualization. Resulting larvae were grown to adulthood while survival and growth were assayed. Embryonic gross development on PCB treatment was abnormal, with defects presenting in a concentration-dependent manner in the liver, pancreas, heart, and blood vessel organization. Polychlorinated biphenyl 126 treatment resulted in the most consistently severe and fatal phenotypes, whereas treatments with PCB 104 and Aroclors resulted in a range of more subtle organ defects. Survival of fish was highly variable although the growth rates of surviving fish were relatively normal, suggesting that maturing PCB-treated fish that survive develop compensatory strategies needed to reach adulthood. Life span analyses of fish from embryogenesis through adulthood, as in the present study, are scarce but important for the field because they help identify foci for further studies. Environ Toxicol Chem 2021;40:187-201. © 2020 SETAC.

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.

Vezf1 regulates cardiac structure and contractile function

Background

Vascular endothelial zinc finger 1 (Vezf1) is a transcription factor previously shown to regulate vasculogenesis and angiogenesis. We aimed to investigate the role of Vezf1 in the postnatal heart.

Methods

The role of Vezf1 in regulating cardiac growth and contractile function was studied in zebrafish and in primary cardiomyocytes.

Findings

We find that expression of Vezf1 is decreased in diseased human myocardium and mouse hearts. Our experimental data shows that knockdown of zebrafish Vezf1 reduces cardiac growth and results in impaired ventricular contractile response to β-adrenergic stimuli. However, Vezf1 knockdown is not associated with dysregulation of cardiomyocyte Ca²⁺ transient kinetics. Gene ontology enrichment analysis indicates that Vezf1 regulates cardiac muscle contraction and dilated cardiomyopathy related genes and we identify cardiomyocyte Myh7/β-MHC as key target for Vezf1. We further identify a key role for an MCAT binding site in the Myh7 promoter regulating the response to Vezf1 knockdown and show that TEAD-1 is a binding partner of Vezf1.

Interpretation

We demonstrate a role for Vezf1 in regulation of compensatory cardiac growth and cardiomyocyte contractile function, which may be relevant in human cardiac disease.

Quantitative Trait Loci (QTL) Mapping

Quantitative trait loci (QTL) are genetic regions that influence phenotypic variation of a complex trait, often through genetic interactions with each other and the environment. These are commonly identified through a statistical genetic analysis known as QTL mapping. Here, I present a step-by-step, practical approach to QTL mapping along with a sample data file. I focus on methods commonly used and discoveries that have been made in fishes, and utilize a multiple QTL mapping (MQM) approach in the free software package R/qtl.

An Interrogation of Shared and Unique Copy Number Variants Across Genetically Distinct Zebrafish Strains

Zebrafish (Danio rerio) are a widely utilized model system for human disorders, but common laboratory strains have distinct behavioral and physiological differences. Accompanying these known strain differences, commonly used "wildtype" zebrafish strains have both shared and unique suites of single nucleotide polymorphisms and copy number variants (CNVs). Despite this, genomic variation is often ignored in study design, and the actual strain used is often not adequately reported. The goal of this study was to assess CNVs across three common laboratory strains of zebrafish-AB, Tubingen (TU), and WIK-and provide these data as a tool for the zebrafish community. Herein we identified 1351 CNV regions within the most recent genome assembly (GRCz11) covering 1.9% of the zebrafish genome (31.7 Mb). CNVs were found across all chromosomes, and 2200 genes (5121 transcripts) lie within ±5 kb of identified CNVs, pointing to likely cis regulatory actions of CNVs on nearby gene neighbors. We have created a Public Session accessible on the UCSC Genome Browser to view CNVs from this study titled "danRer11 zebrafish CNV across strains" as a tool for the zebrafish community.

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.

When evolution is the solution to pollution: Key principles, and lessons from rapid repeated adaptation of killifish (Fundulus heteroclitus) populations

For most species, evolutionary adaptation is not expected to be sufficiently rapid to buffer the effects of human‐mediated environmental changes, including environmental pollution. Here we review how key features of populations, the characteristics of environmental pollution, and the genetic architecture underlying adaptive traits, may interact to shape the likelihood of evolutionary rescue from pollution. Large populations of Atlantic killifish (Fundulus heteroclitus) persist in some of the most contaminated estuaries of the United States, and killifish studies have provided some of the first insights into the types of genomic changes that enable rapid evolutionary rescue from complexly degraded environments. We describe how selection by industrial pollutants and other stressors has acted on multiple populations of killifish and posit that extreme nucleotide diversity uniquely positions this species for successful evolutionary adaptation. Mechanistic studies have identified some of the genetic underpinnings of adaptation to a well‐studied class of toxic pollutants; however, multiple genetic regions under selection in wild populations seem to reflect more complex responses to diverse native stressors and/or compensatory responses to primary adaptation. The discovery of these pollution‐adapted killifish populations suggests that the evolutionary influence of anthropogenic stressors as selective agents occurs widely. Yet adaptation to chemical pollution in terrestrial and aquatic vertebrate wildlife may rarely be a successful “solution to pollution” because potentially adaptive phenotypes may be complex and incur fitness costs, and therefore be unlikely to evolve quickly enough, especially in species with small population sizes.

Sequencing and de novo draft assemblies of a fathead minnow (Pimephales promelas) reference genome

The present study was undertaken to provide the foundation for development of genome-scale resources for the fathead minnow (Pimephales promelas), an important model organism widely used in both aquatic toxicology research and regulatory testing. The authors report on the first sequencing and 2 draft assemblies for the reference genome of this species. Approximately 120× sequence coverage was achieved via Illumina sequencing of a combination of paired-end, mate-pair, and fosmid libraries. Evaluation and comparison of these assemblies demonstrate that they are of sufficient quality to be useful for genome-enabled studies, with 418 of 458 (91%) conserved eukaryotic genes mapping to at least 1 of the assemblies. In addition to its immediate utility, the present work provides a strong foundation on which to build further refinements of a reference genome for the fathead minnow.

Early life exposure to PCB126 results in delayed mortality and growth impairment in the zebrafish larvae

The occurrence of chronic or delayed toxicity resulting from the exposure to sublethal chemical concentrations is an increasing concern in environmental risk assessment. The Fish Embryo Toxicity (FET) test with zebrafish provides a reliable prediction of acute toxicity in adult fish, but it cannot yet be applied to predict the occurrence of chronic or delayed toxicity. Identification of sublethal FET endpoints that can assist in predicting the occurrence of chronic or delayed toxicity would be advantageous. The present study characterized the occurrence of delayed toxicity in zebrafish larvae following early exposure to PCB126, previously described to cause delayed effects in the common sole. The first aim was to investigate the occurrence and temporal profiles of delayed toxicity during zebrafish larval development and compare them to those previously described for sole to evaluate the suitability of zebrafish as a model fish species for delayed toxicity assessment. The second aim was to examine the correlation between the sublethal endpoints assessed during embryonal and early larval development and the delayed effects observed during later larval development. After exposure to PCB126 (3-3000ng/L) until 5 days post fertilization (dpf), larvae were reared in clean water until 14 or 28 dpf. Mortality and sublethal morphological and behavioural endpoints were recorded daily, and growth was assessed at 28 dpf. Early life exposure to PCB126 caused delayed mortality (300 ng/L and 3000 ng/L) as well as growth impairment and delayed development (100 ng/L) during the clean water period. Effects on swim bladder inflation and cartilaginous tissues within 5 dpf were the most promising for predicting delayed mortality and sublethal effects, such as decreased standard length, delayed metamorphosis, reduced inflation of swim bladder and column malformations. The EC50 value for swim bladder inflation at 5 dpf (169 ng/L) was similar to the LC50 value at 8 dpf (188 and 202 ng/L in two experiments). Interestingly, the patterns of delayed mortality and delayed effects on growth and development were similar between sole and zebrafish. This indicates the comparability of critical developmental stages across divergent fish species such as a cold water marine flatfish and a tropical freshwater cyprinid. Additionally, sublethal effects in early embryo-larval stages were found promising for predicting delayed lethal and sublethal effects of PCB126. Therefore, the proposed method with zebrafish is expected to provide valuable information on delayed mortality and delayed sublethal effects of chemicals and environmental samples that may be extrapolated to other species.

Comparative mapping for bighead carp (Aristichthys nobilis) against model and non-model fishes provides insights into the genomic evolution of cyprinids

Comparative mapping provides an efficient method to connect genomes of non-model and model fishes. In this study, we used flanking sequences of the 659 microsatellites on a genetic map of bighead carp (Aristichthys nobilis) to comprehensively study syntenic relationships between bighead carp and nine model and non-model fishes. Of the five model and two food fishes with whole genome data, Cyprinus carpio showed the highest rate of positive BLAST hits (95.3 %) with bighead carp map, followed by Danio rerio (70.9 %), Oreochromis niloticus (21.7 %), Tetraodon nigroviridis (6.4 %), Gasterosteus aculeatus (5.2 %), Oryzias latipes (4.7 %) and Fugu rubripes (3.5 %). Chromosomal syntenic analyses showed that inversion was the basic chromosomal rearrangement during genomic evolution of cyprinids, and the extent of inversions and translocations was found to be positively correlated with evolutionary relationships among fishes studied. Among the five investigated cyprinids, linkage groups (LGs) of bighead carp, Hypophthalmichthys molitrix and Ctenopharyngodon idella exhibited a one-to-one relationship. Besides, LG 9 of bighead carp and homologous LGs of silver carp and grass carp all corresponded to the chromosomes 10 and 22 of zebrafish, suggesting that chromosomal fission may have occurred in the ancestor of zebrafish. On the other hand, LGs of bighead carp and common carp showed an approximate one-to-two relationship with extensive translocations, confirming the occurrence of a 4th whole genome duplication in common carp. This study provides insights into the understanding of genome evolution among cyprinids and would aid in transferring positional and functional information of genes from model fish like zebrafish to non-model fish like bighead carp.

Zebrafish: A marvel of high-throughput biology for 21st century toxicology

The evolutionary conservation of genomic, biochemical and developmental features between zebrafish and humans is gradually coming into focus with the end result that the zebrafish embryo model has emerged as a powerful tool for uncovering the effects of environmental exposures on a multitude of biological processes with direct relevance to human health. In this review, we highlight advances in automation, high-throughput (HT) screening, and analysis that leverage the power of the zebrafish embryo model for unparalleled advances in our understanding of how chemicals in our environment affect our health and wellbeing.

Genetic variation at aryl hydrocarbon receptor (AHR) loci in populations of Atlantic killifish (Fundulus heteroclitus) inhabiting polluted and reference habitats

Background

The non-migratory killifish Fundulus heteroclitus inhabits clean and polluted environments interspersed throughout its range along the Atlantic coast of North America. Several populations of this species have successfully adapted to environments contaminated with toxic aromatic hydrocarbon pollutants such as polychlorinated biphenyls (PCBs). Previous studies suggest that the mechanism of resistance to these and other “dioxin-like compounds” (DLCs) may involve reduced signaling through the aryl hydrocarbon receptor (AHR) pathway. Here we investigated gene diversity and evidence for positive selection at three AHR-related loci (AHR1, AHR2, AHRR) in F. heteroclitus by comparing alleles from seven locations ranging over 600 km along the northeastern US, including extremely polluted and reference estuaries, with a focus on New Bedford Harbor (MA, USA), a PCB Superfund site, and nearby reference sites.

Results

We identified 98 single nucleotide polymorphisms within three AHR-related loci among all populations, including synonymous and nonsynonymous substitutions. Haplotype distributions were spatially segregated and F-statistics suggested strong population genetic structure at these loci, consistent with previous studies showing strong population genetic structure at other F. heteroclitus loci. Genetic diversity at these three loci was not significantly different in contaminated sites as compared to reference sites. However, for AHR2 the New Bedford Harbor population had significant F_ST values in comparison to the nearest reference populations. Tests for positive selection revealed ten nonsynonymous polymorphisms in AHR1 and four in AHR2. Four nonsynonymous SNPs in AHR1 and three in AHR2 showed large differences in base frequency between New Bedford Harbor and its reference site. Tests for isolation-by-distance revealed evidence for non-neutral change at the AHR2 locus.

Conclusion

Together, these data suggest that F. heteroclitus populations in reference and polluted sites have similar genetic diversity, providing no evidence for strong genetic bottlenecks for populations in polluted locations. However, the data provide evidence for genetic differentiation among sites, selection at specific nucleotides in AHR1 and AHR2, and specific AHR2 SNPs and haplotypes that are associated with the PCB-resistant phenotype in the New Bedford Harbor population. The results suggest that AHRs, and especially AHR2, may be important, recurring targets for selection in local adaptation to dioxin-like aromatic hydrocarbon contaminants.

Targeted approach to identify genetic loci associated with evolved dioxin tolerance in Atlantic killifish (Fundulus heteroclitus)

Background

The most toxic aromatic hydrocarbon pollutants are categorized as dioxin-like compounds (DLCs) to which extreme tolerance has evolved independently and contemporaneously in (at least) four populations of Atlantic killifish (Fundulus heteroclitus). Surprisingly, the magnitude and phenotype of DLC tolerance is similar among these killifish populations that have adapted to varied, but highly aromatic hydrocarbon-contaminated urban/industrialized estuaries of the US Atlantic coast. Multiple tolerant and neighboring sensitive killifish populations were compared with the expectation that genetic loci associated with DLC tolerance would be revealed.

Results

Since the aryl hydrocarbon receptor (AHR) pathway partly or fully mediates DLC toxicity in vertebrates, single nucleotide polymorphisms (SNPs) from 42 genes associated with the AHR pathway were identified to serve as targeted markers. Wild fish (N = 36/37) from four highly tolerant killifish populations and four nearby sensitive populations were genotyped using 59 SNP markers. Similar to other killifish population genetic analyses, strong genetic differentiation among populations was detected, consistent with isolation by distance models. When DLC-sensitive populations were pooled and compared to pooled DLC-tolerant populations, multi-locus analyses did not distinguish the two groups. However, pairwise comparisons of nearby tolerant and sensitive populations revealed high differentiation among sensitive and tolerant populations at these specific loci: AHR 1 and 2, cathepsin Z, the cytochrome P450s (CYP1A and 3A30), and the NADH dehydrogenase subunits. In addition, significant shifts in minor allele frequency were observed at AHR2 and CYP1A loci across most sensitive/tolerant pairs, but only AHR2 exhibited shifts in the same direction across all pairs.

Conclusions

The observed differences in allelic composition at the AHR2 and CYP1A SNP loci were identified as significant among paired sensitive/tolerant populations of Atlantic killifish with multiple statistical tests. The genetic patterns reported here lend support to the argument that AHR2 and CYP1A play a role in the adaptive response to extreme DLC contamination. Additional functional assays are required to isolate the exact mechanism of DLC tolerance.

Genetic susceptibility to teratogens: state of the art

There is evidence that the susceptibility to the teratogenic effect of drugs within human populations varies extremely from one individual to another, even after identical exposures. One of the factors that may explain these interindividual differences is the genetic makeup in the pharmacokinetics and pharmacodynamics of the respective drugs. In fact, both maternal and embryonic/fetal genotypes can affect placental transport, absorption, metabolism, distribution and receptor binding of an agent, influencing its teratogenicity. We have reviewed the literature and commented on the reported correlations between genetic factors and drug-induced birth defects. There is still a clear lack of knowledge regarding this issue and the available data are often conflicting. However, the identification of specific polymorphisms associated with predisposition to teratogenesis may allow in the future the development of personalized non-teratogenic therapies for pregnant women.

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.

Common mechanism underlies repeated evolution of extreme pollution tolerance

Human alterations to the environment can exert strong evolutionary pressures, yet contemporary adaptation to human-mediated stressors is rarely documented in wildlife populations. A common-garden experimental design was coupled with comparative transcriptomics to discover evolved mechanisms enabling three populations of killifish resident in urban estuaries to survive normally lethal pollution exposure during development, and to test whether mechanisms are unique or common across populations. We show that killifish populations from these polluted sites have independently converged on a common adaptive mechanism, despite variation in contaminant profiles among sites. These populations are united by a similarly profound desensitization of aryl-hydrocarbon receptor-mediated transcriptional activation, which is associated with extreme tolerance to the lethal effects of toxic dioxin-like pollutants. The rapid, repeated, heritable and convergent nature of evolved tolerance suggests that ancestral killifish populations harboured genotypes that enabled adaptation to twentieth-century industrial pollutants.