The chromosome-level genome of dragon fruit reveals whole-genome duplication and chromosomal co-localization of betacyanin biosynthetic genes

Altitudinal variation of dragon fruit metabolite profiles as revealed by UPLC-MS/MS-based widely targeted metabolomics analysis

BACKGROUND

Geographical factors affect the nutritional, therapeutic and commercial values of fruits. Dragon fruit (Hylocereus spp) is a popular fruit in Asia and a potential functional food with diverse pharmacological attributes. Although it is produced in various localities, the information related to the altitudinal variation of dragon fruit nutrients and active compounds is scarce. Hence, this study aimed to investigate the variations in metabolite profiles of H. polyrhizus (variety Jindu1) fruit pulps from three different altitudes of China, including Wangmo (WM, 650 m), Luodian (LD, 420 m), and Zhenning (ZN, 356 m). Jindu1 is the main cultivated pitaya variety in Guizhou province, China.

RESULTS

The LC-MS (liquid chromatography-mass spectroscopy)-based widely targeted metabolic profiling identified 645 metabolites, of which flavonoids (22.64%), lipids (13.80%), phenolic acids (12.40%), amino acids and derivatives (10.39%), alkaloids (8.84%), and organic acids (8.37%) were dominant. Multivariate analyses unveiled that the metabolite profiles of the fruit differed regarding the altitude. Fruits from WM (highest altitude) were prime in quality, with higher levels of flavonoids, alkaloids, nucleotides and derivatives, amino acids and derivatives, and vitamins. Fruits from LD and ZN had the highest relative content of phenolic acids and terpenoids, respectively. We identified 69 significantly differentially accumulated metabolites across the pulps of the fruits from the three locations. KEGG analysis revealed that flavone and flavonol biosynthesis and isoflavonoid biosynthesis were the most differentially regulated. It was noteworthy that most active flavonoid compounds exhibited an increasing accumulation pattern along with the increase in altitude. Vitexin and isovitexin were the major differentially accumulated flavonoids. Furthermore, we identified two potential metabolic biomarkers (vitexin and kaempferol 3-O-[2-O-β-D-galactose-6-O-a-L-rhamnose]-β-D-glucoside) to discriminate between dragon fruits from different geographical origins.

CONCLUSION

Our findings provide insights into metabolic changes in dragon fruits grown at different altitudes. Furthermore, they show that growing pitaya at high altitudes can produce fruit with higher levels of bioactive compounds, particularly flavonoids.

Advances in the Understanding of Postharvest Physiological Changes and the Storage and Preservation of Pitaya

Highly prized for its unique taste and appearance, pitaya is a tasty, low-calorie fruit. It has a high-water content, a high metabolism, and a high susceptibility to pathogens, resulting in an irreversible process of tissue degeneration or quality degradation and eventual loss of commercial value, leading to economic loss. High quality fruits are a key guarantee for the healthy development of economic advantages. However, the understanding of postharvest conservation technology and the regulation of maturation, and senescence of pitaya are lacking. To better understand the means of postharvest storage of pitaya, extend the shelf life of pitaya fruit and prospect the postharvest storage technology, this paper analyzes and compares the postharvest quality changes of pitaya fruit, preservation technology, and senescence regulation mechanisms. This study provides research directions for the development of postharvest storage and preservation technology.

The American Cherimoya Genome Reveals Insights into the Intra-Specific Divergence, the Evolution of Magnoliales, and a Putative Gene Cluster for Acetogenin Biosynthesis

Annona cherimola (cherimoya) is a species renowned for its delectable fruit and medicinal properties. In this study, we developed a chromosome-level genome assembly for the cherimoya 'Booth' cultivar from the United States. The genome assembly has a size of 794 Mb with a N50 = 97.59 Mb. The seven longest scaffolds account for 87.6% of the total genome length, which corresponds to the seven pseudo-chromosomes. A total of 45,272 protein-coding genes (≥30 aa) were predicted with 92.9% gene content completeness. No recent whole genome duplications were identified by an intra-genome collinearity analysis. Phylogenetic analysis supports that eudicots and magnoliids are more closely related to each other than to monocots. Moreover, the Magnoliales was found to be more closely related to the Laurales than the Piperales. Genome comparison revealed that the 'Booth' cultivar has 200 Mb less repeats than the Spanish cultivar 'Fino de Jete', despite their highly similar (>99%) genome sequence identity and collinearity. These two cultivars were diverged during the early Pleistocene (1.93 Mya), which suggests a different origin and domestication of the cherimoya. Terpene/terpenoid metabolism functions were found to be enriched in Magnoliales, while TNL (Toll/Interleukin-1-NBS-LRR) disease resistance gene has been lost in Magnoliales during evolution. We have also identified a gene cluster that is potentially responsible for the biosynthesis of acetogenins, a class of natural products found exclusively in Annonaceae. The cherimoya genome provides an invaluable resource for supporting characterization, conservation, and utilization of Annona genetic resources.

Gene duplications facilitate C4-CAM compatibility in common purslane

Common purslane (Portulaca oleracea) integrates both C4 and crassulacean acid metabolism (CAM) photosynthesis pathways and is a promising model plant to explore C4-CAM plasticity. Here, we report a high-quality chromosome-level genome of nicotinamide adenine dinucleotide (NAD)-malic enzyme (ME) subtype common purslane that provides evidence for 2 rounds of whole-genome duplication (WGD) with an ancient WGD (P-β) in the common ancestor to Portulacaceae and Cactaceae around 66.30 million years ago (Mya) and another (Po-α) specific to common purslane lineage around 7.74 Mya. A larger number of gene copies encoding key enzymes/transporters involved in C4 and CAM pathways were detected in common purslane than in related species. Phylogeny, conserved functional site, and collinearity analyses revealed that the Po-α WGD produced the phosphoenolpyruvate carboxylase-encoded gene copies used for photosynthesis in common purslane, while the P-β WGD event produced 2 ancestral genes of functionally differentiated (C4- and CAM-specific) beta carbonic anhydrases involved in the C4 + CAM pathways. Additionally, cis-element enrichment analysis in the promoters showed that CAM-specific genes have recruited both evening and midnight circadian elements as well as the Abscisic acid (ABA)-independent regulatory module mediated by ethylene-response factor cis-elements. Overall, this study provides insights into the origin and evolutionary process of C4 and CAM pathways in common purslane, as well as potential targets for engineering crops by integrating C4 or CAM metabolism.

Elucidation of the mescaline biosynthetic pathway in peyote (Lophophora williamsii)

Peyote (Lophophora williamsii) is an entheogenic and medicinal cactus native to the Chihuahuan desert. The psychoactive and hallucinogenic properties of peyote are principally attributed to the phenethylamine derivative mescaline. Despite the isolation of mescaline from peyote over 120 years ago, the biosynthetic pathway in the plant has remained undiscovered. Here, we use a transcriptomics and homology-guided gene discovery strategy to elucidate a near-complete biosynthetic pathway from l-tyrosine to mescaline. We identified a cytochrome P450 that catalyzes the 3-hydroxylation of l-tyrosine to l-DOPA, a tyrosine/DOPA decarboxylase yielding dopamine, and four substrate-specific and regiospecific substituted phenethylamine O-methyltransferases. Biochemical assays with recombinant enzymes or functional analyses performed by feeding putative precursors to engineered yeast (Saccharomyces cerevisiae) strains expressing candidate peyote biosynthetic genes were used to determine substrate specificity, which served as the basis for pathway elucidation. Additionally, an N-methyltransferase displaying broad substrate specificity and leading to the production of N-methylated phenethylamine derivatives was identified, which could also function as an early step in the biosynthesis of tetrahydroisoquinoline alkaloids in peyote.

Plastome variations reveal the distinct evolutionary scenarios of plastomes in the subfamily Cereoideae (Cactaceae)

BACKGROUND

The cactus family (Cactaceae) has been reported to have evolved a minimal photosynthetic plastome size, with the loss of inverted-repeat (IR) regions and NDH gene suites. However, there are very limited genomic data on the family, especially Cereoideae, the largest subfamily of cacti.

RESULTS

In the present study, we assembled and annotated 35 plastomes, 33 of which were representatives of Cereoideae, alongside 2 previously published plastomes. We analyzed the organelle genomes of 35 genera in the subfamily. These plastomes have variations rarely observed in those of other angiosperms, including size differences (with ~ 30 kb between the shortest and longest), dramatic dynamic changes in IR boundaries, frequent plastome inversions, and rearrangements. These results suggested that cacti have the most complex plastome evolution among angiosperms.

CONCLUSION

These results provide unique insight into the dynamic evolutionary history of Cereoideae plastomes and refine current knowledge of the relationships within the subfamily.

Stem transcriptome screen for selection in wild and cultivated pitahaya (Selenicereus undatus): an epiphytic cactus with edible fruit

Dragon fruit, pitahaya or pitaya are common names for the species in the Hylocereus group of Selenicereus that produce edible fruit. These Neotropical epiphytic cacti are considered promising underutilized crops and are currently cultivated around the world. The most important species, S. undatus, has been managed in the Maya domain for centuries and is the focus of this article. Transcriptome profiles from stems of wild and cultivated plants of this species were compared. We hypothesized that differences in transcriptomic signatures could be associated with genes related to drought stress. De novo transcriptome assembly and the analysis of differentially expressed genes (DEGs) allowed us to identify a total of 9,203 DEGs in the Hunucmá cultivar relative of wild Mozomboa plants. Of these, 4,883 represent up-regulated genes and 4,320, down-regulated genes. Additionally, 6,568 DEGs were identified from a comparison between the Umán cultivar and wild plants, revealing 3,286 up-regulated and 3,282 down-regulated genes. Approximately half of the DEGs are shared by the two cultivated plants. Differences between the two cultivars that were collected in the same region could be the result of differences in management. Metabolism was the most representative functional category in both cultivars. The up-regulated genes of both cultivars formed a network related to the hormone-mediated signaling pathway that includes cellular responses to auxin stimulus and to hormone stimulus. These cellular reactions have been documented in several cultivated plants in which drought-tolerant cultivars modify auxin transport and ethylene signaling, resulting in a better redistribution of assimilates.

The genome of the recretohalophyte Limonium bicolor provides insights into salt gland development and salinity adaptation during terrestrial evolution

Halophytes have evolved specialized strategies to cope with high salinity. The extreme halophyte sea lavender (Limonium bicolor) lacks trichomes but possesses salt glands on its epidermis that can excrete harmful ions, such as sodium, to avoid salt damage. Here, we report a high-quality, 2.92-Gb, chromosome-scale L. bicolor genome assembly based on a combination of Illumina short reads, single-molecule, real-time long reads, chromosome conformation capture (Hi-C) data, and Bionano genome maps, greatly enriching the genomic information on recretohalophytes with multicellular salt glands. Although the L. bicolor genome contains genes that show similarity to trichome fate genes from Arabidopsis thaliana, it lacks homologs of the decision fate genes GLABRA3, ENHANCER OF GLABRA3, GLABRA2, TRANSPARENT TESTA GLABRA2, and SIAMESE, providing a molecular explanation for the absence of trichomes in this species. We identified key genes (LbHLH and LbTTG1) controlling salt gland development among classical trichome homologous genes and confirmed their roles by showing that their mutations markedly disrupted salt gland initiation, salt secretion, and salt tolerance, thus offering genetic support for the long-standing hypothesis that salt glands and trichomes may share a common origin. In addition, a whole-genome duplication event occurred in the L. bicolor genome after its divergence from Tartary buckwheat and may have contributed to its adaptation to high salinity. The L. bicolor genome resource and genetic evidence reported in this study provide profound insights into plant salt tolerance mechanisms that may facilitate the engineering of salt-tolerant crops.

Identification of Allergens in White- and Red-Fleshed Pitaya (Selenicereus undatus and Selenicereus costaricensis) Seeds Using Bottom-Up Proteomics Coupled with Immunoinformatics

White-fleshed pitaya (Selenicereus undatus) and red-fleshed pitaya (Selenicereus costaricensis) are becoming increasingly popular because of their nutritional and medicinal benefits. However, in addition to their beneficial properties, allergy to pitaya fruits has occurred in daily life. In this study, we investigated the protein profile of pitaya fruit seeds and focused on the most reactive proteins against immunoglobulin E (IgE) in sera from allergic patients by immunoblotting. A protein band of approximately 20 kDa displayed a clear reaction with the serum IgE. The protein bands of interest were excised, in-gel digested, and analyzed using liquid chromatography–tandem mass spectrometry (LC–MS/MS), followed by data searching against a restricted database (Caryophyllales in UniProtKB) for protein identification. Immunoinformatic tools were used to predict protein allergenicity. The potential allergens included cupin_1 and heat shock protein 70 (HSP70) in white-fleshed pitaya seeds, and cupin_1, heat shock protein 70, and heat shock protein sti1-like in red-fleshed pitaya seeds are potential allergens. The expression of potential allergens was further verified at the transcriptional level in the species of S. undatus and S. costaricensis.

Pitaya Genome and Multiomics Database (PGMD): A Comprehensive and Integrative Resource of Selenicereus undatus

Pitaya (Selenicereus) is a kind of novel fruit with a delicious taste and superior horticulture ornamental value. The potential economic impact of the pitaya lies in its diverse uses not only as agricultural produce and processed foods but also in industrial and medicinal products. It is also an excellent plant material for basic and applied biological research. A comprehensive database of pitaya would facilitate studies of pitaya and the other Cactaceae plant species. Here, we constructed pitaya genome and multiomics database, which is a collection of the most updated and high-quality pitaya genomic assemblies. The database contains various information such as genomic variation, gene expression, miRNA profiles, metabolite and proteomic data from various tissues and fruit developmental stages of different pitaya cultivars. In PGMD, we also uploaded videos on the flowering process and planting tutorials for practical usage of pitaya. Overall, these valuable data provided in the PGMD will significantly facilitate future studies on population genetics, molecular breeding and function research of pitaya.

A target Capture Probe Set Useful for Deep- and Shallow-Level Phylogenetic Studies in Cactaceae

The molecular phylogenies of Cactaceae have enabled us to better understand their systematics, biogeography, and diversification ages. However, most of the phylogenetic relationships within Cactaceae major groups remain unclear, largely due to the lack of an appropriate set of molecular markers to resolve its contentious relationships. Here, we explored the genome and transcriptome assemblies available for Cactaceae and identified putative orthologous regions shared among lineages of the subfamily Cactoideae. Then we developed a probe set, named Cactaceae591, targeting both coding and noncoding nuclear regions for representatives from the subfamilies Pereskioideae, Opuntioideae, and Cactoideae. We also sampled inter- and intraspecific variation to evaluate the potential of this panel to be used in phylogeographic studies. We retrieved on average of 547 orthologous regions per sample. Targeting noncoding nuclear regions showed to be crucial to resolving inter- and intraspecific relationships. Cactaceae591 covers 13 orthologous genes shared with the Angiosperms353 kit and two plastid regions largely used in Cactaceae studies, enabling the phylogenies generated by our panel to be integrated with angiosperm and Cactaceae phylogenies, using these sequences. We highlighted the importance of using coalescent-based species tree approaches on the Cactaceae591 dataset to infer accurate phylogenetic trees in the presence of extensive incomplete lineage sorting in this family.

Breeding an underutilized fruit crop: a long-term program for Hylocereus

This review describes three decades of introduction, agro-technology development, breeding and selection of Hylocereus species, known as pitaya or dragon fruit, as an example of a holistic program aimed to develop the horticultural potential of a perennial underutilized fruit crop. Interspecific homoploid and interploid crosses and embryo rescue procedures produced improved hybrids, some of which have been released to farmers. Molecular tools and morphological and phenological comparisons between the parental species and the resulting hybrids provided valuable information on dominant/recessive traits and on genetic relationships that could be exploited for further hybridizations. In addition, Hylocereus were crossed with species of the closely related genus Selenicereus, producing valuable intergeneric hybrids. In situ chromosome doubling resulted in the production of autopolyploid lines, from which an understanding of the effect of increased ploidy on fruit traits and metabolomic profiles was obtained. Gamete-derived lines were produced, adding to our biobank homozygote lines that were subsequently used for further hybridization. Spontaneous chromosome doubling occurred in haploid gamete-derived Hylocereus monacanthus lines and in interspecific interploid Hylocereus megalanthus × H. undatus hybrids obtained from an embryo rescue procedure, resulting in plants with double the expected ploidy. Challenging technical problems were addressed by the development of protocols for DNA isolation, flow cytometry, in situ chromosome doubling, androgenesis, gynogenesis and embryo rescue following interspecific and interploidy crosses. Current research leading to the development of genomics and molecular tools, including a draft genome of H. undatus, is also presented. Perspectives for further development of Hylocereus species and hybrids are discussed.

Origin and diversification of the saguaro cactus (Carnegiea gigantea): a within-species phylogenomic analysis

Reconstructing accurate historical relationships within a species poses numerous challenges, not least in many plant groups in which gene flow is high enough to extend well beyond species boundaries. Nonetheless, the extent of tree-like history within a species is an empirical question on which it is now possible to bring large amounts of genome sequence to bear. We assess phylogenetic structure across the geographic range of the saguaro cactus, an emblematic member of Cactaceae, a clade known for extensive hybridization and porous species boundaries. Using 200 Gb of whole genome resequencing data from 20 individuals sampled from 10 localities, we assembled two data sets comprising 150,000 biallelic single nucleotide polymorphisms (SNPs) from protein coding sequences. From these we inferred within-species trees and evaluated their significance and robustness using five qualitatively different inference methods. Despite the low sequence diversity, large census population sizes, and presence of wide-ranging pollen and seed dispersal agents, phylogenetic trees were well resolved and highly consistent across both data sets and all methods. We inferred that the most likely root, based on marginal likelihood comparisons, is to the east and south of the region of highest genetic diversity, which lies along the coast of the Gulf of California in Sonora, Mexico. Together with striking decreases in marginal likelihood found to the north, this supports hypotheses that saguaro's current range reflects post-glacial expansion from the refugia in the south of its range. We conclude with observations about practical and theoretical issues raised by phylogenomic data sets within species, in which SNP-based methods must be used rather than gene tree methods that are widely used when sequence divergence is higher. These include computational scalability, inference of gene flow, and proper assessment of statistical support in the presence of linkage effects.

Plant genome sequence assembly in the era of long reads: Progress, challenges and future directions

Third-generation long-read sequencing is transforming plant genomics. Oxford Nanopore Technologies and Pacific Biosciences are offering competing long-read sequencing technologies and enable plant scientists to investigate even large and complex plant genomes. Sequencing projects can be conducted by single research groups and sequences of smaller plant genomes can be completed within days. This also resulted in an increased investigation of genomes from multiple species in large scale to address fundamental questions associated with the origin and evolution of land plants. Increased accessibility of sequencing devices and user-friendly software allows more researchers to get involved in genomics. Current challenges are accurately resolving diploid or polyploid genome sequences and better accounting for the intra-specific diversity by switching from the use of single reference genome sequences to a pangenome graph.

Polyphenol utilization proteins in human gut microbiome

Dietary polyphenols can significantly benefit human health, but their bioavailability is metabolically controlled by human gut microbiota. To facilitate the study of polyphenol metabolism for human gut health, we have manually curated experimentally characterized polyphenol utilization proteins (PUPs) from published literature. This resulted in 60 experimentally characterized PUPs (named seeds) with various metadata, such as species and substrate. Further database search found 107,851 homologs of the seeds from UniProt and UHGP (Unified Human Gastrointestinal Protein) databases. All PUP seeds and homologs were classified into protein classes, families and subfamilies based on Enzyme Commission (EC) numbers, Pfam (protein family) domains and sequence similarity networks. By locating PUP homologs in the genomes of UHGP, we have identified 1,074 physically linked PUP gene clusters (PGCs), which are potentially involved in polyphenol metabolism in the human gut. The gut microbiome of Africans was consistently ranked the top in terms of the abundance and prevalence of PUP homologs and PGCs among all geographical continents. This reflects the fact that dietary polyphenols are more commonly consumed by African population than other populations such as Europeans and North Americans. A case study of the Hadza hunter-gatherer microbiome verified the feasibility of using dbPUP to profile metagenomic data for biologically meaningful discovery, suggesting an association between diet and PUP abundance. A Pfam domain enrichment analysis of PGCs identified a number of putatively novel PUP families. Lastly, a user-friendly web interface (https://bcb.unl.edu/dbpup/) provides all the data online to facilitate the research of polyphenol metabolism for improved human health. Importance Long-term consumption of polyphenol-rich foods have been shown to lower the risk of various human diseases such as cardiovascular diseases, cancers, and metabolic diseases. Raw polyphenols are often enzymatically processed by gut microbiome, which encode various polyphenol utilization proteins (PUPs) to produce metabolites with much higher bioaccessibility to gastrointestinal cells. This study delivered dbPUP as an online database for experimentally characterized PUPs and their homologs in human gut microbiome. This work also performed a systematic classification of PUPs into enzyme classes, families, and subfamilies. The signature Pfam domains were identified for PUP families, enabling conserved domain-based PUP annotation. This standardized sequence similarity-based PUP classification system offered a guideline for the future inclusion of new experimentally characterized PUPs and the creation of new PUP families. An in-depth data analysis was further conducted on PUP homologs and physically linked PUP gene clusters (PGCs) in gut microbiomes of different human populations.

Genomic Approaches for Improvement of Tropical Fruits: Fruit Quality, Shelf Life and Nutrient Content

The breeding of tropical fruit trees for improving fruit traits is complicated, due to the long juvenile phase, generation cycle, parthenocarpy, polyploidy, polyembryony, heterozygosity and biotic and abiotic factors, as well as a lack of good genomic resources. Many molecular techniques have recently evolved to assist and hasten conventional breeding efforts. Molecular markers linked to fruit development and fruit quality traits such as fruit shape, size, texture, aroma, peel and pulp colour were identified in tropical fruit crops, facilitating Marker-assisted breeding (MAB). An increase in the availability of genome sequences of tropical fruits further aided in the discovery of SNP variants/Indels, QTLs and genes that can ascertain the genetic determinants of fruit characters. Through multi-omics approaches such as genomics, transcriptomics, metabolomics and proteomics, the identification and quantification of transcripts, including non-coding RNAs, involved in sugar metabolism, fruit development and ripening, shelf life, and the biotic and abiotic stress that impacts fruit quality were made possible. Utilizing genomic assisted breeding methods such as genome wide association (GWAS), genomic selection (GS) and genetic modifications using CRISPR/Cas9 and transgenics has paved the way to studying gene function and developing cultivars with desirable fruit traits by overcoming long breeding cycles. Such comprehensive multi-omics approaches related to fruit characters in tropical fruits and their applications in breeding strategies and crop improvement are reviewed, discussed and presented here.

A Genome-Wide Identification Study Reveals That HmoCYP76AD1, HmoDODAα1 and HmocDOPA5GT Involved in Betalain Biosynthesis in Hylocereus

Betalains are water-soluble nitrogen-containing pigments with multiple bioactivities. Pitayas are the only at large-scale commercially grown fruit containing abundant betalains for consumers. Currently, the key genes involved in betalain biosynthesis remain to be fully elucidated. Moreover, genome-wide analyses of these genes in betalain biosynthesis are not available in betalain-producing plant species. In this study, totally 53 genes related to betalain biosynthesis were identified from the genome data of Hylocereus undatus. Four candidate genes i.e., one cytochrome P-450 R gene (HmoCYP76AD1), two L-DOPA 4,5-dioxygenase genes (HmoDODAα1 and HmoDODAα2), and one cyclo-DOPA 5-O glucosyltransferase gene (HmocDOPA5GT) were initially screened according to bioinformatics and qRT-PCR analyses. Silencing HmoCYP76AD1, HmoDODAα1, HmoDODAα2 or HmocDOPA5GT resulted in loss of red pigment. HmoDODAα1 displayed a high level of L-DOPA 4,5-dioxygenase activity to produce betalamic acid and formed yellow betaxanthin. Co-expression of HmoCYP76AD1, HmoDODAα1 and HmocDOPA5GT in Nicotiana benthamiana and yeast resulted in high abundance of betalain pigments with a red color. These results suggested that HmoCYP76AD1, HmoDODAα1, and HmocDOPA5GT play key roles in betalain biosynthesis in Hylocereus. The results of the present study provide novel genes for molecular breeding programs of pitaya.

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

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.