A pan-cancer analysis of collagen VI family on prognosis, tumor microenvironment, and its potential therapeutic effect

Label-based comparative proteomics of oral mucosal tissue to understand progression of precancerous lesions to oral squamous cell carcinoma

Introduction

Oral squamous cell carcinomas typically arise from precancerous lesions such as leukoplakia and erythroplakia. These lesions exhibit a range of histological changes from hyperplasia to dysplasia and carcinoma in situ, during their transformation to malignancy. The molecular mechanisms driving this multistage transition remain incompletely understood. To bridge this knowledge gap, our current study utilizes label based comparative proteomics to compare protein expression profiles across different histopathological grades of leukoplakia, erythroplakia, and oral squamous cell carcinoma samples, aiming to elucidate the molecular changes underlying lesion evolution.

Methodology

An 8-plex iTRAQ proteomics of 4 biological replicates from 8 clinical phenotypes of leukoplakia and erythroplakia, with hyperplasia, mild dysplasia, moderate dysplasia; along with phenotypes of well differentiated squamous cell carcinoma and moderately differentiated squamous cell carcinoma was carried out using the Orbitrap Fusion Lumos mass spectrometer. Raw files were processed with Maxquant, and statistical analysis across groups was conducted using MetaboAnalyst. Statistical tools such as ANOVA, PLS-DA VIP scoring, and correlation analysis were employed to identify differentially expressed proteins that had a linear expression variation across phenotypes of hyperplasia to cancer. Validation was done using Bioinformatic tools such as ClueGO + Cluepedia plugin in Cytoscape to extract functional annotations from gene ontology and pathway databases.

Results and discussion

A total of 2685 protein groups and 12,397 unique peptides were identified, and 61 proteins consistently exhibited valid reporter ion corrected intensities across all samples. Of these, 6 proteins showed linear varying expression across the analysed sample phenotypes. Collagen type VI alpha 2 chain (COL6A2), Fibrinogen β chain (FGB), and Vimentin (VIM) were found to have increased linear expression across pre-cancer phenotypes of leukoplakia to cancer, while Annexin A7 (ANXA7) was seen to be having a linear decreasing expression. Collagen type VI alpha 2 chain (COL6A2) and Annexin A2 (ANXA2) had increased linear expression across precancer phenotypes of erythroplakia to cancer. The mass spectrometry proteomics data have been deposited to the ProteomeXchanger Consortium via the PRIDE partner repository with the data set identifier PXD054190. These differentially expressed proteins mediate cancer progression mainly through extracellular exosome; collagen-containing extracellular matrix, hemostasis, platelet aggregation, and cell adhesion molecule binding.

Conclusion

Label-based proteomics is an ideal platform to study oral cancer progression. The differentially expressed proteins provide insights into the molecular mechanisms underlying the progression of oral premalignant lesions to malignant phenotypes. The study has translational value for early detection, risk stratification, and potential therapeutic targeting of oral premalignant lesions and in its prevention to malignant forms.

Revealing the structural microenvironment of high metastatic risk uveal melanomas following decellularisation

Uveal melanoma (UM) is a rare aggressive intraocular tumour that spreads most commonly to the liver in tumours with loss of one copy of chromosome 3 (HR-M3); current treatments for metastatic disease remain largely ineffective. Pre-clinical research is increasingly using three-dimensional models that better recapitulate the tumour microenvironment (TME). One aspect of the TME is the acellular extracellular matrix (ECM) that influences cell proliferation, migration and response to therapy. Although commercial matrices are used in culture, the composition and biochemical properties may not be representative of the tumour ECM in vivo. This study identifies UM metastatic risk specific ECM proteins by developing methodology for decellularisation of low- and high- metastatic risk tissue samples (LR-D3 vs. HR-M3). Proteomic analysis revealed a matrisome signature of 34 core ECM and ECM-associated proteins upregulated in HR-M3 UM. Combining additional UM secretome and whole cell iTRAQ proteomic datasets revealed enriched GO and KEGG pathways including 'regulating ECM binding' and 'PI3K/Akt signalling'. Structural analyses of decellularised matrices revealed microarchitecture of differing fibre density and expression differences in collagen 4, collagen 6A1 and nidogen 1, between metastatic risk groups. This approach is a powerful tool for the generation of ECM matrices relevant to high metastatic risk UM.

Matrisomics: Beyond the extracellular matrix for unveiling tumor microenvironment

The matrisome, a group of proteins constituting or interacting with the extracellular matrix (ECM), has garnered attention as a potent regulator of cancer progression. An increasing number of studies have focused on cancer matrisome utilizing diverse -omics approaches. Here, we present diverse patterns of matrisomal populations within cancer tissues, exploring recent -omics studies spanning different '-omics' levels (epigenomics, genomics, transcriptomics, and proteomics), as well as newly developed sequencing techniques such as single-cell RNA sequencing and spatial transcriptomics. Some matrisome genes showed uniform patterns of upregulated or downregulated expression across various cancers, while others displayed different expression patterns according to the cancer types. This matrisomal dysregulation in cancer was further examined according to their originating cell type and spatial location in the tumor tissue. Experimental studies were also collected to demonstrate the identified roles of matrisome genes during cancer progression. Interestingly, many studies on cancer matrisome have suggested matrisome genes as effective biomarkers in cancer research. Although the specific mechanisms and clinical applications of cancer matrisome have not yet been fully elucidated, recent techniques and analyses on cancer matrisomics have emphasized their biological importance in cancer progression and their clinical implications in deciding the efficacy of cancer treatment.

Establishment of homotrimer collagen type I signature and its association with clinical manifestation and tertiary lymphoid structures formation in liver cancer

Background

collagen type I is a fundamental composition of extracellular matrix. Typically it exists in the form of a heterotrimer, consisting of two α1 chains encoded by COL1A1 and one α2 chain encoded by COL1A2. However, in cancer a homotrimeric form of collagen type I comprises three α1 chains encoded by COL1A1 was founded. There is still a lack of transcriptional and histologic methods for detecting homotrimeric collagen type I. Furthermore, a comprehensive analysis of the pan-cancer distribution pattern and clinical relevance of homotrimeric collagen type I is conspicuously absent.

Method

Using transcriptional and immunoflourance method, we established homocol signature, which is able to transcriptionally and histologically detect homotrimeric collagen type I. We investigated the diagnostic and prognostic potential of homocol as a novel cancer biomarker in a pan-cancer cohort. Furthermore, we assessed its association with clinical manifestations in a liver cancer cohort undergoing treatment at our institute.

Result

Homotrimer Collagen Type I is predominantly expressed by cancer cells and is linked to several critical cancer hallmarks, particularly inflammatory response and proliferation. Survival analyses have indicated that a high Homocol expression is correlated with poor outcomes in most types of cancer studied. In terms of cancer detection, Homocol demonstrated strong performance in Receiver Operating Characteristic (ROC) analysis, with an Area Under Curve (AUC) of 0.83 for pan-cancer detection and between 0.72 and 0.99 for individual cancers.In cohorts undergoing PD1 treatment, we noted a higher presence of Homocol in the response group. In a Hepatocellular Carcinoma (HCC) clinical set, high Homocol expression was associated with an increased formation of intra-tumor tertiary lymphoid structures (TLS), larger tumor sizes, more advanced Barcelona Clinic Liver Cancer (BCLC) stages, higher microvascular invasion (MVI) grades, absence of a capsule, and an enriched para-tumor collagen presence.

Conclusion

our research has led to the development of a novel gene signature that facilitates the detection of Homotrimer Collagen Type I. This may greatly assist efforts in cancer detection, prognosis, treatment response prediction, and further research into Homotrimer Collagen Type I.

Construction and analysis of a reliable five-gene prognostic signature for colon adenocarcinoma associated with the wild-type allelic state of the COL6A6 gene

Background

Tumors emerge by acquiring a number of mutations over time. The first mutation provides a selective growth advantage compared to adjacent epithelial cells, allowing the cell to create a clone that can outgrow the cells that surround it. Subsequent mutations determine the risk of the tumor progressing to metastatic cancer. Some secondary mutations may inhibit the aggressiveness of the tumor while still increasing the survival of the clone. Meaningful mutations in genes may provide a strong molecular foundation for developing novel therapeutic strategies for cancer.

Methods

The somatic mutation and prognosis in colon adenocarcinoma (COAD) were analyzed. The copy number variation (CNV) and differentially expressed genes (DEGs) between the collagen type VI alpha 6 chain (COL6A6) mutation (COL6A6-MUT) and the COL6A6 wild-type (COL6A6-WT) subgroups were evaluated. The independent prognostic signatures based on COL6A6-allelic state were determined to construct a Cox model. The biological characteristics and the immune microenvironment between the two risk groups were compared.

Results

COL6A6 was found to be highly mutated in COAD at a frequency of 9%. Patients with COL6A6-MUT had a good overall survival (OS) compared to those with COL6A6-WT, who had a different CNV pattern. Significant differences in gene expression were established for 593 genes between the COL6A6-MUT and COL6A6-WT samples. Among them, MUC16, ASNSP1, PRR18, PEG10, and RPL26P8 were determined to be independent prognostic factors. The internally validated prognostic risk model, constructed using these five genes, demonstrated its value by revealing a significant difference in patient prognosis between the high-risk and low-risk groups. Specifically, patients in the high-risk group exhibited a considerably worse prognosis than did those in the low-risk group. The high-risk group had a significantly higher proportion of patients over 60 years of age and patients in stage III. Moreover, the tumor immune dysfunction and exclusion (TIDE) score and the expression of human leukocyte antigen (HLA) family genes were all higher in the high-risk group than that in the low-risk group.

Conclusions

The allelic state of COL6A6 and the five associated DEGs were identified as novel biomarkers for the diagnosis and prognosis of COAD and may be therapeutic targets in COAD.

Transcript and protein signatures derived from shared molecular interactions across cancers are associated with mortality

BACKGROUND

Characterization of shared cancer mechanisms have been proposed to improve therapy strategies and prognosis. Here, we aimed to identify shared cell-cell interactions (CCIs) within the tumor microenvironment across multiple solid cancers and assess their association with cancer mortality.

METHODS

CCIs of each cancer were identified by NicheNet analysis of single-cell RNA sequencing data from breast, colon, liver, lung, and ovarian cancers. These CCIs were used to construct a shared multi-cellular tumor model (shared-MCTM) representing common CCIs across cancers. A gene signature was identified from the shared-MCTM and tested on the mRNA and protein level in two large independent cohorts: The Cancer Genome Atlas (TCGA, 9185 tumor samples and 727 controls across 22 cancers) and UK biobank (UKBB, 10,384 cancer patients and 5063 controls with proteomics data across 17 cancers). Cox proportional hazards models were used to evaluate the association of the signature with 10-year all-cause mortality, including sex-specific analysis.

RESULTS

A shared-MCTM was derived from five individual cancers. A shared gene signature was extracted from this shared-MCTM and the most prominent regulatory cell type, matrix cancer-associated fibroblast (mCAF). The signature exhibited significant expression changes in multiple cancers compared to controls at both mRNA and protein levels in two independent cohorts. Importantly, it was significantly associated with mortality in cancer patients in both cohorts. The highest hazard ratios were observed for brain cancer in TCGA (HR [95%CI] = 6.90[4.64-10.25]) and ovarian cancer in UKBB (5.53[2.08-8.80]). Sex-specific analysis revealed distinct risks, with a higher mortality risk associated with the protein signature score in males (2.41[1.97-2.96]) compared to females (1.84[1.44-2.37]).

CONCLUSION

We identified a gene signature from a comprehensive shared-MCTM representing common CCIs across different cancers and revealed the regulatory role of mCAF in the tumor microenvironment. The pathogenic relevance of the gene signature was supported by differential expression and association with mortality on both mRNA and protein levels in two independent cohorts.

A pathway-based computational framework for identification of a new modal of multi-omics biomarkers and its application in esophageal cancer

BACKGROUND

The pathway-based strategy has been recently proposed for identifying biomarkers with the advantages of higher biological interpretability and cross-data robustness than the conventional gene-based strategy. However, its utility in clinical applications has been limited due to the high computational complexity and ill-defined performance.

OBJECTIVE

The current study presents a machine learning-based computational framework using multi-omics data for identifying a new modal of biomarkers, called pathway-derived core biomarkers, which have the advantages of both gene-based and pathway-based biomarkers.

METHODS

Machine-learning methods and gene-pathway network were integrated to select the pathway-derived core biomarkers. Multiple machine-learning algorithms were used to construct and validate the diagnostic models of the biomarkers based on more than 1400 multi-omics clinical samples of esophageal squamous cell carcinoma (ESCC).

RESULTS

The results showed that the classifier models based on the new modal biomarkers achieved superior performance in the training datasets with an average AUC/accuracy of 0.98/0.95 and 0.89/0.81 for mRNAs and miRNA, respectively, higher than the currently known classifier models based on the conventional gene-based strategy and pathway-based strategy. In the testing cohorts, the AUC/accuracy increased by 6.1 %/7.3 % than the models based on the native gene-based biomarkers. The improved performance was further confirmed in independent validation cohorts. Specifically, the sensitivity/specificity increased by ∼3 % and the variance significantly decreased by ∼69 % compared with that of the native gene-based biomarkers. Importantly, the pathway-derived core biomarkers also recovered 45 % more previously reported biomarkers than the gene-based biomarkers and are more functionally relevant to the ESCC etiology (involved in 14 versus 7 pathways related with ESCC or other cancer), highlighting the cross-data robustness of this new modal of biomarkers via enhanced functional relevance.

CONCLUSIONS

The results demonstrated that the new modal of biomarkers not only have improved predicting performance and robustness, but also exhibit higher functional interpretability thus leading to the potential application in cancer diagnosis.

Collagen code in tumor microenvironment: Functions, molecular mechanisms, and therapeutic implications

The tumor microenvironment (TME) is crucial in cancer progression, and the extracellular matrix (ECM) is an important TME component. Collagen is a major ECM component that contributes to tumor cell infiltration, expansion, and distant metastasis during cancer progression. Recent studies reported that collagen is deposited in the TME to form a collagen wall along which tumor cells can infiltrate and prevent drugs from working on the tumor cells. Collagen-tumor cell interaction is complex and requires the activation of multiple signaling pathways for biochemical and mechanical signaling interventions. In this review, we examine the effect of collagen deposition in the TME on tumor progression and discuss the interaction between collagen and tumor cells. This review aims to illustrate the functions and mechanisms of collagen in tumor progression in the TME and its role in tumor therapy. The findings indicated collagen in the TME appears to be a better target for cancer therapy.

Inferring spatial transcriptomics markers from whole slide images to characterize metastasis-related spatial heterogeneity of colorectal tumors: A pilot study

Over 150 000 Americans are diagnosed with colorectal cancer (CRC) every year, and annually over 50 000 individuals will die from CRC, necessitating improvements in screening, prognostication, disease management, and therapeutic options. Tumor metastasis is the primary factor related to the risk of recurrence and mortality. Yet, screening for nodal and distant metastasis is costly, and invasive and incomplete resection may hamper adequate assessment. Signatures of the tumor-immune microenvironment (TIME) at the primary site can provide valuable insights into the aggressiveness of the tumor and the effectiveness of various treatment options. Spatially resolved transcriptomics technologies offer an unprecedented characterization of TIME through high multiplexing, yet their scope is constrained by cost. Meanwhile, it has long been suspected that histological, cytological, and macroarchitectural tissue characteristics correlate well with molecular information (e.g., gene expression). Thus, a method for predicting transcriptomics data through inference of RNA patterns from whole slide images (WSI) is a key step in studying metastasis at scale. In this work, we collected tissue from 4 stage-III (pT3) matched colorectal cancer patients for spatial transcriptomics profiling. The Visium spatial transcriptomics (ST) assay was used to measure transcript abundance for 17 943 genes at up to 5000 55-micron (i.e., 1-10 cells) spots per patient sampled in a honeycomb pattern, co-registered with hematoxylin and eosin (H&E) stained WSI. The Visium ST assay can measure expression at these spots through tissue permeabilization of mRNAs, which are captured through spatially (i.e., x-y positional coordinates) barcoded, gene specific oligo probes. WSI subimages were extracted around each co-registered Visium spot and were used to predict the expression at these spots using machine learning models. We prototyped and compared several convolutional, transformer, and graph convolutional neural networks to predict spatial RNA patterns at the Visium spots under the hypothesis that the transformer- and graph-based approaches better capture relevant spatial tissue architecture. We further analyzed the model's ability to recapitulate spatial autocorrelation statistics using SPARK and SpatialDE. Overall, the results indicate that the transformer- and graph-based approaches were unable to outperform the convolutional neural network architecture, though they exhibited optimal performance for relevant disease-associated genes. Initial findings suggest that different neural networks that operate on different scales are relevant for capturing distinct disease pathways (e.g., epithelial to mesenchymal transition). We add further evidence that deep learning models can accurately predict gene expression in whole slide images and comment on understudied factors which may increase its external applicability (e.g., tissue context). Our preliminary work will motivate further investigation of inference for molecular patterns from whole slide images as metastasis predictors and in other applications.

Nidogen-2 (NID2) is a Key Factor in Collagen Causing Poor Response to Immunotherapy in Melanoma

Background

The incidence of cutaneous melanoma continues to rise rapidly and has an extremely poor prognosis. Immunotherapy strategies are the most effective approach for patients who have developed metastases, but not all cases have been successful due to the complex and variable mechanisms of melanoma response to immune checkpoint inhibition.

Methods

We synthesized collagen-coding gene expression data (second-generation and single-cell sequencing) from public Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases. Bioinformatics analysis was performed using R software and several database resources such as Metascape database, Gene Set Cancer Analysis (GSCA) database, and Cytoscape software, etc., to investigate the biological mechanisms that may be related with collagens. Immunofluorescence and immunohistochemical staining were used to validate the expression and localization of Nidogen-2 (NID2).

Results

Melanoma patients can be divided into two collagen clusters. Patients with high collagen levels (C1) had a shorter survival than those with low collagen levels (C2) and were less likely to benefit from immunotherapy. We demonstrated that NID2 is a potential key factor in the collagen phenotype, is involved in fibroblast activation in melanoma, and forms a barrier to limit the proximity of CD8+ T cells to tumor cells.

Conclusion

We clarified the adverse effects of collagen on melanoma patients and identified NID2 as a potential therapeutic target.