sunTILL: a TILLING resource for gene function analysis in sunflower
© Sabetta et al; licensee BioMed Central Ltd. 2011
Received: 9 March 2011
Accepted: 30 June 2011
Published: 30 June 2011
Cultivated sunflower (Helianthus annus L.) is a globally important oilseed crop, subjected to intensive genetic and genomic studies. Although classical mutagenesis has successfully been applied to Helianthus genus in the past, we have developed the first sunflower TILLING resource.
To balance the maximum mutation density with an acceptable plant survival rate, a 'kill curve' analysis was first conducted with different ethylmethanesulfonate (EMS) dosages and different exposure times. According to the germination rate, a treatment with 0.7% EMS for 6 h was chosen. An M2 progeny of 3,651 fertile plants was obtained. Totally, 4.79% of the whole population showed clear aberrant phenotypes. A microsatellite analysis on a representative sample of the original seed stock and mutant lines confirmed the uniformity of the genetic background of plant material. The TILLING procedure was successfully applied to sunflower genome, initially by a Cel I-nuclease mismatch cleavage assay coupled with a DNA-pooling level test. To investigate the efficiency of the mutagenic treatment, a pilot screening was carried out on 1,152 M2 lines focusing on four genes, three involved in the fatty acid biosynthetic pathway and one for downy mildew resistance. A total of 9 mutant lines were identified and confirmed by sequencing; thereby, the estimated overall mutation frequency for the pilot assay resulted to be 1/475 kb.
A first TILLING population for a high throughput identification of EMS-induced point mutations in sunflower genome has been successfully obtained. This represents a powerful tool to a better understanding of gene function in sunflower.
Cultivated sunflower (Helianthus annuus L.) is a globally important oilseed crop subjected to intensive molecular genetic and genomic studies during recent decades . Sunflower belongs to the genus Helianthus, which is native to temperate areas of North America and includes 12 annual and 37 perennial species [2, 3]. The world production of sunflower is estimated to be 23.4 million tons ; seed production in the world increased by 24% (or 5 million metric tons) between 1993 and 2003 . After maize, it is the second largest hybrid crop, and the fifth largest among the oilseed crops, after soybean, rapeseed, cottonseed, and groundnut. The largest producers of sunflower seeds are Russia and Ukraine, followed by the European Union with 4.9 million tons , where this crop is mainly cultivated in France, Hungary, Italy, and Spain, thanks to their ideal climate. Sunflower can be used for different purposes: as an oilseed crop, edible confection, birdseed and, to a much lesser extent, as an ornamental for home gardens and the cut-flower industry.
Wild species of sunflower are characterized by a high genetic diversity as a consequence of their adaptation to a wide range of environments . They harbour a significant variability with respect to a number of traits such as disease and pest resistance, quality and composition of seed compounds. The genetic polymorphism of cultivated sunflower is reduced to 40-50% of the diversity observed in the wild germplasm as a result of domestication . The development of new genetic variability in sunflower is aimed at producing a source of agriculturally useful alleles or new genotypes.
Many efforts in sunflower plant breeding have been proposed for the improvement of desirable traits. Genetic engineering , traditional breeding approaches , in vitro breeding techniques (somaclonal variation) and conventional mutation technologies  have been used to improve yield, oil quality and disease-, salt- and pest-resistance of the sunflower crop. TILLING (Targeting Induced Local Lesions IN Genomes) is a reverse genetic technique that is suitable for most plants . In a TILLING project, the chemical mutagenic treatment providing an easy and cost-effective way to saturate a genome is coupled with a PCR-based mutation detection. By using alkylating agents such as ethylmethanesulfonate (EMS) that cause random point mutations at high density, an allelic series of missense mutations can be discovered by TILLING; short insertion/deletions (INDELs) are reported to be detected by this technology too . Thus, with only a small population, multiple alleles of a specific gene may be obtained regardless of the gene size . Gene regions are targeted for mutation discovery, using PCR and standard SNP (single nucleotide polymorphism) discovery methods. The success of the TILLING approach relies on the construction of high quality DNA mutant libraries, in which DNA sampling and preparation are the most critical steps. An accurate evaluation of the genetic uniformity both of seed stock used for the mutagenic treatment and of the mutagenized material is extremely important. As reported by Wu et al. , genotyping plant material by means of microsatellite markers (SSR) can be an efficient tool to verify the uniformity and purity of the seed stock used in a TILLING project.
Although classical mutagenesis has been successfully employed in sunflower breeding programs over the last 40 years, sunflower TILLING resources have not been reported yet. The objectives of the present study were: 1) the development of the first sunflower TILLING population, where single nucleotide lesions are induced by conventional EMS mutagenesis; 2) the optimization of TILLING procedure for an efficient sunflower SNP detection system and the mutation density evaluation by screening a subset of M2 lines; 3) the phenotyping of the M2 generation for the presence of mutations in comparison to the untreated sunflower.
Plant material and EMS mutagenesis
Seeds from the inbred line GV342 of cultivated sunflower (Helianthus annus) were used for the mutagenesis experiment. To test the 'kill curve', eight batches of 100 seeds each were treated with different EMS concentrations (0.5, 0.7, 1.0 and 1.5%); all concentrations were also tested for two different exposure times (3 and 6 h) at 20°C and applying a gentle shaking. Then seeds were washed with tap water three times for 5 min each and a fourth time for 30 min. Subsequently, they were transferred on wet Whatman paper in Petri dishes in a growth chamber at 20°C and 8 h photoperiod.
Based on the percentage germination, the 6 h-treatment with 0.7% EMS was chosen to mutagenize a batch of ~30,000 seeds. After an overnight air-drying on Whatman paper at room temperature, seeds were mechanically sown in the experimental field close to Bari (Southern Italy) with a cultivation distance between and within rows of 70 and 30 cm respectively. Fertilization and herbicidal treatments were applied once and the field was regularly watered. To allow self-pollination, the floral bottom of each M1 plant was covered with a TNT bag before anthesis. At physiological maturity, the main head from each individual plant was harvested and M2 seeds were collected and stored at 4°C.
Four to ten M2 seeds for each family were sowed in the field using the same conditions as for the M1 generation. At the 3-4 leaf stage, most of the emerged plants were manually thinned to have one plant per family. M2 plants were regularly monitored for the presence of morphological mutants, using the untreated plants as reference, and forced to self-pollinate by bagging the inflorescences. Germination and sterility rates of the M2 generation were also evaluated. Based on a visual characterization of plants at the key development stages, from germination to maturity, a systematic phenotypic scoring of the mutant population was carried out. Seven main traits were chosen from the worldwide sunflower descriptor list promoted by the International Board for Plant Genetic Resources. Phenotypes were organized into the following categories: leaf colour, shape and size; internode length; type of branching; venation colour and seed colour. Phenotypic alterations were annotated under field-growing conditions and referred to untreated plants; a detailed photographic documentation was also collected. M3 seeds from individual plants were harvested, aired and vacuum-stored at 4°C.
DNA extraction and normalization
Leaf material from each M2 plant was collected for DNA isolation. Approximately 100 mg of lyophilised mature leaves were located in 2 ml single tubes and crashed using a Mixer-Mill (Retsch-Muhle MM30, Qiagen, Germany) two times for 45 sec at 30 Hz with two 3.175/III-mm-diameter inox spheres. DNA was extracted according to Li et al.  with one modification: all the centrifugations in the second part of the protocol were elongated up to 2-5 min. Genomic DNA concentration and quality were measured both by means of a NanoDrop™ 1000 Spectrophotometer (Thermo Scientific, USA) and 1% agarose gel. DNAs were transferred into 96-well plates and normalized to a standard concentration of 20 ng/μl by adding 0.1 × TE buffer (10 mM Tris-HCl pH 8.0 and 1 mM EDTA).
Genetic uniformity analysis by microsatellite markers.
Detected alleles (bp)
M 2 lines
245 - 255
320 - 340
126 - 130 - 140
365 - 370 - 380
250 - 260 - 270
380 - 395 - 410
295 - 320
190 - 200
311 - 313 - 315 - 320
405 - 425 - 460
200 - 220 - 235
368 - 370
240 - 250
440 - 450 - 460
290 - 310
SNP detection test
Two sunflower genotypes (i.e. XRQ and PSC8 from INRA, France), known to be polymorphic for a single nucleotide at one specific locus (Dr. P. Vicourt, personal communication), were used as positive control to optimize the sunflower SNP detection system. Two different tests were carried out: a Cel I-nuclease mismatch cleavage assay and a DNA pooling test. In the first test, DNAs from XRQ and PSC8 genotypes were mixed together in the same proportions to create two different 2-fold pools, that were used as template for the amplification and the subsequent preliminary digestion test. Once optimized the Cel I restriction analysis, the best enzymatic conditions were used for the DNA pooling test, which was performed by the combination of DNAs from XRQ and PSC8 genotypes in different proportions simulating 2-, 4-, 8-, 12-, 16-, 20-fold pools. Forward and reverse primers were 5'-end labelled with IRDy671 and IRDy781 dyes respectively (Eurofins MWG Operon, Germany). Mixtures of 3:2 (labelled:unlabelled) ratio for the forward primers and 3:1 (labelled:unlabelled) ratio for the reverse primers were used in the PCR reaction. PCR amplification was carried out in a volume of 20 μl containing 20 ng/μl of pooled DNA, 1 × Phusion HF Buffer (Finnzymes, Finland), 0.2 mM each dNTP, 3% DMSO, 0.3 μM each primer and 0.015 U Phusion High-Fidelity DNA Polymerase (Finnzymes, Finland). The PCR was performed in a thermal cycler (Applied Biosystems 9800 Fast Thermal Cycler) using the following cycling program: 98°C for 30 sec; 35 cycles at 98°C for 10 sec, 52°C for 30 sec, 72°C for 30 sec; then a final extension at 72°C for 10 min. The amplification step was followed by the heteroduplex formation: inactivation at 99°C for 10 min; 23 cycles of re-annealing process for 20 sec at 70°C to 69.4°C, decrementing 0.9°C per cycling .
Amplification of candidate genes
The sequences of kasII, kasIII, fad2-1 and AY490791 genes were identified on SRS-EMBL, NCBI, CGP, TIGR and HeliGene databases and subjected to bioinformatic analysis by CODDLE (http://www.proweb.org/input/) and SoftBerry (http://linux1.softberry.com/berry.phtml) software to predict the gene structures. All the primer pairs were set by using Primer3 (http://frodo.wi.mit.edu/primer3/) and OligoExplorer v. 1.2 programs. PCR reactions were performed in a thermal cycler (Applied Biosystems 9800 Fast Thermal Cycler) using 96-well microtiter plates and carried out in a 20 μl volume consisting of dH2O, 1 × Phusion HF Buffer, 0.2 mM each dNTP, 3% DMSO, 0.3 μM each primer (with a ratio of 3:2 and 3:1 labelled:unlabeled for forward and reverse primer, respectively) (Eurofins MWG Operon, Germany), 0.015 U Phusion High-Fidelity DNA Polymerase (Finnzymes, Finland). The thermocycling conditions were 98°C for 30 sec for initial denaturation, followed by 30 cycles at 98°C for 10 sec, an annealing temperature specific for each primer pair for 30 sec, 72°C for 40 sec, one cycle at 72°C for 10 min and 4°C hold for storage. Five μl of each amplified product were electrophoresed on 1% SeaKem® LE Agarose gel (Lonza, Switzerland) to verify the PCR efficiency and concentration before the subsequent TILLING step.
Cel I nuclease mismatch cleavage assay and purification step
The preliminary Cel I digestion test was carried out using progressive enzyme dilutions and different reaction times. Ten μl of amplification product were incubated for different times (15, 30 and 45 min) at 45°C with differentially diluted (1:2, 1:4, 1:10 and 1:20) Surveyor Cel I enzyme (Transgenomics, Omaha, USA) in a 10 × Buffer (10 mM HEPES, pH 7.0, 10 mM KCl, 10 mM MgCl4 7H20, 0.002% Triton X-100 and 10 μg/ml bovine serum albumen). Cel I digestion was stopped by adding 5 μl of 75 mM EDTA (pH 8.0) and freezing samples (-20°C) for at least 3 h. For the pilot screening, a 1:20 enzyme dilution for 45 min reaction was always used as standard conditions.
Sample purification was performed by precipitation with 5 μl 3 M Na-Acetate (pH 5.2), 75 μl 99.8% EtOH and 10 μl dH2O. After shaking for 15-20 min, the samples were centrifuged at 4,500 rpm for 30 min at 20°C; the supernatant was carefully discarded. Then plates were spin upside down on a filter paper up to 500 rpm, washed with 100 μl 75% EtOH and re-centrifuged at 4,500 rpm for 30 min at 20°C. Again, the supernatant was discarded on filter paper and plates were spin facedown on a filter paper up to 500 rpm. Then, samples were dried for 20 min at room temperature and resuspended in 8 μl formamide loading buffer (33% deionised formamide, 25 mM Tris pH 7.5, 25 mM EDTA and ~0.02% bromophenol blue) by shaking (300 rpm) for 5-10 min.
Li-COR gel electrophoresis
Samples were denatured by incubation at 95°C for 5 min, placed on ice, loaded on gel by a 100 teeth paper comb and finally electrophoresed on Li-COR DNA Analyser. Electrophoresis was performed through a 6.0% LongRanger® polyacrylamide (Biozym, FMC Corporation, Austria), 7 M urea gel in 1 × Tris-Borate-EDTA running buffer. The laser focusing and gel warming to 50°C were obtained by a pre-run at 1.500 V, 40 mA and 40 W for 20 min. Gel images were analysed visually for the presence of cleavage products using GelBuddy (http://www.proweb.org/gelbuddy) and Adobe Photoshop (Adobe System Inc., CA, USA) software.
To estimate the mutation frequency for each screened gene, the total number of analysed base pairs (obtained multiplying the number of screened plants per the fragment size) was divided by the number of identified mutations. Because of priming and PCR artefacts on the top and the bottom of LiCOR gels, the detection of mutations in the terminal 50 bp at each end of the amplicons was difficult and for this reason 100 bp were subtracted from the size of each amplicon to obtain the effective screened window size.
Production of a sunflower mutant population
In the past, a broad spectrum of genetic resources was developed for cultivated sunflower and most of the available germplasm has been accurately characterized by molecular markers . Information about genetic diversity [19–21], molecular mapping [15, 22], gene mapping [23–26], construction of cDNA, BAC and BIBAC libraries [27, 28] are currently available. Many examples of successfully applied traditional mutagenesis have been reported [8, 29, 30], but no sunflower TILLING resource has previously been developed.
Kill curve analysis.
Exposure time (h)
EMS concentration (%)
Thirty thousand M0 seeds were mutagenized with 0.7% EMS for 6 h and grown in field conditions. About 13,000 M1 plants were obtained, but only 50% of them reached the complete maturity and set seeds. Four to 10 M2 seeds per each family were sowed in the field and, at the stage of 3-4 leaves, most of the emerged plants were manually thinned to leave only one M2 plant per family. About 64% of the sown seeds were able to germinate: thereby, 4,211 M2 plants were obtained and regularly monitored for the presence of phenotypic variations in comparison with the untreated wild-type. Finally, since M3 seeds were harvested from 86.7% of fertile plants, an M2 population of 3,651 lines was used for leaf DNA sampling and M3 seed stocking.
Phenotypic classes observed by the screening of M2 generation.
1. Leaf colour
2. Leaf size
3. Leaf shape
4. Internode length
5. Type of branching
Fully branched with central head
Fully branched without central head
6. Venation colour
7. Seed colour
Plants with shorter (Figure 1e) or longer internodes (Figure 1f) than the wild-type in total represented 4% of the population. Regarding the 'plant habit', 98.8% of the population were like the wild-type in not showing lateral shoot branching; only one plant showed top branching, while 18 plants had basal branching (Figure 1g). Ten fully branched plants with the main head (Figure 1h) were observed, while in 15 fully branched plants the main head was missing.
One mutant with an anthocyanic venation colour (Figure 1i) was observed in comparison with the common green colour of the untreated plants. White or anthocyanic coloured seeds, in comparison with the dark wild-type, were grouped in the 'seed colour' category. Eight white and 31 anthocyanic seeds were observed. As a general indication of M3 seed yield, only 13.3% of the mutant population resulted sterile and seeds were collected from 3,651 M2 plants.
In conclusion, the most common observed phenotypic variation was for plant size (146 plants), followed by leaf colour (68 plants); 42 and 21 plants showed variations in leaf shape and size, respectively. Four distinct branching typologies affected a total of 44 plants, while there were up to 39 individual lines with pale coloured seeds.
In total 175 M2 plants, representing 4.8% of the population, showed at least one altered trait. Among these, 138 plants that displayed multiple mutant traits (meaning that they fell into more than one major category) were counted only once. This led to an underestimation of the mutation frequency at the phenotypic level. Moreover, some phenotypic traits (i.e. root or flower characteristics) were not taken into account, therefore the mutation frequency observed in the M2 generation may have been underestimated. Detailed morphological and biochemical characterization of larger numbers of M2 plants would likely have resulted in the scoring of increased numbers of mutant phenotypes.
Genetic uniformity analysis
Optimization of sunflower TILLING procedure
To optimize the TILLING procedure on sunflower genome, two genotypes (i.e. XRQ and PSC8 from INRA, France), carrying a single nucleotide polymorphism at one specific locus in cluster HuCL00001C250 (http://www.heliagene.org), were used as positive control. These SNP markers allowed the validation of each step of the TILLING procedure.
Two alternative DNA extraction protocols were initially tested. The use of a commercial kit (Nucleospin®, Mecherey-Nagel, Germany) provided high quality DNA, but low yields made this protocol unusable for the development of a TILLING platform. Some modifications to the DNA extraction protocol from Li et al.  as reported in Methods, allowed us to obtain high quality and highly concentrated DNA from both wild-type and mutant samples.
Since the efficiency of heteroduplex detection could be influenced by the pool size in terms of DNA amount [10, 33], a DNA pooling assay was performed. DNAs from XRQ and PSC8 genotypes were mixed together in different proportions, simulating 2-, 4-, 8-, 12-, 16- and 20-fold pools. As expected, a progressive reduction of the SNP detection was observed as the pool level increased (Figure 3B): even though the cleaved products were still visible up to the pool ratio of 1:19, the 8-fold pooling size in a two-dimensional format was chosen for the M2 DNAs arrangement.
Candidate genes and reverse genetic pilot screening
To estimate the efficiency of the EMS treatment and to evaluate the mutation frequency, a pilot assay was initially performed on DNAs from 1,152 M2 lines, arranged in a two dimensional 8-fold pooling scheme. In particular, four candidate genes were selected: the kasII, kasIII, fad2-1 and AY490791 genes.
Pilot assay on a sub-set of M2 lines and estimation of the mutation frequencies.
Fragment size (bp)
Number of screened plants
Another important enzyme involved into the fatty acid pathway and responsible for the synthesis of linoleic acid from oleic acid is the isoform 1 of the Δ-12 oleate desaturase (FAD2; EC 126.96.36.199) encoded by the fad2-1 gene . It is a membrane-bound enzyme that catalyzes the first extra-plastidial desaturation in plants, converting oleic acid to linoleic acid. The reaction involves the concomitant reduction of molecular oxygen to water and requires the presence of an electron donor system . To identify the most suitable region for TILLING screening a 3,008 bp genomic sequence (EMBL: AY802989) was subjected to bioinformatic analysis. Based on the CODDLE prediction, a primer pair amplifying a 960 bp fragment was designed that flanked the unique large exon of the gene and used for screening. Three mutations were identified: one resulted in a missense change (F26L), a second caused a silent change (R46 =) and a third was situated in the non-coding region. The mutation frequency for this gene was estimated as 1/332 kb (Table 4).
The AY490791 gene is annotated on NCBI website as a LRR gene involved into Plasmopara halstedii resistance process. This biotrophic oomycete is the causal agent of downy mildew in wild and cultivated sunflower and over 100 host species in the Asteraceae family . This is one of the main diseases in the cultivation of sunflower, causing economic losses of up to 95% . The most promising region of this gene, a highly conserved 1,430 bp sequence, was selected. The molecular screening of 1,152 lines allowed us to identify two individuals with mutations in their intronic sequences, leading to an estimated mutation frequency for this gene of 1/766 kb.
Conventional mutagenesis using chemical and physical mutagens have been widely used as forward genetic approach for crop improvement, including sunflower species [29, 41, 42]. Mutagenic treatments have been successfully applied to improve oil-quality and to induce tolerance to biotic and abiotic stress, such as salt-, disease-, pest-, cold- and metal concentration-resistance . New functional approaches, like RNA interference  or transposon mutagenesis , have been developed in the last decade, as a consequence of the rapid accumulation of large-scale sequence data and the increasing need to determine gene function for crop improvement. Several TILLING populations have been produced for different species, using EMS, MNU, diepoxybutane, NaN3 or a mixture of these chemicals as mutagenic agents [11, 45]. Even though the sunflower genome has not been sequenced yet, a large number of EST sequences deriving from the Composite Genome Project and from several databases are available (http://compgenomics.ucdavis.edu/cgp_sitemap.php); (http://compbio.dfci.harvard.edu/tgi/plant.html); (http://www.heliagene.org). In the present study the first sunflower population for reverse genetic has been generated. To scan for mutations, fluorescent-labeling and double-stranded cutting of mismatches with the Cel I cleavage assay were combined with a gel electrophoresis-based analysis, according to Colbert et al. .
The identification of the optimal dose of a chemical mutagen that maximizes the mutation frequency without extreme toxic effects on the vitality and the fertility of biological tissues is a key factor during the establishment of a TILLING project. In fact, the main problem associated with chemical mutagenesis is the high cell toxicity of the mutagen, which can decrease the germination rate and viability of M1 seedlings to unacceptable levels without providing sufficient mutation density in the survival plants [11, 46]. Moreover, in the next generation, undesirable phenotypes could appear as a result of the random distribution of mutations in the genome. This usually translates into the necessity to perform a preliminary 'kill curve' analysis, to find a compromise between mutagen toxicity, genome mutation saturation and possible accumulation of undesirable phenotypes [47, 48].
Our aim was to assemble a sunflower TILLING population of fertile individuals, so that the progeny of each plant carrying a mutant allele could be directly recovered. To achieve this we initially performed a 'kill curve' analysis and then the most suitable EMS dosage was applied to our seed material without excessively compromising the M1 germination rate (~43%) and the M2 fertility rate (87%). Our data cannot be directly compared with those reported in other TILLING projects, since different negative effects to mutagenic treatment could be displayed according to the genetic background of the treated species or variety [17, 45, 49, 50]. A sunflower TILLING population of 3,651 independent mutant lines was generated in this work, providing an important tool for the identification of interesting phenotypes.
Since considerable efforts are needed for the propagation of mutant lines to the second generation, a phenotypic marker of mutation frequency in the early stages of population development is often desirable. Thereby, seven phenotypic descriptors among those listed and promoted by the International Board for Plant Genetic Resources were chosen for the phenotypic screening of our population. Moreover, because of the out-crossing nature of sunflower and the need to bag the inflorescences for self-pollination, the phenotypic scoring for flowering or flower morphology could not be performed. However, about 5% of the population showed altered morphologies in comparison with the wild-type. Some of the mutant lines (37) showed a single altered trait, while 138 lines displayed multiple mutant traits. Moreover, future screenings of the M3 generation will clarify the nature of the observed sunflower mutant lines by comparison of the phenotypic percentages recorded in two generations, as reported for indica rice IR64  and tomato , where not all the identified M2 mutant lines were confirmed in the next generation.
Information about the genetic uniformity of both the untreated and mutagenized plant material and about the maintenance of seed stocks have been provided. This molecular marker analysis allowed us to detect the existence of any contamination of the population. The original seed stock represented an inbred line (F8), and two more selfing generations were performed to obtain the M2 progeny; thereby a high level of genetic uniformity and homozygosity was expected. However, out-crossing events with unwanted pollen may occur in field conditions, when an open pollinated species is analyzed and obliged to self-fertilize by bagging its inflorescences. The SSR genotyping performed in the present work confirmed the genetic uniformity and purity of both original seed stock and randomly analyzed M2 plants of the TILLING population. Wu et al.  reported this approach to be a cost-effective method of purity control of mutant stocks; they routinely genotyped the identified rice mutants to ensure they originated from the original seed stock, before seed distribution . Moreover, SSR analysis allowed Xin et al.  to identify that 30% of M2 sorghum plants resulted from cross-pollination.
Several enzymes have been used for mismatching specific cleavage, but Cel I is the most common one in TILLING projects [11, 45] for genotyping applications  and detection of heterozygous polymorphism . The Cel I-mismatch cleavage assay allowed us to establish the best experimental conditions by the background reduction from gel images and the increase in SNP resolution power. The results obtained in our test provided a further advantage, since the use of minimal enzyme concentrations translated into a considerable economic saving.
Since DNA sample quality, normalization and pooling level can directly affect the efficiency and cost of mutation discovery [9, 11], the main goal of our DNA pooling test was to maximize the throughput by increasing DNA pooling level while still clearly detecting the expected polymorphism. Although the simulation of a 20-fold pool of individual samples resulted in a successful detection of mutations, the 8-fold pool and the two-dimensional format were chosen as the pooling strategy for our sunflower TILLING population.
The TILLING strategy we developed has been successfully applied here as a pilot assay on 1,152 sunflower M2 lines. Four genes have been subjected to the reverse genetic screening: the kasII and kasIII genes, respectively codifying the isoforms II and III of the β-keto-acyl-ACP-synthetase; the fad2-1 gene, encoding the enzyme responsible of the converting reaction of oleic acid to linoleic acid; the AY490791 gene, involved in P. halstedii resistance. Our interest was first focused on some key enzymes of the fatty acid pathway, because of the interest in increasing the nutritional value of sunflower oil by the reduction of the ratio of saturated to unsaturated fatty acids. Moreover, P. halstedii is one of the most dangerous pathogens that affects sunflower cultivation in the Mediterranean area . Therefore the availability of a stable and effective system, as genetic resistance, for the pest-control results of prime importance.
Since few genomic sequences are publicly available for sunflower, the reverse genetic screening was preceded by an accurate reconstruction of candidate gene models, by the amplification and the subsequent sequencing of short overlapping fragments. We focused our efforts in this first step on identifying for each candidate gene the most promising region for TILLING analysis. In this way, new primer pairs flanking this region could be targeted to the intronic sequences, with the aim of improving the screening efficiency on the coding regions in the pilot assay. In total, nine mutant lines have been identified. Each has been confirmed by sequencing and genotyped by microsatellite markers to exclude any individuals originating from cross-pollination events. The results of this first reverse genetic screening translated into an average mutation frequency of 1/475 kb, which is not substantially different from those observed in TILLING populations of other diploid species such as barley (1/374 kb for cv Morex and 1/500 kb for cv Barke) [55, 56], Arabidopsis (1/300 kb) , Brassica oleracea (1/447 kb) , rice (1/280 kb for cv Nipponbare, 1/135 kb for cv Taichung 65 and 1/500 kb for cv IR64) [12, 58, 59] and sorghum (1/526 kb) .
The establishment of the EMS-mutagenized TILLING population described here represents an important advance in the generation of new genetic variation in sunflower. The development of TILLING technology can lead to the identification of new alleles that may be directly of value for crop improvement. Furthermore, intensive investigation of the role of key genes, becomes a feasible goal, especially important where genomic information is lacking, as it is in sunflower. Moreover, TILLING methodology makes it possible to focus on specific genes or genomic regions, bypassing problems with other functional genomic tools (such as T-DNA knock-outs or RNAi-based gene silencing) that require the generation of transgenic plants. Thus we have developed and established an exciting tool for forward and reverse genetics in sunflower, one that is available for scientific collaborations and that aims to contribute to a global understanding of sunflower gene organization and regulation.
The authors are grateful to Dr. M. Turi and Prof. GP. Vannozzi from the University of Udine (Italy) for kindly supplying the GV342 and hybrid seed materials; to Prof. F. Sunseri from the University of Reggio Calabria (Italy) for supporting the EMS treatment; to Prof. P. Vincour and Dr. D. Rengel from LIPM-INRA/CNRS of Toulouse (France) for kindly providing the sunflower XRQ and PSC8 genotypes used as SNP markers. Moreover, authors are grateful to Dr. Nils Stein and Dr. Sven Gottwald at IPK Institute of Gatersleben (Germany) for their technical support during the implementation of the sunTILL platform. This research was funded by the University of Bari, projects: Ateneo PRIN-COFIN 2006 and 2007, Ateneo 2008 and ALISAL-MIPAF 2009.
- Jan CC, Seiler GJ: Oilseed Crops. Genetic resources, chromosome engineering and crop improvement. Edited by: Singh RJ. 2007, New York: Taylor & Francis Group, 4: 103-168. Chapter 5Google Scholar
- Harter AV, Gardner KA, Falush D, Lentz DL, Bye RA, Rieseberg1 LH: Origin of extant domesticated sunflowers in eastern North America. Nature. 2004, 430: 201-205. 10.1038/nature02710.View ArticlePubMedGoogle Scholar
- Seiler GJ: Utilization of wild sunflower species for the improvement of cultivated sunflower. Field Crops Research. 1992, 30 (Suppl 3-4): 195-230.View ArticleGoogle Scholar
- U.S. Department of Agriculture. [http://www/fas.usda.gov/oilseeds/circular/2005]
- Kleingartner LW: World outlook and future development of sunflower markets around the world. Proceedings of the 16th International Sunflower Conference, Paris. Edited by: Seiler GJ. 2004, International Sunflower Association, 69-77.Google Scholar
- Micic Z, Hahn V, Bauer E, Schon C, Melchinger A: QTL mapping of resistance to Sclerotinia midstalk rot in RIL of sunflower population NDBLOSsel × CM625. Theor Appl Genet. 2005, 110 (Suppl 8): 1490-1498.View ArticlePubMedGoogle Scholar
- Lucas O, Kallerhoff J, Alibert G: Production of stable transgenic sunflowers (Helianthus annuus L.) from wounded immature embryos by particle bombardment and co-cultivation with Agrobacterium tumefaciens. Mol Breeding. 2000, 6: 479-487. 10.1023/A:1026583931327.View ArticleGoogle Scholar
- Encheva J, Shindrova P, Penchev E: Developing mutant sunflower lines (Helianthus annuus L.) through induced mutagenesis. HELIA. 2008, 31 (Suppl 48): 61-72.View ArticleGoogle Scholar
- McCallum CM, Comai L, Greene EA, Henikoff S: Targeting Induced Local Lesions IN Genomes (TILLING) for plant functional genomics. Plant Physiol. 2000, 123: 439-442. 10.1104/pp.123.2.439.PubMed CentralView ArticlePubMedGoogle Scholar
- Greene EA, Codomo CA, Taylor NE, Henikoff JG, Till BJ, Reynolds SH, Enns LC, Burtner C, Johnson JE, Odden AR, Comai L, Henikoff S: Spectrum of chemically induced mutations from a large-scale reverse-genetic screen in Arabidopsis. Genetics. 2003, 164: 731-740.PubMed CentralPubMedGoogle Scholar
- Till BJ, Comai L, Henikoff S: TILLING and EcoTILLING for Crop Improvement. Genomics-Assisted Crop Improvement. Volume: Genomics Approaches and Platforms, Chapter 17. Edited by: Varshney RK, Tuberosa R. 2007, Springer, Netherlands, 333-349.View ArticleGoogle Scholar
- Wu JL, Wu C, Lei C, Baraoidan M, Bordeos A, Madama MRS, Ramos-Pamplona M, Mauleon R, Portugal A, Ulat VJ, Bruskiewich R, Wang G, Leach J, Khush G, Leung H: Chemical- and irradiation-induced mutants of indica rice IR64 for forward and reverse genetics. Plant Mol Biol. 2005, 59: 85-97. 10.1007/s11103-004-5112-0.View ArticlePubMedGoogle Scholar
- International Board for Plant Genetic Resources (IBPGR): Descriptors for cultivated and wild sunflower. IBPGR Secretariat, Rome, Italy. 1985, 37-Google Scholar
- Li JT, Yang J, Chen DC, Zhang XL, Tang ZS: An optimized mini-preparation method to obtain high-quality genomic DNA from mature leaves of sunflower. Genet Mol Res. 2007, 6 (Suppl 4): 1064-1071.PubMedGoogle Scholar
- Tang S, Yu JK, Slabaugh MB, Shintani DK, Knapp SJ: Simple sequence repeat map of the sunflower genome. Theor Appl Genet. 2002, 105: 1124-1136. 10.1007/s00122-002-0989-y.View ArticlePubMedGoogle Scholar
- Paniego N, Echaide M, Munoz M, Fernandez L, Torales S, Faccio P, Fuxan I, Carrera M, Zandomeni R, Suarez EY, Hopp HE: Microsatellite isolation and characterization in sunflower (Helianthus annuus L.). Genome. 2002, 45: 34-43. 10.1139/g01-120.View ArticlePubMedGoogle Scholar
- Colbert T, Till BJ, Tompa R, Reynolds S, Steine MN, Yeung AT, McCallum CM, Comai L, Henikoff S: High-throughput screening for induced point mutations. Plant Physiol. 2001, 126: 480-484. 10.1104/pp.126.2.480.PubMed CentralView ArticlePubMedGoogle Scholar
- Zerr T, Henikoff S: Automated band mapping in electrophoretic gel images using background information. Nucleic Acids Res. 2005, 33: 2806-2812. 10.1093/nar/gki580.PubMed CentralView ArticlePubMedGoogle Scholar
- Chapman MA, Pashley CH, Wenzler J, Hvala J, Tang S, Knapp SJ, Burke JM: A genomic scan for selection reveals candidates for genes involved in the evolution of cultivated sunflower (Helianthus annuus). Plant Cell. 2008, 20 (Suppl 11): 2931-45.PubMed CentralView ArticlePubMedGoogle Scholar
- Iqbal MA, Sadaqat HA, Khan IA: Estimation of genetic diversity among sunflower genotypes through random amplified polymorphic DNA analysis. Genet Mol Res. 2008, 7 (Suppl 4): 1408-13.View ArticlePubMedGoogle Scholar
- Lai Z, Kane NC, Zou Y, Rieseberg LH: Natural Variation in Gene Expression Between Wild and Weedy Populations of Helianthus annuus. Genetics. 2008, 179: 1881-1890. 10.1534/genetics.108.091041.PubMed CentralView ArticlePubMedGoogle Scholar
- Lai Z, Livingstone K, Zou Y, Church SA, Knapp SJ, Andrews J, Rieseberg LH: Identification and mapping of SNPs from ESTs in sunflower. Theor Appl Genet. 2005, 111: 1532-44. 10.1007/s00122-005-0082-4.PubMed CentralView ArticlePubMedGoogle Scholar
- Pérez-Vich B, Berry ST, Velasco L, Fernández-Martínez JM, Gandhi S, Freeman C, Heesacker A, Knapp SJ, Leon AJ: Molecular mapping of nuclear malesterility genes in sunflower. Crop Sci. 2005, 54: 1851-7.View ArticleGoogle Scholar
- Rojas-Barros P, Hu J, Jan CC: Molecular mapping of an apical branching gene of cultivated sunflower (Helianthus annuus L.). Theor Appl Genet. 2008, 117: 19-28. 10.1007/s00122-008-0748-9.View ArticlePubMedGoogle Scholar
- Yue B, Cai X, Vick B, Hu J: Genetic characterization and molecular mapping of a chlorophyll deficiency gene in sunflower (Helianthus annuus). J Plant Physiol. 2009, 166: 644-651. 10.1016/j.jplph.2008.09.008.View ArticlePubMedGoogle Scholar
- Yue B, Vick BA, Yuan W, Hu J: Mapping one of the 2 genes controlling lemon ray flower color in sunflower (Helianthus annuus L.). J Hered. 2008, 99: 564-7. 10.1093/jhered/esn033.View ArticlePubMedGoogle Scholar
- Tamborindeguy C, Ben C, Liboz T, Gentzbittel L: Sequence evaluation of four specific cDNA libraries for developmental genomics of sunflower. Mol Gen Genomics. 2004, 271: 367-375. 10.1007/s00438-004-0989-5.View ArticleGoogle Scholar
- Feng J, Vick BA, Lee MK, Zhang HB, Jan CC: Construction of BAC and BIBAC libraries from sunflower and identification of linkage group-specific clones by overgo hybridization. Theor Appl Genet. 2006, 113: 23-32. 10.1007/s00122-006-0265-7.View ArticlePubMedGoogle Scholar
- Nehnevajova E, Herzig R, Federer G, Erismann KH, Schwitzguébel JP: Chemical mutagenesis - a promising technique to increase metal concentration and extraction in sunflowers. Int J Phytoremediation. 2007, 9 (Suppl 2): 149-65.View ArticlePubMedGoogle Scholar
- Skorić D, Jocić S, Sakac Z, Lecić N: Genetic possibilities for altering sunflower oil quality to obtain novel oils. Can J Physiol Pharmacol. 2008, 86 (Suppl 4): 215-21.View ArticlePubMedGoogle Scholar
- Hvarleva Tz, Bakalova A, Chepinski I, Hristova-Cherbadji M, Hristov M, Atanasov A: Characterization of Bulgarian sunflower cultivars and inbred lines with microsatellite markers. Biotechnol & Biotechnol Eq. 2007, 408-412.Google Scholar
- Solodenko A, Sivolap Y: Genotyping of Helianthus based on microsatellite sequences. HELIA. 2005, 28 (Suppl 42): 19-26.View ArticleGoogle Scholar
- Till BJ, Reynolds SH, Weil C, Springer N, Burtner C, Young K, Bowers E, Codomo CA, Enns LC, Odden AR, Greene EA, Comai L, Henikoff S: Discovery of induced point mutation in maize genes by TILLING. BMC Plant Biol. 2004, 4: 12-10.1186/1471-2229-4-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Carlsson1 AS, LaBrie1 ST, Kinney AJ, von Wettstein-Knowles P, Browse J: A KAS2 cDNA complements the phenotypes of the Arabidopsis fab1 mutant that differs in a single residue bordering the substrate binding pocket. Plant J. 2002, 29 (6): 761-770. 10.1046/j.1365-313X.2002.01253.x.View ArticlePubMedGoogle Scholar
- McGuire KA, McGuire JN, von Wettstein-Knowles P: Acyl carrier protein (ACP) inhibition and other differences between β-ketoacyl-synthase (KAS) I and II. Biochem Soc Transactions. 2000, 28: part 6Google Scholar
- Abbadi A, Brummel M, Spener F: Knockout of the regulatory site of 3-ketoacyl-ACP synthase III enhances short- and medium-chain acyl ACP synthesis. Plant J. 2000, 24 (Suppl 1): 1-9.View ArticlePubMedGoogle Scholar
- Shanklin J, Cahoon EB: Desaturation and related modifications of fatty acids. Annu Rev Plant Physiol Plant Mol Biol. 1998, 49: 611-641. 10.1146/annurev.arplant.49.1.611.View ArticlePubMedGoogle Scholar
- Smith MA, Cross AR, Jones OTG, Griffiths WT, Stymne S, Stobart K: Electron-transport components of the 1-acyl-2-oleoyl-sn glycerol-3-phosphocholine D12-desaturase (D12-desaturase) in microsomal preparations from developing safflower (Carthamus tinctorius L.) cotyledons. Biochem J. 1990, 272: 23-29.PubMed CentralView ArticlePubMedGoogle Scholar
- Radwan O, Mouzeyar S, Venisse JS, Nicolas P, Bouzid MF: Resistance of sunflower to the biotrophic oomycete Plasmopara halstedii is associated with a delayed hypersensitive response within the hypocotyls. Journal of Experimental Botany. 2005, 56 (Suppl 420): 2683-2693.View ArticlePubMedGoogle Scholar
- Sackston WE: Downy mildew of sunflower. The downy mildews. Edited by: Spencer DM. 1981, Academic Press, London, UK, 545-575.Google Scholar
- Gundaev AI: Basic principles of sunflower selection. Genetic Principles of Plant Selection, Nauka, Moscow. 1971, 417-465. (Transl. Department of the Secretary of State, Ottawa, Canada, 1972)Google Scholar
- Ivanov P, Petakov D, Nikolova V, Pentchev E: Sunflower breeding for high palmitic acid content in the oil. Proceedings of the 12th International Sunflower Conference, Paris. 1988, International Sunflower Association, 463-465.Google Scholar
- Voinnet O: RNA silencing: small RNAs as ubiquitous regulators of gene expression. Curr Opin Plant Biol. 2002, 5 (Suppl 5): 444-451.View ArticlePubMedGoogle Scholar
- Tissier AF, Marillonnet S, Klimyuk V, Patel K, Torres MA, Murphy G, Jones JD: Multiple independent defective suppressor-mutator transposon insertions in Arabidopsis: a tool for functional genomics. Plant Cell. 1999, 11 (Suppl 10): 1841-1852.PubMed CentralView ArticlePubMedGoogle Scholar
- Barkley NA, Wang ML: Application of TILLING and EcoTILLING as reverse genetic approaches to elucidate the function of genes in plants and animals. Curr Gen. 2008, 9: 212-226. 10.2174/138920208784533656.View ArticleGoogle Scholar
- Till BJ, Reynolds SH, Greene EA, Codomo CA, Enns LC, Johnson JE, Burter C, Odden AR, Young K, Taylor NE, Henikoff JG, Comai L, Henikoff S: Large-scale discovery of induced point mutations with high-throughput TILLING. Genome Res. 2003, 13: 524-530. 10.1101/gr.977903.PubMed CentralView ArticlePubMedGoogle Scholar
- Comai L, Henikoff S: TILLING: practical single-nucleotide mutation discovery. Plant J. 2006, 45: 684-694. 10.1111/j.1365-313X.2006.02670.x.View ArticlePubMedGoogle Scholar
- Hohmann U, Jacobs G, Jung C: An EMS mutagenesis protocol for sugar beet and isolation of non-bolting mutants. Plant Breeding. 2005, 124: 317-321. 10.1111/j.1439-0523.2005.01126.x.View ArticleGoogle Scholar
- Gottschalk W, Wolf G: Induced mutations in plant breeding. Monograph on Theoretical and Applied Genetics 7. 1983, Springer, Berlin Heidelberg New YorkGoogle Scholar
- Caldwell DG, McCallum N, Shaw P, Muehlbauer GJ, Marshall DF, Waugh R: A structured mutant population for forward and reverse genetics in Barley (Hordeum vulgare L.). Plant J. 2004, 40: 143-50. 10.1111/j.1365-313X.2004.02190.x.View ArticlePubMedGoogle Scholar
- Menda N, Semel Y, Peled D, Eshed Y, Zamir D: In silico screening of a saturated mutation library of tomato. Plant J. 2004, 38: 861-872. 10.1111/j.1365-313X.2004.02088.x.View ArticlePubMedGoogle Scholar
- Xin Z, Wang ML, Barkley NA, Burow G, Franks C, Pederson G, Burke J: Applying genotyping (TILLING) and phenotyping analyses to elucidate gene function in a chemically induced sorghum mutant population. BMC Plant Biol. 2008, 8: 103-10.1186/1471-2229-8-103.PubMed CentralView ArticlePubMedGoogle Scholar
- Oleykowski CA, Mullins CRB, Godwin AK, Yeung AT: Mutation detection using a novel plant endonuclease. Nucleic Acids Res. 1998, 26: 4597-4602. 10.1093/nar/26.20.4597.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang B, Wen X, Kodali NS, Oleykowski CA, Miller CG, Kulinski J, Besack D, Yeung JA, Kowalski D, Yeung AT: Purification, cloning and characterization of the Cel I nuclease. Biochemistry. 2000, 39: 3533-3541. 10.1021/bi992376z.View ArticlePubMedGoogle Scholar
- Talamè V, Bovina R, Sanguineti MC, Tuberosa R, Lundqvist U, Salvi S: TILLMore, a resource for the discovery of chemically induced mutants in barley. Plant Biotechnol J. 2008, 6: 477-485. 10.1111/j.1467-7652.2008.00341.x.View ArticlePubMedGoogle Scholar
- Gottwald S, Bauer P, Komatsuda T, Lundqvist U, Stein N: TILLING in the two-rowed barley cultivar 'Barke' reveals preferred sites of functional diversity in the gene HvHox1. BMC Res Notes. 2009, 2: 258-10.1186/1756-0500-2-258.PubMed CentralView ArticlePubMedGoogle Scholar
- Himelblau E, Gilchrist EJ, Buono K, Bizzell C, Mentzer L, Vogelzang R, Osborn T, Amasino RM, Parkin IAP, Haughn GW: Forward and reverse genetics of rapid-cycling Brassica oleracea. Theor Appl Genet. 2009, 118: 953-961. 10.1007/s00122-008-0952-7.View ArticlePubMedGoogle Scholar
- Till B, Cooper J, Thai TH, Colowit P, Greene EA, Henikoff S, Comai L: Discovery of chemically induced mutations in rice by TILLING. BMC Plant Biol. 2007, 7: 19-10.1186/1471-2229-7-19.PubMed CentralView ArticlePubMedGoogle Scholar
- Suzuki T, Eiguchi M, Kumamaru T, Satoh H, Matsusaka H, Moriguchi K, Nagato Y, Kurata N: MNU-induced mutant pools and high performance TILLING enable finding of any gene mutation in rice. Mol Genet Genomics. 2008, 279: 213-223. 10.1007/s00438-007-0293-2.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.