Study population
Among 262 participants from three memory clinics (Inha University Hospital and Yonsei University Severance Hospital) across South Korea who agreed to provide blood samples, 59 subjects were cognitively normal (CN), 105 subjects were diagnosed as MCI, and 98 subjects were diagnosed as probable AD, according to the criteria for clinical diagnosis as previously described [14]. Plasma samples from 15 subjects (7 CN and 8 AD) out of 262 participants and an additional 8 patients with Parkinson’s disease (PD) were used to discover differentially expressed miRNAs able to discriminate AD patients from healthy controls or PD. To test the diagnostic utility of discovered plasma miRNAs, samples from 100 subjects (51 CN, 22 MCI, and 27 AD) recruited from Inha University Hospital were used for comparing the levels of differentially expressed plasma miRNAs. For comparing levels of the differentially expressed miRNAs in plasma EVs, we added 162 plasma samples from other independent cohorts from Yonsei University Hospital and Inha University Hospital.
Participants underwent amyloid imaging via 18F-flutemetamol positron emission tomography (PET). Amyloid PET positivity was determined independently by clinical neurologists. The Institutional Review Board and Ethical Committee of Inha University Hospital approved the study (approval No.: INHAUH2016-06-010-002). Written informed consent to participate was obtained from all volunteers. The demographic characteristics and clinical information of the subjects are summarized in Table 1.
Table 1
Demographic and clinical information of 262 participants.
Diagnosis | n (M:F) | Aβ-PET (+/-, n)a | Age, yr (mean±S.D.) | Education, yr (mean±S.D.) | MMSE (mean±S.D.) |
CN | 59 (18:41) | 4/49 | 68±6.5 | 9.6±4.3 | 27.5±2.1 |
MCI | 105 (37:68) | 43/62 | 73±6.8*** | 8.3±4.9 | 24.3±3.3*** |
AD | 98 (25:73) | 71/22 | 75±8.5*** | 8.1±4.7 | 16.7±5.9***,### |
P value | n.s. b | <0.001 b | <0.001 | n.s. | <0.001 |
Total | 262 (80:182) | 118/133 | 73±7.9 | 8.5±4.8 | 22.1±6.2 |
a6 of CN and 5 of AD did not undergo Aβ-PET analysis. bChi-square test, ***p<0.001 vs. CN, ###p<0.001 vs. MCI by Dunn’s multiple comparison. n.s., not significant. |
Blood collection and plasma preparation
To exclude the contribution of anticoagulants to the pre-analytical variability in the EV yield [15] or the subsequent quantitative polymerase chain reaction (qPCR) [16], we selected the most commonly used anticoagulant, ethylenediaminetetraacetic acid (EDTA) to prepare the plasma. All participants provided their fasting blood, which was drawn from the antecubital vein into EDTA tubes (BD Vacutainer K2E-EDTA; Beckton and Dickinson, Franklin Lake, NJ, USA). Plasma was prepared by immediate centrifugation at 3,000 ⋅ g for 10 min at room temperature (RT). To remove platelets, cell debris, and large particles, the plasma was centrifuged again at 10,000 ⋅ g for 20 min at RT. The supernatant was stored at –70 °C until EV preparation.
Extraction of plasma small RNAs
We collected 1.5 mL blood in PAXgene Blood RNA tubes (PreAnalytiX) from patients with AD (n = 8), patients with PD (n = 8), and from CN (n = 7). After centrifugation at 15,000 rpm, we collected 0.5 mL plasma and added 0.5 mL of extraction buffer. We added 0.2 mL of chloroform to the mixture of plasma and extraction buffer and incubated for 5 min. After centrifugation for 10 min at 15,000 rpm, we transferred the supernatant to a new tube and added 1 ⋅ volume of isopropanol. By filtering the mixed solution through a MinElute column (Qiagen, Venlo, Netherlands), we collected RNAs that were washed with 0.5 mL of Qiagen washing buffer (RPE buffer), and dissolved in DEPC-treated water.
Plasma miRNAome analysis
For small RNA sequencing, the quality of the constructed library was checked using an Agilent Bioanalyzer (Agilent, Santa Clara, CA, USA) and sequenced on the Illumina Nextseq 500 (Illumina, San Diego, CA, USA) platform. The quality and band size of the libraries were assessed using an Agilent 2100 Bioanalyzer (Agilent). Libraries were quantified by qPCR using the CFX96 Real-Time System (Bio-Rad, Hercules, CA, USA). After normalization, library sequencing was conducted on the Nextseq 500 system (Illumina) with 75 bp single end reads. We obtained 41, 36, and 35 million clean reads on average for patients with AD, patients with PD, and from CN, respectively. Potentially existing sequencing adapters and raw quality bases in the low reads were trimmed using Skewer and FastaQC. The cleaned high-quality reads, after trimming the low-quality bases and sequencing adapters were mapped to the reference genome using QuickMIRseq. To quantify the mapped reads on the reference genome into the miRNA expression values, QuickMIRseq with default options was used. The miRNA and hairpin annotation of the reference genome from miRBase was used and the expression values were calculated in Reads Per Ten Million (RPTM) units. The differentially expressed miRNAs between the two selected biological conditions were analyzed by NOISeq in Bioconductor. The expression levels of miRNAs were calculated by normalizing the miRNA counts with the total number of clean reads in the small RNA libraries.
EV isolation, extraction of RNAs and amplification of specific EV miRNAs
For isolation of EVs in plasma, we modified the manufacturer’s provided protocol, as previously described [17]. Briefly, we used a commercially available PP-based kit (miRCURY Exosome Plasma Kit, #76603, Qiagen, Venlo, Netherlands) with pre-incubation of plasma with thrombin (#36402.01, SERVA, Heidelberg, Germany) and proteinase K (0.5 mg/mL; Ambion, Austin, TX, USA). We used the serum-like supernatant from plasma samples (0.6 mL) pre-incubated with thrombin and centrifuged for defibrinating. Then, the supernatant was mixed with proteinase K and incubated at 37 °C for 10 min, followed by incubation with precipitation buffer at 4 °C for 1 h. RNA pellets were obtained following centrifugation at 1,000 ⋅ g for 5 min. Subsequently, we isolated total RNA using a commercial RNA extraction kit (miRNeasy Plasma Advanced Kit, #217204, Qiagen) following the manufacturer’s instructions. EV pellets were lysed with 0.5 mL RNase-free water and the lysed EVs were mixed with lysis RPL buffer (150 µL) including 0.625 ng synthetic RNA (#339390, Qiagen), followed by precipitation in 50 µL of RPP buffer and centrifugation at 12,000 ⋅ g for 3 min at RT. The supernatants were transferred to new tubes and mixed with an equal volume of isopropanol. The mixtures were eluted through a miRNA spin column (RNeasy UCP MinElute column, Qiagen) after washing with 20 µL of RNase-free water. Using 2 µL of the extracted RNA, cDNA was synthesized with a miRCURY LNA RT Kit (#339340, Qiagen) following the manufacturer’s instructions. cDNA was processed through a reverse transcription step at 42 °C for 60 min and an inactivation reaction step at 95 °C for 5 min on a PCR Thermal Cycler (TP600, TAKARA BIO INC., Japan). After reverse transcription, specific miRNAs were amplified using cDNA with 10-fold dilution. We used primer sets for miRNAs from Qiagen, and carried out qPCR under conditions of the miRCURY LNA SYBR Green PCR kit (#339345, Qiagen) on a CFX Connect Real-Time PCR Detection System (#1855201, Bio-Rad, Hercules, CA, USA). The PCR cycling conditions consisted of PCR initial heat activation at 95 °C for 2 min and 45 cycles of 2-step cycling with denaturation at 95 °C for 10 s and combined annealing/extension at 56 °C for 60 s, followed by melting curve analysis from 60 °C to 95 °C. The quantitative real-time PCR (qRT-PCR) data were normalized to UniSp2 as an exogenous reference gene. The primer sequences are summarized in
miRNA target prediction and pathway enrichment analysis
Using the results from differential miRNA expression data, gene–miRNA interaction was predicted using miRDB, which is an online database for miRNA target prediction and functional annotations (version 6.0) [18]. Candidate transcripts with scores ≥ 80 are presented as predicted targets for further analysis. Pathway enrichment analysis was employed to investigate the function of significantly differentially expressed miRNAs. Metascape (http://metascape.org) was used to carry out gene ontology (GO) term analysis for enriched pathways, providing gene list annotation and analysis source [19]. Pairwise similarities between any two enriched terms are computed based on a Kappa-test score (threshold of a 0.3 similarity). In addition, to explore further the relationship between differentially expressed genes (DEGs) at the protein level, we carried out protein–protein interaction (PPI) network analysis using the following databases: STRING, BioGrid, OmniPath, and InWeb IM. If the network contains between 3 and 1000 proteins with physical interaction, the plug-ins Molecular Complex Detection (MCODE) algorithm was applied to identify densely connected network components.
Cell culture and transfection
H4 cells overexpressing the Swedish FAD mutant APP695 (H4-APPswe) were grown in Dulbecco’s modified Eagle’s medium (Invitrogen) with 10% fetal bovine serum (FBS), 25 mM glucose, 400 µg/mL G418, 100 U/mL penicillin, and 100 µg/mL streptomycin at 37 °C in a humidified atmosphere of 95% air and 5% CO2. For transfection of miRNA, 400,000 cells per well were plated in 6-well plates. The next day, cells were transfected with 100 nM of pre-miRNAs (BIONEER, Daejeon, Korea) or scrambled sequence using Lipofectamine RNAiMAX (Invitrogen, MA, USA) following the manufacturer’s instructions. The transfected cells were grown in complete medium for 24 h and refed with serum-free medium for an additional 24 h. Harvested cells and culture media (CM) were used for immunoblots, qRT-PCR, or ELISA.
Measurement of Aβ in culture media
We collected CM of H4-APPswe cells at 24 h after miRNA transfection. Media were centrifuged to remove particulates and stored at –80 °C until use. The quantity of Aβ40 and Aβ42 was measured using commercially available ELISA kits (Invitrogen, MA, USA) according to the manufacturer’s instructions. The CM was diluted appropriately and transferred into the supplied microplate pre-coated with the target-specific capture antibody. Next, the sandwich was formed by incubating with a detection antibody for 3 h at RT. Enzyme labeled antibody was added and reacted with substrate solution. Absorbance was read at 450 nm.
Amplification of targeted mRNAs
For RNA extraction, whole cells were collected using the Trizol reagent (Ambion, Carlsbad, CA, USA), and subjected to RNA isolation as previously described [20]. Total RNA quality and quantity were measured using a NanoDrop 1000 (Thermo Fisher Scientific, MA, USA), before 20 µL reactions containing RNA template (1 µg) were converted to cDNA using the PrimeScript™ 1st strand cDNA Synthesis Kit (Takara, Japan). The mRNA expression levels of three miRNA target genes (ADAM10, BDNF, and JAG1) were estimated via qRT-PCR. Amplification reactions contained 10 µL of 2X SYBR supermix (#170-8882, Bio-Rad), each forward and reverse primer (1 pM), 4 µL of nuclease-free water, and 2 µL of cDNA. The qRT-PCR was carried out under the following conditions: Initial step at 95 °C for 10 min; 40 cycles of 95 °C 15 s; 58 °C 25 s, and 72 °C 35 s, followed by melting curve analysis from 60 °C to 95 °C. Primers were designed using the NCBI primer blast tool and synthesized by BIONIC INC. (Seoul, ROK). The primer sequences are summarized in Supplementary Table 2. The relative expression level was calculated using the 2−ΔΔCt method (ΔCt = CtTarget gene – CtGAPDH; ΔΔCt = ΔCtTarget miRNA – ΔCtscrambled miRNA), after normalization with the glyceraldehyde-3-phosphate-dehydrogenase (GADPH) gene as a housekeeping gene.
SDS-PAGE and immunoblot analysis
Cells were treated with miRNA for 24 h and collected to evaluate the expression of amyloid precursor protein (APP) and its C-terminal fragments (CTFs). To prepare whole-cell lysates, cells were washed with ice-cold PBS, scraped into radioimmunoprecipitation assay (RIPA) buffer (GE Healthcare Bio-Sciences, Pittsburgh, PA, USA) containing a cocktail of both protease inhibitors and phosphatase inhibitors (Sigma-Aldrich, St. Louis, MO, USA), sonicated to homogenize the cell suspension, and centrifuged to remove cell debris. Total protein extract (10–100 µg) was subjected to western blot analysis with an appropriate antibody. The isolated protein was electrophoresed in a 10% or 15% polyacrylamide-SDS gel and transferred to a nitrocellulose membrane. After blocking for 1 h with 5% skim milk, the membrane was incubated with a primary antibody. After washing, the membrane was then incubated with horseradish peroxidase-conjugated goat anti-rabbit IgG or goat anti-mouse IgG antibody. Immunoreactive signals were detected with an enhanced chemiluminescence detection kit (Merck Millipore, Temecula, CA, USA) and analyzed using a Chemi-Doc System (Bio-Rad).
Statistical analysis
Groups were compared using chi-square tests for categorical variables or the Kruskal–Wallis test followed by post hoc Dunn’s multiple comparisons among three groups or the Mann–Whitney U test to compare two groups for continuous variables. To evaluate the effects of miRNA on the levels of target gene expression or APP protein expression and fragmentation, we used the Mann–Whitney U test or one-sample t tests to compare data from 4 or 6 independent experiments, respectively. Receiver-operating characteristic (ROC) curve analyses were performed to assess the diagnostic performance of various miRNAs to distinguish patients from normal subjects using Prism (v. 8.0; GraphPad Software, San Diego, CA, USA). We considered P < 0.05 as significant.