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

Populus euphratica R2R3-MYB transcription factor RAX2 binds ANN1 promoter to increase cadmium enrichment in Arabidopsis

The expression of R2R3-MYB transcription factor PeRAX2 increased transiently upon CdCl2 exposure (100 μM, 48 h) in leaves and roots of Populus euphratica. We observed that overexpression of PeRAX2 increased Cd2+ concentration in Arabidopsis root cells and Cd2+ amount in whole plant, which was due to the increased Cd2+ influx into root tips. However, the Cd2+ influx facilitated by PeRAX2 overexpression was substantially reduced by LaCl3 (an inhibitor of Ca2+-channels), suggesting that PeRAX2 could promote the Cd2+ entering through PM Ca2+-permeable channels (CaPCs) in the roots. It is noting that the expression of annexin1 (AtANN1), which mediates the influx of divalent cations through the PM calcium channels, was upregulated by Cd2+ in PeRAX2-transgenic Arabidopsis. Bioinformatic analysis revealed that the AtANN1 promoter (AtANN1-pro) contains four cis-elements for MYB binding. The PeRAX2 interaction with AtANN1-pro was validated by LUC reporter assay, EMSA, and Y1H assay. Our data showed that PeRAX2 binds to the AtANN1 promoter region to regulate gene transcription and that AtANN1 mediates the Cd2+ entry through CaPCs in the PM, leading to a Cd2+ enrichment in transgenic plants. The PeRAX2-stimulated Cd2+ enrichment consequently resulted in high H2O2 production in root cells of transgenic plants. The expression of AtSOD and AtPOD and activities of CAT, SOD, POD increased in the transgenic lines under Cd2+ stress. However, the Cd2+-upregulated expression and activity of antioxidative enzymes were less pronounced in the PeRAX2-overexpressed lines, compared to the wildtype and vector controls. As a result, root length and plant growth were more suppressed by Cd2+ in the transgenic lines. Our data suggest that transcriptional regulation of AtANN1 by PeRAX2 can be utilized to improve Cd2+ enrichment and phytoremediation, although the enriched Cd2+ affected antioxidant defense system and plant growth in the model species.

Jasmonate enhances cold acclimation in jojoba by promoting flavonol synthesis

Jojoba is an industrial oil crop planted in tropical arid areas, and its low-temperature sensitivity prevents its introduction into temperate areas. Studying the molecular mechanisms associated with cold acclimation in jojoba is advantageous for developing breeds with enhanced cold tolerance. In this study, metabolomic analysis revealed that various flavonols accumulate in jojoba during cold acclimation. Time-course transcriptomic analysis and weighted correlation network analysis (WGCNA) demonstrated that flavonol biosynthesis and jasmonates (JAs) signaling pathways played crucial roles in cold acclimation. Combining the biochemical and genetic analyses showed that ScMYB12 directly activated flavonol synthase gene (ScFLS). The interaction between ScMYB12 and transparent testa 8 (ScTT8) promoted the expression of ScFLS, but the negative regulator ScJAZ13 in the JA signaling pathway interacted with ScTT8 to attenuate the transcriptional activity of the ScTT8 and ScMYB12 complex, leading to the downregulation of ScFLS. Cold acclimation stimulated the production of JA in jojoba leaves, promoted the degradation of ScJAZ13, and activated the transcriptional activity of ScTT8 and ScMYB12 complexes, leading to the accumulation of flavonols. Our findings reveal the molecular mechanism of JA-mediated flavonol biosynthesis during cold acclimation in jojoba and highlight the JA pathway as a promising means for enhancing cold tolerance in breeding efforts.

A Sweet Potato MYB Transcription Factor IbMYB330 Enhances Tolerance to Drought and Salt Stress in Transgenic Tobacco

MYB transcription factors (TFs) play vital roles in plant growth, development, and response to adversity. Although the MYB gene family has been studied in many plant species, there is still little known about the function of R2R3 MYB TFs in sweet potato in response to abiotic stresses. In this study, an R2R3 MYB gene, IbMYB330 was isolated from sweet potato (Ipomoea batatas). IbMYB330 was ectopically expressed in tobacco and the functional characterization was performed by overexpression in transgenic plants. The IbMYB330 protein has a 268 amino acid sequence and contains two highly conserved MYB domains. The molecular weight and isoelectric point of IbMYB330 are 29.24 kD and 9.12, respectively. The expression of IbMYB330 in sweet potato is tissue-specific, and levels in the root were significantly higher than that in the leaf and stem. It showed that the expression of IbMYB330 was strongly induced by PEG-6000, NaCl, and H2O2. Ectopic expression of IbMYB330 led to increased transcript levels of stress-related genes such as SOD, POD, APX, and P5CS. Moreover, compared to the wild-type (WT), transgenic tobacco overexpression of IbMYB330 enhanced the tolerance to drought and salt stress treatment as CAT activity, POD activity, proline content, and protein content in transgenic tobacco had increased, while MDA content had decreased. Taken together, our study demonstrated that IbMYB330 plays a role in enhancing the resistance of sweet potato to stresses. These findings lay the groundwork for future research on the R2R3-MYB genes of sweet potato and indicates that IbMYB330 may be a candidate gene for improving abiotic stress tolerance in crops.

Identification and Expression Analysis of R2R3-MYB Transcription Factors Associated with Flavonoid Biosynthesis in Panax quinquefolius

Panax quinquefolius L. is an important medicinal plant, and flavonoids are among its main secondary metabolites. The R2R3-MYB transcription factor plays an irreplaceable role in plant growth, development, and secondary metabolism. In our study, we identified 159 R2R3-MYBs and analyzed their physical and chemical properties in P. quinquefolius. The protein length of 159 PqMYBs varied from 107 to 1050 amino acids. The molecular weight ranged from 12.21 to 116.44 kDa. The isoelectric point was between 4.57 and 10.34. We constructed a phylogenetic tree of P. quinquefolius and Arabidopsis thaliana R2R3-MYB family members, and PqMYB members were divided into 33 subgroups. Transcriptome data analysis showed that the expression patterns of PqMYBs in root, leaf, and flower were significantly different. Following the MeJA treatment of seedlings, five candidate PqMYB genes demonstrated a response. A correlation analysis of PqMYBs and candidate flavonoid pathway genes showed that PqMYB2, PqMYB46, and PqMYB72 had correlation coefficients that were higher than 0.8 with PqCHS, PqANS4, and PqCCoAMT10, respectively. Furthermore, a transient expression assay confirmed that the three PqMYBs were localized in the nucleus. We speculated that these three PqMYBs were related to flavonoid biosynthesis in P. quinquefolius. These results provided a theoretical basis and a new perspective for further understanding the R2R3-MYB gene family and the biosynthesis mechanism of secondary metabolites in P. quinquefolius.

Genome-Wide Identification and Characterization of MYB Gene Family and Analysis of Its Sex-Biased Expression Pattern in Spinacia oleracea L

The members of the myeloblastosis (MYB) family of transcription factors (TFs) participate in a variety of biological regulatory processes in plants, such as circadian rhythm, metabolism, and flower development. However, the characterization of MYB genes across the genomes of spinach Spinacia oleracea L. has not been reported. Here, we identified 140 MYB genes in spinach and described their characteristics using bioinformatics approaches. Among the MYB genes, 54 were 1R-MYB, 80 were 2R-MYB, 5 were 3R-MYB, and 1 was 4R-MYB. Almost all MYB genes were located in the 0-30 Mb region of autosomes; however, the 20 MYB genes were enriched at both ends of the sex chromosome (chromosome 4). Based on phylogeny, conserved motifs, and the structure of genes, 2R-MYB exhibited higher conservation relative to 1R-MYB genes. Tandem duplication and collinearity of spinach MYB genes drive their evolution, enabling the functional diversification of spinach genes. Subcellular localization prediction indicated that spinach MYB genes were mainly located in the nucleus. Cis-acting element analysis confirmed that MYB genes were involved in various processes of spinach growth and development, such as circadian rhythm, cell differentiation, and reproduction through hormone synthesis. Furthermore, through the transcriptome data analysis of male and female flower organs at five different periods, ten candidate genes showed biased expression in spinach males, suggesting that these genes might be related to the development of spinach anthers. Collectively, this study provides useful information for further investigating the function of MYB TFs and novel insights into the regulation of sex determination in spinach.

Comprehensive Analysis of WUSCEL-Related Homeobox Gene Family in Ramie (Boehmeria nivea) Indicates Its Potential Role in Adventitious Root Development

A WUSCHEL-related homeobox (WOX) gene family has been implicated in promoting vegetative organs to embryonic transition and maintaining plant embryonic stem cell identity. Using genome-wide analysis, we identified 17 candidates, WOX genes in ramie (Boehmeria nivea). The genes (BnWOX) showed highly conserved homeodomain regions typical of WOX. Based on phylogenetic analysis, they were classified into three distinct groups: modern, intermediate, and ancient clades. The genes displayed 65% and 35% collinearities with their Arabidopsis thaliana and Oryza sativa ortholog, respectively, and exhibited similar motifs, suggesting similar functions. Furthermore, four segmental duplications (BnWOX10/14, BnWOX13A/13B, BnWOX9A/9B, and BnWOX6A/Maker00021031) and a tandem-duplicated pair (BnWOX5/7) among the putative ramie WOX genes were obtained, suggesting that whole-genome duplication (WGD) played a role in WOX gene expansion. Expression profiling analysis of the genes in the bud, leaf, stem, and root of the stem cuttings revealed higher expression levels of BnWOX10 and BnWOX14 in the stem and root and lower in the leaf consistent with the qRT-PCR analysis, suggesting their direct roles in ramie root formation. Analysis of the rooting characteristics and expression in the stem cuttings of sixty-seven different ramie genetic resources showed a possible involvement of BnWOX14 in the adventitious rooting of ramie. Thus, this study provides valuable information on ramie WOX genes and lays the foundation for further research.

Comparative genome and metabolome analyses uncover the evolution and flavonoid biosynthesis between Apocynum venetum and Apocynum hendersonii

Apocynum species have great application prospects in textile and phytoremediation of saline soil, are rich in flavonoids, and possess medicinal significance. Here, we report the draft genomes of Apocynum venetum and Apocynum hendersonii, and elucidate their evolutionary relationship. The high synteny and collinearity between the two suggested that they have experienced the same WGD event. Comparative analysis revealed that flavone 3-hydroxylase (ApF3H) and differentially evolved flavonoid 3-O-glucosyl transferase (ApUFGT) genes are critical for determining natural variation in flavonoid biosynthesis between the species. Overexpression of ApF3H-1 enhanced the total flavonoid content and promoted the antioxidant capacity of transformed plants compared to the wild-type. ApUFGT5 and 6 explained the diversification of flavonoids or their derivatives. These data provide biochemical insight and knowledge on the genetic regulation of flavonoid biosynthesis, supporting the adoption of these genes in breeding programs aimed at the multipurpose utilization of the plants.

Advances in the study of the function and mechanism of the action of flavonoids in plants under environmental stresses

MAIN CONCLUSION

This review summarizes the anti-stress effects of flavonoids in plants and highlights its role in the regulation of polar auxin transport and free radical scavenging mechanism. As secondary metabolites widely present in plants, flavonoids play a vital function in plant growth, but also in resistance to stresses. This review introduces the classification, structure and synthetic pathways of flavonoids. The effects of flavonoids in plant stress resistance were enumerated, and the mechanism of flavonoids in plant stress resistance was discussed in detail. It is clarified that plants under stress accumulate flavonoids by regulating the expression of flavonoid synthase genes. It was also determined that the synthesized flavonoids are transported in plants through three pathways: membrane transport proteins, vesicles, and bound to glutathione S-transferase (GST). At the same time, the paper explores that flavonoids regulate polar auxin transport (PAT) by acting on the auxin export carrier PIN-FORMED (PIN) in the form of ATP-binding cassette subfamily B/P-glycoprotein (ABCB/PGP) transporter, which can help plants to respond in a more dominant form to stress. We have demonstrated that the number and location of hydroxyl groups in the structure of flavonoids can determine their free radical scavenging ability and also elucidated the mechanism by which flavonoids exert free radical removal in cells. We also identified flavonoids as signaling molecules to promote rhizobial nodulation and colonization of arbuscular mycorrhizal fungi (AMF) to enhance plant-microbial symbiosis in defense to stresses. Given all this knowledge, we can foresee that the in-depth study of flavonoids will be an essential way to reveal plant tolerance and enhance plant stress resistance.

Genome-Wide Analysis of MYB Transcription Factor Gene Superfamily Reveals BjPHL2a Involved in Modulating the Expression of BjCHI1 in Brassica juncea

Brassica juncea is an economically important vegetable and oilseed crop. The MYB transcription factor superfamily is one of the largest transcription factor families in plants, and plays crucial roles in regulating the expression of key genes involved in a variety of physiological processes. However, a systematic analysis of the MYB transcription factor genes in Brassica juncea (BjMYB) has not been performed. In this study, a total of 502 BjMYB superfamily transcription factor genes were identified, including 23 1R-MYBs, 388 R2R3-MYBs, 16 3R-MYBs, 4 4R-MYBs, 7 atypical MYBs, and 64 MYB-CCs, which is approximately 2.4-fold larger than that of AtMYBs. Phylogenetic relationship analysis revealed that the MYB-CC subfamily consists of 64 BjMYB-CC genes. The expression pattern of members of PHL2 subclade homologous genes in Brassica juncea (BjPHL2) after Botrytis cinerea infection were determined, and BjPHL2a was isolated from a yeast one-hybrid screen with the promoter of BjCHI1 as bait. BjPHL2a was found to localize mainly in the nucleus of plant cells. An EMSA assay confirmed that BjPHL2a binds to the Wbl-4 element of BjCHI1. Transiently expressed BjPHL2a activates expression of the GUS reporter system driven by a BjCHI1 mini-promoter in tobacco (Nicotiana benthamiana) leaves. Taken together, our data provide a comprehensive evaluation of BjMYBs and show that BjPHL2a, one of the members of BjMYB-CCs, functions as a transcription activator by interacting with the Wbl-4 element in the promoter of BjCHI1 for targeted gene-inducible expression.

Genome-Wide Analysis of AP2/ERF Gene Superfamily in Ramie (Boehmeria nivea L.) Revealed Their Synergistic Roles in Regulating Abiotic Stress Resistance and Ramet Development

AP2/ERF transcription factors (TFs) are one of the largest superfamilies in plants, and play vital roles in growth and response to biotic/abiotic stresses. Although the AP2/ERF family has been extensively characterized in many species, very little is known about this family in ramie (Boehmeria nivea L.). In this study, 138 AP2/ERF TFs were identified from the ramie genome and were grouped into five subfamilies, including the AP2 (19), RAV (5), Soloist (1), ERF (77), and DREB (36). Unique motifs were found in the DREB/ERF subfamily members, implying significance to the AP2/ERF TF functions in these evolutionary branches. Segmental duplication events were found to play predominant roles in the BnAP2/ERF TF family expansion. Light-, stress-, and phytohormone-responsive elements were identified in the promoter region of BnAP2/ERF genes, with abscisic acid response elements (ABRE), methyl jasmonate response elements, and the dehydration response element (DRE) being dominant. The integrated transcriptome and quantitative real-time PCR (qPCR) revealed 12 key BnAP2/ERF genes positively responding to waterlogging. Five of the genes are also involved in ramet development, with two (BnERF-30 and BnERF-32) further showing multifunctional roles. The protein interaction prediction analysis further verified their crosstalk mechanism in coordinating waterlogging resistance and ramet development. Our study provides new insights into the presence of AP2/ERF TFs in ramie, and provides candidate AP2/ERF TFs for further studies on breeding varieties with coupling between water stress tolerance and high yield.

Genome-Wide Identification and Analysis of the R2R3-MYB Gene Family in Theobroma cacao

The MYB gene family is involved in the regulation of plant growth, development and stress responses. In this paper, to identify Theobroma cacao R2R3-MYB (TcMYB) genes involved in environmental stress and phytohormones, we conducted a genome-wide analysis of the R2R3-MYB gene family in Theobroma cacao (cacao). A total of 116 TcMYB genes were identified, and they were divided into 23 subgroups according to the phylogenetic analysis. Meanwhile, the conserved motifs, gene structures and cis-acting elements of promoters were analyzed. Moreover, these TcMYB genes were distributed on 10 chromosomes. We conducted a synteny analysis to understand the evolution of the cacao R2R3-MYB gene family. A total of 37 gene pairs of TcMYB genes were identified through tandem or segmental duplication events. Additionally, we also predicted the subcellular localization and physicochemical properties. All the studies showed that TcMYB genes have multiple functions, including responding to environmental stresses. The results provide an understanding of R2R3-MYB in Theobroma cacao and lay the foundation for a further functional analysis of TcMYB genes in the growth of cacao.

Comparative transcriptomics of drought stress response of taproot meristem region of contrasting purple carrot breeding lines supported by physio-biochemical parameters

Carrot is one of the nutritious vegetable crops sensitive to drought stress resulting in loss of quality and yield. There are a lot of studies on detailed molecular mechanisms of drought stress response of main crops; however, very little information available on vegetables, including carrots. Hence, in this study, we investigated root transcriptome profiles from the meristematic region of two contrasting purple carrot (B7262A, drought tolerant; P1129, drought sensitive) lines under varying stress levels (85% and 70%) by using RNA-Seq technique. The morpho-physiological and biochemical response of B7262A line exhibited tolerance behavior to both DS (85% and 70%). RNA-Seq analysis revealed that 15,839 genes were expressed commonly in both carrot lines. The carrot line B7262A showed regulation of 514 genes in response to 85% DS, whereas P1129 showed differential regulation of 622 genes under 70% DS. The B7262A carrot line showed higher upregulation of transcripts that suggested its resilient behavior contrary to P1129 line. Furthermore, validation of transcript gene by qRT-PCR also confirmed the RNA-Seq analysis resulting in elevated expression levels of MYB48 transcription factor, MAPK mitogen-activated protein kinase ANP1, GER geraniol 8-hydroxylase, ABA ABA-induced in somatic embryo 3, FBOX putative F-box protein, FRO ferric reduction oxidase, and PDR probable disease resistance protein. Current study provided unprecedented insights of purple carrot lines that can be potentially exploited for the screening and development of resilient carrot.

Analysis of WRKY Resistance Gene Family in Boehmeria nivea (L.) Gaudich: Crosstalk Mechanisms of Secondary Cell Wall Thickening and Cadmium Stress

A total of 60 WRKY family genes of ramie were identified in the ramie. The genes were unevenly distributed across 14 chromosomes in the specie and highly concentrated (72%) in the distal telomeric region. Phylogenetic analysis placed these genes into seven distinct subfamilies groups: I, II (a, b, c, d, e), and III, with group IIc containing only the variant of heptapetide sequence (WRKYGKK). Segmental duplication events (41.7%) was found to be the main driver of BnGWRKY evolution. Thirty eight from among the genes showed collinear relationships with WRKY genes from Arabidopsis thaliana, Cannabis sativa, Oryza sativa, and Zea mays. The number and density of stress and hormone responsives cis-acting elements were comparably higher than other elements, with abundant ARE and rare LTR cis-acting elements indicating the long-standing adaptability of ramie to its natural environment. GO and KEGG enrichment analysis of the WRKY target genes revealed their involvement in response to stimuli, immune system processes, transporter protein activity and antioxidant activity. Expression analysis show that most WRKYs were activated by the cadmium stress, more especially the BnGWRKY2, BnGWRKY15, BnGWRKY20, BnGWRKY50 and BnGWRKY58. Combining transcriptome, orthologous gene relationships and qPCR result, we established the possible involvement of BnGWRKY50 and BnGWRKY58 in crosstalk mechanism between secondary cell wall thickening and Cd²⁺ stress. This provided information into the role of BnGWRKY proteins in ramie secondary wall development and cadmium stress response to, and could serve as basis for improvement of the ramie.