The R2R3-MYB transcription factor gene family in maize

MYB4 in Lilium pumilum affects plant saline-alkaline tolerance

Lilium pumilum DC (L. pumilum DC) plays an important role in the rational utilization of salinized soil. To explore the molecular mechanism of salt-tolerant L. pumilum, the LpMYB4 was cloned. LpMYB4 close relationship with Bambusa emeiensis and Zea mays MYB4 throughout the phylogenetic tree construction. LpMYB4 protein was found to be localized in the nucleus. Prokaryotic and eukaryotic bacterial solution resistance experiments proved that the exogenous introduction of LpMYB4 made the overexpression strains obtain better survival ability under saline-alkaline stress. Compared with wild-type plants, tobacco plants overexpressing LpMYB4 had better growth and lower leaf wilting and lodging, the content of chlorophyll was higher, the content of hydrogen peroxide and superoxide anion was lower, the activity of peroxidase and superoxide dismutase was higher and the relative conductivity was lower under saline-alkaline stress. The analysis of seed germination and seedling resistance of transgenic plants under salt stress showed that LpMYB4 transgenic seeds were more tolerant to salt stress during germination and growth. Yeast two-hybrid and two-luciferase complementation experiments showed that LpMYB4 interacted with yeast two-hybrid and LpGPX6. The analysis of the role of LpMYB4 in improving plant saline-alkali resistance is helpful to the transformation of plant germplasm resources and has great significance for agriculture and sustainable development.

Genome-wide identification of R2R3-MYB transcription factors in Betula platyphylla and functional analysis of BpMYB95 in salt tolerance

The Myeloblastosis (MYB) transcription factor (TF) family is one of the largest transcription factor families in plants and plays an important role in various physiological processes. At present, there are few reports on birch (Betula platyphylla Suk.) of R2R3-MYB-TFs, and most BpMYBs still need to be characterized. In this study, 111 R2R3-MYB-TFs with conserved R2 and R3 MYB domains were identified. Phylogenetic tree analysis showed that the MYB family members of Arabidopsis thaliana and birch were divided into 23 and 21 subgroups, respectively. The latter exhibited an uneven distribution across 14 chromosomes. There were five tandem duplication events and 17 segmental duplication events between BpMYBs, and repeat events play an important role in the expansion of the family. In addition, the promoter region of MYBs was rich in various cis-acting elements, and MYB-TFs were involved in plant growth and development, light responses, biotic stress, and abiotic stress. RNA-sequencing (RNA-seq) and quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) results revealed that most R2R3-MYB-TFs in birch responded to salt stress. In particular, the expression of BpMYBs in the S20 subfamily was significantly induced by salt, drought, abscisic acid, and methyl jasmonate stresses. Based on the weighted co-expression network analysis of physiological and RNA-seq data of birch under salt stress, a key MYB-TF BpMYB95 (BPChr12G24087), was identified in response to salt stress, and its expression level was induced by salt stress. BpMYB95 is a nuclear localization protein with transcriptional activation activity in yeast and overexpression of this gene significantly enhanced salt tolerance in Saccharomyces cerevisiae. The qRT-PCR and histochemical staining results showed that BpMYB95 exhibited the highest expression in the roots, young leaves, and petioles of birch plants. Overexpression of BpMYB95 significantly improved salt-induced browning and wilting symptoms in birch leaves and alleviated the degree of PSII photoinhibition caused by salt stress in birch seedlings. In conclusion, most R2R3-MYB-TFs found in birch were involved in the salt stress response mechanisms. Among these, BpMYB95 was a key regulatory factor that significantly enhanced salt tolerance in birch. The findings of this study provide valuable genetic resources for the development of salt-tolerant birch varieties.

Genome-Wide Identification, Characterization, and Expression Analysis of the MYB-R2R3 Gene Family in Black Pepper (Piper nigrum L.)

Black pepper (Piper nigrum L.), a prominent spice crop, known as the "king of spices", originated from India. The growth and development of black pepper are influenced by various environmental conditions. MYB transcription factors play a crucial role in controlling metabolic processes, abiotic stress management, and plant growth and development. In this study, we identified 160 PnMYB transcription factors in the black pepper genome. Phylogenetic analysis was performed using 125 R2R3-MYB proteins from black pepper and Arabidopsis thaliana, resulting in the mapping of 20 groups on the phylogenetic tree, each containing members from both species. Most members of the PnMYB family possess two introns, and motif 3 and motif 4 are conserved in all members. The number of genes on each chromosome ranges from 1 to 10. Collinear analysis indicated the creation of new members through gene fragments and tandem replication. The Ka/Ks ratio indicated that purifying selection and positive selection acted on PnMYB of pepper. The majority of pepper PnMYB family members were in the nucleus. Significant differences in gene expression levels were observed between different species and infection periods when Piper nigrum L. and Piper flaviflorum were infected with Phytophthora capsici. These findings are valuable for future studies on the biological role and molecular mechanism of the PnMYB gene.

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.

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.

Genome-Wide Analysis of Transcription Factor R2R3-MYB Gene Family and Gene Expression Profiles during Anthocyanin Synthesis in Common Walnut (Juglans regia L.)

The R2R3-MYB gene family, encoding plant transcriptional regulators, participates in many metabolic pathways of plant physiology and development, including flavonoid metabolism and anthocyanin synthesis. This study proceeded as follows: the JrR2R3-MYB gene family was analyzed genome-wide, and the family members were identified and characterized using the high-quality walnut reference genome "Chandler 2.0". All 204 JrR2R3-MYBs were established and categorized into 30 subgroups via phylogenetic analysis. JrR2R3-MYBs were unevenly distributed over 16 chromosomes. Most JrR2R3-MYBs had similar structures and conservative motifs. The cis-acting elements exhibit multiple functions of JrR2R3-MYBs such as light response, metabolite response, and stress response. We found that the expansion of JrR2R3-MYBs was mainly caused by WGD or segmental duplication events. Ka/Ks analysis indicated that these genes were in a state of negative purifying selection. Transcriptome results suggested that JrR2R3-MYBs were widely entangled in the process of walnut organ development and differentially expressed in different colored varieties of walnuts. Subsequently, we identified 17 differentially expressed JrR2R3-MYBs, 9 of which may regulate anthocyanin biosynthesis based on the results of a phylogenetic analysis. These genes were present in greater expression levels in 'Zijing' leaves than in 'Lvling' leaves, as revealed by the results of qRT-PCR experiments. These results contributed to the elucidation of the functions of JrR2R3-MYBs in walnut coloration. Collectively, this work provides a foundation for exploring the functional characteristics of the JrR2R3-MYBs in walnuts and improving the nutritional value and appearance quality of walnuts.

An ARF gene mutation creates flint kernel architecture in dent maize

Dent and flint kernel architectures are important characteristics that affect the physical properties of maize kernels and their grain end uses. The genes controlling these traits are unknown, so it is difficult to combine the advantageous kernel traits of both. We found mutation of ARFTF17 in a dent genetic background reduces IAA content in the seed pericarp, creating a flint-like kernel phenotype. ARFTF17 is highly expressed in the pericarp and encodes a protein that interacts with and inhibits MYB40, a transcription factor with the dual functions of repressing PIN1 expression and transactivating genes for flavonoid biosynthesis. Enhanced flavonoid biosynthesis could reduce the metabolic flux responsible for auxin biosynthesis. The decreased IAA content of the dent pericarp appears to reduce cell division and expansion, creating a shorter, denser kernel. Introgression of the ARFTF17 mutation into dent inbreds and hybrids improved their kernel texture, integrity, and desiccation, without affecting yield.

Genome-wide identification of R2R3-MYB family genes and gene response to stress in ginger

Ginger (Zingiber officinale Roscoe) is an important plant used worldwide for medicine and food. The R2R3-MYB transcription factor (TF) family has essential roles in plant growth, development, and stresses resistance, and the number of genes in the family varies greatly among different types of plants. However, genome-wide discovery of ZoMYBs and gene responses to stresses have not been reported in ginger. Therefore, genome-wide analysis of R2R3-MYB genes in ginger was conducted in this study. Protein phylogenetic relations and conserved motifs and chromosome localization and duplication, structure, and cis-regulatory elements were analyzed. In addition, the expression patterns of selected genes were analyzed under two different stresses. A total of 299 candidate ZoMYB genes were discovered in ginger. Based on groupings of R2R3-MYB genes in the model plant Arabidopsis thaliana (L.) Heynh., ZoMYBs were divided into eight groups. Genes were distributed across 22 chromosomes at uneven densities. In gene duplication analysis, 120 segmental duplications were identified in the ginger genome. Gene expression patterns of 10 ZoMYBs in leaves of ginger under abscisic acid (ABA) and low-temperature stress treatments were different. The results will help to determine the exact roles of ZoMYBs in anti-stress responses in ginger.

Mechanistic and structural insight into R2R3-MYB transcription factor in plants: molecular dynamics based binding free energy analysis

The MYB transcription factor (TF) family is essential for various plant growth and development processes, including responses to biotic and abiotic stresses. This study investigated the R2R3-MYB protein structure from five plants, including cereal crops. The R2R3-MYB protein structure was docked with the DNA structure, and the best complexes were selected for two runs of molecular dynamics (MD) simulations to investigate the key interacting residues and conformational changes in the R2R3-MYB proteins caused by DNA binding. The MM/PBSA method calculated the binding free energy for each R2R3-MYB protein-DNA complex, showing strong interaction. Hydrophobic and hydrogen bonds significantly stabilized the R2R3-MYB protein-DNA complexes. The principal component analysis showed high restrictions on the movement of protein atoms in the phase space. A similar MD simulation analysis was performed using the crystal structure of the R2R3-MYB protein-DNA complex from Arabidopsis thaliana, and the generated complexes resembled the X-ray crystal structure. This is the first detailed study on the R2R3-MYB protein-DNA complex in cereal crops, providing a cost-effective solution to identify the key interacting residues and analyze the conformational changes in the MYB domain before and after DNA binding.Communicated by Ramaswamy H. Sarma.

The R2R3 MYB gene TaMYB305 positively regulates anther and pollen development in thermo-sensitive male-sterility wheat with Aegilops kotschyi cytoplasm

MAIN CONCLUSION

The loss of TaMYB305 function down-regulated the expression of jasmonic acid synthesis pathway genes, which may disturb the jasmonic acid synthesis, resulting in abnormal pollen development and reduced fertility. The MYB family, as one of the largest transcription factor families found in plants, regulates plant development, especially the development of anthers. Therefore, it is important to identify potential MYB transcription factors associated with pollen development and to study its role in pollen development. Here, the transcripts of an R2R3 MYB gene TaMYB305 from KTM3315A, a thermo-sensitive cytoplasmic male-sterility line with Aegilops kotschyi cytoplasm (K-TCMS) wheat, was isolated. Quantitative real-time PCR (qRT-PCR) and promoter activity analysis revealed that TaMYB305 was primarily expressed in anthers. The TaMYB305 protein was localized in the nucleus, as determined by subcellular localization analysis. Our data demonstrated that silencing of TaMYB305 was related to abnormal development of stamen, including anther indehiscence and pollen abortion in KAM3315A plants. In addition, TaMYB305-silenced plants exhibited alterations in the transcriptional levels of genes involved in the synthesis of jasmonic acid (JA), indicating that TaMYB305 may regulate the expression of genes related to JA synthesis and play an important role during anther and pollen development of KTM3315A. These results provide novel insight into the function and molecular mechanism of R2R3-MYB genes in pollen development.

Genome-wide identification of the Pyrus R2R3-MYB gene family and PhMYB62 regulation analysis in Pyrus hopeiensis flowers at low temperature

The R2R3-MYB gene family play an important role in plant growth, development and stress responses. In this study, a total of 122 PcoR2R3-MYB genes were identified and grouped into 26 clades in pear. And these PcoMYBs were unevenly distributed among 17 chromosomes. The sequence characteristics, conversed motifs, exon/intron structures, classification, duplication events and cis-acting elements were also investigated. The gene duplication events showed that segmental duplication may play key roles in expansion of the PcoMYB gene family. Pyrus hopeiensis, which is a valuable wild resource, has strong cold resistance. An integrative analyses of miRNA and mRNA showed that PhMYB62 was involved in regulating low-temperature stress in P. hopeiensis flower organs. Subcellular localization analysis showed that PhMYB62 protein was specifically localized to the nucleus. The result of DAP-seq showed that PhMYB62 responded to low-temperature stress in P. hopeiensis by regulating TFs, which were associated with plant stress resistance, and POD, GAUT12, AUX28 and CHS genes. Subsequently, yeast one-hybrid verified that PhMYB62 could bind and activate the promoter of POD gene. The current study would provide a comprehensive information for further functional research on the stress-responsive R2R3-MYB gene candidates in pear, and may help to identify the genes associated with cold resistance for the cultivation of cold-resistant pear varieties.

Genome-Wide Analysis of Cotton MYB Transcription Factors and the Functional Validation of GhMYB in Response to Drought Stress

MYB transcription factors play important roles during abiotic stress responses in plants. However, little is known about the accurate systematic analysis of MYB genes in the four cotton species, Gossypium hirsutum, G. barbadense, G. arboreum and G. raimondii. Herein, we performed phylogenetic analysis and showed that cotton MYBs and Arabidopsis MYBs were clustered in the same subfamilies for each species. The identified cotton MYBs were distributed unevenly on chromosomes in various densities for each species, wherein genome-wide tandem and segment duplications were the main driving force of MYB family expansion. Synteny analysis suggested that the abundant collinearity pairs of MYBs were identified between G. hirsutum and the other three species, and that they might have undergone strong purification selection. Characteristics of conserved motifs, along with their consensus sequence, promoter cis elements and gene structure, revealed that MYB proteins might be highly conserved in the same subgroups for each species. Subsequent analysis of differentially expressed genes and expression patterns indicated that most GhMYBs might be involved in response to drought (especially) and salt stress, which was supported by the expression levels of nine GhMYBs using real-time quantitative PCR. Finally, we performed a workflow that combined virus-induced gene silencing and the heterologous transformation of Arabidopsis, which confirmed the positive roles of GhMYBs under drought conditions, as validated by determining the drought-tolerant phenotypes, damage index and/or water loss rate. Collectively, our findings not only expand our understanding of the relationships between evolution and function of MYB genes, but they also provide candidate genes for cotton breeding.

Genome-Wide Identification and Expression Analysis of the MYB Transcription Factor Family in Salvia nemorosa

The MYB transcription factor gene family is among the most extensive superfamilies of transcription factors in plants and is involved in various essential functions, such as plant growth, defense, and pigment formation. Salvia nemorosa is a perennial herb belonging to the Lamiaceae family, and S. nemorosa has various colors and high ornamental value. However, there is little known about its genome-wide MYB gene family and response to flower color formation. In this study, 142 SnMYB genes (MYB genes of S. nemorosa) were totally identified, and phylogenetic relationships, conserved motifs, gene structures, and expression profiles during flower development stages were analyzed. A phylogenetic analysis indicated that MYB proteins in S. nemorosa could be categorized into 24 subgroups, as supported by the conserved motif compositions and gene structures. Furthermore, according to their similarity with AtMYB genes associated with the control of anthocyanin production, ten SnMYB genes related to anthocyanin biosynthesis were speculated and chosen for further qRT-PCR analyses. The results indicated that five SnMYB genes (SnMYB75, SnMYB90, SnMYB6, SnMYB82, and SnMYB12) were expressed significantly differently in flower development stages. In conclusion, our study establishes the groundwork for understanding the anthocyanin biosynthesis of the SnMYB gene family and has the potential to enhance the breeding of S. nemorosa.

Genome-Wide Analysis of MYB Genes in Primulina eburnea (Hance) and Identification of Members in Response to Drought Stress

Due to periodic water deficiency in karst environments, Primulina eburnea experiences sporadic drought stress in its habitat. Despite being one of the largest gene families and functionally diverse in terms of plant growth and development, MYB transcription factors in P. eburnea have not been studied. Here, a total of 230 MYB genes were identified in P. eburnea, including 67 1R-MYB, 155 R2R3-MYB, six 3R-MYB, and two 4R-MYB genes. The R2R3-type PebMYB genes could be classified into 16 subgroups, while the remaining PebMYB genes (1R-MYB, 3R-MYB, and 4R-MYB genes) were divided into 10 subgroups. Notably, the results of the phylogenetic analysis were further supported by the motif and gene structure analysis, which showed that individuals in the same subgroup had comparable motif and structure organization. Additionally, gene duplication and synteny analyses were performed to better understand the evolution of PebMYB genes, and 291 pairs of segmental duplicated genes were found. Moreover, RNA-seq analysis revealed that the PebMYB genes could be divided into five groups based on their expression characteristics. Furthermore, 11 PebMYB genes that may be involved in drought stress response were identified through comparative analysis with Arabidopsis thaliana. Notably, seven of these genes (PebMYB3, PebMYB13, PebMYB17, PebMYB51, PebMYB142, PebMYB69, and PebMYB95) exhibited significant differences in expression between the control and drought stress treatments, suggesting that they may play important roles in drought stress response. These findings clarified the characteristics of the MYB gene family in P. eburnea, augmenting our comprehension of their potential roles in drought stress adaptation.

Petunia as a model for MYB transcription factor action under salt stress

Salinity is a current and growing problem, affecting crops worldwide by reducing yields and product quality. Plants have different mechanisms to adapt to salinity; some crops are highly studied, and their salinity tolerance mechanisms are widely known. However, there are other crops with commercial importance that still need characterization of their molecular mechanisms. Usually, transcription factors are in charge of the regulation of complex processes such as the response to salinity. MYB-TFs are a family of transcription factors that regulate various processes in plant development, and both central and specialized metabolism. MYB-TFs have been studied extensively as mediators of specialized metabolism, and some are master regulators. The influence of MYB-TFs on highly orchestrated mechanisms, such as salinity tolerance, is an attractive research target. The versatility of petunia as a model species has allowed for advances to be made in multiple fields: metabolomic pathways, quality traits, stress resistance, and signal transduction. It has the potential to be the link between horticultural crops and lab models, making it useful in translating discoveries related to the MYB-TF pathways into other crops. We present a phylogenetic tree made with Petunia axillaris and Petunia inflata R2R3-MYB subfamily sequences, which could be used to find functional conservation between different species. This work could set the foundations to improve salinity resistance in other commercial crops in later studies.

Maize response to endophytic Metarhizium robertsii is altered by water stress

To defend against damage from environmental stress, plants have evolved strategies to respond to stress efficiently. One such strategy includes forming mutualist relationships with endophytes which confer stress-alleviating plant defensive and growth promoting effects. Metarhizium robertsii is an entomopathogen and plant-protective and growth-promoting endophyte. To determine the context dependency of the relationship between M. robertsii and maize, we conducted a greenhouse experiment that imposed stress as deficit and excess soil moisture on maize plants which were inoculated or not inoculated with M. robertsii and measured plant growth and defense indicators. Maize height and endophytic root colonization by M. robertsii were positively correlated in the deficit water treatment, but not in the adequate or excess water treatments. The relative expression of ZmLOX1 in the jasmonic acid (JA) biosynthesis pathway was significantly greater in M. robertsii-inoculated than in non-inoculated plants, but water treatment had no effect. There was significant interaction between M. robertsii and water treatments on foliar concentrations of JA and jasmonoyl isoleucine (JA-ILE), suggesting that water stress impacts M. robertsii as a modulator of plant defense. Water stress, but not inoculation with M. robertsii, had a significant effect on the expression of MYB (p = 0.021) and foliar concentrations of abscisic acid (p<0.001), two signaling molecules associated with abiotic stress response. This study contributes toward understanding the highly sophisticated stress response signaling network and context dependency of endophytic mutualisms in crops.

Genome-Wide Identification and Expression Pattern of MYB Family Transcription Factors in Erianthus fulvus

MYB family genes have many functions and are widely involved in plant abiotic-stress responses. Erianthus fulvus is an important donor material for stress-resistance genes in sugarcane breeding. However, the MYB family genes in E. fulvus have not been systematically investigated. In this study, 133 EfMYB genes, including 48 Ef1R-MYB, 84 EfR2R3-MYB and 1 Ef3R-MYB genes, were identified in the E. fulvus genome. Among them, the EfR2R3-MYB genes were classified into 20 subgroups. In addition, these EfMYB genes were unevenly distributed across 10 chromosomes. A total of 4 pairs of tandemly duplicated EfMYB genes and 21 pairs of segmentally duplicated EfMYB genes were identified in the E. fulvus genome. Protein-interaction analysis predicted that 24 EfMYB proteins had potential interactions with 14 other family proteins. The EfMYB promoter mainly contains cis-acting elements related to the hormone response, stress response, and light response. Expression analysis showed that EfMYB39, EfMYB84, and EfMYB124 could be significantly induced using low-temperature stress. EfMYB30, EfMYB70, EfMYB81, and EfMYB101 responded positively to drought stress. ABA treatment significantly induced EfMYB1, EfMYB30, EfMYB39, EfMYB84, and EfMYB130. All nine genes were induced using MeJA treatment. These results provide comprehensive information on EfMYB genes and can serve as a reference for further studies of gene function.

Molecular mechanism analysis of ZmRL6 positively regulating drought stress tolerance in maize

MYB-related genes, a subclass of MYB transcription factor family, have been documented to play important roles in biological processes such as secondary metabolism and stress responses that affect plant growth and development. However, the regulatory roles of MYB-related genes in drought stress response remain unclear in maize. In this study, we discovered that a 1R-MYB gene, ZmRL6, encodes a 96-amino acid protein and is highly drought-inducible. We also found that it is conserved in both barley (Hordeum vulgare L.) and Aegilops tauschii. Furthermore, we observed that overexpression of ZmRL6 can enhance drought tolerance while knock-out of ZmRL6 by CRISPR-Cas9 results in drought hypersensitivity. DAP-seq analyses additionally revealed the ZmRL6 target genes mainly contain ACCGTT, TTACCAAAC and AGCCCGAG motifs in their promoters. By combining RNA-seq and DAP-seq results together, we subsequently identified eight novel target genes of ZmRL6 that are involved in maize's hormone signal transduction, sugar metabolism, lignin synthesis, and redox signaling/oxidative stress. Collectively, our data provided insights into the roles of ZmRL6 in maize's drought response.

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

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

Genomic-wide identification and expression analysis of R2R3-MYB transcription factors related to flavonol biosynthesis in Morinda officinalis

BACKGROUND

The R2R3-MYB transcription factors are a crucial and extensive gene family in plants, which participate in diverse processes, including development, metabolism, defense, differentiation, and stress response. In the Lingnan region of China, Morinda officinalis is extensively grown and is renowned for its use as both a medicinal herb and food source. However, there are relatively few reports on the R2R3-MYB transcription factor family in M.officinalis.

RESULTS

In this study, we identified 97 R2R3-MYB genes in the genome of Morinda officinalis and classified them into 32 subgroups based on phylogenetic comparison with Arabidopsis thaliana. The lack of recent whole-genome duplication events in M.officinalis may be the reason for the relatively few members of the R2R3-MYB family. We also further analyzed the physical and chemical characteristics, conserved motifs, gene structure, and chromosomal location. Gene duplication events found 21 fragment duplication pairs and five tandem duplication event R2R3-MYB genes in M.officinalis may also affect gene family expansion. Based on phylogenetic analysis, cis-element analysis, co-expression analysis and RT-qPCR, we concluded that MoMYB33 might modulate flavonol levels by regulating the expression of 4-coumarate-CoA ligase Mo4CL2, chalcone isomerase MoCHI3, and flavonol synthase MoFLS4/11/12. MoMYB33 and AtMYB111 showed the highest similarity of 79% and may be involved in flavonol synthase networks by the STRING database. Moreover, we also identified MoMYB genes that respond to methyl Jasmonate (MeJA) and abscisic acid (ABA) stress by RT-qPCR.

CONCLUSIONS

This study offers a thorough comprehension of R2R3-MYB in M.officinalis, which lays the foundation for the regulation of flavonol synthesis and the response of MoMYB genes to phytohormones in M.officinalis.

FaMYB5 Interacts with FaBBX24 to Regulate Anthocyanin and Proanthocyanidin Biosynthesis in Strawberry (Fragaria × ananassa)

MYB and BBX transcription factors play important roles in flavonoid biosynthesis. Here, we obtained transgenic woodland strawberry with stable overexpression of FaMYB5, demonstrating that FaMYB5 can increase anthocyanin and proanthocyanidin content in roots, stems and leaves of woodland strawberry. In addition, bimolecular fluorescence complementation assays and yeast two-hybridization demonstrated that the N-terminal (1-99aa) of FaBBX24 interacts with FaMYB5. Transient co-expression of FaBBX24 and FaMYB5 in cultivated strawberry 'Xiaobai' showed that co-expression strongly promoted the expression of F3'H, 4CL-2, TT12, AHA10 and ANR and then increased the content of anthocyanin and proanthocyanidin in strawberry fruits. We also determined that FaBBX24 is also a positive regulator of anthocyanin and proanthocyanidin biosynthesis in strawberry. The results reveal a novel mechanism by which the FaMYB5-FaBBX24 module collaboratively regulates anthocyanin and proanthocyanidin in strawberry fruit.

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.

An unconventional proanthocyanidin pathway in maize

Proanthocyanidins (PAs), flavonoid polymers involved in plant defense, are also beneficial to human health and ruminant nutrition. To date, there is little evidence for accumulation of PAs in maize (Zea mays), although maize makes anthocyanins and possesses the key enzyme of the PA pathway, anthocyanidin reductase (ANR). Here, we explore whether there is a functional PA biosynthesis pathway in maize using a combination of analytical chemistry and genetic approaches. The endogenous PA biosynthetic machinery in maize preferentially produces the unusual PA precursor (+)-epicatechin, as well as 4β-(S-cysteinyl)-catechin, as potential PA starter and extension units. Uncommon procyanidin dimers with (+)-epicatechin as starter unit are also found. Expression of soybean (Glycine max) anthocyanidin reductase 1 (ANR1) in maize seeds increases the levels of 4β-(S-cysteinyl)-epicatechin and procyanidin dimers mainly using (-)-epicatechin as starter units. Introducing a Sorghum bicolor transcription factor (SbTT2) specifically regulating PA biosynthesis into a maize inbred deficient in anthocyanin biosynthesis activates both anthocyanin and PA biosynthesis pathways, suggesting conservation of the PA regulatory machinery across species. Our data support the divergence of PA biosynthesis across plant species and offer perspectives for future agricultrural applications in maize.

Systematic investigation of the R2R3-MYB gene family in Aquilaria sinensis reveals a transcriptional repressor AsMYB054 involved in 2-(2-phenylethyl)chromone biosynthesis

Trees in the genus Aquilaria produce agarwood, a valuable resin used in medicine, perfumes, and incense. 2-(2-Phenethyl)chromones (PECs) are characteristic components of agarwood; however, molecular mechanisms underlying PEC biosynthesis and regulation remain largely unknown. The R2R3-MYB transcription factors play important regulatory roles in the biosynthesis of various secondary metabolites. In this study, 101 R2R3-MYB genes in Aquilaria sinensis were systematically identified and analyzed at the genome-wide level. Transcriptomic analysis revealed that 19 R2R3-MYB genes were significantly regulated by an agarwood inducer, and showed significant correlations with PEC accumulation. Expression and evolutionary analyses revealed that AsMYB054, a subgroup 4 R2R3-MYB, was negatively correlated with PEC accumulation. AsMYB054 was located in the nucleus and functioned as a transcriptional repressor. Moreover, AsMYB054 could bind to the promoters of the PEC biosynthesis related genes AsPKS02 and AsPKS09, and inhibit their transcriptional activity. These findings suggested that AsMYB054 functions as a negative regulator of PEC biosynthesis via the inhibition of AsPKS02 and AsPKS09 in A. sinensis. Our results provide a comprehensive understanding of the R2R3-MYB subfamily in A. sinensis and lay a foundation for further functional analyses of R2R3-MYB genes in PEC biosynthesis.

Genomic survey of MYB gene family in six pearl millet (Pennisetum glaucum) varieties and their response to abiotic stresses

In addition to their roles in developmental and metabolic processes, MYB transcription factors play crucial roles in plant defense mechanisms and stress responses. A comprehensive analysis of six pearl millet genomes revealed the presence of 1133 MYB genes, which can be classified into four phylogenetically distinct subgroups. The duplication pattern of MYB genes across the pearl millet genomes demonstrates their conserved and similar evolutionary history. Overall, MYB genes were observed to be involved in drought and heat stress responses, with stronger differential expressed observed in root tissues. Multiple analyses indicated that MYB genes mediate abiotic stress responses by modulating abscisic acid-related pathways, circadian rhythms, and histone modification processes. A substantial number of duplicated genes were determined to exhibit differential expression under abiotic stress. The consistent positive expression trend observed in duplicated gene pairs, such as PMA5G04432.1 and PMA2G00728.1, across various abiotic stresses suggests that duplicated MYB genes plays a key role in the evolution of adaptive responses of pearl millet to abiotic stresses.

Genome-Wide Identification and Analysis of R2R3-MYB Genes Response to Saline-Alkali Stress in Quinoa

Soil saline-alkalization inhibits plant growth and development and seriously affects crop yields. Over their long-term evolution, plants have formed complex stress response systems to maintain species continuity. R2R3-MYB transcription factors are one of the largest transcription factor families in plants, widely involved in plant growth and development, metabolism, and stress response. Quinoa (Chenopodium quinoa Willd.), as a crop with high nutritional value, is tolerant to various biotic and abiotic stress. In this study, we identified 65 R2R3-MYB genes in quinoa, which are divided into 26 subfamilies. In addition, we analyzed the evolutionary relationships, protein physicochemical properties, conserved domains and motifs, gene structure, and cis-regulatory elements of CqR2R3-MYB family members. To investigate the roles of CqR2R3-MYB transcription factors in abiotic stress response, we performed transcriptome analysis to figure out the expression file of CqR2R3-MYB genes under saline-alkali stress. The results indicate that the expression of the six CqMYB2R genes was altered significantly in quinoa leaves that had undergone saline-alkali stress. Subcellular localization and transcriptional activation activity analysis revealed that CqMYB2R09, CqMYB2R16, CqMYB2R25, and CqMYB2R62, whose Arabidopsis homologues are involved in salt stress response, are localized in the nucleus and exhibit transcriptional activation activity. Our study provides basic information and effective clues for further functional investigation of CqR2R3-MYB transcription factors in quinoa.

Comparative analysis of the MYB gene family in seven Ipomoea species

The MYB transcription factors regulate plant growth, development, and defense responses. However, information about the MYB gene family in Ipomoea species is rare. Herein, we performed a comprehensive genome-wide comparative analysis of this gene family among seven Ipomoea species, sweet potato (I. batatas), I. trifida, I. triloba, I. nil, I. purpurea, I. cairica, and I. aquatic, and identified 296, 430, 411, 291, 226, 281, and 277 MYB genes, respectively. The identified MYB genes were classified into five types: 1R-MYB (MYB-related), 2R-MYB (R2R3-MYB), 3R-MYB (R1R2R3-MYB), 4R-MYB, and 5R-MYB, and the MYB-related or R2R3-MYB type was the most abundant MYB genes in the seven species. The Ipomoea MYB genes were classed into distinct subgroups based on the phylogenetic topology and the classification of the MYB superfamily in Arabidopsis. Analysis of gene structure and protein motifs revealed that members within the same phylogenetic group presented similar exon/intron and motif organization. The identified MYB genes were unevenly mapped on the chromosomes of each Ipomoea species. Duplication analysis indicated that segmental and tandem duplications contribute to expanding the Ipomoea MYB genes. Non-synonymous substitution (Ka) to synonymous substitution (Ks) [Ka/Ks] analysis showed that the duplicated Ipomoea MYB genes are mainly under purifying selection. Numerous cis-regulatory elements related to stress responses were detected in the MYB promoters. Six sweet potato transcriptome datasets referring to abiotic and biotic stresses were analyzed, and MYB different expression genes' (DEGs') responses to stress treatments were detected. Moreover, 10 sweet potato MYB DEGs were selected for qRT-PCR analysis. The results revealed that four responded to biotic stress (stem nematodes and Ceratocystis fimbriata pathogen infection) and six responded to the biotic stress (cold, drought, and salt). The results may provide new insights into the evolution of MYB genes in the Ipomoea genome and contribute to the future molecular breeding of sweet potatoes.

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.

Overexpression of a Fragaria vesca 1R-MYB Transcription Factor Gene (FvMYB114) Increases Salt and Cold Tolerance in Arabidopsis thaliana

The MYB (v-MYB avian myeloblastosis viral oncogene homolog) transcription factor (TF) family has numerous members with complex and diverse functions, which play an indispensable role in regulating the response of plants to stress. In this study, a new 1R-MYB TF gene was obtained from Fragaria vesca (a diploid strawberry) by cloning technology and given a new name, FvMYB114. According to the subcellular localization results, FvMYB114 protein was a nuclear localization protein. Overexpression of FvMYB114 greatly enhanced the adaptability and tolerance of Arabidopsis thaliana to salt and low temperature. Under salt and cold stress, the transgenic plants had greater proline and chlorophyll contents and higher activities of superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT) than the wild-type (WT) and unloaded-line (UL) A. thaliana. However, malondialdehyde (MDA) was higher in the WT and UL lines. These results suggested that FvMYB114 may be involved in regulating the response of A. thaliana to salt stress and cold stress. FvMYB114 can also promote the expression of genes, such as the genes AtSOS1/3, AtNHX1 and AtLEA3 related to salt stress and the genes AtCCA1, AtCOR4 and AtCBF1/3 related to cold stress, further improving the tolerance of transgenic plants to salt and cold stress.

CgMYB1, an R2R3-MYB transcription factor, can alleviate abiotic stress in an annual halophyte Chenopodium glaucum

MYB transcription factors (TFs) are important regulators of the stress response in plants. In the present study, we characterized the CgMYB1 gene in Chenopodium glaucum, a member of the R2R3-MYB TF family. CgMYB1 was located in the nucleus with an activating domain at the C terminus. The CgMYB1 gene could be induced by salt and cold stress in C. glaucum. Overexpressing CgMYB1 in Arabidopsis significantly enhanced salt and cold tolerance, probably by improving physiological performance and stress-related gene expression. Further analysis suggests that the positive response of CgMYB1 to abiotic stress may partially be attributed to the interaction between CgMYB1 and the CgbHLH001 promoter followed by activation of downstream stress-responsive genes, which mediates stress tolerance. Our findings should contribute to further understanding of the function of R2R3 MYB TF in response to abiotic stress.

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 identification and expression pattern analysis of R2R3-MYB transcription factor gene family involved in puerarin biosynthesis and response to hormone in Pueraria lobata var. thomsonii

BACKGROUND

R2R3-MYB transcription factors regulate secondary metabolism, stress responses and development in various plants. Puerarin is a bioactive ingredient and most abundant secondary metabolite isolated from Pueraria lobata. The biosynthesis of puerarin proceeds via the phenylpropanoid pathway and isoflavonoids pathway, in which 9 key enzymes are involved. The expression of these structural genes is under control of specific PtR2R3-MYB genes in different plant tissues. However, how PtR2R3-MYB genes regulates structural genes in puerarin biosynthesis remains elusive. This study mined the PtR2R3-MYB genes involved in puerarin biosynthesis and response to hormone in Pueraria lobata var. thomsonii.

RESULTS

A total of 209 PtR2R3-MYB proteins were identified, in which classified into 34 subgroups based on the phylogenetic topology and the classification of the R2R3-MYB superfamily in Arabidopsis thaliana. Furtherly physical and chemical characteristics, gene structure, and conserved motif analysis were also used to further analyze PtR2R3-MYBs. Combining puerarin content and RNA-seq data, speculated on the regulated puerarin biosynthesis of PtR2R3-MYB genes and structural genes, thus 21 PtR2R3-MYB genes and 25 structural genes were selected for validation gene expression and further explore its response to MeJA and GSH treatment by using qRT-PCR analysis technique. Correlation analysis and cis-acting element analysis revealed that 6 PtR2R3-MYB genes (PtMYB039, PtMYB057, PtMYB080, PtMYB109, PtMYB115 and PtMYB138) and 7 structural genes (PtHID2, PtHID9, PtIFS3, PtUGT069, PtUGT188, PtUGT286 and PtUGT297) were directly or indirectly regulation of puerarin biosynthesis in ZG11. It is worth noting that after MeJA and GSH treatment for 12-24 h, the expression changes of most candidate genes were consistent with the correlation of puerarin biosynthesis, which also shows that MeJA and GSH have the potential to mediate puerarin biosynthesis by regulating gene expression in ZG11.

CONCLUSIONS

Overall, this study provides a comprehensive understanding of the PtR2R3-MYB and will paves the way to reveal the transcriptional regulation of puerarin biosynthesis and response to phytohormone of PtR2R3-MYB genes in Pueraria lobata var. thomsonii.

Genome-wide analysis of MYB family in Nicotiana benthamiana and the functional role of the key members in resistance to Bemisia tabaci

MYB transcription factors (TFs) play a key role in plant resistance to abiotic and biotical stresses. However, little is currently known about their involvement in the plant defense to piercing-sucking insects. Here, we studied the MYB TFs that responded to and resisted Bemisia tabaci whitefly in the model plant Nicotiana benthamiana. Firstly, a total of 453 NbMYB TFs in N. benthamiana genome were identified and 182 R2R3-MYB TFs were analyzed for molecular characteristics, phylogenetic analysis, genetic structure, motif composition, and cis-elements. Then, six stress-related NbMYB genes were selected for further study. The expression pattern shows they were highly expressed in mature leaves and intensively induced upon whitefly attack. Combined with bioinformatic analysis, overexpression, β-Glucuronidase (GUS) assay, and virus-induced silencing tests, we determined the transcriptional regulation of these NbMYBs on the genes in lignin biosynthesis and SA-signaling pathways. Meanwhile, we tested the performance of whitefly on plants with increased or silenced NbMYB genes expression and found that NbMYB42, NbMYB107, NbMYB163, and NbMYB423 were resistant to whitefly. Our results contribute to a comprehensive understanding of the MYB TFs in N. benthamiana. Furthermore, our findings will facilitate further studies on the role of MYB TFs in the interaction between plants and piercing-sucking insects.

Identification and Phylogenetic Analysis of the R2R3-MYB Subfamily in Brassica napus

The R2R3-MYB sub-family proteins are composed of most members of MYB (v-Myb avian myeloblastosis viral oncogene homolog) protein, a plant-specific transcription factor (TF) that is classified into four classes depending on the number of MYB repeats. R2R3-MYB TFs are involved in physiological and biochemical processes. However, the functions of the Brassica napus R2R3-MYB genes are still mainly unknown. In this study, 35 Brassica napus MYB (BnaMYB) genes were screened in the genome of Brassica napus, and details about their physical and chemical characteristics, evolutionary relationships, chromosome locations, gene structures, three-dimensional protein structures, cis-acting promoter elements, and gene duplications were uncovered. The BnaMYB genes have undergone segmental duplications and positive selection pressure, according to evolutionary studies. The same subfamilies have similar intron-exon patterns and motifs, according to the genes' structure and conserved motifs. Additionally, through cis-element analysis, many drought-responsive and other stress-responsive cis-elements have been found in the promoter regions of the BnaMYB genes. The expression of the BnaMYB gene displays a variety of tissue-specific patterns. Ten lignin-related genes were chosen for drought treatment. Our research screened four genes that showed significant upregulation under drought stress, and thus may be important drought-responsive genes. The findings lay a new foundation for understanding the complex mechanisms of BnaMYB in multiple developmental stages and pathways related to drought stress in rapeseed.

MbMYBC1, a M. baccata MYB transcription factor, contribute to cold and drought stress tolerance in transgenic Arabidopsis

Cold and drought stress considerably suppress the development of plants. In this study, a new MYB (v-myb avian myeloblastosis viral)TF gene, MbMYBC1, was isolated from the M. baccata and located in nucleus. MbMYBC1 has a positive response to low temperature and drought stress. After being introduced into Arabidopsis thaliana, the physiological indicators of transgenic Arabidopsis had corresponding changes under these two stresses, the activities of catalase (CAT), peroxidase (POD) and superoxide dismutase (SOD) increased, electrolyte leakage rate (EL) and the content of proline increased, but the content of chlorophyll decreased. In addition, its overexpression can also activate the downstream expression of AtDREB1A, AtCOR15a, AtERD10B and AtCOR47 related to cold stress and AtSnRK2.4, AtRD29A, AtSOD1and AtP5CS1 related to drought stress. Based on these results, we speculate that MbMYBC1 can respond to cold and hydropenia signals, and can be used in transgenic technology to improve plant tolerance to low temperature and drought stress.

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.

Genome-Wide Identification and Characterization of R2R3-MYB Provide Insight into Anthocyanin Biosynthesis Regulation Mechanism of Ananas comosus var. bracteatus

The R2R3-MYB proteins comprise the largest class of MYB transcription factors, which play an essential role in regulating anthocyanin synthesis in various plant species. Ananas comosus var. bracteatus is an important colorful anthocyanins-rich garden plant. The spatio-temporal accumulation of anthocyanins in chimeric leaves, bracts, flowers, and peels makes it an important plant with a long ornamental period and highly improves its commercial value. We conducted a comprehensive bioinformatic analysis of the R2R3-MYB gene family based on genome data from A. comosus var. bracteatus. Phylogenetic analysis, gene structure and motif analysis, gene duplication, collinearity, and promoter analysis were used to analyze the characteristics of this gene family. In this work, a total of 99 R2R3-MYB genes were identified and classified into 33 subfamilies according to phylogenetic analysis, and most of them were localized in the nucleus. We found these genes were mapped to 25 chromosomes. Gene structure and protein motifs were conserved among AbR2R3-MYB genes, especially within the same subfamily. Collinearity analysis revealed four pairs of tandem duplicated genes and 32 segmental duplicates in AbR2R3-MYB genes, indicating that segmental duplication contributed to the amplification of the AbR2R3-MYB gene family. A total of 273 ABRE responsiveness, 66 TCA elements, 97 CGTCA motifs, and TGACG motifs were the main cis elements in the promoter region under response to ABA, SA, and MEJA. These results revealed the potential function of AbR2R3-MYB genes in response to hormone stress. Ten R2R3-MYBs were found to have high homology to MYB proteins reported to be involved in anthocyanin biosynthesis from other plants. RT-qPCR results revealed the 10 AbR2R3-MYB genes showed tissue-specific expression patterns, six of them expressed the highest in the flower, two genes in the bract, and two genes in the leaf. These results suggested that these genes may be the candidates that regulate anthocyanin biosynthesis of A. comosus var. bracteatus in the flower, leaf, and bract, respectively. In addition, the expressions of these 10 AbR2R3-MYB genes were differentially induced by ABA, MEJA, and SA, implying that these genes may play crucial roles in hormone-induced anthocyanin biosynthesis. Our study provided a comprehensive and systematic analysis of AbR2R3-MYB genes and identified the AbR2R3-MYB genes regulating the spatial-temporal anthocyanin biosynthesis in A. comosus var. bracteatus, which would be valuable for further study on the anthocyanin regulation mechanism of A. comosus var. bracteatus.

MYB Transcription Factor Family in Pearl Millet: Genome-Wide Identification, Evolutionary Progression and Expression Analysis under Abiotic Stress and Phytohormone Treatments

Transcription factors (TFs) are the regulatory proteins that act as molecular switches in controlling stress-responsive gene expression. Among them, the MYB transcription factor family is one of the largest TF family in plants, playing a significant role in plant growth, development, phytohormone signaling and stress-responsive processes. Pearl millet (Pennisetum glaucum L.) is one of the most important C4 crop plants of the arid and semi-arid regions of Africa and Southeast Asia for sustaining food and fodder production. To explore the evolutionary mechanism and functional diversity of the MYB family in pearl millet, we conducted a comprehensive genome-wide survey and identified 279 MYB TFs (PgMYB) in pearl millet, distributed unevenly across seven chromosomes of pearl millet. A phylogenetic analysis of the identified PgMYBs classified them into 18 subgroups, and members of the same group showed a similar gene structure and conserved motif/s pattern. Further, duplication events were identified in pearl millet that indicated towards evolutionary progression and expansion of the MYB family. Transcriptome data and relative expression analysis by qRT-PCR identified differentially expressed candidate PgMYBs (PgMYB2, PgMYB9, PgMYB88 and PgMYB151) under dehydration, salinity, heat stress and phytohormone (ABA, SA and MeJA) treatment. Taken together, this study provides valuable information for a prospective functional characterization of the MYB family members of pearl millet and their application in the genetic improvement of crop plants.

Genome-wide analysis of R2R3-MYB transcription factors in Boehmeria nivea (L.) gaudich revealed potential cadmium tolerance and anthocyanin biosynthesis genes

Gene family, especially MYB as one of the largest transcription factor family in plants, the study of its subfunctional characteristics is a key step in the study of plant gene function. The sequencing of ramie genome provides a good opportunity to study the organization and evolutionary characters of the ramie MYB gene at the whole genome level. In this study, a total of 105 BnGR2R3-MYB genes were identified from ramie genome and subsequently grouped into 35 subfamilies according to phylogeny divergence and sequences similarity. Chromosomal localization, gene structure, synteny analysis, gene duplication, promoter analysis, molecular characteristics and subcellular localization were accomplished using several bioinformatics tools. Collinearity analysis showed that the segmental and tandem duplication events is the dominant form of the gene family expansion, and duplications prominent in distal telomeric regions. Highest syntenic relationship was obtained between BnGR2R3-MYB genes and that of Apocynum venetum (88). Furthermore, transcriptomic data and phylogenetic analysis revealed that BnGMYB60, BnGMYB79/80 and BnGMYB70 might inhibit the biosynthesis of anthocyanins, and UPLC-QTOF-MS data further supported the results. qPCR and phylogenetic analysis revealed that the six genes (BnGMYB9, BnGMYB10, BnGMYB12, BnGMYB28, BnGMYB41, and BnGMYB78) were cadmium stress responsive genes. Especially, the expression of BnGMYB10/12/41 in roots, stems and leaves all increased more than 10-fold after cadmium stress, and in addition they may interact with key genes regulating flavonoid biosynthesis. Thus, a potential link between cadmium stress response and flavonoid synthesis was identified through protein interaction network analysis. The study thus provided significant information into MYB regulatory genes in ramie and may serve as a foundation for genetic enhancement and increased productivity.

Genome-Wide Characterization and Analysis of R2R3-MYB Genes Related to Fruit Ripening and Stress Response in Banana (Musa acuminata L. AAA Group, cv. 'Cavendish')

MYB is an important type of transcription factor in eukaryotes. It is widely involved in a variety of biological processes and plays a role in plant morphogenesis, growth and development, primary and secondary metabolite synthesis, and other life processes. In this study, bioinformatics methods were used to identify the R2R3-MYB transcription factor family members in the whole Musa acuminata (DH-Pahang) genome, one of the wild ancestors of banana. A total of 280 MaMYBs were obtained, and phylogenetic analysis indicated that these MaMYBs could be classified into 33 clades with MYBs from Arabidopsis thaliana. The amino acid sequences of the R2 and R3 Myb-DNA binding in all MaMYB protein sequences were quite conserved, especially Arg-12, Arg-13, Leu-23, and Leu-79. Distribution mapping results showed that 277 MaMYBs were localized on the 11 chromosomes in the Musa acuminata genome. The MaMYBs were distributed unevenly across the 11 chromosomes. More than 40.0% of the MaMYBs were located in collinear fragments, and segmental duplications likely played a key role in the expansion of the MaMYBs. Moreover, the expression profiles of MaMYBs in different fruit development and ripening stages and under various abiotic and biotic stresses were investigated using available RNA-sequencing data to obtain fruit development, ripening-specific, and stress-responsive candidate genes. Weighted gene co-expression network analysis (WGCNA) was used to analyze transcriptome data of banana from the above 11 samples. We found MaMYBs participating in important metabolic biosynthesis pathways in banana. Collectively, our results represent a comprehensive genome-wide study of the MaMYB gene family, which should be helpful in further detailed studies on MaMYBs functions related to fruit development, postharvest ripening, and the seedling response to stress in an important banana cultivar.

Genome-Wide Analysis of Specific PfR2R3-MYB Genes Related to Paulownia Witches' Broom

Paulownia witches' broom (PaWB), caused by phytoplasmas, is the most devastating infectious disease of Paulownia. R2R3-MYB transcription factors (TF) have been reported to be involved in the plant's response to infections caused by these pathogens, but a comprehensive study of the R2R3-MYB genes in Paulownia has not been reported. In this study, we identified 138 R2R3-MYB genes distributed on 20 chromosomes of Paulownia fortunei. These genes were classified into 27 subfamilies based on their gene structures and phylogenetic relationships, which indicated that they have various evolutionary relationships and have undergone rich segmental replication events. We determined the expression patterns of the 138 R2R3-MYB genes of P. fortunei by analyzing the RNA sequencing data and found that PfR2R3-MYB15 was significantly up-regulated in P. fortunei in response to phytoplasma infections. PfR2R3-MYB15 was cloned and overexpressed in Populus trichocarpa. The results show that its overexpression induced branching symptoms. Subsequently, the subcellular localization results showed that PfR2R3-MYB15 was located in the nucleus. Yeast two-hybrid and bimolecular fluorescence complementation experiments showed that PfR2R3-MYB15 interacted with PfTAB2. The analysis of the PfR2R3-MYB15 gene showed that it not only played an important role in plant branching, but also might participate in the biosynthesis of photosystem elements. Our results will provide a foundation for future studies of the R2R3-MYB TF family in Paulownia and other plants.

A Novel R2R3-MYB Transcription Factor SbMYB12 Positively Regulates Baicalin Biosynthesis in Scutellaria baicalensis Georgi

Scutellaria baicalensis Georgi is an annual herb from the Scutellaria genus that has been extensively used as a traditional medicine for over 2000 years in China. Baicalin and other flavonoids have been identified as the principal bioactive ingredients. The biosynthetic pathway of baicalin in S. baicalensis has been elucidated; however, the specific functions of R2R3-MYB TF, which regulates baicalin synthesis, has not been well characterized in S. baicalensis to date. Here, a S20 R2R3-MYB TF (SbMYB12), which encodes 263 amino acids with a length of 792 bp, was expressed in all tested tissues (mainly in leaves) and responded to exogenous hormone methyl jasmonate (MeJA) treatment. The overexpression of SbMYB12 significantly promoted the accumulation of flavonoids such as baicalin and wogonoside in S. baicalensis hairy roots. Furthermore, biochemical experiments revealed that SbMYB12 is a nuclear-localized transcription activator that binds to the SbCCL7-4, SbCHI-2, and SbF6H-1 promoters to activate their expression. These results illustrate that SbMYB12 positively regulates the generation of baicalin and wogonoside. In summary, this work revealed a novel S20 R2R3-MYB regulator and enhances our understanding of the transcriptional and regulatory mechanisms of baicalin biosynthesis, as well as sheds new light on metabolic engineering in S. baicalensis.

Investigation and Expression Analysis of R2R3-MYBs and Anthocyanin Biosynthesis-Related Genes during Seed Color Development of Common Bean (Phaseolus vulgaris)

Anthocyanins are responsible for the coloration of common bean seeds, and their accumulation is positively correlated with the expression level of anthocyanin biosynthetic genes. The MBW (MYB-bHLH-WD40) complex is thought to regulate the expression of these genes, and MYB proteins, which are a key factor in activating anthocyanin pathway genes, have been identified in several plants. This study demonstrated gene structures, chromosomal placements, gene duplications of R2R3-MYBs, miRNAs associated with R2R3-MYBs, and the interaction of these genes with other flavonoid regulatory genes. qRT-PCR was used to investigate the role of specific R2R3-MYBs and flavonoid genes in common bean seed color development. As a result of a comprehensive analysis with the help of in silico tools, we identified 160 R2R3-MYB genes in the common bean genome. We divided these genes into 16 classes on the basis of their intron-exon and motif structures. Except for three, the rest of the common bean R2R3-MYB members were distributed to all chromosomes with different densities, primarily located on chromosomes 3 and 8. We identified a total of 44 duplicated gene pairs dispersed across 11 chromosomes and evolved under purifying selection (Ka/Ks  <  1), 19 of which were derived from a whole-genome duplication. Our research uncovered 25 putative repressor PvMYB proteins that contain the EAR motif. Additionally, fifty different cis-regulatory elements regulated by light, stress, and hormone were identified. Within the genome of the common bean, we discovered a total of 36 microRNAs that target a total of 72 R2R3-MYB transcripts. The effect of 16 R2R3-MYB genes and 16 phenylpropanoid pathway genes, selected on the basis of their interaction in the protein-protein interaction map, playing role in the regulation of seed coat color development was evaluated using qRT-PCR in 5 different tissues at different developmental stages. The results revealed that these specific genes have different expression levels during different developmental periods, with higher levels in the pod filling and early pod stages than in the rest of the developmental periods. Furthermore, it was shown that PvTT8 (bHLH), PvTT2 (PvMYB42), PvMYB113, PvTTG1, and PvWD68 genes have effects on the regulation of seed coat color. The findings of this study, which is the first to use whole-genome analysis to identify and characterize the R2R3-MYB genes in common bean, may serve as a reference for future functional research in the legume.

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.

Genome-Wide Identification and Characterization of the SBP Gene Family in Passion Fruit (Passiflora edulis Sims)

The SQUAMOSA promoter binding proteins (SBPs) gene family plays important roles in plant growth and development. The SBP gene family has been identified and reported in many species, but it has not been well studied in passion fruit. In this study, a total of 14 SBP genes were identified in passion fruit and named from PeSBP1 to PeSBP14 based on their chromosomal distribution. The phylogenetic tree, gene structure, conserved motifs, collinearity analysis, and expression patterns of the identified SBP members were analyzed. We classified the PeSBP genes into eight groups (I to VIII) according to the phylogenetic tree, gene structure, and conserved motifs. Synteny analysis found that 5 homologous gene pairs existed in PeSBP genes and 11 orthologous gene pairs existed between passion fruit and Arabidopsis. Synonymous nucleotide substitution analysis showed that the PeSBP genes were under strong negative selection. The expression pattern of PeSBP genes in seed, root, leaf, and flower showed that nine of the PeSBP genes displayed high expression in the leaf and the flower. The expression patterns of PeSBP3/6/8/9/10 were further detected by qRT-PCR. In addition, differences in the expression levels occurred for each gene in the different flower organs and at the different developmental stages. There were large differences among SBPs based on transcriptional levels under cold, heat, salt, and osmotic stress conditions. Altogether, this study provides an overview of SBP genes in passion fruit and lays the foundation for further functional analysis.

Deciphering the roles of tobacco MYB transcription factors in environmental stress tolerance

The MYB members play important roles in development, metabolism, and stress tolerance in plants. In the current study, a total of 246 tobacco R2R3-MYB transcription factors were identified and systemically analyzed from the latest genome annotation. The newly identified tobacco members were divided into 33 subgroups together with the Arabidopsis members. Furthermore, 44 NtMYB gene pairs were identified to arise from duplication events, which might lead to the expansion of tobacco MYB genes. The expression patterns were revealed by transcriptomic analysis. Notably, the results from phylogenetic analysis, synthetic analysis, and expression analysis were integrated to predict the potential functions of these members. Particularly, NtMYB102 was found to act as the homolog of AtMYB70 and significantly induced by drought and salt treatments. The further assays revealed that NtMYB102 had transcriptional activities, and the overexpression of the encoding gene enhanced the drought and salt stress tolerance in transgenic tobacco. The results of this study may be relevant for future functional analyses of the MYB genes in tobacco.

Genome-wide analysis of R2R3-MYB transcription factors in Japanese morning glory

The R2R3-MYB transcription factor is one of the largest transcription factor families in plants. R2R3-MYBs play a variety of functions in plants, such as cell fate determination, organ and tissue differentiations, primary and secondary metabolisms, stress and defense responses and other physiological processes. The Japanese morning glory (Ipomoea nil) has been widely used as a model plant for flowering and morphological studies. In the present study, 127 R2R3-MYB genes were identified in the Japanese morning glory genome. Information, including gene structure, protein motif, chromosomal location and gene expression, were assigned to the InR2R3-MYBs. Phylogenetic tree analysis revealed that the 127 InR2R3-MYBs were classified into 29 subfamilies (C1-C29). Herein, physiological functions of the InR2R3-MYBs are discussed based on the functions of their Arabidopsis orthologues. InR2R3-MYBs in C9, C15, C16 or C28 may regulate cell division, flavonol biosynthesis, anthocyanin biosynthesis or response to abiotic stress, respectively. C16 harbors the known anthocyanin biosynthesis regulator, InMYB1 (INIL00g10723), and putative anthocyanin biosynthesis regulators, InMYB2 (INIL05g09650) and InMYB3 (INIL05g09651). In addition, INIL05g09649, INIL11g40874 and INIL11g40875 in C16 were suggested as novel anthocyanin biosynthesis regulators. We organized the R2R3-MYB transcription factors in the morning glory genome and assigned information to gene and protein structures and presuming their functions. Our study is expected to facilitate future research on R2R3-MYB transcription factors in Japanese morning glory.

Comparative analysis of R2R3-MYB transcription factors in the flower of Iris laevigata identifies a novel gene regulating tobacco cold tolerance

Breeding for flower cold resistance is a priority for flower breeding research in northern China. The identification of cold resistance genes will not only provide genetic resources for cold resistance breeding, but also form a basis for the study of plant cold resistance mechanisms. Based on the flower transcriptome of Iris laevigata, 20 R2R3-MYBs were identified and comprehensive analysis, including conservative domain, phylogenetic analyses and functional distribution, were performed for R2R3-MYBs. Expression patterns of the abiotic stress genes under cold stress were detected, the upregulated gene was genetically transformed into tobacco, and the related physiological indicators of the transgenic tobacco were measured. A novel cold resistance gene, IlMYB306, was obtained. qRT-PCR indicated that IlMYB306 was dramatically induced by cold stress and was significantly upregulated in roots. The free proline content, MDA, SOD and POD activity of the transgenic tobacco improved after cold stress, and the chlorophyll content decreased slowly. In addition, overexpression of IlMYB306 improved cold resistance of the seeds. SEM results showed leaves of transgenic tobacco had obvious folds, more grooves and bulges on the lower leaf surface. Overall, we report a novel cold resistance R2R3-MYB gene, IlMYB306, in the flower of I. laevigata, which could improve tobacco cold stress tolerance by thickening the waxy layer, increasing antioxidant activity and the content of proline.

An R2R3-MYB transcription factor VyMYB24, isolated from wild grape Vitis yanshanesis J. X. Chen., regulates the plant development and confers the tolerance to drought

In grapevines, the MYB transcription factors play an important regulatory role in the phenylpropanoid pathway including proanthocyanidin, anthocyanin, and flavonoid biosynthesis. However, the role of MYB in abiotic stresses is not clear. In this study, an R2R3-MYB transcription factor, VyMYB24, was isolated from a high drought-tolerant Chinese wild Vitis species V. yanshanesis. Our findings demonstrated that it was involved in plant development and drought tolerance. VyMYB24 is a nuclear protein and is significantly induced by drought stress. When over-expressed in tobacco, VyMYB24 caused plant dwarfing including plant height, leaf area, flower size, and seed weight. The GA1+3 content in transgenic plants was reduced significantly, and spraying exogenous gibberellin could recover the dwarf phenotype of VyMYB24 transgenic plants, suggesting that VyMYB24 might inhibit plant development by the regulation of gibberellin (GA) metabolism. Under drought stress, the VyMYB24 transgenic plants improved their tolerance to drought with a lower wilting rate, lower relative electrical conductivity, and stronger roots. Compared to wild-type tobacco plants, VyMYB24 transgenic plants accumulated less reactive oxygen, accompanied by increased antioxidant enzyme activity and upregulated gene expression levels of superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) genes. In addition, transgenic plants accumulated more proline, and their related synthetic genes NtP5CR and NtP5CS genes were significantly upregulated when exposed to drought. Besides, abiotic stress-responsive genes, NtDREB, NtERD10C, NtERD10D, and NtLEA5, were upregulated significantly in VyMYB24 transgenic plants. These results indicate that VyMYB24 plays a positive regulatory role in response to drought stress and also regulates plant development, which provides new evidence to further explore the molecular mechanism of drought stress of the MYB gene family.

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.