Characterizing Relevant MicroRNA Editing Sites in Parkinson’s Disease

MicroRNAs (miRNAs) are extensively edited in human brains. However, the functional relevance of the miRNA editome is largely unknown in Parkinson’s disease (PD). By analyzing small RNA sequencing profiles of brain tissues of 43 PD patients and 88 normal controls, we found that the editing levels of five A-to-I and two C-to-U editing sites are significantly correlated with the ages of normal controls, which is disrupted in PD patients. We totally identified 362 miRNA editing sites with significantly different editing levels in prefrontal cortices of PD patients (PD-PC) compared to results of normal controls. We experimentally validated that A-to-I edited miR-497-5p, with significantly higher expression levels in PD-PC compared to normal controls, directly represses OPA1 and VAPB. Furthermore, overexpression of A-to-I edited miR-497-5p downregulates OPA1 and VAPB in two cell lines, and inhibits proliferation of glioma cells. These results suggest that the hyperediting of miR-497-5p in PD contributes to enhanced progressive neurodegeneration of PD patients. Our results provide new insights into the mechanistic understanding, novel diagnostics, and therapeutic clues of PD.

To comprehensively characterize miRNA editing sites in brain tissues of PD patients, we analyzed 43 and 88 small RNA sequencing profiles of PD and control brain tissues, respectively. We identified 362 miRNA editing sites that had significantly different editing levels in the prefrontal cortices of PD patients (PD-PC) compared with the control group. The editing levels of five A-to-I and two C-to-U editing sites were significantly correlated with the ages of normal controls, which were disrupted in PD patients. One A-to-I editing site in miR-497-5p had significantly higher editing levels compared to normal controls in PD-PC samples and edited miR-497-5p directly repressed OPA1 and VAPB. We also demonstrated that overexpression of A-to-I edited miR-497-5p downregulated the expression of OPA1 and VAPB, and repressed proliferation of glioma cells. These results demonstrate that miRNA editing is severely disturbed and relevant in PD and offer novel insights into the etiology of PD. The abnormal miRNA editing patterns in PD might be used to develop novel diagnostic biomarkers and/or therapeutic targets for PD.

The Small RNA Sequencing Profiles Used
To comprehensively identify miRNA mutation and editing (M/E) sites in PD, we collected 131 sRNA-seq profiles of postmortem PD patients and normal people from the NCBI SRA Database (Supplementary Table S1). These profiles included 29 prefrontal cortex samples of PD patients (PD-PC), 14 amygdalae samples of PD patients (PD-Am), 36 prefrontal cortex samples of normal controls (PC), 14 amygdalae samples of normal controls (Am), 6 frontal cortex samples of normal controls (FC), 6 corpus callosum samples of normal controls (CC), 2 inferior parietal lobe samples of normal controls (IPL), 2 temporal neocortex gray matter samples of normal controls (NG), 3 astrocyte cell lines of normal controls (As), and 19 unknown brain regions of normal controls (Unknown). The PD-PC and PC samples were from Brodmann Area 9 (BA9).

The Gene and Protein Expression Profiles Used
To understand the potential function of edited miRNAs, we examined the deregulated genes from brain samples of PD patients and normal controls. As listed in Supplementary Table S1, we selected 107 and 122 gene expression profiles of PD patients and normal controls, respectively, from seven cohorts, including two cohorts for the prefrontal cortex (Brodmann Area 9, briefly BA9) (GSE20168 and GSE68719), two cohorts for the substantia nigra (GSE7621 and GSE20292), one cohort for the putamen (GSE20291), one cohort for the cingulate gyrus (GSE110716), and one cohort for the subventricular zone (GSE130752).
A previous study reported proteomics profiles of prefrontal cortex (BA9) samples in 12 PD patients and 12 healthy people [33]. After comparing the proteomics profiles of the 12 PD patients to those of the 12 healthy people, 283 deregulated proteins were detected (limma package, corrected p < 0.05) and used in our analysis.

Genome and Annotation of miRNAs Used
The unmasked human genomic sequence (hg38, GRCh38) was downloaded from the UCSC Genome Browser [34]. The bowtie-build program in the Bowtie package [35] was used to generate index files of the human genome. The pre-miRNA sequences and genomic positions of miRNAs in gff3 format were downloaded from miRBase (release 21) [36].

Analysis of Small RNA Sequencing Profiles
The selected sRNA-seq profiles were analyzed using the MiRME pipeline [37] with the default settings. Briefly, the raw reads whose sequencing scores of the first 25 nucleotides from the 5 end had sequencing scores of 30 or higher were kept as qualified reads. Then, the unique sequences of the remaining reads were obtained and counts of unique reads with more than 18 nucleotides were calculated. Next, the unique reads were aligned to human pre-miRNAs using NCBI BLASTN [38] with the options of "-S 1 -m 8 -e 0.01" and the reads mapped to human pre-miRNAs were retrieved. Next, these reads mapped to pre-miRNAs were aligned to the genome using Bowtie (v1.0.0) [35] with the options of "-a -best -S -v 1". Then, the alignments of reads to the genome were examined using the cross-mapping correction method [39] to adjust the weights or percentages of a unique read at each of its genomic loci. In the main step, the MiRME algorithm with the default parameters was used to identify M/E sites in miRNAs from the sequences and structures of pre-miRNAs, the alignment of reads to the genome generated by Bowtie, the reads mapped to pre-miRNAs, the alignments of reads to pre-miRNAs generated by BLASTN, and the results of the cross-mapping correction method [39].
The following criteria were used to define M/E sites with significant editing levels: (i) the relative level of editing was at least 5%; (ii) at least 10 reads supported the editing event; (iii) the score threshold of sequencing reads was 30; (iv) a multiple-test corrected P-value (using the Benjamini and Hochberg method [40]) smaller than 0.05; and (v) only M/E sites that had significant editing levels in 10% of all samples were kept for further analysis. Then, the obtained results of different samples were combined by a separate program in the MiRME package (see details in [37]). Based on the positions of M/E sites in miRNAs and mutations in dbSNP, the identified M/E sites were classified into nine different editing types, i.e., A-to-I, C-to-U, 3 -A, 3 -U, 3 -Other, 5 -editing, Other, SNP, and Pseudo [37].
To remove M/E sites due to random sequencing errors, 860 M/E sites that had significant editing levels in at least 10% of the selected samples (13 samples) were kept for further analysis. The 860 identified M/E sites were compared to reported editing sites in miRNAs in the DARNED database [41], the RADAR database [42], and in the literature [19,20,37,[43][44][45]. Finally, the predicted M/E sites that belonged to A-to-I, C-to-U, and Other were manually examined.
All identified M/E sites were named by the names of the pre-miRNAs, positions of the M/E sites in the pre-miRNAs, the nucleotides from the reference pre-miRNA sequences in upper cases, and the edited/mutated nucleotides at the sites in lower cases. Additionally, edited miRNA was named by the pre-miRNA name, the position of the M/E site in pre-miRNA, and the edited/mutated nucleotide in lower case.

Comparing the M/E Sites to Reported SNPs
The identified M/E sites were compared to known SNPs in miRNAs organized in [46] (which was based on the dbSNP v137) and reported SNPs in dbSNP (v151). Only sites that satisfied the following criteria were regarded as SNPs: (i) they had the same genomic positions as the SNPs, (ii) had the same nucleotides as the alleles of the SNPs for both the original and changed nucleotides, and (iii) had editing levels of 100% in at least one of the 131 samples selected.

Identifying Conserved Editing Sites in miRNAs
The A-to-I and C-to-U editing sites were compared to their counterparts in Macaca mulatta [47] and Mus musculus [18,19]. The editing sites of the same editing types that were located on the same positions of mature miRNAs of at least two different species were considered as conserved editing sites.

Identifying Age-Related miRNA Editing Sites
The Pearson correlation, as well as its p-value, between editing level (in %) of each of the 860 editing sites in Supplementary Table S2 and the age of death (y) was calculated with the corr function in MatLab (Mathworks, MA, USA) for the 29 PD-PC and 36 PC samples, respectively. The obtained p-values were corrected using the Benjamini and Horchberg method [40] using the mafdr function in MatLab (Mathworks, MA, USA). The editing sites with corrected p-values smaller than 0.05 were regarded as age-related.

Identifying M/E Sites with Significantly Different Editing Levels in PD
The editing levels of 860 M/E sites in the PD-PC and PC, and PC-Am and AM samples were compared, respectively, using the Mann-Whitney U-tests. The obtained p-values were corrected with the Benjamini and Hochberg method [40]. M/E sites with multiple test corrected p-values smaller than 0.05 were regarded as having significantly different editing levels in the PD-PC or PC-Am samples.

Clustering and Principle Component Analysis Using the Editing Levels of Selected M/E Sites
The editing levels of the 276 editing sites that had significantly different editing levels in PD-PC samples when compared to PC samples were used to perform hierarchical clustering with the hclust function in R. The editing levels were multiplied times 100 and added plus one, were then log2-scaled and used to calculate the correlation coefficient between samples, and one minus correlation coefficient values were used as distances between samples. The Ward D2 method was used in clustering analysis. The editing levels times 100 plus one of these 276 editing sites were log2-scaled and used to perform principle component analysis with the prcomp function in R.

Identifying Targets for Original and Edited miRNAs
M/E sites that satisfied the following criteria were chosen to identify the targets of original and edited miRNAs: (i) they had significantly different editing levels in PD-PC samples compared to PC samples; and (ii) were located in the seed regions of mature miRNAs. The targets of original and edited miRNAs were predicted using the MiCPAR algorithm [48] with its default parameters.
As listed in Supplementary Table S1.3, seven PAR-CLIP sequencing profiles prepared from HEK293 cells stably expressing FLAG/HA-tagged AGO proteins (AGO1, AGO2, AGO3, and AGO4) [49] were downloaded from the NCBI SRA database using the series accession number SRP002487. Four PAR-CLIP sequencing profiles prepared from HEK293 cell lines stably expressing HIS/FLAG/HA-tagged AGO1 or AGO2 [50] were downloaded from the NCBI SRA database using the series accession number SRP018015. Raw reads in these profiles were filtered to make sure that the first 25 nucleotides of the qualified reads had sequencing scores of 30 or higher. The 3 adapters were cut for qualified reads. The remaining reads in these 11 profiles were combined and used in the identification of miRNA targets with the MiCPAR algorithm [48]. The annotation of NCBI RefSeq genes in the GTF file, the mRNA sequences of NCBI RefSeq genes (refMrna.fa.gz, version hg38), and soft-masked human genome sequences (version GRCh38) were downloaded from the UCSC Genome Browser [51] and used as inputs of the MiCPAR algorithm. The targets with at least one PAR-CLIP read with T-to-C variation were kept for further analysis.

GO and Pathway Analysis for the Original and Edited miRNAs
The GO (Gene Ontology) term and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment of the genes that were only targeted by the original or edited miRNAs were analyzed with KOBAS2, respectively [52]. The significantly enriched GO terms (with a multiple test corrected p-values smaller than 0.05) were divided into three major categories, i.e., Biological Process, Cellular Component, and Molecular Function. Then, the enriched GO terms and KEGG pathways of the original and edited miRNAs were compared.

Identifying Meaningful Targets for Edited miRNAs
A previous study [33] identified 1095 deregulated genes using the Wald tests (in DESeq2 [53]) by analyzing RNA-seq data sets for 19 PD and 24 control samples of the prefrontal cortex (Brodmann Area 9, BA9). This study also examined the protein abundance levels of BA9 samples of 12 PD patients and 12 normal people using MS3 proteomics and identified 283 proteins with significantly different levels (limma package, corrected p < 0.05) in PD [33]. Another study [54] generated 30 gene expression profiles of BA9 with microarray, with 15 profiles for both PD and PD-PC. The deregulated genes were identified with the limma package.
Because these RNA-seq and proteomics profiles were also from the BA9 regions, which were the same as the regions of PD-PC and PC sRNA-seq profiles in this study, these deregulated genes (mRNAs) and proteins were used to identify meaningful targets for edited miRNAs. For edited miRNAs that had higher editing levels and normalized abundances in the PD-PC samples, we chose their targets, identified with the MiCPAR algorithm, which showed lower expression levels in PD-PC samples, and vise versa.
A human 293T cell line was cultured in DMEM high glucose medium with 10% FBS and 1% NEAA. A total of 500 µL medium was inoculated to each well of a 24-well plate. Totals of 2 µg p497_25g and 2 µg pGLO-OPA1 were added to one well as a group of co-transfection. Similarly, one pGLO plasmid and one pCDNA were co-transfected to 293T cells, respectively. Then, 2 µL of DNA-INVI DNA Transfection Reagent (Invigentech, Carlsbad, CA 92008, USA) was added and mixed into each well. In the experiments, each group had three biological replicates. The TransDetect double-luciferase Reporter Assay Kit (Beijing TransGen Biotech, Beijing, China) was used to detect luciferase activities for each of the biological replicates with three technical repeats. The mean value of the three technical repeats was used as the value of one biological replicate. The luciferase activities of different groups were then compared with two-tailed t-tests.

Conservation Analysis of Edited miR-497-5p Complementary Sites
The conservation scores (PhyloP scores [55] for 30 mammals) of the A-to-I edited miR-497-5p complementary sites on OPA1 and VAPB were downloaded from the UCSC Genome Browser [34]. The sequences of the A-to-I edited miR-497-5p complementary sites on OPA1 and VAPB in 30 mammals were downloaded from the UCSC Genome Browser and used to build phylogenetic trees with ClustalX (v2.1) [56]. The obtained phylogenetic trees were visualized with TreeView (v1.6.6) [57].

Cell Culture
Human glioma cell lines (U-118-MG and Hs683) were purchased from the Cell Bank of the Kunming Institute of Zoology, Chinese Academy of Sciences (Kunming, China). All cells were cultured in Dulbecco modified Eagle's medium (DMEM) (Corning), supplied with 10% FBS and 1% penicillin/streptomycin solution. The cell culture was maintained in a 37 • C incubator with 5% CO 2 .

Cell Proliferation Assay
To assess the cell proliferation activities, cell groups were gathered, resuspended in 10% FBS medium (2 × 10 4 cells/mL), seeded in a 96-well plate and then incubated in a 37 • C incubator with 5% CO 2 . Then, cells were transiently transfected with the p497 and p497_25g plasmids, using Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's protocol. Then, after 24, 48, and 72 h of culture, each well received 10 µL of the Cell Counting Kit-8 solution (Beijing TransGen Biotech, Beijing, China), and was subsequently incubated at 37 • C for 1 h. The proliferation activities of cells were measured using an absorbance microplate reader (PHERAstar FS, BMG LABTECH), at wavelengths of 490 and 450 nm. The cell proliferation activities of different groups were then compared with two-tailed t-tests.

Rna Isolation and Quantitative Real-Time Polymerase Chain Reaction
Total RNAs were extracted from the cells 48 h after p497_25g plasmid treatment using the TriQuick total RNA extraction Reagent Kit (Solarbio, Beijing, China). The first strand complementary DNAs (cDNAs) were generated with the Servicebio RT First Strand cDNA Synthesis Kit (Servicebio, Wuhan, China) for OPA1 and VAPB, according to the manufacture's protocols. Real-time PCR evaluation of mRNAs was performed using the 2xSYBR Green qPCR Master mix PCR Kit (Servicebio, Wuhan, China). The qRT-PCR was performed on a Step one plus Real-time PCR System (Applied Biosystems, ABI). The PCR primer sequences are provided in Supplementary Table S6. The relative RNA expression was examined using the 2 −∆∆Ct method [58]. The expression levels of OPA1 and VAPB in different groups were then compared using two-tailed t-tests.

Role of the Funding Source
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

An Overview of Identified Editing Sites in miRNAs
We used the MiRME pipeline [37] to analyze the 131 sRNA-seq profiles selected with the default settings. Totally, we found 860 significant M/E sites (as shown in Figure 1 and Supplementary Table S2) by employing the criteria of at least 10 reads supporting the M/E sites, multiple test-corrected p-values smaller than 0.05, and the M/E sites being identified in at least 10% of the selected sRNA-seq profiles. The largest category of these M/E sites was 3 -A, accounting for 43.8% (377 sites), followed by 3 -U (41.7%), 3 -Other (5.5%), and 5 -editing (2.7%) (Figure 1a and Supplementary Table S7). A-to-I and C-to-U editing sites accounted for 2.7% (23 sites) and 0.5% (4 sites), respectively (Figure 1a). We found 14 SNPs after comparing the identified M/E sites to SNPs in dbSNP (Figure 1a). The Other type of editing sites accounted for 0.5% (four sites) (Figure 1a).
The A-to-I, C-to-U, and Other sites were further classified based on their types of variations (Supplementary Figure S1a and Supplementary Table S12). Similar to previous work [37], there could be several 3 -editing sites in one pre-miRNA, but most pre-miRNAs only had 1 or 2 central sites and 5 sites in our selected samples (Supplementary Figure S1b). However, some miRNAs may have had a few editing events at the 5' end in our previous study [37].  Figure S1 and Table S7 for source data.

Age-Related Editing Sites in PD and Normal Controls
Because previous studies reported that A-to-I editing level was increasing during the developmental process [21,32], we examined the Spearman correlation (ρ) between editing levels of the 860 editing sites and the ages (at death) of 36 PC and 29 PD-PC samples, respectively. There were 227 sites with significant ρ in PC, but only 60 sites with significant ρ in PD-PC samples (Figure 1b and Supplementary Table S2). When these sites were separated as ones with positive and negative ρ values, the numbers of sites with significant ρ values in PD-PC samples were also much smaller than those in PC samples (Figure 1b and Supplementary Figure S1c). The numbers of types of editing sites were decreased in PD-PC samples too. For example, there were no A-to-I sites with positive ρ values in the PD-PC samples (Supplementary Figure S1c). Furthermore, the Kullback-Leibler divergences of the distributions of different types of editing sites in PC and PD-PC samples were 1.81, 2.48, and 3.5 for all, sites with positive ρ values, and sites with negative ρ values, respectively (Supplementary Figure S1c). These results suggest that the gradually increasing or decreasing editing levels of miRNAs during the aging procedure were severely disturbed in PD.
The ρ values of 23 A-to-I and 4 C-to-U editing sites are carefully examined in Figure 1c and Supplementary Table S7. In the PC samples, there were three and one A-to-I sites with significant positive and negative ρ values, respectively (Figure 1c-e). In comparison, only one A-to-I site had a significant negative ρ value in PD-PC (Figure 1c,d). The median value of ρ of A-to-I sites in PC was slightly larger than 0 (Figure 1c), which is consistent with the ever-increasing A-to-I editing level in the developmental process noticed previously [21,32]. However, the median ρ value of A-to-I sites in PD-PC samples was smaller than 0 (Figure 1c), further suggesting the disrupted A-to-I editing of miRNAs in PD. By contrast, the C-to-U editing sites had a weak negative median ρ value in PC (Figure 1c), suggesting that C-to-U editing sites of miRNAs generally had decreasing editing levels when people were aging. Two C-to-U editing sites (see Figure 1c,d,f) had significant ρ values in PC.

A-To-I Editing Sites
We found 23 significant A-to-I editing sites in total from the selected samples (Supplementary Figure S2a and Table S13). At least 14 of these 23 sites were conserved in other mammals [19,47] (Supplementary Figure S2a). Similar to previous results [19,20,37,47], the 5 and 3 sides of these A-to-I editing sites preferred to be U and G, respectively (Supplementary Figure S2b). For example, two identified A-to-I editing sites are shown in Supplementary Figure S2c

C-To-U Editing Sites
We carefully explored the four C-to-U editing sites (Supplementary Figure S3a and Supplementary Table S14). Three of these four C-to-U editing sites were conserved in primates after being compared with results of monkey [47] (see Supplementary Figure S3a). The neighboring nucleotides of these C-to-U sites preferred to be C on both the 5 and 3 sides (Supplementary Figure S3b), consistent with the CCC motif of APOBEC3G [59]. Three examples of C-to-U editing sites are shown in Supplementary Figure S3c

Relevant miRNA Editing Sites in PD
To identify editing sites that were relevant in PD, we compared the editing levels of 860 significant editing sites in PD-PC and PC, and PD-Am and Am samples. We obtained 276 M/E sites that had significantly different editing levels in PD-PC samples compared to PC samples (Supplementary Figure S5a and Table S3). Most of these 276 sites (63.04% or 174 sites) had increased editing levels in PD-PC samples, and 36.96% of these sites had decreased editing levels in PD-PC samples (Supplementary Figure S5b). Three A-to-I editing sites (hsa-mir-497_25_A_g, hsa-mir-1251_10_A_g, and hsa-mir-1301_52_A_g) had increased editing levels in PD-PC samples, and two A-to-I editing sites (hsa-let-7a-2_28_A_g and hsa-mir-411_20_A_g) showed decreased editing levels in PD-PC samples compared to PC samples (Figure 2a and Supplementary Table S8). Three C-to-U and one Other editing sites had decreased and increased editing levels in PD-PC samples compared to PC samples, respectively (Figure 2b). There were no M/E sites with significantly different editing levels in the comparison between PD-Am and Am samples.
We used the editing levels of these 276 M/E sites to perform principle component analysis and hierarchical clustering for the PD-PC and PC samples. The results show that although PD-PC samples were grouped together, 12 PC samples were clustered together with the PD-PC samples (Figure 2c,d). As indicated by the arrow in Figure 2c, some PC samples seem to be gradually falling to the cluster of PD-PC samples. These results suggest that although these 12 PC samples were clinically normal, their miRNAs showed similar editing patterns to those of PD patients. Therefore, we treated these 12 PC samples as PD-PC samples (named as PD-PC-like samples). Seven of the nine editing sites shown in Supplementary Figure S6a,b and Table S17 demonstrated significantly different editing levels (corrected p < 0.05, Mann-Whitney U-tests) when the 12 PD-PC-like samples were compared to other PC samples, suggesting that these 12 PD-PC-like samples had different miRNA editing patterns compared to other PC samples (see Supplementary Figure S6a,b).
To check the age difference of PD-PC-like samples, we carefully examined the editing levels of hsa-mir-497_25_A_g (see Figure 2e) and ages of PC, PD-PC-like, and PD-PC samples. As in Figure 2e, PD-PC-like samples had higher editing levels for hsa-mir-497_25_A_g (corrected p = 1.9 × 10 −4 , Mann-Whitney U-test, see Supplementary Table S17) and larger ages than other PC samples (p = 1.1 × 10 −3 , Student's t-test, Figure 2f). In comparison, the ages of PD-PC-like and PD-PC samples had no significant difference (p = 0.69, Student's t-test, Figure 2f).  Table S8 for source data.

Target Analysis of A-to-I Editing Sites in Seeds
We selected four A-to-I sites that located in the seed regions of mature miRNAs (Supplementary Table S3) and identified the targets of A-to-I edited and original miRNAs with the MiCPAR pipeline [48], as listed in Supplementary Table S4. In the MiCPAR analysis, 11 PAR-CLIP sequencing profiles were used and only targets with at least 1 PAR-CLIP read with T-to-C variation were kept for further analysis (see Materials and Methods for details). After comparing the targets of A-to-I edited and original miRNAs, we found that these A-to-I editing events severely changed the target sets of these miRNAs (Supplementary Figure S5c), as well as the enriched GO terms and KEGG pathways of these targets (Supplementary Figure S5c and Table S5).
Since A-to-I edited hsa-miR-497-5p (hsa-mir-497_25g) had a much higher expression level than the other three A-to-I edited miRNAs (Supplementary Figure S5d and Table S16), we next carefully examined the targets of hsa-mir-497_25g. After comparing the predicted targets of edited miRNAs with the deregulated genes and proteins in prefrontal cortex (BA9) samples (Figure 3a and Supplementary Table S9, see details in Materials and Methods), we found that hsa-mir-497_25g targeted OPA1 mitochondrial dynamin-like GTPase (OPA1) and VAMP-associated protein B and C (VAPB) (Figure 3b). Presumably due to the upregulation of hsa-mir-497_25g (Figure 3c) in PD-PC samples, OPA1 and VAPB were significantly downregulated in PD-PC samples (Figure 3d,e). The protein expression levels of OPA1 and VAPB were also significantly downregulated in PD-PC samples (corrected p < 0.05, as in [33]).  Table S9 for source data.
Furthermore, OPA1 was significantly downregulated in the substantia nigra (SN) (Supplementary Figure S7a,b and Table S18), and was mildly downregulated in the putamen of PD patients (Supplementary Figure S7c). However, in the subventricular zone (SZ) and cingulate gyrus, the expression of OPA1 did not change significantly (Supplementary Figure S7d,e). In SN, VAPB had a weak decreasing trend in PD patients (Supplementary Figure S8a and Table S19, but had no change in another cohort (Supplementary Figure S8b) and in the putamen, SZ and cingulate gyrus (Supplementary Figure S8c to S8e, respectively).
Importantly, consistent with the significant positive correlation between the editing level of hsa-mir-497_25_A_g and ages of individuals in PC (Figure 1d Figures S7h and S8h, respectively), and were positive in SZ (ρ = 0.29 and 0.60 in Supplementary Figures S7i and S8i, respectively). In PD, these correlations were generally weakened ( Supplementary Figures S7j-l,n and S8j-l, S8n for OPA1 and VAPB, respectively), except in putamen, where the negative correlations were slightly enhanced ( Supplementary Figures S7m and S8m for OPA1 and VAPB, respectively).

Hsa-Mir-497_25g Directly Represses OPA1 and VAPB
The complementary sites of hsa-mir-497_25g on OPA1 and VAPB were only partially conserved in mammals (Figure 4a-d and Table S10). For both OPA1 and VAPB, in the 8mer opposite to the seed of edited miR-497-5p (red rectangles in Figure 4a,b), a few nucleotides had negative PhyloP conservation scores [55], indicating that these two sites were fastevolving in primates. Even within primates, four and seven species did not have conserved sequences complementary to the seed of hsa-mir-497_25g for OPA1 and VAPB, respectively (Figure 4c,d, respectively).
The complementary sites of hsa-mir-497_25g on OPA1 and VAPB were only conserved in some primates including monkey, but non-conserved in rodents (Figure 4a-d). As expected, A-to-I edited miR-497-5p could repress OPA1 and VAPB of monkey (Supplementary Figure S9a,b and Table S20), but could not repress Opa1 and Vapb of mouse (Figure 4g-h).
To summarize, hsa-mir-497_25g directly suppressed OPA1 and VAPB in human and in some non-human primates, but not in mouse and those primates with non-conserved complementary sites of hsa-mir-497_25g. The non-conserved complementary sites of hsamir-497_25g might contribute to the advanced functions and fast evolving of prefrontal cortex of human.

hsa-mir-497_25g Suppresses the Proliferation of Glioma Cells
To further validate that hsa-mir-497_25g repressed OPA1 and VAPB, we transfected p497_25g into two glioma cell lines, i.e., U-118-MG and Hs683, and examined the expression levels of OPA1 and VAPB with quantitative RT-PCR (qRT-PCR) experiments. The results of qRT-PCR experiments showed that transfection of p497_25g plasmid downregulated the expression levels of OPA1 and VAPB in the two cell lines, respectively (Figure 5a,d and Supplementary Table S11).
To explore the role of hsa-mir-497_25g in glioma cell line growth, we transfected U-118-MG and Hs683 with p497 and p497_25g plasmids, respectively. Then, the cell lines were subjected to CCK-8 assay of proliferation. We measured the proliferation activities of the cells after transfection of the plasmids at 24, 48, and 72 h. U-118-MG cells without transfection of plasmids and treated with Lipofectamine 2000 were also measured as controls. The results showed that miR-497_25g had a noticeable inhibitory effect on the proliferation of U-118-MG cells when compared to untreated cells, cells treated with Lipofectamine 2000, and cells transfected with p497 (Figure 5e-g). As shown in Figure 5h,i, miR-497_25g also had a noticeable inhibitory effect on the proliferation of Hs683 cells when compared to untreated cells and cells transfected with p497. Snapshots of the cells as controls and cells transfected with p497_25g are also shown in Figure 5k,l, respectively.
In summary, these results indicate that transfection of the miR-497_25g plasmid significantly represses the proliferation of U-118-MG cells.  Supplementary Table S9. The three species with blue names, i.e., human (hsa, Homo sapiens), monkey (mfa, Macaca fascicularis), and mouse (mus, Mus musculus) were selected for luciferase experiments. (e-h) The luciferase activities when co-transfecting a pGLO plasmid of human and mouse 3 UTRs in Part (a,b) and a pCDNA plasmid containing original pre-hsa-mir-497 (p497) or prehsa-mir-497_25g (p497_25g), respectively. The values shown were mean values ± SDs. p-values were based on two-tailed t-tests. See also Supplementary Figure S9 for monkey results and Supplementary Table S10 for source data.

Discussion
By analyzing small RNA sequencing profiles of brain tissues of 43 PD patients and 88 normal controls, we identified 362 miRNA editing sites with significantly different editing levels in the prefrontal cortices of PD patients (PD-PC) compared to results of normal controls. Although a few studies reported A-to-I editing sites of miRNAs in brains of normal people, our results provide a comprehensive view of the miRNA editome with different types of editing events in brain tissues of PD for the first time, which significantly increases our knowledge about miRNA editing in PD.
The ever increasing editing levels of A-to-I editing during developmental process has been noticed in several previous studies [21,32]. Similarly, we found that the editing levels of three A-to-I editing sites in miRNAs were significantly positively correlated with ages in normal PC samples and these significant correlations were disturbed in PD-PC samples (Figure 1d). Furthermore, our results demonstrate that the editing levels of two C-to-U editing sites (see Figure 1c,f) were negatively correlated with the ages of normal controls, which also disappeared in PD patients. These results suggest potential functional relevances of these C-to-U editing sites in PD.
OPA1 localizes to the inner mitochondrial membrane and plays roles in regulating mitochondrial stability and energy output. Mutations in OPA1 have been associated with optic atrophy type 1 [60,61]. Recently, the association between OPA1 and parkinsonism was also noticed [62,63]. One recent study found that three members of a family had autosomal dominant optic atrophy caused by OPA1 mutations and two of them developed nonsyndromic PD [64]. OPA1 was downregulated in PD-PC (BA9) (Figure 3d,e) and in the substantia nigra of PD patients (Supplementary Figure S7a,b). PD patients carrying the G2019S mutation showed decreased levels of mature OPA1 [65]. OPA1 haploinsufficiency led to progressive loss of iPSC-derived dopaminergic neurons [66]. Autosomal dominant LRRK2 mutations were associated with both familial and sporadic PD [65,67]. LRRK2 directly interacted with OPA1, presumably to regulate membrane dynamics important for mediating endocytosis and mitochondrial function [65].
Neuronal synaptic dysfunction is one of the key features of Parkinson's disease [68,69]. VAPB protein is a membrane protein found in plasma and intracellular vesicle membranes, and is involved in vesicle trafficking. VAPB was downregulated in BA9 of PD-PC (Figure 3d,e). Downregulation of VAPB and its interacting protein PTPIP51 reduced dentritic spine numbers and synaptic activity [70]. A missense mutation (P56S) in VAPB was reported to be associated with motorneuron degeneration in affected amyotrophic lateral sclerosis (ALS) patients [71]. Transgenic mice that expressed human P56S VAPB developed progressive hyperactivities and other motor abnormalities, and showed progressive loss of corticospinal motor neurons [72]. By contrast, overexpression of VAPB in mouse ALS model slowed the motor impairment and promoted the survival of spinal motor neurons [73].
α-synuclein is a PD-related protein that localizes to the nerve terminal [74] and directly binds to VAPB [74] and VAMP2 [75]. Overexpression wild-type and familial Parkinson's disease mutant α-synuclein disrupt the interaction between VAPB and PTPIP51 to loosen ER-mitochondria associations [74]. The downregulation of VAPB in PD patients may exaggerate the toxic effect of α-synuclein.
Furthermore, we showed that transfection of p497_25g plasmid significantly downregulated the expression levels of OPA1 and VAPB in two cell lines (Figure 5a-d). Additionally, tranfection of miR-497_25g had significant inhibitory effects on the proliferation of U-118-MG and Hs683 cells compared to controls (Figure 5e-j).
In summary, hsa-mir-497 was edited at a normal level and A-to-I edited miR-497-5p repressed OPA1 and VAPB to induce neorodegeneration at a normal pace in healthy PC samples (Figure 6a). In comparison, the editing level of hsa-mir-497_25_A_g was significantly increased in PD-PC samples and the A-to-I edited miR-497-5p dominantly repressed OPA1 and VAPB, which led to significant downregulation of OPA1 and VAPB in PD-PC (Figure 6b). Consequently, the significant repression of OPA1 and VAPB contributed to the enhanced neurodegeneration in PD. These results suggest that novel therapies of PD might be designed by either over-expressing OPA1/VAPB or repressing the editing level of hsa-mir-497_25_A_g. Furthermore, the significant correlation between ages of normal people and editing level of hsa-mir-497_25_A_g, as well as expression levels of OPA1 and VAPB, and the significantly increased editing level of hsa-mir-497_25_A_g in PD might be used to develop novel diagnostic methods of PD. In (a-j), p-values were based on two-tailed t-tests. Figure 6. The role of hsa-mir-497_25_A_g in PD. (a) In PC tissue, the A-to-I editing of hsa-mir-497 is performed at a normal level. The miR-497_25g with a normal expression level represses OPA1 and VAPB to induce neurodegeneration at a normal pace. (b) The A-to-I editing of hsa-mir-497 increases in PD-PC tissue. The higher expression level of miR-497_25g leads to a significant downregulation of OPA1 and VAPB, which consequently contributes to enhanced neurodegeneration in PD-PC tissue.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/cells12010075/s1, Figure S1: The distribution of different types of editing sites in miRNAs; Figure S2: The details of 23 identified A-to-I editing sites in miRNAs; Figure S3: The details of 4 identified C-to-U editing sites in miRNAs; Figure S4: The details of 14 identified SNPs in miRNAs; Figure S5: The distribution of PD relevant editing sites; Figure S6: The editing levels of five A-to-I and four other types of M/E sites in PC samples and PD-PC-like samples; Figure S7: The expression of OPA1 in different brain regions of PD patients and normal controls; Figure S8: The expression of VAPB in different brain regions of PD patients and normal controls; Figure S9: Validation of monkey OPA1 and VAPB as targets of A-to-I edited miR-497-5p; Table S1: Data sets analyzed or used in this study; Table S2: The 860 editing sites that have significant editing levels in the selected sRNA-seq profiles; Table S3: The 276 editing sites whose editing levels in PC and PD-PC samples are significantly different; Table S4: Targets of four A-to-I edited miRNAs and their corresponding original miRNAs; Table S5: GO and KEGG pathway analysis of targets of four A-to-I edited miRNAs and their corresponding original miRNAs; Table S6: Information of primers used in this study; Table S7: Source data of Figure 1; Table S8: Source data of Figure 2; Table S9: Source data of Figure 3; Table S10: Source data of Figure 4; Table S11: Source data of Figure 5; Table S12: Source data of Figure S1; Table S13: Source data of Figure S2; Table S14: Source data of Figure S3; Table S15: Source data of Figure S4; Table S16: Source data of Figure S5; Table S17: Source data of Figure S6; Table S18: Source data of Figure S7; Table S19: Source data of Figure S8; and Table S20: Source data of Figure S9.

Conflicts of Interest:
The authors declare that they have no competing interests.