Elucidating the role of tumor-associated ALOX5+ mast cells with transformative function in cervical cancer progression via single-cell RNA sequencing.

Cervical cancer (CC) is the fourth most common malignancy among women globally and serves as the main cause of cancer-related deaths among women in developing countries. The early symptoms of CC are often not apparent, with diagnoses typically made at advanced stages, which lead to poor clinical prognoses. In recent years, numerous studies have shown that there is a close relationship between mast cells (MCs) and tumor development. However, research on the role MCs played in CC is still very limited at that time. Thus, the study conducted a single-cell multi-omics analysis on human CC cells, aiming to explore the mechanisms by which MCs interact with the tumor microenvironment in CC. The goal was to provide a scientific basis for the prevention, diagnosis, and treatment of CC, with the hope of improving patients' prognoses and quality of life. The present study acquired single-cell RNA sequencing data from ten CC tumor samples in the ArrayExpress database. Slingshot and AUCcell were utilized to infer and assess the differentiation trajectory and cell plasticity of MCs subpopulations. Differential expression analysis of MCs subpopulations in CC was performed, employing Gene Ontology, gene set enrichment analysis, and gene set variation analysis. CellChat software package was applied to predict cell communication between MCs subpopulations and CC cells. Cellular functional experiments validated the functionality of TNFRSF12A in HeLa and Caski cell lines. Additionally, a risk scoring model was constructed to evaluate the differences in clinical features, prognosis, immune infiltration, immune checkpoint, and functional enrichment across various risk scores. Copy number variation levels were computed using inference of copy number variations. The obtained 93,524 high-quality cells were classified into ten cell types, including T_NK cells, endothelial cells, fibroblasts, smooth muscle cells, epithelial cells, B cells, plasma cells, MCs, neutrophils, and myeloid cells. Furthermore, a total of 1,392 MCs were subdivided into seven subpopulations: C0 CTSG+ MCs, C1 CALR+ MCs, C2 ALOX5+ MCs, C3 ANXA2+ MCs, C4 MGP+ MCs, C5 IL32+ MCs, and C6 ADGRL4+ MCs. Notably, the C2 subpopulation showed close associations with tumor-related MCs, with Slingshot results indicating that C2 subpopulation resided at the intermediate-to-late stage of differentiation, potentially representing a crucial transition point in the benign-to-malignant transformation of CC. CNVscore and bulk analysis results further confirmed the transforming state of the C2 subpopulation. CellChat analysis revealed TNFRSF12A as a key receptor involved in the actions of C2 ALOX5+ MCs. Moreover, in vitro experiments indicated that downregulating the TNFRSF12A gene may partially inhibit the development of CC. Additionally, a prognosis model and immune infiltration analysis based on the marker genes of the C2 subpopulation provided valuable guidance for patient prognosis and clinical intervention strategies. We first identified the transformative tumor-associated MCs subpopulation C2 ALOX5+ MCs within CC, which was at a critical stage of tumor differentiation and impacted the progression of CC. In vitro experiments confirmed the inhibitory effect of knocking down the TNFRSF12A gene on the development of CC. The prognostic model constructed based on the C2 ALOX5+MCs subset demonstrated excellent predictive value. These findings offer a fresh perspective for clinical decision-making in CC.

Zhao F ，Hong J ，Zhou G ，Huang T ，Lin Z ，Zhang Y ，Liang L ，Tang H ... -  《Frontiers in Immunology》

Single-cell analysis uncovers high-proliferative tumour cell subtypes and their interactions in the microenvironment of gastric cancer.

Gastric cancer (GC) remains a prominent malignancy that poses a significant threat to human well-being worldwide. Despite advancements in chemotherapy and immunotherapy, which have effectively augmented patient survival rates, the mortality rate associated with GC remains distressingly high. This can be attributed to the elevated proliferation and invasive nature exhibited by GC. Our current understanding of the drivers behind GC cell proliferation remains limited. Hence, in order to reveal the molecular biological mechanism behind the swift advancement of GC, we employed single-cell RNA-sequencing (scRNA-seq) to characterize the tumour microenvironment in this study. The scRNA-seq data of 27 patients were acquired from the Gene Expression Omnibus database. Differential gene analysis, Gene Ontology, Kyoto Encyclopedia of Genes and Genomes and Gene Set Enrichment Analysis were employed to investigate 38 samples. The copy number variation level exhibited by GC cells was determined using InferCNV. The CytoTRACE, Monocle and Slingshot analysis were used to discern the cellular stemness and developmental trajectory of GC cells. The CellChat package was utilized for the analysis of intercellular communication crosstalk. Moreover, the findings of the data analysis were validated through cellular functional tests conducted on the AGS cell line and SGC-7901 cell line. Finally, this study constructed a risk scoring model to evaluate the differences of different risk scores in clinical characteristics, immune infiltration, immune checkpoints, functional enrichment, tumour mutation burden and drug sensitivity. Within the microenvironment of GC, we identified the presence of 8 cell subsets, encompassing NK_T cells, B_Plasma cells, epithelial cells, myeloid cells, endothelial cells, mast cells, fibroblasts, pericytes. By delving deeper into the characterization of GC cells, we identified 6 specific tumour cell subtypes: C0 PSCA+ tumour cells, C1 CLDN7+ tumour cells, C2 UBE2C+ tumour cells, C3 MUC6+ tumour cells, C4 CHGA+ tumour cells and C5 MUC2+ tumour cells. Notably, the C2 UBE2C+ tumour cells demonstrated a close association with cell mitosis and the cell cycle, exhibiting robust proliferative capabilities. Our findings were fortified through enrichment analysis, pseudotime analysis and cell communication analysis. Meanwhile, knockdown of the transcription factor CREB3, which is highly active in UBE2C+ tumour cells, significantly impedes the proliferation, migration and invasion of GC cells. And the prognostic score model constructed with CREB3-related genes showcased commendable clinical predictive capacity, thus providing valuable guidance for patients' prognosis and clinical treatment decisions. We have identified a highly proliferative cellular subgroup C2 UBE2C+ tumour cells in GC for the first time. The employment of a risk score model, which is based on genes associated with UBE2C expression, exhibits remarkable proficiency in predicting the prognosis of GC patients. In our investigation, we observed that the knockdown of the transcription factor CREB3 led to a marked reduction in cellular proliferation, migration and invasion in GC cell line models. Implementing a stratified treatment approach guided by this model represents a judicious and promising methodology.

Zhang W ，Wang X ，Dong J ，Wang K ，Jiang W ，Fan C ，Liu H ，Fan L ，Zhao L ，Li G ... -  《-》

Unraveling the ecological landscape of mast cells in esophageal cancer through single-cell RNA sequencing.

Esophageal cancer (EC) is a major health issue, ranking seventh in incidence and sixth in mortality worldwide. Despite advancements in multidisciplinary treatment approaches, the 5-year survival rate for EC remains low at 21%. Challenges in EC treatment arise from late-stage diagnosis, high malignancy, and poor prognosis. Understanding the tumor microenvironment is critical, as it includes various cellular and extracellular components that influence tumor behavior and treatment response. Mast cells (MCs), as tissue-resident immune cells, play dual roles in tumor dynamics. High-throughput single-cell RNA sequencing offers a powerful tool for analyzing tumor heterogeneity and immune interactions, although its application in EC is limited. In this study, we investigated the immune microenvironment of EC using single-cell RNA sequencing and established a comprehensive immune profile. We also performed analysis of upstream transcription factors and downstream pathway enrichment to further comprehensively decipher MCs in EC. Besides, we performed knockdown experiments to explore the role of epidermal growth factor receptor (EGFR) signaling pathway in MCs-tumor cell interactions, highlighting its potential as a prognostic marker. Finally, we constructed a prognostic model for EC, which provided valuable suggestions for the diagnosis and prognosis of EC. Our analysis identified 11 major cell types, of which MCs were particularly present in pericarcinoma tissues. Further grouping of the 5,001 MCs identified 8 distinct subtypes, including SRSF7-highly expressed MCs, which showed strong tumor preference and potential tumor-promoting properties. Moreover, we identified the key signaling receptor EGFR and validated it by in vitro knockdown experiments, demonstrating its cancer-promoting effects. In addition, we established an independent prognostic indicator, SRSF7+ MCs risk score (SMRS), which showed a correlation between high SMRS group and poor prognosis. These findings illuminate the complex interactions within the tumor microenvironment of EC and suggest that targeting specific MCs subtypes, particularly via the EGFR signaling pathway, may present novel therapeutic strategies. This study establishes a comprehensive immune map of EC, offering insights for improved treatment approaches.

Zhang S ，Zhang X ，Xiahou Z ，Zuo S ，Xue J ，Zhang Y ... -  《-》

HPV-driven heterogeneity in cervical cancer: study on the role of epithelial cells and myofibroblasts in the tumor progression based on single-cell RNA sequencing analysis.

Cervical cancer (CC) is a neoplasia with a high heterogeneity. We aimed to explore the characteristics of tumor microenvironment (TME) for CC treatment. HPV positive (+) and negative (-) samples from cervical cancer (CC) patients were sourced from the Gene Expression Omnibus (GEO) database. The single-cell RNA sequencing (scRNA-seq) data were processed and annotated for cell types utilizing the Seurat package. Following this, the expression levels and biological roles of the marker genes were analyzed applying real-time PCR (RT-PCR) and transwell assays. Furthermore, the enrichment of genes with significantly differential expressions and copy number variations was assessed by the ClusterProlifer and inferCNV software packages. Seven main cell clusters were classified based on a total of 12,431 cells. The HPV- CC samples exhibited a higher immune cell infiltration level, while epithelial cells and myofibroblasts had higher proportion in the HPV+ CC samples with extensive heterogeneity. Immune pathways including antigen treatment and presentation, immunoglobulin production and T cell mediated immunity were significantly activated in the HPV- CC group with lower cell cycle and proliferation activity. However, the anti-tumor immunity of these cells was inhibited in HPV+ CC group with higher cell proliferation activity. Moreover, the amplification and loss of CNVs also supported that these cells in HPV- CC samples were prone to anti-tumor activation. Further cell validation results showed that except GZMA, the levels of APOC1, CEACAM6, FOXP3, SFRP4 and TFF3 were all higher in CC cells Hela, and that silencing TFF3 could inhibit the migration and invasion of CC cells in-vitro. This study highlighted the critical role of HPV infection in CC progression, providing a novel molecular basis for optimizing the current preventive screening and personalized treatment for the cancer.

Zhang Y ，Zhang Y ，Pan C ，Wang W ，Yu Y ... -  《-》

Single-cell RNA sequencing reveals that MYBL2 in malignant epithelial cells is involved in the development and progression of ovarian cancer.

Ovarian carcinoma (OC) is a prevalent gynecological malignancy associated with high recurrence rates and mortality, often diagnosed at advanced stages. Despite advances in immunotherapy, immune exhaustion remains a significant challenge in achieving optimal tumor control. However, the exploration of intratumoral heterogeneity of malignant epithelial cells and the ovarian cancer tumor microenvironment is still limited, hindering our comprehensive understanding of the disease. Utilizing single-cell RNA sequencing (scRNA-seq), we comprehensively investigated the cellular composition across six ovarian cancer patients with omental metastasis. Our focus centered on analysis of the malignant epithelial cells. Employing CytoTRACE and slingshot pseudotime analyses, we identified critical subpopulations and explored associated transcription factors (TFs) influencing ovarian cancer progression. Furthermore, by integrating clinical factors from a large cohort of bulk RNA sequencing data, we have established a novel prognostic model to investigate the impact of the tumor immune microenvironment on ovarian cancer patients. Furthermore, we have investigated the condition of immunological exhaustion. Our study identified a distinct and highly proliferative subgroup of malignant epithelial cells, known as C2 TOP2A+ TCs. This subgroup primarily consisted of patients who hadn't received neoadjuvant chemotherapy. Ovarian cancer patients with elevated TOP2A expression exhibited heightened sensitivity to neoadjuvant chemotherapy (NACT). Moreover, the transcription factor MYBL2 in this subgroup played a critical role in ovarian cancer development. Additionally, we developed an independent prognostic indicator, the TOP2A TCs Risk Score (TTRS), which revealed a correlation between the High TTRS Group and unfavorable outcomes. Furthermore, immune infiltration and drug sensitivity analyses demonstrated increased responsiveness to Paclitaxel, Cisplatin, and Gemcitabine in the Low TTRS Group. This research deepens our understanding of malignant epithelial cells in ovarian cancer and enhances our knowledge of the ovarian cancer immune microenvironment and immune exhaustion. We have revealed the heightened susceptibility of the C2 TOP2A+ TCs subgroup to neoadjuvant chemotherapy and emphasized the role of MYBL2 within the C2 subgroup in promoting the occurrence and progression of ovarian cancer. These insights provide valuable guidance for the management of ovarian cancer treatment.

Shao W ，Lin Z ，Xiahou Z ，Zhao F ，Xu J ，Liu X ，Cai P ... -  《Frontiers in Immunology》