GCNA-Cluster: A Gene Co-Expression Network Alignment to Cluster Cancer Patients Algorithm for Identifying Subtypes of Pancreatic Ductal Adenocarcinoma.

Cancer heterogeneity makes it necessary to use different treatment strategies for patients with the same pathological features. Accurate identification of cancer subtypes is a crucial step in this approach. The current studies of pancreatic ductal adenocarcinoma (PDAC) subtypes mainly focus on single genes and ignore the synergistic effects of genes. Here we proposed a network alignment algorithm GCNA-cluster to cluster patients based on gene co-expression networks. We constructed weighted gene co-expression networks for patients and aligned the networks of two patients to estimate the similarity of patients and their cancer subtypes. A scoring function is defined to measure the network alignment result and the score can indicate the similarity between patients. Then, the patients are clustered based on their similarities. We validated the accuracy of the algorithm on the GEO-PDAC dataset with real labels, and the experimental results show that the GCNA-cluster algorithm has better results than classical cancer subtyping algorithms. In addition, the GCNA-cluster algorithm applied to the TCGA-PDAC dataset identified two subtypes based on the Silhouette Coefficient. Biomarkers identified for the PDAC subtypes hint to cell growth, cell cycle or apoptosis as targets for new therapeutic strategies.

RSC-based differential model with correlation removal for improving multi-omics clustering.

Multi-omics clustering plays an important role in cancer subtyping. However, the data of different kinds of omics are often related, these correlations may reduce the clustering algorithm performance. It is crucial to eliminate the unexpected redundant information caused by these correlations between different omics. We proposed RSC-based differential model with correlation removal for improving multi-omics clustering (RSC-MCR). This method first introduced RSC to calculate the pairwise correlations of all features, and decomposed it to obtain the pairwise correlations of different omics features, thus built the connection between different omics based on the pairwise correlations of different omics features. Then, to remove the redundant correlation, we designed a differential model to calculate the degree of difference between the original feature matrix and the correlation matrix which contained the most relevant information between different omics. We compared the performance of RSC-MCR with decorrelation methods on different clustering methods (CC, FCM, SNF, NMF, LRAcluster). The experimental results on five cancer datasets show the efficiency of the RSC-MCR as well as improvements over other decorrelation methods.

A 70­RNA model based on SVR and RFE for predicting the pancreatic cancer clinical prognosis.

Researches on the prognosis of pancreatic cancer is of great significance to improve the patient treatment effect and survival. Current researches mainly focus on the prediction of the survival status and the determination of prognostic markers. Each patient has its own characteristics, there is no report about the prediction of survival time. However, accurate prediction of survival time is critical for personalized medicine. In this paper, a hybrid algorithm of Support Vector Regression (SVR) and Recursive Feature Elimination (RFE) was used to construct a quantitative prediction model of Overall Survival (OS) for pancreatic cancer patients, 70 RNAs related to OS were determined, including 33 mRNAs, 28 lncRNAs, and 9 miRNAs. The results of 10-fold cross-validation (R2 is 0.9693) and the generalization ability (R2 is 0.9666) showed that the model has reliable predictive performance and these 70 RNAs are important factors influencing the OS of pancreatic cancer patients. To further study the relationship between RNA-RNA interaction and the survival, competitive endogenous RNA (ceRNA) regulation network was constructed. Degree centrality, betweenness centrality and closeness centrality of nodes in the ceRNA network showed that hsa-mir-570, hsa-mir-944, hsa-mir-6506, hsa-mir-3136, MMP16, PLGLB2, HPGD, FUT1, MFSD2A, SULT1E1, SLC13A5, ZNF488, F2RL2, TNFRSF8, TNFSF11, FHDC1, ISLR2 and THSD7B are hub nodes, which are key RNAs closely determining the OS of pancreatic cancer patients.

Identification of subtype specific biomarkers of clear cell renal cell carcinoma using random forest and greedy algorithm.

Suitable biomarkers can be good indicator for cancer subtype. To find biomarkers that can accurately distinguish clear cell renal cell carcinoma (ccRCC) subtypes, we first determined ccRCC subtypes based on the expression of mRNA, miRNA and lncRNA, named clear cell type 1 (ccluster1) and 2 (ccluster2), using three unsupervised clustering algorithms. Besides being associated with the expression pattern derived from the single type of RNA, the differences between subtypes are relevant to the interactions between RNAs. Then, based on ceRNA network, the optimal combination features are selected using random forest and greedy algorithm. Further, in survival-related sub-ceRNA, competing gene pairs centering on miR-106a, miR-192, miR-193b, miR-454, miR-32, miR-98, miR-143, miR-145, miR-204, miR-424 and miR-1271 can also well identify ccluster1 and ccluster2 with prediction accuracy over 92%. These subtype-specific features potentially enhance the accuracy with which machine learning methods predict specific ccRCC subtypes. Simultaneously, the changes of miR-106 and OIP5-AS1 affect cell proliferation and the prognosis of ccluster1. The changes of miR-145 and FAM13A-AS1 in ccluster2 have an effect on cell invasion, apoptosis, migration and metabolism function. Here miR-192 displays a unique characteristic in both subtypes. Two subtypes also display notable differences in diverse pathways. Tumors belonging to ccluster1 are characterized by Fc gamma R-mediated phagocytosis pathway that affects tissue remodeling and repair, whereas those belonging to ccluster2 are characterized by EGFR tyrosine kinase inhibitor resistance pathway that participates in regulation of cell homeostasis. In conclusion, identifying these gene pairs can shed light on therapeutic mechanisms of ccRCC subtypes.