Advances in molecular breeding in potato have been limited by its complex biological system, which includes vegetative propagation, autotetraploidy, and extreme heterozygosity. The availability of the potato genome and accompanying gene complement with corresponding gene structure, location, and functional annotation are powerful resources for understanding this complex plant and advancing molecular breeding efforts. Here, we report a reference for the potato transcriptome using 32 tissues and growth conditions from the doubled monoploid Solanum tuberosum Group Phureja clone DM1-3 516R44 for which a genome sequence is available. Analysis of greater than 550 million RNA-Seq reads permitted the detection and quantification of expression levels of over 22,000 genes. Hierarchical clustering and principal component analyses captured the biological variability that accounts for gene expression differences among tissues suggesting tissue-specific gene expression, and genes with tissue or condition restricted expression. Using gene co-expression network analysis, we identified 18 gene modules that represent tissue-specific transcriptional networks of major potato organs and developmental stages. This information provides a powerful resource for potato research as well as studies on other members of the Solanaceae family.
References
[1]
Papademetriou MK (2008) (ed.) In: RAP Publication (FAO), no. 2008/07; Workshop to Commemorate the International Year of Potato, Bangkok (Thailand), 6 May 2008/FAO, Bangkok (Thailand). Regional Office for Asia and the Pacific, 2008. 84 p. RAP publication (FAO).
[2]
Mendiburu AO, Peloquin SJ (1977) The significance of 2N gametes in potato breeding. Theoretical and Applied Genetics 49: 53–61.
[3]
The Potato Genome Sequencing Consortium (2011) Genome sequence and analysis of the tuber crop potato. Nature 475: 189–195.
[4]
Crookshanks M, Emmersen J, Welinder KG, Lehmann Nielsen K (2001) The potato tuber transcriptome: analysis of 6077 expressed sequence tags. FEBS Letters 506: 123–126.
[5]
Flinn B, Rothwell C, Griffiths R, L？gue M, DeKoeyer D, et al. (2005) Potato Expressed Sequence Tag generation and analysis using standard and unique cDNA Libraries. Plant Molecular Biology 59: 407–433.
[6]
Li XQ (2007) EST Sequencing and analysis from cold-stored and reconditioned potato tubers. Acta Horticulturae 745: 6.
[7]
Rensink W, Hart A, Liu J, Ouyang S, Zismann V, et al. (2005) Analyzing the potato abiotic stress transcriptome using expressed sequence tags. Genome 48: 598–605.
[8]
Ronning CM, Stegalkina SS, Ascenzi RA, Bougri O, Hart AL, et al. (2003) Comparative analyses of potato Expressed Sequence Tag libraries. Plant Physiol 131: 419–429.
[9]
Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, et al. (2007) NCBI GEO: mining tens of millions of expression profiles–database and tools update. Nucleic Acids Research 35: D760–D765.
[10]
Parkinson H, Sarkans U, Shojatalab M, Abeygunawardena N, Contrino S, et al. (2005) ArrayExpress–a public repository for microarray gene expression data at the EBI. Nucleic Acids Research 33: D553–D555.
[11]
Rensink W, Iobst S, Hart A, Stegalkina S, Liu J, et al. (2005) Gene expression profiling of potato responses to cold, heat, and salt stress. Functional Integrative Genomics 5: 201–207.
[12]
Bachem C, van der Hoeven R, Lucker J, Oomen R, Casarini E, et al. (2000) Functional genomic analysis of potato tuber life-cycle. Potato Research 43: 297–312.
[13]
Restrepo S, Myers KL, del Pozo O, Martin GB, Hart AL, et al. (2005) Gene Profiling of a compatible interaction Between Phytophthora infestans and Solanum tuberosum suggests a role for carbonic anhydrase. Molecular Plant-Microbe Interactions 18: 913–922.
[14]
Evers D, Lefevre I, Legay S, Lamoureux D, Hausman J-F, et al. (2010) Identification of drought-responsive compounds in potato through a combined transcriptomic and targeted metabolite approach. J Exp Bot 61: 2327–2343.
[15]
Ginzberg I, Barel G, Ophir R, Tzin E, Tanami Z, et al. (2009) Transcriptomic profiling of heat-stress response in potato periderm. J Exp Bot 60: 4411–4421.
[16]
Kloosterman B, De Koeyer D, Griffiths R, Flinn B, Steuernagel B, et al. (2008) Genes driving potato tuber initiation and growth: identification based on transcriptional changes using the POCI array. Functional Integrative Genomics 8: 329–340.
[17]
Wang Z, Gerstein M, Snyder M (2009) RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet 10: 57–63.
[18]
Morin RD, Bainbridge M, Fejes A, Hirst A, Krzywinski M, et al. (2008) Profiling the HeLa S3 transcriptome using randomly primed cDNA and massively parallel short-read sequencing. Bio Techniques 45: 14.
[19]
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, et al. (2010) Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature Biotechology 28: 511–515.
[20]
Nicot N, Hausman J-Fo, Hoffmann L, Evers Dl (2005) Housekeeping gene selection for real-time RT-PCR normalization in potato during biotic and abiotic stress. Journal of Experimental Botany 56: 2907–2914.
[21]
Zhang B, Horvath S (2005) A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol 4: Article 17.
[22]
Suzuki M, Wang HH-Y, McCarty DR (2007) Repression of the LEAFY COTYLEDON 1/B3 regulatory network in plant embryo development by VP1/ABSCISIC ACID INSENSITIVE 3-LIKE B3 Genes. PLANT PHYSIOLOGY 143: 902–911.
Altschul SF, Madden TL, Sch？ffer AA, Zhang J, Zhang Z, et al. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Research 25: 3389–3402.
[25]
Zdobnov EM, Apweiler R (2001) InterProScan–an integration platform for the signature-recognition methods in InterPro. Bioinformatics 17: 847–848.
[26]
Langfelder P, Horvath S (2008) WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 9: 559.
[27]
Langfelder P, Horvath S (2007) Eigengene networks for studying the relationships between co-expression modules. BMC Systems Biology 1: 54.
