Expansion and diversification of the Populus R2R3-MYB family of transcription factors

RNA-Seq Analysis of Prickled and Prickle-Free Epidermis Provides Insight into the Genetics of Prickle Development in Red Raspberry (Rubus ideaus L.)

Red raspberry (Rubus idaeus L.) is a globally commercialized specialty crop with growing demand worldwide. The presence of prickles on the stems, petioles and undersides of the leaves complicates both the field management and harvesting of raspberries. An RNA sequencing analysis was used to identify differentially expressed genes in the epidermal tissue of prickled “Caroline” and prickle-free “Joan J.” and their segregating progeny. Expression patterns of differentially expressed genes (DEGs) in prickle-free plants revealed the downregulation of some vital development-related transcription factors (TFs), including a MIXTA-like R2R3-MYB family member; MADS-box; APETALA2/ETHYLENE RESPONSIVE FACTOR (AP2/ERF) and NAM, ATAF1/2 and CUC2 (NAC) in prickle-free epidermis tissue. The downregulation of these TFs was confirmed by qRT-PCR analysis, indicating a key regulatory role in prickle development. This study adds to the understanding of prickle development mechanisms in red raspberries needed for utilizing genetic engineering strategies for developing prickle-free raspberry cultivars and, possibly, other Rubus species, such as blackberry (Rubus sp.) and black raspberry (R. occidentalis L.).

In Silico Identification of MYB and bHLH Families Reveals Candidate Transcription Factors for Secondary Metabolic Pathways in Cannabis sativa L

Plant secondary metabolic pathways are finely regulated by the activity of transcription factors, among which members of the bHLH and MYB subfamilies play a main role. Cannabis sativa L. is a unique officinal plant species with over 600 synthesized phytochemicals having diverse scale-up industrial and pharmaceutical usage. Despite comprehensive knowledge of cannabinoids' metabolic pathways, very little is known about their regulation, while the literature on flavonoids' metabolic pathways is still scarce. In this study, we provide the first genome-wide analysis of bHLH and MYB families in C. sativa reference cultivar CBDRx and identification of candidate coding sequences for these transcription factors. Cannabis sativa bHLHs and MYBs were then classified into functional subfamilies through comparative phylogenetic analysis with A. thaliana transcription factors. Analyses of gene structure and motif distribution confirmed that CsbHLHs and CsMYBs belonging to the same evolutionary clade share common features at both gene and amino acidic level. Candidate regulatory genes for key metabolic pathways leading to flavonoid and cannabinoid synthesis in Cannabis were also retrieved. Furthermore, a candidate gene approach was used to identify structural enzyme-coding genes for flavonoid and cannabinoid synthesis. Taken as a whole, this work represents a valuable resource of candidate genes for further investigation of the C. sativa cannabinoid and flavonoid metabolic pathways for genomic studies and breeding programs.

Structure and expression analysis of seven salt-related ERF genes of Populus

Ethylene response factors (ERFs) are plant-specific transcription factors (TFs) that play important roles in plant growth and stress defense and have received a great amount of attention in recent years. In this study, seven ERF genes related to abiotic stress tolerance and response were identified in plants of the Populus genus. Systematic bioinformatics, including sequence phylogeny, genome organisation, gene structure, gene ontology (GO) annotation, etc. were detected. Expression-pattern of these seven ERF genes were analyzed using RT-qPCR and cross validated using RNA-Seq. Data from a phylogenetic tree and multiple alignment of protein sequences indicated that these seven ERF TFs belong to three subfamilies and contain AP2, YRG, and RAYD conserved domains, which may interact with downstream target genes to regulate the plant stress response. An analysis of the structure and promoter region of these seven ERF genes showed that they have multiple stress-related motifs and cis-elements, which may play roles in the plant stress-tolerance process through a transcriptional regulation mechanism; moreover, the cellular_component and molecular_function terms associated with these ERFs determined by GO annotation supported this hypothesis. In addition, the spatio-temporal expression pattern of these seven ERFs, as detected using RT-qPCR and RNA-seq, suggested that they play a critical role in mediating the salt response and tolerance in a dynamic and tissue-specific manner. The results of this study provide a solid basis to explore the functions of the stress-related ERF TFs in Populus abiotic stress tolerance and development process.

Identification of the wheat C3H gene family and expression analysis of candidates associated with seed dormancy and germination

C3H zinc finger transcription factors play important roles in managing various biotic/abiotic stresses in Aarabidopsis, rice, and maize. The functions of these factors in wheat, however, remain largely unclear. We identified 88 TaC3H genes that were divided into four subfamilies in this analysis. Gene structure and conserved domain analyses indicate that most members of the same subfamily have similar structures. A total of 76 paralogous and 48 orthologous pairs were identified and Ka/Ks values were used to analyze replication relationships amongst wheat, rice, and Arabidopsis. Gene ontology (GO) annotation analysis showed that most TaC3H genes possessed molecular functions, while transcriptome results showed that the 88 TaC3H genes responded to water imbibition. Microarray data for 53 TaC3H genes were obtained and heat maps were generated; these results indicate that these genes are expressed in 13 wheat tissues. Subcellular localization prediction analysis indicates that most TaC3H genes are located in the nucleus. Promoter analysis indicates that most TaC3H genes contained cis-elements including ABRE, GARE-motif, and MBS, indicating that these can respond to various biotic/abiotic stresses. Transcriptome data and quantitative real-time PCR analysis of wheat cultivars with contrasting seed dormancy phenotypes show that five genes TaC3H4/-18/-37/-51/-72 were very likely involved in seed dormancy and germination. Exogenous ABA treatment further indicated that these five genes were responsive to ABA, suggesting that there may be a crosstalk between these genes and ABA signaling pathway in controlling seed dormancy and germination. These results provide a theoretical basis for subsequent studies on TaC3H gene function and also contribute to studies on the C3H gene in other species.

MYB Transcription Factor161 Mediates Feedback Regulation of Secondary wall-associated NAC-Domain1 Family Genes for Wood Formation

Wood formation is a complex process that involves cell differentiation, cell expansion, secondary wall deposition, and programmed cell death. We constructed a four-layer wood formation transcriptional regulatory network (TRN) in Populus trichocarpa (black cottonwood) that has four Secondary wall-associated NAC-Domain1 (PtrSND1) transcription factor (TF) family members as the top-layer regulators. We characterized the function of a MYB (PtrMYB161) TF in this PtrSND1-TRN, using transgenic P trichocarpa cells and whole plants. PtrMYB161 is a third-layer regulator that directly transactivates five wood formation genes. Overexpression of PtrMYB161 in P. trichocarpa (OE-PtrMYB161) led to reduced wood, altered cell type proportions, and inhibited growth. Integrative analysis of wood cell-based chromatin-binding assays with OE-PtrMYB161 transcriptomics revealed a feedback regulation system in the PtrSND1-TRN, where PtrMYB161 represses all four top-layer regulators and one second-layer regulator, PtrMYB021, possibly affecting many downstream TFs in, and likely beyond, the TRN, to generate the observed phenotypic changes. Our data also suggested that the PtrMYB161's repressor function operates through interaction of the base PtrMYB161 target-binding system with gene-silencing cofactors. PtrMYB161 protein does not contain any known negative regulatory domains. CRISPR-based mutants of PtrMYB161 in P. trichocarpa exhibited phenotypes similar to the wild type, suggesting that PtrMYB161's activator functions are redundant among many TFs. Our work demonstrated that PtrMYB161 binds to multiple sets of target genes, a feature that allows it to function as an activator as well as a repressor. The balance of the two functions may be important to the establishment of regulatory homeostasis for normal growth and development.

MYB134-RNAi poplar plants show reduced tannin synthesis in leaves but not roots, and increased susceptibility to oxidative stress

The importance of the poplar MYB134 gene in controlling condensed tannin (CT) biosynthesis was tested by suppressing its expression using RNA interference (RNAi). MYB134-RNAi plants grew normally but showed reduced accumulation of stress-induced CTs in leaves. RNA-seq analysis indicated that flavonoid- and CT-related genes, as well as additional CT regulators, were strongly and specifically down-regulated by MYB134 suppression. This confirmed that the primary MYB134 target is the leaf flavonoid and CT pathway. Root CT accumulation was not impacted by MYB suppression, suggesting that additional CT regulators are active in roots and emphasizing the complexity of the regulation of CTs in poplar. To test the effect of CT down-regulation on oxidative stress resistance, leaves of MYB134-RNAi and control plants were exposed to the reactive oxygen species generator methyl viologen. MYB134-RNAi leaves sustained significantly more photosystem II damage, as seen in reduced chlorophyll fluorescence, compared with wild-type leaves. MYB134-RNAi leaves also contained more hydrogen peroxide, a reactive oxygen species, compared with the wild type. Our data thus corroborate the hypothesis that CT can act as an antioxidant in vivo and protect against oxidative stress. Overall, MYB134 was shown to be a central player in the regulation of CT synthesis in leaves.

The R2R3-MYB gene family in banana (Musa acuminata): Genome-wide identification, classification and expression patterns

The R2R3-MYB genes comprise one of the largest transcription factor gene families in plants, playing regulatory roles in plant-specific developmental processes, defense responses and metabolite accumulation. To date MYB family genes have not yet been comprehensively identified in the major staple fruit crop banana. In this study, we present a comprehensive, genome-wide analysis of the MYB genes from Musa acuminata DH-Pahang (A genome). A total of 285 R2R3-MYB genes as well as genes encoding three other classes of MYB proteins containing multiple MYB repeats were identified and characterised with respect to structure and chromosomal organisation. Organ- and development-specific expression patterns were determined from RNA-Seq data. For 280 M. acuminata MYB genes for which expression was found in at least one of the analysed samples, a variety of expression patterns were detected. The M. acuminata R2R3-MYB genes were functionally categorised, leading to the identification of seven clades containing only M. acuminata R2R3-MYBs. The encoded proteins may have specialised functions that were acquired or expanded in Musa during genome evolution. This functional classification and expression analysis of the MYB gene family in banana establishes a solid foundation for future comprehensive functional analysis of MaMYBs and can be utilized in banana improvement programmes.

R2R3-MYB transcription factor PpMYB17 positively regulates flavonoid biosynthesis in pear fruit

PpMYB17 positively regulates flavonoid biosynthesis in pear fruit by activating PpCHS, PpCHI, PpF3H, and PpFLS in the flavonoid biosynthesis pathway independently of bHLH or WD40 cofactors in the MBW complex. Flavonoids are important secondary metabolites in plants. The flavonoid biosynthesis pathway is regulated by various transcription factors, with MYB transcription factors considered to be the key regulators. However, the regulation of flavonoid biosynthesis in the pear fruit has not been fully characterized. The R2R3-MYB transcription factor PpMYB17 was isolated from 'Red Zaosu' pear fruit and functionally characterized. An exposure to light upregulated PpMYB17 expression in the pear fruit. A phylogenetic analysis indicated PpMYB17 is related to the flavonol regulators. A subcellular localization assay suggested that PpMYB17 is a nuclear protein. Overexpression of PpMYB17 increased the flavonoid content of pear calli and Arabidopsis via the upregulated expression of structural genes in the flavonoid biosynthesis pathway, especially FLS. The LC-MS/MS analysis revealed most of the differentially accumulated flavonols, flavanones, flavones, isoflavones, and anthocyanins were significantly more abundant in PpMYB17-overexpressing calli than in wild-type calli. Moreover, PpMYB17 did not interact with PpbHLH3, PpbHLH33, or PpWD40 in a yeast system. Dual-luciferase assays demonstrated that PpMYB17 strongly activates the promoters of PpCHS, PpCHI, PpF3H, PpFLS, and PpUFGT which are key downstream genes in the flavonoid biosynthesis pathway, independently of the PpbHLH3 cofactor. These gene expression changes may enhance flavonoid biosynthesis in pear fruit. The data presented may be useful for further elucidating the flavonoid biosynthesis regulatory network, potentially leading to the development of new pear cultivars that produce fruits with increased flavonoid contents.

Genome-wide identification and expression analyses of R2R3-MYB transcription factor genes from two Orchid species

MYB transcription factors play important roles in different plant biological processes during plant growth, development and stress response. In this study, 101 (DoMYB1-101) and 99 (PaMYB1-99) R2R3-MYB genes were identified in the genomes of Dendrobium officinale and Phalaenopsis aphrodite, respectively. To classify the isolated candidate genes, the R2R3-MYB genes from A. thaliana were selected as references. As a result, all identified DoMYB and PaMYB genes were classified into 22 subfamilies. Phylogenetic analysis revealed that S21 had the largest number of members of all the subfamilies. The numbers of introns, exons and conserved sequences in all of the identified genes are different. In addition, 20 DoMYB genes from six subfamilies were selected for further analysis of tissue-specific expression and responses to various abiotic stresses treatments. The results showed that all of the DoMYB genes in S4 and S19 subfamilies exhibited the highest relative expression levels in flowers. And five DoMYB genes including DoMYB31, DoMYB40, DoMYB49, DoMYB52 and DoMYB54, responded to the stress response. These results may provide useful information for further studies of the R2R3-MYB gene family.

The Phoebe genome sheds light on the evolution of magnoliids

Lauraceae includes the genus Phoebe, and the family is linked to the evolution of magnoliids. We sequenced the genome of Phoebe bournei Nanmu. The assembled genome size was 989.19 Mb, with a contig N50 value of 2.05 Mb. A total of 28,198 protein-coding genes were annotated in P. bournei. Whole-genome duplication (WGD) analysis showed that Lauraceae has experienced two WGD events; the older WGD event occurred just before the divergence of Lauraceae and Magnoliales, and the more recent WGD was shared by all lineages of Lauraceae. The phylogenetic tree showed that magnoliids form a sister clade to monocots and eudicots. We also identified 63 MADS-box genes, including AGL12-like genes that may be related to the regulation of P. bournei roots and FIN219-like genes encoding GH3 proteins, which are involved in photomorphogenesis. SAUR50-like genes involved in light signal-mediated pedicel or stem development were also identified. Four ATMYB46- and three PtrEPSP-homologous genes related to lignin biosynthesis were identified. These genes may be associated with the formation of straight trunks in P. bournei. Overall, the P. bournei reference genome provides insight into the origin, evolution, and diversification of Phoebe and other magnoliids.

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

MYB proteins constitute one of the largest transcription factor families in plants, members of which are involved in various plant physiological and biochemical processes. Japanese plum (Prunus salicina) is one of the important stone fruit crops worldwide. To date, no comprehensive study of the MYB family in Japanese plum has been reported. In this study, we performed genome-wide analysis of MYB genes in Japanese plum including the phylogeny, gene structures, protein motifs, chromosomal locations, collinearity and expression patterns analysis. A total of 96 Japanese plum R2R3-MYB (PsMYB) genes were characterized and distributed on 8 chromosomes at various densities. Collinearity analysis indicated that the segmental duplication events played a crucial role in the expansion of PsMYB genes, and the interspecies synteny analysis revealed the orthologous gene pairs between Japanese plum and other four selected Rosaceae species. The 96 PsMYB genes could be classified into 27 subgroups based on phylogenetic topology, as supported by the conserved gene structures and motif compositions. Further comparative phylogenetic analysis revealed the functional divergence of MYB gene family during evolution, and three subgroups which included only Rasaceae MYB genes were identified. Expression analysis revealed the distinct expression profiles of the PsMYB genes, and further functional predictions found some of them might be associated with the plum fruit quality traits. Our researches provide a global insight into the organization, phylogeny, evolution and expression patterns of the PsMYB genes, and contribute to the greater understanding of their functional roles in Japanese plum.

Xylem systems genetics analysis reveals a key regulator of lignin biosynthesis in Populus deltoides

Despite the growing resources and tools for high-throughput characterization and analysis of genomic information, the discovery of the genetic elements that regulate complex traits remains a challenge. Systems genetics is an emerging field that aims to understand the flow of biological information that underlies complex traits from genotype to phenotype. In this study, we used a systems genetics approach to identify and evaluate regulators of the lignin biosynthesis pathway in Populus deltoides by combining genome, transcriptome, and phenotype data from a population of 268 unrelated individuals of P. deltoides The discovery of lignin regulators began with the quantitative genetic analysis of the xylem transcriptome and resulted in the detection of 6706 and 4628 significant local- and distant-eQTL associations, respectively. Among the locally regulated genes, we identified the R2R3-MYB transcription factor MYB125 (Potri.003G114100) as a putative trans-regulator of the majority of genes in the lignin biosynthesis pathway. The expression of MYB125 in a diverse population positively correlated with lignin content. Furthermore, overexpression of MYB125 in transgenic poplar resulted in increased lignin content, as well as altered expression of genes in the lignin biosynthesis pathway. Altogether, our findings indicate that MYB125 is involved in the control of a transcriptional coexpression network of lignin biosynthesis genes during secondary cell wall formation in P. deltoides.

Assessments of fine-scale spatial patterns of SNPs in an old-growth beech forest

The spatial patterns of non-neutral genetic variations at fine spatial scales and their possible associations with microenvironments have not been well-documented for tree populations. Based on 25-32 SNP markers, we examine whether non-neutral SNPs and their associations with microenvironments can be detected in FcMYB1603, a gene homologous to that encoding a protein induced by drought stress in Arabidopsis thaliana for the 166 adult trees in a 1-ha plot in a mature population of Fagus crenata. In the 83 individuals of a younger cohort of below canopy trees, the nonsynonymous SNP at locus FcMYB1603_684 exhibited a spatial signature representing a departure from the expected spatial patterns of neutral genetic variation. Evaluations of non-neutrality for this locus were robust against the potential risks of false positives due to the low number of SNP loci, a low criterion set for minor allele frequency, and any edge effect on the trees' spatial structure. An older cohort exhibited no signal of the existence of non-neutral genetic variation, suggesting that temporal fluctuation in the microenvironmental conditions on the forest floor may have exposed different cohorts to different magnitudes of selection pressure. Although genotypes of the locus showed a spatial association with a microenvironmental variable potentially related to soil moisture, the present study was subject to a limitation due to the generally low polymorphism of nonsynonymous loci within the single plot, which suggests that it will be important to replicate the study design in order to carry out research on fine-scale non-neutral genetic variations.

MYB Superfamily in Brassica napus: Evidence for Hormone-Mediated Expression Profiles, Large Expansion, and Functions in Root Hair Development

MYB proteins are involved in diverse important biological processes in plants. Herein, we obtained the MYB superfamily from the allotetraploid Brassica napus, which contains 227 MYB-related (BnMYBR/Bn1R-MYB), 429 R2R3-MYB (Bn2R-MYB), 22 R1R2R3-MYB (Bn3R-MYB), and two R1R2R2R1/2-MYB (Bn4R-MYB) genes. Phylogenetic analysis classified the Bn2R-MYBs into 43 subfamilies, and the BnMYBRs into five subfamilies. Sequence characteristics and exon/intron structures within each subfamily of the Bn2R-MYBs and BnMYBRs were highly conserved. The whole superfamily was unevenly distributed on 19 chromosomes and underwent unbalanced expansion in B. napus. Allopolyploidy between B. oleracea and B. rapa mainly contributed to the expansion in their descendent B. napus, in which B. rapa-derived genes were more retained. Comparative phylogenetic analysis of 2R-MYB proteins from nine Brassicaceae and seven non-Brassicaceae species identified five Brassicaceae-specific subfamilies and five subfamilies that are lacking from the examined Brassicaceae species, which provided an example for the adaptive evolution of the 2R-MYB gene family alongside angiosperm diversification. Ectopic expression of four Bn2R-MYBs under the control of the viral CaMV35S and/or native promoters could rescue the lesser root hair phenotype of the Arabidopsis thaliana wer mutant plants, proving the conserved negative roles of the 2R-MYBs of the S15 subfamily in root hair development. RNA-sequencing data revealed that the Bn2R-MYBs and BnMYBRs had diverse transcript profiles in roots in response to the treatments with various hormones. Our findings provide valuable information for further functional characterizations of B. napusMYB genes.

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

As one of the largest families of transcription factors (TFs) in plants, R2R3-MYB proteins play crucial roles in regulating a series of plant-specific biological processes. Although the diversity of plant R2R3-MYB TFs has been studied previously, the processes and mechanisms underlying the expansion of these proteins remain unclear. Here, we performed evolutionary analyses of plant R2R3-MYB TFs with dense coverage of streptophyte algae and embryophytes. Our analyses revealed that ancestral land plants exhibited 10 subfamilies of R2R3-MYB proteins, among which orthologs of seven subfamilies were present in chlorophytes and charophycean algae. We found that asymmetric gene duplication events in different subfamilies account for the expansion of R2R3-MYB proteins in embryophytes. We further discovered that the largest subfamily of R2R3-MYBs in land plants, subfamily VIII, emerged in the common ancestor of Zygnematophyceae and embryophytes. During plant terrestrialization, six duplication events gave rise to seven clades of subfamily VIII. Subsequently, this TF subfamily showed a tendency for expansion in bryophytes, lycophytes, and ferns and extensively diversified in ancestral gymnosperms and angiosperms in clades VIII-A-1, VIII-D, and VIII-E. In contrast to subfamily VIII, other subfamilies of R2R3-MYB TFs have remained less expanded across embryophytes. The findings regarding phylogenetic analyses, auxiliary motifs, and DNA-binding specificities provide insight into the evolutionary history of plant R2R3-MYB TFs and shed light on the mechanisms underlying the extensive expansion and subsequent sub- and neofunctionalization of these proteins.

The Role of Stress-Responsive Transcription Factors in Modulating Abiotic Stress Tolerance in Plants

Abiotic stresses, such as drought, high temperature, and salinity, affect plant growth and productivity. Furthermore, global climate change may increase the frequency and severity of abiotic stresses, suggesting that development of varieties with improved stress tolerance is critical for future sustainable crop production. Improving stress tolerance requires a detailed understanding of the hormone signaling and transcriptional pathways involved in stress responses. Abscisic acid (ABA) and jasmonic acid (JA) are key stress-response hormones in plants, and some stress-responsive transcription factors such as ABFs and MYCs function as direct components of ABA and JA signaling, playing a pivotal role in plant tolerance to abiotic stress. In addition, extensive studies have identified other stress-responsive transcription factors belonging to the NAC, AP2/ERF, MYB, and WRKY families that mediate plant response and tolerance to abiotic stress. These suggest that transcriptional regulation of stress-responsive genes is an essential step to determine the mechanisms underlying plant stress responses and tolerance to abiotic stress, and that these transcription factors may be important targets for development of crops with enhanced abiotic stress tolerance. In this review, we briefly describe the mechanisms underlying plant abiotic stress responses, focusing on ABA and JA metabolism and signaling pathways. We then summarize the diverse array of transcription factors involved in plant responses to abiotic stress, while noting their potential applications for improvement of stress tolerance.

MYB transcription factor PdMYB118 directly interacts with bHLH transcription factor PdTT8 to regulate wound-induced anthocyanin biosynthesis in poplar

BACKGROUND

R2R3-MYB transcription factors (TFs) play important roles in plant growth and development, and response to biotic and abiotic stresses. However, their regulatory mechanisms in wound-induced anthocyanin biosynthesis in woody plants are largely unknown.

RESULTS

In this work, we report that expression of anthocyanin biosynthesis genes (ABGs) were activated by PdMYB118, a MYB TF encoding gene from Populus deltoids, and the activation of PdMYB118 was significantly enhanced by PdTT8, a bHLH protein, through its direct interaction with PdMYB118. PdMYB118 and some ABGs were evidently induced by wound induction and methyl jasmonate (MeJA) treatment. Overexpression of PdMYB118 promoted anthocyanin accumulation in transgenic poplar upon wound induction. Furthermore, a poplar JASMONATE ZIM-domain (JAZ) protein, PtrJAZ1, repressed the transcriptional function of PdMYB118/PdTT8 complex by binding to PdTT8, and wound stimulated the biosynthesis of jasmonic acid (JA) and the degradation of PtrJAZ1.

CONCLUSIONS

Based on these observations, we proposed that PtrJAZ1 degradation triggered the expression of ABGs, leading to increased biosynthesis of anthocyanins in the wounded leaves of transgenic poplar. Therefore, our findings not only illustrate the crucial role of PdMYB118 in wound-induced anthocyanin biosynthesis in poplar, but also provide a molecular basis for the genetic engineering of colorful tree species.

CmMYB8 encodes an R2R3 MYB transcription factor which represses lignin and flavonoid synthesis in chrysanthemum

R2R3-MYB transcription factors are important regulators of the growth and development of plants. Here, CmMYB8 a chrysanthemum gene encoding an R2R3-MYB transcription factor, was isolated and functionally characterized. The gene was transcribed throughout the plant, but most strongly in the stem. When CmMYB8 was over-expressed, a number of genes encoding components of lignin synthesis were down-regulated, and the plants' lignin content was reduced. The composition of the lignin in the transgenic plants was also altered, and its S/G ratio was reduced. A further consequence of the over-expression of CmMYB8 was to lessen the transcript abundance of key genes involved in flavonoid synthesis, resulting in a reduced accumulation of flavonoids. The indication is that the CmMYB8 protein participates in the negative regulation of both lignin and flavonoid synthesis.

Genome-Wide Identification, Classification, and Expression Profiling Reveals R2R3-MYB Transcription Factors Related to Monoterpenoid Biosynthesis in Osmanthus fragrans

Osmanthus fragrans is widely grown for the purpose of urban greening and the pleasant aroma emitted from its flowers. The floral scent is determined by several monoterpenoid volatiles, such as linalool and its oxides, which are a few of the most common volatiles and the main components of the essential oils in most sweet osmanthus cultivars. In addition, the relative contents of cis- and trans-linalool oxide (furan) may affect the aromas and quality of the essential oils. MYB proteins represent the largest family of transcription factors in plants and participate in regulating secondary metabolites. Several cis-elements, especially AC-rich regions, are known to be bound by 2R-MYBs and could be found in the promoter of the enzyme genes in the terpenoid metabolic pathway. However, there has to date been no investigation into the 2R-MYB family genes involved in regulating terpenoid biosynthesis in O. fragrans. Here, 243 non-redundant 2R-MYB proteins were grouped into 33 clusters based on the phylogeny and exon-intron distribution. These genes were unevenly distributed on 23 chromosomes. Ka/Ks analysis showed that the major mode of 2R-MYB gene evolution was purifying selection. Expression analysis indicated that 2R-MYB genes in O. fragrans exhibited varied expression patterns. A total of 35 OfMYBs representing the highest per kilobase per million mapped reads in the flower were selected for quantitative real-time PCR analysis. The correlation analysis between the expression level and the contents of fragrant compounds at different flowering stages suggested that OfMYB19/20 exhibited remarkably positive correlation with the accumulation of cis-linalool oxides. OfMYB51/65/88/121/137/144 showed significantly negative correlations with one or more linalool oxides. Characterization of these proteins revealed that OfMYB19 and OfMYB137 were localized in the nuclei, but did not show transcriptional activation in the yeast system, which suggested that they may be bound to other transcription factors to exert regulatory functions. These findings provide useful information for further functional investigation of the 2R-MYBs and offer a foundation for clarifying the 2R-MYB transcription factors involved in the molecular mechanism of the regulation of monoterpenoid biosynthesis in Osmanthus fragrans.

Transcriptional Regulation of Sorghum Stem Composition: Key Players Identified Through Co-expression Gene Network and Comparative Genomics Analyses

Most sorghum biomass accumulates in stem secondary cell walls (SCW). As sorghum stems are used as raw materials for various purposes such as feed, energy and fiber reinforced polymers, identifying the genes responsible for SCW establishment is highly important. Taking advantage of studies performed in model species, most of the structural genes contributing at the molecular level to the SCW biosynthesis in sorghum have been proposed while their regulatory factors have mostly not been determined. Validation of the role of several MYB and NAC transcription factors in SCW regulation in Arabidopsis and a few other species has been provided. In this study, we contributed to the recent efforts made in grasses to uncover the mechanisms underlying SCW establishment. We reported updated phylogenies of NAC and MYB in 9 different species and exploited findings from other species to highlight candidate regulators of SCW in sorghum. We acquired expression data during sorghum internode development and used co-expression analyses to determine groups of co-expressed genes that are likely to be involved in SCW establishment. We were able to identify two groups of co-expressed genes presenting multiple evidences of involvement in SCW building. Gene enrichment analysis of MYB and NAC genes provided evidence that while NAC SECONDARY WALL THICKENING PROMOTING FACTOR NST genes and SECONDARY WALL-ASSOCIATED NAC DOMAIN PROTEIN gene functions appear to be conserved in sorghum, NAC master regulators of SCW in sorghum may not be as tissue compartmentalized as in Arabidopsis. We showed that for every homolog of the key SCW MYB in Arabidopsis, a similar role is expected for sorghum. In addition, we unveiled sorghum MYB and NAC that have not been identified to date as being involved in cell wall regulation. Although specific validation of the MYB and NAC genes uncovered in this study is needed, we provide a network of sorghum genes involved in SCW both at the structural and regulatory levels.

The R2R3-MYB transcription factor family in Taxus chinensis: identification, characterization, expression profiling and posttranscriptional regulation analysis

The MYB transcription factor family is one of the largest gene families playing regulatory roles in plant growth and development. The MYB family has been studied in a variety of plant species but has not been reported in Taxus chinensis. Here we identified 72 putative R2R3-MYB genes in T. chinensis using a comprehensive analysis. Sequence features, conversed domains and motifs were characterized. The phylogenetic analysis showed TcMYBs and AtMYBs were clustered into 36 subgroups, of which 24 subgroups included members from T. chinensis and Arabidopsis thaliana, while 12 subgroups were specific to one species. This suggests the conservation and specificity in structure and function of plant R2R3-MYBs. The expression of TcMYBs in various tissues and different ages of xylem were investigated. Additionally, miRNA-mediated posttranscriptional regulation analysis revealed that TcMYBs were the targets of miR858, miR159 and miR828, suggesting the posttranscriptional regulation of MYBs is highly conserved in plants. The results provide a basis for further study the role of TcMYBs in the regulation of secondary metabolites of T. chinensis.

Genome-Wide Analysis of the MYB Transcription Factor Superfamily in Physcomitrella patens

MYB transcription factors (TFs) are one of the largest TF families in plants to regulate numerous biological processes. However, our knowledge of the MYB family in Physcomitrella patens is limited. We identified 116 MYB genes in the P. patens genome, which were classified into the R2R3-MYB, R1R2R3-MYB, 4R-MYB, and MYB-related subfamilies. Most R2R3 genes contain 3 exons and 2 introns, whereas R1R2R3 MYB genes contain 10 exons and 9 introns. N3R-MYB (novel 3RMYB) and NR-MYBs (novel RMYBs) with complicated gene structures appear to be novel MYB proteins. In addition, we found that the diversity of the MYB domain was mainly contributed by domain shuffling and gene duplication. RNA-seq analysis suggested that MYBs exhibited differential expression to heat and might play important roles in heat stress responses, whereas CCA1-like MYB genes might confer greater flexibility to the circadian clock. Some R2R3-MYB and CCA1-like MYB genes are preferentially expressed in the archegonium and during the transition from the chloronema to caulonema stage, suggesting their roles in development. Compared with that of algae, the numbers of MYBs have significantly increased, thus our study lays the foundation for further exploring the potential roles of MYBs in the transition from aquatic to terrestrial environments.

Genome-wide analysis and expression profiles of the StR2R3-MYB transcription factor superfamily in potato (Solanum tuberosum L.)

MYB transcription factors comprise one of the largest families in plant kingdom, which play a variety of functions in plant developmental processes and defence responses, the R2R3-MYB members are the predominant form found in higher plants. In the present study, a total of 111 StR2R3-MYB transcription factors were identified and further phylogenetically classified into 31 subfamilies, as supported by highly conserved gene structures and motifs. Collinearity analysis showed that the segmental duplication events played a crucial role in the expansion of StR2R3-MYB gene family. Synteny analysis indicated that 37 and 13 StR2R3-MYB genes were orthologous to Arabidopsis and wheat (Triticum aestivum), respectively, and these gene pairs have evolved under strong purifying selection. RNA-seq data from different tissues and abiotic stresses revealed tissue-preferential and abiotic stress-responsive StR2R3-MYB genes. We further analyzed StR2R3-MYB genes might be involved in anthocyanin biosynthesis and drought stress by using RNA-seq data of pigmented tetraploid potato cultivars and drought-sensitive and -tolerant tetraploid potato cultivars under drought stress, respectively. Moreover, EAR motifs were found in 21 StR2R3-MYB proteins and 446 pairs of proteins were predicted to interact with 21 EAR motif-containing StR2R3-MYB proteins by constructing the interaction network with medium confidence (0.4). Additionally, Gene Ontology (GO) analysis of the 21 EAR motif-containing StR2R3-MYB proteins was performed to further investigate their functions. This work will facilitate future biologically functional studies of potato StR2R3-MYB transcription factors and enrich the knowledge of MYB superfamily genes in plant species.

Identification and characterization of circular RNAs during wood formation of poplars in acclimation to low nitrogen availability

Circular RNA (circRNA) identification and expression profiles, and construction of circRNAs-miRNAs-mRNAs networks indicates that circRNAs are involved in wood formation of poplars in acclimation to low nitrogen availability. Circular RNAs (circRNAs) are covalently closed non-coding RNAs that play pivotal roles in various biological processes. However, circRNAs' roles in wood formation of poplars in acclimation to low nitrogen (N) availability are currently unknown. Here, we undertook a systematic identification and characterization of circRNAs in the wood of Populus × canescens exposed to either 50 (low N) or 500 (normal N) µM NH₄NO₃ using rRNA-depleted RNA-sequencing. A total of 2,509 unique circRNAs were identified, and 163 (ca. 6.5%) circRNAs were significantly differentially expressed (DE) under low N condition. We observed a positive correlation between the expression patterns of DE circRNAs and their hosting protein-coding genes. Moreover, circRNAs-miRNAs-mRNAs' networks were identified in the wood of poplars under low N availability. For instance, upregulated several circRNAs, such as circRNA1226, circRNA 1732, and circRNA392 induced increases in nuclear factor Y, subunit A1-A (NFYA1-A), NFYA1-B, and NFYA10 transcript levels via the mediation of miR169b members, which is in line with reduced xylem width and cell layers of the xylem in the wood of low N-supplied poplars. Upregulation of circRNA1006, circRNA1344, circRNA1941, circRNA901, and circRNA146 caused increased transcript level of MYB61 via the mediation of a miR5021 member, corresponding well to the higher lignin concentration in the wood of low N-treated poplars. Overall, these results indicated that DE circRNAs play an essential role in regulating gene expression via circRNAs-miRNAs-mRNAs' networks to modulate wood anatomical and chemical properties of poplars in acclimation to low N availability.

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.

Identification of functional single nucleotide polymorphism of Populus trichocarpa PtrEPSP-TF and determination of its transcriptional effect

In plants, the phenylpropanoid pathway is responsible for the synthesis of a diverse array of secondary metabolites that include lignin monomers, flavonoids, and coumarins, many of which are essential for plant structure, biomass recalcitrance, stress defense, and nutritional quality. Our previous studies have demonstrated that Populus trichocarpa PtrEPSP-TF, an isoform of 5-enolpyruvylshikimate 3-phosphate (EPSP) synthase, has transcriptional activity and regulates phenylpropanoid biosynthesis in Populus. In this study, we report the identification of single nucleotide polymorphism (SNP) of PtrEPSP-TF that defines its functionality. Populus natural variants carrying this SNP were shown to have reduced lignin content. Here, we demonstrated that the SNP-induced substitution of 142nd amino acid (PtrEPSP-TFD142E) dramatically impairs the DNA-binding and transcriptional activity of PtrEPSP-TF. When introduced to a monocot species rice (Oryza sativa) in which an EPSP synthase isoform with the DNA-binding helix-turn-helix (HTH) motif is absent, the PtrEPSP-TF, but not PtrEPSP-TFD142E, activated genes in the phenylpropanoid pathway. More importantly, heterologous expression of PtrEPSP-TF uncovered five new transcriptional regulators of phenylpropanoid biosynthesis in rice. Collectively, this study identifies the key amino acid required for PtrEPSP-TF functionality and provides a strategy to uncover new transcriptional regulators in phenylpropanoid biosynthesis.

The poplar R2R3 MYB transcription factor PtrMYB94 coordinates with abscisic acid signaling to improve drought tolerance in plants

In plants, R2R3 MYB transcription factors (TFs) consist of one large gene family and are involved in the regulation of many developmental processes and various stresses. However, the functions of most of MYB TFs in woody plants remain unknown. Here, PtrMYB94, an R2R3 MYB TF from Populus trichocarpa, is characterized to be involved in the regulation of drought responses and abscisic acid (ABA) signaling. PtrMYB94 encodes a nuclear-localized R2R3 MYB TF. RT-PCR results showed that the PtrMYB94 transcripts were relatively abundant in leaves and stems, and were induced rapidly in response to dehydration stress. Overexpression of PtrMYB94 improved plant drought responses, suggesting that this MYB TF may functionally regulate poplar adaptability to drought stress. Furthermore, the analysis of transcriptional expression and PtrMYB94 promoter: GUS activity showed that PtrMYB94 responded to ABA induction. PtrMYB94-overexpressing plants exhibited the inhibition of seed germination compared with the wild-type (WT) control under ABA exposure condition. The ABA content was evidently increased in the PtrMYB94-overexpressing plants relative to the WT plants. In addition, transcript levels of several ABA- and drought-responsive genes, such as ABA1 and DREB2B, were up-regulated. Taken together, our results suggest that PtrMYB94 is involved in an ABA-dependent drought stress regulation in Populus.

CRISPR/Cas9 Genome Editing Technology: A Valuable Tool for Understanding Plant Cell Wall Biosynthesis and Function

For the past 5 years, clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) technology has appeared in the molecular biology research spotlight. As a game-changing player in genome editing, CRISPR/Cas9 technology has revolutionized animal research, including medical research and human gene therapy as well as plant science research, particularly for crop improvement. One of the most common applications of CRISPR/Cas9 is to generate genetic knock-out mutants. Recently, several multiplex genome editing approaches utilizing CRISPR/Cas9 were developed and applied in various aspects of plant research. Here we summarize these approaches as they relate to plants, particularly with respect to understanding the biosynthesis and function of the plant cell wall. The plant cell wall is a polysaccharide-rich cell structure that is vital to plant cell formation, growth, and development. Humans are heavily dependent on the byproducts of the plant cell wall such as shelter, food, clothes, and fuel. Genes involved in the assembly of the plant cell wall are often highly redundant. To identify these redundant genes, higher-order knock-out mutants need to be generated, which is conventionally done by genetic crossing. Compared with genetic crossing, CRISPR/Cas9 multi-gene targeting can greatly shorten the process of higher-order mutant generation and screening, which is especially useful to characterize cell wall related genes in plant species that require longer growth time. Moreover, CRISPR/Cas9 makes it possible to knock out genes when null T-DNA mutants are not available or are genetically linked. Because of these advantages, CRISPR/Cas9 is becoming an ideal and indispensable tool to perform functional studies in plant cell wall research. In this review, we provide perspectives on how to design CRISPR/Cas9 to achieve efficient gene editing and multi-gene targeting in plants. We also discuss the recent development of the virus-based CRISPR/Cas9 system and the application of CRISPR/Cas9 to knock in genes. Lastly, we summarized current progress on using CRISPR/Cas9 for the characterization of plant cell wall-related genes.

The BpMYB4 Transcription Factor From Betula platyphylla Contributes Toward Abiotic Stress Resistance and Secondary Cell Wall Biosynthesis

The MYB (v-myb avian myeloblastosis viral oncogene homolog) family is one of the largest transcription factor families in plants, and is widely involved in the regulation of plant metabolism. In this study, we show that a MYB4 transcription factor, BpMYB4, identified from birch (Betula platyphylla Suk.) and homologous to EgMYB1 from Eucalyptus robusta Smith and ZmMYB31 from Zea mays L. is involved in secondary cell wall synthesis. The expression level of BpMYB4 was higher in flowers relative to other tissues, and was induced by artificial bending and gravitational stimuli in developing xylem tissues. The expression of this gene was not enriched in the developing xylem during the active season, and showed higher transcript levels in xylem tissues around sprouting and near the dormant period. BpMYB4 also was induced express by abiotic stress. Functional analysis indicated that expression of BpMYB4 in transgenic Arabidopsis (Arabidopsis thaliana) plants could promote the growth of stems, and result in increased number of inflorescence stems and shoots. Anatomical observation of stem sections showed lower lignin deposition, and a chemical contents test also demonstrated increased cellulose and decreased lignin content in the transgenic plants. In addition, treatment with 100 mM NaCl and 200 mM mannitol resulted in the germination rate of the over-expressed lines being higher than that of the wild-type seeds. The proline content in transgenic plants was higher than that in WT, but MDA content was lower than that in WT. Further investigation in birch using transient transformation techniques indicated that overexpression of BpMYB4 could scavenge hydrogen peroxide and O₂ ^.- and reduce cell damage, compared with the wild-type plants. Therefore, we believe that BpMYB4 promotes stem development and cellulose biosynthesis as an inhibitor of lignin biosynthesis, and has a function in abiotic stress resistance.

Genome-wide identification and characterization of R2R3-MYB genes in Medicago truncatula

MYB is a large family of plant transcription factors. Its function has been identified in several plants, while there are few reports in Medicago truncatula. In this study, we used RNA-seq data to analyze and identify R2R3-MYB genes in the genome of Medicago truncatula. Phylogenetic analysis classified 150 MtMYB genes into 21 subfamilies with homologs. Out of the 150 MtMYB genes, 139 were distributed among 8 chromosomes, with tandem duplications (TD) and segment duplications (SD). Microarray data were used for functional analysis of the MtMYB genes during growth and developmental processes providing evidence for a role in tissues differentiation, seed development processes, and especially the nodulation process. Furthermore, we investigated the expression of MtMYB genes in response to abiotic stresses using RNA-seq data, which confirmed the critical roles in signal transduction and regulation processes under abiotic stress. We used quantitative real-time PCR (qRT-PCR) to validate expression profiles. The expression pattern of M. truncatula MYB genes under different abiotic stress conditions suggest that some may play a major role in cross-talk among different signal transduction pathways in response to abiotic stresses. Our study will serve as a foundation for future research into the molecular function of M. truncatula R2R3-MYB genes.

Dissecting the Genome-Wide Evolution and Function of R2R3-MYB Transcription Factor Family in Rosa chinensis

Rosa chinensis, an important ancestor species of Rosa hybrida, the most popular ornamental plant species worldwide, produces flowers with diverse colors and fragrances. The R2R3-MYB transcription factor family controls a wide variety of plant-specific metabolic processes, especially phenylpropanoid metabolism. Despite their importance for the ornamental value of flowers, the evolution of R2R3-MYB genes in plants has not been comprehensively characterized. In this study, 121 predicted R2R3-MYB gene sequences were identified in the rose genome. Additionally, a phylogenomic synteny network (synnet) was applied for the R2R3-MYB gene families in 35 complete plant genomes. We also analyzed the R2R3-MYB genes regarding their genomic locations, Ka/Ks ratio, encoded conserved motifs, and spatiotemporal expression. Our results indicated that R2R3-MYBs have multiple synteny clusters. The RcMYB114a gene was included in the Rosaceae-specific Cluster 54, with independent evolutionary patterns. On the basis of these results and an analysis of RcMYB114a-overexpressing tobacco leaf samples, we predicted that RcMYB114a functions in the phenylpropanoid pathway. We clarified the relationship between R2R3-MYB gene evolution and function from a new perspective. Our study data may be relevant for elucidating the regulation of floral metabolism in roses at the transcript level.

Identification and Analysis of the GASR Gene Family in Common Wheat (Triticum aestivum L.) and Characterization of TaGASR34, a Gene Associated With Seed Dormancy and Germination

Seed dormancy and germination are important agronomic traits in wheat (Triticum aestivum L.) because they determine pre-harvest sprouting (PHS) resistance and thus affect grain production. These processes are regulated by Gibberellic Acid-Stimulated Regulator (GASR) genes. In this study, we identified 37 GASR genes in common wheat, which were designated TaGASR1-37. Moreover, we identified 40 pairs of paralogous genes, of which only one had a Ka/Ks value greater than 1, indicating that most TaGASR genes have undergone negative selection. Chromosomal location and duplication analysis revealed 25 pairs of segmentally duplicated genes and seven pairs of tandemly duplicated genes, suggesting that large-scale duplication events may have contributed to the expansion of TaGASR gene family. Microarray analysis of the expression of 18 TaGASR genes indicated that these genes play diverse roles in different biological processes. Using wheat varieties with contrasting seed dormancy phenotypes, we investigated the expression patterns of TaGASR genes and the corresponding seed germination index phenotypes in response to water imbibition, exogenous ABA and GA treatment, and low- and high-temperature treatment. Based on these data, we identified the TaGASR34 gene as potentially associated with seed dormancy and germination. Further, we used a SNP mutation of the TaGASR34 promoter (-16) to develop the CAPS marker GS34-7B, which was then used to validate the association of TaGASR34 with seed dormancy and germination by evaluating two natural populations across environments. Notably, the frequency of the high-dormancy GS34-7Bb allele was significantly lower than that of the low-dormancy GS34-7Ba allele, implying that the favorable GS34-7Bb allele has not previously been used in wheat breeding. These results provide valuable information for further functional analysis of TaGASR genes and present a useful gene and marker combination for future improvement of PHS resistance in wheat.

Phosphorylation of LTF1, an MYB Transcription Factor in Populus, Acts as a Sensory Switch Regulating Lignin Biosynthesis in Wood Cells

Lignin is specifically deposited in plant secondary cell walls, and initiation of lignin biosynthesis is regulated by a variety of developmental and environmental signals. However, the mechanisms governing the regulation of lignin biosynthesis remain to be elucidated. In this study, we identified a lignin biosynthesis-associated transcription factor (LTF) from Populus, LTF1, which binds the promoter of a key lignin biosynthetic gene encoding 4-coumarate-CoA ligase (4CL). We showed that LTF1 in its unphosphorylated state functions as a regulator restraining lignin biosynthesis. When LTF1 becomes phosphorylated by PdMPK6 in response to external stimuli such as wounding, it undergoes degradation through a proteasome pathway, resulting in activation of lignification. Expression of a phosphorylation-null mutant version of LTF1 led to stable protein accumulation and persistent attenuation of lignification in wood cells. Taken together, our study reveals a mechanism whereby LTF1 phosphorylation acts as a sensory switch to regulate lignin biosynthesis in response to environmental stimuli. The discovery of novel modulators and mechanisms modifying lignin biosynthesis has important implications for improving the utilization of cell-wall biomass.

Genome-wide identification, expression and functional analysis of Populus xylogen-like genes

As an extracellular arabinogalactan protein (AGP) containing a non-specific lipid transfer protein (nsLTP) domain, xylogen mediates the local intercellular communication required for tracheary element (TE) differentiation in Zinnia cell culture. Although XYLP (xylogen-like protein) gene families have been reported in Arabidopsis and rice, no comprehensive analysis has been performed in woody plants. In this work, 31 XYLP genes in five phylogenetic groups were identified from Populus trichocarpa genome and a comprehensive bioinformatic analysis including gene and protein structures, chromosomal locations and duplication events were conducted. In-silico data and qRT-PCR results indicated that PtXYLP1 is predominantly expressed in poplar apex, young leaves and roots, while PtXYLP2 is uniformly expressed across a variety of tissues with a low abundance. Analysis on PtXYLP1pro:GUS and PtXYLP2pro:GUS in Arabidopsis revealed their differential expression patterns during seed germination and specific inductions by exogenously applied phytohormones including auxin, cytokinin and GA. When overexpressed in Arabidopsis, PtXYLP1 but not PtXYLP2 resulted in cotyledons with defective venation patterns and interrupted secondary (2°) vein loops, which phenotype was underpinned by the down-regulation of genes indispensably required by embryonic venation development at procambium and/or vessel level.

Ectopic Expression of a R2R3-MYB Transcription Factor Gene LjaMYB12 from Lonicera japonica Increases Flavonoid Accumulation in Arabidopsis thaliana

Lonicera japonica Thunb. is a widely used medicinal plant and is rich in a variety of active ingredients. Flavonoids are one of the important components in L. japonica and their content is an important indicator for evaluating the quality of this herb. To study the regulation of flavonoid biosynthesis in L. japonica, an R2R3-MYB transcription factor gene LjaMYB12 was isolated and characterized. Bioinformatics analysis indicated that LjaMYB12 belonged to the subgroup 7, with a typical R2R3 DNA-binding domain and conserved subgroup 7 motifs. The transcriptional level of LjaMYB12 was proportional to the total flavonoid content during the development of L. japonica flowers. Subcellular localization analysis revealed that LjaMYB12 localized to the nucleus. Transactivation activity assay indicated that LjaMYB12 was a transcriptional activator. Then, ectopic expression of LjaMYB12 in Arabidopsis could increase PAL activity and flavonoid content and promote transcription of a range of flavonoid biosynthetic genes. Interestingly, the fold changes of downstream genes in the flavonoid biosynthetic pathway were significantly higher than that of the upstream genes, which suggested that LjaMYB12 may have different regulatory patterns for the upstream and downstream pathways of flavonoid biosynthesis. The results provided here will effectively facilitate the study of subgroup 7 MYBs and transcriptional regulation of flavonoid biosynthesis in L. japonica.

Expression Analysis of the NAC Transcription Factor Family of Populus in Response to Salt Stress

Research Highlights: Sequence phylogeny, genome organisation, gene structure, conserved motifs, promoter cis-element and expression profiling of poplar NACs related to salt stress were detected. In addition, expression of two salt-induced NACs was analysed. Background and Objectives: NAC transcription factor (TF) proteins are involved in a wide range of functions during plant development and stress-related endurance processes. To understand the function of Populus NAC TFs in salt stress tolerance, we characterised the structure and expression profile of a total of 289 NAC members. Materials and Methods: Sequence phylogeny, genome organisation, gene structure, motif composition and promoter cis-element were detected using bioinformatics. The expression pattern of Populus NAC TFs under salt stress was also detected using RNA-Seq and RT-qPCR. Results: Synteny analysis showed that 46 and 37 Populus NAC genes were involved in whole-genome duplication and tandem duplication events, respectively. The expression pattern of Populus NAC TFs under salt stress showed the expression of the 289 PtNACs of 84K poplar was induced. Similar expression trends of NACs were found in Populus simonii × P. nigra T. S. Hwang et Liang and Arabidopsis thaliana (L.) Heynh. Conclusions: The correlation analysis showed that the expression of two differentially expressed NAC genes PtNAC024 and PtNAC182 was significantly associated with most of the 63 differentially expressed genes tested. The expression of PtNAC024 and PtNAC182 in different tissues was also analysed in silico and different expression patterns were found. Together, this study provides a solid basis to explore stress-related NAC TF functions in Populus salt tolerance and development.

Investigation of the effect of UV-B light on Arabidopsis MYB4 (AtMYB4) transcription factor stability and detection of a putative MYB4-binding motif in the promoter proximal region of AtMYB4

Here, we have investigated the possible effect of UV-B light on the folding/unfolding properties and stability of Arabidopsis thaliana MYB4 (AtMYB4) transcription factor in vitro by using biophysical approaches. Urea-induced equilibrium unfolding analyses have shown relatively higher stability of the wild-type recombinant AtMYB4 protein than the N-terminal deletion forms after UV-B exposure. However, as compared to wild-type form, AtMYB4Δ2 protein, lacking both the two N-terminal MYB domains, showed appreciable alteration in the secondary structure following UV-B exposure. UV-B irradiated AtMYB4Δ2 also displayed higher propensity of aggregation in light scattering experiments, indicating importance of the N-terminal modules in regulating the stability of AtMYB4 under UV-B stress. DNA binding assays have indicated specific binding activity of AtMYB4 to a putative MYB4 binding motif located about 212 bp upstream relative to transcription start site of AtMYB4 gene promoter, while relatively weak DNA binding activity was detected for another putative MYB4 motif located at -908 bp in AtMYB4 promoter. Gel shift and fluorescence anisotropy studies have shown increased binding affinity of UV-B exposed AtMYB4 to the promoter proximal MYB4 motif. ChIP assay has revealed binding of AtMYB4 to the promoter proximal (-212 position) MYB4 motif (ACCAAAC) in vivo. Docking experiments further revealed mechanistic detail of AtMYB4 interaction with the putative binding motifs. Overall, our results have indicated that the N-terminal 62-116 amino acid residues constituting the second MYB domain plays an important role in maintaining the stability of the C-terminal region and the overall stability of the protein, while a promoter proximal MYB-motif in AtMYB4 promoter may involve in the regulation of its own expression under UV-B light.

R2R3-MYB transcription factor MYB6 promotes anthocyanin and proanthocyanidin biosynthesis but inhibits secondary cell wall formation in Populus tomentosa

The secondary cell wall is an important carbon sink in higher plants and its biosynthesis requires coordination of metabolic fluxes in the phenylpropanoid pathway. In Arabidopsis (Arabidopsis thaliana), MYB75 and the KNOX transcription factor KNAT7 form functional complexes to regulate secondary cell wall formation in the inflorescence stem. However, the molecular mechanism by which these transcription factors control different branches of the phenylpropanoid pathway remains poorly understood in woody species. We isolated an R2R3-MYB transcription factor MYB6 from Populus tomentosa and determined that it was expressed predominately in young leaves. Overexpression of MYB6 in transgenic poplar upregulated flavonoid biosynthetic gene expression, resulting in significantly increased accumulation of anthocyanin and proanthocyanidins. MYB6-overexpression plants showed reduced secondary cell wall deposition, accompanied by repressed expression of secondary cell wall biosynthetic genes. We further showed that MYB6 interacted physically with KNAT7 and formed functional complexes that acted to repress secondary cell wall development in poplar and Arabidopsis. The results provide an insight into the transcriptional mechanisms involved in the regulation of the metabolic fluxes between the flavonoid and lignin biosynthetic pathways in poplar.

PdMYB118, isolated from a red leaf mutant of Populus deltoids, is a new transcription factor regulating anthocyanin biosynthesis in poplar

A new anthocyanin biosynthesis transcription factor PdMYB118, which could be used for the genetic engineering of colorful tree species, was indentified from a red leaf mutant of Populus deltoids. In higher plants, the biosynthesis of anthocyanins is regulated by several classes of transcription factors (TFs), including R2R3-MYB, bHLH and WD-repeat proteins. In this work, we isolated an MYB gene regulating anthocyanin biosynthesis from a red leaf mutant of Populus deltoids, which accumulated more anthocyanins in the leaves and showed higher expression levels of anthocyanin biosynthesis genes than did the wild type. Gene expression analyses of all TFs regulating anthocyanin biosynthesis demonstrated that only a MYB118 homologous gene, PdMYB118, was up-regulated in the mutant compared with the wide type. Subcellular localization analyses in poplar leaf mesophyll protoplasts showed that PdMYB118-YFP fusion protein was specifically located in nucleus. When transiently expressed in poplar leaf protoplasts, PdMYB118 specifically promoted the expression of anthocyanidin biosynthesis genes. Dual-luciferase assays revealed that PdMYB118 can directly activate the promoters of these genes. When overexpressed in Shanxin Yang (P. davidiana × P. bolleana), a hybrid clone commercially grown for landscaping in the northern part of China, transgenic plants overexpressing PdMYB118 produced more anthocyanins in the leaves and turned their color into redness when grown in both greenhouse and field. Consistently, transcripts of some important anthocyanidin biosynthesis genes were significantly increased in the leaves of transgenic plants. All these results indicate that PdMYB118 functions as an essential transcription factor regulating anthocyanin biosynthesis in poplar and could be used for the genetic engineering of colorful tree species.

MADS-box genes underground becoming mainstream: plant root developmental mechanisms

Plant growth is largely post-embryonic and depends on meristems that are active throughout the lifespan of an individual. Developmental patterns rely on the coordinated spatio-temporal expression of different genes, and the activity of transcription factors is particularly important during most morphogenetic processes. MADS-box genes constitute a transcription factor family in eukaryotes. In Arabidopsis, their proteins participate in all major aspects of shoot development, but their role in root development is still not well characterized. In this review we synthetize current knowledge pertaining to the function of MADS-box genes highly expressed in roots: XAL1, XAL2, ANR1 and AGL21, as well as available data for other MADS-box genes expressed in this organ. The role of Trithorax group and Polycomb group complexes on MADS-box genes' epigenetic regulation is also discussed. We argue that understanding the role of MADS-box genes in root development of species with contrasting architectures is still a challenge. Finally, we propose that MADS-box genes are key components of the gene regulatory networks that underlie various gene expression patterns, each one associated with the distinct developmental fates observed in the root. In the case of XAL1 and XAL2, their role within these networks could be mediated by regulatory feedbacks with auxin.

The R2R3 MYB transcription factor MYB189 negatively regulates secondary cell wall biosynthesis in Populus

Secondary cell wall (SCW) biosynthesis during wood formation in trees is controlled by a multilevel regulatory network that coordinates the expression of substantial genes. However, few transcription factors involved in the negative regulation of secondary wall biosynthesis have been characterized in tree species. In this study, we isolated an R2R3 MYB transcription factor MYB189 from Populus trichocarpa, which is expressed predominantly in secondary vascular tissues, especially in the xylem. A novel repression motif was identified in the C-terminal region of MYB189, which indicates this factor was a transcriptional repressor. Overexpression (OE) of MYB189 in Arabidopsis and poplar resulted in a significant reduction in the contents of lignin, cellulose and hemicelluloses. Vascular development in stems of MYB189 OE lines was markedly inhibited, leading to a dramatic decrease in SCW thickness of xylem fibers. Gene expression analyses showed that most of the structural genes involved in the biosynthesis of lignin, cellulose and xylans were significantly downregulated in MYB189-overexpressing poplars compared with the wild-type control. Chromatin immunoprecipitation-quantitative real-time polymerase chain reaction and transient expression assays revealed that MYB189 could directly bind to the promoters of secondary wall biosynthetic genes to repress their expression. Together, these data suggest that MYB189 acts as a repressor to regulate SCW biosynthesis in poplar.

Overexpression of CmMYB15 provides chrysanthemum resistance to aphids by regulating the biosynthesis of lignin

MYB transcription factors are widely involved in the development of and physiological processes in plants. Here, we isolated the chrysanthemum R2R3-MYB family transcription factor CmMYB15, a homologous gene of AtMYB15. It was demonstrated that CmMYB15 expression was induced by aphids and that CmMYB15 could bind to AC elements, which usually exist in the promoter of lignin biosynthesis genes. Overexpression of CmMYB15 in chrysanthemum enhanced the resistance of aphids. Additionally, the content of lignin and the expression of several lignin biosynthesis genes increased. In summary, the results indicate that CmMYB15 regulates lignin biosynthesis genes that enhance the resistance of chrysanthemum to aphids.

Genome-wide characterization of aspartic protease (AP) gene family in Populus trichocarpa and identification of the potential PtAPs involved in wood formation

BACKGROUND

Aspartic protease (AP) is one of four large proteolytic enzyme families that are involved in plant growth and development. Little is known about the AP gene family in tree species, although it has been characterized in Arabidopsis, rice and grape. The AP genes that are involved in tree wood formation remain to be determined.

RESULTS

A total of 67 AP genes were identified in Populus trichocarpa (PtAP) and classified into three categories (A, B and C). Chromosome mapping analysis revealed that two-thirds of the PtAP genes were located in genome duplication blocks, indicating the expansion of the AP family by segmental duplications in Populus. The microarray data from the Populus eFP browser demonstrated that PtAP genes had diversified tissue expression patterns. Semi-qRT-PCR analysis further determined that more than 10 PtAPs were highly or preferentially expressed in the developing xylem. When the involvement of the PtAPs in wood formation became the focus, many SCW-related cis-elements were found in the promoters of these PtAPs. Based on PtAP_promoter::GUS techniques, the activities of PtAP66 promoters were observed only in fiber cells, not in the vessels of stems as the xylem and leaf veins developed in the transgenic Populus tree, and strong GUS signals were detected in interfascicular fiber cells, roots, anthers and sepals of PtAP17_promoter::GUS transgenic plants. Intensive GUS activities in various secondary tissues implied that PtAP66 and PtAP17 could function in wood formation. In addition, most of the PtAP proteins were predicted to contain N- and (or) O-glycosylation sites, and the integration of PNGase F digestion and western blotting revealed that the PtAP17 and PtAP66 proteins were N-glycosylated in Populus.

CONCLUSIONS

Comprehensive characterization of the PtAP genes suggests their functional diversity during Populus growth and development. Our findings provide an overall understanding of the AP gene family in trees and establish a better foundation to further describe the roles of PtAPs in wood formation.

Genome-wide identification and expression analysis of R2R3, 3R- and 4R-MYB transcription factors during lignin biosynthesis in flax (Linum usitatissimum)

MYB transcription factors (TFs) have vital roles in regulating lignin or fiber development. Flax (Linum usitatissimum) is known as one of the plants with high fiber production capacity. However, no studies have been conducted to identify and characterize MYB TFs in the flax genome. Results showed that flax genome harbours 167 R2R3, seven 3R, and one 4R-type MYB TFs. 22 MYB genes (%13) were estimated to be tandem duplicated dated around 13.3-86.98 Mya. 130 flax MYB members have apparent orthologous with Arabidopsis, in which 17 R2R3 MYBs are associated with lignin biosynthesis. MYB062, MYB072, MYB096, MYB141, and MYB146 genes were up-regulated in tissues having higher lignin production capacity. In opposite, MYB012 and MYB113 genes were down-regulated which points out the involvement of those genes in the lignin biosynthesis mechanism. This comprehensive study can provide a basis for understanding the role of MYBs in fiber or lignin production in flax.

A novel MYB transcription factor regulates ascorbic acid synthesis and affects cold tolerance

Dehydroascorbate reductase (DHAR) plays an important role in stress responses, but the transcriptional regulation of DHAR in response to abiotic stress is still poorly understood. In this study, we isolated a novel R2R3-type MYB transcription factor from Pyrus betulaefolia by yeast one-hybrid screening, designated as PbrMYB5. PbrMYB5 was localized in the nucleus and could bind specifically to the promoter of PbrDHAR2. PbrMYB5 was greatly induced by cold and salt but slightly by dehydration. Overexpression of PbrMYB5 in tobacco conferred enhanced tolerance to chilling stresses, whereas down-regulation of PbrMYB5 in P. betulaefolia by virus-induced gene silencing resulted in elevated chilling sensitivity. Transgenic tobacco exhibited higher expression levels of NtDHAR2 and accumulated larger amount of ascorbic acid (AsA) than the wild-type plants. Virus-induced gene silencing of PbrMYB5 in P. betulaefolia down-regulated PbrDHAR2 abundance and decreased AsA level, accompanied by an increased sensitivity to the chilling stress. Taken together, these results demonstrated that PbrMYB5 was an activator of AsA biosynthesis and may play a positive role in chilling tolerance, at least in part, due to the modulation of AsA synthesis by regulating the PbrDHAR2 expression.

Aquaporins Respond to Chilling in the Phloem by Altering Protein and mRNA Expression

Previous experiments using heat exchangers (liquid cooled blocks) to chill a portion of plant stem have shown a transient stoppage in phloem translocation and an increase in measured phloem pressure. Although a chilled-induced stoppage of phloem transport has been known for over 100 years, the mechanism of this phenomenon is still poorly understood. Recently, work has highlighted that aquaporins occur within the plasma membrane of the sieve tubes along the entire source-to-sink pathway, and that isoforms of these water channel proteins may change dynamically. Aquaporins show regulatory roles in controlling tissue and cellular water status in response to environmental hardships. Thus, we tested if protein localization and mRNA transcript abundance changes occur in response to chilling in balsam poplar (Populus balsamifera) using immunohistochemistry and qrtPCR. The results of the immunolocalization experiments show that the labeling intensity of the sieve elements treated for only 2 min of chill time significantly increased for PIP2. After 10 min of chilling, this signal declined significantly to lower than that of the pre-chilled sieve elements. Overall, the abundance of mRNA transcript increased for the tested PIP2s following cold application. We discuss the implication that aquaporins are responsible for the alleviation of sieve tube pressure and the resumption of flow following a cold-induced blockage event.