- Open Access
Unveiling sex-based differences in Parkinson's disease: a comprehensive meta-analysis of transcriptomic studies
Biology of Sex Differences volume 13, Article number: 68 (2022)
In recent decades, increasing longevity (among other factors) has fostered a rise in Parkinson's disease incidence. Although not exhaustively studied in this devastating disease, the impact of sex represents a critical variable in Parkinson’s disease as epidemiological and clinical features differ between males and females.
To study sex bias in Parkinson’s disease, we conducted a systematic review to select sex-labeled transcriptomic data from three relevant brain tissues: the frontal cortex, the striatum, and the substantia nigra. We performed differential expression analysis on each study chosen. Then we summarized the individual differential expression results with three tissue-specific meta-analyses and a global all-tissues meta-analysis. Finally, results from the meta-analysis were functionally characterized using different functional profiling approaches.
The tissue-specific meta-analyses linked Parkinson’s disease to the enhanced expression of MED31 in the female frontal cortex and the dysregulation of 237 genes in the substantia nigra. The global meta-analysis detected 15 genes with sex-differential patterns in Parkinson’s disease, which participate in mitochondrial function, oxidative stress, neuronal degeneration, and cell death. Furthermore, functional analyses identified pathways, protein–protein interaction networks, and transcription factors that differed by sex. While male patients exhibited changes in oxidative stress based on metal ions, inflammation, and angiogenesis, female patients exhibited dysfunctions in mitochondrial and lysosomal activity, antigen processing and presentation functions, and glutamic and purine metabolism. All results generated during this study are readily available by accessing an open web resource (http://bioinfo.cipf.es/metafun-pd/) for consultation and reuse in further studies.
Our in silico approach has highlighted sex-based differential mechanisms in typical Parkinson Disease hallmarks (inflammation, mitochondrial dysfunction, and oxidative stress). Additionally, we have identified specific genes and transcription factors for male and female Parkinson Disease patients that represent potential candidates as biomarkers to diagnosis.
Females show a significant increase in the expression of MED31 in the frontal cortex. This gene is involved in lipid metabolism and neural diseases.
We found 237 genes having sex-based significantly differential expression in substantia nigra.
Functional profiling of these genes reveals a differential sex-related behavior in PD regarding their biological functions, protein-protein interaction networks, and transcription factors activation.
There are remarkable sex based differential mechanisms in typical PD hallmarks: inflammation, mitochondrial dysfunction, and oxidative stress.
Studies on sex differences in PD are needed to improve more targeted interventions.
Parkinson's disease (PD) is the second most common progressive neurodegenerative disease after Alzheimer's disease and the most rapidly expanding related disorder in the elderly population, affecting 1% of individuals over 65 and 2% of individuals over 85 . Importantly, PD also occurs in younger patients, with a prevalence of 0.004% in individuals between the ages of 40 and 50 . PD affects different areas of the brain, with the frontal cortex (FC—involved in mental tasks), the striatum tissue (ST—coordinates multiple aspects of cognition including motor and action planning, decision-making, motivation, reinforcement, and reward perception), and the substantia nigra (SN—controls eye movement, motor planning, reward-seeking, learning, and addiction) the most affected. The selective degeneration of SN dopaminergic neurons that innervate the ST (nigrostriatal dopaminergic neurons—NSDA) and the formation of Lewy bodies (α-synuclein accumulation in neurons) in several different brain regions  prompt motor symptoms (e.g., resting tremor, muscle stiffness, bradykinesia, and postural instability) and non-motor symptoms such as cognitive and mood impairment. In advanced-stage PD, a third of patients develop dementia . Oxidative stress, excitotoxicity, and neuroinflammation influence neuron death in PD [5, 6], while nitric oxide and other reactive nitrogen species [7, 8] promote PD progression. Mitochondrial dysfunction also plays a significant role in PD, as in other neurodegenerative diseases . While the etiology of PD remains unclear, studies have provided evidence that genetic factors and environmental triggers, including genetic variation, sex, age, and exposure to environmental toxins (e.g., herbicides and pesticides), all play roles in pathology development .
Importantly, male and female PD patients display specific sex-based differences . While male PD patients suffer from a greater risk of PD , female PD patients suffer from higher mortality rates and require earlier documented placements in nursing care . Other sex-based differences in PD include nigrostriatal degeneration, time of symptom onset, motor and non-motor symptoms, REM sleep behavior disorder, treatment outcomes, and disease mechanisms [14,15,16,17,18]. Inflammation and mitochondrial function also display sex-based alterations in PD patients [19, 20]. Genetic factors, sex chromosome genes, hormones, and neuroactive steroids represent the leading causes for these sex-based differences [15, 18, 21].
While a minority of PD cases benefit from genetic testing, routine clinical methods for PD diagnosis at early disease stages do not exist, complicating disease management. There exist therapeutic options such as l-DOPA (levodopa or l-3,4-dihydroxyphenylalanine) and deep brain stimulation, which improve symptoms and patient quality-of-life; however, these treatments do not delay neurodegeneration. Thus, identifying PD biomarkers represents an unmet clinical task that may allow for early diagnosis and the creation of improved therapies.
To evaluate sex-based differences in PD, we now present a systematic review and 4 transcriptomic meta-analyses on different brain regions followed by a functional characterization. These results demonstrate significant differences in gene expression and biological functions when comparing male and female PD patients and brain regions.
All bioinformatics and statistical analysis were performed using R software v.3.6.3 . Additional file 2: Table S2 details R packages and versions.
Study search and selection
Available datasets were collected from the Gene Expression Omnibus (GEO)  and ArrayExpress  public repositories. A systematic search of all published studies in public repositories (2002–2020) was conducted between June and September of 2020, following the preferred reporting items for systematic reviews and meta-analyses (PRISMA) guidelines . Keywords employed in the search were "Parkinson”, "Parkinson's Disease”, and "PD”. We applied the following inclusion criteria:
Transcriptomics studies on Homo sapiens.
Control and PD-affected subjects included.
Sex, disease status, and brain region variables registered.
RNA extracted directly from post-mortem brain tissues (no cell lines or cultures).
Brain tissues from either FC, ST, or SN.
Sample size > 3 for case and control groups in both sexes.
Finally, normalized gene expression data of 8 array PD datasets (E-MTAB-1194, E-MEXP-1416, GSE28894, GSE8397, GSE20295, GSE20159, GSE7621, and GSE20146) were retrieved using R packages GEOquery  and ArrayExpress .
Individual transcriptomics analysis
Individual transcriptomics analysis was performed on every selected study, which comprised 2 steps: preprocessing and differential expression analysis.
Data preprocessing included the standardization of the nomenclature for the clinical variables in each study, the homogenization of gene annotation, and an exploratory data analysis. The normalization methods performed by the original authors were assessed for each dataset, and data matrices were log2 transformed when necessary. All probe sets were annotated to HUGO gene symbols  using the biomaRt R package . When dealing with duplicated probe-to-symbol mappings, the median of their expression values was calculated. The exploratory analysis included unsupervised clustering and PCA to detect expression patterns between samples and genes and the presence of batch effects in each study. At this point, the GSE20159 study was excluded for presenting a strong batch effect with a critical impact on differential expression analysis.
Differential expression analyses were performed using the limma R package  to detect sex-based differentially expressed genes. To achieve this goal, the following comparison was applied:
This comparison allows the detection of genes with sex-based differential behavior in the development of PD. Genes with a Log2 fold change (LFC) greater than 0 show either a higher increase or a lesser decrease in expression in females when comparing the effect of the disease between sexes. For simplicity, these genes are referred to as increased in females. On the contrary, genes with an LFC lower than 0 have a higher increase or a lesser decrease in expression in males when comparing the effect of the disease between sexes. For simplicity, these genes are referred to as increased in males.
This comparison was applied to each brain region separately and all 3 regions together. When necessary, the batch effect was included as a categorical variable on the limma linear model to reduce its impact on data. P-values were corrected using the Benjamini–Hochberg procedure  and considered significant below a threshold of 0.05.
Gene expression meta-analysis
Differential gene expression results were integrated into a single meta-analysis  for each brain region (FC, ST, and SN) and a fourth meta-analysis for all regions. Meta-analyses were implemented with the R package metafor , under the DerSimonian and Laird random-effects model , considering individual study heterogeneity. This model considers the variability of individual studies by increasing the weights of studies with less variability when computing meta-analysis results. Thus, the most robust functions between studies are highlighted.
P-values, corrected p-values, LFC, LFC standard error (SE), and LFC 95% confidence intervals (CI) were calculated for each evaluated gene. Functions and pathways with corrected p-values of < 0.05 were considered significant, and both funnel and forest plots were computed for each. These representations were evaluated to assess for possible biased results, where LFC represents the effect size of a gene, and the SE of the LFC serves as a study precision measure . Sensitivity analysis (leave-one-out cross-validation) was conducted for each significant gene to verify alterations in the results due to the inclusion of any study. The Open Targets platform  (release 21.06) was used to explore the associations of the significant genes to PD.
Sex-based functional signature in the SN
Gene meta-analysis of SN data revealed sets of differentially expressed genes between male and female PD patients. Several analyses were carried out to identify the functional implications of these differences.
Over-Representation Analysis (ORA)  through the R package clusterProfiler  was first used to determine the biological functions and pathways overrepresented in the following gene sets: (i) differentially expressed "LFC > 0" genes; (ii) differentially expressed "LFC < 0" genes, and (iii) all differentially expressed genes. p-values and corrected p-values were calculated for each GO term from the 3 GO ontologies  and each KEGG pathway . Every function and pathway with a corrected p-value < 0.05 was labeled as overrepresented in each gene set.
Protein–protein interaction (PPI) networks were then calculated using the STRING web tool for each subset of genes . The total number of edges was examined, and PPI enrichment was assessed using the default parameters for each network.
Finally, a VIPER  analysis with human regulons obtained from the DoRothEA R package  was performed to estimate transcription factor activity using the consensus LFC of each gene evaluated in meta-analysis as gene expression signature. Regulons with a confidence level of A, B, C, or D were selected, excluding those with less than 25 genes (n = 217). The i-values were corrected with the Benjamini and Hochberg method. Normalized enrichment scores (NES) were calculated by VIPER as a measure of relative transcription factor activity.
Metafun-PD Web Tool
All data and results generated in the different steps of the meta-analysis are available in the Metafun-PD web tool , which is freely accessible to any user and allows the confirmation of the results described in this manuscript and the exploration of other results of interest. The front-end was developed using the Angular Framework, the interactive graphics used in this web resource have been implemented with plotly , and the exploratory analysis cluster plot was generated with the ggplot2 R package .
This easy-to-use resource is divided into 7 sections: (1) summary of analysis results in each phase. Then, for each of the studies, the detailed results of the (2) exploratory analysis and (3) differential expression. (4) The gene meta-analysis results of the four different performed meta-analyses. The user can interact with the web tool through graphics and tables and search for specific information for a gene or function. Finally, Sects. 5–7 provide the detailed tables and figures corresponding to the results of the 3 functional profiling methods (ORA, PPI, and VIPER).
To evaluate sex-based differences in PD, we performed a systematic review and 4 meta-analyses of transcriptomic studies that included information on patient sex from the GEO  and ArrayExpress databases . The meta-analyses comprised 1 for each of the primary brain regions affected by PD pathogenesis—FC (3 studies), ST (3 studies), and SN (f studies)—and a fourth global meta-analysis for the 3 brain regions combined (7 studies). Finally, we explored the biological implications of the SN meta-analysis results by applying 3 different functional profiling methods: over-representation analysis (ORA), the generation of protein–protein interactions (PPI) networks, and the transcription factors activity analysis (Fig. 1).
Study search and selection
The systematic review identified 83 non-duplicated studies, of which 39 (47%) included both male and female PD patients. We selected 8 comparable studies after applying inclusion and exclusion criteria (Methods, Fig. 2); however, we discarded 1 study after the exploratory analysis. Thus, we analyzed 7 studies that included 267 samples (132 controls and 135 PD cases) from the FC, ST, and SN brain regions (Table 1). Figure 3 describes sex distribution by study and region (overall, 59% males and 41% females); the median age was 75. Table 1 and Figure 2 contain further information regarding the selected studies and their sample clinicopathological characteristics.
We carried out exploratory and processing steps on the datasets to ensure their comparability in subsequent analyses. We applied log2 transformation to studies GSE28894, GSE20295, GSE7621, and GSE20146 to homogenize magnitude order and filtered out samples from regions different to FC, SN, and ST from all studies. Exploratory analysis revealed an anomalous behavior (outlier) by the A04022 sample from the E-MTAB-1194 study, which we excluded from further analysis.
The differential expression results for each study provided only a small number of significantly affected genes in the SN tissues from the GSE8397 study (Table 2).
Gene expression meta-analysis
We performed a gene expression meta-analysis in the FC, ST, and SN analyzing 3 studies for FC, 3 studies for ST, and 4 for SN. Finally, we performed a global meta-analysis that integrated all 7 studies from these brain regions. In the meta-analysis by regions, we found no significant genes in the FC, 1 significant gene (MED31—mediator of RNA polymerase II transcription subunit 31) in the ST, and 237 significant genes in the SN, reflecting the well-known impact of PD in this region. Of those, 75 genes exhibited increased expression in males (LFC < 0, see Methods) and 162 in females (LFC > 0) (Table 3). 16 and 44 of these genes in male and female PD patients, respectively, had a documented association with PD in the OpenTargets  database.
The global meta-analysis integrating the 7 studies and the 3 primary brain regions affected by PD revealed 15 differentially expressed genes by sex, 4 significant genes increased in females and 11 in males (Table 4). 1 gene of each group had already known associations with PD in the OpenTargets database.
Regarding those genes increased in males, we highlight AC073283.7 (a long non-coding RNA), ATP5J2-PTCD1 (a locus that represents naturally occurring read-through transcription between the genes ATP5J2 [ATP synthase, H + transporting, mitochondrial Fo complex, subunit F2] and PTCD1 [pentatricopeptide repeat domain 1]), GPR4 (G Protein-Coupled Receptor 4; may mediate central respiratory sensitivity to CO2 in the brain), PNMT (pentalenolactone synthase; cerebral disease and intellectual disability, hypotonia and mitochondrial disease, already linked to PD), PRDM12 (PR domain zinc finger protein 12; inflammation or degeneration of the sensory nerves), ITGA5 (Integrin Subunit Alpha 5; paraneoplastic neurologic syndrome), MEGF6 (Multiple EGF Like Domains 6 intracranial hemorrhage), LAMA5 (Laminin Subunit Alpha 5; cerebral diseases of vascular origin with epilepsy), MAP3K6 (Mitogen-Activated Protein Kinase 6, musculoskeletal or connective tissue disease genetic), IL15RA (Interleukin 15 Receptor Subunit Alpha; inflammatory brain disease), and RASIP1 (Ras Interacting Protein 1; juvenile primary lateral sclerosis [JPLS], a very rare motor neuron disease characterized by progressive upper motor neuron dysfunction leading to loss of the ability to walk; mitochondrial oxidative phosphorylation disorder with no known mechanism).
The significant genes with increased expression in females included OPTN (Optineurin; amyotrophic lateral sclerosis [ALS], a degenerative disorder affecting upper motor neurons in the brain, already linked to PD), UBA6 (Ubiquitin Like Modifier Activating Enzyme 6; developmental and epileptic encephalopathy), BZW1 (Basic Leucine Zipper And W2 Domains 1; Rolandic epilepsy—speech dyspraxia; intellectual disability), and TMX2 (Thioredoxin Related Transmembrane Protein 2; regulates mitochondrial activity; abnormality in the process of thought including the ability to process information; cognitive abilities or memory anomaly, intellectual disability in general).
Additional file 1: Table S1 details the results for all significant genes from all meta-analyses, including the adjusted p-value, the LFC and its 95% confidence interval (CI), and the standard error (SE) of the LFC.
Over-represented functions in the SN
ORA by R package clusterProfiler  on the SN significant genes revealed that enriched biological functions in the genes increased in males related to oxidative stress, including detoxification of inorganic compounds and stress response to metal ions (Fig. 4). In the genes increased in females, ORA overrepresented functions relate to antigen processing and presentation, pH regulation, proton transmembrane transport, mitochondrial functions, autophagy, cytoskeleton, and microtubule transport. These functions may create a neurodegenerative microenvironment that results in neuronal death. In addition, purine metabolism and ribose metabolism also appeared as enriched in this gene set.
Protein–protein interaction networks for SN genes
Using the STRING web tool, we created PPI networks for significant SN genes with positive and negative LFC . We found high connectivity in the PPI network generated for the 75 genes with negative LFC (p-value 3.87e−05) (Fig. 5a). Within the PPI analysis, we encountered the following protein clusters:
Angiogenesis cluster (in red): HSPG2, CSPG4, ITGA5, LAMA5, ENG, PDGFRG, APOC2, APOL1, and IGFBP7. A central cluster of proteins participates in angiogenic processes.
Metallothioneins cluster (in pearl): MT2A, MT1E, and MT1M are linked to the interferon-induced transmembrane proteins IFITM2 and IFITM3 through HLA-E.
STOM, TMBIM1, CYBA cluster (in green). These proteins have roles in redox homeostasis, response to oxidative stress, inflammation, and blood vessel remodeling.
We also encountered 3 independent pairs of interacting proteins: POLR1E-GTF2H4, MICALL2-RAB13, and ARPC1B-PDGFRB.
We also encountered a highly connected PPI network generated from the 162 genes with positive LFC (p-value < 1.0e−16) (Fig. 5b). Within this PPI network, we encountered the following protein clusters:
Microtubule and transport cluster (in green). This cluster includes proteins related to microtubule reorganization and transport (tubulins alpha: TUBA1B, TUBA3C, and beta TUBB4B, Myosin, and Actin).
V-ATPase and Secretogranin clusters (in blue). Formed by V-ATPases subunits involved in neuronal diseases and members of the chromogranin/secretogranin family of neuroendocrine secretory proteins.
Protein folding, recycling, and degradation (in red). Composed by proteasome subunits PSMB5, PSMB7; Ubiquitin UBB, Ubiquitin-protein hydrolase protein UCHL1, Chaperone HSPA8, and clathrin and adaptins (E-L system).
Mitochondria function and metabolism cluster (in yellow). Proteins such as GOT1 and GOT2 (glutamic metabolism) or VDAC2 and VDAC3, which are voltage-dependent anion-selective channel proteins involved in mitochondrial outer membrane permeabilization, negative regulation of intrinsic apoptotic signaling pathway, and negative regulation of protein polymerization, among other functions.
Transcription factor activity in the SN
The TF activity analysis of the significant SN genes performed with Dorothea  reported 40 differentially activated TFs (adjusted p-value < 0.05), with 29 more activated in male PD patients (NES < 0) and 11 more activated in female PD patients (NES > 0) (Fig. 6 and Additional file 3: Table S3). Among the 29 TFs activated in male PD patients, ADNP, JUN, MBD3, PRDM14, and ESR1 have known associations with neurodegenerative disorders such as Alzheimer's disease, mental retardation, schizophrenia, and autism. Additionally, a subset of specific TFs comprising BATF, CEBPA, ETS1, KLF6, LYL1, NFKB1, RELA, SOX13, SP3, SPI1, STAT1, STAT2, and STAT3 has documented involvement in immune and neuroinflammatory functions (Additional file 3: Table S3).
Specific TFs activated with an NES > 0 (female PD patients), such as ATF3, BCL6, or PCGF2, have been associated with neurodegenerative diseases or cognitive disabilities. TFs such as ZBTB7A, ZC3H8, or NME2 may be involved in processes related to neuroinflammation. These TFs also play additional roles in other processes, including immune response, hemopoiesis, and regulating the differentiation and activation of T and B lymphocytes (Additional file 3: Table S3). To the best of our knowledge, the remaining significant TFs remain poorly characterized or have not been studied in the context of PD.
Metafun-PD web tool
The Metafun-PD web tool  contains information about the 7 studies and 267 samples involved in this study. This resource includes statistical indicators of each performed analysis for each study, which users can explore to identify profiles of interest.
We carried out a total of 4 meta-analyses. For each of the significantly altered genes, Metafun-PD depicts the global activation level for all studies and each study's specific LFC, confidence interval of LFC, and p-value, and shows graphical representations by gene as forest and funnel plots. This open resource hopes to contribute to data sharing between researchers, the elaboration of innovative studies, and the discovery of new findings.
Meta-analysis of transcriptomic profiles
We carried out 3 gene expression meta-analyses in the 3 brain regions most affected by PD (FC, ST, and SN) and a fourth that integrated studies from all regions to better understand sex-based differences in PD. We identified the MED31 gene as significantly increased in females in the ST, 237 differentially expressed genes in the SN (75 and 162 significantly increased in males and females, respectively), a number that probably reflects the known impact of PD on the SN , but failed to find differences in the FC.
Functions associated with MED31 (mediator of RNA polymerase II transcription subunit 31) include the regulation of lipid metabolism by peroxisome proliferator-activated receptor alpha (PPARα) and gene expression [55, 56]. Dysregulated lipid metabolism supports the accumulation of α-synuclein and the formation of Lewy bodies in the brainstem, limbic system, and cortical areas . Long-chain fatty acids, ceramides, and lipopolysaccharides can also induce cellular stress and inflammatory responses in the brain . MED31 has been associated with neural diseases such as ALS types 3 and 5 , which affect motor neurons; however, this study represents the first association of MED31 with PD.
The genes increased in the SN of male PD patients can be clustered into several families:
Metallothioneins, including genes such as MT2A, MT1E, and MT1M. Metallothioneins are small cysteine-rich proteins that play essential roles in metal homeostasis and toxicity, DNA damage, and oxidative stress. Oxidative stress and inflammation (hallmarks of neurodegenerative diseases) influence the regulation of metallothionein expression. Metallothioneins have been described as potential markers of neurologic disease processes and treatment response in lysosomal storage disorders . Moreover, studies have encountered increased metallothionein expression in the SN and FC of PD patients .
Apolipoproteins—APOL1 and APOC2. These proteins control the binding of lipids to form lipoproteins, whose primary function involves lipid transport. Other lipoproteins, such as ApoE, have been related to aging and neurodegenerative diseases [62, 63], while lipid metabolism has strong links to neurotrophic disorders and correlates with symptoms in PD patients .
Interferon-induced transmembrane proteins—IFITM2 and IFITM3. These proteins play roles in the modulation of innate immunity and have been previously related to neurodegenerative diseases. In Alzheimer's disease, inflammatory cytokines induce IFITM3 expression in neurons and astrocytes; IFTM3 then binds to γ-secretase to upregulate activity, increasing the production of amyloid-β (the main component of the amyloid plaques found in Alzheimer's disease patients) [65, 66].
In summary, male PD patients present alterations in pathways related to oxidative stress, inflammation, and innate immune response, which represent hallmarks of neurodegeneration.
Genes with a higher expression in the SN of female PD patients can also be clustered into several families:
Secretogranins—SCG2, SCG3, SCG5, and CHGB. Members of this family control the delivery of peptides and neurotransmitters. Alterations in the granin family have been associated with PD .
V-ATPase subunits—ATP5B, ATP6AP1, ATP6V0D, ATP6V1B2, ATP6V1E1, and ATP6V1G2. Dysfunction of V-ATPase affects lysosomal acidification, which disrupts substrate clearance and leads to many disorders, including neurodegenerative diseases [68, 69].
TIM22 complex subunits—TIMM8B and TIMM10B form a complex involved in mitochondrial protein import. Alterations in TIM complexes such as TIM23 have been suggested as relevant mechanisms in neurodegenerative diseases .
Axonal transport and cytoskeleton stability: Dyneins—DYNC1LI1, DYNC1I1, and DYNLL1, Dynactin subunits—DCTN3, DCTN6, and ACTR10; Actin-related proteins—ACTR3B and ACTR10, and Tubulin subunits—TUBA1B, TUBA3C, and TUBB4B. Alterations in these complexes disrupt axonal transport, prompt the accumulation of misfolded proteins and motor neuron diseases, and are described in disorders such as ALS and other neurodegenerative processes [71,72,73].
Mitochondrial porins—VDAC2 and VDAC3 are responsible for voltage-dependent anion channels and mitochondrial dysfunction and contribute to neurodegenerative diseases .
Glutamic-oxaloacetic transaminase—cytoplasmic and mitochondrial forms GOT1 and GOT2. Metabolic processes also contribute to PD progression  and glutamate metabolism associated with excitotoxicity and neuron death .
Serine/threonine phosphatase subunits—PPP3CB and PPP2CA. PP3CB is a calcium-dependent calmodulin-stimulated protein phosphatase that plays an essential role in the transduction of intracellular Ca2+-mediated signals . PP3CB is related to biological processes such as axon extension, learning, locomotion, lymphangiogenesis, memory, regulation of synaptic plasticity and synaptic vesicle endocytosis, and cytokine and T cell responses. PPP2CA is the major phosphatase for microtubule-associated proteins (MAPs), which can affect GABA receptor binding, tau protein binding, regulation of apoptotic processes, and microtubule-binding as indicated in UniProtKB .
The features altered in the SN of female PD patients highlight the importance of acidification, microtubule stability, mitochondrial and lysosomal dysfunction, glutamic metabolism, and neurotoxicity to neurodegeneration and neuronal death in PD. The potential dysfunction of mitochondrial and lysosomal activity in PD remains of particular interest, with both systems playing a vital role in cellular redox homeostasis. Mitochondrial dysfunction can promote a decline in energy production, the increased generation of reactive oxygen species, and the induction of stress-induced apoptosis. Meanwhile, lysosomes participate in the turnover and degradation of organelles and proteins; targets such as the mitochondria and alpha-synuclein aggregates, respectively, may have relevance to PD [78,79,80].
Notably, the functional profiles inferred in male and female PD patients appear similar (oxidative stress, inflammation, and neurodegeneration); however, we found differences in the protein clusters defined by differential gene expression and the underlying mechanisms. While male PD patients present an environment related to metal homeostasis, lipid metabolism, and immunity, the female PD patients' environment exhibits mitochondrial and lysosomal dysfunction and alterations to cytoskeletal proteins and glutamic metabolism. Importantly, these data may help further understand PD development and guide personalized, sex-specific therapeutic interventions for PD patients.
Finally, we identified 15 differentially expressed genes in the global meta-analysis. Eleven of these genes appeared significantly increased in male PD patients (AC073283.7, ATP5J2-PTCD1, GPR4, IL15RA, ITGA5, LAMA5, MAP3K6, MEGF6, PNTM, PRDM12, and RASIP1), with the majority associated with oxidative stress, inflammation, and cerebral disorders . The remaining 4 genes displaying a significant increase in female PD patients (OPTN, UBA6, BZW1, and TMX2) have been linked to apoptosis, ubiquitination, and mitochondrial activity. OPTN interacts with adenovirus E3-14.7K protein and mediates apoptosis, inflammation, and vasoconstriction through tumor necrosis factor-alpha or Fas-ligand pathways . Several viruses have been related to the etiology of neural diseases, although the underlying mechanisms remain incompletely understood [83, 84]. Ubiquitin dysregulation can affect normal PINK1 and Parkin function, which govern mitochondrial quality control and mitophagy . In female subjects, the alteration of these genes results in mitochondrial dysfunction, stress condition, energy depletion, and necrotic cell death [86, 87]. Even though the genes involved differ, we note similar overall results in male and female PD patients: the overrepresentation of mitochondrial dysfunction, oxidative stress, and inflammation, which translates into neuronal death and cognitive/intellectual disorders. This study represents the first description of a sex-based association of these genes in PD, highlighting them as candidates for future studies.
We also note that many genes identified in the SN and global/three-region meta-analysis have been linked to other neural, intellectual, and cognitive disorders in the Open Targets database. In particular, specific genes have been linked to PD, which provides confidence to the meta-analysis results. Further exploration of these genes may open new perspectives for biomarker identification, early diagnostics, and therapeutic approaches in PD and related disorders.
SN sex-based functional profiling
Based on the differential gene expression analysis in SN, we applied different approaches to analyze the functional scenario. ORA demonstrated significantly more increased genes in male PD patients associated with metal ion detoxification and stress responses to several ions. Detailed studies have encountered metal ions in protein aggregates in PD brains [88, 89] and have demonstrated the contribution of metal ions to oxidative stress, toxicity, and degeneration of dopaminergic neurons in PD . In correlation with ORA results, PPI analyses identified metallothioneins, redox homeostasis, oxidative stress, and inflammation as the most prominent protein interaction clusters in male PD patients and angiogenesis. As interferon-stimulated genes [91, 92], metallothioneins may combine with Interferon-gamma to upregulate microglia genes, a mechanism observed in multiple system atrophy [93, 94]. The MICAL-like protein 2 Ras-related protein Rab is an effector of small Rab GTPs such as Rab13, which modulate alpha-synuclein levels, aggregation, and toxicity . The ARPC1B cluster may be involved in the cytolytic activity of CD8 cytotoxic T lymphocytes , which is related to brain inflammation. Several studies support the role of STOM/TMBIM1/CYBA proteins in redox homeostasis and oxidative stress, inflammation, and blood vessel remodeling [77, 97]. Finally, increased angiogenesis has been described in SN post-mortem tissues from PD patients .
In female PD patients, overrepresented functions correspond to microtubule and cytoskeletal transport, mitochondrial transport along the microtubule, macroautophagy, and several functions associated with antigen processing and presentation. α-Synuclein is a primary trigger of the immune response in PD  and may activate both the innate and adaptive immune system since the degree of microglial activation directly correlates with α-synuclein load in post-mortem brains . Purine metabolism is also overrepresented in the "LFC > 0" gene set, with levels of urate and other purines known to correlate with PD severity . Regarding purine-ribose dysregulation observed in females, a reported link to microglial activation (neuroinflammation) and neurodegeneration  could inspire the design of specific therapeutic strategies in male and female PD patients. As in males, PPI analyses in female PD patients identified clusters closely related to the gene families and ORA results. These clusters provided evidence for the critical role of mitochondria and metabolism in PD, as reported for other neurodegenerative diseases . Glutamate metabolism is vital for neuronal excitability, playing a critical role in memory, synaptic plasticity, and neuronal development ; however, glutamate overstimulation is also implicated in toxicity and neurodegeneration . V-ATPase subunits (ATP-dependent proton pumps present in both intracellular compartments and the plasma membrane) have significant involvement in neurodegenerative diseases. Proton pumps contribute to defective lysosomal acidification in lysosomal storage disorders and common neurodegenerative conditions such as Alzheimer's disease and PD . Supporting our results, the dysregulation of secretogranins in PD has also been reported . Finally, protein folding, recycling, and degradation could be related to the misfolding and aggregation of toxic α-synuclein [103, 104].
Overall, the PPI and ORA results correlated well; furthermore, sex-related differences in inflammation, mitochondrial dysfunction, and oxidative stress mechanisms have been documented in several other studies [14, 18,19,20], which provides additional support to our results. Importantly, we identify differences in the protein networks involved in these mechanisms between male and female PD patients, which could be crucial to designing personalized therapeutic strategies.
Transcription factors represent critical modulators of gene expression and pathways controlling responses to stimuli such as oxidative stress, microglial activation, chronic inflammation, neurotoxins, and DNA damage . A wide range of transcription factors have been linked to PD pathogenesis, and as with other cells, the development, maintenance, and survival of neurons depend on the precise control of gene expression . Activated transcription factors in male PD patients included ADNP, JUN, MBD3, PRDM14, and ESR1. In the OpenTargets database , alterations in ADNP have been associated with autism spectrum disorders, intellectual disability, dysmorphic features, and hypotonia , JUN with Alzheimer's disease, and PRDM14 with schizophrenia. Meanwhile, MBD3 is a critical neurodevelopmental transcription factor [108, 109] associated with neurogenesis and connectivity [108, 110]. Alterations in MBD3 have been related to PD  and to Neuropathy, Hereditary Sensory, Type Ie in the GeneCards database . Polymorphisms in ESR1 may contribute to increased PD susceptibility . Male PD patients also displayed the altered activity of SP1, which regulates the expression of LRRK2 , a contributing factor to PD pathogenesis; furthermore, a study has suggested that SP1 inhibition may provide beneficial effects in PD models . SMAD3 also presents higher activity in male PD patients and plays a vital role in PD, with protein deficiency known to reduce neurogenesis significantly. SMAD3 dysfunction leads to the formation of α-synuclein aggregates and a reduction in the number of dopaminergic axons and dendrites . Other transcription factors altered in male PD patients (i.e., BATF, CEBPA, ETS1, KLF6, LYL1, NFKB1, RELA, SOX13, SP3, SPI1, STAT1, STAT2, and STAT3) may participate in neuroinflammatory processes and regulate various functions of the immune response.
Of the transcription factors altered in female PD patients, ATF3, BCL6, and PCGF2 have been previously associated with neurodegenerative diseases or cognitive disabilities. Alterations in ATF3 in response to neurological damage and reactive oxygen species production have been demonstrated in a PD model [116, 117]. The BCL6 transcriptional repressor targets the ITM2B gene , which has important links to neurodegenerative diseases such as AD. The ITM2B protein binds amyloid precursor protein and inhibits processing, thereby reducing the secretion and accumulation of beta-amyloid peptides. PCGF2 is associated with phenotypes such as intellectual disability, global developmental delay, and mental retardation in the OpenTargets database . Additionally, the ZBTB7A, ZC3H8, or NME2 transcription factors have been linked to neuroinflammation processes. These transcription factors play roles in multiple processes, including the immune response, modulating hemopoiesis, and regulating the differentiation and activation of T and B lymphocytes (Additional file 3: Table S3). The transcription factors identified in this study may represent biomarkers for the detection of neuron degeneration, although this will require more in-depth validation studies.
Strengths and limitations
We performed an in silico approach to study sex-based differences in PD. In silico strategies using computational models represent a powerful tool to evaluate and integrate data, and the results obtained possess greater consensus and statistical power. In this approach, the sample size increases with an augmentation of the number of studies integrated into the meta-analysis; therefore, more subtle effects can be detected. Additionally, we based our in silico analysis of PD with a sex perspective on FAIR data (Findable, Accessible, Interoperable, Reusable) , which we believe to be critically important. Indeed, we strongly believe that research data should be legally sharable and reusable in new research. In silico integrative approaches to analyzing sex-based differences in PD have been carried out before; for example, Mariani et al.  systematically meta-analyzed SN microarray data using the Transcriptome Mapper software (TRAM version 1.2). Our study also includes 2 other primarily affected brain regions (FC and ST) and has been carried out using the limma R package for differential expression analyses  and the metafor R package for gene expression meta-analysis . This allowed us to analyze data following our sex differences comparison, which includes 4 experimental groups, and identifies sex differences in PD, considering the inherent variability among males and females in healthy conditions.
Sex-based differences influence progression, symptoms, treatment response, and mortality in many diseases, including PD; however, the segregation of data by sex does not always occur in research studies, not even in those carried out to explore diagnostic or prognostic factors. For example, we excluded 44 (53%) studies from our systematic review due to the absence of sex information. We highlight the need to include sex information in research studies and databases, as they are relevant to health.
Perspectives and significance
The novel strategy presented in this work and the results obtained, help to clarify the molecular mechanisms associated with sex differences in PD. The identification of consensus expression profiles specific to men and women is relevant because it generates useful knowledge for diagnosis, prognosis and drug repositioning that take sex information into account. Likewise, the identification of differential biological processes between men and women with PD allows other researchers to compare and evaluate their results in order to understand the functioning of PD in their own work. This work opens up numerous research perspectives that will have an impact on improving the treatment of patients with PD.
In conclusion, our in silico approach has highlighted sex-based differential mechanisms in typical PD hallmarks (inflammation, mitochondrial dysfunction, and oxidative stress). Additionally, we have identified specific genes and transcription factors for male and female PD patients that represent potential candidates as biomarkers to diagnosis. Further studies that consider the sex perspective are urgently required to better understand PD and develop tailored interventions that consider the distinct requirements of male and female PD patients. Finally, we underscore the essential nature of sharing data and using open platforms for scientific advancement.
The data used for the analyses described in this work are publicly available at GEO  and ArrayExpress . The accession numbers of the GEO datasets downloaded are GSE7621, GSE8397, GSE20146, GSE20295, and GSE28894. The accession numbers of the ArrayExpress datasets downloaded are E-MTAB-1194 and E-MEXP-1416.
The code developed for the analyses described in this work has been made publicly available at GitHub (https://github.com/adlpecer/PD_GeneMetaanalysis). All software and versions used are detailed in Additional file 2: Table S2.
Gene Expression Omnibus
Log2 fold change
Normalized enrichment score
Preferred Reporting Items for Systematic Reviews and Meta-Analyses
Benito-Leon J. Epidemiology of Parkinson’s disease in Spain and its contextualisation in the world. Rev Neurol. 2018;66(4):125–34.
Poewe W, Seppi K, Tanner CM, Halliday GM, Brundin P, Volkmann J, et al. Parkinson disease. Nat Rev Dis Primer. 2017;3(1):17013.
Heller J, Dogan I, Schulz JB, Reetz K. Evidence for gender differences in cognition, emotion and quality of life in Parkinson’s disease? Aging Dis. 2014;5(1):63–75.
Hanagasi HA, Tufekcioglu Z, Emre M. Dementia in Parkinson’s disease. J Neurol Sci. 2017;15(374):26–31.
Leszek J, Barreto GE, Gąsiorowski K, Koutsouraki E, Ávila-Rodrigues M, Aliev G. Inflammatory Mechanisms and Oxidative Stress as Key Factors Responsible for Progression of Neurodegeneration: Role of Brain Innate Immune System. CNS Neurol Disord Drug Targets. 2016;15(3):329–36.
Cabezas R, Ãvila M, Gonzalez J, El-Bach RS, Bez E, GarcaSegura LM, et al. Astrocytic modulation of blood brain barrier: perspectives on Parkinson disease. Front Cell Neurosci. 2014. https://doi.org/10.3389/fncel.2014.00211/abstract.
Fukae J, Mizuno Y, Hattori N. Mitochondrial dysfunction in Parkinson’s disease. Mitochondrion. 2007;7(1–2):58–62.
Lee Mosley R, Benner EJ, Kadiu I, Thomas M, Boska MD, Hasan K, et al. Neuroinflammation, oxidative stress, and the pathogenesis of Parkinson’s disease. Clin Neurosci Res. 2006;6(5):261–81.
Misiak M, Beyer C, Arnold S. Gender-specific role of mitochondria in the vulnerability of 6-hydroxydopamine-treated mesencephalic neurons. Biochim Biophys Acta BBA - Bioenerg. 2010;1797(6–7):1178–88.
Tysnes OB, Storstein A. Epidemiology of Parkinson’s disease. J Neural Transm Vienna Austria. 2017;124(8):901–5.
Picillo M, Nicoletti A, Fetoni V, Garavaglia B, Barone P, Pellecchia MT. The relevance of gender in Parkinson’s disease: a review. J Neurol. 2017;264(8):1583–607.
Hirsch L, Jette N, Frolkis A, Steeves T, Pringsheim T. The incidence of Parkinson’s disease: a systematic review and meta-analysis. Neuroepidemiology. 2016;46(4):292–300.
Dahodwala N, Pei Q, Schmidt P. Sex differences in the clinical progression of Parkinson’s Disease. J Obstet Gynecol Neonatal Nurs. 2016;45(5):749–56.
Colombo D, Abbruzzese G, Antonini A, Barone P, Bellia G, Franconi F, et al. The “Gender Factor” in Wearing-Off among Patients with Parkinson’s Disease: A Post Hoc Analysis of DEEP Study. Sci World J. 2015;2015:1–10.
Bjørnarå KA, Dietrichs E, Toft M. REM sleep behavior disorder in Parkinson’s disease – Is there a gender difference? Parkinsonism Relat Disord. 2013;19(1):120–2.
Betancourt E, Wachtel J, Michaelos M, Haggerty M, Conforti J, Kritzer MF. The impact of biological sex and sex hormones on cognition in a rat model of early, pre-motor Parkinson’s disease. Neuroscience. 2017;345:297–314.
Martinez-Martin P, Rodriguez-Blazquez C, Kurtis MM, Chaudhuri KR. The impact of non-motor symptoms on health-related quality of life of patients with Parkinson’s disease: Nms and HRQ O L in Parkinson’s Disease. Mov Disord. 2011;26(3):399–406.
Haaxma CA, Bloem BR, Borm GF, Oyen WJG, Leenders KL, Eshuis S, et al. Gender differences in Parkinson’s disease. J Neurol Neurosurg Amp Psychiatry. 2007;78(8):819–24.
Loke H, Harley V, Lee J. Biological factors underlying sex differences in neurological disorders. Int J Biochem Cell Biol. 2015;65:139–50.
Mitra S, Chakrabarti N, Dutta SS, Ray S, Bhattacharya P, Sinha P, et al. Gender-specific brain regional variation of neurons, endogenous estrogen, neuroinflammation and glial cells during rotenone-induced mouse model of Parkinson’s disease. Neuroscience. 2015;292:46–70.
Gillies GE, Pienaar IS, Vohra S, Qamhawi Z. Sex differences in Parkinson’s disease. Front Neuroendocrinol. 2014;35(3):370–84.
R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing. https://www.R-project.org
Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. 2012;41(D1):D991–5.
Athar A, Füllgrabe A, George N, Iqbal H, Huerta L, Ali A, et al. ArrayExpress update - from bulk to single-cell expression data. Nucleic Acids Res. 2019;47(D1):D711–5.
Moher D, Liberati A, Tetzlaff J, Altman DG. Preferred Reporting Items for Systematic Reviews and Meta-Analyses: The PRISMA Statement. PLoS Med. 2009;6(7):e1000097.
Davis S, Meltzer PS. GEOquery: a bridge between the Gene Expression Omnibus (GEO) and BioConductor. Bioinforma Oxf Engl. 2007;23(14):1846–7.
Kauffmann A, Rayner TF, Parkinson H, Kapushesky M, Lukk M, Brazma A, et al. Importing ArrayExpress datasets into R/Bioconductor. Bioinforma Oxf Engl. 2009;25(16):2092–4.
Tweedie S, Braschi B, Gray K, Jones TEM, Seal RL, Yates B, et al. Genenames.org: the HGNC and VGNC resources in 2021. Nucleic Acids Res. 2021;49(1):D939–46.
Durinck S, Spellman PT, Birney E, Huber W. Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt. Nat Protoc. 2009;4(8):1184–91.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47–e47.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995;57(1):289–300.
García-García F. Methods of functional enrichment analysis in genomic studies [PhD Thesis]. 2016.
Viechtbauer W. Conducting Meta-Analyses in R with the metafor Package. J Stat Softw. 2010;36:3.
DerSimonian R, Laird N. Meta-analysis in clinical trials. Control Clin Trials. 1986;7(3):177–88.
Sterne JAC, Egger M. Funnel plots for detecting bias in meta-analysis. J Clin Epidemiol. 2001;54(10):1046–55.
Ochoa D, Hercules A, Carmona M, Suveges D, Gonzalez-Uriarte A, Malangone C, et al. Open Targets Platform: supporting systematic drug–target identification and prioritisation. Nucleic Acids Res. 2021;49(D1):D1302–10.
Boyle EI, Weng S, Gollub J, Jin H, Botstein D, Cherry JM, et al. GO::TermFinder–open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics. 2004;20(18):3710–5.
Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z, et al. clusterProfiler 4.0: A universal enrichment tool for interpreting omics data. Innovation. 2021;100141:78.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene Ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25–9.
Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 1999;27(1):29–34.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.
Alvarez MJ, Shen Y, Giorgi FM, Lachmann A, Ding BB, Ye BH, et al. Functional characterization of somatic mutations in cancer using network-based inference of protein activity. Nat Genet. 2016;48(8):838–47.
Garcia-Alonso L, Holland CH, Ibrahim MM, Turei D, Saez-Rodriguez J. Benchmark and integration of resources for the estimation of human transcription factor activities. Genome Res. 2019;29(8):1363–75.
López-Cerdán A, Andreu Z, Hidalgo MR, Grillo-Risco R, Català-Senent JF, Soler-Sáez I, et al. Parkinson disease. http://bioinfo.cipf.es/metafun-pd/
Sievert C. Interactive Web-Based Data Visualization with R, plotly, and shiny. Chapman and Hall/CRC; 2020. https://plotly-r.com
Wickham H. ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York; 2016. https://ggplot2.tidyverse.org
Garcia-Esparcia P, Schlüter A, Carmona M, Moreno J, Ansoleaga B, Torrejón-Escribano B, et al. Functional genomics reveals dysregulation of cortical olfactory receptors in Parkinson disease: novel putative chemoreceptors in the human brain. J Neuropathol Exp Neurol. 2013;72(6):524–39.
Cantuti-Castelvetri I, Keller-McGandy C, Bouzou B, Asteris G, Clark TW, Frosch MP, et al. Effects of gender on nigral gene expression and parkinson disease. Neurobiol Dis. 2007;26(3):606–14.
Moran LB, Duke DC, Deprez M, Dexter DT, Pearce RKB, Graeber MB. Whole genome expression profiling of the medial and lateral substantia nigra in Parkinson’s disease. Neurogenetics. 2006;7(1):1–11.
Duke DC, Moran LB, Pearce RKB, Graeber MB. The medial and lateral substantia nigra in Parkinson’s disease: mRNA profiles associated with higher brain tissue vulnerability. Neurogenetics. 2007;8(2):83–94.
Zhang Y, James M, Middleton FA, Davis RL. Transcriptional analysis of multiple brain regions in Parkinson’s disease supports the involvement of specific protein processing, energy metabolism, and signaling pathways, and suggests novel disease mechanisms. Am J Med Genet Part B Neuropsychiatr Genet Off Publ Int Soc Psychiatr Genet. 2005;137B(1):5–16.
Zheng B, Liao Z, Locascio JJ, Lesniak KA, Roderick SS, Watt ML, et al. PGC-1α, a potential therapeutic target for early intervention in Parkinson’s disease. Sci Transl Med. 2010;2(52):52ra73.
Lesnick TG, Papapetropoulos S, Mash DC, Ffrench-Mullen J, Shehadeh L, de Andrade M, et al. A genomic pathway approach to a complex disease: axon guidance and Parkinson disease. PLoS Genet. 2007;3(6): e98.
Viveros-Paredes JM, Gonzalez-Castañeda RE, Escalante-Castañeda A, Tejeda-Martínez AR, Castañeda-Achutiguí F, Flores-Soto ME. Effect of inhibition of fatty acid amide hydrolase on MPTP-induced dopaminergic neuronal damage. Neurol Barc Spain. 2019;34(3):143–52.
Beadle EP, Straub JA, Bunnell BA, Newman JJ. MED31 involved in regulating self-renewal and adipogenesis of human mesenchymal stem cells. Mol Biol Rep. 2018;45(5):1545–50.
Schneider M, Hellerschmied D, Schubert T, Amlacher S, Vinayachandran V, Reja R, et al. The Nuclear Pore-Associated TREX-2 complex employs mediator to regulate gene expression. Cell. 2015;162(5):1016–28.
Alecu I, Bennett SAL. Dysregulated Lipid Metabolism and Its Role in α-Synucleinopathy in Parkinson’s Disease. Front Neurosci. 2019;13:328.
Mullins CA, Gannaban RB, Khan MS, Shah H, Siddik MAB, Hegde VK, et al. Neural underpinnings of obesity: the role of oxidative stress and inflammation in the brain. Antioxid Basel Switz. 2020;9(10):E1018.
Stelzer G, Rosen N, Plaschkes I, Zimmerman S, Twik M, Fishilevich S, et al. The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses. Curr Protoc Bioinforma. 2016. https://doi.org/10.1002/cpbi.5.
Cesani M, Cavalca E, Macco R, Leoncini G, Terreni MR, Lorioli L, et al. Metallothioneins as dynamic markers for brain disease in lysosomal disorders. Ann Neurol. 2014;75(1):127–37.
Michael GJ, Esmailzadeh S, Moran LB, Christian L, Pearce RKB, Graeber MB. Up-regulation of metallothionein gene expression in Parkinsonian astrocytes. Neurogenetics. 2011;12(4):295–305.
Yin Y, Wang Z. ApoE and Neurodegenerative Diseases in Aging. In: Wang Z, editor. Aging and Aging-Related Diseases. Singapore: Springer; 2018. p. 77–92.
Krasemann S, Madore C, Cialic R, Baufeld C, Calcagno N, El Fatimy R, et al. The TREM2-APOE pathway drives the transcriptional phenotype of dysfunctional microglia in neurodegenerative diseases. Immunity. 2017;47(3):566-581.e9.
Dong MX, Wei YD, Hu L. The disturbance of lipid metabolism is correlated with neuropsychiatric symptoms in patients with Parkinson’s disease. Chem Phys Lipids. 2021;239: 105112.
Hur JY, Frost GR, Wu X, Crump C, Pan SJ, Wong E, et al. The innate immunity protein IFITM3 modulates γ-secretase in Alzheimer’s disease. Nature. 2020;586(7831):735–40.
Yao AY, Yan R. Activity of Alzheimer’s γ-secretase is linked to changes of interferon-induced transmembrane proteins (IFITM) in innate immunity. Mol Neurodegener. 2020;15(1):69.
Rotunno MS, Lane M, Zhang W, Wolf P, Oliva P, Viel C, et al. Cerebrospinal fluid proteomics implicates the granin family in Parkinson’s disease. Sci Rep. 2020;10(1):2479.
Song Q, Meng B, Xu H, Mao Z. The emerging roles of vacuolar-type ATPase-dependent Lysosomal acidification in neurodegenerative diseases. Transl Neurodegener. 2020;9(1):17.
Bagh MB, Peng S, Chandra G, Zhang Z, Singh SP, Pattabiraman N, et al. Misrouting of v-ATPase subunit V0a1 dysregulates lysosomal acidification in a neurodegenerative lysosomal storage disease model. Nat Commun. 2017;8(1):14612.
Yablonska S, Ganesan V, Ferrando LM, Kim J, Pyzel A, Baranova OV, et al. Mutant huntingtin disrupts mitochondrial proteostasis by interacting with TIM23. Proc Natl Acad Sci. 2019;116(33):16593–602.
Tanaka F, Waza M, Yamamoto M, Sobue G. Exploration of pathogenesis and therapy development for ALS employing sporadic disease model. Rinsho Shinkeigaku. 2009;49(11):811–3.
Tanaka F, Ikenaka K, Yamamoto M, Sobue G. Neuropathology and omics in motor neuron diseases: Neuropathology and Omics in MND. Neuropathology. 2012;32(4):458–62.
Eira J, Silva CS, Sousa MM, Liz MA. The cytoskeleton as a novel therapeutic target for old neurodegenerative disorders. Prog Neurobiol. 2016;141:61–82.
Caterino M, Ruoppolo M, Mandola A, Costanzo M, Orrù S, Imperlini E. Protein–protein interaction networks as a new perspective to evaluate distinct functional roles of voltage-dependent anion channel isoforms. Mol Biosyst. 2017;13(12):2466–76.
Venkatesan D, Iyer M, Narayanasamy A, Siva K, Vellingiri B. Kynurenine pathway in Parkinson’s disease—An update. Neurological Sci. 2020;21:100270.
Iovino L, Tremblay ME, Civiero L. Glutamate-induced excitotoxicity in Parkinson’s disease: The role of glial cells. J Pharmacol Sci. 2020;144(3):151–64.
Bateman A, Martin MJ, Orchard S, Magrane M, Agivetova R, et al. UniProt: the universal protein knowledgebase in 2021. Nucleic Acids Res. 2021;49(1):D480–9.
Kao AW, McKay A, Singh PP, Brunet A, Huang EJ. Progranulin, lysosomal regulation and neurodegenerative disease. Nat Rev Neurosci. 2017;18(6):325–33.
Mahajan UV, Varma VR, Griswold ME, Blackshear CT, An Y, Oommen AM, et al. Dysregulation of multiple metabolic networks related to brain transmethylation and polyamine pathways in Alzheimer disease: A targeted metabolomic and transcriptomic study. PLoS Med. 2020;17(1): e1003012.
Ryter SW, Bhatia D, Choi ME. Autophagy: a lysosome-dependent process with implications in cellular redox homeostasis and human disease. Antioxid Redox Signal. 2019;30(1):138–59.
Rocha EM, De Miranda B, Sanders LH. Alpha-synuclein: Pathology, mitochondrial dysfunction and neuroinflammation in Parkinson’s disease. Neurobiol Dis. 2018;109(Pt B):249–57.
Klingseisen L, Ehrenschwender M, Heigl U, Wajant H, Hehlgans T, Schütze S, et al. E3–14.7K is recruited to TNF-receptor 1 and blocks TNF cytolysis independent from interaction with optineurin. PLoS ONE. 2012;7(6):e38348.
You L, Chen J, Liu W, Xiang Q, Luo Z, Wang W, et al. Enterovirus 71 induces neural cell apoptosis and autophagy through promoting ACOX1 downregulation and ROS generation. Virulence. 2020;11(1):537–53.
He Z, An S, Chen J, Zhang S, Tan C, Yu J, et al. Neural progenitor cell pyroptosis contributes to Zika virus-induced brain atrophy and represents a therapeutic target. Proc Natl Acad Sci U S A. 2020;117(38):23869–78.
Pickrell AM, Youle RJ. The roles of PINK1, parkin, and mitochondrial fidelity in Parkinson’s disease. Neuron. 2015;85(2):257–73.
Kramer NJ, Haney MS, Morgens DW, Jovičić A, Couthouis J, Li A, et al. CRISPR–Cas9 screens in human cells and primary neurons identify modifiers of C9ORF72 dipeptide-repeat-protein toxicity. Nat Genet. 2018;50(4):603–12.
Fiesel FC, Moussaud-Lamodière EL, Ando M, Springer W. Select E2 enzymes differentially regulate parkin activation and mitophagy. J Cell Sci. 2014;45:147520.
Moons R, Konijnenberg A, Mensch C, Van Elzen R, Johannessen C, Maudsley S, et al. Metal ions shape α-synuclein. Sci Rep. 2020;10(1):16293.
Albertson BD, Frederick KL, Maronian NC, Feuillan P, Schorer S, Dunn JF, et al. The effect of ketoconazole on steroidogenesis: I. Leydig cell enzyme activity in vitro. Res Commun Chem Pathol Pharmacol. 1988;61(1):17–26.
Dias V, Junn E, Mouradian MM. The role of oxidative stress in Parkinson’s Disease. J Park Dis. 2013;3(4):461–91.
Schneider WM, Chevillotte MD, Rice CM. Interferon-stimulated genes: a complex web of host defenses. Annu Rev Immunol. 2014;32(1):513–45.
Boehm U, Klamp T, Groot M, Howard JC. Cellular responses to interferon-γ. Annu Rev Immunol. 1997;15(1):749–95.
Moran LB, Duke DC, Graeber MB. The microglial gene regulatory network activated by interferon-gamma. J Neuroimmunol. 2007;183(1–2):1–6.
Inoue Y, Ayaki T, Ishimoto T, Yamakado H, Maki T, Matsuzawa S, et al. The stimulator of interferon genes (STING) pathway is upregulated in striatal astrocytes of patients with multiple system atrophy. Neurosci Lett. 2021;13(757): 135972.
Gonçalves SA, Outeiro TF. Traffic jams and the complex role of α-Synuclein aggregation in Parkinson disease. Small GTPases. 2017;8(2):78–84.
Randzavola LO, Strege K, Juzans M, Asano Y, Stinchcombe JC, Gawden-Bone CM, et al. Loss of ARPC1B impairs cytotoxic T lymphocyte maintenance and cytolytic activity. J Clin Invest. 2019;129(12):5600–14.
Gong F, Gu J, Wang H. Up regulated Tmbim1 activation promotes high fat diet (HFD)-induced cardiomyopathy by enhancement of inflammation and oxidative stress. Biochem Biophys Res Commun. 2018;504(4):797–804.
Desai Bradaric B, Patel A, Schneider JA, Carvey PM, Hendey B. Evidence for angiogenesis in Parkinson’s disease, incidental Lewy body disease, and progressive supranuclear palsy. J Neural Transm. 2012;119(1):59–71.
Allen Reish HE, Standaert DG. Role of α-synuclein in inducing innate and adaptive immunity in Parkinson Disease. J Park Dis. 2015;5(1):1–19.
Von Seggern M, Szarowicz C, Swanson M, Cavotta S, Pike ST, Lamberts JT. Purine molecules in Parkinson’s disease: Analytical techniques and clinical implications. Neurochem Int. 2020;139: 104793.
Zhang Z, Zhang S, Fu P, Zhang Z, Lin K, Ko JKS, et al. Roles of Glutamate Receptors in Parkinson’s Disease. Int J Mol Sci. 2019;20(18):E4391.
Ferraguti F, Crepaldi L, Nicoletti F. Metabotropic glutamate 1 receptor: current concepts and perspectives. Pharmacol Rev. 2008;60(4):536–81.
Ebrahimi-Fakhari D, Saidi LJ, Wahlster L. Molecular chaperones and protein folding as therapeutic targets in Parkinson’s disease and other synucleinopathies. Acta Neuropathol Commun. 2013;5(1):79.
Gonçalves SA, Macedo D, Raquel H, Simões PD, Giorgini F, Ramalho JS, et al. shRNA-Based Screen Identifies Endocytic Recycling Pathway Components That Act as Genetic Modifiers of Alpha-Synuclein Aggregation, Secretion and Toxicity. PLoS Genet. 2016;12(4): e1005995.
Tao K, Wang B, Feng D, Zhang W, Lu F, Lai J, et al. Salidroside Protects Against 6-Hydroxydopamine-Induced Cytotoxicity by Attenuating ER Stress. Neurosci Bull. 2016;32(1):61–9.
Wang R, Yang S, Nie T, Zhu G, Feng D, Yang Q. Transcription factors: potential cell death markers in Parkinson’s Disease. Neurosci Bull. 2017;33(5):552–60.
Arnett AB, Rhoads CL, Hoekzema K, Turner TN, Gerdts J, Wallace AS, et al. The autism spectrum phenotype in ADNP syndrome. Autism Res Off J Int Soc Autism Res. 2018;11(9):1300–10.
Knock E, Pereira J, Lombard PD, Dimond A, Leaford D, Livesey FJ, et al. The methyl binding domain 3/nucleosome remodelling and deacetylase complex regulates neural cell fate determination and terminal differentiation in the cerebral cortex. Neural Develop. 2015;10(1):13.
Hoffmann A, Spengler D. Chromatin Remodeling Complex NuRD in neurodevelopment and neurodevelopmental disorders. Front Genet. 2019;24(10):682.
Moon BS, Yun HM, Chang WH, Steele BH, Cai M, Choi SH, et al. Smek promotes corticogenesis through regulating Mbd3’s stability and Mbd3/NuRD complex recruitment to genes associated with neurogenesis. PLoS Biol. 2017;15(5): e2001220.
Liu X, Cheng R, Ye X, Verbitsky M, Kisselev S, Mejia-Santana H, et al. Increased rate of sporadic and recurrent rare genic copy number variants in Parkinson’s disease among Ashkenazi Jews. Mol Genet Genomic Med. 2013;1(3):142–54.
Chung SJ, Armasu SM, Biernacka JM, Lesnick TG, Rider DN, Cunningham JM, et al. Variants in estrogen-related genes and risk of Parkinson’s disease. Mov Disord Off J Mov Disord Soc. 2011;26(7):1234–42.
Wang J, Song W. Regulation of LRRK2 promoter activity and gene expression by Sp1. Mol Brain. 2016;9(1):33.
Yao L, Dai X, Sun Y, Wang Y, Yang Q, Chen X, et al. Inhibition of transcription factor SP1 produces neuroprotective effects through decreasing MAO B activity in MPTP/MPP+ Parkinson’s disease models. J Neurosci Res. 2018;96(10):1663–76.
Muñoz M, de la Fuente N, Sánchez-Capelo A. TGF-β/Smad3 Signalling Modulates GABA neurotransmission: implications in Parkinson’s Disease. Int J Mol Sci. 2020;21(2):590.
Holtz WA, Turetzky JM, Jong YJI, O’Malley KL. Oxidative stress-triggered unfolded protein response is upstream of intrinsic cell death evoked by parkinsonian mimetics. J Neurochem. 2006;99(1):54–69.
Zhao Q, Yang X, Cai D, Ye L, Hou Y, Zhang L, et al. Echinacoside Protects Against MPP+-Induced Neuronal Apoptosis via ROS/ATF3/CHOP Pathway Regulation. Neurosci Bull. 2016;32(4):349–62.
Baron BW, Pytel P. Expression Pattern of the BCL6 and ITM2B proteins in normal human brains and in Alzheimer Disease. Appl Immunohistochem Mol Morphol. 2017;25(7):489–96.
Wilkinson MD, Dumontier M, Aalbersberg IJJ, Appleton G, Axton M, Baak A, et al. The FAIR Guiding Principles for scientific data management and stewardship. Sci Data. 2016;15(3): 160018.
Mariani E, Lombardini L, Facchin F, Pizzetti F, Frabetti F, Tarozzi A, et al. Sex-specific transcriptome differences in substantia nigra tissue: a meta-analysis of Parkinson’s Disease data. Genes. 2018;9(6):275.
The authors thank the Principe Felipe Research Center (CIPF) for providing access to the cluster, which is co-funded by European Regional Development Funds (ERDF) in Valencian Community 2014-2020. The authors also thank Stuart P. Atkinson for reviewing the manuscript.
This research was supported by and partially funded by the Institute of Health Carlos III (project IMPaCT-Data, exp. IMP/00019), co-funded by the European Union, European Regional Development Fund (ERDF, “A way to make Europe”), and PID2021-124430OA-I00 funded by MCIN/AEI/10.13039/501100011033/FEDER,UE (“A way to make Europe”).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
Results for all significant genes from all meta-analyses. It includes the adjusted p-value, the LFC and its 95% confidence interval (CI), and the standard error (SE) of the LFC.
Additional file 2: Table S2.
List of R packages and versions used in this work.
Additional file 3: Table S3.
VIPER results for each TF evaluated. It includes the adjusted p-value, the NES, the size of the regulons and their functional annotation (From Gene Ontology Biological Processes).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
López-Cerdán, A., Andreu, Z., Hidalgo, M.R. et al. Unveiling sex-based differences in Parkinson's disease: a comprehensive meta-analysis of transcriptomic studies. Biol Sex Differ 13, 68 (2022). https://doi.org/10.1186/s13293-022-00477-5
- Parkinson disease
- Sex differences