Evolution of the 3R-MYB Gene Family in Plants

MYB transcription factors in plants: A comprehensive review of their discovery, structure, classification, functional diversity and regulatory mechanism

The MYB transcription factor (TF) family is one of the largest families in plants and performs highly diverse regulatory functions, particularly in relation to pathogen/pest resistance, nutrient/noxious substance absorption, drought/salt resistance, trichome growth, stamen development, leaf senescence, and flavonoid/terpenoid biosynthesis. Owing to their vital role in various biological regulatory processes, the mechanisms of MYB TFs have been extensively studied. Notably, MYB TFs not only directly regulate targets, such as phytohormones, reactive oxygen species signaling and secondary cell wall formation, but also serve as crucial points of crosstalk between these signaling networks. Here, we have comprehensively described the structures, classifications, and biological functions of MYB TFs, with a specific focus on their roles and mechanisms in the response to biotic and abiotic stresses, plant morphogenesis, and secondary metabolite biosynthesis. Different from other reported reviews, this review provides comprehensive knowledge on plant MYB TFs and will provide valuable insights in understanding regulatory networks and associated functions of plant MYB TFs to apply in resistance breeding and crop improvement.

Comparative transcriptional analysis of Persea americana MYB, WRKY and AP2/ERF transcription factors following Phytophthora cinnamomi infection

Plant cells undergo extensive transcriptional reprogramming following pathogen infection, with these reprogramming patterns becoming more complex when pathogens, such as hemibiotrophs, exhibit different lifestyles. These transcriptional changes are often orchestrated by MYB, WRKY and AP2/ERF transcription factors (TFs), which modulate both growth and defence-related gene expression. Transcriptional analysis of defence-related genes in avocado (Persea americana) infected with Phytophthora cinnamomi indicated differential immune response activation when comparing a partially resistant and susceptible rootstock. This study identified 226 MYB, 82 WRKY, and 174 AP2/ERF TF-encoding genes in avocado, using a genome-wide approach. Phylogenetic analysis revealed substantial sequence conservation within TF groups underscoring their functional significance. RNA-sequencing analysis in a partially resistant and susceptible avocado rootstock infected with P. cinnamomi was indicative of an immune response switch occurring in either rootstock after 24 and 6 h post-inoculation, respectively. Different clusters of co-expressed TF genes were observed at these times, suggesting the activation of necrotroph-related immune responses at varying intervals between the two rootstocks. This study aids our understanding of avocado immune response activation following P. cinnamomi infection, and the role of the TFs therein, elucidating the transcriptional reprogramming disparities between partially resistant and susceptible rootstocks.

The multifaceted roles of Myb domain-containing proteins in apicomplexan parasites

Apicomplexan parasites are a large and diverse clade of protists responsible for significant diseases of humans and animals. Central to the ability of these parasites to colonize their host and evade immune responses is an expanded repertoire of gene-expression programs that requires the coordinated action of complex transcriptional networks. DNA-binding proteins and chromatin regulators are essential orchestrators of apicomplexan gene expression that often act in concert. Although apicomplexan genomes encode various families of putative DNA-binding proteins, most remain functionally and mechanistically unexplored. This review highlights the versatile role of myeloblastosis (Myb) domain-containing proteins in apicomplexan parasites as transcription factors and chromatin regulators. We explore the diversity of Myb domain structure and use phylogenetic analysis to identify common features across the phylum. This provides a framework to discuss functional heterogeneity and regulation of Myb domain-containing proteins particularly emphasizing their role in parasite differentiation.

Characteristics and Functions of MYB (v-Myb avivan myoblastsis virus oncogene homolog)-Related Genes in Arabidopsis thaliana

The MYB (v-Myb avivan myoblastsis virus oncogene homolog) transcription factor family is one of the largest families of plant transcription factors which plays a vital role in many aspects of plant growth and development. MYB-related is a subclass of the MYB family. Fifty-nine Arabidopsis thaliana MYB-related (AtMYB-related) genes have been identified. In order to understand the functions of these genes, in this review, the promoters of AtMYB-related genes were analyzed by means of bioinformatics, and the progress of research into the functions of these genes has been described. The main functions of these AtMYB-related genes are light response and circadian rhythm regulation, root hair and trichome development, telomere DNA binding, and hormone response. From an analysis of cis-acting elements, it was found that the promoters of these genes contained light-responsive elements and plant hormone response elements. Most genes contained elements related to drought, low temperature, and defense and stress responses. These analyses suggest that AtMYB-related genes may be involved in A. thaliana growth and development, and environmental adaptation through plant hormone pathways. However, the functions of many genes do not occur independently but instead interact with each other through different pathways. In the future, the study of the role of the gene in different pathways will be conducive to a comprehensive understanding of the function of the gene. Therefore, gene cloning and protein functional analyses can be subsequently used to understand the regulatory mechanisms of AtMYB-related genes in the interaction of multiple signal pathways. This review provides theoretical guidance for the follow-up study of plant MYB-related genes.

Identification of the passion fruit (Passiflora edulis Sims) MYB family in fruit development and abiotic stress, and functional analysis of PeMYB87 in abiotic stresses

Environmental stresses are ubiquitous in agricultural cultivation, and they affect the healthy growth and development of edible tissues in passion fruit. The study of resistance mechanisms is important in understanding the adaptation and resistance of plants to environmental stresses. In this work, two differently resistant passion fruit varieties were selected, using the expression characteristics of the transcription factor MYB, to explore the resistance mechanism of the MYB gene under various environmental stresses. A total of 174 MYB family members were identified using high-quality passion fruit genomes: 98 2R-MYB, 5 3R-MYB, and 71 1R-MYB (MYB-relate). Their family information was systematically analyzed, including subcellular localization, physicochemical properties, phylogeny at the genomic level, promoter function, encoded proteins, and reciprocal regulation. In this study, bioinformatics and transcriptome sequencing were used to identify members of the PeMYB genes in passion fruit whole-genome data, and biological techniques, such as qPCR, gene clone, and transient transformation of yeast, were used to determine the function of the passion fruit MYB genes in abiotic stress tolerance. Transcriptomic data were obtained for differential expression characteristics of two resistant and susceptible varieties, three expression patterns during pulp development, and four induced expression patterns under abiotic stress conditions. We further focused on the resistance mechanism of PeMYB87 in environmental stress, and we selected 10 representative PeMYB genes for quantitative expression verification. Most of the genes were differentially induced by four abiotic stresses, among which PeMYB87 responded significantly to high-temperature-induced expression and overexpression of the PeMYB87 gene in the yeast system. The transgenic PeMYB87 in yeast showed different degrees of stress resistance under exposure to cold, high temperatures, drought, and salt stresses. These findings lay the foundation for further analysis of the biological functions of PeMYBs involved in stress resistance in passion fruit.

Phylogenetic Analysis of R2R3-MYB Family Genes in Tetrastigma hemsleyanum Diels et Gilg and Roles of ThMYB4 and ThMYB7 in Flavonoid Biosynthesis

Tetrastigma hemsleyanum Diels et Gilg (T. hemsleyanum) is an extensively used Chinese folk herb with multiple bioactivities. Among these bioactivities, flavonoids are recognized as the representative active ingredients. We previously found an elevated accumulation of flavonoids in T. hemsleyanum under water stress; however, the mechanism remains unclear. R2R3-MYB transcription factors play vital roles in the plant response to environmental stress and the regulation of secondary metabolites. Herein, a systematic transcriptome identification of R2R3-MYB family genes under water stress in T. hemsleyanum was performed to explore their potential function in the biosynthesis of flavonoids. A total of 26 R2R3-MYB genes were identified, most of which were clustered into functional branches of abiotic stress. ThMYB4 and ThMYB7 were then screened out to be associated with the biosynthesis of flavonoids through a protein-protein interaction prediction. An expression correlation analysis based on RNA-seq further confirmed that ThMYB4 and ThMYB7 were positively related to the flavonoid biosynthetic pathway genes of T. hemsleyanum. In ThMYB4- and ThMYB7-overexpression hairy roots, it was found that the expression of ThCHS and ThCHI was significantly increased, suggesting that ThMYB4 and ThMYB7 may act as regulators in flavonoid biosynthesis. This will shed new light on the promotion of flavonoid production and the medicinal value of T. hemsleyanum by manipulating transcription factors.

Rapid evolution of T2/S-RNase genes in Fragaria linked to multiple transitions from self-incompatibility to self-compatibility

The T2/RNase gene family is widespread in eukaryotes, and particular members of this family play critical roles in the gametophytic self-incompatibility (GSI) system in plants. Wild diploid strawberry (Fragaria) species have diversified their sexual systems via self-incompatible and self-compatible traits, yet how these traits evolved in Fragaria remains elusive. By integrating the published and de novo assembled genomes and the newly generated RNA-seq data, members of the RNase T2 gene family were systematically identified in six Fragaria species, including three self-incompatible species (Fragaria nipponica, Fragaria nubicola, and Fragaria viridis) and three self-compatible species (Fragaria nilgerrensis, Fragaria vesca, and Fragaria iinumae). In total, 115 RNase T2 genes were identified in the six Fragaria genomes and can be classified into three classes (I-III) according to phylogenetic analysis. The identified RNase T2 genes could be divided into 22 homologous gene sets according to amino acid sequence similarity and phylogenetic and syntenic relationships. We found that extensive gene loss and pseudogenization coupled with small-scale duplications mainly accounted for variations in the RNase T2 gene numbers in Fragaria. Multiple copies of homologous genes were mainly generated from tandem and segmental duplication events. Furthermore, we newly identified five S-RNase genes in three self-incompatible Fragaria genomes, including two in F. nipponica, two in F. viridis, and one in F. nubicola, which fit for typical features of a pistil determinant, including highly pistil-specific expression, highly polymorphic proteins and alkaline isoelectric point (pI), while no S-RNase genes were found in all three self-compatible Fragaria species. Surprisingly, these T2/S-RNase genes contain at least one large intron (>10 kb). This study revealed that the rapid evolution of T2/S-RNase genes within the Fragaria genus could be associated with its sexual mode, and repeated evolution of the self-compatible traits in Fragaria was convergent via losses of S-RNase.

MYB Transcription Factor Family in Pearl Millet: Genome-Wide Identification, Evolutionary Progression and Expression Analysis under Abiotic Stress and Phytohormone Treatments

Transcription factors (TFs) are the regulatory proteins that act as molecular switches in controlling stress-responsive gene expression. Among them, the MYB transcription factor family is one of the largest TF family in plants, playing a significant role in plant growth, development, phytohormone signaling and stress-responsive processes. Pearl millet (Pennisetum glaucum L.) is one of the most important C4 crop plants of the arid and semi-arid regions of Africa and Southeast Asia for sustaining food and fodder production. To explore the evolutionary mechanism and functional diversity of the MYB family in pearl millet, we conducted a comprehensive genome-wide survey and identified 279 MYB TFs (PgMYB) in pearl millet, distributed unevenly across seven chromosomes of pearl millet. A phylogenetic analysis of the identified PgMYBs classified them into 18 subgroups, and members of the same group showed a similar gene structure and conserved motif/s pattern. Further, duplication events were identified in pearl millet that indicated towards evolutionary progression and expansion of the MYB family. Transcriptome data and relative expression analysis by qRT-PCR identified differentially expressed candidate PgMYBs (PgMYB2, PgMYB9, PgMYB88 and PgMYB151) under dehydration, salinity, heat stress and phytohormone (ABA, SA and MeJA) treatment. Taken together, this study provides valuable information for a prospective functional characterization of the MYB family members of pearl millet and their application in the genetic improvement of crop plants.

Comparative Transcriptome Analysis Identifies Key Defense Genes and Mechanisms in Mulberry (Morus alba) Leaves against Silkworms (Bombyx mori)

As a consequence of long-term coevolution and natural selection, the leaves of mulberry (Morus alba) trees have become the best food source for silkworms (Bombyx mori). Nevertheless, the molecular and genomic basis of defense response remains largely unexplored. In the present study, we assessed changes in the transcriptome changes of mulberry in response to silkworm larval feeding at 0, 3, and 6 h. A total of 4709 (up = 2971, down = 1738) and 3009 (up = 1868, down = 1141) unigenes were identified after 3 and 6 h of silkworm infestation, respectively. MapMan enrichment analysis results show structural traits such as leaf surface wax, cell wall thickness and lignification form the first physical barrier to feeding by the silkworms. Cluster analysis revealed six unique temporal patterns of transcriptome changes. We predicted that mulberry promoted rapid changes in signaling and other regulatory processes to deal with mechanical damage, photosynthesis impairment, and other injury caused by herbivores within 3–6 h. LRR-RK coding genes (THE1, FER) was predicted participated in perception of cell wall perturbation in mulberry responding to silkworm feeding. Ca²⁺ signal sensors (CMLs), ROS (OST1, SOS3), RBOHD/F, CDPKs, and ABA were part of the regulatory network after silkworm feeding. Jasmonic acid (JA) signal transduction was predicted to act in silkworm feeding response, 10 JA signaling genes (such as OPR3, JAR1, and JAZ1) and 21 JA synthesis genes (such as LOX2, AOS, and ACX1) were upregulated after silkworm feeding for 3 h. Besides, genes of “alpha-Linolenic acid metabolism” and “phenylpropanoid biosynthesis” were activated in 3 h to reprogram secondary metabolism. Collectively, these findings provided valuable insights into silkworm herbivory-induced regulatory and metabolic processes in mulberry, which might help improve the coevolution of silkworm and mulberry.

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.

Comprehensive analysis of 1R- and 2R-MYBs reveals novel genic and protein features, complex organisation, selective expansion and insights into evolutionary tendencies

Myeloblastosis (MYB) family, the largest plant transcription factor family, has been subcategorised based on the number and type of repeats in the MYB domain. In spite of several reports, evolution of MYB genes and repeats remains enigmatic. Brassicaceae members are endowed with complex genomes, including dysploidy because of its unique history with multiple rounds of polyploidisation, genomic fractionations and rearrangements. The present study is an attempt to gain insights into the complexities of MYB family diversity, understand impacts of genome evolution on gene families and develop an evolutionary framework to understand the origin of various subcategories of MYB gene family. We identified and analysed 1129 MYBs that included 1R-, 2R-, 3R- and atypical-MYBs across sixteen species representing protists, fungi, animals and plants and exclude MYB identified from Brassicaceae except Arabidopsis thaliana; in addition, a total of 1137 2R-MYB genes from six Brassicaceae species were also analysed. Comparative analysis revealed predominance of 1R-MYBs in protists, fungi, animals and lower plants. Phylogenetic reconstruction and analysis of selection pressure suggested ancestral nature of R1-type repeat containing 1R-MYBs that might have undergone intragenic duplication to form multi-repeat MYBs. Distinct differences in gene structure between 1R-MYB and 2R-MYBs were observed regarding intron number, the ratio of gene length to coding DNA sequence (CDS) length and the length of exons encoding the MYB domain. Conserved as well as novel and lineage-specific intron phases were identified. Analyses of physicochemical properties revealed drastic differences indicating functional diversification in MYBs. Phylogenetic reconstruction of 1R- and 2R-MYB genes revealed a shared structure-function relationship in clades which was supported when transcriptome data was analysed in silico. Comparative genomics to study distribution pattern and mapping of 2R-MYBs revealed congruency and greater degree of synteny and collinearity among closely related species. Micro-synteny analysis of genomic segments revealed high conservation of genes that are immediately flanking the surrounding tandemly organised 2R-MYBs along with instances of local duplication, reorganisations and genome fractionation. In summary, polyploidy, dysploidy, reshuffling and genome fractionation were found to cause loss or gain of 2R-MYB genes. The findings need to be supported with functional validation to understand gene structure-function relationship along the evolutionary lineage and adaptive strategies based on comparative functional genomics in plants.

Analysis of Plant-Plant Interactions Reveals the Presence of Potent Antileukemic Compounds

A method to identify anticancer compounds in plants was proposed based on the hypothesis that these compounds are primarily present in plants to provide them with an ecological advantage over neighboring plants and other competitors. According to this view, identifying plants that contain compounds that inhibit or interfere with the development of other plant species may facilitate the discovery of novel anticancer agents. The method was developed and tested using Magnolia grandiflora, Gynoxys verrucosa, Picradeniopsis oppositifolia, and Hedyosmum racemosum, which are plant species known to possess compounds with cytotoxic activities. Plant extracts were screened for growth inhibitory activity, and then a thin-layer chromatography bioautography assay was conducted. This located the major antileukemic compounds 1, 2, 4, and 5 in the extracts. Once the active compounds were located, they were extracted and purified, and their structures were determined. The growth inhibitory activity of the purified compounds showed a significant correlation with their antileukemic activity. The proposed approach is rapid, inexpensive, and can easily be implemented in areas of the world with high biodiversity but with less access to advanced facilities and biological assays.

The SINGLE FLOWER (SFL) gene encodes a MYB transcription factor that regulates the number of flowers produced by the inflorescence of chickpea

Legumes usually have compound inflorescences, where flowers/pods develop from secondary inflorescences (I2), formed laterally at the primary inflorescence (I1). Number of flowers per I2, characteristic of each legume species, has important ecological and evolutionary relevance as it determines diversity in inflorescence architecture; moreover, it is also agronomically important for its potential impact on yield. Nevertheless, the genetic network controlling the number of flowers per I2 is virtually unknown. Chickpea (Cicer arietinum) typically produces one flower per I2 but single flower (sfl) mutants produce two (double-pod phenotype). We isolated the SFL gene by mapping the sfl-d mutation and identifying and characterising a second mutant allele. We analysed the effect of sfl on chickpea inflorescence ontogeny with scanning electron microscopy and studied the expression of SFL and meristem identity genes by RNA in situ hybridisation. We show that SFL corresponds to CaRAX1/2a, which codes a MYB transcription factor specifically expressed in the I2 meristem. Our findings reveal SFL as a central factor controlling chickpea inflorescence architecture, acting in the I2 meristem to regulate the length of the period for which it remains active, and therefore determining the number of floral meristems that it can produce.

Genome-wide Identification and Expression Analysis of RcMYB Genes in Rhodiola crenulata

Modern research has proved that the main medicinal component of Rhodiola crenulata, which has a wide range of medicinal value, is its secondary metabolite salidroside. The MYB transcription factor family is widely involved in biosynthesis of second metabolism and other roles in the stress response in plants, so a genome-wide identification and analysis for this family in R. crenulata is worth conducting. In this research, genome-wide analysis identified 139 MYB genes based on conserved domains in the R. crenulata genome, and 137 genes were used to construct a phylogenetic tree and modified with expression files to reveal evolutionary characteristics. Physical and chemical characteristics, gene structure, and conserved motif analysis were also used to further analyze RcMYBs. Additionally, cis-acting elements related to transcription, hormone, and MYB binding were found in the promoter region of the selected RcMYBs. Four RcMYBs were cloned, sequenced, and their gene expression pattern was analyzed for further analysis of their functions. The research results lay the foundation for further research on the function of RcMYB and R. crenulata.

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.

Independent Evolution of the MYB Family in Brown Algae

Myeloblastosis (MYB) proteins represent one of the largest families of eukaryotic transcription factors and regulate important processes in growth and development. Studies on MYBs have mainly focused on animals and plants; however, comprehensive analysis across other supergroups such as SAR (stramenopiles, alveolates, and rhizarians) is lacking. This study characterized the structure, evolution, and expression of MYBs in four brown algae, which comprise the biggest multicellular lineage of SAR. Subfamily 1R-MYB comprised heterogeneous proteins, with fewer conserved motifs found outside the MYB domain. Unlike the SHAQKY subgroup of plant 1R-MYB, THAQKY comprised the largest subgroup of brown algal 1R-MYBs. Unlike the expansion of 2R-MYBs in plants, brown algae harbored more 3R-MYBs than 2R-MYBs. At least ten 2R-MYBs, fifteen 3R-MYBs, and one 6R-MYB orthologs existed in the common ancestor of brown algae. Phylogenetic analysis showed that brown algal MYBs had ancient origins and a diverged evolution. They showed strong affinity with stramenopile species, while not with red algae, green algae, or animals, suggesting that brown algal MYBs did not come from the secondary endosymbiosis of red and green plastids. Sequence comparison among all repeats of the three types of MYB subfamilies revealed that the repeat of 1R-MYBs showed higher sequence identity with the R3 of 2R-MYBs and 3R-MYBs, which supports the idea that 1R-MYB was derived from loss of the first and second repeats of the ancestor MYB. Compared with other species of SAR, brown algal MYB proteins exhibited a higher proportion of intrinsic disordered regions, which might contribute to multicellular evolution. Expression analysis showed that many MYB genes are responsive to different stress conditions and developmental stages. The evolution and expression analyses provided a comprehensive analysis of the phylogeny and functions of MYBs in brown algae.

Plant protein-coding gene families: Their origin and evolution

Steady advances in genome sequencing methods have provided valuable insights into the evolutionary processes of several gene families in plants. At the core of plant biodiversity is an extensive genetic diversity with functional divergence and expansion of genes across gene families, representing unique phenomena. The evolution of gene families underpins the evolutionary history and development of plants and is the subject of this review. We discuss the implications of the molecular evolution of gene families in plants, as well as the potential contributions, challenges, and strategies associated with investigating phenotypic alterations to explain the origin of plants and their tolerance to environmental stresses.

Evolution of 14-3-3 Proteins in Angiosperm Plants: Recurring Gene Duplication and Loss

14-3-3 proteins are key regulatory factors in plants and are involved in a broad range of physiological processes. We addressed the evolutionary history of 14-3-3s from 46 angiosperm species, including basal angiosperm Amborella and major lineage of monocotyledons and eudicotyledons. Orthologs of Arabidopsis isoforms were detected. There were several rounds of duplication events in the evolutionary history of the 14-3-3 protein family in plants. At least four subfamilies (iota, epsilon, kappa, and psi) formed as a result of ancient duplication in a common ancestor of angiosperm plants. Recent duplication events followed by gene loss in plant lineage, among others Brassicaceae, Fabaceae, and Poaceae, further shaped the high diversity of 14-3-3 isoforms in plants. Coexpression data showed that 14-3-3 proteins formed different functional groups in different species. In some species, evolutionarily related groups of 14-3-3 proteins had coexpressed together under certain physiological conditions, whereas in other species, closely related isoforms expressed in the opposite manner. A possible explanation is that gene duplication and loss is accompanied by functional plasticity of 14-3-3 proteins.

Cis-acting elements involved in the G2/M-phase-specific transcription of the cyclin B gene in the unicellular alga Cyanidioschyzon merolae

M-specific activator (MSA) cis-acting elements have been determined to be involved in the regulation of G2/M-phase-specific transcription in spermatophytes. In this study, the involvement of MSA-core elements in G2/M-phase-specific transcription was examined in the unicellular red alga Cyanidioschyzon merolae. In the C. merolae genome, MSA-core elements do not accumulate specifically in the upstream of mitosis-specific transcriptional genes. Mutations of the four MSA-core elements of the cyclin B gene, which encodes a central factor of the G2-to-M-phase transition, have resulted in the abolishment of transcription or permission of transcription even in the G1 phase. These results suggest that all four MSA-core elements located in the upstream region of cyclin B are involved in G2/M-phase-specific transcription in C. merolae; however, the nature of the involvement of MSA-core elements in G2/M-phase-specific transcription differed among the four elements. Thus, MSA-core-element-mediated G2/M-phase-specific transcription in C. merolae seems to be regulated by a complex mechanism.

Identification of two tobacco genes encoding MYB3R proteins with repressor function and showing cell cycle-regulated transcript accumulation

MYB3R family transcription factors play a central role in the regulation of G2/M-specific gene transcription in Arabidopsis thaliana. Among the members of this family, MYB3R3 and MYB3R5 are structurally closely related and are involved in the transcriptional repression of target genes in both proliferating and quiescent cells. This type of MYB3R repressor is widespread in plants; however, apart from the studies on MYB3Rs in Arabidopsis thaliana, little information about them is available. Here we isolated tobacco cDNA clones encoding two closely related MYB3R proteins designated as NtmybC1 and NtmybC2 and determined the nucleotide sequences of the entire coding regions. Phylogenetic analysis suggested that NtmybC1 and NtmybC2 can be grouped into a conserved subfamily of plant MYB3Rs that also contains MYB3R3 and MYB3R5. When transiently expressed in protoplasts prepared from tobacco BY-2 cells, NtmybC1 and NtmybC2 repressed the activity of target promoters and blocked promoter activation mediated by NtmybA2, a MYB3R activator from tobacco. Unlike MYB3R3 and MYB3R5, NtmybC1 and NtmybC2 showed cell cycle-regulated transcript accumulation. In synchronized cultures of BY-2 cells, mRNAs for both NtmybC1 and NtmybC2 were preferentially expressed during the G2 and M phases, coinciding with the expression of NtmybA2 and G2/M-specific target genes. These results not only broadly confirm our fundamental view that this type of MYB3R protein acts as transcriptional repressor of G2/M-specific genes but also suggest a possible divergence of MYB3R repressors in terms of the mechanisms of their action and regulation.

Regulation of the Plant Cell Cycle in Response to Hormones and the Environment

Developmental and environmental signals converge on cell cycle machinery to achieve proper and flexible organogenesis under changing environments. Studies on the plant cell cycle began 30 years ago, and accumulated research has revealed many links between internal and external factors and the cell cycle. In this review, we focus on how phytohormones and environmental signals regulate the cell cycle to enable plants to cope with a fluctuating environment. After introducing key cell cycle regulators, we first discuss how phytohormones and their synergy are important for regulating cell cycle progression and how environmental factors positively and negatively affect cell division. We then focus on the well-studied example of stress-induced G2 arrest and view the current model from an evolutionary perspective. Finally, we discuss the mechanisms controlling the transition from the mitotic cycle to the endocycle, which greatly contributes to cell enlargement and resultant organ growth in plants.

Comprehensive analysis of SSRs and database construction using all complete gene-coding sequences in major horticultural and representative plants

Simple sequence repeats (SSRs) are one of the most important genetic markers and widely exist in most species. Here, we identified 249,822 SSRs from 3,951,919 genes in 112 plants. Then, we conducted a comprehensive analysis of these SSRs and constructed a plant SSR database (PSSRD). Interestingly, more SSRs were found in lower plants than in higher plants, showing that lower plants needed to adapt to early extreme environments. Four specific enriched functional terms in the lower plant Chlamydomonas reinhardtii were detected when it was compared with seven other higher plants. In addition, Guanylate_cyc existed in more genes of lower plants than of higher plants. In our PSSRD, we constructed an interactive plotting function in the chart interface, and users can easily view the detailed information of SSRs. All SSR information, including sequences, primers, and annotations, can be downloaded from our database. Moreover, we developed Web SSR Finder and Batch SSR Finder tools, which can be easily used for identifying SSRs. Our database was developed using PHP, HTML, JavaScript, and MySQL, which are freely available at http://www.pssrd.info/ . We conducted an analysis of the Myb gene families and flowering genes as two applications of the PSSRD. Further analysis indicated that whole-genome duplication and whole-genome triplication played a major role in the expansion of the Myb gene families. These SSR markers in our database will greatly facilitate comparative genomics and functional genomics studies in the future.

MYB-Mediated Regulation of Anthocyanin Biosynthesis

Anthocyanins are natural water-soluble pigments that are important in plants because they endow a variety of colors to vegetative tissues and reproductive plant organs, mainly ranging from red to purple and blue. The colors regulated by anthocyanins give plants different visual effects through different biosynthetic pathways that provide pigmentation for flowers, fruits and seeds to attract pollinators and seed dispersers. The biosynthesis of anthocyanins is genetically determined by structural and regulatory genes. MYB (v-myb avian myeloblastosis viral oncogene homolog) proteins are important transcriptional regulators that play important roles in the regulation of plant secondary metabolism. MYB transcription factors (TFs) occupy a dominant position in the regulatory network of anthocyanin biosynthesis. The TF conserved binding motifs can be combined with other TFs to regulate the enrichment and sedimentation of anthocyanins. In this study, the regulation of anthocyanin biosynthetic mechanisms of MYB-TFs are discussed. The role of the environment in the control of the anthocyanin biosynthesis network is summarized, the complex formation of anthocyanins and the mechanism of environment-induced anthocyanin synthesis are analyzed. Some prospects for MYB-TF to modulate the comprehensive regulation of anthocyanins are put forward, to provide a more relevant basis for further research in this field, and to guide the directed genetic modification of anthocyanins for the improvement of crops for food quality, nutrition and human health.

MYB3R-mediated active repression of cell cycle and growth under salt stress in Arabidopsis thaliana

Under environmental stress, plants are believed to actively repress their growth to save resource and alter its allocation to acquire tolerance against the stress. Although a lot of studies have uncovered precise mechanisms for responding to stress and acquiring tolerance, the mechanisms for regulating growth repression under stress are not as well understood. It is especially unclear which particular genes related to cell cycle control are involved in active growth repression. Here, we showed that decreased growth in plants exposed to moderate salt stress is mediated by MYB3R transcription factors that have been known to positively and negatively regulate the transcription of G2/M-specific genes. Our genome-wide gene expression analysis revealed occurrences of general downregulation of G2/M-specific genes in Arabidopsis under salt stress. Importantly, this downregulation is significantly and universally mitigated by the loss of MYB3R repressors by mutations. Accordingly, the growth performance of Arabidopsis plants under salt stress is significantly recovered in mutants lacking MYB3R repressors. This growth recovery involves improved cell proliferation that is possibly due to prolonging and accelerating cell proliferation, which were partly suggested by enlarged root meristem and increased number of cells positive for CYCB1;1-GUS. Our ploidy analysis further suggested that cell cycle progression at the G2 phase was delayed under salt stress, and this delay was recovered by loss of MYB3R repressors. Under salt stress, the changes in expression of MYB3R activators and repressors at both the mRNA and protein levels were not significant. This observation suggests novel mechanisms underlying MYB3R-mediated growth repression under salt stress that are different from the mechanisms operating under other stress conditions such as DNA damage and high temperature.

TRICHOME AND ARTEMISININ REGULATOR 2 positively regulates trichome development and artemisinin biosynthesis in Artemisia annua

Glandular secretory trichomes (GSTs) are regarded as biofactories for synthesizing, storing, and secreting artemisinin. It is necessary to figure out the initiation and development regulatory mechanism of GSTs to cultivate high-yielding Artemisia annua. Here, we identified an MYB transcription factor, AaTAR2, from bioinformatics analysis of the A. annua genome database and Arabidopsis trichome development-related genes. AaTAR2 is mainly expressed in young leaves and located in the nucleus. Repression and overexpression of AaTAR2 resulted in a decrease and increase, respectively, in the GSTs numbers, leaf biomass, and the artemisinin content in transgenic plants. Furthermore, the morphological characteristics changed obviously in trichomes, suggesting AaTAR2 plays a key role in trichome formation. In addition, the expression of flavonoid biosynthesis genes and total flavonoid content increased dramatically in AaTAR2-overexpressing transgenic plants. Owing to flavonoids possibly counteracting emerging resistance to artemisinin in Plasmodium species, AaTAR2 is a potential target to improve the effect of artemisinin in clinical therapy. Taken together, AaTAR2 positively regulates trichome development and artemisinin and flavonoid biosynthesis. A better understanding of this 'multiple functions' transcription factor may enable enhanced artemisinin and flavonoids yield. AaTAR2 is a potential breeding target for cultivating high-quality A. annua.

Jasmonate induced alternative splicing responses in Arabidopsis

Jasmonate is an essential phytohormone regulating plant growth, development, and defense. Alternative splicing (AS) in jasmonate ZIM-domain (JAZ) repressors is well-characterized and plays an important role in jasmonate signaling regulation. However, it is unknown whether other genes in the jasmonate signaling pathway are regulated by AS. We explore the potential for AS regulation in three Arabidopsis genotypes (WT, jaz2, jaz7) in response to methyl jasmonate (MeJA) treatment with respect to: (a) differential AS, (b) differential miRNA targeted AS, and (c) AS isoforms with novel functions. AS events identified from transcriptomic data were validated with proteomic data. Protein interaction networks identified two genes, SKIP and ALY4 whose products have both DNA- and RNA-binding affinities, as potential key regulators mediating jasmonate signaling and AS regulation. We observed cases where AS alone, or AS and transcriptional regulation together, can influence gene expression in response to MeJA. Twenty-one genes contain predicted miRNA target sites subjected to AS, which implies that AS is coupled to miRNA regulation. We identified 30 cases where alternatively spliced isoforms may have novel functions. For example, AS of bHLH160 generates an isoform without a basic domain, which may convert it from an activator to a repressor. Our study identified potential key regulators in AS regulation of jasmonate signaling pathway. These findings highlight the importance of AS regulation in the jasmonate signaling pathway, both alone and in collaboration with other regulators.

SIGNIFICANCE STATEMENT

By exploring alternative splicing, we demonstrate its regulation in the jasmonate signaling pathway alone or in collaboration with other posttranscriptional regulations such as nonsense and microRNA-mediated decay. A signal transduction network model for alternative splicing in jasmonate signaling pathway was generated, contributing to our understanding for this important, prevalent, but relatively unexplored regulatory mechanism in plants.

The alcohol dehydrogenase gene family in sugarcane and its involvement in cold stress regulation

BACKGROUND

Alcohol dehydrogenases (ADHs) in plants are encoded by a multigene family. ADHs participate in growth, development, and adaptation in many plant species, but the evolution and function of the ADH gene family in sugarcane is still unclear.

RESULTS

In the present study, 151 ADH genes from 17 species including 32 ADH genes in Saccharum spontaneum and 6 ADH genes in modern sugarcane cultivar R570 were identified. Phylogenetic analysis demonstrated two groups of ADH genes and suggested that these genes underwent duplication during angiosperm evolution. Whole-genome duplication (WGD)/segmental and dispersed duplications played critical roles in the expansion of ADH family in S. spontaneum and R570, respectively. ScADH3 was cloned and preferentially expressed in response to cold stress. ScADH3 conferred improved cold tolerance in E. coli cells. Ectopic expression showed that ScADH3 can also enhance cold tolerance in transgenic tobacco. The accumulation of reactive oxygen species (ROS) in leaves of transgenic tobacco was significantly lower than in wild-type tobacco. The transcript levels of ROS-related genes in transgenic tobacco increased significantly. ScADH3 seems to affect cold tolerance by regulating the ROS-related genes to maintain the ROS homeostasis.

CONCLUSIONS

This study depicted the size and composition of the ADH gene family in 17 species, and investigated their evolution pattern. Comparative genomics analysis among the ADH gene families of S. bicolor, R570 and S. spontaneum revealed their close evolutionary relationship. Functional analysis suggested that ScADH3, which maintained the steady state of ROS by regulating ROS-related genes, was related to cold tolerance. These findings will facilitate research on evolutionary and functional aspects of the ADH genes in sugarcane, especially for the understanding of ScADH3 under cold stress.

Phylogenomic Analysis of R2R3 MYB Transcription Factors in Sorghum and their Role in Conditioning Biofuel Syndrome

Background

Large scale cultivation of sorghum for food, feed, and biofuel requires concerted efforts for engineering multipurpose cultivars with optimised agronomic traits. Due to their vital role in regulating the biosynthesis of phenylpropanoid-derived compounds, biomass composition, biotic, and abiotic stress response, R2R3-MYB family transcription factors are ideal targets for improving environmental resilience and economic value of sorghum.

Methods

We used diverse computational biology tools to survey the sorghum genome to identify R2R3-MYB transcription factors followed by their structural and phylogenomic analysis. We used in-house generated as well as publicly available high throughput expression data to analyse the R2R3 expression patterns in various sorghum tissue types.

Results

We have identified a total of 134 R2R3-MYB genes from sorghum and developed a framework to predict gene functions. Collating information from the physical location, duplication, structural analysis, orthologous sequences, phylogeny, and expression patterns revealed the role of duplications in clade-wise expansion of the R2R3-MYB family as well as intra-clade functional diversification. Using publicly available and in-house generated RNA sequencing data, we provide MYB candidates for conditioning biofuel syndrome by engineering phenylpropanoid biosynthesis and sugar signalling pathways in sorghum.

Conclusion

The results presented here are pivotal to prioritize MYB genes for functional validation and optimize agronomic traits in sorghum.

Genome-Wide Identification and Comparative Analysis of MYB Transcription Factor Family in Musa acuminata and Musa balbisiana

MYB transcription factors (TFs) make up one of the most important TF families in plants. These proteins play crucial roles in processes related to development, metabolism, and stimulus-response; however, very few studies have been reported for the characterization of MYB TFs from banana. The current study identified 305 and 251 MYB genes from Musa acuminata and Musa balbisiana, respectively. Comprehensive details of MYBs are reported in terms of gene structure, protein domain, chromosomal localization, phylogeny, and expression patterns. Based on the exon-intron arrangement, these genes were classified into 12 gene models. Phylogenetic analysis of MYBs involving both species of banana, Oryza sativa, and Arabidopsis thaliana distributed these genes into 27 subfamilies. This highlighted not only the conservation, but also the gain/loss of MYBs in banana. Such genes are important candidates for future functional investigations. The MYB genes in both species exhibited a random distribution on chromosomes with variable densities. Estimation of gene duplication events revealed that segmental duplications represented the major factor behind MYB gene family expansion in banana. Expression profiles of MYB genes were also explored for their potential involvement in acetylene response or development. Collectively, the current comprehensive analysis of MYB genes in both species of banana will facilitate future functional studies.

Characterization of the Poplar R2R3-MYB Gene Family and Over-Expression of PsnMYB108 Confers Salt Tolerance in Transgenic Tobacco

The MYB, one of the largest transcription factor families in plants, is related to various biological processes. For an example, the R2R3-MYB family plays an important role in regulation of primary and secondary metabolism, plant growth and development, and responses to hormones and stresses. However, functional studies on the poplar R2R3-MYB genes are limited. In this study, we identified 207 poplar R2R3-MYB genes that are unevenly distributed on the 19 chromosomes of poplar, followed by characterization of their conserved domains. On the basis of phylogenetic analysis, these genes can be divided into 23 groups. Evidence from synteny analyses indicated that the poplar R2R3-MYB gene family is featured by tandem and segmental duplication events. On the basis of RNA-Seq data, we investigated salt responsive genes and explored their expression patterns. Furthermore, we cloned the PsnMYB108 gene from poplar, which is significantly up-regulated in roots and leaves in response to salt stress. To validate its function, we developed transgenic tobacco plants that over-express the PsnMYB108 gene. It appears that the transgenic lines are more tolerant to salt stress than the wild type does. Evidence from physiological analyses demonstrated that over-expression of PsnMYB108 may improve tobacco salt stress tolerance by increasing the reactive oxygen species scavenging ability and the accumulation of proline. These results laid the foundation for future analysis and functional studies of poplar R2R3-MYB family members, and revealed that PsnMYB108 plays an important role in improving plant salt stress tolerance.

Whole-Genome and Transposed Duplication Contributes to the Expansion and Diversification of TLC Genes in Maize

TLC (TRAM/LAG/CRN8) proteins play important roles in ceramide metabolism and mycotoxin resistance. Herein a comparative genomics analysis of TLCs was performed in 31 plant and 3 species from other kingdoms, with an emphasis mainly on maize. TLCs were conserved across kingdoms and expanded in angiosperms, largely due to whole-genome/segmental duplication (WGD/SD) under purifying selection. Phylogeny reconstruction by maximum-likelihood method uncovered five TLC clades, subsequently named as TRAM/LAG, CLN8, PS-TLC, TM136 and TLCD clades. Each clade of TLCs shared specific transmembrane regions and motif composition. Divisions of conserved motifs to subunits may have occurred in TM136-type TLCs. Focusing on maize, five WGD and two DNA-mediated transposed duplication (TD) pairs were discovered, accounting for 61.11% ZmTLCs. Combined with further expression analysis, significant divergence was found in expression patterns between most maize WGD pairs, indicating subfunctionalization or/and neofunctionalization. Moreover, ZmTLC5, a deduced parental copy in a TD pair, was highly induced under FB₁ and fungus pathogen injection and exhibited potential capacity to respond to environmental stimuli. Additionally, population genetics analysis showed that ZmTLC10 in the CLN8-clade may have experienced significant positive selection and differentiated between wild and inbred maize populations. Overall, our results help to decipher the evolutionary history of TLCs in maize and plants, facilitating further functional analysis of them.

Identification and characterization of GmMYB118 responses to drought and salt stress

BACKGROUND

Abiotic stress severely influences plant growth and development. MYB transcription factors (TFs), which compose one of the largest TF families, play an important role in abiotic stress responses.

RESULT

We identified 139 soybean MYB-related genes; these genes were divided into six groups based on their conserved domain and were distributed among 20 chromosomes (Chrs). Quantitative real-time PCR (qRT-PCR) indicated that GmMYB118 highly responsive to drought, salt and high temperature stress; thus, this gene was selected for further analysis. Subcellular localization revealed that the GmMYB118 protein located in the nucleus. Ectopic expression (EX) of GmMYB118 increased tolerance to drought and salt stress and regulated the expression of several stress-associated genes in transgenic Arabidopsis plants. Similarly, GmMYB118-overexpressing (OE) soybean plants generated via Agrobacterium rhizogenes (A. rhizogenes)-mediated transformation of the hairy roots showed improved drought and salt tolerance. Furthermore, compared with the control (CK) plants, the clustered, regularly interspaced, short palindromic repeat (CRISPR)-transformed plants exhibited reduced drought and salt tolerance. The contents of proline and chlorophyll in the OE plants were significantly greater than those in the CK plants, whose contents were greater than those in the CRISPR plants under drought and salt stress conditions. In contrast, the reactive oxygen species (ROS) and malondialdehyde (MDA) contents were significantly lower in the OE plants than in the CK plants, whose contents were lower than those in the CRISPR plants under stress conditions.

CONCLUSIONS

These results indicated that GmMYB118 could improve tolerance to drought and salt stress by promoting expression of stress-associated genes and regulating osmotic and oxidizing substances to maintain cell homeostasis.

TRANSPARENT TESTA GLABRA 1-Dependent Regulation of Flavonoid Biosynthesis

The flavonoid composition of various tissues throughout plant development is of biological relevance and particular interest for breeding. Arabidopsis thaliana TRANSPARENT TESTA GLABRA 1 (AtTTG1) is an essential regulator of late structural genes in flavonoid biosynthesis. Here, we provide a review of the regulation of the pathway's core enzymes through AtTTG1-containing R2R3-MYELOBLASTOSIS-basic HELIX-LOOP-HELIX-WD40 repeat (MBW(AtTTG1)) complexes embedded in an evolutionary context. We present a comprehensive collection of A. thalianattg1 mutants and AtTTG1 orthologs. A plethora of MBW(AtTTG1) mechanisms in regulating the five major TTG1-dependent traits is highlighted.

Divergence and Conservative Evolution of XTNX Genes in Land Plants

The Toll-interleukin-1 receptor (TIR) and Nucleotide-binding site (NBS) domains are two major components of the TIR-NBS-leucine-rich repeat family plant disease resistance genes. Extensive functional and evolutionary studies have been performed on these genes; however, the characterization of a small group of genes that are composed of atypical TIR and NBS domains, namely XTNX genes, is limited. The present study investigated this specific gene family by conducting genome-wide analyses of 59 green plant genomes. A total of 143 XTNX genes were identified in 51 of the 52 land plant genomes, whereas no XTNX gene was detected in any green algae genomes, which indicated that XTNX genes originated upon emergence of land plants. Phylogenetic analysis revealed that the ancestral XTNX gene underwent two rounds of ancient duplications in land plants, which resulted in the formation of clades I/II and clades IIa/IIb successively. Although clades I and IIb have evolved conservatively in angiosperms, the motif composition difference and sequence divergence at the amino acid level suggest that functional divergence may have occurred since the separation of the two clades. In contrast, several features of the clade IIa genes, including the absence in the majority of dicots, the long branches in the tree, the frequent loss of ancestral motifs, and the loss of expression in all detected tissues of Zea mays, all suggest that the genes in this lineage might have undergone pseudogenization. This study highlights that XTNX genes are a gene family originated anciently in land plants and underwent specific conservative pattern in evolution.