Comparative sequence analysis of patient-matched primary colorectal cancer, metastatic, and recurrent metastatic tumors after adjuvant FOLFOX chemotherapy

Exploring microRNA patterns as biomarkers of FOLFOX chemotherapy-induced peripheral neuropathy in patients with colorectal cancer

FOLFOX is a combination of chemotherapeutic agents (5-fluorouracil, leucovorin, and oxaliplatin) and is used to treat advanced colorectal cancer (CRC) but induces various side effects. Chemotherapy-induced peripheral neuropathy (CIPN) is one of the most critical side effects that compromise the quality of life of patients with CRC undergoing FOLFOX chemotherapy. This study aimed to evaluate circulating miRNA, cortisol and catecholamine as potential biomarkers that can predict FOLFOX-CIPN symptoms. High-throughput microRNA (miRNA) sequencing was performed on the RNA circulating in the plasma of eight patients with CRC who underwent FOLFOX chemotherapy. miRNA expression profiles were evaluated according to two groups: those who underwent ≤3 cycles and those who underwent ≥6 cycles of FOLFOX chemotherapy. The identified miRNAs were validated in 27 patients with CRC who underwent FOLFOX chemotherapy using quantitative reverse transcription polymerase chain reaction. Target genes were predicted using bioinformatics and functional analyses. Cortisol and catecholamine concentrations in peripheral plasma were measured using an enzyme-linked immunosorbent assay. miR-3184-5p was differentially expressed when miRNA expression was compared between the groups that underwent ≤3 and ≥6 cycles of FOLFOX chemotherapy. Cortisol levels were significantly higher in the group that underwent ≥6 cycles of FOLFOX chemotherapy than in the group that underwent ≤3 cycles. This study suggests that miR-3184-5p may be a potential marker for predicting CIPN.

Identification of Hub Genes Associated with Gastric Cancer via Bioinformatics Analysis and Validation Studies

Introduction

Hub genes related to the development of gastric cancer (GC) were identified based on bioinformatics methods. This study aimed to identify GC hub genes, explore the expression of genes in GC and their correlation with prognosis, so as to provide strategies for GC diagnosis and targeted therapy.

Methods

Two messenger RNA (mRNA) microarray datasets were downloaded from GEO database. These data were combined with TCGA database to obtain common DEGs between GC tissues and normal tissues. GO and KEGG pathway enrichment analysis was performed. Visualized PPI network analysis was performed by Cytoscape to further identify hub genes. GEPIA database was used to evaluate the prognostic value of hub genes. The online software Ualcan was applied to analyze the expression of the prognosis-related genes in cancer tissues and normal tissues from different perspectives of primary GC, TNM stage, nodal metastasis status and tumor grade. Immunohistochemical staining of GC tissues and normal tissues was performed to evaluate the expression of signature genes in GC.

Results

Eighty-four common differentially expressed genes (DEGs) in GC were identified. These genes were closely related to the P13K-Akt signal pathway and other signaling pathways. Ten hub genes were identified. Collagen type I alpha 1 (COL1A1) and collagen type IV alpha 1 (COL4A1) were significantly associated with poor prognosis of GC and were all positively correlated with T stage, distant metastasis, and TNM stage of GC. Immunohistochemistry revealed that the expression of these 2 genes was upregulated in GC tissues. These 2 genes expression was negatively related with 5-year survival rate of GC patients.

Conclusion

Ten highly expressed hub genes in GC tissue were mined by bioinformatics method. COL1A1 and COL4A1 were significantly associated with the prognosis of GC. This study provided a theoretical basis for the pathogenesis, clinical diagnosis and therapeutic targets of GC.

Glutathione peroxidase 2: A key factor in the development of microsatellite instability in colon cancer

BACKGROUND

Much research has focused on detecting microsatellite instability (MSI), which is frequently employed in the diagnosis and treatment of patients with colon cancer. However, the causes and progression of MSI in colon cancer have not yet been thoroughly elucidated. In this study, we screened and validated the genes associated with MSI in colorectal adenocarcinoma (COAD) using bioinformatics analysis.

METHODS

MSI-related genes of COAD were obtained from the Gene Expression Omnibus dataset, Search Tool for the Retrieval of Interaction Gene/Proteins, Gene Set Enrichment Analysis, and Human Protein Atlas. The function, prognostic value, and immune connection of MSI-related genes in COAD were examined using Cytoscape 3.9.1, the Human Gene Database, and the Tumor IMmune Estimation Resource. Key genes were verified using The Cancer Genome Atlas database and immunohistochemistry of clinical tumor samples.

RESULTS

We identified 59 MSI-related genes in patients with colon cancer. The protein interaction network of these genes was developed, and numerous functional modules associated with MSI were discovered. Pathways related to MSI were identified using KEGG enrichment analysis, and these included chemokine signaling, thyroid hormone synthesis, cytokine receptor interaction, estrogen signaling, and Wnt signaling pathways. Further analyses were used to identify the MSI-related gene, glutathione peroxidase 2 (GPX2), which was closely related to the occurrence of COAD and tumor immunity.

CONCLUSIONS

In COAD, GPX2 may be crucial for the establishment of MSI and tumor immunity, and its deficiency may result in MSI and immune cell infiltration in colon cancer.

A novel prognostic gene set for colon adenocarcinoma relative to the tumor microenvironment, chemotherapy, and immune therapy

Background: Colon adenocarcinoma (COAD) is a common aggressive malignant tumor. Heterogeneity in tumorigenesis and therapy response leads to an unsatisfactory overall survival of colon adenocarcinoma patients. Our study aimed to identify tools for a better prediction of colon adenocarcinoma prognosis, bolstering the development of a better personalized treatment and management. Method: We used the least absolute shrinkage and selection operator (LASSO) Cox model to analyze the prognosis-related gene datasets from the Gene Expression Omnibus (GEO) database and verified them using The Cancer Genome Atlas (TCGA) database. The area under the curve (AUC) was calculated using the receiver operating characteristic (ROC) curve to evaluate the predictive ability of the risk score model. Gene Set Enrichment Analysis (GSEA) was used to identify the significantly enriched and depleted biological processes. The tumor immune dysfunction and exclusion (TIDE) algorithm was taken to explore the relationship between the risk score and immunotherapy. The observations collectively helped us construct a nomogram to predict prognosis. Finally, the correlation between drug sensitivity and prognostic gene sets was conducted based on the Cancer Therapeutics Response Portal (CTRP) analyses. Results: We constructed a scoring model to assess the significance of the prognosis risk-related gene signatures, which was relative to common tumor characteristics and tumor mutational burdens. Patients with a high-risk score had higher tumor stage and poor prognosis (p< 0.05). Moreover, the expressions of these genes were in correlation with changes in the tumor microenvironment (TME). The risk score is an independent prognostic factor for COAD (p< 0.05). The accuracy of the novel nomogram model with a risk score and TNM-stage prediction prognosis in the predicting prognosis was higher than that of the TNM stage. Further analysis showed that a high-risk score was associated with tumor immune rejection. Patients with a low-risk score have a better prognosis with chemotherapy than those with a high-risk score. Compared to patients in the high-risk group, patients in the low-risk group had a significant survival advantage after receiving chemotherapy. In addition, the prognostic gene sets aid the assessment of drug sensitivity. Conclusion: This study establishes a new prognostic model to better predict the clinical outcome and TME characteristics of colon adenocarcinoma. We believe, our model also serves as a useful clinical tool to strengthen the functioning of chemotherapy, immunotherapy, and targeted drugs.

Construction and Validation of an Oxaliplatin-Resistant Gene Signature in Colorectal Cancer Patients Who Underwent Chemotherapy

Aberrant expression of genes contributes to the chemoresistance of colorectal cancer (CRC) treatment. This study aimed to identify genes associated with the chemoresistance of oxaliplatin-based chemotherapy in CRC patients and to construct a signature. Oxaliplatin resistance-related genes were screened by analyzing the gene profiles of cell lines and tissue samples that underwent oxaliplatin-based treatment. Oxaliplatin resistance-related genes were used to establish a signature. The association of the signature had clinical significance, so the prognostic value of the signature was analyzed. Independent cohorts and CRC cell lines were used to validate the value of the gene signature and the oxaliplatin-resistant genes. There were 64 oxaliplatin resistance-related genes identified after overlapping the genes from the dataset of oxaliplatin-treated CRC cells and the dataset of patients treated with oxaliplatin-based chemotherapy. A gene signature based on five oxaliplatin resistance-related genes was established. This gene signature effectively predicted the prognosis of CRC patients who underwent chemotherapy. No significant associations were found between the gene mutations and survival of the patients; however, two genes were associated with microsatellite instability status. Two external independent cohorts and CRC cell line experiments validated the prognostic values of the signature and expression of the genes after oxaliplatin treatment. In conclusion, the oxaliplatin resistance-related gene signature involving five genes was a novel biomarker for the prediction of the chemotherapy response and prognosis of CRC patients who underwent oxaliplatin-based chemotherapy.

Tanshinone IIA reverses oxaliplatin resistance in colorectal cancer through microRNA-30b-5p/AVEN axis

This research aims to explore the role of Tanshinone IIA (Tan IIA) and microRNA (miR)-30b-5p in chemoresistance of colorectal cancer (CRC). The expression levels of miR-30b-5p and apoptosis and caspase activation inhibitor (AVEN) was detected by reverse transcription-quantitative polymerase chain reaction assay. The cell proliferation and apoptosis were examined by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) and flow cytometry assays. The target relationship between miR-30b-5p and AVEN was confirmed by Dual-luciferase reporter assay. Transwell assay was performed to assess CRC cells’ metastasis. Western blot was carried out to measure the apoptosis-related protein. The results showed that miR-30b-5p was lowly expressed in oxaliplatin-resistance CRC cells SW480 (SW480/R) compared to SW480 cells. Overexpression of miR-30b-5p significantly suppressed the malignant biological behaviors of SW480/R cells and significantly promoted the sensitivity of SW480/R cells to oxaliplatin by down-regulated AVEN expression. Besides, Tan IIA treatment upregulated miR-30b-5p expression in SW480/R cells. Moreover, miR-30b-5p upregulation strengthened the promoting effect of Tan IIA on the sensitivity of SW480/R cells to oxaliplatin. In conclusion, Tan IIA and miR-30b-5p could reverse oxaliplatin resistance of CRC cells and may thus be potential treatment strategies for treating patients with CRC.

Effects of recombinant human adenovirus type 5 combined with transarterial chemoembolization on postoperative metastasis and recurrence of hepatocellular carcinoma patients

BACKGROUND

Postoperative recurrence is currently the main factor affecting the long-term survival of hepatocellular carcinoma (HCC) patients. The folinic acid, fluorouracil, and oxaliplatin (FOLFOX) regimen with transarterial chemoembolization (TACE) is a commonly used postoperative chemotherapy strategy, but its effect is still limited. The aim of this study was to analyze the effects of recombinant human adenovirus type 5 (rhAd5) combined with TACE on postoperative metastasis and recurrence of HCC.

METHODS

Patients with HCC undergoing surgical treatment were collected and divided into the rhAd5 group and control group according to whether rhAd5 was performed. The rhAd5 group was combined with rhAd5 treatment based on TACE. The recurrence and metastasis rates of the two participant groups were compared. The changes of liver function, kidney function, blood routine, and adverse reactions during treatment were analyzed.

RESULTS

The basic data of the two groups were not significantly different (P>0.05). The recurrence and metastasis rates of rhAd5 group participants were significantly lower than those of the control group (P<0.05). There was no significant difference in the incidence of adverse reactions between the rhAd5 group and control group (P>0.05). There were no significant differences in the incidences of adverse reactions between the rhAd5 group and control group (P>0.05).

CONCLUSIONS

The combination of FOLFOX and rhAd5 after surgery can significantly inhibit the occurrence of metastasis and recurrence of HCC patients, improve progression free survival, and has certain safety.

N1-methyladenosine methylation in tRNA drives liver tumourigenesis by regulating cholesterol metabolism

Hepatocellular carcinoma (HCC) accounts for the majority of primary liver cancers and is characterized by high recurrence and heterogeneity, yet its mechanism is not well understood. Here we show that N1-methyladenosine methylation (m1A) in tRNA is remarkably elevated in hepatocellular carcinoma (HCC) patient tumour tissues. Moreover, m1A methylation signals are increased in liver cancer stem cells (CSCs) and are negatively correlated with HCC patient survival. TRMT6 and TRMT61A, forming m1A methyltransferase complex, are highly expressed in advanced HCC tumours and are negatively correlated with HCC survival. TRMT6/TRMT61A-mediated m1A methylation is required for liver tumourigenesis. Mechanistically, TRMT6/TRMT61A elevates the m1A methylation in a subset of tRNA to increase PPARδ translation, which in turn triggers cholesterol synthesis to activate Hedgehog signaling, eventually driving self-renewal of liver CSCs and tumourigenesis. Finally, we identify a potent inhibitor against TRMT6/TRMT61A complex that exerts effective therapeutic effect on liver cancer.

Key Candidate Genes - VSIG2 of Colon Cancer Identified by Weighted Gene Co-Expression Network Analysis

Background

Colon adenocarcinoma (COAD) is one of the most common malignancies. To identify candidate genes that may be involved in colon adenocarcinoma development and progression, weighted gene co-expression network analysis (WGCNA) was used to construct gene co-expression networks to explore associations between gene sets and clinical features and to identify candidate biomarkers. Moreover, we intend to make a preliminary exploration on it.

Methods

Gene expression profiles and clinical information were collected from The Cancer Genome Atlas COAD database for analysis. The gene expression profiles of GSE106582 and GSE110224 were screened from the Gene Expression Omnibus database for verification. WGCNA analysis, functional pathway enrichment analysis, and prognosis analysis were performed on three databases. Target genes were selected from the key genes for experimental verification and research.

Results

Key genes obtained by WGCNA analysis were mainly enriched in key functions and pathways such as drug metabolism, steroid hormones, and retinol metabolism. A total of four prognostic genes were screened out: SELENBP1, NAT2, VSIG2, and CES2. VSIG2 was selected as the target gene for experimental verification, and its encoded protein was found to be mainly expressed in immune cells. Its expression was positively correlated with immune infiltration.

Conclusions

VSIG2 was shown to be associated with immune invasion and antigen presentation in COAD, suggesting it plays an important role in COAD development and progression. It could be used as a potential biomarker or therapeutic target for COAD.

Longitudinal change of genetic variations in cetuximab-treated metastatic colorectal cancer

Recurrent gene mutations and copy number alterations in cancer patients are presumably associated with resistance to targeted therapy. In the present study, we assessed the gene mutations and copy number alterations that recurrently occurred in cetuximab-treated patients with metastatic colorectal cancer (mCRC). Targeted next-generation sequencing was performed in the tumor samples obtained pre- and postcetuximab treatment to assess the variations that occurred during cetuximab treatment. Moreover, we identified the emergent gene mutations (CDK6, EPHA3, ERCC2, MYC, PCMTD1, PIK3CA, PRIM2, RICTOR, and ZNRF3) and copy number alterations (ARAF, BCL2, BRCA2, EGFR, MYC, and SMAD4) that were recurrently observed only in postprogression samples and not in pretreatment or posttreatment samples from patients revealing clinical response. Furthermore, to identify the feasible candidate variations implicated in treatment resistance, we examined the variants with clonal expansion during treatment and discovered PCBP1 as a variant associated with posttreatment progression. Various recurrent mutations were enriched in the TGF-beta signaling pathway. Collectively, we identified recurrent variations in mCRC samples exhibiting post-cetuximab progression. Additionally, future studies are required to evaluate the therapeutic potential of these variations.

Comprehensive genomic analysis contrasting primary colorectal cancer and matched liver metastases

Recent studies have revealed that colorectal cancer (CRC) displays intratumor genetic heterogeneity, and that the cancer microenvironment plays an important role in the proliferation, invasion and metastasis of CRC. The present study performed genomic analysis on paired primary CRC and synchronous colorectal liver metastasis (CRLM) tissues collected from 22 patients using whole-exome sequencing, cancer gene panels and microarray gene expression profiling. In addition, immunohistochemical analysis was used to confirm the protein expression levels of genes identified as highly expressed in CRLM by DNA microarray analysis. The present study identified 10 genes that were highly expressed in CRLM compared with in CRC, from 36,022 probes obtained from primary CRC, CRLM and normal liver tissues by gene expression analysis with DNA microarrays. Of the 10 genes identified, five were classified as encoding 'matricellular proteins' [(osteopontin, periostin, thrombospondin-2, matrix Gla protein (MGP) and glycoprotein nonmetastatic melanoma protein B (GPNMB)] and were selected for immunohistochemical analysis. Osteopontin was strongly expressed in CRLM (6 of 22 cases: 27.3%), but not in CRC (0 of 22: 0%; P=0.02). Periostin also exhibited strong immunoreactivity in CRLM (17 of 22: 68.2%) compared with in CRC (7 of 22: 31.8%; P=0.006). Thrombospondin-2 exhibited strong immunoreactivity in both CRC and CRLM (54.5% in CRC, 45.5% in CRLM; P=0.55). GPNMB and MGP were rarely positive for both CRC and CRLM. A comparison of immunoreactive positive factors for these five genes revealed the complexities of gene expression in CRLM. Of the cases examined, 16 (72.7%) cases of CRC showed zero or only one positive immunoreactive factor. By contrast, CRLM showed more frequent and multiple immunoreactive factors; for example, 16 cases (72.7%) shared two or more factors, which was statistically more frequent than in CRC (P=0.007). The present study revealed the genomic heterogeneity between paired primary CRC and CRLM, in terms of cancer cell microenvironment. This finding may lead to novel diagnostic and therapeutic targets in the era of genome-guided personalized cancer treatment.

Copy number and transcriptome alterations associated with metastatic lesion response to treatment in colorectal cancer

BACKGROUND

Therapeutic resistance is the main cause of death in metastatic colorectal cancer. To investigate genomic plasticity, most specifically of metastatic lesions, associated with response to first-line systemic therapy, we collected longitudinal liver metastatic samples and characterized the copy number aberration (CNA) landscape and its effect on the transcriptome.

METHODS

Liver metastatic biopsies were collected prior to treatment (pre, n = 97) and when clinical imaging demonstrated therapeutic resistance (post, n = 43). CNAs were inferred from whole exome sequencing and were correlated with both the status of the lesion and overall patient progression-free survival (PFS). We used RNA sequencing data from the same sample set to validate aberrations as well as independent datasets to prioritize candidate genes.

RESULTS

We identified a significantly increased frequency gain of a unique CN, in liver metastatic lesions after first-line treatment, on chr18p11.32 harboring 10 genes, including TYMS, which has not been reported in primary tumors (GISTIC method and test of equal proportions, FDR-adjusted p = 0.0023). CNA lesion profiles exhibiting different treatment responses were compared and we detected focal genomic divergences in post-treatment resistant lesions but not in responder lesions (two-tailed Fisher's Exact test, unadjusted p ≤ 0.005). The importance of examining metastatic lesions is highlighted by the fact that 15 out of 18 independently validated CNA regions found to be associated with PFS in this study were only identified in the metastatic lesions and not in the primary tumors.

CONCLUSION

This investigation of genomic-phenotype associations in a large colorectal cancer liver metastases cohort identified novel molecular features associated with treatment response, supporting the clinical importance of collecting metastatic samples in a defined clinical setting.

Clonal Evolution and Timing of Metastatic Colorectal Cancer

Colorectal cancer (CRC) is the third most frequently diagnosed cancer worldwide, where ~50% of patients develop metastasis, despite current improved management. Genomic characterisation of metastatic CRC, and elucidating the effects of therapy on the metastatic process, are essential to help guide precision medicine. Multi-region whole-exome sequencing was performed on 191 sampled tumour regions of patient-matched therapy-naïve and treated CRC primary tumours (n = 92 tumour regions) and metastases (n = 99 tumour regions), in 30 patients. Somatic variants were analysed to define the origin, composition, and timing of seeding in the metastatic progression of therapy-naïve and treated metastatic CRC. High concordance, with few genomic differences, was observed between primary CRC and metastases. Most cases supported a late dissemination model, via either monoclonal or polyclonal seeding. Polyclonal seeding appeared more common in therapy-naïve metastases than in treated metastases. Whereby, treatment prompted for the selection of distinct resistant clones, through monoclonal seeding to distant metastatic sites. Overall, this study reinforces the importance of early clinical detection and surgical excision of the CRC tumour, whilst further highlighting the clinical challenges for metastatic CRC with increased intratumour heterogeneity (either due to early dissemination or polyclonal metastatic spread) and the underlying risk of future therapeutic resistance in treated patients.

Tanshinone IIA Reverses Oxaliplatin Resistance In Human Colorectal Cancer Via Inhibition Of ERK/Akt Signaling Pathway

Background

Oxaliplatin (OXA)-based chemotherapy is generally used to treat human cancers, whereas OXA resistance is a main obstacle for the treatment of colorectal cancer (CRC). Evidence has shown that tanshinone IIA (Tan IIA) could induce apoptosis in CRC cells. However, the role of combination of OXA and Tan IIA on OXA-resistance CRC cells remains unknown. Thus, this study aimed to investigate the effects of Tan IIA in combination with OXA on OXA-resistance CRC cells.

Methods

MTT assay, Ki67 immunofluorescence staining and flow cytometry were used to detect viability, proliferation and apoptosis in OXA-resistant cell line SW480/OXA, respectively. The expressions of Bcl-2, Bax, active caspase 3, p-Akt and p-ERK in SW480/OXA cells were detected with Western blot. In vivo animal study was performed finally.

Results

In this study, the inhibitory effects of OXA on the proliferation and invasion of SW480/OXA cells were significantly enhanced by Tan IIA. In addition, Tan IIA obviously enhanced the anti-apoptosis effects of OXA on SW480/OXA cells via decreasing the levels of Bcl-2, p-Akt and p-ERK, and increasing the levels of Bax and active caspase 3. In vivo experiments confirmed that Tan IIA enhanced OXA sensitivity in SW480/OXA xenograft model.

Conclusion

We found that Tan IIA could reverse OXA resistance in OXA-resistance CRC cells. Therefore, OXA combined with Tan IIA might be considered as a therapeutic approach for the treatment of OXA-resistant CRC.

Identification of key candidate genes for colorectal cancer by bioinformatics analysis

Colorectal cancer (CRC) is one of the most common cancers of the digestive tract. Although numerous studies have been conducted to elucidate the cause of CRC, the exact mechanism of CRC development remains to be determined. To identify candidate genes that may be involved in CRC development and progression, the microarray datasets GSE41657, GSE77953 and GSE113513 were downloaded from the Gene Expression Omnibus database. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes were used for functional enrichment analysis of differentially expressed genes (DEGs). A protein-protein interaction network was constructed, and the hub genes were subjected to module analysis and identification using Search Tool for the Retrieval of Interacting Genes/Proteins and Cytoscape. A total of 142 DEGs were identified, with enriched functions and pathways in the ‘cell cycle’, ‘cell proliferation’, ‘the mitotic cell cycle’ and ‘one-carbon metabolic process’. In addition, 10 hub genes were identified, and functional analysis indicated that these genes are mainly enriched in ‘cell division’, ‘cell cycle’ and functions associated with nucleotide binding processes. Survival analysis demonstrated that DNA topoisomerase II α, cyclin-dependent kinase 1 and CDC28 protein kinase regulatory subunit 2 may be involved in cancer invasion or recurrence. The DEGs identified in the present study may help explain the molecular mechanisms of CRC development and progression.

DNA Copy Number Variations as Markers of Mutagenic Impact

DNA copy number variation (CNV) occurs due to deletion or duplication of DNA segments resulting in a different number of copies of a specific DNA-stretch on homologous chromosomes. Implications of CNVs in evolution and development of different diseases have been demonstrated although contribution of environmental factors, such as mutagens, in the origin of CNVs, is poorly understood. In this review, we summarize current knowledge about mutagen-induced CNVs in human, animal and plant cells. Differences in CNV frequencies induced by radiation and chemical mutagens, distribution of CNVs in the genome, as well as adaptive effects in plants, are discussed. Currently available information concerning impact of mutagens in induction of CNVs in germ cells is presented. Moreover, the potential of CNVs as a new endpoint in mutagenicity test-systems is discussed.