MELOGEN: an EST database for melon functional genomics

Genetic Dissection of ToLCNDV Resistance in Resistant Sources of Cucumis melo

Tomato leaf curl New Delhi virus (ToLCNDV) is a begomovirus causing significant melon (Cucumis melo) crop losses globally. This study aims to map the ToLCNDV resistance in the PI 414723 melon accession, previously identified and characterized through phenotypic studies, thereby exploring shared genomic regions with the established resistant source WM-7. In the present study, WM-7 and PI 414723 were crossed with the susceptible accessions 'Rochet' and 'Blanco' respectively, to generate F1 hybrids. These hybrids were self-pollinated to generate the populations for mapping the ToLCNDV resistance region and designing markers for marker-assisted selection. Disease evaluation included visual symptom scoring, viral-load quantification and tissue printing. Genotyping-by-sequencing and SNP markers were used for quantitative trait loci (QTL) mapping. For genetic analysis, qPCR and bulked segregant RNA-seq (BSR-seq) were performed. Gene expression was assessed using RNA-seq, and qRT-PCR was used for confirmation. The research narrows the candidate region for resistance in WM-7 and identifies overlapping QTLs on chromosome 11 in PI 414723, found in the region of the DNA primase large subunit. BSR-seq and expression analyses highlight potential regulatory roles of chromosome 2 in conferring resistance. Differential expression was confirmed for six genes in the candidate region on chromosome 2. This study confirms the existence of common resistance genes in PI 414723 and WM-7.

Insights into the early transcriptomic response against watermelon mosaic virus in melon

BACKGROUND

Watermelon mosaic virus (WMV) is one of the most prevalent viruses affecting melon worldwide. Recessive resistance to WMV in melon has previously been reported in the African accession TGR-1551. Moreover, the genomic regions associated to the resistance have also been described. Nevertheless, the transcriptomic response that might infer the resistance to this potyvirus has not been explored.

RESULTS

We have performed a comparative transcriptomic analysis using mock and WMV-inoculated plants of the susceptible cultivar "Bola de oro" (BO) and a resistant RIL (Recombinant inbred line) derived from the initial cross between "TGR-1551" and BO. In total, 616 genes were identified as differentially expressed and the weighted gene co-expression network analysis (WGCNA) detected 19 gene clusters (GCs), of which 7 were differentially expressed for the genotype x treatment interaction term. SNPs with a predicted high impact on the protein function were detected within the coding regions of most of the detected DEGs. Moreover, 3 and 16 DEGs were detected within the QTL regions previously described in chromosomes 11 and 5, respectively. In addition to these two specific genomic regions, we also observde large transcriptomic changes from genes spread across the genome in the resistant plants in response to the virus infection. This early response against WMV implied genes involved in plant-pathogen interaction, plant hormone signal transduction, the MAPK signaling pathway or ubiquitin mediated proteolysis, in detriment to the photosynthetic and basal metabolites pathways. Moreover, the gene MELO3C021395, which coded a mediator of RNA polymerase II transcription subunit 33A (MED33A), has been proposed as the candidate gene located on chromosome 11 conferring resistance to WMV.

CONCLUSIONS

The comparative transcriptomic analysis presented here showed that, even though the resistance to WMV in TGR-1551 has a recessive nature, it triggers an active defense response at a transcriptomic level, which involves broad-spectrum resistance mechanisms. Thus, this study represents a step forward on our understanding of the mechanisms underlaying WMV resistance in melon. In addition, it sheds light into a broader topic on the mechanisms of recessive resistances.

Development and Application of EST-SSR Markers in Cephalotaxus oliveri From Transcriptome Sequences

Cephalotaxus oliveri is an endemic conifer of China, which has medicinal and ornamental value. However, the limited molecular markers and genetic information are insufficient for further genetic studies of this species. In this study, we characterized and developed the EST-SSRs from transcriptome sequences for the first time. The results showed that a total of 5089 SSRs were identified from 36446 unigenes with a density of one SSR per 11.1 kb. The most common type was trinucleotide repeats, excluding mononucleotide repeats, followed by dinucleotide repeats. AAG/CTT and AT/AT exhibited the highest frequency in the trinucleotide and dinucleotide repeats, respectively. Of the identified SSRs, 671, 1125, and 1958 SSRs were located in CDS, 3′UTR, and 5′UTR, respectively. Functional annotation showed that the SSR-containing unigenes were involved in growth and development with various biological functions. Among successfully designed primer pairs, 238 primer pairs were randomly selected for amplification and validation of EST-SSR markers and 47 primer pairs were identified as polymorphic. Finally, 28 high-polymorphic primers were used for genetic analysis and revealed a moderate level of genetic diversity. Seven natural C. oliveri sampling sites were divided into two genetic groups. Furthermore, the 28 EST-SSRs had 96.43, 71.43, and 78.57% of transferability rate in Cephalotaxus fortune, Ametotaxus argotaenia, and Pseudotaxus chienii, respectively. These markers developed in this study lay the foundation for further genetic and adaptive evolution studies in C. oliveri and related species.

RNA-Seq Transcriptome Analysis Provides Candidate Genes for Resistance to Tomato Leaf Curl New Delhi Virus in Melon

Tomato leaf curl New Delhi virus (ToLCNDV) emerged in the Mediterranean Basin in 2012 as the first DNA bipartite begomovirus (Geminiviridae family), causing severe yield and economic losses in cucurbit crops. A major resistance locus was identified in the wild melon accession WM-7 (Cucumis melo kachri group), but the mechanisms involved in the resistant response remained unknown. In this work, we used RNA-sequencing to identify disease-associated genes that are differentially expressed in the course of ToLCNDV infection and could contribute to resistance. Transcriptomes of the resistant WM-7 genotype and the susceptible cultivar Piñonet Piel de Sapo (PS) (C. melo ibericus group) in ToLCNDV and mock inoculated plants were compared at four time points during infection (0, 3, 6, and 12 days post inoculation). Different gene expression patterns were observed over time in the resistant and susceptible genotypes in comparison to their respective controls. Differentially expressed genes (DEGs) in ToLCNDV-infected plants were classified using gene ontology (GO) terms, and genes of the categories transcription, DNA replication, and helicase activity were downregulated in WM-7 but upregulated in PS, suggesting that reduced activity of these functions reduces ToLCNDV replication and intercellular spread and thereby contributes to resistance. DEGs involved in the jasmonic acid signaling pathway, photosynthesis, RNA silencing, transmembrane, and sugar transporters entail adverse consequences for systemic infection in the resistant genotype, and lead to susceptibility in PS. The expression levels of selected candidate genes were validated by qRT-PCR to corroborate their differential expression upon ToLCNDV infection in resistant and susceptible melon. Furthermore, single nucleotide polymorphism (SNPs) with an effect on structural functionality of DEGs linked to the main QTLs for ToLCNDV resistance have been identified. The obtained results pinpoint cellular functions and candidate genes that are differentially expressed in a resistant and susceptible melon line in response to ToLCNDV, an information of great relevance for breeding ToLCNDV-resistant melon cultivars.

A mutation in the melon Vacuolar Protein Sorting 41prevents systemic infection of Cucumber mosaic virus

In the melon exotic accession PI 161375, the gene cmv1, confers recessive resistance to Cucumber mosaic virus (CMV) strains of subgroup II. cmv1 prevents the systemic infection by restricting the virus to the bundle sheath cells and impeding viral loading to the phloem. Here we report the fine mapping and cloning of cmv1. Screening of an F2 population reduced the cmv1 region to a 132 Kb interval that includes a Vacuolar Protein Sorting 41 gene. CmVPS41 is conserved among plants, animals and yeast and is required for post-Golgi vesicle trafficking towards the vacuole. We have validated CmVPS41 as the gene responsible for the resistance, both by generating CMV susceptible transgenic melon plants, expressing the susceptible allele in the resistant cultivar and by characterizing CmVPS41 TILLING mutants with reduced susceptibility to CMV. Finally, a core collection of 52 melon accessions allowed us to identify a single amino acid substitution (L348R) as the only polymorphism associated with the resistant phenotype. CmVPS41 is the first natural recessive resistance gene found to be involved in viral transport and its cellular function suggests that CMV might use CmVPS41 for its own transport towards the phloem.

De novo comparative transcriptome analysis of genes involved in fruit morphology of pumpkin cultivars with extreme size difference and development of EST-SSR markers

The genetic basis of fruit size and shape was investigated for the first time in Cucurbita species and genetic loci associated with fruit morphology have been identified. Although extensive genomic resources are available at present for tomato (Solanum lycopersicum), cucumber (Cucumis sativus), melon (Cucumis melo) and watermelon (Citrullus lanatus), genomic databases for Cucurbita species are limited. Recently, our group reported the generation of pumpkin (Cucurbita pepo) transcriptome databases from two contrasting cultivars with extreme fruit sizes. In the current study we used these databases to perform comparative transcriptome analysis in order to identify genes with potential roles in fruit morphology and fruit size. Differential Gene Expression (DGE) analysis between cv. 'Munchkin' (small-fruit) and cv. 'Big Moose' (large-fruit) revealed a variety of candidate genes associated with fruit morphology with significant differences in gene expression between the two cultivars. In addition, we have set the framework for generating EST-SSR markers, which discriminate different C. pepo cultivars and show transferability to related Cucurbitaceae species. The results of the present study will contribute to both further understanding the molecular mechanisms regulating fruit morphology and furthermore identifying the factors that determine fruit size. Moreover, they may lead to the development of molecular marker tools for selecting genotypes with desired morphological traits.

Analysis of codon usage patterns in Morus notabilis based on genome and transcriptome data

Codons play important roles in regulating gene expression levels and mRNA half-lives. However, codon usage and related studies in multicellular organisms still lag far behind those in unicellular organisms. In this study, we describe for the first time genome-wide patterns of codon bias in Morus notabilis (mulberry tree), and analyze genome-wide codon usage in 12 other species within the order Rosales. The codon usage of M. notabilis was affected by nucleotide composition, mutation pressure, nature selection, and gene expression level. Translational selection optimal codons were identified and highly expressed genes of M. notabilis tended to use the optimal codons. Genes with higher expression levels have shorter coding region and lower amino acid complexity. Housekeeping genes showed stronger translational selection, which, notably, was not caused by the large differences between the expression level of housekeeping genes and other genes.

Transcriptomic profiling of Melon necrotic spot virus-infected melon plants revealed virus strain and plant cultivar-specific alterations

Background

Viruses are among the most destructive and difficult to control plant pathogens. Melon (Cucumis melo L.) has become the model species for the agriculturally important Cucurbitaceae family. Approaches that take advantage of recently developed genomic tools in melon have been extremely useful for understanding viral pathogenesis and can contribute to the identification of target genes for breeding new resistant cultivars. In this work, we have used a recently described melon microarray for transcriptome profiling of two melon cultivars infected with two strains of Melon necrotic spot virus (MNSV) that only differ on their 3′-untranslated regions.

Results

Melon plant tissues from the cultivars Tendral or Planters Jumbo were locally infected with either MNSV-Mα5 or MNSV-Mα5/3’264 and analysed in a time-course experiment. Principal component and hierarchical clustering analyses identified treatment (healthy vs. infected) and sampling date (3 vs. 5 dpi) as the primary and secondary variables, respectively. Out of 7566 and 7074 genes deregulated by MNSV-Mα5 and MNSV-Mα5/3’264, 1851 and 1356, respectively, were strain-specific. Likewise, MNSV-Mα5/3’264 specifically deregulated 2925 and 1618 genes in Tendral and Planters Jumbo, respectively. The GO categories that were significantly affected were clearly different for the different virus/host combinations. Grouping genes according to their patterns of expression allowed for the identification of two groups that were specifically deregulated by MNSV-Mα5/3’264 with respect to MNSV-Mα5 in Tendral, and one group that was antagonistically regulated in Planters Jumbo vs. Tendral after MNSV-Mα5/3’264 infection. Genes in these three groups belonged to diverse functional classes, and no obvious regulatory commonalities were identified. When data on MNSV-Mα5/Tendral infections were compared to equivalent data on cucumber mosaic virus or watermelon mosaic virus infections, cytokinin-O-glucosyltransferase2 was identified as the only gene that was deregulated by all three viruses, with infection dynamics correlating with the amplitude of transcriptome remodeling.

Conclusions

Strain-specific changes, as well as cultivar-specific changes, were identified by profiling the transcriptomes of plants from two melon cultivars infected with two MNSV strains. No obvious regulatory features shared among deregulated genes have been identified, pointing toward regulation through differential functional pathways.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2772-5) contains supplementary material, which is available to authorized users.

Evolutionary Consequences of DNA Methylation in a Basal Metazoan

Gene body methylation (gbM) is an ancestral and widespread feature in Eukarya, yet its adaptive value and evolutionary implications remain unresolved. The occurrence of gbM within protein-coding sequences is particularly puzzling, because methylation causes cytosine hypermutability and hence is likely to produce deleterious amino acid substitutions. We investigate this enigma using an evolutionarily basal group of Metazoa, the stony corals (order Scleractinia, class Anthozoa, phylum Cnidaria). We show that patterns of coral gbM are similar to other invertebrate species, predicting wide and active transcription and slower sequence evolution. We also find a strong correlation between gbM and codon bias, resulting from systematic replacement of CpG bearing codons. We conclude that gbM has strong effects on codon evolution and speculate that this may influence establishment of optimal codons.

Transcriptome analysis and molecular marker discovery in Solanum incanum and S. aethiopicum, two close relatives of the common eggplant (Solanum melongena) with interest for breeding

Background

Solanum incanum is a close wild relative of S. melongena with high contents of bioactive phenolics and drought tolerance. S. aethiopicum is a cultivated African eggplant cross-compatible with S. melongena. Despite their great interest in S. melongena breeding programs, the genomic resources for these species are scarce.

Results

RNA-Seq was performed with NGS from pooled RNA of young leaf, floral bud and young fruit tissues, generating more than one hundred millions raw reads per species. The transcriptomes were assembled in 83,905 unigenes for S. incanum and in 87,084 unigenes for S. aethiopicum with an average length of 696 and 722 bp, respectively. The unigenes were structurally and functionally annotated based on comparison with public databases by using bioinformatic tools. The single nucleotide variant calling analysis (SNPs and INDELs) was performed by mapping our S. incanum and S. aethiopicum reads, as well as reads from S. melongena and S. torvum available on NCBI database (National Center for Biotechnology Information), against the eggplant genome. Both intraspecific and interspecific polymorphisms were identified and subsets of molecular markers were created for all species combinations. 36 SNVs were selected for validation in the S. incanum and S. aethiopicum accessions and 96 % were correctly amplified confirming the polymorphisms. In addition, 976 and 1,278 SSRs were identified in S. incanum and S. aethiopicum transcriptomes respectively, and a set of them were validated.

Conclusions

This work provides a broad insight into gene sequences and allelic variation in S. incanum and S. aethiopicum. This work is a first step toward better understanding of target genes involved in metabolic pathways relevant for eggplant breeding. The molecular markers detected in this study could be used across all the eggplant genepool, which is of interest for breeding programs as well as to perform marker-trait association and QTL analysis studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2631-4) contains supplementary material, which is available to authorized users.

De Novo Transcriptome Analysis of Cucumis melo L. var. makuwa

Oriental melon (Cucumis melo L. var. makuwa) is one of six subspecies of melon and is cultivated widely in East Asia, including China, Japan, and Korea. Although oriental melon is economically valuable in Asia and is genetically distinct from other subspecies, few reports of genome-scale research on oriental melon have been published. We generated 30.5 and 36.8 Gb of raw RNA sequence data from the female and male flowers, leaves, roots, and fruit of two oriental melon varieties, Korean landrace (KM) and Breeding line of NongWoo Bio Co. (NW), respectively. From the raw reads, 64,998 transcripts from KM and 100,234 transcripts from NW were de novo assembled. The assembled transcripts were used to identify molecular markers (e.g., single-nucleotide polymorphisms and simple sequence repeats), detect tissue-specific expressed genes, and construct a genetic linkage map. In total, 234 single-nucleotide polymorphisms and 25 simple sequence repeats were screened from 7,871 and 8,052 candidates, respectively, between the KM and NW varieties and used for construction of a genetic map with 94 F2 population specimens. The genetic linkage map consisted of 12 linkage groups, and 248 markers were assigned. These transcriptome and molecular marker data provide information useful for molecular breeding of oriental melon and further comparative studies of the Cucurbitaceae family.

Genome Sequence and Analysis of Peptoclostridium difficile Strain ZJCDC-S82

Peptoclostridium difficile (Clostridium difficile) is the major pathogen associated with infectious diarrhea in humans. Concomitant with the increased incidence of C. difficile infection worldwide, there is an increasing concern regarding this infection type. This study reports a draft assembly and detailed sequence analysis of C. difficile strain ZJCDC-S82. The de novo assembled genome was 4.19 Mb in size, which includes 4,013 protein-coding genes, 41 rRNA genes, and 84 tRNA genes. Along with the nuclear genome, we also assembled sequencing information for a single plasmid consisting of 11,930 nucleotides. Comparative genomic analysis of C. difficile ZJCDC-S82 and two other previously published strains, such as M120 and CD630, showed extensive similarity. Phylogenetic analysis revealed that genetic diversity among C. difficile strains was not influenced by geographic location. Evolutionary analysis suggested that four genes encoding surface proteins exhibited positive selection in C. difficile ZJCDC-S82. Codon usage analysis indicated that C. difficile ZJCDC-S82 had high codon usage bias toward A/U-ended codons. Furthermore, codon usage patterns in C. difficile ZJCDC-S82 were predominantly affected by mutation pressure. Our results provide detailed information pertaining to the C. difficile genome associated with a strain from mainland China. This analysis will facilitate the understanding of genomic diversity and evolution of C. difficile strains in this region.

Assessment of Suitable Reference Genes for Quantitative Gene Expression Studies in Melon Fruits

Melon (Cucumis melo L.) is an attractive model plant for investigating fruit development because of its morphological, physiological, and biochemical diversity. Quantification of gene expression by quantitative reverse transcription polymerase chain reaction (qRT-PCR) with stably expressed reference genes for normalization can effectively elucidate the biological functions of genes that regulate fruit development. However, the reference genes for data normalization in melon fruits have not yet been systematically validated. This study aims to assess the suitability of 20 genes for their potential use as reference genes in melon fruits. Expression variations of these genes were measured in 24 samples that represented different developmental stages of fertilized and parthenocarpic melon fruits by qRT-PCR analysis. GeNorm identified ribosomal protein L (CmRPL) and cytosolic ribosomal protein S15 (CmRPS15) as the best pair of reference genes, and as many as five genes including CmRPL, CmRPS15, TIP41-like family protein (CmTIP41), cyclophilin ROC7 (CmCYP7), and ADP ribosylation factor 1 (CmADP) were required for more reliable normalization. NormFinder ranked CmRPS15 as the best single reference gene, and RAN GTPase gene family (CmRAN) and TATA-box binding protein (CmTBP2) as the best combination of reference genes in melon fruits. Their effectiveness was further validated by parallel analyses on the activities of soluble acid invertase and sucrose phosphate synthase, and expression profiles of their respective encoding genes CmAIN2 and CmSPS1, as well as sucrose contents during melon fruit ripening. The validated reference genes will help to improve the accuracy of gene expression studies in melon fruits.

Comparative transcriptional profiling analysis of developing melon (Cucumis melo L.) fruit from climacteric and non-climacteric varieties

Background

In climacteric fruit-bearing species, the onset of fruit ripening is marked by a transient rise in respiration rate and autocatalytic ethylene production, followed by rapid deterioration in fruit quality. In non-climacteric species, there is no increase in respiration or ethylene production at the beginning or during fruit ripening. Melon is unusual in having climacteric and non-climacteric varieties, providing an interesting model system to compare both ripening types. Transcriptomic analysis of developing melon fruits from Védrantais and Dulce (climacteric) and Piel de sapo and PI 161375 (non-climacteric) varieties was performed to understand the molecular mechanisms that differentiate the two fruit ripening types.

Results

Fruits were harvested at 15, 25, 35 days after pollination and at fruit maturity. Transcript profiling was performed using an oligo-based microarray with 75 K probes. Genes linked to characteristic traits of fruit ripening were differentially expressed between climacteric and non-climacteric types, as well as several transcription factor genes and genes encoding enzymes involved in sucrose catabolism. The expression patterns of some genes in PI 161375 fruits were either intermediate between. Piel de sapo and the climacteric varieties, or more similar to the latter. PI 161375 fruits also accumulated some carotenoids, a characteristic trait of climacteric varieties.

Conclusions

Simultaneous changes in transcript abundance indicate that there is coordinated reprogramming of gene expression during fruit development and at the onset of ripening in both climacteric and non-climacteric fruits. The expression patterns of genes related to ethylene metabolism, carotenoid accumulation, cell wall integrity and transcriptional regulation varied between genotypes and was consistent with the differences in their fruit ripening characteristics. There were differences between climacteric and non-climacteric varieties in the expression of genes related to sugar metabolism suggesting that they may be potential determinants of sucrose content and post-harvest stability of sucrose levels in fruit. Several transcription factor genes were also identified that were differentially expressed in both types, implicating them in regulation of ripening behaviour. The intermediate nature of PI 161375 suggested that classification of melon fruit ripening behaviour into just two distinct types is an over-simplification, and that in reality there is a continuous spectrum of fruit ripening behaviour.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1649-3) contains supplementary material, which is available to authorized users.

A Comprehensive Analysis of Codon Usage Patterns in Blunt Snout Bream (Megalobrama amblycephala) Based on RNA-Seq Data

Blunt snout bream (Megalobrama amblycephala) is an important fish species for its delicacy and high economic value in China. Codon usage analysis could be helpful to understand its codon biology, mRNA translation and vertebrate evolution. Based on RNA-Seq data for M. amblycephala, high-frequency codons (CUG, AGA, GUG, CAG and GAG), as well as low-frequency ones (NUA and NCG codons) were identified. A total of 724 high-frequency codon pairs were observed. Meanwhile, 14 preferred and 199 avoided neighboring codon pairs were also identified, but bias was almost not shown with one or more intervening codons inserted between the same pairs. Codon usage bias in the regions close to start and stop codons indicated apparent heterogeneity, which even occurs in the flanking nucleotide sequence. Codon usage bias (RSCU and SCUO) was related to GC₃ (GC content of 3rd nucleotide in codon) bias. Six GO (Gene ontology) categories and the number of methylation targets were influenced by GC₃. Codon usage patterns comparison among 23 vertebrates showed species specificities by using GC contents, codon usage and codon context analysis. This work provided new insights into fish biology and new information for breeding projects.

CmMDb: a versatile database for Cucumis melo microsatellite markers and other horticulture crop research

Cucumis melo L. that belongs to Cucurbitaceae family ranks among one of the highest valued horticulture crops being cultivated across the globe. Besides its economical and medicinal importance, Cucumis melo L. is a valuable resource and model system for the evolutionary studies of cucurbit family. However, very limited numbers of molecular markers were reported for Cucumis melo L. so far that limits the pace of functional genomic research in melon and other similar horticulture crops. We developed the first whole genome based microsatellite DNA marker database of Cucumis melo L. and comprehensive web resource that aids in variety identification and physical mapping of Cucurbitaceae family. The Cucumis melo L. microsatellite database (CmMDb: http://65.181.125.102/cmmdb2/index.html) encompasses 39,072 SSR markers along with its motif repeat, motif length, motif sequence, marker ID, motif type and chromosomal locations. The database is featured with novel automated primer designing facility to meet the needs of wet lab researchers. CmMDb is a freely available web resource that facilitates the researchers to select the most appropriate markers for marker-assisted selection in melons and to improve breeding strategies.

The 2-C-methylerythritol 4-phosphate pathway in melon is regulated by specialized isoforms for the first and last steps

The 2-C-methyl-d-erythritol-4-phosphate (MEP) pathway provides the precursors for the biosynthesis of plastidial isoprenoids, which include the carotenoid pigments of many fruits. We have analysed the genes encoding the seven enzymes of the MEP pathway in melon (Cucumis melo L.) and determined that the first one, 1-deoxyxylulose 5-phosphate synthase (DXS), and the last one, 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase (HDR), are represented in the genome as a small gene family and paralogous pair, respectively. In the case of DXS, three genes encode functional DXS activities which fall into previously established type I (CmDXS1) and II (CmDXS2a and CmDXS2b) categories, while a fourth DXS-like gene belonging to the type III group did not encode a protein with DXS activity. Their expression patterns and phylogenies suggest that CmDXS1 is functionally specialized for developmental and photosynthetic processes, while CmDXS2a and CmDXS2b are induced in flowers and ripening fruit of orange- (but not white-) fleshed varieties, coinciding with β-carotene accumulation. This is the first instance connecting type II DXS genes to specialized isoprenoid biosynthesis in the fruit of an agronomically important species. Two HDR paralogues were shown to encode functional enzymes, although only CmHDR1 was highly expressed in the tissues and developmental stages tested. Phylogenetic analysis showed that in cucurbits such as melon, these HDR paralogues probably arose through individual gene duplications in a common angiosperm ancestor, mimicking a prior division in gymnosperms, while other flowering plants, including apple, soy, canola, and poplar, acquired HDR duplicates recently as homoeologues through large-scale genome duplications. We report the influence of gene duplication history on the regulation of the MEP pathway in melon and the role of specialized MEP-pathway isoforms in providing precursors for β-carotene production in orange-fleshed melon varieties.

Fatty acid profile and unigene-derived simple sequence repeat markers in tung tree (Vernicia fordii)

Tung tree (Vernicia fordii) provides the sole source of tung oil widely used in industry. Lack of fatty acid composition and molecular markers hinders biochemical, genetic and breeding research. The objectives of this study were to determine fatty acid profiles and develop unigene-derived simple sequence repeat (SSR) markers in tung tree. Fatty acid profiles of 41 accessions showed that the ratio of α-eleostearic acid was increasing continuously with a parallel trend to the amount of tung oil accumulation while the ratios of other fatty acids were decreasing in different stages of the seeds and that α-eleostearic acid (18∶3) consisted of 77% of the total fatty acids in tung oil. Transcriptome sequencing identified 81,805 unigenes from tung cDNA library constructed using seed mRNA and discovered 6,366 SSRs in 5,404 unigenes. The di- and tri-nucleotide microsatellites accounted for 92% of the SSRs with AG/CT and AAG/CTT being the most abundant SSR motifs. Fifteen polymorphic genic-SSR markers were developed from 98 unigene loci tested in 41 cultivated tung accessions by agarose gel and capillary electrophoresis. Genbank database search identified 10 of them putatively coding for functional proteins. Quantitative PCR demonstrated that all 15 polymorphic SSR-associated unigenes were expressed in tung seeds and some of them were highly correlated with oil composition in the seeds. Dendrogram revealed that most of the 41 accessions were clustered according to the geographic region. These new polymorphic genic-SSR markers will facilitate future studies on genetic diversity, molecular fingerprinting, comparative genomics and genetic mapping in tung tree. The lipid profiles in the seeds of 41 tung accessions will be valuable for biochemical and breeding studies.

Codon usage patterns in Chinese bayberry (Myrica rubra) based on RNA-Seq data

BACKGROUND

Codon usage analysis has been a classical topic for decades and has significances for studies of evolution, mRNA translation, and new gene discovery, etc. While the codon usage varies among different members of the plant kingdom, indicating the necessity for species-specific study, this work has mostly been limited to model organisms. Recently, the development of deep sequencing, especial RNA-Seq, has made it possible to carry out studies in non-model species.

RESULT

RNA-Seq data of Chinese bayberry was analyzed to investigate the bias of codon usage and codon pairs. High frequency codons (AGG, GCU, AAG and GAU), as well as low frequency ones (NCG and NUA codons) were identified, and 397 high frequency codon pairs were observed. Meanwhile, 26 preferred and 141 avoided neighboring codon pairs were also identified, which showed more significant bias than the same pairs with one or more intervening codons. Codon patterns were also analyzed at the plant kingdom, organism and gene levels. Changes during plant evolution were evident using RSCU (relative synonymous codon usage), which was even more significant than GC3s (GC content of 3rd synonymous codons). Nine GO categories were differentially and independently influenced by CAI (codon adaptation index) or GC3s, especially in 'Molecular function' category. Within a gene, the average CAI increased from 0.720 to 0.785 in the first 50 codons, and then more slowly thereafter. Furthermore, the preferred as well as avoided codons at the position just following the start codon AUG were identified and discussed in relation to the key positions in Kozak sequences.

CONCLUSION

A comprehensive codon usage Table and number of high-frequency codon pairs were established. Bias in codon usage as well as in neighboring codon pairs was observed, and the significance of this in avoiding DNA mutation, increasing protein production and regulating protein synthesis rate was proposed. Codon usage patterns at three levels were revealed and the significance in plant evolution analysis, gene function classification, and protein translation start site predication were discussed. This work promotes the study of codon biology, and provides some reference for analysis and comprehensive application of RNA-Seq data from other non-model species.

Genetic and chemical characterization of an EMS induced mutation in Cucumis melo CRTISO gene

In order to broaden the available genetic variation of melon, we developed an ethyl methanesulfonate mutation library in an orange-flesh 'Charentais' type melon line that accumulates β-carotene. One mutagenized M2 family segregated for a novel recessive trait, a yellow-orange fruit flesh ('yofI'). HPLC analysis revealed that 'yofI' accumulates pro-lycopene (tetra-cis-lycopene) as its major fruit pigment. The altered carotenoid composition of 'yofI' is associated with a significant change of the fruit aroma since cleavage of β-carotene yields different apocarotenoids than the cleavage of pro-lycopene. Normally, pro-lycopene is further isomerized by CRTISO (carotenoid isomerase) to yield all-trans-lycopene, which is further cyclized to β-carotene in melon fruit. Cloning and sequencing of 'yofI' CRTISO identified two mRNA sequences which lead to truncated forms of CRTISO. Sequencing of the genomic CRTISO identified an A-T transversion in 'yofI' which leads to a premature STOP codon. The early carotenoid pathway genes were up regulated in yofI fruit causing accumulation of other intermediates such as phytoene and ζ-carotene. Total carotenoid levels are only slightly increased in the mutant. Mutants accumulating pro-lycopene have been reported in both tomato and watermelon fruits, however, this is the first report of a non-lycopene accumulating fruit showing this phenomenon.

Assessment of a 1H high-resolution magic angle spinning NMR spectroscopy procedure for free sugars quantification in intact plant tissue

In most plants, sucrose is the primary product of photosynthesis, the transport form of assimilated carbon, and also one of the main factors determining sweetness in fresh fruits. Traditional methods for sugar quantification (mainly sucrose, glucose and fructose) require obtaining crude plant extracts, which sometimes involve substantial sample manipulation, making the process time-consuming and increasing the risk of sample degradation. Here, we describe and validate a fast method to determine sugar content in intact plant tissue by using high-resolution magic angle spinning nuclear magnetic resonance spectroscopy (HR-MAS NMR). The HR-MAS NMR method was used for quantifying sucrose, glucose and fructose in mesocarp tissues from melon fruits (Cucumis melo var. reticulatus and Cucumis melo var. cantalupensis). The resulting sugar content varied among individual melons, ranging from 1.4 to 7.3 g of sucrose, 0.4-2.5 g of glucose; and 0.73-2.83 g of fructose (values per 100 g fw). These values were in agreement with those described in the literature for melon fruit tissue, and no significant differences were found when comparing them with those obtained using the traditional, enzymatic procedure, on melon tissue extracts. The HR-MAS NMR method offers a fast (usually <30 min) and sensitive method for sugar quantification in intact plant tissues, it requires a small amount of tissue (typically 50 mg fw) and avoids the interferences and risks associated with obtaining plant extracts. Furthermore, this method might also allow the quantification of additional metabolites detectable in the plant tissue NMR spectrum.

SNP genotyping in melons: genetic variation, population structure, and linkage disequilibrium

Novel sequencing technologies were recently used to generate sequences from multiple melon (Cucumis melo L.) genotypes, enabling the in silico identification of large single nucleotide polymorphism (SNP) collections. In order to optimize the use of these markers, SNP validation and large-scale genotyping are necessary. In this paper, we present the first validated design for a genotyping array with 768 SNPs that are evenly distributed throughout the melon genome. This customized Illumina GoldenGate assay was used to genotype a collection of 74 accessions, representing most of the botanical groups of the species. Of the assayed loci, 91 % were successfully genotyped. The array provided a large number of polymorphic SNPs within and across accessions. This set of SNPs detected high levels of variation in accessions from this crop's center of origin as well as from several other areas of melon diversification. Allele distribution throughout the genome revealed regions that distinguished between the two main groups of cultivated accessions (inodorus and cantalupensis). Population structure analysis showed a subdivision into five subpopulations, reflecting the history of the crop. A considerably low level of LD was detected, which decayed rapidly within a few kilobases. Our results show that the GoldenGate assay can be used successfully for high-throughput SNP genotyping in melon. Since many of the genotyped accessions are currently being used as the parents of breeding populations in various programs, this set of mapped markers could be used for future mapping and breeding efforts.

De novo sequence assembly and characterization of Lycoris aurea transcriptome using GS FLX titanium platform of 454 pyrosequencing

Background

Lycoris aurea, also called Golden Magic Lily, is an ornamentally and medicinally important species of the Amaryllidaceae family. To date, the sequencing of its whole genome is unavailable as a non-model organism. Transcriptomic information is also scarce for this species. In this study, we performed de novo transcriptome sequencing to produce the first comprehensive expressed sequence tag (EST) dataset for L. aurea using high-throughput sequencing technology.

Methodology and Principal Findings

Total RNA was isolated from leaves with sodium nitroprusside (SNP), salicylic acid (SA), or methyl jasmonate (MeJA) treatment, stems, and flowers at the bud, blooming, and wilting stages. Equal quantities of RNA from each tissue and stage were pooled to construct a cDNA library. Using 454 pyrosequencing technology, a total of 937,990 high quality reads (308.63 Mb) with an average read length of 329 bp were generated. Clustering and assembly of these reads produced a non-redundant set of 141,111 unique sequences, comprising 24,604 contigs and 116,507 singletons. All of the unique sequences were involved in the biological process, cellular component and molecular function categories by GO analysis. Potential genes and their functions were predicted by KEGG pathway mapping and COG analysis. Based on our sequence analysis and published literatures, many putative genes involved in Amaryllidaceae alkaloids synthesis, including PAL, TYDC OMT, NMT, P450, and other potentially important candidate genes, were identified for the first time in this Lycoris. Furthermore, 6,386 SSRs and 18,107 high-confidence SNPs were identified in this EST dataset.

Conclusions

The transcriptome provides an invaluable new data for a functional genomics resource and future biological research in L. aurea. The molecular markers identified in this study will provide a material basis for future genetic linkage and quantitative trait loci analyses, and will provide useful information for functional genomic research in future.

Generation and analysis of expressed sequence tags in the extreme large genomes Lilium and Tulipa

BACKGROUND

Bulbous flowers such as lily and tulip (Liliaceae family) are monocot perennial herbs that are economically very important ornamental plants worldwide. However, there are hardly any genetic studies performed and genomic resources are lacking. To build genomic resources and develop tools to speed up the breeding in both crops, next generation sequencing was implemented. We sequenced and assembled transcriptomes of four lily and five tulip genotypes using 454 pyro-sequencing technology.

RESULTS

Successfully, we developed the first set of 81,791 contigs with an average length of 514 bp for tulip, and enriched the very limited number of 3,329 available ESTs (Expressed Sequence Tags) for lily with 52,172 contigs with an average length of 555 bp. The contigs together with singletons covered on average 37% of lily and 39% of tulip estimated transcriptome. Mining lily and tulip sequence data for SSRs (Simple Sequence Repeats) showed that di-nucleotide repeats were twice more abundant in UTRs (UnTranslated Regions) compared to coding regions, while tri-nucleotide repeats were equally spread over coding and UTR regions. Two sets of single nucleotide polymorphism (SNP) markers suitable for high throughput genotyping were developed. In the first set, no SNPs flanking the target SNP (50 bp on either side) were allowed. In the second set, one SNP in the flanking regions was allowed, which resulted in a 2 to 3 fold increase in SNP marker numbers compared with the first set. Orthologous groups between the two flower bulbs: lily and tulip (12,017 groups) and among the three monocot species: lily, tulip, and rice (6,900 groups) were determined using OrthoMCL. Orthologous groups were screened for common SNP markers and EST-SSRs to study synteny between lily and tulip, which resulted in 113 common SNP markers and 292 common EST-SSR. Lily and tulip contigs generated were annotated and described according to Gene Ontology terminology.

CONCLUSIONS

Two transcriptome sets were built that are valuable resources for marker development, comparative genomic studies and candidate gene approaches. Next generation sequencing of leaf transcriptome is very effective; however, deeper sequencing and using more tissues and stages is advisable for extended comparative studies.

Root transcriptional responses of two melon genotypes with contrasting resistance to Monosporascus cannonballus (Pollack et Uecker) infection

BACKGROUND

Monosporascus cannonballus is the main causal agent of melon vine decline disease. Several studies have been carried out mainly focused on the study of the penetration of this pathogen into melon roots, the evaluation of symptoms severity on infected roots, and screening assays for breeding programs. However, a detailed molecular view on the early interaction between M. cannonballus and melon roots in either susceptible or resistant genotypes is lacking. In the present study, we used a melon oligo-based microarray to investigate the gene expression responses of two melon genotypes, Cucumis melo 'Piel de sapo' ('PS') and C. melo 'Pat 81', with contrasting resistance to the disease. This study was carried out at 1 and 3 days after infection (DPI) by M. cannonballus.

RESULTS

Our results indicate a dissimilar behavior of the susceptible vs. the resistant genotypes from 1 to 3 DPI. 'PS' responded with a more rapid infection response than 'Pat 81' at 1 DPI. At 3 DPI the total number of differentially expressed genes identified in 'PS' declined from 451 to 359, while the total number of differentially expressed transcripts in 'Pat 81' increased from 187 to 849. Several deregulated transcripts coded for components of Ca2+ and jasmonic acid (JA) signalling pathways, as well as for other proteins related to defence mechanisms. Transcriptional differences in the activation of the JA-mediated response in 'Pat 81' compared to 'PS' suggested that JA response might be partially responsible for their observed differences in resistance.

CONCLUSIONS

As a result of this study we have identified for the first time a set of candidate genes involved in the root response to the infection of the pathogen causing melon vine decline. This information is useful for understanding the disease progression and resistance mechanisms few days after inoculation.

Why assembling plant genome sequences is so challenging

In spite of the biological and economic importance of plants, relatively few plant species have been sequenced. Only the genome sequence of plants with relatively small genomes, most of them angiosperms, in particular eudicots, has been determined. The arrival of next-generation sequencing technologies has allowed the rapid and efficient development of new genomic resources for non-model or orphan plant species. But the sequencing pace of plants is far from that of animals and microorganisms. This review focuses on the typical challenges of plant genomes that can explain why plant genomics is less developed than animal genomics. Explanations about the impact of some confounding factors emerging from the nature of plant genomes are given. As a result of these challenges and confounding factors, the correct assembly and annotation of plant genomes is hindered, genome drafts are produced, and advances in plant genomics are delayed.

Melon RNA interference (RNAi) lines silenced for Cm-eIF4E show broad virus resistance

Efficient and sustainable control of plant viruses may be achieved using genetically resistant crop varieties, although resistance genes are not always available for each pathogen; in this regard, the identification of new genes that are able to confer broad-spectrum and durable resistance is highly desirable. Recently, the cloning and characterization of recessive resistance genes from different plant species has pointed towards eukaryotic translation initiation factors (eIF) of the 4E family as factors required for the multiplication of many different viruses. Thus, we hypothesized that eIF4E may control the susceptibility of melon (Cucumis melo L.) to a broad range of viruses. To test this hypothesis, Cm-eIF4E knockdown melon plants were generated by the transformation of explants with a construct that was designed to induce the silencing of this gene, and the plants from T2 generations were genetically and phenotypically characterized. In transformed plants, Cm-eIF4E was specifically silenced, as identified by the decreased accumulation of Cm-eIF4E mRNA and the appearance of small interfering RNAs derived from the transgene, whereas the Cm-eIF(iso)4E mRNA levels remained unaffected. We challenged these transgenic melon plants with eight agronomically important melon-infecting viruses, and identified that they were resistant to Cucumber vein yellowing virus (CVYV), Melon necrotic spot virus (MNSV), Moroccan watermelon mosaic virus (MWMV) and Zucchini yellow mosaic virus (ZYMV), indicating that Cm-eIF4E controls melon susceptibility to these four viruses. Therefore, Cm-eIF4E is an efficient target for the identification of new resistance alleles able to confer broad-spectrum virus resistance in melon.

Transcriptome sequencing for SNP discovery across Cucumis melo

Background

Melon (Cucumis melo L.) is a highly diverse species that is cultivated worldwide. Recent advances in massively parallel sequencing have begun to allow the study of nucleotide diversity in this species. The Sanger method combined with medium-throughput 454 technology were used in a previous study to analyze the genetic diversity of germplasm representing 3 botanical varieties, yielding a collection of about 40,000 SNPs distributed in 14,000 unigenes. However, the usefulness of this resource is limited as the sequenced genotypes do not represent the whole diversity of the species, which is divided into two subspecies with many botanical varieties variable in plant, flowering, and fruit traits, as well as in stress response. As a first step to extensively document levels and patterns of nucleotide variability across the species, we used the high-throughput SOLiD™ system to resequence the transcriptomes of a set of 67 genotypes that had previously been selected from a core collection representing the extant variation of the entire species.

Results

The deep transcriptome resequencing of all of the genotypes, grouped into 8 pools (wild African agrestis, Asian agrestis and acidulus, exotic Far Eastern conomon, Indian momordica and Asian dudaim and flexuosus, commercial cantalupensis, subsp. melo Asian and European landraces, Spanish inodorus landraces, and Piel de Sapo breeding lines) yielded about 300 M reads. Short reads were mapped to the recently generated draft genome assembly of the DHL line Piel de Sapo (inodorus) x Songwhan Charmi (conomon) and to a new version of melon transcriptome. Regions with at least 6X coverage were used in SNV calling, generating a melon collection with 303,883 variants. These SNVs were dispersed across the entire C. melo genome, and distributed in 15,064 annotated genes. The number and variability of in silico SNVs differed considerably between pools. Our finding of higher genomic diversity in wild and exotic agrestis melons from India and Africa as compared to commercial cultivars, cultigens and landraces from Eastern Europe, Western Asia and the Mediterranean basin is consistent with the evolutionary history proposed for the species. Group-specific SNVs that will be useful in introgression programs were also detected. In a sample of 143 selected putative SNPs, we verified 93% of the polymorphisms in a panel of 78 genotypes.

Conclusions

This study provides the first comprehensive resequencing data for wild, exotic, and cultivated (landraces and commercial) melon transcriptomes, yielding the largest melon SNP collection available to date and representing a notable sample of the species diversity. This data provides a valuable resource for creating a catalog of allelic variants of melon genes and it will aid in future in-depth studies of population genetics, marker-assisted breeding, and gene identification aimed at developing improved varieties.

Transcriptome sequencing for SNP discovery across Cucumis melo

BACKGROUND

Melon (Cucumis melo L.) is a highly diverse species that is cultivated worldwide. Recent advances in massively parallel sequencing have begun to allow the study of nucleotide diversity in this species. The Sanger method combined with medium-throughput 454 technology were used in a previous study to analyze the genetic diversity of germplasm representing 3 botanical varieties, yielding a collection of about 40,000 SNPs distributed in 14,000 unigenes. However, the usefulness of this resource is limited as the sequenced genotypes do not represent the whole diversity of the species, which is divided into two subspecies with many botanical varieties variable in plant, flowering, and fruit traits, as well as in stress response. As a first step to extensively document levels and patterns of nucleotide variability across the species, we used the high-throughput SOLiD™ system to resequence the transcriptomes of a set of 67 genotypes that had previously been selected from a core collection representing the extant variation of the entire species.

RESULTS

The deep transcriptome resequencing of all of the genotypes, grouped into 8 pools (wild African agrestis, Asian agrestis and acidulus, exotic Far Eastern conomon, Indian momordica and Asian dudaim and flexuosus, commercial cantalupensis, subsp. melo Asian and European landraces, Spanish inodorus landraces, and Piel de Sapo breeding lines) yielded about 300 M reads. Short reads were mapped to the recently generated draft genome assembly of the DHL line Piel de Sapo (inodorus) x Songwhan Charmi (conomon) and to a new version of melon transcriptome. Regions with at least 6X coverage were used in SNV calling, generating a melon collection with 303,883 variants. These SNVs were dispersed across the entire C. melo genome, and distributed in 15,064 annotated genes. The number and variability of in silico SNVs differed considerably between pools. Our finding of higher genomic diversity in wild and exotic agrestis melons from India and Africa as compared to commercial cultivars, cultigens and landraces from Eastern Europe, Western Asia and the Mediterranean basin is consistent with the evolutionary history proposed for the species. Group-specific SNVs that will be useful in introgression programs were also detected. In a sample of 143 selected putative SNPs, we verified 93% of the polymorphisms in a panel of 78 genotypes.

CONCLUSIONS

This study provides the first comprehensive resequencing data for wild, exotic, and cultivated (landraces and commercial) melon transcriptomes, yielding the largest melon SNP collection available to date and representing a notable sample of the species diversity. This data provides a valuable resource for creating a catalog of allelic variants of melon genes and it will aid in future in-depth studies of population genetics, marker-assisted breeding, and gene identification aimed at developing improved varieties.

Diversity in expression of phosphorus (P) responsive genes in Cucumis melo L

BACKGROUND

Phosphorus (P) is a major limiting nutrient for plant growth in many soils. Studies in model species have identified genes involved in plant adaptations to low soil P availability. However, little information is available on the genetic bases of these adaptations in vegetable crops. In this respect, sequence data for melon now makes it possible to identify melon orthologues of candidate P responsive genes, and the expression of these genes can be used to explain the diversity in the root system adaptation to low P availability, recently observed in this species.

METHODOLOGY AND FINDINGS

Transcriptional responses to P starvation were studied in nine diverse melon accessions by comparing the expression of eight candidate genes (Cm-PAP10.1, Cm-PAP10.2, Cm-RNS1, Cm-PPCK1, Cm-transferase, Cm-SQD1, Cm-DGD1 and Cm-SPX2) under P replete and P starved conditions. Differences among melon accessions were observed in response to P starvation, including differences in plant morphology, P uptake, P use efficiency (PUE) and gene expression. All studied genes were up regulated under P starvation conditions. Differences in the expression of genes involved in P mobilization and remobilization (Cm-PAP10.1, Cm-PAP10.2 and Cm-RNS1) under P starvation conditions explained part of the differences in P uptake and PUE among melon accessions. The levels of expression of the other studied genes were diverse among melon accessions, but contributed less to the phenotypical response of the accessions.

CONCLUSIONS

This is the first time that these genes have been described in the context of P starvation responses in melon. There exists significant diversity in gene expression levels and P use efficiency among melon accessions as well as significant correlations between gene expression levels and phenotypical measurements.

Microarray analysis shows that recessive resistance to Watermelon mosaic virus in melon is associated with the induction of defense response genes

Resistance to Watermelon mosaic virus (WMV) in melon (Cucumis melo L.) accession TGR-1551 is characterized by a significant reduction in virus titer, and is inherited as a recessive, loss-of-susceptibility allele. We measured virus RNA accumulation in TGR-1551 plants and a susceptible control ('Tendral') by real-time quantitative polymerase chain reaction, and also profiled the expression of 17,443 unigenes represented on a melon microarray over a 15-day time course. The virus accumulated to higher levels in cotyledons of the resistant variety up to 9 days postinoculation (dpi) but, thereafter, levels increased in the susceptible variety while those in the resistant variety declined. Microarray experiments looking at the early response to infection (1 and 3 dpi), as well as responses after 7 and 15 dpi, revealed more profound transcriptomic changes in resistant plants than susceptible ones. The gene expression profiles revealed deep and extensive transcriptome remodeling in TGR-1551 plants, often involving genes with pathogen response functions. Overall, our data suggested that resistance to WMV in TGR-1551 melon plants is associated with a defense response, which contrasts with the recessive nature of the resistance trait.

Determination of the melon chloroplast and mitochondrial genome sequences reveals that the largest reported mitochondrial genome in plants contains a significant amount of DNA having a nuclear origin

Background

The melon belongs to the Cucurbitaceae family, whose economic importance among vegetable crops is second only to Solanaceae. The melon has a small genome size (454 Mb), which makes it suitable for molecular and genetic studies. Despite similar nuclear and chloroplast genome sizes, cucurbits show great variation when their mitochondrial genomes are compared. The melon possesses the largest plant mitochondrial genome, as much as eight times larger than that of other cucurbits.

Results

The nucleotide sequences of the melon chloroplast and mitochondrial genomes were determined. The chloroplast genome (156,017 bp) included 132 genes, with 98 single-copy genes dispersed between the small (SSC) and large (LSC) single-copy regions and 17 duplicated genes in the inverted repeat regions (IRa and IRb). A comparison of the cucumber and melon chloroplast genomes showed differences in only approximately 5% of nucleotides, mainly due to short indels and SNPs. Additionally, 2.74 Mb of mitochondrial sequence, accounting for 95% of the estimated mitochondrial genome size, were assembled into five scaffolds and four additional unscaffolded contigs. An 84% of the mitochondrial genome is contained in a single scaffold. The gene-coding region accounted for 1.7% (45,926 bp) of the total sequence, including 51 protein-coding genes, 4 conserved ORFs, 3 rRNA genes and 24 tRNA genes. Despite the differences observed in the mitochondrial genome sizes of cucurbit species, Citrullus lanatus (379 kb), Cucurbita pepo (983 kb) and Cucumis melo (2,740 kb) share 120 kb of sequence, including the predicted protein-coding regions. Nevertheless, melon contained a high number of repetitive sequences and a high content of DNA of nuclear origin, which represented 42% and 47% of the total sequence, respectively.

Conclusions

Whereas the size and gene organisation of chloroplast genomes are similar among the cucurbit species, mitochondrial genomes show a wide variety of sizes, with a non-conserved structure both in gene number and organisation, as well as in the features of the noncoding DNA. The transfer of nuclear DNA to the melon mitochondrial genome and the high proportion of repetitive DNA appear to explain the size of the largest mitochondrial genome reported so far.

Towards a TILLING platform for functional genomics in Piel de Sapo melons

BACKGROUND

The availability of genetic and genomic resources for melon has increased significantly, but functional genomics resources are still limited for this crop. TILLING is a powerful reverse genetics approach that can be utilized to generate novel mutations in candidate genes. A TILLING resource is available for cantalupensis melons, but not for inodorus melons, the other main commercial group.

RESULTS

A new ethyl methanesulfonate-mutagenized (EMS) melon population was generated for the first time in an andromonoecious non-climacteric inodorus Piel de Sapo genetic background. Diverse mutant phenotypes in seedlings, vines and fruits were observed, some of which were of possible commercial interest. The population was first screened for mutations in three target genes involved in disease resistance and fruit quality (Cm-PDS, Cm-eIF4E and Cm-eIFI(iso)4E). The same genes were also tilled in the available monoecious and climacteric cantalupensis EMS melon population. The overall mutation density in this first Piel de Sapo TILLING platform was estimated to be 1 mutation/1.5 Mb by screening four additional genes (Cm-ACO1, Cm-NOR, Cm-DET1 and Cm-DHS). Thirty-three point mutations were found for the seven gene targets, six of which were predicted to have an impact on the function of the protein. The genotype/phenotype correlation was demonstrated for a loss-of-function mutation in the Phytoene desaturase gene, which is involved in carotenoid biosynthesis.

CONCLUSIONS

The TILLING approach was successful at providing new mutations in the genetic background of Piel de Sapo in most of the analyzed genes, even in genes for which natural variation is extremely low. This new resource will facilitate reverse genetics studies in non-climacteric melons, contributing materially to future genomic and breeding studies.

Analysis of the melon (Cucumis melo) small RNAome by high-throughput pyrosequencing

BACKGROUND

Melon (Cucumis melo L.) is a commercially important fruit crop that is cultivated worldwide. The melon research community has recently benefited from the determination of a complete draft genome sequence and the development of associated genomic tools, which have allowed us to focus on small RNAs (sRNAs). These are short, non-coding RNAs 21-24 nucleotides in length with diverse physiological roles. In plants, they regulate gene expression and heterochromatin assembly, and control protection against virus infection. Much remains to be learned about the role of sRNAs in melon.

RESULTS

We constructed 10 sRNA libraries from two stages of developing ovaries, fruits and photosynthetic cotyledons infected with viruses, and carried out high-throughput pyrosequencing. We catalogued and analysed the melon sRNAs, resulting in the identification of 26 known miRNA families (many conserved with other species), the prediction of 84 melon-specific miRNA candidates, the identification of trans-acting siRNAs, and the identification of chloroplast, mitochondrion and transposon-derived sRNAs. In silico analysis revealed more than 400 potential targets for the conserved and novel miRNAs.

CONCLUSION

We have discovered and analysed a large number of conserved and melon-specific sRNAs, including miRNAs and their potential target genes. This provides insight into the composition and function of the melon small RNAome, and paves the way towards an understanding of sRNA-mediated processes that regulate melon fruit development and melon-virus interactions.

EuroPineDB: a high-coverage web database for maritime pine transcriptome

Background

Pinus pinaster is an economically and ecologically important species that is becoming a woody gymnosperm model. Its enormous genome size makes whole-genome sequencing approaches are hard to apply. Therefore, the expressed portion of the genome has to be characterised and the results and annotations have to be stored in dedicated databases.

Description

EuroPineDB is the largest sequence collection available for a single pine species, Pinus pinaster (maritime pine), since it comprises 951 641 raw sequence reads obtained from non-normalised cDNA libraries and high-throughput sequencing from adult (xylem, phloem, roots, stem, needles, cones, strobili) and embryonic (germinated embryos, buds, callus) maritime pine tissues. Using open-source tools, sequences were optimally pre-processed, assembled, and extensively annotated (GO, EC and KEGG terms, descriptions, SNPs, SSRs, ORFs and InterPro codes). As a result, a 10.5× P. pinaster genome was covered and assembled in 55 322 UniGenes. A total of 32 919 (59.5%) of P. pinaster UniGenes were annotated with at least one description, revealing at least 18 466 different genes. The complete database, which is designed to be scalable, maintainable, and expandable, is freely available at: http://www.scbi.uma.es/pindb/. It can be retrieved by gene libraries, pine species, annotations, UniGenes and microarrays (i.e., the sequences are distributed in two-colour microarrays; this is the only conifer database that provides this information) and will be periodically updated. Small assemblies can be viewed using a dedicated visualisation tool that connects them with SNPs. Any sequence or annotation set shown on-screen can be downloaded. Retrieval mechanisms for sequences and gene annotations are provided.

Conclusions

The EuroPineDB with its integrated information can be used to reveal new knowledge, offers an easy-to-use collection of information to directly support experimental work (including microarray hybridisation), and provides deeper knowledge on the maritime pine transcriptome.

Analysis of expressed sequence tags generated from full-length enriched cDNA libraries of melon

BACKGROUND

Melon (Cucumis melo), an economically important vegetable crop, belongs to the Cucurbitaceae family which includes several other important crops such as watermelon, cucumber, and pumpkin. It has served as a model system for sex determination and vascular biology studies. However, genomic resources currently available for melon are limited.

RESULT

We constructed eleven full-length enriched and four standard cDNA libraries from fruits, flowers, leaves, roots, cotyledons, and calluses of four different melon genotypes, and generated 71,577 and 22,179 ESTs from full-length enriched and standard cDNA libraries, respectively. These ESTs, together with ~35,000 ESTs available in public domains, were assembled into 24,444 unigenes, which were extensively annotated by comparing their sequences to different protein and functional domain databases, assigning them Gene Ontology (GO) terms, and mapping them onto metabolic pathways. Comparative analysis of melon unigenes and other plant genomes revealed that 75% to 85% of melon unigenes had homologs in other dicot plants, while approximately 70% had homologs in monocot plants. The analysis also identified 6,972 gene families that were conserved across dicot and monocot plants, and 181, 1,192, and 220 gene families specific to fleshy fruit-bearing plants, the Cucurbitaceae family, and melon, respectively. Digital expression analysis identified a total of 175 tissue-specific genes, which provides a valuable gene sequence resource for future genomics and functional studies. Furthermore, we identified 4,068 simple sequence repeats (SSRs) and 3,073 single nucleotide polymorphisms (SNPs) in the melon EST collection. Finally, we obtained a total of 1,382 melon full-length transcripts through the analysis of full-length enriched cDNA clones that were sequenced from both ends. Analysis of these full-length transcripts indicated that sizes of melon 5' and 3' UTRs were similar to those of tomato, but longer than many other dicot plants. Codon usages of melon full-length transcripts were largely similar to those of Arabidopsis coding sequences.

CONCLUSION

The collection of melon ESTs generated from full-length enriched and standard cDNA libraries is expected to play significant roles in annotating the melon genome. The ESTs and associated analysis results will be useful resources for gene discovery, functional analysis, marker-assisted breeding of melon and closely related species, comparative genomic studies and for gaining insights into gene expression patterns.

Distinct colonization patterns and cDNA-AFLP transcriptome profiles in compatible and incompatible interactions between melon and different races of Fusarium oxysporum f. sp. melonis

BACKGROUND

Fusarium oxysporum f. sp. melonis Snyd. & Hans. (FOM) causes Fusarium wilt, the most important infectious disease of melon (Cucumis melo L.). The four known races of this pathogen can be distinguished only by infection on appropriate cultivars. No molecular tools are available that can discriminate among the races, and the molecular basis of compatibility and disease progression are poorly understood. Resistance to races 1 and 2 is controlled by a single dominant gene, whereas only partial polygenic resistance to race 1,2 has been described. We carried out a large-scale cDNA-AFLP analysis to identify host genes potentially related to resistance and susceptibility as well as fungal genes associated with the infection process. At the same time, a systematic reisolation procedure on infected stems allowed us to monitor fungal colonization in compatible and incompatible host-pathogen combinations.

RESULTS

Melon plants (cv. Charentais Fom-2), which are susceptible to race 1,2 and resistant to race 1, were artificially infected with a race 1 strain of FOM or one of two race 1,2 w strains. Host colonization of stems was assessed at 1, 2, 4, 8, 14, 16, 18 and 21 days post inoculation (dpi), and the fungus was reisolated from infected plants. Markedly different colonization patterns were observed in compatible and incompatible host-pathogen combinations. Five time points from the symptomless early stage (2 dpi) to obvious wilting symptoms (21 dpi) were considered for cDNA-AFLP analysis. After successful sequencing of 627 transcript-derived fragments (TDFs) differentially expressed in infected plants, homology searching retrieved 305 melon transcripts, 195 FOM transcripts expressed in planta and 127 orphan TDFs. RNA samples from FOM colonies of the three strains grown in vitro were also included in the analysis to facilitate the detection of in planta-specific transcripts and to identify TDFs differentially expressed among races/strains.

CONCLUSION

Our data suggest that resistance against FOM in melon involves only limited transcriptional changes, and that wilting symptoms could derive, at least partially, from an active plant response.We discuss the pathogen-derived transcripts expressed in planta during the infection process and potentially related to virulence functions, as well as transcripts that are differentially expressed between the two FOM races grown in vitro. These transcripts provide candidate sequences that can be further tested for their ability to distinguish between races.Sequence data from this article have been deposited in GenBank, Accession Numbers: HO867279-HO867981.

Transcriptome characterization and high throughput SSRs and SNPs discovery in Cucurbita pepo (Cucurbitaceae)

Background

Cucurbita pepo belongs to the Cucurbitaceae family. The "Zucchini" types rank among the highest-valued vegetables worldwide, and other C. pepo and related Cucurbita spp., are food staples and rich sources of fat and vitamins. A broad range of genomic tools are today available for other cucurbits that have become models for the study of different metabolic processes. However, these tools are still lacking in the Cucurbita genus, thus limiting gene discovery and the process of breeding.

Results

We report the generation of a total of 512,751 C. pepo EST sequences, using 454 GS FLX Titanium technology. ESTs were obtained from normalized cDNA libraries (root, leaves, and flower tissue) prepared using two varieties with contrasting phenotypes for plant, flowering and fruit traits, representing the two C. pepo subspecies: subsp. pepo cv. Zucchini and subsp. ovifera cv Scallop. De novo assembling was performed to generate a collection of 49,610 Cucurbita unigenes (average length of 626 bp) that represent the first transcriptome of the species. Over 60% of the unigenes were functionally annotated and assigned to one or more Gene Ontology terms. The distributions of Cucurbita unigenes followed similar tendencies than that reported for Arabidopsis or melon, suggesting that the dataset may represent the whole Cucurbita transcriptome. About 34% unigenes were detected to have known orthologs of Arabidopsis or melon, including genes potentially involved in disease resistance, flowering and fruit quality. Furthermore, a set of 1,882 unigenes with SSR motifs and 9,043 high confidence SNPs between Zucchini and Scallop were identified, of which 3,538 SNPs met criteria for use with high throughput genotyping platforms, and 144 could be detected as CAPS. A set of markers were validated, being 80% of them polymorphic in a set of variable C. pepo and C. moschata accessions.

Conclusion

We present the first broad survey of gene sequences and allelic variation in C. pepo, where limited prior genomic information existed. The transcriptome provides an invaluable new tool for biological research. The developed molecular markers are the basis for future genetic linkage and quantitative trait loci analysis, and will be essential to speed up the process of breeding new and better adapted squash varieties.

Engineering melon plants with improved fruit shelf life using the TILLING approach

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.

Sequencing of 6.7 Mb of the melon genome using a BAC pooling strategy

BACKGROUND

Cucumis melo (melon) belongs to the Cucurbitaceae family, whose economic importance among horticulture crops is second only to Solanaceae. Melon has a high intra-specific genetic variation, morphologic diversity and a small genome size (454 Mb), which make it suitable for a great variety of molecular and genetic studies. A number of genetic and genomic resources have already been developed, such as several genetic maps, BAC genomic libraries, a BAC-based physical map and EST collections. Sequence information would be invaluable to complete the picture of the melon genomic landscape, furthering our understanding of this species' evolution from its relatives and providing an important genetic tool. However, to this day there is little sequence data available, only a few melon genes and genomic regions are deposited in public databases. The development of massively parallel sequencing methods allows envisaging new strategies to obtain long fragments of genomic sequence at higher speed and lower cost than previous Sanger-based methods.

RESULTS

In order to gain insight into the structure of a significant portion of the melon genome we set out to perform massive sequencing of pools of BAC clones. For this, a set of 57 BAC clones from a double haploid line was sequenced in two pools with the 454 system using both shotgun and paired-end approaches. The final assembly consists of an estimated 95% of the actual size of the melon BAC clones, with most likely complete sequences for 50 of the BACs, and a total sequence coverage of 39x. The accuracy of the assembly was assessed by comparing the previously available Sanger sequence of one of the BACs against its 454 sequence, and the polymorphisms found involved only 1.7 differences every 10,000 bp that were localized in 15 homopolymeric regions and two dinucleotide tandem repeats. Overall, the study provides approximately 6.7 Mb or 1.5% of the melon genome. The analysis of this new data has allowed us to gain further insight into characteristics of the melon genome such as gene density, average protein length, or microsatellite and transposon content. The annotation of the BAC sequences revealed a high degree of collinearity and protein sequence identity between melon and its close relative Cucumis sativus (cucumber). Transposon content analysis of the syntenic regions suggests that transposition activity after the split of both cucurbit species has been low in cucumber but very high in melon.

CONCLUSIONS

The results presented here show that the strategy followed, which combines shotgun and BAC-end sequencing together with anchored marker information, is an excellent method for sequencing specific genomic regions, especially from relatively compact genomes such as that of melon. However, in agreement with other results, this map-based, BAC approach is confirmed to be an expensive way of sequencing a whole plant genome. Our results also provide a partial description of the melon genome's structure. Namely, our analysis shows that the melon genome is highly collinear with the smaller one of cucumber, the size difference being mainly due to the expansion of intergenic regions and proliferation of transposable elements.

Generation and analysis of ESTs from strawberry (Fragaria xananassa) fruits and evaluation of their utility in genetic and molecular studies

BACKGROUND

Cultivated strawberry is a hybrid octoploid species (Fragaria xananassa Duchesne ex. Rozier) whose fruit is highly appreciated due to its organoleptic properties and health benefits. Despite recent studies on the control of its growth and ripening processes, information about the role played by different hormones on these processes remains elusive. Further advancement of this knowledge is hampered by the limited sequence information on genes from this species, despite the abundant information available on genes from the wild diploid relative Fragaria vesca. However, the diploid species, or one ancestor, only partially contributes to the genome of the cultivated octoploid. We have produced a collection of expressed sequence tags (ESTs) from different cDNA libraries prepared from different fruit parts and developmental stages. The collection has been analysed and the sequence information used to explore the involvement of different hormones in fruit developmental processes, and for the comparison of transcripts in the receptacle of ripe fruits of diploid and octoploid species. The study is particularly important since the commercial fruit is indeed an enlarged flower receptacle with the true fruits, the achenes, on the surface and connected through a network of vascular vessels to the central pith.

RESULTS

We have sequenced over 4,500 ESTs from Fragaria xananassa, thus doubling the number of ESTs available in the GenBank of this species. We then assembled this information together with that available from F. xananassa resulting a total of 7,096 unigenes. The identification of SSRs and SNPs in many of the ESTs allowed their conversion into functional molecular markers. The availability of libraries prepared from green growing fruits has allowed the cloning of cDNAs encoding for genes of auxin, ethylene and brassinosteroid signalling processes, followed by expression studies in selected fruit parts and developmental stages. In addition, the sequence information generated in the project, jointly with previous information on sequences from both F. xananassa and F. vesca, has allowed designing an oligo-based microarray that has been used to compare the transcriptome of the ripe receptacle of the diploid and octoploid species. Comparison of the transcriptomes, grouping the genes by biological processes, points to differences being quantitative rather than qualitative.

CONCLUSIONS

The present study generates essential knowledge and molecular tools that will be useful in improving investigations at the molecular level in cultivated strawberry (F. xananassa). This knowledge is likely to provide useful resources in the ongoing breeding programs. The sequence information has already allowed the development of molecular markers that have been applied to germplasm characterization and could be eventually used in QTL analysis. Massive transcription analysis can be of utility to target specific genes to be further studied, by their involvement in the different plant developmental processes.

Generation of a BAC-based physical map of the melon genome

BACKGROUND

Cucumis melo (melon) belongs to the Cucurbitaceae family, whose economic importance among horticulture crops is second only to Solanaceae. Melon has high intra-specific genetic variation, morphologic diversity and a small genome size (450 Mb), which make this species suitable for a great variety of molecular and genetic studies that can lead to the development of tools for breeding varieties of the species. A number of genetic and genomic resources have already been developed, such as several genetic maps and BAC genomic libraries. These tools are essential for the construction of a physical map, a valuable resource for map-based cloning, comparative genomics and assembly of whole genome sequencing data. However, no physical map of any Cucurbitaceae has yet been developed. A project has recently been started to sequence the complete melon genome following a whole-genome shotgun strategy, which makes use of massive sequencing data. A BAC-based melon physical map will be a useful tool to help assemble and refine the draft genome data that is being produced.

RESULTS

A melon physical map was constructed using a 5.7 x BAC library and a genetic map previously developed in our laboratories. High-information-content fingerprinting (HICF) was carried out on 23,040 BAC clones, digesting with five restriction enzymes and SNaPshot labeling, followed by contig assembly with FPC software. The physical map has 1,355 contigs and 441 singletons, with an estimated physical length of 407 Mb (0.9 x coverage of the genome) and the longest contig being 3.2 Mb. The anchoring of 845 BAC clones to 178 genetic markers (100 RFLPs, 76 SNPs and 2 SSRs) also allowed the genetic positioning of 183 physical map contigs/singletons, representing 55 Mb (12%) of the melon genome, to individual chromosomal loci. The melon FPC database is available for download at http://melonomics.upv.es/static/files/public/physical_map/.

CONCLUSIONS

Here we report the construction of the first physical map of a Cucurbitaceae species described so far. The physical map was integrated with the genetic map so that a number of physical contigs, representing 12% of the melon genome, could be anchored to known genetic positions. The data presented is already helping to improve the quality of the melon genomic sequence available as a result of a project currently being carried out in Spain, adopting a whole genome shotgun approach based on 454 sequencing data.

Computational annotation of genes differentially expressed along olive fruit development

BACKGROUND

Olea europaea L. is a traditional tree crop of the Mediterranean basin with a worldwide economical high impact. Differently from other fruit tree species, little is known about the physiological and molecular basis of the olive fruit development and a few sequences of genes and gene products are available for olive in public databases. This study deals with the identification of large sets of differentially expressed genes in developing olive fruits and the subsequent computational annotation by means of different software.

RESULTS

mRNA from fruits of the cv. Leccino sampled at three different stages [i.e., initial fruit set (stage 1), completed pit hardening (stage 2) and veraison (stage 3)] was used for the identification of differentially expressed genes putatively involved in main processes along fruit development. Four subtractive hybridization libraries were constructed: forward and reverse between stage 1 and 2 (libraries A and B), and 2 and 3 (libraries C and D). All sequenced clones (1,132 in total) were analyzed through BlastX against non-redundant NCBI databases and about 60% of them showed similarity to known proteins. A total of 89 out of 642 differentially expressed unique sequences was further investigated by Real-Time PCR, showing a validation of the SSH results as high as 69%. Library-specific cDNA repertories were annotated according to the three main vocabularies of the gene ontology (GO): cellular component, biological process and molecular function. BlastX analysis, GO terms mapping and annotation analysis were performed using the Blast2GO software, a research tool designed with the main purpose of enabling GO based data mining on sequence sets for which no GO annotation is yet available. Bioinformatic analysis pointed out a significantly different distribution of the annotated sequences for each GO category, when comparing the three fruit developmental stages. The olive fruit-specific transcriptome dataset was used to query all known KEGG (Kyoto Encyclopaedia of Genes and Genomes) metabolic pathways for characterizing and positioning retrieved EST records. The integration of the olive sequence datasets within the MapMan platform for microarray analysis allowed the identification of specific biosynthetic pathways useful for the definition of key functional categories in time course analyses for gene groups.

CONCLUSION

The bioinformatic annotation of all gene sequences was useful to shed light on metabolic pathways and transcriptional aspects related to carbohydrates, fatty acids, secondary metabolites, transcription factors and hormones as well as response to biotic and abiotic stresses throughout olive drupe development. These results represent a first step toward both functional genomics and systems biology research for understanding the gene functions and regulatory networks in olive fruit growth and ripening.

An oligo-based microarray offers novel transcriptomic approaches for the analysis of pathogen resistance and fruit quality traits in melon (Cucumis melo L.)

Background

Melon (Cucumis melo) is a horticultural specie of significant nutritional value, which belongs to the Cucurbitaceae family, whose economic importance is second only to the Solanaceae. Its small genome of approx. 450 Mb coupled to the high genetic diversity has prompted the development of genetic tools in the last decade. However, the unprecedented existence of a transcriptomic approaches in melon, highlight the importance of designing new tools for high-throughput analysis of gene expression.

Results

We report the construction of an oligo-based microarray using a total of 17,510 unigenes derived from 33,418 high-quality melon ESTs. This chip is particularly enriched with genes that are expressed in fruit and during interaction with pathogens. Hybridizations for three independent experiments allowed the characterization of global gene expression profiles during fruit ripening, as well as in response to viral and fungal infections in plant cotyledons and roots, respectively. Microarray construction, statistical analyses and validation together with functional-enrichment analysis are presented in this study.

Conclusion

The platform validation and enrichment analyses shown in our study indicate that this oligo-based microarray is amenable for future genetic and functional genomic studies of a wide range of experimental conditions in melon.

A set of EST-SNPs for map saturation and cultivar identification in melon

BACKGROUND

There are few genomic tools available in melon (Cucumis melo L.), a member of the Cucurbitaceae, despite its importance as a crop. Among these tools, genetic maps have been constructed mainly using marker types such as simple sequence repeats (SSR), restriction fragment length polymorphisms (RFLP) and amplified fragment length polymorphisms (AFLP) in different mapping populations. There is a growing need for saturating the genetic map with single nucleotide polymorphisms (SNP), more amenable for high throughput analysis, especially if these markers are located in gene coding regions, to provide functional markers. Expressed sequence tags (ESTs) from melon are available in public databases, and resequencing ESTs or validating SNPs detected in silico are excellent ways to discover SNPs.

RESULTS

EST-based SNPs were discovered after resequencing ESTs between the parental lines of the PI 161375 (SC) x 'Piel de sapo' (PS) genetic map or using in silico SNP information from EST databases. In total 200 EST-based SNPs were mapped in the melon genetic map using a bin-mapping strategy, increasing the map density to 2.35 cM/marker. A subset of 45 SNPs was used to study variation in a panel of 48 melon accessions covering a wide range of the genetic diversity of the species. SNP analysis correctly reflected the genetic relationships compared with other marker systems, being able to distinguish all the accessions and cultivars.

CONCLUSION

This is the first example of a genetic map in a cucurbit species that includes a major set of SNP markers discovered using ESTs. The PI 161375 x 'Piel de sapo' melon genetic map has around 700 markers, of which more than 500 are gene-based markers (SNP, RFLP and SSR). This genetic map will be a central tool for the construction of the melon physical map, the step prior to sequencing the complete genome. Using the set of SNP markers, it was possible to define the genetic relationships within a collection of forty-eight melon accessions as efficiently as with SSR markers, and these markers may also be useful for cultivar identification in Occidental melon varieties.

Dissection of the oligogenic resistance to Cucumber mosaic virus in the melon accession PI 161375

Resistance to Cucumber mosaic virus (CMV) in the exotic melon accession PI 161375, cultivar "Sonwang Charmi" (SC) had previously been described as oligogenic, recessive and quantitative, with a major QTL residing in linkage group XII (LGXII). We have used a collection of near isogenic lines (NILs) with introgressions of SC into the genome of the susceptible accession Piel de Sapo (PS) to further characterise this resistance. Infection of NILs carrying introgressions on LGXII showed that only NIL SC12-1 was resistant to CMV strains P9 and P104.82, but not to strains M6 and TL. Further mapping of this region showed that the resistance, named cmv1 maps in an area of 2.2 cM, between markers CMN61_44 and CMN21_55. Moreover, cmv1 confers total resistance to strains P9 and P104.82, indicating that in these cases it is not quantitative and that cmv1 is sufficient to confer full resistance to these CMV strains. Candidate gene mapping of ten translation initiation factors in the melon genome failed to find any of them in the interval between markers CMN61_44 and CMN21_55. All these results suggest that the resistance to CMV present in SC is oligogenic, where different loci confer resistance to different CMV strains, but not necessarily quantitative, since at least one of these genes (cmv1) confers total resistance, similar to that of the parental SC, and does not need the contribution of other loci.

Identification of QTLs for resistance to powdery mildew and SSR markers diagnostic for powdery mildew resistance genes in melon (Cucumis melo L.)

Powdery mildew caused by Podosphaera xanthii is an important foliar disease in melon. To find molecular markers for marker-assisted selection, we constructed a genetic linkage map of melon based on a population of 93 recombinant inbred lines derived from crosses between highly resistant AR 5 and susceptible 'Earl's Favourite (Harukei 3)'. The map spans 877 cM and consists of 167 markers, comprising 157 simple sequence repeats (SSRs), 7 sequence characterized amplified region/cleavage amplified polymorphic sequence markers and 3 phenotypic markers segregating into 20 linkage groups. Among them, 37 SSRs and 6 other markers were common to previous maps. Quantitative trait locus (QTL) analysis identified two loci for resistance to powdery mildew. The effects of these QTLs varied depending on strain and plant stage. The percentage of phenotypic variance explained for resistance to the pxA strain was similar between QTLs (R (2) = 22-28%). For resistance to pxB strain, the QTL on linkage group (LG) XII was responsible for much more of the variance (41-46%) than that on LG IIA (12-13%). The QTL on LG IIA was located between two SSR markers. Using an independent population, we demonstrated the effectiveness of these markers. This is the first report of universal and effective markers linked to a gene for powdery mildew resistance in melon.