收藏切换
Plasma Extracellular Vesicle MicroRNA Analysis of Alzheimer’s Disease Reveals Dysfunction of a Neural Correlation Network
收藏切换
PDF
Yuzhe Sun1, 2, 3, *, Zhen Hefu2, 3, Benchao Li1, Wang Lifang2, 3, Song Zhijie2, 4, Li Zhou1, Yan Deng1, Liu Zhili2, 4, Jiahong Ding2, 3, Tao Li2, Wenwei Zhang2, 3, *, Nie Chao2, 3, *, Shuang Rong1, *
Research. Vol 6 Article ID 0114
Less
收藏切换
Research. Vol 6 Article ID 0114
• Research Article •
Plasma Extracellular Vesicle MicroRNA Analysis of Alzheimer’s Disease Reveals Dysfunction of a Neural Correlation Network
Full
Yuzhe Sun1, 2, 3, *, Zhen Hefu2, 3, Benchao Li1, Wang Lifang2, 3, Song Zhijie2, 4, Li Zhou1, Yan Deng1, Liu Zhili2, 4, Jiahong Ding2, 3, Tao Li2, Wenwei Zhang2, 3, *, Nie Chao2, 3, *, Shuang Rong1, *
Affiliations
  • 1Department of Nutrition and Food Hygiene, School of Public Health, Medical College, Wuhan University of Science and Technology, Wuhan, China.
  • 2BGI-Shenzhen, Shenzhen, China.
  • 3Shenzhen Key Laboratory of Neurogenomics, BGI-Shenzhen, Shenzhen 518120, China.
  • 4College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China.
About Author:

† These authors contributed equally to this work.

Published: 2023-04-20 doi: 10.34133/research.0114
Outline
收藏切换

Small extracellular vesicle (sEV) is an emerging source of potential biomarkers of Alzheimer's disease (AD), but the role of microRNAs (miRNAs) in sEV is not well understood. In this study, we conducted a comprehensive analysis of sEV-derived miRNAs in AD using small RNA sequencing and coexpression network analysis. We examined a total of 158 samples, including 48 from AD patients, 48 from patients with mild cognitive impairment (MCI), and 62 from healthy controls. We identified an miRNA network module (M1) that was strongly linked to neural function and showed the strongest association with AD diagnosis and cognitive impairment. The expression of miRNAs in the module was decreased in both AD and MCI patients compared to controls. Conservation analysis revealed that M1 was highly preserved in the healthy control group but dysfunctional in the AD and MCI groups, suggesting that changes in the expression of miRNAs in this module may be an early response to cognitive decline prior to the appearance of AD pathology. We further validated the expression levels of the hub miRNAs in M1 in an independent population. The functional enrichment analysis showed that 4 hub miRNAs might interact with a GDF11-centered network and play a critical role in the neuropathology of AD. In summary, our study provides new insights into the role of sEV-derived miRNAs in AD and suggests that M1 miRNAs may serve as potential biomarkers for the early diagnosis and monitoring of AD.

Yuzhe Sun, Zhen Hefu, Benchao Li, Wang Lifang, Song Zhijie, Li Zhou, Yan Deng, Liu Zhili, Jiahong Ding, Tao Li, Wenwei Zhang, Nie Chao, Shuang Rong. Plasma Extracellular Vesicle MicroRNA Analysis of Alzheimer’s Disease Reveals Dysfunction of a Neural Correlation Network[J]. Research, 2023 , 6 (4) : 0114 - . DOI: 10.34133/research.0114
Neurodegeneration-derived Alzheimer's disease (AD) is the second most prevalent affliction among elderly populations [1]. Due to the prolonged pathological progression of AD and its latency, the clinical presentation of AD is classified into 3 categories: cognitively normal, mild cognitive impairment (MCI), and full-blown AD [2]. MCI due to AD refers to an intermediate state between normal cognitive aging and symptomatic AD. It is characterized by memory loss, which is inconsistent with age but has not yet reached the clinical diagnosis standard of dementia [3–5]. Longitudinal cohort studies suggest that slowly progressing neuropathological changes, such as amyloid β (Aβ) deposition, begin in the early stages of MCI [6,7]. It is plausible that extending the early stages of MCI may delay the onset of AD pathology, underscoring the need to identify biomarkers of pathological and molecular changes during this stage [6]. However, early diagnosis of AD is challenging because neurodegenerative phenotypes are frequently accompanied by several cognitive disorders that are nearly indistinguishable from age-related cognitive decline [4,7]. Despite efforts to identify early-stage AD biomarkers, the most important clinical indicators remain Aβ and Tau [8,9].
Small extracellular vesicle (sEV), which contains proteins and RNAs that are transported between cells over distances that could trigger changes in gene expression and cellular function, exhibits a unique mode of cell communication [10–12]. As sEVs are capable of crossing the blood–brain barrier, they represent a valuable asset in the molecular diagnostics of neurological diseases [13]. In the central nervous system (CNS), 70% of microRNAs (miRNAs) are released from human brain cells and possibly regulate transcription of more than one-third of all genes [14,15]. Tissue-specific miRNAs can be isolated from enriched extracellular vesicles such as exosomes, the most studied subtype of EVs. In brain disorders, neuronal cells release sEVs that are enriched with miRNAs and other cellular components, which spread from the cerebrospinal fluid (CSF) into the peripheral blood system [16]. The presence of sEVs in the CSF, detectable in peripheral blood, makes them a promising pool of potential biomarker for neurodegenerative disease [17].
To date, several potential miRNA biomarkers for neurodegenerative disorders (NDs), including AD, Parkinson's disease, Huntington's disease, and other memory disorders, have been examined in blood, plasma, and serum by either quantitative polymerase chain reaction (qPCR) or small RNA (sRNA) sequencing, including miR-149-5p, miR-124-3p, miR-9-3p/5p, and miR-125b-5p [18–21]. However, the complexity of miRNA regulation presents challenges in fully understanding gene interaction networks. Weighted gene correlation network analysis (WGCNA) was developed to identify highly correlated expression patterns and calculate relationships between selected modules and external sample traits [22]. Gene expression-based network analysis has demonstrated its usefulness in the prognosis of various ND subtypes [23,24]. The identification of these disorders at an early stage through molecular characterization holds the potential to enhance the precision and prediction of diagnosis. Aberrantly expressed genes play a modular role in the progression of disease symptoms, and the deactivation and activation of coexpression networks may provide insight into the underlying mechanisms of NDs [23].
In this study, we aimed to explore the diagnostic potential of sEVs and their contained miRNA in the pathological classification of AD and MCI. To this end, we isolated sEVs from plasma of AD, MCI, and NC (normal control) and performed a comprehensive comparison of sEV miRNA expression between the 3 groups. Our results demonstrated the feasibility of using differentially expressed (DE) miRNAs as a diagnostic tool for differentiating between MCI and AD. Additionally, we performed WGCNA to uncover highly correlated patterns of gene expression and to assess the relationships between these patterns and the external traits of the samples. Our findings revealed a network that was strongly correlated with AD diagnosis and was conserved in NC but not in MCI and AD, implying its dysfunction in the latter 2 groups. Conversely, the majority of miRNA patterns were found to be similar between NC and MCI, suggesting that this network could serve as a biomarker of response to cognitive impairment. Tissue-specific analysis showed that most hub miRNAs in this network were expressed predominantly in the brain, providing evidence for the hypothesis that these miRNAs may originate from the CNS and serve as markers of posttranscriptional regulation in this organ. These results provide a basis for future studies aimed at uncovering the mechanisms underlying the progression of NDs and developing novel strategies for early diagnosis and intervention.
To investigate the miRNA expression profile of AD patients, we collected 158 plasma samples from 48 AD cases, 48 sporadic MCI cases, and 62 NC cases (Fig. 1A). The diagnosis of AD and MCI was made according to the clinical diagnostic criteria set forth by the National Institute on Aging and Alzheimer's Association (NIA-AA) and Petersen's criteria, respectively [3,25,26]. To evaluate cognitive function, patients were subjected to various cognitive tests, including mini-mental state examination (MMSE), clinical dementia rating (CDR), Hachinski ischemic score (HIS), and geriatric depression scale (GDS). Detailed demography of patients is shown in Table.
Demographic analysis showed that there were no significant differences in age or sex among the AD, MCI, and NC groups (Table). However, the MMSE scores were significantly lower in both the AD and MCI groups compared to the NC group, while the Montreal Cognitive Assessment (MoCA) scores were significantly lower in the MCI group compared to the NC group (Table and Table S1). These findings indicate a continuous decline in cognitive function across the 3 groups. Body mass index (BMI) was significantly lower in the AD group than in the MCI and NC groups. To gain a comprehensive understanding of the participants, disease history and lifestyle were collected via detailed initial health and demographic questionnaires. Education and marriage differed among 3 groups, which was consistent with the previous studies [27]. Furthermore, our findings suggested that preexisting conditions of hypertension and stroke might predispose individuals to developing AD. Lifestyle factors such as smoking, tea drinking, entertainment, reading frequency, exercise, and communication with neighbors and children varied significantly among the 3 groups (Table).
We isolated sEVs from plasma samples of AD, MCI, and NC groups, using a previously described precipitation method [17]. Transmission electron microscopy and Western blot were performed to confirm successfully isolated sEVs. Additionally, Brownian motion-based nanoparticle tracking analysis (NTA) was applied to measure distribution of nanoparticle size in the 3 groups (Fig. S1).
Following the extraction and sequencing of sRNA, we identified 2,026 miRNAs from 158 sRNA libraries (Fig. 1A and Table S2). Using R packages DESeq2 and EdgeR, 11 and 16 miRNAs were found to be significantly down-regulated in AD and MCI, respectively, compared to NC. MiR-6891-5p and miR-7975 were significantly up-regulated in AD compared to NC, and no up-regulated miRNA was observed in MCI, pointing toward a general decline of EV-derived miRNAs in both AD and MCI (Fig. 1B and Table S3). The down-regulation of miR-9-3p and miR-9-5p was also observed, which was consistent with previous research in brain tissues [28,29]. Interestingly, the expression levels of 14 miRNAs were significantly higher and 2 miRNAs were significantly lower in AD as compared to MCI (Table S3).
Target prediction and gene enrichment analysis showed that the DE miRNAs in AD were enriched in pathways related to axonogenesis, forebrain development, axon guidance, and Rap1 signaling, consistent with the widely observed dysfunction of synapses in AD patients (Fig. 1C) [30]. DE miRNAs in MCI were enriched in pathways related to reproductive system/structure development, regulation of cell morphogenesis, histone modification, and covalent chromatin modification (Fig. 1C). We also identified 8 miRNAs commonly suppressed in both AD and MCI, enriched in axon- and synapse-related functions (Fig. 1D).
Our study further examined the predictive powers of DE miRNAs as biomarkers, constructing 7-miRNA models of AD-NC and 3-miRNA models of MCI-NC, which showed area under the curve (AUC) values of 0.80 (P < 0.001) and 0.71 (P < 0.001), respectively. Besides, we downloaded sEV miRNA profiles of previous studies, and using the same miRNAs and criteria, AUC of the AD-NC model was 0.61 and 0.59 (Fig. 1E and F) [31,32]. However, when compared to a previous study (AUC = 0.91, P < 0.001) that combined plasma P-tau217, memory, executive function, and APOE, our results indicated that plasma sEV miRNA biomarkers were not more advantageous in AD diagnosis than traditional biomarkers [33].
To further investigate the biological function of sEV miRNA, we used the WGCNA algorithm to generate a coexpression network from the 1,000 most abundant miRNAs. The resulting network consisted of 13 miRNA coexpression modules of similar expression patterns across the 158 cases analyzed (Fig. S2). The modules ranged in size from 202 miRNAs (M12, turquoise) to 12 miRNAs (M4, salmon). GO analysis of miRNA module members revealed putative functions of 13 modules, out of which 11 modules had a clear ontology, encompassing a diverse mix of biological processes and functions (Fig. 2A and Table S4). To assess whether a given coexpression module was related to either AD or MCI phenotypes, we correlated module eigengene values (the first principal component of the module miRNA expression level) to hallmarks of demography, personal health, and habits. We also correlated module eigengene values to cognitive function as assessed by the MMSE to capture module–AD/MCI relationships (Fig. 2A). We observed 3 modules that were significantly correlated with AD diagnosis and MMSE: M1 neural function, M3 transcription repressor, and M7 guanosine triphosphatase (GTPase) binding. The M1 neural function module exhibited the strongest AD trait correlations (AD diagnosis, P = 9 × 10−5; MMSE, P = 0.004), indicating that expression levels of M1 miRNAs might regulate AD neuropathology. Besides, M1 was correlated to education, smoking, and reading, which are important factors in AD development [34–36]. The M3 transcription repressor and M7 GTPase binding were correlated with education and numerous personal habits. Education level and reading behavior were associated with more than half of the modules (Fig. 2A). These findings suggested that biological regulation of M1, M3, and M7 might be altered during disease progression. It is also worth noting that 4 modules (M2 SMAD binding, M5 calcium transporting, M6 transcription coregulator, and M12 proximal promoter) were correlated with MMSE (P < 0.05), indicating that the regulatory functions of these modules are possibly affected by cognition decline.
To investigate the relationship between diagnostic classification and coexpression modules, we compared the module eigengene values by case status. Five modules exhibited significant differences among 3 diagnostic categories, including M1 neural function, M3 transcription repressor, M5 calcium transporting, M7 GTPase binding, and M8 regulation by MECP2 modules (Fig. 2B). Compared to NC, eigengene values of AD and MCI in M1 were significantly lower (P = 3.43 × 10−5), eigengene values of M3 and M5 decreased in AD but increased in MCI (P = 0.00028 and 0.0392, respectively), while eigengene values of M7 and M8 increased in AD but decreased in MCI (P = 0.015 and 0.0212, respectively) (Fig. 2B). An intramodular analysis of gene significance and module membership in the 5 modules was performed. AD diagnosis and module membership of M1, M3, and M8 exhibit significant correlation, and M1 was the most relevant module with AD (Fig. 2C). In summary, module eigengenes in AD when compared to NC were conversely altered in MCI, except M1. This suggested that M1 was the module especially related to AD, and miRNAs in M1 were potential biomarkers for AD staging.
The coexpression networks including M1, M3, M5, M7, and M8 contained 105, 25, 86, 20, and 25 miRNAs, respectively (Table S5). In order to determine the key miRNAs within each network, we calculated the membership of each miRNA (kME) and defined the miRNAs with the top 25% kME values in each coexpression module as “hub miRNAs.” Subsequently, we analyzed the beta diversity (sample dissimilarity) based on weighted distances, which was calculated using expression levels of hub miRNAs by case status. Unlike the previously discussed eigengene values, our results showed a significant decrease in diversity of AD and MCI compared to NC in modules M1, M3, M5, and M8. These findings suggest that the expression levels of hub miRNAs in the respective coexpression modules may have a shared pattern of change in the presence of AD or MCI (Fig. S3).
In our study, we identified 26 hub miRNAs in the M1 network and discovered that 13 of these miRNAs were previously reported as AD markers (Table S6). Then, we queried the M1 hub miRNAs in the Human miRNA tissue atlas (https://ccb-web.cs.uni-saarland.de/tissueatlas/) and found that 24 of 26 M1 hub miRNAs were specifically expressed in brain tissues (Table S6), indicating a potential role for these miRNAs in nervous system function [37].
To determine the status of the M1 network in AD, MCI, and NC, we evaluated the preservation of the model for each group (Fig. 3). We found that 12 of 13 network modules were preserved in AD and NC, and 10 of 13 modules were preserved in MCI. The preservation Zsummary values of 13 modules in MCI and NC were relatively close but changed greatly in AD, indicating that the miRNA coexpression networks did not change greatly in MCI, but rather in AD. Remarkably, M1 neural function was highly preserved in NC (Zsummary = 12.9) but not preserved in MCI (Zsummary = 0.2) and AD (Zsummary = 1.3), suggesting that the dysfunction of M1 miRNA network may occur at the MCI stage prior to the onset of AD. Additionally, M12 and M13 proximal promoters were highly preserved (Zsummary > 10) in all 3 groups, potentially regulating housekeeping functions that do not change during cognitive decline (Fig. 3).
We examined the relationship between M1 hub miRNAs and 38 plasma proteins that were reported associated with incident dementia [38]. By using a comprehensive approach that included target predictions from 8 different databases, we obtained 49 pairs of correspondences in at least 2 databases (Table S7). After removing the proteins that were only linked to one miRNA, the association of 11 genes with 10 M1 hub miRNAs was retained. Our analysis revealed that hsa-miR-125b-5p, hsa-miR-9-3p, hsa-miR-9-5p, and hsa-miR-218-5p emerged as key miRNAs, with growth differentiation factor 11 (GDF11) being the most highly connected gene (Fig. 4A). We also observed that the immunologically relevant cellular adhesion protein, Sushi, von Willebrand factor type A, epidermal growth factor, and pentraxin domain-containing protein 1 (SVEP1), which has been implicated in brain atrophy and Alzheimer's pathology, was targeted by hsa-miR-338-3p and involved in the transforming growth factor-β (TGF-β) signaling pathway. These findings suggest that the M1 miRNAs may contribute to Alzheimer's pathology through their involvement in a GDF11-centered network (Fig. 4B).
To validate the sRNA sequencing analysis, the abundance of selected miRNAs (miR-9-5p, miR338-3p, and miR-125b-5p) was validated by quantitative reverse transcription PCR (qRT-PCR) assay (Fig. 5A). We used 60 samples in the experiment, 20 AD, 20 MCI, and 20 NC samples. The qRT-PCR validation results of 3 miRNA expression levels corroborated those of the sequencing analyses (Fig. 5B). To reinforce the validity of the findings, the miRNA expression levels were also measured in a separate set of 40 samples (20 MCI and 20 NC). There were no significant differences in age, BMI, education, marriage, hypertension, diabetes, smoking, or alcohol consumption between MCI and NC samples. The MMSE scores of the MCI group were significantly lower than the NC group (Table S8). The MCI group had lower levels of miR-9-5p and miR-338-3p expressions than the NC group, consistent with our qRT-PCR results (Fig. 5C). The MCI group had insignificantly lower expression levels of miR-125b-5p than the NC group (Fig. 5C). These results provide additional support for the reliability of the sequencing analysis.
Aβ and Tau proteins have been widely studied as biomarkers for AD, and they have shown potential for the early diagnosis and monitoring of the disease. Expression of miRNA and accumulation of Aβ and Tau are different components of the biology of AD, and their potential for clinical diagnosis is also different. Growing evidence suggests that miRNAs contribute to the onset and development of NDs through synaptic plasticity and various signaling pathways [39–42]. MiRNAs are enveloped in sEVs and released from cells to the environment or circulation system, which can be recruited by cells at a distance or in another tissue where they can affect gene expression. However, enrichment of miRNA varies widely in plasma, blood immune cells, CSF, and brain tissue, and the biological functions of miRNA in sEVs of peripheral blood remain elusive [43,44].
In this study, we measured expression levels of sRNA from plasma sEVs in AD, MCI, and NC groups. The miRNA coexpression networks were constructed, and M1 neural function was presumed dysfunctional in AD and MCI (Fig. 3). Of the M1 hub miRNAs, both hsa-miR-9-5p and hsa-miR-9-3p are proved to be enriched in the CNS and regulate essential neuronal processes such as neuronal differentiation and synaptic plasticity [45–47]. In particular, as the potential target of miR-9-5p and miR-9-3p, GDF11 is a member of the TGF-β superfamily and is reportedly beneficial in preventing age-related degeneration in CNS, enhancing cognitive function, and stimulating tissue regeneration [48,49]. Besides hsa-miR-9-5p and has-miR-9-3p, hsa-miR-125b-5p negatively regulates GDF11 expression, which triggers canonical signal transduction mediated by R-SMAD proteins (Fig. 4B). A previous study has indicated that a reduction in hsa-miR-125b-5p levels leads to an increase in GDF11 expression and attenuates neuronal injury [50]. This supported the hypothesis that miRNAs in M1 may be involved in the underlying mechanisms of AD.
In the intricate dance of neural biology, Beta-secretase 1 (BACE1) steps to the forefront as the leading performer, responsible for the creation of Aβ peptides [51]. Yet, recent studies have uncovered a perturbation in this delicate balance, and the decreased expression levels of hsa-miR-9-5p, hsa-miR-9-3p, and hsa-miR-338-3p lead to Aβ42 increase by targeting BACE1 and exacerbating the disease [52,53]. Besides BACE1, hsa-miR-338-3p also targets growth hormone (GH) [54], whose receptor (GHR) was targeted by hsa-miR-9-3p and hsa-miR-9-5p, indicating that M1 miRNAs might be involved in the GH-related circadian rhythm (Fig. 4). SVEP1 is a newly reported plasma biomarker, high levels of which present a causal relationship with AD [38]. Combining hsa-miR-338-3p’s regulatory relation with BACE1 and GH and its putative target SVEP1, hsa-miR-338-3p plays an essential role in the development of neuropathology in AD.
This study has some inherent limitations. First, the number of AD samples is relatively small owing to its low incidence. Another reason is that the period of AD diagnosis is usually more extended, and AD lacks effective treatment methods as well, which will affect patients' enthusiasm for participating in this study. Second, this community-based case–control study precluded us from making any temporal association between the miRNA expression and the AD diagnosis or cognitive decline. Next, the significance of M1 comes from the relationship between coexpression indices of interaction profiles within networks and clinical data of the AD patients. This relationship mainly indicates the biological mechanism of the network and the disease, rather than diagnostic power. Although we have validated our findings in another independent MCI and NC population, we cannot conduct an effective assessment of AD early diagnostic biomarkers from miRNA expression, and we would like to measure the miRNA expression by constructing a longitudinal cohort in the future.
In summary, we found a completely altered network preservation pattern in the AD group, corroborating the clinical pathology of AD cognitive impairment and MMSE score (Table). Our findings, backed by clinical pathology and cognitive impairment scores, paint a picture of a disrupted network and a dysfunction in the M1 neural function. The discovery of eight miRNAs with decreased expression, all found within the M1 module, and confirmed through tissue-specific analysis, provides further evidence of a breakdown in the delicate balance of posttranscriptional regulation. The M1 module could emerge as a potential source of biomarkers for the early detection of AD. Our investigation represents a critical step toward understanding the intricacies of the fundamental mechanisms of AD.
In this study, 158 individuals were enrolled, including 48 AD patients, 48 MCI patients, and 62 NC volunteers. Participants were from communities of Wuhan, Huangshi, and Jingmen Cities in Hubei province, China. Cases were frequency-matched by age and sex. Cognitive tests performed were MMSE, MoCA, CDR, HIS, and GDS. The diagnosis of AD was based on the criteria of the NIA-AA [25]. Briefly, the AD was diagnosed in clinical settings with the following criteria: (a) meeting the criteria for dementia, (b) excluding vascular dementia by distinguishing cerebrovascular diseases by computed tomography scans of the brain, (c) excluding the patients with prominent features of Lewy body dementia or frontotemporal dementia, and (d) excluding the active neurological diseases or other medication-induced cognitive disorders. The diagnosis of MCI was based on the Petersen's criteria from previous studies [3,26]. All participants' enrollment was approved by the Institutional Review Board of BGI (BG-IRB 20149) and the Medical Ethics Committee of Wuhan University of Science and Technology (No. 049). The standardized questionnaire was used to obtain demographic characteristics, history of diseases, and lifestyle factors. Written consent was obtained from all subjects.
To isolate sEVs from plasma, the commercial kit System Biosciences ExoQuick Kit was used in accordance with the manufacturer's instructions [12]. Briefly, plasma samples were first thawed and mixed by vortexing. Subsequently, 50 μl of plasma sample was subjected to centrifugation at 3,000g for 15 min to separate cells and cellular debris from the supernatant. The supernatant was then transferred to a sterile vessel, mixed with 63 μl of ExoQuick Precipitation Solution (63 μl), and refrigerated for 30 min. After refrigeration, the mixture was centrifuged at 1,500g for 30 min followed by a second centrifugation at 1,500g for 5 min to remove any residual ExoQuick solution. The resulting pellet was then resuspended in 25 μl of nuclease-free water.
sRNA was extracted from the resuspended sEV solution with the miRNeasy Serum/Plasma Kit (catalog no. 217084, QIAGEN, Germany) as per the manufacturer's instructions. Using Agilent's small RNA kit, sRNA samples were subsequently analyzed on a Bioanalyzer 2100 system (Agilent, Santa Clara, USA) to determine RNA integrity numbers (RINs). The concentration of sRNA was measured using NanoDrop 2000 (Thermo Scientific). sRNA sequencing libraries were prepared from total sEV sRNA with the MGIEasy Small RNA Library Prep Kit (MGI, Shenzhen, China). Pooled libraries were loaded and sequenced on a BGISEQ-500 platform (BGI), and more than 2.0 GB of reads were obtained from each library. Generated sequenced reads were deposited into the CNGB Sequence Archive (CNSA) of China National GeneBank DataBase (CNGBdb) under accession number CNP0001975 [55,56].
We generated an average of 2.2 Gb raw data per library from sRNA sequencing. Raw reads were processed by trimming adaptor sequences and then culling low-quality reads using SOAPnuke v1.5.0 (with the setting -Q 2 -q -c 0). High-quality reads that ranged from 70 to 500 Mb per sample were mapped against the reference genome (hg19) using Bowtie 2 (with the setting -q -L 16 -p 6 --phred64 --rdg 1,10 --rfg 1,10) [57] with allowance for only one mismatch. The matched reads were aligned to mature miRNAs in miRbase version 20 [58], and miRNA counts were calculated using perfectly matched reads. The remaining reads were mapped to the Rfam database [59] to predict novel miRNAs by identification and removal of protein-coding genes, transfer RNA (tRNA), and ribosomal RNA (rRNA). We used miRDeep2 to predict the novel miRNA from the retained reads [60].
Analysis of DE miRNA was performed using EdgeR and DESeq2 [61,62]. The thresholds for significant miRNA expression changes were >±1-fold (log2) and false discovery rate (FDR) ≤ 0.05, and P values were corrected to FDR using the Benjamin and Hochberg (BH) method. R package MultimiR was used to predict the targets of miRNA [63]. Gene ontology (GO) analysis was performed on putative target genes using clusterProfiler [64]. DE miRNAs were selected to calculate the diagnostic power. R package Caret was used to combine the predictive powers of DE miRNAs. Logistic regression with internal 10-fold cross-validation was used to develop the model. The AUC was calculated to assess the model's power by R package pROC.
The coexpression network analysis was performed using the R package WGCNA [22]. To construct a weighted coexpression network, expression levels of 158 sRNA-seq samples were used, including AD, MCI, and NC groups. The WGCNA network function was used with the following settings: soft threshold power β = 5, deepSplit = 4, minimum module size of 12, merge cut height of 0.07. Pearson correlations between every miRNA and module eigengene were performed. After the initial network construction, 13 modules consisting of 12 to 202 miRNAs were detected. The first principal component (eigengene value) was calculated and considered as representative of each module. We discarded the gray module, which cannot be merged into any other modules, and then correlated the eigengene values of modules with phenotype traits shared by AD, MCI, and NC samples using the function corPvalueStudent() of the Pearson method. MiRNAs with the highest eigengene values were selected and predicted the targets, and functional enrichment analysis was performed using Metascape. Based on the GO annotation of miRNA targets, we then summarized the potential function of each network. Using the above networks as the template, Zsummary composite preservation scores were calculated for each target group, with 500 permutations. The modulePreservation() function in the WGCNA package was used.
In this study, 60 samples were selected out of 158 for qRT-PCR validation (20 samples from each of the 3 groups). The sEV miRNA was converted into cDNA using the TaqMan Advanced miRNA cDNA Synthesis Kit (Applied Biosystems, #A28007), and qRT-PCR was performed on a StepOnePlus Real-Time PCR System using miRNA assays in a 96-well format (TaqMan microRNA assays, 20×, Applied Biosystems, #A25576). Has-miR-186-5p was used as an endogenous control. Data were tested for normality using the Shapiro test. Additionally, expression levels of miR-9-5p, miR338-3p, and miR-125b-5p were assessed in an independent set of samples composed of 20 MCI and 20 NC samples.
1.
Schneider JA, Wilson RS, Bienias JL, Evans DA, Bennett DA. Cerebral infarctions and the likelihood of dementia from Alzheimer disease pathology. Neurology. 2004;62(7):1148–1155.
2.
Villemagne VL, Burnham S, Bourgeat P, Brown B, Ellis KA, Salvado O, Szoeke C, Macaulay SL, Martins R, Maruff P, et al.Amyloid β deposition, neurodegeneration, and cognitive decline in sporadic Alzheimer's disease: A prospective cohort study. Lancet Neurol. 2013;12(4):357–367.
3.
Petersen RC. Mild cognitive impairment as a diagnostic entity. J Intern Med. 2004;256(3):183–194.
4.
Langa KM, Levine DA. The diagnosis and management of mild cognitive impairment: A clinical review. JAMA. 2014;312(23):2551–2561.
5.
Albert MS, DeKosky ST, Dickson D, Dubois B, Feldman HH, Fox NC, Gamst A, Holtzman DM, Jagust WJ, Petersen RC, et al.The diagnosis of mild cognitive impairment due to Alzheimer's disease: Recommendations from the National Institute on Aging-Alzheimer's Association workgroups on diagnostic guidelines for Alzheimer's disease. Alzheimers Dement. 2011;7(3):270–279.
6.
Winblad B, Palmer K, Kivipelto M, Jelic V, Fratiglioni L, Wahlund L-O, Nordberg A, Bäckman L, Albert M, Almkvist O, et al.Mild cognitive impairment—Beyond controversies, towards a consensus: Report of the international working group on mild cognitive impairment. J Intern Med. 2004;256(3):240–246.
7.
Sims R, Hill M, Williams J. The multiplex model of the genetics of Alzheimer's disease. Nat Neurosci. 2020;23(3):311–322.
8.
Blennow K, Zetterberg H. Biomarkers for Alzheimer's disease: Current status and prospects for the future. J Intern Med. 2018;284(6):643–663.
9.
Masters CL, Bateman R, Blennow K, Rowe CC, Sperling RA, Cummings JL. Alzheimer's disease. Nat Rev Dis Primers. 2015;1:15056.
10.
Thomou T, Mori MA, Dreyfuss JM, Konishi M, Sakaguchi M, Wolfrum C, Rao TN, Winnay JN, Garcia-Martin R, Grinspoon SK, et al.Adipose-derived circulating miRNAs regulate gene expression in other tissues. Nature. 2017;542(7642):450–455.
11.
Crewe C, Joffin N, Rutkowski JM, Kim M, Zhang F, Towler DA, Gordillo R, Scherer PE. An endothelial-to-adipocyte extracellular vesicle axis governed by metabolic state. Cell. 2018;175(3):695–708.e13.
12.
Saugstad JA, Lusardi TA, van Keuren-Jensen KR, Phillips JI, Lind B, Harrington CA, McFarland TJ, Courtright AL, Reiman RA, Yeri AS, et al.Analysis of extracellular RNA in cerebrospinal fluid. J Extracell Vesicles. 2017;6(1):Article 1317577.
13.
Hill AF. Extracellular vesicles and neurodegenerative diseases. J Neurosci. 2019;39(47):9269–9273.
14.
Nowak JS, Michlewski G. miRNAs in development and pathogenesis of the nervous system. Biochem Soc Trans. 2013;41(4):815–820.
15.
Kosik KS. The neuronal microRNA system. Nat Rev Neurosci. 2006;7(12):911–920.
16.
Cheng L, Vella LJ, Barnham KJ, McLean C, Masters CL, Hill AF. Small RNA fingerprinting of Alzheimer's disease frontal cortex extracellular vesicles and their comparison with peripheral extracellular vesicles. J Extracell Vesicles. 2020;9(1):Article 1766822.
17.
Jia L, Zhu M, Kong C, Pang Y, Zhang H, Qiu Q, Wei C, Tang Y, Wang Q, Li Y, et al.Blood neuro-exosomal synaptic proteins predict Alzheimer's disease at the asymptomatic stage. Alzheimers Dement. 2021;17(1):49–60.
18.
Chen F, Chen H, Jia Y, Lu H, Tan Q, Zhou X. miR-149-5p inhibition reduces Alzheimer's disease β-amyloid generation in 293/APPsw cells by upregulating H4K16ac via KAT8. Exp Ther Med. 2020;20(5):88.
19.
Zhou Y, Deng J, Chu X, Zhao Y, Guo Y. Role of post-transcriptional control of calpain by miR-124-3p in the development of Alzheimer's disease. J Alzheimers Dis. 2019;67(2):571–581.
20.
Coolen M, Katz S, Bally-Cuif L. miR-9: A versatile regulator of neurogenesis. Front Cell Neurosci. 2013;7:220.
21.
Jin Y, Tu Q, Liu M. MicroRNA125b regulates Alzheimer's disease through SphK1 regulation. Mol Med Rep. 2018;18(2):2373–2380.
22.
Langfelder P, Horvath S. WGCNA: An R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.
23.
Johnson ECB, Dammer EB, Duong DM, Ping L, Zhou M, Yin L, Higginbotham LA, Guajardo A, White B, Troncoso JC, et al.Large-scale proteomic analysis of Alzheimer's disease brain and cerebrospinal fluid reveals early changes in energy metabolism associated with microglia and astrocyte activation. Nat Med. 2020;26(5):769–780.
24.
Kakati T, Bhattacharyya DK, Barah P, Kalita JK. Comparison of methods for differential co-expression analysis for disease biomarker prediction. Comput Biol Med. 2019;113:Article 103380.
25.
McKhann GM, Knopman DS, Chertkow H, Hyman BT, Jack CR Jr, Kawas CH, Klunk WE, Koroshetz WJ, Manly JJ, Mayeux R, et al.The diagnosis of dementia due to Alzheimer's disease: Recommendations from the National Institute on Aging-Alzheimer's Association workgroups on diagnostic guidelines for Alzheimer's disease. Alzheimers Dement. 2011;7(3):263–269.
26.
Lu J, Li D, Li F, Zhou A, Wang F, Zuo X, Jia XF, Song H, Jia J. Montreal cognitive assessment in detecting cognitive impairment in Chinese elderly individuals: A population-based study. J Geriatr Psychiatry Neurol. 2011;24(4):184–190.
27.
Xu W, Tan L, Wang HF, Jiang T, Tan MS, Tan L, Zhao QF, Li JQ, Wang J, Yu JT. Meta-analysis of modifiable risk factors for Alzheimer's disease. J Neurol Neurosurg Psychiatry. 2015;86(12):1299–1306.
28.
Geekiyanage H, Jicha GA, Nelson PT, Chan C. Blood serum miRNA: Non-invasive biomarkers for Alzheimer's disease. Exp Neurol. 2012;235(2):491–496.
29.
Geekiyanage H, Chan C. MicroRNA-137/181c regulates serine palmitoyltransferase and in turn amyloid β, novel targets in sporadic Alzheimer's disease. J Neurosci. 2011;31(41):14820–14830.
30.
Mota SI, Ferreira IL, Rego AC. Dysfunctional synapse in Alzheimer's disease—A focus on NMDA receptors. Neuropharmacology. 2014;76(Pt A):16–26.
31.
Cheng L, Doecke JD, Sharples RA, Villemagne VL, Fowler CJ, Rembach A, Martins RN, Rowe CC, Macaulay SL, Masters CL, et al.Prognostic serum miRNA biomarkers associated with Alzheimer's disease shows concordance with neuropsychological and neuroimaging assessment. Mol Psychiatry. 2015;20(10):1188–1196.
32.
Lugli G, Cohen AM, Bennett DA, Shah RC, Fields CJ, Hernandez AG, Smalheiser NR. Plasma Exosomal miRNAs in persons with and without Alzheimer disease: Altered expression and prospects for biomarkers. PLoS One. 2015;10(10):Article e0139233.
33.
Palmqvist S, Tideman P, Cullen N, Zetterberg H, Blennow K, the Alzheimer’s Disease Neuroimaging Initiative, Dage JL, Stomrud E, Janelidze S, Mattsson-Carlgren N, et al.Prediction of future Alzheimer's disease dementia using plasma phospho-tau combined with other accessible measures. Nat Med. 2021;27(6):1034–1042.
34.
Bennett DA, Wilson RS, Schneider JA, Evans DA, Mendes de Leon CF, Arnold SE, Barnes LL, Bienias JL. Education modifies the relation of AD pathology to level of cognitive function in older persons. Neurology. 2003;60(12):1909–1915.
35.
Aggarwal NT, Bienias JL, Bennett DA, Wilson RS, Morris MC, Schneider JA, Shah RC, Evans DA. The relation of cigarette smoking to incident Alzheimer's disease in a biracial urban community population. Neuroepidemiology. 2006;26(3):140–146.
36.
Friedman RB, Ferguson S, Robinson S, Sunderland T. Dissociation of mechanisms of reading in Alzheimer's disease. Brain Lang. 1992;43(3):400–413.
37.
Ludwig N, Leidinger P, Becker K, Backes C, Fehlmann T, Pallasch C, Rheinheimer S, Meder B, Stähler C, Meese E, et al.Distribution of miRNA expression across human tissues. Nucleic Acids Res. 2016;44(8):3865–3877.
38.
Walker KA, Chen J, Zhang J, Fornage M, Yang Y, Zhou L, Grams ME, Tin A, Daya N, Hoogeveen RC, et al.Large-scale plasma proteomic analysis identifies proteins and pathways associated with dementia risk. Nat Aging. 2021;1(5):473–489.
39.
Qureshi IA, Mehler MF. Non-coding RNA networks underlying cognitive disorders across the lifespan. Trends Mol Med. 2011;17(6):337–346.
40.
Colantuoni C, Lipska BK, Ye T, Hyde TM, Tao R, Leek JT, Colantuoni EA, Elkahloun AG, Herman MM, Weinberger DR, et al.Temporal dynamics and genetic control of transcription in the human prefrontal cortex. Nature. 2011;478(7370):519–523.
41.
Gao J, Wang WY, Mao YW, Gräff J, Guan JS, Pan L, Mak G, Kim D, Su SC, Tsai LH. A novel pathway regulates memory and plasticity via SIRT1 and miR-134. Nature. 2010;466(7310):1105–1109.
42.
Peng D, Wang Y, Xiao Y, Peng M, Mai W, Hu B, Jia Y, Chen H, Yang Y, Xiang Q, et al.Extracellular vesicles derived from astrocyte-treated with haFGF(14-154) attenuate Alzheimer phenotype in AD mice. Theranostics. 2022;12(8):3862–3881.
43.
Takousis P, Sadlon A, Schulz J, Wohlers I, Dobricic V, Middleton L, Lill CM, Perneczky R, Bertram L. Differential expression of microRNAs in Alzheimer's disease brain, blood, and cerebrospinal fluid. Alzheimers Dement. 2019;15(11):1468–1477.
44.
Yoon S, Kim SE, Ko Y, Jeong GH, Lee KH, Lee J, Solmi M, Jacob L, Smith L, Stickley A, et al.Differential expression of MicroRNAs in Alzheimer's disease: A systematic review and meta-analysis. Mol Psychiatry. 2022;27(5):2405–2413.
45.
Madelaine R, Sloan SA, Huber N, Notwell JH, Leung LC, Skariah G, Halluin C, Paşca SP, Bejerano G, Krasnow MA, et al.MicroRNA-9 couples brain neurogenesis and angiogenesis. Cell Rep. 2017;20(7):1533–1542.
46.
Giusti SA, Vogl AM, Brockmann MM, Vercelli CA, Rein ML, Trümbach D, Wurst W, Cazalla D, Stein V, Deussing JM, et al.MicroRNA-9 controls dendritic development by targeting REST. Elife. 2014;3:Article e02755.
47.
Yuva-Aydemir Y, Simkin A, Gascon E, Gao FB. MicroRNA-9: Functional evolution of a conserved small regulatory RNA. RNA Biol. 2011;8(4):557–564.
48.
Katsimpardi L, Litterman NK, Schein PA, Miller CM, Loffredo FS, Wojtkiewicz GR, Chen JW, Lee RT, Wagers AJ, Rubin LL. Vascular and neurogenic rejuvenation of the aging mouse brain by young systemic factors. Science. 2014;344(6184):630–634.
49.
Sinha M, Jang YC, Oh J, Khong D, Wu EY, Manohar R, Miller C, Regalado SG, Loffredo FS, Pancoast JR, et al.Restoring systemic GDF11 levels reverses age-related dysfunction in mouse skeletal muscle. Science. 2014;344(6184):649–652.
50.
Chen W, Wang H, Feng J, Chen L. Overexpression of circRNA circUCK2 attenuates cell apoptosis in cerebral ischemia-reperfusion injury via miR-125b-5p/GDF11 signaling. Mol Ther Nucleic Acids. 2020;22:673–683.
51.
Yuan P, Zhang M, Tong L, Morse TM, McDougal RA, Ding H, Chan D, Cai Y, Grutzendler J. PLD3 affects axonal spheroids and network defects in Alzheimer's disease. Nature. 2022;612(7939):328–337.
52.
Xie H, Zhao Y, Zhou Y, Liu L, Liu Y, Wang D, Zhang S, Yang M. MiR-9 regulates the expression of BACE1 in dementia induced by chronic brain hypoperfusion in rats. Cell Physiol Biochem. 2017;42(3):1213–1226.
53.
Qian Q, Zhang J, He FP, Bao WX, Zheng TT, Zhou DM, Pan HY, Zhang H, Zhang XQ, He X, et al.Down-regulated expression of microRNA-338-5p contributes to neuropathology in Alzheimer's disease. FASEB J. 2019;33(3):4404–4417.
54.
Lee YJ, Cho JM, Moon JH, Ku CR, Kim J, Kim SH, Lee EJ. Increased miR-338-3p expression correlates with invasiveness of GH-producing pituitary adenomas. Endocrine. 2017;58(1):184–189.
55.
Guo X, Chen F, Gao F, Li L, Liu K, You L, Hua C, Yang F, Liu W, Peng C, et al.CNSA: A data repository for archiving omics data. Database (Oxford). 2020;2020:Article baaa055.
56.
Chen FZ, You LJ, Yang F, Wang LN, Guo XQ, Gao F, Hua C, Tan C, Fang L, Shan RQ, et al.CNGBdb: China National GeneBank DataBase. Yi Chuan. 2020;42(8):799–809.
57.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–359.
58.
Kozomara A, Griffiths-Jones S. miRBase: Annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014;42(Database issue):D68–D73.
59.
Burge SW, Daub J, Eberhardt R, Tate J, Barquist L, Nawrocki EP, Eddy SR, Gardner PP, Bateman A. Rfam 11.0: 10 years of RNA families. Nucleic Acids Res. 2013;41(Database issue):D226–D232.
60.
Friedländer MR, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012;40(1):37–52.
61.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–140.
62.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
63.
Ru Y, Kechris KJ, Tabakoff B, Hoffman P, Radcliffe RA, Bowler R, Mahaffey S, Rossi S, Calin GA, Bemis L, et al.The multiMiR R package and database: Integration of microRNA-target interactions along with their disease and drug associations. Nucleic Acids Res. 2014;42(17):e133.
64.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: An R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–287.
Year 2023 volume 6 Issue 4
PDF
92
49
Cite this Article
BibTeX
Article Info
doi: 10.34133/research.0114
  • Receive Date:2022-12-31
  • Online Date:2026-10-28
  • Published:2023-04-20
Article Data
Affiliations
History
  • Received:2022-12-31
  • Accepted:2023-03-20
Affiliations
    1Department of Nutrition and Food Hygiene, School of Public Health, Medical College, Wuhan University of Science and Technology, Wuhan, China.
    2BGI-Shenzhen, Shenzhen, China.
    3Shenzhen Key Laboratory of Neurogenomics, BGI-Shenzhen, Shenzhen 518120, China.
    4College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China.

Corresponding:

* Address correspondence to: (Y.S.)
References
Share
https://castjournals.cast.org.cn/joweb/research/EN/10.34133/research.0114
Share to
QR

Scan QR to access full text

Cite this article
BibTeX
Citations
表12种不同金属材料的力学参数
科
Family
属数
Number of
genus
种数
Number of
species
占总种数比例
Percentage of
total species (%)
属
Genus
种数
Number of
species
占总种数比例
Percentage of total
species (%)
鹅膏菌科Amanitaceae 2 11 5.26 鹅膏菌属 Amanita 10 4.78
小菇科 Mycenaceae 2 12 5.74 丝盖伞属 Inocybe 5 2.39
多孔菌科 Polyporaceae 8 14 6.70 蜡蘑属 Laccaria 5 2.39
红菇科 Russulaceae 3 23 11.00 小皮伞属 Marasmius 6 2.87
小菇属 Mycena 11 5.26
光柄菇属 Pluteus 5 2.39
红菇属 Russula 17 8.13
栓菌属 Trametes 5 2.39
关闭全屏
  • BibTeX
  • EndNote
  • RefWorks
  • TxT