Generation of EMS mutagenized population

We developed TILLING populations in tetraploid and hexaploid wheat using EMS as a chemical mutagen. For the tetraploid population we mutagenized seeds of the Desert durum^® variety 'Kronos', which was developed by Arizona Plant Breeders from a male sterile population (selection D03–21). For the hexaploid TILLING population we used the Hard Red Spring common wheat breeding line 'UC1041+Gpc-B1/Yr36'. UC1041 is a short stature breeding line developed by the University of California from the cross Tadinia/Yecora Rojo. 'UC1041+Gpc-B1/Yr36' was developed later by backcrossing for six generations a 6BS chromosome segment from T. turgidum ssp. dicoccoides that carries the high grain protein gene Gpc-B1 [31] and the partial stripe rust resistance gene Yr36 [32]. The EMS concentrations used to mutagenize the populations were 0.7 to 0.75% (57 to 60 mM) for Kronos and 0.9 to 1.0% (73 to 80 mM) for 'UC1041+Gpc-B1/Yr36'. Similar EMS concentrations have been used previously to create TILLING population in wheat [18, 21]. Germination rates for EMS-treated seeds were ~50–60% (results not shown). We extracted DNA from single M₂ plants and collected their M₃ seeds to have independent and non-redundant mutations in our libraries. DNAs from a total of 1,368 M₂ (tetraploid) and 1,536 M₂ (hexaploid) plants were pooled in groups of four DNAs and organized into four 96-well plates for convenient screening (342 and 384 4× pools in the tetraploid and hexaploid populations, respectively). The tetraploid TILLING population is currently being expanded to 1,536 lines.

Development of genome specific primers

We characterized the TILLING libraries by screening for mutations in the two Starch Branching Enzyme II genes, SBEIIa and SBEIIb, for tetraploid wheat and for mutations in the Wheat Kinase Start (WKS) 1, WKS2 and SBEIIa genes in hexaploid wheat. WKS1 and WKS2 are single copy genes on chromosome arm 6BS [32], so there was no need to develop genome specific primers. SBEIIa and SBEIIb map to chromosome 2 and have homoeologous loci in each of the different wheat genomes. To screen for mutations in each of the homoeologous copies we designed primers specific for each copy taking advantage of polymorphic indels and SNPs between the different homoeologues. We designed primers complementary to intron sequences flanking the target exons and positioned approximately ~200-bp from the sequence of interest. The genome specificity of the SBEII primers was validated using nulli-tetrasomic lines of chromosome 2 (N2AT2B, N2BT2D, N2DT2A), DNAs from BAC clones from the A and B genomes obtained from a tetraploid BAC library [8] and Aegilops tauschii genomic DNA.

Wheat TILLING platform using the Non-denaturing Polyacrylamide Detection Method

The screening of a TILLING population includes three fundamental steps: an initial screen of DNA pools to identify those that contain mutant individuals; a second screen to identify the individual within each pool that contains a putative mutation; and lastly, confirmation of the individual mutations by sequencing of the PCR products. Although some modifications exist, such as directly sequencing all individuals from a positive pool, most TILLING approaches follow this general framework.

We developed a TILLING platform that uses a non-denaturing polyacrylamide detection method to perform the two rounds of screening. After digestion of mismatches in heteroduplexes with celery juice extract (CJE), the fragments resulting from dsDNA cuts at mismatched sites are separated in native polyacrylamide gels and visualized through the fluorescence of bound ethidium bromide. For the initial step, the screen identifies the 4× pools with mutant individuals (Figure 1A). Targets are amplified by PCR from the four genomic DNAs that comprise each positive pool (a, b, c, d in Figure 1B). PCR products are then combined in four two-fold pools, heated and annealed to achieve heteroduplex formation, and finally digested with CJE (Figure 1C). Depending on the banding pattern (Figure 1C), the mutation is assigned to one of the four individual DNAs (Figure 1D, only one banding pattern example is shown). The mutant individual is then sequenced and the identity of the mutation is established (Figure 1E). A more detailed description can be found in the Methods section.

This method provides an independent validation of the mutation, identifies its location within the target region, and determines which individual from the pool carries the mutation. The paired pooling (Figure 1B) is necessary to detect homozygous mutations in the M₂ plants since combining two samples allows heteroduplex formation and detection. This strategy also reduces the number of false positive errors as a true mutation should be observed in two separate gel lanes (Figure 1C).

In an M₂ population, 33% of the mutations are expected to be found in homozygous state. For both populations we had a slight bias towards homozygous mutations, 37% in hexaploid and 42% in tetraploid, although these percentages were not significantly different from the expected 33% (hexaploid χ² = 1.18, P = 0.28; tetraploid χ² = 3.09, P = 0.08). Sequencing also confirmed that over 99% of the mutations were G to A or C to T transitions as expected from alkylation by EMS, with only one exception in the SBEIIa A genome target which was a C to G transversion (T6-2312).

Using the CODDLe (Choose codons to Optimize the Detection of Deleterious Lesions) program [33], we predicted the effect of EMS mutations in the different amplicons. In the hexaploid library we identified a total of 186 mutations of which 40% were missense and 4.3% were truncations (nonsense or splice junction mutations). The predicted effects by CODDLE were 35% missense and 4.5% truncations, very close to the observed values (Figure 2). In the tetraploid screen we identified 93 mutations of which 28% were missense and 5.4% were truncations, whereas CODDLE predicted 22% missense mutations and 3.9% truncations. For both libraries, the distribution of silent, missense and truncation mutations were not significantly different from those predicted by CODDLE (tetraploid χ² = 2.37, P = 0.31; hexaploid χ² = 2.66, P = 0.26).

Several of the mutations in the WKS and SBE genes were identified in more than one independent individual. In the hexaploid library, 11.8% of the mutations were found in duplicate or triplicate (8 duplicate and 2 triplicate mutations, or 22 mutations in 186), which was higher (P < 0.01) than the expected 6.5% calculated using a Poisson distribution and the number of potential GC sites in the screened region. In the tetraploid library 24.7% of the mutations were found in more than one individual (8 duplicate, 1 triplicate, 1 quadruple, or 23 mutations in 93). These numbers were again higher (P < 0.01) than the expected 3.4% predicted by a Poisson distribution.

Comparison between LI-COR and Non-denaturing Polyacrylamide Detection Method

Laser detection of fluorescently labeled DNA fragments, using a LI-COR genotyping platform, is the most widely used detection method to screen TILLING populations for mutations. To evaluate an alternative to this detection method, we screened two regions in both WKS1 and WKS2 using a non-denaturing 3% polyacrylamide set-up (Figure 1) and compared the results with an established LI-COR platform.

TILLING of the same four WKS targets in 768 M₂ individuals from the hexaploid library revealed that these two methods detect comparable number of mutations (Table 4). We used the method of Greene et al. [34] to estimate mutation densities (cumulative length of sequence screened divided by total number of mutants). We adjusted the total length of each target as mutations in the regions closest to the primers are not readily detected by either method. For the LI-COR screen we subtracted 160-bp [34], whereas for the 3% polyacrylamide screen we subtracted 10% of the target region at each end. This value was determined empirically as mutations were only detected in the central 80% of the target sequence (effective target region) for all genes (Figure 3).

		LI-COR			Polyacrylamide/Ethidium bromide
Gene	Region	Sequence screened (kb)	Mut.	Mutation Frequency	Sequence screened (kb)	Mut.	Mutation Frequency
WKS1	Kinase	930.0	18	1/52 kb	842.3	20	1/42 kb
	START	852.5	28	1/30 kb	390.1a	15	1/26 kb
WKS2	Kinase	998.4	25	1/40 kb	897.0	25	1/36 kb
	START	1053.7	25	1/42 kb	941.3	14	1/67 kb
Total/mean		3835	96	1/39.9 kb	3071	74	1/41.5 kb

^a 384 M₂ plants were screened

Four target regions where examined in the same 768 M₂ plants from the hexaploid TILLING population and the total sequence screened was adjusted according to each method (see text).

Using the LI-COR detection technology we estimated a mutation frequency of 1 mutation per 40 kb (96 mutants in 3.84 Mb screened), whereas using the polyacrylamide setup we estimated 1 mutation per 41.5 kb (74 mutants in 3.07 Mb screened). Overall, each method detected two or three mutants not detected by the other method, but the majority of the mutations were detected by both. The exception to this was the WKS2 START domain target region in which the LI-COR screen identified eleven additional mutations (Table 4). Despite this, the average mutation frequencies for the four targets were almost identical. This suggests that when using fourfold pools (eight-fold dilution for mutations in heterozygous individuals), the polyacrylamide/ethidium bromide set-up has similar sensitivity to detect SNPs compared to the LI-COR platform, although there may be slight differences depending on the target region.

Characterization of the EMS mutagenized populations

The use of reverse genetic approaches to determine gene function will become increasingly important as large amounts of sequence information become available in wheat. In an effort to address this, we created EMS-induced TILLING libraries in tetraploid and hexaploid wheat. We use the tetraploid TILLING population to generate mutants for basic research projects because it is easier and faster to generate complete null mutants. A single generation of crosses between A and B genome mutations, followed by selection of homozygous double mutants in the F₂ populations is sufficient to generate null mutants. However, when a targeted mutant has important commercial applications (e.g. the sbeIIa mutants with predicted high amylose phenotype [35]) we screen the hexaploid TILLING population for mutations, because hexaploid wheat represents most of the wheat grown around the world (~95%) [36].

Characterization of these populations through the screening of several targets revealed mutation densities of at least one mutation per 49.4-kb and 68-kb in the hexaploid and tetraploid libraries, respectively. These mutation densities are lower than those found by Slade et al. [18] in wheat using similar EMS concentrations (one mutation per 24-kb and 40-kb screened in hexaploid and tetraploid libraries, respectively). The difference observed in these two studies is likely dependent on the different GC content of the target regions employed in the two studies, because EMS mutagenesis acts predominantly on GC residues [34]. Slade et al. [18] characterized their libraries by screening for mutants in the waxy genes. The regions included in their work have an unusually high GC content [Wx-A1 (56.4%), Wx-B1 (59.6%), Wx-D1 (55.4%)] whereas regions targeted in our study have an average GC content of 37.3% and 38.8% in the tetraploid and hexaploid targets, respectively. These values must be taken into account when estimating mutation densities because they represent the maximum number of mutations that can be found in those particular targets. For example, adjusting our reported mutation densities to the average GC content in Slade et al. [18] (55.9% hexaploid, 58.0% tetraploid) would yield new densities of one mutation per 34-kb and 44-kb for the hexaploid and tetraploid libraries, respectively. For future studies involving species where EMS mutagenesis is limited to GC>AT changes, it would be beneficial to report mutation densities corrected for a 50% GC content, or specify the GC content of the target regions, to allow for more meaningful comparisons. Applying this criterion, our mutation densities would be one mutation per 38-kb and 51-kb for the hexaploid and tetraploid libraries, respectively, in a target with 50% GC content.

Independent of the GC content, our reported mutation densities were lower than those reported by Slade et al. [18], by 43% in the hexaploid and 9% in the tetraploid population. The response of different hexaploid genetic backgrounds to EMS could account for some differences between the hexaploid libraries. This explanation cannot be applied to the tetraploid libraries since both studies used the same cultivar Kronos. Slight differences in EMS concentration and treatment conditions, environmental effects and experimental differences in the detection methods for both studies (for both LI-COR and 3% non-denaturing polyacrylamide) could account for the remaining variation.

Wheat is especially well suited for TILLING because of the tolerance of recently evolved polyploid species to high mutation densities [36]. The vast majority of greenhouse grown plants was fertile and displayed no apparent mutant phenotype. The mutation frequencies for wheat reported here and by Slade et al. [18] are five to ten times higher than mutation rates found in diploids such as barley, pea and Arabidopsis [25, 37, 38]. This high mutation frequency facilitates the identification of large allelic series in target genes using relatively small TILLING populations. For example, by screening 1,536 lines for a 1625-bp target region (1300-bp effectively screened, 50% GC content) we would expect to recover approximately 52 mutant alleles in the hexaploid library and 39 mutant alleles in the tetraploid library. Analysis of the mutations obtained in this study confirmed that the frequencies predicted by CODDLE were accurate and can be used to estimate the expected proportion of the different types of mutations to be recovered.

In an average TILLING fragment, truncation mutations are expected in 4 to 5% of the cases. Therefore, the large number of mutant alleles expected when TILLING an effective 1.3 kb region (1625-bp total target, 50% GC content) provides a high probability (>90% in hexaploid and 84% in tetraploid; P = [1-(1-0.045)^{number of alleles}]) of obtaining at least one truncation mutation. In our screen using targets with lower GC content (<40%) we found truncations for 71.4% and 75% of the targets in the hexaploid and tetraploid libraries, respectively. This probability will vary according to GC content and can be improved by increasing the size of the target region if the gene is large enough (for example TILLING two regions of the same gene).

Strategies for TILLING in polyploid genomes

The high probability of identifying truncation mutants is very important in a polyploid species, such as wheat, where the phenotype of a single mutant may be masked by the wild-type homoeologue present in another genome. Because of gene redundancy, it is generally necessary to cross single mutants in the A and B genome homoeologues to obtain a functional knockout in tetraploid wheat or create the triple A/B/D mutant in hexaploid wheat. Employing missense mutations in these lengthy genetic schemes is risky because if one of the mutations is not effective, it may be sufficient to limit the effect of the combined mutations on function. For the A and B genomes, the search for nonsense or splice junction mutations can involve both tetraploid and hexaploid TILLING populations, since mutations can be transferred by crossing. Hybridization of Kronos and UC1041 produces a pentaploid F₁, which can be used as a female in subsequent backcrossing until fertility is restored.

Bioinformatics algorithms such as SIFT [39] and PSSM can be used to prioritize mutants for phenotypic evaluation, as reported before for WKS1 [32]. All five mutations with significant PSSM and SIFT scores were loss-of-function mutants of the resistance gene that led to susceptibility to the causal agent of stripe rust, Puccinia striiformis f. sp. tritici (Table 2). The only mutant line that remained resistant was T6-89 that had a non-significant SIFT score (0.46) and a borderline significant PSSM score (10.4). Despite this successful example for the use of SIFT and PSSM, the decision of using a missense mutation should be weighed against the amount of time and work that would be invested in producing double and triple mutants. The optimum strategy will depend on the objective and the gene being studied, as in some cases homoeologues are naturally deleted (as was the case for WKS1) or are not expressed.

The high mutation density in our libraries also implies that any given individual is predicted to carry between 260,000 (tetraploid) and 415,000 (hexaploid) mutations. Since most of the wheat genome (>83%) is represented by highly repetitive elements, and likely less than 3% of the wheat DNA encodes for genes (assuming a similar gene space per genome as Arabidopsis), most of the mutations will be outside the genes. Even after correcting for repetitive regions, for coding sequence space within a gene, and for the proportion of mutations that result in missense or truncations, each individual from the TILLING population is expected to carry thousands of missense mutations and hundreds of truncations. A simple way to reduce this large amount of background mutations is to backcross the mutants to the non-mutagenized recurrent parent for two to three generations. These backcross generations are essential when the mutations are being used for wheat breeding because the background mutations can reduce the average performance of the populations generated directly from crosses with the original mutants. The mutant SNP can be tracked in the backcrossing scheme by direct sequencing of the genome-specific amplicon in each generation. Alternatively, many SNPs lead to polymorphisms in restriction sites which can be used to develop Cleavage Amplified Polymorphism (CAPs) markers. Alternatively, derived CAP (dCAP) markers can also be designed [40]. Ultimately, the most effective strategy will depend on the costs of sequencing and restriction enzymes for each lab.

For research projects with a clear phenotype, selecting for sister lines homozygous for the presence and absence of the mutation is an effective strategy. Sister lines with and without the mutations share many of the same background mutations, thus serving as a better control than the wild-type line with no background mutations. This approach is especially powerful when multiple sets of independent sister lines are examined [32] as the probability of finding mutations by chance in a linked gene is extremely low. For example, with the mutation densities of the hexaploid population (1 mutation per 38 kb), the probability of finding at least one an amino acid change in any 1.5 kb coding region is 2.5% ({1- [(1-(1/38000)]¹⁵⁰⁰}*0.66); as 66% of GC>AT codon changes are non-synonymous]. If two or three independent lines were examined, then this probability drops dramatically (P < 0.0007 and P < 0.00002, respectively).

Primer design in polyploid species

Primer design is an important aspect of TILLING in polyploid species. Genome specificity needs to be combined with a high yielding PCR product for proper mutant detection. The first step in designing genome specific primers is the sequencing of the different homoeologous copies. These sequences can be obtained for highly expressed genes by a bioinformatics characterization of available wheat ESTs. Alternatively, genome sequences can be generated by screening the BAC libraries and sequencing from individual BACs or by sequencing the diploid donors of the different wheat genomes. We routinely use T. urartu for the A genome, Ae. tauschii for the D genome (accession AL8/78 closely related to the D genome of wheat) and T. speltoides as the best approximation to the B genome. If better sequences for the B genome are required, a fast strategy is to clone and sequence several clones from PCR products obtained from tetraploid wheat.

The optimal target regions were defined by the CODDLE program using the following criteria: a) mutations close to primers (~10% of target sequence) are not readily detected, particularly in large amplification products b) maximize exons and/or intron-exon splice junctions and c) maximize regions encoding for conserved domains within the protein. Primers are usually designed in the introns or 5' or 3' UTRs flanking the target exons as these regions are more polymorphic (important for genome specificity).

Different strategies can be used to generate genome-specific primers (Figure 4) [41, 42]. If one of the primers can overlap a unique in/del or multiple intergenomic SNPs, this is usually sufficient to generate genome specificity (e.g. SBEIIa A and D genomes, SBEIIb B genome). In other cases, where there is lower polymorphism between homoeologues, both primers can be designed such that the first nucleotide from the 3' end of the primers aligns to genome-specific SNPs (e.g. SBEIIb A genome). In these cases, increased specificity can be attained by introducing a mismatch in the primer at the third or fourth position from the 3' end. Although this generates a mismatch between the target sequence and the primer (at the third or fourth position from the 3' end), the two mismatches with the other homoeologues increase the probability of genome-specific amplification. These strategies can be combined as in the SBEIIa B genome primers (3' end SNP, unique in/del overlap, and introduced mismatch; Figure 4) and used in conjunction with touchdown PCR to increase specificity.

Non-denaturing Polyacrylamide Detection Method

We report here the use of a modified screening technique that can be used to detect mutations in TILLING populations. We found equivalent mutations using the LI-COR and 3% non-denaturing polyacrylamide, suggesting that this system represents a viable low-cost alternative to the current technologies. Our 3% non-denaturing polyacrylamide system is based on ethidium bromide staining, eliminating the need for the genotyping instrument and fluorescently labeled primers. This is especially relevant in polyploid genomes as the fluorescent label attached to primers can reduce their genome specificity, requiring additional PCR optimization. Samples can be loaded directly after stopping the CJE digestion reaction with 0.225 M EDTA, eliminating the subsequent steps of sample purification and volume reduction required in LI-COR screens. We also found that by extending and optimizing the CJE digestion time (determined empirically, results not shown) we observed digestion of both strands, despite CelI being a single strand mismatch-specific endonuclease. This additional CJE activity eliminates the need of denaturing polyacrylamide gels which are more time consuming and technically more difficult than a non-denaturing system. Although we are able to increase the size of our target regions to over 1.5-kb, we were unable to find mutations in the first 10% of the sequences adjacent to each primer. Therefore the total sequence screened is roughly similar to the LI-COR method (1.3 to 1.4-kb) with the disadvantage that more sequence information is needed in our method to accommodate the larger distance between the primers and the region where mutations are effectively detected.

An additional advantage of the LI-COR system, is that the use of different dyes for each primer allows a precise estimation for the location of the mutant. In our set-up, two possible locations are estimated since we have no information as to whether the estimated distance is from the forward or the reverse primer. This implies an additional cost for sequencing a larger number of mutants to identify EMS-induced polymorphism in the desired regions. Another possible drawback of this set-up is the need for manual analysis of the gels since no software has been developed for this system. Despite this, gel image analysis requires approximately 20–30 minutes, similar to the time required with GelBuddy or similar gel analysis programs.

For both libraries we found a number of duplicate mutations, as well as a few triple and quadruple mutations, that were higher than expected by chance. These mutations are likely residual polymorphisms in the mutagenized seed, originated from residual heterozygous alleles in some of the plants used for the production of breeder's seed of Kronos or the seed stock for UC1041+Gpc-B1/Yr36. For example, we confirmed that the original seed stock of UC1041+Gpc-B1/Yr36 is polymorphic for a known 1-bp deletion in the coding region of the VRN-D3 allele [43]. These observations also suggest that the polyacrylamide detection method should be amenable for EcoTILLING [44].

The polyacrylamide detection method is especially relevant for species, such as wheat, that have no central TILLING service available. Even if a central service for wheat becomes available, individual researchers may need to TILL genotypes carrying specific alleles (such as for disease resistance) that may be absent in available TILLING populations. Although several alternative detection methods have been published, most rely on expensive equipment (sequencers, HPLC, gene analyzers) that precludes many individual laboratories from performing TILLING. The development of a non-denaturing polyacrylamide detection system makes TILLING more accessible to a larger set of researchers and breeding programs and may facilitate the development of multiple wheat TILLING populations. We plan to make the DNAs of our TILLING lines available on a cost recovery basis for other research groups to screen [45]. This should enable different research groups to screen for mutations in their gene of interest and expand their capabilities for wheat functional genomics.

We have pursued further characterization of over 20 mutants for the WKS [32] and SBEII genes. Although the phenotypic characterization of the SBEII mutants is beyond the scope of this work, we have successfully confirmed each WKS and SBEII mutation in its corresponding M₃ seed stock. This is an important final validation step, suggesting that the two libraries will be an effective reverse genetics resource.

Other groups are also currently developing TILLING resources in wheat. These include tetraploid libraries in variety 'Cham 1'[46] as well as hexaploid libraries in several varieties ('Alpowa', 'Louise' and 'Jagger' in the USA (Camille M. Steber, personal communication), Cadenza [46] and 'QAL2000' and 'Ventura' in Australia [21]). The generation of multiple TILLING libraries in wheat will allow an even greater flexibility and robustness to wheat TILLING as any targets which are missing in one population could be screened for in others.

Mutagenesis and population development and growth conditions

For the EMS treatments, 10–20 g of seeds were placed in 1 L plastic buckets with 40 mL 10% Tween 20 solution. The buckets were placed on an orbital shaker for 15 minutes at 112.5 rpm after which the solution was discarded. 250 mL of tap water were added to each bucket to wash off excess Tween 20 and they were then placed on the orbital shaker at 112.5 rpm for 5 minutes. The water was discarded and this washing step was repeated three more time. A final volume of 250 mL of tap water was used for the EMS treatment and the buckets were placed on the orbital shaker at 112.5 rpm for 18 hours at room temperature. After EMS treatment, seeds were thoroughly washed with tap water for 3 hours and then placed at 4°C for 5 days before being transferred to room temperature.

The EMS-treated seeds (M₁) were sown in individual cones with soil and the M₂ seeds from each individual plant were collected and labeled. Designators 'T6' and 'T4' were used for the hexaploid and tetraploid lines, respectively. Five M₂ seeds per line were sown in one pot and they were later thinned to leave a single M₂ plants per line.

DNA isolation and preparation of plates for pooling

Genomic DNA was isolated from 1,386 and 1,536 M₂ plants for the tetraploid and hexaploid populations, respectively, using a previously published protocol [47]. This large-scale DNA extraction protocol yields approximately 1–2 mg of DNA and was selected to generate abundant DNA of excellent quality for distribution. DNA concentrations were measured using a spectrophotometer (Nanodrop, Thermo Scientific, Wilmington, DE, USA) and standardized. Using equivalent amounts of DNA from individual plants, samples were pooled fourfold and organized into 96-well format. A total of four 96-well plates were built for each TILLING population, representing the 1,386 and 1,536 M₂ plants.

Primer design

Screening technique and two-step strategy

The protocol for the LI-COR detection method has been previously published [48]. The protocol for the targets visualized through the polyacrylamide detection methods includes a two-step screening approach. The first PCR screen of the complete set of DNA pools was carried out in a 25-μl reaction volume using 50–100 ng of pooled DNA (a large amount of DNA is required given the large genome size of wheat), 1 U of Taq polymerase and the following cycling conditions: initial denaturation at 95°C for 2 min, followed by 14 cycles of touchdown at 94°C for 20 s, from 64 to 57 for 30 s (0.5°C decrease per cycle) and extension at 72°C for 75 s. This touchdown cycle is followed by 37 cycles of 94°C for 20 s, 57°C for 30 s and 72°C for 75 s. A denaturing and re-annealing step is included at the end of the PCR reaction (99°C for 10 min, 70 cycles of 70°C for 20 s decreasing 0.3°C per cycle) to allow the formation of heteroduplexes if a mutation is present in the pool.

After PCR amplification, 12 μl of sample (~500 ng) was digested with celery juice extract (CJE) which was obtained using the protocol described by Till et al. [48]. Due to the inherent variability of different celery juice extracts, the optimal amount of CJE for heteroduplex-digestion was determined empirically [48] using targets with known mutations. This was done by performing TILLING reactions at 45°C for 30 minutes, but with varying amounts of CJE. Using too little CJE results in a strong top band (corresponding to the PCR reaction), but no visible lower bands (corresponding to the cut fragments). Too much CJE leads to complete digestion of the original PCR product resulting in a smear along the gel. A good signal-to-noise ratio is depicted in Figure 1A and 1C. Digestions for this study included: 12 μl PCR product, 0.12 μl CJE, 1.7 μl 10× digestion buffer [48] and 3.18 μl dH₂O for a final volume of 17 μl. The digestion was carried out at 45°C for 30 minutes and stopped immediately by adding 5 μl of 0.225 M EDTA per sample and mixing thoroughly. Two micro liters of bromophenol blue loading dye were added and the complete volume (24 μl) was loaded on the gel.

Samples were visualized on a 3% polyacrylamide gel (19:1 Acrylamide:bis ratio) in 0.5× TBE running buffer with ethidium bromide. We use a 100-lane vertical electrophoresis system (gel size: 22 cm tall by 1.5 mm thick; CBS Scientific, Del Mar, CA, USA) which was run at 350 V for 45–60 minutes. Gel images were analyzed manually on PowerPoint (Microsoft Corp., Seattle, WA, USA) and positive pools were identified by the presence of cleaved products whose combined size was similar to the original PCR product.

The second screen was performed using individuals DNAs only from the DNA pools that showed cleaved products in the first screen. PCR amplifications were as described before with the exception that the heteroduplex formation step was postponed. First, 6 μl from each PCR sample were pooled in four pairs (a+b, c+d, a+c and b+d) following the diagram in Figure 1B. The heteroduplex formation step was then performed on the mixed pairs and the samples followed the same detection protocol as described above. The pooling step is necessary to detect homozygous mutations in the M₂ plants since combining two samples allows heteroduplex formation and detection which would otherwise go undetected in a single homozygous sample. In addition to the identification of the individual from the pool that carries the mutation, this step also provides an independent validation of the mutations and a better estimation of its location within the target region.

Sequence analysis

To characterize the individual mutations, the residual PCR product from the selected individuals was used as a template for sequencing. PCR products (2 μl/sample) were cleaned using ExoSAP-IT (USB Corp., Cleveland, OH, USA) according to manufacturer's instructions, and subsequently sequenced using BigDye Terminator sequencing kit and an ABI -3730 DNA Sequencer (Applied Biosystems, Foster City, CA, USA).