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

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.