A chromosome-scale genome sequence of pitaya (Hylocereus undatus) provides novel insights into the genome evolution and regulation of betalain biosynthesis

Abstract

Pitaya (Hylocereus) is the most economically important fleshy-fruited tree of the Cactaceae family that is grown worldwide, and it has attracted significant attention because of its betalain-abundant fruits. Nonetheless, the lack of a pitaya reference genome significantly hinders studies focused on its evolution, as well as the potential for genetic improvement of this crop. Herein, we employed various sequencing approaches, namely, PacBio-SMRT, Illumina HiSeq paired-end, 10× Genomics, and Hi-C (high-throughput chromosome conformation capture) to provide a chromosome-level genomic assembly of ‘GHB’ pitaya (H. undatus, 2n = 2x = 22 chromosomes). The size of the assembled pitaya genome was 1.41 Gb, with a scaffold N50 of ~127.15 Mb. In total, 27,753 protein-coding genes and 896.31 Mb of repetitive sequences in the H. undatus genome were annotated. Pitaya has undergone a WGT (whole-genome triplication), and a recent WGD (whole-genome duplication) occurred after the gamma event, which is common to the other species in Cactaceae. A total of 29,328 intact LTR-RTs (~696.45 Mb) were obtained in H. undatus, of which two significantly expanded lineages, Ty1/copia and Ty3/gypsy, were the main drivers of the expanded genome. A high-density genetic map of F1 hybrid populations of ‘GHB’ × ‘Dahong’ pitayas (H. monacanthus) and their parents were constructed, and a total of 20,872 bin markers were identified (56,380 SNPs) for 11 linkage groups. More importantly, through transcriptomic and WGCNA (weighted gene coexpression network analysis), a global view of the gene regulatory network, including structural genes and the transcription factors involved in pitaya fruit betalain biosynthesis, was presented. Our data present a valuable resource for facilitating molecular breeding programs of pitaya and shed novel light on its genomic evolution, as well as the modulation of betalain biosynthesis in edible fruits.

References

[1] Bellec, F. L., Vaillant, F. & Imbert, E. Pitahaya (Hylocereus spp.), a new fruit crop, a market with a future. Fruits 61, 237-250 (2006).

[2] Ibrahim, S. R. M., Mohamed, G. A., Khedr, A. I. M., Zayed, M. F. & El-Kholy, A. A. S. Genus Hylocereus: beneficial phytochemicals, nutritional importance, and biological relevance-A review. J. Food Biochem. 42, e12491 (2018).

[3] Lee, E. J. et al. Betalain and betaine composition of greenhouse- or field-produced beetroot (Beta vulgaris L.) and inhibition of HepG2 cell proliferation. J. Agric. Food Chem. 62, 1324-1331 (2014).

[4] Suh, D. H. et al. Metabolite profiling of red and white pitayas (Hylocereus polyrhizus and Hylocereus undatus) for comparing betalain biosynthesis and antioxidant activity. J. Agric. Food Chem. 62, 8764-8771 (2014).

[5] Teixeira Da Silva, D. V. et al. Betanin, a natural food additive: stability, bioavailability, antioxidant and preservative ability assessments. Molecules 24, 458 (2019).

[6] Mizrahi, Y., Nerd, A. & Nobel, P. S. Cacti as crops. Hortic. Rev. 18, 291-319 (1997).

[7] Mizrahi, Y., Mouyal, J., Nerd, A. & Sitrit, Y. Metaxenia in the vine cacti Hylocereus polyrhizus and Selenicereus spp. Ann. Bot. 93, 469-472 (2004).

[8] Tel-Zur, N. et al. Phenotypic and genomic characterization of vine cactus collection (Cactaceae). Genet. Res. Crop Evol. 58, 1075-1085 (2011).

[9] Copetti, D. et al. Extensive gene tree discordance and hemiplasy shaped the genomes of North American columnar cacti. Proc. Natl Acad. Sci. USA 114, 12003-12008 (2017).

[10] Zheng, J. et al. The chromosome-level genome of dragon fruit reveals whole-genome duplication and chromosomal co-localization of betacyanin biosynthetic genes. Hortic. Res. 8, 63 (2021).

[11] Gandia-Herrero, F. & Garcia-Carmona, F. The dawn of betalains. N. Phytol. 227, 664-666 (2020).

[12] Khan, M. I. & Giridhar, P. Plant betalains: chemistry and biochemistry. Phytochem 117, 267-295 (2015).

[13] Gandia-Herrero, F., Escribano, J. & Garcia-Carmona, F. Biological activities of plant pigments betalains. Crit. Rev. Food Sci. Nutr. 56, 937-945 (2016).

[14] Fernandez-Lopez, J. A., Fernandez-Lledo, V. & Angosto, J. M. New insights into red plant pigments: more than just natural colorants. RSC Adv. 10, 24669-24682 (2020).

[15] Polturak, G. & Aharoni, A. Advances and future directions in betalain metabolic engineering. N. Phytol. 224, 1472-1478 (2019).

[16] Lopez-Nieves, S. et al. Relaxation of tyrosine pathway regulation underlies the evolution of betalain pigmentation in Caryophyllales. N. Phytol. 217, 896-908 (2018).

[17] Gandía-Herrero, F. & García-Carmona, F. Characterization of recombinant Beta vulgaris 4,5-DOPA-extradiol-dioxygenase active in the biosynthesis of betalains. Planta 236, 91-100 (2012).

[18] Hatlestad, G. J. et al. The beet R locus encodes a new cytochrome P450 required for red betalain production. Nat. Genet. 44, 816-820 (2012).

[19] Tanaka, Y., Sasaki, N. & Ohmiya, A. Biosynthesis of plant pigments: anthocyanins, betalains and carotenoids. Plant J. 54, 733-749 (2008).

[20] Heuer, S., Vogt, T., Boehm, H. & Strack, D . Partial purification and characterization of UDP-glucose: betanidin 5-O- and 6-O-glucosyltransferases from cell suspension cultures of Dorotheanthus bellidiformis (Burm. F.) N.E.Br. Planta 199, 244-250 (1996).

[21] Hans, J., Brandt, W. & Vogt, T. Site-directed mutagenesis and protein 3D-homology modelling suggest a catalytic mechanism for UDP-glucose-dependent betanidin 5-O-glucosyltransferase from Dorotheanthus bellidiformis. Plant J. 39, 319-333 (2004).

[22] Sasaki, N. et al. Isolation and characterization of cDNAs encoding an enzyme with glucosyltransferase activity for cyclo-DOPA from four o’clocks and feather cockscombs. Plant Cell Physiol. 46, 666-670 (2005).

[23] Polturak, G. & Aharoni, A. “la vie en rose”: Biosynthesis, sources, and applications of betalain pigments. Mol. Plant 11, 7-22 (2018).

[24] Imamura, T. et al. Isolation of amaranthin synthetase from Chenopodium quinoa and construction of an amaranthin production system using suspension-cultured tobacco BY-2 cells. Plant Biotechnol. J. 17, 969-981 (2019).

[25] Hatlestad, G. J. et al. The beet Y locus encodes an anthocyanin MYB-like protein that activates the betalain red pigment pathway. Nat. Genet. 47, 92-96 (2015).

[26] Hua, Q. Transcriptomic analysis reveals key genes related to betalain biosynthesis in pulp coloration of Hylocereus polyrhizus. Front. Plant Sci. 6, 1179 (2016).

[27] Xie, F. et al. Transcriptomics-based identification and characterization of glucosyltransferases involved in betalain biosynthesis in Hylocereus megalanthus. Plant Physiol. Biochem. 152, 112-124 (2020).

[28] Cheng, M. et al. The WRKY transcription factor HpWRKY44 regulates CytP450-like1 expression in red pitaya fruit (Hylocereus polyrhizus). Hortic. Res. 4, 17039 (2017).

[29] Parra, G., Bradnam, K. & Korf, I. CEGMA: A pipeline to accurately annotate core genes in eukaryotic genornes. Bioinformatics 23, 1061-1067 (2007).

[30] Simao, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V. & Zdobnov, E. M. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210-3212 (2015).

[31] Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754-1760 (2009).

[32] Altschul, S. F. et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389-3402 (1997).

[33] Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. J. Mol. Biol. 215, 403-410 (1990).

[34] Nawrocki, E. P. & Eddy, S. R. Infernal 1.1: 100-Fold faster RNA homology searches. Bioinformatics 29, 2933-2935 (2013).

[35] Chen, F. et al. The sequenced angiosperm genomes and genome databases. Front. Plant Sci. 9, 418 (2018).

[36] Chen, F. et al. Genome sequences of horticultural plants: past, present, and future. Hortic. Res. 6, 112 (2019).

[37] Edgar, R. C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792-1797 (2004).

[38] Han, M. V., Thomas, G. W. C., Lugo-Martinez, J. & Hahn, M. W. Estimating gene gain and loss rates in the presence of error in genome assembly and annotation using CAFE 3. Mol. Biol. Evol. 30, 1987-1997 (2013).

[39] Xu, Z. & Wang, H. LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Res. 35, W265-W268 (2007).

[40] Ellinghaus, D., Kurtz, S. & Willhoeft, U. LTRharvest, an efficient and flexible software for de novo detection of LTR retrotransposons. BMC Bioinforma. 9, 18 (2008).

[41] Ossowski, S. et al. The rate and molecular spectrum of spontaneous mutations in Arabidopsis thaliana. Science. 327, 92-94 (2010).

[42] Steinbiss, S., Willhoeft, U., Gremme, G. & Kurtz, S. Fine-grained annotation and classification of de novo predicted LTR retrotransposons. Nucleic Acids Res. 37, 7002-7013 (2009).

[43] Vilella, A. J. et al. Ensembl compara genetrees: complete, duplication-aware phylogenetic trees in vertebrates. Genome Res. 19, 327-335 (2009).

[44] Macas, J., Neumann, P. & Navratilova, A . Repetitive DNA in the pea (Pisum sativum L.) genome: comprehensive characterization using 454 sequencing and comparison to soybean and Medicago truncatula. BMC Genomics 8, 427 (2007).

[45] Salse, J. Ancestors of modern plant crops. Curr. Opin. Plant Biol. 30, 134-142 (2016).

[46] Murat, F., Armero, A., Pont, C., Klopp, C. & Salse, J. Reconstructing the genome of the most recent common ancestor of flowering plants. Nat. Genet. 49, 490-496 (2017).

[47] Robinson, M. D., McCarthy, D. J. & Smyth, G. K. EdgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139-140 (2010).

[48] Langfelder, P. & Horvath, S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinforma. 9, 559 (2008).

[49] Shannon, P. et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498-2504 (2003).

[50] Dohm, J. C. et al. The genome of the recently domesticated crop plant sugar beet (Beta vulgaris). Nature 505, 546-549 (2014).

[51] Guo, Y. Gene family evolution in green plants with emphasis on the origination and evolution of Arabidopsis thaliana genes. Plant J. 73, 941-951 (2013).

[52] Guan, R. et al. Draft genome of the living fossil Ginkgo biloba. Gigascience 5, 49 (2016).

[53] De Bie, T., Cristianini, N., Demuth, J. P. & Hahn, M. W. CAFE: a computational tool for the study of gene family evolution. Bioinformatics 22, 1269-1271 (2006).

[54] Nystedt, B. et al. The Norway spruce genome sequence and conifer genome evolution. Nature 497, 579-584 (2013).

[55] Xu, C. et al. Draft genome of spinach and transcriptome diversity of 120 Spinacia accessions. Nat. Commun. 8, 15275 (2017).

[56] Yagi, M. et al. Sequence analysis of the genome of carnation (Dianthus caryophyllus L.). DNA Res. 21, 231-241 (2014).

[57] Zhou, Z. et al. Combined transcriptome and metabolome analysis of pitaya fruit unveiled the mechanisms underlying peel and pulp color formation. BMC Genomics 21, 734 (2020).

[58] Strack, D., Vogt, T. & Schliemann, W. Recent advances in betalain research. Phytochem 62, 247-269 (2003).

[59] Gandia-Herrero, F. & Garcia-Carmona, F. Biosynthesis of betalains: yellow and violet plant pigments. Trends Plant Sci. 18, 334-343 (2013).

[60] Timoneda, A. et al. Redirecting primary metabolism to boost production of tyrosine-derived specialised metabolites in planta. Sci. Rep. 8, 17256 (2018).

[61] Chung, H. et al. Characterisation of betalain biosynthesis in Parakeelya flowers identifies the key biosynthetic gene DOD as belonging to an expanded LigB gene family that is conserved in betalain-producing species. Front. Plant Sci. 6, 499 (2015).

[62] Sheehan, H. et al. Evolution of L-DOPA 4,5-dioxygenase activity allows for recurrent specialisation to betalain pigmentation in Caryophyllales. N. Phytol. 227, 914-929 (2020).

[63] Xu, W., Dubos, C. & Lepiniec, L. Transcriptional control of flavonoid biosynthesis by MYB-bHLH-WDR complexes. Trends Plant Sci. 20, 176-185 (2015).

[64] Ma, D. & Constabel, C. P. MYB repressors as regulators of phenylpropanoid metabolism in plants. Trends Plant Sci. 24, 275-289 (2019).

[65] Lloyd, A. et al. Advances in the MYB-bHLH-WD repeat (MBW) pigment regulatory model: addition of a WRKY factor and co-option of an anthocyanin MYB for betalain regulation. Plant Cell Physiol. 58, 1431-1441 (2017).

[66] Timoneda, A. et al. The evolution of betalain biosynthesis in Caryophyllales. N. Phytol. 224, 71-85 (2019).