Transcriptome Analysis of Differentially Expressed Genes Involved in Proanthocyanidin Accumulation in the Rhizomes of Fagopyrum dibotrys and an Irradiation-Induced Mutant

Genomics and resequencing of Fagopyrum dibotrys from different geographic regions reveals species evolution and genetic diversity

Fagopyrum dibotrys, belonging to the family Polygonaceae and genus Fagopyrum, is used in traditional Chinese medicine and is rich in beneficial components, such as flavonoids. As its abundant medicinal value has become increasingly recognized, its excessive development poses a considerable challenge to wild germplasm resources, necessitating artificial cultivation and domestication. Considering these factors, a high-quality genome of F. dibotrys was assembled and the evolutionary relationships within Caryophyllales were compared, based on which 58 individual samples of F. dibotrys were re-sequenced. We found that the samples could be categorized into three purebred populations and regions distributed at distinct elevations. Our varieties were cultivated from the parental populations of the subpopulation in central Yunnan. F. dibotrys is speculated to have originated in the high-altitude Tibetan Plateau region, and that its combination with flavonoids can protect plants against ultraviolet radiation; this infers a subpopulation with a high accumulation of flavonoids. This study assembled a high-quality genome and provided a theoretical foundation for the future introduction, domestication, and development of cultivated varieties of F. dibotrys.

Gamma irradiation delays tomato (Solanum lycopersicum) ripening by inducing transcriptional changes

BACKGROUND

Tomato (Solanum lycopersicum) has a relatively short shelf life as a result of rapid ripening, limiting its transportability and marketability. Recently, gamma irradiation has emerged as a viable method for delaying tomato fruit ripening. Although few studies have shown that gamma irradiation delays the ripening of tomatoes, the underlying mechanism remains unknown. Therefore, the present study aimed to examine the effects of gamma irradiation on tomato fruit ripening and the underlying mechanisms using transcriptomics.

RESULTS

Following gamma irradiation, the total microbial count, weight loss, and decay rate of tomatoes significantly reduced during storage. Furthermore, the redness (a*), color change (∆E), and lycopene content of gamma-irradiated tomatoes decreased in a dose-dependent manner during storage. Moreover, gamma irradiation significantly upregulated the expression levels of genes associated with DNA, chloroplast, and oxidative damage repairs, whereas those of ethylene and auxin signaling-, ripening-, and cell wall metabolism-related, as well as carotenoid genes, were downregulated.

CONCLUSION

Gamma irradiation effectively delayed ripening by downregulating the expression of ripening-related genes and inhibiting microbial growth, which prevented decay and prolonged the shelf life of tomatoes. © 2023 The Authors. Journal of The Science of Food and Agriculture published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Succession of endophytic fungi and rhizosphere soil fungi and their correlation with secondary metabolites in Fagopyrum dibotrys

Golden buckwheat (Fagopyrum dibotrys, also known as F. acutatum) is a traditional edible herbal medicinal plant with a large number of secondary metabolites and is considered to be a source of therapeutic compounds. Different ecological environments have a significant impact on their compound content and medicinal effects. However, little is known about the interactions between soil physicochemical properties, the rhizosphere, endophytic fungal communities, and secondary metabolites in F. dibotrys. In this study, the rhizosphere soil and endophytic fungal communities of F. dibotrys in five different ecological regions in China were identified based on high-throughput sequencing methods. The correlations between soil physicochemical properties, active components (total saponins, total flavonoids, proanthocyanidin, and epicatechin), and endophytic and rhizosphere soil fungi of F. dibotrys were analyzed. The results showed that soil pH, soil N, OM, and P were significantly correlated with the active components of F. dibotrys. Among them, epicatechin, proanthocyanidin, and total saponins were significantly positively correlated with soil pH, while proanthocyanidin content was significantly positively correlated with STN, SAN, and OM in soil, and total flavone content was significantly positively correlated with P in soil. In soil microbes, Mortierella, Trechispora, Exophiala, Ascomycota_unclassified, Auricularia, Plectosphaerella, Mycena, Fungi_unclassified, Agaricomycetes_unclassified, Coprinellus, and Pseudaleuria were significantly related to key secondary metabolites of F. dibotrys. Diaporthe and Meripilaceae_unclassified were significantly related to key secondary metabolites in the rhizome. This study presents a new opportunity to deeply understand soil-plant-fungal symbioses and secondary metabolites in F. dibotrys, as well as provides a scientific basis for using biological fertilization strategies to improve the quality of F. dibotrys.

Comparative transcriptomic and metabolomic analyses reveal differences in flavonoid biosynthesis between PCNA and PCA persimmon fruit

The fruit of the persimmon (Diospyros kaki.) has high economic and nutritional value and is rich in flavonoids. Flavonoids are essential secondary metabolisms in plants. The association between persimmon astringency and changes in the proanthocyanidins (a flavonoid subclass) content is well-known. However, information on the relationships between different astringency types and other flavonoid subclasses and biosynthetic genes is more limited. In this study, an initial correlation analysis between total flavonoids and fruit astringency type, and KEGG analysis of metabolites showed that flavonoid-related pathways were linked to differences between mature pollination-constant non-astringent (PCNA) varieties ('Jiro' and 'Yohou') and pollination-constant astringent (PCA) fruit varieties ('Zhongshi5' and 'Huojing'). Based on these findings, variations in the expression of genes and metabolites associated with flavonoid biosynthesis were investigated between typical PCNA ('Jiro') and PCA ('Huojing') persimmons during fruit development. The flavonoid concentration in 'Huojing' fruit was significantly higher than that of 'Jiro' fruit, especially, in levels of proanthocyanin precursor epicatechin and anthocyanin cyanidin derivatives. Combined WGCNA and KEGG analyses showed that genes such as PAL, C4H, CHI, CHS, F3H, F3'5'H, FLS, DFR, ANR, ANS, and UF3GT in the phenylpropanoid and flavonoid biosynthesis pathways may be significant factors impacting the proanthocyanin precursor and anthocyanin contents. Moreover, interactions between the R2R3MYB (evm.TU.contig7272.598) and WD40 (evm.TU.contig3208.5) transcription factors were found to be associated with the above structural genes. These findings provide essential information on flavonoid biosynthesis and its regulation in the persimmon and lay a foundation for further investigation into how astringency types affect flavor components in PCNA and PCA persimmons.

Comparison of buckwheat genomes reveals the genetic basis of metabolomic divergence and ecotype differentiation

Golden buckwheat (Fagopyrum dibotrys or Fagopyrum cymosum) and Tartary buckwheat (Fagopyrum tataricum) belong to the Polygonaceae and the Fagopyrum genus is rich in flavonoids. Golden buckwheat is a wild relative of Tartary buckwheat, yet golden buckwheat is a traditional Chinese herbal medicine and Tartary buckwheat is a food crop. The genetic basis of adaptive divergence between these two buckwheats is poorly understood. Here, we assembled a high-quality chromosome-level genome of golden buckwheat and found a one-to-one syntenic relationship with the chromosomes of Tartary buckwheat. Two large inversions were identified that differentiate golden buckwheat and Tartary buckwheat. Metabolomic and genetic comparisons of golden buckwheat and Tartary buckwheat indicate an amplified copy number of FdCHI, FdF3H, FdDFR, and FdLAR gene families in golden buckwheat, and a parallel increase in medicinal flavonoid content. Resequencing of 34 wild golden buckwheat accessions across the two morphologically distinct ecotypes identified candidate genes, including FdMYB44 and FdCRF4, putatively involved in flavonoid accumulation and differentiation of plant architecture, respectively. Our comparative genomic study provides abundant genomic resources of genomic divergent variation to improve buckwheat with excellent nutritional and medicinal value.

The spatiotemporal regulations of epicatechin biosynthesis under normal flowering and the continuous inflorescence removal treatment in Fagopyrum dibotrys

BACKGROUND

Flowering is a critical physiological change that interferes with not only biomass yield but also secondary metabolism, such as the biosynthesis of flavonoids, in rhizome/root plants. The continuous inflorescence removal (CIR) treatment is frequently conducted to weaken this effect. Fagopyrum dibotrys (D.Don) H.Hara (Golden buckwheat) is a kind of rhizome medicinal plant rich in flavonoids and is widely used for the treatment of lung diseases. The CIR treatment is usually conducted in F. dibotrys because of its excessive reproductive growth. To uncover the molecular mechanisms, comprehensive analysis was performed using metabolome and transcriptome data obtained from normally bloomed and the CIR treated plants.

RESULTS

Metabolome results demonstrated that in the rhizomes of F. dibotrys, its bioactive compound called epicatechin has higher amount than most of the detected precursors. Compared with the normally bloomed plants, the level of epicatechin in the rhizomes of the CIR group increased by 25% at the withering stage. Based on 96 samples of the control and the CIR groups at 4 flowering stages for 4 tissues, RNA-Seq results revealed a 3 ~ 5 times upregulations of all the key enzyme genes involved in the biosynthesis of epicatechin in both time (from the bud stage to the withering stage) and spatial dimensions (from the top of branch to rhizome) under the CIR treatment compared to normal flowering. Integrated analysis of LC-MS/MS and transcriptome revealed the key roles of several key enzyme genes besides anthocyanidin reductase (ANR). A total of 93 transcription factors were identified to co-expressed with the genes in epicatechin biosynthetic pathway. The flowering activator SQUAMOSA promoter-binding protein like (SPLs) exhibited opposite spatiotemporal expression patterns to that of the epicatechin pathway genes; SPL3 could significantly co-express with all the key enzyme genes rather than the flowering repressor DELLA. Weighted gene co-expression network analysis (WGCNA) further confirmed the correlations among chalcone synthases (CHSs), chalcone isomerases (CHIs), ANRs, SPLs and other transcription factors.

CONCLUSIONS

SPL3 might dominantly mediate the effect of normal flowering and the CIR treatment on the biosynthesis of epicatechin in rhizomes mainly through the negative regulations of its key enzyme genes including CHS, CHI and ANR.

Transcriptome and HPLC Analysis Reveal the Regulatory Mechanisms of Aurantio-Obtusin in Space Environment-Induced Senna obtusifolia Lines

Senna obtusifolia is a famous medicinal plant that is widely used in Asian countries. Its seed plays an important role in the treatment of many diseases because it contains various anthraquinones and flavonoids. Our previous studies have indicated that three space environment-induced S. obtusifolia lines (SP-lines) i.e., QC10, QC29, and QC46, have higher seed yield and aurantio-obtusin (AO) content. However, the underlying mechanism of higher AO content in SP-lines is still unknown. Herein, transcriptome sequencing and HPLC were employed to analyze the differences between SP-lines and ground control (GC3) and elucidate the regulatory mechanisms of AO accumulation in SP-lines. The results show that 4002 differentially expressed genes (DEGs) were identified in SP-lines versus (vs.) GC3. DEGs in the QC10 vs. GC3, QC29 vs. GC3, and QC46 vs. GC3 comparisons were classified into 28, 36, and 81 GO terms and involved in 63, 74, and 107 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. KEGG pathway and gene expression analysis revealed that DEGs involved in anthraquinone pathways were significantly elevated in QC10 and QC46. Integrating the results of GO annotation, KEGG enrichment, and gene expression analysis, we propose that the elevated genes such as DAHPS, DHQS, and MenB enhance the metabolic flux in the anthraquinone pathway and promote AO content in QC10 and QC46. Taken together, this study elucidated the mechanism of AO content in SP-lines and provides valuable genetic information for S. obtusifolia. In addition, to the best of our knowledge, this study presents the first transcriptome analysis of environment-induced medicinal plants and paves the way to select elite S. obtusifolia varieties in the future.

Integrated Transcriptomics and Widely Targeted Metabolomics Analyses Provide Insights Into Flavonoid Biosynthesis in the Rhizomes of Golden Buckwheat (Fagopyrum cymosum)

Golden buckwheat (Fagopyrum cymosum) is used in Traditional Chinese Medicine. It has received attention because of the high value of its various medicinal and nutritional metabolites, especially flavonoids (catechin and epicatechin). However, the metabolites and their encoding genes in golden buckwheat have not yet been identified in the global landscape. This study performed transcriptomics and widely targeted metabolomics analyses for the first time on rhizomes of golden buckwheat. As a result, 10,191 differentially expressed genes (DEGs) and 297 differentially regulated metabolites (DRMs) were identified, among which the flavonoid biosynthesis pathway was enriched in both transcriptome and metabolome. The integration analyses of the transcriptome and the metabolome revealed a network related to catechin, in which four metabolites and 14 genes interacted with each other. Subsequently, an SG5 R2R3-MYB transcription factor, named FcMYB1, was identified as a transcriptional activator in catechin biosynthesis, as it was positively correlated to eight flavonoid biosynthesis genes in their expression patterns and was directly bound to the promoters of FcLAR2 and FcF3'H1 by yeast one hybrid analysis. Finally, a flavonoid biosynthesis pathway was proposed in the rhizomes of golden buckwheat, including 13 metabolites, 11 genes encoding 9 enzymes, and 1 MYB transcription factor. The expression of 12 DEGs were validated by qRT-PCR, resulting in a good agreement with the Pearson R ranging from 0.83 to 1. The study provided a comprehensive flavonoid biosynthesis and regulatory network of golden buckwheat.

Towards a better understanding of Fagopyrum dibotrys: a systematic review

Fagopyrum dibotrys (F. dibotrys) (D.Don) H.Hara is a well-known edible herbal medicine in Asian countries. It has been widely used for the treatment of lung diseases, swelling, etc., and is also an important part of many Chinese medicine prescriptions. At present, more than 100 compounds have been isolated and identified from F. dibotrys, and these compounds can be primarily divided into flavonoids, phenols, terpenes, steroids, and fatty acids. Flavonoids and phenolic compounds are considered to be the main active ingredients of F. dibotrys. Previous pharmacological studies have shown that F. dibotrys possesses anti-inflammatory, anti-cancer, anti-oxidant, anti-bacterial, and anti-diabetic activities. Additional studies on functional genes have led to a better understanding of the metabolic pathways and regulatory factors related with the flavonoid active ingredients in F. dibotrys. In this paper, we systemically reviewed the research advances on the phytochemistry and pharmacology of F. dibotrys, as well as the functional genes related to the synthesis of active ingredients, aiming to promote the development and utilization of F. dibotrys.

De novo transcriptome analysis unravels tissue-specific expression of candidate genes involved in major secondary metabolite biosynthetic pathways of Plumbago zeylanica: implication for pharmacological potential

KEY MESSAGE

The present study provides comparative transcriptome analysis, besides identifying functional secondary metabolite genes of Plumbago zeylanica with pharmacological potential for future functional genomics, and metabolomic engineering of secondary metabolites from this plant towards diversified biomedical applications.

ABSTRACT

Plumbago zeylanica is a widely used medicinal plant of the traditional Indian system of medicine with wide pharmacological potential to treat several disorders. The present study aimed to carry out comparative transcriptome analysis in leaf and root tissue of P. zeylanica using Illumina paired end sequencing to identify tissue-specific functional genes involved in the biosynthesis of secondary metabolites, contributing to its therapeutic efficacy. De novo sequencing assembly resulted in the identification of 62,321 "Unigenes" transcripts with an average size of 1325 bp. Functional annotation using BLAST2GO resulted in the identification of 50,301 annotated transcripts (80.71%) and GO assigned to 18,814 transcripts. KEGG pathway annotation of the "Unigenes" revealed that 2465 transcripts could be assigned to 242 KEGG pathway maps wherein the number of transcripts involved in secondary metabolism was distinct in root and leaf transcriptome. Among the secondary metabolite biosynthesis pathways, the cluster of "Unigenes" encoding enzymes of 'Phenylpropanoid biosynthesis pathway' represents the largest group (84 transcripts) followed by 'Terpenoid Backbone biosynthesis' (48 transcripts). The transcript levels of the candidate unigenes encoding key enzymes of phenylpropanoid (PAL, TAL) and flavanoid biosynthesis (CHS, ANS, FLS) pathways were up-regulated in root, while the expression levels of candidate "Unigenes" transcript for monoterpenoid (DXS, ISPF), diterpenoid biosynthesis (SPS, SDS) and indole alkaloid pathways (STR) were significantly higher in leaf of P. zeylanica. Interestingly, validation of differential gene expression profile by qRT-PCR also confirmed that candidate "Unigenes" enzymes of phenylpropanoid and flavonoid biosynthesis were highly expressed in the root, while the key regulatory enzymes of terpenoid and indole alkaloid compounds were up-regulated in the leaf, suggesting that (differences in) the levels of these functional genes could be attributed to the (differential) pharmacological activity (between root and leaf) in tissues of P. zeylanica.

Complete chloroplast genome sequence of Fagopyrum dibotrys: genome features, comparative analysis and phylogenetic relationships

Fagopyrum dibotrys, belongs to Polygonaceae family, is one of national key conserved wild plants of China with important medicinal and economic values. Here, the complete chloroplast (cp) genome sequence of F. dibotrys is reported. The cp genome size is 159,919 bp with a typical quadripartite structure and consisting of a pair of inverted repeat regions (30,738 bp) separated by large single copy region (85,134 bp) and small single copy region (13,309 bp). Sequencing analyses indicated that the cp genome encodes 131 genes, including 80 protein-coding genes, 28 tRNA genes and 4 rRNA genes. The genome structure, gene order and codon usage are typical of angiosperm cp genomes. We also identified 48 simple sequence repeats (SSR) loci, fewer of them are distributed in the protein-coding sequences compared to the noncoding regions. Comparison of F. dibotrys cp genome to other Polygonaceae cp genomes indicated the inverted repeats (IRs) and coding regions were more conserved than single copy and noncoding regions, and several variation hotspots were detected. Coding gene sequence divergence analyses indicated that five genes (ndhK, petL rpoC2, ycf1, ycf2) were subject to positive selection. Phylogenetic analysis among 42 species based on cp genomes and 50 protein-coding genes indicated a close relationship between F. dibotrys and F. tataricum. In summary, the complete cp genome sequence of F. dibotrys reported in this study will provide useful plastid genomic resources for population genetics and pave the way for resolving phylogenetic relationships of order Caryophyllales.

Gene Profiling in Late Blight Resistance in Potato Genotype SD20

Late blight caused by the oomycete fungus Phytophthora infestans (Pi) is the most serious obstacle to potato (Solanum tuberosum) production in the world. A super race isolate, CN152, which was identified from Sichuan Province, China, could overcome nearly all known late blight resistance genes and caused serious damage in China. The potato genotype SD20 was verified to be highly resistant to CN152; however, the molecular regulation network underlying late blight resistance pathway remains unclear in SD20. Here, we performed a time-course experiment to systematically profile the late blight resistance response genes using RNA-sequencing in SD20. We identified 3354 differentially expressed genes (DEGs), which mainly encoded transcription factors and protein kinases, and also included four NBS-LRR genes. The late blight responsive genes showed time-point-specific induction/repression. Multi-signaling pathways of salicylic acid, jasmonic acid, and ethylene signaling pathways involved in resistance and defense against Pi in SD20. Gene Ontology and KEGG analyses indicated that the DEGs were significantly enriched in metabolic process, protein serine/threonine kinase activity, and biosynthesis of secondary metabolites. Forty-three DEGs were involved in immune response, of which 19 were enriched in hypersensitive response reaction, which could play an important role in broad-spectrum resistance to Pi infection. Experimental verification confirmed the induced expression of the responsive genes in the late blight resistance signaling pathway, such as WRKY, ERF, MAPK, and NBS-LRR family genes. Our results provided valuable information for understanding late blight resistance mechanism of potato.

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

Buckwheat is an important minor crop with pharmaceutical functions due to rutin enrichment in the seed. Seeds of common buckwheat cultivars (Fagopyrum esculentum, Fes) usually have much lower rutin content than tartary buckwheat (F. tartaricum, Ft). We previously found a wild species of common buckwheat (F. esculentum ssp. ancestrale, Fea), with seeds that are high in rutin, similar to Ft. In the present study, we investigated the mechanism by which rutin production varies among different buckwheat cultivars, Fea, a Ft variety (Xide) and a Fes variety (No.2 Pingqiao) using RNA sequencing of filling stage seeds. Sequencing data generated approximately 43.78-Gb of clean bases, all these data were pooled together and assembled 180,568 transcripts, and 109,952 unigenes. We established seed gene expression profiles of each buckwheat sample and assessed genes involved in flavonoid biosynthesis, storage proteins production, CYP450 family, starch and sucrose metabolism, and transcription factors. Differentially expressed genes between Fea and Fes were further analyzed due to their close relationship than with Ft. Expression levels of flavonoid biosynthesis gene FLS1 (Flavonol synthase 1) were similar in Fea and Ft, and much higher than in Fes, which was validated by qRT-PCR. This suggests that FLS1 transcript levels may be associated with rutin accumulation in filling stage seeds of buckwheat species. Further, we explored transcription factors by iTAK, and multiple gene families were identified as being involved in the coordinate regulation of metabolism and development. Our extensive transcriptomic data sets provide a complete description of metabolically related genes that are differentially expressed in filling stage buckwheat seeds and suggests that FLS1 is a key controller of rutin synthesis in buckwheat species. FLS1 can effectively convert dihydroflavonoids into flavonol products. These findings provide a basis for further studies of flavonoid biosynthesis in buckwheat breeding to help accelerate flavonoid metabolic engineering that would increase rutin content in cultivars of common buckwheat.

Neuraminidase Inhibitory Activity and Constituent Characterization of Fagopyrum dibotrys

This study aimed to identify a new biological activity of the widely distributed species Fagopyrum dibotrys. Four F. dibotrys extracts (ethyl acetate (EA), petroleum ether (P), ethanol (E), and water (W)) were explored for their anti-neuraminidase (NA) activity. A total of 32 compounds were identified using UHPLC-Q-Exactive Orbitrap HRMS in the EA extract, which had the best NA inhibitory effects. We used the docking data for supporting compounds’ anti-neuraminidase activity. Among them, five compounds including one flavonoid, three organic acids, and one glucoside were discovered for the first time in F. dibotrys. Docking studies and NA activity assay revealed the remarkable NA inhibitory activity of eight components in EA extract, especially rutin, hesperidin, procyanidin B₂, and quercitrin. Therefore, F. dibotrys could be used to develop anti-influenza drugs.

Biorhizome: A Biosynthetic Platform for Colchicine Biomanufacturing

Colchicine is one of the oldest plant-based medicines used to treat gout and one of the most important alkaloid-based antimitotic drugs with anticancer potential, which is commercially extracted from Gloriosa superba. Clinical trials suggest that colchicine medication could prevent atrial fibrillation recurrence after cardiac surgery. In addition, therapeutic colchicine is undergoing clinical trials to treat non-diabetic metabolic syndrome and diabetic nephropathy. However, the industrial-scale biomanufacturing of colchicine have not yet been established. Clearly, further studies on detailed biorhizome-specific transcriptome analysis, gene expression, and candidate gene validation are required before uncover the mechanism of colchicine biosynthesis and biorhizome-based colchicine biomanufacturing. Annotation of 32312 assembled multiple-tissues transcripts of G. superba represented 15088 unigenes in known plant specific gene ontology. This could help understanding colchicine biosynthesis in G. superba. This review highlights the biorhizomes, rhizome specific genes or gene what expressed with high level in rhizomes, and deep fluid dynamics in a bioreactor specifically for the biomanufacture of colchicine.