Map-based cloning and CRISPR/Cas9-based editing uncover BoNA1 as the causal gene for the no-anthocyanin-accumulation phenotype in curly kale (Brassica oleracea var. sabellica)

Abstract

Brassica oleracea comprises several important vegetable and ornamental crops, including curly kale, ornamental kale, cabbage, broccoli, and others. The accumulation of anthocyanins, important secondary metabolites valuable to human health, in these plants varies widely and is responsible for their pink to dark purple colors. Some curly kale varieties lack anthocyanins, making these plants completely green. The genetic basis of this trait is still unknown. We crossed the curly kale inbred line BK2019 (without anthocyanins) with the cabbage inbred line YL1 (with anthocyanins) and the Chinese kale inbred line TO1000 (with anthocyanins) to generate segregating populations. The no-anthocyanin trait was genetically controlled by a recessive gene, bona1. We generated a linkage map and mapped bona1 to a 256-kb interval on C09. We identified one candidate gene, Bo9g058630, in the target genomic region; this gene is homologous to AT5G42800, which encodes a dihydroflavonol-4-reductase-like (DFR-like) protein in Arabidopsis. In BK2019, a 1-bp insertion was observed in the second exon of Bo9g058630 and directly produced a stop codon. To verify the candidate gene function, CRISPR/Cas9 gene editing technology was applied to knock out Bo9g058630. We generated three bona1 mutants, two of which were completely green with no anthocyanins, confirming that Bo9g058630 corresponds to BoNA1. Different insertion/deletion mutations in BoNA1 exons were found in all six of the other no-anthocyanin kale varieties examined, supporting that independent disruption of BoNA1 resulted in no-anthocyanin varieties of B. oleracea. This study improves the understanding of the regulation mechanism of anthocyanin accumulation in B. oleracea subspecies.

Acknowledgements

This work was funded by the Science and Technology Innovation Program of the Chinese Academy of Agricultural Sciences (CAAS- ASTIP-IVFCAAS) and the China Agriculture Research System of MOF and MARA (CARS-23).

Author contributions

H.L. designed the work. K.Y. and X.Z. carried out the experiments. K.Y. and X.Z. wrote and revised the manuscript. F.H. helped with data analyses. W.S., L.Y., Y.Z., Y.W., J.J., and Z.F. assisted with the experiments.

Data availability

All the data supporting this work are presented in the manuscript and its supplementary files.

Conflict of interest statement

The authors declare no conflict of interest.

References

[1] Tanaka Y, Sasaki N, Ohmiya A . Biosynthesis of plant pigments: anthocyanins, betalains and carotenoids. Plant J. 2008; 54: 733-49

[2] Fang J . Bioavailability of anthocyanins. Drug Metab Rev. 2014; 46: 508-20

[3] Dixon R, Achnine L, Kota P et al. The phenylpropanoid pathway and plant defence-a genomics perspective. Mol Plant Pathol. 2002; 3: 371-90

[4] Naing A, Kim C . Abiotic stress-induced anthocyanins in plants: Their role in tolerance to abiotic stresses. Physiol Plant. 2021; 172: 1711-23

[5] Francavilla A, Joye I . Anthocyanins in whole grain cereals and their potential effect on health. Nutrients. 2020; 12: 2922

[6] Shi M, Xie D . Biosynthesis and metabolic engineering of anthocyanins in Arabidopsis thaliana. Recent Pat Biotechnol. 2014; 8: 47-60

[7] Xu W, Grain D, Bobet S et al. Complexity and robustness of the flavonoid transcriptional regulatory network revealed by comprehensive analyses of MYB-bHLH-WDR complexes and their targets in Arabidopsis seed. New Phytol. 2014; 202: 132-44

[8] Albert N, Davies KM, Lewis DH et al. A conserved network of transcriptional activators and repressors regulates anthocyanin pigmentation in eudicots. Plant Cell. 2014; 26: 962-80

[9] Gonzalez A, Zhao M, Leavitt J et al. Regulation of the anthocyanin biosynthetic pathway by the TTG1/bHLH/Myb transcriptional complex in Arabidopsis seedlings . Plant J. 2008; 53: 814-27

[10] Dubos C, le Gourrierec J, Baudry A et al. MYBL2 is a new regulator of flavonoid biosynthesis in Arabidopsis thaliana. Plant J. 2008; 55: 940-53

[11] Gómez-Campo C . Morphology and morpho-taxonomy of the tribe Brassiceae. In: Brassica Crops and Wild Allies. 1980, 3-31

[12] Liu S, Liu Y, Yang X et al. The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes . Nat Commun. 2014; 5: 1-11

[13] Bayer P, Golicz AA, Tirnaz S et al. Variation in abundance of predicted resistance genes in the Brassica oleracea pangenome . Plant Biotechnol J. 2019; 17: 789-800

[14] Zhang Y, Feng X, Liu Y et al. A single-base insertion in BoDFR1 results in loss of anthocyanins in green-leaved ornamental kale . Theor Appl Genet. 2022; 135: 1855-65

[15] Tang Q, Tian M, An G et al. Rapid identification of the purple stem (Ps) gene of Chinese kale (Brassica oleracea var. alboglabra) in a segregation distortion population by bulked segregant analysis and RNA sequencing . Mol Breed. 2017; 37: 1-9

[16] Liu X, Gao B-z, Han F-q et al. Genetics and fine mapping of a purple leaf gene, BoPr, in ornamental kale (Brassica oleracea L. var. acephala). BMC Genomics. 2017; 18: 1-9

[17] Guo N, Han S, Zong M et al. Identification and differential expression analysis of anthocyanin biosynthetic genes in leaf color variants of ornamental kale. BMC Genomics. 2019; 20: 1-17

[18] Han F, Zhang X, Yang L et al. Genome-wide characterization and analysis of the anthocyanin biosynthetic genes in Brassica oleracea. Planta. 2021; 254: 1-14

[19] Chiu L, Zhou X, Burke S et al. The purple cauliflower arises from activation of a MYB transcription factor. Plant Physiol. 2010; 154: 1470-80

[20] Yan C, An G, Zhu T et al. Independent activation of the BoMYB2 gene leading to purple traits in Brassica oleracea. Theor Appl Genet. 2019; 132: 895-906

[21] Song H, Yi H, Lee M et al. Purple Brassica oleracea var. capitata F. rubra is due to the loss of BoMYBL2-1 expression . BMC Plant Biol. 2018; 18: 1-16

[22] Zhu P, Cheng M, Feng X et al. Mapping of Pi, a gene conferring pink leaf in ornamental kale (Brassica oleracea L. var. acephala DC). Euphytica. 2016; 207: 377-85

[23] Feng X, Zhang Y, Wang H et al. The dihydroflavonol 4-reductase BoDFR1 drives anthocyanin accumulation in pink-leaved ornamental kale . Theor Appl Genet. 2021; 134: 159-69

[24] Zhang B, Hu Z, Zhang Y et al. A putative functional MYB transcription factor induced by low temperature regulates anthocyanin biosynthesis in purple kale (Brassica Oleracea var. acephala f. tricolor). Plant Cell Rep. 2012; 31: 281-9

[25] Ahn J, Jung Y, Song H et al. Alleles disrupting LBD37-like gene by an 136bp insertion show different distributions between green and purple cabbages (Brassica oleracea var. capitata). Genes & Genomics. 2021; 43: 679-88

[26] Rahim M, Afrin KS, Jung HJ et al. Molecular analysis of anthocyanin biosynthesis-related genes reveal BoTT8 associated with purple hypocotyl of broccoli (Brassica oleracea var. italica L.). Genome. 2019; 62: 253-66

[27] Parkin I, Koh C, Tang H et al. Transcriptome and methylome profiling reveals relics of genome dominance in the mesopolyploid Brassica oleracea. Genome Biol. 2014; 15: R77-18

[28] Yi D, Cui L, Wang L et al. Pyramiding of Bt cry1Ia8 and cry1Ba3 genes into cabbage (Brassica oleracea L. var. capitata) confers effective control against diamondback moth . Plant Cell Tissue Organ Cult. 2013; 115: 419-28

[29] Hong H, Phan A, O’Hare T . Temperature and maturity stages affect anthocyanin development and phenolic and sugar content of purple-pericarp supersweet sweetcorn during storage. J Agric Food Chem. 2021; 69: 922-31

[30] Lila M, Burton-Freeman B, Grace M et al. Unraveling anthocyanin bioavailability for human health. Annu Rev Food Sci Technol. 2016; 7: 375-93

[31] Zhang Y, Butelli E, Martin C . Engineering anthocyanin biosynthesis in plants. Curr Opin Plant Biol. 2014; 19: 81-90

[32] Song H, Park JI, Hwang BH et al. SNP in DFR1 coding sequence is tightly associated with anthocyanin accumulation in cabbage (B. oleracea var. capitata f. alba) at low temperature . Agronomy. 2020; 10: 602

[33] Zhang F, Wen Y, Guo X . CRISPR/Cas9 for genome editing: progress, implications and challenges. Hum Mol Genet. 2014; 23: R40-6

[34] Brooks C, Nekrasov V, Lippman Z et al. Efficient gene editing in tomato in the first generation using the clustered regularly interspaced short palindromic repeats/CRISPR-associated9 system. Plant Physiol. 2014; 166: 1292-7

[35] Ran F, Hsu PD, Wright J et al. Genome engineering using the CRISPR-Cas9 system. Nat Protoc. 2013; 8: 2281-308

[36] Gupta D, Bhattacharjee O, Mandal D et al. CRISPR-Cas9 system: a new-fangled dawn in gene editing. Life Sci. 2019; 232: 116636

[37] Mishra R, Joshi R, Zhao K . Genome editing in rice: recent advances, challenges, and future implications. Front Plant Sci. 2018; 9: 1361

[38] Rönspies M, Schindele P, Wetzel R et al. CRISPR-Cas9-mediated chromosome engineering in Arabidopsis thaliana. Nat Protoc. 2022; 17: 1332-58

[39] Agarwal A, Yadava P, Kumar K et al. Insights into maize genome editing via CRISPR/Cas9. Physiol Mol Biol Plants. 2018; 24: 175-83

[40] Sun B, Zheng A, Jiang M et al. CRISPR/Cas9-mediated mutagenesis of homologous genes in Chinese kale. Sci Rep. 2018; 8: 1-10

[41] LaManna L, Parulekar M, Maliga P . Multiple sgRNAs for onestep inactivation of the duplicated acetyl-coenzyme A carboxylase 2 (ACC2) genes in Brassica napus. Plant Physiol. 2022; 189: 178-87

[42] Han F, Yang C, Fang ZY et al. Inheritance and InDel markers closely linked to petal color gene (cpc-1) in Brassica oleracea. Mol Breed. 2015; 35: 1-8

[43] Martin M . Cutadapt removes adapter sequences from highthroughput sequencing reads. Embnet journal. 2011; 17: 10-2

[44] Li H, Durbin R . Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009; 25: 1754-60

[45] Kumar S, Stecher G, Tamura K . MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016; 33: 1870-4

[46] Gómez-Lobato M, Hasperué JH, Civello PM et al. Effect of 1-MCP on the expression of chlorophyll degrading genes during senescence of broccoli (Brassica oleracea L.). Sci Hortic. 2012; 144: 208-11

[47] Han F, Yuan K, Kong C et al. Fine mapping and candidate gene identification of the genic male-sterile gene ms3 in cabbage 51S . Theor Appl Genet. 2018; 131: 2651-61

[48] Ji J, Yang LM, Fang ZY et al. Recessive male sterility in cabbage (Brassica oleracea var. capitata) caused by loss of function of BoCYP704B1 due to the insertion of a LTR-retrotransposon . Theor Appl Genet. 2017; 130: 1441-51

[49] Livak K, Schmittgen T . Analysis of relative gene expression data using real-time quantitative PCR and the 2 ^−ΔΔCT method. Methods. 2001; 25: 402-8