Whole genome sequencing analysis identifies recurrent structural alterations in esophageal squamous cell carcinoma

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.

Immuno-genomic analysis reveals eosinophilic feature and favorable prognosis of female non-smoking esophageal squamous cell carcinomas

Most of esophageal squamous cell carcinoma (ESCC) develop in smoking males in Japan, but the genomic etiology and immunological characteristics of rare non-smoking female ECSS remain unclear. To elucidate the genomic and immunological features of ESCC in non-smoking females, we analyzed whole-genome or transcriptome sequencing data from 94 ESCCs, including 20 rare non-smoking female cases. In addition, 31,611 immune cells were extracted from four ESCC tissues and subject to single-cell RNA-seq. We compared their immuno-genomic and microbiome profiles between non-smoking female and smoking ESCCs. Non-smoking females showed much better prognosis. Whole-genome sequencing analysis showed no significant differences in driver genes or copy number alterations depending on smoking status. The mutational signatures specifically observed in non-smoking females ESCC could be attributed to aging. Immune profiling from RNA-seq revealed that ESCC in non-smoking females had high tumor microenvironment signatures and a high abundance of eosinophils with a favorable prognosis. Single-cell RNA-sequencing of intratumor immune cells revealed gender differences of eosinophils and their activation in female cases. ESCCs in non-smoking females have age-related mutational signatures and gender-specific tumor immune environment with eosinophils, which is likely to contribute to their favorable prognosis.

Comparing Genomic Landscapes of Oral and Cutaneous Squamous Cell Carcinoma of the Head and Neck: Quest for Novel Diagnostic Markers

Squamous cell carcinoma is the most common head and neck malignancy arising from the oral mucosa and the skin. The histologic and immunohistochemical features of oral squamous cell carcinoma (OSCC) and head and neck cutaneous squamous cell carcinoma (HNcSCC) are similar, making it difficult to identify the primary site in cases of metastases. With the advent of immunotherapy, reliable distinction of OSCC and HNcSCC at metastatic sites has important treatment and prognostic implications. Here, we investigate and compare the genomic landscape of OSCC and HNcSCC to identify diagnostically useful biomarkers. Whole-genome sequencing data from 57 OSCC and 41 HNcSCC patients were obtained for tumor and matched normal samples. Tumor mutation burden (TMB), Catalogue of Somatic Mutations in Cancer (COSMIC) mutational signatures, frequent chromosomal alterations, somatic single nucleotide, and copy number variations were analyzed. The median TMB of 3.75 in primary OSCC was significantly lower (P < .001) than that of 147.51 mutations/Mb in primary HNcSCC. The COSMIC mutation signatures were significantly different (P < .001) between OSCC and HNcSCC. OSCC showed COSMIC single-base substitution (SBS) mutation signature 1 and AID/APOBEC activity-associated signature 2 and/or 13. All except 1 HNcSCC from hair-bearing scalp showed UV damage-associated COSMIC SBS mutation signature 7. Both OSCC and HNcSCC demonstrated a predominance of tumor suppressor gene mutations, predominantly TP53. The most frequently mutated oncogenes were PIK3CA and MUC4 in OSCC and HNcSCC, respectively. The metastases of OSCC and HNcSCC demonstrated TMB and COSMIC SBS mutation signatures similar to their primary counterparts. The combination of high TMB and UV signature in a metastatic keratinizing squamous cell carcinoma suggests HNcSCC as the primary site and may also facilitate decisions regarding immunotherapy. HNcSCC and OSCC show distinct genomic profiles despite histologic and immunohistochemical similarities. Their genomic characteristics may underlie differences in behavior and guide treatment decisions in recurrent and metastatic settings.

Modeling tissue-specific breakpoint proximity of structural variations from whole-genomes to identify cancer drivers

Structural variations (SVs) in cancer cells often impact large genomic regions with functional consequences. However, identification of SVs under positive selection is a challenging task because little is known about the genomic features related to the background breakpoint distribution in different cancers. We report a method that uses a generalized additive model to investigate the breakpoint proximity curves from 2,382 whole-genomes of 32 cancer types. We find that a multivariate model, which includes linear and nonlinear partial contributions of various tissue-specific features and their interaction terms, can explain up to 57% of the observed deviance of breakpoint proximity. In particular, three-dimensional genomic features such as topologically associating domains (TADs), TAD-boundaries and their interaction with other features show significant contributions. The model is validated by identification of known cancer genes and revealed putative drivers in cancers different than those with previous evidence of positive selection.

Identifying structural variants (SVs) under positive selection in cancer is challenging. Here, the authors develop CSVDriver, a method that computes SV breakpoint proximity and the contribution of elements such as topologically associating domains, and identifies loci that show signs of positive selection and contain known and putative drivers.

Immuno-genomic profiling of biopsy specimens predicts neoadjuvant chemotherapy response in esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is one of the most aggressive cancers and is primarily treated with platinum-based neoadjuvant chemotherapy (NAC). Some ESCCs respond well to NAC. However, biomarkers to predict NAC sensitivity and their response mechanism in ESCC remain unclear. We perform whole-genome sequencing and RNA sequencing analysis of 141 ESCC biopsy specimens before NAC treatment to generate a machine-learning-based diagnostic model to predict NAC reactivity in ESCC and analyzed the association between immunogenomic features and NAC response. Neutrophil infiltration may play an important role in ESCC response to NAC. We also demonstrate that specific copy-number alterations and copy-number signatures in the ESCC genome are significantly associated with NAC response. The interactions between the tumor genome and immune features of ESCC are likely to be a good indicator of therapeutic capability and a therapeutic target for ESCC, and machine learning prediction for NAC response is useful.

•

Four different immune subtypes from RNA-seq of ESCC biopsy specimen

•

Neutrophils within tumors are associated with tumor sensitivity to NAC

•

Specific copy-number changes and signatures in ESCC are associated with NAC response

•

Machine learning prediction for NAC response using immunogenomics of ESCC is useful

CSMD1 suppresses cancer progression by inhibiting proliferation, epithelial-mesenchymal transition, chemotherapy-resistance and inducing immunosuppression in esophageal squamous cell carcinoma

Human CUB and Sushi multiple domains (CSMD1) is considered a crucial role in cancer progression, but the specific function in esophageal squamous cell carcinoma (ESCC) is not clear. Understanding the role of CSMD1 in ESCC progression may lead to a novel strategy for ESCC treatment. Here, we found that both CSMD1 mRNA and protein levels were downregulated in ESCC tissues. Reduced CSMD1 expression was correlated with a poor prognosis in ESCC patients. CSMD1 expression inhibited proliferation, migration and invasion in ESCC cell lines in vitro. CSMD1 deficiency in established xenografted tumors increases tumor size and weight. We further found that CSMD1-overexpression cells are more sensitive to chemotherapy. Moreover, we addressed the role of CSMD1 in the CD8⁺ T cell immune response. An in vitro killing assay showed that the cytotoxicity of CD8⁺ T cells was inhibited in CSMD1-overexpression tumor cells. In vivo, in CSMD1 deficiency tumor-bearing mice activation and expansion of CD8⁺ T cells were increased. Further investigation showed that CSMD1 expression on tumor cells was positively correlated with CD8⁺ T cells infiltration and cytokines secretion. These findings highlight that CSMD1 is a tumor suppressor gene in ESCC patients and a positive regulator of CD8⁺ T cells expansion and activation, and could increase cytokines secretion, indicating that tumor cell-associated CSMD1 might be a target for ESCC.

Integrative Genomic Analyses of 1,145 Patient Samples Reveal New Biomarkers in Esophageal Squamous Cell Carcinoma

Due to the lack of effective diagnostic markers and therapeutic targets, esophageal squamous cell carcinoma (ESCC) shows a poor 5 years survival rate of less than 30%. To explore the potential therapeutic targets of ESCC, we integrated and reanalyzed the mutation data of WGS (whole genome sequencing) or WES (whole exome sequencing) from a total of 1,145 samples in 7 large ESCC cohorts, including 270 ESCC gene expression data. Two new mutation signatures and 20 driver genes were identified in our study. Among them, AP3S1, MUC16, and RPS15 were reported for the first time. We also discovered that the KMT2D was associated with the multiple clinical characteristics of ESCC, and KMT2D knockdown cells showed enhanced cell migration and cell invasion. Furthermore, a few neoantigens were shared between ESCC patients. For ESCC, compared to TMB, neoantigen might be treated as a better immunotherapy biomarker. Our research expands the understanding of ESCC mutations and helps the identification of ESCC biomarkers, especially for immunotherapy biomarkers.

Chromothripsis is a frequent event and underlies typical genetic changes in early T-cell precursor lymphoblastic leukemia in adults

Chromothripsis is a mitotic catastrophe that arises from multiple double strand breaks and incorrect re-joining of one or a few chromosomes. We report on incidence, distribution, and features of chromothriptic events in T-cell acute lymphoblastic leukemias (T-ALL). SNP array was performed in 103 T-ALL (39 ETP/near ETP, 59 non-ETP, and 5 with unknown stage of differentiation), including 38 children and 65 adults. Chromothripsis was detected in 11.6% of all T-ALL and occurred only in adult cases with an immature phenotype (12/39 cases; 30%). It affected 1 to 4 chromosomes, and recurrently involved chromosomes 1, 6, 7, and 17. Abnormalities of genes typically associated with T-ALL were found at breakpoints of chromothripsis. In addition, it gave rise to new/rare alterations, such as, the SFPQ::ZFP36L2 fusion, reported in pediatric T-ALL, deletions of putative suppressors, such as IKZF2 and CSMD1, and amplification of the BCL2 gene. Compared to negative cases, chromothripsis positive T-ALL had a significantly higher level of MYCN expression, and a significant downregulation of RGCC, which is typically induced by TP53 in response to DNA damage. Furthermore we identified mutations and/or deletions of DNA repair/genome stability genes in all cases, and an association with NUP214 rearrangements in 33% of cases.

LncRNA FOXP4-AS1 Promotes the Progression of Esophageal Squamous Cell Carcinoma by Interacting With MLL2/H3K4me3 to Upregulate FOXP4

Malignant tumors are a grave threat to human health. Esophageal squamous cell carcinoma (ESCC) is a common gastrointestinal malignant tumor. China has a high incidence of ESCC, and its morbidity and mortality are higher than the global average. Increasingly, studies have shown that long noncoding RNAs (lncRNAs) play a vital function in the occurrence and development of tumors. Although the biological function of FOXP4-AS1 has been demonstrated in various tumors, the potential molecular mechanism of FOXP4-AS1 in ESCC is still poorly understood. The expression of FOXP4 and FOXP4-AS1 was detected in ESCC by quantitative real-time PCR (qRT–PCR) or SP immunohistochemistry (IHC). shRNA was used to silence gene expression. Apoptosis, cell cycle, MTS, colony formation, invasion and migration assays were employed to explore the biological functions of FOXP4 and FOXP4-AS1. The potential molecular mechanism of FOXP4-AS1 in ESCC was determined by dual-luciferase reporter, RNA immunoprecipitation (RIP) and chromatin immunoprecipitation (ChIP). Here, we demonstrated that FOXP4-AS1 was significantly increased in ESCC tissues and cell lines, associated with lymph node metastasis and TNM staging. Cell function experiments showed that FOXP4-AS1 promoted the proliferation, invasion and migration ability of ESCC cells. The expression of FOXP4-AS1 and FOXP4 in ESCC tissues was positively correlated. Further research found that FOXP4-AS1, upregulated in ESCC, promotes FOXP4 expression by enriching MLL2 and H3K4me3 in the FOXP4 promoter through a “molecular scaffold”. Moreover, FOXP4, a transcription factor of β-catenin, promotes the transcription of β-catenin and ultimately leads to the malignant progression of ESCC. Finally, FOXP4-AS1 may be a new therapeutic target for ESCC.

Liquid Biopsy to Detect Minimal Residual Disease: Methodology and Impact

One reason why some patients experience recurrent disease after a curative-intent treatment might be the persistence of residual tumor cells, called minimal residual disease (MRD). MRD cannot be identified by standard radiological exams or clinical evaluation. Tumor-specific alterations found in the blood indirectly diagnose the presence of MRD. Liquid biopsies thus have the potential to detect MRD, allowing, among other things, the detection of circulating tumor DNA (ctDNA), circulating tumor cells (CTC), or tumor-specific microRNA. Although liquid biopsy is increasingly studied, several technical issues still limit its clinical applicability: low sensitivity, poor standardization or reproducibility, and lack of randomized trials demonstrating its clinical benefit. Being able to detect MRD could give clinicians a more comprehensive view of the risk of relapse of their patients and could select patients requiring treatment escalation with the goal of improving cancer survival. In this review, we are discussing the different methodologies used and investigated to detect MRD in solid cancers, their respective potentials and issues, and the clinical impacts that MRD detection will have on the management of cancer patients.

Integrative genome-wide analyses reveal the transcriptional aberrations in Japanese esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is a malignant disease. At present, the genomic profiles of ESCC are known to a considerable extent, and DNA methylation and gene expression profiles have been mainly used for the classification of ESCC subtypes, but integrative genomic, transcriptomic, and epigenomic analyses remain insufficient. Therefore, we performed integrative analyses using whole-exome sequencing, DNA methylation, and RNA sequencing (RNA-seq) analyses of Japanese patients with ESCC. In cancer-related genes, such as NOTCH family genes, RTK/PI3K pathway genes, and NFE2L2 pathway genes, variants and copy number amplification were detected frequently. Japanese ESCC cases were clustered into two mutational signatures: an APOBEC-associated signature and an age-related signature. In imprinted genes, DNA methylation was aberrant in gene promoter regions and correlated well with gene expression profiles. Nonsynonymous single-nucleotide variants and allelic expression imbalance were detected frequently in FAT family genes. Our integrative genome-wide analyses, including DNA methylation and allele-specific gene expression profiles, revealed altered gene regulation of imprinted genes and FAT family genes in ESCC.

Amplification of 8p11.23 in cancers and the role of amplicon genes

Copy number alterations are widespread in cancer genomes and are part of the genomic instability underlying the pathogenesis of neoplastic diseases. Recurrent copy number alterations of specific chromosomal loci may result in gains of oncogenes or losses of tumor suppressor genes and become entrenched in the genomic framework of certain types of cancers. The locus at chromosome 8p11.23 presents recurrent amplifications most commonly in squamous lung carcinomas, breast cancers, squamous esophageal carcinomas, and urothelial carcinomas. Amplification is rare in other cancers. The amplified segment involves several described oncogenes that may promote cancer cell survival and proliferation, as well as less well characterized genes that could also contribute to neoplastic processes. Genes proposed to be "drivers" in 8p11.23 amplifications include ZNF703, FGFR1 and PLPP5. Additional genes in the locus that could be functionally important in neoplastic networks include co-chaperone BAG4, lysine methyltransferase NSD3, ASH2L, a member of another methyltransferase complex, MLL and the mRNA processing and translation regulators LSM1 and EIF4EBP1. In this paper, genes located in the amplified segment of 8p11.23 will be examined for their role in cancer and data arguing for their importance for cancers with the amplification will be presented.