Redundancy, Feedback, and Robustness in the Arabidopsis thaliana BZR/BEH Gene Family

Genome-wide identification and expression analysis of the BZR gene family in Zanthoxylum armatum DC and functional analysis of ZaBZR1 in drought tolerance

MAIN CONCLUSION

In this study, six ZaBZRs were identified in Zanthoxylum armatum DC, and all the ZaBZRs were upregulated by abscisic acid (ABA) and drought. Overexpression of ZaBZR1 enhanced the drought tolerance of transgenic Nicotiana benthamian. Brassinosteroids (BRs) are a pivotal class of sterol hormones in plants that play a crucial role in plant growth and development. BZR (brassinazole resistant) is a crucial transcription factor in the signal transduction pathway of BRs. However, the BZR gene family members have not yet been identified in Zanthoxylum armatum DC. In this study, six members of the ZaBZR family were identified by bioinformatic methods. All six ZaBZRs exhibited multiple phosphorylation sites. Phylogenetic and collinearity analyses revealed a closest relationship between ZaBZRs and ZbBZRs located on the B subgenomes. Expression analysis revealed tissue-specific expression patterns of ZaBZRs in Z. armatum, and their promoter regions contained cis-acting elements associated with hormone response and stress induction. Additionally, all six ZaBZRs showed upregulation upon treatment after abscisic acid (ABA) and polyethylene glycol (PEG), indicating their participation in drought response. Subsequently, we conducted an extensive investigation of ZaBZR1. ZaBZR1 showed the highest expression in the root, followed by the stem and terminal bud. Subcellular localization analysis revealed that ZaBZR1 is present in the cytoplasm and nucleus. Overexpression of ZaBZR1 in transgenic Nicotiana benthamiana improved seed germination rate and root growth under drought conditions, reducing water loss rates compared to wild-type plants. Furthermore, ZaBZR1 increased proline content (PRO) and decreased malondialdehyde content (MDA), indicating improved tolerance to drought-induced oxidative stress. The transgenic plants also showed a reduced accumulation of reactive oxygen species. Importantly, ZaBZR1 up-regulated the expression of drought-related genes such as NbP5CS1, NbDREB2A, and NbWRKY44. These findings highlight the potential of ZaBZR1 as a candidate gene for enhancing drought resistance in transgenic N. benthamiana and provide insight into the function of ZaBZRs in Z. armatum.

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.

Generation and Characterization of Engineered Ubiquitin Variants to Modulate the Ubiquitin Signaling Cascade

The ubiquitin signaling cascade plays a crucial role in human cells. Consistent with this, malfunction of ubiquitination and deubiquitination is implicated in the initiation and progression of numerous human diseases, including cancer. Therefore, the development of potent and specific modulators of ubiquitin signal transduction has been at the forefront of drug development. In the past decade, a structure-based combinatorial protein-engineering approach has been used to generate ubiquitin variants (UbVs) as protein-based modulators of multiple components in the ubiquitin-proteasome system. Here, we review the design and generation of phage-displayed UbV libraries, including the processes of binder selection and library improvement. We also provide a comprehensive overview of the general in vitro and cellular methodologies involved in characterizing UbV binders. Finally, we describe two recent applications of UbVs for developing molecules with therapeutic potential.

Identification and Characterization of the BZR Transcription Factor Genes Family in Potato (Solanum tuberosum L.) and Their Expression Profiles in Response to Abiotic Stresses

Brassinazole resistant (BZR) genes act downstream of the brassinosteroid signaling pathway regulating plant growth and development and participating in plant stress responses. However, the BZR gene family has not systematically been characterized in potato. We identified eight BZR genes in Solanum tuberosum, which were distributed among seven chromosomes unequally and were classified into three subgroups. Potato and tomato BZR proteins were shown to be closely related with high levels of similarity. The BZR gene family members in each subgroup contained similar conserved motifs. StBZR genes exhibited tissue-specific expression patterns, suggesting their functional differentiation during evolution. StBZR4, StBZR7, and StBZR8 were highly expressed under white light in microtubers. StBZR1 showed a progressive up-regulation from 0 to 6 h and a progressive down-regulation from 6 to 24 h after drought and salt stress. StBZR1, StBZR2, StBZR4, StBZR5, StBZR6, StBZR7 and StBZR8 were significantly induced from 0 to 3 h under BR treatment. This implied StBZR genes are involved in phytohormone and stress response signaling pathways. Our results provide a theoretical basis for understanding the functional mechanisms of BZR genes in potato.

LTP2 hypomorphs show genotype-by-environment interaction in early seedling traits in Arabidopsis thaliana

Isogenic individuals can display seemingly stochastic phenotypic differences, limiting the accuracy of genotype-to-phenotype predictions. The extent of this phenotypic variation depends in part on genetic background, raising questions about the genes involved in controlling stochastic phenotypic variation. Focusing on early seedling traits in Arabidopsis thaliana, we found that hypomorphs of the cuticle-related gene LIPID TRANSFER PROTEIN 2 (LTP2) greatly increased variation in seedling phenotypes, including hypocotyl length, gravitropism and cuticle permeability. Many ltp2 hypocotyls were significantly shorter than wild-type hypocotyls while others resembled the wild-type. Differences in epidermal properties and gene expression between ltp2 seedlings with long and short hypocotyls suggest a loss of cuticle integrity as the primary determinant of the observed phenotypic variation. We identified environmental conditions that reveal or mask the increased variation in ltp2 hypomorphs and found that increased expression of its closest paralog LTP1 is necessary for ltp2 phenotypes. Our results illustrate how decreased expression of a single gene can generate starkly increased phenotypic variation in isogenic individuals in response to an environmental challenge.

The caleosin RD20/CLO3 regulates lateral root development in response to abscisic acid and regulates flowering time in conjunction with the caleosin CLO7

The caleosins are encoded by multi-gene families in Arabidopsis thaliana and other plant species. This work investigates the role of two family members, RD20/CLO3 and CLO7, in flowering transition and in root development in response to ABA treatment. Gene expression of the caleosin RD20/CLO3 is induced by ABA in the root tissues and RD20/CLO3 has a negative affect on the total number of lateral roots as well as the length of the lateral roots in response to ABA treatment. The rd20/clo3 mutant has more and longer lateral roots in response to ABA treatment compared to the wild-type, showing that RD20/CLO3 plays a role in the ABA signaling pathway affecting this trait. In contrast, the caleosin CLO7 is not expressed in the roots and does not affect root architecture in response to ABA treatment. The disruption of both RD20/CLO3 and CLO7 together causes a dramatic early-flowering phenotype under long-day conditions, whereas single mutations in these genes do not affect flowering time under these conditions. Both yeast two-hybrid and bimolecular fluorescence complementation showed that both RD20/CLO3 and CLO7 interact with each other and can form homodimers and heterodimers. Taken together, these findings suggest that members of the caleosin gene family play both different and redundant roles in plant development.

LTP2 hypomorphs show genotype-by-environment interaction in early seedling traits in Arabidopsis thaliana

Isogenic individuals can display seemingly stochastic phenotypic differences, limiting the accuracy of genotype-to-phenotype predictions. The extent of this phenotypic variation depends in part on genetic background, raising questions about the genes involved in controlling stochastic phenotypic variation. Focusing on early seedling traits in Arabidopsis thaliana, we found that hypomorphs of the cuticle-related gene LTP2 greatly increased variation in seedling phenotypes, including hypocotyl length, gravitropism and cuticle permeability. Many ltp2 hypocotyls were significantly shorter than wild-type hypocotyls while others resembled the wild type. Differences in epidermal properties and gene expression between ltp2 seedlings with long and short hypocotyls suggest a loss of cuticle integrity as the primary determinant of the observed phenotypic variation. We identified environmental conditions that reveal or mask the increased variation in ltp2 hypomorphs, and found that increased expression of its closest paralog LTP1 is necessary for ltp2 phenotypes. Our results illustrate how decreased expression of a single gene can generate starkly increased phenotypic variation in isogenic individuals in response to an environmental challenge.

Plant biochemical genetics in the multiomics era

Our understanding of plant biology has been revolutionized by modern genetics and biochemistry. However, biochemical genetics can be traced back to the foundation of Mendelian genetics; indeed, one of Mendel's milestone discoveries of seven characteristics of pea plants later came to be ascribed to a mutation in a starch branching enzyme. Here, we review both current and historical strategies for the elucidation of plant metabolic pathways and the genes that encode their component enzymes and regulators. We use this historical review to discuss a range of classical genetic phenomena including epistasis, canalization, and heterosis as viewed through the lens of contemporary high-throughput data obtained via the array of approaches currently adopted in multiomics studies.

BRI1 EMS SUPPRESSOR1 genes regulate abiotic stress and anther development in wheat (Triticum aestivum L.)

BRI1 EMS SUPPRESSOR1 (BES1) family members are crucial downstream regulators that positively mediate brassinosteroid signaling, playing vital roles in the regulation of plant stress responses and anther development in Arabidopsis. Importantly, the expression profiles of wheat (Triticum aestivum L.) BES1 genes have not been analyzed comprehensively and systematically in response to abiotic stress or during anther development. In this study, we identified 23 BES1-like genes in common wheat, which were unevenly distributed on 17 out of 21 wheat chromosomes. Phylogenetic analysis clustered the BES1 genes into four major clades; moreover, TaBES1-3A2, TaBES1-3B2 and TaBES1-3D2 belonged to the same clade as Arabidopsis BES1/BZR1 HOMOLOG3 (BEH3) and BEH4, which participate in anther development. The expression levels of 23 wheat BES1 genes were assessed using real-time quantitative PCR under various abiotic stress conditions (drought, salt, heat, and cold), and we found that most TaBES1-like genes were downregulated under abiotic stress, particularly during drought stress. We therefore used drought-tolerant and drought-sensitive wheat cultivars to explore TaBES1 expression patterns under drought stress. TaBES1-3B2 and TaBES1-3D2 expression was high in drought-tolerant cultivars but substantially repressed in drought-sensitive cultivars, while TaBES1-6D presented an opposite pattern. Among genes preferentially expressed in anthers, TaBES1-3B2 and TaBES1-3D2 expression was substantially downregulated in thermosensitive genic male-sterile wheat lines compared to common wheat cultivar under sterile conditions, while we detected no obvious differences under fertile conditions. This result suggests that TaBES1-3B2 and TaBES1-3D2 might not only play roles in regulating drought tolerance, but also participate in low temperature-induced male sterility.

Genome-wide identification and expression analysis of BZR gene family and associated responses to abiotic stresses in cucumber (Cucumis sativus L.)

BACKGROUND

BRASSINAZOLE-RESISTANT (BZR) is a class of specific transcription factor (TFs) involved in brassinosteroid (BR) signal transduction. The regulatory mechanism of target genes mediated by BZR has become one of the key research areas in plant BR signaling networks. However, the functions of the BZR gene family in cucumber have not been well characterized.

RESULTS

In this study, six CsBZR gene family members were identified by analyzing the conserved domain of BES1 N in the cucumber genome. The size of CsBZR proteins ranges from 311 to 698 amino acids and are mostly located in the nucleus. Phylogenetic analysis divided CsBZR genes into three subgroups. The gene structure and conserved domain showed that the BZR genes domain in the same group was conserved. Cis-acting element analysis showed that cucumber BZR genes were mainly involved in hormone response, stress response and growth regulation. The qRT-PCR results also confirmed CsBZR response to hormones and abiotic stress.

CONCLUSION

Collectively, the CsBZR gene is involved in regulating cucumber growth and development, particularly in hormone response and response to abiotic stress. These findings provide valuable information for understanding the structure and expression patterns of BZR genes.

MBNL1 regulates programmed postnatal switching between regenerative and differentiated cardiac states

Discovering determinants of cardiomyocyte maturity and the maintenance of differentiated states is critical to both understanding development and potentially reawakening endogenous regenerative programs in adult mammalian hearts as a therapeutic strategy. Here, the RNA binding protein Muscleblind-like 1 (MBNL1) was identified as a critical regulator of cardiomyocyte differentiated states and their regenerative potential through transcriptome-wide control of RNA stability. Targeted MBNL1 overexpression early in development prematurely transitioned cardiomyocytes to hypertrophic growth, hypoplasia, and dysfunction, whereas loss of MBNL1 function increased cardiomyocyte cell cycle entry and proliferation through altered cell cycle inhibitor transcript stability. Moreover, MBNL1-dependent stabilization of the estrogen-related receptor signaling axis was essential for maintaining cardiomyocyte maturity. In accordance with these data, modulating MBNL1 dose tuned the temporal window of cardiac regeneration, where enhanced MBNL1 activity arrested myocyte proliferation, and MBNL1 deletion promoted regenerative states with prolonged myocyte proliferation. Collectively these data suggest MBNL1 acts as a transcriptome-wide switch between regenerative and mature myocyte states postnatally and throughout adulthood.

AGO1 and HSP90 buffer different genetic variants in Arabidopsis thaliana

Argonaute 1 (AGO1), the principal protein component of microRNA-mediated regulation, plays a key role in plant growth and development. AGO1 physically interacts with the chaperone HSP90, which buffers cryptic genetic variation in plants and animals. We sought to determine whether genetic perturbation of AGO1 in Arabidopsis thaliana would also reveal cryptic genetic variation, and if so, whether AGO1-dependent loci overlap with those dependent on HSP90. To address these questions, we introgressed a hypomorphic mutant allele of AGO1 into a set of mapping lines derived from the commonly used Arabidopsis strains Col-0 and Ler. Although we identified several cases in which AGO1 buffered genetic variation, none of the AGO1-dependent loci overlapped with those buffered by HSP90 for the traits assayed. We focused on 1 buffered locus where AGO1 perturbation uncoupled the traits days to flowering and rosette leaf number, which are otherwise closely correlated. Using a bulk segregant approach, we identified a nonfunctional Ler hua2 mutant allele as the causal AGO1-buffered polymorphism. Introduction of a nonfunctional hua2 allele into a Col-0 ago1 mutant background recapitulated the Ler-dependent ago1 phenotype, implying that coupling of these traits involves different molecular players in these closely related strains. Taken together, our findings demonstrate that even though AGO1 and HSP90 buffer genetic variation in the same traits, these robustness regulators interact epistatically with different genetic loci, suggesting that higher-order epistasis is uncommon. Plain Language Summary Argonaute 1 (AGO1), a key player in plant development, interacts with the chaperone HSP90, which buffers environmental and genetic variation. We found that AGO1 buffers environmental and genetic variation in the same traits; however, AGO1-dependent and HSP90-dependent loci do not overlap. Detailed analysis of a buffered locus found that a nonfunctional HUA2 allele decouples days to flowering and rosette leaf number in an AGO1-dependent manner, suggesting that the AGO1-dependent buffering acts at the network level.

Far from home: the role of glial mRNA localization in synaptic plasticity

Neurons and glia are highly polarized cells, whose distal cytoplasmic functional subdomains require specific proteins. Neurons have axonal and dendritic cytoplasmic extensions containing synapses whose plasticity is regulated efficiently by mRNA transport and localized translation. The principles behind these mechanisms are equally attractive for explaining rapid local regulation of distal glial cytoplasmic projections, independent of their cell nucleus. However, in contrast to neurons, mRNA localization has received little experimental attention in glia. Nevertheless, there are many functionally diverse glial subtypes containing extensive networks of long cytoplasmic projections with likely localized regulation that influence neurons and their synapses. Moreover, glia have many other neuron-like properties, including electrical activity, secretion of gliotransmitters and calcium signaling, influencing, for example, synaptic transmission, plasticity and axon pruning. Here, we review previous studies concerning glial transcripts with important roles in influencing synaptic plasticity, focusing on a few cases involving localized translation. We discuss a variety of important questions about mRNA transport and localized translation in glia that remain to be addressed, using cutting-edge tools already available for neurons.

A glutathione-independent DJ-1/PfpI domain-containing tomato glyoxalaseIII2, SlGLYIII2, confers enhanced tolerance under salt and osmotic stresses

In plants, glyoxalase enzymes are activated under stress conditions to mitigate the toxic effects of hyperaccumulated methylglyoxal (MG), a highly reactive carbonyl compound. Until recently, a glutathione-dependent bi-enzymatic pathway involving glyoxalase I (GLYI) and glyoxalase II (GLYII) was considered the primary MG-detoxification system. Recently, a new glutathione-independent glyoxalase III (GLYIII) mediated direct route was also reported in plants. However, the physiological significance of this new pathway remains to be elucidated across plant species. This study identified the full complement of 22 glyoxalases in tomato. Based on their strong induction under multiple abiotic stresses, SlGLYI4, SlGLYII2 and SlGLYIII2 were selected candidates for further functional characterisation. Stress-inducible overexpression of both glutathione-dependent (SlGLYI4 + SlGLYII2) and independent (SlGLYIII2) pathways led to enhanced tolerance in both sets of transgenic plants under abiotic stresses. However, SlGLYIII2 overexpression (OE) plants outperformed the SlGLYI4 + SlGLYII2 OE counterparts for their stress tolerance under abiotic stresses. Further, knockdown of SlGLYIII2 resulted in plants with exacerbated stress responses than those silenced for both SlGLYI4 and SlGLYII2. The superior performance of SlGLYIII2 OE tomato plants for better growth and yield under salt and osmotic treatments could be attributed to better GSH/GSSG ratio, lower reactive oxygen species levels, and enhanced antioxidant potential, indicating a prominent role of GLYIII MG-detoxification pathway in abiotic stress mitigation in this species.

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.

Genome-Wide Identification of Brassicaceae Hormone-Related Transcription Factors and Their Roles in Stress Adaptation and Plant Height Regulation in Allotetraploid Rapeseed

Phytohormone-related transcription factors (TFs) are involved in regulating stress responses and plant growth. However, systematic analysis of these TFs in Brassicaceae is limited, and their functions in stress adaptation and plant height (PH) regulation remain unclear. In this study, 2115 hormone-related TFs were identified in nine Brassicaceae species. Specific domains were found in several Brassicaceae hormone-related TFs, which may be associated with diverse functions. Syntenic analysis indicated that expansion of these genes was mainly caused by segmental duplication, with whole-genome duplication occurring in some species. Differential expression analysis and gene co-expression network analysis identified seven phytohormone-related TFs (BnaWRKY7, 21, 32, 38, 52, BnaGL3-4, and BnaAREB2-5) as possible key genes for cadmium (Cd) toxicity, salinity stress, and potassium (K) and nitrogen (N) deficiencies. Furthermore, BnaWRKY42 and BnaARR21 may play essential roles in plant height. Weighted gene co-expression network analysis (WGCNA) identified 15 phytohormone-related TFs and their potential target genes regulating stress adaptation and plant height. Among the above genes, BnaWRKY56 and BnaWRKY60 responded to four different stresses simultaneously, and BnaWRKY42 was identified in two dwarf rapeseeds. In summary, several candidate genes for stress resistance (BnaWRKY56 and BnaWRKY60) and plant height (BnaWRKY42) were identified. These findings should help elucidate the biological roles of Brassicaceae hormone-related TFs, and the identified candidate genes should provide a genetic resource for the potential development of stress-tolerant and dwarf oilseed plants.

GhBZR3 suppresses cotton fiber elongation by inhibiting very-long-chain fatty acid biosynthesis

The BRASSINAZOLE-RESISTANT (BZR) transcription factor is a core component of brassinosteroid (BR) signaling and is involved in the development of many plant species. BR is essential for the initiation and elongation of cotton fibers. However, the mechanism of BR-regulating fiber development and the function of BZR is poorly understood in Gossypium hirsutum L. (cotton). Here, we identified a BZR family transcription factor protein referred to as GhBZR3 in cotton. Overexpression of GhBZR3 in Arabidopsis caused shorter root hair length, hypocotyl length, and hypocotyl cell length, indicating that GhBZR3 negatively regulates cell elongation. Pathway enrichment analysis from VIGS-GhBZR3 cotton plants found that fatty acid metabolism and degradation might be the regulatory pathway that is primarily controlled by GhBZR3. Silencing GhBZR3 expression in cotton resulted in taller plant height as well as longer fibers. The very-long-chain fatty acid (VLCFA) content was also significantly increased in silenced GhBZR3 plants compared with the wild type. The GhKCS13 promoter, a key gene for VLCFA biosynthesis, contains two GhBZR3 binding sites. The results of yeast one-hybrid, electrophoretic mobility shift, and luciferase assays revealed that GhBZR3 directly interacted with the GhKCS13 promoter to suppress gene expression. Taken together, these results indicate that GhBZR3 negatively regulates cotton fiber development by reducing VLCFA biosynthesis. This study not only deepens our understanding of GhBZR3 function in cotton fiber development, but also highlights the potential of improving cotton fiber length and plant growth using GhBZR3 and its related genes in future cotton breeding programs.

Organ-specific COP1 control of BES1 stability adjusts plant growth patterns under shade or warmth

Under adverse conditions such as shade or elevated temperatures, cotyledon expansion is reduced and hypocotyl growth is promoted to optimize plant architecture. The mechanisms underlying the repression of cotyledon cell expansion remain unknown. Here, we report that the nuclear abundance of the BES1 transcription factor decreased in the cotyledons and increased in the hypocotyl in Arabidopsis thaliana under shade or warmth. Brassinosteroid levels did not follow the same trend. PIF4 and COP1 increased their nuclear abundance in both organs under shade or warmth. PIF4 directly bound the BES1 promoter to enhance its activity but indirectly reduced BES1 expression. COP1 physically interacted with the BES1 protein, promoting its proteasome degradation in the cotyledons. COP1 had the opposite effect in the hypocotyl, demonstrating organ-specific regulatory networks. Our work indicates that shade or warmth reduces BES1 activity by transcriptional and post-translational regulation to inhibit cotyledon cell expansion.

Genome-Wide Prediction, Functional Divergence, and Characterization of Stress-Responsive BZR Transcription Factors in B. napus

BRASSINAZOLE RESISTANT (BZR) are transcriptional factors that bind to the DNA of targeted genes to regulate several plant growth and physiological processes in response to abiotic and biotic stresses. However, information on such genes in Brassica napus is minimal. Furthermore, the new reference Brassica napus genome offers an excellent opportunity to systematically characterize this gene family in B. napus. In our study, 21 BnaBZR genes were distributed across 19 chromosomes of B. napus and clustered into four subgroups based on Arabidopsis thaliana orthologs. Functional divergence analysis among these groups evident the shifting of evolutionary rate after the duplication events. In terms of structural analysis, the BnaBZR genes within each subgroup are highly conserved but are distinctive within groups. Organ-specific expression analyses of BnaBZR genes using RNA-seq data and quantitative real-time polymerase chain reaction (qRT-PCR) revealed complex expression patterns in plant tissues during stress conditions. In which genes belonging to subgroups III and IV were identified to play central roles in plant tolerance to salt, drought, and Sclerotinia sclerotiorum stress. The insights from this study enrich our understanding of the B. napus BZR gene family and lay a foundation for future research in improving rape seed environmental adaptability.

Transcription factor BES1 interacts with HSFA1 to promote heat stress resistance of plants

Heat stress is a major environmental stress type that can limit plant growth and development. To survive sudden temperature increases, plants utilize the heat shock response, an ancient signaling pathway. Initial results had suggested a role for brassinosteroids (BRs) in this response. Brassinosteroids are growth-promoting steroid hormones whose activity is mediated by transcription factors of the BES1/BZR1 subfamily. Here, we provide evidence that BES1 can contribute to heat stress signaling. In response to heat, BES1 is activated even in the absence of BRs and directly binds to heat shock elements (HSEs), known binding sites of heat shock transcription factors (HSFs). HSFs of the HSFA1 type can interact with BES1 and facilitate its activity in HSE binding. These findings lead us to propose an extended model of the heat stress response in plants, in which the recruitment of BES1 is a means of heat stress signaling cross-talk with a central growth regulatory pathway.

Brassinosteroids regulate petal spur length in Aquilegia by controlling cell elongation

BACKGROUND AND AIMS

Aquilegia produce elongated, three-dimensional petal spurs that fill with nectar to attract pollinators. Previous studies have shown that the diversity of spur length across the Aquilegia genus is a key innovation that is tightly linked with its recent and rapid diversification into new ranges, and that evolution of increased spur lengths is achieved via anisotropic cell elongation. Previous work identified a brassinosteroid response transcription factor as being enriched in the early developing spur cup. Brassinosteroids are known to be important for cell elongation, suggesting that brassinosteroid-mediated response may be an important regulator of spur elongation and potentially a driver of spur length diversity in Aquilegia. In this study, we investigated the role of brassinosteroids in the development of the Aquilegia coerulea petal spur.

METHODS

We exogenously applied the biologically active brassinosteroid brassinolide to developing petal spurs to investigate spur growth under high hormone conditions. We used virus-induced gene silencing and gene expression experiments to understand the function of brassinosteroid-related transcription factors in A. coerulea petal spurs.

KEY RESULTS

We identified a total of three Aquilegia homologues of the BES1/BZR1 protein family and found that these genes are ubiquitously expressed in all floral tissues during development, yet, consistent with the previous RNAseq study, we found that two of these paralogues are enriched in early developing petals. Exogenously applied brassinosteroid increased petal spur length due to increased anisotropic cell elongation as well as cell division. We found that targeting of the AqBEH genes with virus-induced gene silencing resulted in shortened petals, a phenotype caused in part by a loss of cell anisotropy.

CONCLUSIONS

Collectively, our results support a role for brassinosteroids in anisotropic cell expansion in Aquilegia petal spurs and highlight the brassinosteroid pathway as a potential player in the diversification of petal spur length in Aquilegia.

Gene co-expression network analysis identifies BEH3 as a stabilizer of secondary vascular development in Arabidopsis

In plants, vascular stem cells located in the cambium continuously undergo self-renewal and differentiation during secondary growth. Recent advancements in cell sorting techniques have enabled access to the transcriptional regulatory framework of cambial cells. However, mechanisms underlying the robust control of vascular stem cells remain unclear. Here, we identified a new cambium-related regulatory module through co-expression network analysis using multiple transcriptome datasets obtained from an ectopic vascular cell transdifferentiation system using Arabidopsis cotyledons, Vascular cell Induction culture System Using Arabidopsis Leaves (VISUAL). The cambium gene list included a gene encoding the transcription factor BES1/BZR1 Homolog 3 (BEH3), whose homolog BES1 negatively affects vascular stem cell maintenance. Interestingly, null beh3 mutant alleles showed a large variation in their vascular size, indicating that BEH3 functions as a stabilizer of vascular stem cells. Genetic analysis revealed that BEH3 and BES1 perform opposite functions in the regulation of vascular stem cells and the differentiation of vascular cells in the context of the VISUAL system. At the biochemical level, BEH3 showed weak transcriptional repressor activity and functioned antagonistically to other BES/BZR members by competing for binding to the brassinosteroid response element. Furthermore, mathematical modeling suggested that the competitive relationship between BES/BZR homologs leads to the robust regulation of vascular stem cells.

Topological clustering of regulatory genes confers pathogenic tolerance to cassava brown streak virus (CBSV) in cassava

Robustness, a naïve property of biological systems, enables organisms to maintain functions during perturbation and is crucial for improving the resilience of crops to prevailing stress conditions and diseases, guaranteeing food security. Most studies of robustness in crops have focused on genetic superiority based upon individual genes, overlooking the collaborative actions of multiple responsive genes and the regulatory network topology. This research aims to uncover patterns of gene cooperation leading to organismal robustness by studying the topology of gene co-expression networks (GCNs) of both CBSV virus resistant and susceptible cassava cultivars. The resulting GCNs show higher topological clustering of cooperative genes in the resistant cultivar, suggesting that the network architecture is central to attaining robustness. Despite a reduction in the number of hub genes in the resistant cultivar following the perturbation, essential biological functions contained in the network were maintained through neighboring genes that withstood the shock. The susceptible cultivar seemingly coped by inducing more gene actions in the network but could not maintain the functions required for plant growth. These findings underscore the importance of regulatory network architecture in ensuring phenotypic robustness and deepen our understanding of transcriptional regulation.

Integrative inference of transcriptional networks in Arabidopsis yields novel ROS signalling regulators

Gene regulation is a dynamic process in which transcription factors (TFs) play an important role in controlling spatiotemporal gene expression. To enhance our global understanding of regulatory interactions in Arabidopsis thaliana, different regulatory input networks capturing complementary information about DNA motifs, open chromatin, TF-binding and expression-based regulatory interactions were combined using a supervised learning approach, resulting in an integrated gene regulatory network (iGRN) covering 1,491 TFs and 31,393 target genes (1.7 million interactions). This iGRN outperforms the different input networks to predict known regulatory interactions and has a similar performance to recover functional interactions compared to state-of-the-art experimental methods. The iGRN correctly inferred known functions for 681 TFs and predicted new gene functions for hundreds of unknown TFs. For regulators predicted to be involved in reactive oxygen species (ROS) stress regulation, we confirmed in total 75% of TFs with a function in ROS and/or physiological stress responses. This includes 13 ROS regulators, previously not connected to any ROS or stress function, that were experimentally validated in our ROS-specific phenotypic assays of loss- or gain-of-function lines. In conclusion, the presented iGRN offers a high-quality starting point to enhance our understanding of gene regulation in plants by integrating different experimental data types.

BES1/BZR1 Homolog 3 cooperates with E3 ligase AtRZF1 to regulate osmotic stress and brassinosteroid responses in Arabidopsis

Proline (Pro) metabolism plays important roles in protein synthesis, redox balance, and abiotic stress response. However, it is not known if cross-talk occurs between proline and brassinosteroid (BR) signaling pathways. Here, an Arabidopsis intergenic enhancer double mutant, namely proline content alterative 41 (pca41), was generated by inserting a T-DNA tag in the Arabidopsis thaliana ring zinc finger 1 (atrzf1 ) mutant background. pca41 had a T-DNA inserted at the site of the gene encoding BES1/BZR1 Homolog 3 (BEH3). pca41 has a drought-insensitive phenotype that is stronger than atrzf1 under osmotic stress, including high Pro accumulation and decreased amounts of reactive oxygen species. Analysis of physiological, genetic, and molecular networks revealed that negative regulation of BEH3 during abiotic stress was linked to the BR signaling pathway. Our data also suggest that AtRZF1, an E3 ubiquitin ligase, might control osmotic stress, abscisic acid, and BR responses in a BEH3-dependent manner. Under darkness, pca41 displays a long hypocotyl phenotype, which is similar to atrzf1 and beh3, suggesting that BEH3 acts in the same pathway as AtRZF1. Overexpression of BEH3 results in an osmotic stress-sensitive phenotype, which is reversed by exogenous BR application. Taken together, our results indicate that AtRZF1 and BEH3 may play important roles in the osmotic stress response via ubiquitination and BR signaling.

Evolution and Plasticity of the Transcriptome Under Temperature Fluctuations in the Fungal Plant Pathogen Zymoseptoria tritici

Most species live in a variable environment in nature. Yet understanding the evolutionary processes underlying molecular adaptation to fluctuations remains a challenge. In this study we investigate the transcriptome of the fungal wheat pathogen Zymoseptoria tritici after experimental evolution under stable or fluctuating temperature, by comparing ancestral and evolved populations simultaneously. We found that temperature regimes could have a large and pervasive effect on the transcriptome evolution, with as much as 38% of the genes being differentially expressed between selection regimes. Although evolved lineages showed different changes of gene expression based on ancestral genotypes, we identified a set of genes responding specifically to fluctuation. We found that transcriptome evolution in fluctuating conditions was repeatable between parallel lineages initiated from the same genotype for about 60% of the differentially expressed genes. Further, we detected several hotspots of significantly differentially expressed genes in the genome, in regions known to be enriched in repetitive elements, including accessory chromosomes. Our findings also evidenced gene expression evolution toward a gain of robustness (loss of phenotypic plasticity) associated with the fluctuating regime, suggesting robustness is adaptive in changing environment. This work provides valuable insight into the role of transcriptional rewiring for rapid adaptation to abiotic changes in filamentous plant pathogens.

Multifaceted roles of HEAT SHOCK PROTEIN 90 molecular chaperones in plant development

HEAT SHOCK PROTEINS 90 (HSP90s) are molecular chaperones that mediate correct folding and stability of many client proteins. These chaperones act as master molecular hubs involved in multiple aspects of cellular and developmental signalling in diverse organisms. Moreover, environmental and genetic perturbations affect both HSP90s and their clients, leading to alterations of molecular networks determining respectively plant phenotypes and genotypes and contributing to a broad phenotypic plasticity. Although HSP90 interaction networks affecting the genetic basis of phenotypic variation and diversity have been thoroughly studied in animals, such studies are just starting to emerge in plants. Here, we summarize current knowledge and discuss HSP90 network functions in plant development and cellular homeostasis.

Analysis of Centranthera grandiflora Benth Transcriptome Explores Genes of Catalpol, Acteoside and Azafrin Biosynthesis

Cardiovascular diseases (CVDs) are a major cause of health loss in the world. Prevention and treatment of this disease by traditional Chinese medicine is a promising method. Centranthera grandiflora Benth is a high-value medicinal herb in the prevention and treatment of CVDs; its main medicinal components include iridoid glycosides, phenylethanoid glycosides, and azafrin in roots. However, biosynthetic pathways of these components and their regulatory mechanisms are unknown. Furthermore, there are no genomic resources of this herb. In this article, we provide sequence and transcript abundance data for the root, stem, and leaf transcriptome of C. grandiflora Benth obtained by the Illumina Hiseq2000. More than 438 million clean reads were obtained from root, stem, and leaf libraries, which produced 153,198 unigenes. Based on databases annotation, a total of 557, 213, and 161 unigenes were annotated to catalpol, acteoside, and azafrin biosynthetic pathways, respectively. Differentially expressed gene analysis identified 14,875 unigenes differentially enriched between leaf and root with 8,054 upregulated genes and 6,821 downregulated genes. Candidate MYB transcription factors involved in catalpol, acteoside, and azafrin biosynthesis were also predicated. This work is the first transcriptome analysis in C. grandiflora Benth which will aid the deciphering of biosynthesis pathways and regulatory mechanisms of active components.

Auxin promotion of seedling growth via ARF5 is dependent on the brassinosteroid-regulated transcription factors BES1 and BEH4

Seedlings must continually calibrate their growth in response to the environment. Auxin and brassinosteroids (BRs) are plant hormones that work together to control growth responses during photomorphogenesis. We used our previous analysis of promoter architecture in an auxin and BR target gene to guide our investigation into the broader molecular bases and biological relevance of transcriptional co-regulation by these hormones. We found that the auxin-regulated transcription factor Auxin Responsive Factor 5 (ARF5) and the brassinosteroid-regulated transcription factor BRI1-EMS-Suppressor 1/Brassinazole Resistant 2 (BES1) co-regulated a subset of growth-promoting genes via conserved bipartite cis-regulatory elements. Moreover, ARF5 binding to DNA could be enriched by increasing BES1 levels. The evolutionary loss of bipartite elements in promoters results in loss of hormone responsiveness. We also identified another member of the BES1/BZR1 family called BEH4 that acts partially redundantly with BES1 to regulate seedling growth. Double mutant analysis showed that BEH4 and not BZR1 were required alongside BES1 for normal auxin response during early seedling development. We propose that an ARF5-BES1/BEH4 transcriptional module acts to promote growth via modulation of a diverse set of growth-associated genes.

Dynamics of Gene Expression in Single Root Cells of Arabidopsis thaliana

Single cell RNA sequencing can yield high-resolution cell-type-specific expression signatures that reveal new cell types and the developmental trajectories of cell lineages. Here, we apply this approach to Arabidopsis (Arabidopsis thaliana) root cells to capture gene expression in 3,121 root cells. We analyze these data with Monocle 3, which orders single cell transcriptomes in an unsupervised manner and uses machine learning to reconstruct single cell developmental trajectories along pseudotime. We identify hundreds of genes with cell-type-specific expression, with pseudotime analysis of several cell lineages revealing both known and novel genes that are expressed along a developmental trajectory. We identify transcription factor motifs that are enriched in early and late cells, together with the corresponding candidate transcription factors that likely drive the observed expression patterns. We assess and interpret changes in total RNA expression along developmental trajectories and show that trajectory branch points mark developmental decisions. Finally, by applying heat stress to whole seedlings, we address the longstanding question of possible heterogeneity among cell types in the response to an abiotic stress. Although the response of canonical heat-shock genes dominates expression across cell types, subtle but significant differences in other genes can be detected among cell types. Taken together, our results demonstrate that single cell transcriptomics holds promise for studying plant development and plant physiology with unprecedented resolution.

Dynamics of Gene Expression in Single Root Cells of Arabidopsis thaliana

Single cell RNA sequencing can yield high-resolution cell-type-specific expression signatures that reveal new cell types and the developmental trajectories of cell lineages. Here, we apply this approach to Arabidopsis (Arabidopsis thaliana) root cells to capture gene expression in 3,121 root cells. We analyze these data with Monocle 3, which orders single cell transcriptomes in an unsupervised manner and uses machine learning to reconstruct single cell developmental trajectories along pseudotime. We identify hundreds of genes with cell-type-specific expression, with pseudotime analysis of several cell lineages revealing both known and novel genes that are expressed along a developmental trajectory. We identify transcription factor motifs that are enriched in early and late cells, together with the corresponding candidate transcription factors that likely drive the observed expression patterns. We assess and interpret changes in total RNA expression along developmental trajectories and show that trajectory branch points mark developmental decisions. Finally, by applying heat stress to whole seedlings, we address the longstanding question of possible heterogeneity among cell types in the response to an abiotic stress. Although the response of canonical heat-shock genes dominates expression across cell types, subtle but significant differences in other genes can be detected among cell types. Taken together, our results demonstrate that single cell transcriptomics holds promise for studying plant development and plant physiology with unprecedented resolution.

On Zebrafish Disease Models and Matters of the Heart

Coronary artery disease (CAD) is the leading form of cardiovascular disease (CVD), which is the primary cause of mortality worldwide. It is a complex disease with genetic and environmental risk factor contributions. Reports in human and mammalian models elucidate age-associated changes in cardiac function. The diverse mechanisms involved in cardiac diseases remain at the center of the research interest to identify novel strategies for prevention and therapy. Zebrafish (Danio rerio) have emerged as a valuable vertebrate model to study cardiovascular development over the last few decades. The facile genetic manipulation via forward and reverse genetic approaches combined with noninvasive, high-resolution imaging and phenotype-based screening has provided new insights to molecular pathways that orchestrate cardiac development. Zebrafish can recapitulate human cardiac pathophysiology due to gene and regulatory pathways conservation, similar heart rate and cardiac morphology and function. Thus, generations of zebrafish models utilize the functional analysis of genes involved in CAD, which are derived from large-scale human population analysis. Here, we highlight recent studies conducted on cardiovascular research focusing on the benefits of the combination of genome-wide association studies (GWAS) with functional genomic analysis in zebrafish. We further summarize the knowledge obtained from zebrafish studies that have demonstrated the architecture of the fundamental mechanisms underlying heart development, homeostasis and regeneration at the cellular and molecular levels.

Redundancy, Feedback, and Robustness in the Arabidopsis thaliana BZR/BEH Gene Family

Organismal development is remarkably robust, tolerating stochastic errors to produce consistent, so-called canalized adult phenotypes. The mechanistic underpinnings of developmental robustness are poorly understood, but recent studies implicate certain features of genetic networks such as functional redundancy, connectivity, and feedback. Here, we examine the BZR/BEH gene family, whose function contributes to embryonic stem development in the plant Arabidopsis thaliana, to test current assumptions on functional redundancy and trait robustness. Our analyses of BZR/BEH gene mutants and mutant combinations revealed that functional redundancy among these gene family members is not necessary for trait robustness. Connectivity is another commonly cited determinant of robustness; however, we found no correlation between connectivity among gene family members or their connectivity with other transcription factors and effects on developmental robustness. Instead, our data suggest that BEH4, the earliest diverged family member, modulates developmental robustness. We present evidence indicating that regulatory cross-talk among gene family members is integrated by BEH4 to promote wild-type levels of developmental robustness. Further, the chaperone HSP90, a known determinant of developmental robustness, appears to act via BEH4 in maintaining robustness of embryonic stem length. In summary, we demonstrate that even among closely related transcription factors, trait robustness can arise through the activity of a single gene family member, challenging common assumptions about the molecular underpinnings of robustness.