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

Abstract

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.

References

[1] Ferriol M, Picó B . Pumpkin and winter squash. In: Prohens J, Nuez F (ed.) Vegetables I. Vol. 1. New York: Springer, 2008: 317-349.

[2] Azevedo-Meleiro CH, Rodriguez-Amaya DB . Qualitative and quantitative differences in carotenoid composition among Cucurbita moschata, Cucurbita maxima, and Cucurbita pepo . J Agric Food Chem 2007; 55: 4027-4033.

[3] Cazzonelli CI, Pogson BJ . Source to sink: regulation of carotenoid biosynthesis in plants. Trends Plant Sci 2010; 15: 266-274.

[4] Idouraine A, Kohlhepp EA, Weber CW, Warid WA, MartinezTellez JJ . Nutrient constituents from eight lines of naked seed squash (Cucurbita pepo L) . J Agric Food Chem 1996; 44: 721-724.

[5] Murkovic M, Mulleder U, Neunteufl H . Carotenoid content in different varieties of pumpkins. J Food Compos Anal 2002; 15: 633-638.

[6] Yoshida H, Shougaki Y, Hirakawa Y, Tomiyama Y, Mizushina Y . Lipid classes, fatty acid composition and triacylglycerol molecular species in the kernels of pumpkin (Cucurbita spp.) seeds . J Sci Food Agric 2004; 84: 158-163.

[7] Robinson RW, Decker-Walters DS . Cucurbits. New York: CAB International, 1997.

[8] Gajewski M, Radzanowska J, Danilcenko H, Jariene E, Cerniauskiene J . Quality of pumpkin cultivars in relation to sensory characteristics. Not Bot Hort Agrobot Cluj 2008; 36: 73-79.

[9] Itle RA, Kabelka EA . Correlation between L*a*b* color space values and carotenoid content in pumpkins and squash (Cucurbita spp.) . HortScience 2009; 44: 633-637.

[10] Cumarasamy R, Corrigan V, Hurst P, Bendall M . Cultivar differences in New Zealand ‘‘Kabocha’’ (buttercup squash, Cucurbita maxima) . NZ J Crop Hort Sci 2002; 30: 197-208.

[11] Corrigan VK, Hurst PL, Potter JF . Winter squash (Cucurbita maxima) texture: sensory, chemical, and physical measures . NZ J Crop Hort Sci 2001; 29: 111-124.

[12] Hurst PL, Corrigan VK, Koolaard J . Genetic analysis of sweetness and textural attributes in winter squash (Cucurbita maxima) . NZ J Crop Hort Sci 2006; 34: 359-367.

[13] Irving DE, Hurst PL, Ragg JS . Changes in carbohydrates and carbohydrate metabolizing enzymes during the development, maturation, and ripening of buttercup squash (Cucurbita maxima D ‘Delica’) . J Am Soc Hort Sci 1997; 122: 310-314.

[14] Lu S, Li L . Carotenoid metabolism: biosynthesis, regulation, and beyond. J Integr Plant Biol 2008; 50: 778-785.

[15] Huang SW, Li RQ, Zhang ZH et al. The genome of the cucumber, Cucumis sativus L. Nat Genet 2009; 41: 1275-1281.

[16] Clepet C, Joobeur T, Zheng Y et al. Analysis of expressed sequence tags generated from full-length enriched cDNA libraries of melon. BMC Genomics 2011; 12: 1-12.

[17] Diaz A, Fergany M, Formisano G et al. A consensus linkage map for molecular markers and Quantitative Trait Loci associated with economically important traits in melon (Cucumis melo L.) . BMC Plant Biol 2011; 11: 1-14.

[18] Garcia-Mas J, Benjak A, Sanseverino W et al. The genome of melon (Cucumis melo L.) . Proc Natl Acad Sci USA 2012; 109: 11872-11877.

[19] Ren Y, Zhao H, Kou Q et al. A high resolution genetic map anchoring scaffolds of the sequenced watermelon genome. PLoS One 2012; 7: e29453.

[20] Zhang WW, Pan JS, He HL et al. Construction of a high density integrated genetic map for cucumber (Cucumis sativus L.) . Theor Appl Genet 2012; 124: 249-259.

[21] Guo SG, Zhang JG, Sun HH et al. The draft genome of watermelon (Citrullus lanatus) and resequencing of 20 diverse accessions . Nat Genet 2013; 45: 51-58.

[22] Arumuganathan K, Earle ED . Nuclear DNA content of some important plant species. Plant Mol Biol Rep 1991; 9: 208-218.

[23] Weeden F. Isozyme studies indicate that the genus Cucurbita is an ancient tetraploid . Cucurbit Genet Cooper Rep 1984; 7: 84-85.

[24] Zraidi A, Stift G, Pachner M, Shojaeiyan A, Gong L, Lelley T . A consensus map for Cucurbita pepo . Mol Breeding 2007; 20: 375-388.

[25] Gong L, Stift G, Kofler R, Pachner M, Lelley T . Microsatellites for the genus Cucurbita and an SSR-based genetic linkage map of Cucurbita pepo L. Theor Appl Genet 2008; 117: 37-48.

[26] Gong L, Pachner M, Kalai K, Lelley T . SSR-based genetic linkage map of Cucurbita moschata and its synteny with Cucurbita pepo . Genome 2008; 51: 878-887.

[27] Esteras C, Gomez P, Monforte AJ et al. High-throughput SNP genotyping in Cucurbita pepo for map construction and quantitative trait loci mapping . BMC Genomics 2012; 13: 1-21.

[28] Blanca J, Canizares J, Roig C, Ziarsolo P, Nuez F, Pico B . Transcriptome characterization and high throughput SSRs and SNPs discovery in Cucurbita pepo (Cucurbitaceae) . BMC Genomics 2011; 12: 1-15.

[29] Wu T, Luo S, Wang R et al. The first Illumina-based de novo transcriptome sequencing and analysis of pumpkin (Cucurbita moschata Duch.) and SSR marker development . Mol Breeding 2014; 34: 1437-1447.

[30] Ando K, Grumet R . Transcriptional profiling of rapidly growing cucumber fruit by 454-pyrosequencing analysis. J Am Soc Hort Sci 2010; 135: 291-302.

[31] Ando K, Carr KM, Grumet R . Transcriptome analyses of early cucumber fruit growth identifies distinct gene modules associated with phases of development. BMC Genomics 2012; 13: 1-16.

[32] Zenoni S, Ferrarini A, Giacomelli E et al. Characterization of transcriptional complexity during berry development in Vitis vinifera using RNA-Seq . Plant Physiol 2010; 152: 1787-1795.

[33] Rohrmann J, Tohge T, Alba R et al. Combined transcription factor profiling, microarray analysis and metabolite profiling reveals the transcriptional control of metabolic shifts occurring during tomato fruit development. Plant J 2011; 68: 999-1013.

[34] Loy JB . Morpho-physiological aspects of productivity and quality in squash and pumpkins (Cucurbita spp.) . Crit Rev Plant Sci 2004; 23: 337-363.

[35] Grabherr MG, Haas BJ, Yassour M et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol 2011; 29: 644-652.

[36] Zheng Y, Zhao LJ, Gao JP, Fei ZJ . iAssembler: a package for de novo assembly of Roche-454/Sanger transcriptome sequences. BMC Bioinformatics 2011; 12: 1-8.

[37] Schmieder R, Edwards R . Fast identification and removal of sequence contamination from genomic and metagenomic datasets. PLoS One 2011; 6: e17288.

[38] Altschul SF, Madden TL, Schaffer AA et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997; 25: 3389-3402.

[39] Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M . Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005; 21: 3674-3676.

[40] Li B, Dewey CN . RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics 2011; 12: 1-16.

[41] Iseli C, Jongeneel C, Bucher P . ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. Proc Int Conf Intell Syst Mol Biol. 1999; 1999: 138-148.

[42] Apweiler R, Martin MJ, O’Donovan C et al. Update on activities at the Universal Protein Resource (UniProt) in 2013. Nucleic Acids Res 2013; 41: D43-D47.

[43] Zdobnov EM, Apweiler R . InterProScan-an integration platform for the signature-recognition methods in InterPro. Bioinformatics 2001; 17: 847-848.

[44] Camon E, Magrane M, Barrell D et al. The gene ontology annotation (GOA) project: Implementation of GO in SWISS-PROT, TrEMBL, and InterPro. Genome Res 2003; 13: 662-672.

[45] Sanchez A, Salicru M, Ocana J . Statistical methods for the analysis of high-throughput data based on functional profiles derived from the Gene Ontology. J Stat Plann Infer 2007; 137: 3975-3989.

[46] Carlson M, Falcon S, Pages H, Li N . GO.db: a set of annotation maps describing the entire Gene Ontology. R package version 2.6.1 ed.

[47] Ashburner M, Ball CA, Blake JA et al. Gene Ontology: tool for the unification of biology. Nat Genet 2000; 25: 25-29.

[48] Murphy TM, Gao M . Multiple forms of formamidopyrimidine-DNA glycosylase produced by alternative splicing in Arabidopsis thaliana . J Photochem Photobiol B-Biol 2001; 61: 87-93.

[49] Li JG, Li XJ, Guo L et al. A subgroup of MYB transcription factor genes undergoes highly conserved alternative splicing in Arabidopsis and rice . J Exp Botany 2006; 57: 1263-1273.

[50] Lois LM, Rodriguez-Concepcion M, Gallego F, Campos N, Boronat A . Carotenoid biosynthesis during tomato fruit development: regulatory role of 1-deoxy-D-xylulose 5-phosphate synthase. Plant J 2000; 22: 503-513.

[51] Hirschberg J. Carotenoid biosynthesis in flowering plants. Curr Opin Plant Biol 2001; 4: 210-218.

[52] Harjes CE, Rocheford TR, Bai L et al. Natural genetic variation in lycopene epsilon cyclase tapped for maize biofortification. Science 2008; 319: 330-333.

[53] Lu S, van Eck J, Zhou X et al. The cauliflower or gene encodes a DnaJ cysteine-rich domain-containing protein that mediates high levels of beta-carotene accumulation. Plant Cell 2006; 18: 3594-3605.

[54] Huber SC, Huber JL . Role and regulation of sucrose-phosphate synthase in higher plants. Annu Rev Plant Physiol Plant Mol Biol 1996; 47: 431-444.

[55] Zhang HP, Wu JY, Qin GH et al. The role of sucrose-metabolizing enzymes in pear fruit that differ in sucrose accumulation. Acta Physiol Plant 2014; 36: 71-77.

[56] Hubbard NL, Huber SC, Pharr DM . Sucrose phosphate synthase and acid invertase as determinants of sucrose concentration in developing muskmelon (Cucumis melo L.) fruits . Plant Physiol 1989; 91: 1527-1534.

[57] Geigenberger P, Stitt M, Fernie AR . Metabolic control analysis and regulation of the conversion of sucrose to starch in growing potato tubers. Plant Cell Environ 2004; 27: 655-673.

[58] Irving DE, Shingleton GJ, Hurst PL . Starch degradation in buttercup squash (I). J Am Soc Hort Sci 1999; 124: 587-590.