Maarten van Iterson

Disease variants alter transcription factor levels and methylation of their binding sites

Most disease-associated genetic variants are noncoding, making it challenging to design experiments to understand their functional consequences. Identification of expression quantitative trait loci (eQTLs) has been a powerful approach to infer the downstream effects of disease-associated variants, but most of these variants remain unexplained. The analysis of DNA methylation, a key component of the epigenome, offers highly complementary data on the regulatory potential of genomic regions. Here we show that disease-associated variants have widespread effects on DNA methylation in trans that likely reflect differential occupancy of trans binding sites by cis-regulated transcription factors. Using multiple omics data sets from 3,841 Dutch individuals, we identified 1,907 established trait-associated SNPs that affect the methylation levels of 10,141 different CpG sites in trans (false discovery rate (FDR) < 0.05). These included SNPs that affect both the expression of a nearby transcription factor (such as NFKB1, CTCF and NKX2-3) and methylation of its respective binding site across the genome. Trans methylation QTLs effectively expose the downstream effects of disease-associated variants.

Genetic and environmental influences interact with age and sex in shaping the human methylome

The methylome is subject to genetic and environmental effects. Their impact may depend on sex and age, resulting in sex- and age-related physiological variation and disease susceptibility. Here we estimate the total heritability of DNA methylation levels in whole blood and estimate the variance explained by common single nucleotide polymorphisms at 411,169 sites in 2,603 individuals from twin families, to establish a catalogue of between-individual variation in DNA methylation. Heritability estimates vary across the genome (mean=19%) and interaction analyses reveal thousands of sites with sex-specific heritability as well as sites where the environmental variance increases with age. Integration with previously published data illustrates the impact of genome and environment across the lifespan at methylation sites associated with metabolic traits, smoking and ageing. These findings demonstrate that our catalogue holds valuable information on locations in the genome where methylation variation between people may reflect disease-relevant environmental exposures or genetic variation.

Identification of context-dependent expression quantitative trait loci in whole blood

Genetic risk factors often localize to noncoding regions of the genome with unknown effects on disease etiology. Expression quantitative trait loci (eQTLs) help to explain the regulatory mechanisms underlying these genetic associations. Knowledge of the context that determines the nature and strength of eQTLs may help identify cell types relevant to pathophysiology and the regulatory networks underlying disease. Here we generated peripheral blood RNA–seq data from 2,116 unrelated individuals and systematically identified context-dependent eQTLs using a hypothesis-free strategy that does not require previous knowledge of the identity of the modifiers. Of the 23,060 significant cis-regulated genes (false discovery rate (FDR) ≤ 0.05), 2,743 (12%) showed context-dependent eQTL effects. The majority of these effects were influenced by cell type composition. A set of 145 cis-eQTLs depended on type I interferon signaling. Others were modulated by specific transcription factors binding to the eQTL SNPs.

The Effects of Low Levels of Dystrophin on Mouse Muscle Function and Pathology

Duchenne muscular dystrophy (DMD) is a severe progressive muscular disorder caused by reading frame disrupting mutations in the DMD gene, preventing the synthesis of functional dystrophin. As dystrophin provides muscle fiber stability during contractions, dystrophin negative fibers are prone to exercise-induced damage. Upon exhaustion of the regenerative capacity, fibers will be replaced by fibrotic and fat tissue resulting in a progressive loss of function eventually leading to death in the early thirties. With several promising approaches for the treatment of DMD aiming at dystrophin restoration in clinical trials, there is an increasing need to determine more precisely which dystrophin levels are sufficient to restore muscle fiber integrity, protect against muscle damage and improve muscle function. To address this we generated a new mouse model (mdx-Xist Δhs) with varying, low dystrophin levels (3–47%, mean 22.7%, stdev 12.1, n = 24) due to skewed X-inactivation. Longitudinal sections revealed that within individual fibers, some nuclei did and some did not express dystrophin, resulting in a random, mosaic pattern of dystrophin expression within fibers. Mdx-Xist Δhs, mdx and wild type females underwent a 12 week functional test regime consisting of different tests to assess muscle function at base line, or after chronic treadmill running exercise. Overall, mdx-Xist Δhs mice with 3–14% dystrophin outperformed mdx mice in the functional tests. Improved histopathology was observed in mice with 15–29% dystrophin and these levels also resulted in normalized expression of pro-inflammatory biomarker genes, while for other parameters >30% of dystrophin was needed. Chronic exercise clearly worsened pathology, which needed dystrophin levels >20% for protection. Based on these findings, we conclude that while even dystrophin levels below 15% can improve pathology and performance, levels of >20% are needed to fully protect muscle fibers from exercise-induced damage.

Blood lipids influence DNA methylation in circulating cells

BackgroundCells 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.ResultsThis 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 ABCG1 (for HDL cholesterol), which are all key regulators of lipid metabolism.ConclusionsOur analysis suggests a role for epigenetic priming in end-product feedback control of lipid metabolism and highlights Mendelian randomization as an effective tool to infer causal relationships in integrative genomics data.

MethylAid: visual and interactive quality control of large Illumina 450k datasets

UNLABELLED The Illumina 450k array is a frequently used platform for large-scale genome-wide DNA methylation studies, i.e. epigenome-wide association studies. Currently, quality control of 450k data can be performed with Illuminas GenomeStudio and is part of a limited number 450k analysis pipelines. However, GenomeStudio cannot handle large-scale studies, and existing pipelines provide limited options for quality control and neither support interactive exploration by the user. To aid the detection of bad-quality samples in large-scale genome-wide DNA methylation studies as flexible and transparent as possible, we have developed MethylAid; a visual and interactive Web application using RStudios shiny package. Bad-quality samples are detected using sample-dependent and sample-independent quality control probes present on the array and user-adjustable thresholds. In-depth exploration of bad-quality samples can be performed using several interactive diagnostic plots. Furthermore, plots can be annotated with user-provided metadata, for example, to identify outlying batches. Our new tool makes quality assessment of 450k array data interactive, flexible and efficient and is, therefore, expected to be useful for both data analysts and core facilities. AVAILABILITY AND IMPLEMENTATION MethylAid is implemented as an R/Bioconductor package (www.bioconductor.org/packages/3.0/bioc/html/MethylAid.html). A demo application is available from shiny.bioexp.nl/MethylAid.

Single Molecule Sequencing of Free DNA from Maternal Plasma for Noninvasive Trisomy 21 Detection

BACKGROUND Noninvasive fetal aneuploidy detection by use of free DNA from maternal plasma has recently been shown to be achievable by whole genome shotgun sequencing. The high-throughput next-generation sequencing platforms previously tested use a PCR step during sample preparation, which results in amplification bias in GC-rich areas of the human genome. To eliminate this bias, and thereby experimental noise, we have used single molecule sequencing as an alternative method. METHODS For noninvasive trisomy 21 detection, we performed single molecule sequencing on the Helicos platform using free DNA isolated from maternal plasma from 9 weeks of gestation onwards. Relative sequence tag density ratios were calculated and results were directly compared to the previously described Illumina GAII platform. RESULTS Sequence data generated without an amplification step show no GC bias. Therefore, with the use of single molecule sequencing all trisomy 21 fetuses could be distinguished more clearly from euploid fetuses. CONCLUSIONS This study shows for the first time that single molecule sequencing is an attractive and easy to use alternative for reliable noninvasive fetal aneuploidy detection in diagnostics. With this approach, previously described experimental noise associated with PCR amplification, such as GC bias, can be overcome.

Comparison of skeletal muscle pathology and motor function of dystrophin and utrophin deficient mouse strains

The genetic defect of mdx mice resembles that of Duchenne muscular dystrophy, although their functional performance and life expectancy is nearly normal. By contrast, mice lacking utrophin and dystrophin (mdx/utrn -/-) are severely affected and die prematurely. Mice with one utrophin allele (mdx/utrn +/-) are more severely affected than mdx mice, but outlive mdx/utrn -/- mice. We subjected mdx/utrn +/+, +/-, -/- and wild type males to a 12week functional test regime of four different functional tests. Mdx/utrn +/+ and +/- mice completed the regime, while mdx/utrn -/- mice died prematurely. Mdx/utrn +/- mice performed significantly worse compared to mdx/utrn +/+ mice in functional tests. Creatine kinase levels, percentage of fibrotic/necrotic tissue, morphology of neuromuscular synapses and expression of biomarker genes were comparable, whereas mdx/utrn +/- and -/- mice had increased levels of regenerating fibers. This makes mdx/utrn +/- mice valuable for testing the benefit of potential therapies on muscle function parameters.

Age-related accrual of methylomic variability is linked to fundamental ageing mechanisms

BackgroundEpigenetic change is a hallmark of ageing but its link to ageing mechanisms in humans remains poorly understood. While DNA methylation at many CpG sites closely tracks chronological age, DNA methylation changes relevant to biological age are expected to gradually dissociate from chronological age, mirroring the increased heterogeneity in health status at older ages.ResultsHere, we report on the large-scale identification of 6366 age-related variably methylated positions (aVMPs) identified in 3295 whole blood DNA methylation profiles, 2044 of which have a matching RNA-seq gene expression profile. aVMPs are enriched at polycomb repressed regions and, accordingly, methylation at those positions is associated with the expression of genes encoding components of polycomb repressive complex 2 (PRC2) in trans. Further analysis revealed trans-associations for 1816 aVMPs with an additional 854 genes. These trans-associated aVMPs are characterized by either an age-related gain of methylation at CpG islands marked by PRC2 or a loss of methylation at enhancers. This distinct pattern extends to other tissues and multiple cancer types. Finally, genes associated with aVMPs in trans whose expression is variably upregulated with age (733 genes) play a key role in DNA repair and apoptosis, whereas downregulated aVMP-associated genes (121 genes) are mapped to defined pathways in cellular metabolism.ConclusionsOur results link age-related changes in DNA methylation to fundamental mechanisms that are thought to drive human ageing.

Controlling bias and inflation in epigenome- and transcriptome-wide association studies using the empirical null distribution

We show that epigenome- and transcriptome-wide association studies (EWAS and TWAS) are prone to significant inflation and bias of test statistics, an unrecognized phenomenon introducing spurious findings if left unaddressed. Neither GWAS-based methodology nor state-of-the-art confounder adjustment methods completely remove bias and inflation. We propose a Bayesian method to control bias and inflation in EWAS and TWAS based on estimation of the empirical null distribution. Using simulations and real data, we demonstrate that our method maximizes power while properly controlling the false positive rate. We illustrate the utility of our method in large-scale EWAS and TWAS meta-analyses of age and smoking.