- Research article
- Open Access
- Open Peer Review
The DNA methylation drift of the atherosclerotic aorta increases with lesion progression
BMC Medical Genomicsvolume 8, Article number: 7 (2015)
Atherosclerosis severity-independent alterations in DNA methylation, a reversible and highly regulated DNA modification, have been detected in aortic atheromas, thus supporting the hypothesis that epigenetic mechanisms participate in the pathogenesis of atherosclerosis. One yet unaddressed issue is whether the progression of atherosclerosis is associated with an increase in DNA methylation drift in the vascular tissue. The purpose of the study was to identify CpG methylation profiles that vary with the progression of atherosclerosis in the human aorta.
We interrogated a set of donor-matched atherosclerotic and normal aortic samples ranging from histological grade III to VII, with a high-density (>450,000 CpG sites) DNA methylation microarray.
We detected a correlation between histological grade and intra-pair differential methylation for 1,985 autosomal CpGs, the vast majority of which drifted towards hypermethylation with lesion progression. The identified CpG loci map to genes that are regulated by known critical transcription factors involved in atherosclerosis and participate in inflammatory and immune responses. Functional relevance was corroborated by crossing the DNA methylation profiles with expression data obtained in the same human aorta sample set, by a transcriptome-wide analysis of murine atherosclerotic aortas and from available public databases.
Our work identifies for the first time atherosclerosis progression-specific DNA methylation profiles in the vascular tissue. These findings provide potential novel markers of lesion severity and targets to counteract the progression of the atheroma.
DNA methylation is an important epigenetic mechanism of transcriptional regulation [1,2]. The DNA methylome undergoes programmed changes during cellular differentiation, but it can be also modified by exogenous stimuli such as the diet and environmental factors, some of which are of the same typology as known atherosclerosis risk factors [3,4]. Therefore, it has been proposed that the latter may act by imposing aberrant, proatherogenic DNA methylation patterns [5,6]. Indeed, differentially methylated sites mapping to genes participating in atherogenesis have been identified by candidate gene-based studies and in a recent epigenome-wide survey of normal and atherosclerotic human aortas, which showed a genome-wide increase in DNA methylation in atherosclerosis [7-9]. Accordingly, biochemical inhibition of DNA methyltransferase activity decreases vascular lesion size in mice [10,11]. Furthermore, epigenome-wide association studies (EWAS) have linked peripheral blood cell DNA methylation profiles of specific loci to hyperlipidemia, hyperglycemia and obesity, thus uncovering candidate circulating epigenetic markers of atherosclerosis and metabolic conditions predisposing to atherosclerosis [12-14].
By exploiting a high-coverage DNA methylation microarray  that determines the methylation status of more than 450,000 CpG sites in the human genome, we previously identified a set of CpG loci (dm-CpGs) that are differentially methylated between donor-matched atherosclerotic and normal human aortas . dm-CpG methylation profiles are independent of lesion histological grade in the grade III-VII range, suggesting that they are established relatively early in the evolution of the atheroma. One question left unanswered in our previous study is whether DNA methylation changes at any CpG loci during the progression of aortic atherosclerosis. To address that issue, we searched the original array data  for DNA methylation profiles that are significantly correlated with the histological grade of the atheroma.
Human vascular samples
The post mortem donor-matched atherosclerotic and non-atherosclerotic portions of human aortas (referred to as “donor-matched sample pairs”) were previously described .
Aortic samples were obtained at the Bellvitge Hospital and the Clinic Hospital in Barcelona, Spain according to a protocol approved by the ethics committee of the Bellvitge Hospital (authorization no. PR311/11) and following signed consent by the relatives. Post-mortem time was between 3 and 26 h. Atherosclerotic and normal (A and N, respectively) portions were macroscopically identified by a trained Pathology Technician. Approximately 3-cm long segments of each portion were transferred to RNAlater (Ambion) and refrigerated. A fragment was obtained from the middle portion of each segment (either A or N) and the severity of its histological atherosclerosis was scored on the AHA scale by a qualified Pathologist. The remaining samples were transferred to −80°C within 1 day of dissection. Atherosclerotic lesions and the underlying media were used in all the procedures described here. The relevant donor (age, sex) and sample information for the aortic specimens is shown in Additional file 1: Table S1 in Zaina et al. .
A diagram describing the sample origin and workflow is presented in Figure 1.
DNA methylation analysis
The microarray-based DNA methylation analysis was conducted with the Infinium HumanMethylation450 BeadChip (450 k array). DNA quality checks, bisulfite modification, hybridization, data normalization, statistical filtering and Beta value calculation were carried out as described elsewhere [15,16]. SNP-harboring probes (11,858) obtained from the list provided by Illumina (support.illumina.com/downloads/infinium_hd_methylation_snp_list.html) were excluded from the analysis. The array data for the samples presented here are extracted from our previously validated and published data set , freely available at the Gene Expression Omnibus (GEO) database: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=zfephwswcyoeotw&acc=GSE46401. The methylation level for each cytosine was expressed as a Beta value calculated as the fluorescence intensity ratio of the methylated to unmethylated versions of the probes. Beta values ranged between 0 (unmethylated) and 1 (methylated). The annotation relating to CGIs uses the following categorization: “shore”, each of the 2 kb-sequences flanking a CGI; “shelf”, each of the 2 kb-sequences next to a shore; “open sea”, DNA not included in any of the previous sequences or in CGIs2. TSS200 orTSS1500 indicate the region between position −200 bp or −1,500 bp from the Transcription Start Site (TSS), respectively. Functional gene annotation was performed with the DAVID database (http://david.abcc.ncifcrf.gov). Clustering was performed by Ward’s hierarchical clustering method, using simultaneously samples and CpGs.
Whole genome bisulfite sequencing data of one atherosclerotic (grade VII)/normal human aortic tissue pair were previously published . The complete DNA methylation sequence at the single nucleotide resolution is freely available at the GEO database: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=zfephwswcyoeotw&acc=GSE46401.
RNA-seq of mouse aortas
The mouse experimental protocol was approved by the ethics committee of the Department of Medical Sciences, University of Guanajuato, León, Mexico. Apolipoprotein E-null (APOE-null) mice homozygous for the disrupted Apoe allele  and sex- and age-matched WT siblings were obtained by mating heterozygous mutants that have been maintained in the C57BL/6 background for >100 generations. Mice were fed a normal chow diet (LabDiet 5001) until sacrifice at the age of 1.6 years. After culling by decapitation under anesthesia with Isoflurane (Forane, Baxter), tissues were dissected and kept in RNAlater (Ambion) at −80°C. Aortas were flushed in situ with PBS by cardiac puncture after removal of the thoracic and abdominal organs, and carefully cleaned of any fat tissue under the microscope. Trizol-extracted RNA Integrity Number (RIN) and concentration was measured using the Agilent 2100 Bioanalyzer system. After verifying that each sample presented a RIN number > = 7, Illumina RNA-seq libraries were prepared using the TruSeq RNA Sample Preparation Kit, following the vendor instructions and the resulting cDNA libraries were sequenced using the Illumina Genome Analyzer IIx platform with a Paired End configuration of 144 cycles (72 bp read length pairs). The resulting reads were mapped using the BWA aligner (version 0.7.5a-r405) with default parameters  to a compilation of CDS sequences from the mouse reference genome (Mus musculus strain C57BL/6 J GRCm38.p2 C57BL/6 J) downloaded from ENSEMBL database using BioMart (http://www.ensembl.org/biomart). Duplicates were discarded by using the Picard MarkDuplicate tool (http://picard.sourceforge.net) and from the filtered alignment, the read count for each gene was calculated. The gene differential expression analysis between the wild type and mutant mouse aortas was performed using the NOIseq R package  with TMM normalization and a cutoff line of 0.9 (ranking score) for differential expression significance, which means that a certain gene is about 20 times more likely to be differentially expressed than non-differentially expressed. RNA-seq data are available in the NCBI BioProject database (Project accession PRJNA262445).
In the case of 450 k arrays, Beta values were converted to M values by a Logit transformation . Intra-pair differences in M values were correlated (Spearman’s correlation test) with the histological grade (III-VII range) converted to numeral values. The correlation test was used as a first filter to identify CpGs that exhibit a significant methylation change with histological grade and to determine the direction of that change (i.e. progressive hypermethylation or hypomethylation). The significance of the identified CpGs was corroborated by ANOVA adjusted for age, sex and post mortem time. The significance threshold in whole 450 k array data analysis was a Bonferroni-corrected p < 10−7, corresponding to an uncorrected p < 0.05. As for the statistical power, the paired analysis of the same sample (n = 15 pairs) yielded a sizeable number of CpG sites with significant differential methylation after Bonferroni correction . Therefore, we assumed that the sample size would provide the sufficient power for the statistical analysis presented here. As Beta values are a more intuitive representation of methylation profiles, they will be used to describe the data in the Results and Discussion sections. Percentages were compared by using the Chi-square test. The Pearson’s correlation test was used to compare RNA-seq-derived gene expression level with the direction of methylation change (i.e. the correlation r between the delta-M and histological grade). All tests were performed using the STATISTICA software.
Aberrant DNA methylation profile of atherosclerotic lesions becomes more frequent with histological grade
We sought DNA methylation profiles that changed with the progression of atherosclerosis. To this end, we conducted a correlation analysis between the methylation M value difference within donor-matched atherosclerotic (A)/normal (N) aortic sample pairs (n = 15) and the respective lesion grade for the >450,000 CpG sites profiled by the Infinium HumanMethylation450 BeadChip (450 k array). Aortic sample histological grade ranged from III to VII. Significantly correlated CpG sites will be referred to as grade-CpGs. We set the following criteria to assign a grade-CpG status: 1) absolute difference >0.05 between the average intra-pair Delta-Betas of the extreme histological grade subsets - i.e. grade III and grade VII. 2) Significant (p < 0.05) Spearman’s correlation rho between histological grade and intra-pair Delta-M, and no significant correlation with donor’s age, sex and post-mortem time. The significant of the correlation test result was further corroborated by submitting the identified CpG set to ANOVA with donor’s age, sex and post-mortem time as covariates (significance threshold p < 0.05). 3) No significant difference between M values of grade VII and grade III N samples, to exclude possible confounding factors that might over- or underestimate the effect of histological grade.
We identified 1,985 autosomal CpG sites (corresponding to 1,206 genes) that matched the grade-CpG criteria (Additional file 1: Table S1). Spearman’s correlation r-values ranged from 0.863 (p = 3.4×10−5) to −0.558 (p = 0.031), thus no grade-CpG displayed a significance level below the Bonferroni-corrected p-value. Yet, clustering analysis of the 1,985 grade-CpGs showed a clear grouping by low (III-IV) and high (V-VII) histological grade, with a tendency for hypermethylation in high-histological grade samples (Figure 2A). Accordingly, the vast majority of the grade-CpG methylation profiles (1,631 or 82.2%) correlated positively with atherosclerosis grade, suggesting a gradient change of DNA methylation with disease progression, predominantly towards hypermethylation in advanced lesions (Figure 2B). On average, grade-CpGs displayed weak differential DNA methylation (average absolute grade III versus grade VII A sample Delta-Beta = 0.069; range: 0.050 to 0.522). Notably, 111 genes contained both grade-CpGs and the previously described histological grade-independent dm-CpGs . Thus, grade-CpGs reinforced the DNA methylation drift represented by dm-CpG. As examples, two relevant genes are shown, one involved in vascular smooth muscle cell (VSMC) proliferation (platelet-derived growth factor alpha polypeptide, PDGFA)  and the other participating in the induction of repressive chromatin, in adipose tissue homeostasis and identified as part of the gene expression signature of stable atherosclerotic plaques (PR domain containing 16, PRDM16) [22-24] (Figure 2C).
As for genomic distribution, grade-CpGs mapped preferentially to gene bodies and were underrepresented in promoters (Chi-square test, p = 0.029, compared with the probe distribution on the 450 k array), whereas their distribution relative to CGIs was only marginally different from expected (p = 0.059) (Figure 2D). When grade-CpGs were grouped according to the direction of methylation change with histological grade, noticeable differences in the genomic distribution were observed. First, CGIs were significantly underrepresented within the hypomethylated grade-CpG set in comparison with the hypermethylated counterpart (6.7% and 21.5%, respectively, Chi-square p = 2.39 × 10−5), whereas only a marginal difference in distribution was observed relative to gene compartments (p = 0.07). Second, CpGs undergoing hypomethylation were initially highly methylated - i.e. in the N (baseline) sample, while the hypermethylated counterparts were evenly distributed across the 0–1 Beta range in the N sample (average Beta 0.63 ± 0.20 and 0.48 ± 0.26 for the hypomethylated and hypermethylated grade-CpGs, respectively; t-test p = 5.35 × 10−31) (Additional file 2: Figure S1A,B). Accordingly, when grade-CpGs of either set were further fractioned by N sample methylation level (high: Beta > 0.75; intermediate: 0.75 ≥ Beta ≥ 0.25; low: Beta < 0.25), a significant underrepresentation of the low methylation fraction within hypomethylated grade-CpGs was observed (7.6% and 26.3% in hypometylated and hypermethylated grade-CpGs, respectively; Chi-square p = 2.3 × 10−5) (Additional file 2: Figure S1C). Thus, the histological grade-related hypomethylation and hypermethylation differed, as the former targeted highly methylated CpGs preferentially outside CGIs, in comparison with the latter.
Functionally, grade-CpG-harboring genes were significantly enriched in biological processes involving immune and inflammatory response genes (Fc-gamma receptor-mediated phagocytosis, peroxisome proliferator-activated receptor (PPAR) signaling, T-cell receptor signaling, immune cell-specific expression; Fisher’s exact test, p < 0.05) (Table 1). Noticeably, genes regulated by the transcriptional factors ZEB1 (AREB6) and PPAR-gamma, involved in the cellular response to oxidized lipoproteins and macrophage activation [25,26], respectively, were significantly enriched (949 and 924 genes, p = 4.2 × 10−10 and p = 2.4 × 10−10) (Table 1). Eighteen grade-CpGs mapped to 9 genes previously identified as differentially methylated in atherosclerosis by candidate gene-based or epigenomics studies (Table 2). Furthermore, 24 grade-CpGs were located within 250 kb from 16 cardiovascular disease-associated single nucleotide polymorphisms (SNPs) , suggesting an interplay between the risk genotype and differential DNA methylation as methylation quantitative trait loci (Additional file 3: Table S2) . When grouped according to the direction of methylation change with histological grade, hypermethylated grade-CpG-harboring genes were enriched for the same functional categories as the bulk of grade-CpGs, while no significant enrichment was observed for the hypomethylated counterpart (not shown). Based on our data, we cannot conclude whether the latter result reflects an underlying biological property or is due to a sample size effect.
We performed a technical validation of 450 k data by comparing the promoter grade-CpG methylation profiles with our previous whole genome bisulfite sequencing (WGBS) data for one of the donor-paired aortic samples presented here (sample 93, displaying a grade VII lesion) . Eighteen promoter grade-CpGs mapped to promoter differentially methylated regions (DMRs) identified by WGBS. The correlation between the average intra-pair Delta-Beta of grade VII aortic samples (n = 9) and the intra-pair difference in methylation for the single WGBS aortic sample was significant (r = 0.540, p = 0.021, n = 18), thus ruling out technique-related confounding factors. As for biological validation, the same 450 k array data from which the grade-CpG profiles were obtained, were previously validated in 15 discovery and 24 validation paired aortic samples , including the grade-CpG-harboring genes PDGFA and C9orf3/MIR23b.
A grade-CpG subset is differentially expressed in human atheromas and in a heterologous model of advanced atherosclerosis
To further corroborate the pathobiological relevance of grade-CpGs, we crossed their methylation profiles with expression data obtained experimentally or from public databases. For the PDGFA and C9orf3/MIR23b loci, we previously showed that differential methylation coincided with significant changes in expression in aortic samples . Furthermore, 105 grade-CpG-harboring genes overlapped with a human atherosclerotic plaque-specific genes expression signature , and 87 overlapped with the list of atherosclerosis-related genes in the HuGENet™ database. The full list of the overlapping genes is shown in Additional file 4: Table S3. Overlapping genes are involved in relevant functions such as lipid homeostasis (ALOX12, CETP, LPL, LRP1, LRP5, PLA2G2A), chemokine signaling (CCR5), VSMC proliferation (PDGFA), atheroprotection (ESR1), local inflammatory response (VCAM1, C3, ITGB2) and T-cell receptor signaling (genes listed in Table 1).
Next, we asked whether any grade-CpG-harboring gene was differentially expressed in a mouse model of advanced atherosclerosis. We reasoned that a strong correlation between human aortic DNA methylation profiles and gene expression in a heterologous model will identify grade-CpGs with potential pathobiological relevance and suitability for comparative mouse model studies. To this end, we performed a RNA-seq-based gene expression analysis in aortic atherosclerotic lesions of one 1.6 years old male APOE-null mouse, a model of hyperlipidemia-induced atherosclerosis  and one congenic, sex- and age-matched WT sibling. As expected in a mutant mouse of advanced age, the thoracic aorta showed extensive atheromas (Additional file 2: Figure S2). The relevant descriptive statistics for the RNA-seq analysis of WT and mutant aortic samples is shown in Additional file 5: Table S4. After comparing the samples using the NOIseq R package, among the genes showing a >2-fold expression difference between mutant and control mice, 53 up-regulated and 40 down-regulated genes in the APOE-null mouse harbored grade-CpGs. A total of 99 grade-CpGs mapped to those 93 differentially expressed mouse genes (mgrade-CpGs). The mouse genes and corresponding grade-CpGs are listed in Additional file 6: Table S5. The 23 promoter mgrade-CpGs showed a significant inverse correlation between the histological grade-associated methylation change (i.e. the Spearman’s correlation rho between the methylation status and histological grade) in human aortas and the differential expression of the corresponding mouse genes - i.e. APOE-null relative to WT (Pearson’s r = −0.71, p = 1.5×10−4). Thus, the aberrant methylation of promoter grade-CpGs during atherosclerosis progression is associated with an expected expression change of the corresponding gene in a heterologous advanced lesion. A significant inverse correlation for open sea grade-CpG was also observed (r = −0.384, p = 0.009, n = 46), reflecting the abundance of promoter and 5′ gene portions (5′UTR, first exon) mapping to open sea in the mgrade-CpG set (21 out of 46, or 45.6%). The methylation status of gene-body mgrade-CpGs did not show any significant correlation with expression in accordance with the published evidence that intragenic methylation may enhance or repress transcription in a context-specific manner [34,35] nor did CpG island-, shore- or shelf-mapping mgrade-CpGs (data not shown).
As for gene function, 29 of the mgrade-CpG-harboring up-regulated genes in the APOE-null aorta were targets of the master regulator of inflammation nuclear factor kappa B (NFKB) , (p = 1.8 × 10−4). Interestingly, pathway analysis (BioCarta) revealed enrichment in local inflammatory response and transcription regulation by the arginine methyltransferase CARM1 genes (PRKAR1B, RARA) (p = 0.006 and p = 0.05, respectively). These results are consistent with the participation of CARM1 in NFKB transcriptional activation . Of the down-regulated genes, 32 were regulated by myocyte enhancer factor 2A (MEF2A), a transcription factor implicated in coronary artery disease (p = 3.5 × 10−7) . Notably, the observation that loss of MEF2A leads to VSMC dedifferentiation is consistent with the enrichment of its targets in the down-regulated gene set .
We detected a significant correlation between intra-aortic pair DNA methylation Delta-Beta and lesion histological grade for 1,985 autosomal CpGs (grade-CpGs). Although no grade-CpG lies below the multiple testing-corrected significance level, the results are of potential interest due the convergence of multiple independent supporting evidence. First, we observed a progressive increase in methylation of macrophage function and immunity genes with atherosclerosis progression, reflecting the well-characterized migration of immune and inflammatory cells into the vascular wall and lesion mass during the natural history of atherosclerosis . Second, the hypermethylation of the majority of grade-CpGs with lesion severity represents an exacerbation of the hypermethylated profiles of previously described atherosclerosis-specific dm-CpGs . dm-Cpgs are lesion progression-independent profiles and therefore represent an epigenetic drift of the vascular genome that is already present in low-grade lesions and is amplified by grade-CpGs during lesion progression. Although overlapping in selected loci, the bulk of dm-CpGs and grade-CpGs map to functionally distinct genes, i.e. to SMC-related and immune/inflammatory genes, respectively, pointing to the specific regulation of different pathways in early and progressing atherosclerosis . Noticeably, the hypermethylation of the majority of grade-CpGs in high-grade lesions matches the observation that high expression of the DNA methyltransferase DNMT1 is a signature of inflamed, advanced human plaques . Interestingly, histological grade-related hypermethylation indiscriminately targeted low, intermediate and high methylation CpGs, whereas hypomethylation showed a preference for initially (N sample baseline) highly methylated CpGs in CG-poor regions. Further studies will be necessary to fully elucidate the mechanisms and the biological significance of this observation. Third, the genomic distribution of grade-CpGs generally reflects the one of dm-CpGs, suggesting a common instructive mechanism of DNA methylation drift, particularly the preference for gene-body hypermethylation, despite the functional differences mentioned above. Fourth, several grade-CpGs map to previously identified differentially methylated genes in atherosclerosis. In particular, both a previous epigenomics study carried out in peripheral blood samples  and our work in the aorta , identify differential methylation in the C1QL4, CTNNA3 and IMMT genes. These loci undergo an epigenetic drift towards hypermethylation in independent cardiovascular sample sets representing different tissues and are therefore potential important functional markers. Lastly, we crossed the methylation profiles with available public human databases and transcriptome-wide expression data for a mouse aortic lesion. Although the expression data lack replicates, we used an analysis package (NOIseq) that was developed for differential expression analysis conditions with no replicates as a non-parametric algorithm that evaluates the significance of gene expression changes by comparison against data noise simulation . Hence, the strong inverse association between the methylation status of human promoter grade-CpGs and the expression levels in that heterologous model strongly suggest a general pathobiological relevance for grade-CpGs.
We provide for the first time a list of CpG loci in the vascular lesion genome that undergo a DNA methylation drift with the progression of the lesion. The identified grade-CpGs represent promising markers of atheroma progression or targets to counteract the progressive DNA hypermethylation associated with the natural history of atherosclerosis.
Availability of supporting data
The DNA methylation array data have been deposited in the Gene Expression Omnibus (GEO) database: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=zfephwswcyoeotw&acc=GSE46401.
The RNA-seq data have been deposited in the NCBI BioProject database (www.ncbi.nlm.nih.gov/bioproject) with project accession: PRJNA262445.
Reik W. Stability and flexibility of epigenetic gene regulation in mammalian development. Nature. 2007;447:425–32.
Szulwach KE, Jin P. Integrating DNA methylation dynamics into a framework for understanding epigenetic codes. Bioessays. 2014;36:107–17.
Apostolou E, Hochedlinger K. Chromatin dynamics during cellular reprogramming. Nature. 2013;502:462–71.
Waterland RA, Jirtle RL. Transposable elements: targets for early nutritional effects on epigenetic gene regulation. Mol Cell Biol. 2003;23:5293–300.
Newman PE. Can reduced folic acid and vitamin B12 levels cause deficient DNA methylation producing mutations which initiate atherosclerosis? Med Hypotheses. 1999;53:421–4.
Ordovás JM, Smith CE. Epigenetics and cardiovascular disease. Nat Rev Cardiol. 2010;7:510–9.
Laukkanen MO, Mannermaa S, Hiltunen MO, Aittomäki S, Airenne K, Jänne J, et al. Local hypomethylation in atherosclerosis found in rabbit ec-sod gene. Arterioscler Thromb Vasc Biol. 1999;19:2171–8.
Zaina S. Unraveling the DNA methylome of atherosclerosis. Curr Opin Lipidol. 2014;25:148–53.
Zaina S, Heyn H, Carmona FJ, Varol N, Sayols S, Condom E, et al. A DNA methylation map of human atherosclerosis. Circ Cardiovasc Genet. 2014;7:692–700.
Dunn J, Qiu H, Kim S, Jjingo D, Hoffman R, Kim CW, et al. Flow-dependent epigenetic DNA methylation regulates endothelial gene expression and atherosclerosis. J Clin Invest. 2014;124:3187–99.
Cao Q, Wang X, Jia L, Mondal AK, Diallo A, Hawkins GA, et al. Inhibiting DNA Methylation by 5-aza-2′-deoxycytidine ameliorates atherosclerosis through suppressing macrophage inflammation. Endocrinology. 2014. [Epub ahead of print].
Irvin MR, Zhi D, Joehanes R, Mendelson M, Aslibekyan S, Claas SA, et al. Epigenome-wide association study of fasting blood lipids in the genetics of lipid lowering drugs and diet network study. Circulation. 2014. [Epub ahead of print].
Hidalgo B, Irvin MR, Sha J, Zhi D, Aslibekyan S, Absher D, et al. Epigenome-wide association study of fasting measures of glucose, insulin, and HOMA-IR in the genetics of lipid lowering drugs and diet network study. Diabetes. 2014;63:801–7.
Dick KJ, Nelson CP, Tsaprouni L, Sandling JK, Aïssi D, Wahl S, et al. DNA methylation and body-mass index: a genome-wide analysis. Lancet. 2014;383:1990–8.
Sandoval J, Heyn H, Moran S, Serra-Musach J, Pujana MA, Bibikova M, et al. Validation of a DNA methylation microarray for 450,000 CpG sites in the human genome. Epigenetics. 2011;6:692–702.
Heyn H, Li N, Ferreira HJ, Moran S, Pisano DG, Gomez A, et al. Distinct DNA methylomes of newborns and centenarians. Proc Natl Acad Sci U S A. 2012;109:10522–7.
Zhang SH, Reddick RL, Piedrahita JA, Maeda N. Spontaneous hypercholesterolemia and arterial lesions in mice lacking apolipoprotein E. Science. 1992;258:468–71.
Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25:1754–60.
Tarazona S, Garcia-Alcalde F, Dopazo J, Ferrer A, Conesa A. Differential expression in RNA-seq: a matter of depth. Genome Res. 2011;21:4436.
Du P, Zhang X, Huang CC, Jafari N, Kibbe WA, Hou L, et al. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinformatics. 2010;11:587.
George SJ, Williams A, Newby AC. An essential role for platelet-derived growth factor in neointima formation in human saphenous vein in vitro. Atherosclerosis. 1996;120:227–40.
Pinheiro I, Margueron R, Shukeir N, Eisold M, Fritzsch C, Richter FM, et al. Prdm3 and Prdm16 are H3K9me1 methyltransferases required for mammalian heterochromatin integrity. Cell. 2012;150:948–60.
Hondares E, Rosell M, Díaz-Delfín J, Olmos Y, Monsalve M, Iglesias R, et al. Peroxisome proliferator-activated receptor α (PPARα) induces PPARγ coactivator 1α (PGC-1α) gene expression and contributes to thermogenic activation of brown fat: involvement of PRDM16. J Biol Chem. 2011;286:43112–22.
Puig O, Yuan J, Stepaniants S, Zieba R, Zycband E, Morris M, et al. A gene expression signature that classifies human atherosclerotic plaque by relative inflammation status. Circ Cardiovasc Genet. 2011;4:595–604.
Karagiannis GS, Weile J, Bader GD, Minta J. Integrative pathway dissection of molecular mechanisms of moxLDL-induced vascular smooth muscle phenotype transformation. BMC Cardiovasc Disord. 2013;13:4.
Ricote M, Li AC, Willson TM, Kelly CJ, Glass CK. The peroxisome proliferator-activated receptor-gamma is a negative regulator of macrophage activation. Nature. 1998;391:79–82.
Stylianou IM, Bauer RC, Reilly MP, Rader DJ. Genetic basis of atherosclerosis: insights from mice and humans. Circ Res. 2012;110:337–55.
Gibbs JR, van der Brug MP, Hernandez DG, Traynor BJ, Nalls MA, Lai SL, et al. Abundant quantitative trait loci exist for DNA methylation and gene expression in human brain. PLoS Genet. 2010;6:e1000952.
Sharma P, Kumar J, Garg G, Kumar A, Patowary A, Karthikeyan G, et al. Detection of altered global DNA methylation in coronary artery disease patients. DNA Cell Biol. 2008;27:357–65.
Huang YS, Zhi YF, Wang SR. Hypermethylation of estrogen receptor-alpha gene in atheromatosis patients and its correlation with homocysteine. Pathophysiology. 2009;16:259–65.
Zawada AM, Rogacev KS, Hummel B, Grün OS, Friedrich A, Rotter B, et al. SuperTAG methylation-specific digital karyotyping reveals uremia-induced epigenetic dysregulation of atherosclerosis-related genes. Circ Cardiovasc Genet. 2012;5:611–20.
Zhang D, Chen Y, Xie X, Liu J, Wang Q, Kong W, et al. Homocysteine activates vascular smooth muscle cells by DNA demethylation of platelet-derived growth factor in endothelial cells. J Mol Cell Cardiol. 2012;53:487–96.
Zhu S, Goldschmidt-Clermont PJ, Dong C. Inactivation of monocarboxylate transporter MCT3 by DNA methylation in atherosclerosis. Circulation. 2005;112:1353–61.
Lorincz MC, Dickerson DR, Schmitt M, Groudine M. Intragenic DNA methylation alters chromatin structure and elongation efficiency in mammalian cells. Nat Struct Mol Biol. 2004;11:1068–75.
Flanagan JM, Wild L. An epigenetic role for noncoding RNAs and intragenic DNA methylation. Genome Biol. 2007;8:307.
Shimizu H, Mitomo K, Watanabe T, Okamoto S, Yamamoto K. Involvement of a NF-kappa B-like transcription factor in the activation of the interleukin-6 gene by inflammatory lymphokines. Mol Cell Biol. 1990;10:561–8.
Covic M, Hassa PO, Saccani S, Buerki C, Meier NI, Lombardi C, et al. Arginine methyltransferase CARM1 is a promoter-specific regulator of NF-kappaB-dependent gene expression. EMBO J. 2005;24:85–96.
Wang L, Fan C, Topol SE, Topol EJ, Wang Q. Mutation of MEF2A in an inherited disorder with features of coronary artery disease. Science. 2003;302:1578–81.
Zhao W, Zhao SP, Peng DQ. The effects of myocyte enhancer factor 2A gene on the proliferation, migration and phenotype of vascular smooth muscle cells. Cell Biochem Funct. 2012;30:108–13.
Ross R. Atherosclerosis - an inflammatory disease. N Engl J Med. 1999;340:115–26.
We thank Esteve González (Department of Pathology, Bellvitge University Hospital) and Abel Muñoz (Department of Anatomy and Pathology, Hospital Clinic) for technical assistance. We also thank the “Unidad Universitaria de Secuenciacion Masiva de ADN” and the “Unidad Universitaria de Apoyo Bioinformatico” for the support in generating the Illumina RNA-seq data and in bioinformatics analysis, respectively; Veronica Jimenez-Jancinto and Jerome Verleyen for data organization and submission to the NCBI repository.
This work was supported by the European Research Council (ERC) grant EPINORC under agreement n° 268626, the MICINN Project – SAF2011-22803, the Cellex Foundation, the Botin Foundation, the European Community’s Seventh Framework Programme (FP7/2007-2013) from grant HEALTH-F5-2011-282510 – BLUEPRINT and the Health and Science Departments of the Generalitat de Catalunya. S.Z. was supported by a CONACYT (Mexico) Sabbatical Fellowship no. 166058 and by the University of Guanajuato (DAIP) “2014 Interdisciplinary Projects” grant no. 420. H.H. is a Juan de la Cierva Postdoctoral Researcher. M.E. is an Institució Catalana de Recerca i Estudis Avançats (ICREA) Research Professor.
The authors declare that they have no competing interests.
MPV-M co-designed the mouse experiments, performed mouse dissections and RNA extraction, participated in the RNA-seq and mouse bioinformatics work; SZ co-designed the human and mouse studies, extracted DNA, analyzed the normalized array data and wrote the manuscript draft; HH and AG participated in the human bioinformatics analysis; FJC and NV performed pyrosequencing; SS performed the bioinformatics analysis of WGBS data; EC organized the collection of human samples and performed the histological grading of human lesions; JR-R organized the collection of human samples; SM performed the human array work; GL critically revised the manuscript; DR-R participated in the mouse molecular work; GL-G, MR-N and CR-M organized mouse breeding and genotyping; AS-F performed the RNA-seq experiments and the majority of the bioinformatics analysis of RNA-seq data; ME co-designed the study. All authors read and approved the manuscript.
María del Pilar Valencia-Morales and Silvio Zaina contributed equally to this work.
List of autosomal grade-CpGs.
Differential distribution of grade-CpGs that undergo hypermethylation or hypomethylation with histological grade, among N sample Beta classes. A,B: N sample Beta of hypermethylated and hypomethylated CpGs (rho>0 and rho<0, respectively). CpGs are ordered for decreasing methylation/grade Spearman’s rho (left to right). C, Distribution of hypermethylated (solid bars) and hypomethylated (open bars) grade-CpGs among high, intermediate and low methylation fractions (left to right). Notice the near-symmetrical distribution of hypermethylated grade-CpGs and the relative underrepresentation of the low methylation fraction among the hypomethylated counterpart. Figure S2: Aortic atherosclerosis in the APOE-null mouse aorta analyzed by RNA-seq. The heart and thoracic aorta of sex-matched, 1.6 years old APOE-null (left) and wt (right) mice are shown. Notice the abundant lesions in the aortic root and arch (arrows) and in the descending aorta in the APOE-null mouse.
List of CVD-related SNPs mapping close to grade-CpGs.
Grade-CpG-harboring genes overlapping with publicly available CVD gene expression signatures.
Grade-CpG-harboring genes that are differentially expressed in APOE-null and WT mice.