Genome-wide identification of genes regulating DNA methylation using genetic anchors for causal inference

BACKGROUND: DNA methylation is a key epigenetic modification in human development and disease, yet there is limited understanding of its highly coordinated regulation. Here, we identify 818 genes that affect DNA methylation patterns in blood using large-scale population genomics data. RESULTS: By employing genetic instruments as causal anchors, we establish directed associations between gene expression and distant DNA methylation levels, while ensuring specificity of the associations by correcting for linkage disequilibrium and pleiotropy among neighboring genes. The identified genes are enriched for transcription factors, of which many consistently increased or decreased DNA methylation levels at multiple CpG sites. In addition, we show that a substantial number of transcription factors affected DNA methylation at their experimentally determined binding sites. We also observe genes encoding proteins with heterogenous functions that have widespread effects on DNA methylation, e.g., NFKBIE, CDCA7(L), and NLRC5, and for several examples, we suggest plausible mechanisms underlying their effect on DNA methylation. CONCLUSION: We report hundreds of genes that affect DNA methylation and provide key insights in the principles underlying epigenetic regulation

Blood lipids influence DNA methylation in circulating cells

Background: Cells can be primed by external stimuli to obtain a long-term epigenetic memory. We hypothesize that long-term exposure to elevated blood lipids can prime circulating immune cells through changes in DNA methylation, a process that may contribute to the development of atherosclerosis. To interrogate the causal relationship between triglyceride, low-density lipoprotein (LDL) cholesterol, and high-density lipoprotein (HDL) cholesterol levels and genome-wide DNA methylation while excluding confounding and pleiotropy, we perform a stepwise Mendelian randomization analysis in whole blood of 3296 individuals. Results: This analysis shows that differential methylation is the consequence of inter-individual variation in blood lipid levels and not vice versa. Specifically, we observe an effect of triglycerides on DNA methylation at three CpGs, of LDL cholesterol at one CpG, and of HDL cholesterol at two CpGs using multivariable Mendelian randomization. Using RNA-seq data available for a large subset of individuals (N = 2044), DNA methylation of these six CpGs is associated with the expression of CPT1A and SREBF1 (for triglycerides), DHCR24 (for LDL cholesterol) and

Genomic and phenotypic insights from an atlas of genetic effects on DNA methylation

Characterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. In the present study, we describe results of DNAm quantitative trait locus (mQTL) analyses on 32,851 participants, identifying genetic variants associated with DNAm at 420,509 DNAm sites in blood. We present a database of >270,000 independent mQTLs, of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15–17% of the additive genetic variance of DNAm. We show that the genetic architecture of DNAm levels is highly polygenic. Using shared genetic control between distal DNAm sites, we constructed networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic variants are associated with both DNAm levels and complex diseases, but only in a minority of cases do these associations reflect causal relationships from DNAm to trait or vice versa, indicating a more complex genotype–phenotype map than previously anticipated

Genomic and phenotypic insights from an atlas of genetic effects on DNA methylation.

Characterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. In the present study, we describe results of DNAm quantitative trait locus (mQTL) analyses on 32,851 participants, identifying genetic variants associated with DNAm at 420,509 DNAm sites in blood. We present a database of >270,000 independent mQTLs, of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15–17% of the additive genetic variance of DNAm. We show that the genetic architecture of DNAm levels is highly polygenic. Using shared genetic control between distal DNAm sites, we constructed networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic variants are associated with both DNAm levels and complex diseases, but only in a minority of cases do these associations reflect causal relationships from DNAm to trait or vice versa, indicating a more complex genotype–phenotype map than previously anticipated.C.L.R., G.D.S., G.S., J.L.M., K.B., M. Suderman, T.G.R. and T.R.G. are supported by the UK Medical Research Council (MRC) Integrative Epidemiology Unit at the University of Bristol (MC_UU_00011/1, MC_UU_00011/4, MC_UU_00011/5). C.L.R. receives support from a Cancer Research UK Programme grant (no. C18281/A191169). G.H. is funded by the Wellcome Trust and the Royal Society (208806/Z/17/Z). E.H. and J.M. were supported by MRC project grants (nos. MR/K013807/1 and MR/R005176/1 to J.M.) and an MRC Clinical Infrastructure award (no. MR/M008924/1 to J.M.). B.T.H. is supported by the Netherlands CardioVascular Research Initiative (the Dutch Heart Foundation, Dutch Federation of University Medical Centres, the Netherlands Organisation for Health Research and Development, and the Royal Netherlands Academy of Sciences) for the GENIUS project ‘Generating the best evidence-based pharmaceutical targets for atherosclerosis’ (CVON2011-19, CVON2017-20). J.T.B. was supported by the Economic and Social Research Council (grant no. ES/N000404/1). The present study was also supported by JPI HDHL-funded DIMENSION project (administered by the BBSRC UK, grant no. BB/S020845/1 to J.T.B., and by ZonMW the Netherlands, grant no. 529051021 to B.T.H). A.D.B. has been supported by a Wellcome Trust PhD Training Fellowship for Clinicians and the Edinburgh Clinical Academic Track programme (204979/Z/16/Z). J. Klughammer was supported by a DOC fellowship of the Austrian Academy of Sciences. Cohort-specific acknowledgements and funding are presented in the Supplementary Note

Genome-wide identification of directed gene networks using large-scale population genomics data

Identification of causal drivers behind regulatory gene networks is crucial in understanding gene function. Here, we develop a method for the large-scale inference of gene–gene interactions in observational population genomics data that are both directed (using local genetic instruments as causal anchors, akin to Mendelian Randomization) and specific (by controlling for linkage disequilibrium and pleiotropy). Analysis of genotype and whole-blood RNA-sequencing data from 3072 individuals identified 49 genes as drivers of downstream transcriptional changes (Wald P < 7 × 10−10), among which transcription factors were overrepresented (Fisher’s P = 3.3 × 10−7). Our analysis suggests new gene functions and targets, including for SENP7 (zinc-finger genes involved in retroviral repression) and BCL2A1 (target genes possibly involved in auditory dysfunction). Our work highlights the utility of population genomics data in deriving directed gene expression networks. A resource of trans-effects for all 6600 genes with a genetic instrument can be explored individually using a web-based browser

Upper Extremity Muscle Strength in Children With Unilateral Spastic Cerebral Palsy:A Bilateral Problem?

Objective. The objective was to investigate whether muscle strength in the nonaffected and affected upper extremities (UEs) in children (7-12 years) with unilateral spastic cerebral palsy (USCP) differs from that in children with typical development (TD). Methods. A cross-sectional study design was used. Isometric arm strength (wrist flexion, wrist extension with flexed and extended fingers, elbow flexion/extension) was assessed in 72 children (mean age = 9.3 [SD = 1.9] years) with USCP, and isometric grip/pinch strength was assessed in 86 children (mean age = 9.3 [SD = 1.8] years) with USCP. Arm/grip/pinch strength was assessed in 120 children (mean age = 9.5 [SD = 1.7] years) with TD. Arm strength was measured with a hand-held dynamometer, and grip/pinch strength was measured with a calibrated, modified (digitized) grip dynamometer and a pinch meter. The nonaffected UE of children with USCP was compared with the preferred UE of children with TD because both sides represent the preferred UE. The affected UE was compared with the nonpreferred UE of children with TD, as both sides represent the nonpreferred UE. Results. In all measurements except for grip strength of the preferred UE, children with USCP were weaker than children with TD. Conclusions. In children with USCP, muscle strength weakness exists in both UEs. Impact. When unimanual or bimanual ability limitations are present in children with unilateral cerebral palsy, investigation of the muscle strength of the nonaffected UE should be part of the assessment

Upper Extremity Muscle Strength in Children With Unilateral Spastic Cerebral Palsy: A Bilateral Problem?

Objective. The objective was to investigate whether muscle strength in the nonaffected and affected upper extremities (UEs) in children (7-12 years) with unilateral spastic cerebral palsy (USCP) differs from that in children with typical development (TD). Methods. A cross-sectional study design was used. Isometric arm strength (wrist flexion, wrist extension with flexed and extended fingers, elbow flexion/extension) was assessed in 72 children (mean age = 9.3 [SD = 1.9] years) with USCP, and isometric grip/pinch strength was assessed in 86 children (mean age = 9.3 [SD = 1.8] years) with USCP. Arm/grip/pinch strength was assessed in 120 children (mean age = 9.5 [SD = 1.7] years) with TD. Arm strength was measured with a hand-held dynamometer, and grip/pinch strength was measured with a calibrated, modified (digitized) grip dynamometer and a pinch meter. The nonaffected UE of children with USCP was compared with the preferred UE of children with TD because both sides represent the preferred UE. The affected UE was compared with the nonpreferred UE of children with TD, as both sides represent the nonpreferred UE. Results. In all measurements except for grip strength of the preferred UE, children with USCP were weaker than children with TD. Conclusions. In children with USCP, muscle strength weakness exists in both UEs. Impact. When unimanual or bimanual ability limitations are present in children with unilateral cerebral palsy, investigation of the muscle strength of the nonaffected UE should be part of the assessment

Lipid-induced transcriptomic changes in blood link to lipid metabolism and allergic response

Immune cell function can be altered by lipids in circulation, a process potentially relevant to lipid-associated inflammatory diseases including atherosclerosis and rheumatoid arthritis. To gain further insight in the molecular changes involved, we here perform a transcriptome-wide association analysis of blood triglycerides, HDL cholesterol, and LDL cholesterol in 3229 individuals, followed by a systematic bidirectional Mendelian randomization analysis to assess the direction of effects and control for pleiotropy. Triglycerides are found to induce transcriptional changes in 55 genes and HDL cholesterol in 5 genes. The function and cell-specific expression pattern of these genes implies that triglycerides downregulate both cellular lipid metabolism and, unexpectedly, allergic response. Indeed, a Mendelian randomization approach based on GWAS summary statistics indicates that several of these genes, including interleukin-4 (IL4) and IgE receptors (FCER1A, MS4A2), affect the incidence of allergic diseases. Our findings highlight the interplay between triglycerides and immune cells in allergic disease