Alternatively spliced BobCAL transcripts alter curd morphotypes in a collection of Chinese cauliflower accessions

Abstract

The curd of cauliflower (Brassica oleracea L. var. botrytis) is a modified inflorescence that is consumed as a vegetable. Curd formation is proposed to be due to a mutation in the BobCAULIFLOWER (BobCAL) gene, but the genetic relationship between BobCAL variation and curd morphotypes remains obscure. To address this question, we collected and classified a collection of 78 cauliflower accessions into four subpopulations according to curd surface features: smooth, coarse, granular, and hairy curd morphotypes. Through the cDNA sequencing of BobCAL alleles, we showed that smooth and coarse accessions characterized by inflorescence meristem arrest presented a strong association with the 451T SNP (BobCAL_T), whereas granular and hairy accessions marked with floral organ arrest presented an association with 451G (BobCAL_G). Interestingly, all BobCAL alleles were alternatively spliced, resulting in a total of four alternative splice (AS) variants due to the retention of the fourth and/or seventh introns. Among accessions with BobCAL_G alleles, the total expression of all these AS variants in granular plants was almost equal to that in hairy plants; however, the expression of the individual AS variants encoding intact proteins relative to those encoding truncated proteins differed. Hairy accessions showed relatively high expression of the individual variants encoding intact proteins, whereas granular accessions displayed relatively low expression. In smooth cauliflower, the overexpression of the BobCAL_Ga variant caused an alteration in the curd morphotype from smooth to hairy, concurrent with an increase in the expression levels of downstream floral identity genes. These results reveal that alternative splicing of BobCAL transcripts is involved in the determination of cauliflower curd morphotypes.

References

[1] Smyth, D. R. Flower development: origin of the cauliflower. Curr. Biol. 5, 361-363 (1995).

[2] Kieffer, M., Fuller, M. & Jellings, A. Mathematical model of cauliflower curd architecture based on biometrical analysis. Acta Hortic. 407, 361-368 (1996).

[3] Sadik, S. Morphology of the curd of cauliflower. Am. J. Bot. 49, 290-297 (1962).

[4] De Candolle, A. Prodromus (Treuttel and Wurtz, 1824).

[5] Masters, M. T. On the structure of the flower in the genus Napoleona, &c. Botanical J. Linn. Soc. 10, 492-504 (1869).

[6] Sadik, S. & Ozbun, J. The association of carbohydrate changes in the shoot tip of cauliflower with flowering. Plant Physiol. 43, 1696-1698 (1968).

[7] Anthony, R. G., James, P. E. & Jordan, B. R. Cauliflower (Brassica oleracea var. botrytis L.) curd development: the expression of meristem identity genes. J. Exp. Bot. 47, 181-188 (1996).

[8] Carr, S. M. & Irish, V. F. Floral homeotic gene expression defines developmental arrest stages in Brassica oleracea L. vars. botrytis and italica. Planta 201, 179-188 (1997).

[9] Bowman, J. L. et al. Control of flower development in Arabidopsis thaliana by APETALA1 and interacting genes. Development 119, 721-743 (1993).

[10] Kempin, S. A., Savidge, B. & Yanofsky, M. F. Molecular basis of the cauliflower phenotype in Arabidopsis. Science 267, 522-525 (1995).

[11] Ferrandiz, C., Gu, Q., Martienssen, R. & Yanofsky, M. F. Redundant regulation of meristem identity and plant architecture by FRUITFULL, APETALA1 and CAULIFLOWER. Development 127, 725-734 (2000).

[12] Alvarez-Buylla, E. R., Garcia-Ponce, B. & Garay-Arroyo, A. Unique and redundant functional domains of APETALA1 and CAULIFLOWER, two recently duplicated Arabidopsis thaliana floral MADS-box genes. J. Exp. Bot. 57, 3099-3107 (2006).

[13] Pelaz, S. et al. B and C floral organ identity functions require SEPALLATA MADS-box genes. Nature 405, 200-203 (2000).

[14] Pelaz, S. et al. APETALA1 and SEPALLATA3 interact to promote flower development. Plant J. Cell Mol. Biol. 26, 385-394 (2001).

[15] Samach, A. et al. Distinct roles of CONSTANS target genes in reproductive development of Arabidopsis. Science 288, 1613-1616 (2000).

[16] Hartmann, U. et al. Molecular cloning of SVP: a negative regulator of the floral transition in Arabidopsis. Plant J. Cell Mol. Biol. 21, 351-360 (2000).

[17] Krizek, B. A. & Meyerowitz, E. M. Mapping the protein regions responsible for the functional specificities of the Arabidopsis MADS domain organ-identity proteins. Proc. Natl Acad. Sci. USA 93, 4063-4070 (1996).

[18] Riechmann, J. L., Krizek, B. A. & Meyerowitz, E. M. Dimerization specificity of Arabidopsis MADS domain homeotic proteins APETALA1, APETALA3, PISTILLATA, and AGAMOUS. Proc. Natl Acad. Sci. USA 93, 4793-4798 (1996).

[19] Fan, H. Y., Hu, Y., Tudor, M. & Ma, H. Specific interactions between the K domains of AG and AGLs, members of the MADS domain family of DNA binding proteins. Plant J. 12, 999-1010 (1997).

[20] Mizukami, Y. et al. Functional domains of the floral regulator AGAMOUS: characterization of the DNA binding domain and analysis of dominant negative mutations. Plant Cell 8, 831-845 (1996).

[21] Moon, Y.-H. et al. Determination of the motif responsible for interaction between the rice APETALA1/AGAMOUS-LIKE9 family proteins using a yeast two-hybrid system. Plant Physiol. 120, 1193-1204 (1999).

[22] Yang, Y. & Jack, T. Defining subdomains of the K domain important for protein-protein interactions of plant MADS proteins. Plant Mol. Biol. 55, 45 (2004).

[23] Huang, H. et al. The Arabidopsis MADS-box gene AGL3 is widely expressed and encodes a sequence-specific DNA-binding protein. Plant Mol. Biol. 28, 549-567 (1995).

[24] Smith, L. B. & King, G. J. The distribution of BoCAL-a alleles in Brassica oleracea is consistent with a genetic model for curd development and domestication of the cauliflower. Mol. Breed. 6, 603-613 (2000).

[25] Purugganan, M. D., Boyles, A. L. & Suddith, J. I. Variation and selection at the CAULIFLOWER floral homeotic gene accompanying the evolution of domesticated Brassica oleracea. Genetics 155, 855-862 (2000).

[26] Labate, J. A., Robertson, L. D., Baldo, A. M. & Björkman, T. Inflorescence identity gene alleles are poor predictors of inflorescence type in broccoli and cauliflower. J. Am. Soc. Hortic. Sci. 131, 667-673 (2006).

[27] Ratcliffe, O. J., Bradley, D. J. & Coen, E. S. Separation of shoot and floral identity in Arabidopsis. Development 126, 1109-1120 (1999).

[28] Kieffer, M. & Davies, B. Developmental programmes in floral organ formation. Semin. Cell Dev. Biol. 12, 373-380 (2001).

[29] Pidkowich, M. S., Klenz, J. E. & Haughn, G. W. The making of a flower: control of floral meristem identity in Arabidopsis. Trends Plant Sci. 4, 64-70 (1999).

[30] Zhao, Z. et al. Identification of candidate genes involved in curd riceyness in cauliflower. Int. J. Mol. Sci. 21, 1999 (2020).

[31] Yousef, E. A. A., Lampei, C. & Schmid, K. J. Evaluation of cauliflower genebank accessions under organic and conventional cultivation in Southern Germany. Euphytica 201, 389-400 (2015).

[32] Labate, J. A., Robertson, L. D. & Baldo, A. M. Inflorescence identity gene alleles are poor predictors of inflorescence type in broccoli and cauliflower. J. Am. Soc. Hortic. Sci. 131, 667-673 (2006).

[33] Chen, H., Nelson, R. & Sherwood, J. Enhanced recovery of transformants of Agrobacterium tumefaciens after freeze-thaw transformation and drug selection. Biotechniques 16, 664-668, 670 (1994).

[34] Sambrook, J., Fritsch, E. & Maniatis, T. Molecular Cloning: A Laboratory Manual. (Cold Spring Harbor Press, New York, 1989).

[35] Porebski, S., Bailey, L. G. & Baum, B. R. Modification of a CTAB DNA extraction protocol for plants containing high polysaccharide and polyphenol components. Plant Mol. Biol. Rep. 15, 8-15 (1997).

[36] Thompson, J. D. et al. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 25, 4876-4882 (1997).

[37] Tukey, J. W. Exploratory Data Analysis (Addison-Wesley, Reading, 1977).

[38] Kruskal, W. H. & Wallis, W. A. Use of ranks in one-criterion variance analysis. J. Am. Stat. Assoc. 47, 583-621 (1952).