An LTR retrotransposon in the promoter of a PsMYB10.2 gene associated with the regulation of fruit flesh color in Japanese plum

Abstract

Japanese plums exhibit wide diversity of fruit coloration. The red to black hues are caused by the accumulation of anthocyanins, while their absence results in yellow, orange or green fruits. In Prunus, MYB10 genes are determinants for anthocyanin accumulation. In peach, QTLs for red plant organ traits map in an LG3 region with three MYB10 copies (PpMYB10.1, PpMYB10.2 and PpMYB10.3). In Japanese plum the gene copy number in this region differs with respect to peach: there are at least three copies of PsMYB10.1, with the expression of one of them (PsMYB10.1a) correlating with fruit skin color. The objective of this study was to determine a possible role of LG3-PsMYB10 genes in the natural variability of the flesh color trait and to develop a molecular marker for marker-assisted selection (MAS). We explored the variability within the LG3-PsMYB10 region using long-range sequences obtained in previous studies through CRISPR-Cas9 enrichment sequencing. We found that the PsMYB10.2 gene was only expressed in red flesh fruits. Its role in promoting anthocyanin biosynthesis was validated by transient overexpression in Japanese plum fruits. The analysis of long-range sequences identified an LTR retrotransposon in the promoter of the expressed PsMYB10.2 gene that explained the trait in 93.1% of the 145 individuals analyzed. We hypothesize that the LTR retrotransposon may promote the PsMYB10.2 expression and activate the anthocyanin biosynthesis pathway. We propose for the first time the PsMYB10.2 gene as candidate for the flesh color natural variation in Japanese plum and provide a molecular marker for MAS.

References

[1] Saigo T, Wang T, Watanabe M et al. Diversity of anthocyanin and proanthocyanin biosynthesis in land plants. Curr Opin Plant Biol. 2020; 55: 93-9.

[2] Winkel-Shirley B . Biosynthesis of flavonoids and effects of stress. Curr Opin Plant Biol. 2002; 5: 218-23.

[3] Khoo HE, Azlan A, Tang ST et al. Anthocyanidins and anthocyanins: colored pigments as food, pharmaceutical ingredients, and the potential health benefits. Food Nutr Res. 2017; 61: 1361779.

[4] Espley RV, Hellens RP, Putterill J et al. Red colouration in apple fruit is due to the activity of the MYB transcription factor, MdMYB10. Plant J. 2007; 49: 414-27.

[5] 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.

[6] Tuan PA, Bai S, Yaegaki H et al. The crucial role of PpMYB10.1 in anthocyanin accumulation in peach and relationships between its allelic type and skin color phenotype. BMC Plant Biol. 2015; 15: 280.

[7] Feng S, Wang Y, Yang S et al. Anthocyanin biosynthesis in pears is regulated by a R2R3-MYB transcription factor PyMYB10 . Planta. 2010; 232: 245-55.

[8] Jin W, Wang H, Li M et al. The R2R3 MYB transcription factor PavMYB10.1 involves in anthocyanin biosynthesis and determines fruit skin colour in sweet cherry (Prunus avium L.) . Plant Biotechnol J. 2016; 14: 2120-33.

[9] Zhou H, Lin-Wang K, Wang H et al. Molecular genetics of blood-fleshed peach reveals activation of anthocyanin biosynthesis by NAC transcription factors. Plant J. 2015;82:105-21.

[10] Fiol A, Howad W, Surya A et al. Development of molecular markers for fruit skin color in Japanese plum (Prunus salicina Lindl.) . Acta Hortic. 2021; 1307: 221-6.

[11] Albert NW, 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.

[12] Stracke R, Werber M, Weisshaar B . The R2R3-MYB gene family in Arabidopsis thaliana. Curr Opin Plant Biol. 2001; 4: 447-56.

[13] Lin-Wang K, Bolitho K, Grafton K et al. An R2R3 MYB transcription factor associated with regulation of the anthocyanin biosynthetic pathway in Rosaceae. BMC Plant Biol. 2010; 10: 50.

[14] Koes R, Verweij W, Quattrocchio F . Flavonoids: a colorful model for the regulation and evolution of biochemical pathways. Trends Plant Sci. 2005; 10: 236-42.

[15] Heppel SC, Jaffé FW, Takos AM et al. Identification of key amino acids for the evolution of promoter target specificity of anthocyanin and proanthocyanidin regulating MYB factors. Plant Mol Biol. 2013; 82: 457-71.

[16] Lai Y, Li H, Yamagishi M . A review of target gene specificity of flavonoid R2R3-MYB transcription factors and a discussion of factors contributing to the target gene selectivity. Front Biol. 2013; 8: 577-98.

[17] Ravaglia D, Espley RV, Henry-Kirk RA et al. Transcriptional regulation of flavonoid biosynthesis in nectarine (Prunus persica) by a set of R2R3 MYB transcription factors . BMC Plant Biol. 2013; 13: 68-8.

[18] Zhou H, Lin-Wang K, Wang F et al. Activator-type R2R3-MYB genes induce a repressor-type R2R3-MYB gene to balance anthocyanin and proanthocyanidin accumulation. New Phytol. 2019; 221: 1919-34.

[19] Yan H, Pei X, Zhang H et al. MYB-mediated regulation of anthocyanin biosynthesis. Int J Mol Sci. 2021; 22: 3103.

[20] Espley RV, Brendolise C, Chagné D et al. Multiple repeats of a promoter segment causes transcription factor autoregulation in red apples. Plant Cell. 2009; 21: 168-83.

[21] 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-3.

[22] Pierantoni L, Dondini L, de Franceschi P et al. Mapping of an anthocyanin-regulating MYB transcription factor and its expression in red and green pear, Pyrus communis . Plant Physiol Bioch. 2010; 48: 1020-6.

[23] Bretó M, Cantín C, Iglesias I et al. Mapping a major gene for red skin color suppression (highlighter) in peach. Euphytica. 2017; 213: 14.

[24] Starkevič P, Paukštytė J, Kazanavičiūtė V et al. Expression and anthocyanin biosynthesis-modulating potential of sweet cherry (Prunus avium L.) MYB10 and bHLH genes. PLoS One. 2015; 10: e0126991.

[25] Fiol A, García-Gómez BE, Jurado-Ruiz F et al. Characterization of Japanese plum (Prunus salicina) PsMYB10 alleles reveals structural variation and polymorphisms correlating with fruit skin color . Front Plant Sci. 2021; 12: 655267.

[26] Zhang L et al. A high-quality apple genome assembly reveals the association of a retrotransposon and red fruit colour. Nat Commun. 2019; 10: 1-13.

[27] Cheng Y, Liu L, Yuan C et al. Molecular characterization of ethylene-regulated anthocyanin biosynthesis in plums during fruit ripening. Plant Mol Biol Rep. 2016; 34: 777-85.

[28] Takos AM, Jaffé FW, Jacob SR et al. Light-induced expression of a MYB gene regulates anthocyanin biosynthesis in red apples. Plant Physiol. 2006; 142: 1216-32.

[29] Niu J, Zhang G, Zhang W et al. Anthocyanin concentration depends on the counterbalance between its synthesis and degradation in plum fruit at high temperature. Sci Rep. 2017; 7: 1-16.

[30] Fang Z, Lin-Wang K, Jiang C et al. Postharvest temperature and light treatments induce anthocyanin accumulation in peel of ‘Akihime’ plum (Prunus salicina Lindl.) via transcription factor PsMYB10.1 . Postharvest Biol Technol. 2021; 179: 111592.

[31] Sooriyapathirana SS, Khan A, Sebolt AM et al. QTL analysis and candidate gene mapping for skin and flesh color in sweet cherry fruit (Prunus avium L.) . Tree Genet Genomes. 2010; 6: 821-32.

[32] Donoso JM, Picañol R, Serra O et al. Exploring almond genetic variability useful for peach improvement: mapping major genes and QTLs in two interspecific almond x peach populations. Mol Breeding. 2016; 36: 1-17.

[33] Salazar JA, Pacheco I, Shinya P et al. Genotyping by sequencing for SNP-based linkage analysis and identification of QTLs linked to fruit quality traits in Japanese plum (Prunus salicina Lindl.) . Front Plant Sci. 2017; 8: 476.

[34] Rahim MA, Busatto N, Trainotti L . Regulation of anthocyanin biosynthesis in peach fruits. Planta. 2014; 240: 913-29.

[35] Shirasawa K, Isuzugawa K, Ikenaga M et al. The genome sequence of sweet cherry (Prunus avium) for use in genomics-assisted breeding . DNA Res. 2017; 24: 499-508.

[36] Jiang F, Zhang J, Wang S et al. The apricot (Prunus armeniaca L.) genome elucidates Rosaceae evolution and beta-carotenoid synthesis . Hortic Res. 2019; 6: 128.

[37] Huang Z, Shen F, Chen Y et al. Chromosome-scale genome assembly and population genomics provide insights into the adaptation, domestication, and flavonoid metabolism of Chinese plum. Plant J. 2021; 108: 1174-92.

[38] Werner DJ, Okie W . A history and description of the Prunus persica: plant introduction collection. HortScience. 1998; 33: 787-93.

[39] Shen Z, Confolent C, Lambert P et al. Characterization and genetic mapping of a new blood-flesh trait controlled by the single dominant locus DBF in peach. Tree Genet Genomes. 2013; 9: 1435-46.

[40] Yamamoto T, Yamaguchi M, Hayashi T . An integrated genetic linkage map of peach by SSR, STS, AFLP and RAPD. J Jpn Soc Hort Sci. 2005; 74: 204-13.

[41] Beckman T, Alcazar JR, Sherman W et al. Evidence for qualitative suppression of red skin color in peach. HortScience. 2005; 40: 523-4.

[42] Guo J, Cao K, Deng C et al. An integrated peach genome structural variation map uncovers genes associated with fruit traits. Genome Biol. 2020; 21: 1-19.

[43] Salazar JA, Pacheco I, Silva C et al. Development and applicability of GBS approach for genomic studies in Japanese plum (Prunus salicina Lindl.) . J Hortic Sci Biotechnol. 2019; 94: 284-94.

[44] Acuña CV et al. Characterization of genetic diversity in accessions of prunus salicina lindl: keeping fruit flesh color ideotype while adapting to water stressed environments. Agronomy. 2019; 9: 487.

[45] Fiol A, Jurado-Ruiz F, Lopez-Girona E et al. An efficient CRISPR-Cas9 enrichment sequencing strategy for characterizing complex and highly duplicated genomic regions. A case study in the Prunus salicina LG3-MYB10 genes cluster. Plant Methods. 2022; 18: 105.

[46] Zhou H, Liao L, Xu S et al. Two amino acid changes in the R3 repeat cause functional divergence of two clustered MYB10 genes in peach. Plant Mol Biol. 2018; 98: 169-83.

[47] Yang J, Song HD, Chen Y et al. A single amino acid substitution in the R2R3 conserved domain of the BrPAP1a transcription factor impairs anthocyanin production in turnip (Brassica rapa subsp. rapa) . Plant Physiol Biochem. 2021; 162: 124-36.

[48] Verde I, Jenkins J, Dondini L et al. The peach v2.0 release: high-resolution linkage mapping and deep resequencing improve chromosome-scale assembly and contiguity. BMC Genomics. 2017; 18: 225.

[49] Liu C, Feng C, Peng W et al. Chromosome-level draft genome of a diploid plum (Prunus salicina) . GigaScience. 2020; 9: giaa130.

[50] N’Diaye A, Haile JK, Cory AT et al. Single marker and haplotype-based association analysis of semolina and pasta colour in elite durum wheat breeding lines using a high-density consensus map. PLoS One. 2017; 12: e0170941.

[51] Maldonado C, Mora F, Scapim CA et al. Genome-wide haplotype-based association analysis of key traits of plant lodging and architecture of maize identifies major determinants for leaf angle: hap LA4. PLoS One. 2019; 14: e0212925.

[52] Calle A, Wünsch A . Multiple-population QTL mapping of maturity and fruit-quality traits reveals LG4 region as a breeding target in sweet cherry (Prunus avium L.) . Hortic Res. 2020; 7: 1-13.

[53] Fang Z-Z, Lin-Wang K, Zhou D-R et al. Activation of PsMYB10.2 transcription causes anthocyanin accumulation in flesh of the red-fleshed mutant of ‘Sanyueli’ (Prunus salicina Lindl.) . Front Plant Sci. 2021; 12: 1167.

[54] Fang Z-Z, Zhou D-R, Ye X-F et al. Identification of candidate anthocyanin-related genes by transcriptomic analysis of ‘Furongli’ plum (Prunus salicina Lindl.) during fruit ripening using RNA-seq . Front Plant Sci. 2016; 7: 1338.

[55] Amarasinghe SL, Su S, Dong X et al. Opportunities and challenges in long-read sequencing data analysis. Genome Biol. 2020; 21: 1-16.

[56] Wei L, Cao X . The effect of transposable elements on phenotypic variation: insights from plants to humans. Sci China Life Sci. 2016; 59: 24-37.

[57] Jung S, Venkatesh J, Kang M-Y et al. A non-LTR retrotransposon activates anthocyanin biosynthesis by regulating a MYB transcription factor in Capsicum annuum . Plant Sci. 2019; 287: 110181.

[58] 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.

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

[60] Lallemand T, Leduc M, Landès C et al. An overview of duplicated gene detection methods: why the duplication mechanism has to be accounted for in their choice. Genes. 2020; 11: 1046.

[61] Otto SP, Yong P . The evolution of gene duplicates. Adv Genet. 2002; 46: 451-83.

[62] Prince VE, Pickett FB . Splitting pairs: the diverging fates of duplicated genes. Nat Rev Genet. 2002; 3: 827-37.

[63] Okie W. Prunus domestica (European Plum)- Prunus salicina (Japanese plum) . In: Janick J, Paull RE, eds. The Enciclopedia of Fruit and Nuts.Cambridge University press, 2008, 694-705.

[64] Salazar JA, Pacheco I, Zapata P et al. Identification of loci controlling phenology, fruit quality and post-harvest quantitative parameters in Japanese plum (Prunus salicina Lindl.) . Postharvest Biol Technol. 2020; 169: 111292.

[65] Valderrama-Soto D, Salazar J, Sepúlveda-González A et al. Detection of quantitative trait loci controlling the content of phenolic compounds in an Asian plum (Prunus salicina L.) F1 population . Front Plant Sci. 2021; 12: 1331.

[66] Crisosto CH, Garner D, Crisosto GM et al. Increasing ‘Blackamber’ plum (Prunus salicina Lindell) consumer acceptance . Postharvest Biol Technol. 2004; 34: 237-44.

[67] Manganaris GA, Vicente AR, Crisosto CH et al. Effect of delayed storage and continuous ethylene exposure on flesh reddening of ‘Royal Diamond’ plums. J Sci Food Agric. 2008; 88: 2180-5.

[68] Wang L, Sang W, Xu R et al. Alteration of flesh color and enhancement of bioactive substances via the stimulation of anthocyanin biosynthesis in ‘friar’ plum fruit by low temperature and the removal. Food Chem. 2020; 310: 125862.

[69] Doyle JJ, Doyle JL . A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochemistry. 1987; 19: 11-15.

[70] Kim H-Y et al. RNA-seq analysis of spatiotemporal gene expression patterns during fruit development revealed reference genes for transcript normalization in plums. Plant Mol Biol Rep. 2015; 33: 1634-49.

[71] Artimo P, Jonnalagedda M, Arnold K et al. ExPASy: SIB bioinformatics resource portal. Nucleic Acids Res. 2012; 40: W597-603.

[72] Sievers F, Wilm A, Dineen D et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal omega. Mol Syst Biol. 2011; 7: 539.

[73] Waterhouse AM, Procter JB, Martin DM et al. Jalview version 2-a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009; 25: 1189-91.

[74] Canli F, Tian L . Assessment of regeneration and transient expression factors for agrobacterium-mediated transformation of Prunus salicina Lindl. Eur J Hortic Sci. 2009; 74: 66.

[75] Xi W, Feng J, Liu Y et al. The R2R3-MYB transcription factor PaMYB10 is involved in anthocyanin biosynthesis in apricots and determines red blushed skin. BMC Plant Biol. 2019; 19: 287.

[76] Spolaore S, Trainotti L, Casadoro G . A simple protocol for transient gene expression in ripe fleshy fruit mediated by agrobacterium. J Exp Bot. 2001; 52: 845-50.

[77] Zhao Y, Mao W, Chen Y et al. Optimization and standardization of transient expression assays for gene functional analyses in strawberry fruits. Hortic Res. 2019; 6: 1-13.

[78] Altschul SF, Gish W, Miller W et al. Basic local alignment search tool. J Mol Biol. 1990; 215: 403-10.

[79] Ma K, Zhang Q, Cheng T et al. Identification of transposons near predicted lncRNA and mRNA pools of Prunus mume using an integrative transposable element database constructed from Rosaceae plant genomes. Mol Gen Genomics. 2018; 293: 1301-16.

[80] Xu Z, Wang H . LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Res. 2007; 35: W265-8.

[81] Kõressaar T, Lepamets M, Kaplinski L et al. Primer3_masker: integrating masking of template sequence with primer design software. Bioinformatics. 2018; 34: 1937-8.