ce-Subpathway: Identification of ceRNA-mediated subpathways via joint power of ceRNAs and pathway topologies

Identification and Functional Analysis of Individual-Specific Subpathways in Lung Adenocarcinoma

Small molecular networks within complex pathways are defined as subpathways. The identification of patient-specific subpathways can reveal the etiology of cancer and guide the development of personalized therapeutic strategies. The dysfunction of subpathways has been associated with the occurrence and development of cancer. Here, we propose a strategy to identify aberrant subpathways at the individual level by calculating the edge score and using the Gene Set Enrichment Analysis (GSEA) method. This provides a novel approach to subpathway analysis. We applied this method to the expression data of a lung adenocarcinoma (LUAD) dataset from The Cancer Genome Atlas (TCGA) database. We validated the effectiveness of this method in identifying LUAD-relevant subpathways and demonstrated its reliability using an independent Gene Expression Omnibus dataset (GEO). Additionally, survival analysis was applied to illustrate the clinical application value of the genes and edges in subpathways that were associated with the prognosis of patients and cancer immunity, which could be potential biomarkers. With these analyses, we show that our method could help uncover subpathways underlying lung adenocarcinoma.

Illuminating lncRNA Function Through Target Prediction

Most of the transcribed human genome codes for noncoding RNAs (ncRNAs), and long noncoding RNAs (lncRNAs) make for the lion's share of the human ncRNA space. Despite growing interest in lncRNAs, because there are so many of them, and because of their tissue specialization and, often, lower abundance, their catalog remains incomplete and there are multiple ongoing efforts to improve it. Consequently, the number of human lncRNA genes may be lower than 10,000 or higher than 200,000. A key open challenge for lncRNA research, now that so many lncRNA species have been identified, is the characterization of lncRNA function and the interpretation of the roles of genetic and epigenetic alterations at their loci. After all, the most important human genes to catalog and study are those that contribute to important cellular functions-that affect development or cell differentiation and whose dysregulation may play a role in the genesis and progression of human diseases. Multiple efforts have used screens based on RNA-mediated interference (RNAi), antisense oligonucleotide (ASO), and CRISPR screens to identify the consequences of lncRNA dysregulation and predict lncRNA function in select contexts, but these approaches have unresolved scalability and accuracy challenges. Instead-as was the case for better-studied ncRNAs in the past-researchers often focus on characterizing lncRNA interactions and investigating their effects on genes and pathways with known functions. Here, we focus most of our review on computational methods to identify lncRNA interactions and to predict the effects of their alterations and dysregulation on human disease pathways.

Construction and Analysis of the Dysregulated ceRNA Network and Identification of Risk Long Noncoding RNAs in Breast Cancer

Breast cancer (BRCA) is the second leading cause of cancer-related mortality in women worldwide. However, the molecular mechanism involved in the development of BRCA is not fully understood. In this study, based on the miRNA-mediated long non-coding RNA (lncRNA)-protein coding gene (PCG) relationship and lncRNA-PCG co-expression information, we constructed and analyzed a specific dysregulated lncRNA-PCG co-expression network in BRCA. Then, we performed the random walk with restart (RWR) method to prioritize BRCA-related lncRNAs through comparing their RWR score and significance. As a result, we identified 30 risk lncRNAs for BRCA, which can distinguish normal and tumor samples. Moreover, through gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis, we found that these risk lncRNAs mainly synergistically exerted functions related to cell cycle and DNA separation and replication. At last, we developed a four-lncRNA prognostic signature (including AP000851.1, LINC01977, MAFG-DT, SIAH2-AS1) and assessed the survival accuracy of the signature by performing time-dependent receiver operating characteristic (ROC) analysis. The areas under the ROC curve for 1, 3, 5, and 10 years of survival prediction were 0.68, 0.61, 0.62, and 0.63, respectively. The multivariable Cox regression results verified that the four-lncRNA signature could be used as an independent prognostic biomarker in BRCA. In summary, these results have important reference value for the study of diagnosis, treatment, and prognosis evaluation of BRCA.

CCR7 and its related molecules may be potential biomarkers of pulmonary arterial hypertension

Pulmonary arterial hypertension (PAH) is a chronic progressive cardiovascular disease characterized by vascular remodeling and leading to right‐heart failure. The purpose of this research was to further study the pathogenesis of PAH and to detect potential prognostic signatures. Differentially expressed genes (DEGs) selected from GSE38267 were mostly enriched in inflammation‐related pathways, suggesting inflammation may be involved in the occurrence and development of PAH. Through the prediction and verification of related miRNAs and long noncoding RNAs using online databases and Gene Expression Omnibus (GEO) datasets, CCR7 and its related molecules, including hsa‐let‐7e‐5p and SNHG12, were identified as possible targets. The expression levels of CCR7, hsa‐let‐7e‐5p and SNHG12 were then verified by quantitative RT‐PCR in vivo and in vitro. Further study showed that silencing of SNHG12 decreased the expression of CCR7 under hypoxia treatment in vitro. Dual‐luciferase reporter assays were used to verify the relationship between hsa‐let‐7e‐5p and SNHG12. Collectively, our research reveals that a long noncoding RNA–miRNA–mRNA interaction network may be involved in the pathogenesis of PAH and suggests SNHG12, hsa‐let‐7e‐5p and CCR7 as potential biomarkers for PAH.

MiRNA-Mediated Subpathway Identification and Network Module Analysis to Reveal Prognostic Markers in Human Pancreatic Cancer

Background

Pancreatic cancer (PC) remains one of the most lethal cancers. In contrast to the steady increase in survival for most cancers, the 5-year survival remains low for PC patients.

Methods

We describe a new pipeline that can be used to identify prognostic molecular biomarkers by identifying miRNA-mediated subpathways associated with PC. These modules were then further extracted from a comprehensive miRNA-gene network (CMGN). An exhaustive survival analysis was performed to estimate the prognostic value of these modules.

Results

We identified 105 miRNA-mediated subpathways associated with PC. Two subpathways within the MAPK signaling and cell cycle pathways were found to be highly related to PC. Of the miRNA-mRNA modules extracted from CMGN, six modules showed good prognostic performance in both independent validated datasets.

Conclusions

Our study provides novel insight into the mechanisms of PC. We inferred that six miRNA-mRNA modules could serve as potential prognostic molecular biomarkers in PC based on the pipeline we proposed.

Replication Study: Coding-independent regulation of the tumor suppressor PTEN by competing endogenous mRNAs

As part of the Reproducibility Project: Cancer Biology, we published a Registered Report (Phelps et al., 2016) that described how we intended to replicate selected experiments from the paper 'Coding-independent regulation of the tumor suppressor PTEN by competing endogenous mRNAs' (Tay et al., 2011). Here, we report the results. We found depletion of putative PTEN competing endogenous mRNAs (ceRNAs) in DU145 cells did not impact PTEN 3'UTR regulation using a reporter, while the original study reported decreased activity when SERINC1, VAPA, and CNOT6L were depleted (Figure 3C; Tay et al., 2011). Using the same reporter, we found decreased activity when ceRNA 3'UTRs were overexpressed, while the original study reported increased activity (Figure 3D; Tay et al., 2011). In HCT116 cells, ceRNA depletion resulted in decreased PTEN protein levels, a result similar to the findings reported in the original study (Figure 3G,H; Tay et al., 2011); however, while the original study reported an attenuated ceRNA effect in microRNA deficient (Dicer^Ex5) HCT116 cells, we observed increased PTEN protein levels. Further, we found depletion of the ceRNAs VAPA or CNOT6L did not statistically impact DU145, wild-type HCT116, or Dicer^Ex5 HCT116 cell proliferation. The original study reported increased DU145 and wild-type HCT116 cell proliferation when these ceRNAs were depleted, which was attenuated in the Dicer^Ex5 HCT116 cells (Figure 5B; Tay et al., 2011). Differences between the original study and this replication attempt, such as variance between biological repeats, are factors that might have influenced the results. Finally, we report meta-analyses for each result.

Identification of differentially expressed subpathways via a bilevel consensus scoring of network topology and gene expression

Identifying differentially expressed subpathways connected to the emergence of a disease that can be considered as candidates for pharmacological intervention, with minimal off-target effects, is a daunting task. In this direction, we present a bilevel subpathway analysis method to identify differentially expressed subpathways that are connected with an experimental condition, while taking into account potential crosstalks between subpathways which arise due to their connectivity in a combined multi-pathway network. The efficacy of the method is demonstrated on a hematopoietic stem cell aging dataset, with findings corroborated using recent literature.

Exploring prognostic potential of long noncoding RNAs in colorectal cancer based on a competing endogenous RNA network

BACKGROUND

Colorectal cancer (CRC) is one of the most prevalent tumors worldwide. Recently, long noncoding RNAs (lncRNAs) have been shown to influence tumorigenesis and tumor progression by acting as competing endogenous RNAs (ceRNAs). It is difficult to extract prognostic lncRNAs and useful bioinformation from most ceRNA networks constructed previously.

AIM

To construct a prognostic related ceRNA regulatory network and lncRNA related signature based on risk score in CRC.

METHODS

RNA transcriptome profile and clinical information of 506 CRC patients were downloaded from the Cancer Genome Atlas database. R packages and Perl program were used for data processing. Cox regression analysis was used for prognostic model construction. Quantitative real-time polymerase chain reaction was used to detect the expression of lncRNAs.

RESULTS

A prognostic-related ceRNA network was constructed, including 9 lncRNAs, 44 mRNAs, and 30 miRNAs. In addition, a four-lncRNA model was constructed using multivariate Cox regression analysis, which could be an independent prognostic model in CRC. The risk score for each patient was calculated, and the 506 patients were divided into high and low-risk groups (253 for each group) based on the median risk score. The results of the survival analysis showed that patients with a high-risk score had a poor survival rate. Furthermore, the predictive value of the four-lncRNA model was evaluated in GSE38832. Patient survival probabilities could be better predicted when combing the risk score and clinical features. Gene Set Enrichment Analysis results verified that a number of cancer-related signaling pathways were enriched with a high-risk score in CRC. Finally, we validated a novel lncRNA (LINC00488) using quantitative real-time polymerase chain reaction in 22 paired CRC patient tumor tissues compared to adjacent non-tumor tissues.

CONCLUSION

The four-lncRNA model could give better predictive value for CRC patients. Our understanding of the lncRNA-related ceRNA regulatory mechanism could provide a potential diagnostic indicator for CRC patients.

ce-Subpathway: Identification of ceRNA-mediated subpathways via joint power of ceRNAs and pathway topologies

Competing endogenous RNAs (ceRNAs) represent a novel mechanism of gene regulation that may mediate key subpathway regions and contribute to the altered activities of pathways. However, the classical methods used to identify pathways fail to specifically consider ceRNAs within the pathways and key regions impacted by them. We proposed a powerful strategy named ce-Subpathway for the identification of ceRNA-mediated functional subpathways. It provided an effective level of pathway analysis via integrating ceRNAs, differentially expressed (DE) genes and their key regions within the given pathways. We respectively analysed one pulmonary arterial hypertension (PAH) and one myocardial infarction (MI) data sets and demonstrated that ce-Subpathway could identify many subpathways whose corresponding entire pathways were ignored by those non-ceRNA-mediated pathway identification methods. And these pathways have been well reported to be associated with PAH/MI-related cardiovascular diseases. Further evidence showed reliability of ceRNA interactions and robustness/reproducibility of the ce-Subpathway strategy by several data sets of different cancers, including breast cancer, oesophageal cancer and colon cancer. Survival analysis was finally applied to illustrate the clinical application value of the ceRNA-mediated functional subpathways using another data sets of pancreatic cancer. Comprehensive analyses have shown the power of a joint ceRNAs/DE genes and subpathway strategy based on their topologies.