CSF2 Overexpression Is Associated with STAT5 Phosphorylation and Poor Prognosis in Patients with Urothelial Carcinoma

The potential mechanism of WT1-associated protein-induced N-6-methyladenosine modification of colony-stimulating factor 2 in the progression of oral squamous cell carcinoma by JAK/STAT3 pathway regulation

Colony-stimulating factor 2 (CSF2) plays a regulatory role in numerous cancers. However, there is needed to investigate the role of CSF2 in oral squamous cell carcinoma (OSCC) malignant phenotype and the specific mechanisms of CSF2 N-6-methyladenosine (m6A) modification. Therefore, we investigated the regulatory mechanism of m6A-modified CSF2 by WT1-associated protein (WTAP) in OSCC via qRT-PCR, western blot, WTAP and CSF2 overexpression in OSCC. In a panel of OSCCs, Kaplan-Meier plot analysis indicated that high expression of CSF2 was associated with poorer prognosis. Cell functional experiments revealed that enrichment of CSF2 promoted the proliferation and migration of OSCC cells by activating the JAK/STAT3 pathway, whereas the reduced expression of CSF2 resulted in the malignant decline of OSCC cells by blocking the JAK/STAT3 pathway. This study also confirmed that WTAP enhanced the m6A level of CSF2 and facilitated the expression of CSF2 and that CSF2 silencing blocked the invasive phenotype of OSCC cells and reversed the malignancy induced by WTAP overexpression. Overall, this study demonstrated that WTAP mediates the m6A modification of CSF2 and the JAK/STAT3 pathway, which plays an oncogenic role in the development of OSCC and can be a target for the treatment of patients with OSCC.

CSF2 Impairs Nrf2 Signaling through the Akt/Mtor Pathway in the Development of Bladder Cancer

Bladder Cancer (BCa) is one of the most common cancers of the urinary system. Colony-stimulating factor 2 (CSF2) is involved in many cancers, but not BCa. We investigated the effect of CSF2 on BCa in this study and the underlying molecular mechanisms. CSF2 mRNA levels in BCa were analyzed using the Cancer Genome Atlas (TCGA) database. Western blot was conducted to verify CSF2 expression in BCa tissue samples and cell lines. The effect of CSF2 on the growth of BCa cells was assessed by CCK8 and colony formation. To determine the migration and invasion capabilities of BCa cells, transwell analysis and wound healing assays were conducted. Next, western blot was used to explore the underlying mechanism. In the end, a xenografted BCa mouse model was established to examine the effects of CSF2 on tumorigenesis in vivo. Results showed that CSF2 mRNA was upregulated in BCa samples. Knocking down CSF2 significantly inhibited the proliferation and tumorigenesis of BCa cells in vitro and in vivo. Mechanism analysis revealed that CSF2 knockdown inhibited the proliferation and invasion of BCa cells via AKT/mTOR signaling. Based on these results, CSF2 promotes the proliferation and tumorigenesis of BCa.

Analysis of Zinc and Stromal Immunity in Disuse Osteoporosis: Mendelian Randomization and Transcriptomic Analysis

OBJECTIVE

Disuse osteoporosis is known to be primarily caused by a lack of exercise. However, the causal relationships between zinc and immunity and disuse osteoporosis remain unknown. This study investigated these relationships and their potential mechanisms.

METHODS

This study was an integrative study combining genome-wide association studies and transcriptomics. Two-sample Mendelian randomization analysis (MR) was used to analyze the causal relationships between exposures (zinc, immunity, physical activity) and the outcome (osteoporosis) with the aid of single-nucleotide polymorphisms (SNPs) as instrumental variables (IVs). Four models, MR-Egger, inverse variance weighted, weighted median and MR-Pleiotrophy RESidual Sum and Outlier (MRPRESSO), were used to calculate odds ratio values. Sensitivity and heterogeneity analyses were also performed using MRPRESSO and MR-Egger methods. The mRNA transcriptomic analysis was subsequently conducted. Zinc metabolism scores were acquired through single-sample Gene Set Enrichment Analysis algorithms. Stromal scores were obtained using the R Package "estimate" algorithms. Important Kyoto Encyclopedia of Genes and Genomes and Gene Ontology pathways were also derived through gene set variation analysis. Cytoscape software helped construct the transcription factor (TF)-mRNA-microRNA (miRNA) network. Virtual screening and molecular docking were performed. Polymerase chain reaction validation was also carried out in vivo.

RESULTS

Causal relationships were demonstrated between zinc and exercise (95% confidence interval [CI] = 1.30-2.95, p = 0.001), exercise and immunity (95% CI = 0.36-0.80, p = 0.002), exercise and osteoporosis (95% CI = 0.97-0.99, p = 0.0007), and immunity disorder and osteoporosis (95% CI = 1.30-2.03, p = 0.00002). One hundred and seventy-nine mRNAs in important modules were screened. Combining the differential expressional genes (DEGs) and the Boruta selection, six DEGs were screened (AHNAK, CSF2, ADAMTS12, SRA1, RUNX2, and SLC39A14). TF HOXC10 and miRNA hsa-miR-204 were predicted. Then, the TF-mRNA-miRNA network was successfully constructed. RUNX2 and SLC39A14 were identified as hub mRNAs in the TF-mRNA-miRNA network. Eventually, the novel small drug C6O4NH5 was designed according to the pharmacophore structure of SLC39A14. The docking energy for the novel drug was -5.83 kcal/mol. SLC39A14 and RUNX2 were downregulated (of statistical significance p-value < 0.05) in our animal experiment.

CONCLUSION

This study revealed that zinc had a protective causal relationship with disuse osteoporosis by promoting exercise and immunity. SLC39A14 and RUNX2 mRNA participated in this zinc-related mechanism.

Unravelling the Role of Iroquois Homeobox 4 and its Interplay with Androgen Receptor in Prostate Cancer

Genome-wide association studies have linked Iroquois-Homeobox 4 (IRX4) as a robust expression quantitative-trait locus associated with prostate cancer (PCa) risk. However, the intricate mechanism and regulatory factors governing IRX4 expression in PCa remain poorly understood. Here, we unveil enrichment of androgen-responsive gene signatures in metastatic prostate tumors exhibiting heightened IRX4 expression. Furthermore, we uncover a novel interaction between IRX4 and the androgen receptor (AR) co-factor, FOXA1, suggesting that IRX4 modulates PCa cell behavior through AR cistrome alteration. Remarkably, we identified a distinctive short insertion-deletion polymorphism (INDEL), upstream of the IRX4 gene that differentially regulates IRX4 expression through the disruption of AR binding. This INDEL emerges as the most significant PCa risk-associated variant within the 5p15 locus, in a genetic analysis involving 82,591 PCa cases and 61,213 controls and was associated with PCa survival in patients undergoing androgen-deprivation therapy. These studies suggest the potential of this INDEL as a prognostic biomarker for androgen therapy in PCa and IRX4 as a potential therapeutic target in combination with current clinical management.

Establishment of two oxaliplatin-resistant gallbladder cancer cell lines and comprehensive analysis of dysregulated genes

Acquired resistance to chemotherapeutic drugs in gallbladder cancer (GBC) results in therapy failure. This study is aimed to establish oxaliplatin (OXA)-resistant GBC cell lines and uncover their gene expression profiles. First, two OXA-resistant GBC cell lines (GBC-SD/OXA and SGC996/OXA) were established by gradually increasing the drug concentration, and the resistance index was 4-5. The two resistant cell lines showed slower proliferation and higher stemness, colony formation, and migration abilities. Epithelial mesenchymal transformation and increased levels of P-glycoprotein were also detected. Next RNA-sequence analysis identified 4,675 dysregulated genes (DGs) in resistant cells, and most of the 12 randomly selected DGs were verified to be consistent with the sequence results. Kyoto Encyclopedia of Genes and Genomes analysis indicated that several DGs were involved in resistance- and phenotype-related pathways, of which the activations of PD-L1 and ERK1/2 were both verified in resistant cell lines. In conclusion, this study is the first to report the gene expression profile of OXA-resistant GBC cells and provides a useful database for target development.

Lactoferrin attenuated urban particulate matter-induced nephrotoxicity by regulating the CSF2/CENPE axis

Several epidemiological studies regarding the adverse effect of air pollution have notably accelerated in recent years. Urban particulate matter (PM) gains access to the respiratory system and translocates into the circulation to affect several tissues, such as the liver and kidneys. Lactoferrin is a substance belonging to the non-heme iron-binding glycoprotein which is present in breast milk and other exocrine fluids. Lactoferrin is protective against many pathophysiological conditions. In the present study, we explored the potential influence of lactoferrin on PM-induced nephrotoxicity. We found that lactoferrin rescued PM-induced cell death but did not affect apoptosis in human kidney cells. Lactoferrin decreased necroptosis and fibrosis but increased autophagy in human kidney cells. Furthermore, the gene expression profiles of PM and lactoferrin were analyzed by RNA sequencing. The transcriptional profiles were uploaded and analyzed by ingenuity pathway analysis software and gene set enrichment analysis. The results showed that the crucial role of the CSF2/CENPE pathway was involved in human kidney cells treated with PM and lactoferrin. In a mouse model, lactoferrin ameliorates PM-induced nephrotoxicity by regulating necroptosis, fibrosis, autophagy and the CSF2/CENPE axis. In summary, these findings showed that lactoferrin could be a novel therapeutic or preventive agent for renal disorders caused by airborne PM pollution.

A Novel Lymphangiogenesis-Related Gene Signature can Predict Prognosis and Immunosuppressive Microenvironment in Patients with Clear Cell Renal Cell Carcinoma

Background: Lymphangiogenesis represents a key event in the progression and metastasis of patients with clear cell renal cell carcinoma (ccRCC). Nevertheless, the prognostic value of lymphangiogenesis-related genes (LRGs) in ccRCC patients remains unknown. Method: Differential analyses were performed to identify differentially expressed LRGs between normal and tumor tissues. A univariate Cox analysis was performed to identify differently expressed LRGs associated with overall survival (OS). LASSO and multivariate Cox analyses were performed to construct and optimize the LRG signature. To further explore the molecular characterization of the LRG signature, a functional enrichment analysis, immune signature, somatic mutations, and drug sensitivity were assessed. Immunohistochemistry (IHC) and immunofluorescence staining were performed to validate the relationship between lymphangiogenesis and immunity using our ccRCC samples. Results: Four candidate genes (IL4, CSF2, PROX1, and TEK) were eventually available to construct the LRG signature in the training set. Patients in the high-risk group had a shorter survival than those in the low-risk group. The LRG signature was an independent prognostic factor of OS. These results were confirmed in the validation group. The LRG signature was correlated with immunosuppressive cell infiltration, T cell exhaustion markers, somatic mutations, and drug sensitivity. The IHC and immunofluorescence staining results confirmed the correlation between lymphangiogenesis and CD163⁺ macrophages, exhausted CD8⁺PD-1⁺, and CD8⁺ LAG3⁺ T cells. Conclusion: A novel prognostic signature based on LRGs could provide insight into the prognostic evaluation and treatment of ccRCC patients.

Identifying and characterizing drug sensitivity-related lncRNA-TF-gene regulatory triplets

Recently, many studies have shown that lncRNA can mediate the regulation of TF-gene in drug sensitivity. However, there is still a lack of systematic identification of lncRNA-TF-gene regulatory triplets for drug sensitivity. In this study, we propose a novel analytic approach to systematically identify the lncRNA-TF-gene regulatory triplets related to the drug sensitivity by integrating transcriptome data and drug sensitivity data. Totally, 1570 drug sensitivity-related lncRNA-TF-gene triplets were identified, and 16 307 relationships were formed between drugs and triplets. Then, a comprehensive characterization was performed. Drug sensitivity-related triplets affect a variety of biological functions including drug response-related pathways. Phenotypic similarity analysis showed that the drugs with many shared triplets had high similarity in their two-dimensional structures and indications. In addition, Network analysis revealed the diverse regulation mechanism of lncRNAs in different drugs. Also, survival analysis indicated that lncRNA-TF-gene triplets related to the drug sensitivity could be candidate prognostic biomarkers for clinical applications. Next, using the random walk algorithm, the results of which we screen therapeutic drugs for patients across three cancer types showed high accuracy in the drug-cell line heterogeneity network based on the identified triplets. Besides, we developed a user-friendly web interface-DrugSETs (http://bio-bigdata.hrbmu.edu.cn/DrugSETs/) available to explore 1570 lncRNA-TF-gene triplets relevant with 282 drugs. It can also submit a patient’s expression profile to predict therapeutic drugs conveniently. In summary, our research may promote the study of lncRNAs in the drug resistance mechanism and improve the effectiveness of treatment.

Perineural Invasion is a Powerful Prognostic Factor for Upper Tract Urothelial Carcinoma Following Radical Nephroureterectomy

BACKGROUND

Taiwan has the highest incidence of upper tract urothelial carcinoma (UTUC) worldwide. Although many pathological factors can predict the prognosis of UTUC, previous studies have rarely discussed perineural invasion (PNI). Therefore, we aimed to investigate the effect of PNI on a well-established cohort of patients with UTUC.

METHODS

This retrospective study included 803 patients with non-metastatic UTUC who underwent radical nephroureterectomy between June 2000 and August 2019. Demographic and clinicopathological parameters, including PNI, were collected for analysis. Using the Kaplan-Meier method and Cox proportional hazards model, we evaluated the significance of PNI with respect to progression-free survival (PFS), cancer-specific survival (CSS), and overall survival (OS).

RESULTS

The median follow-up was 30.9 months, and there were 83 cases of PNI (10.3%). PNI-positive patients had unfavorable pathological features, including high pT stage, positive lymph node involvement, high tumor grade, and more lymphovascular invasion (all p < 0.001). Kaplan-Meier analysis showed that PNI was significantly associated with PFS, CSS, and OS (all p < 0.00001), and when combined with lymphovascular invasion, patients could be divided into groups with distinct survival rates (all p < 0.00001). In multivariate analysis, PNI was an independent factor leading to worse PFS (hazard ratio [HR] 1.72, 95% confidence interval [CI] 1.19-2.50; p = 0.004), CSS (HR 2.54, 95% CI 1.58-4.10; p = 0.0001), and OS (HR 1.78, 95% CI 1.19-2.65; p = 0.005).

CONCLUSIONS

We demonstrated an association between PNI and the prognosis of UTUC. Routine assessment of PNI in UTUC with standardized protocols may help achieve better risk stratification and subject selection for perioperative treatment.

Construction of a Two-Gene Immunogenomic-Related Prognostic Signature in Lung Squamous Cell Carcinoma

Lung cancer has the highest tumor incidence in China. Lung squamous cell carcinoma (LUSC) is the most common type, accounting for 40–51% of primary lung cancers. LUSC is slow in growth and late in metastasis. Immune-related genes (IRGs) and immune infiltrating cells play a vital role in the clinical outcomes of LUSC. It is important to systematically study its immune gene map to help the prognosis of cancer patients. In this study, we combined the prognostic landscape and expression status of IRGs downloaded from the TCGA and InnatedDB databases and systematically analyzed the prognostic information of LUSC patients to obtain IRGs. After systematically exploring the survival analysis, prognosis-related genes were found, and the PPI network revealed that a total of 11 genes were hub genes. A two-gene prognosis risk model was established by multivariate Cox analysis. Two IRGs were closely correlated with the prognosis of LUSC. Based on these two genes, a new independent prognostic risk model was established, and this model was further verified in the GEO database. Moreover, the risk score of the model was correlated with sex, survival status, and lymphatic metastasis in LUSC patients, and the predictive risk of the prognostic risk model was significantly positively correlated with five kinds of immune cells (CD4 T cells, CD8 T cells, neutrophils, macrophages, and dendritic cells). This study comprehensively analyzed immunogenomics and presented immune-related prognostic biomarkers for LUSC.

A Neutrophil Extracellular Traps Signature Predicts the Clinical Outcomes and Immunotherapy Response in Head and Neck Squamous Cell Carcinoma

Background: Neutrophil extracellular traps (NETs) play an important role in the occurrence, metastasis and immune escape of cancers. This study aimed to investigate NET-related genes, their clinical prognostic value and their correlation with immunotherapy and anticancer drugs in patients with head and neck squamous cell carcinoma (HNSCC).

Methods: Differentially expressed NET-related genes in HNSCC were identified based on multiple public databases. To improve the clinical practicability and avoid overfitting, univariable, least absolute shrinkage and selection operator (LASSO) and multivariable Cox algorithms were used to construct a prognostic risk model. A nomogram was further used to explore the clinical value of the model. Internal and external validation were conducted to test the model. Furthermore, the immune microenvironment, immunophenoscore (IPS) and sensitivity to anticancer drugs in HNSCC patients with different prognostic risks were explored.

Results: Six NET-related genes were screened to construct the risk model. In the training cohort, Kaplan–Meier (K-M) analysis showed that the overall survival (OS) of low-risk HNSCC patients was significantly better than that of high-risk HNSCC patients (p < 0.001). The nomogram also showed a promising prognostic value with a better C-index (0.726 vs 0.640) and area under the curve (AUC) (0.743 vs 0.706 at 3 years, 0.743 vs 0.645 at 5 years) than those in previous studies. Calibration plots and decision curve analysis (DCA) also showed the satisfactory predictive capacity of the nomogram. Internal and external validation further strengthened the credibility of the clinical prognostic model. The level of tumor mutational burden (TMB) in the high-risk group was significantly higher than that in the low-risk group (p = 0.017), and the TMB was positively correlated with the risk score (R = 0.11; p = 0.019). Moreover, the difference in immune infiltration was significant in HNSCC patients with different risks (p < 0.05). Furthermore, the IPS analysis indicated that anti-PD-1 (p < 0.001), anti-CTLA4 (p < 0.001) or combining immunotherapies (p < 0.001) were more beneficial for low-risk HNSCC patients. The response to anticancer drugs was also closely correlated with the expression of NET-related genes (p < 0.001).

Conclusion: This study identified a novel prognostic model that might be beneficial to develop personalized treatment for HNSCC patients.

Low Tumor Infiltrating Mast Cell Density Reveals Prognostic Benefit in Cervical Carcinoma

Objectives: Research on the role of mast cells (MCs) in cervical tumor immunity is more limited. Therefore, our study aimed to evaluate the prognostic value of MCs and their correlation with the immune microenvironment of cervical carcinoma (CC). Methods: The Cancer Genome Atlas (TCGA) data was utilized to obtain the degree of immune infiltration of MCs in CC. Meanwhile, this study retrospectively collected patient clinical characteristic data and tissue specimens to further verify the relevant conclusions. Mast cell density (MCD) was measured by the CIBERSORT algorithm in TCGA data and immunohistochemical staining of tryptase in CC tissues. Finally, differentially expressed genes (DEGs) of TCGA data were performed using "limma" packages and key gene modules were identified using the MCODE application in Cytoscape. Results: The results showed MCs were diffusely distributed in CC tissues. Moreover, we found that low tumor-infiltrating MCD was beneficial for overall survival (OS) in the TCGA cohort. Consistent conclusions were also obtained in a clinical cohort. In addition, a total of 305 DEGs were analyzed between the high tumor-infiltrating MCD and low tumor-infiltrating MCD group. Seven key modules, a total of 34 genes, were screened through the MCODE plug-in, which was mainly related to inflammatory response and immune response and closely correlated with cytokines including CSF2, CCL20, IL1A, IL1B, and CXCL8. Conclusion: In short, high tumor-infiltration MCs in CC tissue was associated with worse OS in patients. Furthermore, MCs were closely related to cytokines in the tumor microenvironment, suggesting that they collectively played a role in the immune response of the tumor. Therefore, MCD may be a potential prognostic indicator and immunotherapy target of CC.

HIST1H1B Promotes Basal-Like Breast Cancer Progression by Modulating CSF2 Expression

Background

Basal-like breast cancer (BLBC) is associated with a poor clinical outcome; however, the mechanism of BLBC aggressiveness is still unclear. It has been shown that a linker histone functions as either a positive or negative regulator of gene expression in tumors. Here, we aimed to investigate the possible involvement and mechanism of HIST1H1B in BLBC progression.

Experimental design

We analyzed multiple gene expression datasets to determine the relevance of HIST1H1B expression with BLBC. We employed quantitative real-time PCR, transwell assay, colony formation assay, and mammosphere assay to dissect the molecular events associated with the expression of HIST1H1B in human breast cancer. We studied the association of HIST1H1B with CSF2 by ChIP assay. Using tumorigenesis assays, we determine the effect of HIST1H1B expression on tumorigenicity of BLBC cells.

Results

Here, we show that the linker histone HIST1H1B is dramatically elevated in BLBC due to HIST1H1B copy number amplification and promoter hypomethylation. HIST1H1B upregulates colony-stimulating factor 2 (CSF2) expression by binding the CSF2 promoter. HIST1H1B expression promotes, whereas knockdown of HIST1H1B expression suppresses tumorigenicity. In breast cancer patients, HIST1H1B expression is positively correlated with large tumor size, high grade, metastasis and poor survival.

Conclusion

HIST1H1B contributes to basal-like breast cancer progression by modulating CSF2 expression, indicating a potential prognostic marker and therapeutic target for this disease.

Transcriptomic Profiling Predicts Multiple Pathways and Molecules Associated With the Metastatic Phenotype of Oral Cancer Cells

BACKGROUND/AIM

Metastasis to cervical lymph nodes of oral squamous cell carcinoma (OSCC) leads to a poor prognosis. The present study aimed at investigating the pathways and molecules associated with OSCC metastasis.

MATERIALS AND METHODS

The transcriptome between HSC-3 cells and their highly metastatic subline, HSC-3-M3 cells, was examined using gene expression microarray. Gene enrichment analyses and Ingenuity Pathway Analysis were performed. Kaplan-Meier plot analysis using a publicly available dataset was conducted to assess whether candidate molecules are prognosticators.

RESULTS

A total of 1,018 genes were differentially expressed, and the inflammatory pathway and NF-kB were predicted to be activated in HSC-3-M3 cells. CSF2 was suggested to be an indicator of poor prognosis in head and neck cancers.

CONCLUSION

Inflammation and NF-kB may be involved in the metastasis of OSCC, and CSF2 is a promising diagnostic and therapeutic molecule. Moreover, HSC-3-M3 cells are a useful cell line model for studying OSCC progression.

Identification of prognostic marker genes in head and neck squamous cell carcinoma: A study based on The Cancer Genome Atlas database and experimental validation

BACKGROUND

Early detection and prognostic prediction are crucial in improving the survival of patients with head and neck squamous cell carcinoma (HNSCC). Therefore, we provided potential molecular markers in this study for early diagnosis and prognosis of this cancer based on The Cancer Genome Atlas (TCGA) database analysis and experimental validations.

METHODS

Differentially expressed genes (DEGs) between HNSCC tumor and normal samples were identified by TCGA database-based analyses. Univariate and multivariate Cox regression analyses were applied, respectively, to identify survival-related DEGs and independent prognostic factors in HNSCC. Further, RT-qPCR was employed to verify expression of DEGs in cancer and adjacent tissues from HNSCC patients recruited in our hospital, in which we also clarified the correlation between candidate genes and clinicopathological characteristics and prognosis of HNSCC patients.

RESULTS

TCGA data analyses yielded 59 DEGs. Cox analyses identified 13 candidate genes closely related to prognosis of HNSCC patients and established a five-gene signature comprising AC103702.2, LINC00941, RPL29, FOXL2, and CCL11. This five-gene signature could classify patients into high- and low-risk groups. The survival rate of the high-risk group was significantly lower than that of the low-risk group. Clinical tissue experiments further confirmed that AC103702.2, LINC00941, CCL11, and RPL29P19 genes were inversely associated with the prognosis of HNSCC patients, while CCL11 gene was positively associated. We also found that high-risk HNSCC patients presented a higher incidence of lymph node metastasis.

CONCLUSION

Five prognostic marker genes (AC103702.2, LINC00941, CCL11, RPL29P19, and FOXL2) as a gene cluster may serve as prognostic marker genes in HNSCC.

Integrated analysis of DNA methylation and mRNA expression profiles to identify key genes in head and neck squamous cell carcinoma

DNA methylation has been demonstrated to play significant roles in the etiology and pathogenesis of head and neck squamous cell carcinoma (HNSCC). In the present study, methylation microarray dataset (GSE87053) and gene expression microarray dataset (GSE23558) were downloaded from GEO database and analyzed through R language. A total of 255 hypermethylated-downregulated genes and 114 hypomethylated-upregulated genes were finally identified. Functional enrichment analyses were performed and a comprehensive protein–protein interaction (PPI) network was constructed. Subsequently, the top ten hub genes selected by Cytoscape software were subjected to further analyses. It was illustrated that the expression level of CSF2, CTLA4, ETS1, PIK3CD, and CFTR was intimately associated with HNSCC. Survival analysis suggested that CTLA4 and FGFR2 could serve as effective independent prognostic biomarkers for HNSCC patients. Overall, our study lay a groundwork for further investigation into the underlying molecular mechanisms in HNSCC carcinogenesis, providing potential biomarkers and therapeutic targets for HNSCC.

A Novel Nine Apoptosis-Related Genes Signature Predicting Overall Survival for Kidney Renal Clear Cell Carcinoma and its Associations with Immune Infiltration

Background: This study was designed to establish a sensitive prognostic model based on apoptosis-related genes to predict overall survival (OS) in patients with clear cell renal cell carcinoma (ccRCC). Methods: Obtaining the expression of apoptosis-related genes and associated clinical parameters from online datasets (The Cancer Genome Atlas, TCGA), their biological function analyses were performed through differently expressed genes. By means of LASSO, unadjusted and adjusted Cox regression analyses, this predictive signature was constructed and validated by internal and external databases (both TCGA and ArrayExpress). Results: A total of nine apoptosis-related genes (SLC27A2, TNFAIP2, IFI44, CSF2, IL4, MDK, DOCK8, WNT5A, APP) were ultimately screened as associated hub genes and utilized to construct a prognosis model. Then our constructed riskScore model significantly passed the validation in both the internal and external datasets of OS (all p < 0.05) and verified their expressions by qRT-PCR. Moreover, we conducted the Receiver Operating Characteristic (ROC), finding the area under the ROC curves (AUCs) were all above 0.70 which indicated that riskScore was a stable independent prognostic factor (p < 0.05). Furthermore, prognostic nomograms were established to figure out the relationship between 1-, 3- and 5-year OS and individual parameters for ccRCC patients. Additionally, survival analyses indicated that our riskScore worked well in predicting OS in subgroups of age, gender, grade, stage, T, M, N0, White (all p < 0.05), except for African, Asian and N1 (p > 0.05). We also explored its association with immune infiltration and applied cMap database to seek out highly correlated small molecule drugs. Conclusion: Our study successfully constructed a prognostic model containing nine hub apoptosis-related genes for ccRCC, helping clinicians predict patients' OS and making the prognostic assessment more standardized. Future prospective studies are required to validate our findings.

Identification of intestinal flora-related key genes and therapeutic drugs in colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is a multifactorial tumor and a leading cause of cancer-specific deaths worldwide. Recent research has shown that the alteration of intestinal flora contributes to the development of CRC. However, the molecular mechanism by which intestinal flora influences the pathogenesis of CRC remains unclear. This study aims to explore the key genes underlying the effect of intestinal flora on CRC and therapeutic drugs for CRC.

METHODS

Intestinal flora-related genes were determined using text mining. Based on The Cancer Genome Atlas database, differentially expressed genes (DEGs) between CRC and normal samples were identified with the limma package of the R software. Then, the intersection of the two gene sets was selected for enrichment analyses using the tool Database for Annotation, Visualization and Integrated Discovery. Protein interaction network analysis was performed for identifying the key genes using STRING and Cytoscape. The correlation of the key genes with overall survival of CRC patients was analyzed. Finally, the key genes were queried against the Drug-Gene Interaction database to find drug candidates for treating CRC.

RESULTS

518 genes associated with intestinal flora were determined by text mining. Based on The Cancer Genome Atlas database, we identified 48 DEGs associated with intestinal flora, including 25 up-regulated and 23 down-regulated DEGs in CRC. The enrichment analyses indicated that the selected genes were mainly involved in cell-cell signaling, immune response, cytokine-cytokine receptor interaction, and JAK-STAT signaling pathway. The protein-protein interaction network was constructed with 13 nodes and 35 edges. Moreover, 8 genes in the significant cluster were considered as the key genes and chemokine (C-X-C motif) ligand 8 (CXCL8) correlated positively with the overall survival of CRC patients. Finally, a total of 24 drugs were predicted as possible drugs for CRC treatment using the Drug-Gene Interaction database.

CONCLUSIONS

These findings of this study may provide new insights into CRC pathogenesis and treatments. The prediction of drug-gene interaction is of great practical significance for exploring new drugs or novel targets for existing drugs.

Weighted Gene Coregulation Network Analysis of Promoter DNA Methylation on All-Cause Mortality in Old-Aged Birth Cohorts Finds Modules of High-Risk Associated Biomarkers

Overall or all-cause mortality is a key measure of health in a population. Multiple epigenome-wide association studies have been conducted on all-cause mortality with limited significant findings and low replication. To elucidate the coregulated DNA methylation patterns associated with all-cause mortality, we conducted a weighted DNA methylation coregulation network analysis on whole-blood samples of 1,425 older individuals from the Lothian Birth Cohorts of 1921 and 1936. Our network-based analysis defined coregulated DNA methylation patterns in gene promoters into clusters or modules whose correlation with all-cause mortality was assessed by survival analysis. We found two significant modules or gene clusters associated with all-cause mortality in LBC1921 based on their eigengenes; one negatively correlated (p = 8.14E-03, 698 genes) and one positively correlated (p = 4.26E-02, 1,431 genes) with the risk of death. The two modules were replicated in LBC1936 with the same directions of correlation (p = 6.35E-02 and p = 3.64E-02, respectively). Furthermore, the modules revealed 32 genes associated with all-cause mortality (FDR < 0.05) linked to various diseases, including cancer and diabetes. Additionally, we performed pathway analysis and found 22 pathways (FDR < 0.05), including a pathway for taste transduction, which has been shown to be associated with poor prognosis in acutely hospitalized patients, and several pathways were linked to different types of cancer. The results from our network analysis show that DNA methylation of multiple genes could have been coregulated in an association with the overall risk of death. The identified epigenetic markers might help with our understanding of the molecular basis of all-cause mortality and general health.

Targeting Stat3 signaling impairs the progression of bladder cancer in a mouse model

Bladder cancer is the fourth most commonly diagnosed malignancy in men worldwide and has one of the highest recurrence rates of all cancers. This cancer type is unique because chronic inflammation caused by Schistosoma haematobium can cause bladder cancer, while inflammation induced by Bacillus Calmette Guerin is the therapeutic cornerstone for this cancer type. Activation of proinflammatory IL-6/Stat3 axis promotes the development of different cancers by acting on cancer cells as well as by modulating cancer microenvironment. Using a genetic and pharmacological approach in a mouse model, we demonstrated the importance of IL-6 and Stat3 signaling in bladder cancer. Our findings show that pharmacological inhibition of Stat3 with WP1066 effectively delays progression and invasiveness of bladder cancer in N-butyl-N-(4-hydroxybutyl) nitrosamine-induced mouse model. Moreover, either IL-6 blockade or Stat3 inhibition sensitized bladder cancer to anti-PD-L1 immune therapy. Taken together, our study demonstrates an important role of IL-6/Stat3 signaling in bladder cancer and creates a rationale for testing the therapeutic potential of Stat3 inhibitors in human MIBC both alone or in combination with anti-PD-L1 and anti-IL-6 therapy.

Identification of novel genes in testicular cancer microenvironment based on ESTIMATE algorithm-derived immune scores

Testicular cancer is the most common solid malignancy among young men. We downloaded data of testicular cancer patients from The Cancer Genome Atlas database to find novel genes in the testicular cancer microenviroment based on ESTIMATE algorithm-derived immune scores. A total of 156 cases of testicular cancer were included in this study and 165 cases of normal testicular tissues were used. We divided the testicular cancer patients into high- and low-score groups based on their immune scores. We identified 1,226 differentially expressed genes (fold change > 2, false discovery rate < 0.05), including 688 downregulated genes and 538 upregulated genes, between these two groups. The top Gene Ontology terms were involved in the immune response-regulating cell surface receptor signaling pathway, immune response-activating cell surface receptor signaling pathway, external side of the plasma membrane, and receptor ligand activity. By performing the Kyoto Encyclopedia of Genes and Genomes analysis, we demonstrated that cAMP signaling pathway was highly enriched among these differentially expressed genes. High expression of LINC01564, LINC02208, ODAM, RNA5SP111, and RNU6-196P were found to be associated with poor overall survival. The expression of genes was further validated by the Human Protein Atlas and only ALB and IFNG were demonstrated to be differentially expressed between testis tissue and testicular cancer tissue.

An optimal prognostic model based on gene expression for clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the most prevalent type of RCC; however, prognostic prediction tools for ccRCC are scant. Developing mRNA or long non-coding RNA (lncRNA)-based risk assessment tools may improve the prognosis in patients with ccRCC. RNA-sequencing and prognostic data from patients with ccRCC were downloaded from The Cancer Genome Atlas and the European Bioinformatics Institute Array database at the National Center for Biotechnology Information. Differentially expressed (DE) RNAs (DERs) and prognostic DERs were screened between less favorable and favorable prognoses using the limma package in R 3.4.1, and analyzed using univariate and multivariate Cox regression analyses, respectively. Risk score models were constructed using optimal combinations of DEmRNAs and DElncRNAs identified using the Least Absolute Shrinkage And Selection Operator Cox regression model of the penalized package. Associations between risk score models and overall survival time were evaluated. Independent prognostic clinical factors were screened using univariate and multivariate Cox regression analyses, and nomogram models were constructed. Gene Ontology biological processes and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were conducted using the clusterProfiler package in R3.4.1. A total of 451 DERs were identified, including 404 mRNAs and 47 lncRNAs, between less favorable and favorable prognoses, and 269 DERs, including 233 mRNAs and 36 lncRNAs, were identified as independent prognostic factors. Optimal combinations including 10 DEmRNAs or 10 DElncRNAs were screened using four risk score models based on the status or expression levels of the 10 DEmRNAs or 10 DElncRNAs. The model based on the expression levels of the 10 DEmRNAs had the highest prognostic power. These prognostic DEmRNAs may be involved in biological processes associated with the inflammatory response, complement and coagulation cascades and neuroactive ligand-receptor interaction pathways. The present validated risk assessment tool based on the expression levels of these 10 DEmRNAs may help to identify patients with ccRCC at a high risk of mortality. These 10 DEmRNAs in optimal combinations may serve as prognostic biomarkers and help to elucidate the pathogenesis of ccRCC.

Bioinformatics Analysis of the Expression of Key Long Intergenic Non-Protein Coding RNA Genes in Bladder Cancer

Background

Evidence indicates that there is an important role for long non-coding RNAs (lncRNA) in numerous cellular processes and that lncRNAs dysregulation contributes to tumor progression. Improved insight into the molecular characteristics of bladder cancer is required to predict outcomes and to develop a new rationale for targeted therapeutic strategies. Bioinformatics methods, including functional enrichment and network analysis combined with survival analysis, are required to process a large volume of data to obtain further information about differentially expressed genes (DEGs) in bladder cancer. This study aimed to explore the role of lncRNAs and their regulation network in bladder cancer.

Material/Methods

We analyzed bladder cancer data by The Cancer Genome Atlas profiling to identify differentially expressed lncRNAs in bladder cancer. The genes involved in the circlncRNAnet database were evaluated using Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), evolutionary relationship analysis, and protein-protein interaction (PPI) networks.

Results

Two new lncRNAs, ADAMTS9-AS1 and LINC00460, were shown to be differentially expressed in bladder cancer. Patients were divided into 2 groups (high expression and low expression) according to their median expression values. The overall survival and disease-free survival of patients with high ADAMTS9-AS1 bladder cancer were significantly shorter; the expression of LINC00460 had no significant correlation with survival. GO and KEGG analysis of the 2 lncRNA-related genes revealed that these lncRNAs played a vital role in tumorigenesis. Bioinformatics analysis showed that key genes related to LINC00460, including CXCL, CCL, and CSF2, may be related to the development of bladder cancer. The low expression of ADAMTS9-AS1 may influence the survival rate of bladder cancer with the hub gene as a target.

Conclusions

LncRNA, including LINC00460 and ADAMTS9-AS1, might play a crucial role in the biosynthesis network of bladder cancer. Differential expression results of ADAMTS9-AS1 suggests it may be correlated with a worse prognosis and a shorter survival time. We outlined the biosynthesis network that regulates lncRNAs in bladder cancer. Further experimental data is needed to validate our results.

The Effect of Baicalin on the Expression Profiles of Long Non-Coding RNAs and mRNAs in Porcine Aortic Vascular Endothelial Cells Infected with Haemophilus parasuis

Haemophilus parasuis can elicit serious inflammatory responses, which contribute to huge economic losses to the swine industry. However, the pathogenic mechanisms underlying inflammation-related damage induced by H. parasuis remain unclear. Accumulating evidence indicates that long non-coding RNAs (lncRNAs) have important functions in the regulation of autoimmune disorders. Baicalin has been shown to have anti-inflammatory, anti-microbial, and anti-oxidant activities. In this study, we investigated whether lncRNAs were involved in the vascular injury or inflammation triggered by H. parasuis and whether baicalin regulated the lncRNA profiles of porcine aortic vascular endothelial cells (PAVECs) infected with H. parasuis. The results showed that the lncRNA and mRNA expression profiles of PAVECs were changed by H. parasuis. Important functions of lncRNAs and mRNAs were predicted. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses demonstrated that the targets of differentially expressed lncRNAs of H. parasuis infected PAVECs were mainly involved in the tumor necrosis factor (TNF) signaling pathway, apoptosis, and N-glycan biosynthesis; whereas nicotinate and nicotinamide metabolism, the cytosolic DNA-sensing pathway, the TNF signaling pathway, and the nuclear factor (NF)-kappa B signaling pathway were enriched in PAVECs pretreated with baicalin. In addition, top hub genes and lncRNAs were identified and validated by quantitative polymerase chain reaction. CCL5, GBP1, and SAMHD1 were significantly upregulated after H. parasuis infection, whereas they were significantly downregulated with baicalin pretreatment. LncRNA ALDBSSCT0000001677, ALDBSSCT0000001353, MSTRG.10724.2, and ALDBSSCT0000010434 had the same expression pattern. Collectively, these data suggested that baicalin could modify changes to the lncRNAs profiles or regulate lncRNAs that participate in inflammation-related signaling pathways, thereby alleviating tissue damage or inflammatory responses induced by H. parasuis. To our best knowledge, this is the first article of H. parasuis stimulating changes to the lncRNA profiles of PAVECs and the capability of baicalin to regulate lncRNA changes in PAVECs infected with H. parasuis, which might provide a novel therapeutic target for the control of H. parasuis infection.

Novel methylation gene panel in adjacent normal tissues predicts poor prognosis of colorectal cancer in Taiwan

BACKGROUND

It is evident that current clinical criteria are suboptimal to accurately estimate patient prognosis. Studies have identified epigenetic aberrant changes as novel prognostic factors for colorectal cancer (CRC).

AIM

To estimate whether a methylation gene panel in different clinical stages can reflect a different prognosis.

METHODS

We enrolled 120 CRC patients from Tri-Service General Hospital in Taiwan and used the candidate gene approach to select six genes involved in carcinogenesis pathways. Patients were divided into two groups based on the methylation status of the six evaluated genes, namely, the < 3 aberrancy group and ≥ 3 aberrancy group. Various tumor stages were divided into two subgroups (local and advanced stages) on the basis of the pathological type of the following tissues: Tumor and adjacent normal tissues (matched normal). We assessed DNA methylation in tumors and adjacent normal tissues from CRC patients and analyzed the association between DNA methylation with different cancer stages and the prognostic outcome including time to progression (TTP) and overall survival.

RESULTS

We observed a significantly increasing trend of hazard ratio as the number of hypermethylated genes increased both in normal tissue and tumor tissue. The 5-year TTP survival curves showed a significant difference between the ≥ 3 aberrancy group and the < 3 aberrancy group. Compared with the < 3 aberrancy group, a significantly shorter TTP was observed in the ≥ 3 aberrancy group. We further analyzed the interaction between CRC prognosis and different cancer stages (local and advanced) according to the methylation status of the selected genes in both types of tissues. There was a significantly shorter 5-year TTP for tumors at advanced stages with the promoter methylation status of selected genes than for those with local stages. We found an interaction between cancer stages and the promoter methylation status of selected genes in both types of tissues.

CONCLUSION

Our data provide a significant association between the methylation markers in normal tissues with advanced stage and prognosis of CRC. We recommend using these novel markers to assist in clinical decision-making.

Multiple gene promoter methylation and clinical stage in adjacent normal tissues: Effect on prognosis of colorectal cancer in Taiwan

This study provide an insight that the panel genes methylation status in different clinical stage tended to reflect a different prognosis even in matched normal tissues, to clinical recommendation. We enrolled 153 colorectal cancer patients from a medical center in Taiwan and used the candidate gene approach to select five genes involved in carcinogenesis pathways. We analyzed the relationship between DNA methylation with different cancer stages and the prognostic outcome. There were significant trends of increasing risk of 5-year time to progression and event-free survival of subjects with raising number of hypermethylation genes both in normal tissue and tumor tissue. The group with two or more genes with aberrant methylation in the advanced cancer stages (Me/advanced) had lower 5-year event-free survival among patients with colorectal cancer in either normal or tumor tissue. The adjusted hazard ratios in the group with two or more genes with aberrant methylation with advanced cancer stages (Me/advanced) were 8.04 (95% CI, 2.80–23.1; P for trend <0.01) and 8.01 (95% CI, 1.92–33.4; P for trend <0.01) in normal and tumor tissue, respectively. DNA methylation status was significantly associated with poor prognosis outcome. This finding in the matched normal tissues of colorectal cancer patients could be an alternative source of prognostic markers to assist clinical decision making.

Demethylation and Overexpression of CSF2 are Involved in Immune Response, Chemotherapy Resistance, and Poor Prognosis in Colorectal Cancer

Purpose

This study aimed to evaluate the role of colony-stimulating factor 2 (CSF2) in chemotherapy resistance, prognosis, and immune response and to identify its possible mechanisms underlying drug resistance.

Methods

Drug-resistant cell lines were obtained by successively increasing drug concentration. RNA-Seq was performed to screen hub genes. CSF2 expression was analyzed via immunohistochemistry. Moreover, The Cancer Genome Atlas (TCGA), Tumor Immune Estimation Resource (TIMER) dataset, and R2 platform were used to explore the correlations among CSF2 expression, prognosis, and immune response.

Results

RNA-Seq indicated that microRNAs in cancer, P53 signaling pathway, and cell cycle were associated with FOLFOX chemotherapy resistance. Protein-protein interaction (PPI), molecular complex detection (MOCDE), and qRT-PCR analysis verified CSF2 as the hub gene in chemotherapy resistance. Moreover, CSF2 expression was lower in the normal tissue than in the cancerous tissue (P<0.05). Higher expression of CSF2 was associated with poor OS and DFS in colon cancer patients (P<0.05). We further found similar results in the Oncomine database and R2 platform (P<0.05). A higher expression of CSF2 in the CRC tissue may be caused by demethylation, which was verified using the TCGA datasets. Moreover, GSEA demonstrated that CSF2 was associated with immune response, which was consistent with results reported using TIMER datasets.

Conclusion

CSF2 is a novel biomarker and a prognostic factor for the survival of CRC patients affecting the immune response, and an overexpression of CSF2 in CRC patients may be caused by DNA demethylation.

PDHA1 Gene Knockout In Human Esophageal Squamous Cancer Cells Resulted In Greater Warburg Effect And Aggressive Features In Vitro And In Vivo

Background

One of the remarkable metabolic characteristics of cancer cells is that they prefer glycolysis rather than oxidative phosphorylation (OXPHOS). Pyruvate dehydrogenase E1 alpha subunit (PDHA1) is an important prerequisite for OXPHOS. Our previous studies have shown that low level of PDHA1 protein expression in esophageal squamous cell cancer (ESCC) was correlated with poor prognosis. However, the effect of PDHA1 inhibition on metabolism and biological behavior of esophageal cancer cells remains unclear.

Methods And Results

In this study, a KYSE450 PDHA1 knockout (KO) cell line of esophageal cancer was established by CRISPR/Cas9 technology. Then, the glycose metabolism, cell proliferation and migration abilities, chemotherapeutic tolerance and angiogenesis of the PDHA1 KO cells were investigated in vitro and in vivo. In the PDHA1 KO cells, the glycolysis and the consumption of glucose and glutamine were significantly enhanced, while the OXPHOS was significantly suppressed, implying Warburg effect in the PDHA1 KO cells. Furthermore, it was also proved in vitro experiments that the PDHA1 KO cell obtained proliferation advantage, as well as significantly greater chemotherapy tolerance and migration ability. Xenograft experiments discovered not only larger tumors but also increased angiogenesis in the PDHA1 KO cell group.

Conclusion

Inhibition of PDHA1 gene expression in human ESCC leads to metabolic reprogramming of Warburg effect and increased malignancies. Targeting ESCC metabolic reprogramming may become a potential therapeutic target.

Knockdown of RAD18 inhibits glioblastoma development

This study aimed at investigating the role of RAD18 in the regulation of glioblastoma development as well as the underlying mechanisms. The human glioblastoma U251 and U87MG cells were transfected with siRNAs specifically targeting RAD18, and the effects of knockdown of RAD18 on the viability, apoptosis, migration, and invasion of U251 and U87MG cells were investigated. Transcriptome sequencing of the siRNA-RAD18-tranfected and siRNA-NC-transfected U251 cells was performed, followed by bioinformatic analyses for sequencing data. The results showed that knockdown of RAD18 significantly inhibited cell viability, promoted apoptosis, and suppressed migration and invasion of U251 and U87MG cells. Bioinformatic analyses of sequencing data identified 1,051 differentially expressed genes (DEGs) (369 up- and 682 downregulated genes) in the siRNA-RAD18-transfected U251 cells compared with siRNA-NC-transfected U251 cells. Eleven DEGs, including nerve growth factor (NGF), colony-stimulating factor 2 (CSF2), matrix metallopeptidase 1 (MMP1), platelet-derived growth factor receptor α (PDGFRA), and heme oxygenase 1 (HMOX1), were identified as the hub nodes in protein-protein interaction (PPI) network. Moreover, the aforementioned 11 hub genes were significantly enriched in PI3K-Akt signaling pathway and GO functions associated with the extracellular region. Notably, quantitative real-time polymerase chain reaction further confirmed that the expression levels of NGF, CSF2, HMOX1, and MMP1 were significantly downregulated, while that of PDGFRA was markedly upregulated in the siRNA-RAD18-transfected U251 cells than in the siRNA-NC cells. In conclusion, the knockdown of RAD18 may inhibit glioblastoma development by regulating the expression of the aforementioned key DEGs.

MCM10 overexpression implicates adverse prognosis in urothelial carcinoma

Urothelial carcinoma (UC) occurs in the upper urinary tract (UTUC) and the urinary bladder (UBUC). The molecular pathogenesis of UC has not been fully elucidated. Through data mining of a published transcriptome of UBUC (GSE31684), we identified Minichromosome Maintenance Complex Component 2 (MCM2) and MCM10 as the two most significantly upregulated genes in UC progression among the MCM gene family, the key factors for the initiation of DNA replication. To validate the clinical significance of MCM2 and MCM10, immunohistochemistry, evaluated by H-score, was used in a pilot study of 50 UTUC and 50 UBUC samples. Only a high expression level of MCM10 predicted worse disease-specific survival (DSS) and inferior metastasis-free survival (MeFS) for both UTUC and UBUC. Correspondingly, evaluation of MCM10 mRNA expression in 36 UTUCs and 30 UBUCs showed significantly upregulated levels in high stage UC, suggesting its role in tumor progression. Evaluation of 340 UTUC and 296 UBUC tissue samples, respectively, demonstrated that high MCM10 immunoexpression was significantly associated with advanced primary tumors, nodal status, and the presence of vascular invasion in both groups of UCs. In multivariate Cox regression analyses, adjusted for standard clinicopathological features, MCM10 overexpression was independently associated with DSS (UTUC hazard ratio [HR]=2.401, P = 0.013; UBUC HR=4.323, P=0.001) and with MeFS (UTUC HR=3.294, P<0.001; UBUC HR=1.972, P=0.015). In vitro, knockdown of MCM10 gene significantly suppressed cell proliferation in both J82 and TCCSUP cells. In conclusion, MCM10 overexpression was associated with unfavorable clinicopathological characteristics and independent negative prognostic effects, justifying its potential theranostic value in UC.

Determination of Genes Related to Uveitis by Utilization of the Random Walk with Restart Algorithm on a Protein-Protein Interaction Network

Uveitis, defined as inflammation of the uveal tract, may cause blindness in both young and middle-aged people. Approximately 10–15% of blindness in the West is caused by uveitis. Therefore, a comprehensive investigation to determine the disease pathogenesis is urgent, as it will thus be possible to design effective treatments. Identification of the disease genes that cause uveitis is an important requirement to achieve this goal. To begin to answer this question, in this study, a computational method was proposed to identify novel uveitis-related genes. This method was executed on a large protein–protein interaction network and employed a popular ranking algorithm, the Random Walk with Restart (RWR) algorithm. To improve the utility of the method, a permutation test and a procedure for selecting core genes were added, which helped to exclude false discoveries and select the most important candidate genes. The five-fold cross-validation was adopted to evaluate the method, yielding the average F1-measure of 0.189. In addition, we compared our method with a classic GBA-based method to further indicate its utility. Based on our method, 56 putative genes were chosen for further assessment. We have determined that several of these genes (e.g., CCL4, Jun, and MMP9) are likely to be important for the pathogenesis of uveitis.

High Expression of EphA4 Predicted Lesser Degree of Tumor Regression after Neoadjuvant Chemoradiotherapy in Rectal Cancer

Background: Numerous transmembrane receptor tyrosine kinase pathways have been found to play an important role in tumor progression in some cancers. This study was aimed to evaluate the clinical impact of Eph receptor A4 (EphA4) in patients with rectal cancer treated with neoadjuvant concurrent chemoradiotherapy (CCRT) combined with mesorectal excision, with special emphasis on tumor regression.

Methods: Analysis of the publicly available expression profiling dataset of rectal cancer disclosed that EphA4 was the top-ranking, significantly upregulated, transmembrane receptor tyrosine kinase pathway-associated gene in the non-responders to CCRT, compared with the responders. Immunohistochemical study was conducted to assess the EphA4 expression in pre-treatment biopsy specimens from 172 rectal cancer patients without distant metastasis. The relationships between EphA4 expression and various clinicopathological factors or survival were statistically analyzed.

Results: EphA4 expression was significantly associated with vascular invasion (P=0.015), post-treatment depth of tumor invasion (P=0.006), pre-treatment and post-treatment lymph node metastasis (P=0.004 and P=0.011, respectively). More importantly, high EphA4 expression was significantly predictive for lesser degree of tumor regression after CCRT (P=0.031). At univariate analysis, high EphA4 expression was a negative prognosticator for disease-specific survival (P=0.0009) and metastasis-free survival (P=0.0001). At multivariate analysis, high expression of EphA4 still served as an independent adverse prognostic factor for disease-specific survival (HR, 2.528; 95% CI, 1.131-5.651; P=0.024) and metastasis-free survival (HR, 3.908; 95% CI, 1.590-9.601; P=0.003).

Conclusion: High expression of EphA4 predicted lesser degree of tumor regression after CCRT and served as an independent negative prognostic factor in patients with rectal cancer.

High Expression of 17β-hydroxysteroid Dehydrogenase Type 2 is Associated with a Better Prognosis in Urothelial Carcinoma of the Urinary Tract

Object: To identify associations between 17β-hydroxysteroid dehydrogenase type 2 (HSD17B2) expression and clinicopathological variables and prognoses in patients with urothelial carcinoma of the urinary tract.

Materials and Methods: We analyzed one public transcriptome dataset (GSE31684) from the Gene Expression Omnibus. HSD17B2 showed the highest log2-transformed fold-change in expression, and it was therefore further analyzed. We enrolled 340 upper urinary tract and 295 urinary bladder tissue samples that were obtained from patients with urothelial carcinoma between 1996 and 2004 to evaluate the expression of HSD17B2 using immunohistochemistry. The endpoints were disease-specific survival and metastasis-free survival. Univariate and multivariate analyses were performed to assess the relationships between HSD17B2, survival and clinicopathological parameters.

Results: High expression of HSD17B2 was significantly associated with better clinicopathological parameters, including the following parameters in urothelial carcinoma of upper urinary tracts: Ta (non-invasive papillary carcinoma) and T1 (invaded subepithelial connective tissue) disease, without nodal metastasis, without vascular invasion, without perineal invasion and with a low histological grade, and the following parameters specifically in urothelial carcinoma of the urinary bladder: Ta (noninvasive papillary carcinoma) and T1 (invaded into the lamina propria but not into the muscularis propria) disease, without nodal metastasis (all P<0.05). Additionally, HSD17B2 high expression predicted a better prognosis, including improved disease-specific survival and metastasis-free survival in urothelial carcinomas of the urinary tract system.

Conclusions: High expression of HSD17B2 is associated with a better prognosis and is therefore a prognostic biomarker that can be used to predict favorable survival in patients with urothelial carcinoma of the urinary tract system.

Downregulation of RNF128 Predicts Progression and Poor Prognosis in Patients with Urothelial Carcinoma of the Upper Tract and Urinary Bladder

Background: The TP53 tumor suppressor gene plays a crucial role in the carcinogenesis of many malignancies, including urothelial carcinoma (UC). Overexpression of p53 is associated with poor prognosis in UC. Recently, RING finger protein 128 (RNF128) was shown to be involved in p53-induced apoptosis, forming a negative feedback loop. However, the significance of RNF128 in patients with UC remains unknown. In this study, our aim was to evaluate the expression of RNF128 in UC and to assess its predictive and prognostic value in a well-established cohort. Methods: Through data mining from a published transcriptome (GSE31684), RNF128 was identified as the most differentially expressed gene in UC among those associated with negative regulation of the cytokine biosynthetic process (GO:0042036). Its immunoexpression was further evaluated using the H-scores of 340 patients with upper urinary tract UC (UTUC) and 295 with urinary bladder UC (UBUC). The scores were correlated with clinicopathological features, disease-specific survival (DSS) and metastasis-free survival (MeFS). We also used Western blot analysis to evaluate RNF128 protein expression in human urothelial cell (HUC) lines. Results: Downregulation of RNF128 expression was significantly associated with advanced pT stage (p<0.001), high histological grade (UTUC, p<0.001; UBUC, p=0.035), nodal metastasis (UTUC, p<0.001; UBUC, p=0.001), vascular invasion (UTUC, p<0.001; UBUC, p=0.008) and high mitotic rate (UTUC, p=0.003; UBUC, p=0.023). Low expression of RNF128 was an adverse prognosticator for DSS (UTUC, p<0.0001; UBUC, p<0.0001) and MeFS (UTUC, p<0.0001; UBUC, p=0.0002). Moreover, low expression was predictive of poor DSS (UTUC, p=0.006; UBUC, p=0.003) and MeFS (UTUC, p=0.009; UBUC, p=0.036) in multivariate comparisons. Western blot analysis showed that the RNF128 protein was downregulated in invasive urothelial cancer cell lines. Conclusion: Our findings showed that downregulation of RNF128 was correlated with cancer invasiveness and metastasis as well as reduced survival in patients with UTUC and UBUC, identifying RNF128 as a prognostic factor in UC.

Predictive models and prognostic factors for upper tract urothelial carcinoma: a comprehensive review of the literature

In the context of customized patient care for upper tract urothelial carcinoma (UTUC), decision-making could be facilitated by risk assessment and prediction tools. The aim of this study was to provide a critical overview of existing predictive models and to review emerging promising prognostic factors for UTUC. A literature search of articles published in English from January 2000 to June 2016 was performed using PubMed. Studies on risk group stratification models and predictive tools in UTUC were selected, together with studies on predictive factors and biomarkers associated with advanced-stage UTUC and oncological outcomes after surgery. Various predictive tools have been described for advanced-stage UTUC assessment, disease recurrence and cancer-specific survival (CSS). Most of these models are based on well-established prognostic factors such as tumor stage, grade and lymph node (LN) metastasis, but some also integrate newly described prognostic factors and biomarkers. These new prediction tools seem to reach a high level of accuracy, but they lack external validation and decision-making analysis. The combinations of patient-, pathology- and surgery-related factors together with novel biomarkers have led to promising predictive tools for oncological outcomes in UTUC. However, external validation of these predictive models is a prerequisite before their introduction into daily practice. New models predicting response to therapy are urgently needed to allow accurate and safe individualized management in this heterogeneous disease.

Prognostic Impact of Thrombospodin-2 (THBS2) Overexpression on Patients with Urothelial Carcinomas of Upper Urinary Tracts and Bladders

Purpose: Urothelial carcinoma (UC) is a type of tumor, especially of the urinary bladder, that affects people worldwide. Clarification of its detailed tumor biology and discovery of potential targets for developing treatment strategies are imperative because of frequent recurrences and poor prognosis of advanced UCs. By data mining a published dataset of UC of bladder (UCB) transcriptome (GSE31684) from Gene Expression Omnibus, National Center of Biotechnology Information (GEO, NCBI), we identified that THBS2 was the most significantly upregulated gene among those related to structural molecule activity (GO:0005198). Therefore, we evaluated the clinical significance and prognostic impact of thrombospondin-2 (THBS2) protein, A.K.A. TSP2, which encoded by THBS2 gene.

Materials and Methods: THBS2 immunostaining was performed in 340 UCs of upper urinary tract (UC-UUTs) and 295 UCBs; subsequently, both groups were dichotomized into high- and low-expression subgroups. Moreover, statistical analyses were performed to correlate the association between THBS2 expression and clinicopathological parameters with two survival indexes: disease-specific survival (DSS) and metastasis-free survival (MeFS).

Results: High THBS2 immunoexpression was significantly associated with advanced primary tumor status, nodal metastasis, and vascular invasion in both UC-UUT and UCB groups (all P ≤ .001). In addition, THBS2 overexpression was linked to adverse DSS and MeFS in univariate analyses and served as an independent prognosticator indicating poor outcomes in both groups in multivariate analyses.

Conclusion: THBS2 may play a crucial role in UC progression and may be a novel prognostic marker. Additional investigations to elucidate the molecular pathway are necessary for developing potential THBS2-targeted therapies for UCs.