ePlant: Visualizing and Exploring Multiple Levels of Data for Hypothesis Generation in Plant Biology

20 years of the Bio-Analytic Resource for Plant Biology

The Bio-Analytic Resource for Plant Biology ('the BAR', at https://bar.utoronto.ca) is celebrating its 20th year in operation in 2025. The BAR encompasses and provides visualization tools for large 'omics data sets from plants. The BAR covers data from Arabidopsis, tomato, wheat, barley and 29 other plant species (with data for 2 others to be released soon). These data include nucleotide and protein sequence data, gene expression data, protein-protein and protein-DNA interactions, protein structures, subcellular localizations, and polymorphisms. The data are stored in more than 200 relational databases holding 186 GB of data and are presented to the researchers via web apps. These web apps provide data analysis and visualization tools. Some of the most popular tools are eFP ('electronic fluorescent pictograph') Browsers, ePlants and ThaleMine (an Arabidopsis-specific instance of InterMine). The BAR was designated a Global Core Biodata Resource in 2023. Like other GCBRs, the BAR has excellent operational stability, provides access without login requirement, and provides an API for researchers to be able to access BAR data programmatically. We present in this update a new overarching search tool called Gaia that permits easy access to all BAR data, powered by machine learning and artificial intelligence.

Shade-induced ROS/NO reinforce COP1-mediated diffuse cell growth

Canopy shade enhances the activity of PHYTOCHROME INTERACTING FACTORs (PIFs) to boost auxin synthesis in the cotyledons. Auxin, together with local PIFs and their positive regulator CONSTITUTIVELY PHOTOMORPHOGENIC 1 (COP1), promotes hypocotyl growth to facilitate access to light. Whether shade alters the cellular redox status thereby affecting growth responses, remains unexplored. Here, we show that, under shade, high auxin levels increased reactive oxygen species and nitric oxide accumulation in the hypocotyl of Arabidopsis. This nitroxidative environment favored the promotion of hypocotyl growth by COP1 under shade. We demonstrate that COP1 is S-nitrosylated, particularly under shade. Impairing this redox regulation enhanced COP1 degradation by the proteasome and diminished the capacity of COP1 to interact with target proteins and to promote hypocotyl growth. Disabling this regulation also generated transversal asymmetries in hypocotyl growth, indicating poor coordination among different cells, which resulted in random hypocotyl bending and predictably low ability to compete with neighbors. These findings highlight the significance of redox signaling in the control of diffuse growth during shade avoidance.

THI1 Gene Evolutionary Trends: A Comprehensive Plant-Focused Assessment via Data Mining and Large-Scale Analysis

Molecular evolution analysis typically involves identifying selection pressure and reconstructing evolutionary trends. This process usually requires access to specific data related to a target gene or gene family within a particular group of organisms. While recent advancements in high-throughput sequencing techniques have resulted in the rapid accumulation of extensive genomics and transcriptomics data and the creation of new databases in public repositories, extracting valuable insights from such vast data sets remains a significant challenge for researchers. Here, we elucidated the evolutionary history of THI1, a gene responsible for encoding thiamine thiazole synthase. The thiazole ring is a precursor for vitamin B1 and a crucial cofactor in primary metabolic pathways. A thorough search of complete genomes available within public repositories reveals 702 THI1 homologs of Archaea and Eukarya. Throughout its diversification, the plant lineage has preserved the THI1 gene by incorporating the N-terminus and targeting the chloroplasts. Likewise, evolutionary pressures and lifestyle appear to be associated with retention of TPP riboswitch sites and consequent dual posttranscriptional regulation of the de novo biosynthesis pathway in basal groups. Multicopy retention of THI1 is not a typical plant pattern, even after successive genome duplications. Examining cis-regulatory sites in plants uncovers two shared motifs across all plant lineages. A data mining of 484 transcriptome data sets supports the THI1 homolog expression under a light/dark cycle response and a tissue-specific pattern. Finally, the work presented brings a new look at public repositories as an opportunity to explore evolutionary trends to THI1.

Analysis of Ion Transport Properties of Glycine max HKT Transporters and Identifying a Regulation of GmHKT1;1 by the Non-Functional GmHKT1;4

High-affinity potassium transporters (HKTs) play an important role in plants responding to salt stress, but the transport properties of the soybean HKT transporters at the molecular level are still unclear. Here, using Xenopus oocyte as a heterologous expression system and two-electrode voltage-clamp technique, we identified four HKT transporters, GmHKT1;1, GmHKT1;2, GmHKT1;3 and GmHKT1;4, all of which belong to type I subfamily, but have distinct ion transport properties. While GmHKT1;1, GmHKT1;2 and GmHKT1;3 function as Na+ transporters, GmHKT1;1 is less selective against K+ than the two other transporters. Astonishingly, GmHKT1;4, which lacks transmembrane segments and has no ion permeability, is significantly expressed, and its gene expression pattern is different from the other three GmHKTs under salt stress. Interestingly, GmHKT1;4 reduced the Na+/K+ currents mediated by GmHKT1;1. Further study showed that the transport ability of GmHKT1;1 regulated by GmHKT1;4 was related to the structural differences in the first intracellular domain and the fourth repeat domain. Overall, we have identified one unique GmHKT member, GmHKT1;4, which modulates the Na+ and K+ transport ability of GmHKT1;1 via direct interaction. Thus, we have revealed a new type of HKT interaction model for altering their ion transport properties.

Evolutionary analysis of tonoplast intrinsic proteins (TIPs) unraveling the role of TIP3s in plant seed development

Tonoplast intrinsic proteins (TIPs) are crucial in facilitating the transportation of water and various small solutes across biological membranes. The evolutionary path and functional roles of TIPs is poorly understood in plants. In the present study, a total of 976 TIPs were identified in 104 diverse species and subsequently studied to trace their lineage-specific evolutionary path and tissue-specific function. Interestingly, TIPs were found to be absent in lower forms such as algae and fungi and they evolved later in primitive plants like bryophytes. Bryophytes possess a distant class of TIPs, denoted as TIP6, which is not found in higher plants. The aromatic/arginine (ar/R) selectivity filter found in TIP6 of certain liverworts share similarity with hybrid intrinsic protein (HIP), suggesting an evolutionary kinship. As plants evolved to more advanced forms, TIPs diversified into five different sub-groups (TIP1 to TIP5). Notably, TIP5 is a sub-group unique to angiosperms. The evolutionary history of the TIP subfamily reveals an interesting observation that the TIP3 subgroup has evolved within seed-bearing Spermatophyta. Further, TIPs exhibit tissue-specific expression that is conserved within various plant species. Specifically, the TIP3s were found to be exclusively expressed in seeds. Quantitative PCR analysis of TIP3s showed gradually increasing expression in soybean seed developmental stages. The expression of TIP3s in different plant species was also found to be gradually increasing during seed maturation. The results presented here address the knowledge gap concerning the evolutionary background of TIPs, specifically TIP3 in plants, and provide valuable insights for a deeper comprehension of the functions of TIPs in plants.

Molecular simulation reveals that pathogenic mutations in BTB/ANK domains of Arabidopsis thaliana NPR1 circumscribe the EDS1-mediated immune regulation

The NPR1 (nonexpressor of pathogenesis-related genes 1) is a key regulator of the salicylic-acid-mediated immune response caused by pathogens in Arabidopsis thaliana. Mutations C150Y and H334Y in the BTB/ANK domains of NPR1 inhibit the defense response, and transcriptional co-activity with enhanced disease susceptibility 1 (EDS1) has been revealed experimentally. This study examined the conformational changes and reduced NPR1-EDS1 interaction upon mutation using a molecular dynamics simulation. Initially, BTBC150YNPR1 and ANKH334YNPR1 were categorized as pathological mutations rather than others based on sequence conservation. A distant ortholog was used to map the common residues shared among the wild-type because the mutations were highly conserved. Overall, 179 of 373 residues were determining the secondary structures and fold versatility of conformations. In addition, the mutational hotspots Cys150, Asp152, Glu153, Cys155, His157, Cys160, His334, Arg339 and Lys370 were crucial for oligomer-to-monomer exchange. Subsequently, the atomistic simulations with free energy (MM/PB(GB)SA) calculations predicted structural displacements engaging in the N-termini α5133-178α7 linker connecting the central ANK regions (α13260-290α14 and α18320-390α22), where prominent long helices (α516) and short helices (α310) replaced with β-turns and loops disrupting hydrogen bonds and salt bridges in both mutants implicating functional regulation and activation. Furthermore, the mutation repositions the intact stability of multiple regions (L13C149-N356α20BTB/ANK-α17W301-E357α21N-ter/coiled-coil) compromising a dynamic interaction of NPR1-EDS1. By unveiling the transitions between the distinct functions of mutational perception, this study paves the way for future investigation to orchestrate additive host-adapted transcriptional reprogramming that controls defense-related regulatory mechanisms of NPR1s in plants.

ggPlantmap: an open-source R package for the creation of informative and quantitative ggplot maps derived from plant images

As plant research generates an ever-growing volume of spatial quantitative data, the need for decentralized and user-friendly visualization tools to explore large and complex datasets becomes crucial. Existing resources, such as the Plant eFP (electronic Fluorescent Pictograph) viewer, have played a pivotal role on the communication of gene expression data across many plant species. However, although widely used by the plant research community, the Plant eFP viewer lacks open and user-friendly tools for the creation of customized expression maps independently. Plant biologists with less coding experience can often encounter challenges when attempting to explore ways to communicate their own spatial quantitative data. We present 'ggPlantmap' an open-source R package designed to address this challenge by providing an easy and user-friendly method for the creation of ggplot representative maps from plant images. ggPlantmap is built in R, one of the most used languages in biology, to empower plant scientists to create and customize eFP-like viewers tailored to their experimental data. Here, we provide an overview of the package and tutorials that are accessible even to users with minimal R programming experience. We hope that ggPlantmap can assist the plant science community, fostering innovation, and improving our understanding of plant development and function.

Sugar signaling modulates SHOOT MERISTEMLESS expression and meristem function in Arabidopsis

In plants, development of all above-ground tissues relies on the shoot apical meristem (SAM) which balances cell proliferation and differentiation to allow life-long growth. To maximize fitness and survival, meristem activity is adjusted to the prevailing conditions through a poorly understood integration of developmental signals with environmental and nutritional information. Here, we show that sugar signals influence SAM function by altering the protein levels of SHOOT MERISTEMLESS (STM), a key regulator of meristem maintenance. STM is less abundant in inflorescence meristems with lower sugar content, resulting from plants being grown or treated under limiting light conditions. Additionally, sucrose but not light is sufficient to sustain STM accumulation in excised inflorescences. Plants overexpressing the α1-subunit of SUCROSE-NON-FERMENTING1-RELATED KINASE 1 (SnRK1) accumulate less STM protein under optimal light conditions, despite higher sugar accumulation in the meristem. Furthermore, SnRK1α1 interacts physically with STM and inhibits its activity in reporter assays, suggesting that SnRK1 represses STM protein function. Contrasting the absence of growth defects in SnRK1α1 overexpressors, silencing SnRK1α in the SAM leads to meristem dysfunction and severe developmental phenotypes. This is accompanied by reduced STM transcript levels, suggesting indirect effects on STM. Altogether, we demonstrate that sugars promote STM accumulation and that the SnRK1 sugar sensor plays a dual role in the SAM, limiting STM function under unfavorable conditions but being required for overall meristem organization and integrity under favorable conditions. This highlights the importance of sugars and SnRK1 signaling for the proper coordination of meristem activities.

HIGH PLOIDY2-mediated SUMOylation of transcription factor ARR1 controls two-component signaling in Arabidopsis

Cytokinins regulate plant growth, development, and responses to environmental stresses such as cold via phosphorelay from cytokinin receptors to the ARABIDOPSIS RESPONSE REGULATORs (ARRs). However, the molecular mechanisms underlying the activation of type-B ARR transcriptional activity in Arabidopsis (Arabidopsis thaliana) remain unclear. Here, we show that the E3 SUMO ligase HIGH PLOIDY2 SUMOylates ARR1, a type-B ARR, at K236, triggering its activation. Cold- or cytokinin-induced phosphorylation of ARR1 at D89 is crucial for its interaction with HPY2. Lysine 236 is critical for ARR1's transactivation without compromising its DNA-binding ability, while D89 is crucial for ARR1's binding to target gene promoters. Cytokinin enhances ARR1's chromatin binding, but cold does not. ARR1 K236 plays a critical role in promoting histone H3 acetylation in response to both cytokinin and cold without affecting chromatin binding. The K236R mutation in ARR1 reduces target gene expression and alters cytokinin and cold response phenotypes. This study unveils a mechanism of ARR1 activation wherein phosphorylated ARR1 interacts with HPY2 and binds to chromatin in response to cytokinin. Cold triggers a phosphorelay targeting chromatin-bound ARR1. HPY2 then catalyzes ARR1 SUMOylation at K236, enhancing histone H3 acetylation and leading to transcriptional activation of ARR1 in response to both cold and cytokinin.

Characterization of Arabidopsis aldolases AtFBA4, AtFBA5, and their inhibition by morin and interaction with calmodulin

Fructose bisphosphate aldolases (FBAs) catalyze the reversible cleavage of fructose 1,6-bisphosphate into dihydroxyacetone phosphate and glyceraldehyde 3-phosphate. We analyzed two previously uncharacterized cytosolic Arabidopsis FBAs, AtFBA4 and AtFBA5. Based on a recent report, we examined the interaction of AtFBA4 with calmodulin (CaM)-like protein 11 (AtCML11). AtFBA4 did not bind AtCML11; however, we found that CaM bound AtFBA5 in a Ca2+-dependent manner with high specificity and affinity (KD ~ 190 nm) and enhanced its stability. AtFBA4 and AtFBA5 exhibited Michaelis-Menten kinetics with Km and Vmax values of 180 μm and 4.9 U·mg-1 for AtFBA4, and 6.0 μm and 0.30 U·mg-1 for AtFBA5, respectively. The flavonoid morin inhibited both isozymes. Our study suggests that Ca2+ signaling and flavanols may influence plant glycolysis/gluconeogenesis.

Characterization of the small Arabidopsis thaliana GTPase and ADP-ribosylation factor-like 2 protein TITAN 5

Small GTPases switch between GDP- and GTP-bound states during cell signaling. The ADP-ribosylation factor (ARF) family of small GTPases is involved in vesicle trafficking. Although evolutionarily well conserved, little is known about ARF and ARF-like GTPases in plants. We characterized biochemical properties and cellular localization of the essential small ARF-like GTPase TITAN 5 (TTN5; also known as HALLIMASCH, ARL2 and ARLC1) from Arabidopsis thaliana, and two TTN5 proteins with point mutants in conserved residues, TTN5T30N and TTN5Q70L, that were expected to be unable to perform nucleotide exchange and GTP hydrolysis, respectively. TTN5 exhibited very rapid intrinsic nucleotide exchange and remarkably low GTP hydrolysis activity, functioning as a non-classical small GTPase being likely present in a GTP-loaded active form. We analyzed signals from YFP-TTN5 and HA3-TTN5 by in situ immunolocalization in Arabidopsis seedlings and through use of a transient expression system. Colocalization with endomembrane markers and pharmacological treatments suggests that TTN5 can be present at the plasma membrane and that it dynamically associates with membranes of vesicles, Golgi stacks and multivesicular bodies. Although TTN5Q70L mirrored wild-type TTN5 behavior, the TTN5T30N mutant differed in some aspects. Hence, the unusual rapid nucleotide exchange activity of TTN5 is linked with its membrane dynamics, and TTN5 likely has a role in vesicle transport within the endomembrane system.

Sugar sensing in C4 source leaves: a gap that needs to be filled

Plant growth depends on sugar production and export by photosynthesizing source leaves and sugar allocation and import by sink tissues (grains, roots, stems, and young leaves). Photosynthesis and sink demand are tightly coordinated through metabolic (substrate, allosteric) feedback and signalling (sugar, hormones) mechanisms. Sugar signalling integrates sugar production with plant development and environmental cues. In C3 plants (e.g. wheat and rice), it is well documented that sugar accumulation in source leaves, due to source-sink imbalance, negatively feeds back on photosynthesis and plant productivity. However, we have a limited understanding about the molecular mechanisms underlying those feedback regulations, especially in C4 plants (e.g. maize, sorghum, and sugarcane). Recent work with the C4 model plant Setaria viridis suggested that C4 leaves have different sugar sensing thresholds and behaviours relative to C3 counterparts. Addressing this research priority is critical because improving crop yield requires a better understanding of how plants coordinate source activity with sink demand. Here we review the literature, present a model of action for sugar sensing in C4 source leaves, and suggest ways forward.

Application of orthology and network biology to infer gene functions in non-model plants

Approximately 60% of the genes and gene products in the model species Arabidopsis thaliana have been functionally characterized. In non-model plant species, the functional annotation of the gene space is largely based on homology, with the assumption that genes with shared common ancestry have conserved functions. However, the wide variety in possible morphological, physiological, and ecological differences between plant species gives rise to many species- and clade-specific genes, for which this transfer of knowledge is not possible. Other complications, such as difficulties with genetic transformation, the absence of large-scale mutagenesis methods, and long generation times, further lead to the slow characterization of genes in non-model species. Here, we discuss different resources that integrate plant gene function information. Different approaches that support the functional annotation of gene products, based on orthology or network biology, are described. While sequence-based tools to characterize the functional landscape in non-model species are maturing and becoming more readily available, easy-to-use network-based methods inferring plant gene functions are not as prevalent and have limited functionality.

P-type calcium ATPases play important roles in biotic and abiotic stress signaling

MAIN CONCLUSION

Knowledge of Ca2+-ATPases is imperative for improving crop quality/ food security, highly threatened due to global warming. Ca2+-ATPases modulates calcium, essential for stress signaling and modulating growth, development, and immune activities. Calcium is considered a versatile secondary messenger and essential for short- and long-term responses to biotic and abiotic stresses in plants. Coordinated transport activities from both calcium influx and efflux channels are required to generate cellular calcium signals. Various extracellular stimuli cause an induction in cytosolic calcium levels. To cope with such stresses, it is important to maintain intracellular Ca2+ levels. Plants need to evolve efficient efflux mechanisms to maintain Ca2+ ion homeostasis. Plant Ca2+-ATPases are members of the P-type ATPase superfamily and localized in the plasma membrane and endoplasmic reticulum (ER). They are required for various cellular processes, including plant growth, development, calcium signaling, and even retorts to environmental stress. These ATPases play an essential role in Ca2+ homeostasis and are actively involved in Ca2+ transport. Plant Ca2+-ATPases are categorized into two major classes: type IIA and type IIB. Although these two classes of ATPases share similarities in protein sequence, they differ in their structure, cellular localization, and sensitivity to inhibitors. Due to the emerging role of Ca2+-ATPase in abiotic and biotic plant stress, members of this family may help promote agricultural improvement under stress conditions. This review provides a comprehensive overview of P-type Ca2+-ATPase, and their role in Ca2+ transport, stress signaling, and cellular homeostasis focusing on their classification, evolution, ion specificities, and catalytic mechanisms. It also describes the main aspects of the role of Ca2+-ATPase in transducing signals during plant biotic and abiotic stress responses and its role in plant development and physiology.

Arabidopsis CaLB1 undergoes phase separation with the ESCRT protein ALIX and modulates autophagosome maturation

Autophagy is relevant for diverse processes in eukaryotic cells, making its regulation of fundamental importance. The formation and maturation of autophagosomes require a complex choreography of numerous factors. The endosomal sorting complex required for transport (ESCRT) is implicated in the final step of autophagosomal maturation by sealing of the phagophore membrane. ESCRT-III components were shown to mediate membrane scission by forming filaments that interact with cellular membranes. However, the molecular mechanisms underlying the recruitment of ESCRTs to non-endosomal membranes remain largely unknown. Here we focus on the ESCRT-associated protein ALG2-interacting protein X (ALIX) and identify Ca2+-dependent lipid binding protein 1 (CaLB1) as its interactor. Our findings demonstrate that CaLB1 interacts with AUTOPHAGY8 (ATG8) and PI(3)P, a phospholipid found in autophagosomal membranes. Moreover, CaLB1 and ALIX localize with ATG8 on autophagosomes upon salt treatment and assemble together into condensates. The depletion of CaLB1 impacts the maturation of salt-induced autophagosomes and leads to reduced delivery of autophagosomes to the vacuole. Here, we propose a crucial role of CaLB1 in augmenting phase separation of ALIX, facilitating the recruitment of ESCRT-III to the site of phagophore closure thereby ensuring efficient maturation of autophagosomes.

Localization of four class I glutaredoxins in the cytosol and the secretory pathway and characterization of their biochemical diversification

Class I glutaredoxins (GRXs) are catalytically active oxidoreductases and considered key proteins mediating reversible glutathionylation and deglutathionylation of protein thiols during development and stress responses. To narrow in on putative target proteins, it is mandatory to know the subcellular localization of the respective GRXs and to understand their catalytic activities and putative redundancy between isoforms in the same compartment. We show that in Arabidopsis thaliana, GRXC1 and GRXC2 are cytosolic proteins with GRXC1 being attached to membranes through myristoylation. GRXC3 and GRXC4 are identified as type II membrane proteins along the early secretory pathway with their enzymatic function on the luminal side. Unexpectedly, neither single nor double mutants lacking both GRXs isoforms in the cytosol or the ER show phenotypes that differ from wild-type controls. Analysis of electrostatic surface potentials and clustering of GRXs based on their electrostatic interaction with roGFP2 mirrors the phylogenetic classification of class I GRXs, which clearly separates the cytosolic GRXC1 and GRXC2 from the luminal GRXC3 and GRXC4. Comparison of all four studied GRXs for their oxidoreductase function highlights biochemical diversification with GRXC3 and GRXC4 being better catalysts than GRXC1 and GRXC2 for the reduction of bis(2-hydroxyethyl) disulfide. With oxidized roGFP2 as an alternative substrate, GRXC1 and GRXC2 catalyze the reduction faster than GRXC3 and GRXC4, which suggests that catalytic efficiency of GRXs in reductive reactions depends on the respective substrate. Vice versa, GRXC3 and GRXC4 are faster than GRXC1 and GRXC2 in catalyzing the oxidation of pre-reduced roGFP2 in the reverse reaction.

Class I TCP in fruit development: much more than growth

Fruit development can be viewed as the succession of three main steps consisting of the fruit initiation, growth and ripening. These processes are orchestrated by different factors, notably the successful fertilization of flowers, the environmental conditions and the hormones whose action is coordinated by a large variety of transcription factors. Among the different transcription factor families, TEOSINTE BRANCHED 1, CYCLOIDEA, PROLIFERATING CELL FACTOR (TCP) family has received little attention in the frame of fruit biology despite its large effects on several developmental processes and its action as modulator of different hormonal pathways. In this respect, the comprehension of TCP functions in fruit development remains an incomplete puzzle that needs to be assembled. Building on the abundance of genomic and transcriptomic data, this review aims at collecting available TCP expression data to allow their integration in the light of the different functional genetic studies reported so far. This reveals that several Class I TCP genes, already known for their involvement in the cell proliferation and growth, display significant expression levels in developing fruit, although clear evidence supporting their functional significance in this process remains scarce. The extensive expression data compiled in our study provide convincing elements that shed light on the specific involvement of Class I TCP genes in fruit ripening, once these reproductive organs acquire their mature size. They also emphasize their putative role in the control of specific biological processes such as fruit metabolism and hormonal dialogue.

Characterization of the small Arabidopsis thaliana GTPase and ADP-ribosylation factor-like 2 protein TITAN 5

Small GTPases function by conformational switching ability between GDP- and GTP-bound states in rapid cell signaling events. The ADP-ribosylation factor (ARF) family is involved in vesicle trafficking. Though evolutionarily well conserved, little is known about ARF and ARF-like GTPases in plants. Here, we characterized biochemical properties and cellular localization of the essential small ARF-like GTPase TITAN 5/HALLIMASCH/ARL2/ARLC1 (hereafter termed TTN5) from Arabidopsis thaliana. Two TTN5 variants were included in the study with point mutations at conserved residues, suspected to be functional for nucleotide exchange and GTP hydrolysis, TTN5T30N and TTN5Q70L. We found that TTN5 had a very rapid intrinsic nucleotide exchange capacity with a conserved nucleotide switching mechanism. TTN5 acted as a non-classical small GTPase with a remarkably low GTP hydrolysis activity, suggesting it is likely present in GTP-loaded active form in the cell. We analyzed signals from yellow fluorescent protein (YFP)-tagged TTN5 and from in situ immunolocalization of hemagglutine-tagged HA3-TTN5 in Arabidopsis seedlings and in a transient expression system. Together with colocalization using endomembrane markers and pharmacological treatments the microscopic analysis suggests that TTN5 can be present at the plasma membrane and dynamically associated with membranes of vesicles, Golgi stacks and multivesicular bodies. While the TTN5Q70L variant showed similar GTPase activities and localization behavior as wild-type TTN5, the TTN5T30N mutant differed in some aspects. Hence, the unusual capacity of rapid nucleotide exchange activity of TTN5 is linked with cell membrane dynamics, likely associated with vesicle transport pathways in the endomembrane system.

A method for mining condition-specific co-expressed genes in Camellia sinensis based on k-means clustering

BACKGROUND

As one of the world's most important beverage crops, tea plants (Camellia sinensis) are renowned for their unique flavors and numerous beneficial secondary metabolites, attracting researchers to investigate the formation of tea quality. With the increasing availability of transcriptome data on tea plants in public databases, conducting large-scale co-expression analyses has become feasible to meet the demand for functional characterization of tea plant genes. However, as the multidimensional noise increases, larger-scale co-expression analyses are not always effective. Analyzing a subset of samples generated by effectively downsampling and reorganizing the global sample set often leads to more accurate results in co-expression analysis. Meanwhile, global-based co-expression analyses are more likely to overlook condition-specific gene interactions, which may be more important and worthy of exploration and research.

RESULTS

Here, we employed the k-means clustering method to organize and classify the global samples of tea plants, resulting in clustered samples. Metadata annotations were then performed on these clustered samples to determine the "conditions" represented by each cluster. Subsequently, we conducted gene co-expression network analysis (WGCNA) separately on the global samples and the clustered samples, resulting in global modules and cluster-specific modules. Comparative analyses of global modules and cluster-specific modules have demonstrated that cluster-specific modules exhibit higher accuracy in co-expression analysis. To measure the degree of condition specificity of genes within condition-specific clusters, we introduced the correlation difference value (CDV). By incorporating the CDV into co-expression analyses, we can assess the condition specificity of genes. This approach proved instrumental in identifying a series of high CDV transcription factor encoding genes upregulated during sustained cold treatment in Camellia sinensis leaves and buds, and pinpointing a pair of genes that participate in the antioxidant defense system of tea plants under sustained cold stress.

CONCLUSIONS

To summarize, downsampling and reorganizing the sample set improved the accuracy of co-expression analysis. Cluster-specific modules were more accurate in capturing condition-specific gene interactions. The introduction of CDV allowed for the assessment of condition specificity in gene co-expression analyses. Using this approach, we identified a series of high CDV transcription factor encoding genes related to sustained cold stress in Camellia sinensis. This study highlights the importance of considering condition specificity in co-expression analysis and provides insights into the regulation of the cold stress in Camellia sinensis.

Nuclear lamina component KAKU4 regulates chromatin states and transcriptional regulation in the Arabidopsis genome

BACKGROUND

The nuclear lamina links the nuclear membrane to chromosomes and plays a crucial role in regulating chromatin states and gene expression. However, current knowledge of nuclear lamina in plants is limited compared to animals and humans.

RESULTS

This study mainly focused on elucidating the mechanism through which the putative nuclear lamina component protein KAKU4 regulates chromatin states and gene expression in Arabidopsis leaves. Thus, we constructed a network using the association proteins of lamin-like proteins, revealing that KAKU4 is strongly associated with chromatin or epigenetic modifiers. Then, we conducted ChIP-seq technology to generate global epigenomic profiles of H3K4me3, H3K27me3, and H3K9me2 in Arabidopsis leaves for mutant (kaku4-2) and wild-type (WT) plants alongside RNA-seq method to generate gene expression profiles. The comprehensive chromatin state-based analyses indicate that the knockdown of KAKU4 has the strongest effect on H3K27me3, followed by H3K9me2, and the least impact on H3K4me3, leading to significant changes in chromatin states in the Arabidopsis genome. We discovered that the knockdown of the KAKU4 gene caused a transition between two types of repressive epigenetics marks, H3K9me2 and H3K27me3, in some specific PLAD regions. The combination analyses of epigenomic and transcriptomic data between the kaku4-2 mutant and WT suggested that KAKU4 may regulate key biological processes, such as programmed cell death and hormone signaling pathways, by affecting H3K27me3 modification in Arabidopsis leaves.

CONCLUSIONS

In summary, our results indicated that KAKU4 is directly and/or indirectly associated with chromatin/epigenetic modifiers and demonstrated the essential roles of KAKU4 in regulating chromatin states, transcriptional regulation, and diverse biological processes in Arabidopsis.

NKS1/ELMO4 is an integral protein of a pectin synthesis protein complex and maintains Golgi morphology and cell adhesion in Arabidopsis

Adjacent plant cells are connected by specialized cell wall regions, called middle lamellae, which influence critical agricultural characteristics, including fruit ripening and organ abscission. Middle lamellae are enriched in pectin polysaccharides, specifically homogalacturonan (HG). Here, we identify a plant-specific Arabidopsis DUF1068 protein, called NKS1/ELMO4, that is required for middle lamellae integrity and cell adhesion. NKS1 localizes to the Golgi apparatus and loss of NKS1 results in changes to Golgi structure and function. The nks1 mutants also display HG deficient phenotypes, including reduced seedling growth, changes to cell wall composition, and tissue integrity defects. These phenotypes are comparable to qua1 and qua2 mutants, which are defective in HG biosynthesis. Notably, genetic interactions indicate that NKS1 and the QUAs work in a common pathway. Protein interaction analyses and modeling corroborate that they work together in a stable protein complex with other pectin-related proteins. We propose that NKS1 is an integral part of a large pectin synthesis protein complex and that proper function of this complex is important to support Golgi structure and function.

Profiling of the gene expression and alternative splicing landscapes of Eucalyptus grandis

Eucalyptus is a widely planted hardwood tree species due to its fast growth, superior wood properties and adaptability. However, the post-transcriptional regulatory mechanisms controlling tissue development and stress responses in Eucalyptus remain poorly understood. In this study, we performed a comprehensive analysis of the gene expression profile and the alternative splicing (AS) landscape of E. grandis using strand-specific RNA-Seq, which encompassed 201 libraries including different organs, developmental stages, and environmental stresses. We identified 10 416 genes (33.49%) that underwent AS, and numerous differentially expressed and/or differential AS genes involved in critical biological processes, such as primary-to-secondary growth transition of stems, adventitious root formation, aging and responses to phosphorus- or boron-deficiency. Co-expression analysis of AS events and gene expression patterns highlighted the potential upstream regulatory role of AS events in multiple processes. Additionally, we highlighted the lignin biosynthetic pathway to showcase the potential regulatory functions of AS events in the KNAT3 and IRL3 genes within this pathway. Our high-quality expression atlas and AS landscape serve as valuable resources for unravelling the genetic control of woody plant development, long-term adaptation, and understanding transcriptional diversity in Eucalyptus. Researchers can conveniently access these resources through the interactive ePlant browser (https://bar.utoronto.ca/eplant_eucalyptus).

Rice requires a chromatin remodeler for Polymerase IV-small interfering RNA production and genomic immunity

Transgenes are often spontaneously silenced, which hinders the application of genetic modifications to crop breeding. While gene silencing has been extensively studied in Arabidopsis (Arabidopsis thaliana), the molecular mechanism of transgene silencing remains elusive in crop plants. We used rice (Oryza sativa) plants silenced for a 35S::OsGA2ox1 (Gibberellin 2-oxidase 1) transgene to isolate five elements mountain (fem) mutants showing restoration of transgene expression. In this study, we isolated multiple fem2 mutants defective in a homolog of Required to Maintain Repression 1 (RMR1) of maize (Zea mays) and CLASSY (CLSY) of Arabidopsis. In addition to failing to maintain transgene silencing, as occurs in fem3, in which mutation occurs in NUCLEAR RNA POLYMERASE E1 (OsNRPE1), the fem2 mutant failed to establish transgene silencing of 35S::OsGA2ox1. Mutation in FEM2 eliminated all RNA POLYMERASE IV (Pol-IV)-FEM1/OsRDR2 (RNA-DEPENDENT RNA POLYMERASE 2)-dependent small interfering RNAs (siRNAs), reduced DNA methylation on genome-wide scale in rice seedlings, caused pleiotropic developmental defects, and increased disease resistance. Simultaneous mutation in 2 FEM2 homologous genes, FEM2-Like 1 (FEL1) and FEL2, however, did not affect DNA methylation and rice development and disease resistance. The predominant expression of FEM2 over FEL1 and FEL2 in various tissues was likely caused by epigenetic states. Overexpression of FEL1 but not FEL2 partially rescued hypomethylation of fem2, indicating that FEL1 maintains the cryptic function. In summary, FEM2 is essential for establishing and maintaining gene silencing; moreover, FEM2 is solely required for Pol IV-FEM1 siRNA biosynthesis and de novo DNA methylation.

Rice OseIF6.1 encodes a eukaryotic translation initiation factor and is essential for the development of grain and anther

The eIF6 proteins are distributed extensively in eukaryotes and play diverse and essential roles. The bona fide eIF6 protein in Arabidopsis, At-eIF6;1, is essential for embryogenesis. However, the role of eIF6 proteins in rice growth and development remains elusive and requires further investigation. Here, we characterized the functions of OseIF6.1, which is homologous to At-eIF6;1. OseIF6.1 encodes an eukaryotic translation initiation factor with a conserved eIF6 domain. The knockdown of OseIF6.1 resulted in a decrease in grain length and pollen sterility, whereas the overexpression of OseIF6.1 displayed opposite phenotypes. Further studies revealed that OseIF6.1 regulates grain shape by influencing cell expansion and proliferation. In addition, OseIF6.1 interacts with OsNMD3, which is a nuclear export adaptor for the 60S ribosomal subunit. The knockdown of OsNMD3 in plants exhibited reduced fertility and seed setting. Therefore, our findings have significantly enriched the current understanding of the role of OseIF6.1 in rice growth and development.

Chalkiness and premature controlled by energy homeostasis in OsNAC02 Ko-mutant during vegetative endosperm development

BACKGROUND

Chalkiness is a common phenotype induced by various reasons, such as abiotic stress or the imbalance of starch synthesis and metabolism during the development period. However, the reason mainly for one gene losing its function such as NAC (TFs has a large family in rice) which may cause premature is rarely known to us.

RESULTS

The Ko-Osnac02 mutant demonstrated an obviously early maturation stage compared to the wild type (WT) with 15 days earlier. The result showed that the mature endosperm of Ko-Osnac02 mutant exhibited chalkiness, characterized by white-core and white-belly in mature endosperm. As grain filling rate is a crucial factor in determining the yield and quality of rice (Oryza sativa, ssp. japonica), it's significant that mutant has a lower amylose content (AC) and higher soluble sugar content in the mature endosperm. Interestingly among the top DEGs in the RNA sequencing of N2 (3DAP) and WT seeds revealed that the OsBAM2 (LOC_Os10g32810) expressed significantly high in N2 mutant, which involved in Maltose up-regulated by the starch degradation. As Prediction of Protein interaction showed in the chalky endosperm formation in N2 seeds (3 DAP), seven genes were expressed at a lower-level which should be verified by a heatmap diagrams based on DEGs of N2 versus WT. The Tubulin genes controlling cell cycle are downregulated together with the MCM family genes MCM4 ( ↓), MCM7 ( ↑), which may cause white-core in the early endosperm development. In conclusion, the developing period drastically decreased in the Ko-Osnac02 mutants, which might cause the chalkiness in seeds during the early endosperm development.

CONCLUSIONS

The gene OsNAC02 which controls a great genetic co-network for cell cycle regulation in early development, and KO-Osnac02 mutant shows prematurity and white-core in endosperm.

KAKU4 regulates leaf senescence through modulation of H3K27me3 deposition in the Arabidopsis genome

Lamins are the major components of the nuclear lamina, which regulate chromatin structure and gene expression. KAKU4 is a unique nuclear lamina component in the nuclear periphery, modulates nuclear shape and size in Arabidopsis. The knowledge about the regulatory role of KAKU4 in leaf development remains limited. Here we found that knockdown of KAKU4 resulted in an accelerated leaf senescence phenotype, with elevated levels of H2O2 and hormones, particularly SA, JA, and ABA. Our results demonstrated the importance of KAKU4 as a potential negative regulator in age-triggered leaf senescence in Arabidopsis. Furthermore, we conducted combination analyses of transcriptomic and epigenomic data for the kaku4 mutant and WT leaves. The knockdown of KAKU4 lowered H3K27me3 deposition in the up-regulated genes associated with hormone pathways, programmed cell death, and leaf senescence, including SARD1, SAG113/HAI1, PR2, and so forth. In addition, we found the functional crosstalks between KAKU4 and its associated proteins (CRWN1/4, PNET2, GBPL3, etc.) through comparing multiple transcriptome datasets. Overall, our results indicated that KAKU4 may inhibit the expression of a series of genes related to hormone signals and H2O2 metabolism by affecting the deposition of H3K27me3, thereby suppressing leaf senescence.

spatialHeatmap: visualizing spatial bulk and single-cell assays in anatomical images

Visualizing spatial assay data in anatomical images is vital for understanding biological processes in cell, tissue, and organ organizations. Technologies requiring this functionality include traditional one-at-a-time assays, and bulk and single-cell omics experiments, including RNA-seq and proteomics. The spatialHeatmap software provides a series of powerful new methods for these needs, and allows users to work with adequately formatted anatomical images from public collections or custom images. It colors the spatial features (e.g. tissues) annotated in the images according to the measured or predicted abundance levels of biomolecules (e.g. mRNAs) using a color key. This core functionality of the package is called a spatial heatmap plot. Single-cell data can be co-visualized in composite plots that combine spatial heatmaps with embedding plots of high-dimensional data. The resulting spatial context information is essential for gaining insights into the tissue-level organization of single-cell data, or vice versa. Additional core functionalities include the automated identification of biomolecules with spatially selective abundance patterns and clusters of biomolecules sharing similar abundance profiles. To appeal to both non-expert and computational users, spatialHeatmap provides a graphical and a command-line interface, respectively. It is distributed as a free, open-source Bioconductor package (https://bioconductor.org/packages/spatialHeatmap) that users can install on personal computers, shared servers, or cloud systems.

Single-cell transcriptomics is revolutionizing the improvement of plant biotechnology research: recent advances and future opportunities

Single-cell approaches are a promising way to obtain high-resolution transcriptomics data and have the potential to revolutionize the study of plant growth and development. Recent years have seen the advent of unprecedented technological advances in the field of plant biology to study the transcriptional information of individual cells by single-cell RNA sequencing (scRNA-seq). This review focuses on the modern advancements of single-cell transcriptomics in plants over the past few years. In addition, it also offers a new insight of how these emerging methods will expedite advance research in plant biotechnology in the near future. Lastly, the various technological hurdles and inherent limitations of single-cell technology that need to be conquered to develop such outstanding possible knowledge gain is critically analyzed and discussed.

Plant biotechnology research with single-cell transcriptome: recent advancements and prospects

KEY MESSAGE

Single-cell transcriptomic techniques have emerged as powerful tools in plant biology, offering high-resolution insights into gene expression at the individual cell level. This review highlights the rapid expansion of single-cell technologies in plants, their potential in understanding plant development, and their role in advancing plant biotechnology research. Single-cell techniques have emerged as powerful tools to enhance our understanding of biological systems, providing high-resolution transcriptomic analysis at the single-cell level. In plant biology, the adoption of single-cell transcriptomics has seen rapid expansion of available technologies and applications. This review article focuses on the latest advancements in the field of single-cell transcriptomic in plants and discusses the potential role of these approaches in plant development and expediting plant biotechnology research in the near future. Furthermore, inherent challenges and limitations of single-cell technology are critically examined to overcome them and enhance our knowledge and understanding.

The root-knot nematode effector Mi2G02 hijacks a host plant trihelix transcription factor to promote nematode parasitism

Root-knot nematodes (RKNs) cause huge agricultural losses every year. They secrete a repertoire of effectors to facilitate parasitism through the induction of plant-derived giant feeding cells, which serve as their sole source of nutrients. However, the mode of action of these effectors and their targeted host proteins remain largely unknown. In this study, we investigated the role of the effector Mi2G02 in Meloidogyne incognita parasitism. Host-derived Mi2G02 RNA interference in Arabidopsis thaliana affected giant cell development, whereas ectopic expression of Mi2G02 promoted root growth and increased plant susceptibility to M. incognita. We used various combinations of approaches to study the specific interactions between Mi2G02 and A. thaliana GT-3a, a trihelix transcription factor. GT-3a knockout in A. thaliana affected feeding-site development, resulting in production of fewer egg masses, whereas GT-3a overexpression in A. thaliana increased susceptibility to M. incognita and also root growth. Moreover, we demonstrated that Mi2G02 plays a role in maintaining GT-3a protein stabilization by inhibiting the 26S proteasome-dependent pathway, leading to suppression of TOZ and RAD23C expression and thus promoting nematode parasitism. This work enhances our understanding of how a pathogen effector manipulates the role and regulation of a transcription factor by interfering with a proteolysis pathway to reprogram gene expression for development of nematode feeding cells.

Phototropin2 3'UTR overlaps with the AT5G58150 gene encoding an inactive RLK kinase

BACKGROUND

This study examines the biological implications of an overlap between two sequences in the Arabidopsis genome, the 3'UTR of the PHOT2 gene and a putative AT5G58150 gene, encoded on the complementary strand. AT5G58150 is a probably inactive protein kinase that belongs to the transmembrane, leucine-rich repeat receptor-like kinase family. Phot2 is a membrane-bound UV/blue light photoreceptor kinase. Thus, both proteins share their cellular localization, on top of the proximity of their loci.

RESULTS

The extent of the overlap between 3'UTR regions of AT5G58150 and PHOT2 was found to be 66 bp, using RACE PCR. Both the at5g58150 T-DNA SALK_093781C (with insertion in the promoter region) and 35S::AT5G58150-GFP lines overexpress the AT5G58150 gene. A detailed analysis did not reveal any substantial impact of PHOT2 or AT5G58150 on their mutual expression levels in different light and osmotic stress conditions. AT5G58150 is a plasma membrane protein, with no apparent kinase activity, as tested on several potential substrates. It appears not to form homodimers and it does not interact with PHOT2. Lines that overexpress AT5G58150 exhibit a greater reduction in lateral root density due to salt and osmotic stress than wild-type plants, which suggests that AT5G58150 may participate in root elongation and formation of lateral roots. In line with this, mass spectrometry analysis identified proteins with ATPase activity, which are involved in proton transport and cell elongation, as putative interactors of AT5G58150. Membrane kinases, including other members of the LRR RLK family and BSK kinases (positive regulators of brassinosteroid signalling), can also act as partners for AT5G58150.

CONCLUSIONS

AT5G58150 is a membrane protein that does not exhibit measurable kinase activity, but is involved in signalling through interactions with other proteins. Based on the interactome and root architecture analysis, AT5G58150 may be involved in plant response to salt and osmotic stress and the formation of roots in Arabidopsis.

Mutation of BAM2 rescues the sunn hypernodulation phenotype in Medicago truncatula, suggesting that a signaling pathway like CLV1/BAM in Arabidopsis affects nodule number

The unique evolutionary adaptation of legumes for nitrogen-fixing symbiosis leading to nodulation is tightly regulated by the host plant. The autoregulation of nodulation (AON) pathway negatively regulates the number of nodules formed in response to the carbon/nitrogen metabolic status of the shoot and root by long-distance signaling to and from the shoot and root. Central to AON signaling in the shoots of Medicago truncatula is SUNN, a leucine-rich repeat receptor-like kinase with high sequence similarity with CLAVATA1 (CLV1), part of a class of receptors in Arabidopsis involved in regulating stem cell populations in the root and shoot. This class of receptors in Arabidopsis includes the BARELY ANY MERISTEM family, which, like CLV1, binds to CLE peptides and interacts with CLV1 to regulate meristem development. M. truncatula contains five members of the BAM family, but only MtBAM1 and MtBAM2 are highly expressed in the nodules 48 hours after inoculation. Plants carry mutations in individual MtBAMs, and several double BAM mutant combinations all displayed wild-type nodule number phenotypes. However, Mtbam2 suppressed the sunn-5 hypernodulation phenotype and partially rescued the short root length phenotype of sunn-5 when present in a sunn-5 background. Grafting determined that bam2 suppresses supernodulation from the roots, regardless of the SUNN status of the root. Overexpression of MtBAM2 in wild-type plants increases nodule numbers, while overexpression of MtBAM2 in some sunn mutants rescues the hypernodulation phenotype, but not the hypernodulation phenotypes of AON mutant rdn1-2 or crn. Relative expression measurements of the nodule transcription factor MtWOX5 downstream of the putative bam2 sunn-5 complex revealed disruption of meristem signaling; while both bam2 and bam2 sunn-5 influence MtWOX5 expression, the expression changes are in different directions. We propose a genetic model wherein the specific root interactions of BAM2/SUNN are critical for signaling in nodule meristem cell homeostasis in M. truncatula.

Arabidopsis clathrin adaptor EPSIN1 but not MODIFIED TRANSPORT TO THE VACOULE1 contributes to effective plant immunity against pathogenic Pseudomonas bacteria

In eukaryotes, EPSINs are Epsin N-terminal Homology (ENTH) domain-containing proteins that serve as monomeric clathrin adaptors at the plasma membrane (PM) or the trans-Golgi Network (TGN)/early endosomes (EE). The model plant Arabidopsis thaliana encodes for seven ENTH proteins, of which so far, only AtEPSIN1 (AtEPS1) and MODIFIED TRANSPORT TO THE VACUOLE1 (AtMTV1) localize to the TGN/EE and contribute to cargo trafficking to both the cell surface and the vacuole. However, relatively little is known about role(s) of any plant EPSIN in governing physiological responses. We have recently shown that AtEPS1 is a positive modulator of plant immune signaling and pattern-triggered immunity against flagellated Pseudomonas syringae pv. tomato (Pto) DC3000 bacteria. In eps1 mutants, impaired immune responses correlate with reduced accumulation of the receptor FLAGELLIN SENSING2 (AtFLS2) and the convergent immune co-receptor BRASSINOSTEROID INSENTIVE1-ASSOCIATED RECEPTOR KINASE1 (AtBAK1) in the PM. Here, we report that in contrast to AtEPS1, the TGN/EE-localized AtMTV1 did not contribute significantly to immunity against pathogenic Pto DC3000 bacteria. We also compared the amino acid sequences, peptide motif structures and in silico tertiary structures of the ENTH domains of AtEPS1 and AtMTV1 in more detail. We conclude that despite sharing the classical tertiary alpha helical ENTH-domain structure and clathrin-binding motifs, the overall low amino acid identity and differences in peptide motifs may explain their role(s) in trafficking of some of the same as well as distinct cargo components to their site of function, with the latter potentially contributing to differences in physiological responses.

Xylan clustering on the pollen surface is required for exine patterning

Xylan is a crosslinking polymer that plays an important role in the assembly of heterogeneous cell wall structures in plants. The pollen wall, a specialized cell wall matrix, exhibits diverse sculpted patterns that serve to protect male gametophytes and facilitate pollination during plant reproduction. However, whether xylan is precisely anchored into clusters and its influence on pollen wall patterning remain unclear. Here, we report xylan clustering on the mature pollen surface in different plant species that is indispensable for the formation of sculpted exine patterns in dicot and monocot plants. Chemical composition analyses revealed that xylan is generally present at low abundance in the mature pollen of flowering plants and shows plentiful variations in terms of substitutions and modifications. Consistent with the expression profiles of their encoding genes, genetic characterization revealed IRREGULAR XYLEM10-LIKE (IRX10L) and its homologous proteins in the GT47 family of glycosyltransferases as key players in the formation of these xylan micro-/nano-compartments on the pollen surface in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa). A deficiency in xylan biosynthesis abolished exine patterning on pollen and compromised male fertility. Therefore, our study outlines a mechanism of exine patterning and provides a tool for manipulating male fertility in crop breeding.

Differing structures of galactoglucomannan in eudicots and non-eudicot angiosperms

The structures of cell wall mannan hemicelluloses have changed during plant evolution. Recently, a new structure called β-galactoglucomannan (β-GGM) was discovered in eudicot plants. This galactoglucomannan has β-(1,2)-Gal-α-(1,6)-Gal disaccharide branches on some mannosyl residues of the strictly alternating Glc-Man backbone. Studies in Arabidopsis revealed β-GGM is related in structure, biosynthesis and function to xyloglucan. However, when and how plants acquired β-GGM remains elusive. Here, we studied mannan structures in many sister groups of eudicots. All glucomannan structures were distinct from β-GGM. In addition, we searched for candidate mannan β-galactosyltransferases (MBGT) in non-eudicot angiosperms. Candidate AtMBGT1 orthologues from rice (OsGT47A-VII) and Amborella (AtrGT47A-VII) did not show MBGT activity in vivo. However, the AtMBGT1 orthologue from rice showed MUR3-like xyloglucan galactosyltransferase activity in complementation analysis using Arabidopsis. Further, reverse genetic analysis revealed that the enzyme (OsGT47A-VII) contributes to proper root growth in rice. Together, gene duplication and diversification of GT47A-VII in eudicot evolution may have been involved in the acquisition of mannan β-galactosyltransferase activity. Our results indicate that β-GGM is likely to be a eudicot-specific mannan.

Laser Capture Microdissection Transcriptome Reveals Spatiotemporal Tissue Gene Expression Patterns of Medicago truncatula Roots Responding to Rhizobia

We report a public resource for examining the spatiotemporal RNA expression of 54,893 Medicago truncatula genes during the first 72 h of response to rhizobial inoculation. Using a methodology that allows synchronous inoculation and growth of more than 100 plants in a single media container, we harvested the same segment of each root responding to rhizobia in the initial inoculation over a time course, collected individual tissues from these segments with laser capture microdissection, and created and sequenced RNA libraries generated from these tissues. We demonstrate the utility of the resource by examining the expression patterns of a set of genes induced very early in nodule signaling, as well as two gene families (CLE peptides and nodule specific PLAT-domain proteins) and show that despite similar whole-root expression patterns, there are tissue differences in expression between the genes. Using a rhizobial response dataset generated from transcriptomics on intact root segments, we also examined differential temporal expression patterns and determined that, after nodule tissue, the epidermis and cortical cells contained the most temporally patterned genes. We circumscribed gene lists for each time and tissue examined and developed an expression pattern visualization tool. Finally, we explored transcriptomic differences between the inner cortical cells that become nodules and those that do not, confirming that the expression of 1-aminocyclopropane-1-carboxylate synthases distinguishes inner cortical cells that become nodules and provide and describe potential downstream genes involved in early nodule cell division. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Genomic analysis and characterization of new loci associated with seed protein and oil content in soybeans

Breeding for increased protein without a reduction in oil content in soybeans [Glycine max (L.) Merr.] is a challenge for soybean breeders but an expected goal. Many efforts have been made to develop new soybean varieties with high yield in combination with desirable protein and/or oil traits. An elite line, R05-1415, was reported to be high yielding, high protein, and low oil. Several significant quantitative trait loci (QTL) for protein and oil were reported in this line, but many of them were unstable across environments or genetic backgrounds. Thus, a new study under multiple field environments using the Infinium BARCSoySNP6K BeadChips was conducted to detect and confirm stable genomic loci for these traits. Genetic analyses consistently detected a single major genomic locus conveying these two traits with remarkably high phenotypic variation explained (R2 ), varying between 24.2% and 43.5%. This new genomic locus is located between 25.0 and 26.7 Mb, distant from the previously reported QTL and did not overlap with other commonly reported QTL and the recently cloned gene Glyma.20G085100. Homolog analysis indicated that this QTL did not result from the paracentric chromosome inversion with an adjacent genomic fragment that harbors the reported QTL. The pleiotropic effect of this QTL could be a challenge for improving protein and oil simultaneously; however, a further study of four candidate genes with significant expressions in the seed developmental stages coupled with haplotype analysis may be able to pinpoint causative genes. The functionality and roles of these genes can be determined and characterized, which lay a solid foundation for the improvement of protein and oil content in soybeans.

Coordinated regulation of the entry and exit steps of aromatic amino acid biosynthesis supports the dual lignin pathway in grasses

Vascular plants direct large amounts of carbon to produce the aromatic amino acid phenylalanine to support the production of lignin and other phenylpropanoids. Uniquely, grasses, which include many major crops, can synthesize lignin and phenylpropanoids from both phenylalanine and tyrosine. However, how grasses regulate aromatic amino acid biosynthesis to feed this dual lignin pathway is unknown. Here we show, by stable-isotope labeling, that grasses produce tyrosine >10-times faster than Arabidopsis without compromising phenylalanine biosynthesis. Detailed in vitro enzyme characterization and combinatorial in planta expression uncovered that coordinated expression of specific enzyme isoforms at the entry and exit steps of the aromatic amino acid pathway enables grasses to maintain high production of both tyrosine and phenylalanine, the precursors of the dual lignin pathway. These findings highlight the complex regulation of plant aromatic amino acid biosynthesis and provide novel genetic tools to engineer the interface of primary and specialized metabolism in plants.

Model-based inference of a plant-specific dual role for HOPS in regulating guard cell vacuole fusion

Stomata are the pores on a leaf surface that regulate gas exchange. Each stoma consists of two guard cells whose movements regulate pore opening and thereby control CO2 fixation and water loss. Guard cell movements depend in part on the remodeling of vacuoles, which have been observed to change from a highly fragmented state to a fused morphology during stomata opening. This change in morphology requires a membrane fusion mechanism that responds rapidly to environmental signals, allowing plants to respond to diurnal and stress cues. With guard cell vacuoles being both large and responsive to external signals, stomata represent a unique system in which to delineate mechanisms of membrane fusion. Fusion of vacuole membranes is a highly conserved process in eukaryotes, with key roles played by two multi-subunit complexes: HOPS (homotypic fusion and vacuolar protein sorting) and SNARE (soluble NSF attachment protein receptor). HOPS is a vacuole tethering factor that is thought to chaperone SNAREs from apposing vacuole membranes into a fusion-competent complex capable of rearranging membranes. To resolve a counter-intuitive observation regarding the role of HOPS in regulating plant vacuole morphology, we derived a quantitative model of vacuole fusion dynamics and used it to generate testable predictions about HOPS-SNARE interactions. We derived our model by applying simulation-based inference to integrate prior knowledge about molecular interactions with limited, qualitative observations of emergent vacuole phenotypes. By constraining the model parameters to yield the emergent outcomes observed for stoma opening - as induced by two distinct chemical treatments - we predicted a dual role for HOPS and identified a stalled form of the SNARE complex that differs from phenomena reported in yeast. We predict that HOPS has contradictory actions at different points in the fusion signaling pathway, promoting the formation of SNARE complexes, but limiting their activity.

Comprehensive and evolutionary analysis of Spodoptera litura-inducible Cytochrome P450 monooxygenase gene family in Glycine max elucidate their role in defense

Plants being sessile organisms and lacking both circulating phagocytic cells and somatic adaptive immune response, have thrived on various defense mechanisms to fend off insect pests and invasion of pathogens. CYP450s are the versatile enzymes, which thwart plants against insect pests by ubiquitous biosynthesis of phytohormones, antioxidants, and secondary metabolites, utilizing them as feeding deterrents and direct toxins. Therefore, a comprehensive analysis of biotic stress-responsive CYPs from Glycine max was performed to ascertain their function against S. litura-infestation. Phylogenetic analysis and evolutionary studies on conserved domains and motifs disclosed the evolutionary correspondence of these GmCYPs with already characterized members of the CYP450 superfamily and close relatedness to Medicago truncatula. These GmCYPs were mapped on 13 chromosomes; they possess 1-8 exons; they have evolved due to duplication and are localized in endoplasmic reticulumn. Further, identification of methyl-jasmonate, salicylic acid, defense responsive and flavonoid biosynthesis regulating cis-acting elements, their interaction with biotic stress regulating proteins and their differential expression in diverse types of tissues, and during herbivory, depicted their responsiveness to biotic stress. Three-dimensional homology modelling of GmCYPs, docking with heme cofactor required for their catalytic activity and enzyme-substrate interactions were performed to understand the functional mechanism of their action. Moreover, to gain insight into their involvement in plant defense, gene expression analysis was evaluated, which revealed differential expression of 11 GmCYPs upon S. litura-infestation, 12 GmCYPs on wounding while foliar spray of ethylene, methyl-jasmonate and salicylic acid differentially regulated 11 GmCYPs, 6 GmCYPs, and 10 GmCYPs respectively. Our study comprehensively analysed the underlying mechanism of GmCYPs function during S. litura-infestation, which can be further utilized for functional characterization to develop new strategies for enhancing soybean resistance to insect pests.

The trehalose 6-phosphate phosphatase family in plants

Trehalose 6-phosphate (Tre6P), the intermediate of trehalose biosynthesis, is an essential signalling metabolite linking plant growth and development to carbon metabolism. While recent work has focused predominantly on the enzymes that produce Tre6P, little is known about the proteins that catalyse its degradation, the trehalose 6-phosphate phosphatases (TPPs). Often occurring in large protein families, TPPs exhibit cell-, tissue- and developmental stage-specific expression patterns, suggesting important regulatory functions in controlling local levels of Tre6P and trehalose as well as Tre6P signalling. Furthermore, growing evidence through gene expression studies and transgenic approaches shows that TPPs play an important role in integrating environmental signals with plant metabolism. This review highlights the large diversity of TPP isoforms in model and crop plants and identifies how modulating Tre6P metabolism in certain cell types, tissues, and at different developmental stages may promote stress tolerance, resilience and increased crop yield.

Identification of Reference Genes for Precise Expression Analysis during Germination in Chenopodium quinoa Seeds under Salt Stress

Chenopodium quinoa Willd. (quinoa), a member of the Amaranthaceae family, is an allotetraploid annual plant, endemic to South America. The plant of C. quinoa presents significant ecological plasticity with exceptional adaptability to several environmental stresses, including salinity. The resilience of quinoa to several abiotic stresses, as well as its nutritional attributes, have led to significant shifts in quinoa cultivation worldwide over the past century. This work first defines germination sensu stricto in quinoa where the breakage of the pericarp and the testa is followed by endosperm rupture (ER). Transcriptomic changes in early seed germination stages lead to unstable expression levels in commonly used reference genes that are typically stable in vegetative tissues. Noteworthy, no suitable reference genes have been previously identified specifically for quinoa seed germination under salt stress conditions. This work aims to identify these genes as a prerequisite step for normalizing qPCR data. To this end, germinating seeds from UDEC2 and UDEC4 accessions, with different tolerance to salt, have been analyzed under conditions of absence (0 mM NaCl) and in the presence (250 mM NaCl) of sodium chloride. Based on the relevant literature, six candidate reference genes, Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), Monensin sensitivity1 (MON1), Polypyrimidine tract-binding protein (PTB), Actin-7 (ACT7), Ubiquitin-conjugating enzyme (UBC), and 18S ribosomal RNA (18S), were selected and assessed for stability using the RefFinder Tool encompassing the statistical algorithms geNorm, NormFinder, BestKeeper, and ΔCt in the evaluation. The data presented support the suitability of CqACT7 and CqUBC as reference genes for normalizing gene expression during seed germination under salinity stress. These recommended reference genes can be valuable tools for consistent qPCR studies on quinoa seeds.

Multi-omics atlas of combinatorial abiotic stress responses in wheat

Field-grown crops rarely experience growth conditions in which yield can be maximized. Environmental stresses occur in combination, with advancements in crop tolerance further complicated by its polygenic nature. Strategic targeting of causal genes is required to meet future crop production needs. Here, we employed a systems biology approach in wheat (Triticum aestivum L.) to investigate physio-metabolic adjustments and transcriptome reprogramming involved in acclimations to heat, drought, salinity and all combinations therein. A significant shift in magnitude and complexity of plant response was evident across stress scenarios based on the agronomic losses, increased proline concentrations and 8.7-fold increase in unique differentially expressed transcripts (DETs) observed under the triple stress condition. Transcriptome data from all stress treatments were assembled into an online, open access eFP browser for visualizing gene expression during abiotic stress. Weighted gene co-expression network analysis revealed 152 hub genes of which 32% contained the ethylene-responsive element binding factor-associated amphiphilic repression (EAR) transcriptional repression motif. Cross-referencing against the 31 DETs common to all stress treatments isolated TaWRKY33 as a leading candidate for greater plant tolerance to combinatorial stresses. Integration of our findings with available literature on gene functional characterization allowed us to further suggest flexible gene combinations for future adaptive gene stacking in wheat. Our approach demonstrates the strength of robust multi-omics-based data resources for gene discovery in complex environmental conditions. Accessibility of such datasets will promote cross-validation of candidate genes across studies and aid in accelerating causal gene validation for crop resiliency.

Cella: 3D data visualization for plant single-cell transcriptomics in Blender

AIMS

Recent advancements in single-cell transcriptomics have facilitated the possibility of acquiring vast amounts of data at single-cell resolution. This development has provided a broader and more comprehensive understanding of complex biological processes. The growing datasets require a visualization tool that transforms complex data into an intuitive representation. To address this challenge, we have utilized an open-source 3D software Blender to design Cella, a cell atlas visualization tool, which transforms data into 3D heatmaps that can be rendered into image libraries. Our tool is designed to support especially research on plant development.

DATA RESOURCES GENERATED

To validate our method, we have created a 3D model representing the Arabidopsis thaliana root meristem and mapped an existing single-cell RNA-seq dataset into the 3D model. This provided a user-friendly visual representation of the expression profiles of 21,489 genes from two perspectives (42,978 images).

UTILITY OF THE RESOURCE

This approach is not limited to single-cell RNA-seq data of the Arabidopsis root meristem. We provide detailed step-by-step instructions to generate 3D models and a script that can be customized to project data onto different tissues.

KEY RESULTS

Our tool provides a proof-of-concept method for how increasingly complex single-cell RNA-seq datasets can be visualized in a simple and cohesive manner.

Wheat Grains as a Sustainable Source of Protein for Health

Protein deficiency is recognized among the major global health issues with an underestimation of its importance. Genetic biofortification is a cost-effective and sustainable strategy to overcome global protein malnutrition. This study was designed to focus on protein-dense grains of wheat (Triticum aestivum L.) and identify the genes governing grain protein content (GPC) that improve end-use quality and in turn human health. Genome-wide association was applied using the 90k iSELECT Infinium and 35k Affymetrix arrays with GPC quantified by using a proteomic-based technique in 369 wheat genotypes over three field-year trials. The results showed significant natural variation among bread wheat genotypes that led to detecting 54 significant quantitative trait nucleotides (QTNs) surpassing the false discovery rate (FDR) threshold. These QTNs showed contrasting effects on GPC ranging from -0.50 to +0.54% that can be used for protein content improvement. Further bioinformatics analyses reported that these QTNs are genomically linked with 35 candidate genes showing high expression during grain development. The putative candidate genes have functions in the binding, remobilization, or transport of protein. For instance, the promising QTN AX-94727470 on chromosome 6B increases GPC by +0.47% and is physically located inside the gene TraesCS6B02G384500 annotated as Trehalose 6-phosphate phosphatase (T6P), which can be employed to improve grain protein quality. Our findings are valuable for the enhancement of protein content and end-use quality in one of the major daily food resources that ultimately improve human nutrition.

Dynamic DNA methylation turnover in gene bodies is associated with enhanced gene expression plasticity in plants

BACKGROUND

In several eukaryotes, DNA methylation occurs within the coding regions of many genes, termed gene body methylation (GbM). Whereas the role of DNA methylation on the silencing of transposons and repetitive DNA is well understood, gene body methylation is not associated with transcriptional repression, and its biological importance remains unclear.

RESULTS

We report a newly discovered type of GbM in plants, which is under constitutive addition and removal by dynamic methylation modifiers in all cells, including the germline. Methylation at Dynamic GbM genes is removed by the DRDD demethylation pathway and added by an unknown source of de novo methylation, most likely the maintenance methyltransferase MET1. We show that the Dynamic GbM state is present at homologous genes across divergent lineages spanning over 100 million years, indicating evolutionary conservation. We demonstrate that Dynamic GbM is tightly associated with the presence of a promoter or regulatory chromatin state within the gene body, in contrast to other gene body methylated genes. We find Dynamic GbM is associated with enhanced gene expression plasticity across development and diverse physiological conditions, whereas stably methylated GbM genes exhibit reduced plasticity. Dynamic GbM genes exhibit reduced dynamic range in drdd mutants, indicating a causal link between DNA demethylation and enhanced gene expression plasticity.

CONCLUSIONS

We propose a new model for GbM in regulating gene expression plasticity, including a novel type of GbM in which increased gene expression plasticity is associated with the activity of DNA methylation writers and erasers and the enrichment of a regulatory chromatin state.

Understanding and applying biological resilience, from genes to ecosystems

The natural world is under unprecedented and accelerating pressure. Much work on understanding resilience to local and global environmental change has, so far, focussed on ecosystems. However, understanding a system's behaviour requires knowledge of its component parts and their interactions. Here we call for increased efforts to understand 'biological resilience', or the processes that enable components across biological levels, from genes to communities, to resist or recover from perturbations. Although ecologists and evolutionary biologists have the tool-boxes to examine form and function, efforts to integrate this knowledge across biological levels and take advantage of big data (e.g. ecological and genomic) are only just beginning. We argue that combining eco-evolutionary knowledge with ecosystem-level concepts of resilience will provide the mechanistic basis necessary to improve management of human, natural and agricultural ecosystems, and outline some of the challenges in achieving an understanding of biological resilience.

Differential effects of arsenite and arsenate on rice (Oryza sativa) plants differing in glutathione S-transferase gene expression

Contamination of paddy soils with arsenic (As) can cause phytotoxicity in rice and increase the accumulation of arsenic in grains. The uptake and accumulation of As in rice depends on the different As species present in the soil. Plants detoxify As by conjugating and sequestering xenobiotic compounds into vacuoles using various enzymes. However, the severity of damage induced by arsenite (As(III)) and arsenate (As(V)), as well as the roles of glutathione S-transferase in detoxifying these As species in rice, are not fully understood. In this study, we developed plant materials overexpressing a glutathione S-transferase gene OsGSTU40 under the control of the maize UBIL promoter. Through systematic investigations of both wild-type Nipponbare (Oryza sativa L., ssp. japonica) and OsGSTU40 overexpression lines under chronic or acute stress of As, we aimed to understand the toxic effects of both As(III) and As(V) on rice plants at the vegetative growth stage. We hypothesized that (i) As(III) and As(V) have different toxic effects on rice plants and (ii) OsGSTU40 played positive roles in As toxicity tolerance. Our results showed that As(III) was more detrimental to plant growth than As(V) in terms of plant growth, biomass, and lipid peroxidation in both chronic and acute exposure. Furthermore, overexpression of OsGSTU40 led to better plant growth even though uptake of As(V), but not As(III), into shoots was enhanced in transgenic plants. In acute As(III) stress, transgenic plants exhibited a lower level of lipid peroxidation than wild-type plants. The element composition of plants was dominated by the different As stress treatments rather than by the genotype, while the As concentration was negatively correlated with phosphorus and silicon. Overall, our findings suggest that As(III) is more toxic to plants than As(V) and that glutathione S-transferase OsGSTU40 differentially affects plant reactions and tolerance to different species of arsenic.

TM3 and STM3 Promote Flowering Together with FUL2 and MBP20, but Act Antagonistically in Inflorescence Branching in Tomato

The moment at which a plant transitions to reproductive development is paramount to its life cycle and is strictly controlled by many genes. The transcription factor SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) plays a central role in this process in Arabidopsis. However, the role of SOC1 in tomato (Solanum lycopersicum) has been sparsely studied. Here, we investigated the function of four tomato SOC1 homologs in the floral transition and inflorescence development. We thoroughly characterized the SOC1-like clade throughout the Solanaceae and selected four tomato homologs that are dynamically expressed upon the floral transition. We show that of these homologs, TOMATO MADS 3 (TM3) and SISTER OF TM3 (STM3) promote the primary and sympodial transition to flowering, while MADS-BOX PROTEIN 23 (MBP23) and MBP18 hardly contribute to flowering initiation in the indeterminate cultivar Moneyberg. Protein-protein interaction assays and whole-transcriptome analysis during reproductive meristem development revealed that TM3 and STM3 interact and share many targets with FRUITFULL (FUL) homologs, including cytokinin regulators. Furthermore, we observed that mutating TM3/STM3 affects inflorescence development, but counteracts the inflorescence-branching phenotype of ful2 mbp20. Collectively, this indicates that TM3/STM3 promote the floral transition together with FUL2/MBP20, while these transcription factors have opposite functions in inflorescence development.