Multi-omics analyses reveal MdMYB10 hypermethylation being responsible for a bud sport of apple fruit color

Abstract

Apple bud sports offer a rich resource for clonal selection of numerous elite cultivars. The accumulation of somatic mutations as plants develop may potentially impact the emergence of bud sports. Previous studies focused on somatic mutation in the essential genes associated with bud sports. However, the rate and function of genome-wide somatic mutations that accumulate when a bud sport arises remain unclear. In this study, we identified a branch from a 10-year-old tree of the apple cultivar ‘Oregon Spur II’ as a bud sport. The mutant branch showed reduced red coloration on fruit skin. Using this plant material, we assembled a high-quality haplotype reference genome consisting of 649.61 Mb sequences with a contig N50 value of 2.04 Mb. We then estimated the somatic mutation rate of the apple tree to be 4.56 × 10⁻⁸ per base per year, and further identified 253 somatic single-nucleotide polymorphisms (SNPs), including five non-synonymous SNPs, between the original type and mutant samples. Transcriptome analyses showed that 69 differentially expressed genes between the original type and mutant fruit skin were highly correlated with anthocyanin content. DNA methylation in the promoter of five anthocyanin-associated genes was increased in the mutant compared with the original type as determined using DNA methylation profiling. Among the genetic and epigenetic factors that directly and indirectly influence anthocyanin content in the mutant apple fruit skin, the hypermethylated promoter of MdMYB10 is important. This study indicated that numerous somatic mutations accumulated at the emergence of a bud sport from a genome-wide perspective, some of which contribute to the low coloration of the bud sport.

References

[1] Deng XX, Wang LR, Li SH et al. Retrospection and prospect of fruit breeding for last four decades in China. J Fruit Sci. 2019; 36: 514-20.

[2] Foster TM, Aranzana MJ . Attention sports fans! The far-reaching contributions of bud sport mutants to horticulture and plant biology. Hortic Res. 2018; 5: 44.

[3] Zotta LA, Logan C . 200 Years and Growing: The Story of Stark Bro’s Nurseries & Orchards Co . Stark Bro’s Nurseries & Orchards Company: Louisiana, 2015.

[4] Xin P, Xiang F . Classification of and relation for the native, introduced and bred apple cultivars in China. J Shandong Agric Univ. 1998; 29: 189-200.

[5] Velasco R, Zharkikh A, Affourtit J et al. The genome of the domesticated apple (Malus × domestica Borkh.) . Nat Genet. 2010; 42: 833-9.

[6] Zhang L, Hu J, Han X et al. A high-quality apple genome assembly reveals the association of a retrotransposon and red fruit colour. Nat Commun. 2019; 10: 1494.

[7] Daccord N, Celton JM, Linsmith G et al. High-quality de novo assembly of the apple genome and methylome dynamics of early fruit development. Nat Genet. 2017; 49: 1099-106.

[8] Sun X, Jiao C, Schwaninger H et al. Phased diploid genome assemblies and pan-genomes provide insights into the genetic history of apple domestication. Nat Genet. 2020; 52: 1423-32.

[9] Muller HJ . The relation of recombination to mutational advance. Mutat Res. 1964; 106: 2-9.

[10] Wang L, Ji Y, Hu Y et al. The architecture of intra-organism mutation rate variation in plants. PLoS Biol. 2019; 17: 1-29.

[11] Wang L, Huang Y, Liu Z et al. Somatic variations led to the selection of acidic and acidless orange cultivars. Nature Plants. 2021; 7: 954-65.

[12] Kobayashi S, Goto-Yamamoto N, Hirochika H . Retrotransposon-induced mutations in grape skin color. Science. 2004; 304: 982.

[13] El-Sharkawy I, Liang D, Xu K . Transcriptome analysis of an apple (Malus × domestica) yellow fruit somatic mutation identifies a gene network module highly associated with anthocyanin and epigenetic regulation. J Exp Bot. 2015; 66: 7359-76.

[14] Xu Y, Feng S, Jiao Q et al. Comparison of MdMYB1 sequences and expression of anthocyanin biosynthetic and regulatory genes between Malus domestica Borkh. cultivar ‘Ralls’ and its blushed sport . Euphytica. 2012; 185: 157-70.

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

[16] Qian M, Sun Y, Allan AC et al. The red sport of ‘Zaosu’ pear and its red-striped pigmentation pattern are associated with demethylation of the PyMYB10 promoter. Phytochemistry. 2014; 107: 16-23.

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

[18] Yang T, Ma H, Zhang J et al. Systematic identification of long noncoding RNAs expressed during light-induced anthocyanin accumulation in apple fruit. Plant Biotechnol J. 2019; 100: 572-90.

[19] Winkel-Shirley B. Flavonoid biosynthesis. A colorful model for genetics, biochemistry, cell biology, and biotechnology. Plant Physiol. 2001; 126: 485-93.

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

[21] Fischer TC, Halbwirth H, Meisel B et al. Molecular cloning, substrate specificity of the functionally expressed dihydroflavonol 4-reductases from Malus domestica and Pyrus communis cultivars and the consequences for flavonoid metabolism . Arch Biochem Biophys. 2003; 412: 223-30.

[22] Chen Z, Yu L, Liu W et al. Research progress of fruit color development in apple (Malus domestica Borkh.) . Plant Physiol Biochem. 2021; 162: 267-79.

[23] Tao R, Yu W, Gao Y et al. Light-induced basic/helix-loop-helix64 enhances anthocyanin biosynthesis and undergoes CONSTITUTIVELY PHOTOMORPHOGENIC1-mediated degradation in pear. Plant Physiol. 2020; 184: 1684-701.

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

[25] Mehrtens F, Kranz H, Bednarek P et al. The Arabidopsis transcription factor MYB12 is a flavonol-specific regulator of phenylpropanoid biosynthesis. Plant Physiol. 2005; 138: 1083-96.

[26] Stracke R, Ishihara H, Huep G et al. Differential regulation of closely related R2R3-MYB transcription factors controls flavonol accumulation in different parts of the Arabidopsis thaliana seedling. Plant J. 2007; 50: 660-77.

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

[28] Ban Y, Honda C, Hatsuyama Y et al. Isolation and functional analysis of a MYB transcription factor gene that is a key regulator for the development of red coloration in apple skin. Plant Cell Physiol. 2007; 48: 958-70.

[29] 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-0.

[30] Jiang S, Chen M, He N et al. MdGSTF6, activated by MdMYB1, plays an essential role in anthocyanin accumulation in apple. Hortic Res. 2019; 6: 40.

[31] Allan AC, Hellens RP, Laing WA . MYB transcription factors that colour our fruit. Trends Plant Sci. 2008; 13: 99-102.

[32] Jiang SH, Sun QG, Chen M et al. Methylome and transcriptome analyses of apple fruit somatic mutations reveal the difference of red phenotype. BMC Genomics. 2019; 20: 117.

[33] Tian Y, Thrimawithana A, Ding T et al. Transposon insertions regulate genome-wide allele-specific expression and underpin flower colour variations in apple (Malus spp.) . Plant Biotechnol J. 2022.

[34] Migicovsky Z, Gardner KM, Richards C et al. Genomic consequences of apple improvement. Hortic Res. 2021; 8: 9.

[35] Vondras AM, Minio A, Blanco-Ulate B et al. The genomic diversification of grapevine clones. BMC Genomics. 2019; 20: 972.

[36] Schoen DJ, Schultz ST . Somatic mutation and evolution in plants. Annu Rev Ecol Evol Syst. 2019; 50: 49-73.

[37] Cheng J, Liao L, Zhou H et al. A small indel mutation in an anthocyanin transporter causes variegated colouration of peach flowers. J Exp Bot. 2015; 66: 7227-39.

[38] Chen P, Li Z, Zhang D et al. Insights into the effect of human civilization on Malus evolution and domestication. Plant Biotechnol J. 2021; 19: 2206-20.

[39] An JP, Zhang XW, Bi SQ et al. The ERF transcription factor MdERF38 promotes drought stress-induced anthocyanin biosynthesis in apple. Plant J. 2020; 101: 573-89.

[40] Xing L, Zhang D, Song X et al. Genome-wide sequence variation identification and floral-associated trait comparisons based on the re-sequencing of the ’Nagafu No. 2′ and ’Qinguan’ varieties of apple (Malus domestica Borkh.) . Front Plant Sci. 2016; 7: 908-8.

[41] Jiao F, Luo R, Dai X et al. Chromosome-level reference genome and population genomic analysis provide insights into the evolution and improvement of domesticated mulberry (Morus alba) . Mol Plant. 2020; 13: 1001-12.

[42] Gonçalves B, Silva AP, Moutinho-Pereira J et al. Effect of ripeness and postharvest storage on the evolution of colour and anthocyanins in cherries (Prunus avium L.) . Food Chem. 2007; 103: 976-84.

[43] Zhou JY, Sun CD, Zhang LL et al. Preferential accumulation of orange-colored carotenoids in Ponkan (Citrus reticulata) fruit peel following postharvest application of ethylene or ethephon . Sci Hortic. 2010; 126: 229-35.

[44] Karagiannis E, Michailidis M, Tanou G et al. Decoding altitude-activated regulatory mechanisms occurring during apple peel ripening. Hortic Res. 2020; 7: 120.

[45] An JP, Qu FJ, Yao JF et al. The bZIP transcription factor MdHY5 regulates anthocyanin accumulation and nitrate assimilation in apple. Hortic Res. 2017; 4: 17023.

[46] Ranjbar S, Ramezanian A, Rahemi M . Nano-calcium and its potential to improve ‘Red Delicious’ apple fruit characteristics. Hortic Environ Biotechnol. 2020; 61: 23-30.

[47] Chen Y, Nie F, Xie SQ et al. Efficient assembly of nanopore reads via highly accurate and intact error correction. Nat Commun. 2021; 12: 60.

[48] Hu J, Fan J, Sun Z et al. NextPolish: a fast and efficient genome polishing tool for long-read assembly. Bioinformatics. 2020; 36: 2253-5.

[49] Guan D, McCarthy SA, Wood J et al. Identifying and removing haplotypic duplication in primary genome assemblies. Bioinformatics. 2020; 36: 2896-8.

[50] Alonge M, Soyk S, Ramakrishnan S et al. RaGOO: fast and accurate reference-guided scaffolding of draft genomes. Genome Biol. 2019; 20: 224.

[51] Tang H, Zhang X, Miao C et al. ALLMAPS: robust scaffold ordering based on multiple maps. Genome Biol. 2015; 16: 3.

[52] Howard NP, Troggio M, Durel CE et al. Integration of Infinium and Axiom SNP array data in the outcrossing species Malus × domestica and causes for seemingly incompatible calls . BMC Genomics. 2021; 22: 246.

[53] Tarailo-Graovac M, Chen N . Using RepeatMasker to identify repetitive elements in genomic sequences. Curr Protoc Bioinformatics. 2009; 25: 4.10.1-14.

[54] Yandell M, Ence D . A beginner’s guide to eukaryotic genome annotation. Nat Rev Genet. 2012; 13: 329-42.

[55] Pratas MI, Aguiar B, Vieira J et al. Inferences on specificity recognition at the Malus × domestica gametophytic self-incompatibility system . Sci Rep. 2018; 8: 1717.

[56] Zhu L, Ni W, Liu S et al. Transcriptomics analysis of apple leaves in response to Alternaria alternata apple pathotype infection. Front Plant Sci. 2017; 8: 22-2.

[57] Stanke M, Keller O, Gunduz I et al. AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res. 2006; 34: W435-9.

[58] Brůna T, Hoff KJ, Lomsadze A et al. BRAKER2: automatic eukaryotic genome annotation with GeneMark-EP+ and AUGUSTUS supported by a protein database. NAR Genom Bioinform. 2021; 3: lqaa108.

[59] Kim D, Langmead B, Salzberg SL . HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015; 12: 357-60.

[60] Pertea M, Pertea GM, Antonescu CM et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015; 33: 290-5.

[61] Simão FA, Waterhouse RM, Ioannidis P et al. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015; 31: 3210-2.

[62] Quevillon E, Silventoinen V, Pillai S et al. InterProScan: protein domains identifier. Nucleic Acids Res. 2005; 33: W116-20.

[63] Duan N, Bai Y, Sun H et al. Genome re-sequencing reveals the history of apple and supports a two-stage model for fruit enlargement. Nat Commun. 2017; 8: 249.

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

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

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

[67] Chen X, Schulz-Trieglaff O, Shaw R et al. Manta: rapid detection of structural variants and indels for germline and cancer sequencing applications. Bioinformatics. 2016; 32: 1220-2.

[68] Layer RM, Chiang C, Quinlan AR et al. LUMPY: a probabilistic framework for structural variant discovery. Genome Biol. 2014; 15: R84.

[69] Rausch T, Zichner T, Schlattl A et al. DELLY: structural variant discovery by integrated paired-end and split-read analysis. Bioinformatics. 2012; 28: i333-9.

[70] Jeffares DC, Jolly C, Hoti M et al. Transient structural variations have strong effects on quantitative traits and reproductive isolation in fission yeast. Nat Commun. 2017; 8: 14061.

[71] Sedlazeck FJ, Rescheneder P, Smolka M et al. Accurate detection of complex structural variations using single-molecule sequencing. Nat Methods. 2018; 15: 461-8.

[72] Robinson JT, Thorvaldsdóttir H, Winckler W et al. Integrative genomics viewer. Nat Biotechnol. 2011; 29: 24-6.

[73] Chang CC, Chow CC, Tellier LC et al. Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience. 2015; 4: 7.

[74] Manichaikul A, Mychaleckyj J, Rich S et al. Robust relationship inference in genome wide association studies. Bioinformatics. 2010; 26: 2867-73.

[75] Retief JD . Phylogenetic analysis using PHYLIP. In: Misener S, Krawetz SA, eds. Bioinformatics Methods and Protocols.Humana Press: New Jersey, 243-58.

[76] Liao Y, Smyth GK, Shi W . featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014; 30: 923- 30.

[77] Love MI, Huber W, Anders S . Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15: 550-0.

[78] Langfelder P, Horvath S . WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 2008; 9: 559-9.

[79] Saito R, Smoot ME, Ono K et al. A travel guide to Cytoscape plugins. Nat Methods. 2012; 9: 1069-76.

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

[81] Krueger F, Andrews SR . Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 2011; 27: 1571-2.

[82] Guo W, Zhu P, Pellegrini M et al. CGmapTools improves the precision of heterozygous SNV calls and supports allele-specific methylation detection and visualization in bisulfite-sequencing data. Bioinformatics. 2018; 34: 381-7.

[83] Gaspar JM, Hart RP . DMRfinder: efficiently identifying differentially methylated regions from MethylC-seq data. BMC Bioinformatics. 2017; 18: 528.

[84] Yu G, Wang LG, He QY . ChIPseeker: an R/Bioconductor package for ChIP peak annotation, comparison and visualization. Bioinformatics. 2015; 31: 2382-3.