Skip to main content

Single-cell sequencing analysis of peripheral blood in patients with moyamoya disease

Abstract

Background

At present, the etiology of moyamoya disease is not clear, and it is necessary to explore the mechanism of its occurrence and development. Although some bulk sequencing data have previously revealed transcriptomic changes in Moyamoya disease, single-cell sequencing data has been lacking.

Methods

Two DSA(Digital Subtraction Angiography)-diagnosed patients with moyamoya disease were recruited between January 2021 and December 2021. Their peripheral blood samples were single-cell sequenced. CellRanger(10 x Genomics, version 3.0.1) was used to process the raw data, demultiplex cellular barcodes, map reads to the transcriptome, and dowm-sample reads(as required to generate normalized aggregate data across samples). There were 4 normal control samples, including two normal samples GSM5160432 and GSM5160434 of GSE168732, and two normal samples of GSE155698, namely GSM4710726 and GSM4710727. Weighted co-expression network analysis was used to explore the gene sets associated with moyamoya disease. GO analysis and KEGG analysis were used to explore gene enrichment pathways. Pseudo-time series analysis and cell interaction analysis were used to explore cell differentiation and cell interaction.

Results

For the first time, we present a peripheral blood single cell sequencing landscape of Moyamoya disease, revealing cellular heterogeneity and gene expression heterogeneity. In addition, by combining with WGCNA analysis in public database and taking intersection, the key genes in moyamoya disease were obtained. namely PTP4A1, SPINT2, CSTB, PLA2G16, GPX1, HN1, LGALS3BP, IFI6, NDRG1, GOLGA2, LGALS3. Moreover, pseudo-time series analysis and cell interaction analysis revealed the differentiation of immune cells and the relationship between immune cells in Moyamoya disease.

Conclusions

Our study can provide information for the diagnosis and treatment of moyamoya disease.

Introduction

Moyamoya disease is a cerebrovascular disease characterized by chronic progressive stenosis or occlusion at the end of bilateral internal carotid artery and the beginning of the anterior cerebral artery and middle cerebral artery, with secondary abnormal vascular network formation at the skull base [1]. Moyamoya disease was first proposed and named by Japanese scholars in 1957 [2]. “Moyamoya,“ which means “smoke” in Japanese, graphically describes the patient’s network of compensatory abnormal blood vessels swirling like smoke [3, 4]. In recent years, both the incidence and prevalence of moyamoya disease have gradually increased, which may depend on advances in diagnostic techniques [5, 6]. Globally, people in Asia have a higher incidence of moyamoya disease, especially in Japan and South Korea [7]. Notably, two peaks were found in the age of onset of Moyamoya disease: 5–10 years and 25–49 years [8,9,10]. The main clinical manifestation of Moyamoya disease is ischemic or hemorrhagic stroke, which has a high mortality rate [11]. However, the pathogenesis of moyamoya disease is still unknown. Previous genomics studies have identified the RNF213 gene as a susceptibility gene in moyamoya disease [12, 13]. However, cellular heterogeneity, transcriptome changes, signaling pathway changes, and intercellular communication in Moyamoya disease remain unclear. In-depth exploration is urgently needed to provide reference for the pathogenesis and treatment of moyamoya disease.

At present, single cell sequencing technology is a revolution in life science research [14]. Although the rapid development of next generation sequencing (NGS) technology and the third generation sequencing (TGS) technology has caused a huge change in the field of biological research, but the sequencing results are a representation of the cells as a whole [15, 16]. However, due to cell heterogeneity, the genetic information of cells with the same phenotype may be significantly different, and much information with low abundance will be lost in the overall characterization [17]. In order to make up for the limitations of traditional high-throughput sequencing, single-cell sequencing technology emerged. It performs high-throughput sequencing of the genome, transcriptome, and epigenome at the single cell level. It can reveal the gene structure and gene expression state of a single cell and reflect the heterogeneity between cells. It plays an important role in the fields of tumor, developmental biology, microbiology, neuroscience and so on, and is becoming the focus of life science research [18,19,20]. However, there is still a blank in the field of sequencing peripheral blood single cells in human patients with Moyamoya disease, which makes it difficult for us to deeply understand the transcriptome expression of various cells and the immune function in the microenvironment of moyamoya disease.

To explore the pathogenesis of moyamoya disease, single cell sequencing was performed in patients with moyamoya disease and normal controls. We performed weighted co-expression network analysis, enrichment analysis, protein-protein interaction analysis, pseudo-time series analysis, and intercellular communication analysis, in order to provide a reference for diagnosis and treatment of moyamoya disease and identification of markers.

Methods

Single cell data acquisition and processing

Two DSA(Digital Subtraction Angiography)-diagnosed patients with moyamoya disease were recruited between January 2021 and December 2021. They signed informed consent and peripheral blood samples were used for single-cell sequencing. The Ethics Committee of the First Affiliated Hospital of Nanjing Medical University approved the study (No. 2021-SR-165). CellRanger(10× Genomics, version 3.0.1) was used to process the raw data, demultiplex cellular barcodes, map reads to the transcriptome, and dowm-sample reads(as required to generate normalized aggregate data across samples). There were 4 normal control samples, including two normal samples GSM5160432 and GSM5160434 of GSE168732, and two normal samples of GSE155698, namely GSM4710726 and GSM4710727.

Quality control of single cell sequencing data

The “Seurat” package was used to process and analyze single cell sequencing data [21]. The criteria for inclusion of genes for subsequent analysis are as follows: (1) Genes expressed in at least 3 cells. (2) Genes whose expression level is above 200. The inclusion criteria were as follows: (1) The gene expression level was greater than 200 but less than 4000. (2) The percentage of mitochondrial gene is less than 25. (3) The total gene expression is less than 15,000.

Processing of data

First, the “LogNormalize” method is used to standardize data. Subsequently, the number of highly variable genes was set to 3000, and the method was “vst” to search for highly variable genes. The SCTransform function of the “Seurat” package was used to integrate the data and eliminate interference from the cell cycle and mitochondria. Set dims to 20 and use the tSNE method to reduce the dimensions. The “KNN” method was used to cluster cells, and K.Paam was set as 20, dims = 20, resolution = 0.8, random.seed = 2021. The “SingleR” package was used to annotate the cells.

Weighted gene co-expression Network Analysis (WGCNA)

The “sample” function was used to randomly select 2000 cells, including both Moyamoya disease and normal samples, for subsequent analysis. In this study, genes with variance fluctuation in the top 50% were selected for subsequent analysis. If the soft threshold ranges from 1 to 10, the step is (1) From 12 to 20, the step size is (2) The “WGNCA” package pickSoftThresgold function is used to explore the soft threshold changes. Set deepSplit = 2, the minimum number of module genes was 100, and the genes were clustered into different modules. Finally, each module was associated with the phenotype to find the module gene most associated with moyamoya disease.

Gene Ontology(GO) enrichment analysis

Gene Ontology(GO) is a database established by the Gene Ontology Consortium, which aims to establish a database applicable to various species to define and describe the functions of genes and proteins [22]. GO enrichment analysis is divided into three categories, namely Biological Process (BP), Cellular Component (CC) and Molecular Function (MF). We used GO enrichment analysis to explore the enrichment function of gene sets.

Kyoto Encyclopedia of genes and genomes (KEGG) enrichment analysis

KEGG(Kyoto Encyclopedia of Genes and Genomes) database is a database that systematically analyzes gene functions, genomic information and functional information, including metabolic pathways database, hierarchical classification database, gene database, genome database, etc [23]. By comparing the studied genes with the pathway gene set in the database, the enriched pathways were obtained.

Pseudo-time series analysis and cell interaction analysis

The “monocle2” package was used for cell differentiation trajectory analysis. Genes with average expression greater than 0.1 were selected for subsequent analysis. Set max_components = 2 and method = DDRTree to reduce the sample dimensions. Cell differentiation was shown by tSNE diagram. The “CellChat” packet was used for cell interaction analysis. The database “Secreted Signaling” was selected and the minimum number of cells was set to 10.

Statistical analysis

The method used to search for genes associated with moyamoya disease is WGCNA. The Chi-squared test and Fisher’s exact test were performed to compare the clinical variables. The “Seurat” package was used to process and analyze single cell sequencing data. The ClusterProfile package was used for enrichment analysis. The “monocle2” package was used for pseudo temporal analysis. The software used was R software, version 4.0.5. Without specific description, p < 0.05 was considered statistically significant.

Results

Figure 1 showed our workflow.

Fig. 1
figure 1

The flow chart

Single cell sequencing data analysis

In this study, two samples of Moyamoya disease and four normal samples were analyzed first. Through quality control, a total of 15,459 cells of Moyamoya disease and 20,079 normal control cells were obtained, as shown in Fig. 2A. Red represents Moyamoya disease samples and yellow represents normal samples. As shown in Fig. 2B, no significant batch effect was observed between the two samples of Moyamoya disease. Similarly, no significant batch effect was observed among the four samples of the normal control, while there was heterogeneity between the moyamoya disease and normal samples. Through dimensionality reduction and cluster analysis, all cells were coclustered into 25 clusters (Fig. 2C). A total of 7 cell types have been annotated, as shown in Fig. 2D: NK cells, T cells, B cells, HSC-G-CSF cells, Monocyte cells, Macrophage cells, Epithelial cells. In order to further explore the pathway activated in moyamoya disease, we first obtained significantly up-regulated genes in moyamoya disease by differential expressed gene analysis. The first five up-regulated genes were RPL11,PARK7, ENO1, RPL22 and DDOST. The five genes that were significantly down-regulated were EFHD2, SPEN, AL02AA55.5, PRDM2 and MICOS10 (Fig. 3A). In Fig. 3B-C, GO and KEGG enrichment analysis showed that functional pathways associated with moyamoya disease include RNA splicing, ribonucleoprotein complex biogenesis, and pathways of neurodegeneration multiple diseases, amyotrophic lateral sclerosis, etc. Then, we screened all moyamoya cells and explored the functional pathways. As shown in Fig. 3D, we found that most of the pathways were highly activated in Monocyte cells, Macrophage cells and Epithelial cells. However, the activation degree of NK cells, T cells and B cells was lower.

Fig. 2
figure 2

Single cell sequencing data analysis. (A) Cells obtained from Moyamoya disease (MMD) and normal controls (N). (B) Batch effect testing. (C) All cells were clustered into 25 clusters. (D) Cell annotation

Fig. 3
figure 3

Differentially expressed genes and their enrichment analysis. (A) Analysis of differentially expressed genes between moyamoya disease and normal controls. (B) GO enrichment analysis of differentially expressed genes. (C) KEGG enrichment analysis of differentially expressed genes. (D) Pathway enrichment of different immune cells

Weighted co-expression Network Analysis (WGCNA), functional enrichment and search for key genes

Then we used the “sample” function to randomly select 2000 samples of Moyamoya disease cells and performed WGCNA analysis. As shown in Fig. 4A, the optimal soft threshold was 2, and the data conform to the power law distribution. As shown in Fig. 4B-C, genes were aggregated into 8 non-gray modules, among which the red module was most closely related to the occurrence of moyamoya disease. As shown in Fig. 4D, module membership of genes in the red module was positively correlated with gene significance for body weight (cor = 0.51, p < 0.001). We then investigated the function of these red module genes, as shown in Fig. 5A-B. These genes are associated with gland development, epithelial cell development, cell-cell junction organization, focal adhesion, and ECM-receptor interaction. Then, we intersected the differential genes obtained from single-cell analysis with the red module genes, and a total of 11 genes were obtained. They were respectively PTP4A1, SPINT2, CSTB, PLA2G16, GPX1, HN1, LGALS3BP, IFI6, NDRG1, GOLGA2, LGALS3(Fig. 5C).

Fig. 4
figure 4

Weighted Coexpression network analysis (WGCNA). (A) Determination of optimal soft thresholds. (B-C) Cluster integration of all genes. The red module is considered the critical module. (D) Correlation of red modules

Fig. 5
figure 5

Enrichment analysis of red module genes obtained from WGCNA. (A) GO enrichment analysis. (B) KEGG enrichment analysis. (C) The intersection of red module genes of WGCNA and differentially expressed genes obtained by single-cell sequencing analysis resulted in 11 genes

Pseudo-time series analysis

Subsequently, immune cells were extracted from samples of Moyamoya disease for pseudo-time series analysis. As shown in Fig. 6A, darker blue means that cells differentiate earlier, from dark blue to light blue. These cells had a total of seven differentiation states (Fig. 6B). NK cells differentiated early, corresponding to differentiation state 1, and corresponding to two main clusters (Fig. 6C and D). Moreover, we found that the differentiation status of the two samples of Moyamoya disease was consistent (Fig. 6E). As shown in Fig. 6F, we found that the expression of 8 of the above 11 important genes was changed during the differentiation of immune cells. Among them, the expression of genes CSTB, LGALS3 and IFI6 increased first and then decreased, while the expression of GOLGA2 and NDRG1 decreased first and then increased, while the expression of SPINT2, LGALS3BP and PTP4A1 increased gradually.

Fig. 6
figure 6

Pseudo-time series analysis. (A) Darker blue indicates early cell differentiation, from dark blue to light blue. (B) The cells have a total of seven differentiated states. (C, D) NK cells differentiated early, corresponding to state 1, corresponding to two major clusters (E) The differentiation states of the two Moyamoya disease samples were consistent. (F) The expression of 8 of the 11 key genes was changed. Among them, the expressions of CSTB, LGALS3 and IFI6 genes were firstly increased and then decreased, the expressions of GOLGA2 and NDRG1 were firstly decreased and then increased, and the expressions of SPINT2, LGALS3BP and PTP4A1 were gradually increased

Cell interaction analysis

Cell interaction analysis was performed to further explore the connections between various cells in Moyamoya disease. Figure 7 A showed the correlation between the five cell types and to what extent. The thicker the line, the closer the connection. As shown in Fig. 7B, B cells were mainly associated with T cells, NK cells and Monocytes, of which Monocyte was the most closely associated. As shown in Fig. 7C, the BTLA-TNFRSF14 pathway was involved in the connection between B cells and other cells, suggesting that this pathway may play an important role. As shown in Fig. 7D, in the BTLA pathway, B cells mainly send and affect signal transmission, while Monocyte cells mainly receive signals. As shown in Fig. 7E, according to the degree of cell communication, cells were divided into four modes, distributed as parten1-4. Among them, B cells were mainly in parten3, and the main gene expressed in parten3 is BTLA. As shown in Fig. 7F, Sankey diagram similarly showed that B cells correspond to parten3, and BTLA was mainly expressed in parten3.

Fig. 7
figure 7

Cell interaction analysis. (A) The correlation between the five cell types. (B) B cells are mainly associated with T cells, NK cells and Monocyte cells, of which Monocyte is the most closely associated. (C) The BTLA-TNFRSF14 pathway is involved in the connection between B cells and other cells, suggesting that this pathway may play an important role. (D) In the BTLA pathway, B cells mainly send and affect signal transmission, while Monocyte cells mainly receive signals. (E) Cells are divided into parten1-4. Among them, B cells are mainly in parten3, and the main gene expressed in parten3 is BTLA. (F) Sankey diagram similarly shows that B cells correspond to parten3, and BTLA is mainly expressed in parten3

Discussion

Stroke is one of the main causes of disability or death in middle-aged and elderly people, with high morbidity, mortality and disability rates [24]. Although most strokes are caused by high blood pressure, they can also be caused by an easily overlooked disease: Moyamoya disease [25]. Stroke caused by Moyamoya disease is associated with progressive narrowing of the intracranial internal carotid artery and its proximal branch [26]. Although vascular reconstructive surgery is capable of improving the progression of Moyamoya disease, identification of the disease at an early stage is critical in order to achieve the best results in patients [27]. However, the pathogenesis of moyamoya disease is still unclear at present, and multi-omics studies such as genome and transcriptome are urgently needed to reveal the potential mechanism of moyamoya disease formation, so as to provide reference data for the early diagnosis and treatment of Moyamoya disease.

For the first time, we explored the single-cell sequencing map of moyamoya disease peripheral blood using a variety of algorithms, including dimension reduction clustering, enrichment analysis, differential expression gene analysis, weighted coexpression network analysis, and intercellular communication analysis. Our results for the first time identified cellular heterogeneity in Moyamoya disease peripheral blood and identified potential key genes PTP4A1, SPINT2, CSTB, PLA2G16, GPX1, HN1, LGALS3BP, IFI6, NDRG1, GOLGA2, LGALS3. Moreover, the interaction between cells was clearly demonstrated, which provides a reference for us to understand the pathogenesis and the crosstalk between cells and pathways. This not only provided reference for us to understand the pathogenesis of moyamoya disease, but also provided us with a target pathway.

This study was the first single-cell sequencing analysis of Moyamoya disease. Compared with previous bulk RNA sequencing alone, our study was the first to cluster Moyamoya disease samples into different cell types through single-cell sequencing analysis, greatly reducing the heterogeneity of sequencing data, and studying the interaction mode between different cells. The BTLA-TNFRSF14 pathway identified by us may be a potential target for the treatment of Moyamoya disease. However, there are limitations to our study. First of all, due to the rarity of Moyamoya disease, we can obtain a small number of samples. Second, our normal controls came from a public database. We will expand the sample size and conduct further in-depth analysis in the future.

Conclusion

This study was the first single-cell sequencing analysis of Moyamoya disease. We explored the single-cell sequencing map of moyamoya disease peripheral blood using a variety of algorithms, including dimension reduction clustering, enrichment analysis, differential expression gene analysis, weighted coexpression network analysis, and intercellular communication analysis. Our results for the first time identified cellular heterogeneity in Moyamoya disease peripheral blood and identified potential key genes.

Availability of data and materials

All the data can be found in the article. The normal sample was downloaded from GEO, accession number GSE168732(https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE168732).

References

  1. Scott RM, Smith ER. Moyamoya disease and moyamoya syndrome. N Engl J Med. 2009 Mar 19;360(12):1226-37. doi: https://doi.org/10.1056/NEJMra0804622. PMID: 19297575.

  2. Shang S, Zhou D, Ya J, Li S, Yang Q, Ding Y, Ji X, Meng R. Progress in moyamoya disease. Neurosurg Rev. 2020 Apr;43(2):371–82. https://doi.org/10.1007/s10143-018-0994-5. Epub 2018 Jun 18. PMID: 29911252.

  3. Ihara M, Yamamoto Y, Hattori Y, Liu W, Kobayashi H, Ishiyama H, Yoshimoto T, Miyawaki S, Clausen T, Bang OY, Steinberg GK, Tournier-Lasserve E, Koizumi A. Moyamoya disease: diagnosis and interventions. Lancet Neurol. 2022 Aug;21(8):747–58. https://doi.org/10.1016/S1474-4422(22)00165-X. Epub 2022 May 20. PMID: 35605621.

  4. Hishikawa T, Sugiu K, Date I. Moyamoya Disease: A Review of Clinical Research. Acta Med Okayama. 2016 Aug;70(4):229 – 36. doi: https://doi.org/10.18926/AMO/54497. PMID: 27549666.

  5. Mertens R, Graupera M, Gerhardt H, Bersano A, Tournier-Lasserve E, Mensah MA, Mundlos S, Vajkoczy P. The genetic basis of Moyamoya Disease. Transl Stroke Res 2022 Feb;13(1):25–45. doi: https://doi.org/10.1007/s12975-021-00940-2. Epub 2021 Sep 16. PMID: 34529262; PMCID: PMC8766392.

  6. Asselman C, Hemelsoet D, Eggermont D, Dermaut B, Impens F. Moyamoya disease emerging as an immune-related angiopathy. Trends Mol Med. 2022 Nov;28(11):939–50. https://doi.org/10.1016/j.molmed.2022.08.009. Epub 2022 Sep 14. PMID: 36115805.

  7. Velo M, Grasso G, Fujimura M, Torregrossa F, Longo M, Granata F, Pitrone A, Vinci SL, Ferra첫 L, La Spina P. Moyamoya Vasculopathy: Cause, Clinical Manifestations, Neuroradiologic Features, and Surgical Management. World Neurosurg. 2022 Mar;159:409–425. doi: 10.1016/j.wneu.2021.11.026. PMID: 35255640.

  8. Zhang X, Xiao W, Zhang Q, Xia D, Gao P, Su J, Yang H, Gao X, Ni W, Lei Y, Gu Y. Progression in Moyamoya Disease: clinical features, neuroimaging evaluation, and treatment. Curr Neuropharmacol. 2022;20(2):292–308. https://doi.org/10.2174/1570159X19666210716114016. PMID: 34279201; PMCID: PMC9413783.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Pandey P, Steinberg GK. Neurosurgical advances in the treatment of moyamoya disease. Stroke. 2011 Nov;42(11):3304-10. doi: https://doi.org/10.1161/STROKEAHA.110.598565. Epub 2011 Oct 6. PMID: 21980214.

  10. Wan M, Duan L. Recent progress in hemorrhagic moyamoya disease. Br J Neurosurg. 2015 Apr;29(2):189 – 91. doi: 10.3109/02688697.2014.976177. Epub 2014 Nov 3. PMID: 25365664.

  11. Liu E, Zhao H, Liu C, Tan X, Luo C, Yang S. Research progress of moyamoya disease combined with renovascular hypertension. Front Surg. 2022 Aug;26:9:969090. https://doi.org/10.3389/fsurg.2022.969090. PMID: 36090342; PMCID: PMC9458923.

  12. Bang OY, Chung JW, Kim DH, Won HH, Yeon JY, Ki CS, Shin HJ, Kim JS, Hong SC, Kim DK, Koizumi A. Moyamoya Disease and Spectrums of RNF213 Vasculopathy. Transl Stroke Res 2020 Aug;11(4):580–9. doi: https://doi.org/10.1007/s12975-019-00743-6. Epub 2019 Oct 24. PMID: 31650369.

  13. Mineharu Y, Miyamoto S. RNF213 and GUCY1A3 in Moyamoya Disease: key regulators of metabolism, inflammation, and Vascular Stability. Front Neurol. 2021 Jul;26:12:687088. https://doi.org/10.3389/fneur.2021.687088. PMID: 34381413; PMCID: PMC8350054.

  14. Ziegenhain C, Vieth B, Parekh S, Reinius B, Guillaumet-Adkins A, Smets M, Leonhardt H, Heyn H, Hellmann I, Enard W. Comparative Analysis of Single-Cell RNA Sequencing Methods. Mol Cell. 2017 Feb 16;65(4):631–643.e4. doi: https://doi.org/10.1016/j.molcel.2017.01.023. PMID: 28212749.

  15. Gao S. Data Analysis in Single-Cell Transcriptome Sequencing. Methods Mol Biol. 2018;1754:311–326. doi: https://doi.org/10.1007/978-1-4939-7717-8_18. PMID: 29536451.

  16. van Dijk EL, Jaszczyszyn Y, Naquin D, Thermes C. The third revolution in sequencing technology. Trends Genet. 2018 Sep;34(9):666–81. https://doi.org/10.1016/j.tig.2018.05.008. Epub 2018 Jun 22. PMID: 29941292.

  17. Daniszewski M, Senabouth A, Nguyen QH, Crombie DE, Lukowski SW, Kulkarni T, Sluch VM, Jabbari JS, Chamling X, Zack DJ, Pébay A, Powell JE, Hewitt AW. Single cell RNA sequencing of stem cell-derived retinal ganglion cells. Sci Data 2018 Feb 13;5:180013. doi: https://doi.org/10.1038/sdata.2018.13. PMID: 29437159; PMCID: PMC5810423.

  18. Li CL, Li KC, Wu D, Chen Y, Luo H, Zhao JR, Wang SS, Sun MM, Lu YJ, Zhong YQ, Hu XY, Hou R, Zhou BB, Bao L, Xiao HS, Zhang X. Somatosensory neuron types identified by high-coverage single-cell RNA-sequencing and functional heterogeneity. Cell Res. 2016 Jan;26(1):83–102. doi: https://doi.org/10.1038/cr.2015.149. Epub 2015 Dec 22. Erratum in: Cell Res. 2016 Aug;26(8):967. PMID: 26691752; PMCID: PMC4816136.

  19. Sklavenitis-Pistofidis R, Getz G, Ghobrial I. Single-cell RNA sequencing: one step closer to the clinic. Nat Med. 2021 Mar;27(3):375–376. doi: https://doi.org/10.1038/s41591-021-01276-y. PMID: 33664491.

  20. Zhang X, Li T, Liu F, Chen Y, Yao J, Li Z, Huang Y, Wang J. Comparative Analysis of Droplet-Based Ultra-High-Throughput Single-Cell RNA-Seq Systems. Mol Cell. 2019 Jan 3;73(1):130–142.e5. doi: https://doi.org/10.1016/j.molcel.2018.10.020. Epub 2018 Nov 21. PMID: 30472192.

  21. Pereira WJ, Almeida FM, Conde D, Balmant KM, Triozzi PM, Schmidt HW, Dervinis C, Pappas GJ Jr, Kirst M. Asc-Seurat: analytical single-cell Seurat-based web application. BMC Bioinformatics. 2021 Nov 18;22(1):556. doi: https://doi.org/10.1186/s12859-021-04472-2. PMID: 34794383; PMCID: PMC8600690.

  22. Gene Ontology Consortium. Gene Ontology Consortium: going forward. Nucleic Acids Res. 2015 Jan;43(Database issue):D1049-56. doi: https://doi.org/10.1093/nar/gku1179. Epub 2014 Nov 26. PMID: 25428369; PMCID: PMC4383973.

  23. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000 Jan 1;28(1):27–30. doi: https://doi.org/10.1093/nar/28.1.27. PMID: 10592173; PMCID: PMC102409.

  24. Feigin VL, Brainin M, Norrving B, Martins S, Sacco RL, Hacke W, Fisher M, Pandian J, Lindsay P. World Stroke Organization (WSO): Global Stroke Fact Sheet 2022. Int J Stroke. 2022 Jan;17(1):18–29. doi: 10.1177/17474930211065917. Erratum in: Int J Stroke. 2022 Apr;17(4):478. PMID: 34986727.

  25. Fluss R, Ligas BA, Chan AW, Ellis JA, Ortiz RA, Langer DJ, Rahme R. Moyamoya-Related stroke risk during pregnancy: an evidence-based reappraisal. World Neurosurg. 2019 Sep;129:e582–5. https://doi.org/10.1016/j.wneu.2019.05.225. Epub 2019 Jun 6. PMID: 31176061.

  26. Fujimura M, Bang OY, Kim JS. Moyamoya Disease. Front Neurol Neurosci. 2016;40:204–220. doi: 10.1159/000448314. Epub 2016 Dec 2. PMID: 27960175.

  27. Gatti JR, Torriente AG, Sun LR. Clinical presentation and stroke incidence Differ by Moyamoya etiology. J Child Neurol. 2021 Mar;36(4):272–80. https://doi.org/10.1177/0883073820967160. Epub 2020 Nov 6. PMID: 33155871.

Download references

Acknowledgements

We thank the GEO database.

Funding

This work was supported by the Jiangsu Province Hospital(the First Affiliated Hospital with Nanjing Medical University) Clinical Capacity Enhancement Project (JSPH-MA-2022-8) and Research Project of Jiangsu Health Commission (No.ZD2022054).

Author information

Authors and Affiliations

Authors

Contributions

QT was responsible for the design of this study. QT, WL, JH, CM, YT and QZ were responsible for all data analysis and manuscript writing. QT, YW, JL, JX, YL and XM were responsible for writing the manuscript. WW provided funding. All the authors agreed on the final version of the article.

Corresponding author

Correspondence to Wei Wu.

Ethics declarations

Ethics approval and consent to participate

All procedures performed were in accordance with the declaration of the ethical standards of the institutional research committee and with the 1964 Helsinki 387 Declaration and its later amendments. The ethics committee has approved this study of the First Affiliated Hospital of Nanjing Medical University(No. 2021-SR-165).

Consent for publication

All authors agreed for publication.

Competing interests

All authors declare that no conflict of interest exists.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tang, Q., Li, W., Huang, J. et al. Single-cell sequencing analysis of peripheral blood in patients with moyamoya disease. Orphanet J Rare Dis 18, 174 (2023). https://doi.org/10.1186/s13023-023-02781-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13023-023-02781-8

Keywords