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

Profiles of phenolics and their synthetic pathways in Asparagus officinalis L

The synthetic pathways of some phenolics compounds in asparagus have been reported, however, the diversified phenolics compounds including their modification and transcription regulation remains unknown. Thus, multi-omics strategies were applied to detect the phenolics profiles, contents, and screen the key genes for phenolics biosynthesis and regulation in asparagus. A total of 437 compounds, among which 204 phenolics including 105 flavonoids and 82 phenolic acids were detected with fluctuated concentrations in roots (Rs), spears (Ss) and flowering twigs (Fs) of the both green and purple cultivars. Based on the detected phenolics profiles and contents correlated to the gene expressions of screened synthetic enzymes and regulatory TFs, a full phenolics synthetic pathway of asparagus was proposed for the first time, essential for future breeding of asparagus and scaled healthy phenolics production using synthetic biological strategies.

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.

The quality evaluation of 30 Asparagus officinalis L. varieties

Asparagus, a vital economic contributor, is a well-liked vegetable grown around the globe, and some secondary metabolites in its spear are beneficial to human health. Asparagus spears possess a significant quantity of nutrients and phytochemicals; however, the difference in these chemical compositions among various varieties has not been sufficiently studied. This work aimed to detect the chemical compositions of 30 varieties of asparagus and to assess them by principal component analysis (PCA). The results showed that the contents of these chemical compositions varied in varieties. Selenium (Se, 1.12-2.9 μg/100 g dry-weight [DW]) was abundant in asparagus, with an average dry matter content of 8.25%. Free amino acids (5.60-9.98 g/100 g DW) and polyphenols (6.34-8.67 mg/g DW) were both present in high amounts, along with flavonoids (4.218-8.22 mg/g DW) and protodioscin (0.44-1.96 mg/g DW). Correlation analysis, PCA, and hierarchical cluster analysis were used to conduct a comprehensive evaluation of asparagus. Atlas, Appolo, Jinggang 111, Jingke 2, and WS-1 were the top five varieties with comprehensive scores. This study provided valuable data for the breeding, quality improvement, processing, and utilization of asparagus varieties in the future.

Steroidal saponin profiles and their key genes for synthesis and regulation in Asparagus officinalis L. by joint analysis of metabolomics and transcriptomics

BACKGROUND

Asparagus officinalis L. is a worldwide cultivated vegetable enrichened in both nutrient and steroidal saponins with multiple pharmacological activities. The upstream biosynthetic pathway of steroidal saponins (USSP) for cholesterol (CHOL) synthesis has been studied, while the downstream pathway of steroidal saponins (DSSP) starting from cholesterol and its regulation in asparagus remains unknown.

RESULTS

Metabolomics, Illumina RNAseq, and PacBio IsoSeq strategies were applied to different organs of both cultivated green and purple asparagus to detect the steroidal metabolite profiles & contents and to screen their key genes for biosynthesis and regulation. The results showed that there is a total of 427 compounds, among which 18 steroids were detected with fluctuated concentrations in roots, spears and flowering twigs of two garden asparagus cultivars. The key genes of DSSP include; steroid-16-hydroxylase (S16H), steroid-22-hydroxylase (S22H) and steroid-22-oxidase-16-hydroxylase (S22O-16H), steroid-26-hydroxylase (S26H), steroid-3-β-glycosyltransferase (S3βGT) and furostanol glycoside 26-O-beta-glucosidases (F26GHs) which were correlated with the contents of major steroidal saponins were screened, and the transcriptional factors (TFs) co-expressing with the resulted from synthetic key genes, including zinc fingers (ZFs), MYBs and WRKYs family genes were also screened.

CONCLUSIONS

Based on the detected steroidal chemical structures, profiles and contents which correlated to the expressions of screened synthetic and TFs genes, the full steroidal saponin synthetic pathway (SSP) of asparagus, including its key regulation networks was proposed for the first time.

Rutin alleviated acrolein-induced cytotoxicity in Caco-2 and GES-1 cells by forming a cyclic hemiacetal product

Acrolein (ACR), an α, β-unsaturated aldehyde, is a toxic compound formed during food processing, and the use of phenolics derived from dietary materials to scavenge ACR is a hot spot. In this study, rutin, a polyphenol widely present in various dietary materials, was used to investigate its capacity to scavenge ACR. It was shown that more than 98% of ACR was eliminated under the conditions of reaction time of 2 h, temperature of 80 °C, and molar ratio of rutin/ACR of 2/1. Further structural characterization of the formed adduct revealed that the adduct of rutin to ACR to form a cyclic hemiacetal compound (RAC) was the main scavenging mechanism. Besides, the stability of RAC during simulated in vitro digestion was evaluated, which showed that more than 83.61% of RAC was remained. Furthermore, the cytotoxicity of RAC against Caco-2 and GES-1 cells was significantly reduced compared with ACR, where the IC50 values of ACR were both below 20 μM while that of RAC were both above 140 μM. And the improvement of the loss of mitochondrial membrane potential (MMP) by RAC might be one of the detoxification pathways. The present study indicated that rutin was one of the potential ACR scavengers among natural polyphenols.

Effects of Mulching on Early-spring Green Asparagus Yield and Quality under Cultivation in Plastic Tunnels

Mulching significantly increases the crop yield and quality by positively affecting the physical features of the soil. The effect of multiple mulching treatments on soil temperature, yield, and nutrient composition of green asparagus grown in a plastic tunnel was assessed. Two mulch materials: transparent plastic film (PF) and rice husk (RH), were applied and compared with non-mulching treatment (CK). The soil temperatures at the soil surface and 10 cm below it were generally higher in the PF mulch than in the CK during the spring. PF mulch accelerated early spear emergence and growth, which led to harvesting 16 days earlier than in the CK. Under the PF mulch, the early yield of spears increased by 26.6% from January to the end of March, and the annual gross income by 14.8% because of the higher price resulting from the significantly higher marketable spear length and diameter; however, they were reduced in the RH mulch. Most nutrient compositions, such as soluble sugar, ascorbic acid, rutin, flavonoid content, and total antioxidant activity, were significantly increased in the PF mulch treatment. PF mulch might benefit green asparagus production during early spring because of its ability to promote early spear emergence and growth.

22R- but not 22S-hydroxycholesterol is recruited for diosgenin biosynthesis

Diosgenin is an important compound in the pharmaceutical industry and it is biosynthesized in several eudicot and monocot species, herein represented by fenugreek (a eudicot), and Dioscorea zingiberensis (a monocot). Formation of diosgenin can be achieved by the early C22,16-oxidations of cholesterol followed by a late C26-oxidation. This study reveals that, in both fenugreek and D. zingiberensis, the early C22,16-oxygenase(s) shows strict 22R-stereospecificity for hydroxylation of the substrates. Evidence against the recently proposed intermediacy of 16S,22S-dihydroxycholesterol in diosgenin biosynthesis was also found. Moreover, in contrast to the eudicot fenugreek, which utilizes a single multifunctional cytochrome P450 (TfCYP90B50) to perform the early C22,16-oxidations, the monocot D. zingiberensis has evolved two separate cytochrome P450 enzymes, with DzCYP90B71 being specific for the 22R-oxidation and DzCYP90G6 for the C16-oxidation. We suggest that the DzCYP90B71/DzCYP90G6 pair represent more broadly conserved catalysts for diosgenin biosynthesis in monocots.

Metabolomics Reveals Heterogeneity in the Chemical Composition of Green and White Spears of Asparagus (A. officinalis)

Green and white asparagus are quite different crops but can be harvested from the same plant. They have distinct morphological differences due to their mode of cultivation and they are characterised by having contrasting appearance and flavour. Significant chemical differences are therefore expected. Spears from three varieties of both green and white forms, harvested in two consecutive seasons were analysed using headspace GC-MS and LC-MS with an untargeted metabolomic workflow. Mainly C5 and C8 alcohols and aldehydes, and phenolic compounds were more abundant in green spears, whereas benzenoids, monoterpenes, unsaturated aldehydes and steroidal saponins were more abundant in white ones. Previously reported key asparagus volatiles and non-volatiles were detected at similar or not significantly different levels in the two asparagus types. Spatial metabolomics revealed also that many volatiles with known positive aroma attributes were significantly more abundant in the upper parts of the spears and showed a decreasing trend towards the base. These findings provide valuable insights into the metabolome of raw asparagus, the contrasts between green and white spears as well as the different chemical distributions along the stem.

Identification of key genes associated with secondary metabolites biosynthesis by system network analysis in Valeriana officinalis

Valeriana officinalis is a medicinal plant, a source of bioactive chemical compounds and secondary metabolites which are applied in pharmaceutical industries. The advent of ethnomedicine has provided alternatives for disease treatment and has increased demands for natural products and bioactive compounds. A set of preliminary steps to answers for such demands can include integrative omics for systems metabolic engineering, as an approach that contributes to the understanding of cellular metabolic status. There is a growing trend of this approach for genetically engineering metabolic pathways in plant systems, by which natural and synthetic compounds can be produced. As in the case of most medicinal plants, there are no sufficient information about molecular mechanisms involved in the regulation of metabolic pathways in V. officinalis. In this research, systems biology was performed on the RNA-seq transcriptome and metabolome data to find key genes that contribute to the synthesis of major secondary metabolites in V. officinalis. The R Package Weighted Gene Co-Expression Network Analysis (WGCNA) was employed to analyze the data. Based on the results, some major modules and hub genes were identified to be associated with the valuable secondary metabolites. In addition, some TF-encoding genes, including AP2/ERF-ERF, WRKY and NAC TF families, as well as some regulatory factors including protein kinases and transporters were identified. The results showed that several novel hub genes, such as PCMP-H24, RPS24B, ANX1 and PXL1, may play crucial roles in metabolic pathways. The current findings provide an overall insight into the metabolic pathways of V. officinalis and can expand the potential for engineering genome-scale pathways and systems metabolic engineering to increase the production of bioactive compounds by plants.

Tissue-specific signatures of metabolites and proteins in asparagus roots and exudates

Comprehensive untargeted and targeted analysis of root exudate composition has advanced our understanding of rhizosphere processes. However, little is known about exudate spatial distribution and regulation. We studied the specific metabolite signatures of asparagus root exudates, root outer (epidermis and exodermis), and root inner tissues (cortex and vasculature). The greatest differences were found between exudates and root tissues. In total, 263 non-redundant metabolites were identified as significantly differentially abundant between the three root fractions, with the majority being enriched in the root exudate and/or outer tissue and annotated as 'lipids and lipid-like molecules' or 'phenylpropanoids and polyketides'. Spatial distribution was verified for three selected compounds using MALDI-TOF mass spectrometry imaging. Tissue-specific proteome analysis related root tissue-specific metabolite distributions and rhizodeposition with underlying biosynthetic pathways and transport mechanisms. The proteomes of root outer and inner tissues were spatially very distinct, in agreement with the fundamental differences between their functions and structures. According to KEGG pathway analysis, the outer tissue proteome was characterized by a high abundance of proteins related to 'lipid metabolism', 'biosynthesis of other secondary metabolites' and 'transport and catabolism', reflecting its main functions of providing a hydrophobic barrier, secreting secondary metabolites, and mediating water and nutrient uptake. Proteins more abundant in the inner tissue related to 'transcription', 'translation' and 'folding, sorting and degradation', in accord with the high activity of cortical and vasculature cell layers in growth- and development-related processes. In summary, asparagus root fractions accumulate specific metabolites. This expands our knowledge of tissue-specific plant cell function.

Determination of the Geographical Origin of Asparagus officinalis L. by 1H NMR Spectroscopy

Food authenticity concerning the geographical origin becomes increasingly important for consumers, food industries, and food authorities. In this study, nontargeted ¹H NMR metabolomics combined with machine learning methodologies was applied to successfully distinguish the geographical origin of 237 samples of white asparagus from Germany, Poland, The Netherlands, Spain, Greece, and Peru. Support vector classification of the geographical origin achieved an accuracy of 91.5% for the entire sample set and 87.8% after undersampling the majority class. Important regions of the spectra could be identified and assigned to potential chemical markers. A subset of samples was compared to isotope-ratio mass spectrometry (IRMS), an established method for the determination of origin of white asparagus in Germany. Here, SVM classification led to accuracies of 79.4% for NMR and 70.9% for IRMS. Finally, the classification of asparagus from different German regions was evaluated, and the influence of year and variety was analyzed.

Genome-wide transcriptome variation landscape in Ruta chalepensis organs revealed potential genes responsible for rutin biosynthesis

Ruta chalepensis L., most commonly known as 'fringed rue,' is an excellent and valuable bioactive plant that produces a range of complex flavonoids, of which rutin is the major compound present in this plant of great pharmaceutical and medicinal significance. The present study is a pioneering attempt to examine the changes in the transcriptomic landscape of leaf, stem, and root tissues and correlate this with rutin quantity in each tissue in order to identify the candidate genes responsible for rutin biosynthesis and to increase genomic resources in fringed rue. Comparative transcriptome sequencing of leaves, stems and roots were performed using the NovaSeq 6000 platform. The de novo transcriptome assembly generated 254,685 transcripts representing 154,018 genes with GC content of 42.60 % and N50 of 2280 bp. Searching assembled transcripts against UniRef90 and SwissProt databases annotated 79.7 % of them as protein coding. The leaf tissues had the highest rutin content followed by stems and roots. Several differentially expressed genes and transcripts relating to rutin biosynthesis were identified in leaves comparing with roots or stems comparing with roots. All the genes known to be involved in rutin biosynthesis showed up-regulation in leaves as compared with roots. These results were confirmed by gene ontology (GO) and pathway enrichment analyses. Up-regulated genes in leaves as compared with roots enriched GO terms with relation to rutin biosynthesis e.g. action of flavonol synthase, biosynthetic mechanism of malonyl-CoA, and action of monooxygenase. Phylogenetic analysis of the rhamnosyltransferase (RT) gene showed that it was highly homologues with RT sequence from Citrus species and all were located in the same clade. This transcriptomic dataset will serve as an important public resource for future genomics and transcriptomic studies in R. chalepensis and will act as a benchmark for the identification and genetic modification of genes involved in the biosynthesis of secondary metabolites.

Comparative Metabolome and Transcriptome Analyses of Susceptible Asparagus officinalis and Resistant Wild A. kiusianus Reveal Insights into Stem Blight Disease Resistance

Phomopsis asparagi is one of the most serious fungal pathogens, which causes stem blight disease in Asparagus officinalis (AO), adversely affecting its production worldwide. Recently, the development of novel asparagus varieties using wild Asparagus genetic resources with natural P. asparagi resistance has become a priority in Japan due to the lack of resistant commercial AO cultivars. In this study, comparative metabolome and transcriptome analyses of susceptible AO and resistant wild Asparagus kiusianus (AK) 24 and 48 h postinoculated (AOI_24 hpi, AOI_48 hpi, AKI_24 hpi and AKI_48 hpi, respectively) with P. asparagi were conducted to gain insights into metabolic and expression changes associated with AK species. Following infection, the resistant wild AK showed rapid metabolic changes with increased levels of flavonoids and steroidal saponins and decreased asparagusic acid glucose ester content, compared with the susceptible AO plants. Transcriptome data revealed a total of 21  differentially expressed genes (DEGs) as the core gene set that displayed upregulation in the resistant AK versus susceptible AO after infection with P. asparagi. Kyoto Encyclopedia of Genes and Genomes pathway analysis of these DEGs identified 11 significantly enriched pathways, including flavonoid biosynthesis and primary metabolite metabolism, in addition to plant signaling and defense-related pathways. In addition, comparative single-nucleotide polymorphism and Indel distributions in susceptible AO and resistant AK plants were evaluated using the latest AO reference genome Aspof.V1. The data generated in this study are important resources for advancing Asparagus breeding programs and for investigations of genetic linkage mapping, phylogenetic diversity and plant defense-related genes.

Development of an Oriental Medicine Discrimination Method through Analysis of Steroidal Saponins in Dioscorea nipponica Makino and Their Anti-Osteosarcoma Effects

To prevent confusing Dioscorea nipponica (DN), an Oriental medicine, with Dioscorea quinquelobata (DQ) and Dioscorea septemloba (DS), a simple and accurate quantitative analysis method using HPLC combined with ultraviolet (UV) detection was developed and verified with UPLC-QTOF/MS through identification of five saponin glycosides: protodioscin (1), protogracillin (2), pseudoprotodioscin (3), dioscin (4), and gracillin (5). The newly developed analysis method showed sufficient reproducibility (<1.91%) and accuracy (92.1%–102.6%) and was able to identify DN based on the presence of compound 3 (13.821 ± 0.037 mg/mL) and the absence of 5. Compound 1, which is present in DN at a relatively high level (159.983 ± 0.064 mg/mL), was also an important marker for identification. Among the three species, DN showed the strongest activation of apoptotic signaling in osteosarcoma cells, while the four compounds detected in DN showed IC₅₀ values of 6.43 (1), 10.61 (2), 10.48 (3), and 6.90 (4). In conclusion, the strong inhibitory effect of DN against osteosarcoma was confirmed to be associated with 1 and 4, which is also related to the quantitative results. Therefore, the results of this study might provide important information for quality control related to Oriental medicine.