Genome-wide identification of Pistacia R2R3-MYB gene family and function characterization of PcMYB113 during autumn leaf coloration in Pistacia chinensis

Molecular mapping of candidate genes in determining red color of perilla leaf

Perilla frutescens is a traditional medicinal plant and functional food in Asian communities, characterized by distinct red and green leaf types that have significant phenotypic and medicinal implications. However, the genetic mechanisms controlling anthocyanin synthesis in this species remain unclear. Genetic analysis serves as a powerful tool for investigating the pivotal genes and regulatory mechanisms governing anthocyanin accumulation in red and green perilla. In this study, an F2 segregation population was constructed from a hybrid of red and green perilla, and representative samples were subjected to mix-sequencing using BSA-seq and BSR-seq. A 6.0 Mb candidate region on chromosome 8 was identified, pinpointing PfMYB113b, PfC4H1, and PfF3H as key genes involved in anthocyanin biosynthesis. The insertion of a repeat sequence in the promoter of PfMYB113b leads to alterations in gene expression levels. Furthermore, PfMYB113b regulates the transcription of PfC4H1 and PfF3H, thereby influencing anthocyanin synthesis. These findings enhance our understanding of the genetic regulatory mechanisms underlying leaf coloration in perilla.

Lipid metabolism during seed germination of Pistacia chinensis and its response to gibberellic acid

Lipid metabolism may play a critical role in fueling seed germination, but the knowledge of lipid metabolism during germination is still ambiguous. Here, we hypothesize that gibberellic acid (GA) promotes germination by means of enhancing lipid mobilization in Chinese pistachio (Pistacia chinensis Bunge), a species belonging to Anacardiaceae with high oil content in its seeds. A multi-omics approach has been applied to measure lipid mobilization during seed germination, and to identify the key regulators involved in GA-mediated lipid metabolism. The results indicated that GA contents increased, while IAA, ABA and JA contents decreased during seed germination. GA3 increased significantly in the two germination stages (i.e. imbibition and radicle protrusion), and it was more abundant than GA1 and GA4. In addition, the relative content of most lipids decreased during germination, and the differentially changed metabolites were significantly enriched in lipid metabolic pathways based on KEGG analysis. WGCNA indicated that GA3 was correlated with more genes in lipid metabolic pathways. Transcriptomic analysis further revealed that differentially expressed genes (DEGs) related to fatty acid biosynthesis, glycerolipid metabolism, glycerophospholipid metabolism and starch and sucrose metabolism were upregulated under GA3 application, such as the acetyl-CoA carboxylase biotin carboxyl carrier protein (ACCB), fatty acyl-ACP thioesterase B (FATB), diacylglycerol acyltransferase (DGAT) and DEFECTIVE IN ANTHER DEHISCENCE 1 (DAD1). Therefore, our study supports the hypothesis that GA promotes seed germination in P. chinensis by enhancing lipid mobilization. This study proposes a novel mechanism of lipid responses to exogenous GA, which contributes to a deep understanding of germination of oleaginous seeds.

MYB transcription factors in Peucedanum Praeruptorum Dunn: the diverse roles of the R2R3-MYB subfamily in mediating coumarin biosynthesis

BACKGROUND

The MYB superfamily (v-myb avian myeloblastosis viral oncogene homolog) plays a role in plant growth and development, environmental stress defense, and synthesis of secondary metabolites. Little is known about the regulatory function of MYB genes in Peucedanum praeruptorum Dunn, although many MYB family members, especially R2R3-MYB genes, have been extensively studied in model plants.

RESULTS

A total of 157 R2R3-MYB transcription factors from P. praeruptorum were identified using bioinformatics analysis. Comprehensive analyses including chromosome location, microsynteny, gene structure, conserved motif, phylogenetic tree, and conserved domain were further performed. The length of the 157 transcription factors ranged from 120 to 1,688 amino acids (molecular weight between 14.21 and 182.69 kDa). All proteins were hydrophilic. Subcellular localization predictions showed that 155 PpMYB proteins were localized in the nucleus, with PpMYB12 and PpMYB157 localized in the chloroplasts and mitochondria, respectively. Ten conserved motifs were identified in the PpMYBs, all of which contained typical MYB domains. Transcriptome analysis identified 47,902 unigenes. Kyoto Encyclopedia of Genes and Genomes analysis revealed 136 pathways, of which 524 genes were associated with the phenylpropanoid pathway. Differential expressed genes (DEGs) before and after bolting showed that 11 genes were enriched in the phenylpropanoid pathway. Moreover, the expression patterns of transcription genes were further verified by qRT-PCR. With high-performance liquid chromatography (HPLC), 8 coumarins were quantified from the root, stem, and leaf tissue samples of P. praeruptorum at different stages. Praeruptorin A was found in both roots and leaves before bolting, whereas praeruptorin B was mainly concentrated in the roots, and the content of both decreased in the roots and stems after bolting. Praeruptorin E content was highest in the leaves and increased with plant growth. The correlation analysis between transcription factors and coumarin content showed that the expression patterns of PpMYB3 and PpMYB103 in roots align with the accumulation trends of praeruptorin A, praeruptorin B, praeruptorin E, scopoletin, and isoscopoletin, which declined in content after bolting, suggesting that these genes may positively regulate the biosynthesis of coumarins. Eleven distinct metabolites and 48 DEGs were identified. Correlation analysis revealed that the expression of all DEGs were significantly related to the accumulation of coumarin metabolites, indicating that these genes are involved in the regulation of coumarin biosynthesis.

CONCLUSIONS

R2R3-MYB transcription factors may be involved in the synthesis of coumarin. Our findings provide basic data and a rationale for future an in-depth studies on the role of R2R3-MYB transcription factors in the growth and regulation of coumarin synthesis.

Catechin promotes the germination of Pistacia chinensis seeds via GA biosynthesis

BACKGROUND AND AIMS

Chinese pistachio (Pistacia chinensis), an important horticultural plant species, holds great ornamental value with beautiful leaves and fruits. Seedling propagation of this tree species is restricted by its erratic seed germination; however, the germination mechanism is ambiguous. The aim of this study was to determine the germination mechanism from a novel perspective based on the multi-omics data.

METHODS

The multi-omics technique combined with hormone content measurement was applied to seed germination of Chinese pistachio.

KEY RESULTS

Due to its great accumulation during seed germination, catechin stood out from the identified metabolites in a broadly targeted metabolomic analysis. Exogenous catechin at 10 mg L-1 significantly improved the germination of Chinese pistachio seeds. An interesting result of hormone analysis was that the improving effect of catechin could be attributed to an increase in gibberellic acid 3 (GA3) content rather than a decrease in abscisic acid (ABA) content before germination. Treatments with paclobutrazol (PAC, a GA biosynthesis inhibitor) and PAC + catechin also showed that the promoting effect of catechin on seed germination depends on GA biosynthesis. Transcriptome analysis and qRT‒PCR further revealed that catechin induced the expression of PcGA20ox5 to activate GA biosynthesis. Several transcription factors were induced by catechin and GA treatments, such as TCP, bZIP and C3H, which may play an important regulatory role in GA biosynthesis in a catechin-mediated way.

CONCLUSIONS

Catechin promotes seed germination via GA biosynthesis in Chinese pistachios. This study proposes a novel mechanism by which catechin promotes seed germination via the GA pathway, which provides new insight into a comprehensive understanding of seed dormancy and germination.

For a Colorful Life: Recent Advances in Anthocyanin Biosynthesis during Leaf Senescence

Leaf senescence is the last stage of leaf development, and it is accompanied by a leaf color change. In some species, anthocyanins are accumulated during leaf senescence, which are vital indicators for both ornamental and commercial value. Therefore, it is essential to understand the molecular mechanism of anthocyanin accumulation during leaf senescence, which would provide new insight into autumn coloration and molecular breeding for more colorful plants. Anthocyanin accumulation is a surprisingly complex process, and significant advances have been made in the past decades. In this review, we focused on leaf coloration during senescence. We emphatically discussed several networks linked to genetic, hormonal, environmental, and nutritional factors in regulating anthocyanin accumulation during leaf senescence. This paper aims to provide a regulatory model for leaf coloration and to put forward some prospects for future development.

The RING zinc finger protein LbRZF1 promotes salt gland development and salt tolerance in Limonium bicolor

The recretohalophyte Limonium bicolor thrives in high-salinity environments because salt glands on the above-ground parts of the plant help to expel excess salt. Here, we characterize a nucleus-localized C3HC4 (RING-HC)-type zinc finger protein of L. bicolor named  RING  ZINC  FINGER PROTEIN  1 (LbRZF1). LbRZF1 was expressed in salt glands and in response to NaCl treatment. LbRZF1 showed no E3 ubiquitin ligase activity. The phenotypes of overexpression and knockout lines for LbRZF1 indicated that LbRZF1 positively regulated salt gland development and salt tolerance in L. bicolor. lbrzf1 mutants had fewer salt glands and secreted less salt than did the wild-type, whereas LbRZF1-overexpressing lines had opposite phenotypes, in keeping with the overall salt tolerance of these plants. A yeast two-hybrid screen revealed that LbRZF1 interacted with LbCATALASE2 (LbCAT2) and the transcription factor LbMYB113, leading to their stabilization. Silencing of LbCAT2 or LbMYB113 decreased salt gland density and salt tolerance. The heterologous expression of LbRZF1 in Arabidopsis thaliana conferred salt tolerance to this non-halophyte. We also identified the transcription factor LbMYB48 as an upstream regulator of LbRZF1 transcription. The study of LbRZF1 in the regulation network of salt gland development also provides a good foundation for transforming crops and improving their salt resistance.

Abscisic acid-induced transcription factor PsMYB306 negatively regulates tree peony bud dormancy release

Bud dormancy is a crucial strategy for perennial plants to withstand adverse winter conditions. However, the regulatory mechanism of bud dormancy in tree peony (Paeonia suffruticosa) remains largely unknown. Here, we observed dramatically reduced and increased accumulation of abscisic acid (ABA) and bioactive gibberellins (GAs) GA1 and GA3, respectively, during bud endodormancy release of tree peony under prolonged chilling treatment. An Illumina RNA sequencing study was performed to identify potential genes involved in the bud endodormancy regulation in tree peony. Correlation matrix, principal component, and interaction network analyses identified a downregulated MYB transcription factor gene, PsMYB306, the expression of which positively correlated with 9-CIS-EPOXYCAROTENOID DIOXYGENASE 3 (PsNCED3) expression. Protein modeling analysis revealed 4 residues within the R2R3 domain of PsMYB306 to possess DNA binding capability. Transcription of PsMYB306 was increased by ABA treatment. Overexpression of PsMYB306 in petunia (Petunia hybrida) inhibited seed germination and plant growth, concomitant with elevated ABA and decreased GA contents. Silencing of PsMYB306 accelerated cold-triggered tree peony bud burst and influenced the production of ABA and GAs and the expression of their biosynthetic genes. ABA application reduced bud dormancy release and transcription of ENT-KAURENOIC ACID OXIDASE 1 (PsKAO1), GA20-OXIDASE 1 (PsGA20ox1), and GA3-OXIDASE 1 (PsGA3ox1) associated with GA biosynthesis in PsMYB306-silenced buds. In vivo and in vitro binding assays confirmed that PsMYB306 specifically transactivated the promoter of PsNCED3. Silencing of PsNCED3 also promoted bud break and growth. Altogether, our findings suggest that PsMYB306 negatively modulates cold-induced bud endodormancy release by regulating ABA production.

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

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

PSEUDO-ETIOLATION IN LIGHT proteins reduce greening by binding GLK transcription factors

Knocking out genes encoding proteins that downregulate the accumulation of pigments may lead to increases in crop quality and yield. PSEUDO-ETIOLATION IN LIGHT 1 (PEL1) downregulates the accumulation of carotenoids in carrot and chlorophyll in Arabidopsis and rice and may inhibit GOLDEN 2-LIKE (GLK) transcription factors. PEL1 belongs to a previously unstudied gene family found only in plants. We used CRISPR/Cas9 technology to knock out each member of the 4-member PEL gene family and both GLK genes in Arabidopsis. In pel mutants, chlorophyll levels were elevated in seedlings; after flowering, chloroplasts increased in size, and anthocyanin levels increased. Although the chlorophyll-deficient phenotype of glk1 glk2 was epistatic to pel1 pel2 pel3 pel4 in most of our experiments, glk1 glk2 was not epistatic to pel1 pel2 pel3 pel4 for the accumulation of anthocyanins in most of our experiments. The pel alleles attenuated growth, altered the accumulation of nutrients in seeds, disrupted an abscisic acid-inducible inhibition of seedling growth response that promotes drought tolerance, and affected the expression of genes associated with diverse biological functions, such as stress responses, cell wall metabolism hormone responses, signaling, growth, and the accumulation of phenylpropanoids and pigments. We found that PEL proteins specifically bind 6 transcription factors that influence the accumulation of anthocyanins, GLK2, and the carboxy termini of GLK1 and Arabidopsis thaliana myeloblastosis oncogene homolog 4 (AtMYB4). Our data indicate that the PEL proteins influence the accumulation of chlorophyll and many other processes, possibly by inhibiting GLK transcription factors and via other mechanisms, and that multiple mechanisms downregulate chlorophyll content.

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

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

Genome-wide identification, classification and expression analysis of MYB gene family in coconut (Cocos nucifera L.)

MYB transcription factors regulate the growth, development, and secondary metabolism of plant species. To investigate the origin of color variations in coconut pericarp, we identified and analyzed the MYB gene family present in coconut. According to the sequence of MYB genes in Arabidopsis thaliana, homologous MYB gene sequences were found in the whole genome database of coconut, the conserved sequence motifs within MYB proteins were analyzed by Motif Elicitation (MEME) tool, and the sequences without conservative structure were eliminated. Additionally, we employed RNA-seq technology to generate gene expression signatures of the R2R3-MYB genes across distinctive coconut parts exhibiting diverse colors. To validate these profiles, we conducted quantitative PCR (qPCR). Through comprehensive genome-wide screening, we successfully identified a collection of 179 MYB genes in coconut. Subsequent phylogenetic analysis categorized these 179 coconut MYB genes into 4-subfamilies: 124 R2R3-MYB, 4 3R-MYB types, 4 4R-MYB type, and 47 unknown types. Furthermore, these genes were further divided into 34 subgroups, with 28 of these subgroups successfully classified into known subfamilies found in Arabidopsis thaliana. By mapping the CnMYB genes onto the 16 chromosomes of the coconut genome, we unveiled a collinearity association between them. Moreover, a preservation of gene structure and motif distribution was observed across the CnMYB genes. Our research encompassed a thorough investigation of the R2R3-MYB genes present in the coconut genome, including the chromosomal localization, gene assembly, conserved regions, phylogenetic associations, and promoter cis-acting elements of the studied genes. Our findings revealed a collection of 12 R2R3-MYB candidate genes, namely CnMYB8, CnMYB15, CnMYB27, CnMYB28, CnMYB61, CnMYB63, CnMYB68, CnMYB94, CnMYB101, CnMYB150, CnMYB153, and CnMYB164. These genes showed differential expressions in diverse tissues and developmental stages of four coconut species, such as CnMYB68, CnMYB101, and CnMYB28 exhibited high expression in majority of tissues and coconut species, while CnMYB94 and CnMYB164 showed lower expression. These findings shed light on the crucial functional divergence of CnMYB genes across various coconut tissues, suggesting these genes as promising candidate genes for facilitating color development in this important crop.

Involvement of a MYB Transcription Factor in Anthocyanin Biosynthesis during Chinese Bayberry (Morella rubra) Fruit Ripening

Anthocyanin is a class of water-soluble flavonoids found in Chinese bayberry (Morella rubra) that is not only responsible for the variety of colors visible in nature but also has numerous health-promoting benefits in humans. Through comparative transcriptomics, we isolated and identified a transcription factor (TF) of the R2R3-MYB type, MrMYB9, in order to explore the anthocyanin biosynthesis pathway in red and white Chinese bayberries. MrMYB9 transcript was positively correlated with anthocyanin level and anthocyanin biosynthetic gene expression during Chinese bayberry fruit maturation (R-values in the range 0.54-0.84, p < 0.05). Sequence analysis revealed that MrMYB9 shared a similar R2R3 domain with MYB activators of anthocyanin biosynthesis in other plants. MrMYB9 substantially transactivated promoters of anthocyanin biosynthesis-related EBGs (MrCHI, MrF3'H, and MrANS) and LBGs (MrUFGT) upon co-expression of the AtEGL3 gene. Our findings indicated that MrMYB9 may positively modulate anthocyanin accumulation in Chinese bayberry.

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

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

Transcriptome-wide expression analysis of MYB gene family leads to functional characterization of flavonoid biosynthesis in fruit coloration of Ziziphus Mill

The Ziziphus mauritiana Lam. and Z. jujuba Mill. are the two most economically important members of the genus Ziziphus. The fruit color of Z. mauritiana remains green throughout fruit development in the majority of commercial cultivars, whereas its close relative, Z. jujuba Mill. turns from green to red in all cultivars. However, the lack of transcriptomic and genomic information confines our understanding of the molecular mechanisms underlying fruit coloration in Z. mauritiana (Ber). In the present study, we performed the transcriptome-wide analysis of MYB transcription factors (TFs) genes in Z. mauritiana and Z. jujuba, and identified 56 ZmMYB and 60 ZjMYB TFs in Z. mauritiana and Z. jujuba, respectively. Through transcriptomic expression analysis, four similar MYB genes (ZmMYB/ZjMYB13, ZmMYB/ZjMYB44, ZmMYB/ZjMYB50, and ZmMYB/ZjMYB56) from Z. mauritiana and Z. jujuba were selected as candidate key genes regulating flavonoid biosynthesis. Among these genes, the ZjMYB44 gene was transiently highly expressed in fruit, and flavonoid content accumulation also increased, indicating that this gene can influence flavonoid content during the period of fruit coloration in Z. jujuba. The current study adds to our understanding of the classification of genes, motif structure, and predicted functions of the MYB TFs, as well as identifying MYBs that regulate flavonoid biosynthesis in Ziziphus (Z. mauritiana and Z. jujuba). Based on this information, we concluded that MYB44 is involved in the flavonoids biosynthesis pathway during the fruit coloring of Ziziphus. Our research results provide an important understanding of the molecular mechanism of flavonoid biosynthesis resulting in fruit coloration and laying a foundation for further genetic improvement of fruit color in Ziziphus.

Identification and functional characterization of AcMYB113 in anthocyanin metabolism of Aesculus chinensis Bunge var. chinensis leaves

Anthocyanins can be induced by environmental factors such as low-temperature and play essential roles in plant color formation. In this study, leaves of Aesculus chinensis Bunge var. chinensis with different colors under natural low-temperature in autumn were collected and grouped into green leaf (GL) and red leaf (RL). To reveal the underlying mechanism of color formation in RL, a combined analysis of the metabolome and transcriptome was conducted with GL and RL. Metabolic analyses revealed that total anthocyanin content and primary anthocyanin components were increased RL relative to GL and cyanidin was the main anthocyanin compound in RL. Transcriptome analysis provided a total of 18720 differentially expressed genes (DEGs), of which 9150 DEGs were upregulated and 9570 DEGs were downregulated in RL relative to GL. KEGG analysis showed that DEGs were mainly enriched in flavonoid biosynthesis, phenylalanine metabolism, and phenylpropanoid biosynthesis. Furthermore, co-expression network analysis indicated that 56 AcMYB transcription factors were highly expressed in RL compared with GL, among which AcMYB113 (an R2R3-MYB TF) had a strong correlation with anthocyanins. Overexpression of AcMYB113 in apple resulted in dark-purple transgenic calluses. In addition, the transient expression experiment showed that AcMYB113 enhanced anthocyanin synthesis by activating pathways of anthocyanin biosynthesis in leaves of Aesculus chinensis Bunge var. chinensis. Taken together, our findings reveal new insights into the molecular mechanism of anthocyanin accumulation in RL and provide candidate genes for the breeding of anthocyanin-rich cultivars.

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 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.

Genome-wide identification and expression analysis of MYB gene family under nitrogen stress in Panax notoginseng

The myeloblastosis (MYB) gene family, involved in regulating many important physiological and biochemical processes, is one of the largest transcript factor superfamilies in plants. Since the identification of genome sequencing of Panax notoginseng has been completed, there was little known about the whole genome of its specific MYB gene family and the response to abiotic stresses, in consideration of the excessive application of nitrogen fertilizers in P. notoginseng. In this study, 123 PnMYB genes (MYB genes of P. notoginseng) have been identified and divided into 3 subfamilies by the phylogenetic analysis. These PnMYB genes were unevenly located on 12 chromosomes. Meanwhile, the gene structure and protein conserved domain were established by MEME Suite. The analysis of collinear relationships reflected that there were 121 homologous genes between P. notoginseng and Arabidopsis and 30 between P. notoginseng and rice. Moreover, cis-acting elements of PnMYB gene promoters were predicted which indicated that PnMYBs are involved in biotic, abiotic stress, and hormone induction. The expressions of PnMYB transcription factors in its roots, flowers, and leaves were detected by qRT-PCR and they had tissue-specific expressions and related to the growth of different tissues. Under nitrogen stress, MYB transcription factors had great feedback. Ten R2R3-MYB subfamily genes were significantly induced and indicated the possible function of protecting P. notoginseng from excess nitrogen. With further knowledge on identification of PnMYB gene related to tissue selectivity and abiotic stresses, this study laid the foundation for the functional development of PnMYB gene family and improved the cultivation of P. notoginseng.

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.

Anthocyanin Biosynthesis Associated with Natural Variation in Autumn Leaf Coloration in Quercus aliena Accessions

Quercus aliena is an economically important tree species and one of the dominant native oak species in China. Although its leaves typically turn yellow in autumn, we observed natural variants with red leaves. It is important to understand the mechanisms involved in leaf color variation in this species. Therefore, we compared a Q. aliena tree with yellow leaves and three variants with red leaves at different stages of senescence in order to determine the causes of natural variation. We found that the accumulation of anthocyanins such as cyanidin 3-O-glucoside and cyanidin 3-O-sambubiglycoside had a significant effect on leaf coloration. Gene expression analysis showed upregulation of almost all genes encoding enzymes involved in anthocyanin synthesis in the red-leaved variants during the early and main discoloration stages of senescence. These findings are consistent with the accumulation of anthocyanin in red variants. Furthermore, the variants showed significantly higher expression of transcription factors associated with anthocyanin synthesis, such as those encoded by genes QaMYB1 and QaMYB3. Our findings provide new insights into the physiological and molecular mechanisms involved in autumn leaf coloration in Q. aliena, as well as provide genetic resources for further development and cultivation of valuable ornamental variants of this species.

Systematic Analysis and Functional Characterization of R2R3-MYB Genes in Scutellaria baicalensis Georgi

R2R3-MYB transcription factors participate in multiple critical biological processes, particularly as relates to the regulation of secondary metabolites. The dried root of Scutellaria baicalensis Georgi is a traditional Chinese medicine and possesses various bioactive attributes including anti-inflammation, anti-HIV, and anti-COVID-19 properties due to its flavonoids. In the current study, a total of 95 R2R3-MYB genes were identified in S. baicalensis and classified into 34 subgroups, as supported by similar exon–intron structures and conserved motifs. Among them, 93 R2R3-SbMYBs were mapped onto nine chromosomes. Collinear analysis revealed that segmental duplications were primarily responsible for driving the evolution and expansion of the R2R3-SbMYB gene family. Synteny analyses showed that the ortholog numbers of the R2R3-MYB genes between S. baicalensis and other dicotyledons had a higher proportion compared to that which is found from the monocotyledons. RNA-seq data indicated that the expression patterns of R2R3-SbMYBs in different tissues were different. Quantitative reverse transcriptase-PCR (qRT-PCR) analysis showed that 36 R2R3-SbMYBs from different subgroups exhibited specific expression profiles under various conditions, including hormone stimuli treatments (methyl jasmonate and abscisic acid) and abiotic stresses (drought and cold shock treatments). Further investigation revealed that SbMYB18/32/46/60/70/74 localized in the nucleus, and SbMYB18/32/60/70 possessed transcriptional activation activity, implying their potential roles in the regulatory mechanisms of various biological processes. This study provides a comprehensive understanding of the R2R3-SbMYBs gene family and lays the foundation for further investigation of their biological function.

The UV-B-Induced Transcription Factor HY5 Regulated Anthocyanin Biosynthesis in Zanthoxylum bungeanum

Pericarp color is an important economic characteristic of Zanthoxylum bungeanum. Anthocyanins are the main reason for the pericarp's red appearance in Z. bungeanum. In this study, through the combined analysis of the metabolome and transcriptome, HY5, whose expression is highly correlated to changes in the anthocyanin content, was screened and identified. Under natural ripening conditions, the Z. bungeanum fruit gradually changed in color from green to red, while bagging resulted in the fruit maintaining its green color. After unbagging, the fruit gradually turned red, and the ZbHY5 expression and anthocyanin content increased. In addition, the leaves changed from green to red after exposure to UV-B radiation, and the ZbHY5 expression and anthocyanin content increased. The transient overexpression of ZbHY5 deepened the redness of the Z. bungeanum leaves and promoted the expression of ZbHY5 and ZbMYB113 as well as anthocyanin accumulation. Bimolecular fluorescence complementation (BIFC) showed that there was an interaction between ZbHY5 and ZbMYB113. These results revealed that under UV-B irradiation, ZbHY5 might regulate the expression levels of the structural genes related to anthocyanin biosynthesis through combination with ZbMYB113, thereby affecting anthocyanin accumulation. This finding provides useful insights for further studies focusing on UV-B-induced anthocyanin accumulation in Z. bungeanum.

Genome-Wide Comparative Analysis of the R2R3-MYB Gene Family in Five Solanaceae Species and Identification of Members Regulating Carotenoid Biosynthesis in Wolfberry

The R2R3-MYB is a large gene family involved in various plant functions, including carotenoid biosynthesis. However, this gene family lacks a comprehensive analysis in wolfberry (Lycium barbarum L.) and other Solanaceae species. The recent sequencing of the wolfberry genome provides an opportunity for investigating the organization and evolutionary characteristics of R2R3-MYB genes in wolfberry and other Solanaceae species. A total of 610 R2R3-MYB genes were identified in five Solanaceae species, including 137 in wolfberry. The LbaR2R3-MYB genes were grouped into 31 subgroups based on phylogenetic analysis, conserved gene structures, and motif composition. Five groups only of Solanaceae R2R3-MYB genes were functionally divergent during evolution. Dispersed and whole duplication events are critical for expanding the R2R3-MYB gene family. There were 287 orthologous gene pairs between wolfberry and the other four selected Solanaceae species. RNA-seq analysis identified the expression level of LbaR2R3-MYB differential gene expression (DEGs) and carotenoid biosynthesis genes (CBGs) in fruit development stages. The highly expressed LbaR2R3-MYB genes are co-expressed with CBGs during fruit development. A quantitative Real-Time (qRT)-PCR verified seven selected candidate genes. Thus, Lba11g0183 and Lba02g01219 are candidate genes regulating carotenoid biosynthesis in wolfberry. This study elucidates the evolution and function of R2R3-MYB genes in wolfberry and the four Solanaceae species.

Identification of the Regulatory Genes of UV-B-Induced Anthocyanin Biosynthesis in Pepper Fruit

Fruit peels of certain pepper (Capsicum annum L.) varieties accumulate a large amount of anthocyanins and exhibit purple color under medium-wave ultraviolet (UV-B) conditions, which severely impacts the commodity value of peppers. However, the regulatory mechanism of the above process has not been well studied so far. To explore which key genes are involved in this regulatory mechanism, pepper variety 19Q6100, the fruit peels of which turn purple under UV-B conditions, was investigated in this study. Transcription factors with expression levels significantly impacted by UV-B were identified by RNA-seq. Those genes may be involved in the regulation of UV-B-induced anthocyanin biosynthesis. Yeast one-hybrid results revealed that seven transcription factors, CabHLH143, CaMYB113, CabHLH137, CaMYBG, CaWRKY41, CaWRKY44 and CaWRKY53 directly bound to the putative promotor regions of the structural genes in the anthocyanin biosynthesis pathway. CaMYB113 was found to interact with CabHLH143 and CaHY5 by yeast two-hybrid assay, and those three genes may participate collaboratively in UV-B-induced anthocyanin biosynthesis in pepper fruit. Virus-induced gene silencing (VIGS) indicated that fruit peels of CaMYB113-silenced plants were unable to turn purple under UV-B conditions. These findings could deepen our understanding of UV-B-induced anthocyanin biosynthesis in pepper.