Search

Sugar composition not only affects fruit flavor but is also an important determinant of fruit taste and consumer preference. In this study, changes in the sugar content and sugar-metabolizing enzymes were investigated from different sections of various fruit development phases of ‘Hongshuijing’ pitaya (Hylocereus monacanthus). Genes related to sugar metabolism were also screened by transcriptome analyses. The results indicated that glucose was the major sugar in mature pitaya fruit, and was mainly regulated by vacuolar acid invertase (VAI) and sucrose synthase (SS) (degradative direction). Sugar accumulation varied in pulp between different sections of the pitaya fruit. VAI, neutral invertase (NI) and SS (degradative direction) are crucial enzymes for sugar accumulation in pitaya. The expression of 17 genes related to sucrose metabolism and obtained from seven databases [NCBI non-redundant protein database (Nr), NCBI non-redundant nucleotide sequence database (Nt), EuKaryotic Orthologous Groups (KOG), The Protein Families (Pfam), Kyoto Encyclopedia of Genes and Genomes (KEGG), Swiss-prot, and Gene Ontology (GO)] were analyzed in different pitaya pulp sections. HpVAI1 had the highest relative expression level on the 29th day after pollination (DAP). Positive correlations were found between HpVAI1 expression and VAI activity; HpNI4 and NI activity; HpSS2, HpSS5, and SS activity (synthetic direction), indicating that HpVAI1, HpNI4, and HpSS2 and HpSS5 were involved in the regulation of VAI, NI, and SS (synthetic direction), respectively. HpVAI1 and HpNI4 regulated sucrose degradation and the accumulation of glucose and fructose, while HpSS2 and HpSS5 regulated sucrose synthesis. These results suggest that HpVAI1 plays a key role in sugar metabolism during fruit development of ‘Hongshuijing’ pitaya. The results of this study provide new information about sugar metabolism in pitaya fruit that could help improve fruit quality and the breeding of new cultivars.

Sucrose synthase (SUS, EC 2.4.1.13) is widely considered as a key enzyme involved in plant sucrose metabolism, and the gene family encoding different SUS isozymes has been identified and characterized in several plant species. However, to date scant information about the SUS genes is available in Litchi chinensis Sonn. Here, we identified five SUS genes in litchi. These LcSUSs shared high levels of similarity in both nucleotide and amino acid sequences. Their gene structure, phylogenetic relationships, and expression profiles were characterized. Gene structure analysis indicated that the LcSUSs have similar exon-intron structures. Phylogenetic analysis revealed that the five members could be classified into three groups (LcSUS1 and LcSUS2 in SUS II, LcSUS4 and LcSUS5 in SUS III, and LcSUS3 in SUS I), demonstrating evolutionary conservation in the SUS family across litchi and other plant species. The expression levels of LcSUSs were investigated via real-time PCR in various tissues and different developmental stages of aril. For tissues and organs, LcSUSs exhibited distinct but partially redundant expression profiles in litchi, being predominantly expressed in young leaves (sink). During aril development, the expression pattern of LcSUS1 was consistent with the trend of sugar accumulation, indicating it may play important roles in determination of sink strength in aril. Moreover, transcript levels of LcSUS2, LcSUS4, and LcSUS5 varied between cultivars with different hexose/sucrose ratios, which may regulate the sugar composition in aril. Our results provide insights into physiological functions of SUS genes in litchi, especially roles in regulating sugar accumulation in aril.

Litchi chinensis Sonn is widely cultivated in subtropical regions and has an important economic value. A high-density genetic map is a valuable tool for mapping quantitative trait loci (QTL) and marker-assisted breeding programs. In this study, a single nucleotide polymorphism (SNP)-based high-density linkage map was constructed by a genotyping-by-sequencing (GBS) protocol using an F1 population of 178 progenies between two commercial litchi cultivars, ‘Ziniangxi’ (dwarf) and ‘Feizixiao’ (vigorous). The genetic map consisted of 3027 SNP markers with a total length of 1711.97 cM in 15 linkage groups (LGs) and an average marker distance of 0.57 cM. Based on this high-density linkage map and three years of phenotyping, a total of 37 QTLs were detected for eight dwarf-related traits, including length of new branch (LNB), diameter of new branch (DNB), length of common petiole (LCP), diameter of common petiole (DCP), length of internode (LI), length of single leaf (LSL), width of single leaf (WSL), and plant height (PH). These QTLs could explain 8.0 to 14.7% (mean=9.7%) of the phenotypic variation. Among them, several QTL clusters were observed, particularly on LG04 and LG11, which might show enrichment for genes regulating the dwarf-related traits in litchi. There were 126 candidate genes identified within the QTL regions, 55 of which are differentially expressed genes by RNA-seq analysis between ‘Ziniangxi’ and ‘Feizixiao’. These DEGs were found to participate in the regulation of cell development, material transportation, signal transduction, and plant morphogenesis, so they might play important roles in regulating plant dwarf-related traits. The high-density genetic map and QTLs identification related to dwarf traits can provide a valuable genetic resource and a basis for marker-assisted selection and genomic studies of litchi.

The

Pitaya (

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.

Background Litchi chinensis Sonn. is an economically important fruit tree in tropical and subtropical regions. However, litchi functional genomics is severely hindered due to its recalcitrance to regeneration and stable transformation. Agrobacterium rhizogenes-mediated hairy root transgenic system provide an alternative to study functional genomics in woody plants. However, the hairy root transgenic system has not been established in litchi. Results In this study, we report a rapid and highly efficient A. rhizogenes-mediated co-transformation system in L. chinensis using Green Fluorescent Protein (GFP) gene as a marker. Both leaf discs and stem segments of L. chinensis cv. ‘Fenhongguiwei’ seedlings were able to induce transgenic hairy roots. The optimal procedure involved the use of stem segments as explants, infection by A. rhizogenes strain MSU440 at optical density (OD600) of 0.7 for 10 min and co-cultivation for 3 days, with a co-transformation efficiency of 9.33%. Furthermore, the hairy root transgenic system was successfully used to validate the function of the key anthocyanin regulatory gene LcMYB1 in litchi. Over-expression of LcMYB1 produced red hairy roots, which accumulated higher contents of anthocyanins, proanthocyanins, and flavonols. Additionally, the genes involving in the flavonoid pathway were strongly activated in the red hairy roots. Conclusion We first established a rapid and efficient transformation system for the study of gene function in hairy roots of litchi using A. rhizogenes strain MSU440 by optimizing parameters. This hairy root transgenic system was effective for gene function analysis in litchi using the key anthocyanin regulator gene LcMYB1 as an example.

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.

Background SWEETs (Sugar Will Eventually be Exported transporters) function as sugar efflux transporters that perform diverse physiological functions, including phloem loading, nectar secretion, seed filling, and pathogen nutrition. The SWEET gene family has been identified and characterized in a number of plant species, but little is known about in Litchi chinensis, which is an important evergreen fruit crop. Results In this study, 16 LcSWEET genes were identified and nominated according to its homologous genes in Arabidopsis and grapevine. Multiple sequence alignment showed that the 7 alpha-helical transmembrane domains (7-TMs) were basically conserved in LcSWEETs. The LcSWEETs were divided into four clades (Clade I to Clade IV) by phylogenetic tree analysis. A total of 8 predicted motifs were detected in the litchi LcSWEET genes. The 16 LcSWEET genes were unevenly distributed in 9 chromosomes and there was one pairs of segmental duplicated events by synteny analysis. The expression patterns of the 16 LcSWEET genes showed higher expression levels in reproductive organs. The temporal and spatial expression patterns of LcSWEET2a and LcSWEET3b indicated they play central roles during early seed development. Conclusions The litchi genome contained 16 SWEET genes, and most of the genes were expressed in different tissues. Gene expression suggested that LcSWEETs played important roles in the growth and development of litchi fruits. Genes that regulate early seed development were preliminarily identified. This work provides a comprehensive understanding of the SWEET gene family in litchi, laying a strong foundation for further functional studies of LcSWEET genes and improvement of litchi fruits.

Background A suitable reference gene is an important prerequisite for guarantying accurate and reliable results in quantitative real-time PCR (qRT-PCR) analyses. However, there is no absolute universality in reference genes among different species. It’s hard to find an ideal reference gene to fit for different tissues and growth periods. Pitaya (Hylocereus) is commercially produced as a new fruit crop at a large scale in tropical and subtropical regions. To date, there is no report on the identification of the most reliable reference genes for qRT-PCR normalization in pitaya. Results In this study, six candidate reference genes i.e. Actin(1), GAPDH, UBC(1), UBC(2) EF1-α(1) and histone(1) were selected from thirty-nine typical candidate reference genes to determine the most stable reference genes for qRT-PCR normalization in different tissues, temperature stresses and fruit developmental stages of pitaya. Among the six candidate reference genes, Actin(1) and EF1-α(1) were the most stable gene according to calculations of three statistical methods (GeNorm, NormFinder and BestKeeper) while UBC(1) and UBC(2) showed the lowest expression stability. The six candidate reference genes were further validated by comparing expression profiles of key genes related to betalain biosynthesis at flesh coloration stages of Guanhuahong (Hylocereus monacanthus) and Guanhuabai (H. undatus) pitayas. Actin(1) was recommended the best reference gene for accurate normalization of qRT-PCR data. Conclusions In this study, the stability of the selected reference genes for normalizing the qRT-PCR data were identified from pitaya. Actin(1) was the most stably expressed genes in different tissues and fruit developmental stages in pitaya. The present work provides the first data of reference gene identification for pitaya and will facilitate further studies in molecular biology and gene function on Hylocereus and other closely related species. Electronic supplementary material The online version of this article (10.1186/s13007-019-0455-3) contains supplementary material, which is available to authorized users.

Anthocyanin accumulation is one of the remarkable physiological changes during fruit ripening. In plants, anthocyanin synthesis is regulated by MYB activators, but the MYB repressors has been recognized recently. Here, we isolated a repressor of anthocyanin synthesis, LcMYBx, from Litchi chinensis Sonn. LcMYBx encoded a typical R3-MYB protein and contained a conserved [D/E]Lx2[R/K]x3Lx6Lx3R motif for interacting with bHLH proteins. Overexpression of LcMYBx in tobacco suppressed anthocyanin accumulation resulting in faded petals from pale-pink to almost white. Gene expression analysis showed the strong down-regulation of endogenous anthocyanin structural and regulatory genes by LcMYBx overexpression. Yeast two-hybrid and bimolecular fluorescence complementation assays indicated that LcMYBx could interact with the transcription factors LcbHLH1 and LcbHLH3. Transient promoter activation assays showed that LcMYBx could inhibit the activation capacity of LcMYB1-LcbHLH3 complex for LcDFR gene. These results suggest that LcMYBx competed with LcMYB1 to LcbHLHs, thus preventing the activation of LcDFR by LcMYB1-LcbHLHs complex and negatively controlling anthocyanin biosynthesis.

Betalains are water-soluble nitrogen-containing pigments with multiple bioactivities. Pitaya is the only large-scale commercially grown fruit containing abundant betalains for consumers. However, the upstream regulators in betalain biosynthesis are still not clear. In this study, HmoWRKY40, a novel WRKY transcription factor, was obtained from the transcriptome data of pitaya (Hylocereus monacanthus). HmoWRKY40 is a member of the Group IIa WRKY family, containing a conserved WRKY motif, and it is located in the nucleus. The betalain contents and expression levels of HmoWRKY40 increased rapidly during the coloration of pitaya and reached their maximums on the 23rd day after artificial pollination (DAAP). Yeast one-hybrid and transient expression assays showed that HmoWRKY40 could bind and activate the promoter of HmoCYP76AD1. Silencing the HmoWRKY40 gene resulted in a significant reduction of betacyanin contents. These results indicate that HmoWRKY40 transcriptionally activates HmoCYP76AD, which is involved in the regulation of pitaya betalain biosynthesis. The results of the present study provide new regulatory networks related to betalain biosynthesis in pitaya.

S-ribonuclease (S-RNase)-based self-incompatibility (SI) mechanisms have been extensively studied in Solanaceae, Rosaceae and Plantaginaceae. S-RNase-based SI is controlled by two closely related genes, S-RNase and S-locus F-box (SLF), located at a polymorphic S-locus. In the SI system, the SCF-type (SKP1-CUL1-F-box-RBX1) complex functions as an E3 ubiquitin ligase complex for ubiquitination of non-self S-RNase. Pummelo (Citrus grandis) and several mandarin cultivars are suggested to utilize an S-RNase-based SI system. However, the molecular mechanism of the non-S-factors involved in the SI reaction is not straightforward in Citrus. To investigate the SCF-type E3 complex responsible for the SI reaction in mandarin, SLF, SKP1-like and CUL1 candidates potentially involved in the SI reaction of ‘Wuzishatangju’ (Citrus reticulata Blanco) were identified based on the genome-wide identification and expression analyses. Sixteen pollen-specific F-box genes (CrFBX1-CrFBX16), one pollen-specific SKP1-like gene (CrSKP1-e) and two CUL1 genes (CrCUL1A and CrCUL1B) were identified and cloned from ‘Wuzishatangju’. Yeast two-hybrid (Y2H) and in vitro binding assays showed that five CrFBX proteins could bind to CrSKP1-e, which is an ortholog of SSK1 (SLF-interacting-SKP1-like), a non-S-factor responsible for the SI reaction. Luciferase complementation imaging (LCI) and in vitro binding assays also showed that CrSKP1-e interacts with the N-terminal region of both CrCUL1A and CrCUL1B. These results indicate that CrSKP1-e may serve as a functional member of the SCF-type E3 ubiquitin ligase complex in ‘Wuzishatangju’.

Background MicroRNAs (miRNAs) and their regulatory functions in anthocyanin, carotenoid, and chlorophyll accumulation have been extensively characterized in many plant species. However, the miRNA regulatory mechanism in betalain biosynthesis remains mostly unknown. Results In this study, 126 conserved miRNAs and 41 novel miRNAs were first isolated from Hylocereus monacanthus, among which 95 conserved miRNAs belonged to 53 miRNA families. Thirty-four candidate miRNAs related to betalain biosynthesis were differentially expressed. The expression patterns of those differential expressed miRNAs were analyzed in various pitaya tissues by RT-qPCR. A significantly negative correlation was detected between the expression levels of half those miRNAs and corresponding target genes. Target genes of miRNAs i.e. Hmo-miR157b-HmSPL6-like, Hmo-miR160a-Hpcyt P450-like3, Hmo-miR6020-HmCYP71A8-like, Hmo-novel-2-HmCYP83B1-like, Hmo-novel-15-HmTPST-like, Hmo-miR828a-HmTT2-like, Hmo-miR858-HmMYB12-like, Hmo-miR858-HmMYBC1-like and Hmo-miR858-HmMYB2-like were verified by 5′RACE and transient expression system in tobacco. Conclusions Hmo-miR157b, Hmo-miR160a, Hmo-miR6020 Hmo-novel-2, Hmo-novel-15, Hmo-miR828a and Hmo-miR858 play important roles in pitaya fruit coloration and betalain accumulation. Our findings provide new insights into the roles of miRNAs and their target genes of regulatory functions involved in betalain biosynthesis of pitaya.

Longan (Dimocarpus longan) is a typical southern subtropical fruit tree species that is sensitive to cold stress. C-repeat binding factors (CBFs), as transcription factors, are crucial components involved in the molecular regulation of the plant response to cold stress. However, the role of CBF homologs in the cold response regulation of longan remains largely unknown. Here, three novel CBF genes, DlCBF1, DlCBF2, and DlCBF3, were cloned from longan. DlCBF1 and DlCBF2 contain an AP2 domain and PKKPAGR and DSAWR CBF signature motifs, while DlCBF3 has mutations within these conserved signature motifs. DlCBF1/2/3 were mainly localized in the nucleus and specifically bound to CRT/DRE cis-elements, resulting in strong transcriptional activation. DlCBF1/2 exhibited tissue expression specificity, and their expression was induced by low temperature, while DlCBF3 had no tissue specificity and barely responded to low temperature. DlCBF1, DlCBF2, and DlCBF3 overexpression in Arabidopsis-enhanced cold tolerance by increasing proline accumulation and reducing reactive oxygen species (ROS) content, accompanied by upregulated expression of cold-responsive genes (AtRD29A, AtCOR15A, AtCOR47, and AtKIN1) in the CBF cold stress response signaling pathway. In conclusion, the biological functions of DlCBF1/2/3 were somewhat conserved, but slow expression of DlCBF1/2 and low expression of DlCBF3 may partly cause the cold sensitivity of longan. Collectively, these results indicated that differences exist in the expression and function of CBF orthologs in the cold-sensitive plant species longan, and these findings may help to improve the understanding of the cold response regulation mechanism and provide important theoretical support for cold-tolerant breeding of longan.

Additional file 9: Table S2. The information of miRNAs in pitaya.

Additional file 17: Table S10. cDNA sequences of 17 target genes.

Additional file 8: Table S1. Statistics of miRNA sequences in pitaya.

Additional file 5: Figure S5. The sequences of 17 target genes from the two pitaya cultivars. ‘-W’, ‘Guanhuabai’ pitaya. ‘-R’, ‘Guanhuahong’ pitaya. Red lines indicate sequence of target genes primers for real-time PCR. Green lines indicate outer specific primers for nested PCR. Blue lines indicate inner specific primers for nested PCR.

Additional file 6: Figure S6. Transient expressions of miRNAs and their target genes (all). A1-A49, bright field; B1-B49, eGFP; C1-C49, merge; A1-C1, 35 s::eGFP; A2-C2, WT; A3-C3, 35 s::Pre-miR164; A4-C4, 35 s::Pre-Hmo-miR6020; A5-C5, 35 s::HmCYP71A8-like::eGFP; A6-C6, 35 s::positive control HmCYP71A8-like::eGFP; A7-C7, 35 s::negative control HmCYP71A8-like::eGFP; A8-C8, 35 s::HmCYP71A8-like::eGFP+ 35 s::Pre-miR164; A9-C9, 35 s::HmCYP71A8-like::eGFP+ 35 s::Pre-Hmo-miR6020; A10-C10, 35 s::positive control HmCYP71A8-like::eGFP+ 35 s::Pre-Hmo-miR6020; A11-C11, 35 s::negative control HmCYP71A8-like::eGFP+ 35 s::Pre-Hmo-miR6020; A12-C12, 35 s::Pre-Hmo-miR858; A13-C13, 35 s::HmMYB12-like::eGFP; A14-C14, 35 s::positive control HmMYB12-like::eGFP; A15-C15, 35 s::negative control HmMYB12-like::eGFP; A16-C16, 35 s::HmMYB12-like::eGFP+ 35 s::Pre-miR164; A17-C17, 35 s::HmMYB12-like::eGFP+ 35 s::Pre-Hmo-miR858; A18-C18, 35 s::positive control HmMYB12-like::eGFP+ 35 s::Pre-Hmo-miR858; A19-C19, 35 s::negative control HmMYB12-like::eGFP+ 35 s::Pre-Hmo-miR858; A20-C20, 35 s::HmMYBC1-like::eGFP; A21-C21, 35 s::HmMYBC1-like::eGFP+ 35 s::Pre-miR164; A22-C22, 35 s::HmMYBC1-like::eGFP+ 35 s::Pre-Hmo-miR858; A23-C23, 35 s::HmMYB2-like::eGFP; A24-C24, 35 s::HmMYB2-like::eGFP+ 35 s::Pre-miR164; A25-C25, 35 s::HmMYB2-like::eGFP+ 35 s::Pre-Hmo-miR858; A26-C26, 35 s::Pre-Hmo-miR160a; A27-C27, 35 s::Hpcyt P450-like2::eGFP; A28-C28, 35 s::positive control Hpcyt P450-like2::eGFP; A29-C29, 35 s::negative control Hpcyt P450-like2::eGFP; A30-C30, 35 s::Hpcyt P450-like2::eGFP+ 35 s::Pre-miR164; A31-C31, 35 s::Hpcyt P450-like2::eGFP+ 35 s::Pre-Hmo-miR160a; A32-C32, 35 s::positive control Hpcyt P450-like2::eGFP+ 35 s::Pre-Hmo-miR160a; A33-C33, 35 s::negative control Hpcyt P450-like2::eGFP+ 35 s::Pre-Hmo-miR160a; A34-C34, 35 s::Pre-Hmo-novel-2; A35-C35, 35 s::HmCYP83B1-like::eGFP; A36-C36, 35 s::positive control HmCYP83B1-like::eGFP; A37-C37, 35 s::negative control HmCYP83B1-like::eGFP; A38-C38, 35 s::HmCYP83B1-like::eGFP+ 35 s::Pre-miR164; A39-C39, 35 s::HmCYP83B1-like::eGFP+35 s::Pre-Hmo-novel-2; A40-C40, 35 s::positive control HmCYP83B1-like::eGFP+ 35 s::Pre-Hmo-novel-2; A41-C41, 35 s::negative control HmCYP83B1-like::eGFP+ 35 s::Pre-Hmo-novel-2; A42-C42, 35 s::Pre-Hmo-novel-15; A43-C43, 35 s::HmTPST-like::eGFP; A44-C44, 35 s::positive control HmTPST-like::eGFP; A45-C45, 35 s::negative control HmTPST-like::eGFP; A46-C46, 35 s::HmTPST-like::eGFP+ 35 s::Pre-miR164; A47-C47, 35 s::HmTPST-like::eGFP+ 35 s::Pre-Hmo-novel-15; A48-C48, 35 s::positive control HmTPST-like::eGFP+ 35 s::Pre-Hmo-novel-15; A49-C49, 35 s::negative control HmTPST-like::eGFP+ 35 s::Pre-Hmo-novel-15.

Additional file 7: Figure S7. Different fruit developmental stages of ‘Guanhuahong’ (A) and ‘Guanhuabai’ (B) pitayas. A1 and B1, 13 d; A2 and B2, 16 d; A3 and B3, 19 d; A4 and B4, 23 d; A5 and B5, 25 d; A6 and B6, 27 d; A7 and B7, 29 d. Bar = 2.0 cm.

Additional file 4: Figure S4. Number of miRNA family members in pitaya pulp.

Additional file 3: Figure S3. GC content distributions of assembled transcripts and genes. A, GC content distributions of transcripts; B, GC content distributions of genes.

Additional file 2: Figure S2. The length distributions of assembled transcripts and genes. A, The length distributions of transcripts; B, The length distributions of genes.

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 H2O2 and O2− 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.

During litchi (Litchi chinensis Sonn.) fruit ripening, two major physiological changes, degreening (Chl degradation) and pigmentation (anthocyanin biosynthesis), are visually apparent. However, the specific factor triggering this important transition is still unclear. In the present study, we found that endogenous ABA content increased sharply when Chl breakdown was initiated and the ABA level peaked just before the onset of anthocyanin accumulation, suggesting that ABA plays an important role during litchi fruit pigmentation. We characterized three ABSCISIC ACID RESPONSE ELEMENT-BINDING FACTORs (LcABF1/2/3) belonging to group A of the basic leucine zipper (bZIP) transcription factors previously shown to be involved in ABA signaling under abiotic stress. LcABF1 transcripts increased at the onset of Chl degradation, and the expression of LcABF3 accumulated in parallel with anthocyanin biosynthesis. In addition, dual luciferase and yeast one-hybrid assays indicated that LcABF1/2 recognized ABA-responsive elements in the promoter region of Chl degradation-related genes (PAO and SGR), while LcABF2/3 bound the promoter region of LcMYB1 and anthocyanin biosynthesis-related structural genes. Indeed, Nicotiana benthamiana leaves transiently expressing LcABF1/2 showed a senescence phenomenon with Chl degradation, and LcABF3 overexpression increased the accumulation of anthocyanin via activation of LcMYB1, which is the key determinant of anthocyanin biosynthesis. These data indicate that LcABF1/2/3 are important transcriptional regulators of ABA-dependent litchi fruit ripening involved in both Chl degradation and anthocyanin biosynthesis.