Combination of genetic analysis and ancient literature survey reveals the divergence of traditional Brassica rapa varieties from Kyoto, Japan

Abstract

Since ancient times, humans have bred several plants that we rely on today. However, little is known about the divergence of most of these plants. In the present study, we investigated the divergence of Mibuna (Brassica rapa L. subsp. nipposinica L. H. Bailey), a traditional leafy vegetable in Kyoto (Japan), by combining genetic analysis and a survey of ancient literature. Mibuna is considered to have been bred 200 years ago from Mizuna, another traditional leafy vegetable in Kyoto. Mibuna has simple spatulate leaves, whereas Mizuna has characteristic serrated leaves. The quantitative trait loci (QTL) and gene expression analyses suggested that the downregulation of BrTCP15 expression contributed to the change in the leaf shape from serrated to simple spatulate. Interestingly, the SNP analysis indicated that the genomic region containing the BrTCP15 locus was transferred to Mibuna by introgression. Furthermore, we conducted a survey of ancient literature to reveal the divergence of Mibuna and found that hybridization between Mizuna and a simple-leaved turnip might have occurred in the past. Indeed, the genomic analysis of multiple turnip cultivars showed that one of the cultivars, Murasakihime, has almost the same sequence in the BrTCP15 region as Mibuna. These results suggest that the hybridization between Mizuna and turnip has resulted in the establishment of Mibuna.

References

[1] Nakamura, Y. et al. Antimutagenic and anticarcinogenic properties of Kyo-yasai, heirloom vegetables in Kyoto. Genes Environ. 30, 41-47 (2008).

[2] Nakamura, T. et al. Protection of Kyo-yasai (heirloom vegetables in Kyoto) from extinction: a case of Sabaka-daikon (Japan’s heirloom white radish, Raphanus sativus) in Maizuru, Japan. J. Ethn. Foods 4, 103-109 (2017).

[3] Zhao, J. et al. Genetic relationships within Brassica rapa as inferred from AFLP fingerprints. Theor. Appl. Genet. 110, 1301-1314 (2005).

[4] Takuno, S., Kawahara, T. & Ohnishi, O. Phylogenetic relationships among cultivated types of Brassica rapa L. em. Metzg. as revealed by AFLP analysis. Genet. Resour. Crop Evol. 54, 279-285 (2007).

[5] Artemyeva, A. M., Budahn, H., Klocke, E. & Chesnokov, Yu. V. Use of CACTA mobile genetic elements for revision of phylogenetic relationships in Brassica rapa L. species. Russ. J. Genet. Appl. Res. 1, 577-586 (2011).

[6] Qi, X. et al. Genomic inferences of domestication events are corroborated by written records in Brassica rapa. Mol. Ecol. 26, 3373-3388 (2017).

[7] Lou, P. et al. Genetic and genomic resources to study natural variation in Brassica rapa. Plant Direct 4, 1-12 (2020).

[8] Kimura, S. & Kawakatsu, Y. Evolutionary history of leaf shape diversification found in Mizuna and Mibuna. Kyoto Sangyo Univ. Essays Humanit. Ser. No. 49 (2016).

[9] Martín-Trillo, M. & Cubas, P. TCP genes: a family snapshot ten years later. Trends Plant Sci. 15, 31-39 (2010).

[10] Kieffer, M., Master, V., Waites, R. & Davies, B. TCP14 and TCP15 affect internode length and leaf shape in Arabidopsis. Plant J. 68, 147-158 (2011).

[11] Cubas, P., Lauter, N., Doebley, J. & Coen, E. The TCP domain: a motif found in proteins regulating plant growth and development. Plant J. 18, 215-222 (1999).

[12] Lucero, L. E. et al. TCP15 modulates cytokinin and auxin responses during gynoecium development in Arabidopsis. Plant J. 84, 267-282 (2015).

[13] Ori, N. et al. Regulation of LANCEOLATE by miR319 is required for compound-leaf development in tomato. Nat. Genet. 39, 787-791 (2007).

[14] Seki, K. et al. A CIN-like TCP transcription factor (LsTCP4) having retrotransposon insertion associates with a shift from Salinas type to Empire type in crisphead lettuce (Lactuca sativa L.). Hortic. Res. 7, 15 (2020).

[15] Crawford, B. C. W., Nath, U., Carpenter, R. & Coen, E. S. CINCINNATA controls both cell differentiation and growth in petal lobes and leaves of Antirrhinum. Plant Physiol. 135, 244-253 (2004).

[16] Nath, U., Crawford, B. C. W., Carpenter, R. & Coen, E. Genetic control of surface curvature. Science 299, 1404-1407 (2003).

[17] Hay, A., Barkoulas, M. & Tsiantis, M. ASYMMETRIC LEAVES1 and auxin activities converge to repress BREVIPEDICELLUS expression and promote leaf development in Arabidopsis. Development 133, 3955-3961 (2006).

[18] Iwakawa, H. et al. The ASYMMETRIC LEAVES2 Gene of Arabidopsis thaliana, required for formation of a symmetric flat leaf lamina, encodes a member of a novel family of proteins characterized by cysteine repeats and a leucine zipper. Plant Cell Physiol. 43, 467-478 (2002).

[19] Shuai, B., Reynaga-Peña, C. G. & Springer, P. S. The lateral organ boundaries gene defines a novel, plant-specific gene family. Plant Physiol. 129, 747-761 (2002).

[20] Tiwari, S. B., Hagen, G. & Guilfoyle, T. The roles of auxin response factor domains in auxin-responsive transcription. Plant Cell 15, 533-543 (2003).

[21] Pekker, I., Alvarez, J. P. & Eshed, Y. Auxin response factors mediate Arabidopsis organ asymmetry via modulation of KANADI activity. Plant Cell 17, 2899-2910 (2005).

[22] Yamaguchi, T., Nukazuka, A. & Tsukaya, H. Leaf adaxial-abaxial polarity specification and lamina outgrowth: evolution and development. Plant Cell Physiol. 53, 1180-1194 (2012).

[23] Hanson, J., Johannesson, H. & Engström, P. Sugar-dependent alterations in cotyledon and leaf development in transgenic plants expressing the HDZhdip gene ATHB13. Plant Mol. Biol. 45, 247-262 (2001).

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

[25] Mao, Y. et al. microRNA319a-targeted Brassica rapa ssp. pekinensis TCP genes modulate head shape in Chinese cabbage by differential cell division arrest in leaf regions. Plant Physiol. 164, 710-720 (2014).

[26] Sun, X. et al. Genome-wide transcriptome analysis reveals molecular pathways involved in leafy head formation of Chinese cabbage (Brassica rapa). Hortic. Res. 6, 1-13 (2019).

[27] Yu, X. et al. QTL Mapping of leafy heads by genome resequencing in the RIL population of Brassica rapa. PLoS ONE 8, e76059 (2013).

[28] 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).

[29] Li, X. et al. Quantitative trait loci mapping in Brassica rapa revealed the structural and functional conservation of genetic loci governing morphological and yield component traits in the A, B, and C subgenomes of Brassica species. DNA Res. 20, 1-16 (2013).

[30] Li, F., Kitashiba, H., Inaba, K. & Nishio, T. A Brassica rapa linkage map of EST-based SNP markers for identification of candidate genes controlling flowering time and leaf morphological traits. DNA Res. 16, 311-323 (2009).

[31] Kawakatsu, Y. et al. A GLABRA1 ortholog on LG A9 controls trichome number in the Japanese leafy vegetables Mizuna and Mibuna (Brassica rapa L. subsp. nipposinica L. H. Bailey): evidence from QTL analysis. J. Plant Res. 130, 539-550 (2017).

[32] Zhang, J. et al. Map-based cloning and characterization of a gene controlling hairiness and seed coat color traits in Brassica rapa. Plant Mol. Biol. 69, 553-563 (2009).

[33] Wang, Y. et al. Identification and validation of a major QTL controlling the presence/absence of leaf lobes in Brassica rapa L. Euphytica 205, 761-771 (2015).

[34] Ni, X., Huang, J., Ali, B., Zhou, W. & Zhao, J. Genetic analysis and fine mapping of the LOBED-LEAF 1 (BnLL1) gene in rapeseed (Brassica napus L.). Euphytica 204, 29-38 (2015).

[35] Baird, N. A. et al. Rapid SNP discovery and genetic mapping using sequenced RAD markers. PLoS ONE 3, e3376 (2008).

[36] Takagi, H. et al. QTL-seq: rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations. Plant J. 74, 174-183 (2013).

[37] Abe, A. et al. Genome sequencing reveals agronomically important loci in rice using MutMap. Nat. Biotechnol. 30, 174-178 (2012).

[38] Palatnik, J. F. et al. Control of leaf morphogenesis by microRNAs. Nature 425, 257-263 (2003).

[39] Uberti-Manassero, N. G., Lucero, L. E., Viola, I. L., Vegetti, A. C. & Gonzalez, D. H. The class I protein AtTCP15 modulates plant development through a pathway that overlaps with the one affected by CIN-like TCP proteins. J. Exp. Bot. 63, 809-823 (2012).

[40] Uchiyama, Y., Fujihira, Y., Matsuoka, H. & Kohsaka, R. Tradition and Japanese vegetables: history, locality, geography, and discursive ambiguity. J. Ethn. Foods 4, 198-203 (2017).

[41] Shengchu W., Christopher J. B. and Zhao-Bang Z. Windows QTL Cartographer 2.5. Department of Statistics (North Carolina State University, 2012).

[42] Kim, D. et al. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14, R36 (2013).

[43] Anders, S., Pyl, P. T. & Huber, W. HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics 31, 166-169 (2015).

[44] Sun, J., Nishiyama, T., Shimizu, K. & Kadota, K. TCC: an R package for comparing tag count data with robust normalization strategies. BMC Bioinform. 14, 219 (2013).

[45] Langmead, B., Trapnell, C., Pop, M. & Salzberg, S. L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10, R25 (2009).

[46] Walker, B. J. et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE 9, e112963 (2014).

[47] McKenna, A. et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297-1303 (2010).

[48] Kumar, S., Stecher, G., Li, M., Knyaz, C. & Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547-1549 (2018).