We thank Dr. Agnes Viale and the MSKCC Genomics Core Laboratory for technical assistance. This protocol was developed with support from the Geoffrey Beene Cancer Research Center and the Farmer Family Foundation.

1. DNA and Reagent Preparation
Note: This protocol describes the simultaneous processing and analysis of 24 samples (e.g. 12 tumor/normal pairs) but can be adapted for smaller and larger batches. DNA samples may derive from FFPE or fresh frozen tissue, cytological specimens, or blood. Typically, both tumor and normal tissue from the same patient will be profiled together in order to distinguish somatic mutations from inherited polymorphisms. The protocol begins immediately following DNA extraction.
<ol>
	<li>Aliquot 50-250 ng (250 ng recommended) of extracted DNA per sample, diluted in 1x Tris-EDTA (pH 8.0) buffer to a final volume of 50 &#181;l, into separate Covaris round bottom tubes.</li>
	<li>Shear the DNA on the Covaris E220 instrument at the following shear settings: 360 sec at 10 duty factor, 175 PIP, and 200 cycles per burst.</li>
	<li>Thaw AMPure XP Beads (Table 2) to room temperature. Thaw all buffers and enzymes on ice prior to the appropriate step.</li>
</ol>
2. Library Preparation - End Repair Module
<ol>
	<li>Prepare the end repair master mix in a sterile microcentrifuge tube using the following volumes per sample: 10 &#181;l of 10x NEBNext End Repair Reaction Buffer (Table 2), 5 &#181;l NEBNext End Repair Enzyme Mix (Table 2), 35 &#181;l sterile nuclease free H2O.</li>
	<li>Aliquot the 50 &#181;l of each sheared DNA into separate wells of a 96-well plate. Add 50 &#181;l of the prepared end repair master mix to each reaction. Incubate in a thermocycler for 30 min at 20 &#176;C.</li>
</ol>
3. Post End Repair Clean-up
<ol>
	<li>Add 2x volume (200 &#181;l) of AMPure XP Beads (Table 2) to each sample and gently pipette the entire volume up and down 10x using a multichannel pipettor. Incubate at RT for 15 min.</li>
	<li>Place on the magnetic stand (Table 2) at RT for 15 min until the sample appears clear.</li>
	<li>Gently remove and discard clear supernatant taking care not to disturb beads. Some liquid may remain in wells.</li>
	<li>With tubes on the stand, gently add 200 &#181;l of freshly prepared 80% ethanol to each sample and incubate plate at room temperature for 30 sec. Carefully, remove ethanol by pipette.</li>
	<li>Repeat step 5, for a total of 2 ethanol washes. Ensure all ethanol has been removed. Remove from the magnetic stand and let dry at RT for 5 min.</li>
	<li>Resuspend dried beads with 44.5 &#181;l sterile H2O. Gently, pipette entire volume up and down 10x mixing thoroughly. Ensure beads are no longer attached to the side of the well.</li>
	<li>Incubate resuspended beads at RT for 2 min. Place on magnetic stand for 5 min until the sample appears clear.</li>
	<li>Gently transfer 42 &#181;l of clear supernatant containing the sample to a new well. 
	The procedure may be safely stopped at this step and samples stored at -20 &#176;C for up to 7 days. To restart, thaw frozen samples on ice before proceeding.</li>
</ol>
4. Library Preparation - dA-Tailing Module
<ol>
	<li>Prepare the dA-Tailing master mix in a sterile microcentrifuge tube using the following volumes per sample: 5 &#181;l of 10x NEBNext dA-Tailing Reaction Buffer (Table 2), 3 &#181;l (3&#39;-5&#39; exo-) Klenow Fragment (Table 2).</li>
	<li>Add 8 &#181;l of the prepared dA-Tailing master mix to each well containing 42 &#181;l of end-repaired DNA. Incubate in a thermocycler for 30 min at 37 &#176;C.</li>
	<li>Perform a post dA-Tailing clean-up: Repeat steps 3.1-3.8 using 2x volume (100 &#181;l) of beads, but resuspend in sterile H2O to a final volume of 33.75 &#181;l. 
	The procedure may be safely stopped at this step and samples stored at -20 &#176;C for up to 7 days. To restart, thaw frozen samples on ice before proceeding.</li>
</ol>
5. Library Preparation - Adapter Ligation Module
<ol>
	<li>Prepare the Ligation master mix in a sterile microcentrifuge tube using the following volumes per sample: 10 &#181;l of 5x NEBNext Quick Ligation Reaction Buffer (New England Biolabs, Table 2), 5 &#181;l Quick T4 DNA Ligase (Table 2).</li>
	<li>Add 15 &#181;l of the Ligation master mix to each well containing the dA-Tailed DNA.</li>
	<li>Add 1.25 &#181;l of the appropriate 25 &#181;M NEXTflex barcoded adapter (Table 2) to each well. Each sample should receive a separate barcoded adapter. Incubate in a thermocycler for 15 min at 20 &#176;C.</li>
	<li>Perform a post adapter ligation clean-up: Repeat steps 3.1-3.8 using 1x volume (50 &#181;l) of beads, but resuspend in sterile H2O to a final volume of 50 &#181;l.</li>
	<li>Repeat steps 3.1-3.8 for a second 1x volume (50 &#181;l) clean-up, but resuspend in sterile H2O to a final volume of 23 &#181;l. 
	The procedure may be safely stopped at this step and samples stored at -20 &#176;C for up to 7 days. To restart, thaw frozen samples on ice before proceeding.</li>
</ol>
6. Library Amplification
<ol>
	<li>Prepare the KAPA HiFi master mix in a sterile microcentrifuge tube using the following volumes per sample: 25 &#181;l 2x KAPA HiFi HS RM (Table 2), 1 &#181;l each of 100 &#181;M Primer 1 and Primer 2 (Table 3).</li>
	<li>Add 27 &#181;l of the master mix to each well containing the ligation product. Set pipette to 50 &#181;l and gently pipette the entire volume up and down 10x.</li>
	<li>Place in the thermocycler for the following PCR cycles: initial denaturation of 45 sec at 98 &#176;C, 10 cycles of 15 sec at 98 &#176;C, 30 sec at 60 &#176;C, 30 sec at 72 &#176;C, and final extension of 1 min at 72 &#176;C.</li>
	<li>Perform a post amplification clean-up: Repeat steps 3.1-3.8 using 1x volume (50 &#181;l) of AMPure XP Beads, and resuspend in sterile water to a final volume of 30 &#181;l.</li>
	<li>Quantify the concentration of each library using Qubit Broad Range Assay (Life Technologies, Table 2) according to the manufacturer&#39;s instructions.</li>
</ol>
7. Roche Nimblegen SeqCAP EZ Library Hybridization, Wash, and Amplification
<ol>
	<li>Thaw universal and index oligonucleotide blockers (Table 3), SeqCap EZ Library, and Nimblegen capture components (COT DNA, 2x hybridization buffer, and hybridization component A, Table 2) on ice. All oligonucleotide blockers are stored in working stocks of 1 mM. Our SeqCap EZ Library is a custom-designed pool of capture probes encompassing all protein-coding exons of 279 cancer-associated genes, though other custom and catalog designs may be used.</li>
	<li>Prepare a 1 mM pool of index oligonucleotide blockers corresponding to the specific barcoded adaptor sequences used in library preparation. Combine 1 &#181;l of each blocker, and vortex.</li>
	<li>In a new sterile 1.5 ml microcentrifuge tube, add 5 &#181;l of 1mg/ml COT DNA, 2 &#181;l of 1 mM universal blocker and 2 &#181;l of the 1 mM pool of index oligonucleotide blockers.</li>
	<li>Pool 24 barcoded libraries into a single reaction. Add a total of 1-3 &#181;g of pooled, barcoded sequence library (prepared above) to the capture mixture. For 24 samples, 100 ng per sample is recommended. When pooling tumors and matched normal samples, we recommended adding more tumor library (e.g. 100 ng per tumor and 50 ng per normal) so that they are sequenced to greater depth of coverage.</li>
	<li>Close the tube&#39;s cap and make 15-20 holes in the top of the cap with an 18-20 G or smaller needle. Speed Vac the amplified sample library/COT DNA/oligos in a DNA vacuum concentrator on high heat (60 &#176;C).</li>
	<li>Once completely dried down, rehydrate with 7.5 &#181;l of 2x hybridization buffer and 3 &#181;l of hybridization component A. The tube should now contain the following: 5 &#181;g COT DNA, variable &#181;g amplified sample library, 1000 pmol of Universal and Index Oligos, 7.5 &#181;l 2x hybridization buffer, 3 &#181;l hybridization component A in a total volume of 10.5 &#181;l.</li>
	<li>Cover the holes in the cap with a small piece of laboratory tape. Vortex the sample for 10 sec and centrifuge at maximum speed for 10 sec.</li>
	<li>Incubate the sample at 95 &#176;C for 10 min to denature the DNA, followed by centrifugation at maximum speed for 10 sec at RT.</li>
	<li>In a new 0.2 ml PCR strip tube, prepare 2.25 &#181;l aliquot of SeqCap EZ Library capture probes and 2.25 &#181;l of sterile nuclease free H2O and add the contents of step 7.7. Gently pipette the entire volume up and down 10x.</li>
	<li>Incubate in a thermal cycler at 47 &#176;C for 48-96 hr. Set and maintain the thermal cycler&#39;s lid temperature to 57 &#176;C.</li>
	<li>Follow steps in Chapter 6 of the NimbleGen SeqCap EX Library SR User&#39;s Guide v3.0 (Table 2) for washing and recovery of the captured DNA.</li>
	<li>Follows steps in Chapter 7 of the Roche NimbleGen protocol for amplification of the captured DNA with the following modifications:
	<ul type="disc">
		<li>Use 2x KAPA HiFi Polymerase (Table 2) instead of Phusion High-Fidelity PCR Master Mix</li>
		<li>Use 100 &#181;M Primer 1 and Primer 2 listed in Table 3.</li>
		<li>Run the following PCR conditions: initial denaturation at 98 &#176;C for 30 sec, 11 cycles of 98 &#176;C for 10 sec, 60 &#176;C for 30 sec, 72 &#176;C for 30 sec, and final extension at 72 &#176;C for 5 min. Hold at 4 &#176;C.</li>
	</ul>
	</li>
	<li>Quantify the concentration of the amplified captured DNA using Qubit High Sensitivity Assay (Table 2) according to the manufacturer&#39;s instructions.</li>
	<li>Analyze the captured DNA quality using Agilent Bioanalyzer DNA HS chip (Table 2) according to the manufacturer&#39;s instructions.
	<ul type="disc">
		<li>The PCR yield should be at least 100 ng (range 100-1,000 ng).</li>
		<li>The average fragment length should be between 150-400 base pairs.</li>
	</ul>
	</li>
</ol>
8. Illumina Hi-Seq Sequencing
<ol>
	<li>Dilute the amplified captured DNA to 2 pM, and load the sample on a single lane of an Illumina HiSeq 2000 flow cell. Sequence paired-end 75 base pair reads according to the manufacturer&#39;s instructions. In our experience, one lane is sufficient to produce 500-700x unique sequence coverage per sample across 279 genes for a pool of 24 samples.</li>
	<li>Illumina software (Real Time Analysis) will convert high-resolution images to cluster intensities to base calls and quality scores. Use CASAVA to demultiplex data according to barcode identity and generate individual FASTQ files for each sample.</li>
</ol>
9. Data Analysis
<ol>
	<li>Trim adapter sequences from the 3&#39; end of sequence reads in FASTQ files using Cutadapt (<a href="http://code.google.com/p/cutadapt/">http://code.google.com/p/cutadapt/</a>). The minimum overlap length (-O) is 3 with an error rate (-e) of 1%. The minimum base length for retention of paired reads after trimming is 25.</li>
	<li>Align sequence reads to the reference human genome hg19 using the Burrows-Wheeler Aligner tool (BWA)12. Perform post-processing steps including duplicate marking, local multiple sequence realignment, and base quality score recalibration using the Genome Analysis Toolkit (GATK) according to their standard best practices (<a href="http://www.broadinstitute.org/gatk/guide/topic?name=best-practices">http://www.broadinstitute.org/gatk/guide/topic?name=best-practices</a>)13. When profiling multiple samples from the same patient (e.g.&#160;matched tumor and normal), local multiple sequence realignment should be performed jointly. The output of these post-processing steps is a separate BAM file for each sample, which contains all sequence, quality, and alignment information and can be loaded directly into the Integrative Genomics Viewer (IGV) for visualization of reads and sequence variants14.</li>
	<li>Calculate sequence performance metrics using Picard tools (<a href="http://picard.sourceforge.net/">http://picard.sourceforge.net/</a>). Informative metrics include alignment rate, fragment size distribution, GC-content associated coverage bias, on-target capture specificity, PCR duplicate rate, library complexity, and mean target coverage. Representative values are displayed in Figure 1. Allele frequencies at common polymorphism sites calculated using DepthOfCoverage from GATK are used to monitor contamination from unrelated DNA and to ensure that matched samples came from the same patients. 
	The following computational analyses depend on the presence of matched tumor and normal samples for the detection of somatic genetic alterations. Unmatched tumors may also be analyzed using a separate normal control sample, but most somatic mutations cannot easily be distinguished from inherited sequence variants.</li>
	<li>Call somatic single nucleotide variants (SNVs) for each tumor-normal pair using a somatic mutation caller, such as muTect15, SomaticSniper16, or Strelka17. We use muTect with modified filtering criteria, allowing for variants with low (non-zero) allele counts in the normal sample as long as the allele frequency is at least 5x&#160;greater in the tumor. Variants recurrently called in a panel of normal DNAs are removed as likely systematic artifacts of sequencing or alignment. Each SNV passing filters is then carefully reviewed using IGV.</li>
	<li>Call somatic indels for each tumor-normal pair using algorithms such as SomaticIndelDetector13 , Dindel18, or SOAPindel19. We use SomaticIndelDetector with the following criteria: at least 10% allele frequency with at least 3 supporting reads, or at least 2% allele frequency with at least 10 supporting reads. We allow for indels with low (non-zero) allele counts in the normal sample as long as the allele frequency is at least 5x&#160;greater in the tumor. Indels called in a panel of normal DNAs are removed, as above. Each indel passing filters is then carefully reviewed using IGV.</li>
	<li>Extrapolate copy number status from the sequence coverage at target exons. For each individual sample, perform a Loess normalization of mean sequence coverage for all target exons, regression over GC percentage. Adjust the unnormalized coverage values by subtracting the Loess fit and adding the sample-wide median coverage. Calculate the ratio in normalized sequence coverage across all target exons for tumor-normal pairs. In the absence of a matched normal sample and/or to decrease noise, other gender-matched diploid normal controls lacking germline copy number variants may be used. Determine somatic copy number gains and losses based on increases or decreases in the coverage ratio.</li>
	<li>Call somatic rearrangements involving cancer genes where at least one genomic breakpoint falls within or near a targeted interval of the genome. Suggested algorithms include GASV20, Breakdancer21, CREST22 , and dRanger23. Intrachromosomal inversions, deletions, and tandem duplications may be distinguished based on the relative position and orientation of supporting reads in discordant pairs. All candidate rearrangements should be manually reviewed using IGV. We recommend experimental validation by PCR and Sanger sequencing across the putative rearrangement breakpoints.</li>
</ol>

<ol>
	<li>Stratton, M. R., Campbell, P. J., Futreal, P. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+cancer+genome.">The cancer genome.</a> Nature. 458 (7239), 719-724 (2009).</li><li>Meyerson, M., Gabriel, S., Getz, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Advances+in+understanding+cancer+genomes+through+second-generation+sequencing.">Advances in understanding cancer genomes through second-generation sequencing.</a> Nat. Rev. Genet. 11 (10), 685-696 (2010).</li><li>Thomas, R. K., Baker, A. C., Debiasi, R. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=High-throughput+oncogene+mutation+profiling+in+human+cancer.">High-throughput oncogene mutation profiling in human cancer.</a> Nat. Genet. 39 (3), 347-351 (2007).</li><li>Macconaill, L. E., Campbell, C. D., Kehoe, S. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Profiling+critical+cancer+gene+mutations+in+clinical+tumor+samples.">Profiling critical cancer gene mutations in clinical tumor samples.</a> PLoS One. 4 (11), e7887 (2009).</li><li>Dias-Santagata, D., Akhavanfard, S., David, S. S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Rapid+targeted+mutational+analysis+of+human+tumours:+a+clinical+platform+to+guide+personalized+cancer+medicine.">Rapid targeted mutational analysis of human tumours: a clinical platform to guide personalized cancer medicine.</a> EMBO Mol. Med. 2 (5), 146-158 (2010).</li><li>Macconaill, L. E., Van Hummelen, P., Meyerson, M., Hahn, W. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Clinical+implementation+of+comprehensive+strategies+to+characterize+cancer+genomes:+opportunities+and+challenges.">Clinical implementation of comprehensive strategies to characterize cancer genomes: opportunities and challenges.</a> Cancer Discov. 1 (4), 297-311 (2011).</li><li>Taylor, B. S., Ladanyi, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Clinical+cancer+genomics:+how+soon+is+now.">Clinical cancer genomics: how soon is now.</a> J. Pathol. 223 (2), 318-326 (2011).</li><li>Gnirke, A., Melnikov, A., Maguire, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Solution+hybrid+selection+with+ultra-long+oligonucleotides+for+massively+parallel+targeted+sequencing.">Solution hybrid selection with ultra-long oligonucleotides for massively parallel targeted sequencing.</a> Nat. Biotechnol. 27 (2), 182-189 (2009).</li><li>Mamanova, L., Coffey, A. J., Scott, C. E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Target-enrichment+strategies+for+next-generation+sequencing.">Target-enrichment strategies for next-generation sequencing.</a> Nat. Methods. 7 (2), 111-118 (2010).</li><li>Craig, D. W., Pearson, J. V., Szelinger, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Identification+of+genetic+variants+using+bar-coded+multiplexed+sequencing.">Identification of genetic variants using bar-coded multiplexed sequencing.</a> Nat. Methods. 5 (10), 887-893 (2008).</li><li>Wagle, N., Berger, M. F., Davis, M. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=High-throughput+detection+of+actionable+genomic+alterations+in+clinical+tumor+samples+by+targeted,+massively+parallel+sequencing.">High-throughput detection of actionable genomic alterations in clinical tumor samples by targeted, massively parallel sequencing.</a> Cancer Discov. 2 (1), 82-93 (2012).</li><li>Li, H., Durbin, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Fast+and+accurate+short+read+alignment+with+Burrows-Wheeler+transform.">Fast and accurate short read alignment with Burrows-Wheeler transform.</a> Bioinformatics. 25 (14), 1754-1760 (2009).</li><li>Depristo, M. A., Banks, E., Poplin, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+framework+for+variation+discovery+and+genotyping+using+next-generation+DNA+sequencing+data.">A framework for variation discovery and genotyping using next-generation DNA sequencing data.</a> Nat. Genet. 43 (5), 491-498 (2011).</li><li>Robinson, J. T., Thorvaldsdottir, H., Winckler, W., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Integrative+genomics+viewer.">Integrative genomics viewer.</a> Nat. Biotechnol. 29 (1), 24-26 (2011).</li><li>Cibulskis, K., Lawrence, M. S., Carter, S. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Sensitive+detection+of+somatic+point+mutations+in+impure+and+heterogeneous+cancer+samples.">Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples.</a> Nat. Biotechnol. 31 (3), 213-219 (2013).</li><li>Larson, D. E., Harris, C. C., Chen, K., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=SomaticSniper:+identification+of+somatic+point+mutations+in+whole+genome+sequencing+data.">SomaticSniper: identification of somatic point mutations in whole genome sequencing data.</a> Bioinformatics. 28 (3), 311-317 (2012).</li><li>Saunders, C. T., Wong, W. S., Swamy, S., Becq, J., Murray, L. J., Cheetham, R. K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Strelka:+Accurate+somatic+small-variant+calling+from+sequenced+tumor-normal+sample+pairs.">Strelka: Accurate somatic small-variant calling from sequenced tumor-normal sample pairs.</a> Bioinformatics. 28 (14), 1811-1817 (2012).</li><li>Albers, C. A., Lunter, G., Macarthur, D. G., Mcvean, G., Ouwehand, W. H., Din del Durbin, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Accurate+indel+calls+from+short-read+data.">Accurate indel calls from short-read data.</a> Genome Res. 21 (6), 961-973 (2011).</li><li>Li, S., Li, R., Li, H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Efficient+identification+of+indels+from+short+paired+reads.">Efficient identification of indels from short paired reads.</a> Genome Res. 23 (1), 195-200 (2013).</li><li>Sindi, S., Helman, E., Bashir, A., Raphael, B. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+geometric+approach+for+classification+and+comparison+of+structural+variants.">A geometric approach for classification and comparison of structural variants.</a> Bioinformatics. 25 (12), 222-230 (2009).</li><li>Chen, K., Wallis, J. W., Mclellan, M. D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=BreakDancer:+an+algorithm+for+high-resolution+mapping+of+genomic+structural+variation.">BreakDancer: an algorithm for high-resolution mapping of genomic structural variation.</a> Nat. Methods. 6 (9), 677-681 (2009).</li><li>Wang, J., Mullighan, C. G., Easton, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=CREST+maps+somatic+structural+variation+in+cancer+genomes+with+base-pair+resolution.">CREST maps somatic structural variation in cancer genomes with base-pair resolution.</a> Nat. Methods. 8 (8), 652-654 (2011).</li><li>Drier, Y., Lawrence, M. S., Carter, S. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Somatic+rearrangements+across+cancer+reveal+classes+of+samples+with+distinct+patterns+of+DNA+breakage+and+rearrangement-induced+hypermutability.">Somatic rearrangements across cancer reveal classes of samples with distinct patterns of DNA breakage and rearrangement-induced hypermutability.</a> Genome Res. 23 (2), 228-235 (2012).</li><li>Wang, Z., Gerstein, M., Snyder, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=RNA-Seq:+a+revolutionary+tool+for+transcriptomics.">RNA-Seq: a revolutionary tool for transcriptomics.</a> Nat. Rev. Genet. 10 (1), 57-63 (2009).</li><li>Roychowdhury, S., Iyer, M. K., Robinson, D. R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Personalized+oncology+through+integrative+high-throughput+sequencing:+a+pilot+study.">Personalized oncology through integrative high-throughput sequencing: a pilot study.</a> Sci. Transl. Med. 3 (111), 111ra121 (2011).</li><li>Kerick, M., Isau, M., Timmermann, B., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Targeted+high+throughput+sequencing+in+clinical+cancer+settings:+formaldehyde+fixed-paraffin+embedded+(FFPE)+tumor+tissues,+input+amount+and+tumor+heterogeneity.">Targeted high throughput sequencing in clinical cancer settings: formaldehyde fixed-paraffin embedded (FFPE) tumor tissues, input amount and tumor heterogeneity.</a> BMC Med. Genomics. 4, 68 (2011).</li></ol>

Dr. Berger has received consulting fees and research funding from Foundation Medicine, a cancer diagnostics company.

Our IMPACT assay produces a high alignment rate, high on-target rate, high target coverage, and high sensitivity for detecting mutations, indels, and copy number alterations. We have demonstrated the capability of our IMPACT assay to sequence DNA from both fresh frozen and archived FFPE samples of low DNA input. By performing targeted exon sequencing of key cancer-associated genes, one can achieve very deep sequence coverage for the exons of these most critical genes thereby maximizing the ability to detect low frequency...

The identification of &#34;driver&#34; tumor genetic events in key oncogenes and tumor suppressor genes plays an essential role in the diagnosis and treatment of many cancers1. Large-scale research efforts utilizing massively parallel &#34;next-generation&#34; sequencing have enabled the identification of many such cancer-associated genes in recent years2. However, these sequencing platforms typically require large quantities of DNA isolated from fresh frozen tissues, thus posing a major limitation in characterizing and analyzing DNA mutations from preserved tissues, such as formalin-fixed paraffin embedded (FFPE) tumor samples. Improved efforts to efficiently and reliably characterize &#34;actionable&#34; genomic information from FFPE tumor samples will enable the retrospective analysis of previously-banked specimens and further encourage individualized approaches to cancer management.
Traditionally, molecular diagnostic laboratories have relied on time-consuming, low-throughput methodologies such as Sanger sequencing and real-time PCR for DNA mutation profiling. More recently, higher-throughput methods utilizing multiplexed PCR or mass spectrometric genotyping have been developed to investigate recurrent somatic mutations in key cancer genes3-5. These approaches, however, are limited in that only predesignated &#34;hotspot&#34; mutations are assayed, making them unsuitable for detecting inactivating mutations in tumor suppressor genes. Massively parallel sequencing offers several advantages over these strategies including the ability to interrogate entire exons for both common and rare mutations, the ability to reveal additional classes of genomic alterations such as copy number gains and losses, and greater detection sensitivity in heterogeneous samples6, 7. Whole genome sequencing represents the most comprehensive approach for mutation discovery, though it is relatively expensive and incurs large computational demands for data analysis and storage.
For clinical applications, where only a small fraction of the genome may be of clinical interest, two particular innovations in sequencing technology have been transformative. First, through hybridization-based exon capture, one can isolate DNA corresponding to key cancer-associated genes for targeted mutation profiling8, . Second, through ligation of molecular barcodes (i.e. DNA sequences 6-8 nucleotides in length), one can pool hundreds of samples per sequencing run and fully take advantage of the ever-increasing capacity of massively parallel sequencing instruments10. When combined, these innovations enable tumors to be profiled for lower cost and at higher throughput, with smaller computational requirements11. Further, by redistributing sequence coverage to only those genes most critical to the particular application, one can achieve greater sequencing depth for higher detection sensitivity for low allele frequency events.
Here we describe our IMPACT assay (Integrated Mutation Profiling of Actionable Cancer Targets), which utilizes exon capture on barcoded sequence library pools by hybridization using custom oligonucleotides to capture all protein-coding exons and select introns of 279 key cancer-associated genes (Table 1). This strategy enables the identification of mutations, indels, copy number alterations, and select structural rearrangements involving these 279 genes. Our method is compatible with DNA isolated from both fresh frozen and FFPE tissue as well as fine needle aspirates and other cytology specimens.

<table border="0" cellpadding="0" cellspacing="0" width="812">
	<tbody> 
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">NEBNext End Repair Module</td>
			<td style="width:247px;">New England Biolabs</td>
			<td style="width:137px;">E6050L</td>
			<td style="width:79px;"> </td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">NEBNext dA-Tailing Module</td>
			<td style="width:247px;">New England Biolabs</td>
			<td style="width:137px;">E6053L</td>
			<td style="width:79px;"> </td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">NEBNext Quick Ligation Module</td>
			<td style="width:247px;">New England Biolabs</td>
			<td style="width:137px;">E6056L</td>
			<td style="width:79px;"></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Agencourt AMPure XP</td>
			<td style="width:247px;">Beckman Coulter Genomics</td>
			<td style="width:137px;"> </td>
			<td style="width:79px;"> </td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">NEXTflex PCR-Free Barcodes &#8211; 24</td>
			<td style="width:247px;">Bioo Scientific</td>
			<td style="width:137px;">514103</td>
			<td style="width:79px;"> </td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">HiFi Library Amplification Kit</td>
			<td style="width:247px;">KAPA Biosystems</td>
			<td style="width:137px;">KK2612</td>
			<td style="width:79px;"> </td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">COT Human DNA, Fluorometric Grade</td>
			<td style="width:247px;">Roche Diagnostics</td>
			<td style="width:137px;">05 480 647 001</td>
			<td style="width:79px;"> </td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">NimbleGen SeqCap EZ Hybridization and Wash kit</td>
			<td style="width:247px;">Roche NimbleGen</td>
			<td style="width:137px;">05 634 261 001</td>
			<td style="width:79px;"> </td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">SeqCap EZ Library Baits</td>
			<td style="width:247px;">Roche NimbleGen</td>
			<td style="width:137px;"> </td>
			<td style="width:79px;"> </td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">QIAquick PCR Purification Kit</td>
			<td style="width:247px;">Qiagen</td>
			<td style="width:137px;">28104</td>
			<td style="width:79px;"> </td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Qubit dsDNA Broad Range (BR) Assay Kit</td>
			<td style="width:247px;">Life Technologies</td>
			<td style="width:137px;">Q32850</td>
			<td style="width:79px;"> </td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Qubit dsDNA High Sensitivity (HS) Assay Kit</td>
			<td style="width:247px;">Life Technologies</td>
			<td style="width:137px;">Q32851</td>
			<td style="width:79px;"> </td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Agilent DNA HS Kit</td>
			<td style="width:247px;">Agilent Technologies</td>
			<td style="width:137px;">5067-4626, 4627</td>
			<td style="width:79px;"> </td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Agilent 2100 Bioanalyzer</td>
			<td style="width:247px;">Agilent Technologies</td>
			<td style="width:137px;"> </td>
			<td style="width:79px;"> </td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Covaris E220</td>
			<td style="width:247px;">Covaris</td>
			<td style="width:137px;"> </td>
			<td style="width:79px;"> </td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Magnetic Stand-96</td>
			<td style="width:247px;">Ambion</td>
			<td style="width:137px;">AM10027</td>
			<td style="width:79px;"> </td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Illumina Hi-Seq 2000</td>
			<td style="width:247px;">Illumina</td>
			<td style="width:137px;"> </td>
			<td style="width:79px;"> </td>
		</tr>
	</tbody>
</table>

detecting somatic genetic alterations in tumor specimens by exon capture and massively parallel sequencing

Efforts to detect and investigate key oncogenic mutations have proven valuable to facilitate the appropriate treatment for cancer patients. The establishment of high-throughput, massively parallel "next-generation" sequencing has aided the discovery of many such mutations. To enhance the clinical and translational utility of this technology, platforms must be high-throughput, cost-effective, and compatible with formalin-fixed paraffin embedded (FFPE) tissue samples that may yield small amounts of degraded or damaged DNA. Here, we describe the preparation of barcoded and multiplexed DNA libraries followed by hybridization-based capture of targeted exons for the detection of cancer-associated mutations in fresh frozen and FFPE tumors by massively parallel sequencing. This method enables the identification of sequence mutations, copy number alterations, and select structural rearrangements involving all targeted genes. Targeted exon sequencing offers the benefits of high throughput, low cost, and deep sequence coverage, thus conferring high sensitivity for detecting low frequency mutations.

One pool of 24 barcoded sequence libraries (12 tumor-normal pairs) was captured using probes corresponding to all protein-coding exons of 279 cancer genes and sequenced as 2 x 75 bp reads on a single lane of a HiSeq 2000 flow cell. Tumor and normal libraries were pooled in a 2:1 ratio. Sample performance metrics for a pool of frozen tumor DNA samples are shown in Figure 1, including alignment rate, fragment size distribution, on-target capture specificity, and mean target coverage. Example somatic mutati...

Watch this Scientific Journal Video about Detecting Somatic Genetic Alterations in Tumor Specimens by Exon Capture and Massively Parallel Sequencing at JoVE.com

Detecting Somatic Genetic Alterations in Tumor Specimens by Exon Capture and Massively Parallel Sequencing

We describe the preparation of barcoded DNA libraries and subsequent hybridization-based exon capture for detection of key cancer-associated mutations in clinical tumor specimens by massively parallel "next generation" sequencing. Targeted exon sequencing offers the benefits of high throughput, low cost, and deep sequence coverage, thus yielding high sensitivity for detecting low frequency mutations.

Efforts  to detect and investigate key oncogenic mutations have proven valuable to  facilitate the appropriate treatment for cancer patients. The ...

detecting-somatic-genetic-alterations-tumor-specimens-exon-capture

Research

JoVE Journal

Biology

19.3K Views.  Memorial Sloan-Kettering Cancer Center. We describe the preparation of barcoded DNA libraries and subsequent hybridization-based exon capture for detection of key cancer-associated mutations in clinical tumor specimens by massively parallel "next generation" sequencing. Targeted exon sequencing offers the benefits of high throughput, low cost, and deep sequence coverage, thus yielding high sensitivity for detecting low frequency mutations.