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

Genome-wide identification of R2R3-MYB transcription factors in Betula platyphylla and functional analysis of BpMYB95 in salt tolerance

The Myeloblastosis (MYB) transcription factor (TF) family is one of the largest transcription factor families in plants and plays an important role in various physiological processes. At present, there are few reports on birch (Betula platyphylla Suk.) of R2R3-MYB-TFs, and most BpMYBs still need to be characterized. In this study, 111 R2R3-MYB-TFs with conserved R2 and R3 MYB domains were identified. Phylogenetic tree analysis showed that the MYB family members of Arabidopsis thaliana and birch were divided into 23 and 21 subgroups, respectively. The latter exhibited an uneven distribution across 14 chromosomes. There were five tandem duplication events and 17 segmental duplication events between BpMYBs, and repeat events play an important role in the expansion of the family. In addition, the promoter region of MYBs was rich in various cis-acting elements, and MYB-TFs were involved in plant growth and development, light responses, biotic stress, and abiotic stress. RNA-sequencing (RNA-seq) and quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) results revealed that most R2R3-MYB-TFs in birch responded to salt stress. In particular, the expression of BpMYBs in the S20 subfamily was significantly induced by salt, drought, abscisic acid, and methyl jasmonate stresses. Based on the weighted co-expression network analysis of physiological and RNA-seq data of birch under salt stress, a key MYB-TF BpMYB95 (BPChr12G24087), was identified in response to salt stress, and its expression level was induced by salt stress. BpMYB95 is a nuclear localization protein with transcriptional activation activity in yeast and overexpression of this gene significantly enhanced salt tolerance in Saccharomyces cerevisiae. The qRT-PCR and histochemical staining results showed that BpMYB95 exhibited the highest expression in the roots, young leaves, and petioles of birch plants. Overexpression of BpMYB95 significantly improved salt-induced browning and wilting symptoms in birch leaves and alleviated the degree of PSII photoinhibition caused by salt stress in birch seedlings. In conclusion, most R2R3-MYB-TFs found in birch were involved in the salt stress response mechanisms. Among these, BpMYB95 was a key regulatory factor that significantly enhanced salt tolerance in birch. The findings of this study provide valuable genetic resources for the development of salt-tolerant birch varieties.

Selection and Validation of Reliable Reference Genes for Liquidambar formosana Leaves with Different Leaf Colors

Liquidambar formosana Hance is renowned for its rich leaf color and possesses notable advantages, such as robust adaptability, strong resistance to diseases and pests, and rapid growth, making it a preferred choice for urban greening and carbon sequestration forest initiatives. The completion of whole-genome sequencing of L. formosana has spurred an increased interest in exploring the molecular mechanisms underlying seasonal changes in leaf color, marking a significant focus in L. formosana breeding research. However, there is currently a lack of stable reference genes suitable for analyzing the expression patterns of functional genes in L. formosana exhibiting varying leaf colors. This study selected five L. formosana varieties with significant differences in leaf colors. Through the RT-qPCR analysis, and evaluation using BestKeeper, geNorm, NormFinder, Delta Ct, and RefFinder, the expression stability of 14 candidate reference genes was examined. Consequently, two reference genes (LifEF1-α and LifACT) with stable expression, suitable for RT-qPCR of L. formosana with diverse leaf colors, were identified. The stability of these selected reference genes was further validated by examining the LifbHLH137 gene, which promoted the biosynthesis of anthocyanins. This advancement facilitated molecular biology and genetic breeding investigations of L. formosana, providing essential data for the precise quantification of functional genes associated with leaf color variation.

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

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

Transcription Factors and Their Regulatory Roles in the Male Gametophyte Development of Flowering Plants

Male gametophyte development in plants relies on the functions of numerous genes, whose expression is regulated by transcription factors (TFs), non-coding RNAs, hormones, and diverse environmental stresses. Several excellent reviews are available that address the genes and enzymes associated with male gametophyte development, especially pollen wall formation. Growing evidence from genetic studies, transcriptome analysis, and gene-by-gene studies suggests that TFs coordinate with epigenetic machinery to regulate the expression of these genes and enzymes for the sequential male gametophyte development. However, very little summarization has been performed to comprehensively review their intricate regulatory roles and discuss their downstream targets and upstream regulators in this unique process. In the present review, we highlight the research progress on the regulatory roles of TF families in the male gametophyte development of flowering plants. The transcriptional regulation, epigenetic control, and other regulators of TFs involved in male gametophyte development are also addressed.

Integrated metabolome, full-length sequencing, and transcriptome analyses unveil the molecular mechanisms of color formation of the canary yellow and red bracts of Bougainvillea × buttiana 'Chitra'

Bougainvillea is a typical tropical flower of great ornamental value due to its colorful bracts. The molecular mechanism behind color formation is not well-understood. Therefore, this research conducted metabolome analysis, transcriptome analysis, and multi-flux full-length sequencing in two color bracts of Bougainvillea × buttiana 'Chitra' to investigate the significantly different metabolites (SDMs) and differentially expressed genes (DEGs). Overall, 261 SDMs, including 62 flavonoids and 26 alkaloids, were detected, and flavonols and betalains were significantly differentially accumulated among the two bracts. Furthermore, the complete-length transcriptome of Bougainvillea × buttiana was also developed, which contained 512 493 non-redundant isoforms. Among them, 341 210 (66.58%) displayed multiple annotations in the KOG, GO, NR, KEGG, Pfam, Swissprot, and NT databases. RNA-seq findings revealed that 3610 DEGs were identified between two bracts. Co-expression analysis demonstrated that the DEGs and SDMs involved in flavonol metabolism (such as CHS, CHI, F3H, FLS, CYP75B1, kaempferol, and quercetin) and betacyanin metabolism (DODA, betanidin, and betacyanins) were the main contributors for the canary yellow and red bract formation, respectively. Further investigation revealed that several putative transcription factors (TFs) might interact with the promoters of the genes mentioned above. The expression profiles of the putative TFs displayed that they may positively and negatively regulate the structural genes' expression profiles. The data revealed a potential regulatory network between important genes, putative TFs, and metabolites in the flavonol and betacyanin biosynthesis of Bougainvillea × buttiana 'Chitra' bracts. These findings will serve as a rich genetic resource for future studies that could create new color bracts.

The Over-Expression of Two R2R3-MYB Genes, PdMYB2R089 and PdMYB2R151, Increases the Drought-Resistant Capacity of Transgenic Arabidopsis

The R2R3-MYB genes in plants play an essential role in the drought-responsive signaling pathway. Plenty of R2R3-MYB S21 and S22 subgroup genes in Arabidopsis have been implicated in dehydration conditions, yet few have been covered in terms of the role of the S21 and S22 subgroup genes in poplar under drought. PdMYB2R089 and PdMYB2R151 genes, respectively belonging to the S21 and S22 subgroups of NL895 (Populus deltoides × P. euramericana cv. 'Nanlin895'), were selected based on the previous expression analysis of poplar R2R3-MYB genes that are responsive to dehydration. The regulatory functions of two target genes in plant responses to drought stress were studied and speculated through the genetic transformation of Arabidopsis thaliana. PdMYB2R089 and PdMYB2R151 could promote the closure of stomata in leaves, lessen the production of malondialdehyde (MDA), enhance the activity of the peroxidase (POD) enzyme, and shorten the life cycle of transgenic plants, in part owing to their similar conserved domains. Moreover, PdMYB2R089 could strengthen root length and lateral root growth. These results suggest that PdMYB2R089 and PdMYB2R151 genes might have the potential to improve drought adaptability in plants. In addition, PdMYB2R151 could significantly improve the seed germination rate of transgenic Arabidopsis, but PdMYB2R089 could not. This finding provides a clue for the subsequent functional dissection of S21 and S22 subgroup genes in poplar that is responsive to drought.

Function and molecular mechanism of a poplar placenta limited MIXTA gene in regulating differentiation of plant epidermal cells

The placenta in fruits of most plants either desiccate and shrink as the fruits mature or develop further to form the fleshy tissues. In poplars, placental epidermal cells protrude collectively to produce catkin fibers. In this study, three carpel limited MIXTA genes, PdeMIXTA02, PdeMIXTA03, PdeMIXTA04, were find to specifically expressed in carpel immediately after pollination. Heterologous expression of the three genes in Arabidopsis demonstrated that PdeMIXTA04 significantly promoted trichomes density and could restore trichomes in the trichomeless mutant. By contrast, such functions were not observed with PdeMIXTA02, PdeMIXTA03. In situ hybridization revealed that PdeMIXTA04 was explicitly expressed in poplar placental epidermal cells. We also confirmed trichome-specific expression of the PdeMIXTA04 promoter. Multiple experimental proofs have confirmed the interaction between PdeMIXTA04, PdeMYC and PdeWD40, indicating PdeMIXTA04 functioned through the MYB-bHLH-WD40 ternary complex. Our work provided distinctive understanding of the molecular mechanism triggering differentiation of poplar catkins.

Systematic investigation of the R2R3-MYB gene family in Aquilaria sinensis reveals a transcriptional repressor AsMYB054 involved in 2-(2-phenylethyl)chromone biosynthesis

Trees in the genus Aquilaria produce agarwood, a valuable resin used in medicine, perfumes, and incense. 2-(2-Phenethyl)chromones (PECs) are characteristic components of agarwood; however, molecular mechanisms underlying PEC biosynthesis and regulation remain largely unknown. The R2R3-MYB transcription factors play important regulatory roles in the biosynthesis of various secondary metabolites. In this study, 101 R2R3-MYB genes in Aquilaria sinensis were systematically identified and analyzed at the genome-wide level. Transcriptomic analysis revealed that 19 R2R3-MYB genes were significantly regulated by an agarwood inducer, and showed significant correlations with PEC accumulation. Expression and evolutionary analyses revealed that AsMYB054, a subgroup 4 R2R3-MYB, was negatively correlated with PEC accumulation. AsMYB054 was located in the nucleus and functioned as a transcriptional repressor. Moreover, AsMYB054 could bind to the promoters of the PEC biosynthesis related genes AsPKS02 and AsPKS09, and inhibit their transcriptional activity. These findings suggested that AsMYB054 functions as a negative regulator of PEC biosynthesis via the inhibition of AsPKS02 and AsPKS09 in A. sinensis. Our results provide a comprehensive understanding of the R2R3-MYB subfamily in A. sinensis and lay a foundation for further functional analyses of R2R3-MYB genes in PEC biosynthesis.

Genome-Wide Analysis of the Expansin Gene Family in Populus and Characterization of Expression Changes in Response to Phytohormone (Abscisic Acid) and Abiotic (Low-Temperature) Stresses

Expansins are a group of cell wall enzyme proteins that help to loosen cell walls by breaking hydrogen bonds between cellulose microfibrils and hemicellulose. Expansins are essential plant proteins that are involved in several key processes, including seed germination, the growth of pollen tubes and root hairs, fruit ripening and abscission processes. Currently, there is a lack of knowledge concerning the role of expansins in woody plants. In this study, we analyzed expansin genes using Populus genome as the study target. Thirty-six members of the expansin gene family were identified in Populus that were divided into four subfamilies (EXPA, EXPB, EXLA and EXLB). We analyzed the molecular structure, chromosome localization, evolutionary relationships and tissue specificity of these genes and investigated expression changes in responses to phytohormone and abiotic stresses of the expansin genes of Populus tremula L. (PtEXs). Molecular structure analysis revealed that each PtEX protein had several conserved motifs and all of the PtEXs genes had multiple exons. Chromosome structure analysis showed that the expansin gene family is distributed on 14 chromosomes. The PtEXs gene family expansion patterns showed segmental duplication. Transcriptome data of Populus revealed that 36 PtEXs genes were differently expressed in different tissues. Cis-element analysis showed that the PtEXs were closely associated with plant development and responses to phytohormone and abiotic stress. Quantitative real-time PCR showed that abscisic acid (ABA) and low-temperature treatment affected the expression of some PtEXs genes, suggesting that these genes are involved in responses to phytohormone and abiotic stress. This study provides a further understanding of the expansin gene family in Populus and forms a basis for future functional research studies.

Comprehensive Genome-Wide Analyses of Poplar R2R3-MYB Transcription Factors and Tissue-Specific Expression Patterns under Drought Stress

R2R3-type MYB transcription factors are implicated in drought stress, which is a primary factor limiting the growth and development of woody plants. The identification of R2R3-MYB genes in the Populus trichocarpa genome has been previously reported. Nevertheless, the diversity and complexity of the conserved domain of the MYB gene caused inconsistencies in these identification results. There is still a lack of drought-responsive expression patterns and functional studies of R2R3-MYB transcription factors in Populus species. In this study, we identified a total of 210 R2R3-MYB genes in the P. trichocarpa genome, of which 207 genes were unevenly distributed across all 19 chromosomes. These poplar R2R3-MYB genes were phylogenetically divided into 23 subgroups. Collinear analysis demonstrated that the poplar R2R3-MYB genes underwent rapid expansion and that whole-genome duplication events were a dominant factor in the process of rapid gene expansion. Subcellular localization assays indicated that poplar R2R3-MYB TFs mainly played a transcriptional regulatory role in the nucleus. Ten R2R3-MYB genes were cloned from P. deltoides × P. euramericana cv. Nanlin895, and their expression patterns were tissue-specific. A majority of the genes showed similar drought-responsive expression patterns in two out of three tissues. This study provides a valid cue for further functional characterization of drought-responsive R2R3-MYB genes in poplar and provides support for the development of new poplar genotypes with elevated drought tolerance.

Genome-wide characterization and functional identification of MYB genes in Malus sieversii infected by Valsa mali

Among the most important transcription factors in plants, the v-myb avian myeloblastosis viral oncogene homolog (MYB) regulates the expression network of response genes under stresses such as fungal infection. In China, the canker disease Valsa mali threatens the survival of Malus sieversii, an ancestor of cultivated apples. Using the M. sieversii genome, we identified 457 MsMYB and 128 R2R3-MsMYB genes that were randomly distributed across 17 chromosomes. Based on protein sequence and structure, the R2R3-MsMYB genes were phylogenetically divided into 29 categories, and 26 conserved motifs were identified. We further predicted cis-elements in the 2000-kb promoter region of R2R3-MsMYBs based on the genome. Transcriptome analysis of M. sieversii under V. mali infection showed that 27 R2R3-MsMYBs were significantly differentially expressed, indicating their key role in the response to V. mali infection. Using transient transformation, MsMYB14, MsMYB24, MsMYB39, MsMYB78, and MsMYB108, which were strongly induced by V. mali infection, were functionally identified. Among the five MsMYBs, MsMYB14 and MsMYB78 were both important in enhancing resistance to diseases, whereas MsMYB24 inhibited resistance. Based on the results of this study, we gained a better understanding of the MsMYB transcription factor family and laid the foundation for a future research program on disease prevention strategies in M. sieversii.

Genome-Wide Characterization and Evolutionary Expansion of Poplar NAC Transcription Factors and Their Tissue-Specific Expression Profiles under Drought

The NAC (NAM, ATAF1/2 and CUC2) is a large gene family of plant-specific transcription factors that play a pivotal role in various physiological processes and abiotic stresses. Due to the lack of genome-wide characterization, intraspecific and interspecific synteny, and drought-responsive expression pattern of NAC genes in poplar, the functional characterization of drought-related NAC genes have been scarcely reported in Populus species. Here, we identified a total of 170 NAC domain-containing genes in the P. trichocarpa genome, 169 of which were unevenly distributed on its nineteen chromosomes. These NAC genes were phylogenetically divided into twenty subgroups, some of which exhibited a similar pattern of exon-intron architecture. The synteny and Ka/Ks analysis indicated that the expansion of NAC genes in poplar was mainly due to gene duplication events occurring before and after the divergence of Populus and Salix. Ten PdNAC (P. deltoids × P. euramericana cv.'Nanlin895') genes were randomly selected and cloned. Their drought-responsive expression profiles showed a tissue-specific pattern. The transcription factor PdNAC013 was verified to be localized in the nucleus. Our research results provide genomic information for the expansion of NAC genes in the poplar genome, and for further characterizing putative poplar NAC genes associated with water-deficit.

Genome-Wide Analysis of Specific PfR2R3-MYB Genes Related to Paulownia Witches' Broom

Paulownia witches' broom (PaWB), caused by phytoplasmas, is the most devastating infectious disease of Paulownia. R2R3-MYB transcription factors (TF) have been reported to be involved in the plant's response to infections caused by these pathogens, but a comprehensive study of the R2R3-MYB genes in Paulownia has not been reported. In this study, we identified 138 R2R3-MYB genes distributed on 20 chromosomes of Paulownia fortunei. These genes were classified into 27 subfamilies based on their gene structures and phylogenetic relationships, which indicated that they have various evolutionary relationships and have undergone rich segmental replication events. We determined the expression patterns of the 138 R2R3-MYB genes of P. fortunei by analyzing the RNA sequencing data and found that PfR2R3-MYB15 was significantly up-regulated in P. fortunei in response to phytoplasma infections. PfR2R3-MYB15 was cloned and overexpressed in Populus trichocarpa. The results show that its overexpression induced branching symptoms. Subsequently, the subcellular localization results showed that PfR2R3-MYB15 was located in the nucleus. Yeast two-hybrid and bimolecular fluorescence complementation experiments showed that PfR2R3-MYB15 interacted with PfTAB2. The analysis of the PfR2R3-MYB15 gene showed that it not only played an important role in plant branching, but also might participate in the biosynthesis of photosystem elements. Our results will provide a foundation for future studies of the R2R3-MYB TF family in Paulownia and other plants.

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.

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

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

A Selection of Reliable Reference Genes for Gene Expression Analysis in the Female and Male Flowers of Salix suchowensis

Salix is a dioecious plant. Research on the molecular regulation mechanism of male and female inflorescence differentiation and development is necessary to analyze sex differentiation in the willow and the underlying mechanisms of unisexual flower development. However, at present, there are no reference genes suitable for stable expression in the process of willow inflorescence development. In this study, Salix suchowensis was used as the research material, nine candidate reference genes (α-TUB1, α-TUB2, ACT, H2A, DnaJ, CDC2, GAPDH, TIP41, β-TUB) were selected, and qRT-PCR technology was used to detect the expression of each candidate reference gene in female and male flowers at different developmental stages and using five algorithms (geNorm, Normfinder, Delta Ct, BestKeeper, and RefFinder) to comprehensively evaluate the stability of candidate reference genes. The results showed that ACT and DnaJ were stably expressed in all samples and could be used as reference genes. In addition, the reliability of the screening results was further verified via an expression pattern analysis of the CFS gene that encodes flower specific transcription factor in different samples. The stable reference genes selected in this study provide the basis for future research on the expression analysis of functional genes related to the development of male and female flowers of S. suchowensis.