Transcriptomics based multi-dimensional characterization and drug screen in esophageal squamous cell carcinoma

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.

Genetic and molecular characterization of metabolic pathway-based clusters in esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is one of the most aggressive types of squamous cell carcinoma and represents a significant proportion of esophageal cancer. Metabolic reprogramming plays a key role in the occurrence and development of ESCC. Unsupervised clustering analysis was employed to stratify ESCC samples into three clusters: MPC1-lipid type, MPC2-amino acid type, and MPC3-energy type, based on the enrichment scores of metabolic pathways extracted from the Reactome database. The MPC3 cluster exhibited characteristics of energy metabolism, with heightened glycolysis, cofactors, and nucleotide metabolism, showing a trend toward increased aggressiveness and poorer survival rates. On the other hand, MPC1 and MPC2 primarily involved lipid and amino acid metabolism, respectively. In addition, liquid chromatography‒mass spectrometry-based metabolite profiles and potential therapeutic agents were explored and compared among ESCC cell lines with different MPCs. MPC3 amplified energy metabolism markers, especially carnitines. In contrast, MPC1 and MPC2 predominantly had elevated levels of lipids (primarily triacylglycerol) and amino acids, respectively. Furthermore, MPC3 demonstrated a suboptimal clinical response to PD-L1 immunotherapy but showed increased sensitivity to the doramapimod chemotherapy regimen, as evident from drug sensitivity evaluations. These insights pave the way for a more personalized therapeutic approach, potentially enhancing treatment precision for ESCC patients.

UPP1 promotes lung adenocarcinoma progression through the induction of an immunosuppressive microenvironment

The complexity of the tumor microenvironment (TME) is a crucial factor in lung adenocarcinoma (LUAD) progression. To gain deeper insights into molecular mechanisms of LUAD, we perform an integrative single-cell RNA sequencing (scRNA-seq) data analysis of 377,574 cells from 117 LUAD patient samples. By linking scRNA-seq data with bulk gene expression data, we identify a cluster of prognostic-related UPP1high tumor cells. These cells, primarily situated at the invasive front of tumors, display a stronger association with the immunosuppressive components in the TME. Our cytokine array analysis reveals that the upregulation of UPP1 in tumor cells leads to the increased release of various immunosuppressive cytokines, with TGF-β1 being particularly prominent. Furthermore, this UPP1 upregulation also elevates the expression of PD-L1 through the PI3K/AKT/mTOR pathway, which contributes to the suppression of CD8 + T cells. Cytometry by time-of-flight (CyTOF) analysis provides additional evidence of the role of UPP1 in shaping the immunosuppressive nature of the TME. Using patient-derived organoids (PDOs), we discover that UPP1high tumors exhibit relatively increased sensitivity to Bosutinib and Dasatinib. Collectively, our study highlights the immunosuppressive role of UPP1 in LUAD, and these findings may provide insights into the molecular features of LUAD and facilitate the development of personalized treatment strategies.

Integrative Analyses of Pyrimidine Salvage Pathway-Related Genes Revealing the Associations Between UPP1 and Tumor Microenvironment

Background

The pyrimidine salvage pathway plays a critical role in tumor progression and patient outcomes. The roles of pyrimidine salvage pathway-related genes (PSPGs) in cancer, however, are not fully understood. This study aims to depict the characteristics of PSPGs across various cancers.

Methods

An integrative pan-cancer analysis of six PSPGs (CDA, UCK1, UCK2, UCKL1, UPP1, and UPP2) was conducted using TCGA data, single-cell RNA sequencing datasets, and patient samples. Single-cell transcriptome analysis and RT-qPCR were used to validate the relation between UPP1 and cytokines. Flow cytometry was performed to validate the role of UPP1 in immune checkpoint regulation. The correlation between UPP1 and tumor associated neutrophils (TAN) were investigated and validated by single-cell transcriptome analysis and tissue microarrays (TMAs).

Results

PSPGs showed low mutation rates but significant copy number variations, particularly amplifications in UCKL1, UPP1, and UCK2 across various cancers. DNA methylation patterns varied, with notable negative correlations between methylation and gene expression in UPP1. PSPGs were broadly up-regulated in multiple cancers, with correlations to clinical staging and prognosis. Proteomic data further confirmed these findings. Functional analysis revealed PSPGs' associations with tumor proliferation, metastasis, and various signaling pathways. UPP1 showed strong correlations with the tumor microenvironment (TME), particularly with cytokines, immune checkpoints, and various immune cells. Single-cell transcriptome analysis confirmed these associations, highlighting UPP1's influence on cytokine expression and immune checkpoint regulation. In esophageal squamous cell carcinoma (ESCC), UPP1-high tumor cells were significantly associated with immunosuppressive cells in the TME. Spatial analysis using TMAs revealed that UPP1+ tumor cells were predominantly located at the invasive margin and closely associated with neutrophils, correlating with poorer patient prognosis.

Conclusion

Our study depicted the multi-dimensional view of PSPGs in cancer, with a particular focus on UPP1's role in the TME. Targeting UPP1 holds promise as a potential strategy for cancer therapy.

IGH repertoire analysis at scale: deciphering the complexity of B cell infiltration and migration in esophageal squamous cell carcinoma

Tumor-infiltrating B-lineage cells have become predictors of prognosis and immunotherapy responses in various cancers. However, limited knowledge about their infiltration and migration patterns has hindered the understanding of their anti-tumor functions. Here, we examined the immunoglobulin heavy chain (IGH) repertoires in 496 multi-regional tumor, 107 normal tissue, and 48 metastatic lymph node samples obtained from 107 patients with esophageal squamous cell carcinoma (ESCC). Our study revealed higher IgG-type B-lineage cells infiltration in tumors than in healthy tissue, which was associated with improved patient outcomes. Genes such as ACTN1, COL6A5, and pathways like focal adhesion, which shapes the physical structure of tumors, could affect B-lineage cell infiltration. Notably, the IGH sequence was used as an identity-tag to monitor B cell migration, and their infiltration schema within the tumor were depicted based on our multi-regional tumor specimens. This analysis revealed an escalation in B cell clones overlapped between metastatic lymph nodes and tumors. Therefore, the Lymph Node Activation Index was defined, which could predict the outcomes of patients with lymph node metastasis. This research introduces a novel framework for probing B cell infiltration and migration within the tumor microenvironment using large-scale transcriptome data, while simultaneously providing fresh perspectives on B cell immunology within ESCC.

HDS screening with patient-derived primary cells guided individualized therapy for esophageal squamous cell carcinoma-in vivo and vitro

Objective

To analyze and evaluate the role of the High-throughput Drug Sensitivity (HDS) screening strategy in identifying highly sensitive drugs against esophageal squamous cell carcinoma (ESCC).

Methods

A total of 80 patients with progressive ESCC were randomly divided into the observation (40 cases) and the control groups (40 cases). In the observation group, primary ESCC cells were isolated from the tumor tissues with a gastroscope, and drug sensitivity screening was performed on cells derived from the 40 ESCC cases using the HDS method, followed by verification in a patient-derived tumor xenograft (PDX) mouse model. Finally, the differences in the therapeutic efficacy (levels of CEA, CYFRA21-1, SCCA after chemotherapy and the rates of overall survival, local progression, and distant metastasis at 12 months and 18 months time points after chemotherapy) were compared between the observation group (Screened drug-treated) and the control group (Paclitaxel combined with cisplatin regimen-treated).

Results

Forty ESCC patients were screened for nine different high-sensitive chemotherapeutics, with the majority showing sensitivity to Bortezomib. Experiments on animal models revealed that the tumor tissue mass of PDX mice treated with the HDS-screened drug was significantly lower than that of the Paclitaxel-treated mice (p < 0.05), and the therapeutic efficacy of the observation group was better than the control group (p < 0.05).

Conclusion

HDS screening technology can be beneficial in screening high-efficacy anticancer drugs for advanced-stage ESCC patients, thereby minimizing adverse drug toxicity in critically ill patients. Moreover, this study provides a new avenue for treating advanced ESCC patients with improved outcomes.

Roles of circRNA dysregulation in esophageal squamous cell carcinoma tumor microenvironment

Esophageal squamous cell carcinoma (ESCC) is the most prevalent histological esophageal cancer characterized by advanced diagnosis, metastasis, resistance to treatment, and frequent recurrence. In recent years, numerous human disorders such as ESCC, have been linked to abnormal expression of circular RNAs (circRNAs), suggesting that they are fundamental to the intricate system of gene regulation that governs ESCC formation. The tumor microenvironment (TME), referring to the area surrounding the tumor cells, is composed of multiple components, including stromal cells, immune cells, the vascular system, extracellular matrix (ECM), and numerous signaling molecules. In this review, we briefly described the biological purposes and mechanisms of aberrant circRNA expression in the TME of ESCC, including the immune microenvironment, angiogenesis, epithelial-to-mesenchymal transition, hypoxia, metabolism, and radiotherapy resistance. As in-depth research into the processes of circRNAs in the TME of ESCC continues, circRNAs are promising therapeutic targets or delivery systems for cancer therapy and diagnostic and prognostic indicators for ESCC.

Transcriptomic analysis of esophageal cancer reveals hub genes and networks involved in cancer progression

Esophageal carcinoma (ESCA) has a 5-year survival rate of fewer than 20%. The study aimed to identify new predictive biomarkers for ESCA through transcriptomics meta-analysis to address the problems of ineffective cancer therapy, lack of efficient diagnostic tools, and costly screening and contribute to developing more efficient cancer screening and treatments by identifying new marker genes. Nine GEO datasets of three kinds of esophageal carcinoma were analyzed, and 20 differentially expressed genes were detected in carcinogenic pathways. Network analysis revealed four hub genes, namely RAR Related Orphan Receptor A (RORA), lysine acetyltransferase 2B (KAT2B), Cell Division Cycle 25B (CDC25B), and Epithelial Cell Transforming 2 (ECT2). Overexpression of RORA, KAT2B, and ECT2 was identified with a bad prognosis. These hub genes modulate immune cell infiltration. These hub genes modulate immune cell infiltration. Although this research needs lab confirmation, we found interesting biomarkers in ESCA that may aid in diagnosis and treatment.

A novel basement membrane-related gene signature for prognosis of lung adenocarcinomas

BACKGROUND

Lung adenocarcinoma (LUAD) remains a global health concern with its poor prognosis and high mortality. Whether tumor cells invade through the basement membrane (BM) is the key factor to determine the prognosis of LUAD. This study aimed to identify the BM-related gene signatures to improve the overall prognosis of LUAD.

MATERIALS & METHODS

A series of bioinformatics analyses were conducted based on TCGA and GEO datasets. Unsupervised consistent cluster analysis was performed, and 500 LUAD patients were assigned to two different groups according to expressions of 222 BM-related genes. The differentially expressed genes (DEGs) between the two clusters were identified, and Lasso regression, ROC curve, univariate and multivariate Cox regression analyses and enrichment analysis were conducted. Besides, ssGSEA, CIBERSORT and ESTIMATE algorithmwere were employed to understand the relationship between the tumor microenvironment (TME) and risk scores. Moreover, single cell clustering and trajectory analyses were performed to further understand the significance of BM-related genes. Finally, qRT-PCR was used to verify the prognosis model.

RESULTS

A total of 31 prognostic BM-related genes were determined for LUAD, and a novel 17-mRNA prognostic model named BMsocre was successfully established to predict the overall survival of LUAD patients. The high BMscore group indicated worse prognosis. Seventeen DEGs were enriched mainly in metabolism, ECM-receptor interaction and immune response. In addition, the high-risk group showed higher TMB and lower immune score. The low-risk group had a better immunotherapeutic response where immune escape was less likely. The BMscore model was verified in our patient cohort. Furthermore, NELL2 was mainly expressed in clusters of T cells, and was identified to play a critical role in T-cell differentiation.

CONCLUSIONS

A novel BMscore model was successfully established and might be effective for providing guidance to LUAD therapy.

Artificial Intelligence-Assisted Transcriptomic Analysis to Advance Cancer Immunotherapy

The emergence of immunotherapy has dramatically changed the cancer treatment paradigm and generated tremendous promise in precision medicine. However, cancer immunotherapy is greatly limited by its low response rates and immune-related adverse events. Transcriptomics technology is a promising tool for deciphering the molecular underpinnings of immunotherapy response and therapeutic toxicity. In particular, applying single-cell RNA-seq (scRNA-seq) has deepened our understanding of tumor heterogeneity and the microenvironment, providing powerful help for developing new immunotherapy strategies. Artificial intelligence (AI) technology in transcriptome analysis meets the need for efficient handling and robust results. Specifically, it further extends the application scope of transcriptomic technologies in cancer research. AI-assisted transcriptomic analysis has performed well in exploring the underlying mechanisms of drug resistance and immunotherapy toxicity and predicting therapeutic response, with profound significance in cancer treatment. In this review, we summarized emerging AI-assisted transcriptomic technologies. We then highlighted new insights into cancer immunotherapy based on AI-assisted transcriptomic analysis, focusing on tumor heterogeneity, the tumor microenvironment, immune-related adverse event pathogenesis, drug resistance, and new target discovery. This review summarizes solid evidence for immunotherapy research, which might help the cancer research community overcome the challenges faced by 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.

A more novel and robust gene signature predicts outcome in patients with esophageal squamous cell carcinoma

BACKGROUND

Esophageal squamous cell carcinoma (ESCC) is a life-threatening thoracic tumor with a poor prognosis. The tumor microenvironment (TME) mainly comprises tumor cells and tumor-infiltrating immune cells mixed with stromal components. The latest research has displayed that tumor immune cell infiltration (ICI) is closely connected with the ESCC patients' clinical prognosis. This study was designed to construct a gene signature based on the ICI of ESCC to predict prognosis.

METHODS

Based on the selection criteria we set, the eligible ESCC cases from the GSE53625 and TCGA-ESCA datasets were chosen for the training cohort and the validation cohort, respectively. Unsupervised clustering detailed grouped ESCC cases of the training cohort based on the ICI profile. We determined the differential expression genes (DEGs) between the ICI clusters, and, subsequently, we adopted the univariate Cox analysis to recognize DEGs with prognostic potential. These screened DEGs underwent a Lasso regression, which then generated a gene signature. The harvested signature's predictive ability was further examined by the Kaplan-Meier analysis, Cox analysis, ROC, IAUC, and IBS. More importantly, we listed similar studies in the most recent year and compared theirs with ours. We performed the functional annotation, immune relevant signature correlation analysis, and immune infiltrating analysis to thoroughly understand the functional mechanism of the signature and the immune cells' roles in the gene signature's predicting capacity.

RESULTS

A sixteen-gene signature (ARSD, BCAT1, BIK, CLDN11, DLEU7-AS1, GGH, IGFBP2, LINC01037, LINC01446, LINC01497, M1AP, PCSK2, PCSK5, PPP2R2A, TIGD7, and TMSB4X) was generated from the Lasso model. We then confirmed the signature as having solid and stable prognostic capacity by several statistical methods. We revealed the superiority of our signature after comparing it to our predecessors, and the GSEA uncovered the specifically mechanism of action related to the gene signature. Two immune relevant signatures, including GZMA and LAG3 were identified associating with our signature. The immune-infiltrating analysis identified crucial roles of resting mast cells, which potentially support the sixteen-gene signature's prognosis ability.

CONCLUSIONS

We discovered a robust sixteen-gene signature that can accurately predict ESCC prognosis. The immune relevant signatures, GZMA and LAG3, and resting mast cells infiltrating were closely linked to the sixteen-gene signature's ability.

Drug repositioning for esophageal squamous cell carcinoma

Esophageal cancer (EC) remains a significant challenge globally, having the 8th highest incidence and 6th highest mortality worldwide. Esophageal squamous cell carcinoma (ESCC) is the most common form of EC in Asia. Crucially, more than 90% of EC cases in China are ESCC. The high mortality rate of EC is likely due to the limited number of effective therapeutic options. To increase patient survival, novel therapeutic strategies for EC patients must be devised. Unfortunately, the development of novel drugs also presents its own significant challenges as most novel drugs do not make it to market due to lack of efficacy or safety concerns. A more time and cost-effective strategy is to identify existing drugs, that have already been approved for treatment of other diseases, which can be repurposed to treat EC patients, with drug repositioning. This can be achieved by comparing the gene expression profiles of disease-states with the effect on gene-expression by a given drug. In our analysis, we used previously published microarray data and identified 167 differentially expressed genes (DEGs). Using weighted key driver analysis, 39 key driver genes were then identified. These driver genes were then used in Overlap Analysis and Network Analysis in Pharmomics. By extracting drugs common to both analyses, 24 drugs are predicted to demonstrate therapeutic effect in EC patients. Several of which have already been shown to demonstrate a therapeutic effect in EC, most notably Doxorubicin, which is commonly used to treat EC patients, and Ixazomib, which was recently shown to induce apoptosis and supress growth of EC cell lines. Additionally, our analysis predicts multiple psychiatric drugs, including Venlafaxine, as repositioned drugs. This is in line with recent research which suggests that psychiatric drugs should be investigated for use in gastrointestinal cancers such as EC. Our study shows that a drug repositioning approach is a feasible strategy for identifying novel ESCC therapies and can also improve the understanding of the mechanisms underlying the drug targets.

Computational exploration of cellular communication in skin from emerging single-cell and spatial transcriptomic data

Tissue development and homeostasis require coordinated cell–cell communication. Recent advances in single-cell sequencing technologies have emerged as a revolutionary method to reveal cellular heterogeneity with unprecedented resolution. This offers a great opportunity to explore cell–cell communication in tissues systematically and comprehensively, and to further identify signaling mechanisms driving cell fate decisions and shaping tissue phenotypes. Using gene expression information from single-cell transcriptomics, several computational tools have been developed for inferring cell–cell communication, greatly facilitating analysis and interpretation. However, in single-cell transcriptomics, spatial information of cells is inherently lost. Given that most cell signaling events occur within a limited distance in tissues, incorporating spatial information into cell–cell communication analysis is critical for understanding tissue organization and function. Spatial transcriptomics provides spatial location of cell subsets along with their gene expression, leading to new directions for leveraging spatial information to develop computational approaches for cell–cell communication inference and analysis. These computational approaches have been successfully applied to uncover previously unrecognized mechanisms of intercellular communication within various contexts and across organ systems, including the skin, a formidable model to study mechanisms of cell–cell communication due to the complex interactions between the different cell populations that comprise it. Here, we review emergent cell–cell communication inference tools using single-cell transcriptomics and spatial transcriptomics, and highlight the biological insights gained by applying these computational tools to exploring cellular communication in skin development, homeostasis, disease and aging, as well as discuss future potential research avenues.