A R2R3-MYB transcription factor, GmMYB12B2, affects the expression levels of flavonoid biosynthesis genes encoding key enzymes in transgenic Arabidopsis plants

Deciphering the regulatory networks involved in mild and severe salt stress responses in the roots of wild grapevine Vitis vinifera spp. sylvestris

Transcriptional regulatory networks are pivotal components of plant's response to salt stress. However, plant adaptation strategies varied as a function of stress intensity, which is mainly modulated by climate change. Here, we determined the gene regulatory networks based on transcription factor (TF) TF_gene co-expression, using two transcriptomic data sets generated from the salt-tolerant "Tebaba" roots either treated with 50 mM NaCl (mild stress) or 150 mM NaCl (severe stress). The analysis of these regulatory networks identified specific TFs as key regulatory hubs as evidenced by their multiple interactions with different target genes related to stress response. Indeed, under mild stress, NAC and bHLH TFs were identified as central hubs regulating nitrogen storage process. Moreover, HSF TFs were revealed as a regulatory hub regulating various aspects of cellular metabolism including flavonoid biosynthesis, protein processing, phenylpropanoid metabolism, galactose metabolism, and heat shock proteins. These processes are essentially linked to short-term acclimatization under mild salt stress. This was further consolidated by the protein-protein interaction (PPI) network analysis showing structural and plant growth adjustment. Conversely, under severe salt stress, dramatic metabolic changes were observed leading to novel TF members including MYB family as regulatory hubs controlling isoflavonoid biosynthesis, oxidative stress response, abscisic acid signaling pathway, and proteolysis. The PPI network analysis also revealed deeper stress defense changes aiming to restore plant metabolic homeostasis when facing severe salt stress. Overall, both the gene co-expression and PPI network provided valuable insights on key transcription factor hubs that can be employed as candidates for future genetic crop engineering programs.

Integration of genome-wide association studies, metabolomics, and transcriptomics reveals phenolic acid- and flavonoid-associated genes and their regulatory elements under drought stress in rapeseed flowers

Introduction

Biochemical and metabolic processes help plants tolerate the adverse effects of drought. In plants accumulating bioactive compounds, understanding the genetic control of the biosynthesis of biochemical pathways helps the discovery of candidate gene (CG)-metabolite relationships.

Methods

The metabolic profile of flowers in 119 rapeseed (Brassica napus) accessions was assessed over two irrigation treatments, one a well-watered (WW) condition and the other a drought stress (DS) regime. We integrated information gained from 52,157 single-nucleotide polymorphism (SNP) markers, metabolites, and transcriptomes to identify linked SNPs and CGs responsible for the genetic control of flower phenolic compounds and regulatory elements.

Results

In a genome-wide association study (GWAS), of the SNPs tested, 29,310 SNPs were qualified to assess the population structure and linkage disequilibrium (LD), of which several SNPs for radical scavenging activity (RSA) and total flavanol content (TFLC) were common between the two irrigation conditions and pleiotropic SNPs were found for chlorogenic and coumaric acids content. The principal component analysis (PCA) and stepwise regression showed that chlorogenic acid and epicatechin in WW and myricetin in DS conditions were the most important components for RSA. The hierarchical cluster analysis (HCA) showed that vanillic acid, myricetin, gallic acid, and catechin were closely associated in both irrigation conditions. Analysis of GWAS showed that 60 CGs were identified, of which 18 were involved in stress-induced pathways, phenylpropanoid pathway, and flavonoid modifications. Of the CGs, PAL1, CHI, UGT89B1, FLS3, CCR1, and CYP75B137 contributed to flavonoid biosynthetic pathways. The results of RNA sequencing (RNA-seq) revealed that the transcript levels of PAL, CHI, and CYP75B137 known as early flavonoid biosynthesis-related genes and FLS3, CCR1, and UGT89B1 related to the later stages were increased during drought conditions. The transcription factors (TFs) NAC035 and ERF119 related to flavonoids and phenolic acids were upregulated under drought conditions.

Discussion

These findings expand our knowledge on the response mechanisms to DS, particularly regarding the regulation of key phenolic biosynthetic genes in rapeseed. Our data also provided specific linked SNPs for marker-assisted selection (MAS) programs and CGs as resources toward realizing metabolomics-associated breeding of rapeseed.

Integration of transcriptomic and metabolomic analyses provides insights into response mechanisms to nitrogen and phosphorus deficiencies in soybean

Nitrogen (N) and phosphorus (P) are two essential plant macronutrients that can limit plant growth by different mechanisms. We aimed to shed light on how soybean respond to low nitrogen (LN), low phosphorus (LP) and their combined deficiency (LNP). Generally, these conditions triggered changes in gene expression of the same processes, including cell wall organization, defense response, response to oxidative stress, and photosynthesis, however, response was different in each condition. A typical primary response to LN and LP was detected also in soybean, i.e., the enhanced uptake of N and P, respectively, by upregulation of genes for the corresponding transporters. The regulation of genes involved in cell wall organization showed that in LP roots tended to produce more casparian strip, in LN more secondary wall biosynthesis occurred, and in LNP reduction in expression of genes involved in secondary wall production accompanied by cell wall loosening was observed. Flavonoid biosynthesis also showed distinct pattern of regulation in different conditions: more anthocyanin production in LP, and more isoflavonoid production in LN and LNP, which we confirmed also on the metabolite level. Interestingly, in soybean the nutrient deficiencies reduced defense response by lowering expression of genes involved in defense response, suggesting a role of N and P nutrition in plant disease resistance. In conclusion, we provide detailed information on how LN, LP, and LNP affect different processes in soybean roots on the molecular and physiological levels.

GsMYB7 encoding a R2R3-type MYB transcription factor enhances the tolerance to aluminum stress in soybean (Glycine max L.)

Background

MYB transcription factor (TF) is one of the largest families of TFs in plants and play essential roles in plant growth and development, and is involved in responses to biological and abiotic stress. However, there are few reports on GsMYB7 gene in soybean under aluminum acid stress, and its regulatory mechanism remains unclear.

Results

The GsMYB7 protein is localized in the nucleus and has transcriptional activation ability. Quantitative real-time PCR (qRT-PCR) results showed that GsMYB7 held a constitutive expression pattern rich in roots. When AlCl₃ concentration was 25 µM, the total root surface area (SA) of GsMYB7 transgenic lines were 34.97% higher than that of wild-type Huachun 6 (HC6). While the accumulation of Al³⁺ in root tip of transgenic plants after aluminum treatment was 17.39% lower than that of wild-type. RNA-sequencing analysis indicated that over 1181 genes were regulated by GsMYB7 and aluminum stress. Among all the regulated genes, the expression levels of glutathione peroxidase, protein kinase, cytochrome and other genes in the transgenic lines were significantly higher than those in wild type by acidic aluminum stress. The bioinformatics and qRT-PCR results showed that 9 candidate genes were induced under the treatments of acidic aluminum stress which were indirectly and/or directly regulated by GsMYB7. After AlCl₃ treatments, the transcripts of these genes in GsMYB7 transgenic seedlings were significantly higher than those of wide-type HC6.

Conclusions

The results suggested that GsMYB7 may enhance soybean tolerance to acidic aluminum stress by regulating the downstream genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08744-w.

A comparative study on photosynthetic characteristics and flavonoid metabolism between Camellia petelotii (Merr.) Sealy and Camellia impressinervis Chang &Liang

Camellia petelotii (Merr.) Sealy and Camellia impressinervis Chang & Liang belong to the golden subgroup of Camellia (Theaceae). This subgroup contains the yellow-flowering species of the genus, which have high medicinal and ornamental value and a narrow geographical distribution. These species differ in their tolerance to high light intensity. This study aimed to explore the differences in their light-stress responses and light damage repair processes, and the effect of these networks on secondary metabolite synthesis. Two-year-old plants of both species grown at 300 µmol·m^-2·s^-1 photosynthetically active radiation (PAR) were shifted to 700 µmol·m^-2·s^-1 PAR for 5 days shifting back to 300 µmol·m^-2·s^-1 PAR for recovery for 5 days. Leaf samples were collected at the start of the experiment and 2 days after each shift. Data analysis included measuring photosynthetic indicators, differential transcriptome expression, and quantifying plant hormones, pigments, and flavonoids. Camellia impressinervis showed a weak ability to recover from photodamage that occurred at 700 µmol·m^-2·s^-1 compared with C. petelotii. Photodamage led to decreased photosynthesis, as shown by repressed transcript abundance for photosystem II genes psbA, B, C, O, and Q, photosystem I genes psaB, D, E, H, and N, electron transfer genes petE and F, and ATP synthesis genes ATPF1A and ATPF1B. High-light stress caused more severe damage to C. impressinervis, which showed a stronger response to reactive oxygen species than C. petelotii. In addition, high-light stress promoted the growth and development of high zeatin signalling and increased transcript abundance of adenylate dimethylallyl transferase (IPT) and histidine-containing phosphotransferase (AHP). The identification of transcriptional differences in the regulatory networks that respond to high-light stress and activate recovery of light damage in these two rare species adds to the resources available to conserve them and improve their value through molecular breeding.

Transcriptome analysis of MYB transcription factors family and PgMYB genes involved in salt stress resistance in Panax ginseng

BACKGROUND

As the king of all herbs, the medicinal value of ginseng is self-evident. The perennial nature of ginseng causes its quality to be influenced by various factors, one of which is the soil environment. During plant growth and development, MYB transcription factors play an important role in responding to abiotic stresses and regulating the synthesis of secondary metabolites. However, there are relatively few reports on the MYB transcription factor family in Panax ginseng.

RESULTS

This study identified 420 PgMYB transcripts under 117 genes ID in the Jilin ginseng transcriptome database. Phylogenetic analysis showed that PgMYB transcripts in Jilin ginseng were classified into 19 functional subclasses. The GO annotation result indicated that the functional differentiation of PgMYB transcripts was annotated to 11 functional nodes at GO Level 2 in ginseng. Expression pattern analysis of PgMYB transcripts based on the expression data (TPM) that PgMYB transcripts were revealed spatiotemporally specific in expression patterns. We performed a weighted network co-expression network analysis on the expression of PgMYB transcripts from different samples. The co-expression network containing 51 PgMYB transcripts was formed under a soft threshold of 0.85, revealing the reciprocal relationship of PgMYB in ginseng. Treatment of adventitious roots of ginseng with different concentrations of NaCl revealed four up-regulated expression of PgMYB transcripts that can candidate genes for salt resistance studies in ginseng.

CONCLUSIONS

The present findings provide data resources for the subsequent study of the functions of MYB transcription factor family members in ginseng, and provide an experimental basis for the anti-salt functions of MYB transcription factors in Panax ginseng.

AgMYB12, a novel R2R3-MYB transcription factor, regulates apigenin biosynthesis by interacting with the AgFNS gene in celery

Overexpression of AgMYB12 in celery improved the accumulation of apigenin by interacting with the AgFNS gene. Celery is a common vegetable, and its essential characteristic is medicine food homology. A natural flavonoid and a major pharmacological component in celery, apigenin plays an important role in human health. In this study, we isolated a novel R2R3-MYB transcription factor that regulates apigenin accumulation from the celery cultivar 'Jinnan Shiqin' through yeast one-hybrid screening and designated it as AgMYB12. The AgMYB12 protein was located in the nucleus. It showed transcriptional activation activity and bound specifically to the promoter of AgFNS, a gene involved in apigenin biosynthesis. Phylogenetic tree analysis demonstrated that AgMYB12 belongs to the flavonoid branch. It contains two flavonoid-related motifs, SG7 and SG7-2, and shared a highly conserved R2R3 domain with flavonoid-related MYBs. The homologous overexpression of AgMYB12 induced the up-regulation of AgFNS gene expression and accumulation of apigenin and luteolin in celery. Additionally, the expression levels of apigenin biosynthesis-related genes, including AgPAL, AgCHI, AgCHS, Ag4CL, and AgC4H, increased in transgenic celery plants. These results indicated that AgMYB12 acted as a positive regulator of apigenin biosynthesis and activated the expression of AgFNS gene. The current study provides new information about the regulation mechanism of apigenin metabolism in celery and offers a strategy for cultivating the plants with high apigenin content.

Transcriptome and Metabolome Analysis Unveil Anthocyanin Metabolism in Pink and Red Testa of Peanut (Arachis hypogaea L.)

Peanut (Arachis hypogaea L.) is an important source of oil and food around the world, and the testa color affects its appearance and commercial value. However, few studies focused on the mechanism of pigment formation in peanut testa. In this study, cultivars Shanhua 15 with pink testa and Zhonghua 12 with red testa were used as materials to perform the combined analysis of transcriptome and metabolome. A total of 198 flavonoid metabolites were detected, among which petunidin 3-O-glucoside and cyanidin O-acetylhexoside in Zhonghua12 were 15.23 and 14.72 times higher than those of Shanhua 15 at the R7 stage, revealing the anthocyanins underlying the red testa. Transcriptome analysis showed that there were 6059 and 3153 differentially expressed genes between Shanhua 15 and Zhonghua 12 in different growth periods, respectively. These differentially expressed genes were significantly enriched in the flavonoid biosynthesis, biosynthesis of secondary metabolites, and metabolic pathways. Integrated analysis of transcriptome and metabolome indicated CHS gene (arahy.CM90T6), F3'H genes (arahy. 8F7PE4 and arahy. K8H9R8), and DFR genes (arahy. LDV9QN and arahy. X8EVF3) may be the key functional genes controlling the formation of pink and red testa in peanut. Transcription factors MYB (arahy.A2IWKV, arahy.US2SKM, arahy.SJGE27, arahy.H8DJRL, and arahy.PR7AYB), bHLH (arahy.26781N, arahy.HM1IVV, and arahy.MP3D3D), and WD40 (arahy.L6JJW9) in the biosynthetic pathway of anthocyanin were significantly upregulated in Zhonghua 12 which may be the key regulatory genes in testa pigment formation. This is a comprehensive analysis on flavonoid metabolites and related genes expression in peanut testa, providing reference for revealing the regulatory mechanism of pigment accumulation in peanut testa.

Full-Length Transcriptome Analysis Reveals Candidate Genes Involved in Terpenoid Biosynthesis in Artemisia argyi

Artemisia argyi is an important medicinal plant widely utilized for moxibustion heat therapy in China. The terpenoid biosynthesis process in A. argyi is speculated to play a key role in conferring its medicinal value. However, the molecular mechanism underlying terpenoid biosynthesis remains unclear, in part because the reference genome of A. argyi is unavailable. Moreover, the full-length transcriptome of A. argyi has not yet been sequenced. Therefore, in this study, de novo transcriptome sequencing of A. argyi's root, stem, and leaf tissues was performed to obtain those candidate genes related to terpenoid biosynthesis, by combining the PacBio single-molecule real-time (SMRT) and Illumina sequencing NGS platforms. And more than 55.4 Gb of sequencing data and 108,846 full-length reads (non-chimeric) were generated by the Illumina and PacBio platform, respectively. Then, 53,043 consensus isoforms were clustered and used to represent 36,820 non-redundant transcripts, of which 34,839 (94.62%) were annotated in public databases. In the comparison sets of leaves vs roots, and leaves vs stems, 13,850 (7,566 up-regulated, 6,284 down-regulated) and 9,502 (5,284 up-regulated, 4,218 down-regulated) differentially expressed transcripts (DETs) were obtained, respectively. Specifically, the expression profile and KEGG functional enrichment analysis of these DETs indicated that they were significantly enriched in the biosynthesis of amino acids, carotenoids, diterpenoids and flavonoids, as well as the metabolism processes of glycine, serine and threonine. Moreover, multiple genes encoding significant enzymes or transcription factors related to diterpenoid biosynthesis were highly expressed in the A. argyi leaves. Additionally, several transcription factor families, such as RLK-Pelle_LRR-L-1 and RLK-Pelle_DLSV, were also identified. In conclusion, this study offers a valuable resource for transcriptome information, and provides a functional genomic foundation for further research on molecular mechanisms underlying the medicinal use of A. argyi leaves.

SlMYB14 promotes flavonoids accumulation and confers higher tolerance to 2,4,6-trichlorophenol in tomato

Flavonoids are small molecular secondary metabolites, which have a variety of biological functions. Transcriptional regulations of key enzyme genes play critical roles in the flavonoid biosynthesis. In this study, an R2R3-MYB transcription factor gene, SlMYB14, was isolated from tomato and characterized. The nucleus-localized SlMYB14 functions as a transcriptional activator in yeast. The expression of SlMYB14 could be induced by methyl jasmonic acid, wounding and ABA. SlMYB14 works downstream of SlMYC2 in the jasmonate signaling pathway. Overexpression of SlMYB14 under the control of CaMV35S promoter in tomato led to increased accumulation of flavonoids. RNA-sequencing analysis revealed that the transcript levels of several structural genes associated with flavonoid biosynthesis were up-regulated in transgenic tomato plants. Gel-shift assays confirmed that SlMYB14 protein could bind to the promoter regions of SlPAL genes. It was also found that overexpression of SlMYB14 improved the tolerance of transgenic plants to 2,4,6-trichlorophenol (2,4,6-TCP), an environmental organic pollutant which could cause serious oxidative damage to plant. These results suggest that SlMYB14 participates in the regulation of flavonoid biosynthesis and might play a role in maintaining reactive oxygen species homeostasis in plant. SlMYB14 gene also has the potential to contribute to the phytoremediation of 2,4,6-TCP-contaminated soils.

Identification of Isoflavonoid Biosynthesis-Related R2R3-MYB Transcription Factors in Callerya speciosa (Champ. ex Benth.) Schot Using Transcriptome-Based Gene Coexpression Analysis

The R2R3-MYB family is one of the largest plant transcription factor (TF) families playing vital roles in defense, plant growth, and secondary metabolism biosynthesis. Although this gene family has been studied in many species, isoflavonoid biosynthesis-related R2R3-MYB TFs in Callerya speciosa (Champ. ex Benth.) Schot, a traditional Chinese medicinal herb, are poorly understood. Here, a total of 101 R2R3-MYB TFs were identified from C. speciosa transcriptome dataset. 25 clades divided into five functional groups were clustered based on the sequence similarity and phylogenetic tree. Conserved motifs and domain distribution, expression patterns, and coexpression networks were also employed to identify the potential R2R3-MYB TFs in the regulation of isoflavonoid biosynthesis. In silico evaluation showed that the deduced R2R3-CsMYB proteins contain highly conserved R2R3 repeat domain at the N-terminal region, that is the signature motif of R2R3-type MYB TFs. Eight potential TFs (CsMYB17, CsMYB36, CsMYB41, CsMYB44, CsMYB45, CsMYB46, CsMYB72, and CsMYB81) had high degrees of coexpression with four key isoflavonoid biosynthetic genes (CsIFS, CsCHS7, CsHID-1, and CsCHI3), in which CsMYB36 as a potential regulator possessed the highest degree. HPLC analysis showed that formononetin and maackiain contents were significantly increased during the development of tuberous roots, which might be controlled by both related R2R3-CsMYBs and structural genes involved in the isoflavonoid biosynthesis pathway. The transcriptome data were further validated by reverse transcription real-time PCR (RT-qPCR) analysis, and similar expression profiles between TFs and key structural genes were identified. This study was the first step toward the understanding of the R2R3-MYB TFs regulating isoflavonoid biosynthesis in C. speciosa. The results will provide information for further functional analysis and quality improvement through genetic manipulation of these potential R2R3-CsMYB genes in C. speciosa.

Integrated transcriptomic and metabolomic data reveal the flavonoid biosynthesis metabolic pathway in Perilla frutescens (L.) leaves

Perilla frutescens (L.) is an important medicinal and edible plant in China with nutritional and medical uses. The extract from leaves of Perilla frutescens contains flavonoids and volatile oils, which are mainly used in traditional Chinese medicine. In this study, we analyzed the transcriptomic and metabolomic data of the leaves of two Perilla frutescens varieties: JIZI 1 and JIZI 2. A total of 9277 differentially expressed genes and 223 flavonoid metabolites were identified in these varieties. Chrysoeriol, apigenin, malvidin, cyanidin, kaempferol, and their derivatives were abundant in the leaves of Perilla frutescens, which were more than 70% of total flavonoid contents. A total of 77 unigenes encoding 15 enzymes were identified as candidate genes involved in flavonoid biosynthesis in the leaves of Perilla frutescens. High expression of the CHS gene enhances the accumulation of flavonoids in the leaves of Perilla frutescens. Our results provide valuable information on the flavonoid metabolites and candidate genes involved in the flavonoid biosynthesis pathways in the leaves of Perilla frutescens.

Flavonoids and Isoflavonoids Biosynthesis in the Model Legume Lotus japonicus; Connections to Nitrogen Metabolism and Photorespiration

Phenylpropanoid metabolism represents an important metabolic pathway from which originates a wide number of secondary metabolites derived from phenylalanine or tyrosine, such as flavonoids and isoflavonoids, crucial molecules in plants implicated in a large number of biological processes. Therefore, various types of interconnection exist between different aspects of nitrogen metabolism and the biosynthesis of these compounds. For legumes, flavonoids and isoflavonoids are postulated to play pivotal roles in adaptation to their biological environments, both as defensive compounds (phytoalexins) and as chemical signals in symbiotic nitrogen fixation with rhizobia. In this paper, we summarize the recent progress made in the characterization of flavonoid and isoflavonoid biosynthetic pathways in the model legume Lotus japonicus (Regel) Larsen under different abiotic stress situations, such as drought, the impairment of photorespiration and UV-B irradiation. Emphasis is placed on results obtained using photorespiratory mutants deficient in glutamine synthetase. The results provide different types of evidence showing that an enhancement of isoflavonoid compared to standard flavonol metabolism frequently occurs in Lotus under abiotic stress conditions. The advance produced in the analysis of isoflavonoid regulatory proteins by the use of co-expression networks, particularly MYB transcription factors, is also described. The results obtained in Lotus japonicus plants can be also extrapolated to other cultivated legume species, such as soybean, of extraordinary agronomic importance with a high impact in feeding, oil production and human health.

Transcriptome analysis of Pueraria candollei var. mirifica for gene discovery in the biosyntheses of isoflavones and miroestrol

BACKGROUND

Pueraria candollei var. mirifica, a Thai medicinal plant used traditionally as a rejuvenating herb, is known as a rich source of phytoestrogens, including isoflavonoids and the highly estrogenic miroestrol and deoxymiroestrol. Although these active constituents in P. candollei var. mirifica have been known for some time, actual knowledge regarding their biosynthetic genes remains unknown.

RESULTS

Miroestrol biosynthesis was reconsidered and the most plausible mechanism starting from the isoflavonoid daidzein was proposed. A de novo transcriptome analysis was conducted using combined P. candollei var. mirifica tissues of young leaves, mature leaves, tuberous cortices, and cortex-excised tubers. A total of 166,923 contigs was assembled for functional annotation using protein databases and as a library for identification of genes that are potentially involved in the biosynthesis of isoflavonoids and miroestrol. Twenty-one differentially expressed genes from four separate libraries were identified as candidates involved in these biosynthetic pathways, and their respective expressions were validated by quantitative real-time reverse transcription polymerase chain reaction. Notably, isoflavonoid and miroestrol profiling generated by LC-MS/MS was positively correlated with expression levels of isoflavonoid biosynthetic genes across the four types of tissues. Moreover, we identified R2R3 MYB transcription factors that may be involved in the regulation of isoflavonoid biosynthesis in P. candollei var. mirifica. To confirm the function of a key-isoflavone biosynthetic gene, P. candollei var. mirifica isoflavone synthase identified in our library was transiently co-expressed with an Arabidopsis MYB12 transcription factor (AtMYB12) in Nicotiana benthamiana leaves. Remarkably, the combined expression of these proteins led to the production of the isoflavone genistein.

CONCLUSIONS

Our results provide compelling evidence regarding the integration of transcriptome and metabolome as a powerful tool for identifying biosynthetic genes and transcription factors possibly involved in the isoflavonoid and miroestrol biosyntheses in P. candollei var. mirifica.

Identification of novel MYB transcription factors involved in the isoflavone biosynthetic pathway by using the combination screening system with agroinfiltration and hairy root transformation

Soybean isoflavones are functionally important secondary metabolites that are mainly accumulated in seeds. Their biosynthetic processes are regulated coordinately at the transcriptional level; however, screening systems for key transcription factors (TFs) are limited. Here we developed a combination screening system comprising a simple agroinfiltration assay and a robust hairy root transformation assay. First, we screened for candidate MYB TFs that could activate the promoters of the chalcone synthase (CHS) gene GmCHS8 and the isoflavone synthase (IFS) genes GmIFS1 and GmIFS2 in the isoflavone biosynthetic pathway. In the agroinfiltration assay, we co-transformed a LjUbi (Lotus japonicus polyubiquitin gene) promoter-fused MYB gene with target promoter-fused GUS (β-glucuronidase) gene constructs, and identified three genes (GmMYB102, GmMYB280, and GmMYB502) as candidate regulators of isoflavone biosynthesis. We then evaluated the functional regulatory role of identified three MYB genes in isoflavone biosynthesis using hairy roots transformation assay in soybean for the accumulation of isoflavones. Three candidate MYB genes showed an increased accumulation of total isoflavones in hairy root transgenic lines. Accumulation of total isoflavones in the three MYB-overexpressing lines was approximately 2-to 4-folds more than that in the vector control, confirming their possible role to regulate isoflavone biosynthesis. However, the significant accumulation of authentic GmCHS8, GmIFS1, and GmIFS2 transcripts could not be observed except for the GmMYB502-overexpressing line. Therefore, the analysis of isoflavone accumulation in transgenic hairy root was effective for evaluation of transactivation activity of MYB TFs for isoflavone biosynthetic genes. Our results demonstrate a simple and robust system that can potentially identify the function of orphan TFs in diverse plant metabolic pathways.

Enhanced Salt Tolerance of Rhizobia-inoculated Soybean Correlates with Decreased Phosphorylation of the Transcription Factor GmMYB183 and Altered Flavonoid Biosynthesis

Soybean (Glycine max (L.) Merrill) is an important component of the human diet and animal feed, but soybean production is limited by abiotic stresses especially salinity. We recently found that rhizobia inoculation enhances soybean tolerance to salt stress, but the underlying mechanisms are unaddressed. Here, we used quantitative phosphoproteomic and metabonomic approaches to identify changes in phosphoproteins and metabolites in soybean roots treated with rhizobia inoculation and salt. Results revealed differential regulation of 800 phosphopeptides, at least 32 of these phosphoproteins or their homologous were reported be involved in flavonoid synthesis or trafficking, and 27 out of 32 are transcription factors. We surveyed the functional impacts of all these 27 transcription factors by expressing their phospho-mimetic/ablative mutants in the roots of composite soybean plants and found that phosphorylation of GmMYB183 could affect the salt tolerance of the transgenic roots. Using data mining, ChIP and EMSA, we found that GmMYB183 binds to the promoter of the soybean GmCYP81E11 gene encoding for a Cytochrome P450 monooxygenase which contributes to the accumulation of ononin, a monohydroxy B-ring flavonoid that negatively regulates soybean tolerance to salinity. Phosphorylation of GmMYB183 was inhibited by rhizobia inoculation; overexpression of GmMYB183 enhanced the expression of GmCYP81E11 and rendered salt sensitivity to the transgenic roots; plants deficient in GmMYB183 function are more tolerant to salt stress as compared with wild-type soybean plants, these results correlate with the transcriptional induction of GmCYP81E11 by GmMYB183 and the subsequent accumulation of ononin. Our findings provide molecular insights into how rhizobia enhance salt tolerance of soybean plants.

Conserved miRNAs modulate the expression of potential transcription factors of isoflavonoid biosynthetic pathway in soybean seeds

Despite the significant importance of soybean isoflavone, the regulatory mechanism of miRNAs during its biosynthesis is highly unexplored. In the present work, nine existing miRNAs along with their ten corresponding target genes were identified and validated in soybean for their possible role during isoflavonoid biosynthesis and accumulation. Temporal expression analysis at four key stages of seed development (35, 45, 55 and 65DAF) of all the miRNA-target pairs showed varying degree of differential accumulation in two soybean genotypes (NRC37: high isoflavone; and NRC7: low isoflavone). Differential expression of MYB65-Gma-miR159, MYB96-Gma-miRNA1534, MYB176-Gma-miRNA5030, SPL9-Gma-miRNA156, TCP3, TCP4-Gma-miRNA319, WD40-Gma-miRNA162, UDP-glucose: flavonoid 3-O-glucosyltransferase-Gma-miRNA396, and CHI3-Gma-miRNA5434 showed an important relationship with their targets in both the soybean genotypes across all the stages. Therefore, the finding of the present work would certainly increase our understanding of molecular regulation of isoflavone biosynthetic pathway mediated by the miRNA which would guide molecular breeder to develop isoflavone rich soybean cultivars.

Comprehensive transcriptome analysis reveals genes potentially involved in isoflavone biosynthesis in Pueraria thomsonii Benth

Pueraria thomsonii Benth is an important medicinal plant. Transcriptome sequencing, unigene assembly, the annotation of transcripts and the study of gene expression profiles play vital roles in gene function research. However, the full-length transcriptome of P. thomsonii remains unknown. Here, we obtained 44,339 nonredundant transcripts of P. thomsonii by using the PacBio RS II Isoform and Illumina sequencing platforms, of which 43,195 were annotated genes. Compared with the expression levels in the plant roots, those of transcripts with a |fold change| ≥ 4 and FDR < 0.01 in the leaves or stems were assigned as differentially expressed transcripts (DETs). In total, we found 9,225 DETs, 32 of which came from structural genes that were potentially involved in isoflavone biosynthesis. The expression profiles of 8 structural genes from the RNA-Seq data were validated by qRT-PCR. We identified 437 transcription factors (TFs) that were positively or negatively correlated with at least 1 of the structural genes involved in isoflavone biosynthesis using Pearson correlation coefficients (r) (r > 0.8 or r < -0.8). We also identified a total of 32 microRNAs (miRNAs), which targeted 805 transcripts. These miRNAs caused enriched function in ‘ATP binding’, ‘defense response’, ‘ADP binding’, and ‘signal transduction’. Interestingly, MIR156a potentially promoted isoflavone biosynthesis by repressing SBP, and MIR319 promoted isoflavone biosynthesis by repressing TCP and HB-HD-ZIP. Finally, we identified 2,690 alternative splicing events, including that of the structural genes of trans-cinnamate 4-monooxygenase and pullulanase, which are potentially involved in the biosynthesis of isoflavone and starch, respectively, and of three TFs potentially involved in isoflavone biosynthesis. Together, these results provide us with comprehensive insight into the gene expression and regulation of P. thomsonii.

GmMYB176 Regulates Multiple Steps in Isoflavonoid Biosynthesis in Soybean

Isoflavonoids are a group of plant natural compounds synthesized almost exclusively by legumes, and are abundant in soybean seeds and roots. They play important roles in plant-microbial interactions and the induction of nod gene expression in Rhizobia that form nitrogen-fixing nodules on soybean roots. Isoflavonoids also contribute to the positive health effects associated with soybean consumption by humans and animals. An R1 MYB transcription factor GmMYB176 regulates isoflavonoid biosynthesis by activating chalcone synthase (CHS) 8 gene expression in soybean. Using a combination of transcriptomic and metabolomic analyses of GmMYB176-RNAi silenced (GmMYB176-Si), GmMYB176-overexpressed (GmMYB176-OE), and control soybean hairy roots, we identified a total of 33 differentially expressed genes (DEGs) and 995 differentially produced metabolite features (DPMF) in GmMYB176-Si hairy roots, and 5727 DEGs and 149 DPMFs in GmMYB176-OE hairy roots. By a targeted approach, 25 isoflavonoid biosynthetic genes and 6 metabolites were identified as differentially regulated in GmMYB176-OE and GmMYB176-Si soybean hairy roots. Taken together, our results demonstrate the complexity of isoflavonoid biosynthesis in soybean roots and suggest that a coordinated expression of pathway genes, substrate flux and product threshold level may contribute to the dynamic of the pathway regulation.

GmMYB58 and GmMYB205 are seed-specific activators for isoflavonoid biosynthesis in Glycine max

GmMYB58 and GmMYB205 are key positive regulators that are involved in isoflavonoid biosynthesis in seeds of Glycine max, and they activate the expression of several structural genes in the isoflavonoid pathway. MYB transcription factors (TFs) are major regulators involved in flavonoid/isoflavonoid biosynthesis in many plant species. However, functions of most MYB TFs remain unknown in flavonoid/isoflavonoid pathway in Glycine max. In this study, we identified 321 MYB TFs by genome-wide searching, and further isolated and functionally characterized two MYB TFs, GmMYB58 and GmMYB205. The deduced GmMYB58 and GmMYB205 proteins contain highly conserved R2R3 repeat domain at the N-terminal region that is the signature motif of R2R3-type MYB TFs. GmMYB58 and GmMYB205 were highly expressed in early seed development stages than in the other tested organs. GmMYB58 and GmMYB205 GFP fusion proteins were found to be localized in the nucleus when they were transiently expressed in Arabidopsis thaliana mesophyll protoplast. Both GmMYB58 and GmMYB205 can activate the promoter activities of GmCHS, GmIFS2, and GmHID in the transient trans-activation assays, and the activation of GmHID by both GmMYB58 and GmMYB205 was further confirmed by yeast one-hybrid assay. In addition, over-expression of GmMYB58 and GmMYB205 resulted in significant increases in expression levels of several pathway genes in soybean hairy roots, in particular, IFS2 by more than fivefolds in GmMYB205-over-expressing lines. Moreover, isoflavonoid contents were remarkably enhanced in the GmMYB58 and GmMYB205 over-expressing hairy roots than in the control. Our results suggest that GmMYB58 and GmMYB205 are seed-specific TFs, and they can enhance isoflavonoid biosynthesis mainly through the regulation of GmIFS2 and GmHID in G. max.

Isolation and characterization of GmMYBJ3, an R2R3-MYB transcription factor that affects isoflavonoids biosynthesis in soybean

Isoflavonoids are secondary metabolites that play a variety of roles in plant-microbe interactions and plant defenses against abiotic stresses. Here we report a new MYB transcription factor (TF) gene, GmMYBJ3, that is involved in the isoflavonoids biosynthesis. The GmMYBJ3 gene is 1,002 bp long and encodes a protein of 333 amino acids. Amino acid sequence analysis showed that GmMYBJ3 is a typical R2R3 MYB TF. Yeast expression experiment demonstrated that GmMYBJ3 has its transcription activity in the nucleus and is transiently expressed in onion epidermal cells. The GmMYBJ3 gene was transformed into soybean and the expression activity of the GmMYBJ3 gene was significantly positively correlated with total isoflavonoid accumulation in soybean. Transient expression assays indicated that GmMYBJ3 can activate CHS8 expression. Furthermore, we analyzed the expressions of several genes known involved in the isoflavonoid biosynthesis, including CHS8, CHI1A, PAL1, IFS2 and F3H, in the GmMYBJ3 transgenic plants. The results showed that the expression levels of CHS8 and CHI1A were significantly increased in the transgenic plants compared to wild-type plants, but those of PAL1, IFS2 and F3H remained similar between the transgenic and wild-type plants. These results suggest that GmMYBJ3 participates in the isoflavonoid biosynthesis through regulation of CHS8 and CHI1A in soybean.

BnLATE, a Cys2/His2-Type Zinc-Finger Protein, Enhances Silique Shattering Resistance by Negatively Regulating Lignin Accumulation in the Silique Walls of Brassica napus

Silique shattering resistance is one of the most important agricultural traits in oil crop breeding. Seed shedding from siliques prior to and during harvest causes devastating losses in oilseed yield. Lignin biosynthesis in the silique walls is thought to affect silique-shattering resistance in oil crops. Here, we identified and characterized B. napus LATE FLOWERING (BnLATE), which encodes a Cys2/His2-type zinc-finger protein. Heterologous expression of BnLATE under the double enhanced CaMV 35S promoter (D35S) in wild-type Arabidopsis plants resulted in a marked decrease in lignification in the replum, valve layer (carpel) and dehiscence zone. pBnLATE::GUS activity was strong in the yellowing silique walls of transgenic lines. Furthermore, the expression pattern of BnLATE and the lignin content gradient in the silique walls at 48 days after pollination (DAP) of 73290, a B. napus silique shattering-resistant line, are similar to those in transgenic Arabidopsis lines expressing BnLATE. Transcriptome sequencing of the silique walls revealed that genes encoding peroxidases, which polymerize monolignols and lignin in the phenylpropanoid pathway, were down-regulated at least two-fold change in the D35S::BnLATE transgenic lines. pBnLATE::BnLATE transgenic lines were further used to identify the function of BnLATE, and the results showed that lignification in the carpel and dehiscence zone of yellowing silique also remarkably decreased compared with the wild-type control, the silique shattering-resistance and expression pattern of peroxidase genes are very similar to results with D35S::BnLATE. These results suggest that BnLATE is a negative regulator of lignin biosynthesis in the yellowing silique walls, and promotes silique-shattering resistance in B. napus through restraining the polymerization of monolignols and lignin.

Transcriptomics, Targeted Metabolomics and Gene Expression of Blackberry Leaves and Fruits Indicate Flavonoid Metabolic Flux from Leaf to Red Fruit

Blackberries (Rubus spp.) are among the high added value food products relevant for human health due to the increasing evidence of the beneficial effects of polyphenols, which are very abundant in these fruits. Interestingly, these compounds also play a role on plant physiology, being especially relevant their role in plant defense against biotic and abiotic stress. Hence, we hypothesize that since blackberry fruits have high amounts of flavonols and anthocyanins, leaves would also have high amounts of these compounds, and can be studied as a source of active molecules; furthermore, leaf synthesis would support their high contents in fruits. To explore this hypothesis, the present study reports a de novo transcriptome analysis on field grown blackberry leaves and fruits at the same time point, to establish the metabolic relationship of these compounds in both organs. Transcripts were aligned against Fragaria vesca genome, and genes were identified and annotated in different databases; tissue expression pattern showed 20,463 genes common to leaves and fruits, while 6,604 genes were significantly overexpressed only in fruits, while another 6,599 genes were significantly overexpressed in leaves, among which flavonol-anthocyanin transporter genes were present. Bioactives characterization indicated that total phenolics in leaves were three-fold, and flavonols were six-fold than in fruits, while concentration of anthocyanins was higher in fruits; HPLC-MS analysis indicated different composition in leaves and fruits, with cyanidin-3-glucoside as the only common compound identified. Next, RT-qPCR of the core genes in the flavonol anthocyanin pathway and regulatory MYB genes were carried out. Interestingly, genes in the flavonol-anthocyanin pathway and flavonol-transport families were overexpressed in leaves, consistent with the higher bioactive levels. On the other hand, transcription factors were overexpressed in fruits anticipating an active anthocyanin biosynthesis upon ripening. This suggests that, in addition to the biosynthesis taking place in the fruits during ripening, translocation of flavonols from leaves to fruits contributes to the high amounts of bioactives starting to accumulate in fruits.

Isolation and characterization of a R2R3-MYB transcription factor gene related to anthocyanin biosynthesis in the spathes of Anthurium andraeanum (Hort.)

KEY MESSAGE

A R2R3-MYB gene AaMYB2 was isolated from Anthurium andraeanum (Hort.) and was functionally characterized to be a positive transcriptional regulator for anthocyanin biosynthesis in the spathes and leaves. Spathe coloration is an important Anthurium andraeanum (Hort.) characteristic, which is mainly contributed by anthocyanins. R2R3-MYB transcription factors (TFs) are important regulators of anthocyanin biosynthesis in plants. Here we describe the identification and characterization of AaMYB2, a member in subgroup 6 of the R2R3-MYB TFs family, which correlated with anthocyanin biosynthesis in A. andraeanum. AaMYB2 was a nuclear-localization protein with positive transcriptional activity, and prominently expressed in the red spathes. Ectopic expression of AaMYB2 in tobacco led to anthocyanin accumulation and up-regulation of the early and late anthocyanin pathway genes, particularly NtDFR, NtANS, and NtUFGT, and the endogenous TF genes NtAn2 and NtAn1 in leaves. In the developing red spathes of 'Tropical' and 'Vitara', the expression of AaMYB2 was closely linked to anthocyanin accumulation, and co-expressed with AaCHS, AaF3H, and AaANS, the latter two of which were regarded as the potential targets of the R locus encoding a TF controlling spathe colors inheritance in anthurium. In addition, the transcription level of AaMYB2 in various cultivars with different color phenotypes showed that AaMYB2 was drastically expressed in the spathes from the red, pink, and purple cultivars, but hardly detected in the spathes from the white and green ones. Besides, AaMYB2 also showed higher expression in newly developmental leaves when anthocyanin was actively biosynthesized. Taken together, AaMYB2 positively related to anthocyanin biosynthesis in anthurium spathes and leaves, and appeared to regulate the expression of AaF3H, AaANS, and possibly AaCHS.

Application of Pseudomonas fluorescens to Blackberry under Field Conditions Improves Fruit Quality by Modifying Flavonoid Metabolism

Application of a plant growth promoting rhizobacterium (PGPR), Pseudomonas fluorescens N21.4, to roots of blackberries (Rubus sp.) is part of an optimised cultivation practice to improve yields and quality of fruit throughout the year in this important fruit crop. Blackberries are especially rich in flavonoids and therefore offer potential benefits for human health in prevention or amelioration of chronic diseases. However, the phenylpropanoid pathway and its regulation during ripening have not been studied in detail, in this species. PGPR may trigger flavonoid biosynthesis as part of an induced systemic response (ISR) given the important role of this pathway in plant defence, to cause increased levels of flavonoids in the fruit. We have identified structural genes encoding enzymes of the phenylpropanoid and flavonoid biosynthetic pathways catalysing the conversion of phenylalanine to the final products including flavonols, anthocyanins and catechins from blackberry, and regulatory genes likely involved in controlling the activity of pathway branches. We have also measured the major flavonols, anthocyanins and catechins at three stages during ripening. Our results demonstrate the coordinated expression of flavonoid biosynthetic genes with the accumulation of anthocyanins, catechins, and flavonols in developing fruits of blackberry. Elicitation of blackberry plants by treatment of roots with P.fluorescens N21.4, caused increased expression of some flavonoid biosynthetic genes and an accompanying increase in the concentration of selected flavonoids in fruits. Our data demonstrate the physiological mechanisms involved in the improvement of fruit quality by PGPR under field conditions, and highlight some of the genetic targets of elicitation by beneficial bacteria.

Identification and expression analyses of MYB and WRKY transcription factor genes in Papaver somniferum L

Papaver somniferum L. is an herbaceous, annual and diploid plant that is important from pharmacological and strategic point of view. The cDNA clones of two putative MYB and WRKY genes were isolated (GeneBank accession numbers KP411870 and KP203854, respectively) from this plant, via the nested-PCR method, and characterized. The MYB transcription factor (TF) comprises 342 amino acids, and exhibits the structural features of the R2R3MYB protein family. The WRKY TF, a 326 amino acid-long polypeptide, falls structurally into the group II of WRKY protein family. Quantitative real-time PCR (qRT-PCR) analyses indicate the presence of these TFs in all organs of P. somniferum L. and Papaver bracteatum L. Highest expression levels of these two TFs were observed in the leaf tissues of P. somniferum L. while in P. bracteatum L. the espression levels were highest in the root tissues. Promoter analysis of the 10 co-expressed gene clustered involved in noscapine biosynthesis pathway in P. somniferum L. suggested that not only these 10 genes are co-expressed, but also share common regulatory motifs and TFs including MYB and WRKY TFs, and that may explain their common regulation.

The soybean R2R3 MYB transcription factor GmMYB100 negatively regulates plant flavonoid biosynthesis

Soybean flavonoids, a group of important signaling molecules in plant-environment interaction, ubiquitously exist in soybean and are tightly regulated by many genes. Here we reported that GmMYB100, a gene encoding a R2R3 MYB transcription factor, is involved in soybean flavonoid biosynthesis. GmMYB100 is mainly expressed in flowers, leaves and immature embryo, and its level is decreased after pod ripening. Subcellular localization assay indicates that GmMYB100 is a nuclear protein. GmMYB100 has transactivation ability revealed by a yeast functional assay; whereas bioinformatic analysis suggests that GmMYB100 has a negative function in flavonoid biosynthesis. GmMYB100-overexpression represses the transcript levels of flavonoid-related genes in transgenic soybean hairy roots and Arabidopsis, and inhibits isoflavonoid (soybean) and flavonol (Arabidopsis) production in transgenic plants. Furthermore, the transcript levels of six flavonoid-related genes and flavonoid (isoflavonoid and flavone aglycones) accumulation are elevated in the GmMYB100-RNAi transgenic hairy roots. We also demonstrate that GmMYB100 protein depresses the promoter activities of soybean chalcone synthase and chalcone isomerase. These findings indicate that GmMYB100 is a negative regulator in soybean flavonoid biosynthesis pathway.

MYB Transcription Factors as Regulators of Phenylpropanoid Metabolism in Plants

Phenylpropanoid-derived compounds represent a diverse family of secondary metabolites that originate from phenylalanine. These compounds have roles in plant growth and development, and in defense against biotic and abiotic stress. Many of these compounds are also beneficial to human health and welfare. V-myb myeloblastosis viral oncogene homolog (MYB) proteins belong to a large family of transcription factors and are key regulators of the synthesis of phenylpropanoid-derived compounds. This review summarizes the current understanding of MYB proteins and their roles in the regulation of phenylpropanoid metabolism in plants.