Robust validation of methylation levels association at CPT1A locus with lipid plasma levels

There is increasing enthusiasm regarding the use of bio-banked whole blood DNA as a model to discover methylation marks associated with biological phenotypes and generate novel mechanistic hypotheses (1Murphy T.M. Mill J. Epigenetics in health and disease: heralding the EWAS era.The Lancet. 2014; : 60269-60275Google Scholar, 2Osório J. Obesity: Looking at the epigenetic link between obesity and its consequences-the promise of EWAS.Nat. Rev. Endocrinol. 2014; 10: 249Crossref PubMed Scopus (4) Google Scholar, 3Callaway E. Epigenomics starts to make its mark.Nature. 2014; 508: 22Crossref PubMed Scopus (22) Google Scholar). DNA methylation has a critical role in cell functions and is cell-type specific. Such cell specificity makes DNA methylation particularly challenging for epidemiological epigenetic investigations because disease relevant cell types might not be accessible due to practical issues such as availability, ethics, and cost associated with more complex specimen collection. Recent work suggests that agnostic methylation-wide association scan (MWAS) in peripheral blood can reflect phenotype-associated methylation marks in other tissues and cell types, with effects detected in established effector cells much stronger than effects detected in blood (4Dick K.J. Nelson C.P. Tsaprouni L. Sandling J.K. Aïssi D. Wahl S. Meduri E. Morange P.E. Gagnon F. Grallert H. et al.DNA methylation and body-mass index: a genome-wide analysis.The Lancet. 2014; : 62674-62682Google Scholar). These observations suggest that marks detected in blood are associated with functions in effector cells. The Illumina HumanMethylation450 (HM450K) array is a robust assay to measure DNA methylation across the genome (4Dick K.J. Nelson C.P. Tsaprouni L. Sandling J.K. Aïssi D. Wahl S. Meduri E. Morange P.E. Gagnon F. Grallert H. et al.DNA methylation and body-mass index: a genome-wide analysis.The Lancet. 2014; : 62674-62682Google Scholar, 5Zeilinger S. Kühnel B. Klopp N. Baurecht H. Kleinschmidt A. Gieger C. Weidinger S. Lattka E. Adamski J. Peters A. et al.Tobacco smoking leads to extensive genome-wide changes in DNA methylation.PLoS ONE. 2013; 8: e63812Crossref PubMed Scopus (517) Google Scholar, 6Bell J.T. Tsai P-C. Yang T-P. Pidsley R. Nisbet J. Glass D. Mangino M. Zhai G. Zhang F. Valdes A. et al.Epigenome-wide scans identify differentially methylated regions for age and age-related phenotypes in a healthy ageing population.PLoS Genet. 2012; 8: e1002629Crossref PubMed Scopus (514) Google Scholar, 7Dayeh T. Volkov P. Salö S. Hall E. Nilsson E. Olsson A.H. Kirkpatrick C.L. Wollheim C.B. Eliasson L. Rönn T. et al.Genome-wide DNA methylation analysis of human pancreatic islets from type 2 diabetic and non-diabetic donors identifies candidate genes that influence insulin secretion.PLoS Genet. 2014; 10: e1004160Crossref PubMed Scopus (324) Google Scholar). For any high-throughput technologies, and in particular for a novel assay such as the HM450K, rigorous quality control procedures are warranted and robustness of findings must be validated through independent replication to avoid reporting spurious associations. In the current issue of the Journal, Frazier-Wood et al. (8Frazier-Wood A.C. Aslibekyan S. Absher D.M. Hopkins P.H. Sha J. Tsai M.Y. Tiwari H.K. Waite L.L. Zhi D. Arnett D.K. Methylation at CPT1A locus is associated with lipoprotein subfraction profiles.J. Lipid Res. 2014; 55: 1324-1330Abstract Full Text Full Text PDF PubMed Scopus (55) Google Scholar) reported the novel findings of significant negative correlations between methylation levels at two CpG sites in the CPT1A locus and plasma levels of VLDL and LDL. Methylation levels were assessed in CD4+ T-cells isolated from peripheral blood DNA using the HM450K array. Given that no independent study samples were available for replication, to circumvent this challenge, the authors adopted an internal validation method by splitting the whole sample into “discovery and replication subsamples”. This strategy provides arguments in favor of the discovered associations but does not provide evidence of robustness against spurious findings due to sampling or confounding biases or any other undetected biases present in the study sample. A robust and thorough validation strategy implies the use of independent study samples and variation in the study designs (9Rosenbaum P.R. Replicating effects and biases.Am. Stat. 2001; 55: 223-227Crossref Scopus (55) Google Scholar, 10Kraft P. Zeggini E. Ioannidis J.P.A. Replication in genome-wide association studies.Stat. Sci. 2009; 24: 561-573Crossref PubMed Scopus (193) Google Scholar). The validation phase is of particular importance in MWAS, as this technique is particularly subjected to confounders (3Callaway E. Epigenomics starts to make its mark.Nature. 2014; 508: 22Crossref PubMed Scopus (22) Google Scholar). Thus, we undertook to test for associations the two CTP1A CpG sites found associated with lipid-related traits by Frazier-Wood et al. using two independent study samples with considerable variations in their respective study design and with the design of the Frazier-Wood study. The studies had differences in sampling scheme, DNA methylation specimen, and array preprocessing approaches. The notable differences in the design and sample characteristics between the three studies are shown in Table 1. Most notable is the method for lipid measurement, nuclear magnetic resonance spectroscopy in Frazier-Wood et al. and spectrophotometry in our studies. In addition to sampling variation, and of particular interest for MWAS studies, Frazier-Wood et al. assessed DNA methylation in isolated CD4+ T-cells, while we assessed methylation in peripheral whole blood, which includes CD4+ T-cell (<30%) and several other leukocyte subtypes. Finally, different normalization procedures were used: we applied the SWAN methodology (4Dick K.J. Nelson C.P. Tsaprouni L. Sandling J.K. Aïssi D. Wahl S. Meduri E. Morange P.E. Gagnon F. Grallert H. et al.DNA methylation and body-mass index: a genome-wide analysis.The Lancet. 2014; : 62674-62682Google Scholar, 11Maksimovic J. Gordon L. Oshlack A. SWAN: Subset-quantile within array normalization for Illumina Infinium HumanMethylation450 BeadChips.Genome Biol. 2012; 13: R44Crossref PubMed Scopus (558) Google Scholar) to globally normalize β values from the Infinium I and II probes, while separate normalization by probe type was applied by Frazier-Wood et al.TABLE 1Main design and sample characteristics of the three MWAS studies on lipidsStudy nameMARTHAF5L-PedigreesGOLDEN (Frazier-Wood et al.)Study designUnrelated individualsExtended pedigreesExtended pedigreesSubjects OriginCaucasians from Marseille area (South of France)French-Canadians from Ottawa area (Canada)European descent from Minneapolis (Minnesota) and Salt Lake City (Utah)DiscoveryValidationN327199663331Age44.1 (14.23)39.6 (16.9)48.6 (16.4)47.7 (16.6)Sex (% male)21.746.74749.2Total cholesterol5.452 (1.019) (g/L)4.896 (1.079) (g/L)NANAHDL-cholesterol1.476 (0.435) (g/L)1.359 (0.353) (g/L)40.0 (5.6) (nmol/L)37.0 (5.8) (nmol/L)LDL-cholesterol3.647 (0.980) (g/L)3.111 (0.901) (g/L)1393.8 (460.0) (nmol/L)1369.5 (1369.5) (nmol/L)Triglycerides1.058 (0.772) (mmol/L)1.487 (0.905) (mmol/L)NANALipid measurement technologySpectrophotometry except for LDL that was derived from the Friedewald's formulaNuclear Magnetic Resonance spectroscopyBlood collectionFastingFasting and non smokingFastingDNA specimenWhole bloodIsolated CD4+ T-cellsMedicationNo exclusionExclusion if on medicationAsked to discontinue the use of lipid lowering drugs and over-the-counter medication that could affect lipid levels.HumanMethylation450k NormalizationNoobaTriche et al. 2013. Low-level processing of Illumina Infinium DNA Methylation BeadArrays. Nucl. Acids Res. 41: e90. and SWANbRef. 11.Separately normalized probes from the Infinium I and II using ComBatcJohnson et al. 2007. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 8: 118–127.AdjustmentAge, sex, batch effect, chip effect, cell type composition, dyslipidemiaAge, sex, batch effect, chip effect, cell type composition,dIn MARTHA, specific measured biological counts of lymphocytes, monocytes, neutrophils, eosinophils and basophils were used to characterize leukocytes composition. In F5L-pedigrees, adjustment for cell type composition was handled by the methods described in Houseman et al. (BMC Bioinformatics 2012;13:86) family structureAge, sex, study site, T-cell purity (based on the first 4 principal components), family structurea Triche et al. 2013. Low-level processing of Illumina Infinium DNA Methylation BeadArrays. Nucl. Acids Res. 41: e90.b Ref. 11.c Johnson et al. 2007. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 8: 118–127.d In MARTHA, specific measured biological counts of lymphocytes, monocytes, neutrophils, eosinophils and basophils were used to characterize leukocytes composition. In F5L-pedigrees, adjustment for cell type composition was handled by the methods described in Houseman et al. (BMC Bioinformatics 2012;13:86) Open table in a new tab Despite the nontrivial differences between these studies, we observed strong statistical evidence for a negative association between the two CTP1A CpG sites (cg00574958 and cg17058475) identified by Frazier-Wood et al. and plasma levels of both LDL and triglycerides (TG) in two independent studies, the MARTHA (4Dick K.J. Nelson C.P. Tsaprouni L. Sandling J.K. Aïssi D. Wahl S. Meduri E. Morange P.E. Gagnon F. Grallert H. et al.DNA methylation and body-mass index: a genome-wide analysis.The Lancet. 2014; : 62674-62682Google Scholar, 12Oudot-Mellakh T. Cohen W. Germain M. Saut N. Kallel C. Zelenika D. Lathrop M. Trégouët D-A. Morange P-E. Genome wide association study for plasma levels of natural anticoagulant inhibitors and protein C anticoagulant pathway: the MARTHA project.Br. J. Haematol. 2012; 157: 230-239Crossref PubMed Scopus (42) Google Scholar) and F5L-pedigree studies (13Antoni G. Morange P-E. Luo Y. Saut N. Burgos G. Heath S. Germain M. Biron-Andreani C. Schved J.F. Pernod G. et al.A multi-stage multi-design strategy provides strong evidence that the BAI3 locus is associated with early-onset venous thromboembolism.J. Thromb. Haemost. 2010; 8: 2671-2679Crossref PubMed Scopus (37) Google Scholar). In our two samples totaling 526 individuals, increased DNA methylation levels at CPTA1 CpG sites were associated with both decreased LDL and TG (Table 2). A 1% increase in cg00574958 DNA methylation levels was associated with a 0.057 ± 0.011 decrease in log TG levels (P = 5.71 10−8). Corresponding values for a 1% increase in cg17058475 levels were 0.030 ± 0.008 (P = 9.83 10−5).TABLE 2Association of cg00574958 and cg17058475 CPT1A CpG variability with plasma TG and LDL levels in the MARTHA and F5L-pedigreesTG (log)LDLcg00574958MARTHA−0.059 (0.013) P = 8.28 10−6−0.023 (0.040) P = 0.57F5L-pedigrees−0.054 (0.018) P = 3.28 10−3−0.046 (0.019) P = 0.12CombinedaResults of the MARTHA and F5L-pedigrees studies were combined into a random effect meta-analysis based on the inverse-variance weighting method.−0.057 (0.011) P = 5.71 10−8−0.038 (0.024) P = 0.11cg17058475MARTHA−0.025 (0.009) P = 8.86 10−3−0.051 (0.029) P = 8.57 10−2F5L-pedigrees−0.041 (0.013) P = 2.88 10−3−0.037 (0.022) P = 9.36 10−2Combined−0.030 (0.008) P = 9.83 10−5−0.042 (0.017) P = 1.7 10−2Association was tested using a linear regression model (mixed linear model in F5L-Pedigrees) where log(TG) (LDL, resp.) was the outcome and the CpG site the predictor variable. Analyses were adjusted for age, sex, cell type, batch and chip effects. Reported coefficients (standard error) represent the increase in outcome value associated with a 1% increase in CpG site variability. In MARTHA, TG and LDL phenotypes were measured in 327 and 180 individuals, respectively. In the F5L-pedigrees study, lipid phenotypes were measured in 199 individuals.a Results of the MARTHA and F5L-pedigrees studies were combined into a random effect meta-analysis based on the inverse-variance weighting method. Open table in a new tab Association was tested using a linear regression model (mixed linear model in F5L-Pedigrees) where log(TG) (LDL, resp.) was the outcome and the CpG site the predictor variable. Analyses were adjusted for age, sex, cell type, batch and chip effects. Reported coefficients (standard error) represent the increase in outcome value associated with a 1% increase in CpG site variability. In MARTHA, TG and LDL phenotypes were measured in 327 and 180 individuals, respectively. In the F5L-pedigrees study, lipid phenotypes were measured in 199 individuals. Of note, cg00574958 and cg17058475 were highly correlated (ρspearman = 0.67 in both studies, P < 10−16); adjusting for cg00574958 in the model abolished the effect observed for cg17058475 on log TG levels. Finally, after adjustment for key covariates (age, sex, BMI, cell type composition, batch, and chip effects), cg00574958 explained ∼4% of log TG plasma levels, both in MARTHA and F5L-pedigrees. Negative association was also observed between plasma LDL levels and cg17058475 (P = 1.7 10−2) but not with cg00574958 (P = 0.11). No association was observed with HDL-cholesterol levels (P = 0.96 for cg00574958 and P = 0.75 for cg17058475), nor with total cholesterol levels (P = 0.16 for cg00574958 and P = 0.53 for cg17058475). The CPT1A protein is essential for fatty acid oxidation (a multistep process that metabolizes fats and converts them into energy) and is expressed in the liver and glandular tissues (14Uhlen M. Oksvold P. Fagerberg L. Lundberg E. Jonasson K. Forsberg M. Zwahlen M. Kampf C. Wester K. Hober S. et al.Towards a knowledge-based Human Protein Atlas.Nat. Biotechnol. 2010; 28: 1248-1250Crossref PubMed Scopus (1696) Google Scholar). This pivotal role in fatty acid metabolism makes CPT1A DNA methylation marks relevant to many metabolic disorders (from lipids to glucose homeostasis). The lipid-related DNA methylation probes in this study (cg00574958 and cg17058475) are designated as falling in a single “CpG shore”, and are flanked by two CpG islands. Human ENCODE HM450K studies performed on over 40 cell lines suggest these two probes show more variable methylation levels than the two CpG islands that flank them. The uncoupled methylation levels at these probes versus the flanking islands suggest that the observed variation is more likely to be regulatory. This region also shows evidence of open chromatin through DNase I hypersensitivity assays (15Thurman R.E. Rynes E. Humbert R. Vierstra J. Maurano M.T. Haugen E. Sheffield N.C. Stergachis A.B Wang H. Vernot B. Garg K. et al.The accessible chromatin landscape of the human genome.Nature. 2012; 489: 75-82Crossref PubMed Scopus (1790) Google Scholar) and gene regulatory potential through chromatin immunoprecipitation sequencing of the epigenetic modification H3K27ac (16ENCODE Project Consortium Bernstein B.E. Birney E. Dunham I. Green E.D. Gunter C. Snyder M. An integrated encyclopedia of DNA elements in the human genome.Nature. 2012; 489: 57-74Crossref PubMed Scopus (11021) Google Scholar). More work is needed to understand the functional impact of DNA methylation on CPT1A gene regulation. Three important conclusions emerge from this validation study. First, despite limitations in the Frazier-Wood et al. replication approach, the published results are robust to variation in sample, study design, normalization procedures, and even DNA blood specimen type. Second, inter-individual variation in lipid-related traits appears to be under the influence of DNA methylation regulation at the CPT1A locus. This epidemiological evidence now requires technical validation and functional work to confirm that these methylation marks are causes rather than consequences of lipid levels variation. Given that DNA methylation marks are potentially reversible, evidence for their role in the regulation of such a key enzyme is of great interest as it could lead to new therapeutic approaches (e.g., drug and/or diet supplementation) to modulate CPT1A expression. Finally, and of major importance for MWAS studies, peripheral whole blood DNA methylation marks were detected in an enzyme gene expressed in the liver and glandular tissues, suggesting that such marks could serve as surrogates for methylation at more closely-related effector cells, such as hepatocytes. The latter adds to the recent paper by Dick et al. (4Dick K.J. Nelson C.P. Tsaprouni L. Sandling J.K. Aïssi D. Wahl S. Meduri E. Morange P.E. Gagnon F. Grallert H. et al.DNA methylation and body-mass index: a genome-wide analysis.The Lancet. 2014; : 62674-62682Google Scholar) also supporting the value of peripheral whole blood DNA methylation marks as biomarkers of methylation in other tissues. We thank Dr. Michael D. Wilson for his judicious comments on the manuscript and for the many fruitful discussions about epigenetic regulation.

Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.