This study aimed to identified the key genes and sequencing metrics for predicting prognosis and efficacy of neoadjuvant chemotherapy (nCT) in rectal cancer (RC) based on genomic DNA sequencing in samples with different origin and multi-omics association database.
We collected 16 RC patients and obtained DNA sequencing data from cancer tissues and plasma cell-free DNA before and after nCT. Various gene variations were analyzed, including single nucleotide variants (SNV), copy number variation (CNV), tumor mutation burden (TMB), copy number instability (CNI) and mutant-allele tumor heterogeneity (MATH). We also identified genes by which CNV level can differentiate the response to nCT. The Cancer Genome Atlas database and the Clinical Proteomic Tumor Analysis Consortium database were used to further evaluate the specific role of therapeutic relevant genes and screen out the key genes in multi-omics levels. After the intersection of the screened genes from differential expression analysis, survival analysis and principal components analysis dimensionality reduction cluster analysis, the key genes were finally identified.
The genes CNV level of principal component genes in baseline blood and cancer tissues could significantly distinguish the two groups of patients. The CNV of HSP90AA1, EGFR, SRC, MTOR, etc. were relatively gained in the better group compared with the poor group in baseline blood. The CNI and TMB was significantly different between the two groups. The increased expression of HSP90AA1, EGFR, and SRC was associated with increased sensitivity to multiple chemotherapeutic drugs. The nCT predictive score obtained by therapeutic relevant genes could be a potential prognostic indicator, and the combination with TMB could further refine prognostic prediction for patients. After a series of analysis in multi-omics association database, EGFR and HSP90AA1 with significant differences in multiple aspects were identified as the key predictive genes related to prognosis and the sensitivity of nCT.
This work revealed that effective combined application and analysis in multi-omics data are critical to search for predictive biomarkers. The key genes EGFR and HSP90AA1 could serve as an effective biomarker to predict prognose and neoadjuvant chemosensitivity.