Identifying Mutations by High Resolution Melting in a TILLING Population of Rice

Summary

In this article, we present the protocol that is described as high-resolution melting analysis (HRM)-based Target Induced Local Lesions In Genomes (TILLING). This method utilizes fluorescence changes during the melting of the DNA duplex and is suitable for high-throughput screening of both insertion/deletion (Indel) and single base substition (SBS).

Abstract

Target Induced Local Lesions In Genomes (TILLING) is a strategy of reverse genetics for the high-throughput screening of induced mutations. However, the TILLING system has less applicability for insertion/deletion (Indel) detection and traditional TILLING needs more complex steps, like CEL I nuclease digestion and gel electrophoresis. To improve the throughput and selection efficiency, and to make the screening of both Indels and single base substitions (SBSs) possible, a new high-resolution melting (HRM)-based TILLING system is developed. Here, we present a detailed HRM-TILLING protocol and show its application in mutation screening. This method can analyze the mutations of PCR amplicons by measuring the denaturation of double-stranded DNA at high temperatures. HRM analysis is directly performed post-PCR without additional processing. Moreover, a simple, safe and fast (SSF) DNA extraction method is integrated with HRM-TILLING to identify both Indels and SBSs. Its simplicity, robustness and high throughput make it potentially useful for mutation scanning in rice and other crops.

Introduction

Mutants are important genetic resources for plant functional genomics research and breeding of new varieties. A forward genetics approach (i.e. from mutant selection to gene cloning or variety development) used to be the main and sole method for the use of induced mutations about 20 years ago. The development of a novel reverse genetics method, TILLING (Targeting Induced Local Lesions In Genomes) by McCallum et al.¹ opened a new paradigm and it has since been applied in a great number of animal and plant species². TILLING is particularly useful for breeding traits that are technically difficult or costly to be determined (e.g., disease resistance, mineral content).

TILLING was initially developed for screening point mutations induced by chemical mutagens (e.g., EMS^1,3). It includes the following steps: the establishment of a TILLING population(s); DNA preparation and pooling of individual plants; PCR amplification of target DNA fragment; heteroduplexes formation by denaturation and annealing of PCR amplicons and cleavage by CEL I nuclease; and identification of mutant individuals and their specific molecular lesions^3,4. However, this method is still relatively complex, time consuming, and low-throughput. To make it more efficient and with higher throughput, many modified TILLING methods have been developed, such as deletion TILLING (De-TILLING) (Table 1)^{1,3,5,6,7,8,9,10,11,12}.

HRM curve analysis, which is based on fluorescence changes during the melting of the DNA duplex, is a simple, cost-effective, and high-throughput method for mutation screening and genotyping¹³. HRM has already been widely used in plant research including HRM based TILLING (HRM-TILLING) for screening SBS mutations induced by EMS mutagenesis¹⁴. Here, we presented detailed HRM-TILLING protocols for screening of mutations (both Indel and SBS) induced by gamma (γ) rays in rice.

Protocol

1. Preparations

1. Development of γ-rays mutagenized populations
   1. Treat about 20,000 dried rice seeds (with moisture content of ca. 14%) of a japonica rice line (e.g., DS552) with ¹³⁷Cs gamma rays at 100 Gy (1 Gy/min) in a γ irradiation facility (e.g., gamma cell).
      NOTE: Seeds used for treatment should have a high viability (e.g., with a germination rate >85%). The irradiation dose for indica rice could be increased to 150 Gy.
   2. Sow the irradiated seeds after germination on a seedling bed and transplant seedlings individually to a paddy field and grow into the M₁ population.
      NOTE: Direct seeding of M₁ plants sparsely could also be applied to save labor cost. Prevent outcrossing of M₁ plants with other rice varieties using physical or biological isolation means.
   3. Bulk-harvest M₂ seeds from the M₁ plants, with 1-2 seeds from each M₁ panic, to form an M₂ population.
      NOTE: In practice and for simplicity, harvest all seeds of M₁ plants and after fully mixing, a portion is sampled to form an M₂ population.
   4. Soak about 5,000 M₂ seeds in water for 24 (indica rice)-36 (japonica rice) h at room temperature. Then let the seeds to germinate at 37 °C for 2 days on moist filter paper in Petri dishes.
      NOTE: More M₂ seeds can be germinated for analysis to increase the probability of identifying mutants.
   5. Place the germinated seeds to seeding panels with small holes and grow them hydroponically for 3-4 weeks in a culture solution modified from Yoshida et al.¹⁵ in a glasshouse with a 12 h photoperiod [daytime (30±2) °C and night (24±2) °C].
2. Sampling of leaf tissues: Cut one disk (Φ ~2 mm) from the fully extended leaf of each seeding at the same position using a hole puncher.
3. Preparation of DNA extraction solutions
   1. Buffer A: Add 2 mL of 5 M NaOH and 10 mL of 20% Tween 20 to make the final volume of 50 mL. Freshly prepare buffer A before DNA extraction.
   2. Buffer B: Add 20 mL of 1 M Tris-HCl (pH 8.0) and 80 μL of 0.5 M EDTA to make a final volume of 100 mL.
4. PCR primers
   1. Primers for HRM analysis: Design primers for amplification of the target sequence using software (e.g., Primer Premier5) and synthesize by a commercial company.
      NOTE: Because HRM is less applicable to the analysis of long fragments, and hence, the amplicons should be less than 400 bp. Fragments with too high (>75%) or too low (<25%) GC content are also not good for HRM analysis.
   2. Primers for quality control: Use the 24 SSR markers distributed on the 12 rice chromosomes from Peng et al.¹⁶.

2. DNA Extraction

1. Place 4 leaf discs into each well of a 96-well PCR plate, add buffer A solution (50 mL/well). Freeze the plate in a -80 °C freezer for 10 min.
2. Defreeze the plate at room temperature, and then incubate at 95 °C for 10 min.
3. Add 50 μL of buffer B to each well and mix well by vortexing.
4. Centrifuge the plate for 1 min at 1,500 x g. The supernatant is ready for PCR.

3. PCR Amplification

1. PCR Optimization
   1. Add the following reagents to each PCR well: 1 μL of DNA (the supernatant in step 2.4), 5 μL of 2x Master Mix (containing 2x PCR buffer, 4 mmol/L MgCl₂, 0.4 mmol/L 2'-deoxyribonucleoside triphosphates (dNTPs), and 50 U/mL Taq DNA polymerase, 0.2 μL each of 10 μmol/L primers, and make a final volume up to 10 μL using nuclease free water.
   2. Use a gradient-capable thermal block to determine the optimal annealing temperature for each target fragment by using the following PCR program: 5 min at 94 °C, followed by 40 cycles of 30 s at 94 °C, 30 s at 52-62 °C (gradient temperature), and 30 s at 72 °C, with a final extension at 72 °C for 8 min and a hold at 16 °C.
   3. Examine amplicons on 1% agarose gels for determination of the optimal annealing temperature.
      NOTE: An optimal annealing temperature should enable specific amplification of the target fragment, without nonspecific amplification and primer dimerization.
2. PCR for HRM analysis
   NOTE: HRM compatible plates and DNA of M₂ seedlings extracted as described in step 2.4 are used for PCR.
   1. Perform PCRs in a final volume of 10 μL, with 1 μL of DNA (supernatant), 5 μL of 2x Master Mix, 0.2 μL each of 10 μmol/L primers, and 1 μL of 10x fluorescence dye. In each plate include one wild type (WT) parented sample and one negative (without DNA) control.
   2. Add a drop of mineral oil to each well to prevent evaporation.
   3. Seal the plate with adhesive film and centrifuge at 1,000 x g for 1 min.
   4. Run the PCR using the optimized annealing temperature.
      NOTE: In a few cases, fluorescence dye may affect PCR amplification, hence the dye is added after completion of PCR. In such cases, the dye is incorporated into DNA strands by post-PCR denaturing and annealing.

4. HRM Scanning and Mutation Confirmation

1. Remove the adhesive film from plate and insert the plate into an HRM machine.
2. Select New Run from the file menu or press the Run button at the top of the screen.
3. Specify the starting and ending temperatures for the melt from 55 °C to 95 °C.
   NOTE: After the first run, the range of melt temperature can be determined for a particular fragment; hence, the melt temperature can be adjusted for subsequent analysis of the same fragment to save time.
4. Select samples for high resolution melting analysis, exclude samples similar to the negative control.
5. Normalize the melting curves to have the same beginning and ending fluorescence. Visually confirm that the Lower Min and Lower Max temperature cursors are in a region of the curves.
6. Keep the ∆F (difference of fluorescence) level at the default setting of 0.05.
7. Select the Common versus Variant from the Standards selection list and choose the Normal sensitivity.
8. Use the WT as the control; samples with a ∆F value of ≥0.05 from WT are considered to contain mutant plant.
9. Identification and confirmation of mutant plants.
   1. Identify the four plants, of which each mutant pool was made.
   2. Extract DNA from each of these plants using a CTAB method according to Allen et al.¹⁷ with some modifications.
      NOTE: DNase-free RNase and NaAc are not used when extracting DNA using the CTAB method described by Allen et al.¹⁷, as the DNA quality is good enough for further PCR amplification. Adjust the DNA to a final concentration of ~25 ng/μL after quantification using a spectrophotometer.
   3. Amplify the target fragment using PCR primers and program the same as for HRM analysis.
   4. Identify specific molecular lesions by Sanger sequencing of the amplicons.
      NOTE: Once the completion of PCR amplification, send the amplicons to a company for Sanger sequencing. The molecular lesion can be identified by comparing the sequences between the M₂ plant and the WT.

5. Quality Control of Selected Mutants with Molecular Markers

1. Perform PCR for the 24 SSR markers in a final volume of 10 μL with 25 ng of genomic DNA (extracted using the CTAB protocol), 5 µL of 2x master mix, 0.2 µL of each of 10 μM SSR primers.
2. Run PCR using the following program: 5 min at 94 °C, followed by 30 cycles of 30 s at 94 °C, 30 s at 55 °C and 30 s at 72 °C, with a final extension at 72 °C for 7 min.
3. Separate the amplicons on 8% polyacrylamide gels and reveal polymorphism of amplified fragments by silver staining¹⁸.
4. Compare the SSR haplotypes of selected variants and the WT.
   NOTE: Induced mutants often have SSR haplotypes identical to their WT, if more than one SSR markers are different between a variant and the WT, the variant is high likely to be a genetic contaminant (e.g., a mixture or outcrossed plant) rather than an induced mutant¹⁸.

Results

HRM Scanning and Analysis

In total, 1,140 pooled DNA samples from 4,560 M₂ seedlings were produced and subjected to PCR amplification. Two fragments with the size of 195 bp and 259 bp were amplified for OsLCT1 and SPDT, respectively (Table 2). Most samples had melting curves not significantly different from the WT (ΔF < 0.05). HRM curves significantly different from the WT (ΔF > 0.05) were grouped with color(s) di...

Discussion

TILLING has proved to be a powerful reverse genetic tool for identifying induced mutations for gene functional analysis and crop breeding. For some traits not easily observed or determined, TILLING with high-throughput PCR-based mutation detection can be a useful method to obtain mutants for different genes. HRM-TILLING method has been used in EMS-mutagenized populations of tomato¹², wheat¹¹ and grapevine²⁰ for mutation screening. In this paper, a si...

Materials

Name	Company	Catalog Number	Comments
2× Taq plus PCR Master Mix	Tiangen, China	KT201	PCR buffer, dNTP and polymerase for PCR amplification
96-well plate	Bio-rad, America	MSP-9651	Specific plate for PCR in HRM analysis
Mastercycler nexus	Eppendorf, German	6333000073	PCR amplification
LightScanner	Idaho Technology, USA	LCSN-ASY-0011	For fluorescence sampling and processing
CALL-IT 2.0	Idaho Technology, USA		For analysis of the fluorescence change
EvaGreen	Biotium, USA	31000-T	Fluorescence dye of HRM
Nanodrop 2000	Thermo Scientific, USA	ND2000	For DNA quantification