The high-quality genome of pummelo provides insights into the tissue-specific regulation of citric acid and anthocyanin during domestication

Abstract

Citric acid and anthocyanin contents were co-selected during Citrus domestication. Pummelo is a founding species in the Citrus genus, but the domestication of pummelo has not been well studied. Here, we compared the citric acid and anthocyanin contents of a low citric acid pummelo (Citrus maxima LCA) and its high citric acid variety (HCA) from the same cultivation area in China. Our study revealed that, unlike the LCA type, the HCA variety accumulated anthocyanin in the pericarp early in fruit development. To investigate the genetic basis of acid and anthocyanin enrichment in HCA pulp and pericarp, respectively, we generated a chromosome-scale HCA genome using long-read sequence reads and Hi-C sequencing data. Transcriptome analysis and transient overexpression assays showed that the accumulation of citric acid and anthocyanin was associated with high expression of CgANTHOCYANIN1 (CgAN1), and two different MYBs transcription factors (CgPH4 and CgRuby1), respectively. Moreover, the CgAN1 promoter was more methylated in the LCA pulp than in the HCA pulp. Treatment with a DNA methylation inhibitor, 5-azacytidine, alleviated the CgAN1 promoter hypermethylation in the LCA pulp, leading to increased CgAN1 expression and citric acid content. This study provides a new high-quality pummelo genome and insight into the molecular mechanism behind the change in tissue-specific citric acid and anthocyanin accumulation during pummelo domestication and provides a conceptual basis for precise genetic manipulation in fruit flavor breeding.

References

[1] Butelli E, Licciardello C, Ramadugu C et al. Noemi controls production of flavonoid pigments and fruit acidity and illustrates the domestication routes of modern citrus varieties. Curr Biol. 2019; 29: 158-164.e2.

[2] Hu DG, Sun CH, Ma QJ et al. MdMYB1 regulates anthocyanin and malate accumulation by directly facilitating their transport into vacuoles in apples. Plant Physiol. 2016; 170: 1315-30.

[3] Strazzer P, Spelt CE, Li S et al. Hyperacidification of citrus fruits by a vacuolar proton-pumping P-ATPase complex. Nat Commun. 2019; 10: 744.

[4] Spelt C, Quattrocchio F, Mol JN et al. ANTHOCYANIN1 of petunia encodes a basic helix-loop-helix protein that directly activates transcription of structural ANTHOCYANIN genes. Plant Cell. 2000; 12: 1619-31.

[5] Huang D, Wang X, Tang Z et al. Subfunctionalization of the Ruby2-Ruby1 gene cluster during the domestication of citrus. Nat Plants. 2018; 4: 930-41.

[6] Spelt C, Quattrocchio F, Mol J et al. ANTHOCYANIN1 of petunia controls pigment synthesis, vacuolar pH, and seed coat development by genetically distinct mechanisms. Plant Cell. 2002; 14: 2121-35.

[7] Zhang Y, Ye J, Liu C et al. Citrus PH4-Noemi regulatory complex is involved in proanthocyanidin biosynthesis via a positive feedback loop. J Exp Bot. 2020; 71: 1306-21.

[8] He JX, Xu Y, Huang D et al. TRIPTYCHON-LIKE regulates aspects of both fruit flavor and color in citrus. J Exp Bot. 2022; 73: 3610-24.

[9] Telias A, Lin-Wang K, Stevenson DE et al. Apple skin patterning is associated with differential expression of MYB10. BMC Plant Biol. 2011; 11: 93.

[10] Jiang S, Wang N, Chen M et al. Methylation of MdMYB1 locus mediated by RdDM pathway regulates anthocyanin biosynthesis in apple. Plant Biotechnol J. 2020; 18: 1736-48.

[11] Jia D, Li Z, Dang Q et al. Anthocyanin biosynthesis and methylation of the MdMYB10 promoter are associated with the red blushed-skin mutant in the red striped-skin "Changfu 2" apple. J Agric Food Chem. 2020; 68: 4292-304.

[12] Shao SY, Xu WJ, Tao J et al. Glycemic index, glycemic load, and glycemic response to pomelo in patients with type 2 diabetes. J Huazhong Univ Sci Technol Med Sci. 2017; 37: 711-8.

[13] Wang X, Xu Y, Zhang S et al. Genomic analyses of primitive, wild and cultivated citrus provide insights into asexual reproduction. Nat Genet. 2017; 49: 765-72.

[14] Castillejo C, Waurich V, Wagner H et al. Allelic variation of MYB10 is the major force controlling natural variation in skin and flesh color in strawberry (Fragaria spp.). Fruit Plant Cell. 2020; 32: 3723-49.

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

[16] Chen J, Yuan Z, Zhang H et al. Cit1,2RhaT and two novel Cit-dGlcTs participate in flavor-related flavonoid metabolism during citrus fruit development. J Exp Bot. 2019; 70: 2759-71.

[17] Yu HW, Yang X, Guo F et al. Genetic diversity and population structure of pummelo (Citrus maxima) germplasm in China. Tree Genet Genomes. 2017; 13, 58.

[18] Butelli E, Licciardello C, Zhang Y et al. Retrotransposons control fruit-specific, cold-dependent accumulation of anthocyanins in blood oranges. Plant Cell. 2012; 24: 1242-55.

[19] Butelli E, Garcia-Lor A, Licciardello C et al. Changes in anthocyanin production during domestication of citrus. Plant Physiol. 2017; 173: 2225-42.

[20] Rodrigues JA, Espley RV, Allan AC . Genomic analysis uncovers functional variation in the C-terminus of anthocyanin-activating MYB transcription factors. Hortic Res. 2021; 8: 77.

[21] Quattrocchio F, Verweij W, Kroon A et al. PH4 of petunia is an R2R3 MYB protein that activates vacuolar acidification through interactions with basic-helix-loop-helix transcription factors of the anthocyanin pathway. Plant Cell. 2006; 18: 1274-91.

[22] Zhu F, Wen W, Fernie AR . Finding Noemi: the transcription factor mutations underlying trait differentiation amongst citrus. Trends Plant Sci. 2019; 24: 384-6.

[23] Song Q, Zhang T, Stelly DM et al. Epigenomic and functional analyses reveal roles of epialleles in the loss of photoperiod sensitivity during domestication of allotetraploid cottons. Genome Biol. 2017; 18: 99.

[24] Shi J, Lai J . Patterns of genomic changes with crop domestication and breeding. Curr Opin Plant Biol. 2015; 24: 47-53.

[25] Ding M, Chen ZJ . Epigenetic perspectives on the evolution and domestication of polyploid plant and crops. Curr Opin Plant Biol. 2018; 42: 37-48.

[26] Sicilia A, Scialo E, Puglisi I et al. Anthocyanin biosynthesis and DNA methylation dynamics in sweet orange fruit [Citrus sinensis L. (Osbeck)] under cold stress. J Agric Food Chem. 2020; 68: 7024-31.

[27] Wang Z, Meng D, Wang A et al. The methylation of the PcMYB10 promoter is associated with green-skinned sport in max red Bartlett pear. Plant Physiol. 2013; 162: 885-96.

[28] Manning K, Tör M, Poole M et al. A naturally occurring epigenetic mutation in a gene encoding an SBP-box transcription factor inhibits tomato fruit ripening. Nat Genet. 2006; 38: 948-52.

[29] Mancinelli AL, Rabino I . Photocontrol of anthocyanin synthesis: IV. Dose dependence and reciprocity relationships in anthocyanin synthesis. Plant Physiol. 1975; 56: 351-5.

[30] Liu Q, Xu J, Liu Y et al. A novel bud mutation that confers abnormal patterns of lycopene accumulation in sweet orange fruit (Citrus sinensis L. Osbeck). J Exp Bot. 2007; 58: 4161-71.

[31] Tan FQ, Tu H, Liang WJ et al. Comparative metabolic and transcriptional analysis of a doubled diploid and its diploid citrus rootstock (C. junos cv. Ziyang Xiangcheng) suggests its potential value for stress resistance improvement. BMC Plant Biol. 2015; 15: 89.

[32] Emms DM, Kelly S . OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol. 2015; 16: 157.

[33] Stamatakis A . RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014; 30: 1312-3.

[34] McKenna A, Hanna M, Banks E et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010; 20: 1297-303.

[35] Cingolani P, Platts A, Wang LL et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly. 2012; 6: 80-92.

[36] Kumar S, Stecher G, Li M et al. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018; 35: 1547-9.

[37] Purcell S, Neale B, Todd-Brown K et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007; 81: 559-75.

[38] Trapnell C, Williams BA, Pertea G et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010; 28: 511-5.

[39] Kim D, Pertea G, Trapnell C et al. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013; 14: R36.

[40] Chen C, Chen H, Zhang Y et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020; 13: 1194-202.

[41] Yu G, Wang LG, Han Y et al. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012; 16: 284-7.

[42] Jefferson RA, Kavanagh TA, Bevan MW . GUS fusions: beta-glucuronidase as a sensitive and versatile gene fusion marker in higher plants. EMBO J. 1987; 6: 3901-7.

[43] Lu S, Zhang Y, Zheng X et al. Isolation and functional characterization of a lycopene β-cyclase gene promoter from citrus. Front Plant Sci. 2016; 7: 1367.

[44] Huang H, Liu R, Niu Q et al. Global increase in DNA methylation during orange fruit development and ripening. Proc Natl Acad Sci USA. 2019; 116: 1430-6.

[45] Cheng Y, de Vicente MC, Meng H et al. A set of primers for analyzing chloroplast DNA diversity in citrus and related genera. Tree Physiol. 2005; 25: 661-72.