TNFRSF10C copy number variation is associated with metastatic colorectal cancer

Prognostic costimulatory molecule-related signature risk model correlates with immunotherapy response in colon cancer

Costimulatory molecules can promote the activation and proliferation of T cells and play an essential role in immunotherapy. However, their role in the prognosis of colon adenocarcinoma remains elusive. In this study, the expression data of costimulatory molecules and clinicopathological information of 429 patients with colon adenocarcinoma were obtained from The Cancer Genome Atlas database. The patients were divided into training and verification cohorts. Correlation, Cox regression, and Lasso regression analyses were performed to identify costimulatory molecules related to prognosis. After mentioning the construction of the risk mode, a nomogram integrating the clinical characteristics and risk scores of patients was constructed to predict prognosis. Eventually, three prognostic costimulatory molecules were identified and used for constructing a risk model. High expression of these three molecules indicated a poor prognosis. The predictive accuracy of the risk model was verified in the GSE17536 dataset. Subsequently, multivariate regression analysis showed that the signature based on the three costimulatory molecules was an independent risk factor in the training cohort (HR = 2.12; 95% CI = 1.26, 3.56). Based on the risk model and clinicopathological data, the AUC values for predicting the 1-, 3-, and 5-year survival probability of patients with colon adenocarcinoma were 0.77, 0.77, and 0.71, respectively. To the best of our knowledge, this study is the first to report a risk signature constructed based on the costimulatory molecules TNFRSF10c, TNFRSF13c, and TNFRSF11a. This risk signature can serve as a prognostic biomarker for colon adenocarcinoma and is related to the immunotherapeutic response of patients.

Identification of immune-related ferroptosis prognostic marker and in-depth bioinformatics exploration of multi-omics mechanisms in thyroid cancer

Background: Factors such as variations in thyroid carcinoma (THCA) gene characteristics could influence the clinical outcome. Ferroptosis and immunity have been verified to play an essential role in various cancers, and could affect the cancer patients’ prognosis. However, their relationship to the progression and prognosis of many types of THCA remains unclear.

Methods: First, we extracted prognosis-related immune-related genes and ferroptosis-related genes from 2 databases for co-expression analysis to obtain prognosis-related differentially expressed immune-related ferroptosis genes (PR-DE-IRFeGs), and screened BID and CDKN2A for building a prognostic model. Subsequently, multiple validation methods were used to test the model’s performance and compare its performance with other 4 external models. Then, we explored the mechanism of immunity and ferroptosis in the occurrence, development and prognosis of THCA from the perspectives of anti-tumor immunity, CDKN2A-related competitive endogenous RNA regulatory, copy number variations and high frequency gene mutation. Finally, we evaluated this model’s clinical practice value.

Results: BID and CDKN2A were identified as prognostic risk and protective factors, respectively. External data and qRT-PCR experiment also validated their differential expression. The model’s excellent performance has been repeatedly verified and outperformed other models. Risk scores were significantly associated with most immune cells/functions. Risk score/2 PR-DE-IRFeGs expression was strongly associated with BRAF/NRAS/HRAS mutation. Single copy number deletion of CDKN2A is associated with upregulation of CDKN2A expression and worse prognosis. The predicted regulatory network consisting of CYTOR, hsa-miRNA-873-5p and CDKN2A was shown to significantly affect prognosis. The model and corresponding nomogram have been shown to have excellent clinical practice value.

Conclusion: The model can effectively predict the THCA patients’ prognosis and guide clinical treatment. Ferroptosis and immunity may be involved in the THCA’s progression through antitumor immunity and BRAF/NRAS/HRAS mutation. CYTOR-hsa-miRNA-873-5p-CDKN2A regulatory networks and single copy number deletion of CDKN2A may also affect THCA′ progression and prognosis.

Copy Number Alterations as Novel Biomarkers and Therapeutic Targets in Colorectal Cancer

In colorectal cancer, somatic mutations have played an important role as prognostic and predictive biomarkers, with some also functioning as therapeutic targets. Another genetic aberration that has shown significance in colorectal cancer is copy number alterations (CNAs). CNAs occur when a change to the DNA structure propagates gain/amplification or loss/deletion in sections of DNA, which can often lead to changes in protein expression. Multiple techniques have been developed to detect CNAs, including comparative genomic hybridization with microarray, low pass whole genome sequencing, and digital droplet PCR. In this review, we summarize key findings in the literature regarding the role of CNAs in the pathogenesis of colorectal cancer, from adenoma to carcinoma to distant metastasis, and discuss the roles of CNAs as prognostic and predictive biomarkers in colorectal cancer.

Novel Immune-Related Ferroptosis Signature in Esophageal Cancer: An Informatics Exploration of Biological Processes Related to the TMEM161B-AS1/hsa-miR-27a-3p/GCH1 Regulatory Network

Background: Considering the role of immunity and ferroptosis in the invasion, proliferation and treatment of cancer, it is of interest to construct a model of prognostic-related differential expressed immune-related ferroptosis genes (PR-DE-IRFeGs), and explore the ferroptosis-related biological processes in esophageal cancer (ESCA).

Methods: Four ESCA datasets were used to identify three PR-DE-IRFeGs for constructing the prognostic model. Validation of our model was based on analyses of internal and external data sets, and comparisons with past models. With the biological-based enrichment analysis as a guide, exploration for ESCA-related biological processes was undertaken with respect to the immune microenvironment, mutations, competing endogenous RNAs (ceRNA), and copy number variation (CNV). The model’s clinical applicability was measured by nomogram and correlation analysis between risk score and gene expression, and also immune-based and chemotherapeutic sensitivity.

Results: Three PR-DE-IRFeGs (DDIT3, SLC2A3, and GCH1), risk factors for prognosis of ESCA patients, were the basis for constructing the prognostic model. Validation of our model shows a meaningful capability for prognosis prediction. Furthermore, many biological functions and pathways related to immunity and ferroptosis were enriched in the high-risk group, and the role of the TMEM161B-AS1/hsa-miR-27a-3p/GCH1 network in ESCA is supported. Also, the KMT2D mutation is associated with our risk score and SLC2A3 expression. Overall, the prognostic model was associated with treatment sensitivity and levels of gene expression.

Conclusion: A novel, prognostic model was shown to have high predictive value. Biological processes related to immune functions, KMT2D mutation, CNV and the TMEM161B-AS1/hsa-miR-27a-3p/GCH1 network were involved in ESCA progression.

A Genome-Wide Investigation of Effects of Aberrant DNA Methylation on the Usage of Alternative Promoters in Hepatocellular Carcinoma

BACKGROUND

The alternative usage of promoters provides a way to regulate gene expression, has a significant influence on the transcriptome, and contributes to the cellular transformation of cancer. However, the function of alternative promoters (APs) in hepatocellular carcinoma (HCC) has not been systematically studied yet. In addition, the potential mechanism of regulation to the usage of APs remains unclear. DNA methylation, one of the most aberrant epigenetic modifications in cancers, is known to regulate transcriptional activity. Whether DNA methylation regulates the usage of APs needs to be explored. Here, we aim to investigate the effects of DNA methylation on usage of APs in HCC.

METHODS

Promoter activities were calculated based on RNA-seq data. Functional enrichment analysis was implemented to conduct GO terms. Correlation tests were used to detect the correlation between promoter activity and methylation status. The LASSO regression model was used to generate a diagnostic model. Kaplan-Meier analysis was used to compare the overall survival between high and low methylation groups. RNA-seq and whole-genome bisulfite sequencing (WGBS) in HCC samples were performed to validate the correlation of promoter activity and methylation.

RESULTS

We identified 855 APs in total, which could be well used to distinguish cancer from normal samples. The correlation of promoter activity and DNA methylation in APs was observed, and the APs with negative correlation were defined as methylation-regulated APs (mrAPs). Six mrAPs were identified to generate a diagnostic model with good performance (AUC = 0.97). Notably, the majority of mrAPs had CpG sites that could be used to predict clinical outcomes by methylation status. Finally, we verified 85.6% of promoter activity variation and 92.3% of methylation changes in our paired RNA-seq and WGBS samples, respectively. The negative correlation between promoter activity and methylation status was further confirmed in our HCC samples.

CONCLUSION

The aberrant methylation status plays a critical role in the precision usage of APs in HCC, which sheds light on the mechanism of cancer development and provides a new insight into cancer screening and treatment.

CNV Analysis of the Correlation between Preoperative Lymph Node Metastasis and Prognosis of Early Tongue Cancer

Objective: To investigate the relationship between preoperative concealed lymph node metastasis (LNM) and prognosis in patients with early tongue cancer.

Methods: According to preoperative lymph node dissection, 41 patients with tongue cancer were divided into the LNM group (n=19) and the non-LNM group (n=22). Analysis of single nucleotide variation (SNV), tumor mutation burden (TMB), mutant allele tumor heterogeneity (MATH), aneuploidy and copy number variation (CNV) was performed to identify differentially expressed genes (DEGs) related to LNM. While KEGG analysis was conducted to reveal the CNV differentially expressed genes and main enriched pathways, the correlation between pathways and genes was analyzed by hierarchical clustering. The clinical information of LNM and data on overall survival (OS) rate were obtained from The Cancer Genome Atlas (TCGA), and survival analysis was performed based on combined LNM information.

Results: We observed significant correlations for the mTOR signaling pathway, Hippo signaling pathway and Wnt signaling pathway with the LNM group, while leukocyte transendothelial migration and cytokine-cytokine receptor interaction were markedly correlated with the non-LNM group. Moreover, TNFRSF10C was identified as the main DEG in the leukocyte transendothelial migration pathway. TCGA-based validation revealed that the disease-free survival (DFS) of the non-LNM group was significantly higher than that of the LNM group (P<0.005) when TNFRSF10C CNV was set to a log2 ratio>-0.163.

Conclusion: The differential expression of TNFRSF10C in leukocyte transendothelial migration, an immune-related pathway, is associated with LNM and DFS. The TNFRSF10C CNV log2 ratio could potentially serve as an indicator of good prognosis for tongue cancer patients with non-LNM after neck dissection.

New insight into long non-coding RNAs associated with bone metastasis of breast cancer based on an integrated analysis

Background

Bone is the most common site of metastatic breast cancer, and it is a leading cause of breast cancer-related death. This study aimed to explore bone metastasis-related long non-coding RNAs (lncRNAs) in breast cancer.

Methods

Four mRNA datasets and two lncRNA datasets of bone metastasis, lung metastasis and liver metastasis of breast cancer were downloaded from Gene Expression Omnibus (GEO) database. The differentially expressed mRNAs (DEmRNAs) and lncRNAs (DElncRNAs) in group of bone metastasis vs lung metastasis and bone metastasis vs liver metastasis, as well as the overlap of the two groups, were identified. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and protein–protein interaction (PPI) network construction of DEmRNAs were conducted. The cis nearby-targeted DEmRNAs of DElncRNAs were obtained. Quantitative real-time polymerase chain reactions (qRT-PCR) was used to detect the expression levels of selected DEmRNAs and DElncRNAs. LOC641518-lymphoid enhancer-binding factor 1 (LEF1) pair was selected to verify its role in migration and invasion capability of breast cancer cells by wounding healing assay and transwell invasion assay.

Results

A total of 237 DEmRNAs were obtained in bone metastasis compared with both lung metastasis and liver metastasis. A total of three DElncRNAs in bone metastasis compared with both lung metastasis and liver metastasis were obtained. A total of seven DElncRNA-nearby-targeted DEmRNA pairs and 15 DElncRNA-nearby-targeted DEmRNA pairs in group of bone metastasis vs lung metastasis and bone metastasis vs liver metastasis, were detected, respectively. Four cis LncRNA-mRNA interaction pairs were identified, which are LOC641518-LEF1, FLJ35024-Very Low Density Lipoprotein Receptor (VLDLR), LOC285972-Retinoic Acid Receptor Responder 2 (RARRES2) and LOC254896-TNF receptor superfamily member 10c (TNFRSF10C). qRT-PCR using clinical samples from our hospital confirms the bioinformatics prediction. siRNA knocking down LOC641518 down-regulates LEF1 mRNA expression, and reduces the migration and invasion capability of breast cancer cells.

Conclusions

We concluded that four LncRNA-mRNA pairs, including LOC641518-LEF1, may play a central role in breast cancer bone metastasis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02068-7.

Copy number variation in triple negative breast cancer samples associated with lymph node metastasis

Triple negative breast cancer (TNBC) is a highly metastatic and aggressive subtype of breast cancer and cases presenting with lymph node involvement have worse outcomes. This study aimed to determine the regions of copy number variation (CNV) associated with lymph node metastasis in TNBC patients. CNV analyses were performed in a study cohort of 23 invasive ductal carcinomas (IDCs), 12 lymph node metastases (LNmets), and 7 normal adjacent tissues (NATs); as well as in an independent cohort containing 70 TNBC IDCs and the same 7 NATs. CNV-associated genes were analyzed using GO-enrichment and Pathway analysis. The prognostic role for genes showing CNV-based changes in messenger RNA expression was determined using the Kaplan-Meier plotter database. For the IDCs, there were a number of variations that were common in both the study and independent cohorts in the amplified regions of 1q, 8q, 19 (p and q), 2p, 5p and the deleted regions in 8p followed by 5q, and 19p. The most frequently amplified regions in the LNmets of the study cohort were 4q28.3, 2p, 3q24, 1q21.2, 10p, 12p11.1, 8q, 20p11.22-20p11.21, 21q22.13, 6p22.1 and the most frequently deleted regions were in 1p36.23, 4q21.1 and 5q. A total of 686 (441 amplified and 245 deleted) genes were associated with LNmets. The LNmet-associated genes were highly enriched for “regulation of complement activation,” “regulation of protein activation cascade,” “regulation of humoral immune response,” “oxytocin signalling pathway,” and “TRAIL binding” pathways. Moreover, 6/686 LNmet-associated genes showed CNV-based changes in their mRNA expression of which, high expression of ASPM and KIF14 was significantly associated with worse relapse-free survival. This study has identified several CNV regions in TNBC that could play a major role in metastasis to the lymph node.

Inhibition of Autophagy Amplifies Baicalein-Induced Apoptosis in Human Colorectal Cancer

Baicalein is a Chinese herbal compound extracted from Scutellaria baicalensis that has anti-tumor properties. The aim of this study was to elucidate the mechanisms of action of baicalein against human colorectal cancer cell lines and to assess whether the anti-proliferative effects of baicalein may be amplified with autophagy inhibition. Human colon cancer cell lines (HT-29, HCT-116, SW480, and SW620) were treated with baicalein alone and in combination with the autophagy inhibitor chloroquine (CQ). Baicalein reduced cell viability in all four colon cancer lines in a dose-dependent fashion. Combination treatment of baicalein and the autophagy inhibitor CQ significantly decreased cell viability compared with baicalein alone in HT-29 and HCT-116 cell lines. Western blot analysis of the HCT-116 cell line treated with both baicalein and CQ demonstrated increased expression of LC3-II, a component of autophagy. The combination of baicalein with CQ culminated in activation of caspase-3-mediated apoptosis. These findings demonstrate that inhibition of autophagy enhanced apoptotic cell death induced by baicalein treatment in colon cancer cell lines. Future work will assess other targetable apoptotic pathways activated by baicalein and autophagy inhibition.

Development of a risk scoring system for evaluating the prognosis of patients with Her2-positive breast cancer

Background

As one of the many breast cancer subtypes, human epidermal growth factor receptor 2 (Her2)-positive breast cancer has higher invasiveness and poor prognosis, although the advent of anti-Her2 drugs has brought good news to patients. However, the emergence of drug resistance still limits its clinical efficacy, so there is an urgent need to explore new targets and develop a risk scoring system to improve treatments and evaluate patient prognosis.

Methods

Differentially expressed mRNAs associated with Her2-positive breast cancer were screened from a TCGA cohort. The prognostic risk scoring system was constructed according to univariate and Lasso Cox regression model analyses and combined with clinical factors (such as age and TNM) for univariate and multivariate analyses to verify the specificity and sensitivity of the risk scoring system. Finally, based on correlation and CNV mutation analyses, we explored the research value of the mRNAs involved in the system as key genes of the model.

Results

In this study, six mRNAs were screened and identified to construct a prognostic risk scoring system, including four up-regulated mRNA (RDH16, SPC25, SPC24, and SCUBE3) and two down-regulated mRNA (DGAT2 and CCDC69). The risk scoring system can divide Her2-positive breast cancer samples into high-risk and low-risk groups to evaluate patient prognosis. In addition, whether through the time-dependent receiver operating characteristics curve or compared with clinical factors, the risk scoring system showed high predictive sensitivity and specificity. Moreover, some CNV mutations in mRNA increase patient risk by influencing expression levels.

Conclusion

The risk scoring system constructed in this study is helpful to improve the screening of high-risk patients with Her2-positive breast cancer and is beneficial for implementing early diagnosis and personalized treatment. It is suggested that these mRNAs may play an important role in the progression of Her2-positive breast cancer.

Identification and meta-analysis of copy number variation-driven circadian clock genes for colorectal cancer

Both copy number variation (CNV) and circadian clock genes play a critical role in the etiology and pathogenesis of colorectal cancer (CRC); however, a comprehensive analysis of CNV-driven circadian clock genes is urgently required. The present study aimed to investigate the systematic associations between somatic cell CNVs and circadian clock gene expression in patients with CRC. Using somatic CNV, legacy clinical information and gene expression data from The Cancer Genome Atlas, 295 genes that were significantly differentially expressed and with significantly different CNV were obtained, and the expression of the genes, among which 15 were circadian clock genes, was significantly associated with CNV. Further analysis revealed that aryl hydrocarbon receptor nuclear translocator-like 2 (ARNTL2) expression and CNV in these circadian clock genes were significantly associated with survival time in patients with CRC, and the expression of ARNTL2 was also significantly associated with the pathological stage of CRC. Gene set enrichment analysis found that ARNTL2 is enriched for gene sets associated with CRC pathogenesis such as the p53 signaling pathway. These results suggest that ARNTL2 may be a promising prognostic biomarker for patients with CRC, and that circadian clock genes play an important role in CRC through CNV.

Effect of FOLFOX6 chemotherapy on serum VEGF expression in advanced colorectal cancer patients

OBJECTIVE

To explore the effect of FOLFOX6 chemotherapy on serum vascular endothelial growth factor (VEGF) expression in advanced colorectal cancer patients.

METHODS

A retrospective analysis of 81 patients with advanced colorectal cancer who visited our hospital from March 2014 to February 2016 was performed. All the patients were treated with FOLFOX6 chemotherapy. On day 1, patients received oxaliplatin 100 mg/m2 ivgtt (2h), calcium folinate 200 mg/m2 ivgtt (2h), 5 fluorouracil 400 mg/m2 iv bolus and 5 fluorouracil 2500 mg/m2 ivgtt (5h). The treatment course was 2 weeks, and 4 treatment courses were required. The changes in the levels of VEGF and CRP and quality of life before and after 4 courses of chemotherapy were observed and therapeutic effects and adverse reactions after chemotherapy were evaluated.

RESULTS

After treatment, the total efficiency of chemotherapy was 82.72% (67/81) with 24 cases in complete remission, 25 cases in partial response, 18 cases in stable disease and 14 cases in progressive disease. The levels of CRP and VEGF after the treatment were significantly lower than those before treatment (5.69±0.77) mg/L vs. (7.99±1.36) mg/L; (443.26±21.55) pg/mL vs. (542.83±20.44) pg/mL] (P<0.05). The KPS grade after treatment was significantly higher than that before treatment (57.84±4.6) point vs. (50.99±3.73) point] (P<0.05). Among them, 3 cases developed a rash, 5 cases experienced hair loss, and 9 cases developed nausea and vomiting.

CONCLUSION

FOLFOX6 chemotherapy can decrease serum VEGF expression in patients with advanced colorectal cancer and enhance the curative effect with high safety, which is good for the improvement of patients' survival.

TNFRSF10C methylation is a new epigenetic biomarker for colorectal cancer

Background

Abnormal methylation of TNFRSF10C was found to be associated with different types of cancers, excluding colorectal cancer (CRC). In this paper, the performance of TNFRSF10C methylation in CRC was studied in two stages.

Method

The discovery stage was involved with 38 pairs of CRC tumor and paired adjacent non-tumor tissues, and 69 pairs of CRC tumor and paired adjacent non-tumor tissues were used for the validation stage. Quantitative methylation specific PCR (qMSP) method and percentage of methylated reference (PMR) were used to test and represent the methylation level of TNFRSF10C, respectively. A dual-luciferase reporter gene experiment was conducted to evaluate the promoter activity of TNFRSF10C fragment.

Results

A significant association of TNFRSF10C promoter hypermethylation with CRC was found and validated (discovery stage: 24.67 ± 7.52 vs. 3.36 ± 0.89; P = 0.003; validation stage: 31.21 ± 12.48 vs. 4.52 ± 1.47; P = 0.0005). Subsequent analyses of TCGA data among 46 pairs of CRC samples further confirmed our findings (cg23965061: P = 4E - 6; cg14015044: P = 1E - 7). Dual-luciferase reporter gene assay revealed that TNFRSF10C fragment was able to significantly promote gene expression (Fold change = 2.375, P = 0.013). Our data confirmed that TNFRSF10C promoter hypermethylation can predict shorter overall survival of CRC patients (P = 0.032). Additionally, bioinformatics analyses indicated that TNFRSF10C hypermethylation was significantly associated with lower TNFRSF10C expression.

Conclusion

Our work suggested that TNFRSF10C hypermethylation was significantly associated with the risk of CRC.

Studying the mechanism of PLAGL2 overexpression and its carcinogenic characteristics based on 3'-untranslated region in colorectal cancer

Pleomorphic adenoma gene like-2 (PLAGL2) is a zinc finger protein transcription factor, which is upregulated and serves an oncogenic function in multiple human malignancies, including colorectal cancer (CRC). First, reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was used to detect the expression levels of PLAGL2 in CRC tissues and normal tissues. Then, bioinformatics analysis, RT-qPCR, western blotting, luciferase reporter assays and RNA-binding protein immunoprecipitation assays were performed to explore whether the underlying mechanisms, including copy number variation (CNV), microRNAs (miRNAs/miRs) and RNA-binding proteins (RBPs) led to the abnormal expression of PLAGL2. Finally, cell counting kit-8 assays, Transwell assays and xenograft models were used to detect carcinogenesis-associated characteristics based on the 3′-untranslated region (3′-UTR) of PLAGL2. In the present study, PLAGL2 was revealed to be upregulated in CRC tissues compared with normal CRC tissues. CNV was one of the causes leading to the upregulation of PLAGL2. miRNA, including downregulated miR-486-5p, and RBPs, including upregulated human antigen R (HuR), were other key underlying causes. In addition, PLAGL2 3′-UTR was revealed to promote the progression of CRC in vitro and in vivo, and to regulate the expression of C-MYC and CD44. To conclude, these results suggested that high expression of PLAGL2 in CRC was associated with CNV, miR-486-5p and HuR expression, whose 3′-UTR may promote colon carcinogenesis and serve as a novel potential biomarker for CRC therapies.