Divergence of three BRX homoeologs in Brassica rapa and its effect on leaf morphology

Abstract

The leafy head characteristic is a special phenotype of Chinese cabbage resulting from artificial selection during domestication and breeding. BREVIS RADIX (BRX) has been suggested to control root elongation, shoot growth, and tiller angle in Arabidopsis and rice. In Brassica rapa, three BrBRX homoeologs have been identified, but only BrBRX.1 and BrBRX.2 were found to be under selection in leaf-heading accessions, indicating their functional diversification in leafy head formation. Here, we show that these three BrBRX genes belong to a plant-specific BRX gene family but that they have significantly diverged from other BRX-like members on the basis of different phylogenetic classifications, motif compositions and expression patterns. Moreover, although the expression of these three BrBRX genes differed, compared with BrBRX.3, BrBRX.1, and BrBRX.2 displayed similar expression patterns. Arabidopsis mutant complementation studies showed that only BrBRX.1 could rescue the brx root phenotype, whereas BrBRX.2 and BrBRX.3 could not. However, overexpression of each of the three BrBRX genes in Arabidopsis resulted in similar pleiotropic leaf phenotypes, including epinastic leaf morphology, with an increase in leaf number and leaf petiole length and a reduction in leaf angle. These leaf traits are associated with leafy head formation. Further testing of a SNP (T/C) in BrBRX.2 confirmed that this allele in the heading accessions was strongly associated with the leaf-heading trait of B. rapa. Our results revealed that all three BrBRX genes may be involved in the leaf-heading trait, but they may have functionally diverged on the basis of their differential expression.

References

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

[2] Liu, S. et al. The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes. Nat. Commun. 5, 3930 (2014).

[3] Yang, J. et al. The genome sequence of allopolyploid Brassica juncea and analysis of differential homoeolog gene expression influencing selection. Nat. Genet. 48, 1225-1232 (2016).

[4] Cui, Y. et al. Segmental translocation contributed to the origin of the Brassica S-locus. Hortic. Plant J. 6, 167-178 (2020).

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

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

[7] Briggs, G. C., Mouchel, C. F. & Hardtke, C. S. Characterization of the plant-specific BREVIS RADIX gene family reveals limited genetic redundancy despite high sequence conservation. Plant Physiol. 140, 1306-1316 (2006).

[8] Beuchat, J. et al. A hyperactive quantitative trait locus allele of Arabidopsis BRX contributes to natural variation in root growth vigor. Proc. Natl Acad. Sci. USA. 107, 8475-8480 (2010).

[9] Mouchel, C. F., Briggs, G. C. & Hardtke, C. S. Natural genetic variation in Arabidopsis identifies BREVIS RADIX, a novel regulator of cell proliferation and elongation in the root. Genes Dev. 18, 700-714 (2004).

[10] Liu, J., Liang, D., Song, Y. & Xiong, L. Systematic identification and expression analysis of BREVIS RADIX-like homologous genes in rice. Plant Sci. 178, 183-191 (2010).

[11] Beuchat, J. et al. BRX promotes Arabidopsis shoot growth. N. Phytol. 188, 23-29 (2010).

[12] Li, Z. et al. OsBRXL4 regulates shoot gravitropism and rice tiller angle through affecting LAZY1 nuclear localization. Mol. Plant 12, 1143-1156 (2019).

[13] Blanc, G. & Wolfe, K. H. Functional divergence of duplicated genes formed by polyploidy during Arabidopsis evolution. Plant Cell 16, 1679-1691 (2004).

[14] Casneuf, T., De Bodt, S., Raes, J., Maere, S. & Van de Peer, Y. Nonrandom divergence of gene expression following gene and genome duplications in the flowering plant Arabidopsis thaliana. Genome Biol. 7, R13 (2006).

[15] Duarte, J. M. et al. Expression pattern shifts following duplication indicative of subfunctionalization and neofunctionalization in regulatory genes of Arabidopsis. Mol. Biol. Evol. 23, 469-478 (2006).

[16] Cheng, F. et al. BRAD, the genetics and genomics database for Brassica plants. BMC Plant Biol. 11, 1-6 (2011).

[17] Bowers, J. E., Chapman, B. A., Rong, J. & Paterson, A. H. Unravelling angiosperm genome evolution by phylogenetic analysis of chromosomal duplication events. Nature 422, 433-438 (2003).

[18] Zhang, J. et al. Three genes encoding AOP2, a protein involved in aliphatic glucosinolate biosynthesis, are differentially expressed in Brassica rapa. J. Exp. Bot. 66, 6205-6218 (2015).

[19] Tong, C., Wang, X., Yu, J., Wu, J. & Liu, S. Comprehensive analysis of RNA-seq data reveals the complexity of the transcriptome in Brassica rapa. BMC Genomics 14, 689 (2013).

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

[21] Wang, Y., Wu, F., Bai, J. & He, Y. BrpSPL9 (Brassica rapa ssp. pekinensis SPL9) controls the earliness of heading time in Chinese cabbage. Plant Biotechnol. J. 12, 312-321 (2014).

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

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

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

[25] Romano, C. P., Robson, P. R. H., Smith, H., Estelle, M. & Klee, H. Transgene-mediated auxin overproduction in Arabidopsis: hypocotyl elongation phenotype and interactions with the hy6-1 hypocotyl elongation and axr1 auxin-resistant mutants. Plant Mol. Biol. 27, 1071-1083 (1995).

[26] Kidner, C. A. & Timmermans, M. C. Mixing and matching pathways in leaf polarity. Curr. Opin. Plant Biol. 10, 13-20 (2007).

[27] Wu, F. et al. The N-terminal double-stranded RNA binding domains of Arabidopsis HYPONASTIC LEAVES1 are sufficient for pre-microRNA processing. Plant Cell 19, 914-925 (2007).

[28] Han, M. H., Goud, S., Song, L. & Fedoroff, N. The Arabidopsis double-stranded RNA-binding protein HYL1 plays a role in microRNA-mediated gene regulation. Proc. Natl Acad. Sci. USA. 101, 1093-1098 (2004).

[29] Liu, Z., Jia, L., Wang, H. & He, Y. HYL1 regulates the balance between adaxial and abaxial identity for leaf flattening via miRNA-mediated pathways. J. Exp. Bot. 62, 4367-4381 (2011).

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

[31] Liang, J., Liu, B., Wu, J., Cheng, F. & Wang, X. Genetic variation and divergence of genes involved in leaf adaxial-abaxial polarity establishment in Brassica rapa. Front. Plant Sci. 7, 94 (2016).

[32] Mouchel, C. F., Osmont, K. S. & Hardtke, C. S. BRX mediates feedback between brassinosteroid levels and auxin signalling in root growth. Nature 443, 458-461 (2006).

[33] Scacchi, E. et al. Dynamic, auxin-responsive plasma membrane-to-nucleus movement of Arabidopsis BRX. Development 136, 2059-2067 (2009).

[34] Marhava, P. et al. A molecular rheostat adjusts auxin flux to promote root protophloem differentiation. Nature 558, 297-300 (2018).

[35] Scacchi, E. et al. Spatio-temporal sequence of cross-regulatory events in root meristem growth. Proc. Natl Acad. Sci. USA 107, 22734-22739 (2010).

[36] Rodrigues, A. et al. The short-rooted phenotype of the brevis radix mutant partly reflects root abscisic acid hypersensitivity. Plant Physiol. 149, 1917-1928 (2009).

[37] Kliver, S. et al. Assembly of the Boechera retrofracta genome and evolutionary analysis of apomixis-associated genes. Genes 9, 185 (2018).

[38] Gan, X. et al. The Cardamine hirsuta genome offers insight into the evolution of morphological diversity. Nat. Plants 2, 16167 (2016).

[39] Zhang, T. et al. Genome of Crucihimalaya himalaica, a close relative of Arabidopsis, shows ecological adaptation to high altitude. Proc. Natl Acad. Sci. USA 116, 7137-7146 (2019).

[40] Dorn, K. M., Fankhauser, J. D., Wyse, D. L. & Marks, M. D. A draft genome of field pennycress (Thlaspi arvense) provides tools for the domestication of a new winter biofuel crop. DNA Res. 22, 121-131 (2015).