An EMT‐related gene signature for the prognosis of human bladder cancer

The transition from non–muscle‐invasive bladder cancer (NMIBC) to muscle‐invasive bladder cancer (MIBC) is detrimental to bladder cancer (BLCA) patients. Here, we aimed to study the underlying mechanism of the subtype transition. Gene set variation analysis (GSVA) revealed the epithelial‐mesenchymal transition (EMT) signalling pathway with the most positive correlation in this transition. Then, we built a LASSO Cox regression model of an EMT‐related gene signature in BLCA. The patients with high risk scores had significantly worse overall survival (OS) and disease‐free survival (DFS) than those with low risk scores. The EMT‐related gene signature also performed favourably in the accuracy of prognosis and in the subtype survival analysis. Univariate and multivariate Cox regression analyses demonstrated that the EMT‐related gene signature, pathological N stage and age were independent prognostic factors for predicting survival in BLCA patients. Furthermore, the predictive nomogram model was able to effectively predict the outcome of BLCA patients by appropriately stratifying the risk score. In conclusion, we developed a novel EMT‐related gene signature that has tumour‐promoting effects, acts as a negative independent prognostic factor and might facilitate personalized counselling and treatment in BLCA.     

An outcome model for human bladder cancer: A comprehensive study based on weighted gene co‐expression network analysis

The precision evaluation of prognosis is crucial for clinical treatment decision of bladder cancer (BCa). Therefore, establishing an effective prognostic model for BCa has significant clinical implications. We performed WGCNA and DEG screening to initially identify the candidate genes. The candidate genes were applied to construct a LASSO Cox regression analysis model. The effectiveness and accuracy of the prognostic model were tested by internal/external validation and pan‐cancer validation and time‐dependent ROC. Additionally, a nomogram based on the parameter selected from univariate and multivariate cox regression analysis was constructed. Eight genes were eventually screened out as progression‐related differentially expressed candidates in BCa. LASSO Cox regression analysis identified 3 genes to build up the outcome model in E‐MTAB‐4321 and the outcome model had good performance in predicting patient progress free survival of BCa patients in discovery and test set. Subsequently, another three datasets also have a good predictive value for BCa patients' OS and DFS. Time‐dependent ROC indicated an ideal predictive accuracy of the outcome model. Meanwhile, the nomogram showed a good performance and clinical utility. In addition, the prognostic model also exhibits good performance in pan‐cancer patients. Our outcome model was the first prognosis model for human bladder cancer progression prediction via integrative bioinformatics analysis, which may aid in clinical decision‐making.     

STAT3及下游基因Cyclin D1在膀胱癌中的表达及意义

目的 探讨STAT3及下游基因CyclinD1在膀胱癌中的表达及临床意义.为膀胱癌的基因治疗提供理论依据.方法 收集武汉大学人民医院病理科2004-2006年膀胱移行细胞癌存档蜡块50例,其中男性38例,女性12例,平均年龄50.3岁,另取5例癌旁组织做对照.标本均经HE染色,光镜观察诊断,按WHO病理分级标准,G1:30例,G2:15例,G3:5例;运用免疫组织化学染色方法检测STAT3和CyclinD1在以上各组中的表达.利用HPIAS-2000图像分析系统测定STAT3和CyclinD1在以上各组中表达的平均光密度和平均阳性面积率.结果 1.STAT3的表达:STAT3在膀胱移行细胞癌中呈高表达;STAT3在癌旁组织中呈低表达.膀胱移行细胞癌与癌旁组织比较,STAT3的表达差异有显著性(P＜0.05);2.CyclinD1的表达:CyclinD1在膀胱移行细胞癌中呈高表达;CyclinD1在癌旁组织中呈低表达.膀胱移行细胞癌与癌旁组织比较,CyclinD1的表达差异有显著性(P＜0.05).结论 1.STAT3可能参与了膀胱癌细胞的增殖分化、细胞周期等的调节;2.STAT3的活化促进了抗凋亡基因CyclinD1的表达;3.STAT3及其下游基因CyclinD1在膀胱癌的形成中起了一定作用.     

JAK‐STAT core cancer pathway: An integrative cancer interactome analysis

Through a comprehensive review and in silico analysis of reported data on STAT‐linked diseases, we analysed the communication pathways and interactome of the seven STATs in major cancer categories and proposed rational targeting approaches for therapeutic intervention to disrupt critical pathways and addictions to hyperactive JAK/STAT in neoplastic states. Although all STATs follow a similar molecular activation pathway, STAT1, STAT2, STAT4 and STAT6 exert specific biological profiles associated with a more restricted pattern of activation by cytokines. STAT3 and STAT5A as well as STAT5B have pleiotropic roles in the body and can act as critical oncogenes that promote many processes involved in cancer development. STAT1, STAT3 and STAT5 also possess tumour suppressive action in certain mutational and cancer type context. Here, we demonstrated member‐specific STAT activity in major cancer types. Through systems biology approaches, we found surprising roles for EGFR family members, sex steroid hormone receptor ESR1 interplay with oncogenic STAT function and proposed new drug targeting approaches of oncogenic STAT pathway addiction.     

IGF2BP3 facilitates cell proliferation and tumorigenesis via modulation of JAK/STAT signalling pathway in human bladder cancer

Insulin‐like growth factor‐2 messenger RNA‐binding protein 3 (IGF2BP3) has been reported to contribute to tumorigenesis in several human cancers. However, the biological functions of IGF2BP3 in bladder cancer are poorly understood. We investigated the relation between IGF2BP3 expression and prognosis of bladder cancer patients. Cell proliferation, cell cycle and cell apoptosis assays were performed to assess IGF2BP3 functions. The results showed that IGF2BP3 was overexpressed in bladder cancer tissues compared with that in normal bladder tissues, and its higher expression was closely correlated with poor prognosis in bladder cancer patients. Overexpression of IGF2BP3 markedly promoted cell proliferation and cell cycle progression and inhibited cell apoptosis, while knockdown of IGF2BP3 notably suppressed the proliferation, promoted cell apoptosis and induced cell cycle arrest at the G0/G1 phase. Mechanistically, we revealed that IGF2BP3 promotes the activation of the JAK/STAT pathway in bladder cancer cells. Moreover, the JAK/STAT inhibitor dramatically blocked the tumour‐promoting activity of IGF2BP3. Tumour growth in vivo was also suppressed by knocking down of IGF2BP3. Hence, IGF2BP3 facilitated bladder cancer cell proliferation by activating the JAK/STAT signalling pathway. These findings suggest that IGF2BP3 exhibits an oncogenic effect in human bladder cancer progression.     

Comprehensive analysis of microRNA-messenger RNA regulatory network in gemcitabine-resistant bladder cancer cells

Chemotherapy is still a standard treatment of unresectable bladder cancer or distant metastases. The chemotherapy resistance always occurs after a period of treatment indicating poor prognosis. The current study aimed to explore the molecular mechanism of chemoresistance in bladder cancer cells. The gene expression profiles of GSE77883, including three untreated T24 cells samples and three gemcitabine-resistant T24 cells samples, was downloaded from Gene Expression Omnibus database. The screening of differentially expressed genes (DEGs), gene function analysis, and interaction prediction between microRNAs (miRNAs) and DEGs were performed by R software. The protein-protein interaction (PPI) and miRNA-DEGs networks were constructed and visualized by Cytoscape software. Then, the small molecules, with potential synergistic or antagonistic effects to gemcitabine resistance, were identified using the Connectivity Map database. Finally, gemcitabine-resistant T24 cell line was established and key genes were validated by quantitative real-time polymerase chain reaction (qRT-PCR). In total, 536 upregulated and 513 downregulated genes were screened and mainly enriched in oxidative stress response and signaling pathways related to extracellular matrix–receptor interaction and focal adhesion. PPI network showed interleukin 6, tumor necrosis factor, kinesin family member 11, and BUB1 mitotic checkpoint serine/threonine kinase B were key genes. The miRNA-DEGs regulatory networks included 18 miRNAs and 185 DEGs, including miR-182-5p, miR-590-3p, miR-320a and serum- and glucocorticoid-regulated kinase 1 (SGK1). Then, the related key genes and miRNAs were confirmed by qRT-PCR. Furthermore, 81 small molecules with antagonistic or synergistic effect to GEM were screened. We have investigated the molecular mechanisms driving GEM-resistance in bladder cancer cells that would contribute to the development of chemotherapy for advanced bladder cancer.     

Identification of a prognostic gene signature based on an immunogenomic landscape analysis of bladder cancer

Cancer immune plays a critical role in cancer progression. Tumour immunology and immunotherapy are one of the exciting areas in bladder cancer research. In this study, we aimed to develop an immune‐related gene signature to improve the prognostic prediction of bladder cancer. Firstly, we identified 392 differentially expressed immune‐related genes (IRGs) based on TCGA and ImmPort databases. Functional enrichment analysis revealed that these genes were enriched in inflammatory and immune‐related pathways, including in ‘regulation of signaling receptor activity’, ‘cytokine‐cytokine receptor interaction’ and ‘GPCR ligand binding’. Then, we separated all samples in TCGA data set into the training cohort and the testing cohort in a ratio of 3:1 randomly. Data set {"type":"entrez-geo","attrs":{"text":"GSE13507","term_id":"13507"}}GSE13507 was set as the validation cohort. We constructed a prognostic six‐IRG signature with LASSO Cox regression in the training cohort, including AHNAK, OAS1, APOBEC3H, SCG2, CTSE and KIR2DS4. Six IRGs reflected the microenvironment of bladder cancer, especially immune cell infiltration. The prognostic value of six‐IRG signature was further validated in the testing cohort and the validation cohort. The results of multivariable Cox regression and subgroup analysis revealed that six‐IRG signature was a clinically independent prognostic factor for bladder cancer patients. Further, we constructed a nomogram based on six‐IRG signature and other clinicopathological risk factors, and it performed well in predict patients'' survival. Finally, we found six‐IRG signature showed significant difference in different molecular subtypes of bladder cancer. In conclusions, our research provided a novel immune‐related gene signature to estimate prognosis for patients'' survival with bladder cancer.     

Integrated analysis of ceRNA network and tumor-infiltrating immune cells in esophageal cancer

Background: Esophageal cancer (ESCA) is one of the most commonly diagnosed cancers in the world. Tumor immune microenvironment is closely related to tumor prognosis. The present study aimed at analyzing the competing endogenous RNA (ceRNA) network and tumor-infiltrating immune cells in ESCA.Methods: The expression profiles of mRNAs, lncRNAs, and miRNAs were downloaded from the Cancer Genome Atlas database. A ceRNA network was established based on the differentially expressed RNAs by Cytoscape. CIBERSORT was applied to estimate the proportion of immune cells in ESCA. Prognosis-associated genes and immune cells were applied to establish prognostic models basing on Lasso and multivariate Cox analyses. The survival curves were constructed with Kaplan–Meier method. The predictive efficacy of the prognostic models was evaluated by the receiver operating characteristic (ROC) curves.Results: The differentially expressed mRNAs, lncRNAs, and miRNAs were identified. We constructed the ceRNA network including 23 lncRNAs, 19 miRNAs, and 147 mRNAs. Five key molecules (HMGB3, HOXC8, HSPA1B, KLHL15, and RUNX3) were identified from the ceRNA network and five significant immune cells (plasma cells, T cells follicular helper, monocytes, dendritic cells activated, and neutrophils) were selected via CIBERSORT. The ROC curves based on key genes and significant immune cells all showed good sensitivity (AUC of 3-year survival: 0.739, AUC of 5-year survival: 0.899, AUC of 3-year survival: 0.824, AUC of 5-year survival: 0.876). There was certain correlation between five immune cells and five key molecules.Conclusion: The present study provides an effective bioinformatics basis for exploring the potential biomarkers of ESCA and predicting its prognosis.     

The AP2M1 gene expression is a promising biomarker for predicting survival of patients with hepatocellular carcinoma

There is a growing need for the discovery of new prognostic factors for cases where the scoring and staging system of hepatocellular carcinoma (HCC) does not result in a clear definition. We analyzed whether AP-2 complex subunit mu (AP2M1) expression could be a new prognostic marker for HCC based on the roles of AP2M1 in influencing hepatocyte growth factor (HGF) promoter regulation and hepatitis C virus (HCV) assembly. Patient data were extracted from cohorts of the Gene Expression Omnibus (GSE10186), International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). Differential expression value between matched cancer and normal liver was identified using ICGC cohort. Subsequently, we compared AP2M1 expression as a prognostic gene with other well-known prognostic genes for HCC, using the time-dependent area under the curve (AUC) of the Uno's C-index, the AUC value of the receiver operating characteristics at 5 years, Kaplan-Meier survival curve, and multivariate analysis. Particularly, TCGA and GSE10186 patients were divided into subgroups based on alcohol intake, hepatitis B, and C viral infections, and analyzed in the same methods. The AP2M1 expression values in patients with cancer were much higher than matched normal liver. The AP2M1 level showed excellent prognosis predictions in comparison with existing markers in the three independent cohorts (n = 647). In particular, it was more predictive of prognosis than other markers in alcohol intake and HCV infections. In conclusion, we were confident that AP2M1 provides sufficient value as a new prognostic marker for HCC especially patients with HCV infection and/or alcohol intake.     

Type I IFN induced IL1-Ra expression in hepatocytes is mediated by activating STAT6 through the formation of STAT2: STAT6 heterodimer

The biological activities of type I interferons (IFNs) are mediated by their binding to a heterodimer receptor complex (IFNAR1 and IFNAR2), resulting in the activation of the JAK (JAK1 and TYK2)-STAT (1, 2, 3, 5 isotypes) signalling pathway. Although several studies have indicated that IFN-alpha and IFN-beta can activate complexes containing STAT6, the biological role of this activation is still unknown. We found that exposure of hepatoma cells (HuH7 and Hep3B) to IFN-alpha or IFN-beta led to the activation of STAT6. Activated STAT6 in turn induced the formation of STAT2: STAT6 complexes, which led to the secretion of IL-1Ra. The activation of STAT6 by type I IFN in hepatocytes was mediated by JAK1 and Tyk2. In addition, IFN-alpha or IFN-beta significantly enhanced the stimulatory effect of IL-1beta on production of IL-1Ra. The present study suggests a novel function of IFN-alpha and IFN-beta signalling in human hepatocytes. Our results provide evidence for the mechanism how IFN-alpha and IFN-beta modulate inflammatory responses through activation of STAT6 and production of secreted IL-1Ra.     

An analysis of genetic factors related to risk of inflammatory bowel disease and colon cancer

Background and aimsPatients with inflammatory bowel disease (IBD) have a higher risk of developing colorectal cancer than the general population. Genome-wide association studies have identified and replicated several loci associated with risk of IBD; however, it is currently unknown whether these loci are also associated with colon cancer risk.MethodsWe selected 15 validated SNPs associated with risk of either Crohn's disease, ulcerative colitis, or both in previous GWAS and tested whether these loci were also associated with colon cancer risk in a two-stage study design.ResultsWe found that rs744166 in STAT3 was associated with colon cancer risk in two studies; however, the direction of the observation was reversed in TP53 mutant tumors possibly due to a nullification of the effect by mutant p53. The SNP, which lies within intron 1 of the STAT3 gene, was associated with lower expression of STAT3 mRNA in TP53 wild-type, but not mutant, tumors.ConclusionsThese data suggest that the STAT3 locus is associated with both IBD and cancer. Further understanding the function of this variant in relation to TP53 could possibly explain the role of this gene in autoimmunity and cancer. Furthermore, an analysis of this locus, specifically in a population with IBD, could help to resolve the relationship between this SNP and cancer.     

Role of Jagged1/STAT3 signalling in platinum‐resistant ovarian cancer

Jagged1, the essential ligand of the Notch signalling pathway, is highly expressed in metastatic prostate cancer, and its high expression in breast cancer is linked to poor survival rates. However, the mechanism of Jagged1′s involvement in platinum‐resistant ovarian cancer has not been thoroughly elucidated to date. The purpose of the present study was to investigate the roles of Jagged1 in the platinum resistance of ovarian cancer and its possible mechanisms. Compared with a platinum responsive group of ovarian epithelial cell carcinomas, we found the positive staining intensity of Notch1, Notch2, Jagged1, STAT3 and Epithelial‐mesenchymal transition (EMT) proteins were lower in a platinum‐resistant group. The DDP‐resistant ovarian cancer cell line (C13K) had a higher IC50 of DDP than its parental cell line (OV2008) (P < 0.05) and acquired an EMT phenotype and invasive characteristics. Inhibiting or knockdown of Jagged1 expression could not only reduce its capacity of migration and invasion but also reverse EMT and down‐regulate the expression of serine 727‐phosphorylated STAT3 (pS727) at the protein level but not total STAT3 or tyrosine 705‐phosphorylated STAT3 (pY705) in C13K cells. Furthermore, it was found that crosstalk between the Jagged1/Notch and JAK/STAT3 signalling pathways were involved in Jagged1‐promoting EMT in C13K cells. Experiments in vivo showed a reduced micrometastatic tumour burden in the lung, liver and spleen of mice implanted with C13K cells with knocked‐down Jagged1 compared with mice implanted with control cells. All of these results demonstrate that Jagged1 can crosstalk with the JAK/STAT3 pathway, and they all cooperate to promote the aberrant occurrence of EMT, further reinforcing the abilities of invasion and migration of platinum‐resistant ovarian cancer in vivo and in vitro.     

Identification of functional lncRNAs in gastric cancer by integrative analysis of GEO and TCGA data

Gastric cancer (GC) is a prevalent malignant cancer of digestive system, identification of novel diagnostic and prognostic biomarkers for GC is urgently demanded. The aim of this study was to determine potential long noncoding RNAs (lncRNAs) associated with the pathogenesis and prognosis of GC. Raw noncoding RNA microarray data (GSE53137, GSE70880, and GSE99417) was downloaded from Gene Expression Omnibus (GEO) database. Differentially expressed genes between GC and adjacent normal gastric tissue samples were screened by an integrated analysis of multiple gene expression profile after gene reannotation and batch normalization. Differentially expressed genes were further confirmed by The Cancer Genome Atlas (TCGA) database. Competing endogenous RNA (ceRNA) network, Gene Ontology term and Kyoto Encyclopedia of Genes and Genomes pathway, survival analysis were extensively applied to identify hub lncRNAs and discover potential biomarkers related to diagnosis and prognosis of GC. In total of 246 integrated differential genes including 15 lncRNAs and 241 messenger RNAs (mRNAs) were obtained after intersections of differential genes between GEO and TCGA database. ceRNA network comprised of three lncRNAs (UCA1, HOTTIP, and HMGA1P4), 26 microRNAs (miRNAs) and 72 mRNAs. Functional analysis revealed that three lncRNAs were mainly dominated in cell cycle and cellular senescence. Survival analysis showed that HMGA1P4 was statistically related to the overall survival rate. For the first time, we identified that HMGA1P4, a target of miR-301b/miR-508, is involved in cell cycle and senescence process by regulating CCNA2 in GC. Finally, the expression levels of three lncRNAs were validated to be upregulated in GC tissues. Thus, three lncRNAs including UCA1, HOTTIP, and HMGA1P4 may contribute to GC development and their potential functions might be associated with the prognosis of GC.     

CDO1 promoter methylation is associated with gene silencing and is a prognostic biomarker for biochemical recurrence-free survival in prostate cancer patients

Molecular biomarkers may facilitate the distinction between aggressive and clinically insignificant prostate cancer (PCa), thereby potentially aiding individualized treatment. We analyzed cysteine dioxygenase 1 (CDO1) promoter methylation and mRNA expression in order to evaluate its potential as prognostic biomarker. CDO1 methylation and mRNA expression were determined in cell lines and formalin-fixed paraffin-embedded prostatectomy specimens from a first cohort of 300 PCa patients using methylation-specific qPCR and qRT-PCR. Univariate and multivariate Cox proportional hazards and Kaplan-Meier analyses were performed to evaluate biochemical recurrence (BCR)-free survival. Results were confirmed in an independent second cohort comprising 498 PCa cases. Methylation and mRNA expression data from the second cohort were generated by The Cancer Genome Atlas (TCGA) Research Network by means of Infinium HumanMethylation450 BeadChip and RNASeq. CDO1 was hypermethylated in PCa compared to normal adjacent tissues and benign prostatic hyperplasia (P < 0.001) and was associated with reduced gene expression (ρ = ?0.91, P = 0.005). Using two different methodologies for methylation quantification, high CDO1 methylation as continuous variable was associated with BCR in univariate analysis (first cohort: HR = 1.02, P = 0.002, 95% CI [1.01–1.03]; second cohort: HR = 1.02, P = 0.032, 95% CI [1.00–1.03]) but failed to reach statistical significance in multivariate analysis. CDO1 promoter methylation is involved in gene regulation and is a potential prognostic biomarker for BCR-free survival in PCa patients following radical prostatectomy. Further studies are needed to validate CDO1 methylation assays and to evaluate the clinical utility of CDO1 methylation for the management of PCa.