Comparative Transcriptome Analysis Combining SMRT- and Illumina-Based RNA-Seq Identifies Potential Candidate Genes Involved in Betalain Biosynthesis in Pitaya Fruit

Multi-omic applications for understanding and enhancing tropical fruit flavour

Consumer trends towards nutrient-rich foods are contributing to global increasing demand for tropical fruit. However, commercial cultivars in the breeding pipeline that are tailored to meet market demand are at risk of possessing reduced fruit flavour qualities. This stems from recurrent prioritised selection for superior agronomic traits and not fruit flavour, which may in turn reduce consumer satisfaction. There is realisation that fruit quality traits, inclusive of flavour, must be equally selected for; but currently, there are limited tools and resources available to select for fruit flavour traits, particularly in tropical fruit species. Although sugars, acids, and volatile organic compounds are known to define fruit flavour, the specific combinations of these, that result in defined consumer preferences, remain unknown for many tropical fruit species. To define and include fruit flavour preferences in selective breeding, it is vital to determine the metabolites that underpin them. Then, objective quantitative analysis may be implemented instead of solely relying on human sensory panels. This may lead to the development of selective genetic markers through integrated omics approaches that target biosynthetic pathways of flavour active compounds. In this review, we explore progress in the development of tools to be able to strategically define and select for consumer-preferred flavour profiles in the breeding of new cultivars of tropical fruit species.

Effect of betanin synthesis on photosynthesis and tyrosine metabolism in transgenic carrot

BACKGROUND

Betalain is a natural pigment with important nutritional value and broad application prospects. Previously, we produced betanin biosynthesis transgenic carrots via expressing optimized genes CYP76AD1S, cDOPA5GTS and DODA1S. Betanin can accumulate throughout the whole transgenic carrots. But the effects of betanin accumulation on the metabolism of transgenic plants and whether it produces unexpected effects are still unclear.

RESULTS

The accumulation of betanin in leaves can significantly improve its antioxidant capacity and induce a decrease of chlorophyll content. Transcriptome and metabolomics analysis showed that 14.0% of genes and 33.1% of metabolites were significantly different, and metabolic pathways related to photosynthesis and tyrosine metabolism were markedly altered. Combined analysis showed that phenylpropane biosynthesis pathway significantly enriched the differentially expressed genes and significantly altered metabolites.

CONCLUSIONS

Results showed that the metabolic status was significantly altered between transgenic and non-transgenic carrots, especially the photosynthesis and tyrosine metabolism. The extra consumption of tyrosine and accumulation of betanin might be the leading causes.

Transcriptome analysis reveals modulation of differentially expressed genes in LPS-treated mouse macrophages (RAW264.7 cells) by grouper (Epinephelus coioides) Epinecidin-1

The marine antimicrobial peptide Epinecidin (Epi)-1 has been shown to exert direct antimicrobial and immunomodulatory actions in teleost, mammalian and avian organisms. For instance, Epi-1 can suppress bacterial endotoxin lipolysachcharide (LPS)-induced proinflammatory cytokines in RAW264.7 murine macrophages. However, it remains unknown how Epi-1 might broadly affect non-activated and LPS-activated macrophages. To address this question, we performed a comparative transcriptomic analysis of non-treated and LPS-treated RAW264.7 cells in the presence and absence of Epi-1. Gene enrichment analysis was conducted on filtered reads, followed by GO and KEGG analyses. The results showed that Epi-1 treatment modulated pathways and genes associated with nucleoside binding, intramolecular oxidoreductase activity, GTPase activity, peptide antigen binding, GTP binding, ribonucleoside/nucleotide binding, phosphatidylinositol binding and phosphatidylinositol-4-phosphate binding. Based on the GO analysis results, we performed real-time PCR at different treatment times to compare expression levels of selected proinflammatory cytokines, anti-inflammatory cytokines, MHC, proliferation and differentiation genes. Epi-1 decreased expression of the proinflammatory cytokines, TNF-α, IL-6 and IL-1β, and it increased the anti-inflammatory cytokine TGFβ and Sytx1. MHC-associated genes, GM7030, Arfip1, Gpb11 and Gem, were induced by Epi-1, which is expected to enhance the immune response against LPS. Immunoglobulin-associated Nuggc was also upregulated by Epi-1. Finally, we found that Epi-1 downregulated the expression of host defense peptides CRAMP, Leap2 and BD3. Taken together, these findings suggest that Epi-1 treatment induces orchestrated changes in the transcriptome of LPS-stimulated RAW264.7 cells.

Genome-Wide Analysis of microRNAs and Their Target Genes in Dongxiang Wild Rice (Oryza rufipogon Griff.) Responding to Salt Stress

Rice (Oryza sativa) is a staple food for more than half of the world's population, and its production is critical for global food security. Moreover, rice yield decreases when exposed to abiotic stresses, such as salinity, which is one of the most detrimental factors for rice production. According to recent trends, as global temperatures continue to rise due to climate change, more rice fields may become saltier. Dongxiang wild rice (Oryza rufipogon Griff., DXWR) is a progenitor of cultivated rice and has a high tolerance to salt stress, making it useful for studying the regulatory mechanisms of salt stress tolerance. However, the regulatory mechanism of miRNA-mediated salt stress response in DXWR remains unclear. In this study, miRNA sequencing was performed to identify miRNAs and their putative target genes in response to salt stress in order to better understand the roles of miRNAs in DXWR salt stress tolerance. A total of 874 known and 476 novel miRNAs were identified, and the expression levels of 164 miRNAs were found to be significantly altered under salt stress. The stem-loop quantitative real-time PCR (qRT-PCR) expression levels of randomly selected miRNAs were largely consistent with the miRNA sequencing results, suggesting that the sequencing results were reliable. The gene ontology (GO) analysis indicated that the predicted target genes of salt-responsive miRNAs were involved in diverse biological pathways of stress tolerance. This study contributes to our understanding of DXWR salt tolerance mechanisms regulated by miRNAs and may ultimately improve salt tolerance in cultivated rice breeding using genetic methods in the future.

Persistent Inflammation Initiated by TORCH Infections and Dysbiotic Microbiome in Autism Spectrum Disorders: A Prospect for Future Interventions

Autism spectrum disorders (ASD) are a range of neurodevelopmental conditions that are clinically present early in childhood with the symptoms of social withdrawal and repetitive behavior. Despite an extensive research on ASD, no commonly accepted theory on the disease etiology exists. Hence, we reviewed several scientific publications, including reviews, preclinical and clinical investigations, and published hypotheses to analyze various opinions on the nature and cause of the disorder. Many studies suggest that infections and inflammation during pregnancy play a significant role in genetic and epigenetic changes in the developing fetus, resulting in an autistic phenotype in a child. Still, there is a lack of comprehensive literature about the multitude of autism inducing factors. Therefore, this article reviews and discusses available scientific evidence on the roles of viral, bacterial, fungal, and parasitic infections, overactivation of the immune system, and intestinal microflora in the pathogenesis and clinical manifestation of ASD. The overview of the scientific publications, including our own studies, suggests that TORCH infections, imbalanced microbiome, and persistent inflammation are significantly associated with the disruption of the social domain in ASD children. The ASD-related changes begin prenatally as maternal-to-fetal immune activation triggered by infection. It results in continuous low-grade inflammation and oxidative stress in a fetus, causing germline and somatic genetic changes in the developing brain and the establishment of the dysregulated immune system. These changes and dysregulations result in central and peripheral nervous systems dysfunctions as well as other comorbid conditions found in autistic children.

Gene Co-expression Network Analysis of the Comparative Transcriptome Identifies Hub Genes Associated With Resistance to Aspergillus flavus L. in Cultivated Peanut (Arachis hypogaea L.)

Cultivated peanut (Arachis hypogaea L.), a cosmopolitan oil crop, is susceptible to a variety of pathogens, especially Aspergillus flavus L., which not only vastly reduce the quality of peanut products but also seriously threaten food safety for the contamination of aflatoxin. However, the key genes related to resistance to Aspergillus flavus L. in peanuts remain unclear. This study identifies hub genes positively associated with resistance to A. flavus in two genotypes by comparative transcriptome and weighted gene co-expression network analysis (WGCNA) method. Compared with susceptible genotype (Zhonghua 12, S), the rapid response to A. flavus and quick preparation for the translation of resistance-related genes in the resistant genotype (J-11, R) may be the drivers of its high resistance. WGCNA analysis revealed that 18 genes encoding pathogenesis-related proteins (PR10), 1-aminocyclopropane-1-carboxylate oxidase (ACO1), MAPK kinase, serine/threonine kinase (STK), pattern recognition receptors (PRRs), cytochrome P450, SNARE protein SYP121, pectinesterase, phosphatidylinositol transfer protein, and pentatricopeptide repeat (PPR) protein play major and active roles in peanut resistance to A. flavus. Collectively, this study provides new insight into resistance to A. flavus by employing WGCNA, and the identification of hub resistance-responsive genes may contribute to the development of resistant cultivars by molecular-assisted breeding.

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.

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.

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.

Selection and validation of reference genes for qRT-PCR analysis during fruit ripening of red pitaya (Hylocereus polyrhizus)

Red pitaya (Hylocereus polyrhizus) is widely cultivated in southern and southwestern China. To provide a basis for studying the molecular mechanisms of the ripening of this fruit, we carried out RNA sequencing (RNA-seq) analysis to identify differentially and stably expressed unigenes. The latter may serve as a resource of potential reference genes for normalization of target gene expression determined using quantitative real-time PCR (qRT-PCR). We selected 11 candidate reference genes from our RNA-seq analysis of red pitaya fruit ripening (ACT7, EF-1α, IF-4α, PTBP, PP2A, EF2, Hsp70, GAPDH, DNAJ, TUB and CYP), as well as β-ACT, which has been used as a reference gene for pitayas in previous studies. We then comprehensively evaluated their expression stability during fruit ripening using four statistical methods (GeNorm, NormFinder, BestKeeper and DeltaCt) and merged the four outputs using the online tool RefFinder for the final ranking. We report that PTBP and DNAJ showed the most stable expression patterns, whereas CYP and ACT7 showed the least stable expression patterns. The relative gene expression of red pitaya sucrose synthase and 4, 5-dihydroxyphenylalanine-extradiol-dioxygenase as determined by quantitative real-time PCR and normalized to PTBP and DNAJ was consistent with the RNA-seq results, suggesting that PTBP and DNAJ are suitable reference genes for studies of red pitaya fruit ripening.

Winemaking Characteristics of Red-Fleshed Dragon Fruit from Three Locations in Guizhou Province, China

The aim of this study was to identify the locations and harvest months in Guizhou province, China, producing the most suitable red dragon fruit (Hylocereus polyrhizus) for winemaking. Fruit from Guanling, Luodian and Zhenfeng counties was harvested separately from successive fruit cycles in August, September and October, respectively. The key traits measured were fruit weight, pulp yield, soluble solids content, and titratable acid. Wine characteristics measured were alcohol content, total carbohydrates, titratable acidity, volatile acidity, and betacyanin content. The overall suitability of fruit from each location for winemaking was evaluated using a multi-factor, unweighted, scorecard. On that basis, fruit from Guanling county harvested in August was the most suitable. Fruit from Luodian, and Zhenfeng was most suitable when harvested in August and September, and September, respectively. These results provide a preliminary guide for the sourcing of red dragon fruit from Guizhou for wine production.

Gene co-expression network analysis of the heat-responsive core transcriptome identifies hub genes in Brassica rapa

Gene co-expression network analysis of the heat-responsive core transcriptome in two contrasting Brassica rapa accessions reveals the main metabolic pathways, key modules and hub genes, are involved in long-term heat stress. Brassica rapa is a widely cultivated and economically important vegetable in Asia. High temperature is a common stress that severely impacts leaf head formation in B. rapa, resulting in reduced quality and production. The purpose of this study was thus to identify candidate heat tolerance genes by comparative transcriptome analysis of two contrasting B. rapa accessions in response to long-term heat stress. Two B. rapa accessions, '268' and '334', which showed significant differences in heat tolerance, were used for RNA sequencing analysis. We identified a total of 11,055 and 8921 differentially expressed genes (DEGs) in '268' and '334', respectively. Functional enrichment analyses of all of the identified DEGs, together with the genes identified from weighted gene co-expression network analyses (WGCNA), revealed that the autophagy pathway, glutathione metabolism, and ribosome biogenesis in eukaryotes were significantly up-regulated, whereas photosynthesis was down-regulated, in the heat resistance of B. rapa '268'. Furthermore, when B. rapa '334' was subjected to long-term high-temperature stress, heat stress caused significant changes in the expression of certain functional genes linked to protein processing in the endoplasmic reticulum and plant hormone signal transduction pathways. Autophagy-related genes might have been induced by persistent heat stress and remained high during recovery. Several hub genes like HSP17.6, HSP17.6B, HSP70-8, CLPB1, PAP1, PYR1, ADC2, and GSTF11 were discussed in this study, which may be potential candidates for further analyses of the response to long-term heat stress. These results should help elucidate the molecular mechanisms of heat stress adaptation in B. rapa.

Identification and expression profiling of genes involved in circadian clock regulation in red dragon fruit (Hylocereus polyrhizus) by full-length transcriptome sequencing

Crassulacean acid metabolism (CAM) plants fix CO2 at night, exhibiting a reversed regulatory pattern of metabolomic pathways compared with most model plants, which have C3 and C4 pathways. In this study, we used a valuable tropic fruit, red dragon fruit (Hylocereus polyrhizus), as model plant to identify and analyze the circadian regulation genes. Due to the absence of red dragon fruit's whole-genome dataset, we established a full-length transcriptome dataset using single-molecule real-time (SMRT) sequencing method. A 7.66-Gb dataset with 4,552,474 subreads was generated, with an average length of 1,683 bp and an N50 of 2,446 bp. Using this dataset, we identified center oscillator genes: CCA1 (CIRCADIAN CLOCK ASSOCIATED1), ELF3 (EARLY FLOWERING 3), GI (GIGANTEA), LHY (LATE ELONGATED HYPOCOTYL), LNK1 (NIGHT LIGHT-INDUCIBLE AND CLOCK-REGULATED 1), and TOC1 (TIMING OF CAB EXPRESSION 1); a gene for the input pathway: CRY1 (CRYPTOCHROME); a gene for the output pathway: CO (CONSTANS); and genes related to the CAM pathway: MDH (MALATE DEHYDROGENASE), ME (MALIC ENZYMES), and PPDK (PYRUVATE PHOSPHATE DIKINASE). We further established the 24-h rhythmic expression pattern of these genes and classified these into three groups: HpCCA1, HpELF3, HpLHY, HpLNK1, and HpGI have expression peaks during the day; HpTOC1, HpCO, and HpCRY1 have highest expression levels at night; The genes involved in the CAM pathways, namely, HpMDH, HpME1, and HpPPDK, have double expression peaks in the day and night. Comparison of these expression patterns between red dragon fruit and model plants could provide clues in understanding the circadian clock regulation and the activity of the CAM pathways in cactus plants.

Full-Length Transcriptome Analyses of Genes Involved in Triterpenoid Saponin Biosynthesis of Psammosilene tunicoides Hairy Root Cultures With Exogenous Salicylic Acid

Triterpenoid saponins constitute a diverse class of bioactive compounds in medicinal plants. Salicylic acid (SA) is an efficient elicitor for secondary metabolite production, but a transcriptome-wide regulatory network of SA-promoted triterpenoid saponin biosynthesis remains little understood. In the current study, we described the establishment of the hairy root culture system for Psammosilene tunicoides, a triterpenoid saponin-producing medicinal herb in China, using genetic transformation by Agrobacterium rhizogenes. Compared to controls, we found that total saponin content was dramatically increased (up to 2.49-fold) by the addition of 5 mg/L SA in hairy roots for 1 day. A combination of single-molecule real-time (SMRT) and next-generation sequencing (Illumina RNA-seq) was generated to analyze the full-length transcriptome data for P. tunicoides, as well as the transcript profiles in treated (8 and 24 h) and non-treated (0 h) groups with 5 mg/L SA in hairy roots. A total of 430,117 circular consensus sequence (CCS) reads, 16,375 unigenes and 4,678 long non-coding RNAs (lncRNAs) were obtained. The average length of unigenes (2,776 bp) was much higher in full-length transcriptome than that derived from single RNA-seq (1,457 bp). The differentially expressed genes (DEGs) were mainly enriched in the metabolic process. SA up-regulated the unigenes encoding SA-binding proteins and antioxidant enzymes in comparison with controls. Additionally, we identified 89 full-length transcripts encoding enzymes putatively involved in saponin biosynthesis. The candidate transcription factors (WRKY, NAC) and structural genes (AACT, DXS, SE, CYP72A) might be the key regulators in SA-elicited saponin accumulation. Their expression was further validated by quantitative real-time PCR (qRT-PCR). These findings preliminarily elucidate the regulatory mechanisms of SA on triterpenoid saponin biosynthesis in the transcriptomic level, laying a foundation for SA-elicited saponin augmentation in P. tunicoides.

Comparative transcriptome analysis reveals ectopic delta-5 and delta-6 desaturases enhance protective gene expression upon Vibrio vulnificus challenge in Tilapia (Oreochromis niloticus)

Background

Tilapia (Oreochromis niloticus) cultures are frequently infected by Vibrio vulnificus, causing major economic losses to production units. Previously, tilapia expressing recombinant delta-5 desaturase and delta-6 desaturase (D56) were found to be resistant to V. vulnificus infection. In this report, we profile the D56-mediated molecular changes underlying this resistance in tilapia. A comparative transcriptome analysis was performed on V. vulnificus-infected wild-type and D56-transgenic tilapia using Illumina’s sequencing-by-synthesis approach. Gene enrichment analysis on differentially expressed unigenes was performed, and the expression patterns were validated by real-time PCR.

Results

Comparative transcriptome analysis was performed on RNA-sequence profiles obtained from wild-type and D56-transgenic tilapia at 0, 6 and 24 h post-infection with V. vulnificaus. GO and KEGG gene enrichment analyses showed that D56 regulates several pathways and genes, including fatty acid (FA) metabolism associated, and inflammatory and immune response. Expression of selected FA metabolism-associated, inflammatory and immune responsive genes was validated by qPCR. The inflammatory and immune responsive genes that are modulated by FA-associated D56 likely contribute to the enhanced resistance against V. vulnificus infection in Tilapia.

Conclusions

Transcriptome profiling and filtering for two-fold change variation showed that 3795 genes were upregulated and 1839 genes were downregulated in D56-transgenic tilapia. These genes were grouped into pathways, such as FA metabolism, FA elongation, FA biosynthesis, biosynthesis of unsaturated FA, FA degradation, inflammation, immune response, and chemokines. FA-associated genes and immune-related genes were modulated by D56 at 6 h and 24 h post infection with V. vulnificus. The expression patterns of FA-related genes, inflammatory genes, antimicrobial peptide genes and immune responsive genes at 0, 3, 6, 12, 24 and 48 h post-infection suggests these genes are involved in the enhanced resistance of D56 transgenic tilapia to V. vulnificus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07521-5.

Elucidation of the core betalain biosynthesis pathway in Amaranthus tricolor

Amaranthus tricolor L., a vegetable Amaranthus species, is an economically important crop containing large amounts of betalains. Betalains are natural antioxidants and can be classified into betacyanins and betaxanthins, with red and yellow colors, respectively. A. tricolor cultivars with varying betalain contents, leading to striking red to green coloration, have been commercially produced. However, the molecular differences underlying betalain biosynthesis in various cultivars of A. tricolor remain largely unknown. In this study, A. tricolor cultivars with different colors were chosen for comparative transcriptome analysis. The elevated expression of AmCYP76AD1 in a red-leaf cultivar of A. tricolor was proposed to play a key role in producing red betalain pigments. The functions of AmCYP76AD1, AmDODAα1, AmDODAα2, and AmcDOPA5GT were also characterized through the heterologous engineering of betalain pigments in Nicotiana benthamiana. Moreover, high and low L-DOPA 4,5-dioxygenase activities of AmDODAα1 and AmDODAα2, respectively, were confirmed through in vitro enzymatic assays. Thus, comparative transcriptome analysis combined with functional and enzymatic studies allowed the construction of a core betalain biosynthesis pathway of A. tricolor. These results not only provide novel insights into betalain biosynthesis and evolution in A. tricolor but also provide a basal framework for examining genes related to betalain biosynthesis among different species of Amaranthaceae.

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

Dragon fruits are tropical fruits economically important for agricultural industries. As members of the family of Cactaceae, they have evolved to adapt to the arid environment. Here we report the draft genome of Hylocereus undatus, commercially known as the white-fleshed dragon fruit. The chromosomal level genome assembly contains 11 longest scaffolds corresponding to the 11 chromosomes of H. undatus. Genome annotation of H. undatus found ~29,000 protein-coding genes, similar to Carnegiea gigantea (saguaro). Whole-genome duplication (WGD) analysis revealed a WGD event in the last common ancestor of Cactaceae followed by extensive genome rearrangements. The divergence time between H. undatus and C. gigantea was estimated to be 9.18 MYA. Functional enrichment analysis of orthologous gene clusters (OGCs) in six Cactaceae plants found significantly enriched OGCs in drought resistance. Fruit flavor-related functions were overrepresented in OGCs that are significantly expanded in H. undatus. The H. undatus draft genome also enabled the discovery of carbohydrate and plant cell wall-related functional enrichment in dragon fruits treated with trypsin for a longer storage time. Lastly, genes of the betacyanin (a red-violet pigment and antioxidant with a very high concentration in dragon fruits) biosynthetic pathway were found to be co-localized on a 12 Mb region of one chromosome. The consequence may be a higher efficiency of betacyanin biosynthesis, which will need experimental validation in the future. The H. undatus draft genome will be a great resource to study various cactus plants.