Next Article in Journal
Convergent Evolution of the Seed Shattering Trait
Next Article in Special Issue
PVTree: A Sequential Pattern Mining Method for Alignment Independent Phylogeny Reconstruction
Previous Article in Journal
Writing Histone Monoubiquitination in Human Malignancy—The Role of RING Finger E3 Ubiquitin Ligases
Previous Article in Special Issue
A Multi-Label Supervised Topic Model Conditioned on Arbitrary Features for Gene Function Prediction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Random Walk Based Cluster Ensemble Approach for Data Integration and Cancer Subtyping

1
College of Computer Science and Technology, Anhui University, Hefei 230601, Anhui, China
2
School of Software Engineering, Qufu Normal University, Qufu 273165, Shandong, China
3
Co-Innovation Center for Information Supply & Assurance Technology, Anhui University, Hefei 230601, Anhui, China
*
Author to whom correspondence should be addressed.
Submission received: 27 November 2018 / Revised: 11 January 2019 / Accepted: 14 January 2019 / Published: 18 January 2019

Abstract

:
Availability of diverse types of high-throughput data increases the opportunities for researchers to develop computational methods to provide a more comprehensive view for the mechanism and therapy of cancer. One fundamental goal for oncology is to divide patients into subtypes with clinical and biological significance. Cluster ensemble fits this task exactly. It can improve the performance and robustness of clustering results by combining multiple basic clustering results. However, many existing cluster ensemble methods use a co-association matrix to summarize the co-occurrence statistics of the instance-cluster, where the relationship in the integration is only encapsulated at a rough level. Moreover, the relationship among clusters is completely ignored. Finding these missing associations could greatly expand the ability of cluster ensemble methods for cancer subtyping. In this paper, we propose the RWCE (Random Walk based Cluster Ensemble) to consider similarity among clusters. We first obtained a refined similarity between clusters by using random walk and a scaled exponential similarity kernel. Then, after being modeled as a bipartite graph, a more informative instance-cluster association matrix filled with the aforementioned cluster similarity was fed into a spectral clustering algorithm to get the final clustering result. We applied our method on six cancer types from The Cancer Genome Atlas (TCGA) and breast cancer from the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC). Experimental results show that our method is competitive against existing methods. Further case study demonstrates that our method has the potential to find subtypes with clinical and biological significance.

1. Introduction

With the efforts of the large-scale projects such as The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) [1,2,3], a wealth of genome-scale molecular data are available and easy to access. The multiple types of omics data from genomes, transcriptomes, proteome, and epigenomes enable researchers to embrace great opportunities and possibilities to explore a more comprehensive view into cancer informatics, such as drug target prediction [4,5], diver gene identification [6,7,8], and so on.
One essential topic in oncology is cancer subtyping, whereby tumors are divided into clinically and biologically relevant subtypes, which could offer insight into tumor progression and provide personalized treatment. However, applying traditional clustering algorithms on a single data type—like gene expression data—does not obtain satisfactory results such as deriving subtypes associated with clinical phenotype [9]. These unsatisfactory results indicate the limitations of expression-based analysis for cancer subtyping. Since different types of molecular data contain information in various aspects that may complement each other, it is beneficial for leveraging different types of omics data simultaneously [10]. Several integrative frameworks have been proposed and gained success [9,10,11,12,13,14].
A promising method for cancer subtyping is cluster ensemble [11,15,16]. It can merge individual clusterings (clustering results obtained from running diverse clustering algorithms, running different types of omics data, etc.) to a consensus to form one robust unit. More importantly, cluster ensemble can naturally be applied on multiple data types as an integrative method. However, traditional cluster ensemble mostly merges different clusterings using a co-association matrix, which measures the frequency of two instances clustering together [17]. In this coarse way, some important information—such as relations among clusters—may be lost after merging the base clusterings. Link-based cluster ensemble (LCE) [15] tries to solve this problem by considering the relationship among clusters in terms of the triplet. The triplet is a subgraph containing three vertices and two non-zero edges. The similarity between two clusters is measured based on the count of all triples between them. However, this only captures local structure since the triplet measures the similarity in a local range.
In this paper, we proposed a new method named Random Walk based Cluster Ensemble (RWCE) to deal with these problems (Figure 1). We first obtained a refined cluster-cluster similarity by using random walk on a network of clusters constructed with Jaccard similarity and applied a scaled exponential similarity kernel, which provided a more global view from the whole cluster network. We then generated a more informative instance-cluster association matrix by filling in the refined cluster-cluster similarity. A bipartite graph was modeled on this resulting matrix in which spectral clustering [18] was used to obtain the final partition. Experiments on six cancer type datasets from the TCGA and the Molecular Taxonomy of Breast Cancer International (METABRIC) breast cancer data set [19] showed that our RWCE was competitive compared with other methods. Further case study demonstrated that our method also had the power to find clinically and biologically relevant subtypes. The source code of RWCE can be found in supplementary File 1.

2. Materials and Methods

2.1. Datasets

In order to show the effectiveness of our method, we used six TCGA cancer types: kidney renal clear cell carcinoma (KIRC), glioblastoma multiforme (GBM), lung squamous cell carcinoma (LUSC), breast invasive carcinoma (BRCA), acute myeloid leukemia (LAML), and colon adenocarcinoma (COAD). The data were processed by PINS (perturbation clustering for data integration and disease subtyping) [20], which is an integrative clustering framework for cancer subtyping. Each cancer type had three molecular omics data types, namely mRNA expression, miRNA expression, and DNA methylation. In addition, the METABRIC breast cancer data set [19] was used for survival analysis. The METABRIC data set included a discovery cohort (997 patients) and a validation cohort (995 patients). Each of them had two molecular data types: mRNA expression and copy number variation data, which were downloaded from the European Genome-Phenome Archive [21] (https://ega-archive.org/).

2.2. Competitive Methods

To show the effectiveness of our method, we compared it with a traditional cluster ensemble method called consensus clustering (CC) [17] as the baseline, and three state-of-the-art methods called link-based cluster ensemble (LCE) [15], perturbation clustering for data integration and disease subtyping (PINS) [20], and entropy-based consensus clustering (ECC) [11].

2.3. Evaluation Metrics

We used cox log-rank p-value [22] for measuring the significance of the difference of survival distributions between subtypes. Normally, a p-value < 0.05 indicates statistical significance, and a lower p-value indicates a more significant difference. We also used the silhouette value to measure consistency within subtypes. The mean value of the silhouette was used as a measure of how tightly grouped all the data in the cluster were. Higher values of the silhouette indicate a well-divided clustering structure.
For survival analysis, we also used concordance index (CI) [23], which measures the consistency between the estimated risk and the real survival time. Higher CI value indicates better performance for survival analysis.

2.4. Methodology Overview of RWCE

Here, we sum up RWCE for cancer subtyping. Suppose we have three data types to use for clustering. There are three steps in the RWCE pipeline. Step 1: For each data type, M basic clusterings are generated using K-means with a number of clusters randomly chosen from 2 to n , where n is the number of instances (Figure 1A). Note that in this step, we can use any clustering method, and in this paper, we fixed it to K-means. Step 2: These M basic clusterings are combined into a consensus clustering by RWCE refinement (Figure 1B), which we introduce in detail later. Step 3: Each data type follows the same operation as in Step 1 and Step 2, and we then have three consensus clusterings— π 1 * , π 2 * , π 3 * . At last, we use RWCE refinement again to combine each data type’s consensus clustering to get the final clustering result π * (Figure 1C).

2.5. Cluster Ensemble

Let X denote an omics data set such as gene expression data with n instances (or conditions, experiments, patients, and so on) and let m denote genes (or biomarkers and so on). A cluster ensemble is a set of M basic clustering solutions generated by different clustering algorithms or a single clustering algorithm with different parameters, which is represented as Π = { π ( 1 ) , π ( 2 ) , π ( M ) } . Each clustering π ( m ) partitions X into K m crisp clusters, represented as π ( m ) = { C 1 ( m ) , C 2 ( m ) , C K m ( m ) } , with C k ( m ) C k ( m ) = ,   k k , and k = 1 K m C k ( m ) = X . The cluster ensemble method then takes these clusterings Π as input and combines these solutions to produce the consensus clustering π * as the output. There are diverse ways for combing [24,25]. One can derive an instance-cluster binary (IC) matrix with ‘1’, indicating that instance belongs to that cluster, otherwise it is indicated as ‘0’. Then, a clustering algorithm or graph segmentation algorithm could be used on this matrix to get the consensus clustering solution [17].

2.6. RWCE Refinement

2.6.1. Generating a Refined Instance-Cluster Association (RIC) Matrix

A problem remains when leveraging an IC matrix or another similarity matrix since they only summarize information at a coarse level. For example, in an IC matrix, only one element is ‘1’ for one instance in each clustering, and others are ‘0’. This may lead to sparsity and does not favor the similarity-based clustering algorithm. Accordingly, in LCE [15], an improved variation of the original IC matrix, refined cluster-association matrix (RM) is generated by modifying the zero entries of the IC matrix with the cluster-cluster similarity discovered by the link-based similarity algorithm. The results show that the refinement is helpful and works better than using the original matrix. However, the algorithm used for measuring the similarity among clusters through focusing on triple is limited in local view. In response, we proposed RWCE, which has a more global view in discovering cluster-cluster similarity. We put forward a refined instance-cluster association (RIC) matrix as a more informative variation of the original IC matrix. It is designed to replace the value of those hidden associations (‘0’) of the IC matrix with the refined cluster similarity. For each clustering π ( m ) ,   m = 1 M and their corresponding clusters C 1 ( m ) , C 2 ( m ) , C K m ( m ) (where K m is the number of clusters in clustering π ( m ) ), the association R I C ( x i , C ) 0   a n d 1 between instance x i X and cluster C { C 1 ( m ) , C 2 ( m ) , C K m ( m ) } is measured as follows:
R I C ( x i , C ) = { 1   i f   C = C * ( m ) ( x i ) s i m ( C , C * ( m ) ( x i ) ) C π ( m ) Ȭ C C * ( m ) ( x i ) s i m ( C , C * ( m ) ( x i ) ) × d c   o t h e r w i s e
where C * ( m ) ( x i ) is a cluster label to which the instance x i belongs in clustering π ( m ) . Moreover, s i m ( C x , C y ) [ 0 , 1 ] measures the similarity between any two clusters C x , C y , which can be calculated using the random-walk based similarity algorithms listed in Section 2.6.2, and d c is a hyperparameter that we empirically set as 1 (performance is robust to dc, thus we fixed it to 1 for the sake of explanation). In this way, we fill in the zero entries of the IC matrix with the normalized similarity between the clusters by using the following random-walk based similarity algorithm.

2.6.2. Random-Walk Based Similarity Algorithm

We first constructed an original cluster-cluster similarity network by using the Jaccard index as follows:
J x y = | L x L y | | L x L y |
where J x y is an edge of the above similarity network between cluster C x and C y , and L x and L y denote the set of samples of clusters C x and C y , respectively. On this initial network, we applied random walk with restart:
F t + 1 = α F t A + ( 1 α ) F 0
where A is the adjacency matrix of the above-mentioned similarity network and F 0 is the IC matrix. ( 1 α ) is the restart probability that the random walker may choose to teleport to the initial node. The random walk process runs iteratively until F t + 1 converges ( | F t + 1 F t | < 1   × 10 6 ). In consequence, the resulted F t + 1 is a real-valued instance-by-cluster association matrix instead of a binary value, on which we can measure the refined similarity between clusters using the scaled exponential similarity kernel:
s i m ( C i , C j ) = e x p ( ρ 2 ( z i , z j ) 2 σ 2 )
where z i and z j are i-th column and j-th column of F t + 1 , representing clusters C i and C j , respectively, ρ 2 ( z i , z j ) denotes the squared Euclidean distance between cluster C i and C j , and σ is a parameter we set to 1.

2.6.3. Applying Spectral Clustering to RIC

As a result, we obtained a refined and informative instance-cluster (RIC) matrix; R I C ( i , j ) 0   a n d   1 is a degree that instance i belongs to cluster j . We then modeled a bipartite graph G = ( V , W ) based on RIC, where V = V C V I . V C is the set of vertices, where each vertex corresponds to a cluster from ensemble Π ; V I is the set of vertices, where each vertex corresponds to an instance from data set X . W denotes a set of weighted edges that can be defined as follows:
W ( i , j ) = 0   i f   v e r t i c e s   v i , v j V C W ( i , j ) = 0   i f   v e r t i c e s   v i , v j V I W ( i , j ) = W ( j , i ) = R I C ( i , j )   i f   v e r t i c e s   v i V C   a n d   v j V I
Note that W can be written as W = [ 0 R I C T R I C 0 ] equivalently. Given such a graph, spectral graph partitioning (SPEC) [18] was then used to generate the final partition of X , denoted as π * . SPEC with normalized cut is simply described as follows. Given graph G = ( V , W ) , it first calculated the degree matrix D with degrees of each node on the diagonal. It then computed the Laplacian matrix L = D W . Next, the normalized Laplacian matrix D 1 2 L D 1 2 , with its K smallest eigenvalues λ 1 λ k and their corresponding eigenvectors u 1 , u 2 u K , were obtained. Then, a matrix U = [ u 1 , u 2 u K ] was formalized after being row normalized. At last, SPEC generated the final clustering result using K-means on U. More details can be found in [18]. We selected the number of clusters k = a r g   i > 1 m a x e i g e n g a p ( i ) , where e i g e n g a p ( i ) = λ i + 1 λ i . To sum up, we called the process of operating on ensemble Π and getting the final clustering π * as RWCE refinement.

2.7. Integrating Multiple Types of Omics Data for Subtyping

Suppose we had T types of omics data to integrate. For each type of omics data X t ,   t = 1 T , we obtained the corresponding clustering result π t * , t = 1…T. Then, we treated these clustering results as a new ensemble Π *   = { π 1 * , π 2 * , π T * } . Finally, we used RWCE refinement again to Π *   to get the final clustering result π * across all T data types.

3. Results

3.1. Evaluation on TCGA Cancer Data Sets

For each cancer type, we counted the number of significant survival analyss results based on three single molecular data types and the integration of the three data types.
According to Figure 2, our method outperformed other methods on both single data type and the integration. By integrating the three molecular data types, our method attained significant subtypes (p-value < 0.05) for all six cancer types (Table 1). This indicates the potential of leveraging multiple data types simultaneously for identifying meaningful subtypes and the power of RWCE as an integrative method.
Table 1 shows the cox log-rank p-value of RWCE on three molecular data types and their integration across six cancer types from TCGA. It indicates that RWCE is a good integrative method for combining multiple omics data for cancer subtype discovery.
In terms of silhouette value (Figure 3), our method still outperformed other methods, indicating good clustering performance at the data level.

3.2. A Case Study: Glioblastoma Multiforme

Our method found three GBM subtypes. The survival curves of them are shown in Figure 4.
From Figure 4, subtype 1 had a bad prognosis while subtype 3 had a favorable prognosis. Moreover, Figure 5 shows that patients from subtype 1 had a favorable response to temozolomide (TMZ), a drug commonly used to treat GBM, and subtype 3 consisted of slightly younger patients.

3.3. Evaluation on METABRIC Data Set

We also tested the performance of survival analysis on the METABRIC breast cancer data set. As seen in Table 2, our method outperformed other clustering methods and was comparable with the PAM50 analysis (a standard breast cancer signature). This indicates the potential of our method for finding subtypes with differential survival profiles.

4. Discussion and Conclusions

In this paper, a new cluster ensemble method named RWCE was introduced for clustering and integrating multiple omics data to discover meaningful cancer subtypes. A novel RIC matrix is used in RWCE that considers relationships among clusters, which contributes to a superior clustering performance in terms of silhouette value and cox log-rank p-value.
Moreover, RWCE can also be utilized as an integrative method to make use of diverse types of omics data together for identifying subtypes with differential survival profiles. Further case study on the GBM subtypes that RWCE generated showed that RWCE could find subtypes with differential drug reactions and age distributions.
Taken together, RWCE provides a new way of thinking by combining basic clusterings in the cluster ensemble method and integrating multiple data types. We hope RWCE can generalize well to identify meaningful subtypes in more cancer types for the improvement of diagnostic and therapeutic intervention, and this is what we will investigate in further work.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4425/10/1/66/s1, File 1: Source code (ZIP, 7 KB).

Author Contributions

C.Y. carried out the experiments, analyses presented in this work and wrote the manuscript. Y.T.W. carried out the data analysis. Y.T.W. and C.H.Z. helped with project design, edited the manuscript and provided guidance and feedback throughout. All authors read and approved the final manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Nos. 61873001, 61872220, 61672037, 61861146002 and 61732012), the Key Project of Anhui Provincial Education Department (No. KJ2017ZD01).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. The International Cancer Genome Consortium. International network of cancer genome projects. Nature 2010, 464, 993. [Google Scholar] [CrossRef] [PubMed]
  2. Levine, D.A. The Cancer Genom Atlas Research Network. Integrated genomic characterization of endometrial carcinoma. Nature 2013, 497, 67. [Google Scholar]
  3. The Cancer Genom Atlas Research. Integrated genomic analyses of ovarian carcinoma. Nature 2011, 474, 609. [Google Scholar] [CrossRef] [PubMed]
  4. Emig, D.; Ivliev, A.; Pustavalova, O.; Lanchasire, L.; Bureeva, S.; Nikolsky, Y.; Bessarabova, M. Drug target prediction and repositioning using an integrated network-based approach. PLoS ONE 2013, 8, e60618. [Google Scholar] [CrossRef] [PubMed]
  5. Yamanishi, Y.; Araki, M.; Gutteridge, A.; Honda, W.; Kanehisa, M. Prediction of drug–target interaction networks from the integration of chemical and genomic spaces. Bioinformatics 2008, 24, i232–i240. [Google Scholar] [CrossRef] [PubMed]
  6. Bashashati, A.; Haffari, G.; Ding, J.; Ha, G.; Lui, K.; Rosner, J.; Huntsman, D.G.; Caldas, C.; Aparico, S.A.; Shah, S.P. DriverNet: Uncovering the impact of somatic driver mutations on transcriptional networks in cancer. Genome Biol. 2012, 13, R124. [Google Scholar] [CrossRef] [PubMed]
  7. Cho, A.; Shim, J.E.; Kim, E.; Supek, F.; Lehner, B.; Lee, I. MUFFINN: Cancer gene discovery via network analysis of somatic mutation data. Genome Biol. 2016, 17, 129. [Google Scholar] [CrossRef]
  8. Hou, J.P.; Ma, J. DawnRank: Discovering personalized driver genes in cancer. Genome Med. 2014, 6, 56. [Google Scholar] [CrossRef]
  9. Hofree, M.; Shen, J.P.; Carter, H.; Gross, A.; Ideker, T. Network-based stratification of tumor mutations. Nat. Methods 2013, 10, 1108. [Google Scholar] [CrossRef]
  10. Shen, R.; Olshen, A.B.; Ladanyi, M. Integrative clustering of multiple genomic data types using a joint latent variable model with application to breast and lung cancer subtype analysis. Bioinformatics 2009, 25, 2906–2912. [Google Scholar] [CrossRef] [Green Version]
  11. Liu, H.; Zhao, R.; Fang, H.; Cheng, F.; Fu, Y.; Liu, Y.Y. Entropy-based consensus clustering for patient stratification. Bioinformatics 2017, 33, 2691–2698. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, B.; Mezlini, A.M.; Demir, F.; Fiume, M.; Tu, Z.; Brudno, M.; Haibe-Kins, B.; Goldenberg, A. Similarity network fusion for aggregating data types on a genomic scale. Nat. Methods 2014, 11, 333. [Google Scholar] [CrossRef] [PubMed]
  13. Liu, J.X.; Gao, Y.L.; Zheng, C.H.; Xu, Y.; Yu, J. Block-constraint robust principal component analysis and its application to integrated analysis of TCGA Data. IEEE Trans. Nanobiosci. 2016, 15, 510–516. [Google Scholar] [CrossRef]
  14. Liu, J.X.; Xu, Y.; Zheng, C.H.; Kong, H.; Lai, Z.H. RPCA-based tumor classification using gene expression data. IEEE/ACM Trans. Comput. Biol. Bioinform. 2015, 12, 964–970. [Google Scholar] [CrossRef]
  15. Iam-On, N.; Boongoen, T.; Garrett, S. LCE: A link-based cluster ensemble method for improved gene expression data analysis. Bioinformatics 2010, 26, 1513–1519. [Google Scholar] [CrossRef] [PubMed]
  16. Lock, E.F.; Dunson, D.B. Bayesian consensus clustering. Bioinformatics 2013, 29, 2610–2616. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Monti, S.; Tamayo, P.; Mesirov, J.; Golub, T. Consensus clustering: A resampling-based method for class discovery and visualization of gene expression microarray data. Mach. Learn. 2003, 52, 91–118. [Google Scholar] [CrossRef]
  18. Ng, A.Y.; Jordan, M.I.; Weiss, Y. On spectral clustering: Analysis and an algorithm. In Proceedings of the Advances in Neural Information Processing Systems, Vancouver, BC, Canada, 9–14 December 2002; pp. 849–856. [Google Scholar]
  19. Curtis, C.; Shah, S.B.; Chin, S.F.; Gulisa, T.; Rueda, O.M.; Dunning, M.J.; Speed, D.; Lynch, A.G.; Samarajiwa, S.; Yuan, Y.; et al. The genomic and transcriptomic architecture of 2000 breast tumours reveals novel subgroups. Nature 2012, 486, 346. [Google Scholar] [CrossRef]
  20. Nguyen, T.; Tagett, R.; Diaz, D.; Draghici, S. A novel approach for data integration and disease subtyping. Genome Res. 2017, 27, 2025. [Google Scholar] [CrossRef]
  21. Lappalainen, I.; Almeida-King, J.; Kumanduri, V.; Senf, A.; Spalding, J.D.; Ur-Rehman, S.; Saunders, G.; Kandasamy, J.; Caccamo, M.; Leinonen, R. The European genome-phenome archive of human data consented for biomedical research. Nat. Genet. 2015, 47, 692. [Google Scholar] [CrossRef]
  22. Hosmer, D.W.; Lemeshow, S.; May, S. Applied survival analysis: Regression modeling of time-to-event data, second edition. J. Stat. Plan. Inference 2000, 91, 173–175. [Google Scholar]
  23. Pencina, M.J.; D’Agostino, R.B. Overall C as a measure of discrimination in survival analysis: Model specific population value and confidence interval estimation. Stat. Med. 2004, 23, 2109–2123. [Google Scholar] [CrossRef] [PubMed]
  24. Strehl, A.; Ghosh, J. Cluster ensembles—A knowledge reuse framework for combining multiple partitions. JMLR 2003, 3, 583–617. [Google Scholar]
  25. Topchy, A.; Jain, A.K.; Punch, W. Clustering ensembles: Models of consensus and weak partitions. IEEE Trans. Pattern. Anal. Mach. Intell. 2005, 27, 1866–1881. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic diagram of the Random Walk based Cluster Ensemble (RWCE) pipeline: (A) traditional clustering algorithm (here we used K-means) was applied to each molecular data type to obtain M basic clusterings. For each basic clustering, the cluster number was randomly chosen from 2 to √n; (B) each data type’s M clusterings were fused into one consensus clustering by RWCE refinement; (C) all data types’ consensus clusterings were fused into one final clustering using RWCE refinement again.
Figure 1. Schematic diagram of the Random Walk based Cluster Ensemble (RWCE) pipeline: (A) traditional clustering algorithm (here we used K-means) was applied to each molecular data type to obtain M basic clusterings. For each basic clustering, the cluster number was randomly chosen from 2 to √n; (B) each data type’s M clusterings were fused into one consensus clustering by RWCE refinement; (C) all data types’ consensus clusterings were fused into one final clustering using RWCE refinement again.
Genes 10 00066 g001
Figure 2. Stacked histogram displaying, for each clustering method (PINS: perturbation clustering for data integration and disease subtyping; ECC: entropy-based consensus clustering; LCE: link-based cluster ensemble; CC: consensus clustering; RWCE: random walk based cluster ensemble), the times it passed the significant tests (p-value < 0.05) of survival analysis on several molecular data types: mRNA expression data (mRNA), DNA methylation data (Methy), miRNA expression data (miRNA) and an integration of all three data types (integration).
Figure 2. Stacked histogram displaying, for each clustering method (PINS: perturbation clustering for data integration and disease subtyping; ECC: entropy-based consensus clustering; LCE: link-based cluster ensemble; CC: consensus clustering; RWCE: random walk based cluster ensemble), the times it passed the significant tests (p-value < 0.05) of survival analysis on several molecular data types: mRNA expression data (mRNA), DNA methylation data (Methy), miRNA expression data (miRNA) and an integration of all three data types (integration).
Genes 10 00066 g002
Figure 3. The heatmap for silhouette value on six TCGA datasets of different methods. KIRC-mRNA indicates mRNA expression data in KIRC was used. The same as the others.
Figure 3. The heatmap for silhouette value on six TCGA datasets of different methods. KIRC-mRNA indicates mRNA expression data in KIRC was used. The same as the others.
Genes 10 00066 g003
Figure 4. The survival curves for TCGA glioblastoma multiforme (GBM) subtypes generated by RWCE.
Figure 4. The survival curves for TCGA glioblastoma multiforme (GBM) subtypes generated by RWCE.
Genes 10 00066 g004
Figure 5. (AC) Survival analysis of GBM patients for treatment with temozolomide (TMZ) in different subtypes generated by RWCE; (D) age distribution of GBM subtypes generated by RWCE.
Figure 5. (AC) Survival analysis of GBM patients for treatment with temozolomide (TMZ) in different subtypes generated by RWCE; (D) age distribution of GBM subtypes generated by RWCE.
Genes 10 00066 g005
Table 1. Performance of RWCE on three molecular data types and their integration across six cancer types from The Cancer Genome Atlas (TCGA).
Table 1. Performance of RWCE on three molecular data types and their integration across six cancer types from The Cancer Genome Atlas (TCGA).
mRNAMethylationmiRNAIntegration
KIRC0.008(2)0.79397(3)0.52883(2)0.00671(2)
GBM0.19041(2)0.00629(2)0.96568(2)0.00343 (3)
LAML0.00272(8)0.58721(2)0.00119(8)0.00158(2)
LUSC0.40747(3)0.04761(7)0.01666(2)0.00827(4)
BRCA0.04193(2)0.58412(2)0.15534(2)0.03006(2)
COAD0.01058(2)0.68703(2)0.81886(6)0.02818(3)
KIRC (kidney renal clear cell carcinoma); GBM (glioblastoma multiforme); LAML (acute myeloid leukemia); LUSC (lung squamous cell carcinoma); BRCA (breast invasive carcinoma); COAD (colon adenocarcinoma). p < 0.05 is highlighted in bold.
Table 2. Cox p-value and concordance index (CI) of subtypes discovered by PAM50, perturbation clustering for data integration and disease subtyping (PINS), consensus clustering (CC), entropy-based consensus clustering (ECC), link-based cluster ensemble (LCE), and our method on METABRIC data. For each discovery and validation cohort, we calculated the p-value and CI with respect to disease free survival (DFS) and overall survival of the patients. For each row, the best p-value (most significant) and the best CI (highest) are in red. The number of clusters in discovery and validation cohort are shown after the name of the clustering methods.
Table 2. Cox p-value and concordance index (CI) of subtypes discovered by PAM50, perturbation clustering for data integration and disease subtyping (PINS), consensus clustering (CC), entropy-based consensus clustering (ECC), link-based cluster ensemble (LCE), and our method on METABRIC data. For each discovery and validation cohort, we calculated the p-value and CI with respect to disease free survival (DFS) and overall survival of the patients. For each row, the best p-value (most significant) and the best CI (highest) are in red. The number of clusters in discovery and validation cohort are shown after the name of the clustering methods.
PAM50 (5, 5)PINS (14, 7)CC (10, 8)ECC (10, 10)LCE (10, 8)RWCE (6, 6)
Discoveryp-valueDFS3.00 × 10−116.50 × 10−102.50 × 10−51.39 × 10−19.50 × 10−11.69 × 10−9
Overall8.50 × 10−51.90 × 10−68.10 × 10−65.59 × 10−24.42 × 10−14.16 × 10−12
CIDFS0.6200.6340.5980.5210.5060.594
Overall0.5780.5980.5720.5290.5080.641
Validationp-valueDFS3.10 × 10−94.30 × 10−51.20 × 10−22.61 × 10−18.44 × 10−29.12 × 10−5
Overall2.90 × 10−5033.80 × 10−37.90 × 10−31.66 × 10−13.53 × 10−29.13 × 10−7
CIDFS0.6360.5890.5720.5210.5200.560
Overall0.5610.5450.5380.5190.5140.607

Share and Cite

MDPI and ACS Style

Yang, C.; Wang, Y.-T.; Zheng, C.-H. A Random Walk Based Cluster Ensemble Approach for Data Integration and Cancer Subtyping. Genes 2019, 10, 66. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10010066

AMA Style

Yang C, Wang Y-T, Zheng C-H. A Random Walk Based Cluster Ensemble Approach for Data Integration and Cancer Subtyping. Genes. 2019; 10(1):66. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10010066

Chicago/Turabian Style

Yang, Chao, Yu-Tian Wang, and Chun-Hou Zheng. 2019. "A Random Walk Based Cluster Ensemble Approach for Data Integration and Cancer Subtyping" Genes 10, no. 1: 66. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10010066

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop