Inference of patient‐specific subpathway activities reveals a functional signature associated with the prognosis of patients with breast cancer

Breast cancer is one of the most deadly forms of cancer in women worldwide. Better prediction of breast cancer prognosis is essential for more personalized treatment. In this study, we aimed to infer patient‐specific subpathway activities to reveal a functional signature associated with the prognosis of patients with breast cancer. We integrated pathway structure with gene expression data to construct patient‐specific subpathway activity profiles using a greedy search algorithm. A four‐subpathway prognostic signature was developed in the training set using a random forest supervised classification algorithm and a prognostic score model with the activity profiles. According to the signature, patients were classified into high‐risk and low‐risk groups with significantly different overall survival in the training set (median survival of 65 vs 106 months, P = 1.82e‐13) and test set (median survival of 75 vs 101 months, P = 4.17e‐5). Our signature was then applied to five independent breast cancer data sets and showed similar prognostic values, confirming the accuracy and robustness of the subpathway signature. Stratified analysis suggested that the four‐subpathway signature had prognostic value within subtypes of breast cancer. Our results suggest that the four‐subpathway signature may be a useful biomarker for breast cancer prognosis.     

Integrating high-throughput microRNA and mRNA expression data to identify risk mRNA signature for pancreatic cancer prognosis

Pancreatic cancer is a malignancy of the digestive system characterized by poor prognosis. A number of prognostic messenger RNA (mRNA) signatures have been identified by using the high-throughput expression profiles. MicroRNAs (miRNA) play a critical role in regulating multiple cellular functions. However, no such integrated analysis of miRNAs and mRNAs for studying the prognostic mechanisms of pancreatic cancer has been reported. In this study, we first identified prognostic mRNAs and miRNAs based on The Cancer Genome Atlas datasets, and then performed an enrichment analysis to explore the underlying biological mechanisms involved in pancreatic cancer prognosis at the mRNA level. Furthermore, we performed an integrated analysis of mRNAs and miRNAs to identify prognostic subpathways, which were closely associated with pancreatic cancer genes and tumor hallmarks and involved in hypoxia, oxidative phosphyorylation and xenobiotic metabolisms. Meanwhile, we performed a random walk algorithm based on global network, prognostic mRNAs and miRNAs, and identified top risk mRNAs as the prognostic signature. Finally, an independent testing set was used to confirm the predictive power of the top mRNA signature, and most of these genes involved were known oncogenes. In conclusion, we performed a series of integrated analyses by comprehensively exploring pancreatic cancer prognosis and systematically optimized the prognostic signature for clinical use.     

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.     

Identification of stably expressed reference small non‐coding RNAs for microRNA quantification in high‐grade serous ovarian carcinoma tissues

MicroRNAs (miRNAs) belong to a family of small non‐coding RNAs (sncRNAs) playing important roles in human carcinogenesis. Multiple investigations reported miRNAs aberrantly expressed in several cancers, including high‐grade serous ovarian carcinoma (HGS‐OvCa). Quantitative PCR is widely used in studies investigating miRNA expression and the identification of reliable endogenous controls is crucial for proper data normalization. In this study, we aimed to experimentally identify the most stable reference sncRNAs for normalization of miRNA qPCR expression data in HGS‐OvCa. Eleven putative reference sncRNAs for normalization (U6, SNORD48, miR‐92a‐3p, let‐7a‐5p, SNORD61, SNORD72, SNORD68, miR‐103a‐3p, miR‐423‐3p, miR‐191‐5p, miR‐16‐5p) were analysed on a total of 75 HGS‐OvCa and 30 normal tissues, using a highly specific qPCR. Both the normal tissues considered to initiate HGS‐OvCa malignant transformation, namely ovary and fallopian tube epithelia, were included in our study. Stability of candidate endogenous controls was evaluated using an equivalence test and validated by geNorm and NormFinder algorithms. Combining results from the three different statistical approaches, SNORD48 emerged as stably and equivalently expressed between malignant and normal tissues. Among malignant samples, considering groups based on residual tumour, miR‐191‐5p was identified as the most equivalent sncRNA. On the basis of our results, we support the use of SNORD48 as best reference sncRNA for relative quantification in miRNA expression studies between HGS‐OvCa and normal controls, including the first time both the normal tissues supposed to be HGS‐OvCa progenitors. In addition, we recommend miR‐191‐5p as best reference sncRNA in miRNA expression studies with prognostic intent on HGS‐OvCa tissues.     

The implications of relationships between human diseases and metabolic subpathways

One of the challenging problems in the etiology of diseases is to explore the relationships between initiation and progression of diseases and abnormalities in local regions of metabolic pathways. To gain insight into such relationships, we applied the "k-clique" subpathway identification method to all disease-related gene sets. For each disease, the disease risk regions of metabolic pathways were then identified and considered as subpathways associated with the disease. We finally built a disease-metabolic subpathway network (DMSPN). Through analyses based on network biology, we found that a few subpathways, such as that of cytochrome P450, were highly connected with many diseases, and most belonged to fundamental metabolisms, suggesting that abnormalities of fundamental metabolic processes tend to cause more types of diseases. According to the categories of diseases and subpathways, we tested the clustering phenomenon of diseases and metabolic subpathways in the DMSPN. The results showed that both disease nodes and subpathway nodes displayed slight clustering phenomenon. We also tested correlations between network topology and genes within disease-related metabolic subpathways, and found that within a disease-related subpathway in the DMSPN, the ratio of disease genes and the ratio of tissue-specific genes significantly increased as the number of diseases caused by the subpathway increased. Surprisingly, the ratio of essential genes significantly decreased and the ratio of housekeeping genes remained relatively unchanged. Furthermore, the coexpression levels between disease genes and other types of genes were calculated for each subpathway in the DMSPN. The results indicated that those genes intensely influenced by disease genes, including essential genes and tissue-specific genes, might be significantly associated with the disease diversity of subpathways, suggesting that different kinds of genes within a disease-related subpathway may play significantly differential roles on the diversity of diseases caused by the corresponding subpathway.     

Reconstruction and analysis of the aberrant lncRNA‐miRNA‐mRNA network based on competitive endogenous RNA in CESC

A growing body of studies has demonstrated that long non‐coding RNA (lncRNA) are regarded as the primary section of the ceRNA network. This is thought to be the case owing to its regulation of protein‐coding gene expression by functioning as miRNA sponges. However, functional roles and regulatory mechanisms of lncRNA‐mediated ceRNA in cervical squamous cell carcinoma (CESC), as well as their use for potential prediction of CESC prognosis, remains unknown. The aberrant expression profiles of mRNA, lncRNA, and miRNA of 306 cervical squamous cancer tissues and three adjacent cervical tissues were obtained from the TCGA database. A lncRNA‐mRNA‐miRNA ceRNA network in CESC was constructed. Meanwhile, Gene Ontology (GO) and KEGG pathway analysis were performed using Cytoscape plug‐in BinGo and DAVID database. We identified a total of 493 lncRNA, 70 miRNA, and 1921 mRNA as differentially expressed profiles. An aberrant lncRNA‐mRNA‐miRNA ceRNA network was constructed in CESC, it was composed of 50 DElncRNA, 18 DEmiRNA, and 81 DEmRNA. According to the overall survival analysis, 3 out of 50 lncRNA, 10 out of 81 mRNA, and 1 out of 18 miRNA functioned as prognostic biomarkers for patients with CESC (P value < 0.05). We extracted the sub‐network in the ceRNA network and found that two novel lncRNA were recognized as key genes. These included lncRNA MEG3 and lncRNA ADAMTS9‐AS2. The present study provides a new insight into a better understanding of the lncRNA‐related ceRNA network in CESC, and the novel recognized ceRNA network will help us to improve our understanding of lncRNA‐mediated ceRNA regulatory mechanisms in the pathogenesis of CESC.     

Comprehensive analysis of the functional microRNA–mRNA regulatory network identifies miRNA signatures associated with glioma malignant progression

Glioma is the most common and fatal primary brain tumour with poor prognosis; however, the functional roles of miRNAs in glioma malignant progression are insufficiently understood. Here, we used an integrated approach to identify miRNA functional targets during glioma malignant progression by combining the paired expression profiles of miRNAs and mRNAs across 160 Chinese glioma patients, and further constructed the functional miRNA–mRNA regulatory network. As a result, most tumour-suppressive miRNAs in glioma progression were newly discovered, whose functions were widely involved in gliomagenesis. Moreover, three miRNA signatures, with different combinations of hub miRNAs (regulations≥30) were constructed, which could independently predict the survival of patients with all gliomas, high-grade glioma and glioblastoma. Our network-based method increased the ability to identify the prognostic biomarkers, when compared with the traditional method and random conditions. Hsa-miR-524-5p and hsa-miR-628-5p, shared by these three signatures, acted as protective factors and their expression decreased gradually during glioma progression. Functional analysis of these miRNA signatures highlighted their critical roles in cell cycle and cell proliferation in glioblastoma malignant progression, especially hsa-miR-524-5p and hsa-miR-628-5p exhibited dominant regulatory activities. Therefore, network-based biomarkers are expected to be more effective and provide deep insights into the molecular mechanism of glioma malignant progression.     

Dissection of human MiRNA regulatory influence to subpathway

The global insight into the relationships between miRNAs and their regulatory influences remains poorly understood. And most of complex diseases may be attributed to certain local areas of pathway (subpathway) instead of the entire pathway. Here, we reviewed the studies on miRNA regulations to pathways and constructed a bipartite miRNAs and subpathways network for systematic analyzing the miRNA regulatory influences to subpathways. We found that a small fraction of miRNAs were global regulators, environmental information processing pathways were preferentially regulated by miRNAs, and miRNAs had synergistic effect on regulating group of subpathways with similar function. Integrating the disease states of miRNAs, we also found that disease miRNAs regulated more subpathways than nondisease miRNAs, and for all miRNAs, the number of regulated subpathways was not in proportion to the number of the related diseases. Therefore, the study not only provided a global view on the relationships among disease, miRNA and subpathway, but also uncovered the function aspects of miRNA regulations and potential pathogenesis of complex diseases. A web server to query, visualize and download for all the data can be freely accessed at http://bioinfo.hrbmu.edu.cn/miR2Subpath.     

Systematic construction and validation of an immune prognostic model for lung adenocarcinoma

Lung adenocarcinoma (LUAD), the most common non‐small‐cell lung cancer, is characterized by a dense lymphocytic infiltrate, which indicates that the immune system plays an active role in the development and growth of this cancer. However, no investigations to date have proposed robust models for predicting survival outcome for patients with LUAD in terms of tumour immunology. A total of 761 LUAD patients were included in this study, in which the database of The Cancer Genome Atlas (TCGA) was utilized for discovery, and the Gene Expression Omnibus (GEO) database was utilized for validation. Bioinformatics analysis and R language tools were utilized to construct an immune prognostic model and annotate biological functions. Lung adenocarcinoma showed a weakened immune phenotype compared with adjacent normal tissues. Immune‐related gene sets were profiled, an immune prognostic model based on 2 immune genes (ANLN and F2) was developed with the TCGA database to distinguish cases as having a low or high risk of unfavourable prognosis, and the model was verified with the GEO database. The model was prognostically significant in stratified cohorts, including stage I‐II, stage III‐IV and epidermal growth factor receptor (EGFR) mutant subsets, and was considered to be an independent prognostic factor for LUAD. Furthermore, the low‐ and high‐risk groups showed marked differences in tumour‐infiltrating leucocytes, tumour mutation burden, aneuploidy and PD‐L1 expression. In conclusion, an immune prognostic model was proposed for LUAD that is capable of independently identifying patients at high risk for poor survival, suggesting a relationship between local immune status and prognosis.     

Integrative analysis from multi‐centre studies identifies a function‐derived personalized multi‐gene signature of outcome in colorectal cancer

Colorectal cancer (CRC) is highly heterogeneous leading to variable prognosis and treatment responses. Therefore, it is necessary to explore novel personalized and reproducible prognostic signatures to aid clinical decision‐making. The present study combined large‐scale gene expression profiles and clinical data of 1828 patients with CRC from multi‐centre studies and identified a personalized gene prognostic signature consisting of 46 unique genes (called function‐derived personalized gene signature [FunPGS]) from an integrated statistics and function‐derived perspective. In the meta‐training and multiple independent validation cohorts, the FunPGS effectively discriminated patients with CRC with significantly different prognosis at the individual level and remained as an independent factor upon adjusting for clinical covariates in multivariate analysis. Furthermore, the FunPGS demonstrated superior performance for risk stratification with respect to other recently reported signatures and clinical factors. The complementary value of the molecular signature and clinical factors was further explored, and it was observed that the composite signature called IMCPS greatly improved the predictive performance of survival estimation relative to molecular signatures or clinical factors alone. With further prospective validation in clinical trials, the FunPGS may become a promising and powerful personalized prognostic tool for stratifying patients with CRC in order to achieve an optimal systemic therapy.     

Integrated analysis of co‐expression and ceRNA network identifies five lncRNAs as prognostic markers for breast cancer

Long non‐coding RNAs (lncRNAs), which competitively bind miRNAs to regulate target mRNA expression in the competing endogenous RNAs (ceRNAs) network, have attracted increasing attention in breast cancer research. We aim to find more effective therapeutic targets and prognostic markers for breast cancer. LncRNA, mRNA and miRNA expression profiles of breast cancer were downloaded from TCGA database. We screened the top 5000 lncRNAs, top 5000 mRNAs and all miRNAs to perform weighted gene co‐expression network analysis. The correlation between modules and clinical information of breast cancer was identified by Pearson's correlation coefficient. Based on the most relevant modules, we constructed a ceRNA network of breast cancer. Additionally, the standard Kaplan‐Meier univariate curve analysis was adopted to identify the prognosis of lncRNAs. Ultimately, a total of 23 and 5 modules were generated in the lncRNAs/mRNAs and miRNAs co‐expression network, respectively. According to the Green module of lncRNAs/mRNAs and Blue module of miRNAs, our constructed ceRNA network consisted of 52 lncRNAs, 17miRNAs and 79 mRNAs. Through survival analysis, 5 lncRNAs (AL117190.1, COL4A2‐AS1, LINC00184, MEG3 and MIR22HG) were identified as crucial prognostic factors for patients with breast cancer. Taken together, we have identified five novel lncRNAs related to prognosis of breast cancer. Our study has contributed to the deeper understanding of the molecular mechanism of breast cancer and provided novel insights into the use of breast cancer drugs and prognosis.     

Characterizing the Network of Drugs and Their Affected Metabolic Subpathways

A fundamental issue in biology and medicine is illustration of the overall drug impact which is always the consequence of changes in local regions of metabolic pathways (subpathways). To gain insights into the global relationship between drugs and their affected metabolic subpathways, we constructed a drug–metabolic subpathway network (DRSN). This network included 3925 significant drug–metabolic subpathway associations representing drug dual effects. Through analyses based on network biology, we found that if drugs were linked to the same subpathways in the DRSN, they tended to share the same indications and side effects. Furthermore, if drugs shared more subpathways, they tended to share more side effects. We then calculated the association score by integrating drug-affected subpathways and disease-related subpathways to quantify the extent of the associations between each drug class and disease class. The results showed some close drug–disease associations such as sex hormone drugs and cancer suggesting drug dual effects. Surprisingly, most drugs displayed close associations with their side effects rather than their indications. To further investigate the mechanism of drug dual effects, we classified all the subpathways in the DRSN into therapeutic and non-therapeutic subpathways representing drug therapeutic effects and side effects. Compared to drug side effects, the therapeutic effects tended to work through tissue-specific genes and these genes tend to be expressed in the adrenal gland, liver and kidney; while drug side effects always occurred in the liver, bone marrow and trachea. Taken together, the DRSN could provide great insights into understanding the global relationship between drugs and metabolic subpathways.     

Identification of prognostic signature in cancer based on DNA methylation interaction network

Background

The identification of prognostic biomarkers for cancer patients is essential for cancer research. These days, DNA methylation has been proved to be associated with cancer prognosis. However, there are few methods which identify the prognostic markers based on DNA methylation data systematically, especially considering the interaction among DNA methylation sites.

Methods

In this paper, we first evaluated the stabilities of microRNA, mRNA, and DNA methylation data in prognosis of cancer. After that, a rank-based method was applied to construct a DNA methylation interaction network. In this network, nodes with the largest degrees (10% of all the nodes) were selected as hubs. Cox regression was applied to select the hubs as prognostic signature. In this prognostic signature, DNA methylation levels of each DNA methylation site are correlated with the outcomes of cancer patients. After obtaining these prognostic genes, we performed the survival analysis in the training group and the test group to verify the reliability of these genes.

Results

We applied our method in three cancers (ovarian cancer, breast cancer and Glioblastoma Multiforme).In all the three cancers, there are more common ones of prognostic genes selected from different samples in DNA methylation data, compared with gene expression data and miRNA expression data, which indicates the DNA methylation data may be more stable in cancer prognosis. Power-law distribution fitting suggests that the DNA methylation interaction networks are scale-free. And the hubs selected from the three networks are all enriched by cancer related pathways. The gene signatures were obtained for the three cancers respectively, and survival analysis shows they can distinguish the outcomes of tumor patients in both the training data sets and test data sets, which outperformed the control signatures.

Conclusions

A computational method was proposed to construct DNA methylation interaction network and this network could be used to select prognostic signatures in cancer.