Genome-wide identification and expression analysis of the MYB transcription factor in Japanese plum (Prunus salicina)

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

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

Comparative genomic studies on the TGF-β superfamily in blue whale

TGF-β supergene family has a wide range of physiological functions including cell adhesion, motility, proliferation, apoptosis, and differentiation. We systematically analyzed and characterized the TGF-β gene superfamily from the whole blue whale (Balaenoptera musculus) genome, using comparative genomic and evolutionary analysis. We identified 30 TGF-β genes and were split into two subgroups, BMP-like and TGF-like. All TGF-β proteins demonstrating a basic nature, with the exception of BMP1, BMP2, BMP10, GDF2, MSTN, and NODAL modulator, had acidic characteristics. All the blue whale (B. musculus) TGF-β proteins, excluding BMP1, are thermostable based on aliphatic index. The instability index showed all proteins except the NODAL modulator was unstable. TGF-β proteins showed a hydrophilic character, with the exception of GDF1 and INHBC. Moreover, all the detected TGF-β genes showed evolutionary conserved nature. A segmental duplication was indicated by TGF-β gene family, and the Ka/Ks ratio showed that the duplicated gene pairs were subjected to selection pressure, indicating both purifying and positive selection pressure. Two possible recombination breakpoints were also predicted. This study provides insights into the genetic characterization and evolutionary aspects of the TGF-β superfamily in blue whales (B. musculus).

Comparative genomic analysis of the Growth-Regulating Factors-Interacting Factors (GIFs) in six Salicaceae species and functional analysis of PeGIF3 reveals their regulatory role in Populus heteromorphic leaves

BACKGROUND

The growth-regulating factor-interacting factor (GIF) gene family plays a vital role in regulating plant growth and development, particularly in controlling leaf, seed, and root meristem homeostasis. However, the regulatory mechanism of heteromorphic leaves by GIF genes in Populus euphratica as an important adaptative trait of heteromorphic leaves in response to desert environment remains unknown.

RESULTS

This study aimed to identify and characterize the GIF genes in P. euphratica and other five Salicaceae species to investigate their role in regulating heteromorphic leaf development. A total of 27 GIF genes were identified and characterized across six Salicaceae species (P. euphratica, Populus pruinose, Populus deltoides, Populus trichocarpa, Salix sinopurpurea, and Salix suchowensis) at the genome-wide level. Comparative genomic analysis among these species suggested that the expansion of GIFs may be derived from the specific Salicaceae whole-genome duplication event after their divergence from Arabidopsis thaliana. Furthermore, the expression data of PeGIFs in heteromorphic leaves, combined with functional information on GIF genes in Arabidopsis, indicated the role of PeGIFs in regulating the leaf development of P. euphratica, especially PeGIFs containing several cis-acting elements associated with plant growth and development. By heterologous expression of the PeGIF3 gene in wild-type plants (Col-0) and atgif1 mutant of A. thaliana, a significant difference in leaf expansion along the medial-lateral axis, and an increased number of leaf cells, were observed between the overexpressed plants and the wild type.

CONCLUSION

PeGIF3 enhances leaf cell proliferation, thereby resulting in the expansion of the central-lateral region of the leaf. The findings not only provide global insights into the evolutionary features of Salicaceae GIFs but also reveal the regulatory mechanism of PeGIF3 in heteromorphic leaves of P. euphratica.

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.

Comparative analysis of the MYB gene family in seven Ipomoea species

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

Genome-wide analysis of R2R3-MYB transcription factors reveals their differential responses to drought stress and ABA treatment in desert poplar (Populus euphratica)

The R2R3-MYB transcription factors are widely involved in the regulation of plant growth, biotic and abiotic stress responses. Meanwhile, seed germination, which is stimulated by internal and external environments, is a critical stage in the plant life cycle. However, the identification, characterization, and expression profiling of the Populus euphratica R2R3-MYB family in drought response during seed germination have been unknown. Our study attempted to identify and characterize the R2R3-MYB genes in P. euphratica (PeR2R3-MYBs) and explore how R2R3-MYBs trigger the drought and abscisic acid (ABA) response mechanism in its seedlings. Based on the analysis of comparative genomics, 174 PeR2R3-MYBs were identified and expanded driven by whole genome duplication or segment duplication events. The analysis of Ka/Ks ratios showed that, in contrast to most PeR2R3-MYBs, the other PeR2R3-MYBs were subjected to positive selection in P. euphratica. Further, the expression data of PeR2R3-MYBs under drought stress and ABA treatment, together with available functional data for Arabidopsis thaliana MYB genes, supported the hypothesis that PeR2R3-MYBs involved in response to drought are dependent or independent on ABA signaling pathway during seed germination, especially PeR2R3-MYBs with MYB binding sites (MBS) cis-element and/or tandem duplication. This study is the first report on the genome-wide analysis of PeR2R3-MYBs, as well as the other two Salicaceae species. The duplication events and differential expressions of PeR2R3-MYBs play important roles in enhancing the adaptation to drought desert environment. Our results provide a reference for prospective functional studies of R2R3-MYBs of poplars and lay the foundation for new breeding strategies to improve the drought tolerance of P. euphratica.

VrMYB90 Functions Synergistically with VrbHLHA and VrMYB3 to Regulate Anthocyanin Biosynthesis in Mung Bean

Mung bean is an important grain-legume crop and its sprout is an economical and nutrient vegetable for the public, but the genetic regulation of anthocyanin production, which is an antioxidant in mung bean, remains elusive. In our study, we characterized a subgroup (SG) 6 R2R3-MYB anthocyanin activator VrMYB90 and a SG 4 R2R3-MYB anthocyanin repressor VrMYB3, which synergistically function in regulating anthocyanin synthesis with VrbHLHA transcription factor. The overexpressed VrMYB90 protein activates the expression of VrMYB3 and VrbHLHA in mung bean hair roots, and also promotes VrDFR and VrANS transcript levels by directly binding to the corresponding promoters at specific motifs (CAACTG and CCGTTG). VrMYB90 interacts with VrbHLHA to enhance its regulatory activities on VrDFR and VrANS. Furthermore, the interaction between VrMYB3 with VrMYB90 and VrbHLHA could result in the restriction of anthocyanin synthesis to prevent excessive anthocyanin accumulation. Our results demonstrate that the VrMYB90 protein, in conjunction with VrMYB3 and VrbHLHA, forms a key regulatory module to fine-tune anthocyanin synthesis in mung bean.

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.

An Upgraded, Highly Saturated Linkage Map of Japanese Plum (Prunus salicina Lindl.), and Identification of a New Major Locus Controlling the Flavan-3-ol Composition in Fruits

Japanese plum fruits are rich in phenolic compounds, such as anthocyanins and flavan-3-ols, whose contents vary significantly among cultivars. Catechin (C) and epicatechin (EC) are flavan-3-ol monomers described in the fruits of this species and are associated with bitterness, astringency, antioxidant capacity, and susceptibility to enzymatic mesocarp browning. In this study, we aimed to identify quantitative trait loci (QTL) associated with the content of flavan-3-ol in Japanese plum fruits. We evaluated the content of C and EC in the mesocarp and exocarp of samples from 79 and 64 seedlings of an F1 progeny (<'98-99' × 'Angeleno'>) in the first and second seasons, respectively. We also constructed improved versions of linkage maps from '98-99' and 'Angeleno,' presently called single-nucleotide polymorphisms (SNPs) after mapping the already available GBS reads to Prunus salicina Lindl. cv. 'Sanyueli' v2.0 reference genome. These data allowed for describing a cluster of QTLs in the cultivar, 'Angeleno,' associated with the flavan-3-ol composition of mesocarp and exocarp, which explain up to 100% of the C/EC ratio. Additionally, we developed a C/EC metabolic marker, which was mapped between the markers with the highest log of odds (LOD) scores detected by the QTL analysis. The C/EC locus was located in the LG1, at an interval spanning 0.70 cM at 108.30-108.90 cM. Our results suggest the presence of a novel major gene controlling the preferential synthesis of C or EC in the Japanese plum fruits. This study is a significant advance in understanding the regulation of synthesizing compounds associated with fruit quality, postharvest, and human health promotion.

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.

Genome wide characterization of R2R3 MYB transcription factor from Apocynum venetum revealed potential stress tolerance and flavonoid biosynthesis genes

MYB transcription factors are crucial in regulating stress tolerance and expression of major genes involved in flavonoid biosynthesis. The functions of MYBs is well explored in a number of plants, yet no studies is reported in Apocynum venetum. We identified a total of 163 MYB candidates, that comprised of 101 (61.96%) R2R3, 6 3R, 1 4R and 55 1R. Syntenic analysis of A. venetum R2R3 (AvMYB) showed highest orthologous pairs with Vitis vinifera MYBs followed by Arabidopsis thaliana among the four species evaluated. Thirty segmental duplications and 6 tandem duplications were obtained among AvMYB gene pairs signifying their role in the MYB gene family expansion. Nucleotide substitution analysis (Ka/Ks) showed the AvMYBs to be under the influence of strong purifying selection. Expression analysis of selected AvMYB under low temperature and cadmium stresses resulted in the identification of AvMYB48, AvMYB97, AvMYB8,AvMYB4 as potential stress responsive genes and AvMYB10 and AvMYB11 in addition, proanthocyanidin biosynthesis regulatory genes which is consistent with their annotated homologues in Arabidopsis. Tissue specific expression profile analysis of AvMYBs further supported the qPCR analysis result. MYBs with higher transcript levels in root, stem and leaf like AvMYB4 forexample, was downregulated under the stresses and such with low transcript level such as AvMYB48 which had low transcript in the leaf was upregulated under both stresses. Transcriptome and phylogenetic analysis suggested AvMYB42 as a potential regulator of anthocyanin biosynthesis. Thus, this study provided valuable information on AvR2R3-MYB gene family with respect to stress tolerance and flavonoid biosynthesis.

Genomic Analysis of the Glutathione S-Transferase Family in Pear (Pyrus communis) and Functional Identification of PcGST57 in Anthocyanin Accumulation

Anthocyanin accumulation in vacuoles results in red coloration in pear peels. Glutathione S-transferase (GST) proteins have emerged as important regulators of anthocyanin accumulation. Here, a total of 57 PcGST genes were identified in the European pear ‘Bartlett’ (Pyrus communis) through comprehensive genomic analysis. Phylogenetic analysis showed that PcGST genes were divided into 10 subfamilies. The gene structure, chromosomal localization, collinearity relationship, cis-elements in the promoter region, and conserved motifs of PcGST genes were analyzed. Further research indicated that glutamic acid (Glu) can significantly improve anthocyanin accumulation in pear peels. RNA sequencing (RNA-seq) analysis showed that Glu induced the expression of most PcGST genes, among which PcGST57 was most significantly induced. Further phylogenetic analysis indicated that PcGST57 was closely related to GST genes identified in other species, which were involved in anthocyanin accumulation. Transcript analysis indicated that PcGST57 was expressed in various tissues, other than flesh, and associated with peel coloration at different developmental stages. Silencing of PcGST57 by virus-induced gene silencing (VIGS) inhibited the expression of PcGST57 and reduced the anthocyanin content in pear fruit. In contrast, overexpression of PcGST57 improved anthocyanin accumulation. Collectively, our results demonstrated that PcGST57 was involved in anthocyanin accumulation in pear and provided candidate genes for red pear breeding.

Evaluation of reference genes and characterization of the MYBs in xylem radial change of Chinese fir stem

The radial change (RC) of tree stem is the process of heartwood formation involved in complex molecular mechanism. Chinese fir (Cunninghamia lanceolata (Lamb.) Hook.), an evergreen species, is an important fast-growing timber tree in southern China. In this study, the top four stable genes (IDH, UBC2, RCA and H2B) were selected in RC tissues of 15 years old Chinese fir stem (RC15) and the genes (H2B, 18S, TIP41 and GAPDH) were selected in RC tissues of 30 years old Chinese fir stem (RC30). The stability of the reference genes is higher in RC30 than in RC15. Sixty-one MYB transcripts were obtained on the PacBio Sequel platform from woody tissues of one 30 years old Chinese fir stem. Based on the number of MYB DNA-binding domain and phylogenetic relationships, the ClMYB transcripts contained 21 transcripts of MYB-related proteins (1R-MYB), 39 transcripts of R2R3-MYB proteins (2R-MYB), one transcript of R1R2R3-MYB protein (3R-MYB) belonged to 18 function-annotated clades and two function-unknown clades. In RC woody tissues of 30 years old Chinese fir stem, ClMYB22 was the transcript with the greatest fold change detected by both RNA-seq and qRT-PCR. Reference genes selected in this study will be helpful for further verification of transcript abundance patterns during the heartwood formation of Chinese fir.

Functional Characterization of Hedychium coronarium J. Koenig MYB132 Confers the Potential Role in Floral Aroma Synthesis

The R2R3-MYB transcription factors (TFs) play several key roles in numerous plant biological processes. Hedychium coronarium is an important ornamental plant well-known for its elegant flower shape and abundant aroma type. The floral aroma of H. coronarium is due to the presence of a large amount of terpenes and benzenoids. However, less is known about the role of R2R3-MYB TFs in the regulatory mechanism of floral aroma production in this breed. Herein, we isolate and functionally characterize the R2R3-MYB TF HcMYB132, which is potentially involved in regulating floral aroma synthesis. Sequence alignment analysis revealed that it includes a nuclear localization signal NLS(s) and a 2R, 3R motif signature in the sequences. A subcellular localization assay revealed that HcMYB132 protein localizes to the nucleus. Real-time qPCR assays showed that HcMYB132 is specifically expressed in flowers and its expression pattern correlates with the emission of floral volatile compounds. In HcMYB132-silenced flowers, the levels of floral volatile compounds were significantly reduced, and the expression of key structural volatile synthesis genes was downregulated compared to control. Collectively, these results suggest that HcMYB132 might play a significant role in the regulation of terpenoid biosynthesis in H. coronarium.

Genome-Wide Identification and Expression Analysis of MYB Transcription Factor Superfamily in Dendrobium catenatum

Dendrobium catenatum is an important traditional Chinese medicine and naturally grows on tree trunks and cliffs, where it can encounter diverse environmental stimuli. MYB transcription factors are widely involved in response to abiotic stresses. However, the MYB gene family has not yet been systematically cataloged in D. catenatum. In this study, a total of 133 MYB proteins were identified in D. catenatum, including 32 MYB-related, 99 R2R3-MYB, 1 3R-MYB, and 1 4R-MYB proteins. Phylogenetic relationships, conserved motifs, gene structures, and expression profiles in response to abiotic stresses were then analyzed. Phylogenetic analysis revealed MYB proteins in D. catenatum could be divided into 14 subgroups, which was supported by the conserved motif compositions and gene structures. Differential DcMYB gene expression and specific responses were analyzed under drought, heat, cold, and salt stresses using RNA-seq and validated by qRT-PCR. Forty-two MYB genes were differentially screened following exposure to abiotic stresses. Five, 12, 11, and 14 genes were specifically expressed in response to drought, heat, cold, and salt stress, respectively. This study identified candidate MYB genes with possible roles in abiotic tolerance and established a theoretical foundation for molecular breeding of D. catenatum.

Genome-wide analysis of MYB transcription factors and their responses to salt stress in Casuarina equisetifolia

BACKGROUND

MYB transcription factors are a kind of DNA binding protein that can specifically interact with the promoter region. Members of MYB TFs are widely involved in plant growth and development, secondary metabolism, stress response, and hormone signal transduction. However, there is no report of comprehensive bioinformatics analysis on the MYB family of Casuarina equisetifolia.

RESULTS

In this study, bioinformatics methods were used to screen out 182 MYB transcription factors from the Casuarina equisetifolia genome database, including 69 1R-MYB, 107 R2R3-MYB, 4 R1R2R3-MYB, and 2 4R-MYB. The C. equisetifolia R2R3-MYB genes were divided into 29 groups based on the phylogenetic topology and the classification of the MYB superfamily in Arabidopsis thaliana, while the remaining MYB genes (1R-MYB, R1R2R3-MYB, and 4R-MYB) was divided into 19 groups. Moreover, the conserved motif and gene structure analysis shown that the members of the CeqMYBs were divided into the same subgroups with mostly similar gene structures. In addition, many conserved amino acids in the R2 and R3 domains of CeqMYBs by WebLogo analysis, especially tryptophan residues (W), with 3 conserved W in R2 repeat and 2 conserved W in R3 repeat. Combining promoter and GO annotation analysis, speculated on the various biological functions of CeqMYBs, thus 32 MYB genes were selected to further explore its response to salt stress by using qPCR analysis technique. Most CeqMYB genes were differentially regulated following multiple salt treatments.

CONCLUSIONS

Seven genes (CeqMYB164, CeqMYB4, CeqMYB53, CeqMYB32, CeqMYB114, CeqMYB71 and CeqMYB177) were assigned to the "response to salt stress" by GO annotation. Among them, the expression level of CeqMYB4 was up-regulated under various salt treatments, indicating CeqMYB4 might participated in the response to salt stress. Our results provide important information for the biological function of C. equisetifolia, as well as offer candidate genes for further study of salt stress mechanism.

Genome-Wide Analysis Reveals the Potential Role of MYB Transcription Factors in Floral Scent Formation in Hedychium coronarium

The MYB gene family is one of the largest groups of transcription factors (TFs) playing diverse roles in several biological processes. Hedychium coronarium (white ginger lily) is a renowned ornamental plant both in tropical and subtropical regions due to its flower shape and strong floral scent mainly composed of terpenes and benzenoids. However, there is no information available regarding the role of the MYB gene family in H. coronarium. In the current study, the MYB gene family was identified and extensively analyzed. The identified 253 HcMYB genes were unevenly mapped on 17 chromosomes at a different density. Promoter sequence analysis showed numerous phytohormones related to cis-regulatory elements. The majority of HcMYB genes contain two to three introns and motif composition analysis showed their functional conservation. Phylogenetic analysis revealed that HcMYBs could be classified into 15 distinct clades, and the segmental duplication events played an essential role in the expansion of the HcMYB gene family. Tissue-specific expression patterns of HcMYB genes displayed spatial and temporal expression. Furthermore, seven HcMYB (HcMYB7/8/75/79/145/238/248) were selected for further investigation. Through RT-qPCR, the response of candidates HcMYB genes toward jasmonic acid methyl ester (MeJA), abscisic acid (ABA), ethylene, and auxin was examined. Yeast one-hybrid (Y1H) assays revealed that candidate genes directly bind to the promoter of bottom structural volatile synthesis genes (HcTPS1, HcTPS3, HcTPS10, and HcBSMT2). Moreover, yeast two-hybrid (Y2H) assay showed that HcMYB7/8/75/145/248 interact with HcJAZ1 protein. In HcMYB7/8/79/145/248-silenced flowers, the floral volatile contents were decreased and downregulated the expression of key structural genes, suggesting that these genes might play crucial roles in floral scent formation in H. coronarium by regulating the expression of floral scent biosynthesis genes. Collectively, these findings indicate that HcMYB genes might be involved in the regulatory mechanism of terpenoids and benzenoid biosynthesis in H. coronarium.

Genome-Wide Analysis of MYB Gene Family in Chinese Bayberry (Morella rubra) and Identification of Members Regulating Flavonoid Biosynthesis

Chinese bayberry (Morella rubra), the most economically important fruit tree in the Myricaceae family, is a rich source of natural flavonoids. Recently the Chinese bayberry genome has been sequenced, and this provides an opportunity to investigate the organization and evolutionary characteristics of MrMYB genes from a whole genome view. In the present study, we performed the genome-wide analysis of MYB genes in Chinese bayberry and identified 174 MrMYB transcription factors (TFs), including 122 R2R3-MYBs, 43 1R-MYBs, two 3R-MYBs, one 4R-MYB, and six atypical MYBs. Collinearity analysis indicated that both syntenic and tandem duplications contributed to expansion of the MrMYB gene family. Analysis of transcript levels revealed the distinct expression patterns of different MrMYB genes, and those which may play important roles in leaf and flower development. Through phylogenetic analysis and correlation analyses, nine MrMYB TFs were selected as candidates regulating flavonoid biosynthesis. By using dual-luciferase assays, MrMYB12 was shown to trans-activate the MrFLS1 promoter, and MrMYB39 and MrMYB58a trans-activated the MrLAR1 promoter. In addition, overexpression of 35S:MrMYB12 caused a significant increase in flavonol contents and induced the expression of NtCHS, NtF3H, and NtFLS in transgenic tobacco leaves and flowers and significantly reduced anthocyanin accumulation, resulting in pale-pink or pure white flowers. This indicates that MrMYB12 redirected the flux away from anthocyanin biosynthesis resulting in higher flavonol content. The present study provides valuable information for understanding the classification, gene and motif structure, evolution and predicted functions of the MrMYB gene family and identifies MYBs regulating different aspects of flavonoid biosynthesis in Chinese bayberry.

Characterization of the Liriodendron Chinense MYB Gene Family and Its Role in Abiotic Stress Response

Liriodendron chinense (Lchi) is a Magnoliaceae plant, which is a basic angiosperm left behind by the Pleistocene and mainly distributed in the south of the Yangtze River. Liriodendron hybrids has good wood properties and is widely used in furniture and in other fields. It is not clear if they can adapt to different environmental conditions, such as drought and high and low temperatures, and the molecular mechanisms for this adaptation are unknown. Among plant transcription factors (TFs), the MYB gene family is one of the largest and is often involved in stress or adversity response signaling, growth, and development. Therefore, studying the role of MYBTFs in regulating abiotic stress signaling, growth, and development in Lchi is helpful to promote afforestation in different environments. In our research, a genome-wide analysis of the LchiMYB gene family was performed, including the phylogenetic relationship tree, gene exon-intron structure, collinearity, and chromosomal position. According to the evolutionary tree, 190 LchiMYBs were divided into three main branches. LchiMYBs were evenly distributed across 19 chromosomes, with their collinearity, suggesting that segment duplication events may have contributed to LchiMYB gene expansion. Transcriptomes from eight tissues, 11 stages of somatic embryogenesis, and leaves after cold, heat, and drought stress were used to analyze the function of the MYB gene family. The results of tissue expression analysis showed that most LchiMYB genes regulated bark, leaf, bud, sepal, stigma, and stamen development, as well as the four important stages (ES3, ES4, ES9, and PL) of somatic embryogenesis. More than 60 LchiMYBs responded to heat, cold, and drought stress; some of which underwent gene duplication during evolution. LchiMYB3 was highly expressed under all three forms of stress, while LchiMYB121 was strongly induced by both cold and heat stress. Eight genes with different expression patterns were selected and verified by quantitative real-time PCR (qRT-PCR) experiments. The results suggested that these LchiMYBs may regulate Lchi growth development and resistance to abiotic stress. This study shows the cross-regulatory function of LchiMYBs in the growth and development, asexual reproduction, and abiotic resistance of Lchi. This information will prove pivotal to directing further studies on the biological function of Lchi MYBTFs in genetic improvement and abiotic stress response.