SHP2 inhibitors maintain TGFß signalling through SMURF2 inhibition.

Despite the promising antitumor activity of SHP2 inhibitors in RAS-dependent tumours, overall responses have been limited by their narrow therapeutic window. Like with all MAPK pathway inhibitors, this is likely the result of compensatory pathway activation mechanisms. However, the underlying mechanisms of resistance to SHP2 inhibition remain unknown. The E3 ligase SMURF2 limits TGFß activity by ubiquitinating and targeting the TGFß receptor for proteosome degradation. Using a functional RNAi screen targeting all known phosphatases, we identify that the tyrosine phosphatase SHP2 is a critical regulator of TGFß activity. Specifically, SHP2 dephosphorylates two key residues on SMURF2, resulting in activation of the enzyme. Conversely, SHP2 depletion maintains SMURF2 in an inactive state, resulting in the maintenance of TGFß activity. Furthermore, we demonstrate that depleting SHP2 has significant implications on TGFß-mediated migration, senescence, and cell survival. These effects can be overcome through the use of TGFß-targeted therapies. Consequently, our findings provide a rationale for combining SHP2 and TGFß inhibitors to enhance tumour responses leading to improved patient outcomes.

Increased chromatin accessibility facilitates intron retention in specific cell differentiation states.

Dynamic intron retention (IR) in vertebrate cells is of widespread biological importance. Aberrant IR is associated with numerous human diseases including several cancers. Despite consistent reports demonstrating that intrinsic sequence features can help introns evade splicing, conflicting findings about cell type- or condition-specific IR regulation by trans-regulatory and epigenetic mechanisms demand an unbiased and systematic analysis of IR in a controlled experimental setting. We integrated matched mRNA sequencing (mRNA-Seq), whole-genome bisulfite sequencing (WGBS), nucleosome occupancy methylome sequencing (NOMe-Seq) and chromatin immunoprecipitation sequencing (ChIP-Seq) data from primary human myeloid and lymphoid cells. Using these multi-omics data and machine learning, we trained two complementary models to determine the role of epigenetic factors in the regulation of IR in cells of the innate immune system. We show that increased chromatin accessibility, as revealed by nucleosome-free regions, contributes substantially to the retention of introns in a cell-specific manner. We also confirm that intrinsic characteristics of introns are key for them to evade splicing. This study suggests an important role for chromatin architecture in IR regulation. With an increasing appreciation that pathogenic alterations are linked to RNA processing, our findings may provide useful insights for the development of novel therapeutic approaches that target aberrant splicing.

New insights into the genetic etiology of Alzheimer's disease and related dementias.

Characterization of the genetic landscape of Alzheimer's disease (AD) and related dementias (ADD) provides a unique opportunity for a better understanding of the associated pathophysiological processes. We performed a two-stage genome-wide association study totaling 111,326 clinically diagnosed/'proxy' AD cases and 677,663 controls. We found 75 risk loci, of which 42 were new at the time of analysis. Pathway enrichment analyses confirmed the involvement of amyloid/tau pathways and highlighted microglia implication. Gene prioritization in the new loci identified 31 genes that were suggestive of new genetically associated processes, including the tumor necrosis factor alpha pathway through the linear ubiquitin chain assembly complex. We also built a new genetic risk score associated with the risk of future AD/dementia or progression from mild cognitive impairment to AD/dementia. The improvement in prediction led to a 1.6- to 1.9-fold increase in AD risk from the lowest to the highest decile, in addition to effects of age and the APOE ε4 allele.

Early life affects late-life health through determining DNA methylation across the lifespan: A twin study.

BACKGROUND: Previous findings for the genetic and environmental contributions to DNA methylation variation were for limited age ranges only. We investigated the lifespan contributions and their implications for human health for the first time. METHODS: 1,720 monozygotic twin (MZ) pairs and 1,107 dizygotic twin (DZ) pairs aged 0-92 years were included. Familial correlations (i.e., correlations between twins) for 353,681 methylation sites were estimated and modelled as a function of twin pair cohabitation history. FINDINGS: The methylome average familial correlation was around zero at birth (MZ pair: -0.01; DZ pair: -0.04), increased with the time of twins living together during childhood at rates of 0.16 (95%CI: 0.12-0.20) for MZ pairs and 0.13 (95%CI: 0.07-0.20) for DZ pairs per decade, and decreased with the time of living apart during adulthood at rates of 0.026 (95%CI: 0.019-0.033) for MZ pairs and 0.027 (95%CI: 0.011-0.043) for DZ pairs per decade. Neither the increasing nor decreasing rate differed by zygosity (both P>0.1), consistent with cohabitation environment shared by twins, rather than genetic factors, influencing the methylation familial correlation changes. Familial correlations for 6.6% (23,386/353,681) sites changed with twin pair cohabitation history. These sites were enriched for high heritability, proximal promoters, and epigenetic/genetic associations with various early-life factors and late-life health conditions. INTERPRETATION: Early life strongly influences DNA methylation variation across the lifespan, and the effects are stronger for heritable sites and sites biologically relevant to the regulation of gene expression. Early life could affect late-life health through influencing DNA methylation. FUNDING: Victorian Cancer Agency, Cancer Australia, Cure Cancer Foundation.

Diagnostic and prognostic plasma biomarkers for preclinical Alzheimer's disease.

INTRODUCTION: This study involved a parallel comparison of the diagnostic and longitudinal monitoring potential of plasma glial fibrillary acidic protein (GFAP), total tau (t-tau), phosphorylated tau (p-tau181 and p-tau231), and neurofilament light (NFL) in preclinical Alzheimer's disease (AD). METHODS: Plasma proteins were measured using Simoa assays in cognitively unimpaired older adults (CU), with either absence (Aß-) or presence (Aß+) of brain amyloidosis. RESULTS: Plasma GFAP, t-tau, p-tau181, and p-tau231 concentrations were higher in Aß+ CU compared with Aß- CU cross-sectionally. GFAP had the highest effect size and area under the curve (AUC) in differentiating between Aß+ and Aß- CU; however, no statistically significant differences were observed between the AUCs of GFAP, p-tau181, and p-tau231, but all were significantly higher than the AUC of NFL, and the AUC of GFAP was higher than the AUC of t-tau. The combination of a base model (BM), comprising the AD risk factors, age, sex, and apolipoprotein E gene (APOE) ε4 status with GFAP was observed to have a higher AUC (>90%) compared with the combination of BM with any of the other proteins investigated in the current study. Longitudinal analyses showed increased GFAP and p-tau181 in Aß+ CU and increased NFL in Aß- CU, over a 12-month duration. GFAP, p-tau181, p-tau231, and NFL showed significant correlations with cognition, whereas no significant correlations were observed with hippocampal volume. DISCUSSION: These findings highlight the diagnostic and longitudinal monitoring potential of GFAP and p-tau for preclinical AD.

Identification of DNA methylation biomarkers with potential to predict response to neoadjuvant chemotherapy in triple-negative breast cancer.

Neoadjuvant chemotherapy (NAC) is used to treat triple-negative breast cancer (TNBC) prior to resection. Biomarkers that accurately predict a patient's response to NAC are needed to individualise therapy and avoid chemotoxicity from unnecessary chemotherapy. We performed whole-genome DNA methylation profiling on diagnostic TNBC biopsy samples from the Sequential Evaluation of Tumours Undergoing Preoperative (SETUP) NAC study. We found 9 significantly differentially methylated regions (DMRs) at diagnosis which were associated with response to NAC. We show that 4 of these DMRs are associated with TNBC overall survival (P < 0.05). Our results highlight the potential of DNA methylation biomarkers for predicting NAC response in TNBC.

DNA methylation is required to maintain both DNA replication timing precision and 3D genome organization integrity.

DNA replication timing and three-dimensional (3D) genome organization are associated with distinct epigenome patterns across large domains. However, whether alterations in the epigenome, in particular cancer-related DNA hypomethylation, affects higher-order levels of genome architecture is still unclear. Here, using Repli-Seq, single-cell Repli-Seq, and Hi-C, we show that genome-wide methylation loss is associated with both concordant loss of replication timing precision and deregulation of 3D genome organization. Notably, we find distinct disruption in 3D genome compartmentalization, striking gains in cell-to-cell replication timing heterogeneity and loss of allelic replication timing in cancer hypomethylation models, potentially through the gene deregulation of DNA replication and genome organization pathways. Finally, we identify ectopic H3K4me3-H3K9me3 domains from across large hypomethylated domains, where late replication is maintained, which we purport serves to protect against catastrophic genome reorganization and aberrant gene transcription. Our results highlight a potential role for the methylome in the maintenance of 3D genome regulation.

Investigating Olfactory Gene Variation and Odour Identification in Older Adults.

Ageing is associated with a decrease in odour identification. Additionally, deficits in olfaction have been linked to age-related disease and mortality. Heritability studies suggest genetic variation contributes to olfactory identification. The olfactory receptor (OR) gene family is the largest in the human genome and responsible for overall odour identification. In this study, we sought to find olfactory gene family variants associated with individual and overall odour identification and to examine the relationships between polygenic risk scores (PRS) for olfactory-related phenotypes and olfaction. Participants were Caucasian older adults from the Sydney Memory and Ageing Study and the Older Australian Twins Study with genome-wide genotyping data (n = 1395, mean age = 75.52 ± 6.45). The Brief-Smell Identification Test (BSIT) was administered in both cohorts. PRS were calculated from independent GWAS summary statistics for Alzheimer's disease (AD), white matter hyperintensities (WMH), Parkinson's disease (PD), hippocampal volume and smoking. Associations with olfactory receptor genes (n = 967), previously identified candidate olfaction-related SNPs (n = 36) and different PRS with BSIT scores (total and individual smells) were examined. All of the relationships were analysed using generalised linear mixed models (GLMM), adjusted for age and sex. Genes with suggestive evidence for odour identification were found for 8 of the 12 BSIT items. Thirteen out of 36 candidate SNPs previously identified from the literature were suggestively associated with several individual BSIT items but not total score. PRS for smoking, WMH and PD were negatively associated with chocolate identification. This is the first study to conduct genetic analyses with individual odorant identification, which found suggestive olfactory-related genes and genetic variants for multiple individual BSIT odours. Replication in independent and larger cohorts is needed.

Genome-wide association study of circulating interleukin 6 levels identifies novel loci.

Interleukin 6 (IL-6) is a multifunctional cytokine with both pro- and anti-inflammatory properties with a heritability estimate of up to 61%. The circulating levels of IL-6 in blood have been associated with an increased risk of complex disease pathogenesis. We conducted a two-staged, discovery and replication meta genome-wide association study (GWAS) of circulating serum IL-6 levels comprising up to 67 428 (ndiscovery = 52 654 and nreplication = 14 774) individuals of European ancestry. The inverse variance fixed effects based discovery meta-analysis, followed by replication led to the identification of two independent loci, IL1F10/IL1RN rs6734238 on chromosome (Chr) 2q14, (Pcombined = 1.8 × 10-11), HLA-DRB1/DRB5 rs660895 on Chr6p21 (Pcombined = 1.5 × 10-10) in the combined meta-analyses of all samples. We also replicated the IL6R rs4537545 locus on Chr1q21 (Pcombined = 1.2 × 10-122). Our study identifies novel loci for circulating IL-6 levels uncovering new immunological and inflammatory pathways that may influence IL-6 pathobiology.

Replication timing and epigenome remodelling are associated with the nature of chromosomal rearrangements in cancer.

DNA replication timing is known to facilitate the establishment of the epigenome, however, the intimate connection between replication timing and changes to the genome and epigenome in cancer remain largely uncharacterised. Here, we perform Repli-Seq and integrated epigenome analyses and demonstrate that genomic regions that undergo long-range epigenetic deregulation in prostate cancer also show concordant differences in replication timing. A subset of altered replication timing domains are conserved across cancers from different tissue origins. Notably, late-replicating regions in cancer cells display a loss of DNA methylation, and a switch in heterochromatin features from H3K9me3-marked constitutive to H3K27me3-marked facultative heterochromatin. Finally, analysis of 214 prostate and 35 breast cancer genomes reveal that late-replicating regions are prone to cis and early-replication to trans chromosomal rearrangements. Together, our data suggests that the nature of chromosomal rearrangement in cancer is related to the spatial and temporal positioning and altered epigenetic states of early-replicating compared to late-replicating loci.

A Meta-Analysis of Genome-Wide Association Studies of Growth Differentiation Factor-15 Concentration in Blood.

Blood levels of growth differentiation factor-15 (GDF-15), also known as macrophage inhibitory cytokine-1 (MIC-1), have been associated with various pathological processes and diseases, including cardiovascular disease and cancer. Prior studies suggest genetic factors play a role in regulating blood MIC-1/GDF-15 concentration. In the current study, we conducted the largest genome-wide association study (GWAS) to date using a sample of ∼5,400 community-based Caucasian participants, to determine the genetic variants associated with MIC-1/GDF-15 blood concentration. Conditional and joint (COJO), gene-based association, and gene-set enrichment analyses were also carried out to identify novel loci, genes, and pathways. Consistent with prior results, a locus on chromosome 19, which includes nine single nucleotide polymorphisms (SNPs) (top SNP, rs888663, p = 1.690 × 10-35), was significantly associated with blood MIC-1/GDF-15 concentration, and explained 21.47% of its variance. COJO analysis showed evidence for two independent signals within this locus. Gene-based analysis confirmed the chromosome 19 locus association and in addition, a putative locus on chromosome 1. Gene-set enrichment analyses showed that the"COPI-mediated anterograde transport" gene-set was associated with MIC-1/GDF15 blood concentration with marginal significance after FDR correction (p = 0.067). In conclusion, a locus on chromosome 19 was associated with MIC-1/GDF-15 blood concentration with genome-wide significance, with evidence for a new locus (chromosome 1). Future studies using independent cohorts are needed to confirm the observed associations especially for the chromosomes 1 locus, and to further investigate and identify the causal SNPs that contribute to MIC-1/GDF-15 levels.

Enduring epigenetic landmarks define the cancer microenvironment.

The growth and progression of solid tumors involves dynamic cross-talk between cancer epithelium and the surrounding microenvironment. To date, molecular profiling has largely been restricted to the epithelial component of tumors; therefore, features underpinning the persistent protumorigenic phenotype of the tumor microenvironment are unknown. Using whole-genome bisulfite sequencing, we show for the first time that cancer-associated fibroblasts (CAFs) from localized prostate cancer display remarkably distinct and enduring genome-wide changes in DNA methylation, significantly at enhancers and promoters, compared to nonmalignant prostate fibroblasts (NPFs). Differentially methylated regions associated with changes in gene expression have cancer-related functions and accurately distinguish CAFs from NPFs. Remarkably, a subset of changes is shared with prostate cancer epithelial cells, revealing the new concept of tumor-specific epigenome modifications in the tumor and its microenvironment. The distinct methylome of CAFs provides a novel epigenetic hallmark of the cancer microenvironment and promises new biomarkers to improve interpretation of diagnostic samples.

Acetylated histone variant H2A.Z is involved in the activation of neo-enhancers in prostate cancer.

Acetylation of the histone variant H2A.Z (H2A.Zac) occurs at active promoters and is associated with oncogene activation in prostate cancer, but its role in enhancer function is still poorly understood. Here we show that H2A.Zac containing nucleosomes are commonly redistributed to neo-enhancers in cancer resulting in a concomitant gain of chromatin accessibility and ectopic gene expression. Notably incorporation of acetylated H2A.Z nucleosomes is a pre-requisite for activation of Androgen receptor (AR) associated enhancers. H2A.Zac nucleosome occupancy is rapidly remodeled to flank the AR sites to initiate the formation of nucleosome-free regions and the production of AR-enhancer RNAs upon androgen treatment. Remarkably higher levels of global H2A.Zac correlate with poorer prognosis. Altogether these data demonstrate the novel contribution of H2A.Zac in activation of newly formed enhancers in prostate cancer.

Quantitative Performance Evaluator for Proteomics (QPEP): Web-based Application for Reproducible Evaluation of Proteomics Preprocessing Methods.

Tandem mass spectrometry is one of the most popular techniques for quantitation of proteomes. There exists a large variety of options in each stage of data preprocessing that impact the bias and variance of the summarized protein-level values. Using a newly released data set satisfying a replicated Latin squares design, a diverse set of performance metrics has been developed and implemented in a web-based application, Quantitative Performance Evaluator for Proteomics (QPEP). QPEP has the flexibility to allow users to apply their own method to preprocess this data set and share the results, allowing direct and straightforward comparison of new methodologies. Application of these new metrics to three case studies highlights that (i) the summarization of peptides to proteins is robust to the choice of peptide summary used, (ii) the differences between iTRAQ labels are stronger than the differences between experimental runs, and (iii) the commercial software ProteinPilot performs equivalently well at between-sample normalization to more complicated methods developed by academics. Importantly, finding (ii) underscores the benefits of using the principles of randomization and blocking to avoid the experimental measurements being confounded by technical factors. Data are available via ProteomeXchange with identifier PXD003608.

Novel genetic loci associated with hippocampal volume.

The hippocampal formation is a brain structure integrally involved in episodic memory, spatial navigation, cognition and stress responsiveness. Structural abnormalities in hippocampal volume and shape are found in several common neuropsychiatric disorders. To identify the genetic underpinnings of hippocampal structure here we perform a genome-wide association study (GWAS) of 33,536 individuals and discover six independent loci significantly associated with hippocampal volume, four of them novel. Of the novel loci, three lie within genes (ASTN2, DPP4 and MAST4) and one is found 200 kb upstream of SHH. A hippocampal subfield analysis shows that a locus within the MSRB3 gene shows evidence of a localized effect along the dentate gyrus, subiculum, CA1 and fissure. Further, we show that genetic variants associated with decreased hippocampal volume are also associated with increased risk for Alzheimer's disease (rg=-0.155). Our findings suggest novel biological pathways through which human genetic variation influences hippocampal volume and risk for neuropsychiatric illness.

Novel genetic loci underlying human intracranial volume identified through genome-wide association.

Intracranial volume reflects the maximally attained brain size during development, and remains stable with loss of tissue in late life. It is highly heritable, but the underlying genes remain largely undetermined. In a genome-wide association study of 32,438 adults, we discovered five previously unknown loci for intracranial volume and confirmed two known signals. Four of the loci were also associated with adult human stature, but these remained associated with intracranial volume after adjusting for height. We found a high genetic correlation with child head circumference (ρgenetic = 0.748), which indicates a similar genetic background and allowed us to identify four additional loci through meta-analysis (Ncombined = 37,345). Variants for intracranial volume were also related to childhood and adult cognitive function, and Parkinson's disease, and were enriched near genes involved in growth pathways, including PI3K-AKT signaling. These findings identify the biological underpinnings of intracranial volume and their link to physiological and pathological traits.

Increased Insulin following an Oral Glucose Load, Genetic Variation near the Melatonin Receptor MTNR1B, but No Biochemical Evidence of Endothelial Dysfunction in Young Asian Men and Women.

AIM: To identify biochemical and genetic variation relating to increased risk of developing type 2 diabetes mellitus and cardiovascular disease in young, lean male and female adults of different ethnicities. METHOD: Fasting blood and urine and non-fasting blood following oral glucose intake were analysed in 90 Caucasians, South Asians and South East/East Asians. RESULTS: There were no differences in age, birthweight, blood pressure, body mass index, percent body fat, total energy, percentage of macronutrient intake, microalbumin, leptin, cortisol, adrenocorticotropic hormone, nitric oxide metabolites, C-reactive protein, homocysteine, tumor necrosis factor-α, interleukin-6, von Willebrand factor, vascular cell adhesion molecule-1, plasminogen activator inhibitor-1, and tissue plasminogen activator. Fasting total cholesterol (P = .000), triglycerides (P = .050), low density lipoprotein (P = .009) and non-fasting blood glucose (15 min) (P = .024) were elevated in South Asians compared with Caucasians, but there was no significant difference in glucose area under curve (AUC). Non-fasting insulin in South Asians (15-120 min), in South East/East Asians (60-120 min), and insulin AUC in South Asians and South East/East Asians, were elevated compared with Caucasians (P≤0.006). The molar ratio of C-peptide AUC/Insulin AUC (P = .045) and adiponectin (P = .037) were lower in South Asians compared with Caucasians. A significant difference in allele frequency distributions in Caucasians and South Asians was found for rs2166706 (P = 0.022) and rs10830963 (P = 0.009), which are both near the melatonin receptor MTNR1B. CONCLUSIONS: Elevated non-fasting insulin exists in young South Asians of normal fasting glucose and insulin. Hepatic clearance of insulin may be reduced in South Asians. No current biochemical evidence exists of endothelial dysfunction at this stage of development. MTNR1B signalling may be a useful therapeutic target in Asian populations in the prevention of type 2 diabetes mellitus.

MiR-score: a novel 6-microRNA signature that predicts survival outcomes in patients with malignant pleural mesothelioma.

BACKGROUND: Prognosis of malignant pleural mesothelioma (MPM) is poor, and predicting the outcomes of treatment is difficult. Here we investigate the potential of microRNA expression to estimate prognosis of MPM patients. METHODS: Candidate microRNAs from microarray profiling of tumor samples from 8 long (median: 53.7 months) and 8 short (median: 6.4 months) survivors following extrapleural pneumonectomy (EPP) were validated by RT-qPCR in 48 additional EPP samples. Kaplan-Meier log ranking was used to further explore the association between microRNA expression and overall survival (OS). Binary logistic regression was used to construct a microRNA signature (miR-Score) that was able to predict an OS of ≥20 months. Performance of the miR-Score was evaluated by receiver operating characteristic (ROC) curve analysis and validated in a series of 43 tumor samples from patients who underwent palliative surgery [pleurectomy/decortication (P/D)]. RESULTS: The miR-Score, using expression data of six microRNAs (miR-21-5p, -23a-3p, -30e-5p, -221-3p, -222-3p, and -31-5p), enabled prediction of long survival with an accuracy of 92.3% for EPP and 71.9% for palliative P/D. Hazard ratios for score-negative patients were 4.12 (95% CI: 2.03-8.37) for EPP and 1.93 (95% CI: 1.01-3.69) for P/D. Importantly, adding the miR-Score to a set of clinical selection criteria (histology, age, gender) increased predictive accuracy in the independent validation set from 76.3% for clinical factors only to 87.3%. CONCLUSIONS: This study has identified a novel 6-microRNA signature (miR-Score) that can accurately predict prognosis of MPM patients.

Surface antigen profiles of leukocytes and melanoma cells in lymph node metastases are associated with survival in AJCC stage III melanoma patients.

There is an urgent need to identify more accurate prognostic biomarkers in melanoma patients, particularly in those with metastatic disease. This study aimed to identify melanoma and leukocyte surface antigens predictive of survival in a prospective series of AJCC stage IIIb/c melanoma patients (n = 29). Live cell suspensions were prepared from melanoma metastases within lymph nodes (LN). The suspensions were immuno-magnetically separated into CD45(+) (leukocyte) and CD45(-) (non-hematopoietic, enriched melanoma cell) fractions. Surface antigens on CD45(-) and CD45(+) cell populations were profiled using DotScan™ microarrays (Medsaic Pty. Ltd.) and showed differential abundance levels for 52 and 78 antigens respectively. Associations of the surface profiles with clinicopathologic and outcome data (median follow-up 35.4 months post LN resection) were sought using univariate (log-rank test) and multivariate (Wald's test; modelled with patient's age, gender and AJCC staging at LN recurrence) survival models. CD9 (p = 0.036), CD39 (p = 0.004) and CD55 (p = 0.005) on CD45(+) leukocytes were independently associated with distant metastasis-free survival using multivariate analysis. Leukocytes with high CD39 levels were also significantly associated with increased overall survival (OS) in multivariate analysis (p = 0.016). LNs containing leukocytes expressing CD11b (p = 0.025), CD49d (p = 0.043) and CD79b (p = 0.044) were associated with reduced OS on univariate analysis. For enriched melanoma cells (CD45(-) cell populations), 11 surface antigens were significantly correlated with the disease-free interval (DFI) between diagnosis of culprit primary melanoma and LN metastasis resection. Nine antigens on CD45(+) leukocytes also correlated with DFI. Following validation in independent datasets, surface markers identified here should enable more accurate determination of prognosis in stage III melanoma patients and provide better risk stratification of patients entering clinical trials.

DNA methylation of the MAPT gene in Parkinson's disease cohorts and modulation by vitamin E in vitro.

Parkinson's disease (PD) is a neurodegenerative disorder for which environmental factors influence disease risk and may act via an epigenetic mechanism. The microtubule-associated protein tau (MAPT) is a susceptibility gene for idiopathic PD. Methylation levels were determined by pyrosequencing of bisulfite-treated DNA in a leukocyte cohort (358 PD patients and 1084 controls) and in two brain cohorts (Brain1, comprising 69 cerebellum controls; and Brain2, comprising 3 brain regions from 28 PD patients and 12 controls). In vitro assays involved the transfection of methylated promoter-luciferase constructs or treatment with an exogenous micronutrient. In normal leukocytes, the MAPT H1/H2 diplotype and sex were predictors of MAPT methylation. Haplotype-specific pyrosequencing confirmed that the H1 haplotype had higher methylation than the H2 haplotype in normal leukocytes and brain tissues. MAPT methylation was negatively associated with MAPT expression in the Brain1 cohort and in transfected cells. Methylation levels differed between three normal brain regions (Brain2 cohort, putamen < cerebellum < anterior cingulate cortex). In PD samples, age at onset was positively associated with MAPT methylation in leukocytes. Moreover, there was hypermethylation in the cerebellum and hypomethylation in the putamen of PD patients compared with controls (Brain2 cohort). Finally, leukocyte methylation status was positively associated with blood vitamin E levels, and the effect was more significant in H2 haplotype carriers; this result was confirmed in cells that were exposed to 100 µM vitamin E. The significant effects of sex, diplotype, and brain region suggest that hypermethylation of the MAPT gene is neuroprotective by reducing MAPT expression. The effect of vitamin E on MAPT represents a possible gene-environment interaction.