No abstract

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.

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.

No abstract