Fifteen hub genes associated with progression and prognosis of clear cell renal cell carcinoma identified by coexpression analysis

Explorations of novel MDR-related hub genes and the potential roles TRIM9 played in drug-resistant hepatocellular carcinoma

Current chemotherapeutic efficacy is limited by the rapid development of multidrug resistance (MDR) in hepatocellular carcinoma (HCC). In this study, 66 MDR-related hub genes in drug-resistant HCC were identified through combined analysis of differential expressed genes (DEGs), gene functional enrichment, Cox proportional regression, weighted gene co-expression network analysis (WGCNA) and protein-protein interaction (PPI) network construction. A prognostic risk model was established through the LASSO-Cox regression analysis. Based on the comparison of gene mutation frequency, tumor mutation burden (TMB) and immune infiltration in high- and low-risk groups, we explored the relationships between the MDR-related hub genes and immune regulation. The competitive endogenous RNA (ceRNA) network and associated non-coding RNAs (ncRNAs) were predicted to investigate the potential mechanisms. Five MDR-related hub genes in drug-resistant HCC were finally confirmed, namely ABCB6, FLNC, MCC, NAV3 and TRIM9. TRIM9 was identified as the most significant gene enhancing MDR. Inhibiting TRIM9 caused a decrease in the IC50 of doxorubicin (DOX), and significant increases in the intracellular uptake, retention and absorption of DOX in HepG2/ADR cells. These findings may provide new insights into the mechanism of MDR development. The MDR-related hub genes, especially TRIM9 may be targeted therapeutically to enhance the prognosis of patients with drug-resistant HCC.

Identification of SV2C and DENR as Key Biomarkers for Parkinson's Disease Based on Bioinformatics, Machine Learning, and Experimental Verification

The objective of this study is to investigate the potential biomarkers and therapeutic target genes for Parkinson's disease (PD). We analyzed four datasets (GSE8397, GSE20292, GSE20163, GSE20164) from the Gene Expression Omnibus database. We employed weighted gene co-expression network analysis and differential expression analysis to select genes and perform functional analysis. We applied three algorithms, namely, random forest, support vector machine recursive feature elimination, and least absolute shrinkage and selection operator, to identify hub genes, perform functional analysis, and assess their clinical diagnostic potential using receiver operating characteristic (ROC) curve analysis. We employed the xCell website to evaluate differences in the composition patterns of immune cells in the GEO datasets. We also collected serum samples from PD patients and established PD cell model to validate the expression of hub genes using enzyme-linked immunosorbent assay and quantitative real-time polymerase chain reaction. Our findings identified SV2C and DENR as two hub genes for PD and decreased in PD brain tissue compared with controls. ROC analysis showed effectively value of SV2C and DENR to diagnose PD, and they were downregulated in the serum of PD patients and cell model. Functional analysis revealed that dopamine vesicle transport and synaptic vesicle recycling are crucial pathways in PD. Besides, the differences in the composition of immune cells, especially basophils and T cells, were discovered between PD and controls. In summary, our study identifies SV2C and DENR as potential biomarkers for diagnosing PD and provides a new perspective for exploring the molecular mechanisms of PD.

Multiomics characterization and verification of clear cell renal cell carcinoma molecular subtypes to guide precise chemotherapy and immunotherapy

Clear cell renal cell carcinoma (ccRCC) is a heterogeneous tumor with different genetic and molecular alterations. Schemes for ccRCC classification system based on multiomics are urgent, to promote further biological insights. Two hundred and fifty-five ccRCC patients with paired data of clinical information, transcriptome expression profiles, copy number alterations, DNA methylation, and somatic mutations were collected for identification. Bioinformatic analyses were performed based on our team's recently developed R package "MOVICS." With 10 state-of-the-art algorithms, we identified the multiomics subtypes (MoSs) for ccRCC patients. MoS1 is an immune exhausted subtype, presented the poorest prognosis, and might be caused by an exhausted immune microenvironment, activated hypoxia features, but can benefit from PI3K/AKT inhibitors. MoS2 is an immune "cold" subtype, which represented more mutation of VHL and PBRM1, favorable prognosis, and is more suitable for sunitinib therapy. MoS3 is the immune "hot" subtype, and can benefit from the anti-PD-1 immunotherapy. We successfully verified the different molecular features of the three MoSs in external cohorts GSE22541, GSE40435, and GSE53573. Patients that received Nivolumab therapy helped us to confirm that MoS3 is suitable for anti-PD-1 therapy. E-MTAB-3267 cohort also supported the fact that MoS2 patients can respond more to sunitinib treatment. We also confirm that SETD2 is a tumor suppressor in ccRCC, along with the decreased SETD2 protein level in advanced tumor stage, and knock-down of SETD2 leads to the promotion of cell proliferation, migration, and invasion. In summary, we provide novel insights into ccRCC molecular subtypes based on robust clustering algorithms via multiomics data, and encourage future precise treatment of ccRCC patients.

Prognostic and therapeutic implications of iron-related cell death pathways in acute myeloid leukemia

Acute myeloid leukemia (AML) is a blood cancer that is diverse in terms of its molecular abnormalities and clinical outcomes. Iron homeostasis and cell death pathways play crucial roles in cancer pathogenesis, including AML. The objective of this study was to examine the clinical significance of genes involved in iron-related cell death and apoptotic pathways in AML, with the intention of providing insights that could have prognostic implications and facilitate the development of targeted therapeutic interventions. Gene expression profiles, clinical information, and molecular alterations were integrated from multiple datasets, including TCGA-LAML and GSE71014. Our analysis identified specific molecular subtypes of acute myeloid leukemia (AML) displaying varying outcomes, patterns of immune cell infiltration, and profiles of drug sensitivity for targeted therapies based on the expression of genes involved in iron-related apoptotic and cell death pathways. We further developed a risk model based on four genes, which demonstrated promising prognostic value in both the training and validation cohorts, indicating the potential of this model for clinical decision-making and risk stratification in AML. Subsequently, Western blot analysis showed that the expression levels of C-Myc and CyclinD1 were significantly reduced after CD4 expression levels were knocked down. The findings underscore the potential of iron-related cell death pathways as prognostic biomarkers and therapeutic targets in AML, paving the way for further research aimed at understanding the molecular mechanisms underlying the correlation between iron balance, apoptosis regulation, and immune modulation in the bone marrow microenvironment.

Screening of Prognostic Markers for Hepatocellular Carcinoma Patients Based on Multichip Combined Analysis

Methods

GSE (14520, 36376, 57957, 76427) datasets were accessed from GEO database. 55 differential mRNAs (DEGs) were obtained by differential analysis based on the datasets. GO and KEGG analysis results indicated that the DEGs were enriched in xenobiotic metabolic process and other pathways. Expression profiles and clinical data of TCGA-LIHC mRNAs were from TCGA database. We established a prognostic model of HCC through univariate and multivariate Cox risk regression analyses. ROC curve analysis was used to examine the prognostic model performance. GSEA analysis was performed between the high- and low-risk score sample groups.

Results

A 4-gene HCC prognostic model was constructed, in which the gene expressions correlated to HCC patients' survival. The AUC value presented 0.734 in the ROC analysis for the prognostic model.

Conclusion

The four-gene model could be introduced as an independent prognostic factors to assess HCC patients' survival status.

The Hippo Pathway Effector Transcriptional Co-activator With PDZ-Binding Motif Correlates With Clinical Prognosis and Immune Infiltration in Colorectal Cancer

The transcriptional co-activator with PDZ-binding motif (TAZ) is a downstream effector of the Hippo pathway. It has been identified as an oncogene in certain tumor types; however, the function and role of TAZ in colorectal cancer (CRC) has not been illustrated. Here, we aimed to analyze the expression and role of TAZ in CRC. In this study, we investigated the expression level of TAZ in 127 CRC and matched adjacent normal tissues by immunohistochemistry (IHC) and analyzed its correlation with clinicopathological characteristics in CRC. Moreover, we further analyzed the role of TAZ in the CRC-associated immunology using integrative bioinformatic analyses. The cBioPortal and WebGestalt database were used to analyze the co-expressed genes and related pathways of TAZ in CRC by gene ontology (GO) and KEGG enrichment analyses. Meanwhile, the correlations between TAZ and the infiltrating immune cells and gene markers were analyzed by TIMER database. Our study revealed that TAZ expression is higher in CRC tissues than in matched adjacent non-tumor tissues. In addition, CRC patients with higher TAZ expression demonstrated poor overall survival (OS) and recurrent-free survival rates as compared to CRC patients with lower expression of TAZ. Furthermore, the TAZ expression was identified to closely associate with the immune infiltration of CD4 + T, CD8 + T, and B cells. Taken together, our findings suggest that TAZ may serve as a promising prognostic biomarker and therapeutic target in CRC.

Identification of Four Novel Prognostic Biomarkers and Construction of Two Nomograms in Adrenocortical Carcinoma: A Multi-Omics Data Study via Bioinformatics and Machine Learning Methods

Background: Adrenocortical carcinoma (ACC) is an orphan tumor which has poor prognoses. Therefore, it is of urgent need for us to find candidate prognostic biomarkers and provide clinicians with an accurate method for survival prediction of ACC via bioinformatics and machine learning methods.

Methods: Eight different methods including differentially expressed gene (DEG) analysis, weighted correlation network analysis (WGCNA), protein-protein interaction (PPI) network construction, survival analysis, expression level comparison, receiver operating characteristic (ROC) analysis, and decision curve analysis (DCA) were used to identify potential prognostic biomarkers for ACC via seven independent datasets. Linear discriminant analysis (LDA), K-nearest neighbor (KNN), support vector machine (SVM), and time-dependent ROC were performed to further identify meaningful prognostic biomarkers (MPBs). Cox regression analyses were performed to screen factors for nomogram construction.

Results: We identified nine hub genes correlated to prognosis of patients with ACC. Furthermore, four MPBs (ASPM, BIRC5, CCNB2, and CDK1) with high accuracy of survival prediction were screened out, which were enriched in the cell cycle. We also found that mutations and copy number variants of these MPBs were associated with overall survival (OS) of ACC patients. Moreover, MPB expressions were associated with immune infiltration level. Two nomograms [OS-nomogram and disease-free survival (DFS)-nomogram] were established, which could provide clinicians with an accurate, quick, and visualized method for survival prediction.

Conclusion: Four novel MPBs were identified and two nomograms were constructed, which might constitute a breakthrough in treatment and prognosis prediction of patients with ACC.

Novel Diagnostic Value of Driver Gene Transcription Signatures to Characterise Clear Cell Renal Cell Carcinoma, ccRCC

Routine molecular tumour diagnostics are augmented by DNA-based qualitative and quantitative molecular techniques detecting mutations of DNA. However, in the past decade, it has been unravelled that the phenotype of cancer, as it’s an extremely complex disease, cannot be fully described and explained by single or multiple genetic variants affecting only the coding regions of the genes. Moreover, studying the manifestation of these somatic mutations and the altered transcription programming—driven by genomic rearrangements, dysregulation of DNA methylation and epigenetic landscape—standing behind the tumorigenesis and detecting these changes could provide a more detailed characterisation of the tumour phenotype. Consequently, novel comparative cancer diagnostic pipelines, including DNA- and RNA-based approaches, are needed for a global assessment of cancer patients. Here we report, that by monitoring the expression patterns of key tumour driver genes by qPCR, the normal and the tumorous samples can be separated into distinct categories. Furthermore, we also prove that by examining the transcription signatures of frequently affected genes at 3p25, 3p21 and 9p21.3 genomic regions, the ccRCC (clear cell renal cell carcinoma) and non-tumorous kidney tissues can be distinguished based on the mRNA level of the selected genes. Our results open new diagnostics possibilities where the mRNA signatures of tumour drivers can supplement the DNA-based approaches providing a more precise diagnostics opportunity leading to determine more precise therapeutic protocols.

Downregulation of FXYD2 Is Associated with Poor Prognosis and Increased Regulatory T Cell Infiltration in Clear Cell Renal Cell Carcinoma

BACKGROUND

FXYD2, a gene coding for the γ subunit of Na⁺/K⁺-ATPase, was demonstrated to involve in carcinogenesis recently. However, the specific role of FXYD2 in clear cell renal cell carcinoma (ccRCC) remains unknown. The current study was conducted to investigate the expression, biological function, and potentially immune-related mechanisms of FXYD2 in ccRCC. Materials and methods. The data from TCGA-KIRC, ICGC, GEO, Oncomine, ArrayExpress, TIMER, HPA datasets, and our clinical samples were used to determine and validate the expression level, prognostic roles, and potentially immune-related mechanisms in ccRCC. Cell function assays were performed to investigate the biological role of FXYD2 in vitro.

RESULTS

FXYD2 was identified to be downregulated in ccRCC tissue compared to normal tissue, which was confirmed by our RT-PCR, WB, and IHC analyses. Kaplan-Meier survival analysis and Cox regression analysis suggested that downregulated FXYD2 could independently predict poor survival of ccRCC patients. Through the ESTIMATE algorithm, ssGSEA algorithm, CIBERSORT algorithm, TIMER database, and our laboratory experiment, FXYD2 was found to correlate with the immune landscape, especially regulatory T cells (Treg), in ccRCC. Gain-of-function experiment revealed that FXYD2 could restrain cell proliferation, migration, and invasion in vitro. Functional enrichment analysis illustrated that TGF-β-SMAD2/3, Notch, and PI3K-Akt-mTOR signaling pathways may be potential signaling pathways of FXYD2 in ccRCC.

CONCLUSIONS

Downregulation of FXYD2 is associated with ccRCC tumorigenesis, poor prognosis, and increased Treg infiltration in ccRCC, which may be related to TGF-β-SMAD2/3, Notch, and PI3K-Akt-mTOR signaling pathways. This will probably provide a novel prognostic marker and potential therapeutic target for ccRCC.

Identification and Validation of CYBB, CD86, and C3AR1 as the Key Genes Related to Macrophage Infiltration of Gastric Cancer

Gastric cancer (GC) is rampant around the world. Most of the GC cases are detected in advanced stages with poor prognosis. The identification of marker genes for early diagnosis is of great significance. Studying the tumor environment is helpful to acknowledge the process of tumorigenesis, development, and metastasis. Twenty-two kinds of immune cells were calculated by CIBERSORT from Gene Expression Omnibus (GEO) database. Subsequently, higher infiltration of macrophages M0 was discovered in GC compared with normal tissues. WGCNA was utilized to construct the network and then identify key modules and genes related to macrophages in TCGA. Finally, 18 hub genes were verified. In the PPI bar chart, the top 3 genes were chosen as hub genes involved in most pathways. On the TIMER and THPA websites, it is verified that the expression levels of CYBB, CD86, and C3AR1 genes in tumor tissues were higher than those in normal tissues. These genes may work as biomarkers or targets for accurate diagnosis and treatment of GC in the future. Our findings may be a new strategy for the treatment of GC.

Identification of the Potential Prognosis Biomarkers in Hepatocellular Carcinoma: An Analysis Based on WGCNA and PPI

Aim

This study was done to determine biomarkers for the prognostic prediction of hepatocellular carcinoma (HCC).

Materials and Methods

In the Gene Expression Omnibus, the gene expression profiles of HCC were downloaded. Biomarkers were identified by weighted gene co-expression network analysis and protein–protein interaction network analysis.

Results

There were 24 modules, which were characterized by the high correlation with HCC. Meanwhile, through enrichment analysis, differentially expressed genes were largely participated in the ubiquitination and autophagy processes. Moreover, PRC1, TOP2A and CKAP2L may be the hub genes involved in HCC tumorigenesis, and their biomarker roles were further demonstrated via Gene Expression Profiling Interactive Analysis (GEPIA) and Oncomine databases. In addition, the levels of PRC1, TOP2A and CKAP2L were obviously up-regulated in the sera of HCC patients.

Conclusion

PRC1, TOP2A and CKAP2L may serve as biomarkers for the prognostic prediction of HCC patients.

Screening Hub Genes of Hepatocellular Carcinoma Based on Public Databases

Tumor recurrence and metastasis often occur in HCC patients after surgery, and the prognosis is not optimistic. Hence, searching effective biomarkers for prognosis of is of great importance. Firstly, HCC-related data was acquired from the TCGA and GEO databases. Based on GEO data, 256 differentially expressed genes (DEGs) were obtained firstly. Subsequently, to clarify function of DEGs, clusterProfiler package was used to conduct functional enrichment analyses on DEGs. Protein-protein interaction (PPI) network analysis screened 20 key genes. The key genes were filtered via GEPIA database, by which 11 hub genes (F9, CYP3A4, ASPM, AURKA, CDC20, CDCA5, NCAP, PRC1, PTTG1, TOP2A, and KIFC1) were screened out. Then, univariate Cox analysis was applied to construct a prognostic model, followed by a prediction performance validation. With the risk score calculated by the model and common clinical features, univariate and multivariate analyses were carried out to assess whether the prognostic model could be used independently for prognostic prediction. In conclusion, the current study screened HCC prognostic gene signature based on public databases.

Immunogenomic Analyses of the Prognostic Predictive Model for Patients With Renal Cancer

Background

Renal cell carcinoma (RCC) is associated with poor prognostic outcomes. The current stratifying system does not predict prognostic outcomes and therapeutic benefits precisely for RCC patients. Here, we aim to construct an immune prognostic predictive model to assist clinician to predict RCC prognosis.

Methods

Herein, an immune prognostic signature was developed, and its predictive ability was confirmed in the kidney renal clear cell carcinoma (KIRC) cohorts based on The Cancer Genome Atlas (TCGA) dataset. Several immunogenomic analyses were conducted to investigate the correlations between immune risk scores and immune cell infiltrations, immune checkpoints, cancer genotypes, tumor mutational burden, and responses to chemotherapy and immunotherapy.

Results

The immune prognostic signature contained 14 immune-associated genes and was found to be an independent prognostic factor for KIRC. Furthermore, the immune risk score was established as a novel marker for predicting the overall survival outcomes for RCC. The risk score was correlated with some significant immunophenotypic factors, including T cell infiltration, antitumor immunity, antitumor response, oncogenic pathways, and immunotherapeutic and chemotherapeutic response.

Conclusions

The immune prognostic, predictive model can be effectively and efficiently used in the prediction of survival outcomes and immunotherapeutic responses of RCC patients.

Analysis and verification of N6-methyladenosine-modified genes as novel biomarkers for clear cell renal cell carcinoma

N⁶-methyladenosine (m⁶A) has been involved in diverse biological processes in cancer, but its function and clinical value in clear cell renal cell carcinoma (ccRCC) remain largely unknown. In this study, we found that 1453 m⁶A-modified differentially expressed genes (DEGs) of ccRCC were mainly enriched in cell cycle, PI3K-AKT, and p53 signaling pathways. Then we constructed a co-expression network of the 1453 m⁶A-modified DEGs and identified a most clinically relevant module, where NUF2, CDCA3, CKAP2L, KIF14, and ASPM were hub genes. NUF2, CDCA3, and KIF14 could combine with a major RNA m⁶A methyltransferase METTL14, serving as biomarkers for ccRCC. Real-time quantitative PCR assay confirmed that NUF2, CDCA3, and KIF14 were highly expressed in ccRCC cell lines and ccRCC tissues. Furthermore, these three genes were modified by m⁶A and negatively regulated by METTL14. This study revealed that NUF2, CDCA3, and KIF14 were m⁶A-modified biomarkers, representing a potential diagnostic, prognostic, and therapeutic target for ccRCC.

Abbreviations: m⁶A: N⁶-methyladenosine; ccRCC: clear cell renal cell carcinoma; DEGs: differentially expressed genes; NUF2: NUF2 component of NDC80 kinetochore complex; CDCA3: cell division cycle associated 3; CKAP2L: cytoskeleton associated protein 2 like; KIF14: kinesin family member 14; ASPM: assembly factor for spindle microtubules; METTL14: methyltransferase 14; OS: overall survival; FPKM: fragments per kilobase million; GEO: gene expression omnibus; TCGA: the Cancer Genome Atlas; RMA: robust multi-array average expression measure; WGCNA: weighted gene co-expression network analysis; GO: gene ontology; KEGG: kyoto encyclopedia of genes and genomes; ROC: receiver operating characteristic curve; AUC: area under the curve; RIP: RNA immunoprecipitation; qPCR: real-time quantitative PCR.

Transcriptome research identifies four hub genes related to primary myelofibrosis: a holistic research by weighted gene co-expression network analysis

Objectives: This study aimed to identify specific diagnostic as well as predictive targets of primary myelofibrosis (PMF).

Methods: The gene expression profiles of GSE26049 were obtained from Gene Expression Omnibus (GEO) dataset, WGCNA was constructed to identify the most related module of PMF. Subsequently, Gene Ontology (GO), Kyoto Encyclopedia Genes and Genomes (KEGG), Gene Set Enrichment Analysis (GSEA) and Protein-Protein interaction (PPI) network were conducted to fully understand the detailed information of the interested green module. Machine learning, Principal component analysis (PCA), and expression pattern analysis including immunohistochemistry and immunofluorescence of genes and proteins were performed to validate the reliability of these hub genes.

Results: Green module was strongly correlated with PMF disease after WGCNA analysis. 20 genes in green module were identified as hub genes responsible for the progression of PMF. GO, KEGG revealed that these hub genes were primarily enriched in erythrocyte differentiation, transcription factor binding, hemoglobin complex, transcription factor complex and cell cycle, etc. Among them, EPB42, CALR, SLC4A1 and MPL had the most correlations with PMF. Machine learning, Principal component analysis (PCA), and expression pattern analysis proved the results in this study.

Conclusions: EPB42, CALR, SLC4A1 and MPL were significantly highly expressed in PMF samples. These four genes may be considered as candidate prognostic biomarkers and potential therapeutic targets for early stage of PMF. The effects are worth expected whether in the diagnosis at early stage or as therapeutic target.

CENPA promotes clear cell renal cell carcinoma progression and metastasis via Wnt/β-catenin signaling pathway

Clear cell renal cell carcinoma (ccRCC) is the most common malignant tumor of the kidney. New and reliable biomarkers are in urgent need for ccRCC diagnosis and prognosis. The CENP family is overexpressed in many types of cancers, but its functions in ccRCC have not been fully clarified. In this paper, we found that several CENP family members were highly expressed in ccRCC tissues. Also, CENPA expression level was related to clinicopathological grade and prognosis by weighted gene co-expression network analysis (WGCNA). CENPA served as a representative CENP family member as a ccRCC biomarker. Further in vitro experiments verified that overexpression of CENPA promoted ccRCC proliferation and metastasis by accelerating the cell cycle and activating the Wnt/β-catenin signaling pathway. The elevated β-catenin led by CENPA overexpression translocated to nucleus for downstream effect. Functional recovery experiment confirmed that Wnt/β-catenin pathway was essential for ccRCC progression and metastasis. Developing selective drugs targeting CENPA may be a promising direction for cancer treatment.

Development and Validation of a Radiomic Nomogram for Predicting the Prognosis of Kidney Renal Clear Cell Carcinoma

Purpose

The present study aims to comprehensively investigate the prognostic value of a radiomic nomogram that integrates contrast-enhanced computed tomography (CECT) radiomic signature and clinicopathological parameters in kidney renal clear cell carcinoma (KIRC).

Methods

A total of 136 and 78 KIRC patients from the training and validation cohorts were included in the retrospective study. The intraclass correlation coefficient (ICC) was used to assess reproducibility of radiomic feature extraction. Univariate Cox analysis and least absolute shrinkage and selection operator (LASSO) as well as multivariate Cox analysis were utilized to construct radiomic signature and clinical signature in the training cohort. A prognostic nomogram was established containing a radiomic signature and clinicopathological parameters by using a multivariate Cox analysis. The predictive ability of the nomogram [relative operating characteristic curve (ROC), concordance index (C-index), Hosmer–Lemeshow test, and calibration curve] was evaluated in the training cohort and validated in the validation cohort. Patients were split into high- and low-risk groups, and the Kaplan–Meier (KM) method was conducted to identify the forecasting ability of the established models. In addition, genes related with the radiomic risk score were determined by weighted correlation network analysis (WGCNA) and were used to conduct functional analysis.

Results

A total of 2,944 radiomic features were acquired from the tumor volumes of interest (VOIs) of CECT images. The radiomic signature, including ten selected features, and the clinical signature, including three selected clinical variables, showed good performance in the training and validation cohorts [area under the curve (AUC), 0.897 and 0.712 for the radiomic signature; 0.827 and 0.822 for the clinical signature, respectively]. The radiomic prognostic nomogram showed favorable performance and calibration in the training cohort (AUC, 0.896, C-index, 0.846), which was verified in the validation cohort (AUC, 0.768). KM curves indicated that the progression-free interval (PFI) time was dramatically shorter in the high-risk group than in the low-risk group. The functional analysis indicated that radiomic signature was significantly associated with T cell activation.

Conclusions

The nomogram combined with CECT radiomic and clinicopathological signatures exhibits excellent power in predicting the PFI of KIRC patients, which may aid in clinical management and prognostic evaluation of cancer patients.

Loss of Setd2 associates with aberrant microRNA expression and contributes to inflammatory bowel disease progression in mice

Both SETD2-mediated H3K36me3 and miRNAs play critical epigenetic roles in inflammatory bowel disease (IBD) and involve in the dysfunctional intestinal barrier. However, little is known about cross-talk between these two types of regulators in IBD progression. We performed small RNA sequencing of Setd2 epithelium-specific knockout mice (Setd2^Vil-KO) and wild-type controls, both with DSS-induced colitis, and designed a framework for integrative analysis. Firstly, we integrated the downloaded ChIP-seq data with miRNA expression profiles and identified a significant intersection of pre-miRNA expression and H3K36me3 modification. A significant inverse correlation was detected between changes of H3K36me3 modification and expression of the 171 peak-covered miRNAs. We further integrated RNA-seq data with predicted miRNA targets to screen negatively regulated miRNA-mRNA pairs and found the H3K36me3-associated differentially expressed microRNAs significantly enriched in cell-cell junction and signaling pathways. Using network analysis, we identified ten hub miRNAs, among which six are H3K36me3-associated, suggesting therapeutic targets for IBD patients with SETD2-deficiency.

The Choice of Candidates in Survival Markers Based on Coordinated Gene Expression in Renal Cancer

We aimed to identify and investigate genes that are essential for the development of clear cell renal cell carcinoma (ccRCC) and sought to shed light on the mechanisms of its progression and create prognostic markers for the disease. We used real-time PCR to study the expression of 20 genes that were preliminarily selected based on their differential expression in ccRCC, in 68 paired tumor/normal samples. Upon ccRCC progression, seven genes that showed an initial increase in expression showed decreased expression. The genes whose expression levels did not significantly change during progression were associated mainly with metabolic and inflammatory processes. The first group included CA9, NDUFA4L2, EGLN3, BHLHE41, VWF, IGFBP3, and ANGPTL4, whose expression levels were coordinately decreased during tumor progression. This expression coordination and gene function is related to the needs of tumor development at different stages. Specifically, the high correlation coefficient of EGLN3 and NDUFA4L2 expression may indicate the importance of the coordinated regulation of glycolysis and mitochondrial metabolism. A panel of CA9, EGLN3, BHLHE41, and VWF enabled the prediction of survival for more than 3.5 years in patients with ccRCC, with a probability close to 90%. Therefore, a coordinated change in the expression of a gene group during ccRCC progression was detected, and a new panel of markers for individual survival prognosis was identified.

Targeted Luminescent Probes for Precise Upconversion/NIR II Luminescence Diagnosis of Lung Adenocarcinoma

In this research, the antibody of the searched hub genes has been proposed to combine with a rare-earth composite for an upconversion luminescence (UCL) and downconversion (DCL) NIR-II imaging strategy for the diagnosis of lung adenocarcinoma (LUAD). Weighted gene co-expression network analysis is used to search the most relevant hub genes, and the required top genes that contribute to tumorigenesis (negative: CLEC3B, MFAP4, PECAM1, and FHL1; positive: CCNB2, CDCA5, HMMR, and TOP2A) are identified and validated by survival analysis and transcriptional and translational results. Meanwhile, fluorescence imaging probes (NaYF₄:Yb,Er,Eu@NaYF₄:Nd, denoted as NYF:Eu NPs) with multimodal optical imaging properties of downconversion and upconversion luminescence in the visible region and luminescence in the near infrared II region are designed with various uniform sizes and enhanced penetration and sensitivity. Finally, when the NYF:Eu NP probe is combined with antibodies of these chosen positive hub genes (such as, TOP2A and CCNB2), the in vitro and in vivo animal experiments (flow cytometry, cell counting kit-8 assay using A549 cells, and in vivo immunohistochemistry IHC microscopy images of LUAD from patient cases) indicate that the designed nanoprobes can be excellently used as a targeted optical probe for future accurate diagnosis and surgery navigation of LUAD in contrast with other cancer cells and normal cells. This strategy of antibodies combined with optical probes provides a dual-modal luminescence imaging method for precise medicine.

Identifying Clear Cell Renal Cell Carcinoma Coexpression Networks Associated with Opioid Signaling and Survival

While opioids constitute the major component of perioperative analgesic regimens for surgery in general, a variety of evidence points to an association between perioperative opioid exposure and longer term oncologic outcomes. The mechanistic details underlying these effects are not well understood. In this study, we focused on clear cell renal cell carcinoma (ccRCC) and utilized RNA sequencing and outcome data from both The Cancer Genome Atlas, as well as a local patient cohort to identify survival-associated gene coexpression networks. We then projected drug-induced transcriptional profiles from in vitro cancer cells to predict drug effects on these networks and recurrence-free, cancer-specific, and overall survival. The opioid receptor agonist, leu-enkephalin, was predicted to have antisurvival effects in ccRCC, primarily through Th2 immune- and NRF2-dependent macrophage networks. Conversely, the antagonist, naloxone, was predicted to have prosurvival effects, primarily through angiogenesis, fatty acid metabolism, and hemopoesis pathways. Eight coexpression networks associated with survival endpoints in ccRCC were identified, and master regulators of the transition from the normal to disease state were inferred, a number of which are linked to opioid pathways. These results are the first to suggest a mechanism for opioid effects on cancer outcomes through modulation of survival-associated coexpression networks. While we focus on ccRCC, this methodology may be employed to predict opioid effects on other cancer types and to personalize analgesic regimens in patients with cancer for optimal outcomes. SIGNIFICANCE: This study suggests a possible molecular mechanism for opioid effects on cancer outcomes generally, with implications for personalization of analgesic regimens.

Identification of KIF4A and its effect on the progression of lung adenocarcinoma based on the bioinformatics analysis

BACKGROUND

Lung adenocarcinoma (LUAD) is the most frequent histological type of lung cancer, and its incidence has displayed an upward trend in recent years. Nevertheless, little is known regarding effective biomarkers for LUAD.

METHODS

The robust rank aggregation method was used to mine differentially expressed genes (DEGs) from the gene expression omnibus (GEO) datasets. The Search Tool for the Retrieval of Interacting Genes (STRING) database was used to extract hub genes from the protein-protein interaction (PPI) network. The expression of the hub genes was validated using expression profiles from TCGA and Oncomine databases and was verified by real-time quantitative PCR (qRT-PCR). The module and survival analyses of the hub genes were determined using Cytoscape and Kaplan-Meier curves. The function of KIF4A as a hub gene was investigated in LUAD cell lines.

RESULTS

The PPI analysis identified seven DEGs including BIRC5, DLGAP5, CENPF, KIF4A, TOP2A, AURKA, and CCNA2, which were significantly upregulated in Oncomine and TCGA LUAD datasets, and were verified by qRT-PCR in our clinical samples. We determined the overall and disease-free survival analysis of the seven hub genes using GEPIA. We further found that CENPF, DLGAP5, and KIF4A expressions were positively correlated with clinical stage. In LUAD cell lines, proliferation and migration were inhibited and apoptosis was promoted by knocking down KIF4A expression.

CONCLUSION

We have identified new DEGs and functional pathways involved in LUAD. KIF4A, as a hub gene, promoted the progression of LUAD and might represent a potential therapeutic target for molecular cancer therapy.

Construction of a Prognostic Risk Prediction Model for Obesity Combined With Breast Cancer

The improvement in the quality of life is accompanied by an accelerated pace of living and increased work-related pressures. Recent decades has seen an increase in the proportion of obese patients, as well as an increase in the prevalence of breast cancer. More and more evidences prove that obesity may be one of a prognostic impact factor in patients with breast cancer. Obesity presents unique diagnostic and therapeutic challenges in the population of breast cancer patients. Therefore, it is essential to have a better understanding of the relationship between obesity and breast cancer. This study aims to construct a prognostic risk prediction model combining obesity and breast cancer. In this study, we obtained a breast cancer sample dataset from the GEO database containing obesity data [determined by the body mass index (BMI)]. A total of 1174 genes that were differentially expressed between breast cancer samples of patients with and without obesity were screened by the rank-sum test. After weighted gene co-expression network analysis (WGCNA), 791 related genes were further screened. Relying on single-factor COX regression analysis to screen the candidate genes to 30, these 30 genes and another set of TCGA data were intersected to obtain 24 common genes. Finally, lasso regression analysis was performed on 24 genes, and a breast cancer prognostic risk prediction model containing 6 related genes was obtained. The model was also found to be related to the infiltration of immune cells. This study provides a new and accurate prognostic model for predicting the survival of breast cancer patients with obesity.

Non-SMC condensin I complex subunit H (NCAPH), a regulator of cell cycle, predicts poor prognosis in lung adenocarcinoma patients: a study mainly based on TCGA and GEO database

Background

Lung adenocarcinoma (LUAD) is the main sub-type of lung cancer, which is a major disease of human death. However, the role of non-SMC condensin I complex subunit H (NCAPH) in LUAD and its possible upstream regulation microRNAs (miRNAs) remains unclearly.

Methods

In this study, we analyzed the NCAPH mRNA and protein expression in normal and cancer tissues mainly based on Human Protein Atlas (HPA) database, The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. With the help of the Kaplan Meier plotter, we explored the prognosis role in LUAD. Furtherly, the co-expressed genes of NCAPH in LUAD were obtained by using cBioPortal, GEPIA and UALCAN database. Then, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of co-expression genes of NCAPH was conducted by DAVID, while the protein-protein interaction (PPI) network was constructed with STRING and hub genes were identified and visualized by Cytoscape software. We also investigated the miRNAs and chemicals that may downregulated the NCAPH expression.

Results

The results showed that NCAPH expression level was elevated in LUAD tissue compared with normal lung tissue and predicted poor prognosis. GO and KEGG pathway enriched analysis of co-expressed genes suggested that NCAPH may play an important role in cell cycle in LUAD. Nine top hub co-expressed genes were all negatively related to the LUAD prognosis. Lastly, 8 miRNAs and 5 chemicals were identified to have the potential to down-regulate the NCAPH expression.

Conclusions

Our study indicated that NCAPH expression in LUAD is a poor prognostic indicator, which may be the potential therapeutic target in the future.

Influence of gene expression on survival of clear cell renal cell carcinoma

Approximately 10%‐20% of patients with clinically localized clear cell renal cell carcinoma (ccRCC) at time of surgery will subsequently experience metastatic progression. Although considerable progression was seen in the systemic treatment of metastatic ccRCC in last 20 years, once ccRCC spreads beyond the confines of the kidney, 5‐year survival is less than 10%. Therefore, significant clinical advances are urgently needed to improve overall survival and patient care to manage the growing number of patients with localized ccRCC. We comprehensively evaluated expression of 388 candidate genes related with survival of ccRCC by using TCGA RNAseq (n = 515), Total Cancer Care (TCC) expression array data (n = 298), and a well characterized Moffitt RCC cohort (n = 248). We initially evaluated all 388 genes for association with overall survival using TCGA and TCC data. Eighty‐one genes were selected for further analysis and tested on Moffitt RCC cohort using NanoString expression analysis. Expression of nine genes (AURKA, AURKB, BIRC5, CCNE1, MK167, MMP9, PLOD2, SAA1, and TOP2A) was validated as being associated with poor survival. Survival prognostic models showed that expression of the nine genes and clinical factors predicted the survival in ccRCC patients with AUC value: 0.776, 0.821 and 0.873 for TCGA, TCC and Moffitt data set, respectively. Some of these genes have not been previously implicated in ccRCC survival and thus potentially offer insight into novel therapeutic targets. Future studies are warranted to validate these identified genes, determine their biological mechanisms and evaluate their therapeutic potential in preclinical studies.

Co-expression Network Analysis Identifies Fourteen Hub Genes Associated with Prognosis in Clear Cell Renal Cell Carcinoma

Renal cancer is a common genitourinary malignance, of which clear cell renal cell carcinoma (ccRCC) has high aggressiveness and leads to most cancer-related deaths. Identification of sensitive and reliable biomarkers for predicting tumorigenesis and progression has great significance in guiding the diagnosis and treatment of ccRCC. Here, we identified 2397 common differentially expressed genes (DEGs) using paired normal and tumor ccRCC tissues from GSE53757 and The Cancer Genome Atlas (TCGA). Then, we performed weighted gene co-expression network analysis and protein-protein interaction network analysis, 17 candidate hub genes were identified. These candidate hub genes were further validated in GSE36895 and Oncomine database and 14 real hub genes were identified. All the hub genes were up-regulated and significantly positively correlated with pathological stage and histologic grade of ccRCC. Survival analysis showed that the higher expression level of each hub gene tended to predict a worse clinical outcome. ROC analysis showed that all the hub genes can accurately distinguish between tumor and normal samples, and between early stage and advanced stage ccRCC. Moreover, all the hub genes were positively associated with distant metastasis, lymph node infiltration, tumor recurrence and the expression of MKi67, suggesting these genes might promote tumor proliferation, invasion and metastasis. Furthermore, the functional annotation demonstrated that most genes were enriched in cell-cycle related biological function. In summary, our study identified 14 potential biomarkers for predicting tumorigenesis and progression, which might contribute to early diagnosis, prognosis prediction and therapeutic intervention.

Weighted correlation network bioinformatics uncovers a key molecular biosignature driving the left-sided heart failure

BACKGROUND

Left-sided heart failure (HF) is documented as a key prognostic factor in HF. However, the relative molecular mechanisms underlying left-sided HF is unknown. The purpose of this study is to unearth significant modules, pivotal genes and candidate regulatory components governing the progression of left-sided HF by bioinformatical analysis.

METHODS

A total of 319 samples in GSE57345 dataset were used for weighted gene correlation network analysis (WGCNA). ClusterProfiler package in R was used to conduct functional enrichment for genes uncovered from the modules of interest. Regulatory networks of genes were built using Cytoscape while Enrichr database was used for identification of transcription factors (TFs). The MCODE plugin was used for identifying hub genes in the modules of interest and their validation was performed based on GSE1869 dataset.

RESULTS

A total of six significant modules were identified. Notably, the blue module was confirmed as the most crucially associated with left-sided HF, ischemic heart disease (ISCH) and dilated cardiomyopathy (CMP). Functional enrichment conveyed that genes belonging to this module were mainly those driving the extracellular matrix-associated processes such as extracellular matrix structural constituent and collagen binding. A total of seven transcriptional factors, including Suppressor of Zeste 12 Protein Homolog (SUZ12) and nuclear factor erythroid 2 like 2 (NFE2L2), adrenergic receptor (AR), were identified as possible regulators of coexpression genes identified in the blue module. A total of three key genes (OGN, HTRA1 and MXRA5) were retained after validation of their prognostic value in left-sided HF. The results of functional enrichment confirmed that these key genes were primarily involved in response to transforming growth factor beta and extracellular matrix.

CONCLUSION

We uncovered a candidate gene signature correlated with HF, ISCH and CMP in the left ventricle, which may help provide better prognosis and therapeutic decisions and in HF, ISCH and CMP patients.

CCNB2, TOP2A, and ASPM Reflect the Prognosis of Hepatocellular Carcinoma, as Determined by Weighted Gene Coexpression Network Analysis

Background

Hepatocellular carcinoma (HCC) is characterized by increased mortality and poor prognosis. We aimed to identify potential prognostic markers by weighted gene coexpression network analysis (WGCNA), to assist clinical outcome prediction and improve treatment decisions for HCC patients.

Methods

Prognosis-related gene modules were first established by WGCNA. Venn diagrams obtained intersection genes of module genes and differentially expressed genes. The Kaplan-Meier overall survival curves and disease-free survival curves of intersection genes were further analyzed on the Gene Expression Profiling Interactive Analysis website. Chi-square tests were performed to explore the associations between prognostic gene expressions and clinicopathological features.

Results

CCNB2, TOP2A, and ASPM were identified as both prognosis-related genes and differentially expressed genes. TOP2A (HR: 1.7, P = 0.003) and ASPM (HR: 1.8, P < 0.001) exhibited a significant difference between the high- and low-expression groups in the overall survival analysis, while CCNB2 (HR: 1.4, P = 0.052) was not statistically significant. CCNB2 (HR: 1.5, P = 0.006), TOP2A (HR: 1.7, P < 0.001), and ASPM (HR: 1.6, P = 0.003) were all statistically significant in the disease-free survival analysis. All three genes were significantly associated with race and fetoprotein values (P < 0.05). CCNB2 expression was associated with tumor stage (P = 0.01), and ASPM expression was associated with new tumor events (P = 0.03).

Conclusion

Overexpression of CCNB2, TOP2A, and ASPM are associated with poor prognosis, and these genes could serve as potential prognostic markers and therapeutic targets for HCC.

Identification of 13 Key Genes Correlated With Progression and Prognosis in Hepatocellular Carcinoma by Weighted Gene Co-expression Network Analysis

Hepatocellular carcinoma (HCC) remains hard to diagnose early and cure due to a lack of accurate biomarkers and effective treatments. Hence, it is necessary to explore the tumorigenesis and tumor progression of HCC to discover new biomarkers for clinical treatment. We performed weighted gene co-expression network analysis (WGCNA) to explore hub genes that have high correlation with clinical information. In this study, we found 13 hub genes (GTSE1, PLK1, NCAPH, SKA3, LMNB2, SPC25, HJURP, DEPDC1B, CDCA4, UBE2C, LMNB1, PRR11, and SNRPD2) that have high correlation with histologic grade in HCC by analyzing TCGA LIHC dataset. All of these 13 hub genes could be used to effectively distinguish high histologic grade from low histologic grade of HCC through analysis of the ROC curve. The overall survival and disease-free survival information showed that high expression of these 13 hub genes led to poor prognosis. Meanwhile, these 13 hub genes had significantly different expression in HCC tumor and non-tumor tissues. We downloaded GSE6764, which contains corresponding clinical information, to validate the expression of these 13 hub genes. At the same time, we performed quantitative real-time PCR to validate the differences in the expression tendencies of these 13 hub genes between HCC tumor tissues and non-tumor tissues and high histologic grade and low histologic grade. We also explored mutation and methylation information of these 13 hub genes for further study. In summary, 13 hub genes correlated with the progression and prognosis of HCC were discovered by WGCNA in our study, and these hub genes may contribute to the tumorigenesis and tumor progression of HCC.

Identification and validation of key modules and hub genes associated with the pathological stage of oral squamous cell carcinoma by weighted gene co-expression network analysis

Background

Oral squamous cell carcinoma (OSCC) is a major lethal malignant cancer of the head and neck region, yet its molecular mechanisms of tumourigenesis are still unclear.

Patients and methods

We performed weighted gene co-expression network analysis (WGCNA) on RNA-sequencing data with clinical information obtained from The Cancer Genome Atlas (TCGA) database. The relationship between co-expression modules and clinical traits was investigated by Pearson correlation analysis. Furthermore, the prognostic value and expression level of the hub genes of these modules were validated based on data from the TCGA database and other independent datasets from the Gene Expression Omnibus (GEO) database and the Human Protein Atlas database. The significant modules and hub genes were also assessed by functional analysis and gene set enrichment analysis (GSEA).

Results

We found that the turquoise module was strongly correlated with pathologic T stage and significantly enriched in critical functions and pathways related to tumourigenesis. PPP1R12B, CFD, CRYAB, FAM189A2 and ANGPTL1 were identified and statistically validated as hub genes in the turquoise module and were closely implicated in the prognosis of OSCC. GSEA indicated that five hub genes were significantly involved in many well-known cancer-related biological functions and signaling pathways.

Conclusion

In brief, we systematically discovered a co-expressed turquoise module and five hub genes associated with the pathologic T stage for the first time, which provided further insight that WGCNA may reveal the molecular regulatory mechanism involved in the carcinogenesis and progression of OSCC. In addition, the five hub genes may be considered candidate prognostic biomarkers and potential therapeutic targets for the precise early diagnosis, clinical treatment and prognosis of OSCC in the future.

Identifying a ten-microRNA signature as a superior prognosis biomarker in colon adenocarcinoma

Background

Increasing studies have suggested that aberrant expression of microRNAs might play essential roles in the progression of cancers. In this study, we sought to construct a high-specific and superior microRNAs signature to improve the survival prediction of colon adenocarcinoma (COAD) patients.

Methods

The genome-wide miRNAs, mRNA and lncRNA expression profiles and corresponding clinical information of COAD were collected from the TCGA database. Differential expression analysis, Kaplan–Meier curve and time-dependent ROC curve were calculated and performed using R software and GraphPad Prism7. Univariate and multivariate Cox analysis was performed to evaluate the prognostic ability of signature. Functional enrichment analysis was analyzed using STRING database.

Results

We identified ten prognosis-related microRNAs, including seven risky factors (hsa-miR-197, hsa-miR-32, hsa-miR-887, hsa-miR-3199-2, hsa-miR-4999, hsa-miR-561, hsa-miR-210) and three protective factors (hsa-miR-3917, hsa-miR-3189, hsa-miR-6854). The Kaplan–Meier survival analysis showed that the patients with high risk score had shorter overall survival (OS) in test series. And the similar results were observed in both validation and entire series. The time-dependent ROC curve suggested this signature have high accuracy of OS for COAD. The Multivariate Cox regression analysis and stratification analysis suggested that the ten-microRNA signature was an independent factor after being adjusted with other clinical characteristics. In addition, we also found microRNA signature have higher AUC than other signature. Furthermore, we identified some miRNA-target genes that affect lymphatic metastasis and invasion of COAD patients.

Conclusion

In this study, we established a ten-microRNA signature as a potentially reliable and independent biomarker for survival prediction of COAD patients.

Bioinformatic analysis identifies potentially key differentially expressed genes in oncogenesis and progression of clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is one of the most common and lethal types of cancer within the urinary system. Great efforts have been made to elucidate the pathogeny. However, the molecular mechanism of ccRCC is still not well understood. The aim of this study is to identify key genes in the carcinogenesis and progression of ccRCC. The mRNA microarray dataset GSE53757 was downloaded from the Gene Expression Omnibus database. The GSE53757 dataset contains tumor and matched paracancerous specimens from 72 ccRCC patients with clinical stage I to IV. The linear model of microarray data (limma) package in R language was used to identify differentially expressed genes (DEGs). The protein–protein interaction (PPI) network of the DEGs was constructed using the search tool for the retrieval of interacting genes (STRING). Subsequently, we visualized molecular interaction networks by Cytoscape software and analyzed modules with MCODE. A total of 1,284, 1,416, 1,610 and 1,185 up-regulated genes, and 932, 1,236, 1,006 and 929 down-regulated genes were identified from clinical stage I to IV ccRCC patients, respectively. The overlapping DEGs among the four clinical stages contain 870 up-regulated and 645 down-regulated genes. The enrichment analysis of DEGs in the top module was carried out with DAVID. The results showed the DEGs of the top module were mainly enriched in microtubule-based movement, mitotic cytokinesis and mitotic chromosome condensation. Eleven up-regulated genes and one down-regulated gene were identified as hub genes. Survival analysis showed the high expression of CENPE, KIF20A, KIF4A, MELK, NCAPG, NDC80, NUF2, TOP2A, TPX2 and UBE2C, and low expression of ACADM gene could be involved in the carcinogenesis, invasion or recurrence of ccRCC. Literature retrieval results showed the hub gene NDC80, CENPE and ACADM might be novel targets for the diagnosis, clinical treatment and prognosis of ccRCC. In conclusion, the findings of present study may help us understand the molecular mechanisms underlying the carcinogenesis and progression of ccRCC, and provide potential diagnostic, therapeutic and prognostic biomarkers.

Bladder Cancer Exhibiting High Immune Infiltration Shows the Lowest Response Rate to Immune Checkpoint Inhibitors

Background: Bladder urothelial cancer (BLCA) treatment using immune checkpoint inhibitors (IMCIs) can result in long-lasting clinical benefits. However, only a fraction of patients respond to such treatment. In this study, we aimed to identify the relationships between immune cell infiltration levels (ICILs) and IMCIs and identify markers for ICILs.

Methods: ICILs were estimated based on single-sample gene set enrichment analysis. The response rates of different ICILs to IMCIs were calculated by combining the ICILs of molecular subtypes in BLCA with the response rates of different molecular subtypes of IMvigor 210 trials to a programmed cell death ligand-1 inhibitor. Weighted gene co-expression network analysis was used to identify modules of interest with ICILs. Functional enrichment analysis was performed to functionally annotate the modules. Screening of key genes and unsupervised clustering were used to identify candidate biomarkers. Tumor IMmune Estimation Resource was used to validate the relationships between the biomarkers and ICILs. Finally, we verified the expression of key genes in molecular subtypes of different response rates for IMCIs.

Findings: The basal squamous subtype and luminal infiltrated subtype, which showed low response rates for IMCIs, had the highest levels of immune infiltration. The neuronal subtypes, which showed the highest response rates to IMCIs, had low ICILs. The modules of interest and key genes were determined based on topological overlap measurement, clustering results, and inclusion criteria. Modules highly correlated with ICILs were mainly enriched in immune responses and epithelial–mesenchymal transition. After screening the key genes in the modules, five candidate biomarkers (CD48, SEPT1, ACAP1, PPP1R16B, and IL16) were selected by unsupervised clustering. The key genes were inversely associated with tumor purity and were mostly expressed in the basal squamous subtype and luminal infiltrated subtypes.

Interpretation: Patients with high ICILs may benefit the least from treatment with IMCIs. Five key genes could predict ICILs in BLCA, and their high expression suggested that the response rate to IMCIs may decrease.

GABRQ expression is a potential prognostic marker for patients with clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the most common type of kidney cancer. Novel biomarkers of ccRCC may provide crucial information on tumor features and prognosis. The present study aimed to determine whether the expression of γ-aminobutyric acid (GABA) A receptor subunit θ (GABRQ) could serve as a novel prognostic marker of ccRCC. GABA is the main inhibitory neurotransmitter in the brain that activates the receptor GABA_A, which is comprised of three subunit isoforms: GABRA3, GABRB3 and GABRQ. A recent study reported that GABRQ is involved in the initiation and progression of hepatocellular carcinoma; however, the role of GABRQ in ccRCC remains unknown. In the present study, clinical and transcriptomic data were obtained from cohorts of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). Differential GABRQ expression levels among early (TI and II), late (TIII and IV), nonmetastatic (M0) and metastatic (M1, primary tumor) stages of ccRCC samples were then identified. Furthermore, the use of GABRQ as a prognostic gene was analyzed using Uno's C-index based on the time-dependent area under the curve (AUC), the AUC of the receiver operating characteristic curve at 5 years, the Kaplan-Meier survival curve and multivariate analysis. The survival curve analysis revealed that low GABRQ mRNA expression was significantly associated with a poor prognosis of ccRCC (P<0.001 and P=0.0012 for TCGA and ICGC data, respectively). In addition, analyses of the C-index and AUC values further supported this discriminatory power. Furthermore, the prognostic value of GABRQ mRNA expression was confirmed by multivariate Cox regression analysis. Taken together, these results suggested that GABRQ mRNA expression may be considered as a novel prognostic biomarker of ccRCC.

Network-based analysis reveals novel gene signatures in the peripheral blood of patients with sporadic nonsyndromic thoracic aortic aneurysm

Thoracic aortic aneurysm (TAA), a serious cardiovascular disease that causes morbidity and mortality worldwide. At present, few biomarkers can accurately diagnose the appearance of TAA before dissection or rupture. Our research has the intention to investigate the developing applicable biomarkers for TAA promising clinically diagnostic biomarkers or probable regulatory targets for TAA. In our research, we built correlation networks utilizing the expression profile of peripheral blood mononuclear cell obtained from a public microarray data set (GSE9106). Furthermore, we chose the turquoise module, which has the strongest significance with TAA and was further analyzed. Fourteen genes that overlapped with differentially expressed proteins in the medial aortic layer were obtained. Subsequently, we verified the results applying quantitative polymerase chain reaction (Q-PCR) to our clinical specimen. In general, the Q-PCR results coincide with the majority of the expression profile. Fascinatingly, a notable change occurred in CLU, DES, MYH10, and FBLN5. In summary, using weighted gene coexpression analysis, our study indicates that CLU, DES, MYH10, and FBLN5 were identified and validated to be related to TAA and might be candidate biomarkers or therapeutic targets for TAA.

Prognostic role of the beta-2 adrenergic receptor in clear cell renal cell carcinoma

The beta-2 adrenergic receptor (ADRB2) regulates the proliferation, apoptosis, angiogenesis, migration, and metastasis of cancer cells. However, its function in the progression of clear cell renal cell carcinoma (ccRCC) is unknown. Here, we report that ADRB2 can be a novel prognostic factor for patients with ccRCC. The differential expression of ADRB2 in low-stage (stages I and II), high-stage (stages III and IV), low-grade (grades I and II), and high-grade (grades III and IV) ccRCC was identified in cohorts of patients from The Cancer Genome Atlas and the International Cancer Genome Consortium. We evaluated ADRB2 expression as a prognostic factor using the Kaplan-Meier survival curve, multivariate analysis, time-dependent area under the curve (AUC) of Uno’s C-index, and AUC of the receiver operating characteristics (ROC) at five years. Kaplan-Meier analysis revealed that reduced ADRB2 expression is associated with poor prognosis in ccRCC patients. Analysis of C-indices and AUC-ROC further confirmed this result. Moreover, multivariate analysis confirmed the prognostic significance of ADRB2 expression. Collectively, these findings suggest that ADRB2 is a potential prognostic factor for ccRCC.

Prognostic Role of TMED3 in Clear Cell Renal Cell Carcinoma: A Retrospective Multi-Cohort Analysis

Transmembrane p24 trafficking protein 3 (TMED3) is a metastatic suppressor in colon cancer and hepatocellular carcinoma. However, its function in the progression of clear cell renal cell carcinoma (ccRCC) is unknown. Here, we report that TMED3 could be a new prognostic marker for ccRCC. Patient data were extracted from cohorts in the Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC). Differential expression of TMED3 was observed between the low stage (Stage I and II) and high stage (Stage III and IV) patients in the TCGA and ICGC cohorts and between the low grade (Grade I and II) and high grade (Grade III and IV) patients in the TCGA cohort. Further, we evaluated TMED3 expression as a prognostic gene using Kaplan-Meier survival analysis, multivariate analysis, the time-dependent area under the curve (AUC) of Uno’s C-index, and the AUC of the receiver operating characteristics at 5 years. The Kaplan-Meier analysis revealed that TMED3 overexpression was associated with poor prognosis for ccRCC patients. Analysis of the C-indices and area under the receiver operating characteristic curve further supported this. Multivariate analysis confirmed the prognostic significance of TMED3 expression levels (P = 0.005 and 0.006 for TCGA and ICGC, respectively). Taken together, these findings demonstrate that TMED3 is a potential prognostic factor for ccRCC.