MYB Transcription Factors as Regulators of Secondary Metabolism in Plants

Identification and Functional Characterization of the SaMYB113 Gene in Solanum aculeatissimum

The MYB transcription factors (TFs) have substantial functions in anthocyanin synthesis as well as being widely associated with plant responses to various adversities. In the present investigation, we found an unreported MYB TF from Solanum aculeatissimum (a wild relative of eggplant) and named it SaMYB113 in reference to its homologous gene. Bioinformatics analysis demonstrated that the open reading frame of SaMYB113 was 825 bp in length, encoding 275 amino acids, with a typical R2R3-MYB gene structure, and predicted subcellular localization in the nucleus. Analysis of the tissue-specific expression pattern through qRT-PCR showed that the SaMYB113 was expressed at a high level in young stems as well as leaves of S. aculeatissimum. Transgenic Arabidopsis and tobacco plants overexpressing SaMYB113 pertinent to the control of the 35S promoter exhibited a distinct purple color trait, suggesting a significant change in their anthocyanin content. Furthermore, we obtained three tobacco transgenic lines with significant differences in anthocyanin accumulation and analyzed the differences in anthocyanin content by LC-MS/MS. The findings demonstrated that overexpression of SaMYB113 caused tobacco to have considerably raised levels of several anthocyanin components, with the most significant increases in delphinidin-like anthocyanins and cyanidin-like anthocyanins. The qRT-PCR findings revealed significant differences in the expression levels of structural genes for anthocyanin synthesis among various transgenic lines. In summary, this study demonstrated that the SaMYB113 gene has a substantial impact on anthocyanin synthesis, and overexpression of the SaMYB113 gene leads to significant modifications to the expression levels of a variety of anthocyanin-synthesizing genes, which leads to complex changes in anthocyanin content and affects plant phenotypes. This present research offers the molecular foundation for the research of the mechanism of anthocyanin formation within plants, as well as providing some reference for the improvement of traits in solanum crops.

The MYB family and their response to abiotic stress in ginger (Zingiber officinale Roscoe)

BACKGROUND

Zingiber officinale Roscoe, colloquially known as ginger, is a crop of significant medicinal and culinary value that frequently encounters adversity stemming from inhospitable environmental conditions. The MYB transcription factors have garnered recognition for their pivotal role in orchestrating a multitude of plant biological pathways. Nevertheless, the enumeration and characterization of the MYBs within Z. officinale Roscoe remains unknown. This study embarks on a genome-wide scrutiny of the MYB gene lineage in ginger, with the aim of cataloging all ZoMYB genes implicated in the biosynthesis of gingerols and curcuminoids, and elucidating their potential regulatory mechanisms in counteracting abiotic stress, thereby influencing ginger growth and development.

RESULTS

In this study, we identified an MYB gene family comprising 231 members in ginger genome. This ensemble comprises 74 singular-repeat MYBs (1R-MYB), 156 double-repeat MYBs (R2R3-MYB), and a solitary triple-repeat MYB (R1R2R3-MYB). Moreover, a comprehensive analysis encompassing the sequence features, conserved protein motifs, phylogenetic relationships, chromosome location, and gene duplication events of the ZoMYBs was conducted. We classified ZoMYBs into 37 groups, congruent with the number of conserved domains and gene structure analysis. Additionally, the expression profiles of ZoMYBs during development and under various stresses, including ABA, cold, drought, heat, and salt, were investigated in ginger utilizing both RNA-seq data and qRT-PCR analysis.

CONCLUSION

This work provides a comprehensive understanding of the MYB family in ginger and lays the foundation for the future investigation of the potential functions of ZoMYB genes in ginger growth, development and abiotic stress tolerance of ginger.

R2R3-MYB transcription factor CsMYB60 controls mature fruit skin color by regulating flavonoid accumulation in cucumber

Skin color is an important trait that determines the cosmetic appearance and quality of fruits. In cucumber, the skin color ranges from white to brown in mature fruits. However, the genetic basis for this important trait remains unclear. We conducted a genome-wide association study of natural cucumber populations, along with map-based cloning techniques, on an F2 population resulting from a cross between Pepino (with yellow-brown fruit skin) and Zaoer-N (with creamy fruit skin). We identified CsMYB60 as a candidate gene responsible for skin coloration in mature cucumber fruits. In cucumber accessions with white to pale yellow skin color, a premature stop mutation (C to T) was found in the second exon region of CsMYB60, whereas light yellow cucumber accessions exhibited splicing premature termination caused by an intronic mutator-like element insertion in CsMYB60. Transgenic CsMYB60c cucumber plants displayed a yellow-brown skin color by promoting accumulation of flavonoids, especially hyperoside, a yellow-colored flavonol. CsMYB60c encodes a nuclear protein that primarily acts as a transcriptional activator through its C-terminal activation motif. RNA sequencing and DNA affinity purification sequencing assays revealed that CsMYB60c promotes skin coloration by directly binding to the YYTACCTAMYT motif in the promoter regions of flavonoid biosynthetic genes, including CsF3'H, which encodes flavonoid 3'-hydroxylase. The findings of our study not only offer insight into the function of CsMYB60 as dominantly controlling fruit coloration, but also highlight that intronic DNA mutations can have a similar phenotypic impact as exonic mutations, which may be valuable in future cucumber breeding programs.

Regulatory microRNAs and phasiRNAs of paclitaxel biosynthesis in Taxus chinensis

Paclitaxel (trade name Taxol) is a rare diterpenoid with anticancer activity isolated from Taxus. At present, paclitaxel is mainly produced by the semi-synthetic method using extract of Taxus tissues as raw materials. The studies of regulatory mechanisms in paclitaxel biosynthesis would promote the production of paclitaxel through tissue/cell culture approaches. Here, we systematically identified 990 transcription factors (TFs), 460 microRNAs (miRNAs), and 160 phased small interfering RNAs (phasiRNAs) in Taxus chinensis to explore their interactions and potential roles in regulation of paclitaxel synthesis. The expression levels of enzyme genes in cone and root were higher than those in leaf and bark. Nearly all enzyme genes in the paclitaxel synthesis pathway were significantly up-regulated after jasmonate treatment, except for GGPPS and CoA Ligase. The expression level of enzyme genes located in the latter steps of the synthesis pathway was significantly higher in female barks than in male. Regulatory TFs were inferred through co-expression network analysis, resulting in the identification of TFs from diverse families including MYB and AP2. Genes with ADP binding and copper ion binding functions were overrepresented in targets of miRNA genes. The miRNA targets were mainly enriched with genes in plant hormone signal transduction, mRNA surveillance pathway, cell cycle and DNA replication. Genes in oxidoreductase activity, protein-disulfide reductase activity were enriched in targets of phasiRNAs. Regulatory networks were further constructed including components of enzyme genes, TFs, miRNAs, and phasiRNAs. The hierarchical regulation of paclitaxel production by miRNAs and phasiRNAs indicates a robust regulation at post-transcriptional level. Our study on transcriptional and posttranscriptional regulation of paclitaxel synthesis provides clues for enhancing paclitaxel production using synthetic biology technology.

Heavy metal stress and mitogen activated kinase transcription factors in plants: Exploring heavy metal-ROS influences on plant signalling pathways

Due to their stationary nature, plants are exposed to a diverse range of biotic and abiotic stresses, of which heavy metal (HM) stress poses one of the most detrimental abiotic stresses, targeting diverse plant processes. HMs instigate the overproduction of reactive oxygen species (ROS), and to mitigate the adverse effects of ROS, plants induce multiple defence mechanisms. Besides the negative implications of overproduction of ROS, these molecules play a multitude of signalling roles in plants, acting as a central player in the complex signalling network of cells. One of the ROS-associated signalling mechanisms is the mitogen-activated protein kinase (MAPK) cascade, a signalling pathway which transduces extracellular stimuli into intracellular responses. Plant MAPKs have been implicated in signalling involved in stress response, phytohormone regulation, and cell cycle cues. However, the influence of various HMs on MAPK activation has not been well documented. In this review, we address and summarise several aspects related to various HM-induced ROS signalling. Additionally, we touch on how these signals activate the MAPK cascade and the downstream transcription factors that influence plant responses to HMs. Moreover, we propose a workflow that could characterise genes associated with MAPKs and their roles during plant HM stress responses.

Comparative transcriptional analysis of Persea americana MYB, WRKY and AP2/ERF transcription factors following Phytophthora cinnamomi infection

Plant cells undergo extensive transcriptional reprogramming following pathogen infection, with these reprogramming patterns becoming more complex when pathogens, such as hemibiotrophs, exhibit different lifestyles. These transcriptional changes are often orchestrated by MYB, WRKY and AP2/ERF transcription factors (TFs), which modulate both growth and defence-related gene expression. Transcriptional analysis of defence-related genes in avocado (Persea americana) infected with Phytophthora cinnamomi indicated differential immune response activation when comparing a partially resistant and susceptible rootstock. This study identified 226 MYB, 82 WRKY, and 174 AP2/ERF TF-encoding genes in avocado, using a genome-wide approach. Phylogenetic analysis revealed substantial sequence conservation within TF groups underscoring their functional significance. RNA-sequencing analysis in a partially resistant and susceptible avocado rootstock infected with P. cinnamomi was indicative of an immune response switch occurring in either rootstock after 24 and 6 h post-inoculation, respectively. Different clusters of co-expressed TF genes were observed at these times, suggesting the activation of necrotroph-related immune responses at varying intervals between the two rootstocks. This study aids our understanding of avocado immune response activation following P. cinnamomi infection, and the role of the TFs therein, elucidating the transcriptional reprogramming disparities between partially resistant and susceptible rootstocks.

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.

Integrative systems biology of wheat susceptibility to Fusarium graminearum uncovers a conserved gene regulatory network and identifies master regulators targeted by fungal core effectors

BACKGROUND

Plant diseases are driven by an intricate set of defense mechanisms counterbalanced by the expression of host susceptibility factors promoted through the action of pathogen effectors. In spite of their central role in the establishment of the pathology, the primary components of plant susceptibility are still poorly understood and challenging to trace especially in plant-fungal interactions such as in Fusarium head blight (FHB) of bread wheat. Designing a system-level transcriptomics approach, we leveraged the analysis of wheat responses from a susceptible cultivar facing Fusarium graminearum strains of different aggressiveness and examined their constancy in four other wheat cultivars also developing FHB.

RESULTS

In this study, we describe unexpected differential expression of a conserved set of transcription factors and an original subset of master regulators were evidenced using a regulation network approach. The dual-integration with the expression data of pathogen effector genes combined with database mining, demonstrated robust connections with the plant molecular regulators and identified relevant candidate genes involved in plant susceptibility, mostly able to suppress plant defense mechanisms. Furthermore, taking advantage of wheat cultivars of contrasting susceptibility levels, a refined list of 142 conserved susceptibility gene candidates was proposed to be necessary host's determinants for the establishment of a compatible interaction.

CONCLUSIONS

Our findings emphasized major FHB determinants potentially controlling a set of conserved responses associated with susceptibility in bread wheat. They provide new clues for improving FHB control in wheat and also could conceivably leverage further original researches dealing with a broader spectrum of plant pathogens.

Genome-Wide Comprehensive Identification and In Silico Characterization of Lectin Receptor-Like Kinase Gene Family in Barley (Hordeum vulgare L.)

Lectin receptor-like kinases (LecRLKs) are a significant subgroup of the receptor-like kinases (RLKs) protein family. They play crucial roles in plant growth, development, immune responses, signal transduction, and stress tolerance. However, the genome-wide identification and characterization of LecRLK genes and their regulatory elements have not been explored in a major cereal crop, barley (Hordeum vulgare L.). Therefore, in this study, integrated bioinformatics tools were used to identify and characterize the LecRLK gene family in barley. Based on the phylogenetic tree and domain organization, a total of 113 LecRLK genes were identified in the barley genome (referred to as HvlecRLK) corresponding to the LecRLK genes of Arabidopsis thaliana. These putative HvlecRLK genes were classified into three groups: 62 G-type LecRLKs, 1 C-type LecRLK, and 50 L-type LecRLKs. They were unevenly distributed across eight chromosomes, including one unknown chromosome, and were predominantly located in the plasma membrane (G-type HvlecRLK (96.8%), C-type HvlecRLK (100%), and L-type HvlecRLK (98%)). An analysis of motif composition and exon-intron configuration revealed remarkable homogeneity with the members of AtlecRLK. Notably, most of the HvlecRLKs (27 G-type, 43 L-type) have no intron, suggesting their rapid functionality. The Ka/Ks and syntenic analysis demonstrated that HvlecRLK gene pairs evolved through purifying selection and gene duplication was the major factor for the expansion of the HvlecRLK gene family. Exploration of gene ontology (GO) enrichment indicated that the identified HvlecRLK genes are associated with various cellular processes, metabolic pathways, defense mechanisms, kinase activity, catalytic activity, ion binding, and other essential pathways. The regulatory network analysis identified 29 transcription factor families (TFFs), with seven major TFFs including bZIP, C2H2, ERF, MIKC_MADS, MYB, NAC, and WRKY participating in the regulation of HvlecRLK gene functions. Most notably, eight TFFs were found to be linked to the promoter region of both L-type HvleckRLK64 and HvleckRLK86. The promoter cis-acting regulatory element (CARE) analysis of barley identified a total of 75 CARE motifs responsive to light responsiveness (LR), tissue-specific (TS), hormone responsiveness (HR), and stress responsiveness (SR). The maximum number of CAREs was identified in HvleckRLK11 (25 for LR), HvleckRLK69 (17 for TS), and HvleckRLK80 (12 for HR). Additionally, HvleckRLK14, HvleckRLK16, HvleckRLK33, HvleckRLK50, HvleckRLK52, HvleckRLK56, and HvleckRLK110 were predicted to exhibit higher responses in stress conditions. In addition, 46 putative miRNAs were predicted to target 81 HvlecRLK genes and HvlecRLK13 was the most targeted gene by 8 different miRNAs. Protein-protein interaction analysis demonstrated higher functional similarities of 63 HvlecRLKs with 7 Arabidopsis STRING proteins. Our overall findings provide valuable information on the LecRLK gene family which might pave the way to advanced research on the functional mechanism of the candidate genes as well as to develop new barley cultivars in breeding programs.

Physiological and Molecular Response of Liriodendron chinense to Varying Stand Density

Stand density affects the potentially superior productivity of forest ecosystems directly by regulating the light and nutrient availability of trees. Understanding how stand density influences the growth and development of trees is crucial for supporting forest management in the context of climate change. We focused on Liriodendron chinense in experimental plantations created in 2003, with planting densities ranging from 277 to 10,000 trees per hectare at six plots. The leaf structure and photosynthetic capacity of L. chinense changed significantly under different stand densities, which had a negative impact on their biomass (leaf mass) and nutrient (total carbohydrate content) accumulation. Transcriptional differences were observed among samples from plots with different planting densities. The expression of 1784 genes was negatively dependent on stand density, participating mainly in the biological processes of "circadian rhythm", "carbon metabolism", and "amino acid biosynthesis". Furthermore, we identified a photosynthesis-related module and constructed a gene regulatory network to discover that the transcription factors of MYB and bHLH may have important roles in the transcriptional regulation of photosynthesis biosynthesis by activating or repressing the expression of petA (Litul.15G096200), psbE (Litul.10G033900), and petD (Litul.17G061600) at different stand densities. Our study quantified the impact of stand density on tree growth at physiological and molecular levels. Our observations provide theoretical support for plantation establishment of L. chinense.

Overexpression of a 'Beta' MYB Factor Gene, VhMYB15, Increases Salinity and Drought Tolerance in Arabidopsis thaliana

'Beta' is a hybrid of Vitis riparia L. and V. labrusca and has a strong ability to adapt to adverse growth environments and is mainly cultivated and used as a resistant rootstock. At present, the most extensively studied MYB TFs are R2R3-type, which have been found to be involved in plant growth, development, and stress response processes. In the present research, VhMYB15, a key transcription factor for abiotic stress tolerance, was screened by bioinformatics in 'Beta' rootstock, and its function under salinity and drought stresses was investigated. VhMYB15 was highly expressed in roots and mature leave under salinity and drought stresses. Observing the phenotype and calculating the survival rate of plants, it was found that VhMYB15-overexpressing plants exhibited relatively less yellowing and wilting of leaves and a higher survival rate under salinity and drought stresses. Consistent with the above results, through the determination of stress-related physiological indicators and the expression analysis of stress-related genes (AtSOS2, AtSOS3, AtSOS1, AtNHX1, AtSnRK2.6, AtNCED3, AtP5CS1, and AtCAT1), it was found that transgenic Arabidopsis showed better stress tolerance and stronger adaptability under salinity and drought stresses. Based on the above data, it was preliminarily indicated that VhMYB15 may be a key factor in salinity and drought regulation networks, enhancing the adaptability of 'Beta' to adverse environments.

Molecular and metabolic insights into purplish leaf coloration through the investigation of two mulberry (Morus alba) genotypes

BACKGROUND

Leaf coloration in plants, attributed to anthocyanin compounds, plays a crucial role in various physiological functions, and also for pharmaceutical and horticultural uses. However, the molecular mechanisms governing leaf coloration and the physiological significance of anthocyanins in leaves remain poorly understood.

RESULTS

In this study, we investigated leaf color variation in two closely related mulberry genotypes, one with purplish-red young leaves (EP) and another with normal leaf color (EW). We integrated transcriptomic and metabolomic approaches to gain insights into the metabolic and genetic basis of purplish-red leaf development in mulberry. Our results revealed that flavonoid biosynthesis, particularly the accumulation of delphinidin-3-O-glucoside, is a key determinant of leaf color. Additionally, the up-regulation of CHS genes and transcription factors, including MYB family members, likely contributes to the increased flavonoid content in purplish-red leaves.

CONCLUSION

These findings enhance our understanding of the molecular mechanisms responsible for the purplish coloration observed in mulberry leaves and also offer supporting evidence for the hypothesis that anthocyanins serve a protective function in plant tissues until the processes of light absorption and carbon fixation reach maturity, thereby ensuring a balanced equilibrium between energy capture and utilization.

Comparative transcriptome profiling reveals differential defense responses among Alternaria brassicicola resistant Sinapis alba and susceptible Brassica rapa

Alternaria blight is a devastating disease that causes significant crop losses in oilseed Brassicas every year. Adoption of conventional breeding to generate disease-resistant varieties has so far been unsuccessful due to the lack of suitable resistant source germplasms of cultivated Brassica spp. A thorough understanding of the molecular basis of resistance, as well as the identification of defense-related genes involved in resistance responses in closely related wild germplasms, would substantially aid in disease management. In the current study, a comparative transcriptome profiling was performed using Illumina based RNA-seq to detect differentially expressed genes (DEGs) specifically modulated in response to Alternaria brassicicola infection in resistant Sinapis alba, a close relative of Brassicas, and the highly susceptible Brassica rapa. The analysis revealed that, at 48 hpi (hours post inoculation), 3396 genes were upregulated and 23239 were downregulated, whereas at 72 hpi, 4023 genes were upregulated and 21116 were downregulated. Furthermore, a large number of defense response genes were detected to be specifically regulated as a result of Alternaria infection. The transcriptome data was validated using qPCR-based expression profiling for selected defense-related DEGs, that revealed significantly higher fold change in gene expression in S. alba when compared to B. rapa. Expression of most of the selected genes was elevated across all the time points under study with significantly higher expression towards the later time point of 72 hpi in the resistant germplasm. S. alba activates a stronger defense response reaction against the disease by deploying an array of genes and transcription factors involved in a wide range of biological processes such as pathogen recognition, signal transduction, cell wall modification, antioxidation, transcription regulation, etc. Overall, the study provides new insights on resistance of S. alba against A. brassicicola, which will aid in devising strategies for breeding resistant varieties of oilseed Brassica.

Identification of MYB gene family in medicinal tea tree Melaleuca alternifolia (Maiden and Betche) cheel and analysis of members regulating terpene biosynthesis

BACKGROUND

The tea tree (Melaleuca alternifolia) is renowned for its production of tea tree oil, an essential oil primarily composed of terpenes extracted from its shoot. MYB transcription factors, which are one of the largest TF families, play a crucial role in regulating primary and secondary metabolite synthesis. However, knowledge of the MYB gene family in M. alternifolia is limited.

METHODS AND RESULTS

Here, we conducted a comprehensive genome-wide analysis of MYB genes in M. alternifolia, referred to as MaMYBs, including phylogenetic relationships, structures, promoter regions, and GO annotations. Our findings classified 219 MaMYBs into four subfamilies: one 5R-MYB, four 3R-MYBs, sixty-one MYB-related, and the remaining 153 are all 2R-MYBs. Seven genes (MYB189, MYB146, MYB44, MYB29, MYB175, MYB162, and MYB160) were linked to terpenoid synthesis based on GO annotation. Phylogenetic analysis with Arabidopsis homologous MYB genes suggested that MYB193 and MYB163 may also be involved in terpenoid synthesis. Additionally, through correlation analysis of gene expression and metabolite content, we identified 42 MYB genes associated with metabolite content.

CONCLUSION

The results provide valuable insights into the importance of MYB transcription factors in essential oil production in M. alternifolia. These findings lay the groundwork for a better understanding of the MYB regulatory network and the development of novel strategies to enhance essential oil synthesis in M. alternifolia.

Interactions between plant-beneficial microorganisms in a consortium: Streptomyces microflavus and Trichoderma harzianum

The construction of microbial consortia is challenging due to many variables to be controlled, including the cross-compatibility of the selected strains and their additive or synergistic effects on plants. In this work, we investigated the interactions in vitro, in planta, and at the molecular level of two elite biological control agents (BCAs), that is Streptomyces microflavus strain AtB-42 and Trichoderma harzianum strain M10, to understand their attitude to cooperate in a consortium. In vitro, we observed a strong cross-antagonism between AtB-42 and M10 in agar plates due to diffusible metabolites and volatile organic compounds. In liquid co-cultures, M10 hindered the growth of AtB-42 very likely because of secondary metabolites and strong competition for the nutrients. The interaction in the co-culture induced extensive transcriptional reprogramming in both strains, especially in the pathways related to ribosomes, protein synthesis, and oxidoreductase activity, suggesting that each strain recognized the counterpart and activated its defence responses. The metabolome of both strains was also significantly affected. In contrast, in the soil, M10 growth was partially contrasted by AtB-42. The roots of tomato seedlings inoculated with the consortium appeared smaller than the control and single-strain-inoculated plants, indicating that plants diverted some energy from the development to defence activation, as evidenced by the leaf transcriptome. The consortium induced a stronger transcriptional change compared to the single inoculants, as demonstrated by a higher number of differentially expressed genes. Although the cross-antagonism observed in vitro, the two strains exerted a synergistic effect on tomato seedlings by inducing resistance responses stronger than the single inoculants. Our observations pose a question on the usefulness of the sole in vitro assays for selecting BCAs to construct a consortium. In vivo experiments should be preferred, and transcriptomics may greatly help to elucidate the activity of the BCAs beyond the phenotypic effects on the plant.

Stress combinations and their interactions in plants database: a one-stop resource on combined stress responses in plants

We have developed a compendium and interactive platform, named Stress Combinations and their Interactions in Plants Database (SCIPDb; http://www.nipgr.ac.in/scipdb.php), which offers information on morpho-physio-biochemical (phenome) and molecular (transcriptome and metabolome) responses of plants to different stress combinations. SCIPDb is a plant stress informatics hub for data mining on phenome, transcriptome, trait-gene ontology, and data-driven research for advancing mechanistic understanding of combined stress biology. We analyzed global phenome data from 939 studies to delineate the effects of various stress combinations on yield in major crops and found that yield was substantially affected under abiotic-abiotic stresses. Transcriptome datasets from 36 studies hosted in SCIPDb identified novel genes, whose roles have not been earlier established in combined stress. Integretome analysis under combined drought-heat stress pinpointed carbohydrate, amino acid, and energy metabolism pathways as the crucial metabolic, proteomic, and transcriptional components in plant tolerance to combined stress. These examples illustrate the application of SCIPDb in identifying novel genes and pathways involved in combined stress tolerance. Further, we showed the application of this database in identifying novel candidate genes and pathways for combined drought and pathogen stress tolerance. To our knowledge, SCIPDb is the only publicly available platform offering combined stress-specific omics big data visualization tools, such as an interactive scrollbar, stress matrix, radial tree, global distribution map, meta-phenome analysis, search, BLAST, transcript expression pattern table, Manhattan plot, and co-expression network. These tools facilitate a better understanding of the mechanisms underlying plant responses to combined stresses.

Characteristics and Functions of MYB (v-Myb avivan myoblastsis virus oncogene homolog)-Related Genes in Arabidopsis thaliana

The MYB (v-Myb avivan myoblastsis virus oncogene homolog) transcription factor family is one of the largest families of plant transcription factors which plays a vital role in many aspects of plant growth and development. MYB-related is a subclass of the MYB family. Fifty-nine Arabidopsis thaliana MYB-related (AtMYB-related) genes have been identified. In order to understand the functions of these genes, in this review, the promoters of AtMYB-related genes were analyzed by means of bioinformatics, and the progress of research into the functions of these genes has been described. The main functions of these AtMYB-related genes are light response and circadian rhythm regulation, root hair and trichome development, telomere DNA binding, and hormone response. From an analysis of cis-acting elements, it was found that the promoters of these genes contained light-responsive elements and plant hormone response elements. Most genes contained elements related to drought, low temperature, and defense and stress responses. These analyses suggest that AtMYB-related genes may be involved in A. thaliana growth and development, and environmental adaptation through plant hormone pathways. However, the functions of many genes do not occur independently but instead interact with each other through different pathways. In the future, the study of the role of the gene in different pathways will be conducive to a comprehensive understanding of the function of the gene. Therefore, gene cloning and protein functional analyses can be subsequently used to understand the regulatory mechanisms of AtMYB-related genes in the interaction of multiple signal pathways. This review provides theoretical guidance for the follow-up study of plant MYB-related genes.

Polyphenols in Plants: Structure, Biosynthesis, Abiotic Stress Regulation, and Practical Applications (Review)

Phenolic compounds or polyphenols are among the most common compounds of secondary metabolism in plants. Their biosynthesis is characteristic of all plant cells and is carried out with the participation of the shikimate and acetate-malonate pathways. In this case, polyphenols of various structures are formed, such as phenylpropanoids, flavonoids, and various oligomeric and polymeric compounds of phenolic nature. Their number already exceeds 10,000. The diversity of phenolics affects their biological activity and functional role. Most of their representatives are characterized by interaction with reactive oxygen species, which manifests itself not only in plants but also in the human body, where they enter through food chains. Having a high biological activity, phenolic compounds are successfully used as medicines and nutritional supplements for the health of the population. The accumulation and biosynthesis of polyphenols in plants depend on many factors, including physiological-biochemical, molecular-genetic, and environmental factors. In the review, we present the latest literature data on the structure of various classes of phenolic compounds, their antioxidant activity, and their biosynthesis, including their molecular genetic aspects (genes and transfactors). Since plants grow with significant environmental changes on the planet, their response to the action of abiotic factors (light, UV radiation, temperature, and heavy metals) at the level of accumulation and composition of these secondary metabolites, as well as their metabolic regulation, is considered. Information is given about plant polyphenols as important and necessary components of functional nutrition and pharmaceutically valuable substances for the health of the population. Proposals on promising areas of research and development in the field of plant polyphenols are presented.

A genome‑wide approach to the systematic and comprehensive analysis of LIM gene family in sorghum (Sorghum bicolor L.)

The LIM domain-containing proteins are dominantly found in plants and play a significant role in various biological processes such as gene transcription as well as actin cytoskeletal organization. Nevertheless, genome-wide identification as well as functional analysis of the LIM gene family have not yet been reported in the economically important plant sorghum (Sorghum bicolor L.). Therefore, we conducted an in silico identification and characterization of LIM genes in S. bicolor genome using integrated bioinformatics approaches. Based on phylogenetic tree analysis and conserved domain, we identified five LIM genes in S. bicolor (SbLIM) genome corresponding to Arabidopsis LIM (AtLIM) genes. The conserved domain, motif as well as gene structure analyses of the SbLIM gene family showed the similarity within the SbLIM and AtLIM members. The gene ontology (GO) enrichment study revealed that the candidate LIM genes are directly involved in cytoskeletal organization and various other important biological as well as molecular pathways. Some important families of regulating transcription factors such as ERF, MYB, WRKY, NAC, bZIP, C2H2, Dof, and G2-like were detected by analyzing their interaction network with identified SbLIM genes. The cis-acting regulatory elements related to predicted SbLIM genes were identified as responsive to light, hormones, stress, and other functions. The present study will provide valuable useful information about LIM genes in sorghum which would pave the way for the future study of functional pathways of candidate SbLIM genes as well as their regulatory factors in wet-lab experiments.

Genome-Wide Identification and Expression Analysis of the RADIALIS-like Gene Family in Camellia sinensis

The RADIALIS-like (RL) proteins are v-myb avian myeloblastosis viral oncogene homolog (MYB)-related transcription factors (TFs), and are involved in many biological processes, including metabolism, development, and response to biotic and abiotic stresses. However, the studies on the RL genes of Camellia sinensis are not comprehensive enough. Therefore, we undertook this study and identified eight CsaRLs based on the typical conserved domain SANT Associated domain (SANT) of RL. These genes have low molecular weights and theoretical pI values ranging from 5.67 to 9.76. Gene structure analysis revealed that six CsaRL genes comprise two exons and one intron, while the other two contain a single exon encompassing motifs 1 and 2, and part of motif 3. The phylogenetic analysis divided one hundred and fifty-eight RL proteins into five primary classes, in which CsaRLs clustered in Group V and were homologous with CssRLs of the Shuchazao variety. In addition, we selected different tissue parts to analyze the expression profile of CsaRLs, and the results show that almost all genes displayed variable expression levels across tissues, with CsaRL1a relatively abundant in all tissues. qRT-PCR (real-time fluorescence quantitative PCR) was used to detect the relative expression levels of the CsaRL genes under various abiotic stimuli, and it was found that CsaRL1a expression levels were substantially higher than other genes, with abscisic acid (ABA) causing the highest expression. The self-activation assay with yeast two-hybrid system showed that CsaRL1a has no transcriptional activity. According to protein functional interaction networks, CsaRL1a was well connected with WIN1-like, lysine histidine transporter-1-like, β-amylase 3 chloroplastic-like, carbonic anhydrase-2-like (CA2), and carbonic anhydrase dnaJC76 (DJC76). This study adds to our understanding of the RL family and lays the groundwork for further research into the function and regulatory mechanisms of the CsaRLs gene family in Camellia sinensis.

Identification and functional analysis of ZmDLS associated with the response to biotic stress in maize

Terpenes are the main class of secondary metabolites produced in response to pest and germ attacks. In maize (Zea mays L.), they are the essential components of the herbivore-induced plant volatile mixture, which functioned as a direct or indirect defense against pest and germ attacks. In this study, 43 maize terpene synthase gene (ZmTPS) family members were systematically identified and analyzed through the whole genomes of maize. Nine genes, including Zm00001d032230, Zm00001d045054, Zm00001d024486, Zm00001d004279, Zm00001d002351, Zm00001d002350, Zm00001d053916, Zm00001d015053, and Zm00001d015054, were isolated for their differential expression pattern in leaves after corn borer (Ostrinia nubilalis) bite. Additionally, six genes (Zm00001d045054, Zm00001d024486, Zm00001d002351, Zm00001d002350, Zm00001d015053, and Zm00001d015054) were significantly upregulated in response to corn borer bite. Among them, Zm00001d045054 was cloned. Heterologous expression and enzyme activity assays revealed that Zm00001d045054 functioned as d-limonene synthase. It was renamed ZmDLS. Further analysis demonstrated that its expression was upregulated in response to corn borer bites and Fusarium graminearum attacks. The mutant of ZmDLS downregulated the expressions of Zm00001d024486, Zm00001d002351, Zm00001d002350, Zm00001d015053, and Zm00001d015054. It was more attractive to corn borer bites and more susceptible to F. graminearum infection. The yeast one-hybrid assay and dual-luciferase assay showed that ZmMYB76 and ZmMYB101 could upregulate the expression of ZmDLS by binding to the promoter region. This study may provide a theoretical basis for the functional analysis and transcriptional regulation of terpene synthase genes in crops.

Transcriptomic and Metabolomic Investigation on Leaf Necrosis Induced by ZmWus2 Transient Overexpression in Nicotiana benthamiana

WUSCHEL (WUS) is a crucial transcription factor in regulating plant stem cell development, and its expression can also improve genetic transformation. However, the ectopic expression of WUS always causes pleiotropic effects during genetic transformation, making it important to understand the regulatory mechanisms underlying these phenomena. In our study, we found that the transient expression of the maize WUS ortholog ZmWus2 caused severe leaf necrosis in Nicotiana benthamiana. We performed transcriptomic and non-target metabolomic analyses on tobacco leaves during healthy to wilted states after ZmWus2 transient overexpression. Transcriptomic analysis revealed that ZmWus2 transformation caused active metabolism of inositol trisphosphate and glycerol-3-phosphate, while also upregulating plant hormone signaling and downregulating photosystem and protein folding pathways. Metabolomic analysis mainly identified changes in the synthesis of phenylpropanoid compounds and various lipid classes, including steroid synthesis. In addition, transcription factors such as ethylene-responsive factors (ERFs), the basic helix-loop-helix (bHLH) factors, and MYBs were found to be regulated by ZmWus2. By integrating these findings, we developed a WUS regulatory model that includes plant hormone accumulation, stress responses, lipid remodeling, and leaf necrosis. Our study sheds light on the mechanisms underlying WUS ectopic expression causing leaf necrosis and may inform the development of future genetic transformation strategies.

A novel transcriptional repressor complex MYB22-TOPLESS-HDAC1 promotes rice resistance to brown planthopper by repressing F3'H expression

The brown planthopper (BPH) is the most destructive pest of rice. The MYB transcription factors are vital for rice immunity, but most are activators. Although MYB22 positively regulates rice resistance to BPH and has an EAR motif associated with active repression, it remains unclear whether it is a transcriptional repressor affecting rice-BPH interaction. Genetic analyses revealed that MYB22 regulates rice resistance to BPH via its EAR motif. Several biochemical experiments (e.g. transient transcription assay, Y2H, LCA, and BiFC) indicated that MYB22 is a transcriptional repressor that interacts with the corepressor TOPLESS via its EAR motif and recruits HDAC1 to form a tripartite complex. Flavonoid-3'-hydroxylase (F3'H) is a flavonoid biosynthesis pathway-related gene that negatively regulates rice resistance to BPH. Based on a bioinformatics analysis and the results of EMSA and transient transcription assays, MYB22 can bind directly to the F3'H promoter and repress gene expression along with TOPLESS and HDAC1. We revealed a transcriptional regulatory mechanism influencing the rice-BPH interaction that differs from previously reported mechanisms. Specifically, MYB22-TOPLESS-HDAC1 is a novel transcriptional repressor complex with components that synergistically and positively regulate rice resistance to BPH through the transcriptional repression of F3'H.

Multilayered regulation of secondary metabolism in medicinal plants

Medicinal plants represent a huge reservoir of secondary metabolites (SMs), substances with significant pharmaceutical and industrial potential. However, obtaining secondary metabolites remains a challenge due to their low-yield accumulation in medicinal plants; moreover, these secondary metabolites are produced through tightly coordinated pathways involving many spatiotemporally and environmentally regulated steps. The first regulatory layer involves a complex network of transcription factors; a second, more recently discovered layer of complexity in the regulation of SMs is epigenetic modification, such as DNA methylation, histone modification and small RNA-based mechanisms, which can jointly or separately influence secondary metabolites by regulating gene expression. Here, we summarize the findings in the fields of genetic and epigenetic regulation with a special emphasis on SMs in medicinal plants, providing a new perspective on the multiple layers of regulation of gene expression.

PsERF1B-PsMYB10.1-PsbHLH3 module enhances anthocyanin biosynthesis in the flesh-reddening of amber-fleshed plum (cv. Friar) fruit in response to cold storage

Flesh-reddening usually occurs in the amber-fleshed plum (Prunus salicina Lindl.) fruit during cold storage but not during ambient storage direct after harvest. It is not clear how postharvest cold signal is mediated to regulate the anthocyanin biosynthesis in the forming of flesh-reddening yet. In this study, anthocyanins dramatically accumulated and ethylene produced in the 'Friar' plums during cold storage, in comparison with plums directly stored at ambient temperature. Expression of genes associated with anthocyanin biosynthesis, as well as transcription factors of PsMYB10.1, PsbHLH3, and PsERF1B were strongly stimulated to upregulated in the plums in the period of cold storage. Suppression of ethylene act with 1-methylcyclopropene greatly suppressed flesh-reddening and downregulated the expression of these genes. Transient overexpression and virus-induced gene silencing assays in plum flesh indicated that PsMYB10.1 encodes a positive regulator of anthocyanin accumulation. The transient overexpression of PsERF1B, coupled with PsMYB10.1 and PsbHLH3, could further prompt the anthocyanin biosynthesis in a tobacco leaf system. Results from yeast two-hybrid and luciferase complementation assays verified that PsERF1B directly interacted with PsMYB10.1. PsERF1B and PsMYB10.1 enhanced the activity of the promoter of PsUFGT individually, and the enhancement was prompted by the co-action of PsERF1B and PsMYB10.1. Overall, the stimulation of the PsERF1B-PsMYB10.1-PsbHLH3 module mediated cold signal in the transcriptomic supervision of the anthocyanin biosynthesis in the 'Friar' plums. The results thereby revealed the underlying mechanism of the postharvest alteration of the flesh phenotype of 'Friar' plums subjected to low temperature.

Genome-Wide Identification of the MYB and bHLH Families in Carnations and Expression Analysis at Different Floral Development Stages

Carnations are one of the most popular ornamental flowers in the world with varied flower colors that have long attracted breeders and consumers alike. The differences in carnation flower color are mainly the result of the accumulation of flavonoid compounds in the petals. Anthocyanins are a type of flavonoid compound that produce richer colors. The expression of anthocyanin biosynthetic genes is mainly regulated by MYB and bHLH transcription factors. However, these TFs have not been comprehensively reported in popular carnation cultivars. Herein, 106 MYB and 125 bHLH genes were identified in the carnation genome. Gene structure and protein motif analyses show that members of the same subgroup have similar exon/intron and motif organization. Phylogenetic analysis combining the MYB and bHLH TFs from Arabidopsis thaliana separates the carnation DcaMYBs and DcabHLHs into 20 subgroups each. Gene expression (RNAseq) and phylogenetic analysis shows that DcaMYB13 in subgroup S4 and DcabHLH125 in subgroup IIIf have similar expression patterns to those of DFR, ANS, and GT/AT, which regulate anthocyanin accumulation, in the coloring of carnations, and in red-flowered and white-flowered carnations, DcaMYB13 and DcabHLH125 are likely the key genes responsible for the formation of red petals in carnations. These results lay a foundation for the study of MYB and bHLH TFs in carnations and provide valuable information for the functional verification of these genes in studies of tissue-specific regulation of anthocyanin biosynthesis.

Genome-Wide Identification and Expression Profiling of Pathogenesis-Related Protein 1 (PR-1) Genes in Durum Wheat (Triticum durum Desf.)

Pathogen-related proteins (PRs) are diversified proteins with a low molecular weight implicated in plant response to biotic and abiotic stress as well in regulating different functions in plant maturation. Interestingly, no systematical study has been conducted in durum wheat (Triticum turgidum subsp. durum). In the present study, 12 PR-1 genes encoding a CAP superfamily domain were identified in the genome of Triticum turgidum subsp. durum, which is an important cereal, using in silico approaches. Additionally, phylogenetic analysis showed that the PR-1 genes were classified into three groups based on their isoelectric point and the conserved motif domain. Moreover, our analysis showed that most of the TdPR-1 proteins presented an N-terminal signal peptide. Expression patterns analysis showed that the PR-1 gene family presented temporal and spatial specificity and was induced by different abiotic stresses. This is the first report describing the genome-scale analysis of the durum wheat PR-1 gene family, and these data will help further study the roles of PR-1 genes during stress responses, leading to crop improvement.

BcMYB111 Responds to BcCBF2 and Induces Flavonol Biosynthesis to Enhance Tolerance under Cold Stress in Non-Heading Chinese Cabbage

Flavonols have been shown to respond to a variety of abiotic stresses in plants, including cold stress. Higher total flavonoid content was found in non-heading Chinese cabbage (NHCC, Brassica campestris (syn. Brassica rapa) ssp. chinensis) after cold stress. A non-targeted metabolome analysis showed a significant increase in flavonol content, including that of quercetin and kaempferol. Here, we found that an R2R3-MYB transcription factor, BcMYB111, may play a role in this process. BcMYB111 was up-regulated in response to cold treatment, with an accompanying accumulation of flavonols. Then, it was found that BcMYB111 could regulate the synthesis of flavonols by directly binding to the promoters of BcF3H and BcFLS1. In the transgenic hairy roots of NHCC or stable transgenic Arabidopsis, overexpression of BcMYB111 increased flavonol synthesis and accumulation, while these were reduced in virus-induced gene silencing lines in NHCC. After cold stress, the higher proline content and lower malondialdehyde (MDA) content showed that there was less damage in transgenic Arabidopsis than in the wild-type (WT). The BcMYB111 transgenic lines performed better in terms of antioxidant capacity because of their lower H₂O₂ content and higher superoxide dismutase (SOD) and peroxidase (POD) enzyme activities. In addition, a key cold signaling gene, BcCBF2, could specifically bind to the DRE element and activate the expression of BcMYB111 in vitro and in vivo. The results suggested that BcMYB111 played a positive role in enhancing the flavonol synthesis and cold tolerance of NHCC. Taken together, these findings reveal that cold stress induces the accumulation of flavonols to increase tolerance via the pathway of BcCBF2-BcMYB111-BcF3H/BcFLS1 in NHCC.

Metabolome and transcriptome analyses reveal the molecular mechanisms of LcMYB1 regulating anthocyanin accumulation in litchi hairy roots

Agrobacterium rhizogenes-mediated hairy root culture offer a promising approach for gene function analysis and production of plant secondary metabolites. Here, we obtained red litchi hairy roots using A. rhizogenes-mediated LcMYB1 transformation. Using high performance liquid chromatography, the main anthocyanins in the red hairy roots were determined to be cyanidin 3-rutinoside and cyanidin 3-glucoside. A total of 164 metabolites were significantly upregulated or downregulated in the red hairy roots, which were mostly involved in flavone and flavonol pathway, and flavonoid pathway. The transcriptome analysis revealed 472 differentially expressed genes (DEGs). Up-regulated genes were considerably enriched in anthocyanin, flavone and flavonol biosynthesis. Integrative metabolite profiling and transcriptome analyses showed that LcF3'H, LcUFGT1, and LcGST4 were key structural genes in anthocyanin biosynthesis. However, the expression of Cinnamyl-alcohol dehydrogenase (CAD) and Peroxidase (POD) leading to the production of lignin were significantly down-regulated, suggesting flavonoids and lignin compete with each other in the phenylpropanoid pathway. A total of 52 DEGs were identified as transcription factors. Correlation analysis showed that 8 transcription factors were positively correlated with LcUFGT1, and LcGST4, involving in anthocyanin biosynthesis. These findings clarify the molecular mechanisms of LcMYB1 regulating anthocyanin accumulation in litchi hairy roots.

The spatiotemporal profile of Dendrobium huoshanense and functional identification of bHLH genes under exogenous MeJA using comparative transcriptomics and genomics

Introduction

Alkaloids are one of the main medicinal components of Dendrobium species. Dendrobium alkaloids are mainly composed of terpene alkaloids. Jasmonic acid (JA) induce the biosynthesis of such alkaloids, mainly by enhancing the expression of JA-responsive genes to increase plant resistance and increase the content of alkaloids. Many JA-responsive genes are the target genes of bHLH transcription factors (TFs), especially the MYC2 transcription factor.

Methods

In this study, the differentially expressed genes involved in the JA signaling pathway were screened out from Dendrobium huoshanense using comparative transcriptomics approaches, revealing the critical roles of basic helix-loop-helix (bHLH) family, particularly the MYC2 subfamily.

Results and discussion

Microsynteny-based comparative genomics demonstrated that whole genome duplication (WGD) and segmental duplication events drove bHLH genes expansion and functional divergence. Tandem duplication accelerated the generation of bHLH paralogs. Multiple sequence alignments showed that all bHLH proteins included bHLH-zip and ACT-like conserved domains. The MYC2 subfamily had a typical bHLH-MYC_N domain. The phylogenetic tree revealed the classification and putative roles of bHLHs. The analysis of cis-acting elements revealed that promoter of the majority of bHLH genes contain multiple regulatory elements relevant to light response, hormone responses, and abiotic stresses, and the bHLH genes could be activated by binding these elements. The expression profiling and qRT-PCR results indicated that bHLH subgroups IIIe and IIId may have an antagonistic role in JA-mediated expression of stress-related genes. DhbHLH20 and DhbHLH21 were considered to be the positive regulators in the early response of JA signaling, while DhbHLH24 and DhbHLH25 might be the negative regulators. Our findings may provide a practical reference for the functional study of DhbHLH genes and the regulation of secondary metabolites.

The Role of Plant Transcription Factors in the Fight against Plant Viruses

Plants are vulnerable to the challenges of unstable environments and pathogen infections due to their immobility. Among various stress conditions, viral infection is a major threat that causes significant crop loss. In response to viral infection, plants undergo complex molecular and physiological changes, which trigger defense and morphogenic pathways. Transcription factors (TFs), and their interactions with cofactors and cis-regulatory genomic elements, are essential for plant defense mechanisms. The transcriptional regulation by TFs is crucial in establishing plant defense and associated activities during viral infections. Therefore, identifying and characterizing the critical genes involved in the responses of plants against virus stress is essential for the development of transgenic plants that exhibit enhanced tolerance or resistance. This article reviews the current understanding of the transcriptional control of plant defenses, with a special focus on NAC, MYB, WRKY, bZIP, and AP2/ERF TFs. The review provides an update on the latest advances in understanding how plant TFs regulate defense genes expression during viral infection.

Identification of the passion fruit (Passiflora edulis Sims) MYB family in fruit development and abiotic stress, and functional analysis of PeMYB87 in abiotic stresses

Environmental stresses are ubiquitous in agricultural cultivation, and they affect the healthy growth and development of edible tissues in passion fruit. The study of resistance mechanisms is important in understanding the adaptation and resistance of plants to environmental stresses. In this work, two differently resistant passion fruit varieties were selected, using the expression characteristics of the transcription factor MYB, to explore the resistance mechanism of the MYB gene under various environmental stresses. A total of 174 MYB family members were identified using high-quality passion fruit genomes: 98 2R-MYB, 5 3R-MYB, and 71 1R-MYB (MYB-relate). Their family information was systematically analyzed, including subcellular localization, physicochemical properties, phylogeny at the genomic level, promoter function, encoded proteins, and reciprocal regulation. In this study, bioinformatics and transcriptome sequencing were used to identify members of the PeMYB genes in passion fruit whole-genome data, and biological techniques, such as qPCR, gene clone, and transient transformation of yeast, were used to determine the function of the passion fruit MYB genes in abiotic stress tolerance. Transcriptomic data were obtained for differential expression characteristics of two resistant and susceptible varieties, three expression patterns during pulp development, and four induced expression patterns under abiotic stress conditions. We further focused on the resistance mechanism of PeMYB87 in environmental stress, and we selected 10 representative PeMYB genes for quantitative expression verification. Most of the genes were differentially induced by four abiotic stresses, among which PeMYB87 responded significantly to high-temperature-induced expression and overexpression of the PeMYB87 gene in the yeast system. The transgenic PeMYB87 in yeast showed different degrees of stress resistance under exposure to cold, high temperatures, drought, and salt stresses. These findings lay the foundation for further analysis of the biological functions of PeMYBs involved in stress resistance in passion fruit.

Genome-Wide Identification and Analysis of MYB Transcription Factor Family in Hibiscus hamabo

MYB transcription factors constitute one of the largest gene families in plants and play essential roles in the regulation of plant growth, responses to stress, and a wide variety of physiological and biochemical processes. In this study, 204 MYB proteins (HhMYBs) were identified in the Hibiscus hamabo Sieb. et Zucc (H. hamabo) genome and systematically analyzed based on their genomic sequence and transcriptomic data. The candidate HhMYB proteins and MYBs of Arabidopsis thaliana were divided into 28 subfamilies based on the analysis of their phylogenetic relationships and their motif patterns. Expression analysis using RNA-seq and quantitative real-time PCR (qRT-PCR) indicated that most HhMYBs are differentially regulated under drought and salt stresses. qRT-PCR analysis of seven selected HhMYBs suggested that the HhMYB family may have regulatory roles in the responses to stress and hormones. This study provides a framework for a more comprehensive analysis of the role of MYBs in the response to abiotic stress in H. hamabo.

Large-scale analysis of putative Euphorbiaceae R2R3-MYB transcription factors identifies a MYB involved in seed oil biosynthesis

BACKGROUND

MYB transcription factors are widely distributed in the plant kingdom and play key roles in regulatory networks governing plant metabolism and biochemical and physiological processes.

RESULTS

Here, we first determined the R2R3-MYB genes in five Euphorbiaceae genomes. The three Trp (W) residues from the first MYB domain (R2) were absolutely conserved, whereas the first W residue from the second MYB domain (R3) was preferentially mutated. The R2R3-MYBs were clustered into 48 functional subfamilies, of which 34 had both R2R3-MYBs of Euphorbiaceae species and AtMYBs, and four contained only Euphorbiaceae R2R3-MYBs. The whole-genome duplication (WGD) and/or segmental duplication (SD) played key roles in the expansion of the R2R3-MYB family. Unlike paralogous R2R3-MYB family members, orthologous R2R3-MYB members contained a higher selective pressure and were subject to a constrained evolutionary rate. VfMYB36 was specifically expressed in fruit, and its trend was consistent with the change in oil content, indicating that it might be involved in oil biosynthesis. Overexpression experiments showed that VfMYB36 could significantly provide linolenic acid (C18:3) content, which eventually led to a significant increase in oil content.

CONCLUSION

Our study first provides insight into understanding the evolution and expression of R2R3-MYBs in Euphorbiaceae species, and also provides a target for the production of biomass diesel and a convenient way for breeding germplasm resources with high linolenic acid content in the future.

Genome-wide characterization and expression analysis of MYB transcription factors in Chrysanthemum nankingense

BACKGROUND

Chrysanthemum is a popular ornamental plant worldwide. MYB (v-myb avian myeloblastosis viral oncogene homolog) transcription factors play an important role in everything from stress resistance to plant growth and development. However, the MYB family of chrysanthemums has not been the subject of a detailed bioinformatics and expression investigation.

RESULTS

In this study, we examined 324 CnMYB transcription factors from Chrysanthemum nankingense genome data, which contained 122 Cn1R-MYB, 183 CnR2R3-MYB, 12 Cn3R-MYB, 2 Cn4R-MYB, and 5 atypical CnMYB. The protein motifs and classification of CnMYB transcription factors were analyzed. Among them, motifs 1, 2, 3, and 4 were found to encode the MYB DNA-binding domain in R2R3-MYB proteins, while in other-MYB proteins, the motifs 1, 2, 3, 4, 5, 6, 7, and 8 encode the MYB DNA-binding domain. Among all CnMYBs, 44 genes were selected due to the presence of CpG islands, while methylation is detected in three genes, including CnMYB9, CnMYB152, and CnMYB219. We analyzed the expression levels of each CnMYB gene in ray floret, disc floret, flower bud, leaf, stem, and root tissues. Based on phylogenetic analysis and gene expression analysis, three genes appeared likely to control cellulose and lignin synthesis in stem tissue, and 16 genes appeared likely to regulate flowering time, anther, pollen development, and flower color. Fifty-one candidate genes that may be involved in stress response were identified through phylogenetic, stress-responseve motif of promoter, and qRT-PCR analyses. According to genes expression levels under stress conditions, six CnMYB genes (CnMYB9, CnMYB172, CnMYB186, CnMYB199, CnMYB219, and CnMYB152) were identified as key stress-responsive genes.

CONCLUSIONS

This research provides useful information for further functional analysis of the CnMYB gene family in chrysanthemums, as well as offers candidate genes for further study of cellulose and lignin synthesis, flowering traits, salt and drought stress mechanism.

MdMYB305-MdbHLH33-MdMYB10 regulates sugar and anthocyanin balance in red-fleshed apple fruits

Sugar and anthocyanin are important indicators of fruit quality, and understanding the mechanism underlying their accumulation is essential for breeding high-quality fruit. We identified an R2R3-MYB transcription factor MdMYB305 in the red-fleshed apple progeny, which was positively correlated with fruit sugar content but negatively correlated with anthocyanin content. Transient injection, stable expression [overexpressing and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9)], and heterologous transformation of tomato confirmed that MdMYB305 promotes the accumulation of sugar and inhibits the synthesis of anthocyanin. A series of molecular experiments (such as electrophoretic mobility shift and luciferase assays) confirmed that MdMYB305 combines with sugar-related genes (MdCWI1/MdVGT3/MdTMT2) and anthocyanin-related genes (MdF3H/MdDFR/MdUFGT), promoting and inhibiting their activities, and finally regulating the sugar and anthocyanin content of fruits. In addition, the study also found that MdMYB305 competes with MdMYB10 for the MdbHLH33 binding site to balance sugar and anthocyanin accumulation in the fruits, which provides a reference value for exploring more functions of the MYB-bHLH-MYB complex and the balance relationship between sugar and anthocyanin in the future.

An R2R3-MYB network modulates stem strength by regulating lignin biosynthesis and secondary cell wall thickening in herbaceous peony

Stem strength is an important agronomic trait affecting plant lodging, and plays an essential role in the quality and yield of plants. Thickened secondary cell walls in stems provide mechanical strength that allows plants to stand upright, but the regulatory mechanism of secondary cell wall thickening and stem strength in cut flowers remains unclear. In this study, first, a total of 11 non-redundant Paeonia lactiflora R2R3-MYBs related to stem strength were identified and isolated from cut-flower herbaceous peony, among which PlMYB43, PlMYB83 and PlMYB103 were the most upregulated differentially expressed genes. Then, the expression characteristics revealed that these three R2R3-MYBs were specifically expressed in stems and acted as transcriptional activators. Next, biological function verification showed that these P. lactiflora R2R3-MYBs positively regulated stem strength, secondary cell wall thickness and lignin deposition. Furthermore, yeast-one-hybrid and dual luciferase reporter assays demonstrated that they could bind to the promoter of caffeic acid O-methyltransferase gene (PlCOMT2) and/or laccase gene (PlLAC4), two key genes involved in lignin biosynthesis. In addition, the function of PlLAC4 in increasing lignin deposition was confirmed by virus-induced gene silencing and overexpression. Moreover, PlMYB83 could also act as a transcriptional activator of PlMYB43. The findings of the study propose a regulatory network of R2R3-MYBs modulating lignin biosynthesis and secondary cell wall thickening for improving stem lodging resistance, and provide a resource for molecular genetic engineering breeding of cut flowers.

Molecular module of CmMYB15-like-Cm4CL2 regulating lignin biosynthesis of chrysanthemum (Chrysanthemum morifolium) in response to aphid (Macrosiphoniella sanborni) feeding

Lignin is a major component of plant cell walls and a conserved basic defense mechanism in higher plants deposited in response to aphid infection. However, the molecular mechanisms of lignin biosynthesis in response to aphid infection and the effect of lignin on aphid feeding behavior remain unclear. We report that 4-Coumarate:coenzyme A ligase 2 (Cm4CL2), a gene encoding a key enzyme in the lignin biosynthesis pathway, is induced by aphid feeding, resulting in lignin deposition and reduced aphid attack. Upstream regulator analysis showed that the expression of Cm4CL2 in response to aphid feeding was directly upregulated by CmMYB15-like, an SG2-type R2R3-MYB transcription factor. CmMYB15-like binds directly to the AC cis-element in the promoter region of Cm4CL2. Genetic validation demonstrated that CmMYB15-like was induced by aphid infection and contributed to lignin deposition and cell wall thickening, which consequently enhanced aphid resistance in a Cm4CL2-dependent manner. This study is the first to show that the CmMYB15-like-Cm4CL2 module regulates lignin biosynthesis in response to aphid feeding.

Transcriptome and metabolome analyses reveal novel insights into the seed germination of Michelia chapensis, an endangered species in China

Michelia chapensis Dandy, a well-known medicinal woody plant endemic to China, is endangered and seriously constricted by seed dormancy-induced low-regeneration in natural conditions. Cold stratification can effectively reduce seed dormancy and promote the seed germination of M. chapensis. However, the molecular events and systematic changes that occurred during seed germination in M. chapensis remain largely unknown. In this study, we carried out transcriptomic and metabolomic analyses to elucidate the potential molecular mechanisms underlying seed germination in M. chapensis under cold stratification. The results showed that the embryo cells became bigger and looser with increasing stratification time. Moreover, the endosperm appeared reduced due to the consumption of nutrients. Seventeen phytohormones were examined by the metabolome targeted for hormones. Compared with the ES (no stratification), the levels of indole-3-acetic acid (IAA) and gibberellin A₃ (GA₃) were increased in the MS (stratification for 45 days), while the abscisic acid (ABA) was downregulated in both MS and LS (stratification for 90 days). The transcriptome profiling identified 24975 differentially expressed genes (DEGs) in the seeds during germination. The seed germination of M. chapensis was mainly regulated by the biological pathways of plant hormone signal transduction, energy supply, secondary metabolite biosynthesis, photosynthesis-related metabolism, and transcriptional regulation. This study reveals the biological evidence of seed germination at the transcriptional level and provides a foundation for unraveling molecular mechanisms regulating the seed germination of M. chapensis.

A MYB-related transcription factor from peanut, AhMYB30, improves freezing and salt stress tolerance in transgenic Arabidopsis through both DREB/CBF and ABA-signaling pathways

Abiotic stresses such as salinity and low temperature have serious impact on peanut growth and yield. The present work investigated the function of a MYB-related transcription factor gene AhMYB30 obtained from peanut under salt and low temperature stresses by transgenic methods. The results indicated that the overexpression of AhMYB30 in Arabidopsis could enhance the resistance of transgenic plants to freezing and salt stresses. The expression of stress-response genes RD29A (Response-to-Dehydration 29A), COR15A (Cold-Regulated 15A), KIN1 (Kinesin 1) and ABI2 (Abscisic acid Insensitive 2) increased in transgenic plants compared with in wild-type. Subcellular localization and transcriptional autoactivation validation demonstrated that AhMYB30 has essential features of transcription factors. Therefore, AhMYB30 may increase salt and freezing stress tolerance as the transcription factor (TF) in Arabidopsis through both DREB/CBF and ABA-signaling pathways. Our results lay the theoretical foundation for exploring stress resistance mechanisms of peanut and offering novel genetic resources for molecular breeding.

Genome-wide analysis of R2R3-MYB genes in cultivated peanut (Arachis hypogaea L.): Gene duplications, functional conservation, and diversification

The cultivated Peanut (Arachis hypogaea L.), an important oilseed and edible legume, are widely grown worldwide. The R2R3-MYB transcription factor, one of the largest gene families in plants, is involved in various plant developmental processes and responds to multiple stresses. In this study we identified 196 typical R2R3-MYB genes in the genome of cultivated peanut. Comparative phylogenetic analysis with Arabidopsis divided them into 48 subgroups. The motif composition and gene structure independently supported the subgroup delineation. Collinearity analysis indicated polyploidization, tandem, and segmental duplication were the main driver of the R2R3-MYB gene amplification in peanut. Homologous gene pairs between the two subgroups showed tissue specific biased expression. In addition, a total of 90 R2R3-MYB genes showed significant differential expression levels in response to waterlogging stress. Furthermore, we identified an SNP located in the third exon region of AdMYB03-18 (AhMYB033) by association analysis, and the three haplotypes of the SNP were significantly correlated with total branch number (TBN), pod length (PL) and root-shoot ratio (RS ratio), respectively, revealing the potential function of AdMYB03-18 (AhMYB033) in improving peanut yield. Together, these studies provide evidence for functional diversity in the R2R3-MYB genes and will contribute to understanding the function of R2R3-MYB genes in peanut.

SbMYB3 transcription factor promotes root-specific flavone biosynthesis in Scutellaria baicalensis

Scutellaria baicalensis Georgi produces abundant root-specific flavones (RSFs), which provide various benefits to human health. We have elucidated the complete biosynthetic pathways of baicalein and wogonin. However, the transcriptional regulation of flavone biosynthesis in S. baicalensis remains unclear. We show that the SbMYB3 transcription factor functions as a transcriptional activator involved in the biosynthesis of RSFs in S. baicalensis. Yeast one-hybrid and transcriptional activation assays showed that SbMYB3 binds to the promoter of flavone synthase II-2 (SbFNSII-2) and enhances its transcription. In S. baicalensis hairy roots, RNAi of SbMYB3 reduced the accumulation of baicalin and wogonoside, and SbMYB3 knockout decreased the biosynthesis of baicalein, baicalin, wogonin, and wogonoside, whereas SbMYB3 overexpression enhanced the contents of baicalein, baicalin, wogonin, and wogonoside. Transcript profiling by qRT-PCR demonstrated that SbMYB3 activates SbFNSII-2 expression directly, thus leading to more abundant accumulation of RSFs. This study provides a potential target for metabolic engineering of RSFs.

Genome-wide characterization and functional identification of MYB genes in Malus sieversii infected by Valsa mali

Among the most important transcription factors in plants, the v-myb avian myeloblastosis viral oncogene homolog (MYB) regulates the expression network of response genes under stresses such as fungal infection. In China, the canker disease Valsa mali threatens the survival of Malus sieversii, an ancestor of cultivated apples. Using the M. sieversii genome, we identified 457 MsMYB and 128 R2R3-MsMYB genes that were randomly distributed across 17 chromosomes. Based on protein sequence and structure, the R2R3-MsMYB genes were phylogenetically divided into 29 categories, and 26 conserved motifs were identified. We further predicted cis-elements in the 2000-kb promoter region of R2R3-MsMYBs based on the genome. Transcriptome analysis of M. sieversii under V. mali infection showed that 27 R2R3-MsMYBs were significantly differentially expressed, indicating their key role in the response to V. mali infection. Using transient transformation, MsMYB14, MsMYB24, MsMYB39, MsMYB78, and MsMYB108, which were strongly induced by V. mali infection, were functionally identified. Among the five MsMYBs, MsMYB14 and MsMYB78 were both important in enhancing resistance to diseases, whereas MsMYB24 inhibited resistance. Based on the results of this study, we gained a better understanding of the MsMYB transcription factor family and laid the foundation for a future research program on disease prevention strategies in M. sieversii.

Overexpression of GmPHR1 Promotes Soybean Yield through Global Regulation of Nutrient Acquisition and Root Development

MYB-CC transcription factors (TFs) are essential for plant growth and development. Members of the MYB-CC subfamily with long N terminal domains, such as phosphate starvation response 1 (PHR1) or PHR1-like TFs, have well documented functions, while those with short N terminal domains remain less understood. In this study, we identified a nodule specific MYB-CC transcription factor 1 (GmPHR1) in soybean that is different from other canonical PHR family genes in that GmPHR1 harbors a short N terminal ahead of its MYB-CC domain and was highly induced by rhizobium infection. The overexpression of GmPHR1 dramatically increased the ratio of deformed root hairs, enhanced subsequent soybean nodulation, and promoted soybean growth in pot experiments. The growth promotion effects of GmPHR1 overexpression were further demonstrated in field trails in which two GmPHR1-OE lines yielded 10.78% and 8.19% more than the wild type line. Transcriptome analysis suggested that GmPHR1 overexpression led to global reprogramming, with 749 genes upregulated and 279 genes downregulated, especially for genes involved in MYB transcription factor activities, root growth, and nutrient acquisition. Taken together, we conclude that GmPHR1 is a key gene involved in the global regulation of nodulation, root growth, and nutrient acquisition in soybeans, and is thus a promising candidate gene to target for soybean yield enhancement.

MYB transcription factors-master regulators of phenylpropanoid biosynthesis and diverse developmental and stress responses

Phenylpropanoids, the largest class of natural products including flavonoids, anthocyanins, monolignols and tannins perform multiple functions ranging from photosynthesis, nutrient uptake, regulating growth, cell division, maintenance of redox homeostasis and biotic and abiotic stress responses. Being sedentary life forms, plants possess several regulatory modules that increase their performance in varying environments by facilitating activation of several signaling cascades upon perception of developmental and stress signals. Of the various regulatory modules, those involving MYB transcription factors are one of the extensive groups involved in regulating the phenylpropanoid metabolic enzymes in addition to other genes. R2R3 MYB transcription factors are a class of plant-specific transcription factors that regulate the expression of structural genes involved in anthocyanin, flavonoid and monolignol biosynthesis which are indispensable to several developmental pathways and stress responses. The aim of this review is to present the regulation of the phenylpropanoid pathway by MYB transcription factors via Phospholipase D/phosphatidic acid signaling, downstream activation of the structural genes, leading to developmental and/or stress responses. Specific MYB transcription factors inducing or repressing specific structural genes of anthocyanin, flavonoid and lignin biosynthetic pathways are discussed. Further the roles of MYB in activating biotic and abiotic stress responses are delineated. While several articles have reported the role of MYB's in stress responses, they are restricted to two or three specific MYB factors. This review is a consolidation of the diverse roles of different MYB transcription factors involved both in induction and repression of anthocyanin, flavonoid, and lignin biosynthesis.

Comparative transcriptome profiling and weighted gene co-expression network analysis to identify core genes in maize (Zea mays L.) silks infected by multiple fungi

Maize (Zea mays L.) is the third most popular Poaceae crop after wheat and rice and used in feed and pharmaceutical sectors. The maize silk contains bioactive components explored by traditional Chinese herbal medicine for various pharmacological activities. However, Fusarium graminearum, Fusarium verticillioides, Trichoderma atroviride, and Ustilago maydis can infect the maize, produce mycotoxins, hamper the quantity and quality of silk production, and further harm the primary consumer's health. However, the defense mechanism is not fully understood in multiple fungal infections in the silk of Z. mays. In this study, we applied bioinformatics approaches to use the publicly available transcriptome data of Z. mays silk affected by multiple fungal flora to identify core genes involved in combatting disease response. Differentially expressed genes (DEGs) were identified among intra- and inter-transcriptome data sets of control versus infected Z. mays silks. Upon further comparison between up- and downregulated genes within the control of datasets, 4,519 upregulated and 5,125 downregulated genes were found. The DEGs have been compared with genes in the modules of weighted gene co-expression network analysis to relevant specific traits towards identifying core genes. The expression pattern of transcription factors, carbohydrate-active enzymes (CAZyme), and resistance genes was analyzed. The present investigation is supportive of our findings that the gene ontology, immunity stimulus, and resistance genes are upregulated, but physical and metabolic processes such as cell wall organizations and pectin synthesis were downregulated respectively. Our results are indicative that terpene synthase TPS6 and TPS11 are involved in the defense mechanism against fungal infections in maize silk.

High-resolution genome-wide association study of a large Czech collection of sweet cherry (Prunus avium L.) on fruit maturity and quality traits

In sweet cherry (Prunus avium L.), quantitative trait loci have been identified for fruit maturity, colour, firmness, and size to develop markers for marker-assisted selection. However, resolution is usually too low in those analyses to directly target candidate genes, and some associations are missed. In contrast, genome-wide association studies are performed on broad collections of accessions, and assemblies of reference sequences from Tieton and Satonishiki cultivars enable identification of single nucleotide polymorphisms after whole-genome sequencing, providing high marker density. Two hundred and thirty-five sweet cherry accessions were sequenced and phenotyped for harvest time and fruit colour, firmness, and size. Genome-wide association studies were used to identify single nucleotide polymorphisms associated with each trait, which were verified in breeding material consisting of 64 additional accessions. A total of 1 767 106 single nucleotide polymorphisms were identified. At that density, significant single nucleotide polymorphisms could be linked to co-inherited haplotype blocks (median size ~10 kb). Thus, markers were tightly associated with respective phenotypes, and individual allelic combinations of particular single nucleotide polymorphisms provided links to distinct phenotypes. In addition, yellow-fruit accessions were sequenced, and a ~ 90-kb-deletion on chromosome 3 that included five MYB10 transcription factors was associated with the phenotype. Overall, the study confirmed numerous quantitative trait loci from bi-parental populations using high-diversity accession populations, identified novel associations, and genome-wide association studies reduced the size of trait-associated loci from megabases to kilobases and to a few candidate genes per locus. Thus, a framework is provided to develop molecular markers and evaluate and characterize genes underlying important agronomic traits.

Interplay of transcription factors orchestrating the biosynthesis of plant alkaloids

Plants produce a range of secondary metabolites primarily as defence molecules. A plant has to invest considerable energy to synthesise alkaloids, and sometimes they are even toxic to themselves. Hence, the biosynthesis of alkaloids is a spatiotemporally regulated process under quantitative feedback regulation which is accomplished by the signal reception, transcriptional/translational regulation, transport, storage and accumulation. The transcription factors (TFs) initiate the biosynthesis of alkaloids after appropriate cues. The present study recapitulates last decade understanding of the role of TFs in alkaloid biosynthesis. The present review discusses TF families, viz. AP2/ERF, bHLH, WRKY, MYB involved in the biosynthesis of various types of alkaloids. It also highlights the role of the jasmonic acid cascade and post-translational modifications of TF proteins. A thorough understanding of TFs will help us to decide a strategy to facilitate successful pathway manipulation and in vitro production.

Dynamic transcriptome and network-based analysis of yellow leaf mutant Ginkgo biloba

BACKGROUND

Golden leaf in autumn is a prominent feature of deciduous tree species like Ginkgo biloba L., a landscape tree widely cultivated worldwide. However, little was known about the molecular mechanisms of leaf yellowing, especially its dynamic regulatory network. Here, we performed a suite of comparative physiological and dynamic transcriptional analyses on the golden-leaf cultivar and the wild type (WT) ginkgo to investigate the underlying mechanisms of leaf yellowing across different seasons.

RESULTS

In the present study, we used the natural bud mutant cultivar with yellow leaves "Wannianjin" (YL) as materials. Physiological analysis revealed that higher ratios of chlorophyll a to chlorophyll b and carotenoid to chlorophyll b caused the leaf yellowing of YL. On the other hand, dynamic transcriptome analyses showed that genes related to chlorophyll metabolism played key a role in leaf coloration. Genes encoding non-yellow coloring 1 (NYC1), NYC1-like (NOL), and chlorophyllase (CLH) involved in the degradation of chlorophyll were up-regulated in spring. At the summer stage, down-regulated HEMA encoding glutamyl-tRNA reductase functioned in chlorophyll biosynthesis, while CLH involved in chlorophyll degradation was up-regulated, causing a lower chlorophyll accumulation. In carotenoid metabolism, genes encoding zeaxanthin epoxidase (ZEP) and 9-cis-epoxy carotenoid dioxygenase (NCED) showed significantly different expression levels in the WT and YL. Moreover, the weighted gene co-expression network analysis (WGCNA) suggested that the most associated transcriptional factor, which belongs to the AP2/ERF-ERF family, was engaged in regulating pigment metabolism. Furthermore, quantitative experiments validated the above results.

CONCLUSIONS

By comparing the golden-leaf cultivar and the wide type of ginkgo across three seasons, this study not only confirm the vital role of chlorophyll in leaf coloration of YL but also provided new insights into the seasonal transcriptome landscape and co-expression network. Our novel results pinpoint candidate genes for further wet-bench experiments in tree species.

Genome-wide identification and analysis of anthocyanin synthesis-related R2R3-MYB genes in Cymbidium goeringii

The MYB gene family plays a vital regulatory role in plant metabolism, stress response, and floral color. The R2R3-MYB gene family of C. goeringii was identified, and its expression was analyzed using bioinformatics in this article. The R2R3-MYB genes of Arabidopsis thaliana were used as a reference to determine 104 CgMYB genes and categorize them into 22 subfamilies. Exon/intron organizations and conserved motif analysis revealed that the majority of CgMYB genes were highly conserved, and chromosome localization and collinearity analysis provided evidence of tandem duplication and segmental duplication events, indicating the phenomenon of gene family expansion and contraction. The function of CgMYB genes was analyzed by cis-acting element and gene ontology (GO) enrichment. In addition, we selected CgMYB91 and CgMYB32 for RT-qPCR, suggesting that CgMYB91 and CgMYB32 are associated with anthocyanin formation. In short, this study provides a comprehensive and specific function of the R2R3-MYB transcription factors (TFs) in orchids.