Integrated transcriptomic and WGCNA analyses reveal candidate genes regulating mainly flavonoid biosynthesis in Litsea coreana var. sinensis.

Litsea coreana Levl. var. sinensis (Allen) Yang et P. H. Huang is a popular ethnic herb and beverage plant known for its high flavonoid content, which has been linked to a variety of pharmacological benefits and crucial health-promoting impacts in humans. The progress in understanding the molecular mechanisms of flavonoid accumulation in this plant has been hindered due to the deficiency of genomic and transcriptomic resources. We utilized a combination of Illumina and Oxford Nanopore Technology (ONT) sequencing to generate a de novo hybrid transcriptome assembly. In total, 126,977 unigenes were characterized, out of which 107,977 were successfully annotated in seven public databases. Within the annotated unigenes, 3,781 were categorized into 58 transcription factor families. Furthermore, we investigated the presence of four valuable flavonoids-quercetin-3-O-ß-D-galactoside, quercetin-3-O-ß-D-glucoside, kaempferol-3-O-ß-D-galactoside, and kaempferol-3-O-ß-D-glucoside in 98 samples, using high-performance liquid chromatography. A weighted gene co-expression network analysis identified two co-expression modules, MEpink and MEturquoise, that showed strong positive correlation with flavonoid content. Within these modules, four transcription factor genes (R2R3-MYB, NAC, WD40, and ARF) and four key enzyme-encoding genes (CHI, F3H, PAL, and C4H) emerged as potential hub genes. Among them, the R2R3-MYB (LcsMYB123) as a homologous gene to AtMYB123/TT2, was speculated to play a significant role in flavonol biosynthesis based on phylogenetic analysis. Our findings provided a theoretical foundation for further research into the molecular mechanisms of flavonoid biosynthesis. Additionally, The hybrid transcriptome sequences will serve as a valuable molecular resource for the transcriptional annotation of L. coreana var. sinensis, which will contribute to the improvement of high-flavonoid materials.

Evolutionary analysis of MADS-box genes in buckwheat species and functional study of FdMADS28 in flavonoid metabolism.

The MADS-box gene family is a transcription factor family that is widely expressed in plants. It controls secondary metabolic processes in plants and encourages the development of tissues like roots and flowers. However, the phylogenetic analysis and evolutionary model of MADS-box genes in Fagopyrum species has not been reported yet. This study identified the MADS-box genes of three buckwheat species at the whole genome level, and conducted systematic evolution and physicochemical analysis. The results showed that these genes can be divided into four subfamilies, with fragment duplication being the main way for the gene family expansion. During the domestication process from golden buckwheat to tartary buckwheat and the common buckwheat, the Ka/Ks ratio indicated that most members of the family experienced strong purification selection pressure, and with individual gene pairs experiencing positive selection. In addition, we combined the expression profile data of the MADS genes, mGWAS data, and WGCNA data to mine genes FdMADS28/48/50 that may be related to flavonoid metabolism. The results also showed that overexpression of FdMADS28 could increase rutin content by decreasing Kaempferol pathway content in hairy roots, and increase the resistance and growth of hairy roots to PEG and NaCl. This study systematically analyzed the evolutionary relationship of MADS-box genes in the buckwheat species, and elaborated on the expression patterns of MADS genes in different tissues under biotic and abiotic stresses, laying an important theoretical foundation for further elucidating their role in flavonoid metabolism.

Genomic and transcriptomic studies on flavonoid biosynthesis in Lagerstroemia indica.

BACKGROUND: Lagerstroemia indica is a widely cultivated ornamental woody shrub/tree of the family Lythraceae that is used as a traditional medicinal plant in East Asia and Egypt. However, unlike other ornamental woody plants, its genome is not well-investigated, which hindered the discovery of the key genes that regulate important traits and the synthesis of bioactive compounds. RESULTS: In this study, the genomic sequences of L. indica were determined using several next-generation sequencing technologies. Altogether, 324.01 Mb sequences were assembled and 98.21% (318.21 Mb) of them were placed in 24 pseudo-chromosomes. The heterozygosity, repeated sequences, and GC residues occupied 1.65%, 29.17%, and 38.64% of the genome, respectively. In addition, 28,811 protein-coding gene models, 327 miRNAs, 552 tRNAs, 214 rRNAs, and 607 snRNAs were identified. The intra- and interspecies synteny and Ks analysis revealed that L. indica exhibits a hexaploidy. The co-expression profiles of the genes involved in the phenylpropanoid (PA) and flavonoid/anthocyanin (ABGs) pathways with the R2R3 MYB genes (137 members) showed that ten R2R3 MYB genes positively regulate flavonoid/anthocyanin biosynthesis. The colors of flowers with white, purple (PB), and deep purplish pink (DPB) petals were found to be determined by the levels of delphinidin-based (Dp) derivatives. However, the substrate specificities of LiDFR and LiOMT probably resulted in the different compositions of flavonoid/anthocyanin. In L. indica, two LiTTG1s (LiTTG1-1 and LiTTG1-2) were found to be the homologs of AtTTG1 (WD40). LiTTG1-1 was found to repress anthocyanin biosynthesis using the tobacco transient transfection assay. CONCLUSIONS: This study showed that the ancestor L. indica experienced genome triplication approximately 38.5 million years ago and that LiTTG1-1 represses anthocyanin biosynthesis. Furthermore, several genes such as LiDFR, LiOMTs, and R2R3 LiMYBs are related to anthocyanin biosynthesis. Further studies are required to clarify the mechanisms and alleles responsible for flower color development.

Integrative transcriptome and metabolome analysis reveals the discrepancy in the accumulation of active ingredients between Lycium barbarum cultivars.

MAIN CONCLUSION: The combined analysis of transcriptome and metabolome provided molecular insight into the dynamics of multiple active ingredients biosynthesis and accumulation across different cultivars of Lycium barbarum. Lycium barbarum L. has a high concentration of active ingredients and is well known in traditional Chinese herbal medicine for its therapeutic properties. However, there are many Lycium barbarum cultivars, and the content of active components varies, resulting in inconsistent quality between Lycium barbarum cultivars. At present, few research has been conducted to reveal the difference in active ingredient content among different cultivars of Lycium barbarum at the molecular level. Therefore, the transcriptome of 'Ningqi No.1' and 'Qixin No.1' during the three development stages (G, T, and M) was constructed in this study. A total of 797,570,278 clean reads were obtained. Between the two types of wolfberries, a total of 469, 2394, and 1531 differentially expressed genes (DEGs) were obtained in the 'G1 vs. G10,' 'T1 vs. T10,' and 'M1 vs. M10,' respectively, and were annotated with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology identifiers. Using these transcriptome data, most DEGs related to the metabolism of the active ingredients in 'Ningqi No.1' and 'Qixin No.1' were identified. Moreover, a widely targeted metabolome analysis of the metabolites of 'Ningqi 1' and 'Qixin 1' fruits at the maturity stage revealed 1,135 differentially expressed metabolites (DEMs) in 'M1 vs. M10,' and many DEMs were associated with active ingredients such as flavonoids, alkaloids, terpenoids, and so on. We further quantified the flavonoid, lignin, and carotenoid contents of the two Lycium barbarum cultivars during the three developmental stages. The present outcome provided molecular insight into the dynamics of multiple active ingredients biosynthesis and accumulation across different cultivars of Lycium barbarum, which would provide the basic data for the formation of Lycium barbarum fruit quality and the breeding of outstanding strains.

Integrated Transcriptome and Metabolome Analyses Reveal Bamboo Culm Color Formation Mechanisms Involved in Anthocyanin Biosynthetic in Phyllostachys nigra.

Phyllostachys nigra has green young culms (S1) and purple black mature culms (S4). Anthocyanins are the principal pigment responsible for color presentation in ornamental plants. We employ a multi-omics approach to investigate the regulatory mechanisms of anthocyanins in Ph. nigra. Firstly, we found that the pigments of the culm of Ph. nigra accumulated only in one to four layers of cells below the epidermis. The levels of total anthocyanins and total flavonoids gradually increased during the process of bamboo culm color formation. Metabolomics analysis indicated that the predominant pigment metabolites observed were petunidin 3-O-glucoside and malvidin O-hexoside, exhibiting a significant increase of up to 9.36-fold and 13.23-fold, respectively, during pigmentation of Ph. nigra culm. Transcriptomics sequencing has revealed that genes involved in flavonoid biosynthesis, phenylpropanoid biosynthesis, and starch and sucrose metabolism pathways were significantly enriched, leading to color formation. A total of 62 differentially expressed structural genes associated with anthocyanin synthesis were identified. Notably, PnANS2, PnUFGT2, PnCHI2, and PnCHS1 showed significant correlations with anthocyanin metabolites. Additionally, certain transcription factors such as PnMYB6 and PnMYB1 showed significant positive or negative correlations with anthocyanins. With the accumulation of sucrose, the expression of PnMYB6 is enhanced, which in turn triggers the expression of anthocyanin biosynthesis genes. Based on these findings, we propose that these key genes primarily regulate the anthocyanin synthesis pathway in the culm and contribute to the accumulation of anthocyanin, ultimately resulting in the purple-black coloration of Ph. nigra.

Haplotype-resolved genome assembly of Phanera championii reveals molecular mechanisms of flavonoid synthesis and adaptive evolution.

Phanera championii is a medicinal liana plant that has successfully adapted to hostile karst habitats. Despite extensive research on its medicinal components and pharmacological effects, the molecular mechanisms underlying the biosynthesis of critical flavonoids and its adaptation to karst habitats remain elusive. In this study, we performed high-coverage PacBio and Hi-C sequencing of P. championii, which revealed its high heterozygosity and phased the genome into two haplotypes: Hap1 (384.60 Mb) and Hap2 (383.70 Mb), encompassing a total of 58 612 annotated genes. Comparative genomes analysis revealed that P. championii experienced two whole-genome duplications (WGDs), with approximately 59.59% of genes originating from WGD events, thereby providing a valuable genetic resource for P. championii. Moreover, we identified a total of 112 genes that were strongly positively selected. Additionally, about 81.60 Mb of structural variations between the two haplotypes. The allele-specific expression patterns suggested that the dominant effect of P. championii was the elimination of deleterious mutations and the promotion of beneficial mutations to enhance fitness. Moreover, our transcriptome and metabolome analysis revealed alleles in different tissues or different haplotypes collectively regulate the synthesis of flavonoid metabolites. In summary, our comprehensive study highlights the significance of genomic and morphological adaptation in the successful adaptation of P. championii to karst habitats. The high-quality phased genomes obtained in this study serve as invaluable genomic resources for various applications, including germplasm conservation, breeding, evolutionary studies, and elucidation of pathways governing key biological traits of P. championii.

Identification of One O-Methyltransferase Gene Involved in Methylated Flavonoid Biosynthesis Related to the UV-B Irradiation Response in Euphorbia lathyris.

Flavonoids are ubiquitous polyphenolic compounds that play a vital role in plants' defense response and medicinal efficacy. UV-B radiation is a vital environmental regulator governing flavonoid biosynthesis in plants. Many plants rapidly biosynthesize flavonoids as a response to UV-B stress conditions. Here, we investigated the effects of flavonoid biosynthesis via UV-B irradiation in Euphorbia lathyris. We found that exposure of the E. lathyris callus to UV-B radiation sharply increased the level of one O-methyltransferase (ElOMT1) transcript and led to the biosynthesis of several methylated flavonoids. The methyltransferase ElOMT1 was expressed heterologously in E. coli, and we tested the catalytic activity of recombinant ElOMT1 with possible substrates, including caffeic acid, baicalin, and luteolin, in vitro. ElOMT1 could efficiently methylate when the hydroxyl groups were contained in the core nucleus of the flavonoid. This molecular characterization identifies a methyltransferase responsible for the chemical modification of the core flavonoid structure through methylation and helps reveal the mechanism of methylated flavonoid biosynthesis in Euphorbiaceae. This study identifies the O-methyltransferase that responds to UV-B irradiation and helps shed light on the mechanism of flavonoid biosynthesis in Euphorbia lathyris.

Depicting the genetic and metabolic panorama of chemical diversity in the tea plant.

As a frequently consumed beverage worldwide, tea is rich in naturally important bioactive metabolites. Combining genetic, metabolomic and biochemical methodologies, here, we present a comprehensive study to dissect the chemical diversity in tea plant. A total of 2837 metabolites were identified at high-resolution with 1098 of them being structurally annotated and 63 of them were structurally identified. Metabolite-based genome-wide association mapping identified 6199 and 7823 metabolic quantitative trait loci (mQTL) for 971 and 1254 compounds in young leaves (YL) and the third leaves (TL), respectively. The major mQTL (i.e., P < 1.05 × 10-5, and phenotypic variation explained (PVE) > 25%) were further interrogated. Through extensive annotation of the tea metabolome as well as network-based analysis, this study broadens the understanding of tea metabolism and lays a solid foundation for revealing the natural variations in the chemical composition of the tea plant. Interestingly, we found that galloylations, rather than hydroxylations or glycosylations, were the largest class of conversions within the tea metabolome. The prevalence of galloylations in tea is unusual, as hydroxylations and glycosylations are typically the most prominent conversions of plant specialized metabolism. The biosynthetic pathway of flavonoids, which are one of the most featured metabolites in tea plant, was further refined with the identified metabolites. And we demonstrated the further mining and interpretation of our GWAS results by verifying two identified mQTL (including functional candidate genes CsUGTa, CsUGTb, and CsCCoAOMT) and completing the flavonoid biosynthetic pathway of the tea plant.

A chromosome-scale genome sequence of Aeonium(Aeonium arboreum 'Velour') provides novel insights into the evolution of anthocyanin synthesis.

Anthocyanin glycoside is a water-soluble flavonoid compound that colors plants and aids in stress resistance. The mechanism driving the evolution of the anthocyanin synthesis pathway in plants remains unclear. Aeonium plants are highly regarded as model organisms for studying adaptive evolution. These plants can be categorized into various types, each distinguished by the content and distribution of anthocyanins in their leaves. The categories include red leaves, green leaves, black leaves, yellow leaves, and a classification known as the 'spot brocade series. In this study, we successfully assembled and annotated the genome of cultivar 'Aeonium arboreum 'Velour'' at chromosomal level. The genome size is 1,334.85 Mb containing 18 chromosomes in a single set, with a contig N50 of 23.47 Mb and a Scaffold N50 of 25.07 Mb. Through homology prediction, de novo prediction, and transcriptome prediction, we identified 166,228 coding genes, 161,656 of which were successfully annotated in the database. Comparative genomic analysis revealed that Aeonium arboreum 'Velour' underwent an independent genome-wide replication event after differentiating from Sedum album, Kalanchoe laxiflora, and Kalanchoe fedtschenkoi. It also shared a genome-wide replication event with Sedum album and Kalanchoe laxiflora. Aeonium arboreum 'Velour' exhibits a higher number of multi-copy gene families compared to other species. A total of 5,129 gene families unique to Aeonium arboreum 'Velour' were identified, primarily enriched in various metabolic pathways, including monoterpenoid biosynthesis, sesquiterpene and triterpene biosynthesis, cyanamide acid metabolism, flavonoid and flavonol biosynthesis, phosphonate and phosphinate metabolism, fatty acid degradation, biosynthesis of unsaturated fatty acid, ether lipid metabolism, tyrosine metabolism, and isoflavone biosynthesis according to the KEGG pathway analysis. Aeonium arboreum 'Velour' and Sedum album diversion dates back to approximately 43.11 million years ago during the Paleogene period, marked by the expansion of 2,807 gene families. In contrast, the divergence from Kalanchoe laxiflora and Kalanchoe fedtschenkoi began around 57.28 million years ago, with 219 gene families expanding. GO analysis highlighted that most of the expansion or contraction gene families were predominantly enriched in flower organs, leaf organ development, anthocyanin metabolism regulation, and light energy absorption and utilization. Remarkably, anthocyanin metabolism regulation is enriched to 80 expanded genes, including 36 bHLH transcription factors, possibly functioning as photosensitive pigment interaction factors (PIFs). We speculate that flavonoids play a pivotal role in the adaptation of Aeonium arboreum 'Velour' to environmental stress. Moreover, the evolution of the anthocyanin synthesis pathway is potentially driven by the plant's capability to absorb and utilize light energy, especially in high CO2 and high-temperature settings characteristic of the early Paleogene.

Full-length transcriptome sequencing provides new insights into the complexity of flavonoid biosynthesis in Glechoma longituba.

Glechoma longituba has been frequently used in treating urolithiasis and cholelithiasis due to the presence of flavonoids, which are its major bioactive constituents. However, research on the molecular background of flavonoid biosynthesis in G. longituba is limited. In this study, we used single-molecule real-time combined with next-generation sequencing technologies to construct the complete transcriptome of G. longituba. We identified 404,648 non-redundant transcripts, including 249,697 coding sequences, 197,811 simple sequence repeats, 174,846 long noncoding RNA, and 176,554 coding RNA. Moreover, we functionally annotated 346,218 isoforms (85.56%) and identified 86,528 differentially expressed genes. We also identified 55 non-redundant full-length isoforms related to the flavonoid biosynthetic pathway. Pearson correlation analysis revealed that the expression levels of some key genes of the flavonoid biosynthesis pathway were significantly positively correlated with the flavonoid metabolites. Furthermore, we performed bioinformatics analysis (sequence and structural) of isoform_47029 (encoding flavanone 3-hydroxylase) and isoform_53692 (encoding flavonol synthase) to evaluate their potential biological functions. Finally, we validated gene expression levels of 12 flavonoid-related key enzyme genes using quantitative real-time PCR. Overall, this study provides full-length transcriptome information on G. longituba for the first time and valuable molecular resources for further research on the medicinal properties of this plant.

De novo biosynthesis of anticarcinogenic icariin in engineered yeast.

Icariin (ICA) has wide applications in nutraceuticals and medicine with strong anticancer activities. However, the structural complexity and low abundance in plants of ICA lead to the unsustainable and high-cost supply from chemical synthesis and plant extraction. Here, the whole biosynthesis pathway of ICA was elucidated, then was constructed in Saccharomyces cerevisiae, including a 13-step heterologous ICA pathway from eleven kinds of plants as well as deletions or overexpression of ten yeast endogenous genes. Spatial regulation of 8-C-prenyltransferase to mitochondria and three-stage sequential control of 4'-O-methyltransferase, 3-OH rhamnosyltransferase, and 7-OH glycosyltransferase expression successfully achieved the de novo synthesis of ICA with a titer of 130 µg/L under shake-flask culture. The ICA synthesis from glucose represents the longest reconstructed pathway of flavonoid in microbe so far. This study provides a potential choice for the sustainable microbial production of number of complex flavonoids.

Phylogenomic Analysis of Cytochrome P450 Gene Superfamily and Their Association with Flavonoids Biosynthesis in Peanut (Arachis hypogaea L.).

Cytochrome P450s (CYPs) constitute extensive enzyme superfamilies in the plants, playing pivotal roles in a multitude of biosynthetic and detoxification pathways essential for growth and development, such as the flavonoid biosynthesis pathway. However, CYPs have not yet been systematically studied in the cultivated peanuts (Arachis hypogaea L.), a globally significant cash crop. This study addresses this knowledge deficit through a comprehensive genome-wide analysis, leading to the identification of 589 AhCYP genes in peanuts. Through phylogenetic analysis, all AhCYPs were systematically classified into 9 clans, 43 gene families. The variability in the number of gene family members suggests specialization in biological functions. Intriguingly, both tandem duplication and fragment duplication events have emerged as pivotal drivers in the evolutionary expansion of the AhCYP superfamily. Ka/Ks analysis underscored the substantial influence of strong purifying selection on the evolution of AhCYPs. Furthermore, we selected 21 genes encoding 8 enzymes associated with the flavonoid pathway. The results of quantitative real-time PCR (qRT-PCR) experiments unveiled stage-specific expression patterns during the development of peanut testa, with discernible variations between pink and red testa. Importantly, we identified a direct correlation between gene expression levels and the accumulation of metabolites. These findings offer valuable insights into elucidating the comprehensive functions of AhCYPs and the underlying mechanisms governing the divergent accumulation of flavonoids in testa of different colors.

Genome-Wide Association Study Reveals the Genetic Basis of Total Flavonoid Content in Brown Rice.

Flavonoids have anti-inflammatory, antioxidative, and anticarcinogenic effects. Breeding rice varieties rich in flavonoids can prevent chronic diseases such as cancer and cardio-cerebrovascular diseases. However, most of the genes reported are known to regulate flavonoid content in leaves or seedlings. To further elucidate the genetic basis of flavonoid content in rice grains and identify germplasm rich in flavonoids in grains, a set of rice core collections containing 633 accessions from 32 countries was used to determine total flavonoid content (TFC) in brown rice. We identified ten excellent germplasms with TFC exceeding 300 mg/100 g. Using a compressed mixed linear model, a total of 53 quantitative trait loci (QTLs) were detected through a genome-wide association study (GWAS). By combining linkage disequilibrium (LD) analysis, location of significant single nucleotide polymorphisms (SNPs), gene expression, and haplotype analysis, eight candidate genes were identified from two important QTLs (qTFC1-6 and qTFC9-7), among which LOC_Os01g59440 and LOC_Os09g24260 are the most likely candidate genes. We also analyzed the geographic distribution and breeding utilization of favorable haplotypes of the two genes. Our findings provide insights into the genetic basis of TFC in brown rice and could facilitate the breeding of flavonoid-rich varieties, which may be a prevention and adjuvant treatment for cancer and cardio-cerebrovascular diseases.

Characterization and expression analysis of chalcone synthase gene family members suggested their roles in the male sterility of a wheat temperature-sensitive genic male sterile (TGMS) line.

Chalcone synthase (CHS), also known as the plants-specific type III polyketide synthases (PKSs), catalyzes the first key step in the biosynthesis of plant flavonoids. Flavonoids are one of the most important secondary metabolites which participate in flower pigmentation and pollen fertility. Recent reports have demonstrated the role of the CHS family in plant pollen exine formation. This study focused on the potential roles of CHS in the pollen exine formation of wheat. In the present study, a genome-wide investigation of the CHS family was carried out, and 87 CHS genes were identified in wheat. TaCHS3, TaCHS10, and TaCHS13 are wheat orthologs of Arabidopsis LESS ADHESIVE POLLEN (LAP5); TaCHS58, TaCHS64, and TaCHS67 are wheat orthologs of AtLAP6. TaCHS3, TaCHS10, and TaCHS67 showed anther-specific patterns. The expression of TaCHS3, TaCHS10, and TaCHS67 was positively co-expressed with sporopollenin biosynthetic genes, including TaCYP703A2, TaCYP704B1, TaDRL1, TaTKPR2, and TaMS2. Coincidently, the expression of TaCHS3, TaCHS10, and TaCHS67, together with those sporopollenin biosynthetic genes, were repressed at the tetrads and uninucleate stages in the temperature-sensitive genic male-sterile (TGMS) line BS366 under sterile conditions. Wheat anther-specific CHS genes might participate in the exine formation of BS366 through co-expressing with sporopollenin biosynthetic genes, which will undoubtedly provide knowledge of the roles of CHS in wheat pollen development.

Integrative analysis of the metabolome and transcriptome reveals the molecular regulatory mechanism of isoflavonoid biosynthesis in Ormosia henryi Prain.

Flavonoids are important components of many phytopharmaceuticals, however, most studies on flavonoids and isoflavonoids have been conducted on herbaceous plants of the family Leguminosae, such as soybean, and less attention has been paid to woody plants. To fill this gap, we characterized the metabolome and transcriptome of five plant organs of Ormosia henryi Prain (OHP), a woody Leguminosae plant with great pharmaceutical value. Our results indicate that OHP possesses a relatively high content of isoflavonoids as well as significant diversity, with greater diversity of isoflavonoids in the roots. Combined with transcriptome data, the pattern of isoflavonoid accumulation was found to be highly correlated with differential expression genes. Furthermore, the use of trait-WGCNA network analysis identified OhpCHSs as a probable hub enzyme that directs the downstream isoflavonoid synthesis pathway. Transcription factors, such as MYB26, MYB108, WRKY53, RAV1 and ZFP3, were found to be involved in the regulation of isoflavonoid biosynthesis in OHP. Our findings will be beneficial for the biosynthesis and utilization of woody isoflavonoids.

Full-length transcriptome analysis provides insights into flavonoid biosynthesis in Ranunculus japonicus.

Ranunculus japonicus Thunb. is a traditional Chinese herb. Plants in the genus Ranunculus are generally rich in flavonoids, which have antibacterial, anti-infective, and other pharmacological effects. However, owing to the lack of reference genomes, little is known about the flavonoid biosynthetic pathway in R. japonicus. In this study, PacBio isoform sequencing (PacBio iso-seq) and DNA nanoball sequencing (DNB-seq) were combined to build a full-length transcriptome database for three different tissues of R. japonicus. A total of 395,402 full-length transcripts were obtained, of which 308,474 were successfully annotated. A Kyoto Encyclopedia of Genes and Genomes analysis identified 29 differentially expressed genes encoding nine key enzymes for flavonoid biosynthesis. Correlation analysis indicated that flavanone 3-hydroxylase and flavonol synthase genes might have key roles in the accumulation of flavonoid substances in the different tissues of R. japonicus. The structures of chalcone synthase and chalcone isomerase enzymes were spatially modeled. Reverse-transcription quantitative PCR was used to verify gene expression levels of key enzymes associated with flavonoid biosynthesis. In addition, 22 MYB transcription factors involved in flavonoid biosynthesis and phenylpropanoid biosynthesis were discovered. The reliable transcriptomic data from this study provide genetic information about R. japonicus as well as insights into the molecular mechanism of flavonoid biosynthesis. The results also provide a basis for developing the medicinal value R. japonicus.

The transcription factors SbMYB45 and SbMYB86.1 regulate flavone biosynthesis in Scutellaria baicalensis.

Scutellaria baicalensis Georgi is an important Chinese medicinal plant that is rich in the flavones baicalin, wogonoside, and wogonin, providing it with anti-cancer, anti-inflammatory, and antibacterial properties. However, although the biosynthetic pathways of baicalin and its derivates have been elucidated, the regulation of flavone biosynthesis in S. baicalensis is poorly understood. Here, we found that the contents of baicalin and its derivates increased and that baicalin biosynthetic pathway genes were induced in response to light, and baicalin and baicalein are not exclusively produced in the roots of S. baicalensis. Based on the fact that MYB transcription factors are known to play important roles in flavone biosynthesis, we identified SbMYB45 and SbMYB86.1 in S. baicalensis and determined that they bind to the promoter of the flavone biosynthesis gene SbCHI to enhance its transcription. Moreover, overexpressing SbMYB45 and SbMYB86.1 enhanced the accumulation of baicalin in S. baicalensis leaves. We demonstrate that SbMYB45 and SbMYB86.1 bind to the cis-acting element MBSII in the promoter of CHI to redundantly induce its expression upon light exposure. These findings indicate that SbMYB45 and SbMYB86.1 transcriptionally activate SbCHI in response to light and enhance flavone contents in S. baicalensis.

Construction of the first high-density genetic linkage map and QTL mapping of flavonoid and leaf-size related traits in Epimedium.

BACKGROUND: Leaves are the main medicinal organ in Epimedium herbs, and leaf flavonoid content is an important criterion of Epimedium herbs. However, the underlying genes that regulate leaf size and flavonoid content are unclear, which limits the use of breeding for Epimedium development. This study focuses on QTL mapping of flavonoid and leaf-size related traits in Epimedium. RESULTS: We constructed the first high-density genetic map (HDGM) using 109 F1 hybrids of Epimedium leptorrhizum and Epimedium sagittatum over three years (2019-2021). Using 5,271 single nucleotide polymorphism (SNP) markers, an HDGM with an overall distance of 2,366.07 cM and a mean gap of 0.612 cM was generated by utilizing genotyping by sequencing (GBS) technology. Every year for three years, 46 stable quantitative trait loci (QTLs) for leaf size and flavonoid contents were discovered, including 31 stable loci for Epimedin C (EC), one stable locus for total flavone content (TFC), 12 stable loci for leaf length (LL), and two stable loci for leaf area (LA). For flavonoid content and leaf size, the phenotypic variance explained for these loci varied between 4.00 and 16.80% and 14.95 and 17.34%, respectively. CONCLUSIONS: Forty-six stable QTLs for leaf size and flavonoid content traits were repeatedly detected over three years. The HDGM and stable QTLs are laying the basis for breeding and gene investigation in Epimedium and will contribute to accelerating the identification of desirable genotypes for Epimedium breeding.

Novel players in organogenesis and flavonoid biosynthesis in cucumber glandular trichomes.

Glandular trichomes (GTs) are outgrowths of plant epidermal cells that secrete and store specialized secondary metabolites that protect plants against biotic and abiotic stresses and have economic importance for human use. While extensive work has been done to understand the molecular mechanisms of trichome organogenesis in Arabidopsis (Arabidopsis thaliana), which forms unicellular, nonglandular trichomes (NGTs), little is known about the mechanisms of GT development or regulation of secondary metabolites in plants with multicellular GTs. Here, we identified and functionally characterized genes associated with GT organogenesis and secondary metabolism in GTs of cucumber (Cucumis sativus). We developed a method for effective separation and isolation of cucumber GTs and NGTs. Transcriptomic and metabolomic analyses showed that flavonoid accumulation in cucumber GTs is positively associated with increased expression of related biosynthesis genes. We identified 67 GT development-related genes, the functions of 7 of which were validated by virus-induced gene silencing. We further validated the role of cucumber ECERIFERUM1 (CsCER1) in GT organogenesis by overexpression and RNA interference transgenic approaches. We further show that the transcription factor TINY BRANCHED HAIR (CsTBH) serves as a central regulator of flavonoid biosynthesis in cucumber GTs. Work from this study provides insight into the development of secondary metabolite biosynthesis in multicellular GTs.

Cloning and functional analysis of GhDFR1, a key gene of flavonoid synthesis pathway in naturally colored cotton.

BACKGROUND: The naturally colored brown cotton fiber is the most widely used environmentally friendly textile material, which primarily contains proanthocyanidins and their derivatives. Many structural genes in the flavonoid synthesis pathway are known to improve the genetic resources of naturally colored cotton. Among them, DFR is a crucial late enzyme to synthesis both anthocyanins and proanthocyanidins in the plant flavonoid pathway. METHODS: The protein sequences of GhDFRs were analyzed using bioinformatic tools. The expression levels of GhDFRs in various tissues and organs of upland cotton Zongxu1 (ZX1), were analyzed by quantitative real-time PCR, and the expression pattern of GhDFR1 during fiber development of white cotton and brown cotton was analyzed further. The function of GhDFR1 in NCC ZX1 was preliminarily analyzed by virus induced gene silencing (VIGS) technology. RESULTS: Bioinformatic analysis revealed that GhDFRs sequences in upland cotton genome were extremely conserved. Furthermore, evolutionary tree analysis revealed that the functions of GhDFR1 and GhDFR2, and GhDFR3 and GhDFR4, presented different and shared some similarities. Our study showed GhDFR1 and GhDFR2 were specifically expressed in fibers, while GhDFR3 and GhDFR4 were specifically expressed in petals. GhDFR1 was exclusively expressed in brown cotton fiber at various stages of development and progressively increased with the growth of fiber, but the trend of expression in white cotton was quite the opposite. We silenced GhDFR1 expression in brown cotton fiber using VIGS technology, and observed the VIGS-interference plants. After reducing the expression level of GhDFR1, the period for significant GhDFR1 expression in the developing fibers changed, reducing the content of anthocyanins, and lightening the color of mature cotton fibers. CONCLUSION: GhDFR1 was preferentially expressed in brown cotton during fiber development. The timing of GhDFR1 expression for flavonoid synthesis altered, resulting in anthocyanin contents reduced and the fiber color of the GhDFR1i lines lightened. These findings showed the role of GhDFR1 in fiber coloration of NCC and provided a new candidate for NCC genetic improvement.