Characterizing relevant microRNA editing sites in Parkinson’s disease
=======================================================================

* Chenyu Lu
* Shuchao Ren
* Zhigang Zhao
* Xingwang Wu
* Angbaji Suo
* Nan Zhou
* Jun Yang
* Shuai Wu
* Tianqing Li
* Chao Peng
* Yun Zheng

## Abstract

MicroRNAs (miRNAs) are extensively edited in human brains. However, the functional relevance of miRNA editome is largely unknown in Parkinson’s disease (PD). By analyzed small RNA sequencing profiles of brain tissues of 43 PD patients and 88 normal controls, we totally identified 421 miRNA editing sites with significantly different editing levels in prefrontal cortices of PD patients (PD-PC). A-to-I edited miR-497-5p has significantly higher expression levels in PD-PC compared to normal controls and directly represses OPA1 and VAPB, which potentially contributes to the progressive neurodegeneration of PD patients. These results provide new insights into mechanistic understanding, novel diagnostic and therapeutic clues of PD.

Key words
*   RNA editing
*   Parkinson’s disease
*   hsa-miR-497-5p
*   OPA1
*   VAPB

## Introduction

MicroRNAs (miRNAs) are small non-coding RNAs with about 22 nucleotides that normally repress their target mRNAs at post-transcriptional level [1]. Some miRNAs are edited, such as Adenosine-to-Inosine (A-to-I) editing [2–7] performed by adenosine deaminase (ADAR) enzymes and C-to-U editing performed by apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC) enzymes [8], during their biogenesis processes. Altered editing of miRNAs lead to human diseases, such as cancer [9–11]. Parkinson’s disease is a neurodegenetative disorder affecting 2-3% of people over 65 years of age [12]. However the functional relevance of miRNA editing in PD is largely unknown, although A-to-I editing is prevalent in brain [13–16] and the editing level is gradually increasing in the developmental procedure [5,16].

To comprehensively characterize miRNA editing sites in brain tissues of PD, we analyzing 43 and 88 small RNA sequencing profiles of PD and control brain tissues, respectively. We identified 421 miRNA editing sites that have significantly different editing levels in prefrontal cortices of PD patients (PD-PC) compared with control group. The editing levels of 6 A-to-I and 3 C-to-U editing sites are significantly correlated with the ages of normal controls which is disrupted in PD patients. One A-to-I editing site in miR-497-5p has significantly higher editing levels compared to normal controls in PD-PC samples and edited miR-497-5p directly represses OPA1 and VAPB, which potentially contributes to the progressive neurodegeneration of PD patients. These results demonstrate that miRNA editing is severely disturbed and relevant in PD and offer novel insight into the etiology of PD. The editing of miRNAs might be used to develop novel diagnostic biomarkers and/or therapeutic targets for PD.

## Results

### Summary of data sets 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 NCBI SRA Database (Supplementary Table S1.1). These profiles include 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). To understand the potential function of edited miRNAs, we examined the deregulated genes and proteins from the BA9 and several other brain regions of PD patients (Table S1.2).

### An overview of identified editing sites in miRNAs

We used the MiRME pipeline [17] to analyze the 131 sRNA-seq profiles selected with the default settings. Totally, we found 2456 significant M/E sites (as shown in Figure 1 and Table S2.1), by the criteria of at least 10 reads and multiple test corrected *P* -values smaller than 0.05. The largest category of these M/E sites are 3’-A, accounting for 37.3% (916 sites), then followed by 3’-U (28.6%), 3’-Other (15.5%) and 5’-editing (5.4%) (Fig. 1a). A-to-I and C-to-U editing sites account for 1.9% (46 sites) and 1.1% (26 sites), respectively (Fig. 1a). We found 60 SNPs after comparing the identified M/E sites to SNPs in dbSNP (Fig. 1a). Other types of editing sites consist of 6.7% (164 sites) (Fig. 1a).

![Fig. 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/12/04/2020.04.06.20054791/F1.medium.gif)

[Fig. 1.](http://medrxiv.org/content/early/2020/12/04/2020.04.06.20054791/F1)

Fig. 1. 
A summary of the identified miRNA mutation and editing sites. (a) The categories of significant M/E sites in miRNAs. (b) The numbers of different types of editing sites that have significant Spearman correlation (*ρ*) between editing levels and the ages of individuals in the PC and PD-PC group. (c) The distributions of *ρ* for A-to-I and C-to-U sites in the PC and PD-PC samples, respectively. (d) Ten selected editing sites with significant *ρ* values in either PC or PD-PC samples. *: *P* < 0.05; **: corrected *P* < 0.05. (e) to (f) The plots of ages of PC samples against the editing levels of two selected editing sites. See also Extended Data Fig. 1.

The A-to-I, C-to-U and Other sites were further classified based on their types of variations (Extended Data: Fig. 1a). Similar to previous work [17], there could be several 3’-editing sites in one pre-miRNA, but most pre-miRNAs only have 1 or 2 central sites in our selected samples (ED: Fig. 1b). However, some miRNAs may have a few editing events at 5’ end in our previous study [17].

### Age-related editing sites in PD and normal controls

Because previous studies reported that A-to-I editing level is increasing during the development process [5,16], we examined the Spearman correlation (*ρ*) between editing levels of the 2456 editing sites and the ages (at death) of 36 PC and 29 PD-PC samples, respectively. There are 411 sites with significant *ρ* in PC, but only 134 sites with significant *ρ* in PD-PC samples (Fig. 1b and Table S2.2-2.3). When these sites were separated as ones with positive and negative *ρ* values, the numbers of sites with significant *ρ* values in PD-PC samples are also much smaller than those in PC samples (Fig. 1b and ED: Fig. 1c). The numbers of types of editing sites are decreased in PD-PC samples too. For example, there are no A-to-I sites with positive *ρ* values in the PD-PC samples (ED: Fig. 1c). Furthermore, the Kullback-Leibler divergences of the distributions of different types of editing sites in PC and PD-PC samples are 0.6, 1.9 and 1.7 for all, sites with positive *ρ* values and negative *ρ* values, respectively (ED: Fig. 1c). These results suggest that the gradually increasing or decreasing editing levels of miRNAs during the aging procedure are severely disturbed in PD.

The *ρ* values of 47 A-to-I and 26 C-to-U editing sites are carefully examined in Fig. 1c. In PC samples, there are 4 and 2 A-to-I sites with significant positive and negative *ρ* values, respectively (Fig. 1c-1e). In comparison, only one A-to-I site has significant negative *ρ* in PD-PC (Fig. 1c-1d). The median value of *ρ* in PC is slightly larger than 0, which is consistent with ever-increasing A-to-I editing level along age noticed previously [5,16]. But the median *ρ* value in PD-PC samples is smaller than 0, further suggesting the disrupted A-to-I editing of miRNAs in PD. In contrast, the C-to-U editing sites have a weak negative median *ρ* value in PC, suggesting C-to-U editing sites of miRNAs generally have decreasing editing levels when people are aging. Three C-to-U editing sites (see Fig. 1c, 1d and 1f) have significant *ρ* value in PC.

### A-to-I editing sites

We totally found 47 significant A-to-I editing sites from the selected samples (ED:Figure 2a and Table S2.4). Among them, 14 sites are newly identified when compared with reported sites (ED: Fig. 2a). At least 19 of these 47 sites are conserved in other mammals [3,18] (ED: Fig. 2a). Similar to previous results [3,4,17,18], the 5’ and 3’ side of these A-to-I editing sites preferred to be U and G, respectively (ED: Fig. 2b). For examples, two A-to-I editing sites were shown in ED: Fig. 2c to 2f.

![Fig. 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/12/04/2020.04.06.20054791/F2.medium.gif)

[Fig. 2.](http://medrxiv.org/content/early/2020/12/04/2020.04.06.20054791/F2)

Fig. 2. 
M/E sites that have significantly different editing levels in PD-PC. (a) Five A-to-I editing sites with significantly different editing levels in PD-PC. (b) Eight editing sites with significantly different editing levels in PD-PC. *: corrected *P* < 0.05, **: corrected *P* < 0.01, \***|: corrected *P* < 0.001, Mann-Whitney *U* -test. See also Extended Data Fig. 6a-6b.

### C-to-U editing sites

We carefully explored the 26 C-to-U editing sites (ED: Fig. 3a and Table S2.5). Five of these 26 C-to-U editing sites are conserved in primates [18] (see ED: Fig. 3a). The neighboring nucleotides of these C-to-U sites prefer to be C on both 5’ and 3’ sides (ED: Fig. 3b), consistent with the CCC motif of APOBEC3G [19]. Three examples of C-to-U editing sites are shown in ED: Fig. 3c-3h. There are multiple C-to-U editing sites on two pre-miRNAs, i.e., hsa-mir-219a-2 and hsa-mir-1249 (ED: Fig. 4 and Table S2.5). Interestingly, seven continuous cytosines on hsa-mir-1249 were found to be C-to-U edited (ED: Fig. 3a and Fig. 4).

![Fig. 3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/12/04/2020.04.06.20054791/F3.medium.gif)

[Fig. 3.](http://medrxiv.org/content/early/2020/12/04/2020.04.06.20054791/F3)

Fig. 3. 
The target analysis of A-to-I edited hsa-miR-497-5p. (a) Integrated analysis of targets of edited hsa-miR-497-5p (hsa-mir-497_25g) and the downregulated genes and proteins in PD-PC (BA9). (b) The complementary sites of hsa-mir-497_25g on OPA1 and VAPB, and the PAR-CLIP sequencing reads from these sites. (c) Comparison of abundances of hsa-mir-497_25g in PC and PD-PC. (d) - (e) Comparison of abundances of OPA1 and VAPB in PC and PD-PC (BA9). *: *P* < 0.05; **: corrected *P* < 0.05; \***|: corrected *P* < 0.001, DESeq2 and limma package, respectively. (f) - (g) The plots of the ages of PC samples against the abundances of OPA1 and VAPB. See also Extended Data Fig. 7 and 8.

![Fig. 4.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/12/04/2020.04.06.20054791/F4.medium.gif)

[Fig. 4.](http://medrxiv.org/content/early/2020/12/04/2020.04.06.20054791/F4)

Fig. 4. 
Validating selected targets of A-to-I edited miR-497-5p. (a) - (b) The segment of 3’ UTR of OPA1 and VAPB, mutated segment of 3’ UTR of OPA1 and VAPB, segment of 3’ UTR of mouse Opa1 and Vapb, and conservation scores for 30 mammals in these region calculated with PhyloP, resepctively. (c) - (d) The luciferase activities when co-transfecting a pGLO plasmid of 3’ UTRs in Part (a) - (b) and a pCDNA plasmid containing original pre-hsa-mir-497 (p497) or pre-hsa-mir-497_25g (p497 25g), respectively. The values shown are mean values *±* SDs. *P* -values were based on two-tailed *t*-tests. See also Extended Data Fig. 9.

### Identified SNPs in miRNAs

By comparing the M/E sites to SNPs reported in dbSNP and examining their editing levels, 60 M/E sites were regarded as SNPs (as shown in ED: Fig. 5a and Table S2.6). Three of the 60 SNP sites were shown in ED: Fig. 5b to 5g. The editing levels of these three sites are 100% in the selected samples (ED: Fig. 5e to 5g).

### Relevant miRNA editing sites in PD

To identify editing sites that are relevant in PD, we compared the editing levels of 2456 significant editing sites in PD-PC and PC; PD-Am and Am samples, respectively. We obtained 421 M/E sites that have significantly different editing levels in PD-PC samples compared to PC samples (ED: Fig. 6a and Table S3). Most of 421 sites (65.8% or 277 sites) have increased editing levels in PD-PC samples, and 34.2% of these sites have decreased editing levels in PD-PC samples (ED: Fig. 6b). 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) have 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) show decreased editing levels in PD-PC samples (Fig. 2a). Four C-to-U and four Other editing sites have decreased editing levels in PD-PC samples (Fig. 2b). There are no M/E sites with significantly different editing levels in the comparison between PD-Am and Am samples.

### Target analysis of A-to-I editing sites in seeds

We selected 4 M/E sites that locate in the seed regions of mature miRNAs (Table S3) and predicted their targets with the MiCPAR pipeline [20], as listed in Table S4.1. After comparing the results to those of original miRNAs (in Table S4.2), we found that these editing events severely changed the target sets of these miRNAs (ED: Fig. 6c), as well as the enriched GO terms and KEGG pathways of these targets (ED: Fig. 6c and Table S5).

Since A-to-I edited hsa-miR-497-5p (hsa-mir-497_25g) has much higher expression level than the other three A-to-I edited miRNAs (ED: Fig. 6d), 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 (Fig. 3a, see details in Methods), we found that hsa-mir-497_25g targets OPA1 mitochondrial dynamin like GT-Pase (OPA1) and VAMP associated protein B and C (VAPB) (Fig. 3b). Presumably due to the upregulation of hsa-mir-497_25g (Fig. 3b) in PD-PC, OPA1 and VAPB are significantly downregulated in PD-PC (Fig. 3d-3e). Furthermore, the protein expression level of OPA1 and VAPB is also significantly downregulated in PD-PC samples (corrected *P* < 0.05, as in [21]).

Furthermore, OPA1 is significantly downregulated in substantia nigra (SN) (ED: Fig. 7a-7b), is mildly decreasing in putamen of PD patients (ED: Fig. 7c). But in subventricular zone (SZ) and cingulate gyrus, the expression of OPA1 does not change significantly (ED: Fig. 7d-7e). In SN, VAPB has a weak decreasing trend in PD (ED: Fig. 7a), but has no change in another cohort (ED: Fig. 7b), and cingulate gyrus, SZ, putamen (ED: Fig. 7c to 7e, respectively).

Importantly, consistent with the significant positive correlation between editing level of hsa-mir-497\_25\_A_g and ages of individuals in PC (Fig. 2f), the expression levels of OPA1 and VAPB are significantly negatively correlated with the ages of individuals in PC (Fig. 3f, 3g and ED: Fig. 7f/8f) and potentially in SN as well (*ρ* = *−*0.36 and *−*0.47, respectively, see ED: Fig. 7g and 8g, respectively). The Spearman correlation *ρ* between expression of OPA1 and VAPB and ages in putamen are weak and insignificant in putamen (ED: Fig. 7h and 8h, respectively), and are positive in SZ (*ρ* = 0.29 and 0.60 in ED: Fig. 7i and 8i, respectively). In PD, these correlations were generally weakened (ED: Fig. 7j-7l/7n and 8j-8l/8n for OPA1 and VAPB, respectively), except in putamen where the negative correlations are slightly enhanced (ED: Fig. 7m and 8m for OPA1 and VAPB, respectively).

### hsa-mir-497_25g directly represses OPA1 and VAPB

We first examined the conservation of the complementary site of hsa-mir-497_25g on VAPB. The 7mer opposite to the seed of edited miR-497-5p is only partially conserved in mammals (Fig. 4a). Especially, the 3rd nucleotide from 5’ end (as pointed by a red arrow in Fig. 4a) has a conservation score of -1.47, indicating that the site is fast-evolving in primates. Even within primates, 7 species did not have conserved sequences complementary to the seed of hsa-mir-497_25g (Fig. 4b).

To verify that hsa-mir-497\_25g directly represses OPA1 and VAPB, we performed luciferase assays in 293T cells for hsa-mir-497\_25g, original hsa-mir-497, OPA1/VAPB, mutated OPA1/VAPB (see Fig. 4a and 4b, respectively). As in Fig. 4c and 4d, hsa-mir-497_25g induced significant decreases of luciferase activities for OPA1 and VAPB, indicating that hsa-mir-497_25g directly represses OPA1 and VAPB, respectively. In comparison, when being co-transfected with mutated OPA1 and mutated VAPB, hsa-mir-497_25g could not repress mutated OPA1 and mutated VAPB (Fig. 4c and 4d, respectively). The original miR-497-5p could not repress OPA1 and VAPB too (Fig. 4c and 4d, respectively). The complementary sites of hsa-mir-497_25g on OPA1 and VAPB have non-conserved seed regions (see red rectangles in Fig. 4a-4b). Actually, these two hsa-mir-497_25g complementary sites are only conserved in some primates including monkey, but non-conserved in rodents (ED: Fig. 9a-9d). As expected, A-to-I edited miR-497-5p could repress monkey OPA1 and VAPB, but could not repress mouse Opa1 and Vapb (ED: Fig. 9e-9h).

To summarize, hsa-mir-497\_25g directly represses OPA1 and VAPB in human, and putatively in some primates, but not in mouse and those primates with non-conserved complementary sites of hsa-mir-497\_25g, which might contributes to the advanced functions and fast evolving of prefrontal cortices in human.

## Discussion

*OPA1* localizes to the inner mitochondrial membrane and plays roles in regulating mito-chondrial stability and energy output. Mutations in OPA1 have been associated with optic atrophy type 1 [22,23]. Recently, the association between OPA1 and parkinsonism is also being noticed [24,25]. One recent study found that 3 members of a family had autosomal dominant optic atrophy caused by OPA1 mutations and 2 of them developed nonsyndromic PD [26]. OPA1 is downregulated in PD-PC (BA9) (Fig. 2d-2e) and in substantia nigra of PD patients (ED: Fig. 7b-7c). PD patients carrying the G2019S mutation showed decreased levels of mature OPA1 [27]. OPA1 haploinsufficinecy leads to progressive loss of iPSC-derived dopaminergic neurons [28]. Autosomal dominant LRRK2 mutations are associated with both familial and sporadic PD [27,29]. LRRK2 directly interacts with OPA1, presumably to regulate membrane dynamics important for mediating endocytosis and mitochondrial function [27].

Neuronal synaptic dysfunction is one of the key features of Parkinsons disease [30,31]. VAPB protein is a membrane protein found in plasma and intracellular vesicle membranes, and involved in vesicle trafficking. VAPB is downregulated in BA9 of PD-PC (Fig. 3d-3e). Down-regulation of VAPB and its interacting protein PTPIP51 reduces dentritic spine numbers and synaptic activity [32]. A missense mutation (P56S) in VAPB was reported to associated with motorneuron degeneration in affected amyotrophic lateral sclerosis (ALS) patients [33]. Transgenic mice that expressed human P56S *VAPB* developed progressive hyperactivities and other motor abnormalities, and showed progressive loss of corticospinal motor neurons [34]. In contrast, overexpression of VAPB in mouse ALS model slows the motor impairment and promotes the survival of spinal motor neurons [35].

α-synuclein is a PD-related protein and localizes to the nerve terminal [36] and directly binds to VAPB [36] and VAMP2 [37]. Disruption of the interaction between VAPB and PTPIP5 has been suggested to mediated the neurotoxic effect of α-synuclein. The downregulation of VAPB in PD patients might exaggerate the toxic effect of α-synuclein.

In summary, our results suggest that the editing level of hsa-mir-497\_25\_A\_g increases in PD-PC samples and the A-to-I edited miR-497-5p dominantly represses OPA1 and VAPB in PD-PC samples, which leads to significant downregulation of OPA1 and VAPB in PD-PC. Consequently, the reduced OPA1 and VAPB contributes to the progressive neurodegeneration in PD. These results suggest that novel therapies of PD might be designed by either over-expressing OPA1/VAPB or repressing 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 and therapies of PD.

## Methods

### The small RNA sequencing profiles used

As summarized in Supplementary Table S1.1, we collected 131 sRNA-seq profiles of post-mortem brain samples of PD patients and normal controls, including 29 prefrontal cortex samples of PD patients, 14 amygdalae samples of PD patients, 36 prefrontal cortex samples of normal controls, 15 amygdalae samples of normal controls, 6 frontal cortex samples of normal controls, 6 corpus callosum samples of normal controls, 3 astrocyte cell line samples of normal controls, 4 inferior parietal lobe samples of normal controls, 2 temporal neocortex gray matter samples of normal controls, and 17 unknown region samples of normal controls.

### 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 summarized in Table S1.2, we chose 107 and 122 gene expression profiles of PD patients and normal controls, respectively, from 7 cohorts, including 2 cohorts for prefrontal cortex (Brodmann Area 9, briefly BA9) (GSE20168 and GSE68719), 2 cohorts for substantia nigra (GSE7621 and GSE20292), 1 cohort for putamen (GSE20291), 1 cohort for cingulate gyrus (GSE110716), and 1 cohort for subventricular zone (GSE130752).

We also checked 283 deregulated proteins in prefrontal contex (BA9) samples of PD by comparing the MS3 proteomics profiles of 12 PD patients to those of 12 normal people reported previously [21].

### Genome and annotation of miRNAs used

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

### Analysis of small RNA sequencing profiles

The selected sRNA-seq profiles were analyzed using the MiRME pipeline [17] with the default settings. The following criteria were used to define M/E sites with significant editing levels: (i) the relative level of editing is at least 5%; (ii) at least 10 reads support the editing event; (iii) the score threshold of sequencing reads is 30; and (iv) a multiple-test corrected *P* -value (using the Benjamini and Hochberg method [41]) of smaller than 0.05. Then, the obtained results of different samples were combined by a separate program in the MiRME package (see details in [17]). The identified M/E sites were compared to reported editing sites in miRNAs in the DARNED database [42], the RADAR database [43] and literature [3,44,4,45,46,17]. 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 sites, the nucleotides from the reference pre-miRNA sequences and the edited/mutated nucleotide at the sites.

### Comparing the M/E sites to reported SNPs

The identified M/E sites were compared to known SNPs in miRNAs organized in [47] (which was based on the dbSNP v137) and reported SNPs in the dbSNP (v151). Only sites that satisfied the following criteria were regarded as SNPs, (i) 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* [18] and *Mus musculus* [2,3]. The editing sites of the same editing types that locate 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 2456 editing sites in Table S1 and the age of death (y) was calculated with the corr function in MatLab (Mathworks, MA) for the 29 PD-PC and 36 PC samples, respectively. The editing sites with *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 2456 M/E sites in the PD-PC and PC; PC-Am and AM samples were compared, respectively, with the Mann-Whitney *U* -test. The obtained *P* -values were corrected with the Benjamini-Hochberg correction method [41]. M/E sites with multiple test corrected *P* -value smaller than 0.05 were regarded as having significantly different editing levels in PD-PC or PC-Am samples.

### Identifying targets for original and edited miRNAs

M/E sites that satisfy the following criteria were chosen to identify the targets of original and edited miRNAs, (i) have significantly different editing levels in PD-PC samples compared to PC samples; and (ii) locate in the seed regions of mature miRNAs. The targets of original and edited miRNAs were predicted using the MiCPAR algorithm [20] with its default parameters. The targets with at least 1 PAR-CLIP read with T-to-C variation were kept in further analysis.

As listed in Table S1.3, seven PAR-CLIP sequencing profiles prepared from HEK293 cells stably expressing FLAG/HA-tagged AGO proteins (AGO1, AGO2, AGO3, and AGO4) [48] 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 [49] 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 have 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 [20]. 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 genome sequences of human (version GRCh38) were downloaded from the UCSC Genome Browser [50] and used as inputs of the MiCPAR algorithm. The targets having at least 1 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 term and KEGG pathway enrichment of the targets that were only targeted by the original or edited miRNAs were analyzed with KOBAS2, respectively [51]. The significantly enriched GO terms (with multiple test corrected *P* -values smaller than 0.05) were divided into 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 [21] identified 1095 deregulated genes with the Wald tests (in DESeq2 [52]) by analyzing RNA-seq data sets for 19 PD and 24 control samples of 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 [21]. Another study [53] generated 30 gene expression profiles of BA9 with microarray, with 15 profiles for PD and PD-PC respectively. The deregulated genes were identified with the limma package.

Because these RNA-seq and proteomics profiles were also from BA9 which is 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 have 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.

### Validating selected targets of edited hsa-miR-497-5p

The 3’ UTR segments of human OPA1/VAPB, mutated OPA1/VAPB, crab-eating monkey (*Macaca fascicularis*) OPA1/VAPB, mutated monkey OPA1/VAPB and mouse (*Mus musculus*) Opa1/Vapb containing the hsa-mir-497_25g complementary sites (*∼* 130 nt with 70 nt on both sides) were synthesized and cloned into the plasmid pmirGLO (Promega, Madison, WI, US) with XholI and SalI sites, and named as pGLO-OPA1/VAPB, pGLO-OPA1m/VAPBm, pGLO-mfaOPA1/mfaVAPB, pGLO-mfaOPA1m/mfaVAPBm, and pGLO-musOpa1/musVapb, respectively. The genomic regions of hsa-pre-mir-497 and hsa-pre-mir-497 25g were synthesized and cloned into the plasmid pCDNA 3.1(+) (Invirtrogen, Carlsbad, CA, USA) with HindIII and BamHI sites, and named as p497 and p497 25g, respectively. All the plasmids were validated by Sanger sequencing (ABI3730, Thermo Fisher, Waltham, MA, US).

Human 293T cells was cultured in DMEM high glucose medium with 10% FBS, 1% NEAA. 500 µL medium was inoculated to each well with 24 well plate. 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 added and mixed 2 µL of DNA-INVI DNA Transfection Reagent (Invigentech, Carlsbad, CA 92008, USA) in each well.

In the experiments, each group had three biological replicates. 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 for 30 mammals) of the A-to-I edited miR-497-5p complementary sites on OPA1 and VAPB were downloaded from UCSC Genome Browser [38]. 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) [54], respectively. The obtained phylogenetic trees were visualized with TreeView (v1.6.6) [55].

## Supporting information

Extended Data: Fig. 1 to Fig. 9 [[supplements/054791_file02.pdf]](pending:yes)

Supplementary Table S9 [[supplements/054791_file03.xlsx]](pending:yes)

Supplementary Table S18 [[supplements/054791_file04.xlsx]](pending:yes)

Supplementary Table S1 [[supplements/054791_file05.xlsx]](pending:yes)

Supplementary Table S2 [[supplements/054791_file06.xlsx]](pending:yes)

Supplementary Table S3 [[supplements/054791_file07.xlsx]](pending:yes)

Supplementary Table S4 [[supplements/054791_file08.xlsx]](pending:yes)

Supplementary Table S5 [[supplements/054791_file09.xlsx]](pending:yes)

Supplementary Table S6 [[supplements/054791_file10.xlsx]](pending:yes)

Supplementary Table S7 [[supplements/054791_file11.xlsx]](pending:yes)

Supplementary Table S8 [[supplements/054791_file12.xlsx]](pending:yes)

Supplementary Table S10 [[supplements/054791_file13.xlsx]](pending:yes)

Supplementary Table S11 [[supplements/054791_file14.xlsx]](pending:yes)

Supplementary Table S12 [[supplements/054791_file15.xlsx]](pending:yes)

Supplementary Table S13 [[supplements/054791_file16.xlsx]](pending:yes)

Supplementary Table S14 [[supplements/054791_file17.xlsx]](pending:yes)

Supplementary Table S15 [[supplements/054791_file18.xlsx]](pending:yes)

Supplementary Table S16 [[supplements/054791_file19.xlsx]](pending:yes)

Supplementary Table S17 [[supplements/054791_file20.xlsx]](pending:yes)

## Data Availability

The data used in this study were available at the NCBI SRA and GEO databases at https://www.ncbi.nlm.nih.gov/sra and https://www.ncbi.nlm.nih.gov/geo/.

## Data Availability

The 131 sRNA-seq profiles, 7 cohorts of gene expression profiles, and 11 PAR-CLIP sequencing profiles were available at the NCBI GEO and SRA databases under the accession numbers listed in Table S1.

## Authors’ contributions

YZ conceived and designed the research. CL, SR, ZZ, XW, AS, NZ, JY, CP, and YZ analyzed the data and organized the results. CL, XW, SW and TL performed the luciferase experiments. YZ developed the computational pipelines for analyzing the data sets. YZ wrote the manuscript. All authors had read and approved this paper in the current form.

## Competing interests

The authors declare that they have no competing interests.

## Additional information

Supplementary information is available for this paper at the website of the journal. Correspondence and requests for materials should be addressed to YZ.

## Acknowledgements

The research was supported in part by a grant (No. 31760314) of National Natural Science Foundation of China ([http://www.nsfc.gov.cn/](http://www.nsfc.gov.cn/)) to YZ and a grant (No. 2018YFA0108502) of the Ministry of Science and Technology of China to YZ. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

## Footnotes

*   Add a new author, Shuai Wu, was added. Update name of author Xingwang Wu. Fig. 4 was updated. A new figure Fig. 9 was added in Extended Data. The related parts in Methods were also updated. Table S9 was updated. A new Supplementary Table, Table S19, was added.

*   Received April 6, 2020.
*   Revision received December 3, 2020.
*   Accepted December 4, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  [1].Bartel, D. P. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 116, 281–297 (2004).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0092-8674(04)00045-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14744438&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000188541700017&link_type=ISI) 

2.  [2].Kawahara, Y., Zinshteyn, B., Sethupathy, P., Iizasa, H., Hatzigeorgiou, A. G., and Nishikura, K. Redirection of silencing targets by Adenosine-to-Inosine editing of miRNAs. Science 315(5815), 1137–1140 (2007).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMTUvNTgxNS8xMTM3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMDQvMjAyMC4wNC4wNi4yMDA1NDc5MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

3.  [3].Kawahara, Y., Megraw, M., Kreider, E., Iizasa, H., Valente, L., Hatzigeorgiou, A., and Nishikura, K. Frequency and fate of microRNA editing in human brain. Nucleic Acids Res. 36(16), 5270–5280 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkn479&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18684997&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000259205600021&link_type=ISI) 

4.  [4].Alon, S., Mor, E., Vigneault, F., Church, G. M., Locatelli, F., Galeano, F., Gallo, A., Shomron, N., and Eisenberg, E. Systematic identification of edited microRNAs in the human brain. Genome Res. 22(8), 1533–1540 (2012).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjk6IjIyLzgvMTUzMyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzEyLzA0LzIwMjAuMDQuMDYuMjAwNTQ3OTEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

5.  [5].Ekdahl, Y., Farahani, H., Behm, M., Lagergren, J., and Öhman, M. A-to-I editing of microRNAs in the mammalian brain increases during development. Genome Res. 22(8), 1477–1487 (2012).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjk6IjIyLzgvMTQ3NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzEyLzA0LzIwMjAuMDQuMDYuMjAwNTQ3OTEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

6.  [6].Nishikura, K. A-to-I editing of coding and non-coding RNAs by ADARs. Nature Reviews. Molecular cell biology 17, 83–96 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrm.2015.4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26648264&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 

7.  [7].Behm, M. and Ohman, M. Rna editing: A contributor to neuronal dynamics in the mammalian brain. Trends in Genetics 32, 165–175, 01 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tig.2015.12.005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26803450&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 

8.  [8].Negi, V., Paul, D., Das, S., Bajpai, P., Singh, S., Mukhopadhyay, A., Agrawal, A., and Ghosh, B. Altered expression and editing of miRNA-100 regulates iTreg differentiation. Nucleic Acids Research 43(16) (2015).
    
    

9.  [9].Shoshan, E., Mobley, A. K., Braeuer, R. R., Kamiya, T., Huang, L., Vasquez, M. E., Salameh, A., Lee, H. J., Kim, S. J., Ivan, C., Velazquez-Torres, G., Nip, K. M., Zhu, K., Brooks, D., Jones, S. J. M., Birol, I., Mosqueda, M., Wen, Y.-y., Eterovic, A. K., Sood, A. K., Hwu, P., Gershenwald, J. E., Robertson, A. G., Calin, G. A., Markel, G., Fidler, I. J., and Bar-Eli, M. Reduced adenosine-to-inosine miR-455-5p editing promotes melanoma growth and metastasis. Nature Cell Biol. 17(3), 311–321, February (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncb3110&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25686251&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 

10. [10].Pinto, Y., Buchumenski, I., Levanon, E. Y., and Eisenberg, E. Human cancer tissues exhibit reduced A-to-I editing of miRNAs coupled with elevated editing of their targets. Nucleic Acids Research 46(1), 71–82, 11 (2017).
    
    

11. [11].Wang, Y., Xu, X., Yu, S., Jeong, K. J., Zhou, Z., Han, L., Tsang, Y. H., Li, J., Chen, H., Mangala, L. S., Yuan, Y., Eterovic, A. K., Lu, Y., Sood, A. K., Scott, K. L., Mills, G. B., and Liang, H. Systematic characterization of A-to-I RNA editing hotspots in microRNAs across human cancers. Genome Research 27(7), 1112–1125 (2017).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjk6IjI3LzcvMTExMiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzEyLzA0LzIwMjAuMDQuMDYuMjAwNTQ3OTEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

12. [12].Williams-Gray, C. H. and Worth, P. F. Parkinson’s disease. Medicine 44(9), 542–546 (2016).
    
    

13. [13].Eisenberg, E., Nemzer, S., Kinar, Y., Sorek, R., Rechavi, G., and Levanon, E. Is abundant A-to-I RNA editing primate-specific? Trends in genetics : TIG 21, 77–81, 03 (2005).
    
    

14. [14].Paz, N., Levanon, E., Nevo, E., Kinar, Y., Harmelin, A., Jacob-Hirsch, J., Amariglio, N., Eisenberg, E., and Rechavi, G. Adenosine-to-inosine RNA editing shapes transcriptome diversity in primates. Proceedings of the National Academy of Sciences of the United States of America 107, 12174–9, 07 (2010).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTA3LzI3LzEyMTc0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMDQvMjAyMC4wNC4wNi4yMDA1NDc5MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

15. [15].Li, J. B. and Church, G. M. Deciphering the functions and regulation of brain-enriched A-to-I RNA editing. Nat Neurosci 16, 1518–1522 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nn.3539&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24165678&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 

16. [16].Zaidan, H., Ramaswami, G., Golumbic, Y., Sher, N., Malik, A., Barak, M., Galiani, D., Dekel, N., Li, J., and Gaisler-Salomon, I. A-to-I RNA editing in the rat brain is age-dependent, region-specific and sensitive to environmental stress across generations. BMC Genomics 19, 12 (2018).
    
    

17. [17].Zheng, Y., Ji, B., Song, R., Wang, S., Li, T., Zhang, X., Chen, K., Li, T., and Li, J. Accurate detection for a wide range of mutation and editing sites of microRNAs from small RNA high-throughput sequencing profiles. Nucleic Acids Research 44(14), e123 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkw471&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27229138&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 

18. [18].Wang, Q., Zhao, Z., Zhang, X., Lu, C., Ren, S., Li, S., Guo, J., Liao, P., Jiang, B., and Zheng, Y. Identifying microRNAs and their editing sites in Macaca mulatta. Cells 8(7) (2019).
    
    

19. [19].Chen, J. and MacCarthy, T. The preferred nucleotide contexts of the AID/APOBEC cytidine deaminases have differential effects when mutating retrotransposon and virus sequences compared to host genes. PLoS computational biology 13, e1005471, 03 (2017).
    
    

20. [20].Zheng, Y. Computational Non-coding RNA Biology. Academic Press, Cambridge, MA, 1st edition, (2018).
    
    

21. [21].Dumitriu, A., Golji, J., Labadorf, A., Gao, B.G. Beach, T., Myers, R.A. Longo, K., and Latourelle, J. Integrative analyses of proteomics and RNA transcriptomics implicate mitochondrial processes, protein folding pathways and GWAS loci in Parkinson disease. BMC Medical Genomics 9, 12 (2015).
    
    

22. [22].Hudson, G., Amati-Bonneau, P., Blakely, E. L., Stewart, J. D., He, L., Schaefer, A. M., Griffiths, P. G., Ahlqvist, K., Suomalainen, A., Reynier, P., McFarland, R., Turnbull, D. M., Chinnery, P. F., and Taylor, R. W. Mutation of OPA1 causes dominant optic atrophy with external ophthalmoplegia, ataxia, deafness and multiple mitochondrial DNA deletions: a novel disorder of mtDNA maintenance. Brain 131(2), 329–337, 12 (2007).
    
    

23. [23].Amati-Bonneau, P., Valentino, M. L., Reynier, P., Gallardo, M. E., Bornstein, B., Boissire, A., Campos, Y., Rivera, H., de la Aleja, J.G., Carroccia, R., Iommarini, L., Labauge, P., Figarella-Branger, D., Marcorelles, P., Furby, A., Beauvais, K., Letournel, F., Liguori, R., La Morgia, C., Montagna, P., Liguori, M., Zanna, C., Rugolo, M., Cossarizza, A., Wissinger, B., Verny, C., Schwarzenbacher, R., Martn, M. n., Arenas, J., Ayuso, C., Garesse, R., Lenaers, G., Bonneau, D., and Carelli, V. OPA1 mutations induce mitochondrial DNA instability and optic atrophy plus phenotypes. Brain 131(2), 338–351, 12 (2007).
    
    

24. [24].Liskova, P., Ulmanov, O., Tesina, P., Melsova, H., Diblik, P., Hansikova, H., Tesaov, M., and Votruba, M. Novel OPA1 missense mutation in a family with optic atrophy and severe widespread neurological disorder. Acta ophthalmologica 91, 02 (2013).
    
    

25. [25].Carelli, V., Musumeci, O., Caporali, L., Zanna, C., Morgia, C., Del Dotto, V., Porcelli, A., Rugolo, M., Valentino, M., Iommarini, L., Maresca, A., Barboni, P., Carbonelli, M., Trombetta, C., Valente, E. M., Patergnani, S., Giorgi, C., Pinton, P., Rizzo, G., and Zeviani, M. Syndromic parkinsonism and dementia associated with OPA1 missense mutations. Annals of Neurology 78, 21–38, 03 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ana.24410&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25820230&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 

26. [26].Lynch, D., Loh, S., Harley, J., Noyce, A., Martins, L. M., Wood, N., Houlden, H., and Plun-Favreau, H. Nonsyndromic Parkinson disease in a family with autosomal dominant optic atrophy due to OPA1 mutations. Neurology Genetics 3, e188, 10 (2017).
    
    

27. [27].Stafa, K., Tsika, E., Moser, R., Musso, A., Glauser, L., Jones, A., Biskup, S., Xiong, Y., Bandopadhyay, R., Dawson, V. L., Dawson, T. M., and Moore, D. J. Functional interaction of Parkinson’s disease-associated LRRK2 with members of the dynamin GTPase superfamily. Human Molecular Genetics 23(8), 2055–2077 (Stafa2014LRRK2).
    
    

28. [28].Jonikas, M., Madill, M., Mathy, A., Zekoll, T., Zois, C., Wigfield, S., Kurzawa-Akanbi, M., Browne, C., Sims, D., Chinnery, P., Cowley, S., and Tofaris, G. Stem cell modeling of mitochondrial parkinsonism reveals key functions of OPA1. Annals of Neurology 83, 03 (2018).
    
    

29. [29].Sironi, L., Restelli, L. M., Tolnay, M., Neutzner, A., and Frank, S. Dysregulated interorganellar crosstalk of mitochondria in the pathogenesis of Parkinson’s disease. Cells 9(1), 233 (2020).
    
    

30. [30].Herms, J. and Dorostkar, M. Dendritic spine pathology in neurodegenerative diseases. Annual Review of Pathology: Mechanisms of Disease 11, 02 (2016).
    
    

31. [31].Bae, J. and Kim, S. H. Synapses in neurodegenerative diseases. BMB reports 50, 03 (2017).
    
    

32. [32].Gómez-Suaga, P., Pérez-Nievas, B., Glennon, L., Lau, D., Paillusson, S., Mórotz, G., Calì, T., Pizzo, P., Noble, W., and Miller, C. The VAPB-PTPIP51 endoplasmic reticulum-mitochondria tethering proteins are present in neuronal synapses and regulate synaptic activity. Acta Neuropathologica Communications 7, 35, 12 (2019).
    
    

33. [33].Nishimura, A., Mitne-Neto, M., Silva, H., Richieri-Costa, A., Middleton, S., Cascio, D., Kok, F., Oliveira, J., Gillingwater, T., Webb, J., Skehel, P., and Zatz, M. A mutation in the vesicle-trafficking protein VAPB causes late-onset spinal muscular atrophy and amyotrophic lateral sclerosis. American journal of human genetics 75, 822–31, 12 (2004).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/425287&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15372378&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000224303500007&link_type=ISI) 

34. [34].Aliaga, L., Lai, C., Yu, J., Chub, N., Shim, H., Sun, L., Xie, C., Yang, W.-J., Lin, X., O’Donovan, M., and Cai, H. Amyotrophic lateral sclerosis-related VAPB P56S mutation differentially affects the function and survival of corticospinal and spinal motor neurons. Human molecular genetics 22, 06 (2013).
    
    

35. [35].Kim, J.-Y., Jang, A., Reddy, R., Yoon, W., and Jankowsky, J. Neuronal overexpression of human VAPB slows motor impairment and neuromuscular denervation in a mouse model of ALS. Human Molecular Genetics 25, ddw294, 08 (2016).
    
    

36. [36].Paillusson, S., Gmez-Suaga, P., Stoica, R., Little, D., Gissen, P., Devine, M., Noble, W., Hanger, D., and Miller, C. α-synuclein binds to the ER-mitochondria tethering protein VAPB to disrupt Ca2+ homeostasis and mitochondrial ATP production. Acta Neuropathologica 134, 03 (2017).
    
    

37. [37].Burré, J., Sharma, M., Tsetsenis, T., Buchman, V., Etherton, M., and Südhof, T. Alpha-synuclein promotes SNARE-complex assembly in vivo and in vitro. Science 329, 1663–7, 09 (2010).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMjkvNTk5OS8xNjYzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMDQvMjAyMC4wNC4wNi4yMDA1NDc5MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

38. [38].Rosenbloom, K. R., Armstrong, J., Barber, G. P., Casper, J., Clawson, H., Diekhans, M., Dreszer, T. R., Fujita, P. A., Guruvadoo, L., Haeussler, M., Harte, R. A., Heitner, S., Hickey, G., Hinrichs, A. S., Hubley, R., Karolchik, D., Learned, K., Lee, B. T., Li, C. H., Miga, K. H., Nguyen, N., Paten, B., Raney, B. J., Smit, A. F. A., Speir, M. L., Zweig, A. S., Haussler, D., Kuhn, R. M., and Kent, W. J. The UCSC Genome Browser database: 2015 update. Nucleic Acids Res. 43(D1), D670–D681 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gku1177&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25428374&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 

39. [39].Langmead, B., Trapnell, C., Pop, M., and Salzberg, S. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biology 10(3), R25–10 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/gb-2009-10-3-r25&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19261174&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 

40. [40].Kozomara, A. and Griffiths-Jones, S. miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Research 42(D1), D68–D73 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkt1181&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24275495&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000331139800011&link_type=ISI) 

41. [41].Benjamini, Y. and Hochberg, Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society. Series B (Methodological) 57(1), 289–300 (1995).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/2346101&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:A1995QE4&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995QE45300017&link_type=ISI) 

42. [42].Kiran, A. and Baranov, P. V. DARNED: a DAtabase of RNa EDiting in humans. Bioinformatics 26(14), 1772–1776 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btq285&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20547637&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000279474400014&link_type=ISI) 

43. [43].Ramaswami, G. and Li, J. B. RADAR: a rigorously annotated database of A-to-I RNA editing. Nucleic Acids Research 42(D1), D109 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkt996&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24163250&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000331139800018&link_type=ISI) 

44. [44].Mizuguchi, Y., Mishima, T., Yokomuro, S., Arima, Y., Kawahigashi, Y., Shigehara, K., Kanda, T., Yoshida, H., Uchida, E., Tajiri, T., et al. Sequencing and bioinformatics-based analyses of the microRNA transcriptome in Hepatitis B–related hepatocellular carcinoma. PloS ONE 6(1), e15304 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0015304&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21283620&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 

45. [45].Peng, Z., Cheng, Y., Tan, B. C., Kang, L., Tian, Z., Zhu, Y., Zhang, W., Liang, Y., Hu, X., Tan, X., Guo, J., Dong, Z., Liang, Y., Bao, L., and Wang, J. Comprehensive analysis of RNA-seq data reveals extensive RNA editing in a human transcriptome. Nat Biotech 30(3), 253–260, March (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nbt.2122&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22327324&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 

46. [46].Gong, J., Wu, Y., Zhang, X., Liao, Y., Sibanda, V. L., Liu, W., and Guo, A.-Y. Comprehensive analysis of human small RNA sequencing data provides insights into expression profiles and miRNA editing. RNA Biol 11(11), 1375–1385 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/15476286.2014.996465&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25692236&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 

47. [47].Han, M. and Zheng, Y. Comprehensive analysis of single nucleotide polymorphisms in human microRNAs. PLoS ONE 8(11), e78028 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0078028&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24223755&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 

48. [48].Hafner, M., Landthaler, M., Burger, L., Khorshid, M., Hausser, J., Berninger, P., Rothballer, A., Ascano, M., Jungkamp, A.-C., Munschauer, M., et al. Transcriptome-wide identification of RNA-binding protein and microRNA target sites by PAR-CLIP. Cell 141(1), 129–141 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2010.03.009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20371350&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000276211100018&link_type=ISI) 

49. [49].Memczak, S., Jens, M., Elefsinioti, A., Torti, F., Krueger, J., Rybak, A., Maier, L., Mackowiak, S. D., Gregersen, L. H., Munschauer, M., et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature 495(7441), 333–338 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature11928&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23446348&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 

50. [50].Fujita, P. A., Rhead, B., Zweig, A. S., Hinrichs, A. S., Karolchik, D., Cline, M. S., Goldman, M., Barber, G. P., Clawson, H., Coelho, A., et al. The UCSC Genome Browser database: update 2011. Nucleic Acids Research 39(suppl 1), D876–D882 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkq963&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20959295&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000285831700137&link_type=ISI) 

51. [51].Xie, C., Mao, X., Huang, J., Ding, Y., Wu, J., Dong, S., Kong, L., Gao, G., Li, C.-Y., and Wei, L. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic acids research 39(suppl 2), W316–W322 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkr483&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21715386&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000292325300051&link_type=ISI) 

52. [52].I Love, M., Huber, W., and Anders, S. Moderated estimation of fold change and dispersion for RNA-Seq data with DESeq2. Genome biology 15, 550, 12 (2014).
    
    

53. [53].Zhang, Y., James, M., Middleton, F. A., and Davis, R. L. Transcriptional analysis of multiple brain regions in parkinson’s disease supports the involvement of specific protein processing, energy metabolism, and signaling pathways, and suggests novel disease mechanisms. American Journal of Medical Genetics Part B Neuropsychiatric Genetics 137b(2005).
    
    

54. [54].Larkin, M. A., Blackshields, G., Brown, N., Chenna, R., McGettigan, P. A., McWilliam, H., Valentin, F., Wallace, I. M., Wilm, A., Lopez, R., et al. Clustal W and Clustal X version 2.0. Bioinformatics 23(21), 2947–2948 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btm404&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17846036&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000251197700021&link_type=ISI) 

55. [55].Page, R. D. TreeView: an application to display phylogenetic trees on personal computers. Computer applications in the biosciences: CABIOS 12(4), 357–358, August (1996).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8902363&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.04.06.20054791.atom)