SNP Discovery

BCM-HGSC Mutation Discovery Pipeline

Single Nucleotide Mutations

Sequencing reads are compared with their respective amplicon reference sequences using a modification of SNPDetector, which employs a relaxed Het peak ratio threshold to compensate for possible heterogeneity of the tumor tissue sample. We use Polyphred 6.0b as a backup discovery method and capture any high-scoring variation missed (very rarely) by SNPDetector. Special mention should be made of our collaboration with Dr. Jinghui Zhang in the laboratory of Dr Ken Beautow, on the development and calibration of SNPdetector. This software has consistently outperformed other routines for the direct discovery of heterozygotes.

Putative polymorphisms accumulated from this analysis are annotated with the following information:

  1. Chromosome and global coordinates

  2. Coincidence with known variation (dbSNP current build, as well as local databases of newly identified variants)

  3. Functional information:

    1. gene compartment (intron, exon, splice junction)

    2. non-synonymous amino acid change if any

    3. position of non-synonymous amino acid in protein

    4. BLOSSUM62 score of variant amino acid compared to reference

Novel SNPs with recognizable functional potential (e.g., non-synonymous SNP or splice junctional variants) are further evaluated. First, they are visually inspected at the trace level and those that are not clearly noise are passed on to experimental validation, currently pyro-sequencing. We plan to resequence with Sanger reads the matched normal tissue in patients with mutations passing pyro-sequencing validation.

All putative genotypes of each individual at each mutation position, along with annotation and validation status will be stored in local databases. Reports are formatted for submission to common data repositories according to protocols jointly established.

Structural Variation Discovery

We are in the process of evaluating Polyphred 6.0b and a new module for SNPDetector designed for detecting intra-exonic indels in biallelic resequencing traces. One or both of these will be used for indel discovery and characterization. Genotype frequencies of constitutional variants (i.e., known SNPs) will be tracked since they might reveal commonly deleted genes or gene segments (LOH) through departures fro, Hardy-Weinberg equilibrium.

Quality Control

Sequencing coverage is a critical factor leading to variation discovery. We track coverage using the SNPDetector program rather than by a single base quality measure. Bases are judged to be covered in a given read if SNPDetector is able to make a call at any given position regardless of their Phred quality score (although there is a high correlation between Phred quality score and SNPDetector coverage).

Related Publications

Liu Y, Qin X, Song X-ZHenry, Jiang H, Shen Y, K Durbin J, et al. Bos taurus genome assembly. BMC Genomics. 2009 ;10:180.

Gibbs RA, Taylor JF, Van Tassell CP, Barendse W, Eversole KA, Gill CA, et al. Genome-wide survey of SNP variation uncovers the genetic structure of cattle breeds. Science. 2009 ;324(5926):528-32.

You M, Wang D, Liu P, Vikis H, James M, Lu Y, et al. Fine mapping of chromosome 6q23-25 region in familial lung cancer families reveals RGS17 as a likely candidate gene. Clin Cancer Res. 2009 ;15(8):2666-74.

Maheshwari M, Shi J, Badner JA, Skol A, Willour VL, Muzny DM, et al. Common and rare variants of DAOA in bipolar disorder. Am J Med Genet B Neuropsychiatr Genet. 2009 ;150B(7):960-6.

Meyer TE, Shiffman D, Morrison AC, Rowland CM, Louie JZ, Bare LA, et al. GOSR2 Lys67Arg is associated with hypertension in whites. Am J Hypertens. 2009 ;22(2):163-8.

Li Y, Wang H, Peng J, Gibbs RA, Lewis RAlan, Lupski JR, et al. Mutation survey of known LCA genes and loci in the Saudi Arabian population. Invest Ophthalmol Vis Sci. 2009 ;50(3):1336-43.

Chung CC, Shimmin L, Natarajan S, Hanis CL, Boerwinkle E, Hixson JE. Glucocorticoid receptor gene variant in the 3' untranslated region is associated with multiple measures of blood pressure. J Clin Endocrinol Metab. 2009 ;94(1):268-76.