MELOGEN: an EST database for melon functional genomics
© Gonzalez-Ibeas et al; licensee BioMed Central Ltd. 2007
Received: 08 May 2007
Accepted: 03 September 2007
Published: 03 September 2007
Melon (Cucumis melo L.) is one of the most important fleshy fruits for fresh consumption. Despite this, few genomic resources exist for this species. To facilitate the discovery of genes involved in essential traits, such as fruit development, fruit maturation and disease resistance, and to speed up the process of breeding new and better adapted melon varieties, we have produced a large collection of expressed sequence tags (ESTs) from eight normalized cDNA libraries from different tissues in different physiological conditions.
We determined over 30,000 ESTs that were clustered into 16,637 non-redundant sequences or unigenes, comprising 6,023 tentative consensus sequences (contigs) and 10,614 unclustered sequences (singletons). Many potential molecular markers were identified in the melon dataset: 1,052 potential simple sequence repeats (SSRs) and 356 single nucleotide polymorphisms (SNPs) were found. Sixty-nine percent of the melon unigenes showed a significant similarity with proteins in databases. Functional classification of the unigenes was carried out following the Gene Ontology scheme. In total, 9,402 unigenes were mapped to one or more ontology. Remarkably, the distributions of melon and Arabidopsis unigenes followed similar tendencies, suggesting that the melon dataset is representative of the whole melon transcriptome. Bioinformatic analyses primarily focused on potential precursors of melon micro RNAs (miRNAs) in the melon dataset, but many other genes potentially controlling disease resistance and fruit quality traits were also identified. Patterns of transcript accumulation were characterised by Real-Time-qPCR for 20 of these genes.
The collection of ESTs characterised here represents a substantial increase on the genetic information available for melon. A database (MELOGEN) which contains all EST sequences, contig images and several tools for analysis and data mining has been created. This set of sequences constitutes also the basis for an oligo-based microarray for melon that is being used in experiments to further analyse the melon transcriptome.
Melon (Cucumis melo L.) is an important horticultural crop grown in temperate, subtropical and tropical regions worldwide. Melon is among the most important fleshy fruits for fresh consumption, its total production in 2004 exceeding 874 million metric tons, of which 72.5% are produced in Asia, 11.7% in Europe, 8.4% in America and 6.1% in Africa, being a significant component of fresh fruit traded internationally . Melon belongs to the Cucurbitaceae family, which comprises up to 750 different species distributed in 90 genera. Species in this family include watermelon, cucumber, squash and marrow, all of them cultivated essentially because of their fruits, but this family also includes species of interest for other reasons, as, for example, their contents in potentially therapeutic compounds (e.g. Momordica charantia) . Melon is a diploid species, with a basic number of chromosomes x = 12 (2x = 2n = 24) and an estimated genome size of 450 to 500 Mb , similar in size to the rice genome (419 Mb) [4, 5] and about three times the size of the Arabidopsis genome (125 Mb) . Melon has been classified into two subspecies, C. melo ssp. agrestis and C. melo ssp. melo with India and Africa being their centres of origin, respectively [7, 8].
Melon has a great potential for becoming a model for understanding important traits in fruiting crops. Melon fruits have wide morphological, physiological and biochemical diversity [7, 9] which can be exploited to dissect biological processes of great technological importance, among them flavour development and textural changes that occur during fruit ripening. The contemporary melon cultivars can be divided into two groups, climacteric and nonclimacteric, according to their ripening patterns . Climacteric fruits are characterized by rapid and profound changes during ripening associated to increased levels of respiration and release of ethylene, whereas the nonclimacteric varieties do not produce ethylene and have long shelf-life. Analyses of climacteric and nonclimacteric melons have illustrated the process of aroma formation [11–14] and the temporal sequence of cell wall disassembly [15–17]. Melon can be also a very useful experimental system to analyse other aspects of fundamental plant biology. For example, melon and other cucurbits have been used to analyse the development of the plant vasculature and the transportation of macromolecules through it [18–20], and different interactions between melon and pests and pathogens have been characterised with varying depths [21–27].
Important genetic tools have been described for melon, as for example linkage genetic maps [28, 29] and the development of a genomic library of near isogenic lines (NILs) from an exotic accession ; also, biotechnology is feasible in melon [31–33]. However, the great majority of genes involved in the aforementioned traits are yet to be identified in melon. Partial sequencing of cDNA inserts of expressed sequence tags (ESTs) have been used as an effective method for gene discovery. By sequencing clones derived from RNA from different sources, and/or by normalizing cDNA libraries, the total set of genes sampled can be maximized. Bioinformatic analysis, annotation and clustering of sequences could yield databases which mining can be used to select candidate genes implicated in traits of interest. EST collections can also serve to construct microarrays useful for identifying sets of plant genes expressed during different developmental stages and/or responding to environmental stimuli [34, 35]. In addition, EST collections are good sources of simple sequence repeats (SSRs) and single-nucleotide polymorphisms (SNPs) that can be used for creating saturated genetic maps [36, 37]. Thus, EST collections have been generated for many plant species, being the most comprehensive those of Arabidopsis and rice . Fruit crops have been less extensively surveyed, but important collections are publicly available for several species, including tomato , apple , grape  and citrus .
Despite the importance of the family Cucurbitaceae, relatively little EST information is currently available: only 16,039 nucleotide sequences have been annotated from the whole Cucurbitaceae family in the publicly accessible GenBank database as of November 2006; out of these, 12,180 correspond to the Cucumis genus and 6,061 to melon. These numbers are in sharp contrast with the data available for families composed of other important food crops like Solanaceae (1,020,102 sequences), Fabaceae (1,466,518 sequences), Brassicaceae (1,010,148 sequences excluding Arabidopsis), Vitaceae (449,478) and Rosaceae (390,066 sequences). Here we describe a public EST sequencing project in melon. We report the determination and analysis of 30,675 high-quality melon ESTs, sequenced from eight normalized cDNA libraries corresponding to different tissues in different physiological conditions. We have classified the sequences into functional categories and described SSRs and SNPs of potential use in genetic maps and marker-assisted breeding programs. A database which contains all EST sequences, contig images and several tools for analysis and data mining has been created. In addition, we have analyzed the EST melon dataset to identify candidate genes potentially coding microRNAs or involved in fruit maturation processes and pathogen defence. The pattern of transcript accumulation in different physiological conditions has been characterised by Real-Time-qPCR for 20 of these candidate genes.
EST Sequencing and Clustering
Description of cDNA libraries
Ssp. melo cv. "Piel de Sapo" T-111
Fruit 15 days after pollination
Ssp. melo cv. "Piel de Sapo" T-111
Fruit 46 days after pollination
Ssp. agrestis accession pat81
Ssp. agrestis accession pat81
Roots infected with M. cannonballus
Ssp. melo var. cantaloupe accession C-35
Photosynthetic cotyledons infected with CMV
Ssp. melo var. cantaloupe accession C-35
Ssp. melo cv. "Piel de Sapo" Piñonet torpedo
Ssp. melo cv. "Piel de Sapo" Piñonet torpedo
Roots infected with M. cannonballus
608.1 ± 175.2
583.0 ± 161.1
700.0 ± 185.4
756.3 ± 137.1
651.4 ± 205.7
669.3 ± 171.1
679.9 ± 198.7
749.3 ± 156.2
ESTs showing significant similarity with fungal sequences
Number of ESTs
SSRs and SNPs
Simple sequence repeats (SSRs) statistics*
Number of di-pSSR
Number of tri-pSSR
Number of tetra-pSSR
Localization of simple sequence repeats (SSRs) with respect to putative initiation and termination codons in the melon sequence dataset*
All SSRs analyzed
Single nucleotide polymorphisms (SNPs) statistics*
Number of pSNP
Number of pSCH
Functional annotation statistics
A. Number of unigenes with BLAST hits
Number of unigenes
B. Number of unigenes with HMMER hits
Number of unigenes
C. Number of unigenes with orthologue
Number of unigenes
We have also identified 6,673 (40.1%) melon unigenes with an ortholog in the Arabidopsis database, and a HMMER motif has been assigned to 4,655 (28.0%) unigenes by comparisons with the Pfam database  (Table 7). All these results are compiled in the MELOGEN database, which also contains direct links to the databases used to carry out analyses.
Genes potentially encoding microRNAs
Potential melon miRNAs
Potential mature miRNA sequence (5'->3')*
Precursor folding free energy (kcal mol-1)§
MFEI (kcal mol-1)Δ
Potential target in melon
Genes potentially encoding pathogen resistance and fruit quality traits
Pathogens affect severely the productivity of melon crops. Three of the cDNA libraries sequenced here correspond to pathogen-infected tissues and, thus, should contain transcripts from genes whose expression is induced in response infection. We have carried out a bioinformatics search for homologs of genes involved in pathogen resistance response (see  for a review) and virus susceptibility [59–61], finding among them at least one melon ortholog to the Arabidopsis FLS2 receptor , several unigenes potentially encoding disease resistance proteins as well as mitogen-activated protein kinases, homologs to translation initiation factors constituting potential virus susceptibility factors, etc. [see Additional file 1].
Fruit development and ripening are the most important processes determining the fruit quality traits of fleshy fruits like melon. At present most of the molecular and genetic data available about fruit development and ripening come from tomato [63, 64] and Arabidopsis [65, 66]. In recent years, several genes and quantitative trait loci controlling fruit quality traits have been described in melon [67, 68]. As for developmental processes, homologs to genes involved in melon fruit development, ripening and quality have been found in the melon dataset. These include several MADS-box genes, homologs to the fw2.2 and ovate QTLs [69, 70], several homologs to members of the SBP-box gene family to which the major tomato ripening gene COLORLESS NON-RIPENING belongs , several ACC synthase and ACC oxidase genes, unigenes from several cell wall-metabolism enzymes, etc [see Additional file 1].
Expression analysis of selected ESTs by Real-Time-qPCR
Transcripts selected from the database for gene expression analysis by Real Time qPCR
Sequence length (bp)
Aminoacid similarity (%)
Annotation (HMMR domain)
Chitinase-like protein 1, similar to class I chitinase (Glyco_hydro_19)
Peptidyl-prolyl cis-trans isomerase, cyclophilin
Eukaryotic translation initiation factor 4A-2 (DEAD)
Eukaryotic translation initiation factor 4E 1 (IF4E)
Ethylene receptor (Response_reg)
Gibberellin 2-oxidase (2OG-FeII_Oxy)
Heat shock protein 101 (AAA_2)
Heat shock protein 70/HSC70-5 (HSP70)
Auxin-responsive protein/indoleacetic acid- induced protein 9 (AUX_IAA)
Small nuclear ribonucleoprotein/snRNP, putative/Sm protein, putative, similar to U6 snRNA-associated Sm-like protein
Lycopene epsilon cyclase (Lycopene_cycl)
Novel cap-binding protein (IF4E)
Short vegetative phase protein, MADS box transcription factor related cluster (SRF-TF)
Band 7 family protein, strong similarity to hypersensitive-induced response protein (Band_7)
Xyloglucan:xyloglucosyl transferase/xyloglucan endotransglycosylase/endo- xyloglucan transferase/TCH4
Tonoplast intrinsic protein (MIP)
Transmembrane protein-related/TOM1 (DUF1084)
Senescence-associated family protein
Tobamovirus multiplication protein 3/THH1 (DUF1084)
UDP-glucose 4-epimerase/UDP-galactose 4- epimerase/Galactowaldenase (Epimerase)
WRKY family transcription factor, DNA-binding protein (WRKY)
The response of selected genes in roots inoculated with M. cannonballus was analysed in melon genotypes known to be susceptible (cultivar "Piel de sapo"; Fig. 4B) and partially resistant (accession pat81 of C. melo L. ssp. agrestis; Fig. 4C) to the infection by this fungus. The patterns of transcript accumulation resulted clearly different for both genotypes. For pat81 (resistant), transcription factors WRKY70 and SVP increased their expression between 2 and 3 times after inoculation; other stress-inducible genes (HSP101, HSP70) showed only a moderate increase (Fig. 4C). For "Piel de sapo" (susceptible), accumulation of WRKY70 and SVP transcripts only increased about 1.5 times after inoculation whereas the expression of HSP101 showed a marked increase (Fig. 4B). It is also worth noting the differential response of the GA2OX1 gene in the two genotypes. Expression of GA2OX1 increased about 1.5 times in pat81 roots after the M. cannonballus attack, whereas it decreased in "Piel de sapo" roots after fungal infection (compare Figs. 4B and 4C).
Comparison of patterns of transcript accumulation at two stages of fruit development showed increased levels of gene expression for 9 of the analysed genes. This was particularly evident for HSP70, TOM2A, TOM3, EIN4 and IAA9. In contrast, decreased levels of transcript accumulation were observed for the other 11 genes.
In this paper we provide an initial platform for functional genomics of melon by the identification of more than 16,000 unigenes assembled from almost 30,000 ESTs sequenced from 8 melon cDNA libraries. It is probably premature to estimate the proportion of melon genes represented in this dataset, but based on available data for other plant species (i.e. Arabidopsis and rice), it is likely that the melon unigene set characterised here represents approximately between half and one-third of the number of expressed, protein coding genes of melon. Libraries were constructed from various tissue types, but with a bias towards fruit development and pathogen-infected tissues. Data from these libraries will become a useful resource of genes for experiments aimed at understanding important processes involved in fruit development and resistance to viral and fungal pathogens. Also, data presented here provide an important tool for generating markers to saturate melon genetic maps.
In contrast to typical EST gene-sampling strategies reported previously, we have found a low degree of redundancy in the sequences determined. The process of clustering reduced the number of sequences to 56%, from 29,604 good quality ESTs to 6,023 contigs and 10,614 singletons. Contigs with more than 8 ESTs were scarce, the majority of them being formed by 3 or 2 ESTs. Redundancy of the sequences derived from each library ranged from 13% to 20%, with singletons constituting approximately one third of the unigenes determined per library. This low redundancy is probably due to the success of the normalization process, responsible for the suppression of superabundant transcripts specific for a given tissue or condition. Normalization precludes in silico analysis of gene expression, but greatly increases the number of unigenes that can be determined by reducing redundancy . Here we have used a recently described normalization protocol which is based on the cleavage of DNA or DNA-RNA duplexes by a specific DNase ; this process, in our hands, has proven simple, reproducible and efficient. Another factor that has contributed to the low redundancy values obtained has been the sequencing of libraries from very distinct tissues. Thus, the number of library specific unigenes was about one half of the total number of unigenes contributed by each library, suggesting that further sequencing of the libraries still has the potential to provide a good number of new, non-redundant sequences.
cDNA sequences are a useful source of SSRs, which are excellent molecular markers due to their high degree of polymorphism. A common feature of cDNA sequences obtained from plants is the high frequency of SSRs that they contain . We have identified more than 1,000 potential SSRs in the melon dataset, with approximately 6% of the melon unigenes containing di-, tri- or tetranucleotide repeats. A clear bias toward AG and AAG repeats existed, that account for 67% of the SSRs. In contrast, the GC repeat was not found in the melon dataset. A similar bias toward AG and against CG repeats has been identified in Arabidopsis and other plant species [40, 78]. As proposed at least in one other instance , this may be due to the tendency of CpG sequences to be methylated , which potentially might inhibit transcription. Another interesting feature of melon SSRs relates to their pattern of localization with respect to putative initiation and termination codons. It is known that the UTRs of transcribed sequences are richer in SSRs than coding regions, particularly at the 5'-UTRs [36, 40]. However, in the melon dataset, a high proportion of SSRs (29.5%) were found in ORFs. An analysis of the localization of di-, tri- and tetranucleotide repeats separately showed that di- and tetranucleotides were preferentially located in UTRs, whereas trinucleotides localised in both, UTRs and ORFs, consistently with maintenance of the ORFs coding capacity. Thus, the prevalence of trinucleotide repeats in the melon dataset (71%) explains this result.
We identified in the melon sequence dataset 356 high-quality SNPs. Since non-redundant sequences analysed here encompassed 4.5 Mb, one SNP was found every 12,000 pb of sequence. This small figure is probably due to the limited number of melon genotypes used and the low redundancy found among libraries. In fact, when the frequency of SNPs is computed in relation to the length and number of contigs containing SNPs, the corresponding value (one SNP in every 616 bp of sequence) is of the same order of magnitude as values previously calculated for melon (441 bp; ) and other plant species . With the advent of high-throughput detection systems, the SSRs and SNPs identified here will constitute an important resource for mapping and marker-assisted breeding in melon and closely related crops.
As an approach to the function of melon unigenes, we carried out a bioinformatics analysis based on BLASTX and matches with the Pfam database . The proportion of melon unigenes with no similar sequences in databases was quite high, suggesting that the melon dataset may encompass an important number of melon-specific sequences. However, the proportion specific sequences might be overestimated because blasting has been made with unigene sequences, which in many cases do not cover the complete length of the transcript. We performed a functional classification of the unigenes following the Gene Ontology scheme, which is one of the more versatile and complete systems for functional classification . A comparison of the distributions of melon and Arabidopsis unigenes in GO categories showed that both followed similar tendencies, suggesting that the melon dataset is representative of the whole melon transcriptome. This is remarkable, as the number of different libraries sequenced has been relatively small; again, this is probably due to the success of the normalization process. We have also carried out specific searches for genes involved in pathways of particular relevance in melon, such are resistance response and fruit development, identifying a remarkable number of melon candidates. For example, an ortholog of the flagellin receptor FLS2 from Arabidopsis  has been identified, together with 163 candidate RLKs that may have critical roles in pathogen recognition or diverse signalling processes. Similarly, up to 8 MADS-box gene homologs with potential roles in development have been found in the melon dataset. Moreover, a bioinformatics approach [52, 53] allowed the identification of potential precursors of melon miRNAs together with several potential targets in the melon dataset. This finding opens the door to biotechnology approaches based on the use of artificial miRNAs to specifically silence melon genes [81, 82].
The transcript accumulation analysis for the 20 selected genes revealed important changes in gene expression associated with pathogen infection and fruit development. For virus infection, the accumulation of transcripts remained unaltered for 12 genes, but showed a significant increase for 6 genes and a decrease for 2. Among the set of genes analysed, TOM2A and EIF4E were known to code for virus susceptibility factors [27, 83]; the expression of TOM2A was increased after CMV infection, consistently with its requirement by the virus, but this was not the case for EIF4E. Different hypotheses can explain this result: since EIF4E is an abundant, housekeeping protein, increased expression may not be essential for virus multiplication; alternatively, CMV may not use EIF4E in melon or may use a factor coded by a different member of the 4E family; it may also be that timing of the sampling for this experiment was not appropriate to detect such an effect, as requirement of EIF4E might occur very early during virus multiplication. In the cases of infection of susceptible ("Piel de sapo") and resistant (pat81) melons by M. cannonballus, more extensive alterations in gene expression seemed to occur in the susceptible than in the resistant accession. Significantly, for the susceptible accession, stress responsive genes (e.g. HSP101) appeared to be maximally induced, whereas for the resistant accession, a gene encoding a WRKY70 transcription factor, potentially involved in resistance response, was induced to high levels. Significantly, expression of GA2OX1 increased about 1.5 times in pat81 after the M. cannonballus attack, whereas it decreased in "Piel de sapo". GA2ox is a major gibberellin (GA) catabolic enzyme, with an important role in controlling GA levels in plants. Hormones control many plant developmental processes, and strong evidence indicates that hormone signalling is involved in the regulation of root growth and architecture [84, 85]. The differential response of the GA2OX1 gene in the two melon genotypes is consistent with an enhanced root growth in pat81 after infection . Notably, other genes involved in hormone-mediated signalling pathways, such as the IAA9 gene, did not show such differential response to M. cannonballus infection in both genotypes. In the case of fruit development, differences in the expression of selected genes between immature and ripening fruits appeared to be even sharper than in the cases of healthy and pathogen-infected tissues. Specific roles during fruit development for HSP70, TOM2A and TOM3 have not been identified, though an increased expression has been shown at least in the case of TOM2A in tomato . The ethylene receptor gene EIN4 showed a two-fold increase in expression. EIN4 is the ortholog of Arabidopsis EIN4 and tomato LeETR4[88, 89]. In tomato, LeETR4 is also highly expressed in ripening fruit, suggesting that it responds by modulating ethylene signalling during ripening . The MADS-box gene (SVP) showed about a four-fold decrease in expression. This gene is the ortholog of tomato JOINTLESS, which specifies the abscission zone in tomato. In tomato fruit microarray hybridizations, the expression of JOINTLESS also decreased from 7 to 57 DAP , in agreement with our data for melon. The lycopene epsilon cyclase (LUT2) and xyloglucan endotransglycosylase (TCH4) genes showed an approximately four-fold decrease in expression during melon fruit development. These findings fit with the patterns of expression of these genes in tomato, where their transcript levels decrease to a non-detectable level in the ripe fruits [90, 91].
In summary, this collection of ESTs represents a substantial increase on the information available for melon. The dataset contains SSR and SNP markers that can be used for breeding, as well as a significant number of candidate genes that can be experimentally tested for their roles in various important processes. This set of genes constitutes also the basis for a microarray for melon that is being used in experiments to further analyse fruit development and maturation and responses to pathogen infections.
The cDNA libraries were prepared using material from four different melon genotypes: the line T-111 (Semillas Fitó, Barcelona, Spain), which corresponds to a Piel de Sapo breeding line, the Piel de Sapo cultivar "Piñonet torpedo" (Semillas Batlle, Barcelona, Spain), the accession C-35 of the germplasm collection of La Mayora-CSIC (EELM-CSIC, Málaga, Spain), which corresponds to a cantaloupe-type of melon, and the accession pat81 of C. melo L. ssp. agrestis (Naud.) Pangalo maintained at the germplasm bank of COMAV (COMAV-UPV, Valencia, Spain) (Table 1). Seeds of line T-111 were germinated at 30°C for two days and plants were grown in a greenhouse in peat bags, drip irrigated, with 0.25-m spacing between plants. Fruits of 15 and 46 days after pollination were collected and mesocarp tissues were recovered and used for RNA extractions. Root samples were from Piel de sapo and pat81 plants, both healthy and inoculated with M. cannonballus. Piel de sapo is fully susceptible to the infection by this fungus whereas pat81 has been shown to be partially resistant [92, 93]. Seeds were pre-germinated in Petri dishes. After 4 days, seedlings were transplanted to 0.5-l pots filled with sterile soil substrate and grown in a greenhouse (20–35°C, 60–85% relative humidity). Inoculations were carried out by adding 50 colony-forming units (CFU) of M. cannonballus per gram of sterile soil as described by Iglesias et al. . Fourteen days after inoculation, healthy and inoculated roots were collected for RNA extraction. The presence of the fungus and the infection levels were assessed by real-time quantitative PCR as described by Picó et al. . CMV infected cotyledons were collected from plants of the C-35 accession. In this case, seeds were pregerminated in Petri dishes for 24 h at 28°C in the dark, planted in 0.5-l pots and maintained in an insect-proof green house (20–28°C, 45 to 85% relative humidity) for 6 to 7 days, until the first true leaf started emerging. At this stage, cotyledons were mechanically inoculated with CMV following standard procedures . Inoculated cotyledons were harvested 4 days after inoculation and used for RNA extractions. Dot-blot hybridisation  was used to check infection by CMV. Plants of the C-35 accession were also used for collecting healthy leaves. Plants were maintained in the greenhouse for 21 days, and second and third leaves above cotyledons were harvested for RNA extractions.
Construction of cDNA libraries and EST sequencing
Total RNA was prepared as described by Aranda et al. . Poly(A+) RNA from total RNA was purified using MicroPoly(A+) Purist (Ambion, Austin, TX, USA), a cellulose-oligo(dT)-based method. Integrity and quality of both total and poly(A+) RNA were tested by gel electrophoresis. cDNA libraries were constructed with the SMART cDNA Library Construction kit (Clontech, Mountain View, CA, USA), using a modified primer to include a Sfi I enzyme restriction site. A normalization step was carried out with TRIMMER kit (Evrogen, Moscow, Russia). After normalization, a cDNA fractionation step was performed with SizeSep 400 Spun Columns (Amersham Biosciences, Buckinghamshire, England). cDNA was digested with Sfi I, generating Sfi IA-Sfi IB cohesive ends for directional cloning into a modified version of BlueScript SK plasmid vector (Stratagene, La Jolla, CA, USA). Ligation products were transformed into E. coli electrocompetent cells DH10B (Invitrogen, Carlsbad, CA, USA) by electroporation. The titer of the libraries was evaluated by plating an aliquot on LB agar plates with ampicillin at 100 μg ml-1. Only libraries of 105 cfu ml-1 or more were considered as acceptable. Prior to large scale sequencing, the average insert size was estimated by restriction analyses of 24 plasmid DNA minipreps per library from randomly picked colonies.
Sequencing was carried out from the 5'-end of the inserts without library amplification using the universal M13 reverse primer. An external custom service was contracted for this task (Macrogen Inc., Seoul, Korea). Approximately 6,000 clones were sequenced from the CI library, and 3,500 clones were sequenced from each of the other libraries (Table 2). Sequences obtained in this work can be found in GenBank [accession numbers AM713476 to AM743079] and MELOGEN .
EST sequences were automatically trimmed, clustered and annotated using the EST2uni analysis pipeline . EST2uni compromises the analysis pipeline written in PERL , a database (MySQL)  and a web site to browse the results coded in PHP . Thus, for the EST pre-processing step, base calling was performed with Phred , low quality regions and vector sequences were trimmed with Lucy , and repeats and low complexity regions were masked with RepeatMasker  and Seqclean . Further vector contamination was also eliminated with Seqclean using NCBI's UniVec . High-quality EST sequences were then assembled to obtain the unigene set using Tgicl .
Detection of SSRs was performed using Sputnik . Putative SNPs were annotated when the least represented allele was present in two EST sequences or more. ORFs were predicted in the ESTs with the aid of the ESTScan software .
For functional annotation, comparisons against the Arabidopsis (TAIR)  and Uniref  databases were carried out using BLASTN or BLASTX for nucleotide or protein sequences, respectively. Functional domains were searched with HMMPFAM  using the Pfam database . The Gene Ontology (GO) classification  was derived from the BLASTN results against the Arabidopsis proteome. Also, a bi-directional BLASTN comparison was performed in order to obtain a set of putative orthologs with Arabidopsis. Finally, a set of superunigenes was obtained grouping different unigenes with the same expected mRNA target, as judged by extensive sequence overlapping.
To assess codon usage, we generated a set of melon sequences predicted to contain full-length coding regions. These sequences were subjected to BLASTX and, after manual inspection, sequences showing a high similarity to Arabidopsis proteins were selected to ensure that no sequences containing frame-shift errors were included in the analysis. From this smaller dataset, which included 588 sequences, ORFs were defined and a codon usage table was created. Codon usage was calculated from sequences using the GCUA program . All codons were found in the dataset, with the least frequent codon represented 134 times.
To identify potential melon miRNAs, the 33,292 melon ESTs were subjected to a BLAST search against mature sequences of known miRNAs from the miRNA Registry Database (released January 2007)  using BLASTN . ESTs with only 0–1 mismatched nucleotides with known miRNAs were considered. Selected ESTs were subjected to a BLAST search against protein databases in order to remove potential protein-coding sequences. ESTs pertaining to the same melon unigene of the MELOGEN database were grouped. The secondary structures of the unigenes encoding potential miRNA precursors were predicted with the web-based tool MFOLD , using default parameters. In each case, only the lower energy structure was selected for visual inspection, as previously described [52, 114]. In order to select unigenes with perfect or near-perfect secondary foldback hairpin structures, only sequences with a maximum size of 3 nucleotides for a bulge in the miRNA sequence and with at least 16 paired nucleotides between the mature sequence and the opposite arm were considered as potential miRNA candidates. In addition, the minimal folding free energy index (MFEI) for each sequence was calculated following Zhang et al. .
Gene expression analyses
Primer sequences for Real Time-qPCR analysis of transcript accumulation
Primer sequence (5'->3') Forward
This work was supported by grants from Ministerio de Educación y Ciencia (Spain) (GEN2003-20237-C06) and Consejería de Educación y Cultura (Región de Murcia, Spain) (BIO2005/04-6436). Wim Deleu, Cristina Roig and Daniel Gonzalez-Ibeas are recipient of a postdoctoral fellowship from the Centre de Recerca en Agrigenòmica CSIC-IRTA (Spain), a Juan de la Cierva grant from Ministerio de Educación y Ciencia (Spain) and a predoctoral fellowship from Ministerio de Educación y Ciencia (Spain), respectively.
- FAOSTAT Agriculture data. [http://faostat.fao.org/default.aspx]
- Jayasooriya AP, Sakono M, Yukizaki C, Kawano M, Yamamoto K, Fukuda N: Effects of Momordica charantia powder on serum glucose levels and various lipid parameters in rats fed with cholesterol-free and cholesterol-enriched diets. J Ethnopharmacol. 2000, 72: 331-336. 10.1016/S0378-8741(00)00259-2.PubMedView ArticleGoogle Scholar
- Arumuganathan K, Earle ED: Nuclear DNA content of some important plant species. Plant Mol Biol Rep. 1991, 9: 208-218.View ArticleGoogle Scholar
- Goff SA, Ricke D, Lan TH, Presting G, Wang R, Dunn M, Glazebrook J, Sessions A, Oeller P, Varma H, others: A Draft Sequence of the Rice Genome (Oryza sativa L. ssp. japonica). Science. 2002, 296: 92-100. 10.1126/science.1068275.PubMedView ArticleGoogle Scholar
- Yu J, Hu S, Wang J, Wong GK-S, Li S, Liu B, Deng Y, Dai L, Zhou Y, Zhang X, others: A Draft Sequence of the Rice Genome (Oryza sativa L. ssp. indica). Science. 2002, 296: 79-92. 10.1126/science.1068037.PubMedView ArticleGoogle Scholar
- Arabidopsis Genome Initiative: The Arabidopsis Information Resource (TAIR): a comprehensive database and web-based information retrieval, analysis, and visualization system for a model plant. Nucleic Acids Res. 2001, 29: 102-105. 10.1093/nar/29.1.102.View ArticleGoogle Scholar
- Kirkbride JH: Biosystematic monograph of the genus Cucumis (Cucurbitaceae). 1993, Boone, North Carolina: Parkway PublishersGoogle Scholar
- Garcia-Mas J, Monforte AJ, Arus P: Phylogenetic relationships among Cucumis species based on the ribosomal internal transcribed spacer sequence and microsatellite markers. Plant Syst Evol. 2004, 248: 191-203. 10.1007/s00606-004-0170-y.View ArticleGoogle Scholar
- Liu L, Kakihara F, Kato M: Characterization of six varieties of Cucumis melo L. based on morphological and physiological characters, including shelf-life of fruit. Euphytica. 2004, 135: 305-313. 10.1023/B:EUPH.0000013330.66819.6f.View ArticleGoogle Scholar
- Miccolis V, Saltveit ME: Morphological and physiological changes during fruit growth and maturation of seven melon cultivars. J Am Soc Hort Sci. 1991, 116: 1025-1029.Google Scholar
- Shalit M, Katzir N, Tadmor Y, Larkov O, Burger Y, Shalekhet F, Lastochkin E, Ravid U, Amar O, Edelstein M, others: Acetyl-CoA: Alcohol acetyltransferase activity and aroma formation in ripening melon fruits. J Agric Food Chem. 2001, 49: 794-799. 10.1021/jf001075p.PubMedView ArticleGoogle Scholar
- Bauchot AD, Mottram DS, Dodson AT, John P: Effect of Aminocyclopropane-1-carboxylic Acid Oxidase Antisense Gene on the Formation of Volatile Esters in Cantaloupe Charentais Melon (Cv. Védrandais). J Agric Food Chem. 1998, 46: 4787-4792. 10.1021/jf980692z.View ArticleGoogle Scholar
- Flores F, El Yahyaoui F, de Billerbeck G, Romojaro F, Latche A, Bouzayen M, Pech JC, Ambid C: Role of ethylene in the biosynthetic pathway of aliphatic ester aroma volatiles in Charentais Cantaloupe melons. J Exp Bot. 2002, 53: 201-206. 10.1093/jexbot/53.367.201.PubMedView ArticleGoogle Scholar
- Yahyaoui FEL, Wongs-Aree C, Latche A, Hackett R, Grierson D, Pech JC: Molecular and biochemical characteristics of a gene encoding an alcohol acyl-transferase involved in the generation of aroma volatile esters during melon ripening. FEBS J. 2002, 269: 2359-2366.View ArticleGoogle Scholar
- Hadfield KA, Bennett AB: Polygalacturonases: Many Genes in Search of a Function. Plant Physiol. 1998, 117: 337-343. 10.1104/pp.117.2.337.PubMed CentralPubMedView ArticleGoogle Scholar
- Rose JKC, Hadfield KA, Labavitch JM, Bennett AB: Temporal sequence of cell wall disassembly in rapidly ripening melon fruit. Plant Physiol. 1998, 117: 345-361. 10.1104/pp.117.2.345.PubMed CentralPubMedView ArticleGoogle Scholar
- Bennett AB: Biochemical and genetic determinants of cell wall disassembly in ripening fruit: A general model. Hortscience. 2002, 37: 447-450.Google Scholar
- Haritatos E, Keller F, Turgeon R: Raffinose oligosaccharide concentrations measured in individual cell and tissue types in Cucumis melo L leaves: Implications for phloem loading. Planta. 1996, 198: 614-622. 10.1007/BF00262649.View ArticleGoogle Scholar
- Volk GM, Turgeon R, Beebe DU: Secondary plasmodesmata formation in the minor-vein phloem of Cucumis melo L and Cucurbita pepo L. Planta. 1996, 199: 425-432. 10.1007/BF00195735.View ArticleGoogle Scholar
- Gomez G, Torres H, Pallas V: Identification of translocatable RNA-binding phloem proteins from melon, potential components of the long-distance RNA transport system. Plant J. 2005, 41: 107-116. 10.1111/j.1365-313X.2004.02278.x.PubMedView ArticleGoogle Scholar
- Chen JQ, Rahbé Y, Delobel B, Sauvion N, Guillaud J, Febvay G: resistance to the aphid Aphis gossypii: behavioural analysis and chemical correlations with nitrogenous compounds. Entomol Exp Appl. 1997, V85: 33-44. 10.1023/A:1003041228333.View ArticleGoogle Scholar
- Luo MZ, Wang YH, Frisch D, Joobeur T, Wing RA, Dean RA: Melon bacterial artificial chromosome (BAC) library construction using improved methods and identification of clones linked to the locus conferring resistance to melon Fusarium wilt (Fom-2). Genome. 2001, 44: 154-162. 10.1139/gen-44-2-154.PubMedView ArticleGoogle Scholar
- Klingler J, Powell G, Thompson GA, Isaacs R: Phloem specific aphid resistance in Cucumis melo line AR 5: effects on feeding behaviour and performance of Aphis gossypii. Entomol Exp Appl. 1998, 86: 79-88. 10.1023/A:1003153410049.View ArticleGoogle Scholar
- Marco CF, Aguilar JM, Abad J, Gomez-Guillamon ML, Aranda MA: Melon resistance to Cucurbit yellow stunting disorder virus is characterized by reduced virus accumulation. Phytopathology. 2003, 93: 844-852. 10.1094/PHYTO.2003.93.7.844.PubMedView ArticleGoogle Scholar
- Diaz JA, Nieto C, Moriones E, Truniger V, Aranda MA: Molecular characterization of a Melon necrotic spot virus strain that overcomes the resistance in melon and nonhost plants. Mol Plant-Microbe Interact. 2004, 17: 668-675. 10.1094/MPMI.2004.17.6.668.PubMedView ArticleGoogle Scholar
- Klingler J, Creasy R, Gao LL, Nair RM, Calix AS, Jacob HS, Edwards OR, Singh KB: Aphid resistance in Medicago truncatula involves antixenosis and phloem-specific, inducible antibiosis, and maps to a single locus flanked by NBS-LRR resistance gene analogs. Plant Physiol. 2005, 137: 1445-1455. 10.1104/pp.104.051243.PubMed CentralPubMedView ArticleGoogle Scholar
- Nieto C, Morales M, Orjeda G, Clepet C, Monfort A, Sturbois B, Puigdomenech P, Pitrat M, Caboche M, Dogimont C, others: An eIF4E allele confers resistance to an uncapped and non-polyadenylated RNA virus in melon. Plant J. 2006, 48: 452-462. 10.1111/j.1365-313X.2006.02885.x.PubMedView ArticleGoogle Scholar
- Perin C, Hagen L, Conto VD, Katzir N, Danin-Poleg Y, Portnoy V, Baudracco-Arnas S, Chadoeuf J, Dogimont C, Pitrat M: A reference map of Cucumis melo based on two recombinant inbred line populations. Theor Appl Genet. 2002, 104: 1017-1034. 10.1007/s00122-002-0864-x.PubMedView ArticleGoogle Scholar
- Gonzalo MJ, Oliver M, Garcia-Mas J, Monfort A, Dolcet-Sanjuan R, Katzir N, Arus P, Monforte A: Simple-sequence repeat markers used in merging linkage maps of melon (Cucumis melo L.). Theor Appl Genet. 2005, 110: 802-811. 10.1007/s00122-004-1814-6.PubMedView ArticleGoogle Scholar
- Eduardo I, Arus P, Monforte AJ: Development of a genomic library of near isogenic lines (NILs) in melon (Cucumis melo L.) from the exotic accession PI161375. Theor Appl Genet. 2005, 112: 139-148. 10.1007/s00122-005-0116-y.PubMedView ArticleGoogle Scholar
- Ayub R, Guis M, BenAmor M, Gillot L, Roustan JP, Latche A, Bouzayen M, Pech JC: Expression of ACC oxidase antisense gene inhibits ripening of cantaloupe melon fruits. Nat Biotechnol. 1996, 14: 862-866. 10.1038/nbt0796-862.PubMedView ArticleGoogle Scholar
- Guis M, Roustan JP, Dogimont C, Pitrat M, Pech JC: Melon biotechnology. Biotechnol Genet Engng Rev. 1998, 15: 289-311.View ArticleGoogle Scholar
- Gaba V, Zelcer A, Gal-On A: Cucurbit biotechnology – The importance of virus resistance. In Vitro Cell Dev Biol Plant. 2004, 40: 346-358. 10.1079/IVP2004554.View ArticleGoogle Scholar
- Rudd S: Expressed sequence tags: alternative or complement to whole genome sequences?. Trends Plant Sci. 2003, 8: 321-329. 10.1016/S1360-1385(03)00131-6.PubMedView ArticleGoogle Scholar
- Alba R, Fei Z, Payton P, Liu Y, Moore SL, Debbie P, Cohn J, D'Ascenzo M, Gordon JS, Rose JKC, others: ESTs, cDNA microarrays, and gene expression profiling: tools for dissecting plant physiology and development. Plant J. 2004, 39: 697-714. 10.1111/j.1365-313X.2004.02178.x.PubMedView ArticleGoogle Scholar
- Morgante M, Hanafey M, Powell W: Microsatellites are preferentially associated with nonrepetitive DNA in plant genomes. Nat Genet. 2002, 30: 194-200. 10.1038/ng822.PubMedView ArticleGoogle Scholar
- Rafalski JA: Novel genetic mapping tools in plants: SNPs and LD-based approaches. Plant Sci. 2002, 162: 329-333. 10.1016/S0168-9452(01)00587-8.View ArticleGoogle Scholar
- Ouyang S, Zhu W, Hamilton J, Lin H, Campbell M, Childs K, Thibaud-Nissen F, Malek RL, Lee Y, Zheng L, others: The TIGR Rice Genome Annotation Resource: improvements and new features. Nucleic Acids Res. 2007, 35: D883-D887. 10.1093/nar/gkl976.PubMed CentralPubMedView ArticleGoogle Scholar
- Fei ZJ, Tang XM, Alba R, Giovannoni J: Tomato Expression Database (TED): a suite of data presentation and analysis tools. Nucleic Acids Res. 2006, 34: D766-D770. 10.1093/nar/gkj110.PubMed CentralPubMedView ArticleGoogle Scholar
- Newcomb RD, Crowhurst RN, Gleave AP, Rikkerink EHA, Allan AC, Beuning LL, Bowen JH, Gera E, Jamieson KR, Janssen BJ, others: Analyses of expressed sequence tags from apple. Plant Physiol. 2006, 141: 147-166. 10.1104/pp.105.076208.PubMed CentralPubMedView ArticleGoogle Scholar
- Goes da Silva F, Iandolino A, Al Kayal F, Bohlmann MC, Cushman MA, Lim H, Ergul A, Figueroa R, Kabuloglu EK, Osborne C, others: Characterizing the Grape Transcriptome. Analysis of Expressed Sequence Tags from Multiple Vitis Species and Development of a Compendium of Gene Expression during Berry Development. Plant Physiol. 2005, 139: 574-597. 10.1104/pp.105.065748.View ArticleGoogle Scholar
- Forment J, Gadea J, Huerta L, Abizanda L, Agusti J, Alamar S, Alos E, Andres F, Arribas R, Beltran JP, others: Development of a citrus genome-wide EST collection and cDNA microarray as resources for genomic studies. Plant Mol Biol. 2005, 57: 375-391. 10.1007/s11103-004-7926-1.PubMedView ArticleGoogle Scholar
- EST2uni. [http://www.melogen.upv.es/genomica/web_estpipe/index.php]
- MELOGEN database. [http://www.melogen.upv.es]
- Sterky F, Bhalerao RR, Unneberg P, Segerman B, Nilsson P, Brunner AM, Charbonnel-Campaa L, Lindvall JJ, Tandre K, Strauss SH, others: A Populus EST resource for plant functional genomics. Proc Natl Acad Sci USA. 2004, 101: 13951-13956. 10.1073/pnas.0401641101.PubMed CentralPubMedView ArticleGoogle Scholar
- Hsiang T, Goodwin PH: Distinguishing plant and fungal sequences in ESTs from infected plant tissues. J Microbiol Methods. 2003, 54: 339-351. 10.1016/S0167-7012(03)00067-8.PubMedView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralPubMedView ArticleGoogle Scholar
- The Gene Ontology Consortium: Gene Ontology: tool for the unification of biology. Nat Genet. 2000, 25: 25-29. 10.1038/75556.PubMed CentralView ArticleGoogle Scholar
- Finn RD, Mistry J, Schuster-Bockler B, Griffiths-Jones S, Hollich V, Lassmann T, Moxon S, Marshall M, Khanna A, Durbin R, others: Pfam: clans, web tools and services. Nucleic Acids Res. 2006, 34: D247-D251. 10.1093/nar/gkj149.PubMed CentralPubMedView ArticleGoogle Scholar
- Carrington JC, Ambros V: Role of microRNAs in plant and animal development. Science. 2003, 301: 336-338. 10.1126/science.1085242.PubMedView ArticleGoogle Scholar
- Bartel DP: MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.PubMedView ArticleGoogle Scholar
- Zhang BH, Pan XP, Cannon CH, Cobb GP, Anderson TA: Conservation and divergence of plant microRNA genes. Plant J. 2006, 46: 243-259. 10.1111/j.1365-313X.2006.02697.x.PubMedView ArticleGoogle Scholar
- Zhang BH, Pan XP, Cox SB, Cobb GP, Anderson TA: Evidence that miRNAs are different from other RNAs. Cell Mol Life Sci. 2006, 63: 246-254. 10.1007/s00018-005-5467-7.PubMedView ArticleGoogle Scholar
- Griffiths-Jones S: The microRNA Registry. Nucleic Acids Res. 2004, 32: D109-D111. 10.1093/nar/gkh023.PubMed CentralPubMedView ArticleGoogle Scholar
- Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31: 3406-3415. 10.1093/nar/gkg595.PubMed CentralPubMedView ArticleGoogle Scholar
- Palatnik JF, Allen E, Wu XL, Schommer C, Schwab R, Carrington JC, Weigel D: Control of leaf morphogenesis by microRNAs. Nature. 2003, 425: 257-263. 10.1038/nature01958.PubMedView ArticleGoogle Scholar
- Gustafson AM, Allen E, Givan S, Smith D, Carrington JC, Kasschau KD: ASRP: the Arabidopsis Small RNA Project Database. Nucleic Acids Res. 2005, 33: D637-D640. 10.1093/nar/gki127.PubMed CentralPubMedView ArticleGoogle Scholar
- Chisholm ST, Coaker G, Day B, Staskawicz BJ: Host-microbe interactions: Shaping the evolution of the plant immune response. Cell. 2006, 124: 803-814. 10.1016/j.cell.2006.02.008.PubMedView ArticleGoogle Scholar
- Kushner DB, Lindenbach BD, Grdzelishvili VZ, Noueiry AO, Paul SM, Ahlquist P: Systematic, genome-wide identification of host genes affecting replication of a positive-strand RNA virus. Proc Natl Acad Sci USA. 2003, 100: 15764-15769. 10.1073/pnas.2536857100.PubMed CentralPubMedView ArticleGoogle Scholar
- Diaz-Pendon JA, Truniger V, Nieto C, Garcia-Mas J, Bendahmane A, Aranda MA: Advances in understanding recessive resistance to plant viruses. Mol Plant Pathol. 2004, 5: 223-233. 10.1111/j.1364-3703.2004.00223.x.PubMedView ArticleGoogle Scholar
- Robaglia C, Caranta C: Translation initiation factors: a weak link in plant RNA virus infection. Trends Plant Sci. 2006, 11: 40-45. 10.1016/j.tplants.2005.11.004.PubMedView ArticleGoogle Scholar
- Gomez-Gomez L, Boller T: FLS2: An LRR receptor-like kinase involved in the perception of the bacterial elicitor flagellin in Arabidopsis. Mol Cell. 2000, 5: 1003-1011. 10.1016/S1097-2765(00)80265-8.PubMedView ArticleGoogle Scholar
- Giovannoni JJ: Genetic regulation of fruit development and ripening. Plant Cell. 2004, 16: S170-S180. 10.1105/tpc.019158.PubMed CentralPubMedView ArticleGoogle Scholar
- Tanksley SD: The genetic, developmental, and molecular bases of fruit size and shape variation in tomato. Plant Cell. 2004, 16: S181-S189. 10.1105/tpc.018119.PubMed CentralPubMedView ArticleGoogle Scholar
- Pinyopich A, Ditta GS, Savidge B, Liljegren SJ, Baumann E, Wisman E, Yanofsky MF: Assessing the redundancy of MADS-box genes during carpel and ovule development. Nature. 2003, 424: 85-88. 10.1038/nature01741.PubMedView ArticleGoogle Scholar
- Dinneny JR, Weigel D, Yanofsky MF: A genetic framework for fruit patterning in Arabidopsis thaliana. Development. 2005, 132: 4687-4696. 10.1242/dev.02062.PubMedView ArticleGoogle Scholar
- Pitrat M: 2002 gene list for melon. Cucurbit Genet Coop Rep. 2002, 25: 76-93.Google Scholar
- Monforte AJ, Oliver M, Gonzalo MJ, Alvarez JM, Dolcet-Sanjuan R, Arus P: Identification of quantitative trait loci involved in fruit quality traits in melon (Cucumis melo L.). Theor Appl Genet. 2004, 108: 750-758. 10.1007/s00122-003-1483-x.PubMedView ArticleGoogle Scholar
- Frary A, Nesbitt TC, Frary A, Grandillo S, van der Knaap E, Cong B, Liu JP, Meller J, Elber R, Alpert KB, others: fw2.2: A quantitative trait locus key to the evolution of tomato fruit size. Science. 2000, 289: 85-88. 10.1126/science.289.5476.85.PubMedView ArticleGoogle Scholar
- Liu JP, Van Eck J, Cong B, Tanksley SD: A new class of regulatory genes underlying the cause of pear-shaped tomato fruit. Proc Natl Acad Sci USA. 2002, 99: 13302-13306. 10.1073/pnas.162485999.PubMed CentralPubMedView ArticleGoogle Scholar
- Manning K, Tor M, Poole M, Hong Y, Thompson AJ, King GJ, Giovannoni JJ, Seymour GB: A naturally occurring epigenetic mutation in a gene encoding an SBP-box transcription factor inhibits tomato fruit ripening. Nat Genet. 2006, 38: 948-952. 10.1038/ng1841.PubMedView ArticleGoogle Scholar
- Aranda MA, Escaler M, Wang D, Maule AJ: Induction of HSP70 and polyubiquitin expression associated with plant virus replication. Proc Natl Acad Sci USA. 1996, 93: 15289-15293. 10.1073/pnas.93.26.15289.PubMed CentralPubMedView ArticleGoogle Scholar
- Escaler M, Aranda MA, Roberts IM, Thomas CL, Maule AJ: A comparison between virus replication and abiotic stress (heat) as modifiers of host gene expression in pea. Mol Plant Pathol. 2000, 1: 159-167. 10.1046/j.1364-3703.2000.00020.x.PubMedView ArticleGoogle Scholar
- Whitham SA, Yang CL, Goodin MM: Global impact: Elucidating plant responses to viral infection. Mol Plant-Microbe Interact. 2006, 19: 1207-1215. 10.1094/MPMI-19-1207.PubMedView ArticleGoogle Scholar
- Aranda M, Maule A: Virus-induced host gene shutoff in animals and plants. Virology. 1998, 243: 261-267. 10.1006/viro.1998.9032.PubMedView ArticleGoogle Scholar
- Bonaldo MDF, Lennon G, Soares MB: Normalization and subtraction: Two approaches to facilitate gene discovery. Genome Res. 1996, 6: 791-806. 10.1101/gr.6.9.791.PubMedView ArticleGoogle Scholar
- Zhulidov PA, Bogdanova EA, Shcheglov AS, Vagner LL, Khaspekov GL, Kozhemyako VB, Matz MV, Meleshkevitch E, Moroz LL, Lukyanov SA, others: Simple cDNA normalization using kamchatka crab duplex-specific nuclease. Nucleic Acids Res. 2004, 32: e37-10.1093/nar/gnh031.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhang LD, Yuan DJ, Yu SW, Li ZG, Cao YF, Miao ZQ, Qian HM, Tang KX: Preference of simple sequence repeats in coding and non-coding regions of Arabidopsis thaliana. Bioinformatics. 2004, 20: 1081-1086. 10.1093/bioinformatics/bth043.PubMedView ArticleGoogle Scholar
- Finnegan EJ, Genger RK, Peacock WJ, Dennis ES: DNA methylation in plants. Annu Rev Plant Physiol Plant Mol Biol. 1998, 49: 223-247. 10.1146/annurev.arplant.49.1.223.PubMedView ArticleGoogle Scholar
- Morales M, Roig E, Monforte AJ, Arus P, Garcia-Mas J: Single-nucleotide polymorphisms detected in expressed sequence tags of melon (Cucumis melo L.). Genome. 2004, 47: 352-360.PubMedView ArticleGoogle Scholar
- Niu QW, Lin SS, Reyes JL, Chen KC, Wu HW, Yeh SD, Chua NH: Expression of artificial microRNAs in transgenic Arabidopsis thaliana confers virus resistance. Nat Biotechnol. 2006, 24: 1420-1428. 10.1038/nbt1255.PubMedView ArticleGoogle Scholar
- Schwab R, Ossowski S, Riester M, Warthmann N, Weigel D: Highly specific gene silencing by artificial microRNAs in Arabidopsis. Plant Cell. 2006, 18: 1121-1133. 10.1105/tpc.105.039834.PubMed CentralPubMedView ArticleGoogle Scholar
- Yamanaka T, Imai T, Satoh R, Kawashima A, Takahashi M, Tomita K, Kubota K, Meshi T, Naito S, Ishikawa M: Complete inhibition of tobamovirus multiplication by simultaneous mutations in two homologous host genes. J Virol. 2002, 76: 2491-2497. 10.1128/jvi.76.5.2491-2497.2002.PubMed CentralPubMedView ArticleGoogle Scholar
- Malamy JE: Intrinsic and environmental response pathways that regulate root system architecture. Plant Cell Environ. 2005, 28: 67-77. 10.1111/j.1365-3040.2005.01306.x.PubMedView ArticleGoogle Scholar
- Radi A, Dina P, Guy A: Expression of sarcotoxin IA gene via a root-specific tob promoter enhanced host resistance against parasitic weeds in tomato plants. Plant Cell Rep. 2006, 25: 297-303. 10.1007/s00299-005-0052-y.PubMedView ArticleGoogle Scholar
- Dias RDS, Pico B, Espinos A, Nuez F: Resistance to melon vine decline derived from Cucumis melo ssp agrestis: genetic analysis of root structure and root response. Plant Breeding. 2004, 123: 66-72. 10.1046/j.1439-0523.2003.00944.x.View ArticleGoogle Scholar
- Tomato expression database. [http://ted.bti.cornell.edu]
- Sakai H, Hua J, Chen QHG, Chang CR, Medrano LJ, Bleecker AB, Meyerowitz EM: ETR2 is an ETR1-like gene involved in ethylene signaling in Arabidopsis. Proc Natl Acad Sci USA. 1998, 95: 5812-5817. 10.1073/pnas.95.10.5812.PubMed CentralPubMedView ArticleGoogle Scholar
- Tieman DM, Klee HJ: Differential expression of two novel members of the tomato ethylene-receptor family. Plant Physiol. 1999, 120: 165-172. 10.1104/pp.120.1.165.PubMed CentralPubMedView ArticleGoogle Scholar
- Ronen G, Cohen M, Zamir D, Hirschberg J: Regulation of carotenoid biosynthesis during tomato fruit development: Expression of the gene for lycopene epsilon-cyclase is down-regulated during ripening and is elevated in the mutant Delta. Plant J. 1999, 17: 341-351. 10.1046/j.1365-313X.1999.00381.x.PubMedView ArticleGoogle Scholar
- Catala C, Rose JKC, Bennett AB: Auxin-regulated genes encoding cell wall-modifying proteins are expressed during early tomato fruit growth. Plant Physiol. 2000, 122: 527-534. 10.1104/pp.122.2.527.PubMed CentralPubMedView ArticleGoogle Scholar
- Esteva J, Nuez F: Field resistance to melon dieback in Cucumis melo L. Cucurbit Genet Coop. 1994, 17: 76-77.Google Scholar
- Iglesias A, Pico B, Nuez F: A temporal genetic analysis of disease resistance genes: resistance to melon vine decline derived from Cucumis melo var. agrestis. Plant Breeding. 2000, 119: 329-334. 10.1046/j.1439-0523.2000.00507.x.View ArticleGoogle Scholar
- Iglesias A, Picó B, Nuez F: Artificial inoculation methods and selection criteria for breeding melons against vine decline. Acta Hortic. 2000, 510: 155-162.View ArticleGoogle Scholar
- Picó B, Roig C, Fita A, Nuez F: Detección de Monosporascus cannonballus en raíces de melón mediante PCR cuantitativa en tiempo real. Acta Port Hortic. 2005, 7: 169-176.Google Scholar
- Hull S: Matthews's plant Virology. 2002, San Diego: Academic PressGoogle Scholar
- Sambrook J, Rusell DW: Molecular cloning. 2001, A laboratory manual. New York: CSHL PRESSGoogle Scholar
- Aranda MA, Escaler M, Thomas CL, Maule AJ: A heat shock transcription factor in pea is differentially controlled by heat and virus replication. Plant J. 1999, 20: 153-161. 10.1046/j.1365-313x.1999.00586.x.PubMedView ArticleGoogle Scholar
- PERL. [http://www.perl.org]
- MySQL. [http://www.mysql.com]
- PHP. [http://www.php.net]
- Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using Phred. I. Accuracy assessment. Genome Res. 1998, 8: 175-185.PubMedView ArticleGoogle Scholar
- Chou HH, Holmes MH: DNA sequence quality trimming and vector removal. Bioinformatics. 2001, 17: 1093-1104. 10.1093/bioinformatics/17.12.1093.PubMedView ArticleGoogle Scholar
- RepeatMasker. [http://www.repeatmasker.org]
- DFCI Gene Indices Software Tools. [http://compbio.dfci.harvard.edu/tgi/software/]
- NCBI's UniVec. [http://www.ncbi.nlm.nih.gov/VecScreen/UniVec.html]
- Sputnik. [http://capb.dbi.udel.edu/main/ssr-proj/ssr-file/description.html]
- ESTScan software. [http://www.isrec.isb-sib.ch/ftp-server/ESTScan/]
- TAIR: the Arabidopsis information resource. [http://www.arabidopsis.org/]
- Uniref. [http://www.ebi.ac.uk/uniref/]
- Eddy SR: Profile hidden Markov models. Bioinformatics. 1998, 14: 755-763. 10.1093/bioinformatics/14.9.755.PubMedView ArticleGoogle Scholar
- Pfam database. [http://www.sanger.ac.uk/Software/Pfam/]
- McInerney JO: GCUA (General Codon Usage Analysis). Bioinformatics. 1998, 14: 372-373. 10.1093/bioinformatics/14.4.372.PubMedView ArticleGoogle Scholar
- Reinhart BJ, Weinstein EG, Rhoades MW, Bartel B, Bartel DP: MicroRNAs in plants. Genes Dev. 2002, 16: 1616-1626. 10.1101/gad.1004402.PubMed CentralPubMedView ArticleGoogle Scholar
- Apweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R, Magrane M, others: UniProt: the Universal Protein knowledgebase. Nucleic Acids Res. 2004, 32: D115-D119. 10.1093/nar/gkh131.PubMed CentralPubMedView ArticleGoogle 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.