Insights into the Diversification and Evolution of R2R3-MYB Transcription Factors in Plants

Cloning, characterization and functional analysis of NtMYB306a gene reveals its role in wax alkane biosynthesis of tobacco trichomes and stress tolerance

Trichomes are specialized hair-like organs found on epidermal cells of many terrestrial plants, which protect plant from excessive transpiration and numerous abiotic and biotic stresses. However, the genetic basis and underlying mechanisms are largely unknown in Nicotiana tabacum (common tobacco), an established model system for genetic engineering and plant breeding. In present study, we identified, cloned and characterized an unknown function transcription factor NtMYB306a from tobacco cultivar K326 trichomes. Results obtained from sequence phylogenetic tree analysis showed that NtMYB306a-encoded protein belonged to S1 subgroup of the plants' R2R3-MYB transcription factors (TFs). Observation of the green fluorescent signals from NtMYB306a-GFP fusion protein construct exhibited that NtMYB306a was localized in nucleus. In yeast transactivation assays, the transformed yeast containing pGBKT7-NtMYB306a construct was able to grow on SD/-Trp-Ade+X-α-gal selection media, signifying that NtMYB306a exhibits transcriptional activation activity. Results from qRT-PCR, in-situ hybridization and GUS staining of transgenic tobacco plants revealed that NtMYB306a is primarily expressed in tobacco trichomes, especially tall glandular trichomes (TGTs) and short glandular trichomes (SGTs). RNA sequencing (RNA-seq) and qRT-PCR analysis of the NtMYB306a-overexpressing transgenic tobacco line revealed that NtMYB306a activated the expression of a set of key target genes which were associated with wax alkane biosynthesis. Gas Chromatography-Mass Spectrometry (GC-MS) exhibited that the total alkane contents and the contents of n-C28, n-C29, n-C31, and ai-C31 alkanes in leaf exudates of NtMYB306a-OE lines (OE-3, OE-13, and OE-20) were significantly greater when compared to WT. Besides, the promoter region of NtMYB306a contained numerous stress-responsive cis-acting elements, and their differential expression towards salicylic acid and cold stress treatments reflected their roles in signal transduction and cold-stress tolerance. Together, these results suggest that NtMYB306a is necessarily a positive regulator of alkane metabolism in tobacco trichomes that does not affect the number and morphology of tobacco trichomes, and that it can be used as a candidate gene for improving stress resistance and the quality of tobacco.

DhMYB22 and DhMYB60 regulate pigment intensity and floral organ shape in Dendrobium hybrid

BACKGROUND

Flower pigment and shape are determined by the coordinated expression of a set of structural genes during flower development. R2R3-MYB transcription factors are known regulators of structural gene expression. The current study focused on two members of this large family of transcription factors that were predicted to have roles in pigment biosynthesis and organ shape development in orchids.

METHODS

Phylogenetic analysis was used to identify candidate Dendrobium catenatum R2R3-MYB (DcaMYB) sequences associated with pigment and cell shape development. Gene silencing of candidate DhMYBs in Dendrobium hybrid by direct application of dsRNA to developing flowers was followed by observation of gene expression level and flower phenotypes. Silencing of the structural gene chalcone synthase was used as a comparative control.

KEY RESULTS

Ten candidate flower-associated DcaMYBs were identified. Flowers treated with dsRNA of DhMYB22 and DhMYB60 sequences were less pigmented and had relatively low expression of anthocyanin biosynthetic genes (F3'H and DFR), lower total anthocyanin concentration and markedly lower levels of cyanidin-3-glucoside and cyanidin-3-rutinoside. Petals of DhMYB22-treated flowers and sepals of DhMYB60-treated flowers showed the greatest colour difference relative to the same organs in untreated flowers. DhMYB22-treated flowers had relatively narrow and constricted lips, while DhMYB60-treated flowers had narrow and constricted sepals. No significant difference in shape was observed for DhCHS-treated or untreated flowers.

CONCLUSIONS

Our results demonstrate that DhMYB22 and DhMYB60 regulate pigment intensity and floral organ shape in Dendrobium. This is a first report of MYB regulation of floral organ shape in orchids.

A genome for Cissus illustrates features underlying its evolutionary success in dry savannas

Cissus is the largest genus in Vitaceae and is mainly distributed in the tropics and subtropics. Crassulacean acid metabolism (CAM), a photosynthetic adaptation to the occurrence of succulent leaves or stems, indicates that convergent evolution occurred in response to drought stress during species radiation. Here we provide the chromosomal level assembly of Cissus rotundifolia (an endemic species in Eastern Africa) and a genome-wide comparison with grape to understand genome divergence within an ancient eudicot family. Extensive transcriptome data were produced to illustrate the genetics underpinning C. rotundifolia's ecological adaption to seasonal aridity. The modern karyotype and smaller genome of C. rotundifolia (n = 12, 350.69 Mb/1C), which lack further whole-genome duplication, were mainly derived from gross chromosomal rearrangements such as fusions and segmental duplications, and were sculpted by a very recent burst of retrotransposon activity. Bias in local gene amplification contributed to its remarkable functional divergence from grape, and the specific proliferated genes associated with abiotic and biotic responses (e.g. HSP-20, NBS-LRR) enabled C. rotundifolia to survive in a hostile environment. Reorganization of existing enzymes of CAM characterized as diurnal expression patterns of relevant genes further confer the ability to thrive in dry savannas.

Regulation of Anthocyanin Biosynthesis by Drought and UV-B Radiation in Wild Tomato (Solanum peruvianum) Fruit

Anthocyanins are plant pigments derived from the phenylpropanoid pathway which are produced in many different species, contributing to defense against stresses by their antioxidant properties. Cultivated tomatoes cannot synthesize flavonoids; however, wild tomatoes such as Solanum chilense and Solanum lycopersicoides have anthocyanin pigmented skin. Other wild tomato species such as Solanum peruvianum have been poorly studied concerning anthocyanin accumulation in the fruit. This research is the first to address the regulation of anthocyanin biosynthesis mediated by drought stress and light radiation in S. peruvianum fruit. Transcript accumulation of SpAN2, encoding for a key MYB type transcription factor for the regulation of anthocyanin biosynthesis, was induced in the fruit of plants exposed to drought treatment. In addition, fruit peel accumulates a greater anthocyanin content in water deficit-treated plants. The expression of SpAN2 was also regulated according to sunlight exposure, reaching a higher expression during maximal daily UV radiation and under controlled UV-B treatments. Similar results were observed for the expression of the late flavonoid biosynthetic gene dihydroflavonol 4-reductase (SpDFR). These results suggest that SpAN2 and SpDFR are involved in anthocyanin biosynthesis under drought stress and UV radiation in S. peruvianum.

Analysis of flavonol regulator evolution in the Brassicaceae reveals MYB12, MYB111 and MYB21 duplications and MYB11 and MYB24 gene loss

BACKGROUND

Flavonols are the largest subgroup of flavonoids, possessing multiple functions in plants including protection against ultraviolet radiation, antimicrobial activities, and flower pigmentation together with anthocyanins. They are of agronomical and economical importance because the major off-taste component in rapeseed protein isolates is a flavonol derivative, which limits rapeseed protein use for human consumption. Flavonol production in Arabidopsis thaliana is mainly regulated by the subgroup 7 (SG7) R2R3-MYB transcription factors MYB11, MYB12, and MYB111. Recently, the SG19 MYBs MYB21, MYB24, and MYB57 were shown to regulate flavonol accumulation in pollen and stamens. The members of each subgroup are closely related, showing gene redundancy and tissue-specific expression in A. thaliana. However, the evolution of these flavonol regulators inside the Brassicaceae, especially inside the Brassiceae, which include the rapeseed crop species, is not fully understood.

RESULTS

We studied the SG7 and SG19 MYBs in 44 species, including 31 species of the Brassicaceae, by phylogenetic analyses followed by synteny and gene expression analyses. Thereby we identified a deep MYB12 and MYB111 duplication inside the Brassicaceae, which likely occurred before the divergence of Brassiceae and Thelypodieae. These duplications of SG7 members were followed by the loss of MYB11 after the divergence of Eruca vesicaria from the remaining Brassiceae species. Similarly, MYB21 experienced duplication before the emergence of the Brassiceae tribe, where the gene loss of MYB24 is also proposed to have happened. The members of each subgroup revealed frequent overlapping spatio-temporal expression patterns in the Brassiceae member B. napus, which are assumed to compensate for the loss of MYB11 and MYB24 in the analysed tissues.

CONCLUSIONS

We identified a duplication of MYB12, MYB111, and MYB21 inside the Brassicaceae and MYB11 and MYB24 gene loss inside the tribe Brassiceae. We propose that polyploidization events have shaped the evolution of the flavonol regulators in the Brassicaceae, especially in the Brassiceae.

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.

Characteristics and molecular identification of glyceraldehyde-3-phosphate dehydrogenases in poplar

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), an essential enzyme of the glycolysis metabolic pathway, plays a vital role in carbon metabolism, plant development, and stress resistance. As a kind of woody plant, poplars are widely cultivated for afforestation. Although the whole genome data of poplars have been published, little information is known about the GAPDH family of genes in poplar. This study performed a genome-wide identification of the poplar GAPDH family, and 13 determined PtGAPDH genes were identified from poplar genome. Phylogenetic tree showed that the PtGAPDH members were divided into PtGAPA/B, PtGAPC, PtGAPCp, and PtGAPN groups. A total of 13 PtGAPDH genes were distributed on eight chromosomes, 13 gene pairs belonging to segmented replication events were detected in poplar, and 23 collinearity gene pairs were determined between poplar and willow. The PtGAPDHcis-acting elements associated with growth and development as well as stress resistance revealed that PtGAPDHs might be involved in these processes. The phosphoglycerate kinase (PGK) and triose-phosphate isomerase (TPI) were predicted as the putative interaction proteins of PtGAPDHs. Gene ontology (GO) analysis showed that PtGAPDHs play a crucial role in the oxidation and reduction processes. PtGAPDH expression levels were induced by NaCl and PEG treatments, which implied that PtGAPDHs might be involved in stress response. Overexpression of PtGAPC1 significantly changed the contents of lipid and carbohydrate metabolites, which indicated that PtGAPC1 plays an essential role in metabolic regulation. This study highlights the characterizations and profiles of PtGAPDHs and reveals that PtGAPC1 is involved in the loop of lipid and carbohydrate metabolisms.

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.

Comparative transcriptomics of tropical woody plants supports fast and furious strategy along the leaf economics spectrum in lianas

Lianas, climbing woody plants, influence the structure and function of tropical forests. Climbing traits have evolved multiple times, including ancestral groups such as gymnosperms and pteridophytes, but the genetic basis of the liana strategy is largely unknown. Here, we use a comparative transcriptomic approach for 47 tropical plant species, including ten lianas of diverse taxonomic origins, to identify genes that are consistently expressed or downregulated only in lianas. Our comparative analysis of full-length transcripts enabled the identification of a core interactomic network common to lianas. Sets of transcripts identified from our analysis reveal features related to functional traits pertinent to leaf economics spectrum in lianas, include upregulation of genes controlling epidermal cuticular properties, cell wall remodeling, carbon concentrating mechanism, cell cycle progression, DNA repair and a large suit of downregulated transcription factors and enzymes involved in ABA-mediated stress response as well as lignin and suberin synthesis. All together, these genes are known to be significant in shaping plant morphologies through responses such as gravitropism, phyllotaxy and shade avoidance.

Summary: The full-length fraction of liana transcriptomes mapped on a protein–protein interactome revealed the nature of their convergence through distinct sets of expressed and downregulated genes not observed in free-standing plants.

The Transcription Factor MYB37 Positively Regulates Photosynthetic Inhibition and Oxidative Damage in Arabidopsis Leaves Under Salt Stress

MYB transcription factors (TFs) mediate plant responses and defenses to biotic and abiotic stresses. The effects of overexpression of MYB37, an R2R3 MYB subgroup 14 transcription factors in Arabidopsis thaliana, on chlorophyll content, chlorophyll fluorescence parameters, reactive oxygen species (ROS) metabolism, and the contents of osmotic regulatory substances were studied under 100 mM NaCl stress. Compared with the wild type (Col-0), MYB37 overexpression significantly alleviated the salt stress symptoms in A. thaliana plants. Chlorophyll a (Chl a) and chlorophyll b (Chl b) contents were significantly decreased in OE-1 and OE-2 than in Col-0. Particularly, the Chl a/b ratio was also higher in OE-1 and OE-2 than in Col-0 under NaCl stress. However, MYB37 overexpression alleviated the degradation of chlorophyll, especially Chl a. Salt stress inhibited the activities of PSII and PSI in Arabidopsis leaves, but did not affect the activity of PSII electron donor side oxygen-evolving complex (OEC). MYB37 overexpression increased photosynthesis in Arabidopsis by increasing PSII and PSI activities. MYB37 overexpression also promoted the transfer of electrons from Q _A to Q _B on the PSII receptor side of Arabidopsis under NaCl stress. Additionally, MYB37 overexpression increased Y(II) and Y(NPQ) of Arabidopsis under NaCl stress and decreased Y(NO). These results indicate that MYB37 overexpression increases PSII activity and regulates the proportion of energy dissipation in Arabidopsis leaves under NaCl stress, thus decreasing the proportion of inactivated reaction centers. Salt stress causes excess electrons and energy in the photosynthetic electron transport chain of Arabidopsis leaves, resulting in the release of reactive oxygen species (ROS), such as superoxide anion and hydrogen peroxide, leading to oxidative damage. Nevertheless, MYB37 overexpression reduced accumulation of malondialdehyde in Arabidopsis leaves under NaCl stress and alleviated the degree of membrane lipid peroxidation caused by ROS. Salt stress also enhanced the accumulation of soluble sugar (SS) and proline (Pro) in Arabidopsis leaves, thus reducing salt stress damage to plants. Salt stress also degraded soluble protein (SP). Furthermore, the accumulation of osmoregulation substances SS and Pro in OE-1 and OE-2 was not different from that in Col-0 since MYB37 overexpression in Arabidopsis OE-1, and OE-2 did not significantly affect plants under NaCl stress. However, SP content was significantly higher in OE-1 and OE-2 than in Col-0. These results indicate that MYB37 overexpression can alleviate the degradation of Arabidopsis proteins under NaCl stress, promote plant growth and improve salt tolerance.

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.

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.

MaMYB4 is a negative regulator and a substrate of RING-type E3 ligases MaBRG2/3 in controlling banana fruit ripening

Fruit ripening is a complex developmental process, which is modulated by both transcriptional and post-translational events. Control of fruit ripening is important in maintaining moderate quality traits and minimizing postharvest deterioration. In this study, we discovered that the transcription factor MaMYB4 acts as a negative regulator of fruit ripening in banana. The protein levels of MaMYB4 decreased gradually with banana fruit ripening, paralleling ethylene production, and decline in firmness. DNA affinity purification sequencing combined with RNA-sequencing analyses showed that MaMYB4 preferentially binds to the promoters of various ripening-associated genes including ethylene biosynthetic and cell wall modifying genes. Furthermore, ectopic expression of MaMYB4 in tomato delayed tomato fruit ripening, which was accompanied by downregulation of ethylene biosynthetic and cell wall modifying genes. Importantly, two RING finger E3 ligases MaBRG2/3, whose protein accumulation increased progressively with fruit ripening, were found to interact with and ubiquitinate MaMYB4, contributing to decreased accumulation of MaMYB4 during fruit ripening. Transient overexpression of MaMYB4 and MaBRG2/3 in banana fruit ripening delayed or promoted fruit ripening by inhibiting or stimulating ethylene biosynthesis, respectively. Taken together, we demonstrate that MaMYB4 negatively modulates banana fruit ripening, and that MaMYB4 abundance could be regulated by protein ubiquitination, thus providing insights into the role of MaMYB4 in controlling fruit ripening at both transcriptional and post-translational levels.

RLM1, Encoding an R2R3 MYB Transcription Factor, Regulates the Development of Secondary Cell Wall in Rice

Leaf morphology is an important component of rice ideal plant type. To date, many regulatory genes influencing leaf morphology in rice have been cloned, and their underlying molecular regulatory mechanism has been preliminarily clarified. However, the fine regulation relationship of leaf morphogenesis and plant type remains largely elusive. In this study, a rolling-leaf mutant, named rlm1-D, was obtained and controlled by a pair of dominant nuclear genes. Cytological observations revealed that the rlm1 was mainly caused by abnormal deposition of secondary cell walls. Molecular evidence showed ectopic expression of a MYB-type transcription factor LOC_Os05g46610 was responsible for the phenotype of rlm1-D. A series of experiments, including the transcription factor-centered technology, DNA-binding assay, and electrophoretic mobility shift assay, verified that RLM1 can bind to the promoter of OsCAD2, a key gene responsible for lignin biosynthesis in rice. An interacting partner of RLM1, OsMAPK10, was identified. Multiple biochemical assays confirmed that OsMAPK10 interacted with RLM1. OsMAPK10 positively regulated the lignin content in the leaves and stems of rice. Moreover, OsMAPK10 contributes to RLM1 activation of downstream target genes. In particular, RLM1 is exclusively expressed in the stems at the mature plant stage. The yield of RLM1 knockdown lines increased by over 11% without other adverse agricultural trait penalties, indicating great practical application value. A MAPK-MYB-OsCAD2 genetic regulatory network controlling SCW was proposed, providing a theoretical significance and practical value for shaping the ideal plant type and improving rice yield.

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.

Transcription factors evolve faster than their structural gene targets in the flavonoid pigment pathway

Dissecting the relationship between gene function and substitution rates is key to understanding genome-wide patterns of molecular evolution. Biochemical pathways provide powerful systems for investigating this relationship because the functional role of each gene is often well characterized. Here, we investigate the evolution of the flavonoid pigment pathway in the colorful Petunieae clade of the tomato family (Solanaceae). This pathway is broadly conserved in plants, both in terms of its structural elements and its MYB, basic helix–loop–helix, and WD40 transcriptional regulators, and its function has been extensively studied, particularly in model species of petunia. We built a phylotranscriptomic data set for 69 species of Petunieae to infer patterns of molecular evolution across pathway genes and across lineages. We found that transcription factors exhibit faster rates of molecular evolution (dN/dS) than their targets, with the highly specialized MYB genes evolving fastest. Using the largest comparative data set to date, we recovered little support for the hypothesis that upstream enzymes evolve slower than those occupying more downstream positions, although expression levels do predict molecular evolutionary rates. Although shifts in floral pigmentation were only weakly related to changes affecting coding regions, we found a strong relationship with the presence/absence patterns of MYB transcripts. Intensely pigmented species express all three main MYB anthocyanin activators in petals, whereas pale or white species express few or none. Our findings reinforce the notion that pathway regulators have a dynamic history, involving higher rates of molecular evolution than structural components, along with frequent changes in expression during color transitions.

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

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

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.

Origin and adaptive evolution of UV RESISTANCE LOCUS 8-mediated signaling during plant terrestrialization

UV RESISTANCE LOCUS 8 (UVR8) mediates photomorphogenic responses and acclimation to UV-B radiation by regulating the transcription of a series of transcription factors (TFs). However, the origin and evolution of UVR8-mediated signaling pathways remain largely unknown. In this study, we investigated the origin and evolution of the major components of the UVR8-mediated signaling pathway (UVR8, REPRESSOR OF UV-B PHOTOMORPHOGENESIS [RUP], BRI1-EMS-SUPPRESSOR1 [BES1], BES1-INTERACTING MYC-LIKE 1 (BIM1), WRKY DNA-BINDING PROTEIN 36 (WRKY36), MYB DOMAIN PROTEIN 73/77/13 [MYB73/MYB77/MYB13], and PHYTOCHROME INTERACTING FACTOR 4/5 [PIF4 and PIF5]) using comparative genomics and phylogenetic approaches. We showed that the central regulator UVR8 presented a conservative evolutionary route during plant evolution, and the evolutionary history of downstream negative regulators and TFs was different from that of green plant phylogeny. The canonical UVR8-CONSTITUTIVELY PHOTOMORPHOGENIC 1(COP1)/SUPPRESSOR OF PHYA-105 (SPA)-ELONGATED HYPOCOTYL 5 (HY5)-RUP signaling pathway originated in chlorophytes and conferred green algae the additional ability to cope with UV-B radiation. Moreover, the emergence of multiple UVR8-mediated signaling pathways in charophytes laid the foundations for the cross-talk between UV-B signals and endogenous hormone responses. Importantly, we observed signatures that reflect plant adaptations to high UV-B irradiance in subaerial/terrestrial environments, including positive selection in UVR8 and RUPs and increased copy number of some vital TFs. These results revealed that green plants not only experienced adaptive modifications in the canonical UVR8-COP1/SPA-HY5-RUP signaling pathway, but also diversified their UV-B signal transduction mechanisms through increasing cross-talk with other pathways, such as those associated with brassinosteroids and auxin. This study greatly expands our understanding of molecular evolution and adaptive mechanisms underlying plant UV-B acclimation.

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.

Genome-wide identification of myeloblastosis gene family and its response to cadmium stress in Ipomoea aquatica

The myeloblastosis (MYB) proteins perform key functions in mediating cadmium (Cd) tolerance of plants. Ipomoea aquatica has strong adaptability to Cd Stress, while the roles of the I. aquatica MYB gene family with respect to Cd stress are still unclear. Here, we identified a total of 183 MYB genes in the I. aquatica genome (laMYB), which were classified into 66 1R-type IaMYB, 112 2R-type IaMYB, four 3R-type IaMYB, and one 4R-type IaMYB based on the number of the MYB repeat in each gene. The analysis of phylogenetic tree indicated that most of IaMYB genes are associated with the diverse biological processes including defense, development and metabolism. Analysis of sequence features showed that the IaMYB genes within identical subfamily have the similar patterns of the motif distributions and gene structures. Analysis of gene duplication events revealed that the dispersed duplication (DSD) and whole-genome duplication (WGD) modes play vital roles in the expansion of the IaMYB gene family. Expression profiling manifests that approximately 20% of IaMYB genes had significant role in the roots of I. aquatica under Cd stress. Promoter profiling implied that the differentially expressed genes might be induced by environmental factors or inherent hormones and thereby execute their function in Cd response. Remarkably, the 2R-type IaMYB157 with abundant light-responsive element G-box and ABA-responsive element ABRE in its promoter region exhibited very strong response to Cd stress. Taken together, our findings provide an important candidate IaMYB gene for further deciphering the molecular regulatory mechanism in plant with respect to Cd stress.

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.

A cellular mechanism underlying the restoration of thermo/photoperiod-sensitive genic male sterility

During anther development, the transformation of the microspore into mature pollen occurs under the protection of first the tetrad wall and later the pollen wall. Mutations in genes involved in this wall transition often lead to microspore rupture and male sterility; some such mutants, such as the reversible male sterile (rvms) mutant, are thermo/photoperiod-sensitive genic male sterile (P/TGMS) lines. Previous studies have shown that slow development is a general mechanism of P/TGMS fertility restoration. In this study, we identified restorer of rvms-2 (res2), which is an allele of QUARTET 3 (QRT3) encoding a polygalacturonase that shows delayed degradation of the tetrad pectin wall. We found that MS188, a tapetum-specific transcription factor essential for pollen wall formation, can activate QRT3 expression for pectin wall degradation, indicating a non-cell-autonomous pathway involved in the regulation of the cell wall transition. Further assays showed that a delay in degradation of the tetrad pectin wall is responsible for the fertility restoration of rvms and other P/TGMS lines, whereas early expression of QRT3 eliminates low temperature restoration of rvms-2 fertility. Taken together, these results suggest a likely cellular mechanism of fertility restoration in P/TGMS lines, that is, slow development during the cell wall transition of P/TGMS microspores may reduce the requirement for their wall protection and thus support their development into functional pollens, leading to restored fertility.

Transcriptomic and Proteomic Profiling Reveal the Key Role of AcMYB16 in the Response of Pseudomonas syringae pv. actinidiae in Kiwifruit

Kiwifruit bacterial canker caused by Pseudomonas syringae pv. actinidiae (Psa), is an important disease of kiwifruit (Actinidia Lind.). Plant hormones may induce various secondary metabolites to resist pathogens via modulation of hormone-responsive transcription factors (TFs), as reported in past studies. In this study, we showed that JA accumulated in the susceptible cultivar Actinidia chinensis 'Hongyang' but decreased in the resistant cultivar of A. chinensis var. deliciosa 'Jinkui' in response to Psa. Integrated transcriptomic and proteomic analyses were carried out using the resistant cultivar 'Jinkui'. A total of 5,045 differentially expressed genes (DEGs) and 1,681 differentially expressed proteins (DEPs) were identified after Psa infection. Two pathways, 'plant hormone signal transduction' and 'phenylpropanoid biosynthesis,' were activated at the protein and transcript levels. In addition, a total of 27 R2R3-MYB transcription factors (TFs) were involved in the response to Psa of 'Jinkui,' including the R2R3-MYB TF subgroup 4 gene AcMYB16, which was downregulated in 'Jinkui' but upregulated in 'Hongyang.' The promoter region of AcMYB16 has a MeJA responsiveness cis-acting regulatory element (CRE). Transient expression of the AcMYB16 gene in the leaves of 'Jinkui' induced Psa infection. Together, these data suggest that AcMYB16 acts as a repressor to regulate the response of kiwifruit to Psa infection. Our work will help to unravel the processes of kiwifruit resistance to pathogens and will facilitate the development of varieties with resistance against bacterial pathogens.

Transcriptome Analysis Revealed the Mechanism by Which Exogenous ABA Increases Anthocyanins in Blueberry Fruit During Veraison

Blueberry (Vaccinium spp.) is a popular healthy fruit worldwide. The health value of blueberry is mainly because the fruit is rich in anthocyanins, which have a strong antioxidant capacity. However, because blueberry is a non-model plant, little is known about the structural and regulatory genes involved in anthocyanin synthesis in blueberries. Previous studies have found that spraying 1,000 mg/L abscisic acid at the late green stage of "Jersey" highbush blueberry fruits can increase the content of anthocyanins. In this experiment, the previous results were verified in "Brightwell" rabbiteye blueberry fruits. Based on the previous results, the anthocyanin accumulation process in blueberry can be divided into six stages from the late green stage to the mature stage, and the transcriptome was used to systematically analyze the blueberry anthocyanin synthesis process. Combined with data from previous studies on important transcription factors regulating anthocyanin synthesis in plants, phylogenetic trees were constructed to explore the key transcription factors during blueberry fruit ripening. The results showed that ABA increased the anthocyanin content of blueberry fruits during veraison. All structural genes and transcription factors (MYB, bHLH, and WD40) involved in the anthocyanin pathway were identified, and their spatiotemporal expression patterns were analyzed. The expression of CHS, CHI, DFR, and LDOX/ANS in ABA-treated fruits was higher in the last two stages of maturity, which was consistent with the change in the anthocyanin contents in fruits. In general, six MYB transcription factors, one bHLH transcription factor and four WD40 transcription factors were found to change significantly under treatment during fruit ripening. Among them, VcMYBA plays a major role in the regulation of anthocyanin synthesis in ABA signaling. This result preliminarily explained the mechanism by which ABA increases the anthocyanin content and improves the efficiency of the industrial use of blueberry anthocyanins.

MYBA and MYBPA transcription factors co-regulate anthocyanin biosynthesis in blue-coloured berries

The regulatory network of R2R3 MYB transcription factors in anthocyanin biosynthesis is not fully understood in blue-coloured berries containing delphinidin compounds. We used blue berries of bilberry (Vaccinium myrtillus) to comprehensively characterise flavonoid-regulating R2R3 MYBs, which revealed a new type of co-regulation in anthocyanin biosynthesis between members of MYBA-, MYBPA1- and MYBPA2-subgroups. VmMYBA1, VmMYBPA1.1 and VmMYBPA2.2 expression was elevated at berry ripening and by abscisic acid treatment. Additionally, VmMYBA1 and VmMYBPA1.1 expression was strongly downregulated in a white berry mutant. Complementation and transient overexpression assays confirmed VmMYBA1 and VmMYBA2 to induce anthocyanin accumulation. Promoter activation assays showed that VmMYBA1, VmMYBPA1.1 and VmMYBPA2.2 had similar activity towards dihydroflavonol 4-reductase (DFR) and anthocyanidin synthase (ANS), but differential regulation activity for UDP-glucose flavonoid 3-O-glucosyltransferase (UFGT) and flavonoid 3'5'-hydroxylase (F3'5'H) promoters. Silencing of VmMYBPA1.1 in berries led to the downregulation of key anthocyanin and delphinidin biosynthesis genes. Functional analyses of other MYBPA regulators, and a member of novel MYBPA3 subgroup, associated them with proanthocyanidin biosynthesis and F3'5'H expression. The existence of 18 flavonoid-regulating MYBs indicated gene duplication, which may have enabled functional diversification among MYBA, MYBPA1 and MYBPA2 subgroups. Our results provide new insights into the intricate regulation of the complex anthocyanin profile found in blue-coloured berries involving regulation of both cyanidin and delphinidin branches.

The Streptochaeta Genome and the Evolution of the Grasses

In this work, we sequenced and annotated the genome of Streptochaeta angustifolia, one of two genera in the grass subfamily Anomochlooideae, a lineage sister to all other grasses. The final assembly size is over 99% of the estimated genome size. We find good collinearity with the rice genome and have captured most of the gene space. Streptochaeta is similar to other grasses in the structure of its fruit (a caryopsis or grain) but has peculiar flowers and inflorescences that are distinct from those in the outgroups and in other grasses. To provide tools for investigations of floral structure, we analyzed two large families of transcription factors, AP2-like and R2R3 MYBs, that are known to control floral and spikelet development in rice and maize among other grasses. Many of these are also regulated by small RNAs. Structure of the gene trees showed that the well documented whole genome duplication at the origin of the grasses (ρ) occurred before the divergence of the Anomochlooideae lineage from the lineage leading to the rest of the grasses (the spikelet clade) and thus that the common ancestor of all grasses probably had two copies of the developmental genes. However, Streptochaeta (and by inference other members of Anomochlooideae) has lost one copy of many genes. The peculiar floral morphology of Streptochaeta may thus have derived from an ancestral plant that was morphologically similar to the spikelet-bearing grasses. We further identify 114 loci producing microRNAs and 89 loci generating phased, secondary siRNAs, classes of small RNAs known to be influential in transcriptional and post-transcriptional regulation of several plant functions.

Concerted control of the LaRAV1-LaCDKB1;3 module by temperature during dormancy release and reactivation of larch

Dormancy release and reactivation of temperate-zone trees involve the temperature-modulated expression of cell-cycle genes. However, information on the detailed regulatory mechanism is limited. Here, we compared the transcriptomes of the stems of active and dormant larch trees, emphasizing the expression patterns of cell-cycle genes and transcription factors and assessed their relationships and responses to temperatures. Twelve cell-cycle genes and 31 transcription factors were strongly expressed in the active stage. Promoter analysis suggested that these 12 genes might be regulated by transcription factors from 10 families. Altogether, 73 cases of regulation between 16 transcription factors and 12 cell-cycle genes were predicted, while the regulatory interactions between LaMYB20 and LaCYCB1;1, and LaRAV1 and LaCDKB1;3 were confirmed by yeast one-hybrid and dual-luciferase assays. Last, we found that LaRAV1 and LaCDKB1;3 had almost the same expression patterns during dormancy release and reactivation induced naturally or artificially by temperature, indicating that the LaRAV1-LaCDKB1;3 module functions in the temperature-modulated dormancy release and reactivation of larch trees. These results provide new insights into the link between temperature and cell-cycle gene expression, helping to understand the temperature control of tree growth and development in the context of climate change.

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.

The Welwitschia genome reveals a unique biology underpinning extreme longevity in deserts

The gymnosperm Welwitschia mirabilis belongs to the ancient, enigmatic gnetophyte lineage. It is a unique desert plant with extreme longevity and two ever-elongating leaves. We present a chromosome-level assembly of its genome (6.8 Gb/1 C) together with methylome and transcriptome data to explore its astonishing biology. We also present a refined, high-quality assembly of Gnetum montanum to enhance our understanding of gnetophyte genome evolution. The Welwitschia genome has been shaped by a lineage-specific ancient, whole genome duplication (~86 million years ago) and more recently (1-2 million years) by bursts of retrotransposon activity. High levels of cytosine methylation (particularly at CHH motifs) are associated with retrotransposons, whilst long-term deamination has resulted in an exceptionally GC-poor genome. Changes in copy number and/or expression of gene families and transcription factors (e.g. R2R3MYB, SAUR) controlling cell growth, differentiation and metabolism underpin the plant's longevity and tolerance to temperature, nutrient and water stress.

R2R3-MYB transcription factor family in tea plant (Camellia sinensis): Genome-wide characterization, phylogeny, chromosome location, structure and expression patterns

MYB transcription factors play essential roles in many biological processes and environmental stimuli. However, the functions of the MYB transcription factor family in tea plants have not been elucidated. Here, a total of 122 CsR2R3-MYB genes were identified from the chromosome level genome of tea plant (Camellia sinensis). The CsR2R3-MYB genes were phylogenetically classified into 25 groups. Results from the structure analysis of the gene, conserved motifs, and chromosomal distribution supported the relative conservation of the R2R3-MYB genes family in the tea plant. Synteny analysis indicated that 122, 34, and 112 CsR2R3-MYB genes were orthologous to Arabidopsis thaliana, Oryza sativa and C. sinensis var. 'huangdan' (HD), respectively. Tissue-specific expression showed that all CsR2R3-MYB genes had different expression patterns in the tea plant tissues, indicating that these genes may perform diverse functions. The expression patterns of representative R2R3-MYB genes and the regulatory network of the main anthocyanin components were analyzed, which suggested that CsMYB17 may played a key role in the regulation of cya-3-O-gal, del-3-O-gal, cya-3-O-glu and pel-3-O-glu. Results from the qRT-PCR validation of selected genes suggested that CsR2R3-MYB genes were induced in response to drought, cold, GA, and ABA treatments. Overall, this study provides comprehensive and systematic information for research on the function of R2R3-MYB genes in tea plants.

Genome-Wide Identification and Analysis of the MYB Transcription Factor Gene Family in Chili Pepper (Capsicum spp.)

The MYB transcription factor family is very large and functionally diverse in plants, however, only a few members of this family have been reported and characterized in chili pepper (Capsicum spp.). In the present study, we performed genome-wide analyses of the MYB family in Capsicum annuum, including phylogenetic relationships, conserved domain, gene structure organization, motif protein arrangement, chromosome distribution, chemical properties predictions, RNA-seq expression, and RT-qPCR expression assays. A total of 235 non-redundant MYB proteins were identified from C. annuum, including R2R3-MYB, 3R-MYB, atypical MYB, and MYB-related subclasses. The sequence analysis of CaMYBs compared with other plant MYB proteins revealed gene conservation, but also potential specialized genes. Tissue-specific expression profiles showed that CaMYB genes were differentially expressed, suggesting that they are functionally divergent. Furthermore, the integration of our data allowed us to propose strong CaMYBs candidates to be regulating phenylpropanoid, lignin, capsaicinoid, carotenoid, and vitamin C biosynthesis, providing new insights into the role of MYB transcription factors in secondary metabolism. This study adds valuable knowledge about the functions of CaMYB genes in various processes in the Capsicum genus.

The land plant-specific MIXTA-MYB lineage is implicated in the early evolution of the plant cuticle and the colonization of land

The evolution of a lipid-based cuticle on aerial plant surfaces that protects against dehydration is considered a fundamental innovation in the colonization of the land by the green plants. However, key evolutionary steps in the early regulation of cuticle synthesis are still poorly understood, owing to limited studies in early-diverging land plant lineages. Here, we characterize a land plant specific subgroup 9 R2R3 MYB transcription factor MpSBG9, in the early-diverging land plant model Marchantia polymorpha, that is homologous to MIXTA proteins in vascular plants. The MpSBG9 functions as a key regulator of cuticle biosynthesis by preferentially regulating expression of orthologous genes for cutin formation, but not wax biosynthesis genes. The MpSBG9 also promotes the formation of papillate cells on the adaxial surface of M. polymorpha, which is consisitent with its canonical role in vascular plants. Our observations imply conserved MYB transcriptional regulation in the control of the cutin biosynthesis pathway as a core genetic network in the common ancestor of all land plants, implicating the land plant-specific MIXTA MYB lineage in the early origin and evolution of the cuticle.

Genome-Wide Mining of MYB Transcription Factors in the Anthocyanin Biosynthesis Pathway of Gossypium Hirsutum

The MYB family, one of the largest transcription factor (TF) families, plays an important role in plant growth, development, and stress response. Although genome-wide analysis of the MYB family has been performed in many species based on sequence similarity, predicting the potential functions of the MYB genes and classifying the regulators into specific metabolic pathways remains difficult. In this study, using a hidden Markov model search and co-expression regulatory network analysis, we demonstrated a process to screen and identify potential MYB TFs in the anthocyanin biosynthesis pathway of Gossypium hirsutum. As a result, we identified 617 and 784 MYB genes (812 in total) from the previously reported and recently released genomes, respectively. Using 126 structural genes involved in the anthocyanin biosynthesis pathway as targets for several co-expression network analyses, we sorted out 31 R2R3-MYB genes, which are potential regulators in the specific pathway. Phylogenetic and collinearity analyses indicated that 83.9% of the 31 MYB genes originated from whole genome duplication or polyploidization. In addition, we revealed relatively specific regulatory relationships between the MYB TFs and their target structural genes. Approximately, 71% of the MYBs could regulate only a single anthocyanin-related structural gene. Moreover, we found that the A- and D- subgenome homoeologs of MYB TFs in G. hirsutum rarely co-regulate the same target gene. The current study not only demonstrated an easy method to rapidly predict potential TFs in a specific metabolic pathway, but also enhanced our understanding of the evolution, gene characteristics, expression, and regulatory pattern of MYB TFs in G. hirsutum.

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.

The Expression of IbMYB1 Is Essential to Maintain the Purple Color of Leaf and Storage Root in Sweet Potato [Ipomoea batatas (L.) Lam]

IbMYB1 was one of the major anthocyanin biosynthesis regulatory genes that has been identified and utilized in purple-fleshed sweet potato breeding. At least three members of this gene, namely, IbMYB1-1, -2a, and -2b, have been reported. We found that IbMYB1-2a and -2b are not necessary for anthocyanin accumulation in a variety of cultivated species (hexaploid) with purple shoots or purplish rings/spots of flesh. Transcriptomic and quantitative reverse transcription PCR (RT-qPCR) analyses revealed that persistent and vigorous expression of IbMYB1 is essential to maintain the purple color of leaves and storage roots in this type of cultivated species, which did not contain IbMYB1-2 gene members. Compared with IbbHLH2, IbMYB1 is an early response gene of anthocyanin biosynthesis in sweet potato. It cannot exclude the possibility that other MYBs participate in this gene regulation networks. Twenty-two MYB-like genes were identified from 156 MYBs to be highly positively or negatively correlated with the anthocyanin content in leaves or flesh. Even so, the IbMYB1 was most coordinately expressed with anthocyanin biosynthesis genes. Differences in flanking and coding sequences confirm that IbMYB2s, the highest similarity genes of IbMYB1, are not the members of IbMYB1. This phenomenon indicates that there may be more members of IbMYB1 in sweet potato, and the genetic complementation of these members is involved in the regulation of anthocyanin biosynthesis. The 3' flanking sequence of IbMYB1-1 is homologous to the retrotransposon sequence of TNT1-94. Transposon movement is involved in the formation of multiple members of IbMYB1. This study provides critical insights into the expression patterns of IbMYB1, which are involved in the regulation of anthocyanin biosynthesis in the leaf and storage root. Notably, our study also emphasized the presence of a multiple member of IbMYB1 for genetic improvement.

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.

Azolla ferns testify: seed plants and ferns share a common ancestor for leucoanthocyanidin reductase enzymes

Questions about in vivo substrates for proanthocyanidin (PA) biosynthesis and condensation have not been resolved and wide gaps in the understanding of transport and biogenesis in 'tannosomes' persist. Here we examined the evolution of PA biosynthesis in ferns not previously reported, asking what PAs are synthesised and how. Chemical and gene-expression analyses were combined to characterise PA biosynthesis, leveraging genome annotation from the floating fern Azolla filiculoides. In vitro assay and phylogenomics of PIP-dehydrogenases served to infer the evolution of leucoanthocyanidin reductase (LAR). Sporophyte-synthesised (epi)catechin polymers, averaging only seven subunits, accumulated to 5.3% in A. filiculoides, and 8% in A. pinnata biomass dry weight. Consistently, a LAR active in vitro was highly expressed in A. filiculoides. LAR, and paralogous fern WLAR-enzymes with differing substrate binding sites, represent an evolutionary innovation of the common ancestor of fern and seed plants. The specific ecological niche of Azolla ferns, a floating plant-microbe mat massively fixing CO₂ and N₂ , shaped their metabolism in which PA biosynthesis predominates and employs novel fern LAR enzymes. Characterisation of in vivo substrates of these LAR, will help to shed light on the recently assigned and surprising dual catalysis of LAR from seed plants.

Far-Red Light-Induced Azolla filiculoides Symbiosis Sexual Reproduction: Responsive Transcripts of Symbiont Nostoc azollae Encode Transporters Whilst Those of the Fern Relate to the Angiosperm Floral Transition

Water ferns of the genus Azolla and the filamentous cyanobacteria Nostoc azollae constitute a model symbiosis that enabled the colonization of the water surface with traits highly desirable for the development of more sustainable crops: their floating mats capture CO₂ and fix N₂ at high rates using light energy. Their mode of sexual reproduction is heterosporous. The regulation of the transition from the vegetative phase to the spore forming phase in ferns is largely unknown, yet a prerequisite for Azolla domestication, and of particular interest as ferns represent the sister lineage of seed plants. Sporocarps induced with far red light could be crossed so as to verify species attribution of strains from the Netherlands but not of the strain from the Anzali lagoon in Iran; the latter strain was assigned to a novel species cluster from South America. Red-dominated light suppresses the formation of dissemination stages in both gametophyte- and sporophyte-dominated lineages of plants, the response likely is a convergent ecological strategy to open fields. FR-responsive transcripts included those from MIKC^C homologues of CMADS1 and miR319-controlled GAMYB transcription factors in the fern, transporters in N. azollae, and ycf2 in chloroplasts. Loci of conserved microRNA (miRNA) in the fern lineage included miR172, yet FR only induced miR529 and miR535, and reduced miR319 and miR159. Phylogenomic analyses of MIKC^C TFs suggested that the control of flowering and flower organ specification may have originated from the diploid to haploid phase transition in the homosporous common ancestor of ferns and seed plants.

Origins and Stepwise Expansion of R2R3-MYB Transcription Factors for the Terrestrial Adaptation of Plants

The R2R3-MYB transcription factors play critical roles in various processes in embryophytes (land plants). Here, we identified genes encoding R2R3-MYB proteins from rhodophytes, glaucophytes, Chromista, chlorophytes, charophytes, and embryophytes. We classified the R2R3-MYB genes into three subgroups (I, II, and III) based on their evolutionary history and gene structure. The subgroup I is the most ancient group that includes members from all plant lineages. The subgroup II was formed before the divergence of charophytes and embryophytes. The subgroup III genes form a monophyletic group and only comprise members from land plants with conserved exon-intron structure. Each subgroup was further divided into multiple clades. The subgroup I can be divided into I-A, I-B, I-C, and I-D. The I-A, I-B, and I-C are the most basal clades that have originated before the divergence of Archaeplastida. The I-D with the II and III subgroups form a monophyletic group, containing only green plants. The II and III subgroups form another monophyletic group with Streptophyta only. Once on land, the subgroup III genes have experienced two rounds of major expansions. The first round occurred before the origin of land plants, and the second round occurred after the divergence of land plants. Due to significant gene expansion, the subgroup III genes have become the predominant group of R2R3-MYBs in land plants. The highly unbalanced pattern of birth and death evolution of R2R3-MYB genes indicates their important roles in the successful adaptation and massive radiation of land plants to occupy a multitude of terrestrial environments.