METHODS article

Front. Genet., 05 February 2019

Sec. Computational Genomics

Volume 10 - 2019 | https://doi.org/10.3389/fgene.2019.00020

Discovering Cancer Subtypes via an Accurate Fusion Strategy on Multiple Profile Data

  • LJ

    Limin Jiang 1†

  • YX

    Yongkang Xiao 2†

  • YD

    Yijie Ding 3

  • JT

    Jijun Tang 1,4*

  • FG

    Fei Guo 1*

  • 1. School of Computer Science and Technology, College of Intelligence and Computing, Tianjin University, Tianjin, China

  • 2. School of Chemical Engineering and Technology, Tianjin University, Tianjin, China

  • 3. School of Electronic and Information Engineering, Suzhou University of Science and Technology, Suzhou, China

  • 4. Department of Computer Science and Engineering, University of South Carolina, Columbia, SC, United States

Abstract

Discovering cancer subtypes is useful for guiding clinical treatment of multiple cancers. Progressive profile technologies for tissue have accumulated diverse types of data. Based on these types of expression data, various computational methods have been proposed to predict cancer subtypes. It is crucial to study how to better integrate these multiple profiles of data. In this paper, we collect multiple profiles of data for five cancers on The Cancer Genome Atlas (TCGA). Then, we construct three similarity kernels for all patients of the same cancer by gene expression, miRNA expression and isoform expression data. We also propose a novel unsupervised multiple kernel fusion method, Similarity Kernel Fusion (SKF), in order to integrate three similarity kernels into one combined kernel. Finally, we make use of spectral clustering on the integrated kernel to predict cancer subtypes. In the experimental results, the P-values from the Cox regression model and survival curve analysis can be used to evaluate the performance of predicted subtypes on three datasets. Our kernel fusion method, SKF, has outstanding performance compared with single kernel and other multiple kernel fusion strategies. It demonstrates that our method can accurately identify more accurate subtypes on various kinds of cancers. Our cancer subtype prediction method can identify essential genes and biomarkers for disease diagnosis and prognosis, and we also discuss the possible side effects of therapies and treatment.

1. Introduction

Cancer is a heterogeneous disease caused by chemical, physical, or genetic factors (Mager, 2006; Liu and Chu, 2014). The development of high-throughput genome analysis techniques on the research of cancer subtypes plays an important role in the analysis and clinical treatment of various kinds of cancers (Kruijf et al., 2013; Prat et al., 2015; Thanki et al., 2017). In recent years, much expression data, including genomes, transcriptome and epigenomes, has accumulated and been stored in various databases. The Cancer Genome Atlas (TCGA) (Katarzyna et al., 2015) is a large-scale project including over 34 cancers and 15 expression data sets. We can conveniently obtain genome-scale molecular data, which contributes to the development of computational methods for discovering cancer subtypes.

Until now, massive computational methods were proposed to discover cancer subtypes. Some methods are based on single expression data, including gene expression data (Nguyen and Rocke, 2002; Brunet et al., 2004; Finnegan and Carey, 2007; Teschendorff et al., 2007) and copy number (Wong et al., 2012) and DNA methylation (Zhang et al., 2017). Gao and Church (2005) employed sparse non-negative matrix factorization (SNMF) and gene expression data to identify subtypes of three cancers. Also, various kinds of expression data (Wei et al., 2017, 2018a,b) and several types of similarity strategies (Zeng et al., 2016; Ding et al., 2017a,b; Pan et al., 2017, 2018; Guo F. et al., 2018; Song et al., 2018) can be applied in many other biological prediction problems.

Generally, we desire a comprehensive view of one disease with a cohort of patients. We cannot analyze just one kind of data, but must separately abstract information from different types of data (Xu et al., 2017). Therefore, many methods improve the robustness of clustering by focusing on data processing (Ren et al., 2015). Wang et al. (2014) proposed the Similarity Network Fusion (SNF) approach for accurately clustering caner subtypes. This method first collects three types of genome-wide data including gene, methylation and miRNA expression. Then, it constructs the networks of samples (e.g., patients) by using three types of expression data, and fuses these networks into one network by using SNF representing the full spectrum of underlying data. Finally, it employs spectral clustering on an integrated network to predict caner subtypes. Ma and Zhang (2017) developed an improved SNF, Affinity Network Fusion (ANF), to integrate multiple similarity networks. Xu et al. (2016) proposed Weighted Similarity Network Fusion (WSNF) to identify cancer subtypes. This method constructs similarity of patients by integrating associations between miRNA, mRNA, and transcription factors. It is applied to two cancer types to demonstrate performance.

Furthermore, the effective models of clustering that we usually use have strong data sensitivity, such as k-means and hierarchical clustering. Today, many clustering methods have been developed to identify cancer subtypes. Le et al. (2016) developed the SRF algorithm, which identifies subtypes by combining mutational and expression information. It diffuses mutation information over an interaction network on the basis of each sample and eliminates scale differences by applying a rank-based transformation based on mutation and expression data. Then, rank matrix factorization is used to jointly factorize the transformed data into a number of ranked factors, and the subtypes are defined as the combination of ranked factors. This method obtains excellent performance, but some of the patients cannot be identified. Shen et al. (2009) proposed the iCluster method, which is based on the Gaussian latent variable model, to discover caner subtypes. This method was tested on breast cancer and lung cancer by using copy number and gene expression data types. Speicher and Pfeifer (2015) pointed out that iCluster has high computational complexity and proposed a dimensionality reduction method to integrate multiple similarity kernels. This method is evaluated by using five cancer types. Ge et al. (2017) developed the Scluster method, which integrates different types of data and maps them into an effective low-dimensional subspace. First, Scluster uses adaptive sparse reduced-rank regression (S-rrr) to map the original data into the principal subspaces. Next, a fused patient-by-patient network is abstracted for these subgroups by a scaled exponential similarity kernel method. It can then obtain the cancer subtypes by spectral clustering.

In this paper, we first collect multiple profile data on The Cancer Genome Atlas (TCGA), including five cancers (lung cancer, kidney cancer, stomach cancer, breast cancer, and colon cancer) and their three types of expression data (gene expression, isoform expression, and miRNA expression). Then, we construct three similarity kernels for all patients of the same cancer by using the three types of expression data. We then propose a novel unsupervised multiple kernel fusion method, Similarity Kernel Fusion (SKF), in order to integrate three similarity kernels into one combined kernel. Compared with SNF, SKF not only keeps the original information of each type of similarity kernel, but also gets rid of the noise in the integrated kernel. Finally, we make use of spectral clustering on the integrated kernel to predict cancer subtypes. To test the effectiveness and robustness of this novel approach, P- values from a Cox regression model and survival curve analysis can be used to evaluate the performance of our method on cancer subtype prediction. We compare the integrated kernel with the single kernel and other fusion methods, and also analyze the survival curve of the clinical data.

2. Materials and Methods

In this paper, we first extract five cancer datasets from The Cancer Genome Atlas (TCGA). For a particular cancer, we construct three patient similarity kernels by using the expression data. Then, we combine these similarity kernels into one similarity kernel by using Similarity Kernel Fusion (SKF). Finally, we employ spectral clustering on the integrated kernel to divide all patients into multiple clusters. The flowchart of our method is shown in Figure 1.

Figure 1

2.1. Dataset

We collect five cancer datasets from the TCGA website, including stomach cancer, lung cancer, kidney cancer, breast cancer, and colon cancer. For each cancer, we extract three kinds of expression data respectively, including gene expression, miRNA expression, and isoform level. Our dataset is denoted as Dataset No.1 in this paper. In addition, we employ anther dataset to evaluate the performance of our method. The second dataset is provided in Wang et al. (2014), which includes lung cancer, kidney cancer, breast cancer, colon cancer, and glioblastoma multiforme (GBM). For each tumor, gene expression, methylation expression, and miRNA expression from TCGA are used to analyze cancer subtypes. We denote this dataset as Dataset No.2. Since genes could be categorized into multiple groups, we selected 18222 coding genes from Dataset No.1, formed as Dataset No.3 . A summary of the three datasets is shown in Table 1. It is clear that Dataset No.1 and Dataset No.3 have more patients and expression factors than Dataset No.2 .

Table 1

DiseasesPatientsGenesIsoformmiRNAs
No.1 DatasetBreast1071604831831881
Colon426604831861881
Kidney868604831761881
Lung981604831741881
Stomach377604832111881
DiseasesPatientsGenesCpG sitesmiRNAs
No.2 DatasetBreast1051781423094354
Colon921781423088312
Kidney1221789924960329
Lung1061204223074352
GBM215120421305534
DiseasesPatientsGenesIsoformmiRNAs
No.3 DatasetBreast1071182221831881
Colon426182221861881
Kidney868182221761881
Lung981182221741881
Stomach377182222111881

Description of three datasets from TCGA.

2.2. Similarity Kernel Construction

A special expression dataset is denoted as E∈Rn×m, where m is the number of expression factors and n is the number of patients. We first normalize E by using Equation (1).

where x is an element of E, x′ is corresponding elements of E after standardization, is the mean of E and S is standard deviation of E. And, we denote normalized expression data as E′.

Based on the processed expression data E′, we construct similarity kernel K∈Rm×m for patients. Here, the similarity between two patients is defined as Equation (2) (Chen et al., 2018b; Zhao et al., 2018a,b).

where Ki,j is the similarity between i-th patient and j-th patient, and is i-th column and j-th column of E′, respectively.

Finally, we get three similarity kernels for a special disease, including similarity kernel by using gene expression, similarity kernel by using miRNA expression, and similarity kernel by using isoform expression.

2.3. Similarity Kernel Fusion

We constructed three similarity kernels for patients in the above section. We propose Similarity Kernel Fusion (SKF) to combine these kernels into one kernel K*∈Rm×m. First, we construct two kernels P∈Rm×m and S∈Rm×m for each similarity kernel by using Equations (3, 4), where P is a normalized kernel and S is a sparse kernel that eliminates weak similarity.

where P satisfies .

where S satisfies ; Ni is a set of all neighbors of the i-th patient, including itself.

Second, we discover more information by using multiple iterations as Equation (5).

where (l = 1, 2, 3) is the status of the l-th kernel after t iterations, α is a coefficient and satisfies α∈[0, 1], (r = 1, 2, 3) represents the initial status of Pr.

After t + 1 iterations, the overall kernel can be computed as Equation (6).

Finally, based on the integrated kernel, we construct a weight matrix to eliminate noise in the integrated kernel as Equation (7).

where Ni is a set of all neighbors of the i-th patient, including itself, and Nj is a set of all neighbors of the j-th patient, including itself.

The final similarity kernel can be obtained as Equation (8).

where K* is the final integrated similarity kernel by using SKF.

2.4. Mining Subtypes Using Spectral Clustering

In this section, we employ spectral clustering (Ng et al., 2001) on the integrated similarity kernel to divide all patients into multiple clusters. Many previous studies, including CSPRV (Guo Y. et al., 2018), Scluster (Ge et al., 2017), and SNF(Wang et al., 2014), have constructed similarity kernels for patients and used spectral clustering to discover cancer subtypes. These methods have achieved excellent performance by using spectral clustering. Additionally, Luxburg (2007) have pointed out that spectral clustering is effective in capturing the global structure of the graph. Therefore, we use spectral clustering to identify cancer subtypes. Then, we will introduce the processes of spectral clustering in detail. We define a matrix Y∈{0, 1}k×n to represent the result of a cluster, where Y(i, j) = 1 if patient pj belongs to i-th cluster, otherwise Y(i, j) = 0. We also use Equation (9) as the optimal question to solve Y.

where , , D is a diagonal matrix whose diagonal element is the sum of the row elements of K*.

3. Results

In this section, we discuss the performance of our method in a variety of ways. First, we introduce an evaluation criteria and a verification method that are used to evaluate the performance significance of the cancer subtype predictions. Second, we analyze the performance of SKF with different parameters α on Dataset No.1 . Third, we discuss the performance of SKF on the three datasets. Fourth, we compare SKF with two other fusion methods on the three datasets. Finally, we analyze the survival probability curves of the predicted subtypes for four cancers.

3.1. Evaluation Criteria and Verification Method

In this paper, we employ the P-value from the Cox regression model to evaluate the performance of our method, where a lower P-value indicates higher significance for performance. When the P-value is less than 0.05, it is of significance to the performance of the model. When the P-value is less than 0.01, the performance of the model is highly significant. Here, we use 0.05 as the threshold for significance. The meaning of the P-value is significance in the difference of survival profiles between cancer subtypes. Moreover, we also use survival analysis to evaluate the performance of the clustering results. The survival curve represents the change in survival rate over time, and it is a monotone decreasing curve without any fluctuation. In the survival curve, we can find that different subtypes have different survival rates. We can analyze some subtypes that have a higher risk of death.

3.2. Parameter Selection for SKF

Particularly, α is an important parameter in the process of SKF. A lower α value represents keeping more initial information in the integrated kernel. A higher α value represents keeping more information after multiple iterations. In the three datasets, we take α from 0 to 1 with a step of 0.1 to find the optimal α for the five cancers. Results are shown in Figure 2, with the X axis representing the α value and the Y axis representing the −log10(Pvalue). A lower P-value is represented by a higher value of −log10(Pvalue). In Figure 2, the P-value maintains clear fluctuation in the range between 0 and 1. It demonstrates that SKF is sensitive to changes in α. We get the optimal P-value when α is equal to 1 for the four cancers except lung cancer on Dataset No.3. From the results of Dataset No.2, we can see that keeping more initial information is necessary for many of the datasets.

Figure 2

3.3. Performance of SKF in Difference Datasets

In this paper, we obtain Dataset No.1 from TCGA. For a specific disease, we extract all 60483 gene expression data points on Dataset No.1 . We employ the three datasets to evaluate the performance of SKF. For each dataset, we compare the performance of SKF with single kernel by using the optimal number of clusters. In Table 2, we can see that SKF achieves outstanding performance compared with single kernel in 12 cases. We also find that the same kernels with different numbers of clusters have different P-values. Therefore, we need to adjust the number of clusters to obtain optimal clustering results. Although P-values do not achieve significant performance for GBM cancer in Dataset No.2 or Kidney cancer in Datasets No.2 and No.3 after SKF, these P-values get remarkable promotion compared to single kernel. Moreover, it is clear that the P-value of Dataset No.3 is better than that in Dataset No.1 , which shows that coding genes play an important role in the clustering of cancer subtypes.

Table 2

DatasetsCancersGenemiRNAIsoformSKF
expressionexpressionexpression
Dataset No.1Stomach (C=12)0.7030.0270.5488.86 × 10−14
Lung (C = 8)0.6210.1370.8293.81 × 10−4
Kidney (C = 3)0.2280.6420.3580.120
Breast (C = 5)0.5160.2810.2819.79 × 10−6
Colon (C = 6)0.0450.7260.1330.025
CancersGeneDNAmiRNASKF
expressionmethylationexpression
Dataset No.2GBM (C = 5)0.1590.0010.4360.037
Lung (C = 3)8.25 × 10−40.0090.2896.66 × 10−5
Kidney (C = 5)0.01770.4670.3680.0372
Breast (C = 5)0.0090.001641.38 × 10−42.7 × 10−7
Colon (C = 5)0.5870.0840.7021.81 × 10−3
CancersGenemiRNAIsoformSKF
expressionexpressionexpression
Dataset No.3Stomach (C = 9)0.05380.4380.6210.003
Lung (C = 3)0.3520.1710.3980.005
Kidney (C = 8)0.0480.00180.7790.101
Breast (C = 7)0.5970.03430.864 × 10−81.06 × 10−34
Colon (C = 7)0.04650.6260.1343.66 × 10−4

Comparison results between SKF and single kernel on three datasets.

3.4. Comparing With Other Fusion Methods

Several multiple kernel fusion strategies have been developed, including similarity network fusion (SNF) (Wang et al., 2014) and unsupervised multiple kernel learning (UMKL) (Mariette and Villavialaneix, 2018). We compared the performance of SKF with these two strategies to find better subtypes for a particular cancer. We tested the three strategies on the three datasets to compare the performance of different fusion methods. All results are found in the Supplementary Table 1. The graphical results are shown in Figure 3, with the X axis representing the number of clusters and the Y axis representing the value of −log10(Pvalue). The blue lines represent the change of SKF, the red lines represent the change of SNF, the green lines represent the change of UMKL and the black dashed lines show the P-value equal to 0.05. In Figure 3, we find that SKF achieved a remarkable level of performance for the clustering of breast and colon cancer subtypes in the three datasets. Additionally, SKF achieved better performance than other kernel fusion strategies for the clustering of lung cancer subtypes in Datasets No.1 and No.3. We also found that SNF performed well for the clustering of kidney cancer subtypes in the three datasets and UMKL reached the best level of performance for the clustering of lung cancer subtypes in Dataset No.2 and stomach cancer subtypes in Dataset No.1. It demonstrates that SKF obtained a significant level performance for discovering subtypes of a particular cancer, and also that the cluster results can be used for guiding clinical treatment.

Figure 3

3.5. Survival Analysis

In this paper, we analyzed the performance of SKF based on six cancers, including breast, lung, kidney, colon, stomach, and GBM cancers. However, since the P-values for the clustering of kidney and GBM cancer subtypes were larger than 0.05, we showed survival probability curves for the four other cancers. We analyzed these cancer subtypes by using Dataset No.3 . In Figure 4, we find that subtype 3 for stomach cancer has a higher death rate. These patients with subtype 3 need more attention to be paid to them. The average survival time of subtype 2 for colon cancer is longer than the other subtypes. Similarly, subtype 3 for other cancers tends to be more aggressive than other subtypes. We also found that the average survival time for breast cancer and lung cancer are longer than for stomach and colon cancer. It demonstrates that the cluster results of SKF can be used to guide clinical treatment.

Figure 4

4. Conclusions

In this paper, we proposed an accurate model for predicting cancer subtypes. First, we extracted a novel dataset with three expression data types (gene expression, miRNA expression, and isoform expression) and five cancers (breast, lung, kidney, colon, and stomach cancers) from the TCGA website. Second, we constructed three similarity kernels by using the three types of expression data for each cancer. Then, we proposed Similarity Kernel Fusion (SKF) to integrate the three kernels into one combined kernel. Finally, we used spectral clustering on integrated kernel to discover cancer subtypes.

We used an evaluation criteria (P-value) and a verification method (survival analysis) to evaluate the performance of SKF for the discovery of cancer subtypes. We compared SKF with single kernel and two kernel fusion strategies (SNF and UMKL) in three datasets. Results showed that SKF obtains a significant level of performance on P-value, and the survival curve of the subtypes was consistent with the clinical data. It demonstrates that SKF is an accurate computational tool for guiding clinical treatment.

Our method also has some limitations that require some attention. Since spectral clustering is a widely used and accepted cluster method, we are attempting to find an improved method to discover cancer subtypes more accurately. We will consider various machine learning methods and constructing kernel methods to predict cancer subtypes (Zeng et al., 2017; Ding et al., 2018; Zhang et al., 2018a,b,c; Zou et al., 2018). We also consider the potential possibility of developing computational models for cancer subtype identification based on microRNA information (Chen and Huang, 2017; Chen et al., 2017, 2018a,b; Hu et al., 2018).

Statements

Data availability statement

The results and codes for this study can be found at the following address: https://github.com/guofei-tju/Cancer-subtypes.

Author contributions

FG and LJ conceived and designed the experiments. LJ and YX performed the experiments and analyzed the data. FG and YX wrote the paper. FG, YD, and JT supervised the experiments and reviewed the manuscript.

Funding

This work is supported by a grant from the National Natural Science Foundation of China (NSFC 61772362), the Tianjin Research Program of Application Foundation and Advanced Technology (16JCQNJC00200) and the National Key R&D Program of China (2018YFC0910405, 2017YFC0908400).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00020/full#supplementary-material

References

  • 1

    BrunetJ. P.TamayoP.GolubT. R.MesirovJ. P. (2004). Metagenes and molecular pattern discovery using matrix factorization. Proc. Natl. Acad. Sci. U.S.A.101, 4164–4169. 10.1073/pnas.0308531101

  • 2

    ChenX.HuangL. (2017). LRSSLMDA: laplacian regularized sparse subspace learning for miRNA-disease association prediction. PLoS Comput. Biol.13:e1005912. 10.1371/journal.pcbi.1005912

  • 3

    ChenX.HuangL.XieD.ZhaoQ. (2018a). EGBMMDA: extreme gradient boosting machine for miRNA-disease association prediction. Cell Death Dis.9:3. 10.1038/s41419-017-0003-x

  • 4

    ChenX.XieD.WangL.ZhaoQ.YouZ.LiuH. (2018b). BNPMDA: bipartite network projection for miRNA–disease association prediction. Bioinformatics34, 3178–3186. 10.1093/bioinformatics/bty333

  • 5

    ChenX.XieD.ZhaoQ.YouZ. H. (2017). MicroRNAs and complex diseases: from experimental results to computational models. Brief. Bioinform.18:558. 10.1093/bib/bbx130

  • 6

    DingY.TangJ.GuoF. (2017a). Identification of drug-target interactions via multiple information integration. Inform. Sci.418, 546–560. 10.1016/j.ins.2017.08.045

  • 7

    DingY.TangJ.GuoF. (2017b). Identification of protein-ligand binding sites by sequence information and ensemble classifier. J. Chem. Inform. Model.57, 3149–3161. 10.1021/acs.jcim.7b00307

  • 8

    DingY.TangJ.GuoF. (2018). Identification of drug-side effect association via multiple information integration with centered kernel alignment. Neurocomputing. 10.1016/j.neucom.2018.10.028

  • 9

    FinneganT. J.CareyL. A. (2007). Gene-expression analysis and the basal-like breast cancer subtype. Future Oncol.3, 55–63. 10.2217/14796694.3.1.55

  • 10

    GaoY.ChurchG. (2005). Improving molecular cancer class discovery through sparse non-negative matrix factorization. Bioinformatics21:3970–3975. 10.1093/bioinformatics/bti653

  • 11

    GeS. G.XiaJ.ShaW.ZhengC. H. (2017). Cancer subtype discovery based on integrative model of multigenomic data. IEEE ACM Trans. Comput. Biol. Bioinform.14, 1115–1121. 10.1109/TCBB.2016.2621769

  • 12

    GuoF.WangD.WangL. (2018). Progressive approach for snp calling and haplotype assembly using single molecular sequencing data. Bioinformatics34, 2012–2018. 10.1093/bioinformatics/bty059

  • 13

    GuoY.QiY.LiZ.ShangX. (2018). Improvement of cancer subtype prediction by incorporating transcriptome expression data and heterogeneous biological networks, in Genome Informatics Workshop, (Kunming).

  • 14

    HuH.ZhangL.AiH.ZhangH.FanY.ZhaoQ.et al. (2018). HLPI-ensemble: prediction of human lncRNA-protein interactions based on ensemble strategy. RNA Biol.15:1. 10.1080/15476286.2018.1457935

  • 15

    KatarzynaT.PatrycjaC.MaciejW. (2015). The cancer genome atlas (TCGA): an immeasurable source of knowledge. Contemp. Oncol.19, 68–77. 10.5114/wo.2014.47136

  • 16

    KruijfE. M.EngelsC. C.van de WaterW.BastiaannetE.SmitV. T.van de VeldeC. J.et al. (2013). Tumor immune subtypes distinguish tumor subclasses with clinical implications in breast cancer patients. Breast Cancer Res. Treat.142, 355–364. 10.1007/s10549-013-2752-2

  • 17

    Le VanT.van LeeuwenM.Carolina FierroA.De MaeyerD.Van den EyndenJ.VerbekeL.et al. (2016). Simultaneous discovery of cancer subtypes and subtype features by molecular data integration. Bioinformatics32:i445. 10.1093/bioinformatics/btw434

  • 18

    LiuX.ChuK. M. (2014). E-cadherin and gastric cancer: cause, consequence, and applications. Biomed. Res. Int.2014:637308. 10.1155/2014/637308

  • 19

    LuxburgU. V. (2007). A tutorial on spectral clustering. Stat. Comput.17, 395–416. 10.1007/s11222-007-9033-z

  • 20

    MaT.ZhangA. (2017). Integrate multi-omic data using affinity network fusion (anf) for cancer patient clustering, in IEEE International Conference on Bioinformatics and Biomedicine, (Kansas City, MO) 398–403.

  • 21

    MagerD. L. (2006). Bacteria and cancer: cause, coincidence or cure. J. Trans. Med.4, 1–18. 10.1186/1479-5876-4-14

  • 22

    MarietteJ.VillavialaneixN. (2018). Unsupervised multiple kernel learning for heterogeneous data integration. Bioinformatics34, 1009–1015. 10.1093/bioinformatics/btx682

  • 23

    NgA. Y.JordanM. I.WeissY. (2001). On spectral clustering: analysis and an algorithm, in International Conference on Neural Information Processing Systems: Natural and Synthetic, (Vancouver, BC) 849–856.

  • 24

    NguyenD. V.RockeD. M. (2002). Multi-class cancer classification via partial least squares with gene expression profiles. Bioinformatics18, 1216–1226. 10.1093/bioinformatics/18.9.1216

  • 25

    PanG.JiangL.TangJ.GuoF. (2018). A novel computational method for detecting DNA methylation sites with DNA sequence information and physicochemical properties. Int. J. Mol. Sci.19:511. 10.3390/ijms19020511

  • 26

    PanG.TangJ.GuoF. (2017). Analysis of co-associated transcription factors via ordered adjacency differences on motif distribution. Sci. Rep.7:43597. 10.1038/srep43597

  • 27

    PratA.PinedaE.AdamoB.GalvnP.FernndezA.GabaL.et al. (2015). Clinical implications of the intrinsic molecular subtypes of breast cancer. Breast24(Suppl. 2), S26–S35. 10.1016/j.breast.2015.07.008

  • 28

    RenX.FuH.JinQ. (2015). Integrating heterogeneous genomic data to accurately identify disease subtypes. BMC Med. Genomics8:78. 10.1186/s12920-015-0154-5

  • 29

    ShenR.OlshenA. B.LadanyiM. (2009). Integrative clustering of multiple genomic data types using a joint latent variable model with application to breast and lung cancer subtype analysis. Bioinformatics25:2906. 10.1093/bioinformatics/btp543

  • 30

    SongJ.TangJ.GuoF. (2018). Identification of inhibitors of mmps enzymes via a novel computational approach:. Int. J. Biol. Sci.14, 863–871. 10.7150/ijbs.24588

  • 31

    SpeicherN. K.PfeiferN. (2015). Integrating different data types by regularized unsupervised multiple kernel learning with application to cancer subtype discovery. Bioinformatics31, i268–i275. 10.1093/bioinformatics/btv244

  • 32

    TeschendorffA. E.MiremadiA.PinderS. E.EllisI. O.CaldasC. (2007). An immune response gene expression module identifies a good prognosis subtype in estrogen receptor negative breast cancer. Genome Biol.8:R157. 10.1186/gb-2007-8-8-r157

  • 33

    ThankiK.NichollsM. E.GajjarA.SenagoreA. J.QiuS.SzaboC.et al. (2017). Consensus molecular subtypes of colorectal cancer and their clinical implications. Int. Biol. Biomed. J.3, 105–111.

  • 34

    WangB.MezliniA. M.DemirF.FiumeM.TuZ.BrudnoM.et al. (2014). Similarity network fusion for aggregating data types on a genomic scale. Nat. Methods11, 333–337. 10.1038/nmeth.2810

  • 35

    WeiL.LuanS.NagaiL. A. E.SuR.ZouQ. (2018a). Exploring sequence-based features for the improved prediction of DNA n4-methylcytosine sites in multiple species. Bioinformatics. [Epub ahead of print]. 10.1093/bioinformatics/bty824

  • 36

    WeiL.XingP.ShiG.JiZ.-L.ZouQ. (2017). Fast prediction of protein methylation sites using a sequence-based feature selection technique. IEEE/ACM Trans. Comput. Biol. Bioinform.34, 4007–401610.1109/TCBB.2017.2670558

  • 37

    WeiL.ZhouC.ChenH.SongJ.SuR. (2018b). ACPred-FL: a sequence-based predictor based on effective feature representation to improve the prediction of anti-cancer peptides. Bioinformatics34, 4007–4016. 10.1093/bioinformatics/bty451

  • 38

    WongG.LeckieC.KowalczykA. (2012). FSR: feature set reduction for scalable and accurate multi-class cancer subtype classification based on copy number. Bioinformatics28:151. 10.1093/bioinformatics/btr644

  • 39

    XuT.LeT. D.LiuL.SuN.WangR.SunB.et al. (2017). Cancersubtypes: an r/bioconductor package for molecular cancer subtype identification, validation and visualization. Bioinformatics33, 3131–3133. 10.1093/bioinformatics/btx378

  • 40

    XuT.LeT. D.LiuL.WangR.SunB.LiJ. (2016). Identifying cancer subtypes from miRNA-tf-mRNA regulatory networks and expression data. PLoS ONE11:e0152792. 10.1371/journal.pone.0152792

  • 41

    ZengX.DingN.Rodrguez-PatnA.QuanZ. (2017). Probability-based collaborative filtering model for predicting gene–disease associations. BMC Med. Genomics10(Suppl. 5):76. 10.1186/s12920-017-0313-y

  • 42

    ZengX.LiaoY.LiuY.ZouQ. (2016). Prediction and validation of disease genes using hetesim scores. IEEE ACM Trans. Comput. Biol. Bioinform.14, 687–695. 10.1109/TCBB.2016.2520947

  • 43

    ZhangW.FengH.WuH.ZhengX. (2017). Accounting for tumor purity improves cancer subtype classification from DNA methylation data. Bioinformatics33:2651. 10.1093/bioinformatics/btx303

  • 44

    ZhangW.LiuX.ChenY.WuW.WangW.LiX. (2018a). Feature-derived graph regularized matrix factorization for predicting drug side effects. Neurocomputing287, 154–162. 10.1016/j.neucom.2018.01.085

  • 45

    ZhangW.YangW.LuX.HuangF.LuoF. (2018b). The bi-direction similarity integration method for predicting microbe-disease associations. IEEE Access6, 38052–38061. 10.1109/ACCESS.2018.2851751

  • 46

    ZhangW.YueX.LinW.WuW.LiuR.HuangF.LiuF. (2018c). Predicting drug-disease associations by using similarity constrained matrix factorization. BMC Bioinform.19:233. 10.1186/s12859-018-2220-4

  • 47

    ZhaoQ.YuH.MingZ.HuH.RenG.LiuH. (2018a). The bipartite network projection-recommended algorithm for predicting long non-coding RNA-protein interactions. Mol. Ther. Nucleic Acids13, 464–471. 10.1016/j.omtn.2018.09.020

  • 48

    ZhaoQ.ZhangY.HuH.RenG.ZhangW.LiuH. (2018b). Irwnrlpi: integrating random walk and neighborhood regularized logistic matrix factorization for lncRNA-protein interaction prediction. Front. Genet.9:239. 10.3389/fgene.2018.00239

  • 49

    ZouQ.LinG.JiangX.LiuX.ZengX. (2018). Sequence clustering in bioinformatics: an empirical study. Brief. Bioinform. [Epub ahead of print]. 10.1093/bib/bby090

Summary

Keywords

cancer subtypes prediction, similarity kernel fusion, spectral clustering, sparse matrix, The Cancer Genome Atlas

Citation

Jiang L, Xiao Y, Ding Y, Tang J and Guo F (2019) Discovering Cancer Subtypes via an Accurate Fusion Strategy on Multiple Profile Data. Front. Genet. 10:20. doi: 10.3389/fgene.2019.00020

Received

30 October 2018

Accepted

15 January 2019

Published

05 February 2019

Volume

10 - 2019

Edited by

Arun Kumar Sangaiah, VIT University, India

Reviewed by

Qi Zhao, Liaoning University, China; Xianwen Ren, Peking University, China

Updates

Copyright

*Correspondence: Jijun Tang Fei Guo

This article was submitted to Bioinformatics and Computational Biology, a section of the journal Frontiers in Genetics

†These authors have contributed equally to this work

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics