Genomic analysis of oesophageal squamous-cell carcinoma identifies alcohol drinking-related mutation signature and genomic alterations

ZNF468-mediated epigenetic upregulation of VEGF-C facilitates lymphangiogenesis and lymphatic metastasis in ESCC via PI3K/Akt and ERK1/2 signaling pathways

PURPOSE

Dysfunctional lymphangiogenesis is pivotal for various pathological processes including tumor lymph node metastasis which is a crucial cause of therapeutic failure for ESCC. In this study, we aim to elucidate the molecular mechanisms and clinical relevance of Zinc-finger protein ZNF468 in lymphangiogenesis and lymphatic metastasis in ESCC.

METHODS

Immunohistochemistry, Western blot, Kaplan-Meier plotter analysis and Gene Set Enrichment Analysis were preformed to detect the association of ZNF468 with lymphangiogenesis and poor prognosis in ESCC patients. Foot-pads lymph node metastasis model, tube formation assay, 3D-culture assay and invasion assay were preformed to verify the effect of ZNF468 on lymphangiogenesis and lymph node metastasis. CUT&Tag analysis, immunoprecipitation and mass spectrometry analysis and ChIP-PCR assay were preformed to study the molecular mechanisms of ZNF468 in lymphangiogenesis.

RESULTS

We found that ectopic expression of ZNF468 was correlated with higher microlymphatic vessel density in ESCC tissues, leading to poorer prognosis of ESCC patients. ZNF468 enhanced the capacity of lymphangiogenesis and promoted lymphatic metastasis in ESCC both in vitro and in vivo. However, silencing ZNF468 reversed these phenotypes in ESCC. Mechanically, we demonstrated that ZNF468 recruits the histone modification factors (PRMT1/HAT1) to increase the levels of H4R2me2a and H3K9ac, which then leads to the recruitment of the transcription initiation complex on the VEGF-C promoter, ultimately promoting the upregulation of VEGF-C transcription. Strikingly, the promoting effect of lymphatic metastasis induced by ZNF468 in ESCC was abrogated by targeting PRMT1 using Arginine methyltransferase inhibitor-1 or silencing VEGF-C. Furthermore, we found that the activation of PI3K/AKT and ERK1/2 signaling is required for ZNF468-medicated lymphatic metastasis in ESCC. Importantly, the clinical relevance between ZNF468 and VEGF-C were confirmed not only in ESCC samples and but also in multiple cancer types.

CONCLUSION

Our results identified a precise mechanism underlying ZNF468-induced epigenetic upregulation of VEGF-C in facilitating lymphangiogenesis and lymph node metastasis of ESCC, which might provide a novel prognostic biomarker and potential therapeutic for ESCC patients.

The essential roles of lncRNAs/PI3K/AKT axis in gastrointestinal tumors

The role of long noncoding RNA (lncRNA) in tumors, particularly in gastrointestinal tumors, has gained significant attention. Accumulating evidence underscores the interaction between various lncRNAs and diverse molecular pathways involved in cancer progression. One such pivotal pathway is the PI3K/AKT pathway, which serves as a crucial intracellular mechanism maintaining the balance among various cellular physiological processes for normal cell growth and survival. Frequent dysregulation of the PI3K/AKT pathway in cancer, along with aberrant activation, plays a critical role in driving tumorigenesis. LncRNAs modulate the PI3K/AKT signaling pathway through diverse mechanisms, primarily by acting as competing endogenous RNA to regulate miRNA expression and associated genes. This interaction significantly influences fundamental biological behaviors such as cell proliferation, metastasis, and drug resistance. Abnormal expression of numerous lncRNAs in gastrointestinal tumors often correlates with clinical outcomes and pathological features in patients with cancer. Additionally, these lncRNAs influence the sensitivity of tumor cells to chemotherapy in multiple types of gastrointestinal tumors through the abnormal activation of the PI3K/AKT pathway. These findings provide valuable insights into the mechanisms underlying gastrointestinal tumors and potential therapeutic targets. However, gastrointestinal tumors remain a significant global health concern, with increasing incidence and mortality rates of gastrointestinal tumors over recent decades. This review provides a comprehensive summary of the latest research on the interactions of lncRNA and the PI3K/AKT pathway in gastrointestinal tumor development. Additionally, it focuses on the functions of lncRNAs and the PI3K/AKT pathway in carcinogenesis, exploring expression profiles, clinicopathological characteristics, interaction mechanisms with the PI3K/AKT pathway, and potential clinical applications.

The role of distinct APOBEC/ADAR mRNA levels in mutational signatures linked to aging and ultraviolet radiation

The APOBEC/AID family is known for its mutator activity, and recent evidence also supports the potential impact of ADARs. Furthermore, the mutator impacts of APOBEC/ADAR mutations have not yet been investigated. Assessment of pancancer TCGA exomes identified enriched somatic variants among exomes with nonsynonymous APOBEC1, APOBEC3B, APOBEC3C, ADAR, and ADARB1 mutations, compared to exomes with synonymous ones. Principal component (PC) analysis reduced the number of potential players to eight in cancer exomes/genomes, and to five in cancer types. Multivariate regression analysis was used to assess the impact of the PCs on each COSMIC mutational signature among pancancer exomes/genomes and particular cancers, identifying several novel links, including SBS17b, SBS18, and ID7 mainly determined by APOBEC1 mRNA levels; SBS40, ID1, and ID2 by age; SBS3 and SBS16 by APOBEC3A/APOBEC3B mRNA levels; ID5 and DBS9 by DNA repair/replication (DRR) defects; and SBS7a-d, SBS38, ID4, ID8, ID13, and DBS1 by ultraviolet (UV) radiation/ADARB1 mRNA levels. APOBEC/ADAR mutations appeared to potentiate the impact of DRR defects on several mutational signatures, and some factors seemed to inversely affect certain signatures. These findings potentially implicate certain APOBEC/ADAR mutations/mRNA levels in distinct mutational signatures, particularly APOBEC1 mRNA levels in aging-related signatures and ADARB1 mRNA levels in UV radiation-related signatures.

ALDH2 mutations and defense against genotoxic aldehydes in cancer and inherited bone marrow failure syndromes

Reactive aldehydes, for instance, formaldehyde and acetaldehyde, are important endogenous or environmental mutagens by virtue of their abilities to produce a DNA lesion called interstrand crosslink (ICL). Aldehyde-metabolizing enzymes such as aldehyde dehydrogenases (ALDHs) and the Fanconi anemia (FA) pathway constitute the main defense lines against aldehyde-induced genotoxicity. Biallelic mutations of genes in any one of the FA complementation groups can impair the ICL repair mechanism and cause FA, a heterogeneous disorder manifested by bone marrow failure (BMF), congenital abnormality and a strong predisposition to cancer. The defective ALDH2 polymorphism rs671 (ALDH2*2) is a known risk and prognostic factor for alcohol drinking-associated cancers. Recent studies suggest that it also promotes BMF and cancer development in FA, and its combination with alcohol dehydrogenase 5 (ADH5) mutations causes aldehyde degradation deficiency syndrome (ADDS), also known by its symptoms as aplastic anemia, mental retardation, and dwarfism syndrome. ALDH2*2 and another pathogenic variant in the alcohol-metabolizing pathway, ADH1B1*1, is prevalent among East Asians. Also, other ALDH2 genotypes with disease-modifying potentials have lately been identified in different populations. Therefore, it would be appropriate to summarize current knowledge of genotoxic aldehydes and defense mechanisms against them to shed new light on the pathogenic effects of ALDH2 variants together with other genetic and environmental modifiers on cancer and inherited BMF syndromes. Lastly, we also presented potential treatment strategies for FA, ADDS and cancer based on the manipulation of aldehyde-induced genotoxicity.

Prognostic value of genomic mutation signature associated with immune microenvironment in southern Chinese patients with esophageal squamous cell carcinoma

The genomic landscape of esophageal squamous cell cancer (ESCC), as well as its impact on the regulation of immune microenvironment, is not well understood. Thus, tumor samples from 92 patients were collected from two centers and subjected to targeted-gene sequencing. We identified frequently mutated genes, including TP53, KMT2C, KMT2D, LRP1B, and FAT1. The most frequent mutation sites were ALOX12B (c.1565C > T), SLX4 (c.2786C > T), LRIG1 (c.746A > G), and SPEN (c.6915_6917del) (6.5%). Pathway analysis revealed dysregulation of cell cycle regulation, epigenetic regulation, PI3K/AKT signaling, and NOTCH signaling. A 17-mutated gene-related risk model was constructed using random survival forest analysis and showed significant prognostic value in both our cohort and the validation cohort. Based on the Estimation of Stromal and Immune cells in Malignant Tumor tissues using Expression (ESTIMATE) algorithm, the Tumor Immune Dysfunction and Exclusion (TIDE) algorithm, and the MCPcounter algorithm, we found that the risk score calculated by the risk model was significantly correlated with stimulatory immune checkpoints (TNFSF4, ITGB2, CXCL10, CXCL9, and BTN3A1; p < 0.05). Additionally, it was significantly associated with markers that are important in predicting response to immunotherapy (CD274, IFNG, and TAMM2; p < 0.05). Furthermore, the results of immunofluorescence double staining showed that patients with high risk scores had a significantly higher level of M2 macrophage than those with low risk scores (p < 0.05). In conclusion, our study provides insights into the genomic landscape of ESCC and highlights the prognostic value of a genomic mutation signature associated with the immune microenvironment in southern Chinese patients with ESCC.

Systematical identification of key genes and regulatory genetic variants associated with prognosis of esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) stands as a highly lethal malignancy characterized by pronounced recurrence and metastasis, resulting in a bleak 5-year survival rate. Despite extensive investigations, encompassing genome-wide association studies, the identification of robust prognostic markers has remained elusive. In this study, leveraging four independent data sets comprising 404 ESCC patients, we conducted a systematic analysis to unveil pivotal genes influencing overall survival. our meta-analysis identified 278 genes significantly associated with ESCC prognosis. Further exploration of the prognostic landscape involved an examination of expression quantitative trait loci for these genes, leading to the identification of six tag single nucleotide polymorphisms predictive of overall survival in a cohort of 904 ESCC patients. Notably, functional annotation spotlighted rs11227223, residing in the enhancer region of nuclear paraspeckle assembly transcript 1 (NEAT1), as a crucial variant likely exerting a substantive biological role. Through a series of biochemistry experiments, we conclusively demonstrated that the rs11227223-T allele, indicative of a poorer prognosis, augmented NEAT1 expression. Our results underscore the substantive role of NEAT1 and its regulatory variant in prognostic predictions for ESCC. This comprehensive analysis not only advances our comprehension of ESCC prognosis but also unveils a potential avenue for targeted interventions, offering promise for enhanced clinical outcomes.

Antitumour effects of artesunate via cell cycle checkpoint controls in human oesophageal squamous carcinoma cells

Artesunate (ART), an effective antimalarial semisynthetic derivative of artemisinin, exhibits antitumour properties, but the mechanism(s) involved remain elusive. In this study, we investigated the antitumour effects of ART on human oesophageal squamous cell carcinoma (ESCC) cell lines. Treatment of ESCC cell lines with ART resulted in the production of excessive reactive oxygen species (ROS) that induced DNA damage, reduced cell proliferation and inhibited clonogenicity via G1-S cell cycle arrest and/or apoptosis in vitro. The administration of ART to nude mice with ESCC cell xenografts inhibited tumour formation in vivo. However, the cytotoxicity of ART strongly differed among the ESCC cell lines tested. Transcriptomic profiling revealed that although the expression of large numbers of genes in ESCC cell lines was affected by ART treatment, these genes could be functionally clustered into pathways involved in regulating cell cycle progression, DNA metabolism and apoptosis. We revealed that p53 and Cdk4/6-p16-Rb cell cycle checkpoint controls were critical determinants required for mediating ART cytotoxicity in ESCC cell lines. Specifically, KYSE30 cells with p53Mut/p16Mut were the most sensitive to ART, KYSE150 and KYSE180 cells with p53Mut/p16Nor exhibited intermediate responses to ART, and Eca109 cells with p53Nor/p16Nor exhibited the most resistance to ATR. Consistently, perturbation of p53 expression using RNA interference (RNAi) and/or Cdk4/6 activity using the inhibitor palbociclib altered ART cytotoxicity in KYSE30 cells. Given that the p53 and Cdk4/6-cyclin D1-p16-Rb genes are commonly mutated in ESCC, our results potentially shed new light on neoadjuvant chemotherapy strategies for ESCC.

Immune-tumor interaction dictates spatially directed evolution of esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is a poor-prognostic cancer type with extensive intra- and inter-patient heterogeneity in both genomic variations and tumor microenvironment (TME). However, the patterns and drivers of spatial genomic and microenvironmental heterogeneity of ESCC remain largely unknown. Here, we generated a spatial multi-omic atlas by whole-exome, transcriptome, and methylome sequencing of 507 tumor samples from 103 patients. We identified a novel tumor suppressor PREX2, accounting for 22% of ESCCs with frequent somatic mutations or hyper-methylation, which promoted migration and invasion of ESCC cells in vitro. Analysis of the TME and quantification of subclonal expansion indicated that ESCCs undergo spatially directed evolution, where subclones mostly originated from the tumor center but had a biased clonal expansion to the upper direction of the esophagus. Interestingly, we found upper regions of ESCCs often underwent stronger immunoediting with increased selective fitness, suggesting more stringent immune selection. In addition, distinct TMEs were associated with variable genomic and clinical outcomes. Among them, hot TME was associated with high immune evasion and subclonal heterogeneity. We also found that immunoediting, instead of CD8+ T cell abundance, acts as an independent prognostic factor of ESCCs. Importantly, we found significant heterogeneity in previously considered potential therapeutic targets, as well as BRCAness characteristics in a subset of patients, emphasizing the importance of focusing on heterogeneity in ESCC targeted therapy. Collectively, these findings provide novel insights into the mechanisms of the spatial evolution of ESCC and inform precision therapeutic strategies.

COL8A1 Regulates Esophageal Squamous Carcinoma Proliferation and Invasion Through PI3K/AKT Pathway

PURPOSE

Esophageal squamous carcinoma (ESCC) is a gastrointestinal malignancy with a high mortality, but the tumorigenesis is still unclear, restricting the target therapy development of ESCC. We explored the role of COL8A1 in ESCC development.

METHODS

Tissue microarrays were used to investigate the expression level of COL8A1 in ESCC tissues. The association between COL8A1 and the overall survival of ESCC patients was assessed. The effect of differential COL8A1 expression on tumor growth was investigated by the xenograft model. The regulation of COL8A1 on tumor growth, migration, and invasion was studied by using ESCC cell lines. The signal transduction pathways involved in COL8A1 were bioinformatically profiled and validated.

RESULTS

The COL8A1 was significantly expressed in cancerous tissues and was associated with poor prognosis in patients with ESCC. In vivo, the tumor growth obviously declined after inhibition of the COL8A1 expression. The abilities of cell proliferation and invasion were both decreased when the expression of COL8A1 was knockdown in ESCC cell line. Furthermore, we found the inactivation of the PI3K/AKT pathway that was mediated by knockdown of COL8A1 in ESCC cells, which was reversed with COL8A1 overexpression, whereas the cell proliferation and invasion ability were restored.

CONCLUSIONS

This is the first report that COL8A1 promote ESCC progression, which hopefully will provide a theoretical basis for clinical targeting of ESCC.

Characterization of macrophages in head and neck squamous cell carcinoma and development of MRG-based risk signature

Macrophages are immune cells in the TME that can not only inhibit angiogenesis, extracellular matrix remodeling, cancer cell proliferation, and metastasis but also mediate the phagocytosis and killing of cancer cells after activation, making them key targets in anti-tumor immunotherapy. However, there is little research on macrophages and their relation to disease prognosis in HNSCC. Initially, we collected scRNA-seq, bulk RNA-seq, and clinical data. Subsequently, we identified macrophages and distinguished MRGs. Using the K-means algorithm, we performed consensus unsupervised clustering. Next, we used ssGSEA analysis to assess immune cell infiltration in MRG clusters. A risk model was established using multivariate Cox analysis. Then, Kaplan-Meier, ROC curves, univariate and multivariate COX analyses, and C-index was used to validate the predictive power of the signature. The TIDE method was applied to assess the response to immunotherapy in patients diagnosed with HNSCC. In addition, drug susceptibility predictions were made for the GDSC database using the calcPhenotype function. We found that 8 MRGs had prognostic potential. Patients in the MRG group A had a higher probability of survival, and MRG clusters A and B had different characteristics. Cluster A had a higher degree of expression and infiltration in MRG, indicating a closer relationship with MRG. The accuracy of the signature was validated using univariate and multivariate Cox analysis, C-index, and nomogram. Immune landscape analysis found that various immune functions were highly expressed in the low-risk group, indicating an improved response to immunotherapy. Finally, drugs with high sensitivity to HNSCC (such as 5-Fluorouracil, Temozolomide, Carmustine, and EPZ5676) were explored and analyze the malignant characteristics of HNSCC. We constructed a prognostic model using multivariate Cox analysis, consisting of 8 MRGs (TGM2, STC1, SH2D3C, PIK3R3, MAP3K8, ITGA5, ARHGAP4, and AQP1). Patients in the low-risk group may have a higher response to immunotherapy. The more prominent drugs for drug selection are 5-fluorouracil, temozolomide and so on. Malignant features associated with HNSCC include angiogenesis, EMT, and the cell cycle. This study has opened up new prospects for the prognosis, prediction, and clinical treatment strategy of HNSCC.

An effective two-stage NMBzA-induced rat esophageal tumor model revealing that the FAT-Hippo-YAP1 axis drives the progression of ESCC

Rat model of N-nitrosomethylbenzylamine (NMBzA)-induced esophageal squamous cell carcinoma (ESCC) is routinely used to study ESCC initiation, progression and new therapeutic strategies. However, the model is time-consuming and malignant tumor incidences are low. Here, we report the usage of multi-kinase inhibitor sorafenib as a tumor promoter to establish an efficient two-stage NMBzA-induced rat ESCC carcinogenesis model, resulting in increments of tumor incidences and shortened tumor formation times. By establishing the model and applying whole-genome sequencing, we discover that benign papillomas and malignant ESCCs harbor most of the "driver" events found in rat ESCCs (e.g. recurrent mutations in Ras family, the Hippo and Notch pathways and histone modifier genes) and the mutational landscapes of rat and human ESCCs overlap extensively. We generate tumor cell lines derived from NMBzA-induced papillomas and ESCCs, showing that papilloma cells retain more characteristics of normal epithelial cells than carcinoma cells, especially their exhibitions of normal rat cell karyotypes and inabilities of forming tumors in immunodeficient mice. Three-dimensional (3-D) organoid cultures and single cell RNA sequencing (scRNA-seq) indicate that, when compared to control- and papilloma-organoids, ESCC-organoids display salient abnormalities at tissue and single-cell levels. Multi-omic analyses indicate that NMBzA-induced rat ESCCs are accompanied by progressive hyperactivations of the FAT-Hippo-YAP1 axis and siRNA or inhibitors of YAP1 block the growth of rat ESCCs. Taken together, these studies provide a framework of using an effective rat ESCC model to investigate multilevel functional genomics of ESCC carcinogenesis, which justify targeting YAP1 as a therapeutic strategy for ESCC.

Single-cell analysis identified POSTN+ cells associated with the aggressive phenotype and risk of esophageal squamous cell carcinoma

Tumors are intricate and heterogeneous systems characterized by mosaic cancer cell populations with diverse expression profiles. Leveraging single-cell technologies, we employed the Scissor algorithm to delineate an epithelial subpopulation associated with the aggressive phenotype in esophageal squamous cell carcinoma (ESCC). This identified subpopulation exhibited elevated expression of genes involved in critical pathways, such as epithelial-mesenchymal transition and PI3K-Akt. Key signature genes within this subpopulation, namely CAV1, COL3A1, COL6A1, POSTN, and TAGLN, demonstrated significant upregulation concomitant with both tumorigenesis and tumor progression across independent single-cell datasets. Furthermore, we selected 1,450 expression quantitative trait loci of the top 62 signature genes of this cell subpopulation to investigate their potential in predicting ESCC risk. The results showed that the POSTN loci were predominantly associated with ESCC susceptibility. Through functional annotation and replication analyses, we identified that the rs1028728 in the POSTN promoter was significantly associated with increased ESCC risk in 7,049 ESCC cases and 8,063 controls (odds ratio = 1.29, 95% confidence interval: 1.18-1.42, p = 4.03 × 10-8). Subsequent biochemical experiments showed that the rs1028728[T] allele enhanced POSTN expression by affecting the binding of PRRX1 in the POSTN promoter. In summary, our meticulous single-cell analysis delineates an invasive epithelial subpopulation in ESCC, with POSTN emerging as an important marker for the aggressive phenotype. These findings offer more insights into potential strategies for the prevention and intervention of ESCC, enriching our understanding of this complex cancer landscape.

Links between melanoma germline risk loci, driver genes and comorbidities: insight from a tissue-specific multi-omic analysis

Genome-wide association studies (GWAS) have associated 76 loci with the risk of developing melanoma. However, understanding the molecular basis of such associations has remained a challenge because most of these loci are in non-coding regions of the genome. Here, we integrated data on epigenomic markers, three-dimensional (3D) genome organization, and expression quantitative trait loci (eQTL) from melanoma-relevant tissues and cell types to gain novel insights into the mechanisms underlying melanoma risk. This integrative approach revealed a total of 151 target genes, both near and far away from the risk loci in linear sequence, with known and novel roles in the etiology of melanoma. Using protein-protein interaction networks, we identified proteins that interact-directly or indirectly-with the products of the target genes. The interacting proteins were enriched for known melanoma driver genes. Further integration of these target genes into tissue-specific gene regulatory networks revealed patterns of gene regulation that connect melanoma to its comorbidities. Our study provides novel insights into the biological implications of genetic variants associated with melanoma risk.

Proteogenomic characterization of primary colorectal cancer and metastatic progression identifies proteome-based subtypes and signatures

Metastatic progression of colorectal adenocarcinoma (CRC) remains poorly understood and poses significant challenges for treatment. To overcome these challenges, we performed multiomics analyses of primary CRC and liver metastases. Genomic alterations, such as structural variants or copy number alterations, were enriched in oncogenes and tumor suppressor genes and increased in metastases. Unsupervised mass spectrometry-based proteomics of 135 primary and 123 metastatic CRCs uncovered distinct proteomic subtypes, three each for primary and metastatic CRCs, respectively. Integrated analyses revealed that hypoxia, stemness, and immune signatures characterize these 6 subtypes. Hypoxic CRC harbors high epithelial-to-mesenchymal transition features and metabolic adaptation. CRC with a stemness signature shows high oncogenic pathway activation and alternative telomere lengthening (ALT) phenotype, especially in metastatic lesions. Tumor microenvironment analysis shows immune evasion via modulation of major histocompatibility complex (MHC) class I/II and antigen processing pathways. This study characterizes both primary and metastatic CRCs and provides a large proteogenomics dataset of metastatic progression.

Genomic alterations in oral multiple primary cancers

Oral squamous cell carcinoma (OSCC) is the predominant type of oral cancer, while some patients may develop oral multiple primary cancers (MPCs) with unclear etiology. This study aimed to investigate the clinicopathological characteristics and genomic alterations of oral MPCs. Clinicopathological data from patients with oral single primary carcinoma (SPC, n = 202) and oral MPCs (n = 34) were collected and compared. Copy number alteration (CNA) analysis was conducted to identify chromosomal-instability differences among oral MPCs, recurrent OSCC cases, and OSCC patients with lymph node metastasis. Whole-exome sequencing was employed to identify potential unique gene mutations in oral MPCs patients. Additionally, CNA and phylogenetic tree analyses were used to gain preliminary insights into the molecular characteristics of different primary tumors within individual patients. Our findings revealed that, in contrast to oral SPC, females predominated the oral MPCs (70.59%), while smoking and alcohol use were not frequent in MPCs. Moreover, long-term survival outcomes were poorer in oral MPCs. From a CNA perspective, no significant differences were observed between oral MPCs patients and those with recurrence and lymph node metastasis. In addition to commonly mutated genes such as CASP8, TP53 and MUC16, in oral MPCs we also detected relatively rare mutations, such as HS3ST6 and RFPL4A. Furthermore, this study also demonstrated that most MPCs patients exhibited similarities in certain genomic regions within individuals, and distinct differences of the similarity degree were observed between synchronous and metachronous oral MPCs.

DNA damage response alterations in clear cell renal cell carcinoma: clinical, molecular, and prognostic implications

BACKGROUND

DNA damage repair (DDR) pathways modulate cancer risk, progression, and therapeutic responses. Nonetheless, the characteristics and significance of DDR alterations in clear cell renal cell carcinoma (ccRCC) remain undefined. This study aimed to explore the predictive role, molecular mechanism, and tumor immune profile of DDR genes in ccRCC.

METHODS

We prospectively sequenced 757 tumors and matched blood DNA samples from Chinese patients with ccRCC using next-generation sequencing (NGS) and analyzed data from 537 patients from The Cancer Genome Atlas (TCGA). A comprehensive analysis was performed.

RESULTS

Fifty-two percent of Chinese patients with ccRCC harbored DDR gene mutations and 57% of TCGA patients. The immunotherapy treatment prognosis of patients with DDR gene mutations was superior to that of patients without DDR gene mutations (p = 0.047). DDR gene mutations were associated with more gene mutations and a higher tumor mutation load (TMB, p < 0.001). Moreover, patients with DDR gene mutations have a distinct mutational signature compared with those with wild-type DDR. Furthermore, the DDR-mut group had elevated neoantigen load (including single-nucleotide variants (SNV) and indel neoantigen load, p = 0.037 and p = 0.002, respectively), TCR Shannon (p = 0.025), and neutrophils (p = 0.010). DDR gene mutations exhibited a distinct immune profile with significantly higher expression levels of TNFSF9, CD70, ICAM1, and indoleamine-2,3-dioxygenase (IDO) and lower expression levels of VTCN1 and IL12A.

CONCLUSIONS

Our data suggest that the detection of somatic mutations in DDR genes can predict the efficacy of immunotherapy in patients with ccRCC. Furthermore, we revealed the unique molecular and immune mechanisms underlying ccRCC with DDR gene mutations.

Molecular landscapes of oral cancers of unknown etiology

The incidence of the mobile tongue cancer in young patients has been rising. This oral cancer (OC) type has no identified risk factors (NIRF), no established molecular markers and is not yet recognized as a distinct clinical entity. To understand this emerging malignancy, we innovatively analyzed the public head and neck cancer multi-omics data. We identified mutational signatures that successfully stratified 307 OC and 109 laryngeal cancer cases according to their clinico-pathological characteristics. The NIRF OCs exhibited significantly increased activities of endogenous clock-like and APOBEC-associated mutagenesis, alongside specific cancer driver gene mutations, distinct methylome patterns and prominent antimicrobial transcriptomic responses. Furthermore, we show that mutational signature SBS16 in OCs reflects the combined effects of alcohol drinking and tobacco smoking. Our study characterizes the unique disease histories and molecular programs of the NIRF OCs revealing that this emerging cancer subtype is likely driven by increased endogenous mutagenesis correlated with responses to microbial insults.

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.

Ethanol induces replication fork stalling and membrane stress in immortalized laryngeal cells

Although ethanol is a class I carcinogen and is linked to more than 700,000 cancer incidences, a clear understanding of the molecular mechanisms underlying ethanol-related carcinogenesis is still lacking. Further understanding of ethanol-related cell damage can contribute to reducing or treating alcohol-related cancers. Here, we investigated the effects of both short- and long-term exposure of human laryngeal epithelial cells to different ethanol concentrations. RNA sequencing shows that ethanol altered gene expression patterns in a time- and concentration-dependent way, affecting genes involved in ribosome biogenesis, cytoskeleton remodeling, Wnt signaling, and transmembrane ion transport. Additionally, ethanol induced a slower cell proliferation, a delayed cell cycle progression, and replication fork stalling. In addition, ethanol exposure resulted in morphological changes, which could be associated with membrane stress. Taken together, our data yields a comprehensive view of molecular changes associated with ethanol stress in epithelial cells of the upper aerodigestive tract.

Genomic alterations driving precancerous to cancerous lesions in esophageal cancer development

Esophageal squamous cell carcinoma (ESCC) develops through a series of increasingly abnormal precancerous lesions. Previous studies have revealed the striking differences between normal esophageal epithelium and ESCC in copy number alterations (CNAs) and mutations in genes driving clonal expansion. However, due to limited data on early precancerous lesions, the timing of these transitions and which among them are prerequisites for malignant transformation remained unclear. Here, we analyze 1,275 micro-biopsies from normal esophagus, early and late precancerous lesions, and esophageal cancers to decipher the genomic alterations at each stage. We show that the frequency of TP53 biallelic inactivation increases dramatically in early precancerous lesion stage while CNAs and APOBEC mutagenesis substantially increase at late stages. TP53 biallelic loss is the prerequisite for the development of CNAs of genes in cell cycle, DNA repair, and apoptosis pathways, suggesting it might be one of the earliest steps initiating malignant transformation.

The mutational pattern of homologous recombination repair genes in urothelial carcinoma and its correlation with immunotherapeutic response

BACKGROUND

The mutational pattern of homologous recombination repair (HRR)-associated gene alterations in Chinese urothelial carcinoma (UC) necessitates comprehensive sequencing efforts, and the clinical implications of HRR gene mutations in UC remain to be elucidated.

MATERIALS AND METHODS

We delineated the mutational landscape of 343 Chinese UC patients from West China Hospital and 822 patients from The Cancer Genome Atlas (TCGA) using next-generation sequencing (NGS). Data from 182 metastatic UC patients from MSK-IMPACT cohort were used to assess the association between HRR mutations and immunotherapy efficacy. Comprehensive transcriptomic analysis was performed to explore the impact of HRR mutations on tumor immune microenvironment.

RESULTS

Among Chinese UC patients, 34% harbored HRR gene mutations, with BRCA2, ATM, BRCA1, CDK12, and RAD51C being the most prevalently mutated genes. Mutational signatures contributing to UC differed between patients with and without HRR mutations. Signature 22 for exposure to aristolochic acid was only observed in Chinese UC patients. The presence of HRR mutations was correlated with higher tumor mutational burden, neoantigen burden, and PD-L1 expression. Importantly, patients with HRR mutations exhibited significantly improved prognosis following immunotherapy compared to those without HRR mutations.

CONCLUSIONS

Our findings provide valuable insights into the genomic landscape of Chinese UC patients and underscore the molecular rationale for utilizing immunotherapy in UC patients with HRR mutations.

Progenitor-like exhausted SPRY1+CD8+ T cells potentiate responsiveness to neoadjuvant PD-1 blockade in esophageal squamous cell carcinoma

Neoadjuvant immune checkpoint blockade (ICB) demonstrates promise in operable esophageal squamous cell carcinoma (ESCC), but lacks available efficacy biomarkers. Here, we perform single-cell RNA-sequencing of tumors from patients with ESCC undergoing neoadjuvant ICB, revealing a subset of exhausted CD8+ T cells expressing SPRY1 (CD8+ Tex-SPRY1) that displays a progenitor exhausted T cell (Tpex) phenotype and correlates with complete response to ICB. We validate CD8+ Tex-SPRY1 cells as an ICB-specific predictor of improved response and survival using independent ICB-/non-ICB cohorts and demonstrate that expression of SPRY1 in CD8+ T cells enforces Tpex phenotype and enhances ICB efficacy. Additionally, CD8+ Tex-SPRY1 cells contribute to proinflammatory phenotype of macrophages and functional state of B cells, which thereby promotes antitumor immunity by enhancing CD8+ T cell effector functions. Overall, our findings unravel progenitor-like CD8+ Tex-SPRY1 cells' role in effective responses to ICB for ESCC and inform mechanistic biomarkers for future individualized immunotherapy.

Multi-omics analyses reveal spatial heterogeneity in primary and metastatic oesophageal squamous cell carcinoma

BACKGROUND

Biopsies obtained from primary oesophageal squamous cell carcinoma (ESCC) guide diagnosis and treatment. However, spatial intra-tumoral heterogeneity (ITH) influences biopsy-derived information and patient responsiveness to therapy. Here, we aimed to elucidate the spatial ITH of ESCC and matched lymph node metastasis (LNmet ).

METHODS

Primary tumour superficial (PTsup ), deep (PTdeep ) and LNmet subregions of patients with locally advanced resectable ESCC were evaluated using whole-exome sequencing (WES), whole-transcriptome sequencing and spatially resolved digital spatial profiling (DSP). To validate the findings, immunohistochemistry was conducted and a single-cell transcriptomic dataset was analysed.

RESULTS

WES revealed 15.72%, 5.02% and 32.00% unique mutations in PTsup , PTdeep and LNmet , respectively. Copy number alterations and phylogenetic trees showed spatial ITH among subregions both within and among patients. Driver mutations had a mixed intra-tumoral clonal status among subregions. Transcriptome data showed distinct differentially expressed genes among subregions. LNmet exhibited elevated expression of immunomodulatory genes and enriched immune cells, particularly when compared with PTsup (all P < .05). DSP revealed orthogonal support of bulk transcriptome results, with differences in protein and immune cell abundance between subregions in a spatial context. The integrative analysis of multi-omics data revealed complex heterogeneity in mRNA/protein levels and immune cell abundance within each subregion.

CONCLUSIONS

This study comprehensively characterised spatial ITH in ESCC, and the findings highlight the clinical significance of unbiased molecular classification based on multi-omics data and their potential to improve the understanding and management of ESCC. The current practices for tissue sampling are insufficient for guiding precision medicine for ESCC, and routine profiling of PTdeep and/or LNmet should be systematically performed to obtain a more comprehensive understanding of ESCC and better inform treatment decisions.

Noncoding RNA mutations in cancer

Cancer is driven by both germline and somatic genetic changes. Efforts have been devoted to characterizing essential genetic variations in cancer initiation and development. Most attention has been given to mutations in protein-coding genes and associated regulatory elements such as promoters and enhancers. The development of sequencing technologies and in silico and experimental methods has allowed further exploration of cancer predisposition variants and important somatic mutations in noncoding RNAs, mainly for long noncoding RNAs and microRNAs. Association studies including GWAS have revealed hereditary variations including SNPs and indels in lncRNA or miRNA genes and regulatory regions. These mutations altered RNA secondary structures, expression levels, and target recognition and then conferred cancer predisposition to carriers. Whole-exome/genome sequencing comparing cancer and normal tissues has revealed important somatic mutations in noncoding RNA genes. Mutation hotspots and somatic copy number alterations have been identified in various tumor-associated noncoding RNAs. Increasing focus and effort have been devoted to studying the noncoding region of the genome. The complex genetic network of cancer initiation is being unveiled. This article is categorized under: RNA in Disease and Development > RNA in Disease.

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.

Low expression of PINK1 and PARK2 predicts poor prognosis in patients with esophageal squamous cell carcinoma

BACKGROUND

The Parkinson's disease (PD) gene family expression is strongly linked to tumor development and progression; PINK1 and PARK2 are essential members of the PD gene family. However, the relationship between PINK1 and PARK2 and esophageal squamous cell carcinoma (ESCC) remains unknown. This research aims to clarify the prognostic value of PINK1 and PARK2 in ESCC.

METHODS

PINK1 and PARK2 protein levels in 232 ESCC specimens, and 125 matched adjacent normal tissues were detected by immunohistochemistry. The relationship between PINK1 and PARK2 protein expression and clinicopathological features were analyzed. Kaplan-Meier survival analysis was performed to estimate the prognostic value of the PINK1 and PARK2 proteins in patients. Cox univariate and multivariate analyses were used to assess the risk factors affecting the OS for patients with ESCC.

RESULTS

PINK1 and PARK2 had low expression in ESCC. Patients with low PINK1 had worse differentiation and advanced T and TNM stages. Lower PARK2 expression was linked to lymph node metastases and an advanced TNM stage. Furthermore, reduced PINK1 and PARK2 levels were associated with a poor prognosis for ESCC. Cox univariate and multivariate analyses revealed that PINK1, PARK2, and tumor size were closely associated with the prognosis of patients with ESCC, and PARK2 was an independent risk factor for patients with ESCC. Finally, the PINK1 and PARK2 proteins were closely related and shared the same signal pathway.

CONCLUSIONS

PINK1 and PARK2 could work as tumor suppressors in ESCC and are likely to become new treatment targets for ESCC.

Integrative analysis of tertiary lymphoid structures and immune microenvironment in patients with esophageal carcinoma

BACKGROUND

Esophageal squamous cell carcinoma (ESCC) is one of the most common upper gastrointestinal malignancies worldwide. Tertiary lymphoid structures (TLS) are tumor-infiltrating immune cells aggregates coupled with stromal cells which are similar to secondary lymphoid organs. The objective of this study is to explore the predictive effects of two common genes associated with TLS models on prognosis and immunotherapy effects in ESCC patients.

METHODS

Clinical information for ESCC patients in the TCGA(The Cancer Genome Altas) cohort and GSE 53625 were collected. All of the samples were classified as either high score group or low score group based on two TLS signatures, and the association between TLS signatures and survival, clinical indicators, genomic burden, stemness indices analysis, tumor microenvironment and immunotherapy response were performed. Furthermore, the mature TLS was also assessed in ESCC tissue microarray.

RESULTS

In our study, we quantified the score of TLS_9 and TLS_12, respectively, reflecting the different statuses of TLS (TLS_9 = B and T cells in TLSs; TLS_12 = neogenesis of TLSs). Subsequently, we explored the effect of TLS score on ESCC tumor microenvironment quantified by multiple algorithms. We found that a correlation analysis indicated that TLS_9 and TLS_12 were all positively correlated with CD8+ T cell, NK cells, CD4+ T cells, M1 macrophages and so on. Meanwhile, some cells present a different correlation pattern of TLS_9 and TLS_12, including activated CD4+ memory T cells and Tgd cells. Immune-related analysis revealed that the TLS_12 and TLS_9 scores were all positively correlated with immune dysfunction, yet negatively correlated with immune exclusion. Following this, the biological roles of TLS_9 and TLS_12 scores were investigated. Also, we noticed that the TLS score could significantly affect the CAFs infiltration and be associated with the genomic burden and tumor stemness. In addition, we explored the prognostic value of mature TLS through tissue microarray (TMA). Our result displayed ESCC patients with the presence of mature TLS had a better prognosis than ESCC patients without it.

CONCLUSIONS

Our study indicated that ESCC patients with the presence of TLS had better outcomes and an inflamed immune microenvironment. In addition, both TLS-9 and TLS-12 gene signatures could be used as potential biomarkers for the immunotherapy of ESCC patients.

Bioinformatics and network pharmacology identify promotional effects and potential mechanisms of ethanol on esophageal squamous cell carcinoma and experimental validation

Ethanol is an important risk factor for esophageal squamous cell carcinoma (ESCC); however, the molecular mechanisms behind how ethanol promotes ESCC development remain poorly understood. In this study, ethanol-ESCC-associated target genes were constructed and screened using network pharmacology and subjected to Kyoto Encyclopedia of Genes and Genomes (KEGG) and bioinformatics analysis. A mouse ethanol-exposed esophageal cancer model was constructed with 4-nitroquinoline-1-oxide (4-NQO) to assess its survival and tumor lesion status, and the mechanism of ethanol-promoted ESCC lesions was verified by qRT-PCR and Western blotting. The results showed that 126 ethanol-ESCC crossover genes were obtained, which were significantly enriched in the PI3K/AKT signaling pathway. Bioinformatics results showed that the target genes TNF, IL6, IL1β and JUN were highly expressed in esophageal tumor samples and positively correlated with tumor proliferation and apoptosis genes, and the genetic information of these genes was mutated to different degrees. Animal model experiments showed that ethanol decreased the survival rate and aggravated the occurrence of esophageal cancer in mice. qRT-PCR showed that ethanol promoted the expression of TNF, IL6, IL1β and JUN mRNA in mouse esophageal tumor tissues, and Western blotting showed that ethanol promoted p-PI3K and p-AKT protein expression in mouse esophageal tumor tissues. In conclusion, ethanol promotes esophageal carcinogenesis by increasing the expression of TNF, IL6, IL1β and JUN and activating the PI3K/AKT signaling pathway.

Circular RNA circ-FIRRE interacts with HNRNPC to promote esophageal squamous cell carcinoma progression by stabilizing GLI2 mRNA

Increasing evidence has shown that circular RNAs (circRNAs) interact with RNA-binding proteins (RBPs) and promote cancer progression. However, the function and mechanism of the circRNA/RBP complex in esophageal squamous cell carcinoma (ESCC) are still largely unknown. Herein, we first characterized a novel oncogenic circRNA, circ-FIRRE, by RNA sequencing (Ribo-free) profiling of ESCC samples. Furthermore, we observed marked circ-FIRRE overexpression in ESCC patients with high TNM stage and poor overall survival. Mechanistic studies indicated that circ-FIRRE, as a platform, interacts with the heterogeneous nuclear ribonucleoprotein C (HNRNPC) protein to stabilize GLI2 mRNA by directly binding to its 3'-UTR in the cytoplasm, thereby resulting in elevated GLI2 protein expression and subsequent transcription of its target genes MYC, CCNE1, and CCNE2, ultimately contributing to ESCC progression. Moreover, HNRNPC overexpression in circ-FIRRE knockdown cells notably abolished circ-FIRRE knockdown-mediated Hedgehog pathway inhibition and ESCC progression impairment in vitro and in vivo. Clinical specimen results showed that circ-FIRRE and HNRNPC expression was positively correlated with GLI2 expression, which reveals the clear significance of the circ-FIRRE/HNRNPC-GLI2 axis in ESCC. In summary, our results indicate that circ-FIRRE could serve as a valuable biomarker and potential therapeutic target for ESCC and highlight a novel mechanism of the circ-FIRRE/HNRNPC complex in ESCC progression regulation.

Topography of mutational signatures in human cancer

The somatic mutations found in a cancer genome are imprinted by different mutational processes. Each process exhibits a characteristic mutational signature, which can be affected by the genome architecture. However, the interplay between mutational signatures and topographical genomic features has not been extensively explored. Here, we integrate mutations from 5,120 whole-genome-sequenced tumors from 40 cancer types with 516 topographical features from ENCODE to evaluate the effect of nucleosome occupancy, histone modifications, CTCF binding, replication timing, and transcription/replication strand asymmetries on the cancer-specific accumulation of mutations from distinct mutagenic processes. Most mutational signatures are affected by topographical features, with signatures of related etiologies being similarly affected. Certain signatures exhibit periodic behaviors or cancer-type-specific enrichments/depletions near topographical features, revealing further information about the processes that imprinted them. Our findings, disseminated via the COSMIC (Catalog of Somatic Mutations in Cancer) signatures database, provide a comprehensive online resource for exploring the interactions between mutational signatures and topographical features across human cancer.

Comparative transcriptome characterization of esophageal squamous cell carcinoma and adenocarcinoma

Background

Esophageal cancers are primarily categorized as esophageal squamous cell carcinoma (ESCC) and esophageal adenocarcinoma (EAC). While various (epi) genomic alterations associated with tumor development in ESCC and EAC have been documented, a comprehensive comparison of the transcriptomes in these two cancer subtypes remains lacking.

Methods

We collected 551 gene expression profiles from publicly available sources, including normal, ESCC, and EAC tissues or cell lines. Subsequently, we conducted a systematic analysis to compare the transcriptomes of these samples at various levels, including gene expression, promoter activity, alternative splicing (AS), alternative polyadenylation (APA), and gene fusion.

Results

Seven distinct cluster gene expression patterns were identified among the differentially expressed genes in normal, ESCC, and EAC tissues. These patterns were enriched in the PI3K-Akt signaling pathway and the activation of extracellular matrix organization and exhibited repression of epidermal development. Notably, we observed additional genes or unique expression levels enriched in these shared pathways and biological processes related to tumor development and immune activation. In addition to the differentially expressed genes, there was an enrichment of lncRNA co-expression networks and downregulation of promoter activity associated with the repression of epidermal development in both ESCC and EAC. This indicates a common feature between these two cancer subtypes. Furthermore, differential AS and APA patterns in ESCC and EAC appear to partially affect the expression of host genes associated with bacterial or viral infections in these subtypes. No gene fusions were observed between ESCC and EAC, thus highlighting the distinct molecular mechanisms underlying these two cancer subtypes.

Conclusions

We conducted a comprehensive comparison of ESCC and EAC transcriptomes and uncovered shared and distinct transcriptomic signatures at multiple levels. These findings suggest that ESCC and EAC may exhibit common and unique mechanisms involved in tumorigenesis.

Mutational pattern off homologous recombination repair (HRR)-related genes in upper tract urothelial carcinoma

BACKGROUND

Homologous recombination (HR) repair (HRR) has been indicated to be a biomarker for immunotherapy, chemotherapy, and poly-ADP ribose polymerase inhibitors inhibitors (PARPis). Nonetheless, their molecular correlates in upper tract urothelial carcinoma (UTUC) have not been well studied. This study aimed to explore the molecular mechanism and tumor immune profile of HRR genes and the relevance of their prognostic value in patients with UTUC.

MATERIALS AND METHODS

One hundred and ninety-seven tumors and matched blood samples from Chinese UTUC were subjected to next-generation sequencing. A total of 186 patients from The Cancer Genome Atlas were included. Comprehensive analysis was performed.

RESULTS

In Chinese patients with UTUC, 5.01% harbored germline HRR gene mutations, and 1.01% had Lynch syndrome-related genes. A total of 37.6% (74/197) of patients carried somatic or germline HRR gene mutations. There was marked discrepancy in the mutation landscapes, genetic interactions, and driver genes between the HRR-mut cohorts and HRR-wt cohorts. Aristolochic acid signatures and defective DNA mismatch repair signatures only existed in individuals in the HRR-mut cohorts. Inversely, the unknown signature (signature A) and signature SBS55 only existed in patients in the HRR-wt cohorts. HRR gene mutations regulated immune activities by NKT cells, plasmacytoid dendritic cells, hematopoietic stem cell, and M1 macrophages. In patients with local recurrence, patients with HRR gene mutations had poorer DFS rates than patients with wild-type HRR genes.

CONCLUSIONS

Our results imply that the detection of HRR gene mutations can predict recurrence in patients with UC. In addition, this study provides a path to explore the role of HRR-directed therapies, including PARPis, chemotherapy, and immunotherapy.

Aldehyde-Associated Mutagenesis─Current State of Knowledge

Aldehydes are widespread in the environment, with multiple sources such as food and beverages, industrial effluents, cigarette smoke, and additives. The toxic effects of exposure to several aldehydes have been observed in numerous studies. At the molecular level, aldehydes damage DNA, cross-link DNA and proteins, lead to lipid peroxidation, and are associated with increased disease risk including cancer. People genetically predisposed to aldehyde sensitivity exhibit severe health outcomes. In various diseases such as Fanconi's anemia and Cockayne syndrome, loss of aldehyde-metabolizing pathways in conjunction with defects in DNA repair leads to widespread DNA damage. Importantly, aldehyde-associated mutagenicity is being explored in a growing number of studies, which could offer key insights into how they potentially contribute to tumorigenesis. Here, we review the genotoxic effects of various aldehydes, focusing particularly on the DNA adducts underlying the mutagenicity of environmentally derived aldehydes. We summarize the chemical structures of the aldehydes and their predominant DNA adducts, discuss various methodologies, in vitro and in vivo, commonly used in measuring aldehyde-associated mutagenesis, and highlight some recent studies looking at aldehyde-associated mutation signatures and spectra. We conclude the Review with a discussion on the challenges and future perspectives of investigating aldehyde-associated mutagenesis.

Comprehensive Analyses Reveal Effects on Tumor Immune Infiltration and Immunotherapy Response of APOBEC Mutagenesis and Its Molecular Mechanisms in Esophageal Squamous Cell Carcinoma

The apolipoprotein B mRNA editing enzyme catalytic polypeptide (APOBEC) mutagenesis is prevalent in esophageal squamous cell carcinoma (ESCC). However, the functional role of APOBEC mutagenesis has yet to be fully delineated. To address this, we collect matched multi-omics data of 169 ESCC patients and evaluate characteristics of immune infiltration using multiple bioinformatic approaches based on bulk and single-cell RNA sequencing (scRNA-seq) data and verified by functional assays. We find that APOBEC mutagenesis prolongs overall survival (OS) of ESCC patients. The reason for this outcome is probably due to high anti-tumor immune infiltration, immune checkpoints expression and immune related pathway enrichment, such as interferon (IFN) signaling, innate and adaptive immune system. The elevated AOBEC3A (A3A) activity paramountly contributes to the footprints of APOBEC mutagenesis and is first discovered to be transactivated by FOSL1. Mechanistically, upregulated A3A exacerbates cytosolic double-stranded DNA (dsDNA) accumulation, thus stimulating cGAS-STING pathway. Simultaneously, A3A is associated with immunotherapy response which is predicted by TIDE algorithm, validated in a clinical cohort and further confirmed in mouse models. These findings systematically elucidate the clinical relevance, immunological characteristics, prognostic value for immunotherapy and underlying mechanisms of APOBEC mutagenesis in ESCC, which demonstrate great potential in clinical utility to facilitate clinical decisions.

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.

An immunogenic and oncogenic feature-based classification for chemotherapy plus PD-1 blockade in advanced esophageal squamous cell carcinoma

Although chemotherapy plus PD-1 blockade (chemo+anti-PD-1) has become the standard first-line therapy for advanced esophageal squamous cell carcinoma (ESCC), reliable biomarkers for this regimen are lacking. Here we perform whole-exome sequencing on tumor samples from 486 patients of the JUPITER-06 study and develop a copy number alteration-corrected tumor mutational burden that depicts immunogenicity more precisely and predicts chemo+anti-PD-1 efficacy. We identify several other favorable immunogenic features (e.g., HLA-I/II diversity) and risk oncogenic alterations (e.g., PIK3CA and TET2 mutation) associated with chemo+anti-PD-1 efficacy. An esophageal cancer genome-based immuno-oncology classification (EGIC) scheme incorporating these immunogenic features and oncogenic alterations is established. Chemo+anti-PD-1 achieves significant survival improvements in EGIC1 (immunogenic feature-favorable and oncogenic alteration-negative) and EGIC2 (either immunogenic feature-favorable or oncogenic alteration-negative) subgroups, but not the EGIC3 subgroup (immunogenic feature-unfavorable and oncogenic alteration-positive). Thus, EGIC may guide future individualized treatment strategies and inform mechanistic biomarker research for chemo+anti-PD-1 treatment in patients with advanced ESCC.

PSMD2 contributes to the progression of esophageal squamous cell carcinoma by repressing autophagy

BACKGROUND

The ubiquitin-proteasome and autophagy-lysosomal systems collaborate in regulating the levels of intracellular proteins. Dysregulation of protein homeostasis is a central feature of malignancy. The gene encoding 26S proteasome non-ATPase regulatory subunit 2 (PSMD2) of the ubiquitin-proteasome system is an oncogene in various types of cancer. However, the detailed role of PSMD2 in autophagy and its relationship to tumorigenesis in esophageal squamous cell carcinoma (ESCC) remain unknown. In the present study, we have investigated the tumor-promoting roles of PSMD2 in the context of autophagy in ESCC.

METHODS

Molecular approaches including DAPgreen staining, 5-Ethynyl-2'-deoxyuridine (EdU), cell counting kit 8 (CCK8), colony formation, transwell assays, and cell transfection, xenograft model, immunoblotting and Immunohistochemical analysis were used to investigate the roles of PSMD2 in ESCC cells. Data-independent acquisition (DIA) quantification proteomics analysis and rescue experiments were used to study the roles of PSMD2 in ESCC cells.

RESULTS

We demonstrate that the overexpression of PSMD2 promotes ESCC cell growth by inhibiting autophagy and is correlated with tumor progression and poor prognosis of ESCC patients. DIA quantification proteomics analysis shows a significant positive correlation between argininosuccinate synthase 1 (ASS1) and PSMD2 levels in ESCC tumors. Further studies indicate that PSMD2 activates the mTOR pathway by upregulating ASS1 to inhibit autophagy.

CONCLUSIONS

PSMD2 plays an important role in repressing autophagy in ESCC, and represents a promising biomarker to predict prognosis and a therapeutic target of ESCC patients.

Computational Methods Summarizing Mutational Patterns in Cancer: Promise and Limitations for Clinical Applications

Since the rise of next-generation sequencing technologies, the catalogue of mutations in cancer has been continuously expanding. To address the complexity of the cancer-genomic landscape and extract meaningful insights, numerous computational approaches have been developed over the last two decades. In this review, we survey the current leading computational methods to derive intricate mutational patterns in the context of clinical relevance. We begin with mutation signatures, explaining first how mutation signatures were developed and then examining the utility of studies using mutation signatures to correlate environmental effects on the cancer genome. Next, we examine current clinical research that employs mutation signatures and discuss the potential use cases and challenges of mutation signatures in clinical decision-making. We then examine computational studies developing tools to investigate complex patterns of mutations beyond the context of mutational signatures. We survey methods to identify cancer-driver genes, from single-driver studies to pathway and network analyses. In addition, we review methods inferring complex combinations of mutations for clinical tasks and using mutations integrated with multi-omics data to better predict cancer phenotypes. We examine the use of these tools for either discovery or prediction, including prediction of tumor origin, treatment outcomes, prognosis, and cancer typing. We further discuss the main limitations preventing widespread clinical integration of computational tools for the diagnosis and treatment of cancer. We end by proposing solutions to address these challenges using recent advances in machine learning.

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.

Multiancestry genomic and transcriptomic analysis of gastric cancer

Gastric cancer is among the most common malignancies worldwide, characterized by geographical, epidemiological and histological heterogeneity. Here, we report an extensive, multiancestral landscape of driver events in gastric cancer, involving 1,335 cases. Seventy-seven significantly mutated genes (SMGs) were identified, including ARHGAP5 and TRIM49C. We also identified subtype-specific drivers, including PIGR and SOX9, which were enriched in the diffuse subtype of the disease. SMGs also varied according to Epstein-Barr virus infection status and ancestry. Non-protein-truncating CDH1 mutations, which are characterized by in-frame splicing alterations, targeted localized extracellular domains and uniquely occurred in sporadic diffuse-type cases. In patients with gastric cancer with East Asian ancestry, our data suggested a link between alcohol consumption or metabolism and the development of RHOA mutations. Moreover, mutations with potential roles in immune evasion were identified. Overall, these data provide comprehensive insights into the molecular landscape of gastric cancer across various subtypes and ancestries.

M1 macrophage predicted efficacy of neoadjuvant camrelizumab combined with chemotherapy vs chemotherapy alone for locally advanced ESCC: A pilot study

Introduction

The efficacy and safety of immunotherapy have been widely recognized in gastrointestinal-related cancers. However, the efficacy of neoadjuvant camrelizumab for locally advanced esophageal squamous cell carcinoma (ESCC) has not been firmly established. This study compared the efficacy of camrelizumab in combination with neoadjuvant DCF (docetaxel, cisplatin and fluorouracil), with DCF alone for ESCC, and exploring biomarkers related to immune infiltration of the ESCC immunotherapy response.

Methods

We enrolled and randomly assigned patients with stage II-IVa ESCC to two study treatments: camrelizumab combined with docetaxel, cisplatin and fluorouracil (DCF) regimen and DCF regimen alone. The tissue for multiplex immunofluorescence (mIF) was obtained before and after neoadjuvant therapy. The Response Evaluation Criteria in Solid Tumors RECIST Version 1.1 (RECIST 1.1) and Tumor Regression Grade (TRG) was used to evaluate efficacy.

Results

A total of 30 patients were enrolled in the study. Following neoadjuvant camrelizumab, the objective response rate (ORR) and the disease control rate (DCR) were 46.7% (7/15) and 95.7% (14/15), respectively. No patients reported complete remission, while ORR and DCR in the chemotherapy group were 26.7% (4/15) and 86.7% (13/15), respectively. R0 resection after neoadjuvant treatment was achieved in 3 out of 15 patients in the combined group and in all patients (15/15) in the chemotherapy group. In the combined group, M1-type tumor-associated macrophages and CD56dim NK cells were more abundant in responders than in non-responders (p &lt; 0.05). A higher M1/M2 ratio was observed in responders (p &lt; 0.05). With respect to the NGS, among the copy number amplified genes, the 11q13 amplicon (CCND1/FGF19/FGF4/FGF3) showed the highest frequency (47%, 7/15).

Conclusions

Neoadjuvant camrelizumab combined with chemotherapy improved ORR in locally advanced ESCC. M1-type tumor-associated macrophages and CD56dim NK cells might be utilized to predict camrelizumab efficacy.

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.

The Molecular Characterization of Genetic Abnormalities in Esophageal Squamous Cell Carcinoma May Foster the Development of Targeted Therapies

Esophageal cancer is among the most common tumors in the world and is associated with poor outcomes, with a 5-year survival rate of about 10-20%. Two main histological subtypes are observed: esophageal squamous cell carcinoma (ESCC), more frequent among Asian populations, and esophageal adenocarcinoma (EAC), the predominant type in Western populations. The development of molecular analysis techniques has led to the definition of the molecular alterations observed in ESCC, consistently differing from those observed in EAC. The genetic alterations observed are complex and heterogeneous and involve gene mutations, gene deletions and gene amplifications. However, despite the consistent progress in the definition of the molecular basis of ESCC, precision oncology for these patients is still virtually absent. The recent identification of molecular subtypes of ESCC with clinical relevance may foster the development of new therapeutic strategies. It is estimated that about 40% of the genetic alterations observed in ESCC are actionable. Furthermore, the recent introduction of solid tumor immunotherapy with immune checkpoint inhibitors (ICIs) showed that a minority of ESCC patients are responsive, and the administration of ICIs, in combination with standard chemotherapy, significantly improves overall survival over chemotherapy in ESCC patients with advanced disease.

Potent molecular-targeted therapies for advanced esophageal squamous cell carcinoma

Esophageal cancer (EC) remains a public health concern with a high mortality and disease burden worldwide. Esophageal squamous cell carcinoma (ESCC) is a predominant histological subtype of EC that has unique etiology, molecular profiles, and clinicopathological features. Although systemic chemotherapy, including cytotoxic agents and immune checkpoint inhibitors, is the main therapeutic option for recurrent or metastatic ESCC patients, the clinical benefits are limited with poor prognosis. Personalized molecular-targeted therapies have been hampered due to the lack of robust treatment efficacy in clinical trials. Therefore, there is an urgent need to develop effective therapeutic strategies. In this review, we summarize the molecular profiles of ESCC based on the findings of pivotal comprehensive molecular analyses, highlighting potent therapeutic targets for establishing future precision medicine for ESCC patients, with the most recent results of clinical trials.

Esophageal Squamous Cancer from 4NQO-Induced Mice Model: CNV Alterations

Squamous esophageal carcinoma is a common pathological type of esophageal carcinoma around the world. The prognosis of esophageal carcinoma is usually poor and diagnosed at late stages. Recently, research suggested that genomic instability occurred in esophageal cells during the development of esophageal squamous cell carcinoma (ESCC). Identifying prognostic and specific genomic characteristics, especially at the early hyperplasia stage, is critical. Mice were given 4-nitroquinoline 1-oxide (4NQO) with drinking water to induce esophageal cancer. The immortalized human esophageal epithelial cell line (NE2) was also treated with 4NQO. We performed histologic analyses, immunofluorescence, and immunohistochemical staining to detect DNA damage at different time points. Whole-exome sequencing was accomplished on the esophagus tissues at different pathological stages to detect single-nucleotide variants and copy number variation (CNV) in the genome. Our findings indicate that all mice were tumor-forming, and a series of changes from simple hyperplasia (ESSH) to intraepithelial neoplasia (IEN) to esophageal squamous cell carcinoma (ESCC) was seen at different times. The expression of γ-H2AX increased from ESSH to ESCC. In addition, mutations of the Muc4 gene were detected throughout the pathological stages. Furthermore, CNV burden appeared in the esophageal tissues from the beginning of ESSH and accumulated more in cancer with the deepening of the lesions. This study demonstrates that mutations caused by the early appearance of DNA damage may appear in the early stage of malignant tissue before the emergence of atypia. The detection of CNV and mutations of the Muc4 gene may be used as an ultra-early screening indicator for esophageal cancer.

Effect of the Age-Adjusted Charlson Comorbidity Index on the Survival of Esophageal Squamous Cell Carcinoma Patients after Radical Esophagectomy

We aimed to investigate whether the age-adjusted Charlson comorbidity index (ACCI) can predict the postoperative overall survival (OS) and cancer-specific survival (CSS) of esophageal squamous cell carcinoma (ESCC) patients. Between 1 July 2015 and 31 July 2021, a retrospective cohort study was conducted among patients with primary ESCC who underwent radical esophagectomy. A total of 352 patients were included, with median age of 63.00 (IQR (interquartile range) 56.00-68.00). The patients were divided into low (n = 300) and high (n = 52) ACCI groups based on the optimal cut-off value of 5 points. Chronic pulmonary disease (38.4%) was the most common comorbidity. The results of the multivariate Cox regression showed that the ACCI (HR = 1.63, 95%CI: 1.04-2.56), tumor size (HR = 1.67, 95%CI: 1.05-2.66), pTNM (II vs. I, HR = 4.74, 95%CI: 1.82-12.32; III vs. I, HR = 6.08, 95%CI: 2.37-15.60), and postoperative chemotherapy (HR = 0.60, 95%CI: 0.40-0.91) were significantly associated with the OS. Furthermore, the ACCI, tumor size, pTNM, and postoperative chemotherapy were also significantly associated with the CSS. Interactions were identified between the ACCI and postoperative chemotherapy, pTNM stage, and tumor size in relation to the OS and CSS. In conclusion, the ACCI may be an independent prognostic factor affecting the long-term prognosis of patients after radical esophagectomy.

Identification of natural product 3, 5-diiodotyrosine as APOBEC3B inhibitor to prevent somatic mutation accumulation and cancer progression

BACKGROUND

The development of cancer is largely dependent on the accumulation of somatic mutations, indicating the potential to develop cancer chemoprevention agents targeting mutation drivers. However, ideal cancer chemoprevention agents that can effectively inhibit the mutation drivers have not been identified yet.

METHODS

The somatic mutation signatures and expression analyses of APOBEC3B were performed in patient with pan-cancer. The computer-aided screening and skeleton-based searching were performed to identify natural products that can inhibit the activity of APOBEC3B. 4-nitroquinoline-1-oxide (4-NQO)-induced spontaneous esophageal squamous cell carcinoma (ESCC) and azoxymethane/dextran sulfate sodium (AOM/DSS)-induced spontaneous colon cancer mouse models were conducted to investigate the influences of APOBEC3B inhibitor on the prevention of somatic mutation accumulation and cancer progression.

RESULTS

Here, we discovered that the cytidine deaminase APOBEC3B correlated somatic mutations were widely observed in a variety of cancers, and its overexpression indicated poor survival. SMC247 (3, 5-diiodotyrosine), as a source of kelp iodine without side effects, could strongly bind APOBEC3B (K_D=65 nM) and effectively inhibit its deaminase activity (IC₅₀=1.69 µM). Interestingly, 3, 5-diiodotyrosine could significantly reduce the clusters of mutations, prevent the precancerous lesion progression, and prolong the survival in 4-NQO-induced spontaneous ESCC and AOM/DSS-induced spontaneous colon cancer mouse models. Furthermore, 3, 5-diiodotyrosine could reduce colitis, increase the proportion and function of T lymphocytes via IL-15 in tumor microenvironment. The synergistic cancer prevention effects were observed when 3, 5-diiodotyrosine combined with PD-1/PD-L1 blockade.

CONCLUSIONS

This is the first prove-of-concept study to elucidate that the natural product 3, 5-diiodotyrosine could prevent somatic mutation accumulation and cancer progression through inhibiting the enzymatic activity of APOBEC3B. In addition, 3, 5-diiodotyrosine could reduce the colitis and increase the infiltration and function of T lymphocytes via IL-15 in tumor microenvironment. 3, 5-diiodotyrosine combined with PD-1/PD-L1 blockade could elicit synergistic cancer prevention effects, indicating a novel strategy for both prevent the somatic mutation accumulation and the immune-suppressive microenvironment exacerbation.

Prognostic value of circulating tumour DNA during post-radiotherapy surveillance in locally advanced esophageal squamous cell carcinoma

BACKGROUND

The potential of circulating tumour DNA (ctDNA) as a reliable biomarker for relapse/metastasis early detection and prognosis in esophageal squamous cell carcinoma (ESCC) after radiotherapy/chemoradiotherapy (RT/CRT) initiation requires comprehensive investigation.

METHODS

Treatment-naive locally advanced ESCC patients with available baseline plasma samples were prospectively enrolled from November 2018 to January 2020. RT/CRT was delivered with a simultaneous integrated boost of radiation dose. Serial plasma samples were collected at baseline (T0 ), week 4 of RT/CRT (T1 ), 1-3 (T2 ) and 3-6 months post-RT/CRT (T3 ). ctDNA was analysed using next-generation sequencing of 474 cancer-relevant genes.

RESULTS

A total of 128 plasma samples from 40 eligible patients were analysed (median age: 64 [range: 40-78], 88% males, 95% stage III/IV), and the median follow-up time was 20.6 months (range: 12.2-33.3). During the post-RT/CRT surveillance including 36 patients, radiological progression was observed in 16 patients, and 69% (11/16) had detectable post-RT/CRT ctDNA prior to radiological progression, with a median lead time of 4.4 months compared with radiological imaging. ctDNA positivity at T1 (hazard ratio, HR: 3.60, 95% confidence interval, CI: 1.30-10.01) or T2 (HR: 5.45, 95% CI: 1.72-17.26) indicated inferior progression-free survival (PFS). ctDNA clearance between T0 -T1 (HR: 0.31, 95% CI: 0.08-1.13) or T0 -T2 (HR: 0.11; 95% CI: 0.02-0.61) was associated with relatively favourable PFS. Similar results were obtained when focusing on patients without esophagectomy after RT/CRT. Notably, detectable ctDNA at T1 was a potential indicator of high local recurrence risks (HR: 4.43, 95% CI: 1.31-15.04).

CONCLUSIONS

ctDNA was identified as a robust biomarker for early detection of disease progression and post-RT/CRT prognosis stratification in ESCC. Detectable ctDNA at week 4 of RT/CRT might indicate higher local recurrence risks, implying the potential clinical utility of ctDNA tests in guiding post-RT/CRT treatments for locoregional control in ESCC.

Characterization of somatic structural variations in 528 Chinese individuals with Esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) demonstrates high genome instability. Here, we analyze 528 whole genomes to investigate structural variations' mechanisms and biological functions. SVs show multi-mode distributions in size, indicating distinct mutational processes. We develop a tool and define five types of complex rearrangements with templated insertions. We highlight a type of fold-back inversion, which is associated with poor outcomes. Distinct rearrangement signatures demonstrate variable genomic metrics such as replicating time, spatial proximity, and chromatin accessibility. Specifically, fold-back inversion tends to occur near the centrosome; TD-c2 (Tandem duplication-cluster2) is significantly enriched in chromatin-accessibility and early-replication region compared to other signatures. Analyses of TD-c2 signature reveal 9 TD hotspots, of which we identify a hotspot consisting of a super-enhancer of PTHLH. We confirm the oncogenic effect of the PTHLH gene and its interaction with enhancers through functional experiments. Finally, extrachromosomal circular DNAs (ecDNAs) are present in 14% of ESCCs and have strong selective advantages to driver genes.

Single-Cell Transcriptomics of Endothelial Cells in Upper and Lower Human Esophageal Squamous Cell Carcinoma

Esophageal squamous cell carcinoma (ESCC) is a type of progressive and distant metastatic tumor. Targeting anti-angiogenic genes could effectively hinder ESCC development and metastasis, whereas ESCC locating on the upper or the lower esophagus showed different response to the same clinical treatment, suggesting ESCC location should be taken into account when exploring new therapeutic targets. In the current study, to find novel anti-angiogenic therapeutic targets, we identified endothelial cell subsets in upper and lower human ESCC using single-cell RNA sequencing (scRNA-seq), screened differentially expressed genes (DEGs), and performed gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The results showed that common DEGs shared in the upper and the lower endothelial cells mainly are involved in vessel development, angiogenesis, and cell motility of endothelial cells by regulating PI3K-AKT, Rap1, Ras, TGF-beta, and Apelin signaling pathways. The critical regulatory genes were identified as ITGB1, Col4A1, Col4A2, ITGA6, LAMA4, LAMB1, LAMC1, VWF, ITGA5, THBS1, PDGFB, PGF, RHOC, and CTNNB1. Cell metabolism-relevant genes, e.g., MGST3, PNP, UPP1, and HYAL2 might be the prospective therapeutic targets. Furthermore, we found that DEGs only in the upper endothelial cells, such as MAPK3, STAT3, RHOA, MAPK11, HIF1A, FGFR1, GNG5, GNB1, and ARHGEF12, mainly regulated cell adhesion, structure morphogenesis, and motility through Phospholipase D, Apelin, and VEGF signaling pathways. Moreover, DEGs only in the lower endothelial cells, for instance PLCG2, EFNA1, CALM1, and RALA, mainly regulated cell apoptosis and survival by targeting calcium ion transport through Rap1, Ras, cAMP, Phospholipase D, and Phosphatidylinositol signaling pathways. In addition, the upper endothelial cells showed significant functional diversity such as cytokine-responsive, migratory, and proliferative capacity, presenting a better angiogenic capacity and making it more sensitive to anti-angiogenic therapy compared with the lower endothelial cells. Our study has identified the potential targeted genes for anti-angiogenic therapy for both upper and lower ESCC, and further indicated that anti-angiogenic therapy might be more effective for upper ESCC, which still need to be further examined in the future.