Gene network inherent in genomic big data improves the accuracy of prognostic prediction for cancer patients

PKI: A bioinformatics method of quantifying the importance of nodes in gene regulatory network via a pseudo knockout index

BACKGROUND

Gene regulatory network (GRN) is a model that characterizes the complex relationships between genes and thereby provides an informatics environment to measure the importance of nodes. The evaluation of important nodes in a GRN can effectively refer to their functional implications severing as key players in particular biological processes, such as master regulator and driver gene. Currently, it is mainly based on network topological parameters and focuses only on evaluating a single node individually. However, genes and products play their functions by interacting with each other. It is worth noting that the effects of gene combinations in GRN are not simply additive. Key combinations discovery is of significance in revealing gene sets with important functions. Recently, with the development of single-cell RNA-sequencing (scRNA-seq) technology, we can quantify gene expression profiles of individual cells that provide the potential to identify crucial nodes in gene regulations regarding specific condition, e.g., stem cell differentiation.

RESULTS

In this paper, we propose a bioinformatics method, called Pseudo Knockout Importance (PKI), to quantify the importance of node and node sets in a specific GRN structure using time-course scRNA-seq data. First, we construct ordinary differential equations to approach the gene regulations during cell differentiation. Then we design gene pseudo knockout experiments and define PKI score evaluation criteria based on the coefficient of determination. The importance of nodes can be described as the influence on the ODE system of removing variables. For key gene combinations, PKI is derived as a combinatorial optimization problem of quantifying the in silico gene knockout effects.

CONCLUSIONS

Here, we focus our analyses on the specific GRN of embryonic stem cells with time series gene expression profile. To verify the effectiveness and advantage of PKI method, we compare its node importance rankings with other twelve kinds of centrality-based methods, such as degree and Latora closeness. For key node combinations, we compare the results with the method based on minimum dominant set. Moreover, the famous combinations of transcription factors in induced pluripotent stem cell are also employed to verify the vital gene combinations identified by PKI. These results demonstrate the reliability and superiority of the proposed method.

Detecting prognostic biomarkers of breast cancer by regularized Cox proportional hazards models

BACKGROUND

The successful identification of breast cancer (BRCA) prognostic biomarkers is essential for the strategic interference of BRCA patients. Recently, various methods have been proposed for exploring a small prognostic gene set that can distinguish the high-risk group from the low-risk group.

METHODS

Regularized Cox proportional hazards (RCPH) models were proposed to discover prognostic biomarkers of BRCA from gene expression data. Firstly, the maximum connected network with 1142 genes by mapping 956 differentially expressed genes (DEGs) and 677 previously BRCA-related genes into the gene regulatory network (GRN) was constructed. Then, the 72 union genes of the four feature gene sets identified by Lasso-RCPH, Enet-RCPH, [Formula: see text]-RCPH and SCAD-RCPH models were recognized as the robust prognostic biomarkers. These biomarkers were validated by literature checks, BRCA-specific GRN and functional enrichment analysis. Finally, an index of prognostic risk score (PRS) for BRCA was established based on univariate and multivariate Cox regression analysis. Survival analysis was performed to investigate the PRS on 1080 BRCA patients from the internal validation. Particularly, the nomogram was constructed to express the relationship between PRS and other clinical information on the discovery dataset. The PRS was also verified on 1848 BRCA patients of ten external validation datasets or collected cohorts.

RESULTS

The nomogram highlighted that the importance of PRS in guiding significance for the prognosis of BRCA patients. In addition, the PRS of 301 normal samples and 306 tumor samples from five independent datasets showed that it is significantly higher in tumors than in normal tissues ([Formula: see text]). The protein expression profiles of the three genes, i.e., ADRB1, SAV1 and TSPAN14, involved in the PRS model demonstrated that the latter two genes are more strongly stained in tumor specimens. More importantly, external validation illustrated that the high-risk group has worse survival than the low-risk group ([Formula: see text]) in both internal and external validations.

CONCLUSIONS

The proposed pipelines of detecting and validating prognostic biomarker genes for BRCA are effective and efficient. Moreover, the proposed PRS is very promising as an important indicator for judging the prognosis of BRCA patients.

Novel Insights into the Prognosis and Immunological Value of the SLC35A (Solute Carrier 35A) Family Genes in Human Breast Cancer

According to statistics 2020, female breast cancer (BRCA) became the most commonly diagnosed malignancy worldwide. Prognosis of BRCA patients is still poor, especially in population with advanced or metastatic. Particular functions of each members of the solute carrier 35A (SLC35A) gene family in human BRCA are still unknown regardless of awareness that they play critical roles in tumorigenesis and progression. Using integrated bioinformatics analyses to identify therapeutic targets for specific cancers based on transcriptomics, proteomics, and high-throughput sequencing, we obtained new information and a better understanding of potential underlying molecular mechanisms. Leveraging BRCA dataset that belongs to The Cancer Genome Atlas (TCGA), which were employed to clarify SLC35A gene expression levels. Then we used a bioinformatics approach to investigate biological processes connected to SLC35A family genes in BRCA development. Beside that, the Kaplan-Meier estimator was leveraged to explore predictive values of SLC35A family genes in BCRA patients. Among individuals of this family gene, expression levels of SLC35A2 were substantially related to poor prognostic values, result from a hazard ratio of 1.3 (with 95 percent confidence interval (95% CI: 1.18-1.44), the p for trend (ptrend) is 3.1 × 10-7). Furthermore, a functional enrichment analysis showed that SLC35A2 was correlated with hypoxia-inducible factor 1A (HIF1A), heat shock protein (HSP), E2 transcription factor (E2F), DNA damage, and cell cycle-related signaling. Infiltration levels observed in specific types of immune cell, especially the cluster of differentiation found on macrophages and neutrophils, were positively linked with SLC35A2 expression in multiple BRCA subclasses (luminal A, luminal B, basal, and human epidermal growth factor receptor 2). Collectively, SLC35A2 expression was associated with a lower recurrence-free survival rate, suggesting that it could be used as a biomarker in treating BRCA.

Identification and Complete Validation of Prognostic Gene Signatures for Human Papillomavirus-Associated Cancers: Integrated Approach Covering Different Anatomical Locations

Human papillomavirus (HPV) infects squamous epithelium and is a major cause of cervical cancer (CC) and a subset of head and neck cancers (HNC). Virus-induced tumorigenesis, molecular alterations, and related prognostic markers are expected to be similar between the two cancers, but they remain poorly understood. We present integrated molecular analysis of HPV-associated tumors from TCGA and GEO databases and identify prognostic biomarkers. Analysis of gene expression profiles identified common upregulated genes and pathways of DNA replication and repair in the HPV-associated tumors. We established 34 prognostic gene signatures with a universal cutoff value in TCGA-CC using Elastic Net Cox regression analysis. We were able to externally validate our results in the TCGA-HNC and several GEO data sets, and demonstrated prognostic power in HPV-associated HNC, but not in HPV-negative cancers. The HPV-related prognostic and predictive indicator did not discriminate other cancers, except bladder urothelial carcinoma. These results identify and completely validate a highly selective prognostic system and its cross-usefulness in HPV-associated cancers, regardless of the tumor's anatomical subsite.IMPORTANCE Persistent infection with high-risk HPV interferes with cell function regulation and causes cell mutations, which accumulate over the long term and eventually develop into cancer. Results of pathway enrichment analysis presumably showed this accumulation of intracellular damage during the chronic HPV-infected state. We used highly advanced statistical methods to identify the most appropriate genes and coefficients and developed the HPV-related prognostic and predictive indicator (HPPI) risk scoring system. We applied the same cutoff value to training and validation sets and demonstrated good prognostic performance in both data sets, and confirmed a consistent trend in external validation. Moreover, HPPI presented significant validation results for bladder cancer suspected to be related to HPV. This suggested that our risk scoring system based on the prognostic gene signature could play an important role in the development of treatment strategies for patients with HPV-related cancer.

Pathway-Based Drug-Repurposing Schemes in Cancer: The Role of Translational Bioinformatics

Cancer is a set of complex pathologies that has been recognized as a major public health problem worldwide for decades. A myriad of therapeutic strategies is indeed available. However, the wide variability in tumor physiology, response to therapy, added to multi-drug resistance poses enormous challenges in clinical oncology. The last years have witnessed a fast-paced development of novel experimental and translational approaches to therapeutics, that supplemented with computational and theoretical advances are opening promising avenues to cope with cancer defiances. At the core of these advances, there is a strong conceptual shift from gene-centric emphasis on driver mutations in specific oncogenes and tumor suppressors-let us call that the silver bullet approach to cancer therapeutics-to a systemic, semi-mechanistic approach based on pathway perturbations and global molecular and physiological regulatory patterns-we will call this the shrapnel approach. The silver bullet approach is still the best one to follow when clonal mutations in driver genes are present in the patient, and when there are targeted therapies to tackle those. Unfortunately, due to the heterogeneous nature of tumors this is not the common case. The wide molecular variability in the mutational level often is reduced to a much smaller set of pathway-based dysfunctions as evidenced by the well-known hallmarks of cancer. In such cases "shrapnel gunshots" may become more effective than "silver bullets". Here, we will briefly present both approaches and will abound on the discussion on the state of the art of pathway-based therapeutic designs from a translational bioinformatics and computational oncology perspective. Further development of these approaches depends on building collaborative, multidisciplinary teams to resort to the expertise of clinical oncologists, oncological surgeons, and molecular oncologists, but also of cancer cell biologists and pharmacologists, as well as bioinformaticians, computational biologists and data scientists. These teams will be capable of engaging on a cycle of analyzing high-throughput experiments, mining databases, researching on clinical data, validating the findings, and improving clinical outcomes for the benefits of the oncological patients.

A User-Friendly, Web-Based Integrative Tool (ESurv) for Survival Analysis: Development and Validation Study

BACKGROUND

Prognostic genes or gene signatures have been widely used to predict patient survival and aid in making decisions pertaining to therapeutic actions. Although some web-based survival analysis tools have been developed, they have several limitations.

OBJECTIVE

Taking these limitations into account, we developed ESurv (Easy, Effective, and Excellent Survival analysis tool), a web-based tool that can perform advanced survival analyses using user-derived data or data from The Cancer Genome Atlas (TCGA). Users can conduct univariate analyses and grouped variable selections using multiomics data from TCGA.

METHODS

We used R to code survival analyses based on multiomics data from TCGA. To perform these analyses, we excluded patients and genes that had insufficient information. Clinical variables were classified as 0 and 1 when there were two categories (for example, chemotherapy: no or yes), and dummy variables were used where features had 3 or more outcomes (for example, with respect to laterality: right, left, or bilateral).

RESULTS

Through univariate analyses, ESurv can identify the prognostic significance for single genes using the survival curve (median or optimal cutoff), area under the curve (AUC) with C statistics, and receiver operating characteristics (ROC). Users can obtain prognostic variable signatures based on multiomics data from clinical variables or grouped variable selections (lasso, elastic net regularization, and network-regularized high-dimensional Cox-regression) and select the same outputs as above. In addition, users can create custom gene signatures for specific cancers using various genes of interest. One of the most important functions of ESurv is that users can perform all survival analyses using their own data.

CONCLUSIONS

Using advanced statistical techniques suitable for high-dimensional data, including genetic data, and integrated survival analysis, ESurv overcomes the limitations of previous web-based tools and will help biomedical researchers easily perform complex survival analyses.

FAM213A is linked to prognostic significance in acute myeloid leukemia through regulation of oxidative stress and myelopoiesis

Accurate prediction of malignancies is important in choosing therapeutic strategies. Although there are many genetic and cytogenetic prognostic factors for acute myeloid leukemia (AML), prognosis is difficult to predict because of the heterogeneity of AML. Prognostic factors, including messenger RNA (mRNA) expression, have been determined for other malignancies, but not for AML. A total of 402 patients from The Cancer Genome Atlas, GSE12417 (GPL96, 97), and GSE12417 (GPL570) were included in this study. In Kaplan-Meier curve analyses, high expression of family with sequence similarity 213 member A (FAM213A), which activates antioxidant proteins, was associated with worse prognosis of AML. Similar to the results of the survival curve, C-indices and area under the curve values were high. Current prognostic factors of AML, unlike those of other cancers, do not take mRNA expression into consideration. Thus, the development of mRNA-based prognostic factors would be beneficial for accurate prediction of the survival of AML patients. Additionally, in vivo validation using zebrafish revealed that fam213a is important for myelopoiesis at the developmental stage and is a negative regulator of the p53 tumor suppressor gene. The findings implicate fam213a as a novel prognostic factor for AML patients.

Prognostic Role of Zinc Finger Homeobox 4 in Ovarian Serous Cystadenocarcinoma

Introduction: The zinc finger homeobox 4 (ZFHX4) protein is a crucial molecular regulator of tumor-initiating stem cell-like functions. Objective: This study aimed to determine the role of ZFHX4 in the progression of ovarian serous cystadenocarcinoma (OSC). Methods: Differential gene expression ZFHX4 among 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) OSC patients was identified using four independent cohorts from the Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO). We compared ZFHX4 expression as a prognostic factor using Kaplan-Meier survival curves, multivariate analysis, the time-dependent area under the curve (AUC) of Uno's C-index, and the AUC of the receiver operating characteristics at 4 years post diagnosis. Results: ZFHX4 gene expression in high-stage tumors is significantly higher than in low-stage tumors (TCGA, p = 0.007; GSE9891, p = 0.001). A Kaplan-Meier analysis revealed that elevated expression of ZFHX4 was associated with a poor prognosis in OSC patients for all cohorts, regardless of stage and grade (TCGA, p = 1e-04; GSE9891, p = 0.0044; GSE13876, p = 0.00078; GSE26712, p = 0.039). Analysis of C-indices and the area under the receiver operating characteristic curve further supported this result (C-index: TCGA, 0.599; GSE9891, 0.642; GSE13876, 0.585; GSE26712, 0.597). Moreover, univariate and multivariate Cox hazards analyses confirmed the prognostic significance of ZFHX4 levels. Conclusion: Collectively, these findings suggest that ZFHX4 is a prognostic factor for OSC.

Incorporating genetic networks into case-control association studies with high-dimensional DNA methylation data

Background

In human genetic association studies with high-dimensional gene expression data, it has been well known that statistical selection methods utilizing prior biological network knowledge such as genetic pathways and signaling pathways can outperform other methods that ignore genetic network structures in terms of true positive selection. In recent epigenetic research on case-control association studies, relatively many statistical methods have been proposed to identify cancer-related CpG sites and their corresponding genes from high-dimensional DNA methylation array data. However, most of existing methods are not designed to utilize genetic network information although methylation levels between linked genes in the genetic networks tend to be highly correlated with each other.

Results

We propose new approach that combines data dimension reduction techniques with network-based regularization to identify outcome-related genes for analysis of high-dimensional DNA methylation data. In simulation studies, we demonstrated that the proposed approach overwhelms other statistical methods that do not utilize genetic network information in terms of true positive selection. We also applied it to the 450K DNA methylation array data of the four breast invasive carcinoma cancer subtypes from The Cancer Genome Atlas (TCGA) project.

Conclusions

The proposed variable selection approach can utilize prior biological network information for analysis of high-dimensional DNA methylation array data. It first captures gene level signals from multiple CpG sites using data a dimension reduction technique and then performs network-based regularization based on biological network graph information. It can select potentially cancer-related genes and genetic pathways that were missed by the existing methods.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-3040-x) contains supplementary material, which is available to authorized users.

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.

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 scoring system for osteosarcoma using network-regularized high-dimensional Cox-regression analysis and potential therapeutic targets

With the recent emphasis on the importance of personalized genomic medicine, studies have performed prognostic stratification using gene signatures in cancers. However, these studies have not considered gene networks with clinical data. Therefore, this study aimed to develop a novel prognostic score using grouped variable selection for patients with osteosarcoma. We assessed messenger RNA (mRNA) expression and clinical data from Gene Expression Omnibus to develop a novel prognostic scoring system for patients with osteosarcoma. Variable selection using Network-Regularized high-dimensional Cox-regression analysis with information regarding gene networks obtained from six large pathway databases was performed. We determined the risk score on the linear combination of regression coefficients and mRNA expression values. Log-rank test, UNO's c-index, and area under the curve (AUC) values were determined to evaluate the discriminatory power between the low- and high-risk groups. A recently reported next-generation Connectivity Map was used to identify future therapeutic targets for osteosarcoma. Our novel model had significantly high discriminatory power in predicting overall survival. An optimal c-index of 0.967 was obtained and time-dependent receiver operating characteristic analysis revealed an acceptable predictive value of AUC between 0.953 and 1.000. Knockdown of BACE2 or ING2 and linifanib treatment may improve the prognosis of patients with osteosarcoma. Herein, this novel prognostic scoring system would not only facilitate a more accurate prediction of patient prognosis, but also contribute to the selection of suitable therapeutic alternatives for osteosarcoma patients.

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.

Prognostic significance of EIF4G1 in patients with pancreatic ductal adenocarcinoma

Background

Advances in genomics have greatly improved the survival rate in cancer patients. However, due to genetic heterogeneity, pancreatic ductal adenocarcinoma (PDAC) is still difficult to diagnose early, and its survival rate is extremely low. Therefore, we identified biomarkers that predict the prognosis of PDAC patients using independent cohort data.

Materials and methods

To develop a novel prognostic biomarker, we used the gene expression and clinical data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). Kaplan-Meier survival curve using median values of genes as cutoff showed that EIF4G1 was the only statistically significant gene in the 3 cohorts. We analyzed the prognostic significance of EIF4G1 using the time-dependent area under the curve (AUC) of Uno's C-index, the AUC value of the receiver operating characteristics (ROC) at 3 years, and multivariate Cox analysis. We also compared EIF4G1 levels between tumors and matched non-tumor tissues.

Results

EIF4G1 is the only prognostic gene in patients with PDAC, which was selected by Kaplan-Meier survival analysis. The survival curve showed that high expression of EIF4G1 was associated with poor prognosis of PDAC with a good discriminative ability in 3 independent cohorts. The risk stratifying ability of EIF4G1 was demonstrated by analyzing C-indices and AUC values. Multivariate Cox regression confirmed its prognostic significance. EIF4G1 expression was significantly higher in PDAC tissues than in the matched normal tissues.

Conclusion

EIF4G1 could be used as a novel prognostic marker for PDAC and to determine suitable treatment options.

Development of a risk scoring system for patients with papillary thyroid cancer

As the importance of personalized therapeutics in aggressive papillary thyroid cancer (PTC) increases, accurate risk stratification is required. To develop a novel prognostic scoring system for patients with PTC (n = 455), we used mRNA expression and clinical data from The Cancer Genome Atlas. We performed variable selection using Network-Regularized high-dimensional Cox-regression with gene network from pathway databases. The risk score was calculated using a linear combination of regression coefficients and mRNA expressions. The risk score and clinical variables were assessed by several survival analyses. The risk score showed high discriminatory power for the prediction of event-free survival as well as the presence of metastasis. In multivariate analysis, the risk score and presence of metastasis were significant risk factors among the clinical variables that were examined together. In the current study, we developed a risk scoring system that will help to identify suitable therapeutic options for PTC.

SAC3D1: a novel prognostic marker in hepatocellular carcinoma

Centrosome-associated proteins are recognized as prognostic factors in many cancers because centrosomes are critical structures for the cell cycle progression and genomic stability. SAC3D1, however, is associated with centrosome abnormality, although its prognostic potential has not been evaluated in hepatocellular carcinoma (HCC). In this study, 3 independent cohorts (GSE10186, n = 80; TCGA, n = 330 and ICGC, n = 237) were used to assess SAC3D1 as a biomarker, which demonstrated SAC3D1 overexpression in HCC tissues when compared to the matched normal tissues. Kaplan-Meier survival analysis also showed that its overexpression was associated with poor prognosis of HCC with good discriminative ability in 3 independent cohorts (GSE10186, P = 0.00469; TCGA, P = 0.0000413 and ICGC, P = 0.0000114). Analysis of the C-indices and AUC values further supported its discriminative ability. Finally, multivariate analysis confirmed its prognostic significance (GSE10186, P = 0.00695; TCGA, P = 0.0000289 and ICGC, P = 0.0000651). These results suggest a potential of SAC3D1 as a biomarker for HCC.