Wild-type Blocking PCR Combined with Direct Sequencing as a Highly Sensitive Method for Detection of Low-Frequency Somatic Mutations

Summary

Wild-type blocking PCR followed by direct sequencing offers a highly sensitive method of detection for low frequency somatic mutations in a variety of sample types.

Abstract

Accurate detection and identification of low frequency mutations can be problematic when assessing residual disease after therapy, screening for emerging resistance mutations during therapy, or when patients have few circulating tumor cells. Wild-type blocking PCR followed by sequencing analysis offers high sensitivity, flexibility, and simplicity as a methodology for detecting these low frequency mutations. By adding a custom designed locked nucleic acid oligonucleotide to a new or previously established conventional PCR based sequencing assay, sensitivities of approximately 1 mutant allele in a background of 1,000 WT alleles can be achieved (1:1,000). Sequencing artifacts associated with deamination events commonly found in formalin fixed paraffin embedded tissues can be partially remedied by the use of uracil DNA glycosylase during extraction steps. The optimized protocol here is specific for detecting MYD88 mutation, but can serve as a template to design any WTB-PCR assay. Advantages of the WTB-PCR assay over other commonly utilized assays for the detection of low frequency mutations including allele specific PCR and real-time quantitative PCR include fewer occurrences of false positives, greater flexibility and ease of implementation, and the ability to detect both known and unknown mutations.

Introduction

Sanger sequencing has traditionally been the gold standard in testing for both known and unknown somatic mutations. One of the limitations of Sanger sequencing is its limit of detection (~ 10 - 20% mutant allele in a background of WT)¹. This level of sensitivity is inappropriate for detecting low level somatic mutations that may be present in samples from premalignant tissues or patients with few circulating tumor cells, or when bone marrow (BM) is patchy. This also makes assessing residual disease after therapy or detecting emerging resistance mutations during therapy difficult by conventional sequencing alone². By replacing conventional PCR with locked nucleic acid (LNA)-mediated wild-type blocking PCR (WTB-PCR) in Sanger sequencing, sensitivities of up to 0.1% mutant allele in a background of WT can be achieved^2,3,4. In WTB-PCR, enrichment for mutant alleles is achieved via the addition of a short (~ 10 - 14 NT) blocking (LNA) oligonucleotide that binds preferentially to WT DNA thereby preventing amplification of WT DNA. The mutant enriched WTB-PCR product can then be sequenced. By blocking WT DNA rather than selecting for specific mutations WTB-PCR allows for enrichment of both known and unknown mutations present in minute cell fractions.

Multiple methods are currently used for detecting mutations in small cells fractions. This includes allele-specific PCR, amplification-refractory mutation system (ARMS), denaturing high performance liquid chromatography (DHPLC), beads, emulsions, amplification, and magnetics (BEAMing), electric field-induced release and measurement (EFIRM), high resolution melting point and others. However, most of these methods are limited by false-positives and the ability to only detect one mutation that the assay was designed for⁴. WTB-PCR, however, allows the user to visualize sequencing traces which enables the detection of multiple mutation types and can aid in ruling out false positives due to artifacts or deamination events. Next-generation sequencing (NGS) may offer a suitable alternative to conventional sequencing, however, substantially greater costs, complexity, and longer assay time render it an unnecessary option for many disease types with few distinct molecular markers or for monitoring patients on therapy for emerging resistance mutations. Furthermore, high false positive rates when detecting variants with mutant allele frequencies of less than 5% can pose a problem for amplicon-based NGS^5,6.

Here we demonstrate the increase in sensitivity achieved by WTB-PCR in screening for mutations in the myeloid differentiation factor 88 gene as described by Albitar et al.³ MYD88 mutations are important diagnostic and prognostic factors in Waldenström's Macroglobulinemia (WM), IgM monoclonal gammopathy of unknown significance (IgM-MGUS), splenic marginal zone lymphoma (SMZL), and diffuse large B-cell lymphoma (DLBCL). MYD88 mutations are found in almost all cases of WM and approximately 50% of patients with Immunoglobulin M (IgM)-secreting MGUS. Contrastingly, MYD88 mutations are found in only 0-6% of patients with SMZL and are absent in multiple myeloma^7,8. Because overlapping morphologic, immunophenotypic, cytogenetic, and clinical characteristics between WM and SMZL or IgM-multiple myeloma can often complicate differential diagnoses, the presence of a MYD88 mutation may serve as a useful identifying factor⁹. MYD88 mutations have also been associated with greater disease burden in patients with DLBCL and poor overall survival following therapy^7,10. Additionally, because MYD88 mutations are more frequently found in activated B-cell-like (ABC) DLBCL than germinal center B-cell-like (GCB) DLBCL or primary mediastinal B-cell lymphoma (PMBL), MYD88 mutation status may serve as a surrogate marker for the ABC subtype^11,12.

The detailed protocol provided here serves as a template from which new assays can be developed or most existing sequencing assays can be easily adapted to accurately detect low frequency mutations in various sample types. The approach can also be used for monitoring and detecting resistant mutations that may develop in tumors or even bacteria that may develop while patients are being treated with targeted therapy or antibiotics. Furthermore, it addresses and remedies many of the issues associated with mutation enrichment, particularly in formalin-fixed paraffin-embedded (FFPE) tissue.

Protocol

Ethics Statement: All testing of human samples was performed after obtaining Institutional Review Board (IRB) approval.

1. DNA Extraction from FFPE Tissue, Peripheral Blood, and Bone Marrow Aspirate

1. For bone marrow FFPE tissue with DNA FFPE extraction kit
   1. Begin with FFPE tissue from unstained slides (5 - 10 sections at 5 - 10 µm thickness).
      NOTE: If beginning with tissue shavings, use 3 - 6 sections at 5 - 10 µm thickness and skip to step 1.1.6.
   2. Place slides in a slide basket and prepare four wash reservoirs (two for Xylene and two for 100% Alcohol). Add a minimum volume of 600 mL of solution to each basket.
   3. Deparaffinize the slides by doing a 5 min xylene wash in the first tray. Transfer the slides to the second xylene wash reservoir for an additional 5 min.
   4. After deparaffinization, wash the slides with 100% alcohol for 5 min. Transfer the slides to the second alcohol wash reservoir for an additional 5 min.
   5. Allow the slides to dry completely under a hood before scraping them with a razor blade into a microcentrifuge tube.
      NOTE: If an H&E slide with tumor region indicated is available, align the slides and scrape only the tumor region.
   6. Proceed with instructions from manufacturer handout included in kit.
2. For BM aspirate and Peripheral Blood
   1. Extract according to the manufacturer's instructions handout with these specifications.
      NOTE: There are numerous commercially available kits and methods for DNA extraction from FFPE tissues and cells. Practically, all can be used. Use a method that is established in the laboratory.
      1. Use 200 µL peripheral blood (PB) or 100 µL BM + 100 µL PBS.
      2. Use 4 µL RNase A stock solution.
      3. Elute with 100 µL elution buffer.
3. DNA Quantification
   1. Measure DNA concentrations using a spectrophotometer ensuring a 260 nm/280 nm ratio of approximately 1.8 (for pure DNA). If the ratio is appreciably lower, it may indicate the presence of protein, phenol, or other contaminants that may interfere with downstream applications.
   2. Adjust DNA concentrations to approximately 50 - 100 ng/µL with appropriate elution buffer.

2. Wild-Type Blocking PCR

1. Primer Design
   1. Design and obtain primers for genes of interest according to previously described general guidelines of PCR primer design¹³. Include M13-forward and reverse universal sequencing primers in PCR primers.
      NOTE: The MYD88 assay was developed to amplify exon 5 of MYD88 (G259 - N291) and 110 nucleotides located in the 5' intronic region to cover the splice site and part of intron 4. The forward and reverse primers were designed with a 5'-M13 sequence (M13-forward: tgt aaa acg acg gcc agt; M13-reverse: cag gaa aca gct atg acc) to allow for annealing of complementary sequencing primers (see Materials Table).
2. Locked Nucleic Acid Oligonucleotide Design
   1. Design the blocking oligonucleotide to be approximately 10 - 15 bases in length and complementary to the WT template where mutant enrichment is desired.
      NOTE: A shorter oligo will improve mismatch discrimination. To achieve high target specificity, it is important not to use too much blocking nucleotides since this can result in a very "sticky" oligonucleotide. The content and the length and blocking nucleotide should to be optimize according to the specific nucleotide that need to be blocked. It is important to achieve high binding specificity without compromising secondary structure and risking self-complementarity.
   2. Design the blocking oligo to have a T_m 10 - 15 °C above the extension temperature during thermocycling. Adjust the T_m by adding or removing blocking bases or by substituting blocking bases for DNA while avoiding long stretches (3 - 4 bases) of blocking G or C bases.
      1. Navigate to the oligo tools website (see Table of Materials) and select the "blocking Oligo Tm Prediction" tool.
      2. Paste the sequence of the WT template that is desired to be blocked into the "Oligo Sequence" box. Add "+" in front of bases to indicate blocker bases.
      3. Select the "Calculate" button to determine the approximate T_m of the DNA-Blocker hybrid.
   3. Design the blocking oligo to avoid secondary structure formation or self-dimers.
      1. Navigate to the oligo tools website (see Table of Materials) and select the "blocking Oligo Optimizer" tool.
      2. Paste the sequence of the WT template that is desired to be blocked into the box. Add "+" in front of bases to indicate blocking bases.
      3. Select the two boxes for "Secondary Structure" and "Self Only". Press the Analyze button to review the oligo for potentially troublesome secondary structures or self-dimers.
         NOTE: The scores represent very rough estimates of the T_m's of the secondary structures and self-dimers. Lower scores are therefore optimal and can be achieved by limiting blocker-blocker pairing. The blocking oligonucleotide for MYD88 [MYD88blocker (TCAGA+AG+C+G+A+C+T+G+A+T+CC/invdT/)] was designed to cover amino acids Q262-I266 and features a 3'-inverted dT to inhibit both extension by DNA polymerase and degradation by 3'-exonuclease. This specific blocking oligo is a 17mer with 10 blocking bases which are denoted as "+N". The remaining 7 bases are ordinary DNA nucleotides.
3. WTB-PCR Setup and Thermocycling
   1. Prepare a WTB-PCR master mix (MMX) using 2.5 µL PCR reaction buffer 10x w/ 20 mM MgCl₂, 250 µM dNTPs, 0.2 µM forward and reverse primer, 1.2 µM MYD88blocker, and DNAse, RNAse-free, ultra-pure H₂O to create a final solution volume of 21.75 µL per reaction.
      NOTE: Working concentrations are for the final reaction volume of 25 µL (after addition of DNA template and polymerase). Conventional PCR MMX is prepared by simply omitting the addition of blocker. All protocol steps remain identical for both WTB and conventional PCR. This can be used in validation and determining enrichment achieved by addition of blocker.
      1. When calculating the amount of MMX to prepare make sure to account for 3 additional reactions (positive and negative controls and a non-template control to check for contamination) and at least 1 additional reaction for pipetting error.
   2. Vortex MMX thoroughly. The MMX can be stored at -80 °C for up to a year.
   3. Add 0.25 µL Taq DNA polymerase per reaction to the MMX and invert to mix. Once the polymerase has been added to the MMX, keep it on ice.
   4. Put a new 96 well PCR plate on a cold plate and pipette 22 µL of MMX with polymerase to each reaction well.
   5. Add 3 µL genomic DNA (50 - 100 ng/µL to each of the wells containing the MMX with polymerase.
   6. Seal the plate and centrifuge briefly to ensure the solution reaches the bottom of each well.
   7. Load the PCR plate on a thermocycler.
      1. Set the thermocycling conditions for the WTB-PCR reaction as follows: initial denaturation at 95 °C for 6 min; 40 cycles of denaturation at 95 °C for 30 s, primer annealing at 56 °C for 30 s, and extension at 72 °C for 1 min 20 s; final extension at 72 °C for 10 min.
         NOTE: Best practice involves completing the remainder of the protocol in a physically separate area to avoid amplicon contamination in future setups.
4. Perform gel electrophoresis according to previously described protocols¹⁴ in order to determine if WTB-PCR was successful and to confirm a lack of amplification in the non-template control.
5. Purification of PCR Product by magnetic beads
   1. Remove magnetic beads from 4 °C and bring to room temperature.
   2. Transfer 10 µL PCR product to a new PCR plate.
   3. Vortex the magnetic beads vigorously to fully resuspend the magnetic particles.
   4. Add 18 µL magnetic beads to each well on the new plate. Pipette mix 10 times.
   5. Incubate at room temperature for 5 min.
   6. Place the PCR plate onto the side-skirted magnet plate for 2 min to separate beads from solution.
   7. Aspirate the supernatant with a multichannel pipette, avoiding the bead pellet.
   8. Dispense 150 µL 70% ethanol into each well and incubate at room temperature for at least 30 s. Aspirate the ethanol with a multichannel pipette and discard tips. Repeat this wash procedure once more.
   9. Using a 20 µL multichannel pipette aspirate the remaining ethanol from each well and discard tips.
   10. Once wells are dry (~ 10 min), remove from magnet plate and add 40 µL nuclease free water to each well and pipette mix 15 times or vortex gently.
   11. Incubate at room temperature for 2 min.
   12. Place the PCR plate onto magnet plate for 1 min to separate beads from solution.
   13. Transfer 35 µL of purified product to a new PCR plate. This is purified PCR product is now ready to proceed with bi-directional sequencing or can be stored at -20 °C until needed.
       NOTE: When developing a new assay, quantification of purified PCR product may be required. 1 - 3 ng/µL amplicon DNA is optimal for bi-directional sequencing. If concentration is consistently low, increase PCR product and magnetic beads proportionally (1:1.8). If concentration is consistently high, elute with a greater volume of water.

3. Sequencing of WTB-PCR Product

1. Bi-Directional Sequencing
   NOTE: The following protocol is a modified form of the manufacturer's instructions that has been optimized to use fewer reagents.
   1. Prepare forward and reverse sequencing reaction mixes with 0.25 µL Ready Reaction mix, 1.88 µL Sequencing Buffer, 1.78 µM M13-F or M13-R sequencing primers and add and DNAse, RNAse-free, ultra-pure H₂O to create a final solution volume of 9 µL per reaction. This sequencing reaction mix can be stored at -20 °C for up to a year.
   2. Pipette 9 µL of the forward sequencing reaction mix into each well of a new 96 well PCR plate for the forward sequencing reaction. Repeat on a separate plate for the reverse sequencing reaction.
   3. Add 1 µL of the purified WTB-PCR product to each well on both the forward and reverse plates.
   4. Seal the plate and centrifuge briefly to ensure the solution reaches the bottom of each well.
   5. Load the PCR plate on a thermocycler.
      1. Set the thermocycling conditions for the sequencing reaction as follows: 96 °C for 1 min; 30 cycles of 96 °C for 10 s, 50 °C for 5 s, 60 °C for 4 min; hold at 4 °C until ready for purification.
2. Purify Sequencing Products
   1. Prepare a fresh 1:25 solution of 3 M sodium acetate (pH 5.2) and 100% ethanol.
   2. Prepare a fresh 70% ethanol solution.
   3. Add 30 µL of sodium acetate/100% ethanol solution to each well of both the forward and reverse sequencing plates and pipette mix 5 times.
   4. Reseal the plate and incubate in the dark at room temperature for 20 min.
   5. Centrifuge the plate at 2,250 x g for 15 min.
   6. Remove the plate sealer and invert over a waste container.
      NOTE: Invert only once or the pellet may loosen from the well bottoms.
   7. Place the inverted plate on a clean paper towel and centrifuge at 150 x g for 1 min.
   8. Add 150 µL 70% ethanol to each well.
   9. Reseal the plate and spin at 2,250 x g for 5 min.
   10. Repeat steps 3.2.6 and 3.2.7.
   11. If the wells are not completely dry, allow them to air dry at room temperature. Make sure the samples are protected from light while drying.
   12. Add 10 µL Formamide to each well and pipette mix 10 times. Reseal the plate.
   13. Denature on thermocycler at 95 °C for 3 min followed by 4 °C for 5 min.
   14. After denaturing, replace the plate sealer with a septa and sequence on sequencing platform according to manufacturer's instructions¹⁵.

4. Analysis of Sequencing Results

1. Visualize sequencing traces using a sequencing data analysis software according to manufacturer's instructions¹⁶ and align to the respective reference sequences. Align MYD88 to NCBI reference sequence "NM_002468".

Results

A conceptual overview of WTB-PCR during extension is presented in Figure 1. Because a single nucleotide mismatch in the blocker-DNA hybrid greatly decreases its melting temperature (ΔT_m=20 - 30 °C), amplification of the WT allele is blocked while mutant template DNA is free to complete extension¹⁷. In this manner, mutant DNA is amplified exponentially while WT DNA is amplified linearly.

Discussion

The WTB-PCR assay described here uses a generic set of primers with a blocking oligo designed to block amplification of WT DNA during extension (Figure 1). The WTB-PCR product is then sequenced for mutational analysis. The utility of WTB-PCR/Sanger lies in its simplicity, high-sensitivity, and high-throughput. Using the guidelines described here, most existing Sanger based assays can be simply modified via the addition of a blocking oligonucleotide to greatly increase sensitivity. In the example assay pr...

Materials

Name	Company	Catalog Number	Comments
1.5 or 2 mL Safe-Lock microcentrifuge tubes	Eppendorf	05-402-25
100% alcohol	VWR	89370-084	Histology grade; 91.5% Ethanol, 5% Isopropyl alcohol, 4.5% Methyl alcohol
3730XL sequencer	ABI		or equivalent
Agencourt AMPure XP	Beckman Coulter	A63881	For magnetic bead PCR purification
Aluminum sealing foils	GeneMate	T-2451-1	For PCR and cold storage
BigDye Terminator v3.1 Cycle sequencing kit	Life Technologies	4337455	For bi-directional sequencing. With 5x Sequencing Buffer
Centrifuge 5804 Series	Eppendorf		A-2-DWP rotor (for PCR plate)
Cold plate for 96 well plates	Eppendorf	Z606634
DNAse, RNAse-free, ultra-pure water
dNTPs (100 mM)	Invitrogen	10297-117
DynaMag-96 Side-Skirted Magnet	Thermo Fisher Scientific	12027	For use in PCR Purification.
Ethanol Absolute	Sigma	E7023	200 proof, for molecular biology
Exiqon website Oligo Tools			www.exiqon.com/oligo-tools
FastStart Taq DNA polymerase (5 U/µL)	Roche	12032937001	With10x concentrated PCR reaction buffer, with 20 mM MgCl2
Gel electrophoresis apparatus			2% agarose gel
GeneRead DNA FFPE extraction Kit	Qiagen	180134	Contains uracil DNA glycosylase necessary for reducing sequencing artifacts
Hi-Di Formamide	ABI	4311320	For sequencing.
LNA oligonucleotide	Exiqon	500100	5'-TCAGA+AG+C+G+A+C+T+G+A+T+CC/invdT/ (+N = LNA bases)
M13-F Sequencing Primer	ABI		5'-tgt aaa acg acg gcc agt
M13-R Sequencing Primer	ABI		5'-cag gaa aca gct atg acc
Mastercycler Pro S Thermocycler	Eppendorf	E950030020
Microcentrifuge Model 5430	Eppendorf		FA-45-30-11 rotor (for 1.5/2 mL microcentrifuge tubes)
NanoDrop 2000 Spectrophotometer	Thermo Fisher Scientific
PCR forward primer	IDT		5'-tgt aaa acg acg gcc agt TGC CAG GGG TAC TTA GAT GG
PCR reverse primer	IDT		5'-cag gaa aca gct atg acc GGT TGG TGT AGT CGC AGA CA
PCR plates	GeneMate	T-3107-1
Pipettors			20, 200, 1,000 µL
Plate septa, 96 well	ABI	4315933
QIAamp DNA Mini Kit	Qiagen	51304	For BM aspirate and peripheral blood
SeqScape Sortware v3.0	ABI	4474978	For sequencing analysis
Slide basket
Sodium Acetate (3 M, pH 5.2)	Sigma	S7899
Sterile filtered pipette tips			20, 200, 1,000 µL
Thermomixer C	Eppendorf	5382000023
Vortex genie	Scientific Industries	SI-0236
Wash reservoir			~ 1,000 mL
Xylene	VWR	89370-088	Histology grade