Genetic variation influencing DNA methylation provides insights into molecular mechanisms regulating genomic function
Hawe JS., Wilson R., Schmid KT., Zhou L., Lakshmanan LN., Lehne BC., Kühnel B., Scott WR., Wielscher M., Yew YW., Baumbach C., Lee DP., Marouli E., Bernard M., Pfeiffer L., Matías-García PR., Autio MI., Bourgeois S., Herder C., Karhunen V., Meitinger T., Prokisch H., Rathmann W., Roden M., Sebert S., Shin J., Strauch K., Zhang W., Tan WLW., Hauck SM., Merl-Pham J., Grallert H., Barbosa EGV., Ahmadi KR., Ainali C., Barrett A., Bataille V., Bell JT., Buil A., Dermitzakis ET., Dimas AS., Durbin R., Glass D., Grundberg E., Hassanali N., Hedman ÅK., Ingle C., Knowles D., Krestyaninova M., Lindgren CM., Lowe CE., McCarthy MI., Meduri E., di Meglio P., Min JL., Montgomery SB., Nestle FO., Nica AC., Nisbet J., O’Rahilly S., Parts L., Potter S., Sandling J., Sekowska M., Shin SY., Small KS., Soranzo N., Spector TD., Surdulescu G., Travers ME., Tsaprouni L., Tsoka S., Wilk A., Yang TP., Zondervan KT., Illig T., Peters A., Paus T., Pausova Z., Deloukas P., Foo RSY., Jarvelin MR., Kooner JS., Loh M., Heinig M., Gieger C., Waldenberger M., Chambers JC.
We determined the relationships between DNA sequence variation and DNA methylation using blood samples from 3,799 Europeans and 3,195 South Asians. We identify 11,165,559 SNP–CpG associations (methylation quantitative trait loci (meQTL), P < 10−14), including 467,915 meQTL that operate in trans. The meQTL are enriched for functionally relevant characteristics, including shared chromatin state, High-throuhgput chromosome conformation interaction, and association with gene expression, metabolic variation and clinical traits. We use molecular interaction and colocalization analyses to identify multiple nuclear regulatory pathways linking meQTL loci to phenotypic variation, including UBASH3B (body mass index), NFKBIE (rheumatoid arthritis), MGA (blood pressure) and COMMD7 (white cell counts). For rs6511961, chromatin immunoprecipitation followed by sequencing (ChIP–seq) validates zinc finger protein (ZNF)333 as the likely trans acting effector protein. Finally, we used interaction analyses to identify population- and lineage-specific meQTL, including rs174548 in FADS1, with the strongest effect in CD8+ T cells, thus linking fatty acid metabolism with immune dysregulation and asthma. Our study advances understanding of the potential pathways linking genetic variation to human phenotype.