Ribosome Loading Research Articles

Machine Learning for Designing Next-Generation mRNA Therapeutics.

Over just the last 2 years, mRNA therapeutics and vaccines have undergone a rapid transition from an intriguing concept to real-world impact. However, whereas some aspects of mRNA therapeutics, such as the use of chemical modifications to increase stability and reduce immunogenicity, have been extensively optimized for over two decades, other aspects, particularly the selection and design of the noncoding leader and trailer sequences which control translation efficiency and stability, have received comparably less attention. In practice, such 5' and 3' untranslated regions (UTRs) are often borrowed from highly expressed human genes with few or no modifications, as in the case for the Pfizer/BioNTech Covid vaccine. Focusing on the 5'UTR, we here argue that model-driven design is a promising alternative that provides unprecedented control over 5'UTR function. We review recent work that combines synthetic biology with machine learning to build quantitative models that relate ribosome loading, and thus translation efficiency, to the 5'UTR sequence. We first introduce an experimental approach that uses polysome profiling and high-throughput sequencing to quantify ribosome loading for hundreds of thousands of 5'UTRs in parallel. We apply this approach to measure ribosome loading in synthetic RNA libraries with a random sequence inserted into the 5'UTR. We then review Optimus 5-Prime, a convolutional neural network model trained on the experimental data. We highlight that very accurate models of biological regulation can be learned from synthetic data sets with degenerate 5'UTRs. We validate model predictions not only on held-out data sets from our random library but also on a large library of over 30 000 human 5'UTR fragments and using translation reporter data collected independently by other groups. Both the experiment and model are compatible with commonly used chemically modified nucleosides, in particular, pseudouridine (Ψ) and 1-methyl-pseudouridine (m1Ψ). We find that, in general, 5'UTRs have very similar impacts when combined with different protein-coding sequences and even in the context of different chemical modifications. We demonstrate that Optimus 5-Prime can be combined with design algorithms to generate de novo sequences with precisely defined translation efficiencies. We emphasize recent developments in design algorithms that rely on activation maximization and generative modeling to improve both the fitness and diversity of designed sequences. Compared with prior approaches such as genetic algorithms, we show that these approaches are not only faster but also less likely to get stuck in local sequence optima. Finally, we discuss how the approach reviewed here can be generalized to other gene regions and applications.

Combining an optimized mRNA template with a double purification process allows strong expression of in vitro transcribed mRNA

mRNA is a blooming technology for vaccination and has gained global attention during the SARS-CoV-2 pandemic. However, the recent clinical trials have highlighted increased reactogenicity when using high mRNA doses. Intending to increase the potency of mRNA therapeutics and to decrease the therapeutic dose, we designed a mRNA backbone and optimized the mRNA purification process. We used the enhanced green fluorescent protein (eGFP) reporter gene flanked by one 5′ untranslated region (UTR) and two 3′ UTRs of the human β-globin as a reference mRNA and identified the most promising mRNA sequence using in vitro and in vivo models. First, we assessed the impact of different poly(A) sizes on translation and selected the most optimal sequence. Then, we selected the best 5′ UTR among synthetic sequences displaying a high ribosome loading. Finally, we evaluated the transfection efficiency of our standard mRNA template after two capping strategies and purification using either double-stranded RNA (dsRNA) depletion or dephosphorylation of 5′PPP RNA or both combined. Double purification was shown to give the best results. Altogether, the use of a newly defined 5′ UTR coupled to post-transcriptional treatments will be of great interest in the mRNA vaccine field, by limiting the amount of the antigen-coding transcript and subsequently the formulation components needed for an efficient vaccination.

Utr.annotation: a tool for annotating genomic variants that could influence post-transcriptional regulation.

Whole genome sequencing of patient populations is identifying thousands of new variants in untranslated regions (UTRs). While the consequences of UTR mutations are not as easily predicted from primary sequence as coding mutations are, there are some known features of UTRs that modulate their function. utr.annotation is an R package that can be used to annotate potential deleterious variants in the UTR regions for both human and mouse species. Given a CSV or VCF format variant file, utr.annotation provides information of each variant on whether and how it alters known translational regulators including upstream open reading frames, upstream Kozak sequences, polyA signals, Kozak sequences at the annotated translation start site, start codons and stop codons, conservation scores in the variant position, and whether and how it changes ribosome loading based on a model derived from empirical data. utr.annotation is freely available on Bitbucket (https://bitbucket.org/jdlabteam/utr.annotation/src/master/) and CRAN (https://cran.r-project.org/web/packages/utr.annotation/index.html). Supplementary data are available at Bioinformatics online.

Nascent Ribo-Seq measures ribosomal loading time and reveals kinetic impact on ribosome density.

In general, mRNAs are assumed to be loaded with ribosomes instantly upon entry into the cytoplasm. To measure ribosome density (RD) on nascent mRNA, we developed nascent Ribo-Seq by combining Ribo-Seq with progressive 4-thiouridine labeling. In mouse macrophages, we determined experimentally the lag between the appearance of nascent mRNA and its association with ribosomes, which was calculated to be 20-22 min for bulk mRNA. In mouse embryonic stem cells, nRibo-Seq revealed an even stronger lag of 35-38 min in ribosome loading. After stimulation of macrophages with lipopolysaccharide, the lag between cytoplasmic and translated mRNA leads to uncoupling between input and ribosome-protected fragments, which gives rise to distorted RD measurements under conditions where mRNA amounts are far from steady-state expression. As a result, we demonstrate that transcriptional changes affect RD in a passive way.

Predicting mean ribosome load for 5'UTR of any length using deep learning.

The 5' untranslated region plays a key role in regulating mRNA translation and consequently protein abundance. Therefore, accurate modeling of 5'UTR regulatory sequences shall provide insights into translational control mechanisms and help interpret genetic variants. Recently, a model was trained on a massively parallel reporter assay to predict mean ribosome load (MRL)-a proxy for translation rate-directly from 5'UTR sequence with a high degree of accuracy. However, this model is restricted to sequence lengths investigated in the reporter assay and therefore cannot be applied to the majority of human sequences without a substantial loss of information. Here, we introduced frame pooling, a novel neural network operation that enabled the development of an MRL prediction model for 5'UTRs of any length. Our model shows state-of-the-art performance on fixed length randomized sequences, while offering better generalization performance on longer sequences and on a variety of translation-related genome-wide datasets. Variant interpretation is demonstrated on a 5'UTR variant of the gene HBB associated with beta-thalassemia. Frame pooling could find applications in other bioinformatics predictive tasks. Moreover, our model, released open source, could help pinpoint pathogenic genetic variants.

95233 Analysis of 5'UTR Variation in Rare Disease Patients Reveals Variants of Potential Disease Relevance

ABSTRACT IMPACT: This work sheds diagnostic insight on patients with idiopathic rare disease and has the potential to further their care and treatment as a result. OBJECTIVES/GOALS: Correct diagnosis is imperative to treating patients with idiopathic, suspected genetic conditions, yet sequencing approaches leave up to 70% of these patients undiagnosed. We sought to improve diagnosis rates for a cohort patients referred for sequencing by characterizing deleterious variants within the 5’UTR. METHODS/STUDY POPULATION: We retrospectively analyzed whole exome sequencing (WES) data from 472 unsolved rare disease patients within the Mayo Clinic Center for Individualized Medicine to identify variants within the 5’UTR that affect the presence of upstream open reading frames (uORFs). uORFs are short regions (typically 30bp - 600bp) that typically influence downstream gene translation by sequestering ribosomes. We specifically searched for variants with the potential to disrupt existing uORFs or introduce new uORFs within the 5’UTR, and developed a pipeline to annotate these variants with information including GnomAD allele frequency and gene loss of function intolerance (pLI) score. To aid in variant interpretation, we applied two deep learning tools to predict variant impacts on transcript ribosome load (TITER and FramePool). RESULTS/ANTICIPATED RESULTS: Our pipeline identified a median of 21 variants per patient that were predicted to have a deleterious impact on the translational efficiency of protein coding transcripts, primarily by introducing new start codons within the 5’UTR or by altering the Kozak consensus of existing start codons. A median of 10 of these variants occur upstream of haploinsufficient genes with an existing disease association. We also identified a subset of variants that are predicted to introduce translationally active N-terminal extensions to protein coding transcripts, with the potential to disrupt protein localization and processing. DISCUSSION/SIGNIFICANCE OF FINDINGS: This work demonstrates that analysis of 5’UTR variants can be incorporated into existing WES pipelines, and identifies a group of variants with potential significance to patient disease. Further experimental evidence is necessary to ascertain the pathogenicity of these variants.

Live-cell imaging reveals kinetic determinants of quality control triggered by ribosome stalling

Translation of problematic mRNA sequences induces ribosome stalling, triggering quality-control events, including ribosome rescue and nascent polypeptide degradation. To define the timing and regulation of these processes, we developed a SunTag-based reporter to monitor translation of a problematic sequence (poly[A]) in real time on single mRNAs. Although poly(A)-containing mRNAs undergo continuous translation over the timescale of minutes to hours, ribosome load is increased by ∼3-fold compared to a control, reflecting long queues of ribosomes extending far upstream of the stall. We monitor the resolution of these queues in real time and find that ribosome rescue is very slow compared to both elongation and termination. Modulation of pause strength, collision frequency, and the collision sensor ZNF598 reveals how the dynamics of ribosome collisions and their recognition facilitate selective targeting for quality control. Our results establish that slow clearance of stalled ribosomes allows cells to distinguish between transient and deleterious stalls.

Preferential Ribosome Loading on the Stress-Upregulated mRNA Pool Shapes the Selective Translation under Stress Conditions.

The reprogramming of gene expression is one of the key responses to environmental stimuli, whereas changes in mRNA do not necessarily bring forth corresponding changes of the protein, which seems partially due to the stress-induced selective translation. To address this issue, we systematically compared the transcriptome and translatome using self-produced and publicly available datasets to decipher how and to what extent the coordination and discordance between transcription and translation came to be in response to wounding (self-produced), dark to light transition, heat, hypoxia, Pi starvation and the pathogen-associated molecular pattern (elf18) in Arabidopsis. We found that changes in total mRNAs (transcriptome) and ribosome-protected fragments (translatome) are highly correlated upon dark to light transition or heat stress. However, this close correlation was generally lost under other four stresses analyzed in this study, especially during immune response, which suggests that transcription and translation are differentially coordinated under distinct stress conditions. Moreover, Gene Ontology (GO) enrichment analysis showed that typical stress responsive genes were upregulated at both transcriptional and translational levels, while non-stress-specific responsive genes were changed solely at either level or downregulated at both levels. Taking wounding responsive genes for example, typical stress responsive genes are generally involved in functional categories related to dealing with the deleterious effects caused by the imposed wounding stress, such as response to wounding, response to water deprivation and response to jasmonic acid, whereas non-stress-specific responsive genes are often enriched in functional categories like S-glycoside biosynthetic process, photosynthesis and DNA-templated transcription. Collectively, our results revealed the differential as well as targeted coordination between transcriptome and translatome in response to diverse stresses, thus suggesting a potential model wherein preferential ribosome loading onto the stress-upregulated mRNA pool could be a pacing factor for selective translation.

Ribosome collisions trigger cis-acting feedback inhibition of translation initiation.

Translation of aberrant mRNAs can cause ribosomes to stall, leading to collisions with trailing ribosomes. Collided ribosomes are specifically recognised by ZNF598 to initiate protein and mRNA quality control pathways. Here we found using quantitative proteomics of collided ribosomes that EDF1 is a ZNF598-independent sensor of ribosome collisions. EDF1 stabilises GIGYF2 at collisions to inhibit translation initiation in cis via 4EHP. The GIGYF2 axis acts independently of the ZNF598 axis, but each pathway's output is more pronounced without the other. We propose that the widely conserved and highly abundant EDF1 monitors the transcriptome for excessive ribosome density, then triggers a GIGYF2-mediated response to locally and temporarily reduce ribosome loading. Only when collisions persist is translation abandoned to initiate ZNF598-dependent quality control. This tiered response to ribosome collisions would allow cells to dynamically tune translation rates while ensuring fidelity of the resulting protein products.

The translational activator Sov1 coordinates mitochondrial gene expression with mitoribosome biogenesis.

Mitoribosome biogenesis is an expensive metabolic process that is essential to maintain cellular respiratory capacity and requires the stoichiometric accumulation of rRNAs and proteins encoded in two distinct genomes. In yeast, the ribosomal protein Var1, alias uS3m, is mitochondrion-encoded. uS3m is a protein universally present in all ribosomes, where it forms part of the small subunit (SSU) mRNA entry channel and plays a pivotal role in ribosome loading onto the mRNA. However, despite its critical functional role, very little is known concerning VAR1 gene expression. Here, we demonstrate that the protein Sov1 is an in bona fide VAR1 mRNA translational activator and additionally interacts with newly synthesized Var1 polypeptide. Moreover, we show that Sov1 assists the late steps of mtSSU biogenesis involving the incorporation of Var1, an event necessary for uS14 and mS46 assembly. Notably, we have uncovered a translational regulatory mechanism by which Sov1 fine-tunes Var1 synthesis with its assembly into the mitoribosome.

Polysomes Bypass a 50-Nucleotide Coding Gap Less Efficiently Than Monosomes Due to Attenuation of a 5′ mRNA Stem–Loop and Enhanced Drop-off

Efficient translational bypassing of a 50-nt non-coding gap in a phage T4 topoisomerase subunit gene (gp60) requires several recoding signals. Here we investigate the function of the mRNA stem–loop 5′ of the take-off codon, as well as the importance of ribosome loading density on the mRNA for efficient bypassing. We show that polysomes are less efficient at mediating bypassing than monosomes, both in vitro and in vivo, due to their preventing formation of a stem–loop 5′ of the take-off codon and allowing greater peptidyl-tRNA drop off. A ribosome profiling analysis of phage T4-infected Escherichia coli yielded protected mRNA fragments within the normal size range derived from ribosomes stalled at the take-off codon. However, ribosomes at this position also yielded some 53-nucleotide fragments, 16 longer. These were due to protection of the nucleotides that form the 5′ stem–loop. NMR shows that the 5′ stem–loop is highly dynamic. The importance of different nucleotides in the 5′ stem–loop is revealed by mutagenesis studies. These data highlight the significance of the 5′ stem–loop for the 50-nt bypassing and further enhance appreciation of relevance of the extent of ribosome loading for recoding.

Light-Dependent Activation of the GCN2 Kinase Under Cold and Salt Stress Is Mediated by the Photosynthetic Status of the Chloroplast.

Regulation of cytosolic mRNA translation is a key node for rapid adaptation to environmental stress conditions. In yeast and animals, phosphorylation of the α-subunit of eukaryotic translation initiation factor eIF2 is the most thoroughly characterized event for regulating global translation under stress. In plants, the GCN2 kinase (General Control Nonderepressible-2) is the only known kinase for eIF2α. GCN2 is activated under a variety of stresses including reactive oxygen species (ROS). Here, we provide new evidence that the GCN2 kinase in Arabidopsis is also activated rapidly and in a light-dependent manner by cold and salt treatments. These treatments alone did not repress global mRNA ribosome loading in a major way. The activation of GCN2 was accompanied by a more oxidative environment and was attenuated by inhibitors of photosynthetic electron transport, suggesting that it is gated by the redox poise or the reactive oxygen status of the chloroplast. In keeping with these results, gcn2 mutant seedlings were more sensitive than wild type to both cold and salt in a root elongation assay. These data suggest that cold and salt stress may both affect the status of the cytosolic translation apparatus via the conserved GCN2-eIF2α module. The potential role of the GCN2 kinase pathway in the global repression of translation under abiotic stress is discussed.

The RNA-binding protein DAZL functions as repressor and activator of mRNA translation during oocyte maturation

Deleted in azoospermia-like (DAZL) is an RNA-binding protein critical for gamete development. In full-grown oocytes, the DAZL protein increases 4-fold during reentry into the meiotic cell cycle. Here, we have investigated the functional significance of this accumulation at a genome-wide level. Depletion of DAZL causes a block in maturation and widespread disruption in the pattern of ribosome loading on maternal transcripts. In addition to decreased translation, DAZL depletion also causes translational activation of a distinct subset of mRNAs both in quiescent and maturing oocytes, a function recapitulated with YFP-3′UTR reporters. DAZL binds to mRNAs whose translation is both repressed and activated during maturation. Injection of recombinant DAZL protein in DAZL-depleted oocytes rescues the translation and maturation to MII. Mutagenesis of putative DAZL-binding sites in these mRNAs mimics the effect of DAZL depletion. These findings demonstrate that DAZL regulates translation of maternal mRNAs, functioning both as the translational repressor and activator during oocyte maturation.

Small RNAs OmrA and OmrB promote class III flagellar gene expression by inhibiting the synthesis of anti-Sigma factor FlgM

ABSTRACT Bacteria can move by a variety of mechanisms, the best understood being flagella-mediated motility. Flagellar genes are organized in a three-tiered cascade allowing for temporally regulated expression that involves both transcriptional and post-transcriptional control. The class I operon encodes the master regulator FlhDC that drives class II gene transcription. Class II genes include fliA and flgM, which encode the Sigma factor σ28, required for class III transcription, and the anti-Sigma factor FlgM, which inhibits σ28 activity, respectively. The flhDC mRNA is regulated by several small regulatory RNAs (sRNAs). Two of these, the sequence-related OmrA and OmrB RNAs, inhibit FlhD synthesis. Here, we report on a second layer of sRNA-mediated control downstream of FhlDC in the flagella pathway. By mutational analysis, we confirm that a predicted interaction between the conserved 5ʹ seed sequences of OmrA/B and the early coding sequence in flgM mRNA reduces FlgM expression. Regulation is dependent on the global RNA-binding protein Hfq. In vitro experiments support a canonical mechanism: binding of OmrA/B prevents ribosome loading and decreases FlgM protein synthesis. Simultaneous inhibition of both FlhD and FlgM synthesis by OmrA/B complicated an assessment of how regulation of FlgM alone impacts class III gene transcription. Using a combinatorial mutation strategy, we were able to uncouple these two targets and demonstrate that OmrA/B-dependent inhibition of FlgM synthesis liberates σ28 to ultimately promote higher expression of the class III flagellin gene fliC.

FttA is a CPSF73 homologue that terminates transcription in Archaea.

Regulated gene expression is achieved in large part by controlling the activities of essential, multi-subunit RNA polymerase transcription elongation complexes (TECs). The extreme stability required of TECs to processively transcribe large genomic regions necessitates robust mechanisms to terminate transcription. Efficient transcription termination is particularly critical for gene-dense bacterial and archaeal genomes1-3 wherein continued transcription would necessarily transcribe immediately adjacent genes, result in conflicts between the transcription and replication apparatuses4-6 and the coupling of transcription and translation7,8 would permit loading of ribosomes onto aberrant transcripts. Only select sequences or transcription termination factors can disrupt the otherwise extremely stable TEC and we demonstrate that one of the last universally conserved archaeal proteins with unknown biological function is the Factor that terminates transcription in Archaea (FttA). FttA resolves the dichotomy of a prokaryotic gene structure (operons and polarity) and eukaryotic molecular homology (general transcription apparatus) observed in Archaea. This missing-link between prokaryotic and eukaryotic transcription regulation provides the most parsimonious link to the evolution of the processing activities involved in RNA 3’-end formation in Eukarya.

Light Activates the Translational Regulatory Kinase GCN2 via Reactive Oxygen Species Emanating from the Chloroplast.

Cytosolic mRNA translation is subject to global and mRNA-specific controls. Phosphorylation of the translation initiation factor eIF2α anchors a reversible regulatory switch that represses cytosolic translation globally. The stress-responsive GCN2 kinase is the only known kinase for eIF2α serine 56 in Arabidopsis (Arabidopsis thaliana). Here, we show that conditions that generate reactive oxygen species (ROS) in the chloroplast, including dark-light transitions, high light, and the herbicide methyl viologen, rapidly activated GCN2 kinase, whereas mitochondrial and endoplasmic reticulum stress did not. GCN2 activation was light dependent and mitigated by photosynthesis inhibitors and ROS quenchers. Accordingly, the seedling growth of multiple Arabidopsis gcn2 mutants was retarded under excess light conditions, implicating the GCN2-eIF2α pathway in responses to light and associated ROS. Once activated, GCN2 kinase preferentially suppressed the ribosome loading of mRNAs for functions such as mitochondrial ATP synthesis, the chloroplast thylakoids, vesicle trafficking, and translation. The gcn2 mutant overaccumulated transcripts functionally related to abiotic stress, including oxidative stress, as well as innate immune responses. Accordingly, gcn2 displayed defects in immune priming by the fungal elicitor, chitin. Therefore, we provide evidence that reactive oxygen species produced by the photosynthetic apparatus help activate the highly conserved GCN2 kinase, leading to eIF2α phosphorylation and thus affecting the status of the cytosolic protein synthesis apparatus.

MRNA Metabolism in Cardiac Development and Disease: Life After Transcription.

The central dogma of molecular biology illustrates the importance of mRNAs as critical mediators between genetic information encoded at the DNA level and proteomes/metabolomes that determine the diverse functional outcome at the cellular and organ levels. Although the total number of protein-producing (coding) genes in the mammalian genome is ~20,000, it is evident that the intricate processes of cardiac development and the highly regulated physiological regulation in the normal heart, as well as the complex manifestation of pathological remodeling in a diseased heart, would require a much higher degree of complexity at the transcriptome level and beyond. Indeed, in addition to an extensive regulatory scheme implemented at the level of transcription, the complexity of transcript processing following transcription is dramatically increased. RNA processing includes post-transcriptional modification, alternative splicing, editing and transportation, ribosomal loading, and degradation. While transcriptional control of cardiac genes has been a major focus of investigation in recent decades, a great deal of progress has recently been made in our understanding of how post-transcriptional regulation of mRNA contributes to transcriptome complexity. In this review, we highlight some of the key molecular processes and major players in RNA maturation and post-transcriptional regulation. In addition, we provide an update to the recent progress made in the discovery of RNA processing regulators implicated in cardiac development and disease. While post-transcriptional modulation is a complex and challenging problem to study, recent technological advancements are paving the way for a new era of exciting discoveries and potential clinical translation in the context of cardiac biology and heart disease.

Local translation of yeast ERG4 mRNA at the endoplasmic reticulum requires the brefeldin A resistance protein Bfr1.

Brefeldin A resistance factor 1 (Bfr1p) is a nonessential RNA-binding protein and multicopy suppressor of brefeldin A sensitivity in Saccharomyces cerevisiae. Deletion of BFR1 leads to multiple defects, including altered cell shape and size, change in ploidy, induction of P-bodies and chromosomal missegregation. Bfr1p has been shown to associate with polysomes, binds to several hundred mRNAs, and can target some of them to P-bodies. Although this implies a role of Bfr1p in translational control of mRNAs, its molecular function remains elusive. In the present study, we show that mutations in RNA-binding residues of Bfr1p impede its RNA-dependent colocalization with ER, yet do not mimic the known cellular defects seen upon BFR1 deletion. However, a Bfr1 RNA-binding mutant is impaired in binding to ERG4 mRNA, which encodes an enzyme required for the final step of ergosterol biosynthesis. Consistently, bfr1Δ strains show a strong reduction in Erg4p protein levels, most likely because of degradation of misfolded Erg4p. Polysome profiling of bfr1Δ or bfr1 mutant strains reveals a strong shift of ERG4 mRNA to polysomes, consistent with a function of Bfr1p in elongation or increased ribosome loading. Collectively, our data reveal that Bfr1 has at least two separable functions: one in RNA binding and cotranslational protein translocation into the ER and one in ploidy control or chromosome segregation.

Human 5' UTR design and variant effect prediction from a massively parallel translation assay.

The ability to predict the impact of cis-regulatory sequences on gene expression would facilitate discovery in fundamental and applied biology. Here we combine polysome profiling of a library of 280,000 randomized 5' untranslated regions (UTRs) with deep learning to build a predictive model that relates human 5' UTR sequence to translation. Together with a genetic algorithm, we use the model to engineer new 5' UTRs that accurately direct specified levels of ribosome loading, providing the ability to tune sequences for optimal protein expression. We show that the same approach can be extended to chemically modified RNA, an important feature for applications in mRNA therapeutics and synthetic biology. We test 35,212 truncated human 5' UTRs and 3,577 naturally occurring variants and show that the model predicts ribosome loading of these sequences. Finally, we provide evidence of 45 single-nucleotide variants (SNVs) associated with human diseases that substantially change ribosome loading and thus may represent a molecular basis for disease.

Submergence and Waterlogging Stress in Plants: A Review Highlighting Research Opportunities and Understudied Aspects.

Soil flooding creates composite and complex stress in plants known as either submergence or waterlogging stress depending on the depth of the water table. In nature, these stresses are important factors dictating the species composition of the ecosystem. On agricultural land, they cause economic damage associated with long-term social consequences. The understanding of the plant molecular responses to these two stresses has benefited from research studying individual components of the stress, in particular low-oxygen stress. To a lesser extent, other associated stresses and plant responses have been incorporated into the molecular framework, such as ion and ROS signaling, pathogen susceptibility, and organ-specific expression and development. In this review, we aim to highlight known or suspected components of submergence/waterlogging stress that have not yet been thoroughly studied at the molecular level in this context, such as miRNA and retrotransposon expression, the influence of light/dark cycles, protein isoforms, root architecture, sugar sensing and signaling, post-stress molecular events, heavy-metal and salinity stress, and mRNA dynamics (splicing, sequestering, and ribosome loading). Finally, we explore biotechnological strategies that have applied this molecular knowledge to develop cultivars resistant to flooding or to offer alternative uses of flooding-prone soils, like bioethanol and biomass production.