Screening of Parkinson’s Differential MicroRNA Based on GEO Database and Its Clinical Verification

Jiang, Xuping; Xiao, Lili; Jiang, Xumei; Li, Guangsheng; Lu, Zhijuan

doi:https://doi.org/10.1155/2021/8171236

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Data Availability Conflicts of Interest References Copyright Related Articles

Special Issue

Accelerating Precision Medicine through Integrated Analysis of Multiscale Biological Data

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 8171236 | https://doi.org/10.1155/2021/8171236

Screening of Parkinson’s Differential MicroRNA Based on GEO Database and Its Clinical Verification

Xuping Jiang,¹Lili Xiao,¹Xumei Jiang,²Guangsheng Li,¹and Zhijuan Lu¹

Academic Editor: Jun Yang

Received02 Jul 2021

Revised18 Oct 2021

Accepted19 Oct 2021

Published13 Nov 2021

Abstract

Objective. This study is set out to explore the potential difference of miR in PD through GEO data and provide diagnostic indicators for clinical practice. Methods. In this study, differential miR was screened through the Gene Expression Omnibus (GEO) database, 68 PD patients treated in our hospital from May 2017 to March 2018 were collected as the research group (RG), and 50 normal subjects who underwent physical examination in our hospital during the same period were collected as the control group (CG). Quantitative real-time polymerase chain reaction (qRT-PCR) was used to detect the expression and diagnostic value of miR-374a-5p in serum of patients. The potential target genes of miR-374a-5p were predicted, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis and Gene Ontology Consortium (GO) were carried out. Results. GEO2R analysis revealed that 193 miRs are expressed differentially, of which 78 were highly expressed and 115 were poorly expressed. The miR-374a-5p expression in the serum of the RG was reduced markedly and had a diagnostic value. Targetscan and miRDB online websites were used to predict their target genes, with 415 common target genes. miR-374a-5p may participate in 27 functional pathways and 8 signal pathways. Conclusion. miR-335-5p has low expression in PD and is expected to be a potential diagnostic indicator.

1. Introduction

Parkinson’s disease (PD) is the second clinically common neurodegenerative disease after Alzheimer’s disease, characterized by gradual loss of dopaminergic neurons in the substantia nigra [1, 2]. Statistics reveal that [3, 4] patients over 60 years old have one PD patient per 800-1000 people, and the incidence of that line is higher than that of women. At present, the relevant mechanism of PD is not clear. Studies have found that [3] PD may be caused by deletion of SNCA and PINK1 genes. However, this disease is sporadic, and its pathogenesis involves environment, epigenetic, etc. In addition, it has a long latent period, the early stage of the disease has no great influence on patients, and the lack of medical knowledge leads to their admission to enter the middle and late stages [5, 6]. With the continuous improvement of the medical level in recent years, PD prevention has been popularized, but patients’ daily life will be seriously affected in the middle and late stages of the disease [7].

MicroRNA (miR) is a short-chain noncoding RNA with a length of 20-25 nucleotides [8]. Studies show that [9, 10] miR can regulate downstream target genes by regulating the 3 untranslated region of downstream target gene mRNA. Previous studies have found that miR is involved in the occurrence and development of various diseases such as tumors, cardiovascular diseases, and neurological diseases [11–13]. Some studies have revealed that there are differential expression levels of various miRs in PD. For example, some studies have found that the miR-221 expression in the serum of [14] PD patients reduces, which is expected to become a potential diagnostic indicator of PD. Other studies have found that [15] miR-155 can regulate the inflammatory response of the Parkinson’s disease model induced by α-synuclein. The GEO database is one of the largest public chip databases in the world [16], which contains gene chips for various diseases. We have found that there are many differences in miR through analysis of the chips in the database this time, but it has not been clinically verified.

Therefore, in this research, we explored the differential genes in expression profile dataset GSE16658 and carried out clinical verification to find potential observation indicators for clinical use.

2. Materials and Methods

2.1. Patient Data

Sixty-eight PD patients treated in our hospital from May 2017 to March 2018 were collected as the research group (RG), and another 50 normal patients who underwent physical examination in our hospital during the same period were collected as the control group (CG). There was no statistical difference between the two groups in age and gender. This study was approved by the Medical Ethics Committee of our hospital. Inclusion criteria include patients who met Parkinson’s diagnostic criteria and those who met Hoehn and Yahr (H&Y) staging criteria [17]. Patients and their families were informed, and they all signed an informed consent form. Exclusion criteria include patients who were secondary or iatrogenic PD, who passed the mini-mental state examination (MMSE) score, and who had renal insufficiency, multiple tumors, and stroke.

MMSE refers to the test that measures global cognitive functioning on domains that include memory, attention, language, praxis, and visuospatial ability. The summed scores of MMSE ranged from 0 to 30, with higher values denoting better cognitive function.

2.2. GEO Chip Analysis

We checked Parkinson’s disease microRNA-related chips by logging into the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/gds) and selected GSE16658 for analysis. The miRNA expression dataset GSE16658 was downloaded from the National Center for Biotechnology Information GEO database (http://www.ncbi.nlm.nih.gov/geo). In the GSE16658 dataset, expression profiles were obtained from peripheral blood mononuclear cells from patients with PD and normal controls and were quantified using the miRCURY LNA microRNA Array, v.10.0-hsa, mmu & rno (Exiqon A/S, Vedbæk, Denmark). The analysis was performed by GEO2R software built in the GEO database. miRs with the top 10 differences were selected for heat mapping. The chip data are shown in Table 1.

2.3. Collection of Test Samples

We collected 5 mL of peripheral venous blood from people in the two groups and then subpacked it to the blood collection tube (5 mL) of inert separation gel and coagulant. Then, it was centrifuged at 3000 rpm for 10 min at 24°C to collect part of serum for subsequent experiments, and the surplus part was placed in an EP tube without RNA enzyme for later use at -80°C.

2.4. Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) Detection

Total RNA was extracted from the collected serum by the TRIzol kit (Invitrogen Company, USA), and its purity, concentration, and integrity were detected by ultraviolet spectrophotometer and agarose gel electrophoresis. Subsequently, reverse transcription was carried out using the TaqMan™ reverse transcription kit (Invitrogen Company, USA), and the transcription steps were strictly operated according to the kit instructions. The obtained cDNA was subjected to subsequent research. PCR amplification was carried out using the PrimeScript RT Master Mix kit (Takarabo Company, Japan). The amplification system was as below: 10 μL SYBR qPCR Mix, 0.8 μL for upstream and downstream primers, 2 μL cDNA product, 0.4 μL 50x ROX reference dye, and RNase-free water supplemented to 20 μL. PCR reaction conditions were as below: 95°C predenaturation for 60 s, 95°C denaturation for 30 s, and 60°C annealing extension for 40 s, with a total of 40 cycles. In the experiment, three parallel repeating holes were designed, and all specimens were repeatedly tested 3 times. miR used U6 as the internal reference and used 2^-ΔΔct to analyze the data [18]. The PCR instrument was 7500PCR instrument from ABI Company.

2.5. Bioinformatics Analysis

The online target gene prediction website was used to predict the target gene, and then, the clusterProfiler package in R software was used to carry out Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis and Gene Ontology Consortium (GO) and draw pictures.

2.6. Statistical Analysis

In this study, the SPSS20.0 software package was used to carry out statistical analysis on the collected data. GraphPad Prism 7 was used to draw the data picture. The usage (%) of the counting data was confirmed by the chi-squared test and expressed by . The K-S test was employed to analyze the data distribution. The measurement data were expressed by the . The comparison of normal distribution data between the two groups was conducted by the independent-samples -test, which was expressed by . ROC was used to analyze the diagnostic value of miR-374a-5p in PD. A value lower than 0.05 was statistically different.

3. Results

3.1. Analysis of Chip Results

Through analysis by online analysis software GEO2R, we found 193 miRs with differences, including 78 with high expression and 115 with low expression. We selected the top 10 miRs for display and combined the references and adjusted differences. Finally, we selected miR-335-5p and miR-374a-5p for further research. More details are shown in Figure 1 and Table 2.

(a)

(b)

3.2. Expression of miR-335-5p and miR-374a-5p in PD Patients

We detected the miR-335 and miR-374a-5p expression in the serum of patients. The results revealed that the miR-335-5p expression in the serum of the RG was not different from that in the CG, while the miR-374a-5p expression in the serum of patients reduced remarkably. Furthermore, we found that miR-374a-5p had a certain clinical value in the diagnosis of PD by drawing ROC curve analysis, and the area under the curve was 0.820. More details are shown in Figure 2.

(a)

(b)

(c)

3.3. Relationship between miR-374a-5p and H&Y Stage

We divided patients into three stages according to the H&Y stage, namely, the first stage (H&Y stage: 1-2, ), the second stage (H&Y stage: 3-4, ), and the third stage (H&Y stage: 5, ). We further detected the miR-374a-5p expression in patients of different stages and explored that miR-374a-5p was differentially expressed in the three stages. The ROC curve analysis displayed that miR-374a-5p had high clinical value in distinguishing patients at different stages. More details are shown in Figure 3 and Table 3.

(a)

(b)

3.4. Bioinformatics Analysis

In order to further explore the relevant mechanism of miR-374a-5p, we predicted its target genes through the online websites of Targetscan and miRDB. The results demonstrated that 415 common target genes existed. To further determine its potential mechanism, we employed the R software clusterProfiler package for GO and KEGG enrichment analysis. The results exhibited that miR-374a-5p might participate in 27 functional pathways and 8 signal pathways. More details are shown in Figure 4 and Tables 4 and 5.

(a)

(b)

(c)

4. Discussion

PD is a serious nervous system disease. The gradual loss of dopaminergic neurons in substantia nigra is one of the primary causes of PD [19]. What is more, due to the long incubation period of PD, patients are not clear about its early symptoms, which leads to PD patients being in the middle and late stages, thus causing them to miss the best treatment period [20]. Therefore, we urgently need to explore the relevant mechanism of PD.

In recent years, more and more studies have found a close relationship between miR and PD. For instance, Ma and others [21] pointed out that serum miR-221 could be used as a biomarker of Parkinson’s disease. Other studies have discovered that [22] inhibition of miR-34b and miR-34c enhances the α-synuclein expression in Parkinson’s disease, thus improving the disease condition. All the above studies signify that miR is involved in the occurrence and development of PD. In order to better screen potential PD differences of miR, we analyzed based on PD chips in the GEO database and found significant differences between miR-335-5p and miR-374a-5p, and the differences between the two after correction were the same, so we chose them for research.

Previous studies on miR-335-5p, miR-374a-5p, and PD are very few. To verify the expression and value of the two miRs in PD, we tested them. As a result, we found that there was no difference in the miR-335-5p expression in the serum of PD patients, while the miR-374a-5p slashed. Moreover, we also plotted the ROC curve and verified that miR-374a-5p had a high clinical value in distinguishing PD from normal people. Previous studies have found that miR-374a-5p is expressed in breast cancer, esophageal cancer, colon cancer, and other tumors and participates in the occurrence and development of tumors through relevant mechanisms [23–25]. This study was the first time we found that miR-374a-5p was differentially expressed in PD. To further determine the clinical value of miR-374a-5p in PD, we also analyzed the diagnostic value of miR-374a-5p in different H&Y stages. The H&Y stage is an important clinical score used to distinguish PD, which is divided into 5 stages. In this study, we further divided the patients into 3 stages based on the H&Y stage. Previously, Zhang and others divided them into groups according to this scheme. After grouping patients, we further detected the miR-374a-5p expression in patients’ serum. At last, we confirmed that miR-374a-5p was differentially expressed in those at different stages through Figure 3(a), and its expression decreased with the increase of stages, which indicated that PD was tied to the miR-374a-5p expression in the serum of severe patients. In addition, through ROC curve analysis, we found that miR-335-5p had a high clinical value in distinguishing patients of different stages, and the area under the curve was more than 0.8. Through the above research, we preliminarily confirmed the clinical value of miR-335-5p in PD, but we are still unclear about its relevant mechanism.

To further determine the relevant mechanism of miR-335-5p, we predicted its downstream target genes. Through the joint prediction of two online websites, we found 415 potential target genes downstream of miR-335-5p. Further enrichment analysis discovered that miR-335-5p participated in many functional pathways and signal pathways. Among them, the TNF signaling pathway, MAPK signaling pathway, and TGF-beta signaling pathway are important signal pathways for PD occurrence [26–28]. For instance, previous studies have found that [29] TNF signal inhibition in the central nervous system has an impact on normal brain function and neurodegenerative diseases. Other studies have verified that [30] activation of the MAPK signal pathway aggravates the condition of PD patients. Besides, studies have found that [31] TGF-β plays a role in the development, maintenance, and neuroprotection of dopamine neurons. Nevertheless, the analysis of this life letter provides us with the direction for future research.

There are still some limitations in this study. First of all, we only collected patient serum samples, and the samples were single. Some studies have found that differential expression of miR has also been detected in cerebrospinal fluid and peripheral blood mononuclear cells (PBMC). Moreover, the number of samples is relatively small, which may affect the data. Last but not least, we have not conducted in-depth research on the mechanism of miR-335-5p. Although we have conducted the credibility analysis to provide ideas for our research, we have not conducted experiments to confirm it. Hence, we hope to add more experiments and samples in future research to improve our research results.

To summarize, miR-335-5p has low expression in PD and is expected to become a potential diagnostic indicator.

Data Availability

The data used during the present study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

S. Sveinbjornsdottir, “The clinical symptoms of Parkinson's disease,” Journal of Neurochemistry, vol. 139, Supplement 1, pp. 318–324, 2016.
View at: Publisher Site | Google Scholar
R. B. Schneider, J. Iourinets, and I. H. Richard, “Parkinson's disease psychosis: presentation, diagnosis and management,” Neurodegenerative Disease Management, vol. 7, no. 6, pp. 365–376, 2017.
View at: Publisher Site | Google Scholar
L. M. Bekris, I. F. Mata, and C. P. Zabetian, “The genetics of Parkinson disease,” Journal of Geriatric Psychiatry and Neurology, vol. 23, no. 4, article 20938043, pp. 228–242, 2010.
View at: Publisher Site | Google Scholar
M. K. Lin and M. J. Farrer, “Genetics and genomics of Parkinson's disease,” Genome Medicine, vol. 6, no. 6, p. 48, 2014.
View at: Publisher Site | Google Scholar
A. Ascherio and M. A. Schwarzschild, “The epidemiology of Parkinson's disease: risk factors and prevention,” Lancet Neurology, vol. 15, no. 12, pp. 1257–1272, 2016.
View at: Publisher Site | Google Scholar
M. Grayson, “Parkinson's disease,” Nature, vol. 538, no. 7626, article S1, 2016.
View at: Publisher Site | Google Scholar
S. Lotankar, K. S. Prabhavalkar, and L. K. Bhatt, “Biomarkers for Parkinson's disease: recent advancement,” Neuroscience Bulletin, vol. 33, no. 5, pp. 585–597, 2017.
View at: Publisher Site | Google Scholar
A. M. Mohr and J. L. Mott, “Overview of microRNA biology,” Seminars in Liver Disease, vol. 35, no. 1, pp. 3–11, 2015.
View at: Publisher Site | Google Scholar
V. Agarwal, G. W. Bell, J. W. Nam, and D. P. Bartel, “Predicting effective microRNA target sites in mammalian mRNAs,” eLife, vol. 4, 2015.
View at: Publisher Site | Google Scholar
J. Fernandes, A. S. Vieira, J. C. Kramer-Soares et al., “Hippocampal microRNA-mRNA regulatory network is affected by physical exercise,” Biochimica et Biophysica Acta - General Subjects, vol. 1862, no. 8, pp. 1711–1720, 2018.
View at: Publisher Site | Google Scholar
C. K. Chou, R. T. Liu, and H. Y. Kang, “MicroRNA-146b: a novel biomarker and therapeutic target for human papillary thyroid cancer,” International Journal of Molecular Sciences, vol. 18, no. 3, 2017.
View at: Publisher Site | Google Scholar
A. Wojciechowska, A. Braniewska, and K. Kozar-Kaminska, “MicroRNA in cardiovascular biology and disease,” Advances in Clinical and Experimental Medicine, vol. 26, no. 5, pp. 865–874, 2017.
View at: Publisher Site | Google Scholar
D. S. Greenberg and H. Soreq, “MicroRNA therapeutics in neurological disease,” Current Pharmaceutical Design, vol. 20, no. 38, pp. 6022–6027, 2014.
View at: Publisher Site | Google Scholar
S. E. Oh, H. J. Park, L. He, C. Skibiel, E. Junn, and M. M. Mouradian, “The Parkinson's disease gene product DJ-1 modulates miR-221 to promote neuronal survival against oxidative stress,” Redox Biology, vol. 19, pp. 62–73, 2018.
View at: Publisher Site | Google Scholar
A. D. Thome, A. S. Harms, L. A. Volpicelli-Daley, and D. G. Standaert, “MicroRNA-155 regulates alpha-synuclein-induced inflammatory responses in models of Parkinson disease,” The Journal of Neuroscience, vol. 36, no. 8, pp. 2383–2390, 2016.
View at: Publisher Site | Google Scholar
T. Barrett, S. E. Wilhite, P. Ledoux et al., “NCBI GEO: archive for functional genomics data sets--update,” Nucleic Acids Research, vol. 41, no. Database issue, pp. D991–D995, 2013.
View at: Publisher Site | Google Scholar
X. Bai, Y. Tang, M. Yu et al., “Downregulation of blood serum microRNA 29 family in patients with Parkinson's disease,” Scientific Reports, vol. 7, no. 1, p. 5411, 2017.
View at: Publisher Site | Google Scholar
K. J. Livak and T. D. Schmittgen, “Analysis of relative gene expression data using real-time quantitative PCR and the 2(-delta delta C(T)) method,” Methods, vol. 25, no. 4, pp. 402–408, 2001.
View at: Publisher Site | Google Scholar
P. M. Antony, N. J. Diederich, R. Kruger, and R. Balling, “The hallmarks of Parkinson's disease,” The FEBS Journal, vol. 280, no. 23, pp. 5981–5993, 2013.
View at: Publisher Site | Google Scholar
A. Schrag, L. Horsfall, K. Walters, A. Noyce, and I. Petersen, “Prediagnostic presentations of Parkinson's disease in primary care: a case-control study,” Lancet Neurology, vol. 14, no. 1, pp. 57–64, 2015.
View at: Publisher Site | Google Scholar
W. Ma, Y. Li, C. Wang, F. Xu, M. Wang, and Y. Liu, “Serum miR-221 serves as a biomarker for Parkinson's disease,” Cell Biochemistry and Function, vol. 34, no. 7, pp. 511–515, 2016.
View at: Publisher Site | Google Scholar
S. Kabaria, D. C. Choi, A. D. Chaudhuri, M. M. Mouradian, and E. Junn, “Inhibition of miR-34b and miR-34c enhances alpha-synuclein expression in Parkinson's disease,” FEBS Letters, vol. 589, no. 3, pp. 319–325, 2015.
View at: Publisher Site | Google Scholar
D. Son, Y. Kim, S. Lim et al., “miR-374a-5p promotes tumor progression by targeting ARRB1 in triple negative breast cancer,” Cancer Letters, vol. 454, pp. 224–233, 2019.
View at: Publisher Site | Google Scholar
W. Chen, Y. Zhang, H. Wang, T. Pan, Y. Zhang, and C. Li, “LINC00473/miR-374a-5p regulates esophageal squamous cell carcinoma via targeting SPIN1 to weaken the effect of radiotherapy,” Journal of Cellular Biochemistry, vol. 120, no. 9, pp. 14562–14572, 2019.
View at: Publisher Site | Google Scholar
M. L. Slattery, A. J. Pellatt, F. Y. Lee et al., “Infrequently expressed miRNAs influence survival after diagnosis with colorectal cancer,” Oncotarget, vol. 8, no. 48, pp. 83845–83859, 2017.
View at: Publisher Site | Google Scholar
Q. S. Zhang, Y. Heng, Y. H. Yuan, and N. H. Chen, “Pathological alpha-synuclein exacerbates the progression of Parkinson's disease through microglial activation,” Toxicology Letters, vol. 265, pp. 30–37, 2017.
View at: Publisher Site | Google Scholar
Y. Dong, L. L. Han, and Z. X. Xu, “Suppressed microRNA-96 inhibits iNOS expression and dopaminergic neuron apoptosis through inactivating the MAPK signaling pathway by targeting CACNG5 in mice with Parkinson's disease,” Molecular Medicine, vol. 24, no. 1, p. 61, 2018.
View at: Publisher Site | Google Scholar
J. H. Koo, Y. C. Jang, D. J. Hwang et al., “Treadmill exercise produces neuroprotective effects in a murine model of Parkinson's disease by regulating the TLR2/MyD88/NF-kappaB signaling pathway,” Neuroscience, vol. 356, pp. 102–113, 2017.
View at: Publisher Site | Google Scholar
P. C. Tiwari and R. Pal, “The potential role of neuroinflammation and transcription factors in Parkinson disease,” Dialogues in Clinical Neuroscience, vol. 19, no. 1, pp. 71–80, 2017.
View at: Google Scholar
S. K. Jha, N. K. Jha, R. Kar, R. K. Ambasta, and P. Kumar, “p38 MAPK and PI3K/AKT signalling cascades inParkinson's disease,” International Journal of Molecular and Cellular Medicine, vol. 4, no. 2, pp. 67–86, 2015.
View at: Google Scholar
R. Calvello, A. Cianciulli, G. Nicolardi et al., “Vitamin D treatment attenuates neuroinflammation and dopaminergic neurodegeneration in an animal model of Parkinson's disease, shifting M1 to M2 microglia responses,” Journal of NeuroImmune Pharmacology, vol. 12, no. 2, pp. 327–339, 2017.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Xuping Jiang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies