Impact Factor 4.504

The Frontiers in Neuroscience journal series is the 1st most cited in Neurosciences

Original Research ARTICLE

Front. Aging Neurosci., 18 June 2018 | https://doi.org/10.3389/fnagi.2018.00178

Integrated Analysis and Identification of Novel Biomarkers in Parkinson’s Disease

Jieshan Chi1,2, Qizhi Xie1,2, Jingjing Jia1, Xiaoma Liu1, Jingjing Sun1, Yuanfei Deng1,3 and Li Yi1,3*
  • 1Department of Neurology, Peking University Shenzhen Hospital, Shenzhen, China
  • 2Department of Clinical Medicine, Shantou University Medical College, Shantou, China
  • 3National Clinical Research Center for Geriatric Diseases Shenzhen Center, Peking University Shenzhen Hospital, Shenzhen, China

Parkinson’s disease (PD) is a quite common neurodegenerative disorder with a prevalence of approximately 1:800–1,000 in subjects over 60 years old. The aim of our study was to determine the candidate target genes in PD through meta-analysis of multiple gene expression arrays datasets and to further combine mRNA and miRNA expression analyses to identify more convincing biological targets and their regulatory factors. Six included datasets were obtained from the Gene Expression Omnibus database by systematical search, including five mRNA datasets (150 substantia nigra samples in total) and one miRNA dataset containing 32 peripheral blood samples. A chip meta-analysis of five microarray data was conducted by using the metaDE package and 94 differentially expressed (DE) mRNAs were comprehensively obtained. And 19 deregulated DE miRNAs were obtained through the analysis of one miRNAs dataset by Qlucore Omics Explorer software. An interaction network formed by DE mRNAs, DE miRNAs, and important pathways was discovered after we analyzed the functional enrichment, protein–protein interactions, and miRNA targetome prediction analysis. In conclusion, this study suggested that five significantly downregulated mRNAs (MAPK8, CDC42, NDUFS1, COX4I1, and SDHC) and three significantly downregulated miRNAs (miR-126-5p, miR-19-3p, and miR-29a-3p) were potentially useful diagnostic markers in clinic, and lipid metabolism (especially non-alcoholic fatty liver disease pathway) and mitochondrial dysregulation may be the keys to biochemically detectable molecular defects. However, the role of these new biomarkers and molecular mechanisms in PD requires further experiments in vivo and in vitro and further clinical evidence.

Introduction

Parkinson’s disease (PD) is a quite common neurodegenerative disorder with a prevalence of approximately 1:800–1,000 in subjects over 60 years of age (Bekris et al., 2010; Lin and Farrer, 2014). The pathological hallmarks of PD are the degeneration of dopaminergic neurons in the substantia nigra (SN) pars compacta and Lewy bodies formed by the accumulation of α-synuclein (SNCA) in the remaining neurons (Spillantini et al., 1997; Braak et al., 2003). PD is a progressive disease characterized by motor signs such as resting tremor, bradykinesia, muscle rigidity, and posture abnormalities, but though in some patients, non-motor symptoms may appear, including depression, anxiety, sleep disturbance, fatigue, cognitive disorders, and gastrointestinal and sexual dysfunction (Schneider and Obeso, 2015).

Parkinson’s disease was originally thought to be a sporadic disease caused by environmental factors such as exposure to 1-methyl-4-phenyl-1,2,3,6-tetrahyropyridine (MPTP) or paraquat (Langston et al., 1983; Ascherio and Schwarzschild, 2016). However, epidemiological investigation in recent decades found that approximately 10% of patients had a family history of PD, and approximately 90% of patients had sporadic PD (Klein and Westenberger, 2012; Ascherio and Schwarzschild, 2016). It is generally accepted that PD is a multifactorial disease caused by the interaction between environmental and genetic factors, and approximately 5–10% of PD cases have a mutation in one of the several disease-associated genes according to the report (Farrer, 2006; Kalia and Lang, 2015). Several genes have been shown to be closely related to PD, such as α-synuclein (SNCA), leucine-rich repeat kinase 2 (LRRK2), Parkin RBRE3 ubiquitin protein ligase (PRKN), Parkinsonism associated deglycase (PARK7), PTEN-induced putative kinase 1 (PINK1), and ATPase 13A2 (ATP13A2) (Bekris et al., 2010; Nuytemans et al., 2010; Klein and Westenberger, 2012; Heman-Ackah et al., 2013). However because these mutations do not cause disease directly, it is more likely that they make people more susceptible to the PD when in cooperation with other risk factors (Lesage and Brice, 2009). In sporadic PD patients, SNCA and mitochondrial dysfunction are the predominant components of Lewy bodies, and complex I is found to be defective in the cytoplasm in the SN (Henchcliffe and Beal, 2008). However, what we know about the potential mechanism of PD is just the tip of the iceberg. At the genomic level, it is rather difficult to analyze the disease, especially sporadic PD, from a single gene mutation, and it is more the imbalance in gene expression and phenotypic changes were caused by a variety of regulatory mechanisms which act as mediators between genotype and phenotype. microRNA (miRNA) is considered as one of these mediators. Some deregulatory miRNAs were detected in the humans that lead to mitochondrial dysfunction, altered mitochondrial dynamics, oxidative stress, excitotoxicity, and the accumulation of SNCA, consequently resulting in neurodegeneration (Martin, 2010). Since miRNAs and its target genes are not in a one-to-one but one-to-many relationship (Qiu et al., 2015), the mutual regulation network between miRNAs and target genes can provide more insight into the disease and possible new treatments. Although the exploration of PD has been explored for nearly a century, the pathogenesis of PD is not yet clear.

Massively parallel microarray analysis allows for the global assessment of more credible relationships between gene expression and clinical manifestations in unbiased, and reveals the etiology of such complex diseases by identifying abnormalities in genes or pathways (Schadt et al., 2005; Scherzer, 2009). Large-scale genetic data can be categorized and detected based on phenotypic characteristics and produce hypotheses about the mechanisms of disease, which may have an underestimated role in the decoding of complex diseases (Scherzer et al., 2008). The utility of genome-wide expression data is often subjected to typical inconsistent analysis, non-replication, and small sample effects in practice.

To effectively reducing the bias of small sample studies and nominating deregulated genes, this integrated analysis consolidated data information on multiple datasets from different platforms. The aim of our work was to determine the candidate target genes in PD through meta-analysis and bioinformatic analysis of multiple datasets in gene expression and to further combine mRNA and miRNA expression analyses to seek out more convincing biological targets and their regulatory factors. In this study, we found several new potential mRNAs and miRNAs as well as one pathway in combination with the mRNA and miRNA microarray analysis and have mapped the intrinsic roles of genes and pathways in PD.

Materials and Methods

Materials and Data Pre-processing

Candidate microarrays associated with PD were acquired by retrieving the human gene expression studies deposited in the Gene Expression Omnibus (GEO) database1. Datasets related to PD were researched with the term “PD” of Homo sapiens (target). By September 14, 2017, a total of 1418 datasets were retrieved, including different types of samples and various types of expression data. The inclusion criteria were as follows: (1) original experimental studies; (2) human SN sample; (3) mRNA expression profile; (4) can obtain the unprocessed raw data (CEL files). The exclusion criteria were as follows: (1) repeated reports from the same institute or hospital; (2) non-human SN sample; (3) a non-expression gene chip, and the unprocessed raw data (CEL files) of these datasets were acquired from GEO database. All the included studies obtain relative ethics approval. All datasets were pre-processed individually (including background adjustment, normalization, summarization) on the base-2 logarithm by robust multi-array average (RMA) and annotated by converting different probe IDs to gene IDs by using R language. We use the Bioconductor software to compute RMA expression measures. Loaded the appropriate software with library(affy) to read all the CEL files in the current working directory. After loading the data, we compute the RMA expression measure. For miRNA microarray, the unqualified chips would be retrieved, and the samples for RNA detection were human peripheral blood. The miRNA dataset was imported into the Qlucore Omics Explorer (QOE) software for data pre-processing (mean = 0, SD = 1).

Integrated Analysis of Gene Expression Datasets

Included gene expression datasets were loaded into R language for objective quality control by MetaQC package, which intended to identify whether the included chip was qualified for genomic meta-analysis. The MetaQC package provides four quantitative quality control indexes, including internal quality control (IQC), external quality control (EQC), accuracy quality control of differentially expressed (DE) gene detection (AQCg), or pathway identification (AQCp) and consistency quality control in genes (CQCg) or pathways (CQCp). A principal component analysis (PCA) was performed to further visualize the quality control results. Eligible microarrays were subjected to threshold screening individually to obtain DE genes in PD under specific conditions using the Linear Models for Microarray (LIMMA) package. For integrated analysis, we further carried out the genomic meta-analysis. Considering the feasibility of the methods, a modified two-sample t-test by adding a fudging parameter was used to extrapolate the P-values and Fisher’s method was implemented for statistical analysis of significance (Tseng et al., 2012). A corrected P-value (P < 0.05) was considered statistically significant for the DE mRNAs. For visualization, DE mRNAs based on specific fault discovery rate (FDR) values (FDR < 0.0001) were plotted by the MetaDE package (heatmap.sig.genes).

Functional Analysis of PD-Related DE Genes

To access the prospective functions of PD-related DE genes found in the meta-analysis, online tools such as the Database for Annotation, Visualization and Integrated Discovery (DAVID) (Huang da et al., 2009a,b) were used. The functional categories Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and Gene Ontology (GO) terms were analyzed. For narrated KEGG pathways and GO terms that enriched the target genes, P-values less than 0.05 were defined as the cut-off criterion. The critical assessment and integration of protein–protein interactions (PPIs) of PD-related DE genes, including direct as well as indirect associations, were analyzed by the STRING database (Szklarczyk et al., 2015). The PPI networks were visualized by Cytoscape software (Shannon et al., 2003).

Analysis of miRNA Expression Dataset and Target Prediction

Qlucore Omics Explorer software was primarily capable of screening for DE miRNAs (Jimbo et al., 2011). Primarily, QOE software performed the normalized processing of mean = 0 and var = 1 for the original data, and further calculated the corresponding standard error σ for each variable. We filtered out the variables with small differences among samples according to the statistics by σ/σmax. After data preprocessing, the comparison between two groups was performed and specific filter variable parameters were selected for statistical analysis of significance (P < 0.001, q < 0.01, fold change > 2). The predicted target genes of DE miRNAs were identified by three different target prediction algorithms: miRDB2 (Wong and Wang, 2015) TargetScan 7.13 (Agarwal et al., 2015) and microT_CDS of Diana Tools4 (Reczko et al., 2012; Paraskevopoulou et al., 2013). Unique genes with target sites on 3′-UTR were incorporated. To reduce the false positive rate and improve persuasion, we obtained the overlap target genes of the three algorithms mentioned above.

Results

Integrated Analysis of PD Gene Expression Datasets

Six primary datasets with available mRNA expression data for SN samples in PD patients and controls were identified by searching the GEO database (GSE20186, GSE8397, GSE20141, GSE20333, GSE7621, and GSE20295). One of these eligible datasets (GSE20295) was excluded after the quality control using the MetaQC package and visualization by PCA (Figure 1A, additional data were given in Online Resource Supplementary Table S1). The remaining five original datasets (accession numbers GSE20186, GSE8397, GSE20141, GSE20333, and GSE7621), six microarrays in total (dataset GSE8397 contained two mRNA chips: 39 individual tissue samples were tested using one A and one B mRNA chip per sample), were incorporated in our meta-analysis, for a total of 150 independent SN samples (80 PD patients and 70 controls, Table 1). Each dataset was analyzed independently to obtain its own DE mRNAs before performing pooled analysis. Approximately 193 DE mRNAs (hereinafter referred to as overlap DE mRNAs) were confirmed to co-exist in all eligible datasets. After merging the datasets, we identified 94 mRNAs (hereinafter referred to as meta-DE mRNAs) showing consistent DE patterns using a penalized t-test by adding a fudging parameter, and the maximum P-value and Fisher’s method by summarizing -log (P-value) across studies were chosen to eliminate the significant influence of the large number of samples (Lu et al., 2010). The heatmap of the meta-DE mRNAs was plotted and visualized by setting a particular FDR (<0.00001, Figure 1B and additional data were given in Online Resource Supplementary Figure S1), which demonstrated that meta-DE mRNAs were all downregulated (meta-DE mRNAs expression in each dataset are shown in Online Resource Supplementary Table S2, and the details of each meta-DE mRNA are given in Online Resource Supplementary Table S3).

FIGURE 1
www.frontiersin.org

FIGURE 1. (A) Quality control result. (B) DE genes through the method of meta-DE. (For the sake of space, only a part of the figure is shown here. The full size figure is displayed in Online Resource Supplementary Figure S1.

TABLE 1
www.frontiersin.org

TABLE 1. Baseline characteristics of datasets.

Enrichment Analysis of the DE Gene

Through the enrichment analysis of KEGG pathways, meta-DE genes were mainly involved in the Huntington’s disease pathway, PD pathway, Alzheimer’s disease pathway, non-alcoholic fatty liver disease (NAFLD) pathway, and the oxidative phosphorylation pathway (Figure 2D). The results of enrichment analysis using three categories of GO were as follows: (1) biological processes: mitochondrial ATP synthesis-coupled electron transport, mitochondrial membrane organization, modulation by virus of host process, branched-chain amino acid catabolic process, aerobic respiration, etc. (Figure 2A); (2) molecular functions: oxidoreductase activity acting on NAD(P)H, oxidoreductase activity of a heme group of donors, SMAD binding, and oxidoreductase activity of a sulfur group of donors (Figure 2B); and (3) cellular component: mitochondrial respiratory chain, respiratory chain complex, organelle envelope lumen, and rough endoplasmic reticulum (Figure 2C). The enrichment analysis by GO revealed that DE genes were mainly involved in biological processes including the ATP metabolic process and electron transfer, and the main organelles in these processes were mitochondria, which were consistent with the results of GO enrichment with cellular components.

FIGURE 2
www.frontiersin.org

FIGURE 2. Functional enrichment analysis of meta-DE genes. (A) Cellular components of GO enrichment analysis. (B) Biological processes of GO enrichment analysis. (C) Molecular functions of GO enrichment analysis. (D) KEGG pathway enrichment analysis.

In addition, we established PPI networks using the STRING database to further investigate the functional partnership and interaction of overlapping DE genes (Figure 3A) and meta-DE genes (Figure 3B), respectively, and the results were visualized by Cytoscape software. Interestingly, MAPK8, CDC42, NDUFS1, COX4I1, RAD23B, and SDHC, respectively, were in the core position of two PPI networks. The significant level of these six core genes was showed in Table 2. Once these hub nodes were expurgated, both PPI networks became inattentive. Of these core genes, NDUFS1, COX4I1, and SDHC, were significantly linked to mitochondrial function which means these three DE genes were potentially related to PD. Furthermore, in addition to the known pathways, we found a pathway (NAFLD pathway) that was not previously reported to be closely associated with PD. Of note, two of these core genes (MAPK8, CDC42) were significantly involved in the NAFLD pathway.

FIGURE 3
www.frontiersin.org

FIGURE 3. PPI networks. (A) PPI network of overlap DE genes. (B) PPI network of meta-DE genes. The size and color of map nodes are determined by the degree value, which renders a gradual process by setting the small size with a low degree in green, large size with a high degree in orange.

TABLE 2
www.frontiersin.org

TABLE 2. The significant level of six core genes.

Analysis of the PD miRNA Expression Dataset

We identified the GSE16658 dataset with available miRNA expression data containing 32 peripheral blood samples of PD patients (19 samples) and controls (13 samples) by searching the GEO database. Nineteen DE miRNAs were identified using QOE software; among them, the top five most significant DE miRNAs were miR-199-3p (miR-199a-3p/miR-199b-3p), miR-126-5p, miR-29a-3p, miR-19b-3p, and miR-301a-3p (Table 3). The miRNA target genes were obtained from experimentally supported databases with experimentally verified and different prediction algorithms. The consensus of the target genes of the top five miRNAs was summarized (Figure 4, additional data are provided in Online Resource Supplementary Table S4). Here, we identified 11 genes coexisted in the union set of all target genes and meta-DE mRNAs (NDUFS1, MAPK8, CDC42, SNCA, VAPA, USP13, TIMM8B, KIF3A, KPNA6, MTCH2, and SUB1).

TABLE 3
www.frontiersin.org

TABLE 3. The DE miRNAs.

FIGURE 4
www.frontiersin.org

FIGURE 4. Summary results of the consensus of targets genes of top five miRNAs.

Comprehensive Analysis of DE Genes and miRNAs

To further investigate the relationship between top 5 DE miRNAs and 11 DE mRNAs mentioned above, we plotted their regulatory networks by Cytoscape. The results showed that the interaction between miRNAs and genes was formed with miR-126-5p, miR-29a-3p, and miR-19b-3p as the center (Figure 5A). Each miRNA interacted with several meta-DE genes, including miR-126-5p with MTCH2, VAVP, SNCA, and NDUFS1; miR-29a-3p with CDC42, TIMM8B, and SUB1; and miR-19b-3p with MAPK8, USP13, and KPNA6. Here, we outline the relationship analyzed above, including genes and pathways, particularly the main underlying pathway related to PD clustering by meta-DE genes (Figure 5B). Altogether, we extracted the most potential biomarkers, five mRNAs and three miRNAs, as well as one new pathway from the whole analysis as a matter of priority (Table 4).

FIGURE 5
www.frontiersin.org

FIGURE 5. The interaction between miRNAs and genes. (A) The relationship of gene–gene and gene–miRNA. (B) The outline of the interaction between the significant KEGG pathways of meta-DE genes and miRNAs.

TABLE 4
www.frontiersin.org

TABLE 4. Priority genes, related miRNAs, and functions.

Discussion

We utilized a careful method to standardize and unify the cross-platform datasets of gene profiling in PD for further integrated analyses. The integrated analysis method described here aimed to reduce the bias of small sample studies and nominating deregulated genes effectively. We have included six gene expression data point from five different platforms in PD and then performed meta-analysis using the MetaDE package in R language for merging and gene filtering of gene expressions (Wang et al., 2012). We identified 94 down regulated DE mRNAs in SN samples of PD and 19 deregulated miRNAs in PB samples of PD. Of these DE mRNAs and miRNAs, we deemed five genes (NDUFS1, COX4I1, SDHC, CDC42, and MAPK8) and three miRNAs (miR-126-5p, miR-29a-3p, and miR-19b-3p) worth exploring due to their relationship with the NAFLD pathway or the biological process of mitochondrial function.

Previous studies have suggested that several miRNAs, such as miR-133b (Kim et al., 2007), miR-16 (Zhang and Cheng, 2014), miR-153 (Doxakis, 2010), miR-205 (Cho et al., 2013), miR-7 (Junn et al., 2009; Doxakis, 2010), miR-64, miR-65, and the let-7 family (Asikainen et al., 2010), were related to PD. In our research, 19 significantly DE miRNAs were identified including let-7f and let-7g (Asikainen et al., 2010), which were confirmed previously and several new miRNAs. To further verify the results of our integrated analysis in a more comprehensive way, we compared the target genes predicted by the top five DE miRNAs with the analyzed DE mRNAs, and we initially conjectured that miR-126-5p, miR-19-3p, and miR-29a-3p were more likely to participate in the pathogenesis of PD.

The enrichment analysis by GO revealed that DE genes were mainly involved in biological processes including ATP metabolism process and electron transfer, and the main organelles in these processes were mitochondria, which was consistent with the results of GO enrichment with cellular components. Mitochondrial dysfunction is closely linked to the occurrence and development of PD (Stoessl, 1999; Koyano et al., 2013). Damage to the electron transport chain, the inhibition of complex I activity, mishandling of calcium, enhanced sensitivity to mitochondrial toxins, oxidative stress, and dysfunctional mitochondrial dynamics could contribute to neuronal dysfunction and participate in the pathogenesis of PD (Schapira et al., 1990; McCoy and Cookson, 2012). As an important energy-producing organelle in cells, mitochondria oxidize carbon and further produced ATP through oxidative phosphorylation (Wallace et al., 2010). Accordingly, damage to mitochondria in any process may disrupt the energy balance and promote the occurrence of the disease. Combining the results of two PPI networks, the core genes linked to mitochondrial function included three mitochondrial-localized genes (NDUFS1, COX4I1, and SDHC).

Both NDUFS1 and COX4I1 are located at the mitochondrial inner membrane and play an integral role in the electron transport chain. It was noted that pyrroloquinoline quinone can protect SH-SY5Y cells from the cytotoxicity induced by a mitochondrial complex I inhibitor (Zhang et al., 2014). However, silencing NDUFS1 in midbrain neurons or SH-SY5Y cells reduced the neuroprotective effect of pyrroloquinoline quinone (Zhang et al., 2016), which means that the down regulation of NDUFS1 may be detrimental to mitochondrial function. As a core subunit of ubiquinone oxidoreductase, the deregulation of NDUFS1 might result from the pathogenesis of PD by undermining mitochondrial function. In addition, production coded by COX4I1 is the terminal enzyme of the mitochondrial respiratory chain. A large-scale study of PD gene expression profiling demonstrated that COX4I1 was downregulated in whole-blood in patients with PD (Shamir et al., 2017). The consensus of deregulation both in blood and SN indicated that COX4I1 may be a potential biomarker in PD. There is currently no significant evidence that SDHC is involved in the pathogenesis of PD. However, its coding product as a member of nuclear-encoded subunits that comprises succinate dehydrogenase (also known as mitochondrial complex II) is a key enzyme complex of the tricarboxylic acid cycle and aerobic respiratory chains of mitochondria, which shows that the relationship between SDHC and PD deserves further study.

Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis for meta-DE genes shows that the genes relate to the neurodegenerative disease signaling pathway, including PD, Alzheimer’s disease and Huntington’s disease, and the NAFLD pathway. Neurodegenerative diseases share fundamental processes, such as mitochondrial anomalies and oxidative stress, in spite of their distinct pathological and clinical features (Di Carlo et al., 2012). Simultaneously, cholinergic deficit is a common pathogenesis among these three diseases and leads to analogous clinical manifestations such as dyskinesia (Pepeu and Grazia Giovannini, 2017). These prevailing conclusions support the effectiveness and validity of our integrated analysis.

Here, we proposed that the NAFLD pathway may probably be a new PD-related pathway that has been reportedly linked to neurodegeneration, especially AD. The study of APP-Tg mice indicates that NAFLD offers the opportunity to accelerate the symptoms of AD (Kim et al., 2016). A long-term high-fat diet induces systemic inflammation that stimulates the central nervous system to cause neurodegeneration (Kim et al., 2016). Studies that explore the association of NAFLD with PD are currently lacking, however, evidence from clinical and experimental investigations is mounting that alterations in lipid metabolism participate in the pathogenesis through direct crosstalk between lipids and SNCA (Ruiperez et al., 2010). Noticeably, two meta-DE genes (MAPK8 and CDC42) (Sharma et al., 2012; Wang et al., 2016), including their related miRNAs (miR-19b-3p and miR-29a-3p, respectively), participate in the NAFLD pathway. An in vivo and in vitro study confirms that miR-29 can significantly inhibit HMGCR expression by targeting 3′-UTR of HMGCR mRNA and participate in the accumulation of free cholesterol in the livers of mice with non-alcoholic steatohepatitis (Liu et al., 2017). Therefore, miR-29 may be an important regulator of hepatic cholesterol homeostasis and a potential therapeutic target for the treatment of NAFLD and other liver diseases associated with free cholesterol accumulation. According to our integrated analysis, it is reasonable to assume that miR-29a-3p is involved in the development and progression of PD through the expression of MAPK8 in the NAFLD pathway. Paradoxically, the inhibition of the MAPK8 (also named JNK1) was reported to be a neuroprotective factor (Busquets et al., 2017) which seems contradictory to our analytic results. The reasons for this phenomenon are obscure and need further verification and exploration. Some reports have illustrated the role of CDC42 in the pathophysiology of inherited PD (Musilli et al., 2016), and similarly, CDC42 is deregulated in NAFLD (Wang et al., 2016) which may significantly indicate the pivotal role of CDC42 in PD especially though the potential NAFLD pathway. Its related miRNA (miR-19b-3p) may participate in neurodegeneration by modulating neural cell apoptosis (Wang et al., 2016). We hereby propose that it is meaningful to conduct in-depth follow-up studies on the relationship between the NAFLD pathway and PD to elucidate the pathogenesis of PD and find new therapeutic targets or preventive measures.

This study has several limitations. One is the limited number of included datasets, especially the miRNA datasets. The other is the limited sample size. In the preliminary protocol, we intended to have brain, blood, and cerebrospinal fluid samples included in the mRNA and miRNA microarray. However, some of these sample types were not included in the current study due to the insufficiency of datasets and the inaccessibility of the raw data. Therefore, we will monitor this progress in PD. Further investigations are required to confirm whether our new biomarkers are potential prognostic predictors or therapeutic targets in PD.

In summary, our integrated analysis of PD genomics provides us with a wealth of resources to explore the role of target genes and miRNAs in PD. Five significantly downregulated genes (MAPK8, CDC42, NDUFS1, COX4I1, and SDHC) and three significantly downregulated miRNAs (miR-126-5p, miR-19-3p, and miR-29a-3p) were potentially useful clinical diagnostic markers. Lipid metabolism (especially NAFLD pathway) and mitochondrial dysregulation may underlie biochemically detectable deeper molecular defects. Future work should be focused on these two aspects to reveal the pathogenesis of PD and to develop new therapeutic targets for the clinical treatment of PD.

Author Contributions

JC and LY designed the protocol. JC, QX, and JJ collected the data. JS and XL standardized the raw data. JC and YD analyzed the experimental results. JC and LY wrote the manuscript.

Funding

This project was supported by grants from the Shenzhen Science and Technology Project (No. 201401027).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2018.00178/full#supplementary-material

Footnotes

  1. ^https://www.ncbi.nlm.nih.gov/geo/
  2. ^http://www.mirdb.org
  3. ^http://www.targetscan.org/vert_71/
  4. ^http://diana.imis.athena-innovation.gr/DianaTools

References

Agarwal, V., Bell, G. W., Nam, J. W., and Bartel, D. P. (2015). Predicting effective microRNA target sites in mammalian mRNAs. eLife 4:e05005. doi: 10.7554/eLife.05005

PubMed Abstract | CrossRef Full Text | Google Scholar

Ascherio, A., and Schwarzschild, M. A. (2016). The epidemiology of Parkinson’s disease: risk factors and prevention. Lancet Neurol. 15, 1257–1272. doi: 10.1016/s1474-4422(16)30230-7

CrossRef Full Text | Google Scholar

Asikainen, S., Rudgalvyte, M., Heikkinen, L., Louhiranta, K., Lakso, M., Wong, G., et al. (2010). Global microRNA expression profiling of Caenorhabditis elegans Parkinson’s disease models. J. Mol. Neurosci. 41, 210–218. doi: 10.1007/s12031-009-9325-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Bekris, L. M., Mata, I. F., and Zabetian, C. P. (2010). The genetics of Parkinson disease. J. Geriatr. Psychiatry Neurol. 23, 228–242. doi: 10.1177/0891988710383572

PubMed Abstract | CrossRef Full Text | Google Scholar

Braak, H., Del Tredici, K., Rub, U., de Vos, R. A., Jansen Steur, E. N., and Braak, E. (2003). Staging of brain pathology related to sporadic Parkinson’s disease. Neurobiol. Aging 24, 197–211.

Google Scholar

Busquets, O., Ettcheto, M., Verdaguer, E., Dario, R., Auladell, C., and Beas-Zarate, C. (2017). JNK1 inhibition by Licochalcone A leads to neuronal protection against excitotoxic insults derived of kainic acid. Neuropharmacology 62, 628–637. doi: 10.1016/j.neuropharm.2017.10.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Cho, H. J., Liu, G., Jin, S. M., Parisiadou, L., Xie, C., Yu, J., et al. (2013). MicroRNA-205 regulates the expression of Parkinson’s disease-related leucine-rich repeat kinase 2 protein. Hum. Mol. Genet. 22, 608–620. doi: 10.1093/hmg/dds470

PubMed Abstract | CrossRef Full Text | Google Scholar

Di Carlo, M., Giacomazza, D., Picone, P., Nuzzo, D., and San Biagio, P. L. (2012). Are oxidative stress and mitochondrial dysfunction the key players in the neurodegenerative diseases? Free Radic Res. 46, 1327–1338. doi: 10.3109/10715762.2012.714466

PubMed Abstract | CrossRef Full Text | Google Scholar

Doxakis, E. (2010). Post-transcriptional regulation of alpha-synuclein expression by mir-7 and mir-153. J. Biol. Chem. 285, 12726–12734. doi: 10.1074/jbc.M109.086827

PubMed Abstract | CrossRef Full Text | Google Scholar

Farrer, M. J. (2006). Genetics of Parkinson disease: paradigm shifts and future prospects. Nat. Rev. Genet. 7, 306–318. doi: 10.1038/nrg1831

PubMed Abstract | CrossRef Full Text | Google Scholar

Heman-Ackah, S. M., Hallegger, M., Rao, M. S., and Wood, M. J. (2013). RISC in PD: the impact of microRNAs in Parkinson’s disease cellular and molecular pathogenesis. Front. Mol. Neurosci. 6:40. doi: 10.3389/fnmol.2013.00040

PubMed Abstract | CrossRef Full Text | Google Scholar

Henchcliffe, C., and Beal, M. F. (2008). Mitochondrial biology and oxidative stress in Parkinson disease pathogenesis. Nat. Clin. Pract. Neurol. 4, 600–609. doi: 10.1038/ncpneuro0924

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang da, W., Sherman, B. T., and Lempicki, R. A. (2009a). Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 37, 1–13. doi: 10.1093/nar/gkn923

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang da, W., Sherman, B. T., and Lempicki, R. A. (2009b). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57. doi: 10.1038/nprot.2008.211

PubMed Abstract | CrossRef Full Text | Google Scholar

Jimbo, R., Fernandez-Rodriguez, J., Sul, Y. T., and Johansson, C. B. (2011). Principal component analysis: a novel analysis to evaluate the characteristics of osseointegration of different implant surfaces. Implant Dent. 20, 364–368. doi: 10.1097/ID.0b013e31822b9c39

PubMed Abstract | CrossRef Full Text | Google Scholar

Junn, E., Lee, K. W., Jeong, B. S., Chan, T. W., Im, J. Y., and Mouradian, M. M. (2009). Repression of alpha-synuclein expression and toxicity by microRNA-7. Proc. Natl. Acad. Sci. U.S.A. 106, 13052–13057. doi: 10.1073/pnas.0906277106

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalia, L. V., and Lang, A. E. (2015). Parkinson’s disease. Lancet 386, 896–912. doi: 10.1016/s0140-6736(14)61393-3

CrossRef Full Text

Kim, D. G., Krenz, A., Toussaint, L. E., Maurer, K. J., Robinson, S. A., Yan, A., et al. (2016). Non-alcoholic fatty liver disease induces signs of Alzheimer’s disease (AD) in wild-type mice and accelerates pathological signs of AD in an AD model. J. Neuroinflammation 13:1. doi: 10.1186/s12974-015-0467-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, J., Inoue, K., Ishii, J., Vanti, W. B., Voronov, S. V., Murchison, E., et al. (2007). A MicroRNA feedback circuit in midbrain dopamine neurons. Science 317, 1220–1224. doi: 10.1126/science.1140481

PubMed Abstract | CrossRef Full Text | Google Scholar

Klein, C., and Westenberger, A. (2012). Genetics of Parkinson’s disease. Cold Spring Harb. Perspect. Med. 2:a008888. doi: 10.1101/cshperspect.a008888

PubMed Abstract | CrossRef Full Text | Google Scholar

Koyano, F., Okatsu, K., Ishigaki, S., Fujioka, Y., Kimura, M., Sobue, G., et al. (2013). The principal PINK1 and Parkin cellular events triggered in response to dissipation of mitochondrial membrane potential occur in primary neurons. Genes Cells 18, 672–681. doi: 10.1111/gtc.12066

PubMed Abstract | CrossRef Full Text | Google Scholar

Langston, J. W., Ballard, P., Tetrud, J. W., and Irwin, I. (1983). Chronic Parkinsonism in humans due to a product of meperidine-analog synthesis. Science 219, 979–980.

PubMed Abstract | Google Scholar

Lesage, S., and Brice, A. (2009). Parkinson’s disease: from monogenic forms to genetic susceptibility factors. Hum. Mol. Genet. 18, R48–R59. doi: 10.1093/hmg/ddp012

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, M. K., and Farrer, M. J. (2014). Genetics and genomics of Parkinson’s disease. Genome Med. 6:48. doi: 10.1186/gm566

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, M. X., Gao, M., Li, C. Z., Yu, C. Z., Yan, H., Peng, C., et al. (2017). Dicer1/miR-29/HMGCR axis contributes to hepatic free cholesterol accumulation in mouse non-alcoholic steatohepatitis. Acta Pharmacol. Sin. 38, 660–671. doi: 10.1038/aps.2016.158

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, S., Li, J., Song, C., Shen, K., and Tseng, G. C. (2010). Biomarker detection in the integration of multiple multi-class genomic studies. Bioinformatics 26, 333–340. doi: 10.1093/bioinformatics/btp669

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, L. J. (2010). Mitochondrial and cell death mechanisms in neurodegenerative diseases. Pharmaceuticals 3, 839–915. doi: 10.3390/ph3040839

PubMed Abstract | CrossRef Full Text | Google Scholar

McCoy, M. K., and Cookson, M. R. (2012). Mitochondrial quality control and dynamics in Parkinson’s disease. Antioxid. Redox Signal. 16, 869–882. doi: 10.1089/ars.2011.4019

PubMed Abstract | CrossRef Full Text | Google Scholar

Musilli, M., Ciotti, M. T., Pieri, M., Martino, A., Borrelli, S., Dinallo, V., et al. (2016). Therapeutic effects of the Rho GTPase modulator CNF1 in a model of Parkinson’s disease. Neuropharmacology 109, 357–365. doi: 10.1016/j.neuropharm.2016.06.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Nuytemans, K., Theuns, J., Cruts, M., and Van Broeckhoven, C. (2010). Genetic etiology of Parkinson disease associated with mutations in the SNCA, PARK2, PINK1, PARK7, and LRRK2 genes: a mutation update. Hum. Mutat. 31, 763–780. doi: 10.1002/humu.21277

PubMed Abstract | CrossRef Full Text | Google Scholar

Paraskevopoulou, M. D., Georgakilas, G., Kostoulas, N., Vlachos, I. S., Vergoulis, T., Reczko, M., et al. (2013). DIANA-microT web server v5.0: service integration into miRNA functional analysis workflows. Nucleic Acids Res. 41, W169–W173. doi: 10.1093/nar/gkt393

PubMed Abstract | CrossRef Full Text | Google Scholar

Pepeu, G., and Grazia Giovannini, M. (2017). The fate of the brain cholinergic neurons in neurodegenerative diseases. Brain Res. 1670, 173–184. doi: 10.1016/j.brainres.2017.06.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, L., Tan, E. K., and Zeng, L. (2015). microRNAs and neurodegenerative diseases. Adv. Exp. Med. Biol. 888, 85–105. doi: 10.1007/978-3-319-22671-2_6

PubMed Abstract | CrossRef Full Text | Google Scholar

Reczko, M., Maragkakis, M., Alexiou, P., Grosse, I., and Hatzigeorgiou, A. G. (2012). Functional microRNA targets in protein coding sequences. Bioinformatics 28, 771–776. doi: 10.1093/bioinformatics/bts043

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruiperez, V., Darios, F., and Davletov, B. (2010). Alpha-synuclein, lipids and Parkinson’s disease. Prog. Lipid Res. 49, 420–428. doi: 10.1016/j.plipres.2010.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Schadt, E. E., Lamb, J., Yang, X., Zhu, J., Edwards, S., Guhathakurta, D., et al. (2005). An integrative genomics approach to infer causal associations between gene expression and disease. Nat. Genet. 37, 710–717. doi: 10.1038/ng1589

PubMed Abstract | CrossRef Full Text | Google Scholar

Schapira, A. H., Cooper, J. M., Dexter, D., Clark, J. B., Jenner, P., and Marsden, C. D. (1990). Mitochondrial complex I deficiency in Parkinson’s disease. J. Neurochem. 54, 823–827.

Google Scholar

Scherzer, C. R. (2009). Chipping away at diagnostics for neurodegenerative diseases. Neurobiol. Dis. 35, 148–156. doi: 10.1016/j.nbd.2009.02.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Scherzer, C. R., Grass, J. A., Liao, Z., Pepivani, I., Zheng, B., Eklund, A. C., et al. (2008). GATA transcription factors directly regulate the Parkinson’s disease-linked gene alpha-synuclein. Proc. Natl. Acad. Sci. U.S.A. 105, 10907–10912. doi: 10.1073/pnas.0802437105

PubMed Abstract | CrossRef Full Text | Google Scholar

Schneider, S. A., and Obeso, J. A. (2015). Clinical and pathological features of Parkinson’s disease. Curr. Top. Behav. Neurosci. 22, 205–220. doi: 10.1007/7854_2014_317

PubMed Abstract | CrossRef Full Text | Google Scholar

Shamir, R., Klein, C., Amar, D., Vollstedt, E. J., Bonin, M., Usenovic, M., et al. (2017). Analysis of blood-based gene expression in idiopathic Parkinson disease. Neurology 89, 1676–1683. doi: 10.1212/wnl.0000000000004516

PubMed Abstract | CrossRef Full Text | Google Scholar

Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, M., Urano, F., and Jaeschke, A. (2012). Cdc42 and Rac1 are major contributors to the saturated fatty acid-stimulated JNK pathway in hepatocytes. J. Hepatol. 56, 192–198. doi: 10.1016/j.jhep.2011.03.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Spillantini, M. G., Schmidt, M. L., Lee, V. M., Trojanowski, J. Q., Jakes, R., and Goedert, M. (1997). Alpha-synuclein in Lewy bodies. Nature 388, 839–840. doi: 10.1038/42166

PubMed Abstract | CrossRef Full Text | Google Scholar

Stoessl, A. J. (1999). Etiology of Parkinson’s disease. Can. J. Neurol. Sci. 26(Suppl. 2), S5–S12.

Google Scholar

Szklarczyk, D., Franceschini, A., Wyder, S., Forslund, K., Heller, D., Huerta-Cepas, J., et al. (2015). STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 43, D447–D452. doi: 10.1093/nar/gku1003

PubMed Abstract | CrossRef Full Text | Google Scholar

Tseng, G. C., Ghosh, D., and Feingold, E. (2012). Comprehensive literature review and statistical considerations for microarray meta-analysis. Nucleic Acids Res. 40, 3785–3799. doi: 10.1093/nar/gkr1265

PubMed Abstract | CrossRef Full Text | Google Scholar

Wallace, D. C., Fan, W., and Procaccio, V. (2010). Mitochondrial energetics and therapeutics. Annu. Rev. Pathol. 5, 297–348. doi: 10.1146/annurev.pathol.4.110807.092314

CrossRef Full Text | Google Scholar

Wang, R., Wang, X., and Zhuang, L. (2016). Gene expression profiling reveals key genes and pathways related to the development of non-alcoholic fatty liver disease. Ann. Hepatol. 15, 190–199. doi: 10.5604/16652681.1193709

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Kang, D. D., Shen, K., Song, C., Lu, S., Chang, L. C., et al. (2012). An R package suite for microarray meta-analysis in quality control, differentially expressed gene analysis and pathway enrichment detection. Bioinformatics 28, 2534–2536. doi: 10.1093/bioinformatics/bts485

PubMed Abstract | CrossRef Full Text | Google Scholar

Wong, N., and Wang, X. (2015). miRDB: an online resource for microRNA target prediction and functional annotations. Nucleic Acids Res. 43, D146–D152. doi: 10.1093/nar/gku1104

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Q., Chen, S., Yu, S., Qin, J., Zhang, J., Cheng, Q., et al. (2016). Neuroprotective effects of pyrroloquinoline quinone against rotenone injury in primary cultured midbrain neurons and in a rat model of Parkinson’s disease. Neuropharmacology 108, 238–251. doi: 10.1016/j.neuropharm.2016.04.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Q., Zhang, J., Jiang, C., Qin, J., Ke, K., and Ding, F. (2014). Involvement of ERK1/2 pathway in neuroprotective effects of pyrroloquinoline quinine against rotenone-induced SH-SY5Y cell injury. Neuroscience 270, 183–191. doi: 10.1016/j.neuroscience.2014.04.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., and Cheng, Y. (2014). miR-16-1 promotes the aberrant alpha-synuclein accumulation in Parkinson disease via targeting heat shock protein 70. ScientificWorldJournal 2014:938348. doi: 10.1155/2014/938348

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Parkinson’s disease, bioinformatics, biomarkers, genes, miRNA, NAFD pathway

Citation: Chi J, Xie Q, Jia J, Liu X, Sun J, Deng Y and Yi L (2018) Integrated Analysis and Identification of Novel Biomarkers in Parkinson’s Disease. Front. Aging Neurosci. 10:178. doi: 10.3389/fnagi.2018.00178

Received: 16 January 2018; Accepted: 24 May 2018;
Published: 18 June 2018.

Edited by:

Nibaldo C. Inestrosa, Pontificia Universidad Católica de Chile, Chile

Reviewed by:

Ming Li, Lieber Institute for Brain Development, United States
Richard Camicioli, University of Alberta, Canada

Copyright © 2018 Chi, Xie, Jia, Liu, Sun, Deng and Yi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Li Yi, yilitj@hotmail.com