DNA methylation patterns in the frontal lobe white matter of multiple system atrophy, Parkinson's disease, and progressive supranuclear palsy: a cross-comparative investigation.

Multiple system atrophy (MSA) is a rare neurodegenerative disease characterized by neuronal loss and gliosis, with oligodendroglial cytoplasmic inclusions (GCIs) containing α-synuclein being the primary pathological hallmark. Clinical presentations of MSA overlap with other parkinsonian disorders, such as Parkinson's disease (PD), dementia with Lewy bodies (DLB), and progressive supranuclear palsy (PSP), posing challenges in early diagnosis. Numerous studies have reported alterations in DNA methylation in neurodegenerative diseases, with candidate loci being identified in various parkinsonian disorders including MSA, PD, and PSP. Although MSA and PSP present with substantial white matter pathology, alterations in white matter have also been reported in PD. However, studies comparing the DNA methylation architectures of white matter in these diseases are lacking. We therefore aimed to investigate genome-wide DNA methylation patterns in the frontal lobe white matter of individuals with MSA (n = 17), PD (n = 17), and PSP (n = 16) along with controls (n = 15) using the Illumina EPIC array, to identify shared and disease-specific DNA methylation alterations. Genome-wide DNA methylation profiling of frontal lobe white matter in the three parkinsonian disorders revealed substantial commonalities in DNA methylation alterations in MSA, PD, and PSP. We further used weighted gene correlation network analysis to identify disease-associated co-methylation signatures and identified dysregulation in processes relating to Wnt signaling, signal transduction, endoplasmic reticulum stress, mitochondrial processes, RNA interference, and endosomal transport to be shared between these parkinsonian disorders. Our overall analysis points toward more similarities in DNA methylation patterns between MSA and PD, both synucleinopathies, compared to that between MSA and PD with PSP, which is a tauopathy. Our results also highlight several shared DNA methylation changes and pathways indicative of converging molecular mechanisms in the white matter contributing toward neurodegeneration in all three parkinsonian disorders.

Epigenetic Age Acceleration in Frontotemporal Lobar Degeneration: A Comprehensive Analysis in the Blood and Brain.

Frontotemporal lobar degeneration (FTLD) includes a heterogeneous group of disorders pathologically characterized by the degeneration of the frontal and temporal lobes. In addition to major genetic contributors of FTLD such as mutations in MAPT, GRN, and C9orf72, recent work has identified several epigenetic modifications including significant differential DNA methylation in DLX1, and OTUD4 loci. As aging remains one of the major risk factors for FTLD, we investigated the presence of accelerated epigenetic aging in FTLD compared to controls. We calculated epigenetic age in both peripheral blood and brain tissues of multiple FTLD subtypes using several DNA methylation clocks, i.e., DNAmClockMulti, DNAmClockHannum, DNAmClockCortical, GrimAge, and PhenoAge, and determined age acceleration and its association with different cellular proportions and clinical traits. Significant epigenetic age acceleration was observed in the peripheral blood of both frontotemporal dementia (FTD) and progressive supranuclear palsy (PSP) patients compared to controls with DNAmClockHannum, even after accounting for confounding factors. A similar trend was observed with both DNAmClockMulti and DNAmClockCortical in post-mortem frontal cortex tissue of PSP patients and in FTLD cases harboring GRN mutations. Our findings support that increased epigenetic age acceleration in the peripheral blood could be an indicator for PSP and to a smaller extent, FTD.

Brain DNA methylomic analysis of frontotemporal lobar degeneration reveals OTUD4 in shared dysregulated signatures across pathological subtypes.

Frontotemporal lobar degeneration (FTLD) is an umbrella term describing the neuropathology of a clinically, genetically and pathologically heterogeneous group of diseases, including frontotemporal dementia (FTD) and progressive supranuclear palsy (PSP). Among the major FTLD pathological subgroups, FTLD with TDP-43 positive inclusions (FTLD-TDP) and FTLD with tau-positive inclusions (FTLD-tau) are the most common, representing about 90% of the cases. Although alterations in DNA methylation have been consistently associated with neurodegenerative diseases, including Alzheimer's disease and Parkinson's disease, little is known for FTLD and its heterogeneous subgroups and subtypes. The main goal of this study was to investigate DNA methylation variation in FTLD-TDP and FTLD-tau. We used frontal cortex genome-wide DNA methylation profiles from three FTLD cohorts (142 FTLD cases and 92 controls), generated using the Illumina 450K or EPIC microarrays. We performed epigenome-wide association studies (EWAS) for each cohort followed by meta-analysis to identify shared differentially methylated loci across FTLD subgroups/subtypes. In addition, we used weighted gene correlation network analysis to identify co-methylation signatures associated with FTLD and other disease-related traits. Wherever possible, we also incorporated relevant gene/protein expression data. After accounting for a conservative Bonferroni multiple testing correction, the EWAS meta-analysis revealed two differentially methylated loci in FTLD, one annotated to OTUD4 (5'UTR-shore) and the other to NFATC1 (gene body-island). Of these loci, OTUD4 showed consistent upregulation of mRNA and protein expression in FTLD. In addition, in the three independent co-methylation networks, OTUD4-containing modules were enriched for EWAS meta-analysis top loci and were strongly associated with the FTLD status. These co-methylation modules were enriched for genes implicated in the ubiquitin system, RNA/stress granule formation and glutamatergic synaptic signalling. Altogether, our findings identified novel FTLD-associated loci, and support a role for DNA methylation as a mechanism involved in the dysregulation of biological processes relevant to FTLD, highlighting novel potential avenues for therapeutic development.

Epigenetic age acceleration is associated with oligodendrocyte proportions in MSA and control brain tissue.

AIMS: Epigenetic clocks are widely applied as surrogates for biological age in different tissues and/or diseases, including several neurodegenerative diseases. Despite white matter (WM) changes often being observed in neurodegenerative diseases, no study has investigated epigenetic ageing in white matter. METHODS: We analysed the performances of two DNA methylation-based clocks, DNAmClockMulti and DNAmClockCortical , in post-mortem WM tissue from multiple subcortical regions and the cerebellum, and in oligodendrocyte-enriched nuclei. We also examined epigenetic ageing in control and multiple system atrophy (MSA) (WM and mixed WM and grey matter), as MSA is a neurodegenerative disease comprising pronounced WM changes and α-synuclein aggregates in oligodendrocytes. RESULTS: Estimated DNA methylation (DNAm) ages showed strong correlations with chronological ages, even in WM (e.g., DNAmClockCortical , r = [0.80-0.97], p < 0.05). However, performances and DNAm age estimates differed between clocks and brain regions. DNAmClockMulti significantly underestimated ages in all cohorts except in the MSA prefrontal cortex mixed tissue, whereas DNAmClockCortical tended towards age overestimations. Pronounced age overestimations in the oligodendrocyte-enriched cohorts (e.g., oligodendrocyte-enriched nuclei, p = 6.1 × 10-5 ) suggested that this cell type ages faster. Indeed, significant positive correlations were observed between estimated oligodendrocyte proportions and DNAm age acceleration estimated by DNAmClockCortical (r > 0.31, p < 0.05), and similar trends were obtained with DNAmClockMulti . Although increased age acceleration was observed in MSA compared with controls, no significant differences were detected upon adjustment for possible confounders (e.g., cell-type proportions). CONCLUSIONS: Our findings show that oligodendrocyte proportions positively influence epigenetic age acceleration across brain regions and highlight the need to further investigate this in ageing and neurodegeneration.

Neurodegenerative movement disorders: An epigenetics perspective and promise for the future.

Neurodegenerative movement disorders (NMDs) are age-dependent disorders that are characterised by the degeneration and loss of neurons, typically accompanied by pathological accumulation of different protein aggregates in the brain, which lead to motor symptoms. NMDs include Parkinson's disease, multiple system atrophy, progressive supranuclear palsy, and Huntington's disease, among others. Epigenetic modifications are responsible for functional gene regulation during development, adult life and ageing and have progressively been implicated in complex diseases such as cancer and more recently in neurodegenerative diseases, such as NMDs. DNA methylation is by far the most widely studied epigenetic modification and consists of the reversible addition of a methyl group to the DNA without changing the DNA sequence. Although this research field is still in its infancy in relation to NMDs, an increasing number of studies point towards a role for DNA methylation in disease processes. This review addresses recent advances in epigenetic and epigenomic research in NMDs, with a focus on human brain DNA methylation studies. We discuss the current understanding of the DNA methylation changes underlying these disorders, the potential for use of these DNA modifications in peripheral tissues as biomarkers in early disease detection, classification and progression as well as a promising role in future disease management and therapy.

Increased brain expression of GPNMB is associated with genome wide significant risk for Parkinson's disease on chromosome 7p15.3.

Genome wide association studies (GWAS) for Parkinson's disease (PD) have previously revealed a significant association with a locus on chromosome 7p15.3, initially designated as the glycoprotein non-metastatic melanoma protein B (GPNMB) locus. In this study, the functional consequences of this association on expression were explored in depth by integrating different expression quantitative trait locus (eQTL) datasets (Braineac, CAGEseq, GTEx, and Phenotype-Genotype Integrator (PheGenI)). Top risk SNP rs199347 eQTLs demonstrated increased expressions of GPNMB, KLHL7, and NUPL2 with the major allele (AA) in brain, with most significant eQTLs in cortical regions, followed by putamen. In addition, decreased expression of the antisense RNA KLHL7-AS1 was observed in GTEx. Furthermore, rs199347 is an eQTL with long non-coding RNA (AC005082.12) in human tissues other than brain. Interestingly, transcript-specific eQTLs in immune-related tissues (spleen and lymphoblastoid cells) for NUPL2 and KLHL7-AS1 were observed, which suggests a complex functional role of this eQTL in specific tissues, cell types at specific time points. Significantly increased expression of GPNMB linked to rs199347 was consistent across all datasets, and taken in combination with the risk SNP being located within the GPNMB gene, these results suggest that increased expression of GPNMB is the causative link explaining the association of this locus with PD. However, other transcript eQTLs and subsequent functional roles cannot be excluded. This highlights the importance of further investigations to understand the functional interactions between the coding genes, antisense, and non-coding RNA species considering the tissue and cell-type specificity to understand the underlying biological mechanisms in PD.

High-resolution arrays reveal burden of copy number variations on Parkinson disease genes associated with increased disease risk in random cohorts.

BACKGROUND: Parkinson disease (PD) is a neurological disease responsible for a considerable rate of mortality and morbidity in the society. Since the symptoms of the disease appear much later than the actual onset of neuron degeneration, a majority of cases remain undiagnosed until the manifestation of the symptoms. OBJECTIVES: In order to investigate the existence of such susceptibility in the population, we analyzed Copy Number Variation (CNV) influences on PD genes in 1715 individuals from 12 different populations. RESULTS: Overall, 16 CNV-PD genes, 3 known to be causal and 13 associated, were found to be significantly enriched. PARK2, was under heavy burden with ~1% of the population containing CNV in the exonic region. The impact of these genes on the genome and disease pathway was analyzed using several genome analysis tools. Protein interaction network of CNV-PD genes revealed a complex interaction of molecules forming a major hub by the α-Synuclein, whose direct interactors, LRRK2, PARK2 and ATP13A2 are under CNV influence. CONCLUSIONS: We hypothesize that CNVs may not be the initiating event in the pathogenesis of PD and remain latent until additional secondary hits are acquired and also propose novel genes that may fall under the PD pathway which contribute in pathogenesis.

Type 2 diabetes mellitus disease risk genes identified by genome wide copy number variation scan in normal populations.

AIMS: To identify the role of copy number variations (CNVs) on disease risk genes and its effect on disease phenotypes in type 2 diabetes mellitus (T2DM) in 12 random populations using high throughput arrays. METHODS: CNV analysis was carried out on a total of 1715 individuals from 12 populations, from ArrayExpress Archive of the European Bioinformatics Institute along with our subjects using Affymetrix Genome Wide SNP 6.0 array. CNV effect on T2DM genes were analyzed using several bioinformatics tools and a molecular protein interaction network was constructed to identify the disease mechanism altered by the CNVs. RESULTS: Analysis showed 34.4% of the total population to be under CNV burden for T2DM, with 83 disease causal and associated genes being under CNV influence. Hotspots were identified on chromosomes 22, 12, 6, 19 and 11.Overlap studies with case cohorts revealed significant disease risk genes such as EGFR, E2F1, PPP1R3A, HLA and TSPAN8. CONCLUSIONS: CNVs play a significant role in predisposing T2DM in normal cohorts and contribute to the phenotypic effects. Thus, CNVs should be considered as one of the major contributors in predisposition of the disease.

Global patterns of large copy number variations in the human genome reveal complexity in chromosome organization.

Global patterns of copy number variations (CNVs) in chromosomes are required to understand the dynamics of genome organization and complexity. For this study, analysis was performed using the Affymetrix Genome-Wide Human SNP Array 6.0 chip and CytoScan High-Density arrays. We identified a total of 44 109 CNVs from 1715 genomes with a mean of 25 CNVs in an individual, which established the first drafts of population-specific CNV maps providing a rationale for prioritizing chromosomal regions. About 19 905 ancient CNVs were identified across all chromosomes and populations at varying frequencies. CNV count, and sometimes CNV size, contributed to the bulk CNV size of the chromosome. Population specific lengthening and shortening of chromosomal length was observed. Sex bias for CNV presence was largely dependent on ethnicity. Lower CNV inheritance rate was observed for India, compared to YRI and CEU. A total of 33 candidate CNV hotspots from 5382 copy number (CN) variable region (CNVR) clusters were identified. Population specific CNV distribution patterns in p and q arms disturbed the assumption that CNV counts in the p arm are less common compared to long arms, and the CNV occurrence and distribution in chromosomes is length independent. This study unraveled the force of independent evolutionary dynamics on genome organization and complexity across chromosomes and populations.

Global spectrum of copy number variations reveals genome organizational plasticity and proposes new migration routes.

Global spectrum of CNVs is required to catalog variations to provide a high-resolution on the dynamics of genome-organization and human migration. In this study, we performed genome-wide genotyping using high-resolution arrays and identified 44,109 CNVs from 1,715 genomes across 12 populations. The study unraveled the force of independent evolutionary dynamics on genome-organizational plasticity across populations. We demonstrated the use of CNV tool to study human migration and identified a second major settlement establishing new migration routes in addition to existing ones.

Copy number variations burden on miRNA genes reveals layers of complexities involved in the regulation of pathways and phenotypic expression.

MicroRNAs are involved in post-transcriptional down-regulation of gene expression. Variations in miRNA genes can severely affect downstream-regulated genes and their pathways. However, population-specific burden of CNVs on miRNA genes and the complexities created towards the phenotype is not known. From a total of 44109 CNVs investigated from 1715 individuals across 12 populations using high-throughput arrays, 4007 miRNA-CNVs (∼ 9%) consisting 6542 (∼ 5%) miRNA genes with a total of 333 (∼ 5%) singleton miRNA genes were identified. We found miRNA-CNVs across the genomes of individuals showing multiple hits in many targets, co-regulated under the same pathway. This study proposes four mechanisms unraveling the many complexities in miRNA genes, targets and co-regulated miRNA genes towards establishment of phenotypic diversity.

Impact of copy number variations burden on coding genome in humans using integrated high resolution arrays.

Copy number variations (CNVs) alter the transcriptional and translational levels of genes by disrupting the coding structure and this burden of CNVs seems to be a significant contributor to phenotypic variations. Therefore it was necessary to assess the complexities of CNV burden on the coding genome. A total of 1715 individuals from 12 populations were used for CNV analysis in the present investigation. Analysis was performed using Affymetrix Genome-Wide Human SNP Array 6·0 chip and CytoScan High-Density arrays. CNVs were more frequently observed in the coding region than in the non-coding region. CNVs were observed vastly more frequently in the coding region than the non-coding region. CNVs were found to be enriched in the regions containing functional genes (83-96%) compared with the regions containing pseudogenes (4-17%). CNVs across the genome of an individual showed multiple hits across many genes, whose proteins interact physically and function under the same pathway. We identified varying numbers of proteins and degrees of interactions within protein complexes of single individual genomes. This study represents the first draft of a population-specific CNV genes map as well as a cross-populational map. The complex relationship of CNVs on genes and their physically interacting partners unravels many complexities involved in phenotype expression. This study identifies four mechanisms contributing to the complexities caused by the presence of multiple CNVs across many genes in the coding part of the genome.

Unravelling the complexity of human olfactory receptor repertoire by copy number analysis across population using high resolution arrays.

Olfactory receptors (OR), responsible for detection of odor molecules, belong to the largest family of genes and are highly polymorphic in nature having distinct polymorphisms associated with specific regions around the globe. Since there are no reports on the presence of copy number variations in OR repertoire of Indian population, the present investigation in 43 Indians along with 270 HapMap and 31 Tibetan samples was undertaken to study genome variability and evolution. Analysis was performed using Affymetrix Genome-Wide Human SNP Array 6.0 chip, Affymterix CytoScan(®) High-Density array, HD-CNV, and MAFFT program. We observed a total of 1527 OR genes in 503 CNV events from 81.3% of the study group, which includes 67.6% duplications and 32.4% deletions encompassing more of genes than pseudogenes. We report human genotypic variation in functional OR repertoire size across populations and it was found that the combinatorial effect of both "orthologous obtained from closely related species" and "paralogous derived sequences" provide the complexity to the continuously occurring OR CNVs.