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

Rapid sequence and functional diversification of a miRNA superfamily targeting calcium signaling components in seed plants

MicroRNA (miRNA)-directed posttranscriptional gene silencing (miR-PTGS) is an integral component of gene regulatory networks governing plant development and responses to the environment. The sequence homology between Sly-miR4376, a miRNA common to Solanaceae and reported to target autoinhibited Ca²⁺ -ATPase 10 (ACA10) messenger RNA (mRNA) in tomato, and Arabidopsis miR391 (Ath-miR391), previously annotated as a nonconserved member of the deeply conserved miR390 family, has prompted us to revisit the function of Ath-miR391, as well as its regulatory conservation. A combination of genetic, molecular, and bioinformatic analyses revealed a hidden conservation for miR-PTGS of ACA10 homologs in spermatophytes. We found that the Arabidopsis ACA10 mRNA undergoes miR391-directed cleavage in vivo. Furthermore, transgenic overexpression of miR391 recapitulated the compact inflorescence (cif) phenotypes characteristic of ACA10 loss-of-function mutants, due to miR391-directed PTGS of ACA10. Significantly, comprehensive data mining revealed robust evidence for widespread PTGS of ACA10 homologs directed by a superfamily of related miRNAs sharing a conserved sequence core. Intriguingly, the ACA-targeting miRNAs in Poaceae also direct PTGS for calmodulin-like proteins which are putative Ca²⁺ sensors. The PTGS of ACA10 homologs is therefore directed by a miRNA superfamily that is of ancient origin and has undergone rapid sequence diversification associated with functional innovation.

Starch branching enzymes as putative determinants of postharvest quality in horticultural crops

Starch branching enzymes (SBEs) are key determinants of the structure and amount of the starch in plant organs, and as such, they have the capacity to influence plant growth, developmental, and fitness processes, and in addition, the industrial end-use of starch. However, little is known about the role of SBEs in determining starch structure-function relations in economically important horticultural crops such as fruit and leafy greens, many of which accumulate starch transiently. Further, a full understanding of the biological function of these types of starches is lacking. Because of this gap in knowledge, this minireview aims to provide an overview of SBEs in horticultural crops, to investigate the potential role of starch in determining postharvest quality. A systematic examination of SBE sequences in 43 diverse horticultural species, identified SBE1, 2 and 3 isoforms in all species examined except apple, olive, and Brassicaceae, which lacked SBE1, but had a duplicated SBE2. Among our findings after a comprehensive and critical review of published data, was that as apple, banana, and tomato fruits ripens, the ratio of the highly digestible amylopectin component of starch increases relative to the more digestion-resistant amylose fraction, with parallel increases in SBE2 transcription, fruit sugar content, and decreases in starch. It is tempting to speculate that during the ripening of these fruit when starch degradation occurs, there are rearrangements made to the structure of starch possibly via branching enzymes to increase starch digestibility to sugars. We propose that based on the known action of SBEs, and these observations, SBEs may affect produce quality, and shelf-life directly through starch accumulation, and indirectly, by altering sugar availability. Further studies where SBE activity is fine-tuned in these crops, can enrich our understanding of the role of starch across species and may improve horticulture postharvest quality.

Identification and Characterization of Stress-Responsive TAS3-Derived TasiRNAs in Melon

Small interfering RNAs (siRNA) are key regulators of gene expression that play essential roles in diverse biological processes. Trans-acting siRNAs (tasiRNAs) are a class of plant-endogenous siRNAs that lead the cleavage of nonidentical transcripts. TasiRNAs are usually involved in fine-tuning development. However, increasing evidence supports that tasiRNAs may be involved in stress response. Melon is a crop of great economic importance extensively cultivated in semiarid regions frequently exposed to changing environmental conditions that limit its productivity. However, knowledge of the precise role of siRNAs in general, and of tasiRNAs in particular, in regulating the response to adverse environmental conditions is limited. Here, we provide the first comprehensive analysis of computationally inferred melon-tasiRNAs responsive to two biotic (viroid-infection) and abiotic (cold treatment) stress conditions. We identify two TAS3-loci encoding to length (TAS3-L) and short (TAS3-S) transcripts. The TAS candidates predicted from small RNA-sequencing data were characterized according to their chromosome localization and expression pattern in response to stress. The functional activity of cmTAS genes was validated by transcript quantification and degradome assays of the tasiRNA precursors and their predicted targets. Finally, the functionality of a representative cmTAS3-derived tasiRNA (TAS3-S) was confirmed by transient assays showing the cleavage of ARF target transcripts.

Inferring the regulatory network of the miRNA-mediated response to biotic and abiotic stress in melon

BACKGROUND

MiRNAs have emerged as key regulators of stress response in plants, suggesting their potential as candidates for knock-in/out to improve stress tolerance in agricultural crops. Although diverse assays have been performed, systematic and detailed studies of miRNA expression and function during exposure to multiple environments in crops are limited.

RESULTS

Here, we present such pioneering analysis in melon plants in response to seven biotic and abiotic stress conditions. Deep-sequencing and computational approaches have identified twenty-four known miRNAs whose expression was significantly altered under at least one stress condition, observing that down-regulation was preponderant. Additionally, miRNA function was characterized by high scale degradome assays and quantitative RNA measurements over the intended target mRNAs, providing mechanistic insight. Clustering analysis provided evidence that eight miRNAs showed a broad response range under the stress conditions analyzed, whereas another eight miRNAs displayed a narrow response range. Transcription factors were predominantly targeted by stress-responsive miRNAs in melon. Furthermore, our results show that the miRNAs that are down-regulated upon stress predominantly have as targets genes that are known to participate in the stress response by the plant, whereas the miRNAs that are up-regulated control genes linked to development.

CONCLUSION

Altogether, this high-resolution analysis of miRNA-target interactions, combining experimental and computational work, Illustrates the close interplay between miRNAs and the response to diverse environmental conditions, in melon.

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.

Haemodynamics Regulate Fibronectin Assembly via PECAM

Fibronectin (FN) assembly and fibrillogenesis are critically important in both development and the adult organism, but their importance in vascular functions is not fully understood. Here we identify a novel pathway by which haemodynamic forces regulate FN assembly and fibrillogenesis during vascular remodelling. Induction of disturbed shear stress in vivo and in vitro resulted in complex FN fibril assembly that was dependent on the mechanosensor PECAM. Loss of PECAM also inhibited the cell-intrinsic ability to remodel FN. Gain- and loss-of-function experiments revealed that PECAM-dependent RhoA activation is required for FN assembly. Furthermore, PECAM^-/- mice exhibited reduced levels of active β1 integrin that were responsible for reduced RhoA activation and downstream FN assembly. These data identify a new pathway by which endothelial mechanotransduction regulates FN assembly and flow-mediated vascular remodelling.

Whole Genome Re-Sequencing and Characterization of Powdery Mildew Disease-Associated Allelic Variation in Melon

Powdery mildew is one of the most common fungal diseases in the world. This disease frequently affects melon (Cucumis melo L.) and other Cucurbitaceous family crops in both open field and greenhouse cultivation. One of the goals of genomics is to identify the polymorphic loci responsible for variation in phenotypic traits. In this study, powdery mildew disease assessment scores were calculated for four melon accessions, 'SCNU1154', 'Edisto47', 'MR-1', and 'PMR5'. To investigate the genetic variation of these accessions, whole genome re-sequencing using the Illumina HiSeq 2000 platform was performed. A total of 754,759,704 quality-filtered reads were generated, with an average of 82.64% coverage relative to the reference genome. Comparisons of the sequences for the melon accessions revealed around 7.4 million single nucleotide polymorphisms (SNPs), 1.9 million InDels, and 182,398 putative structural variations (SVs). Functional enrichment analysis of detected variations classified them into biological process, cellular component and molecular function categories. Further, a disease-associated QTL map was constructed for 390 SNPs and 45 InDels identified as related to defense-response genes. Among them 112 SNPs and 12 InDels were observed in powdery mildew responsive chromosomes. Accordingly, this whole genome re-sequencing study identified SNPs and InDels associated with defense genes that will serve as candidate polymorphisms in the search for sources of resistance against powdery mildew disease and could accelerate marker-assisted breeding in melon.

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.

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.

Combined use of genetic and genomics resources to understand virus resistance and fruit quality traits in melon

The availability of the genome sequence of many crop species during the past few years has opened a new era in plant biology, allowing for the performance of massive genomic studies in plant species other than the classical models Arabidopsis and rice. One of these crop species is melon (Cucumis melo), a cucurbit of high economic value that has become an interesting model for the study of biological processes such as fruit ripening, sex determination and phloem transport. The recent availability of the melon genome sequence, together with a number of genetic and genomic resources, provides powerful tools that can be used to assist in the main melon breeding targets, namely disease resistance and fruit quality. In this review, we will describe recent data obtained combining the use of a melon near isogenic line (NIL) population and genomic resources to gain insight into agronomically important traits as fruit ripening, resistance to Cucumber Mosaic virus (CMV) and the accumulation of sugars in fruits.

Characterization of the cork oak transcriptome dynamics during acorn development

BACKGROUND

Cork oak (Quercus suber L.) has a natural distribution across western Mediterranean regions and is a keystone forest tree species in these ecosystems. The fruiting phase is especially critical for its regeneration but the molecular mechanisms underlying the biochemical and physiological changes during cork oak acorn development are poorly understood. In this study, the transcriptome of the cork oak acorn, including the seed, was characterized in five stages of development, from early development to acorn maturation, to identify the dominant processes in each stage and reveal transcripts with important functions in gene expression regulation and response to water.

RESULTS

A total of 80,357 expressed sequence tags (ESTs) were de novo assembled from RNA-Seq libraries representative of the several acorn developmental stages. Approximately 7.6 % of the total number of transcripts present in Q. suber transcriptome was identified as acorn specific. The analysis of expression profiles during development returned 2,285 differentially expressed (DE) transcripts, which were clustered into six groups. The stage of development corresponding to the mature acorn exhibited an expression profile markedly different from other stages. Approximately 22 % of the DE transcripts putatively code for transcription factors (TF) or transcriptional regulators, and were found almost equally distributed among the several expression profile clusters, highlighting their major roles in controlling the whole developmental process. On the other hand, carbohydrate metabolism, the biological pathway most represented during acorn development, was especially prevalent in mid to late stages as evidenced by enrichment analysis. We further show that genes related to response to water, water deprivation and transport were mostly represented during the early (S2) and the last stage (S8) of acorn development, when tolerance to water desiccation is possibly critical for acorn viability.

CONCLUSIONS

To our knowledge this work represents the first report of acorn development transcriptomics in oaks. The obtained results provide novel insights into the developmental biology of cork oak acorns, highlighting transcripts putatively involved in the regulation of the gene expression program and in specific processes likely essential for adaptation. It is expected that this knowledge can be transferred to other oak species of great ecological value.

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.

Systems approach for exploring the intricate associations between sweetness, color and aroma in melon fruits

BACKGROUND

Melon (Cucumis melo) fruits exhibit phenotypic diversity in several key quality determinants such as taste, color and aroma. Sucrose, carotenoids and volatiles are recognized as the key compounds shaping the above corresponding traits yet the full network of biochemical events underlying their synthesis have not been comprehensively described. To delineate the cellular processes shaping fruit quality phenotypes, a population of recombinant inbred lines (RIL) was used as a source of phenotypic and genotypic variations. In parallel, ripe fruits were analyzed for both the quantified level of 77 metabolic traits directly associated with fruit quality and for RNA-seq based expression profiles generated for 27,000 unigenes. First, we explored inter-metabolite association patterns; then, we described metabolites versus gene association patterns; finally, we used the correlation-based associations for predicting uncharacterized synthesis pathways.

RESULTS

Based on metabolite versus metabolite and metabolite versus gene association patterns, we divided metabolites into two key groups: a group including ethylene and aroma determining volatiles whose accumulation patterns are correlated with the expression of genes involved in the glycolysis and TCA cycle pathways; and a group including sucrose and color determining carotenoids whose accumulation levels are correlated with the expression of genes associated with plastid formation.

CONCLUSIONS

The study integrates multiple processes into a genome scale perspective of cellular activity. This lays a foundation for deciphering the role of gene markers associated with the determination of fruit quality traits.

An acorn squash (Cucurbita pepo ssp. ovifera) fruit and seed transcriptome as a resource for the study of fruit traits in Cucurbita

Acorn squash (Cucurbita pepo) is an iconic fall vegetable in the United States, known for its unique fruit shape and also prized for its culinary properties. Little is known about the metabolism that underlies the development of fruit quality attributes such as color, sweetness, texture and nutritional qualities in acorn squash, or any other winter squash grown worldwide. To provide insight into winter squash fruit and seed development and add to the genomic resources in the Cucurbita genus, RNA sequencing was used to generate an acorn squash fruit and seed transcriptome from the cultivar Sweet REBA at critical points throughout fruit development. 141 838 600 high-quality paired-end Illumina reads were assembled into 55 949 unigenes. 85% of unigenes with predicted open reading frames had homology with previously identified genes and over 62% could be functionally annotated. Comparison with the watermelon and cucumber genomes provided confirmation that the unigenes are full-length and comprehensive, covering an average of 90% of the coding sequence of their homologs and 72% of the cucumber and watermelon exomes. Key candidate genes associated with carotenoid and carbohydrate metabolism were identified toward a resource for winter squash fruit quality trait dissection. This transcriptome represents a major advance in C. pepo genomics, providing significant new sequence information and revealing the repertoire of genes expressed throughout winter squash fruit and seed development. Future studies on the genetic basis of fruit quality and future breeding efforts will be enhanced by tools and insights developed from this resource.

Use of targeted SNP selection for an improved anchoring of the melon (Cucumis melo L.) scaffold genome assembly

Background

The genome of the melon (Cucumis melo L.) double-haploid line DHL92 was recently sequenced, with 87.5 and 80.8% of the scaffold assembly anchored and oriented to the 12 linkage groups, respectively. However, insufficient marker coverage and a lack of recombination left several large, gene rich scaffolds unanchored, and some anchored scaffolds unoriented. To improve the anchoring and orientation of the melon genome assembly, we used resequencing data between the parental lines of DHL92 to develop a new set of SNP markers from unanchored scaffolds.

Results

A high-resolution genetic map composed of 580 SNPs was used to anchor 354.8 Mb of sequence, contained in 141 scaffolds (average size 2.5 Mb) and corresponding to 98.2% of the scaffold assembly, to the 12 melon chromosomes. Over 325.4 Mb (90%) of the assembly was oriented. The genetic map revealed regions of segregation distortion favoring SC alleles as well as recombination suppression regions coinciding with putative centromere, 45S, and 5S rDNA sites. New chromosome-scale pseudomolecules were created by incorporating to the previous v3.5 version an additional 38.3 Mb of anchored sequence representing 1,837 predicted genes contained in 55 scaffolds. Using fluorescent in situ hybridization (FISH) with BACs that produced chromosome-specific signals, melon chromosomes that correspond to the twelve linkage groups were identified, and a standardized karyotype of melon inbred line T111 was developed.

Conclusions

By utilizing resequencing data and targeted SNP selection combined with a large F2 mapping population, we significantly improved the quantity of anchored and oriented melon scaffold genome assembly. Using genome information combined with FISH mapping provided the first cytogenetic map of an inodorus melon type. With these results it was possible to make inferences on melon chromosome structure by relating zones of recombination suppression to centromeres and 45S and 5S heterochromatic regions. This study represents the first steps towards the integration of the high-resolution genetic and cytogenetic maps with the genomic sequence in melon that will provide more information on genome organization and allow for the improvement of the melon genome draft sequence.

Electronic supplementary material

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

Mapping of the gynoecy in bitter gourd (Momordica charantia) using RAD-seq analysis

Momordica charantia is a monoecious plant of the Cucurbitaceae family that has both male and female unisexual flowers. Its unique gynoecious line, OHB61-5, is essential as a maternal parent in the production of F₁ cultivars. To identify the DNA markers for this gynoecy, a RAD-seq (restriction-associated DNA tag sequencing) analysis was employed to reveal genome-wide DNA polymorphisms and to genotype the F₂ progeny from a cross between OHB61-5 and a monoecious line. Based on a RAD-seq analysis of F₂ individuals, a linkage map was constructed using 552 co-dominant markers. In addition, after analyzing the pooled genomic DNA from monoecious or gynoecious F₂ plants, several SNP loci that are genetically linked to gynoecy were identified. GTFL-1, the closest SNP locus to the putative gynoecious locus, was converted to a conventional DNA marker using invader assay technology, which is applicable to the marker-assisted selection of gynoecy in M. charantia breeding.

Screening suitable reference genes for normalization in reverse transcription quantitative real-time PCR analysis in melon

Melon (Cucumis melo. L) is not only an economically important cucurbitaceous crop but also an attractive model for studying many biological characteristics. Screening appropriate reference genes is essential to reverse transcription quantitative real-time PCR (RT-qPCR), which is key to many studies involving gene expression analysis. In this study, 14 candidate reference genes were selected, and the variations in their expression in roots and leaves of plants subjected to biotic stress, abiotic stress, and plant growth regulator treatment were assessed by RT-qPCR. The stability of the expression of the selected genes was determined and ranked using geNorm and NormFinder. geNorm identified the two most stable genes for each set of conditions: CmADP and CmUBIep across all samples, CmUBIep and CmRPL in roots, CmRAN and CmACT in leaves, CmADP and CmRPL under abiotic stress conditions, CmTUA and CmACT under biotic stress conditions, and CmRAN and CmACT under plant growth regulator treatments. NormFinder determined CmRPL to be the best reference gene in roots and under biotic stress conditions and CmADP under the other experimental conditions. CmUBC2 and CmPP2A were not found to be suitable under many experimental conditions. The catalase family genes CmCAT1, CmCAT2, and CmCAT3 were identified in melon genome and used as target genes to validate the reliability of identified reference genes. The catalase family genes showed the most upregulation 3 days after inoculation with Fusarium wilt in roots, after which they were downregulated. Their levels of expression were significantly overestimated when the unsuitable reference gene was used for normalization. These results not only provide guidelines for the selection of reference genes for gene expression analyses in melons but may also provide valuable information for studying the functions of catalase family genes in stress responses.

Screening suitable reference genes for normalization in reverse transcription quantitative real-time PCR analysis in melon

Melon (Cucumis melo. L) is not only an economically important cucurbitaceous crop but also an attractive model for studying many biological characteristics. Screening appropriate reference genes is essential to reverse transcription quantitative real-time PCR (RT-qPCR), which is key to many studies involving gene expression analysis. In this study, 14 candidate reference genes were selected, and the variations in their expression in roots and leaves of plants subjected to biotic stress, abiotic stress, and plant growth regulator treatment were assessed by RT-qPCR. The stability of the expression of the selected genes was determined and ranked using geNorm and NormFinder. geNorm identified the two most stable genes for each set of conditions: CmADP and CmUBIep across all samples, CmUBIep and CmRPL in roots, CmRAN and CmACT in leaves, CmADP and CmRPL under abiotic stress conditions, CmTUA and CmACT under biotic stress conditions, and CmRAN and CmACT under plant growth regulator treatments. NormFinder determined CmRPL to be the best reference gene in roots and under biotic stress conditions and CmADP under the other experimental conditions. CmUBC2 and CmPP2A were not found to be suitable under many experimental conditions. The catalase family genes CmCAT1, CmCAT2, and CmCAT3 were identified in melon genome and used as target genes to validate the reliability of identified reference genes. The catalase family genes showed the most upregulation 3 days after inoculation with Fusarium wilt in roots, after which they were downregulated. Their levels of expression were significantly overestimated when the unsuitable reference gene was used for normalization. These results not only provide guidelines for the selection of reference genes for gene expression analyses in melons but may also provide valuable information for studying the functions of catalase family genes in stress responses.

Leaf-, panel- and latex-expressed sequenced tags from the rubber tree (Hevea brasiliensis) under cold-stressed and suboptimal growing conditions: the development of gene-targeted functional markers for stress response

Hevea brasiliensis is a native species of the Amazon Basin of South America and the primary source of natural rubber worldwide. Due to the occurrence of South American Leaf Blight disease in this area, rubber plantations have been extended to suboptimal regions. Rubber tree breeding is time-consuming and expensive, but molecular markers can serve as a tool for early evaluation, thus reducing time and costs. In this work, we constructed six different cDNA libraries with the aim of developing gene-targeted molecular markers for the rubber tree. A total of 8,263 reads were assembled, generating 5,025 unigenes that were analyzed; 912 expressed sequence tags (ESTs) represented new transcripts, and two sequences were highly up-regulated by cold stress. These unigenes were scanned for microsatellite (SSR) regions and single nucleotide polymorphisms (SNPs). In total, 169 novel EST-SSR markers were developed; 138 loci were polymorphic in the rubber tree, and 98 % presented transferability to six other Hevea species. Locus duplication was observed in H. brasiliensis and other species. Additionally, 43 SNP markers in 13 sequences that showed similarity to proteins involved in stress response, latex biosynthesis and developmental processes were characterized. cDNA libraries are a rich source of SSR and SNP markers and enable the identification of new transcripts. The new markers developed here will be a valuable resource for linkage mapping, QTL identification and other studies in the rubber tree and can also be used to evaluate the genetic variability of other Hevea species, which are valuable assets in rubber tree breeding.

Comprehensive analysis of expressed sequence tags from cultivated and wild radish (Raphanus spp.)

BACKGROUND

Radish (Raphanus sativus L., 2n = 2× = 18) is an economically important vegetable crop worldwide. A large collection of radish expressed sequence tags (ESTs) has been generated but remains largely uncharacterized.

RESULTS

In this study, approximately 315,000 ESTs derived from 22 Raphanus cDNA libraries from 18 different genotypes were analyzed, for the purpose of gene and marker discovery and to evaluate large-scale genome duplication and phylogenetic relationships among Raphanus spp. The ESTs were assembled into 85,083 unigenes, of which 90%, 65%, 89% and 89% had homologous sequences in the GenBank nr, SwissProt, TrEMBL and Arabidopsis protein databases, respectively. A total of 66,194 (78%) could be assigned at least one gene ontology (GO) term. Comparative analysis identified 5,595 gene families unique to radish that were significantly enriched with genes related to small molecule metabolism, as well as 12,899 specific to the Brassicaceae that were enriched with genes related to seed oil body biogenesis and responses to phytohormones. The analysis further indicated that the divergence of radish and Brassica rapa occurred approximately 8.9-14.9 million years ago (MYA), following a whole-genome duplication event (12.8-21.4 MYA) in their common ancestor. An additional whole-genome duplication event in radish occurred at 5.1-8.4 MYA, after its divergence from B. rapa. A total of 13,570 simple sequence repeats (SSRs) and 28,758 high-quality single nucleotide polymorphisms (SNPs) were also identified. Using a subset of SNPs, the phylogenetic relationships of eight different accessions of Raphanus was inferred.

CONCLUSION

Comprehensive analysis of radish ESTs provided new insights into radish genome evolution and the phylogenetic relationships of different radish accessions. Moreover, the radish EST sequences and the associated SSR and SNP markers described in this study represent a valuable resource for radish functional genomics studies and breeding.

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.

Transcriptomic events involved in melon mature-fruit abscission comprise the sequential induction of cell-wall degrading genes coupled to a stimulation of endo and exocytosis

BACKGROUND

Mature-fruit abscission (MFA) in fleshy-fruit is a genetically controlled process with mechanisms that, contrary to immature-fruit abscission, has not been fully characterized. Here, we use pyrosequencing to characterize the transcriptomes of melon abscission zone (AZ) at three stages during AZ-cell separation in order to understand MFA control at an early stage of AZ-activation.

PRINCIPAL FINDINGS

The results show that by early induction of MFA, the melon AZ exhibits major gene induction, while by late induction of MFA, melon AZ shows major gene repression. Although some genes displayed similar regulation in both early and late induction of abscission, such as EXT1-EXT4, EGase1, IAA2, ERF1, AP2D15, FLC, MADS2, ERAF17, SAP5 and SCL13 genes, the majority had different expression patterns. This implies that time-specific events occur during MFA, and emphasizes the value of characterizing multiple time-specific abscission transcriptomes. Analysis of gene-expression from these AZs reveal that a sequential induction of cell-wall-degrading genes is associated with the upregulation of genes involved in endo and exocytosis, and a shift in plant-hormone metabolism and signaling genes during MFA. This is accompanied by transcriptional activity of small-GTPases and synthaxins together with tubulins, dynamins, V-type ATPases and kinesin-like proteins potentially involved in MFA signaling. Early events are potentially controlled by down-regulation of MADS-box, AP2/ERF and Aux/IAA transcription-factors, and up-regulation of homeobox, zinc finger, bZIP, and WRKY transcription-factors, while late events may be controlled by up-regulation of MYB transcription-factors.

SIGNIFICANCE

Overall, the data provide a comprehensive view on MFA in fleshy-fruit, identifying candidate genes and pathways associated with early induction of MFA. Our comprehensive gene-expression profile will be very useful for elucidating gene regulatory networks of the MFA in fleshy-fruit.

Meta-analytical biomarker search of EST expression data reveals three differentially expressed candidates

Background

Researches have been conducted for the identification of differentially expressed genes (DEGs) by generating and mining of cDNA expressed sequence tags (ESTs) for more than a decade. Although the availability of public databases make possible the comprehensive mining of DEGs among the ESTs from multiple tissue types, existing studies usually employed statistics suitable only for two categories. Multi-class test has been developed to enable the finding of tissue specific genes, but subsequent search for cancer genes involves separate two-category test only on the ESTs of the tissue of interest. This constricts the amount of data used. On the other hand, simple pooling of cancer and normal genes from multiple tissue types runs the risk of Simpson's paradox. Here we presented a different approach which searched for multi-cancer DEG candidates by analyzing all pertinent ESTs in all categories and narrowing down the cancer biomarker candidates via integrative analysis with microarray data and selection of secretory and membrane protein genes as well as incorporation of network analysis. Finally, the differential expression patterns of three selected cancer biomarker candidates were confirmed by real-time qPCR analysis.

Results

Seven hundred and twenty three primary DEG candidates (p-value < 0.05 and lower bound of confidence interval of odds ratio ≧ 1.65) were selected from a curated EST database with the application of Cochran-Mantel-Haenszel statistic (CMH). GeneGO analysis results indicated this set as neoplasm enriched. Cross-examination with microarray data further narrowed the list down to 235 genes, among which 96 had membrane or secretory annotations. After examined the candidates in protein interaction network, public tissue expression databases, and literatures, we selected three genes for further evaluation by real-time qPCR with eight major normal and cancer tissues. The higher-than-normal tissue expression of COL3A1, DLG3, and RNF43 in some of the cancer tissues is in agreement with our in silico predictions.

Conclusions

Searching digitized transcriptome using CMH enabled us to identify multi-cancer differentially expressed gene candidates. Our methodology demonstrated simultaneously analysis for cancer biomarkers of multiple tissue types with the EST data. With the revived interest in digitizing the transcriptomes by NGS, cancer biomarkers could be more precisely detected from the ESTs. The three candidates identified in this study, COL3A1, DLG3, and RNF43, are valuable targets for further evaluation with a larger sample size of normal and cancer tissue or serum samples.

Inhibition of a basal transcription factor 3-like gene Osj10gBTF3 in rice results in significant plant miniaturization and typical pollen abortion

BTF3, which was originally recognized as a basal transcription factor, has been known to be involved in transcription initiation, translational regulation and protein localization in many eukaryotic organisms. However, its function remains largely unknown in plant species. In the present study, we analyzed a BTF3-related sequence in Oryza sativa L. subsp. japonica, which shares the conserved domain of a nascent polypeptide-associated complex with human BTF3, and was referred to as Osj10gBTF3. The expression of Osj10gBTF3 was primarily constitutive and generally modulated by salt, high temperature and exogenous phytohormone stress. The Osj10gBTF3::EGFP (enhanced green fluorescence protein) fusion protein was localized in both the nucleus and cytoplasmic membrane system. Inhibition of Osj10gBTF3 led to significant morphological changes in all detected tissues and organs, with a reduced size of between 25% and 52%. Furthermore, the pollen that developed was completely sterile, which was correlated with the altered expression of two Rf (fertility restorer)-like genes that encode pentatricopeptide repeat-containing proteins OsPPR676 and OsPPR920, translational initiation factors OseIF3e and OseIF3h, and the heat shock protein OsHSP82. These findings were verified through a yeast two-hybrid assay using a Nipponbare callus cDNA library as bait followed by the reverse transcription-PCR analysis of total leaf or anther RNAs. Our demonstration of the important role of Osj10gBTF3 in rice growth and development provides new insights showing that more complex regulatory functions are associated with BTF3 in plants.

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.

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.

The genome of melon (Cucumis melo L.)

We report the genome sequence of melon, an important horticultural crop worldwide. We assembled 375 Mb of the double-haploid line DHL92, representing 83.3% of the estimated melon genome. We predicted 27,427 protein-coding genes, which we analyzed by reconstructing 22,218 phylogenetic trees, allowing mapping of the orthology and paralogy relationships of sequenced plant genomes. We observed the absence of recent whole-genome duplications in the melon lineage since the ancient eudicot triplication, and our data suggest that transposon amplification may in part explain the increased size of the melon genome compared with the close relative cucumber. A low number of nucleotide-binding site-leucine-rich repeat disease resistance genes were annotated, suggesting the existence of specific defense mechanisms in this species. The DHL92 genome was compared with that of its parental lines allowing the quantification of sequence variability in the species. The use of the genome sequence in future investigations will facilitate the understanding of evolution of cucurbits and the improvement of breeding strategies.

Co-mapping studies of QTLs for fruit acidity and candidate genes of organic acid metabolism and proton transport in sweet melon (Cucumis melo L.)

Sweet melon cultivars contain a low level of organic acids and, therefore, the quality and flavor of sweet melon fruit is determined almost exclusively by fruit sugar content. However, genetic variability for fruit acid levels in the Cucumis melo species exists and sour fruit accessions are characterized by acidic fruit pH of <5, compared to the sweet cultivars that are generally characterized by mature fruit pH values of >6. In this paper, we report results from a mapping population based on recombinant inbred lines (RILs) derived from the cross between the non-sour 'Dulce' variety and the sour PI 414323 accession. Results show that a single major QTL for pH co-localizes with major QTLs for the two predominant organic acids in melon fruit, citric and malic, together with an additional metabolite which we identified as uridine. While the acidic recombinants were characterized by higher citric and malic acid levels, the non-acidic recombinants had a higher uridine content than did the acidic recombinants. Additional minor QTLs for pH, citric acid and malic acid were also identified and for these the increased acidity was unexpectedly contributed by the non-sour parent. To test for co-localization of these QTLs with genes encoding organic acid metabolism and transport, we mapped the genes encoding structural enzymes and proteins involved in organic acid metabolism, transport and vacuolar H+ pumps. None of these genes co-localized with the major pH QTL, indicating that the gene determining melon fruit pH is not one of the candidate genes encoding this primary metabolic pathway. Linked markers were tested in two additional inter-varietal populations and shown to be linked to the pH trait. The presence of the same QTL in such diverse segregating populations suggests that the trait is determined throughout the species by variability in the same gene and is indicative of a major role of the evolution of this gene in determining the important domestication trait of fruit acidity within the species.

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.

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.

Characterization of transcriptome dynamics during watermelon fruit development: sequencing, assembly, annotation and gene expression profiles

Background

Cultivated watermelon [Citrullus lanatus (Thunb.) Matsum. & Nakai var. lanatus] is an important agriculture crop world-wide. The fruit of watermelon undergoes distinct stages of development with dramatic changes in its size, color, sweetness, texture and aroma. In order to better understand the genetic and molecular basis of these changes and significantly expand the watermelon transcript catalog, we have selected four critical stages of watermelon fruit development and used Roche/454 next-generation sequencing technology to generate a large expressed sequence tag (EST) dataset and a comprehensive transcriptome profile for watermelon fruit flesh tissues.

Results

We performed half Roche/454 GS-FLX run for each of the four watermelon fruit developmental stages (immature white, white-pink flesh, red flesh and over-ripe) and obtained 577,023 high quality ESTs with an average length of 302.8 bp. De novo assembly of these ESTs together with 11,786 watermelon ESTs collected from GenBank produced 75,068 unigenes with a total length of approximately 31.8 Mb. Overall 54.9% of the unigenes showed significant similarities to known sequences in GenBank non-redundant (nr) protein database and around two-thirds of them matched proteins of cucumber, the most closely-related species with a sequenced genome. The unigenes were further assigned with gene ontology (GO) terms and mapped to biochemical pathways. More than 5,000 SSRs were identified from the EST collection. Furthermore we carried out digital gene expression analysis of these ESTs and identified 3,023 genes that were differentially expressed during watermelon fruit development and ripening, which provided novel insights into watermelon fruit biology and a comprehensive resource of candidate genes for future functional analysis. We then generated profiles of several interesting metabolites that are important to fruit quality including pigmentation and sweetness. Integrative analysis of metabolite and digital gene expression profiles helped elucidating molecular mechanisms governing these important quality-related traits during watermelon fruit development.

Conclusion

We have generated a large collection of watermelon ESTs, which represents a significant expansion of the current transcript catalog of watermelon and a valuable resource for future studies on the genomics of watermelon and other closely-related species. Digital expression analysis of this EST collection allowed us to identify a large set of genes that were differentially expressed during watermelon fruit development and ripening, which provide a rich source of candidates for future functional analysis and represent a valuable increase in our knowledge base of watermelon fruit biology.

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.