An expression module of WIPF1-coexpressed genes identifies patients with favorable prognosis in three tumor types

CEBPB-mediated upregulation of SERPINA1 promotes colorectal cancer progression by enhancing STAT3 signaling

Colorectal cancer (CRC) is a highly malignant carcinoma associated with poor prognosis, and metastasis is one of the most common causes of death in CRC. Serpin Family A Member 1 (SERPINA1) is a serine protease inhibitor from the Serpin family. Till now, the function and mechanism of SERPINA1 in CRC progression have not been fully illustrated. We established highly metastatic colorectal cancer cells named as RKO-H and Caco2-H by mice liver metastasis model. By integrative bioinformatic approaches, we analyzed the prognostic value and clinical significance of SERPINA1 in CRC, and predicted potential transcription factors. Colony formation, EDU, MTS, Transwell and wound healing assay were performed to evaluate the biological functions of SERPINA1 in CRC in vitro. Experiments in vivo were conducted to explore the effects of SERPINA1 on liver metastasis of CRC. ChIP and luciferase reporter gene assays were performed to identify the transcriptional regulatory mechanism of SERPINA1 by CEBPB. Our results show that SERPINA1 is highly expressed in CRC and correlated with poor clinical outcomes. SERPINA1 promotes the proliferation, migration by activating STAT3 pathway. Mechanistically, CEBPB binds SERPINA1 gene promoter sequence and promotes the transcription of SERPINA1. SERPINA1 drives CEBPB-induced tumor cell growth and migration via augmenting STAT3 signaling. Our results suggest that SERPINA1 is a potential prognostic marker and may serve as a novel treatment target for CRC.

The role of lncRNA HCG18 in human diseases

A substantial number of long noncoding RNAs (lncRNAs) have been identified as potent regulators of human disease. Human leukocyte antigen complex group 18 (HCG18) is a new type of lncRNA that has recently been proven to play an important role in the occurrence and development of various diseases. Studies have found that abnormal expression of HCG18 is closely related to the clinicopathological characteristics of many diseases. More importantly, HCG18 was also found to promote disease progression by affecting a series of cell biological processes. This article mainly discusses the expression characteristics, clinical characteristics, biological effects and related regulatory mechanisms of HCG18 in different human diseases, providing a scientific theoretical basis for its early clinical application.

WIPF1 promotes gastric cancer progression by regulating PI3K/Akt signaling in a myocardin-dependent manner

Wiskott-Aldrich syndrome protein-interacting protein family member 1 (WIPF1) is associated with malignant tumor progression. However, molecular links between WIPF1 and gastric cancer (GC) remain elusive. The expression of WIPF1 was detected in GC tissues and cells. WIPF1 was overexpressed in GC tissues and cells and high expression of WIPF1 was an independent risk factor for a poor prognosis in patients with GC. Further experiments indicated that WIPF1 promoted the proliferation, invasion, and migration of GC cells in vivo and in vitro. WIPF1-regulated genes were closely related to cell proliferation and migration in GC, and silencing WIPF1 significantly repressed PI3K/AKT signaling pathway activation. WIPF1 was activated by myocardin (MYOCD) translation. Rescue experiments confirmed that MYOCD promotes the proliferation, invasion, and migration of GC cells in a WIPF1-dependent manner and activates the PI3K/AKT signaling pathway. MYOCD may transactivate WIPF1 and facilitate GC cell growth and metastasis by activating the PI3K/AKT signaling pathway.

Elesclomol, a copper-transporting therapeutic agent targeting mitochondria: from discovery to its novel applications

Copper (Cu) is an essential element that is involved in a variety of biochemical processes. Both deficiency and accumulation of Cu are associated with various diseases; and a high amount of accumulated Cu in cells can be fatal. The production of reactive oxygen species (ROS), oxidative stress, and cuproptosis are among the proposed mechanisms of copper toxicity at high concentrations. Elesclomol (ELC) is a mitochondrion-targeting agent discovered for the treatment of solid tumors. In this review, we summarize the synthesis of this drug, its mechanisms of action, and the current status of its applications in the treatment of various diseases such as cancer, tuberculosis, SARS-CoV-2 infection, and other copper-associated disorders. We also provide some detailed information about future directions to improve its clinical performance.

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.

Establishment and external verification of an oxidative stress-related gene signature to predict clinical outcomes and therapeutic responses of colorectal cancer

Objective: Accumulated evidence highlights the biological significance of oxidative stress in tumorigenicity and progression of colorectal cancer (CRC). Our study aimed to establish a reliable oxidative stress-related signature to predict patients' clinical outcomes and therapeutic responses. Methods: Transcriptome profiles and clinical features of CRC patients were retrospectively analyzed from public datasets. LASSO analysis was used to construct an oxidative stress-related signature to predict overall survival, disease-free survival, disease-specific survival, and progression-free survival. Additionally, antitumor immunity, drug sensitivity, signaling pathways, and molecular subtypes were analyzed between different risk subsets through TIP, CIBERSORT, oncoPredict, etc. approaches. The genes in the signature were experimentally verified in the human colorectal mucosal cell line (FHC) along with CRC cell lines (SW-480 and HCT-116) through RT-qPCR or Western blot. Results: An oxidative stress-related signature was established, composed of ACOX1, CPT2, NAT2, NRG1, PPARGC1A, CDKN2A, CRYAB, NGFR, and UCN. The signature displayed an excellent capacity for survival prediction and was linked to worse clinicopathological features. Moreover, the signature correlated with antitumor immunity, drug sensitivity, and CRC-related pathways. Among molecular subtypes, the CSC subtype had the highest risk score. Experiments demonstrated that CDKN2A and UCN were up-regulated and ACOX1, CPT2, NAT2, NRG1, PPARGC1A, CRYAB, and NGFR were down-regulated in CRC than normal cells. In H₂O₂-induced CRC cells, their expression was notably altered. Conclusion: Altogether, our findings constructed an oxidative stress-related signature that can predict survival outcomes and therapeutic response in CRC patients, thus potentially assisting prognosis prediction and adjuvant therapy decisions.

The c-MYC-WDR43 signalling axis promotes chemoresistance and tumour growth in colorectal cancer by inhibiting p53 activity

Oxaliplatin chemoresistance is a major challenge in the clinical treatment of colorectal cancer (CRC), which is one of the most common malignancies worldwide. In this study, we identified the tryptophan-aspartate repeat domain 43 (WDR43) as a potentially critical oncogenic factor in CRC pathogenesis through bioinformatics analysis. It was found that WDR43 is highly expressed in CRC tissues, and WDR43 overexpression is associated with poor prognosis of CRC patients. WDR43 knockdown significantly inhibits cell growth by arresting cell cycle and enhancing the effect of oxaliplatin chemotherapy both in vitro and in vivo. Mechanistically, upon oxaliplatin stimulation, c-MYC promotes the transcriptional regulation and expression of WDR43. WDR43 enhances the ubiquitination of p53 by MDM2 through binding to RPL11, thereby reducing the stability of the p53 protein, which induces proliferation and chemoresistance of CRC cells. Thus, the overexpression of WDR43 promotes CRC progression, and could be a potential therapeutic target of chemoresistance in CRC.

Gene co-expression analyses of health(span) across multiple species

Health(span)-related gene clusters/modules were recently identified based on knowledge about the cross-species genetic basis of health, to interpret transcriptomic datasets describing health-related interventions. However, the cross-species comparison of health-related observations reveals a lot of heterogeneity, not least due to widely varying health(span) definitions and study designs, posing a challenge for the exploration of conserved healthspan modules and, specifically, their transfer across species. To improve the identification and exploration of conserved/transferable healthspan modules, here we apply an established workflow based on gene co-expression network analyses employing GEO/ArrayExpress data for human and animal models, and perform a comprehensive meta-study of the resulting modules related to health(span), yielding a small set of literature backed health(span) candidate genes. For each experiment, WGCNA (weighted gene correlation network analysis) was used to infer modules of genes which correlate in their expression with a 'health phenotype score' and to determine the most-connected (hub) genes (and their interactions) for each such module. After mapping these hub genes to their human orthologs, 12 health(span) genes were identified in at least two species (ACTN3, ANK1, MRPL18, MYL1, PAXIP1, PPP1CA, SCN3B, SDCBP, SKIV2L, TUBG1, TYROBP, WIPF1), for which enrichment analysis by g:profiler found an association with actin filament-based movement and associated organelles, as well as muscular structures. We conclude that a meta-study of hub genes from co-expression network analyses for the complex phenotype health(span), across multiple species, can yield molecular-mechanistic insights and can direct experimentalists to further investigate the contribution of individual genes and their interactions to health(span).

Systematic Pan-Cancer Analysis Identifies CDK1 as an Immunological and Prognostic Biomarker

Cyclin-dependent kinase 1 (CDK1) plays an important role in cancer development, progression, and the overall process of tumorigenesis. However, no pan-cancer analysis has been reported for CDK1, and the predictive role of CDK1 in immune checkpoint inhibitors (ICIs) therapy response remains unexplored. Thus, in this study, we first investigated the potential oncogenic role of CDK1 in 33 tumors by multidimensional bioinformatics analysis based on The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets. Bioinformatic analysis and immunohistochemical experiments confirmed that CDK1 is significantly upregulated in most common cancers and is strongly associated with prognosis. Further analysis indicated that CDK1 may influence tumor immunity mainly by mediating the degree of tumor infiltration of immune-associated cells, and the effect of CDK1 on immunity is diverse across tumor types in tumor microenvironment. CDK1 was also positively correlated with tumor mutational burden (TMB) and microsatellite instability (MSI) in certain cancer types, linking its expression to the assessment of possible treatment response. The results of the pan-cancer analysis study showed that the CDK1 gene was positively associated with the expression of three classes of RNA methylation regulatory proteins, and affects RNA function through multiple mechanisms of action and plays an important role in the posttranscriptional regulation of the tumor microenvironment. These findings shed light on the role of the CDK1 gene in cancer progression and provide information to further study the CDK1 gene as a potential target for pan-cancer.

Repurposing of ingenol mebutate for treating human colorectal cancer by targeting S100 calcium-binding protein A4 (S100A4)

Colorectal cancer (CRC) is the world's second most common cause of cancer-related death. Novel treatments are still urgently needed. S100 calcium-binding protein A4 (S100A4) was demonstrated to be an anticancer therapeutic target. Herein, we found that higher S100A4 expression was associated with a poorer prognosis in publicly available cohorts and a Taiwanese CRC patient cohort. To identify repurposed S100A4 inhibitors, we mined the Connectivity Map (CMap) database for clinical drugs mimicking the S100A4-knockdown gene signature. Ingenol mebutate, derived from the sap of the plant Euphorbia peplus, is approved as a topical treatment for actinic keratosis. The CMap analysis predicted ingenol mebutate as a potent S100A4 inhibitor. Indeed, both messenger RNA and protein levels of S100A4 were attenuated by ingenol mebutate in human CRC cells. In addition, CRC cells with higher S100A4 expressions and/or the wild-type p53 gene were more sensitive to ingenol mebutate, and their migration and invasion were inhibited by ingenol mebutate. Therefore, our results suggest the repurposing of ingenol mebutate for treating CRC by targeting S100A4.

Biomarkers of connective tissue disease-associated interstitial lung disease in bronchoalveolar lavage fluid: A label-free mass spectrometry-based relative quantification study

BACKGROUND

The pathogenesis of connective tissue disease-associated interstitial lung disease (CTD-ILD) is unclear. This study aims to identify differentially expressed proteins (DEPs) in CTD-ILD to determine the potential role of these DEPs that may play in the pathogenesis of CTD-ILD and to offer potential therapeutic targets.

METHODS

Bronchoalveolar lavage fluid (BALF) samples were collected from four patients with CTD-ILD and four patients without CTD-ILD. Label-free mass spectrometry-based relative quantification was used to identify the DEPs. Bioinformatics were used to determine the potential biological processes and signaling pathways associated with these DEPs.

RESULTS

We found 65 upregulated DEPs including SFTPD, CADM1, ACSL4, TSTD1, CD163, LUM, SIGLEC1, CPB2, TGFBI and HGD, and 67 downregulated DEPs including SGSH, WIPF1, SIL1, RAB20, OAS3, GMPR2, PLBD1, DNAJC3, RNASET2 and OAS2. The results of GO functional annotation for the DEPs showed that the DEPS were mainly enriched in the binding, cellular anatomical entity, cellular processes, and biological regulation GO terms. The results of KEGG analyses showed that the pathways most annotated with the DEPs were complement and coagulation cascades, metabolic pathways, pathways in cancer, and PPAR signaling pathway. COG analyses further informed the functions associated with these DEPs, with most focused on signal transduction mechanisms; posttranslational modification, protein turnover, chaperones; intracellular trafficking, secretion, and vesicular transport; amino acid transport and metabolism; and lipid transport and metabolism.

CONCLUSIONS

DEPs identified between patients with vs. without CTD-ILD may play important roles in the development of CTD-ILD and are potential new biomarkers for early diagnosis of CTD-ILD.

Identification of Diagnostic Markers for Breast Cancer Based on Differential Gene Expression and Pathway Network

Background: Breast cancer is the second largest cancer in the world, the incidence of breast cancer continues to rise worldwide, and women's health is seriously threatened. Therefore, it is very important to explore the characteristic changes of breast cancer from the gene level, including the screening of differentially expressed genes and the identification of diagnostic markers. Methods: The gene expression profiles of breast cancer were obtained from the TCGA database. The edgeR R software package was used to screen the differentially expressed genes between breast cancer patients and normal samples. The function and pathway enrichment analysis of these genes revealed significant enrichment of functions and pathways. Next, download these pathways from KEGG website, extract the gene interaction relations, construct the KEGG pathway gene interaction network. The potential diagnostic markers of breast cancer were obtained by combining the differentially expressed genes with the key genes in the network. Finally, these markers were used to construct the diagnostic prediction model of breast cancer, and the predictive ability of the model and the diagnostic ability of the markers were verified by internal and external data. Results: 1060 differentially expressed genes were identified between breast cancer patients and normal controls. Enrichment analysis revealed 28 significantly enriched pathways (p < 0.05). They were downloaded from KEGG website, and the gene interaction relations were extracted to construct the gene interaction network of KEGG pathway, which contained 1277 nodes and 7345 edges. The key nodes with a degree greater than 30 were extracted from the network, containing 154 genes. These 154 key genes shared 23 genes with differentially expressed genes, which serve as potential diagnostic markers for breast cancer. The 23 genes were used as features to construct the SVM classification model, and the model had good predictive ability in both the training dataset and the validation dataset (AUC = 0.960 and 0.907, respectively). Conclusion: This study showed that the difference of gene expression level is important for the diagnosis of breast cancer, and identified 23 breast cancer diagnostic markers, which provides valuable information for clinical diagnosis and basic treatment experiments.

Elesclomol induces copper-dependent ferroptosis in colorectal cancer cells via degradation of ATP7A

Cancer cells reprogram their copper metabolism to adapt to adverse microenvironments, such as oxidative stress. The copper chelator elesclomol has been reported to have considerable anticancer efficacy, but the underlying mechanisms remain largely unknown. In this study, we found that elesclomol-mediated copper overload inhibits colorectal cancer (CRC) both in vitro and in vivo. Elesclomol alone promotes the degradation of the copper transporter copper-transporting ATPase 1 (ATP7A), which retards the proliferation of CRC cells. This property distinguishes it from several other copper chelators. Combinational treatment of elesclomol and copper leads to copper retention within mitochondria due to ATP7A loss, leading to reactive oxygen species accumulation, which in turn promotes the degradation of SLC7A11, thus further enhancing oxidative stress and consequent ferroptosis in CRC cells. This effect accounts for the robust antitumour activity of elesclomol against CRC, which can be reversed by the administration of antioxidants and ferroptosis inhibitors, as well as the overexpression of ATP7A. In summary, our findings indicate that elesclomol-induced copper chelation inhibits CRC by targeting ATP7A and regulating ferroptosis.

Analysis of SYK Gene as a Prognostic Biomarker and Suggested Potential Bioactive Phytochemicals as an Alternative Therapeutic Option for Colorectal Cancer: An In-Silico Pharmaco-Informatics Investigation

BACKGROUND

SYK gene regulates the expression of SYK kinase (Spleen tyrosine kinase), an important non-receptor protein-tyrosine kinase for immunological receptor-mediated signaling, which is also considered a tumor growth metastasis initiator. An onco-informatics analysis was adopted to evaluate the expression and prognostic value of the SYK gene in colorectal cancer (CRC), the third most fatal cancer type; of late, it may be a biomarker as another targeted site for CRC. In addition, identify the potential phytochemicals that may inhibit the overexpression of the SYK kinase protein and minimize the human CRC.

MATERIALS & METHODS

The differential expression of the SYK gene was analyzed using several transcriptomic databases, including Oncomine, UALCAN, GENT2, and GEPIA2. The server cBioPortal was used to analyze the mutations and copy number alterations, whereas GENT2, Gene Expression Profiling Interactive Analysis (GEPIA), Onco-Lnc, and PrognoScan were used to examine the survival rate. The protein-protein interaction network of SYK kinase and its co-expressed genes was conducted via Gene-MANIA. Considering the SYK kinase may be the targeted site, the selected phytochemicals were assessed by molecular docking using PyRx 0.8 packages. Molecular interactions were also observed by following the Ligplot+ version 2.2. YASARA molecular dynamics simulator was applied for the post-validation of the selected phytochemicals.

RESULTS

Our result reveals an increased level of mRNA expression of the SYK gene in colorectal adenocarcinoma (COAD) samples compared to those in normal tissues. A significant methylation level and various genetic alterations recurrence of the SYK gene were analyzed where the fluctuation of the SYK alteration frequency was detected across different CRC studies. As a result, a lower level of SYK expression was related to higher chances of survival. This was evidenced by multiple bioinformatics platforms and web resources, which demonstrated that the SYK gene can be a potential biomarker for CRC. In this study, aromatic phytochemicals, such as kaempferol and glabridin that target the macromolecule (SYK kinase), showed higher stability than the controls, and we have estimated that these bioactive potential phytochemicals might be a useful option for CRC patients after the clinical trial.

CONCLUSIONS

Our onco-informatics investigation suggests that the SYK gene can be a potential prognostic biomarker of CRC. On the contrary, SYK kinase would be a major target, and all selected compounds were validated against the protein using in-silico drug design approaches. Here, more in vitro and in vivo analysis is required for targeting SYK protein in CRC.

Identification of prognostic immune-related gene signature associated with tumor microenvironment of colorectal cancer

Background

The tumor microenvironment (TME) has significantly correlation with tumor occurrence and prognosis. Our study aimed to identify the prognostic immune-related genes (IRGs)in the tumor microenvironment of colorectal cancer (CRC).

Methods

Transcriptome and clinical data of CRC cases were downloaded from TCGA and GEO databases. Stromal score, immune score, and tumor purity were calculated by the ESTIMATE algorithm. Based on the scores, we divided CRC patients from the TCGA database into low and high groups, and the differentially expressed genes (DEGs) were identified. Immune-related genes (IRGs) were selected by venn plots. To explore underlying pathways, protein-protein interaction (PPI) networks and functional enrichment analysis were used. After utilizing LASSO Cox regression analysis, we finally established a multi-IRGs signature for predicting the prognosis of CRC patients. A nomogram consists of the thirteen-IRGs signature and clinical parameters was developed to predict the overall survival (OS). We investigated the association between prognostic validated IRGs and immune infiltrates by TIMER database.

Results

Gene expression profiles and clinical information of 1635 CRC patients were collected from the TCGA and GEO databases. Higher stromal score, immune score and lower tumor purity were observed positive correlation with tumor stage and poor OS. Based on stromal score, immune score and tumor purity, 1517 DEGs, 1296 DEGs, and 1892 DEGs were identified respectively. The 948 IRGs were screened by venn plots. A thirteen-IRGs signature was constructed for predicting survival of CRC patients. Nomogram with a C-index of 0.769 (95%CI, 0.717–0.821) was developed to predict survival of CRC patients by integrating clinical parameters and thirteen-IRGs signature. The AUC for 1-, 3-, and 5-year OS were 0.789, 0.783 and 0.790, respectively. Results from TIMER database revealed that CD1B, GPX3 and IDO1 were significantly related with immune infiltrates.

Conclusions

In this study, we established a novel thirteen immune-related genes signature that may serve as a validated prognostic predictor for CRC patients, thus will be conducive to individualized treatment decisions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08629-3.

Suppression MGP inhibits tumor proliferation and reverses oxaliplatin resistance in colorectal cancer

Matrix Gla protein (MGP), an extracellular matrix protein, has been widely reported to participate in the tumorigenic process and is abnormally expressed in several tumors. However, the role of MGP in colorectal cancer (CRC) remains unknown. Chemotherapy resistance represents a significant limitation in the treatment of CRC. Here, a comprehensive bioinformatics analysis revealed that MGP, which is overexpressed in CRC, might act as one of the critical genes conferring resistance to oxaliplatin (OXA). Furthermore, we found that MGP overexpression in tumor tissue might be correlated with cancer stage and patient prognosis, consistent with the bioinformatics analysis. The upregulation of MGP may act as an independent risk factor for CRC. The knockdown of MGP or inhibition of MGP expression significantly increased the sensitivity of the CRC cell lines to OXA. Suppression of MGP may reverse OXA resistance by upregulating copper transporter 1 (CTR1) and downregulating ATP7A and ATP7B. When used in combination with OXA, the inhibition of MGP reduced cancer cell proliferation, invasion, and migration and increased cell apoptosis in vitro. Suppression of MGP or OXA treatment alone significantly inhibited tumor growth in the CRC mouse model. Additionally, we found that OXA might promote the antitumor immune response in vivo. In summary, our study is the first to provide evidence that MGP expression confers OXA chemotherapy resistance in CRC and provides novel strategies to overcome chemotherapy resistance in CRC.

Comprehensive analysis of prognostic immune-related genes in the tumor microenvironment of colorectal cancer

In this study, we used the ESTIMATE algorithm to analyze clinical data and transcriptome profiles of 1635 colorectal cancer (CRC) samples from the Gene Expression Omnibus and The Cancer Genome Atlas databases and identify prognostic immune-related genes (IRGs). We identified 941 differentially expressed (4 downregulated and 937 upregulated) genes by comparing samples with high and low immune, stromal scores and tumor purity. LASSO Cox regression analyses showed that the risk score based on a ten-IRG signature was an independent prognostic factor in CRC. The nomogram with pathological stages (TNM) and the ten-IRG signature showed a C-index of 0.769 (95% CI, 0.717-0.821), and area under ROC curve values of 0.788, 0.782 and 0.789 for 1-, 3-, and 5-year OS, respectively. TIMER database analysis showed positive correlation between the ten prognostic IRGs and the levels of tumor-infiltrated immune cells, including CD4⁺ and CD8⁺ T cells, macrophages, neutrophils, and dendritic cells. These findings demonstrate that this novel ten-IRG signature correlates with the pathological stages and the levels of multiple tumor-infiltrated immune cell types. This makes the ten-IRG signature a potential prognostic factor for CRC patients.

RNAi Screening Identifies that TEX10 Promotes the Proliferation of Colorectal Cancer Cells by Increasing NF-κB Activation

Colorectal cancer (CRC) has become a predominant cancer worldwide. To understand the process of carcinogenesis, a short hairpin RNA library screening is employed to search for candidate genes that promote proliferation in the CRC cell line HT29. The candidate genes overlap with differentially expressed genes in 32 CRC tumor tissues in the GEO dataset GSE8671. The seventh-ranked testis expressed 10 (TEX10) is upregulated in CRC and its knockdown decreases cell proliferation. The TEX10 high-expression group exhibits worse overall survival (P = 0.003) and progression-free survival (P = 0.001) than the TEX10 low-expression group. TEX10 depletion decreases the growth of CRC cells in vitro and in vivo. Gene set enrichment analysis indicates that the nuclear factor-kappa B pathway is significantly enriched in the genes downregulated by TEX10 knockdown. Mechanistically, TEX10 interacts with RELA and increases its nuclear localization. TEX10 promotes RELA occupancy at gene promoters and regulates the expression of a subset of RELA-targeted genes, including TNFAIP8, SAT1, and IL6ST. Taken together, this study identifies that TEX10 promotes the proliferation of CRC cells in an RELA-dependent manner. In addition, high TEX10 expression is associated with poor prognosis in CRC patients.

Revisiting the Concept of Stress in the Prognosis of Solid Tumors: A Role for Stress Granules Proteins?

Simple Summary

Stress Granules (SGs) were discovered in 1999 and while the first decade of research has focused on some fundamental questions, the field is now investigating their role in human pathogenesis. Since then, evidences of a link between SGs and cancerology are accumulating in vitro and in vivo. In this work we summarized the role of SGs proteins in cancer development and their prognostic values. We find that level of expression of protein involved in SGs formation (and not mRNA level) could serve a prognostic marker in cancer. With this review we strongly suggest that SGs (proteins) could be targets of choice to block cancer development and counteract resistance to improve patients care.

Abstract

Cancer treatments are constantly evolving with new approaches to improve patient outcomes. Despite progresses, too many patients remain refractory to treatment due to either the development of resistance to therapeutic drugs and/or metastasis occurrence. Growing evidence suggests that these two barriers are due to transient survival mechanisms that are similar to those observed during stress response. We review the literature and current available open databases to study the potential role of stress response and, most particularly, the involvement of Stress Granules (proteins) in cancer. We propose that Stress Granule proteins may have prognostic value for patients.

WIP Modulates Oxidative Stress through NRF2/KEAP1 in Glioblastoma Cells

Due to their high metabolic rate, tumor cells produce exacerbated levels of reactive oxygen species that need to be under control. Wiskott-Aldrich syndrome protein (WASP)-interacting protein (WIP) is a scaffold protein with multiple yet poorly understood functions that participates in tumor progression and promotes cancer cell survival. However, its participation in the control of oxidative stress has not been addressed yet. We show that WIP depletion increases the levels of reactive oxygen species and reduces the levels of transcription factor NRF2, the master regulator of redox homeostasis. We found that WIP stabilizes NRF2 by restraining the activity of its main NRF2 repressor, the E3 ligase adapter KEAP1, because the overexpression of a NRF2ΔETGE mutant that is resistant to targeted proteasome degradation by KEAP1 or the knock-down of KEAP1 maintains NRF2 levels in the absence of WIP. Mechanistically, we show that the increased KEAP1 activity in WIP-depleted cells is not due to the protection of KEAP1 from autophagic degradation, but is dependent on the organization of the Actin cytoskeleton, probably through binding between KEAP1 and F-Actin. Our study provides a new role of WIP in maintaining the oxidant tolerance of cancer cells that may have therapeutic implications.

Long noncoding RNA HCG18 up-regulates the expression of WIPF1 and YAP/TAZ by inhibiting miR-141-3p in gastric cancer

Background

Accumulating works show that lncRNAs play critical roles in the development of gastric cancer (GC). LncRNA HLA complex group 18 (HCG18) was implicated in the progression of bladder cancer and glioma, but its role in GC is unknown.

Methods

RT‐PCR was used to detect HCG18 and miR‐141‐3p expression in GC specimen. GC cell lines (AGS and MKN‐28) were exploited as cell model. The biological effect of HCG18 on cancer cells was probed by CCK‐8, colony formation, flow cytometry, Transwell and wound‐healing experiments in vitro, and subcutaneous xenotransplanted tumor model and tail vein injection model in vivo. Interaction between HCG18 and miR‐141‐3p was determined by bioinformatics analysis, RT‐PCR, and luciferase reporter experiments. Downstream gene expression of miR‐141‐3p, including Wiskott–Aldrich syndrome protein interacting protein family member 1 (WIPF1), Yes associated protein 1 (YAP), and tafazzin (TAZ) were detected using Western blot.

Results

HCG18 was markedly up‐regulated in GC specimens, while miR‐141‐3p was markedly down‐regulated. Down‐regulation of HCG18 inhibited viability, migration, and invasion of GC cells, while miR‐141‐3p transfection led to opposite effect. HCG18 could down‐regulate miR‐141‐3p through adsorbing it, and a negative association between HCG18 and miR‐141‐3p was found in GC specimens. HCG18 promoted WIPF1, YAP and TAZ expression, nonetheless, such influence was reversed by co‐transfecting with miR‐141‐3p.

Conclusion

HCG18 was aberrantly up‐regulated in GC tissues, and it indirectly regulated the activity of Hippo signaling through counteracting miR‐141‐3p expression.

Dynamic gene regulatory network reconstruction and analysis based on clinical transcriptomic data of colorectal cancer

Inferring dynamic regulatory networks that rewire at different stages is a reasonable way to understand the mechanisms underlying cancer development. In this study, we reconstruct the stage-specific gene regulatory networks (GRNs) for colorectal cancer to understand dynamic changes of gene regulations along different disease stages. We combined multiple sets of clinical transcriptomic data of colorectal cancer patients and employed a supervised approach to select initial gene set for network construction. We then developed a dynamical system-based optimization method to infer dynamic GRNs by incorporating mutual information-based network sparsification and a dynamic cascade technique into an ordinary differential equations model. Dynamic GRNs at four different stages of colorectal cancer were reconstructed and analyzed. Several important genes were revealed based on the rewiring of the reconstructed GRNs. Our study demonstrated that reconstructing dynamic GRNs based on clinical transcriptomic profiling allows us to detect the dynamic trend of gene regulation as well as reveal critical genes for cancer development which may be important candidates of master regulators for further experimental test.

Prognostic and predictive values of CXCL10 in colorectal cancer

BACKGROUND

The role of CXCL10 in progression and prognosis of colorectal cancer (CRC) has been studied for years, yet results remain controversial.

AIM

This study aims to explore the relationship between CXCL10 and CRC progression and prognosis.

METHODS

We evaluated plasma CXCL10 in CRC patients using ELISA. We also performed a meta-analysis of the associations between CXCL10 and overall survival (OS), disease-free survival (DFS), disease-specific survival (DSS), relapse-free survival (RFS), and clinicopathological features. Finally, correlations between CXCL10 and methylation or immune infiltration were performed using TCGA data.

RESULTS

ELISA analysis showed that CXCL10 was associated with age, red blood cells, blood platelets, and blood urea nitrogen. A separate analysis of 3,763 patients from 24 studies revealed that there were significant associations between low CXCL10 expression and OS (HR 1.25, 95% CI 1.01-1.53), DFS (HR 1.65, 95% CI 1.17-2.34), and RFS (HR 1.43, 95% CI 1.20-1.71) in CRC. Additionally, downregulated CXCL10 expression was significantly correlated with age [odds ratio (OR) 1.31, 95% CI 1.13-1.52], metastasis (OR 1.34, 95% CI 1.11-1.63), recurrence (OR 1.46, 95% CI 1.16-1.83), tumor location (OR 1.88, 95% CI 1.58-2.24), differentiation (OR 0.57, 95% CI 0.35-0.93), microsatellite instability (OR 0.23, 95% CI 0.15-0.35), BRAF mutation (OR 1.62, 95% CI 1.25-2.08), p53 mutation (OR 0.28, 95% CI 0.16-0.47), and CIMP (OR 0.27, 95% CI 0.17-0.43). Furthermore, significant associations were observed between CXCL10 and methylation and immune infiltration.

CONCLUSIONS

The study suggests that CXCL10 might be a potential target for the treatment of CRC.

TRIAL REGISTRATION

NCT03189992. Registered 4 June 2017, https://www.clinicaltrials.gov/ct2/show/study/NCT03189992?term=NCT03189992&rank=1 .

Tumour-infiltrating CD8+ lymphocytes and colorectal cancer recurrence by tumour and nodal stage

BACKGROUND

Intratumoural T-cell infiltrate intensity cortes wrelaith clinical outcome in stage II/III colorectal cancer (CRC). We aimed to determine whether this association varies across this heterogeneous group.

METHODS

We performed a pooled analysis of 1804 CRCs from the QUASAR2 and VICTOR trials. Intratumoural CD8+ and CD3+ densities were quantified by immunohistochemistry in tissue microarray (TMA) cores, and their association with clinical outcome analysed by Cox regression. We validated our results using publicly available gene expression data in a pooled analysis of 1375 CRCs from seven independent series.

RESULTS

In QUASAR2, intratumoural CD8+ was a stronger predictor of CRC recurrence than CD3+ and showed similar discriminative ability to both markers in combination. Pooled multivariable analysis of both trials showed increasing CD8+ density was associated with reduced recurrence risk independent of confounders including DNA mismatch repair deficiency, POLE mutation and chromosomal instability (multivariable hazard ratio [HR] for each two-fold increase = 0.92, 95%CI = 0.87-0.97, P = 3.6 × 10-3). This association was not uniform across risk strata defined by tumour and nodal stage: absent in low-risk (pT3,N0) cases (HR = 1.03, 95%CI = 0.87-1.21, P = 0.75), modest in intermediate-risk (pT4,N0 or pT1-3,N1-2) cases (HR = 0.92, 95%CI = 0.86-1.0, P = 0.046) and strong in high-risk (pT4,N1-2) cases (HR = 0.87, 95%CI = 0.79-0.97, P = 9.4 × 10-3); PINTERACTION = 0.090. Analysis of tumour CD8A expression in the independent validation cohort revealed similar variation in prognostic value across risk strata (PINTERACTION = 0.048).

CONCLUSIONS

The prognostic value of intratumoural CD8+ cell infiltration in stage II/III CRC varies across tumour and nodal risk strata.

OMICfpp: a fuzzy approach for paired RNA-Seq counts

Background

RNA sequencing is a widely used technology for differential expression analysis. However, the RNA-Seq do not provide accurate absolute measurements and the results can be different for each pipeline used. The major problem in statistical analysis of RNA-Seq and in the omics data in general, is the small sample size with respect to the large number of variables. In addition, experimental design must be taken into account and few tools consider it.

Results

We propose OMICfpp, a method for the statistical analysis of RNA-Seq paired design data. First, we obtain a p-value for each case-control pair using a binomial test. These p-values are aggregated using an ordered weighted average (OWA) with a given orness previously chosen. The aggregated p-value from the original data is compared with the aggregated p-value obtained using the same method applied to random pairs. These new pairs are generated using between-pairs and complete randomization distributions. This randomization p-value is used as a raw p-value to test the differential expression of each gene. The OMICfpp method is evaluated using public data sets of 68 sample pairs from patients with colorectal cancer. We validate our results through bibliographic search of the reported genes and using simulated data set. Furthermore, we compared our results with those obtained by the methods edgeR and DESeq2 for paired samples. Finally, we propose new target genes to validate these as gene expression signatures in colorectal cancer. OMICfpp is available at http://www.uv.es/ayala/software/OMICfpp_0.2.tar.gz.

Conclusions

Our study shows that OMICfpp is an accurate method for differential expression analysis in RNA-Seq data with paired design. In addition, we propose the use of randomized p-values pattern graphic as a powerful and robust method to select the target genes for experimental validation.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5496-5) contains supplementary material, which is available to authorized users.

Transposon mutagenesis screen in mice identifies TM9SF2 as a novel colorectal cancer oncogene

New therapeutic targets for advanced colorectal cancer (CRC) are critically needed. Our laboratory recently performed an insertional mutagenesis screen in mice to identify novel CRC driver genes and, thus, potential drug targets. Here, we define Transmembrane 9 Superfamily 2 (TM9SF2) as a novel CRC oncogene. TM9SF2 is an understudied protein, belonging to a well conserved protein family characterized by their nine putative transmembrane domains. Based on our transposon screen we found that TM9SF2 is a candidate progression driver in digestive tract tumors. Analysis of The Cancer Genome Atlas (TCGA) data revealed that approximately 35% of CRC patients have elevated levels of TM9SF2 mRNA, data we validated using an independent set of CRC samples. RNAi silencing of TM9SF2 reduced CRC cell growth in an anchorage-independent manner, a hallmark of cancer. Furthermore, CRISPR/Cas9 knockout of TM9SF2 substantially diminished CRC tumor fitness in vitro and in vivo. Transcriptome analysis of TM9SF2 knockout cells revealed a potential role for TM9SF2 in cell cycle progression, oxidative phosphorylation, and ceramide signaling. Lastly, we report that increased TM9SF2 expression correlates with disease stage and low TM9SF2 expression correlate with a more favorable relapse-free survival. Taken together, this study provides evidence that TM9SF2 is a novel CRC oncogene.

Continuity of transcriptomes among colorectal cancer subtypes based on meta-analysis

BACKGROUND

Previous approaches to defining subtypes of colorectal carcinoma (CRC) and other cancers based on transcriptomes have assumed the existence of discrete subtypes. We analyze gene expression patterns of colorectal tumors from a large number of patients to test this assumption and propose an approach to identify potentially a continuum of subtypes that are present across independent studies and cohorts.

RESULTS

We examine the assumption of discrete CRC subtypes by integrating 18 published gene expression datasets and > 3700 patients, and contrary to previous reports, find no evidence to support the existence of discrete transcriptional subtypes. Using a meta-analysis approach to identify co-expression patterns present in multiple datasets, we identify and define robust, continuously varying subtype scores to represent CRC transcriptomes. The subtype scores are consistent with established subtypes (including microsatellite instability and previously proposed discrete transcriptome subtypes), but better represent overall transcriptional activity than do discrete subtypes. The scores are also better predictors of tumor location, stage, grade, and times of disease-free survival than discrete subtypes. Gene set enrichment analysis reveals that the subtype scores characterize T-cell function, inflammation response, and cyclin-dependent kinase regulation of DNA replication.

CONCLUSIONS

We find no evidence to support discrete subtypes of the CRC transcriptome and instead propose two validated scores to better characterize a continuity of CRC transcriptomes.

Low tumor infiltrating mast cell density confers prognostic benefit and reflects immunoactivation in colorectal cancer

The role of mast cells (MCs) in colorectal cancer (CRC) progression was controversial. Thus, our study was designed to evaluate the prognostic value of MCs as well as their correlation with immune microenvironment. A retrospective cohort of CRC patients of stages I-IV was enrolled in our study. Consecutive patients (854) were divided into training set (427 patients) and validation set (427 patients) randomly. The findings were further validated in a GEO cohort, GSE39582 (556 patients). The mast cell density (MCD) was measured by immunohistochemical staining of tryptase or by CIBERSORT algorithm. Low MCD predicted prolonged overall survival (OS) in training and validation set. Moreover, MCD was identified as an independent prognostic indicator in both sets. Better stratification for CRC prognosis can be achieved by building a MCD based nomogram. The prognostic role of MCD was further validated in GSE39582. In addition, MCD predicted improved survival in stages II and III CRC patients receiving adjuvant chemotherapy (ACT). Multiple immune pathways were enriched in low MCD group while cytokines/chemokines promoting anti-tumor immunity were highly expressed in such group. Furthermore, MCD was negatively correlated with CD8+ T cells infiltration. In conclusion, MCD was identified as an independent prognostic factor, as well as a potential biomarker for ACT benefit in stages II and III CRC. Better stratification of CRC prognosis could be achieved by building a MCD based nomogram. Moreover, immunoactivation in low MCD tumors may contributed to improved prognosis.

WIPF1 antagonizes the tumor suppressive effect of miR-141/200c and is associated with poor survival in patients with PDAC

Background

Aberrant expression of Wiskott–Aldrich syndrome protein interacting protein family member 1 (WIPF1) contributes to the invasion and metastasis of several malignancies. However, the role of WIPF1 in human pancreatic ductal adenocarcinoma (PDAC) remains poorly understood.

Methods

Human pancreatic cancer samples from PDAC patients were collected for methylation analysis. Bioinformatic prediction program and luciferase reporter assay were used to identify microRNAs regulating WIPF1 expression. The association between WIPF1 expression and the overall survival (OS) of patients with PDAC was evaluated by using The Cancer Genome Atlas (TCGA) database. The roles of miR-141/200c and WIPF1 on the invasion and metastasis of PDAC cells were investigated both in vitro and in vivo.

Results

We found that compared to the surrounding non-cancerous tissues, there was significantly increased methylation of miR-200c and miR-141 in human PDAC tissues that resulted in their silencing, whereas the members of the other cluster of miR-200 family including miR-200a, miR-200b and miR-429 were hypomethylated. Our data show that forced expression of miR-141 or miR-200c suppressed invasion and metastasis of PDAC cells both in vitro and in xenograft and identified WIPF1 as a direct target of miR-141 and miR-200c. Both miR-141 and miR-200c inhibit WIPF1 by directly interacting with its 3′-untranslated region. Remarkably, silencing of WIPF1 blocked PDAC growth and metastasis both in vitro and in vivo, whereas forced WIPF1 overexpression antagonized the tumor suppressive effect of miR-141/200c. Additionally, by using TCGA database we showed that high expression of WIPF1 correlated with poor survival in patients with PDAC. Moreover, we show that miR-141 and miR-200c blocked YAP/TAZ expression by suppressing WIPF1.

Conclusions

We have identified WIPF1 as an oncoprotein in PDAC and a direct target of miR-141/miR-200c. We have also defined the miR-141/200c-WIPF1-YAP/TAZ as a novel signaling pathway that is involved in the regulation of the invasion and metastasis of human PDAC cells.

Electronic supplementary material

The online version of this article (10.1186/s13046-018-0848-6) contains supplementary material, which is available to authorized users.

Prolonged survival of patients with colorectal cancer is associated with a higher regucalcin gene expression: Overexpression of regucalcin suppresses the growth of human colorectal carcinoma cells in vitro

Regucalcin plays a crucial role as a regulator of transcriptional signaling activity, and its decreased expression or activity may contribute to the promotion of human carcinogenesis. A higher regucalcin expression in the tumor tissues has been demonstrated to prolong the survival of patients with various types of cancer, including pancreatic cancer, breast cancer, liver cancer and lung adenocarcinoma. The involvement of regucalcin in human colorectal cancer was investigated in the current study. Regucalcin gene expression and the survival data of 62 patients with colorectal cancer were obtained though the Gene Expression Omnibus (GEO) database (GSE12945) for outcome analysis. The data of gene expression revealed that the prolonged survival of patients with colorectal cancer was associated with a higher regucalcin gene expression in tumor tissues. The overexpression of regucalcin suppressed colony formation and proliferation, and induced the death of human colorectal carcinoma RKO cells cultured in a medium containing fetal bovine serum in vitro. Mechanistically, the overexpression of regucalcin induced the G1 and G2/M phase cell cycle arrest of the RKO cells through the suppression of multiple signaling pathways, including Ras, Akt, mitogen-activated protein (MAP) kinase and SAPK/JNK. Of note, the overexpression of regucalcin induced an increase in the levels of the tumor suppressors, p53 and Rb, and the cell cycle inhibitor, p21. Moreover, the levels of the transcription factors, c‑fos, c‑jun, nuclear factor (NF)‑κB p65, β-catenin and signal transducer and activator of transcription 3 (Stat3), were suppressed by the overexpression of regucalcin. On the whole, the findings of this study suggest that regucalcin plays a crucial role as a suppressor in human colorectal cancer, and that the suppressed expression of the regucalcin gene may predispose patients to the promotion of colorectal cancer. The overexpression of regucalcin by gene delivery may thus prove to be a novel therapeutic strategy for colorectal cancer.

WIP/ITSN1 complex is involved in cellular vesicle trafficking and formation of filopodia-like protrusions

WIP (WASP interacting protein) together with N-WASP (neural Wiskott-Aldrich syndrome protein) regulates actin polymerization that is crucial for invadopodia and filopodia formation. Recently, we reported the WIP interaction with ITSN1 which is highly implicated in endo-/exocytosis, apoptosis, mitogenic signaling and cytoskeleton rearrangements. Here we demonstrate that the WIP/ITSN1 complex is involved in the transferrin receptor recycling and partially co-localizes with a marker of the fast recycling endosomes, RAB4. Moreover, ITSN1 recruits WIP to RAB4-positive vesicles upon overexpression. Our data indicate that WIP enhances the interaction of N-WASP with ITSN1 and promotes ITSN1/β-actin association. Moreover, the WIP/ITSN1-L complex facilitates formation of filopodia-like protrusions in MCF-7 cells. Thus, WIP/ITSN1 complex is involved in the cellular vesicle trafficking and actin-dependent membrane processes.

Epigenetically upregulated WIPF1 plays a major role in BRAF V600E-promoted papillary thyroid cancer aggressiveness

How the BRAF V600E mutation promotes the pathogenesis and aggressiveness of papillary thyroid cancer (PTC) is not completely understood. Here we explored a novel mechanism involving WASP interacting protein family member 1 (WIPF1). In PTC tumors, compared with the wild-type BRAF, BRAF V600E was associated with over-expression and hypomethylation of the WIPF1 gene. In thyroid cancer cell lines with wild-type BRAF, WIPF1 expression was robustly upregulated upon introduced expression of BRAF V600E (P=0.03) whereas the opposite was seen upon BRAF knockdown or treatment with BRAF V600E or MEK inhibitors in cells harboring BRAF V600E. Methylation of a functionally critical region of the WIPF1 promoter was decreased by expressing BRAF V600E in cells harboring the wild-type BRAF and increased by BRAF knockdown or treatment with BRAF V600E or MEK inhibitors in cells harboring BRAF V600E mutation. Under-expression and hypermethylation of WIPF1 induced by stable BRAF knockdown was reversed by DNA demethylating agent 5′-azadeoxycytidine. Knockdown of WIPF1 robustly inhibited anchorage-independent colony formation, migration, and invasion of thyroid cancer cells and suppressed xenograft thyroid cancer tumor growth and vascular invasion, mimicking the effects of BRAF knockdown. In human PTC tumors, WIPF1 expression was associated with extrathyroidal invasion (P=0.01) and lymph node metastasis (P=2.64E-05). In summary, BRAF V600E-activated MAP kinase pathway causes hypomethylation and overexpression of WIPF1; WIPF1 then functions like an oncoprotein to robustly promote aggressive cellular and tumor behaviors of PTC. This represents a novel mechanism in BRAF V600E-promoted PTC aggressiveness and identifies WIPF1 as a novel therapeutic target for thyroid cancer.

AQP8 inhibits colorectal cancer growth and metastasis by down-regulating PI3K/AKT signaling and PCDH7 expression

Colorectal cancer (CRC) is one of the most common malignant cancers and the leading cause of cancer-related deaths in worldwide. Although the monoclonal antibody therapy is prescribed for CRC, the metastasis resistant to therapy is the major cause of death of patients with CRC, which indicating the urgent demands for new therapeutic targets discovery. Aquaporin 8 (AQP8) has been identified alter expressed in several cancers including breast cancer, lung cancer and prostatic carcinoma. Our study demonstrated the functional significance of AQP8 in CRC cells growth and metastasis. Over-expression of AQP8 remarkably decreased growth, aggressiveness and colony formation in the CRC SW480 and HT-29 cells. Mechanistically, AQP8 over-expression inhibited tumorigenic phenotype by inactivating PI3K/AKT signaling and inhibiting PCDH7 expression. Furthermore, in vivo studies using nude mice xenograft and metastasis model identified the pivotal role of AQP8 in CRC cells growth and metastasis. Taken together, the present study verifies the vital role of the endogenous AQP8 in colorectal cancer progression.

Prediction of novel target genes and pathways involved in bevacizumab-resistant colorectal cancer

Bevacizumab combined with cytotoxic chemotherapy is the backbone of metastatic colorectal cancer (mCRC) therapy; however, its treatment efficacy is hampered by therapeutic resistance. Therefore, understanding the mechanisms underlying bevacizumab resistance is crucial to increasing the therapeutic efficacy of bevacizumab. The Gene Expression Omnibus (GEO) database (dataset, GSE86525) was used to identify the key genes and pathways involved in bevacizumab-resistant mCRC. The GEO2R web tool was used to identify differentially expressed genes (DEGs). Functional and pathway enrichment analyses of the DEGs were performed using the Database for Annotation, Visualization, and Integrated Discovery(DAVID). Protein–protein interaction (PPI) networks were established using the Search Tool for the Retrieval of Interacting Genes/Proteins database(STRING) and visualized using Cytoscape software. A total of 124 DEGs were obtained, 57 of which upregulated and 67 were downregulated. PPI network analysis showed that seven upregulated genes and nine downregulated genes exhibited high PPI degrees. In the functional enrichment, the DEGs were mainly enriched in negative regulation of phosphate metabolic process and positive regulation of cell cycle process gene ontologies (GOs); the enriched pathways were the phosphoinositide 3-kinase-serine/threonine kinase signaling pathway, bladder cancer, and microRNAs in cancer. Cyclin-dependent kinase inhibitor 1A(CDKN1A), toll-like receptor 4 (TLR4), CD19 molecule (CD19), breast cancer 1, early onset (BRCA1), platelet-derived growth factor subunit A (PDGFA), and matrix metallopeptidase 1 (MMP1) were the DEGs involved in the pathways and the PPIs. The clinical validation of the DEGs in mCRC (TNM clinical stages 3 and 4) revealed that high PDGFA expression levels were associated with poor overall survival, whereas high BRCA1 and MMP1 expression levels were associated with favorable progress free survival(PFS). The identified genes and pathways can be potential targets and predictors of therapeutic resistance and prognosis in bevacizumab-treated patients with mCRC.

miR-196b-5p Regulates Colorectal Cancer Cell Migration and Metastases through Interaction with HOXB7 and GALNT5

Purpose: miR-196b-5p has been previously implicated in malignant transformation; however, its role in colorectal cancer has not been fully explored. In this study, we examine the clinical and biological relevance of miR-196b-5p, and the molecular pathways regulated by miR-196b-5p in colorectal cancer.Experimental Design: miR-196b-5p expression was quantitated by qRT-PCR in 2 independent cohorts composed of 292 patients with colorectal cancer in total, to explore its biomarker potential. Transient and stable gain- and loss-of-function experiments were conducted in a panel of colorectal cancer cell lines and mice, to evaluate the impact of miR-196b-5p on proliferation, chemosensitivity, migration/invasion, and metastases formation in vitro and in vivo The molecular pathways influenced by miR-196b-5p were characterized using whole transcriptome profiling, in silico target prediction tools, luciferase interaction assays, and phenocopy/rescue gene knockdown experiments.Results: Low miR-196b-5p expression was significantly associated with metastases and poor outcomes in 2 independent colorectal cancer patient cohorts (P < 0.05, log-rank test). miR-196b-5p inhibition led to significantly increased colorectal cancer cell migration/invasion and metastases formation in mice, whereas ectopic overexpression showed the opposite phenotype. Molecular profiling and target confirmation identified an interaction between miR-196b-5p and HOXB7 and GALNT5, which in turn regulated colorectal cancer cell migration.Conclusions: The association of low levels of miR-196b-5p and poor prognosis in patients with colorectal cancer can be explained by its influence on cancer cell migration and metastases formation. miR-196b-5p has an impact on colorectal cancer progression pathways through direct interaction with genes involved in cancer cell migration. Clin Cancer Res; 23(17); 5255-66. ©2017 AACR.

Standardising RNA profiling based biomarker application in cancer-The need for robust control of technical variables

Histopathology-based staging of colorectal cancer (CRC) has utility in assessing the prognosis of patient subtypes, but as yet cannot accurately predict individual patient's treatment response. Transcriptomics approaches, using array based or next generation sequencing (NGS) platforms, of formalin fixed paraffin embedded tissue can be harnessed to develop multi-gene biomarkers for predicting both prognosis and treatment response, leading to stratification of treatment. While transcriptomics can shape future biomarker development, currently <1% of published biomarkers become clinically validated tests, often due to poor study design or lack of independent validation. In this review of a large number of CRC transcriptional studies, we identify recurrent sources of technical variability that encompass collection, preservation and storage of malignant tissue, nucleic acid extraction, methods to quantitate RNA transcripts and data analysis pipelines. We propose a series of defined steps for removal of these confounding issues, to ultimately aid in the development of more robust clinical biomarkers.

Tumour suppressor EP300, a modulator of paclitaxel resistance and stemness, is downregulated in metaplastic breast cancer

PURPOSE

We have previously described a novel pathway controlling drug resistance, epithelial-to-mesenchymal transition (EMT) and stemness in breast cancer cells. Upstream in the pathway, three miRs (miR-106b, miR-93 and miR-25) target EP300, a transcriptional activator of E-cadherin. Upregulation of these miRs leads to the downregulation of EP300 and E-cadherin with initiation of an EMT. However, miRs regulate the expression of many genes, and the contribution to EMT by miR targets other than EP300 cannot be ruled out.

METHODS

We used lentiviruses expressing EP300-targeting shRNA to downregulate its expression in MCF-7 cells as well as an EP300-knocked-out colon carcinoma cell line. An EP300-expression plasmid was used to upregulate its expression in basal-like CAL51 and MDA-MB-231 breast cancer cells. Drug resistance was determined by short-term proliferation and long-term colony formation assays. Stemness was determined by tumour sphere formation in both soft agar and liquid cultures as well as by the expression of CD44/CD24/ALDH markers. Gene expression microarray analysis was performed in MCF-7 cells lacking EP300. EP300 expression was analysed by immunohistochemistry in 17 samples of metaplastic breast cancer.

RESULTS

Cells lacking EP300 became more resistant to paclitaxel whereas EP300 overexpression increased their sensitivity to the drug. Expression of cancer stem cell markers, as well as tumour sphere formation, was also increased in EP300-depleted cells, and was diminished in EP300-overexpressing cells. The EP300-regulated gene signature highlighted genes associated with adhesion (CEACAM5), cytoskeletal remodelling (CAPN9), stemness (ABCG2), apoptosis (BCL2) and metastasis (TGFB2). Some genes in this signature were also validated in a previously generated EP300-depleted model of breast cancer using minimally transformed mammary epithelial cells. Importantly, two key genes in apoptosis and stemness, BCL2 and ABCG2, were also upregulated in EP300-knockout colon carcinoma cells and their paclitaxel-resistant derivatives. Immunohistochemical analysis demonstrated that EP300 expression was low in metaplastic breast cancer, a rare, but aggressive form of the disease with poor prognosis that is characterized by morphological and physiological features of EMT.

CONCLUSIONS

EP300 plays a major role in the reprogramming events, leading to a more malignant phenotype with the acquisition of drug resistance and cell plasticity, a characteristic of metaplastic breast cancer.

PML promotes metastasis of triple-negative breast cancer through transcriptional regulation of HIF1A target genes

Elucidating the molecular basis of tumor metastasis is pivotal for eradicating cancer-related mortality. Triple-negative breast cancer (TNBC) encompasses a class of aggressive tumors characterized by high rates of recurrence and metastasis, as well as poor overall survival. Here, we find that the promyelocytic leukemia protein PML exerts a prometastatic function in TNBC that can be targeted by arsenic trioxide. We found that, in TNBC patients, constitutive HIF1A activity induces high expression of PML, along with a number of HIF1A target genes that promote metastasis at multiple levels. Intriguingly, PML controls the expression of these genes by binding to their regulatory regions along with HIF1A. This mechanism is specific to TNBC cells and does not occur in other subtypes of breast cancer where PML and prometastatic HIF1A target genes are underexpressed. As a consequence, PML promotes cell migration, invasion, and metastasis in TNBC cell and mouse models. Notably, pharmacological inhibition of PML with arsenic trioxide, a PML-degrading agent used to treat promyelocytic leukemia patients, delays tumor growth, impairs TNBC metastasis, and cooperates with chemotherapy by preventing metastatic dissemination. In conclusion, we report identification of a prometastatic pathway in TNBC and suggest clinical development toward the use of arsenic trioxide for TNBC patients.

NQO1 inhibits proteasome-mediated degradation of HIF-1α

Overexpression of NQO1 is associated with poor prognosis in human cancers including breast, colon, cervix, lung and pancreas. Yet, the molecular mechanisms underlying the pro-tumorigenic capacities of NQO1 have not been fully elucidated. Here we show a previously undescribed function for NQO1 in stabilizing HIF-1α, a master transcription factor of oxygen homeostasis that has been implicated in the survival, proliferation and malignant progression of cancers. We demonstrate that NQO1 directly binds to the oxygen-dependent domain of HIF-1α and inhibits the proteasome-mediated degradation of HIF-1α by preventing PHDs from interacting with HIF-1α. NQO1 knockdown in human colorectal and breast cancer cell lines suppresses HIF-1 signalling and tumour growth. Consistent with this pro-tumorigenic function for NQO1, high NQO1 expression levels correlate with increased HIF-1α expression and poor colorectal cancer patient survival. These results collectively reveal a function of NQO1 in the oxygen-sensing mechanism that regulates HIF-1α stability in cancers.

WIP promotes in-vitro invasion ability, anchorage independent growth and EMT progression of A549 lung adenocarcinoma cells by regulating RhoA levels

Cancer cell migration and invasion involves actin cytoskeleton reorganization, which is regulated by the WASP (Wiskott Aldrich Syndrome Protein) family of proteins such as WASP, N-WASP (Neural-WASP) and WASP interacting protein (WIP). In this study, we found that the expression of WIP was significantly upregulated in metastatic A5-RT3 cells compared to its parental non-tumorigenic HaCaT cells. Using A549 human lung adenocarcinoma cell line as the model system, we found that WIP regulates cell invasion, proliferation and anchorage-independent growth. Expression of WIP was enhanced during TGF-β1 induced epithelial-mesenchymal transition (EMT) and overexpression of WIP accelerated EMT while knocking down WIP attenuated EMT associated morphological changes. Knocking down WIP expression in A549 cells significantly reduced RhoA levels and WIP was found to interact with RhoA suggesting that WIP might be executing its function by regulating RhoA. Acquisition of invasive, proliferative properties and anoikis resistance is the central step in metastasis indicating a novel function of WIP in cancer progression.

Epigenetic silencing of NTSR1 is associated with lateral and noninvasive growth of colorectal tumors

Our aim was to identify DNA methylation changes associated with the growth pattern and invasiveness of colorectal cancers (CRCs). Comparison of the methylation statuses of large (≥20 mm in diameter along the colonic surface) noninvasive tumors (NTs) and small (<20 mm in diameter along the colonic surface) invasive tumors (ITs) using CpG island microarray analysis showed neurotensin receptor 1 (NTSR1) to be hypermethylated in large NTs. Quantitative bisulfite pyrosequencing revealed that NTSR1 is frequently methylated in colorectal tumors, with large NTs exhibiting the highest methylation levels. The higher NTSR1 methylation levels were associated with better prognoses. By contrast, NTSR1 copy number gains were most frequent among small ITs. Methylation of NTSR1 was associated with the gene's silencing in CRC cell lines, whereas ectopic expression of NTSR1 promoted proliferation and invasion by CRC cells. Analysis of primary tumors composed of adenomatous and malignant portions revealed that NTSR1 is frequently methylated in the adenomatous portion, while methylation levels are generally lower in the cancerous portions. These results suggest that NTSR1 methylation is associated with lateral and noninvasive growth of colorectal tumors, while low levels of methylation may contribute to the malignant potential through activation of NTSR1. Our data also indicate that NTSR1 methylation may be a prognostic biomarker in CRC.

Lymphoma Secondary to Congenital and Acquired Immunodeficiency Syndromes at a Turkish Pediatric Oncology Center

The prevalence of lymphoma in primary immunodeficiency cases and autoimmune diseases, as well as on a background of immunodeficiency following organ transplants, is increasing. The lymphoma treatment success rate is known to be a low prognosis. Our study aimed to emphasize the low survival rates in immunodeficient vs. immunocompetent lymphoma patients and also to investigate the effect of rituximab in patients with ataxia telangiectasia and other immunodeficiencies. We summarized the clinical characteristics and treatment results of 17 cases with primary immunodeficiency that developed non-Hodgkin lymphoma (NHL) and Hodgkin lymphoma (HL) retrospectively. Seven patients were diagnosed with ataxia-telangiectasia, two with common variable immunodeficiency, two with selective IgA deficiency, one with X-related lymphoproliferative syndrome, one with Wiskott-Aldrich syndrome, one with Epstein-Barr virus-related lymphoproliferative syndrome, one with interleukin-2-inducible T-cell kinase (ITK) deficiency, and one with lymphoma developing after autoimmune lymphoproliferative syndrome (ALPS). One patient underwent a renal transplant. Of the nine males and eight females (aged 3-12 years, median = 7) that developed lymphoma, seven were diagnosed with HL and ten with NHL (seven B-cell, three T-cell). The NHL patients were started on the Berlin-Frankfurt-Münster, POG9317, LMB-96, or R-CHOP treatment protocols with reduced chemotherapy dosages. HL cases were started on the doxorubicin, bleomycin, vinblastine, and dacarbazine (ABVD) and/or cyclophosphamide, vincristine, procarbazine, and prednisone (COPP) protocol, also with modified dosages. Importantly, all seven cases of HL are alive and in remission, while six of the ten NHL patients have died. Primary immunodeficiency is a strong predisposing factor for developing lymphoma. Low treatment success rates relative to other lymphomas and difficulties encountered during treatment indicate that new treatment agents are needed. While some success has been achieved by combining rituximab with lymphoma treatment protocols in B-NHL cases with primary immunodeficiency, the need for new treatment approaches for these patients remains critical.

Gene signatures of estrogen and progesterone receptor pathways predict the prognosis of colorectal cancer

The associations of estrogen receptor (ER) and progesterone receptor (PR) pathways with the prognosis of colorectal cancer (CRC) are still controversial. The aim of this study was to readdress these issues by introducing a gene signature-based approach to semiquantitate pathway activity. In this approach, the ER and PR pathway activities in CRC were computed based on the expression profiles of the signature genes of ER and PR pathways, respectively. The results showed that the ER pathway activity was progressively significantly decreased from normal colorectal mucosa, colorectal adenoma to CRC. ER pathway signaling was a favorable factor for the presence of microsatellite stability (MSS) in CRC in seven cohorts tested, while was an unfavorable factor for cancer recurrence in all four CRC cohorts tested (n = 1122; overall HR: 0.311, 95% CI: 0.199-0.488, P < 0.001). Subset stratification in stage II patients showed that ER pathway remained significantly inversely associated with recurrence. PR pathway was also suppressed in colorectal tumors and inversely associated with recurrence of CRC, but to a much lesser extent than ER pathway. Moreover, the inverse association of PR pathway with cancer recurrence was more likely observed in CRC with high ER pathway activity, suggesting the interactions between the two pathways. PR pathway was not associated with MSS in CRC, but it was more significant than ER pathway associated with advance cancer stages and cancer response to adjuvant chemotherapy. These results suggested the potential application of the gene signatures of ER and PR pathways, especially the former, as novel markers for prognosis and management of CRC.

WIP and WICH/WIRE co-ordinately control invadopodium formation and maturation in human breast cancer cell invasion

Cancer cells form actin-rich degradative protrusions (invasive pseudopods and invadopodia), which allows their efficient dispersal during metastasis. Using biochemical and advanced imaging approaches, we demonstrate that the N-WASP-interactors WIP and WICH/WIRE play non-redundant roles in cancer cell invasion. WIP interacts with N-WASP and cortactin and is essential for invadopodium assembly, whereas WICH/WIRE regulates N-WASP activation to control invadopodium maturation and degradative activity. Our data also show that Nck interaction with WIP and WICH/WIRE modulates invadopodium maturation; changes in WIP and WICH/WIRE levels induce differential distribution of Nck. We show that WIP can replace WICH/WIRE functions and that elevated WIP levels correlate with high invasiveness. These findings identify a role for WICH/WIRE in invasiveness and highlight WIP as a hub for signaling molecule recruitment during invadopodium generation and cancer progression, as well as a potential diagnostic biomarker and an optimal target for therapeutic approaches.

ColoFinder: a prognostic 9-gene signature improves prognosis for 871 stage II and III colorectal cancer patients

Colorectal cancer (CRC) is a heterogeneous disease with a high mortality rate and is still lacking an effective treatment. Our goal is to develop a robust prognosis model for predicting the prognosis in CRC patients. In this study, 871 stage II and III CRC samples were collected from six gene expression profilings. ColoFinder was developed using a 9-gene signature based Random Survival Forest (RSF) prognosis model. The 9-gene signature recurrence score was derived with a 5-fold cross validation to test the association with relapse-free survival, and the value of AUC was gained with 0.87 in GSE39582(95% CI [0.83-0.91]). The low-risk group had a significantly better relapse-free survival (HR, 14.8; 95% CI [8.17-26.8]; P < 0.001) than the high-risk group. We also found that the 9-gene signature recurrence score contributed more information about recurrence than standard clinical and pathological variables in univariate and multivariate Cox analyses when applied to GSE17536(p = 0.03 and p = 0.01 respectively). Furthermore, ColoFinder improved the predictive ability and better stratified the risk subgroups when applied to CRC gene expression datasets GSE14333, GSE17537, GSE12945and GSE24551. In summary, ColoFinder significantly improves the risk assessment in stage II and III CRC patients. The 9-gene prognostic classifier informs patient prognosis and treatment response.

Structure-Function Analysis of the Mcl-1 Protein Identifies a Novel Senescence-regulating Domain

Unlike other antiapoptotic Bcl-2 family members, Mcl-1 also mediates resistance to cancer therapy by uniquely inhibiting chemotherapy-induced senescence (CIS). In general, Bcl-2 family members regulate apoptosis at the level of the mitochondria through a common prosurvival binding groove. Through mutagenesis, we determined that Mcl-1 can inhibit CIS even in the absence of its apoptotically important mitochondrion-localizing domains. This finding prompted us to generate a series of Mcl-1 deletion mutants from both the N and C termini of the protein, including one that contained a deletion of all of the Bcl-2 homology domains, none of which impacted anti-CIS capabilities. Through subsequent structure-function analyses of Mcl-1, we identified a previously uncharacterized loop domain responsible for the anti-CIS activity of Mcl-1. The importance of the loop domain was confirmed in multiple tumor types, two in vivo models of senescence, and by demonstrating that a peptide mimetic of the loop domain can effectively inhibit the anti-CIS function of Mcl-1. The results from our studies appear to be highly translatable because we discerned an inverse relationship between the expression of Mcl-1 and of various senescence markers in cancerous human tissues. In summary, our findings regarding the unique structural properties of Mcl-1 provide new approaches for targeted cancer therapy.

Overexpression of Transcription Termination Factor 1 is Associated with a Poor Prognosis in Patients with Colorectal Cancer

BACKGROUND

RNA polymerase 1 transcription termination factor (TTF1) mediates the transcription of ribosomal RNA (rRNA). In the current study, we investigated the clinical and biological significance of the TTF1 gene in colorectal cancer (CRC).

METHODS

The expression of TTF1 messenger RNA (mRNA) in tumor and normal tissues from 136 patients with CRC was examined by quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR). We also performed in vitro cell proliferation and migration assays in TTF1-expressing CRC cells. The biological role of TTF1 in CRC was further elucidated using gene set enrichment analysis (GSEA) with CRC samples.

RESULTS

TTF1 expression was significantly higher in tumor tissues than in corresponding normal tissues (p = 0.016). In clinicopathological analysis, the high-TTF1 expression group showed a higher incidence of liver metastasis and lymphatic invasion than the low-TTF1 expression group (p < 0.05), and tended to have more frequent venous invasion than the low-TTF1 expression group. Furthermore, the high-TTF1 expression group had a significantly poorer prognosis than the low-TTF1 expression group (p = 0.011). Moreover, overexpression of TTF1 enhanced the proliferation and migration capacity of CRC cells in vitro. GSEA revealed that TTF1 was significantly associated with the RAS and MYC pathways, and this observation was confirmed in samples from 136 patients with CRC.

CONCLUSION

TTF1 was involved in cancer progression via the RAS and MYC pathways in CRC, suggesting that TTF1 may be a prognostic indicator and therapeutic target in CRC.

RasGRP1 opposes proliferative EGFR-SOS1-Ras signals and restricts intestinal epithelial cell growth

The character of EGFR signals can influence cell fate but mechanistic insights into intestinal EGFR-Ras signalling are limited. Here we show that two distinct Ras nucleotide exchange factors, RasGRP1 and SOS1, lie downstream of EGFR but act in functional opposition. RasGRP1 is expressed in intestinal crypts where it restricts epithelial growth. High RasGRP1 expression in colorectal cancer (CRC) patient samples correlates with a better clinical outcome. Biochemically, we find that RasGRP1 creates a negative feedback loop that limits proliferative EGFR-SOS1-Ras signals in CRC cells. Genetic Rasgrp1 depletion from mice with either an activating mutation in KRas or with aberrant Wnt signalling due to a mutation in Apc resulted in both cases in exacerbated Ras-ERK signalling and cell proliferation. The unexpected opposing cell biological effects of EGFR-RasGRP1 and EGFR-SOS1 signals in the same cell shed light on the intricacy of EGFR-Ras signalling in normal epithelium and carcinoma.