24-Hour Rhythms of DNA Methylation and Their Relation with Rhythms of RNA Expression in the Human Dorsolateral Prefrontal Cortex

Circadian rhythms modulate the biology of many human tissues, including brain tissues, and are driven by a near 24-hour transcriptional feedback loop. These rhythms are paralleled by 24-hour rhythms of large portions of the transcriptome. The role of dynamic DNA methylation in influencing these rhythms is uncertain. While recent work in Neurospora suggests that dynamic site-specific circadian rhythms of DNA methylation may play a role in modulating the fungal molecular clock, such rhythms and their relationship to RNA expression have not, to our knowledge, been elucidated in mammalian tissues, including human brain tissues. We hypothesized that 24-hour rhythms of DNA methylation exist in the human brain, and play a role in driving 24-hour rhythms of RNA expression. We analyzed DNA methylation levels in post-mortem human dorsolateral prefrontal cortex samples from 738 subjects. We assessed for 24-hour rhythmicity of 420,132 DNA methylation sites throughout the genome by considering methylation levels as a function of clock time of death and parameterizing these data using cosine functions. We determined global statistical significance by permutation. We then related rhythms of DNA methylation with rhythms of RNA expression determined by RNA sequencing. We found evidence of significant 24-hour rhythmicity of DNA methylation. Regions near transcription start sites were enriched for high-amplitude rhythmic DNA methylation sites, which were in turn time locked to 24-hour rhythms of RNA expression of nearby genes, with the nadir of methylation preceding peak transcript expression by 1–3 hours. Weak ante-mortem rest-activity rhythms were associated with lower amplitude DNA methylation rhythms as were older age and the presence of Alzheimer's disease. These findings support the hypothesis that 24-hour rhythms of DNA methylation, particularly near transcription start sites, may play a role in driving 24-hour rhythms of gene expression in the human dorsolateral prefrontal cortex, and may be affected by age and Alzheimer's disease

Modeling the cumulative genetic risk for multiple sclerosis from genome-wide association data

Background: Multiple sclerosis (MS) is the most common cause of chronic neurologic disability beginning in early to middle adult life. Results from recent genome-wide association studies (GWAS) have substantially lengthened the list of disease loci and provide convincing evidence supporting a multifactorial and polygenic model of inheritance. Nevertheless, the knowledge of MS genetics remains incomplete, with many risk alleles still to be revealed. Methods: We used a discovery GWAS dataset (8,844 samples, 2,124 cases and 6,720 controls) and a multi-step logistic regression protocol to identify novel genetic associations. The emerging genetic profile included 350 independent markers and was used to calculate and estimate the cumulative genetic risk in an independent validation dataset (3,606 samples). Analysis of covariance (ANCOVA) was implemented to compare clinical characteristics of individuals with various degrees of genetic risk. Gene ontology and pathway enrichment analysis was done using the DAVID functional annotation tool, the GO Tree Machine, and the Pathway-Express profiling tool. Results: In the discovery dataset, the median cumulative genetic risk (P-Hat) was 0.903 and 0.007 in the case and control groups, respectively, together with 79.9% classification sensitivity and 95.8% specificity. The identified profile shows a significant enrichment of genes involved in the immune response, cell adhesion, cell communication/ signaling, nervous system development, and neuronal signaling, including ionotropic glutamate receptors, which have been implicated in the pathological mechanism driving neurodegeneration. In the validation dataset, the median cumulative genetic risk was 0.59 and 0.32 in the case and control groups, respectively, with classification sensitivity 62.3% and specificity 75.9%. No differences in disease progression or T2-lesion volumes were observed among four levels of predicted genetic risk groups (high, medium, low, misclassified). On the other hand, a significant difference (F = 2.75, P = 0.04) was detected for age of disease onset between the affected misclassified as controls (mean = 36 years) and the other three groups (high, 33.5 years; medium, 33.4 years; low, 33.1 years). Conclusions: The results are consistent with the polygenic model of inheritance. The cumulative genetic risk established using currently available genome-wide association data provides important insights into disease heterogeneity and completeness of current knowledge in MS genetics

Elevated DNA methylation across a 48-kb region spanning the HOXA gene cluster is associated with Alzheimer's disease neuropathology

Introduction Alzheimer's disease is a neurodegenerative disorder that is hypothesized to involve epigenetic dysregulation of gene expression in the brain. Methods We performed an epigenome-wide association study to identify differential DNA methylation associated with neuropathology in prefrontal cortex and superior temporal gyrus samples from 147 individuals, replicating our findings in two independent data sets (N = 117 and 740). Results We identify elevated DNA methylation associated with neuropathology across a 48-kb region spanning 208 CpG sites within the HOXA gene cluster. A meta-analysis of the top-ranked probe within the HOXA3 gene (cg22962123) highlighted significant hypermethylation across all three cohorts (P = 3.11 × 10−18). Discussion We present robust evidence for elevated DNA methylation associated with Alzheimer's disease neuropathology spanning the HOXA gene cluster on chromosome 7. These data add to the growing evidence highlighting a role for epigenetic variation in Alzheimer's disease, implicating the HOX gene family as a target for future investigation

O4–01–03: Genome‐wide association study and admixture mapping of age‐related cognitive decline in African‐Americans

Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/152592/1/alzjjalz201304336.pd

Complex relation of HLA-DRB1*1501, age at menarche, and age at multiple sclerosis onset

Objective: To examine the relationship between 2 markers of early multiple sclerosis (MS) onset, 1 genetic (HLA-DRB1*1501) and 1 experiential (early menarche), in 2 cohorts. Methods: We included 540 white women with MS or clinically isolated syndrome (N = 156 with genetic data available) and 1,390 white women without MS but with a first-degree relative with MS (Genes and Environment in Multiple Sclerosis [GEMS]). Age at menarche, HLA-DRB1*1501 status, and age at MS onset were analyzed. Results: In both cohorts, participants with at least 1 HLA-DRB1*1501 allele had a later age at menarche than did participants with no risk alleles (MS: mean difference = 0.49, 95% confidence interval [CI] = [0.03–0.95], p = 0.036; GEMS: mean difference = 0.159, 95% CI = [0.012–0.305], p = 0.034). This association remained after we adjusted for body mass index at age 18 (available in GEMS) and for other MS risk alleles, as well as a single nucleotide polymorphism near the HLA-A region previously associated with age of menarche (available in MS cohort). Confirming previously reported associations, in our MS cohort, every year decrease in age at menarche was associated with a 0.65-year earlier MS onset (95% CI = [0.07–1.22], p = 0.027, N = 540). Earlier MS onset was also found in individuals with at least 1 HLA-DRB1*1501 risk allele (mean difference = −3.40 years, 95% CI = [−6.42 to −0.37], p = 0.028, N = 156). Conclusions: In 2 cohorts, a genetic marker for earlier MS onset (HLA-DRB1*1501) was inversely related to earlier menarche, an experiential marker for earlier symptom onset. This finding warrants broader investigations into the association between the HLA region and hormonal regulation in determining the onset of autoimmune disease

Alzheimer’s loci: epigenetic associations and interaction with genetic factors

Objective: We explore the role of DNA methylation in Alzheimer’s disease (AD). To elucidate where DNA methylation falls along the causal pathway linking risk factors to disease, we examine causal models to assess its role in the pathology of AD. Methods: DNA methylation profiles were generated in 740 brain samples using the Illumina HumanMet450K beadset. We focused our analysis on CpG sites from 11 AD susceptibility gene regions. The primary outcome was a quantitative measure of neuritic amyloid plaque (NP), a key early element of AD pathology. We tested four causal models: (1) independent associations, (2) CpG mediating the association of a variant, (3) reverse causality, and (4) genetic variant by CpG interaction. Results: Six genes regions (17 CpGs) showed evidence of CpG associations with NP, independent of genetic variation – BIN1 (5), CLU (5), MS4A6A (3), ABCA7 (2), CD2AP (1), and APOE (1). Together they explained 16.8% of the variability in NP. An interaction effect was seen in the CR1 region for two CpGs, cg10021878 (P = 0.01) and cg05922028 (P = 0.001), in relation to NP. In both cases, subjects with the risk allele rs6656401AT/AA display more methylation being associated with more NP burden, whereas subjects with the rs6656401TT protective genotype have an inverse association with more methylation being associated with less NP. Interpretation These observations suggest that, within known AD susceptibility loci, methylation is related to pathologic processes of AD and may play a largely independent role by influencing gene expression in AD susceptibility loci

Genetic epistasis regulates amyloid deposition in resilient aging

AbstractIntroduction The brain-derived neurotrophic factor (BDNF) interacts with important genetic Alzheimer's disease (AD) risk factors. Specifically, variants within the SORL1 gene determine BDNF's ability to reduce amyloid β (Aβ) in vitro. We sought to test whether functional BDNF variation interacts with SORL1 genotypes to influence expression and downstream AD-related processes in humans. Methods We analyzed postmortem brain RNA sequencing and neuropathological data for 441 subjects from the Religious Orders Study/Memory and Aging Project and molecular and structural neuroimaging data for 1285 subjects from the Alzheimer's Disease Neuroimaging Initiative. Results We found one SORL1 RNA transcript strongly regulated by SORL1-BDNF interactions in elderly without pathological AD and showing stronger associations with diffuse than neuritic Aβ plaques. The same SORL1-BDNF interactions also significantly influenced Aβ load as measured with [18F]Florbetapir positron emission tomography. Discussion Our results bridge the gap between risk and resilience factors for AD, demonstrating interdependent roles of established SORL1 and BDNF functional genotypes

Modeling Disease Severity in Multiple Sclerosis Using Electronic Health Records

Objective: To optimally leverage the scalability and unique features of the electronic health records (EHR) for research that would ultimately improve patient care, we need to accurately identify patients and extract clinically meaningful measures. Using multiple sclerosis (MS) as a proof of principle, we showcased how to leverage routinely collected EHR data to identify patients with a complex neurological disorder and derive an important surrogate measure of disease severity heretofore only available in research settings. Methods: In a cross-sectional observational study, 5,495 MS patients were identified from the EHR systems of two major referral hospitals using an algorithm that includes codified and narrative information extracted using natural language processing. In the subset of patients who receive neurological care at a MS Center where disease measures have been collected, we used routinely collected EHR data to extract two aggregate indicators of MS severity of clinical relevance multiple sclerosis severity score (MSSS) and brain parenchymal fraction (BPF, a measure of whole brain volume). Results: The EHR algorithm that identifies MS patients has an area under the curve of 0.958, 83% sensitivity, 92% positive predictive value, and 89% negative predictive value when a 95% specificity threshold is used. The correlation between EHR-derived and true MSSS has a mean R[superscript 2] = 0.38±0.05, and that between EHR-derived and true BPF has a mean R[superscript 2] = 0.22±0.08. To illustrate its clinical relevance, derived MSSS captures the expected difference in disease severity between relapsing-remitting and progressive MS patients after adjusting for sex, age of symptom onset and disease duration (p = 1.56×10[superscript −12]). Conclusion: Incorporation of sophisticated codified and narrative EHR data accurately identifies MS patients and provides estimation of a well-accepted indicator of MS severity that is widely used in research settings but not part of the routine medical records. Similar approaches could be applied to other complex neurological disorders.National Institute of General Medical Sciences (U.S.) (NIH U54-LM008748

Modeling the cumulative genetic risk for multiple sclerosis from genome-wide association data

BACKGROUND: Multiple sclerosis (MS) is the most common cause of chronic neurologic disability beginning in early to middle adult life. Results from recent genome-wide association studies (GWAS) have substantially lengthened the list of disease loci and provide convincing evidence supporting a multifactorial and polygenic model of inheritance. Nevertheless, the knowledge of MS genetics remains incomplete, with many risk alleles still to be revealed. METHODS: We used a discovery GWAS dataset (8,844 samples, 2,124 cases and 6,720 controls) and a multi-step logistic regression protocol to identify novel genetic associations. The emerging genetic profile included 350 independent markers and was used to calculate and estimate the cumulative genetic risk in an independent validation dataset (3,606 samples). Analysis of covariance (ANCOVA) was implemented to compare clinical characteristics of individuals with various degrees of genetic risk. Gene ontology and pathway enrichment analysis was done using the DAVID functional annotation tool, the GO Tree Machine, and the Pathway-Express profiling tool. RESULTS: In the discovery dataset, the median cumulative genetic risk (P-Hat) was 0.903 and 0.007 in the case and control groups, respectively, together with 79.9% classification sensitivity and 95.8% specificity. The identified profile shows a significant enrichment of genes involved in the immune response, cell adhesion, cell communication/signaling, nervous system development, and neuronal signaling, including ionotropic glutamate receptors, which have been implicated in the pathological mechanism driving neurodegeneration. In the validation dataset, the median cumulative genetic risk was 0.59 and 0.32 in the case and control groups, respectively, with classification sensitivity 62.3% and specificity 75.9%. No differences in disease progression or T2-lesion volumes were observed among four levels of predicted genetic risk groups (high, medium, low, misclassified). On the other hand, a significant difference (F = 2.75, P = 0.04) was detected for age of disease onset between the affected misclassified as controls (mean = 36 years) and the other three groups (high, 33.5 years; medium, 33.4 years; low, 33.1 years). CONCLUSIONS: The results are consistent with the polygenic model of inheritance. The cumulative genetic risk established using currently available genome-wide association data provides important insights into disease heterogeneity and completeness of current knowledge in MS genetics