Background {#Sec1}
==========

Parkinson's disease (PD) is one of the most common neurodegenerative diseases. The main symptoms of PD include shaking, bradypragia and postural instability \[[@CR1]\]. PD is characterized by the loss of dopaminergic neurons in the substantia nigra (SN) in brain, which involves the increase of microglia, and the presence of Lewy bodies \[[@CR2]\].

Many studies have indicated that PD may be the result of a combination of genetic and environmental factors, such as exposure to pesticides, metals, solvents, and other toxicants \[[@CR3]\]. Approximately 15% of individuals affected by PD have a family member with the disease \[[@CR4]\] and 5--10% of people with PD are affected by a single specific gene mutation \[[@CR5]\]. Previous studies have indicated that mutations in several specific genes cause PD, including genes encoding alpha-synuclein (SNCA), leucine-rich repeat kinase 2 (LRRK2), glucocerebrosidase (GBA), parkin (PRKN), PTEN-induced putative kinase 1 (PINK1), parkinson disease protein 7 (PARK7), vacuolar protein sorting-associated protein 35 (VPS35), eukaryotic translation initiation factor 4 gamma 1 (EIF4G1), dnaJ heat shock protein family (Hsp40) member C13 (DNAJC13) and coiled-coil-helix-coiled-coil-helix domain containing 2 (CHCHD2) \[[@CR6]\]. The SNCA gene plays an important role in PD, because the encoded protein is the main component of Lewy bodies, which accumulate in the brains of people with PD \[[@CR7]\]. Mutations in the SNCA gene have been found in different groups with sporadic (non-familial) PD and familial PD \[[@CR5]\]. Mutations in the LRRK2 gene, which encodes a protein called dardarin, are associated with many familial and sporadic PD patients \[[@CR8]\]. A recent study showed that carriers of the G2385R variant in the LRRK2 gene showed more of a tendency towards fatigue than non-carriers in PD patients \[[@CR9]\]. A mutation in GBA is proposed to be the greatest genetic risk of PD through its effects in increasing the levels of SNCA \[[@CR10]\]. The GBA variant E326K (rs2230288) was significantly more frequent in PD patients compared with controls, indicating that this variant is a susceptibility allele for PD \[[@CR10]\]. Mutations in PINK1, PRKN, and PARK7 genes may cause mitochondrial dysfunction, which are observed in PD \[[@CR11]\]. Some evidence has shown that low concentrations of urate in the blood serum increase the risk of PD \[[@CR12]\].

These above studies improve our understanding of the molecular mechanism in PD. However, discordance among those studie make the combination of results from many similar studies difficult. However, meta-analysis uses statistical techniques to get combining information from multiple existing studies can increase the power and generalizability of results in PD \[[@CR13]\]. A limited number of meta studies has been conducted on gene expression profiles to identify novel biomarkers for PD. However, different microarray data combination, different sizes of each individual study combination, and different statistical methods are among the factors contributed to the different results. For instance, Elisa Mariani et al. \[[@CR14]\] used TRAM Software to conduct a meta-analysis of PD transcriptome data. Chang D et al. \[[@CR15]\] carried out a meta-analysis with a recent study of over 13,000 PD cases and 95,000 controls with estimating the h2 value for PD, and then identified 17 new PD risk loci. Wang Q et al. \[[@CR16]\] performed a meta-analysis with 9 microarray datasets of PD studies in every brain region using the RankProd method \[[@CR17]\].

Our meta-analysis was conducted on PD patients and healthy controls microarray datasets obtained from the SN brain region and blood using Vote counting generic ways. In this study, five independent SN microarrays (GSE7621, GSE8397-GPL96, GSE8397-GPL97, GSE20163, and GSE20141) from PD patients and healthy controls were integrated and analyzed to screen common differentially expressed genes (DEGs). We also performed a meta-analysis of three independent blood microarrays (GSE99039, GSE6613, and GSE72267) from PD patients and healthy controls to screen common DEGs.

MicroRNAs (miRNAs) are small non-coding RNA molecules (\~ 22 nt long) that are found in plants, animals and some viruses and have been proven to play crucial roles in gene expression. Long non-coding RNAs (lncRNAs), longer than 200 nt, can control gene expression by interacting with the miRNA pathways \[[@CR18]\]. LncRNAs binds to miRNAs through miRNA response elements or binding sites to regulate miRNA target gene expression \[[@CR18]\]. Several studies have revealed that non-coding RNAs, such as miRNAs and lncRNAs, have been implicated in PD pathogenesis \[[@CR13]\].

Finally, we examined the regulatory network involving genes, miRNAs, transcription factors (TFs) and lncRNA in PD progression to better understand the molecular mechanisms involved in this disease. In our study, we use some same bioinformatics software and similar methodological workflow as the recently published paper \[[@CR13]\].

Methods {#Sec2}
=======

Our overarching goal was to identify candidate biomarkers of PD (the focus of this study) using the follow workflow (Fig. [1](#Fig1){ref-type="fig"}).Fig. 1Workflow of the methodology used in our study

Data collection {#Sec3}
---------------

We used "Parkinson's disease" as keywords to search for genome-wide expression studies in NCBI-GEO ([http://www.ncbi.nlm.nih.gov/geo/](http://www.ncbi.nlm.nih.gov/geo)) and EMBL-EBI ArrayExpress databases ([https://www.ebi.ac.uk/arrayexpress/](https://www.ebi.ac.uk/arrayexpress)). Only original experimental studies that screened for different genes between PD and healthy humans were the first choice for inclusion. Additional inclusion criteria included the following: (1) the study type was expression profiling by array; (2) Studies which comprised of cell intensity file (CEL) raw files were available; and (3) studies about blood or brain SN in PD were used in the analysis. With the inclusion criteria, three datasets from blood and five datasets from SN in PD were screened. We then performed meta-analysis of three datasets from blood studies (GSE99039, GSE6613, and GSE72267) and five datasets from SN studies (GSE7621, GSE8397-GPL96, GSE8397-GPL97, GSE20163, and GSE20141). Details of the microarray datasets are provided in Table [1](#Tab1){ref-type="table"}.Table 1Datasets used in the meta-analysisTissueGEO AccessionSample Size(control/PD)PlatformPMIDBloodGSE99039PD = 205; HC = 233GPL570: Affymetrix Human Genome U133 Plus 2.0 Array28,916,538 \[[@CR90]\]GSE6613PD = 50; HC = 22GPL96: Affymetrix Human Genome U133A Array17,215,369 \[[@CR91]\]GSE72267PD = 40; HC = 19GPL571: Affymetrix Human Genome U133A 2.0 Array26,510,930 \[[@CR92]\]SNGSE7621PD = 16; HC = 9GPL570: Affymetrix Human Genome U133 Plus 2.0 Array17,571,925 \[[@CR93]\]GSE8397-GPL96PD = 24; HC = 15GPL96: Affymetrix Human Genome U133A Array;16,344,956 \[[@CR94]\]GSE8397-GPL97PD = 24; HC = 15GPL97: Affymetrix Human Genome U133B ArrayGSE20141PD = 10; HC = 8GPL570: Affymetrix Human Genome U133 Plus 2.0 Array20,926,834 \[[@CR95]\]GSE20163PD = 8; HC = 9GPL96: Affymetrix Human Genome U133A Array20,926,834 \[[@CR95]\]*PD* Parkinson's disease, *HC* healthy control

Searches were executed up to August 2017.

Inclusion criteria for DEGs {#Sec4}
---------------------------

GEO2R is an interactive web tool that is based on GEOquery and limma R packages from the Bioconductor project that is used to identify DEGs by comparing different groups of GEO Series samples \[[@CR19]\]. In this study, each individual dataset was processed using the GEO2R online software to compare healthy control samples and PD patient samples to screen DEGs. The t test and Benjamini and Hochberg method were used to calculate the *P* values and FDR, respectively \[[@CR20]\]. Genes that were differentially expressed between PD and healthy controls were screened with the threshold value of *P* \< 0.05 and \|log fold change\| \> 0.1. The probes with no gene annotation or matched multiple gene symbols were removed, and when multiple probes matched to the same gene symbol, we selected the probe with the highest \|log fold change\| value.

Meta-analysis of DEGs of gene expression microarray datasets {#Sec5}
------------------------------------------------------------

Here, Meta-analysis was performed on the three datasets from blood studies and five datasets from brain SN independently using Vote counting generic ways of combining information \[[@CR21]\]. The results were visualized using a Venn diagram that were plotted using the OmicShare tools, a free online platform for data analysis (<http://www.omicshare.com/tools>).

Next, we mainly used GCBI online software to identify common DEGs associated with PD as PD-specific genes, and to identify other neurodegenerative diseases as other neurodegenerative diseases related genes.

Functional enrichment analysis {#Sec6}
------------------------------

BiNGO is a tool to determine the Gene Ontology (GO) categories that are statistically overrepresented in a set of genes or a subgraph of a biological network \[[@CR22]\]. BiNGO was used to analyze GO enrichment of DEGs in blood tissues and in brain SN of PD. GO enrichment with *p* \< 0.05 and false discovery rate (FDR) \< 0.05 was regarded as statistically significant.

Protein-protein interaction (PPI) network construction for common DEGs {#Sec7}
----------------------------------------------------------------------

The selected common genes from SN tissue and blood were subjected to STRING v.10.5 database \[[@CR23]\] analysis to construct PPI networks. The threshold value was a score of 0.4. PPI networks were visualized by Cytoscape software v. 3.4.0 \[[@CR24]\] and analyzed using the Network Analyzer tool based on degree. The degree indicates the number of interactions of a particular protein. The size of the node is proportional to the degree in the interaction network. The more the nodes connected to the node, the larger the node will be. Large nodes indicate bigger degrees. In this network, we selected some genes that are with a high node degree. At last, the top three nodes with degree values above the average network degree value were identified.

Biological regulatory interaction and networks for common PD-specific genes from our study {#Sec8}
------------------------------------------------------------------------------------------

Experimentally validated miRNA and common PD-specific gene interactions were analyzed using DIANA-Tarbase \[[@CR25]\]. We also examined lncRNA and miRNA interactions using DIANA-LncBase Experimental v.2 \[[@CR26]\]. This database contains more than 70,000 low and high-throughput, (in)direct miRNA: lncRNA experimentally supported interactions, derived from manually curated publications and the analysis of 153 AGO CLIP-Seq libraries \[[@CR26]\]. Experimentally validated (prediction score ≥ 0.90) lncRNAs were selected in our study.

To study molecular regulator interactions, we built a regulatory network comprising common PD-specific genes, TFs associated with the genes, miRNAs associated with the genes and TFs associated with the miRNAs. Information on TFs binding site that were associated with the genes was obtained from TRANSFAC \[[@CR27]\] based on the Match™ algorithm. The TRANSFAC database contains data from a wide variety of eukaryotic organisms, ranging from human to yeast and comprising data on transcription factors, their target genes and regulatory binding sites \[[@CR27]\]. In our study, the search algorithm uses two score values, the matrix similarity score (MSS) and the core similarity score (CSS), to estimate the result. Finally, the inclusion criteria was score = 1 for the MSS and CSS. The gene-miRNA interaction information was analyzed using DIANA-Tarbase \[[@CR25]\]. Regulatory information on the TFs associated with these miRNAs was obtained from the TransmiR database \[[@CR28]\]. We then generated the network using Cytoscape software v. 3.4.0 \[[@CR24]\].

PD-specific and other neurodegenerative diseases related miRNAs identified from miRNAs associated with the PD-specific genes were filtered from miRdSNP v.11.03 online database \[[@CR29]\].

The analysis of lncRNAs and mTFs associated with these PD-specific miRNAs related to the PD-specific genes {#Sec9}
----------------------------------------------------------------------------------------------------------

To identify possible lncRNA-mediated regulation of the PD-specific miRNAs associated with common PD-specific genes of our study, we collected regulatory information in the lncbase module of DIANA-LncBase v.2 \[[@CR26]\]. The experimentally validated lncRNAs (score ≥ 0.90) were selected in our study. Regulatory information on the TFs that were associated with these miRNAs was obtained from the TransmiR v.1.2 database \[[@CR26]\]. We then got a regulatory pattern involving the TFs, lncRNAs and common PD-specific genes associated with the PD-specific miRNAs.

Regulatory analysis of the common non-PD-specific genes {#Sec10}
-------------------------------------------------------

To study the molecular regulator interactions with the common non-PD-specific genes, we perform regulatory analysis including TFs associated with the genes, miRNAs associated with the genes and TFs associated with the miRNAs. Information on the TFs that were associated with the genes was obtained from TRANSFAC \[[@CR27]\] and the inclusion criteria was score = 1for the MSS and CSS. The gene-miRNA interaction information was analyzed using DIANA-Tarbase v.7.0 \[[@CR25]\], miRWalk database v.3.0 \[[@CR20]\] and TargetScan human v.7.1 \[[@CR30]\]. Regulatory information on the TFs that were associated with these miRNAs was obtained from the TransmiR v.1.2 database \[[@CR28]\].

PD-specific miRNAs identified from miRNAs associated with the common non-PD-specific genes were filtered from miRdSNP v.11.03 online database \[[@CR29]\].

SNP analysis of the common PD-specific and other neurodegenerative disease-specific genes {#Sec11}
-----------------------------------------------------------------------------------------

To identify PD-associated SNPs, the common genes in our study were subjected to SNP analysis. SNPs corresponding to these genes were identified from the MirSNP online database ([http://bioinfo.bjmu.edu.cn/mirsnp/search/](http://bioinfo.bjmu.edu.cn/mirsnp/search)) \[[@CR31]\]. We obtained a large number of SNPs related to the genes. To identify the PD-specific or neurodegenerative disease-specific SNPs from this large number of SNPs, we screened the SNPs in miRdSNP v.11.03 \[[@CR29]\] and LincSNP v.2.0 \[[@CR32]\]. Chromosome locus and allele gene information corresponding to each of the SNPs were searched using the dbSNP database(<https://www.ncbi.nlm.nih.gov/snp/?term=>).

Results {#Sec12}
=======

Meta-analysis of DEGs and PD-specific or other neurodegenerative diseases related DEGs identification {#Sec13}
-----------------------------------------------------------------------------------------------------

To identify DEGs between PD and healthy controls, we used three datasets from blood and five datasets from the brain SN region to perform meta-analyses and detected common genes across these datasets. The analyses revealed that 28 downregulated genes and 8 upregulated genes were differentially expressed in the three blood studies (Fig. [2a, b](#Fig2){ref-type="fig"} and Table [2](#Tab2){ref-type="table"}). Out of the total 36 DEGs in the blood studies, 19 were previously known to be associated with several other neurodegenerative diseases, such as Alzheimer's disease and four common genes (HSPA6, MAP2K6, SRPK2 and NOL7) were also previously associated with PD (named PD-specific genes)(identified using GCBI online software). We also identified 11 downregulated and 6 upregulated genes as differentially expressed in the five brain SN studies (Fig. [2c, d](#Fig2){ref-type="fig"} and Table [3](#Tab3){ref-type="table"}). Among the 17 total DEGs in the brain SN studies, nine genes have been previously associated with neurodegenerative diseases, and one gene (SNCA) was identified as PD-specific genes(identified using GCBI online software). Overall, the number of downregulated genes was greater than the number of upregulated genes.Fig. 2The number of common genes obtained from blood and substantia nigra (SN) expression profiling datasets visualized by a Venn diagram. **a** The number of downregulated genes in expression profiling datasets from blood; **b** The number of upregulated genes in expression profiling datasets from blood; **c** The number of downregulated genes in expression profiling datasets from SN; and **d** The number of upregulated genes in expression profiling datasets from SNTable 2Common differentially expressed genes identified from blood in PD patients and healthy controls and that PD-specific and other neurodegenerative diseases related genes identificationExpressionGene symbolGene nameCorresponding neurodegenerative disease other than PDDown-regulatedABCA1^a^ATP binding cassette subfamily A member 1Alzheimer's diseaseABHD5abhydrolase domain containing 5ADGRG3adhesion G protein-coupled receptor G3AKAP13^a^A-kinase anchoring protein 13Alzheimer's diseaseAPMAP^a^adipocyte plasma membrane associated proteinMultiple sclerosis; Alzheimer's diseaseARG1arginase 1BAZ1A^a^bromodomain adjacent to zinc finger domain 1AHuntington diseaseBMX^a^BMX non-receptor tyrosine kinaseAlzheimer's diseaseCDKL5^a^cyclin dependent kinase like 5CEBPD^a^CCAAT/enhancer binding protein deltaAlzheimer's diseaseCSF2RA^a^colony stimulating factor 2 receptor alpha subunitAlzheimer's diseaseCTBP2^a^C-terminal binding protein 2Alzheimer's diseaseFAM120A^a^family with sequence similarity 120AMultiple sclerosisHLA-C^a^major histocompatibility complex, class I, CMultiple sclerosisHSPA6^b^heat shock protein family A (Hsp70) member 6IRS2insulin receptor substrate 2LILRB1leukocyte immunoglobulin like receptor B1LRRFIP1^a^LRR binding FLII interacting protein 1Multiple sclerosis; Huntington diseaseMAP2K6^b^mitogen-activated protein kinase kinase 6MAP2K7mitogen-activated protein kinase kinase 7MGAMmaltase-glucoamylaseNCAM1neural cell adhesion molecule 1NCOA3nuclear receptor coactivator 3PML^a^promyelocytic leukemiaPolyglutamine diseasesSRPK2^b^SRSF protein kinase 2SUPT20H^a^SPT20 homolog, SAGA complex componentCerebellar Purkinje cell degenerationTHOC2^a^THO complex 2Multiple sclerosisTMX4^a^thioredoxin related transmembrane protein 4Motor neuron diseaseUp-regulatedATMATM serine/threonine kinaseBCL2^a^BCL2, apoptosis regulatorAlzheimer's diseaseFAM102A^a^family with sequence similarity 102 member ANeurodegenerative diseaseLRRN3leucine rich repeat neuronal 3NOL7^b^nucleolar protein 7TCF3transcription factor 3TP73-AS1TP73 antisense RNA 1YME1L1^a^YME1 like 1 ATPaseAlzheimer's disease^a^Previously associated with several neurodegenerative diseases (identified using GCBI online software), but not PD^b^Previously associated with PD (identified using GCBI online software)Table 3Common differentially expressed genes identified from substantia nigra in PD patients and healthy controls, and PD-specific and other neurodegenerative diseases related genes identificationExpressionGene symbolGene nameCorresponding neurodegenerative disease other than PDDown-regulatedCARHSP1^a^calcium regulated heat stable protein 1Multiple sclerosisGART^a^phosphoribosylglycinamide formyltransferase phosphoribosylglycinamide synthetase, phosphoribosylaminoimidazole synthetaseAlzheimer's diseaseMUC4mucin 4, cell surface associatedNINnineinNRF1^a^nuclear respiratory factor 1Alzheimer's diseasePRELP^a^proline and arginine rich end leucine rich repeat proteinNervous system disease; Neurodegenerative diseaseRGS12regulator of G-protein signaling 12RNF130^a^ring finger protein 130Motor neuron diseaseSNAP23^a^synaptosome associated protein 23Nervous system diseaseSNTB2syntrophin beta 2TBX1T-box 1Up-regulatedACSL6^a^acyl-CoA synthetase long-chain family member 6Nervous system diseaseATP5S^a^ATP synthase, H+ transporting, mitochondrial Fo complex subunit s (factor B)Alzheimer's diseaseCADPScalcium dependent secretion activatorDCLK1^a^doublecortin like kinase 1Alzheimer's diseasePPFIA2PTPRF interacting protein alpha 2SNCA^b^synuclein alpha^a^Previously associated with several neurodegenerative diseases, other than PD (identified using GCBI online software)^b^Previously associated with PD (PD-specific) (identified using GCBI online software)

GO functional enrichment analysis of all common DEGs {#Sec14}
----------------------------------------------------

The common DEGs from blood and brain SN were then subjected to enrichment analysis in BiNGO. The common DEGs in blood were enriched in biological processes (BP) of protein phosphorylation, DNA damage-induced protein phosphorylation, phosphorylation, response to gamma radiation and B cell lineage commitment, among other processes (Additional file [1](#MOESM1){ref-type="media"}: Table S1 and Fig. [3](#Fig3){ref-type="fig"}). DEGs from brain SN were enriched in BP including cellular nitrogen compound biosynthetic process, negative regulation of G-protein coupled receptor protein signaling pathway, pigment biosynthetic process, regulation of acyl-CoA biosynthetic process, negative regulation of neurotransmitter uptake and negative regulation of catecholamine uptake involved in synaptic transmission (Additional file [1](#MOESM1){ref-type="media"}: Table S2 and Fig. [4](#Fig4){ref-type="fig"}). In blood, these processes such as protein amino acid phosphorylation (9 genes), regulation of neuron maturation (1 gene), positive regulation of neuron maturation (1 gene) are highly related to neuron. SRPK2 and MAP2K6 are found to be involved in PD \[[@CR33], [@CR34]\]. BCL2 was reported in Alzheimer's Disease \[[@CR35]\]. In brain SN, negative regulation of neurotransmitter uptake (1 gene) and negative regulation of catecholamine uptake involved in synaptic transmission (1 gene) are highly related to neuron. SNCA is reported to be in PD \[[@CR36]\]. All these information validate our finding regarding the association of these genes in PD.Fig. 3GO biological processes network of diferentially expressed genes in blood of Parkinson's Disease from BiNGO software. Large nodes indicate more genes involved in. Yellow nodes: *P*-value \< 0.05 and FDR \< 0.05Fig. 4GO biological processes network of diferentially expressed genes in substantia nigra of Parkinson's Disease from BiNGO software. Large nodes indicate more genes involved in. Yellow nodes: P-value \< 0.05 and FDR \< 0.05

PPI network construction for common DEGs {#Sec15}
----------------------------------------

PPI networks were constructed using STRING v.10.5 database and then visualized and analyzed by Cytoscape software v.3.4.0. The PPI information obtained from STRING online software for common DEGs in blood is demonstrated in Additional file [1](#MOESM1){ref-type="media"}: Table S3 and the PPI network is shown in Fig. [5a](#Fig5){ref-type="fig"}. Here, the 29 genes exhibited a wide degree distribution with the highest degree of 18 and lowest degree of 1.The average degree value was found to be 6.14. Finally, the top three hub genes(i.e. top 10% of the total nodes) with higher degree were chosen, including ATM, MGAM and BCL2 genes. Out of the three hubs, one hub (BCL2) was already found to be related to Alzheimer's Disease \[[@CR35]\], which was further studied for their association in human PD. The PPI information obtained from STRING online software for common DEGs in brain SN is exhibited in Additional file [1](#MOESM1){ref-type="media"}: Table S4 and the PPI network is shown in Fig. [5b](#Fig5){ref-type="fig"}. We found that highest degree value was 3 and the lowest was 1 with an average of 1.69. In this analysis, the top three hub genes with higher degree included GART, SNCA and NIN genes. SNCA hub node was already found to be associated with human PD \[[@CR36]\], and GART hub node has been verified in other neurodegenerative diseases such as Alzheimer's Disease \[[@CR37]\], but not in PD.Fig. 5PPI networks obtained from Cytoscape software. **a** PPI network for common genes from blood; **b** PPI network for common genes from the brain substantia nigra region. Red nodes represent upregulated genes, and green nodes represent downregulated genes. Large nodes indicate bigger degrees

Regulatory miRNAs associated with the common PD-specific genes and the analysis of lncRNAs associated with these miRNAs {#Sec16}
-----------------------------------------------------------------------------------------------------------------------

To investigate the regulation of the PD-specific genes identified in our study by non-coding RNAs, miRNA associated with these genes were first retrieved from DIANA-Tarbase v.7.0 database, and then lncRNAs associated with these miRNAs were selected from DIANA-LncBase Experimentally v.2 (prediction score ≥ 0.90). The biological regulatory interaction of PD-specific genes identified in blood is exhibited in Table [4](#Tab4){ref-type="table"}, which shows that four PD-specific genes were associated with several non-coding RNAs. Out of the miRNAs, 13 were previously known to be associated with several other neurodegenerative diseases, such as Alzheimer's disease and nine were also previously associated with PD. The biological regulatory interaction of PD-specific genes identified in brain SN is showed in Table [5](#Tab5){ref-type="table"}, which demonstrates that one PD-specific gene was associated with several non-coding RNAs. Out of the miRNAs, one was previously known to be associated with other neurodegenerative diseases and five were also previously associated with PD. Out of these five, four PD-specific genes (HSPA6, SRPK2, NOL7 and SNCA) were associated with miRNAs that were previously associated with PD. Moreover, most of these miRNAs associated with PD-specific genes were in turn regulated by lncRNAs. Among the miRNAs, 14 has been previously associated with neurodegenerative diseases, and 14 were previously associated with PD.Table 4Regulatory miRNAs associated with the common PD-specific genes identified in blood and the analysis of lncRNAs associated with these miRNAsPD specific genesmiRNAs associated with geneslncRNAs associated with the miRNAsHSPA6hsa-miR-204-3p^b^LINC00999hsa-miR-548o-3pCASC7, GABPB1-AS1, NEAT1, XISTMAP2K6hsa-miR-33a-5p^a^KCNQ1OT1, MCF2L-AS1hsa-miR-590-3p^a^CASC7, CTA-292E10.9, CTB-89H12.4, HCG11, LINC00657, LOC100507577, NEAT1, OIP5-AS1, OTUD6B-AS1, RP11-834C11.4, XISThsa-miR-145-5p^a^KCNQ1OT1,TUG1hsa-miR-425-5pC1orf132, CTD-3025 N20.3, KCNQ1OT1, RP11-15H20.7, SNHG14, TTTY15, ZNRD1-AS1hsa-miR-1306-5pKCNQ1OT1, SENP3-EIF4A1hsa-miR-148a-3p^a^CASC7, KCNQ1OT1, OIP5-AS1, SNHG14hsa-miR-130a-3p^a^CASC7, H19, SNHG14hsa-miR-148b-3p^a^CASC7, OIP5-AS1, SLMO2-ATP5E, SNHG14SRPK2hsa-miR-3200-3pBLOC1S5-TXNDC5, KCNQ1OT1, XISThsa-miR-1292hsa-miR-155-5pCTD-2561 J22.5, RP11-175O19.4, XISThsa-miR-15b-3p^a^hsa-miR-17-3p^b^RP3-323A16.1, XISThsa-miR-181a-5p^b^AC000403.4, CASC7, CTB-89H12.4, KCNQ1OT1, LINC00506, N4BP2L2-IT2, RP11-10E18.7, RP11-1134I14.8, RP11-147 L13.14, RP11-314B1.2, RP11-361F15.2, RP11-707A18.1, RP1-309I22.2, ZNF883, ZNRD1-AS1hsa-miR-181b-5p^b^CASC7, CTB-89H12.4, KCNQ1OT1, RP11-1134I14.8, XISThsa-miR-181c-5p^b^CTB-89H12.4, KCNIP4-IT1, KCNQ1OT1, RP11-1134I14.8hsa-miR-181d^b^hsa-miR-183-3p^a^hsa-miR-1976HNRNPUL2-BSCL2, KCNQ1OT1, LOC100129917, NEAT1, TSIXhsa-miR-19a-3p^b^FAM201A, H19, KCNA3, KCNQ1OT1, LINC00094, RP11-337C18.8, RP11-523G9.3, SNHG14hsa-miR-19b-3p^b^CASC7, FAM201A, H19, KCNA3, KCNQ1OT1, LINC00094, RP11-337C18.8, RP11-523G9.3, SNHG14hsa-miR-21-5p^a^hsa-miR-320aALMS1-IT1, CASC7, CTB-36H16.2, CTB-89H12.4, KCNIP4-IT1, LINC00663, LMCD1-AS1, MALAT1, NEAT1, RP11-145P16.3, XISThsa-miR-3685KCNQ1OT1hsa-miR-3689a-3pTTTY15hsa-miR-4518hsa-miR-944NOL7hsa-miR-199b-3p^a^CTB-89H12.4, ERVK3--1, XISThsa-miR-199a-3p^b^CTB-89H12.4, ERVK3--1, XISThsa-miR-328hsa-miR-129-5p^a^CASC7, ERVK3--1, KCNQ1OT1, MALAT1, NEAT1,hsa-miR-374a-5p^a^CTC-444 N24.11, CTD-2561 J22.5, RP11-613D13.5, TRG-AS1, XIST, ZNRD1-AS1hsa-miR-744-5pFLJ16779,hsa-miR-374b-5p^a^CTA-292E10.9, CTC-444 N24.11, OIP5-AS1, RP11-221 J22.1, RP11-38P22.2, XIST,hsa-miR-548o-3pCASC7, GABPB1-AS1, NEAT1, XIST,^a^Previously associated with several neurodegenerative diseases, but not PD^b^Previously associated with PDTable 5Regulatory miRNAs associated with the common PD-specific genes identified in substantia nigra and the analysis of lncRNAs associated with these miRNAsGenes associated with PDmiRNAs associated with geneslncRNAs associated with the miRNAsSNCAhsa-miR-93-3p^b^AC012065.7,KCNQ1OT1, LINC00342hsa-miR-153^b^hsa-miR-23b-3p^b^CASC7,CTC-459F4.3, KCNQ1OT1, RP11-215G15.5, SNHG14, TOB1-AS1, XIST, ZNRD1-AS1hsa-miR-500a-5pSNHG22hsa-miR-34a-5p^b^AC004951.6, AC092535.3, KCNQ1OT1, LINC00662, PCBP2-OT1, RP11-693 J15.5hsa-miR-29a-3p^a^AC005154.6, H19, KCNQ1OT1, LINC00674, MIR4697HG, NEAT1, RP11-314B1.2, RP11-582E3.6, RP4-630A11.3, THUMPD3-AS1, TTTY15, TUG1hsa-miR-7-5p^b^AC005154.6, DLX6-AS1, KCNQ1OT1, LINC01233, LINC01314, RP11-679B19.1, XISThsa-miR-181a-2-3pKCNQ1OT1, NEAT1^a^Previously associated with several neurodegenerative diseases, but not PD^b^Previously associated with PD

mTF-miRNA-gene-gTF regulatory network construction of PD-specific genes {#Sec17}
-----------------------------------------------------------------------

To further study the regulatory mechanism of these common PD-specific genes (the genes from Tables [2](#Tab2){ref-type="table"} and [3](#Tab3){ref-type="table"} that were defined as PD-specific by authors) in PD pathogenesis, we constructed a regulatory network comprising common PD-specific genes, the TFs associated with these genes (gTFs), miRNAs associated with these genes and the TFs associated with these miRNAs (mTFs). Additional file [1](#MOESM1){ref-type="media"}: Table S5 and Fig. [6](#Fig6){ref-type="fig"} show the mTF-miRNA-gene-gTF regulatory networks identified in blood. Additional file [1](#MOESM1){ref-type="media"}: Table S6 and Fig. [7](#Fig7){ref-type="fig"} show the mTF-miRNA-gene-gTF regulatory networks identified in brain SN.Fig. 6The mTF-miRNA-gene-gTF regulatory network of PD-specific genes identified in blood datasets obtained from Cytoscape software. The diamond-shaped magenta nodes represent miRNAs, the triangle-shaped magenta nodes represent transcription factors (TFs) associated with these miRNAs (mTFs), and the round rectangle-shaped red nodes represent the PD-specific genes. The V-shaped red nodes represent the TFs associated with these genes (gTFs)Fig. 7The mTF-miRNA-gene-gTF regulatory network of PD-specific genes identified in substantia nigra datasets obtained from Cytoscape software. The diamond-shaped magenta nodes represent miRNAs, the triangle-shaped magenta nodes represent transcription factors (TFs) associated with these miRNAs (mTFs), and the round rectangle-shaped red nodes represent the PD-specific genes. The V-shaped red nodes represent the TFs associated with these genes (gTFs)

IncRNAs and mTFs associated with these PD-specific miRNAs related to the PD-specific genes {#Sec18}
------------------------------------------------------------------------------------------

The 14 PD-specific miRNAs shown in Tables [4](#Tab4){ref-type="table"} and [5](#Tab5){ref-type="table"} (hsa-miR-204-3p, hsa-miR-17-3p, hsa-miR-181a-5p, hsa-miR-181b-5p, hsa-miR-181c-5p, hsa-miR-181d, hsa-miR-19a-3p, hsa-miR-19b-3p, hsa-miR-93-3p, hsa-miR-153, hsa-miR-23b-3p, hsa-miR-34a-5p, hsa-miR-9-5p and hsa-miR-7-5p) were further analyzed in the DIANA-LncBase Experimental v.2 database to examine the regulatory information of miRNAs-lncRNAs (Table [6](#Tab6){ref-type="table"}) and in the TransmiR v.1.2 database to examine the regulatory information of miRNAs-mTFs. The 14 miRNAs were found to be associated with 45 lncRNAs. These 14 miRNAs control four PD-specific genes in our study, and these miRNAs are in turn regulated by 31mTFs (Table [6](#Tab6){ref-type="table"}). By combining all the regulatory information obtained from our results, we constructed a regulatory pattern involving the PD-specific miRNAs, lncRNAs, PD-specific genes and mTFs.Table 6The lncRNA-mediated PD-specific miRNAs associated the PD-specific mRNAs identified in our study regulatory networkTissuePD-specific miRNAlncRNAs associated the miRNAsthe PD-specific mRNAs associated the miRNATFs associated with the miRNAsBloodhsa-miR-204-3pLINC00999HSPA6hsa-miR-17-3pRP3-323A16.1, XISTSRPK2Activation:CCND1,E2F1,MYC, MYCN, NKX2--5, TLX1, TLX3, TNF, ERS1, STAT5;\
Repression: NFKB1, SPI1; Neutral: MAPK14hsa-miR-181a-5pAC000403.4, CASC7, CTB-89H12.4, KCNQ1OT1, LINC00506, N4BP2L2-IT2, RP11-10E18.7, RP11-1134I14.8, RP11-147 L13.14, RP11-314B1.2, RP11-361F15.2, RP11-707A18.1, RP1-309I22.2, ZNF883, ZNRD1-AS1SRPK2Neutral: TP63hsa-miR-181b-5pCASC7, CTB-89H12.4, KCNQ1OT1, RP11-1134I14.8, XISTSRPK2hsa-miR-181c-5pCTB-89H12.4, KCNIP4-IT1, KCNQ1OT1, RP11-1134I14.8SRPK2Activation: AKT1hsa-miR-181dSRPK2hsa-miR-19a-3pFAM201A, H19, KCNA3, KCNQ1OT1, LINC00094, RP11-337C18.8, RP11-523G9.3, SNHG14SRPK2Activation: E2F1, MYC, MYCN, NKX2--5, TLX1, TLX3, ERS1, STAT5; Repression:SPI1; Neutral:PTEN, MAPK14hsa-miR-19b-3pCASC7, FAM201A, H19, KCNA3, KCNQ1OT1, LINC00094, RP11-337C18.8, RP11-523G9.3, SNHG14SRPK2Activation: E2F1, MYC, MYCN, NKX2--5, TLX1, TLX3, ERS1;\
Neutral: MAPK14hsa-miR-199a-3pCTB-89H12.4, ERVK3--1, XISTNOL7Neutral: TWIST1Substantia nigrahsa-miR-93-3pAC012065.7,KCNQ1OT1, LINC00342SNCARepression: MYC;\
Neutral: E2F1hsa-miR-153SNCAhsa-miR-23b-3pCASC7,CTC-459F4.3, KCNQ1OT1, RP11-215G15.5, SNHG14, TOB1-AS1, XIST, ZNRD1-AS1SNCARepression: MYC;\
Neutral: ESR2hsa-miR-34a-5pAC005154.6, H19, KCNQ1OT1, LINC00674, MIR4697HG, NEAT1, RP11-314B1.2, RP11-582E3.6, RP4-630A11.3, THUMPD3-AS1, TTTY15, TUG1SNCAActivation: HMGA1, AP-1;\
Repression: MYC, NFKB1, YY1, IL-4, PDGF-B, TGFB1;\
Neutral: CEBPAhsa-miR-7-5pCTB-89H12.4, KCNQ1OT1, RP11-273G15.2, RP11-314B1.2, RP11-793H13.8, SNHG14, TSNAX-DISC1, TUG1SNCAActivation: NFKB1, IL1B, TLR2, TLR4, TLR7,\
TLR8, TNF, MYC, NFKB1;\
Repression: TLX

mTF-miRNA-gene-gTF regulatory analysis of the common non-PD-specific genes {#Sec19}
--------------------------------------------------------------------------

To study the regulatory interaction of the common non-PD-specific genes, we got mTF-miRNA-gene-gTF regulatory patterns from blood (Additional file [1](#MOESM1){ref-type="media"}: Table S7) and brain SN (Additional file [1](#MOESM1){ref-type="media"}: Table S8) respectively. We obtained 359 gTFs (score = 1) that were related to 26 common DEGs from blood datasets. The TCF3, BAZ1A and THOC2 genes were regulated by the highest numbers of gTFs (105, 95 and 95 gTFs, respectively). Previous reports showed TCF3 involvement in PD \[[@CR38]\]. TCF3, a transcriptional repressor \[[@CR39]\], is counteracted at multiple levels by Wnt signaling \[[@CR40]\]. Alteration in normal Wnt signaling pathway has been already found in gene expression studies of blood of PD patients \[[@CR40]\]. It was found that BAZ1A depletion influence the expression of genes important for nervous system development and function \[[@CR41]\]. THOC2, an abundant protein in the developing and mature human \[[@CR42], [@CR43]\] and adult mouse brains \[[@CR44]\]. THOC2 depletion has been already reported to stimulate neurite outgrowth in primary rat hippocampal neurons \[[@CR44]\]. It has been identified also to be involved in neoplasms and multiple sclerosis \[[@CR45], [@CR46]\].

The common non-PD-specific genes were associated with 100 experimentally validated miRNAs. Out of these miRNAs, 25 were previously known to be associated with PD (Table [7](#Tab7){ref-type="table"}). The NCOA3, AKAP13 and BCL2 genes were regulated by the highest numbers of miRNAs (21, 19 and 14 miRNAs, respectively). NCOA3, a nuclear receptor coactivator, has been reported to play an important role in some biological processes, such as cell proliferation and apoptosis \[[@CR47]\]. Previous study indicate that high expression of NCOA3 can suppress cells apoptosis induced by histone deacetylase inhibitor in breast cancer \[[@CR47]\]. BCL2, an anti-apoptotic gene, regulate cell death (apoptosis), by either inducing or inhibiting apoptosis. The observed upregulation of BCL2 by GDNF suggests that BCL2 protects neurons, possibly by regulating cell apoptosis. Previous studys have shown its alteration in Alzheimer's Disease \[[@CR35]\]. AKAP13 protein functions as an anchor protein for the regulatory subunit of protein kinase A which is reported to increase tau phosphorylation. AKAP13 are likely more specifically involved in tau phosphorylation pathways that is highly related to Alzheimer's disease \[[@CR48], [@CR49]\]. Out of these three genes, NCOA3 was associated with hsa-let-7 g-5p, hsa-miR-106b-5p, hsa-miR-17-5p, hsa-miR-181c-5p, hsa-miR-20a-5p, hsa-miR-25-3p and hsa-miR-27b-3p that were reported to be previously known to be associated with PD (Table [7](#Tab7){ref-type="table"}). AKAP13 was associated with hsa-miR-34a-5p, and hsa-miR-34c-5p that were reported to be previously known to be associated with PD (Table [7](#Tab7){ref-type="table"}). BCL2 was associated with hsa-miR-132--5p, hsa-miR-182--5p, hsa-miR-23b-3p and hsa-miR-34a-5p that were reported to be previously known to be associated with PD (Table [7](#Tab7){ref-type="table"}). Overall, the 100 miRNAs were identified to be associated with 98 mTFs.Table 7miRNAs that are previously known to be associated with PD associated with non-PD-specific genes identified in our studymiRNAs from bloodmRNA targets associated with the miRNAsmiRNAs from substantia nigramRNA targets associated with the miRNAshsa-let-7 g-5pNCOA3hsa-let-7a-5pSNAP23hsa-let-7i-5pIRS2hsa-let-7e-5pSNAP23hsa-miR-101-5pTMX4hsa-miR-106a-5pNINhsa-miR-106a-5pNCOA3hsa-miR-106b-5pNINhsa-miR-106b-5pABHD5, NCOA3hsa-miR-125a-5pNINhsa-miR-125a-5pMAP2K7hsa-miR-17-5pSNTB2hsa-miR-132-3pATMhsa-miR-182--5pSNTB2hsa-miR-132--5pBCL2hsa-miR-20a-5pNINhsa-miR-143-3pFAM102Ahsa-miR-23b-3pNIN, SNTB2hsa-miR-17-5pABHD5, NCOA3hsa-miR-25-3pSNTB2hsa-miR-181c-5pNCOA3hsa-miR-302d-5pSNTB2hsa-miR-182--5pBCL2hsa-miR-30a-5pSNTB2hsa-miR-18a-5pATMhsa-miR-30e-5pSNTB2hsa-miR-20a-5pABHD5, NCOA3hsa-miR-93-5pSNTB2hsa-miR-212-3pLRRFIP1hsa-miR-9-5pSNTB2hsa-miR-214-3pFAM120Ahsa-miR-23b-3pBCL2hsa-miR-25-3pNCOA3hsa-miR-27a-3pABCA1, FAM102A, FAM120Ahsa-miR-27b-3pNCOA3hsa-miR-30e-5pIRS2hsa-miR-34a-5pAKAP13, ATM, BCL2hsa-miR-34c-5pAKAP13hsa-miR-363-3pIRS2hsa-miR-93-5pABCA1, ABHD5

In brain SN, 333 gTFs were associated with 14 common genes. The NIN, PPFIA2 and ATP5S genes were regulated by the highest numbers of gTFs (139, 102 and 93 gTFs, respectively). NIN, a large coiled-coil protein, is essential for neurogenesis, angiogenesis and stem cell fate \[[@CR50]\]. Previous study showed that alternative splicing in NIN gene appears sufficient for neural stem cells differentiation into neurons \[[@CR50]\]. PPFIA2 is a member of the leukocyte common antigen-related PTP interacting protein family liprin \[[@CR47]\]. It is known to be downregulated by androgens in prostate cancer cell lines \[[@CR47]\] and is found to play crucial functions in the pathophysiology of schizophrenia and bipolar disorder \[[@CR51]\]. ATP5S gene encodes a subunit of mitochondrial ATP synthase which utilizes the electro chemical gradient to synthetize ATP from ADP in inner mitochondrial membrane by oxidative phosphorylation \[[@CR52]\]. Mitochondrial dysfunction is reported to be a pathological pathway associated with PD \[[@CR53]\]. But little is known about the involvement of ATP5S in PD.

The common non-PD-specific genes were associated with 50 experimentally validated miRNAs. Out of these miRNAs, 16 were previously known to be associated with PD (Table [7](#Tab7){ref-type="table"}). It was found that NIN and SNTB2 were regulated by the highest numbers of miRNAs (18 and 25 miRNAs, respectively), and that NIN and SNTB2 were associated with 16 miRNAs (such as miR-9-5p) (Table [7](#Tab7){ref-type="table"}) that were previously associated with PD. For example, in study, miR-9-5p increased by more than three times in treated PD patients compared with those of controls \[[@CR54]\]. Moreover, miR-9-5p were in turn regulated by mTFs of interleukin 1 beta (IL1B) and nuclear factor kappa B subunit 1 (NFKB1). Genetic variation in the proinflammatory cytokine gene IL1B can contribute to risk of developing PD. These finding supported that miR-9-5p may contributes to the pathogenesis of sporadic PD. Previous finding reported SNTB2 involved in atherosclerosis \[[@CR55]\], and that its risk factors in atherosclerosis include higher plasma urate level \[[@CR56]\] which has linked to lower risk of PD in men \[[@CR57]\]. Overall, the 50 miRNAs were identified to be associated with 63 mTFs.

These data support the finding that these genes may be significant factors in PD, but furture studys are needed.

SNP analysis of the common PD-specific and other neurodegenerative disease-specific genes {#Sec20}
-----------------------------------------------------------------------------------------

We next obtained SNPs corresponding to the common DEGs (both the PD-specific and the neurodegenerative disease-specific genes) of our study from the online databases MirSNP, miRdSNP and LincSNP 2.0. Six miRNAs were found to be associated with two PD-specific SNPs corresponding to the PD-specific gene SNCA (Table [8](#Tab8){ref-type="table"}). Two SNPs were associated with these six miRNAs (Table [8](#Tab8){ref-type="table"}). The chromosome locus and allele gene for these two SNPs were then searched using the dbSNP database. These two SNPs were located on chromosome 4.Table 8SNPs in PD with associated PD-specific genes and miRNAsMicroRNAsSNPsAlleleChromosomeGeneMinorAlleleCountExperimentally_confirmedhsa-miR-141rs17016074A/G4:89726127SNCAA = 0.0503/252hsa-miR-153rs17016074A/GSNCAYeshsa-miR-223rs17016074A/GSNCAhsa-miR-499-3prs17016074A/GSNCAhsa-miR-504rs17016074A/GSNCAhsa-miR-7rs17016074A/GSNCAYeshsa-miR-141rs356165A/G4:89725735SNCAA = 0.4842/2425hsa-miR-153rs356165A/GSNCAYeshsa-miR-223rs356165A/GSNCAhsa-miR-504rs356165A/GSNCAhsa-miR-7rs356165A/GSNCAYes

Twenty-three miRNAs were found to be associated with seven neurodegenerative disease-specific SNPs corresponding to four neurodegenerative disease-specific genes (Table [9](#Tab9){ref-type="table"}). Seven SNPs were associated with these 23 miRNAs. The chromosome locus and allele gene for these seven SNPs were then searched using the dbSNP database. The chromosome loci of these SNPs include 6, 14, 16 and 18.Table 9SNPs with their associated miRNAs and genes in neurodegenerative diseases other than PDMicroRNAsSNPsAlleleChromosomeGeneMinorAlleleCountRegionSNPs related diseases by experimentally_confirmedhsa-miR-4428rs1049853C/T6:31269123HLA-CA = 0.0903/4523'UTRAlzheimer's disease, Rheumatoid arthritishsa-miR-27a-5prs1016860A/G18:63127841BCL2*T* = 0.1166/5843'UTRMultiple sclerosishsa-miR-3127-5prs1016860A/G18:63127841BCL2T = 0.1166/5843'UTRhsa-miR-3158-3prs1016860A/G18:63127841BCL2T = 0.1166/5843'UTRhsa-miR-4720-3prs1016860A/G18:63127841BCL2T = 0.1166/5843'UTRhsa-miR-4789-5prs1016860A/G18:63127841BCL2T = 0.1166/5843'UTRhsa-miR-629-5prs1016860A/G18:63127841BCL2T = 0.1166/5843'UTRhsa-miR-1224-5prs1058929A/G/T16:8853394CARHSP1C = 0.4740/23743'UTRAlzheimer's diseasehsa-miR-1265rs1058929A/G/T16:8853394CARHSP1C = 0.4740/23743'UTRhsa-miR-3605-5prs1058929A/G/T16:8853394CARHSP1C = 0.4740/23743'UTRhsa-miR-3915rs1058929A/G/T16:8853394CARHSP1C = 0.4740/23743'UTRhsa-miR-4710rs1058929A/G/T16:8853394CARHSP1C = 0.4740/23743'UTRhsa-miR-5585-5prs1058929A/G/T16:8853394CARHSP1C = 0.4740/23743'UTRhsa-miR-193a-5prs9953G/T16:8853271CARHSP1G = 0.4507/22573'UTRAlzheimer's disease, glioblastomahsa-miR-3190-3prs9953G/T16:8853271CARHSP1G = 0.4507/22573'UTRhsa-miR-3615rs9953G/T16:8853271CARHSP1G = 0.4507/22573'UTRhsa-miR-335-5prs1058967A/G16:8853179CARHSP1C = 0.4744/23763'UTRAlzheimer's diseasehsa-miR-451brs1058967A/G16:8853179CARHSP1C = 0.4744/23763'UTRhsa-miR-4533rs1058967A/G16:8853179CARHSP1C = 0.4744/23763'UTRhsa-miR-4797-5prs1058967A/G16:8853179CARHSP1C = 0.4744/23763'UTRhsa-miR-3679-3prs2447924G/T16:8854616CARHSP1C = 0.4754/23813'UTRAlzheimer's diseasehsa-miR-4286rs2447924G/T16:8854616CARHSP1C = 0.4754/23813'UTRhsa-miR-4517rs8017316C/G14:50322860ATP5SG = 0.4511/2259IntronAlzheimer's disease

Discussion {#Sec21}
==========

Meta-analysis of common genes {#Sec22}
-----------------------------

In the past decades, microarray has been widely used to identify DEGs and pathways underlying PD pathogenesis. In these studies, most microarray datas are mainly from the brain regions or blood in PD. The analysis of this brain regions from PD may only highlight genes associated with changes in cellular composition \[[@CR58]\], However, brain tissues are not easily obtained. In recent years, there is evidence to prove that low concentrations of urate in the blood serum increased the risk of PD \[[@CR12]\], and that there is a growing interest in the discovery of blood biomarkers for PD. Some genes have been proved to take part in important neurodegeneration molecular pathways in PD patients.

A large number of microarray gene expression studies on PD have been performed, but some of the results showed low consistency among involved genes and pathways. Meta-analysis has suggested new biological insights, as well as identification of greater consistent genes and pathways potential in PD pathogenesis \[[@CR13], [@CR15], [@CR16]\]. Numbers of meta-analysis on PD have been conducted, with few focused on understanding the regulatory network involving genes, miRNAs, transcription factors (TFs) and lncRNA in PD progression. Moreover, few meta-analysis on PD is focused on comparative analysis of microarray datases profiles from brain and blood samples.

In our study, meta-analysis merges several datasets from PD brain SN or blood into a single analysis to obtain more meaningful set of DEGs and these DEGs are analyzed in various biogenetic databases to get related non-coding RNAs and gTFs, among others. In this study of DEGs in PD patients compared with healthy controls, we identified 36 common DEGs in blood studies and 17 common DEGs in five brain SN studies using meta-analysis technique-Vote counting generic ways. Of the total common genes, 28 genes were previously reported as associated with other neurodegenerative diseases (see Tables [2](#Tab2){ref-type="table"} and [3](#Tab3){ref-type="table"}), but not previously known to be associated with PD, and five genes were previously shown to be associated with PD. Among the 28 genes associated with other neurodegenerative diseases, 22 were downregulated and six were upregulated in PD compared with controls. We further studied these 28 genes for their association with PD. Out of the five genes previously associated with PD, three genes (HSPA6, MAP2K6 and SRPK2 genes) were downregulated and two genes (NOL7 and SNCA genes) were upregulated in PD compared with controls. A recent study reported the role of HSPA6 in SH-SY5Y neuronal cells which are used as a model system for neurodegenerative diseases such as Alzheimer's disease and PD \[[@CR59]\]. HSPA6, aslo known as Hsp70B', is a member of heat shock proteins (Hsps) 70 family and the Hsp70 family of heat shock proteins has been well known for its roles in cytoprotective effects against cell death and implicated in neuroprotection \[[@CR59]\]. It was also found that HSPA6 has previously been examined in human neurodegenerative diseases and was proposed as a potential treatment strategy to counter PD \[[@CR60]\]. Studies have found that MAP2K6 binds to and regulates the expression of the PD-related protein LRRK2, which is a prevalent cause of sporadic PD \[[@CR33]\]. SRPK2, a serine/arginine protein-specific kinase, is highly expressed in the brain. Increased amounts of SRPK2 can lead to the hyperphosphorylation of serine-arginine-rich proteins, which in turn induces changes in alternative pre-mRNA splicing observed in PD \[[@CR34]\]. NOL7 is a candidate cancer suppressor that localizes to 6p23, a segment with frequent loss of heterozygosity (LOH) in many tumors \[[@CR61]\]. NOL7 interacts with amyloid precursor protein (APP) protein which accumulates in mitochondrial membrane in PD, and that APP interacts with LRRK2 and then is phosphorylated at Thr668 within its intracellular domain to promote neurotoxicity in PD \[[@CR62]\]. SNCA was reported to be associated with PD, and different SNPs in the SNCA gene had a correlation with increased or decreased risk of PD \[[@CR63]\]. All these data help validate our results showing an association of these genes in PD.

Moreover, our study screen the overlapped genes among common DEGs obtained from two different tissues (brain SN and blood), but did not get any common genes. And It was found that the common DEGs from blood and brain SN are enriched in different BP.

The analysis of lncRNAs and mTFs associated with these PD-specific miRNAs related to the PD-specific genes {#Sec23}
----------------------------------------------------------------------------------------------------------

The biological regulations of the common PD-specific genes are shown in Tables [4](#Tab4){ref-type="table"} and [5](#Tab5){ref-type="table"}. We examined the relevance of experimentally validated miRNAs and lncRNAs with these five PD genes. Our analysis showed that four out of the five PD-specific genes (HSPA6, SRPK2, NOL7 and SNCA) are regulated by 14 PD-specific miRNAs (Tables [4](#Tab4){ref-type="table"} and [5](#Tab5){ref-type="table"}). Interestingly, SRPK2 (obtained from blood) was regulated by the highest number of PD-specific miRNAs (hsa-miR-17-3p, hsa-miR-181a-5p, hsa-miR-181b-5p, hsa-miR-181c-5p, hsa-miR-181d, hsa-miR-19a-3p and hsa-miR-19b-3p) and these seven PD-specific miRNAs are regulated by 25 lncRNAs and 16 mTFs (Table [6](#Tab6){ref-type="table"}).

mTF-miRNA-gene-gTF regulatory network analysis {#Sec24}
----------------------------------------------

In this study, we screened the overlapped gTFs of PD-specific and non-PD-specific genes mTF-miRNA-gene-gTF regulatory network obtained from blood studies. Result showed that 135 gTFs are found to be involved in these two networks. It was found that core promoter element-binding protein (CPBP), Spi-B transcription factor (SPIB) and KID3 are the gTFs, which regulated a maximum number of genes, 28, 27 and 24 genes respectively. The CPBP was associated with four PD-specific genes (HSPA6, MAP2K6, SRPK2 and NOL7) and 24 non-PD-specific genes as obtained from TRANSFAC database \[[@CR27]\]. SPIB was found to be associated with four PD-specific genes(HSPA6, MAP2K6, SRPK2 and NOL7) and 23 non-PD-specific genes. KID3 was found to be associated with four PD-specific genes(HSPA6, MAP2K6, SRPK2 and NOL7) and 20 non-PD-specific genes in GeneCards database (<http://www.genecards.org>/). CPBP is reported to be with lung cancer \[[@CR64]\] and SPIB is found to be associated with leukemia cells \[[@CR65]\]. Previous study reported that KID3 was to be highly expressed in adult brain \[[@CR66]\]. These three gTFs are not reported previously to be associated with PD. But they regulate the four PD-specific genes and many non-PD-specific genes from blood.

In PD-specific and non-PD-specific genes mTF-miRNA-gene-gTF regulatory network obtained from brain SN studies, 59 overlapped gTFs is found to be involved in these two networks. CPBP, caudal type homeobox 1(CDX1) and SPIB are the gTFs, which regulated maximum genes i.e. 13, 12 and 12 genes, respectively. All three gTFs were found to regulate the PD-specific gene SNCA and 12, 11, and 11 non-PD-specific genes respectively.

The overlapped gTFs were filtered out in PD-specific regulatory patterns from blood and brain SN and 42 gTFs were screened. CPBP and SPIB are the gTFs, which regulated a maximum number of genes, 5 PD-specific genes(HSPA6, MAP2K6, SRPK2, NOL7 and SNCA) respectively.

Finally, we screened the overlapped gTFs in non-PD-specific regulatory patterns from blood and brain SN, 236 gTFs were identified. CPBP, KID3 and SPIB are the gTFs, which regulated maximum genes i.e. 36, 34 and 34 genes.

Combination of above results, we supposed that CPBP and SPIB may play significant roles in PD. But there is no verification of experiments. So further researches are needed.

Next, analysis of these miRNAs associated with PD-specific genes and non-PD-specific genes resulted in the identification of common PD-specific miRNAs (hsa-miR-181c-5p, hsa-miR-23b-3p and hsa-miR-34a-5p). Hsa-miR-181c-5p was shown to regulate PD-specific gene SRPK2 from blood (shown in Table [4](#Tab4){ref-type="table"}) and regulate non-PD-specific gene NCOA3 from blood (Table [7](#Tab7){ref-type="table"}). Hsa-miR-23b-3p is found to regulate PD-specific gene SNCA from brain SN (shown in Table [4](#Tab4){ref-type="table"}), and regulate non-PD-specific gene BCL2 from blood and NIN and SNTB2 from brain SN (Table [7](#Tab7){ref-type="table"}). It was found that hsa-miR-34a-5p regulate PD-specific gene SNCA from brain SN (shown in Table [4](#Tab4){ref-type="table"}), and regulate non-PD-specific gene AKAP13, ATM and BCL2 from blood and SNTB2 from brain SN (Table [7](#Tab7){ref-type="table"}). Maybe, these three common miRNAs play an important role in PD, but further study is needed.

The biological regulations of the common PD-specific genes are shown in Tables [4](#Tab4){ref-type="table"} and [5](#Tab5){ref-type="table"}. We examined the relevance of experimentally validated miRNAs and lncRNAs with these five PD genes. Our analysis showed that four out of the five PD-specific genes (HSPA6, SRPK2, NOL7 and SNCA) are regulated by 14 PD-specific miRNAs (Tables [4](#Tab4){ref-type="table"} and [5](#Tab5){ref-type="table"}). Interestingly, SRPK2 was regulated by the highest number of PD-specific miRNAs (hsa-miR-17-3p, hsa-miR-181a-5p, hsa-miR-181b-5p, hsa-miR-181c-5p, hsa-miR-181d, hsa-miR-19a-3p and hsa-miR-19b-3p) and these seven PD-specific miRNAs are regulated by 25 lncRNAs and 16 mTFs (Table [6](#Tab6){ref-type="table"}).

Identification of feed forward loops from mTF-miRNA-gene-gTF regulatory networks {#Sec25}
--------------------------------------------------------------------------------

Analysis of the regulatory network of PD-specific genes identified in blood datasets revealed the presence of two interesting feed-forward loops (FFLs), in which a TF regulates a miRNA and in turn they both regulate a target gene. One FFL was detected among the SPRK2 gene, hsa-miR-19a-3p and Spi-1 proto-oncogene (SPI1). The TransmiR information indicated that hsa-miR-19a-3p is inhibited by the SPI1 TF. By combing the TransmiR and DIANA-Tarbase data, we found that SPI1 and hsa-miR-19a-3p both regulate the expression of its target gene SPRK2. Studies have found that SPI1 plays an important role in the regulation of genes related to the function of microglia, which accumulate in PD \[[@CR67]\]. Another study revealed that lower SPI1 expression reduces risk for AD by regulating myeloid gene expression and cell function \[[@CR68]\]. Another study revealed a function of miR-19a in multiple system atrophy, which is a sporadic neurodegenerative disease, identifying lower expression of this miRNA in multiple system atrophy patients than controls \[[@CR69]\]. Moreover, previous study reported the association of downregulated levels of miR-19a in idiopathic Parkinson's disease (IPD) \[[@CR70]\]. This information provides functional insights into the finding of downregulation of hsa-miR-19a-3p by SPI1. However, the FFL gene SPRK2 was identified downregulated in our study.

The other FFL is among the SPRK2 gene, hsa-miR-17-3p and SPI1. The TransmiR information indicated that hsa-miR-17-3p is inhibited by the SPI1 TF. By combing the TransmiR and DIANA-Tarbase data, we found that SPI1 and hsa-miR-17-3p both regulate the expression of its target gene SPRK2. Another study showed that miR-17 is upregulated in multiple system atrophy \[[@CR71]\] and reduced during neurogenesis \[[@CR72]\]. By combining this information, we elucidated possible relationships among the genes, miRNAs and TFs in FFLs in neurodegenerative disease by using Cytoscape software v. 3.4.0 (Fig. [8](#Fig8){ref-type="fig"}). Further study on these FFLs may help us understand the molecular biology of PD.Fig. 8The feed forward loop from mTF-miRNAgene-gTF regulatory networks obtained from Cytoscape software. In this network, the diamond shaped magenta nodes represent miRNAs, the triangle-shaped magenta nodes represent TFs, and the round rectangle-shaped red nodes represent the PD-specific genes

Significance of IncRNAs and mTFs associated with these PD-specific miRNAs related to the PD-specific genes {#Sec26}
----------------------------------------------------------------------------------------------------------

The network showed the relationship among the PD-specific miRNAs, mTFs and lncRNAs. In this regulatory network, we identified 45 lncRNAs associated with these PD-specific miRNAs. The regulations of 31mTFs associated with 14 PD-specific miRNAs showed that most of these interactions had an activating effect (Table [6](#Tab6){ref-type="table"}). The regulation of hsa-miR-204-3p of this network is noteworthy. It is associated with one lncRNAs but has no mTFs associated with it (Table [6](#Tab6){ref-type="table"}). This miRNA represses HSPA6 which is known to be involved in PD \[[@CR59]\]. In contrast to the above findings, hsa-miR-181a-5p has a maximum number of lncRNAs (15 lncRNAs) associated with it (Table [6](#Tab6){ref-type="table"}). This miRNA is in turn repress one gene namely SRPK2. This gene is known to be involved in PD \[[@CR34]\]. However, the miRNA is associated with only one mTFs (TP63) (Table [6](#Tab6){ref-type="table"}). In order to find out the functional role of these 45 lncRNAs, we further analyzed them in EVLncRNAs database ([http://biophy.dzu.edu.cn/EVLncRNAs/](http://biophy.dzu.edu.cn/EVLncRNAs)) \[[@CR73]\]. Two lncRNAs (NEAT1 and TUG1) were found to be associated with neurodegenerative disease. Nuclear paraspeckle assembly transcript 1 (NEAT1) is known to be involved in neurological disorders including amyotrophic lateral sclerosis \[[@CR74]\] and Huntington's disease \[[@CR75]\]. Taurine up-regulated 1 (TUG1) is known to be involved in neurological disorders including Huntington's disease \[[@CR76]\]. But their functions in PD were not identified by previous studies.Since these two lncRNAs of this regulatory network are associated with PD-specific miRNAs of our study, they might play an important role in PD.

Significance of SNP analysis of the common PD-specific and other neurodegenerative disease-specific genes {#Sec27}
---------------------------------------------------------------------------------------------------------

Two SNPs were identified to be associated with the SNCA gene (shown in Table [8](#Tab8){ref-type="table"}), which is in turn controlled by miRNAs including the PD-specific miRNA hsa-miR-153. The dbSNP database showed that SNCA is located on chromosome 4. SNPs rs17016074 and rs356165 are location in the SNCA gene, which is related to PD \[[@CR77], [@CR78]\]. Table [8](#Tab8){ref-type="table"} shows that these two SNPs are identified to be associated with SNCA gene which is in turn controlled by PD-specific miRNA hsa-miR-153 and this finding was experimentally validated \[[@CR77], [@CR78]\]. This strengthens the association of these two SNPs in PD. In addition, SNPs rs17016074 and rs356165 in SNCA gene associated with hsa-miR-7 is involved in PD and it was experimentally validated. However, hsa-miR-7 was not previously linked with PD. Therefore, hsa-miR-7 may be a significant PD epigenetic biomarker that requires further study.

We also analyzed SNPs associated with several neurodegenerative disease-specific genes that are in turn controlled by several miRNAs(shown in Table [9](#Tab9){ref-type="table"}). Twenty-three SNPs were identified from four neurodegenerative disease-specific genes. HLA-C was found to be located on chromosome 6 and several previous studies have found the function of HLA-C in multiple sclerosis that was activated by killer cell immunoglobulin like receptor, two Ig domains and long cytoplasmic tail 2(KIR2DL2) \[[@CR79]\]. The SNP rs1049853 in HLA-C was identified in multiple sclerosis patients and was controlled by the miRNA hsa-miR-4428. BCL2 is located on chromosome 18; SNP rs1016860 is located in the BCL2 gene, and the two possible nucleotide variations (A or G) are considered alleles for this base position. This SNP is regulated by several miRNAs: hsa-miR-27a-5p, hsa-miR-3127-5p, hsa-miR-3158-3p, hsa-miR-4720-3p, hsa-miR-4789-5p and hsa-miR-629-5p. MiR-27a was previously shown to be involved in neurodegenerative diseases, such as Huntington's disease and Alzheimer disease \[[@CR80], [@CR81]\]. MiR-3127 and miR-629 were shown to play a role in cancer \[[@CR82], [@CR83]\]. CARHSP1 is located on chromosome 16 and there are four SNPs (rs1058929, rs9953, rs1058967 and rs2447924) at this specific base position; three possible nucleotide variations (A, G or T) are considered alleles for the rs1058929 base position, two possible nucleotide variations (G or T) are considered alleles for rs9953 and rs2447924 base positions, and two possible nucleotide variations (G or A) are considered alleles for the rs1058967 base position. These four SNPs may be regulated by several miRNAs: hsa-miR-1224-5p, hsa-miR-1265, hsa-miR-3605-5p, hsa-miR-3915, hsa-miR-4710, hsa-miR-5585-5p, hsa-miR-193a-5p, hsa-miR-3190-3p, hsa-miR-3615, hsa-miR-335-5p, hsa-miR-451b, hsa-miR-4533, hsa-miR-4797-5p, hsa-miR-3679-3p and hsa-miR-4286. MiR-1224 was found the fountion in PD targeting LRRK2 \[[@CR84]\]. MiR-1265, miR-193a, miR-3615, miR-335 and miR-4286 were shown to play a role in cancer \[[@CR85]--[@CR88]\]. Two SNPs, rs6430498 and rs12512664, in miR-3679 were significantly associated with femoral neck bone mineral density \[[@CR89]\]. Therefore, future studies should pursue the function of these SNPs and associated miRNAs in PD disease.

Conclusion {#Sec28}
==========

In this study, we performed a meta-analysis with three microarray datasets from PD blood studies and five from PD SN studies to study DEGs, gene regulatory networks and lncRNA-mediated regulatory networks. The meta-analysis identified 36 common DEGs in PD blood studies and 17 common DEGs in PD SN studies. Further, we identified five PD-specific genes in our study: HSPA6, MAP2K6, SRPK2, NOL7 and SNCA genes. Analysis of the regulatory miRNAs associated with the common PD-specific genes resulted in the identification of PD-specific miRNAs (hsa-miR-204-3p, hsa-miR-17-3p, hsa-miR-181a-5p, hsa-miR-181b-5p, hsa-miR-181c-5p, hsa-miR-181d, hsa-miR-19a-3p, hsa-miR-19b-3p, hsa-miR-93-3p, hsa-miR-153, hsa-miR-23b-3p, hsa-miR-34a-5p, hsa-miR-9-5p and hsa-miR-7-5p). Analysis of the mTF-miRNA-gene-gTF network also led to the identification of two FFLs: one FFL between the SPRK2 gene, hsa-miR-19a-3p and SPI1 and the other between the SPRK2 gene, hsa-miR-17-3p and SPI. In the lncRNA-mediated regulatory network, 45 lncRNAs were associated with known PD-specific miRNAs in our study. Notably, these were not identified in previous studies and warrant further investigation as they might be important epigenetic regulators in PD. Moreover, SNP analysis identified two significant SNPs associated with PD-specific genes and six regulatory miRNAs, and seven significant SNPs associated with other neurodegenerative disease-specific genes and 23 regulatory miRNAs. These SNPs could be considered as latent risk factors for further validation. Thus, the findings of our study should be explored in further investigation of PD.

Additional file
===============

 {#Sec29}

Additional file 1:**Table S1**. GO biological processes analysis of diferentially expressed genes in blood of Parkinson's Disease. **Table S2**. GO biological processes analysis of diferentially expressed genes in substantia nigra of Parkinson's Disease. **Table S3**. The PPI information for common DEGs in blood obtained from STRING online software. **Table S4**. The PPI information for common DEGs in substantia nigra obtained from STRING online software. **Table S5**. The mTF-miRNA-gene-gTF regulatory network of common PD-specific genes identified in blood datasets. **Table S6**. The mTF-miRNA-gene-gTF regulatory network of common PD-specific genes identified in brain substantia nigra. **Table S7**. The mTF-miRNA-gene-gTF regulatory network of common genes identified in blood datasets(other than PD-specific genes). **Table S8**. The mTF-miRNA-gene-gTF regulatory network of common genes identified in substantia nigra datasets(other than PD-specific genes). (XLS 1121 kb)

BP

:   Biological processes

DE

:   Differentially expressed

FFL

:   Feed forward loop

gTF

:   Gene transcription factor

HC

:   Healthy control

lncRNA

:   Long non-coding RNA

miRNA

:   microRNA

mTF

:   microRNA transcription factor

PD

:   Parkinson's disease

PPI

:   Protein protein interaction

SN

:   Substantia nigra

SNP

:   Single nucleotide polymorphism

TF

:   Transcription factor

**Electronic supplementary material**

The online version of this article (10.1186/s12920-018-0357-7) contains supplementary material, which is available to authorized users.

We thank Edanz Group ([www.edanzediting.com/ac](http://www.edanzediting.com/ac)) for editing a draft of this manuscript.

Funding {#FPar1}
=======

This work was supported by the Key Project of Hebei North University (No. 120177), science and technology bureau research development plan of zhangjiakou city in Hebei (No. 0911021D-4) and the Science and Technology Research Project of Hebei Province Department Institutions of Higher Learning (No. Z2015047), China.

Availability of data and materials {#FPar2}
==================================

All the data supporting the findings of this study are contained in the manuscript and the corresponding supplementary files.

LS and HW designed the study, LS, CW, CZ, HW and XS performed the data analysis, and LS wrote the manuscript. HW supervised this work. All authors read and approved the final manuscript.

Ethics approval and consent to participate {#FPar3}
==========================================

Not applicable.

Consent for publication {#FPar4}
=======================

Not applicable.

Competing interests {#FPar5}
===================

The authors declare that they have no competing interests.

Publisher's Note {#FPar6}
================

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
