Genetic analysis of the “head top shape” quality trait of Chinese cabbage and its association with rosette leaf variation

Abstract

The agricultural and consumer quality of Chinese cabbage is determined by its shape. The shape is defined by the folding of the heading leaves, which defines the head top shape (HTS). The overlapping HTS, in which the heading leaves curve inward and overlap at the top, is the shape preferred by consumers. To understand the genetic regulation of HTS, we generated a large segregating F₂ population from a cross between pak choi and Chinese cabbage, with phenotypes ranging from nonheading to heading with either outward curving or inward curving overlapping heading leaves. HTS was correlated with plant height, outer/rosette leaf length, and petiole length. A high-density genetic map was constructed. Quantitative trait locus (QTL) analysis resulted in the identification of 22 QTLs for leafy head-related traits, which included five HTS QTLs. Bulked segregant analysis (BSA) was used to confirm HTS QTLs and identify candidate genes based on informative single-nucleotide polymorphisms. Interestingly, the HTS QTLs colocalized with QTLs for plant height, outer/rosette leaf, and petiole length, consistent with the observed phenotypic correlations. Combined QTL analysis and BSA laid a foundation for molecular marker-assisted breeding of Chinese cabbage HTS and directions for further research on the genetic regulation of this trait.

References

[1] Xu, J. B. et al. Guidelines for the Conduct of Tests for Distinctness Uniformity and Stability-Chinese Cabbage (The National Standard of the People’s Republic of China, 2004).

[2] Horiguchi, G., Ferjani, A., Fujikura, U. & Tsukaya, H. Coordination of cell proliferation and cell expansion in the control of leaf size in Arabidopsis thaliana. J. Plant Res. 119, 37-42 (2006).

[3] Tsukaya, H. Organ shape and size: a lesson from studies of leaf morphogenesis. Curr. Opin. Plant Biol. 6, 57-62 (2003).

[4] Horiguchi, G., Fujikura, U., Ferjani, A., Ishikawa, N. & Tsukaya, H. Large-scale histological analysis of leaf mutants using two simple leaf observation methods: identification of novel genetic pathways governing the size and shape of leaves. Plant J. 48, 638-644 (2006).

[5] Sandalio, L. M., Rodríguez-Serrano, M. & Romero-Puertas, M. C. Leaf epinasty and auxin: a biochemical and molecular overview. Plant Sci. 253, 187-193 (2016).

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

[7] Yu, X. et al. Cloning and structural and expressional characterization of BcpLH gene preferentially expressed in folding leaf of Chinese cabbage. Sci. China Ser. C 43, 321-329 (2000).

[8] He, Y. K., Xue, W. X., Sun, Y. D., Yu, X. H. & Liu, P. L. Leafy head formation of the progenies of transgenic plants of Chinese cabbage with exogenous auxin genes. Cell Res. 10, 151-160 (2000).

[9] Sun, X. et al. Genetic analysis of Chinese cabbage reveals correlation between rosette leaf and leafy head variation. Front. Plant Sci. 9, 1455 (2018).

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

[11] Inoue, T., Kubo, N., Kondo, T. & Hirai, M. Detection of quantitative trait loci for heading traits in Brassica rapa using different heading types of Chinese cabbage. J. Hortic. Sci. Biotechnol. 90, 311-317 (2015).

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

[13] Ge, Y. et al. Mapping quantitative trait loci for leaf and heading-related traits in Chinese cabbage (Brassica rapa L. ssp. pekinesis). Hortic. Environ. Biotechnol. 52, 494 (2011).

[14] Gu, A. et al. Coupling Seq-BSA and RNA-Seq analyses reveal the molecular pathway and genes associated with heading type in Chinese cabbage. Front. Genet. 8, 176 (2017).

[15] Li, R. et al. Conjunctive analyses of BSA-Seq and BSR-Seq to reveal the molecular pathway of leafy head formation in Chinese cabbage. Plants 8, 603 (2019).

[16] Kim, M. S. et al. Identification and characterization of the leaf specific networks of inner and rosette leaves in Brassica rapa. Biochem. Biophys. Res. Commun. 490, 821-826 (2017).

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

[18] Wang, F. et al. Transcriptome analysis of rosette and folding leaves in Chinese cabbage using high-throughput RNA sequencing. Genomics 99, 299-307 (2012).

[19] Gao, L. W. et al. Genome-wide analysis of auxin transport genes identifies the hormone responsive patterns associated with leafy head formation in Chinese cabbage. Sci. Rep. 7, 1-13 (2017).

[20] Ren, W. et al. BcpLH organizes a specific subset of microRNAs to form a leafy head in Chinese cabbage (Brassica rapa ssp. pekinensis). Hortic. Res. 7, 1-13 (2020).

[21] Himelblau, E. & Amasino, R. M. Nutrients mobilized from leaves of Arabidopsis thaliana during leaf senescence. J. Plant Physiol. 158, 1317-1323 (2001).

[22] Wingler, A., Marès, M. & Pourtau, N. Spatial patterns and metabolic regulation of photosynthetic parameters during leaf senescence. N. Phytologist 161, 781-789 (2004).

[23] Shen, Y., Xiang, Y., Xu, E., Ge, X. & Li, Z. Major co-localized QTL for plant height, branch initiation height, stem diameter, and flowering time in an alien introgression derived Brassica napus DH population. Front. Plant Sci. 9, 390 (2018).

[24] Luo, H. et al. Co-localization of major quantitative trait loci for pod size and weight to a 3.7 cM interval on chromosome A05 in cultivated peanut (Arachis hypogaea L). BMC Genomics 18, 58 (2017).

[25] Weng, X. et al. Grain number, plant height, and heading date7 is a central regulator of growth, development, and stress response. Plant Physiol. 164, 735-747 (2014).

[26] Kalve, S., De Vos, D. & Beemster, G. T. S. Leaf development: a cellular perspective. Front. Plant Sci. 5, 362 (2014).

[27] Kaplan, D. R. & Hagemann, W. The relationship of cell and organism in vascular plants. Bioscience 41, 693-703 (1991).

[28] Garg, R., Singh, V. K., Rajkumar, M. S., Kumar, V. & Jain, M. Global transcriptome and coexpression network analyses reveal cultivar-specific molecular signatures associated with seed development and seed size/weight determination in chickpea. Plant J. 91, 1088-1107 (2017).

[29] Cao, X., Liu, B. & Zhang, Y. SEA: a software package of segregation analysis of quantitative traits in plants. Journal of Nanjing Agricultural University (2013).

[30] Elshire, R. J. et al. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS ONE 6(5), e19379 (2011).

[31] Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754-1760 (2009).

[32] Li, H. et al. The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078-2079 (2009).

[33] Stam, P. Construction of integrated genetic linkage maps by means of a new computer package: Join Map. Plant J. 3, 739-744 (1993).

[34] Tang, H. et al. ALLMAPS: robust scaffold ordering based on multiple maps. Genome Biol. 16, 3 (2015).

[35] Broman, K. W., Wu, H., Sen, Ś. & Churchill, G. A. R/qtl: QTL mapping in experimental crosses. Bioinformatics 19, 889-890 (2003).

[36] Churchill, G. A. & Doerge, R. W. Empirical threshold values for quantitative trait mapping. Genetics 138, 963-971 (1994).