<figure><img src="image_1.png" /></figure>

HAL
open science

# Development and evaluation of a cucumber TILLING population

Rina Fraenkel, Irina Kovalski, Christelle Troadec, Abdelhafid Bendahmane, Rafael Perl-Treves

<figure><img src="image_2.png" /></figure>

## To cite this version:

Rina Fraenkel, Irina Kovalski, Christelle Troadec, Abdelhafid Bendahmane, Rafael Perl-Treves. Development and evaluation of a cucumber TILLING population. BMC Research Notes, 2014, 7, pp.846. 10.1186/1756-0500-7-846. doi:10.1186/1756-0500-7-846

HAL Id: hal-02630077

https://hal.inrae.fr/hal-02630077

Submitted on 27 May 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_3.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Development and evaluation of a cucumber TILLING population

Rina Fraenkel<fnref n="1" />, Irina Kovalski<fnref n="1" />, Christelle Troadec², Abdelhafid Bendahmane² and Rafael Perl-Treves<fnref n="1" />*

## Abstract

**Background:** Ordered collections of mutants serve as invaluable tools in biological research. TILLING (Targeting Induced Local Lesions IN Genomes) provides an efficient method to discover, in mutagenized populations, the possible phenotypes controlled by gene sequences whose function is unknown. This method can replace transgenic techniques for the functional validation of cloned genes, especially in the case of transformation-recalcitrant plants such as cucumber.

**Results:** We report the development of a TILLING cucumber population, generated by EMS mutagenesis in the Poinsett76 genetic background. The population was evaluated by screening for morphological mutations, and a range of developmental, pigmentation and spontaneous lesion mutants were observed. Suitability for detecting single nucleotide polymorphism in selected genes has been tested by screening a sample of amplicons, with detection rate of 1 SNP in ~1 Mbp.

**Conclusion:** The population described in this Research Note represents a useful asset in cucumber research, to be exploited for forward genetic screens and functional genomics purposes.

**Keywords:** *Cucumis sativus*, EMS mutagenesis, Mutant screen, TILLING

## Background

Reverse genetic strategies are required to attribute a function, or phenotype, to the wealth of genes whose sequences are known, but their precise role in metabolism or development is unknown. The most common reverse genetic methods in plants involve T-DNA insertions or RNAi silencing, however, both methods require an efficient plant transformation platform and long term efforts to generate large arrays of transgenic lines. TILLING (Targeting Induced Local Lesions IN Genomes) is a relatively inexpensive, non-transgenic technique for interrogating the function of a given sequence; it is suitable for any species, regardless of genome size [1-3]. To generate a saturated TILLING population, wild-type seeds are typically treated with EMS (ethyl methanesulfonate), a mutagen that saturates the genome with point mutations, mostly G/C to A/T transitions [4]. The treated seeds are germinated, and the resulting M₁ plants self pollinated to collect M₂ seeds. DNA samples from mutagenized seed pools are

screened by PCR amplification of the gene of interest, followed by digestion with an endonuclease that cleaves at mismatched sites. If a particular M₂ family carries a point mutation in the amplified fragment, heteroduplex DNA carrying a mismatch at the mutated site will form during replication of mutant and wild-type alleles present together in the particular family-sample. Genetic analysis of the segregating family, *i.e.*, sequencing individuals from the suspect family and their phenotypic inspection, allows the discovery of mutant phenotypes and sheds light on the gene's function.

In a TILLING population, mutations are spread randomly in the genome; a spectrum of mutated alleles, including weak ones that are desirable for studying essential genes, can be obtained, and can even be used to infer structure-function relations of the encoded protein [5,6]. The method has been refined by Bendahmane and co-workers for different crops, including *Cucumis melo*, the melon [7,8]; see also [9], pea [10] and tomato [5]. The use of an Arabidopsis endonuclease, ENDO-1, rather than the commonly used CEL-1 nuclease from celery, rendered the screening of molecular mismatches particularly efficient [7,11,12].

<footnote>* Correspondence: rafi.perl@gmail.com</footnote>

<footnote n="1">¹The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan, Israel</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
Cucumber (*Cucumis sativus*) is a species of economic importance, whose genome has been sequenced recently [13-15]. However, being transformation-recalcitrant, functional genomics studies in cucumber have lagged behind. Here we report on a TILLING population that could enhance the cucumber genomic toolkit for future research. Studying a sample of the population, we show the range of morphological mutations that is seen, and demonstrate the occurrence of nucleotide substitutions in a few gene fragments.

## Results and discussion

### Mutagenesis and multiplication

To select the optimal EMS concentration that is likely to produce many mutations but will not heavily affect germination and fertility, we performed two calibration experiments, in which seeds of cucumber 'Poinsett76' were exposed to 0-2.5% EMS. Figure 1 displays the results of an experiment, in which 100 M₀-seed aliquots were treated with five concentrations (0, 1, 1.5, 2, 2.5% EMS, respectively). Germination rates were recorded and ranged between 95% (untreated seeds) and 80%, being moderately affected by the mutagen. We also recorded the rate of seedlings exhibiting somatic mutations, i.e., leaf distortion and mosaicism (dark/light variegation) on the first and second true leaves. Such rates increased gradually from 3% total ratio of apparent abnormalities in untreated seedlings, to 97% and 100% in the 2% and 2.5% treatments, respectively. We selected 1.5%-2% EMS as optimal concentrations that mildly affected germination but had a

substantial proportion of visible somatic damage. Subsamples of these plants were grown further and we noted that in most plants the new leaves that developed were normal, and the mosaic symptoms recorded in the first leaves did not extend to subsequent leaves. Apparently, the somatic mutations did not affect the shoot apical meristem but only the first and second leaf-primordia in the mutagenized seed; only a small proportion displayed persistent mutations at the M₁ generation. A similar recovery of melon M₁ plants from growth inhibition following 1-2% EMS mutagenesis was reported by Dahmani-Mardas et al. [8]. For pea, however, a lower 0.25% concentration (20 mM) was selected [10], with less than 30% of the plants setting seeds at higher concentrations; in Arabidopsis, high mutation rates were achieved by 0.25-0.5% EMS [16] and the authors suggested that DNA repair mechanisms could be responsible for inter-specific differences in mutation yields.

Thirty seedlings from the 1%, 1.5% and 2% treatments were transferred to the greenhouse for further growth and fertility assessment. We concluded that fertility following self pollination was similarly good at all three EMS concentrations, with 87% of the plants treated with either 1.5% or 2% mutagen producing >100 seeds. This agrees with Dahmani-Mardas et al. [8] who reported that 3% EMS, but not 1-2%, significantly affected melon fertility.

To produce the TILLING population, we applied 1.5 or 2% EMS to Poinsett76 cucumber seeds in three different batches. A total of 1200 M₁ plants were grown to full maturity in a farmer's net-house in Netiv Ha'asara, Israel Southern coastal region, and in Bar-Ilan University net-house, and self-pollinated to collect M₂ seeds. This resulted in ~1000 M₂ families with adequate yield of 50-300 seeds. Six seedlings per M₂ family were germinated, and a mixed DNA sample from four of them was prepared. The first 768 families were arranged in 96 well plate wells (see Methods) and used to evaluate the population.

## Phenotypic evaluation of the population

We utilized the M₂ seedling samples (6 individuals per family) grown for DNA extraction to record mutant phenotypes that were observable at the cotyledon and first leaf stages. Whereas we cannot exclude that some phenotypes could be due to environmental variation or seedling physiology – especially differences in size or germination ability - in many cases, the phenotype was apparent in two individuals of the same family, demonstrating the likely heritability of the putative mutation. Table 1 summarizes the phenotypes that we recorded and Figure 2 provides several examples. The most frequent seedling phenotypes included post-germination lethality, dwarfism, and spontaneous lesions in the cotyledons. The latter class was unlikely to result from

<figure><img src="image_4.png" /></figure>

**Figure 1 Calibration of EMS treatment.** Aliquots of 100 seeds were treated with 0, 1, 1.5, 2 and 2.5% EMS and sown in germination trays. Percent germination was recorded after 9 and 16 days, and seedlings exhibiting somatic mutations in the cotyledons, first and second true leaf were counted and expressed as percent of the fully germinated seedlings. Mutations included smaller or distorted cotyledons or leaves, as well as dark-light leaf patterns. Standard errors of the percentage of the final germination or mutated seedlings rates was computed as 100× [b × (1-b)/n]⁰.⁵, where b is the proportion of a given class, and n – the total number of plants [17].
<table><caption>Table 1 Major mutant phenotypes recovered by visual screening of the TILLING population at the seedling stage</caption><thead><tr><td>Phenotype</td><td>No. of families displaying phenotype</td><td>Description</td><td>Representative family</td></tr></thead><tbody><tr><td>Seedling lethality</td><td>19</td><td>Seedling dies, or fails to develop a root after germination</td><td>144, 157</td></tr><tr><td>Dwarf</td><td>10</td><td>Cotyledons smaller, short hypocotyl</td><td>422</td></tr><tr><td>Necrotic lesions</td><td>15</td><td>Spontaneous necrotic spots appear on cotyledons</td><td>38</td></tr><tr><td>Albino</td><td>3</td><td>White cotyledons</td><td>164</td></tr><tr><td>Yellow leaf</td><td>7</td><td>True leaves pale-green or yellow, or mosaic green and yellow</td><td>218</td></tr><tr><td>Tall seedlings</td><td>5</td><td>Hypocotyl >2 cm taller than wild type</td><td>180</td></tr><tr><td>Glabrous</td><td>1</td><td>Cotyledons lack trichomes</td><td>16</td></tr><tr><td>Non-serrated leaf</td><td>3</td><td>Leaf edge appears smooth</td><td>189</td></tr><tr><td>Dark, narrow cotyledon</td><td>1</td><td>Dark green, narrow cotyledons</td><td>411</td></tr><tr><td>Fused cotyledons</td><td>1</td><td>Cotyledons fused together</td><td>714</td></tr><tr><td>Small leaf</td><td>7</td><td>True leaf very small</td><td>17</td></tr><tr><td>Distorted leaf/cotyledon</td><td>5</td><td>Irregular organ shape</td><td>160, 213</td></tr><tr><td>Total mutants/ families screened</td><td>77/768</td><td></td><td></td></tr></tbody></table>

Six seedlings per M₂ family were sown in trays and inspected at the cotyledon-first true leaf stage. The number of families and examples of specific families that segregate for a given phenotype-class are indicated. About 10% of the families exhibited morphological alterations at the seedling stage.

response to pathogens in the growth chamber, as it appeared sporadically, often with two seedlings per family showing the lesions while the rest of the sample and the adjacent seedlings in the tray lacked them. Such mutations could identify genes involved in disease-resistance signaling. Mutants with distorted cotyledons and/or leaves, albino cotyledons and yellow or mosaic-yellow leaves were also apparent, as well as taller seedlings, dark-narrow and fused cotyledons. In total, 10% of the families exhibited a morphological alteration.

Several seedlings that exhibited prominent abnormalities were grown to maturity, and outstanding phenotypes could be seen among the mature plants as well. These included persistent virescent or yellow-leaf character (Family 424), a fasciated plant with floral organ abnormalities and organ-fusions (Family 176), and more; Figure 3 provides a few examples. This demonstrates that the population is a rich source of morphological and developmental mutations that could be tapped using forward-genetic schemes. Since each M₁ plant and the descendant M₂ family harbors multiple point mutations, genetic analysis will be required to discern them and correlate a specific mutant phenotype with a single molecular event.

## Molecular evaluation of test-genes

To check whether we could interrogate our population with a given gene sequence and recover point mutations in selected gene fragments, we chose six genes that yielded visible phenotypes when mutated in other plants. Primers were designed to amplify the more conserved

parts of the coding sequence; in some cases, intronic sequences were also included. Figure 4A depicts the amplification scheme of four genes, for which point mutations were recovered. Table 2 specifies the primers that were used.

Cucumber *phytoene desaturase-3* (*Pds-3*), gene accession number Csa002881, encodes a carotene biosynthetic enzyme. Inactive alleles could give rise to an albino phenotype [18]. We screened three amplicons and identified four mutated families, two in amplicon B and two in D. We sequenced and verified the mutations in families 53 and 254, and each carried a different C-to-T substitution. In family 254, three plants were heterozygous and two had the wild type allele; in family 53, two were wild-type, two homozygous for the mutation and four were heterozygous (Table 3, Figure 4). All four mutations mapped to introns present in the amplicons, and no mutant phenotypes were observed. This underscores the importance of targeting exons and avoiding introns when performing a TILLING screen.

The *Female* sex-determining gene encodes an ACC synthase enzyme that controls the differentiation of female flowers and abolishes male flowers [19]. We detected a single nucleotide substitution in our screen that segregated in the expected 1:2:1 ratio, and mapped it to intron 2 of the gene (Table 3, Figure 4A). As expected from its intronic location, the mutation had no phenotypic effect.

*Ramosus-3* and *Ramosus-4* are genes that affect apical dominance and were extensively studied in pea, as well as other plant species [20]; *rms* mutants often show extensive branching. We designed and screened amplicons for each
<figure><img src="image_5.png" /><figcaption>**Figure 2** Selection of morphological mutations segregating in M₂ mutated families at the seedling stage. Family 164: albino cotyledons, green leaf. Another mutation for pointed cotyledons seems to segregate as well. Family 411: dark, seemingly anthocyanin-enriched, narrow cotyledons. Family 38: spontaneous necrosis of cotyledons and leaf. Family 218: yellow-light green pigmentation. Family 147: dwarf phenotype, small dark cotyledons.</figcaption></figure>

gene, and recovered one mutant in *rms-3* and two in *rms-4*. All three substitutions mapped to the coding sequence (Table 3, Figure 4A). Two of them involved silent mutations and no phenotype was observable. The third mutation, in the single large exon of the gene (Family 53) caused a threonine-to-leucine substitution. In the M₂ generation, we did not observe progeny with increased branching, or other phenotype that we could associate with certainty with the Thr256Ile substitution. We have analyzed the

respective protein sequence using SIFT (Sorting Intolerant From Tolerant, [21]) and found that substituting Thr to Ile at this position is predicted to be tolerated and thus, unlikely to affect the phenotype. In addition to the above four genes, we screened also a cucumber homolog of the tomato *self-pruning* gene for determinate growth habit (cucumber *sp*, Csa010707) and a MADS box homolog, *CUM-1* (Csa00068; Tables 2 and 3), but no mutations were recovered for these two genes.
<figure><img src="image_6.png" /><figcaption>**Figure 3** Selection of morphological mutants at the mature plant stage. Family 38: deeper-lobed leaf (right, compared to wild-type leaf on the left), and smaller male and female flowers (bottom, compared to wt flowers on top). Family 424: pale green plant. Family 164: darker leaf with shallow lobes (left, compared to wt leaf on the right). Family 176: fasciated, sessile inflorescences with reiterated organs, multiple petalled flowers and a branched ovary. Family FX: "cauliflower" mutant with arrested-development, reiterating, inflorescence with dense trichomes and a lanceolate leaf.</figcaption></figure>

We have recovered, by screening six genes, in a 768 families-sample of the population (of the ~1000 presently available), a total of 8 mutants, 6 of which were confirmed by sequencing the M₂ progeny; the other two were estimated to be in the intron, judging from their endonuclease product sizes, and were not analyzed further. The mutations were recovered by screening 13 amplicons that totaled 10,766 nucleotides (Table 3). The average mutation frequency in the collection is, therefore, 8 mutations per 10,766 × 768 bp of sequence screened across the population, i.e., 1 in 1.03 Mbp. If we ignore amplicon borders, where endonuclease digestion would pass undetected, frequency would be 1 mutation in 859 kb. Such density is

fair, albeit lower than the density observed by for pea (1/200 kbp, [10]), *Arabidopsis* (1/300 kbp, [16] and melon (1/573 bp, [8]); it is similar to the observed rate in barley (1/1 Mbp, [22]) and allows for Reverse Genetics screening of the population. In a parallel study, another TILLING population has been prepared in the cucumber 'Beit-Alpha' background, and a similar mutation density (1/1.14 Mbp) was recovered [23].

In the amplicons that we screened, we could compare SNP rates in exons (comprising 60% of the amplicons) vs introns (40%). Of the eight mutations, five were located in introns and three in exons. Of the six sequenced mutations, five were C to T, and one was a G to A transition,
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 4** (See legend on next page.)</figcaption></figure>
(See figure on previous page.)

**Figure 4 Screen for nucleotide substitutions in selected gene fragments.** A. Gene model and amplification schemes of the cucumber *phytoene desaturase-3* gene, *Female* ACC synthase, *ramosus-3* and *ramosus-4* homologous genes. The approximate positions of the amplicons screened by TILLING are indicated, each delimited by two pairs of nested primers (arrows indicate the internal amplicons). The mutations that were verified by sequencing (from families 53, 254, 48, 540, 928) are indicated; in family 53, two independent point mutations were recovered in two of the genes, respectively. B. Chromatogram of the C to T mutation (read by a reverse primer as G to A), discovered in the *PDS-3* gene, in family 53. Top: wild type sequence in cultivar Poinsett76. Below: plant 53–13 is homozygous for the mutation, plant 53–20 is heterozygous.

which are the two common products of EMS mutagenesis [16]. The ratio of intronic mutations exceeded the ratio of exonic ones by a factor of 2.5 (= 60/40 × 5/3). Moreover, two of three exonic mutations were silent; it has been estimated that about half of the non silent (missense) mutations are likely to affect protein activity [24]. Thus, if one wishes to recover mutated alleles with a higher

probability of affecting the phenotype, population size should be increased to ~3000 plants and inclusion of introns in the amplicons screened should be minimized, although, when exons are small, introns are hard to avoid. The present SNP sampling was not aimed at recovering mutant phenotypes, the goal being a mere assessment of DNA-level mutation rates.

<table><caption>**Table 2 Primer pairs used to generate amplicons for the TILLING screen of six genes, as detailed in the methods**</caption><thead><tr><th>Gene</th><th>Amplicon</th><th>External primers</th><th>Internal primers (M13 tag bold)</th></tr></thead><tbody><tr><td rowspan="5">PDS-3</td><td rowspan="2">pdsB</td><td>Pds3for1, CACAGATGACATTCTTCCAAT</td><td>pds3for3, <b>CACGACGTTGTAAAACGACT</b>CTAACATACCCATAGG</td></tr><tr><td>Pds3rev2, CCTAGTTCTACCCTTTGTTCTTGG</td><td>pds3rev4, <b>ATAACAATTTCACACAGG</b>GTITTCATGCTGGCTGCC</td></tr><tr><td rowspan="2">pdsD</td><td>Pds3for5, GGAAATTTGTCTCAACATGTGTGC</td><td>Pds3for7, <b>CACGACGTTGTAAAACGACT</b>CTCCAAACTAGTGAC</td></tr><tr><td>Pds3rev6, CTTGTGCCACATGGCTAGAATAG</td><td>Pds3rev8, <b>GATAACAATTTCACACAGG</b>CTGCCGGTGATGCTGG</td></tr><tr><td>pdsF</td><td>Pds3for9, AAGGGGCTCGACTGTTCAGAAA</td><td>Pds3for11, <b>CACGACGTTGTAAAACGAC</b>CACTAGAAGACTCAGC</td></tr><tr><td rowspan="4">Female</td><td rowspan="2">acsA</td><td>Pds3rev10, CTGGTAGTGATTCTCGGTTTCA</td><td>Pds3rev12, <b>GATAACAATTTCACACAGG</b>CATCCACCGAAGTAGA</td></tr><tr><td>AcsFor1, GAACTATCTACCATATTCCAACC</td><td>AcsFor2, <b>CACGACGTTGTAAAACGAC</b>GTACCTATATACCTCACCTCAACAT</td></tr><tr><td rowspan="2">acsB</td><td>AcsRev3, CACCAACTCGAAAACCTGGGAGCC</td><td>AcsRev1, <b>GGATAACAATTTCACACAGG</b>CCTCGTCTTCGTTACTCCTCTCCT</td></tr><tr><td>AcsFor5, TTGATAGAGATTTGAAATGGAGA</td><td>AcsFfor6, <b>CACGACGTTGTAAAACGAC</b>CCGGAGTTGAGATTGTGCCAATTC</td></tr><tr><td rowspan="5">RMS-3</td><td rowspan="2">rms3A</td><td>AcsRev7F, CCGAGTGCACTTTTCTTTTTTC</td><td>AcsRev5, <b>GGATAACAATTTCACACAGG</b>ATTCCCCCAAATATGGATGATG</td></tr><tr><td>Rms3for1, GTGTCACCGTGCATGCAATTGCCG</td><td>Rms3for2, GTACACTTCAAATCATAAACGGCTG</td></tr><tr><td rowspan="2">rms3B</td><td>Rms3Rev2, TTCATCTCAGTTTTCCTACCTAAT</td><td>Rms3Rev1, GGGTTACAAACGCTGGCCTTC</td></tr><tr><td>Rms3for6, ATTGCAAACATAGCCATCAAAATC</td><td>Rms3for7, GTCAAAATTATCATTTCTACGCAGG</td></tr><tr><td rowspan="5">RMS-4</td><td rowspan="2">rms4A</td><td>Rms3rev6, CAAGTAAAAACACAGCTCTCAACCT</td><td>Rms3Rev5, GAGTGGATGCTATTCTTTTCGATG</td></tr><tr><td>Rms4for5, CTCTCTCCGTTGCTAAGACAAACCC</td><td>Rms4for1, TCCGATTACTGTATCTTCCTGCTCG</td></tr><tr><td rowspan="2">rms4B</td><td>Rms4rev6, CAAACTCACCATTGTTCTCAAACCC</td><td>Rms4rev1, CGACAGCGATTCAAGCCCTTGACAA</td></tr><tr><td>Rms4for4, TTCTCGTGGCCAATCCTCTGAC</td><td>Rms4for2, CCATTACCGAGGCTTGCCCTAACCT</td></tr><tr><td rowspan="3">rmsC</td><td>Rms4rev5, GCCCCACAAATCATTACCACTGCAT</td><td>Rms4rev2, GCCCCACAAATCATTACCACTGCAT</td></tr><tr><td>Rms4G3for7, GGATGGAACTATGGTGGATATA</td><td>Rms4for3, CCAGGTACTCCACCGACGCTGATTG</td></tr><tr><td>Rms4G3rev7, CCCATGAAAGATTGTGAAATCACA</td><td>Rms4rev3, GTGATACAGCTAATCTCAAAGTAAC</td></tr><tr><td rowspan="3">Cum-1</td><td rowspan="2">Cum1B</td><td>Cum1for1, GCTCTTTTCCTCATCAGGTTAGTG</td><td>Cum1for3, <b>ACGACGTTGTAAAACGACT</b>CTCCTGCTCATTCCAC</td></tr><tr><td>Cum1rev2, GTATACACCAAACTGAGAACCAG</td><td>Cum1rev4, <b>GATAACAATTTCACACAGG</b>TACAACCCAGATTCCC</td></tr><tr><td>Cum1D</td><td>Cum1for5, GATATCAATTAAACCATGCGGGC</td><td>Cum1for7, <b>CACGACGTTGTAAAACGACT</b>TGGTATGAAATGGGGG</td></tr><tr><td rowspan="4">SP</td><td rowspan="2">spB</td><td>Cum1rev6, GATTATCGGTTTCATCTCCATGG</td><td>Cum1rev8, <b>GATAACAATTTCACACAGG</b>CCCAATTCAGACCTTC</td></tr><tr><td>SPfor1, GGACAGCACAAGAAAAGGTCAC</td><td>SPfor2, <b>CACGACGTTGTAAAACGACT</b>TGCCCTCTCTCTGCT</td></tr><tr><td rowspan="2">spC</td><td>SPrev1, CACATCATTTCTTGCCAATTGTC</td><td>SPrev2, <b>GATAACAATTTCACACAGG</b>GGGCATCTTTTGCAAC</td></tr><tr><td>SPfor1, GGACAGCACAAGAAAAGGTCAC</td><td>SPfor3, <b>CACGACGTTGTAAAACGAC</b>CAACCAATTCCCCAAAC</td></tr><tr><td></td><td></td><td>SPrev1, CACATCATTTCTTGCCAATTGTC</td><td>SPrev3, <b>GATAACAATTTCACACAGG</b>CACTCTCTCTCACCAG</td></tr></tbody></table>

The bold portion of the internal primers corresponds to the M13 tail that binds a third pair of universal primers that were fluorescently labeled. The corresponding internal primers of *Rms-3* and *Rms-4* were directly labeled.
<table><caption>Table 3 Nucleotide substitutions recovered in query genes by TILLING screen</caption><thead><tr><th>Gene</th><th>Amplicon</th><th>Amplicon size, bp</th><th>% GC content</th><th>% Exons</th><th>Mutated family</th><th>SNP</th><th>Diagnostic CAPS</th><th>Location/ substitution</th><th>M3 Progeny analyzed: WT/WT: Heteroz.: mut/mut</th></tr></thead><tbody><tr><td rowspan="4">Phytoene desaturase-3 (PDS-3), Csa002881</td><td rowspan="2">pdsB</td><td rowspan="2">622</td><td rowspan="2">36.5</td><td rowspan="2">55</td><td>53</td><td>C5310T</td><td></td><td>Intron 8</td><td>2: 4: 2</td></tr><tr><td>188</td><td>nd</td><td></td><td>Intron 8</td><td>nd</td></tr><tr><td rowspan="2">pdsD</td><td rowspan="2">1198</td><td rowspan="2">34.5</td><td rowspan="2">23</td><td>254</td><td>C3997T</td><td>-</td><td>Intron 6</td><td>2: 3: 0</td></tr><tr><td>70</td><td>nd</td><td></td><td>Intron 6</td><td>nd</td></tr><tr><td rowspan="2">ACC synthase (F), Csa012150</td><td>acsA</td><td>1111</td><td>39</td><td>70</td><td>none</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>acsB</td><td>1025</td><td>41</td><td>96</td><td>48</td><td>C594T</td><td>Dral</td><td>Intron 2</td><td>4: 9: 5</td></tr><tr><td rowspan="2">Ramosus-3 (RMS-3), Csa010158</td><td>Rms3A</td><td>513</td><td>47</td><td>72</td><td>none</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Rms3B</td><td>639</td><td>43</td><td>70</td><td>540</td><td>G1376A</td><td>-</td><td>Gln191Gln</td><td>3: 4: 0</td></tr><tr><td rowspan="4">Ramosus-4 (RMS-4), Csa003326</td><td>Rms4A</td><td>1222</td><td>51</td><td>86</td><td>53</td><td>C767T</td><td>-</td><td>Thr256Ile</td><td>3: 10: 3</td></tr><tr><td>Rms4B</td><td>929</td><td>46</td><td>100</td><td>928</td><td>C1258T</td><td>Sacl / Xbal</td><td>Leu420Leu</td><td>7: 7: 0</td></tr><tr><td>Rms4C</td><td>798</td><td>45</td><td>84</td><td>none</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Cucumber MADS1 (Cum1), Csa000681</td><td>Cum1B</td><td>575</td><td>35</td><td>41</td><td>none</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="2">Self-pruning (sp), Csa010707</td><td>Cum1D</td><td>702</td><td>34</td><td>30</td><td>none</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>SpB</td><td>704</td><td>29.8</td><td>39</td><td>none</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Total</td><td></td><td>10,766</td><td>39.7</td><td>63</td><td>8</td><td></td><td></td><td></td><td></td></tr></tbody></table>

Genes are indicated by name and by accession numbers (cucumber genome project, http://cucumber.genomics.org.cn). Amplicon size (calculated between internal primer pairs) is shown, as well as the GC composition and the ratio of coding sequence (exons) to total amplicon length. Mutated position is determined according to the genomic sequence, from the start codon (ATG), and its location, either in an intron or in the protein coding-sequence, is indicated. Nd - non determined. CAPS marker in ACC synthase gene: the wild type amplicon is digested by Dral, mutant amplicon is uncut. CAPS marker in RMS4: wild type amplicon is digested by Sacl (and also Xbal), no restriction in the mutant. A small number of M2 progeny was genotyped to demonstrate inheritance of the nucleotide substitution. The total sequence screened (10,766 bp) was calculated by summing up all the internal amplicons, and detracting the overlapping regions found between the Acs (F) and Rms-4 amplicons (see Figure 4).

By sequencing a small sample of progeny of mutated families in the M2 generation we confirmed the Mendelian inheritance of the nucleotide substitutions. In two cases we facilitated screening by developing a diagnostic CAPS marker. In three cases we did not recover the homozygous-mutant class (Table 3). Since we did not expect the respective mutations to be lethal (none of them affects the protein product of the gene), this segregation could be due to the small sample of progeny, or to pollination of a heterozygous mutant plant with wild-type pollen; this, in turn, could result from outcrossing with a neighbor plant, or from a chimeric plant composition following seed mutagenesis. In such cases, self pollinating of heterozygous individuals can readily provide the desired homozygous mutants.

## Conclusions

A cucumber mutant population in the Poinsett76 background was successfully constructed using 1.5-2% EMS treatment. It provides a rich source for morphological mutant phenotypes that can be recovered by forward genetic screens. It can also be efficiently screened by the reverse genetics TILLING approach, to recover mutations in target genes. With the entire cucumber genome sequence available, our study thus provides a valuable tool for functional genomics in cucumber. Data and seeds for collaborative

research can be obtained by agreement from the corresponding author.

## Methods

### Preparation of the TILLING population

Poinsett76 cucumber seeds were incubated at room temperature for 15 hrs, in 3 volumes of freshly-prepared EMS solution (Sigma M0880), in 0.2 M phosphate buffer (pH = 7), with gentle stirring. After extensive washing (5 × 30 min in tap water at room temperature with stirring), seeds were sown in “Speedling” germination trays and transplanted in 10 L pots in the greenhouse at the 2–3 leaf stage. Plants were grown under standard agro-nomic conditions, self-pollinated, and M2 seeds were harvested (50–300 per plant). From each M2 family, six seeds were sown, and after two weeks, 4 young leaf discs were sampled from four individuals of the family (8 mm diameter, one disc per individual) and pooled for DNA extraction by the CTAB method [25]. Uniform DNA concentrations were obtained by suspending in 200 µl water and samples were kept at –80°C. Individual family-samples of the first 768 families were arrayed in eight 96 well plates, and DNA aliquots were diluted 10 times (to obtain a 400 µl volume, ~10 ng DNA/µl) and kept at –20°C. DNA pools were prepared (8 families /pool) and arrayed in a single 96-well pooling-plate representing
768 families. Such sampling strategy was possible due to
the highly sensitive TILLING screen that allows detection
of mutant alleles in the eight-family x four individuals
mixed sample.

Target gene amplification and screening

Amplification included two steps of nested PCR. The
first involved a pair of external unlabeled primers,
22–25 nt-long, 40-50% recommended GC content, planned
to generate a fragment up to 1.2 kb in size. Primers were
reacted with 2 µL of the pooled genomic DNA (~20 ng) for
30 PCR cycles. The second round of PCR was performed
using 1 µL of the first reaction products as template,
with two pairs of internal primers present together in
the same reaction. The first primer pair in the second
round comprised a 3' region that specifically binds the
PCR products of the first round (having a 40-50%
recommended GC content), while the 5' region of the
primer represents a universal M13 sequence “tail”, resulting
in a 35 nt-long primer. The last primer pair (M13F700, 5'
CACGACGTTGTAAAACGAC and M13R800, 5' GGATA
ACATTTCACACAGG) is fluorescently marked (IRDye™
700 and IRDye™ 800) and matches the M13 sequence tags
of the second primer pair. For the *rms-3* and *rms-4* genes,
the second pair of primers was fluorescently labeled and
M13-labeled primers were not required. The reaction in-
cluded 35 cycles, the first 10 performed at the recom-
mended annealing temperature of the specific primer pair,
then 25 cycles at a lower annealing temperature, 50°C,
required for the universal primers [3,12]. Finally, heterodu-
plex molecules were formed by a temperature gradient of
94°C to 8°C (–0.1°C/sec).

The products were digested for 20 min at 45°C with
EndoI endonuclease and separated on a polyacrylamide
sequencing gel using a LICOR 4300 machine. The
pooled 96 DNA samples were screened, and positive
family-pools displaying amplicon digestion were decon-
voluted to identify the positive M2 family within the
pool. The family's DNA sample (that had been prepared
by mixing four individual progeny) was sequenced, to
look for a mixed peak, indicating a mutation at a pos-
ition predicted by the endonuclease digestion pattern. A
few progeny of the family were sown and their amplicon
sequenced, to confirm the mutation and correlate it with
a possible phenotype. In a few cases, the mutation gen-
erated a restriction site polymorphism; in that case,
genotyping was performed by restriction digestion of the
PCR product, followed by agarose gel electrophoresis
and ethidium bromide staining.

Abbreviations

CAPS: Cleaved amplified polymorphic sequence; EMS: Ethyl
methanesulfonate; PCR: Polymerase chain reaction; SNP: Single nucleotide
polymorphism; TILLING: Targeting Induced Local Lesions IN Genomes.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

RF prepared and organized the mutagenized population and performed the morphological screen and part of the molecular screen. IK assisted in growing and sampling of the populations and validated the mutations. CT supervised and performed the TILLING screen of six genes. AB provided methodologies and supervision of molecular screen. RPT supervised the preparation and screening of the population and wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgements

R. Fraenkel was supported by an M.Sc. Fellowship of Bar-Ilan University.
We thank Mr. David Levy for technical assistance.

Author details

¹The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan, Israel. ²INRA-CNRS, UMR1165, Unité de Recherche en Génomique Végétale, Evry, France.

Received: 20 July 2014 Accepted: 18 November 2014

Published: 26 November 2014

References

1. Comai L, Henikoff S: TILLING: practical single-nucleotide mutation discovery. *Plant J* 2006, 45(4):684–694.

2. McCallum CM, Comai L, Greene EA, Henikoff S: Targeting induced local lesions in genomes (TILLING) for plant functional genomics. *Plant Physiol* 2000, 123(2):439–442.

3. 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(2):480–484.

4. Anderson P: *Mutagenesis. Methods Cell Biol* 1995, 48:31–58.

5. Jones MO, Piron-Prunier F, Marcel F, Piednoir-Barbeau E, Alsadon AA, Wahb-Allah MA, Al-Doss AA, Bowler C, Bramley PM, Fraser PD, Bendahmane A: Characterisation of alleles of tomato light signalling genes generated by TILLING. *Phytochemistry* 2012, 79:78–86.

6. Perry J, Brachmann A, Welham T, Binder A, Charpentier M, Groth M, Haage K, Markmann K, Wang TL, Parniske M: TILLING in Lotus japonicus identified large allelic series for symbiosis genes and revealed a bias in functionally defective ethyl methanesulfonate alleles toward glycine replacements. *Plant Physiol* 2009, 151(3):1281–1291.

7. Nieto C, Piron F, Dalmais M, Marco CF, Moriones E, Gomez-Guillamon ML, Truniger V, Gomez P, Garcia-Mas J, Aranda MA, Bendahmane A: EcoTILLING for the identification of allelic variants of melon eIF4E, a factor that controls virus susceptibility. *BMC Plant Biol* 2007, 7. Article Number: 34, doi:10.1186/1471-2229-7-34.

8. Dahmani-Mardas F, Troadec C, Boualem A, Leveque S, Alsadon AA, Aldoss AA, Dogimont C, Bendahmane A: Engineering melon plants with improved fruit shelf life using the TILLING approach. *PLoS One* 2010, 5(12). Article Number: e15776, doi:10.1371/journal.pone.0015776.

9. Tadmor Y, Katzir N, Meir A, Yaniv-Yaakov A, Sa'ar U, Baumkoler F, Lavee T, Lewinsohn E, Schaffer A, Burger J: Induced mutagenesis to augment the natural genetic variability of melon (*Cucumis melo* L.). *Isr J Plant Sci* 2007, 55(2):159–169.

10. Dalmais M, Schmidt J, Le Signor C, Moussy F, Burstin J, Savois V, Aubert G, Brunaud V, de Oliveira Y, Guichard C, Thompson R, Bendahmane A: *UTILLdb*, a *Pisum sativum* in silico forward and reverse genetics tool. *Genome Biol* 2008, 9(2). Article Number: R43, doi:10.1186/gb-2008-9-2-r43.

11. Triques K, Sturbois B, Gallais S, Dalmais M, Chauvin S, Clepet C, Aubourg S, Rameau C, Caboche M, Bendahmane A: Characterization of *Arabidopsis thaliana* mismatch specific endonucleases: application to mutation discovery by TILLING in *pea*. *Plant J* 2007, 51(6):1116–1125.

12. Triques K, Piednoir E, Dalmais M, Schmidt J, Le Signor C, Sharkey M, Caboche M, Sturbois B, Bendahmane A: Mutation detection using ENDO1: application to disease diagnostics in humans and TILLING and Eco-TILLING in plants. *BMC Mol Biol* 2008, 9. Article Number: 42, doi:10.1186/1471-2199-9-42.

13. Huang S, Li R, Zhang Z, Li L, Gu X, Fan W, Lucas WJ, Wang X, Xie B, Ni P, Ren Y, Zhu H, Li J, Lin K, Jin W, Fei Z, Li G, Staub J, Kilian A, van der Vossen EAG, Wu Y, Guo J, He J, Jia Z, Ren Y, Tian G, Lu Y, Ruan J, Qian W, Wang M,
et al: The genome of the cucumber, Cucumis sativus L. *Nat Genet* 2009, **41**(12):1275–U1229.

14. Yang L, Koo D-H, Li Y, Zhang X, Luan F, Havey MJ, Jiang J, Weng Y: Chromosome rearrangements during domestication of cucumber as revealed by high-density genetic mapping and draft genome assembly. *Plant J* 2012, **71**(6):895–906.

15. Woycicki R, Witkowicz J, Gawronski P, Dabrowska J, Lomsadze A, Pawelkowicz M, Siedlecka E, Yagi K, Plader W, Seroczynska A, Smiech M, Gutman W, Niemirowicz-Szczytt K, Bartoszewski G, Tagashira N, Hoshi Y, Borodovsky M, Karpinski S, Malepszy S, Przybecki Z: *The genome sequence of the North-European Cucumber (Cucumis sativus L.) unravels evolutionary adaptation mechanisms in plants*. *PLoS One* 2011, **6**(7). Article Number: e22728, doi:10.1371/journal.pone.0022728.

16. 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**(2):731–740.

17. Downie NM, Heath RH: *Basic Statistical Methods*. 5th edition. New-York and London: Harper and Row Publishers; 1983.

18. Qin G, Gu H, Ma L, Peng Y, Deng XW, Chen Z, Qu L-J: *Disruption of phytoene desaturase gene results in albino and dwarf phenotypes in Arabidopsis by impairing chlorophyll, carotenoid, and gibberellin biosynthesis*. *Cell Res* 2007, **17**(5):471–482.

19. Trebitsh T, Staub JE, Oneill SD: Identification of a 1-aminocyclopropane-1-carboxylic acid synthase gene linked to the female (F) locus that enhances female sex expression in cucumber. *Plant Physiol* 1997, **113**(3):987–995.

20. Johnson X, Brcich T, Dun EA, Goussot M, Haurogne K, Beveridge CA, Rameau C: Branching genes are conserved across species. *Genes controlling a novel signal in pea are coregulated by other long-distance signals*. *Plant Physiol* 2006, **142**(3):1014–1026.

21. Ng PC, Henikoff S: SIFT: predicting amino acid changes that affect protein function. *Nucleic Acids Res* 2003, **31**:3812–3814.

22. 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**(1):143–150.

23. Boualem A, Fleurier S, Troadec C, Audigier P, Kumar AP, Chatterjee M, Alsadon AA, Sadder MT, Wahb-Allah MA, Al-Doss AA, Bendahmane: Development of a *Cucumis sativus* TILLING platform for forward and reverse genetics. *PLoS One* 2014, **9**(5):e97963.

24. Markiewicz P, Kleina LG, Cruz C, Ehret S, Miller JH: Genetic-sudies of the Lac repressor.14. Analysis of 4000 altered Escherichia-coli Lac repressors reveals essential and nonessential residues, as well as spacers which do not require a specific sequence. *J Mol Biol* 1994, **40**(5):421–433.

25. Rogers SO, Bendich AJ: Extraction of DNA from milligram amounts of fresh, herbarium and mummified plant-tissues. *Plant Mol Biol* 1985, **5**(2):69–76.

doi:10.1186/1756-0500-7-846
Cite this article as: Fraenkel et al.: Development and evaluation of a cucumber TILLING population. BMC Research Notes 2014 7:846.

Submit your next manuscript to BioMed Central
and take full advantage of:

* Convenient online submission

* Thorough peer review

* No space constraints or color figure charges

* Immediate publication on acceptance

* Inclusion in PubMed, CAS, Scopus and Google Scholar

* Research which is freely available for redistribution

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central