Unravelling biology and shifting paradigms in cancer with single-cell sequencing

Single-cell multi-omics in the study of digestive system cancers

Digestive system cancers are prevalent diseases with a high mortality rate, posing a significant threat to public health and economic burden. The diagnosis and treatment of digestive system cancer confront conventional cancer problems, such as tumor heterogeneity and drug resistance. Single-cell sequencing (SCS) emerged at times required and has developed from single-cell RNA-seq (scRNA-seq) to the single-cell multi-omics era represented by single-cell spatial transcriptomics (ST). This article comprehensively reviews the advances of single-cell omics technology in the study of digestive system tumors. While analyzing and summarizing the research cases, vital details on the sequencing platform, sample information, sampling method, and key findings are provided. Meanwhile, we summarize the commonly used SCS platforms and their features, as well as the advantages of multi-omics technologies in combination. Finally, the development trends and prospects of the application of single-cell multi-omics technology in digestive system cancer research are prospected.

Single-cell mass cytometry in immunological skin diseases

Immune-related skin diseases represent a collective of dermatological disorders intricately linked to dysfunctional immune system processes. These conditions are primarily characterized by an immoderate activation of the immune system or deviant immune responses, involving diverse immune components including immune cells, antibodies, and inflammatory mediators. However, the precise molecular dysregulation underlying numerous individual cases of these diseases and unique subsets respond under disease conditions remains elusive. Comprehending the mechanisms and determinants governing the homeostasis and functionality of diseases could offer potential therapeutic opportunities for intervention. Mass cytometry enables precise and high-throughput quantitative measurement of proteins within individual cells by utilizing antibodies labeled with rare heavy metal isotopes. Imaging mass cytometry employs mass spectrometry to obtain spatial information on cell-to-cell interactions within tissue sections, simultaneously utilizing more than 40 markers. The application of single-cell mass cytometry presents a unique opportunity to conduct highly multiplexed analysis at the single-cell level, thereby revolutionizing our understanding of cell population heterogeneity and hierarchy, cellular states, multiplexed signaling pathways, proteolysis products, and mRNA transcripts specifically in the context of many autoimmune diseases. This information holds the potential to offer novel approaches for the diagnosis, prognostic assessment, and monitoring responses to treatment, thereby enriching our strategies in managing the respective conditions. This review summarizes the present-day utilization of single-cell mass cytometry in studying immune-related skin diseases, highlighting its advantages and limitations. This technique will become increasingly prevalent in conducting extensive investigations into these disorders, ultimately yielding significant contributions to their accurate diagnosis and efficacious therapeutic interventions.

Advancements in 3D In Vitro Models for Colorectal Cancer

The process of drug discovery and pre-clinical testing is currently inefficient, expensive, and time-consuming. Most importantly, the success rate is unsatisfactory, as only a small percentage of tested drugs are made available to oncological patients. This is largely due to the lack of reliable models that accurately predict drug efficacy and safety. Even animal models often fail to replicate human-specific pathologies and human body's complexity. These factors, along with ethical concerns regarding animal use, urge the development of suitable human-relevant, translational in vitro models.

Combining single-cell and bulk RNA sequencing, NK cell marker genes reveal a prognostic and immune status in pancreatic ductal adenocarcinoma

The NK cell is an important component of the tumor microenvironment of pancreatic ductal adenocarcinoma (PDAC), also plays a significant role in PDAC development. This study aimed to explore the relationship between NK cell marker genes and prognosis, immune response of PDAC patients. By scRNA-seq data, we found the proportion of NK cells were significantly downregulated in PDAC and 373 NK cell marker genes were screened out. By TCGA database, we enrolled 7 NK cell marker genes to construct the signature for predicting prognosis in PDAC patients. Cox analysis identified the signature as an independent factor for pancreatic cancer. Subsequently, the predictive power of signature was validated by 6 GEO datasets and had an excellent evaluation. Our analysis of relationship between the signature and patients' immune status revealed that the signature has a strong correlation with immunocyte infiltration, inflammatory reaction, immune checkpoint inhibitors (ICIs) response. The NK cell marker genes are closely related to the prognosis and immune capacity of PDAC patients, and they have potential value as a therapeutic target.

Single-Cell RNA Sequencing: Technological Progress and Biomedical Application in Cancer Research

Single-cell RNA-seq (scRNA-seq) is a revolutionary technology that allows for the genomic investigation of individual cells in a population, allowing for the discovery of unusual cells associated with cancer and metastasis. ScRNA-seq has been used to discover different types of cancers with poor prognosis and medication resistance such as lung cancer, breast cancer, ovarian cancer, and gastric cancer. Besides, scRNA-seq is a promising method that helps us comprehend the biological features and dynamics of cell development, as well as other disorders. This review gives a concise summary of current scRNA-seq technology. We also explain the main technological steps involved in implementing the technology. We highlight the present applications of scRNA-seq in cancer research, including tumor heterogeneity analysis in lung cancer, breast cancer, and ovarian cancer. In addition, this review elucidates potential applications of scRNA-seq in lineage tracing, personalized medicine, illness prediction, and disease diagnosis, which reveals that scRNA-seq facilitates these events by producing genetic variations on the single-cell level.

Single-cell analysis in rheumatic and allergic diseases: insights for clinical practice

Since the advent of single-cell RNA sequencing (scRNA-seq) methodology, single-cell analysis has become a powerful tool for exploration of cellular networks and dysregulated immune responses in disease pathogenesis. Advanced bioinformatics tools have enabled the combined analysis of scRNA-seq data and information on various cell properties, such as cell surface molecular profiles, chromatin accessibility and spatial information, leading to a deeper understanding of pathology. This Review provides an overview of the achievements in single-cell analysis applied to clinical samples of rheumatic and allergic diseases, including rheumatoid arthritis, systemic lupus erythematosus, systemic sclerosis, allergic airway diseases and atopic dermatitis, with an expanded scope beyond peripheral blood cells to include local diseased tissues. Despite the valuable insights that single-cell analysis has provided into disease pathogenesis, challenges remain in translating single-cell findings into clinical practice and developing personalized treatment strategies. Beyond understanding the atlas of cellular diversity, we discuss the application of data obtained in each study to clinical practice, with a focus on identifying biomarkers and therapeutic targets.

ScRNA-seq revealed the tumor microenvironment heterogeneity related to the occurrence and metastasis in upper urinary tract urothelial carcinoma

Metastasis is the greatest clinical challenge for UTUCs, which may have distinct molecular and cellular characteristics from earlier cancers. Herein, we provide single-cell transcriptome profiles of UTUC para cancer normal tissue, primary tumor lesions, and lymphatic metastases to explore possible mechanisms associated with UTUC occurrence and metastasis. From 28,315 cells obtained from normal and tumor tissues of 3 high-grade UTUC patients, we revealed the origin of UTUC tumor cells and the homology between metastatic and primary tumor cells. Unlike the immunomicroenvironment suppression of other tumors, we found no immunosuppression in the tumor microenvironment of UTUC. Moreover, it is imperative to note that stromal cells are pivotal in the advancement of UTUC. This comprehensive single-cell exploration enhances our comprehension of the molecular and cellular dynamics of metastatic UTUCs and discloses promising diagnostic and therapeutic targets in cancer-microenvironment interactions.

Progress of single-cell RNA sequencing combined with spatial transcriptomics in tumour microenvironment and treatment of pancreatic cancer

In recent years, single-cell analyses have revealed the heterogeneity of the tumour microenvironment (TME) at the genomic, transcriptomic, and proteomic levels, further improving our understanding of the mechanisms of tumour development. Single-cell RNA sequencing (scRNA-seq) technology allow analysis of the transcriptome at the single-cell level and have unprecedented potential for exploration of the characteristics involved in tumour development and progression. These techniques allow analysis of transcript sequences at higher resolution, thereby increasing our understanding of the diversity of cells found in the tumour microenvironment and how these cells interact in complex tumour tissue. Although scRNA-seq has emerged as an important tool for studying the tumour microenvironment in recent years, it cannot be used to analyse spatial information for cells. In this regard, spatial transcriptomics (ST) approaches allow researchers to understand the functions of individual cells in complex multicellular organisms by understanding their physical location in tissue sections. In particular, in related research on tumour heterogeneity, ST is an excellent complementary approach to scRNA-seq, constituting a new method for further exploration of tumour heterogeneity, and this approach can also provide unprecedented insight into the development of treatments for pancreatic cancer (PC). In this review, based on the methods of scRNA-seq and ST analyses, research progress on the tumour microenvironment and treatment of pancreatic cancer is further explained.

Integration of single-cell sequencing and bulk RNA-seq to identify and develop a prognostic signature related to colorectal cancer stem cells

The prognosis for patients with colorectal cancer (CRC) remains worse than expected due to metastasis, recurrence, and resistance to chemotherapy. Colorectal cancer stem cells (CRCSCs) play a vital role in tumor metastasis, recurrence, and chemotherapy resistance. However, there are currently no prognostic markers based on CRCSCs-related genes available for clinical use. In this study, single-cell transcriptome sequencing was employed to distinguish cancer stem cells (CSCs) in the CRC microenvironment and analyze their properties at the single-cell level. Subsequently, data from TCGA and GEO databases were utilized to develop a prognostic risk model for CRCSCs-related genes and validate its diagnostic performance. Additionally, functional enrichment, immune response, and chemotherapeutic drug sensitivity of the relevant genes in the risk model were investigated. Lastly, the key gene RPS17 in the risk model was identified as a potential prognostic marker and therapeutic target for further comprehensive studies. Our findings provide new insights into the prognostic treatment of CRC and offer novel perspectives for a systematic and comprehensive understanding of CRC development.

Raman cell sorting for single-cell research

Cells constitute the fundamental units of living organisms. Investigating individual differences at the single-cell level facilitates an understanding of cell differentiation, development, gene expression, and cellular characteristics, unveiling the underlying laws governing life activities in depth. In recent years, the integration of single-cell manipulation and recognition technologies into detection and sorting systems has emerged as a powerful tool for advancing single-cell research. Raman cell sorting technology has garnered attention owing to its non-labeling, non-destructive detection features and the capability to analyze samples containing water. In addition, this technology can provide live cells for subsequent genomics analysis and gene sequencing. This paper emphasizes the importance of single-cell research, describes the single-cell research methods that currently exist, including single-cell manipulation and single-cell identification techniques, and highlights the advantages of Raman spectroscopy in the field of single-cell analysis by comparing it with the fluorescence-activated cell sorting (FACS) technique. It describes various existing Raman cell sorting techniques and introduces their respective advantages and disadvantages. The above techniques were compared and analyzed, considering a variety of factors. The current bottlenecks include weak single-cell spontaneous Raman signals and the requirement for a prolonged total cell exposure time, significantly constraining Raman cell sorting technology's detection speed, efficiency, and throughput. This paper provides an overview of current methods for enhancing weak spontaneous Raman signals and their associated advantages and disadvantages. Finally, the paper outlines the detailed information related to the Raman cell sorting technology mentioned in this paper and discusses the development trends and direction of Raman cell sorting.

Spatially resolved analysis of microenvironmental gradient impact on cancer cell phenotypes

Despite the physiological and pathophysiological significance of microenvironmental gradients, e.g., for diseases such as cancer, tools for generating such gradients and analyzing their impact are lacking. Here, we present an integrated microfluidic-based workflow that mimics extracellular pH gradients characteristic of solid tumors while enabling high-resolution live imaging of, e.g., cell motility and chemotaxis, and preserving the capacity to capture the spatial transcriptome. Our microfluidic device generates a pH gradient that can be rapidly controlled to mimic spatiotemporal microenvironmental changes over cancer cells embedded in a 3D matrix. The device can be reopened allowing immunofluorescence analysis of selected phenotypes, as well as the transfer of cells and matrix to a Visium slide for spatially resolved analysis of transcriptional changes across the pH gradient. This workflow is easily adaptable to other gradients and multiple cell types and can therefore prove invaluable for integrated analysis of roles of microenvironmental gradients in biology.

Integrating single-cell and spatial transcriptomics to elucidate the crosstalk between cancer-associated fibroblasts and cancer cells in hepatocellular carcinoma with spleen-deficiency syndrome

Background and aim

Most patients with hepatocellular carcinoma (HCC) in China have been diagnosed with spleen deficiency syndrome (SDS), which accelerates the progression of HCC by disrupting the tumor microenvironment homeostasis. This study aimed to investigate the intercellular crosstalk in HCC with SDS.

Experimental procedure

An HCC-SDS mouse model was established using orthotopic HCC transplantation based on reserpine-induced SDS. Single-cell data analysis and cancer cell prediction were conducted using Seurat and CopyKAT package, respectively. Intercellular interactions were explored using CellPhoneDB and CellChat and subsequently validated using co-culture assays, ELISA and histological staining. We performed pathway activity analysis using gene set variation analysis and the Seurat package. The extracellular matrix (ECM) remodeling was assessed using a gel contraction assay, atomic force microscopy, and Sirius red staining. The deconvolution of the spatial transcriptomics data using the "CARD" package based on single-cell data.

Results and conclusion

We successfully established the HCC-SDS mouse model. Twenty-nine clusters were identified. The interactions between cancer cells and cancer-associated fibroblasts (CAFs) were significantly enhanced via platelet-derived growth factor (PDGF) signaling in HCC-SDS. CAFs recruited in HCC-SDS lead to ECM remodeling and the activation of TGF-β signaling pathway. Deconvolution of the spatial transcriptome data revealed that CAFs physically surround cancer cells in HCC-SDS. This study reveals that the crosstalk of CAFs-cancer cells is crucial for the tumor-promoting effect of SDS. CAFs recruited by HCC via PDGFA may lead to ECM remodeling through activation of the TGF-β pathway, thereby forming a physical barrier to block immune cell infiltration under SDS.

Integrating single-cell and spatial analysis reveals MUC1-mediated cellular crosstalk in mucinous colorectal adenocarcinoma

BACKGROUND

Mucinous colorectal adenocarcinoma (MCA) is a distinct subtype of colorectal cancer (CRC) with the most aggressive pattern, but effective treatment of MCA remains a challenge due to its vague pathological characteristics. An in-depth understanding of transcriptional dynamics at the cellular level is critical for developing specialised MCA treatment strategies.

METHODS

We integrated single-cell RNA sequencing and spatial transcriptomics data to systematically profile the MCA tumor microenvironment (TME), particularly the interactome of stromal and immune cells. In addition, a three-dimensional bioprinting technique, canonical ex vivo co-culture system, and immunofluorescence staining were further applied to validate the cellular communication networks within the TME.

RESULTS

This study identified the crucial intercellular interactions that engaged in MCA pathogenesis. We found the increased infiltration of FGF7+/THBS1+ myofibroblasts in MCA tissues with decreased expression of genes associated with leukocyte-mediated immunity and T cell activation, suggesting a crucial role of these cells in regulating the immunosuppressive TME. In addition, MS4A4A+ macrophages that exhibit M2-phenotype were enriched in the tumoral niche and high expression of MS4A4A+ was associated with poor prognosis in the cohort data. The ligand-receptor-based intercellular communication analysis revealed the tight interaction of MUC1+ malignant cells and ZEB1+ endothelial cells, providing mechanistic information for MCA angiogenesis and molecular targets for subsequent translational applications.

CONCLUSIONS

Our study provides novel insights into communications among tumour cells with stromal and immune cells that are significantly enriched in the TME during MCA progression, presenting potential prognostic biomarkers and therapeutic strategies for MCA.

KEY POINTS

Tumour microenvironment profiling of MCA is developed. MUC1+ tumour cells interplay with FGF7+/THBS1+ myofibroblasts to promote MCA development. MS4A4A+ macrophages exhibit M2 phenotype in MCA. ZEB1+ endotheliocytes engage in EndMT process in MCA.

Identification and validation of a novel predictive signature based on hepatocyte-specific genes in hepatocellular carcinoma by integrated analysis of single-cell and bulk RNA sequencing

BACKGROUND

Hepatocellular carcinoma represents a significant global burden in terms of cancer-related mortality, posing a substantial risk to human health. Despite the availability of various treatment modalities, the overall survival rates for patients with hepatocellular carcinoma remain suboptimal. The objective of this study was to explore the potential of novel biomarkers and to establish a novel predictive signature utilizing multiple transcriptome profiles.

METHODS

The GSE115469 and CNP0000650 cohorts were utilized for single cell analysis and gene identification. The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) datasets were utilized in the development and evaluation of a predictive signature. The expressions of hepatocyte-specific genes were further validated using the GSE135631 cohort. Furthermore, immune infiltration results, immunotherapy response prediction, somatic mutation frequency, tumor mutation burden, and anticancer drug sensitivity were analyzed based on various risk scores. Subsequently, functional enrichment analysis was performed on the differential genes identified in the risk model. Moreover, we investigated the expression of particular genes in chronic liver diseases utilizing datasets GSE135251 and GSE142530.

RESULTS

Our findings revealed hepatocyte-specific genes (ADH4, LCAT) with notable alterations during cell maturation and differentiation, leading to the development of a novel predictive signature. The analysis demonstrated the efficacy of the model in predicting outcomes, as evidenced by higher risk scores and poorer prognoses in the high-risk group. Additionally, a nomogram was devised to forecast the survival rates of patients at 1, 3, and 5 years. Our study demonstrated that the predictive model may play a role in modulating the immune microenvironment and impacting the anti-tumor immune response in hepatocellular carcinoma. The high-risk group exhibited a higher frequency of mutations and was more likely to benefit from immunotherapy as a treatment option. Additionally, we confirmed that the downregulation of hepatocyte-specific genes may indicate the progression of hepatocellular carcinoma and aid in the early diagnosis of the disease.

CONCLUSION

Our research findings indicate that ADH4 and LCAT are genes that undergo significant changes during the differentiation of hepatocytes into cancer cells. Additionally, we have created a unique predictive signature based on genes specific to hepatocytes.

Methylated lncRNAs suppress apoptosis of gastric cancer stem cells via the lncRNA-miRNA/protein axis

BACKGROUND

Long noncoding RNAs (lncRNAs) play essential roles in the tumorigenesis of gastric cancer. However, the influence of lncRNA methylation on gastric cancer stem cells (GCSCs) remains unclear.

METHODS

The N6-methyladenosine (m6A) levels of lncRNAs in gastric cancer stem cells were detected by methylated RNA immunoprecipitation sequencing (MeRIP-seq), and the results were validated by MeRIP-quantitative polymerase chain reaction (qPCR). Specific sites of m6A modification on lncRNAs were detected by single-base elongation- and ligation-based qPCR amplification (SELECT). By constructing and transfecting the plasmid expressing methyltransferase-like 3 (METTL3) fused with catalytically inactivated Cas13 (dCas13b) and guide RNA targeting specific methylation sites of lncRNAs, we obtained gastric cancer stem cells with site-specific methylation of lncRNAs. Reverse transcription (RT)-qPCR and Western blot were used for detecting the stemness of treated gastric cancer stem cells.

RESULTS

The site-specific methylation of PSMA3-AS1 and MIR22HG suppressed apoptosis and promoted stemness of GCSCs. LncRNA methylation enhanced the stability of PSMA3-AS1 and MIR22HG to suppress apoptosis of GCSCs via the PSMA3-AS1-miR-411-3p- or MIR22HG-miR-24-3p-SERTAD1 axis. Simultaneously, the methylated lncRNAs promoted the interaction between PSMA3-AS1 and the EEF1A1 protein or MIR22HG and the LRPPRC protein, stabilizing the proteins and leading to the suppression of apoptosis. The in vivo data revealed that the methylated PSMA3-AS1 and MIR22HG triggered tumorigenesis of GCSCs.

CONCLUSIONS

Our study revealed the requirement for site-specific methylation of lncRNAs in the tumorigenesis of GCSCs, contributing novel insights into cancer development.

Advances in Engineered Macrophages: A New Frontier in Cancer Immunotherapy

Macrophages, as pivotal cells within the tumour microenvironment, significantly influence the impact of and reactions to treatments for solid tumours. The rapid evolution of bioengineering technology has revealed the vast potential of engineered macrophages in immunotherapy, disease diagnosis, and tissue engineering. Given this landscape, the goal of harnessing and innovating macrophages as a novel strategy for solid tumour immunotherapy cannot be overstated. The diverse strategies for engineered macrophages in the realm of cancer immunotherapy, encompassing macrophage drug delivery systems, chimeric antigen receptor macrophage therapy, and synergistic treatment approaches involving bacterial outer membrane vesicles and macrophages, are meticulously examined in this review. These methodologies are designed to enhance the therapeutic efficacy of macrophages against solid tumours, particularly those that are drug-resistant and metastatic. Collectively, these immunotherapies are poised to supplement and refine current solid tumour treatment paradigms, thus heralding a new frontier in the fight against malignant tumours.

Identification of Treg-related prognostic molecular subtypes and individualized characteristics in clear cell renal cell carcinoma through single-cell transcriptomes and bulk RNA sequencing

BACKGROUND

In clear cell renal cell carcinoma (ccRCC), the role of Regulatory T cells (Treg cells) as prognostic and immunotherapy response predictors is not fully explored.

METHODS

Analyzing renal clear cell carcinoma datasets from TISCH, TCGA, and GEO, we focused on 8 prognostic Treg genes to study patient subtypes in ccRCC. We assessed Treg subtypes in relation to patient prognosis, tumor microenvironment, metabolism. Using Cox regression and principal component analysis, we devised Treg scores for individual patient characterization and explored the molecular role of C1QL1, a critical gene in the Treg model, through in vivo and in vitro studies.

RESULTS

Eight Treg-associated prognostic genes were identified, classifying ccRCC patients into cluster A and B. Cluster A patients showed poorer prognosis with distinct clinical and molecular profiles, potentially benefiting more from immunotherapy. Low Treg scores correlated with worse outcomes and clinical progression. Low scores also suggested that patients might respond better to immunotherapy and targeted therapies. In ccRCC, C1QL1 knockdown reduced tumor proliferation and invasion via NF-kb-EMT pathways and decreased Treg cell infiltration, enhancing immune efficacy.

CONCLUSIONS

The molecular subtype and Treg score in ccRCC, based on Treg cell marker genes, are crucial in personalizing ccRCC treatment and underscore C1QL1's potential as a tumor biomarker and target for immunotherapy.

Integrative analysis of bulk and single-cell RNA-seq reveals the molecular characterization of the immune microenvironment and oxidative stress signature in melanoma

Background

The immune microenvironment and oxidative stress of melanoma show significant heterogeneity, which affects tumor growth, invasion and treatment response. Single-cell and bulk RNA-seq data were used to explore the heterogeneity of the immune microenvironment and oxidative stress of melanoma.

Methods

The R package Seurat facilitated the analysis of the single-cell dataset, while Harmony, another R package, was employed for batch effect correction. Cell types were classified using Uniform Manifold Approximation and Projection (UMAP). The Secreted Signaling algorithm from CellChatDB.human was applied to elucidate cell-to-cell communication patterns within the single-cell data. Consensus clustering analysis for the skin cutaneous melanoma (SKCM) samples was executed with the R package ConsensusClusterPlus. To quantify immune infiltrating cells, we utilized CIBERSORT, ESTIMATE, and TIMERxCell algorithms provided by the R package Immuno-Oncology Biological Research (IOBR). Single nucleotide variant (SNV) analysis was conducted using Maftools, an R package specifically designed for this purpose. Subsequently, the expression levels of PXDN and PAPSS2 genes were assessed in melanoma tissues compared to adjacent normal tissues. Furthermore, in vitro experiments were conducted to evaluate the proliferation and reactive oxygen species expression in melanoma cells following transfection with siRNA targeting PXDN and PAPSS2.

Results

Malignant tumor cell populations were reclassified based on a comprehensive single-cell dataset analysis, which yielded six distinct tumor subsets. The specific marker genes identified for these subgroups were then used to interrogate the Cancer Genome Atlas Skin Cutaneous Melanoma (TCGA-SKCM) cohort, derived from bulk RNA sequencing data, resulting in the delineation of two immune molecular subtypes. Notably, patients within the cluster2 (C2) subtype exhibited a significantly more favorable prognosis compared to those in the cluster1 (C1) subtype. An alignment of immune characteristics was observed between the C2 subtype and unique immune functional tumor cell subsets. Genes differentially expressed across these subtypes were subsequently leveraged to construct a predictive risk model. In vitro investigations further revealed elevated expression levels of PXDN and PAPSS2 in melanoma tissue samples. Functional assays indicated that modulation of PXDN and PAPSS2 expression could influence the production of reactive oxygen species (ROS) and the proliferative capacity of melanoma cells.

Conclusion

The constructed six-gene signature can be used as an immune response and an oxidative stress marker to guide the clinical diagnosis and treatment of melanoma.

Interpreting single-cell messages in normal and aberrant hematopoiesis with the Cell Marker Accordion

Single-cell technologies offer a unique opportunity to explore cellular heterogeneity in hematopoiesis, reveal malignant hematopoietic cells with clinically significant features and measure gene signatures linked to pathological pathways. However, reliable identification of cell types is a crucial bottleneck in single-cell analysis. Available databases contain dissimilar nomenclature and non-concurrent marker sets, leading to inconsistent annotations and poor interpretability. Furthermore, current tools focus mostly on physiological cell types, lacking extensive applicability in disease. We developed the Cell Marker Accordion, a user-friendly platform for the automatic annotation and biological interpretation of single-cell populations based on consistency weighted markers. We validated our approach on peripheral blood and bone marrow single-cell datasets, using surface markers and expert-based annotation as the ground truth. In all cases, we significantly improved the accuracy in identifying cell types with respect to any single source database. Moreover, the Cell Marker Accordion can identify disease-critical cells and pathological processes, extracting potential biomarkers in a wide variety of contexts in human and murine single-cell datasets. It characterizes leukemia stem cell subtypes, including therapy-resistant cells in acute myeloid leukemia patients; it identifies malignant plasma cells in multiple myeloma samples; it dissects cell type alterations in splicing factor-mutant cells from myelodysplastic syndrome patients; it discovers activation of innate immunity pathways in bone marrow from mice treated with METTL3 inhibitors. The breadth of these applications elevates the Cell Marker Accordion as a flexible, faithful and standardized tool to annotate and interpret hematopoietic populations in single-cell datasets focused on the study of hematopoietic development and disease.

Single-cell transcriptome sequencing provides insight into multiple chemotherapy resistance in a patient with refractory DLBCL: a case report

Relapsed and refractory diffuse large B-cell lymphoma (DLBCL) is associated with poor prognosis. As such, a comprehensive analysis of intratumoral components, intratumoral heterogeneity, and the immune microenvironment is essential to elucidate the mechanisms driving the progression of DLBCL and to develop new therapeutics. Here, we used single-cell transcriptome sequencing and conventional bulk next-generation sequencing (NGS) to understand the composite tumor landscape of a single patient who had experienced multiple tumor recurrences following several chemotherapy treatments. NGS revealed several key somatic mutations that are known to contribute to drug resistance. Based on gene expression profiles at the single-cell level, we identified four clusters of malignant B cells with distinct transcriptional signatures, showing high intra-tumoral heterogeneity. Among them, heterogeneity was reflected in activating several key pathways, human leukocyte antigen (HLA)-related molecules' expression, and key oncogenes, which may lead to multi-drug resistance. In addition, FOXP3+ regulatory CD4+ T cells and exhausted cytotoxic CD8+ T cells were identified, accounted for a significant proportion, and showed highly immunosuppressive properties. Finally, cell communication analysis indicated complex interactions between malignant B cells and T cells. In conclusion, this case report demonstrates the value of single-cell RNA sequencing for visualizing the tumor microenvironment and identifying potential therapeutic targets in a patient with treatment-refractory DLBCL. The combination of NGS and single-cell RNA sequencing may facilitate clinical decision-making and drug selection in challenging DLBCL cases.

Patterns of metastatic spread and tumor burden in unselected cancer patients using PET imaging: Implications for the oligometastatic spectrum theory

Introduction and background

Metastatic disease has been proposed as a continuum, with no clear cut-off between oligometastatic and polymetastatic disease. This study aims to quantify tumor burden and patterns of spread in unselected metastatic cancer patients referred for PET-based staging, response assessment of restaging.

Materials and methods

All oncological fluorodeoxyglucose (FDG-) and prostate-specific membrane antigen (PSMA-) positron emission tomography (PET) scans conducted at a single academic center in 2020 were analyzed. Imaging reports of all patients with metastatic disease were reviewed and assessed.

Results

For this study, 7,000 PET scans were screened. One third of PET scans (n = 1,754; 33 %) from 1,155 unique patients showed presence of metastatic disease from solid malignancies, of which 601 (52 %) and 554 (48 %) were classified as oligometastatic (maximum 5 metastases) and polymetastatic (>5 metastases), respectively. Lung and pleural cancer, skin cancer, and breast cancer were the most common primary tumor histologies with 132 (23.8 %), 88 (15.9 %), and 72 (13.0 %) cases, respectively. Analysis of the number of distant metastases showed a strong bimodal distribution of the metastatic burden with 26 % of patients having one solitary metastasis and 43 % of patients harboring >10 metastases. Yet, despite 43 % of polymetastatic patients having >10 distant metastases, their pattern of distribution was restricted to one or two organs in about two thirds of patients, and there was no association between the number of distant metastases and the number of involved organs.

Conclusion

The majority of metastatic cancer patients are characterized by either a solitary metastasis or a high tumor burden with >10 metastases, the latter was often associated with affecting a limited number of organs. These findings support both the spectrum theory of metastasis and the seed and soil hypothesis and can support in designing the next generation of clinical trials in the field of oligometastatic disease.

The role of microglia heterogeneity in synaptic plasticity and brain disorders: Will sequencing shed light on the discovery of new therapeutic targets?

Microglia play a crucial role in interacting with neuronal synapses and modulating synaptic plasticity. This function is particularly significant during postnatal development, as microglia are responsible for removing excessive synapses to prevent neurodevelopmental deficits. Dysregulation of microglial synaptic function has been well-documented in various pathological conditions, notably Alzheimer's disease and multiple sclerosis. The recent application of RNA sequencing has provided a powerful and unbiased means to decipher spatial and temporal microglial heterogeneity. By identifying microglia with varying gene expression profiles, researchers have defined multiple subgroups of microglia associated with specific pathological states, including disease-associated microglia, interferon-responsive microglia, proliferating microglia, and inflamed microglia in multiple sclerosis, among others. However, the functional roles of these distinct subgroups remain inadequately characterized. This review aims to refine our current understanding of the potential roles of heterogeneous microglia in regulating synaptic plasticity and their implications for various brain disorders, drawing from recent sequencing research and functional studies. This knowledge may aid in the identification of pathogenetic biomarkers and potential factors contributing to pathogenesis, shedding new light on the discovery of novel drug targets. The field of sequencing-based data mining is evolving toward a multi-omics approach. With advances in viral tools for precise microglial regulation and the development of brain organoid models, we are poised to elucidate the functional roles of microglial subgroups detected through sequencing analysis, ultimately identifying valuable therapeutic targets.

scCircle-seq unveils the diversity and complexity of extrachromosomal circular DNAs in single cells

Extrachromosomal circular DNAs (eccDNAs) have emerged as important intra-cellular mobile genetic elements that affect gene copy number and exert in trans regulatory roles within the cell nucleus. Here, we describe scCircle-seq, a method for profiling eccDNAs and unraveling their diversity and complexity in single cells. We implement and validate scCircle-seq in normal and cancer cell lines, demonstrating that most eccDNAs vary largely between cells and are stochastically inherited during cell division, although their genomic landscape is cell type-specific and can be used to accurately cluster cells of the same origin. eccDNAs are preferentially produced from chromatin regions enriched in H3K9me3 and H3K27me3 histone marks and are induced during replication stress conditions. Concomitant sequencing of eccDNAs and RNA from the same cell uncovers the absence of correlation between eccDNA copy number and gene expression levels, except for a few oncogenes, including MYC, contained within a large eccDNA in colorectal cancer cells. Lastly, we apply scCircle-seq to one prostate cancer and two breast cancer specimens, revealing cancer-specific eccDNA landscapes and a higher propensity of eccDNAs to form in amplified genomic regions. scCircle-seq is a scalable tool that can be used to dissect the complexity of eccDNAs across different cell and tissue types, and further expands the potential of eccDNAs for cancer diagnostics.

Single-Cell RNA Sequencing Revealed That the Enrichment of TPI1+ Malignant Hepatocytes Was Linked to HCC Metastasis and Immunosuppressive Microenvironment

Background

Tumor metastasis is the leading cause of high mortality in hepatocellular carcinoma (HCC). The metastasis-related HCC microenvironment is characterized by high heterogeneity. Single-cell RNA sequencing (scRNA-seq) may aid in determining specific cell clusters involved in regulating the immune microenvironment of HCC.

Methods

The scRNA-seq data of 10 HCC samples were collected from the Gene Expression Omnibus (GEO) database GSE124395. Correlations between key gene expression and clinicopathological data were determined using public databases. HCC tissues and matched tumor-adjacent and normal tissue samples were obtained by surgical resection at Sichuan Cancer Hospital. Immune cell infiltration analysis was performed and verified by immunohistochemistry and immunofluorescent staining.

Results

Nine malignant hepatocyte clusters with different marker genes and biological functions were identified. C3_Hepatocyte-SERF2 and C6_Hepatocyte-IL13RA2 were mainly involved in the regulation of the immune microenvironment, which was also a significant pathway in regulating HCC metastasis. Key genes in malignant hepatocyte clusters that associated with HCC metastasis were further screened by LASSO regression analysis. TPI1, a key gene in C6_Hepatocyte-IL13RA2 and HCC metastasis, could participate in regulating the HCC immune microenvironment in The Cancer Genome Atlas (TCGA) and Tumor Immune Estimation Resource (TIMER) databases. Moreover, immunohistochemistry analysis demonstrated that TPI1 expression was positively correlated with HCC metastasis and poor prognosis, while negatively correlated with CD8+ T cell infiltration. The negative correlation between TPI1 expression and CD8+ T cell infiltration was further confirmed by immunofluorescence staining.

Conclusion

In summary, a cluster of TPI1+ malignant hepatocytes was associated with the suppression of CD8+ T cell infiltration and HCC metastasis, providing novel insights into potential biomarkers for immunotherapy in HCC.

Single-cell RNA sequencing of anaplastic ependymoma and H3K27M-mutant diffuse midline glioma

BACKGROUND

Anaplastic ependymoma and H3K27M-mutant diffuse midline glioma are two common subtypes of brain tumors with poor long-term prognosis. The present study analyzed and compared the differences in cell types between two tumors by single-cell RNA sequencing (scRNA-seq) technology.

METHODS

ScRNA-seq was performed to profile cells from cancer tissue from anaplastic ependymoma patient and H3K27M-mutant diffuse midline glioma patient. Cell clustering, marker gene identification, cell type annotation, copy number variation analysis and function analysis of differentially expressed genes were then performed.

RESULTS

A total of 11,219 cells were obtained from anaplastic ependymoma and H3K27M mutant diffuse midline glioma, and these cells categorized into 12 distinct clusters. Each cell cluster could be characterized with specific cell markers to indicate cellular heterogeneity. Five cell types were annotated in each sample, including astrocyte, oligodendrocytes, microglial cell, neural progenitor cell and immune cell. The cluster types and proportion of cell types were not consistent between the two brain tumors. Functional analyses suggest that these cell clusters are involved in tumor-associated pathways, with slight differences in the cells of origin between the two tumors. In addition, cell communication analysis showed that the NRG3-ERBB4 pair is a key Ligand-receptor pair for anaplastic ependymoma, while in H3K27M-mutant diffuse midline glioma it is the PTN-PTPRZ1 pair that establishes contact with other cells.

CONCLUSION

There was intratumor heterogeneity in anaplastic ependymoma and H3K27M mutant diffuse midline glioma, and that the subtype differences may be due to differences in the origin of the cells.

A statistical learning method for simultaneous copy number estimation and subclone clustering with single-cell sequencing data

The availability of single-cell sequencing (SCS) enables us to assess intra-tumor heterogeneity and identify cellular subclones without the confounding effect of mixed cells. Copy number aberrations (CNAs) have been commonly used to identify subclones in SCS data using various clustering methods, as cells comprising a subpopulation are found to share a genetic profile. However, currently available methods may generate spurious results (e.g., falsely identified variants) in the procedure of CNA detection, thereby diminishing the accuracy of subclone identification within a large, complex cell population. In this study, we developed a subclone clustering method based on a fused lasso model, referred to as FLCNA, which can simultaneously detect CNAs in single-cell DNA sequencing (scDNA-seq) data. Spike-in simulations were conducted to evaluate the clustering and CNA detection performance of FLCNA, benchmarking it against existing copy number estimation methods (SCOPE, HMMcopy) in combination with commonly used clustering methods. Application of FLCNA to a scDNA-seq data set of breast cancer revealed different genomic variation patterns in neoadjuvant chemotherapy-treated samples and pretreated samples. We show that FLCNA is a practical and powerful method for subclone identification and CNA detection with scDNA-seq data.

Paradigm Shift: A Comprehensive Review of Ovarian Cancer Management in an Era of Advancements

Ovarian cancer (OC) is the female genital malignancy with the highest lethality. Patients present a poor prognosis mainly due to the late clinical presentation allied with the common acquisition of chemoresistance and a high rate of tumour recurrence. Effective screening, accurate diagnosis, and personalised multidisciplinary treatments are crucial for improving patients' survival and quality of life. This comprehensive narrative review aims to describe the current knowledge on the aetiology, prevention, diagnosis, and treatment of OC, highlighting the latest significant advancements and future directions. Traditionally, OC treatment involves the combination of cytoreductive surgery and platinum-based chemotherapy. Although more therapeutical approaches have been developed, the lack of established predictive biomarkers to guide disease management has led to only marginal improvements in progression-free survival (PFS) while patients face an increasing level of toxicity. Fortunately, because of a better overall understanding of ovarian tumourigenesis and advancements in the disease's (epi)genetic and molecular profiling, a paradigm shift has emerged with the identification of new disease biomarkers and the proposal of targeted therapeutic approaches to postpone disease recurrence and decrease side effects, while increasing patients' survival. Despite this progress, several challenges in disease management, including disease heterogeneity and drug resistance, still need to be overcome.

Single-cell transcriptomic analysis reveals tumor cell heterogeneity and immune microenvironment features of pituitary neuroendocrine tumors

BACKGROUND

Pituitary neuroendocrine tumors (PitNETs) are one of the most common types of intracranial tumors. Currently, the cellular characteristics of normal pituitary and various other types of PitNETs are still not completely understood.

METHODS

We performed single-cell RNA sequencing (scRNA-seq) on 4 normal samples and 24 PitNET samples for comprehensive bioinformatics analysis. Findings regarding the function of PBK in the aggressive tumor cells were validated by siRNA knockdown, overexpression, and transwell experiments.

RESULTS

We first constructed a reference cell atlas of the human pituitary. Subsequent scRNA-seq analysis of PitNET samples, representing major tumor subtypes, shed light on the intrinsic cellular heterogeneities of the tumor cells and tumor microenvironment (TME). We found that the expression of hormone-encoding genes defined the major variations of the PIT1-lineage tumor cell transcriptomic heterogeneities. A sub-population of TPIT-lineage tumor cells highly expressing GZMK suggested a novel subtype of corticotroph tumors. In immune cells, we found two clusters of tumor-associated macrophages, which were both highly enriched in PitNETs but with distinct functional characteristics. In PitNETs, the stress response pathway was significantly activated in T cells. While a majority of these tumors are benign, our study unveils a common existence of aggressive tumor cells in the studied samples, which highly express a set of malignant signature genes. The following functional experiments confirmed the oncogenic role of selected up-regulated genes. The over-expression of PBK could promote both tumor cell proliferation and migration, and it was also significantly associated with poor prognosis in PitNET patients.

CONCLUSIONS

Our data and analysis manifested the basic cell types in the normal pituitary and inherent heterogeneity of PitNETs, identified several features of the tumor immune microenvironments, and found a novel epithelial cell sub-population with aggressive signatures across all the studied cases.

Single-cell sequencing reveals SATB2/NOTCH1 signaling promotes the progression of malignancy of epithelial cells from papillary thyroid cancer

Although most papillary thyroid cancers (PTCs) are deemed to have a favorable clinical course and outcome, some develop an aggressive biological behavior at diagnosis or during treatment. Single-cell RNA sequencing (scRNA-seq), which is based on quantifying the features of individual cells to resolve tumor tissue heterogeneity, was used to uncover gene regulatory relationships and trace the transcriptional trajectories underlying the malignant transformation. In this study, we performed single-cell sequencing on samples from four PTC patients and one benign thyroid tumor patient. These included two papillary thyroid microcarcinoma cancers (PTMC) patients, two age-matched advanced PTC patients with invading surrounding tissues, and one patient undergoing surgical treatment due to a benign thyroid tumor. We constructed a new PTC RNA spectrum using single-cell sequencing. Single-cell sequencing analysis indicated that there was a highly invasive subgroup in the PTC epithelial cells, the expression of SATB2 (special AT-rich binding protein-2) was related to the prognosis and clinical progress of PTC, and SATB2 could promote the proliferation, migration, and invasion of PTC cells. We found that NOTCH1 was the key target gene of SATB2, and the activation of the SATB2/NOTCH1 pathway was one of the reasons for the high carcinogenicity of this subgroup.

scInterpreter: a knowledge-regularized generative model for interpretably integrating scRNA-seq data

BACKGROUND

The rapid emergence of single-cell RNA-seq (scRNA-seq) data presents remarkable opportunities for broad investigations through integration analyses. However, most integration models are black boxes that lack interpretability or are hard to train.

RESULTS

To address the above issues, we propose scInterpreter, a deep learning-based interpretable model. scInterpreter substantially outperforms other state-of-the-art (SOTA) models in multiple benchmark datasets. In addition, scInterpreter is extensible and can integrate and annotate atlas scRNA-seq data. We evaluated the robustness of scInterpreter in a variety of situations. Through comparison experiments, we found that with a knowledge prior, the training process can be significantly accelerated. Finally, we conducted interpretability analysis for each dimension (pathway) of cell representation in the embedding space.

CONCLUSIONS

The results showed that the cell representations obtained by scInterpreter are full of biological significance. Through weight sorting, we found several new genes related to pathways in PBMC dataset. In general, scInterpreter is an effective and interpretable integration tool. It is expected that scInterpreter will bring great convenience to the study of single-cell transcriptomics.

Unraveling the enigma of B cells in diffuse large B-cell lymphoma: unveiling cancer stem cell-like B cell subpopulation at single-cell resolution

Background

Diffuse large B-cell lymphoma (DLBCL) represents the most prevalent form of aggressive non-Hodgkin lymphoma. Despite receiving standard treatment, a subset of patients undergoes refractory or recurrent cases, wherein the involvement of cancer stem cells (CSCs) could be significant.

Methods

We comprehensively characterized B cell subpopulations using single-cell RNA sequencing data from three DLBCL samples and one normal lymph tissue. The CopyKat R package was employed to assess the malignancy of B cell subpopulations based on chromosomal copy number variations. CIBERSORTx software was utilized to estimate the proportions of B cell subpopulations in 230 DLBCL tissues. Furthermore, we employed the pySCENIC to identify key transcription factors that regulate the functionality of B cell subpopulations. By employing CellphoneDB, we elucidated the interplay among tumor microenvironment components within the B cell subpopulations. Finally, we validated our findings through immunofluorescence experiments.

Results

Our analysis revealed a specific cancer stem cell-like B cell subpopulation exhibiting self-renewal and multilineage differentiation capabilities based on the exploration of B cell subpopulations in DLBCL and normal lymph tissues at the single-cell level. Notably, a high infiltration of cancer stem cell-like B cells correlated with a poor prognosis, potentially due to immune evasion mediated by low expression of major histocompatibility complex molecules. Furthermore, we identified key transcription factor regulatory networks regulated by HMGB3, SAP30, and E2F8, which likely played crucial roles in the functional characterization of the cancer stem cell-like B cell subpopulation. The existence of cancer stem cell-like B cells in DLBCL was validated through immunofluorescent staining. Finally, cell communication between B cells and tumor-infiltrating T cell subgroups provided further insights into the functional characterization of the cancer stem cell-like B cell subpopulation.

Conclusions

Our research provides a systematic description of a specific cancer stem cell-like B cell subpopulation associated with a poor prognosis in DLBCL. This study enhances our understanding of CSCs and identifies potential therapeutic targets for refractory or recurrent DLBCL patients.

Single cell multi-omics reveal intra-cell-line heterogeneity across human cancer cell lines

Human cancer cell lines have long served as tools for cancer research and drug discovery, but the presence and the source of intra-cell-line heterogeneity remain elusive. Here, we perform single-cell RNA-sequencing and ATAC-sequencing on 42 and 39 human cell lines, respectively, to illustrate both transcriptomic and epigenetic heterogeneity within individual cell lines. Our data reveal that transcriptomic heterogeneity is frequently observed in cancer cell lines of different tissue origins, often driven by multiple common transcriptional programs. Copy number variation, as well as epigenetic variation and extrachromosomal DNA distribution all contribute to the detected intra-cell-line heterogeneity. Using hypoxia treatment as an example, we demonstrate that transcriptomic heterogeneity could be reshaped by environmental stress. Overall, our study performs single-cell multi-omics of commonly used human cancer cell lines and offers mechanistic insights into the intra-cell-line heterogeneity and its dynamics, which would serve as an important resource for future cancer cell line-based studies.

scMuffin: an R package to disentangle solid tumor heterogeneity by single-cell gene expression analysis

INTRODUCTION

Single-cell (SC) gene expression analysis is crucial to dissect the complex cellular heterogeneity of solid tumors, which is one of the main obstacles for the development of effective cancer treatments. Such tumors typically contain a mixture of cells with aberrant genomic and transcriptomic profiles affecting specific sub-populations that might have a pivotal role in cancer progression, whose identification eludes bulk RNA-sequencing approaches. We present scMuffin, an R package that enables the characterization of cell identity in solid tumors on the basis of a various and complementary analyses on SC gene expression data.

RESULTS

scMuffin provides a series of functions to calculate qualitative and quantitative scores, such as: expression of marker sets for normal and tumor conditions, pathway activity, cell state trajectories, Copy Number Variations, transcriptional complexity and proliferation state. Thus, scMuffin facilitates the combination of various evidences that can be used to distinguish normal and tumoral cells, define cell identities, cluster cells in different ways, link genomic aberrations to phenotypes and identify subtle differences between cell subtypes or cell states. We analysed public SC expression datasets of human high-grade gliomas as a proof-of-concept to show the value of scMuffin and illustrate its user interface. Nevertheless, these analyses lead to interesting findings, which suggest that some chromosomal amplifications might underlie the invasive tumor phenotype and the presence of cells that possess tumor initiating cells characteristics.

CONCLUSIONS

The analyses offered by scMuffin and the results achieved in the case study show that our tool helps addressing the main challenges in the bioinformatics analysis of SC expression data from solid tumors.

Integrative Analysis of Single-Cell and Bulk RNA Sequencing Reveals Prognostic Characteristics of Macrophage Polarization-Related Genes in Lung Adenocarcinoma

Background

Lung adenocarcinoma (LUAD) is a group of cancers with poor prognosis. The combination of single-cell RNA sequencing (scRNA-seq) and bulk RNA sequencing (RNA-seq) can identify important genes involved in cancer development and progression from a broader perspective.

Methods

The scRNA-seq data and bulk RNA-seq data of LUAD were downloaded from the Gene Expression Omnibus (GEO) database and the Cancer Genome Atlas (TCGA) database. Analyzing scRNA-seq for core cells in the GSE131907 dataset, and the uniform manifold approximation and projection (UMAP) was used for dimensionality reduction and cluster identification. Macrophage polarization-associated subtypes were acquired from the TCGA-LUAD dataset after analysis, followed by further identification of differentially expressed genes (DEGs) in the TCGA-LUAD dataset (normal/LUAD tissue samples, two subtypes). Venn diagrams were utilized to visualize differentially expressed and highly variable macrophage polarization-related genes. Subsequently, a prognostic risk model for LUAD patients was constructed by univariate Cox and Least Absolute Shrinkage and Selection Operator (LASSO), and the model was investigated for stability in the external data GSE72094. After analyzing the correlation between the trait genes and significantly mutated genes, the immune infiltration between the high/low-risk groups was then examined. The Monocle package was applied to analyze the pseudo-temporal trajectory analysis of different cell clusters in macrophage clusters. Subsequently, cell clusters of data macrophages were selected as key cell clusters to explore the role of characteristic genes in different cell populations and to identify transcription factors (TFs) that affect signature genes. Finally, qPCR were employed to validate the expression levels of prognosis signature genes in LUAD.

Results

424 macrophage highly variable genes, 3920 DEGs, and 9561 DEGs were obtained from macrophage clusters, the macrophage polarization-related subtypes, and normal/LUAD tissue samples, respectively. Twenty-eight differentially expressed and highly mutated MPRGs were obtained. A prognostic risk model with 7 DE-MPRGs (RGS13, ADRB2, DDIT4, MS4A2, ALDH2, CTSH, and PKM) was constructed. This prognostic model still has a good prediction effect in the GSE72094 dataset. ZNF536 and DNAH9 were mutated in the low-risk group, while COL11A1 was mutated in the high-risk group, and they were highly correlated with the characteristic genes. A total of 11 immune cells were significantly different in the high/low-risk groups. Five cell types were again identified in the macrophage cluster, and then NK cells: CD56hiCD62L+ differentiated earlier and were present mainly on 2 branches. While macrophages were present on 2 branches and differentiated later. It was found that the expression levels of BCLAF1 and MAX were higher in cluster 1, which might be the TFs affecting the expression of the characteristic genes. Moreover, qPCR confirmed that the expression of the prognosis genes was generally consistent with the results of the bioinformatic analysis.

Conclusion

Seven MPRGs (RGS13, ADRB2, DDIT4, MS4A2, ALDH2, CTSH, and PKM) were identified as prognostic genes for LUAD and revealed the mechanisms of MPRGs at the single-cell level.

Application of single-cell sequencing to the research of tumor microenvironment

Single-cell sequencing is a technique for detecting and analyzing genomes, transcriptomes, and epigenomes at the single-cell level, which can detect cellular heterogeneity lost in conventional sequencing hybrid samples, and it has revolutionized our understanding of the genetic heterogeneity and complexity of tumor progression. Moreover, the tumor microenvironment (TME) plays a crucial role in the formation, development and response to treatment of tumors. The application of single-cell sequencing has ushered in a new age for the TME analysis, revealing not only the blueprint of the pan-cancer immune microenvironment, but also the heterogeneity and differentiation routes of immune cells, as well as predicting tumor prognosis. Thus, the combination of single-cell sequencing and the TME analysis provides a unique opportunity to unravel the molecular mechanisms underlying tumor development and progression. In this review, we summarize the recent advances in single-cell sequencing and the TME analysis, highlighting their potential applications in cancer research and clinical translation.

Spatial downregulation of CD74 signatures may drive invasive component development in part-solid lung adenocarcinoma

Pulmonary nodules with part-solid imaging features manifest during the progression from preinvasive to invasive lung adenocarcinoma. To define the spatial composition and evolutionary trajectories of early-stage lung adenocarcinoma, we combined spatial transcriptomics (ST) and pathological annotations from 20 part-solid nodules (PSNs), four of which were matched with single-cell RNA sequencing. Two malignant cell populations (MC1 and MC2) were identified, and a linear evolutionary relationship was observed. Compared to MC2, the pre-existing malignant MC1 exhibited a lower metastatic signature, corresponding to the preinvasive component (lepidic) on pathology and the ground glass component on PSN imaging. Higher immune infiltration was observed among MC1 regions in ST profiles, and further analysis revealed that macrophages may be involved in this process through the CD74 axis. This work provides deeper insights into the evolutionary process and spatial immune cell composition behind PSNs and highlights the mechanisms of immune escape behind this adenocarcinoma trajectory.

High-dimensional deconstruction of pancreatic cancer identifies tumor microenvironmental and developmental stemness features that predict survival

Numerous cell states are known to comprise the pancreatic ductal adenocarcinoma (PDAC) tumor microenvironment (TME). However, the developmental stemness and co-occurrence of these cell states remain poorly defined. Here, we performed single-cell RNA sequencing (scRNA-seq) on a cohort of treatment-naive PDAC time-of-diagnosis endoscopic ultrasound-guided fine needle biopsy (EUS-FNB) samples (n = 25). We then combined these samples with surgical resection (n = 6) and publicly available samples to increase statistical power (n = 80). Following annotation into 25 distinct cell states, cells were scored for developmental stemness, and a customized version of the Ecotyper tool was used to identify communities of co-occurring cell states in bulk RNA-seq samples (n = 268). We discovered a tumor microenvironmental community comprised of aggressive basal-like malignant cells, tumor-promoting SPP1+ macrophages, and myofibroblastic cancer-associated fibroblasts associated with especially poor prognosis. We also found a developmental stemness continuum with implications for survival that is present in both malignant cells and cancer-associated fibroblasts (CAFs). We further demonstrated that high-dimensional analyses predictive of survival are feasible using standard-of-care, time-of-diagnosis EUS-FNB specimens. In summary, we identified tumor microenvironmental and developmental stemness characteristics from a high-dimensional gene expression analysis of PDAC using human tissue specimens, including time-of-diagnosis EUS-FNB samples. These reveal new connections between tumor microenvironmental composition, CAF and malignant cell stemness, and patient survival that could lead to better upfront risk stratification and more personalized upfront clinical decision-making.

Unveiling the Role of Tumor-Infiltrating T Cells and Immunotherapy in Hepatocellular Carcinoma: A Comprehensive Review

Hepatocellular carcinoma (HCC) is a rapidly rising global health concern, ranking as the third-leading cause of cancer-related mortality. Despite medical advancements, the five-year survival rate remains a dismal 18%, with a daunting 70% recurrence rate within a five-year period. Current systematic treatments, including first-line sorafenib, yield an overall response rate (ORR) below 10%. In contrast, immunotherapies have shown promise by improving ORR to approximately 30%. The IMbravel150 clinical trial demonstrates that combining atezolizumab and bevacizumab surpasses sorafenib in terms of median progression-free survival (PFS) and overall survival (OS). However, the therapeutic efficacy for HCC patients remains unsatisfactory, highlighting the urgent need for a comprehensive understanding of antitumor responses and immune evasion mechanisms in HCC. In this context, understanding the immune landscape of HCC is of paramount importance. Tumor-infiltrating T cells, including cytotoxic T cells, regulatory T cells, and natural killer T cells, are key components in the antitumor immune response. This review aims to shed light on their intricate interactions within the immunosuppressive tumor microenvironment and explores potential strategies for revitalizing dysfunctional T cells. Additionally, current immune checkpoint inhibitor (ICI)-based trials, ICI-based combination therapies, and CAR-T- or TCR-T-cell therapies for HCC are summarized, which might further improve OS and transform the management of HCC in the future.

Mapping pediatric brain tumors to their origins in the developing cerebellum

BACKGROUND

Distinguishing the cellular origins of childhood brain tumors is key for understanding tumor initiation and identifying lineage-restricted, tumor-specific therapeutic targets. Previous strategies to map the cell-of-origin typically involved comparing human tumors to murine embryonal tissues, which is potentially limited due to species-specific differences. The aim of this study was to unravel the cellular origins of the 3 most common pediatric brain tumors, ependymoma, pilocytic astrocytoma, and medulloblastoma, using a developing human cerebellar atlas.

METHODS

We used a single-nucleus atlas of the normal developing human cerebellum consisting of 176 645 cells as a reference for an in-depth comparison to 4416 bulk and single-cell transcriptome tumor datasets, using gene set variation analysis, correlation, and single-cell matching techniques.

RESULTS

We find that the astroglial cerebellar lineage is potentially the origin for posterior fossa ependymomas. We propose that infratentorial pilocytic astrocytomas originate from the oligodendrocyte lineage and MHC II genes are specifically enriched in these tumors. We confirm that SHH and Group 3/4 medulloblastomas originate from the granule cell and unipolar brush cell lineages. Radiation-induced gliomas stem from cerebellar glial lineages and demonstrate distinct origins from the primary medulloblastoma. We identify tumor genes that are expressed in the cerebellar lineage of origin, and genes that are tumor specific; both gene sets represent promising therapeutic targets for future study.

CONCLUSION

Based on our results, individual cells within a tumor may resemble different cell types along a restricted developmental lineage. Therefore, we suggest that tumors can arise from multiple cellular states along the cerebellar "lineage of origin."

Neoadjuvant chemotherapy for breast cancer: an evaluation of its efficacy and research progress

Neoadjuvant chemotherapy (NAC) for breast cancer is widely used in the clinical setting to improve the chance of surgery, breast conservation and quality of life for patients with advanced breast cancer. A more accurate efficacy evaluation system is important for the decision of surgery timing and chemotherapy regimen implementation. However, current methods, encompassing imaging techniques such as ultrasound and MRI, along with non-imaging approaches like pathological evaluations, often fall short in accurately depicting the therapeutic effects of NAC. Imaging techniques are subjective and only reflect macroscopic morphological changes, while pathological evaluation is the gold standard for efficacy assessment but has the disadvantage of delayed results. In an effort to identify assessment methods that align more closely with real-world clinical demands, this paper provides an in-depth exploration of the principles and clinical applications of various assessment approaches in the neoadjuvant chemotherapy process.

The prognostic and biology of tumour-infiltrating lymphocytes in the immunotherapy of cancer

Tumour immunotherapy has achieved remarkable clinical success in many different types of cancer in the past two decades. The outcome of immune checkpoint inhibitors in cancer patients has been linked to the quality and magnitude of T cell, NK cell, and more recently, B cell within the tumour microenvironment, suggesting that the immune landscape of a tumour is highly connected to patient response and prognosis. It is critical to understanding tumour immune microenvironments for identifying immune modifiers of cancer progression and developing cancer immunotherapies. The infiltration of solid tumours by immune cells with anti-tumour activity is both a strong prognostic factor and a therapeutic goal. Recent approaches and applications of new technologies, especially single-cell mRNA analysis in dissecting tumour microenvironments have brought important insights into the biology of tumour-infiltrating immune cells, revealed a remarkable degree of cellular heterogeneity and distinct patterns of immune response. In this review, we will discuss recent advances in the understanding of tumour infiltrated lymphocytes, their prognostic benefit, and predictive value for immunotherapy.

Follicular lymphoma microenvironment: insights provided by single-cell analysis

Follicular lymphoma (FL) is the most frequent indolent lymphoma and is characterized by the abundant infiltration of tumor microenvironment (TME) cells. The activity of TME cells reportedly plays an important role in the biology of FL. TME cells that reside within neoplastic follicles, such as T-follicular helper cells and follicular dendritic cells, have been shown to aid in FL development and progression through interactions with malignant B cells, whereas regulatory T cells have unexpectedly shown an apparently favorable prognostic impact in FL. Unfortunately, the understanding of the FL TME, particularly regarding minor cell subsets, has been hampered by unknown cell heterogeneity. As with other solid and hematologic cancers, novel single-cell analysis technologies have recently been applied to FL research and have uncovered previously unrecognized heterogeneities, not only in malignant B cells but also in TME cells. These reports have greatly increased the resolution of our understanding of the FL TME and, at the same time, raised questions about newly identified TME cells. This review provides an overview of the unique aspects of FL TME cells with a clinical viewpoint and highlights recent discoveries from single-cell analysis, while also suggesting potential future directions.

Development and validation of cuproptosis-related lncRNAs associated with pancreatic cancer immune microenvironment based on single-cell

Background

Cuproptosis, a novel mode of cell death associated with the tricarboxylic acid (TCA) cycle, is relevant to the development of cancer. However, the impact of single-cell-based Cuproptosis-associated lncRNAs on the Tumor immune microenvironment (TIME) of Pancreatic adenocarcinoma (PAAD) and its potential value for individualized immunotherapy has not been clarified.

Methods

14 immune-related CRGs were screened by exploring the interaction between differentially expressed Immune-Related Genes (IRGs) and Cuproptosis-Related Genes (CRGs) in PAAD. Next, the expression amount and expression distribution of CRGs in single-cell samples were analyzed by focusing on 7-CRGs with significant expressions. On the one hand, MAP2K2, SOD1, and VEGFA, which were significantly differentially expressed between PAAD sites and normal tissues adjacent to them, were subjected to immunohistochemical validation and immune landscape analysis. On the other hand, from these 7-CRGs, prognostic signatures of lncRNAs were established by co-expression and LASSO-COX regression analysis, and their prognostic value and immune relevance were assessed. In addition, this study not only validated the hub CRGs and the lncRNAs constituting the signature in a PAAD animal model treated with immunotherapy-based combination therapy using immunohistochemistry and qRT-PCR but also explored the potential value of the combination of targeted, chemotherapy and immunotherapy.

Results

Based on the screening of 7-CRGs significantly expressed in a PAAD single-cell cohort and their co-expressed Cuproptosis-Related lncRNAs (CRIs), this study constructed a prognostic signature of 4-CRIs named CIR-score. A Nomogram integrating the CIR-score and clinical risk factors was constructed on this basis to predict the individualized survival of patients. Moreover, high and low-risk groups classified according to the median of signatures exhibited significant differences in clinical prognosis, immune landscape, bioenrichment, tumor burden, and drug sensitivity. And the immunohistochemical and qRT-PCR results of different mouse PAAD treatment strategies were consistent with the trend of inter-group variability in drug sensitivity of hub CRGs and CIR-score. The combination of immunotherapy, targeted therapy, and chemotherapy exhibited a better tumor suppression effect.

Conclusion

CIR-score, as a Cuproptosis-related TIME-specific prognostic signature based on PAAD single cells, not only predicts the prognosis and immune landscape of PAAD patients but also provides a new strategy for individualized immunotherapy-based combination therapy.

Microfluidic Isolation of Disseminated Tumor Cells from the Bone Marrow of Breast Cancer Patients

Disseminated tumor cells (DTCs) in the bone marrow (BM) of breast cancer (BC) patients are putative precursors of metastatic disease, and their presence is associated with an adverse clinical outcome. To achieve the personalization of therapy on a clinical routine level, the characterization of DTCs and in vitro drug testing on DTCs are of great interest. Therefore, biobanking methods, as well as novel approaches to DTC isolation, need to be developed. In this study, we established a protocol for the biobanking of BM samples and evaluated a microfluidic-based separation system (Parsortix®) for the enrichment of cryopreserved DTCs. We were able to successfully isolate viable DTCs after the prior cryopreservation of BM samples. We calculated a significant increase of up to 90-fold in harvested DTCs with the proposed method compared to the current standard techniques, opening up new analysis possibilities for DTCs. Our advanced method further presents options for 3D DTC cultures, enabling the individualized testing of targeted therapies for BC patients. In conclusion, we present a novel approach for DTC enrichment, with possibilities for future clinical implications.

CellAnn: a comprehensive, super-fast, and user-friendly single-cell annotation web server

MOTIVATION

Single-cell sequencing technology has become a routine in studying many biological problems. A core step of analyzing single-cell data is the assignment of cell clusters to specific cell types. Reference-based methods are proposed for predicting cell types for single-cell clusters. However, the scalability and lack of preprocessed reference datasets prevent them from being practical and easy to use.

RESULTS

Here, we introduce a reference-based cell annotation web server, CellAnn, which is super-fast and easy to use. CellAnn contains a comprehensive reference database with 204 human and 191 mouse single-cell datasets. These reference datasets cover 32 organs. Furthermore, we developed a cluster-to-cluster alignment method to transfer cell labels from the reference to the query datasets, which is superior to the existing methods with higher accuracy and higher scalability. Finally, CellAnn is an online tool that integrates all the procedures in cell annotation, including reference searching, transferring cell labels, visualizing results, and harmonizing cell annotation labels. Through the user-friendly interface, users can identify the best annotation by cross-validating with multiple reference datasets. We believe that CellAnn can greatly facilitate single-cell sequencing data analysis.

AVAILABILITY AND IMPLEMENTATION

The web server is available at www.cellann.io, and the source code is available at https://github.com/Pinlyu3/CellAnn_shinyapp.

Automated microarray platform for single-cell sorting and collection of lymphocytes following HIV reactivation

A promising strategy to cure HIV-infected individuals is to use latency reversing agents (LRAs) to reactivate latent viruses, followed by host clearance of infected reservoir cells. However, reactivation of latent proviruses within infected cells is heterogeneous and often incomplete. This fact limits strategies to cure HIV which may require complete elimination of viable virus from all cellular reservoirs. For this reason, understanding the mechanism(s) of reactivation of HIV within cellular reservoirs is critical to achieve therapeutic success. Methodologies enabling temporal tracking of single cells as they reactivate followed by sorting and molecular analysis of those cells are urgently needed. To this end, microraft arrays were adapted to image T-lymphocytes expressing mCherry under the control of the HIV long terminal repeat (LTR) promoter, in response to the application of LRAs (prostratin, iBET151, and SAHA). In response to prostratin, iBET151, and SAHA, 30.5%, 11.2%, and 12.1% percentage of cells, respectively. The arrays enabled large numbers of single cells (>25,000) to be imaged over time. mCherry fluorescence quantification identified cell subpopulations with differing reactivation kinetics. Significant heterogeneity was observed at the single-cell level between different LRAs in terms of time to reactivation, rate of mCherry fluorescence increase upon reactivation, and peak fluorescence attained. In response to prostratin, subpopulations of T lymphocytes with slow and fast reactivation kinetics were identified. Single T-lymphocytes that were either fast or slow reactivators were sorted, and single-cell RNA-sequencing was performed. Different genes associated with inflammation, immune activation, and cellular and viral transcription factors were found.

Barcoded multiple displacement amplification for high coverage sequencing in spatial genomics

Determining mutational landscapes in a spatial context is essential for understanding genetically heterogeneous cell microniches. Current approaches, such as Multiple Displacement Amplification (MDA), offer high genome coverage but limited multiplexing, which hinders large-scale spatial genomic studies. Here, we introduce barcoded MDA (bMDA), a technique that achieves high-coverage genomic analysis of low-input DNA while enhancing the multiplexing capabilities. By incorporating cell barcodes during MDA, bMDA streamlines library preparation in one pot, thereby overcoming a key bottleneck in spatial genomics. We apply bMDA to the integrative spatial analysis of triple-negative breast cancer tissues by examining copy number alterations, single nucleotide variations, structural variations, and kataegis signatures for each spatial microniche. This enables the assessment of subclonal evolutionary relationships within a spatial context. Therefore, bMDA has emerged as a scalable technology with the potential to advance the field of spatial genomics significantly.

How single-cell techniques help us look into lung cancer heterogeneity and immunotherapy

Lung cancer patients tend to have strong intratumoral and intertumoral heterogeneity and complex tumor microenvironment, which are major contributors to the efficacy of and drug resistance to immunotherapy. From a new perspective, single-cell techniques offer an innovative way to look at the intricate cellular interactions between tumors and the immune system and help us gain insights into lung cancer and its response to immunotherapy. This article reviews the application of single-cell techniques in lung cancer, with focuses directed on the heterogeneity of lung cancer and the efficacy of immunotherapy. This review provides both theoretical and experimental information for the future development of immunotherapy and personalized treatment for the management of lung cancer.

Spatial transcriptomics: recent developments and insights in respiratory research

The respiratory system's complex cellular heterogeneity presents unique challenges to researchers in this field. Although bulk RNA sequencing and single-cell RNA sequencing (scRNA-seq) have provided insights into cell types and heterogeneity in the respiratory system, the relevant specific spatial localization and cellular interactions have not been clearly elucidated. Spatial transcriptomics (ST) has filled this gap and has been widely used in respiratory studies. This review focuses on the latest iterative technology of ST in recent years, summarizing how ST can be applied to the physiological and pathological processes of the respiratory system, with emphasis on the lungs. Finally, the current challenges and potential development directions are proposed, including high-throughput full-length transcriptome, integration of multi-omics, temporal and spatial omics, bioinformatics analysis, etc. These viewpoints are expected to advance the study of systematic mechanisms, including respiratory studies.

An image-guided microfluidic system for single-cell lineage tracking

Cell lineage tracking is a long-standing and unresolved problem in biology. Microfluidic technologies have the potential to address this problem, by virtue of their ability to manipulate and process single-cells in a rapid, controllable and efficient manner. Indeed, when coupled with traditional imaging approaches, microfluidic systems allow the experimentalist to follow single-cell divisions over time. Herein, we present a valve-based microfluidic system able to probe the decision-making processes of single-cells, by tracking their lineage over multiple generations. The system operates by trapping single-cells within growth chambers, allowing the trapped cells to grow and divide, isolating sister cells after a user-defined number of divisions and finally extracting them for downstream transcriptome analysis. The platform incorporates multiple cell manipulation operations, image processing-based automation for cell loading and growth monitoring, reagent addition and device washing. To demonstrate the efficacy of the microfluidic workflow, 6C2 (chicken erythroleukemia) and T2EC (primary chicken erythrocytic progenitors) cells are tracked inside the microfluidic device over two generations, with a cell viability rate in excess of 90%. Sister cells are successfully isolated after division and extracted within a 500 nL volume, which was demonstrated to be compatible with downstream single-cell RNA sequencing analysis.