This pipeline computes the correlation between significantly recurrent gene mutations and selected clinical features.
Testing the association between mutation status of 29 genes and 10 clinical features across 417 patients, 2 significant findings detected with Q value < 0.25.
-
BAP1 mutation correlated to 'NEOPLASM.DISEASESTAGE' and 'PATHOLOGY.T.STAGE'.
Clinical Features |
Time to Death |
AGE |
NEOPLASM DISEASESTAGE |
PATHOLOGY T STAGE |
PATHOLOGY N STAGE |
PATHOLOGY M STAGE |
GENDER |
KARNOFSKY PERFORMANCE SCORE |
RACE | ETHNICITY | ||
nMutated (%) | nWild-Type | logrank test | Wilcoxon-test | Fisher's exact test | Fisher's exact test | Fisher's exact test | Fisher's exact test | Fisher's exact test | Wilcoxon-test | Fisher's exact test | Fisher's exact test | |
BAP1 | 41 (10%) | 376 |
0.00204 (0.529) |
0.597 (1.00) |
0.0002 (0.0526) |
0.00027 (0.0707) |
0.672 (1.00) |
0.0116 (1.00) |
0.0016 (0.417) |
0.682 (1.00) |
1 (1.00) |
|
SETD2 | 44 (11%) | 373 |
0.169 (1.00) |
0.285 (1.00) |
0.186 (1.00) |
0.107 (1.00) |
0.617 (1.00) |
0.198 (1.00) |
0.181 (1.00) |
0.834 (1.00) |
0.707 (1.00) |
|
PTEN | 18 (4%) | 399 |
0.358 (1.00) |
0.807 (1.00) |
0.0473 (1.00) |
0.245 (1.00) |
0.11 (1.00) |
0.329 (1.00) |
0.129 (1.00) |
0.714 (1.00) |
0.618 (1.00) |
1 (1.00) |
PBRM1 | 136 (33%) | 281 |
0.454 (1.00) |
0.473 (1.00) |
0.946 (1.00) |
0.206 (1.00) |
0.515 (1.00) |
1 (1.00) |
0.23 (1.00) |
0.909 (1.00) |
0.242 (1.00) |
0.231 (1.00) |
KDM5C | 26 (6%) | 391 |
0.0501 (1.00) |
0.0881 (1.00) |
0.988 (1.00) |
0.972 (1.00) |
0.584 (1.00) |
1 (1.00) |
0.00224 (0.578) |
0.756 (1.00) |
0.377 (1.00) |
|
VHL | 166 (40%) | 251 |
0.762 (1.00) |
0.777 (1.00) |
0.11 (1.00) |
0.154 (1.00) |
1 (1.00) |
0.222 (1.00) |
0.463 (1.00) |
0.283 (1.00) |
0.114 (1.00) |
0.818 (1.00) |
FAM200A | 5 (1%) | 412 |
0.481 (1.00) |
0.733 (1.00) |
0.751 (1.00) |
0.575 (1.00) |
1 (1.00) |
1 (1.00) |
0.053 (1.00) |
1 (1.00) |
1 (1.00) |
|
TP53 | 9 (2%) | 408 |
0.00125 (0.327) |
0.332 (1.00) |
0.137 (1.00) |
0.0575 (1.00) |
1 (1.00) |
0.162 (1.00) |
0.726 (1.00) |
1 (1.00) |
1 (1.00) |
|
NEFH | 5 (1%) | 412 |
0.685 (1.00) |
0.522 (1.00) |
0.115 (1.00) |
0.182 (1.00) |
1 (1.00) |
0.585 (1.00) |
0.661 (1.00) |
1 (1.00) |
1 (1.00) |
|
CCDC120 | 4 (1%) | 413 |
0.129 (1.00) |
0.9 (1.00) |
0.0789 (1.00) |
0.0911 (1.00) |
1 (1.00) |
1 (1.00) |
1 (1.00) |
1 (1.00) |
0.205 (1.00) |
|
MTOR | 23 (6%) | 394 |
0.116 (1.00) |
0.149 (1.00) |
0.271 (1.00) |
0.379 (1.00) |
0.172 (1.00) |
0.236 (1.00) |
0.378 (1.00) |
0.713 (1.00) |
1 (1.00) |
|
GUSB | 4 (1%) | 413 |
0.852 (1.00) |
0.0799 (1.00) |
0.205 (1.00) |
0.486 (1.00) |
1 (1.00) |
1 (1.00) |
1 (1.00) |
1 (1.00) |
1 (1.00) |
|
PIK3CA | 12 (3%) | 405 |
0.587 (1.00) |
0.898 (1.00) |
0.523 (1.00) |
0.545 (1.00) |
1 (1.00) |
1 (1.00) |
0.029 (1.00) |
1 (1.00) |
0.575 (1.00) |
|
PCK1 | 5 (1%) | 412 |
0.179 (1.00) |
0.322 (1.00) |
0.142 (1.00) |
0.253 (1.00) |
1 (1.00) |
1 (1.00) |
1 (1.00) |
1 (1.00) |
1 (1.00) |
|
ARID1A | 12 (3%) | 405 |
0.34 (1.00) |
0.343 (1.00) |
0.0846 (1.00) |
0.151 (1.00) |
1 (1.00) |
0.7 (1.00) |
0.122 (1.00) |
1 (1.00) |
0.162 (1.00) |
|
ATM | 11 (3%) | 406 |
0.79 (1.00) |
0.51 (1.00) |
0.149 (1.00) |
0.0479 (1.00) |
0.39 (1.00) |
0.692 (1.00) |
0.527 (1.00) |
1 (1.00) |
0.37 (1.00) |
|
DPCR1 | 6 (1%) | 411 |
0.407 (1.00) |
0.541 (1.00) |
0.1 (1.00) |
1 (1.00) |
1 (1.00) |
0.0553 (1.00) |
1 (1.00) |
0.273 (1.00) |
1 (1.00) |
|
FGFR3 | 4 (1%) | 413 |
0.942 (1.00) |
0.621 (1.00) |
0.286 (1.00) |
0.381 (1.00) |
0.0588 (1.00) |
1 (1.00) |
0.302 (1.00) |
1 (1.00) |
1 (1.00) |
|
PTCH1 | 7 (2%) | 410 |
0.0607 (1.00) |
0.272 (1.00) |
0.38 (1.00) |
0.0843 (1.00) |
1 (1.00) |
1 (1.00) |
0.7 (1.00) |
1 (1.00) |
1 (1.00) |
|
RBMX | 4 (1%) | 413 |
0.731 (1.00) |
0.205 (1.00) |
1 (1.00) |
1 (1.00) |
0.114 (1.00) |
0.505 (1.00) |
0.302 (1.00) |
1 (1.00) |
1 (1.00) |
|
ARAP3 | 3 (1%) | 414 |
0.769 (1.00) |
0.521 (1.00) |
0.818 (1.00) |
0.732 (1.00) |
0.41 (1.00) |
0.281 (1.00) |
0.146 (1.00) |
1 (1.00) |
||
GPR172B | 4 (1%) | 413 |
0.56 (1.00) |
0.838 (1.00) |
0.882 (1.00) |
0.792 (1.00) |
0.167 (1.00) |
1 (1.00) |
0.126 (1.00) |
1 (1.00) |
0.0734 (1.00) |
|
EGFR | 7 (2%) | 410 |
0.0201 (1.00) |
0.554 (1.00) |
0.0187 (1.00) |
0.00594 (1.00) |
0.217 (1.00) |
0.313 (1.00) |
0.43 (1.00) |
1 (1.00) |
1 (1.00) |
|
IL1RAP | 4 (1%) | 413 |
0.0718 (1.00) |
0.0629 (1.00) |
0.531 (1.00) |
0.381 (1.00) |
1 (1.00) |
0.505 (1.00) |
1 (1.00) |
0.188 (1.00) |
1 (1.00) |
|
BCL6 | 5 (1%) | 412 |
0.638 (1.00) |
0.442 (1.00) |
0.686 (1.00) |
0.697 (1.00) |
1 (1.00) |
0.585 (1.00) |
1 (1.00) |
1 (1.00) |
1 (1.00) |
|
NFAT5 | 6 (1%) | 411 |
0.856 (1.00) |
0.678 (1.00) |
0.707 (1.00) |
0.748 (1.00) |
1 (1.00) |
0.595 (1.00) |
0.67 (1.00) |
1 (1.00) |
1 (1.00) |
|
LARP1 | 5 (1%) | 412 |
0.396 (1.00) |
0.278 (1.00) |
0.391 (1.00) |
0.471 (1.00) |
1 (1.00) |
0.585 (1.00) |
0.167 (1.00) |
1 (1.00) |
1 (1.00) |
|
KIAA1751 | 6 (1%) | 411 |
0.869 (1.00) |
0.208 (1.00) |
0.573 (1.00) |
1 (1.00) |
0.264 (1.00) |
0.595 (1.00) |
0.67 (1.00) |
1 (1.00) |
1 (1.00) |
|
INSRR | 4 (1%) | 413 |
0.0254 (1.00) |
0.112 (1.00) |
1 (1.00) |
0.793 (1.00) |
1 (1.00) |
0.505 (1.00) |
0.614 (1.00) |
1 (1.00) |
1 (1.00) |
P value = 2e-04 (Fisher's exact test), Q value = 0.053
nPatients | STAGE I | STAGE II | STAGE III | STAGE IV |
---|---|---|---|---|
ALL | 197 | 40 | 113 | 67 |
BAP1 MUTATED | 7 | 7 | 15 | 12 |
BAP1 WILD-TYPE | 190 | 33 | 98 | 55 |
P value = 0.00027 (Fisher's exact test), Q value = 0.071
nPatients | T1 | T2 | T3 | T4 |
---|---|---|---|---|
ALL | 202 | 49 | 160 | 6 |
BAP1 MUTATED | 8 | 10 | 23 | 0 |
BAP1 WILD-TYPE | 194 | 39 | 137 | 6 |
-
Mutation data file = transformed.cor.cli.txt
-
Clinical data file = KIRC-TP.merged_data.txt
-
Number of patients = 417
-
Number of significantly mutated genes = 29
-
Number of selected clinical features = 10
-
Exclude genes that fewer than K tumors have mutations, K = 3
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 binary or multi-class clinical features (nominal or ordinal), 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.