Mutation Analysis (MutSigCV v0.6)
Colon Adenocarcinoma (Primary solid tumor)
22 February 2013  |  analyses__2013_02_22
Maintainer Information
Citation Information
Maintained by Dan DiCara (Broad Institute)
Cite as Broad Institute TCGA Genome Data Analysis Center (2013): Mutation Analysis (MutSigCV v0.6). Broad Institute of MIT and Harvard. doi:10.7908/C1BZ647W
Overview
Introduction

This report serves to describe the mutational landscape and properties of a given individual set, as well as rank genes and genesets according to mutational significance. MutSigCV v0.6 was used to generate the results found in this report.

  • Working with individual set: COAD-TP

  • Number of patients in set: 155

Input

The input for this pipeline is a set of individuals with the following files associated for each:

  1. An annotated .maf file describing the mutations called for the respective individual, and their properties.

  2. A .wig file that contains information about the coverage of the sample.

Summary
Results
Target Coverage for Each Individual

The x axis represents the samples. The y axis represents the exons, one row per exon, and they are sorted by average coverage across samples. For exons with exactly the same average coverage, they are sorted next by the %GC of the exon. (The secondary sort is especially useful for the zero-coverage exons at the bottom).

Figure 1. 

Distribution of Mutation Counts, Coverage, and Mutation Rates Across Samples

Figure 2.  Patients counts and rates file used to generate this plot: COAD-TP.patients.counts_and_rates.txt

CoMut Plot

Figure 3.  Get High-res Image The matrix in the center of the figure represents individual mutations in patient samples, color-coded by type of mutation, for the significantly mutated genes. The rate of synonymous and non-synonymous mutations is displayed at the top of the matrix. The barplot on the left of the matrix shows the number of mutations in each gene. The percentages represent the fraction of tumors with at least one mutation in the specified gene. The barplot to the right of the matrix displays the q-values for the most significantly mutated genes. The purple boxplots below the matrix (only displayed if required columns are present in the provided MAF) represent the distributions of allelic fractions observed in each sample. The plot at the bottom represents the base substitution distribution of individual samples, using the same categories that were used to calculate significance.

Significantly Mutated Genes

Column Descriptions:

  • nnon = number of (nonsilent) mutations in this gene across the individual set

  • npat = number of patients (individuals) with at least one nonsilent mutation

  • nsite = number of unique sites having a non-silent mutation

  • nflank = number of noncoding mutations from this gene's flanking region, across the individual set

  • nsil = number of silent mutations in this gene across the individual set

  • p = p-value (overall)

  • q = q-value, False Discovery Rate (Benjamini-Hochberg procedure)

Table 1.  Get Full Table A Ranked List of Significantly Mutated Genes. Number of significant genes found: 8. Number of genes displayed: 35. Click on a gene name to display its stick figure depicting the distribution of mutations and mutation types across the chosen gene (this feature may not be available for all significant genes).

gene Nnon Nsil Nflank nnon npat nsite nsil nflank nnei fMLE p score time q
APC 1034160 289850 291 163 107 108 4 0 4 1 1.7e-15 500 0.11 3e-11
TP53 146475 42780 206 77 75 50 1 0 4 1 3.4e-15 260 0.11 3.1e-11
KRAS 93775 22940 106 59 59 9 0 0 1 0 5.6e-15 160 0.11 3.4e-11
NRAS 71455 18910 82 15 15 7 0 0 20 0.59 1.5e-11 48 0.1 6.8e-08
FBXW7 301010 82615 239 33 29 21 2 0 20 1.5 1.4e-10 81 0.11 5.2e-07
SMAD4 205220 57195 222 22 18 18 0 0 20 1.2 2.3e-07 54 0.11 0.0007
FAM123B 383780 113770 98 19 19 17 1 0 20 2 4.2e-07 77 0.11 0.0011
SOX9 146010 42935 52 9 9 9 0 0 20 1.5 1e-05 47 0.11 0.024
PIK3CA 402845 103075 396 33 26 18 1 0 20 1.2 0.000052 53 0.11 0.11
ACVR2A 192820 50840 222 10 9 8 1 0 10 1.1 0.0001 37 0.1 0.19
C15orf60 73625 21390 70 4 4 4 0 0 8 0 0.00018 17 0.092 0.3
ACVR1B 176080 50685 162 13 13 13 0 0 7 0.31 0.00024 33 0.1 0.36
BRAF 267685 76415 339 21 20 3 0 0 19 1.2 0.00029 42 0.1 0.4
TXNDC3 224905 54715 294 13 10 13 2 0 20 1.1 0.00048 36 0.11 0.62
SMAD2 173135 48360 199 11 10 8 1 0 20 1.4 0.00062 35 0.1 0.76
NAALAD2 279155 75330 381 13 11 13 0 0 4 0 0.00071 32 0.1 0.81
CASP8 225680 52855 213 11 10 10 0 0 20 0.88 0.00083 32 0.1 0.88
SAT1 64790 15500 117 4 4 2 0 0 20 0.91 0.00087 20 0.097 0.88
ESR1 188635 53475 146 12 11 12 0 0 20 0.67 0.0011 29 0.1 1
TMEM105 41230 14570 35 3 3 3 0 0 19 0.4 0.0014 14 0.09 1
TPX2 281015 75020 322 7 5 7 0 0 6 0 0.0014 23 0.1 1
MGC26647 92070 25420 24 7 7 7 0 0 13 1.4 0.0015 25 0.1 1
CRTC1 201190 63550 274 5 5 3 1 0 20 0.45 0.0015 26 0.1 1
ZNF593 25885 7750 43 3 3 3 0 0 20 0.75 0.0016 14 0.094 1
RNF17 608530 157015 672 11 10 11 1 0 5 0.2 0.0017 31 0.1 1
CCDC160 60295 14570 5 7 5 7 0 0 20 1.2 0.0022 18 0.095 1
GPHN 284580 85095 460 11 11 10 1 0 3 0.25 0.0023 30 0.1 1
C6orf211 164610 41850 100 5 5 5 0 0 19 0.41 0.0023 20 0.096 1
RP9 67115 15190 102 3 3 2 0 0 20 0.38 0.003 15 0.091 1
WBSCR17 214055 61690 211 17 17 16 2 0 18 1.5 0.003 39 0.11 1
ELF3 137795 37820 158 5 5 5 0 0 20 0.82 0.0034 24 0.099 1
TBC1D10C 86800 28830 92 4 4 2 0 0 20 1.1 0.0034 23 0.099 1
ATM 1157385 297445 1233 28 21 28 1 0 2 0.4 0.0034 56 0.11 1
MIER3 206460 53165 236 10 10 10 0 0 20 1.1 0.004 29 0.1 1
TCF7L2 230175 64945 273 13 11 13 3 0 20 1.4 0.0041 37 0.1 1
Methods & Data
Methods

In brief, we tabulate the number of mutations and the number of covered bases for each gene. The counts are broken down by mutation context category: four context categories that are discovered by MutSig, and one for indel and 'null' mutations, which include indels, nonsense mutations, splice-site mutations, and non-stop (read-through) mutations. For each gene, we calculate the probability of seeing the observed constellation of mutations, i.e. the product P1 x P2 x ... x Pm, or a more extreme one, given the background mutation rates calculated across the dataset. [1]

Download Results

This is an experimental feature. The full results of the analysis summarized in this report can be downloaded from the TCGA Data Coordination Center.

References
[1] TCGA, Integrated genomic analyses of ovarian carcinoma, Nature 474:609 - 615 (2011)