The Evolutionary History of R2R3-MYB Proteins Across 50 Eukaryotes: New Insights Into Subfamily Classification and Expansion

Genome-wide identification of R2R3-MYB transcription factor subfamily genes involved in salt stress in rice (Oryza sativa L.)

BACKGROUND

R2R3-MYB transcription factors belong to one of the largest gene subfamilies in plants, and they are involved in diverse biological processes. However, the role of R2R3-MYB transcription factor subfamily genes in the response of rice (Oryza sativa L.) to salt stress has been rarely reported.

RESULTS

In this study, we performed a genome-wide characterization and expression identification of rice R2R3-MYB transcription factor subfamily genes. We identified a total of 117 R2R3-MYB genes in rice and characterized their gene structure, chromosomal location, and cis-regulatory elements. According to the phylogenetic relationships and amino acid sequence homologies, the R2R3-MYB genes were divided into four groups. qRT-PCR of the R2R3-MYB genes showed that the expression levels of 10 genes significantly increased after 3 days of 0.8% NaCl treatment. We selected a high expression gene OsMYB2-115 for further analysis. OsMYB2-115 was highly expressed in the roots, stem, leaf, and leaf sheath. OsMYB2-115 was found to be localized in the nucleus, and the yeast hybrid assay showed that OsMYB2-115 has transcriptional activation activity.

CONCLUSION

This result provides important information for the functional analyses of rice R2R3-MYB transcription factor subfamily genes related to the salt stress response and reveals that OsMYB2-115 may be an important gene associated with salt tolerance in rice.

R2R3-MYB Gene Family in Coptis teeta Wall.: Genome-Wide Identification, Phylogeny, Evolutionary Expansion, and Expression Analyses during Floral Development

The R2R3-MYB gene family represents a widely distributed class of plant transcription factors. This gene family plays an important role in many aspects of plant growth and development. However, the characterization of R2R3-MYB genes present in the genome of Coptis teeta has not been reported. Here, we describe the bioinformatic identification and characterization of 88 R2R3-MYB genes in this species, and the identification of members of the R2R3-MYB gene family in species within the order Ranales most closely related to Coptis teeta. The CteR2R3-MYB genes were shown to exhibit a higher degree of conservation compared to those of A. thaliana, as evidenced by phylogeny, conserved motifs, gene structure, and replication event analyses. Cis-acting element analysis confirmed the involvement of CteR2R3-MYB genes in a variety of developmental processes, including growth, cell differentiation, and reproduction mediated by hormone synthesis. In addition, through homology comparisons with the equivalent gene family in A. thaliana, protein regulatory network prediction and transcriptome data analysis of floral organs across three time periods of flower development, 17 candidate genes were shown to exhibit biased expression in two floral phenotypes of C. teeta. This suggests their potential involvement in floral development (anther development) in this species.

MaMYBR30, a Novel 1R-MYB, Plays Important Roles in Plant Development and Abiotic Stress Resistance

The V-myb myeloblastosis viral oncogene homolog (MYB) family participate in various bioprocesses including development and abiotic stress responses. In the present study, we first report a 1R SHAQKYF-class MYB, MaMYBR30, in mulberry. Subcellular localization and sequence analysis indicated MaMYBR30 is located in the nucleus and belongs to a CCA-like subgroup with a conserved SHAQKYF motif. Expression profile analysis showed that MaMYBR30 is expressed in leaves and can be induced by drought and salt stress. The down-regulation of MaMYBR30 using virus-induced gene silence (VIGS) in mulberry and the overexpression of MaMYBR30 in Arabidopsis were induced to explore the function of MaMYBR30. The functional characterization of MaMYBR30 in vivo indicated that MaMYBR30 can positively regulate the resistance of mulberry to drought while negatively regulating the resistance of mulberry to salt stress. In addition, MaMYBR30 also affects flower development and reproductive growth, especially after exposure to salt stress. Weighted gene co-expression network analysis (WGCNA) primarily revealed the possible genes and signal pathways that are regulated by MaMYBR30. Our results also imply that complex molecular mechanisms mediated by MaMYBR30, including crosstalk of ion toxicity, phytohormone signal transduction, flowering development, and epigenetic modification, need to be further explored in the future.

SmuMYB113 is the determinant of fruit color in pepino (Solanum muricatum)

Pepino (Solanum muricatum) is an herbaceous crop phylogenetically related to tomato and potato. Pepino fruit vary in color, size and shape, and are eaten fresh. In this study, we use pepino as a fruit model to understand the transcriptional regulatory mechanisms controlling fruit quality. To identify the key genes involved in anthocyanin biosynthesis in pepino, two genotypes were studied that contrasted in foliar and fruit pigmentation. Anthocyanin profiles were analyzed, as well as the expression of genes that encode enzymes for anthocyanin biosynthesis and transcriptional regulators using both RNA-seq and quantitative PCR. The differential expression of the transcription factor genes R2R3 MYB SmuMYB113 and R3MYB SmuATV suggested their association with purple skin and foliage phenotype. Functional analysis of these genes in both tobacco and pepino showed that SmuMYB113 activates anthocyanins, while SmuATV suppresses anthocyanin accumulation. However, despite elevated expression in all tissues, SmuMYB113 does not significantly elevate flesh pigmentation, suggesting a strong repressive background in fruit flesh tissue. These results will aid understanding of the differential regulation controlling fruit quality aspects between skin and flesh in other fruiting species.

Identification of a candidate gene for the I locus determining the dominant white bulb color in onion (Allium cepa L.)

KEY MESSAGE

Through a map-based cloning approach, a gene coding for an R2R3-MYB transcription factor was identified as a causal gene for the I locus controlling the dominant white bulb color in onion. White bulb colors in onion (Allium cepa L.) are determined by either the C or I loci. The causal gene for the C locus was previously isolated, but the gene responsible for the I locus has not been identified yet. To identify candidate genes for the I locus, an approximately 7-Mb genomic DNA region harboring the I locus was obtained from onion and bunching onion (A. fistulosum) whole genome sequences using two tightly linked molecular markers. Within this interval, the AcMYB1 gene, known as a positive regulator of anthocyanin production, was identified. No polymorphic sequences were found between white and red AcMYB1 alleles in the 4,860-bp full-length genomic DNA sequences. However, a 4,838-bp LTR-retrotransposon was identified in the white allele, in the 79-bp upstream coding region from the stop codon. The insertion of this LTR-retrotransposon created a premature stop codon, resulting in the replacement of 26 amino acids with seven different residues. A molecular marker was developed based on the insertion of this LTR-retrotransposon to genotype the I locus. A perfect linkage between bulb color phenotypes and marker genotypes was observed among 5,303 individuals of segregating populations. The transcription of AcMYB1 appeared to be normal in both red and white onions, but the transcription of CHS-A, which encodes chalcone synthase and is involved in the first step of the anthocyanin biosynthesis pathway, was inactivated in the white onions. Taken together, an aberrant AcMYB1 protein produced from the mutant allele might be responsible for the dominant white bulb color in onions.

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

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

An R2R3-MYB Transcriptional Factor LuMYB314 Associated with the Loss of Petal Pigmentation in Flax (Linum usitatissimum L.)

The loss of anthocyanin pigments is one of the most common evolutionary transitions in petal color, yet the genetic basis for these changes in flax remains largely unknown. In this study, we used crossing studies, a bulk segregant analysis, genome-wide association studies, a phylogenetic analysis, and transgenic testing to identify genes responsible for the transition from blue to white petals in flax. This study found no correspondence between the petal color and seed color, refuting the conclusion that a locus controlling the seed coat color is associated with the petal color, as reported in previous studies. The locus controlling the petal color was mapped using a BSA-seq analysis based on the F2 population. However, no significantly associated genomic regions were detected. Our genome-wide association study identified a highly significant QTL (BP4.1) on chromosome 4 associated with flax petal color in the natural population. The combination of a local Manhattan plot and an LD heat map identified LuMYB314, an R2R3-MYB transcription factor, as a potential gene responsible for the natural variations in petal color in flax. The overexpression of LuMYB314 in both Arabidopsis thaliana and Nicotiana tabacum resulted in anthocyanin deposition, indicating that LuMYB314 is a credible candidate gene for controlling the petal color in flax. Additionally, our study highlights the limitations of the BSA-seq method in low-linkage genomic regions, while also demonstrating the powerful detection capabilities of GWAS based on high-density genomic variation mapping. This study enhances our genetic insight into petal color variations and has potential breeding value for engineering LuMYB314 to develop colored petals, bast fibers, and seeds for multifunctional use in flax.

Jasmonic acid (JA)-mediating MYB transcription factor1, JMTF1, coordinates the balance between JA and auxin signalling in the rice defence response

The plant hormone jasmonic acid (JA) is a signalling compound involved in the regulation of cellular defence and development in plants. In this study, we investigated the roles of a JA-responsive MYB transcription factor, JMTF1, in the JA-regulated defence response against rice bacterial blight caused by Xanthomonas oryzae pv. oryzae (Xoo). JMTF1 did not interact with any JASMONATE ZIM-domain (JAZ) proteins. Transgenic rice plants overexpressing JMTF1 showed a JA-hypersensitive phenotype and enhanced resistance against Xoo. JMTF1 upregulated the expression of a peroxidase, OsPrx26, and monoterpene synthase, OsTPS24, which are involved in the biosynthesis of lignin and antibacterial monoterpene, γ-terpinene, respectively. OsPrx26 was mainly expressed in the vascular bundle. Transgenic rice plants overexpressing OsPrx26 showed enhanced resistance against Xoo. In addition to the JA-hypersensitive phenotype, the JMTF1-overexpressing rice plants showed a typical auxin-related phenotype. The leaf divergence and shoot gravitropic responses were defective, and the number of lateral roots decreased significantly in the JMTF1-overexpressing rice plants. JMTF1 downregulated the expression of auxin-responsive genes but upregulated the expression of OsIAA13, a suppressor of auxin signalling. The rice gain-of-function mutant Osiaa13 showed high resistance against Xoo. Transgenic rice plants overexpressing OsEXPA4, a JMTF1-downregulated auxin-responsive gene, showed increased susceptibility to Xoo. JMTF1 is selectively bound to the promoter of OsPrx26 in vivo. These results suggest that JMTF1 positively regulates disease resistance against Xoo by coordinating crosstalk between JA- and auxin-signalling in rice.

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 and Expression Analysis of the MYB Transcription Factor Family in Salvia nemorosa

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

Genome-Wide Identification and Characterization of MYB Gene Family and Analysis of Its Sex-Biased Expression Pattern in Spinacia oleracea L

The members of the myeloblastosis (MYB) family of transcription factors (TFs) participate in a variety of biological regulatory processes in plants, such as circadian rhythm, metabolism, and flower development. However, the characterization of MYB genes across the genomes of spinach Spinacia oleracea L. has not been reported. Here, we identified 140 MYB genes in spinach and described their characteristics using bioinformatics approaches. Among the MYB genes, 54 were 1R-MYB, 80 were 2R-MYB, 5 were 3R-MYB, and 1 was 4R-MYB. Almost all MYB genes were located in the 0-30 Mb region of autosomes; however, the 20 MYB genes were enriched at both ends of the sex chromosome (chromosome 4). Based on phylogeny, conserved motifs, and the structure of genes, 2R-MYB exhibited higher conservation relative to 1R-MYB genes. Tandem duplication and collinearity of spinach MYB genes drive their evolution, enabling the functional diversification of spinach genes. Subcellular localization prediction indicated that spinach MYB genes were mainly located in the nucleus. Cis-acting element analysis confirmed that MYB genes were involved in various processes of spinach growth and development, such as circadian rhythm, cell differentiation, and reproduction through hormone synthesis. Furthermore, through the transcriptome data analysis of male and female flower organs at five different periods, ten candidate genes showed biased expression in spinach males, suggesting that these genes might be related to the development of spinach anthers. Collectively, this study provides useful information for further investigating the function of MYB TFs and novel insights into the regulation of sex determination in spinach.

Plant-TFClass: a structural classification for plant transcription factors

Transcription factors (TFs) bind DNA at specific sequences to regulate gene expression. This universal process is achieved via their DNA-binding domain (DBD). In mammals, the vast diversity of DBD structural conformations and the way in which they contact DNA has been used to organize TFs in the TFClass hierarchical classification. However, the numerous DBD types present in plants but absent from mammalian genomes were missing from this classification. We reviewed DBD 3D structures and models available for plant TFs to classify most of the 56 recognized plant TF types within the TFClass framework. This extended classification adds eight new classes and 37 new families corresponding to DBD structures absent in mammals. Plant-TFClass provides a unique resource for TF comparison across families and organisms.

A suberized exodermis is required for tomato drought tolerance

Plant roots integrate environmental signals with development using exquisite spatiotemporal control. This is apparent in the deposition of suberin, an apoplastic diffusion barrier, which regulates flow of water, solutes and gases, and is environmentally plastic. Suberin is considered a hallmark of endodermal differentiation but is absent in the tomato endodermis. Instead, suberin is present in the exodermis, a cell type that is absent in the model organism Arabidopsis thaliana. Here we demonstrate that the suberin regulatory network has the same parts driving suberin production in the tomato exodermis and the Arabidopsis endodermis. Despite this co-option of network components, the network has undergone rewiring to drive distinct spatial expression and with distinct contributions of specific genes. Functional genetic analyses of the tomato MYB92 transcription factor and ASFT enzyme demonstrate the importance of exodermal suberin for a plant water-deficit response and that the exodermal barrier serves an equivalent function to that of the endodermis and can act in its place.

The Characterization of R2R3-MYB Genes in Ammopiptanthus nanus Uncovers That the miR858-AnaMYB87 Module Mediates the Accumulation of Anthocyanin under Osmotic Stress

R2R3-MYB transcription factors (TFs) participate in the modulation of plant development, secondary metabolism, and responses to environmental stresses. Ammopiptanthus nanus, a leguminous dryland shrub, tolerates a high degree of environmental stress, including drought and low-temperature stress. The systematic identification, structural analysis, evolutionary analysis, and gene profiling of R2R3-MYB TFs under cold and osmotic stress in A. nanus were performed. Up to 137 R2R3-MYB TFs were identified and clustered into nine clades, with most A. nanus R2R3-MYB members belonging to clade VIII. Tandem and segmental duplication events drove the expansion of the A. nanus R2R3-MYB family. Expression profiling revealed that multiple R2R3-MYB genes significantly changed under osmotic and cold stress conditions. MiR858 and miR159 targeted 88 R2R3-MYB genes. AnaMYB87, an miR858-targeted clade VIII R2R3-MYB TF, was up-regulated under both osmotic and cold stress. A transient expression assay in apples showed that the overexpression of AnaMYB87 promoted anthocyanin accumulation. A luciferase reporter assay in tobacco demonstrated that AnaMYB87 positively affected the transactivation of the dihydroflavonol reductase gene, indicating that the miR858-MYB87 module mediates anthocyanin accumulation under osmotic stress by regulating the dihydroflavonol reductase gene in A. nanus. This study provides new data to understand the roles of R2R3-MYB in plant stress responses.

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

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

Function of the R2R3-MYB Transcription Factors in Dalbergia odorifera and Their Relationship with Heartwood Formation

R2R3-MYB transcription factors (TFs) form one of the most important TF families involved in regulating various physiological functions in plants. The heartwood of Dalbergia odorifera is a kind of high-grade mahogany and valuable herbal medicine with wide application. However, the role of R2R3-MYB genes in the growth and development of D. odorifera, especially their relevance to heartwood formation, has not been revealed. A total of 126 R2R3-MYBs were screened from the D. odorifera genome and named DodMYB1-126 based on their location on 10 chromosomes. The collinearity results showed that purification selection was the main driving force for the evolution of the R2R3-MYB TFs family, and whole genome/fragment replication event was the main form for expanding the R2R3-MYB family, generating a divergence of gene structure and function. Comparative phylogenetic analysis classified the R2R3-MYB TFs into 33 subfamilies. S3-7,10,12-13,21 and N4-7 were extensively involved in the metabolic process; S9,13,16-19,24-25 and N1-3,8 were associated with the growth and development of D. odorifera. Based on the differential transcriptional expression levels of R2R3-MYBs in different tissues, DodMYB32, DodMYB55, and DodMYB89 were tentatively screened for involvement in the regulatory process of heartwood. Further studies have shown that the DodMYB89, localized in the nucleus, has transcriptional activation activity and is involved in regulating the biosynthesis of the secondary metabolites of heartwood by activating the promoters of the structural genes DodI2'H and DodCOMT. This study aimed to comprehensively analyze the functions of the R2R3-MYB TFs and screen for candidate genes that might be involved in heartwood formation of D. odorifera.

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.

Genome-Wide Comparative Analysis of the R2R3-MYB Gene Family in Six Ipomoea Species and the Identification of Anthocyanin-Related Members in Sweet Potatoes

Sweet potatoes (Ipomoea batatas) are one of the important tuberous root crops cultivated worldwide, and thier storage roots are rich in antioxidants, such as anthocyanins. R2R3-MYB is a large gene family involved in various biological processes, including anthocyanin biosynthesis. However, few reports about the R2R3-MYB gene family of sweet potatoes have been released to date. In the present study, a total of 695 typical R2R3-MYB genes were identified in six Ipomoea species, including 131 R2R3-MYB genes in sweet potatoes. A maximum likelihood phylogenetic analysis divided these genes into 36 clades, referring to the classification of 126 R2R3-MYB proteins of Arabidopsis. Clade C25(S12) has no members in six Ipomoea species, whereas four clades (i.e., clade C21, C26, C30, and C36), including 102 members, had no members in Arabidopsis, and they were identified as Ipomoea-specific clades. The identified R2R3-MYB genes were unevenly distributed on all chromosomes in six Ipomoea species genomes, and the collinearity analysis among hexaploid I. batatas and another five diploid Ipomoea species suggested that the sweet potato genome might have undergone a larger chromosome rearrangement during the evolution process. Further analyses of gene duplication events showed that whole-genome duplication, transposed duplication, and dispersed duplication events were the primary forces driving the R2R3-MYB gene family expansion of Ipomoea plants, and these duplicated genes experienced strong purifying selection because of their Ka/Ks ratio, which is less than 1. Additionally, the genomic sequence length of 131 IbR2R3-MYBs varied from 923 bp to ~12.9 kb with a mean of ~2.6 kb, and most of them had more than three exons. The Motif 1, 2, 3, and 4 formed typical R2 and R3 domains and were identified in all IbR2R3-MYB proteins. Finally, based on multiple RNA-seq datasets, two IbR2R3-MYB genes (IbMYB1/g17138.t1 and IbMYB113/g17108.t1) were relatively highly expressed in pigmented leaves and tuberous root flesh and skin, respectively; thus, they were identified to regulate tissue-specific anthocyanin accumulation in sweet potato. This study provides a basis for the evolution and function of the R2R3-MYB gene family in sweet potatoes and five other Ipomoea species.

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

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

Genome-wide analysis of R2R3-MYB genes in cultivated peanut (Arachis hypogaea L.): Gene duplications, functional conservation, and diversification

The cultivated Peanut (Arachis hypogaea L.), an important oilseed and edible legume, are widely grown worldwide. The R2R3-MYB transcription factor, one of the largest gene families in plants, is involved in various plant developmental processes and responds to multiple stresses. In this study we identified 196 typical R2R3-MYB genes in the genome of cultivated peanut. Comparative phylogenetic analysis with Arabidopsis divided them into 48 subgroups. The motif composition and gene structure independently supported the subgroup delineation. Collinearity analysis indicated polyploidization, tandem, and segmental duplication were the main driver of the R2R3-MYB gene amplification in peanut. Homologous gene pairs between the two subgroups showed tissue specific biased expression. In addition, a total of 90 R2R3-MYB genes showed significant differential expression levels in response to waterlogging stress. Furthermore, we identified an SNP located in the third exon region of AdMYB03-18 (AhMYB033) by association analysis, and the three haplotypes of the SNP were significantly correlated with total branch number (TBN), pod length (PL) and root-shoot ratio (RS ratio), respectively, revealing the potential function of AdMYB03-18 (AhMYB033) in improving peanut yield. Together, these studies provide evidence for functional diversity in the R2R3-MYB genes and will contribute to understanding the function of R2R3-MYB genes in peanut.

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

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

Chromosome-level reference genome of tetraploid Isoetes sinensis provides insights into evolution and adaption of lycophytes

BACKGROUND

The Lycophyta species are the extant taxa most similar to early vascular plants that were once abundant on Earth. However, their distribution has greatly diminished. So far, the absence of chromosome-level assembled lycophyte genomes has hindered our understanding of evolution and environmental adaption of lycophytes.

FINDINGS

We present the reference genome of the tetraploid aquatic quillwort, Isoetes sinensis, a lycophyte. This genome represents the first chromosome-level assembled genome of a tetraploid seed-free plant. Comparison of genomes between I. sinensis and Isoetestaiwanensis revealed conserved and different genomic features between diploid and polyploid lycophytes. Comparison of the I. sinensis genome with those of other species representing the evolutionary lineages of green plants revealed the inherited genetic tools for transcriptional regulation and most phytohormones in I. sinensis. The presence and absence of key genes related to development and stress responses provide insights into environmental adaption of lycophytes.

CONCLUSIONS

The high-quality reference genome and genomic analysis presented in this study are crucial for future genetic and environmental studies of not only I. sinensis but also other lycophytes.

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.

MYB Transcription Factors Becoming Mainstream in Plant Roots

The function of the root system is crucial for plant survival, such as anchoring plants, absorbing nutrients and water from the soil, and adapting to stress. MYB transcription factors constitute one of the largest transcription factor families in plant genomes with structural and functional diversifications. Members of this superfamily in plant development and cell differentiation, specialized metabolism, and biotic and abiotic stress processes are widely recognized, but their roles in plant roots are still not well characterized. Recent advances in functional studies remind us that MYB genes may have potentially key roles in roots. In this review, the current knowledge about the functions of MYB genes in roots was summarized, including promoting cell differentiation, regulating cell division through cell cycle, response to biotic and abiotic stresses (e.g., drought, salt stress, nutrient stress, light, gravity, and fungi), and mediate phytohormone signals. MYB genes from the same subfamily tend to regulate similar biological processes in roots in redundant but precise ways. Given their increasing known functions and wide expression profiles in roots, MYB genes are proposed as key components of the gene regulatory networks associated with distinct biological processes in roots. Further functional studies of MYB genes will provide an important basis for root regulatory mechanisms, enabling a more inclusive green revolution and sustainable agriculture to face the constant changes in climate and environmental conditions.

The Pea R2R3-MYB Gene Family and Its Role in Anthocyanin Biosynthesis in Flowers

Pea (Pisum sativum L.) is one of the most important legume crops in the world, and it has attracted great attention for its high nutritive values. Recently, the crop breeding program has been focused on the crop metabolic engineering (i.e., color, flavor, nutrition) to improve the quality of crop. As a major group of transcription factors forming the ternary MYB–bHLH–WD repeat protein (MBW) complex to regulate the anthocyanin biosynthesis pathway, members of R2R3-MYB gene family have always been the focus of research targets to improve the valuable metabolic product of crops. Until now, few report about the R2R3-MYB gene family of pea has been released. In this study, we identified 119 R2R3-MYB genes in the assembled pea genome (Version 1a), of which 111 were distributed across 14 chromosomes. Combining with the 126 R2R3-MYB protein sequences of Arabidopsis, we categorized 245 R2R3-MYB proteins into 36 subgroups according to sequence similarity and phylogenetic relationships. There was no member from subgroup 12, 15 and 29 existing in pea genome, whereas three novel subgroups were found in pea and named as N1-N3. Further analyses of conserved domains and Motifs, gene structures, and chromosomal locations showed that the typical R2 and R3 domains were present across all R2R3-MYB proteins, and Motif 1, 2, and 3 were identified in most members. Most of them had no more than two introns. Additionally, 119 pea R2R3-MYB genes did not experience large-scale duplication events. Finally, we concluded that several candidate genes may be responsible for the spatiotemporal accumulation of anthocyanins in pea petals. PsMYB116 was predominantly expressed in the dorsal petals to presumably activate the anthocyanin biosynthesis pathway, while PsMYB37 and PsMYB32 may positively regulates the anthocyanin accumulation in the lateral petals. This study not only provides a good reference to further characterize the diverse functions of R2R3-MYB genes but also helps researchers to understand the color formation of pea flowers.

Mutations in Rht-B1 Locus May Negatively Affect Frost Tolerance in Bread Wheat

The wheat semi-dwarfing genes Rht (Reduced height) are widely distributed among the contemporary wheat varieties. These genes also exert pleiotropic effects on plant tolerance towards various abiotic stressors. In this work, frost tolerance was studied in three near-isogenic lines of the facultative variety ‘April Bearded’ (AB), carrying the wild type allele Rht-B1a (tall phenotype), and the mutant alleles Rht-B1b (semi-dwarf) and Rht-B1c (dwarf), and was further compared with the tolerance of a typical winter type variety, ‘Mv Beres’. The level of freezing tolerance was decreasing in the order ‘Mv Beres’ > AB Rht-B1a > AB Rht-B1b > AB Rht-B1c. To explain the observed differences, cold acclimation-related processes were studied: the expression of six cold-related genes, the phenylpropanoid pathway, carbohydrates, amino acids, polyamines and compounds in the tricarboxylic acid cycle. To achieve this, a comprehensive approach was applied, involving targeted analyses and untargeted metabolomics screening with the help of gas chromatography/liquid chromatography—mass spectrometry setups. Several cold-related processes exhibited similar changes in these genotypes; indeed, the accumulation of eight putrescine and agmatine derivatives, 17 flavones and numerous oligosaccharides (max. degree of polymerization 18) was associated with the level of freezing tolerance in the ‘April Bearded’ lines. In summary, the mutant Rht alleles may further decrease the generally low frost tolerance of the Rht-B1a, and, based on the metabolomics study, the mechanisms of frost tolerance may differ for a typical winter variety and a facultative variety. Present results point to the complex nature of frost resistance.

In-depth comparative analysis of Tritrichomonas foetus transcriptomics reveals novel genes linked with adaptation to feline host

Tritrichomonas foetus is a flagellated parasite able to infect cattle, cats, and pigs. Despite its prevalence, feline tritrichomonosis has received markedly less attention than venereal infection, and little information about the molecular mechanisms that participate in feline host infection is available. Through a bioinformatics approach, we integrated public transcriptomic data for three T. foetus isolates and explored the differences at transcript level with a focus on pathogenesis and adaptation processes, particularly for the feline isolate. Our analysis revealed higher abundance levels of predicted virulence factors, such as proteases and surface antigens. Additionally, by a comparative and expression analysis of T. foetus genes, we proposed putative virulence factors that could be involved in feline infection. Finally, we identified a great proportion of predicted transcription factors of the MYB protein family and, by a promoter analysis, we revealed that MYB-related proteins could participate in the regulation of gene transcription in T. foetus. In conclusion, this integrated approach is a valuable resource for future studies of host–pathogen interactions and identifying new gene targets for improved feline tritrichomonosis diagnosis and treatment.

Comprehensive analysis of the MYB transcription factor gene family in Morus alba

BACKGROUND

The V-myb myeloblastosis viral oncogene homolog (MYB) family of proteins is large, containing functionally diverse transcription factors. However, MYBs in Morus are still poorly annotated and a comprehensive functional analysis of these transcription factors is lacking.

RESULTS

In the present study, a genome-wide identification of MYBs in Morus alba was performed. In total 166 MaMYBs were identified, including 103 R2R3-MYBs and four 3R-MaMYBs. Comprehensive analyses, including the phylogenetic analysis with putative functional annotation, motif and structure analysis, gene structure organization, promoter analysis, chromosomal localization, and syntenic relationships of R2R3-MaMYBs and 3R-MaMYBs, provided primary characterization for these MaMYBs. R2R3-MaMYBs covered the subgroups reported for R2R3-MYBs in Arabidopsis and Populus, and had two Morus-specific subgroups, indicating the high retention of MYBs in Morus. Motif analysis revealed high conservative residues at the start and end of each helix and residues consisting of the third helix in R2 and R3 repeats. Thirteen intron/exon patterns (a-m) were summarized, and the intron/exon pattern of two introns with phase numbers of 0 and 2 was the prevalent pattern for R2R3-MaMYBs. Various cis-elements in promoter regions were identified, and were mainly related to light response, development, phytohormone response, and abiotic and biotic stress response and secondary metabolite production. Expression patterns of R2R3-MaMYBs in different organs showed that MaMYBs involved in secondary cell wall components and stress responsiveness were preferentially expressed in roots or stems. R2R3-MaMYBs involved in flavonoid biosynthesis and anthocyanin accumulation were identified and characterized based on functional annotation and correlation of their expression levels with anthocyanin contents.

CONCLUSION

Based on a comprehensive analysis, this work provided functional annotation for R2R3-MYBs and an informative reference for further functional dissection of MYBs in Morus.

Genome-Wide Characterization of High-Affinity Nitrate Transporter 2 (NRT2) Gene Family in Brassica napus

Nitrate transporter 2 (NRT2) plays an essential role in Nitrogen (N) uptake, transport, utilization, and stress resistance. In this study, the NRT2 gene family in two sequenced Brassica napus ecotypes were identified, including 31 genes in ‘Zhongshuang11’ (BnaZSNRT2s) and 19 in ‘Darmor-bzh’ (BnaDarNRT2s). The candidate genes were divided into three groups (Group I−III) based on phylogenetic analyses, supported by a conserved intron-exon structure in each group. Collinearity analysis revealed that the large expansion of BnaZSNRT2s attributed to allopolyploidization of ancestors Brassica rapa and Brassica oleracea, and small-scale duplication events in B. napus. Transcription factor (TF) binding site prediction, cis-element analysis, and microRNA prediction suggested that the expressions of BnaZSNRT2s are regulated by multiple factors, and the regulatory pattern is relatively conserved in each group and is tightly connected between groups. Expression assay showed the diverse and differentiated spatial-temporal expression profiles of BnaZSNRT2s in Group I, but conserved patterns were observed in Group II/III; and the low nitrogen (LN) stress up-regulated expression profiles were presented in Group I−III, based on RNA-seq data. RT-qPCR analyses confirmed that BnaZSNRT2.5A-1 and BnaZSNRT2.5C-1 in Group II were highly up-regulated under LN stress in B. napus roots. Our results offer valid information and candidates for further functional BnaZSNRT2s studies.

Genome-wide identification and characterization of PdbHLH transcription factors related to anthocyanin biosynthesis in colored-leaf poplar (Populus deltoids)

Basic helix-loop-helix (bHLH) proteins are transcription factors (TFs) that have been shown to regulate anthocyanin biosynthesis in many plant species. However, the bHLH gene family in Populus deltoids has not yet been reported. In this study, 185 PdbHLH genes were identified in the Populus deltoids genome and were classified into 15 groups based on their sequence similarity and phylogenetic relationships. Analysis of the gene structure, chromosome location and conserved motif of the PdbHLH genes were performed by bioinformatic methods. Gene duplication analyses revealed that 114 PdbHLH were expanded and retained after WGD/segmental and proximal duplication. Investigation of cis-regulatory elements of PdbHLH genes indicated that many PdbHLH genes are involved in the regulation of anthocyanin biosynthesis. The expression patterns of PdbHLHs were obtained from previous data in two colored-leaf poplar (QHP and JHP) and green leaf poplar (L2025). Further analysis revealed that 12 candidate genes, including 3 genes (PdbHLH57, PdbHLH143, and PdbHLH173) from the subgroup III(f) and 9 gene from other groups, were positively associated with anthocyanin biosynthesis. In addition, 4 genes (PdbHLH4, PdbHLH1, PdbHLH18, and PdbHLH164) may be involved in negatively regulating the anthocyanin biosynthesis. These results provide a basis for the functional characterization of bHLH genes and investigations on the molecular mechanisms of anthocyanin biosynthesis in colored-leaf poplar.

Automatic identification and annotation of MYB gene family members in plants

BACKGROUND

MYBs are among the largest transcription factor families in plants. Consequently, members of this family are involved in a plethora of processes including development and specialized metabolism. The MYB families of many plant species were investigated in the last two decades since the first investigation looked at Arabidopsis thaliana. This body of knowledge and characterized sequences provide the basis for the identification, classification, and functional annotation of candidate sequences in new genome and transcriptome assemblies.

RESULTS

A pipeline for the automatic identification and functional annotation of MYBs in a given sequence data set was implemented in Python. MYB candidates are identified, screened for the presence of a MYB domain and other motifs, and finally placed in a phylogenetic context with well characterized sequences. In addition to technical benchmarking based on existing annotation, the transcriptome assembly of Croton tiglium and the annotated genome sequence of Castanea crenata were screened for MYBs. Results of both analyses are presented in this study to illustrate the potential of this application. The analysis of one species takes only a few minutes depending on the number of predicted sequences and the size of the MYB gene family. This pipeline, the required bait sequences, and reference sequences for a classification are freely available on github: https://github.com/bpucker/MYB_annotator .

CONCLUSIONS

This automatic annotation of the MYB gene family in novel assemblies makes genome-wide investigations consistent and paves the way for comparative studies in the future. Candidate genes for in-depth analyses are presented based on their orthology to previously characterized sequences which allows the functional annotation of the newly identified MYBs with high confidence. The identification of orthologs can also be harnessed to detect duplication and deletion events.

Evolution and functional diversification of R2R3-MYB transcription factors in plants

R2R3-MYB genes (R2R3-MYBs) form one of the largest transcription factor gene families in the plant kingdom, with substantial structural and functional diversity. However, the evolutionary processes leading to this amazing functional diversity have not yet been clearly established. Recently developed genomic and classical molecular technologies have provided detailed insights into the evolutionary relationships and functions of plant R2R3-MYBs. Here, we review recent genome-level and functional analyses of plant R2R3-MYBs, with an emphasis on their evolution and functional diversification. In land plants, this gene family underwent a large expansion by whole genome duplications and small-scale duplications. Along with this population explosion, a series of functionally conserved or lineage-specific subfamilies/groups arose with roles in three major plant-specific biological processes: development and cell differentiation, specialized metabolism, and biotic and abiotic stresses. The rapid expansion and functional diversification of plant R2R3-MYBs are highly consistent with the increasing complexity of angiosperms. In particular, recently derived R2R3-MYBs with three highly homologous intron patterns (a, b, and c) are disproportionately related to specialized metabolism and have become the predominant subfamilies in land plant genomes. The evolution of plant R2R3-MYBs is an active area of research, and further studies are expected to improve our understanding of the evolution and functional diversification of this gene family.

MYB gene family in the diatom Phaeodactylum tricornutum revealing their potential functions in the adaption to nitrogen deficiency and diurnal cycle

The MYB transcription factor (TF) family is one of the largest and most important TF families, regulating the growth and response of microalgae to stress. However, the gene structure and characteristics of Phaeodactylum tricornutum MYB TFs, and their functions under nitrogen deficiency, have not been explored yet. To identify all P. tricornutum MYB (PtMYB) genes, the MYB gene family was analyzed at the genome-wide level in this study. A total ofm26 PtMYB genes were identified from the genome of P. tricornutum. These PtMYB genes were divided into 5 subfamilies: 5R-MYB, 4R-MYB, R2R3-MYB, R1R2R3-MYB, and MYB-related proteins. Phylogenetical motif and gene structure analyses of MYB genes indicated that the number and proportion of MYB TFs were species-specific, and MYB genes exhibited a lot of duplication events from microalgae to higher plants. Furthermore, the differentially expressed patterns of all 26 PtMYB TFs implied that PtMYB genes might have functional specificity under nitrogen deficiency. Homology analysis of MYB genes revealed that PtMYB3, PtMYB15, and PtMYB21 might play important roles in the regulation of the diurnal cycle and response to nitrogen stress in P. tricornutum. PtMYB3, PtMYB15, and PtMYB21 genes might be used as potential candidate genes for further studying the regulatory mechanisms of P. tricornutum under nitrogen deficiency. This work provides an important foundation for the future research of the potential functions of PtMYB genes and its diurnal regulatory mechanisms under nitrogen deficiency.

Genome-Wide Identification and Expression Analysis of the R2R3-MYB Transcription Factor Family Revealed Their Potential Roles in the Flowering Process in Longan (Dimocarpus longan)

Longan (Dimocarpus longan Lour.) is a productive fruit crop with high nutritional and medical value in tropical and subtropical regions. The MYB gene family is one of the most widespread plant transcription factor (TF) families participating in the flowering regulation. However, little is known about the MYB TFs involved in the flowering process in longan and its regulatory network. In this study, a total of 119 DlR2R3-MYB genes were identified in the longan genome and were phylogenetically grouped into 28 subgroups. The groupings were supported by highly conserved gene structures and motif composition of DlR2R3-MYB genes in each subgroup. Collinearity analysis demonstrated that segmental replications played a more crucial role in the expansion of the DlR2R3-MYB gene family compared to tandem duplications, and all tandem/segmental duplication gene pairs have evolved under purifying selection. Interspecies synteny analysis among longan and five representative species implied the occurrence of gene duplication events was one of the reasons contributing to functional differentiation among species. RNA-seq data from various tissues showed DlR2R3-MYB genes displayed tissue-preferential expression patterns. The pathway of flower development was enriched with six DlR2R3-MYB genes. Cis-acting element prediction revealed the putative functions of DlR2R3-MYB genes were related to the plant development, phytohormones, and environmental stresses. Notably, the orthologous counterparts between Arabidopsis and longan R2R3-MYB members tended to play conserved roles in the flowering regulation and stress responses. Transcriptome profiling on off-season flower induction (FI) by KClO₃ indicated two up-regulated and four down-regulated DlR2R3-MYB genes involved in the response to KClO₃ treatment compared with control groups. Additionally, qRT-PCR confirmed certain genes exhibited high expression in flowers/flower buds. Subcellular localization experiments revealed that three predicted flowering-associated MYB proteins were localized in the nucleus. Future functional studies on these potential candidate genes involved in the flowering development could further the understanding of the flowering regulation mechanism.

R2R3-MYBs in Durum Wheat: Genome-Wide Identification, Poaceae-Specific Clusters, Expression, and Regulatory Dynamics Under Abiotic Stresses

MYB transcription factors (TFs) represent one of the biggest TF families in plants, being involved in various specific plant processes, such as responses to biotic and abiotic stresses. The implication of MYB TFs in the tolerance mechanisms to abiotic stress is particularly interesting for crop breeding, since environmental conditions can negatively affect growth and productivity. Wheat is a worldwide-cultivated cereal, and is a major source of plant-based proteins in human food. In particular, durum wheat plays an important role in global food security improvement, since its adaptation to hot and dry conditions constitutes the base for the success of wheat breeding programs in future. In the present study, a genome-wide identification of R2R3-MYB TFs in durum wheat was performed. MYB profile search and phylogenetic analyses based on homology with Arabidopsis and rice MYB TFs led to the identification of 233 R2R3-TdMYB (Triticum durum MYB). Three Poaceae-specific MYB clusters were detected, one of which had never been described before. The expression of eight selected genes under different abiotic stress conditions, revealed that most of them responded especially to salt and drought stress. Finally, gene regulatory network analyses led to the identification of 41 gene targets for three TdR2R3-MYBs that represent novel candidates for functional analyses. This study provides a detailed description of durum wheat R2R3-MYB genes and contributes to a deeper understanding of the molecular response of durum wheat to unfavorable climate conditions.

Full-length transcriptomic identification of R2R3-MYB family genes related to secondary cell wall development in Cunninghamia lanceolata (Chinese fir)

BACKGROUND

R2R3-MYB is a class of transcription factor crucial in regulating secondary cell wall development during wood formation. The regulation of wood formation in gymnosperm has been understudied due to its large genome size. Using Single-Molecule Real-Time sequencing, we obtained full-length transcriptomic libraries from the developmental stem of Cunninghamia lanceolata, a perennial conifer known as Chinese fir. The R2R3-MYB of C. lanceolata (hereafter named as ClMYB) associated with secondary wall development were identified based on phylogenetic analysis, expression studies and functional study on transgenic line.

RESULTS

The evolutionary relationship of 52 ClMYBs with those from Arabidopsis thaliana, Eucalyptus grandis, Populus trichocarpa, Oryza sativa, two gymnosperm species, Pinus taeda, and Picea glauca were established by neighbour-joining phylogenetic analysis. A large number of ClMYBs resided in the woody-expanded subgroups that predominated with the members from woody dicots. In contrast, the woody-preferential subgroup strictly carrying the members of woody dicots contained only one candidate. The results suggest that the woody-expanded subgroup emerges before the gymnosperm/angiosperm split, while most of the woody-preferential subgroups are likely lineage-specific to woody dicots. Nine candidates shared the same subgroups with the A. thaliana orthologs, with known function in regulating secondary wall development. Gene expression analysis inferred that ClMYB1/2/3/4/5/26/27/49/51 might participate in secondary wall development, among which ClMYB1/2/5/26/27/49 were significantly upregulated in the highly lignified compression wood region, reinforcing their regulatory role associated with secondary wall development. ClMYB1 was experimentally proven a transcriptional activator that localised in the nucleus. The overexpression of ClMYB1 in Nicotiana benthamiana resulted in an increased lignin deposition in the stems. The members of subgroup S4, ClMYB3/4/5 shared the ERF-associated amphiphilic repression motif with AtMYB4, which is known to repress the metabolism of phenylpropanoid derived compounds. They also carried a core motif specific to gymnosperm lineage, suggesting divergence of the regulatory process compared to the angiosperms.

CONCLUSIONS

This work will enrich the collection of full-length gymnosperm-specific R2R3-MYBs related to stem development and contribute to understanding their evolutionary relationship with angiosperm species.

Genome sequencing of the multicellular alga Astrephomene provides insights into convergent evolution of germ-soma differentiation

Germ-soma differentiation evolved independently in many eukaryotic lineages and contributed to complex multicellular organizations. However, the molecular genetic bases of such convergent evolution remain unresolved. Two multicellular volvocine green algae, Volvox and Astrephomene, exhibit convergent evolution of germ-soma differentiation. The complete genome sequence is now available for Volvox, while genome information is scarce for Astrephomene. Here, we generated the de novo whole genome sequence of Astrephomene gubernaculifera and conducted RNA-seq analysis of isolated somatic and reproductive cells. In Volvox, tandem duplication and neofunctionalization of the ancestral transcription factor gene (RLS1/rlsD) might have led to the evolution of regA, the master regulator for Volvox germ-soma differentiation. However, our genome data demonstrated that Astrephomene has not undergone tandem duplication of the RLS1/rlsD homolog or acquisition of a regA-like gene. Our RNA-seq analysis revealed the downregulation of photosynthetic and anabolic gene expression in Astrephomene somatic cells, as in Volvox. Among genes with high expression in somatic cells of Astrephomene, we identified three genes encoding putative transcription factors, which may regulate somatic cell differentiation. Thus, the convergent evolution of germ-soma differentiation in the volvocine algae may have occurred by the acquisition of different regulatory circuits that generate a similar division of labor.

Capsicum chinense MYB Transcription Factor Genes: Identification, Expression Analysis, and Their Conservation and Diversification With Other Solanaceae Genomes

Myeloblastosis (MYB) genes are important transcriptional regulators of plant growth, development, and secondary metabolic biosynthesis pathways, such as capsaicinoid biosynthesis in Capsicum. Although MYB genes have been identified in Capsicum annuum, no comprehensive study has been conducted on other Capsicum species. We identified a total of 251 and 240 MYB encoding genes in Capsicum chinense MYBs (CcMYBs) and Capsicum baccatum MYBs (CbMYBs). The observation of twenty tandem and 41 segmental duplication events indicated expansion of the MYB gene family in the C. chinense genome. Five CcMYB genes, i.e., CcMYB101, CcMYB46, CcMYB6, CcPHR8, and CcRVE5, and two CaMYBs, i.e., CaMYB3 and CaHHO1, were found within the previously reported capsaicinoid biosynthesis quantitative trait loci. Based on phylogenetic analysis with tomato MYB proteins, the Capsicum MYBs were classified into 24 subgroups supported by conserved amino acid motifs and gene structures. Also, a total of 241 CcMYBs were homologous with 225 C. annuum, 213 C. baccatum, 125 potato, 79 tomato, and 23 Arabidopsis MYBs. Synteny analysis showed that all 251 CcMYBs were collinear with C. annuum, C. baccatum, tomato, potato, and Arabidopsis MYBs spanning over 717 conserved syntenic segments. Using transcriptome data from three fruit developmental stages, a total of 54 CcMYBs and 81 CaMYBs showed significant differential expression patterns. Furthermore, the expression of 24 CcMYBs from the transcriptome data was validated by quantitative real-time (qRT) PCR analysis. Eight out of the 24 CcMYBs validated by the qRT-PCR were highly expressed in fiery hot C. chinense than in the lowly pungent C. annuum. Furthermore, the co-expression analysis revealed several MYB genes clustered with genes from the capsaicinoid, anthocyanin, phenylpropanoid, carotenoid, and flavonoids biosynthesis pathways, and related to determining fruit shape and size. The homology modeling of 126 R2R3 CcMYBs showed high similarity with that of the Arabidopsis R2R3 MYB domain template, suggesting their potential functional similarity at the proteome level. Furthermore, we have identified simple sequence repeat (SSR) motifs in the CcMYB genes, which could be used in Capsicum breeding programs. The functional roles of the identified CcMYBs could be studied further so that they can be manipulated for Capsicum trait improvement.

De novo transcriptome characterization of Iris atropurpurea (the Royal Iris) and phylogenetic analysis of MADS-box and R2R3-MYB gene families

The Royal Irises (section Oncocyclus) are a Middle-Eastern group of irises, characterized by extremely large flowers with a huge range of flower colors and a unique pollination system. The Royal Irises are considered to be in the course of speciation and serve as a model for evolutionary processes of speciation and pollination ecology. However, no transcriptomic and genomic data are available for these plants. Transcriptome sequencing is a valuable resource for determining the genetic basis of ecological-meaningful traits, especially in non-model organisms. Here we describe the de novo transcriptome assembly of Iris atropurpurea, an endangered species endemic to Israel's coastal plain. We sequenced and analyzed the transcriptomes of roots, leaves, and three stages of developing flower buds. To identify genes involved in developmental processes we generated phylogenetic gene trees for two major gene families, the MADS-box and MYB transcription factors, which play an important role in plant development. In addition, we identified 1503 short sequence repeats that can be developed for molecular markers for population genetics in irises. This first reported transcriptome for the Royal Irises, and the data generated, provide a valuable resource for this non-model plant that will facilitate gene discovery, functional genomic studies, and development of molecular markers in irises, to complete the intensive eco-evolutionary studies of this group.

Genome-Wide Characterization of R2R3-MYB Transcription Factors in Pitaya Reveals a R2R3-MYB Repressor HuMYB1 Involved in Fruit Ripening through Regulation of Betalain Biosynthesis by Repressing Betalain Biosynthesis-Related Genes

The MYB (myeloblastosis) superfamily constitutes one of the most abundant transcription factors (TFs) regulating various biological processes in plants. However, the molecular characteristics and functions of MYB TFs in pitaya remain unclear. To date, no genome-wide characterization analysis of this gene family has been conducted in the Cactaceae species. In this study, 105 R2R3-MYB members were identified from the genome data of Hylocereus undatus and their conserved motifs, physiological and biochemical characteristics, chromosome locations, synteny relationship, gene structure and phylogeny were further analyzed. Expression analyses suggested that three up-regulated HuMYBs and twenty-two down-regulated HuMYBs were probably involved in fruit ripening of pitaya. Phylogenetic analyses of R2R3-MYB repressors showed that seven HuMYBs (HuMYB1, HuMYB21, HuMYB48, HuMYB49, HuMYB72, HuMYB78 and HuMYB101) were in clades containing R2R3-MYB repressors. HuMYB1 and HuMYB21 were significantly down-regulated with the betalain accumulation during fruit ripening of ‘Guanhuahong’ pitaya (H. monacanthus). However, only HuMYB1 had R2 and R3 repeats with C1, C2, C3 and C4 motifs. HuMYB1 was localized exclusively to the nucleus and exhibited transcriptional inhibition capacities. Dual luciferase reporter assay demonstrated that HuMYB1 inhibited the expression of betalain-related genes: HuADH1, HuCYP76AD1-1 and HuDODA1. These results suggested that HuMYB1 is a potential repressor of betalain biosynthesis during pitaya fruit ripening. Our results provide the first genome-wide analyses of the R2R3-MYB subfamily involved in pitaya betalain biosynthesis and will facilitate functional analysis of this gene family in the future.

Innovation, conservation, and repurposing of gene function in root cell type development

Plant species have evolved myriads of solutions, including complex cell type development and regulation, to adapt to dynamic environments. To understand this cellular diversity, we profiled tomato root cell type translatomes. Using xylem differentiation in tomato, examples of functional innovation, repurposing, and conservation of transcription factors are described, relative to the model plant Arabidopsis. Repurposing and innovation of genes are further observed within an exodermis regulatory network and illustrate its function. Comparative translatome analyses of rice, tomato, and Arabidopsis cell populations suggest increased expression conservation of root meristems compared with other homologous populations. In addition, the functions of constitutively expressed genes are more conserved than those of cell type/tissue-enriched genes. These observations suggest that higher order properties of cell type and pan-cell type regulation are evolutionarily conserved between plants and animals.

Evolution of the Subgroup 6 R2R3-MYB Genes and Their Contribution to Floral Color in the Perianth-Bearing Piperales

Flavonoids, carotenoids, betalains, and chlorophylls are the plant pigments responsible for floral color. Anthocyanins, a class of flavonoids, are largely responsible for the red, purple, pink, and blue colors. R2R3-MYB genes belonging to subgroup 6 (SG6) are the upstream regulatory factors of the anthocyanin biosynthetic pathway. The canonical members of these genes in Arabidopsis include AtMYB75, AtMYB90, AtMYB113, and AtMYB114. The Aristolochiaceae is an angiosperm lineage with diverse floral groundplans and perianth colors. Saruma henryi exhibits a biseriate perianth with green sepals and yellow petals. All other genera have sepals only, with colors ranging from green (in Lactoris) to a plethora of yellow to red and purple mixtures. Here, we isolated and reconstructed the SG6 R2R3-MYB gene lineage evolution in angiosperms with sampling emphasis in Aristolochiaceae. We found numerous species-specific duplications of this gene lineage in core eudicots and local duplications in Aristolochiaceae for Saruma and Asarum. Expression of SG6 R2R3-MYB genes examined in different developmental stages and plant organs of four Aristolochiaceae species, largely overlaps with red and purple pigments, suggesting a role in anthocyanin and flavonoid synthesis and accumulation. A directed RNA-seq analysis corroborated our RT-PCR analyses, by showing that these structural enzymes activate during perianth development in Aristolochia fimbriata and that the regulatory genes are expressed in correlation with color phenotype. Finally, the reconstruction of the flavonoid and anthocyanin metabolic pathways using predicted peptides from transcriptomic data show that all pivotal enzymes are present in the analyzed species. We conclude that the regulatory genes as well as the biosynthetic pathway are largely conserved across angiosperms. In addition, the Aristolochiaceae emerges as a remarkable group to study the genetic regulatory network for floral color, as their members exhibit an outstanding floral diversity with elaborate color patterns and the genetic complement for SG6 R2R3-MYB genes is simpler than in core eudicot model species.

Functional and evolutionary analysis of the Arabidopsis 4R-MYB protein SNAPc4 as part of the SNAP complex

Transcription initiation of the genes coding for small nuclear RNA (snRNA) has been extensively analyzed in humans and fruit fly, but only a single ortholog of a snRNA-activating protein complex (SNAPc) subunit has so far been characterized in plants. The genome of the model plant Arabidopsis thaliana encodes orthologs of all three core SNAPc subunits, including A. thaliana SNAP complex 4 (AtSNAPc4)-a 4R-MYB-type protein with four-and-a-half adjacent MYB repeat units. We report the conserved role of AtSNAPc4 as subunit of a protein complex involved in snRNA gene transcription and present genetic evidence that AtSNAPc4 is an essential gene in gametophyte and zygote development. We present experimental evidence that the three A. thaliana SNAPc subunits assemble into a SNAP complex and demonstrate the binding of AtSNAPc4 to snRNA promoters. In addition, co-localization studies show a link between AtSNAPc4 accumulation and Cajal bodies, known to aggregate at snRNA gene loci in humans. Moreover, we show the strong evolutionary conservation of single-copy 4R-MYB/SNAPc4 genes in a broad range of eukaryotes and present additional shared protein features besides the MYB domain, suggesting a conservation of the snRNA transcription initiation machinery along the course of the eukaryotic evolution.

Overexpression of GmMYB14 improves high-density yield and drought tolerance of soybean through regulating plant architecture mediated by the brassinosteroid pathway

MYB transcription factors (TFs) have been reported to regulate the biosynthesis of secondary metabolites, as well as to mediate plant adaption to abiotic stresses, including drought. However, the roles of MYB TFs in regulating plant architecture and yield potential remain poorly understood. Here, we studied the roles of the dehydration-inducible GmMYB14 gene in regulating plant architecture, high-density yield and drought tolerance through the brassinosteroid (BR) pathway in soybean. GmMYB14 was shown to localize to nucleus and has a transactivation activity. Stable GmMYB14-overexpressing (GmMYB14-OX) transgenic soybean plants displayed a semi-dwarfism and compact plant architecture associated with decreased cell size, resulting in a decrease in plant height, internode length, leaf area, leaf petiole length and leaf petiole angle, and improved yield in high density under field conditions. Results of the transcriptome sequencing suggested the involvement of BRs in regulating GmMYB14-OX plant architecture. Indeed, GmMYB14-OX plants showed reduced endogenous BR contents, while exogenous application of brassinolide could partly rescue the phenotype of GmMYB14-OX plants. Furthermore, GmMYB14 was shown to directly bind to the promoter of GmBEN1 and up-regulate its expression, leading to reduced BR content in GmMYB14-OX plants. GmMYB14-OX plants also displayed improved drought tolerance under field conditions. GmBEN1 expression was also up-regulated in the leaves of GmMYB14-OX plants under polyethylene glycol treatment, indicating that the GmBEN1-mediated reduction in BR level under stress also contributed to drought/osmotic stress tolerance of the transgenic plants. Our findings provided a strategy for stably increasing high-density yield and drought tolerance in soybean using a single TF-encoding gene.

Genomic analysis uncovers functional variation in the C-terminus of anthocyanin-activating MYB transcription factors

MYB transcription factors regulate diverse aspects of plant development and secondary metabolism, often by partnering in transcriptional regulatory complexes. Here, we harness genomic resources to identify novel MYBs, thereby producing an updated eudicot MYB phylogeny with revised relationships among subgroups as well as new information on sequence variation in the disordered C-terminus of anthocyanin-activating MYBs. BLAST® and hidden Markov model scans of gene annotations identified a total of 714 MYB transcription factors across the genomes of four crops that span the eudicots: apple, grape, kiwifruit and tomato. Codon model-based phylogenetic inference identified novel members of previously defined subgroups, and the function of specific anthocyanin-activating subgroup 6 members was assayed transiently in tobacco leaves. Sequence conservation within subgroup 6 highlighted one previously described and two novel short linear motifs in the disordered C-terminal region. The novel motifs have a mix of hydrophobic and acidic residues and are predicted to be relatively ordered compared with flanking protein sequences. Comparison of motifs with the Eukaryotic Linear Motif database suggests roles in protein-protein interaction. Engineering of motifs and their flanking regions from strong anthocyanin activators into weak activators, and vice versa, affected function. We conclude that, although the MYB C-terminal sequence diverges greatly even within MYB clades, variation within the C-terminus at and near relatively ordered regions offers opportunities for exploring MYB function and developing superior alleles for plant breeding.

Nitric oxide signaling in the plant nucleus: the function of nitric oxide in chromatin modulation and transcription

Nitric oxide (NO) is involved in a vast number of physiologically important processes in plants, such as organ development, stress resistance, and immunity. Transduction of NO bioactivity is generally achieved by post-translational modification of proteins, with S-nitrosation of cysteine residues as the predominant form. While traditionally the subcellular location of the factors involved was of lesser importance, recent studies identified the connection between NO and transcriptional activity and thereby raised the question about the route of NO into the nuclear sphere. Identification of NO-affected transcription factors and chromatin-modifying histone deacetylases implicated the important role of NO signaling in the plant nucleus as a regulator of epigenetic mechanisms and gene transcription. Here, we discuss the relationship between NO and its directly regulated protein targets in the nuclear environment, focusing on S-nitrosated chromatin modulators and transcription factors.

AP2/ERF and R2R3-MYB family transcription factors: potential associations between temperature stress and lipid metabolism in Auxenochlorella protothecoides

BACKGROUND

Both APETALA2/Ethylene Responsive Factor (AP2/ERF) superfamily and R2R3-MYB family were from one of the largest diverse families of transcription factors (TFs) in plants, and played important roles in plant development and responses to various stresses. However, no systematic analysis of these TFs had been conducted in the green algae A. protothecoides heretofore. Temperature was a critical factor affecting growth and lipid metabolism of A. protothecoides. It also remained largely unknown whether these TFs would respond to temperature stress and be involved in controlling lipid metabolism process.

RESULTS

Hereby, a total of six AP2 TFs, six ERF TFs and six R2R3-MYB TFs were identified and their expression profiles were also analyzed under low-temperature (LT) and high-temperature (HT) stresses. Meanwhile, differential adjustments of lipid pathways were triggered, with enhanced triacylglycerol accumulation. A co-expression network was built between these 18 TFs and 32 lipid-metabolism-related genes, suggesting intrinsic associations between TFs and the regulatory mechanism of lipid metabolism.

CONCLUSIONS

This study represented an important first step towards identifying functions and roles of AP2 superfamily and R2R3-MYB family in lipid adjustments and response to temperature stress. These findings would facilitate the biotechnological development in microalgae-based biofuel production and the better understanding of photosynthetic organisms' adaptive mechanism to temperature stress.

Systematic investigation and expression profiles of the GbR2R3-MYB transcription factor family in ginkgo (Ginkgo biloba L.)

As one of the largest families of transcription factors, the R2R3-MYB family plays a significant role in plant growth, development, and response to hormone and environmental stress. To explore its evolutionary mechanism and potential function in Ginkgo biloba, a gymnosperm of great economic and ecological value, we presented a comprehensive analysis of the R2R3-MYB genes in ginkgo. Sixty-nine GbR2R3-MYB genes were identified and these genes could be classified into 33 groups based on the characteristics of the amino acid sequence of the R2R3-MYB domain and gene structure. Syntenic analyses indicated that few tandem and segmental duplications possibly resulted in the contraction of the GbR2R3-MYB gene family. Based on the transcriptome data, expression profiles of eight different tissues and different developmental stages of leaf and kernel showed that GbR2R3-MYB genes had distinct temporal and spatial expression characteristics. Specific expression patterns of the sixteen GbR2R3-MYB genes were also identified in response to different abiotic stresses and hormonal exposures. Further investigation revealed that GbR2R3-MYB19 was located in the nucleus and possessed transcriptional activity, implying its potential roles in the regulation of multiple biological processes. Our findings provide a robust basis for future comprehensive evolutionary and functional analyses of GbR2R3-MYB genes in ginkgo.

Evolution of A bHLH Interaction Motif

Intrinsically disordered proteins and regions with their associated short linear motifs play key roles in transcriptional regulation. The disordered MYC-interaction motif (MIM) mediates interactions between MYC and MYB transcription factors in Arabidopsis thaliana that are critical for constitutive and induced glucosinolate (GLS) biosynthesis. GLSs comprise a class of plant defense compounds that evolved in the ancestor of the Brassicales order. We used a diverse set of search strategies to discover additional occurrences of the MIM in other proteins and in other organisms and evaluate the findings by means of structural predictions, interaction assays, and biophysical experiments. Our search revealed numerous MIM instances spread throughout the angiosperm lineage. Experiments verify that several of the newly discovered MIM-containing proteins interact with MYC TFs. Only hits found within the same transcription factor family and having similar characteristics could be validated, indicating that structural predictions and sequence similarity are good indicators of whether the presence of a MIM mediates interaction. The experimentally validated MIMs are found in organisms outside the Brassicales order, showing that MIM function is broader than regulating GLS biosynthesis.

Genome-Wide Comparative Analysis of R2R3 MYB Gene Family in Populus and Salix and Identification of Male Flower Bud Development-Related Genes

The MYB transcription factor (TF) family is one of the largest plant transcription factor gene family playing vital roles in plant growth and development, including defense, cell differentiation, secondary metabolism, and responses to biotic and abiotic stresses. As a model tree species of woody plants, in recent years, the identification and functional prediction of certain MYB family members in the poplar genome have been reported. However, to date, the characterization of the gene family in the genome of the poplar's sister species willow has not been done, nor are the differences and similarities between the poplar and willow genomes understood. In this study, we conducted the first genome-wide investigation of the R2R3 MYB subfamily in the willow, identifying 216 R2R3 MYB gene members, and combined with the poplar R2R3 MYB genes, performed the first comparative analysis of R2R3 MYB genes between the poplar and willow. We identified 81 and 86 pairs of R2R3 MYB paralogs in the poplar and willow, respectively. There were 17 pairs of tandem repeat genes in the willow, indicating active duplication of willow R2R3 MYB genes. A further 166 pairs of poplar and willow orthologs were identified by collinear and synonymous analysis. The findings support the duplication of R2R3 MYB genes in the ancestral species, with most of the R2R3 MYB genes being retained during the evolutionary process. The phylogenetic trees of the R2R3 MYB genes of 10 different species were drawn. The functions of the poplar and willow R2R3 MYB genes were predicted using reported functional groupings and clustering by OrthoFinder. Identified 5 subgroups in general expanded in woody species, three subgroups were predicted to be related to lignin synthesis, and we further speculate that the other two subgroups also play a role in wood formation. We analyzed the expression patterns of the GAMYB gene of subgroup 18 (S18) related to pollen development in the male flower buds of poplar and willow at different developmental stages by qRT-PCR. The results showed that the GAMYB gene was specifically expressed in the male flower bud from pollen formation to maturity, and that the expression first increased and then decreased. Both the specificity of tissue expression specificity and conservation indicated that GAMYB played an important role in pollen development in both poplar and willow and was an ideal candidate gene for the analysis of male flower development-related functions of the two species.