Identification of biomarkers associated with progression and prognosis in bladder cancer via co-expression analysis

COL6A1 expression as a potential prognostic biomarker for risk stratification of T1 high grade bladder cancer: Unveiling the aggressive nature of a distinct non-muscle invasive subtype

PURPOSE

T1 high grade (T1HG) bladder cancer (BC) is a type of non-muscle invasive BC (NMIBC) that is recognized as an aggressive subtype with a heightened propensity for progression. Current risk stratification methods for NMIBC rely on clinicopathological indicators; however, these approaches do not adequately capture the aggressive nature of T1HG BC. Thus, new, more accurate biomarkers for T1HG risk stratification are needed. Here, we enrolled three different patient cohorts and investigated expression of collagen type VI alpha 1 (COL6A1), a key component of the extracellular matrix, at different stages and grades of BC, with a specific focus on T1HG BC.

MATERIALS AND METHODS

Samples from 298 BC patients were subjected to RNA sequencing and real-time polymerase chain reaction.

RESULTS

We found that T1HG BC and muscle invasive BC (MIBC) exhibited comparable expression of COL6A1, which was significantly higher than that by other NMIBC subtypes. In particular, T1HG patients who later progressed to MIBC had considerably higher expression of COL6A1 than Ta, T1 low grade patients, and patients that did not progress, highlighting the aggressive nature and higher risk of progression associated with T1HG BC. Moreover, Cox and Kaplan-Meier survival analyses revealed a significant association between elevated expression of COL6A1 and poor progression-free survival of T1HG BC patients (multivariate Cox hazard ratio, 16.812; 95% confidence interval, 3.283-86.095; p=0.001 and p=0.0002 [log-rank test]).

CONCLUSIONS

These findings suggest that COL6A1 may be a promising biomarker for risk stratification of T1HG BC, offering valuable insight into disease prognosis and guidance of personalized treatment decisions.

Molecular mechanisms of Thrombospondin-2 modulates tumor vasculogenic mimicry by PI3K/AKT/mTOR signaling pathway

Vasculogenic mimicry (VM) differs from the classical tumor angiogenesis model. VM does not depend on endothelial cells; instead, highly aggressive tumor cells mimic endothelial cells to form a vascular-like channel structure. VM mediated by tumor cells is significantly and positively associated with a poor prognosis and low survival rates in patients with highly aggressive cancer. In the treatment of highly aggressive malignancies, the presence of VM is considered an important reason for the unsatisfactory clinical efficacy of anti-tumor-angiogenesis therapy (e.g., therapy targeting vascular endothelial growth factor A). Many targeted therapeutic drugs based on traditional tumor blood vessels have been used clinically. Although some progress has been made in certain tumors, problems such as drug resistance have restricted the expected therapeutic effects. Thrombospondin 2 (THBS2) is one of the most important genes associated with angiogenesis, and this gene exerts angiogenesis-related functions through the PI3K/AKT signaling pathway. Although the PI3K/AKT/mTOR signaling pathway is closely related to the progression of VM, the mechanism by which the promising biomarker THBS2 participates in and regulates tumor VM by activating the PI3K/AKT/mTOR signaling pathway is unclear. In this review, we analyze the monomer structure and biological activity of THBS2, the structure and potential synthesis mechanisms of VM, and the complex mechanisms between THBS2, the PI3K/AKT/mTOR signaling pathway, and VM.

CCNB1 is involved in bladder cancer pathogenesis and silencing CCNB1 decelerates tumor growth and improves prognosis of bladder cancer

In search of an effective therapeutic target for bladder urothelial carcinoma (BLCA), the present study aimed to investigate the expression of cyclin B1 (CCNB1) and its putative mechanism in BLCA. BLCA sequencing data from Gene Expression Omnibus and The Cancer Genome Atlas were used to analyze expression of CCNB1 mRNA and high CCNB1 expression had a poorer prognosis compared with those with low expression. Immunohistochemistry (IHC) samples collected from the Human Protein Atlas database were analyzed for CCNB1 protein expression. Short hairpin (sh) CCNB1-transfected BLCA T24 and 5637 cells were used to investigate the effects of CCNB1 and inhibit the proliferation, migration and invasion of BLCA cells, affect the cell cycle distribution and promote apoptosis of 5637 cells. A sh-CCNB1 BLCA chicken embryo chorioallantoic membrane (CAM) transplantation model was established to observe the impacts of sh-CCNB1 on the tumorigenesis of BLCA in vivo. Analysis of sequencing data showed that CCNB1 mRNA was significantly elevated in tumor and BLCA compared with normal tissues [standardized mean difference (SMD)=1.21; 95% CI: 0.26-2.15; I²=95.9%]. IHC indicated that CCNB1 protein was localized in the nucleus and cytoplasm and was significantly increased in BLCA tumor tissues. The in vitro tests demonstrated that proliferation of T24 and 5637 cells transfected with sh-CCNB1 was significantly inhibited and cell migration and invasion ability were significantly decreased. sh-CCNB1 decreased the percentage of T24 cells in G0/G1, 5637 cells in the G0/G1 phase and S phase and increased percentage of 5637 cells in the G2/M phase and increased early apoptosis of 5637 cells. The in vivo experiments demonstrated that the mass of transplanted tumors was significantly decreased compared with the control group following silencing of CCNB1. The present results suggested that CCNB1 was involve in the development and prognosis of BLCA and silencing of CCNB1 may be a promising targeted therapy for BLCA.

Bioinformatics Identification of Therapeutic Gene Targets for Gastric Cancer

INTRODUCTION

The global prevalence of gastric cancer (GC) is increasing, and novel chemotherapeutic targets are needed.

METHODS

We searched for potential biomarkers for GC in three microarray data sets within the Gene Expression Omnibus (GEO) database. FunRich (v3.1.3) was used to perform Gene Ontology (GO) analyses and STRUN and Cytoscape (v3.6.0) were employed to construct a protein-protein interaction (PPI) network. To explore hub gene expression and survival, we used Gene Expression Profiling Interactive Analysis (GEPIA) and Kaplan-Meier (KM) plotter. Drugs that were closely related to key genes were screened by the Gene Set Cancer Analysis (GSCA), and relevant correlations were verified experimentally. We validated that the sensitivity of a GC cell line to these drugs was correlated with fibrillin 1 (FBN1) mRNA expression levels.

RESULTS

We identified 83 upregulated and 133 downregulated differentially expressed genes (DEGs) and these were enriched with regards to their cellular component (extracellular and exosomes), molecular function (extracellular matrix structural constituent and catalytic activity), and biological process (cell growth and/or maintenance and metabolism). The biological pathways most prominently involved were epithelial-to-mesenchymal transition (EMT) and β3 integrin cell surface interactions. For the PPI network, we selected 10 hub genes, and 70% of these were significantly connected to poor overall survival (OS) in patients with GC. We found a significant link between the expression of FBN1 and two small molecule drugs (PAC-1 and PHA-793887).

CONCLUSIONS

Overall, we suggest that these hub genes can be used as biomarkers and novel targets for GC. FBN1 may be associated with drug resistance in gastric cancer.

Weighted Gene Co-Expression Network Analysis of Oxymatrine in Psoriasis Treatment

Purpose

Psoriasis is a common, chronic, inflammatory, recurrent, immune-mediated skin disease. Oxymatrine is effective for treating moderate and severe psoriasis. Here, transcriptional changes in skin lesions before and after oxymatrine treatment of patients with psoriasis were identified using full-length transcriptome analysis and then compared with those of normal skin tissues.

Patients and Methods

Co-expression modules were constructed by combining the psoriasis area and severity index (PASI) score with weighted gene co-expression network analysis to explore the action mechanism of oxymatrine in improving clinical PASI. The expression of selected genes was verified using immunohistochemistry, quantitative real-time PCR, and Western blotting.

Results

Kyoto Encyclopedia of Gene and Genome pathway analysis revealed that oxymatrine treatment reversed the abnormal pathways, with an improvement in lesions and a reduction in PASI scores. Gene Ontology (GO) analysis revealed that oxymatrine treatment led to altered GO terms being regulated with a decrease in the PASI score in patients. Therefore, oxymatrine treatment may improve the skin barrier, differentiation of keratinocytes, and alleviate abnormality of organelles such as desmosomes. Protein-protein interaction network interaction analysis revealed that the top five hub genes among many interrelated genes were CNFN, S100A8, SPRR2A, SPRR2D, and SPRR2E, associated with the epidermal differentiation complex (EDC). EDC regulates keratinocyte differentiation. This result indicates that oxymatrine treatment can restore keratinocyte differentiation by regulating the expression of EDC-related genes.

Conclusion

Oxymatrine can improve erythema, scales, and other clinical symptoms of patients with psoriasis by regulating EDC-related genes and multiple pathways, thereby promoting the repair of epithelial tissue and maintaining the dynamic balance of skin keratosis.

A novel pyroptosis-related prognostic lncRNAs signature, tumor immune microenvironment and the associated regulation axes in bladder cancer

Bladder cancer (BC) is the most common malignancy of the urinary system. Pyroptosis is a host programmed cell death. However, the effects of pyroptosis-related lncRNAs (PRLs) on BC have not yet been completely elucidated. In this study, a prognostic PRLs model and two ceRNA networks were established using sufficient bioinformatics analysis and preliminary RT-qPCR validation in vitro. 6 PRLs were identified to construct a prognostic model. Then, the prognostic model risk score was verified to be an effective independent factor (Training cohort: Univariate analysis: HR = 1.786, 95% Cl = 1.416-2.252, p < 0.001; multivariate analysis: HR = 1.664, 95% Cl = 1.308-2.116, p < 0.001; testing cohort: Univariate analysis: HR = 1.268, 95% Cl = 1.144-1.405, p < 0.001; multivariate analysis: HR = 1.141, 95% Cl = 1.018-1.280, p = 0.024). Moreover, ROC and nomogram were performed to assess the accuracy of this signature (1-year-AUC = 0.764, 3-years-AUC = 0.769, 5-years-AUC = 0.738). Consequently, we evaluated the survival curves of these 6 lncRNAs using Kaplan–Meier survival analysis, demonstrating that MAFG-DT was risk lncRNA, while OCIAD1-AS1, SLC25A25-AS1, SNHG18, PSMB8-AS1 and TRM31-AS1 were protective lncRNAs. We found a strong correlation between PRLs and tumor immune microenvironment by Pearson’s correlation analysis. As for sensitivity of anti-tumor drugs, the high-risk group was more sensitive to Sorafenib, Bicalutamide and Cisplatin, while the low-risk group was more sensitive to AKT.inhibitor.VIII, Salubrinal and Lenalidomide, etc. Meanwhile, we identified lncRNA OCIAD1-AS1/miR-141-3p/GPM6B and lncRNA OCIAD1-AS1/miR-200a-3p/AKAP11 regulatory axes, which may play a potential role in the progression of BC.

Tumor suppressive role of microRNA-139-5p in bone marrow mesenchymal stem cells-derived extracellular vesicles in bladder cancer through regulation of the KIF3A/p21 axis

The emerging roles of extracellular vesicles (EVs) in bladder cancer have recently been identified. This study aims to elucidate the role of microRNA-139-5p (miR-139-5p) shuttled by bone marrow mesenchymal stem cells (BMSCs)-derived EVs (BMSCs-EVs) in bladder cancer, with the possible mechanism explored. Expression of miR-139-5p and KIF3A was tested, followed by an analysis of their correlation. EVs were isolated from BMSCs and co-cultured with T24 or BOY-12E cells with miR-139-5p mimic/inhibitor, oe-KIF3A, and/or si-p21 transfected to study the roles of miR-139-5p/KIF3A/p21 in bladder cancer cell functions. A nude mouse model of subcutaneous xenograft tumor was constructed to detect the effect of miR-139-5p in BMSCs-EVs on the tumorigenesis and lung metastasis of bladder cancer cells in vivo. It was identified that miR-139-5p was highly expressed in BMSCs-EVs, but poorly expressed in bladder cancer. BMSCs-EVs transferred miR-139-5p into bladder cancer cells where miR-139-5p inhibited the malignant features of bladder cancer cells in vitro. miR-139-5p in BMSCs-EVs targeted KIF3A and inhibited the expression of KIF3A, thereby activating p21. miR-139-5p in BMSCs-EVs arrested the tumorigenesis and lung metastasis of bladder cancer cells in vivo by modulation of the KIF3A/p21 axis. Altogether, BMSCs-EVs carried miR-139-5p targeted KIF3A to activate p21, thus delaying the occurrence of bladder cancer.

An Integrated Bioinformatics Analysis towards the Identification of Diagnostic, Prognostic, and Predictive Key Biomarkers for Urinary Bladder Cancer

Simple Summary

Bladder cancer is evidently a challenge as far as its prognosis and treatment are concerned. The investigation of potential biomarkers and therapeutic targets is indispensable and still in progress. Most studies attempt to identify differential signatures between distinct molecular tumor subtypes. Therefore, keeping in mind the heterogeneity of urinary bladder tumors, we attempted to identify a consensus gene-related signature between the common expression profile of bladder cancer and control samples. In the quest for substantive features, we were able to identify key hub genes, whose signatures could hold diagnostic, prognostic, or therapeutic significance, but, primarily, could contribute to a better understanding of urinary bladder cancer biology.

Abstract

Bladder cancer (BCa) is one of the most prevalent cancers worldwide and accounts for high morbidity and mortality. This study intended to elucidate potential key biomarkers related to the occurrence, development, and prognosis of BCa through an integrated bioinformatics analysis. In this context, a systematic meta-analysis, integrating 18 microarray gene expression datasets from the GEO repository into a merged meta-dataset, identified 815 robust differentially expressed genes (DEGs). The key hub genes resulted from DEG-based protein–protein interaction and weighted gene co-expression network analyses were screened for their differential expression in urine and blood plasma samples of BCa patients. Subsequently, they were tested for their prognostic value, and a three-gene signature model, including COL3A1, FOXM1, and PLK4, was built. In addition, they were tested for their predictive value regarding muscle-invasive BCa patients’ response to neoadjuvant chemotherapy. A six-gene signature model, including ANXA5, CD44, NCAM1, SPP1, CDCA8, and KIF14, was developed. In conclusion, this study identified nine key biomarker genes, namely ANXA5, CDT1, COL3A1, SPP1, VEGFA, CDCA8, HJURP, TOP2A, and COL6A1, which were differentially expressed in urine or blood of BCa patients, held a prognostic or predictive value, and were immunohistochemically validated. These biomarkers may be of significance as prognostic and therapeutic targets for BCa.

Deoxyschizandrin Inhibits the Proliferation, Migration, and Invasion of Bladder Cancer Cells through ALOX5 Regulating PI3K-AKT Signaling Pathway

Objective

Deoxyschizandrin has a significant inhibitory effect on a variety of tumor cells. However, the effect of Deoxyschizandrin on bladder cancer cells and its mechanism are still unclear.

Methods

Bladder cancer cells were treated with different concentrations of Deoxyschizandrin for 24 h, 48 h, and 72 h. The inhibition rate of cell proliferation was detected by CCK-8 assay. The changes of cell migration and invasion were detected by wound healing and Transwell assay. Based on the structure of Deoxyschizandrin, the protein targets of Deoxyschizandrin were predicted by bioinformatics database and verified by RNA and protein. Then, the expressions of ALOX5 and PI3K-AKT signaling pathway proteins were detected by Western blot in bladder cancer cells treated with Deoxyschizandrin.

Result

Deoxyschizandrin inhibited the proliferation, migration, and invasion of bladder cancer cells in a time- and concentration-dependent manner. Bioinformatics analysis showed that Deoxyschizandrin had 100 protein targets; among them, the score of ALOX5 was the highest, and the mRNA and protein levels of ALOX5 decreased after treatment with different concentrations of Deoxyschizandrin. Western blot results showed that compared with the control group, Deoxyschizandrin could significantly reduce the expression of p-PI3K and p-AKT, and overexpression of ALOX5 could significantly enhance the expression of p-PI3K and p-AKT. Compared with Deoxyschizandrin or overexpression of ALOX5, the expression of p-PI3K and p-AKT of Deoxyschizandrin combined with overexpression of ALOX5 recovered.

Conclusion

Deoxyschizandrin inhibits the proliferation, migration, and invasion of bladder cancer cells through ALOX5 regulating PI3K-AKT signaling pathway.

Proteins in the pathway from high red blood cell width distribution to all-cause mortality

Background

The pathophysiological mechanisms underlying the association between red blood cell distribution width (RDW) and all-cause mortality are unknown. We conducted a data-driven discovery investigation to identify plasma proteins that mediate the association between RDW and time to death in community-dwelling adults.

Methods

At baseline, 962 adults (women, 54·4%; age range, 21–98 years) participated in the InCHIANTI, “Aging in the Chianti Area” study, and proteomics data were generated from their plasma specimens. Of these, 623 participants had proteomics data available at the 9-year follow-up. For each visit, a total of 1301 plasma proteins were measured using SOMAscan technology. Complete data on vital status were available up to the 15-year follow-up period. Protein-specific exponential distribution accelerated failure time, and linear regression analyses adjusted for possible covariates were used for mortality and mediation analyses, respectively (survival data analysis).

Findings

Baseline values of EGFR, GHR, NTRK3, SOD2, KLRF1, THBS2, TIMP1, IGFBP2, C9, APOB, and LRP1B mediated the association between baseline RDW and all-cause mortality. Changes in IGFBP2 and C7 over 9 years mediated the association between changes in RDW and 6-year all-cause mortality.

Interpretation

Cellular senescence may contribute to the association between RDW and mortality.

Funding

This study was funded by grants from the National Institutes of Health (NIH) and the National Institute on Aging (NIA) contract and was supported by the Intramural Research Program of the NIA, NIH. The InCHIANTI study was supported as a ‘targeted project’ by the Italian Ministry of Health and in part by the U.S. NIA.

Identification of glycolysis genes signature for predicting prognosis in malignant pleural mesothelioma by bioinformatics and machine learning

BACKGROUND

Glycolysis-related genes as prognostic markers in malignant pleural mesothelioma (MPM) is still unclear. We hope to explore the relationship between glycolytic pathway genes and MPM prognosis by constructing prognostic risk models through bioinformatics and machine learning.

METHODS

The authors screened the dataset GSE51024 from the GEO database for Gene set enrichment analysis (GSEA), and performed differentially expressed genes (DEGs) of glycolytic pathway gene sets. Then, Cox regression analysis was used to identify prognosis-associated glycolytic genes and establish a risk model. Further, the validity of the risk model was evaluated using the dataset GSE67487 in GEO database, and finally, a specimen classification model was constructed by support vector machine (SVM) and random forest (RF) to further screen prognostic genes.

RESULTS

By DEGs, five glycolysis-related pathway gene sets (17 glycolytic genes) were identified to be highly expressed in MPM tumor tissues. Also 11 genes associated with MPM prognosis were identified in TCGA-MPM patients, and 6 (COL5A1, ALDH2, KIF20A, ADH1B, SDC1, VCAN) of them were included by Multi-factor COX analysis to construct a prognostic risk model for MPM patients, with Area under the ROC curve (AUC) was 0.830. Further, dataset GSE67487 also confirmed the validity of the risk model, with a significant difference in overall survival (OS) between the low-risk and high-risk groups (P < 0.05). The final machine learning screened the five prognostic genes with the highest risk of MPM, in order of importance, were ALDH2, KIF20A, COL5A1, ADH1B and SDC1.

CONCLUSIONS

A risk model based on six glycolytic genes (ALDH2, KIF20A, COL5A1, ADH1B, SDC1, VCAN) can effectively predict the prognosis of MPM patients.

Orthodenticle homeobox OTX1 is a potential prognostic biomarker for bladder cancer

Bladder cancer (BC) is one of the most aggressive tumors worldwide. OTX1 (orthodenticle homeobox 1) is an important transcription factor involved in various diseases, such as cancers. The aim of this study was to further investigate the role of OTX1 in BC. In this study, differentially expressed genes (DEGs) were screened from tumor tissues and para-cancerous tissues by bioinformatics. The expression of protein and RNA was separately detected by western blotting and immunohistochemistry (IHC), and quantitative polymerase chain reaction (qPCR); cell viability and cell growth were determined by 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) and clone formation assays, respectively; cell motility was measured by transwell and wound healing assays; cell cycle was measured by flow cytometry. In this study, 9 DEGs were screened out, and OTX1 was employed as a candidate gene for subsequent study. Results found that OTX1 was highly expressed in BC cells and BC tissues, which was significantly associated with poor prognosis of patients. In addition, OTX1 silencing significantly reduced cell viability, and inhibited cell growth and motility, while OTX1 overexpression got opposite results. Moreover, OTX1 co-expressed genes were enriched in cell cycle-related pathways, suggesting that the role of OTX1 in BC may be related to cell cycle, which was confirmed by flow cytometry analysis. Furthermore, in vivo experiments showed that OTX1 silencing significantly inhibited tumor growth in tumor-bearing mice. Taken together, our findings suggested that OTX1 may play a promotional role in BC progression.

A pancancer overview of FBN1, asprosin and its cognate receptor OR4M1 with detailed expression profiling in ovarian cancer

Ovarian cancer affects >295,000 women worldwide and is the most lethal of gynaecological malignancies. Often diagnosed at a late stage, current research efforts seek to further the molecular understanding of its aetiopathogenesis and the development of novel biomarkers. The present study investigated the expression levels of the glucogenic hormone asprosin [encoded by fibrillin-1 (FBN1)], and its cognate receptor, olfactory receptor 4M1 (OR4M1), in ovarian cancer. A blend of in silico open access The Cancer Genome Atlas data, as well as in vitro reverse transcription-quantitative PCR (RT-qPCR), immunohistochemistry and immunofluorescence data were used. RT-qPCR revealed expression levels of OR4M1 and FBN1 in clinical samples and in ovarian cancer cell lines (SKOV-3, PEO1, PEO4 and MDAH-2774), as well as the normal human ovarian surface epithelial cell line (HOSEpiC). Immunohistochemical staining of a tissue microarray was used to identify the expression levels of OR4M1 and asprosin in ovarian cancer samples of varying histological subtype and grade, including clear cell carcinoma, serous ovarian cancer and mucinous adenocarcinoma. Immunofluorescence analysis revealed asprosin expression in SKOV-3 and HOSEpiC cells. These results demonstrated the expression of both asprosin and OR4M1 in normal and malignant human ovarian tissues. This research invokes further investigation to advance the understanding of the role of asprosin and OR4M1 within the ovarian tumour microenvironment.

Prognostic signatures associated with high infiltration of Tregs in bone metastatic prostate cancer

Metastatic cancer especially bone metastasis (BM) is the lethal end-stage of castration-resistant prostate cancer (CRPC). To understand the possible molecular mechanisms underlying the development of the distant metastasis is of potential clinical value. We sought to identify differentially expressed genes between patient-matched primary and bone metastatic CRPC tumors. Functional enrichment, protein-protein interaction networks, and survival analysis of DEGs were performed. DEGs with a prognostic value considered as candidate genes were evaluated, followed by genetic analysis of tumor infiltrating immune cells based on Wilcoxon test and immunofluorescence identification. Expression profiles analysis showed that 381 overlapping genes were screened as differentially expressed genes (DEGs), of which 16 DEGs were randomly selected to be validated and revealed that most of these genes showed a transcriptional profile similar to that seen in the datasets (Pearson’s r = 0.76). Six core genes were found to be involved in regulation of extracellular matrix receptor interaction and chemotactic activity, and four of them were significantly correlated with the survival of PCa patients with bone metastases. Immune infiltration analysis showed that the expressions levels of COL3A1, RAC1, FN1, and SDC2 in CD4+T cells were significantly higher than those in tumor cells, especially regulatory T cell infiltration was significantly increased in BM tumors. We analyzed gene expression signatures specifically associated with the development of bone metastases of CRPC patients. Characterization of genes associated with BM of mCRPC is critical for identification of predictive biomarkers and potential therapeutic targets.

LINC-PINT Inhibited Malignant Progression of Bladder Cancer by Targeting miR-155-5p

BACKGROUND

This study mainly explored the expression level of LINC-PINT in bladder cancer and its relationship with prognosis. Meanwhile, the effect of LINC-PINT on the biological function of bladder cancer was also explored.

METHODS

The expression levels of LINC-PINT and miR-155-5p were detected by qRT-PCR. The prognostic significance of LINC-PINT in bladder cancer was studied by the Kaplan-Meier curve and Log rank test. CCK-8 and Transwell assays were used to analyze the proliferation, migration, and invasion ability. The targeting relationship between LINC-PINT and miR-155-5p was analyzed using bioinformatics and dual-luciferase reporter assays.

RESULTS

The expression of LINC-PINT was downregulated in bladder cancer tissues and cell lines, and miR-155-5p showed the opposite trend in bladder cancer tissues. Kaplan-Meier curve proved that the patients with low LINC-PINT expression had a lower five-year survival rate and the Log rank test displayed that LINC-PINT was a prognostic factor of BC. CCK-8 and Transwell results showed that LINC-PINT could inhibit the ability of proliferation, migration, and invasion. LINC-PINT was proved to target miR-155-5p in bladder cancer. Dual-luciferase reporter gene assay showed that the relative luciferase activity of overexpression miR-155-5p co-transfected with LINC-PINT-wt was significantly lower. LINC-PINT was negatively correlated with miR-155-5p.

CONCLUSION

LINC-PINT is a potential prognostic marker of bladder cancer, and the up-regulation of Lin-PINT can inhibit the proliferation, invasion, and migration of bladder cancer cells by targeting miR-155-5p.

An Integrated Data Analysis of mRNA, miRNA and Signaling Pathways in Pancreatic Cancer

Although many genes and miRNAs have been reported for various cancers, pancreatic cancer's specific genes or miRNAs have not been studied precisely yet. Therefore, we have analyzed the gene and miRNA expression profile of pancreatic cancer data in the gene expression omnibus (GEO) database. The microarray-derived miRNAs and mRNAs were annotated by gene ontology (GO) and signaling pathway analysis. We also recognized mRNAs that were targeted by miRNA through the mirDIP database. An integrated analysis of the microarray revealed that only 6 out of 43 common miRNAs had significant differences in their expression profiles between the tumor and normal groups (P value < 0.05 and |log Fold Changes (logFC)|> 1). The hsa-miR-210 had upregulation, whereas hsa-miR-375, hsa-miR-216a, hsa-miR-217, hsa-miR-216b and hsa-miR-634 had downregulation in pancreatic cancer (PC). The analysis results also revealed 109 common mRNAs by microarray and mirDIP 4.1 databases. Pathway analysis showed that amoebiasis, axon guidance, PI3K-Akt signaling pathway, absorption and focal adhesion, adherens junction, platelet activation, protein digestion, human papillomavirus infection, extracellular matrix (ECM) receptor interaction, and riboflavin metabolism played important roles in pancreatic cancer. GO analysis revealed the significant enrichment in the three terms of biological process, cellular component, and molecular function, which were identified as the most important processes associated strongly with pancreatic cancer. In conclusion, DTL, CDH11, COL5A1, ITGA2, KIF14, SMC4, VCAN, hsa-mir-210, hsa-mir-217, hsa-mir-216a, hsa-mir-216b, hsa-mir-375 and hsa-mir-634 can be reported as the novel diagnostic or even therapeutic markers for the future studies. Also, the hsa-mir-107 and hsa-mir-125a-5p with COL5A1, CDH11 and TGFBR1 genes can be introduced as major miRNA and genes on the miRNA-drug-mRNA network. The new regulatory network created in our study could give a deeper knowledge of the pancreatic cancer.

Role of COL3A1 and POSTN on Pathologic Stages of Esophageal Cancer

BACKGROUND

Esophageal cancer (EC) is a primary malignant tumor originating from the esophageal of the epithelium. Surgical resection is a potential treatment for EC, but this is only appropriate for patients who have locally resectable lesions suitable for surgery. However, most patients with EC are at a late stage when diagnosed. Therefore, there is an urgent need to further explore the pathogenesis of EC to enable early diagnosis and treatment.

METHODS

Our study downloaded 2 expression spectrum datasets (GSE92396 and GSE100942) in the Gene Expression Omnibus (GEO) database. GEO2 R was used to identify the Differentially expressed genes (DEGs) between the samples of EC and control. Using the DAVID tool to make the Functional enrichment analysis. Constructing A protein-protein interaction (PPI) network. Identifying the Hub genes. The impact of hub gene expression on overall survival and their expression based on immunohistochemistry were analyzed. Associated microRNAs were also predicted.

RESULTS

There were 36 common DEGs identified. The analysis of GO and KEGG results shown that the variations were predominantly concentrated in the extracellular matrix (ECM), ECM organization, DNA binding, platelet activation, and ECM-receptor interactions. COL3A1 and POSTN had high expression in EC tissues which was compared with their expression in healthy tissues. Analysis of pathologic stages showed that when COL3A1 and POSTN were highly expressed, the stage of the pathologic of EC patients was relatively high (P < 0.005).

CONCLUSIONS

COL3A1 and POSTN may play an important role in the advancement and occurrence of EC. These genes could provide some novel ideas and basis for the diagnosis and targeted treatment of EC.

Urinary glycoproteomic profiling of non-muscle invasive and muscle invasive bladder carcinoma patients reveals distinct N-glycosylation pattern of CD44, MGAM, and GINM1

Clinical management of bladder carcinomas (BC) remains a major challenge and demands comprehensive multi-omics analysis for better stratification of the disease. Identification of patients on risk requires identification of signatures predicting prognosis risk of the patients. Understanding the molecular alterations associated with the disease onset and progression could improve the routinely used diagnostic and therapy procedures. In this study, we investigated the aberrant changes in N-glycosylation pattern of proteins associated with tumorigenesis as well as disease progression in bladder cancer. We integrated and compared global N-glycoproteomic and proteomic profile of urine samples from bladder cancer patients at different clinicopathological stages (non-muscle invasive and muscle-invasive patients [n = 5 and 4 in each cohort]) with healthy subjects (n = 5) using SPEG method. We identified 635 N-glycopeptides corresponding to 381 proteins and 543 N-glycopeptides corresponding to 326 proteins in NMIBC and MIBC patients respectively. Moreover, we identified altered glycosylation in 41 NMIBC and 21 MIBC proteins without any significant change in protein abundance levels. In concordance with the previously published bladder cancer cell line N-glycoproteomic data, we also observed dysregulated glycosylation in ECM related proteins. Further, we identified distinct N-glycosylation pattern of CD44, MGAM, and GINM1 between NMIBC and MIBC patients, which may be associated with disease progression in bladder cancer. These aberrant protein glycosylation events would provide a novel approach for bladder carcinoma diagnosis and further define novel mechanisms of tumor initiation and progression.

Identification of prognostic biomarkers associated with stromal cell infiltration in muscle-invasive bladder cancer by bioinformatics analyses

Muscle-invasive bladder cancer (MIBC) is one of the common malignant tumors. Patients with MIBC still have high tumor recurrence and progression rates after surgery. Bioinformatics analysis of stromal infiltration-related genes in the tumor microenvironment (TME) of MIBC patients was performed in this study to determine the major stromal cells types and biomarkers for their progression and poor prognosis. The ESTIMATE algorithm was applied to evaluate the stromal score and immune score of samples from MIBC patients in The Cancer Genome Atlas (TCGA) and found that stromal score was closely related to the clinical characteristics of the patients. The Gene Set Enrichment Analysis (GSEA) further revealed that stromal cells were involved in biological processes such as activation of leukocytes and positive regulation of cell migration during MIBC progression, as well as PI3K-Akt, MAPK, and Rap1 signaling pathways. Five hub genes related to prognosis, including ACTA2, COL5A1, DCN, LUM, and PRRX1 were identified by the Weighted Gene Co-Expression Network Analysis (WGCNA), Protein-Protein Interaction (PPI), survival analysis, and Oncomine, Gene Expression Omnibus (GEO) database validation. Besides, we identified five stromal cell types associated with overall survival time, among which chondrocytes and fibroblasts were identified as the major stromal cell types through correlation analysis. Finally, the Receiver Operating Characteristic (ROC) curve and immunohistochemistry were used to verify the diagnostic value and expression of hub genes in different invasive tumors. In summary, we investigated the biological behavior of stromal cells in the TME of MIBC to promote tumor progression obtained hub genes associated with progression and poor prognosis and identified the main stromal cells types in the TME.

Reduced Expression of METTL3 Promotes Metastasis of Triple-Negative Breast Cancer by m6A Methylation-Mediated COL3A1 Up-Regulation

The abnormal m6A modification caused by m6A modulators is a common feature of various tumors; however, little is known about which m6A modulator plays the most important role in triple-negative breast cancer (TNBC). In this study, when analyzing the influence of m6A modulators (METTL3, METTL14, WTAP, FTO, and ALKBH5) on the prognosis of breast cancer, especially in TNBC using several on-line databases, methyltransferase-like 3 (METTL3) was found to have low expression in breast cancer, and was closely associated with short-distance-metastasis-free survival in TNBC. Further investigation showed that knockdown of METTL3 could enhance the ability of migration, invasion, and adhesion by decreasing m6A level in TNBC cell lines. Collagen type III alpha 1 chain (COL3A1) was identified and verified as a target gene of METTL3. METTL3 could down-regulate the expression of COL3A1 by increasing its m6A methylation, ultimately inhibiting the metastasis of TNBC cells. Finally, with immunohistochemistry staining in breast cancer tissues, it was proved that METTL3 expression was negatively correlated with COL3A1 in TNBC, but not in non-TNBC. This study demonstrated the potential mechanism of m6A modification in metastasis and provided potential targets for treatment in TNBC.

Weighted gene co-expression network analysis can sort cancer-associated fibroblast-specific markers promoting bladder cancer progression

The role of cancer-associated fibroblasts (CAFs) has been thoroughly investigated in tumour microenvironments but not in bladder urothelial carcinoma (BLCA). The cell fraction of CAFs gradually increased with BLCA progression. Weighted gene co-expression network analysis (WGCNA) revealed a specific gene expression module of CAFs that are relevant to cancer progression and survival status. Fifteen key genes of the module were consistent with a fibroblast signature in single-cell RNA sequencing, functionally related to the extracellular matrix, and significant in survival analysis and tumour staging. A comparison of the luminal-infiltrated versus luminal-papillary subtypes and fibroblast versus urothelial carcinoma cell lines and immunohistochemical data analysis demonstrated that the key genes were specifically expressed in CAFs. Moreover, these genes are highly correlated with previously reported CAF markers. In summary, CAFs play a major role in the progression of BLCA, and the 15 key genes act as BLCA-specific CAF markers and can predict CAF changes. WGCNA can, therefore, be used to sort CAF-specific gene set in cancer tissues.

Functional and genomic characterization of three novel cell lines derived from a metastatic gallbladder cancer tumor

Background

Gallbladder cancer (GBC) is the most common tumor of the biliary tract. The incidence of GBC shows a large geographic variability, being particularly frequent in Native American populations. In Chile, GBC represents the second cause of cancer-related death among women. We describe here the establishment of three novel cell lines derived from the ascitic fluid of a Chilean GBC patient, who presented 46% European, 36% Mapuche, 12% Aymara and 6% African ancestry.

Results

After immunocytochemical staining of the primary cell culture, we isolated and comprehensively characterized three independent clones (PUC-GBC1, PUC-GBC2 and PUC-GBC3) by short tandem repeat DNA profiling and RNA sequencing as well as karyotype, doubling time, chemosensitivity, in vitro migration capability and in vivo tumorigenicity assay. Primary culture cells showed high expression of CK7, CK19, CA 19-9, MUC1 and MUC16, and negative expression of mesothelial markers. The three isolated clones displayed an epithelial phenotype and an abnormal structure and number of chromosomes. RNA sequencing confirmed the increased expression of cytokeratin and mucin genes, and also of TP53 and ERBB2 with some differences among the three cells lines, and revealed a novel exonic mutation in NF1. The PUC-GBC3 clone was the most aggressive according to histopathological features and the tumorigenic capacity in NSG mice.

Conclusions

The first cell lines established from a Chilean GBC patient represent a new model for studying GBC in patients of Native American descent.

Plasma proteomic signature of the risk of developing mobility disability: A 9-year follow-up

INTRODUCTION

Mobility disability is a powerful indicator of poor health in older adults. The biological and pathophysiological mechanism underlying the development of mobility disability remains unknown. This study conducted a data-driven discovery phase investigation to identify plasma proteins that predict the incidence of mobility disability in community-dwelling older adults without mobility disability at baseline.

METHODS

We investigated 660 women and men, aged 71.9 ± 6.0 (60-94) years, who participated in the Invecchiare in Chianti, "Aging in the Chianti Area" study and completed the 400-m walk at fast pace (400-m walk) at enrollment. Median follow-up time was 8.57 [interquartile, 3.20-9.08] years. SOMAscan technology was used to measure 1,301 plasma proteins at enrollment. The incident of mobility disability was defined as inability to complete the 400-m walk. Protein-specific Cox proportional hazard model was adjusted for sex, age, and other important covariates.

RESULTS

Plasma levels of 75 proteins predicted mobility disability (p < .05). Significant proteins were enriched for the KEGG "PI3K-Akt signaling," "phagosomes," and "cytokine-cytokine receptor interaction" pathways. After multiple comparison adjustment, plasma cathepsin S (CTSS; hazard ratio [HR] 1.33, 95% CI: 1.17, 1.51, q = 0.007), growth/differentiation factor 15 (GDF15; HR: 1.45, 95% CI: 1.23, 1.72, q = 0.007), and thrombospondin-2 (THBS2; HR: 1.44, 95% CI: 1.22, 1.69, q = 0.007) remained significantly associated with high risk of losing mobility.

CONCLUSION

CTSS, GDF15, and THBS2 are novel blood biomarkers associated with new mobility disability in community-dwelling individuals. Overall, our analysis suggests that cellular senescence and inflammation should be targeted for prevention of mobility disability.

In silico analysis of the immune microenvironment in bladder cancer

BACKGROUND

Infiltrating immune and stromal cells are vital components of the bladder cancer (BC) microenvironment, which can significantly affect BC progression and outcome. However, the contribution of each subset of tumour-infiltrating immune cells is unclear. The objective of this study was to perform cell phenotyping and transcriptional profiling of the tumour immune microenvironment and analyse the association of distinct cell subsets and genes with BC prognosis.

METHODS

Clinical data of 412 patients with BC and 433 transcription files for normal and cancer tissues were downloaded from The Cancer Genome Atlas. The CIBERSORT algorithm was used to determine the relative abundance of 22 immune cell types in each sample and the ESTIMATE algorithm was used to identify differentially expressed genes within the tumour microenvironment of BC, which were subjected to functional enrichment and protein-protein interaction (PPI) analyses. The association of cell subsets and differentially expressed genes with patient survival and clinical parameters was examined by Cox regression analysis and the Kaplan-Meier method.

RESULTS

Resting natural killer cells and activated memory CD4+ and CD8+ T cells were associated with favourable patient outcome, whereas resting memory CD4+ T cells were associated with poor outcome. Differential expression analysis revealed 1334 genes influencing both immune and stromal cell scores; of them, 97 were predictive of overall survival in patients with BC. Among the top 10 statistically significant hub genes in the PPI network, CXCL12, FN1, LCK, and CXCR4 were found to be associated with BC prognosis.

CONCLUSION

Tumour-infiltrating immune cells and cancer microenvironment-related genes can affect the outcomes of patients and are likely to be important determinants of both prognosis and response to immunotherapy in BC.

Construction and validation of an eight-gene signature with great prognostic value in bladder cancer

Bladder cancer (BC) is one of the most common malignancies in urinary system with a common malignancy in urinary system with a high mortality and recurrence rate, so we attempt to construct a gene signature to predict the prognosis of BCs. We initially established a co-expression network by performing WGCNA analysis and further identified magenta module as key module (P = 8e-05, R2 = 0.4). Subsequently, we screened 12 genes associated with survival from the key module, which were selected to construct an eight-gene signature by establishing a LASSO Cox model. Moreover, we reckoned the risk score (RS) of each sample, through which we could divide samples into two groups (the high-risk and low-risk groups) and verify the signature, in the training set and 3 validation sets (internal test set, GSE13507and E-MTAB-4321). This signature could distinguish between the high- and low- risk patients well (survival analysis: P = 0.015; AUC: 0.61 at 1 year, 0.61 at 3 years and 0.61 at 5 years). In the validation sets, this signature also showed good performance, which was consistent with the training test. Furthermore, we plotted a nomogram to predict the possibility of the overall survival (OS) and three calibration curves to predict the effectiveness of the nomogram, which suggested good value and clinical utility of the nomogram. In conclusion, we established an eight-gene signature, which was probably effective in the prediction of prognosis of patients with BC.

Characterization of long non-coding RNA and messenger RNA profiles in laryngeal cancer by weighted gene co-expression network analysis

Laryngeal cancer (LC) is a malignant tumor in the head and neck region. It was recently elucidated that long non-coding RNAs (lncRNAs) participate in the pathogenesis of LC. However, the detailed mechanism of lncRNA in LC and whether long non-coding RNAs serve as effective biomarkers remains unclear. Ribonucleic acid (RNA) sequence data of LC and 11 patient clinical traits were extracted from The Cancer Genome Atlas (TCGA) database and analyzed by weighted gene co-expression network analysis (WGCNA). A total of 9 co-expression modules were identified. The co-expression Pink module significantly correlated with four clinical traits, including history of smoking, lymph node count, tumor status, and the success of follow-up treatment. Based on the co-expression Pink module, lncRNA-microRNA (miRNA)-messenger RNA (mRNA) and lncRNA-RNA binding protein-mRNA networks were constructed. We found that 8 lncRNAs significantly impacted overall survival (OS) in LC patients. These identified lncRNA and hub gene biomarkers were also validated in multiple LC cells in vitro via qPCR. Taken together, this study provided the framework of co-expression gene modules of LC and identified some important biomarkers in LC development and disease progression.