Prioritizing Parkinson's disease genes using population-scale transcriptomic data

Exploring the dynamic pathogenesis of Parkinson's disease by case-control and longitudinal blood transcriptome analyses

The pathogenesis of Parkinson's disease (PD) was recently hypothesized to change along with the disease course. Given the fact that transcriptional changes in blood can provide insightful clues for PD pathogenesis, we performed case-control and longitudinal whole blood transcriptome analyses to identify the signature genes underlying the hypothesized dynamic pathogenesis of PD. In the case-control study, we compared the gene expression patterns in healthy control (N = 189), prodromal (N = 58) and de novo idiopathic PD subjects (N = 390). The results showed that the prodromal subjects were at the tipping-point stage, which is characterized by the abnormal expression patterns of 414 genes associated with oxygen transport and reactive oxygen species metabolic process. We next performed a longitudinal transcriptome analysis on 255 PD patients from the baseline to the third year, and identified 203 genes related to immune and inflammatory responses during disease progression. These findings not just offer deeper insights into the dynamic pathogenesis of PD, but also help to find potential drugs to prevent the early neurodegeneration process.

Computational and functional prioritization identifies genes that rescue behavior and reduce tau protein in fly and human cell models of Alzheimer disease

Genome-wide association studies (GWASs) in Alzheimer disease (AD) have uncovered over 70 loci significantly associated with AD risk, but identifying the true causal gene(s) at these loci requires systematic functional validation that is rarely performed due to limitations of time and cost. Here, we integrate transcriptome-wide association study (TWAS) with colocalization analysis, fine-mapping, and additional annotation of AD GWAS variants to identify 123 genes at known and suggestive AD risk loci. A comparison with human AD brain transcriptome data confirmed that many of these candidate genes are dysregulated in human AD and correlate with neuropathology. We then tested all available orthologs in two well-established Drosophila AD models that express either wild-type tau or secreted β-amyloid (β42). Experimental perturbation of the 60 available candidates pinpointed 46 that modulated neuronal dysfunction in one or both fly models. The effects of 18 of these genes were concordant with the TWAS prediction, such that the direction of misexpression predicted to increase AD risk in humans exacerbated behavioral impairments in the AD fly models. Reversing the aberrant down- or upregulation of 11 of these genes (MTCH2, ELL, TAP2, HDC, DMWD, MYCL, SLC4A9, ABCA7, CSTF1, PTK2B, and CD2AP) proved neuroprotective in vivo. We further studied MTCH2 and found that it regulates steady-state tau protein levels in the Drosophila brain and reduces tau accumulation in human neural progenitor cells. This systematic, integrative approach effectively prioritizes genes at GWAS loci and reveals promising AD-relevant candidates for further investigation as risk factors or targets for therapeutic intervention.

Identification of Common Brain Protein and Genetic Loci Between Parkinson's Disease and Lewy Body Dementia

BACKGROUND

Parkinson's disease (PD) and Lewy body dementia (LBD) have many common features, including clinical manifestations, neurochemistry, and pathology, but little is known about their shared brain proteins and genetic factors.

METHODS

To identify susceptibility-related brain proteins that are shared between PD and LBD patients, proteome-wide association studies (PWASs) were conducted by integrating human brain protein quantitative trait loci (pQTLs) with large-scale genome-wide association studies (GWASs) of both diseases. Subsequently, pleiotropy-informed conditional false discovery rate (pleioFDR) analysis was performed to identify common risk genetic loci between PD and LBD. Finally, the downregulation of these risk genes in different disease states was validated by differential gene expression analysis.

RESULTS

PWASs identified 12 PD risk proteins and nine LBD risk proteins, among which TMEM175 (zPD = -7.25, PPD = 4.12E-13; zLBD = -6.02, PLBD = 1.75E-09) and DOC2A (zPD = -4.13, PPD = 3.71E-05; zLBD = -3.91, PLBD = 9.08E-05) were shared. PleioFDR analysis revealed that five genetic risk loci mapped to eight genes associated with PD and LBD, including the proteome-wide significant risk gene TMEM175 (ConjFDR = 5.74E-03). Differential expression analysis verified that TMEM175 was significantly downregulated in the midbrains of PD patients (p = 1.19E-02), and further exploration revealed that TMEM175 was also dramatically downregulated in the substantia nigra of PD patients (p = 1.16E-02) and incidental Lewy body disease patients (p = 7.52E-03). Moreover, TMEM175 was significantly downregulated in induced pluripotent stem cell-derived dopaminergic neurons from PD patients (p = 4.60E-02).

CONCLUSION

Dysregulation of TMEM175 may confer PD and LBD risk and may be partly responsible for their comorbidity. Our results revealed the common genetic risk factors between PD and LBD, which elucidated the shared genetic basis of these diseases.

GPNMB and ATP6V1A interact to mediate microglia phagocytosis of multiple types of pathological particles

Pronounced elevation of glycoprotein non-metastatic melanoma B (GPNMB) is a common phenomenon in a variety of brain diseases, but the expression patterns, functions, and molecular signaling of GPNMB have not been well studied. Here, we showed that pathological factors, including neuronal degeneration caused by seizures, caspase-3-induced neuronal apoptosis, neuronal debris, and β-amyloid, induced "on-demand" GPNMB expression in hippocampal microglia. Genetic ablation of GPNMB did not affect acute seizures but worsened chronic epileptogenesis. We found that GPNMB functioned in phagocytosis, deficiency of which resulted in defects in both phagocytic engulfment and degradation. GPNMB could be internalized into cells, where it wrapped engulfed pathogenic particles and presented them to lysosomes through interaction with lysosomal vacuolar-type proton ATPase catalytic subunit A (ATP6V1A). Activating ATP6V1A was able to rescue GPNMB-deficiency-caused phagocytosis impairment. Thus, microglial GPNMB-ATP6V1A might be a common treatment target of a batch of chronic neurological disorders, and clearing the degenerative neurons might be more valuable than reserving them to protect the brain.

Long-read RNA sequencing atlas of human microglia isoforms elucidates disease-associated genetic regulation of splicing

Microglia, the innate immune cells of the central nervous system, have been genetically implicated in multiple neurodegenerative diseases. Mapping the genetics of gene expression in human microglia has identified several loci associated with disease-associated genetic variants in microglia-specific regulatory elements. However, identifying genetic effects on splicing is challenging because of the use of short sequencing reads. Here, we present the isoform-centric microglia genomic atlas (isoMiGA), which leverages long-read RNA sequencing to identify 35,879 novel microglia isoforms. We show that these isoforms are involved in stimulation response and brain region specificity. We then quantified the expression of both known and novel isoforms in a multi-ancestry meta-analysis of 555 human microglia short-read RNA sequencing samples from 391 donors, and found associations with genetic risk loci in Alzheimer's and Parkinson's disease. We nominate several loci that may act through complex changes in isoform and splice-site usage.

Single-cell RNA sequencing reveals peripheral immunological features in Parkinson's Disease

Although many researchers of Parkinson's disease (PD) have shifted their focus from the central nervous system (CNS) to the peripheral blood, a significant knowledge gap remains between PD severity and the peripheral immune response. In the current study, we aimed to map the peripheral immunity atlas in peripheral blood mononuclear cells (PBMCs) from PD patients and healthy controls using single-cell RNA sequencing (scRNA-seq). Our study employed scRNA-seq analysis to map the peripheral immunity atlas in PD by profiling PBMCs from PD-Early, PD-Late patients and matched controls. By enlarging the blood sample size, we validated the roles of NK cells in numerous immune-related biological processes. We also detected the infiltration of NK cells into the cerebral motor cortex as the disease progressed, using human brain sections, and elucidated the communication between the periphery and CNS and its implications for PD. As a result, cell subpopulation atlases in PBMCs from PD patients and healthy controls along with differentially expressed genes in NK cells were identified by scRNA-seq analysis, representing 6 major immune cell subsets among which NK cells declined in the progression of PD. We further validated NK cell reduction in increasing samples and found that they participated in numerous immune-related biological processes and infiltration into the cerebral motor cortex as the disease proceeded, evidencingd the close communication between the peripheral immune response and CNS. Strikingly, XCL2 positively correlated with PD severity, with good predictive performance of PD and specific expression in subclusters C2 and C5 of NK cells. All these findings delineated the critical role of peripheral immune response mediated by NK cells in the pathogenesis of PD. NK cell-specific XCL2 could be used as a diagnostic marker for treating PD. The indispensable function of NK cells and NK cell-specific molecular biomarkers highlighted the implication of the peripheral immune response in PD progression. Trial registration: ChiCTR, ChiCTR1900023975. Registered 20 June 2019 - Retrospectively registered, https://www.chictr.org.cn/showproj.html?proj=31035 .

Mediator complex in transcription regulation and DNA repair: Relevance for human diseases

The Mediator complex is an essential coregulator of RNA polymerase II transcription. More recent developments suggest Mediator functions as a link between transcription regulation, genome organisation and DNA repair mechanisms including nucleotide excision repair, base excision repair, and homologous recombination. Dysfunctions of these processes are frequently associated with human pathologies, and growing evidence shows Mediator involvement in cancers, neurological, metabolic and infectious diseases. The detailed deciphering of molecular mechanisms of Mediator functions, using interdisciplinary approaches in different biological models and considering all functions of this complex, will contribute to our understanding of relevant human diseases.

A systematic exploration of unexploited genes for oxidative stress in Parkinson's disease

Human disease-associated gene data are accessible through databases, including the Open Targets Platform, DisGeNET, miRTex, RNADisease, and PubChem. However, missing data entries in such databases are anticipated because of curational errors, biases, and text-mining failures. Additionally, the extensive research on human diseases has led to challenges in registering comprehensive data. The lack of essential data in databases hinders knowledge sharing and should be addressed. Therefore, we propose an analysis pipeline to explore missing entries of unexploited genes in the human disease-associated gene databases. Using this pipeline for genes in Parkinson's disease with oxidative stress revealed two unexploited genes: nuclear protein 1 (NUPR1) and ubiquitin-like with PHD and ring finger domains 2 (UHRF2). This methodology enhances the identification of underrepresented disease-associated genes, facilitating easier access to potential human disease-related functional genes. This study aims to identify unexploited genes for further research and does not include independent experimental validation.

LRRK2 G2019S variant is associated with transcriptional changes in Parkinson's disease human myeloid cells under proinflammatory environment

The G2019S mutation in the leucine-rich repeat kinase 2 (LRRK2) gene is a major risk factor for the development of Parkinson's disease (PD). LRRK2, although ubiquitously expressed, is highly abundant in cells of the innate immune system. Given the importance of central and peripheral immune cells in the development of PD, we sought to investigate the consequences of the G2019S mutation on microglial and monocyte transcriptome and function. We have generated large-scale transcriptomic profiles of isogenic human induced microglial cells (iMGLs) and patient derived monocytes carrying the G2019S mutation under baseline culture conditions and following exposure to the proinflammatory factors IFNγ and LPS. We demonstrate that the G2019S mutation exerts a profound impact on the transcriptomic profile of these myeloid cells, and describe corresponding functional differences in iMGLs. The G2019S mutation led to an upregulation in lipid metabolism and phagolysosomal pathway genes in untreated and LPS/IFNγ stimulated iMGLs, which was accompanied by an increased phagocytic capacity of myelin debris. We also identified dysregulation of cell cycle genes, with a downregulation of the E2F4 regulon. Transcriptomic characterization of human-derived monocytes carrying the G2019S mutation confirmed alteration in lipid metabolism associated genes. Altogether, these findings reveal the influence of G2019S on the dysregulation of the myeloid cell transcriptome under proinflammatory conditions.

Cross-ancestry atlas of gene, isoform, and splicing regulation in the developing human brain

Neuropsychiatric genome-wide association studies (GWASs), including those for autism spectrum disorder and schizophrenia, show strong enrichment for regulatory elements in the developing brain. However, prioritizing risk genes and mechanisms is challenging without a unified regulatory atlas. Across 672 diverse developing human brains, we identified 15,752 genes harboring gene, isoform, and/or splicing quantitative trait loci, mapping 3739 to cellular contexts. Gene expression heritability drops during development, likely reflecting both increasing cellular heterogeneity and the intrinsic properties of neuronal maturation. Isoform-level regulation, particularly in the second trimester, mediated the largest proportion of GWAS heritability. Through colocalization, we prioritized mechanisms for about 60% of GWAS loci across five disorders, exceeding adult brain findings. Finally, we contextualized results within gene and isoform coexpression networks, revealing the comprehensive landscape of transcriptome regulation in development and disease.

Transcriptomic imputation of genetic risk variants uncovers novel whole-blood biomarkers of Parkinson's disease

Blood-based gene expression signatures could potentially be used as biomarkers for PD. However, it is unclear whether genetically-regulated transcriptomic signatures can provide novel gene candidates for use as PD biomarkers. We leveraged on the Genotype-Tissue Expression (GTEx) database to impute whole-blood transcriptomic expression using summary statistics of three large-scale PD GWAS. A random forest classifier was used with the consensus whole-blood imputed gene signature (IGS) to discriminate between cases and controls. Outcome measures included Area under the Curve (AUC) of Receiver Operating Characteristic (ROC) Curve. We demonstrated that the IGS (n = 37 genes) is conserved across PD GWAS studies and brain tissues. IGS discriminated between cases and controls in an independent whole-blood RNA-sequencing study (1176 PD, 254 prodromal, and 860 healthy controls) with mean AUC and accuracy of 64.8% and 69.4% for PD cohort, and 78.8% and 74% for prodromal cohort. PATL2 was the top-performing imputed gene in both PD and prodromal PD cohorts, whose classifier performance varied with biological sex (higher performance for males and females in the PD and prodromal PD, respectively). Single-cell RNA-sequencing studies (scRNA-seq) of healthy humans and PD patients found PATL2 to be enriched in terminal effector CD8+ and cytotoxic CD4+ cells, whose proportions are both increased in PD patients. We demonstrated the utility of GWAS transcriptomic imputation in identifying novel whole-blood transcriptomic signatures which could be leveraged upon for PD biomarker derivation. We identified PATL2 as a potential biomarker in both clinical and prodromic PD.

Prioritizing susceptibility genes for the prognosis of male-pattern baldness with transcriptome-wide association study

BACKGROUND

Male-pattern baldness (MPB) is the most common cause of hair loss in men. It can be categorized into three types: type 2 (T2), type 3 (T3), and type 4 (T4), with type 1 (T1) being considered normal. Although various MPB-associated genetic variants have been suggested, a comprehensive study for linking these variants to gene expression regulation has not been performed to the best of our knowledge.

RESULTS

In this study, we prioritized MPB-related tissue panels using tissue-specific enrichment analysis and utilized single-tissue panels from genotype-tissue expression version 8, as well as cross-tissue panels from context-specific genetics. Through a transcriptome-wide association study and colocalization analysis, we identified 52, 75, and 144 MPB associations for T2, T3, and T4, respectively. To assess the causality of MPB genes, we performed a conditional and joint analysis, which revealed 10, 11, and 54 putative causality genes for T2, T3, and T4, respectively. Finally, we conducted drug repositioning and identified potential drug candidates that are connected to MPB-associated genes.

CONCLUSIONS

Overall, through an integrative analysis of gene expression and genotype data, we have identified robust MPB susceptibility genes that may help uncover the underlying molecular mechanisms and the novel drug candidates that may alleviate MPB.

Joint-tissue integrative analysis identifies high-risk genes for Parkinson's disease

The loss of dopaminergic neurons in the substantia nigra and the abnormal accumulation of synuclein proteins and neurotransmitters in Lewy bodies constitute the primary symptoms of Parkinson's disease (PD). Besides environmental factors, scholars are in the early stages of comprehending the genetic factors involved in the pathogenic mechanism of PD. Although genome-wide association studies (GWAS) have unveiled numerous genetic variants associated with PD, precisely pinpointing the causal variants remains challenging due to strong linkage disequilibrium (LD) among them. Addressing this issue, expression quantitative trait locus (eQTL) cohorts were employed in a transcriptome-wide association study (TWAS) to infer the genetic correlation between gene expression and a particular trait. Utilizing the TWAS theory alongside the enhanced Joint-Tissue Imputation (JTI) technique and Mendelian Randomization (MR) framework (MR-JTI), we identified a total of 159 PD-associated genes by amalgamating LD score, GTEx eQTL data, and GWAS summary statistic data from a substantial cohort. Subsequently, Fisher's exact test was conducted on these PD-associated genes using 5,152 differentially expressed genes sourced from 12 PD-related datasets. Ultimately, 29 highly credible PD-associated genes, including CTX1B, SCNA, and ARSA, were uncovered. Furthermore, GO and KEGG enrichment analyses indicated that these genes primarily function in tissue synthesis, regulation of neuron projection development, vesicle organization and transportation, and lysosomal impact. The potential PD-associated genes identified in this study not only offer fresh insights into the disease's pathophysiology but also suggest potential biomarkers for early disease detection.

Genetic imputation of kidney transcriptome, proteome and multi-omics illuminates new blood pressure and hypertension targets

Genetic mechanisms of blood pressure (BP) regulation remain poorly defined. Using kidney-specific epigenomic annotations and 3D genome information we generated and validated gene expression prediction models for the purpose of transcriptome-wide association studies in 700 human kidneys. We identified 889 kidney genes associated with BP of which 399 were prioritised as contributors to BP regulation. Imputation of kidney proteome and microRNAome uncovered 97 renal proteins and 11 miRNAs associated with BP. Integration with plasma proteomics and metabolomics illuminated circulating levels of myo-inositol, 4-guanidinobutanoate and angiotensinogen as downstream effectors of several kidney BP genes (SLC5A11, AGMAT, AGT, respectively). We showed that genetically determined reduction in renal expression may mimic the effects of rare loss-of-function variants on kidney mRNA/protein and lead to an increase in BP (e.g., ENPEP). We demonstrated a strong correlation (r = 0.81) in expression of protein-coding genes between cells harvested from urine and the kidney highlighting a diagnostic potential of urinary cell transcriptomics. We uncovered adenylyl cyclase activators as a repurposing opportunity for hypertension and illustrated examples of BP-elevating effects of anticancer drugs (e.g. tubulin polymerisation inhibitors). Collectively, our studies provide new biological insights into genetic regulation of BP with potential to drive clinical translation in hypertension.

Machine learning nominates the inositol pathway and novel genes in Parkinson's disease

There are 78 loci associated with Parkinson's disease in the most recent genome-wide association study (GWAS), yet the specific genes driving these associations are mostly unknown. Herein, we aimed to nominate the top candidate gene from each Parkinson's disease locus and identify variants and pathways potentially involved in Parkinson's disease. We trained a machine learning model to predict Parkinson's disease-associated genes from GWAS loci using genomic, transcriptomic and epigenomic data from brain tissues and dopaminergic neurons. We nominated candidate genes in each locus and identified novel pathways potentially involved in Parkinson's disease, such as the inositol phosphate biosynthetic pathway (INPP5F, IP6K2, ITPKB and PPIP5K2). Specific common coding variants in SPNS1 and MLX may be involved in Parkinson's disease, and burden tests of rare variants further support that CNIP3, LSM7, NUCKS1 and the polyol/inositol phosphate biosynthetic pathway are associated with the disease. Functional studies are needed to further analyse the involvements of these genes and pathways in Parkinson's disease.

A distinct class of pan-cancer susceptibility genes revealed by an alternative polyadenylation transcriptome-wide association study

Alternative polyadenylation plays an important role in cancer initiation and progression; however, current transcriptome-wide association studies mostly ignore alternative polyadenylation when identifying putative cancer susceptibility genes. Here, we perform a pan-cancer 3' untranslated region alternative polyadenylation transcriptome-wide association analysis by integrating 55 well-powered (n > 50,000) genome-wide association studies datasets across 22 major cancer types with alternative polyadenylation quantification from 23,955 RNA sequencing samples across 7,574 individuals. We find that genetic variants associated with alternative polyadenylation are co-localized with 28.57% of cancer loci and contribute a significant portion of cancer heritability. We further identify 642 significant cancer susceptibility genes predicted to modulate cancer risk via alternative polyadenylation, 62.46% of which have been overlooked by traditional expression- and splicing- studies. As proof of principle validation, we show that alternative alleles facilitate 3' untranslated region lengthening of CRLS1 gene leading to increased protein abundance and promoted proliferation of breast cancer cells. Together, our study highlights the significant role of alternative polyadenylation in discovering new cancer susceptibility genes and provides a strong foundational framework for enhancing our understanding of the etiology underlying human cancers.

Joint analysis of proteome, transcriptome, and multi-trait analysis to identify novel Parkinson's disease risk genes

Genome-wide association studies (GWAS) have identified multiple risk variants for Parkinson's disease (PD). Nevertheless, how the risk variants confer the risk of PD remains largely unknown. We conducted a proteome-wide association study (PWAS) and summary-data-based mendelian randomization (SMR) analysis by integrating PD GWAS with proteome and protein quantitative trait loci (pQTL) data from human brain, plasma and CSF. We also performed a large transcriptome-wide association study (TWAS) and Fine-mapping of causal gene sets (FOCUS), leveraging joint-tissue imputation (JTI) prediction models of 22 tissues to identify and prioritize putatively causal genes. We further conducted PWAS, SMR, TWAS, and FOCUS using a multi-trait analysis of GWAS (MTAG) to identify additional PD risk genes to boost statistical power. In this large-scale study, we identified 16 genes whose genetically regulated protein abundance levels were associated with Parkinson's disease risk. We undertook a large-scale analysis of PD and correlated traits, through TWAS and FOCUS studies, and discovered 26 casual genes related to PD that had not been reported in previous TWAS. 5 genes (CD38, GPNMB, RAB29, TMEM175, TTC19) showed significant associations with PD at both the proteome-wide and transcriptome-wide levels. Our study provides new insights into the etiology and underlying genetic architecture of PD.

omicSynth: An open multi-omic community resource for identifying druggable targets across neurodegenerative diseases

Treatments for neurodegenerative disorders remain rare, but recent FDA approvals, such as lecanemab and aducanumab for Alzheimer disease (MIM: 607822), highlight the importance of the underlying biological mechanisms in driving discovery and creating disease modifying therapies. The global population is aging, driving an urgent need for therapeutics that stop disease progression and eliminate symptoms. In this study, we create an open framework and resource for evidence-based identification of therapeutic targets for neurodegenerative disease. We use summary-data-based Mendelian randomization to identify genetic targets for drug discovery and repurposing. In parallel, we provide mechanistic insights into disease processes and potential network-level consequences of gene-based therapeutics. We identify 116 Alzheimer disease, 3 amyotrophic lateral sclerosis (MIM: 105400), 5 Lewy body dementia (MIM: 127750), 46 Parkinson disease (MIM: 605909), and 9 progressive supranuclear palsy (MIM: 601104) target genes passing multiple test corrections (pSMR_multi < 2.95 × 10-6 and pHEIDI > 0.01). We created a therapeutic scheme to classify our identified target genes into strata based on druggability and approved therapeutics, classifying 41 novel targets, 3 known targets, and 115 difficult targets (of these, 69.8% are expressed in the disease-relevant cell type from single-nucleus experiments). Our novel class of genes provides a springboard for new opportunities in drug discovery, development, and repurposing in the pre-competitive space. In addition, looking at drug-gene interaction networks, we identify previous trials that may require further follow-up such as riluzole in Alzheimer disease. We also provide a user-friendly web platform to help users explore potential therapeutic targets for neurodegenerative diseases, decreasing activation energy for the community.

Long-read RNA-seq atlas of novel microglia isoforms elucidates disease-associated genetic regulation of splicing

Microglia, the innate immune cells of the central nervous system, have been genetically implicated in multiple neurodegenerative diseases. We previously mapped the genetic regulation of gene expression and mRNA splicing in human microglia, identifying several loci where common genetic variants in microglia-specific regulatory elements explain disease risk loci identified by GWAS. However, identifying genetic effects on splicing has been challenging due to the use of short sequencing reads to identify causal isoforms. Here we present the isoform-centric microglia genomic atlas (isoMiGA) which leverages the power of long-read RNA-seq to identify 35,879 novel microglia isoforms. We show that the novel microglia isoforms are involved in stimulation response and brain region specificity. We then quantified the expression of both known and novel isoforms in a multi-ethnic meta-analysis of 555 human microglia short-read RNA-seq samples from 391 donors, the largest to date, and found associations with genetic risk loci in Alzheimer's disease and Parkinson's disease. We nominate several loci that may act through complex changes in isoform and splice site usage.

Epigenome-wide association study of peripheral immune cell populations in Parkinson's disease

Understanding the contribution of immune mechanisms to Parkinson's disease pathogenesis is an important challenge, potentially of major therapeutic implications. To further elucidate the involvement of peripheral immune cells, we studied epigenome-wide DNA methylation in isolated populations of CD14+ monocytes, CD19+ B cells, CD4+ T cells, and CD8+ T cells from Parkinson's disease patients and healthy control participants. We included 25 patients with a maximum five years of disease duration and 25 controls, and isolated four immune cell populations from each fresh blood sample. Epigenome-wide DNA methylation profiles were generated from 186 samples using the Illumina MethylationEpic array and association with disease status was tested using linear regression models. We identified six differentially methylated CpGs in CD14+ monocytes and one in CD8 + T cells. Four differentially methylated regions were identified in monocytes, including a region upstream of RAB32, a gene that has been linked to LRRK2. Methylation upstream of RAB32 correlated negatively with mRNA expression, and RAB32 expression was upregulated in Parkinson's disease both in our samples and in summary statistics from a previous study. Our epigenome-wide association study of early Parkinson's disease provides evidence for methylation changes across different peripheral immune cell types, highlighting monocytes and the RAB32 locus. The findings were predominantly cell-type-specific, demonstrating the value of isolating purified cell populations for genomic studies.

Expanding causal genes for Parkinson's disease via multi-omics analysis

Genome‑wide association studies (GWASs) have revealed numerous loci associated with Parkinson's disease (PD). However, some potential causal/risk genes were still not revealed and no etiological therapies are available. To find potential causal genes and explore genetically supported drug targets for PD is urgent. By integrating the expression quantitative trait loci (eQTL) and protein quantitative trait loci (pQTL) datasets from multiple tissues (blood, cerebrospinal fluid (CSF) and brain) and PD GWAS summary statistics, a pipeline combing Mendelian randomization (MR), Steiger filtering analysis, Bayesian colocalization, fine mapping, Protein-protein network and enrichment analysis were applied to identify potential causal genes for PD. As a result, GPNMB displayed a robust causal role for PD at the protein level in the blood, CSF and brain, and transcriptional level in the brain, while the protective role of CD38 (in brain pQTL and eQTL) was also identified. We also found inconsistent roles of DGKQ on PD between protein and mRNA levels. Another 9 proteins (CTSB, ARSA, SEC23IP, CD84, ENTPD1, FCGR2B, BAG3, SNCA, FCGR2A) were associated with the risk for PD based on only a single pQTL after multiple corrections. We also identified some proteins' interactions with known PD causative genes and therapeutic targets. In conclusion, this study suggested GPNMB, CD38, and DGKQ may act in the pathogenesis of PD, but whether the other proteins involved in PD needs more evidence. These findings would help to uncover the genes underlying PD and prioritize targets for future therapeutic interventions.

The human microglia responsome: a resource to better understand microglia states in health and disease

Microglia, the immune cells of the brain, are increasingly implicated in neurodegenerative disorders through genetic studies. However, how genetic risk factors for these diseases are related to microglial gene expression, microglial function, and ultimately disease, is still largely unknown. Microglia change rapidly in response to alterations in their cellular environment, which is regulated through changes in transcriptional programs, which are as yet poorly understood. Here, we compared the effects of a set of inflammatory and restorative stimuli (lipopolysaccharide, interferon-gamma, resiquimod, tumor necrosis factor-alpha, adenosine triphosphate, dexamethasone, and interleukin-4) on human microglial cells from 67 different donors (N = 398 samples) at the gene and transcript level. We show that microglia from different anatomical brain regions show distinct responses to inflammatory stimuli. We observed a greater overlap between human stimulated microglia and human monocytes than with mouse microglia. We define specific microglial signatures across conditions which are highly relevant for a wide range of biological functions and complex human diseases. Finally, we used our stimulation signatures to interpret associations from Alzheimer's disease (AD) genetic studies and microglia by integrating our inflammatory gene expression profiles with common genetic variants to map cis -expression QTLs (eQTLs). Together, we provide the most comprehensive transcriptomic database of the human microglia responsome.

Highlights

RNA-sequencing of 398 human microglial samples exposed to six different triggers.Microglia from different anatomical regions show distinct stimulation responses.Responses in human microglia show a greater overlap with human monocytes than murine microglia.Mapping of response Quantitative Trait Loci identifies interactions between genotype and effect of stimulation on gene expression.Our atlas provides a reference map for interpreting microglia signatures in health and disease.

The systemic immune response in Parkinson's disease: focus on the peripheral immune component

During Parkinson's disease (PD), both the central nervous system (CNS) and peripheral nervous system (PNS) are affected. In parallel, innate immune cells respond early to neuronal changes and alpha-synuclein (α-syn) pathology. Moreover, some of the affected neuronal groups innervate organs with a relevant role in immunity. Consequently, not only microglia, but also peripheral immune cells are altered, resulting in a systemic immune response. Innate and adaptive immune cells may participate in the neurodegenerative process by acting peripherally, infiltrating the brain, or releasing mediators that can protect or harm neurons. However, the sequence of the changes and the significance of each immune compartment in the disease remain to be clarified. In this review, we describe current understanding of the peripheral immune response in PD and discuss the road ahead.

A Transcriptome-Wide Analysis of Psoriasis: Identifying the Potential Causal Genes and Drug Candidates

Psoriasis is a chronic inflammatory skin disease characterized by cutaneous eruptions and pruritus. Because the genetic backgrounds of psoriasis are only partially revealed, an integrative and rigorous study is necessary. We conducted a transcriptome-wide association study (TWAS) with the new Genotype-Tissue Expression version 8 reference panels, including some tissue and multi-tissue panels that were not used previously. We performed tissue-specific heritability analyses on genome-wide association study data to prioritize the tissue panels for TWAS analysis. TWAS and colocalization (COLOC) analyses were performed with eight tissues from the single-tissue panels and the multi-tissue panels of context-specific genetics (CONTENT) to increase tissue specificity and statistical power. From TWAS, we identified the significant associations of 101 genes in the single-tissue panels and 64 genes in the multi-tissue panels, of which 26 genes were replicated in the COLOC. Functional annotation and network analyses identified that the genes were associated with psoriasis and/or immune responses. We also suggested drug candidates that interact with jointly significant genes through a conditional and joint analysis. Together, our findings may contribute to revealing the underlying genetic mechanisms and provide new insights into treatments for psoriasis.

omicSynth: an Open Multi-omic Community Resource for Identifying Druggable Targets across Neurodegenerative Diseases

Treatments for neurodegenerative disorders remain rare, although recent FDA approvals, such as Lecanemab and Aducanumab for Alzheimer's Disease, highlight the importance of the underlying biological mechanisms in driving discovery and creating disease modifying therapies. The global population is aging, driving an urgent need for therapeutics that stop disease progression and eliminate symptoms. In this study, we create an open framework and resource for evidence-based identification of therapeutic targets for neurodegenerative disease. We use Summary-data-based Mendelian Randomization to identify genetic targets for drug discovery and repurposing. In parallel, we provide mechanistic insights into disease processes and potential network-level consequences of gene-based therapeutics. We identify 116 Alzheimer's disease, 3 amyotrophic lateral sclerosis, 5 Lewy body dementia, 46 Parkinson's disease, and 9 Progressive supranuclear palsy target genes passing multiple test corrections (pSMR_multi < 2.95×10-6 and pHEIDI > 0.01). We created a therapeutic scheme to classify our identified target genes into strata based on druggability and approved therapeutics - classifying 41 novel targets, 3 known targets, and 115 difficult targets (of these 69.8% are expressed in the disease relevant cell type from single nucleus experiments). Our novel class of genes provides a springboard for new opportunities in drug discovery, development and repurposing in the pre-competitive space. In addition, looking at drug-gene interaction networks, we identify previous trials that may require further follow-up such as Riluzole in AD. We also provide a user-friendly web platform to help users explore potential therapeutic targets for neurodegenerative diseases, decreasing activation energy for the community [https://nih-card-ndd-smr-home-syboky.streamlit.app/].

LRRK2 and Parkinson's disease: from genetics to targeted therapy

LRRK2 variants are implicated in both familial and sporadic PD. LRRK2-PD has a generally benign clinical presentation and variable pathology, with inconsistent presence of Lewy bodies and marked Alzheimer's disease pathology. The mechanisms underlying LRRK2-PD are still unclear, but inflammation, vesicle trafficking, lysosomal homeostasis, and ciliogenesis have been suggested, among others. As novel therapies targeting LRRK2 are under development, understanding the role and function of LRRK2 in PD is becoming increasingly important. Here, we outline the epidemiological, pathophysiological, and clinical features of LRRK2-PD, and discuss the arising therapeutic approaches targeting LRRK2 and possible future directions for research.

Phosphatidylinositol-3,4,5-trisphosphate interacts with alpha-synuclein and initiates its aggregation and formation of Parkinson's disease-related fibril polymorphism

Lipid interaction with α-synuclein (αSyn) has been long implicated in the pathogenesis of Parkinson's disease (PD). However, it has not been fully determined which lipids are involved in the initiation of αSyn aggregation in PD. Here exploiting genetic understanding associating the loss-of-function mutation in Synaptojanin 1 (SYNJ1), a phosphoinositide phosphatase, with familial PD and analysis of postmortem PD brains, we identified a novel lipid molecule involved in the toxic conversion of αSyn and its relation to PD. We first established a SYNJ1 knockout cell model and found SYNJ1 depletion increases the accumulation of pathological αSyn. Lipidomic analysis revealed SYNJ1 depletion elevates the level of its substrate phosphatidylinositol-3,4,5-trisphosphate (PIP₃). We then employed Caenorhabditis elegans model to examine the effect of SYNJ1 defect on the neurotoxicity of αSyn. Mutations in SYNJ1 accelerated the accumulation of αSyn aggregation and induced locomotory defects in the nematodes. These results indicate that functional loss of SYNJ1 promotes the pathological aggregation of αSyn via the dysregulation of its substrate PIP₃, leading to the aggravation of αSyn-mediated neurodegeneration. Treatment of cultured cell line and primary neurons with PIP₃ itself or with PIP₃ phosphatase inhibitor resulted in intracellular formation of αSyn inclusions. Indeed, in vitro protein-lipid overlay assay validated that phosphoinositides, especially PIP₃, strongly interact with αSyn. Furthermore, the aggregation assay revealed that PIP₃ not only accelerates the fibrillation of αSyn, but also induces the formation of fibrils sharing conformational and biochemical characteristics similar to the fibrils amplified from the brains of PD patients. Notably, the immunohistochemical and lipidomic analyses on postmortem brain of patients with sporadic PD showed increased PIP₃ level and its colocalization with αSyn. Taken together, PIP₃ dysregulation promotes the pathological aggregation of αSyn and increases the risk of developing PD, and PIP₃ represents a potent target for intervention in PD.

Neurodegeneration cell per cell

The clinical definition of neurodegenerative diseases is based on symptoms that reflect terminal damage of specific brain regions. This is misleading as it tells little about the initial disease processes. Circuitry failures that underlie the clinical symptomatology are themselves preceded by clinically mostly silent, slowly progressing multicellular processes that trigger or are triggered by the accumulation of abnormally folded proteins such as Aβ, Tau, TDP-43, and α-synuclein, among others. Methodological advances in single-cell omics, combined with complex genetics and novel ways to model complex cellular interactions using induced pluripotent stem (iPS) cells, make it possible to analyze the early cellular phase of neurodegenerative disorders. This will revolutionize the way we study those diseases and will translate into novel diagnostics and cell-specific therapeutic targets, stopping these disorders in their early track before they cause difficult-to-reverse damage to the brain.

Genetic modifiers of synucleinopathies-lessons from experimental models

α-Synuclein is a pleiotropic protein underlying a group of progressive neurodegenerative diseases, including Parkinson's disease and dementia with Lewy bodies. Together, these are known as synucleinopathies. Like all neurological diseases, understanding of disease mechanisms is hampered by the lack of access to biopsy tissues, precluding a real-time view of disease progression in the human body. This has driven researchers to devise various experimental models ranging from yeast to flies to human brain organoids, aiming to recapitulate aspects of synucleinopathies. Studies of these models have uncovered numerous genetic modifiers of α-synuclein, most of which are evolutionarily conserved. This review discusses what we have learned about disease mechanisms from these modifiers, and ways in which the study of modifiers have supported ongoing efforts to engineer disease-modifying interventions for synucleinopathies.

Genome- and transcriptome-wide splicing associations with alcohol use disorder

Genetic mechanisms of alternative mRNA splicing have been shown in the brain for a variety of neuropsychiatric traits, but not substance use disorders. Our study utilized RNA-sequencing data on alcohol use disorder (AUD) in four brain regions (n = 56; ages 40-73; 100% 'Caucasian'; PFC, NAc, BLA and CEA) and genome-wide association data on AUD (n = 435,563, ages 22-90; 100% European-American). Polygenic scores of AUD were associated with AUD-related alternative mRNA splicing in the brain. We identified 714 differentially spliced genes between AUD vs controls, which included both putative addiction genes and novel gene targets. We found 6463 splicing quantitative trait loci (sQTLs) that linked to the AUD differentially spliced genes. sQTLs were enriched in loose chromatin genomic regions and downstream gene targets. Additionally, the heritability of AUD was enriched for DNA variants in and around differentially spliced genes associated with AUD. Our study also performed splicing transcriptome-wide association studies (TWASs) of AUD and other drug use traits that unveiled specific genes for follow-up and splicing correlations across SUDs. Finally, we showed that differentially spliced genes between AUD vs control were also associated with primate models of chronic alcohol consumption in similar brain regions. Our study found substantial genetic contributions of alternative mRNA splicing in AUD.

Post-translational proteomics platform identifies neurite outgrowth impairments in Parkinson's disease GBA-N370S dopamine neurons

Variants at the GBA locus, encoding glucocerebrosidase, are the strongest common genetic risk factor for Parkinson's disease (PD). To understand GBA-related disease mechanisms, we use a multi-part-enrichment proteomics and post-translational modification (PTM) workflow, identifying large numbers of dysregulated proteins and PTMs in heterozygous GBA-N370S PD patient induced pluripotent stem cell (iPSC) dopamine neurons. Alterations in glycosylation status show disturbances in the autophagy-lysosomal pathway, which concur with upstream perturbations in mammalian target of rapamycin (mTOR) activation in GBA-PD neurons. Several native and modified proteins encoded by PD-associated genes are dysregulated in GBA-PD neurons. Integrated pathway analysis reveals impaired neuritogenesis in GBA-PD neurons and identify tau as a key pathway mediator. Functional assays confirm neurite outgrowth deficits and identify impaired mitochondrial movement in GBA-PD neurons. Furthermore, pharmacological rescue of glucocerebrosidase activity in GBA-PD neurons improves the neurite outgrowth deficit. Overall, this study demonstrates the potential of PTMomics to elucidate neurodegeneration-associated pathways and potential drug targets in complex disease models.

LRRK2: Genetic mechanisms vs genetic subtypes

In 2004, the identification of pathogenic variants in the LRRK2 gene across several families with autosomal dominant late-onset Parkinson's disease (PD) revolutionized our understanding of the role of genetics in PD. Previous beliefs that genetics in PD was limited to rare early-onset or familial forms of the disease were quickly dispelled. Currently, we recognize LRRK2 p.G2019S as the most common genetic cause of both sporadic and familial PD, with more than 100,000 affected carriers across the globe. The frequency of LRRK2 p.G2019S is also highly variable across populations, with some regions of Asian or Latin America reporting close to 0%, contrasting to Ashkenazi Jews or North African Berbers reporting up to 13% and 40%, respectively. Patients with LRRK2 pathogenic variants are clinically and pathologically heterogeneous, highlighting the age-related variable penetrance that also characterizes LRRK2-related disease. Indeed, the majority of patients with LRRK2-related disease are characterized by a relatively mild Parkinsonism with less motor symptoms with variable presence of α-synuclein and/or tau aggregates, with pathologic pleomorphism widely described. At a functional cellular level, it is likely that pathogenic variants mediate a toxic gain-of-function of the LRRK2 protein resulting in increased kinase activity perhaps in a cell-specific manner; by contrast, some LRRK2 variants appear to be protective reducing PD risk by decreasing the kinase activity. Therefore, employing this information to define appropriate patient populations for clinical trials of targeted kinase LRRK2 inhibition strategies is very promising and demonstrates a potential future application for PD using precision medicine.

Going broad and deep: sequencing-driven insights into plant physiology, evolution, and crop domestication

Deep sequencing is a term that has become embedded in the plant genomic literature in recent years and with good reason. A torrent of (largely) high-quality genomic and transcriptomic data has been collected and most of this has been publicly released. Indeed, almost 1000 plant genomes have been reported (www.plabipd.de) and the 2000 Plant Transcriptomes Project has long been completed. The EarthBioGenome project will dwarf even these milestones. That said, massive progress in understanding plant physiology, evolution, and crop domestication has been made by sequencing broadly (across a species) as well as deeply (within a single individual). We will outline the current state of the art in genome and transcriptome sequencing before we briefly review the most visible of these broad approaches, namely genome-wide association and transcriptome-wide association studies, as well as the compilation of pangenomes. This will include both (i) the most commonly used methods reliant on single nucleotide polymorphisms and short InDels and (ii) more recent examples which consider structural variants. We will subsequently present case studies exemplifying how their application has brought insight into either plant physiology or evolution and crop domestication. Finally, we will provide conclusions and an outlook as to the perspective for the extension of such approaches to different species, tissues, and biological processes.

Integrative genetic analysis illuminates ALS heritability and identifies risk genes

Amyotrophic lateral sclerosis (ALS) has substantial heritability, in part shared with fronto-temporal dementia (FTD). We show that ALS heritability is enriched in splicing variants and in binding sites of 6 RNA-binding proteins including TDP-43 and FUS. A transcriptome wide association study (TWAS) identified 6 loci associated with ALS, including in NUP50 encoding for the nucleopore basket protein NUP50. Independently, rare variants in NUP50 were associated with ALS risk (P = 3.71.10-03; odds ratio = 3.29; 95%CI, 1.37 to 7.87) in a cohort of 9,390 ALS/FTD patients and 4,594 controls. Cells from one patient carrying a NUP50 frameshift mutation displayed a decreased level of NUP50. Loss of NUP50 leads to death of cultured neurons, and motor defects in Drosophila and zebrafish. Thus, our study identifies alterations in splicing in neurons as critical in ALS and provides genetic evidence linking nuclear pore defects to ALS.

Integrative transcriptomic analysis of the amyotrophic lateral sclerosis spinal cord implicates glial activation and suggests new risk genes

Amyotrophic lateral sclerosis (ALS) is a progressively fatal neurodegenerative disease affecting motor neurons in the brain and spinal cord. In this study, we investigated gene expression changes in ALS via RNA sequencing in 380 postmortem samples from cervical, thoracic and lumbar spinal cord segments from 154 individuals with ALS and 49 control individuals. We observed an increase in microglia and astrocyte gene expression, accompanied by a decrease in oligodendrocyte gene expression. By creating a gene co-expression network in the ALS samples, we identified several activated microglia modules that negatively correlate with retrospective disease duration. We mapped molecular quantitative trait loci and found several potential ALS risk loci that may act through gene expression or splicing in the spinal cord and assign putative cell types for FNBP1, ACSL5, SH3RF1 and NFASC. Finally, we outline how common genetic variants associated with splicing of C9orf72 act as proxies for the well-known repeat expansion, and we use the same mechanism to suggest ATXN3 as a putative risk gene.

MED13 and glycolysis are conserved modifiers of α-synuclein-associated neurodegeneration

α-Synuclein (α-syn) is important in synucleinopathies such as Parkinson's disease (PD). While genome-wide association studies (GWASs) of synucleinopathies have identified many risk loci, the underlying genes have not been shown for most loci. Using Drosophila, we screened 3,471 mutant chromosomes for genetic modifiers of α-synuclein and identified 12 genes. Eleven modifiers have human orthologs associated with diseases, including MED13 and CDC27, which lie within PD GWAS loci. Drosophila Skd/Med13 and glycolytic enzymes are co-upregulated by α-syn-associated neurodegeneration. While elevated α-syn compromises mitochondrial function, co-expressing skd/Med13 RNAi and α-syn synergistically increase the ratio of oxidized-to-reduced glutathione. The resulting neurodegeneration can be suppressed by overexpressing a glycolytic enzyme or treatment with deferoxamine, suggesting that compensatory glycolysis is neuroprotective. In addition, the functional relationship between α-synuclein, MED13, and glycolytic enzymes is conserved between flies and mice. We propose that hypoxia-inducible factor and MED13 are part of a druggable pathway for PD.

Mendelian Randomization Study Using Dopaminergic Neuron-Specific eQTL Nominates Potential Causal Genes for Parkinson's Disease

BACKGROUND

Large-scale genome-wide association studies (GWASs) have reported multiple risk variants for Parkinson's disease (PD). However, little is known about how these reported risk variants confer risk of PD.

OBJECTIVE

To nominate genes whose genetically regulated expression in dopaminergic neurons may have a causal role in PD.

METHODS

We conducted a two-sample Mendelian randomization (MR) study by integrating large-scale genome-wide associations and expression quantitative trait loci (eQTL) data from dopaminergic neurons.

RESULTS

MR analysis nominated 10 risk genes whose genetically regulated expression in dopaminergic neurons may have a causal role in PD. These MR significant genes include FAM200B, NDUFAF2, NUP42, SH3GL2, STX1B, CCDC189, KAT8, PRSS36, VAMP4, and ZSWIM7.

CONCLUSIONS

We report the first MR study of PD by using dopaminergic neuron-specific eQTL and nominate novel risk genes for PD. Further functional characterization of the nominated risk genes will provide mechanistic insights into PD pathogenesis and potential therapeutic targets. © 2022 International Parkinson and Movement Disorder Society.

Bridging the splicing gap in human genetics with long-read RNA sequencing: finding the protein isoform drivers of disease

Aberrant splicing underlies many human diseases, including cancer, cardiovascular diseases and neurological disorders. Genome-wide mapping of splicing quantitative trait loci (sQTLs) has shown that genetic regulation of alternative splicing is widespread. However, identification of the corresponding isoform or protein products associated with disease-associated sQTLs is challenging with short-read RNA-seq, which cannot precisely characterize full-length transcript isoforms. Furthermore, contemporary sQTL interpretation often relies on reference transcript annotations, which are incomplete. Solutions to these issues may be found through integration of newly emerging long-read sequencing technologies. Long-read sequencing offers the capability to sequence full-length mRNA transcripts and, in some cases, to link sQTLs to transcript isoforms containing disease-relevant protein alterations. Here, we provide an overview of sQTL mapping approaches, the use of long-read sequencing to characterize sQTL effects on isoforms, the linkage of RNA isoforms to protein-level functions and comment on future directions in the field. Based on recent progress, long-read RNA sequencing promises to be part of the human disease genetics toolkit to discover and treat protein isoforms causing rare and complex diseases.

A transcriptome-wide association study of uterine fibroids to identify potential genetic markers and toxic chemicals

Uterine fibroid is one of the most prevalent benign tumors in women, with high socioeconomic costs. Although genome-wide association studies (GWAS) have identified several loci associated with uterine fibroid risks, they could not successfully interpret the biological effects of genomic variants at the gene expression levels. To prioritize uterine fibroid susceptibility genes that are biologically interpretable, we conducted a transcriptome-wide association study (TWAS) by integrating GWAS data of uterine fibroid and expression quantitative loci data. We identified nine significant TWAS genes including two novel genes, RP11-282O18.3 and KBTBD7, which may be causal genes for uterine fibroid. We conducted functional enrichment network analyses using the TWAS results to investigate the biological pathways in which the overall TWAS genes were involved. The results demonstrated the immune system process to be a key pathway in uterine fibroid pathogenesis. Finally, we carried out chemical–gene interaction analyses using the TWAS results and the comparative toxicogenomics database to determine the potential risk chemicals for uterine fibroid. We identified five toxic chemicals that were significantly associated with uterine fibroid TWAS genes, suggesting that they may be implicated in the pathogenesis of uterine fibroid. In this study, we performed an integrative analysis covering the broad application of bioinformatics approaches. Our study may provide a deeper understanding of uterine fibroid etiologies and informative notifications about potential risk chemicals for uterine fibroid.

Immunogenetic Determinants of Parkinson’s Disease Etiology

Parkinson’s disease (PD) is increasingly recognised as a systemic disorder in which inflammation might play a causative role rather than being a consequence or an epiphenomenon of the neurodegenerative process. Although growing genetic evidence links the central and peripheral immune system with both monogenic and sporadic PD, our understanding on how the immune system contributes to PD pathogenesis remains a daunting challenge. In this review, we discuss recent literature aimed at exploring the role of known genes and susceptibility loci to PD pathogenesis through immune system related mechanisms. Furthermore, we outline shared genetic etiologies and interrelations between PD and autoimmune diseases and underlining challenges and limitations faced in the translation of relevant allelic and regulatory risk loci to immune-pathological mechanisms. Lastly, with the field of immunogenetics expanding rapidly, we place these insights into a future context highlighting the prospect of immune modulation as a promising disease-modifying strategy.

DNA methylation signatures of Alzheimer's disease neuropathology in the cortex are primarily driven by variation in non-neuronal cell-types

Alzheimer's disease (AD) is a chronic neurodegenerative disease characterized by the progressive accumulation of amyloid-beta and neurofibrillary tangles of tau in the neocortex. We profiled DNA methylation in two regions of the cortex from 631 donors, performing an epigenome-wide association study of multiple measures of AD neuropathology. We meta-analyzed our results with those from previous studies of DNA methylation in AD cortex (total n = 2013 donors), identifying 334 cortical differentially methylated positions (DMPs) associated with AD pathology including methylomic variation at loci not previously implicated in dementia. We subsequently profiled DNA methylation in NeuN+ (neuronal-enriched), SOX10+ (oligodendrocyte-enriched) and NeuN-/SOX10- (microglia- and astrocyte-enriched) nuclei, finding that the majority of DMPs identified in 'bulk' cortex tissue reflect DNA methylation differences occurring in non-neuronal cells. Our study highlights the power of utilizing multiple measures of neuropathology to identify epigenetic signatures of AD and the importance of characterizing disease-associated variation in purified cell-types.

Genetic control of RNA splicing and its distinct role in complex trait variation

Most genetic variants identified from genome-wide association studies (GWAS) in humans are noncoding, indicating their role in gene regulation. Previous studies have shown considerable links of GWAS signals to expression quantitative trait loci (eQTLs) but the links to other genetic regulatory mechanisms, such as splicing QTLs (sQTLs), are underexplored. Here, we introduce an sQTL mapping method, testing for heterogeneity between isoform-eQTLeffects (THISTLE), with improved power over competing methods. Applying THISTLE together with a complementary sQTL mapping strategy to brain transcriptomic (n = 2,865) and genotype data, we identified 12,794 genes with cis-sQTLs at P < 5 × 10⁻⁸, approximately 61% of which were distinct from eQTLs. Integrating the sQTL data into GWAS for 12 brain-related complex traits (including diseases), we identified 244 genes associated with the traits through cis-sQTLs, approximately 61% of which could not be discovered using the corresponding eQTL data. Our study demonstrates the distinct role of most sQTLs in the genetic regulation of transcription and complex trait variation.

A powerful method for splicing quantitative trait loci (sQTL) mapping, THISTLE, is presented and applied to a collection of 2,865 brain samples. Integration with GWAS identifies 244 genes associated via cis-sQTLs, of which 61% were not identified using expression QTLs.

The Protein-Protein Interaction Network of Hereditary Parkinsonism Genes Is a Hierarchical Scale-Free Network

PURPOSE

Hereditary parkinsonism genes consist of causative genes of familial Parkinson's disease (PD) with a locus symbol prefix (PARK genes) and hereditary atypical parkinsonian disorders that present atypical features and limited responsiveness to levodopa (non-PARK genes). Although studies have shown that hereditary parkinsonism genes are related to idiopathic PD at the phenotypic, gene expression, and genomic levels, no study has systematically investigated connectivity among the proteins encoded by these genes at the protein-protein interaction (PPI) level.

MATERIALS AND METHODS

Topological measurements and physical interaction enrichment were performed to assess PPI networks constructed using some or all the proteins encoded by hereditary parkinsonism genes (n=96), which were curated using the Online Mendelian Inheritance in Man database and literature.

RESULTS

Non-PARK and PARK genes were involved in common functional modules related to autophagy, mitochondrial or lysosomal organization, catecholamine metabolic process, chemical synapse transmission, response to oxidative stress, neuronal apoptosis, regulation of cellular protein catabolic process, and vesicle-mediated transport in synapse. The hereditary parkinsonism proteins formed a single large network comprising 51 nodes, 83 edges, and three PPI pairs. The probability of degree distribution followed a power-law scaling behavior, with a degree exponent of 1.24 and a correlation coefficient of 0.92. LRRK2 was identified as a hub gene with the highest degree of betweenness centrality; its physical interaction enrichment score was 1.28, which was highly significant.

CONCLUSION

Both PARK and non-PARK genes show high connectivity at the PPI and biological functional levels.

Integrative transcriptome-wide analysis of atopic dermatitis for drug repositioning

Atopic dermatitis (AD) is one of the most common inflammatory skin diseases, which significantly impact the quality of life. Transcriptome-wide association study (TWAS) was conducted to estimate both transcriptomic and genomic features of AD and detected significant associations between 31 expression quantitative loci and 25 genes. Our results replicated well-known genetic markers for AD, as well as 4 novel associated genes. Next, transcriptome meta-analysis was conducted with 5 studies retrieved from public databases and identified 5 additional novel susceptibility genes for AD. Applying the connectivity map to the results from TWAS and meta-analysis, robustly enriched perturbations were identified and their chemical or functional properties were analyzed. Here, we report the first research on integrative approaches for an AD, combining TWAS and transcriptome meta-analysis. Together, our findings could provide a comprehensive understanding of the pathophysiologic mechanisms of AD and suggest potential drug candidates as alternative treatment options.

Integrative genomic and transcriptomic analyses on publicly available data-sets together with in silico drug repositioning identifies alternative therapeutic options to treat atopic dermatitis.

Opioid Use Disorder and Alternative mRNA Splicing in Reward Circuitry

Opiate/opioid use disorder (OUD) is a chronic relapsing brain disorder that has increased in prevalence in the last two decades in the United States. Understanding the molecular correlates of OUD may provide key insights into the pathophysiology of this syndrome. Using publicly available RNA-sequencing data, our study investigated the possible role of alternative mRNA splicing in human brain tissue (dorsal-lateral prefrontal cortex (dlPFC), nucleus accumbens (NAc), and midbrain) of 90 individuals with OUD or matched controls. We found a total of 788 differentially spliced genes across brain regions. Alternative mRNA splicing demonstrated mostly tissue-specific effects, but a functionally characterized splicing change in the clathrin and AP-2-binding (CLAP) domain of the Bridging Integrator 1 (BIN1) gene was significantly linked to OUD across all brain regions. We investigated two hypotheses that may underlie differential splicing in OUD. First, we tested whether spliceosome genes were disrupted in the brains of individuals with OUD. Pathway enrichment analyses indicated spliceosome perturbations in OUD across brain regions. Second, we tested whether alternative mRNA splicing regions were linked to genetic predisposition. Using a genome-wide association study (GWAS) of OUD, we found no evidence that DNA variants within or surrounding differentially spliced genes were implicated in the heritability of OUD. Altogether, our study contributes to the understanding of OUD pathophysiology by providing evidence of a possible role of alternative mRNA splicing in OUD.

Molecular Quantitative Trait Locus Mapping in Human Complex Diseases

Mapping quantitative trait loci (QTLs) for molecular traits from chromatin to metabolites (i.e., xQTLs) provides insight into the locations and effect modes of genetic variants that influence these molecular phenotypes and the propagation of functional consequences of each variant. xQTL studies indirectly interrogate the functional landscape of the molecular basis of complex diseases, including the impact of non-coding regulatory variants, the tissue specificity of regulatory elements, and their contribution to disease by integrating with genome-wide association studies (GWAS). We summarize a variety of molecular xQTL studies in human tissues and cells. In addition, using the Alzheimer's Disease Sequencing Project (ADSP) as an example, we describe the ADSP xQTL project, a collaborative effort across the ADSP Functional Genomics Consortium (ADSP-FGC). The project's ultimate goal is a reference map of Alzheimer's-related QTLs using existing datasets from multiple omics layers to help us study the consequences of genetic variants identified in the ADSP. xQTL studies enable the identification of the causal genes and pathways in GWAS loci, which will likely aid in the discovery of novel biomarkers and therapeutic targets for complex diseases. © 2022 Wiley Periodicals LLC.

Fine-mapping of Parkinson's disease susceptibility loci identifies putative causal variants

Recent genome-wide association studies have identified 78 loci associated with Parkinson's disease susceptibility but the underlying mechanisms remain largely unclear. To identify likely causal variants for disease risk, we fine-mapped these Parkinson's-associated loci using four different fine-mapping methods. We then integrated multi-assay cell type-specific epigenomic profiles to pinpoint the likely mechanism of action of each variant, allowing us to identify Consensus single nucleotide polymorphism (SNPs) that disrupt LRRK2 and FCGR2A regulatory elements in microglia, an MBNL2 enhancer in oligodendrocytes, and a DYRK1A enhancer in neurons. This genome-wide functional fine-mapping investigation of Parkinson's disease substantially advances our understanding of the causal mechanisms underlying this complex disease while avoiding focus on spurious, non-causal mechanisms. Together, these results provide a robust, comprehensive list of the likely causal variants, genes and cell-types underlying Parkinson's disease risk as demonstrated by consistently greater enrichment of our fine-mapped SNPs relative to lead GWAS SNPs across independent functional impact annotations. In addition, our approach prioritized an average of 3/85 variants per locus as putatively causal, making downstream experimental studies both more tractable and more likely to yield disease-relevant, actionable results. Large-scale studies comparing individuals with Parkinson's disease to age-matched controls have identified many regions of the genome associated with the disease. However, there is widespread correlation between different parts of the genome, making it difficult to tell which genetic variants cause Parkinson's and which are simply co-inherited with causal variants. We therefore applied a suite of statistical models to identify the most likely causal genetic variants (i.e. fine-mapping). We then linked these genetic variants with epigenomic and gene expression signatures across a wide variety of tissues and cell types to identify how these variants cause disease. Therefore, this study provides a comprehensive and robust list of cellular and molecular mechanisms that may serve as targets in the development of more effective Parkinson's therapeutics.

HyMM: hybrid method for disease-gene prediction by integrating multiscale module structure

MOTIVATION

Identifying disease-related genes is an important issue in computational biology. Module structure widely exists in biomolecule networks, and complex diseases are usually thought to be caused by perturbations of local neighborhoods in the networks, which can provide useful insights for the study of disease-related genes. However, the mining and effective utilization of the module structure is still challenging in such issues as a disease gene prediction.

RESULTS

We propose a hybrid disease-gene prediction method integrating multiscale module structure (HyMM), which can utilize multiscale information from local to global structure to more effectively predict disease-related genes. HyMM extracts module partitions from local to global scales by multiscale modularity optimization with exponential sampling, and estimates the disease relatedness of genes in partitions by the abundance of disease-related genes within modules. Then, a probabilistic model for integration of gene rankings is designed in order to integrate multiple predictions derived from multiscale module partitions and network propagation, and a parameter estimation strategy based on functional information is proposed to further enhance HyMM's predictive power. By a series of experiments, we reveal the importance of module partitions at different scales, and verify the stable and good performance of HyMM compared with eight other state-of-the-arts and its further performance improvement derived from the parameter estimation.

CONCLUSIONS

The results confirm that HyMM is an effective framework for integrating multiscale module structure to enhance the ability to predict disease-related genes, which may provide useful insights for the study of the multiscale module structure and its application in such issues as a disease-gene prediction.

N-Acetyl-Cysteine: Modulating the Cysteine Redox Proteome in Neurodegenerative Diseases

In the last twenty years, significant progress in understanding the pathophysiology of age-associated neurodegenerative diseases has been made. However, the prevention and treatment of these diseases remain without clinically significant therapeutic advancement. While we still hope for some potential genetic therapeutic approaches, the current reality is far from substantial progress. With this state of the issue, emphasis should be placed on early diagnosis and prompt intervention in patients with increased risk of neurodegenerative diseases to slow down their progression, poor prognosis, and decreasing quality of life. Accordingly, it is urgent to implement interventions addressing the psychosocial and biochemical disturbances we know are central in managing the evolution of these disorders. Genomic and proteomic studies have shown the high molecular intricacy in neurodegenerative diseases, involving a broad spectrum of cellular pathways underlying disease progression. Recent investigations indicate that the dysregulation of the sensitive-cysteine proteome may be a concurrent pathogenic mechanism contributing to the pathophysiology of major neurodegenerative diseases, opening new therapeutic opportunities. Considering the incidence and prevalence of these disorders and their already significant burden in Western societies, they will become a real pandemic in the following decades. Therefore, we propose large-scale investigations, in selected groups of people over 40 years of age with decreased blood glutathione levels, comorbidities, and/or mild cognitive impairment, to evaluate supplementation of the diet with low doses of N-acetyl-cysteine, a promising and well-tolerated therapeutic agent suitable for long-term use.

LRRK2 mutant knock-in mouse models: therapeutic relevance in Parkinson's disease

Mutations in the leucine-rich repeat kinase 2 gene (LRRK2) are one of the most frequent genetic causes of both familial and sporadic Parkinson's disease (PD). Mounting evidence has demonstrated pathological similarities between LRRK2-associated PD (LRRK2-PD) and sporadic PD, suggesting that LRRK2 is a potential disease modulator and a therapeutic target in PD. LRRK2 mutant knock-in (KI) mouse models display subtle alterations in pathological aspects that mirror early-stage PD, including increased susceptibility of nigrostriatal neurotransmission, development of motor and non-motor symptoms, mitochondrial and autophagy-lysosomal defects and synucleinopathies. This review provides a rationale for the use of LRRK2 KI mice to investigate the LRRK2-mediated pathogenesis of PD and implications from current findings from different LRRK2 KI mouse models, and ultimately discusses the therapeutic potentials against LRRK2-associated pathologies in PD.