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

Integrative genomics unveils basement membrane-related diagnostic markers and therapeutic targets in esophageal squamous cell carcinoma

BACKGROUND

Esophageal squamous cell carcinoma (ESCC) is often diagnosed at advanced stages due to the inherent limitations of current screening methodologies. Central to evaluating tumor invasion and prognostic assessment in ESCC is the integrity of the basement membrane (BM). However, current research on the implications of BM-related genes (BMRGs) in diagnosing ESCC remains sparse.

METHODS

We performed a comprehensive analysis using single-cell RNA-sequencing (scRNA-seq) data from the Gene Expression Omnibus (GEO) database, alongside gene expression profiles acquired from GEO and The Cancer Genome Atlas (TCGA) databases. This identified differentially expressed BMRGs in ESCC. Employing LASSO, RF, and SVM-RFE, we selected potential BM biomarkers and crafted a diagnostic nomogram for ESCC, validated by ROC curves and AUC values. We also explored immune infiltration and biological mechanisms through consensus clustering and GSVA, and utilized single cell trajectory analysis and GSCALite to study gene distributions and pathways. In vitro experiments further elucidated the role of these genes in ESCC carcinogenesis.

RESULTS

Here, we discovered that ESCC cell types exhibited markedly elevated BM-related scores. Our analysis pinpointed seven BM genes upregulated and linked to immune infiltration, showcasing unique gene expression profiles and varying immune cell densities across the BM-related subtypes. Furthermore, a robust positive correlation was observed between these genes expression and EMT activity. The knockdown of BGN significantly suppressed cell proliferation, migration, invasion, while also augmenting cell viability following chemotherapy drug treatment.

CONCLUSION

Our study identified seven key BMRGs (BGN, LAMB3, SPARC, MMP1, LUM, COL4A1, and NELL2) and established a diagnostic nomogram for ESCC. Of noteworthy significance is the discovery of BGN as a promising drug target, indicating a novel strategy for future clinical combination therapies in ESCC.

Molecular patterns and mechanisms of tumorigenesis in HPV-associated and HPV-independent sinonasal squamous cell carcinoma

Mechanisms of tumorigenesis in sinonasal squamous cell carcinoma (SNSCC) remain poorly described due to its rare nature. A subset of SNSCC are associated with the human papillomavirus (HPV); however, it is unknown whether HPV is a driver of HPV-associated SNSCC tumorigenesis or merely a neutral bystander. We hypothesized that performing the first large high-throughput sequencing study of SNSCC would reveal molecular mechanisms of tumorigenesis driving HPV-associated and HPV-independent SNSCC and identify targetable pathways. High-throughput sequencing was performed on 64 patients with HPV-associated and HPV-independent sinonasal carcinomas. Mutation annotation, viral integration, copy number, and pathway-based analyses were performed. Analysis of HPV-associated SNSCC revealed similar mutational patterns observed in HPV-associated cervical and head and neck squamous cell carcinoma, including lack of TP53 mutations and the presence of known hotspot mutations in PI3K and FGFR3. Further similarities included enrichment of APOBEC mutational signature, viral integration at known hotspot locations, and frequent mutations in epigenetic regulators. HPV-associated SNSCC-specific recurrent mutations were also identified including KMT2C , UBXN11 , AP3S1 , MT-ND4 , and MT-ND5 . Mutations in KMT2D and FGFR3 were associated with decreased overall survival. We developed the first known HPV-associated SNSCC cell line and combinatorial small molecule inhibition of YAP/TAZ and PI3K pathways synergistically inhibited tumor cell clonogenicity. In conclusion, HPV-associated SNSCC and HPV-independent SNSCC are driven by molecularly distinct mechanisms of tumorigenesis. Combinatorial blockade of YAP/TAZ and vertical inhibition of the PI3K pathway may be useful in targeting HPV-associated SNSCC whereas targeting MYC and horizontal inhibition of RAS/PI3K pathways for HPV-independent SNSCC.

One Sentence Summary

This study solidifies HPV as a driver of HPV-associated SNSCC tumorigenesis, identifies molecular mechanisms distinguishing HPV-associated and HPV-independent SNSCC, and elucidates YAP/TAZ and PI3K blockade as key targets for HPV-associated SNSCC.

Predictive role of ctDNA in esophageal squamous cell carcinoma receiving definitive chemoradiotherapy combined with toripalimab

The combination of toripalimab (an anti-PD-1 antibody) with definitive chemoradiotherapy (CRT) demonstrated encouraging efficacy against locally advanced esophageal squamous cell carcinoma (ESCC) in the EC-CRT-001 phase II trial (NCT04005170). The primary endpoint of this trial was the clinical complete response rate (cCR), and the secondary endpoints included overall survival (OS), progression-free survival (PFS), duration of response, and quality of life. The exploratory analyses of EC-CRT-001 include exploring the role of circulating tumor DNA (ctDNA) and blood-based tumor mutational burden (bTMB) in predicting the response and survival. In total, 118 blood and 35 tissue samples from 42 enrolled patients were included in the analyses. We found that ctDNA-negative patients achieved a higher cCR compared to those with detectable ctDNA during CRT (83%, 19/23 vs. 39%, 7/18; p = 0.008) or post-CRT (78%, 21/27 vs. 30%, 3/10; p = 0.017). Patients with detectable ctDNA during CRT had shorter PFS (p = 0.014). Similarly, patients with post-CRT detectable ctDNA had a significantly shorter PFS (p = 0.012) and worse OS (p = 0.004). Moreover, patients with high bTMB levels during CRT had prolonged OS (p = 0.027). In conclusion, ctDNA and bTMB have the potential to predict treatment efficacy and survival in ESCC treated with CRT and immunotherapy.

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.

Differences in genome, transcriptome, miRNAome, and methylome in synchronous and metachronous liver metastasis of colorectal cancer

Despite distant metastases being the critical factor affecting patients' survival, they remain poorly understood. Our study thus aimed to molecularly characterize colorectal cancer liver metastases (CRCLMs) and explore whether molecular profiles differ between Synchronous (SmCRC) and Metachronous (MmCRC) colorectal cancer. This characterization was performed by whole exome sequencing, whole transcriptome, whole methylome, and miRNAome. The most frequent somatic mutations were in APC, SYNE1, TP53, and TTN genes. Among the differently methylated and expressed genes were those involved in cell adhesion, extracellular matrix organization and degradation, neuroactive ligand-receptor interaction. The top up-regulated microRNAs were hsa-miR-135b-3p and -5p, and the hsa-miR-200-family while the hsa-miR-548-family belonged to the top down-regulated. MmCRC patients evinced higher tumor mutational burden, a wider median of duplications and deletions, and a heterogeneous mutational signature than SmCRC. Regarding chronicity, a significant down-regulation of SMOC2 and PPP1R9A genes in SmCRC compared to MmCRC was observed. Two miRNAs were deregulated between SmCRC and MmCRC, hsa-miR-625-3p and has-miR-1269-3p. The combined data identified the IPO5 gene. Regardless of miRNA expression levels, the combined analysis resulted in 107 deregulated genes related to relaxin, estrogen, PI3K-Akt, WNT signaling pathways, and intracellular second messenger signaling. The intersection between our and validation sets confirmed the validity of our results. We have identified genes and pathways that may be considered as actionable targets in CRCLMs. Our data also provide a valuable resource for understanding molecular distinctions between SmCRC and MmCRC. They have the potential to enhance the diagnosis, prognostication, and management of CRCLMs by a molecularly targeted approach.

AP3S1 is a Novel Prognostic Biomarker and Correlated With an Immunosuppressive Tumor Microenvironment in Pan-Cancer

Background: Adaptor-related protein complex 3, sigma one subunit (AP3S1) is one of the encoding subunits of the adaptor complex AP-3. However, its role in various tumor types and relationship with the tumor immune microenvironment (TIME) remains unclear.

Methods: AP3S1 expression was analyzed using datasets from The Cancer Genome Atlas, Genotype-Tissue Expression, UALCAN, and HPA databases. Then, we performed a systematic analysis of the genetic alterations, clinical features, and prognostic value of AP3S1 in pan-cancer. Gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA) were used to identify the signaling pathways associated with AP3S1. The correlation between immune cell infiltration and AP3S1 expression was analyzed using immune cell infiltration data from the ImmuCellAI, TIMER2, and a previous study. Finally, we analyzed the association of AP3S1 with tumor mutational burden (TMB), microsatellite instability (MSI), and immune-related genes.

Results: We found AP3S1 overexpression in most tumors and a significant association with low survival rates. GSEA and GSVA results show that AP3S1 is involved in tumor progression and associated with immune pathways in different tumor types. We also found that AP3S1 expression was positively correlated with the level of infiltration of immunosuppressive cells (tumor-associated macrophages, cancer-associated fibroblasts, Tregs) and negatively correlated with immune killer cells, including NK cells and CD8+ T cells, in pan-cancer. The expression of AP3S1 could affect TMB and MSI in various cancers. In addition, AP3S1 was positively correlated with most immunosuppressive genes, including PD-1, PD-L1, CTLA4, LAG3 and TIGIT in most cancer types.

Conclusion: Our study reveals that AP3S1 is a potential pan-cancer oncogene and plays an essential role in tumorigenesis and cancer immunity. Elevated expression of AP3S1 indicates an immunosuppressive microenvironment and can be used as a potential prognostic biomarker and a target for immunotherapy.