MYB Transcription Factor Family Research Articles

The MYB Transcription Factor Family in Eucommia ulmoides: Genome-Wide Identification, Characterization, and Network Analysis in Relation to the Rubber Biosynthetic Genes

The MYB Transcription Factor Family in Eucommia ulmoides: Genome-Wide Identification, Characterization, and Network Analysis in Relation to the Rubber Biosynthetic Genes

Physiological and transcriptomic analysis of the mangrove species Kandelia obovata in response to flooding stress

Flooding stress on mangroves is growing continually with rising sea level. In this study, the physiology and transcriptome of the mangrove species Kandelia obovata under flooding stress were analyzed. With increasing inundation time, malondialdehyde (MDA), superoxide dismutase (SOD), glutathione (GSH), soluble sugar (SS), soluble protein (SP), and proline (Pro) content declined, while peroxidase (POD) and ascorbate peroxidase (APX) activity rose significantly. According to the KEGG pathway enrichment analysis, upregulated differentially expressed genes (DEGs) were enriched in the plant hormone signaling pathway. Furthermore, MYB44 and MYB108 genes from the MYB transcription factor family and RAP2.12, DREB2B, and ERF4 genes from the AP2/ERF family were up-regulated under flooding conditions. A strong correlation was established between the expression levels of 12 DEGs under flooding stress and RNA sequencing data and was verified by qRT-PCR. These results provide new insights into the molecular mechanism of K. obovata in response to flooding stress.

Integrated metabolic profiling and transcriptome analysis of Lonicera japonica flowers for chlorogenic acid, luteolin and endogenous hormone syntheses

The active ingredients of many medicinal plants are the secondary metabolites associated with the growth period. Lonicera japonica Thunb. is an important traditional Chinese medicine, and the flower development stage is an important factor that influences the quality of medicinal ingredients. In this study, transcriptomics and metabolomics were performed to reveal the regulatory mechanism of secondary metabolites during flowering of L. japonica. The results showed that the content of chlorogenic acid (CGA) and luteolin gradually decreased from green bud stage (Sa) to white flower stage (Sc), especially from white flower bud stage (Sb) to Sc. Most of the genes encoding the crucial rate-limiting enzymes, including PAL, C4H, HCT, C3′H, F3′H and FNSII, were down-regulated in three comparisons. Correlation analysis identified some members of the MYB, AP2/ERF, bHLH and NAC transcription factor families that are closely related to CGA and luteolin biosynthesis. Furthermore, differentially expressed genes (DEGs) involved in hormone biosynthesis, signalling pathways and flowering process were analysed in three flower developmental stage.

Interactive effects of light quality and nitrate supply on growth and metabolic processes in two lettuce cultivars (Lactuca sativa L.)

Light quality and nutrient parameters alter plant growth and metabolic pathways in leafy greens, therefore being primarily optimized in indoor precision agriculture. While their individual effects were studied extensively with large-scale phenomics and transcriptomics, the combinatorial effects remain poorly understood. Here, we explored the hidden interactions among light- and nitrate-responsive transcriptomes and their associations with six growth and biochemical traits, using the full factorial experimental design with three lights (red (R), blue (B) and a 1:1 mixture of red and blue (RB)), two nitrate concentrations and two cultivars (Green oak-leaf and Red oak-leaf). Multifactor analysis of 36 RNA-sequencing data revealed their largely independent but additive effects that are associated with MYB, bHLH and GATA family transcription factors. Red light and sufficient nitrate supply additively increased the primary metabolite pathway and bHLH and GATA genes, leading to a high growth rate. Blue and RB lights, together with nitrate deficiency, induced flavonoid pathway and MYB genes. Variety-specific deviations were detected in late anthocyanin biosynthetic genes and specific MYB transcription factors that would explain anthocyanin accumulations in the Red oak-leaf cultivar. Taken together, our large-scale transcriptomic and physiologic results suggest that the two environmental parameters (i.e., light quality and nitrate level) (1) utilize independent transcriptional regulatory mechanisms to achieve similar phenotypic and metabolic changes, and (2) induce the variety-specific anthocyanin accumulations via a small number of genes that show different gene expression profiles between the two cultivars.

In-silico mining and characterization of MYB family genes in wilt-resistant hybrid guava (Psidium guajava × Psidium molle)

BackgroundThe MYB family is one of the most significant groups of transcription factors in plants. However, several MYBs have been linked to secondary metabolism and are important for determining the color of fruit’s peel and pulp. Despite being a substantial fruit crop in tropical and subtropical areas of the world, wilt-resistant hybrid guava (Psidium guajava × Psidium molle; PGPM) has not yet been the subject of a thorough examination. This study’s goal was to assess the expression of MYB in guava fruit pulp, roots, and seeds to predict its function by in silico analysis of the guava root transcriptome data. ResultsIn the current study, we have mined the MYBs family of MYB genes from the transcriptome of the PGPM guava root. We have mined 15 distinct MYB transcription factor genes/transcripts viz MYB3, MYB4, MYB23, MYB86, MYB90, MYB308, MYB5, MYB82, MYB114, MYB6, MYB305, MYB44, MYB51, MYB46, and MYB330. From the analyses, it was found that R2-MYB and R3-MYB domains are conserved in all known guava MYB proteins. The expression of six different MYB TFs was examined using semi-quantitative RT-PCR in “Shweta” pulp (white colour pulp), “Lalit” pulp (red color pulp), “Lalit” root, and “Lalit” seed. ConclusionThere were 15 MYB family members observed in guava. They were unequally distributed across the chromosomes, most likely as a result of gene duplication. Additionally, the expression patterns of the particular MYBs showed that MYB may be involved in the control of wilt, fruit ripening, seed development, and root development. Our results allow for a more thorough functional characterization of the guava MYB family genes and open the door to additional research into one essential MYB transcription factor family of genes and its involvement in the growth and ripening of guava fruit.

MYB112 connects light and circadian clock signals to promote hypocotyl elongation in Arabidopsis.

Ambient light and the endogenous circadian clock play key roles in regulating Arabidopsis (Arabidopsis thaliana) seedling photomorphogenesis. PHYTOCHROME-INTERACTING FACTOR 4 (PIF4) acts downstream of both light and the circadian clock to promote hypocotyl elongation. Several members of the R2R3-MYB transcription factor (TF) family, the most common type of MYB TF family in Arabidopsis, have been shown to be involved in regulating photomorphogenesis. Nonetheless, whether R2R3-MYB TFs are involved in connecting the light and clock signaling pathways during seedling photomorphogenesis remains unknown. Here, we report that MYB112, a member of the R2R3-MYB family, acts as a negative regulator of seedling photomorphogenesis in Arabidopsis. The light signal promotes the transcription and protein accumulation of MYB112. myb112 mutants exhibit short hypocotyls in both constant light and diurnal cycles. MYB112 physically interacts with PIF4 to enhance the transcription of PIF4 target genes involved in the auxin pathway, including YUCCA8 (YUC8), INDOLE-3-ACETIC ACID INDUCIBLE 19 (IAA19), and IAA29. Furthermore, MYB112 directly binds to the promoter of LUX ARRHYTHMO (LUX), the central component of clock oscillators, to repress its expression mainly in the afternoon and relieve LUX-inhibited expression of PIF4. Genetic evidence confirms that LUX acts downstream of MYB112 in regulating hypocotyl elongation. Thus, the enhanced transcript accumulation and transcriptional activation activity of PIF4 by MYB112 additively promotes the expression of auxin-related genes, thereby increasing auxin synthesis and signaling and fine-tuning hypocotyl growth under diurnal cycles.

High-throughput sgRNA testing reveals rules for Cas9 specificity and DNA repair in tomato cells.

CRISPR/Cas9 technology has the potential to significantly enhance plant breeding. To determine the specificity and the mutagenic spectrum of SpCas9 in tomato, we designed 89g(uide) RNAs targeting genes of the tomato MYB transcription factor family with varying predicted specificities. Plasmids encoding sgRNAs and Cas9 were introduced into tomato protoplasts, and target sites as well as 224 predicted off-target sites were screened for the occurrence of mutations using amplicon sequencing. Algorithms for the prediction of efficacy of the sgRNAs had little predictive power in this system. The analysis of mutations suggested predictable identity of single base insertions. Off-target mutations were found for 13 out of 89 sgRNAs and only occurred at positions with one or two mismatches (at 14 and 3 sites, respectively). We found that PAM-proximal mismatches do not preclude low frequency off-target mutations. Off-target mutations were not found at all 138 positions that had three or four mismatches. We compared off-target mutation frequencies obtained with plasmid encoding sgRNAs and Cas9 with those induced by ribonucleoprotein (RNP) transfections. The use of RNPs led to a significant decrease in relative off-target frequencies at 6 out of 17, no significant difference at 9, and an increase at 2 sites. Additionally, we show that off-target sequences with insertions or deletions relative to the sgRNA may be mutated, and should be considered during sgRNA design. Altogether, our data help sgRNA design by providing insight into the Cas9-induced double-strand break repair outcomes and the occurrence of off-target mutations.

Weighted gene coexpression correlation network analysis reveals the potential molecular regulatory mechanism of citrate and anthocyanin accumulation between postharvest ‘Bingtangcheng’ and ‘Tarocco’ blood orange fruit

BackgroundOrganic acids and anthocyanins are the most important compounds for the flavor and nutritional quality of citrus fruit. However, there are few reports on the involvement of co-regulation of citrate and anthocyanin metabolism. Here, we performed a comparative transcriptome analysis to elucidate the genes and pathways involved in both citrate and anthocyanin accumulation in postharvest citrus fruit with ‘Tarocco’ blood orange (TBO; high accumulation) and ‘Bingtangcheng’ sweet orange (BTSO; low accumulation).ResultsA robust core set of 825 DEGs were found to be temporally associated with citrate and anthocyanin accumulation throughout the storage period through transcriptome analysis. Further according to the results of weighted gene coexpression correlation network analysis (WGCNA), the turquoise and brown module was highly positively correlated with both of the content of citrate and anthocyanin, and p-type ATPase (PH8), phosphoenolpyruvate carboxylase kinase (PEPCK), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavonoid 3’-hydroxylase (F3’H) and glutathione S transferase (GST) were considered key structural genes. Moreover, MYB family transcription factor (PH4), Zinc finger PHD-type transcription factor (CHR4, HAC12), Zinc finger SWIM-type transcription factor (FAR1) and Zinc finger C3H1-type transcription factor (ATC3H64) were considered hub genes related to these structural genes. Further qRT-PCR analysis verified that these transcription factors were highly expressed in TBO fruit and their expression profiles were significantly positively correlated with the structural genes of citrate and anthocyanin metabolism as well as the content of citrate and anthocyanin content.ConclusionsThe findings suggest that the CHR4, FAR1, ATC3H64 and HAC12 may be the new transcription regulators participate in controlling the level of citrate and anthocyanin in postharvest TBO fruit in addition to PH4. These results may providing new insight into the regulation mechanism of citrate and anthocyanin accumulation in citrus fruit.

A Medicago truncatula Autoregulation of Nodulation Mutant Transcriptome Analysis Reveals Disruption of the SUNN Pathway Causes Constitutive Expression Changes in Some Genes, but Overall Response to Rhizobia Resembles Wild-Type, Including Induction of TML1 and TML2

Nodule number regulation in legumes is controlled by a feedback loop that integrates nutrient and rhizobia symbiont status signals to regulate nodule development. Signals from the roots are perceived by shoot receptors, including a CLV1-like receptor-like kinase known as SUNN in Medicago truncatula. In the absence of functional SUNN, the autoregulation feedback loop is disrupted, resulting in hypernodulation. To elucidate early autoregulation mechanisms disrupted in SUNN mutants, we searched for genes with altered expression in the loss-of-function sunn-4 mutant and included the rdn1-2 autoregulation mutant for comparison. We identified constitutively altered expression of small groups of genes in sunn-4 roots and in sunn-4 shoots. All genes with verified roles in nodulation that were induced in wild-type roots during the establishment of nodules were also induced in sunn-4, including autoregulation genes TML2 and TML1. Only an isoflavone-7-O-methyltransferase gene was induced in response to rhizobia in wild-type roots but not induced in sunn-4. In shoot tissues of wild-type, eight rhizobia-responsive genes were identified, including a MYB family transcription factor gene that remained at a baseline level in sunn-4; three genes were induced by rhizobia in shoots of sunn-4 but not wild-type. We cataloged the temporal induction profiles of many small secreted peptide (MtSSP) genes in nodulating root tissues, encompassing members of twenty-four peptide families, including the CLE and IRON MAN families. The discovery that expression of TML2 in roots, a key factor in inhibiting nodulation in response to autoregulation signals, is also triggered in sunn-4 in the section of roots analyzed, suggests that the mechanism of TML regulation of nodulation in M. truncatula may be more complex than published models.

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

FvMYB44, a Strawberry R2R3-MYB Transcription Factor, Improved Salt and Cold Stress Tolerance in Transgenic Arabidopsis

MYB (v-MYB avian myeloblastosis viral oncogene homolog) transcription factor (TF) family has numerous members with complex and diverse functions, which perform an integral role in regulating the plant’s response to adversity. This study used cloning to obtain a novel MYB TF gene from the diploid strawberry Fragaria vesca, which was given the designation FvMYB44. Subcellular localization results showed that the protein of FvMYB44 was a nuclear localization protein. The resistance of Arabidopsis thaliana to salt and low temperature stresses was greatly enhanced by the overexpression of FvMYB44. When subjected to salt and temperature stress, transgenic plants showed higher proline and chlorophyll concentrations and higher superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) activities than wild-type (WT) and unloaded line (UL) of A. thaliana. In contrast, WT and UL lines had higher malondialdehyde (MDA) content and reactive oxygen species ROS (O2− and H2O2) content. These findings suggest that FvMYB44 may perform a role in controlling the response of A. thaliana to cold and salt stress.

Mining and identification of members of MYB transcription factor family in Lonicera macranthoides

As a large family of transcription factors, the MYB family plays a vital role in regulating flower development. We studied the MYB family members in Lonicera macranthoides for the first time and identified three sequences of 1R-MYB, 47 sequences of R2R3-MYB, two sequences of 3R-MYB, and one sequence of 4R-MYB from the transcriptome data. Further, their physicochemical properties, conserved domains, phylogenetic relationship, protein structure, functional information, and expression were analyzed. The results show that the 53 MYB transcription factors had different conserved motifs, physicochemical properties, structures, and functions in wild type and 'Xianglei' cultivar of L. macranthoides, indicating their conservation and diversity in evolution. The transcript level of LmMYB was significantly different between the wild type and 'Xianglei' cultivar as well as between flowers and leaves, and some genes were specifically expressed. Forty-three out of 53 LmMYB sequences were expressed in both flowers and leaves, and 9 of the LmMYB members showed significantly different transcript levels between the wild type and 'Xianglei' cultivar, which were up-regulated in the wild type. The results provide a theoretical basis for further studying the specific functional mechanism of the MYB family.

Comprehensive Transcriptome Analysis of Responses during Cold Stress in Wheat (Triticum aestivum L.).

Wheat production is often impacted by pre-winter freezing damage and cold spells in later spring. To study the influences of cold stress on wheat seedlings, unstressed Jing 841 was sampled once at the seedling stage, followed by 4 °C stress treatment for 30 days and once every 10 days. A total of 12,926 differentially expressed genes (DEGs) were identified from the transcriptome. K-means cluster analysis found a group of genes related to the glutamate metabolism pathway, and many genes belonging to the bHLH, MYB, NAC, WRKY, and ERF transcription factor families were highly expressed. Starch and sucrose metabolism, glutathione metabolism, and plant hormone signal transduction pathways were found. Weighted Gene Co-Expression Network Analysis (WGCNA) identified several key genes involved in the development of seedlings under cold stress. The cluster tree diagram showed seven different modules marked with different colors. The blue module had the highest correlation coefficient for the samples treated with cold stress for 30 days, and most genes in this module were rich in glutathione metabolism (ko00480). A total of eight DEGs were validated using quantitative real-time PCR. Overall, this study provides new insights into the physiological metabolic pathways and gene changes in a cold stress transcriptome, and it has a potential significance for improving freezing tolerance in wheat.

Genome-Wide Identification and Expression Analysis of the R2R3-MYB Gene Family in Rubber Trees

The plant MYB transcription factor family featured as highly conserved DNA-binding domains consisting of 1 to 4 imperfect repeats (R). Increasing evidence indicates that MYB genes participates in growth, differentiation, metabolism, and biotic and abiotic stress responses. However, the functions of MYB genes in the rubber tree remain to be deeply elucidated, especially R2R3-MYB gene family. In this study, molecular biology, bioinformatics, and qRT-PCR were used to identify and analyze HbR2R3-MYB gene family members in the rubber tree. A total of 132 members of the R2R3-MYB gene family were identified in the rubber tree based on genome-wide level. Most of the HbR2R3-MYBs were mapped to 17 rubber tree chromosomes except four genes. A phylogenetic analysis divided all the HbR2R3-MYBs into 20 subgroups with Arabidopsis thaliana. MEME analysis showed that the protein of HbR2R3-MYBs was characterized by 9 conserved motifs. Twenty-six representative R2R3 HbMYBs from different subgroups were selected for expression profiles analysis and the results revealed that the HbR2R3-MYBs members showed various expression patterns in different tissues, powdery mildew-infected and ethylene treatment, implying the diversity of their functions in rubber trees. These results provide fundamental knowledge for further studying the response of the HbR2R3-MYB family to stress and regulation latex flow in rubber tree.

Combined transcriptomic and metabolomic analyses reveal the mechanism of debagged ‘Fuji’ apple sunburn

Bagging of apples can improve the quality and commercial value of fruit. However, the sudden exposure of bagged fruit to high temperature and high light can lead to sunburn of the fruit surface. The molecular mechanism underlying the process of apple skin sunburn is poorly characterized. Here, we analyzed transcriptomic and metabolomic data from debagged apples subjected to high temperature and high light stress. Specifically, we studied the differential gene expression and metabolite accumulation in normal apple skin (NM), photooxidized apple skin (PS), and skin with browned (SB). High light and heat stress led to enrichment of differentially expressed genes (DEGs) related to the flavonoid biosynthesis pathway, and screening revealed 32 transcription factors (TFs) and 11 structural genes as targets. Furthermore, we found that MYB family TFs were significantly correlated with most structural genes in the flavonoid metabolism pathway, particularly the MdMYBR9 TF. We injected the MdMYBR9 TF into apple skin, and we found that the apple skin with overexpression of MdMYBR9 TF exhibited less severe sunburn, the apple with silenced MdMYBR9 TF exhibited severe sunburn. These findings provide a theoretical basis for further understanding the transcriptional regulation mechanism related to apple skin sunburn.

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

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

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.

RRS1 shapes robust root system to enhance drought resistance in rice.

A strong root system facilitates the absorption of water and nutrients from the soil, to improve the growth of crops. However, to date, there are still very few root development regulatory genes that can be used in crop breeding for agriculture. In this study, we cloned a negative regulator gene of root development, Robust Root System 1 (RRS1), which encodes an R2R3-type MYB family transcription factor. RRS1 knockout plants showed enhanced root growth, including longer root length, longer lateral root length, and larger lateral root density. RRS1 represses root development by directly activating the expression of OsIAA3 which is involved in the auxin signaling pathway. A natural variation in the coding region of RRS1 changes the transcriptional activity of its protein. RRS1T allele, originating from wild rice, possibly increases root length by means of weakening regulation of OsIAA3. Knockout of RRS1 enhances drought resistance by promoting water absorption and improving water use efficiency. This study provides a new gene resource for improving root systems and cultivating drought-resistant rice varieties with important values in agricultural applications.

LiMYB108 is involved in floral monoterpene biosynthesis induced by light intensity in Lilium 'Siberia'.

We find that the MYB family transcription factor, LiMYB108, has a novel function to regulate the floral fragrance affected by light intensity. Floral fragrance determines the commercial value of flowers and is influenced by many environmental factors, especially light intensity. However, the mechanism by which light intensity affects the release of floral fragrance is unclear. Here, we isolated an R2R3-type MYB transcription factor LiMYB108, the expression of which was induced by light intensity and located in the nucleus. Light of 200 and 600μmolm-1s-1 significantly increased the expression of LiMYB108, which was consistent with the improving trend of monoterpene synthesis under light. Virus-induced gene silencing (VIGS) of LiMYB108 in Lilium not only significantly inhibited the synthesis of ocimene and linalool, but also decreased the expression of LoTPS1; however, transient overexpression of LiMYB108 exerted opposite effects. Furthermore, yeast one-hybrid assays, dual-luciferase assays, and electrophoretic mobility shift assays (EMSA) demonstrated that LiMYB108 directly activated the expression of LoTPS1 by binding to the MYB binding site (MBS) (CAGTTG). Our findings demonstrate that light intensity triggered the high expression of LiMYB108, and then LiMYB108 as a transcription factor to activate the expression of LoTPS1, thus promoting the synthesis of the ocimene and linalool, which are important components of floral fragrance. These results provide new insights into the effects of light intensity on floral fragrance synthesis.

Transcriptome Analysis Reveals the Central Role of the Transcription Factor MYB in Regulating Anthocyanin Accumulation in Economic Grape Species (Vitis vinifera)

To cultivate different grape varieties according to market needs, it is necessary to study the regulation mechanism of color changes in different development stages of grapes. In this study, RNA-sequencing (RNA-Seq) technology was used to compare and analyze the transcriptome data of four grape varieties at the same development stage. Among the annotated differential genes, the anthocyanin synthesis pathway in the flavonoid pathway was mainly studied. Further RT-qPCR analysis of key enzyme genes, in the flavonoid synthesis pathway of the anthocyanin metabolism pathway, showed that the MYB transcription factor family had binding sites at the start of the four enzyme genes. The relative expression of the MYB transcription factor and enzyme gene in the transcriptome data was verified by reverse transcription polymerase chain reaction. Subcellular localization and gene function verification of the transcription factor MYB2 confirmed its regulatory role in anthocyanins.