The R2R3-MYB gene family in Arabidopsis thaliana

Genome-wide characterization of R2R3-MYB gene family in Santalum album and their expression analysis under cold stress

Sandalwood (Santalum album) is a high-value multifunctional tree species that is rich in aromatic substances and is used in medicine and global cosmetics. Due to the scarcity of land resources in tropical and subtropical regions, land in temperate regions is a potential resource for the development of S. album plantations in order to meet the needs of S. album production and medicine. The R2R3-MYB transcription factor family is one of the largest in plants and plays an important role in the response to various abiotic stresses. However, the R2R3-MYB gene family of S. album has not been studied. In this study, 144 R2R3-MYB genes were successfully identified in the assembly genome sequence, and their characteristics and expression patterns were investigated under various durations of low temperature stress. According to the findings, 31 of the 114 R2R3-MYB genes showed significant differences in expression after cold treatment. Combining transcriptome and weighted gene co-expression network analysis (WGCNA) revealed three key candidate genes (SaMYB098, SaMYB015, and SaMYB068) to be significantly involved in the regulation of cold resistance in S. album. The structural characteristics, evolution, and expression pattern of the R2R3-MYB gene in S. album were systematically examined at the whole genome level for the first time in this study. It will provide important information for future research into the function of the R2R3-MYB genes and the mechanism of cold stress response in S. album.

Molecular mechanism of naringenin regulation on flavonoid biosynthesis to improve the salt tolerance in pigeon pea (Cajanus cajan (Linn.) Millsp.)

Flavonoids are important secondary metabolites in the plant growth and development process. As a medicinal plant, pigeon pea is rich in secondary metabolites. As a flavonoid, there are few studies on the regulation mechanism of naringenin in plant stress resistance. In our study, we found that naringenin can increase the pigeon pea's ability to tolerate salt and influence the changes that occur in flavonoids including naringenin, genistein and biochanin A. We analyzed the transcriptome data after 1 mM naringenin treatment, and identified a total of 13083 differentially expressed genes. By analyzing the metabolic pathways of these differentially expressed genes, we found that these differentially expressed genes were enriched in the metabolic pathways of phenylpropanoid biosynthesis, starch and sucrose metabolism and so on. We focused on the analysis of flavonoid biosynthesis related pathways. Among them, the expression levels of enzyme genes CcIFS, CcCHI and CcCHS in the flavonoid biosynthesis pathway had considerably higher expression levels. By counting the number of transcription factors and the binding sites on the promoter of the enzyme gene, we screened the transcription factors CcMYB62 and CcbHLH35 related to flavonoid metabolism. Among them, CcMYB62 has a higher expression level than the others. The hairy root transgene showed that CcMYB62 could induce the upregulation of CcCHI, and promote the accumulation of naringenin, genistein and biochanin A. Our study revealed the molecular mechanism of naringenin regulating flavonoid biosynthesis under salt stress in pigeon pea, and provided an idea for the role of flavonoids in plant resistance to abiotic stresses.

Mdm-miR858 targets MdMYB9 and MdMYBPA1 to participate anthocyanin biosynthesis in red-fleshed apple

Anthocyanins are important secondary metabolites in plants. They are important for human health because of their antioxidant activities and because their dietary intake reduces the incidence of cardiovascular and cerebrovascular diseases and tumors. The biosynthesis of anthocyanins and its regulation in fruits and vegetables is a global research hotspot. Compared with cultivated apples, the red-fleshed apple is a relatively new and popular commodity in the market. Previous studies on red-fleshed apples have focused on the basis for the high anthocyanin content and the transcriptional regulation of anthocyanin synthesis. In the present study, we focused on the mechanism of microRNA-mediated post-transcriptional regulation of anthocyanin synthesis in red-fleshed apples. We identified a microRNA (miRNA), designated mdm-miR858, that is specifically expressed in the flesh of apple fruit. The expression level of miR858 was significantly lower in red-fleshed apples than in white-fleshed apples. The overexpression of mdm-miR858 significantly inhibited anthocyanin accumulation, whereas the silencing of mdm-miR858 promoted anthocyanin synthesis in STTM858 transgenic apple calli. Further analyses showed that mdm-miR858 targets the transcription factor genes MdMYB9 and MdMYBPA1 to participate anthocyanin accumulation in apple. Our results also show that MdHY5, a transcription factor in the light signaling pathway, can bind to the promoter of mdm-miR858 to inhibit its transcription, thereby regulating anthocyanin synthesis. Based on our results, we describe a novel HY5-miR858-MYB loop involved in the modulation of anthocyanin biosynthesis. These findings provide new information about how plant miRNAs regulate anthocyanin anabolism and provide a basis for breeding new anthocyanin-rich, red-fleshed apple varieties.

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.

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.

Fine mapping and characterisation of a PV-PUR mediating anthocyanin synthesis in snap bean (Phaseolus vulgaris L.)

Anthocyanin makes snap bean (Phaseolus vulgaris L.) pods purple, which helps seed dispersal and protects against environmental stress. In this study, we characterised the snap bean purple mutant pv-pur, which has purple cotyledon, hypocotyl, stem, leaf vein, flower and pod tissues. Total anthocyanin, delphinidin and malvidin levels in mutant pods were significantly higher than in wild-type plants. We constructed two populations for fine mapping of the PV-PUR purple mutation gene, located in the 243.9-kb region of chromosome 06. We identified Phvul.006g018800.3, encoding F3'5'H, as a candidate gene for PV-PUR. Six single-base mutations occurred in the coding region of this gene, altering protein structure. PV-PUR and pv-pur genes were transferred into Arabidopsis, respectively. Compared with the wild-type, the leaf base and internode of T-PV-PUR plant were purple, and the phenotype of T-pv-pur plant remained unchanged, which verified the function of the mutant gene. The results demonstrated that PV-PUR is a crucial gene for anthocyanin biosynthesis in snap bean, resulting in purple colouration. The findings lay a foundation for future breeding and improvement of snap bean.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01362-8.

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.

Genome-Wide Analysis of MYB Transcription Factors and Screening of MYBs Involved in the Red Color Formation in Rhododendron delavayi

Flower color is one of the crucial traits of ornamental plants. Rhododendron delavayi Franch. is a famous ornamental plant species distributed in the mountain areas of Southwest China. This plant has red inflorescence and young branchlets. However, the molecular basis of the color formation of R. delavayi is unclear. In this study, 184 MYB genes were identified based on the released genome of R. delavayi. These genes included 78 1R-MYB, 101 R2R3-MYB, 4 3R-MYB, and 1 4R-MYB. The MYBs were divided into 35 subgroups using phylogenetic analysis of the MYBs of Arabidopsis thaliana. The members of the same subgroup in R. delavayi had similar conserved domains and motifs, gene structures, and promoter cis-acting elements, which indicate their relatively conserved function. In addition, transcriptome based on unique molecular identifier strategy and color difference of the spotted petals, unspotted petals, spotted throat, unspotted throat, and branchlet cortex were detected. Results showed significant differences in the expression levels of R2R3-MYB genes. Weighted co-expression network analysis between transcriptome and chromatic aberration values of five types of red samples showed that the MYBs were the most important TFs involved in the color formation, of which seven were R2R3-MYB, and three were 1R-MYB. Two R2R3-MYB (DUH019226.1 and DUH019400.1) had the highest connectivity in the whole regulation network, and they were identified as hub genes for red color formation. These two MYB hub genes provide references for the study of transcriptional regulation of the red color formation of R. delavayi.

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.

AtBBX29 integrates photomorphogenesis and defense responses in Arabidopsis

Light is an environmental signal that modulates plant defenses against attackers. Recent research has focused on the effects of light on defense hormone signaling; however, the connections between light signaling pathways and the biosynthesis of specialized metabolites involved in plant defense have been relatively unexplored. Here, we show that Arabidopsis BBX29, a protein that belongs to the B-Box transcription factor (TF) family, integrates photomorphogenic signaling with defense responses by promoting flavonoid, sinapate and glucosinolate accumulation in Arabidopsis leaves. AtBBX29 transcript levels were up regulated by light, through photoreceptor signaling pathways. Genetic evidence indicated that AtBBX29 up-regulates MYB12 gene expression, a TF known to induce genes related to flavonoid biosynthesis in a light-dependent manner, and MYB34 and MYB51, which encode TFs involved in the regulation of glucosinolate biosynthesis. Thus, bbx29 knockout mutants displayed low expression levels of key genes of the flavonoid biosynthetic pathway, and the opposite was true in BBX29 overexpression lines. In agreement with the transcriptomic data, bbx29 mutant plants accumulated lower levels of kaempferol glucosides, sinapoyl malate, indol-3-ylmethyl glucosinolate (I3M), 4-methylsulfinylbutyl glucosinolate (4MSOB) and 3-methylthiopropyl glucosinolate (3MSP) in rosette leaves compared to the wild-type, and showed increased susceptibility to the necrotrophic fungus Botrytis cinerea and to the herbivore Spodoptera frugiperda. In contrast, BBX29 overexpressing plants displayed increased resistance to both attackers. In addition, we found that AtBBX29 plays an important role in mediating the effects of ultraviolet-B (UV-B) radiation on plant defense against B. cinerea. Taken together, these results suggest that AtBBX29 orchestrates the accumulation of specific light-induced metabolites and regulates Arabidopsis resistance against pathogens and herbivores.

Integrative analysis of metabolome and transcriptome reveals a dynamic regulatory network of potato tuber pigmentation

Potatoes consist of flavonoids that provide health benefits for human consumers. To learn more about how potato tuber flavonoid accumulation and flesh pigmentation are controlled, we analyzed the transcriptomic and metabolomic profile of potato tubers from three colored potato clones at three developmental phases using an integrated approach. From the 72 flavonoids identified in pigmented flesh, differential abundance was noted for anthocyanins, flavonols, and flavones. Weighted gene co-expression network analysis further allowed modules and candidate genes that positively or negatively regulate flavonoid biosynthesis to be identified. Furthermore, an R2R3-MYB repressor StMYB3 and an R3-MYB repressor StMYBATV involved in the modulation of anthocyanin biosynthesis during tuber development were identified. Both StMYB3 and StMYBATV could interact with the cofactor StbHLH1 and repress anthocyanin biosynthesis. Our results indicate a feedback regulatory mechanism of a coordinated MYB activator-repressor network on fine-tuning of potato tuber pigmentation during tuber development.

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.

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.

Reduced glucosinolate content in oilseed rape (Brassica napus L.) by random mutagenesis of BnMYB28 and BnCYP79F1 genes

The presence of anti-nutritive compounds like glucosinolates (GSLs) in the rapeseed meal severely restricts its utilization as animal feed. Therefore, reducing the GSL content to < 18 µmol/g dry weight in the seeds is a major breeding target. While candidate genes involved in the biosynthesis of GSLs have been described in rapeseed, comprehensive functional analyses are missing. By knocking out the aliphatic GSL biosynthesis genes BnMYB28 and BnCYP79F1 encoding an R2R3 MYB transcription factor and a cytochrome P450 enzyme, respectively, we aimed to reduce the seed GSL content in rapeseed. After expression analyses on single paralogs, we used an ethyl methanesulfonate (EMS) treated population of the inbred winter rapeseed 'Express617' to detect functional mutations in the two gene families. Our results provide the first functional analysis by knock-out for the two GSL biosynthesis genes in winter rapeseed. We demonstrate that independent knock-out mutants of the two genes possessed significantly reduced seed aliphatic GSLs, primarily progoitrin. Compared to the wildtype Express617 control plants (36.3 µmol/g DW), progoitrin levels were decreased by 55.3% and 32.4% in functional mutants of BnMYB28 (16.20 µmol/g DW) and BnCYP79F1 (24.5 µmol/g DW), respectively. Our study provides a strong basis for breeding rapeseed with improved meal quality in the future.

Jasmonic acid negatively regulates branch growth in pear

The quality of seedlings is an important factor for development of the pear industry. A strong seedling with few branches and suitable internodes is ideal material as a rootstock for grafting and breeding. Several branching mutants of pear rootstocks were identified previously. In the present study, 'QAU-D03' (Pyrus communis L.) and it's mutants were used to explore the mechanism that affects branch formation by conducting phenotypic trait assessment, hormone content analysis, and transcriptome analysis. The mutant plant (MP) showed fewer branches, shorter 1-year-old shoots, and longer petiole length, compared to original plants (OP), i.e., wild type. Endogenous hormone analysis revealed that auxin, cytokinin, and jasmonic acid contents in the stem tips of MP were significantly higher than those of the original plants. In particular, the jasmonic acid content of the MP was 1.8 times higher than that of the original plants. Transcriptome analysis revealed that PcCOI1, which is a transcriptional regulatory gene downstream of the jasmonic acid signaling pathway, was expressed more highly in the MP than in the original plants, whereas the expression levels of PcJAZ and PcMYC were reduced in the MP compared with that of the original plants. In response to treatment with exogenous methyl jasmonate, the original plants phenotype was consistent with that of the MP in developing less branches. These results indicate that jasmonic acid negatively regulates branch growth of pear trees and that jasmonic acid downstream regulatory genes play a crucial role in regulating branching.

Genome-Wide Identification Analysis of the R2R3-MYB Transcription Factor Family in Cymbidium sinense for Insights into Drought Stress Responses

Cymbidium sinense represents a distinctive Orchidaceae plant that is more tolerant than other terrestrial orchids. Studies have shown that many members of the MYB transcription factor (TF) family, especially the R2R3-MYB subfamily, are responsive to drought stress. This study identified 103 CsMYBs; phylogenetic analysis classified these genes into 22 subgroups with Arabidopsis thaliana. Structural analysis showed that most CsMYB genes contained the same motifs, three exons and two introns, and showed a helix-turn-helix 3D structure in each R repeat. However, the members of subgroup 22 contained only one exon and no intron. Collinear analysis revealed that C. sinense had more orthologous R2R3-MYB genes with wheat than A. thaliana and rice. Ka/Ks ratios indicated that most CsMYB genes were under purifying negative selection pressure. Cis-acting elements analysis revealed that drought-related elements were mainly focused on subgroups 4, 8, 18, 20, 21, and 22, and Mol015419 (S20) contained the most. The transcriptome analysis results showed that expression patterns of most CsMYB genes were upregulated in leaves in response to slight drought stress and downregulated in roots. Among them, members in S8 and S20 significantly responded to drought stress in C. sinense. In addition, S14 and S17 also participated in these responses, and nine genes were selected for the real-time reverse transcription quantitative PCR (RT-qPCR) experiment. The results were roughly consistent with the transcriptome. Our results, thus, provide an important contribution to understanding the role of CsMYBs in stress-related metabolic processes.

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.

Comparative physiological and transcriptomic analyses reveal the mechanisms of CO2 enrichment in promoting the growth and quality in Lactuca sativa

The increase in the concentration of CO2 in the atmosphere has attracted widespread attention. To explore the effect of elevated CO2 on lettuce growth and better understand the mechanism of elevated CO2 in lettuce cultivation, 3 kinds of lettuce with 4 real leaves were selected and planted in a solar greenhouse. One week later, CO2 was applied from 8:00 a.m. to 10:00 a.m. on sunny days for 30 days. The results showed that the growth potential of lettuce was enhanced under CO2 enrichment. The content of vitamin C and chlorophyll in the three lettuce varieties increased, and the content of nitrate nitrogen decreased. The light saturation point and net photosynthetic rate of leaves increased, and the light compensation point decreased. Transcriptome analysis showed that there were 217 differentially expressed genes (DEGs) shared by the three varieties, among which 166 were upregulated, 44 were downregulated, and 7 DEGs were inconsistent in the three materials. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that these DEGs involved mainly the ethylene signaling pathway, jasmonic acid signaling pathway, porphyrin and chlorophyll metabolism pathway, starch and sucrose metabolism pathway, etc. Forty-one DEGs in response to CO2 enrichment were screened out by Gene Ontology (GO) analysis, and the biological processes involved were consistent with KEGG analysis. which suggested that the growth and nutritional quality of lettuce could be improved by increasing the enzyme activity and gene expression levels of photosynthesis, hormone signaling and carbohydrate metabolism. The results laid a theoretical foundation for lettuce cultivation in solar greenhouses and the application of CO2 fertilization technology.

Chromosome-Level Assembly of Flowering Cherry (Prunus campanulata) Provides Insight into Anthocyanin Accumulation

The flowering cherries (genus Prunus, subgenus Cerasus) are popular ornamental trees in China, Japan, Korea, and elsewhere. Prunus campanulata Maxim. is an important species of flowering cherry native to Southern China, which is also distributed in Taiwan, the Ryukyu Islands of Japan, and Vietnam. It produces bell-shaped flowers with colors ranging from bright pink to crimson during the Chinese Spring Festival from January to March each year. We selected the P. campanulata cultivar "Lianmeiren", with only 0.54% of heterozygosity, as the focus of this study, and generated a high-quality chromosome-scale genome assembly of P. campanulata by combining Pacific Biosciences (PacBio) single-molecule sequencing, 10× Genomics sequencing, and high-throughput chromosome conformation capture (Hi-C) technology. We first assembled a 300.48 Mb genome assembly with a contig N50 length of 2.02 Mb. In total, 28,319 protein-coding genes were predicted from the genome, 95.8% of which were functionally annotated. Phylogenetic analyses indicated that P. campanulata diverged from a common ancestor of cherry approximately 15.1 million years ago. Comparative genomic analyses showed that the expanded gene families were significantly involved in ribosome biogenesis, diterpenoid biosynthesis, flavonoid biosynthesis, and circadian rhythm. Furthermore, we identified 171 MYB genes from the P. campanulata genome. Based on the RNA-seq of five organs at three flowering stages, expression analyses revealed that the majority of the MYB genes exhibited tissue-specific expression patterns, and some genes were identified as being associated with anthocyanin accumulation. This reference sequence is an important resource for further studies of floral morphology and phenology, and comparative genomics of the subgenera Cerasus and Prunus.

OsPPR11 encoding P-type PPR protein that affects group II intron splicing and chloroplast development

OsPPR11 belongs to the P-type PPR protein family and can interact with OsCAF2 to regulate Group II intron splicing and affect chloroplast development in rice. Pentatricopeptide repeat (PPR) proteins participate in chloroplasts or mitochondria group II introns splicing in plants. The PPR protein family contains 491 members in rice, but most of their functions are unknown. In this study, we identified a nuclear gene encoding the P-type PPR protein OsPPR11 in chloroplasts. The qRT-PCR analysis demonstrated that OsPPR11 was expressed in all plant tissues, but leaves had the highest expression. The osppr11 mutants had yellowing leaves and a lethal phenotype that inhibited chloroplast development and photosynthesis-related gene expression and reduced photosynthesis-related protein accumulation in seedlings. Moreover, photosynthetic complex accumulation decreased significantly in osppr11 mutants. The OsPPR11 is required for ndhA, and ycf3-1 introns splicing and interact with CRM family protein OsCAF2, suggesting that these two proteins may form splicing complexes to regulate group II introns splicing. Further analysis revealed that OsCAF2 interacts with OsPPR11 through the N-terminus. These results indicate that OsPPR11 is essential for chloroplast development and function by affecting group II intron splicing in rice.

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.

Foxtail millet MYB-like transcription factor SiMYB16 confers salt tolerance in transgenic rice by regulating phenylpropane pathway

R2R3-MYB transcription factors play an important role in the synthesis of phenylpropanoid-derived compounds, which in turn provide salt tolerance in plant. In this study, we found that the expression of foxtail millet R2R3-MYB factor SiMYB16 can be induced by salt and drought. SiMYB16 is localized in the nucleus and acts as a transcriptional activator. Phylogenetic analysis indicates that SiMYB16 belongs to the R2R3-MYB transcription factor family subgroup 24. Transgenic rice expressing SiMYB16 (OX16) had a higher survival rate, lower malondialdehyde content, and heavier fresh weight compared with type (WT) under salt stress conditions. The transgenic plants also had a higher germination rate in salt treatment conditions and higher yield in the field compared with wild-type plants. Transcriptome analysis revealed that the up-regulated differential expression genes in the transgenic rice were mainly involved in phenylpropanoid biosynthesis, fatty acid elongation, phenylalanine metabolism, and flavonoid biosynthesis pathways. Quantitative real-time PCR analysis also showed that the genes encoding the major enzymes in the lignin and suberin biosynthesis pathways had higher expression level in SiMYB16 transgenic plants. Correspondingly, the content of flavonoid and lignin, and the activity of fatty acid synthase increased in SiMYB16 transgenic rice compared with wild-type plants under salt stress treatment. These results indicate that SiMYB16 gene can enhance plant salt tolerance by regulating the biosynthesis of lignin and suberin.

A MYB Transcription Factor Atlas Provides Insights into the Evolution of Environmental Adaptations in Plants

The MYB transcription factor superfamily includes key regulators of plant development and responses to environmental changes. The diversity of lifestyles and morphological characteristics exhibited by plants are potentially associated with the genomic dynamics of the MYB superfamily. With the release of the plant genomes, a comprehensive phylogenomic analysis of the MYB superfamily across Viridiplantae is allowed. The present study performed phylogenetic, phylogenomic, syntenic, horizontal gene transfer, and neo/sub-functionalization analysis of the MYB superfamily to explore the evolutionary contributions of MYB members to species diversification, trait formation, and environmental adaptation in 437 different plant species. We identified major changes in copy number variation and genomic context within subclades across lineages. Multiple MYB subclades showed highly conserved copy number patterns and synteny across flowering plants, whereas others were more dynamic and showed lineage-specific patterns. As examples of lineage-specific morphological divergence, we hypothesize that the gain of a MYB orthogroup associated with flower development and environmental responses and an orthogroup associated with auxin and wax biosynthesis in angiosperms were correlated with the emergence of flowering plants, unbiased neo-/sub-functionalization of gene duplicates contributed to environmental adaptation, and species-specific neo-/sub-functionalization contributed to phenotype divergence between species. Transposable element insertion in promoter regions may have facilitated the sub-/neo-functionalization of MYB genes and likely played a tissue-specific role contributing to sub-/neo-functionalization in plant root tissues. This study provides new insights into the evolutionary divergence of the MYB superfamily across major flowering and non-flowering lineages and emphasizes the need for lineage-/tissue-specific characterization to further understand trait variability and environmental adaptation.

Genome-Wide Identification and Transcriptional Analysis of the MYB Gene Family in Pearl Millet (Pennisetum glaucum)

The MYB gene family widely exists in the plant kingdom and participates in the regulation of plant development and stress response. Pearl millet (Pennisetum glaucum (L.) R. Br.), as one of the most important cereals, is not only considered a good source of protein and nutrients but also has excellent tolerances to various abiotic stresses (e.g., salinity, water deficit, etc.). Although the genome sequence of pearl millet was recently published, bioinformatics and expression pattern analysis of the MYB gene family are limited. Here, we identified 208 PgMYB genes in the pearl millet genome and employed 193 high-confidence candidates for downstream analysis. Phylogenetic and structural analysis classified these PgMYBs into four subgroups. Eighteen pairs of segmental duplications of the PgMYB gene were found using synteny analysis. Collinear analysis revealed pearl millet had the closest evolutionary relationship with foxtail millet. Nucleotide substitution analysis (Ka/Ks) revealed PgMYB genes were under purifying positive selection pressure. Reverse transcription-quantitative PCR analysis of eleven R2R3-type PgMYB genes revealed they were preferentially expressed in shoots and seeds and actively responded to various environment stimuli. Current results provide insightful information regarding the molecular features of the MYB family in pearl millet to support further functional characterizations.

RNA-Seq reveals the key pathways and genes involved in the light-regulated flavonoids biosynthesis in mango (Mangifera indica L.) peel

Introduction

Flavonoids are important water soluble secondary metabolites in plants, and light is one of the most essential environmental factors regulating flavonoids biosynthesis. In the previous study, we found bagging treatment significantly inhibited the accumulation of flavonols and anthocyanins but promoted the proanthocyanidins accumulation in the fruit peel of mango (Mangifera indica L.) cultivar 'Sensation', while the relevant molecular mechanism is still unknown.

Methods

In this study, RNA-seq was conducted to identify the key pathways and genes involved in the light-regulated flavonoids biosynthesis in mango peel.

Results

By weighted gene co-expression network analysis (WGCNA), 16 flavonoids biosynthetic genes were crucial for different flavonoids compositions biosynthesis under bagging treatment in mango. The higher expression level of LAR (mango026327) in bagged samples might be the reason why light inhibits proanthocyanidins accumulation in mango peel. The reported MYB positively regulating anthocyanins biosynthesis in mango, MiMYB1, has also been identified by WGCNA in this study. Apart from MYB and bHLH, ERF, WRKY and bZIP were the three most important transcription factors (TFs) involved in the light-regulated flavonoids biosynthesis in mango, with both activators and repressors. Surprisingly, two HY5 transcripts, which are usually induced by light, showed higher expression level in bagged samples.

Discussion

Our results provide new insights of the regulatory effect of light on the flavonoids biosynthesis in mango fruit peel.

RNAseq-Based Working Model for Transcriptional Regulation of Crosstalk between Simultaneous Abiotic UV-B and Biotic Stresses in Plants

Plants adjust their secondary metabolism by altering the expression of corresponding genes to cope with both abiotic and biotic stresses. In the case of UV-B radiation, plants produce protective flavonoids; however, this reaction is impeded during pattern-triggered immunity (PTI) induced by pathogens. Pathogen attack can be mimicked by the application of microbial associated molecular patterns (e.g., flg22) to study crosstalk between PTI and UV-B-induced signaling pathways. Switching from Arabidopsis cell cultures to in planta studies, we analyzed whole transcriptome changes to gain a deeper insight into crosstalk regulation. We performed a comparative transcriptomic analysis by RNAseq with four distinct mRNA libraries and identified 10778, 13620, and 11294 genes, which were differentially expressed after flg22, UV-B, and stress co-treatment, respectively. Focusing on genes being either co-regulated with the UV-B inducible marker gene chalcone synthase CHS or the flg22 inducible marker gene FRK1 identified a large set of transcription factors from diverse families, such as MYB, WRKY, or NAC. These data provide a global view of transcriptomic reprogramming during this crosstalk and constitute a valuable dataset for further deciphering the underlying regulatory mechanism(s), which appear to be much more complex than previously anticipated. The possible involvement of MBW complexes in this context is discussed.

Betalain biosynthesis in red pulp pitaya is regulated via HuMYB132: a R-R type MYB transcription factor

BACKGROUND

Multiple MYB transcription factors (TFs) are involved in the regulation of plant coloring. Betalain is a kind of natural plant pigment and its biosynthesis is regulated by a number of enzymes. Despite this, little is known about the molecular properties and roles of MYB TFs in pitaya betalain biosynthesis.

RESULTS

In the present study, we identified a 1R-MYB gene, HuMYB132, which is preferentially expressed in red-pulp pitaya at the mature stage. It was clustered with Arabidopsis R-R-type genes and had two DNA-binding domains and a histidine-rich region. The expression assays in N. benthamiana and yeast indicated that HuMYB132 is a nucleus-localized protein with transcriptional activation activity. Dual luciferase reporter assay and electrophoretic mobility shift assays (EMSA) demonstrated that HuMYB132 could promote the transcriptional activities of HuADH1, HuCYP76AD1-1, and HuDODA1 by binding to their promoters. Silencing HuMYB132 reduced betalain accumulation and the expression levels of betalain biosynthetic genes in pitaya pulps.

CONCLUSIONS

According to our findings, HuMYB132, a R-R type member of 1R-MYB TF subfamily, positively regulates pitaya betalain biosynthesis by regulating the expression of HuADH1, HuCYP76AD1-1, and HuDODA1. The present study provides a new theoretical reference for the management of pitaya betalain biosynthesis and also provides an essential basis for future regulation of betalain biosynthesis in Hylocereus.

A Genome-Wide View of the Transcriptome Dynamics of Fresh-Cut Potato Tubers

Fresh fruits and vegetable products are easily perishable during postharvest handling due to enzymatic browning reactions. This phenomenon has contributed to a significant loss of food. To reveal the physiological changes in fresh-cut potato tubers at the molecular level, a transcriptome analysis of potato tubers after cutting was carried out. A total of 10,872, 10,449, and 11,880 differentially expressed genes (DEGs) were identified at 4 h, 12 h and 24 h after cutting, respectively. More than 87.5% of these DEGs were classified into the categories of biological process (BP) and molecular function (MF) based on Gene Ontology (GO) analysis. There was a difference in the response to cutting at different stages after the cutting of potato tubers. The genes related to the phenol and fatty biosynthesis pathways, which are responsible for enzymatic browning and wound healing in potato tubers, were significantly enriched at 0-24 h after cutting. Most genes related to the enzymatic browning of potato tubers were up-regulated in response to cut-wounding. Plant hormone biosynthesis, signal molecular biosynthesis and transduction-related genes, such as gibberelin (GA), cytokinin (CK), ethylene (ET), auxin (IAA), jasmonic acid (JA), salicylic (SA), and Respiratory burst oxidase (Rboh) significantly changed at the early stage after cutting. In addition, the transcription factors involved in the wound response were the most abundant at the early stage after cutting. The transcription factor with the greatest response to injury was MYB, followed by AP2-EREBP, C3H and WRKY. This study revealed the physiological changes at the molecular level of fresh-cut potato tubers after cutting. This information is needed for developing a better approach to enhancing the postharvest shelf life of fresh processed potato and the breeding of potato plants that are resistant to enzymatic browning.

Comparative expression analysis in three Brassicaceae species revealed compensatory changes of the underlying gene regulatory network

Trichomes are regularly distributed on the leaves of Arabidopsis thaliana. The gene regulatory network underlying trichome patterning involves more than 15 genes. However, it is possible to explain patterning with only five components. This raises the questions about the function of the additional components and the identification of the core network. In this study, we compare the relative expression of all patterning genes in A. thaliana, A. alpina and C. hirsuta by qPCR analysis and use mathematical modelling to determine the relative importance of patterning genes. As the involved proteins exhibit evolutionary conserved differential complex formation, we reasoned that the genes belonging to the core network should exhibit similar expression ratios in different species. However, we find several striking differences of the relative expression levels. Our analysis of how the network can cope with such differences revealed relevant parameters that we use to predict the relevant molecular adaptations in the three species.

Characterization of Highbush Blueberry (Vaccinium corymbosum L.) Anthocyanin Biosynthesis Related MYBs and Functional Analysis of VcMYB Gene

As one of the most important transcription factors regulating plant anthocyanin biosynthesis, MYB has attracted great attentions. In this study, we identified fifteen candidate anthocyanin biosynthesis related MYB (ABRM) proteins, including twelve R2R3-MYBs and three 1R-MYBs, from highbush blueberry. The subcellular localization prediction results showed that, with the exception of VcRVE8 (localized in chloroplast and nucleus), all of the blueberry ABRMs were nucleus-localized. The gene structure analysis revealed that the exon numbers of the blueberry ABRM genes varied greatly, ranging between one and eight. There are many light-responsive, phytohormone-responsive, abiotic stress-responsive and plant growth and development related cis-acting elements in the promoters of the blueberry ABRM genes. It is noteworthy that almost all of their promoters contain light-, ABA- and MeJA-responsive elements, which is consistent with the well-established results that anthocyanin accumulation and the expression of MYBs are influenced significantly by many factors, such as light, ABA and JA. The gene expression analysis revealed that VcMYB, VcMYB6, VcMYB23, VcMYBL2 and VcPH4 are expressed abundantly in blueberry fruits, and VcMYB is expressed the highest in the red, purple and blue fruits among all blueberry ABRMs. VcMYB shared high similarity with functionally proven ABRMs from many other plant species. The gene cloning results showed that VcMYB had three variable transcripts, but only the transient overexpression of VcMYB-1 promoted anthocyanin accumulation in the green fruits. Our study can provide a basis for future research on the anthocyanin biosynthesis related MYBs in blueberry.

MYB30 and MYB14 form a repressor-activator module with WRKY8 that controls stilbene biosynthesis in grapevine

When exposed to pathogen infection or ultraviolet (UV) radiation, grapevine (Vitis vinifera) plants rapidly accumulate the stilbenoid resveratrol (Res) with concomitant increase of stilbene synthase (STS), the key enzyme in stilbene biosynthesis. Although a few transcription factors have been shown to regulate STSs, the molecular mechanism governing the regulation of STSs is not well elucidated. Our previous work showed that a VvMYB14-VvWRKY8 regulatory loop fine-tunes stilbene biosynthesis in grapevine through protein-protein interaction; overexpression of VvWRKY8 down-regulates VvMYB14 and VvSTS15/21; and application of exogenous Res up-regulates WRKY8 expression. Here, we identified an R2R3-MYB repressor, VvMYB30, which competes with the activator VvMYB14 for binding to the common binding sites in the VvSTS15/21 promoter. Similar to VvMYB14, VvMYB30 physically interacts with VvWRKY8 through their N-termini, forming a complex that does not bind DNA. Exposure to UV-B/C stress induces VvMYB14, VvWRKY8, and VvSTS15/21, but represses VvMYB30 in grapevine leaves. In addition, MYB30 expression is up-regulated by VvWRKY8-overexpression or exogenous Res. These findings suggest that the VvMYB14-VvWRKY8-VvMYB30 regulatory circuit allows grapevine to respond to UV stress by producing Res and prevents over-accumulation of Res to balance metabolic costs. Our work highlights the stress-mediated induction and feedback inhibition of stilbene biosynthesis through a complex regulatory network involving multiple positive and negative transcriptional regulators.

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.

R2R3 MYB transcription factor SbMYBHv33 negatively regulates sorghum biomass accumulation and salt tolerance

SbMYBHv33 negatively regulated biomass accumulation and salt tolerance in sorghum and Arabidopsis by regulating reactive oxygen species accumulation and ion levels. Salt stress is one of the main types of environmental stress leading to a reduction in crop yield worldwide. Plants have also evolved a variety of corresponding regulatory pathways to resist environmental stress damage. This study aimed to identify a SbMYBHv33 transcription factor that downregulates in salt, drought, and abscisic acid (ABA) in the salt-tolerant inbred line sorghum M-81E. The findings revealed that overexpression of SbMYBHv33 in sorghum significantly reduced sorghum biomass accumulation at the seedling stage and also salinity tolerance. Meanwhile, a heterologous transformation of Arabidopsis with SbMYBHv33 produced a similar phenotype. The loss of function of the Arabidopsis homolog of SbMYBHv33 resulted in longer roots and increased salt tolerance. Under normal conditions, SbMYBHV33 overexpression promoted the expression of ABA pathway genes in sorghum and inhibited growth. Under salt stress conditions, the gene expression of SbMYBHV33 decreased in the overexpressed lines, and the promotion of these genes in the ABA pathway was attenuated. This might be an important reason for the difference in growth and stress resistance between SbMYBHv33-overexpressed sorghum and ectopic expression Arabidopsis. Hence, SbMYBHv33 is an important component of sorghum growth and development and the regulation of salt stress response, and it could negatively regulate salt tolerance and biomass accumulation in sorghum.

The R2R3MYB transcription factors MaMYBF and MaMYB1 regulate flavonoid biosynthesis in grape hyacinth

R2R3 MYBs play vital roles in the regulation of flavonoid biosynthesis. However, the regulatory network of R2R3 MYBs in flavonoid biosynthesis is not fully understood in grape hyacinth (Muscari spp.). Here, we identified two R2R3 MYBs, MaMYBF and MaMYB1, as potential regulators of flavonol and anthocyanin biosynthesis, respectively. MaMYBF and MaMYB1 expression was elevated during flower development and was light-induced, and the expression patterns were related to those of the flavonoid structural genes MaFLS and MaDFR, respectively. The BiFC assay verified that MaMYB1 interacts with MabHLH1, but MaMYBF does not. A dual luciferase assay revealed that MaMYBF alone strongly activated pMaFLS, and its activation was attenuated at reduced doses of MaMYBF in the presence of MabHLH1, MaMybA, and MaMYB1. MaDFR transcription mediated by MaMybA and MabHLH1 was inhibited by MaMYB1. Moreover, overexpression of MaMYBF and MaMYB1 in tobacco reduced flower pigmentation and repressed the expression of flavonoid pathway key structural genes. Therefore, MaMYBF regulates the flavonol pathway independently of cofactors. Whereas MaMYB1 regulates anthocyanin biosynthesis by binding to MabHLH1 and disrupting the MaMybA-bHLH complex in grape hyacinth. Our results offer new insights into the intricate regulatory network of flavonoids in grape hyacinth involving the regulation of both flavonol and anthocyanin.

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.

Ablation of glucosinolate accumulation in the oil crop Camelina sativa by targeted mutagenesis of genes encoding the transporters GTR1 and GTR2 and regulators of biosynthesis MYB28 and MYB29

Camelina sativa is an oil crop with low input costs and resistance to abiotic and biotic stresses. The presence of glucosinolates, plant metabolites with adverse health effects, restricts the use of camelina for human and animal nutrition. Cas9 endonuclease-based targeted mutagenesis of the three homeologs of each of the glucosinolate transporters CsGTR1 and CsGTR2 caused a strong decrease in glucosinolate amounts, highlighting the power of this approach for inactivating multiple genes in a hexaploid crop. Mutagenesis of the three homeologs of each of the transcription factors CsMYB28 and CsMYB29 resulted in the complete loss of glucosinolates, representing the first glucosinolate-free Brassicaceae crop. The oil and protein contents and the fatty acid composition of the csgtr1csgtr2 and csmyb28csmyb29 mutant seeds were not affected. The decrease and elimination of glucosinolates improves the quality of the oil and press cake of camelina, which thus complies with international standards regulating glucosinolate levels for human consumption and animal feeding.

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.

Identification of MYB Transcription Factors Involving in Fruit Quality Regulation of Fragaria × ananassa Duch

The cultivated strawberry (Fragaria × ananassa Duch.) is an important horticultural crop. The economic values of strawberry cultivars are decided by their fruit qualities including taste, color and aroma. The important role of MYB transcription factors in fruit quality regulation is recognized increasingly with the identification of MYB genes involved in metabolism. A total of 407 MYB genes of F. × ananassa (FaMYBs) were identified in the genome-wide scale and named according to subgenome locations. The 407 FaMYBs were clustered into 36 groups based on phylogenetic analysis. According to synteny analysis, whole genome duplication and segmental duplication contributed over 90% of the expansion of the FaMYBs family. A total of 101 FaMYB loci with 1-6 alleles were identified by the homologous gene groups on homologous chromosomes. The differentially expressed FaMYB profiles of three cultivars with different fruit quality and fruit ripe processes provided the 8 candidate loci involved in fruit quality regulation. In this experiment, 7, 5, and 4 FaMYBs were screeded as candidate genes involved in the regulation of metabolism/transportation of anthocyanins, sugars or organic acids and 4-hydroxy-2, 5-dimethyl-3(2H)-furanone, respectively. These results pointed out the key FaMYBs for further functional analysis of gene regulation of strawberry fruit quality and would be helpful in the clarification on ofe roles of MYBs in the metabolism of fruit crops.

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.

The Anthocyanin Accumulation Related ZmBZ1, Facilitates Seedling Salinity Stress Tolerance via ROS Scavenging

Anthocyanins are a class of antioxidants that scavenge free radicals in cells and play an important role in promoting human health and preventing many diseases. Here, we characterized a maize Bronze gene (BZ1) from the purple colored W22 introgression line, which encodes an anthocyanin 3-O-glucosyltransferase, a key enzyme in the anthocyanin synthesis pathway. Mutation of ZmBZ1 showed bronze-colored seeds and reduced anthocyanins in seeds aleurone layer, seedlings coleoptile, and stem of mature plants by comparison with purple colored W22 (WT). Furthermore, we proved that maize BZ1 is an aleurone layer-specific expressed protein and sub-located in cell nucleus. Real-time tracing of the anthocyanins in developing seeds demonstrated that the pigment was visible from 16 DAP (day after pollination) in field condition, and first deposited in the crown part then spread all over the seed. Additionally, it was transferred along with the embryo cell activity during seed germination, from aleurone layer to cotyledon and coleoptile, as confirmed by microscopy and real-time qRT-PCR. Finally, we demonstrated that the ZmBZ1 contributes to stress tolerance, especially salinity. Further study proved that ZmBZ1 participates in reactive oxygen scavenging (ROS) by accumulating anthocyanins, thereby enhancing the tolerance to abiotic stress.

Silencing of SlMYB50 affects tolerance to drought and salt stress in tomato

High salinity and drought stresses often cause plants to produce ROS, including hydrogen peroxide (H₂O₂) and superoxide (O₂^-), which interfere with plant growth and affect crop yield. The transcription factors of the MYB family are involved in responses to biotic and abiotic stresses. Here, we isolated the R2R3-MYB transcription factor gene SlMYB50 and found that silencing of SlMYB50 increased resistance to PEG 6000, mannitol and salt. In addition, the resistance of transgenic tomatoes increased under high salt and drought stress. After stress treatment, the relative water content, chlorophyll content (critical for carbon fixation) and root vitality of the SlMYB50-RNAi lines were higher than those of the wild-type (WT). The opposite was true the water loss rate, relative conductivity, and MDA (as a sign of cell wall disruption). Under drought stress conditions, SlMYB50-silenced lines exhibited less H₂O₂ and less O₂^- accumulation, as well as higher CAT enzyme activity, than were exhibited by the WT. Notably, after stress treatment, the expression levels of chlorophyll-synthesis-related, flavonoid-synthesis-related, carotenoid-related, antioxidant-enzyme-related and ABA-biosynthesis-related genes were all upregulated in SlMYB50-silenced lines compared to those of WT. A dual-luciferase reporter system was used to verify that SlMYB50 could bind to the CHS1 promoter. In summary, this study identified essential roles for SlMYB50 in regulating drought and salt tolerance.

Genome editing of SlMYB3R3, a cell cycle transcription factor gene of tomato, induces elongated fruit shape

Fruit shape is an important trait that attracts consumers, and the regulation of genes related to cell division is crucial for shaping multicellular organs. In Arabidopsis, MYB3R transcription factors, which harbor three imperfect repeats in the N-terminus, control organ growth by regulating cell division. However, the function of MYB3Rs in tomato remains unknown. Here, we characterized tomato SlMYB3R3, which was preferentially expressed in flowers and placed in a subclade with two Arabidopsis cell cycle suppressors (MYB3R3/5). slmyb3r3 knockout mutants were generated using the CRISPR/Cas9 system. Morphological observation of the slmyb3r3 mutants showed that fruits that were elongated and occasionally peanut-like in shape were formed, which was caused by significantly increased cell numbers in the longitudinal direction. Transcriptome and yeast one-hybrid assay results suggested that SlMYB3R3 acted as a suppressor of cell-cycle-related genes by binding to the mitosis-specific activator (MSA) motifs in their promoters. Taken together, knock out of the suppressor SlMYB3R3 leads to elongated fruit, which results from the altered cell division pattern at the ovary stage, by regulating cell-cycle-related genes in an MSA-dependent manner. Our results suggest that SlMYB3R3 and its orthologs have the potential to change fruit shape as part of the molecular breeding of fruit crops.

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.

A Novel R2R3-MYB Transcription Factor FtMYB22 Negatively Regulates Salt and Drought Stress through ABA-Dependent Pathway

Tartary buckwheat (Fagopyrum tataricum Gaertn.) is a coarse cereal with strongly abiotic resistance. The MYB family plays a regulatory role in plant growth, development, and responses to biotic and abiotic stresses. However, the characteristics and regulatory mechanisms of MYB transcription factors in Tartary buckwheat remain unclarified. Here, this study cloned the FtMYB22 gene from Tartary buckwheat, and investigated its involvement in responding to individual water deficit and salt stress in Arabidopsis. Sequence analysis highlighted that the N-termini of FtMYB22 contained two highly conserved SANT domains and one conserved domain from the SG20 subfamily. Nucleus-localized FtMYB22 did not have individual transcriptional activation activity. Water deficiency and salt stress induced the high expression of the GUS gene, which was driven by the promoter of FtMYB22. Yeast stress experiments showed that the overexpression of FtMYB22 significantly reduced the growth activity of transgenic yeast under water deficit or salt stress. Consistently, the overexpression of FtMYB22 reduced the salt and water deficit stress resistance of the transgenic plants. In addition, physiological parameters showed that transgenic plants had lower proline and antioxidant enzyme activity under stress conditions. Compared to the wild-type (WT), transgenic plants accumulated more malondialdehyde (MDA), H2O2, and O2−; they also showed higher ion permeability and water loss rates of detached leaves under stress treatments. Notably, FtMYB22 was involved in plant stress resistance through an ABA-dependent pathway. Under stress conditions, the expression of RD29A, RD29B, PP2CA, KIN1, COR15A, and other genes in response to plant stress in transgenic lines was significantly lower than that in the WT (p < 0.05). Furthermore, yeast two-hybrid assay showed that there was a significant interaction between FtMYB22 and the ABA receptor protein RCAR1/2, which functioned in the ABA signal pathway. Altogether, FtMYB22, as a negative regulator, inhibited a variety of physiological and biochemical reactions, affected gene expression and stomatal closure in transgenic plants through the ABA-dependent pathway, and reduced the tolerance of transgenic Arabidopsis to water deficiency and salt stress. Based on these fundamental verifications, further studies would shed light on the hormone signal response mechanism of FtMYB22.

Genome-Scale Identification, Classification, and Expression Profiling of MYB Transcription Factor Genes in Cinnamomum camphora

The camphor tree (Cinnamomum camphora (L.) Presl.) is the representative species of subtropical evergreen broadleaved forests in eastern Asia and an important raw material for essential oil production worldwide. Although MYBs have been comprehensively characterized and their functions have been partially resolved in many plants, it has not been explored in C. camphora. In this study, 121 CcMYBs were identified on 12 chromosomes in the whole genome of C. camphora and found that CcMYBs were mainly expanded by segmental duplication. They were divided into 28 subgroups based on phylogenetic analysis and gene structural characteristics. In the promoter regions, numerous cis-acting elements were related to biological processes. Analysis of RNA sequencing data from seven tissues showed that CcMYBs exhibited different expression profiles, suggesting that they have various roles in camphor tree development. In addition, combined with the correlation analysis of structural genes in the flavonoid synthesis pathway, we identified CcMYBs from three subgroups that might be related to the flavonoid biosynthesis pathway. This study systematically analyzed CcMYBs in C. camphora, which will set the stage for subsequent research on the functions of CcMYBs during their lifetime and provide valuable insights for the genetic improvement of camphor trees.

A naturally-occurring phenomenon of flower color change during flower development in Xanthoceras sorbifolium

Xanthoceras sorbifolium (yellowhorn) is originated in China and is a unique tree in northern China. Yellowhorn is very popular because of the gradual color change of its flower at different flower developmental stages. During flowering development, the color at the base of yellowhorn flower petals gradually changes from yellow to purple. The mechanism of this miraculous phenomenon is still unclear. Here we show that anthocyanin accumulation during flowering development is the main reason for this color change. RT-PCR results show that the expression level of a variety of anthocyanin biosynthesis genes changes in different flower developmental stages. Realtime results show that the expression changes of these anthocyanin biosynthesis genes are positively regulated by a cluster of R2R3-MYB transcription factor genes, XsMYB113s. Furthermore, the DNA methylation analysis showed that CHH methylation status on the transposon element near the XsMYB113-1 influence its transcript level during flowering development. Our results suggest that dynamic epigenetic regulation of the XsMYB113-1 leads to the accumulation of anthocyanins during yellowhorn flower color change. These findings reemphasize the role of epigenetic regulation in flower development and provide a foundation for further studies of epigenetic regulation in long-lived woody perennials.

Glandular trichome specificity of menthol biosynthesis pathway gene promoters from Mentha × piperita

Several cis-elements including Myb-binding motifs together confer glandular trichome specificity as revealed from heterologous expression and analysis of menthol biosynthesis pathway gene promoters. Glandular Trichomes (GTs) are result of division of epidermal cells that produce diverse metabolites. Species of mint family are important for their essential oil containing many high-value terpenoids, biosynthesized and stored in these GTs. Hence, GTs constitute attractive targets for metabolic engineering and GT-specific promoters are important. In this investigation, the upstream regions of the Mentha × piperita menthol biosynthetic pathway genes (-)-limonene synthase, (-)-P450 limonene-3- hydroxylase, (-)-trans-isopiperitenol dehydrogenase, (-)-Isopiperitenone reductase, ( +)-Pulegone reductase, (-)-Menthone reductase/ (-)-Menthol dehydrogenase and a branched pathway gene ( +)-menthofuran synthase were isolated and characterized. These fragments, fused to β-glucuronidase (GUS) reporter gene of pBI101 binary vector, are able to drive high level gene expression in transgenic tobacco trichomes with strong signals in GTs, except for (-)-Isopiperitenone reductase. The GT-enriched tissue from transformed plants were analysed for GUS enzyme activity and RNA expression which correlates the GUS staining. To characterize the cis-elements responsible for GT-specific expression, a series of 5' deletion constructs for MpPLS and MpPMFS were cloned and analysed in stable transgenic tobacco lines. The specificity of trichome expression was located to -  797 to-  598 bp sequence for (-)-limonene synthase and-  629 to -   530 bp for ( +)-menthofuran synthase promoters containing specific Myb-binding motifs in addition to other unique motifs described for developmental regulation without any defined pattern. All other pathway promoters also recruits specific but different Myb factors as indicated by this analysis.

PeMYB4L interacts with PeMYC4 to regulate anthocyanin biosynthesis in Phalaenopsis orchid

Phalaenopsis spp., one genus of Orchidaceae, have become very popular worldwide for their fascinating flowers with various colors and pigmentation patterns. Several R2R3-MYB transcription factors have been reported to function in anthocyanin accumulation in Phalaenopsis spp. However, its molecular mechanism underlying the detailed regulatory pathway remains poorly understood. In this study, we identified a novel subgroup 2 R2R3-MYB transcription factor PeMYB4L, the expression profile of which was concomitant with red color formation in Phalaenopsis spp. flowers. Virus-induced gene silencing (VIGS) and transient overexpression assay verified that PeMYB4L promotes anthocyanin accumulation in flower tissues. In addition, PeMYB4L could directly regulates the expression of Phalaenopsis spp. chalcone synthase gene (PeCHS) through MYBST1 (GGATA) binding site. It's interesting that the basic-helix-loop-helix (bHLH) protein PeMYC4 shows opposite expression pattern from PeMYB4L in anthocyanin accumulation. Furthermore, PeMYC4 was verified to form MYB-bHLH complex with PeMYB4L, and attenuated the expression of PeCHS and weakened anthocyanin production, indicating a novel regulatory model of MYB-bHLH complex. Our findings uncover the detailed regulatory pathway of MYB-bHLH, and might provide a new insight into the complicated anthocyanin pigmentation in Phalaenopsis spp.