What makes turnips: anatomy, physiology and transcriptome during early stages of its hypocotyl-tuber development

Abstract

Brassica species are characterized by their tremendous intraspecific diversity, exemplified by leafy vegetables, oilseeds, and crops with enlarged inflorescences or above ground storage organs. In contrast to potato tubers that are edible storage organs storing energy as starch and are the vegetative propagation modules, the storage organs of turnips, grown from true seed, are swollen hypocotyls with varying degrees of root and stem that mainly store glucose and fructose. To highlight their anatomical origin, we use the term “hypocotyl-tuber” for these turnip vegetative storage organs. We combined cytological, physiological, genetic and transcriptomic approaches, aiming to identify the initial stages, molecular pathways and regulatory genes for hypocotyl-tuber induction in turnips (B. rapa subsp. rapa). We first studied the development of the hypocotyl zone of turnip and Pak choi and found that 16 days after sowing (DAS) morphological changes occurred in the xylem which indicated the early tuberization stage. Tissue culture experiments showed a clear effect of auxin on hypocotyl-tuber growth. Differentially expressed genes between 1 and 6 weeks after sowing in turnip hypocotyls, located in genomic regions involved in tuber initiation and/or tuber growth defined by QTL and selective sweeps for tuber formation, were identified as candidate genes that were studied in more detail for their role in hypocotyl-tuber formation. This included a Bra-FLOR1 paralogue with increased expression 16 DAS, when the hypocotyl starts swelling, suggesting dual roles for duplicated flowering time genes in flowering and hypocotyl-tuber induction. Bra-CYP735A2 was identified for its possible role in tuber growth via trans-zeatin. Weighted Co-expression Network Analysis (WGCNA) identified 59 modules of co-expressed genes. Bra-FLOR1 and Bra-CYP735A2 were grouped in a module that included several genes involved in carbohydrate transport and metabolism, cell-wall growth, auxin regulation and secondary metabolism that serve as starting points to illuminate the transcriptional regulation of hypocotyl-tuber formation and development.

References

[1] Gowers, S. “ Swedes and Turnips,” in Root and Tuber Crops (ed Bradshaw, J. E.) 245-289 (Springer New York Dordrecht Heidelberg London, 2010).

[2] Ewing, E. E. & Struik, P. C. Tuber formation in potato: induction, initiation, and growth. Hortic. Rev. 14, 89-198 (1992).

[3] Jackson, S. D. Multiple signaling pathways control tuber induction in potato. Plant Physiol. 119, 1-8 (1999).

[4] Abelenda, J. A., Navarro, C. & Prat, S. From the model to the crop: genes controlling tuber formation in potato. Curr. Opin. Biotechnol. 22, 287-292 (2011).

[5] Aksenova, N. P., Konstantinova, T. N., Golyanovskaya, S. A., Sergeeva, L. I. & Romanov, G. A. Hormonal regulation of tuber formation in potato plants. Russ. J. Plant Physiol. 59, 451-466 (2012).

[6] Sonnewald, S. & Sonnewald, U. Regulation of potato tuber sprouting. Planta 239, 27-38 (2014).

[7] Bonnema, G., Del Carpio, D. P. & Zhao, J. J. Diversity analysis and molecular taxonomy of Brassica vegetable crops. in Genetics, Genomics and Breeding of Vegetable Brassicas (ed Sadowski, J. and Kole, C.) 81-124 (CRC Press, New York, 2011).

[8] Zhang, N. W. et al. Morphology, carbohydrate composition and vernalization response in a genetically diverse collection of Asian and European turnips (Brassica rapa subsp. rapa). PLoS One 9, e114241 (2014).

[9] Lu, G., Cao, J. S., Yu, X. L., Xiang, X. & Chen, H. Mapping QTLs for root morphological traits in Brassica rapa L. based on AFLP and RAPD markers. J. Appl. Genet. 49, 23-31 (2008).

[10] Jang, G. et al. Cytokinin-dependent secondary growth determines root biomass in radish (Raphanus sativus L.). J. Exp. Bot. 66, 4607-4619 (2015).

[11] Lou, P. et al. Quantitative trait loci for flowering time and morphological traits in multiple populations of Brassica rapa. J. Exp. Bot. 58, 4005-4016 (2007).

[12] Zhao, J. J. et al. Design of a Brassica rapa core collection for association mapping studies. Genome 53, 884-898 (2010).

[13] Lin, K. et al. Beyond genomic variation-comparison and functional annotation of three Brassica rapa genomes: a turnip, a rapid cycling and a Chinese cabbage. BMC Genom. 15, 250 (2014).

[14] Wang, X. W. et al. The genome of the mesopolyploid crop species Brassica rapa. Nat. Genet. 43, 1035-1039 (2011).

[15] Cheng, F. et al. Subgenome parallel selection is associated with morphotype diversification and convergent crop domestication in Brassica rapa and Brassica oleracea. Nat. Genet. 48, 1218-1224 (2016).

[16] Trick, M. et al. Complexity of genome evolution by segmental rearrangement in Brassica rapa revealed by sequence-level analysis. BMC Genom. 10, 539 (2009).

[17] Cheng, F. et al. Deciphering the diploid ancestral genome of the mesohexaploid Brassica rapa. Plant Cell 25, 1541-1554 (2013).

[18] Cheng, F., Wu, J. & Wang, X. W. Genome triplication drove the diversification of Brassica plants. Hortic. Res. 1, 14024 (2014).

[19] Kitashiba, H. et al. Draft sequences of the radish (Raphanus sativus L.) genome. DNA Res. 21, 481-490 (2014).

[20] Peterson, R. L. Control of cambial activity in roots of turnip (Brassica rapa). Can. J. Bot. 51, 475-480 (1973).

[21] Takahashi, H., Kimura, M., Suge, H. & Saito, T. Interactions between vernalization and photoperiod on the flowering and bolting of different turnip varieties. J. Jpn. Soc. Hortic. Sci. 63, 99-108 (1994).

[22] Kubo, N. et al. Detection of quantitative trait loci controlling morphological traits in Brassica rapa L.. Breed. Sci. 60, 164-171 (2010).

[23] Szklarczyk, D. et al. String v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 43, D447-D452 (2015).

[24] Petrasch, S. Temporal and Spatial Gene Expression During Turnip Formation in Brassica rapa. MSc Thesis, Wageningen University, Wageningen (2016).

[25] Livak, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitative PCR and the 2^−ΔΔCT method . Methods 25, 402-408 (2001).

[26] Basnet, R. K. et al. Genome-wide analysis of coordinated transcript abundance during seed development in different Brassica rapa morphotypes. BMC Genom. 14, 840 (2013).

[27] Smyth, G. K. Limma: linear models for microarray data. in Bioinformatics and Computational Biology Solutions Using R and Bioconductor (ed Gentleman, R., Carey, V., Dudoit, S., Irizarry, R. and Huber, W.) 397-420 (Springer, New York, 2005).

[28] Langfelder, P. & Horvath, S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinf. 9, 559 (2008).

[29] Zhang, B. & Horvath, S. A general framework for weighted gene co-expression network analysis. Stat. Appl. Genet. Mol. Biol. 4, 17 (2005).

[30] Langfelder, P., Zhang, B. & Horvath, S. Defining clusters from a hierarchical cluster tree: the Dynamic Tree Cut package for R. Bioinformatics 24, 719-720 (2008).

[31] Bassetti, N. Genome-wide Transcriptomic Profiling Of Turnip (Brassica rapa ssp. rapa) During Tuber Development. MSc thesis, Wageningen University (2015).

[32] Schmitt, U. et al. The vascular cambium of trees and its involvement in defining xylem anatomy. in Secondary Xylem Biology: Origins, Functions, and Applications (ed Kim, Y. S., Funada, R. & Singh, A. P.) 3-22 (Academic Press, London, 2016).

[33] Lachaud, S., Catesson, A. M. & Bonnemain, J. L. Structure and function of the vascular cambium. C. R. Acad. Sci. Ser. III 322, 633-650 (1999).

[34] Roumeliotis, E., Visser, R. G. F. & Bachem, C. W. B. A crosstalk of auxin and GA during tuber development. Plant Signal. Behav. 7, 1360-1363 (2012).

[35] Gukasyan, I. A. et al. Effect of Rol transgenes, IAA, and kinetin on starch content and the size of starch granules in tubers of in vitro potato plants. Russ. J. Plant Physiol. 52, 809-813 (2005).

[36] Kloosterman, B., Vorst, O., Hall, R. D., Visser, R. G. & Bachem, C. W. Tuber on a chip: differential gene expression during potato tuber development. Plant Biotechnol. J. 3, 505-519 (2005).

[37] Kloosterman, B. et al. Genes driving potato tuber initiation and growth: identification based on transcriptional changes using the POCI array. Funct. Integr. Genom. 8, 329-340 (2008).

[38] Firon, N. et al. Transcriptional profiling of sweetpotato (Ipomoea batatas) roots indicates down-regulation of lignin biosynthesis and up-regulation of starch biosynthesis at an early stage of storage root formation. BMC Genom. 14, 460 (2013).

[39] Mitsui, Y. et al. The radish genome and comprehensive gene expression profile of tuberous root formation and development. Sci. Rep. 5, 10835 (2015).

[40] Sun, X. C. et al. Identification of novel and salt-responsive miRNAs to explore miRNA-mediated regulatory network of salt stress response in radish (Raphanus sativus L.). BMC Genom. 16, 197 (2015).

[41] Parry, G., Ward, S., Cernac, A., Dharmasiri, S. & Estelle, M. The Arabidopsis suppressor of auxin resistance proteins are nucleoporins with an important role in hormone signaling and development. Plant Cell 18, 1590 (2006).

[42] Sampedro, J. & Cosgrove, D. J. The expansin superfamily. Genome Biol. 6, 242 (2005).

[43] Li, X. X., Zhao, J., Walk, T. C. & Liao, H. Characterization of soybean β-expansin genes and their expression responses to symbiosis, nutrient deficiency, and hormone treatment. Appl. Microbiol. Biotechnol. 98, 2805-2817 (2014).

[44] Gupta, A. K., Singh, J. & Kaur, N. Sink development, sucrose metabolising enzymes and carbohydrate status in turnip (Brassica rapa L.). Acta Physiol. Plant. 23, 31-36 (2001).

[45] Farrar, J. Regulation of root weight ratio is mediated by sucrose: Opinion. Plant Soil 185, 13-19 (1996).

[46] Sturm, A. & Tang, G. Q. The sucrose-cleaving enzymes of plants are crucial for development, growth and carbon partitioning. Trends Plant Sci. 4, 401-407 (1999).

[47] Yamada, K. et al. Monosaccharide absorption activity of Arabidopsis roots depends on expression profiles of transporter genes under high salinity conditions. J. Biol. Chem. 286, 43577-43586 (2011).

[48] Acevedo, F. G. et al. FLOR1, a putative interaction partner of the floral homeotic protein AGAMOUS, is a plant-specific intracellular LRR. Plant Sci. 167, 225-231 (2004).

[49] Torti, S. et al. Analysis of the Arabidopsis shoot meristem transcriptome during floral transition identifies distinct regulatory patterns and a leucine-rich repeat protein that promotes flowering. Plant Cell 24, 444-462 (2012).

[50] Kloosterman, B. A. et al. Naturally occurring allele diversity allows potato cultivation in northern latitudes. Nature 495, 246-250 (2013).

[51] Schmitz, R. Y., Skoog, F., Playtis, A. J. & Leonard, N. J. Cytokinins: synthesis and biological activity of geometric and position isomers of zeatin. Plant Physiol. 50, 702-705 (1972).

[52] Takei, K., Yamaya, T. & Sakakibara, H. Arabidopsis CYP735A1 and CYP735A2 encode cytokinin hydroxylases that catalyze the biosynthesis of trans-zeatin. J. Biol. Chem. 279, 41866-41872 (2004).

[53] Kiba, T., Takei, K., Kojima, M. & Sakakibara, H. Side-chain modification of cytokinins controls shoot growth in Arabidopsis. Dev. Cell 27, 452-461 (2013).

[54] Matsumoto-Kitano, M. et al. Cytokinins are central regulators of cambial activity. Proc. Natl. Acad. Sci. USA 105, 20027-20031 (2008).

[55] Ramireddy, E., Chang, L. & Schmülling, T. Cytokinin as a mediator for regulating root system architecture in response to environmental cues. Plant Signal. Behav. 9, 5021-5032 (2014).

[56] Paponov, I. A. et al. Comprehensive transcriptome analysis of auxin responses in Arabidopsis. Mol. Plant 1, 321-337 (2008).

[57] Kathare, P. K. Functional Characterization of SAUR Genes In Plant Growth and Development. PhD Thesis, Texas State University, Texas (2015).