Esophageal squamous cell carcinoma (ESCC) is highly heterogeneous. Our understanding of full molecular and immune landscape of ESCC remains limited, hindering the development of personalised therapeutic strategies. To address this, we perform genomic-transcriptomic characterizations and AI-aided histopathological image analysis of 120 Chinese ESCC patients. Here we show that ESCC can be categorized into differentiated, metabolic, immunogenic and stemness subtypes based on bulk and single-cell RNA-seq, each exhibiting specific molecular and histopathological features based on an amalgamated deep-learning model. The stemness subgroup with signature genes, such as WFDC2, SFRP1, LGR6 and VWA2, has the poorest prognosis and is associated with downregulated immune activities, a high frequency of EP300 mutation/activation, functional mutation enrichment in Wnt signalling and the highest level of intratumoural heterogeneity. The immune profiling by transcriptomics and immunohistochemistry reveals ESCC cells overexpress natural killer cell markers XCL1 and CD160 as immune evasion. Strikingly, XCL1 expression also affects the sensitivity of ESCC cells to common chemotherapy drugs. This study opens avenues for ESCC treatment and provides a valuable public resource to better understand ESCC.

Oridonin, a diterpenoid compound isolated from traditional Chinese medicine Rabdosia rubescens, has shown antitumor effects to esophageal cancer. However, its molecular mechanism is not fully understood, which limits its clinical application. In the present study, we used RNA-seq analysis to check the transcriptome changes after oridonin treatment and we found genes controlling the GSH-ROS system were up-regulated, namely SLC7A11, TXNRD1, TRIM16, SRXN1, GCLM, and GCLC. Furthermore, our data suggest that oridonin significantly increased the production of ROS in EC109 and TE1 cells, which can be inhibited by NAC. Interestingly, oridonin can dramatically reduce intracellular GSH levels in TE1 cells in a concentration and time-dependent manner. In addition, cell death caused by oridonin was strongly inhibited by GSH (1 mM), while GSSG (1 mM) had little effect. At the same time, we also found that oridonin showed selective cytotoxicity to esophageal squamous carcinoma cell with p53 mutation since mut-p53 cells had lower SLC7A11 expression, a component of the cystine/glutamate antiporter. We also found that γ-glutamyl cysteine synthetase inhibitor (BSO) synergizes with oridonin to strongly inhibit EC109 cells at a low dose. These results suggested that the antitumor effects of oridonin are based on its –SH reactivity and glutathione depletion. Esophageal squamous carcinoma cells with p53-mutation showed hypersensitivity to oridonin because of the suppression of SLC7A11 expression by p53 mutation.

Esophageal Squamous Cell carcinomas (ESCC) is a highly heterogeneous malignancy that is among the leading cause of cancer-related death worldwide. B cells play pivotal roles in the immune defense system and cancer progression and regression, yet the repertoire of tumor infiltrating B cells (TIBs) and its association with clinical outcome remains unexplored in ESCC. Here we collected bulk RNA-seq sequencing data from 119 ESCC tumors and matched adjacent normal samples to delineate the B cell repertoire. We found that ESCC is more heavily infiltrated by B cells and plasma cells compared to activated T cells. The immunoglobulin heavy chain variable region (IGHV) gene usage was remarkably biased and IGHV3-74 was under-represented in ESCC tumors. The TIBs showed a more oligoclonal profile along with widespread clonal expansion and IgG subclass switch events (CSRs). Survival analysis revealed several unexpected associations between tumor infiltrating B cells and prognosis. Higher levels of immunoglobulin expression (IGH), CD138 expression, IGH to MS4A1 ratio, CSR events and clone diversity are all associated with better survival. Notably, we found that the abundance of CD20-negative IgG2-producing plasma cells has a strong positive effect on overall survival with a hazard ratio (HR) of 0.40 (log-rank p: 0.002). Combing molecular subtyping, the IgG2-producing plasma cells could stratify high-risk patients more accurately with a HR of 0.253 (log-rank p: 0.0006). The direct link between protective B cell populations and ESCC prognosis provides biomarkers for high-risk patient selection and holds great promise for developing strategies for immunotherapy targeting B cells in ESCC patients.

In this study, we determined whether different subtypes of epidermal growth factor receptor (EGFR) exon19 mutation are associated with the therapeutic effect of EGFR-tyrosine kinase inhibitors (TKIs) on advanced non-small cell lung adenocarcinoma. A total of 122 patients with stage III or IV non-small cell lung adenocarcinoma were retrospectively reviewed. Clinical characteristics of these patients, including progression-free survival (PFS) outcome for EGFR-TKI treatment, were analyzed. According to the mutation pattern, we classified the in-frame deletions detected on EGFR Exon19 into three different types: codon deletion (CD), with a deletion of one or more original codons; codon substitution and skipping (CSS), with a deletion of one or two nucleotides but the residues could be translated into a new amino acid without changing following amino acid sequence; CD or CSS plus single nucleotide variant (SNV) (CD/CSS+SNV), exclude CD or CSS, there’s another SNV nearby the deletion region. The clinical characteristics of three groups were analyzed and as a result, no significant difference was found. By comparing the average number of missing bases and amino acids of the three mutation subtypes, it could be discovered that the number of missing bases and amino acids of the three mutation subtypes is diverse, and group CSS> group CD> group CD/CSS+SNV. Finally, survival analysis was performed between three groups of patients. The median PFS of group CD, group CSS and group CD/CSS+SNV was 11 months, 9 months and 14 months respectively. There was a distinct difference in the PFS between group CSS and group CD/CSS+SNV (P = 0.035<0.05), and the PFS of group CD/CSS+SNV was longer. Different mutation subtypes of EGFR exon19 can predict the therapeutic effect of EGFR-TKIs on advanced non-small cell lung adenocarcinoma.