OsBZR1 turnover mediated by OsSK22-regulated U-box E3 ligase OsPUB24 in rice BR response

The maize GSK3-like kinase ZmSK1 negatively regulates drought tolerance by phosphorylating the transcription factor ZmCPP2

Glycogen synthase kinase 3 (GSK3)-like kinases play important roles in stress responses in plants. However, the mechanism of GSK3-like kinases in drought-induced antioxidant defense is not clear. In this study, we discovered that the GSK3-like kinase SHAGGY-like kinase 1 (ZmSK1) negatively regulates drought tolerance by inhibiting antioxidant defense in maize (Zea mays). Then, we determined that cysteine-rich polycomb-like protein 2 (ZmCPP2) interacts with ZmSK1 and enhances maize drought tolerance by inducing antioxidant defense. ZmCPP2 is phosphorylated at Ser-250 by ZmSK1, which is dependent on ZmSK1 kinase activity and attenuates maize drought tolerance. Furthermore, ZmCPP2 directly binds to the promoter of the superoxide dismutase (SOD) gene ZmSOD4, encoding an antioxidant defense enzyme, and activates its expression. ZmSK1 phosphorylating ZmCPP2 at Ser-250 represses the binding of ZmCPP2 to the ZmSOD4 promoter. Taken together, our results indicate that the phosphorylation of ZmCPP2 by ZmSK1 results in decreased SOD activity and thus reduces drought tolerance in maize. These findings reveal a mechanism of GSK3-like kinases regulating antioxidant defense in the drought stress response.

The BES1/BZR1 transcriptional factor SlBES2 regulates photosynthetic apparatus in tomato fruit

BACKGROUND

Fruit photosynthetic apparatus development comprises a series of biological processes which is essential in determining fruit development and quality formation. However, the understanding of the regulation of fruit photosynthetic apparatus development remains poor.

RESULTS

In this study, we identified a transcriptional factor SlBES2, the closest homolog of BES1 and BZR1 in tomato BES1 family, is highly expressed in fruit at mature green (MG) stage and exhibited transcriptional inhibition activity. Down-regulation of SlBES2 resulted in fruits showing paler fruit than wild type at MG stage, in contrast, SlBES2-overexpressing tomato lines bore deeper green fruits. Notably, chlorophyll content and number of thylakoids per chloroplast in fruit was substantially increased in SlBES2-overexpressing lines, while markedly decreased in SlBES2-suppressing lines. Comparative transcriptome analysis revealed that multiple genes of the photosystem, chloroplast development and chlorophyll metabolism pathways were regulated by SlBES2. Further verification revealed that SlBES2 can significantly repress the transcriptional activity of SlNYC1 and Green-Flesh, and physically interact with protein SlHY5.

CONCLUSIONS

Collectively, this study demonstrated that SlBES2 plays an important role in regulating fruit photosynthetic apparatus development through either transcriptional repression of genes involved in chlorophyll breakdown, or posttranscriptional regulation of proteins associated with plant photomorphogenesis and chloroplast development. Our findings add a new actor to the complex mechanisms underlying photosynthetic apparatus during fruit development.

Identification of a novel locus qGW12/OsPUB23 regulating grain shape and weight in rice (Oryza sativa L.)

KEY MESSAGE

Key message A major quantitative trait locus (qGW12) for grain shape and weight has been isolated in rice, corresponding to LOC_Os12g17900/OsPUB23, and its encoded protein interacts with OsMADS1. Grain shape in rice is an important trait that influences both yield and quality. The primary determinants of grain shape are quantitative trait loci (QTLs) inherited from natural variation in crops. In recent years, much attention has been paid to the molecular role of QTLs in regulating grain shape and weight. In this study, we report the cloning and characterization of qGW12, a major QTL regulating grain shape and weight in rice, using a series of chromosome fragment substitution lines (CSSLs) derived from Oryza sativa indica cultivar 9311 (acceptor) and Oryza rufipogon Griff (donor). One CSSL line, Q187, harboring the introgression of qGW12, exhibited a significant decrease in grain-shape-related traits (including grain length and width) and thousand-grain weight compared to the cultivar 9311. Subsequent backcrossing of Q187 with 9311 resulted in the generation of secondary segregating populations, which were used to fine-map qGW12 to a 24-kb region between markers Seq-44 and Seq-48. Our data indicated that qGW12 encodes a previously unreported U-box type E3 ubiquitin ligase, designated OsPUB23, which exhibited E3 ubiquitin ligase activity. Overexpression of OsPUB23 in rice resulted in higher plant yield than the wild type due to an increase in grain size and weight. Conversely, loss of OsPUB23 function resulted in the opposite tendency. Yeast two-hybrid screening and split luciferase complementation assays revealed that OsPUB23 interacts with OsMADS1. The functional characterization of OsPUB23 provides new genetic resources for improving of grain yield and quality in crops.

Post-translational Regulation of BRI1-EMS Suppressor 1 and Brassinazole-Resistant 1

Brassinosteroid-insensitive 1 (BRI1)-EMS suppressor 1 (BES1) and Brassinazole-resistant 1 (BZR1) are two highly similar master transcription factors of the brassinosteroid (BR) signaling pathway that regulates a variety of plant growth and development processes as well as stress responses. Previous genetic and biochemical analyses have established a complex regulatory network to control the two transcription factors. This network includes coordination with other transcription factors and interactors, multiple post-translational modifications (PTMs) and differential subcellular localizations. In this review, we systematically detail the functions and regulatory mechanisms of various PTMs: phosphorylation/dephosphorylation, ubiquitination/deubiquitination, SUMOylation/deSUMOylation and oxidation/reduction, in regulating the subcellular localization, protein stability and the transcriptional activity of BES1/BZR1. We also discuss the current knowledge about the BES1/BZR1 interactors mediating the dynamic nucleocytoplasmic shuttling of BES1 and BZR1.

Genome-wide identification and expression analysis of the U-box gene family related to biotic and abiotic stresses in Coffea canephora L

Plant U-box genes play an important role in the regulation of plant hormone signal transduction, stress tolerance, and pathogen resistance; however, their functions in coffee (Coffea canephora L.) remain largely unexplored. In this study, we identified 47 CcPUB genes in the C. canephora L. genome, clustering them into nine groups via phylogenetic tree. The CcPUB genes were unevenly distributed across the 11 chromosomes of C. canephora L., with the majority (11) on chromosome 2 and none on chromosome 8. The cis-acting elements analysis showed that CcPUB genes were involved in abiotic and biotic stresses, phytohormone responsive, and plant growth and development. RNA-seq data revealed diverse expression patterns of CcPUB genes across leaves, stems, and fruits tissues. qRT-PCR analyses under dehydration, low temperature, SA, and Colletotrichum stresses showed significant up-regulation of CcPUB2, CcPUB24, CcPUB34, and CcPUB40 in leaves. Furthermore, subcellular localization showed CcPUB2 and CcPUB34 were located in the plasma membrane and nucleus, and CcPUB24 and CcPUB40 were located in the nucleus. This study provides valuable insights into the roles of PUB genes in stress responses and phytohormone signaling in C. canephora L., and provided basis for functional characterization of PUB genes in C. canephora L.

Brassinosteroid biosynthesis and signaling: Conserved and diversified functions of core genes across multiple plant species

Brassinosteroids (BRs) are important regulators that control myriad aspects of plant growth and development, including biotic and abiotic stress responses, such that modulating BR homeostasis and signaling presents abundant opportunities for plant breeding and crop improvement. Enzymes and other proteins involved in the biosynthesis and signaling of BRs are well understood from molecular genetics and phenotypic analysis in Arabidopsis thaliana; however, knowledge of the molecular functions of these genes in other plant species, especially cereal crop plants, is minimal. In this manuscript, we comprehensively review functional studies of BR genes in Arabidopsis, maize, rice, Setaria, Brachypodium, and soybean to identify conserved and diversified functions across plant species and to highlight cases for which additional research is in order. We performed phylogenetic analysis of gene families involved in the biosynthesis and signaling of BRs and re-analyzed publicly available transcriptomic data. Gene trees coupled with expression data provide a valuable guide to supplement future research on BRs in these important crop species, enabling researchers to identify gene-editing targets for BR-related functional studies.

A distinct protein posttranslational modifications-linked OsATL32-OsPPKL2-OsGSK2 loop modulates rice immunity against blast disease

Protein posttranslational modifications play crucial roles in plant immunity through modulating a complicated signaling network mediated by different hormones. We previously demonstrated that OsATL32, an ATL-type E3 ligase, negatively contributes to rice immunity against Magnaporthe oryzae. Here, we show that OsATL32 forms a loop with OsPPKL2 and OsGSK2 through distinct protein posttranslational modifications to modulate rice immunity. OsATL32 ubiquitinates OsPPKL2, a protein phosphatase with Kelch-like repeat domains that exerts positive roles in regulating rice immunity against M. oryzae and chitin-triggered immune responses, for degradation. The glycogen synthase kinase 2 (OsGSK2), which acts as a negative regulator of rice immunity against M. oryzae and chitin-triggered immune responses, phosphorylates OsATL32 to elevate its protein stability and E3 ligase activity on OsPPKL2. Moreover, OsPPKL2 directly dephosphorylates OsGSK2, affecting its kinase activity on substrates including OsATL32 for phosphorylation. Like OsGSK2 as a BR signaling repressor, OsATL32 negatively regulates BR signaling; conversely, OsPPKL2 plays a positive role in BR signaling. These findings provide a molecular mechanism in which OsATL32 serves as a node connecting BR signaling and immunity by associating with OsPPKL2 and OsGSK2, assembling into a distinct protein posttranslational modifications-linked loop that functions in rice BR signaling and immunity.

The TaSnRK1-TabHLH489 module integrates brassinosteroid and sugar signalling to regulate the grain length in bread wheat

Regulation of grain size is a crucial strategy for improving the crop yield and is also a fundamental aspect of developmental biology. However, the underlying molecular mechanisms governing grain development in wheat remain largely unknown. In this study, we identified a wheat atypical basic helix-loop-helix (bHLH) transcription factor, TabHLH489, which is tightly associated with grain length through genome-wide association study and map-based cloning. Knockout of TabHLH489 and its homologous genes resulted in increased grain length and weight, whereas the overexpression led to decreased grain length and weight. TaSnRK1α1, the α-catalytic subunit of plant energy sensor SnRK1, interacted with and phosphorylated TabHLH489 to induce its degradation, thereby promoting wheat grain development. Sugar treatment induced TaSnRK1α1 protein accumulation while reducing TabHLH489 protein levels. Moreover, brassinosteroid (BR) promotes grain development by decreasing TabHLH489 expression through the transcription factor BRASSINAZOLE RESISTANT1 (BZR1). Importantly, natural variations in the promoter region of TabHLH489 affect the TaBZR1 binding ability, thereby influencing TabHLH489 expression. Taken together, our findings reveal that the TaSnRK1α1-TabHLH489 regulatory module integrates BR and sugar signalling to regulate grain length, presenting potential targets for enhancing grain size in wheat.

Interaction of the Transcription Factors BES1/BZR1 in Plant Growth and Stress Response

Bri1-EMS Suppressor 1 (BES1) and Brassinazole Resistant 1 (BZR1) are two key transcription factors in the brassinosteroid (BR) signaling pathway, serving as crucial integrators that connect various signaling pathways in plants. Extensive genetic and biochemical studies have revealed that BES1 and BZR1, along with other protein factors, form a complex interaction network that governs plant growth, development, and stress tolerance. Among the interactome of BES1 and BZR1, several proteins involved in posttranslational modifications play a key role in modifying the stability, abundance, and transcriptional activity of BES1 and BZR1. This review specifically focuses on the functions and regulatory mechanisms of BES1 and BZR1 protein interactors that are not involved in the posttranslational modifications but are crucial in specific growth and development stages and stress responses. By highlighting the significance of the BZR1 and BES1 interactome, this review sheds light on how it optimizes plant growth, development, and stress responses.

The MYC2-PUB22-JAZ4 module plays a crucial role in jasmonate signaling in tomato

Jasmonates (JAs), a class of lipid-derived stress hormones, play a crucial role across an array of plant physiological processes and stress responses. Although JA signaling is thought to rely predominantly on the degradation of specific JAZ proteins by SCFCOI1, it remains unclear whether other pathways are involved in the regulation of JAZ protein stability. Here, we report that PUB22, a plant U-box type E3 ubiquitin ligase, plays a critical role in the regulation of plant resistance against Helicoverpa armigera and other JA responses in tomato. Whereas COI1 physically interacts with JAZ1/2/5/7, PUB22 physically interacts with JAZ1/3/4/6. PUB22 ubiquitinates JAZ4 to promote its degradation via the 26S proteasome pathway. Importantly, we observed that pub22 mutants showreduced resistance to H. armigera, whereas jaz4 single mutants and jaz1 jaz3 jaz4 jaz6 quadruple mutants have enhanced resistance. The hypersensitivity of pub22 mutants to herbivores could be partially rescued by JAZ4 mutation. Moreover, we found that expression of PUB22 can be transcriptionally activated by MYC2, thus forming a positive feedback circuit in JA signaling. We noticed that the PUB22-JAZ4 module also regulates other JA responses, including defense against B. cinerea, inhibition of root elongation, and anthocyanin accumulation. Taken together, these results indicate that PUB22 plays a crucial role in plant growth and defense responses, together with COI1-regulated JA signaling, by targeting specific JAZs.

Alternative polyadenylation profiles of susceptible and resistant rice (Oryza sativa L.) in response to bacterial leaf blight using RNA-seq

BACKGROUND

Alternative polyadenylation (APA) is an important pattern of post-transcriptional regulation of genes widely existing in eukaryotes, involving plant physiological and pathological processes. However, there is a dearth of studies investigating the role of APA profile in rice leaf blight.

RESULTS

In this study, we compared the APA profile of leaf blight-susceptible varieties (CT 9737-613P-M) and resistant varieties (NSIC RC154) following bacterial blight infection. Through gene enrichment analysis, we found that the genes of two varieties typically exhibited distal poly(A) (PA) sites that play different roles in two kinds of rice, indicating differential APA regulatory mechanisms. In this process, many disease-resistance genes displayed multiple transcripts via APA. Moreover, we also found five polyadenylation factors of similar expression patterns of rice, highlighting the critical roles of these five factors in rice response to leaf blight about PA locus diversity.

CONCLUSION

Notably, the present study provides the first dynamic changes of APA in rice in early response to biotic stresses and proposes a possible functional conjecture of APA in plant immune response, which lays the theoretical foundation for in-depth determination of the role of APA events in plant stress response and other life processes.

Optimizing rice grain size by attenuating phosphorylation-triggered functional impairment of a chromatin modifier ternary complex

Histone acetylation affects numerous cellular processes, such as gene transcription, in both plants and animals. However, the posttranslational modification-participated regulatory networks for crop-yield-related traits are largely unexplored. Here, we characterize a regulatory axis for controlling rice grain size and yield, centered on a potent histone acetyltransferase (chromatin modifier) known as HHC4. HHC4 interacts with and forms a ternary complex with adaptor protein ADA2 and transcription factor bZIP23, wherein bZIP23 recruits HHC4 to specific promoters, and ADA2 and HHC4 additively enhance bZIP23 transactivation on target genes. Meanwhile, HHC4 interacts with and is phosphorylated by GSK3-like kinase TGW3. The resultant phosphorylation triggers several functional impairments of the HHC4 ternary complex. In addition, we identify two major phosphorylation sites of HHC4 by TGW3-sites which play an important role in controlling rice grain size. Overall, our findings thus have critical implications for understanding epigenetic basis of grain size control and manipulating the knowledge for higher crop productivity.

Genome-wide analysis of the U-box E3 ligases gene family in potato (Solanum tuberosum L.) and overexpress StPUB25 enhance drought tolerance in transgenic Arabidopsis

BACKGROUND

Plant U-box (PUB) E3 ubiquitin ligases have vital effects on various biological processes. Therefore, a comprehensive and systematic identification of the members of the U-box gene family in potato will help to understand the evolution and function of U-box E3 ubiquitin ligases in plants.

RESULTS

This work identified altogether 74 PUBs in the potato (StPUBs) and examined their gene structures, chromosomal distributions, and conserved motifs. There were seventy-four StPUB genes on ten chromosomes with diverse densities. As revealed by phylogenetic analysis on PUBs within potato, Arabidopsis, tomato (Solanum lycopersicum), cabbage (Brassica oleracea), rice (Oryza sativa), and corn (Zea mays), were clustered into eight subclasses (C1-C8). According to synteny analysis, there were 40 orthologous StPUB genes to Arabidopsis, 58 to tomato, 28 to cabbage, 7 to rice, and 8 to corn. In addition, RNA-seq data downloaded from PGSC were utilized to reveal StPUBs' abiotic stress responses and tissue-specific expression in the doubled-monoploid potato (DM). Inaddition, we performed RNA-seq on the 'Atlantic' (drought-sensitive cultivar, DS) and the 'Qingshu NO.9' (drought-tolerant cultivar, DT) in early flowering, full-blooming, along with flower-falling stages to detect genes that might be involved in response to drought stress. Finally, quantitative real-time PCR (qPCR) was carried out to analyze three candidate genes for their expression levels within 100 mM NaCl- and 10% PEG 6000 (w/v)-treated potato plantlets for a 24-h period. Furthermore, we analyzed the drought tolerance of StPUB25 transgenic plants and found that overexpression of StPUB25 significantly increased peroxidase (POD) activity, reduced ROS (reactive oxygen species) and MDA (malondialdehyde) accumulation compared with wild-type (WT) plants, and enhancing drought tolerance of the transgenic plants.

CONCLUSION

In this study, three candidate genes related to drought tolerance in potato were excavated, and the function of StPUB25 under drought stress was verified. These results should provide valuable information to understand the potato StPUB gene family and investigate the molecular mechanisms of StPUBs regulating potato drought tolerance.

Leveraging brassinosteroids towards the next Green Revolution

The use of gibberellin-related dwarfing genes significantly increased grain yield during the Green Revolution. Brassinosteroids (BRs) play a vital role in regulating agronomic traits and stress resistance. The potential of BR-related genes in crop improvement has been well demonstrated, positioning BRs as crucial targets for the next agricultural biotechnological revolution. However, BRs exert pleiotropic effects on plants, and thus present both opportunities and challenges for their application. Recent research suggests promising strategies for leveraging BR regulatory molecules for crop improvement, such as exploring function-specific genes, identifying beneficial alleles, inducing favorable mutations, and optimizing spatial hormone distribution. Advancing our understanding of the roles of BRs in plants is imperative to implement these strategies effectively.

Brassinosteroid signaling regulator BIM1 integrates brassinolide and jasmonic acid signaling during cold tolerance in apple

Although brassinolide (BR) and jasmonic acid (JA) play essential roles in the regulation of cold stress responses, the molecular basis of their crosstalk remains elusive. Here, we show a key component of BR signaling in apple (Malus × domestica), BR INSENSITIVE1 (BRI1)-EMS-SUPPRESSOR1 (BES1)-INTERACTING MYC-LIKE PROTEIN1 (MdBIM1), increases cold tolerance by directly activating expression of C-REPEAT BINDING FACTOR1 (MdCBF1) and forming a complex with C-REPEAT BINDING FACTOR2 (MdCBF2) to enhance MdCBF2-activated transcription of cold-responsive genes. Two repressors of JA signaling, JAZMONATE ZIM-DOMAIN1 (MdJAZ1) and JAZMONATE ZIM-DOMAIN2 (MdJAZ2), interact with MdBIM1 to integrate BR and JA signaling under cold stress. MdJAZ1 and MdJAZ2 reduce MdBIM1-promoted cold stress tolerance by attenuating transcriptional activation of MdCBF1 expression by MdBIM1 and interfering with the formation of the MdBIM1-MdCBF2 complex. Furthermore, the E3 ubiquitin ligase ARABIDOPSIS TÓXICOS en LEVADURA73 (MdATL73) decreases MdBIM1-promoted cold tolerance by targeting MdBIM1 for ubiquitination and degradation. Our results not only reveal crosstalk between BR and JA signaling mediated by a JAZ-BIM1-CBF module but also provide insights into the posttranslational regulatory mechanism of BR signaling.

Gene expression and expression quantitative trait loci analyses uncover natural variations underlying the improvement of important agronomic traits during modern maize breeding

Maize (Zea mays L.) is a major staple crop worldwide, and during modern maize breeding, cultivars with increased tolerance to high-density planting and higher yield per plant have contributed significantly to the increased yield per unit land area. Systematically identifying key agronomic traits and their associated genomic changes during modern maize breeding remains a significant challenge because of the complexity of genetic regulation and the interactions of the various agronomic traits, with most of them being controlled by numerous small-effect quantitative trait loci (QTLs). Here, we performed phenotypic and gene expression analyses for a set of 137 elite inbred lines of maize from different breeding eras in China. We found four yield-related traits are significantly improved during modern maize breeding. Through gene-clustering analyses, we identified four groups of expressed genes with distinct trends of expression pattern change across the historical breeding eras. In combination with weighted gene co-expression network analysis, we identified several candidate genes regulating various plant architecture- and yield-related agronomic traits, such as ZmARF16, ZmARF34, ZmTCP40, ZmPIN7, ZmPYL10, ZmJMJ10, ZmARF1, ZmSWEET15b, ZmGLN6 and Zm00001d019150. Further, by combining expression quantitative trait loci (eQTLs) analyses, correlation coefficient analyses and population genetics, we identified a set of candidate genes that might have been under selection and contributed to the genetic improvement of various agronomic traits during modern maize breeding, including a number of known key regulators of plant architecture, flowering time and yield-related traits, such as ZmPIF3.3, ZAG1, ZFL2 and ZmBES1. Lastly, we validated the functional variations in GL15, ZmPHYB2 and ZmPYL10 that influence kernel row number, flowering time, plant height and ear height, respectively. Our results demonstrates the effectiveness of our combined approaches for uncovering key candidate regulatory genes and functional variation underlying the improvement of important agronomic traits during modern maize breeding, and provide a valuable genetic resource for the molecular breeding of maize cultivars with tolerance for high-density planting.

Cosuppression of AtGELP22 and AtGELP23, two ubiquitinated target proteins of RING E3 ligase AtAIRP5, increases tolerance to drought stress in Arabidopsis

AtAIRP5 RING E3 ubiquitin ligase was recently identified as a positive regulator of the abscisic acid (ABA)-mediated drought stress response by stimulating the degradation of serine carboxypeptidase-like 1. Here, we identified GDSL-type esterase/lipase 22 (AtGELP22) and AtGELP23 as additional interacting partners of AtAIRP5. Yeast two-hybrid, pull-down, co-immunoprecipitation, and ubiquitination analyses verified that AtGELP22 and AtGELP23 are ubiquitinated target proteins of AtAIRP5. AtGELP22 and AtGELP23 were colocalized with AtAIRP5 to punctate-like structures in the cytosolic fraction, in which PYK10 and NAI2, two ER body marker proteins, are localized. T-DNA insertion atgelp22 and atgelp23 single knockout mutant plants showed phenotypes indistinguishable from those of wild-type plants under ABA treatment. In contrast, RNAi-mediated cosuppression of AtGELP22 and AtGELP23 resulted in hypersensitive ABA-mediated stomatal movements and higher tolerance to drought stress than that of the single mutant and wild-type plants. Taken together, our results suggest that the putative GDSL-type esterases/lipases AtGELP22 and AtGELP23 act as redundant negative regulators of the ABA-mediated drought stress response in Arabidopsis.

Phosphorylation and ubiquitination of OsWRKY31 are integral to OsMKK10-2-mediated defense responses in rice

Mitogen-activated protein kinase (MPK) cascades play vital roles in plant innate immunity, growth, and development. Here, we report that the rice (Oryza sativa) transcription factor gene OsWRKY31 is a key component in a MPK signaling pathway involved in plant disease resistance in rice. We found that the activation of OsMKK10-2 enhances resistance against the rice blast pathogen Magnaporthe oryzae and suppresses growth through an increase in jasmonic acid and salicylic acid accumulation and a decrease of indole-3-acetic acid levels. Knockout of OsWRKY31 compromises the defense responses mediated by OsMKK10-2. OsMKK10-2 and OsWRKY31 physically interact, and OsWRKY31 is phosphorylated by OsMPK3, OsMPK4, and OsMPK6. Phosphomimetic OsWRKY31 has elevated DNA-binding activity and confers enhanced resistance to M. oryzae. In addition, OsWRKY31 stability is regulated by phosphorylation and ubiquitination via RING-finger E3 ubiquitin ligases interacting with WRKY 1 (OsREIW1). Taken together, our findings indicate that modification of OsWRKY31 by phosphorylation and ubiquitination functions in the OsMKK10-2-mediated defense signaling pathway.

Analyses of genes encoding the Glycogen Synthase Kinases in rice and Arabidopsis reveal mechanisms which regulate their expression during development and responses to abiotic stresses

Plant Glycogen Synthase Kinases (GSKs) enable a crosstalk among the brassinosteroid signaling and phytohormonal- and stress-response pathways to regulate various physiological processes. Initial information about regulation of the GSK proteins' activity was obtained, however, mechanisms that modulate expression of the GSK genes during plant development and stress responses remain largely unknown. Taking into account the importance of the GSK proteins, combined with the lack of in-depth knowledge about modulation of their expression, research in this area may provide a significant insight into mechanisms regulating these aspects of plant biology. In the current study, a detailed analysis of the GSK promoters in rice and Arabidopsis was performed, including identification of the CpG/CpNpG islands, tandem repeats, cis-acting regulatory elements, conserved motifs, and transcription factor-binding sites. Moreover, characterization of expression profiles of the GSK genes in different tissues, organs and under various abiotic stress conditions was perfomed. Additionally, protein-protein interactions between products of the GSK genes were predicted. Results of this study provided intriguing information about these aspects and insight into various regulatory mechanisms that influence non-redundant and diverse functions of the GSK genes during development and stress responses.Therefore, they may constitute a reference for future research in other plant species.

Regulatory network of GSK3-like kinases and their role in plant stress response

Glycogen synthase kinase 3 (GSK3) family members are evolutionally conserved Ser/Thr protein kinases in mammals and plants. In plants, the GSK3s function as signaling hubs to integrate the perception and transduction of diverse signals required for plant development. Despite their role in the regulation of plant growth and development, emerging research has shed light on their multilayer function in plant stress responses. Here we review recent advances in the regulatory network of GSK3s and the involvement of GSK3s in plant adaptation to various abiotic and biotic stresses. We also discuss the molecular mechanisms underlying how plants cope with environmental stresses through GSK3s-hormones crosstalk, a pivotal biochemical pathway in plant stress responses. We believe that our overview of the versatile physiological functions of GSK3s and underlined molecular mechanism of GSK3s in plant stress response will not only opens further research on this important topic but also provide opportunities for developing stress-resilient crops through the use of genetic engineering technology.

OsARF4 regulates leaf inclination via auxin and brassinosteroid pathways in rice

Leaf inclination is a vital agronomic trait and is important for plant architecture that affects photosynthetic efficiency and grain yield. To understand the molecular mechanisms underlying regulation of leaf inclination, we constructed an auxin response factor (arf) rice mutant-osarf4-showing increased leaf inclination using CRISPR/Cas9 gene editing technology. OsARF4 encodes a nuclear protein that is expressed in the lamina joint (LJ) at different developmental stages in rice. Histological analysis indicated that an increase in cell differentiation on the adaxial side resulted in increased leaf inclination in the osarf4 mutants; however, OsARF4-overexpressing lines showed a decrease in leaf inclination, resulting in erect leaves. Additionally, a decrease in the content and distribution of indole-3-acetic acid (IAA) in osarf4 mutant led to a greater leaf inclination, whereas the OsARF4-overexpressing lines showed the opposite phenotype with increased IAA content. RNA-sequencing analysis revealed that the expression of genes related to brassinosteroid (BR) biosynthesis and response was different in the mutants and overexpressing lines, suggesting that OsARF4 participates in the BR signaling pathway. Moreover, BR sensitivity assay revealed that OsARF4-overexpressing lines were more sensitive to exogenous BR treatment than the mutants. In conclusion, OsARF4, a transcription factor in auxin signaling, participates in leaf inclination regulation and links auxin and BR signaling pathways. Our results provide a novel insight into l leaf inclination regulation, and have significant implications for improving rice architecture and grain yield.

BES1/BZR1 Family Transcription Factors Regulate Plant Development via Brassinosteroid-Dependent and Independent Pathways

The BES1/BZR1 family is a plant-specific small group of transcription factors possessing a non-canonical bHLH domain. Genetic and biochemical analyses within the last two decades have demonstrated that members of this family are key transcription factors in regulating the expression of brassinosteroid (BR) response genes. Several recent genetic and evolutionary studies, however, have clearly indicated that the BES1/BZR1 family transcription factors also function in regulating several aspects of plant development via BR-independent pathways, suggesting they are not BR specific. In this review, we summarize our current understanding of this family of transcription factors, the mechanisms regulating their activities, DNA binding motifs, and target genes. We selectively discuss a number of their biological functions via BR-dependent and particularly independent pathways, which were recently revealed by loss-of-function genetic analyses. We also highlight a few possible future directions.

E3 ligase AtAIRP5/GARU regulates drought stress response by stimulating SERINE CARBOXYPEPTIDASE-LIKE1 turnover

Ubiquitination is a major mechanism of eukaryotic posttranslational protein turnover that has been implicated in abscisic acid (ABA)-mediated drought stress response. Here, we isolated T-DNA insertion mutant lines in which ABA-insensitive RING protein 5 (AtAIRP5) was suppressed, resulting in hyposensitive ABA-mediated germination compared to wild-type Arabidopsis (Arabidopsis thaliana) plants. A homology search revealed that AtAIRP5 is identical to gibberellin (GA) receptor RING E3 ubiquitin (Ub) ligase (GARU), which downregulates GA signaling by degrading the GA receptor GID1, and thus AtAIRP5 was renamed AtAIRP5/GARU. The atairp5/garu knockout progeny were impaired in ABA-dependent stomatal closure and were markedly more susceptible to drought stress than wild-type plants, indicating a positive role for AtAIRP5/GARU in the ABA-mediated drought stress response. Yeast two-hybrid, pull-down, target ubiquitination, and in vitro and in planta degradation assays identified serine carboxypeptidase-like1 (AtSCPL1), which belongs to the clade 1A AtSCPL family, as a ubiquitinated target protein of AtAIRP5/GARU. atscpl1 single and atairp5/garu-1 atscpl1-2 double mutant plants were more tolerant to drought stress than wild-type plants in an ABA-dependent manner, suggesting that AtSCPL1 is genetically downstream of AtAIRP5/GARU. After drought treatment, the endogenous ABA levels in atscpl1 and atairp5/garu-1 atscpl1-2 mutant leaves were higher than those in wild-type and atairp5/garu leaves. Overall, our results suggest that AtAIRP5/GARU RING E3 Ub ligase functions as a positive regulator of the ABA-mediated drought response by promoting the degradation of AtSCPL1.

OsCPL3 is involved in brassinosteroid signaling by regulating OsGSK2 stability

Glycogen synthase kinase 3 (GSK3) proteins play key roles in brassinosteroid (BR) signaling during plant growth and development by phosphorylating various substrates. However, how GSK3 protein stability and activity are themselves modulated is not well understood. Here, we demonstrate in vitro and in vivo that C-TERMINAL DOMAIN PHOSPHATASE-LIKE 3 (OsCPL3), a member of the RNA Pol II CTD phosphatase-like family, physically interacts with OsGSK2 in rice (Oryza sativa). OsCPL3 expression was widely detected in various tissues and organs including roots, leaves and lamina joints, and was induced by exogenous BR treatment. OsCPL3 localized to the nucleus, where it dephosphorylated OsGSK2 at the Ser-222 and Thr-284 residues to modulate its protein turnover and kinase activity, in turn affecting the degradation of BRASSINAZOLE-RESISTANT 1 (BZR1) and BR signaling. Loss of OsCPL3 function resulted in higher OsGSK2 abundance and lower OsBZR1 levels, leading to decreased BR responsiveness and alterations in plant morphology including semi-dwarfism, leaf erectness and grain size, which are of fundamental importance to crop productivity. These results reveal a previously unrecognized role for OsCPL3 and add another layer of complexity to the tightly controlled BR signaling pathway in plants.

OsNAC016 regulates plant architecture and drought tolerance by interacting with the kinases GSK2 and SAPK8

Ideal plant architecture and drought tolerance are important determinants of yield potential in rice (Oryza sativa). Here, we found that OsNAC016, a rice NAC (NAM, ATAF, and CUC) transcription factor, functions as a regulator in the crosslink between brassinosteroid (BR)-mediated plant architecture and abscisic acid (ABA)-regulated drought responses. The loss-of-function mutant osnac016 exhibited erect leaves and shortened internodes, but OsNAC016-overexpressing plants had opposite phenotypes. Further investigation revealed that OsNAC016 regulated the expression of the BR biosynthesis gene D2 by binding to its promoter. Moreover, OsNAC016 interacted with and was phosphorylated by GSK3/SHAGGY-LIKE KINASE2 (GSK2), a negative regulator in the BR pathway. Meanwhile, the mutant osnac016 had improved drought stress tolerance, supported by a decreased water loss rate and enhanced stomatal closure in response to exogenous ABA, but OsNAC016-overexpressing plants showed attenuated drought tolerance and reduced ABA sensitivity. Further, OSMOTIC STRESS/ABA-ACTIVATED PROTEIN KINASE8 (SAPK8) phosphorylated OsNAC016 and reduced its stability. The ubiquitin/26S proteasome system is an important degradation pathway of OsNAC016 via the interaction with PLANT U-BOX PROTEIN43 (OsPUB43) that mediates the ubiquitination of OsNAC016. Notably, RNA-sequencing analysis revealed global roles of OsNAC016 in promoting BR-mediated gene expression and repressing ABA-dependent drought-responsive gene expression, which was confirmed by chromatin immunoprecipitation quantitative PCR analysis. Our findings establish that OsNAC016 is positively involved in BR-regulated rice architecture, negatively modulates ABA-mediated drought tolerance, and is regulated by GSK2, SAPK8, and OsPUB43 through posttranslational modification. Our data provide insights into how plants balance growth and survival by coordinately regulating the growth-promoting signaling pathway and response under abiotic stresses.

Evolution and Functions of Plant U-Box Proteins: From Protein Quality Control to Signaling

Posttranslational modifications add complexity and diversity to cellular proteomes. One of the most prevalent modifications across eukaryotes is ubiquitination, which is orchestrated by E3 ubiquitin ligases. U-box-containing E3 ligases have massively expanded in the plant kingdom and have diversified into plant U-box proteins (PUBs). PUBs likely originated from two or three ancestral forms, fusing with diverse functional subdomains that resulted in neofunctionalization. Their emergence and diversification may reflect adaptations to stress during plant evolution, reflecting changes in the needs of plant proteomes to maintain cellular homeostasis. Through their close association with protein kinases, they are physically linked to cell signaling hubs and activate feedback loops by dynamically pairing with E2-ubiquitin-conjugating enzymes to generate distinct ubiquitin polymers that themselves act as signals. Here, we complement current knowledgewith comparative genomics to gain a deeper understanding of PUB function, focusing on their evolution and structural adaptations of key U-box residues, as well as their various roles in plant cells.

Leaf angle: a target of genetic improvement in cereal crops tailored for high-density planting

High-density planting is an effective measure for increasing crop yield per unit land area. Leaf angle (LA) is a key trait of plant architecture and a target for genetic improvement of crops. Upright leaves allow better light capture in canopy under high-density planting, thus enhancing photosynthesis efficiency, ventilation and stress resistance, and ultimately higher grain yield. Here, we summarized the latest progress on the cellular and molecular mechanisms regulating LA formation in rice and maize. We suggest several standing out questions for future studies and then propose some promising strategies to manipulate LA for breeding of cereal crops tailored for high-density planting.

OsATL38 mediates mono-ubiquitination of the 14-3-3 protein OsGF14d and negatively regulates the cold stress response in rice

One of the major regulatory pathways that permits plants to convert an external stimulus into an internal cellular response within a short period of time is the ubiquitination pathway. In this study, OsATL38 was identified as a low temperature-induced gene that encodes a rice homolog of Arabidopsis Tóxicos en Levadura RING-type E3 ubiquitin (Ub) ligase, which was predominantly localized to the plasma membrane. OsATL38-overexpressing transgenic rice plants exhibited decreased tolerance to cold stress as compared with wild-type rice plants. In contrast, RNAi-mediated OsATL38 knockdown transgenic progeny exhibited markedly increased tolerance to cold stress relative to that of wild-type plants, which indicated a negative role of OsATL38 in response to cold stress. Yeast two-hybrid, in vitro pull-down, and co-immunoprecipitation assays revealed that OsATL38 physically interacted with OsGF14d, a rice 14-3-3 protein. An in vivo target ubiquitination assay indicated that OsGF14d was mono-ubiquitinated by OsATL38. osgf14d knockout mutant plants were more sensitive to cold stress than wild-type rice plants, indicating that OsGF14d is a positive factor in the response to cold stress. These results provide evidence that the RING E3 Ub ligase OsATL38 negatively regulates the cold stress response in rice via mono-ubiquitination of OsGF14d 14-3-3 protein.

Glycogen synthase kinases in model and crop plants - From negative regulators of brassinosteroid signaling to multifaceted hubs of various signaling pathways and modulators of plant reproduction and yield

Glycogen synthase kinases, also known as SHAGGY-like Kinases (GSKs/SKs), are highly conserved serine/threonine protein kinases present both in animals and plants. Plant genomes contain multiple homologs of the GSK3 genes which participate in various biological processes. Plant GSKs/SKs, and their best known representative in Arabidopsis thaliana - Brassinosteroid Insentisive2 (BIN2/SK21) in particular, were first identified as components of the brassinosteroid (BR) signaling pathway. As phytohormones, BRs regulate a wide range of physiological processes in plants - from germination, cell division, elongation and differentiation to leaf senescence, and response to environmental stresses. The GSKs/SKs proteins belong to a group of several highly conserved components of the BR signaling which evolved early during evolution of this molecular relay. However, recent reports indicated that the GSKs/SKs proteins are also implicated in signaling pathways of other phytohormones and stress-response processes. As a consequence, the GSKs/SKs proteins became hubs of various signaling pathways and modulators of plant development and reproduction. Thus, it is very important to understand molecular mechanisms regulating activity of the GSKs/SKs proteins, but also to get insights into role of the GSKs/SKs proteins in modulation of stability and activity of various substrate proteins which participate in the numerous signaling pathways. Although elucidation of these aspects is still in progress, this review presents a comprehensive and detailed description of these processes and their implications for regulation of development, stress response, and reproduction of model and crop species. The GSKs/SKs proteins and their activity are modulated through phosphorylation and de-phosphorylation reactions which are regulated by various proteins. Importantly, both phosphorylations and de-phosphorylations may have positive and negative effects on the activity of the GSKs/SKs proteins. Additionally, the activity of the GSKs/SKs proteins is positively regulated by reactive oxygen species, whereas it is negatively regulated through ubiquitylation, deacetylation, and nitric oxide-mediated nitrosylation. On the other hand, the GSKs/SKs proteins interact with proteins representing various signaling pathways, and on the basis of the complicated network of interactions the GSKs/SKs proteins differentially regulate various physiological, developmental, stress response, and yield-related processes.

Brassinazole Resistant 1 Activity Is Organ-Specific and Genotype-Dependent in Barley Seedlings

Brassinosteroids (BRs) control many plant developmental processes by regulating different groups of transcription factors, and consequently gene expressions. The most known is BZR1, the main member of the BES1 family. However, to date, it is poorly characterized in crop species. The main goal of the presented study was to identify HvBZR1 and determine its activity in 5-day-old barley (the stage is related to one leaf on the main shoot and a few seminal roots) using two cultivars with different sensitivities to BRs. Using the anti-OsBZR1 antibody, we identified the forms of HvBZR1 transcription factor with different molecular weights, which can be related to different phosphorylated forms of serine/threonine residues. Two phosphorylated forms in the shoots and one dephosphorylated form in the roots were determined. A minor amount of the dephosphorylated form of the HvBZR1 in the Haruna Nijo shoots was also found. The phosphorylated forms gave a higher band intensity for Golden Promise than Haruna Nijo. The bands were similar in their intensity, when two different phosphorylated forms were compared in Golden Promise, while a reduced intensity was detected for the phosphorylated form with a lower molecular weight for Haruna Nijo. Degradation of the phosphorylated forms in the shoots (complete degradation in Golden Promise and significant but not complete in Haruna Nijo) and the presence of the dephosphorylated form in the roots were proven for the etiolated barley. In the case of Haruna Nijo, a wider range of the regulators of the BR biosynthesis and signaling pathways induced the expected effects, 24-EBL (0.001 µM) and bikinin (10 and 50 µM) caused low amount of the phosphorylated forms, and at the same time, a tiny band of dephosphorylated form was detected. However, the expression of genes related to the BR biosynthesis and signaling pathways was not a determinant for the protein amount.

GSK3s: nodes of multilayer regulation of plant development and stress responses

Glycogen synthase kinase 3 (GSK3) family members are highly conserved serine/threonine protein kinases in eukaryotes. Unlike animals, plants have evolved with multiple homologs of GSK3s involved in a diverse array of biological processes. Emerging evidence suggests that GSK3s act as signaling hubs for integrating perception and transduction of diverse signals required for plant development and responses to abiotic and biotic cues. Here we review recent advances in understanding the molecular interactions between GSK3s and an expanding spectrum of their upstream regulators and downstream substrates in plants. We further discuss how GSK3s act as key signaling nodes of multilayer regulation of plant development and stress response through either being regulated at the post-translational level or regulating their substrates via phosphorylation.

The ubiquitin E3 ligase SR1 modulates the submergence response by degrading phosphorylated WRKY33 in Arabidopsis

Oxygen deprivation caused by flooding activates acclimation responses to stress and restricts plant growth. After experiencing flooding stress, plants must restore normal growth; however, which genes are dynamically and precisely controlled by flooding stress remains largely unknown. Here, we show that the Arabidopsis thaliana ubiquitin E3 ligase SUBMERGENCE RESISTANT1 (SR1) regulates the stability of the transcription factor WRKY33 to modulate the submergence response. SR1 physically interacts with WRKY33 in vivo and in vitro and controls its ubiquitination and proteasomal degradation. Both the sr1 mutant and WRKY33 overexpressors exhibited enhanced submergence tolerance and enhanced expression of hypoxia-responsive genes. Genetic experiments showed that WRKY33 functions downstream of SR1 during the submergence response. Submergence induced the phosphorylation of WRKY33, which enhanced the activation of RAP2.2, a positive regulator of hypoxia-response genes. Phosphorylated WRKY33 and RAP2.2 were degraded by SR1 and the N-degron pathway during reoxygenation, respectively. Taken together, our findings reveal that the on-and-off module SR1-WRKY33-RAP2.2 is connected to the well-known N-degron pathway to regulate acclimation to submergence in Arabidopsis. These two different but related modulation cascades precisely balance submergence acclimation with normal plant growth.

OsbHLH98 regulates leaf angle in rice through transcriptional repression of OsBUL1

Leaf angle is an important agronomic trait in cereals that helps determine plant yield by affecting planting density. However, the regulation mechanism of leaf angle remained elusive. Here, we show that OsbHLH98, a rice bHLH transcription factor, negatively regulates leaf angle. osbhlh98 mutant leaves formed a larger leaf angle, whereas transgenic plants overexpressing OsbHLH98 exhibited a slight reduction in leaf angle. We determined that the changes in leaf angle resulted from increased number and size of parenchyma cells on the adaxial side of the lamina joint in osbhlh98 mutants. Experiments using reporter constructs showed that OsbHLH98 is expressed on the adaxial side of lamina joints, consistent with its proposed function in regulating leaf angle. Furthermore, we established by chromatin immunoprecipitation and CUT&RUN that OsBUL1 is a direct downstream target of OsbHLH98. Transactivation assays and reverse transcription quantitative polymerase chain reaction (RT-qPCR) analysis indicated that OsbHLH98 represses OsBUL1 transcription. Our results demonstrate that OsbHLH98 negatively regulates leaf angle by counteracting brassinosteroid-induced cell elongation via the repression of OsBUL1 transcription. The characterization of OsbHLH98 and its role in determining leaf angle will lay the foundation to develop the ideal plant architecture for adaptation to high planting density.

Versatile Physiological Functions of Plant GSK3-Like Kinases

The plant glycogen synthase kinase 3 (GSK3)-like kinases are highly conserved protein serine/threonine kinases that are grouped into four subfamilies. Similar to their mammalian homologs, these kinases are constitutively active under normal growth conditions but become inactivated in response to diverse developmental and environmental signals. Since their initial discoveries in the early 1990s, many biochemical and genetic studies were performed to investigate their physiological functions in various plant species. These studies have demonstrated that the plant GSK3-like kinases are multifunctional kinases involved not only in a wide variety of plant growth and developmental processes but also in diverse plant stress responses. Here we summarize our current understanding of the versatile physiological functions of the plant GSK3-like kinases along with their confirmed and potential substrates.

AtKPNB1, an Arabidopsis importin-β protein, is downstream of the RING E3 ubiquitin ligase AtAIRP1 in the ABA-mediated drought stress response

AtKPNB1, an Arabidopsis importin-β protein, was regulated by AtAIRP1 E3 ubiquitin ligase, which intensified the ABA-mediated drought stress response. As an early step in the abscisic acid (ABA)-mediated drought response, the ABA signal is transduced into the nucleus, and thus the nuclear transport system is crucially involved in the drought stress response. AtKPNB1, an importin-β protein, which is a core component of nuclear transport, was previously reported to be a negative factor in the ABA-mediated drought stress response (Luo et al. Luo et al., Plant J 75:377-389, 2013). Here, we report that AtAIPR1, an Arabidopsis RING-type E3 ubiquitin (Ub) ligase, interacted with and ubiquitinated AtKPNB1. A null mutation of AtKPNB1 suppressed the ABA-insensitive germination phenotype of atairp1 mutant seedlings as compared to that of the wild-type plants. Furthermore, the ABA-insensitive stomatal closure and drought-susceptible phenotypes of atairp1 were rescued in atairp1atkpnb1 double mutant progeny, indicating that AtKPNB1 functions downstream of AtAIRP1. These data suggest that AtAIRP1 regulates the ABA-mediated drought response in Arabidopsis via ubiquitination of AtKPNB1.

Quantitative phosphoproteomics analysis reveals that protein modification and sugar metabolism contribute to sprouting in potato after BR treatment

Brassinosteroids (BRs), a new class of steroid hormones, are involved in the regulation of plant cell elongation and seed germination. Nevertheless, the molecular mechanism of the effect of BRs on tuber sprouting remains largely unknown. In this study, quantitative phosphoproteomics was employed to investigate the protein phosphorylation changes in sprouting induced by BRs. Our results showed that BRs accelerated the conversion of starch into soluble sugar in tubers. A functional enrichment cluster analysis suggested that the "amino acid metabolism pathway" was upregulated and that "plant hormone signal transduction and protein export" were downregulated. BR treatment also changed the phosphorylation of proteins involved in the BR, ABA, starch and sugar signal transduction pathways, such as serine/threonine-protein kinase (BSK), 14-3-3, alpha-glucan water dikinase (GWD), sucrose-phosphate synthase (SPS), sucrose synthase (SS) and alkaline/neutral invertase (A/N-INV). These results shed more light on the pattern of protein phosphorylation in BR promoting potato sprouting.

Exploring the Brassinosteroid Signaling in Monocots Reveals Novel Components of the Pathway and Implications for Plant Breeding

Brassinosteroids (BRs) are a class of steroidal phytohormones which are key regulators of diverse processes during whole life cycle of plants. Studies conducted in the dicot model species Arabidopsis thaliana have allowed identification and characterization of various components of the BR signaling. It is currently known that the BR signaling is interconnected at various stages with other phytohormonal and stress signaling pathways. It enables a rapid and efficient adaptation of plant metabolism to constantly changing environmental conditions. However, our knowledge about mechanism of the BR signaling in the monocot species is rather limited. Thus, identification of new components of the BR signaling in monocots, including cereals, is an ongoing process and has already led to identification of some monocot-specific components of the BR signaling. It is of great importance as disturbances in the BR signaling influence architecture of mutant plants, and as a consequence, the reaction to environmental conditions. Currently, the modulation of the BR signaling is considered as a target to enhance yield and stress tolerance in cereals, which is of particular importance in the face of global climate change.

Genome-Wide Distribution, Expression and Function Analysis of the U-Box Gene Family in Brassica oleracea L

The plant U-box (PUB) protein family plays an important role in plant growth and development. The U-box gene family has been well studied in Arabidopsis thaliana, Brassica rapa, rice, etc., but there have been no systematic studies in Brassica oleracea. In this study, we performed genome-wide identification and evolutionary analysis of the U-box protein family of B. oleracea. Firstly, based on the Brassica database (BRAD) and the Bolbase database, 99 Brassica oleracea PUB genes were identified and divided into seven groups (I-VII). The BoPUB genes are unevenly distributed on the nine chromosomes of B. oleracea, and there are tandem repeat genes, leading to family expansion from the A. thaliana genome to the B. oleracea genome. The protein interaction network, GO annotation, and KEGG pathway enrichment analysis indicated that the biological processes and specific functions of the BoPUB genes may mainly involve abiotic stress. RNA-seq transcriptome data of different pollination times revealed spatiotemporal expression specificity of the BoPUB genes. The differential expression profile was consistent with the results of RT-qPCR analysis. Additionally, a large number of pollen-specific cis-acting elements were found in promoters of differentially expressed genes (DEG), which verified that these significantly differentially expressed genes after self-pollination (SP) were likely to participate in the self-incompatibility (SI) process, including gene encoding ARC1, a well-known downstream protein of SI in B. oleracea. Our study provides valuable information indicating that the BoPUB genes participates not only in the abiotic stress response, but are also involved in pollination.