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

Chromosome-scale genome assembly of Docynia delavayi provides new insights into the α-farnesene biosynthesis

Docynia delavayi is an economically significant fruit species with a high market potential due to the special aroma of its fruit. Here, a 653.34 Mb high-quality genome of D. delavayi was first reported, of which 93.8 % of the sequences (612.98 Mb) could be anchored to 17 chromosomes, containing 48,325 protein-coding genes. Ks analysis proved that two whole genome duplication (WGD) events occurred in D. delavayi, resulting in the expansion of genes associated with terpene biosynthesis, which promoted its fruit-specific aroma production. Combined multi-omics analysis, α-farnesene was detected as the most abundant aroma substance emitted by D. delavayi fruit during storage, meanwhile one α-farnesene synthase gene (AFS) and 15 transcription factors (TFs) were identified as the candidate genes potentially involved in α-farnesene biosynthesis. Further studies for the regulation network of α-farnesene biosynthesis revealed that DdebHLH, DdeERF1 and DdeMYB could activate the transcription of DdeAFS. To our knowledge, it is the first report that MYB TF plays a regulatory role in α-farnesene biosynthesis, which will greatly facilitate future breeding programs for D. delavayi.

Isoform-resolved genome annotation enables mapping of tissue-specific betalain regulation in amaranth

Betalains are coloring pigments produced in some families of the order Caryophyllales, where they replace anthocyanins as coloring pigments. While the betalain pathway itself is well studied, the tissue-specific regulation of the pathway remains mostly unknown. We enhance the high-quality Amaranthus hypochondriacus reference genome and produce a substantially more complete genome annotation, incorporating isoform details. We annotate betalain and anthocyanin pathway genes along with their regulators in amaranth and map the genetic control and tissue-specific regulation of the betalain pathway. Our improved genome annotation allowed us to identify causal mutations that lead to a knock-out of red betacyanins in natural accessions of amaranth. We reveal the tissue-specific regulation of flower color via a previously uncharacterized MYB transcription factor, AhMYB2. Downregulation of AhMYB2 in the flower leads to reduced expression of key betalain enzyme genes and loss of red flower color. Our improved amaranth reference genome represents the most complete genome of amaranth to date and is a valuable resource for betalain and amaranth research. High similarity of the flower betalain regulator AhMYB2 to anthocyanin regulators and a partially conserved interaction motif support the co-option of anthocyanin regulators for the betalain pathway as a possible reason for the mutual exclusiveness of the two pigments.

Genome-wide analysis of the PYL-PP2C-SnRK2s family in the ABA signaling pathway of pitaya reveals its expression profiles under canker disease stress

BACKGROUND

Abscisic acid (ABA) plays a crucial role in seed dormancy, germination, and growth, as well as in regulating plant responses to environmental stresses during plant growth and development. However, detailed information about the PYL-PP2C-SnRK2s family, a central component of the ABA signaling pathway, is not known in pitaya.

RESULTS

In this study, we identified 19 pyrabactin resistance-likes (PYLs), 70 type 2 C protein phosphatases (PP2Cs), and 14 SNF1-related protein kinase 2s (SnRK2s) from pitaya. In pitaya, tandem duplication was the primary mechanism for amplifying the PYL-PP2C-SnRK2s family. Co-linearity analysis revealed more homologous PYL-PP2C-SnRK2s gene pairs located in collinear blocks between pitaya and Beta vulgaris L. than that between pitaya and Arabidopsis. Transcriptome analysis showed that the PYL-PP2C-SnRK2s gene family plays a role in pitaya's response to infection by N. dimidiatum. By spraying ABA on pitaya and subsequently inoculating it with N. dimidiatum, we conducted qRT-PCR experiments to observe the response of the PYL-PP2C-SnRK2s gene family and disease resistance-related genes to ABA. These treatments significantly enhanced pitaya's resistance to pitaya canker. Further protein interaction network analysis helped us identify five key PYLs genes that were upregulated during the interaction between pitaya and N. dimidiatum, and their expression patterns were verified by qRT-PCR. Subcellular localization analysis revealed that the PYL (Hp1879) gene is primarily distributed in the nucleus.

CONCLUSION

This study enhances our understanding of the response of PYL-PP2C-SnRK2s to ABA and also offers a new perspective on pitaya disease resistance.

Unique regulatory network of dragon fruit simultaneously mitigates the effect of vanadium pollutant and environmental factors

Under changing climatic conditions, plants are simultaneously facing conflicting stresses in nature. Plants can sense different stresses, induce systematic ROS signals, and regulate transcriptomic, hormonal, and stomatal responses. We performed transcriptome analysis to reveal the integrative stress response regulatory mechanism underlying heavy metal stress alone or in combination with heat and drought conditions in pitaya (dragon fruit). A total of 70 genes were identified from 31,130 transcripts with conserved differential expression. Furthermore, weighted gene co-expression network analysis (WGCNA) identified trait-associated modules. By integrating information from three modules and protein-protein interaction (PPI) networks, we identified 10 interconnected genes associated with the multifaceted defense mechanism employed by pitaya against co-occurring stresses. To further confirm the reliability of the results, we performed a comparative analysis of 350 genes identified by three trait modules and 70 conserved genes exhibiting their dynamic expression under all treatments. Differential expression pattern of genes and comparative analysis, have proven instrumental in identifying ten putative structural genes. These ten genes were annotated as PLAT/LH2, CAT, MLP, HSP, PB1, PLA, NAC, HMA, and CER1 transcription factors involved in antioxidant activity, defense response, MAPK signaling, detoxification of metals and regulating the crosstalk between the complex pathways. Predictive analysis of putative candidate genes, potentially governing single, double, and multifactorial stress response, by several signaling systems and molecular patterns. These findings represent a valuable resource for pitaya breeding programs, offering the potential to develop resilient "super pitaya" plants.

Effect of BvAgNP on growth, development, and glyoxalase gene expression analysis in Mammillaria bombycina and Selenicereus undatus

BACKGROUND

Silver nanoparticles (AgNPs) have been used in plant tissue culture as growth stimulants, promoting bud initiation, germination, and rooting. In prior studies, AgNPs were synthesized and characterized by green synthesis using extracts from Beta vulgaris var. cicla (BvAgNP), and their functionality as seed disinfectant and antimicrobial was verified. In this study, we evaluated the effect of BvAgNP on the growth and development of Mammillaria bombycina and Selenicereus undatus in vitro, as well as the expression of glyoxalase genes.

METHODS

Explants from M. bombycina and S. undatus in vitro were treated with 25, 50, and 100 mg/L of BvAgNP. After 90 days, morphological characteristics were evaluated, and the expression of glyoxalase genes was analyzed by qPCR.

RESULTS

All treatments inhibited rooting for M. bombycina and no bud initiation was observed. S. undatus, showed a maximum response in rooting and bud generation at 25 mg/L of BvAgNP. Scanning electron microscopy (SEM) results exhibited a higher number of vacuoles in stem cells treated with BvAgNP compared to the control for both species. Expression of glyoxalase genes in M. bombycina increased in all treatments, whereas it decreased for S. undatus, however, increasing in roots.

CONCLUSIONS

This study presents the effects of BvAgNP on the growth and development of M. bombycina and S. undatus, with the aim of proposing treatments that promote in vitro rooting and bud initiation.

A chromosome-level genome reveals genome evolution and molecular basis of anthraquinone biosynthesis in Rheum palmatum

BACKGROUND

Rhubarb is one of common traditional Chinese medicine with a diverse array of therapeutic efficacies. Despite its widespread use, molecular research into rhubarb remains limited, constraining our comprehension of the geoherbalism.

RESULTS

We assembled the genome of Rheum palmatum L., one of the source plants of rhubarb, to elucidate its genome evolution and unpack the biosynthetic pathways of its bioactive compounds using a combination of PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C scaffolding approaches. Around 2.8 Gb genome was obtained after assembly with more than 99.9% sequences anchored to 11 pseudochromosomes (scaffold N50 = 259.19 Mb). Transposable elements (TE) with a continuous expansion of long terminal repeat retrotransposons (LTRs) is predominant in genome size, contributing to the genome expansion of R. palmatum. Totally 30,480 genes were predicted to be protein-coding genes with 473 significantly expanded gene families enriched in diverse pathways associated with high-altitude adaptation for this species. Two successive rounds of whole genome duplication event (WGD) shared by Fagopyrum tataricum and R. palmatum were confirmed. We also identified 54 genes involved in anthraquinone biosynthesis and other 97 genes entangled in flavonoid biosynthesis. Notably, RpALS emerged as a compelling candidate gene for the octaketide biosynthesis after the key residual screening.

CONCLUSION

Overall, our findings offer not only an enhanced understanding of this remarkable medicinal plant but also pave the way for future innovations in its genetic breeding, molecular design, and functional genomic studies.

Functional characterisation of Dof gene family and expression analysis under abiotic stresses and melatonin-mediated tolerance in pitaya (Selenicereus undatus)

DNA binding proteins with one finger (Dof ) transcription factors are essential for seed development and defence against various biotic and abiotic stresses in plants. Genomic analysis of Dof has not been determined yet in pitaya (Selenicereus undatus ). In this study, we have identified 26 Dof gene family members, renamed as HuDof-1 to HuDof-26 , and clustered them into seven subfamilies based on conserved motifs, domains, and phylogenetic analysis. The gene pairs of Dof family members were duplicated by segmental duplications that faced purifying selection, as indicated by the K a /K s ratio values. Promoter regions of HuDof genes contain many cis -acting elements related to phytohormones including abscisic acid, jasmonic acid, gibberellin, temperature, and light. We exposed pitaya plants to different environmental stresses and examined melatonin's influence on Dof gene expression levels. Signifcant expression of HuDof -2 and HuDof -6 were observed in different developmental stages of flower buds, flowers, pericarp, and pulp. Pitaya plants were subjected to abiotic stresses, and transcriptome analysis was carried out to identify the role of Dof gene family members. RNA-sequencing data and reverse transcription quantitative PCR-based expression analysis revealed three putative candidate genes (HuDof -1, HuDof -2, and HuDof -8), which might have diverse roles against the abiotic stresses. Our study provides a theoretical foundation for functional analysis through traditional and modern biotechnological tools for pitaya trait improvement.

The CAM lineages of planet Earth

BACKGROUND AND SCOPE

The growth of experimental studies of crassulacean acid metabolism (CAM) in diverse plant clades, coupled with recent advances in molecular systematics, presents an opportunity to re-assess the phylogenetic distribution and diversity of species capable of CAM. It has been more than two decades since the last comprehensive lists of CAM taxa were published, and an updated survey of the occurrence and distribution of CAM taxa is needed to facilitate and guide future CAM research. We aimed to survey the phylogenetic distribution of these taxa, their diverse morphology, physiology and ecology, and the likely number of evolutionary origins of CAM based on currently known lineages.

RESULTS AND CONCLUSIONS

We found direct evidence (in the form of experimental or field observations of gas exchange, day-night fluctuations in organic acids, carbon isotope ratios and enzymatic activity) for CAM in 370 genera of vascular plants, representing 38 families. Further assumptions about the frequency of CAM species in CAM clades and the distribution of CAM in the Cactaceae and Crassulaceae bring the currently estimated number of CAM-capable species to nearly 7 % of all vascular plants. The phylogenetic distribution of these taxa suggests a minimum of 66 independent origins of CAM in vascular plants, possibly with dozens more. To achieve further insight into CAM origins, there is a need for more extensive and systematic surveys of previously unstudied lineages, particularly in living material to identify low-level CAM activity, and for denser sampling to increase phylogenetic resolution in CAM-evolving clades. This should allow further progress in understanding the functional significance of this pathway by integration with studies on the evolution and genomics of CAM in its many forms.

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.

Complete mitogenome assembly of Selenicereus monacanthus revealed its molecular features, genome evolution, and phylogenetic implications

BACKGROUND

Mitochondria are the powerhouse of the cell and are critical for plant growth and development. Pitaya (Selenicereus or Hylocereus) is the most important economic crop in the family Cactaceae and is grown worldwide, however its mitogenome is unreported.

RESULTS

This study assembled the complete mitogenome of the red skin and flesh of pitaya (Selenicereus monacanthus). It is a full-length, 2,290,019 bp circular molecule encoding 59 unique genes that only occupy 2.17% of the entire length. In addition, 4,459 pairs of dispersed repeats (≥ 50 bp) were identified, accounting for 84.78% of the total length, and three repeats (394,588, 124,827, and 13,437 bp) mediating genomic recombination were identified by long read mapping and Sanger sequencing. RNA editing events were identified in all 32 protein-coding genes (PCGs), among which four sites (nad1-2, nad4L-2, atp9-copy3-223, and ccmFC-1309) were associated with the initiation or termination of PCGs. Seventy-eight homologous fragments of the chloroplast genome were identified in the mitogenome, the longest having 4,523 bp. In addition, evolutionary analyses suggest that S. monacanthus may have undergone multiple genomic reorganization events during evolution, with the loss of at least nine PCGs (rpl2, rpl10, rps2, rps3, rps10, rps11, rps14, rps19, and sdh3).

CONCLUSIONS

This study revealed the genetic basis of the S. monacanthus mitogenome, and provided a scientific basis for further research on phenotypic traits and germplasm resource development.

A chromosome-level genome assembly provides insights into Cornus wilsoniana evolution, oil biosynthesis, and floral bud development

Cornus wilsoniana W. is a woody oil plant with high oil content and strong hypolipidemic effects, making it a valuable species for medicinal, landscaping, and ecological purposes in China. To advance genetic research on this species, we employed PacBio together with Hi-C data to create a draft genome assembly for C. wilsoniana. Based on an 11-chromosome anchored chromosome-level assembly, the estimated genome size was determined to be 843.51 Mb. The N50 contig size and N50 scaffold size were calculated to be 4.49 and 78.00 Mb, respectively. Furthermore, 30 474 protein-coding genes were annotated. Comparative genomics analysis revealed that C. wilsoniana diverged from its closest species ~12.46 million years ago (Mya). Furthermore, the divergence between Cornaceae and Nyssaceae occurred >62.22 Mya. We also found evidence of whole-genome duplication events and whole-genome triplication γ, occurring at ~44.90 and 115.86 Mya. We further inferred the origins of chromosomes, which sheds light on the complex evolutionary history of the karyotype of C. wilsoniana. Through transcriptional and metabolic analysis, we identified two FAD2 homologous genes that may play a crucial role in controlling the oleic to linoleic acid ratio. We further investigated the correlation between metabolites and genes and identified 33 MADS-TF homologous genes that may affect flower morphology in C. wilsoniana. Overall, this study lays the groundwork for future research aimed at identifying the genetic basis of crucial traits in C. wilsoniana.

Pitaya Nutrition, Biology, and Biotechnology: A Review

Pitaya (Hylocereus spp.) is a member of the cactus family that is native to Central and South America but is now cultivated throughout the sub-tropical and tropical regions of the world. It is of great importance due to its nutritional, ornamental, coloring, medicinal, industrial, and high consumption values. In order to effectively utilize and develop the available genetic resources, it is necessary to appreciate and understand studies pertaining to the usage, origin, nutrition, diversity, evaluation, characterization, conservation, taxonomy, and systematics of the genus Hylocereus. Additionally, to gain a basic understanding of the biology of the plant, this review has also discussed how biotechnological tools, such as cell and tissue culture, micropropagation (i.e., somatic embryogenesis, organogenesis, somaclonal variation, mutagenesis, androgenesis, gynogenesis, and altered ploidy), virus-induced gene silencing, and molecular marker technology, have been used to enhance pitaya germplasm.

Identification of HuSWEET Family in Pitaya (Hylocereus undatus) and Key Roles of HuSWEET12a and HuSWEET13d in Sugar Accumulation

The sugar composition and content of fruit have a significant impact on their flavor and taste. In pitaya, or dragon fruit, sweetness is a crucial determinant of fruit taste and consumer preference. The sugars will eventually be exported transporters (SWEETs), a novel group of sugar transporters that have various physiological functions, including phloem loading, seed filling, nectar secretion, and fruit development. However, the role of SWEETs in sugar accumulation in pitaya fruit is not yet clear. Here, we identified 19 potential members (HuSWEET genes) of the SWEET family in pitaya and analyzed their conserved motifs, physiochemical characteristics, chromosomal distribution, gene structure, and phylogenetic relationship. Seven highly conserved α-helical transmembrane domains (7-TMs) were found, and the HuSWEET proteins can be divided into three clades based on the phylogenetic analysis. Interestingly, we found two HuSWEET genes, HuSWEET12a and HuSWEET13d, that showed strong preferential expressions in fruits and an upward trend during fruit maturation, suggesting they have key roles in sugar accumulation in pitaya. This can be further roughly demonstrated by the fact that transgenic tomato plants overexpressing HuSWEET12a/13d accumulated high levels of sugar in the mature fruit. Together, our result provides new insights into the regulation of sugar accumulation by SWEET family genes in pitaya fruit, which also set a crucial basis for the further functional study of the HuSWEETs.

A chromosome-scale and haplotype-resolved genome assembly of carnation (Dianthus caryophyllus) based on high-fidelity sequencing

Dianthus caryophyllus is an economic species often considered excellent cut flowers and is suitable for bouquets and gardens. Here, we assembled the haplotype-resolved genome of D. caryophyllus 'Aili' at the chromosome level for the first time. The total lengths of the two assembled haplotypes of carnation were 584.88 Mb for haplotype genome 1 (hap1) and 578.78 Mb for haplotype genome 2 (hap2), respectively. We predicted a total of 44,098 and 42,425 protein-coding genes, respectively. The remarkable structure variation was identified between two haplotypes. Moreover, we identified 403.80 Mb of transposable elements (TEs) in hap1, which accounted for 69.34% of the genome. In contrast, hap2 had 402.70 Mb of TEs, representing 69.61% of the genome. Long terminal repeats were the predominant transposable elements. Phylogenetic analysis showed that the species differentiation time between carnation and gypsophila was estimated to be ~54.43 MYA. The unique gene families of carnation genomes were identified in 'Aili' and previously published 'Francesco' and 'Scarlet Queen'. The assembled and annotated haplotype-resolved D. caryophyllus genome not only promises to facilitate molecular biology studies but also contributes to genome-level evolutionary studies.

The Pitaya Flower Tissue's Gene Differential Expression Analysis between Self-Incompatible and Self-Compatible Varieties for the Identification of Genes Involved in Self-Incompatibility Regulation

Self-incompatible pitaya varieties have low fruit-setting rates under natural conditions, leading to higher production costs and hindering industrial prosperity. Through transcriptome sequencing, we obtained the 36,900 longest transcripts (including 9167 new transcripts) from 60 samples of flowers. Samples were collected pre- and post-pollination (at 0 h, 0.5 h, 2 h, 4 h, and 12 h) from two varieties of pitaya (self-compatible Jindu No. 1 and self-incompatible Cu Sha). Using the RNA-Seq data and comparison of reference genomes, we annotated 28,817 genes in various databases, and 1740 genes were optimized in their structure for annotation. There were significant differences in the expression of differentially expressed genes (DEGs) in the pitaya stigmas under different pollination types, especially at the late post-pollination stage, where the expression of protease genes increasedal significantly under cross-pollination. We identified DEGs involved in the ribosomal, ubiquitination-mediated, and phyto-signaling pathways that may be involved in pitaya SI regulation. Based on the available transcriptome data and bioinformatics analysis, we tentatively identified HuS-RNase2 as a candidate gynogenetic S gene in the pitaya GSI system.

Genomic Identification of CCCH-Type Zinc Finger Protein Genes Reveals the Role of HuTZF3 in Tolerance of Heat and Salt Stress of Pitaya (Hylocereus polyrhizus)

Pitaya (Hylocereus polyrhizus) is cultivated in a broad ecological range, due to its tolerance to drought, heat, and poor soil. The zinc finger proteins regulate gene expression at the transcriptional and post-transcriptional levels, by interacting with DNA, RNA, and proteins, to play roles in plant growth and development, and stress response. Here, a total of 81 CCCH-type zinc finger protein genes were identified from the pitaya genome. Transcriptomic analysis showed that nine of them, including HuTZF3, responded to both salt and heat stress. RT-qPCR results showed that HuTZF3 is expressed in all tested organs of pitaya, with a high level in the roots and stems, and confirmed that expression of HuTZF3 is induced by salt and heat stress. Subcellular localization showed that HuTZF3 is targeted in the processing bodies (PBs) and stress granules (SGs). Heterologous expression of HuTZF3 could improve both salt and heat tolerance in Arabidopsis, reduce oxidative stress, and improve the activity of catalase and peroxidase. Therefore, HuTZF3 may be involved in post-transcriptional regulation via localizing to PBs and SGs, contributing to both salt and heat tolerance in pitaya.

Betalain biosynthesis in red pulp pitaya is regulated via HuMYB132: a R-R type MYB transcription factor

BACKGROUND

Multiple MYB transcription factors (TFs) are involved in the regulation of plant coloring. Betalain is a kind of natural plant pigment and its biosynthesis is regulated by a number of enzymes. Despite this, little is known about the molecular properties and roles of MYB TFs in pitaya betalain biosynthesis.

RESULTS

In the present study, we identified a 1R-MYB gene, HuMYB132, which is preferentially expressed in red-pulp pitaya at the mature stage. It was clustered with Arabidopsis R-R-type genes and had two DNA-binding domains and a histidine-rich region. The expression assays in N. benthamiana and yeast indicated that HuMYB132 is a nucleus-localized protein with transcriptional activation activity. Dual luciferase reporter assay and electrophoretic mobility shift assays (EMSA) demonstrated that HuMYB132 could promote the transcriptional activities of HuADH1, HuCYP76AD1-1, and HuDODA1 by binding to their promoters. Silencing HuMYB132 reduced betalain accumulation and the expression levels of betalain biosynthetic genes in pitaya pulps.

CONCLUSIONS

According to our findings, HuMYB132, a R-R type member of 1R-MYB TF subfamily, positively regulates pitaya betalain biosynthesis by regulating the expression of HuADH1, HuCYP76AD1-1, and HuDODA1. The present study provides a new theoretical reference for the management of pitaya betalain biosynthesis and also provides an essential basis for future regulation of betalain biosynthesis in Hylocereus.

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.

Genome-Wide Identification and Expression Pattern of the GRAS Gene Family in Pitaya (Selenicereus undatus L.)

The GRAS gene family is one of the most important families of transcriptional factors that have diverse functions in plant growth and developmental processes including axillary meristem patterning, signal-transduction, cell maintenance, phytohormone and light signaling. Despite their importance, the function of GRAS genes in pitaya fruit (Selenicereus undatus L.) remains unknown. Here, 45 members of the HuGRAS gene family were identified in the pitaya genome, which was distributed on 11 chromosomes. All 45 members of HuGRAS were grouped into nine subfamilies using phylogenetic analysis with six other species: maize, rice, soybeans, tomatoes, Medicago truncatula and Arabidopsis. Among the 45 genes, 12 genes were selected from RNA-Seq data due to their higher expression in different plant tissues of pitaya. In order to verify the RNA-Seq data, these 12 HuGRAS genes were subjected for qRT-PCR validation. Nine HuGRAS genes exhibited higher relative expression in different tissues of the plant. These nine genes which were categorized into six subfamilies inlcuding DELLA (HuGRAS-1), SCL-3 (HuGRAS-7), PAT1 (HuGRAS-34, HuGRAS-35, HuGRAS-41), HAM (HuGRAS-37), SCR (HuGRAS-12) and LISCL (HuGRAS-18, HuGRAS-25) might regulate growth and development in the pitaya plant. The results of the present study provide valuable information to improve tropical pitaya through a molecular and conventional breeding program.

Chromosome-scale Amaranthus tricolor genome provides insights into the evolution of the genus Amaranthus and the mechanism of betalain biosynthesis

Amaranthus tricolor is a vegetable and ornamental amaranth, with high lysine, dietary fibre and squalene content. The red cultivar of A. tricolor possesses a high concentration of betalains, which has been used as natural food colorants. Here, we constructed the genome of A. tricolor, the first reference genome for the subgenus Albersia, combining PacBio HiFi, Nanopore ultra-long and Hi-C data. The contig N50 size was 906 kb, and 99.58% of contig sequence was anchored to the 17 chromosomes, totalling 520 Mb. We annotated 27,813 protein-coding genes with an average 1.3 kb coding sequence and 5.3 exons. We inferred that A. tricolor underwent a whole-genome duplication (WGD) and that the WGD shared by amaranths occurred in the last common ancestor of subfamily Amaranthoideae. Moreover, we comprehensively identified candidate genes in betalain biosynthesis pathway. Among them, DODAα1 and CYP76ADα1, located in one topologically associated domain (TAD) of an active (A) compartment on chromosome 16, were more highly expressed in red leaves than in green leaves, and DODAα1 might be the rate-limiting enzyme gene in betalains biosynthesis. This study presents new genome resources and enriches our understanding of amaranth evolution, betalains production, facilitating molecular breeding improvements and the understanding of C4 plants evolution.

Building a foundation for gene family analysis in Rosaceae genomes with a novel workflow: A case study in Pyrus architecture genes

The rapid development of sequencing technologies has led to a deeper understanding of plant genomes. However, direct experimental evidence connecting genes to important agronomic traits is still lacking in most non-model plants. For instance, the genetic mechanisms underlying plant architecture are poorly understood in pome fruit trees, creating a major hurdle in developing new cultivars with desirable architecture, such as dwarfing rootstocks in European pear (Pyrus communis). An efficient way to identify genetic factors for important traits in non-model organisms can be to transfer knowledge across genomes. However, major obstacles exist, including complex evolutionary histories and variable quality and content of publicly available plant genomes. As researchers aim to link genes to traits of interest, these challenges can impede the transfer of experimental evidence across plant species, namely in the curation of high-quality, high-confidence gene models in an evolutionary context. Here we present a workflow using a collection of bioinformatic tools for the curation of deeply conserved gene families of interest across plant genomes. To study gene families involved in tree architecture in European pear and other rosaceous species, we used our workflow, plus a draft genome assembly and high-quality annotation of a second P. communis cultivar, 'd'Anjou.' Our comparative gene family approach revealed significant issues with the most recent 'Bartlett' genome - primarily thousands of missing genes due to methodological bias. After correcting assembly errors on a global scale in the 'Bartlett' genome, we used our workflow for targeted improvement of our genes of interest in both P. communis genomes, thus laying the groundwork for future functional studies in pear tree architecture. Further, our global gene family classification of 15 genomes across 6 genera provides a valuable and previously unavailable resource for the Rosaceae research community. With it, orthologs and other gene family members can be easily identified across any of the classified genomes. Importantly, our workflow can be easily adopted for any other plant genomes and gene families of interest.

Diploid chromosome-level reference genome and population genomic analyses provide insights into Gypenoside biosynthesis and demographic evolution of Gynostemma pentaphyllum (Cucurbitaceae)

Gynostemma pentaphyllum (Thunb.) Makino is a perennial creeping herbaceous plant in the family Cucurbitaceae, which has great medicinal value and commercial potential, but urgent conservation efforts are needed due to the gradual decreases and fragmented distribution of its wild populations. Here, we report the high-quality diploid chromosome-level genome of G. pentaphyllum obtained using a combination of next-generation sequencing short reads, Nanopore long reads, and Hi-C sequencing technologies. The genome is anchored to 11 pseudo-chromosomes with a total size of 608.95 Mb and 26 588 predicted genes. Comparative genomic analyses indicate that G. pentaphyllum is estimated to have diverged from Momordica charantia 60.7 million years ago, with no recent whole-genome duplication event. Genomic population analyses based on genotyping-by-sequencing and ecological niche analyses indicated low genetic diversity but a strong population structure within the species, which could classify 32 G. pentaphyllum populations into three geographical groups shaped jointly by geographic and climate factors. Furthermore, comparative transcriptome analyses showed that the genes encoding enzyme involved in gypenoside biosynthesis had higher expression levels in the leaves and tendrils. Overall, the findings obtained in this study provide an effective molecular basis for further studies of demographic genetics, ecological adaption, and systematic evolution in Cucurbitaceae species, as well as contributing to molecular breeding, and the biosynthesis and biotransformation of gypenoside.

A high-quality Buxus austro-yunnanensis (Buxales) genome provides new insights into karyotype evolution in early eudicots

BACKGROUND

Eudicots are the most diverse group of flowering plants that compromise five well-defined lineages: core eudicots, Ranunculales, Proteales, Trochodendrales, and Buxales. However, the phylogenetic relationships between these five lineages and their chromosomal evolutions remain unclear, and a lack of high-quality genome analyses for Buxales has hindered many efforts to address this knowledge gap.

RESULTS

Here, we present a high-quality chromosome-level genome of Buxus austro-yunnanensis (Buxales). Our phylogenomic analyses revealed that Buxales and Trochodendrales are genetically similar and classified as sisters. Additionally, both are sisters to the core eudicots, while Ranunculales was found to be the first lineage to diverge from these groups. Incomplete lineage sorting and hybridization were identified as the main contributors to phylogenetic discordance (34.33%) between the lineages. In fact, B. austro-yunnanensis underwent only one whole-genome duplication event, and collinear gene phylogeny analyses suggested that separate independent polyploidizations occurred in the five eudicot lineages. Using representative genomes from these five lineages, we reconstructed the ancestral eudicot karyotype (AEK) and generated a nearly gapless karyotype projection for each eudicot species. Within core eudicots, we recovered one common chromosome fusion event in asterids and malvids, respectively. Further, we also found that the previously reported fused AEKs in Aquilegia (Ranunculales) and Vitis (core eudicots) have different fusion positions, which indicates that these two species have different karyotype evolution histories.

CONCLUSIONS

Based on our phylogenomic and karyotype evolution analyses, we revealed the likely relationships and evolutionary histories of early eudicots. Ultimately, our study expands genomic resources for early-diverging eudicots.

Genome-Wide Identification of WRKY Gene Family in Pitaya Reveals the Involvement of HmoWRKY42 in Betalain Biosynthesis

The WRKY gene family is a plant-specific transcription factor (TF) that regulates many physiological processes and (a) biotic stress responses. Despite this, little is known about the molecular properties and roles of WRKY TFs in pitaya betalain biosynthesis. Here we report the identification of 70 WRKY in Hylocereus undatus, their gene structure, locations on each chromosome, systematic phylogenetic analysis, conserved motif analysis, and synteny of HuWRKY genes. HmoWRKY42 is a Group IIb WRKY protein and contains a coiled-coil motif, a WRKY domain and a C₂H₂ zinc-finger motif (CX₅CX₂₃HXH). Results from yeast one-hybrid and transient dual-luciferase assays showed that HmoWRKY42 was a transcriptional repressor and could repress HmocDOPA5GT1 expression by binding to its promoter. Yeast two-hybrid assays showed that HmoWRKY42 could interact with itself to form homodimers. Knocking out the coiled-coil motif of HmoWRKY42 prevented its self-interaction and prevented it from binding to the HmocDOPA5GT1 promoter. Knocking out the WRKY domain and C₂H₂ zinc-finger motif sequence of HmoWRKY42 also prevented it from binding to the HmocDOPA5GT1 promoter. The coiled-coil motif, the WRKY domain and the C₂H₂ zinc finger motif are key motifs for the binding of HmoWRKY42 to the HmocDOPA5GT1 promoter. HmoWRKY42 is localized in the nucleus and possesses trans-activation ability responsible for pitaya betalain biosynthesis by repressing the transcription of HmocDOPA5GT1. As far as we know, no reports are available on the role of HmoWRKY42 in pitaya betalain biosynthesis. The results provide an important foundation for future analyses of the regulation and functions of the HuWRKY gene family.

Prospects and potentials of underutilized leafy Amaranths as vegetable use for health-promotion

Climate change causes environmental variation worldwide, which is one of the most serious threats to global food security. In addition, more than 2 billion people in the world are reported to suffer from serious malnutrition, referred to as 'hidden hunger.' Dependence on only a few crops could lead to the loss of genetic diversity and high fragility of crop breeding in systems adapting to global scale climate change. The exploitation of underutilized species and genetic resources, referred to as orphan crops, could be a useful approach for resolving the issue of adaptability to environmental alteration, biodiversity preservation, and improvement of nutrient quality and quantity to ensure food security. Moreover, the use of these alternative crops will help to increase the human health benefits and the income of farmers in developing countries. In this review, we highlight the potential of orphan crops, especially amaranths, for use as vegetables and health-promoting nutritional components. This review highlights promising diversified sources of amaranth germplasms, their tolerance to abiotic stresses, and their nutritional, phytochemical, and antioxidant values for vegetable purposes. Betalains (betacyanins and betaxanthins), unique antioxidant components in amaranth vegetables, are also highlighted regarding their chemodiversity across amaranth germplasms and their stability and degradation. In addition, we discuss the physiological functions, antioxidant, antilipidemic, anticancer, and antimicrobial activities, as well as the biosynthesis pathway, molecular, biochemical, genetics, and genomic mechanisms of betalains in detail.

The association of HpDof1.7 and HpDof5.4 with soluble sugar accumulation in pitaya fruit by transcriptionally activating sugar metabolic genes

Soluble sugar is one of the important factors affecting fruit flavor and quality. Here, we report the identification of two Dof (DNA-binding with one finger) transcription factors termed HpDof1.7 and HpDof5.4, and their roles in influencing sugar accumulation in pitayas. HpDof1.7 and HpDof5.4 shared a similar expression pattern with sugar metabolism-related genes HpSuSy1 and HpINV2, and sugar transporter genes HpTMT2 and HpSWEET14 during pitayas maturation, and their expression pattern was also consistent with the accumulation of glucose and fructose, which were the predominant sugars in pitayas. HpDof1.7 and HpDof5.4 were both typical nucleus-localized proteins with trans-activation ability. Gel mobility shift assay revealed that HpDof1.7 and HpDof5.4 were bound to promoters of HpSuSy1, HpINV2, HpTMT2 and HpSWEET14. Finally, transient expression assays in tobacco leaves showed that HpDof1.7 and HpDof5.4 increased the activities of HpSuSy1, HpINV2, HpTMT2 and HpSWEET14 promoters, thus facilitating sugar accumulation by transcriptionally enhancing sugar metabolic pathway genes. Our findings provide a new perspective on the regulatory mechanisms of Dof transcription factors in sugar accumulation and pitaya fruit quality formation.

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.

Wheat-Thinopyrum Substitution Lines Imprint Compensation Both From Recipients and Donors

Even frequently used in wheat breeding, we still have an insufficient understanding of the biology of the products via distant hybridization. In this study, a transcriptomic analysis was performed for six Triticum aestivum-Thinopyrum elongatum substitution lines in comparison with the host plants. All the six disomic substitution lines showed much stronger "transcriptomic-shock" occurred on alien genomes with 57.43-69.22% genes changed expression level but less on the recipient genome (2.19-8.97%). Genome-wide suppression of alien genes along chromosomes was observed with a high proportion of downregulated genes (39.69-48.21%). Oppositely, the wheat recipient showed genome-wide compensation with more upregulated genes, occurring on all chromosomes but not limited to the homeologous groups. Moreover, strong co-upregulation of the orthologs between wheat and Thinopyrum sub-genomes was enriched in photosynthesis with predicted chloroplastic localization, which indicates that the compensation happened not only on wheat host genomes but also on alien genomes.

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.

Metabolic Profiling of Sugars and Organic Acids, and Expression Analyses of Metabolism-Associated Genes in Two Yellow-Peel Pitaya Species

Sugar and organic acids are important factors determining pitaya fruit quality. However, changes in sugars and acids, and expressions of metabolism-associated genes during fruit maturation of yellow-peel pitayas are not well-documented. In this study, metabolic and expression analyses in pulps of different fruit developmental stages of ‘Wucihuanglong’ (‘WCHL’, Hylocereus undatus) and ‘Youcihuanglong’ pitaya (‘YCHL’, Hylocereus megalanthus) were used to explore the sugar and organic acid metabolic process. Total phenols and flavonoids were mainly accumulated at S1 in pitaya pulps. Ascorbic acid contents of ‘WCHL’ pitaya were higher than that of ‘YCHL’ pitaya during fruit maturation. Starch was mainly accumulated at early fruit development stages while soluble sugars were rich in late stages. Sucrose, fructose, and glucose were the main sugar components of ‘YCHL’ pitaya while glucose was dominant in ‘WCHL’ pitaya. Malic and citric acids were the main organic acids in ‘WCHL’ and ‘YCHL’ pitayas, respectively. Based on the transcriptome analyses, 118 genes involved in pitaya sugar and organic acid metabolism were obtained. Results from the correlation analyses between the expression profiling of candidate genes and the contents of sugar and organic acid showed that 51 genes had a significant correlation relationship and probably perform key role in pitaya sugar and organic acid metabolism processes. The finding of the present study provides new information for quality regulation of pitayas.

Evolutionary Genetics of Cacti: Research Biases, Advances and Prospects

Here, we present a review of the studies of evolutionary genetics (phylogenetics, population genetics, and phylogeography) using genetic data as well as genome scale assemblies in Cactaceae (Caryophyllales, Angiosperms), a major lineage of succulent plants with astonishing diversity on the American continent. To this end, we performed a literature survey (1992-2021) to obtain detailed information regarding key aspects of studies investigating cactus evolution. Specifically, we summarize the advances in the following aspects: molecular markers, species delimitation, phylogenetics, hybridization, biogeography, and genome assemblies. In brief, we observed substantial growth in the studies conducted with molecular markers in the past two decades. However, we found biases in taxonomic/geographic sampling and the use of traditional markers and statistical approaches. We discuss some methodological and social challenges for engaging the cactus community in genomic research. We also stressed the importance of integrative approaches, coalescent methods, and international collaboration to advance the understanding of cactus evolution.

HuNAC20 and HuNAC25, Two Novel NAC Genes from Pitaya, Confer Cold Tolerance in Transgenic Arabidopsis

NAC transcription factors are one of the largest families of transcriptional regulators in plants, and members of the gene family play vital roles in regulating plant growth and development processes including biotic/abiotic stress responses. However, little information is available about the NAC family in pitaya. In this study, we conducted a genome-wide analysis and a total of 64 NACs (named HuNAC1-HuNAC64) were identified in pitaya (Hylocereus). These genes were grouped into fifteen subgroups with diversities in gene proportions, exon–intron structures, and conserved motifs. Genome mapping analysis revealed that HuNAC genes were unevenly scattered on all eleven chromosomes. Synteny analysis indicated that the segmental duplication events played key roles in the expansion of the pitaya NAC gene family. Expression levels of these HuNAC genes were analyzed under cold treatments using qRT-PCR. Four HuNAC genes, i.e., HuNAC7, HuNAC20, HuNAC25, and HuNAC30, were highly induced by cold stress. HuNAC7, HuNAC20, HuNAC25, and HuNAC30 were localized exclusively in the nucleus. HuNAC20, HuNAC25, and HuNAC30 were transcriptional activators while HuNAC7 was a transcriptional repressor. Overexpression of HuNAC20 and HuNAC25 in Arabidopsis thaliana significantly enhanced tolerance to cold stress through decreasing ion leakage, malondialdehyde (MDA), and H₂O₂ and O₂⁻ accumulation, accompanied by upregulating the expression of cold-responsive genes (AtRD29A, AtCOR15A, AtCOR47, and AtKIN1). This study presents comprehensive information on the understanding of the NAC gene family and provides candidate genes to breed new pitaya cultivars with tolerance to cold conditions through genetic transformation.

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.

Genome-Wide Characterization of R2R3-MYB Transcription Factors in Pitaya Reveals a R2R3-MYB Repressor HuMYB1 Involved in Fruit Ripening through Regulation of Betalain Biosynthesis by Repressing Betalain Biosynthesis-Related Genes

The MYB (myeloblastosis) superfamily constitutes one of the most abundant transcription factors (TFs) regulating various biological processes in plants. However, the molecular characteristics and functions of MYB TFs in pitaya remain unclear. To date, no genome-wide characterization analysis of this gene family has been conducted in the Cactaceae species. In this study, 105 R2R3-MYB members were identified from the genome data of Hylocereus undatus and their conserved motifs, physiological and biochemical characteristics, chromosome locations, synteny relationship, gene structure and phylogeny were further analyzed. Expression analyses suggested that three up-regulated HuMYBs and twenty-two down-regulated HuMYBs were probably involved in fruit ripening of pitaya. Phylogenetic analyses of R2R3-MYB repressors showed that seven HuMYBs (HuMYB1, HuMYB21, HuMYB48, HuMYB49, HuMYB72, HuMYB78 and HuMYB101) were in clades containing R2R3-MYB repressors. HuMYB1 and HuMYB21 were significantly down-regulated with the betalain accumulation during fruit ripening of ‘Guanhuahong’ pitaya (H. monacanthus). However, only HuMYB1 had R2 and R3 repeats with C1, C2, C3 and C4 motifs. HuMYB1 was localized exclusively to the nucleus and exhibited transcriptional inhibition capacities. Dual luciferase reporter assay demonstrated that HuMYB1 inhibited the expression of betalain-related genes: HuADH1, HuCYP76AD1-1 and HuDODA1. These results suggested that HuMYB1 is a potential repressor of betalain biosynthesis during pitaya fruit ripening. Our results provide the first genome-wide analyses of the R2R3-MYB subfamily involved in pitaya betalain biosynthesis and will facilitate functional analysis of this gene family in the future.

Genome-Wide Identification of Aquaporin Gene Family in Pitaya Reveals an HuNIP6;1 Involved in Flowering Process

Aquaporins (AQPs) are essential membrane proteins involved in seed maturation and germination, stomata movement, photosynthesis, and regulation of plant flowering processes. Pitaya flowers are open at night and wither at daybreak, which shows an obvious circadian rhythm. In this study, a comprehensive genome-wide analysis of AQPs in Hylocereus undantus was conducted to screen key genes associated with flowering processes. A total of 33 HuAQP genes were identified from the H. undantus genome. The 33 HuAQPs were grouped into four subfamilies: 10 PIPs, 13 TIPs, 8 NIPs, and 2 SIPs, which were distributed on 9 out of 11 pitaya chromosomes (Chr) (except for Chr7 and Chr10). Results from expression profiles showed that HuNIP6;1 may be involved in pitaya's floral opening. HuNIP6;1 was localized exclusively in the cell membrane. Overexpression of HuNIP6;1 in Arabidopsis thaliana significantly promoted early flowering through regulating negative flowering regulators of MJM30, COL9, and PRR5, suggesting that HuNIP6;1 plays key roles in regulating flowering time. The present study provides the first genome-wide analysis of the AQP gene family in pitaya and valuable information for utilization of HuAQPs.