Dataset Information

Integration of scRNA-Seq and bulk RNA-Seq uncover perturbed immune cell types and pathways of Kawasaki disease.

ABSTRACT:

Introduction

Kawasaki disease (KD) is an acute febrile illness primarily affecting children and characterized by systemic inflammation and vasculitis that can lead to coronary artery complications. The aim of this study was to gain a comprehensive understanding of immune dysregulation in KD.

Methods

To this end, we employed integration of single-cell RNA sequencing (scRNA-Seq) and bulk RNA sequencing (bulk RNA-Seq) data. Furthermore, we conducted flow cytometry analysis for a cohort of 82 KD patients.

Results

Our analysis revealed significant heterogeneity within immune cell populations in KD patients, with distinct clusters of T cells, B cells, and natural killer (NK) cells. Importantly, CD4⁺ naïve T cells in KD patients were found to predominantly differentiate into Treg cells and Th2 cells, potentially playing a role in the excessive inflammation and vascular damage characteristic of the disease. Dysregulated signaling pathways were also identified, including the mTOR signaling pathway, cardiomyopathy pathway, COVID-19 signaling pathway, and pathways involved in bacterial or viral infection.

Discussion

These findings provide insights into the immunopathogenesis of KD, emphasizing the importance of immune cell dysregulation and dysregulated signaling pathways. Integration of scRNA-Seq and bulk RNA-Seq data offers a comprehensive view of the molecular and cellular alterations in KD and highlights potential therapeutic targets for further investigation. Validation and functional studies are warranted to elucidate the roles of the identified immune cell types and pathways in KD pathogenesis and to develop targeted interventions to improve patient outcomes.

SUBMITTER: Cao N

PROVIDER: S-EPMC10568768 | biostudies-literature | 2023

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Integration of scRNA-Seq and bulk RNA-Seq uncover perturbed immune cell types and pathways of Kawasaki disease.

Cao Naixin N Ouyang Huayi H Zhang Xing X Xu Yuanyuan Y Li Jun J Chen Yanfei Y

Frontiers in immunology 20230928

<h4>Introduction</h4>Kawasaki disease (KD) is an acute febrile illness primarily affecting children and characterized by systemic inflammation and vasculitis that can lead to coronary artery complications. The aim of this study was to gain a comprehensive understanding of immune dysregulation in KD.<h4>Methods</h4>To this end, we employed integration of single-cell RNA sequencing (scRNA-Seq) and bulk RNA sequencing (bulk RNA-Seq) data. Furthermore, we conducted flow cytometry analysis for a coho ...[more]

PMID: 37841239

Similar Datasets

Project description:BackgroundColorectal cancer is a malignant tumor of the digestive system originating from abnormal cell proliferation in the colon or rectum, often leading to gastrointestinal symptoms and severe health issues. Nucleotide metabolism, which encompasses the synthesis of DNA and RNA, is a pivotal cellular biochemical process that significantly impacts both the progression and therapeutic strategies of colorectal cancer METHODS: For single-cell RNA sequencing (scRNA-seq), five functions were employed to calculate scores related to nucleotide metabolism. Cell developmental trajectory analysis and intercellular interaction analysis were utilized to explore the metabolic characteristics and communication patterns of different epithelial cells. These findings were further validated using spatial transcriptome RNA sequencing (stRNA-seq). A risk model was constructed using expression profile data from TCGA and GEO cohorts to optimize clinical decision-making. Key nucleotide metabolism-related genes (NMRGs) were functionally validated by further in vitro experiments.ResultsIn both scRNA-seq and stRNA-seq, colorectal cancer (CRC) exhibited unique cellular heterogeneity, with myeloid cells and epithelial cells in tumor samples displaying higher nucleotide metabolism scores. Analysis of intercellular communication revealed enhanced signaling pathways and ligand-receptor interactions between epithelial cells with high nucleotide metabolism and fibroblasts. Spatial transcriptome sequencing confirmed elevated nucleotide metabolism states in the core region of tumor tissue. After identifying differentially expressed NMRGs in epithelial cells, a risk prognostic model based on four genes effectively predicted overall survival and immunotherapy outcomes in patients. High-risk group patients exhibited an immunosuppressive microenvironment and relatively poorer prognosis and responses to chemotherapy and immunotherapy. Finally, based on data analysis and a series of cellular functional experiments, ACOX1 and CPT2 were identified as novel therapeutic targets for CRC.ConclusionIn this study, a comprehensive analysis of NMRGs in CRC was conducted using a combination of single-cell sequencing, spatial transcriptome sequencing, and high-throughput data. The prognostic model constructed with NMRGs shows potential as a standalone prognostic marker for colorectal cancer patients and may significantly influence the development of personalized treatment approaches for CRC.

Project description:Lung adenocarcinoma (LUAD) is the most common type of lung cancer, accounting for approximately 35-40% of lung cancers, and the overall survival time of patients with LUAD is still very poor. B cells are important effector cells of adaptive immunity, and B-cell infiltration increases in various tumors. The role of B cells in LUAD is still largely unknown. Therefore, it is particularly important to clarify the role of B cells in LUAD. GSE164983, GSE50081, GSE37745 and GSE30219 were obtained from the GEO database. The TCGA-LUAD dataset was obtained from the TCGA database. UMAP was used to perform clustering descending and subgroup identification on single-cell RNA-sequencing (scRNA-seq) data to obtain B-cell markers. The TCGA cohort was used to obtain differentially expressed genes (DEGs). B-cell-related differentially expressed genes (BRGs) were identified through the intersection of B-cell markers and DEGs. The LASSO method was used to identify characteristic genes of BRGs and construct a prognostic risk model. LUAD patients were divided into high-risk and low-risk groups based on risk scores, and the immune landscape of the two groups was evaluated. We also analyzed the differences in clinical characteristics, mutations, immunotherapy, and drug sensitivity between the two groups. Thirty BRGs were obtained, and 6 characteristic genes were identified. Based on the characteristic genes, a prognostic risk model was constructed. According to the prognostic risk model, LUAD patients were divided into two groups: high-risk group and low-risk group. Patients in the high-risk group had worse outcomes and shorter survival times. Low-risk patients had better survival, while patients with high TNM stage accounted for a greater proportion of patients in the high-risk group. In addition, high-risk patients had a greater probability of mutation and worse immunotherapy response. Finally, we found different susceptibility profiles between the high-risk and low-risk groups. The prognostic risk model built based on the BRGs had good predictive performance, providing a new perspective on the prognosis and immunotherapy of LUAD patients and a new reference for LUAD research.

Project description:BackgroundConsiderable evidence suggests that the heterogeneity of ovarian cancer (OC) is a major cause of treatment failure. Single-cell RNA sequencing (scRNA-seq) is a powerful tool to analyse the heterogeneity of the tumour at the single-cell level, leading to a better understanding of cell function at the genetic and cellular levels.MethodsOC scRNA-seq data were extracted from the Gene Expression Omnibus (GEO) database and the FindCluster () package used for cell cluster analysis. The GSVA package was used for single-sample gene set enrichment analysis (ssGSEA) analysis to obtain a Hallmark gene set score and bulk RNA-seq data were used to analyse the key genes of OC-associated immune cell subsets. CIBERSORT was used to identify immune scores of cells and the "WGCNA" package for the weighted correlation network analysis (WGCNA). KEGG (Kyoto Encyclopedia of Genes and Genomes) and GO (Gene Ontology) analyses of subtype groups were performed by GSEA. Then, univariate Cox and lasso regression were performed to further establish a signature. Finally, qPCR and immunohistochemistry staining were used to evaluate the expression of signature genes in OC.ResultsTwo scRNA-seq (GSE154600 and GES158937) datasets were integrated to obtain 20 cell clusters. T cells or NK cells (cluster 5, 6, 7, 11), B cells (cluster 16, 19, 20) and myeloid cells (cluster 4, 9, 10) were clustered according to immune cell markers. The ssGSEA revealed that M1- and M2-like myeloid cell-related genes were significantly upregulated in P3 and P4 patients in the GSE154600 data. Immune cell analysis in TCGA-OC showed that a high abundance of M1-like tumour-associated macrophages (TAMS) predicts better survival. WGCNA, univariate Cox and lasso Cox regression established a two-gene signature (RiskScore=-0.059*CXCL13-0.034*IL26). Next, the TCGA-test and TCGA-OC were used to test the risk prediction ability of the signature, showing a good effect in the datasets. Moreover, the qPCR and immunohistochemistry staining revealed that the expression of CXCL13 and IL26 was reduced in OC tissues.ConclusionA two-gene signature prognostic stratification system (CXCL13 and IL26) was developed based on the heterogeneity of OC immune cells to accurately evaluate the prognostic risk.

Project description:BackgroundRecent studies have suggested that cell death may be involved in bone loss or the resolution of inflammation in periodontitis. Immunogenic cell death (ICD), a recently identified cell death pathway, may be involved in the development of this disease.MethodsBy analyzing single-cell RNA sequencing (scRNA-seq) for periodontitis and scoring gene set activity, we identified cell populations associated with ICD, which were further verified by qPCR, enzyme linked immunosorbent assay (ELISA) and immunofluorescence (IF) staining. By combining the bulk transcriptome and applying machine learning methods, we identified several potential ICD-related hub genes, which were then used to build diagnostic models. Subsequently, consensus clustering analysis was performed to identify ICD-associated subtypes, and multiple bioinformatics algorithms were used to investigate differences in immune cells and pathways between subtypes. Finally, qPCR and immunohistochemical staining were performed to validate the accuracy of the models.ResultsSingle-cell gene set activity analysis found that in non-immune cells, fibroblasts had a higher ICD activity score, and KEGG results showed that fibroblasts were enriched in a variety of ICD-related pathways. qPCR, Elisa and IF further verified the accuracy of the results. From the bulk transcriptome, we identified 11 differentially expressed genes (DEGs) associated with ICD, and machine learning methods further identified 5 hub genes associated with ICD. Consensus cluster analysis based on these 5 genes showed that there were differences in immune cells and immune functions among subtypes associated with ICD. Finally, qPCR and immunohistochemistry confirmed the ability of these five genes as biomarkers for the diagnosis of periodontitis.ConclusionFibroblasts may be the main cell source of ICD in periodontitis. Adaptive immune responses driven by ICD may be one of the pathogenesis of periodontitis. Five key genes associated with ICD (ENTPD1, TLR4, LY96, PRF1 and P2RX7) may be diagnostic biomarkers of periodontitis and future therapeutic targets.

Dataset Information

Integration of scRNA-Seq and bulk RNA-Seq uncover perturbed immune cell types and pathways of Kawasaki disease.

Introduction

Methods

Results

Discussion

Publications

Integration of scRNA-Seq and bulk RNA-Seq uncover perturbed immune cell types and pathways of Kawasaki disease.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets