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

Genome-Wide Characterization and Analysis of R2R3-MYB Genes Related to Fruit Ripening and Stress Response in Banana (Musa acuminata L. AAA Group, cv. 'Cavendish')

MYB is an important type of transcription factor in eukaryotes. It is widely involved in a variety of biological processes and plays a role in plant morphogenesis, growth and development, primary and secondary metabolite synthesis, and other life processes. In this study, bioinformatics methods were used to identify the R2R3-MYB transcription factor family members in the whole Musa acuminata (DH-Pahang) genome, one of the wild ancestors of banana. A total of 280 MaMYBs were obtained, and phylogenetic analysis indicated that these MaMYBs could be classified into 33 clades with MYBs from Arabidopsis thaliana. The amino acid sequences of the R2 and R3 Myb-DNA binding in all MaMYB protein sequences were quite conserved, especially Arg-12, Arg-13, Leu-23, and Leu-79. Distribution mapping results showed that 277 MaMYBs were localized on the 11 chromosomes in the Musa acuminata genome. The MaMYBs were distributed unevenly across the 11 chromosomes. More than 40.0% of the MaMYBs were located in collinear fragments, and segmental duplications likely played a key role in the expansion of the MaMYBs. Moreover, the expression profiles of MaMYBs in different fruit development and ripening stages and under various abiotic and biotic stresses were investigated using available RNA-sequencing data to obtain fruit development, ripening-specific, and stress-responsive candidate genes. Weighted gene co-expression network analysis (WGCNA) was used to analyze transcriptome data of banana from the above 11 samples. We found MaMYBs participating in important metabolic biosynthesis pathways in banana. Collectively, our results represent a comprehensive genome-wide study of the MaMYB gene family, which should be helpful in further detailed studies on MaMYBs functions related to fruit development, postharvest ripening, and the seedling response to stress in an important banana cultivar.

A Novel R2R3-MYB Transcription Factor SbMYB12 Positively Regulates Baicalin Biosynthesis in Scutellaria baicalensis Georgi

Scutellaria baicalensis Georgi is an annual herb from the Scutellaria genus that has been extensively used as a traditional medicine for over 2000 years in China. Baicalin and other flavonoids have been identified as the principal bioactive ingredients. The biosynthetic pathway of baicalin in S. baicalensis has been elucidated; however, the specific functions of R2R3-MYB TF, which regulates baicalin synthesis, has not been well characterized in S. baicalensis to date. Here, a S20 R2R3-MYB TF (SbMYB12), which encodes 263 amino acids with a length of 792 bp, was expressed in all tested tissues (mainly in leaves) and responded to exogenous hormone methyl jasmonate (MeJA) treatment. The overexpression of SbMYB12 significantly promoted the accumulation of flavonoids such as baicalin and wogonoside in S. baicalensis hairy roots. Furthermore, biochemical experiments revealed that SbMYB12 is a nuclear-localized transcription activator that binds to the SbCCL7-4, SbCHI-2, and SbF6H-1 promoters to activate their expression. These results illustrate that SbMYB12 positively regulates the generation of baicalin and wogonoside. In summary, this work revealed a novel S20 R2R3-MYB regulator and enhances our understanding of the transcriptional and regulatory mechanisms of baicalin biosynthesis, as well as sheds new light on metabolic engineering in S. baicalensis.

Investigation and Expression Analysis of R2R3-MYBs and Anthocyanin Biosynthesis-Related Genes during Seed Color Development of Common Bean (Phaseolus vulgaris)

Anthocyanins are responsible for the coloration of common bean seeds, and their accumulation is positively correlated with the expression level of anthocyanin biosynthetic genes. The MBW (MYB-bHLH-WD40) complex is thought to regulate the expression of these genes, and MYB proteins, which are a key factor in activating anthocyanin pathway genes, have been identified in several plants. This study demonstrated gene structures, chromosomal placements, gene duplications of R2R3-MYBs, miRNAs associated with R2R3-MYBs, and the interaction of these genes with other flavonoid regulatory genes. qRT-PCR was used to investigate the role of specific R2R3-MYBs and flavonoid genes in common bean seed color development. As a result of a comprehensive analysis with the help of in silico tools, we identified 160 R2R3-MYB genes in the common bean genome. We divided these genes into 16 classes on the basis of their intron-exon and motif structures. Except for three, the rest of the common bean R2R3-MYB members were distributed to all chromosomes with different densities, primarily located on chromosomes 3 and 8. We identified a total of 44 duplicated gene pairs dispersed across 11 chromosomes and evolved under purifying selection (Ka/Ks  <  1), 19 of which were derived from a whole-genome duplication. Our research uncovered 25 putative repressor PvMYB proteins that contain the EAR motif. Additionally, fifty different cis-regulatory elements regulated by light, stress, and hormone were identified. Within the genome of the common bean, we discovered a total of 36 microRNAs that target a total of 72 R2R3-MYB transcripts. The effect of 16 R2R3-MYB genes and 16 phenylpropanoid pathway genes, selected on the basis of their interaction in the protein-protein interaction map, playing role in the regulation of seed coat color development was evaluated using qRT-PCR in 5 different tissues at different developmental stages. The results revealed that these specific genes have different expression levels during different developmental periods, with higher levels in the pod filling and early pod stages than in the rest of the developmental periods. Furthermore, it was shown that PvTT8 (bHLH), PvTT2 (PvMYB42), PvMYB113, PvTTG1, and PvWD68 genes have effects on the regulation of seed coat color. The findings of this study, which is the first to use whole-genome analysis to identify and characterize the R2R3-MYB genes in common bean, may serve as a reference for future functional research in the legume.

MYB transcription factors-master regulators of phenylpropanoid biosynthesis and diverse developmental and stress responses

Phenylpropanoids, the largest class of natural products including flavonoids, anthocyanins, monolignols and tannins perform multiple functions ranging from photosynthesis, nutrient uptake, regulating growth, cell division, maintenance of redox homeostasis and biotic and abiotic stress responses. Being sedentary life forms, plants possess several regulatory modules that increase their performance in varying environments by facilitating activation of several signaling cascades upon perception of developmental and stress signals. Of the various regulatory modules, those involving MYB transcription factors are one of the extensive groups involved in regulating the phenylpropanoid metabolic enzymes in addition to other genes. R2R3 MYB transcription factors are a class of plant-specific transcription factors that regulate the expression of structural genes involved in anthocyanin, flavonoid and monolignol biosynthesis which are indispensable to several developmental pathways and stress responses. The aim of this review is to present the regulation of the phenylpropanoid pathway by MYB transcription factors via Phospholipase D/phosphatidic acid signaling, downstream activation of the structural genes, leading to developmental and/or stress responses. Specific MYB transcription factors inducing or repressing specific structural genes of anthocyanin, flavonoid and lignin biosynthetic pathways are discussed. Further the roles of MYB in activating biotic and abiotic stress responses are delineated. While several articles have reported the role of MYB's in stress responses, they are restricted to two or three specific MYB factors. This review is a consolidation of the diverse roles of different MYB transcription factors involved both in induction and repression of anthocyanin, flavonoid, and lignin biosynthesis.

Genome-Scale Identification, Classification, and Expression Profiling of MYB Transcription Factor Genes in Cinnamomum camphora

The camphor tree (Cinnamomum camphora (L.) Presl.) is the representative species of subtropical evergreen broadleaved forests in eastern Asia and an important raw material for essential oil production worldwide. Although MYBs have been comprehensively characterized and their functions have been partially resolved in many plants, it has not been explored in C. camphora. In this study, 121 CcMYBs were identified on 12 chromosomes in the whole genome of C. camphora and found that CcMYBs were mainly expanded by segmental duplication. They were divided into 28 subgroups based on phylogenetic analysis and gene structural characteristics. In the promoter regions, numerous cis-acting elements were related to biological processes. Analysis of RNA sequencing data from seven tissues showed that CcMYBs exhibited different expression profiles, suggesting that they have various roles in camphor tree development. In addition, combined with the correlation analysis of structural genes in the flavonoid synthesis pathway, we identified CcMYBs from three subgroups that might be related to the flavonoid biosynthesis pathway. This study systematically analyzed CcMYBs in C. camphora, which will set the stage for subsequent research on the functions of CcMYBs during their lifetime and provide valuable insights for the genetic improvement of camphor trees.

Genome-wide analysis of R2R3-MYB transcription factors in Japanese morning glory

The R2R3-MYB transcription factor is one of the largest transcription factor families in plants. R2R3-MYBs play a variety of functions in plants, such as cell fate determination, organ and tissue differentiations, primary and secondary metabolisms, stress and defense responses and other physiological processes. The Japanese morning glory (Ipomoea nil) has been widely used as a model plant for flowering and morphological studies. In the present study, 127 R2R3-MYB genes were identified in the Japanese morning glory genome. Information, including gene structure, protein motif, chromosomal location and gene expression, were assigned to the InR2R3-MYBs. Phylogenetic tree analysis revealed that the 127 InR2R3-MYBs were classified into 29 subfamilies (C1-C29). Herein, physiological functions of the InR2R3-MYBs are discussed based on the functions of their Arabidopsis orthologues. InR2R3-MYBs in C9, C15, C16 or C28 may regulate cell division, flavonol biosynthesis, anthocyanin biosynthesis or response to abiotic stress, respectively. C16 harbors the known anthocyanin biosynthesis regulator, InMYB1 (INIL00g10723), and putative anthocyanin biosynthesis regulators, InMYB2 (INIL05g09650) and InMYB3 (INIL05g09651). In addition, INIL05g09649, INIL11g40874 and INIL11g40875 in C16 were suggested as novel anthocyanin biosynthesis regulators. We organized the R2R3-MYB transcription factors in the morning glory genome and assigned information to gene and protein structures and presuming their functions. Our study is expected to facilitate future research on R2R3-MYB transcription factors in Japanese morning glory.

An R2R3-MYB transcription factor VyMYB24, isolated from wild grape Vitis yanshanesis J. X. Chen., regulates the plant development and confers the tolerance to drought

In grapevines, the MYB transcription factors play an important regulatory role in the phenylpropanoid pathway including proanthocyanidin, anthocyanin, and flavonoid biosynthesis. However, the role of MYB in abiotic stresses is not clear. In this study, an R2R3-MYB transcription factor, VyMYB24, was isolated from a high drought-tolerant Chinese wild Vitis species V. yanshanesis. Our findings demonstrated that it was involved in plant development and drought tolerance. VyMYB24 is a nuclear protein and is significantly induced by drought stress. When over-expressed in tobacco, VyMYB24 caused plant dwarfing including plant height, leaf area, flower size, and seed weight. The GA1+3 content in transgenic plants was reduced significantly, and spraying exogenous gibberellin could recover the dwarf phenotype of VyMYB24 transgenic plants, suggesting that VyMYB24 might inhibit plant development by the regulation of gibberellin (GA) metabolism. Under drought stress, the VyMYB24 transgenic plants improved their tolerance to drought with a lower wilting rate, lower relative electrical conductivity, and stronger roots. Compared to wild-type tobacco plants, VyMYB24 transgenic plants accumulated less reactive oxygen, accompanied by increased antioxidant enzyme activity and upregulated gene expression levels of superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) genes. In addition, transgenic plants accumulated more proline, and their related synthetic genes NtP5CR and NtP5CS genes were significantly upregulated when exposed to drought. Besides, abiotic stress-responsive genes, NtDREB, NtERD10C, NtERD10D, and NtLEA5, were upregulated significantly in VyMYB24 transgenic plants. These results indicate that VyMYB24 plays a positive regulatory role in response to drought stress and also regulates plant development, which provides new evidence to further explore the molecular mechanism of drought stress of the MYB gene family.

The R2R3-MYB gene family in Cicer arietinum: genome-wide identification and expression analysis leads to functional characterization of proanthocyanidin biosynthesis regulators in the seed coat

We identified 119 typical CaMYB encoding genes and reveal the major components of the proanthocyanidin regulatory network. CaPARs emerged as promising targets for genetic engineering toward improved agronomic traits in C. arietinum. Chickpea (Cicer arietinum) is among the eight oldest crops and has two main types, i.e., desi and kabuli, whose most obvious difference is the color of their seeds. We show that this color difference is due to differences in proanthocyanidin content of seed coats. Using a targeted approach, we performed in silico analysis, metabolite profiling, molecular, genetic, and biochemical studies to decipher the transcriptional regulatory network involved in proanthocyanidin biosynthesis in the seed coat of C. arietinum. Based on the annotated C. arietinum reference genome sequence, we identified 119 typical CaMYB encoding genes, grouped in 32 distinct clades. Two CaR2R3-MYB transcription factors, named CaPAR1 and CaPAR2, clustering with known proanthocyanidin regulators (PARs) were identified and further analyzed. The expression of CaPAR genes correlated well with the expression of the key structural proanthocyanidin biosynthesis genes CaANR and CaLAR and with proanthocyanidin levels. Protein-protein interaction studies suggest the in vivo interaction of CaPAR1 and CaPAR2 with the bHLH-type transcription factor CaTT8. Co-transfection analyses using Arabidopsis thaliana protoplasts showed that the CaPAR proteins form a MBW complex with CaTT8 and CaTTG1, able to activate the promoters of CaANR and CaLAR in planta. Finally, transgenic expression of CaPARs in the proanthocyanidin-deficient A. thaliana mutant tt2-1 leads to complementation of the transparent testa phenotype. Taken together, our results reveal main components of the proanthocyanidin regulatory network in C. arietinum and suggest that CaPARs are relevant targets of genetic engineering toward improved agronomic traits.

MicroRNA828 negatively regulates lignin biosynthesis in stem of Populus tomentosa through MYB targets

Lignin biosynthesis in the sclerenchyma cells is strictly controlled by a complex network of genetic and environmental signals. In the last decades, the transcriptional regulation of lignin synthesis in woody species has been established. However, the role of microRNA-mediated post-transcriptional modulation in secondary cell wall biosynthesis remains poorly understood. Here, we identified a microRNA, miR828, involved in the regulation specific to lignin biosynthesis during stem development in Populus tomentosa Carr. miR828 is preferentially expressed in the secondary vascular tissues during stem development. Two MYB genes (MYB171 and MYB011) were validated as direct targets of miR828 by degradome analysis and green fluorescent protein signal detection. Overexpression of miR828 in poplar downregulated genes for lignin biosynthesis, resulting in reduced lignin content in cell walls. Conversely, suppression of miR828 in plants by the short tandem target mimics elevated the expression of lignin biosynthetic genes and increased lignin deposition. We further revealed that poplar MYB171, as the most abundant miR828 target in the stem, is a positive regulator for lignin biosynthesis. Transient expression assays showed that both MYB171 and MYB011 activated PAL1 and CCR2 transcription, whereas the introduction of miR828 significantly suppressed their expression that was induced by MYB171 or MYB011. Collectively, our results demonstrate that the miR828-MYBs module precisely regulates lignin biosynthesis during the stem development in P. tomentosa through transcriptional and post-transcriptional manners.

BcWRKY33A Enhances Resistance to Botrytis cinerea via Activating BcMYB51-3 in Non-Heading Chinese Cabbage

The transcription factor WRKY33 is a vital regulator of the biological process of the necrotrophic fungus Botrytis cinerea (B. cinerea). However, its specific regulatory mechanism remains to be further investigated. In non-heading Chinese cabbage (NHCC, Brassica campestris (syn. Brassica rapa) ssp. Chinensis), our previous study showed that BcWRKY33A is induced not only by salt stress, but also by B. cinerea infection. Here, we noticed that BcWRKY33A is expressed in trichomes and confer plant defense resistance. Disease symptoms and qRT-PCR analyses revealed that BcWRKY33A-overexpressing and -silencing lines were less and more severely impaired, respectively, than wild type upon B. cinerea treatment. Meanwhile, the transcripts’ abundance of indolic glucosinolates’ (IGSs) biosynthetic genes is consistent with plants’ B. cinerea tolerance. Identification and expression pattern analysis of BcMYB51s showed that BcMYB51-3 has a similar trend to BcWRKY33A upon B. cinerea infection. Moreover, BcWRKY33A directly binds to the BcMYB51-3 promoter, which was jointly confirmed by Y1H, dual-LUC, and EMSA assays. The importance of MYB51, the homolog of BcMYB51-3, in the BcWRKY33A-mediated B. cinerea resistance was also verified using the TRV-based VIGS system. Overall, our data concludes that BcWRKY33A directly activates the expression of BcMYB51-3 and downstream IGSs’ biosynthetic genes, thereby improving the B. cinerea tolerance of NHCC plants.

R2R3-MYB gene family: Genome-wide identification provides insight to improve the content of proanthocyanidins in Trifolium repens

The R2R3-MYB family is one of largest transcription factor families in plants playing significant roles in regulating anthocyanin and proanthocyanidin biosynthesis. Proanthocyanidins are one of major objectives to improve the quality of white clover (Trifolium repens L.), which have a beneficial effect on ruminant to prevent the lethal pasture bloat. A total of 133 TrR2R3-MYB genes were identified and distributed on all 16 chromosomes based on the whole genome information of white clover. Also, by exploring the gene structure, motifs and duplication events of TrR2R3-MYBs, as well as the evolutionary relationship with TrR2R3-MYB genes of other species, 10 TrR2R3-MYB genes with the potential to regulate the anthocyanins and proanthocyanidins biosynthesis were screened. These TrR2R3-MYB genes responded significantly to low temperature in white clover. In addition, they have different expression patterns in leaves, petioles and inflorescences of white clover. Importantly, TrMYB116 and TrMYB118 may positively regulate anthocyanin accumulation and low temperature response in white clover. TrMYB118 may also be associated with anthocyanin pigmentation pattern in Purple leaves. This study provides a basis for verifying the function of TrR2R3-MYB and breeding white clover cultivars with high proanthocyanidins.

Assessing the Influencing Factors of Electronic Word-of-Mouth on CSR Issues. A Case of Hospotality Service Industry of China

Corporate social responsibility (CSR) information can be effectively disseminated via social media in a variety of industries, including the hospitality sector. In the same way, the media has a significant impact on CSR because the news media helps companies achieve their CSR goals. Prior research has not examined the main factors that influence electronic word of mouth (eWOM) on media coverage of CSR issues via social networking websites. For the purpose of examining the most significant predictors of intention to share or comment on negative CSR news reported by one media outlet on a specific social networking site (SNS). 677 Wechat users in China were surveyed in order to test the proposed model empirically. According to the findings of the study, eWOM intentions are positively influenced by environmental CSR content, and advertisement related CSR content. It also confirmed that the value of information is positively influenced by the credibility of the source. The variables interpersonal influencer impact and privacy concerns had no significant relationship, nor did they have any significant relationship with the intentions to share and comment on Wechat. Further the study findings suggest the theoretical and managerial policy recommendation for decision makers.

Genome-wide identification and expression analysis of the R2R3-MYB gene family in tobacco (Nicotiana tabacum L.)

BACKGROUND

The R2R3-MYB transcription factor is one of the largest gene families in plants and involved in the regulation of plant development, hormone signal transduction, biotic and abiotic stresses. Tobacco is one of the most important model plants. Therefore, it will be of great significance to investigate the R2R3-MYB gene family and their expression patterns under abiotic stress and senescence in tobacco.

RESULTS

A total of 174 R2R3-MYB genes were identified from tobacco (Nicotiana tabacum L.) genome and were divided into 24 subgroups based on phylogenetic analysis. Gene structure (exon/intron) and protein motifs were especially conserved among the NtR2R3-MYB genes, especially members within the same subgroup. The NtR2R3-MYB genes were distributed on 24 tobacco chromosomes. Analysis of gene duplication events obtained 3 pairs of tandem duplication genes and 62 pairs of segmental duplication genes, suggesting that segmental duplications is the major pattern for R2R3-MYB gene family expansion in tobacco. Cis-regulatory elements of the NtR2R3-MYB promoters were involved in cellular development, phytohormones, environmental stress and photoresponsive. Expression profile analysis showed that NtR2R3-MYB genes were widely expressed in different maturity tobacco leaves, and however, the expression patterns of different members appeared to be diverse. The qRT-PCR analysis of 15 NtR2R3-MYBs confirmed their differential expression under different abiotic stresses (cold, salt and drought), and notably, NtMYB46 was significantly up-regulated under three treatments.

CONCLUSIONS

In summary, a genome-wide identification, evolutionary and expression analysis of R2R3-MYB gene family in tobacco were conducted. Our results provided a solid foundation for further biological functional study of NtR2R3-MYB genes in tobacco.

In silico genome wide identification and expression analysis of the WUSCHEL-related homeobox gene family in Medicago sativa

Alfalfa (Medicago sativa) is an important food and feed crop which rich in mineral sources. The WUSCHEL-related homeobox (WOX) gene family plays important roles in plant development and identification of putative gene families, their structure, and potential functions is a primary step for not only understanding the genetic mechanisms behind various biological process but also for genetic improvement. A variety of computational tools, including MAFFT, HMMER, hidden Markov models, Pfam, SMART, MEGA, ProtTest, BLASTn, and BRAD, among others, were used. We identified 34 MsWOX genes based on a systematic analysis of the alfalfa plant genome spread in eight chromosomes. This is an expansion of the gene family which we attribute to observed chromosomal duplications. Sequence alignment analysis revealed 61 conserved proteins containing a homeodomain. Phylogenetic study sung reveal five evolutionary clades with 15 motif distributions. Gene structure analysis reveals various exon, intron, and untranslated structures which are consistent in genes from similar clades. Functional analysis prediction of promoter regions reveals various transcription binding sites containing key growth, development, and stress-responsive transcription factor families such as MYB, ERF, AP2, and NAC which are spread across the genes. Most of the genes are predicted to be in the nucleus. Also, there are duplication events in some genes which explain the expansion of the family. The present research provides a clue on the potential roles of MsWOX family genes that will be useful for further understanding their functional roles in alfalfa plants.

Multiplex knockout of trichome-regulating MYB duplicates in hybrid poplar using a single gRNA

As the focus for CRISPR/Cas-edited plants moves from proof-of-concept to real-world applications, precise gene manipulation will increasingly require concurrent multiplex editing for polygenic traits. A common approach for editing across multiple sites is to design one guide RNA (gRNA) per target; however, this complicates construct assembly and increases the possibility of off-target mutations. In this study, we utilized one gRNA to target MYB186, a known positive trichome regulator, as well as its paralogs MYB138 and MYB38 at a consensus site for mutagenesis in hybrid poplar (Populus tremula × P. alba INRA 717-1B4). Unexpected duplications of MYB186 and MYB138 resulted in eight alleles for the three targeted genes in the hybrid poplar. Deep sequencing and polymerase chain reaction analyses confirmed editing across all eight targets in nearly all of the resultant glabrous mutants, ranging from small indels to large genomic dropouts, with no off-target activity detected at four potential sites. This highlights the effectiveness of a single gRNA targeting conserved exonic regions for multiplex editing. Additionally, cuticular wax and whole-leaf analyses showed a complete absence of triterpenes in the trichomeless mutants, hinting at a previously undescribed role for the nonglandular trichomes of poplar.

Diverse roles of MYB transcription factors in regulating secondary metabolite biosynthesis, shoot development, and stress responses in tea plants (Camellia sinensis)

Tea (Camellia sinensis) is concocted from tea plant shoot tips that produce catechins, caffeine, theanine, and terpenoids, which collectively determine the rich flavors and health benefits of the infusion. However, little is known about the integrated regulation of shoot tip development and characteristic secondary metabolite biosynthesis in tea plants. Here, we demonstrate that MYB transcription factors (TFs) play key and yet diverse roles in regulating leaf and stem development, secondary metabolite biosynthesis, and environmental stress responses in tea plants. By integrating transcriptomic and metabolic profiling data in different tissues at a series of developmental stages or under various stress conditions, alongside biochemical and genetic analyses, we predicted the MYB TFs involved in regulating shoot development (CsMYB2, 98, 107, and 221), epidermal cell initiation (CsMYB184, 41, 139, and 219), stomatal initiation (CsMYB113 and 153), and the biosynthesis of flavonoids (including catechins, anthocyanins, and flavonols; CsMYB8 and 99), caffeine (CsMYB85 and 86), theanine (CsMYB9 and 49), carotenoids (CsMYB110), mono-/sesquiterpenoid volatiles (CsMYB68, 147, 148, and 193), lignin (CsMYB164 and 192), and indolic compounds (CsMYB139, 162, and 198), as well as the MYB TFs that are likely involved in hormone signaling-mediated environmental stress and defense responses. We characterized the functions of some key MYBs in regulating flavonoid and carotenoid biosynthesis for tea quality and flavor. This study provides a cross-family analysis of MYBs in tea alongside new insights into the coordinated regulation of tea plant shoot development and secondary metabolism, paving the way towards understanding of tea quality trait formation and genetic improvement of quality tea plants.

Dimerization of PtrMYB074 and PtrWRKY19 mediates transcriptional activation of PtrbHLH186 for secondary xylem development in Populus trichocarpa

Wood formation is controlled by transcriptional regulatory networks (TRNs) involving regulatory homeostasis determined by combinations of transcription factor (TF)-DNA and TF-TF interactions. Functions of TF-TF interactions in wood formation are still in the early stages of identification. PtrMYB074 is a woody dicot-specific TF in a TRN for wood formation in Populus trichocarpa. Here, using yeast two-hybrid and bimolecular fluorescence complementation, we conducted a genome-wide screening for PtrMYB074 interactors and identified 54 PtrMYB074-TF pairs. Of these pairs, 53 are novel. We focused on the PtrMYB074-PtrWRKY19 pair, the most highly expressed and xylem-specific interactor, and its direct transregulatory target, PtrbHLH186, the xylem-specific one of the pair's only two direct TF target genes. Using transient and CRISPR-mediated transgenesis in P. trichocarpa coupled with chromatin immunoprecipitation and electrophoretic mobility shift assays, we demonstrated that PtrMYB074 is recruited by PtrWRKY19 and that the PtrMYB074-PtrWRKY19 dimers are required to transactive PtrbHLH186. Overexpressing PtrbHLH186 in P. trichocarpa resulted in retarded plant growth, increased guaiacyl lignin, a higher proportion of smaller stem vessels and strong drought-tolerant phenotypes. Knowledge of the PtrMYB074-PtrWRKY19-PtrbHLH186 regulation may help design genetic controls of optimal growth and wood formation to maximize beneficial wood properties while minimizing negative effects on growth.

Identification and Characterization of SOG1 (Suppressor of Gamma Response 1) Homologues in Plants Using Data Mining Resources and Gene Expression Profiling

SOG1 (Suppressor of the Gamma response 1) is the master-regulator of plant DNA damage response (DDR), a highly coordinated network of DNA damage sensors, transducers, mediators, and effectors, with highly coordinated activities. SOG1 transcription factor belongs to the NAC/NAM protein family, containing the well-conserved NAC domain and five serine-glutamine (SQ) motifs, preferential targets for phosphorylation by ATM and ATR. So far, the information gathered for the SOG1 function comes from studies on the model plant Arabidopsis thaliana. To expand the knowledge on plant-specific DDR, it is opportune to gather information on other SOG1 orthologues. The current study identified plants where multiple SOG1 homologues are present and evaluated their functions by leveraging the information contained in publicly available transcriptomics databases. This analysis revealed the presence of multiple SOG1 sequences in thirteen plant species, and four (Medicago truncatula, Glycine max, Kalankoe fedtschenkoi, Populus trichocarpa) were selected for gene expression data mining based on database availability. Additionally, M. truncatula seeds and seedlings exposed to treatments known to activate DDR pathways were used to evaluate the expression profiles of MtSOG1a and MtSOG1b. The experimental workflow confirmed the data retrieved from transcriptomics datasets, suggesting that the SOG1 homologues have redundant functions in different plant species.

Identification of the Wheat (Triticum aestivum) IQD Gene Family and an Expression Analysis of Candidate Genes Associated with Seed Dormancy and Germination

The IQ67 Domain (IQD) gene family plays important roles in plant developmental processes and stress responses. Although IQDs have been characterized in model plants, little is known about their functions in wheat (Triticum aestivum), especially their roles in the regulation of seed dormancy and germination. Here, we identified 73 members of the IQD gene family from the wheat genome and phylogenetically separated them into six major groups. Gene structure and conserved domain analyses suggested that most members of each group had similar structures. A chromosome positional analysis showed that TaIQDs were unevenly located on 18 wheat chromosomes. A synteny analysis indicated that segmental duplications played significant roles in TaIQD expansion, and that the IQD gene family underwent strong purifying selection during evolution. Furthermore, a large number of hormone, light, and abiotic stress response elements were discovered in the promoters of TaIQDs, implying their functional diversity. Microarray data for 50 TaIQDs showed different expression levels in 13 wheat tissues. Transcriptome data and a quantitative real-time PCR analysis of wheat varieties with contrasting seed dormancy and germination phenotypes further revealed that seven genes (TaIQD4/-28/-32/-58/-64/-69/-71) likely participated in seed dormancy and germination through the abscisic acid-signaling pathway. The study results provide valuable information for cloning and a functional investigation of candidate genes controlling wheat seed dormancy and germination; consequently, they increase our understanding of the complex regulatory networks affecting these two traits.

Identification and Expression Analysis of R2R3-MYB Family Genes Associated with Salt Tolerance in Cyclocarya paliurus

R2R3-MYB transcription factors are most abundant in the MYB superfamily, while the R2R3-MYB genes play an important role in plant growth and development, especially in response to environmental stress. Cyclocarya paliurus is a multifunction tree species, and the existing resources cannot meet the requirement for its leaf production and medical use. Therefore, lands with some environmental stresses would be potential sites for developing C. paliurus plantations. However, the function of R2R3-MYB genes in C.paliurus in response to environmental stress remains unknown. In this study, to identify the roles of R2R3-MYB genes associated with salt stress response, 153 CpaMYB genes and their corresponding protein sequences were identified from the full-length transcriptome. Based on the comparison with MYB protein sequences of Arabidopsis thaliana, 69 R2R3-MYB proteins in C. paliurus were extracted for further screening combined with conserved functional domains. Furthermore, the MYB family members were analyzed from the aspects of protein sequences alignment, evolution, motif prediction, promoter cis-acting element analysis, and gene differential expression under different salt treatments using both a pot experiment and hydroponic experiment. The results showed that the R2R3-MYB genes of C.paliurus conserved functional domains, whereas four R2R3-MYB genes that might respond to salt stress via regulating plant hormone signals were identified in this study. This work provides a basis for further functional characterization of R2R3-MYB TFs in C. paliurus.

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.

Characterization of the wheat VQ protein family and expression of candidate genes associated with seed dormancy and germination

BACKGROUND

Seed dormancy and germination determine wheat resistance to pre-harvest sprouting and thereby affect grain yield and quality. Arabidopsis VQ genes have been shown to influence seed germination; however, the functions of wheat VQ genes have not been characterized.

RESULTS

We identified 65 TaVQ genes in common wheat and named them TaVQ1-65. We identified 48 paralogous pairs, 37 of which had Ka/Ks values greater than 1, suggesting that most TaVQ genes have experienced positive selection. Chromosome locations, gene structures, promoter element analysis, and gene ontology annotations of the TaVQs showed that their structures determined their functions and that structural changes reflected functional diversity. Transcriptome-based expression analysis of 62 TaVQ genes and microarray analysis of 11 TaVQ genes indicated that they played important roles in diverse biological processes. We compared TaVQ gene expression and seed germination index values among wheat varieties with contrasting seed dormancy and germination phenotypes and identified 21 TaVQ genes that may be involved in seed dormancy and germination.

CONCLUSIONS

Sixty-five TaVQ proteins were identified for the first time in common wheat, and bioinformatics analyses were used to investigate their phylogenetic relationships and evolutionary divergence. qRT-PCR data showed that 21 TaVQ candidate genes were potentially involved in seed dormancy and germination. These findings provide useful information for further cloning and functional analysis of TaVQ genes and introduce useful candidate genes for the improvement of PHS resistance in wheat.

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.

TcMYB29a, an ABA-Responsive R2R3-MYB Transcriptional Factor, Upregulates Taxol Biosynthesis in Taxus chinensis

Paclitaxel (Taxol), a highly modified diterpene agent mainly obtained from Taxus species, is the most widely used anticancer drug. Abscisic acid (ABA) is a well-known stress hormone that plays important roles in the secondary metabolism of plants, and it can also induce the accumulation of taxol in Taxus cell suspension cultures. However, the mechanism behind the regulation of taxol biosynthesis by ABA remains largely unknown. In previous research, a R2R3 MYB transcription factor (TF) TcMYB29a was observed to show a significant correlation with taxol biosynthesis, indicative of its potential role in the taxol biosynthesis. In this study, the TcMYB29a encoded by its gene was further characterized. An expression pattern analysis revealed that TcMYB29a was highly expressed in the needles and roots. Overexpression of TcMYB29a in Taxus chinensis cell suspension cultures led to an increased accumulation of taxol, and upregulated expression of taxol-biosynthesis-related genes, including the taxadiene synthase (TS) gene, the taxane 5α-hydroxylase (T5OH) gene, and the 3'-N-debenzoyl-2'-deoxytaxol-N-benzoyltransferase (DBTNBT) gene as compared to the controls. Chromatin immunoprecipitation (ChIP) assays, yeast one-hybrid (Y1H) assays, electrophoretic mobility shift assays (EMSAs), and dual-luciferase reporter assays verified that TcMYB29a could bind and activate the promoter of TcT5OH. Promoter sequence analysis of TcMYB29a revealed that its promoter containing an AERB site from -313 to -319 was a crucial ABA-responsive element. Subsequently, the ABA treatment assay showed that TcMYB29a was strongly upregulated at 6 h after ABA pretreatment. Furthermore, TcMYB29a was strongly suppressed at 3 h after the methyl jasmonate (MeJA) treatment and was depressed to the platform at 12 h. Taken together, these results reveal that TcMYB29a is an activator that improves the accumulation of taxol in Taxus chinensis cells through an ABA-medicated signaling pathway which is different from JA-medicated signaling pathways for the accumulation of taxol. These findings provide new insights into the potential regulatory roles of MYBs on the expression of taxol biosynthetic genes in Taxus.

Identification of Malus halliana R2R3-MYB gene family under iron deficiency stress and functional characteristics of MhR2R3-MYB4 in Arabidopsis thaliana

Iron (Fe) is an essential element for plant growth and development. Fe deficiency can trigger leaf chlorosis and reduce fruit yield. Therefore, it is necessary to explore transcription factors in response to Fe deficiency stress. A total of 29 MhR2R3-MYB transcription factors were identified based on the transcriptome of Malus halliana under Fe deficiency stress. A comprehensive analysis of physical and chemical properties, gene structures, conserved motif composition, evolutionary relationship and chromosome distribution was performed. Subsequently, based on the transcriptome, 14 genes with the most significant expression under Fe deficiency stress were screened for qRT-PCR verification. Among them,the functional characteristics of MhR2R3-MYB4 (MD05G1089600) were further studied in Arabidopsis thaliana. Expression of 13 out of these 14 genes was upregulated, only one was downregulated. Maximum upregulation of MhR2R3-MYB4 under Fe deficiency was 36.39-fold and 58.21-fold compared with day 0 in leaves and roots, respectively. Overexpression of MhR2R3-MYB4 enhanced tolerance to Fe deficiency in A. thaliana and led to multiple biochemical changes: transgenic lines have higher chl a, chl b and Fe²⁺ content, higher enzyme activity (SOD, POD, CAT and FCR) and lower chlorosis than the wild type in Fe deficiency conditions. We suggest that MhR2R3-MYB4 plays an important part in Fe deficiency stress, which may contribute to improve Fe deficiency tolerance of apple in future.

Characterization and Interaction Analysis of the Secondary Cell Wall Synthesis-Related Transcription Factor PmMYB7 in Pinus massoniana Lamb

In vascular plants, the importance of R2R3-myeloblastosis (R2R3-MYB) transcription factors (TFs) in the formation of secondary cell walls (SCWs) has long been a controversial topic due to the lack of empirical evidence of an association between TFs and downstream target genes. Here, we found that the transcription factor PmMYB7, which belongs to the R2R3-MYB subfamily, is involved in lignin biosynthesis in Pinus massoniana. PmMYB7 was highly expressed in lignified tissues and upon abiotic stress. As a bait carrier, the PmMYB7 protein had no toxicity or autoactivation in the nucleus. Forty-seven proteins were screened from the P. massoniana yeast library. These proteins were predicted to be mainly involved in resistance, abiotic stress, cell wall biosynthesis, and cell development. We found that the PmMYB7 protein interacted with caffeoyl CoA 3-O-methyltransferase-2 (PmCCoAOMT2)—which is involved in lignin biosynthesis—but not with beta-1, 2-xylosyltransferase (PmXYXT1) yeast two-hybrid (Y2H) studies. Our in vivo coimmunoprecipitation (Co-IP) assay further showed that the PmMYB7 and PmCCoAOMT2 proteins could interact. Therefore, we concluded that PmMYB7 is an upstream TF that can interact with PmCCoAOMT2 in plant cells. These findings lay a foundation for further research on the function of PmMYB7, lignin biosynthesis and molecular breeding in P. massoniana.

Ectopic expression of Populus MYB10 promotes secondary cell wall thickening and inhibits anthocyanin accumulation

Secondary cell wall (SCW) formation is regulated by a multilevel transcriptional regulatory network, in which MYB transcription factors (TFs) play key roles. In woody plants, hundreds of MYB TFs have been identified, most of which have unknown functions in wood SCW biosynthesis. Here, we characterized the function of a Populus MYB gene, PtoMYB10. PtoMYB10 was found to encode an R2R3-MYB TF and exhibit dominant expression in xylem tissues. PtoMYB10 was determined to be located in the nucleus with the ability to activate transcription. Overexpression of PtoMYB10 in Populus resulted in a drastic increase in SCW thickening in xylem fiber cells as well as ectopic deposition of lignin in cortex cells. The expression of genes associated with lignin biosynthesis was induced in PtoMYB10 overexpressing plants, whereas repressed gene expression was found with the anthocyanin biosynthesis pathway. Lignin and anthocyanin are both produced from metabolites of the phenylpropanoid pathway. Accordingly, the anthocyanin content of Populus overexpressing PtoMYB10 decreased by more than 68%. These results indicate that PtoMYB10 can positively regulate xylary fiber SCW thickening, accompanied by the reprogramming of phenylpropanoid metabolism, which redirects metabolic flux from anthocyanin biosynthesis to monolignol biosynthesis.

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

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

Genome-Wide Identification of R2R3-MYB Transcription Factors: Discovery of a "Dual-Function" Regulator of Gypenoside and Flavonol Biosynthesis in Gynostemma pentaphyllum

The R2R3-MYB gene family participates in several plant physiological processes, especially the regulation of the biosynthesis of secondary metabolites. However, little is known about the functions of R2R3-MYB genes in Gynostemma pentaphyllum (G. pentaphyllum), a traditional Chinese medicinal herb that is an excellent source of gypenosides (a class of triterpenoid saponins) and flavonoids. In this study, a systematic genome-wide analysis of the R2R3-MYB gene family was performed using the recently sequenced G. pentaphyllum genome. In total, 87 R2R3-GpMYB genes were identified and subsequently divided into 32 subgroups based on phylogenetic analysis. The analysis was based on conserved exon-intron structures and motif compositions within the same subgroup. Collinearity analysis demonstrated that segmental duplication events were majorly responsible for the expansion of the R2R3-GpMYB gene family, and Ka/Ks analysis indicated that the majority of the duplicated R2R3-GpMYB genes underwent purifying selection. A combination of transcriptome analysis and quantitative reverse transcriptase-PCR (qRT-PCR) confirmed that Gynostemma pentaphyllum myeloblastosis 81 (GpMYB81) along with genes encoding gypenoside and flavonol biosynthetic enzymes exhibited similar expression patterns in different tissues and responses to methyl jasmonate (MeJA). Moreover, GpMYB81 could bind to the promoters of Gynostemma pentaphyllum farnesyl pyrophosphate synthase 1 (GpFPS1) and Gynostemma pentaphyllum chalcone synthase (GpCHS), the key structural genes of gypenoside and flavonol biosynthesis, respectively, and activate their expression. Altogether, this study highlights a novel transcriptional regulatory mechanism that suggests that GpMYB81 acts as a "dual-function" regulator of gypenoside and flavonol biosynthesis in G. pentaphyllum.

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.

Genome-wide identification and genetic characterization of the CaMYB family and its response to five types of heavy metal stress in hot pepper (Capsicum annuum cv. CM334)

MYB proteins play a crucial role in plant growth and development and stress responses. In this study, 160 members of the MYB gene family from the pepper genome database were used to analyze gene structures, chromosome localization, collinearity, genetic affinity and expression in response to heavy metals. The results identified R2R3-MYB members and further phylogenetically classified them into 35 subgroups based on highly conserved gene structures and motifs. Collinearity analysis showed that segmental duplication events played a crucial role in the functional expansion of the CaMYB gene family by intraspecific collinearity, and at least 12 pairs of CaMYB genes existed between species prior to the differentiation between monocots and dicots. Moreover, the upstream CaMYB genes were mainly localized to the phytohormone elements ABRE and transcription factor elements MYB and MYC. Further analysis revealed that MYB transcription factors were closely associated with a variety of abiotic stress-related proteins (e.g., MAC-complex and SKIP). Under the stress of five metal ions, Cd²⁺, Cu²⁺, Pb²⁺, Zn²⁺, and Fe³⁺, the expression levels of some CaMYB family genes were upregulated. Of these genes, pairing homologous 1 (PH-1), PH-13, and PH-15 in the roots of Capsicum annuum were upregulated to the greatest extent, indicating that these three MYB family members are particularly sensitive to these five metals. This study provides a theoretical reference for the analysis of the molecular regulatory mechanism of MYB family genes in mediating the response to heavy metals in plants. This study reveals the mode of interaction between MYB and a variety of abiotic stress proteins and clarifies the biological functions of CaMYB family members in the regulation of heavy metal stress.

Genome-Wide Identification and Expression Analysis of R2R3-MYB Family Genes Associated with Petal Pigment Synthesis in Liriodendron

The MYB transcription factor family is one of the largest families in plants, and its members have various biological functions. R2R3-MYB genes are involved in the synthesis of pigments that yield petal colors. Liriodendron plants are widely cultivated as ornamental trees owing to their peculiar leaves, tulip-like flowers, and colorful petals. However, the mechanism underlying petal coloring in this species is unknown, and minimal information about MYB genes in Liriodendron is available. Herein, this study aimed to discern gene(s) involved in petal coloration in Liriodendron via genome-wide identification, HPLC, and RT-qPCR assays. In total, 204 LcMYB superfamily genes were identified in the Liriodendron chinense genome, and 85 R2R3-MYB genes were mapped onto 19 chromosomes. Chromosome 4 contained the most (10) R2R3-MYB genes, and chromosomes 14 and 16 contained the fewest (only one). MEME analysis showed that R2R3-MYB proteins in L. chinense were highly conserved and that their exon-intron structures varied. The HPLC results showed that three major carotenoids were uniformly distributed in the petals of L. chinense, while lycopene and β-carotene were concentrated in the orange band region in the petals of Liriodendron tulipifera. Furthermore, the expression profiles via RT-qPCR assays revealed that four R2R3-MYB genes were expressed at the highest levels at the S3P/S4P stage in L. tulipifera. This result combined with the HPLC results showed that these four R2R3-MYB genes might participate in carotenoid synthesis in the petals of L. tulipifera. This work laid a cornerstone for further functional characterization of R2R3-MYB genes in Liriodendron plants.

Genome-wide analysis and expression profiles of PdeMYB transcription factors in colored-leaf poplar (Populus deltoids)

BACKGROUND

MYB transcription factors, comprising one of the largest transcription factor families in plants, play many roles in secondary metabolism, especially in anthocyanin biosynthesis. However, the functions of the PdeMYB transcription factor in colored-leaf poplar remain elusive.

RESULTS

In the present study, genome-wide characterization of the PdeMYB genes in colored-leaf poplar (Populus deltoids) was conducted. A total of 302 PdeMYB transcription factors were identified, including 183 R2R3-MYB, five R1R2R3-MYB, one 4R-MYB, and 113 1R-MYB transcription factor genes. Genomic localization and paralogs of PdeMYB genes mapped 289 genes on 19 chromosomes, with collinearity relationships among genes. The conserved domain, gene structure, and evolutionary relationships of the PdeMYB genes were also established and analyzed. The expression levels of PdeMYB genes were obtained from previous data in green leaf poplar (L2025) and colored leaf poplar (QHP) as well as our own qRT-PCR analysis data in green leaf poplar (L2025) and colored leaf poplar (CHP), which provide valuable clues for further functional characterization of PdeMYB genes.

CONCLUSIONS

The above results provide not only comprehensive insights into the structure and functions of PdeMYB genes but also provide candidate genes for the future improvement of leaf colorization in Populus deltoids.

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

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

Recent Advanced Metabolic and Genetic Engineering of Phenylpropanoid Biosynthetic Pathways

The MYB transcription factors (TFs) are evolving as critical role in the regulation of the phenylpropanoid and tanshinones biosynthetic pathway. MYB TFs relate to a very important gene family, which are involved in the regulation of primary and secondary metabolisms, terpenoids, bioactive compounds, plant defense against various stresses and cell morphology. R2R3 MYB TFs contained a conserved N-terminal domain, but the domain at C-terminal sorts them different regarding their structures and functions. MYB TFs suppressors generally possess particular repressive motifs, such as pdLNLD/ELxiG/S and TLLLFR, which contribute to their suppression role through a diversity of complex regulatory mechanisms. A novel flower specific "NF/YWSV/MEDF/LW" conserved motif has a great potential to understand the mechanisms of flower development. In the current review, we summarize recent advanced progress of MYB TFs on transcription regulation, posttranscriptional, microRNA, conserved motif and propose directions to future prospective research. We further suggest there should be more focus on the investigation for the role of MYB TFs in microalgae, which has great potential for heterologous protein expression system for future perspectives.

Crocus transcription factors CstMYB1 and CstMYB1R2 modulate apocarotenoid metabolism by regulating carotenogenic genes

Two MYB genes have been identified which regulate apocarotenoid metabolism in Crocus sativus. Apocarotenoids like crocin, picrocrocin and safranal are restricted to genus Crocus and are synthesized by oxidative cleavage of zeaxanthin followed by glycosylation reactions. In Crocus sativus, these apocarotenoids are synthesized in stigma part of the flower in developmentally regulated manner. Most of the genes of apocarotenoid pathway are known, however, the mechanism that regulates its tissue and stage specific biosynthesis remains elusive. MYB family was identified as the largest transcription factor family from Crocus transciptome which indicated its possible role in apocarotenoid regulation besides regulating other metabolic pathways. Towards this, we started with identification of 150 MYB genes from Crocus transcriptome databases. The phylogenetic analysis of Crocus MYB genes divided them into 27 clusters. Domain analysis resulted in identification of four groups of MYBs depending upon the number of R repeats present. Expression profiling indicated that 12 MYBs are upregulated in stigma out of which expression of four genes CstMYB1, CstMYB14, CstMYB16 and CstMYB1R2 correlated with crocin accumulation. Transient overexpression of two nuclear localized MYB genes (CstMYB1 and CstMYB1R2) in Crocus confirmed their role in regulating carotenoid metabolism. Yeast-one-hybrid confirmed that CstMYB1 binds to carotenoid cleavage dioxygenase 2 (CCD2) promoter while CstMYB1R2 binds to phytoene synthase (PSY) and CCD2 promoters. Overall, our study established that CstMYB1 and CstMYB1R2 regulate apocarotenoid biosynthesis by directly binding to promoters of pathway genes.

PtrMYB3, a R2R3-MYB Transcription Factor from Poncirus trifoliata, Negatively Regulates Salt Tolerance and Hydrogen Peroxide Scavenging

MYB transcription factors are widely present in plants and play significant roles in abiotic stresses. However, most MYB genes have not been identified in plants and their functions in abiotic stresses are still unknown. In this study, one MYB gene, designated as PtrMYB3, was cloned from trifoliate orange (Poncirus trifoliata (L.) Raf.), and its function in salt tolerance was investigated. PtrMYB3 contains a conserved R2R3-MYB domain, which is the typical property of R2R3-MYB subfamily proteins. Expression profiling under abiotic stresses indicated that PtrMYB3 could be induced by salt, dehydration and cold stresses. PtrMYB3 was found to be localized to the nucleus and possessed transactivation activity. Overexpression of PtrMYB3 by genetic transformation in tobacco impaired its salt tolerance, whereas silencing of PtrMYB3 by VIGS (virus-induced gene silencing) in trifoliate orange conferred significantly enhanced salt tolerance, indicating that PtrMYB3 negatively regulates salt tolerance. Furthermore, a peroxidase gene (PtrPOD) was found to be greatly upregulated in PtrMYB3-silenced trifoliate orange, and a dual LUC (luciferase) assay confirmed that PtrMYB3 could suppress the expression of PtrPOD. The hydrogen peroxide (H₂O₂) accumulation in PtrMYB3 transgenic tobacco plants after salt stress was higher than the wild type (WT), further confirming that overexpression of PtrMYB3 inhibited PtrPOD-mediated H₂O₂ scavenging. Taken together, these results demonstrate that PtrMYB3 negatively regulates salt tolerance, at least in part, due to the excess accumulation of H₂O₂.

Genome-Wide Analysis of the R2R3-MYB Gene Family in Fragaria × ananassa and Its Function Identification During Anthocyanins Biosynthesis in Pink-Flowered Strawberry

The strawberry (Fragaria × ananassa) is an economically important fruit throughout the world. The large R2R3-MYB gene family participates in a variety of plant functions, including anthocyanin biosynthesis. The present study is the first genome-wide analysis of the MYB gene family in the octoploid strawberry and describes the identification and characterization of the family members using the recently sequenced F. × ananassa genome. Specifically, we aimed to identify the key MYBs involved in petal coloration in the pink-flowered strawberry, which increases its ornamental value. A comprehensive, genome-wide analysis of F. × ananassa R2R3-FaMYBs was performed, investigating gene structures, phylogenic relationships, promoter regions, chromosomal locations, and collinearity. A total of 393 R2R3-FaMYB genes were identified in the F. × ananassa genome and divided into 36 subgroups based on phylogenetic analysis. Most genes with similar functions in the same subgroup exhibited similar exon-intron structures and motif compositions. These R2R3-FaMYBs were unevenly distributed over 28 chromosomes. The expansion of the R2R3-FaMYB gene family in the F. × ananassa genome was found to be caused mainly by segmental duplication. The Ka/Ks analysis indicated that duplicated R2R3-FaMYBs mostly experienced purifying selection and showed limited functional divergence after the duplication events. To elucidate which R2R3-FaMYB genes were associated with anthocyanin biosynthesis in the petals of the pink-flowered strawberry, we compared transcriptional changes in different flower developmental stages using RNA-seq. There were 131 differentially expressed R2R3-FaMYB genes identified in the petals, of which three genes, FaMYB28, FaMYB54, and FaMYB576, appeared likely, based on the phylogenetic analysis, to regulate anthocyanin biosynthesis. The qRT-PCR showed that these three genes were more highly expressed in petals than in other tissues (fruit, leaf, petiole and stolon) and their expressions were higher in red compared to pink and white petals. These results facilitate the clarification on the roles of the R2R3-FaMYB genes in petal coloration in the pink-flowered strawberry. This work provides useful information for further functional analysis on the R2R3-FaMYB gene family in F. × ananassa.

Genome-wide analysis of R2R3-MYB transcription factors family in the autopolyploid Saccharum spontaneum: an exploration of dominance expression and stress response

BACKGROUND

Sugarcane (Saccharum) is the most critical sugar crop worldwide. As one of the most enriched transcription factor families in plants, MYB genes display a great potential to contribute to sugarcane improvement by trait modification. We have identified the sugarcane MYB gene family at a whole-genome level through systematic evolution analyses and expression profiling. R2R3-MYB is a large subfamily involved in many plant-specific processes.

RESULTS

A total of 202 R2R3-MYB genes (356 alleles) were identified in the polyploid Saccharum spontaneum genomic sequence and classified into 15 subgroups by phylogenetic analysis. The sugarcane MYB family had more members by a comparative analysis in sorghum and significant advantages among most plants, especially grasses. Collinearity analysis revealed that 70% of the SsR2R3-MYB genes had experienced duplication events, logically suggesting the contributors to the MYB gene family expansion. Functional characterization was performed to identify 56 SsR2R3-MYB genes involved in various plant bioprocesses with expression profiling analysis on 60 RNA-seq databases. We identified 22 MYB genes specifically expressed in the stem, of which RT-qPCR validated MYB43, MYB53, MYB65, MYB78, and MYB99. Allelic expression dominance analysis implied the differential expression of alleles might be responsible for the high expression of MYB in the stem. MYB169, MYB181, MYB192 were identified as candidate C4 photosynthetic regulators by C4 expression pattern and robust circadian oscillations. Furthermore, stress expression analysis showed that MYB36, MYB48, MYB54, MYB61 actively responded to drought treatment; 19 and 10 MYB genes were involved in response to the sugarcane pokkah boeng and mosaic disease, respectively.

CONCLUSIONS

This is the first report on genome-wide analysis of the MYB gene family in sugarcane. SsMYBs probably played an essential role in stem development and the adaptation of various stress conditions. The results will provide detailed insights and rich resources to understand the functional diversity of MYB transcription factors and facilitate the breeding of essential traits in sugarcane.

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

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

VvMYB114 mediated by miR828 negatively regulates trichome development of Arabidopsis

Trichome is a specialized structure differentiated during the morphogenesis of plant leaf epidermal cells. In recent years, with the continuous researches on trichome development of Arabidopsis and other plants, more and more genes related to trichome morphogenesis have been discovered, including R2R3-type MYB genes. In this study, we cloned a R2R3-type MYB family gene from grape, VvMYB114, a target gene of vvi-miR828. qRT-PCR showed that VvMYB114 mRNA accumulated during grape fruit ripening, and VvMYB114 protein had transcriptional activation activity. Heterologous overexpression of VvMYB114 in Arabidopsis reduced the number of trichome on leaves and stems. Mutating the miR828-binding site in VvMYB114 without altering amino-acid sequence had no effect on trichome development in Arabidopsis. The results showed a different role of the regulation of miR828 to VvMYB114 in Arabidopsis from in grape, which indicated the functional divergence of miRNA targeting homoeologous genes in different species played an important roles in evolution and useful trait selection.

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.

Transcriptional Factor-Mediated Regulation of Active Component Biosynthesis in Medicinal Plants

Plants produce thousands of chemically diverse secondary metabolites, many of which have valuable pharmaceutical properties. There is much interest in the synthesis of these pharmaceuticallyvaluable compounds, including the key enzymes and the transcription factors involved. The function and regulatory mechanism of transcription factors in biotic and abiotic stresses have been studied in depth. However, their regulatory roles in the biosynthesis of bioactive compounds, especially in medicinal plants, have only begun. Here, we review what is currently known about how transcription factors contribute to the synthesis of bioactive compounds (alkaloids, terpenoids, flavonoids, and phenolic acids) in medicinal plants. Recent progress has been made in the cloning and characterization of transcription factors in medicinal plants on the genome scale. So far, several large transcription factors have been identified in MYB, WRKY, bHLH, ZIP, AP2/ERF transcription factors. These transcription factors have been predicted to regulate bioactive compound production. These transcription factors positively or negatively regulate the expression of multiple genes encoding key enzymes, and thereby control the metabolic flow through the biosynthetic pathway. Although the research addressing this niche topic is in its infancy, significant progress has been made, and advances in high-throughput sequencing technology are expected to accelerate the discovery of key regulatory transcription factors in medicinal plants. This review is likely to be useful for those interested in the synthesis of pharmaceutically- valuable plant compounds, especially those aiming to breed or engineer plants that produce greater yields of these compounds.

Enzymatic characterization of ancestral/group-IV clade xyloglucan endotransglycosylase/hydrolase enzymes reveals broad substrate specificities

Xyloglucan endotransglycosylase/hydrolase (XTH) enzymes play important roles in cell wall remodelling. Although previous studies have shown a pathway of evolution for XTH genes from bacterial licheninases, through plant endoglucanases (EG16), the order of development within the phylogenetic clades of true XTHs is yet to be elucidated. In addition, recent studies have revealed interesting and potentially useful patterns of transglycosylation beyond the standard xyloglucan-xyloglucan donor/acceptor substrate activities. To study evolutionary relationships and to search for enzymes with useful broad substrate specificities, genes from the 'ancestral' XTH clade of two monocots, Brachypodium distachyon and Triticum aestivum, and two eudicots, Arabidopsis thaliana and Populus tremula, were investigated. Specific activities of the heterologously produced enzymes showed remarkably broad substrate specificities. All the enzymes studied had high activity with the cellulose analogue HEC (hydroxyethyl cellulose) as well as with mixed-link β-glucan as donor substrates, when compared with the standard xyloglucan. Even more surprising was the wide range of acceptor substrates that these enzymes were able to catalyse reactions with, opening a broad range of possible roles for these enzymes, both within plants and in industrial, pharmaceutical and medical fields. Genome screening and expression analyses unexpectedly revealed that genes from this clade were found only in angiosperm genomes and were predominantly or solely expressed in reproductive tissues. We therefore posit that this phylogenetic group is significantly different and should be renamed as the group-IV clade.

Poplar MYB117 promotes anthocyanin synthesis and enhances flavonoid B-ring hydroxylation by up-regulating the flavonoid 3',5'-hydroxylase gene

Flavonoids, such as anthocyanins, proanthocyanidins, and flavonols, are widespread plant secondary metabolites and important for plant adaptation to diverse abiotic and biotic stresses. Flavonoids can be variously hydroxylated and decorated; their biological activity is partly dependent on the degree of hydroxylation of the B-ring. Flavonoid biosynthesis is regulated by MYB transcription factors, which have been identified and characterized in a diversity of plants. Here we characterize a new MYB activator, MYB117, in hybrid poplar (Populus tremula×tremuloides). When overexpressed in transgenic poplar plants, MYB117 enhanced anthocyanin accumulation in all tissues. Transcriptome analysis of MYB117-overexpressing poplars confirmed the up-regulation of flavonoid and anthocyanin biosynthesis genes, as well as two flavonoid 3',5'-hydroxylase (F3'5'H) genes. We also identified up-regulated cytochrome b5 genes, required for full activity of F3'5'H . Phytochemical analysis demonstrated a corresponding increase in B-ring hydroxylation of anthocyanins, proanthocyanidins, and flavonols in these transgenics. Similarly, overexpression of F3'5'H1 directly in hybrid poplar also resulted in increased B-ring hydroxylation, but without affecting overall flavonoid content. However, the overexpression of the cytochrome b5 gene in F3'5'H1-overexpressing plants did not further increase B-ring hydroxylation. Our data indicate that MYB117 regulates the biosynthesis of anthocyanins in poplar, but also enhances B-ring hydroxylation by up-regulating F3'5'H1.

Genome-Wide Identification and Expression Analysis of MYB Transcription Factors and Their Responses to Abiotic Stresses in Woodland Strawberry (Fragaria vesca)

Woodland strawberry (Fragaria vesca) is a diploid strawberry that is widely used as a model of cultivated octoploid strawberry (Fragaria × ananassa). It has also been used as a model for Rosaceae fruits, non-climacteric fruits, and stolons. The MYB superfamily is the largest transcription factor family in plants, and its members play important roles in plant growth and development. However, the complete MYB superfamily in woodland strawberry has not been studied. In this study, a total of 217 MYB genes were identified in woodland strawberry and classified into four groups: one 4R-MYB protein, five 3R-MYB proteins, 113 2R-MYB proteins, and 98 1R-MYB proteins. The phylogenetic relationship of each MYB subgroup was consistent in terms of intron/exon structure and conserved motif composition. The MYB genes in woodland strawberry underwent loss and expansion events during evolution. The transcriptome data revealed that most FveMYB genes are expressed in several organs, whereas 15 FveMYB genes exhibit organ-specific expression, including five genes (FveMYB101, -112, -44, and -8; FveMYB1R81) in roots, two genes (FveMYB62 and -77) in stolon tips, three genes (FveMYB99 and -35; FveMYB1R96) in open flowers, and five genes (FveMYB76 and -100; FveMYB1R4, -5, and -86) in immature fruits. During fruit ripening of woodland strawberry, the expression levels of 84 FveMYB genes were decreased, of which five genes (FveMYB4, -22, -50, and -66; FveMYB1R57) decreased more than 10-fold, whereas those 18 FveMYB genes were increased, especially FveMYB10 and FveMYB74 increased more than 30-fold. In addition, the expression levels of 36, 68, 52, and 62 FveMYB genes were altered by gibberellic acid, abscisic acid, cold, and heat treatments, respectively, and among them, several genes exhibited similar expression patterns for multiple treatments, suggesting possible roles in the crosstalk of multiple signaling pathways. This study provides candidate genes for the study of stolon formation, fruit development and ripening, and abiotic stress responses.

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

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

MYB Transcription Factors and Its Regulation in Secondary Cell Wall Formation and Lignin Biosynthesis during Xylem Development

The secondary wall is the main part of wood and is composed of cellulose, xylan, lignin, and small amounts of structural proteins and enzymes. Lignin molecules can interact directly or indirectly with cellulose, xylan and other polysaccharide molecules in the cell wall, increasing the mechanical strength and hydrophobicity of plant cells and tissues and facilitating the long-distance transportation of water in plants. MYBs (v-myb avian myeloblastosis viral oncogene homolog) belong to one of the largest superfamilies of transcription factors, the members of which regulate secondary cell-wall formation by promoting/inhibiting the biosynthesis of lignin, cellulose, and xylan. Among them, MYB46 and MYB83, which comprise the second layer of the main switch of secondary cell-wall biosynthesis, coordinate upstream and downstream secondary wall synthesis-related transcription factors. In addition, MYB transcription factors other than MYB46/83, as well as noncoding RNAs, hormones, and other factors, interact with one another to regulate the biosynthesis of the secondary wall. Here, we discuss the biosynthesis of secondary wall, classification and functions of MYB transcription factors and their regulation of lignin polymerization and secondary cell-wall formation during wood formation.

Synergies and Entanglement in Secondary Cell Wall Development and Abiotic Stress Response in Trees

A major challenge for sustainable food, fuel, and fiber production is simultaneous genetic improvement of yield, biomass quality, and resilience to episodic environmental stress and climate change. For Populus and other forest trees, quality traits involve alterations in the secondary cell wall (SCW) of wood for traditional uses, as well as for a growing diversity of biofuels and bioproducts. Alterations in wood properties that are desirable for specific end uses can have negative effects on growth and stress tolerance. Understanding of the diverse roles of SCW genes is necessary for the genetic improvement of fast-growing, short-rotation trees that face perennial challenges in their growth and development. Here, we review recent progress into the synergies and antagonisms of SCW development and abiotic stress responses, particularly, the roles of transcription factors, SCW biogenesis genes, and paralog evolution.

Overexpression of PtoCYCD3;3 Promotes Growth and Causes Leaf Wrinkle and Branch Appearance in Populus

D-type cyclin (cyclin D, CYCD), combined with cyclin-dependent kinases (CDKs), participates in the regulation of cell cycle G1/S transition and plays an important role in cell division and proliferation. CYCD could affect the growth and development of herbaceous plants, such as Arabidopsis thaliana, by regulating the cell cycle process. However, its research in wood plants (e.g., poplar) is poor. Phylogenetic analysis showed that in Populus trichocarpa, CYCD3 genes expanded to six members, namely PtCYCD3;1–6. P. tomentosa CYCD3 genes were amplified based on the CDS region of P. trichocarpa CYCD3 genes. PtoCYCD3;3 showed the highest expression in the shoot tip, and the higher expression in young leaves among all members. Therefore, this gene was selected for further study. The overexpression of PtoCYCD3;3 in plants demonstrated obvious morphological changes during the observation period. The leaves became enlarged and wrinkled, the stems thickened and elongated, and multiple branches were formed by the plants. Anatomical study showed that in addition to promoting the differentiation of cambium tissues and the expansion of stem vessel cells, PtoCYCD3;3 facilitated the division of leaf adaxial epidermal cells and palisade tissue cells. Yeast two-hybrid experiment exhibited that 12 PtoCDK proteins could interact with PtoCYCD3;3, of which the strongest interaction strength was PtoCDKE;2, whereas the weakest was PtoCDKG;3. Molecular docking experiments further verified the force strength of PtoCDKE;2 and PtoCDKG;3 with PtoCYCD3;3. In summary, these results indicated that the overexpression of PtoCYCD3;3 significantly promoted the vegetative growth of Populus, and PtoCYCD3;3 may interact with different types of CDK proteins to regulate cell cycle processes.