Crosstalk between heterotrimeric G protein-coupled signaling pathways and WRKY transcription factors modulating plant responses to suboptimal micronutrient conditions

HOMEOBOX2, the paralog of SIX-ROWED SPIKE1/HOMEOBOX1, is dispensable for barley spikelet development

The HD-ZIP class I transcription factor Homeobox 1 (HvHOX1), also known as Vulgare Row-type Spike 1 (VRS1) or Six-rowed Spike 1, regulates lateral spikelet fertility in barley (Hordeum vulgare L.). It was shown that HvHOX1 has a high expression only in lateral spikelets, while its paralog HvHOX2 was found to be expressed in different plant organs. Yet, the mechanistic functions of HvHOX1 and HvHOX2 during spikelet development are still fragmentary. Here, we show that compared with HvHOX1, HvHOX2 is more highly conserved across different barley genotypes and Hordeum species, hinting at a possibly vital but still unclarified biological role. Using bimolecular fluorescence complementation, DNA-binding, and transactivation assays, we validate that HvHOX1 and HvHOX2 are bona fide transcriptional activators that may potentially heterodimerize. Accordingly, both genes exhibit similar spatiotemporal expression patterns during spike development and growth, albeit their mRNA levels differ quantitatively. We show that HvHOX1 delays the lateral spikelet meristem differentiation and affects fertility by aborting the reproductive organs. Interestingly, the ancestral relationship of the two genes inferred from their co-expressed gene networks suggested that HvHOX1 and HvHOX2 might play a similar role during barley spikelet development. However, CRISPR-derived mutants of HvHOX1 and HvHOX2 demonstrated the suppressive role of HvHOX1 on lateral spikelets, while the loss of HvHOX2 does not influence spikelet development. Collectively, our study shows that through the suppression of reproductive organs, lateral spikelet fertility is regulated by HvHOX1, whereas HvHOX2 is dispensable for spikelet development in barley.

Genome-Wide Comparative Analysis of Five Amaranthaceae Species Reveals a Large Amount of Repeat Content

Amaranthus is a genus of C4 dicotyledonous herbaceous plant species that are widely distributed in Asia, Africa, Australia, and Europe and are used as grain, vegetables, forages, and ornamental plants. Amaranth species have gained significant attention nowadays as potential sources of nutritious food and industrial products. In this study, we performed a comparative genome analysis of five amaranth species, namely, Amaranthus hypochondriacus, Amaranthus tuberculatus, Amaranthus hybridus, Amaranthus palmeri, and Amaranthus cruentus. The estimated repeat content ranged from 54.49% to 63.26% and was not correlated with the genome sizes. Out of the predicted repeat classes, the majority of repetitive sequences were Long Terminal Repeat (LTR) elements, which account for about 13.91% to 24.89% of all amaranth genomes. Phylogenetic analysis based on 406 single-copy orthologous genes revealed that A. hypochondriacus is most closely linked to A. hybridus and distantly related to A. cruentus. However, dioecious amaranth species, such as A. tuberculatus and A. palmeri, which belong to the subgenera Amaranthus Acnida, have formed their distinct clade. The comparative analysis of genomic data of amaranth species will be useful to identify and characterize agronomically important genes and their mechanisms of action. This will facilitate genomics-based, evolutionary studies, and breeding strategies to design faster, more precise, and predictable crop improvement programs.

G protein signaling and metabolic pathways as evolutionarily conserved mechanisms to combat calcium deficiency

Calcium (Ca) deficiency causes necrotic symptoms of foliar edges known as tipburn; however, the underlying cellular mechanisms have been poorly understood due to the lack of an ideal plant model and research platform. Using a phenotyping system that quantitates growth and tipburn traits in the model bryophyte Marchantia polymorpha, we evaluated metabolic compounds and the Gβ-null mutant (gpb1) that modulate the occurrence and expansion of the tipburn. Transcriptomic comparisons between wild-type and gpb1 plants revealed the phenylalanine/phenylpropanoid biosynthesis pathway and reactive oxygen species (ROS) important for Ca deficiency responses. gpb1 plants reduced ROS production possibly through transcriptomic regulations of class III peroxidases and induced the expression of phenylpropanoid pathway enzymes without changing downstream lignin contents. Supplementation of intermediate metabolites and chemical inhibitors further confirmed the cell-protective mechanisms of the phenylpropanoid and ROS pathways. Marchantia polymorpha, Arabidopsis thaliana, and Lactuca sativa showed comparable transcriptomic changes where genes related to phenylpropanoid and ROS pathways were enriched in response to Ca deficiency. In conclusion, our study demonstrated unresolved signaling and metabolic pathways of Ca deficiency response. The phenotyping platform can speed up the discovery of chemical and genetic pathways, which could be widely conserved between M. polymorpha and angiosperms.

The interplay of GTP-binding protein AGB1 with ER stress sensors IRE1a and IRE1b modulates Arabidopsis unfolded protein response and bacterial immunity

In eukaryotic cells, the accumulation of unfolded or misfolded proteins in the endoplasmic reticulum (ER) results in ER stress that induces a cascade of reactions called the unfolded protein response (UPR). In Arabidopsis, the most conserved UPR sensor, Inositol-requiring enzyme 1 (IRE1), responds to both abiotic- and biotic-induced ER stress. Guanine nucleotide-binding proteins (G proteins) constitute another universal and conserved family of signal transducers that have been extensively investigated due to their ubiquitous presence and diverse nature of action. Arabidopsis GTP-binding protein β1 (AGB1) is the only G-protein β-subunit encoded by the Arabidopsis genome that is involved in numerous signaling pathways. Mounting evidence suggests the existence of a crosstalk between IRE1 and G protein signaling during ER stress. AGB1 has previously been shown to control a distinct UPR pathway independently of IRE1 when treated with an ER stress inducer tunicamycin. Our results obtained with combinatorial knockout mutants support the hypothesis that both IRE1 and AGB1 synergistically contribute to ER stress responses chemically induced by dithiothreitol (DTT) as well as to the immune responses against a phytopathogenic bacterium Pseudomonas syringae pv. tomato strain DC3000. Our study highlights the crosstalk between the plant UPR transducers under abiotic and biotic stress.

G protein controls stress readiness by modulating transcriptional and metabolic homeostasis in Arabidopsis thaliana and Marchantia polymorpha

The core G protein signaling module, which consists of Gα and extra-large Gα (XLG) subunits coupled with the Gβγ dimer, is a master regulator of various stress responses. In this study, we compared the basal and salt stress-induced transcriptomic, metabolomic and phenotypic profiles in Gα, Gβ, and XLG-null mutants of two plant species, Arabidopsis thaliana and Marchantia polymorpha, and showed that G protein mediates the shift of transcriptional and metabolic homeostasis to stress readiness status. We demonstrated that such stress readiness serves as an intrinsic protection mechanism against further stressors through enhancing the phenylpropanoid pathway and abscisic acid responses. Furthermore, WRKY transcription factors were identified as key intermediates of G protein-mediated homeostatic shifts. Statistical and mathematical model comparisons between A. thaliana and M. polymorpha revealed evolutionary conservation of transcriptional and metabolic networks over land plant evolution, whereas divergence has occurred in the function of plant-specific atypical XLG subunit. Taken together, our results indicate that the shifts in transcriptional and metabolic homeostasis at least partially act as the mechanisms of G protein-coupled stress responses that are conserved between two distantly related plants.

Diversification of heat shock transcription factors expanded thermal stress responses during early plant evolution

The copy numbers of many plant transcription factor (TF) genes substantially increased during terrestrialization. This allowed TFs to acquire new specificities and thus create gene regulatory networks (GRNs) with new biological functions to help plants adapt to terrestrial environments. Through characterizing heat shock factor (HSF) genes MpHSFA1 and MpHSFB1 in the liverwort Marchantia polymorpha, we explored how heat-responsive GRNs widened their functions in M. polymorpha and Arabidopsis thaliana. An interspecies comparison of heat-induced transcriptomes and the evolutionary rates of HSFs demonstrated the emergence and subsequent rapid evolution of HSFB prior to terrestrialization. Transcriptome and metabolome analyses of M. polymorpha HSF-null mutants revealed that MpHSFA1 controls canonical heat responses such as thermotolerance and metabolic changes. MpHSFB1 also plays essential roles in heat responses, as well as regulating developmental processes including meristem branching and antheridiophore formation. Analysis of cis-regulatory elements revealed development- and stress-related TFs that function directly or indirectly downstream of HSFB. Male gametophytes of M. polymorpha showed higher levels of thermotolerance than female gametophytes, which could be explained by different expression levels of MpHSFA1U and MpHSFA1V on sex chromosome. We propose that the diversification of HSFs is linked to the expansion of HS responses, which enabled coordinated multicellular reactions in land plants.

Interplay between ARABIDOPSIS Gβ and WRKY transcription factors differentiates environmental stress responses

Different environmental stresses often evoke similar physiological disorders such as growth retardation; however, specific consequences reported among individual stresses indicate potential mechanisms to distinguish different stress types in plants. Here, we examined mechanisms to differentiate between stress types in Arabidopsis (Arabidopsis thaliana). Gene expression patterns recapitulating several abiotic stress responses suggested abscisic acid (ABA) as a mediator of the common stress response, while stress type-specific responses were related to metabolic adaptations. Transcriptome and metabolome analyses identified Arabidopsis Gβ (AGB1) mediating the common stress-responsive genes and primary metabolisms under nitrogen excess. AGB1 regulated the expressions of multiple WRKY transcription factors. Gene Ontology and mutant analyses revealed different roles among WRKYs: WRKY40 is involved in ABA and common stress responses, while WRKY75 regulates metabolic processes. The AGB1-WRKY signaling module controlled developmental plasticity in roots under nitrogen excess. Signal transmission from AGB1 to a selective set of WRKYs would be essential to evoke unique responses to different types of stresses.

Metal crossroads in plants: modulation of nutrient acquisition and root development by essential trace metals

The metals iron, zinc, manganese, copper, molybdenum, and nickel are essential for the growth and development of virtually all plant species. Although these elements are required at relatively low amounts, natural factors and anthropogenic activities can significantly affect their availability in soils, inducing deficiencies or toxicities in plants. Because essential trace metals can shape root systems and interfere with the uptake and signaling mechanisms of other nutrients, the non-optimal availability of any of them can induce multi-element changes in plants. Interference by one essential trace metal with the acquisition of another metal or a non-metal nutrient can occur prior to or during root uptake. Essential trace metals can also indirectly impact the plant's ability to capture soil nutrients by targeting distinct root developmental programs and hormone-related processes, consequently inducing largely metal-specific changes in root systems. The presence of metal binding domains in many regulatory proteins also enables essential trace metals to coordinate nutrient uptake by acting at high levels in hierarchical signaling cascades. Here, we summarize the known molecular and cellular mechanisms underlying trace metal-dependent modulation of nutrient acquisition and root development, and highlight the importance of considering multi-element interactions to breed crops better adapted to non-optimal trace metal availabilities.

Solving the puzzle of Fe homeostasis by integrating molecular, mathematical, and societal models

To ensure optimal utilization and bioavailability, iron uptake, transport, subcellular localization, and assimilation are tightly regulated in plants. Herein, we examine recent advances in our understanding of cellular responses to Fe deficiency. We then use intracellular mechanisms of Fe homeostasis to discuss how formalizing cell biology knowledge via a mathematical model can advance discovery even when quantitative data is limited. Using simulation-based inference to identify plausible systems mechanisms that conform to known emergent phenotypes can yield novel, testable hypotheses to guide targeted experiments. However, this approach relies on the accurate encoding of domain-expert knowledge in exploratory mathematical models. We argue that this would be facilitated by fostering more "systems thinking" life scientists and that diversifying your research team may be a practical path to achieve that goal.

Comparative genome-wide analysis of WRKY, MADS-box and MYB transcription factor families in Arabidopsis and rice

Transcription factors (TFs) form the major class of regulatory genes and play key roles in multiple plant stress responses. In most eukaryotic plants, transcription factor (TF) families (WRKY, MADS-box and MYB) activate unique cellular-level abiotic and biotic stress-responsive strategies, which are considered as key determinants for defense and developmental processes. Arabidopsis and rice are two important representative model systems for dicot and monocot plants, respectively. A comprehensive comparative study on 101 OsWRKY, 34 OsMADS box and 122 OsMYB genes (rice genome) and, 71 AtWRKY, 66 AtMADS box and 144 AtMYB genes (Arabidopsis genome) showed various relationships among TFs across species. The phylogenetic analysis clustered WRKY, MADS-box and MYB TF family members into 10, 7 and 14 clades, respectively. All clades in WRKY and MYB TF families and almost half of the total number of clades in the MADS-box TF family are shared between both species. Chromosomal and gene structure analysis showed that the Arabidopsis-rice orthologous TF gene pairs were unevenly localized within their chromosomes whilst the distribution of exon–intron gene structure and motif conservation indicated plausible functional similarity in both species. The abiotic and biotic stress-responsive cis-regulatory element type and distribution patterns in the promoter regions of Arabidopsis and rice WRKY, MADS-box and MYB orthologous gene pairs provide better knowledge on their role as conserved regulators in both species. Co-expression network analysis showed the correlation between WRKY, MADs-box and MYB genes in each independent rice and Arabidopsis network indicating their role in stress responsiveness and developmental processes.

Rice G protein γ subunit qPE9-1 modulates root elongation for phosphorus uptake by involving 14-3-3 protein OsGF14b and plasma membrane H+ -ATPase

Heterotrimeric G protein is involved in plant growth and development, while the role of rice (Oryza sativa) G protein γ subunit qPE9-1 in response to low-phosphorus (LP) conditions remains unclear. The gene expression of qPE9-1 was significantly induced in rice roots under LP conditions. Rice varieties carrying the qPE9-1 allele showed a stronger primary root response to LP than the varieties carrying the qpe9-1 allele (mutant of the qPE9-1 allele). Transgenic rice plants with the qPE9-1 allele had longer primary roots and higher P concentrations than those with the qpe9-1 allele under LP conditions. The plasma membrane (PM) H⁺ -ATPase was important for the qPE9-1-mediated response to LP. Furthermore, OsGF14b, a 14-3-3 protein that acts as a key component in activating PM H⁺ -ATPase for root elongation, is also involved in the qPE9-1 mediation. Moreover, the overexpression of OsGF14b in WYJ8 (carrying the qpe9-1 allele) partially increased primary root length under LP conditions. Experiments using R18 peptide (a 14-3-3 protein inhibitor) showed that qPE9-1 is important for primary root elongation and H⁺ efflux under LP conditions by involving the 14-3-3 protein. In addition, rhizosheath weight, total P content, and the rhizosheath soil Olsen-P concentration of qPE9-1 lines were higher than those of qpe9-1 lines under soil drying and LP conditions. These results suggest that the G protein γ subunit qPE9-1 in rice plants modulates root elongation for phosphorus uptake by involving the 14-3-3 protein OsGF14b and PM H⁺ -ATPase, which is required for rice P use.

Arabidopsis G-Protein β Subunit AGB1 Negatively Regulates DNA Binding of MYB62, a Suppressor in the Gibberellin Pathway

Plant G proteins are versatile components of transmembrane signaling transduction pathways. The deficient mutant of heterotrimeric G protein leads to defects in plant growth and development, suggesting that it regulates the GA pathway in Arabidopsis. However, the molecular mechanism of G protein regulation of the GA pathway is not understood in plants. In this study, two G protein β subunit (AGB1) mutants, agb1-2 and N692967, were dwarfed after exogenous application of GA₃. AGB1 interacts with the DNA-binding domain MYB62, a GA pathway suppressor. Transgenic plants were obtained through overexpression of MYB62 in two backgrounds including the wild-type (MYB62/WT Col-0) and agb1 mutants (MYB62/agb1) in Arabidopsis. Genetic analysis showed that under GA₃ treatment, the height of the transgenic plants MYB62/WT and MYB62/agb1 was lower than that of WT. The height of MYB62/agb1 plants was closer to MYB62/WT plants and higher than that of mutants agb1-2 and N692967, suggesting that MYB62 is downstream of AGB1 in the GA pathway. qRT-PCR and competitive DNA binding assays indicated that MYB62 can bind MYB elements in the promoter of GA2ox7, a GA degradation gene, to activate GA2ox7 transcription. AGB1 affected binding of MYB62 on the promoter of GA2ox7, thereby negatively regulating th eactivity of MYB62.

Two Arabidopsis Homologs of Human Lysine-Specific Demethylase Function in Epigenetic Regulation of Plant Defense Responses

Epigenetic marks such as covalent histone modification and DNA methylation are crucial for mitotically and meiotically inherited cellular memory-based plant immunity. However, the roles of individual players in the epigenetic regulation of plant immunity are not fully understood. Here we reveal the functions of two Arabidopsis thaliana homologs of human lysine-specific demethylase1-like1, LDL1 and LDL2, in the maintenance of methyl groups at lysine 4 of histone H3 and in plant immunity to Pseudomonas syringae infection. The growth of virulent P. syringae strains was reduced in ldl1 and ldl2 single mutants compared to wild-type plants. Local and systemic disease resistance responses, which coincided with the rapid, robust transcription of defense-related genes, were more stably expressed in ldl1 ldl2 double mutants than in the single mutants. At the nucleosome level, mono-methylated histone H3K4 accumulated in ldl1 ldl2 plants genome-wide and in the mainly promoter regions of the defense-related genes examined in this study. Furthermore, in silico comparative analysis of RNA-sequencing and chromatin immunoprecipitation data suggested that several WRKY transcription factors, e.g., WRKY22/40/70, might be partly responsible for the enhanced immunity of ldl1 ldl2. These findings suggest that LDL1 and LDL2 control the transcriptional sensitivity of a group of defense-related genes to establish a primed defense response in Arabidopsis.

Strawberry FaWRKY25 Transcription Factor Negatively Regulated the Resistance of Strawberry Fruits to Botrytis cinerea

WRKY genes and jasmonic acid (JA) play a crucial role in plants’ responses against biotic and abiotic stress. However, the regulating mechanism of WRKY genes on strawberry fruits’ resistance against Botrytis cinerea is largely unknown, and few studies have been performed on their effect on the JA-mediated defense mechanism against B. cinerea. This study explored the effect of FaWRKY25 on the JA-mediated strawberry resistance against B. cinerea. Results showed that the JA content decreased significantly as the fruits matured, whereas the FaWRKY25 expression rose substantially, which led to heightened susceptibility to B. cinerea and in strawberries. External JA treatment significantly increased the JA content in strawberries and reduced the FaWRKY25 expression, thereby enhancing the fruits’ resistance against B. cinerea. FaWRKY25 overexpression significantly lowered the fruits’ resistance against B. cinerea, whereas FaWRKY25 silencing significantly increased resistance. Moreover, FaWRKY25 overexpression significantly lowered the JA content, whereas FaWRKY25 silencing significantly increased it. FaWRKY25 expression level substantially affects the expression levels of genes related to JA biosynthesis and metabolism, other members of the WRKY family, and defense genes. Accordingly, FaWRKY25 plays a crucial role in regulating strawberries’ resistance against B. cinerea and may negatively regulate their JA-mediated resistance mechanism against B. cinerea.