Original Research ARTICLE
Developmental Dioxin Exposure Alters the Methylome of Adult Male Zebrafish Gonads
- 1Department of Pharmacology, Wayne State University, Detroit, MI, United States
- 2Institute of Environmental Health Sciences, Wayne State University, Detroit, MI, United States
- 3Applied Genome Technology Center, School of Medicine, Wayne State University, Detroit, MI, United States
2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a persistent environmental toxicant and endocrine disrupting compound with reproductive and developmental effects in humans and model organisms, including zebrafish. Our previous microarray and histological studies found defects in spermatogenesis and fertility of zebrafish in response to acute developmental TCDD exposure. These effects are apparent following exposure during reproductive development, modeling fetal basis of adult-onset disease. Some outcomes of these previous studies (reduced fertility, changes in sex ratio, transcriptomic alterations) are also transgenerational – persisting to unexposed generations – through the male germline. We hypothesized that DNA methylation could be a possible mechanism for these reproductive effects and performed whole genome bisulfite sequencing (WGBS), which identifies whole genome DNA methylation status at the base pair level, on testes of adult zebrafish exposed to TCDD (two separate hour-long exposures to 50 pg/mL TCDD at 3 and 7 weeks post fertilization). In response to TCDD exposure, multiple genes were differentially methylated; many of which are involved in reproductive processes or epigenetic modifications, suggesting a role of DNA methylation in later-life health outcomes. Additionally, several differentially methylated genes corresponded with gene expression changes identified in TCDD-exposed zebrafish testes, indicating a potential link between DNA methylation and gene expression. Ingenuity pathway analysis of WGBS and microarray data revealed genes involved in reproductive processes and development, RNA regulation, the cell cycle, and cellular morphology and development. We conclude that site-specific changes in DNA methylation of adult zebrafish testes occur in response to acute developmental TCDD exposure.
2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD), a by-product of industrial processes such as fossil fuel combustion, waste incineration, and metal production, is a widespread and persistent environmental toxicant. As a well-characterized model aryl-hydrocarbon receptor agonist, this endocrine disrupting compound affects a variety of organ systems, leading to human health outcomes including cancer, type II diabetes, and metabolic disease (Steenland et al., 1999; Manikkam et al., 2012; Warner et al., 2013). People exposed to TCDD suffer decreased semen quality, have longer menstrual cycles, and produce more female than male offspring (i.e., altered sex ratios) (Signorini et al., 2000; Eskenazi et al., 2002; Mocarelli et al., 2008). Similar reproductive outcomes have been demonstrated in model organisms including zebrafish (King Heiden et al., 2009; Manikkam et al., 2012; Baker et al., 2013).
Due to their short generation time, high fecundity, and externally developing embryos, zebrafish are an excellent transgenerational model for studying the reproductive and developmental outcomes of chemical exposure. Our earlier research demonstrated impaired fertility in adult zebrafish in response to TCDD exposure during sexual differentiation (two separate hour-long exposures to 50 pg/mL TCDD at 3 and 7 weeks post fertilization), modeling fetal basis of adult-onset disease (Baker et al., 2013). Several reproductive outcomes (altered sex ratio, decreased egg release, and transcriptional changes in reproductive genes) persisted through two subsequent generations (Baker et al., 2014a,b, 2016; Meyer et al., 2018). Reduced reproductive capacity was characterized by decreases in both egg release (a product of courtship between male and female zebrafish) and percentage of fertilized eggs, which presented across three generations (F0–F2). We determined these effects were male-mediated by outcrossing control and exposed lineage fish from each generation; only males from the exposed lineage demonstrated transgenerational reproductive effects (Baker et al., 2014b). Subsequent microarray and histologic analyses revealed differentially expressed genes involved in reproduction and defects in spermatogenesis, indicating a possible mechanism for the male-mediated decreased fertility (Baker et al., 2016; Meyer et al., 2018).
Parental exposures to environmental chemicals and other stressors, including paternal exposures, can contribute to adverse outcomes in offspring and subsequent generations (Jirtle and Skinner, 2007; Rodgers et al., 2013; Robledo et al., 2015). The mechanisms of adult-onset and transgenerational TCDD effects are unknown but thought to be epigenetic as studies have demonstrated changes in common epigenetic mechanisms such as DNA and histone methylation, as well as differential expression of DNA methyltransferases in response to TCDD exposure (Manikkam et al., 2012; Olsvik et al., 2014; Aluru et al., 2015; Ma et al., 2015; Baker et al., 2016; Sanabria et al., 2016). Because the DNA methylation pattern in zebrafish is passed down paternally through the sperm (Potok et al., 2013; Jiang et al., 2014), inheritance of epimutations in the DNA methylome is a promising mechanism of transgenerational male-mediated reproductive defects resulting from TCDD exposure. Other studies examining DNA methylation in zebrafish in response to TCDD exposure (Olsvik et al., 2014; Aluru et al., 2015) focused on either locus-specific or global assays, both of which have limitations for informing DNA methylation status. Reduced representative bisulfite sequencing improves on these methods, but still only focuses on a subset of CpG islands and promoter regions. In this study, we utilize whole genome bisulfite sequencing (WGBS) which provides base pair level DNA methylation status of the entire genome, eliminating potential masking of bidirectional changes in global analysis and providing more genome wide coverage than locus specific methods (Doherty and Couldrey, 2014).
The purpose of this study is to utilize WGBS to expand the current knowledge of TCDD-induced epigenetic effects, demonstrating that sublethal TCDD exposure during zebrafish gonadal differentiation and maturation causes changes in DNA methylation at specific genome loci which persist to adulthood. As our previous work revealed transcriptomic alterations in zebrafish testes induced by TCDD exposure, we also interpret our findings of changes in the testicular methylome within the context of this data (Baker et al., 2016). To our knowledge, we have performed the first WGBS on zebrafish in response to TCDD exposure. These data will provide valuable insight into the epigenetic effects of environmental endocrine-disrupting compounds as potential mechanisms of adult-onset and multigenerational disease.
Materials and Methods
Zebrafish (AB strain) were kept at 28°C in buffered reverse osmosis water (60 mg/L Instant Ocean Salts; Aquarium Systems, Mentor, OH, United States) with a standard light/dark cycle of 14/10 h and fed Aquatox Fish Diet flakes (Zeigler, PA, United States) twice per day, supplemented with brine shrimp. Fish were raised in beakers with daily water changes of 40–60% at a density of five fish per 400 mL beaker between 3 and 6 weeks post-fertilization (wpf), and five fish per 800 mL beaker between 6 and 9 wpf. Adult fish were raised on a recirculating system at a maximum density of five fish per liter until euthanization at 1 year post fertilization. Fish were euthanized with tricaine methanesulfonate (1.67 mg/mL). Animal use protocols were approved by the Institutional Animal Care and Use Committees at Wayne State University and the University of Wisconsin-Madison, according to the National Institutes of Health Guide to the Care and Use of Laboratory Animals (Protocol No. M00489).
Exposures were performed as previously described in Baker et al. (2013). TCDD (>99% purity; Chemsyn) was used as a 0.4 ng/μL stock solution in dimethyl sulfoxide (DMSO). Zebrafish were exposed at 3 wpf and again at 7 wpf to waterborne TCDD (50 pg/mL) or vehicle (0.1% DMSO) for 1 h in glass beakers with gentle rocking. The number of fish per volume of dosing solution was 1 fish/mL at 3 wpf and 1 fish/2 mL at 7 wpf. All results are derived from three independent TCDD exposure experiments performed in successive blocks. The fish in the current experiment were collected and stored from the same blocks of the previously published microarray and histology studies (Baker et al., 2013, 2016).
Testes were extracted from 1-year-old TCDD- and DMSO vehicle control-treated zebrafish and flash frozen in liquid nitrogen. Samples were kept at -80°C until DNA isolation. DNA isolation was performed on four DMSO control and four TCDD exposed samples using the BioRobot EZ1 workstation following the provided protocol from Qiagen (Hilden, Germany). Concentration and quality of DNA was measured using Qubit (Invitrogen, Carlsbad, CA, United States) and DropSense 96 (Trinean, Gentbrugge, Belgium), respectively (Supplementary Table S1). DNA was collected from different fish than those used for RNA collection and microarray analysis.
Microarray data were previously analyzed according to Baker et al. (2016). Microarray was performed on F0 zebrafish testes. Data were uploaded to NCBI GEO database (GSE77335). Microarray data were analyzed with one-way between subject ANOVA using the Transcriptome Analysis Console (TAC, Affymetrix). Genes uploaded into Ingenuity Pathway Analysis (IPA; Qiagen Bioinformatics; Redwood City, CA, United States) for pathway analysis were defined as significantly altered with a p-value ≤ 0.05 and absolute fold change ≥ 1.5. Gene expression was validated by qPCR of 22 genes of interest. The validation and IPA analysis data was reported in Baker et al. (2016). Although all microarray data were previously uploaded, fold changes for a subset of genes not reported in the previous paper are discussed here for the first time.
Whole Genome Bisulfite Sequencing
Bisulfite conversion of methylated cytosines of DNA from zebrafish testes was performed using the EZ DNA Methylation Kit (Zymo, Irvine, CA, United States). The resulting bisulfite-converted ssDNA was converted to Illumina libraries with the TruSeq DNA Methylation Kit (Illumina, San Diego, CA, United States). DNA libraries were sequenced with 100 bp paired-end reads on an Illumina HiSeq 2500 run in high output mode. One of the TCDD exposed samples did not pass our quality control cutoffs post amplification and was removed from the study. Reads from remaining four control and three TCDD-treated fish were aligned (bismark 0.18.1; bowtie2 2.2.3) to the zebrafish genome (dR10) and differentially methylated sites were determined (methylKit 1.6.3) between conditions (Supplementary Table S2) (Langmead et al., 2009; Krueger and Andrews, 2011; Akalin et al., 2012). Methylated sites (percent methylation change > 5%; p-value ≤ 0.05) were annotated to the nearest gene. 5% methylation change was used as a cutoff after referencing previous studies that used a 5% cutoff for significant methylation change or considered changes less than 5% significant (Yang et al., 2006; Byun et al., 2013; Dayeh et al., 2014). Differentially expressed genes from the microarray data were overlapped with annotations for differentially methylated sites. Individual replicates from each condition were kept separate and not pooled.
Ingenuity Pathway Analysis
Genes associated with differential methylation (methylation change > 5%, p-value ≤ 0.05) were converted to homologous human genes and uploaded into IPA software and analyzed using RefSeq ID as the identifier. Fifty-three molecules were available for pathway analysis of enriched disease and biologic functions (Supplementary Table S4).
All DNA methylation data can be found in Supplementary Table S3. WGBS of control and TCDD-treated zebrafish revealed that 12.38% of total cytosines were methylated in control samples and 12.57% were methylated in TCDD exposed samples. In a CpG context, 86.55% of cytosines were methylated in control samples and 87.57% of cytosines were methylated in TCDD exposed samples. These changes in global methylation are not statistically significant with p-values of 0.48 for total cytosines and 0.54 for cytosines in a CpG context. There were, however, significant alterations in DNA methylation in specific genomic locations between the control and exposed groups.
A total of 397 sites were significantly differentially methylated between control and exposed fish. Of these, 284 were hypomethylated and 113 were hypermethylated (Figure 1). Of these differentially methylated sites, 89 sites were found <100 kb downstream of genes, followed by 97 sites <100 kb upstream of genes. Within genes, more differentially methylated sites were found within introns (42) rather than within exons (15) (Figure 2).
Figure 1. Identification of significantly differentially methylated sites between dioxin treated and untreated zebrafish gonads. Volcano plot depicting percent methylation change and p-value of methylated sites (circles). Significantly changed CpG sites (percent methylation change 5%; p-value ≤ 0.05) are in yellow.
Figure 2. Distribution of significant (percent methylation change ≥5%; p-value ≤ 0.05) differentially methylated sites relative to nearest gene.
Differentially methylated genes of interest related to reproduction, RNA regulation, and epigenetic processes are listed in Table 1. Of the genes related to reproduction, bambia, chrna4, hug, and fmr1 are all involved in spermatogenesis, which is consistent with our previous studies that observed defects in spermatogenesis in response to TCDD exposure. Fmr1 is also involved in RNA regulation as are dhx36, klf15, and pym1, while eya3, wdr5, and zgc:162967 are involved in histone modification.
Table 1. Genes with significant changes in DNA methylation (p-value <0.05, methylation change >5%) that are involved in reproduction, RNA regulation, or epigenetic regulation.
Differentially Methylated Regions
In addition to individual CpG sites, we also discovered differentially methylated regions (DMRs). At a window size of 1,000 bp and a step increment of 500 bp, there were a total of 148 statistically significant (p-value < 0.05) DMRs linked to 64 different genes (several DMRs were linked to the same gene). Maps of individual methylated sites linked to three genes of interest are shown in Figure 3. 41 CpG-linked genes were also linked to at least one DMR, and 23 DMR-linked genes were found that were not linked to any individual CpG sites; for example, the reproductive gene dnaja2 (Figure 3A). Notably, several reproductive genes are linked to both individual CpG sites as well as DMRs: fmr1, ndufc2, and stoml2 (Figure 3C). Pym1 (Figure 3B) and dhx36, which are both involved in RNA regulatory mechanisms, were also linked to both CpG sites and DMRs. Finally, crestin, which is both differentially expressed and linked to a CpG site, is also linked to a DMR.
Figure 3. Maps of individual methylated sites linked to three genes: (A) dnaja2, (B) pym1, and (C) stoml2. Vertical black bars indicate individual DNA methylation sites and height of the bars represents methylation value. Blue dots represent differentially methylated CpG sites, and bolded horizontal black bars represent differentially methylated regions. Each row represents a different sample.
We performed IPA on genes linked to differentially methylated sites in TCDD exposed zebrafish testes and found 405 total associated diseases and annotation pathways. Reproductive system disease (p-values: 4.79E-02 – 1.46E-03), as well as reproductive development and function (p-values: 4.93E-02 – 1.46E-03) were reported under top diseases and biological functions. Two of the top networks affected included gene expression (Gene Expression, Embryonic Development, Nervous System Development and Function, Hypersensitivity) and cell cycle (Cell Cycle, Auditory Disease, Cellular Development). Table 2 shows the diseases or functions from IPA analyses that are involved in reproduction, gene expression regulation, or cellular morphology and development.
Table 2. Diseases or functions from IPA analysis of genes with significant changes in DNA methylation between control and TCDD-exposed groups that are involved in reproduction, gene expression regulation, or cellular morphology and development.
Differentially Expressed Genes Involved in Epigenetic Regulation
Our microarray analysis identified five differentially methylated genes involved in epigenetic regulation (Table 3): wdr5 (fold change: -1.26, p-value: 0.029), setd1bb (fold change: -1.31, p-value: 0.049), dnmt5 (fold change: -1.49, p-value: 0.04), dnmt3 (fold change: -1.44, p-value: 0.03), and dnmt3b (fold change: -1.32, p-value: 0.02). Wdr5 and setd1bb are histone methyltransferases and dnmt5, dnmt3, and dnmt3b are DNA methyltransferases. Wdr5 is likely involved in positioning lysine 4 of histone H3 for methylation and setd1bb trimethylates this site. These DNA methyltransferases are involved in de novo DNA methylation which occurs during development, as opposed to maintenance methylation.
Table 3. Genes from microarray with significant differential expression with roles in epigenetic regulation.
Relationship Between DNA Methylation and Gene Expression
Genes with proximally modified DNA methylation sites and differential expression between TCDD-exposed and control testes are displayed in Table 4. The differentially methylated sites near riox1, crestin, wdr5, and sox21b are upstream of their respective genes while the differentially methylated sites near flrt2, rab11bb, and gabra1 are located downstream of their respective genes.
Table 4. Genes with significant differential expression and proximal changes in DNA methylation status.
Flrt2 and crestin, involved in embryonic development (Luo et al., 2001; Müller et al., 2011), and the transcription factor sox21b, involved in neuronal differentiation (Sandberg et al., 2005), have methylation changes directly related to gene expression. Flrt2 has increased expression and hyper-methylation at three CpG sites; likewise, sox21b is upregulated and hyper-methylated at two associated DMRs. Crestin expression, however, is reduced, with hypo-methylation of both the neighboring CpG site and associated DMR. Gabra1, a GABA receptor subunit (Bhat et al., 2018) has decreased expression while the associated DMR is hyper-methylated, demonstrating a relationship between methylation and gene expression. Riox1 and wdr5, respectively, a histone demethylase and histone methyltransferase (Aravind et al., 2011; Tao et al., 2013), exhibit decreased expression while the endosomal GTPase rab11bb (Clark et al., 2011) has increased expression, but all three are surrounded by a mixture of hyper and hypo-methylated sites (riox and wdr5) and DMRs (rab11bb).
Whole genome methylation analysis revealed both DMR- and CpG-specific changes in the DNA methylation profile of adult zebrafish developmentally exposed to sublethal levels of TCDD. Previously, we found that this exposure paradigm induced transgenerational adult infertility mediated through the male germ line (Baker et al., 2014a,b). As the exposure occurred during reproductive development, the latency to effect (3–10 months) and transmission of these outcomes across generations strongly implicate the presence of persistent epimutations in reproductive pathways as a mechanism for the infertility phenotype. Our promising findings of TCDD-induced differential methylation in reproductive disease, gene expression, and cell morphology and development pathways correspond with our phenotypic and transcriptomic outcomes and are not unexpected, considering recent findings that disruption of DNA methylation during gonad development can lead to changes in reproductive and epigenetic genes, such as increases in dnmt1 in adult zebrafish (Ribas et al., 2017).
About two-thirds of DMR-associated genes overlap with individual CpG-associated genes, including genes involved in mitochondrial function. For example, stoml2, associated with a mix of hypo- and hyper-methylated sites, regulates mitochondrial biogenesis (Christie et al., 2011) along with ndufc2, a member of the electron transport chain (Stroud et al., 2016) associated with hypo-methylated regions. Dysregulation in mitochondrial pathways can impact multiple physiological processes, including reproductive health and fertility; thus, epimutations in these pathways have potential to impair fertility through several distinct mechanisms. Spermatogenesis is a closely controlled process with high energy demands, thus any shifts or deficits in energy balance due to altered quantity and output of gonadal mitochondria may arrest sperm cell maturation (Nakada et al., 2006). Additionally, impairment of ndufc2, as demonstrated in a rodent knock-out model, can induce mitochondrial damage and oxidative stress, which is known to damage the sperm membrane and induce infertility (Tremellen, 2008; Ramalho-Santos et al., 2009; Raffa et al., 2017). Another gene implicated in oxidative stress response, dnaja2, is an hsp70 co-chaperone associated solely with hypo-methylated DMRs. Dnaja2 plays a role in maintaining protein stability and regulating apoptosis; knockout of a closely related DnaJ homolog resulted in disrupted Sertoli-germ cell adherens junctions and severely impaired spermatogenesis, suggesting an important role in sperm development and androgen signaling (Terada et al., 2005; Meccariello et al., 2014).
Analysis of differentially methylated individual CpGs uncovered linked genes associated with several pathways dysregulated in infertility, including spermatogenesis and glucose/lipid metabolism. Bambia, associated with a mix of hypo- and hyper-methylated sites, and chrna4b, associated with a single hypo-methylated site, are implicated in the initiation and progression of spermatogenesis in germ cells and Leydig cells, respectively, suggesting that epimutations in these genes may contribute to the infertility phenotype (Loveland et al., 2003; Hamra et al., 2004; Bray et al., 2005; Ge et al., 2005; Rolland et al., 2009). Glucose and lipid metabolism, two components of energy metabolism, cooperate to regulate sex hormone production and initiation of gametogenesis (Fontana and Della Torre, 2016). CpG methylation changes linked to several metabolism genes, namely klf15, a transcription factor involved in glucose homeostasis, lipid accumulation, and testosterone production (Gray et al., 2002; Du et al., 2009), and nr1h3, which controls cholesterol homeostasis and steroidogenic gene expression (Zhao and Dahlman-Wright, 2010), suggest dysregulation of the tightly regulated feedback loop between the glucose and lipid metabolism pathways and the reproductive system. Klf15 is linked to a hyper-methylated CpG site, whereas nr1h3 is linked to a mixture of hyper- and hypo-methylated sites. An imbalance in these pathways can result in outcomes such as impaired gonadal somatic (Leydig and Sertoli) cell function, decreased sperm quality, and subsequent infertility (Alves et al., 2013; Ding et al., 2015; Parhofer, 2015), aligning with our phenotypic and transcriptomic results of decreased fertilization, a shift in germ cell ratio toward less mature cells, and altered expression of genes linked to spermatogenesis, testicular development, lipid metabolism, and steroidogenesis pathways (Baker et al., 2016).
In addition to reproductive genes, we observed differential methylation of a subset of histone methyltransferases, histone demethylases, miRNAs, and RNA modifying genes. Although histone methylation interacts both independently and reciprocally with DNA methylation in controlling gene expression and chromatin state (Rose and Klose, 2014), few studies have assessed the impact of TCDD exposure on histone methylation. Of particular interest is wdr5, a core component of several histone methyltransferase complexes, including the SET1 complex, which primarily modulates methylation of the classic gene activation mark H3K4 (Barski et al., 2007; Aravind et al., 2011). Along with a crucial role in global gene activation, wdr5 is involved in programming embryonic stem cell differentiation and self-renewal (Ang et al., 2011), as well as epigenetically controlling gluconeogenic pathways (Ravnskjaer et al., 2013), that are implicated in fertility and reproduction. Wdr5 and riox1, a histone demethylase, are both differentially methylated (mixture of hypo- and hyper-methylated CpG sites) and downregulated, suggesting that their methylation state and expression profile are directly linked, and that decreases in histone methyltransferase and demethylase function are contributing to the fertility-linked epimutations present in the genome.
Transcriptomic data also confirms TCDD-induced downregulation of several other critical histone and DNA methyltransferase genes, including the histone methyltransferase setd1bb, involved in H3K4 methylation alongside wdr5, and the de novo DNA methyltransferases dnmt3, dnmt3b, and dnmt5, involved in establishment of DNA methylation patterns during gametogenesis (Campos et al., 2012). As histone and DNA methyltransferases mediate the methylome, toxicant-induced changes in their expression could account for the changes in methylation pattern. Decreases in DNA and histone methyltransferase activity would generally be expected to result in global and/or gene-specific hypomethylation and subsequent increases in gene expression. However, we found no global hypomethylation and very limited site-specific hypomethylation/corresponding upregulation of gene expression. These outcomes align with Aluru et al. (2015), which reported decreases in several dnmt3 genes (dnmt3a1, 3b1, and 3b4) after developmental exposure to TCDD in zebrafish, and both hyper- and hypo-methylation in the promoter region of AhR target genes that overall did not correlate with gene expression. Somm et al. (2013) also reported limited relationships between gene expression and methylation of imprinted genes in the sperm of male mice gestationally exposed to TCDD. Conversely, both Papoutsis et al. (2015) and Sekaran and Jagadeesan (2015) found that gestational exposure of rodents to TCDD and DEHP, respectively, resulted in gene-specific hyper-methylation and corresponding decreases in gene expression, with increased Dnmt levels reported by the latter. These outcomes support a multifaceted interpretation of the methylation-gene expression relationship, suggesting not only a variable, gene-dependent response, but also the possible involvement of multiple epigenetic processes in mediating transcriptomic response to toxicant insult. Results showing that the degree of H3K4 methylation can regulate de novo DNA methyltransferase activity (Rose and Klose, 2014) strengthen the link between DNA and histone methylation as potential mechanisms of epimutation with respect to our methylomic and transcriptomic findings. However, as the location and extent of histone methylation in TCDD-exposed fish remains to be studied, inferring further directional relationships from the complex interplay between DNA methylation, histone methylation, and gene expression remains a challenge. Several miRNAs, including mir2185-1, are repressive epigenetic factors that interact with histone modification and DNA methylation that were also differentially methylated by TCDD exposure, and may thus contribute to modifications of the transcriptome (Silahtaroglu and Stenvang, 2010; Osella et al., 2014). Mir2185-1 is associated with several hyper-methylated DMRs and a mix of hyper- and hypo-methylated individual CpGs. Taken together, TCDD-induced differential methylation of epigenetic genes has potential to alter both genome-wide and pathway-specific patterns of methylation through a variety of inter-related mechanisms, leading to transcriptomic changes and infertility.
Several differentially methylated genes were implicated in the regulation of gene expression. Fmr1, pym1, and dhx36 are associated with both differentially methylated individual CpG sites and DMRs, whereas elavl1 is specifically associated with individual CpG sites. Although not epigenetic factors by definition, these RNA-modifying genes regulate the transcriptome and proteome; therefore, altered methylation and potentially altered expression of these genes can have broad consequences for overall gene transcription and translation (Bächner et al., 1993; Diem et al., 2007; Chi et al., 2011; Ferder et al., 2013; Booy et al., 2014). In some cases, the regulatory function of these genes directly impacts reproductive outcomes; for example, post-transcriptional activity of elavl1 is essential for both meiotic division and spermatid differentiation, thus leading to extensive gonadal cell death and subsequent azoospermia in a mouse knockout model (Chi et al., 2011). Likewise, fmr1 interacts with chromatin to regulate DNA damage during meiosis, resulting in macroorchidism, impaired spermatogenesis, and maturation arrest when downregulated (Slegtenhorst-Eegdeman et al., 1998; Tian et al., 2013; Alpatov et al., 2014). However, the majority of these genes (excepting dhx36) are linked to a mix of hyper- and hypo-methylated sites, posing a challenge in the interpretation of downstream effects. Overall, differential methylation of these genes corresponds with our findings of dysregulation of cell morphology and development pathways, with potential for broad defects in gonadal structure and function leading to infertility outcomes in adulthood (Baker et al., 2016).
No global shift in methylation was detected; however, a substantial number of site-specific changes were present, highlighting the advantage of WGBS as a genome-wide, high-specificity method. EDC exposure has both global and gene-specific effects on the methylome (Anway et al., 2005; Bollati et al., 2007; Dolinoy et al., 2007; Pilsner et al., 2007; Onishchenko et al., 2008), but outcomes vary widely depending on targeted pathways and critical windows of exposure. TCDD-mediated AhR recruitment of methyltransferase enzymes to specific loci (Papoutsis et al., 2015; Alavian-Ghavanini and Rüegg, 2018) may underlie our site-specific findings.
About 71% of differentially methylated genes were hypomethylated in adulthood as a result of developmental TCDD exposure; this corresponds with reports of AhR and nuclear receptor-driven site-specific CpG demethylation (Métivier et al., 2008; Amenya et al., 2016) and aligns with our findings of transcriptomic downregulation of DNA and histone methyltransferase genes. Our work demonstrates that 69% of differentially methylated sites were located within 100 kb of linked genes, with about 9% located within 10 kb. Although promoter methylation is commonly associated with differential gene expression, McGaughey et al. (2014) found the pattern of methylation 10 kb upstream or downstream of a gene to be more predictive of gene expression than promoter methylation status, suggesting a broadened, more complex view of the interplay between methylome and transcriptome.
Despite the overall similarities in affected transcriptomic and methylomic functional pathways in TCDD-exposed fish, seven genes are both differentially expressed and differentially methylated. One gene involved in neural crest development, crestin, is linked to both individual CpGs and DMRs, while the GABA receptor subunit gabra1 and neuronal transcription factor sox21b are associated specifically with DMRs. The remaining four genes are associated only with individual CpG sites: histone-modifying genes wdr5 and riox1 as mentioned above, along with the endosomal GTPase rab11bb and flrt2, involved in embryonic development. This is not wholly unexpected considering prior literature, as other studies exposing zebrafish embryos to TCDD resulted in limited genome-wide alterations in methylation and did not demonstrate an inverse relationship between gene expression and promoter methylation (Olsvik et al., 2014; Aluru et al., 2015). It is possible that other epigenetic or regulatory effects are controlling this interaction. This disparity also suggests a tightly regulated temporal aspect to dysregulation of reproductive pathways: exposure to TCDD during critical periods alters the methylome, which shifts the careful interplay of rapid gene activation and repression during early life (Nagy and Turecki, 2012; Piferrer, 2013). Although such transcriptomic effects may not linger into adulthood, there may be lifelong consequences of disrupting development-specific gene programming, including changes to gonadal histology that present in adult fish or shifts in the sensitivity and response of the hypothalamic–pituitary–gonadal axis. Another factor that may contribute to the minimal overlap between differentially expressed and differentially methylated genes is the assessment of whole testes. As testicular tissue is a heterogeneous mixture of germ cells at various stages of maturity with supporting Leydig and Sertoli cells, taking sperm samples or sorting testicular cell populations prior to analysis may help prevent any potential masking of transcriptomic and methylomic effects. Future directions include single cell analysis and further investigations of these potential relationships by qPCR and pyrosequencing.
As exposure to TCDD occurred during juvenile development, detection of differential methylation in reproductive pathways well into adulthood, taken together with the transgenerational infertility phenotype, suggests persistent epimutations that may be inherited in the F1 or F2 generations. DNA methylation, although a potential mechanism for the phenotypic and transcriptomic effects previously demonstrated, is only one of several epigenetic factors, which also include histone modification and microRNA, that cooperate to regulate the transcriptome (Sun et al., 2013; Alavian-Ghavanini and Rüegg, 2018). Differential methylation and transcription of several DNA methylation, histone methylation, and miRNA genes suggests that multiple epigenetic mechanisms interact with each other and with the transcriptome to mediate toxicant-induced infertility in zebrafish, establishing the need for future studies to investigate these mechanisms and determine promising methods of clinical intervention for adult-onset and transgenerational reproductive disease.
CA and DM wrote the manuscript, generated the tables, and performed the IPA analysis. TB edited the manuscript, performed the experimental exposures, dissections, tissue collection, and analyzed the microarray data. KG analyzed all WGBS data and generated the corresponding figures.
This work was supported by the National Institute of Environmental Health Sciences (T32 ES007015), the National Center for Advancement of Translational Sciences (K01 OD01462 to TB), the WSU Center for Urban Responses to Environmental Stressors (P30 ES020957 to DM and TB), and the National Science Foundation (Grant No. 1735038 to CA). Funding was also provided by the NIH Center grant (P30 CA022453 to KG) to the Karmanos Cancer Institute at Wayne State University and the Perinatology Research Branch of the National Institutes of Child Health and Development at Wayne State University.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We would like to acknowledge the Wayne State University Applied Genomics Technology Center for providing WGBS services and the use of Ingenuity Pathway Analysis software, and the University of Wisconsin Biotechnology Gene Expression Center for providing Affymetrix GeneChip services. We are grateful to Emily Crofts, Jeremy Shields, Bridget Baker, Marcelino Akemann, Anna Boegehold, and the other members of the Baker Laboratory for the time and effort they have dedicated to fish husbandry, laboratory maintenance, and their general advice and support on this manuscript.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2018.00719/full#supplementary-material
TABLE S1 | Concentrations, A260/A230, and A260/A280 of DNA after tissue isolation, bisulfite conversion, and amplification.
TABLE S2 | Mapping efficiency, total C analyzed, percent C methylated, total CpG methylated and percent C methylated in CpG context for control and TCDD-treated DNA.
TABLE S3 | Data for each differentially methylated site including genomic location, p-value, q-value, methylation difference, distance to nearest gene, and location and annotation of nearest gene.
TABLE S4 | Diseases and biological functions results from Ingenuity Pathway Analysis (IPA) of genes with significant changes in DNA methylation between control and TCDD-exposed groups.
Akalin, A., Kormaksson, M., Li, S., Garrett-Bakelman, F. E., Figueroa, M. E., Melnick, A., et al. (2012). MethylKit: a comprehensive R package for the analysis of genome-wide DNA methylation profiles. Genome Biol. 13:R87. doi: 10.1186/gb-2012-13-10-r87s
Alavian-Ghavanini, A., and Rüegg, J. (2018). Understanding epigenetic effects of endocrine disrupting chemicals: from mechanisms to novel test methods. Basic Clin. Pharmacol. Toxicol. 122, 38–45. doi: 10.1111/bcpt.12878
Alpatov, R., Lesch, B. J., Nakamoto-Kinoshita, M., Blanco, A., Chen, S., Stützer, A., et al. (2014). A chromatin-dependent role of the fragile X mental retardation protein FMRP in the DNA damage response. Cell 157, 869–881. doi: 10.1016/j.cell.2014.03.040
Aluru, N., Kuo, E., Helfrich, L. W., Karchner, S. I., Linney, E. A., Pais, J. E., et al. (2015). Developmental exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin alters DNA methyltransferase (dnmt) expression in zebrafish (Danio rerio). Toxicol. Appl. Pharmacol. 284, 142–151. doi: 10.1016/j.taap.2015.02.016
Alves, M. G., Martins, A. D., Rato, L., Moreira, P. I., Socorro, S., and Oliveira, P. F. (2013). Molecular mechanisms beyond glucose transport in diabetes-related male infertility. Biochem. Biophys. Acta Mol. Basis Dis. 1832, 626–635. doi: 10.1016/j.bbadis.2013.01.011
Aravind, L., Abhiman, S., and Iyer, L. M. (2011). Natural history of the eukaryotic chromatin protein methylation system. Prog. Mol. Biol. Transl. Sci. 101, 105–176. doi: 10.1016/B978-0-12-387685-0.00004-4
Bächner, D., Manca, A., Stelnbach, P., Wöhrle, D., Just, W., Vogel, W., et al. (1993). Enhanced expression of the murine FMR1 gene during germ cell proliferation suggests a special function in both the male and the female gonad. Hum. Mol. Genet. 12, 2043–2050. doi: 10.1093/hmg/2.12.2043
Baker, B. B., Yee, J. S., Meyer, D. N., Yang, D., and Baker, T. R. (2016). Histological and transcriptomic changes in male zebrafish testes due to early life exposure to low level 2,3,7,8-tetrachlorodibenzo-p-dioxin. Zebrafish 13, 413–423. doi: 10.1089/zeb.2016.1275
Baker, T. R., King-Heiden, T. C., Peterson, R. E., and Heideman, W. (2014a). Dioxin induction of transgenerational inheritance of disease in zebrafish. Mol. Cell. Endocrinol. 398, 36–41. doi: 10.1016/j.mce.2014.08.011
Barski, A., Cuddapah, S., Cui, K., Roh, T. Y., Schones, D. E., Wang, Z., et al. (2007). High-resolution profiling of histone methylations in the human genome. Cell 129, 823–837. doi: 10.1016/j.cell.2007.05.009
Bhat, M. A., Guru, S. A., Mir, R., Waza, A. A., Zuberi, M., Sumi, M. P., et al. (2018). Association of GABAA receptor gene with epilepsy syndromes. J. Mol. Neurosci. 65, 141–153. doi: 10.1007/s12031-018-1081-7
Bollati, V., Baccarelli, A., Hou, L., Bonzini, M., Fustinoni, S., Cavallo, D., et al. (2007). Changes in DNA methylation patterns in subjects exposed to low-dose benzene. Cancer Res. 67, 876–880. doi: 10.1158/0008-5472.CAN-06-2995
Booy, E. P., Howard, R., Marushchak, O., Ariyo, E. O., Meier, M., Novakowski, S. K., et al. (2014). The RNA helicase RHAU (DHX36) suppresses expression of the transcription factor PITX1. Nucleic Acids Res. 42, 3346–3361. doi: 10.1093/nar/gkt1340
Bray, C., Son, J. H., Kumar, P., and Meizel, S. (2005). Mice deficient in CHRNA7, a subunit of the nicotinic acetylcholine receptor, produce sperm with impaired motility1. Biol. Reprod. 73, 807–814. doi: 10.1095/biolreprod.105.042184
Byun, H., Motta, V., Panni, T., Bertazzi, P. A., Apostoli, P., Hou, L., et al. (2013). Evolutionary age of repetitive element subfamilies and sensitivity of DNA methylation to airborne pollutants. Part. Fibre Toxicol. 10:28. doi: 10.1186/1743-8977-10-28
Campos, C., Valente, L. M., and Fernandes, J. M. (2012). Molecular evolution of zebrafish DNMT3 genes and thermal plasticity of their expression during embryonic development. Gene 500, 93–100. doi: 10.1016/j.gene.2012.03.041
Chi, M. N., Auriol, J., Jegou, B., Kontoyiannis, D. L., Turner, J. M., de Rooij, D. G., et al. (2011). The RNA-binding protein ELAVL1/HuR is essential for mouse spermatogenesis, acting both at meiotic and postmeiotic stages. Mol. Biol. Cell 22, 2875–2885. doi: 10.1091/mbc.E11-03-0212
Christie, D. A., Lemke, C. D., Elias, I. M., Chau, L. A., Kirchhof, M. G., Li, B., et al. (2011). Stomatin-like protein 2 binds cardiolipin and regulates mitochondrial biogenesis and function. Mol. Cell. Biol. 31, 3845–3856. doi: 10.1128/MCB.05393-11
Clark, B. S., Winter, M., Cohen, A. R., and Link, B. A. (2011). Generation of Rab-based transgenic lines for in vivo studies of endosome biology in zebrafish. Dev. Dyn. 240, 2452–2465. doi: 10.1002/dvdy.22758
Dayeh, T., Volkov, P., Salö, S., Hall, E., Nilsson, E., Olsson, A. H., et al. (2014). Genome-wide DNA methylation analysis of human pancreatic islets from type 2 diabetic and non-diabetic donors identifies candidate genes that influence insulin secretion. PLoS Genet. 10:e1004160. doi: 10.1371/journal.pgen.1004160
Diem, M. D., Chan, C. C., Younis, I., and Dreyfuss, G. (2007). PYM binds the cytoplasmic exon-junction complex and ribosomes to enhance translation of spliced mRNAs. Nat. Struct. Mol. Biol. 14, 1173–1179. doi: 10.1038/nsmb1321
Ding, G. L., Liu, Y., Liu, M. E., Pan, J. X., Guo, M. X., Sheng, J. Z., et al. (2015). The effects of diabetes on male fertility and epigenetic regulation during spermatogenesis. Asian J. Androl. 17, 948–953. doi: 10.4103/1008-682X.150844
Dolinoy, D. C., Huang, D., and Jirtle, R. L. (2007). Maternal nutrient supplementation counteracts bisphenol A-induced DNA hypomethylation in early development. Proc. Natl. Acad. Sci. U.S.A. 104, 13056–13061. doi: 10.1073/pnas.0703739104
Du, X., Rosenfield, R. L., and Qin, K. (2009). KLF15 is a transcriptional regulator of the human 17β-hydroxysteroid dehydrogenase type 5 gene. A potential link between regulation of testosterone production and fat stores in women. J. Clin. Endocrinol. Metab. 94, 2594–2601. doi: 10.1210/jc.2009-0139
Eskenazi, B., Mocarelli, P., Warner, M., Samuels, S., Vercellini, P., Olive, D., et al. (2002). Serum dioxin concentrations and endometriosis: a cohort study in Seveso, Italy. Environ. Health Perspect. 110, 629–634. doi: 10.1289/ehp.02110629
Ferder, I., Parborell, F., Sundblad, V., Chiauzzi, V., Gómez, K., Charreau, E. H., et al. (2013). Expression of fragile X mental retardation protein and FMR1 mRNA during folliculogenesis in the rat. Reproduction 145, 335–343. doi: 10.1530/REP-12-0305
Fontana, R., and Della Torre, S. (2016). The deep correlation between energy metabolism and reproduction: a view on the effects of nutrition for women fertility. Nutrients 8:87. doi: 10.3390/nu8020087
Ge, R. S., Dong, Q., Sottas, C. M., Chen, H., Zirkin, B. R., and Hardy, M. P. (2005). Gene expression in rat Leydig cells during development from the progenitor to adult stage: a cluster analysis. Biol. Reprod. 72, 1405–1415. doi: 10.1095/biolreprod.104.037499
Gray, S., Feinberg, M. W., Hull, S., Kuo, C. T., Watanabe, M., Sen, S., et al. (2002). The Krüppel-like factor KLF15 regulates the insulin-sensitive glucose transporter GLUT4. J. Biol. Chem. 277, 34322–34328. doi: 10.1074/jbc.M201304200
Hamra, F. K., Schultz, N., Chapman, K. M., Grellhesl, D. M., Cronkhite, J. T., Hammer, R. E., et al. (2004). Defining the spermatogonial stem cell. Dev. Biol. 269, 393–410. doi: 10.1016/j.ydbio.2004.01.027
Jiang, L., Zhang, J., Wang, J. J., Wang, L., Zhang, L., Li, G., et al. (2014). Sperm, but not oocyte, DNA methylome is inherited by zebrafish early embryos. Cell 153, 773–784. doi: 10.1016/j.cell.2013.04.041
King Heiden, T. C., Spitsbergen, J., Heideman, W., and Peterson, R. E. (2009). Persistent adverse effects on health and reproduction caused by exposure of zebrafish to 2,3,7,8-Tetrachlorodibenzo-p-dioxin during early development and gonad differentiation. Toxicol. Sci. 109, 75–87. doi: 10.1093/toxsci/kfp048
Langmead, B., Trapnell, C., Pop, M., and Salzberg, S. L. (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10:R25. doi: 10.1186/gb-2009-10-3-r25
Loveland, K. L., Bakker, M., Meehan, T., Christy, E., von Schönfeldt, V., Drummond, A., et al. (2003). Expression of Bambi is widespread in juvenile and adult rat tissues and is regulated in male germ cells. Endocrinology 144, 4180–4186. doi: 10.1210/en.2002-0124
Luo, R., An, M. I. N., Arduini, B. L., and Henion, P. D. (2001). Specific pan-neural crest expression of zebrafish crestin. Dev. Dyn. 174, 169–174. doi: 10.1002/1097-0177(2000)9999:9999<::AID-DVDY1097>3.0.CO;2-1
Ma, J., Chen, X., Liu, Y., Xie, Q., Sun, Y., Chen, J., et al. (2015). Ancestral TCDD exposure promotes epigenetic transgenerational inheritance of imprinted gene Igf2: methylation status and DNMTs. Toxicol. Appl. Pharmacol. 289, 193–202. doi: 10.1016/j.taap.2015.09.024
Manikkam, M., Tracey, R., Guerrero-Bosagna, C., and Skinner, M. K. (2012). Dioxin (TCDD) induces epigenetic transgenerational inheritance of adult onset disease and sperm epimutations. PLoS One 7:e46249. doi: 10.1371/journal.pone.0046249
Meccariello, R., Chianese, R., Ciaramella, V., Fasano, S., and Pierantoni, R. (2014). Molecular chaperones, cochaperones, and ubiquitination/deubiquitination system: involvement in the production of high quality spermatozoa. Biomed Res. Int. 2014:561426. doi: 10.1155/2014/561426
Métivier, R., Gallais, R., Tiffoche, C., Le Péron, C., Jurkowska, R. Z., Carmouche, R. P., et al. (2008). Cyclical DNA methylation of a transcriptionally active promoter. Nature 452, 45–50. doi: 10.1038/nature06544
Meyer, D. N., Baker, B. B., and Baker, T. R. (2018). Ancestral TCDD exposure induces multigenerational histologic and transcriptomic alterations in gonads of male zebrafish. Toxicol. Sci. 164, 603–612. doi: 10.1093/toxsci/kfy115
Mocarelli, P., Gerthoux, P. M., Patterson, D. G., Milani, S., Limonta, G., Bertona, M., et al. (2008). Dioxin exposure, from infancy through puberty, produces endocrine disruption and affects human semen quality. Environ. Health Perspect. 116, 70–77. doi: 10.1289/ehp.10399
Müller, P., Schulz, R., Maretto, S., Costello, I., Srinivas, S., Bikoff, E., et al. (2011). The fibronectin leucine-rich repeat transmembrane protein Flrt2 is required in the epicardium to promote heart morphogenesis. Development 138, 1297–1308. doi: 10.1242/dev.059386
Nakada, K., Sato, A., Yoshida, K., Morita, T., Tanaka, H., Inoue, S. I., et al. (2006). Mitochondria-related male infertility. Proc. Natl. Acad. Sci. U.S.A. 103, 15148–15153. doi: 10.1073/pnas.0604641103
Olsvik, P. A., Williams, T. D., Tung, H. S., Mirbahai, L., Sanden, M., Skjaerven, K. H., et al. (2014). Impacts of TCDD and MeHg on DNA methylation in zebrafish (Danio rerio) across two generations. Comp. Biochem. Physiol. Part C Toxicol. Pharmacol. 165, 17–27. doi: 10.1016/j.cbpc.2014.05.004
Onishchenko, N., Karpova, N., Sabri, F., Castrén, E., and Ceccatelli, S. (2008). Long-lasting depression-like behavior and epigenetic changes of BDNF gene expression induced by perinatal exposure to methylmercury. J. Neurochem. 106, 1378–1387. doi: 10.1111/j.1471-4159.2008.05484.x
Osella, M., Riba, A., Testori, A., Corá, D., and Caselle, M. (2014). Interplay of microRNA and epigenetic regulation in the human regulatory network. Front. Genet. 5:345. doi: 10.3389/fgene.2014.00345
Papoutsis, A. J., Selmin, O. I., Borg, J. L., and Romagnolo, D. F. (2015). Gestational exposure to the AhR agonist 2,3,7,8-tetrachlorodibenzo-p-dioxin induces BRCA-1 promoter hypermethylation and reduces BRCA-1 expression in mammary tissue of rat offspring: preventive effects of resveratrol. Mol. Carcinog. 54, 261–269. doi: 10.1002/mc.22095
Pilsner, J. R., Liu, X., Ahsan, H., Ilievski, V., Slavkovich, V., Levy, D., et al. (2007). Genomic methylation of peripheral blood leukocyte DNA: influences of arsenic and folate in Bangladeshi adults. Am. J. Clin. Nutr. 86, 1179–1186. doi: 10.1093/ajcn/86.4.1179
Potok, M. E., Nix, D. A., Parnell, T. J., and Cairns, B. R. (2013). Reprogramming the maternal zebrafish genome after fertilization to match the paternal methylation pattern. Cell 153, 759–772. doi: 10.1016/j.cell.2013.04.030
Raffa, S., Scrofani, C., Valente, S., Micaloni, A., Forte, M., Bianchi, F., et al. (2017). In vitro characterization of mitochondrial function and structure in rat and human cells with a deficiency of the NADH: ubiquinone oxidoreductase NDUFC2 subunit. Hum. Mol. Genet. 26, 4541–4555. doi: 10.1093/hmg/ddx333
Ramalho-Santos, J., Varum, S., Amaral, S., Mota, P. C., Sousa, A. P., and Amaral, A. (2009). Mitochondrial functionality in reproduction: from gonads and gametes to embryos and embryonic stem cells. Hum. Reprod. Update 15, 553–572. doi: 10.1093/humupd/dmp016
Ravnskjaer, K., Hogan, M. F., Lackey, D., Tora, L., Dent, S. Y. R., Olefsky, J., et al. (2013). Glucagon regulates gluconeogenesis through KAT2B- and WDR5-mediated epigenetic effects. J. Clin. Invest. 123, 4318–4328. doi: 10.1172/JCI69035
Ribas, L., Vanezis, K., Imués, M. A., and Piferrer, F. (2017). Treatment with a DNA methyltransferase inhibitor feminizes zebrafish and induces long-term expression changes in the gonads. Epigenetics Chromatin. 10:59. doi: 10.1186/s13072-017-0168-7
Robledo, C. A., Yeung, E., Mendola, P., Sundaram, R., Maisog, J., Sweeney, A. M., et al. (2015). Preconception maternal and paternal exposure to persistent organic pollutants and birth size: the LIFE study. Environ. Health Perspect. 123, 88–94. doi: 10.1289/ehp.1308016
Rodgers, A. B., Morgan, C. P., Bronson, S. L., Revello, S., and Bale, T. L. J. (2013). Paternal stress exposure alters sperm microRNA content and reprograms offspring HPA stress axis regulation. J. Neurosci. 33, 9003–9012. doi: 10.1523/JNEUROSCI.0914-13.2013
Rolland, A. D., Lareyre, J. J., Goupil, A. S., Montfort, J., Ricordel, M. J., Esquerré, D., et al. (2009). Expression profiling of rainbow trout testis development identifies evolutionary conserved genes involved in spermatogenesis. BMC Genomics 10:546. doi: 10.1186/1471-2164-10-546
Sanabria, M., Cucielo, M. S., Guerra, M. T., dos Santos Borges, C., Banzato, T. P., et al. (2016). Sperm quality and fertility in rats after prenatal exposure to low doses of TCDD: a three-generation study. Reprod. Toxicol. 65, 29–38. doi: 10.1016/j.reprotox.2016.06.019
Signorini, S., Gerthoux, P. M., Dassi, C., Cazzaniga, M., Brambilla, P., Vincoli, N., et al. (2000). Environmental exposure to dioxin: the Seveso experience. Andrologia 32, 263–270. doi: 10.1046/j.1439-0272.2000.00394.x
Slegtenhorst-Eegdeman, K. E., de Roojj, D. G., Verhoef-Post, M., van de Kant, H. J., Bakker, C. E., Oostra, B. A., et al. (1998). Macroorchidism in FMR1 knockout mice is caused by increased Sertoli cell proliferation during testicular development. Endocrinology 139, 156–162. doi: 10.1210/endo.139.1.5706
Somm, E., Stouder, C., and Paoloni-Giacobino, A. (2013). Effect of developmental dioxin exposure on methylation and expression of specific imprinted genes in mice. Reprod. Toxicol. 35, 150–155. doi: 10.1016/j.reprotox.2012.10.011
Steenland, K., Piacitelli, L., Deddens, J., and Fingerhut, M. (1999). Cancer, heart disease, and diabetes in workers exposed to 2,3,7,8-tetrachlorodibenzo-p-dioxin. J. Natl. Cancer Inst. 91, 779–786. doi: 10.1093/jnci/91.9.779
Stroud, D. A., Surgenor, E. E., Formosa, L. E., Reljic, B., Frazier, A. E., Dibley, M. G., et al. (2016). Accessory subunits are integral for assembly and function of human mitochondrial complex I. Nature 538, 123–126. doi: 10.1038/nature19754
Sun, X., He, Y., Huang, C., Ma, T. T., and Li, J. (2013). The epigenetic feedback loop between DNA methylation and microRNAs in fibrotic disease with an emphasis on DNA methyltransferases. Cell. Signal. 25, 1870–1876. doi: 10.1016/j.cellsig.2013.05.013
Tao, Y., Wu, M., Zhou, X., Yin, W., Hu, B., Crombrugghe, B., et al. (2013). Structural insights into histone demethylase NO66 in interaction with osteoblast-specific transcription factor osterix and gene repression. J. Biol. Chem. 288, 16430–16437. doi: 10.1074/jbc.M112.446849
Terada, K., Yomogida, K., Imai, T., Kiyonari, H., Takeda, N., Kadomatsu, T., et al. (2005). A type I DnaJ homolog, DjA2, regulates androgen receptor signaling and spermatogenesis. EMBO J. 24, 611–622. doi: 10.1038/sj.emboj.7600549
Tian, H., Cao, Y. X., Zhang, X. S., Liao, W. P., Yi, Y. H., Lian, J., et al. (2013). The targeting and functions of miRNA-383 are mediated by FMRP during spermatogenesis. Cell Death Dis. 4:e617. doi: 10.1038/cddis.2013.138
Warner, M., Mocarelli, P., Brambilla, P., Wesselink, A., Samuels, S., Signorini, S., et al. (2013). Diabetes, metabolic syndrome, and obesity in relation to serum dioxin concentrations: the Seveso Women’s Health Study. Environ. Health Perspect. 121, 906–911. doi: 10.1289/ehp.1206113
Yang, A. S., Doshi, K. D., Choi, S., Mason, J., Mannari, R. K., and Gharybian, V. (2006). DNA methylation changes after 5-Aza-2’-deoxycytidine therapy in patients with leukemia. Cancer Res. 66, 5495–5503. doi: 10.1158/0008-5472
Keywords: 2, 3, 7, 8-tetrachlorodibenzo-p-dioxin, DNA methylation, zebrafish, epigenetics, reproduction, endocrine disrupting compounds
Citation: Akemann C, Meyer DN, Gurdziel K and Baker TR (2019) Developmental Dioxin Exposure Alters the Methylome of Adult Male Zebrafish Gonads. Front. Genet. 9:719. doi: 10.3389/fgene.2018.00719
Received: 23 September 2018; Accepted: 21 December 2018;
Published: 11 January 2019.
Edited by:Jennifer L. Freeman, Purdue University, United States
Reviewed by:Neel Aluru, Woods Hole Oceanographic Institution, United States
Kristine Lee Willett, University of Mississippi, United States
Copyright © 2019 Akemann, Meyer, Gurdziel and Baker. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tracie R. Baker, firstname.lastname@example.org
†These authors have contributed equally to this work