Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Engineering Melon Plants with Improved Fruit Shelf Life Using the TILLING Approach

  • Fatima Dahmani-Mardas ,

    Contributed equally to this work with: Fatima Dahmani-Mardas, Christelle Troadec, Adnane Boualem

    Affiliation Unité de Recherche en Génomique Végétale, UMR1165 ERL8196 INRA-UEVE-CNRS, Evry, France

  • Christelle Troadec ,

    Contributed equally to this work with: Fatima Dahmani-Mardas, Christelle Troadec, Adnane Boualem

    Affiliation Unité de Recherche en Génomique Végétale, UMR1165 ERL8196 INRA-UEVE-CNRS, Evry, France

  • Adnane Boualem ,

    Contributed equally to this work with: Fatima Dahmani-Mardas, Christelle Troadec, Adnane Boualem

    Affiliation Unité de Recherche en Génomique Végétale, UMR1165 ERL8196 INRA-UEVE-CNRS, Evry, France

  • Sylvie Lévêque,

    Affiliation Unité de Génétique et Amélioration des Fruits et Légumes, INRA UR1052, Montfavet, France

  • Abdullah A. Alsadon,

    Affiliation Department of Plant Production, College of Food and Agricultural Sciences, King Saud University, Riyadh, Saudi Arabia

  • Abdullah A. Aldoss,

    Affiliation Department of Plant Production, College of Food and Agricultural Sciences, King Saud University, Riyadh, Saudi Arabia

  • Catherine Dogimont,

    Affiliation Unité de Génétique et Amélioration des Fruits et Légumes, INRA UR1052, Montfavet, France

  • Abdelhafid Bendahmane

    bendahm@evry.inra.fr

    Affiliations Unité de Recherche en Génomique Végétale, UMR1165 ERL8196 INRA-UEVE-CNRS, Evry, France, Department of Plant Production, College of Food and Agricultural Sciences, King Saud University, Riyadh, Saudi Arabia

Abstract

Background

Fruit ripening and softening are key traits that have an effect on food supply, fruit nutritional value and consequently, human health. Since ethylene induces ripening of climacteric fruit, it is one of the main targets to control fruit over ripening that leads to fruit softening and deterioration. The characterization of the ethylene pathway in Arabidopsis and tomato identified key genes that control fruit ripening.

Methodology/Principal Findings

To engineer melon fruit with improved shelf-life, we conducted a translational research experiment. We set up a TILLING platform in a monoecious and climacteric melon line, cloned genes that control ethylene production and screened for induced mutations that lead to fruits with enhanced shelf life. Two missense mutations, L124F and G194D, of the ethylene biosynthetic enzyme, ACC oxidase 1, were identified and the mutant plants were characterized with respect to fruit maturation. The L124F mutation is a conservative mutation occurring away from the enzyme active site and thus was predicted to not affect ethylene production and thus fruit ripening. In contrast, G194D modification occurs in a highly conserved amino acid position predicted, by crystallographic analysis, to affect the enzymatic activity. Phenotypic analysis of the G194D mutant fruit showed complete delayed ripening and yellowing with improved shelf life and, as predicted, the L124F mutation did not have an effect.

Conclusions/Significance

We constructed a mutant collection of 4023 melon M2 families. Based on the TILLING of 11 genes, we calculated the overall mutation rate of one mutation every 573 kb and identified 8 alleles per tilled kilobase. We also identified a TILLING mutant with enhanced fruit shelf life. This work demonstrates the effectiveness of TILLING as a reverse genetics tool to improve crop species. As cucurbits are model species in different areas of plant biology, we anticipate that the developed tool will be widely exploited by the scientific community.

Introduction

Melon (Cucumis melo) belongs to the Cucurbitaceae family that contains about 800 species mainly distributed in tropical and subtropical regions [1]. Cucurbitaceae includes several other economically important cultivated plants, such as cucumber (C. sativus), watermelon (Citrullus lanatus), squash and pumpkin (Cucurbita spp.). Aside tomato (Solanum lycopersicum), watermelons, cucumbers and melons are the most cultivated vegetable species (http://faostat.fao.org). For decades, melon has played key roles in the field of plant molecular biology and plant physiology, serving as an excellent model organism for investigating sex determination [2], [3] and ripening processes [4], [5]. Melon is also a reference model to investigate vascular fluxes, as xylem and phloem saps can be readily collected [6], [7].

Melon is a diploid species (2n = 2x = 24) with an estimated genome size of 450 Mb [8]. Research efforts have been undertaken in order to gather genetic and genomic informations on melon, such as EST sequences [9], BAC libraries [10], [11], oligo-based microarrays [12] and high resolution genetic maps [13][15]. Melon has truly entered in the genomic era by the sequencing of the highly syntenic cucumber genome [16]. With the completion of the melon genome sequencing project in the near future [17], a major challenge is to determine gene functions. In plants, the most common techniques to produce altered or loss of function mutations are T-DNA or transposon insertional mutagenesis [18] and RNA interference [19]. However, because these methods are mainly based on the ability of a given plant to be transformed, their usefulness as general reverse genetics methods is limited to very few plant species. On the other hand, ethyl methanesulfonate (EMS) mutagenesis is a simple method to saturate a genome with mutations [20]. TILLING (Targeting Induced Local Lesions IN Genomes) combines advantages of random chemical mutagenesis and high throughput mutation discovery methods [21] and generates allelic series of the targeted genes which makes it possible to dissect the function of the protein as well as to investigate the role of lethal genes. This technique has been successfully applied to a large variety of organisms including plants and animals [20], [22][33].

To establish a reverse genetics platform in melon and to identify novel alleles of agronomic importance, we have set up a melon TILLING platform and performed a screen for mutations in genes that control fruit ripening. Fruit ripening is a process characterized by a number of biochemical and physiological processes that alter fruit firmness, colour, flavour, aroma, and texture [34], [35]. Ripening patterns are conserved across fleshy fruit species. Characterizations of homologous genes involved in fruit ripening in different species suggest that genetic mechanisms are also conserved [36]. In the fleshy fruit model, tomato, the plant hormone, ethylene, plays a central role in ripening [37][40]. Several fruit ripening mutants have been characterized in tomato [41]. The rin (ripening-inhibitor), nor (non-ripening) and cnr (colorless non-ripening) genes were shown to act upstream the ethylene signaling in ripening [41][43]. Mutations in the ethylene receptor, the Nr (Never-ripe) gene, were also shown to affect fruit ripening [44]. Ethylene biosynthesis requires the conversion of aminocyclopropane-1-carboxylic acid (ACC) to ethylene by the ACC oxidase (ACO). ACC oxidases are encoded by multigene families in plants [45], [46] and have been described to be involved in ripening, growth, and development. In melon, CmACO1 silencing inhibits fruit ripening and extends fruit shelf life [47].

Melon is an attractive model for investigating fruit ripening. Unlike tomato, melon has two different ripening patterns, climacteric and non climacteric. To investigate further the role of ethylene in melon fruit ripening, we have developed a reference EMS mutant population under controlled conditions and established a TILLING platform. Then, we screened for mutations in 11 genes, mainly involved in ripening processes, and characterized CmACO1 TILLING mutants. This work yields a missense mutation in CmACO1 that inhibits fruit ripening and extends fruit storage life. The use of TILLING as a translational research tool is discussed.

Results

Production of CharMono mutant collection

The melon inbred line CharMono is a monoecious climacteric Charentais type cultivar (Cucumis melo L. subsp. melo var cantalupensis) that bears male flowers on the main stem and female flowers on axillary branches. The success of the TILLING approach relies on the construction of high quality mutant libraries. Ideally the mutant population must produce a mutation frequency that is conducive to high-throughput screening but is below a threshold that causes extensive sterility and plant development alteration. To optimize the mutagenesis, we first conducted a “kill-curve” analysis on batches of 100 seeds, using a dose range from 1% to 3% EMS (Table 1). Most M1 treated plants exhibited growth retardation at early seedling stage, but all of them recovered, except for the 3% EMS fraction for which only 40 plants were recovered. The M1 plants were assessed for seed production and seed germination. Low seed production was observed at 3% EMS with 37% of the plants producing less than 40 seeds per fruit. At 2% EMS, only 9% of the plants yielded less than 40 seeds per fruit (Table 1). Thus, the highest EMS doses allowing good seed set and seed germination, 1.5% and 2%, were retained and tested on large batches of seeds. The EMS treated seeds were sown in soil and seedlings were grown to fruit maturity in insect-proof plastic tunnels to avoid cross-pollination. Female flowers were hand pollinated with male flowers from the same plants, bagged and the fruits left to develop to maturity. M2 seeds were harvested from individual M1 plants. We also produced M2 seeds from 617 and 40 M1 plants treated with 1% and 3% EMS, respectively. In total, 4260 M2 seed stocks were collected.

thumbnail
Table 1. Impact of EMS concentration on M2 seed setting and germination.

https://doi.org/10.1371/journal.pone.0015776.t001

To assess the quality of the mutagenesis, we investigated the rate of appearance of depigmentation mutants at the cotyledon stage and developmental alterations at fruit maturity. Albino and chlorotic plants, the most frequently observed phenotype in mutagenized populations [22], [48], occurred at the rate of 1,3% of the M2 families. This rate is similar to the rate reported for other well characterized mutant collections [49][51]. The most commonly observed developmental phenotypes were related to cotyledon number and morphology, dwarfism, shoot branching and plant architecture. At the flowering stage, we observed mutants with flower sex type transition, abnormal flower morphology and architecture. At the fruiting stage, altered fruit shape and size, ovary number and flesh color and thickness variations were also observed. Examples of the observed phenotypes are shown in Figure 1.

thumbnail
Figure 1. Examples of mutant phenotypes observed in CharMono mutant collection.

Plant affected in cotyledon color (A) or number (B), leaf morphology (C), carpel number (E), reduced fruit size (left fruit in F) and fleshless fruit (right fruit in F). Fruit from control wild type plant is shown in (D). Scale bars: 5 cm.

https://doi.org/10.1371/journal.pone.0015776.g001

Melon TILLING platform

To set up the TILLING platform, DNA samples were prepared from 4023 M2 families, each representing an independent M1 family and organized in pools of 8 families as described previously [25]. One key factor in TILLING is the availability of an annotated genomic sequence of the gene to be tilled, which is facilitated by the ongoing sequencing of melon genome (International Cucurbit Genomics Initiative, www.icugi.org), the availability of a large melon EST collection [9] and the 95% sequence similarity over coding regions between melon and the available cucumber genome sequence [16]. Primer design and PCR amplifications were performed as described previously [22]. Mutations were detected, in the amplified targets, using the mismatch-specific endonuclease ENDO1 [21]. Individual mutant lines were identified following a pool deconvolution step, and the mutated bases were identified by sequencing.

To assess the quality of CharMono mutant collection and to estimate the mutation density, we screened for induced mutations in 11 genes related to fruit quality. In total, we identified and confirmed by sequencing 134 induced mutations in 18,3 kb total length of tilled amplicons (Table 2). Most induced mutations were as expected, G/C to A/T transitions [20], with the exception of the three following mutations, G/C to C/G, T/A to A/T and C/G to A/T (Table S1), which have also been reported in other TILLING studies [24], [26], [52], [53].

thumbnail
Table 2. Tilled genes and mutation frequency in CharMono mutant population.

https://doi.org/10.1371/journal.pone.0015776.t002

Sequence analysis of the exonic induced mutations showed that 31,3% were silent, 65,1% were missense, 2,4% were stop and 1,2% were splice junction mutations (Figure S1, Table 3). The number of stop mutations was lower than the CODDLE predicted proportions. In contrast, silent, missense, and splice junction mutations were recovered in an equivalent proportion as predicted (Table 3). As many tilled amplicons harbored non coding segments, some recovered mutations could potentially affect the splicing or the stability of the mRNA, such impacts are unpredictable. Thus, non coding mutations were not characterized further. In contrast, the many non-synonymous mutations discovered were investigated to dissect the function of the proteins as they may lead to gain- or loss-of-function phenotypes. We estimated the mutation density based on Greene et al. [20] to 1 mutation every 573 kb. This observed density is proportional to the EMS dose used, except for the 1.5% EMS fraction showing a lower mutation frequency may be due to some experimental variations (Table 4).

thumbnail
Table 3. Expected and observed frequencies of induced mutation types in tilled gene-coding regions.

https://doi.org/10.1371/journal.pone.0015776.t003

thumbnail
Table 4. Effect of EMS dose on mutation frequency in CharMono mutant collection.

https://doi.org/10.1371/journal.pone.0015776.t004

Characterization of CmACO1 TILLING mutants

Fruit softening is a major factor that determines fruit quality and shelf life. In melon and tomato, silencing of enzymes or regulators of the ethylene biosynthesis pathway inhibits fruit ripening and extend fruit storage life [47], [54][58]. To identify melon lines with improved fruit shelf life, we screened for mutations in the CmACO1 enzyme that catalyses the last step of ethylene biosynthesis (Figure 2A). In this TILLING screen, we identified seven independent point mutations among which two led to L124F and G194D missense mutations. The L124F and G194D changes occurred in a highly conserved amino acid positions and may therefore affect the function of the protein (Figure 2B, [59]). However, X-ray crystallography studies determined that L124 is located in the α-5 helix of the protein away from the active site and thus, it is predicted to not affect the function of the protein (Figure 2C–E, amino acid indicated in green; [60]). In contrast, the residue G194 is located in the β-7 strand, one of the eight β strands (β-4 to β-11) that form the distorted double-stranded β helix (DSBH or jellyroll) core of the active site, common to all members of the 2-OG-oxygenases [60][65]. The residue G194 is located in the inner face of the active site, and thus, a mutation at this position is predicted to affect the function of the protein (Figure 2C-E, amino acid indicated in red; [60]).

thumbnail
Figure 2. Sequence and structural analysis of CmACO1.

(A) Schematic diagram of CmACO1 gene structure. The numbers indicate the size of the 4 exons (filled boxes) and the 3 introns (black lines) in bps. Intronic, silent and the L124F and G194D missense mutations are represented by black, grey, green and red triangles, respectively. (B) Amino acid alignment of CmACO1 (Cucumis melo, Q04644) and homologous proteins from Cucumis sativus (Cs, BAA33377), Arabidopsis thaliana (At, NP_171994), Solanum lycopersicum (Sl, P05116), Petunia x hybrida (Ph, Q08506), Vitis vinifera (Vv, XP_002273430), Medicago truncatula (Mt, AAL35971), Populus trichocarpa (Pt, XP_002320487) and Oryza sativa (Os, NP_001063330). Numbers above the alignment indicate the amino acid positions along the CmACO1 protein. L124F and G194D EMS-induced mutations are shown above the alignment in green and red, respectively. (C–D) 3D structure model of CmACO1. (C) Superposition of the Petunia ACO structure determined by X-ray crystallography [60], indicated in pink and the 3D model of CmACO1, indicated in white. The melon model was determined using the Geno3D server (http://geno3d-pbil.ibcp.fr). (D and E) Zoom in of the ACO1 active site. Yellow sticks represent Ile171, His177, Asp179 and His234 residues involved in the Fe (II) cofactor binding (orange ball). Blue sticks represent the ligating phosphate or sulfate ion. Amino acids L124 and G194 are represented in green and red sticks, respectively. (F–K) Fruit-related traits of the L124F and G194D EMS mutants. Data for the duration between the pollination and the fruit ripening (days, F), firmness (Durofel index, G), shape (ratio length/width in mm, H), soluble sugars (° Brix, I) and flesh (J) and rind (K) color (b value) are shown. G124D mutant and its relative wild type control and L124F and its relative wild type control are indicated by red, black, green and open bars, respectively. Asterisks indicate significant differences.

https://doi.org/10.1371/journal.pone.0015776.g002

To test whether induced mutations in CmACO1 could affect melon fruit ripening, L124F and G194D TILLING mutants were evaluated for different fruit traits (Figure 2F–K). The L124F mutant did not show any difference from its wild type control for the evaluated fruit traits. This is consistent with the conservative L124F mutation, two amino acids with hydrophobic side chains, and the position of the L124 residue away from the active site of the protein (Figure 2F–K, white and green bars). In contrast, the duration between the pollination and the fruit ripening as well as the fruit firmness were significantly increased in the G194D mutant compared to the wild type control (Figure 2F and G, black and red bars). G194D mutation was also associated with a significant change in the fruit rind color leading to a delayed fruit yellowing (Figure 2K). Fruit shape, soluble sugar content and flesh color showed no differences between the G194D mutant line and the wild type control (Figure 2H–J). No differences were also observed for the fruit peduncle abscission layer that is activated at fruit maturity and lead the fruit to drop from the plant.

To confirm the observed phenotypes, G194D heterozygote mutant lines were selfed and 40 homozygote mutant and wild type segregants were assessed for the above fruit traits. Homozygous wild type plants for this locus had no alteration of any fruit traits described above and behaved as the CharMono wild type parent. As predicted, plants homozygous for the G194D mutation presented an increased duration between the pollination and the fruit ripening, a better fruit firmness and a delayed fruit yellowing without any changes in fruit shape, peduncle abscission layer, soluble sugar content and flesh color. These data confirmed that the observed phenotypes resulted from G194D induced mutation in CmACO1. To investigate the robustness of the observed phenotypes, we repeated the phenotypic characterization of the mutant lines in two different locations, Montfavet and Saint Rémy de Provence, and confirmed the observed phenotypes (Figure S2).

Discussion

In the genomic era and in the absence of efficient tools for homologous recombination, TILLING has become an obligate technology to dissect gene function as well as to engineer alleles of agronomic importance in crops. To set up the melon TILLING platform, we first developed a reference EMS mutant collection under controlled conditions. To assess the quality of the mutagenesis, we estimated the rate of occurrence of chlorotic and albino phenotypes. The observed frequency of 1,3%, comparable to well characterized mutant collections, confirmed the quality of the mutagenesis [49][51]. Genomic DNA was prepared from the mutant lines and organized in pools for bulked screening using the mismatch specific endonuclease, Endo1 [21]. The value of this melon mutant collection was then confirmed by screening for mutation in 11 genes. On average, we identified 8 alleles per tilled kilobase, calculated from the tilled 18,5 kb and the 134 identified alleles. We also estimated the overall mutation rate as one mutation every 573 kb in our CharMono mutant collection. This mutation frequency is two fold lower than the rate of one mutation per 200 kb reported for Pisum sativum [22], two fold higher than the rate of one mutation per megabase reported for barley [53] and equivalent to the rate of mutations described in tomato [24], [25]. A much more saturated mutation density has been observed in polyploid species which withstand much higher doses of EMS without obvious impact on plant survival (1/40 kb in tetraploid wheat and 1/24 kb in hexaploid wheat [66]).

Fruit ripening and softening are key traits that determine the shelf life of fleshy climacteric fruits. Post-harvest losses can reach 70% in developing countries because of the lack of post-harvest infrastructure to store and to retail the commodities. A wide list of approaches, ranging from the application of growth regulators to delay ripening, fruit storage in controlled atmospheres and breeding or genetic engineering of leader lines have been used to control fruit ripening and softening. Ethylene has been identified as the major hormone that controls fruit ripening and softening in climacteric fruits. Investigations in the model species as well as in crops have accumulated considerable evidences at the genetic, physiological, biochemical and molecular levels that pointed out the ethylene function at various levels. This includes ethylene biosynthesis, its perception by the target cells through specific receptors, signal transduction pathway involving both positive and negative regulators and finally regulation of target genes expression by transcription factors such as ethylene response factors [46]. The final step of ethylene biosynthesis in plants is catalyzed by the ACC oxidase which converts ACC to ethylene. ACC oxidases belong to a multigene family and members have been shown to be highly expressed during fruit ripening in climacteric fruits. ACC oxidase role in ripening was demonstrated by co-suppression of the fruit ACC oxidases in fleshy fruit plants. The ACC oxidases silenced fruits exhibited non-softening phenotypes and thus, an extended shelf life [47], [54][58].

The TILLING of CmACO1 identified a mutation at the conserved G194 residue, located in the active site of the enzyme. Detailed phenotypic characterization of the TILLING mutants showed that the G194D mutation mimic the phenotype of CmACO1 antisense line [47]. The G194D mutant line showed longer fruit maturation, enhanced fruit firmness and delayed fruit yellowing. As reported for the antisense ACC oxidase melon lines [47], there are also no difference in the fruit sugar content, shape or flesh color. Besides, we noticed that the peduncle abscission layer that is activated at fruit maturity in wild type plant was also observed on the G194D mutant fruits in contrast with the antisense melon lines that failed to drop even at a very late developmental stage (65 days) [47]. This difference could be explained by the antisense strategy that can not, with complete conviction, exclude that other related homologous enzymes are affected in the peduncle abscission layer.

In conclusion, to investigate fruit maturation and fruit softening we tilled a list of candidate genes and isolated key mutants that enhanced melon fruit quality (data not shown). Through the TILLING approach we knocked out a specific gene, CmACO1 and pointed out an important amino acid, Glycine at the position 194, for the fruit maturation, not reported previously. Moreover, we created a new allele for the melon post-harvest management (maturation duration) and transport (firmness), which is an economically attractive trait for breeders. This is even more important regarding the transgenic strategies that face strong concern in public acceptance. The characterization of the other mutants will likely deliver alleles that may be exploited in melon breeding. Cucurbits are also an important model species in many key areas of plant research, including sex determination [2], [3], [67], fruit maturation [4], [5] and the investigation of vascular trafficking of molecules [6], [7]. Hence, by opening the TILLING platform to the scientific community, we hope to fulfill the expectations of both crop breeders and scientists who are using melon as their model of study.

Materials and Methods

Plant material and EMS mutagenesis

Experiments were carried out using the melon inbred line CharMono, a monoecious climacteric Charentais type cultivar (Cucumis melo L. subsp. melo var cantalupensis). Mature seeds from CharMono line were immersed in bottles containing 250 ml of EMS diluted at right concentration in deionized water, and placed on a rotary shaker overnight (16 h) at 24°C in the dark. The EMS treatment was stopped by addition of 100 ml of Na2SO3 at 0,1 M. Treated seeds were then transferred into a new solution of Na2SO3 0,1M for 15 minutes. The EMS solution was then removed and seeds were washed extensively. Treated seeds were sown in soil and grown under insect-proof plastic tunnels according to standard melon agronomic practice. At the end of the fruit ripening phase, M2 seeds were collected from individual M1 plants and stored.

Genomic DNA extraction and pooling

Sixteen melon leaf discs (diameter 10 mm) from eight individual plants per M2 family were collected in 96-well plates containing 2 steel beads (4 mm) per well, and tissues were ground using a bead mill. Genomic DNA was isolated using the Dneasy 96 Plant Kit (Qiagen, Hilden, Germany). DNAs were quantified on a 1% agarose gel using λ DNA (Invitrogen, Carlsbad, USA) as a concentration reference. DNA samples were then diluted ten fold and pooled eight fold in a 96-well format.

PCR amplification and mutation detection

PCR amplification is based on nested-PCR. The first PCR amplification is a standard PCR reaction using target-specific primers (Table S2) and 4 ng of melon genomic DNA. One µl of the first PCR served as a template for the second nested PCR amplification, using combination of specific primers carrying M13 tail and M13 universal primers, M13F700 (5′-CACGACGTTGTAAAACGAC-3′) and M13R800 (5′-GGATAACATTTCACACAGG-3′), labelled at the 5′end with infra-red dyes IRD700 and IRD800 (LI-COR®, Lincoln, Nebraska, USA), respectively. Mutation detection was carried out as described previously [22]. The identity of the mutations was determined by sequencing. TILLING request should be addressed to the corresponding author.

Sequence analysis tools

CODDLE (Codons Optimized to Discover Deleterious Lesions, http://www.proweb.org/coddle/) was used to identify regions of the target gene in which G/C to A/T transitions are most likely to result in deleterious effects on the protein. PARSESNP (Project Aligned Related Sequences and Evaluate SNPs, http://www.proweb.org/parsesnp/) was used to illustrate the distribution of mutations within the gene, and to indicate the nature of each single mutation. SIFT (Sorting Intolerant from Tolerant, http://sift.jcvi.org/www/SIFT_seq_submit2.html) was used to predict the impact of the mutation on the protein. Multiple sequence alignment of full-length protein sequences was performed using ClustalW (http://www.ebi.ac.uk/Tools/clustalw2).

Protein structure modeling

The CmACO1 three-dimensional structures were generated using the Geno3D server (http://geno3d-pbil.ibcp.fr). Superposition of the petunia ACO structure (1W9Y.pdb) determined by X-ray crystallography [60], and our CmACO1 model, was carried out and visualized using Chimera (http://www.cgl.ucsf.edu/chimera).

TILLING mutant phenotyping

Mutant and wild type plants were grown in two locations, Montfavet and Saint Rémy de Provence, in France, under greenhouse conditions. Fruits were harvested at maturity, according to the harvest indicators usually observed in Charentais type fruits, namely: the development of abscission layer, yellow rind and death of the first leaf beside the peduncle. At maturity, the duration between pollination and harvesting was evaluated. Mature fruit shape was measured as the ratio of the polar diameter (peduncle to blossom end), referred to as fruit length, to the equatorial width of the fruit (measured halfway between the peduncle and blossom end). Soluble sugar content of the flesh was measured with a hand refractometer. Firmness was measured in two different areas of the flesh after cutting the fruit lengthwise, using a Durofel electronic penetrometer or with Gullimex penetrometer equipped with an 8 mm diameter probe. Rind and flesh colors were recorded using a Minolta colorimeter at two different sites per fruit. The (b) parameter was selected for monitoring changes in the yellow hue [68].

Supporting Information

Figure S1.

Graphic representations of genes screened for mutations. Black boxes represent the exons. Lanes linking exons indicate introns. Dashed lines in red indicate the genomic regions screened for mutations. Triangles pointing up indicate mutations in coding regions, whereas those pointing down indicate mutations in intron-exon splicing sites. Red, black and grey triangles represent alterations causing truncations, missense and silent mutations, respectively. Mutations in introns, 5′ and 3′UTR are not shown.

https://doi.org/10.1371/journal.pone.0015776.s001

(TIF)

Figure S2.

Multi-location evaluation of fruit traits for the G194D EMS mutant. (A–F) Fruit-related traits tested at INRA Avignon experimental station in Montfavet, France and (G–H) at a second experimental station at Saint Rémy de Provence, France. Data for the duration between the pollination and the fruit ripening (days, A and G), firmness (Durofel index, B and kg/0,5 cm2, H), shape (ratio length/width in mm, C and I), soluble sugars (° Brix, D and J) and flesh (E and K) and rind (F) color (b value) are shown. G124D mutant and its relative wild type control are indicated by red and black bars, respectively. Blue stars indicate significant difference. Values represent the means of measurements obtained for at least 10 fruits from different plants.

https://doi.org/10.1371/journal.pone.0015776.s002

(TIF)

Table S1.

Missense, STOP and splicing induced mutations. Nucleotide substitution, amino acid changes and the predicted impact on the protein function are reported for each EMS-mutant line. Asterisks indicate non conventional EMS mutations. The amino acid substitution is predicted damaging if the score is < or  = 0.05 and tolerated if the score is >0.05.

https://doi.org/10.1371/journal.pone.0015776.s003

(XLS)

Table S2.

Primers used in this study. Asterisk indicate primers labeled with infra-red dyes IRD700 (forward primers) and IRD800 (reverse primers).

https://doi.org/10.1371/journal.pone.0015776.s004

(XLS)

Acknowledgments

The authors thank N. Giovinazzo, D. Besombes, V. Chareyron, P. Audigier, C. Lepage and the INRA Avignon experimental team for the plant care, seed production and technical assistance. We thank F. de Langen and D. Potey (Clause Company, Saint Rémy de Provence, France) for their contribution to the mutants phenotyping, J. Eleblu, B. Lasseur and C. Clepet for comments on the manuscript. We are grateful to the seed companies, ASL, Gautier Semences and Takii France for their help to create the melon mutagenized population.

Author Contributions

Conceived and designed the experiments: AB AB CD AAA AAA. Performed the experiments: FDM CT SL CD. Analyzed the data: AB AB FDM CT SL CD. Contributed reagents/materials/analysis tools: AB AB FDM CT SL CD. Wrote the paper: AB AB FDM.

References

  1. 1. Jeffrey C (1990) Systematics of the Cucurbitaceae: an overview. In: Bates DM RW, Jeffrey C, editors. Biology and utilization of the Cucurbitaceae. New-York, USA: Cornell University. pp. 3–9.
  2. 2. Boualem A, Fergany M, Fernandez R, Troadec C, Martin A, et al. (2008) A conserved mutation in an ethylene biosynthesis enzyme leads to andromonoecy in melons. Science 321: 836–838.
  3. 3. Martin A, Troadec C, Boualem A, Rajab M, Fernandez R, et al. (2009) A transposon-induced epigenetic change leads to sex determination in melon. Nature 461: 1135–1138.
  4. 4. Ezura H, Owino WO (2008) Melon, an alternative model plant for elucidating fruit ripening. Plant Science 175: 121–129.
  5. 5. Pech JC, Bouzayen M, Latché A (2008) Climacteric fruit ripening: Ethylene-dependent and independent regulation of ripening pathways in melon fruit Plant Science 175: 114–120.
  6. 6. Haritatos E, Keller F, Turgeon R (1996) Raffinose oligosaccharide concentrations measured in individual cell and tissue types in Cucumis melo L. leaves: implications for phloem loading. Planta 198: 614–622.
  7. 7. Gomez G, Torres H, Pallas V (2005) Identification of translocatable RNA-binding phloem proteins from melon, potential components of the long-distance RNA transport system. Plant J 41: 107–116.
  8. 8. Arumuganathan K, Earle ED (1991) Nuclear DNA content of some important plant species. Plant Molecular Biology Reporter 9: 208–218.
  9. 9. Gonzalez-Ibeas D, Blanca J, Roig C, Gonzalez-To M, Pico B, et al. (2007) MELOGEN: an EST database for melon functional genomics. BMC Genomics 8: 306.
  10. 10. van Leeuwen H, Monfort A, Zhang HB, Puigdomenech P (2003) Identification and characterization of a melon genomic region containing a resistance gene cluster from a constructed BAC library. Microcolinearity between Cucumis melo and Arabidopsis thaliana. Plant Mol Biol 51: 703–718.
  11. 11. Morales M, Orjeda G, Nieto C, van Leeuwen H, Monfort A, et al. (2005) A physical map covering the nsv locus that confers resistance to Melon necrotic spot virus in melon (Cucumis melo L.). Theor Appl Genet 111: 914–922.
  12. 12. Mascarell-Creus A, Canizares J, Vilarrasa-Blasi J, Mora-Garcia S, Blanca J, et al. (2009) An oligo-based microarray offers novel transcriptomic approaches for the analysis of pathogen resistance and fruit quality traits in melon (Cucumis melo L.). BMC Genomics 10: 467.
  13. 13. Perin C, Gomez-Jimenez M, Hagen L, Dogimont C, Pech JC, et al. (2002) Molecular and genetic characterization of a non-climacteric phenotype in melon reveals two loci conferring altered ethylene response in fruit. Plant Physiol 129: 300–309.
  14. 14. Fernandez-Silva I, Eduardo I, Blanca J, Esteras C, Pico B, et al. (2008) Bin mapping of genomic and EST-derived SSRs in melon (Cucumis melo L.). Theor Appl Genet 118: 139–150.
  15. 15. Deleu W, Esteras C, Roig C, Gonzalez-To M, Fernandez-Silva I, et al. (2009) A set of EST-SNPs for map saturation and cultivar identification in melon. BMC Plant Biol 9: 90.
  16. 16. Huang S, Li R, Zhang Z, Li L, Gu X, et al. (2009) The genome of the cucumber, Cucumis sativus L. Nat Genet 41: 1275–1281.
  17. 17. Gonzalez VM, Garcia-Mas J, Arus P, Puigdomenech P (2010) Generation of a BAC-based physical map of the melon genome. BMC Genomics 11: 339.
  18. 18. Alonso JM, Stepanova AN, Leisse TJ, Kim CJ, Chen H, et al. (2003) Genome-wide insertional mutagenesis of Arabidopsis thaliana. Science 301: 653–657.
  19. 19. Waterhouse PM, Graham MW, Wang MB (1998) Virus resistance and gene silencing in plants can be induced by simultaneous expression of sense and antisense RNA. Proc Natl Acad Sci U S A 95: 13959–13964.
  20. 20. Greene EA, Codomo CA, Taylor NE, Henikoff JG, Till BJ, et al. (2003) Spectrum of chemically induced mutations from a large-scale reverse-genetic screen in Arabidopsis. Genetics 164: 731–740.
  21. 21. Triques K, Sturbois B, Gallais S, Dalmais M, Chauvin S, et al. (2007) Characterization of Arabidopsis thaliana mismatch specific endonucleases: application to mutation discovery by TILLING in pea. Plant J 51: 1116–1125.
  22. 22. Dalmais M, Schmidt J, Le Signor C, Moussy F, Burstin J, et al. (2008) UTILLdb, a Pisum sativum in silico forward and reverse genetics tool. Genome Biol 9: R43.
  23. 23. Le Signor C, Savois V, Aubert G, Verdier J, Nicolas M, et al. (2009) Optimizing TILLING populations for reverse genetics in Medicago truncatula. Plant Biotechnol J 7: 430–441.
  24. 24. Minoia S, Petrozza A, D'Onofrio O, Piron F, Mosca G, et al. (2010) A new mutant genetic resource for tomato crop improvement by TILLING technology. BMC Res Notes 3: 69.
  25. 25. Piron F, Nicolai M, Minoia S, Piednoir E, Moretti A, et al. (2010) An induced mutation in tomato eIF4E leads to immunity to two potyviruses. PLoS One 5: e11313.
  26. 26. Cooper JL, Till BJ, Laport RG, Darlow MC, Kleffner JM, et al. (2008) TILLING to detect induced mutations in soybean. BMC Plant Biol 8: 9.
  27. 27. Till BJ, Reynolds SH, Greene EA, Codomo CA, Enns LC, et al. (2003) Large-scale discovery of induced point mutations with high-throughput TILLING. Genome Res 13: 524–530.
  28. 28. Till BJ, Reynolds SH, Weil C, Springer N, Burtner C, et al. (2004) Discovery of induced point mutations in maize genes by TILLING. BMC Plant Biol 4: 12.
  29. 29. Uauy C, Paraiso F, Colasuonno P, Tran RK, Tsai H, et al. (2009) A modified TILLING approach to detect induced mutations in tetraploid and hexaploid wheat. BMC Plant Biol 9: 115.
  30. 30. Stephenson P, Baker D, Girin T, Perez A, Amoah S, et al. (2010) A rich TILLING resource for studying gene function in Brassica rapa. BMC Plant Biol 10: 62.
  31. 31. Winkler S, Schwabedissen A, Backasch D, Bokel C, Seidel C, et al. (2005) Target-selected mutant screen by TILLING in Drosophila. Genome Res 15: 718–723.
  32. 32. Moens CB, Donn TM, Wolf-Saxon ER, Ma TP (2008) Reverse genetics in zebrafish by TILLING. Brief Funct Genomic Proteomic 7: 454–459.
  33. 33. Gilchrist EJ, O'Neil NJ, Rose AM, Zetka MC, Haughn GW (2006) TILLING is an effective reverse genetics technique for Caenorhabditis elegans. BMC Genomics 7: 262.
  34. 34. Brady CJ (1987) Fruit Ripening. Annual Review of Plant Physiology 38: 155–178.
  35. 35. Brummell DA, Harpster MH (2001) Cell wall metabolism in fruit softening and quality and its manipulation in transgenic plants. Plant Mol Biol 47: 311–340.
  36. 36. Adams-Phillips L, Barry C, Kannan P, Leclercq J, Bouzayen M, et al. (2004) Evidence that CTR1-mediated ethylene signal transduction in tomato is encoded by a multigene family whose members display distinct regulatory features. Plant Mol Biol 54: 387–404.
  37. 37. Giovannoni J (2001) Molecular Biology of Fruit Maturation and Ripening. Annu Rev Plant Physiol Plant Mol Biol 52: 725–749.
  38. 38. Giovannoni JJ (2004) Genetic regulation of fruit development and ripening. Plant Cell. 16. pp. S170–180.
  39. 39. Seymour G, Poole M, Manning K, King GJ (2008) Genetics and epigenetics of fruit development and ripening. Curr Opin Plant Biol 11: 58–63.
  40. 40. Barry C, Giovannoni JJ (2007) Ethylene and Fruit Ripening Journal of Plant Growth Regulation 26: 143–159.
  41. 41. Giovannoni JJ (2007) Fruit ripening mutants yield insights into ripening control. Curr Opin Plant Biol 10: 283–289.
  42. 42. Vrebalov J, Ruezinsky D, Padmanabhan V, White R, Medrano D, et al. (2002) A MADS-box gene necessary for fruit ripening at the tomato ripening-inhibitor (rin) locus. Science 296: 343–346.
  43. 43. Manning K, Tor M, Poole M, Hong Y, Thompson AJ, et al. (2006) A naturally occurring epigenetic mutation in a gene encoding an SBP-box transcription factor inhibits tomato fruit ripening. Nat Genet 38: 948–952.
  44. 44. Wilkinson JQ, Lanahan MB, Yen HC, Giovannoni JJ, Klee HJ (1995) An ethylene-inducible component of signal transduction encoded by never-ripe. Science 270: 1807–1809.
  45. 45. Alexander L, Grierson D (2002) Ethylene biosynthesis and action in tomato: a model for climacteric fruit ripening. J Exp Bot 53: 2039–2055.
  46. 46. Lin Z, Zhong S, Grierson D (2009) Recent advances in ethylene research. J Exp Bot 60: 3311–3336.
  47. 47. Ayub R, Guis M, Ben Amor M, Gillot L, Roustan JP, et al. (1996) Expression of ACC oxidase antisense gene inhibits ripening of cantaloupe melon fruits. Nat Biotechnol 14: 862–866.
  48. 48. Wu JL, Wu C, Lei C, Baraoidan M, Bordeos A, et al. (2005) Chemical- and irradiation-induced mutants of indica rice IR64 for forward and reverse genetics. Plant Mol Biol 59: 85–97.
  49. 49. Chawade A, Sikora P, Brautigam M, Larsson M, Vivekanand V, et al. (2010) Development and characterization of an oat TILLING-population and identification of mutations in lignin and beta-glucan biosynthesis genes. BMC Plant Biol 10: 86.
  50. 50. Gottwald S, Bauer P, Komatsuda T, Lundqvist U, Stein N (2009) TILLING in the two-rowed barley cultivar ‘Barke’ reveals preferred sites of functional diversity in the gene HvHox1. BMC Res Notes 2: 258.
  51. 51. Wang N, Wang Y, Tian F, King GJ, Zhang C, et al. (2008) A functional genomics resource for Brassica napus: development of an EMS mutagenized population and discovery of FAE1 point mutations by TILLING. New Phytol 180: 751–765.
  52. 52. Till BJ, Cooper J, Tai TH, Colowit P, Greene EA, et al. (2007) Discovery of chemically induced mutations in rice by TILLING. BMC Plant Biol 7: 19.
  53. 53. Caldwell DG, McCallum N, Shaw P, Muehlbauer GJ, Marshall DF, et al. (2004) A structured mutant population for forward and reverse genetics in Barley (Hordeum vulgare L.). Plant J 40: 143–150.
  54. 54. Silva JA, Da Costa TS, Lucchetta L, Marini LJ, Zanuzo MR, et al. (2004) Characterization of ripening behavior in transgenic melons expressing an antisense 1-aminocyclopropane-1-carboxylate (ACC) oxidase gene from apple. Postharvest Biology and Technology 32: 263–268.
  55. 55. Nunez-Palenius HG, Cantliffe DJ, Huber DJ, Ciardi J, Klee HJ (2006) Transformation of a muskmelon ‘Galia’ hybrid parental line (Cucumis melo L. var. reticulatus Ser.) with an antisense ACC oxidase gene. Plant Cell Rep 25: 198–205.
  56. 56. Hamilton AJ, Lycett GW, Grierson D (1990) Antisense gene that inhibits synthesis of the hormone ethylene in transgenic plants. Nature 346: 284–287.
  57. 57. Murray AJ, Hobson GE, Schuch W, Bird CR (1993) Reduced ethylene synthesis in EFE antisense tomatoes has differential effects on fruit ripening processes. Postharvest Biology and Technology 2: 301–313.
  58. 58. Xiong AS, Yao QH, Peng RH, Li X, Han PL, et al. (2005) Different effects on ACC oxidase gene silencing triggered by RNA interference in transgenic tomato. Plant Cell Rep 23: 639–646.
  59. 59. Ng , Henikoff (2003) SIFT: Predicting amino acid changes that affect protein function. Nucleic acids research 31: 3812–3814.
  60. 60. Zhang Z, Ren JS, Clifton IJ, Schofield CJ (2004) Crystal structure and mechanistic implications of 1-aminocyclopropane-1-carboxylic acid oxidase–the ethylene-forming enzyme. Chem Biol 11: 1383–1394.
  61. 61. Zhang Z, Ren J, Stammers DK, Baldwin JE, Harlos K, et al. (2000) Structural origins of the selectivity of the trifunctional oxygenase clavaminic acid synthase. Nat Struct Biol 7: 127–133.
  62. 62. Roach PL, Clifton IJ, Fulop V, Harlos K, Barton GJ, et al. (1995) Crystal structure of isopenicillin N synthase is the first from a new structural family of enzymes. Nature 375: 700–704.
  63. 63. Roach PL, Clifton IJ, Hensgens CM, Shibata N, Schofield CJ, et al. (1997) Structure of isopenicillin N synthase complexed with substrate and the mechanism of penicillin formation. Nature 387: 827–830.
  64. 64. Valegard K, van Scheltinga AC, Lloyd MD, Hara T, Ramaswamy S, et al. (1998) Structure of a cephalosporin synthase. Nature 394: 805–809.
  65. 65. Elkins JM, Ryle MJ, Clifton IJ, Dunning Hotopp JC, Lloyd JS, et al. (2002) X-ray crystal structure of Escherichia coli taurine/alpha-ketoglutarate dioxygenase complexed to ferrous iron and substrates. Biochemistry 41: 5185–5192.
  66. 66. Slade AJ, Fuerstenberg SI, Loeffler D, Steine MN, Facciotti D (2005) A reverse genetic, nontransgenic approach to wheat crop improvement by TILLING. Nat Biotechnol 23: 75–81.
  67. 67. Boualem A, Troadec C, Kovalski I, Sari MA, Perl-Treves R, et al. (2009) A conserved ethylene biosynthesis enzyme leads to andromonoecy in two Cucumis species. PLoS One 4: e6144.
  68. 68. Voss HD (1992) Relating colorimeter measurement of plant color to the Royal Horticultural Society Colour Chart. HortScience 27: 1256–1260.