Single-cell transcriptomic analysis in a mouse model deciphers cell transition states in the multistep development of esophageal cancer

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.

Predictive model for acute radiation esophagitis in esophageal carcinoma based on prognostic nutritional index and systemic inflammatory index and its application

Acute radiation esophagitis (ARE) is a common complication in patients with esophageal cancer undergoing radiotherapy. Therefore, it is important to construct an effective ARE risk-prediction model for clinical treatment. The present study performed a retrospective analysis of 225 patients with esophageal cancer who received radiotherapy at the First Affiliated Hospital of Anhui Medical University (Hefei, China) from January 2018 to December 2022. Univariate and logistic regression analyses were performed to screen patients with esophageal cancer after radiotherapy. The results revealed that 147 patients developed radiation esophagitis. Logistic regression analysis results demonstrated that the prognostic nutritional index [odds ratio (OR), 0.864; 95% confidence interval (CI), 0.809-0.924], neutrophil to lymphocyte ratio (OR, 1.795; 95% CI, 1.209-2.667) and platelet to lymphocyte ratio (OR, 1.011; 95% CI, 1.000-1.022) were independent predictors of ARE in patients receiving intensity-modulated conformal radiotherapy for esophagus cancer (P<0.05). A nomogram model for predicting the occurrence of ARE was established based on the three risk factors. The decision curve suggested a high net benefit value when the threshold probability was within 0.25-1.0. External verification confirmed the reproducibility and generalizability of the nomogram model. In general, the calibration curve of this model was close to the ideal curve and had excellent prediction accuracy. Therefore, it may be used as a new tool for early prediction of the ARE risk.

Temporal single-cell RNA sequencing dataset of gastroesophagus development from embryonic to post-natal stages

Gastroesophageal disorders and cancers impose a significant global burden. Particularly, the prevalence of esophageal adenocarcinoma (EAC) has increased dramatically in recent years. Barrett's esophagus, a precursor of EAC, features a unique tissue adaptation at the gastroesophageal squamo-columnar junction (GE-SCJ), where the esophagus meets the stomach. Investigating the evolution of GE-SCJ and understanding dysregulation in its homeostasis are crucial for elucidating cancer pathogenesis. Here, we present the technical quality of the comprehensive single-cell RNA sequencing (scRNA-seq) dataset from mice that captures the transcriptional dynamics during the development of the esophagus, stomach and the GE-SCJ at embryonic, neonatal and adult stages. Through integration with external scRNA-seq datasets and validations using organoid and animal models, we demonstrate the dataset's consistency in identified cell types and transcriptional profiles. This dataset will be a valuable resource for studying developmental patterns and associated signaling networks in the tissue microenvironment. By offering insights into cellular programs during homeostasis, it facilitates the identification of changes leading to conditions like metaplasia and cancer, crucial for developing effective intervention strategies.

Emerging strategies to investigate the biology of early cancer

Early detection and intervention of cancer or precancerous lesions hold great promise to improve patient survival. However, the processes of cancer initiation and the normal-precancer-cancer progression within a non-cancerous tissue context remain poorly understood. This is, in part, due to the scarcity of early-stage clinical samples or suitable models to study early cancer. In this Review, we introduce clinical samples and model systems, such as autochthonous mice and organoid-derived or stem cell-derived models that allow longitudinal analysis of early cancer development. We also present the emerging techniques and computational tools that enhance our understanding of cancer initiation and early progression, including direct imaging, lineage tracing, single-cell and spatial multi-omics, and artificial intelligence models. Together, these models and techniques facilitate a more comprehensive understanding of the poorly characterized early malignant transformation cascade, holding great potential to unveil key drivers and early biomarkers for cancer development. Finally, we discuss how these new insights can potentially be translated into mechanism-based strategies for early cancer detection and prevention.

Landscape of epithelial cell subpopulations in the human esophageal squamous cell carcinoma microenvironment

Aims

We sought to reveal the landscape of epithelial cell subpopulations in the human esophageal squamous cell carcinoma microenvironment and investigate their parts on esophageal squamous carcinoma (ESCC) development.

Background

Epithelial cells play an important role in the occurrence and development of ESCC through multiple mechanisms. While the landscape of epithelial cell subpopulations in ESCC, remains unclear.

Objective

Exploring the role of epithelial cell subpopulations in ESCC progression.

Methods

Seurat R package was used for single-cell RNA sequencing (scRNA-seq) data filtering, dimensionality reduction, clustering and differentially expressed genes analysis. Cellmarker database was adopted for cell cluster annotation. Functional enrichment analysis was carried out by Gene Ontology (GO) analysis. InferCNV package was conducted for copy number variation (CNV) of epithelial cell subpopulations in all chromosomal regions. Pseudotime trajectory analysis was implemented for exploring differentiation trajectory of epithelial cells subgroups during the cancer progression. CellChat analysis was used for probing the interactions between epithelial cells and NK/T cells. cellular experiments were performed using Quantitative Real-Time Polymerase Chain Reaction (RT-qPCR), Wound-Healing Assay and transwell.

Results

11 major cell subpopulations were identified in ESCC and adjunct tissues. Further reclassification of epithelial cells uncovered 4 subpopulations. Enrichment analysis revealed that highly expressed genes in 4 epithelial cell subpopulations were related to cell proliferation, immune response and angiogenesis. CNV analysis found that UBD + epithelial cells and GAS2L3+ epithelial cells had a higher proportion of CNV. Cell differentiation trajectories disclosed that KRT6C+ and GSTA1+ epithelial cells were in an intermediate state of differentiation, while UBD+ and GAS2L3+ epithelial cells are in an end state of differentiation during ESCC progression. Finally, we found that four epithelial cell subpopulations all inhibited NK/T cells through NECTIN2-TIGIT and CLEC2B-KLRB1. Low ATF3 and DDIT3 mRNA expression inhibited ESCC cell migration and invasion.

Conclusion

Here, we obtained a through epithelial cell atlas of ESCC at single-cell resolution, explored the role of epithelial cell in ESCC progression, and unveiled immunosuppressive signals to NK/T cells in promoting ESCC. Our findings expand the comprehension of epithelial cells and offer a theoretical guidance for future anti-epithelial cell treatment of ESCC.

Ethanol exposure exacerbates 4-nitroquinoline-1-oxide induced esophageal carcinogenesis and induces invasive carcinoma with muscularis propria infiltration in a mouse model

Esophageal squamous cell carcinoma (ESCC) is one of the most fatal cancers worldwide. Most ESCC patients are diagnosed at an advanced stage; however, current research on in vivo animal models accurately reflecting their clinical presentation is lacking. Alcohol consumption is a major risk factor for ESCC and has been used in several disease models for disease induction. In this study, we used 4-nitroquinoline-1-oxide in combination with ethanol to induce an in vivo ESCC mouse model. Esophageal tissues were stained with hematoxylin and eosin for histopathological examination and lesion scoring. In cellular experiments, cell adhesion and migration invasion ability were observed using phalloidin staining, cell scratch and transwell assays, respectively, and the expression of epithelial-mesenchymal transition-related markers was detected using quantitative reverse transcription polymerase chain reaction and western blotting. The results showed that ethanol-exposed mice lost more weight and had an increased number of esophageal nodules. Histological examination revealed that the lesion scores of the ethanol-exposed esophageal samples were significantly higher than those of the unexposed esophageal samples. Furthermore, ethanol-exposed esophageal cancer samples had more severe lesions with infiltration of tumor cells into the muscularis propria. In vitro cellular experiments showed that ethanol exposure induced cytoskeletal microfilament formation, promoted cell migration invasion elevated the expression of N-cadherin and Snail, and decreased the expression of E-cadherin. In conclusion, ethanol exposure exacerbates ESCC, promotes tumor cell infiltration into the muscularis propria, and could be an effective agent for establishing innovative models of invasive carcinoma.

Thermosensitive gel-nano system against esophageal cancer via restoring p53 activity and boosting T-cell immunity

In esophageal cancer (EC), clinical specimen testing has uncovered a significant increase in BTB and CNC homolog 1 (BACH1) expression and a shift towards an immunosuppressive environment, alongside a notable decrease in p53 protein expression. Therefore, therapeutic strategies focusing on BACH1 inhibition and p53 upregulation appear promising. Traditional oral treatments for EC lack precision and efficacy. Here, we propose a novel approach employing tumor-targeted nanoparticles (NPs) for drug delivery. However, the formation of a drug reservoir at the esophageal site, crucial for the sustained release of therapeutics, presents significant challenges in nano-delivery systems for EC treatment. To address this, we developed a thermosensitive hydrogel composed of F127 and tannic acid, serving as a vehicle for NP loading. These NPs, synthesized through the emulsion/volatization methods of mPEG-PLGA-PLL-cRGD, facilitate in situ drug delivery. Upon contacting esophageal tissue, the hydrogel transitions to a gel, adhering to the lining and enabling sustained release of encapsulated therapeutics. The formulation encompasses NPs laden with small interfering RNA targeting BACH1 (siBACH1) and the p53 activator PRIMA-1, creating a cohesive gel-nano system. Preliminary biological assessments demonstrate that this injectable, thermosensitive gel-nano system adheres effectively to esophageal tissue and targets EC cells. For better modeling clinical outcomes, a patient-derived organoid xenograft (PDOX) model was innovated, involving transplantation of EC-derived organoids into humanized mice, reconstructed with peripheral blood mononuclear cells (PBMCs). Post-treatment analysis showed substantial EC growth inhibition (89.51% tumor inhibition rate), significant BACH1 level reduction, restored anti-tumor immune responses, and pronounced tumor apoptosis. In summary, our study introduces a thermosensitive gel-nano system for EC treatment via restoring p53 activity and boosting T-cell immunity, with potential for clinical application.

Computational single-cell methods for predicting cancer risk

Despite recent biotechnological breakthroughs, cancer risk prediction remains a formidable computational and experimental challenge. Addressing it is critical in order to improve prevention, early detection and survival rates. Here, I briefly summarize some key emerging theoretical and computational challenges as well as recent computational advances that promise to help realize the goals of cancer-risk prediction. The focus is on computational strategies based on single-cell data, in particular on bottom-up network modeling approaches that aim to estimate cancer stemness and dedifferentiation at single-cell resolution from a systems-biological perspective. I will describe two promising methods, a tissue and cell-lineage independent one based on the concept of diffusion network entropy, and a tissue and cell-lineage specific one that uses transcription factor regulons. Application of these tools to single-cell and single-nucleus RNA-seq data from stages prior to invasive cancer reveal that they can successfully delineate the heterogeneous inter-cellular cancer-risk landscape, identifying those cells that are more likely to turn cancerous. Bottom-up systems biological modeling of single-cell omic data is a novel computational analysis paradigm that promises to facilitate the development of preventive, early detection and cancer-risk prediction strategies.

Alarmin S100A8 imparts chemoresistance of esophageal cancer by reprogramming cancer-associated fibroblasts

Chemotherapy remains the first-line treatment for advanced esophageal cancer. However, durable benefits are achieved by only a limited subset of individuals due to the elusive chemoresistance. Here, we utilize patient-derived xenografts (PDXs) from esophageal squamous-cell carcinoma to investigate chemoresistance mechanisms in preclinical settings. We observe that activated cancer-associated fibroblasts (CAFs) are enriched in the tumor microenvironment of PDXs resistant to chemotherapy. Mechanistically, we reveal that cancer-cell-derived S100A8 triggers the intracellular RhoA-ROCK-MLC2-MRTF-A pathway by binding to the CD147 receptor of CAFs, inducing CAF polarization and leading to chemoresistance. Therapeutically, we demonstrate that blocking the S100A8-CD147 pathway can improve chemotherapy efficiency. Prognostically, we found the S100A8 levels in peripheral blood can serve as an indicator of chemotherapy responsiveness. Collectively, our study offers a comprehensive understanding of the molecular mechanisms underlying chemoresistance in esophageal cancer and highlights the potential value of S100A8 in the clinical management of esophageal cancer.

Tumor cell enrichment by tissue suspension improves sensitivity to copy number variation in diffuse gastric cancer with low tumor content

The detection of copy number variations (CNVs) and somatic mutations in cancer is important for the selection of specific drugs for patients with cancer. In cancers with sporadic tumor cells, low tumor content prevents the accurate detection of somatic alterations using targeted sequencing. To efficiently identify CNVs, we performed tumor cell enrichment using tissue suspensions of formalin-fixed paraffin-embedded (FFPE) tissue sections with low tumor cell content. Tumor-enriched and residual fractions were separated from FFPE tissue suspensions of intestinal and diffuse-type gastric cancers containing sporadic tumor cells, and targeted sequencing was performed on 225 cancer-related genes. Sequencing of a targeted panel of cancer-related genes using tumor-enriched fractions increased the number of detectable CNVs and the copy number of amplified genes. Furthermore, CNV analysis using the normal cell-enriched residual fraction as a reference for CNV scoring allowed targeted sequencing to detect CNV characteristics of diffuse-type gastric cancer with low tumor content. Our approach improves the CNV detection rate in targeted sequencing with tumor enrichment and the accuracy of CNV detection in archival samples without paired blood.

Tumor niche network-defined subtypes predict immunotherapy response of esophageal squamous cell cancer

Despite the promising outcomes of immune checkpoint inhibitors (ICIs), resistance to ICI presents a new challenge. Therefore, selecting patients for specific ICI applications is crucial for maximizing therapeutic efficacy. Herein, we curated 69 human esophageal squamous cell cancer (ESCC) patients' tumor microenvironment (TME) single-cell transcriptomic datasets to subtype ESCC. Integrative analyses of the cellular network and transcriptional signatures of T cells and myeloid cells define distinct ESCC subtypes characterized by T cell exhaustion, and interleukin (IL) and interferon (IFN) signaling. Furthermore, this approach classifies ESCC patients into ICI responders and non-responders, as validated by whole tumor transcriptomes and liquid biopsy-based single-cell transcriptomes of anti-PD-1 ICI responders and non-responders. Our study stratifies ESCC patients based on TME transcriptional network, providing novel insights into tumor niche remodeling and potentially predicting ICI responses in ESCC patients.

Stimuli-responsive incremental DNA machine auto-catalyzed CRISPR-Cas12a feedback amplification permits ultrasensitive molecular diagnosis of esophageal cancer-related microRNA

Development of new diagnostic methods is essential for disease diagnosis and treatment. In this work, we present a stimuli-responsive incremental DNA machine auto-catalyzed CRISPR-Cas12a (SRI-DNA machine/CRISPR-Cas12a) feedback amplification for ultrasensitive molecular detection of miRNA-21, which is an important biomarker related closely to the initiation and development of cancers, such as esophageal cancer. Strategically, the powerful SRI-DNA machine and efficient trans-cleavage activity of the CRISPR-Cas12a system are ingeniously integrated via a rationally designed probe termed as stem-elongated functional hairpin probe (SEF-HP). The SRI-DNA machine begins with the target miRNA, the trigger of the reaction, binding complementarily to the SEF-HP, followed by autonomously performed mechanical strand replication, cleavage, and displacement circuit at multiple sites. This conversion process led to the amplified generation of numerous DNA activators that are complementary with CRISPR RNA (CrRNA). Once formed the DNA activator/CrRNA heteroduplex, the trans-cleavage activity of the CRISPR-Cas12a was activated to nonspecific cleavage of single-stranded areas of a reporter probe for fluorescence emission. Under optimal conditions, the target miRNA can be detected with a wide linear range and an excellent specificity. As a proof-of-concept, this SRI-DNA machine/CRISPR-Cas12a feedback amplification system is adaptable and scalable to higher-order artificial amplification circuits for biomarkers detection, highlighting its promising potential in early diagnosis and disease treatment.

Cathepsin L promotes oesophageal squamous cell carcinoma development and may be associated with tumour-associated macrophages

Background

Oesophageal squamous cell carcinoma (ESCC) is a leading cause of cancer-related deaths worldwide because existing treatments are often insufficient. Therefore, new, reliable biomarkers must be identified. CTSL overexpression is closely associated with tumour progression and poor prognosis. However, the role and mechanism of CTSL as an oncogene in ESCC remain unclear.

Methods

Genome-wide association study (GWAS) data were used for Mendelian randomization analysis to investigate the possible relationships between CTSL and ESCC. The correlation between CTSL expression and prognosis was analysed using GEO, TCGA, and GEPIA data. We compared CTSL expression among the cell types using single-cell sequencing. Correlations between CTSL and the tumour microenvironment, immune cell infiltration, tumour mutational load, immunological checkpoints, and treatment sensitivity in patients with ESCC were investigated. Finally, using mouse models and cellular investigations, we assessed the effects of CTSL on the growth, apoptosis, and metastasis of ESCC tumour cells.

Results

CTSL was overexpressed in ESCC and correlated with prognosis. We also discovered its close association with cell immunity, especially with tumour-associated macrophages and immune checkpoints in the tumour microenvironment. CTSL may play a key role in ESCC development by affecting M2 macrophage polarisation. CTSL and the M2 macrophage marker genes showed significant positive correlations. Mendelian randomization analysis confirmed a relationship between CTSL and ESCC. Finally, our in vitro and in vivo experiments demonstrated that CTSL promoted the proliferation and migration of ESCC cells, validating our bioinformatic analysis.

Conclusion

CTSL emerged as a crucial gene in ESCC that influences patient prognosis and immunity, particularly in association with M2 macrophages. Therefore, targeting or modulating CTSL levels may provide new therapeutic strategies for patients with ESCC.

Transcriptome sequencing reveals novel biomarkers and immune cell infiltration in esophageal tumorigenesis

BACKGROUND

Esophageal squamous cell carcinoma (ESCC) is one of the most common malignancies worldwide, and its development comprises a multistep process from intraepithelial neoplasia (IN) to carcinoma (CA). However, the critical regulators and underlying molecular mechanisms remain largely unknown.

AIM

To explore the genes and infiltrating immune cells in the microenvironment that are associated with the multistage progression of ESCC to facilitate diagnosis and early intervention.

METHODS

A mouse model mimicking the multistage development of ESCC was established by providing warter containing 4-nitroquinoline 1-oxide (4NQO) to C57BL/6 mice. Moreover, we established a control group without 4NQO treatment of mice. Then, transcriptome sequencing was performed for esophageal tissues from patients with different pathological statuses, including low-grade IN (LGIN), high-grade IN (HGIN), and CA, and controlled normal tissue (NOR) samples. Differentially expressed genes (DEGs) were identified in the LGIN, HGIN, and CA groups, and the biological functions of the DEGs were analyzed via Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses. The CIBERSORT algorithm was used to detect the pattern of immune cell infiltration. Immunohistochemistry (IHC) was also conducted to validate our results. Finally, the Luminex multiplex cytokine analysis was utilized to measure the serum cytokine levels in the mice.

RESULTS

Compared with those in the NOR group, a total of 681541, and 840 DEGs were obtained in the LGIN, HGIN, and CA groups, respectively. Using the intersection of the three sets of DEGs, we identified 86 genes as key genes involved in the development of ESCC. Enrichment analysis revealed that these genes were enriched mainly in the keratinization, epidermal cell differentiation, and interleukin (IL)-17 signaling pathways. CIBERSORT analysis revealed that, compared with those in the NOR group, M0 and M1 macrophages in the 4NQO group showed stronger infiltration, which was validated by IHC. Serum cytokine analysis revealed that, compared with those in the NOR group, IL-1β and IL-6 were upregulated, while IL-10 was downregulated in the LGIN, HGIN, and CA groups. Moreover, the expression of the representative key genes, such as S100a8 and Krt6b, was verified in external human samples, and the results of immunohistochemical staining were consistent with the findings in mice.

CONCLUSION

We identified a set of key genes represented by S100a8 and Krt6b and investigated their potential biological functions. In addition, we found that macrophage infiltration and abnormal alterations in the levels of inflammation-associated cytokines, such as IL-1β, IL-6, and IL-10, in the peripheral blood may be closely associated with the development of ESCC.

Decoding spatiotemporal transcriptional dynamics and epithelial fibroblast crosstalk during gastroesophageal junction development through single cell analysis

The gastroesophageal squamocolumnar junction (GE-SCJ) is a critical tissue interface between the esophagus and stomach, with significant relevance in the pathophysiology of gastrointestinal diseases. Despite this, the molecular mechanisms underlying GE-SCJ development remain unclear. Using single-cell transcriptomics, organoids, and spatial analysis, we examine the cellular heterogeneity and spatiotemporal dynamics of GE-SCJ development from embryonic to adult mice. We identify distinct transcriptional states and signaling pathways in the epithelial and mesenchymal compartments of the esophagus and stomach during development. Fibroblast-epithelial interactions are mediated by various signaling pathways, including WNT, BMP, TGF-β, FGF, EGF, and PDGF. Our results suggest that fibroblasts predominantly send FGF and TGF-β signals to the epithelia, while epithelial cells mainly send PDGF and EGF signals to fibroblasts. We observe differences in the ligands and receptors involved in cell-cell communication between the esophagus and stomach. Our findings provide insights into the molecular mechanisms underlying GE-SCJ development and fibroblast-epithelial crosstalk involved, paving the way to elucidate mechanisms during adaptive metaplasia development and carcinogenesis.

Single-cell profiling of response to neoadjuvant chemo-immunotherapy in surgically resectable esophageal squamous cell carcinoma

BACKGROUND

The efficacy of neoadjuvant chemo-immunotherapy (NAT) in esophageal squamous cell carcinoma (ESCC) is challenged by the intricate interplay within the tumor microenvironment (TME). Unveiling the immune landscape of ESCC in the context of NAT could shed light on heterogeneity and optimize therapeutic strategies for patients.

METHODS

We analyzed single cells from 22 baseline and 24 post-NAT treatment samples of stage II/III ESCC patients to explore the association between the immune landscape and pathological response to neoadjuvant anti-PD-1 combination therapy, including pathological complete response (pCR), major pathological response (MPR), and incomplete pathological response (IPR).

RESULTS

Single-cell profiling identified 14 major cell subsets of cancer, immune, and stromal cells. Trajectory analysis unveiled an interesting link between cancer cell differentiation and pathological response to NAT. ESCC tumors enriched with less differentiated cancer cells exhibited a potentially favorable pathological response to NAT, while tumors enriched with clusters of more differentiated cancer cells may resist treatment. Deconvolution of transcriptomes in pre-treatment tumors identified gene signatures in response to NAT contributed by specific immune cell populations. Upregulated genes associated with better pathological responses in CD8 + effector T cells primarily involved interferon-gamma (IFNγ) signaling, neutrophil degranulation, and negative regulation of the T cell apoptotic process, whereas downregulated genes were dominated by those in the immune response-activating cell surface receptor signaling pathway. Natural killer cells in pre-treatment tumors from pCR patients showed a similar upregulation of gene expression in response to IFNγ but a downregulation of genes in the neutrophil-mediated immunity pathways. A decreased cellular contexture of regulatory T cells in ESCC TME indicated a potentially favorable pathological response to NAT. Cell-cell communication analysis revealed extensive interactions between CCL5 and its receptor CCR5 in various immune cells of baseline pCR tumors. Immune checkpoint interaction pairs, including CTLA4-CD86, TIGIT-PVR, LGALS9-HAVCR2, and TNFSF4-TNFRSF4, might serve as additional therapeutic targets for ICI therapy in ESCC.

CONCLUSIONS

This pioneering study unveiled an intriguing association between cancer cell differentiation and pathological response in esophageal cancer patients, revealing distinct subgroups of tumors for which neoadjuvant chemo-immunotherapy might be effective. We also delineated the immune landscape of ESCC tumors in the context of clinical response to NAT, which provides clinical insights for better understanding how patients respond to the treatment and further identifying novel therapeutic targets for ESCC patients in the future.

Recent advances in differential expression analysis for single-cell RNA-seq and spatially resolved transcriptomic studies

Differential expression (DE) analysis is a necessary step in the analysis of single-cell RNA sequencing (scRNA-seq) and spatially resolved transcriptomics (SRT) data. Unlike traditional bulk RNA-seq, DE analysis for scRNA-seq or SRT data has unique characteristics that may contribute to the difficulty of detecting DE genes. However, the plethora of DE tools that work with various assumptions makes it difficult to choose an appropriate one. Furthermore, a comprehensive review on detecting DE genes for scRNA-seq data or SRT data from multi-condition, multi-sample experimental designs is lacking. To bridge such a gap, here, we first focus on the challenges of DE detection, then highlight potential opportunities that facilitate further progress in scRNA-seq or SRT analysis, and finally provide insights and guidance in selecting appropriate DE tools or developing new computational DE methods.

Reciprocal inhibition between TP63 and STAT1 regulates anti-tumor immune response through interferon-γ signaling in squamous cancer

Squamous cell carcinomas (SCCs) are common and aggressive malignancies. Immune check point blockade (ICB) therapy using PD-1/PD-L1 antibodies has been approved in several types of advanced SCCs. However, low response rate and treatment resistance are common. Improving the efficacy of ICB therapy requires better understanding of the mechanism of immune evasion. Here, we identify that the SCC-master transcription factor TP63 suppresses interferon-γ (IFNγ) signaling. TP63 inhibition leads to increased CD8+ T cell infiltration and heighten tumor killing in in vivo syngeneic mouse model and ex vivo co-culture system, respectively. Moreover, expression of TP63 is negatively correlated with CD8+ T cell infiltration and activation in patients with SCC. Silencing of TP63 enhances the anti-tumor efficacy of PD-1 blockade by promoting CD8+ T cell infiltration and functionality. Mechanistically, TP63 and STAT1 mutually suppress each other to regulate the IFNγ signaling by co-occupying and co-regulating their own promoters and enhancers. Together, our findings elucidate a tumor-extrinsic function of TP63 in promoting immune evasion of SCC cells. Over-expression of TP63 may serve as a biomarker predicting the outcome of SCC patients treated with ICB therapy, and targeting TP63/STAT/IFNγ axis may enhance the efficacy of ICB therapy for this deadly cancer.

Autophagy Contributes to Homeostasis in Esophageal Epithelium Where High Autophagic Vesicle Level Marks Basal Cells With Limited Proliferation and Enhanced Self-Renewal Potential

BACKGROUND & AIMS

Autophagy plays roles in esophageal pathologies both benign and malignant. Here, we aim to define the role of autophagy in esophageal epithelial homeostasis.

METHODS

We generated tamoxifen-inducible, squamous epithelial-specific Atg7 (autophagy related 7) conditional knockout mice to evaluate effects on esophageal homeostasis and response to the carcinogen 4-nitroquinoline 1-oxide (4NQO) using histologic and biochemical analyses. We fluorescence-activated cell sorted esophageal basal cells based on fluorescence of the autophagic vesicle (AV)-identifying dye Cyto-ID and then subjected these cells to transmission electron microscopy, image flow cytometry, three-dimensional organoid assays, RNA sequencing, and cell cycle analysis. Three-dimensional organoids were subjected to passaging, single-cell RNA sequencing, cell cycle analysis, and immunostaining.

RESULTS

Genetic autophagy inhibition in squamous epithelium resulted in increased proliferation of esophageal basal cells under homeostatic conditions and also was associated with significant weight loss in mice treated with 4NQO that further displayed perturbed epithelial tissue architecture. Esophageal basal cells with high AV level (Cyto-IDHigh) displayed limited organoid formation capability on initial plating but passaged more efficiently than their counterparts with low AV level (Cyto-IDLow). RNA sequencing suggested increased autophagy in Cyto-IDHigh esophageal basal cells along with decreased cell cycle progression, the latter of which was confirmed by cell cycle analysis. Single-cell RNA sequencing of three-dimensional organoids generated by Cyto-IDLow and Cyto-IDHigh cells identified expansion of 3 cell populations and enrichment of G2/M-associated genes in the Cyto-IDHigh group. Ki67 expression was also increased in organoids generated by Cyto-IDHigh cells, including in basal cells localized beyond the outermost cell layer.

CONCLUSIONS

Autophagy contributes to maintenance of the esophageal proliferation-differentiation gradient. Esophageal basal cells with high AV level exhibit limited proliferation and generate three-dimensional organoids with enhanced self-renewal capacity.

Advanced insights on tumor-associated macrophages revealed by single-cell RNA sequencing: The intratumor heterogeneity, functional phenotypes, and cellular interactions

Single-cell RNA sequencing (scRNA-seq) is an emerging technology used for cellular transcriptome analysis. The application of scRNA-seq has led to profoundly advanced oncology research, continuously optimizing novel therapeutic strategies. Intratumor heterogeneity extensively consists of all tumor components, contributing to different tumor behaviors and treatment responses. Tumor-associated macrophages (TAMs), the core immune cells linking innate and adaptive immunity, play significant roles in tumor progression and resistance to therapies. Moreover, dynamic changes occur in TAM phenotypes and functions subject to the regulation of the tumor microenvironment. The heterogeneity of TAMs corresponding to the state of the tumor microenvironment has been comprehensively recognized using scRNA-seq. Herein, we reviewed recent research and summarized variations in TAM phenotypes and functions from a developmental perspective to better understand the significance of TAMs in the tumor microenvironment.

scRNA-seq revealed high stemness epithelial malignant cell clusters and prognostic models of lung adenocarcinoma

Lung adenocarcinoma (LUAD) is one of the sole causes of death in lung cancer patients. This study combined with single-cell RNA-seq analysis to identify tumor stem-related prognostic models to predict the prognosis of lung adenocarcinoma, chemotherapy agents, and immunotherapy efficacy. mRNA expression-based stemness index (mRNAsi) was determined by One Class Linear Regression (OCLR). Differentially expressed genes (DEGs) were detected by limma package. Single-cell RNA-seq analysis in GSE123902 dataset was performed using Seurat package. Weighted Co-Expression Network Analysis (WGCNA) was built by rms package. Cell differentiation ability was determined by CytoTRACE. Cell communication analysis was performed by CellCall and CellChat package. Prognosis model was constructed by 10 machine learning and 101 combinations. Drug predictive analysis was conducted by pRRophetic package. Immune microenvironment landscape was determined by ESTIMATE, MCP-Counter, ssGSEA analysis. Tumor samples have higher mRNAsi, and the high mRNAsi group presents a worse prognosis. Turquoise module was highly correlated with mRNAsi in TCGA-LUAD dataset. scRNA analysis showed that 22 epithelial cell clusters were obtained, and higher CSCs malignant epithelial cells have more complex cellular communication with other cells and presented dedifferentiation phenomenon. Cellular senescence and Hippo signaling pathway are the major difference pathways between high- and low CSCs malignant epithelial cells. The pseudo-temporal analysis shows that cluster1, 2, high CSC epithelial cells, are concentrated at the end of the differentiation trajectory. Finally, 13 genes were obtained by intersecting genes in turquoise module, Top200 genes in hdWGCNA, DEGs in high- and low- mRNAsi group as well as DEGs in tumor samples vs. normal group. Among 101 prognostic models, average c-index (0.71) was highest in CoxBoost + RSF model. The high-risk group samples had immunosuppressive status, higher tumor malignancy and low benefit from immunotherapy. This work found that malignant tumors and malignant epithelial cells have high CSC characteristics, and identified a model that could predict the prognosis, immune microenvironment, and immunotherapy of LUAD, based on CSC-related genes. These results provided reference value for the clinical diagnosis and treatment of LUAD.

Aberrant epithelial cell interaction promotes esophageal squamous-cell carcinoma development and progression

Epithelial-mesenchymal transition (EMT) and proliferation play important roles in epithelial cancer formation and progression, but what molecules and how they trigger EMT is largely unknown. Here we performed spatial transcriptomic and functional analyses on samples of multistage esophageal squamous-cell carcinoma (ESCC) from mice and humans to decipher these critical issues. By investigating spatiotemporal gene expression patterns and cell-cell interactions, we demonstrated that the aberrant epithelial cell interaction via EFNB1-EPHB4 triggers EMT and cell cycle mediated by downstream SRC/ERK/AKT signaling. The aberrant epithelial cell interaction occurs within the basal layer at early precancerous lesions, which expands to the whole epithelial layer and strengthens along the cancer development and progression. Functional analysis revealed that the aberrant EFNB1-EPHB4 interaction is caused by overexpressed ΔNP63 due to TP53 mutation, the culprit in human ESCC tumorigenesis. Our results shed new light on the role of TP53-TP63/ΔNP63-EFNB1-EPHB4 axis in EMT and cell proliferation in epithelial cancer formation.

Deciphering tumor-infiltrating dendritic cells in the single-cell era

Dendritic cells (DCs) serve as a pivotal link connecting innate and adaptive immunity by processing tumor-derived antigens and activating T cells. The advent of single-cell sequencing has revolutionized the categorization of DCs, enabling a high-resolution characterization of the previously unrecognized diversity of DC populations infiltrating the intricate tumor microenvironment (TME). The application of single-cell sequencing technologies has effectively elucidated the heterogeneity of DCs present in the tumor milieu, yielding invaluable insights into their subpopulation structures and functional diversity. This review provides a comprehensive summary of the current state of knowledge regarding DC subtypes in the TME, drawing from single-cell studies conducted across various human tumors. We focused on the categorization, functions, and interactions of distinct DC subsets, emphasizing their crucial roles in orchestrating tumor-related immune responses. Additionally, we delve into the potential implications of these findings for the identification of predictive biomarkers and therapeutic targets. Enhanced insight into the intricate interplay between DCs and the TME promises to advance our comprehension of tumor immunity and, in turn, pave the way for the development of more efficacious cancer immunotherapies.

Collagen 1-mediated CXCL1 secretion in tumor cells activates fibroblasts to promote radioresistance of esophageal cancer

Esophageal squamous-cell carcinoma (ESCC) is commonly treated with radiotherapy; however, radioresistance hinders its clinical effectiveness, and the underlying mechanism remains elusive. Here, we develop patient-derived xenografts (PDXs) from 19 patients with ESCC to investigate the mechanisms driving radioresistance. Using RNA sequencing, cytokine arrays, and single-cell RNA sequencing, we reveal an enrichment of cancer-associated fibroblast (CAF)-derived collagen type 1 (Col1) and tumor-cell-derived CXCL1 in non-responsive PDXs. Col1 not only promotes radioresistance by augmenting DNA repair capacity but also induces CXCL1 secretion in tumor cells. Additionally, CXCL1 further activates CAFs via the CXCR2-STAT3 pathway, establishing a positive feedback loop. Directly interfering with tumor-cell-derived CXCL1 or inhibiting the CXCL1-CXCR2 pathway effectively restores the radiosensitivity of radioresistant xenografts in vivo. Collectively, our study provides a comprehensive understanding of the molecular mechanisms underlying radioresistance and identifies potential targets to improve the efficacy of radiotherapy for ESCC.

Mammalian esophageal stratified tissue homeostasis is maintained distinctively by the epithelial pluripotent p63+Sox2+ and p63-Sox2+ cell populations

Self-renewing, damage-repair and differentiation of mammalian stratified squamous epithelia are subject to tissue homeostasis, but the regulation mechanisms remain elusive. Here, we investigate the esophageal squamous epithelial tissue homeostasis in vitro and in vivo. We establish a rat esophageal organoid (rEO) in vitro system and show that the landscapes of rEO formation, development and maturation trajectories can mimic those of rat esophageal epithelia in vivo. Single-cell RNA sequencing (scRNA-seq), snapshot immunostaining and functional analyses of stratified "matured" rEOs define that the epithelial pluripotent stem cell determinants, p63 and Sox2, play crucial but distinctive roles for regulating mammalian esophageal tissue homeostasis. We identify two cell populations, p63+Sox2+ and p63-Sox2+, of which the p63+Sox2+ population presented at the basal layer is the cells of origin required for esophageal epithelial stemness maintenance and proliferation, whereas the p63-Sox2+ population presented at the suprabasal layers is the cells of origin having a dual role for esophageal epithelial differentiation (differentiation-prone fate) and rapid tissue damage-repair responses (proliferation-prone fate). Given the fact that p63 and Sox2 are developmental lineage oncogenes and commonly overexpressed in ESCC tissues, p63-Sox2+ population could not be detected in organoids formed by esophageal squamous cell carcinoma (ESCC) cell lines. Taken together, these findings reveal that the tissue homeostasis is maintained distinctively by p63 and/or Sox2-dependent cell lineage populations required for the tissue renewing, damage-repair and protection of carcinogenesis in mammalian esophagi.

Autophagy contributes to homeostasis in esophageal epithelium where high autophagic vesicle content marks basal cells with limited proliferation and enhanced self-renewal potential

Background & Aims

Autophagy has been demonstrated to play roles in esophageal pathologies both benign and malignant. Here, we aim to define the role of autophagy in esophageal epithelium under homeostatic conditions.

Methods

We generated tamoxifen-inducible, squamous epithelial-specific Atg7 (autophagy related 7) conditional knockout mice to evaluate effects on esophageal homeostasis and response to the carcinogen 4-nitroquinoline 1-oxide (4NQO) using histological and biochemical analyses. We FACS sorted esophageal basal cells based upon fluorescence of the autophagic vesicle (AV)-identifying dye Cyto-ID, then subjected these cells to transmission electron microscopy, image flow cytometry, 3D organoid assays, RNA-Sequencing (RNA-Seq), and cell cycle analysis. 3D organoids were subjected to passaging, single cell (sc) RNA-Seq, cell cycle analysis, and immunostaining.

Results

Genetic autophagy inhibition in squamous epithelium resulted in increased proliferation of esophageal basal cells. Esophageal basal cells with high AV level (Cyto-ID High ) displayed limited organoid formation capability upon initial plating but passaged more efficiently than their counterparts with low AV level (Cyto-ID Low ). RNA-Seq suggested increased autophagy in Cyto- ID High esophageal basal cells along with decreased cell cycle progression, the latter of which was confirmed by cell cycle analysis. scRNA-Seq of 3D organoids generated by Cyto-ID Low and Cyto- ID High cells identified expansion of 3 cell populations, enrichment of G2/M-associated genes, and aberrant localization of cell cycle-associated genes beyond basal cell populations in the Cyto- ID High group. Ki67 expression was also increased in organoids generated by Cyto-ID High cells, including in cells beyond the basal cell layer. Squamous epithelial-specific autophagy inhibition induced significant weight loss in mice treated with 4NQO that further displayed perturbed epithelial tissue architecture.

Conclusions

High AV level identifies esophageal epithelium with limited proliferation and enhanced self-renewal capacity that contributes to maintenance of the esophageal proliferation- differentiation gradient in vivo .

Single-cell transcriptomic analysis uncovers the origin and intratumoral heterogeneity of parotid pleomorphic adenoma

Pleomorphic adenoma (PA) is the most common benign tumour in the salivary gland and has high morphological complexity. However, the origin and intratumoral heterogeneity of PA are largely unknown. Here, we constructed a comprehensive atlas of PA at single-cell resolution and showed that PA exhibited five tumour subpopulations, three recapitulating the epithelial states of the normal parotid gland, and two PA-specific epithelial cell (PASE) populations unique to tumours. Then, six subgroups of PASE cells were identified, which varied in epithelium, bone, immune, metabolism, stemness and cell cycle signatures. Moreover, we revealed that CD36+ myoepithelial cells were the tumour-initiating cells (TICs) in PA, and were dominated by the PI3K-AKT pathway. Targeting the PI3K-AKT pathway significantly inhibited CD36+ myoepithelial cell-derived tumour spheres and the growth of PA organoids. Our results provide new insights into the diversity and origin of PA, offering an important clinical implication for targeting the PI3K-AKT signalling pathway in PA treatment.

Key Genetic Determinants Driving Esophageal Squamous Cell Carcinoma Initiation and Immune Evasion

BACKGROUND & AIMS

Despite recent progress in identifying aberrant genetic and epigenetic alterations in esophageal squamous cell carcinoma (ESCC), the mechanism of ESCC initiation remains unknown.

METHODS

Using CRISPR/Cas 9-based genetic ablation, we targeted 9 genes (TP53, CDKN2A, NOTCH1, NOTCH3, KMT2D, KMT2C, FAT1, FAT4, and AJUBA) in murine esophageal organoids. Transcriptomic phenotypes of organoids and chemokine released by organoids were analyzed by single-cell RNA sequencing. Tumorigenicity and immune evasion of organoids were monitored by allograft transplantation. Human ESCC single-cell RNA sequencing data sets were analyzed to classify patients and find subsets relevant to organoid models and immune evasion.

RESULTS

We established 32 genetically engineered esophageal organoids and identified key genetic determinants that drive ESCC initiation. A single-cell transcriptomic analysis uncovered that Trp53, Cdkn2a, and Notch1 (PCN) triple-knockout induces neoplastic features of ESCC by generating cell lineage heterogeneity and high cell plasticity. PCN knockout also generates an immunosuppressive niche enriched with exhausted T cells and M2 macrophages via the CCL2-CCR2 axis. Mechanistically, CDKN2A inactivation transactivates CCL2 via nuclear factor-κB. Moreover, comparative single-cell transcriptomic analyses stratified patients with ESCC and identified a specific subtype recapitulating the PCN-type ESCC signatures, including the high expression of CCL2 and CD274/PD-L1.

CONCLUSIONS

Our study unveils that loss of TP53, CDKN2A, and NOTCH1 induces esophageal neoplasia and immune evasion for ESCC initiation and proposes the CCL2 blockade as a viable option for targeting PCN-type ESCC.

Eosinophils Exert Antitumorigenic Effects in the Development of Esophageal Squamous Cell Carcinoma

BACKGROUND AND AIMS

Eosinophils are present in several solid tumors and have context-dependent function. Our aim is to define the contribution of eosinophils in esophageal squamous cell carcinoma (ESCC), as their role in ESCC is unknown.

METHODS

Eosinophils were enumerated in tissues from 2 ESCC cohorts. Mice were treated with 4-NQO for 8 weeks to induce precancer or 16 weeks to induce carcinoma. The eosinophil number was modified by a monoclonal antibody to interleukin-5 (IL5mAb), recombinant IL-5 (rIL-5), or genetically with eosinophil-deficient (ΔdblGATA) mice or mice deficient in eosinophil chemoattractant eotaxin-1 (Ccl11-/-). Esophageal tissue and eosinophil-specific RNA sequencing was performed to understand eosinophil function. Three-dimensional coculturing of eosinophils with precancer or cancer cells was done to ascertain direct effects of eosinophils.

RESULTS

Activated eosinophils are present in higher numbers in early-stage vs late-stage ESCC. Mice treated with 4-NQO exhibit more esophageal eosinophils in precancer vs cancer. Correspondingly, epithelial cell Ccl11 expression is higher in mice with precancer. Eosinophil depletion using 3 mouse models (Ccl11-/- mice, ΔdblGATA mice, IL5mAb treatment) all display exacerbated 4-NQO tumorigenesis. Conversely, treatment with rIL-5 increases esophageal eosinophilia and protects against precancer and carcinoma. Tissue and eosinophil RNA sequencing revealed eosinophils drive oxidative stress in precancer. In vitro coculturing of eosinophils with precancer or cancer cells resulted in increased apoptosis in the presence of a degranulating agent, which is reversed with NAC, a reactive oxygen species scavenger. ΔdblGATA mice exhibited increased CD4 T cell infiltration, IL-17, and enrichment of IL-17 protumorigenic pathways.

CONCLUSION

Eosinophils likely protect against ESCC through reactive oxygen species release during degranulation and suppression of IL-17.

Spatial transcriptomics analysis of esophageal squamous precancerous lesions and their progression to esophageal cancer

Esophageal squamous precancerous lesions (ESPL) are the precursors of esophageal squamous cell carcinoma (ESCC) including low-grade and high-grade intraepithelial neoplasia. Due to the absence of molecular indicators, which ESPL will eventually develop into ESCC and thus should be treated is not well defined. Indicators, for predicting risks of ESCC at ESPL stages, are an urgent need. We perform spatial whole-transcriptome atlas analysis, which can eliminate other tissue interference by sequencing the specific ESPL regions. In this study, the expression of TAGLN2 significantly increases, while CRNN expression level decreases along the progression of ESCC. Additionally, TAGLN2 protein level significantly increases in paired after-progression tissues compared with before-progression samples, while CRNN expression decreases. Functional studies suggest that TAGLN2 promotes ESCC progression, while CRNN inhibits it by regulating cell proliferation. Taken together, TAGLN2 and CRNN are suggested as candidate indicators for the risk of ESCC at ESPL stages.

IL-1β transgenic mouse model of inflammation driven esophageal and oral squamous cell carcinoma

Chronic inflammation is integral to the development of esophageal adenocarcinoma (EAC) and esophageal squamous cell carcinoma (ESCC), although the latter has not been associated with reflux esophagitis. The L2-IL-1β transgenic mice, expressing human interleukin (IL)-1β in the oral, esophageal and forestomach squamous epithelia feature chronic inflammation and a stepwise development of Barrett's esophagus-like metaplasia, dysplasia and adenocarcinoma at the squamo-columnar junction. However, the functional consequences of IL-1β-mediated chronic inflammation in the oral and esophageal squamous epithelia remain elusive. We report for the first time that in addition to the previously described Barrett's esophagus-like metaplasia, the L2-IL-1β mice also develop squamous epithelial dysplasia with progression to squamous cell carcinoma (SCC) in the esophagus and the tongue. L2-IL-1β showed age-dependent progression of squamous dysplasia to SCC with approximately 40% (n = 49) and 23.5% (n = 17) incidence rates for esophageal and tongue invasive SCC respectively, by 12-15 months of age. Interestingly, SCC development and progression in L2-IL-1β was similar in both Germ Free (GF) and Specific Pathogen Free (SPF) conditions. Immunohistochemistry revealed a T cell predominant inflammatory profile with enhanced expression of Ki67, Sox2 and the DNA double-strand break marker, γ-H2AX, in the dysplastic squamous epithelia of L2-IL-1β mice. Pro-inflammatory cytokines, immunomodulatory players, chemoattractants for inflammatory cells (T cells, neutrophils, eosinophils, and macrophages) and oxidative damage marker, iNOS, were significantly increased in the esophageal and tongue tissues of L2-IL-1β mice. Our recent findings have expanded the translational utility of the IL-1β mouse model to aid in further characterization of the key pathways of inflammation driven BE and EAC as well as ESCC and Oral SCC.

Assessing the Link between Diabetic Metabolic Dysregulation and Breast Cancer Progression

Diabetes mellitus is a burdensome disease that affects various cellular functions through altered glucose metabolism. Several reports have linked diabetes to cancer development; however, the exact molecular mechanism of how diabetes-related traits contribute to cancer progression is not fully understood. The current study aimed to explore the molecular mechanism underlying the potential effect of hyperglycemia combined with hyperinsulinemia on the progression of breast cancer cells. To this end, gene dysregulation induced by the exposure of MCF7 breast cancer cells to hyperglycemia (HG), or a combination of hyperglycemia and hyperinsulinemia (HGI), was analyzed using a microarray gene expression assay. Hyperglycemia combined with hyperinsulinemia induced differential expression of 45 genes (greater than or equal to two-fold), which were not shared by other treatments. On the other hand, in silico analysis performed using a publicly available dataset (GEO: GSE150586) revealed differential upregulation of 15 genes in the breast tumor tissues of diabetic patients with breast cancer when compared with breast cancer patients with no diabetes. SLC26A11, ALDH1A3, MED20, PABPC4 and SCP2 were among the top upregulated genes in both microarray data and the in silico analysis. In conclusion, hyperglycemia combined with hyperinsulinemia caused a likely unique signature that contributes to acquiring more carcinogenic traits. Indeed, these findings might potentially add emphasis on how monitoring diabetes-related metabolic alteration as an adjunct to diabetes therapy is important in improving breast cancer outcomes. However, further detailed studies are required to decipher the role of the highlighted genes, in this study, in the pathogenesis of breast cancer in patients with a different glycemic index.

IL1R1+ cancer-associated fibroblasts drive tumor development and immunosuppression in colorectal cancer

Fibroblasts have a considerable functional and molecular heterogeneity and can play various roles in the tumor microenvironment. Here we identify a pro-tumorigenic IL1R1⁺, IL-1-high-signaling subtype of fibroblasts, using multiple colorectal cancer (CRC) patient single cell sequencing datasets. This subtype of fibroblasts is linked to T cell and macrophage suppression and leads to increased cancer cell growth in 3D co-culture assays. Furthermore, both a fibroblast-specific IL1R1 knockout and IL-1 receptor antagonist Anakinra administration reduce tumor growth in vivo. This is accompanied by reduced intratumoral Th17 cell infiltration. Accordingly, CRC patients who present with IL1R1-expressing cancer-associated-fibroblasts (CAFs), also display elevated levels of immune exhaustion markers, as well as an increased Th17 score and an overall worse survival. Altogether, this study underlines the therapeutic value of targeting IL1R1-expressing CAFs in the context of CRC.

Eosinophils exert direct and indirect anti-tumorigenic effects in the development of esophageal squamous cell carcinoma

Background/Aims

Eosinophils are present in several solid tumors and have context-dependent function. Our aim is to define the contribution of eosinophils in esophageal squamous cell carcinoma (ESCC), since their role in ESCC is unknown.

Methods

Eosinophils were enumerated in tissues from two ESCC cohorts. Mice were treated with 4-nitroquinolone-1-oxide (4-NQO) for 8 weeks to induce pre-cancer or 16 weeks to induce carcinoma. Eosinophil number was modified by monoclonal antibody to IL-5 (IL5mAb), recombinant IL-5 (rIL-5), or genetically with eosinophil-deficient (ΔdblGATA) mice or mice deficient in eosinophil chemoattractant eotaxin-1 ( Ccl11 -/- ). Esophageal tissue and eosinophil specific RNA-sequencing was performed to understand eosinophil function. 3-D co-culturing of eosinophils with pre-cancer or cancer cells was done to ascertain direct effects of eosinophils.

Results

Activated eosinophils are present in higher numbers in early stage versus late stage ESCC. Mice treated with 4-NQO exhibit more esophageal eosinophils in pre-cancer versus cancer. Correspondingly, epithelial cell Ccl11 expression is higher in mice with pre-cancer. Eosinophil depletion using three mouse models ( Ccl11 -/- mice, ΔdblGATA mice, IL5mAb treatment) all display exacerbated 4-NQO tumorigenesis. Conversely, treatment with rIL-5 increases esophageal eosinophilia and protects against pre-cancer and carcinoma. Tissue and eosinophil RNA-sequencing revealed eosinophils drive oxidative stress in pre-cancer. In vitro co-culturing of eosinophils with pre-cancer or cancer cells resulted in increased apoptosis in the presence of a degranulating agent, which is reversed with N-acetylcysteine, a reactive oxygen species (ROS) scavenger. ΔdblGATA mice exhibited increased CD4 T cell infiltration, IL-17, and enrichment of IL-17 pro-tumorigenic pathways.

Conclusion

Eosinophils likely protect against ESCC through ROS release during degranulation and suppression of IL-17.

scRNASequest: an ecosystem of scRNA-seq analysis, visualization, and publishing

BACKGROUND

Single-cell RNA sequencing is a state-of-the-art technology to understand gene expression in complex tissues. With the growing amount of data being generated, the standardization and automation of data analysis are critical to generating hypotheses and discovering biological insights.

RESULTS

Here, we present scRNASequest, a semi-automated single-cell RNA-seq (scRNA-seq) data analysis workflow which allows (1) preprocessing from raw UMI count data, (2) harmonization by one or multiple methods, (3) reference-dataset-based cell type label transfer and embedding projection, (4) multi-sample, multi-condition single-cell level differential gene expression analysis, and (5) seamless integration with cellxgene VIP for visualization and with CellDepot for data hosting and sharing by generating compatible h5ad files.

CONCLUSIONS

We developed scRNASequest, an end-to-end pipeline for single-cell RNA-seq data analysis, visualization, and publishing. The source code under MIT open-source license is provided at https://github.com/interactivereport/scRNASequest . We also prepared a bookdown tutorial for the installation and detailed usage of the pipeline: https://interactivereport.github.io/scRNAsequest/tutorial/docs/ . Users have the option to run it on a local computer with a Linux/Unix system including MacOS, or interact with SGE/Slurm schedulers on high-performance computing (HPC) clusters.

A method for scoring 4-nitroquinoline 1-oxide-induced murine esophageal squamous neoplasia

A mouse model for esophageal squamous cell carcinoma (ESCC) is induced by oral administration of the carcinogen 4-nitroquinoline 1-oxide (4-NQO). There is not an objective method for determining histopathologic severity of disease in this model. We aim to create a clearly defined and easily applied scoring system that can quantify the severity of 4-NQO-induced murine ESCC. Fifteen wild-type C57BL/6J mice were treated with 4-NQO for 8 (n = 8) or 16 (n = 7) weeks, while the rest (n = 9) were treated with vehicle, as 8 weeks of 4-NQO typically results in dysplasia and 16 weeks in carcinoma. We identified histologic abnormalities of the esophagus in this model and developed metrics to grade severity of dysplasia, papillomas, and invasion. Scores were then calculated using quantitative digitized image analysis for measuring depth and extent of each feature within the entire sample. Each feature was also assigned a weight based on its relation to cancer severity. Histology scores were significantly different in the three groups, suggesting that this method can discriminate dysplasia from carcinoma. This model can be applied to any mouse treated with 4-NQO.

Toosendanin targeting eEF2 impedes Topoisomerase I & II protein translation to suppress esophageal squamous cell carcinoma growth

BACKGROUND

Although molecular targets such as HER2, TP53 and PIK3CA have been widely studied in esophageal cancer, few of them were successfully applied for clinical treatment. Therefore, it is urgent to discover novel actionable targets and inhibitors. Eukaryotic translational elongation factor 2 (eEF2) is reported to be highly expressed in various cancers. However, its contribution to the maintenance and progression of cancer has not been fully clarified.

METHODS

In the present study, we utilized tissue array to evaluate eEF2 protein expression and clinical significance in esophageal squamous cell carcinoma (ESCC). Next, we performed knockdown, overexpression, RNA-binding protein immunoprecipitation (RIP) sequence, and nascent protein synthesis assays to explore the molecular function of eEF2. Furthermore, we utilized compound screening, Surface Plasmon Resonance (SPR), Isothermal Titration Calorimetry (ITC) assay, cell proliferation and Patient derived xenograft (PDX) mouse model assays to discover an eEF2 inhibitor and assess its effects on ESCC growth.

RESULTS

We found that eEF2 were highly expressed in ESCC and negatively associated with the prognosis of ESCC patients. Knocking down of eEF2 suppressed the cell proliferation and colony formation of ESCC. eEF2 bond with the mRNA of Topoisomerase II (TOP1) and Topoisomerase II (TOP2) and enhanced the protein biosynthesis of TOP1 and TOP2. We also identified Toosendanin was a novel inhibitor of eEF2 and Toosendanin inhibited the growth of ESCC in vitro and in vivo.

CONCLUSIONS

Our findings show that Toosendanin treatment suppresses ESCC growth through targeting eEF2 and regulating downstream TOP1 and TOP2 biosynthesis. eEF2 could be supplied as a potential therapeutic target in the further clinical studies.

CD8+ T cell exhaustion and cancer immunotherapy

Immunotherapy plays an increasingly important role in the treatment of most malignant tumors, and CD8⁺ T cells are the most important antitumor effector cells in the process of immunotherapy, and their number and functional status largely determine the antitumor effect. However, under continuous antigen exposure and the stimulation of inflammatory factors, CD8⁺ T cells gradually show a weakened proliferation and effector function, accompanied by the expression of a variety of inhibitory receptors. This state is known as CD8⁺ T cell "exhaustion" and often leads to the loss of control and progression of tumors. Recent studies provided us a better understanding of the mechanisms of T cell exhaustion, this review provides an overview of the activation, exhaustion mechanisms and exhaustion characteristics of CD8⁺ T cells. Although immunotherapy can reverse the exhaustion of CD8⁺ T cells and significantly improve the antitumor effects, single immunotherapy often has limitations, and it is difficult to achieve satisfactory antitumor effects, therefore, this review also summarizes up-to-date information related to cancer immunotherapy, and these emerging insights provide promising clues to the future management of malignant tumors.

Cellular states are coupled to genomic and viral heterogeneity in HPV-related oropharyngeal carcinoma

Head and neck squamous cell carcinoma (HNSCC) includes a subset of cancers driven by human papillomavirus (HPV). Here we use single-cell RNA-seq to profile both HPV-positive and HPV-negative oropharyngeal tumors, uncovering a high level of cellular diversity within and between tumors. First, we detect diverse chromosomal aberrations within individual tumors, suggesting genomic instability and enabling the identification of malignant cells even at pathologically negative margins. Second, we uncover diversity with respect to HNSCC subtypes and other cellular states such as the cell cycle, senescence and epithelial-mesenchymal transitions. Third, we find heterogeneity in viral gene expression within HPV-positive tumors. HPV expression is lost or repressed in a subset of cells, which are associated with a decrease in HPV-associated cell cycle phenotypes, decreased response to treatment, increased invasion and poor prognosis. These findings suggest that HPV expression diversity must be considered during diagnosis and treatment of HPV-positive tumors, with important prognostic ramifications.

Epithelial cells activate fibroblasts to promote esophageal cancer development

Esophageal squamous-cell carcinoma (ESCC) develops through multistage epithelial cancer formation, i.e., from normal epithelium, low- and high-grade intraepithelial neoplasia to invasive carcinoma. However, how the precancerous lesions progress to carcinoma remains elusive. Here, we report a comprehensive single-cell RNA sequencing and spatial transcriptomic study of 79 multistage esophageal lesions from 29 patients with ESCC. We reveal a gradual and significant loss of ANXA1 expression in epithelial cells due to its transcription factor KLF4 suppression along the lesion progression. We demonstrate that ANXA1 is a ligand to formyl peptide receptor type 2 (FPR2) on fibroblasts that maintain fibroblast homeostasis. Loss of ANXA1 leads to uncontrolled transformation of normal fibroblasts into cancer-associated fibroblasts (CAFs), which can be enhanced by secreted TGF-β from malignant epithelial cells. Given the role of CAFs in cancer, our study underscores ANXA1/FPR2 signaling as an important crosstalk mechanism between epithelial cells and fibroblasts in promoting ESCC.

Mitochondrial DNA Mutations as Natural Barcodes for Lineage Tracing of Murine Tumor Models

Murine models are indispensable tools for functional genomic studies and preclinical testing of novel therapeutic approaches. Mitochondrial single-cell assay for transposase-accessible chromatin using sequencing (mtscATAC-seq) enables the dissection of cellular heterogeneity and clonal dynamics by capturing chromatin accessibility, copy-number variations (CNV), and mitochondrial DNA (mtDNA) mutations, yet its applicability to murine studies remains unexplored. By leveraging mtscATAC-seq in novel chronic lymphocytic leukemia and Richter syndrome mouse models, we report the detection of mtDNA mutations, particularly in highly proliferative murine cells, alongside CNV and chromatin state changes indicative of clonal evolution upon secondary transplant. This study thus demonstrates the feasibility and utility of multi-modal single-cell and natural barcoding approaches to characterize murine cancer models.

SIGNIFICANCE

mtDNA mutations can serve as natural barcodes to enable lineage tracing in murine cancer models, which can be used to provide new insights into disease biology and to identify therapeutic vulnerabilities.

Single-cell transcriptomic analysis deciphers key transitional signatures associated with oncogenic evolution in human intramucosal oesophageal squamous cell carcinoma

BACKGROUND AND AIMS

The early diagnosis and intervention of oesophageal squamous cell carcinoma (ESCC) are particularly important because of the lack of effective therapies and poor prognosis. Comprehensive research on early ESCC at the single-cell level is rare due to the need for fresh and high-quality specimens obtained from ESD. This study aims to systematically describe the cellular atlas of human intramucosal ESCC.

METHODS

Five paired samples of intramucosal ESCC, para-ESCC oesophageal tissues from endoscopically resected specimens and peripheral blood mononuclear cells were adopted for scRNA-seq analysis. Computational pipeline scMetabolism was applied to quantify the metabolic diversity of single cells.

RESULTS

A total of 164 715 cells were profiled. Epithelial cells exhibited high intra-tumoural heterogeneity and two evolutionary trajectories during ESCC tumorigenesis initiated from proliferative cells, and then through an intermediate state, to two different terminal states of normally differentiated epithelial cells or malignant cells, respectively. The abundance of CD8+ TEX s, Tregs and PD1+ CD4+ T cells suggested an exhausted and suppressive immune microenvironment. Several genes in immune cells, such as CXCL13, CXCR5 and PADI4, were identified as new biomarkers for poor prognosis. A new subcluster of malignant cells associated with metastasis and angiogenesis that appeared at an early stage compared with progressive ESCC was also identified in this study. Intercellular interaction analysis based on ligand-receptor pairs revealed the subcluster of malignant cells interacting with CAFs via the MDK-NCL pathway, which was verified by cell proliferation assay and IHC. This indicates that the interaction may be an important hallmark in the early change of tumour microenvironment and serves as a sign of CAF activation to stimulate downstream pathways for facilitating tumour invasion.

CONCLUSION

This study demonstrates the changes of cell subsets and transcriptional levels in human intramucosal ESCC, which may provide unique insights into the development of novel biomarkers and potential intervention strategies.

Towards artificial intelligence to multi-omics characterization of tumor heterogeneity in esophageal cancer

Esophageal cancer is a unique and complex heterogeneous malignancy, with substantial tumor heterogeneity: at the cellular levels, tumors are composed of tumor and stromal cellular components; at the genetic levels, they comprise genetically distinct tumor clones; at the phenotypic levels, cells in distinct microenvironmental niches acquire diverse phenotypic features. This heterogeneity affects almost every process of esophageal cancer progression from onset to metastases and recurrence, etc. Intertumoral and intratumoral heterogeneity are major obstacles in the treatment of esophageal cancer, but also offer the potential to manipulate the heterogeneity themselves as a new therapeutic strategy. The high-dimensional, multi-faceted characterization of genomics, epigenomics, transcriptomics, proteomics, metabonomics, etc. of esophageal cancer has opened novel horizons for dissecting tumor heterogeneity. Artificial intelligence especially machine learning and deep learning algorithms, are able to make decisive interpretations of data from multi-omics layers. To date, artificial intelligence has emerged as a promising computational tool for analyzing and dissecting esophageal patient-specific multi-omics data. This review provides a comprehensive review of tumor heterogeneity from a multi-omics perspective. Especially, we discuss the novel techniques single-cell sequencing and spatial transcriptomics, which have revolutionized our understanding of the cell compositions of esophageal cancer and allowed us to determine novel cell types. We focus on the latest advances in artificial intelligence in integrating multi-omics data of esophageal cancer. Artificial intelligence-based multi-omics data integration computational tools exert a key role in tumor heterogeneity assessment, which will potentially boost the development of precision oncology in esophageal cancer.

Early development of esophageal squamous cell cancer: Stem cells, cellular origins and early clone evolution

Esophageal squamous cell carcinoma (ESCC), a highly malignant cancer with poor prognosis, is an example of the classical view of cancer development based on stem cell origin and multistep progression. In the past five years, the applications of large-scale sequencing and single-cell sequencing have expanded to human esophageal normal tissues and precancerous lesions, which, coupled with the application of transgenic lineage tracing technology in mouse models, has provided a more comprehensive and detailed understanding of esophageal stem cell heterogeneity and early clonal evolution of ESCC. In this review, we discuss the heterogeneity of esophageal basal-layer stem cells and their potential relationship with cells of ESCC origin. We present evidence that expansion of NOTCH1 mutants may call into play an evolutionarily conserved anti-cancer mechanism and mold the model of early clonal evolution in ESCCs. Finally, we discuss the potential avenues in this context. This review provides a focused understanding of the early development of ESCC, as a background for early tumor detection, intervention, and prevention strategies.

Integrative analysis of bulk and single-cell gene expression profiles to identify tumor-associated macrophage-derived CCL18 as a therapeutic target of esophageal squamous cell carcinoma

BACKGROUND

Esophageal squamous cell carcinoma (ESCC) is a common gastrointestinal malignancy with poor patient prognosis. Current treatment for ESCC, including immunotherapy, is only beneficial for a small subset of patients. Better characterization of the tumor microenvironment (TME) and the development of novel therapeutic targets are urgently needed.

METHODS

In the present study, we hypothesized that integration of single-cell transcriptomic sequencing and large microarray sequencing of ESCC biopsies would reveal the key cell subtypes and therapeutic targets that determine the prognostic and tumorigenesis of ESCC. We characterized the gene expression profiles, gene sets enrichment, and the TME landscape of a microarray cohort including 84 ESCC tumors and their paired peritumor samples. We integrated single-cell transcriptomic sequencing and bulk microarray sequencing of ESCC to reveal key cell subtypes and druggable targets that determine the prognostic and tumorigenesis of ESCC. We then designed and screened a blocking peptide targeting Chemokine C-C motif ligand 18 (CCL18) derived from tumor associated macrophages and validated its potency by MTT assay. The antitumor activity of CCL18 blocking peptide was validated in vivo by using 4-nitroquinoline-1-oxide (4-NQO) induced spontaneous ESCC mouse model.

RESULTS

Comparative gene expression and cell-cell interaction analyses revealed dysregulated chemokine and cytokine pathways during ESCC carcinogenesis. TME deconvolution and cell interaction analyses allow us to identify the chemokine CCL18 secreted by tumor associated macrophages could promote tumor cell proliferation via JAK2/STAT3 signaling pathway and lead to poor prognosis of ESCC. The peptide Pep3 could inhibit the proliferation of EC-109 cells promoted by CCL18 and significantly restrain the tumor progression in 4-NQO-induced spontaneous ESCC mouse model.

CONCLUSIONS

For the first time, we discovered and validated that CCL18 blockade could significantly prevent ESCC progression. Our study revealed the comprehensive cell-cell interaction network in the TME of ESCC and provided novel therapeutic targets and strategies to ESCC treatment.

Tumor Niche Network-Defined Subtypes Predict Immunotherapy Response of Esophageal Squamous Cell Cancer

Despite the promising outcomes of immune checkpoint blockade (ICB), resistance to ICB presents a new challenge. Therefore, selecting patients for specific ICB applications is crucial for maximizing therapeutic efficacy. Herein we curated 69 human esophageal squamous cell cancer (ESCC) patients' tumor microenvironment (TME) single-cell transcriptomic datasets to subtype ESCC. Integrative analyses of the cellular network transcriptional signatures of T cells, myeloid cells, and fibroblasts define distinct ESCC subtypes characterized by T cell exhaustion, Interferon (IFN) a/b signaling, TIGIT enrichment, and specific marker genes. Furthermore, this approach classifies ESCC patients into ICB responders and non-responders, as validated by liquid biopsy single-cell transcriptomics. Our study stratifies ESCC patients based on TME transcriptional network, providing novel insights into tumor niche remodeling and predicting ICB responses in ESCC patients.

4‑Methoxydalbergione inhibits esophageal carcinoma cell proliferation and migration by inactivating NF‑κB

4‑Methoxydalbergione (4‑MD) can inhibit the progression of certain types of cancer; however, its effects on esophageal cancer (EC) remain unclear. The present study aimed to investigate the inhibitory effect of 4‑MD on EC and its molecular mechanism. ECA‑109 and KYSE‑105 cells were treated with or without lipopolysaccharide (LPS) and 4‑MD. Cell Counting Kit‑8 and colony formation assays were used to analyze cell proliferation. Wound healing assay was performed to evaluate cell migration. ELISA and western blotting were performed to measure the expression levels of NF‑κB and inflammatory cytokines. In cells treated with 4‑MD, proliferation and migration were significantly inhibited, the levels of inflammatory cytokines were downregulated and the NF‑κB signaling pathway was inactivated. Notably, proliferation, migration, inflammation and NF‑κB were promoted by LPS, whereas 4‑MD reversed the increases induced by LPS in EC cells. In conclusion, 4‑MD may attenuate the proliferation and migration of EC cells by inactivating the NF‑κB signaling pathway.

A Systemic and Integrated Analysis of p63-Driven Regulatory Networks in Mouse Oral Squamous Cell Carcinoma

Oral squamous cell carcinoma (OSCC) is the most common malignancy of the oral cavity and is linked to tobacco exposure, alcohol consumption, and human papillomavirus infection. Despite therapeutic advances, a lack of molecular understanding of disease etiology, and delayed diagnoses continue to negatively affect survival. The identification of oncogenic drivers and prognostic biomarkers by leveraging bulk and single-cell RNA-sequencing datasets of OSCC can lead to more targeted therapies and improved patient outcomes. However, the generation, analysis, and continued utilization of additional genetic and genomic tools are warranted. Tobacco-induced OSCC can be modeled in mice via 4-nitroquinoline 1-oxide (4NQO), which generates a spectrum of neoplastic lesions mimicking human OSCC and upregulates the oncogenic master transcription factor p63. Here, we molecularly characterized established mouse 4NQO treatment-derived OSCC cell lines and utilized RNA and chromatin immunoprecipitation-sequencing to uncover the global p63 gene regulatory and signaling network. We integrated our p63 datasets with published bulk and single-cell RNA-sequencing of mouse 4NQO-treated tongue and esophageal tumors, respectively, to generate a p63-driven gene signature that sheds new light on the role of p63 in murine OSCC. Our analyses reveal known and novel players, such as COTL1, that are regulated by p63 and influence various oncogenic processes, including metastasis. The identification of new sets of potential biomarkers and pathways, some of which are functionally conserved in human OSCC and can prognosticate patient survival, offers new avenues for future mechanistic studies.

Deciphering functional tumor states at single-cell resolution

Within a tumor, cancer cells exist in different states that are associated with distinct tumor functions, including proliferation, differentiation, invasion, metastasis, and resistance to anti-cancer therapy. The identification of the gene regulatory networks underpinning each state is essential for better understanding functional tumor heterogeneity and revealing tumor vulnerabilities. Here, we review the different studies identifying tumor states by single-cell sequencing approaches and the mechanisms that promote and sustain these functional states and regulate their transitions. We also describe how different tumor states are spatially distributed and interact with the specific stromal cells that compose the tumor microenvironment. Finally, we discuss how the understanding of tumor plasticity and transition states can be used to develop new strategies to improve cancer therapy.

Dietary Inflammatory Nutrients and Esophageal Squamous Cell Carcinoma Risk: A Case-Control Study

We conducted a case-control study (532 cases and 532 control) in Chinese adults to investigate the independent and interactive effects of dietary nutrients (pro- or anti-inflammation) on Esophageal Squamous Cell Carcinoma (ESCC) risk. Dietary data were collected using a food questionnaire survey that included 171 items. Two algorithms, the Least Absolute Shrinkage and Selector Operation (LASSO) and Bayesian Kernel Machine Regression (BKMR) were employed to select indicators and evaluate the interactive effect of nutrients' mixture on ESCC risk. Thirteen nutrients were selected, including three pro-inflammatory nutrients (protein, fat and carbohydrate) and ten anti-inflammatory nutrients (fiber, Vitamin A, riboflavin, niacin, Vitamin C, Fe, Se, MUFA, n-3 PUFA and n-6 PUFA). Single-exposure effects of fat, carbohydrate and fiber significantly contributed to ESCC risk. The pro-inflammatory nutrients' submodel discovered that the combined effect was statistically associated with increased ESCC risk. In addition, a higher fat level was significantly associated with ESCC risk. On the contrary, for fiber and riboflavin, the anti-inflammatory nutrients' submodel delineated a significant negative effect on the risk of ESCC. Our result implies that dietary nutrients and their inflammatory traits significantly impacted ESCC occurrence. Additional studies are warranted to verify our findings.