Single-cell Sequencing Research Articles

Noncaloric monosaccharides induce excessive sprouting angiogenesis in zebrafish via foxo1a-marcksl1a signal.

Artificially sweetened beverages containing noncaloric monosaccharides were suggested as healthier alternatives to sugar-sweetened beverages. Nevertheless, the potential detrimental effects of these noncaloric monosaccharides on blood vessel function remain inadequately understood. We have established a zebrafish model that exhibits significant excessive angiogenesis induced by high glucose, resembling the hyperangiogenic characteristics observed in proliferative diabetic retinopathy (PDR). Utilizing this model, we observed that glucose and noncaloric monosaccharides could induce excessive formation of blood vessels, especially intersegmental vessels (ISVs). The excessively branched vessels were observed to be formed by ectopic activation of quiescent endothelial cells (ECs) into tip cells. Single-cell transcriptomic sequencing analysis of the ECs in the embryos exposed to high glucose revealed an augmented ratio of capillary ECs, proliferating ECs, and a series of upregulated proangiogenic genes. Further analysis and experiments validated that reduced foxo1a mediated the excessive angiogenesis induced by monosaccharides via upregulating the expression of marcksl1a. This study has provided new evidence showing the negative effects of noncaloric monosaccharides on the vascular system and the underlying mechanisms.

Multimodal Identification of Molecular Factors Linked to Severe Diabetic Foot Ulcers Using Artificial Intelligence.

Diabetic foot ulcers (DFUs) are a severe complication of diabetes mellitus (DM), which often lead to hospitalization and non-traumatic amputations in the United States. Diabetes prevalence estimates in South Texas exceed the national estimate and the number of diagnosed cases is higher among Hispanic adults compared to their non-Hispanic white counterparts. San Antonio, a predominantly Hispanic city, reports significantly higher annual rates of diabetic amputations compared to Texas. The late identification of severe foot ulcers minimizes the likelihood of reducing amputation risk. The aim of this study was to identify molecular factors related to the severity of DFUs by leveraging a multimodal approach. We first utilized electronic health records (EHRs) from two large demographic groups, encompassing thousands of patients, to identify blood tests such as cholesterol, blood sugar, and specific protein tests that are significantly associated with severe DFUs. Next, we translated the protein components from these blood tests into their ribonucleic acid (RNA) counterparts and analyzed them using public bulk and single-cell RNA sequencing datasets. Using these data, we applied a machine learning pipeline to uncover cell-type-specific and molecular factors associated with varying degrees of DFU severity. Our results showed that several blood test results, such as the Albumin/Creatinine Ratio (ACR) and cholesterol and coagulation tissue factor levels, correlated with DFU severity across key demographic groups. These tests exhibited varying degrees of significance based on demographic differences. Using bulk RNA-Sequenced (RNA-Seq) data, we found that apolipoprotein E (APOE) protein, a component of lipoproteins that are responsible for cholesterol transport and metabolism, is linked to DFU severity. Furthermore, the single-cell RNA-Seq (scRNA-seq) analysis revealed a cluster of cells identified as keratinocytes that showed overexpression of APOE in severe DFU cases. Overall, this study demonstrates how integrating extensive EHRs data with single-cell transcriptomics can refine the search for molecular markers and identify cell-type-specific and molecular factors associated with DFU severity while considering key demographic differences.

CDT1 is a Potential Therapeutic Target for the Progression of NAFLD to HCC and the Exacerbation of Cancer

Aims: This study aimed to identify potential therapeutic targets in the progression from non-alcoholic fatty liver disease (NAFLD) to hepatocellular carcinoma (HCC), with a focus on genes that could influence disease development and progression. Background: Hepatocellular carcinoma, significantly driven by non-alcoholic fatty liver disease, represents a major global health challenge due to late-stage diagnosis and limited treatment options. This study utilizes bioinformatics to analyze data from GEO and TCGA, aiming to uncover molecular biomarkers that bridge NAFLD to HCC. Through identifying critical genes and pathways, our research seeks to advance early detection and develop targeted therapies, potentially improving prognosis and personalizing treatment for NAFLD-HCC patients. Objective: Identify key genes that differ between NAFLD and HCC; Analyze these genes to understand their roles in disease progression; Validate the functions of these genes in NAFLD to HCC transition. Methods: Initially, we identified a set of genes differentially expressed in both NAFLD and HCC using second-generation sequencing data from the GEO and TCGA databases. We then employed a Cox proportional hazards model and a Lasso regression model, applying machine learning techniques to the large sample data from TCGA. This approach was used to screen for key disease-related genes, and an external dataset was utilized for model validation. Additionally, pseudo-temporal sequencing analysis of single-cell sequencing data was performed to further examine the variations in these genes in NAFLD and HCC. Results: The machine learning analysis identified IGSF3, CENPW, CDT1, and CDC6 as key genes. Furthermore, constructing a machine learning model for CDT1 revealed it to be the most critical gene, with model validation yielding an ROC value greater than 0.80. The single-cell sequencing data analysis confirmed significant variations in the four predicted key genes between the NAFLD and HCC groups. Conclusion: Our study underscores the pivotal role of CDT1 in the progression from NAFLD to HCC. This finding opens new avenues for early diagnosis and targeted therapy of HCC, highlighting CDT1 as a potential therapeutic target.

Identification of a seven-gene prognostic model for renal cell carcinoma associated with CD8+T lymphocyte cell.

CD8+ T lymphocytes are important elements of the tumor microenvironment, hence their involvement in the development and progression of tumors is complex. Data on the precise tumor-infiltrating lymphocytes gene signature in renal cell carcinoma (RCC) remain limited. Therefore, this study created a tumor-infiltrating lymphocytes-related predictive model for patients with RCC using data from The Cancer Genome Atlas. The most important genes associated with CD8 + T lymphocytes were identified using weighted gene co-expression network analysis. Functional categories of important genes were revealed using gene ontology enrichment and Kyoto Encyclopedia of Genes and Genomes signaling pathway analyses. A CD8 + T lymphocyte-related prognostic model with 7 important genes was simultaneously created using the least absolute shrinkage and selection operator, univariate and multivariate Cox regressions, and the 7 genes were expressed particularly in CD8 + T lymphocytes according to single-cell sequencing data obtained from the Gene Expression Omnibus. This study identified a seven-gene prognostic model associated with CD8 + T lymphocytes that may significantly influence risk stratification in patients with RCC. The genes included in the model are apolipoprotein B mRNA editing catalytic polypeptide 3G, CD3 gamma, eomesodermin, protein tyrosine phosphatase, non-receptor type 7, signal regulatory protein gamma, Fas ligand, and T-cell immunoreceptor with Ig and ITIM domains.

NitraTh epitope-based neoantigen vaccines for effective tumor immunotherapy

Neoantigen vaccines represent an emerging and promising strategy in the field of tumor immunotherapy. Despite their potential, designing an effective neoantigen vaccine remains a challenge due to the current limitations in predicting CD4+ T cell epitopes with high accuracy. Here, we introduce a novel approach to neoantigen vaccine design that does not rely on computational prediction of CD4+ T cell epitopes. Utilizing nitrated helper T cell epitope containing p-nitrophenylalanine, termed “NitraTh epitope,” we have successfully engineered a series of tumor neoantigen vaccines capable of eliciting robust neoantigen-specific immune responses. With the help of NitraTh epitope, even mutations with low predicted affinity for MHC class I molecules were successfully induced to elicit neoantigen-specific responses. In H22 cell allograft and patient-derived xenograft (PDX) liver cancer mouse models, the NitraTh epitope-based neoantigen vaccines significantly suppressed tumor progression. More strikingly, through single-cell sequencing we found that the NitraTh epitope-based neoantigen vaccines regulate macrophage reprogramming and modulate macrophages to decrease the levels of the immunosuppressive molecule prostaglandin E2 (PGE2), which in turn reshapes the tumor immunosuppressive microenvironment. In summary, NitraTh epitope-based neoantigen vaccines possess the dual effects of potently activating neoantigen-specific immunity and alleviating immunosuppression, potentially providing a new paradigm for the design of tumor neoantigen vaccines.

CIEC: Cross-tissue Immune Cell Type Enrichment and Expression Map Visualization for Cancer.

Single-cell transcriptome sequencing technology has been applied to decode the cell types and functional states of immune cells, revealing their tissue-specific gene expression patterns and functions in cancer immunity. Comprehensive assessments of immune cells within and across tissues will provide us with a deeper understanding of the tumor immune system in general. Here, we present Cross-tissue Immune cell type or state Enrichment analysis of gene lists for Cancer (CIEC), the first web-based application that integrates database and enrichment analysis to estimate the cross-tissue immune cell type or state. CIEC version 1.0 consists of 480 samples covering primary tumor, adjacent normal tissue, lymph node, metastasis tissue, and peripheral blood from 323 cancer patients. By applying integrative analysis, we constructed an immune cell-type/state map for each context and adopted our previously developed Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthology Based Annotation System (KOBAS) algorithm to estimate the enrichment for context-specific immune cell type/state. In addition, CIEC also provides an easy-to-use online interface for users to comprehensively analyze the immune cell characteristics mapped across multiple tissues, including expression map, correlation, similar genes detection, signature score, and expression comparison. We believe that CIEC will be a valuable resource for exploring the intrinsic characteristics of immune cells in cancer patients and for potentially guiding novel cancer-immune biomarker development and immunotherapy strategies. CIEC is freely accessible at http://ciec.gene.ac/.

Screening and validating genes associated with cuproptosis in systemic lupus erythematosus by expression profiling combined with machine learning.

Cell death has long been a focal point in life sciences research, and recently, scientists have discovered a novel form of cell death induced by copper, termed cuproptosis. This paper aimed to identify genes associated with cuproptosis in Systemic Lupus Erythematosus (SLE) through machine learning, combined with single-cell sequencing (scRNA-seq), to screen and validate related genes. The analytical results were then experimentally verified. Two published microarray gene expression datasets (GSE65391 and GSE61635) from SLE and control peripheral blood samples were downloaded from the GEO database. The GSE65391 dataset was used as the training group, while the GSE61635 dataset served as the validation group. Differentially expressed genes (DEGs) from GSE65391 identified 12 differential genes. Nine diagnostic genes, considered potential biomarkers, were selected using the least absolute shrinkage and selection operator (LASSO) and support vector machine recursive feature elimination (SVM-RFE) analysis. The Receiver Operating Characteristic (ROC) curves for both the training and validation groups were used to calculate the Area Under the Curve (AUC) to assess discriminatory properties. CIBERSORT was used to assess the relationship between these diagnostic genes and a reference set of infiltrating immune cells. Single-cell RNA sequencing data (GSE162577) from SLE patients were also obtained from the GEO database and analyzed. Experimental validation of the most important SLE biomarkers was performed. Twelve significantly different cuproptosis-related genes were identified in the GSE65391 training set. Immune cell analysis revealed 12 immune cell types and identified nine signature genes, including PDHB, GLS, DLAT, LIAS, MTF1, DLST, DLD, LIPT1, and FDX1. In the GSE61635 validation set, seven genes were weakly expressed, and two genes were strongly expressed in the treatment group. According to the ROC curves, PDHB and GLS demonstrated significant diagnostic value. Additionally, correlation analysis was conducted on the nine characteristic genes in relation to immune infiltration. The distribution of key genes in immune cells was determined using single-cell RNA sequencing data. Finally, the mRNA expression of the nine diagnostic genes was validated using qPCR.

Integrin subunit beta 6 is a potential diagnostic marker for acute kidney injury in patients with diabetic kidney disease: a single cell sequencing data analysis

Background Diabetic kidney disease (DKD), a prevalent complication of diabetes mellitus, is often associated with acute kidney injury (AKI). Thus, the development of preventive and therapeutic strategies is crucial for delaying the progression of AKI and DKD. Methods The GSE183276 dataset, comprising the data of 20 healthy controls and 12 patients with AKI, was downloaded from the Gene Expression Omnibus (GEO) database to analyze the AKI group. For analyzing the DKD group, the GSE131822 dataset, comprising the data of 3 healthy controls and 3 patients with DKD, was downloaded from the GEO database. The common differentially expressed genes (DEGs) in renal tubular epithelial cells (TECs) were subjected to enrichment analyses. Next, a protein-protein interaction (PPI) network was constructed using the Search Tool for the Retrieval of Interacting Genes database to analyze gene-related regulatory networks. Finally, the AKI animal models and the DKD and AKI cell models were established, and the reliability of the identified genes was validated using quantitative real-time polymerase chain reaction analysis. Results Functional analysis was performed with 40 common DEGs in TECs. Eight hub genes were identified using the PPI and gene-related networks. Finally, validation experiments with the in vivo animal model and the in vitro cellular model revealed the four common DEGs. Four DEGs that share molecular mechanisms in the pathogenesis of DKD and AKI were identified. In particular, the expression of Integrin Subunit Beta 6(ITGB6), a hub and commonly upregulated gene, was upregulated in the in vitro models. Conclusion ITGB6 may serve as a biomarker for early AKI diagnosis in patients with DKD and as a target for early intervention therapies.

A single amplified genome catalog reveals the dynamics of mobilome and resistome in the human microbiome

BackgroundThe increase in metagenome-assembled genomes (MAGs) has advanced our understanding of the functional characterization and taxonomic assignment within the human microbiome. However, MAGs, as population consensus genomes, often aggregate heterogeneity among species and strains, thereby obfuscating the precise relationships between microbial hosts and mobile genetic elements (MGEs). In contrast, single amplified genomes (SAGs) derived via single-cell genome sequencing can capture individual genomic content, including MGEs.ResultsWe introduce the first substantial SAG dataset (bbsag20) from the human oral and gut microbiome, comprising 17,202 SAGs above medium-quality without co-assembly. This collection unveils a diversity of bacterial lineages across 312 oral and 647 gut species, demonstrating different taxonomic compositions from MAGs. Moreover, the SAGs showed cellular-level evidence of the translocation of oral bacteria to the gut. We also identified broad-host-range MGEs harboring antibiotic resistance genes (ARGs), which were not detected in the MAGs.ConclusionsThe difference in taxonomic composition between SAGs and MAGs indicates that combining both methods would be effective in expanding the genome catalog. By connecting mobilomes and resistomes in individual samples, SAGs could meticulously chart a dynamic network of ARGs on MGEs, pinpointing potential ARG reservoirs and their spreading patterns in the microbial community.DX-3H-wDoVG35Spx9tFU_UVideo

Advances in the application of single-cell sequencing technology in the study of neonatal necrotizing enterocolitis

Advances in the application of single-cell sequencing technology in the study of neonatal necrotizing enterocolitis

The critical involvement of monocytes/macrophages in the pathogenesis of primary Sjögren's syndrome: New evidence from Mendelian randomization and single-cell sequencing

BackgroundPrimary Sjögren's syndrome (pSS) stands as a chronic autoimmune disease characterized by an elusive pathogenesis. The synergy of single-cell RNA sequencing and Mendelian randomization (MR) analysis provides an opportunity to comprehensively unravel the contributory role of monocytes/macrophages in the intricate pathogenesis of pSS. MethodsDifferentially expressed genes (DEGs) of various types of immune cells were analyzed after annotating single-cell RNA sequencing (scRNA-seq) data. MR analysis of expression quantitative trait loci (eQTL) and protein quantitative trait loci (pQTL) was conducted to search for key pathogenic genes and proteins. Cellular localization of pathogenic genes was performed based on scRNA-seq data. Variations in signaling pathways between immune cells were further analyzed. ResultsA total of 1434 significant DEGs were identified. Among these, 60 genes exhibited strong relevance to the occurrence of pSS, of which 32 genes differentially expressed in monocytes/macrophages. CTSS was found to be a significant risk protein with a p-value of 0.001 and an odds ratio of 1.384 (1.147–1.669), showing pronounced expression in monocytes/macrophages. Furthermore, monocytes/macrophages displayed heightened expression levels of MXD1, AMPD2, TNFSF10, FTL, UBXN11, CSF3R, and LILRA5. The analysis of intercellular signaling revealed increased signal intensity in both incoming and outgoing signals in monocytes/macrophages. The signaling interactions between monocytes/macrophages, B cells, and T cells exhibited varying degrees of deviation. ConclusionsThis study highlights the significant involvement of monocytes/macrophages in the pathogenesis of pSS, as evidenced by MR analysis and scRNA-seq analysis. This suggests monocytes/macrophages as a focal point for pathogenesis research and potential therapeutic targeting in pSS.

Single-cell sequencing decodes the secrets of the RAP phenomenon of corticotomy.

Corticotomy-assisted tooth movement is commonly performed in clinics, however, its time-limited efficacy and the fear of surgery among patients significantly limit its clinical application. Hence, researchers have investigated non-invasive methods to accelerate tooth movement. However, the molecular mechanisms underlying corticotomy-assisted tooth movement are not fully understood. Micro-CT and TRAP stain were used to tooth movement and bone resorption. Single-cell RNA sequencing was used to study the transcriptome heterogeneity of macrophages after corticotomy. Transmission electron microscopy and iron ion detection was used to evaluate ferroptosis and iron metabolism. In addition, we carried out immunohistochemistry, quantitative real-time and flow cytometry verify the effect of iron on macrophage polarization. Single-cell RNA sequencing of digested alveolar bone identified a significant increase in iron metabolism-related genes post-corticotomy. Macrophages play a central role in this field. Following the dimensionality reduction of macrophages, we revealed a new developmental state via pseudotime analysis post-corticotomy. SCENIC analysis revealed that Atf3 is a key transcription factor influencing this new state. We found that Atf3+ macrophages were closely associated with osteoclasts. Moreover, cell chat revealed an increase in cellular communication between Atf3+ macrophages and other cell types after corticotomy. These findings suggested that Atf3+ macrophages might play a key role in corticotomy-accelerated tooth movement, thus providing potential targets for drug development.

Processing pipelines and analytical methods for single-cell DNA methylation sequencing data.

Single-cell DNA methylation sequencing technology has seen rapid advancements in recent years, playing a crucial role in uncovering cellular heterogeneity and the mechanisms of epigenetic regulation. As sequencing technologies have progressed, the quality and quantity of single-cell methylation data have also increased, making standardized preprocessing workflows and appropriate analysis methods essential for ensuring data comparability and result reliability. However, a comprehensive data analysis pipeline to guide researchers in mining existing data has yet to be established. This review systematically summarizes the preprocessing steps and analysis methods for single-cell methylation data, introduces relevant algorithms and tools, and explores the application prospects of single-cell methylation technology in neuroscience, hematopoietic differentiation, and cancer research. The aim is to provide guidance for researchers in data analysis and to promote the development and application of single-cell methylation sequencing technology.

Osteoarthritis year in review 2024: Biology

Osteoarthritis (OA) research is a fast-growing and extremely wide field, in which a substantial increase in knowledge has been achieved over the last year. It covers many different topics, however, a PubMed search using the terms 'osteoarthritis' and 'biology' resulted in only a limited number of studies that were published between April 2023 and April 2024. In order to identify OA-relevant studies that focus on mechanistic studies of biological processes at the tissue, cellular, and molecular level, the following keywords were included as search terms: tissue interactions, single cell sequencing, transcriptomics, extracellular matrix, signaling, ion channels, and pain. The final selection of publications presented in this 'year in review' was influenced by the personal preferences of the authors, and eventually three larger key themes emerged: 1) Joint tissue interactions covering meniscus, subchondral bone, fat tissue, synovium, and synovial fluid. 2) Degeneration of the cartilage extracellular matrix and generation of bioactive fragments. 3) Receptors, ion channels, signaling pathways, and cellular metabolism. Many of the studies summarized here identified novel potential targets for OA treatment, and promising results were already obtained addressing these targets in different animal models. It will be exciting to see which findings can be translated into future clinical studies and eventually lead to novel treatment approaches for human OA.

Inferring replication timing and proliferation dynamics from single-cell DNA sequencing data

Dysregulated DNA replication is a cause and a consequence of aneuploidy in cancer, yet the interplay between copy number alterations (CNAs), replication timing (RT) and cell cycle dynamics remain understudied in aneuploid tumors. We developed a probabilistic method, PERT, for simultaneous inference of cell-specific replication and copy number states from single-cell whole genome sequencing (scWGS) data. We used PERT to investigate clone-specific RT and proliferation dynamics in >50,000 cells obtained from aneuploid and clonally heterogeneous cell lines, xenografts and primary cancers. We observed bidirectional relationships between RT and CNAs, with CNAs affecting X-inactivation producing the largest RT shifts. Additionally, we found that clone-specific S-phase enrichment positively correlated with ground-truth proliferation rates in genomically stable but not unstable cells. Together, these results demonstrate robust computational identification of S-phase cells from scWGS data, and highlight the importance of RT and cell cycle properties in studying the genomic evolution of aneuploid tumors.

Spatial Transcriptomics of the Respiratory System.

Over the last decade, single-cell genomics has revealed remarkable heterogeneity and plasticity of cell types in the lungs and airways. The challenge now is to understand how these cell types interact in three-dimensional space to perform lung functions, facilitating airflow and gas exchange while simultaneously providing barrier function to avoid infection. An explosion in novel spatially resolved gene expression technologies, coupled with computational tools that harness machine learning and deep learning, now promise to address this challenge. Here, we review the most commonly used spatial analysis workflows, highlighting their advantages and limitations, and outline recent developments in machine learning and artificial intelligence that will augment how we interpret spatial data. Together these technologies have the potential to transform our understanding of the respiratory system in health and disease, and we showcase studies in lung development, COVID-19, lung cancer, and fibrosis where spatially resolved transcriptomics is already providing novel insights.

Energy Score Analysis in Immune Cells Linked to Atypical Hemolytic Uremic Syndrome Using Single-Cell Sequencing

Energy Score Analysis in Immune Cells Linked to Atypical Hemolytic Uremic Syndrome Using Single-Cell Sequencing

Research progress on single-cell expression quantitative trait loci.

Expression quantitative trait loci (eQTL) represent genetic variants that regulate gene expression levels. eQTL analysis has become a crucial method for identifying the functional roles of disease-associated genetic variants in the post-genome-wide association study (GWAS) era, yielding numerous significant discoveries. Traditional eQTL analysis relies on whole-genome sequencing combined with bulk RNA-seq, which obscures gene expression differences between cells and thus fails to identify cell type- or state-dependent eQTL. This limitation makes it challenging to elucidate the roles of disease-associated genetic variants under specific conditions. In recent years, with the development and widespread application of single-cell RNA sequencing (scRNA-seq) technology, scRNA-seq-based eQTL (sc-eQTL) research has emerged as a focal point. The advantage of this approach lies in its ability to leverage the resolution and granularity of single-cell sequencing to uncover eQTL that are dependent on cell type, cell state, and cellular dynamics. This significantly enhances our ability to analyze genetic variants associated with gene expression. Consequently, it holds substantial significance for advancing our understanding of the formation of complex organs and the mechanisms underlying disease onset, progression, intervention, and treatment. This review comprehensively examines the recent advancements in sc-eQTL studies, focusing on their development, experimental design strategies, modeling approaches, and current challenges. The aim is to offer researchers novel perspectives for identifying disease-associated loci and elucidating gene regulatory mechanisms.

Heterogeneity in extracellular matrix and immune microenvironment of anterior vaginal wall revealed by single-cell sequencing in women with stress urinary incontinence

Stress urinary incontinence (SUI), characterized by involuntary urine leakage during increased abdominal pressure, remains poorly understood regarding its pathophysiology and treatment. In this study, we utilized single-cell sequencing to analyze the transcriptomic profiles of different cell types in anterior vaginal wall of SUI patients, aiming to explore the heterogeneity of the extracellular matrix (ECM) and immune microenvironment in SUI pathogenesis. Our results identified eleven cell types, including connective tissue cells, immune cells, and glial cells. Specifically, fibroblasts, smooth muscle cells, epithelial cells and T cells displayed transcriptional characteristics highly relevant to SUI pathogenesis. We observed that most cell types participate in ECM metabolism and immune-inflammatory responses, indicating a synergistic role of multiple vaginal cell types in SUI. Furthermore, altered intercellular communication, particularly between fibroblasts and T cells, was noted in SUI. This study provides novel single-cell insights into SUI and identifies potential biomarkers and therapeutic targets for future research.

Dissecting the Single-Cell Diversity and Heterogeneity Underlying Cervical Precancerous Lesions and Cancer Tissues.

The underlying cellular diversity and heterogeneity from cervix precancerous lesions to cervical squamous cell carcinoma (CSCC) is investigated. Four single-cell datasets including normal tissues, normal adjacent tissues, precancerous lesions, and cervical tumors were integrated to perform disease stage analysis. Single-cell compositional data analysis (scCODA) was utilized to reveal the compositional changes of each cell type. Differentially expressed genes (DEGs) among cell types were annotated using BioCarta. An assay for transposase-accessible chromatin sequencing (ATAC-seq) analysis was performed to correlate epigenetic alterations with gene expression profiles. Lastly, a logistic regression model was used to assess the similarity between the original and new cohort data (HRA001742). After global annotation, seven distinct cell types were categorized. Eight consensus-upregulated DEGs were identified in B cells among different disease statuses, which could be utilized to predict the overall survival of CSCC patients. Inferred copy number variation (CNV) analysis of epithelial cells guided disease progression classification. Trajectory and ATAC-seq integration analysis identified 95 key transcription factors (TF) and one immunohistochemistry (IHC) testified key-node TF (YY1) involved in epithelial cells from CSCC initiation to progression. The consistency of epithelial cell subpopulation markers was revealed with single-cell sequencing, bulk sequencing, and RT-qPCR detection. KRT8 and KRT15, markers of Epi6, showed progressively higher expression with disease progression as revealed by IHC detection. The logistic regression model testified the robustness of the resemblance of clusters among the various datasets utilized in this study. Valuable insights into CSCC cellular diversity and heterogeneity provide a foundation for future targeted therapy.