An Epigenetic Signature in Peripheral Blood Associated with the Haplotype on 17q21.31, a Risk Factor for Neurodegenerative Tauopathy
Progressive supranuclear palsy (PSP) and frontotemporal dementia (FTD) are two neurodegenerative diseases linked, at the pathologic and genetic level, to the microtubule associated protein tau. We studied epigenetic changes (DNA methylation levels) in peripheral blood from patients with PSP, FTD, and unaffected controls. Analysis of genome-wide methylation patterns revealed significant differentially methylated probes in patients versus unaffected controls. Remarkably, differentially methylated probes in PSP vs. controls were preferentially clustered within the 17q21.31 region, previously known to harbor the major genetic risk factor for PSP. We identified and replicated a dose-dependent effect of the risk-associated H1 haplotype on methylation levels within the region in independent datasets in blood and brain. These data reveal that the H1 haplotype increases risk for tauopathy via differential methylation, indicating a mediating role for methylation in dementia pathophysiology.
Published in the journal:
. PLoS Genet 10(3): e32767. doi:10.1371/journal.pgen.1004211
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pgen.1004211
Summary
Progressive supranuclear palsy (PSP) and frontotemporal dementia (FTD) are two neurodegenerative diseases linked, at the pathologic and genetic level, to the microtubule associated protein tau. We studied epigenetic changes (DNA methylation levels) in peripheral blood from patients with PSP, FTD, and unaffected controls. Analysis of genome-wide methylation patterns revealed significant differentially methylated probes in patients versus unaffected controls. Remarkably, differentially methylated probes in PSP vs. controls were preferentially clustered within the 17q21.31 region, previously known to harbor the major genetic risk factor for PSP. We identified and replicated a dose-dependent effect of the risk-associated H1 haplotype on methylation levels within the region in independent datasets in blood and brain. These data reveal that the H1 haplotype increases risk for tauopathy via differential methylation, indicating a mediating role for methylation in dementia pathophysiology.
Introduction
Epigenetics is one of the most rapidly expanding fields in biology, and is uncovering additional levels of complexity in the human genome, including DNA methylation, histone modifications, and intra- and inter-chromosomal interactions mediated by chromatin proteins [1], [2]. Changes in methylation represent a key area where environmental factors can modify or interact with inherited genetic factors (DNA sequence) to alter the functional output of the genome. Disease-causing genes involved in epigenetic modifications have been identified, most notably for neurodevelopmental disorders such as Rett syndrome [3]. A very limited number of studies have addressed specific epigenetic modifications relevant to neurological diseases and dementia (reviewed in [4]–[7]). Additionally, epigenetic signatures have been reported for different brain regions [8], , for regional brain aging [10], and aging in general [11] further supporting epigenetic studies in patients with neurodegenerative diseases.
Progressive supranuclear palsy (PSP) is a neurodegenerative disease typically characterized by parkinsonism, postural instability, and cognitive impairment [12]. Pathologically, PSP is defined by the accumulation of tau protein in subcortical and cortical regions, (reviewed in Williams and Lees, 2009 [13]), showing substantial overlap with other neurodegenerative diseases characterized by tau accumulation and grouped under the generic name of tauopathies, including approximately one-half of all frontotemporal dementia (FTD) cases and Alzheimer's disease [14]. Both rare [15], [16] and common [17] genetic variation have been shown to mediate risk for tauopathies. The major common variant risk for PSP, a prototypical tauopathy, involves a region surrounding the tau locus [18], but how such genetic variation might mediate risk is not known.
We profiled the methylation status in peripheral blood from patients with two tau-related neurodegenerative conditions, PSP and FTD, using Illumina DNA methylation arrays. We then integrated these methylation data with SNP and gene expression data to identify a mediating role for methylation in genetic risk for PSP. We replicate this finding in independent studies and show that it is conserved in brain, providing the first evidence for a role for DNA methylation in mediating the risk for neurodegenerative dementia.
Results
Differential methylation analysis
We first analyzed methylation profiles in 171 patients with FTD (n = 128) and PSP (n = 43) and compared them with 185 subjects with no evidence of dementia or other neurological conditions using Illumina HumanMethylation 450 k arrays (Table S1). Two datasets were generated in two batches, samples were compared within each dataset to condition out a potential batch effect, and the resulting differentially methylated probes (DMPs) were combined (see Methods).
Differential methylation analysis identified a number of DMPs between affected subjects and controls, with partial overlap between PSP and FTD (Figure 1a–b, complete list of DMPs is in Table S2). DMPs were mostly clustered within CpG islands (defined according to the Illumina annotation), with most being hypermethylated in PSP vs. controls (Figure 1c, Table 1). Gene ontology analysis of DMPs in PSP vs. controls showed overrepresentation of genes involved in a number of pathways, including DNA binding and transcription factor binding (Figure S1). We then assessed the chromosomal distribution of the DMPs, and observed – only in PSP samples vs. controls – an overrepresentation of probes from chromosomes 19 (hypergeometric test p-value = 1.32×10−6), 22 (p = 8.63×10−6), and 17 (p = 5.82×10−5, Figure 1d), with most top DMPs (after filtering for absolute average beta difference (aβD)>0.1) located within the 17q21.31 region (Figure 1e). The most significant DMPs when comparing PSP vs. controls (n = 14, absolute aβD>0.1) are listed in Table 2. Of note, 4 DMPs (all hypomethylated in PSP) were located within the NFYA gene, encoding for a component of a nuclear transcription factor. Importantly, 3 of the 14 significant DMPs are located in 17q21.31 (Figure 1e, p = 2.23×10−7, hypergeometric test). Despite being located in a relatively limited genomic region, these 3 probes were both hypermethylated and hypomethylated in PSP vs. controls, suggesting complex disease-associated patterns of differential methylation in this region.
17q21.31 haplotype and methylation
The location of several DMPs in the 17q21.31 region was intriguing because the 17q21.31 locus contains an established risk factor for neurodegeneration, first reported in 1997 by Conrad et al. [18] for PSP and then confirmed in multiple series (reviewed in Wade-Martins, 2012 [19]). Two main haplotypes (H1 and H2) have been described at this locus. The more common H1 haplotype is over-represented (95% vs. 57%) in PSP vs. normal controls [18]–[20]. The H1/H2 locus spans at least 1.8 Mb and includes multiple genes (>40, many of which are actively transcribed in the brain), notably including MAPT, encoding for the microtubule-associated protein tau [19]. Mutations in MAPT cause FTD and PSP, and hyperphosphorylated tau accumulation is a hallmark in a number of neurodegenerative conditions, including AD, PSP, FTD and others, collectively named ‘tauopathies’.
Consistent with previous reports, the H1 haplotype was overrepresented in our PSP cohort, with a H1 allelic frequency of 97.1% vs. 80.4% in controls (p = 1.86×10−4, Fisher's exact test, Table S1), further confirming – even in this relatively small data set – the H1 haplotype as a risk factor for PSP. We hypothesized that the clustering of DMPs in 17q21.31 in PSP cases vs. controls might be related to the H1 haplotype risk factor. To detect an effect of the 17q21.31 haplotype on methylation levels, we compared samples based on their genotype at this region, independent of disease classification. As for previous analyses, we compared samples within datasets to avoid potential batch effects. Genotype distribution across diseases and datasets is reported in Table S3.
We compared carriers of the risk-associated H1 haplotype (H1/H1 and H1/H2 genotypes) to H2/H2 samples (dominant model) within each dataset and, after filtering DMPs for adjusted p≤0.05, identified two overlapping sets of 57 and 34 DMPs (Figure 2a,b), markedly clustered within the 17q21.31 region (Figure 2d). Similar results (Figure 2a,c) were obtained when comparing H1/H1 samples to H2 carriers (H1/H2 and H2/H2, recessive model), supporting the hypothesis of a strong cis effect of the H1/H2 locus on methylation levels in peripheral blood (Figure 2). After filtering for absolute aβD>0.1, 8 of the top 9 DMPs identified in both datasets were within 17q21.31 (Table 3) in the dominant model. As noted in PSP cases vs. controls, DMPs in this region are both hyper- and hypo-methylated, suggesting a complex cis-regulation of methylation levels (Figure 3a). Scatterplots of the methylation levels for the top DMPs shared between the dominant and recessive models indicate that the H1 haplotype influences methylation levels at these sites in a dose-dependent fashion (Figure 3b), accounting for a majority of methylation variability at these sites (e.g. R-squared = 0.835 and 0.866 in dataset #1 and #2, respectively, for cg22968622). Similar results were obtained when comparing subjects based on their genotype at 17q21.31, but only within controls, FTD, or AD patients (Text S1). The H1 haplotype can be further divided into sub-haplotypes [21]. We obtained sub-haplotype information for 93 H1 carriers in our cohort using the SNPs described in Kauwe et al 2008 [22]. Hierarchical clustering of the methylation signal in the 17q21.31 region and principal component analysis did not reveal a particular clustering of H1 sub-haplotypes (data not shown). These results – although based on a subset of our cohort – suggest that haplotype structure is not the major determinant of 17q21.31 methylation overall.
To test the contribution of haplotype status on PSP-associated DMPs, we repeated the differential methylation analysis only on samples with the H1/H1 haplotype (n = 31 PSP cases and 59 unaffected controls). Of the resulting 341 significant DMPs (after application of the Benjamini-Hochberg procedure, FDR = 0.05), 21 were located in chromosome 17 and 2 were located in the 17q21.31 band. Neither the chromosome nor the region were found to be significantly overrepresented by the hypergeometric test (p = 0.342 and 0.149, respectively). The lack of overrepresentation within 17q21.31 after conditioning on strata defined by the 17q21.31 haplotype suggests that either the previously identified overrepresentation on chromosome 17 was due to the 17q21.31 haplotype effect on methylation levels, or that it could also reflect reduced power due to small sample size. To address the issue that the strata contained too few samples, we also carried out a multivariate regression model analysis that included 17q21.31 haplotype as covariate. Specifically, the methylation level of each of the 3 top PSP-related DMPs located in 17q21.31 (Table 2) was regressed on PSP status, 17q21.31 haplotype, ethnicity, and age using a multivariate linear regression model. We found that, except for cg23758822, other PSP-related DMPs were no longer significant (p = 0.410 on average, Table S4) in a multivariate model once it included the H1 genotype. We also calculated the relative weight of each predictor using the R package relaimpo [23], and determined that the H1 haplotype accounted for the majority of explained variance (78.2±25.9%, Figure S2). Finally, we estimated relative cell count composition in peripheral blood using methylation data [24]–[26]. Correction for inferred cell count did not significantly change our findings (Text S1, Table S5, Figure S11).
Taken together, these findings indicate 1) a strong effect of the 17q21.31 haplotype on methylation levels at 17q21.31, 2) that the risk-associated H1 determines most of the methylation changes observed with confidence in PSP patients vs. controls, and 3) that additional DMPs outside the 17q21.31 region may be at play in determining risk susceptibility for PSP in H1 carriers, though larger sample sizes will be needed to clarify their importance.
Genome-wide methylation QTL analysis confirms a cis methQTL at 17q21.31
Our findings strongly indicate a cis regulation of methylation levels at the 17q21.31 locus. To test whether there were additional potential genetic determinants of methylation levels at 17q21.31 in our dataset, we performed a methylation QTL (methQTL) analysis in a subset of 226 individuals of European descent for whom whole-genome SNP and methylation data were available (Table S6). We assessed association of genetic variants with methylation levels at 3 CpGs within 17q21.31 (cg22968622, cg17117718, cg19832721) in each dataset. We identified on average 110 genome-wide significant signals (Bonferroni-adjusted p≤0.05), all located within the 17q21.31 region (Figure 3c, Table S7). These variants accounted for a proportion of variability ranging between 25.5% and 98.2% (mean R-squared = 0.701, Figure S3, Table S8) further confirming that genetic variants at 17q21.31 are controlling methylation levels in cis in the same region. We focused on Caucasian individuals because of the differences in frequency of the H2 haplotype across populations. In fact, consistent with previous reports [27], we observed that the H2 haplotype occurs more frequently in Caucasians (H2 allelic frequency = 19.2%, Table 4) than in other ethnic groups (H2 allelic frequency in Asians = 1.3%; p = 3.83×10−6, Fisher's exact test). However, similar results were observed when including all the 273 individuals for whom SNP and methylation data were available (Text S1, Figure S4, Tables S9, S10).
Cis methQTL effects at the 17q21.31 locus in additional datasets
To confirm that the 17q21.31 haplotype regulates methylation in cis at this locus in an independent dataset, we downloaded and reanalyzed raw data from a previously published study, for which SNP and methylation data in peripheral blood from 12 samples were publicly available [28]. Using the rs1052553 SNP to call the H1/H2 haplotype and adopting the same statistical thresholds, we compared H1/H1 vs. H1/H2 subjects and identified one hypomethylated probe (cg22968622, adjusted p-value = 2.37×10−8, aβD = −0.42) within 17q21.31, which was also identified in our analysis. We also performed a methQTL analysis for cg22968622 in the same dataset. Of the 310 significant SNPs, 206 were located in the 17q21.31 region (p = 0, hypergeometric test), further supporting the presence of a cis methQTL at this locus.
To provide independent validation of the methylation array assay, we performed reduced representation bisulfite sequencing (RRBS) on a representative set of 7 samples from the study (2 H1/H1 controls, 1 H1/H1 PSP patient, 1 H1/H2 control, 1 H1/H2 PSP patient, and 2 H2/H2 controls). As a sequencing-based approach, RRBS would not suffer from some of the technical biases present in arrays, e.g. due to hybridization. At CpG sites that were covered by both RRBS and array, the methylation measurements were highly correlated (Pearson r>0.9) in all seven samples (Figure S5).
To validate our findings from peripheral blood, we analyzed RRBS data from whole-blood DNA of a separate cohort of 80 healthy subjects (comprising 54 H1/H1, 24 H1/H2, and 2 H2/H2). On average, the methylation level computed from RRBS was highly correlated with the array in both dataset #1 (r = 0.965) and dataset #2 (r = 0.963) (Figure S6). Consistent with the array, we found differences in methylation that were significant even after strict Bonferroni correction for multiple testing, mostly localized to the 17q21.31 cytoband (Table S11). Since local methylation levels are often highly correlated, we would expect the differentially methylated CpGs identified by both the array and the sequencing method to be in close proximity. Indeed, the differentially methylated loci identified by RRBS in the 17q21.31 region are nearby those identified by the Illumina Human Methylation array, overlapping the same genes MAPT and KIAA1267. This degree of overlap is striking, given the large extent of the haplotype inversion (Figure S7). Methylated regions identified by the array but not by RRBS may be a result of the higher power (greater sample numbers) in the array, and differences in coverage. Therefore, we looked at probes that both demonstrated haplotype-specific methylation in 17q21.31 on the Illumina array and were covered by RRBS reads in the additional cohort. One probe, cg08113562, met these criteria. The methylation pattern followed a statistically significant dose-dependent relationship with the H1 versus H2 haplotype, with mean methylation fractions of 0.001 in H1/H1 subjects, 0.022 in H1/H2 subjects, and 0.048 in H2/H2 subjects (two-sided p = 0.03, ANOVA), in the same direction as that reported by the array.
To assess the relevance of our findings to brain tissue, we analyzed 2 published methylation QTL studies in brain involving 150 [29] and 153 [30] subjects, respectively. Gibbs et al. [29] identified 9 SNP-CpG association pairs (out of 52,345 significant methQTL, mean R-squared = 0.232) at the 17q21.31 locus in two (frontal cortex and cerebellum) of the four studied brain regions; Zhang et al. [30] identified in cerebellar samples 122 SNP-CpG pairs (significant in at least one of the three thresholds they used) at the 17q21.31 locus (out of 12,117 significant methQTLs, mean R-squared = 0.136). Together, these data demonstrate that the cis methQTL we identified in our study are present in independent studies in peripheral blood, and are preserved in brain.
Causal inference identifies three methylated regions that may mediate PSP risk
The strong association between the haplotype at 17q21.31 and methylation status raises the question of whether methylation levels mediate the protective or pathogenic effects of haplotype variants. Recent developments in the field of causal inference have yielded quantitative methods to predict the hierarchy of causation given genetic variants [31]–[34]. Network Edge Orienting (NEO), for example, uses structural equation models to choose the best fitting causal model, assuming that the genetic variation is fixed by meiosis and thus “anchors” each model (that is, genotype precedes phenotype) [35]. NEO allows one to evaluate which of five testable causal models (Figure S8) best explains the relationship between genetic variants, methylation levels, and disease status. For instance, the genotype may lead to patterns of methylation that directly contribute to the disease phenotype (Figure S8b). Under this model, the DMPs within 17q21.31 would be the most interesting, as they would correspond to the epigenetic markers mediating the increased risk conferred by the H1 haplotype. Alternatively, the genotype may independently give rise to the methylation and disease phenotype, with neither contributing to the other (Figure S8c). DMPs under this model are only associated with the disease because of the common source of variation due to the 17q21.31 locus.
We applied NEO to calculate a relative fitting index of the “mediation model” for the 9 haplotype-associated DMPs (Table S12) using 35 PSP cases and 184 unaffected controls for whom these data were available. The “mediation model” best explained the methylation pattern in three sites, one of which (cg17117718) was statistically significant (see Methods) (Figure 3d). These results support the hypothesis that methylation status at certain sites likely is a causal mediator of the major known genetic risk related to PSP pathogenesis. Taken together, these results predict – for the first time – a link between epigenetic changes and tauopathies, and will need to be further validated with functional studies.
Haplotype-associated differences in gene expression
Methylation changes have been associated with changes in gene expression. We examined microarray expression data in peripheral blood available for 120 subjects, to test whether the methylation associated with the 17q21.31 haplotype had such an effect. Among 88 healthy subjects with H1H1 haplotype, 24 healthy subjects with the H1H2 haplotype, and 8 healthy subjects with the H2H2 haplotype, we identified three probes significantly differentially expressed in peripheral blood, mapping to MAPK8IP1 (on chromosome 11), LRRC37A4 (located within the 17q21.31 region), and MTFP1 (on chromosome 22, Benjamini-Hochberg adjusted p-values of 2.4×10−20, 8.4×10−5, and 4.0×10−2, respectively). The three probes demonstrated a log fold change of 0.69, −0.25, and 0.10, respectively, in H2 vs. H1 carriers.
The strong haplotype-associated methylation changes identified in our study included DMPs within the MAPT, KIAA1267, ARHGAP27, and DND1 genes. We used linear regression to test whether the 17q21.31 haplotype was associated with differential expression of these genes, and found no correlation between haplotype and gene expression for these transcripts (adjusted R-squared = 0.006, 0.000, −0.008, and −0.011, respectively). Thus, while haplotype was shown to affect mRNA expression of MAPK8IP1 and LRRC37A4 in peripheral blood, there was no detectable correlation between DMP-containing genes and their corresponding expression levels.
Discussion
The goal of this study was to assess whether changes in DNA methylation in peripheral blood are observed in patients with neurodegenerative diseases. By performing microarray-based differential methylation analysis, we identified a methylation signature associated with disease status in PSP and, to a lesser extent, FTD. Using SNP data available in a subset of our series, we showed that a remarkable proportion of the observed changes in methylation status in PSP are associated with a common haplotype at the 17q21.31 locus, strongly suggesting the presence of a cis methylation QTL in this region. Although we included patients with neurodegenerative disorders in our analysis, the observed pattern seems to be related to the haplotype at 17q21.31, independent of disease status. Integrative analyses including SNP and gene expression data support a model whereby genetic variation at the 17q21.31 locus modulates the risk for neurodegenerative tauopathy at least partially via differential methylation.
The H1 haplotype at 17q21.31, a large linkage disequilibrium block due to an inverted chromosomal sequence of ∼970 kb, is the major known risk locus for PSP [36] and other neurodegenerative diseases [17], [37]. Although the genetic contribution of this locus to the risk for neurodegeneration is established and widely replicated, the mechanism by which risk is increased is largely unknown. This region spans from CRHR1 (corticotrophin) to IMP5 (a presenilin homologue) at the centromeric end of LD, while WNT3 and NSF (N-ethylmaleimide-sensitive factor) are at the telomeric end of the LD block [38]; therefore it spans at least 1.8 Mb, including 48 RefSeq genes – many of which actively transcribed in the brain – and constitutes the largest haplotype block in the human genome. Stefansson et al. [39] showed that the complete disequilibrium was due to an inversion occurring in the H2 haplotype relative to the H1 human reference and subsequent absence of recombination between inverted and non-inverted chromosomes. The study of this region to understand susceptibility to neurodegeneration has been mostly focused on one gene, MAPT, encoding for the microtubule-associated protein tau. This focus on tau is well motivated, as hyperphosphorylated tau accumulates within neurofibrillary tangles – the pathological hallmark of AD – and because mutations in MAPT cause FTD, the second most common neurodegenerative dementia. Several in vitro studies have reported alterations of transcription levels in MAPT due to common variants in the region [40], [41], but this finding has not been consistently replicated [42], [43]. More consistent evidence exists for a higher expression of exon 3 in brains from H2 carriers [43], [44] and of exon 10 in H1 carriers [41], [45], suggesting that splicing abnormalities are involved in increasing risk.
Recently, additional genetic evidence has been reported implicating this locus in other neurodegenerative diseases, such as Parkinson's disease [46], essential tremor, and multisystem atrophy [47]. This is important, as these diseases are not typical tauopathies, suggesting that the effect of this risk-associated region may be complex and involve multiple genes [17].
Our results indicate a novel mechanism by which the H1/H2 locus may affect the risk for tauopathies: significant alterations in methylation mediating increased disease susceptibility. Importantly, these methylation changes are not at the MAPT locus only, but are consistently observed in at least 3 neighboring genes as well, suggesting that genes other than MAPT might be at play in increasing disease susceptibility. In addition, DMPs in the region were both hyper and hypo-methylated, suggesting a complex regulation of methylation levels at this locus. Further studies will be needed to understand whether the observed methylation signature at 17q21.31 is increasing susceptibility through a MAPT-dependent or independent mechanism, or both. Interestingly, a recent study focused on rheumatoid arthritis [48] linked a genetic susceptibility region for the disease, the MHC locus, with methylation changes in the same region, supporting the notion that epigenetic changes might mediate complex disease susceptibility induced by genetic risk factors.
This is the first study of DNA methylation levels in blood in PSP and FTD, disorders that mainly affect brain. Although methylation patterns may be tissue specific [49], comparative studies of blood and brain showed both methylation patterns that are tissue-specific and conserved across tissues [50]. We show that the particular H1 haplotype-related methylation pattern identified in blood is at least partially conserved in brain. This is encouraging, since – in contrast to brain – blood is available from living patients, yielding a higher potential for future use as biomarker and the possibility of large-scale studies. We and others have used gene expression in peripheral blood to gain insights into the biology of neurodegenerative disorders [51], [52]. This study supports the notion that a disease-related signature is present in methylation data as well. Finally, we decided to focus on the risk-associated 17q21.31 region as an initial step, but many interesting candidates for further study emerged from the differential methylation analysis in PSP patients vs. controls, namely the nuclear transcription factor NFYA.
Although our analysis of published datasets supports the presence of the 17q21.31-associated methylation signature in brain tissues, further studies focused on brain samples from patients will be needed to test whether methylation changes are contributing to tissue-specific gene expression abnormalities, and ultimately explain the mechanism of action of genetic susceptibility alleles, and the striking regional vulnerability of these disorders.
Materials and Methods
Ethics statement
All subjects and/or their proxies signed informed consents for genetic studies. The research protocol was approved by the University of California San Francisco (UCSF) and Los Angeles (UCLA) University Institutional Review Boards for human research.
Sample description
Patients were enrolled as part of a large genetic study in neurodegenerative dementia (Genetic Investigation in Frontotemporal Dementia, GIFT) at the UCSF Memory and Aging Center (UCSF-MAC) [53]. 371 unrelated subjects were enrolled in the study (Table S1), including patients with neurodegenerative disorders (128 FTD, 43 PSP, and 15 AD), and 185 healthy controls.
Sample preparation
DNA was extracted from peripheral blood using standard methods. No cell sorting or cell selection was conducted, therefore our data measure methylation levels in whole blood. Total RNA was extracted from the same individuals from peripheral blood using Paxgene Blood RNA tubes (Qiagen).
Methylation arrays
Whole-genome methylation patterns were assayed by the Infinium Human Methylation450 BeadChip Kit (96 samples per chip). This work was performed in two stages (each including 2 chips), resulting in a total of 371 samples. Samples were hybridized as follows: Dataset #1: PSP (n = 40), FTD (n = 55), Control (n = 93); Dataset #2: FTD (n = 73) AD (n = 15), Control (n = 92).
Genotyping
Taqman genotyping
Genotype variants at APOE (rs429358 and rs7412) and MAPT H1/H2 (rs1560310) were obtained using Taqman assays. Genome-wide SNP data was obtained using the Illumina HumanOmni1-Quad BeadChip. 17q21.31 sub-haplotypes were obtained by genotyping 6 SNPs as previously reported [22] using Taqman assays.
SNP arrays
High-throughput SNP genotyping data (Illumina HumanOmni1-Quad BeadChip) from a larger dataset containing 702 samples were available for 273 subjects (14 AD, 110 FTD, 15 PSP, and 134 Controls) in this study. SNP genotypes were called and exported from Illumina GenomeStudio (versions 1.6.3 and 1.8.4). Quality control included filtering for 1) SNPs with <95% genotype call rates (n = 157,121), 2) a minor allele frequency <1% (n = 126,502); and 3) with Hardy-Weinberg Equilibrium p-value <1×10−6 in the control group (n = 16,954). A total of 788,694 SNPs were included in the final analysis.
Ethnicity
We inferred ethnicity for 271 samples (out of the 273 for whom SNP data were available), by using SNP clustering compared to Hapmap data. Briefly, MDS analysis was applied on a merged dataset, including 702 samples from our data (273 of which were included in present study) and 1184 subjects from HapMap phase III. MDS plot shows that the first two principal components can cluster samples by ethnicity, and our data had good overlap with HapMap data (Figure S9). We used self-reported ethnicity for an additional 88 samples. For 12 samples ethnicity remained unknown either because we were not able to call ethnicity with certainty using SNP data (n = 2), or because of the lack of SNP data and self-reported ethnicity (n = 10).
Microarray-based gene expression analysis
Microarray expression data (Illumina HumanRef-8 v3.0) were available in 120 subjects with H1/H1 or H1/H2 haplotypes at 17q21.31.
Reduced representation bisulfite sequencing
RRBS was performed on seven samples multiplexed in two lanes of the Illumina HiSeq. Library preparation was performed using the Msp I restriction enzyme as previously described [54]. Read alignment and methylation level calls were performed using BS-Seeker2 [55] (parameters: --aligner = bowtie -m 5 -g hg18.fa -r --low = 50 --up = 500 -a adapter.txt).
Published datasets
SNP data from Heyn et al. [28] reporting the methylome analysis of newborns and centenarians, including 40 samples assayed by Illumina HumanMethylation450 BeadChip and 14 genotyped by Illumina HumanOmni5-Quad BeadChip (both methylation and SNP genotyping data were available for 12 samples), were downloaded from the Gene Expression Omnibus Database (GEO, http://www.ncbi.nlm.nih.gov/geo/, GSE31438). H1/H2 haplotypes were inferred from SNP rs1052553. For methylation data, beta and detection p-values were downloaded from GSE30870, and 65 sites containing missing values were removed.
Statistical analysis
Methylation arrays
In order to avoid potential confounders from batch effects, the two datasets were processed separately. Raw data was processed using the Illumina GenomeStudio software (version 2010.3). Background correction and color normalization were performed using the R package minfi version 1.2.0 [56], and normalization using Subset-quantile Within Array Normalization (SWAN) [57]. Probes were excluded from further analysis if >95% samples had detection p-value >0.01. In summary, 3,027 probes were removed from dataset #1, and 26,306 probes were removed from dataset #2. In order to avoid potential confounders, 66,877 SNP-containing probes were also excluded from further analysis [58]. Beta values (ratio between methylated probe intensity and the overall intensity) were computed using the R package minfi.
Linear models and empirical Bayes methods as implemented in the limma package were used for differential methylation analysis [59]. P-values were adjusted by using the Benjamini-Hochberg (BH) false discovery rate method. Multi-dimensional scaling (MDS) did not show obvious biases between chips within each dataset, but a batch effect could be observed between datasets #1 and #2 (Figure S10), similar to what has been reported in the literature [30], [60]. To avoid potential confounders due to this batch effect, we compared each disease category with the set of controls within the same batch (i.e. conditioning on batch effect). Similar differential methylation results were obtained when we performed a combined analysis across the entire dataset, after correcting for batch effects using ComBat [61] (Text S1). Two filters were applied to conservatively identify differentially methylated sites: 1) a p-value-based filter (BH-adjusted p≤0.05) and 2) an absolute average beta difference (aβD) filter (absolute aβD>0.1). Chromosomal enrichment analysis was performed by using the hypergeometric test as implemented in the R phyper function.
The impact of relative cell counts in peripheral blood was estimated as previously described [1], [2] and based on a subset (n = 385) of the 500 loci whose methylation levels reflect the relative proportions of immune cells in unfractionated whole blood. After estimating the blood cell type distribution for each sample using the methylation level of the 385 loci, we applied a linear mixed-effect model considering (1) main blood cell types distribution as dependent variables; (2) disease status (or 17q21.31 haplotype), age, ethnicity, and gender as fixed effects; and (3) chip number as a random effect (Text S1).
Raw and normalized methylation data were deposited in the Gene Expression Omnibus (GEO, www.ncbi.nlm.nih.gov/geo), accession number: GSE53740.
Methylation QTL analysis (methQTL) was performed by regressing methylation level at select CpG sites of interest on SNP genotypes. Age and the first two principal component generated from MDS analysis were also included in the multivariate regression model as covariates. Linear regression methQTL analysis was performed using PLINK [62].
Microarray-based expression analysis
Raw data were processed as previously described [63]. Briefly, after quantile normalization, batch effects were removed using ComBat [61]. Differential expression analysis was performed using the limma package [59], applying a false discovery rate filter of ≤0.05 and an absolute log fold change filter of >0.1.
Causality analysis
Causality analysis was performed using the software package Network Edge Orienting (NEO) [35], a structural equation modeling software for determining the direction of causality among various phenotypes (e.g. clinical, molecular) given genotype data. Subjects from all batches were included (n = 219). The relative fitting index of the model is estimated by the Single Marker LEO.NB score, defined as the base-10 logarithm of the probability ratio between the mediation model and the next most likely causal model (i.e., a LEO.NB score of 1 means that the fit of the mediation model is 10 times better than that of the next best alternative causal model). The genotype was encoded as the dosage of the minor allele (A) at rs1560310 tagging H2 (i.e., 0 for H1/H1, 1 for H1/H2, and 2 for H2/H2). Gene expression levels were included for differentially expressed probes according to the 17q21.31 haplotype, and were encoded with the ComBat-corrected relative expression levels. Clinical phenotype was encoded as a binary variable (i.e., 0 for unaffected, 1 for affected with PSP). The Single-Marker analysis option was used, and results surpassing the thresholds of a LEO.NB score (a likelihood ratio of model fit) >0.8 and RMSEA index <0.05 were considered significant fits to the mediation model [64].
Supporting Information
Zdroje
1. FeinbergAP (2010) Epigenomics reveals a functional genome anatomy and a new approach to common disease. Nat Biotechnol 28: 1049–1052 doi:10.1038/nbt1010-1049
2. PortelaA, EstellerM (2010) Epigenetic modifications and human disease. Nat Biotechnol 28: 1057–1068 doi:10.1038/nbt.1685
3. ZoghbiHY (2009) Rett syndrome: what do we know for sure? Nat Neurosci 12: 239–240 doi:10.1038/nn0309-239
4. UrdinguioRG, Sanchez-Mut JV, EstellerM (2009) Epigenetic mechanisms in neurological diseases: genes, syndromes, and therapies. Lancet Neurol 8: 1056–1072 doi:10.1016/S1474-4422(09)70262-5
5. JakovcevskiM, AkbarianS (2012) Epigenetic mechanisms in neurological disease. Nat Med 18: 1194–1204 doi:10.1038/nm.2828
6. AkbarianS, BeeriMS, HaroutunianV (2013) Epigenetic Determinants of Healthy and Diseased Brain Aging and Cognition. JAMA Neurol 1–8 doi:10.1001/jamaneurol.2013.1459
7. LuH, LiuX, DengY, QingH (2013) DNA methylation, a hand behind neurodegenerative diseases. Front Aging Neurosci 5: 85 doi:10.3389/fnagi.2013.00085
8. Ladd-AcostaC, PevsnerJ, SabunciyanS, YolkenRH, WebsterMJ, et al. (2007) DNA methylation signatures within the human brain. Am J Hum Genet 81: 1304–1315 doi:10.1086/524110
9. Van EijkKR, de JongS, BoksMP, LangeveldT, ColasF, et al. (2012) Genetic analysis of DNA methylation and gene expression levels in whole blood of healthy human subjects. BMC Genomics 13: 636 doi:10.1186/1471-2164-13-636
10. HernandezDG, NallsMa, GibbsJR, ArepalliS, van der BrugM, et al. (2011) Distinct DNA methylation changes highly correlated with chronological age in the human brain. Hum Mol Genet 20: 1164–1172 doi:10.1093/hmg/ddq561
11. HorvathS (2013) DNA methylation age of human tissues and cell types. Genome Biol 14: R115 doi:10.1186/gb-2013-14-10-r115
12. SteeleJ, RichardsonJ, OlzewskiJ (1964) Progressive supranuclear palsy. Arch Neurol 10: 333–359.
13. WilliamsDR, LeesAJ (2009) Progressive supranuclear palsy: clinicopathological concepts and diagnostic challenges. Lancet Neurol 8: 270–279 doi:10.1016/S1474-4422(09)70042-0
14. BoeveBF (2012) Progressive supranuclear palsy. Parkinsonism Relat Disord 18 Suppl 1: S192–4 doi:10.1016/S1353-8020(11)70060-8
15. HuttonM, LendonCL, RizzuP, BakerM, FroelichS, et al. (1998) Association of missense and 5′-splice-site mutations in tau with the inherited dementia FTDP-17. Nature 393: 702–705.
16. CoppolaG, ChinnathambiS, LeeJJ, DombroskiBa, BakerMC, et al. (2012) Evidence for a role of the rare p.A152T variant in MAPT in increasing the risk for FTD-spectrum and Alzheimer's diseases. Hum Mol Genet 21: 3500–3512 doi:10.1093/hmg/dds161
17. HöglingerGU, MelhemNM, DicksonDW, SleimanPMa, WangL-S, et al. (2011) Identification of common variants influencing risk of the tauopathy progressive supranuclear palsy. Nat Genet 43: 699–705 doi:10.1038/ng.859
18. ConradC, AndreadisA, TrojanowskiJQ, DicksonDW, KangD, et al. (1997) Genetic evidence for the involvement of tau in progressive supranuclear palsy. Ann Neurol 41: 277–281.
19. Wade-MartinsR (2012) Genetics: The MAPT locus-a genetic paradigm in disease susceptibility. Nat Rev Neurol 8: 477–478 doi:10.1038/nrneurol.2012.169
20. KalinderiK, FidaniL, BostantjopoulouS (2009) From 1997 to 2007: a decade journey through the H1 haplotype on 17q21 chromosome. Parkinsonism Relat Disord 15: 2–5 doi:10.1016/j.parkreldis.2008.03.001
21. PittmanAM, MyersAJ, Abou-SleimanP, FungHC, KaleemM, et al. (2005) Linkage disequilibrium fine mapping and haplotype association analysis of the tau gene in progressive supranuclear palsy and corticobasal degeneration. J Med Genet 42: 837–846 doi:10.1136/jmg.2005.031377
22. KauweJSK, CruchagaC, MayoK, FenoglioC, BertelsenS, et al. (2008) Variation in MAPT is associated with cerebrospinal fluid tau levels in the presence of amyloid-beta deposition. Proc Natl Acad Sci U S A 105: 8050–8054 doi:10.1073/pnas.0801227105
23. GrömpingU (2006) Relative Importance for Linear Regression in R: The Package relaimpo. J Stat Softw 17: 1–27.
24. HousemanEA, AccomandoWP, KoestlerDC, ChristensenBC, MarsitCJ, et al. (2012) DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics 13: 86 doi:10.1186/1471-2105-13-86
25. KoestlerDC, ChristensenB, KaragasMR, MarsitCJ, LangevinSM, et al. (2013) Blood-based profiles of DNA methylation predict the underlying distribution of cell types: a validation analysis. Epigenetics 8: 816–826 doi:10.4161/epi.25430
26. LiuY, AryeeMJ, PadyukovL, FallinMD, HesselbergE, et al. (2013) Epigenome-wide association data implicate DNA methylation as an intermediary of genetic risk in rheumatoid arthritis. Nat Biotechnol 31: 142–147 doi:10.1038/nbt.2487
27. EvansW, FungHC, SteeleJ, EerolaJ, TienariP, et al. (2004) The tau H2 haplotype is almost exclusively Caucasian in origin. Neurosci Lett 369: 183–185 doi:10.1016/j.neulet.2004.05.119
28. HeynH, LiN, FerreiraHJ, MoranS, PisanoDG, et al. (2012) Distinct DNA methylomes of newborns and centenarians. Proc Natl Acad Sci U S A 109: 10522–10527 doi:10.1073/pnas.1120658109
29. GibbsJR, van der BrugMP, HernandezDG, TraynorBJ, NallsMa, et al. (2010) Abundant quantitative trait loci exist for DNA methylation and gene expression in human brain. PLoS Genet 6: e1000952 doi:10.1371/journal.pgen.1000952
30. ZhangD, ChengL, BadnerJA, ChenC, ChenQ, et al. (2010) Genetic control of individual differences in gene-specific methylation in human brain. Am J Hum Genet 86: 411–419 doi:10.1016/j.ajhg.2010.02.005
31. SchadtEE, LambJ, YangX, ZhuJ, EdwardsS, et al. (2005) An integrative genomics approach to infer causal associations between gene expression and disease. Nat Genet 37: 710–717 doi:10.1038/ng1589
32. ZhuJ, WienerMC, ZhangC, FridmanA, MinchE, et al. (2007) Increasing the power to detect causal associations by combining genotypic and expression data in segregating populations. PLoS Comput Biol 3: e69 doi:10.1371/journal.pcbi.0030069
33. Pearl J (2009) Causality, 2nd Edition. Cambridge: Cambridge University Press.
34. VansteelandtS, LangeC (2012) Causation and causal inference for genetic effects. Hum Genet 1665–1676 doi:10.1007/s00439-012-1208-9
35. AtenJE, FullerTF, LusisAJ, HorvathS (2008) Using genetic markers to orient the edges in quantitative trait networks: the NEO software. BMC Syst Biol 2: 34 doi:10.1186/1752-0509-2-34
36. BakerM, LitvanI, HouldenH, AdamsonJ, DicksonD, et al. (1999) Association of an extended haplotype in the tau gene with progressive supranuclear palsy. Hum Mol Genet 8: 711–715.
37. CaffreyTM, Wade-MartinsR (2012) The role of MAPT sequence variation in mechanisms of disease susceptibility. Biochem Soc Trans 40: 687–692 doi:10.1042/BST20120063
38. PittmanAM, MyersAJ, DuckworthJ, BrydenL, HansonM, et al. (2004) The structure of the tau haplotype in controls and in progressive supranuclear palsy. Hum Mol Genet 13: 1267–1274 doi:10.1093/hmg/ddh138
39. StefanssonH, HelgasonA, ThorleifssonG, SteinthorsdottirV, MassonG, et al. (2005) A common inversion under selection in Europeans. Nat Genet 37: 129–137 doi:10.1038/ng1508
40. RademakersR, MelquistS, CrutsM, TheunsJ, Del-FaveroJ, et al. (2005) High-density SNP haplotyping suggests altered regulation of tau gene expression in progressive supranuclear palsy. Hum Mol Genet 14: 3281–3292 doi:10.1093/hmg/ddi361
41. MyersAJ, PittmanAM, ZhaoAS, RohrerK, KaleemM, et al. (2007) The MAPT H1c risk haplotype is associated with increased expression of tau and especially of 4 repeat containing transcripts. Neurobiol Dis 25: 561–570 doi:10.1016/j.nbd.2006.10.018
42. HayesmooreJB, BrayNJ, CrossWC, OwenMJ, O'DonovanMC, et al. (2009) The effect of age and the H1c MAPT haplotype on MAPT expression in human brain. Neurobiol Aging 30: 1652–1656 doi:10.1016/j.neurobiolaging.2007.12.017
43. TrabzuniD, WrayS, VandrovcovaJ, RamasamyA, WalkerR, et al. (2012) MAPT expression and splicing is differentially regulated by brain region: relation to genotype and implication for tauopathies. Hum Mol Genet 1–10 doi:10.1093/hmg/dds238
44. CaffreyTM, JoachimC, Wade-MartinsR (2008) Haplotype-specific expression of the N-terminal exons 2 and 3 at the human MAPT locus. Neurobiol Aging 29: 1923–1929 doi:10.1016/j.neurobiolaging.2007.05.002
45. CaffreyTM, JoachimC, ParacchiniS, EsiriMM, Wade-MartinsR (2006) Haplotype-specific expression of exon 10 at the human MAPT locus. Hum Mol Genet 15: 3529–3537 doi:10.1093/hmg/ddl429
46. Simón-SánchezJ, SchulteC, BrasJM, SharmaM, GibbsJR, et al. (2009) Genome-wide association study reveals genetic risk underlying Parkinson's disease. Nat Genet 41: 1308–1312 doi:10.1038/ng.487
47. Vilariño-GüellC, Soto-OrtolazaAI, RajputA, MashDC, PapapetropoulosS, et al. (2011) MAPT H1 haplotype is a risk factor for essential tremor and multiple system atrophy. Neurology 76: 670–672 doi:10.1212/WNL.0b013e31820c30c1
48. LiuY, AryeeMJ, PadyukovL, FallinMD, HesselbergE, et al. (2013) Epigenome-wide association data implicate DNA methylation as an intermediary of genetic risk in rheumatoid arthritis. Nat Biotechnol doi:10.1038/nbt.2487
49. GhoshS, YatesAJ, FrühwaldMC, MiecznikowskiJC, PlassC, et al. (2010) Tissue specific DNA methylation of CpG islands in normal human adult somatic tissues distinguishes neural from non-neural tissues. Epigenetics 5: 527–538.
50. DaviesMN, VoltaM, PidsleyR, LunnonK, DixitA, et al. (2012) Functional annotation of the human brain methylome identifies tissue-specific epigenetic variation across brain and blood. Genome Biol 13: R43 doi:10.1186/gb-2012-13-6-r43
51. CoppolaG, KarydasA, RademakersR, WangQ, BakerM, et al. (2008) Gene expression study on peripheral blood identifies progranulin mutations. Ann Neurol 64: 92–96 doi:10.1002/ana.21397
52. CoppolaG, BurnettR, PerlmanS, VersanoR, GaoF, et al. (2011) A gene expression phenotype in lymphocytes from Friedreich ataxia patients. Ann Neurol 70: 790–804 doi:10.1002/ana.22526
53. CoppolaG, MillerB, ChuiH, VarpetianA, LeveyA, et al. (2007) Genetic Investigation in Frontotemporal Dementia and Alzheimer's Disease: the GIFT Study. Ann Neurol 62.
54. GuH, SmithZD, BockC, BoyleP, GnirkeA, et al. (2011) Preparation of reduced representation bisulfite sequencing libraries for genome-scale DNA methylation profiling. Nat Protoc 6: 468–481 doi:10.1038/nprot.2010.190
55. GuoW, FizievP, YanW, CokusS, SunX, et al. (2013) BS-Seeker2: a versatile aligning pipeline for bisulfite sequencing data. BMC Genomics 14: 774 doi:10.1186/1471-2164-14-774
56. Hansen KD, Aryee M (n.d.) minfi: Analyze Illumina's 450 k methylation arrays. R package version 1.0.0.
57. MaksimovicJ, GordonL, OshlackA (2012) SWAN: Subset-quantile within array normalization for illumina infinium HumanMethylation450 BeadChips. Genome Biol 13: R44 doi:10.1186/gb-2012-13-6-r44
58. ChenY-A, LemireM, ChoufaniS, ButcherDT, GrafodatskayaD, et al. (2013) Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 microarray. Epigenetics 8: 203–209.
59. Smyth GK (2005) Limma: linear models for microarray data. In: Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W, editors. Bioinformatics and Computational Biology Solutions using R and Bioconductor. Springer, New York. pp. 397–420.
60. BellJT, TsaiP-C, YangT-P, PidsleyR, NisbetJ, et al. (2012) Epigenome-Wide Scans Identify Differentially Methylated Regions for Age and Age-Related Phenotypes in a Healthy Ageing Population. PLoS Genet 8: e1002629 doi:10.1371/journal.pgen.1002629
61. JohnsonWE, LiC, RabinovicA (2007) Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics 8: 118–127 doi:10.1093/biostatistics/kxj037
62. PurcellS, NealeB, Todd-BrownK, ThomasL, FerreiraMAR, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575 doi:10.1086/519795
63. CoppolaG (2011) Designing, performing, and interpreting a microarray-based gene expression study. Methods Mol Biol 793: 417–439 doi:_10.1007/978-1-61779-328-8_28
64. Horvath S (2011) Weighted Network Analysis. Springer.
65. KrzywinskiM, ScheinJ, BirolI, ConnorsJ, GascoyneR, HorsmanD, JonesSJ, MarraMA (2009) Circos: an information aesthetic for comparative genomics. Genome res 19: 1639–1645.
Štítky
Genetika Reprodukční medicínaČlánek vyšel v časopise
PLOS Genetics
2014 Číslo 3
- Primární hyperoxalurie – aktuální možnosti diagnostiky a léčby
- Srdeční frekvence embrya může být faktorem užitečným v předpovídání výsledku IVF
- Akutní intermitentní porfyrie
- Vztah užívání alkoholu a mužské fertility
- Šanci na úspěšný průběh těhotenství snižují nevhodné hladiny progesteronu vznikající při umělém oplodnění
Nejčtenější v tomto čísle
- Worldwide Patterns of Ancestry, Divergence, and Admixture in Domesticated Cattle
- Genome-Wide DNA Methylation Analysis of Human Pancreatic Islets from Type 2 Diabetic and Non-Diabetic Donors Identifies Candidate Genes That Influence Insulin Secretion
- Genetic Dissection of Photoreceptor Subtype Specification by the Zinc Finger Proteins Elbow and No ocelli
- GC-Rich DNA Elements Enable Replication Origin Activity in the Methylotrophic Yeast