Transcriptome analysis of filling stage seeds among three buckwheat species with emphasis on rutin accumulation

Rutin distribution in Tartary buckwheat: Identifying prime dietary sources through comparative analysis of post-processing treatments

Rutin is a crucial bioactive compound that determines the nutritional value of Tartary buckwheat (TB). However, the potential of utilizing TB as a dietary source of rutin for human consumption remains largely unexplored. This study aims to address these knowledge gaps by conducting a detailed analysis of rutin content distribution in TB tissues. Our findings revealed a significant variation in rutin content across different plant tissues. Notably, higher levels of rutin were found in embryos and cotyledons compared to other tissues, highlighting them as the primary sites of rutin accumulation in TB seeds and sprouts. Additional research on the processing of TB showed that sprouts and seeds retain high rutin levels even after boiling, steaming, deep-frying, stir-frying, and popping. Comparative analysis of different TB-derived products confirmed that cooked seeds and sprouts can serve as significant dietary sources of rutin. This study offers a foundational framework for the development of future dietary recommendations and applications of TB.

Buckwheat OMICS: present status and future prospects

Buckwheat (Fagopyrum spp.) is an underutilized resilient crop of North Western Himalayas belonging to the family Polygonaceae and is a source of essential nutrients and therapeutics. Common Buckwheat and Tatary Buckwheat are the two main cultivated species used as food. It is the only grain crop possessing rutin, an important metabolite with high nutraceutical potential. Due to its inherent tolerance to various biotic and abiotic stresses and a short life cycle, Buckwheat has been proposed as a model crop plant. Nutritional security is one of the major concerns, breeding for a nutrient-dense crop such as Buckwheat will provide a sustainable solution. Efforts toward improving Buckwheat for nutrition and yield are limited due to the lack of available: genetic resources, genomics, transcriptomics and metabolomics. In order to harness the agricultural importance of Buckwheat, an integrated breeding and OMICS platforms needs to be established that can pave the way for a better understanding of crop biology and developing commercial varieties. This, coupled with the availability of the genome sequences of both Buckwheat species in the public domain, should facilitate the identification of alleles/QTLs and candidate genes. There is a need to further our understanding of the molecular basis of the genetic regulation that controls various economically important traits. The present review focuses on: the food and nutritional importance of Buckwheat, its various omics resources, utilization of omics approaches in understanding Buckwheat biology and, finally, how an integrated platform of breeding and omics will help in developing commercially high yielding nutrient rich cultivars in Buckwheat.

Whole-Genome Sequencing and Analysis of Tumour-Forming Radish (Raphanus sativus L.) Line

Spontaneous tumour formation in higher plants can occur in the absence of pathogen invasion, depending on the plant genotype. Spontaneous tumour formation on the taproots is consistently observed in certain inbred lines of radish (Raphanus sativus var. radicula Pers.). In this paper, using Oxford Nanopore and Illumina technologies, we have sequenced the genomes of two closely related radish inbred lines that differ in their ability to spontaneously form tumours. We identified a large number of single nucleotide variants (amino acid substitutions, insertions or deletions, SNVs) that are likely to be associated with the spontaneous tumour formation. Among the genes involved in the trait, we have identified those that regulate the cell cycle, meristem activity, gene expression, and metabolism and signalling of phytohormones. After identifying the SNVs, we performed Sanger sequencing of amplicons corresponding to SNV-containing regions to validate our results. We then checked for the presence of SNVs in other tumour lines of the radish genetic collection and found the ERF118 gene, which had the SNVs in the majority of tumour lines. Furthermore, we performed the identification of the CLAVATA3/ESR (CLE) and WUSCHEL (WOX) genes and, as a result, identified two unique radish CLE genes which probably encode proteins with multiple CLE domains. The results obtained provide a basis for investigating the mechanisms of plant tumour formation and also for future genetic and genomic studies of radish.

Multi-omics identification of a key glycosyl hydrolase gene FtGH1 involved in rutin hydrolysis in Tartary buckwheat (Fagopyrum tataricum)

Rutin, a flavonoid rich in buckwheat, is important for human health and plant resistance to external stresses. The hydrolysis of rutin to quercetin underlies the bitter taste of Tartary buckwheat. In order to identify rutin hydrolysis genes, a 200 genotypes mini-core Tartary buckwheat germplasm resource was re-sequenced with 30-fold coverage depth. By combining the content of the intermediate metabolites of rutin metabolism with genome resequencing data, metabolite genome-wide association analyses (GWAS) eventually identified a glycosyl hydrolase gene FtGH1, which could hydrolyse rutin to quercetin. This function was validated both in Tartary buckwheat overexpression hairy roots and in vitro enzyme activity assays. Mutation of the two key active sites, which were determined by molecular docking and experimentally verified via overexpression in hairy roots and transient expression in tobacco leaves, exhibited abnormal subcellular localization, suggesting functional changes. Sequence analysis revealed that mutation of the FtGH1 promoter in accessions of two haplotypes might be necessary for enzymatic activity. Co-expression analysis and GWAS revealed that FtbHLH165 not only repressed FtGH1 expression, but also increased seed length. This work reveals a potential mechanism behind rutin metabolism, which should provide both theoretical support in the study of flavonoid metabolism and in the molecular breeding of Tartary buckwheat.

Metabolomics Combined with Correlation Analysis Revealed the Differences in Antioxidant Activities of Lotus Seeds with Varied Cultivars

Functional foods have potential health benefits for humans. Lotus seeds (LS) as functional foods have excellent antioxidant activities. However, the differences in chemical composition of different LS cultivars may affect their antioxidant activities. This study comprehensively analyzed the differences among five LS cultivars based on metabolomics and further revealed the effects of metabolites on antioxidant activities by correlation analysis. A total of 125 metabolites were identified in LS using UPLC-Q/TOF-MS. Then, 15 metabolites were screened as differential metabolites of different LS cultivars by chemometrics. The antioxidant activities of LS were evaluated by DPPH•, FRAP, and ABTS•+ assays. The antioxidant activities varied among different LS cultivars, with the cultivar Taikong 66 showing the highest antioxidant activities. The correlation analysis among metabolites and antioxidant activities highlighted the important contribution of phenolics and alkaloids to the antioxidant activities of LS. Particularly, 11 metabolites such as p-coumaric acid showed significant positive correlation with antioxidant activities. Notably, 6 differential metabolites screened in different LS cultivars showed significant effects on antioxidant activities. These results revealed the important effects of phytochemicals on the antioxidant activities of different LS cultivars. This study provided evidence for the health benefits of different LS cultivars.

Gene Cloning and Characterization of Transcription Factor FtNAC10 in Tartary Buckwheat (Fagopyrum tataricum (L.) Gaertn.)

NAC transcription factors play a significant role in plant stress responses. In this study, an NAC transcription factor, with a CDS of 792 bp encoding 263 amino acids, was cloned from Fagopyrum tataricum (L.) Gaertn. (F. tataricum), a minor cereal crop, which is rich in flavonoids and highly stress resistant. The transcription factor was named FtNAC10 (NCBI accession number: MK614506.1) and characterized as a member of the NAP subgroup of NAC transcriptions factors. The gene exhibited a highly conserved N-terminal, encoding about 150 amino acids, and a highly specific C-terminal. The resulting protein was revealed to be hydrophilic, with strong transcriptional activation activity. FtNAC10 expression occurred in various F. tataricum tissues, most noticeably in the root, and was regulated differently under various stress treatments. The over-expression of FtNAC10 in transgenic Arabidopsis thaliana (A. thaliana) seeds inhibited germination, and the presence of FtNAC10 enhanced root elongation under saline and drought stress. According to phylogenetic analysis and previous reports, our experiments indicate that FtNAC10 may regulate the stress response or development of F. tataricum through ABA-signaling pathway, although the mechanism is not yet known. This study provides a reference for further analysis of the regulatory function of FtNAC10 and the mechanism that underlies stress responses in Tartary buckwheat.

Transcriptome and metabolome analyses reveal the key genes related to grain size of big grain mutant in Tartary Buckwheat (Fagopyrum tartaricum)

Grain size with high heritability and stability is an important selection target during Tartary buckwheat breeding. However, the mechanisms that regulate Tartary buckwheat grain development are unknown. We generated transcriptome and metabolome sequencing from 10 and 15 days past anthesis (DPA) grains of big grain mutant (bg1) and WT, and identified 4108 differentially expressed genes (DEGs) including 93 significantly up-regulated differential genes and 85 significantly down-regulated genes in both stages, simultaneously. Meanwhile, we identified DEGs involved in ubiquitin-proteasome pathway, HAI-KU (IKU) pathway, mitogen-activated protein kinase (MAPK) signaling pathway, plant hormone (auxin, brassinosteroids and cytokinins) transduction pathway and five transcription factor families, including APETALA (AP2), GROWTH-REGULATING FACTORS (GRF), AUXIN RESPONSE FACTOR (ARF), WRKY and MYB. Weighted gene co-expression network analysis (WGCNA) was performed and obtained 9 core DEGs. Conjoint analyses of transcriptome and metabolome sequencing screened out 394 DEGs. Using a combined comprehensive analysis, we identified 24 potential candidate genes that encode E3 ubiquitin-protein ligase HIP1, EMBRYO-DEFECTIVE (EMB) protein, receptor-like protein kinase FERONIA (FER), kinesin-4 protein SRG1, and so on, which may be associated with the big-grain mutant bg1. Finally, a quantitative real-time Polymerase Chain Reaction (qRT-PCR) assay was conducted to validate the identified DEGs. Our results provide additional knowledge for identification and functions of causal candidate genes responsible for the variation in grain size and will be an invaluable resource for the genetic dissection of Tartary buckwheat high-yield molecular breeding.

Six Underutilized Grain Crops for Food and Nutrition in China

Underutilized grain crops are an essential part of the food system that supports humankind. A number of these crops can be found in China, such as barley, buckwheat, broomcorn millet, foxtail millet, oat, and sorghum, which have characteristics such as containing more nutritional elements, being resistant to biotic and abiotic stresses, and having strong adaptability to poor environments. The diversity of these crops provides options for farmers’ livelihoods and healthy food for the population. Although some mentioned crops such as barley, oat, and sorghum are not underutilized crops globally, they could be considered underutilized in China as they were more important in the past and could be revitalized for food and nutrition in the future. This paper reviews current progress in research and development in the areas of germplasm resource conservation, variety improvement, cultivation technologies, processing, and the nutrition and benefits of six underutilized grain crops in China. It is concluded that underutilized grain crops could play a critical role in food and nutritional security in China.

Dynamic transcriptome analysis suggests the key genes regulating seed development and filling in Tartary buckwheat (Fagopyrum tataricum Garetn.)

Tartary buckwheat is highly attractive for the richness of nutrients and quality, yet post-embryonic seed abortion greatly halts the yield. Seed development is crucial for determining grain yield, whereas the molecular basis and regulatory network of Tartary buckwheat seed development and filling is not well understood at present. Here, we assessed the transcriptional dynamics of filling stage Tartary buckwheat seeds at three developmental stages by RNA sequencing. Among the 4249 differentially expressed genes (DEGs), genes related to seed development were identified. Specifically, 88 phytohormone biosynthesis signaling genes, 309 TFs, and 16 expansin genes participating in cell enlargement, 37 structural genes involved in starch biosynthesis represented significant variation and were candidate key seed development genes. Cis-element enrichment analysis indicated that the promoters of differentially expressed expansin genes and starch biosynthesis genes are rich of hormone-responsive (ABA-, AUX-, ET-, and JA-), and seed growth-related (MYB, MYC and WRKY) binding sites. The expansin DEGs showed strong correlations with DEGs in phytohormone pathways and transcription factors (TFs). In total, phytohormone ABA, AUX, ET, BR and CTK, and related TFs could substantially regulate seed development in Tartary buckwheat through targeting downstream expansin genes and structural starch biosynthetic genes. This transcriptome data could provide a theoretical basis for improving yield of Tartary buckwheat.

Variation in Phenolic Compounds and Antioxidant Activity of Various Organs of African Cabbage (Cleome gynandra L.) Accessions at Different Growth Stages

The presence of nutritional and health-benefiting compounds has increased awareness of orphan leafy vegetables such as Cleome gynandra (CG), whose phytochemicals vary among accessions and organs during growth. This study investigated the polyphenol accumulation and antioxidant activities (AOA) of eight CG accessions from the vegetative stage to the seed set stage. Plants were separated into leaves and stem (LS), flowers, and silique organs, and extracts were analyzed for total phenolic content (TPC), total flavonoid content (TFC), rutin and astragalin content, and AOA using 2,2-diphenyl-1-picrylhydrazyl-hydrate (DPPH) and 2,2′-azino-bis(3-ethylbenzothiazoline-6-sulphonic acid) (ABTS). There were significant interaction effects of growth stages and accessions that contributed to changes in compounds content and AOA. TPC accumulated in plant generative parts, whereas flavonoids accumulated in young plant organs. HPLC profiling revealed that rutin was the most abundant compound in all organs, with flowers having the highest levels, while astragalin was only found in flowers. Silique extracts, particularly accession KF-14, recorded the highest TPC, which corresponded to the strongest radical scavenging activity in ABTS and DPPH assays and a strong linear correlation. The germplasm contained accessions with significantly different and varying levels of bioactive compounds and AOA. These findings potentiate the exploitation of CG organs such as siliques for AOA, flowers for rutin and astragalin, and young shoots for flavonoids. Moreover, the significant accumulation of the compounds in particular accessions of the germplasms suggest that such superior accessions may be useful candidates in genetic breeding programs to improve CG vegetable.

Transcriptome analysis reveals important candidate gene families related to oligosaccharides biosynthesis in Morindaofficinalis

Morinda officinalis How (MO) is one of the best-known traditional herbs and is widely cultivated in subtropical and tropical areas for many years, especially in southern China. Oligosaccharides are the major constituents in the roots of MO, which is well known for its therapeutic effects with anti-depression, anti-osteoporosis, memory-enhancing, ect. To date, the main gene families that regulate the biosynthetic pathway of MO oligosaccharides metabolism yet have been published. In our study, six cDNA libraries generated from six plants of MO were sequenced utilizing an Illumina HiSeq 4000 platform. Corresponding totals of more than 132.60 million clean reads were obtained from the six libraries and assembled into 25,812 unigenes with an average length of 1288 bp. Moreover, 6036 unigenes were found to be allocated to 26 pathways maps using several public databases, and 2538 differential expression genes (DEGs) were screened. Among them, 25 genes from three families were selected as the mainly candidate genes related to MO oligosaccharides biosynthesis. Then, the expression patterns of six DEGs closely related to MO oligosaccharides biosynthesis were verified by quantitative real-time PCR (qRT-PCR). Besides, the MO was clustered more closely to Coffea arabica of Rubiaceae. In summary, the transcriptomic analysis was used to investigate the differences in expression genes of oligosaccharides biosynthesis, with the notable outcome that several key gene families were closely linked to oligosaccharides biosynthesis.

Integrated microRNA and transcriptome profiling reveal key miRNA-mRNA interaction pairs associated with seed development in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

Tartary buckwheat seed development is an extremely complex process involving many gene regulatory pathways. MicroRNAs (miRNAs) have been identified as the important negative regulators of gene expression and performed crucial regulatory roles in various plant biological processes. However, whether miRNAs participate in Tartary buckwheat seed development remains unexplored.

RESULTS

In this study, we first identified 26 miRNA biosynthesis genes in the Tartary buckwheat genome and described their phylogeny and expression profiling. Then we performed small RNA (sRNA) sequencing for Tartary buckwheat seeds at three developmental stages to identify the miRNAs associated with seed development. In total, 230 miRNAs, including 101 conserved and 129 novel miRNAs, were first identified in Tartary buckwheat, and 3268 target genes were successfully predicted. Among these miRNAs, 76 exhibited differential expression during seed development, and 1534 target genes which correspond to 74 differentially expressed miRNAs (DEMs) were identified. Based on integrated analysis of DEMs and their targets expression, 65 miRNA-mRNA interaction pairs (25 DEMs corresponding to 65 target genes) were identified that exhibited significantly opposite expression during Tartary buckwheat seed development, and 6 of the miRNA-mRNA pairs were further verified by quantitative real-time polymerase chain reaction (qRT-PCR) and ligase-mediated rapid amplification of 5' cDNA ends (5'-RLM-RACE). Functional annotation of the 65 target mRNAs showed that 56 miRNA-mRNA interaction pairs major involved in cell differentiation and proliferation, cell elongation, hormones response, organogenesis, embryo and endosperm development, seed size, mineral elements transport, and flavonoid biosynthesis, which indicated that they are the key miRNA-mRNA pairs for Tartary buckwheat seed development.

CONCLUSIONS

Our findings provided insights for the first time into miRNA-mediated regulatory pathways in Tartary buckwheat seed development and suggested that miRNAs play important role in Tartary buckwheat seed development. These findings will be help to study the roles and regulatory mechanism of miRNAs in Tartary buckwheat seed development.

Mapping QTLs for 1000-grain weight and genes controlling hull type using SNP marker in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

Tartary buckwheat (Fagopyrum tataricum), an important pseudocereal crop, has high economic value due to its nutritional and medicinal properties. However, dehulling of Tartary buckwheat is difficult owing to its thick and tough hull, which has greatly limited the development of the Tartary buckwheat processing industry. The construction of high-resolution genetic maps serves as a basis for identifying quantitative trait loci (QTLs) and qualitative trait genes for agronomic traits. In this study, a recombinant inbred lines (XJ-RILs) population derived from a cross between the easily dehulled Rice-Tartary type and Tartary buckwheat type was genotyped using restriction site-associated DNA (RAD) sequencing to construct a high-density SNP genetic map. Furthermore, QTLs for 1000-grain weight (TGW) and genes controlling hull type were mapped in multiple environments.

RESULTS

In total, 4151 bin markers comprising 122,185 SNPs were used to construct the genetic linkage map. The map consisted of 8 linkage groups and covered 1444.15 cM, with an average distance of 0.35 cM between adjacent bin markers. Nine QTLs for TGW were detected and distributed on four loci on chromosome 1 and 4. A major locus detected in all three trials was mapped in 38.2-39.8 cM region on chromosome 1, with an LOD score of 18.1-37.0, and explained for 23.6-47.5% of the phenotypic variation. The genes controlling hull type were mapped to chromosome 1 between marker Block330 and Block331, which was closely followed by the major locus for TGW. The expression levels of the seven candidate genes controlling hull type present in the region between Block330 and Block336 was low during grain development, and no significant difference was observed between the parental lines. Six non-synonymous coding SNPs were found between the two parents in the region.

CONCLUSIONS

We constructed a high-density SNP genetic map for the first time in Tartary buckwheat. The mapped major loci controlling TGW and hull type will be valuable for gene cloning and revealing the mechanism underlying grain development and easy dehulling, and marker-assisted selection in Tartary buckwheat.

Comparative transcriptome analysis of Rheum australe, an endangered medicinal herb, growing in its natural habitat and those grown in controlled growth chambers

Rheum australe is an endangered medicinal herb of high altitude alpine region of Himalayas and is known to possess anti-cancerous properties. Unlike many herbs of the region, R. australe has broad leaves. The species thrives well under the environmental extremes in its niche habitat, therefore an understanding of transcriptome of R. australe to environmental cues was of significance. Since, temperature is one of the major environmental variables in the niche of R. australe, transcriptome was studied in the species growing in natural habitat and those grown in growth chambers maintained at 4 °C and 25 °C to understand genes associated with different temperatures. A total of 39,136 primarily assembled transcripts were obtained from 10,17,74,336 clean read, and 21,303 unigenes could match to public databases. An analysis of transcriptome by fragments per kilobase of transcript per million, followed by validation through qRT-PCR showed 22.4% up- and 22.5% down-regulated common differentially expressed genes in the species growing under natural habitat and at 4 °C as compared to those at 25 °C. These genes largely belonged to signaling pathway, transporters, secondary metabolites, phytohormones, and those associated with cellular protection, suggesting their importance in imparting adaptive advantage to R. australe in its niche.

Breeding Buckwheat for Increased Levels of Rutin, Quercetin and Other Bioactive Compounds with Potential Antiviral Effects

Common buckwheat (Fagopyrum esculentum Moench) and Tartary buckwheat (Fagopyrum tataricum (L.) Gaertn.) are sources of many bioactive compounds, such as rutin, quercetin, emodin, fagopyrin and other (poly)phenolics. In damaged or milled grain under wet conditions, most of the rutin in common and Tartary buckwheat is degraded to quercetin by rutin-degrading enzymes (e.g., rutinosidase). From Tartary buckwheat varieties with low rutinosidase activity it is possible to prepare foods with high levels of rutin, with the preserved initial levels in the grain. The quercetin from rutin degradation in Tartary buckwheat grain is responsible in part for inhibition of α-glucosidase in the intestine, which helps to maintain normal glucose levels in the blood. Rutin and emodin have the potential for antiviral effects. Grain embryos are rich in rutin, so breeding buckwheat with the aim of producing larger embryos may be a promising strategy to increase the levels of rutin in common and Tartary buckwheat grain, and hence to improve its nutritional value.

Strategic enhancement of genetic gain for nutraceutical development in buckwheat: A genomics-driven perspective

Buckwheat (Fagopyrum spp.) under the family Polygonaceae is an ancient pseudocereal with stupendous but less studied nutraceutical properties. The gluten free nature of protein, balanced amino acid profile and health promoting bioactive flavonoids make it a golden crop of future. Besides a scanty basic research, not much attention has been paid to the improvement of plant type and breeding of nutraceutical traits. Scanning of scientific literature indicates that adequate genetic variation exists for agronomic and nutritional traits in mainstream and wild gene pool of buckwheat. However, the currently employed conventional approaches together with poorly understood genetic mechanisms restrict effective utilization of the existing genetic variation in nutraceutical breeding of buckwheat. The latest trends in buckwheat genomics, particularly avalilabity of draft genome sequences for both the cultivated species (F. esculentum and F.tataricum) hold immense potential to overcome these limitations. Utilizing the transgenic hairy rot cultures, role of various transcription factors and gene families have been deduced in production and biosynthesis of bioactive flavonoids. Further, the acquisition of high-density genomics data coupled with the next-generation phenotyping will certainly improve our understanding of underlying genetic regulation of nutraceutical traits. The present paper highlights the application of multilayered omics interventions for tailoring a nutrient rich buckwheat cultivar and nutraceutical product development.

Transcriptome Analysis Reveals Key Seed-Development Genes in Common Buckwheat (Fagopyrum esculentum)

Seed development is an essential and complex process, which is involved in seed size change and various nutrients accumulation, and determines crop yield and quality. Common buckwheat (Fagopyrum esculentum Moench) is a widely cultivated minor crop with excellent economic and nutritional value in temperate zones. However, little is known about the molecular mechanisms of seed development in common buckwheat (Fagopyrum esculentum). In this study, we performed RNA-Seq to investigate the transcriptional dynamics and identify the key genes involved in common buckwheat seed development at three different developmental stages. A total of 4619 differentially expressed genes (DEGs) were identified. Based on the results of Gene Ontology (GO) and KEGG analysis of DEGs, many key genes involved in the seed development, including the Ca²⁺ signal transduction pathway, the hormone signal transduction pathways, transcription factors (TFs), and starch biosynthesis-related genes, were identified. More importantly, 18 DEGs were identified as the key candidate genes for seed size through homologous query using the known seed size-related genes from different seed plants. Furthermore, 15 DEGs from these identified as the key genes of seed development were selected to confirm the validity of the data by using quantitative real-time PCR (qRT-PCR), and the results show high consistency with the RNA-Seq results. Taken together, our results revealed the underlying molecular mechanisms of common buckwheat seed development and could provide valuable information for further studies, especially for common buckwheat seed improvement.

Transcriptome analysis of Asparagus officinalis reveals genes involved in the biosynthesis of rutin and protodioscin

Garden asparagus (Asparagus officinalis L.) is a popular vegetable cultivated worldwide. The secondary metabolites in its shoot are helpful for human health. We analyzed A. officinalis transcriptomes and identified differentially expressed genes (DEGs) involved in the biosynthesis of rutin and protodioscin, which are health-promoting functional compounds, and determined their association with stem color. We sequenced the complete mRNA transcriptome using the Illumina high-throughput sequencing platform in one white, three green, and one purple asparagus cultivars. A gene set was generated by de novo assembly of the transcriptome sequences and annotated using a BLASTx search. To investigate the relationship between the contents of rutin and protodioscin and their gene expression levels, rutin and protodioscin were analyzed using high-performance liquid chromatography. A secondary metabolite analysis using high-performance liquid chromatography showed that the rutin content was higher in green asparagus, while the protodioscin content was higher in white asparagus. We studied the genes associated with the biosynthesis of the rutin and protodioscin. The transcriptomes of the five cultivars generated 336 599 498 high-quality clean reads, which were assembled into 239 873 contigs with an average length of 694 bp, using the Trinity v2.4.0 program. The green and white asparagus cultivars showed 58 932 DEGs. A comparison of rutin and protodioscin biosynthesis genes revealed that 12 of the 57 genes associated with rutin and two of the 50 genes associated with protodioscin showed more than four-fold differences in expression. These DEGs might have caused a variation in the contents of these two metabolites between green and white asparagus. The present study is possibly the first to report transcriptomic gene sets in asparagus. The DEGs putatively involved in rutin and protodioscin biosynthesis might be useful for molecular engineering in asparagus.