Correlation between copy number variations of arm-level result and selected clinical features
Uterine Carcinosarcoma (Primary solid tumor)
21 August 2015  |  analyses__2015_08_21
Maintainer Information
Citation Information
Maintained by TCGA GDAC Team (Broad Institute/MD Anderson Cancer Center/Harvard Medical School)
Cite as Broad Institute TCGA Genome Data Analysis Center (2015): Correlation between copy number variations of arm-level result and selected clinical features. Broad Institute of MIT and Harvard. doi:10.7908/C1KP81HD
Overview
Introduction

This pipeline computes the correlation between significant arm-level copy number variations (cnvs) and selected clinical features.

Summary

Testing the association between copy number variation 79 arm-level events and 5 clinical features across 56 patients, no significant finding detected with Q value < 0.25.

  • No arm-level cnvs related to clinical features.

Results
Overview of the results

Table 1.  Get Full Table Overview of the association between significant copy number variation of 79 arm-level events and 5 clinical features. Shown in the table are P values (Q values). Thresholded by Q value < 0.25, no significant finding detected.

Clinical
Features
Time
to
Death
YEARS
TO
BIRTH
RADIATION
THERAPY
HISTOLOGICAL
TYPE
RACE
nCNV (%) nWild-Type logrank test Wilcoxon-test Fisher's exact test Fisher's exact test Fisher's exact test
1p gain 24 (43%) 32 0.739
(1.00)
0.169
(1.00)
1
(1.00)
0.451
(1.00)
0.269
(1.00)
1q gain 31 (55%) 25 0.202
(1.00)
0.509
(1.00)
0.787
(1.00)
0.697
(1.00)
0.149
(1.00)
2p gain 23 (41%) 33 0.0931
(1.00)
0.526
(1.00)
0.0998
(1.00)
0.69
(1.00)
1
(1.00)
2q gain 21 (38%) 35 0.689
(1.00)
0.635
(1.00)
0.0919
(1.00)
0.884
(1.00)
0.873
(1.00)
3p gain 12 (21%) 44 0.633
(1.00)
0.322
(1.00)
1
(1.00)
0.453
(1.00)
0.0943
(1.00)
3q gain 23 (41%) 33 0.655
(1.00)
0.309
(1.00)
1
(1.00)
0.833
(1.00)
0.678
(1.00)
4p gain 9 (16%) 47 0.457
(1.00)
0.738
(1.00)
0.718
(1.00)
0.524
(1.00)
0.163
(1.00)
4q gain 3 (5%) 53 0.935
(1.00)
0.985
(1.00)
0.0927
(1.00)
0.449
(1.00)
1
(1.00)
5p gain 23 (41%) 33 0.68
(1.00)
0.683
(1.00)
0.787
(1.00)
1
(1.00)
0.877
(1.00)
5q gain 8 (14%) 48 0.811
(1.00)
0.833
(1.00)
0.449
(1.00)
0.105
(1.00)
0.239
(1.00)
6p gain 30 (54%) 26 0.899
(1.00)
0.114
(1.00)
0.784
(1.00)
0.384
(1.00)
0.452
(1.00)
6q gain 27 (48%) 29 0.776
(1.00)
0.178
(1.00)
1
(1.00)
0.279
(1.00)
0.599
(1.00)
7p gain 21 (38%) 35 0.66
(1.00)
0.192
(1.00)
0.777
(1.00)
0.175
(1.00)
0.232
(1.00)
7q gain 17 (30%) 39 0.556
(1.00)
0.532
(1.00)
0.557
(1.00)
0.0875
(1.00)
0.867
(1.00)
8p gain 19 (34%) 37 0.0964
(1.00)
0.401
(1.00)
0.249
(1.00)
0.0335
(1.00)
0.306
(1.00)
8q gain 30 (54%) 26 0.129
(1.00)
0.336
(1.00)
0.0132
(1.00)
0.0911
(1.00)
0.45
(1.00)
9p gain 6 (11%) 50 0.953
(1.00)
0.353
(1.00)
1
(1.00)
0.0216
(1.00)
0.7
(1.00)
10p gain 21 (38%) 35 0.411
(1.00)
0.78
(1.00)
0.783
(1.00)
0.884
(1.00)
0.664
(1.00)
10q gain 17 (30%) 39 0.337
(1.00)
0.865
(1.00)
0.566
(1.00)
0.812
(1.00)
0.355
(1.00)
11p gain 5 (9%) 51 0.947
(1.00)
0.752
(1.00)
1
(1.00)
0.838
(1.00)
1
(1.00)
11q gain 7 (12%) 49 0.965
(1.00)
0.457
(1.00)
0.692
(1.00)
0.582
(1.00)
1
(1.00)
12p gain 22 (39%) 34 0.131
(1.00)
0.45
(1.00)
0.576
(1.00)
0.347
(1.00)
0.284
(1.00)
12q gain 11 (20%) 45 0.961
(1.00)
0.584
(1.00)
1
(1.00)
0.688
(1.00)
0.393
(1.00)
13q gain 15 (27%) 41 0.268
(1.00)
0.505
(1.00)
0.535
(1.00)
0.471
(1.00)
0.734
(1.00)
14q gain 7 (12%) 49 0.0259
(1.00)
0.766
(1.00)
1
(1.00)
0.775
(1.00)
0.558
(1.00)
15q gain 4 (7%) 52 0.557
(1.00)
0.166
(1.00)
1
(1.00)
0.669
(1.00)
0.202
(1.00)
16p gain 10 (18%) 46 0.553
(1.00)
0.416
(1.00)
0.736
(1.00)
0.405
(1.00)
0.394
(1.00)
16q gain 6 (11%) 50 0.373
(1.00)
0.853
(1.00)
0.675
(1.00)
0.468
(1.00)
0.488
(1.00)
17p gain 9 (16%) 47 0.302
(1.00)
0.396
(1.00)
0.275
(1.00)
0.269
(1.00)
0.28
(1.00)
17q gain 18 (32%) 38 0.99
(1.00)
0.765
(1.00)
1
(1.00)
0.816
(1.00)
0.75
(1.00)
18p gain 18 (32%) 38 0.384
(1.00)
0.854
(1.00)
0.394
(1.00)
0.152
(1.00)
0.653
(1.00)
18q gain 14 (25%) 42 0.502
(1.00)
0.985
(1.00)
0.766
(1.00)
0.329
(1.00)
0.856
(1.00)
19p gain 24 (43%) 32 0.612
(1.00)
0.637
(1.00)
0.0568
(1.00)
0.885
(1.00)
0.598
(1.00)
19q gain 28 (50%) 28 0.829
(1.00)
0.522
(1.00)
0.00582
(0.974)
0.207
(1.00)
0.592
(1.00)
20p gain 37 (66%) 19 0.637
(1.00)
0.209
(1.00)
0.155
(1.00)
0.721
(1.00)
0.439
(1.00)
20q gain 44 (79%) 12 0.802
(1.00)
0.697
(1.00)
0.755
(1.00)
0.384
(1.00)
0.195
(1.00)
21q gain 18 (32%) 38 0.605
(1.00)
0.752
(1.00)
0.394
(1.00)
0.933
(1.00)
0.257
(1.00)
22q gain 8 (14%) 48 0.0944
(1.00)
0.038
(1.00)
1
(1.00)
0.625
(1.00)
0.27
(1.00)
xp gain 18 (32%) 38 0.832
(1.00)
0.0074
(0.974)
0.25
(1.00)
0.412
(1.00)
1
(1.00)
xq gain 15 (27%) 41 0.943
(1.00)
0.604
(1.00)
0.766
(1.00)
0.735
(1.00)
0.204
(1.00)
1p loss 9 (16%) 47 0.515
(1.00)
0.584
(1.00)
1
(1.00)
0.187
(1.00)
0.614
(1.00)
1q loss 9 (16%) 47 0.32
(1.00)
0.454
(1.00)
0.431
(1.00)
0.129
(1.00)
0.613
(1.00)
3p loss 20 (36%) 36 0.34
(1.00)
0.784
(1.00)
0.249
(1.00)
0.939
(1.00)
0.517
(1.00)
3q loss 14 (25%) 42 0.451
(1.00)
0.302
(1.00)
0.113
(1.00)
0.622
(1.00)
0.731
(1.00)
4p loss 31 (55%) 25 0.286
(1.00)
0.98
(1.00)
1
(1.00)
0.151
(1.00)
1
(1.00)
4q loss 33 (59%) 23 0.185
(1.00)
0.537
(1.00)
0.583
(1.00)
0.141
(1.00)
0.518
(1.00)
5p loss 9 (16%) 47 0.0866
(1.00)
0.118
(1.00)
0.262
(1.00)
0.893
(1.00)
0.611
(1.00)
5q loss 18 (32%) 38 0.534
(1.00)
0.533
(1.00)
0.232
(1.00)
0.933
(1.00)
0.748
(1.00)
6p loss 5 (9%) 51 0.0969
(1.00)
0.508
(1.00)
0.653
(1.00)
0.838
(1.00)
0.403
(1.00)
6q loss 7 (12%) 49 0.254
(1.00)
0.682
(1.00)
1
(1.00)
0.774
(1.00)
0.182
(1.00)
7p loss 13 (23%) 43 0.897
(1.00)
0.734
(1.00)
0.545
(1.00)
0.719
(1.00)
1
(1.00)
7q loss 13 (23%) 43 0.0794
(1.00)
0.923
(1.00)
0.0647
(1.00)
0.921
(1.00)
0.716
(1.00)
8p loss 23 (41%) 33 0.884
(1.00)
0.653
(1.00)
0.783
(1.00)
0.292
(1.00)
0.223
(1.00)
8q loss 9 (16%) 47 0.915
(1.00)
0.577
(1.00)
0.065
(1.00)
0.268
(1.00)
0.284
(1.00)
9p loss 34 (61%) 22 0.447
(1.00)
0.118
(1.00)
0.579
(1.00)
0.0528
(1.00)
0.524
(1.00)
9q loss 39 (70%) 17 0.215
(1.00)
0.178
(1.00)
0.772
(1.00)
0.374
(1.00)
0.412
(1.00)
10p loss 23 (41%) 33 0.624
(1.00)
0.214
(1.00)
0.579
(1.00)
0.831
(1.00)
0.264
(1.00)
10q loss 21 (38%) 35 0.277
(1.00)
0.748
(1.00)
0.573
(1.00)
0.117
(1.00)
0.284
(1.00)
11p loss 26 (46%) 30 0.0291
(1.00)
0.23
(1.00)
0.106
(1.00)
0.194
(1.00)
0.883
(1.00)
11q loss 24 (43%) 32 0.425
(1.00)
0.573
(1.00)
0.585
(1.00)
0.31
(1.00)
0.768
(1.00)
12p loss 13 (23%) 43 0.405
(1.00)
0.907
(1.00)
0.751
(1.00)
1
(1.00)
0.302
(1.00)
12q loss 14 (25%) 42 0.356
(1.00)
0.191
(1.00)
1
(1.00)
0.677
(1.00)
0.731
(1.00)
13q loss 26 (46%) 30 0.834
(1.00)
0.831
(1.00)
0.79
(1.00)
0.889
(1.00)
1
(1.00)
14q loss 27 (48%) 29 0.288
(1.00)
0.533
(1.00)
0.0611
(1.00)
0.0875
(1.00)
0.0824
(1.00)
15q loss 32 (57%) 24 0.565
(1.00)
0.797
(1.00)
0.167
(1.00)
0.695
(1.00)
0.357
(1.00)
16p loss 32 (57%) 24 0.467
(1.00)
0.868
(1.00)
0.583
(1.00)
0.544
(1.00)
0.517
(1.00)
16q loss 37 (66%) 19 0.117
(1.00)
0.965
(1.00)
0.154
(1.00)
0.302
(1.00)
0.327
(1.00)
17p loss 34 (61%) 22 0.214
(1.00)
0.574
(1.00)
0.407
(1.00)
0.371
(1.00)
0.242
(1.00)
17q loss 17 (30%) 39 0.0756
(1.00)
0.562
(1.00)
1
(1.00)
0.464
(1.00)
0.74
(1.00)
18p loss 18 (32%) 38 0.778
(1.00)
0.699
(1.00)
0.142
(1.00)
0.0995
(1.00)
0.258
(1.00)
18q loss 20 (36%) 36 0.276
(1.00)
0.515
(1.00)
0.777
(1.00)
0.171
(1.00)
0.0362
(1.00)
19p loss 15 (27%) 41 0.926
(1.00)
0.505
(1.00)
0.371
(1.00)
1
(1.00)
0.737
(1.00)
19q loss 13 (23%) 43 0.68
(1.00)
0.437
(1.00)
0.188
(1.00)
0.515
(1.00)
0.397
(1.00)
20p loss 7 (12%) 49 0.885
(1.00)
0.823
(1.00)
0.431
(1.00)
0.884
(1.00)
0.557
(1.00)
20q loss 4 (7%) 52 0.824
(1.00)
0.886
(1.00)
0.615
(1.00)
1
(1.00)
1
(1.00)
21q loss 17 (30%) 39 0.145
(1.00)
0.175
(1.00)
0.00219
(0.865)
0.708
(1.00)
1
(1.00)
22q loss 33 (59%) 23 0.363
(1.00)
0.683
(1.00)
0.271
(1.00)
0.615
(1.00)
0.101
(1.00)
xp loss 15 (27%) 41 0.619
(1.00)
0.247
(1.00)
1
(1.00)
0.86
(1.00)
0.862
(1.00)
xq loss 18 (32%) 38 0.914
(1.00)
0.635
(1.00)
0.394
(1.00)
0.816
(1.00)
0.655
(1.00)
Methods & Data
Input
  • Copy number data file = broad_values_by_arm.txt from GISTIC pipeline

  • Processed Copy number data file = /xchip/cga/gdac-prod/tcga-gdac/jobResults/GDAC_Correlate_Genomic_Events_Preprocess/UCS-TP/19782451/transformed.cor.cli.txt

  • Clinical data file = /xchip/cga/gdac-prod/tcga-gdac/jobResults/Append_Data/UCS-TP/19775681/UCS-TP.merged_data.txt

  • Number of patients = 56

  • Number of significantly arm-level cnvs = 79

  • Number of selected clinical features = 5

  • Exclude regions that fewer than K tumors have mutations, K = 3

Survival analysis

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

Fisher's exact test

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

Q value calculation

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.

Download Results

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.

References
[1] Bland and Altman, Statistics notes: The logrank test, BMJ 328(7447):1073 (2004)
[2] Fisher, R.A., On the interpretation of chi-square from contingency tables, and the calculation of P, Journal of the Royal Statistical Society 85(1):87-94 (1922)
[3] Benjamini and Hochberg, Controlling the false discovery rate: a practical and powerful approach to multiple testing, Journal of the Royal Statistical Society Series B 59:289-300 (1995)