Arabidopsis REI-LIKE proteins activate ribosome biogenesis during cold acclimation

Dynamics of ribosome composition and ribosomal protein phosphorylation in immune signaling in Arabidopsis thaliana

In plants, the detection of microbe-associated molecular patterns (MAMPs) induces primary innate immunity by the activation of mitogen-activated protein kinases (MAPKs). We show here that the MAMP-activated MAPK MPK6 not only modulates defense through transcriptional regulation but also via the ribosomal protein translation machinery. To understand the effects of MPK6 on ribosomes and their constituent ribosomal proteins (RPs), polysomes, monosomes and the phosphorylation status of the RPs, MAMP-treated WT and mpk6 mutant plants were analysed. MAMP-activation induced rapid changes in RP composition of monosomes, polysomes and in the 60S ribosomal subunit in an MPK6-specific manner. Phosphoproteome analysis showed that MAMP-activation of MPK6 regulates the phosphorylation status of the P-stalk ribosomal proteins by phosphorylation of RPP0 and the concomitant dephosphorylation of RPP1 and RPP2. These events coincide with a significant decrease in the abundance of ribosome-bound RPP0s, RPP1s and RPP3s in polysomes. The P-stalk is essential in regulating protein translation by recruiting elongation factors. Accordingly, we found that RPP0C mutant plants are compromised in basal resistance to Pseudomonas syringae infection. These data suggest that MAMP-induced defense also involves MPK6-induced regulation of P-stalk proteins, highlighting a new role of ribosomal regulation in plant innate immunity.

Cytosolic RGG RNA-binding proteins are temperature sensitive flowering time regulators in Arabidopsis

mRNA translation is tightly regulated by various classes of RNA-binding proteins (RBPs) during development and in response to changing environmental conditions. In this study, we characterize the arginine-glycine-glycine (RGG) motif containing RBP family of Arabidopsis thaliana representing homologues of the multifunctional translation regulators and ribosomal preservation factors Stm1 from yeast (ScStm1) and human SERBP1 (HsSERBP1). The Arabidopsis genome encodes three RGG proteins named AtRGGA, AtRGGB and AtRGGC. While AtRGGA is ubiquitously expressed, AtRGGB and AtRGGC are enriched in dividing cells. All AtRGGs localize almost exclusively to the cytoplasm and bind with high affinity to ssRNA, while being capable to interact with most nucleic acids, except dsRNA. A protein-interactome study shows that AtRGGs interact with ribosomal proteins and proteins involved in RNA processing and transport. In contrast to ScStm1, AtRGGs are enriched in ribosome-free fractions in polysome profiles, suggesting additional plant-specific functions. Mutant studies show that AtRGG proteins differentially regulate flowering time, with a distinct and complex temperature dependency for each AtRGG protein. In conclusion, we suggest that AtRGGs function in fine-tuning translation efficiency to control flowering time and potentially other developmental processes in response to environmental changes.

Transcriptomic analysis in tomato fruit reveals divergences in genes involved in cold stress response and fruit ripening

Cold storage is widely used to extend the postharvest life of most horticultural crops, including tomatoes, but this practice triggers cold stress and leads to the development of undesirable chilling injury (CI) symptoms. The underlying mechanisms of cold stress response and CI development in fruits remain unclear as they are often intermingled with fruit ripening changes. To gain insight into cold responses in fruits, we examined the effect of the potent ethylene signaling inhibitor 1-methylcyclopropene (1-MCP) on fruit ripening, CI occurrence and gene expression in mature green tomatoes during storage at 20°C and 5°C. 1-MCP treatments effectively inhibited ethylene production and peel color changes during storage at 20°C. Storage at 5°C also inhibited both ethylene production and peel color change; during rewarming at 20°C, 1-MCP treatments inhibited peel color change but failed to inhibit ethylene production. Furthermore, fruits stored at 5°C for 14 d developed CI symptoms (surface pitting and decay) during the rewarming period at 20°C regardless of 1-MCP treatment. Subsequent RNA-Seq analysis revealed that cold stress triggers a large-scale transcriptomic adjustment, as noticeably more genes were differentially expressed at 5°C (8,406) than at 20°C (4,814). More importantly, we have found some important divergences among genes involved in fruit ripening (up- or down-regulated at 20°C; inhibited by 1-MCP treatment) and those involved in cold stress (up- or down-regulated at 5°C; unaffected by 1-MCP treatment). Transcriptomic adjustments unique to cold stress response were associated with ribosome biogenesis, NcRNA metabolism, DNA methylation, chromatin formation/remodeling, and alternative splicing events. These data should foster further research into cold stress response mechanisms in fruits with the ultimate aim of improving tolerance to low temperature and reduction of CI symptoms during cold storage.

Flower transcriptional response to long term hot and cold environments in Antirrhinum majus

Short term experiments have identified heat shock and cold response elements in many biological systems. However, the effect of long-term low or high temperatures is not well documented. To address this gap, we grew Antirrhinum majus plants from two-weeks old until maturity under control (normal) (22/16°C), cold (15/5°C), and hot (30/23°C) conditions for a period of two years. Flower size, petal anthocyanin content and pollen viability obtained higher values in cold conditions, decreasing in middle and high temperatures. Leaf chlorophyll content was higher in cold conditions and stable in control and hot temperatures, while pedicel length increased under hot conditions. The control conditions were optimal for scent emission and seed production. Scent complexity was low in cold temperatures. The transcriptomic analysis of mature flowers, followed by gene enrichment analysis and CNET plot visualization, showed two groups of genes. One group comprised genes controlling the affected traits, and a second group appeared as long-term adaptation to non-optimal temperatures. These included hypoxia, unsaturated fatty acid metabolism, ribosomal proteins, carboxylic acid, sugar and organic ion transport, or protein folding. We found a differential expression of floral organ identity functions, supporting the flower size data. Pollinator-related traits such as scent and color followed opposite trends, indicating an equilibrium for rendering the organs for pollination attractive under changing climate conditions. Prolonged heat or cold cause structural adaptations in protein synthesis and folding, membrane composition, and transport. Thus, adaptations to cope with non-optimal temperatures occur in basic cellular processes.

Solanum tuberosum Microtuber Development under Darkness Unveiled through RNAseq Transcriptomic Analysis

Potato microtuber (MT) development through in vitro techniques are ideal propagules for producing high quality potato plants. MT formation is influenced by several factors, i.e., photoperiod, sucrose, hormones, and osmotic stress. We have previously developed a protocol of MT induction in medium with sucrose (8% w/v), gelrite (6g/L), and 2iP as cytokinin under darkness. To understand the molecular mechanisms involved, we performed a transcriptome-wide analysis. Here we show that 1715 up- and 1624 down-regulated genes were involved in this biological process. Through the protein-protein interaction (PPI) network analyses performed in the STRING database (v11.5), we found 299 genes tightly associated in 14 clusters. Two major clusters of up-regulated proteins fundamental for life growth and development were found: 29 ribosomal proteins (RPs) interacting with 6 PEBP family members and 117 cell cycle (CC) proteins. The PPI network of up-regulated transcription factors (TFs) revealed that at least six TFs-MYB43, TSF, bZIP27, bZIP43, HAT4 and WOX9-may be involved during MTs development. The PPI network of down-regulated genes revealed a cluster of 83 proteins involved in light and photosynthesis, 110 in response to hormone, 74 in hormone mediate signaling pathway and 22 related to aging.

Regulatory dynamics of gene expression in the developing male gametophyte of Arabidopsis

Sexual reproduction in angiosperms requires the production and delivery of two male gametes by a three-celled haploid male gametophyte. This demands synchronized gene expression in a short developmental window to ensure double fertilization and seed set. While transcriptomic changes in developing pollen are known for Arabidopsis, no studies have integrated RNA and proteomic data in this model. Further, the role of alternative splicing has not been fully addressed, yet post-transcriptional and post-translational regulation may have a key role in gene expression dynamics during microgametogenesis. We have refined and substantially updated global transcriptomic and proteomic changes in developing pollen for two Arabidopsis accessions. Despite the superiority of RNA-seq over microarray-based platforms, we demonstrate high reproducibility and comparability. We identify thousands of long non-coding RNAs as potential regulators of pollen development, hundreds of changes in alternative splicing and provide insight into mRNA translation rate and storage in developing pollen. Our analysis delivers an integrated perspective of gene expression dynamics in developing Arabidopsis pollen and a foundation for studying the role of alternative splicing in this model.

Cryo-EM structure and rRNA modification sites of a plant ribosome

Protein synthesis in crop plants contributes to the balance of food and fuel on our planet, which influences human metabolic activity and lifespan. Protein synthesis can be regulated with respect to changing environmental cues via the deposition of chemical modifications into rRNA. Here, we present the structure of a plant ribosome from tomato and a quantitative mass spectrometry analysis of its rRNAs. The study reveals fine features of the ribosomal proteins and 71 plant-specific rRNA modifications, and it re-annotates 30 rRNA residues in the available sequence. At the protein level, isoAsp is found in position 137 of uS11, and a zinc finger previously believed to be universal is missing from eL34, suggesting a lower effect of zinc deficiency on protein synthesis in plants. At the rRNA level, the plant ribosome differs markedly from its human counterpart with respect to the spatial distribution of modifications. Thus, it represents an additional layer of gene expression regulation, highlighting the molecular signature of a plant ribosome. The results provide a reference model of a plant ribosome for structural studies and an accurate marker for molecular ecology.

Gene regulatory network inference methodology for genomic and transcriptomic data acquired in genetically related heterozygote individuals

MOTIVATION

Inferring gene regulatory networks in non-independent genetically related panels is a methodological challenge. This hampers evolutionary and biological studies using heterozygote individuals such as in wild sunflower populations or cultivated hybrids.

RESULTS

First, we simulated 100 datasets of gene expressions and polymorphisms, displaying the same gene expression distributions, heterozygosities and heritabilities as in our dataset including 173 genes and 353 genotypes measured in sunflower hybrids. Secondly, we performed a meta-analysis based on six inference methods [least absolute shrinkage and selection operator (Lasso), Random Forests, Bayesian Networks, Markov Random Fields, Ordinary Least Square and fast inference of networks from directed regulation (Findr)] and selected the minimal density networks for better accuracy with 64 edges connecting 79 genes and 0.35 area under precision and recall (AUPR) score on average. We identified that triangles and mutual edges are prone to errors in the inferred networks. Applied on classical datasets without heterozygotes, our strategy produced a 0.65 AUPR score for one dataset of the DREAM5 Systems Genetics Challenge. Finally, we applied our method to an experimental dataset from sunflower hybrids. We successfully inferred a network composed of 105 genes connected by 106 putative regulations with a major connected component.

AVAILABILITY AND IMPLEMENTATION

Our inference methodology dedicated to genomic and transcriptomic data is available at https://forgemia.inra.fr/sunrise/inference_methods.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Engineering Ribosomes to Alleviate Abiotic Stress in Plants: A Perspective

As the centerpiece of the biomass production process, ribosome activity is highly coordinated with environmental cues. Findings revealing ribosome subgroups responsive to adverse conditions suggest this tight coordination may be grounded in the induction of variant ribosome compositions and the differential translation outcomes they might produce. In this perspective, we go through the literature linking ribosome heterogeneity to plants’ abiotic stress response. Once unraveled, this crosstalk may serve as the foundation of novel strategies to custom cultivars tolerant to challenging environments without the yield penalty.

Phosphorus Acquisition and Utilization in Plants

Tremendous progress has been made on molecular aspects of plant phosphorus (P) nutrition, often without heeding information provided by soil scientists, ecophysiologists, and crop physiologists. This review suggests ways to integrate information from different disciplines. When soil P availability is very low, P-mobilizing strategies are more effective than mycorrhizal strategies. Soil parameters largely determine how much P roots can acquire from P-impoverished soil, and kinetic properties of P transporters are less important. Changes in the expression of P transporters avoid P toxicity. Plants vary widely in photosynthetic P-use efficiency, photosynthesis per unit leaf P. The challenge is to discover what the trade-offs are of different patterns of investment in P fractions. Less investment may save P, but are costs incurred? Are these costs acceptable for crops? These questions can be resolved only by the concerted action of scientists working at both molecular and physiological levels, rather than pursuing these problems independently.

Dynamic Changes in Seed Germination under Low-Temperature Stress in Maize

Low-temperature stress delays seed germination in maize. Different maize inbred lines display various low-temperature resistance, but the dynamic changes in seed germination under low-temperature stress in maize remain unknown, especially at the transcriptome level. In this study, low-temperature-resistant maize (RM) inbred line 04Qun0522-1-1 had a significantly faster germination speed than low-temperature-sensitive maize (SM) line B283-1 under low-temperature stress. Moreover, the total antioxidant capacity, superoxide dismutase, and peroxidase activities were notably higher in the RM line than in the SM line from 3 to 6 d. In contrast, the SM line showed significantly higher malondialdehyde (MDA) content than the RM line at 6 d. Gene ontology (GO) enrichment analysis showed that in 2dvs0d, both SM and RM lines displayed the downregulation of ribosome-related genes. Moreover, photosystem II and heat shock protein binding-related genes were also downregulated in the SM line. In 4dvs2d, the RM line showed a higher degree of upregulation of the ribosome and peroxidase (POD)-related genes than the SM line. In 6dvs4d, POD-related genes were continuously upregulated in both SM and RM lines, but the degree of upregulation of the genes was higher in the SM line than in the RM line. Moreover, vitamin B6-related genes were specifically upregulated in the RM line. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that in 6dvs4d, phenylpropanoid biosynthesis was the most significantly enriched pathway in both SM and RM lines. Moreover, phenylpropanoid biosynthesis was also enriched in the RM line in 4dvs2d. More than half of the differentially expressed genes (DEGs) in phenylpropanoid biosynthesis were peroxidase, and the DEGs were similar to the GO enrichment analysis. The results provide new insights into maize seed germination in response to low-temperature stress.

Pollen Number and Ribosome Gene Expression Altered in a Genome-Editing Mutant of REDUCED POLLEN NUMBER1 Gene

The number of pollen grains varies within and between species. However, little is known about the molecular basis of this quantitative trait, in contrast with the many studies available on cell differentiation in the stamen. Recently, the first gene responsible for pollen number variation, REDUCED POLLEN NUMBER1 (RDP1), was isolated by genome-wide association studies of Arabidopsis thaliana and exhibited the signature of natural selection. This gene encodes a homolog of yeast Mrt4 (mRNA turnover4), which is an assembly factor of the large ribosomal subunit. However, no further data were available to link ribosome function to pollen development. Here, we characterized the RDP1 gene using the standard A. thaliana accession Col-0. The frameshift mutant, rdp1-3 generated by CRISPR/Cas9 revealed the pleiotropic effect of RDP1 in flowering, thus demonstrating that this gene is required for a broad range of processes other than pollen development. We found that the natural Col-0 allele conferred a reduced pollen number against the Bor-4 allele, as assessed using the quantitative complementation test, which is more sensitive than transgenic experiments. Together with a historical recombination event in Col-0, which was identified by sequence alignment, these results suggest that the coding sequence of RDP1 is the candidate region responsible for the natural phenotypic variation. To elucidate the biological processes in which RDP1 is involved, we conducted a transcriptome analysis. We found that genes responsible for ribosomal large subunit assembly/biogenesis were enriched among the differentially regulated genes, which supported the hypothesis that ribosome biogenesis is disturbed in the rdp1-3 mutant. Among the pollen-development genes, three key genes encoding basic helix-loop-helix (bHLH) transcription factors (ABORTED MICROSPORES (AMS), bHLH010, and bHLH089), as well as direct downstream genes of AMS, were downregulated in the rdp1-3 mutant. In summary, our results suggest a specialized function of ribosomes in pollen development through RDP1, which harbors natural variants under selection.

Fine-Tuning Cold Stress Response Through Regulated Cellular Abundance and Mechanistic Actions of Transcription Factors

Inflictions caused by cold stress can result in disastrous effects on the productivity and survival of plants. Cold stress response in plants requires crosstalk between multiple signaling pathways including cold, heat, and reactive oxygen species (ROS) signaling networks. CBF, MYB, bHLH, and WRKY families are among the TFs that function as key players in the regulation of cold stress response at the molecular level. This review discusses some of the latest understanding on the regulation of expression and the mechanistic actions of plant TFs to address cold stress response. It was shown that the plant response consists of early and late responses as well as memory reprogramming for long-term protection against cold stress. The regulatory network can be differentiated into CBF-dependent and independent pathways involving different sets of TFs. Post-transcriptional regulation by miRNAs, control during ribosomal translation process, and post-translational regulation involving 26S proteosomic degradation are processes that affect the cellular abundance of key regulatory TFs, which is an important aspect of the regulation for cold acclimation. Therefore, fine-tuning of the regulation by TFs for adjusting to the cold stress condition involving the dynamic action of protein kinases, membrane ion channels, adapters, and modifiers is emphasized in this review.

[Formula: see text]: ComplexOme-Structural Network Interpreter used to study spatial enrichment in metazoan ribosomes

BACKGROUND

Upon environmental stimuli, ribosomes are surmised to undergo compositional rearrangements due to abundance changes among proteins assembled into the complex, leading to modulated structural and functional characteristics. Here, we present the ComplexOme-Structural Network Interpreter ([Formula: see text]), a computational method to allow testing whether ribosomal proteins (rProteins) that exhibit abundance changes under specific conditions are spatially confined to particular regions within the large ribosomal complex.

RESULTS

[Formula: see text] translates experimentally determined structures into graphs, with nodes representing proteins and edges the spatial proximity between them. In its first implementation, [Formula: see text] considers rProteins and ignores rRNA and other objects. Spatial regions are defined using a random walk with restart methodology, followed by a procedure to obtain a minimum set of regions that cover all proteins in the complex. Structural coherence is achieved by applying weights to the edges reflecting the physical proximity between purportedly contacting proteins. The weighting probabilistically guides the random-walk path trajectory. Parameter tuning during region selection provides the option to tailor the method to specific biological questions by yielding regions of different sizes with minimum overlaps. In addition, other graph community detection algorithms may be used for the [Formula: see text] workflow, considering that they yield different sized, non-overlapping regions. All tested algorithms result in the same node kernels under equivalent regions. Based on the defined regions, available abundance change information of proteins is mapped onto the graph and subsequently tested for enrichment in any of the defined spatial regions. We applied [Formula: see text] to the cytosolic ribosome structures of Saccharomyces cerevisiae, Oryctolagus cuniculus, and Triticum aestivum using datasets with available quantitative protein abundance change information. We found that in yeast, substoichiometric rProteins depleted from translating polysomes are significantly constrained to a ribosomal region close to the tRNA entry and exit sites.

CONCLUSIONS

[Formula: see text] offers a computational method to partition multi-protein complexes into structural regions and a statistical approach to test for spatial enrichments of any given subsets of proteins. [Formula: see text] is applicable to any multi-protein complex given appropriate structural and abundance-change data. [Formula: see text] is publicly available as a GitHub repository https://github.com/MSeidelFed/COSNet_i and can be installed using the python installer pip.

Spatially Enriched Paralog Rearrangements Argue Functionally Diverse Ribosomes Arise during Cold Acclimation in Arabidopsis

Ribosome biogenesis is essential for plants to successfully acclimate to low temperature. Without dedicated steps supervising the 60S large subunits (LSUs) maturation in the cytosol, e.g., Rei-like (REIL) factors, plants fail to accumulate dry weight and fail to grow at suboptimal low temperatures. Around REIL, the final 60S cytosolic maturation steps include proofreading and assembly of functional ribosomal centers such as the polypeptide exit tunnel and the P-Stalk, respectively. In consequence, these ribosomal substructures and their assembly, especially during low temperatures, might be changed and provoke the need for dedicated quality controls. To test this, we blocked ribosome maturation during cold acclimation using two independent reil double mutant genotypes and tested changes in their ribosomal proteomes. Additionally, we normalized our mutant datasets using as a blank the cold responsiveness of a wild-type Arabidopsis genotype. This allowed us to neglect any reil-specific effects that may happen due to the presence or absence of the factor during LSU cytosolic maturation, thus allowing us to test for cold-induced changes that happen in the early nucleolar biogenesis. As a result, we report that cold acclimation triggers a reprogramming in the structural ribosomal proteome. The reprogramming alters the abundance of specific RP families and/or paralogs in non-translational LSU and translational polysome fractions, a phenomenon known as substoichiometry. Next, we tested whether the cold-substoichiometry was spatially confined to specific regions of the complex. In terms of RP proteoforms, we report that remodeling of ribosomes after a cold stimulus is significantly constrained to the polypeptide exit tunnel (PET), i.e., REIL factor binding and functional site. In terms of RP transcripts, cold acclimation induces changes in RP families or paralogs that are significantly constrained to the P-Stalk and the ribosomal head. The three modulated substructures represent possible targets of mechanisms that may constrain translation by controlled ribosome heterogeneity. We propose that non-random ribosome heterogeneity controlled by specialized biogenesis mechanisms may contribute to a preferential or ultimately even rigorous selection of transcripts needed for rapid proteome shifts and successful acclimation.

Membrane-Enriched Proteomics Link Ribosome Accumulation and Proteome Reprogramming With Cold Acclimation in Barley Root Meristems

Due to their sessile nature, plants rely on root systems to mediate many biotic and abiotic cues. To overcome these challenges, the root proteome is shaped to specific responses. Proteome-wide reprogramming events are magnified in meristems due to their active protein production. Using meristems as a test system, here, we study the major rewiring that plants undergo during cold acclimation. We performed tandem mass tag-based bottom-up quantitative proteomics of two consecutive segments of barley seminal root apexes subjected to suboptimal temperatures. After comparing changes in total and ribosomal protein (RP) fraction-enriched contents with shifts in individual protein abundances, we report ribosome accumulation accompanied by an intricate translational reprogramming in the distal apex zone. Reprogramming ranges from increases in ribosome biogenesis to protein folding factors and suggests roles for cold-specific RP paralogs. Ribosome biogenesis is the largest cellular investment; thus, the vast accumulation of ribosomes and specific translation-related proteins during cold acclimation could imply a divergent ribosomal population that would lead to a proteome shift across the root. Consequently, beyond the translational reprogramming, we report a proteome rewiring. First, triggered protein accumulation includes spliceosome activity in the root tip and a ubiquitous upregulation of glutathione production and S-glutathionylation (S-GSH) assemblage machineries in both root zones. Second, triggered protein depletion includes intrinsically enriched proteins in the tip-adjacent zone, which comprise the plant immune system. In summary, ribosome and translation-related protein accumulation happens concomitantly to a proteome reprogramming in barley root meristems during cold acclimation. The cold-accumulated proteome is functionally implicated in feedbacking transcript to protein translation at both ends and could guide cold acclimation.

Genetic Architecture Underlying the Metabolites of Chlorogenic Acid Biosynthesis in Populus tomentosa

Chlorogenic acid (CGA) plays a crucial role in defense response, immune regulation, and the response to abiotic stress in plants. However, the genetic regulatory network of CGA biosynthesis pathways in perennial plants remains unclear. Here, we investigated the genetic architecture for CGA biosynthesis using a metabolite-based genome-wide association study (mGWAS) and expression quantitative trait nucleotide (eQTN) mapping in a population of 300 accessions of Populus tomentosa. In total, we investigated 204 SNPs which were significantly associated with 11 metabolic traits, corresponding to 206 genes, and were mainly involved in metabolism and cell growth processes of P. tomentosa. We identified 874 eQTNs representing 1066 genes, in which the expression and interaction of causal genes affected phenotypic variation. Of these, 102 genes showed significant signatures of selection in three geographical populations, which provided insights into the adaptation of CGA biosynthesis to the local environment. Finally, we constructed a genetic network of six causal genes that coordinately regulate CGA biosynthesis, revealing the multiple regulatory patterns affecting CGA accumulation in P. tomentosa. Our study provides a multiomics strategy for understanding the genetic basis underlying the natural variation in the CGA biosynthetic metabolites of Populus, which will enhance the genetic development of abiotic-resistance varieties in forest trees.