Next Article in Journal
Identification of Novel Core Genes Involved in Malignant Transformation of Inflamed Colon Tissue Using a Computational Biology Approach and Verification in Murine Models
Previous Article in Journal
Lacrimal Gland Epithelial Cells Shape Immune Responses through the Modulation of Inflammasomes and Lipid Metabolism
Previous Article in Special Issue
MicroRNAs Expression in Response to rhNGF in Epithelial Corneal Cells: Focus on Neurotrophin Signaling Pathway
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

circRNA from APP Gene Changes in Alzheimer’s Disease Human Brain

by
Amaya Urdánoz-Casado
1,
Javier Sánchez-Ruiz de Gordoa
1,2,
Maitane Robles
1,
Miren Roldan
1,
Mónica Macías Conde
1,2,
Blanca Acha
1,
Idoia Blanco-Luquin
1 and
Maite Mendioroz
1,2,*
1
Neuroepigenetics Laboratory-Navarrabiomed, Complejo Hospitalario de Navarra-IdiSNA (Navarra Institute for Health Research), Universidad Pública de Navarra (UPNA), Pamplona, 31008 Navarra, Spain
2
Department of Neurology, Complejo Hospitalario de Navarra-IdiSNA (Navarra Institute for Health Research), Pamplona, 31008 Navarra, Spain
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(5), 4308; https://doi.org/10.3390/ijms24054308
Submission received: 22 December 2022 / Revised: 9 February 2023 / Accepted: 18 February 2023 / Published: 21 February 2023
(This article belongs to the Special Issue RNA Regulatory Networks at the Crossroad of Human Diseases 2.0)

Abstract

:
Alzheimer’s disease (AD) is the most common cause of age-related dementia. Amyloid precursor protein (APP) is the precursor of Aβ peptides, and its role in AD has been widely investigated. Recently, it has been reported that a circular RNA (circRNA) originated from APP gene can serve as a template for Aβ synthesis, postulating it as an alternative pathway for the Aβ biogenesis. Moreover, circRNAs play important roles in brain development and in neurological diseases. Therefore, our aim was to study the expression of a circAPP (hsa_circ_0007556) and its linear cognate in AD human entorhinal cortex, a brain region most vulnerable to AD pathology. First, we confirmed the presence of circAPP (hsa_circ_0007556) in human entorhinal cortex samples using RT-PCR and Sanger sequencing of PCR products. Next, a 0.49-fold decrease in circAPP (hsa_circ_0007556) levels was observed in entorhinal cortex of AD cases compared to controls (p-value < 0.05) by qPCR. In contrast, APP mRNA expression did not show changes in the entorhinal cortex between AD cases and controls (Fold-change = 1.06; p-value = 0.81). A negative correlation was found between Aβ deposits and circAPP (hsa_circ_0007556) and APP expression levels (Rho Spearman = −0.56, p-value < 0.001 and Rho Spearman = −0.44, p-values < 0.001, respectively). Finally, by using bioinformatics tools, 17 miRNAs were predicted to bind circAPP (hsa_circ_0007556), and the functional analysis predicted that they were involved in some pathways, such as the Wnt-signaling pathway (p = 3.32 × 10−6). Long-term potentiation (p = 2.86 × 10−5), among others, is known to be altered in AD. To sum up, we show that circAPP (hsa_circ_0007556) is deregulated in the entorhinal cortex of AD patients. These results add to the notion that circAPP (hsa_circ_0007556) could be playing a role in the pathogenesis of AD disease.

1. Introduction

Alzheimer’s disease (AD) is a chronic and irreversible neurodegenerative disease [1]. AD is the leading cause of age-related dementia and is also the most common neurodegenerative disease [2,3]. The main anatomopathological features of AD are the brain deposition of intraneuronal neurofibrillary tangles (NFTs) of hyperphosphorylated tau protein and extracellular plaques of Amyloid beta peptide (Aβ) in the parenchyma and blood vessels.
AD can be classified into two types, familial AD and sporadic AD. The former constitutes 5% of AD cases and is characterized by mutations in the Amyloid Precursor Protein (APP), Presenilin 1 (PSEN1) or Presenilin 2 (PSEN2) genes, all of which are related to the synthesis and processing of Aβ-peptide. Regarding sporadic AD, the cause of the disease remains unknown, although it is considered a multifactorial disease in which genetic and environmental risk factors contribute to its development [2].
APP is a membrane protein which performs several crucial cellular functions. So far, a number of APP isoforms have been described by alternative splicing; one of which is commonly expressed in the brain, where it participates in synaptogenesis and synaptic plasticity [4,5]. APP can be processed through the amyloidogenic pathway, in which the final products are going to be the Aβ40 and Aβ42 peptides, components of amyloid plaque, and the APP intracellular domain (amyloid precursor protein intracellular domain, AICD), following proteolytic cleavage of the APP protein carried out by the enzymes β-secretase (BACE1) and γ-secretase. In the non-amyloidogenic pathway, the enzymes in charge of processing the APP protein are α-secretase (ADAM metallopeptidase domain 10, ADAM10) and γ-secretase, and the final product will be irrelevant peptides for this pathology [1,6].
Although the role of Aβ in the AD pathogenesis is yet not well understood [7], and some authors believe that aberrant Aβ expression may not be the primary cause of all EOAD [8], the amyloid cascade theory still remains the most widely accepted pathogenic model [6]. This theory proposes that Aβ deposition is the first critical event that triggers a cascade of molecular phenomena leading to neurofibrillary deposits, synaptic failure, neuronal death, neurodegeneration and, finally, AD dementia. It seems clear that deposition of Aβ occurs as a consequence of altered processing of APP or Aβ clearance. However, it is still not well understood how the peptide is produced [9].
Recently, it has been described that alternative splicing of APP results in the generation of several circRNAs. Mo et al. even demonstrated that Aβ can be transcribed from a specific circRNA (hsa_circ_0007556) [9]. circRNAs are single-stranded RNA molecules characterized, as their name suggests, by a circular structure. They lack 5′ and 3′ ends since, after transcription, a covalent bond is established between these ends [10,11]. circRNAs can act as miRNA sponges, protein templates or transcriptional regulators, among other described functions [12,13]. circRNAs are evolutionarily conserved and expressed in a large number of body tissues [10,11,14,15]. However, it is in the brain that they are most highly expressed; in fact, 20% of brain genes encode circRNAs [16,17]. Interestingly, neuronal specialization has been related to a high level of alternative splicing [18], and this process could also explain the enormous amount of brain-specific circRNAs [14]. Nevertheless, the expression levels of circRNAs vary from one brain region to another, with the synapse being the site where their expression is highest [19]. Thus, alterations in the expression levels of circRNAs have been described in different neurological diseases, such as multiple sclerosis (MS), Parkinson’s disease (PD), amyotrophic lateral sclerosis (ALS) and AD, among others [20,21,22,23,24].
In the last few years, thanks to advances in massive sequencing and transcriptome analysis, knowledge about circRNAs involvement in AD has increased exponentially. Numerous studies, both in human brain or blood samples and in cellular and mouse models of AD, have shown dysregulation of circRNAs in this disease. A number of circRNAs are now known to play important regulatory functions in neuroinflammation (e.g., circ_0000950), oxidative stress (e.g., mmu_circRNA_013636 and mmu_circRNA_012180) and autophagy (e.g., circNF1-419), as well as Aβ production and degradation (e.g., circHDAC9, CDR1as, circHOMER1 or circCORO1C) [25,26,27,28,29,30,31,32,33]. In addition, significant changes in circRNAs expression have been detected whose genes of origin are closely related to AD pathology, involved in synaptic plasticity and neuronal survival (e.g., HOMER1, DOCK1, NTRK2 and APC) or vesicular trafficking (e.g., DGL1/SAP97, TRAPPC9 and KIF1B) [34]. From the APP and the MAPT genes, which encode the main proteins involved in AD [9,35], or from the APOE gene [36], circRNAs also originate. Likewise, several possible candidate circRNAs for diagnosis, prognosis and disease progression have been postulated, such as circHOMER1, circCOROC1 or hsa_circ_0003391 [25,26,37,38,39,40]. Of all these, circRNAs derived from the APP gene are particularly interesting. Thirty-three circRNAs derived from the APP gene (circAPPs) have been described by in silico analysis of RNA sequencing data, and 17 circAPP have been observed by RT-PCR and Sanger sequencing. Interestingly, only two of them, namely hsa_circ_0007556 and hsa_circ_0115725, were in both lists, with hsa_circ_0007556 being the most abundant one [9,14,41,42]. It is known that this circAPP (hsa_circ_0007556) can be detected in two brain regions, hippocampus and prefrontal lobe, but its expression in the entorhinal cortex is unknown. Therefore, our aim was to study the expression of circAPP (hsa_circ_0007556) and its related linear form in human entorhinal cortex, given that knowledge of the circAPP (hsa_circ_0007556) expression pattern of in both healthy and AD-affected human entorhinal cortex is scarce, and this region constitutes one of the brain areas most vulnerable to the development of AD. Neuronal loss at entorhinal cortex occurs very early in AD [43], and dysfunction in this region is implicated in memory and learning impairment [44]. In addition, entorhinal cortex is the main communication pathway between the hippocampus and neocortex [44]. Furthermore, gene expression profiling changes related with memory and learning functions have been identified in entorhinal cortex from AD and control post-mortem brain samples [45].

2. Results

2.1. circAPP (hsa_circ_0007556) Identification in Entorhinal Human Brain

First, we wanted to test whether circAPP (hsa_circ_0007556) was expressed in the human entorhinal cortex. For this purpose, APP-specific divergent primers were designed to amplify circular but not linear RNA. Then, RT-PCR was performed on RNA samples isolated from the entorhinal cortex of AD patients and controls. After electrophoresis, PCR products were selected from the agarose gel for Sanger sequencing analysis (Figure 1). As a result, we succeeded to amplify circAPP (hsa_circ_0007556) in the human entorhinal cortex in both AD patients and controls.

2.2. Differential Expression of circAPP (hsa_circ_0007556) and APP mRNA in AD Entorhinal Cortex

In order to study the expression levels of circAPP (hsa_circ_0007556) detected in the previous section and to explore whether there were expression differences in entorhinal cortex samples between AD and controls, RT-qPCR technique was performed. A total of 29 entorhinal cortex samples from AD patients and 16 controls were studied. All samples passed the RNA quality criteria.
It should be noted that statistically significant differences in age (mean SD, 55.94 ± 5.31 in controls versus 82.07 ± 1.96 in AD, p-value < 0.001) and gender (% female, 31.25% in controls versus 62.07% in AD, p-value < 0.05) were found between AD samples and controls, so we fitted multivariate linear regression models for each transcript (circAPP (hsa_circ_0007556) and its corresponding linear mRNA), including the presence or absence of the disease as predictor and adjusting for age and gender as potential confounders, and it was found that only the presence of the disease significantly explained differences in the expression of circAPP transcript (p = 0.025). Regression coefficients and standard errors can be found in Table S1. Additionally, we checked whether circAPP and APP mRNA expression is decreased with age (Figure S1) or by gender. We observed no correlation between circAPP expression and age (p-value = 0.306) or gender (p-value = 0.655), nor between APP mRNA expression and age (p-value = 0.158) or gender (p-value = 0.063).
A significant decrease (fold-change (FC) = 0.49, p-value < 0.05) in the circAPP (hsa_circ_0007556) levels was observed in the entorhinal cortex of AD cases compared to controls. However, APP mRNA expression (FC = 0.80, p-value = 0.81) showed no significant changes in the entorhinal cortex region of AD samples compared to controls (Figure 2A). To evaluate whether the magnitude of the differences in the transcript expression levels between AD cases and controls was gender-dependent, we included in the multivariate linear regression models the interaction between the presence/absence of the disease and gender, resulting in non-statistically significant interactions (circAPP (hsa_circ_0007556), p-value = 0.544 and for APP mRNA, p-value = 0.549). These results demonstrated the absence of sex dependence in the differential expression levels of the transcripts between AD cases and controls.
We also wanted to analyze whether circAPP (hsa_circ_0007556) and APP mRNA expression changed across AD neuropathological stages, according to the ABC score. A significant decrease in circAPP (hsa_circ_0007556) expression was observed in the group with higher ABC score values with respect to controls (p-value < 0.05) (Figure 2B). For its part, APP linear transcript showed a significant decrease in expression in the group with high level of ABC score with respect to the intermediate group (p-value < 0.05) and the low group (p-value < 0.001) (Figure 2B).
Next, using a univariate general linear model, where the dependent variable was the RNA expression levels and the independent variables were the specific RNA variants (circAPP (hsa_circ_0007556) and APP mRNA) and the presence of disease (AD or control), it was observed that the expression of the two transcripts together in AD cases was 37.75% lower (p-values < 0.05) with respect to the controls. Considering linear APP as a reference, circAPP (hsa_circ_0007556) expression was 99.93% lower (p-value < 0.0001). However, when the cases (AD and control) and the two RNA variants were studied together in the model, relative expression ratios of the two transcripts (circAPP/APP) were maintained in the AD cases and controls (99.94% and 99.90% lower circAPP (hsa_circ_0007556) relative to APP mRNA, respectively, p-value = 0.100), meaning that the decrease in expression of the 2 RNA variants is proportional within AD cases, although only circAPP (hsa_circ_0007556) was significantly decreased in samples with AD compared to controls (Figure S2).
We repeated the above analysis, but taking into consideration the ABC score instead of the presence of the disease (AD or control). Relative expression ratios of the two transcripts (circAPP/APP) were also maintained in each level of ABC scores and controls (99.95% at the low level, 99.93% at the intermediate level, 99.93% at the high level, 99.90% at the control level and lower circAPP relative to APP mRNA, p-value = 0.455, p-value = 0.955, p-value = 0.362 between groups) (Figure S3).
We can conclude that the decrease in expression of the different variants is proportional within AD or ABC score levels and control samples; in fact, the circAPP/APP ratio is maintained between AD patients or ABC score levels and controls.

2.3. Correlation between circAPP (hsa_circ_0007556) and APP mRNA Expression and Aβ Deposits

Since one of the main pathophysiological features of AD is the deposition of Aβ peptide, which derives from APP protein processing, we sought to study the relationship between quantitative assessment of Aβ deposits and APP-derived RNA transcripts expression levels. Thus, a negative correlation was found between Aβ deposits and both circAPP (hsa_circ_0007556) and linear APP expression (Rho Spearman = −0.48, p-value < 0.05 and Rho Spearman = −0.652, p-value <0.01, respectively) (Table 1, Figure S4).
Additionally, we decided to study the correlation between circAPP (hsa_circ_007556) or APP mRNA expression and the ABC score and each group individually forming ABC score (Method of Thal, Braak and Braak classification, Method of CERAD). We observed a negative correlation between circAPP (hsa_circ_007556) expression and all pathological features analyzed (Method of Thal, Braak and Braak classification, Method of CERAD, ABC score and Global average area of Aβ deposits), while APP mRNA showed a negative correlation with all of them, except the Method of CERAD) (Table 1).

2.4. In Silico Prediction of Biological Function

It has been demonstrated that circAPP (hsa_circ_0007556) could be the template for Aβ peptide transcription, but little is known about the role of circAPP (hsa_circ_0007556) in the brain [9]. Since circRNAs may function as microRNAs sponges, we decided to identify those miRNAs that target circAPP (hsa_circ_0007556). For this purpose, we used the miRNA target sites tool from CircInteractome database [46]. We identified 17 miRNAs that could potentially target circAPP (hsa_circ_0007556) (Table 2). Some of these miRNAs have been described in different neurological disorders, such as major depression, epilepsy, ALS, PD or AD [47,48,49,50,51,52,53,54]. Moreover, we wanted to know the biological pathways in which these miRNAs are involve and tried to approach the biological function of circAPP (hsa_circ_0007556). We employed the microT-CDS tool of the DIANA mirPath v.3 software [55] and used the 17 miRNAs as the input. A significant association between these miRNAs and diverse KEGG pathways were found, as follows: Wnt-signaling pathway (p = 3.32 × 10−6), Long-term potentiation (p = 2.86 × 10−5), Glycosaminoglycan biosynthesis—heparan sulfate/heparin (hsa00534) (p = 3.38 × 10−5), Ubiquitin mediated proteolysis (p = 6.68 × 10−5), Axon guidance (7.01 × 10−5) and Glutamatergic synapse (p = 7.76 × 10−5), among others (Table S2).

3. Discussion

In this study, we observed downregulation of circAPP (hsa_circ_0007556) expression in the entorhinal cortex of AD patients versus controls. Furthermore, taking into account the progression of AD neuropathological change and according to the ABC score, circAPP (hsa_circ_0007556) expression was downregulated in the group with higher ABC score values with respect to controls and APP mRNA showed downregulation in the group with higher ABC score values with respect to intermediate and low groups. In addition, both transcripts showed a negative correlation with Aβ deposits in the brain tissue.
The APP protein is the precursor of Aβ peptides and its role in AD has been extensively investigated [73]. Databases of circRNAs predict up to 33 different circRNAs derived from the APP gene. Some of these circRNAs were detected by massive sequencing in different brain regions such as cerebellum, frontal cortex, diencephalon, occipital lobe, parietal lobe and temporal lobe [14,41,42]. In recent years, the focus has been on studying the expression of circRNAs in different brain regions affected by AD. For example, Lo et al. [40] detected by RNA-seq, in four brain regions, i.e., the anterior prefrontal cortex, superior temporal lobe, inferior frontal lobe and hippocampus, five circRNAs originating from APP, but only one was differentially expressed in the inferior frontal lobe between AD patients and controls. However, they did not observe the circAPP (hsa_circ_0007556) studied in the present work. On the other hand, Dube et al. [25] surveyed the parietal cortex, inferior frontal lobe, frontal pole, superior temporal lobe and parahippocampal lobe also by RNA-seq but found none of the APP-derived circRNAs. That could be because they may not meet the minimum requirements for the subsequent differential analysis.
It is worth noting that, of all the circAPPs revealed by massive sequencing, only hsa_circ_0007556 and hsa_circ_0115725 has been validated by another technique in two brain regions, the hippocampus and prefrontal lobe [9]. Following the results obtained in our study, the expression of hsa_circ_0007556 in the entorhinal cortex of both AD patients and controls can also be confirmed.
Moreover, with this work, knowledge of the differential expression of circAPP (hsa_circ_0007556), which was found to be less expressed in samples with AD compared to controls, is added. A priori, an increase in its expression would have been expected, since the mechanism by which Aβ deposition is increased in AD brain is not entirely clear, and an increase in this circRNA could imply an increase in the production of Aβ [9]. In fact, it has been recently published that this circAPP (hsa_circ_0007556) may serve as a “template” for Aβ synthesis, postulating itself as an alternative pathway for this peptide biogenesis [9]. On the other hand, a negative correlation of peptide Aβ burden and circAPP (hsa_circ_0007556) expression levels is observed here (although when was taking into consideration the ABC score, only the group with higher ABC score values showed a significant decrease in circAPP (hsa_circ_0007556) expression respect to control, but this observation could be due to the reduced sample size by segmenting into ABC stages limits the power of the analysis), leading to hypothesize that it could play a regulatory role, directly or indirectly, on the expression of enzymes in charge of APP protein processing. For example, due to the upregulation of the peptide, a downregulation of circRNA could be induced to compensate for the excess of Aβ that could be produced from it. In fact, the expression of this circRNA seems to be reduced across the progression of the disease. An alternative explanation would be that enzymes involved in Aβ peptide processing may bind to circAPP (hsa_circ_0007556), so that circAPP (hsa_circ_0007556) would control the availability of them; a circAPP (hsa_circ_0007556) dysfunction would lead to the release of the enzymes, resulting in an increase of Aβ. In any case, our study is observational and merely shows a statistical association, not causality. Therefore, to elucidate whether there is a true involvement of circAPP (hsa_circ_0007556) in the generation of the peptide, further experimental studies would be necessary.
In any case, it would be interesting to study several aspects of this circRNA, such as the expression of the protein that originates from circAPP (hsa_circ_0007556), the methylation of circAPP (hsa_circ_0007556), since it has been demonstrated that m6A RNA methylation can efficiently promote translation to proteins [27] or the interaction of this circRNA with other molecules to try to elucidate its function in the brain and its role in AD. It would also be interesting to study the localization of this circRNA at the cellular level, as it is known that the expression of circRNAs is enriched at the synapse [16], and synaptic dysfunction is one of the characteristics of AD [7].
Anyway, it must be taken into consideration that the decrease in the level of circAPP may be due to a number of reasons, such as downregulation of its expression, higher processing or destruction or sequestration of circRNA in the amyloid plaques.
Although circAPP (hsa_circ_0007556) expression is downregulated, the expression of linear APP is not altered on the whole. These findings are in the same direction as others previously shown in the literature [9]. Nevertheless, one would expect that the ratio of circAPP (hsa_circ_0007556) to linear APP between patients and controls may be altered, given that circAPP (hsa_circ_0007556) shows expression changes in the face of AD while APP does not. However, what is observed is that the circAPP (hsa_circ_0007556)/APP ratio is maintained between AD patients and controls. This may be due to a decrease in APP expression, although not statistically significant, in AD patients and this small decrease may be sufficient to maintain the circAPP (hsa_circ_0007556)/APP ratio.
Both RNA transcripts show a negative correlation with Aβ deposits, and this points to other mechanisms being involved in this process, such as those concerning miRNAs. In the CircInteractome database [46], circAPP (hsa_circ_0007556) is predicted to have a binding site for 17 miRNAs. Some of them have been associated with several diseases, but it is worth mentioning that at least five of them have been related with AD. For example, hsa-miR-598, hsa-miRNA-659 and hsa-miR-324-5p have been identified as candidate biomarkers in different fluids such as plasma or cerebrospinal fluid [47,48,49]. In an AD model in primary mouse hippocampal neurons, circ_0004381 was found to regulate PSEN1 expression through miR-647 [67]. has-miR-186 is a strong negative regulator of BACE1 expression, a protein involved in the amyloidogenic pathway [61]. Furthermore, in the review by He et al. [74], it is pointed out that different miRNAs participate in the metabolism of Aβ peptides acting at different levels and on different genes.
Moreover, these 17 miRNAs found in the in silico functional analysis have been significantly associated with multiple pathways, some of which are altered in AD. For instance, Wnt-signaling pathway is implicated in neurodevelopment and neurogenesis and, in AD, is downregulated in several cell types in AD brains [75]. Other relevant pathways, such as synaptic long-term potentiation, axon guidance and dysregulation in ubiquitin-mediated proteolysis, which are implicated in synaptic plasticity and synaptic development, have been related to cognitive impairment observed in AD [76,77,78]. Glycosaminoglycan biosynthesis—heparan sulfate/heparin is associated with the formation of Aβ plaques [79,80,81], and the glutamatergic synapse pathway has been observed altered in AD, while glutamatergic receptors have been suggested as pharmacological targets [82]. However, this is still an in silico study, and other functional studies are required to confirm these predictions.

4. Methods

4.1. Human Entorhinal Cortex Samples

Brain entorhinal cortex samples from 29 AD patients and 16 controls were provided by Navarrabiomed Brain Bank. After death, half-brain specimens from donors were cryopreserved at −80 °C. A neuropathological examination was completed following the usual recommendations [83] and according to the updated National Institute on Aging-Alzheimer’s Association guidelines [84].
Assessment of Aβ deposition was carried out by immunohistochemical staining of paraffin-embedded sections (3–5 μm thick) with a mouse monoclonal (S6 F/3D) anti-Aβ antibody (Leica Biosystems Newcastle Ltd., Newcastle upon Tyne, United Kingdom). Evaluation of neurofibrillary pathology was performed with a mouse monoclonal antibody anti-human PHF-TAU, clone AT-8 (Tau AT8) (Innogenetics, Gent, Belgium), which identifies hyperphosphorylated tau (p-tau) [85]. The reaction product was visualized using an automated slide immunostainer (Leica Bond Max) with Bond Polymer Refine Detection (Leica Biosystems, Newcastle Ltd., UK). Other protein deposits, such as synuclein deposits, were ruled out by a monoclonal antibody against α-synuclein (NCL-L-ASYN; Leica Biosystems, Wetzlar, Germany). The staging of AD was performed by using the ABC score according to the updated National Institute on Aging-Alzheimer’s Association guidelines [84]. ABC score combines histopathologic assessments of Aβ deposits determined by the method of Thal (A) [84], staging of neurofibrillary tangles by Braak and Braak classification (B) [85], and scoring of neuritic plaques by the method of CERAD (Consortium to Establish A Registry for Alzheimer’s Disease) (C) [86] to characterize AD neuropathological changes. Thus, the ABC score shows three levels of AD neuropathological severity: low, intermediate and high. A summary of the characteristics of subjects considered in this study is shown in Table S3.

4.2. RNA Isolation and Reverse Transcription—Polymerase Chain Reaction (RT-PCR)

Total RNA, including small RNA species, was isolated from cells and entorhinal cortex samples with miRNAeasy mini Kit (QIAGEN, Redwood City, CA, USA) following the manufacturer’s instructions. Concentration and purity of RNA were both evaluated with NanoDrop spectrophotometer. Complementary DNA (cDNA) was reverse transcribed from 500 ng total RNA with SuperScript® III First-Strand Synthesis Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA) after priming with random primers. RT-PCR was performed by using GoTaq® DNA polymerase (Promega, Madison, WI, USA) in an Applied Biosystems™ Veriti™ Thermal Cycler, 96-Well (Applied Biosystems, Foster City, CA, USA). PCR conditions were as follows: denaturation at 95 °C for 20 s, extension at 72 °C for 30 s and annealing temperatures 60 °C for 40s and cycles used were 40. Primer3 software was used for divergent primers design (Table S4).

4.3. Candidate Band Selection and Sanger Sequencing

Candidate bands were selected after 1.8% agarose gel electrophoresis of RT-PCR products. Bands purification were made with Wizard® SV Gel and PCR Clean-Up System (Promega, Madison, WI, USA). Next, Sanger sequencing was performed and UCSC (University of California Santa Cruz) Genome Browser software was used for the sequence alignment [87,88].

4.4. Real Time Quantitative PCR (RT-qPCR) Assay

Total RNA were isolated from the entorhinal cortex with RNAeasy Lipid Tissue mini Kit (QIAGEN, Redwood City, CA, USA) following the manufacturer’s instructions. Genomic DNA was removed with recombinant DNase (TURBO DNA-free™ Kit, Ambion, Austin, TX, USA). Concentration and purity of RNA were both evaluated with NanoDrop spectrophotometer. Complementary DNA (cDNA) was reverse transcribed from 500 ng total RNA with SuperScript® III First-Strand Synthesis Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA) after priming with random primers. RT-qPCR reactions were performed in triplicate with Power SYBR Green PCR Master Mix (Invitrogen, Carlsbad, CA, USA) on the QuantStudio 12K Flex real-time PCR system (Applied Biosystems, Foster City, CA, USA) and repeated twice on independent cDNA samples. The sequences of convergent primer pairs for linear RNA detection were designed using the IDT real-time PCR tool (Coralville, IA, USA) and Primer3 software and are listed in Table S4. The relative expression level of mRNA in each sample was calculated using the delta delta-CT method and the geometric mean of GAPDH and ACTB genes was used as a reference to normalize the expression values [89].

4.5. Quantitative Assessment of Aβ Deposits in Brain Tissues

In order to quantitatively assess the Aβ burden for further statistical analysis, we applied a method to quantify protein deposits. This method generates a numeric measurement that represents the extent of Aβ deposition. Sections of the entorhinal cortex were examined after performing immunostaining with anti Aβ antibody as described above in Human Entorhinal Samples. Three pictures were obtained for each immunostained section by using an Olympus BX51 microscope at ×10 magnification power. Focal deposit of Aβ, as described by Braak & Braak (neuritic, immature, and compact plaque) [85], was manually determined and was further edited and analyzed with ImageJ software. Then, the Aβ plaque count, referred to as amyloid plaque score (APS) and total area of Aβ deposition, was automatically measured by ImageJ and averaged for each section (Figure S5).

4.6. Prediction of circAPP (hsa_circ_0007556) Interaction with miRNAs and Biological Function

With the help of public databases containing information about ncRNAs we can predict which genes, proteins and miRNAs candidate circRNAs interact with. For the prediction of miRNAs that could bind to circAPP (hsa_circ_0007556), the miRNA target sites tool from the CircInteractome database was used [46]. For the study of the potential biological pathways that miRNAs might be regulating, microT-CDS tool of the DIANA mirPath v.3 software [55] was employed, and as for input, we included the miRNAs obtained in the previous step.

4.7. Statistical Analysis

Statistical analysis was performed with SPSS 25.0 (IBM, Inc., Chicago, IL, USA). Before performing differential analysis, we checked whether continuous variables follow a normal distribution, as per one-sample Kolgomorov–Smirnov test and the normal quantil–quantil (Q–Q) plots. For the analysis of the differential expression of the distinct APP transcripts, we fitted multivariate linear regression models for each transcript (circAPP and its corresponding linear mRNA), including the presence or absence of the disease as predictor and adjusting for age and gender as potential confounders. All three explanatory variables were included using the Enter method (all variables are entered in a single step). To evaluate the homogeneity of variances, the Levene’s test was used, and the normality of the regression residuals was assessed by visualization of the histograms and Q–Q plots. All models met the aforementioned requirements. In order to analyze differences in the expression levels of the different transcripts studied between the ABC scale groups, a general univariate linear model adjusted for gender and age and the Bonferroni post hoc test were developed. On the other hand, a general linear univariate model was used to determine the proportions between the expression levels of the APP RNA variants in the AD samples versus the control samples and ABC score stages. For each gene, an expression variable was created where the log(expression) of each variant and another categorical variable with the type of transcript was collected, leaving a model where the expression variable was the dependent variable and as fixed factors the case variables ((AD or control) or (ABC score stages (control, low, intermediate and high) and the type of transcript were included. Spearman’s test was used to assess the correlation between the continuous variables circRNA or mRNA expression and Aβ deposition or Method of Thal or Braak and Braak classification or Method of CERAD or ABC score. GraphPad Prism version 9.00 for Windows (GraphPad Software, La Jolla, CA, USA) was used to draw graphs.

5. Conclusions

We observed the expression of circAPP (hsa_circ_0007556) in the human entorhinal cortex, and we also show that circAPP (hsa_circ_0007556) is downregulated in the entorhinal cortex of AD patients compared to controls and at late stages of neuropathological changes. These results add to the notion that circAPP (hsa_circ_0007556) could be playing a role in the pathogenesis of AD disease.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms24054308/s1: Table S1. Regression coefficients and standard errors of each multivariate linear regression model of each transcript circAPP (hsa_circ_0007556) and APP mRNA. Table S2. List of KEGG pathways significantly associated with miRNAs predicted to bind circAPP (hsa_circ_0007556). Table S3. Characteristics of the samples included in the study. Table S4. It shows the sequences of the primers used in the different PCRs carried out in this work. Figure S1. APP RNA variants and Age correlation. Figure S2. Representation of the proportional expression between controls and AD cases of the log (relative expression) of the two APP variants. Figure S3: Representation of the proportional expression between controls and AD neuropathological stages, according to the ABC score of the log(relative expression) of the two APP variants. Figure S4. APP RNA variants and Aβ deposits correlation. Figure S5. Quantitative assessment of Aβ deposits.

Author Contributions

A.U.-C. contributed to the study concept and design, running of the experiments, acquisition of the data, analysis and interpretation of the data, statistical analysis, figure drawing and drafting/revising of the manuscript for content. J.S.-R.d.G. contributed to the analysis and interpretation of the data (amyloid deposits), sorting of the patients into different stages and acquisition of the image data. M.R. (Maitane Robles) and M.R. (Miren Roldan) contributed to the running of the experiments and acquisition of the data. M.M.C. and B.A. contributed to figure drawing and drafting/revising of the manuscript for the content. I.B.-L. contributed to the study concept and design and drafting/revising of the manuscript for the content. M.M. contributed to the drafting/revising of the manuscript for content, study concept and design, analysis and interpretation of the data, study supervision and obtaining the funding. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Spanish Government through grants from the Institute of Health Carlos III (FIS PI20/01701). In addition, A.U.-C. received two grants “Doctorados industriales 2018–2020” and “Contrato predoctoral en investigación en ciencias y tecnologías de la salud en el periodo 2019–2022”, both of them founded by the Government of Navarra. In addition, BA was supported by a PFIS fellowship from the Spanish Government (FI18/00150), MMC was beneficiary of a grant “Río Hortega” from the Spanish Government (CM20/00240) and M.M. received a grant “Programa de intensificación” (LCF/PR/PR15/51100006) funded by Fundación Bancaria “la Caixa” and Fundación Caja-Navarra and “Contrato de intensificación” from the Institute of Health Carlos III (INT19/00029).

Institutional Review Board Statement

The Navarrabiomed Brain Bank board reviewed and approved this study. This study was carried out in accordance with the principles of the Declaration of Helsinki, and handling of the human brain samples was performed according to the current Spanish national legislation (Royal Decree RD1716/2011). The study was approved by the Ethics Committee of the Complejo Hospitalario de Navarra (90/2014) and the Scientific Committee of Navarrabiomed Brain Bank, Spain.

Informed Consent Statement

Written informed consent was obtained by subjects or relatives before donation of brain tissue was performed.

Data Availability Statement

Not applicable.

Acknowledgments

We want to kindly thank Mónica Enguita Germán (Methodology Department of Navarrabiomed, technical support) and Isabel Gil-Aldea (Navarrabiomed Biobank, technical support) for their help. Finally, we are very grateful to the patients and relatives that generously donated the brain tissue to the Navarrabiomed Brain Bank.

Conflicts of Interest

The authors declare that they have no competing interests.

References

  1. Tortajada-Soler, M.; Sánchez-Valdeón, L.; Blanco-Nistal, M.; Benítez-Andrades, J.A.; Liébana-Presa, C.; Bayón-Darkistade, E. Prevalence of Comorbidities in Individuals Diagnosed and Undiagnosed with Alzheimer’s Disease in León, Spain and a Proposal for Contingency Procedures to Follow in the Case of Emergencies Involving People with Alzheimer’s Disease. Int. J. Environ. Res. Public Health 2020, 17, 3398. [Google Scholar] [CrossRef] [PubMed]
  2. Reitz, C.; Mayeux, R. Alzheimer disease: Epidemiology, diagnostic criteria, risk factors and biomarkers. Biochem. Pharmacol. 2014, 88, 640–651. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Alzheimer’s Association. 2021 Alzheimer’s disease facts and figures. Alzheimers Dement. 2021, 17, 327–406. [Google Scholar] [CrossRef]
  4. O’Brien, R.J.; Wong, P.C. Amyloid precursor protein processing and Alzheimer’s disease. Annu. Rev. Neurosci. 2011, 34, 185–204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Gralle, M.; Ferreira, S.T. Structure and functions of the human amyloid precursor protein: The whole is more than the sum of its parts. Prog. Neurobiol. 2007, 82, 11–32. [Google Scholar] [CrossRef] [PubMed]
  6. Niu, H.; Álvarez-Álvarez, I.; Guillén-Grima, F.; Aguinaga-Ontoso, I. Prevalence and incidence of Alzheimer’s disease in Europe: A meta-analysis. Neurologia 2017, 32, 523–532. [Google Scholar] [CrossRef]
  7. Morris, G.P.; Clark, I.A.; Vissel, B. Inconsistencies and controversies surrounding the amyloid hypothesis of Alzheimer’s disease. Acta Neuropathol. Commun. 2014, 2, 135. [Google Scholar] [CrossRef] [Green Version]
  8. Del-Aguila, J.L.; Benitez, B.A.; Li, Z.; Dube, U.; Mihindukulasuriya, K.A.; Budde, J.P.; Farias, F.H.G.; Fernandez, M.V.; Ibanez, L.; Jiang, S.; et al. TREM2 brain transcript-specific studies in AD and TREM2 mutation carriers. Mol. Neurodegener. 2019, 14, 18. [Google Scholar] [CrossRef]
  9. Mo, D.; Li, X.; Raabe, C.A.; Rozhdestvensky, T.S.; Skryabin, B.V.; Brosius, J. Circular RNA Encoded Amyloid Beta peptides-A Novel Putative Player in Alzheimer’s Disease. Cells 2020, 9, 2196. [Google Scholar] [CrossRef]
  10. 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 2013, 495, 333–338. [Google Scholar] [CrossRef]
  11. Xiao, M.S.; Ai, Y.; Wilusz, J.E. Biogenesis and Functions of Circular RNAs Come into Focus. Trends Cell Biol. 2020, 30, 226–240. [Google Scholar] [CrossRef] [PubMed]
  12. Hansen, T.B.; Jensen, T.I.; Clausen, B.H.; Bramsen, J.B.; Finsen, B.; Damgaard, C.K.; Kjems, J. Natural RNA circles function as efficient microRNA sponges. Nature 2013, 495, 384–388. [Google Scholar] [CrossRef] [PubMed]
  13. Piwecka, M.; Glazar, P.; Hernandez-Miranda, L.R.; Memczak, S.; Wolf, S.A.; Rybak-Wolf, A.; Filipchyk, A.; Klironomos, F.; Cerda Jara, C.A.; Fenske, P.; et al. Loss of a mammalian circular RNA locus causes miRNA deregulation and affects brain function. Science 2017, 357, eaam8526. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Rybak-Wolf, A.; Stottmeister, C.; Glazar, P.; Jens, M.; Pino, N.; Giusti, S.; Hanan, M.; Behm, M.; Bartok, O.; Ashwal-Fluss, R.; et al. Circular RNAs in the Mammalian Brain Are Highly Abundant, Conserved, and Dynamically Expressed. Mol. Cell 2015, 58, 870–885. [Google Scholar] [CrossRef] [Green Version]
  15. Barrett, S.P.; Salzman, J. Circular RNAs: Analysis, expression and potential functions. Development 2016, 143, 1838–1847. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Hanan, M.; Soreq, H.; Kadener, S. CircRNAs in the brain. RNA Biol. 2017, 14, 1028–1034. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Chen, W.; Schuman, E. Circular RNAs in Brain and Other Tissues: A Functional Enigma. Trends Neurosci. 2016, 39, 597–604. [Google Scholar] [CrossRef] [PubMed]
  18. Lipscombe, D.; Lopez Soto, E.J. Alternative splicing of neuronal genes: New mechanisms and new therapies. Curr. Opin. Neurobiol. 2019, 57, 26–31. [Google Scholar] [CrossRef] [PubMed]
  19. You, X.; Vlatkovic, I.; Babic, A.; Will, T.; Epstein, I.; Tushev, G.; Akbalik, G.; Wang, M.; Glock, C.; Quedenau, C.; et al. Neural circular RNAs are derived from synaptic genes and regulated by development and plasticity. Nat. Neurosci. 2015, 18, 603–610. [Google Scholar] [CrossRef] [Green Version]
  20. Martínez-Iglesias, O.; Naidoo, V.; Cacabelos, N.; Cacabelos, R. Epigenetic Biomarkers as Diagnostic Tools for Neurodegenerative Disorders. Int. J. Mol. Sci. 2021, 23, 13. [Google Scholar] [CrossRef] [PubMed]
  21. Li, Q.S.; Vasanthakumar, A.; Davis, J.W.; Idler, K.B.; Nho, K.; Waring, J.F.; Saykin, A.J. Association of peripheral blood DNA methylation level with Alzheimer’s disease progression. Clin. Epig. 2021, 13, 191. [Google Scholar] [CrossRef] [PubMed]
  22. Mano, T.; Sato, K.; Ikeuchi, T.; Toda, T.; Iwatsubo, T.; Iwata, A. Peripheral Blood BRCA1 Methylation Positively Correlates with Major Alzheimer’s Disease Risk Factors. J. Prev. Alzheimers Dis. 2021, 8, 477–482. [Google Scholar] [CrossRef] [PubMed]
  23. Jia, L.; Zhu, M.; Yang, J.; Pang, Y.; Wang, Q.; Li, Y.; Li, T.; Li, F.; Wei, Y. Prediction of P-tau/Aβ42 in the cerebrospinal fluid with blood microRNAs in Alzheimer’s disease. BMC Med. 2021, 19, 264. [Google Scholar] [CrossRef] [PubMed]
  24. Yuen, S.C.; Liang, X.; Zhu, H.; Jia, Y.; Leung, S.W. Prediction of differentially expressed microRNAs in blood as potential biomarkers for Alzheimer’s disease by meta-analysis and adaptive boosting ensemble learning. Alzheimers Res. Ther. 2021, 13, 126. [Google Scholar] [CrossRef] [PubMed]
  25. Dube, U.; Del-Aguila, J.L.; Li, Z.; Budde, J.P.; Jiang, S.; Hsu, S.; Ibanez, L.; Fernandez, M.V.; Farias, F.; Norton, J.; et al. An atlas of cortical circular RNA expression in Alzheimer disease brains demonstrates clinical and pathological associations. Nat. Neurosci. 2019, 22, 1903–1912. [Google Scholar] [CrossRef]
  26. Cervera-Carles, L.; Dols-Icardo, O.; Molina-Porcel, L.; Alcolea, D.; Cervantes-Gonzalez, A.; Muñoz-Llahuna, L.; Clarimon, J. Assessing circular RNAs in Alzheimer’s disease and frontotemporal lobar degeneration. Neurobiol. Aging 2020, 92, 7–11. [Google Scholar] [CrossRef]
  27. Zhang, Y.; Zhao, Y.; Liu, Y.; Wang, M.; Yu, W.; Zhang, L. Exploring the regulatory roles of circular RNAs in Alzheimer’s disease. Transl. Neurodegener. 2020, 9, 35. [Google Scholar] [CrossRef]
  28. Shi, Z.; Chen, T.; Yao, Q.; Zheng, L.; Zhang, Z.; Wang, J.; Hu, Z.; Cui, H.; Han, Y.; Han, X.; et al. The circular RNA ciRS-7 promotes APP and BACE1 degradation in an NF-kappaB-dependent manner. FEBS J. 2017, 284, 1096–1109. [Google Scholar] [CrossRef] [Green Version]
  29. Lukiw, W.J. Circular RNA (circRNA) in Alzheimer’s disease (AD). Front. Genet. 2013, 4, 307. [Google Scholar] [CrossRef] [Green Version]
  30. Zhao, Y.; Alexandrov, P.N.; Jaber, V.; Lukiw, W.J. Deficiency in the Ubiquitin Conjugating Enzyme UBE2A in Alzheimer’s Disease (AD) is Linked to Deficits in a Natural Circular miRNA-7 Sponge (circRNA; ciRS-7). Genes 2016, 7, 116. [Google Scholar] [CrossRef] [Green Version]
  31. Lu, Y.; Tan, L.; Wang, X. Circular HDAC9/microRNA-138/Sirtuin-1 Pathway Mediates Synaptic and Amyloid Precursor Protein Processing Deficits in Alzheimer’s Disease. Neurosci. Bull. 2019, 35, 877–888. [Google Scholar] [CrossRef] [PubMed]
  32. Diling, C.; Yinrui, G.; Longkai, Q.; Xiaocui, T.; Yadi, L.; Xin, Y.; Guoyan, H.; Ou, S.; Tianqiao, Y.; Dongdong, W.; et al. Circular RNA NF1-419 enhances autophagy to ameliorate senile dementia by binding Dynamin-1 and Adaptor protein 2 B1 in AD-like mice. Aging 2019, 11, 12002–12031. [Google Scholar] [CrossRef] [PubMed]
  33. Yang, H.; Wang, H.; Shang, H.; Chen, X.; Yang, S.; Qu, Y.; Ding, J.; Li, X. Circular RNA circ_0000950 promotes neuron apoptosis, suppresses neurite outgrowth and elevates inflammatory cytokines levels via directly sponging miR-103 in Alzheimer’s disease. Cell Cycle 2019, 18, 2197–2214. [Google Scholar] [CrossRef] [PubMed]
  34. Cochran, K.R.; Veeraraghavan, K.; Kundu, G.; Mazan-Mamczarz, K.; Coletta, C.; Thambisetty, M.; Gorospe, M.; De, S. Systematic Identification of circRNAs in Alzheimer’s Disease. Genes 2021, 12, 1258. [Google Scholar] [CrossRef]
  35. Welden, J.R.; van Doorn, J.; Nelson, P.T.; Stamm, S. The human MAPT locus generates circular RNAs. Biochim. Biophys Acta 2018, 1864, 2753–2760. [Google Scholar] [CrossRef]
  36. Lee, E.G.; Tulloch, J.; Chen, S.; Leong, L.; Saxton, A.D.; Kraemer, B.; Darvas, M.; Keene, C.D.; Shutes-David, A.; Todd, K.; et al. Redefining transcriptional regulation of the APOE gene and its association with Alzheimer’s disease. PLoS ONE 2020, 15, e0227667. [Google Scholar] [CrossRef]
  37. Liu, L.; Chen, X.; Chen, Y.H.; Zhang, K. Identification of Circular RNA hsa_Circ_0003391 in Peripheral Blood Is Potentially Associated With Alzheimer’s Disease. Front. Aging Neurosci. 2020, 12, 601965. [Google Scholar] [CrossRef]
  38. Li, Y.; Lv, Z.; Zhang, J.; Ma, Q.; Li, Q.; Song, L.; Gong, L.; Zhu, Y.; Li, X.; Hao, Y.; et al. Profiling of differentially expressed circular RNAs in peripheral blood mononuclear cells from Alzheimer’s disease patients. Metab. Brain Dis. 2020, 35, 201–213. [Google Scholar] [CrossRef]
  39. Li, Y.; Fan, H.; Sun, J.; Ni, M.; Zhang, L.; Chen, C.; Hong, X.; Fang, F.; Zhang, W.; Ma, P. Circular RNA expression profile of Alzheimer’s disease and its clinical significance as biomarkers for the disease risk and progression. Int. J. Biochem. Cell Biol. 2020, 123, 105747. [Google Scholar] [CrossRef]
  40. Lo, I.; Hill, J.; Vilhjálmsson, B.J.; Kjems, J. Linking the association between circRNAs and Alzheimer’s disease progression by multi-tissue circular RNA characterization. RNA Biol. 2020, 17, 1789–1797. [Google Scholar] [CrossRef]
  41. Jeck, W.R.; Sorrentino, J.A.; Wang, K.; Slevin, M.K.; Burd, C.E.; Liu, J.; Marzluff, W.F.; Sharpless, N.E. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA 2013, 19, 141–157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Salzman, J.; Chen, R.E.; Olsen, M.N.; Wang, P.L.; Brown, P.O. Cell-type specific features of circular RNA expression. PLoS Genet. 2013, 9, e1003777. [Google Scholar] [CrossRef]
  43. Gómez-Isla, T.; Price, J.L.; McKeel, D.W.; Morris, J.C.; Growdon, J.H.; Hyman, B.T. Profound loss of layer II entorhinal cortex neurons occurs in very mild Alzheimer’s disease. J. Neurosci. 1996, 16, 4491–4500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Bevilaqua, L.R.; Rossato, J.I.; Bonini, J.S.; Myskiw, J.C.; Clarke, J.R.; Monteiro, S.; Lima, R.H.; Medina, J.H.; Cammarota, M.; Izquierdo, I. The role of the entorhinal cortex in extinction: Influences of aging. Neural Plast. 2008, 2008, 595282. [Google Scholar] [CrossRef] [Green Version]
  45. Liang, W.S.; Dunckley, T.; Beach, T.G.; Grover, A.; Mastroeni, D.; Walker, D.G.; Caselli, R.J.; Kukull, W.A.; McKeel, D.; Morris, J.C.; et al. Gene expression profiles in anatomically and functionally distinct regions of the normal aged human brain. Physiol. Genom. 2007, 28, 311–322. [Google Scholar] [CrossRef] [PubMed]
  46. Dudekula, D.B.; Panda, A.C.; Grammatikakis, I.; De, S.; Abdelmohsen, K.; Gorospe, M. CircInteractome: A web tool for exploring circular RNAs and their interacting proteins and microRNAs. RNA Biol. 2016, 13, 34–42. [Google Scholar] [CrossRef] [Green Version]
  47. Lugli, G.; Cohen, A.M.; Bennett, D.A.; Shah, R.C.; Fields, C.J.; Hernandez, A.G.; Smalheiser, N.R. Plasma Exosomal miRNAs in Persons with and without Alzheimer Disease: Altered Expression and Prospects for Biomarkers. PLoS ONE 2015, 10, e0139233. [Google Scholar] [CrossRef] [Green Version]
  48. Riancho, J.; Vázquez-Higuera, J.L.; Pozueta, A.; Lage, C.; Kazimierczak, M.; Bravo, M.; Calero, M.; Gonalezález, A.; Rodríguez, E.; Lleó, A.; et al. MicroRNA Profile in Patients with Alzheimer’s Disease: Analysis of miR-9-5p and miR-598 in Raw and Exosome Enriched Cerebrospinal Fluid Samples. J. Alzheimers Dis. 2017, 57, 483–491. [Google Scholar] [CrossRef]
  49. Guévremont, D.; Tsui, H.; Knight, R.; Fowler, C.J.; Masters, C.L.; Martins, R.N.; Abraham, W.C.; Tate, W.P.; Cutfield, N.J.; Williams, J.M. Plasma microRNA vary in association with the progression of Alzheimer’s disease. Alzheimers Dement. 2022, 14, e12251. [Google Scholar] [CrossRef]
  50. Ravanidis, S.; Bougea, A.; Papagiannakis, N.; Koros, C.; Simitsi, A.M.; Pachi, I.; Breza, M.; Stefanis, L.; Doxakis, E. Validation of differentially expressed brain-enriched microRNAs in the plasma of PD patients. Ann. Clin. Transl. Neurol. 2020, 7, 1594–1607. [Google Scholar] [CrossRef]
  51. Rahmani, S.; Kadkhoda, S.; Ghafouri-Fard, S. Synaptic plasticity and depression: The role of miRNAs dysregulation. Mol. Biol. Rep. 2022, 49, 9759–9765. [Google Scholar] [CrossRef] [PubMed]
  52. Dong, Y.; Xiong, J.; Ji, L.; Xue, X. MiR-421 Aggravates Neurotoxicity and Promotes Cell Death in Parkinson’s Disease Models by Directly Targeting MEF2D. Neurochem. Res. 2021, 46, 299–308. [Google Scholar] [CrossRef] [PubMed]
  53. Cui, H.; Zhang, W. The Neuroprotective Effect of miR-136 on Pilocarpine-Induced Temporal Lobe Epilepsy Rats by Inhibiting Wnt/. Comput. Math. Methods Med. 2022, 2022, 1938205. [Google Scholar] [CrossRef]
  54. Raheja, R.; Regev, K.; Healy, B.C.; Mazzola, M.A.; Beynon, V.; Von Glehn, F.; Paul, A.; Diaz-Cruz, C.; Gholipour, T.; Glanz, B.I.; et al. Correlating serum micrornas and clinical parameters in amyotrophic lateral sclerosis. Muscle Nerve 2018, 58, 261–269. [Google Scholar] [CrossRef] [PubMed]
  55. Vlachos, I.S.; Zagganas, K.; Paraskevopoulou, M.D.; Georgakilas, G.; Karagkouni, D.; Vergoulis, T.; Dalamagas, T.; Hatzigeorgiou, A.G. DIANA-miRPath v3.0: Deciphering microRNA function with experimental support. Nucleic Acids Res. 2015, 43, W460–W466. [Google Scholar] [CrossRef]
  56. Boese, A.S.; Saba, R.; Campbell, K.; Majer, A.; Medina, S.; Burton, L.; Booth, T.F.; Chong, P.; Westmacott, G.; Dutta, S.M.; et al. MicroRNA abundance is altered in synaptoneurosomes during prion disease. Mol. Cell Neurosci. 2016, 71, 13–24. [Google Scholar] [CrossRef]
  57. Ben Halima, S.; Siegel, G.; Rajendran, L. miR-186 in Alzheimer’s disease: A big hope for a small RNA? J. Neurochem. 2016, 137, 308–311. [Google Scholar] [CrossRef] [Green Version]
  58. Rahman, M.R.; Islam, T.; Zaman, T.; Shahjaman, M.; Karim, M.R.; Huq, F.; Quinn, J.M.W.; Holsinger, R.M.D.; Gov, E.; Moni, M.A. Identification of molecular signatures and pathways to identify novel therapeutic targets in Alzheimer’s disease: Insights from a systems biomedicine perspective. Genomics 2020, 112, 1290–1299. [Google Scholar] [CrossRef]
  59. You, J.; Qian, F.; Huang, Y.; Guo, Y.; Lv, Y.; Yang, Y.; Lu, X.; Guo, T.; Wang, J.; Gu, B. lncRNA WT1-AS attenuates hypoxia/ischemia-induced neuronal injury during cerebral ischemic stroke via miR-186-5p/XIAP axis. Open Med. 2022, 17, 1338–1349. [Google Scholar] [CrossRef]
  60. Cai, S.C.; Li, X.P.; Li, X.; Tang, G.Y.; Yi, L.M.; Hu, X.S. Oleanolic Acid Inhibits Neuronal Pyroptosis in Ischaemic Stroke by Inhibiting miR-186-5p Expression. Exp. Neurobiol. 2021, 30, 401–414. [Google Scholar] [CrossRef]
  61. Kim, J.; Yoon, H.; Chung, D.E.; Brown, J.L.; Belmonte, K.C. miR-186 is decreased in aged brain and suppresses BACE1 expression. J. Neurochem. 2016, 137, 436–445. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Tiwari, D.; Brager, D.H.; Rymer, J.K.; Bunk, A.T.; White, A.R.; Elsayed, N.A.; Krzeski, J.C.; Snider, A.; Schroeder Carter, L.M.; Danzer, S.C.; et al. MicroRNA inhibition upregulates hippocampal A-type potassium current and reduces seizure frequency in a mouse model of epilepsy. Neurobiol. Dis. 2019, 130, 104508. [Google Scholar] [CrossRef] [PubMed]
  63. Dong, X.; He, X.; Yang, L.; Li, Q.; Xu, Y. Inhibition of miR-421 Preserves Mitochondrial Function and Protects against Parkinson’s Disease Pathogenesis via Pink1/Parkin-Dependent Mitophagy. Dis. Markers 2022, 2022, 5186252. [Google Scholar] [CrossRef]
  64. Hu, F.; Shao, L.; Zhang, J.; Zhang, H.; Wen, A.; Zhang, P. Knockdown of ZFAS1 Inhibits Hippocampal Neurons Apoptosis and Autophagy by Activating the PI3K/AKT Pathway via Up-regulating miR-421 in Epilepsy. Neurochem. Res. 2020, 45, 2433–2441. [Google Scholar] [CrossRef] [PubMed]
  65. Wen, X.; Han, X.R.; Wang, Y.J.; Wang, S.; Shen, M.; Zhang, Z.F.; Fan, S.H.; Shan, Q.; Wang, L.; Li, M.Q.; et al. MicroRNA-421 suppresses the apoptosis and autophagy of hippocampal neurons in epilepsy mice model by inhibition of the TLR/MYD88 pathway. J. Cell. Physiol. 2018, 233, 7022–7034. [Google Scholar] [CrossRef] [PubMed]
  66. Chen, J.; Zhao, H.; Huang, Y.; Li, Y.; Fan, J.; Wang, R.; Han, Z.; Yang, Z.; Wu, L.; Wu, D.; et al. Dysregulation of Principal Circulating miRNAs in Non-human Primates Following Ischemic Stroke. Front. Neurosci. 2021, 15, 738576. [Google Scholar] [CrossRef]
  67. Li, N.; Zhang, D.; Guo, H.; Yang, Q.; Li, P.; He, Y. Inhibition of circ_0004381 improves cognitive function via miR-647/PSEN1 axis in an Alzheimer disease mouse model. J. Neuropathol. Exp. Neurol. 2022, 82, 84–92. [Google Scholar] [CrossRef]
  68. Li, K.; Shen, L.; Zheng, P.; Wang, Y.; Wang, L.; Meng, X.; Lv, Y.; Xue, Z.; Guo, X.; Zhang, A.; et al. Identification of MicroRNAs as potential biomarkers for detecting ischemic stroke. Genes Genom. 2022, 44, 9–17. [Google Scholar] [CrossRef]
  69. Rademakers, R.; Eriksen, J.L.; Baker, M.; Robinson, T.; Ahmed, Z.; Lincoln, S.J.; Finch, N.; Rutherford, N.J.; Crook, R.J.; Josephs, K.A.; et al. Common variation in the miR-659 binding-site of GRN is a major risk factor for TDP43-positive frontotemporal dementia. Hum. Mol. Genet. 2008, 17, 3631–3642. [Google Scholar] [CrossRef] [Green Version]
  70. Piscopo, P.; Grasso, M.; Fontana, F.; Crestini, A.; Puopolo, M.; Del Vescovo, V.; Venerosi, A.; Calamandrei, G.; Vencken, S.F.; Greene, C.M.; et al. Reduced miR-659-3p Levels Correlate with Progranulin Increase in Hypoxic Conditions: Implications for Frontotemporal Dementia. Front. Mol. Neurosci. 2016, 9, 31. [Google Scholar] [CrossRef] [Green Version]
  71. Ibrahim, H.M.; AlZahrani, A.; Hanieh, H.; Ahmed, E.A.; Thirugnanasambantham, K. MicroRNA-7188-5p and miR-7235 regulates Multiple sclerosis in an experimental mouse model. Mol. Immunol. 2021, 139, 157–167. [Google Scholar] [CrossRef] [PubMed]
  72. Miller-Delaney, S.F.; Bryan, K.; Das, S.; McKiernan, R.C.; Bray, I.M.; Reynolds, J.P.; Gwinn, R.; Stallings, R.L.; Henshall, D.C. Differential DNA methylation profiles of coding and non-coding genes define hippocampal sclerosis in human temporal lobe epilepsy. Brain 2015, 138, 616–631. [Google Scholar] [CrossRef] [Green Version]
  73. Karisetty, B.C.; Bhatnagar, A.; Armour, E.M.; Beaver, M.; Zhang, H.; Elefant, F. Amyloid-β Peptide Impact on Synaptic Function and Neuroepigenetic Gene Control Reveal New Therapeutic Strategies for Alzheimer’s Disease. Front. Mol. Neurosci. 2020, 13, 577622. [Google Scholar] [CrossRef] [PubMed]
  74. He, C.; Huang, Z.S.; Yu, C.C.; Wang, H.H.; Zhou, H.; Kong, L.H. Epigenetic Regulation of Amyloid-beta Metabolism in Alzheimer’s Disease. Curr. Med. Sci. 2020, 40, 1022–1030. [Google Scholar] [CrossRef] [PubMed]
  75. Wang, Q.; Huang, X.; Su, Y.; Yin, G.; Wang, S.; Yu, B.; Li, H.; Qi, J.; Chen, H.; Zeng, W.; et al. Activation of Wnt/β-catenin pathway mitigates blood-brain barrier dysfunction in Alzheimer’s disease. Brain 2022, 145, 4474–4488. [Google Scholar] [CrossRef] [PubMed]
  76. Hegde, A.N.; Smith, S.G.; Duke, L.M.; Pourquoi, A.; Vaz, S. Perturbations of Ubiquitin-Proteasome-Mediated Proteolysis in Aging and Alzheimer’s Disease. Front. Aging Neurosci. 2019, 11, 324. [Google Scholar] [CrossRef]
  77. Zhang, L.; Qi, Z.; Li, J.; Li, M.; Du, X.; Wang, S.; Zhou, G.; Xu, B.; Liu, W.; Xi, S.; et al. Roles and Mechanisms of Axon-Guidance Molecules in Alzheimer’s Disease. Mol. Neurobiol. 2021, 58, 3290–3307. [Google Scholar] [CrossRef]
  78. Skaper, S.D.; Facci, L.; Zusso, M.; Giusti, P. Synaptic Plasticity, Dementia and Alzheimer Disease. CNS Neurol. Disord. Drug Targets 2017, 16, 220–233. [Google Scholar] [CrossRef]
  79. Snow, A.D.; Cummings, J.A.; Lake, T. The Unifying Hypothesis of Alzheimer’s Disease: Heparan Sulfate Proteoglycans/Glycosaminoglycans Are Key as First Hypothesized Over 30 Years Ago. Front. Aging Neurosci. 2021, 13, 710683. [Google Scholar] [CrossRef]
  80. Sandwall, E.; O’Callaghan, P.; Zhang, X.; Lindahl, U.; Lannfelt, L.; Li, J.P. Heparan sulfate mediates amyloid-beta internalization and cytotoxicity. Glycobiology 2010, 20, 533–541. [Google Scholar] [CrossRef] [Green Version]
  81. Beckman, M.; Holsinger, R.M.; Small, D.H. Heparin activates beta-secretase (BACE1) of Alzheimer’s disease and increases autocatalysis of the enzyme. Biochemistry 2006, 45, 6703–6714. [Google Scholar] [CrossRef] [PubMed]
  82. Bukke, V.N.; Archana, M.; Villani, R.; Romano, A.D.; Wawrzyniak, A.; Balawender, K.; Orkisz, S.; Beggiato, S.; Serviddio, G.; Cassano, T. The Dual Role of Glutamatergic Neurotransmission in Alzheimer’s Disease: From Pathophysiology to Pharmacotherapy. Int. J. Mol. Sci. 2020, 21, 7452. [Google Scholar] [CrossRef]
  83. Bell, J.E.; Alafuzoff, I.; Al-Sarraj, S.; Arzberger, T.; Bogdanovic, N.; Budka, H.; Dexter, D.T.; Falkai, P.; Ferrer, I.; Gelpi, E.; et al. Management of a twenty-first century brain bank: Experience in the BrainNet Europe consortium. Acta Neuropathol. 2008, 115, 497–507. [Google Scholar] [CrossRef] [PubMed]
  84. Montine, T.J.; Phelps, C.H.; Beach, T.G.; Bigio, E.H.; Cairns, N.J.; Dickson, D.W.; Duyckaerts, C.; Frosch, M.P.; Masliah, E.; Mirra, S.S.; et al. National Institute on Aging-Alzheimer’s Association guidelines for the neuropathologic assessment of Alzheimer’s disease: A practical approach. Acta Neuropathol. 2012, 123, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Braak, H.; Alafuzoff, I.; Arzberger, T.; Kretzschmar, H.; Del Tredici, K. Staging of Alzheimer disease-associated neurofibrillary pathology using paraffin sections and immunocytochemistry. Acta Neuropathol. 2006, 112, 389–404. [Google Scholar] [CrossRef] [Green Version]
  86. Mirra, S.S.; Heyman, A.; McKeel, D.; Sumi, S.M.; Crain, B.J.; Brownlee, L.M.; Vogel, F.S.; Hughes, J.P.; van Belle, G.; Berg, L. The Consortium to Establish a Registry for Alzheimer’s Disease (CERAD). Part II. Standardization of the neuropathologic assessment of Alzheimer’s disease. Neurology 1991, 41, 479–486. [Google Scholar] [CrossRef]
  87. Kent, W.J.; Sugnet, C.W.; Furey, T.S.; Roskin, K.M.; Pringle, T.H.; Zahler, A.M.; Haussler, D. The human genome browser at UCSC. Genome Res. 2002, 12, 996–1006. [Google Scholar] [CrossRef] [Green Version]
  88. Navarro Gonzalez, J.; Zweig, A.S.; Speir, M.L.; Schmelter, D.; Rosenbloom, K.R.; Raney, B.J.; Powell, C.C.; Nassar, L.R.; Maulding, N.D.; Lee, C.M.; et al. The UCSC Genome Browser database: 2021 update. Nucleic Acids Res. 2021, 49, D1046–D1057. [Google Scholar] [CrossRef]
  89. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
Figure 1. Identification of circAPP (hsa_circ_0007556). (A) Schematic representation of circAPP (hsa_circ_0007556) and APP mRNA, arrows represent primers (blue arrows for circAPP (hsa_circ_0007556) divergent and green arrows for APP mRNA). (B) Electrophoresis agarose gel of circAPP (hsa_circ_0007556), arrow points the circAPP (hsa_circ_0007556) PCR product. (C) Electropherogram showing the circAPP (hsa_circ_0007556) covalent junction.
Figure 1. Identification of circAPP (hsa_circ_0007556). (A) Schematic representation of circAPP (hsa_circ_0007556) and APP mRNA, arrows represent primers (blue arrows for circAPP (hsa_circ_0007556) divergent and green arrows for APP mRNA). (B) Electrophoresis agarose gel of circAPP (hsa_circ_0007556), arrow points the circAPP (hsa_circ_0007556) PCR product. (C) Electropherogram showing the circAPP (hsa_circ_0007556) covalent junction.
Ijms 24 04308 g001
Figure 2. Expression of circAPP (hsa_circ_0007556) and APP mRNA and as a function of the ABC score. (A) The box-plot represents the percentage of relative expression of circAPP (hsa_circ_0007556) and APP mRNA in patients and controls with respect to the geometric mean of the housekeeping genes GAPDH and ACTB. (B) The box-plot represents the percentage of the relative expression of circAPP (hsa_circ_0007556) and APP mRNA according to the ABC score groups, with respect to the geometric mean of the housekeeping genes GAPDH and ACTB. * p-value < 0.05, ** p-value < 0.01 and *** p-value < 0.001.
Figure 2. Expression of circAPP (hsa_circ_0007556) and APP mRNA and as a function of the ABC score. (A) The box-plot represents the percentage of relative expression of circAPP (hsa_circ_0007556) and APP mRNA in patients and controls with respect to the geometric mean of the housekeeping genes GAPDH and ACTB. (B) The box-plot represents the percentage of the relative expression of circAPP (hsa_circ_0007556) and APP mRNA according to the ABC score groups, with respect to the geometric mean of the housekeeping genes GAPDH and ACTB. * p-value < 0.05, ** p-value < 0.01 and *** p-value < 0.001.
Ijms 24 04308 g002
Table 1. Correlation between the Method of Thal, Braak and Braak classification, Method of CERAD, Global average area of Aβ deposits and circAPP or APP mRNA expression levels.
Table 1. Correlation between the Method of Thal, Braak and Braak classification, Method of CERAD, Global average area of Aβ deposits and circAPP or APP mRNA expression levels.
Method of ThalBraak and Braak ClassificationMethod of CERADABC ScoreGlobal Average Area of Aβ Deposits
log(circAPP)Correlation coeficient (Spearman’s Rho)−0.496 **−0.503 **−0.499 **−0.515 **−0.480 *
Sig. (bilateral)0.0010.0010.0010.0000.024
N4344434322
log(APP mRNA)Correlation coeficient (Spearman’s Rho)−0.339 *−0.415 **−0.243−0.384 *−0.652 **
Sig. (bilateral)0.0260.0050.1160.0110.001
N4344434322
* p-value < 0.05, ** p-value < 0.01.
Table 2. miRNAs predicted to bind circAPP (hsa_circ_0007556) with circInteractome tool and their implication in neurological disorders.
Table 2. miRNAs predicted to bind circAPP (hsa_circ_0007556) with circInteractome tool and their implication in neurological disorders.
miRNANeurological Disorders RelatedExpression in ADTissue/Fluid StudiedReferences
hsa-miR-1200-
hsa-miR-1208-
hsa-miR-1270-
hsa-miR-1272-
hsa-miR-136Prion disease, epilepsy, depression, PD [50,51,53,56]
hsa-miR-186AD, ischemic stroke-Neuro-2a cells, 7PA2 cells[57,58,59,60,61]
hsa-miR-324-5pADplasma[49,62]
hsa-miR-421PD, epilepsy, ischemic stroke [63,64,65,66]
hsa-miR-518a-5p-
hsa-miR-527-
hsa-miR-549-
hsa-miR-598ADNo detectedCSF[48]
↑(tendency)CSF exosomes
hsa-miR-620-
hsa-miR-646-
hsa-miR-647AD [67]
hsa-miR-659AD, ischemic stroke, frontotemporal dementia, MShippocampal neurons AD mouse model[47,68,69,70,71]
hsa-miR-876-3pEpilepsyPlasma exosomes[72]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Urdánoz-Casado, A.; Sánchez-Ruiz de Gordoa, J.; Robles, M.; Roldan, M.; Macías Conde, M.; Acha, B.; Blanco-Luquin, I.; Mendioroz, M. circRNA from APP Gene Changes in Alzheimer’s Disease Human Brain. Int. J. Mol. Sci. 2023, 24, 4308. https://doi.org/10.3390/ijms24054308

AMA Style

Urdánoz-Casado A, Sánchez-Ruiz de Gordoa J, Robles M, Roldan M, Macías Conde M, Acha B, Blanco-Luquin I, Mendioroz M. circRNA from APP Gene Changes in Alzheimer’s Disease Human Brain. International Journal of Molecular Sciences. 2023; 24(5):4308. https://doi.org/10.3390/ijms24054308

Chicago/Turabian Style

Urdánoz-Casado, Amaya, Javier Sánchez-Ruiz de Gordoa, Maitane Robles, Miren Roldan, Mónica Macías Conde, Blanca Acha, Idoia Blanco-Luquin, and Maite Mendioroz. 2023. "circRNA from APP Gene Changes in Alzheimer’s Disease Human Brain" International Journal of Molecular Sciences 24, no. 5: 4308. https://doi.org/10.3390/ijms24054308

APA Style

Urdánoz-Casado, A., Sánchez-Ruiz de Gordoa, J., Robles, M., Roldan, M., Macías Conde, M., Acha, B., Blanco-Luquin, I., & Mendioroz, M. (2023). circRNA from APP Gene Changes in Alzheimer’s Disease Human Brain. International Journal of Molecular Sciences, 24(5), 4308. https://doi.org/10.3390/ijms24054308

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop