Next Article in Journal
Evidence-Based Overview of Accelerometer-Measured Physical Activity during School Recess: An Updated Systematic Review
Next Article in Special Issue
Early Life Exposure to Aflatoxin B1 in Rats: Alterations in Lipids, Hormones, and DNA Methylation among the Offspring
Previous Article in Journal
Harm Perceptions and Beliefs about Potential Modified Risk Tobacco Products
Previous Article in Special Issue
Prenatal Multivitamin Use and MTHFR Genotype Are Associated with Newborn Cord Blood DNA Methylation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 18
Issue 2
10.3390/ijerph18020577
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sex-Specific Alterations in Cardiac DNA Methylation in Adult Mice by Perinatal Lead Exposure

by
Laurie K. Svoboda
1,*,
Kai Wang
2,
Tamara R. Jones
1,
Justin A. Colacino
1,2,3,
Maureen A. Sartor
2,4 and
Dana C. Dolinoy
1,3
1
Environmental Health Sciences, University of Michigan School of Public Health, Ann Arbor, MI 48109-2029, USA
2
Department of Computational Medicine and Bioinformatics, University of Michigan Medical School, Ann Arbor, MI 48109-2216, USA
3
Nutritional Sciences, University of Michigan School of Public Health, Ann Arbor, MI 48109-2029, USA
4
Department of Biostatistics, University of Michigan School of Public Health, Ann Arbor, MI 48109-2029, USA
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(2), 577; https://doi.org/10.3390/ijerph18020577
Submission received: 22 November 2020 / Revised: 27 December 2020 / Accepted: 4 January 2021 / Published: 12 January 2021
(This article belongs to the Special Issue Environmental Epigenetics: Implications for Health and Disease)
Browse Figures
Versions Notes

Abstract

:
Environmental factors play an important role in the etiology of cardiovascular diseases. Cardiovascular diseases exhibit marked sexual dimorphism; however, the sex-specific effects of environmental exposures on cardiac health are incompletely understood. Perinatal and adult exposures to the metal lead (Pb) are linked to several adverse cardiovascular outcomes, but the sex-specific effects of this toxicant on the heart have received little attention. Perinatal environmental exposures can lead to disease through disruption of the normal epigenetic programming that occurs during early development. Using a mouse model of human-relevant perinatal environmental exposure, we investigated the effects of exposure to Pb during gestation and lactation on DNA methylation in the hearts of adult offspring mice (n = 6 per sex). Two weeks prior to mating, dams were assigned to control or Pb acetate (32 ppm) water, and exposure continued until offspring were weaned at three weeks of age. Enhanced reduced-representation bisulfite sequencing was used to measure DNA methylation in the hearts of offspring at five months of age. Although Pb exposure stopped at three weeks of age, we discovered hundreds of differentially methylated cytosines (DMCs) and regions (DMRs) in males and females at five months of age. DMCs/DMRs and their associated genes were sex-specific, with a small, but statistically significant subset overlapping between sexes. Pathway analysis revealed altered methylation of genes important for cardiac and other tissue development in males, and histone demethylation in females. Together, these data demonstrate that perinatal exposure to Pb induces sex-specific changes in cardiac DNA methylation that are present long after cessation of exposure, and highlight the importance of considering sex in environmental epigenetics and mechanistic toxicology studies.

1. Introduction

Environmental exposures that occur very early in life can have long-lasting influences on cardiovascular disease (CVD) risk [1], consistent with the Developmental Origins of Health and Disease (DOHaD) hypothesis. In spite of this, the underlying mechanisms by which they do so are incompletely understood. Developmental environments may influence CVD risk by altering the structure of the heart, including the final number of cells in the heart or density of blood vessels [2]. Exposures may also alter the expression of genes and their protein products critical for normal cardiac function, such as those of the mitochondrial electron transport chain [2]. There are considerable sex differences in the incidence and pathogenesis of CVDs [3]. For example, although men and women are both prone to ischemic heart disease, the pathogenesis differs between sexes, with males more frequently exhibiting obstructive coronary artery disease [3]. In contrast, ischemia in women is more often due to non-obstructive coronary artery disease or microvascular dysfunction. Likewise, women are more likely than men to experience lethal arrhythmias as a result of pharmacologic interventions (antibiotics, antidepressants, antihistamines) or genetic syndromes [3,4]. The molecular basis for these differences is unclear but likely involves genetic, epigenetic, and hormonal factors [3]. How these important sex differences interact with environmental exposures to influence cardiac health has not been investigated. One potential mechanism by which toxicant exposures early in life can lead to CVD is by interfering with the widespread epigenetic and transcriptional programming that drive normal cardiac development [1]. DNA methylation, or the addition of a methyl group to cytosine bases on DNA (5-methylcytosine), plays a critical role in normal development, and disruptions in DNA methylation have been linked to environment-induced disease. The development of the heart is characterized by dynamic changes in DNA methylation [5], and alterations in this epigenetic mark have been identified in a variety of cardiovascular disease states [1,5,6], as well as with toxicant exposures [6,7,8]. Notably, sex differences in epigenetic modifications and transcriptional profiles are present in the embryonic stem cell stage and persist throughout cardiac differentiation and into adulthood [9,10]. Thus, sex-specific alterations in DNA methylation by environmental exposures may have important implications for sexually dimorphic cardiovascular health.
Recent incidents of Pb contamination in municipal water systems underscore the ongoing threat that Pb poses to human health in the US and worldwide. In addition to drinking water, common sources of Pb exposure in the US include household dust from Pb-based paint, imported consumer products, and industrial exposures [11]. Worldwide, e-waste recycling, traditional medicines, industrial emissions, and glazed ceramics comprise additional sources of exposure [11]. It is estimated that 815 million children worldwide have blood lead concentrations of at least 5 μg/dL, a level at which adverse neurodevelopmental, psychiatric and neurological outcomes are well-documented [11,12]. In addition to the established effects on the nervous system, the contribution of Pb to CVD mortality in the US is far greater than previously thought [13]. Pb exposure is linked to high blood pressure, myocardial infarction, stroke, and cardiac arrhythmias in humans and animals [13,14,15,16,17,18]. In a rodent model, neonatal Pb exposure leads to enhanced sensitivity to the arrhythmogenic effects of norepinephrine in adulthood, suggesting that early development is a critical window of susceptibility to the effects of Pb on cardiac health [19,20]. In spite of established links between Pb and CVD, sex differences in these outcomes, and the underlying molecular mechanisms, are unclear. Early life exposure to Pb plays a clear role in reprogramming of DNA methylation in several non-cardiac tissues, and sex-specific effects of perinatal Pb exposure on DNA methylation have been reported [21,22]. However, the sex-specific effects of developmental Pb exposure on the epigenome of the heart, and the implications this may have for disease, are unknown. In this work, we hypothesized that Pb exposure during this critical window of development would also affect cardiac DNA methylation. To test this hypothesis, we utilized a mouse model of perinatal environmental exposure to investigate the effects of an environmentally relevant dose of Pb on DNA methylation in the hearts of adult male and female mice.

2. Materials and Methods

2.1. Animal Exposure Paradigm

This work was conducted as part of a larger study under the National Institute of Environmental Health Sciences (NIEHS) Toxicant Exposures and Responses by Genomic and Epigenomic Regulators of Transcription (TaRGET II) Consortium [23,24]. Mice utilized for this study have recently been described [8]. Procedures for Pb preparation and exposure were conducted exactly as outlined previously [24]. Briefly, virgin a/a females (6–8 weeks old) were mated with virgin a/a males (7–9 weeks old), and randomly assigned to receive control or Pb through drinking water. Pb-acetate was mixed with water to result in a Pb concentration of 32 ppm in drinking water, which results in a human-relevant maternal exposure in the 16–60 μg/dL range [22,24]. Pb-supplemented water was made by dissolving Pb (II) acetate trihydrate (Sigma-Aldrich) in a single batch of distilled water, and Pb concentrations were verified using inductively coupled plasma mass spectrometry with a limit of detection of 1.0 µg/L (ICPMS; NSF International, Ann Arbor, MI, USA). Animals were maintained on a phytoestrogen-free modified AIN-93G diet (TD.95092, 7% Corn Oil Diet, (Envigo, Indianapolis, IN USA). Dams were exposed to either control or Pb-supplemented drinking water for two weeks prior to mating, and exposure continued during gestation and lactation. After weaning on postnatal day 21, the resulting pups were weighed and switched to Pb-free drinking water (Figure 1). Approximately 1–2 male and 1–2 female offspring per litter were followed until 5 months of age (n = 6 animals per sex/exposure- control female, control male, Pb female, Pb male). All animals had access to food and drinking water ad libitum throughout the experiment, remained on a 12-h light/dark cycle, and were housed in polycarbonate-free cages. Health checks were carried out daily by lab personnel and the University of Michigan Unit for Laboratory Animal Medicine (ULAM, Ann Arbor, MI, USA). This study protocol was approved by the University of Michigan Institutional Animal Care and Use Committee (IACUC), protocol # PRO00009800.

2.2. Euthanasia and Tissue Collection

Each mouse was weighed on a weekly basis (Mettler Toledo, Columbus, OH, USA) and given regular health checks. Health checks consisted of a general assessment of appearance (fur coat intact and well-groomed, eyes clear, no signs of fight wounds) and behavior (mobility, nest building, etc.). In addition to checks by lab personnel, a designated animal handler from ULAM checked these cues on a daily basis, and a veterinarian assessed the health status of the mice at least once a week. Animals were euthanized at 5 months of age according to protocols established by the TaRGET II Consortium (Figure 1 and [23]). Euthanasia and tissue collection procedures were recently described [8]. Heart samples were immediately snap-frozen in liquid nitrogen and stored at −80 °C until DNA and RNA extraction.

2.3. DNA Extraction and Enhanced Reduced Representation Bisulfite Sequencing

DNA extraction (1–2 male and 1–2 female mice per litter) and sequencing were performed as outlined previously [8]. Enhanced reduced representation bisulfite (ERRBS) was performed at the University of Michigan Epigenomics and Advanced Genomics Cores as described previously [24,25], and each sample met the quality control criteria for next generation sequencing. Bisulfite conversion efficiencies for all samples exceeded 99.8% (Supplementary Table S1). The average mapping efficiency was 58.9% (Supplementary Table S1). Single end, 50 nucleotide sequencing was performed on a HiSeq4000 platform (Illumina, San Diego, CA, USA). Libraries were multiplexed and sequenced over 2 lanes. Library sizes (with adapters) ranged from 200–400 bp, and the average sequencing depth was >118 million reads per sample. On average, this method captured 4.8% of genomic CpGs.

2.4. Bioinformatics Pipeline, Quality Control, and Differential Methylation Analysis

DNA methylation analysis, including quality control, trimming, alignment, and methylation calling were conducted exactly as outlined previously [8]. Briefly, we removed CpGs with read coverage >1000 or <10. Opposite strand CpGs at the same position were combined via destranding. Sex chromosomes were included in this analysis. We performed differential methylation testing on individual CpG sites (DMCs), and differentially methylated regions (DMRs) were identified in 1000 bp tiles using the same process. In order to be included in the analysis, sufficient sequencing coverage for a minimum of 4 samples from the Pb group and 4 samples from the control group was required. Differentially methylated CpGs and regions were identified exactly as recently outlined [8]. Run was included as a covariate in the model to adjust for batch effects. After obtaining p-values, we adjusted for multiple testing using the FDR approach. Sites and regions with FDR < 0.05 and an absolute difference in methylation of >10% were considered significant. To determine the distribution of differentially methylated sites across the genome, we used the annotatR R Bioconductor package (v1.5.9) to annotate the CpGs to the mouse mm10 genome ([26] and methods from reference [8]). To determine whether the proportion of differentially methylated cytosines (DMCs) falling into each annotation was significantly different from the total regions tested, we conducted a Chi Square test [8].

2.5. Pathway Analysis of Differentially Methylated Regions (DMRs)

Poly-enrich analysis was conducted using all DMRs with a p-value of at least 0.1, with the following settings: genesets = GOBP, GOCC, and GOMF; locusdef = 1 kb, min_geneset_size = 15, max_geneset_size = 2000. For GREAT [27] analysis, BED files of DMRs for each sex were uploaded to the GREAT web interface [27] using the mouse mm10 species assembly. Association rule setting “basal plus extension” was utilized with the following parameters: Proximal: 5 kb upstream, 1 kb downstream, plus Distal: up to 1000 kb. STRING network analysis, genes associated with DMRs were analyzed using default parameters: full STRING network and a required interaction score of 0.4.

2.6. Gene Expression Analysis

RNA-seq library preparation and sequencing were performed at the University of Michigan Advanced Genomics Core (N = 6 animals per sex, per condition). Library preparation was carried out using the KAPA mRNA Hyper Prep Kit (Roche, Wilmington, MA, USA) with Dual Indexing Adapters following manufacturer instructions. Quantity and quality of the prepared libraries were confirmed with the Agilent 2200 TapeStation (Agilent, Santa Clara, CA, USA). Sequencing of paired-end 50 base pair reads was carried out on the Illumina NovaSeq 6000 (Illumina, San Diego, CA, USA) in the S2 flow cell. Sequenced reads were trimmed via Trim Galore [28], and quality control was assessed with FastQC [29]. STAR was used for the alignment step [30]. Trimming, QC and alignment were all carried out with default parameters. Normalized read counts for Pb-exposed vs. control samples were obtained for each gene using the TMM method of edgeR [31], stratifying by sex. Statistical analysis was conducted as noted in the next section.

2.7. Statistical Analysis

For heart weights and gene expression analysis (lead vs. control for each gene of interest in a targeted analysis), animals were stratified by sex, and linear mixed-effects regression was carried out using the lme4 and lmerTest packages in R version 3.6.1 [32]. Litter-specific random effects were included to account for within-litter correlation. Statistical analysis of overlapping genes or sites between groups was conducted using a hypergeometric test [33]. p-values and representation factors are reported, where representation factor = the number of overlapping genes divided by the number of genes expected to overlap by chance.

3. Results

3.1. Litter Parameters and Phenotypic Effects

Exposure to Pb during gestation and lactation did not significantly alter litter size, pup mortality, or the percentage of females in each litter, and animal weights at 5 months of age were not significantly different between control and Pb exposed animals [24]. Pb exposure had no significant effect on relative heart weights in either males or females at 5 months of age (Figure 2). This finding may have been due to a relatively small sample size compared to our previous studies in which we observed phenotypic effects of Pb exposure [22,34].

3.2. Genome-Wide Changes in DNA Methylation with Developmental Pb Exposure

In order to investigate the effects of gestational and lactational Pb exposure on DNA methylation, we utilized ERRBS to measure DNA methylation in isolated whole heart tissue from male and female offspring at 5 months of age. Although lead exposure had ceased months before, at 5 months of age we observed >1000 differentially methylated cytosines (DMCs), and several hundred differentially methylated regions (DMRs) in both males and females (Table 1 and Table 2 and Supplementary Tables S2–S5). The total number of cytosines and regions tested was similar between sexes (Table 1 and Table 2). For DMCs, the absolute magnitude of methylation change was as high as 57% in males and 68% in females (Figure 3). For DMRs, we observed a maximum magnitude of 72% and 38% for males and females, respectively (Supplementary Figure S1). We annotated these regions to the mouse mm10 genome and found that, similar to our observations with DEHP and BPA exposed animals [8,35], the majority of DMCs and DMRs fell within open sea and intronic regions, and were depleted from CpG islands and promoters (Figure 4). The top 10 most hyper and hypomethylated DMCs and DMRs, and the genes associated with them, are shown in Table 3 and Table 4.

3.3. Pathway Analysis

In order to determine the pathways enriched among the DMRs, we first performed analysis using Poly-Enrich [36], which has been shown to have a more accurate false positive rate than other pathway enrichment tests. We stratified our analysis by sex and direction of methylation change. In both males and females, we observed enrichment for pathways relevant to cardiac function, including ion channel activity, but the results were not statistically significant (Supplementary Tables S6–S9). As a parallel approach, we utilized the Genomic Regions Enrichment of Annotations (GREAT) tool [27]. Among males, DMRs were enriched for several pathways important for normal heart development and function, including the Notch and hedgehog signaling pathways (smoothened is a component of the hedgehog pathway) [37,38], as well as regulation of cardiac muscle hypertrophy (Figure 5A). Among females, DMRs were enriched for pathways associated with histone demethylation, arginine hydroxylation, and body morphogenesis (Figure 5B).
In order to further understand whether DMRs interacted within common biological networks, we utilized STRING network analysis. We stratified the data by sex and direction of methylation and conducted separate analyses for each. In both males and females, there were fewer than 100 differentially hypomethylated genes, so we included all of them in the analysis (35 hypo DMRs in females and 67 hypo DMRs in males). For differentially hypermethylated regions in both sexes, we included the top 100 genes with the largest changes in methylation. In females, the interactions among DMR-associated genes were not statistically significant for hyper or hypomethylated genes. In males, however, we identified significantly more interactions than would be expected by chance for both hypo and hypermethylated DMRs (enrichment p-values p = 0.003 and p = 0.009 respectively, Figure 6). Many of the genes in these networks are important for normal heart development and function, or are associated with disease. Among hypermethylated regions, they included Afap1, Prkce, Atg5, and Tmod1. Among hypomethylated regions, they included Bcas3, Cux1, and Hnrnpu. Consistent with GREAT analysis, Atg5 [39,40] and Cux1 [41,42] interact with the Notch and Hedgehog signaling pathways during normal development and in cancer.

3.4. Overlap between Sexes

We next investigated the sex specificity of DNA methylation after developmental Pb exposure. To this end, we first determined whether any DMCs or DMRs overlapped directly between males and females. Consistent with sex-specific effects on DNA methylation, only two regions were found in common between males and females among both DMCs and DMRs (Supplementary Tables S10 and S11). The sex-specific overlap of DMCs was not statistically significant (p = 0.43, representation factor = 1.4), and the overlap of DMRs approached statistical significance (p = 0.05, representation factor = 5.4) as calculated by hypergeometric test. We then compared the genes mapping to DMCs and DMRs between males and females to further interrogate sex specificity. In males, DMCs and DMRs mapped to 899 and 171 genes, respectively (Supplementary Tables S2 and S4). In females, we identified 753 and 147 genes associated with DMCs and DMRs, respectively (Supplementary Tables S3 and S5). The vast majority of DMCs and DMRs were sex-specific, with a small subset of genes overlapping between sexes (Figure 7A,B and Supplementary Tables S12 and S13). Although a minority of genes were found to be in common between sexes, the overlaps for DMCs and DMRs were statistically significant (p = 3.4 × 10−30, representation factor = 3.7 and 8.1 × 10−6, representation factor = 8.0, respectively, hypergeometric test). We then determined whether the overlapping genes might represent sex-independent biomarkers of Pb exposure, focusing on the 8 genes in common among DMRs. Of these 8 genes, we identified 3 that were associated with cardiovascular diseases, including Rbfox1, Galnt2, and Pi16 [43,44,45]. Altered DNA methylation at each gene occurred in distinct locations based on sex (Table 5 and Supplementary Table S13). Changes in DNA methylation at Galnt2, and Pi16 occurred in the same direction in both sexes (Table 5). Interrogation of the sex-specific DMRs also revealed cardiovascular disease-relevant genes. Among males, these included Atg5, Tmod1, Smad6, Slc26a6, Prkce, Ank2, Cux1, and Lamp2 [46,47,48,49,50,51,52,53]. Among females, they included Tgfb2, Rbfox2, Vdr, Tlr4, Timp3, Pde4b, Akap1, and Grk5 [54,55,56,57,58,59,60,61].

3.5. Gene Expression Analysis

We then investigated whether the changes in DNA methylation at sex-independent and sex-dependent DMRs were associated with altered gene expression. To this end, we first interrogated RNA-seq read count data for each of the three sex-independent genes from control and Pb-exposed males and females. RNA utilized in the RNA-seq analysis was extracted from the same tissue samples utilized for ERRBS. Gene expression data for Rbfox1, Galnt2, and Pi16 are depicted in Figure 7C–E. Rbfox1 exhibited trends toward increased and decreased expression in males and females, respectively that did not reach statistical significance (Figure 7C). For Galnt2, we observed a significant increase in expression in females but not males (Figure 7D). No significant changes were observed in expression of Pi16 (Figure 7E). Among cardiovascular disease-relevant, differentially methylated genes in males, Atg5, Ank2, Cux1 and Lamp2 exhibited significant increases in gene expression. Expression of Tmod1, Smad6, Slc26a6 and Prkce were not significantly different between Pb and control (Supplementary Figure S2A). In females, expression of Akap1 was significantly increased and Grk5 was significantly decreased. The remaining genes showed no significant changes in expression with Pb exposure (Supplementary Figure S2B).

4. Discussion

Although there are clear sex differences in cardiac physiology and pathophysiology [3], the sex-specific effects of environmental exposures on cardiovascular health are poorly understood. In this work, we demonstrate that Pb exposure during gestation and lactation leads to changes in DNA methylation in the heart that are present in adulthood. The levels of Pb used in this study result in maternal blood Pb levels comparable to those currently observed in women of child-bearing age in poorer countries [62]. Although blood Pb levels have fallen significantly in the US since the 1970s–1980s, levels within this range were not uncommon during the mid-late 20th century [63]. Importantly, children exposed to higher levels of Pb during pregnancy in the mid-late 20th century are now at an age in which CVDs are of significant concern. Our findings are novel for several reasons. To our knowledge, this is the first report demonstrating that Pb exposure during early development leads to altered cardiac epigenetic programming that is present in adulthood. Moreover, sex-specific effects of Pb exposure on the cardiac epigenome have not yet been investigated. Indeed, although the effects of Pb on the nervous system are well-established, this work contributes to a growing body of evidence highlighting potential adverse effects of this metal on the cardiovascular system.

4.1. Genome-Wide Changes in DNA Methylation with Developmental Pb Exposure

In this work, we discovered that exposure to Pb during gestation and lactation resulted in genome-wide changes in cardiac tissue DNA methylation. In both males and females, DMCs and DMRs were enriched in introns and intergenic regions of the genome, consistent with what we observed previously for DEHP [8] and BPA [35] exposures. Intergenic and intronic regions of the genome harbor enhancers and other regulatory elements critical for the regulation of tissue-specific gene expression [64]. Likewise, dynamic DNA methylation at these elements is critically important for regulation of normal tissue differentiation, and is de-regulated in the context of diseases such as cancer [65]. Consistent with this, GREAT pathway analysis revealed enrichment of pathways associated with development and epigenetic regulation. In males, this included the regulation of the Notch pathway, as well as smoothened signaling, a component of the hedgehog pathway. Both pathways play critical roles in normal cardiac development, as well as in adult cardiovascular diseases [37,38,66]. In females, these included pathways associated with lysine (H3K36) demethylation and arginine hydroxylation, epigenetic processes that are important for regulation of chromatin structure and function in the heart and other tissues [67,68]. As a complementary approach, we conducted STRING network analysis and found that, in males, genes associated with the DMRs interacted in networks associated with cardiovascular development and disease, including Afap1, Prkce, Atg5, Tmod1, Bcas3, Cux1, and Hnrnpu [47,50,52,69,70,71,72]. The implications of altered DNA methylation at these development and disease-relevant loci for long-term cardiac health require further investigation.

4.2. Implications of Pb Exposure for Cardiac Function

The effects of Pb-induced changes in DNA methylation on heart function are not yet clear. Pb exists predominately as a divalent cation in the body, and thus interferes with the function of other divalent cations, including calcium. Acute exposure to Pb disrupts cardiac calcium signaling, which is critical for normal heart function, and leads to arrhythmias and impaired cardiac contractility [73,74]. Pb may also exert toxic effects in part through activation of the aryl hydrocarbon receptor [75], and can potentiate the effects of other cardiac toxicants [76]. The effects of developmental Pb exposure on cardiac function, however, are far less clear. In this study, potentially due to the small sample size, we observed no significant changes in heart weights in either males or females, and we did not assess the functional effects of early Pb exposure. Previous work in rats demonstrated that exposure to Pb during the early postnatal period lead to an increased incidence of cardiac arrhythmias in response to norepinephrine in adulthood, and the presence of the toxicant during the early developmental period was necessary for the observed effect [19,20]. These findings suggest that early Pb exposure may affect the sensitivity to additional stressors later in life. Given the extensive epigenetic programming that occurs during this critical window of development, and the observed functional effects during adulthood, it is plausible that developmental Pb exposure may elicit these effects through epigenetic mechanisms. Additional studies designed to investigate this important question are currently underway.

4.3. Potential Mechanisms

The molecular mechanism(s) by which Pb exposure leads to altered cardiac DNA methylation are currently unclear. DNA methylation is mediated by DNA methyltransferases (DNMTs), and active removal of this mark occurs via the activity of TET dioxygenases (TETs). DNMTs require the cofactor S-adenosylmethionine (SAM), which provides the methyl group for cytosine methylation. TETs utilize iron (II), ascorbic acid, and alpha ketoglutarate to catalyze the conversion of 5-methylcytosine to the more oxidized products 5-hydroxymethylcytosine (5hmC), 5-fluorocytosine 5 (fC), and 5-carboxylcytosine (5caC). Conversion of these oxidized species back to unmethylated cytosine occurs through thymine DNA glycosylase (TDG)-mediated base excision repair. Thus, Pb may interfere with DNA methylation by altering the expression of these enzymes and/or the levels of their cofactors. Consistent with the first possibility, developmental Pb exposure leads to reduced activity and expression of DNMTs in vitro and in vivo [77,78]. No effects of Pb exposure on expression of TETs or TDG have been reported; however, Pb does cause alterations in the levels of 5hmC [79]. In support of an effect on cofactors, exogenous SAM mitigates the effects of Pb exposure on acute toxicity and adverse neurodevelopmental outcomes, suggesting that Pb may perturb SAM homeostasis [80,81]. Moreover, Pb interferes with the activity of the TCA cycle enzyme isocitrate dehydrogenase, which generates alpha ketoglutarate [16,82]. Vitamin C has also been shown to attenuate the adverse effects of Pb, although whether this is due to modulation of DNA methylation is unclear [83,84].
The molecular underpinnings of sex-specificity in differential DNA methylation are currently unknown. Sex differences in epigenetic profiles are present during cardiac differentiation and in adulthood, highlighting differential epigenetic regulation of cardiac development and function between males and females [9]. During early development, DNA methylation undergoes two distinct waves of widespread erasure and re-methylation in both males and females; however, the kinetics and underlying mechanisms driving this programming differ between sexes. The first wave occurring post-fertilization results in removal of gamete methylation patterns and establishment of the embryonic methylation program [85]. The second wave occurs in the primordial germ cells of the developing embryo, in which sex-specific methylation patterns are established. In pre-implantation embryos, removal of DNA methylation in females occurs via a slower, passive loss of methylation, while in males, demethylation occurs with faster kinetics via an active mechanism mediated by TETs [85]. In the primordial germ cells, the process of re-methylation is complete before birth in males, but in females, this re-methylation continues until puberty [85]. These differences in the timing of methylation erasure and re-establishment may contribute to the observed sex differences in DNA methylation in response to toxicant exposures. Additional mechanistic studies are necessary to address this question.

4.4. Sex Differences in the Effects of Pb Exposure on the Heart

In this work, we discovered that the effects of gestational and lactational Pb exposure on the heart are sex-specific. Specifically, we observed few DMCs/DMRs that overlapped directly between males and females, and a minority of DMC/DMR-associated genes were found to be in common between sexes. This is consistent with our previous findings for DEHP exposed heart, as well as lead-exposed liver [8,24]. Genome-wide and site-specific changes in cardiac DNA methylation have been reported for other environmental exposures, but studies thus far have largely focused on males [7,86]. Notably, although the epigenetic mechanisms underlying sex differences in CVD are unclear, sexually dimorphic patterns of DNA methylation are associated with various CVDs [87]. Thus, sex-specific perturbations in DNA methylation might differentially affect susceptibility to adverse cardiovascular outcomes after perinatal exposure to Pb. Among the small number of DMRs that overlapped between sexes, we identified several associated with cardiovascular disease, leading us to hypothesize that these genes may be similarly de-regulated by Pb exposure in both sexes. However, Galnt2 was the only sex-independent gene investigated that exhibited significant changes in expression and DNA methylation, and altered expression occurred in females but not males. As the locations of the DMRs within each of the genes differed between males and females, differential effects on gene expression are not surprising. Several sex-specific DMR-associated genes also exhibited significant changes in expression concomitant with differential methylation. Collectively, these data suggest that sex had a profound effect on Pb-induced changes in DNA methylation. As only a minority of studies in cardiovascular epigenetics stratify data by sex [88], in spite of clear sex differences in CVD, our work highlights the necessity of considering sex as a biological variable in cardiovascular environmental health studies.

4.5. Limitations of the Study

Although our findings have important implications for cardiovascular environmental health, there are several limitations to this study. First, because we did not measure the level of Pb in the heart at 5 months of age, we cannot rule out the possibility that there may still be Pb present in the hearts at this time point. However, the majority of Pb is stored in bone, with a relatively small quantity absorbed into soft tissues [89,90]. Moreover, the turnover of Pb is far more rapid in soft tissues compared to bone [89,90]. Although Pb is released from bone during pregnancy and age-related bone loss, the animals examined in this study did not experience pregnancy or old age [91,92]. Thus, in the absence of Pb exposure for over 4 months, we would expect minimal residual Pb in the hearts of adult animals. Whether Pb remains in the body after early developmental exposure or not, our data suggest that Pb exposure during early life causes epigenetic changes that may have important implications for health across the life course. A second limitation to this study is the use of the ERRBS method, which is designed to enrich for CpG-rich regions of the genome. Because of this bias, it is plausible that other regions of biological significance in the genome were altered by Pb exposure but were not detected using this method. In addition, this method uses a traditional bisulfite conversion protocol, which does not permit us to discriminate between 5mC and 5hmC. Thus, some of the changes in DNA methylation may reflect changes in 5hmC rather than 5mC. Importantly, 5hmC is also present in the heart and is important for cardiac development and disease [93,94]. Studies to investigate the effects of Pb exposure on other epigenetic modifications, including 5hmC, are warranted.

5. Conclusions

In conclusion, we have demonstrated, for the first time, that Pb exposure during early development leads to changes in DNA methylation in adulthood that are strongly dependent on sex. This work adds to a growing body of evidence linking the early developmental environment to epigenetic changes in the cardiovascular system. It will be of great interest to determine whether the observed changes in DNA methylation are associated with an increased risk of CVD across the life course.

Supplementary Materials

The following are available online at https://www.mdpi.com/1660-4601/18/2/577/s1, Figure S1: Volcano plots depicting DMRs for Pb exposed compared to control in males (A) and females (B). DMRs in red did not meet the criteria for significance (at least a 10% change in methylation and FDR <0.05). Significantly hypomethylated DMRs are shown in blue, and significantly hypermethylated DMRs are shown in green, Table S1: ERRBS QC data, Table S2: ERRBS QC data, Table S3: DMCs in 5 month female hearts, Table S4: DMRs in 5 month male hearts, Table S5: DMRs in 5 month female hearts, Table S6: Poly-Enrich pathway analysis of hypermethylated DMRs in males, Table S7: Poly-Enrich pathway analysis of hypomethylated DMRs in males, Table S8: Poly-Enrich pathway analysis of hypermethylated DMRs in females, Table S9: Poly-Enrich pathway analysis of hypomethylated DMRs in females, Table S10: Differentially methylated cytosines (DMCs) with direct overlap between male and female hearts, Table S11: Differentially methylated regions (DMRs) with direct overlap between male and female hearts, Table S12: DMC-associated genes that overlap between males and females, Table S13: DMR-associated genes that overlap between males and females.

Author Contributions

All authors contributed to the conceptualization of the study. T.R.J. performed experiments. L.K.S. and K.W. analyzed data. L.K.S., D.C.D., M.A.S. and J.A.C. constructed the study, supervised experiments and data collection, and interpreted the data. L.K.S. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Institute of Environmental Health Sciences (NIEHS) TaRGET II Consortium (ES026697), NIEHS Grant K01 (ES032048), NIEHS Grant R35 (ES031686), the University of Michigan NIEHS/EPA Children’s Environmental Health and Disease Prevention Center P01 (ES022844/RD83543601), the Michigan Lifestage Environmental Exposures and Disease (M-LEEaD) NIEHS Core Center (P30 ES017885), Institutional Training Grant T32 (ES007062), and NIEHS grant R01 (ES028802).

Institutional Review Board Statement

This study protocol was approved by the University of Michigan Institutional Animal Care and Use Committee (IACUC), protocol # PRO00009800.

Informed Consent Statement

Not applicable.

Data Availability Statement

The ERRBS data presented in this study are openly available in GEO, accession number GSE152711. All other data can be found in the Supplementary Excel file.

Acknowledgments

We would like to acknowledge the members of the University of Michigan Epigenomics Core, the Advanced Genomics Core, and the Michigan Lifestage Environmental Exposures and Disease Center (M-LEEaD) Omics and Bioinformatics Core for providing their expertise and facilities for generation and analysis of ERRBS data. We especially thank Claudia Lalancette in the Epigenomics Core for her advice and expertise on sample preparation and data analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Thornburg, K.L. The programming of cardiovascular disease. J. Dev. Orig. Health Dis. 2015, 6, 366–376. [Google Scholar] [CrossRef]
  2. Perrino, C.; Ferdinandy, P.; Botker, H.E.; Brundel, B.; Collins, P.; Davidson, S.M.; den Ruijter, H.M.; Engel, F.B.; Gerdts, E.; Girao, H.; et al. Improving Translational Research in Sex-specific Effects of Comorbidities and Risk Factors in Ischemic Heart Disease and Cardioprotection: Position Paper and Recommendations of the ESC Working Group on Cellular Biology of the Heart. Cardiovasc. Res. 2020, cvaa155. [Google Scholar] [CrossRef]
  3. Regitz-Zagrosek, V.; Kararigas, G. Mechanistic Pathways of Sex Differences in Cardiovascular Disease. Physiol Rev. 2017, 97, 1–37. [Google Scholar] [CrossRef] [Green Version]
  4. Ebert, S.N.; Liu, X.K.; Woosley, R.L. Female gender as a risk factor for drug-induced cardiac arrhythmias: Evaluation of clinical and experimental evidence. J. Womens Health 1998, 7, 547–557. [Google Scholar] [CrossRef]
  5. Gilsbach, R.; Preissl, S.; Gruning, B.A.; Schnick, T.; Burger, L.; Benes, V.; Wurch, A.; Bonisch, U.; Gunther, S.; Backofen, R.; et al. Dynamic DNA methylation orchestrates cardiomyocyte development, maturation and disease. Nat. Commun. 2014, 5, 5288. [Google Scholar] [CrossRef] [Green Version]
  6. Zhou, W.; Tian, D.; He, J.; Yan, X.; Zhao, J.; Yuan, X.; Peng, S. Prolonged exposure to carbon nanoparticles induced methylome remodeling and gene expression in zebrafish heart. J. Appl. Toxicol. 2019, 39, 322–332. [Google Scholar] [CrossRef]
  7. Goodson, J.M.; Weldy, C.S.; MacDonald, J.W.; Liu, Y.; Bammler, T.K.; Chien, W.M.; Chin, M.T. In utero exposure to diesel exhaust particulates is associated with an altered cardiac transcriptional response to transverse aortic constriction and altered DNA methylation. FASEB J. 2017, 31, 4935–4945. [Google Scholar] [CrossRef] [Green Version]
  8. Svoboda, L.K.; Wang, K.; Cavalcante, R.G.; Neier, K.; Colacino, J.A.; Sartor, M.A.; Dolinoy, D.C. Sex-Specific Programming of Cardiac DNA Methylation by Developmental Phthalate Exposure. Epigenet Insights 2020, 13, 2516865720939971. [Google Scholar] [CrossRef] [PubMed]
  9. Deegan, D.F.; Karbalaei, R.; Madzo, J.; Kulathinal, R.J.; Engel, N. The developmental origins of sex-biased expression in cardiac development. Biol. Sex. Differ. 2019, 10, 46. [Google Scholar] [CrossRef] [PubMed]
  10. McCormick, H.; Young, P.E.; Hur, S.S.J.; Booher, K.; Chung, H.; Cropley, J.E.; Giannoulatou, E.; Suter, C.M. Isogenic mice exhibit sexually-dimorphic DNA methylation patterns across multiple tissues. BMC Genom. 2017, 18, 966. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Obeng-Gyasi, E. Sources of lead exposure in various countries. Rev. Environ. Health 2019, 34, 25–34. [Google Scholar] [PubMed]
  12. Burki, T. Report says 815 million children have high blood lead levels. Lancet 2020, 396, 370. [Google Scholar] [CrossRef]
  13. Lanphear, B.P.; Rauch, S.; Auinger, P.; Allen, R.W.; Hornung, R.W. Low-level lead exposure and mortality in US adults: A population-based cohort study. Lancet Public Health 2018, 3, e177–e184. [Google Scholar] [CrossRef]
  14. Kieltucki, J.; Dobrakowski, M.; Pawlas, N.; Sredniawa, B.; Boron, M.; Kasperczyk, S. The analysis of QT interval and repolarization morphology of the heart in chronic exposure to lead. Hum. Exp. Toxicol. 2017, 36, 1081–1086. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, A.; Hu, H.; Sanchez, B.N.; Ettinger, A.S.; Park, S.K.; Cantonwine, D.; Schnaas, L.; Wright, R.O.; Lamadrid-Figueroa, H.; Tellez-Rojo, M.M. Association between prenatal lead exposure and blood pressure in children. Environ. Health Perspect. 2012, 120, 445–450. [Google Scholar] [PubMed]
  16. Basha, D.C.; Basha, S.S.; Reddy, G.R. Lead-induced cardiac and hematological alterations in aging Wistar male rats: Alleviating effects of nutrient metal mixture. Biogerontology 2012, 13, 359–368. [Google Scholar] [PubMed]
  17. Afridi, H.I.; Kazi, T.G.; Kazi, N.; Kandhro, G.A.; Baig, J.A.; Jamali, M.K.; Arain, M.B.; Shah, A.Q.; Shah, F.; Khan, S.; et al. Association of environmental toxic elements in biological samples of myocardial infarction patients at different stages. Biol. Trace Elem. Res. 2011, 141, 26–40. [Google Scholar] [CrossRef]
  18. Navas-Acien, A.; Guallar, E.; Silbergeld, E.K.; Rothenberg, S.J. Lead exposure and cardiovascular disease—A systematic review. Environ. Health Perspect. 2007, 115, 472–482. [Google Scholar] [CrossRef] [Green Version]
  19. Hejtmancik, M., Jr.; Williams, B.J. Time and level of perinatal lead exposure for development of norepinephrine cardiotoxicity. Res. Commun. Chem. Pathol. Pharm. 1979, 24, 367–376. [Google Scholar]
  20. Williams, B.J.; Hejtmancik, M.R., Jr.; Abreu, M. Cardiac effects of lead. Fed. Proc. 1983, 42, 2989–2993. [Google Scholar]
  21. Sobolewski, M.; Varma, G.; Adams, B.; Anderson, D.W.; Schneider, J.S.; Cory-Slechta, D.A. Developmental Lead Exposure and Prenatal Stress Result in Sex-Specific Reprograming of Adult Stress Physiology and Epigenetic Profiles in Brain. Toxicol. Sci. 2018, 163, 478–489. [Google Scholar] [CrossRef] [PubMed]
  22. Faulk, C.; Barks, A.; Liu, K.; Goodrich, J.M.; Dolinoy, D.C. Early-life lead exposure results in dose- and sex-specific effects on weight and epigenetic gene regulation in weanling mice. Epigenomics 2013, 5, 487–500. [Google Scholar] [CrossRef] [Green Version]
  23. Wang, T.; Pehrsson, E.C.; Purushotham, D.; Li, D.; Zhuo, X.; Zhang, B.; Lawson, H.A.; Province, M.A.; Krapp, C.; Lan, Y.; et al. The NIEHS TaRGET II Consortium and environmental epigenomics. Nat. Biotechnol. 2018, 36, 225–227. [Google Scholar] [CrossRef] [PubMed]
  24. Svoboda, L.K.; Neier, K.; Wang, K.; Cavalcante, R.G.; Rygiel, C.A.; Tsai1, Z.; Jones, T.R.; Liu, S.; Goodrich, J.M.; Lalancette, C.; et al. Tissue and sex-specific programming of DNA methylation by perinatal lead exposure: Implications for environmental epigenetics studies. Epigenetics 2020, 1–21. [Google Scholar] [CrossRef] [PubMed]
  25. Garrett-Bakelman, F.E.; Sheridan, C.K.; Kacmarczyk, T.J.; Ishii, J.; Betel, D.; Alonso, A.; Mason, C.E.; Figueroa, M.E.; Melnick, A.M. Enhanced reduced representation bisulfite sequencing for assessment of DNA methylation at base pair resolution. J. Vis. Exp. 2015, 96, e52246. [Google Scholar] [CrossRef] [PubMed]
  26. Cavalcante, R.G.; Sartor, M.A. annotatr: Genomic regions in context. Bioinformatics 2017, 33, 2381–2383. [Google Scholar] [CrossRef]
  27. McLean, C.Y.; Bristor, D.; Hiller, M.; Clarke, S.L.; Schaar, B.T.; Lowe, C.B.; Wenger, A.M.; Bejerano, G. GREAT improves functional interpretation of cis-regulatory regions. Nat. Biotechnol. 2010, 28, 495–501. [Google Scholar] [CrossRef] [Green Version]
  28. Trim Galore. Available online: https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/ (accessed on 19 November 2019).
  29. FastQC: A Quality Control Tool for High Throughput Sequence Data. Available online: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 1 August 2019).
  30. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  31. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [Green Version]
  32. Bates, D.M.M.; Bolker, B.; Walker, S. Fitting Linear Mixed-Effects Models Using lme4. J. Stat. Softw. 2015, 67. [Google Scholar] [CrossRef]
  33. Lund, J. Available online: http://nemates.org/MA/progs/overlap_stats.html (accessed on 8 January 2021).
  34. Faulk, C.; Barks, A.; Sanchez, B.N.; Zhang, Z.; Anderson, O.S.; Peterson, K.E.; Dolinoy, D.C. Perinatal lead (Pb) exposure results in sex-specific effects on food intake, fat, weight, and insulin response across the murine life-course. PLoS ONE 2014, 9, e104273. [Google Scholar] [CrossRef] [PubMed]
  35. Kochmanski, J.J.; Marchlewicz, E.H.; Cavalcante, R.G.; Perera, B.P.U.; Sartor, M.A.; Dolinoy, D.C. Longitudinal Effects of Developmental Bisphenol A Exposure on Epigenome-Wide DNA Hydroxymethylation at Imprinted Loci in Mouse Blood. Environ. Health Perspect. 2018, 126, 077006. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Lee, C.T.; Cavalcante, R.G.; Lee, C.; Qin, T.; Patil, S.; Wang, S.; Tsai, Z.T.Y.; Boyle, A.P.; Sartor, M.A. Poly-Enrich: Count-based methods for gene set enrichment testing with genomic regions. NAR Genom. Bioinform. 2020, 2, lqaa006. [Google Scholar] [CrossRef] [PubMed]
  37. MacGrogan, D.; Munch, J.; de la Pompa, J.L. Notch and interacting signalling pathways in cardiac development, disease, and regeneration. Nat. Rev. Cardiol. 2018, 15, 685–704. [Google Scholar] [CrossRef]
  38. Wiegering, A.; Ruther, U.; Gerhardt, C. The Role of Hedgehog Signalling in the Formation of the Ventricular Septum. J. Dev. Biol. 2017, 5, 17. [Google Scholar] [CrossRef] [Green Version]
  39. Jia, Z.; Wang, J.; Wang, W.; Tian, Y.; XiangWei, W.; Chen, P.; Ma, K.; Zhou, C. Autophagy eliminates cytoplasmic beta-catenin and NICD to promote the cardiac differentiation of P19CL6 cells. Cell Signal. 2014, 26, 2299–2305. [Google Scholar] [CrossRef] [Green Version]
  40. Fan, J.; Zhang, X.; Wang, S.; Chen, W.; Li, Y.; Zeng, X.; Wang, Y.; Luan, J.; Li, L.; Wang, Z.; et al. Regulating autophagy facilitated therapeutic efficacy of the sonic Hedgehog pathway inhibition on lung adenocarcinoma through GLI2 suppression and ROS production. Cell Death Dis. 2019, 10, 626. [Google Scholar] [CrossRef] [Green Version]
  41. Sharma, M.; Fopma, A.; Brantley, J.G.; Vanden Heuvel, G.B. Coexpression of Cux-1 and Notch signaling pathway components during kidney development. Dev. Dyn. 2004, 231, 828–838. [Google Scholar] [CrossRef]
  42. Ellis, T.; Gambardella, L.; Horcher, M.; Tschanz, S.; Capol, J.; Bertram, P.; Jochum, W.; Barrandon, Y.; Busslinger, M. The transcriptional repressor CDP (Cutl1) is essential for epithelial cell differentiation of the lung and the hair follicle. Genes Dev. 2001, 15, 2307–2319. [Google Scholar] [CrossRef] [Green Version]
  43. Gao, C.; Ren, S.; Lee, J.H.; Qiu, J.; Chapski, D.J.; Rau, C.D.; Zhou, Y.; Abdellatif, M.; Nakano, A.; Vondriska, T.M.; et al. RBFox1-mediated RNA splicing regulates cardiac hypertrophy and heart failure. J. Clin. Investig. 2016, 126, 195–206. [Google Scholar] [CrossRef] [Green Version]
  44. Peng, P.; Wang, L.; Yang, X.; Huang, X.; Ba, Y.; Chen, X.; Guo, J.; Lian, J.; Zhou, J. A preliminary study of the relationship between promoter methylation of the ABCG1, GALNT2 and HMGCR genes and coronary heart disease. PLoS ONE 2014, 9, e102265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Frost, R.J.; Engelhardt, S. A secretion trap screen in yeast identifies protease inhibitor 16 as a novel antihypertrophic protein secreted from the heart. Circulation 2007, 116, 1768–1775. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Taneike, M.; Yamaguchi, O.; Nakai, A.; Hikoso, S.; Takeda, T.; Mizote, I.; Oka, T.; Tamai, T.; Oyabu, J.; Murakawa, T.; et al. Inhibition of autophagy in the heart induces age-related cardiomyopathy. Autophagy 2010, 6, 600–606. [Google Scholar] [CrossRef] [PubMed]
  47. Fritz-Six, K.L.; Cox, P.R.; Fischer, R.S.; Xu, B.; Gregorio, C.C.; Zoghbi, H.Y.; Fowler, V.M. Aberrant myofibril assembly in tropomodulin1 null mice leads to aborted heart development and embryonic lethality. J. Cell Biol. 2003, 163, 1033–1044. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Niu, H.M.; Liu, C.L. The aberrant expression of Smad6 and TGF-beta in obesity linked cardiac disease. Eur. Rev. Med. Pharm. Sci. 2017, 21, 138–142. [Google Scholar]
  49. Sirish, P.; Ledford, H.A.; Timofeyev, V.; Thai, P.N.; Ren, L.; Kim, H.J.; Park, S.; Lee, J.H.; Dai, G.; Moshref, M.; et al. Action Potential Shortening and Impairment of Cardiac Function by Ablation of Slc26a6. Circ. Arrhythm. Electrophysiol. 2017, 10. [Google Scholar] [CrossRef]
  50. Scruggs, S.B.; Wang, D.; Ping, P. PRKCE gene encoding protein kinase C-epsilon-Dual roles at sarcomeres and mitochondria in cardiomyocytes. Gene 2016, 590, 90–96. [Google Scholar] [CrossRef] [Green Version]
  51. Roberts, J.D.; Murphy, N.P.; Hamilton, R.M.; Lubbers, E.R.; James, C.A.; Kline, C.F.; Gollob, M.H.; Krahn, A.D.; Sturm, A.C.; Musa, H.; et al. Ankyrin-B dysfunction predisposes to arrhythmogenic cardiomyopathy and is amenable to therapy. J. Clin. Investig. 2019, 129, 3171–3184. [Google Scholar] [CrossRef]
  52. Gan, X.T.; Zhao, G.; Huang, C.X.; Rowe, A.C.; Purdham, D.M.; Karmazyn, M. Identification of fat mass and obesity associated (FTO) protein expression in cardiomyocytes: Regulation by leptin and its contribution to leptin-induced hypertrophy. PLoS ONE 2013, 8, e74235. [Google Scholar] [CrossRef] [Green Version]
  53. Dvornikov, A.V.; Wang, M.; Yang, J.; Zhu, P.; Le, T.; Lin, X.; Cao, H.; Xu, X. Phenotyping an adult zebrafish lamp2 cardiomyopathy model identifies mTOR inhibition as a candidate therapy. J. Mol. Cell Cardiol. 2019, 133, 199–208. [Google Scholar] [CrossRef]
  54. Hanna, A.; Frangogiannis, N.G. The Role of the TGF-beta Superfamily in Myocardial Infarction. Front. Cardiovasc. Med. 2019, 6, 140. [Google Scholar] [CrossRef] [PubMed]
  55. Hu, J.; Gao, C.; Wei, C.; Xue, Y.; Shao, C.; Hao, Y.; Gou, L.T.; Zhou, Y.; Zhang, J.; Ren, S.; et al. RBFox2-miR-34a-Jph2 axis contributes to cardiac decompensation during heart failure. Proc. Natl. Acad. Sci. USA 2019, 116, 6172–6180. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Glenn, D.J.; Cardema, M.C.; Gardner, D.G. Amplification of lipotoxic cardiomyopathy in the VDR gene knockout mouse. J. Steroid. Biochem. Mol. Biol 2016, 164, 292–298. [Google Scholar] [CrossRef] [PubMed]
  57. Liu, L.; Wang, Y.; Cao, Z.Y.; Wang, M.M.; Liu, X.M.; Gao, T.; Hu, Q.K.; Yuan, W.J.; Lin, L. Up-regulated TLR4 in cardiomyocytes exacerbates heart failure after long-term myocardial infarction. J. Cell Mol. Med. 2015, 19, 2728–2740. [Google Scholar] [CrossRef] [PubMed]
  58. Liu, H.; Jing, X.; Dong, A.; Bai, B.; Wang, H. Overexpression of TIMP3 Protects Against Cardiac Ischemia/Reperfusion Injury by Inhibiting Myocardial Apoptosis Through ROS/Mapks Pathway. Cell Physiol. Biochem. 2017, 44, 1011–1023. [Google Scholar] [CrossRef]
  59. Karam, S.; Margaria, J.P.; Bourcier, A.; Mika, D.; Varin, A.; Bedioune, I.; Lindner, M.; Bouadjel, K.; Dessillons, M.; Gaudin, F.; et al. Cardiac Overexpression of PDE4B Blunts beta-Adrenergic Response and Maladaptive Remodeling in Heart Failure. Circulation 2020, 142, 161–174. [Google Scholar] [CrossRef]
  60. Schiattarella, G.G.; Boccella, N.; Paolillo, R.; Cattaneo, F.; Trimarco, V.; Franzone, A.; D’Apice, S.; Giugliano, G.; Rinaldi, L.; Borzacchiello, D.; et al. Loss of Akap1 Exacerbates Pressure Overload-Induced Cardiac Hypertrophy and Heart Failure. Front. Physiol. 2018, 9, 558. [Google Scholar] [CrossRef] [Green Version]
  61. Hullmann, J.E.; Grisanti, L.A.; Makarewich, C.A.; Gao, E.; Gold, J.I.; Chuprun, J.K.; Tilley, D.G.; Houser, S.R.; Koch, W.J. GRK5-mediated exacerbation of pathological cardiac hypertrophy involves facilitation of nuclear NFAT activity. Circ. Res. 2014, 115, 976–985. [Google Scholar] [CrossRef] [Green Version]
  62. Bede-Ojimadu, O.; Amadi, C.N.; Orisakwe, O.E. Blood Lead Levels in Women of Child-Bearing Age in Sub-Saharan Africa: A Systematic Review. Front. Public Health 2018, 6, 367. [Google Scholar] [CrossRef] [Green Version]
  63. Ettinger, A.S.; Egan, K.B.; Homa, D.M.; Brown, M.J. Blood Lead Levels in U.S. Women of Childbearing Age, 1976–2016. Environ. Health Perspect. 2020, 128, 17012. [Google Scholar] [CrossRef] [Green Version]
  64. Bulger, M.; Groudine, M. Enhancers: The abundance and function of regulatory sequences beyond promoters. Dev. Biol. 2010, 339, 250–257. [Google Scholar] [CrossRef] [Green Version]
  65. Clermont, P.L.; Parolia, A.; Liu, H.H.; Helgason, C.D. DNA methylation at enhancer regions: Novel avenues for epigenetic biomarker development. Front. Biosci. 2016, 21, 430–446. [Google Scholar]
  66. Luxan, G.; D’Amato, G.; MacGrogan, D.; de la Pompa, J.L. Endocardial Notch Signaling in Cardiac Development and Disease. Circ. Res. 2016, 118, e1–e18. [Google Scholar] [CrossRef] [PubMed]
  67. Oda, M.; Wakabayashi, S.; Ari Wijetunga, N.; Yuasa, S.; Enomoto, H.; Kaneda, R.; Yoon, S.H.; Mittal, N.; Jing, Q.; Suzuki, M.; et al. Selective modulation of local linkages between active transcription and oxidative demethylation activity shapes cardiomyocyte-specific gene-body epigenetic status in mice. BMC Genom. 2018, 19, 349. [Google Scholar] [CrossRef] [PubMed]
  68. Tsukada, Y. Hydroxylation mediates chromatin demethylation. J. Biochem. 2012, 151, 229–246. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Chen, Z.; Wang, X.; Hou, X.; Ding, F.; Yi, K.; Zhang, P.; You, T. Knockdown of Long Non-Coding RNA AFAP1-AS1 Promoted Viability and Suppressed Death of Cardiomyocytes in Response to I/R In Vitro and In Vivo. J. Cardiovasc. Transl. Res. 2020, 13, 996–1007. [Google Scholar] [CrossRef] [PubMed]
  70. Gurusamy, N.; Das, D.K. Autophagy, redox signaling, and ventricular remodeling. Antioxid. Redox. Signal. 2009, 11, 1975–1988. [Google Scholar] [CrossRef] [Green Version]
  71. Shetty, R.; Joshi, D.; Jain, M.; Vasudevan, M.; Paul, J.C.; Bhat, G.; Banerjee, P.; Abe, T.; Kiyonari, H.; VijayRaghavan, K.; et al. Rudhira/BCAS3 is essential for mouse development and cardiovascular patterning. Sci. Rep. 2018, 8, 5632. [Google Scholar] [CrossRef]
  72. Ye, J.; Beetz, N.; O’Keeffe, S.; Tapia, J.C.; Macpherson, L.; Chen, W.V.; Bassel-Duby, R.; Olson, E.N.; Maniatis, T. hnRNP U protein is required for normal pre-mRNA splicing and postnatal heart development and function. Proc. Natl. Acad. Sci. USA 2015, 112, E3020–E3029. [Google Scholar] [CrossRef] [Green Version]
  73. Silva, M.A.; de Oliveira, T.F.; Almenara, C.C.; Broseghini-Filho, G.B.; Vassallo, D.V.; Padilha, A.S.; Silveira, E.A. Exposure to a Low Lead Concentration Impairs Contractile Machinery in Rat Cardiac Muscle. Biol. Trace Elem. Res. 2015, 167, 280–287. [Google Scholar] [CrossRef]
  74. de Mattos, G.F.; Costa, C.; Savio, F.; Alonso, M.; Nicolson, G.L. Lead poisoning: Acute exposure of the heart to lead ions promotes changes in cardiac function and Cav1.2 ion channels. Biophys. Rev. 2017, 9, 807–825. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Ansari, M.A.; Maayah, Z.H.; Bakheet, S.A.; El-Kadi, A.O.; Korashy, H.M. The role of aryl hydrocarbon receptor signaling pathway in cardiotoxicity of acute lead intoxication in vivo and in vitro rat model. Toxicology 2013, 306, 40–49. [Google Scholar] [CrossRef] [PubMed]
  76. Krishnamoorthy, M.S.; Muthu, P.; Parthiban, N. Influence of pre-, post-, and simultaneous perfusion of elevated calcium on the effect of ascending concentrations of lead on digoxin-induced cardiac arrest in isolated frog heart. Bull. Environ. Contam Toxicol. 1995, 55, 610–617. [Google Scholar] [CrossRef] [PubMed]
  77. Sanchez, O.F.; Lee, J.; Yu King Hing, N.; Kim, S.E.; Freeman, J.L.; Yuan, C. Lead (Pb) exposure reduces global DNA methylation level by non-competitive inhibition and alteration of dnmt expression. Metallomics 2017, 9, 149–160. [Google Scholar] [CrossRef]
  78. Schneider, J.S.; Kidd, S.K.; Anderson, D.W. Influence of developmental lead exposure on expression of DNA methyltransferases and methyl cytosine-binding proteins in hippocampus. Toxicol. Lett. 2013, 217, 75–81. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Sen, A.; Cingolani, P.; Senut, M.C.; Land, S.; Mercado-Garcia, A.; Tellez-Rojo, M.M.; Baccarelli, A.A.; Wright, R.O.; Ruden, D.M. Lead exposure induces changes in 5-hydroxymethylcytosine clusters in CpG islands in human embryonic stem cells and umbilical cord blood. Epigenetics 2015, 10, 607–621. [Google Scholar] [CrossRef]
  80. Paredes, S.R.; Fukuda, H.; Kozicki, P.A.; Rossetti, M.V.; Conti, H.; Batlle, A.M. S-adenosyl-L-methionine and lead intoxication: Its therapeutic effect varying the route of administration. Ecotoxicol. Environ. Saf. 1986, 12, 252–260. [Google Scholar] [CrossRef]
  81. Cao, X.J.; Huang, S.H.; Wang, M.; Chen, J.T.; Ruan, D.Y. S-adenosyl-L-methionine improves impaired hippocampal long-term potentiation and water maze performance induced by developmental lead exposure in rats. Eur. J. Pharm. 2008, 595, 30–34. [Google Scholar] [CrossRef]
  82. Seddik, L.; Bah, T.M.; Aoues, A.; Slimani, M.; Benderdour, M. Elucidation of mechanisms underlying the protective effects of olive leaf extract against lead-induced neurotoxicity in Wistar rats. J. Toxicol. Sci. 2011, 36, 797–809. [Google Scholar] [CrossRef] [Green Version]
  83. Karamian, R.; Komaki, A.; Salehi, I.; Tahmasebi, L.; Komaki, H.; Shahidi, S.; Sarihi, A. Vitamin C reverses lead-induced deficits in hippocampal synaptic plasticity in rats. Brain Res. Bull. 2015, 116, 7–15. [Google Scholar] [CrossRef]
  84. Nam, S.M.; Seo, J.S.; Go, T.H.; Nahm, S.S.; Chang, B.J. Ascorbic Acid Supplementation Prevents the Detrimental Effects of Prenatal and Postnatal Lead Exposure on the Purkinje Cell and Related Proteins in the Cerebellum of Developing Rats. Biol. Trace Elem. Res. 2019, 190, 446–456. [Google Scholar] [CrossRef]
  85. Zeng, Y.; Chen, T. DNA Methylation Reprogramming during Mammalian Development. Genes 2019, 10, 257. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Haddad, R.; Kasneci, A.; Mepham, K.; Sebag, I.A.; Chalifour, L.E. Gestational exposure to diethylstilbestrol alters cardiac structure/function, protein expression and DNA methylation in adult male mice progeny. Toxicol. Appl. Pharm. 2013, 266, 27–37. [Google Scholar] [CrossRef] [PubMed]
  87. Asllanaj, E.; Zhang, X.; Ochoa Rosales, C.; Nano, J.; Bramer, W.M.; Portilla-Fernandez, E.; Braun, K.V.E.; Gonzalez-Jaramillo, V.; Ahrens, W.; Ikram, A.; et al. Sexually dimorphic DNA-methylation in cardiometabolic health: A systematic review. Maturitas 2020, 135, 6–26. [Google Scholar] [CrossRef] [PubMed]
  88. Hartman, R.J.G.; Huisman, S.E.; den Ruijter, H.M. Sex differences in cardiovascular epigenetics-a systematic review. Biol. Sex. Differ. 2018, 9, 19. [Google Scholar] [CrossRef] [Green Version]
  89. Barry, P.S. A comparison of concentrations of lead in human tissues. Br. J. Ind. Med. 1975, 32, 119–139. [Google Scholar] [CrossRef]
  90. Schroeder, H.A. and Tipton, I.H. The human body burden of lead. Arch. Environ. Health 1968, 17, 965–978. [Google Scholar] [CrossRef]
  91. Popovic, M.; McNeill, F.E.; Chettle, D.R.; Webber, C.E.; Lee, C.V.; Kaye, W.E. Impact of occupational exposure on lead levels in women. Environ. Health Perspect. 2005, 113, 478–484. [Google Scholar] [CrossRef] [Green Version]
  92. Berkowitz, G.S.; Wolff, M.S.; Lapinski, R.H.; Todd, A.C. Prospective study of blood and tibia lead in women undergoing surgical menopause. Environ. Health Perspect. 2004, 112, 1673–1678. [Google Scholar] [CrossRef] [Green Version]
  93. Tabish, A.M.; Arif, M.; Song, T.; Elbeck, Z.; Becker, R.C.; Knoll, R.; Sadayappan, S. Association of intronic DNA methylation and hydroxymethylation alterations in the epigenetic etiology of dilated cardiomyopathy. Am. J. Physiol. Heart Circ. Physiol. 2019, 317, H168–H180. [Google Scholar] [CrossRef]
  94. Greco, C.M.; Kunderfranco, P.; Rubino, M.; Larcher, V.; Carullo, P.; Anselmo, A.; Kurz, K.; Carell, T.; Angius, A.; Latronico, M.V.; et al. DNA hydroxymethylation controls cardiomyocyte gene expression in development and hypertrophy. Nat. Commun. 2016, 7, 12418. [Google Scholar] [CrossRef] [PubMed]
Ijerph 18 00577 g001 550
Figure 1. Schematic of experimental design and sample collection. Dams were exposed to Pb 2 weeks prior to mating via drinking water. Maternal (and, in turn, offspring) exposure continued until weaning, when offspring reached 3 weeks of age. 6 males and 6 females per exposure were sacrificed at 5 months of age for tissue collection and ERRBS analysis.
Figure 1. Schematic of experimental design and sample collection. Dams were exposed to Pb 2 weeks prior to mating via drinking water. Maternal (and, in turn, offspring) exposure continued until weaning, when offspring reached 3 weeks of age. 6 males and 6 females per exposure were sacrificed at 5 months of age for tissue collection and ERRBS analysis.
Ijerph 18 00577 g001
Ijerph 18 00577 g002 550
Figure 2. Heart weights, expressed as a percentage of body weight, were assessed for males (A) and females (B) at 5 months of age. Data were analyzed using linear mixed-effects regression with litter-specific random effects to account for within-litter correlation. Black dots represent control animals and red squares depict Pb-treated animals. There were no statistically significant differences in hearts from Pb-exposed animals compared to controls.
Figure 2. Heart weights, expressed as a percentage of body weight, were assessed for males (A) and females (B) at 5 months of age. Data were analyzed using linear mixed-effects regression with litter-specific random effects to account for within-litter correlation. Black dots represent control animals and red squares depict Pb-treated animals. There were no statistically significant differences in hearts from Pb-exposed animals compared to controls.
Ijerph 18 00577 g002
Ijerph 18 00577 g003 550
Figure 3. Volcano plots depicting differentially methylated cytosines (DMCs) for Pb exposed compared to control in males (A) and females (B). DMCs in red did not meet the criteria for significance (at least a 10% change in methylation and FDR < 0.05). Significantly hypomethylated DMCs are shown in blue, and significantly hypermethylated DMCs are shown in green.
Figure 3. Volcano plots depicting differentially methylated cytosines (DMCs) for Pb exposed compared to control in males (A) and females (B). DMCs in red did not meet the criteria for significance (at least a 10% change in methylation and FDR < 0.05). Significantly hypomethylated DMCs are shown in blue, and significantly hypermethylated DMCs are shown in green.
Ijerph 18 00577 g003
Ijerph 18 00577 g004 550
Figure 4. Annotation summary plots depicting the total number of CpGs tested in pink, hypermethylated differentially methylated cytosines (DMCs) in green, and hypomethylated DMCs in blue for each genomic annotation using the R annotatr package for males (A) and females (B). * Denotes statistically significance, p value < 0.05.
Figure 4. Annotation summary plots depicting the total number of CpGs tested in pink, hypermethylated differentially methylated cytosines (DMCs) in green, and hypomethylated DMCs in blue for each genomic annotation using the R annotatr package for males (A) and females (B). * Denotes statistically significance, p value < 0.05.
Ijerph 18 00577 g004
Ijerph 18 00577 g005 550
Figure 5. GREAT pathway analysis of DMRs in males (A) and females (B). Pathways are arranged from most significant (top) to least significant (bottom).
Figure 5. GREAT pathway analysis of DMRs in males (A) and females (B). Pathways are arranged from most significant (top) to least significant (bottom).
Ijerph 18 00577 g005
Ijerph 18 00577 g006 550
Figure 6. Genes associated with differentially methylated regions (DMRs) in males and females were placed into the STRING network interaction tool. DMR-associated genes in females showed no significant network interactions and are not shown in the figure. Panel (A) depicts genes associated with the top 100 most hypermethylated DMRs based on percentage methylation change. Panel (B) depicts all hypomethylated DMRs, since there were fewer than 100. To aid in visualization, only the genes that showed significant network interactivity are shown.
Figure 6. Genes associated with differentially methylated regions (DMRs) in males and females were placed into the STRING network interaction tool. DMR-associated genes in females showed no significant network interactions and are not shown in the figure. Panel (A) depicts genes associated with the top 100 most hypermethylated DMRs based on percentage methylation change. Panel (B) depicts all hypomethylated DMRs, since there were fewer than 100. To aid in visualization, only the genes that showed significant network interactivity are shown.
Ijerph 18 00577 g006
Ijerph 18 00577 g007 550
Figure 7. Venn Diagrams showing overlap in differentially methylated cytosine (DMC) (A) and differentially methylated regions (DMR) (B)-associated genes between males and females. Statistical significance of overlap was determined using a hypergeometric test. In panel B, three DMR-associated genes (Rbfox1, Galnt2, Pi16) that overlap between males and females are highlighted due to their association with cardiovascular disease. Panels (C–E) depict normalized RNA-seq read count data for these three genes in males (top panel) and females (bottom panel) at 5 months of age. Black circles and red squares represent control and Pb-treated animals, respectively. Statistical significance was determined using linear mixed-effects regression, with litter-specific random effects to account for within-litter correlation.
Figure 7. Venn Diagrams showing overlap in differentially methylated cytosine (DMC) (A) and differentially methylated regions (DMR) (B)-associated genes between males and females. Statistical significance of overlap was determined using a hypergeometric test. In panel B, three DMR-associated genes (Rbfox1, Galnt2, Pi16) that overlap between males and females are highlighted due to their association with cardiovascular disease. Panels (C–E) depict normalized RNA-seq read count data for these three genes in males (top panel) and females (bottom panel) at 5 months of age. Black circles and red squares represent control and Pb-treated animals, respectively. Statistical significance was determined using linear mixed-effects regression, with litter-specific random effects to account for within-litter correlation.
Ijerph 18 00577 g007
Table
Table 1. Differentially methylated cytosines (DMCs) in 5 month offspring mouse hearts.
Table 1. Differentially methylated cytosines (DMCs) in 5 month offspring mouse hearts.
ConditionTotal# Hypermethylated (% Total)# Hypomethylated
(% Total)		Total Tested
Female Pb1204850 (71%)354 (29%)1,229,906
Male Pb1441894 (62%)547 (38%)1,242,400
# = number.
Table
Table 2. Differentially methylated regions (DMRs) in 5 month offspring mouse hearts.
Table 2. Differentially methylated regions (DMRs) in 5 month offspring mouse hearts.
ConditionTotal# Hypermethylated (% Total)# Hypomethylated
(% Total)		Total Tested
Female Pb243183 (75%)60 (25%)199,245
Male Pb303200 (66%)103 (34%)200,236
# = number.
Table
Table 3. Top 10 differentially hypo- and hyper-methylated cytosines (DMCs) in each sex, ranked by methylation change.
Table 3. Top 10 differentially hypo- and hyper-methylated cytosines (DMCs) in each sex, ranked by methylation change.
Males
ChromosomeChromosomal CoordinateMethylation
Change		FDRGeneGenomic Annotation
7123,192,72057.240.004Tnrc6aExon, Intron
1527,622,92646.040.009OtulinIntron
1745,811,57445.450.024N/AIntergenic, Intron
5113,144,94945.009.49 × 10−42900026A02RikIntron
Y2,408,70644.973.51 × 10−7N/AIntergenic
575,512,84143.820.009N/AIntergenic
1748,467,93743.360.024Unc5clExon, Intron, 3′UTR
1745,131,08742.981.96 × 10−11N/AIntergenic
1474,943,18642.210.048Lrch1Intron
850,297,70941.540.007N/AIntergenic, Intron
4127,676,072−57.412.6 × 10−4N/AIntergenic
11101,777,136−52.950.026Etv41 to 5 kb, Exon
X151,886,783−51.630.002Huwe1Exon
1858,711,965−42.840.007N/AIntergenic
1946,075,205−41.210.003Nolc1Promoter
515,589,011−41.200.018N/AIntergenic, Intron
560,922,979−41.140.003N/AIntergenic
10126,913,858−40.790.049N/AIntergenic
1482,306,636−40.710.001N/AIntergenic
12112,184,867−40.333.43 × 10−9N/AIntergenic
Females
ChromosomeChromosomal CoordinateMethylation
Change		FDRGeneGenomic Annotation
11100,799,43557.520.016Stat5bIntron
X93,674,70847.422.16 × 10−4Pcyt1bPromoter, Intron
1044,097,41247.170.001Crybg1Intron
11107,243,28245.180.023Pitpnc1Intron
10127,239,42945.160.048Kif5aIntron
2153,872,82243.946.39 × 10−5Bpifb2/Sun51 to 5 kb
5122,761,31342.200.005Camkk2Intron
X138,915,14040.620.001NrkIntron
1049,704,55440.360.008Grik2Intron
1597,844,91640.140.023Hdac7Promoter
1937,685,435−68.117.33 × 10−9Cyp26c1Promoter
5119,687,886−53.300.020Tbx3os2Intron
1573,723,726−45.972.47 × 10−8Ptp4a31 to 5 kb, Intron
1776,107,125−45.740.018N/AIntergenic
175,507,498−45.040.030Zdhhc14Intron
12104,471,048−42.530.047N/A1 to 5 kb, CpG Shore
12104,471,092−42.191.01 × 10−4N/A1 to 5 kb, CpG Shore
170,455,889−41.230.010Spag16Intron
5124,782,320−39.660.021Dnah10Exon
6101,272,100−39.340.050Pdzrn3Intron
FDR = False Discovery Rate.
Table
Table 4. Top 10 differentially hypo- and hyper-methylated regions (DMRs) in each sex, ranked by methylation change.
Table 4. Top 10 differentially hypo- and hyper-methylated regions (DMRs) in each sex, ranked by methylation change.
Males
ChromosomeStartEndMethylation
Change		FDRGeneGenomic Annotation
13119,637,001119,638,00050.970.026Ccl28Intron
1630,745,00130,746,00042.596.93 × 10−4N/AIntergenic
6101,083,001101,084,00041.810.035N/AIntron, Intergenic
536,531,00136,532,00037.560.018Tbc1d14Promoter, 1 to 5 kb, Exon, Intron
834,051,00134,052,00034.430.029N/AIntergenic
575,512,00175,513,00034.130.036N/AIntergenic
814,324,00114,325,00032.020.001Dlgap2Intron
5113,026,001113,027,00031.520.010N/AIntergenic
1786,302,00186,303,00031.150.013PrkceIntron
1286,540,00186,541,00030.638.24 × 10−7N/AIntergenic
1944,248,00144,249,000−71.903.34 × 10−5N/ACpG Island, CpG Shore
142,629,00142,630,000−56.340.010Pantr1Intron
6125,711,001125,712,000−51.240.008Ano2Promoter, Exon, Intron
1185,458,00185,459,000−47.980.022Bcas3Intron
5113,462,001113,463,000−47.456.37 × 10−7N/AIntron
728,925,00128,926,000−43.610.032Actn4Intron
560,922,00160,923,000−41.916.44 × 10−4N/AIntergenic
962,371,00162,372,000−41.690.013Anp32a1 to 5 kb, Exon, Intron
515,589,00115,590,000−41.200.014N/AIntron
1831,921,00131,922,000−40.270.007Lims2Promoter, 1 to 5 kb
Females
ChromosomeStartEndMethylation
Change		FDRGeneGenomic Annotation
1044,097,00144,098,00037.880.044Crybg1Intron
2153,872,001153,873,00037.880.007Sun5Promoter, 1 to 5 kb
2153,872,001153,873,00037.880.007Bpifb21 to 5 kb
1573,544,00173,545,00035.590.027Dennd3Exon, Intron
1254,490,00154,491,00033.840.031N/AIntergenic
1925,422,00125,423,00032.990.041Kank1Promoter, Exon, Intron
1580,000,00180,001,00031.600.024Mir7213Promoter, 1 to 5 kb
1580,000,00180,001,00031.600.024PdgfbExon, Intron
1078,470,00178,471,00031.090.043N/AExon, Intergenic
536,308,00136,309,00031.043.63 × 10−4Sorcs2Intron
659,473,00159,474,000−37.630.002N/AIntergenic
2102,253,001102,254,000−37.550.022N/AIntergenic
734,845,00134,846,000−33.150.037N/AIntergenic
1188,495,00188,496,000−30.700.035Msi2Intron
1240,559,00140,560,000−28.420.022Dock4Intron
745,521,00145,522,000−27.210.044Plekha41 to 5 kb
745,521,00145,522,000−27.210.044Tulp2Exon
1910,904,00110,905,000−26.960.020Prpf19Intron
12102,456,001102,457,000−26.280.044N/AIntergenic, Intron
1873,806,00173,807,000−24.510.036Me2Intron
Table
Table 5. DNA methylation changes at DMRs in cardiovascular disease-associated genes that overlap between males and females.
Table 5. DNA methylation changes at DMRs in cardiovascular disease-associated genes that overlap between males and females.
Males
GeneChrStartEndLead MethCtrl MethChange in MethDirecFDRLocation
Galnt28124,294,001124,295,00030.415.3825.03Hyper9.7 × 10−5Promoter, Intron, 1 to 5 kb
120,995,001120,996,00083.3867.6615.72Hyper0.003Intron
Pi161729,319,00129,320,00042.1315.9226.22Hyper7.8 × 10−9Promoter, 5’UTR, Exon, Intron
Rbfox1167,097,0017,098,00074.7995.40−20.61Hypo0.015Intron
Females
GeneChrStartEndLead MethCtrl MethChange in MethDirecFDRLocation
Galnt28121,760,001121,761,00075.1161.0314.09Hyper0.040Intron
122,085,001122,086,00084.8870.1014.78Hyper0.003Intron
Pi161729,323,00129,324,00064.5647.7116.85Hyper0.024Intron
Rbfox1166,227,0016,228,00049.6039.3610.24Hyper0.014Intron
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Svoboda, L.K.; Wang, K.; Jones, T.R.; Colacino, J.A.; Sartor, M.A.; Dolinoy, D.C. Sex-Specific Alterations in Cardiac DNA Methylation in Adult Mice by Perinatal Lead Exposure. Int. J. Environ. Res. Public Health 2021, 18, 577. https://doi.org/10.3390/ijerph18020577

AMA Style

Svoboda LK, Wang K, Jones TR, Colacino JA, Sartor MA, Dolinoy DC. Sex-Specific Alterations in Cardiac DNA Methylation in Adult Mice by Perinatal Lead Exposure. International Journal of Environmental Research and Public Health. 2021; 18(2):577. https://doi.org/10.3390/ijerph18020577

Chicago/Turabian Style

Svoboda, Laurie K., Kai Wang, Tamara R. Jones, Justin A. Colacino, Maureen A. Sartor, and Dana C. Dolinoy. 2021. "Sex-Specific Alterations in Cardiac DNA Methylation in Adult Mice by Perinatal Lead Exposure" International Journal of Environmental Research and Public Health 18, no. 2: 577. https://doi.org/10.3390/ijerph18020577

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop