Next Article in Journal
Cancer-Associated Fibroblast Heterogeneity in Malignancy with Focus on Oral Squamous Cell Carcinoma
Next Article in Special Issue
Phosphoproteome Microarray Analysis of Extracellular Particles as a Tool to Explore Novel Biomarker Candidates for Alzheimer’s Disease
Previous Article in Journal
Density Functional Theory Calculations: A Useful Tool to Investigate Mechanisms of 1,3-Dipolar Cycloaddition Reactions
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 25
Issue 2
10.3390/ijms25021299
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Network Analysis Performed on Transcriptomes of Parkinson’s Disease Patients Reveals Dysfunction in Protein Translation

by
Simone D’Angiolini
,
Maria Lui
,
Emanuela Mazzon
* and
Marco Calabrò
IRCCS Centro Neurolesi “Bonino-Pulejo”, Via Provinciale Palermo, Contrada Casazza, 98124 Messina, Italy
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(2), 1299; https://doi.org/10.3390/ijms25021299
Submission received: 24 December 2023 / Revised: 17 January 2024 / Accepted: 19 January 2024 / Published: 21 January 2024
(This article belongs to the Special Issue Peripheral Biomarkers in Neurodegenerative Diseases—4th Edition)
Browse Figures
Versions Notes

Abstract

:
Parkinson’s disease (PD) is a prevalent neurodegenerative disorder characterized by the progressive degeneration of dopaminergic neurons in the substantia nigra region of the brain. The hallmark pathological feature of PD is the accumulation of misfolded proteins, leading to the formation of intracellular aggregates known as Lewy bodies. Recent data evidenced how disruptions in protein synthesis, folding, and degradation are events commonly observed in PD and may provide information on the molecular background behind its etiopathogenesis. In the present study, we used a publicly available transcriptomic microarray dataset of peripheral blood of PD patients and healthy controls (GSE6613) to investigate the potential dysregulation of elements involved in proteostasis-related processes at the transcriptomic level. Our bioinformatics analysis revealed 375 differentially expressed genes (DEGs), of which 281 were down-regulated and 94 were up-regulated. Network analysis performed on the observed DEGs highlighted a cluster of 36 elements mainly involved in the protein synthesis processes. Different enriched ontologies were related to translation initiation and regulation, ribosome structure, and ribosome components nuclear export. Overall, this data consistently points to a generalized impairment of the translational machinery and proteostasis. Dysregulation of these mechanics has been associated with PD pathogenesis. Understanding the precise regulation of such processes may shed light on the molecular mechanisms of PD and provide potential data for early diagnosis.

1. Introduction

Parkinson’s disease (PD) ranks as the second most prevalent neurodegenerative disorder after Alzheimer’s disease, reaching ~8.5 × 103 cases worldwide [1]. The clinical phenotype of this pathology is characterized by a distinctive pattern of impairment involving primary motor symptoms, secondary motor symptoms, and non-motor symptoms [2]. The motor syndrome is characterized by movement dysfunctions like slow movement (bradykinesia), reduced movement (hypokinesia), or even an inability to move (akinesia). Another distinctive motor feature of PD is the resting tremor, like trembling or shaking of the hands, but also affecting limbs, legs, feet, the face, and arms [3]. Patients suffering from this multifactorial disease also experience secondary motor symptoms such as gait rigidity, postural instability, dystonic cramps, dysphagia, dysarthria, together with impaired balance and coordination, which drive the progression of the motor disability over time [4]. This progressive neurodegenerative disorder also occurs with a broad spectrum of non-motor symptoms including autonomic dysfunction (hyposmia, constipation, urinary dysfunction, orthostatic hypotension), cognitive/neurobehavioral disorders (memory loss, depression, dementia pain), and sensory and sleep abnormalities [5]. The causative factors at the basis of PD are still under study since it is a pathologically and genetically heterogeneous disease with a multifactorial etiology based on environmental and genetic factors. More than 100 genes have been identified to be involved in PD etiology [6], but most of them are likely to explain only a small portion of PD and about 90% of cases are related to the sporadic form [7]. Considering the previously reported studies, there is an active field of research around PD to better clarify the molecular mechanisms underneath the disease and eventually improve current therapeutic strategies [8]. In this context, the use of PD samples deposited on publicly available datasets could be a useful resource to expand the current knowledge about this neurodegenerative disorder [9]. Nowadays, the large amount of deposited whole blood and peripheral blood mononuclear cell expression data of PD patients have raised great interest in the analysis of the transcription profiles of peripheral tissues. Several studies on expression profiles of PD patients revealed common pathways altered in both brain and blood samples, suggesting the latter as a valid tissue to evaluate the state of the disorder [10]. Indeed, this moved the interest for the investigation of the transcriptome of PD at the peripheral blood level [11], which can represent a valid and non-invasive approach for not only clarifying the disease-leading mechanism, but also for diagnosis and prognosis purposes and monitoring therapeutic responses [12]. Multiple pathways and biological processes have been found to be altered in PD and have been proposed to play a key role underneath the disease. This includes proteostasis [13], whose alteration is often linked to various neurodegenerative diseases. In recent years, numerous studies have been conducted in this direction to elucidate the actual interplay between proteostasis and PD [14,15]. To gain deeper insights into this process, the current study aimed to investigate the dysregulation of elements involved in the proteostasis-related mechanism at the transcriptomic level. This exploration was conducted using a sample of blood tissues from PD patients sourced from a publicly available dataset deposited in the NCBI Gene Expression Omnibus (GEO) database [16].

2. Results

In the present study, we analyzed gene expression data from blood samples of healthy controls (HC) and PD patients using the rank product method to identify differentially expressed genes (DEGs). Starting from the complete dataset that included 15,786 probes, our analysis revealed a total of 375 DEGs, among which 281 identified as down-regulated and 94 as up-regulated. The terms up-regulated refer to the DEGs that have a statistically confirmed increase in the PD sample, wherein the down-regulated DEGs are less expressed in PD compared to HC. The statistical validity of our DEGs was confirmed by adjusting p values using the false discovery rate (FDR) method; this approach allowed us to carefully control the false positive rate. DEGs with an FDR-adjusted p value (q value) < 0.05 were considered statistically significant. In order to attain the most comprehensive perspective possible and to understand potential alterations across entire gene ontologies, no filters based on fold change were applied. The results about the number of up-regulated, down-regulated, and non-significative altered genes are shown in Figure 1.
Scores related to the fold change, p value, and q value of each gene are available in Supplementary Table S1. A comprehensive analysis was undertaken by uploading the entire set of DEGs onto the STRING database to elucidate potential interactions among our identified DEGs. Advanced settings were configured using the highest confidence level, 0.900, and the maximum allowable FDR, 1%. These parameters were selected to ensure a robust depiction of the network, emphasizing connections between nodes that could be deemed highly reliable. By opting for the highest confidence threshold and the strictest FDR stringency, we aimed to attain a secure and trustworthy representation of the network topology, thereby enhancing the credibility of identified interactions. This approach was undertaken to provide a comprehensive and reliable overview of potential relationships within the biological system under investigation. The result is a network composed of 371 nodes and 453 edges. Analysis of the network allowed us to obtain, for each node, information about different parameters. All the information about the nodes of our network can be consulted in the Supplementary Table S2. We used the score of degree to filter the DEGs of our network and isolate those ones with the highest levels. Score of degree indicates the number of edges that connect a node to other nodes; higher degrees correspond to more connected nodes among the network. In Figure 2 is shown the histogram with the distribution of the number of DEGs for each score of degree.
As depicted in Figure 2, it is evident that the majority of DEGs exhibit a degree score near 0. Specifically, the number of isolated DEGs, denoting those with no connections and possessing a degree score of 0, amounts to 207. Notably, an intriguing observation pertains to the group of DEGs related to the higher degree scores. In this context, 21 DEGs fall within a degree value between 18 and 24. Considering the relevant score of degree, these 21 DEGs underwent further scrutiny to assess their potential involvement in specific subnetworks. The 21 DEGs with a degree ≥ 18 are “GSPT1”, “RPL15”, “RPL18A”, “RPL29” “RPL3”, “RPL4”, “RPL5”, “RPL7”, “RPL9”, “RPS15A”, “RPS16”, “RPS17”, “RPS21”, “RPS24”, “RPS26”, “RPS3A”, “RPS4X”, “RPS4Y1”, “RPS5”, “RPS7”, “SERBP1”. Filtering these DEGs in the network, we discover that all of them are connected to each other through 207 edges with an average node degree of 19.7 and an average local clustering coefficient of 0.987. We diffused the subnetwork that included the 21 DEGs to better observe their context and we obtained a network composed of a total of 36 DEGs connected by 266 edges. This last network is shown in Figure 3.
Among the 36 DEGs reported in Figure 3, two were up-regulated and 34 down-regulated. The observed disparity in the number of up and down-regulated genes may indicate a generalized down-regulation linked to processes in which these DEGs could be involved. Prior to exploring the ontologies most pertinent to these DEGs, we assessed their log2 fold change levels, as depicted in Figure 4.
In order to ascertain the predominant gene ontologies (GO) within all the 375 DEGs in our network, we conducted an enrichment analysis for biological process (BP), cellular component (CC), and molecular function (MF). The genes associated with the enriched ontologies were cross-referenced with the 36 DEGs identified in the subnetwork illustrated in Figure 3. In pursuit of a comprehensive investigation, our analysis aimed to scrutinize transcriptional imbalances not only at the level of individual genes but also within the broader context of ontological categories. The enrichment analysis performed using all the 375 DEGs resulted in 276 BP, 63 CC, and 45 MF. We included as enriched all the GO with a q value < 0.05. All the enriched GO with associated DEGs can be retrieved in Supplementary Table S3. In Table 1 we reported all the enriched GO containing at least half of the DEGs from our focal subnetwork showcased in Figure 3.

3. Discussion

The multifactorial nature of PD makes it difficult to pinpoint the genetic background in action in the development of the disease and/or evaluate PD progression. In the present study, we tried to infer the dynamic changes in gene expression using a subset of subjects from the GSE6613 dataset [17,18]. Our aim was to clarify the processes likely altered in this disease and, possibly, suggest potential processes to be targeted to ameliorate PD symptoms or to be used as markers of PD progression. According to our results, 375 genes had a significant different expression in the blood tissues of PD patients compared to controls. The network analyses further identified a cluster of highly connected genes, as shown in Figure 3, which likely concur on the same/similar function. Interestingly, our results consistently point to a generalized impairment of the translational machinery and proteostasis.
As a whole, proteostasis is fundamental for cell viability, maintenance, and response to stressors [19,20,21,22,23,24]. The deregulation of such process can be considered a physiological event observed with aging; a growing body of evidence also highlights that this decline is dramatically accelerated in neurodegenerative diseases, including in models of PD [13,21,25,26,27]. With “proteostasis” we identify the processes related to protein synthesis for accurate protein biogenesis, trafficking, folding, and removal [13]. A balanced proteostasis is essential for neuronal activity as 33–35% of newly synthesized proteins are prone to misfolding and aberrant protein accumulation was proven lethal for neurons [28,29]. Its dysregulation in the central neurons is apparent in PD pathology along with altered signaling events and mechanisms related to translational machinery [28,30,31]. Multiple factors are likely interconnected to preserve this balance, including regulation of protein synthesis, chaperone-assisted protein folding, and protein degradation pathways [32,33].
Our analysis highlighted a group of dysregulated genes mainly involved in the protein synthesis processes, including translation initiation and regulation, ribosome structure, and ribosome components nuclear export (please refer to Supplementary Table S4). In detail, we observed a generalized down-regulation of genes whose products directly participate in the phases of translation (initiation phase, termination phase, and regulation of protein production). In particular, the eukaryotic initiation factors constitute the core element for translation initiation (EIF1AX, EIF3A, EIF4A2, EIF5). The function of these DEGs is finely regulated on multiple levels (isoforms interactions and cleaving), giving the related complexes a key role not only in translation but also in its regulation [34,35]. While no direct associations with PD were described in the literature, the down-regulation of these elements has been linked with the suppression of protein synthesis [36]. We observed the same alteration in our PD sample. The termination and protein trafficking also seem to be influenced. Our results evidenced a down-regulation of SRP, which encodes for a protein involved in the first steps of neo-synthesized proteins transport and correct localization to the endothelial reticulum [37,38,39]. A reduction in the levels of this element may lead to protein accumulation. Hernandez et al. observed that this element is involved in the regulation of αSyn biogenesis [40].
SERBP1, which encodes for a protein whose main role is to bind ribosomes to stabilize their structure and protect them from degradation, was found to be down-regulated in our sample. The ribosome is a cellular structure responsible for protein synthesis, and any disruptions in its components or function can impact the synthesis of proteins crucial for cellular health [21,41]. The down-regulation of SERBP1 likely promotes the degradation of ribosomes, which may ultimately lead to a reduction in protein synthesis [42].
The consequences of SRBP1 down-regulation are further potentiated by the reduction in the genes encoding for ribosome components’ expression (RPL15, RPL18A, RPL29, RPL3, RPL4, RPL5, RPL7, RPL9, RPS15A, RPS16, RPS17, RPS21, RPS24, RPS26, RPS3A, RPS4X, RPS5, RPS7). It should be noted that many of the components we found down-regulated have functions beyond protein biosynthesis that have been suggested in the literature; in particular, elements like RPS15A, RPL3, and RPL7 are involved in the regulation of cell proliferation and cell death [43,44,45,46]. Their altered regulation may contribute to the dopaminergic cell degeneration commonly observed during PD pathogenesis [21,41].
A similar correlation with cell death was found for NPM1, RAN, and XPO1, which are mainly involved with nucleocytoplasmic transport; for HMGN4, KDM5D, and RBBP4, which participate in chromatin-modifying complexes; and for CCND2, CDKN1B, and DDX3X, which are mainly related to cell cycle control.
Regarding nucleocytoplasmic transport, this process comprises the multistep assembly, transport, and disassembly of protein and RNA cargoes entering and exiting nuclear pores. The literature data have implicated dysfunction in nucleocytoplasmic transport in various neurodegenerative diseases, including PD [47,48]. This dysfunction causes the dysregulation in nucleus-cytoplasm protein partitioning and influences proteostasis balance, ultimately promoting intracellular deposits of proteins and cargoes, another hallmark typical of neurodegeneration [49,50,51]. This accumulation is known to be toxic for cells and may represent a trigger for cell death and neurodegeneration [52]. Additionally, nuclear transport elements (NTE), in particular NPM1, which encodes for an NTE that acts as a chaperone, has the function of suppressing the aggregation of misfolded proteins [53]. The down-regulation we found in our sample may promote an increase of pathological proteins aggregation.
Regarding the chromatin-modifying related DEGs evidenced by our analysis, they play a crucial role in the regulation of gene expression, slightly changing the overall accessibility of DNA from the replicative machinery, which ultimately leads to altered levels of proteins synthesis. Although HMGN4, KDM5D, and RBBP4 have not been directly correlated with PD, the dysregulation of chromatin-modifying complexes has been implicated in various neurodegenerative diseases, including PD [13,54,55,56]. It should be noted that, according to the literature data, they seemingly act in regulating the balance between pro-apoptotic and cell survival signals. Data on these elements mainly refer to their role in carcinogenesis [57,58,59]. Nevertheless, some hypotheses can be drafted from this data: the down-regulation of HMGN4 and RBBP4, as well as the up-regulation of KDM5D tend toward pro-apoptotic signaling. Protein encoded by KDM5D seems to act on E2F1 through demethylation, thereby suppressing cell growth and metastasis in some form of colorectal cancer [57]. Also, its suppression hampers cell apoptosis [60]. The overexpression we found in our data may suggest its potential role in triggering apoptotic signals in neuronal cells. Kugler et al. observed that increased expression of HMGN4 leads to a decrease in the levels of various tumor suppressors [58]. Consequently, reducing its expression may trigger an opposite effect, promoting cell cycle arrest and apoptotic signaling. Finally, it has been observed that RBBP4 loss disrupts neural progenitor cell cycle regulation and leads to Tp53 acetylation and apoptosis [59].
The DEGs related to cell cycle control likely exert their influence on multiple processes, including cell survival: CCND2 encodes for a cyclin2, whose main function, in complex with CDK4 or CDK6, is to regulate and control the G1/S transition during the cell cycle. CCND2 down-regulation has been associated with the suppression of cell proliferation [61] and pro-apoptotic mechanics [62,63]. The protein encoded by CDKN1B prevents the activation of cyclin E-CDK2 or cyclin D-CDK4 complexes, and thus controls the cell cycle progression at G1. DDX3X encodes for a member of the large DEAD-box protein family and its roles encompass transcriptional regulation, mRNA assembly, splicing, and export. Chen et al. evidenced how depletion of this element impairs neurite development [64] and, in addition, its down-regulation promotes apoptosis [65,66]. To date, the mechanism of action of these proteins in protein synthesis is not known. Further and more specific investigations are needed as no direct associations are known between them and the inhibition of protein synthesis. However, it should be noted that the promotion of cell cycle arrest and apoptosis usually correlates with a reduction in protein synthesis [67].
The global reduction in protein synthesis is not merely a consequence but potentially a conservative response of cells to stressors [24]. The complexity of this regulation involves a wide range of elements, with the literature data suggesting a significant contribution of PD-related genes to the overall control of protein synthesis. It has been observed that various genes related to PD, both dominant and recessive, such as EIF4G, SNCA, LRRK2, PRKN, PINK1, and PARK7, are connected to protein translation. The most apparent link is with EIF4G1. The protein encoded by EIF4G1 works as a translation initiation factor and acts as a scaffold in the eIF4F translation initiation complex, recruiting ribosomes and tRNAs to mRNA, thereby playing a crucial role in protein synthesis [68]. Its alterations have been linked to PD [69] and play a role in the deregulation of protein translation [70]. EIF4G1 is known to interact with EIF1, which was found down-regulated in our analysis, suggesting a possible common pathway of action. Regarding SNCA, it encodes for α-synuclein, a protein usually observed to form aggregates in the brains of PD patients. It has been shown that this protein may induce transcriptional repression through its interaction with ribosomal RNA and ribosomal proteins, leading to disruptions of ribosome function [71,72]. These interactions may contribute to the impairment of protein synthesis and the accumulation of misfolded proteins. Another notable player is the leucine-rich repeat kinase 2 gene (LRRK2), which is associated with both autosomal-dominant PD and the familial form of the disease [73,74]. The protein encoded by LRRK2 has recently been show to play an active role in the global decrease in protein synthesis [75]. LRRK2 may intervene in regulating protein expression via multiple mechanisms, it has been observed that this protein interacts with several translation proteins, including many translation initiation factors [76] and ribosomal proteins [77], actively regulating protein synthesis. In addition, LRRK2 appears to regulate miRNA machinery [78], through which it further regulates protein expression levels by the inhibition of translation or the degradation of mRNA [79,80]. Another element likely associated with impaired protein synthesis is the protein encoded by PINK1. This protein acts in concert with the protein encoded by PRKN to regulate mitophagy [81,82]. It has been observed that PINK1 is involved in the translation switch from cap-dependent to cap-independent translation during hypoxia, showcasing its role in adapting translation processes to cellular stress [83,84]. Finally, PARK7 (or DJ-1) is a redox-sensitive chaperone protein. When mutated, this protein has been shown to alter its affinity with ribosomal proteins, leading to decreased translation capacity [85]. In summary, our results support the presence of a global reduction in protein synthesis in PD and proteostasis imbalances, which likely translates in synaptic proteins loss, which in turn leads to eventual synaptic failure [41]. Also, the impairment of protein synthesis has been linked with the function of many PD genes, as described before. Many of these elements appear to play a regulative role in cell survival, as they are able to influence apoptotic pathways, and may take part in the pathogenesis of PD. In Figure 5 are reported all the DEGs presented with relative functions.

4. Materials and Methods

4.1. Microarray Dataset Selection

A systematic search has been performed on publicly available databases GEO and ArrayExpress [86] to select a suitable dataset. Keywords used for the search were “Parkinson” and “Blood tissue”. All datasets were then filtered to retain only the ones containing transcriptomic data from microarray or RNAseq sources. To be selected, the datasets should include PD samples and control samples. Additionally, the dataset should comprise a sufficiently large number of samples to enhance the statistical significance of our analysis.

4.2. Dataset Description

Based on the previous section filter, we selected the data samples from the publicly available dataset deposited by Scherzer et al. [17,18] on GEO database with accession number GSE6613. In detail, we selected a subset of 72 samples, 50 PD patients, and 22 age-matched HC. According to the original authors, PD patients were eligible for the inclusion in the dataset if they met the criteria for PD according to the modified United Kingdom Parkinson’s Disease Society Brain Bank [87] clinical diagnostic criteria. To note, the United Kingdom Parkinson’s Disease Society Brain Bank criteria were changed to require at least three cardinal features, three supportive features, and none of the 16 exclusion criteria. Patients showed a mean age of 69.4 ± 8.4 years, with a 78% male ratio. Also, patients included in the dataset did not take any therapeutic treatment for PD. Clinical data are reported in Supplementary Table S5. HC were enrolled from the Partners Parkinson’s and Movement Disorders Center and the Memory Disorders Unit at Massachusetts General Hospital. HC had no personal or family history of neurodegenerative diseases. Mean age of controls was 64.4 ± 10.7 years with a 50% male ratio. Further clinical data are reported in Supplementary Table S5. Exclusion criteria for all study subjects were age < 21 years, hematologic malignancies or coagulopathies, known severe anemia (hematocrit < 30), and known pregnancy. Regarding the sample processing, it has been stated that the procedures outlined below have been reported. The PAXgene Blood RNA Kit (QIAGEN, Germantown, MD, USA) was used for the purification of total RNA from whole blood samples collected in a PAXgene Blood RNA Tube and total RNA was extracted using the PAXgene procedure; four μg of total RNA was eventually reverse-transcribed into cDNA, biotin-labeled, purified, hybridized to the Affymetrix Human Genome U133A Array (GPL96) (Santa Clara, CA, USA), and scanned using standards Affymetrix procedures, for the final production of raw CEL files deposited on the database and used in our analysis. For further details on the cohort data, please refer to the author’s original work [17].

4.3. Bioinformatics Analysis

The pipeline implemented in BioTEA v1.1.0 [88] was run to download, preprocess, and perform Differential gene Expression Analysis (DEAs) with microarray data from the selected GSE6613 dataset. The first steps of the pipeline were composed by the data retrieval of raw CEL files and data preprocessing, which included background-subtraction, quantile-quantile normalization, and quality-checks on top of raw expression data. The analysis commenced with the initiation of the differential expression analysis and the establishment of specific parameters. The variables of interest, delineated by the disease state, and the targeted contrast, comparing PD samples to HC, were defined. To ensure compatibility with data generated using the Affymetrix Human Genome U133A Array (GPL96), a log2 expression threshold of 4 was applied. Genes meeting the criterion of a q value < 0.05 were retained, and their regulation status was determined based on whether the log2 fold change was positive or negative. The options and the expression matrix were eventually read and DEA was performed with the “RankProd” package, which was used to handle DEGs detection under the two experimental conditions [89,90]. Our choice to use this non-parametric statistical technique was moved by previous results obtained using rank products that have already been shown to be more biologically meaningful than the one generated by linear models, especially in studies with noisy data and a low number of replicates. The STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) database [91] was employed to construct the interaction network among DEGs in our study. STRING provides a comprehensive resource for evaluating protein–protein interactions, incorporating both experimental and predicted interaction data. For the analysis, we used as an advanced setting the highest confidence (0.900) for “Required score”, and for “FDR stringency” we used the highest (1%). We used the version 12.0 of STRING accessed on 21 December 2023 through the address https://string-db.org/. The DEGs interactions retrieved using STRING were explored through Cytoscape v. 3.10.1 [92], which allowed us to perform the analysis of the network and obtain data about different network parameters. In addition to network analysis, the list of DEGs was also employed for conducting GO enrichment analysis. This analysis identified gene ontologies that exhibited statistically significant alterations. Enrichment analysis was performed using R v. 4.2.0 (R Core Team) with the package “clusterProfiler” v. 4.6.2 [93] and “biomaRt” v. 2.54.1 [94]. The list of enriched ontologies provided insights into the specific areas where transcriptomic changes were concentrated and which processes were most affected by PD.

5. Conclusions

Our findings support the presence of a global reduction in protein synthesis in PD and proteostasis imbalances. The precise mechanisms that trigger the deregulation of protein synthesis have yet to be fully elucidated. We hypothesize that the repression of protein synthesis may be involved in neurodegeneration through a reduction in synaptic proteins. These effects are in line with what is commonly observed in PD, particularly in dopaminergic neurons, which are particularly vulnerable to disruptions in protein synthesis and cellular stress. Our results highlight protein translation defects in PD patients that are also found in blood. The precise regulation of these events may provide insights into the molecular basis of PD. The dysregulation of proteostasis can be considered a pioneering step in the pathogenesis of PD, and the screening of proteostasis-related biomarkers can be explored for early diagnosis.

Limits and Future Prospects

Our study provides insights into the possible role of translational processes and their reduction in PD patients. It is important to highlight that our analysis was carried out on peripheral blood data. Therefore, future studies correlating the findings with brain tissue would be necessary. Furthermore, extending our analysis to a larger sample would help improve the robustness and generalizability of our findings. By increasing the patient cohort and clinical information, future studies will allow us to better understand the molecular mechanisms underlying the disease.
Our study has provided preliminary data that pave the way for further studies that include other omics data (such as proteomics and metabolomics), and which will allow us to improve the early diagnosis and management framework of PD.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25021299/s1.

Author Contributions

Conceptualization, E.M.; methodology, S.D., M.C. and M.L.; software, M.C., S.D. and M.L.; formal analysis, S.D.; investigation, M.C.; data curation, M.L., S.D. and M.C.; writing—original draft preparation, M.C. and M.L.; writing—review and editing, E.M. and S.D.; visualization, S.D. and M.C.; supervision, E.M.; project administration, E.M. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by Current Research Funds 2023, Ministry of Health, Italy.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available in the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/gds; accessed on 27 November 2023), under the accession number GSE6613.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of this study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Ou, Z.; Pan, J.; Tang, S.; Duan, D.; Yu, D.; Nong, H.; Wang, Z. Global Trends in the Incidence, Prevalence, and Years Lived With Disability of Parkinson’s Disease in 204 Countries/Territories From 1990 to 2019. Front. Public Health 2021, 9, 776847. [Google Scholar] [CrossRef] [PubMed]
  2. Jankovic, J. Parkinson’s disease: Clinical features and diagnosis. J. Neurol. Neurosurg. Psychiatry 2008, 79, 368–376. [Google Scholar] [CrossRef] [PubMed]
  3. Lees, A. The bare essentials: Parkinson’s disease. Pract. Neurol. 2010, 10, 240–246. [Google Scholar] [CrossRef] [PubMed]
  4. Magrinelli, F.; Picelli, A.; Tocco, P.; Federico, A.; Roncari, L.; Smania, N.; Zanette, G.; Tamburin, S. Pathophysiology of Motor Dysfunction in Parkinson’s Disease as the Rationale for Drug Treatment and Rehabilitation. Park. Dis. 2016, 2016, 9832839. [Google Scholar] [CrossRef] [PubMed]
  5. Tolosa, E.; Garrido, A.; Scholz, S.W.; Poewe, W. Challenges in the diagnosis of Parkinson’s disease. Lancet Neurol. 2021, 20, 385–397. [Google Scholar] [CrossRef] [PubMed]
  6. Ye, H.; Robak, L.A.; Yu, M.; Cykowski, M.; Shulman, J.M. Genetics and Pathogenesis of Parkinson’s Syndrome. Annu. Rev. Pathol. 2023, 18, 95–121. [Google Scholar] [CrossRef] [PubMed]
  7. de Lau, L.M.; Breteler, M.M. Epidemiology of Parkinson’s disease. Lancet Neurol. 2006, 5, 525–535. [Google Scholar] [CrossRef]
  8. Ibrahim, R.; Pasic, M.; Yousef, G.M. Omics for personalized medicine: Defining the current we swim in. Expert Rev. Mol. Diagn. 2016, 16, 719–722. [Google Scholar] [CrossRef]
  9. Scholz, S.W.; Mhyre, T.; Ressom, H.; Shah, S.; Federoff, H.J. Genomics and bioinformatics of Parkinson’s disease. Cold Spring Harb. Perspect. Med. 2012, 2, a009449. [Google Scholar] [CrossRef]
  10. Moni, M.A.; Rana, H.K.; Islam, M.B.; Ahmed, M.B.; Xu, H.; Hasan, M.A.M.; Lei, Y.; Quinn, J.M.W. A computational approach to identify blood cell-expressed Parkinson’s disease biomarkers that are coordinately expressed in brain tissue. Comput. Biol. Med. 2019, 113, 103385. [Google Scholar] [CrossRef]
  11. Kurvits, L.; Lattekivi, F.; Reimann, E.; Kadastik-Eerme, L.; Kasterpalu, K.M.; Koks, S.; Taba, P.; Planken, A. Transcriptomic profiles in Parkinson’s disease. Exp. Biol. Med. 2021, 246, 584–595. [Google Scholar] [CrossRef] [PubMed]
  12. Irmady, K.; Hale, C.R.; Qadri, R.; Fak, J.; Simelane, S.; Carroll, T.; Przedborski, S.; Darnell, R.B. Blood transcriptomic signatures associated with molecular changes in the brain and clinical outcomes in Parkinson’s disease. Nat. Commun. 2023, 14, 3956. [Google Scholar] [CrossRef] [PubMed]
  13. Kulkarni, A.; Preeti, K.; Tryphena, K.P.; Srivastava, S.; Singh, S.B.; Khatri, D.K. Proteostasis in Parkinson’s disease: Recent development and possible implication in diagnosis and therapeutics. Ageing Res. Rev. 2023, 84, 101816. [Google Scholar] [CrossRef] [PubMed]
  14. Perez-Carrion, M.D.; Posadas, I.; Solera, J.; Cena, V. LRRK2 and Proteostasis in Parkinson’s Disease. Int. J. Mol. Sci. 2022, 23, 6808. [Google Scholar] [CrossRef]
  15. Bosco, D.A.; LaVoie, M.J.; Petsko, G.A.; Ringe, D. Proteostasis and movement disorders: Parkinson’s disease and amyotrophic lateral sclerosis. Cold Spring Harb. Perspect. Biol. 2011, 3, a007500. [Google Scholar] [CrossRef] [PubMed]
  16. Edgar, R.; Domrachev, M.; Lash, A.E. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30, 207–210. [Google Scholar] [CrossRef] [PubMed]
  17. Scherzer, C.R.; Eklund, A.C.; Morse, L.J.; Liao, Z.; Locascio, J.J.; Fefer, D.; Schwarzschild, M.A.; Schlossmacher, M.G.; Hauser, M.A.; Vance, J.M.; et al. Molecular markers of early Parkinson’s disease based on gene expression in blood. Proc. Natl. Acad. Sci. USA 2007, 104, 955–960. [Google Scholar] [CrossRef]
  18. Scherzer, C.R.; Grass, J.A.; Liao, Z.; Pepivani, I.; Zheng, B.; Eklund, A.C.; Ney, P.A.; Ng, J.; McGoldrick, M.; Mollenhauer, B.; et al. GATA transcription factors directly regulate the Parkinson’s disease-linked gene alpha-synuclein. Proc. Natl. Acad. Sci. USA 2008, 105, 10907–10912. [Google Scholar] [CrossRef]
  19. Taylor, R.C.; Dillin, A. Aging as an event of proteostasis collapse. Cold Spring Harb. Perspect. Biol. 2011, 3, a004440. [Google Scholar] [CrossRef]
  20. Chen, B.; Retzlaff, M.; Roos, T.; Frydman, J. Cellular strategies of protein quality control. Cold Spring Harb. Perspect. Biol. 2011, 3, a004374. [Google Scholar] [CrossRef]
  21. Martin, I. Decoding Parkinson’s Disease Pathogenesis: The Role of Deregulated mRNA Translation. J. Park. Dis. 2016, 6, 17–27. [Google Scholar] [CrossRef]
  22. Somers, J.; Wilson, L.A.; Kilday, J.P.; Horvilleur, E.; Cannell, I.G.; Poyry, T.A.; Cobbold, L.C.; Kondrashov, A.; Knight, J.R.; Puget, S.; et al. A common polymorphism in the 5′ UTR of ERCC5 creates an upstream ORF that confers resistance to platinum-based chemotherapy. Genes Dev. 2015, 29, 1891–1896. [Google Scholar] [CrossRef] [PubMed]
  23. King, H.A.; Cobbold, L.C.; Pichon, X.; Poyry, T.; Wilson, L.A.; Booden, H.; Jukes-Jones, R.; Cain, K.; Lilley, K.S.; Bushell, M.; et al. Remodelling of a polypyrimidine tract-binding protein complex during apoptosis activates cellular IRESs. Cell Death Differ. 2014, 21, 161–171. [Google Scholar] [CrossRef] [PubMed]
  24. Powley, I.R.; Kondrashov, A.; Young, L.A.; Dobbyn, H.C.; Hill, K.; Cannell, I.G.; Stoneley, M.; Kong, Y.W.; Cotes, J.A.; Smith, G.C.; et al. Translational reprogramming following UVB irradiation is mediated by DNA-PKcs and allows selective recruitment to the polysomes of mRNAs encoding DNA repair enzymes. Genes Dev. 2009, 23, 1207–1220. [Google Scholar] [CrossRef] [PubMed]
  25. Lehtonen, S.; Sonninen, T.M.; Wojciechowski, S.; Goldsteins, G.; Koistinaho, J. Dysfunction of Cellular Proteostasis in Parkinson’s Disease. Front. Neurosci. 2019, 13, 457. [Google Scholar] [CrossRef] [PubMed]
  26. Bi, M.; Du, X.; Jiao, Q.; Chen, X.; Jiang, H. Expanding the role of proteasome homeostasis in Parkinson’s disease: Beyond protein breakdown. Cell Death Dis. 2021, 12, 154. [Google Scholar] [CrossRef] [PubMed]
  27. Yang, L.; Cao, Y.; Zhao, J.; Fang, Y.; Liu, N.; Zhang, Y. Multidimensional Proteomics Identifies Declines in Protein Homeostasis and Mitochondria as Early Signals for Normal Aging and Age-associated Disease in Drosophila. Mol. Cell. Proteom. MCP 2019, 18, 2078–2088. [Google Scholar] [CrossRef]
  28. Amrutha, K.; Mishra, A.; Singh, S. Implications of intracellular protein degradation pathways in Parkinson’s disease and therapeutics. J. Neurosci. Res. 2022, 100, 1834–1844. [Google Scholar] [CrossRef]
  29. Mymrikov, E.V.; Daake, M.; Richter, B.; Haslbeck, M.; Buchner, J. The Chaperone Activity and Substrate Spectrum of Human Small Heat Shock Proteins. J. Biol. Chem. 2017, 292, 672–684. [Google Scholar] [CrossRef]
  30. Darios, F.; Stevanin, G. Impairment of Lysosome Function and Autophagy in Rare Neurodegenerative Diseases. J. Mol. Biol. 2020, 432, 2714–2734. [Google Scholar] [CrossRef]
  31. Kip, E.; Parr-Brownlie, L.C. Reducing neuroinflammation via therapeutic compounds and lifestyle to prevent or delay progression of Parkinson’s disease. Ageing Res. Rev. 2022, 78, 101618. [Google Scholar] [CrossRef] [PubMed]
  32. Lazaro, D.F.; Outeiro, T.F. The Interplay Between Proteostasis Systems and Parkinson’s Disease. Adv. Exp. Med. Biol. 2020, 1233, 223–236. [Google Scholar] [CrossRef] [PubMed]
  33. Lim, J.; Yue, Z. Neuronal aggregates: Formation, clearance, and spreading. Dev. Cell 2015, 32, 491–501. [Google Scholar] [CrossRef] [PubMed]
  34. Kumar Deshmukh, F.; Yaffe, D.; Olshina, M.A.; Ben-Nissan, G.; Sharon, M. The Contribution of the 20S Proteasome to Proteostasis. Biomolecules 2019, 9, 190. [Google Scholar] [CrossRef] [PubMed]
  35. Hoshino, S.; Imai, M.; Kobayashi, T.; Uchida, N.; Katada, T. The eukaryotic polypeptide chain releasing factor (eRF3/GSPT) carrying the translation termination signal to the 3’-Poly(A) tail of mRNA. Direct association of erf3/GSPT with polyadenylate-binding protein. J. Biol. Chem. 1999, 274, 16677–16680. [Google Scholar] [CrossRef] [PubMed]
  36. Sehrawat, U.; Koning, F.; Ashkenazi, S.; Stelzer, G.; Leshkowitz, D.; Dikstein, R. Cancer-Associated Eukaryotic Translation Initiation Factor 1A Mutants Impair Rps3 and Rps10 Binding and Enhance Scanning of Cell Cycle Genes. Mol. Cell. Biol. 2019, 39, e00441-18. [Google Scholar] [CrossRef] [PubMed]
  37. Akopian, D.; Shen, K.; Zhang, X.; Shan, S.O. Signal recognition particle: An essential protein-targeting machine. Annu. Rev. Biochem. 2013, 82, 693–721. [Google Scholar] [CrossRef]
  38. Kellogg, M.K.; Miller, S.C.; Tikhonova, E.B.; Karamyshev, A.L. SRPassing Co-translational Targeting: The Role of the Signal Recognition Particle in Protein Targeting and mRNA Protection. Int. J. Mol. Sci. 2021, 22, 6284. [Google Scholar] [CrossRef]
  39. Mary, C.; Scherrer, A.; Huck, L.; Lakkaraju, A.K.; Thomas, Y.; Johnson, A.E.; Strub, K. Residues in SRP9/14 essential for elongation arrest activity of the signal recognition particle define a positively charged functional domain on one side of the protein. RNA 2010, 16, 969–979. [Google Scholar] [CrossRef]
  40. Hernandez, S.M.; Tikhonova, E.B.; Baca, K.R.; Zhao, F.; Zhu, X.; Karamyshev, A.L. Unexpected Implication of SRP and AGO2 in Parkinson’s Disease: Involvement in Alpha-Synuclein Biogenesis. Cells 2021, 10, 2792. [Google Scholar] [CrossRef]
  41. Moreno, J.A.; Radford, H.; Peretti, D.; Steinert, J.R.; Verity, N.; Martin, M.G.; Halliday, M.; Morgan, J.; Dinsdale, D.; Ortori, C.A.; et al. Sustained translational repression by eIF2alpha-P mediates prion neurodegeneration. Nature 2012, 485, 507–511. [Google Scholar] [CrossRef] [PubMed]
  42. Shetty, S.; Hofstetter, J.; Battaglioni, S.; Ritz, D.; Hall, M.N. TORC1 phosphorylates and inhibits the ribosome preservation factor Stm1 to activate dormant ribosomes. EMBO J. 2023, 42, e112344. [Google Scholar] [CrossRef] [PubMed]
  43. Feng, W.; Liang, C.; Wang, C.; Yu, X.; Li, Q.; Yang, H. Knockdown of ribosomal protein S15A inhibits proliferation of breast cancer cells through induction of apoptosis in vitro. Cytotechnology 2018, 70, 1315–1323. [Google Scholar] [CrossRef] [PubMed]
  44. Dai, M.S.; Shi, D.; Jin, Y.; Sun, X.X.; Zhang, Y.; Grossman, S.R.; Lu, H. Regulation of the MDM2-p53 pathway by ribosomal protein L11 involves a post-ubiquitination mechanism. J. Biol. Chem. 2006, 281, 24304–24313. [Google Scholar] [CrossRef]
  45. Dai, M.S.; Challagundla, K.B.; Sun, X.X.; Palam, L.R.; Zeng, S.X.; Wek, R.C.; Lu, H. Physical and functional interaction between ribosomal protein L11 and the tumor suppressor ARF. J. Biol. Chem. 2012, 287, 17120–17129. [Google Scholar] [CrossRef]
  46. Shi, Y.; Zhai, H.; Wang, X.; Han, Z.; Liu, C.; Lan, M.; Du, J.; Guo, C.; Zhang, Y.; Wu, K.; et al. Ribosomal proteins S13 and L23 promote multidrug resistance in gastric cancer cells by suppressing drug-induced apoptosis. Exp. Cell Res. 2004, 296, 337–346. [Google Scholar] [CrossRef] [PubMed]
  47. Argueti-Ostrovsky, S.; Alfahel, L.; Kahn, J.; Israelson, A. All Roads Lead to Rome: Different Molecular Players Converge to Common Toxic Pathways in Neurodegeneration. Cells 2021, 10, 2438. [Google Scholar] [CrossRef]
  48. Ma, K.L.; Song, L.K.; Yuan, Y.H.; Zhang, Y.; Han, N.; Gao, K.; Chen, N.H. The nuclear accumulation of alpha-synuclein is mediated by importin alpha and promotes neurotoxicity by accelerating the cell cycle. Neuropharmacology 2014, 82, 132–142. [Google Scholar] [CrossRef]
  49. Li, N.; Lagier-Tourenne, C. Nuclear pores: The gate to neurodegeneration. Nat. Neurosci. 2018, 21, 156–158. [Google Scholar] [CrossRef]
  50. Bitetto, G.; Di Fonzo, A. Nucleo-cytoplasmic transport defects and protein aggregates in neurodegeneration. Transl. Neurodegener. 2020, 9, 25. [Google Scholar] [CrossRef]
  51. Ferreira, P.A. Nucleocytoplasmic transport at the crossroads of proteostasis, neurodegeneration and neuroprotection. FEBS Lett. 2023, 597, 2567–2589. [Google Scholar] [CrossRef] [PubMed]
  52. Pan, Y.; Zong, Q.; Li, G.; Wu, Z.; Du, T.; Huang, Z.; Zhang, Y.; Ma, K. Nuclear localization of alpha-synuclein affects the cognitive and motor behavior of mice by inducing DNA damage and abnormal cell cycle of hippocampal neurons. Front. Mol. Neurosci. 2022, 15, 1015881. [Google Scholar] [CrossRef] [PubMed]
  53. Frottin, F.; Schueder, F.; Tiwary, S.; Gupta, R.; Korner, R.; Schlichthaerle, T.; Cox, J.; Jungmann, R.; Hartl, F.U.; Hipp, M.S. The nucleolus functions as a phase-separated protein quality control compartment. Science 2019, 365, 342–347. [Google Scholar] [CrossRef] [PubMed]
  54. Sugeno, N.; Jackel, S.; Voigt, A.; Wassouf, Z.; Schulze-Hentrich, J.; Kahle, P.J. alpha-Synuclein enhances histone H3 lysine-9 dimethylation and H3K9me2-dependent transcriptional responses. Sci. Rep. 2016, 6, 36328. [Google Scholar] [CrossRef]
  55. Ostvold, A.C.; Grundt, K.; Wiese, C. NUCKS1 is a highly modified, chromatin-associated protein involved in a diverse set of biological and pathophysiological processes. Biochem. J. 2022, 479, 1205–1220. [Google Scholar] [CrossRef]
  56. Fujioka, S.; Wszolek, Z.K. Clinical aspects of familial forms of frontotemporal dementia associated with parkinsonism. J. Mol. Neurosci. 2011, 45, 359–365. [Google Scholar] [CrossRef]
  57. Liu, M.; Gao, N. KDM5D inhibits the transcriptional activation of FKBP4 by suppressing the expression of E2F1 in colorectal cancer in males. Biochem. Pharmacol. 2021, 194, 114814. [Google Scholar] [CrossRef]
  58. Kugler, J.; Postnikov, Y.V.; Furusawa, T.; Kimura, S.; Bustin, M. Elevated HMGN4 expression potentiates thyroid tumorigenesis. Carcinogenesis 2017, 38, 391–401. [Google Scholar] [CrossRef]
  59. Schultz-Rogers, L.E.; Thayer, M.L.; Kambakam, S.; Wierson, W.A.; Helmer, J.A.; Wishman, M.D.; Wall, K.A.; Greig, J.L.; Forsman, J.L.; Puchhalapalli, K.; et al. Rbbp4 loss disrupts neural progenitor cell cycle regulation independent of Rb and leads to Tp53 acetylation and apoptosis. Dev. Dyn. 2022, 251, 1267–1290. [Google Scholar] [CrossRef]
  60. Jangravi, Z.; Tabar, M.S.; Mirzaei, M.; Parsamatin, P.; Vakilian, H.; Alikhani, M.; Shabani, M.; Haynes, P.A.; Goodchild, A.K.; Gourabi, H.; et al. Two Splice Variants of Y Chromosome-Located Lysine-Specific Demethylase 5D Have Distinct Function in Prostate Cancer Cell Line (DU-145). J. Proteome Res. 2015, 14, 3492–3502. [Google Scholar] [CrossRef]
  61. Jiang, P.; Mao, L.; Lei, X.; Luo, C.; Zhang, Y.; Zhong, X.; Yin, Z.; Xu, X.; Li, D.; Zheng, Q. miR-1297 inhibits osteosarcoma cell proliferation and growth by targeting CCND2. Am. J. Cancer Res. 2022, 12, 3464–3478. [Google Scholar] [PubMed]
  62. Yadav, S.; Pandey, A.; Shukla, A.; Talwelkar, S.S.; Kumar, A.; Pant, A.B.; Parmar, D. miR-497 and miR-302b regulate ethanol-induced neuronal cell death through BCL2 protein and cyclin D2. J. Biol. Chem. 2011, 286, 37347–37357. [Google Scholar] [CrossRef] [PubMed]
  63. Zhang, F.; Zhu, Y.; Fan, G.; Hu, S. MicroRNA-2682-3p inhibits osteosarcoma cell proliferation by targeting CCND2, MMP8 and Myd88. Oncol. Lett. 2018, 16, 3359–3364. [Google Scholar] [CrossRef] [PubMed]
  64. Chen, H.H.; Yu, H.I.; Tarn, W.Y. DDX3 Modulates Neurite Development via Translationally Activating an RNA Regulon Involved in Rac1 Activation. J. Neurosci. 2016, 36, 9792–9804. [Google Scholar] [CrossRef] [PubMed]
  65. Nguyen, C.N.; Nguyen, K.V.A.; Eom, S.; Choi, Y.J.; Kang, L.J.; Lee, J.; Kim, C.; Lee, S.; Lee, S.G.; Lee, J.H. Regulation of p21 expression for anti-apoptotic activity of DDX3 against sanguinarine-induced cell death on intrinsic pathway. Phytomedicine 2019, 65, 153096. [Google Scholar] [CrossRef] [PubMed]
  66. Szappanos, D.; Tschismarov, R.; Perlot, T.; Westermayer, S.; Fischer, K.; Platanitis, E.; Kallinger, F.; Novatchkova, M.; Lassnig, C.; Muller, M.; et al. The RNA helicase DDX3X is an essential mediator of innate antimicrobial immunity. PLoS Pathog. 2018, 14, e1007397. [Google Scholar] [CrossRef] [PubMed]
  67. Jeffrey, I.W.; Bushell, M.; Tilleray, V.J.; Morley, S.; Clemens, M.J. Inhibition of protein synthesis in apoptosis: Differential requirements by the tumor necrosis factor alpha family and a DNA-damaging agent for caspases and the double-stranded RNA-dependent protein kinase. Cancer Res. 2002, 62, 2272–2280. [Google Scholar]
  68. Sonenberg, N.; Dever, T.E. Eukaryotic translation initiation factors and regulators. Curr. Opin. Struct. Biol. 2003, 13, 56–63. [Google Scholar] [CrossRef]
  69. Chartier-Harlin, M.C.; Dachsel, J.C.; Vilarino-Guell, C.; Lincoln, S.J.; Lepretre, F.; Hulihan, M.M.; Kachergus, J.; Milnerwood, A.J.; Tapia, L.; Song, M.S.; et al. Translation initiator EIF4G1 mutations in familial Parkinson disease. Am. J. Hum. Genet. 2011, 89, 398–406. [Google Scholar] [CrossRef]
  70. Dhungel, N.; Eleuteri, S.; Li, L.B.; Kramer, N.J.; Chartron, J.W.; Spencer, B.; Kosberg, K.; Fields, J.A.; Stafa, K.; Adame, A.; et al. Parkinson’s Disease Genes VPS35 and EIF4G1 Interact Genetically and Converge on alpha-Synuclein. Neuron 2023, 111, 138. [Google Scholar] [CrossRef]
  71. Ho, D.H.; Nam, D.; Jeong, S.; Seo, M.K.; Park, S.W.; Seol, W.; Son, I. Expression of transduced nucleolin promotes the clearance of accumulated alpha-synuclein in rodent cells and animal model. Neurobiol. Dis. 2021, 154, 105349. [Google Scholar] [CrossRef] [PubMed]
  72. Ho, D.H.; Kim, H.; Nam, D.; Heo, J.; Son, I. Nuclear alpha-Synuclein-Derived Cytotoxic Effect via Altered Ribosomal RNA Processing in Primary Mouse Embryonic Fibroblasts. Int. J. Mol. Sci. 2023, 24, 2132. [Google Scholar] [CrossRef] [PubMed]
  73. Ohtonen, S.; Giudice, L.; Jantti, H.; Fazaludeen, M.F.; Shakirzyanova, A.; Gomez-Budia, M.; Valimaki, N.N.; Niskanen, J.; Korvenlaita, N.; Fagerlund, I.; et al. Human iPSC-derived microglia carrying the LRRK2-G2019S mutation show a Parkinson’s disease related transcriptional profile and function. Sci. Rep. 2023, 13, 22118. [Google Scholar] [CrossRef]
  74. Sanz Murillo, M.; Villagran Suarez, A.; Dederer, V.; Chatterjee, D.; Alegrio Louro, J.; Knapp, S.; Mathea, S.; Leschziner, A.E. Inhibition of Parkinson’s disease-related LRRK2 by type I and type II kinase inhibitors: Activity and structures. Sci. Adv. 2023, 9, eadk6191. [Google Scholar] [CrossRef] [PubMed]
  75. Deshpande, P.; Flinkman, D.; Hong, Y.; Goltseva, E.; Siino, V.; Sun, L.; Peltonen, S.; Elo, L.L.; Kaasinen, V.; James, P.; et al. Protein synthesis is suppressed in sporadic and familial Parkinson’s disease by LRRK2. FASEB J. 2020, 34, 14217–14233. [Google Scholar] [CrossRef] [PubMed]
  76. Reyniers, L.; Del Giudice, M.G.; Civiero, L.; Belluzzi, E.; Lobbestael, E.; Beilina, A.; Arrigoni, G.; Derua, R.; Waelkens, E.; Li, Y.; et al. Differential protein-protein interactions of LRRK1 and LRRK2 indicate roles in distinct cellular signaling pathways. J. Neurochem. 2014, 131, 239–250. [Google Scholar] [CrossRef] [PubMed]
  77. Martin, I.; Kim, J.W.; Lee, B.D.; Kang, H.C.; Xu, J.C.; Jia, H.; Stankowski, J.; Kim, M.S.; Zhong, J.; Kumar, M.; et al. Ribosomal protein s15 phosphorylation mediates LRRK2 neurodegeneration in Parkinson’s disease. Cell 2014, 157, 472–485. [Google Scholar] [CrossRef]
  78. Gehrke, S.; Imai, Y.; Sokol, N.; Lu, B. Pathogenic LRRK2 negatively regulates microRNA-mediated translational repression. Nature 2010, 466, 637–641. [Google Scholar] [CrossRef]
  79. Bartel, D.P. MicroRNAs: Target recognition and regulatory functions. Cell 2009, 136, 215–233. [Google Scholar] [CrossRef]
  80. Cipolla, G.A. A non-canonical landscape of the microRNA system. Front. Genet. 2014, 5, 337. [Google Scholar] [CrossRef]
  81. Narendra, D.P.; Jin, S.M.; Tanaka, A.; Suen, D.F.; Gautier, C.A.; Shen, J.; Cookson, M.R.; Youle, R.J. PINK1 is selectively stabilized on impaired mitochondria to activate Parkin. PLoS Biol. 2010, 8, e1000298. [Google Scholar] [CrossRef] [PubMed]
  82. Geisler, S.; Holmstrom, K.M.; Skujat, D.; Fiesel, F.C.; Rothfuss, O.C.; Kahle, P.J.; Springer, W. PINK1/Parkin-mediated mitophagy is dependent on VDAC1 and p62/SQSTM1. Nat. Cell Biol. 2010, 12, 119–131. [Google Scholar] [CrossRef] [PubMed]
  83. Lin, W.; Wadlington, N.L.; Chen, L.; Zhuang, X.; Brorson, J.R.; Kang, U.J. Loss of PINK1 attenuates HIF-1alpha induction by preventing 4E-BP1-dependent switch in protein translation under hypoxia. J. Neurosci. 2014, 34, 3079–3089. [Google Scholar] [CrossRef] [PubMed]
  84. Ottone, C.; Galasso, A.; Gemei, M.; Pisa, V.; Gigliotti, S.; Piccioni, F.; Graziani, F.; Verrotti di Pianella, A. Diminution of eIF4E activity suppresses parkin mutant phenotypes. Gene 2011, 470, 12–19. [Google Scholar] [CrossRef] [PubMed]
  85. Kthiri, F.; Gautier, V.; Le, H.T.; Prere, M.F.; Fayet, O.; Malki, A.; Landoulsi, A.; Richarme, G. Translational defects in a mutant deficient in YajL, the bacterial homolog of the parkinsonism-associated protein DJ-1. J. Bacteriol. 2010, 192, 6302–6306. [Google Scholar] [CrossRef] [PubMed]
  86. Parkinson, H.; Kapushesky, M.; Shojatalab, M.; Abeygunawardena, N.; Coulson, R.; Farne, A.; Holloway, E.; Kolesnykov, N.; Lilja, P.; Lukk, M.; et al. ArrayExpress—A public database of microarray experiments and gene expression profiles. Nucleic Acids Res. 2007, 35, D747–D750. [Google Scholar] [CrossRef]
  87. Hughes, A.J.; Daniel, S.E.; Kilford, L.; Lees, A.J. Accuracy of clinical diagnosis of idiopathic Parkinson’s disease: A clinico-pathological study of 100 cases. J. Neurol. Neurosurg. Psychiatry 1992, 55, 181–184. [Google Scholar] [CrossRef]
  88. Visentin, L.; Scarpellino, G.; Chinigo, G.; Munaron, L.; Ruffinatti, F.A. BioTEA: Containerized Methods of Analysis for Microarray-Based Transcriptomics Data. Biology 2022, 11, 1346. [Google Scholar] [CrossRef]
  89. Hong, F.; Breitling, R.; McEntee, C.W.; Wittner, B.S.; Nemhauser, J.L.; Chory, J. RankProd: A bioconductor package for detecting differentially expressed genes in meta-analysis. Bioinformatics 2006, 22, 2825–2827. [Google Scholar] [CrossRef]
  90. Del Carratore, F.; Jankevics, A.; Eisinga, R.; Heskes, T.; Hong, F.; Breitling, R. RankProd 2.0: A refactored bioconductor package for detecting differentially expressed features in molecular profiling datasets. Bioinformatics 2017, 33, 2774–2775. [Google Scholar] [CrossRef]
  91. Szklarczyk, D.; Gable, A.L.; Nastou, K.C.; Lyon, D.; Kirsch, R.; Pyysalo, S.; Doncheva, N.T.; Legeay, M.; Fang, T.; Bork, P.; et al. The STRING database in 2021: Customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021, 49, D605–D612. [Google Scholar] [CrossRef] [PubMed]
  92. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef] [PubMed]
  93. Yu, G.; Wang, L.G.; Han, Y.; He, Q.Y. clusterProfiler: An R package for comparing biological themes among gene clusters. Omics 2012, 16, 284–287. [Google Scholar] [CrossRef] [PubMed]
  94. Kinsella, R.J.; Kahari, A.; Haider, S.; Zamora, J.; Proctor, G.; Spudich, G.; Almeida-King, J.; Staines, D.; Derwent, P.; Kerhornou, A.; et al. Ensembl BioMarts: A hub for data retrieval across taxonomic space. Database 2011, 2011, bar030. [Google Scholar] [CrossRef]
Ijms 25 01299 g001
Figure 1. In the dot plot are reported all the features explored in the analysis. The line that intercepts y axis is related to our threshold of significance of q value (0.05), all the genes above this line are considered differentially expressed. In the x axis is reported the log2 fold change that discriminate up and down-regulated DEGs defined by red or green color respectively.
Figure 1. In the dot plot are reported all the features explored in the analysis. The line that intercepts y axis is related to our threshold of significance of q value (0.05), all the genes above this line are considered differentially expressed. In the x axis is reported the log2 fold change that discriminate up and down-regulated DEGs defined by red or green color respectively.
Ijms 25 01299 g001
Ijms 25 01299 g002
Figure 2. In the histogram are reported, on the x axis, the degree scores from the lowest (0) to the highest (24), while in the y axis is reported the number of DEGs for each degree value.
Figure 2. In the histogram are reported, on the x axis, the degree scores from the lowest (0) to the highest (24), while in the y axis is reported the number of DEGs for each degree value.
Ijms 25 01299 g002
Ijms 25 01299 g003
Figure 3. Subnetwork composed of 36 DEGs connected through 266 edges. In green are reported the down-regulated DEGs in PD, while in red are reported the up-regulated ones. The DEGs with a rectangular shape and thickened edges are the 21 with a value of degree ≥ 18.
Figure 3. Subnetwork composed of 36 DEGs connected through 266 edges. In green are reported the down-regulated DEGs in PD, while in red are reported the up-regulated ones. The DEGs with a rectangular shape and thickened edges are the 21 with a value of degree ≥ 18.
Ijms 25 01299 g003
Ijms 25 01299 g004
Figure 4. In the lollipop plot is reported the list of 36 DEGs shown in Figure 3 on the x axis and the associated log2 fold change on the y axis. Green dots are associated to down-regulated DEGs, and red dots to up-regulated DEGs.
Figure 4. In the lollipop plot is reported the list of 36 DEGs shown in Figure 3 on the x axis and the associated log2 fold change on the y axis. Green dots are associated to down-regulated DEGs, and red dots to up-regulated DEGs.
Ijms 25 01299 g004
Ijms 25 01299 g005
Figure 5. We summarized the processes linked to the DEGs we found in our analyses. In green we reported the down-regulated DEGs in PD patients, while in red we reported the up-regulated ones. (1) Chromatin Modifications (HMGN4 ↓, KDM5D ↑, RBBP4 ↓); (2) Ribosome Structure of 60S (RPL15 ↓, RPL18A ↓, RPL29 ↓, RPL3 ↓, RPL4 ↓, RPL5 ↓, RPL7 ↓, RPL9 ↓) and 40S (RPS15A ↓, RPS16 ↓, RPS17 ↓, RPS21 ↓, RPS24 ↓,RPS26 ↓, RPS3A ↓, RPS4X ↓, RPS4Y ↑, RPS5 ↓, RPS7 ↓) subunits; (3) Nuclear transport (NPM1 ↓, RAN ↓, XPO1 ↓); (4) 43S pre-initiation complex assembly and mRNA Cap Recognition (EIF1AX ↓, EIF3A ↓, EIF4A2 ↓, EIF5 ↓); (5) Ribosome Stabilization (SERBP1 ↓); (6) Signal Sequence recognition and Protein localization (SRP ↓); (7) Cell Cycle Control (CCND2 ↓, CDKN1B ↓, DDX3X ↓, RPS15A ↓, RPL3 ↓, RPL7 ↓). Parts of the figure were drawn by using pictures from Servier Medical Art. Servier Medical Art by Servier is licensed under a Creative Commons Attribution 3.0 Unported License (https://creativecommons.org/licenses/by/3.0/, 19 January 2024).
Figure 5. We summarized the processes linked to the DEGs we found in our analyses. In green we reported the down-regulated DEGs in PD patients, while in red we reported the up-regulated ones. (1) Chromatin Modifications (HMGN4 ↓, KDM5D ↑, RBBP4 ↓); (2) Ribosome Structure of 60S (RPL15 ↓, RPL18A ↓, RPL29 ↓, RPL3 ↓, RPL4 ↓, RPL5 ↓, RPL7 ↓, RPL9 ↓) and 40S (RPS15A ↓, RPS16 ↓, RPS17 ↓, RPS21 ↓, RPS24 ↓,RPS26 ↓, RPS3A ↓, RPS4X ↓, RPS4Y ↑, RPS5 ↓, RPS7 ↓) subunits; (3) Nuclear transport (NPM1 ↓, RAN ↓, XPO1 ↓); (4) 43S pre-initiation complex assembly and mRNA Cap Recognition (EIF1AX ↓, EIF3A ↓, EIF4A2 ↓, EIF5 ↓); (5) Ribosome Stabilization (SERBP1 ↓); (6) Signal Sequence recognition and Protein localization (SRP ↓); (7) Cell Cycle Control (CCND2 ↓, CDKN1B ↓, DDX3X ↓, RPS15A ↓, RPL3 ↓, RPL7 ↓). Parts of the figure were drawn by using pictures from Servier Medical Art. Servier Medical Art by Servier is licensed under a Creative Commons Attribution 3.0 Unported License (https://creativecommons.org/licenses/by/3.0/, 19 January 2024).
Ijms 25 01299 g005
Table 1. DEGs involved in enriched GO also present in our subnetwork.
Table 1. DEGs involved in enriched GO also present in our subnetwork.
OntologyDescriptionTotal DEGsSubnetwork DEGsDEGs Symbol
BPCytoplasmic
translation		26/37521/36EIF3A/EIF4A2/EIF5/RPL15/
RPL18A/RPL29/RPL3/
RPL4/RPL5/RPL7/RPL9/
RPS15A/RPS16/RPS17/
RPS21/RPS24/RPS26/RPS3A/
RPS4X/RPS5/RPS7
CCRibosomal
subunit		21/37520/36NPM1/RPL15/RPL18A/RPL29/
RPL3/RPL4/RPL5/
RPL7/RPL9/RPS15A/RPS16/
RPS17/RPS21/RPS24/
RPS26/RPS3A/RPS4X/RPS4Y1/
RPS5/RPS7
CCRibosome22/37520/36NPM1/RPL15/RPL18A/RPL29/
RPL3/RPL4/RPL5/
RPL7/RPL9/RPS15A/RPS16/
RPS17/RPS21/RPS24/
RPS26/RPS3A/RPS4X/
RPS4Y1/RPS5/RPS7
CCCytosolic
ribosome		20/37519/36RPL15/RPL18A/RPL29/
RPL3/RPL4/RPL5/RPL7/
RPL9/RPS15A/RPS16/
RPS17/RPS21/RPS24/
RPS26/RPS3A/RPS4X/
RPS4Y1/RPS5/RPS7
MFStructural constituent
of ribosome		20/37519/36RPL15/RPL18A/RPL29/
RPL3/RPL4/RPL5/RPL7/
RPL9/RPS15A/RPS16/RPS17/
RPS21/RPS24/
RPS26/RPS3A/RPS4X/
RPS4Y1/RPS5/RPS7
In the first column is reported the type of the enriched GO described in the second column. In the third column are reported the numbers of DEGs included in the GO versus the total number DEGs. In the fourth column are reported the numbers of DEGs included in the GO versus the number of DEGs in the subnetwork. In the last column are reported the symbols for the DEGs related to our subnetwork included in the enriched GO.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

D’Angiolini, S.; Lui, M.; Mazzon, E.; Calabrò, M. Network Analysis Performed on Transcriptomes of Parkinson’s Disease Patients Reveals Dysfunction in Protein Translation. Int. J. Mol. Sci. 2024, 25, 1299. https://doi.org/10.3390/ijms25021299

AMA Style

D’Angiolini S, Lui M, Mazzon E, Calabrò M. Network Analysis Performed on Transcriptomes of Parkinson’s Disease Patients Reveals Dysfunction in Protein Translation. International Journal of Molecular Sciences. 2024; 25(2):1299. https://doi.org/10.3390/ijms25021299

Chicago/Turabian Style

D’Angiolini, Simone, Maria Lui, Emanuela Mazzon, and Marco Calabrò. 2024. "Network Analysis Performed on Transcriptomes of Parkinson’s Disease Patients Reveals Dysfunction in Protein Translation" International Journal of Molecular Sciences 25, no. 2: 1299. https://doi.org/10.3390/ijms25021299

APA Style

D’Angiolini, S., Lui, M., Mazzon, E., & Calabrò, M. (2024). Network Analysis Performed on Transcriptomes of Parkinson’s Disease Patients Reveals Dysfunction in Protein Translation. International Journal of Molecular Sciences, 25(2), 1299. https://doi.org/10.3390/ijms25021299

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop