Three-dimensional mRNA measurements reveal minimal regional heterogeneity in esophageal squamous cell carcinoma

Global comparative transcriptomes uncover novel and population-specific gene expression in esophageal squamous cell carcinoma

BACKGROUND

Esophageal squamous cell carcinoma (ESCC) has a poor prognosis and is one of the deadliest gastrointestinal malignancies. Despite numerous transcriptomics studies to understand its molecular basis, the impact of population-specific differences on this disease remains unexplored.

AIMS

This study aimed to investigate the population-specific differences in gene expression patterns among ESCC samples obtained from six distinct global populations, identify differentially expressed genes (DEGs) and their associated pathways, and identify potential biomarkers for ESCC diagnosis and prognosis. In addition, this study deciphers population specific microbial and chemical risk factors in ESCC.

METHODS

We compared the gene expression patterns of ESCC samples from six different global populations by analyzing microarray datasets. To identify DEGs, we conducted stringent quality control and employed linear modeling. We cross-compared the resulting DEG lists of each populations along with ESCC ATLAS to identify known and novel DEGs. We performed a survival analysis using The Cancer Genome Atlas Program (TCGA) data to identify potential biomarkers for ESCC diagnosis and prognosis among the novel DEGs. Finally, we performed comparative functional enrichment and toxicogenomic analysis.

RESULTS

Here we report 19 genes with distinct expression patterns among populations, indicating population-specific variations in ESCC. Additionally, we discovered 166 novel DEGs, such as ENDOU, SLCO1B3, KCNS3, IFI35, among others. The survival analysis identified three novel genes (CHRM3, CREG2, H2AC6) critical for ESCC survival. Notably, our findings showed that ECM-related gene ontology terms and pathways were significantly enriched among the DEGs in ESCC. We also found population-specific variations in immune response and microbial infection-related pathways which included genes enriched for HPV, Ameobiosis, Leishmaniosis, and Human Cytomegaloviruses. Our toxicogenomic analysis identified tobacco smoking as the primary risk factor and cisplatin as the main drug chemical interacting with the maximum number of DEGs across populations.

CONCLUSION

This study provides new insights into population-specific differences in gene expression patterns and their associated pathways in ESCC. Our findings suggest that changes in extracellular matrix (ECM) organization may be crucial to the development and progression of this cancer, and that environmental and genetic factors play important roles in the disease. The novel DEGs identified may serve as potential biomarkers for diagnosis, prognosis and treatment.

CD39-Expressing CD8+ T Cells as a New Molecular Marker for Diagnosis and Prognosis of Esophageal Squamous Cell Carcinoma

We aimed to explore the effect of CD39 expression on CD8⁺ T cells and on the diagnosis and prognosis of esophageal squamous cell carcinoma (ESCC). The independent prognostic factors for the surgical specimens of the 95 ESCC patients were screened by multivariate Cox regression analysis. Differential gene expression analysis was performed by the NetworkAnalyst platform based on data from the Gene Expression Omnibus (GEO). The expression of CD39 on CD8⁺ T cells in the CK⁺ region was higher in cancer tissue than in paracancerous tissue (p = 0.011), and high CD39-expressing CD8⁺ T cells in the CK⁺ region (HR, 2.587; p = 0.033) and high CD39-expressing CD8⁺ T cells in the CK^- region (HR, 3.090; p = 0.008) were independent risk factors for prognosis in ESCC patients; the expression of ENTPD1 was upregulated in ESCC tissues compared to normal tissues (adjusted p < 0.001; log₂ fold change = 1.99), and its expression was significantly positively correlated with the expression of PDCD1, CTLA4, and HAVCR2. High CD39-expressing CD8⁺ T cells can be used as a new molecular marker for the diagnosis and prognosis of ESCC, and the restoration of partially exhausted CD8⁺ T cells by inhibiting CD39 may be a new strategy for treating ESCC.

A comprehensive analysis of mRNA expression profiles of Esophageal Squamous Cell Carcinoma reveals downregulation of Desmoglein 1 and crucial genomic targets

AIM

Esophageal Squamous Cell Carcinoma (ESCC) is a histological subtype of esophageal cancer that begins in the squamous cells in the esophagus. In only 19% of the ESCC-diagnosed patients, a five-year survival rate has been seen. This necessitates the identification of high-confidence biomarkers for early diagnosis, prognosis, and potential therapeutic targets for the mitigation of ESCC.

METHOD

We performed a meta-analysis of 10 mRNA datasets and identified consistently perturbed genes across the studies. Then, integrated with ESCC ATLAS to segregate 'core' genes to identify consequences of primary gene perturbation events leading to gene-gene interactions and dysregulated molecular signaling pathways. Further, by integrating with toxicogenomics data, inferences were drawn for gene interaction with environmental exposures, trace elements, chemical carcinogens, and drug chemicals. We also deduce the clinical outcomes of candidate genes based on survival analysis using the ESCC related dataset in The Cancer Genome Atlas.

RESULT

We identified 237 known and 18 novel perturbed candidate genes. Desmoglein 1 (DSG1) is one such gene that we found significantly downregulated (Fold Change =-1.89, p-value = 8.2e-06) in ESCC across six different datasets. Further, we identified 31 'core' genes (that either harbor genetic variants or are regulated by epigenetic modifications) and found regulating key biological pathways via adjoining genes in gene-gene interaction networks. Functional enrichment analysis showed dysregulated biological processes and pathways including "Extracellular matrix", "Collagen trimmer" and "HPV infection" are significantly overrepresented in our candidate genes. Based on the toxicogenomic inferences from Comparative Toxicogenomics Database we report the key genes that interacted with risk factors such as tobacco smoking, zinc, nitroso benzylmethylamine, and drug chemicals such as cisplatin, Fluorouracil, and Mitomycin in relation to ESCC. We also point to the STC2 gene that shows a high risk for mortality in ESCC patients.

CONCLUSION

We identified novel perturbed genes in relation to ESCC and explored their interaction network. DSG1 is one such gene, its association with microbiota and a clinical presentation seen commonly with ESCC hints that it is a good candidate for early diagnostic marker. Besides, in this study we highlight candidate genes and their molecular connections to risk factors, biological pathways, drug chemicals, and the survival probability of ESCC patients.

Multi-region sequencing with spatial information enables accurate heterogeneity estimation and risk stratification in liver cancer

BACKGROUND

Numerous studies have used multi-region sampling approaches to characterize intra-tumor heterogeneity (ITH) in hepatocellular carcinoma (HCC). However, conventional multi-region sampling strategies do not preserve the spatial details of samples, and thus, the potential influences of spatial distribution on patient-wise ITH (represents the overall heterogeneity level of the tumor in a given patient) have long been overlooked. Furthermore, gene-wise transcriptional ITH (represents the expression pattern of genes across different intra-tumor regions) in HCC is also under-explored, highlighting the need for a comprehensive investigation.

METHODS

To address the problem of spatial information loss, we propose a simple and easy-to-implement strategy called spatial localization sampling (SLS). We performed multi-region sampling and sequencing on 14 patients with HCC, collecting a total of 75 tumor samples with spatial information and molecular data. Normalized diversity score and integrated heterogeneity score (IHS) were then developed to measure patient-wise and gene-wise ITH, respectively.

RESULTS

A significant correlation between spatial and molecular heterogeneity was uncovered, implying that spatial distribution of sampling sites did influence ITH estimation in HCC. We demonstrated that the normalized diversity score had the ability to overcome sampling location bias and provide a more accurate estimation of patient-wise ITH. According to this metric, HCC tumors could be divided into two classes (low-ITH and high-ITH tumors) with significant differences in multiple biological properties. Through IHS analysis, we revealed a highly heterogenous immune microenvironment in HCC and identified some low-ITH checkpoint genes with immunotherapeutic potential. We also constructed a low-heterogeneity risk stratification (LHRS) signature based on the IHS results which could accurately predict the survival outcome of patients with HCC on a single tumor biopsy sample.

CONCLUSIONS

This study provides new insights into the complex phenotypes of HCC and may serve as a guide for future studies in this field.

LAMP2 as a Biomarker Related to Prognosis and Immune Infiltration in Esophageal Cancer and Other Cancers: A Comprehensive Pan-Cancer Analysis

Esophageal cancer (ESCA) is a common malignant tumor with poor prognosis. Accumulating evidence indicates an important role of lysosomal-associated membrane protein 2 (LAMP2) in the progression and development of various cancers. In this study, we obtained RNA-sequencing raw count data and the corresponding clinical information for ESCA samples from The Cancer Genome Atlas and Gene Expression Omnibus databases. We comprehensively investigated the expression and prognostic significance of LAMP2 and relationships between LAMP2 expression and prognosis, different clinicopathological parameters, and immune cell infiltration in ESCA. We also obtained the differentially expressed genes between the high LAMP2 expression and low LAMP2 expression groups in ESCA and performed a functional enrichment analysis of the 250 linked genes most positively related to LAMP2 expression. Moreover, we performed the pan-cancer analysis of LAMP2 to further analyze the role of LAMP2 in 25 commonly occurring types of human cancer. We also verified and compared the expression of LAMP2 in 40 samples of human ESCA tissue and adjacent tissues. The results indicated that LAMP2 expression was significantly upregulated in ESCA and various human cancers. In addition, LAMP2 expression was associated with certain clinicopathological parameters, prognosis, and immune infiltration in ESCA and the other types of cancer. Our study represents a comprehensive pan-cancer analysis of LAMP2 and supports the potential use of the modulation of LAMP2 in the management of ESCA and various cancers.

Genomic Characterization and Comparison of Multi-Regional and Pooled Tumor Biopsy Specimens

A single tumor biopsy specimen is typically used in cancer genome studies. However, it may represent incompletely the underlying mutational and transcriptional profiles of tumor biology. Multi-regional biopsies have the advantage of increased sensitivity for genomic profiling, but they are not cost-effective. The concept of an alternative method such as the pooling of multiple biopsies is a challenge. In order to determine if the pooling of distinct regions is representative at the genomic and transcriptome level, we performed sequencing of four regional samples and pooled samples for four cancer types including colon, stomach, kidney and liver cancer. Subsequently, a comparative analysis was conducted to explore differences in mutations and gene expression profiles between multiple regional biopsies and pooled biopsy for each tumor. Our analysis revealed a marginal level of regional difference in detected variants, but in those with low allele frequency, considerable discrepancies were observed. In conclusion, sequencing pooled samples has the benefit of detecting many variants with moderate allele frequency that occur in partial regions, but it is not applicable for detecting low-frequency mutations that require deep sequencing.