Two rice MYB transcription factors maintain male fertility in response to photoperiod by modulating sugar partitioning

OsMYB14, an R2R3-MYB transcription factor, regulates plant height through the control of hormone metabolism in rice

Plant growth must be regulated throughout the plant life cycle. The myeloblastosis (MYB) transcription factor (TF) family is one of the largest TF families and is involved in metabolism, lignin biosynthesis, and developmental processes. Here, we showed that OsMYB14, a rice R2R3-MYB TF, was expressed in leaves and roots, especially in rice culm and panicles, and that it localized to the nucleus. Overexpression of OsMYB14 (OsMYB14-ox) in rice resulted in a 30% reduction in plant height compared to that of the wild type (WT), while the height of the osmyb14-knockout (osmyb14-ko) mutant generated using the CRISPR/Cas9 system was not significantly different. Microscopic observations of the first internode revealed that the cell size did not differ significantly among the lines. RNA sequencing analysis revealed that genes associated with plant development, regulation, lipid metabolism, carbohydrate metabolism, and gibberellin (GA) and auxin metabolic processes were downregulated in the OsMYB14-ox line. Hormone quantitation revealed that inactive GA19 accumulated in OsMYB14-ox but not in the WT or knockout plants, suggesting that GA20 generation was repressed. Indole-3-acetic acid (IAA) and IAA-aspartate accumulated in OsMYB14-ox and osmyb14-ko, respectively. Indeed, real-time PCR analysis revealed that the expression of OsGA20ox1, encoding GA20 oxidase 1, and OsGH3-2, encoding IAA-amido synthetase, was downregulated in OsMYB14-ox and upregulated in osmyb14-ko. A protein-binding microarray revealed the presence of a consensus DNA-binding sequence, the ACCTACC-like motif, in the promoters of the OsGA20ox1 and GA20ox2 genes. These results suggest that OsMYB14 may act as a negative regulator of biological processes affecting plant height in rice by regulating GA biosynthesis and auxin metabolism.

The miR159a-DUO1 module regulates pollen development by modulating auxin biosynthesis and starch metabolism in citrus

Achieving seedlessness in citrus varieties is one of the important objectives of citrus breeding. Male sterility associated with abnormal pollen development is an important factor in seedlessness. However, our understanding of the regulatory mechanism underlying the seedlessness phenotype in citrus is still limited. Here, we determined that the miR159a-DUO1 module played an important role in regulating pollen development in citrus, which further indirectly modulated seed development and fruit size. Both the overexpression of csi-miR159a and the knocking out of DUO1 in Hong Kong kumquat (Fortunella hindsii) resulted in small and seedless fruit phenotypes. Moreover, pollen was severely aborted in both transgenic lines, with arrested pollen mitotic I and abnormal pollen starch metabolism. Through additional cross-pollination experiments, DUO1 was proven to be the key target gene for miR159a to regulate male sterility in citrus. Based on DNA affinity purification sequencing (DAP-seq), RNA-seq, and verified interaction assays, YUC2/YUC6, SS4 and STP8 were identified as downstream target genes of DUO1, those were all positively regulated by DUO1. In transgenic F. hindsii lines, the miR159a-DUO1 module down-regulated the expression of YUC2/YUC6, which decreased indoleacetic acid (IAA) levels and modulated auxin signaling to repress pollen mitotic I. The miR159a-DUO1 module reduced the expression of the starch synthesis gene SS4 and sugar transport gene STP8 to disrupt starch metabolism in pollen. Overall, this work reveals a new mechanism by which the miR159a-DUO1 module regulates pollen development and elucidates the molecular regulatory network underlying male sterility in citrus.

OsMADS6-OsMADS32 and REP1 control palea cellular heterogeneity and morphogenesis in rice

Precise regulation of cell proliferation and differentiation is vital for organ morphology. Rice palea, serving as sepal, comprises two distinct regions: the marginal region (MRP) and body of palea (BOP), housing heterogeneous cell populations, which makes it an ideal system for studying organ morphogenesis. We report that the transcription factor (TF) REP1 promotes epidermal cell proliferation and differentiation in the BOP, resulting in hard silicified protrusion cells, by regulating the cyclin-dependent kinase gene, OsCDKB1;1. Conversely, TFs OsMADS6 and OsMADS32 are expressed exclusively in the MRP, where they limit cell division rates by inhibiting OsCDKB2;1 expression and promote endoreduplication, yielding elongated epidermal cells. Furthermore, reciprocal inhibition between the OsMADS6-OsMADS32 complex and REP1 fine-tunes the balance between cell division and differentiation during palea morphogenesis. We further show the functional conservation of these organ identity genes in heterogeneous cell growth in Arabidopsis, emphasizing a critical framework for controlling cellular heterogeneity in organ morphogenesis.

mRNA cleavage by 21-nucleotide phasiRNAs determines temperature-sensitive male sterility in rice

Temperature-sensitive male sterility is one of the core components for hybrid rice (Oryza sativa) breeding based on the 2-line system. We previously found that knockout of ARGONAUTE 1d (AGO1d) causes temperature-sensitive male sterility in rice by influencing phased small interfering RNA (phasiRNA) biogenesis and function. However, the specific phasiRNAs and their targets underlying the temperature-sensitive male sterility in the ago1d mutant remain unknown. Here, we demonstrate that the ago1d mutant displays normal female fertility but complete male sterility at low temperature. Through a multiomics analysis of small RNA (sRNA), degradome, and transcriptome, we found that 21-nt phasiRNAs account for the greatest proportion of the 21-nt sRNA species in rice anthers and are sensitive to low temperature and markedly downregulated in the ago1d mutant. Moreover, we found that 21-nt phasiRNAs are essential for the mRNA cleavage of a set of fertility- and cold tolerance-associated genes, such as Earlier Degraded Tapetum 1 (EDT1), Tapetum Degeneration Retardation (TDR), OsPCF5, and OsTCP21, directly or indirectly determined by AGO1d-mediated gene silencing. The loss of function of 21-nt phasiRNAs can result in upregulation of their targets and causes varying degrees of defects in male fertility and grain setting. Our results highlight the essential functions of 21-nt phasiRNAs in temperature-sensitive male sterility in rice and suggest their promising application in 2-line hybrid rice breeding in the future.

OsUGE2 Regulates Plant Growth through Affecting ROS Homeostasis and Iron Level in Rice

BACKGROUND

The growth and development of rice (Oryza sativa L.) are affected by multiple factors, such as ROS homeostasis and utilization of iron. Here, we demonstrate that OsUGE2, a gene encoding a UDP-glucose 4-epimerase, controls growth and development by regulating reactive oxygen species (ROS) and iron (Fe) level in rice. Knockout of this gene resulted in impaired growth, such as dwarf phenotype, weakened root growth and pale yellow leaves. Biochemical analysis showed that loss of function of OsUGE2 significantly altered the proportion and content of UDP-Glucose (UDP-Glc) and UDP-Galactose (UDP-Gal). Cellular observation indicates that the impaired growth may result from decreased cell length. More importantly, RNA-sequencing analysis showed that knockout of OsUGE2 significantly influenced the expression of genes related to oxidoreductase process and iron ion homeostasis. Consistently, the content of ROS and Fe are significantly decreased in OsUGE2 knockout mutant. Furthermore, knockout mutants of OsUGE2 are insensitive to both Fe deficiency and hydrogen peroxide (H2O2) treatment, which further confirmed that OsUGE2 control rice growth possibly through Fe and H2O2 signal. Collectively, these results reveal a new pathway that OsUGE2 could affect growth and development via influencing ROS homeostasis and Fe level in rice.

MADS8 is indispensable for female reproductive development at high ambient temperatures in cereal crops

Temperature is a major factor that regulates plant growth and phenotypic diversity. To ensure reproductive success at a range of temperatures, plants must maintain developmental stability of their sexual organs when exposed to temperature fluctuations. However, the mechanisms integrating plant floral organ development and temperature responses are largely unknown. Here, we generated barley and rice loss-of-function mutants in the SEPALLATA-like MADS-box gene MADS8. The mutants in both species form multiple carpels that lack ovules at high ambient temperatures. Tissue-specific markers revealed that HvMADS8 is required to maintain floral meristem determinacy and ovule initiation at high temperatures, and transcriptome analyses confirmed that temperature-dependent differentially expressed genes in Hvmads8 mutants predominantly associate with floral organ and meristem regulation. HvMADS8 temperature-responsive activity relies on increased binding to promoters of downstream targets, as revealed by a cleavage under targets and tagmentation (CUT&Tag) analysis. We also demonstrate that HvMADS8 directly binds to 2 orthologs of D-class floral homeotic genes to activate their expression. Overall, our findings revealed a new, conserved role for MADS8 in maintaining pistil number and ovule initiation in cereal crops, extending the known function of plant MADS-box proteins in floral organ regulation.

MYB transcription factors and their roles in the male reproductive development of flowering plants

As one of the largest transcription factor families with complex functional differentiation in plants, the MYB transcription factors (MYB TFs) play important roles in the physiological and biochemical processes of plant growth and development. Male reproductive development, an essential part of sexual reproduction in flowering plants, is undoubtedly regulated by MYB TFs. In this review, we summarize the roles of the MYB TFs involved in the three stages of male reproductive development: pollen grains formation and maturation, filament elongation and anther dehiscence, and fertilization. Also, the potential downstream target genes and upstream regulators of these MYB TFs are discussed. Furthermore, we propose the underlying regulatory mechanisms of these MYB TFs: (1) A complex network of MYB TFs regulates various aspects of male reproductive development; (2) MYB homologous genes in different species may be functionally conserved or differentiated; (3) MYB TFs often form regulatory complexes with bHLH TFs.

Transcriptome analysis provides insight into the regulatory mechanisms underlying pollen germination recovery at normal high ambient temperature in wild banana (Musa itinerans)

Introduction

Cultivated banana are polyploid, with low pollen fertility, and most cultivars are male sterile, which leads to difficulties in banana breeding research. The selection of male parent with excellent resistance and pollen fertility is therefore essential for banana breeding. Wild banana (Musa itinerans) have developed many good characteristics during natural selection and constitute an excellent gene pool for breeding. Therefore, research on wild banana breeding is very important for banana breeding.

Results

In the current analysis, we examined the changes in viability of wild banana pollens at different temperatures by in vitro germination, and found that the germination ability of wild banana pollens cultured at 28°C for 2 days was higher than that of pollens cultured at 23°C (pollens that could not germinate normally under low temperature stress), 24°C (cultured at a constant temperature for 2 days) and 32°C (cultured at a constant temperature for 2 days). To elucidate the molecular mechanisms underlying the germination restoration process in wild banana pollens, we selected the wild banana pollens that had lost its germination ability under low temperature stress (23°C) as the control group (CK) and the wild banana pollens that had recovered its germination ability under constant temperature incubation of 28°C for 2 days as the treatment group (T) for transcriptome sequencing. A total of 921 differentially expressed genes (DEGs) were detected in CK vs T, of which 265 were up-regulated and 656 were down-regulated. The combined analysis of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) revealed that the activation, metabolism of various substances (lipids, sugars, amino acids) play a major role in restoring pollen germination capacity. TCA cycle and the sesquiterpenoid and triterpenoid biosynthetic pathways were also significantly enriched in the KEGG pathway. And we found that some DEGs may be associated with pollen wall formation, DNA methylation and DNA repair. The cysteine content, free fatty acid (FFA) content, H2O2 content, fructose content, and sucrose content of pollen were increased at treatment of 28°C, while D-Golactose content was decreased. Finally, the GO pathway was enriched for a total of 24 DEGs related to pollen germination, of which 16 DEGs received targeted regulation by 14 MYBs.

Discussions

Our study suggests that the balance between various metabolic processes, pollen wall remodelling, DNA methylation, DNA repairs and regulation of MYBs are essential for germination of wild banana pollens.

Transcriptome Co-Expression Network Analysis of Peach Fruit with Different Sugar Concentrations Reveals Key Regulators in Sugar Metabolism Involved in Cold Tolerance

Peach fruits are known to be highly susceptible to chilling injury (CI) during low-temperature storage, which has been linked to the level of sugar concentration in the fruit. In order to better understand the relationship between sugar metabolism and CI, we conducted a study examining the concentration of sucrose, fructose, and glucose in peach fruit with different sugar concentrations and examined their relationship with CI. Through transcriptome sequencing, we screened the functional genes and transcription factors (TFs) involved in the sugar metabolism pathway that may cause CI in peach fruit. Our results identified five key functional genes (PpSS, PpINV, PpMGAM, PpFRK, and PpHXK) and eight TFs (PpMYB1/3, PpMYB-related1, PpWRKY4, PpbZIP1/2/3, and PpbHLH2) that are associated with sugar metabolism and CI development. The analysis of co-expression network mapping and binding site prediction identified the most likely associations between these TFs and functional genes. This study provides insights into the metabolic and molecular mechanisms regulating sugar changes in peach fruit with different sugar concentrations and presents potential targets for breeding high-sugar and cold-tolerant peach varieties.

Over-expression of CcMYB24, encoding a R2R3-MYB transcription factor from a high-leaf-number mutant of Cymbidium, increases the number of leaves in Arabidopsis

Ornamental foliage plants have long been cultivated for their attractive leaves. Variation in leaf traits of ornamental foliage plants is one of the goals in breeding. MYB transcription factors regulate many aspects of leaf development, and thus influence morphological traits of leaves. However, little is known about the function of MYB transcription factors in leaf development of Cymbidium, one of the most economically important ornamental plants in the world. In the present study, a MYB transcription factor, CcMYB24, was identified and the corresponding gene cloned from a new orchid mutant, TRIR-2, which produces more leaves than control plants. The CcMYB24 showed a higher expression level in 'TRIR-2' than in control plants, and the protein was located in the nucleus. The sequence of CcMYB24 showed a high similarity with RAX2-like genes which belong to the R2R3-MYB gene family in other Cymbidium plants. Overexpression of CcMYB24 resulted in a phenotype with an increased number of leaves, elevated chlorophyll content, and decreased contents of carotenoids and flavonoids in Arabidopsis. These results provide functional evidence for the role of CcMYB24 in promoting the production of leaves in 'TRIR-2'. Understanding the role of CcMYB24 in Cymbidium will be beneficial for the molecular breeding of ornamental foliage plants.

Environment-sensitive genic male sterility in rice and other plants

Environment-sensitive genic male sterility is a type of male sterility that is affected by both genetic and environmental factors. Environment-sensitive genic male sterile lines are not only used in two-line hybrid breeding but are also good materials for studying plant-environment interactions. In this study we review the research progress on environment-sensitive genic male sterility in rice from the perspectives of epigenetic, transcriptional, posttranscriptional, posttranslational and metabolic mechanisms as well as signal transduction processes. While significant progress has been made in the genetics, gene cloning and understanding of the molecular mechanisms of environment-sensitive genic male sterility in recent years, the relevant regulatory network is still poorly understood in rice. We therefore also review studies of environment-sensitive genic male sterility in Arabidopsis and other crops, hoping to promote research in this field in rice. Finally, we analyse the challenges posed by environment-sensitive genic male sterility and provide corresponding suggestions. This review will contribute towards an understanding the molecular genetics of environment-sensitive genic male sterility and its application in two-line hybrid breeding in rice and other species.

Knocking Out the Transcription Factor OsNAC092 Promoted Rice Drought Tolerance

Environmental drought stress threatens rice production. Previous studies have reported that related NAC (NAM, ATAF1/2, and CUC) transcription factors play an important role in drought stress. Herein, we identified and characterized OsNAC092, encoding an NAC transcription factor that is highly expressed and induced during drought tolerance. OsNAC092 knockout lines created using the clustered regularly interspaced palindromic repeats (CRISPR)-associated protein 9 (Cas9) system exhibited increased drought resistance in rice. RNA sequencing showed that the knockout of OsNAC092 caused a global expression change, and differential gene expression is chiefly associated with "response to light stimulus," "MAPK signaling pathway," "plant hormone signal transduction," "response to oxidative stress," "photosynthesis," and "water deprivation." In addition, the antioxidants and enzyme activities of the redox response were significantly increased. OsNAC092 mutant rice exhibited a higher ability to scavenge more ROS and maintained a high GSH/GSSG ratio and redox level under drought stress, which could protect cells from oxidant stress, revealing the importance of OsNAC092 in the rice's response to abiotic stress. Functional analysis of OsNAC092 will be useful to explore many rice resistance genes in molecular breeding to aid in the development of modern agriculture.

Source-To-Sink Transport of Sugar and Its Role in Male Reproductive Development

Sucrose is produced in leaf mesophyll cells via photosynthesis and exported to non-photosynthetic sink tissues through the phloem. The molecular basis of source-to-sink long-distance transport in cereal crop plants is of importance due to its direct influence on grain yield-pollen grains, essential for male fertility, are filled with sugary starch, and rely on long-distance sugar transport from source leaves. Here, we overview sugar partitioning via phloem transport in rice, especially where relevant for male reproductive development. Phloem loading and unloading in source leaves and sink tissues uses a combination of the symplastic, apoplastic, and/or polymer trapping pathways. The symplastic and polymer trapping pathways are passive processes, correlated with source activity and sugar gradients. In contrast, apoplastic phloem loading/unloading involves active processes and several proteins, including SUcrose Transporters (SUTs), Sugars Will Eventually be Exported Transporters (SWEETs), Invertases (INVs), and MonoSaccharide Transporters (MSTs). Numerous transcription factors combine to create a complex network, such as DNA binding with One Finger 11 (DOF11), Carbon Starved Anther (CSA), and CSA2, which regulates sugar metabolism in normal male reproductive development and in response to changes in environmental signals, such as photoperiod.

Genome-Wide Identification of R2R3-MYB Transcription Factor and Expression Analysis under Abiotic Stress in Rice

The myeloblastosis (MYB) family comprises a large group of transcription factors (TFs) that has a variety of functions. Among them, the R2R3-MYB type of proteins are the largest group in plants, which are involved in controlling various biological processes such as plant growth and development, physiological metabolism, defense, and responses to abiotic and biotic stresses. In this study, bioinformatics was adopted to conduct genome-wide identification of the R2R3-MYB TFs in rice. We identified 190 MYB TFs (99 R2R3-MYBs), which are unevenly distributed on the 12 chromosomes of rice. Based on the phylogenetic clustering and protein sequence characteristics, OsMYBs were classified into five subgroups, and 59.6% of the Os2R_MYB genes contained two introns. Analysis of cis-acting elements in the 2000 bp upstream region of Os2R_MYB genes showed that all Os2R_MYB genes contained plant hormones-related or stress-responsive elements since 91.9%, 79.8%, 79.8%, and 58.6% of Os2R_MYB genes contain ABRE, TGACG, CGTCA, and MBS motifs, respectively. Protein–protein network analysis showed that the Os2R_MYBs were involved in metabolic process, biosynthetic process, and tissue development. In addition, some genes showed a tissue-specific or developmental-stage-specific expression pattern. Moreover, the transcription levels of 20 Os2R_MYB genes under polyethylene glycol (PEG) and cadmium chloride (CdCl₂) stress inducers were dissected by qRT-PCR. The results indicated genes with an altered expression upon PEG or CdCl₂ stress induction. These results potentially supply a basis for further research on the role that Os2R_MYB genes play in plant development and stress responses.

The regulatory role of CARBON STARVED ANTHER-mediated photoperiod-dependent male fertility in rice

Environmental signals, especially daylength, play important roles in determining fertility in photoperiod-sensitive genic male sterile (PGMS) lines that are critical to sustain production of high-yielding hybrid rice (Oryza sativa) varieties. However, the mechanisms by which PGMS lines perceive changes in photoperiod and transmit those signals to elicit downstream effects are not well understood. In this study, we compared the transcriptomes from the leaves and anthers of carbon starved anther (csa), a PGMS line, to wild-type (WT) tissues under different photoperiods. Components of circadian clock in the leaves, including Circadian Clock-Associated 1 and Pseudo-Response Regulator (PRR95), played vital roles in sensing the photoperiod signals. Photoperiod signals were weakly transduced to anthers, where gene expression was mainly controlled by the CSA allele. CSA played a critical role in regulating sugar metabolism and cell wall synthesis in anthers under short-day conditions, and transcription of key genes inducing csa-directed sterility was upregulated under long-day (LD) conditions though not to WT levels, revealing a mechanism to explain the partial restoration of fertility in rice under LD conditions. Eight direct targets of CSA regulation were identified, all of which were genes involved in sugar metabolism and transport (cell wall invertases, SWEETs, and monosaccharide transporters) expressed only in reproductive tissues. Several hub genes coordinating the effects of CSA regulation were identified as critical elements determining WT male fertility and further analysis of these and related genes will reveal insights into how CSA coordinates sugar metabolism, cell wall biosynthesis, and photoperiod sensing in rice anther development.

VvMYB1 potentially affects VvTOR gene expression by regulating VvTOR promoter and participates in glucose accumulation

MYB (v-myb avian myeloblastosis viral oncogene homolog) transcription factors make up one of the largest protein families in plants. The TOR (target of rapamycin) signaling network plays a pivotal role in sugar metabolism and plant growth. In this article, we utilized grape (Vitis vinifera) calli to explore the relationship between VvMYB1 and VvTOR. By using yeast one-hybrid and dual-luciferase reporter system, we speculated that there may be other proteins that help VvMYB1 and VvTOR promoter bond in grape calli, and the interaction action sites were located between the VvTOR 400-bp promoter fragment and the 1200-bp promoter fragment. The subcellular localization results suggest that VvMYB1 is found in the nucleus. Moreover, the expression level of VvTOR increased in the transgenic calli with overexpression of VvMYB1. These findings provide further evidence that VvMYB1 regulates VvTOR expression. We also found that overexpression of VvMYB1 increased glucose accumulation and affected expression of sugar-related genes. Our results suggest that there is a crosstalk between VvMYB1, VvTOR, and glucose accumulation.

Transcriptomics-Based Identification of Genes Related to Tapetum Degradation and Microspore Development in Lily

Lily is a popular and economically ornamental crop around the world. However, its high production of pollen grains causes serious problems to consumers, including allergies and staining of clothes. During anther development, the tapetum is a crucial step for pollen formation and microspore release. Therefore, it is important to understand the mechanism of tapetum degradation and microspore development in lily where free pollen contamination occurs. Here, we used the cut lily cultivar ‘Siberia’ to characterize the process of tapetum degradation through the use of cytology and transcriptomic methods. The cytological observation indicated that, as the lily buds developed from 4 cm (Lo 4 cm) to 8 cm (Lo 8 cm), the tapetum completed the degradation process and the microspores matured. Furthermore, by comparing the transcriptome profiling among three developmental stages (Lo 4 cm, Lo 6 cm and Lo 8 cm), we identified 27 differentially expressed genes. These 27 genes were classed into 4 groups by function, namely, cell division and expansion, cell-wall morphogenesis, transcription factors, LRR-RLK (leucine-rich repeat receptor-like kinases), plant hormone biosynthesis and transduction. Quantitative real-time PCR was performed as validation of the transcriptome data. These selected genes are candidate genes for the tapetum degradation and microspore development of lily and our work provides a theoretical basis for breeding new lily cultivars without pollen.

Genome wide characterization of R2R3 MYB transcription factor from Apocynum venetum revealed potential stress tolerance and flavonoid biosynthesis genes

MYB transcription factors are crucial in regulating stress tolerance and expression of major genes involved in flavonoid biosynthesis. The functions of MYBs is well explored in a number of plants, yet no studies is reported in Apocynum venetum. We identified a total of 163 MYB candidates, that comprised of 101 (61.96%) R2R3, 6 3R, 1 4R and 55 1R. Syntenic analysis of A. venetum R2R3 (AvMYB) showed highest orthologous pairs with Vitis vinifera MYBs followed by Arabidopsis thaliana among the four species evaluated. Thirty segmental duplications and 6 tandem duplications were obtained among AvMYB gene pairs signifying their role in the MYB gene family expansion. Nucleotide substitution analysis (Ka/Ks) showed the AvMYBs to be under the influence of strong purifying selection. Expression analysis of selected AvMYB under low temperature and cadmium stresses resulted in the identification of AvMYB48, AvMYB97, AvMYB8,AvMYB4 as potential stress responsive genes and AvMYB10 and AvMYB11 in addition, proanthocyanidin biosynthesis regulatory genes which is consistent with their annotated homologues in Arabidopsis. Tissue specific expression profile analysis of AvMYBs further supported the qPCR analysis result. MYBs with higher transcript levels in root, stem and leaf like AvMYB4 forexample, was downregulated under the stresses and such with low transcript level such as AvMYB48 which had low transcript in the leaf was upregulated under both stresses. Transcriptome and phylogenetic analysis suggested AvMYB42 as a potential regulator of anthocyanin biosynthesis. Thus, this study provided valuable information on AvR2R3-MYB gene family with respect to stress tolerance and flavonoid biosynthesis.