ORIGINAL RESEARCH article

Front. Genet., 20 May 2022

Sec. Genetics of Common and Rare Diseases

Volume 13 - 2022 | https://doi.org/10.3389/fgene.2022.816919

MicroRNA Expression in Circulating Leukocytes and Bioinformatic Analysis of Patients With Moyamoya Disease

  • 1. Department of Neurology, Beijing Tiantan Hospital, Capital Medical University, Beijing, China

  • 2. China National Clinical Research Center for Neurological Diseases, Beijing, China

  • 3. Center of Stroke, Beijing Institute for Brain Disorders, Beijing, China

  • 4. Department of Neurosurgery, Beijing Tiantan Hospital, Capital Medical University, Beijing, China

  • 5. Research Unit of Artificial Intelligence in Cerebrovascular Disease, Chinese Academy of Medical Sciences, Beijing, China

Abstract

Objective: MicroRNAs (miRNAs) in exosomes had been implicated differentially expressed in patient with moyamoya disease (MMD), but the miRNAs expression in circulating leukocytes remains unclear. This study was investigated on the differential expression of miRNAs in peripheral leukocytes between MMD patients and healthy adults, and among patients with subtypes of MMD.

Materials and methods: A total of 30 patients with MMD and 10 healthy adults were enrolled in a stroke center from October 2017 to December 2018. The gene microarray was used to detect the differential expression profiles of miRNA in leukocytes between MMD patients and controls, and the differentially expressed miRNAs were verified by the method of real-time PCR. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were used to explore the key signaling pathways and possible pathogenesis of MMD.

Results: The microarray results showed 12 differentially expressed miRNAs in leukocytes of MMD patients compared with controls (fold change >2.0, p < 0.05 and FDR <0.05), of which 8 miRNAs were upregulated (miRNA-142-5p, miRNA-29b-3p, miRNA-424-5p, MiRNA-582-5p, miRNA-6807-5p, miRNA-142-3p, miRNA-340-5p, miRNA-4270), and 4 miRNAs were downregulated (miRNA-144-3p, miRNA-451a, miRNA-486-5p, miRNA-363-3p). The real-time PCR confirmed seven differentially expressed miRNAs (p < 0.05), of which 4 miRNAs (miRNA-29b-3p, miRNA-142-3p, miRNA-340-5p, miRNA-582-5p) were upregulated, and 3 miRNAs (miRNA-363-3p, miRNA-451a and miRNA-486-5p) were downregulated. Both GO and KEGG analysis suggested that the Wnt signaling pathway may be involved in the pathogenesis of MMD. In addition, miRNAs were also differentially expressed among patients with subtypes of MMD.

Conclusion: This study indicated that miRNAs are differentially expressed in peripheral leukocytes between MMD patients and healthy adults, and among patients with subtypes of MMD. The Wnt signaling pathway is probably involved in the pathogenesis of MMD.

Introduction

Moyamoya disease (MMD) is a relatively rare cerebrovascular disease associated with recurrent stroke, characterised by progressive stenosis or occlusion of the circle of Willis with growth of pathological collaterals at the base of the brain (; ). The prevalence of MMD have increased around the world, especially in Japan, Korea, and China (; ; ; ; ; ). Currently, the etiology and pathogenesis of MMD are still unclear, so there is a lack of reliable molecular biological markers and effective drugs (; ; ). The incidence of MMD has been demonstrated as racial propensity and familial clustering, and some genes have been indicated to be abnormally expressed, which suggests that the MMD may be closely associated with gene expressions (; ; ; ; ). MMD can be manifested as ischemic, hemorrhagic, or asymptomatic, and hemorrhagic stroke is one of the main factors leading to acute death and severe disability in patients with MMD (; ). However, the mechanisms and predictors of hemorrhagic stroke associated with MMD are still limited to the vascular or hemodynamic characteristics, while the underlying pathophysiological mechanisms and biological processes are still unclear (; ; ; ).

MicroRNA (miRNA) is a kind of endogenous small RNA with a length of about 18–24 nucleotides, which plays an important regulatory role in the normal metabolism of cells and the development of diseases by inhibiting the translation of messenger RNA (mRNA) into proteins or promoting the degradation of mRNA (; ). Although miRNA has been investigated for several years, its relationship with MMD remains to be elucidated. It has been indicated that miRNAs were abnormally expressed in peripheral blood and cerebrospinal fluid of patients with MMD (; ; ; ). However, previous studies had mostly focused on the differences of miRNAs in exosomes, which cannot fully reflect the actual expression of intracellular miRNAs. Also, studies have shown that some mononuclear cells in peripheral blood can differentiate into vascular endothelial progenitor cells and vascular smooth muscle progenitor cells, and the pathophysiological process of MMD is often accompanied by inflammatory changes involving a variety of cytokines, suggesting that MMD may be closely associated with peripheral leukocytes (; ).

In this study, we investigated the differential expression of miRNA in peripheral leukocytes between MMD patients and healthy adults, and among patients with subtypes (hemorrhagic, ischemic, and asymptomatic) of MMD, to explore the pathogenesis of MMD, and provide the basis for clinical diagnosis, outcome prediction and therapeutic strategy of MMD.

Materials and Methods

Study Population

The study was performed according to the guidelines from the Helsinki Declaration, and was approved by the Research Ethics Committee of the hospital. Written informed consent was obtained from all the subjects or their legally authorized representatives. A total of 30 patients, who were diagnosed with MMD in our stroke center from October 2017 to December 2018, were consecutively included in this study. The inclusion criteria included: 1) the age was 18–65 years; 2) diagnosed as MMD according to the 2012 Japanese Moyamoya disease diagnosis and treatment guidelines (); 3) the patient’s informed consent to the study protocol. The exclusion criteria included: 1) moyamoya syndrome indicated by clinical manifestations, high-resolution MRI, or laboratory examinations; 2) with history of stroke in the past 3 months; 3) patients who had received revascularization surgery; 4) patients with immunological diseases, tumors, or pregnancy. In addition, 10 healthy adults with matching gender and age were included as normal controls.

Sample Collection and RNA Extraction

Peripheral venous blood samples (about 6 ml) were taken from the antecubital vein into an anticoagulant drying tube under fasting conditions. The whole blood was centrifuged at 2,600 rpm for 10 min at room temperature within an hour. The white blood cells were separated from the middle layer, and the mixed red blood cells were removed with red blood cell lysis buffer. A kind of RNA stabilizer (RNAlater, Thermo Fisher Scientific, Baltics UAB, Vilnius, Lithuania) was added to the cryotube containing white blood cells and was stored at 4°C for 24 h, then stored in −80°C refrigerator for standby. All the samples (from MMD patients and controls) were collected and processed in the same way. Trizol Reagent (Invitrogen, Carlsbad, CA, United States) was used for leukocyte denaturation. The miRNeasy Mini Kit (Qiagen p/n 217004) was used for RNA collection and purification according to the manufacturer’s protocol. The quality evaluation of the total RNA was performed using the Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, United States).

MicroRNA Microarray

The Agilent Human (8*60K) V21.0 miRNA microarrays (Agilent, Santa Clara, CA, United States) were performed using a Gene Expression Hybridization Kit (Agilent’s miRNA Complete Labeling and Hyb Kit-p/n5190-0456) according to the manufacturer’s instructions. Slides were washed in staining dishes with a Gene Expression Wash Buffer Kit (Agilent, Santa Clara, CA, United States) and scanned by an Agilent Microarray Scanner with default settings according to the manufacturer’s instructions. The data was extracted using an Agilent Feature Extraction (AFE) software (version 10.7.1.1), and the raw data were normalized by Quantile algorithm using Gene Spring Software 12.6 (Agilent Technologies). Differentially expressed miRNAs were identified through the filtering of fold-change (>2), p value (<0.05) and FDR (<0.05) using R software (version 3.2.3) with the samr package. Also, the differentially expressed miRNAs should not have probes with flag A in at least one group.

Real-Time Polymerase Chain Reaction

Several differentially expressed miRNAs were quantified by using real-time polymerase chain reaction (RT-PCR) based on the Taqman method. A total of 100 ng of purified RNA was reversely transcribed to cDNA using the Taqman MicroRNA Reverse Transcription Kit (ABI, 4366597) and the stem-loop primers in TaqMan MicroRNA Assays (ABI, United States) according to the manufacturer’s instructions. The RT-PCR was performed using QuantStudio 5 Real-Time PCR System (ABI, United States) and QuantStudio™ Design & Analysis Software with the Taqman Universal PCR Master Mix (ABI, 4440049). The program was as follows: 95°C, 10 min; (95°C, 15 s; 60°C, 1 min) 40 cycles. Expression threshold for each miRNA detector was automatically determined. The homogenous miRNA-16 was used as an internal reference for the analysis. The relative expression level of each miRNA was calculated as fold change adopting the 2−ΔΔCt method. All RT-PCR analyses were performed in triplicate to test the reproducibility.

Receiver Operating Characteristic Curve Analysis

To assess the potential diagnostic value of identified miRNAs for MMD and subtypes of MMD, receiver operating characteristic (ROC) curve analysis of individual miRNA and combined model of multiple miRNAs was conducted.

Bioinformatics Analysis

In this study, we used TargetScan (http://www.targetscan.org/, TargetScanHuman 7.2) to identify the targets of differentially expressed miRNAs. We defined the predicted target genes by no less than 10 differentially expressed miRNAs as potential functional targets in MMD. To determine the biological relationship between the potential miRNA target genes, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were performed using DAVID Bioinformatics Resources (http://david.abcc.ncifcrf.gov/, DAVID Bioinformatics Resources 6.8). The significance threshold was set to 0.05 in our enrichment analysis.

Statistical Analysis

The statistical analysis was performed using SPSS (version 22.0, IBM-SPSS, Chicago, IL, United States) and R software (version 3.2.3). Genes with a two-sided p value of <0.05 and fold change >2.0 were regarded as statistically significant genes. ROC curve analysis of individual miRNA and combined model of multiple miRNAs was conducted to assess the potential diagnostic value of identified miRNAs for MMD.

Results

Demographic Information and Clinical Characteristics

The study recruited 30 patients with MMD (including 10 hemorrhagic MMD patients, 12 ischemic MMD patients, and 8 asymptomatic MMD patients), and 10 healthy adults with matched demographic characteristics. The detailed characteristics of the included subjects were shown in Table 1.

TABLE 1

IDSubtypesGenderAge (years)Suzuki stage
LeftRight
HM_01Hemorrhagic MMDFemale3244
HM_02Hemorrhagic MMDFemale4323
HM_03Hemorrhagic MMDFemale2433
HM_04Hemorrhagic MMDFemale4966
HM_05Hemorrhagic MMDMale4963
HM_06Hemorrhagic MMDMale4255
HM_07Hemorrhagic MMDFemale5033
HM_08Hemorrhagic MMDFemale3433
HM_09Hemorrhagic MMDMale5003
HM_10Hemorrhagic MMDFemale4244
IM_01Ischemic MMDMale3433
IM_02Ischemic MMDFemale3631
IM_03Ischemic MMDFemale4633
IM_04Ischemic MMDMale3533
IM_05Ischemic MMDFemale3043
IM_06Ischemic MMDfemale4633
IM_07Ischemic MMDMale4743
IM_08Ischemic MMDFemale2123
IM_09Ischemic MMDFemale2922
IM_10Ischemic MMDMale3954
IM_11Ischemic MMDMale3722
IM_12Ischemic MMDMale3922
AM_01Asymptomatic MMDMale4222
AM_02Asymptomatic MMDMale3453
AM_03Asymptomatic MMDFemale2933
AM_04Asymptomatic MMDFemale3833
AM_05Asymptomatic MMDFemale3444
AM_06Asymptomatic MMDFemale4822
AM_07Asymptomatic MMDMale5033
AM_08Asymptomatic MMDFemale3365
C_01ControlMale42--
C_02ControlMale35--
C_03ControlMale35--
C_04ControlMale34--
C_05ControlFemale34--
C_06ControlFemale32--
C_07ControlFemale33--
C_08ControlFemale35--
C_09ControlFemale31--
C_10ControlFemale33--

Demographic and clinical characteristics of the subjects.

HM, hemorrhagic moyamoya disease; IM, ischemic moyamoya disease; AM, asymptomatic moyamoya disease; C, controls.

MicroRNA Microarray

In analysis of 30 MMD patients compared with 10 controls, 8 miRNAs were found to be differentially expressed, with 5 upregulated (miRNA-142-5p, miRNA-29b-3p, miRNA-424-5p, miRNA-582-5p, miRNA-6807-5p) and 3 downregulated (miRNA-144-3p, miRNA-451a, miRNA-486-5p) (Table 2, Figure 1). In addition, according to the intra-group correlation analysis, there were significant differences in RNA expression in 6 samples (4 in the MMD group and 2 in the control group). After the elimination of these 6 samples, miRNA-142-3p, miRNA-340-5p, miRNA-363-3p, miRNA-4270 were also differentially expressed in MMD patients (miRNA-142-3p, miRNA-340-5p and miRNA-4270 upregulated, while miRNA-363-3p downregulated).

TABLE 2

Systematic namep valuesFDRFold changeRegulation
miRNA-142-5p0.0000200.0022102.020448Up
miRNA-29b-3p0.0000060.0013112.420699Up
miRNA-424-5p0.0000430.0032482.309035Up
miRNA-582-5p0.0012930.0346482.228220Up
miRNA-6807-5p0.0000100.0016282.227829Up
miRNA-142-3pa0.0000180.0012912.160976Up
miRNA-340-5pa0.0000250.0015492.047370Up
miRNA-4270a0.0004780.0112712.110506Up
miRNA-363-3pa0.0001100.0041340.492900Down
miRNA-144-3p0.0009580.0271250.489388Down
miRNA-451a0.0000680.0046950.403396Down
miRNA-486-5p0.0001650.0091260.181111Down

Differentially expressed miRNAs in microarray analysis between MMD patients and controls.

a

The four miRNAs were found differentially expressed in MMD patients after elimination of the 6 samples, in which miRNA expression was significantly different from others according to the intra-group correlation analysis.

FIGURE 1

In analysis of 10 hemorrhagic MMD patients compared with 10 controls, 17 miRNAs were found to be differentially expressed, with 13 upregulated (miRNA-101-3p, miRNA-142-3p, miRNA-142-5p, miRNA-148a-3p, miRNA-29b-3p, miRNA-340-5p, miRNA-374a-5p, miRNA-424-5p, miRNA-450a-5p, miRNA-542-3p, miRNA-582-5p, miRNA-590-5p, miRNA-6807-5p) and 4 downregulated (miRNA-144-3p, miRNA-144-5p, miRNA-451a, miRNA-486-5p). In analysis of 8 asymptomatic MMD patients compared with 10 controls, 7 miRNAs were found to be differentially expressed, with 3 upregulated (miRNA-142-3p, miRNA-29b-3p, miRNA-6515-5p) and 4 downregulated (miRNA-144-3p, miRNA-144-5p, miRNA-451a, miRNA-486-5p). In analysis of 12 ischemic MMD patients compared with 10 controls, 3 miRNAs were found to be differentially expressed (miRNA-29b-3p and miRNA-424-5p upregulated, while miRNA-486-5p downregulated) (Figure 2).

FIGURE 2

In addition, miRNA-486-5p was found to be differentially downregulated in hemorrhagic MMD patients, compared with ischemic MMD patients. Also, miRNA-6515-5p was found to be differentially upregulated in asymptomatic MMD patients, compared with symptomatic (hemorrhagic and ischemic) MMD patients.

Real-Time Polymerase Chain Reaction

The RT-PCR verification of 10 differentially expressed miRNAs was performed between 30 MMD patients and 10 controls, including miRNA-142-5p, miRNA-29b-3p, miRNA-424-5p, miRNA-582-5p, miRNA-144-3p, miRNA-451a, miRNA-486-5p, miRNA-142-3p, miRNA-340-5p, miRNA-363-3p. The other two miRNAs (miRNA-4270 and miRNA-6807-5p) were not validated because the potential functional targets had not been well established. The probe was synthesized by Thermo Fisher Scientific, and the sequences of the miRNAs in RT-PCR is shown in Table 3. All RT-PCR analyses were performed in triplicate, which demonstrated strong reproducibility, and all of the 10 genes were identified successfully. However, two normal control samples had higher overall CT values than other samples, but miRNA-16 (internal reference) was relatively lower, which led to higher overall expression and affected the overall analysis results, so the two samples were removed during the final analysis. The fold change in the microarray and RT-PCR for all the validated miRNAs was shown in Figure 3.

TABLE 3

miRNATarget sequence of the probe
miRNA-16UAG​CAG​CAC​GUA​AAU​AUU​GGC​G
miRNA-29b-3pUAG​CAC​CAU​UUG​AAA​UCA​GUG​UU
miRNA-142-3pUGU​AGU​GUU​UCC​UAC​UUU​AUG​GA
miRNA-142-5pCAU​AAA​GUA​GAA​AGC​ACU​ACU
miRNA-144-3pUAC​AGU​AUA​GAU​GAU​GUA​CU
miRNA-340-5pUUA​UAA​AGC​AAU​GAG​ACU​GAU​U
miRNA-363-3pAAU​UGC​ACG​GUA​UCC​AUC​UGU​A
miRNA-424-5pCAG​CAG​CAA​UUC​AUG​UUU​UGA​A
miRNA-451aAAA​CCG​UUA​CCA​UUA​CUG​AGU​U
miRNA-486-5pUCC​UGU​ACU​GAG​CUG​CCC​CGA​G
miRNA-582-5pUUA​CAG​UUG​UUC​AAC​CAG​UUA​CU

The sequences of the validated miRNAs in real-time PCR.

FIGURE 3

Compared with 10 healthy adults, 7 miRNAs had differential expression in 30 MMD patients (p < 0.05), of which miRNA-29b-3p, miRNA-142-3p, miRNA-340-5p, miRNA-582-5p were upregulated, while miRNA-363-3p, miRNA-451a, miRNA-486-5p were downregulated (Figure 4).

FIGURE 4

Compared with 10 healthy adults, 9 miRNAs had differential expression in 10 hemorrhagic MMD patients (p < 0.05), of which miRNA-29b-3p, miRNA-142-3p, miRNA-142-5p, miRNA-340-5p, miRNA-582-5p were upregulated, while miRNA-451a, miRNA-486-5p, miRNA-144-3p, miRNA-363-3p were downregulated. Compared with 10 healthy adults, 6 miRNAs had differential expression in 8 asymptomatic MMD patients (p < 0.05), of which miRNA-29b-3p, miRNA-142-3p, miRNA-340-5p, miRNA-582-5p were upregulated, while miRNA-363-3p, miRNA-451a were downregulated. Compared with 10 healthy adults, 4 miRNAs had differential expression in 12 ischemic MMD patients (p < 0.05), of which miRNA-29b-3p, miRNA-142-5p, miRNA-340-5p were upregulated, while miRNA-451a was downregulated.

In addition, the miRNA-142-3p was upregulated in hemorrhagic MMD patients and miRNA-340-5p was upregulated in hemorrhagic and the asymptomatic MMD patients compared with ischemic MMD patients (p < 0.05). The miRNA-486-5p was downregulated in hemorrhagic and asymptomatic MMD patients, but there was no statistically significant difference (p > 0.05). There was no significant differential expression between the hemorrhagic and the asymptomatic MMD patients.

Diagnostic Value of Identified MicroRNAs

To determine the diagnostic values of the significantly differentially expressed miRNAs (miRNA-29b-3p, miRNA-142-3p, miRNA-340-5p, miRNA-582-5p, miRNA-363-3p, miRNA-451a and miRNA-486-5p), ROC curve analysis was performed. The details of area under the curve (AUC) and 95% confidence intervals (95% CI), sensitivity, and specificity of individual miRNA and combined model for diagnosis of MMD can be found in Table 4 and Figure 5.

TABLE 4

miRNAsDisease entityAUC (95% CI)Sensitivity (%)Specificity (%)
miRNA-142-3pMMD0.818 (0.665–0.972)77.3100.0
miRNA-29b-3pMMD0.945 (0.873–1.000)77.3100.0
miRNA-340-5pMMD0.968 (0.910–1.000)95.590.0
miRNA-363-3pMMD0.764 (0.507–1.000)86.480.0
miRNA-451aMMD0.882 (0.760–1.000)72.7100.0
miRNA-486-5pMMD0.905 (0.799–1.000)77.3100.0
miRNA-582-5pMMD0.864 (0.728–0.999)81.890.0
Combined model of 7 miRNAsMMD1100100
Combined model of 7 miRNAsHemorrhagic MMD1100100
Combined model of 7 miRNAsAsymptomatic MMD1100100
Combined model of 7 miRNAsIschemic MMD1100100

The details of the AUC, sensitivity, and specificity of each identified miRNA and combined model for diagnosis of MMD.

AUC(95%CI), area under the curve (95% confidence interval).

FIGURE 5

Target Gene Prediction and Bioinformatics Analysis

We used the TargetScan target gene prediction program and miRTarBase program to predict the potential target genes of differentially expressed miRNAs by microarray and RT-PCR. According to the seven differentially expressed miRNAs in peripheral leukocytes of MMD patients, a total of 3,728 target genes were predicted based on the TargetScan program, and 63 target genes were predicted based on TargetScan program and miRTarBase program (Table 5). Of the 3,728 target genes, the biological processes and signal transduction pathways that may be involved in the enrichment analysis are carried out through GO and KEGG. The associated biological processes and signal pathways were screened (according to the criteria: p value ≤0.05, number of target genes ≥10) and sequenced according to the enrichment degree, and the top 30 biological processes and signal pathways are shown in Figure 6. Both the GO and KEGG analysis suggested that the Wnt signaling pathway may be involved in the pathogenesis of MMD (Figures 6, 7).

TABLE 5

miRNAsTarget genes
miRNA-142-3pBOD1, PROM1, PTPN23, HMGA1, ARNTL, LPP, EGR2, CCNT2, RAC1, ROCK2, TGFBR1, LRRC32, TAB2, HSPA1B, APC, HMGB1
miRNA-29b-3pPPP1R13B, AQP4, ITGA6, AKT2, LOX, HDAC4, DUSP2, DNMT3B, LRP6, SERPINH1, NKIRAS2, LAMC1, COL3A1, BACE1, COL4A1, ADAM12, PPIC, DNMT3A
miRNA-340-5pRHOA, IL4, CCNG2
miRNA-363-3pS1PR1, BCL2L11, FBXW7
miRNA-451aADAM10, OXTR, MIF, CPNE3, CAB39, IL6R, MAP3K1, RAB5A, RAB14, CDKN2D, TSC1
miRNA-486-5pSMAD2, FOXO1, PIK3R1, OLFM4, CDK4, DOCK3, PTEN, ARHGAP5, CIT
miRNA-582-5pCASP3, CREB1, RAB27A

Target genes of identified miRNAs based on TargetScan program and miRTarBase program.

FIGURE 6

FIGURE 7

Discussions

The results of this study suggested that miRNA-29b-3p, miRNA-142-3p, miRNA-340-5p, and miRNA-582-5p were significantly upregulated, while miRNA-363-3p, miRNA-451a, miRNA-486-5p were significantly downregulated in peripheral leukocytes of patients with MMD compared with healthy controls. In addition, miRNAs were also differentially expressed among patients with subtypes of MMD.

Previous studies had suggested that miRNA-106b, miRNA-130a, miRNA-126, miRNA-125-3p, miRNA-196a, let-7c were abnormally expressed in the peripheral serum of MMD patients (; ; ). Also, the miRNA-3679-5p, miRNA-6165, miRNA-6760-5p, and miRNA-574-5p had been indicated to be differentially expressed in the cerebrospinal fluid of MMD patients (). The miRNA-196a was suggested to regulate cell proliferation and apoptosis by regulating the expression of ANXA1 gene in endothelial cells and vascular smooth muscle cells, which may be related to the onset of MMD (). It was indicated that the serum let-7c in MMD patients was significantly increased, which can bind to the 3′non-coding region of RNF213 transcribed mRNA, thereby affecting the biological activity of RNF213(21). The RNF213 mutation is currently considered to be a susceptible gene for MMD in Asian populations, which may be involved in the pathogenesis of MMD (). In the present study, however, these miRNAs did not show abnormal expression in peripheral leukocytes of MMD patients compared with controls, which suggested that there may be significant differences in the expression of miRNA between peripheral leukocytes and serum of patients with MMD. The circulating serum miRNAs reflect only the level of extracellular miRNAs secreted by a variety of cells through exosomes, but cannot fully reflect the true expression of specific intracellular miRNAs, which impeded further investigations on the miRNA-related target genes and cell functions.

The pathological characteristics of MMD are mainly manifested as a gradual degeneration of the smooth muscle in the arterial media, and the abnormal proliferation of the smooth muscle in the arterial intima, leading to progressive stenosis and occlusion of the involved vessels, accompanied by the development of collaterals (; ; ). The present study indicated that several biological processes or signal transduction pathways may be involved in the pathophysiological process of MMD, including vascular smooth muscle cell proliferation and metastasis, apoptosis, neovascularization, immune and inflammatory reactions, and thrombosis.

Both the GO analysis and KEGG analysis suggested that the Wnt signaling pathway may be involved in the pathogenesis of MMD (Figures 6, 7). The Wnt signaling pathway is a complex regulatory network consisting of the secreted Wnt protein family, the transmembrane receptor Frizzled family, CK1, Deshevelled, GSK3, APC, Axin, β-catenin, and TCF/Lef family (). When the extracellular ligand binds to the cell surface receptor, the intracellular segment of the surface receptor was activated, which transmits extracellular signals into the cells. Previous studies have shown that the Wnt signaling pathway is activated in some cardiovascular and cerebrovascular diseases, and plays an important role in the proliferation, migration, apoptosis, and differentiation of vascular smooth muscle cells, which is one of the mechanisms of MMD (; ; ; ). Therefore, the role of the Wnt signaling pathway in MMD is worthy of further investigation.

In addition, miRNAs were also differentially expressed among patients with subtypes of MMD. Patients with MMD need to undergo a chronic conversion of cerebral hemodynamics and feeding arteries, from the original internal carotid artery system to the external carotid artery system or vertebrobasilar artery system. Depending on the degree or weight of vascular occlusion and collateral vasodilation, MMD can be manifested as ischemic, asymptomatic, or hemorrhagic (; ). In this study, the microarray results demonstrated an increasing number of differentially expressed miRNAs from ischemic MMD to asymptomatic MMD to hemorrhagic MMD compared with controls. Also, the real-time PCR results suggested that the miRNA-142-3p was upregulated in hemorrhagic MMD patients, and the miRNA-340-5p was upregulated in hemorrhagic and asymptomatic MMD patients compared with ischemic MMD patients. This indicated that miRNA may be involved in the process of arterial occlusion and collateral development.

In this study, we found that the miRNAs are abnormally expressed in peripheral leukocytes of patients with MMD. A couple of genes had been reported to be mutated in patients with MMD, especially RNF213, which has a high mutation rate in Asian MMD patients (; ; ; ; ). However, the causal relationship and associated mechanism between those mutated genes and MMD is still unclear. In our future studies, we will investigate whether the identified miRNAs are correlated with the mutated genes reported previously. After further verification of real-time PCR with larger sample size, and verification of target genes regulation and cell function regulation, these differentially expressed miRNAs are likely to become new molecular biological markers of MMD and associated stroke, and hold a promise to become new potential targets of therapeutic strategy.

Potential limitations of our studies should be mentioned. First, all of the patients in this study were enrolled from a single center, so potential selection bias may be inevitable. We, however, tried our best to reduce in-house selection bias by collecting patients consecutively. Second, the case-control nature of the study necessitates further prospective cohorts to confirm our conclusions. Third, the sample size of this study was relatively small, and further PCR verification with cross-validation is needed in a large population. Fourth, we used target gene prediction and functional enrichment analysis to indirectly obtain the biological processes that may be involved in the pathogenesis, which lacked verification of miRNA’s regulation of target genes and cell functions, and further vitro cytology experiments are necessitated.

Conclusion

This study indicated that miRNAs are differentially expressed in peripheral leukocytes between MMD patients and healthy adults, and among patients with subtypes of MMD. The Wnt signaling pathway is probably involved in the pathogenesis of MMD.

Statements

Data availability statement

The data presented in the study are deposited in the GEO repository, accession number GSE178501. The datasets can also be found in the article/Supplementary Material.

Ethics statement

The studies involving human participants were reviewed and approved by Institutional Review Board (IRB) of Beijing Tiantan Hospital. The patients/participants provided their written informed consent to participate in this study.

Author contributions

KK and YS contributed to the conception/design of the study, the acquisition, analysis and interpretation of data, drafting of the manuscript, final approval of the version to be published. QZ, JL, YJ, RJ, NL, JW, BY, JL, and XL were responsible for the acquisition of data, final approval of the version to be published. XZ and DZ were responsible for the conception/design of the study, revision of the work, critical revision of the manuscript, final approval of the version to be published, and agreement to be accountable for all aspects of the work.

Funding

This study was supported by the National Science and Technology Major Project (2017ZX09304018), Chinese Academy of Medical Sciences Innovation Fund for Medical Sciences (2019-I2M-5-029), Beijing Municipal Committee of Science and Technology (Z201100005620010), National Natural Science Foundation of China (81371293), and National Natural Science Foundation of China (81501001).

Acknowledgments

We thank all of the patients and health care providers who participated in the present study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.816919/full#supplementary-material

References

  • 1

    BangO. Y.FujimuraM.KimS.-K. (2016). The Pathophysiology of Moyamoya Disease: An Update. J. Stroke18 (1), 1220. 10.5853/jos.2015.01760

  • 2

    BaoX.-Y.WangQ.-N.ZhangY.ZhangQ.LiD.-S.YangW.-Z.et al (2019). Epidemiology of Moyamoya Disease in China: Single-Center, Population-Based Study. World Neurosurg.122, e917e923. 10.1016/j.wneu.2018.10.175

  • 3

    ChenP.-C.YangS.-H.ChienK.-L.TsaiI.-J.KuoM.-F. (2014). Epidemiology of Moyamoya Disease in Taiwan. Stroke45 (5), 12581263. 10.1161/strokeaha.113.004160

  • 4

    DaiD.LuQ.HuangQ.YangP.HongB.XuY.et al (2014). Serum miRNA Signature in Moyamoya Disease. PloS one9 (8), e102382. 10.1371/journal.pone.0102382

  • 5

    DolzS.GórrizD.TemblJ. I.SánchezD.ForteaG.ParkhutikV.et al (2017). Circulating MicroRNAs as Novel Biomarkers of Stenosis Progression in Asymptomatic Carotid Stenosis. Stroke48 (1), 1016. 10.1161/strokeaha.116.013650

  • 6

    DuanL.BaoX.-Y.YangW.-Z.ShiW.-C.LiD.-S.ZhangZ.-S.et al (2012). Moyamoya Disease in China. Stroke43 (1), 5660. 10.1161/strokeaha.111.621300

  • 7

    DuanL.WeiL.TianY.ZhangZ.HuP.WeiQ.et al (2018). Novel Susceptibility Loci for Moyamoya Disease Revealed by a Genome-wide Association Study. Stroke49 (1), 1118. 10.1161/strokeaha.117.017430

  • 8

    FoulquierS.DaskalopoulosE. P.LluriG.HermansK. C. M.DebA.BlankesteijnW. M. (2018). WNT Signaling in Cardiac and Vascular Disease. Pharmacol. Rev.70 (1), 68141. 10.1124/pr.117.013896

  • 9

    JungK.-H.ChuK.LeeS.-T.ParkH.-K.KimD.-H.KimJ.-H.et al (2008). Circulating Endothelial Progenitor Cells as a Pathogenetic Marker of Moyamoya Disease. J. Cereb. Blood Flow Metab.28 (11), 17951803. 10.1038/jcbfm.2008.67

  • 10

    KangH.-S.MoonY.-J.KimY.-Y.ParkW.-Y.ParkA. K.WangK.-C.et al (2014). Smooth-muscle Progenitor Cells Isolated from Patients with Moyamoya Disease: Novel Experimental Cell Model. Jns120 (2), 415425. 10.3171/2013.9.jns131000

  • 11

    KangS.LiuX.ZhangD.WangR.ZhangY.ZhangQ.et al (2019). Natural Course of Moyamoya Disease in Patients with Prior Hemorrhagic Stroke. Stroke50 (5), 10601066. 10.1161/strokeaha.118.022771

  • 12

    KimJ. S. (2016). Moyamoya Disease: Epidemiology, Clinical Features, and Diagnosis. J. Stroke18 (1), 211. 10.5853/jos.2015.01627

  • 13

    KimK. M.KimJ. E.ChoW.-S.KangH.-S.SonY.-J.HanM. H.et al (2017). Natural History and Risk Factor of Recurrent Hemorrhage in Hemorrhagic Adult Moyamoya Disease. Neurosurgery81 (2), 289296. 10.1093/neuros/nyw179

  • 14

    KimT.LeeH.BangJ. S.KwonO.-K.HwangG.OhC. W. (2015). Epidemiology of Moyamoya Disease in Korea: Based on National Health Insurance Service Data. J. Korean Neurosurg. Soc.57 (6), 390395. 10.3340/jkns.2015.57.6.390

  • 15

    KuriyamaS.KusakaY.FujimuraM.WakaiK.TamakoshiA.HashimotoS.et al (2008). Prevalence and Clinicoepidemiological Features of Moyamoya Disease in Japan. Stroke39 (1), 4247. 10.1161/strokeaha.107.490714

  • 16

    KurodaS.HoukinK. (2008). Moyamoya Disease: Current Concepts and Future Perspectives. Lancet Neurol.7 (11), 10561066. 10.1016/s1474-4422(08)70240-0

  • 17

    KurodaS.KashiwazakiD.IshikawaT.NakayamaN.HoukinK. (2013). Incidence, Locations, and Longitudinal Course of Silent Microbleeds in Moyamoya Disease. Stroke44 (2), 516518. 10.1161/strokeaha.112.678805

  • 18

    LiuP.HanC.LiD.-S.LvX.-L.LiY.-X.DuanL. (2016). Hemorrhagic Moyamoya Disease in Children. Stroke47 (1), 240243. 10.1161/strokeaha.115.010512

  • 19

    MahajanS. G.FenderA. C.Meyer-KirchrathJ.KurtM.BarthM.SagbanT. A.et al (2012). A Novel Function of FoxO Transcription Factors in Thrombin-Stimulated Vascular Smooth Muscle Cell Proliferation. Thromb. Haemost.108 (1), 148158. 10.1160/TH11-11-0756

  • 20

    MenetR.LecordierS.ElAliA. (2020). Wnt Pathway: An Emerging Player in Vascular and Traumatic Mediated Brain Injuries. Front. Physiol.11, 565667. 10.3389/fphys.2020.565667

  • 21

    MoriM. A.LudwigR. G.Garcia-MartinR.BrandãoB. B.KahnC. R. (2019). Extracellular miRNAs: From Biomarkers to Mediators of Physiology and Disease. Cel Metab.30 (4), 656673. 10.1016/j.cmet.2019.07.011

  • 22

    MoriokaM.HamadaJ.-I.KawanoT.TodakaT.YanoS.KaiY.et al (2003). Angiographic Dilatation and branch Extension of the Anterior Choroidal and Posterior Communicating Arteries Are Predictors of Hemorrhage in Adult Moyamoya Patients. Stroke34 (1), 9095. 10.1161/01.str.0000047120.67507.0d

  • 23

    ParkY. S.JeonY. J.LeeB. E.KimT. G.ChoiJ.-U.KimD.-S.et al (2012). Association of the miR-146aC>G, miR-196a2C>T, and miR-499A>G Polymorphisms with Moyamoya Disease in the Korean Population. Neurosci. Lett.521 (1), 7175. 10.1016/j.neulet.2012.05.062

  • 24

    RcotpatosootcoWillis. (2012). Guidelines for Diagnosis and Treatment of Moyamoya Disease (Spontaneous Occlusion of the circle of Willis). Neurol. Med. Chir (Tokyo)52 (5), 245266. 10.2176/nmc.52.245

  • 25

    ScottR. M.SmithE. R. (2009). Moyamoya Disease and Moyamoya Syndrome. N. Engl. J. Med.360 (12), 12261237. 10.1056/nejmra0804622

  • 26

    ShaoY.ChenJ.FreemanW.DongL.-J.ZhangZ.-H.XuM.et al (2019). Canonical Wnt Signaling Promotes Neovascularization through Determination of Endothelial Progenitor Cell Fate via Metabolic Profile Regulation. Stem Cells37 (10), 13311343. 10.1002/stem.3049

  • 27

    TakahashiJ. C.FunakiT.HoukinK.InoueT.OgasawaraK.NakagawaraJ.et al (2016). Significance of the Hemorrhagic Site for Recurrent Bleeding. Stroke47 (1), 3743. 10.1161/strokeaha.115.010819

  • 28

    TsaousiA.WilliamsH.LyonC. A.TaylorV.SwainA.JohnsonJ. L.et al (2011). Wnt4/β-Catenin Signaling Induces VSMC Proliferation and Is Associated with Intimal Thickening. Circ. Res.108 (4), 427436. 10.1161/circresaha.110.233999

  • 29

    WangG.WenY.FaletiO. D.ZhaoQ.LiuJ.ZhangG.et al (2020). A Panel of Exosome-Derived miRNAs of Cerebrospinal Fluid for the Diagnosis of Moyamoya Disease. Front. Neurosci.14, 548278. 10.3389/fnins.2020.548278

  • 30

    WangX.WangY.NieF.LiQ.ZhangK.LiuM.et al (2020). Association of Genetic Variants with Moyamoya Disease in 13 000 Individuals. Stroke51 (6), 16471655. 10.1161/strokeaha.120.029527

  • 31

    WangY.YangL.WangX.ZengF.ZhangK.ZhangQ.et al (2021). Meta‐analysis of Genotype and Phenotype Studies to Confirm the Predictive Role of the RNF213 p.R4810K Variant for Moyamoya Disease. Eur. J. Neurol.28 (3), 823836. 10.1111/ene.14635

  • 32

    ZhangQ.LiuY.ZhangD.WangR.ZhangY.WangS.et al (2016). RNF213 as the Major Susceptibility Gene for Chinese Patients with Moyamoya Disease and its Clinical Relevance. J. Neurosurg., 126, 11061113. 10.3171/2016.2.JNS152173

  • 33

    ZhaoS.GongZ.ZhangJ.XuX.LiuP.GuanW.et al (2015). Elevated Serum MicroRNA Let-7c in Moyamoya Disease. J. Stroke Cerebrovasc. Dis.24 (8), 17091714. 10.1016/j.jstrokecerebrovasdis.2015.01.041

Summary

Keywords

moyamoya disease, stroke, non-coding RNA, microRNA, wnt signaling pathway

Citation

Kang K, Shen Y, Zhang Q, Lu J, Ju Y, Ji R, Li N, Wu J, Yang B, Lin J, Liang X, Zhang D and Zhao X (2022) MicroRNA Expression in Circulating Leukocytes and Bioinformatic Analysis of Patients With Moyamoya Disease. Front. Genet. 13:816919. doi: 10.3389/fgene.2022.816919

Received

21 December 2021

Accepted

14 April 2022

Published

20 May 2022

Volume

13 - 2022

Edited by

Fan Jin, Zhejiang University, China

Reviewed by

William K.K. Wu, Chinese University of Hong Kong, China

Yue-qiu Tan, Central South University, China

Updates

Copyright

*Correspondence: Xingquan Zhao, ; Dong Zhang,

† These authors have contributed equally to this work

This article was submitted to Genetics of Common and Rare Diseases, a section of the journal Frontiers in Genetics

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics