This pipeline computes the correlation between cancer subtypes identified by different molecular patterns and selected clinical features.
Testing the association between subtypes identified by 8 different clustering approaches and 6 clinical features across 34 patients, 3 significant findings detected with P value < 0.05 and Q value < 0.25.
-
5 subtypes identified in current cancer cohort by 'Copy Number Ratio CNMF subtypes'. These subtypes correlate to 'PATHOLOGY.N.STAGE'.
-
4 subtypes identified in current cancer cohort by 'METHLYATION CNMF'. These subtypes correlate to 'GENDER'.
-
CNMF clustering analysis on sequencing-based mRNA expression data identified 4 subtypes that do not correlate to any clinical features.
-
Consensus hierarchical clustering analysis on sequencing-based mRNA expression data identified 3 subtypes that correlate to 'GENDER'.
-
3 subtypes identified in current cancer cohort by 'MIRSEQ CNMF'. These subtypes do not correlate to any clinical features.
-
3 subtypes identified in current cancer cohort by 'MIRSEQ CHIERARCHICAL'. These subtypes do not correlate to any clinical features.
-
3 subtypes identified in current cancer cohort by 'MIRseq Mature CNMF subtypes'. These subtypes do not correlate to any clinical features.
-
3 subtypes identified in current cancer cohort by 'MIRseq Mature cHierClus subtypes'. These subtypes do not correlate to any clinical features.
Clinical Features |
Time to Death |
AGE |
NEOPLASM DISEASESTAGE |
PATHOLOGY T STAGE |
PATHOLOGY N STAGE |
GENDER |
Statistical Tests | logrank test | ANOVA | Chi-square test | Chi-square test | Fisher's exact test | Fisher's exact test |
Copy Number Ratio CNMF subtypes |
0.0844 (1.00) |
0.196 (1.00) |
0.303 (1.00) |
0.721 (1.00) |
0.000174 (0.00835) |
0.263 (1.00) |
METHLYATION CNMF |
0.0281 (1.00) |
0.806 (1.00) |
0.592 (1.00) |
0.742 (1.00) |
0.18 (1.00) |
0.000689 (0.0324) |
RNAseq CNMF subtypes |
0.00759 (0.342) |
0.293 (1.00) |
0.468 (1.00) |
0.34 (1.00) |
0.0951 (1.00) |
0.042 (1.00) |
RNAseq cHierClus subtypes |
0.0194 (0.816) |
0.26 (1.00) |
0.166 (1.00) |
0.26 (1.00) |
0.014 (0.604) |
0.00391 (0.18) |
MIRSEQ CNMF |
0.0493 (1.00) |
0.273 (1.00) |
0.108 (1.00) |
0.0393 (1.00) |
0.122 (1.00) |
0.021 (0.861) |
MIRSEQ CHIERARCHICAL |
0.536 (1.00) |
0.0754 (1.00) |
0.0454 (1.00) |
0.0848 (1.00) |
0.327 (1.00) |
0.151 (1.00) |
MIRseq Mature CNMF subtypes |
0.119 (1.00) |
0.399 (1.00) |
0.0721 (1.00) |
0.0592 (1.00) |
0.223 (1.00) |
0.0354 (1.00) |
MIRseq Mature cHierClus subtypes |
0.536 (1.00) |
0.0137 (0.602) |
0.0454 (1.00) |
0.0848 (1.00) |
0.327 (1.00) |
0.355 (1.00) |
Cluster Labels | 1 | 2 | 3 | 4 | 5 |
---|---|---|---|---|---|
Number of samples | 10 | 9 | 3 | 8 | 3 |
P value = 0.0844 (logrank test), Q value = 1
nPatients | nDeath | Duration Range (Median), Month | |
---|---|---|---|
ALL | 33 | 8 | 6.9 - 121.2 (29.2) |
subtype1 | 10 | 4 | 8.3 - 50.7 (23.9) |
subtype2 | 9 | 0 | 6.9 - 121.2 (41.3) |
subtype3 | 3 | 1 | 10.2 - 23.3 (17.8) |
subtype4 | 8 | 3 | 11.3 - 63.2 (28.8) |
subtype5 | 3 | 0 | 22.2 - 73.0 (46.8) |
P value = 0.196 (ANOVA), Q value = 1
nPatients | Mean (Std.Dev) | |
---|---|---|
ALL | 33 | 51.1 (14.5) |
subtype1 | 10 | 46.9 (17.6) |
subtype2 | 9 | 55.9 (8.7) |
subtype3 | 3 | 59.7 (6.1) |
subtype4 | 8 | 53.1 (11.9) |
subtype5 | 3 | 36.7 (22.0) |
P value = 0.303 (Chi-square test), Q value = 1
nPatients | STAGE I | STAGE II | STAGE III | STAGE IV |
---|---|---|---|---|
ALL | 4 | 13 | 4 | 8 |
subtype1 | 1 | 3 | 2 | 3 |
subtype2 | 1 | 5 | 2 | 1 |
subtype3 | 0 | 0 | 0 | 3 |
subtype4 | 1 | 3 | 0 | 1 |
subtype5 | 1 | 2 | 0 | 0 |
P value = 0.721 (Chi-square test), Q value = 1
nPatients | T1 | T2+T3 | T4 |
---|---|---|---|
ALL | 4 | 16 | 9 |
subtype1 | 1 | 4 | 4 |
subtype2 | 1 | 6 | 2 |
subtype3 | 0 | 1 | 2 |
subtype4 | 1 | 3 | 1 |
subtype5 | 1 | 2 | 0 |
P value = 0.000174 (Chi-square test), Q value = 0.0084
nPatients | 0 | 1 |
---|---|---|
ALL | 26 | 4 |
subtype1 | 8 | 1 |
subtype2 | 9 | 0 |
subtype3 | 0 | 3 |
subtype4 | 6 | 0 |
subtype5 | 3 | 0 |
P value = 0.263 (Chi-square test), Q value = 1
nPatients | FEMALE | MALE |
---|---|---|
ALL | 17 | 16 |
subtype1 | 8 | 2 |
subtype2 | 3 | 6 |
subtype3 | 1 | 2 |
subtype4 | 4 | 4 |
subtype5 | 1 | 2 |
Cluster Labels | 1 | 2 | 3 | 4 | 5 |
---|---|---|---|---|---|
Number of samples | 6 | 2 | 12 | 11 | 3 |
P value = 0.0281 (logrank test), Q value = 1
nPatients | nDeath | Duration Range (Median), Month | |
---|---|---|---|
ALL | 32 | 8 | 6.9 - 121.2 (28.3) |
subtype1 | 6 | 4 | 8.3 - 33.8 (18.3) |
subtype3 | 12 | 3 | 10.2 - 106.5 (36.9) |
subtype4 | 11 | 1 | 6.9 - 121.2 (29.2) |
subtype5 | 3 | 0 | 12.6 - 22.2 (22.0) |
P value = 0.806 (ANOVA), Q value = 1
nPatients | Mean (Std.Dev) | |
---|---|---|
ALL | 32 | 50.1 (14.4) |
subtype1 | 6 | 52.3 (19.0) |
subtype3 | 12 | 50.8 (12.5) |
subtype4 | 11 | 50.3 (13.5) |
subtype5 | 3 | 42.3 (20.3) |
P value = 0.592 (Chi-square test), Q value = 1
nPatients | STAGE I | STAGE II | STAGE III | STAGE IV |
---|---|---|---|---|
ALL | 4 | 11 | 4 | 8 |
subtype1 | 0 | 2 | 2 | 2 |
subtype3 | 1 | 2 | 1 | 4 |
subtype4 | 2 | 5 | 1 | 2 |
subtype5 | 1 | 2 | 0 | 0 |
P value = 0.742 (Chi-square test), Q value = 1
nPatients | T1 | T2+T3 | T4 |
---|---|---|---|
ALL | 4 | 14 | 9 |
subtype1 | 0 | 3 | 3 |
subtype3 | 1 | 4 | 3 |
subtype4 | 2 | 5 | 3 |
subtype5 | 1 | 2 | 0 |
P value = 0.18 (Fisher's exact test), Q value = 1
nPatients | 0 | 1 |
---|---|---|
ALL | 24 | 4 |
subtype1 | 5 | 1 |
subtype3 | 6 | 3 |
subtype4 | 10 | 0 |
subtype5 | 3 | 0 |
P value = 0.000689 (Fisher's exact test), Q value = 0.032
nPatients | FEMALE | MALE |
---|---|---|
ALL | 15 | 17 |
subtype1 | 6 | 0 |
subtype3 | 2 | 10 |
subtype4 | 4 | 7 |
subtype5 | 3 | 0 |
Cluster Labels | 1 | 2 | 3 | 4 |
---|---|---|---|---|
Number of samples | 9 | 4 | 10 | 11 |
P value = 0.00759 (logrank test), Q value = 0.34
nPatients | nDeath | Duration Range (Median), Month | |
---|---|---|---|
ALL | 34 | 8 | 6.9 - 121.2 (29.8) |
subtype1 | 9 | 5 | 8.3 - 37.1 (18.1) |
subtype2 | 4 | 0 | 10.2 - 41.4 (27.3) |
subtype3 | 10 | 3 | 17.4 - 106.5 (46.6) |
subtype4 | 11 | 0 | 6.9 - 121.2 (29.2) |
P value = 0.293 (ANOVA), Q value = 1
nPatients | Mean (Std.Dev) | |
---|---|---|
ALL | 34 | 50.8 (14.3) |
subtype1 | 9 | 54.1 (15.6) |
subtype2 | 4 | 59.2 (9.2) |
subtype3 | 10 | 51.1 (10.3) |
subtype4 | 11 | 44.8 (16.8) |
P value = 0.468 (Chi-square test), Q value = 1
nPatients | STAGE I | STAGE II | STAGE III | STAGE IV |
---|---|---|---|---|
ALL | 4 | 13 | 4 | 8 |
subtype1 | 0 | 4 | 2 | 3 |
subtype2 | 0 | 1 | 1 | 2 |
subtype3 | 1 | 3 | 1 | 2 |
subtype4 | 3 | 5 | 0 | 1 |
P value = 0.34 (Chi-square test), Q value = 1
nPatients | T1 | T2+T3 | T4 |
---|---|---|---|
ALL | 4 | 16 | 9 |
subtype1 | 0 | 5 | 4 |
subtype2 | 0 | 3 | 1 |
subtype3 | 1 | 3 | 3 |
subtype4 | 3 | 5 | 1 |
P value = 0.0951 (Fisher's exact test), Q value = 1
nPatients | 0 | 1 |
---|---|---|
ALL | 26 | 4 |
subtype1 | 8 | 1 |
subtype2 | 2 | 2 |
subtype3 | 6 | 1 |
subtype4 | 10 | 0 |
P value = 0.042 (Fisher's exact test), Q value = 1
nPatients | FEMALE | MALE |
---|---|---|
ALL | 17 | 17 |
subtype1 | 8 | 1 |
subtype2 | 2 | 2 |
subtype3 | 3 | 7 |
subtype4 | 4 | 7 |
Cluster Labels | 1 | 2 | 3 |
---|---|---|---|
Number of samples | 8 | 17 | 9 |
P value = 0.0194 (logrank test), Q value = 0.82
nPatients | nDeath | Duration Range (Median), Month | |
---|---|---|---|
ALL | 34 | 8 | 6.9 - 121.2 (29.8) |
subtype1 | 8 | 4 | 8.3 - 37.1 (18.3) |
subtype2 | 17 | 1 | 6.9 - 121.2 (29.2) |
subtype3 | 9 | 3 | 10.2 - 63.2 (31.2) |
P value = 0.26 (ANOVA), Q value = 1
nPatients | Mean (Std.Dev) | |
---|---|---|
ALL | 34 | 50.8 (14.3) |
subtype1 | 8 | 54.4 (16.7) |
subtype2 | 17 | 46.8 (14.7) |
subtype3 | 9 | 55.3 (10.1) |
P value = 0.166 (Chi-square test), Q value = 1
nPatients | STAGE I | STAGE II | STAGE III | STAGE IV |
---|---|---|---|---|
ALL | 4 | 13 | 4 | 8 |
subtype1 | 0 | 4 | 2 | 2 |
subtype2 | 4 | 7 | 1 | 2 |
subtype3 | 0 | 2 | 1 | 4 |
P value = 0.26 (Chi-square test), Q value = 1
nPatients | T1 | T2+T3 | T4 |
---|---|---|---|
ALL | 4 | 16 | 9 |
subtype1 | 0 | 5 | 3 |
subtype2 | 4 | 7 | 3 |
subtype3 | 0 | 4 | 3 |
P value = 0.014 (Fisher's exact test), Q value = 0.6
nPatients | 0 | 1 |
---|---|---|
ALL | 26 | 4 |
subtype1 | 7 | 1 |
subtype2 | 15 | 0 |
subtype3 | 4 | 3 |
P value = 0.00391 (Fisher's exact test), Q value = 0.18
nPatients | FEMALE | MALE |
---|---|---|
ALL | 17 | 17 |
subtype1 | 8 | 0 |
subtype2 | 6 | 11 |
subtype3 | 3 | 6 |
Cluster Labels | 1 | 2 | 3 |
---|---|---|---|
Number of samples | 13 | 7 | 14 |
P value = 0.0493 (logrank test), Q value = 1
nPatients | nDeath | Duration Range (Median), Month | |
---|---|---|---|
ALL | 34 | 8 | 6.9 - 121.2 (29.8) |
subtype1 | 13 | 4 | 10.1 - 106.5 (22.0) |
subtype2 | 7 | 3 | 8.3 - 42.5 (18.5) |
subtype3 | 14 | 1 | 6.9 - 121.2 (48.7) |
P value = 0.273 (ANOVA), Q value = 1
nPatients | Mean (Std.Dev) | |
---|---|---|
ALL | 34 | 50.8 (14.3) |
subtype1 | 13 | 51.1 (13.0) |
subtype2 | 7 | 57.9 (15.8) |
subtype3 | 14 | 47.1 (14.4) |
P value = 0.108 (Chi-square test), Q value = 1
nPatients | STAGE I | STAGE II | STAGE III | STAGE IV |
---|---|---|---|---|
ALL | 4 | 13 | 4 | 8 |
subtype1 | 1 | 3 | 2 | 5 |
subtype2 | 0 | 2 | 1 | 3 |
subtype3 | 3 | 8 | 1 | 0 |
P value = 0.0393 (Chi-square test), Q value = 1
nPatients | T1 | T2+T3 | T4 |
---|---|---|---|
ALL | 4 | 16 | 9 |
subtype1 | 1 | 4 | 6 |
subtype2 | 0 | 3 | 3 |
subtype3 | 3 | 9 | 0 |
P value = 0.122 (Fisher's exact test), Q value = 1
nPatients | 0 | 1 |
---|---|---|
ALL | 26 | 4 |
subtype1 | 8 | 3 |
subtype2 | 5 | 1 |
subtype3 | 13 | 0 |
P value = 0.021 (Fisher's exact test), Q value = 0.86
nPatients | FEMALE | MALE |
---|---|---|
ALL | 17 | 17 |
subtype1 | 9 | 4 |
subtype2 | 5 | 2 |
subtype3 | 3 | 11 |
Cluster Labels | 1 | 2 | 3 |
---|---|---|---|
Number of samples | 14 | 5 | 15 |
P value = 0.536 (logrank test), Q value = 1
nPatients | nDeath | Duration Range (Median), Month | |
---|---|---|---|
ALL | 34 | 8 | 6.9 - 121.2 (29.8) |
subtype1 | 14 | 3 | 6.9 - 121.2 (30.6) |
subtype2 | 5 | 1 | 29.2 - 93.4 (46.8) |
subtype3 | 15 | 4 | 8.3 - 106.5 (22.0) |
P value = 0.0754 (ANOVA), Q value = 1
nPatients | Mean (Std.Dev) | |
---|---|---|
ALL | 34 | 50.8 (14.3) |
subtype1 | 14 | 53.9 (14.5) |
subtype2 | 5 | 37.6 (13.9) |
subtype3 | 15 | 52.4 (12.6) |
P value = 0.0454 (Chi-square test), Q value = 1
nPatients | STAGE I | STAGE II | STAGE III | STAGE IV |
---|---|---|---|---|
ALL | 4 | 13 | 4 | 8 |
subtype1 | 1 | 9 | 0 | 3 |
subtype2 | 2 | 1 | 1 | 0 |
subtype3 | 1 | 3 | 3 | 5 |
P value = 0.0848 (Chi-square test), Q value = 1
nPatients | T1 | T2+T3 | T4 |
---|---|---|---|
ALL | 4 | 16 | 9 |
subtype1 | 1 | 9 | 3 |
subtype2 | 2 | 2 | 0 |
subtype3 | 1 | 5 | 6 |
P value = 0.327 (Fisher's exact test), Q value = 1
nPatients | 0 | 1 |
---|---|---|
ALL | 26 | 4 |
subtype1 | 13 | 1 |
subtype2 | 4 | 0 |
subtype3 | 9 | 3 |
P value = 0.151 (Fisher's exact test), Q value = 1
nPatients | FEMALE | MALE |
---|---|---|
ALL | 17 | 17 |
subtype1 | 6 | 8 |
subtype2 | 1 | 4 |
subtype3 | 10 | 5 |
Cluster Labels | 1 | 2 | 3 |
---|---|---|---|
Number of samples | 14 | 7 | 13 |
P value = 0.119 (logrank test), Q value = 1
nPatients | nDeath | Duration Range (Median), Month | |
---|---|---|---|
ALL | 34 | 8 | 6.9 - 121.2 (29.8) |
subtype1 | 14 | 4 | 8.3 - 106.5 (19.9) |
subtype2 | 7 | 3 | 18.1 - 63.2 (33.8) |
subtype3 | 13 | 1 | 6.9 - 121.2 (46.8) |
P value = 0.399 (ANOVA), Q value = 1
nPatients | Mean (Std.Dev) | |
---|---|---|
ALL | 34 | 50.8 (14.3) |
subtype1 | 14 | 51.6 (12.7) |
subtype2 | 7 | 56.1 (16.3) |
subtype3 | 13 | 47.1 (15.0) |
P value = 0.0721 (Chi-square test), Q value = 1
nPatients | STAGE I | STAGE II | STAGE III | STAGE IV |
---|---|---|---|---|
ALL | 4 | 13 | 4 | 8 |
subtype1 | 1 | 3 | 3 | 5 |
subtype2 | 0 | 3 | 0 | 3 |
subtype3 | 3 | 7 | 1 | 0 |
P value = 0.0592 (Chi-square test), Q value = 1
nPatients | T1 | T2+T3 | T4 |
---|---|---|---|
ALL | 4 | 16 | 9 |
subtype1 | 1 | 5 | 6 |
subtype2 | 0 | 3 | 3 |
subtype3 | 3 | 8 | 0 |
P value = 0.223 (Fisher's exact test), Q value = 1
nPatients | 0 | 1 |
---|---|---|
ALL | 26 | 4 |
subtype1 | 9 | 3 |
subtype2 | 5 | 1 |
subtype3 | 12 | 0 |
P value = 0.0354 (Fisher's exact test), Q value = 1
nPatients | FEMALE | MALE |
---|---|---|
ALL | 17 | 17 |
subtype1 | 10 | 4 |
subtype2 | 4 | 3 |
subtype3 | 3 | 10 |
Cluster Labels | 1 | 2 | 3 |
---|---|---|---|
Number of samples | 14 | 6 | 14 |
P value = 0.536 (logrank test), Q value = 1
nPatients | nDeath | Duration Range (Median), Month | |
---|---|---|---|
ALL | 34 | 8 | 6.9 - 121.2 (29.8) |
subtype1 | 14 | 4 | 8.3 - 106.5 (22.7) |
subtype2 | 6 | 1 | 10.1 - 93.4 (39.0) |
subtype3 | 14 | 3 | 6.9 - 121.2 (30.6) |
P value = 0.0137 (ANOVA), Q value = 0.6
nPatients | Mean (Std.Dev) | |
---|---|---|
ALL | 34 | 50.8 (14.3) |
subtype1 | 14 | 54.2 (10.8) |
subtype2 | 6 | 35.8 (13.1) |
subtype3 | 14 | 53.9 (14.5) |
P value = 0.0454 (Chi-square test), Q value = 1
nPatients | STAGE I | STAGE II | STAGE III | STAGE IV |
---|---|---|---|---|
ALL | 4 | 13 | 4 | 8 |
subtype1 | 1 | 3 | 3 | 5 |
subtype2 | 2 | 1 | 1 | 0 |
subtype3 | 1 | 9 | 0 | 3 |
P value = 0.0848 (Chi-square test), Q value = 1
nPatients | T1 | T2+T3 | T4 |
---|---|---|---|
ALL | 4 | 16 | 9 |
subtype1 | 1 | 5 | 6 |
subtype2 | 2 | 2 | 0 |
subtype3 | 1 | 9 | 3 |
P value = 0.327 (Fisher's exact test), Q value = 1
nPatients | 0 | 1 |
---|---|---|
ALL | 26 | 4 |
subtype1 | 9 | 3 |
subtype2 | 4 | 0 |
subtype3 | 13 | 1 |
P value = 0.355 (Fisher's exact test), Q value = 1
nPatients | FEMALE | MALE |
---|---|---|
ALL | 17 | 17 |
subtype1 | 9 | 5 |
subtype2 | 2 | 4 |
subtype3 | 6 | 8 |
-
Cluster data file = ACC-TP.mergedcluster.txt
-
Clinical data file = ACC-TP.merged_data.txt
-
Number of patients = 34
-
Number of clustering approaches = 8
-
Number of selected clinical features = 6
-
Exclude small clusters that include fewer than K patients, K = 3
consensus non-negative matrix factorization clustering approach (Brunet et al. 2004)
Resampling-based clustering method (Monti et al. 2003)
For survival clinical features, the Kaplan-Meier survival curves of tumors with and without gene mutations were plotted and the statistical significance P values were estimated by logrank test (Bland and Altman 2004) using the 'survdiff' function in R
For continuous numerical clinical features, one-way analysis of variance (Howell 2002) was applied to compare the clinical values between tumor subtypes using 'anova' function in R
For multi-class clinical features (nominal or ordinal), Chi-square tests (Greenwood and Nikulin 1996) were used to estimate the P values using the 'chisq.test' function in R
For binary clinical features, two-tailed Fisher's exact tests (Fisher 1922) were used to estimate the P values using the 'fisher.test' function in R
For multiple hypothesis correction, Q value is the False Discovery Rate (FDR) analogue of the P value (Benjamini and Hochberg 1995), defined as the minimum FDR at which the test may be called significant. We used the 'Benjamini and Hochberg' method of 'p.adjust' function in R to convert P values into Q values.
In addition to the links below, the full results of the analysis summarized in this report can also be downloaded programmatically using firehose_get, or interactively from either the Broad GDAC website or TCGA Data Coordination Center Portal.