- Open Access
In utero exposure to diesel exhaust is associated with alterations in neonatal cardiomyocyte transcription, DNA methylation and metabolic perturbation
Particle and Fibre Toxicologyvolume 16, Article number: 17 (2019)
Developmental exposure to particulate matter air pollution is harmful to cardiovascular health, but the mechanisms by which this exposure mediates susceptibility to heart disease is poorly understood. We have previously shown, in a mouse model, that gestational exposure to diesel exhaust (DE) results in increased cardiac hypertrophy, fibrosis and susceptibility to heart failure in the adult offspring following transverse aortic constriction.
In this study, we have analyzed gene expression in neonatal cardiomyocytes after gestational exposure by RNA-sequencing and have identified 300 genes that are dysregulated, including many involved in cardiac metabolism. We subsequently determined that these cardiomyocytes exhibit reduced metabolic activity as measured by Seahorse extracellular flux analysis. We also surveyed for modifications in DNA methylation at global regulatory regions using reduced representation bisulfite sequencing and found hypomethylation of DNA in neonatal cardiomyocytes isolated from in utero DE exposed neonates.
We have demonstrated that in utero exposure to diesel exhaust alters the neonatal cardiomyocyte transcriptional and epigenetic landscapes, as well as the metabolic capability of these cells. Understanding how exposure alters the developing heart through dysregulation of gene expression, metabolism and DNA methylation is vital for identifying therapeutic interventions for air pollution-related heart failure.
Exposure to particulate matter air pollution has been associated with increased incidence of cardiovascular disease [1,2,3], including arrhythmias , myocardial infarction  and heart failure . Prenatal and early-life exposure to PM2.5 and PM10 is associated with deleterious effects, such as decreased placental mitochondrial DNA , decreased fetal growth , reduced birth weight  and increased newborn systolic blood pressure . Currently, understanding the effects of in utero exposure to DE on newborn and early life health has become a public health concern, as multiple studies link maternal exposure to air pollution with increased risk of congenital heart defects [10, 11], childhood cancers such as retinoblastoma and acute lymphoblastic leukemia , and increased incidence of maternal gestational type 2 diabetes mellitus , which can in turn increase infant risk of metabolic and cardiac disorders . Our lab has shown that in utero exposure to DE in mice predisposes the male offspring to heart failure after transverse aortic constriction (TAC) surgery .
There is growing evidence of alterations in epigenetic control due to PM2.5 exposure. PM2.5 and PM10 exposure have been correlated with significantly decreased global DNA methylation in blood cells as indicated by LINE1 and Alu methylation [16,17,18], as well as changes in targeted loci such as decreased methylation in the iNOS promoter  and increased methylation at the p16 promoter . Diesel exhaust exposure specifically has also been shown to cause dysregulation of epigenetic signatures. Human bronchial epithelial cells exposed to diesel exhaust particulates via culture showed dysregulation of > 60% of microRNAs , and BALB/c mice exposed to diesel exhaust and given allergenic stimuli showed hypermethylation of the interferon-gamma promoter and hypomethylation of interleukin-4 in T-cells . We have also previously shown hypomethylation in exon 1 of GM6307 in neonatal cardiomyocytes (NCMs) resulting from in utero exposure to DE, which correlates with dysregulation of the contiguous miR133a-2 gene .
We have previously identified transcriptional dysregulation in the hearts of male adult mice exposed in utero to DE . To hone in on the early changes that might underlie the previously observed later life susceptibility to heart failure, we have performed RNA-sequencing on NCMs isolated from postnatal day 0 (PND0) offspring exposed to either filtered air (FA) or DE in utero. We observed a large number of transcriptional changes associated with in utero exposure to DE, many of which affect genes involved in normal cardiac metabolism. We have also determined the baseline metabolic phenotype of neonatal cardiomyocytes after gestational exposure and have found that DE exposure is associated with metabolic dysregulation. Additionally, we have performed reduced representation bisulfite sequencing (RRBS) on NCMs, and have found diminished methylation at CpG regions in DE-exposed NCMs. This study focuses on male mice based on our previous findings of in utero DE exposure affecting male adult cardiac outcomes [15, 23], and our observation that in utero DE exposure did not result in any change in female adult cardiac outcomes (C. Weldy, Y. Liu and M.T. Chin, unpublished data). This study is the first of its kind to show widespread cardiomyocyte changes in DNA methylation, gene transcription and metabolism as a result of in utero exposure to DE.
In utero exposure to DE results in altered gene expression
To assess alterations in transcription due to in utero exposure to DE (dams were exposed to a time weighted hourly average of 53 μg/m3 diesel exhaust-created PM2.5, 5 days a week from embryonic day 0.5 (E0.5) to E17.5), RNA-seq was performed on isolated NCMs from FA and DE PND0 offspring. No effects on birth weight, litter size or sex were observed (data not shown). We identified 300 transcripts that were significantly differentially expressed between FA and DE exposure. Table 1 shows the top ten genes with differential expression between exposure groups. Table 2 shows the top ten pathway hits using Qiagen’s Ingenuity Pathway Analysis, with arrows indicating the direction of expression change for genes identified as altered. In both, we can see a clear link between exposure and altered transcription of genes that affect cellular metabolism, showing a consistent decrease in expression in the DE NCMs. Additionally, we observe that many genes involving cell cycle control and cardiomyocyte function were affected. Gene expression differences were validated by qPCR for a select number of top scoring genes, shown in Fig. 1.
In utero exposure to DE-induced changes in gene expression mimic changes observed in hypertrophic NCMs
Given our observations of increased cardiac hypertrophy in adult mice following in utero exposure to DE and TAC surgery [15, 23], we wanted to compare our observed changes in NCM gene expression to other experiments examining the changes in NCM gene expression due to other hypertrophic stressors. As shown in Table 3, using a dataset examining changes in hypertrophic NCM gene expression following transgenic overexpression of CHF1/Hey2 and serum stimulation  (GSE14288), we can observe significant changes in the top ten altered pathways identified in Table 2 in the GSE14288 dataset. When accounting for differential expression following the direction of majority gene change noted in Table 2, we can see that the most of pathways we identified share significant changes, including the mitochondrial L-carnitine shuttle pathway, fatty acid β-oxidation and AMPK signaling. When removing the bias of directional change, we observe that all of our top ten pathways show a significant difference in the GSE14288 dataset. This demonstrates that the gene expression changes we observe in our in utero DE NCMs are also present in NCMs that undergo hypertrophic stimuli.
In utero exposure to DE alters neonatal cardiomyocyte oxygen consumption
To determine whether the transcriptional changes identified in Tables 1 and 2 correlated to a functional change in the metabolic capacity of NCMs, we tested the rate of oxygen consumption of DE and FA NCMs using the Seahorse XFe24 Analyzer. However, since the switch from glucose oxidation to fatty acid oxidation in the neonatal heart doesn’t occur until sometime between postnatal day 1–7 [25, 26], we tested the response of NCMs using pyruvate/glutamate to fuel respiration. Figure 2a shows the average oxygen consumption of each group at baseline and after the injections of oligomycin, carbonyl cyanide-p-trifluoromethoxyphenylhydrazone (FCCP), and Antimycin A/Rotenone. After subtracting the non-mitochondrial respiration, we found that both the baseline and maximal respiration are significantly lower in the DE NCMs (Fig. 2b and c). The ATP production and coupling efficiency of the DE NCMs was also found to be significantly lower than the FA NCMs (Fig. 2d and f). DE NCMs also demonstrated a lower spare respiratory capacity, though it was not statistically significant (Fig. 2e). There was no significant difference in extracellular acidification rate between FA and DE NCMs at any of these stages (data not shown). Mitochondrial DNA content was also not significantly different between the two exposure groups (data not shown).
In utero exposure to DE results in altered DNA methylation patterning
The temporal delay between in utero exposure, changes in adult gene expression  and adult onset of disease  suggests that an epigenetic mechanism may mediate the observed effects on adult increased susceptibility to heart failure. We performed reduced representation bisulfite sequencing (RRBS) on isolated NCMs from PND0 mice exposed to in utero FA or DE to uncover alterations in DNA methylation at CpG-rich regions, which are potential regions of genetic control. We identified 63 differentially methylated regions (DMRs), of which only 1 was hypermethylated in the DE samples. The top ten DMRs and their gene calls are represented in Table 4. To determine whether any particular pathways were being selectively altered in their DNA methylation as a result of exposure, we performed pathway analysis using Qiagen Ingenuity Pathway Analysis (Qiagen). Genes located within DMRs were given a hypothesized significant change of expression, with decreased methylation given a call of increased expression and vice versa. The top ten pathways showing altered DNA methylation are presented in Table 5. The finding of significantly reduced methylation in genes involved in pathways such as β-adrenergic, GPCR, and hypertrophic signaling suggests a mechanism by which gestational DE exposure predisposes to adult heart failure.
In utero exposure to DE significantly decreases DNA methylation in regulatory regions
The discovery of 62 DMRs showing decreased methylation and only 1 showing an increase in the DE NCMs suggests a significant loss of methylation in the DE NCMs. This general reduction in overall DNA methylation is further demonstrated in Fig. 3a, where we examine the frequency of each CpG showing a specific percent methylation. We observe a clear decrease in the amount of methylation present in the DE samples, both through a higher proportion of 0–5% methylated CpGs and the lower amount of CpGs that were found to be methylated in > 25% of the sequences. Additionally, the genomic location of the identified DMRs is graphed in Fig. 3b, with the majority of DMRs located in intergenic and intronic regions. The genes Gnas, Gng12 and Pde6H were heavily implicated in many of the altered pathways, and the DMRs associated with these genes were found in the intronic (Gnas) and intergenic (Gng12, Pde6H) regions.
When evaluating the individual CpGs that were identified to have significantly different methylation status between the FA and DE NCMs given a false discovery rate (FDR) of 0.1, including those that were incorporated into DMRs, we identified 12,629 differentially methylated CpGs (564 of which were included in DMRs). Of those, only 293 showed increased methylation in the DE NCMs. Assuming a null hypothesis of differential methylation affecting the samples equally, under the binomial distribution we would reject the null hypothesis with a p-value of 2.2e-16. Taken together, these indicate that in utero exposure to DE significantly decreases DNA methylation in CpG-rich regions.
Comparison between transcriptional and DNA methylation alterations
In order to determine how the observed change in DNA methylation may link to the transcriptional alterations, the datasets were compared by identifying any genes within 1 Mb of a DMR that showed significant changes in expression. Surprisingly, only 14 genes were found to have significantly altered transcription near identified DMRs: Acacb, Cdc45, Chst12, Cox4i2, Cys1, Etv4, Fkbp5, Glipr2, Gng11, Melk, Pibd1, Sh3gl3, Tubb3 and Vwa8. However, only Etv4 and Glipr2 were identified to have associated DMRs in the more stringent assessment of altered DNA methylation discussed above. Of these genes, Acacb and Cox4i2 are involved in cellular metabolism, and Gng11 is involved in GPCR signaling. Acetyl-CoA Carboxylase Beta (Acacb) is interesting due to its role in catalyzing the conversion of acetyl-CoA to malonyl-CoA in the fatty acid cycle, as well as being a key regulator of fatty acid oxidation through the ability of malonyl-CoA to block carnitine-palmitoyl-CoA transferase . Cytochrome c oxidase subunit 4 isoform 2 (Cox4i2) is a subunit of the cytochrome c oxidase complex within the electron transport chain . G protein subunit gamma 11 (Gng11) is a subunit of the GPCR signaling complex, which has also been described as a controller of cellular senescence . While these genes are biologically interesting, we would expect to see a stronger link between changes in DNA methylation and changes in gene expression. It is possible that the altered DNA methylation would not manifest in transcriptional changes until the cells are stressed, or at a later age.
We have shown that the transcriptional profile of NCMs is significantly altered as a result of in utero exposure to DE, in large part affecting genes associated with metabolic pathways (Tables 1, 2). Metabolic dysregulation is a hallmark feature of heart failure , with decreased fatty acid oxidation and increased glycolysis playing an important role in late-stage heart failure [31, 32]. While it is unclear how these findings might affect adult susceptibility to heart failure, multiple studies have discovered a link between exposure to air pollution and metabolic dysregulation [33,34,35]. Combined with our prior findings of increased expression of Ptprf , a transcriptional change associated with decreased sensitivity to insulin signaling , this suggests a potential role for metabolic dysregulation as a causal link between in utero DE exposure and adult susceptibility to heart failure. Additionally, we have shown that respiratory capabilities of NCMs from in utero DE exposed neonates are hindered at both baseline and maximal respiration (Fig. 2b and c). NCMs from in utero DE exposed neonates also demonstrate decreased ATP production and coupling efficiency (Fig. 2d and f), indicative of higher energy loss and ROS production . Given the known role of impaired metabolic activity on heart failure , these findings suggest that metabolic dysregulation during development may increase the susceptibility to heart failure at adulthood in our model. However, the persistence of the observed transcriptional alterations and the observed metabolic defects through adulthood and the causal relationships with the development of heart failure remain to be elucidated. Overall, our findings of transcriptional dysregulation in NCMs as a result of in utero DE exposure will aid in guiding the experimental focus for determining modifications in cardiac function caused by exposure.
While the observed changes in metabolic transcription and function are potentially important mediators of in utero exposure-induced susceptibility to heart failure, other mechanisms involving cardiac remodeling and fetal gene expression may be involved. One gene that showed differential expression was the alpha myosin heavy gene, Myh6. The RNA-seq profile showed a significant reduction in expression of Myh6 in DE NCMs compared to FA. This gene is essential in cardiac development, with ablation experiments showing that Myh6−/− is embryonic lethal, and Myh6+/− mice showing impaired contractility and relaxation, as well as irregular sarcomere structure . MYH6 mutations in humans have been associated with congenital heart defects such as atrial septal defect , as well as in hypertrophic and dilated cardiomyopathy [41, 42]. Further exploration into the dysregulation of Myh6 in this model would aid in understanding the potential developmental defects associated with in utero DE exposure.
Our prior findings of adult cardiovascular disease following in utero exposure suggest a temporally durable mechanism of action. Alterations of epigenetics through DNA hypomethylation have been implicated in the development and progression of cardiovascular disease [43, 44], and due to the double-stranded nature of CpG methylation these changes are incredibly stable . Changes in DNA methylation both at the targeted level and genome-wide have been associated with exposure to air pollution [22, 46]. We are the first to demonstrate a change in DNA methylation at CpG-rich regions of cardiomyocytes as a result of in utero exposure to DE. The discovery of decreased methylation in DMRs associated with genes involved in GPCR and cardiac hypertrophic signaling suggests a predisposition to dysregulation in these pathways. However, there was no observed increase in transcription of Gnas, Gng12, Adra2c or Pde6h in the transcriptional profile of DE-exposed NCMs, though Gng11 was identified to have increased expression. It is possible that the changes in transcription affected by DNA methylation would be stimulated through stress, or will manifest themselves sometime beyond the neonatal period. GPCR signaling also plays a role in the autonomic nervous system control of cardiac regulation. Considering our prior findings of decreased blood pressure at baseline in in utero DE exposed 10 week old mice , it is important to further understand how these molecular changes are affecting physiological outcomes in both the neonatal and adult hearts. Future studies targeting autonomic features such as heart rate variability would likely be informative.
It is unclear whether this change in DNA methylation is occurring through a directed or undirected pathway. Exposure to PM2.5 has been shown to have an associated increase in oxidative stress [47,48,49]. Passive removal of DNA methylation can occur due to oxidative damage from ROS by damaging the DNA to the extent that DNA methlyltransferases (DNMTs) are unable to bind to the lesion area, thereby averting the maintenance of methylation in daughter cells [50, 51]. If this process is occurring early in development during the repatterning of methylation, these changes would be extremely heritable, and likely occurring in multiple cell and organ types . It is also possible that cellular damage due to oxidative stress prevents DNA methylation through sequestration of glutathione . Targeted changes in DNA methylation in response to air pollution have been described previously. Hypermethylation of the p16 promoter in alveolar epithelial cells of mice exposed to PM2.5 in conjunction with increased DNMT1 expression has been reported  and decreased methylation on the first exon of Gm6307, which drives miR133a-2 expression, has also been reported in association with in utero DE exposure . Our studies of transcriptional dysregulation in NCMs did not identify any DNA methyltransferase or TET family enzyme genes as being differentially expressed, though this does not preclude the possibility of altered expression earlier in development leading to decreased methylation that was maintained until birth.
For future studies, refining the contributions of in utero DE exposure and later life cardiac outcomes will be important for focusing future experiments. Meta-analysis of epidemiological studies have demonstrated the increased risk of congenital heart defects associated with maternal exposure to nitrogen dioxide, sulfur dioxide, and particulate matter [53, 54]. Although we favor particulates as being likely responsible for the observed effects, contributions from other components cannot be ruled out. Additionally, it is not clear whether maternal effects, such as placental insufficiency brought on from inflammatory response to exposure, are the main influencers on the offspring cardiac outcomes, or if particulate and/or gaseous matter is able to cross the placental barrier and directly interact with the cardiovascular system of the developing embryos.
Our study is the first to systematically analyze and demonstrate the profound effects of in utero DE exposure on cardiomyocyte gene transcription, DNA methylation and metabolism. It is unclear whether these effects are due to a direct effect of DE exposure or an indirect mechanism via the maternal response to exposure, although our previous observation that maternal exposure leads to placental inflammation  is consistent with an indirect mechanism. A limitation of our study is that the correlation between mRNA expression, DNA methylation and metabolic activity is imperfect. Nevertheless, these findings provide important datasets for mining the potential pathways that mediate the deleterious effects of diesel exhaust on the developing heart that predispose to heart disease. Identifying pathological pathways that promote the development of air pollution-related heart disease is a first step in developing potential therapies.
Diesel exhaust exposure and mice
Housing and exposure facility conditions were carried out as previously described [15, 22, 23]. C57Bl/6 J mice were used in these experiments, housed in specific pathogen free (SPF) conditions, on a 12/12-light/dark cycle. Female and male mice between the ages of 12 to 16 weeks were transferred to our Northlake Diesel Exposure Facility located near the University of Washington (UW) and housed under SPF conditions in Allentown caging systems. Female mice were paired with male mice for timed mating in FA. After observation of a vaginal plug, pregnant mice were put into FA or DE with exposures beginning at E0.5 and lasting until E17.5, at which point pregnant mice were transferred from the Northlake facility to the UW Medicine South Lake Union SPF vivarium for neonatal sample collection.
Diesel exhaust was generated from a single cylinder Yanmar diesel engine (Model YDG5500EV-6EI) operating on 82% load. A detailed analysis of DE particulate components in this system has been previously reported . DE exposures were conducted for 6 h per day (9 am − 3 pm) 5 days a week (Monday – Friday) and DE concentrations were regulated to 300 μg/m3 of PM2.5. A 300 μg/m3 concentration of PM2.5 6 hours a day, 5 days a week equates to a time weighted hourly average of 53 μg/m3. The exposure characteristics detailing gas, particle-bound polycyclic aromatic hydrocarbons (PAH), and particle diameter of this system have been previously reported [55, 56]. In previous measurements, oxides of nitrogen concentrations were 1800 ppb NOx and 60 ppb NO2, carbon monoxide was 2 ppm, and carbon dioxide was 1000 ppm. The mass fraction of particle-bound PAH was 20 ng/mg PM2.5 and the ratio of the organic carbon to elemental carbon mass concentration was 0.10. The mass median aerodynamic diameter of particles was 85 nm (GSD 1.2) and the count median thermodynamic equivalent diameter was 87 nm (GSD 3.0). The exposure approximates real world ambient exposure in an urban setting.
Isolation of NCMs
Upon birth, PND0 hearts were harvested, trimmed of surrounding vascular and atrial tissue, and dissociated as previously described [57, 58]. Males were detected through post-mortem dissection and identification of testes. Dissociation was performed using 1 mg/mL Liberase TH (Roche; Pleasanton, CA, USA) in 1X HBSS by incubating hearts at 37 degrees for 5 min in solution, with pipetting to release cells after incubation. Media containing released NCMs was adjusted to 20% FBS-DMEM, and cellular dissociation was continued until the majority of cells were released. Cells were then filtered using a 70 μm sieve, re-eluted in 20% FBS-DMEM with 20 μm Ara-C and incubated at 37 degrees for 1 h to allow fibroblasts to attach onto the plate. After incubation, media with suspended cardiomyocytes was carefully removed, spun and purified cardiomyocyte pellets were collected and frozen at − 80 °C.
RNA was purified from frozen PND0 isolated NCMs using Trizol (Thermo Fisher Scientific, Waltham, MA, USA), following the manufacturer’s protocol. For each sequencing sample, NCMs from one male neonatal heart were isolated and used for each n (n = 3 for FA, DE). Ribosomal RNA was depleted by poly-A enrichment, and sample libraries were created using the TruSeq Stranded mRNA kit (Illumina, San Diego, CA). Each library was barcoded using the Illumina adapters, and amplified with 13 cycles of PCR. Library concentrations were quantified using the Quant-it dsDNA Assay (Life Technologies, Carlsbad, CA). Libraries were normalized and pooled based on Agilent 2100 Bioanalyzer results (Agilent Technologies, Santa Clara, CA), and the pools were sequenced on an Illumina HiSeq 2500.
To process sequences, Illumina RTA-generated BCL files were converted to FASTQ files and sequence read and base quality were checked using the FASTX-toolkit (http://hannonlab.cshl.edu/fastx_toolkit/) and FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Sequences were aligned to mm10 with reference transcriptome Ensembl v67 using Tophat . Lane level bam data files were merged using the Picard MergeSamFiles tool and suspected PCR duplicates were marked, not removed, in the alignment files using the Picard MarkDuplicates tool (http://broadinstitute.github.io/picard/). Local realignment was performed around indels, and base quality score recalibration was run using GATK tools . Variant detection was performed with the GATK Unified Genotyper version 2.6.5 , and aligned data were used for isoform assembly and quantitation with Cufflinks [59, 62]. Fragment counts were generated using the featureCounts function in the Bioconductor Rsubread package . Comparisons were performed using the Bioconductor edgeR package. Genes were identified as significantly different between groups if they had a 25% or greater difference in expression between the two groups at an FDR adjusted p-value of 0.05 or less. To assess consistency between our data and previous results (GSE14288) we performed self-contained gene set tests using the eleven pathways we identified using IPA (Table 2), testing for evidence that those gene sets were differentially expressed in the previous results as well. We used the limma fry function to perform the gene set tests .
Select genes were validated by qPCR using iTaq Universal Sybr Green Supermix (BioRad). Primers used for validation were as follows:
Acot1-For: 5′-TACGATGACCTCCCCAAGAA-3′, Acot1-Rev: 5′-AGCCCAATTCCAGGTCCTTT-3′, Acot2-For: 5′-ACAACTACGACGACCTCCCC-3′, Acot2-Rev: 5′-AGCCCAATTCCAGGTCCTT-3′, Cpt1a-For: 5′-TCCATGCATACCAAAGTGGAC-3′, Cpt1a-Rev: 5′-CGATGTTCTTCGTCTGGCTT-3′, Cpt1b-For: 5′-AGCCCCCTCATGGTGAAC-3′, Cpt1b-Rev: 5′-AGTTTGCGGCGATACATGA-3′, Scd4-For: 5′-TGGAGCCACCGAACCTATAA-3′, Scd4-Red: 5′-GGCCCATTCATACACGTCAT-3′, 18S-For: 5′-GGACAGGATTGACAGATTGATAG-3′, 18S-Rev: 5′-ATCGCTCCACCAACTAAGAA-3′. Expression was normalized to ribosomal RNA 18S. Samples were run on the 7900HT Fast Real-Time PCR System (ABI). Statistical analysis was carried out using a t-test for comparison of ddCT values.
Metabolic capacity analysis
Cellular respiration was performed on NCMs using the mitochondrial stress test (Agilent; Santa Clara, CA, USA), measured using a Seahorse XFe24 extracellular flux analyzer (Agilent). NCMs were plated on provided 24 well plate coated in fibronectin, plating 3–5 replicates per sample, and cultured for 3 days in 20% FBS-DMEM with 20 μM Ara-C. This time period allowed for the non cardiomyocyte cells to become metabolically inactive after the isolation procedure. On the day of measurement, cells were washed with XF assay medium (Agilent) supplemented with 25 mM glucose, 1 mM sodium pyruvate and 2 mM L-glutamine, and replenished with the same media and incubated at 37 °C without CO2 an hour prior to measurement. Injection scheme is represented in Fig. 1a, with concentrations at 1.5 μM oligomycin, 1 μM FCCP, 1 μM Rotenone and 1 μM antimycin A. Normalization to μg DNA/well was performed by phenol:chloroform extracting DNA from each well, and measuring DNA concentration using the Quant-iT PicoGreen dsDNA assay kit (Thermo Fisher Scientific), following the manufacturers protocol.
Rates were averaged between the 3–5 replicate plated wells per sample, giving one averaged read per sample n. Calculations for basal and maximal respiration, ATP production, spare respiratory capacity, and coupling efficiency were performed as recommended by Agilent:
Non-mitochondrial oxygen consumption (NMOC) = [min. Rate measurement after rotenone/antimycin A injection].
Basal respiration = [last rate measured before first (oligomycin) injection] - NMOC.
Maximal respiration = [max. Rate measured after FCCP injection] – NMOC.
ATP production = [last rate measured before oligomycin] – [min. Rate measurement after oligomycin injection].
Spare respiratory capacity = [maximal respiration] – [basal respiration].
Coupling efficiency = [ATP production rate] / [basal respiration] × 100.
Statistical comparisons were performed using a t-test between the measurements made in male DE vs. FA NCMs (n = 7, n = 8 respectively).
Reduced representation bisulfite sequencing
DNA was isolated from frozen PND0 isolated cardiomyocytes using Qiagen (Germantown; MD, USA) DNeasy Blood and Tissue kit, according to the manufacturer’s protocol. For each sequencing sample, NCMs from one male neonatal heart were isolated and used for each n (n = 4 for FA, DE). 500 ng of DNA was digested overnight with MspI (New England BioLabs Inc., Ipswich, MA, USA) and then directly prepared with the KAPA Hyper Library prep protocol (KAPA Biosystems, Wilmington, MA, USA). Samples were then adapter-ligated using SeqCap Adapters (Roche-Nimblegen, Pleasanton, CA, USA), and cleaned using a 2.5X post-ligation Agencourt AMPure XP bead cleanup (Beckman Coulter, Indianapolis, IN). Size selection was performed to a size range of 160 bp–340 bp on a 2% Pippin Prep gel (Sage Science, Beverly, MA, USA), followed by bisulfite conversion using the Zymo EZ DNA Methylation Lightning kit (Zymo Research, Irvine, CA, USA). The converted DNA was then put back into the KAPA Hyper Library prep protocol at the library amplification step, and amplification was performed with 19 PCR cycles using 2X KAPA Hifi HotStart Uracil+ ReadyMix (KAPA Biosystems, Wilmington, MA, USA). Post-amplification cleanup was performed with 0.8X Agencourt AMPure XP beads (Beckman Coulter). Library size distributions were validated using the Agilent High Sensitivity D1000 ScreenTape run on an Agilent 2200 TapeStation (Agilent Technologies, Inc., Santa Clara, CA, USA). Additional library QC, blending of pooled indexed libraries, and cluster optimization was performed using Life Technologies- Invitrogen Qubit® 2.0 Fluorometer (Life Technologies-Invitrogen, Carlsbad, CA, USA). Sample libraries were sequenced on the Illumina HiSeq 2500 (Illumina; San Diego, CA, USA).
Sequencing reads were processed by discarding the low quality reads using Illumina’s base call quality filter, and raw reads were trimmed of adapter sequences and flanking low quality base reads using the trimmomatic read trimmer , discarding any reads less than 25 nt, after which there were between 13.6–31.4 million reads per sample. Methylation status for each CpG was estimated by aligning reads to the mm10 genome using the Bismark aligner , read into R using the Bioconductor bsseq package , and differential methylation was called using the Bioconductor DSS package . There were between 2.85 and 4.47 million CpGs per sample. Low read coverage CpGs were filtered, based on the criteria that a CpG had to have read coverage > 5 in at least half of the samples, reducing the total number of CpGs to ~ 646 thousand. After this filtering, median read depth per CpG was 8–25. We then selected a set of differentially methylated regions, based on a set of criteria; to be considered a differentially methylated region (DMR), a) the region had to be at least 50 bp in length with at least 3 CpGs present, and b) at least 50% of the CpGs in the region had to be significant, based on a p < 0.001 and show a minimum 10% difference in methylation. Any significant DMRs within 100 bp of each other were merged to a single DMR. Methylation and gene expression data were from different mice, so it wasn’t possible to do any direct analysis between the two. We instead simply looked for any genes within 1 Mb of each DMR that were differentially expressed with an inverse sign (e.g., genes that increase expression as methylation decreases and vice versa).
- Acacb :
acetyl-CoA Carboxylase Beta
- Cox4i2 :
cytochrome c oxidase subunit 4 isoform 2
differentially methylated region
- DNMT :
- Gng11 :
G protein subunit gamma 11
polycyclic aromatic hydrocarbon
- PM10 :
particulate matter with a diameter of 10 μm or less
- PM2.5 :
particulate matter with a diameter of 2.5 μm or less
reduced representation bisulfite sequencing
transverse aortic constriction
Brook R, Rajagopalan S, Pope CA 3rd, Brook JR, Bhatnagar A, Diez-Roux AV, et al. Particulate matter air pollution and cardiovascular disease: an update to the scientific statement from the American Heart Association. Circulation. 2010;121 21:2331–2378.
Krishnan T, Adar SD, Szpiro AA, Jorgensen NW, van Hee VC, Barr RG, et al. Vascular responses to long- and short-term exposure to fine particulate matter: MESA air (multi-ethnic study of atherosclerosis and air pollution). J Am Coll Cardiol. 2012;60(21):2158–66.
Cr P, Muhlestein JB, May HT, Renlund DG, Anderson JL, Horne BD. Ischemic heart disease events triggered by short-term exposure to fine particulate air pollution. Circulation. 2006;114(23):2443–8.
Milojevic A, Wilkinson P, Armstrong B, Bhaskaran K, Smeeth L, Hajat S. Short-term effects of air pollution on a range of cardiovascular events in England and Wales: case-crossover analysis of the MINAP database, hospital admissions and mortality. Heart. 2014;100:1093–8.
To T, Zhu J, Villeneuve PJ, Simatovic J, Feldman L, Gao C, Williams D, Chen J, Weichenthal S, Wall C, Miller AB. Chronic disease prevalence in women and air pollution--a 30-year longitudinal cohort study. Environ Int. 2015;80:26–32.
Janssen B, Munters E, Pieters N, Smeets K, Cox B, Cuypers A, Fierens F, Penders J, Vangronsveld J, Gyselaers W, Nawrot TS. Placental mitochondrial DNA contect and particulate air pollution during in utero life. Environ Health Perspect. 2012;120(9):1346–52.
van den Hooven E, Pierik FH, de Kluizenaar Y, Willemsen SP, Hofman A, van Ratingen SW, et al. Air pollution exposure during pregnancy, ultrasound measures of fetal growth, and adverse birth outcomes: a prospective cohort study. Environ Health Perspect. 2012;120(1):150–6.
Dadvand P, Parker J, Bell ML, Bonzini M, Brauer M, Darrow LA, et al. Maternal exposure to particulate air pollution and term birth weight: a mutli-country evaluation of effect and heterogeneity. Environ Health Perspect. 2013;121(3):367–73.
van Rossem L, Rifas-Shiman SL, Melly SJ, Kloog I, Luttmann-Gibson H, Zanobetti A, et al. Prenatal air pollution exposure and newborn blood pressure. Environ Health Perspect. 2015;123(4):353–9.
Dolk H, Armstrong B, Lachowycz K, Vrijheid M, Rankin J, et al. Ambient air pollution and risk of congenital anomalies in England, 1991-1999. Occup Environ Med. 2010;67(4):223–7.
Gilboa S, Mendola P, Olshan AF, Langlois PH, Savitz DA, et al. Relation between ambient air quality and selected birth defects, seven county study, Texas, 1997-2000. Am J Epidemiol. 2005;162(3):238–52.
Heck J, Wu J, Lombardi C, Qiu J, Meyers TJ, et al. Childhood cancer and traffic-related air pollution exposure in pregnancy and early life. Environ Health Perspect. 2013;121:1385–91.
Eze I, Hemkens LG, Bucher HC, Hoffmann B, Schindler C, et al. Association between ambient air pollution and diabetes mellitus in Europe and North America: systematic review and meta-analysis. Environ Health Perspect. 2015;123(5):381–9.
Mitanchez D, Yzydorczyk C, Simeoni U. What neonatal complications should the pediatrician be aware of in case of maternal gestational diabetes? World J Diabetes. 2015;6(5):734–43.
Weldy C, Liu Y, Liggitt HD, Chin MT. In utero exposure to diesel exhaust air pollution promotes adverse intrauterine conditions, resulting in weight gain, altered blood pressure, and increased susceptibility to heart failure in adult mice. PLoS One. 2014;9(2):e88582.
Baccarelli A, Wright RO, Bollati V, Tarantini L, Litonjua AA, Suh HH, et al. Rapid DNA methylation changes after exposure to traffic particles. Am J Respir Crit Care Med. 2009;179(7):572–8.
Tarantini L, Bonzini M, Apostoli P, Pegoraro V, Bollati V, Marinelli B, et al. Effects of particulate matter on genomic DNA methylation content and iNOS promoter methylation. Environ Health Perspect. 2008;117(2):217–22.
Bellavia A, Urch B, Speck M, Brook RD, Scott JA, et al. DNA hypomethylation, ambient particulate matter, and increased blood pressure: findings from controlled human exposure experiments. J Am Heart Assoc. 2013;2(3):e000212.
Soberanes S, Gonzalez A, Urich D, Chiarella SE, Radigan KA, et al. Particulate matter air pollution induces hypermethylation of the p16 promoter via a mitochondrial ROS-JNK-DNMT1 pathway. Sci Rep. 2012;2:275.
Jardim M, Fry RC, Jaspers I, Dailey L, Diaz-Sanchez D. Disruption of microRNA expression in human airway cells by diesel exhaust particles is linked to tumorigenesis-associated pathways. Environ Health Perspect. 2009;117(11):1745–51.
Liu J, Ballaney M, Al-alem U, Quan C, Jin X, Perera F, Chen L, Miller RL. Combined inhaled diesel exhaust particles and allergen exposure alter methylation of T helper genes and IgE production in vivo. Toxicol Sci. 2008;102(1):76–81.
Goodson J, Weldy CS, MacDonald JW, Liu Y, Bammler TK, et al. 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.
Weldy C, Liu Y, Chang YC, Medvedev IO, Fox JR, Larson TV, et al. In utero and early life exposure to diesel exhaust air pollution increases adult susceptibility to heart failure in mice. Part Fibre Toxicol. 2013;10(1):59.
Yu M, Xiang F, Beyer RP, Farin FM, Bammler TK, Chin MT. Transcription factor CHF1/Hey2 regulates specific pathways in serum stimulated primary cardiac myocytes: implications for cardiac hypertrophy. Curr Genomics. 2010;11(4):287–96.
Taegtmeyer H, Sen S, Vela D. Return to the fetal gene program: a suggested metabolic link to gene expression in the heart. Ann N Y Acad Sci. 2010;1188:191–8.
Onay-Besikci A, Campbell FM, Hopkins TA, Dyck JR, Lopaschuk GD. Relative importance of malonyl CoA and carnitine in maturation of fatty acid oxidation in newborn rabbit heart. Am J Physiol Heart Circ Physiol. 2003;284(1):H283–9.
Abu-Elheiga L, Matzuk MM, Abo-Hashema KA, Wakil SJ. Continuous fatty acid oxidation and reduced fat storage in mice lacking acetyl-CoA carboxylase 2. Science. 2001;291(5513):2613–6.
Hüttemann M, Kadenbach B, Grossman LI. Mammalian subunit IV isoforms of cytochrome c oxidase. Gene. 2001;267(1):111–23.
Hossain MN, Sakemura R, Fujii M, Ayusawa D. G-protein gamma subunit GNG11 strongly regulates cellular senescence. Biochem Biophys Res Commun. 2006;351(3):645–50.
Ingwall J. Energy metabolism in heart failure and remodelling. Cardiovasc Res. 2008;81(3):412–9.
Fürstenwerth H. Rethinking heart failure. Cardiol Res. 2012;3(6):243–57.
Stanley W, Recchia FA, Lopaschuk GD. Myocardial substrate metabolism in the normal and failing heart. Physiol Rev. 2005;85(3):1093–129.
Wei Y, Zhang JJ, Li Z, Gow A, Chung KF, et al. Chronic exposure to air pollution particles increases the risk of obesity and metabolic syndrome: findings from a natural experiment in Beijing. FASEB J. 2016;30(6):2115–22.
Rajagopalan S, Brook RD. Air pollution and type 2 diabetes: mechanistic insights. Diabetes. 2012;61(12):3037–45.
Eze I, Schaffner E, Foraster M, Imboden M, von Eckardstein A, et al. Long-term exposure to ambient air pollution and metabolic syndrome in adults. PLoS One. 2015;10(6):e0130337.
Zabolotny J, Kim YB, Peroni OD, Kim JK, Pani MA, Boss O, et al. Overexpression of the LAR (leukocyte antigen-related) protein-tyrosin phosphatase in muscle causes insulin resistance. PNAS. 2001;98(9):5187–92.
Divakaruni A, Brand MD. The regulation and physiology of mitochondrial proton leak. Physiology. 2011;26(3):192–205.
Ventura-Clapier R, Garnier A, Veksler V. Energy metabolism in heart failure. J Physiol. 2004;555(Pt 1):1–13.
Jones WK, Grupp IL, Doetschman T, Grupp G, Osinska H, Hewett TE, et al. Ablation of the murine alpha myosin heavy chain gene leads to dosage effects and functional deficits in the heart. J Clin Investig. 1996;98(8):1906–17.
Ching YH, Ghosh TK, Cross SJ, Packham EA, Honeyman L, Loughna S, et al. Mutation in myosin heavy chain 6 causes atrial septal defect. Nat Genet. 2005;37(4):423–8.
Carniel E, Taylor MRG, Sinagra G, Di Lenarda A, Ku L, Fain PR, et al. α-Myosin heavy chain: a Sarcomeric gene associated with dilated and hypertrophic phenotypes of cardiomyopathy. Circulation. 2005;112(1):54–9.
Hershberger RE, Norton N, Morales A, Li D, Siegfried JD, Gonzalez-Quintana J. Coding sequence rare variants identified in MYBPC3, MYH6, TPM1, TNNC1, and TNNI3 from 312 patients with familial or idiopathic dilated cardiomyopathy. Circ Cardiovasc Genet. 2010;3(2):155–61.
Khalil C. The emerging role of epigenetics in cardiovascular disease. Therapeutic Advances in Chronic Disease. 2014;5(4):178–87.
Zhong J, Agha G, Baccarelli AA. The role of DNA methylation in cardiovascular risk and disease. Circ Res. 2016;118:119–31.
Bird A. DNA methylation patterns and epigenetic memory. Genes Dev. 2002;16(1):6–21.
Plusquin M, Guida F, Polidoro S, Vermeulen R, Raaschou-Nielsen O, et al. DNA methylation and exposure to ambient air pollution in two prospective cohorts. Environ Int. 2017;108:127–36.
Sorensen M, Daneshvar B, Hansen M, Dragsted LO, Hertel O, Knudsen L, et al. Personal PM2.5 exposure and markers of oxidative stress in blood. Environ Health Perspect. 2003;111(2):161–6.
Riva DR, Magalhaes CB, Lopes AA, Lancas T, Mauad T, Malm O, et al. Low dose of fine particulate matter (PM2.5) can induce acute oxidative stress, inflammation and pulmonary impairment in healthy mice. Inhal Toxicol. 2011;23(5):257–67.
Lodovici M, Bigagli E. Oxidative stress and air pollution exposure. Journal of Toxicology. 2011;2011:1–9.
Turk P, Laayoun A, Smith SS, Weitzman SA. DNA adduct 8-hydroxyl-2′-deoxyguanosine (8-hydroxyguanine) affects function of human DNA methyltransferase. Carcinogenesis. 1995;16(5):1253–5.
Maltseva DV, Baykov AA, Jeltsch A, Gromova ES. Impact of 7,8-dihydro-8-oxoguanine on methylation of the CpG site by Dnmt3a. Biochemistry. 2009;48(6):1361–8.
Campos ACE, Molognoni F, Melo FHM, Galdieri LC, Carneiro CRW, D'Almeida V, et al. Oxidative stress modulates DNA methylation during melanocyte Anchorage blockade associated with malignant transformation. Neoplasia. 2007;9(12):1111–21.
Vrijheid M, Martinez D, Manzanares S, Dadvand P, Schembari A, Rankin J, et al. Ambient air pollution and risk of congenital anomalies: a systematic review and meta-analysis. Environ Health Perspect. 2011;119(5):598–606.
Chen E, Zmirou-Navier D, Padilla C, Deguen S. Effects of air pollution on the risk of congenital anomalies: a systematic review and meta-analysis. Int J Environ Res Public Health. 2014;11(8):7642–68.
Gould T, Larson T, Stewart J, Kaufman JD, Slater D, McEwen N. A controlled inhalation diesel exhaust exposure facility with dynamic feedback control of PM concentration. Inhal Toxicol. 2008;20(1):49–52.
Weldy C, Luttrell IP, White CC, Morgan-Stevenson V, Cox DP, Carosino CM, et al. Glutathione (GSH) and the GSH synthesis gene Gclm modulate plasma redox and vascular responses to acute diesel exhaust inhalation in mice. Inhal Toxicol. 2013;25(8):444–54.
Liu Y, Yu M, Wu L, Chin MT. The bHLH transcription factor CHF1/Hey2 regulates susceptibility to apoptosis and heart failure after pressure overload. Am J Physiol Heart Circ Physiol. 2010;298(6):H2082–92.
Xiang F, Sakata Y, Cui L, Youngblood JM, Nakagami H, Liao JK, et al. Transcription factor CHF1/Hey2 suppresses cardiac hypertrophy through an inhibitory interaction with GATA4. Am J Physiol Heart Circ Physiol. 2009;290(5):H1997–2006.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, et al. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.
McKenna A, Hanna M, Banks E, Sivanchenko A, Cibulskis K, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.
DePristo M, Banks E, Poplin R, Garimella KV, Maguire JR, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43(5):491–8.
Trapnell C, Hendrickson DG, Sauvageau M, Goff L, Rinn JL, et al. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat Biotech. 2013;31:46–53.
Liao Y, Smyth GK, Shi W. The subread aligner: fast, accurate and scalable read mapping by seed-and-vote. Nucleic Acids Res. 2013;41(10):e108.
Wu D, Lim E, Vaillant F, Asselin-Labat ML, Visvader JE, Smyth GK. ROAST: rotation gene set tests for complex microarray experiments. Bioinformatics. 2010;26(17):2176–82.
Bolger A, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Krueger F, Andrews SR. Bismark: a flexible aligner and methylation caller for bisulfite-Seq applications. Bioinformatics. 2011;27(11):1571–2.
Hansen K, Langmead B, Irizarry RA. BSmooth: from whole genome bisulfite sequencing reads to differentially methylated regions. Genome Biol. 2012;13(10):R83.
Feng H, Conneely KN, Wu H. A Bayesian hierarchical model to detect differentially methylated loci from single nucleotide resolution sequencing data. Nucleic Acids Res. 2014;42(8):e69.
We thank Mr. James Stewart for expert technical assistance with diesel exhaust exposures.
JMG was supported by NIH training grant T32 ES007032. TKB, WMC, MTC and the overall project were supported by NIH R21 grant 1R21ES023015-01A1 and American Heart Association grant 15GRNT25390012 to MTC. JWM, TKB and the exposure and genomics aspects of the project were also supported by an NIH grant for the University of Washington Interdisciplinary Center for Exposures, Diseases, Genomics & Environment, P30ES07033.
Availability of data and materials
The datasets generated for RNA-sequencing and RRBS are available in the NCBI GEO Database under the SuperSeries accession number GSE110793.
Ethics approval and consent to participate
All animal studies were approved by the University of Washington Institutional Animal Care and Use Committee under protocol 4134–01. All other methods and experimental protocols were carried out under Institutional, State and Federal guidelines with approval from the appropriate University of Washington Environmental Health and Safety officials overseeing lab safety, chemical safety and biohazards.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.