Subcellular spatial transcriptomics identifies three mechanistically different classes of localizing RNAs

Reshaping transcription and translation dynamics during the awakening of the zygotic genome

During the oocyte-to-embryo transition, the transcriptome and proteome are dramatically reshaped. This transition entails a shift from maternally inherited mRNAs to newly synthesized transcripts, produced during the zygotic genome activation (ZGA). Furthermore, a crucial transcription and translation selectivity is required for early embryonic development. Studies across various model organisms have revealed conserved cis- and trans-regulatory mechanisms dictating the regimes by which mRNA and proteins are produced during this critical phase. In this article, we highlight recent technological and conceptual advances that deepen our understanding of how the tuning of both transcription and translation evolves during ZGA.

In-cell structure and snapshots of copia retrotransposons in intact tissue by cryo-ET

Long terminal repeat (LTR) retrotransposons belong to the transposable elements (TEs), autonomously replicating genetic elements that integrate into the host's genome. Among animals, Drosophila melanogaster serves as an important model organism for TE research and contains several LTR retrotransposons, including the Ty1-copia family, which is evolutionarily related to retroviruses and forms virus-like particles (VLPs). In this study, we use cryo-focused ion beam (FIB) milling and lift-out approaches to visualize copia VLPs in ovarian cells and intact egg chambers, resolving the in situ copia capsid structure to 7.7 Å resolution by cryoelectron tomography (cryo-ET). Although cytoplasmic copia VLPs vary in size, nuclear VLPs are homogeneous and form densely packed clusters, supporting a model in which nuclear import acts as a size selector. Analyzing flies deficient in the TE-suppressing PIWI-interacting RNA (piRNA) pathway, we observe copia's translocation into the nucleus during spermatogenesis. Our findings provide insights into the replication cycle and cellular structural biology of an active LTR retrotransposon.

SpaDiT: diffusion transformer for spatial gene expression prediction using scRNA-seq

The rapid development of spatially resolved transcriptomics (SRT) technologies has provided unprecedented opportunities for exploring the structure of specific organs or tissues. However, these techniques (such as image-based SRT) can achieve single-cell resolution, but can only capture the expression levels of tens to hundreds of genes. Such spatial transcriptomics (ST) data, carrying a large number of undetected genes, have limited its application value. To address the challenge, we develop SpaDiT, a deep learning framework for spatial reconstruction and gene expression prediction using scRNA-seq data. SpaDiT employs scRNA-seq data as an a priori condition and utilizes shared genes between ST and scRNA-seq data as latent representations to construct inputs, thereby facilitating the accurate prediction of gene expression in ST data. SpaDiT enhances the accuracy of spatial gene expression predictions over a variety of spatial transcriptomics datasets. We have demonstrated the effectiveness of SpaDiT by conducting extensive experiments on both seq-based and image-based ST data. We compared SpaDiT with eight highly effective baseline methods and found that our proposed method achieved an 8%-12% improvement in performance across multiple metrics. Source code and all datasets used in this paper are available at https://github.com/wenwenmin/SpaDiT and https://zenodo.org/records/12792074.

Centrocortin potentiates co-translational localization of its mRNA to the centrosome via dynein

Centrosomes rely upon proteins within the pericentriolar material to nucleate and organize microtubules. Several mRNAs also reside at centrosomes, although less is known about how and why they accumulate there. We previously showed that local Centrocortin (Cen) mRNA supports centrosome separation, microtubule organization, and viability in Drosophila embryos. Here, using Cen mRNA as a model, we examine mechanisms of centrosomal mRNA localization. We find that while the Cen N'-terminus is sufficient for protein enrichment at centrosomes, multiple domains cooperate to concentrate Cen mRNA at this location. We further identify an N'-terminal motif within Cen that is conserved among dynein cargo adaptor proteins and test its contribution to RNA localization. Our results support a model whereby Cen protein enables the accumulation of its own mRNA to centrosomes through a mechanism requiring active translation, microtubules, and the dynein motor complex. Taken together, our data uncover the basis of translation-dependent localization of a centrosomal RNA required for mitotic integrity.

Orchestrated centers for the production of proteins or "translation factories"

The mechanics of how proteins are generated from mRNA is increasingly well understood. However, much less is known about how protein production is coordinated and orchestrated within the crowded intracellular environment, especially in eukaryotic cells. Recent studies suggest that localized sites exist for the coordinated production of specific proteins. These sites have been termed "translation factories" and roles in protein complex formation, protein localization, inheritance, and translation regulation have been postulated. In this article, we review the evidence supporting the translation of mRNA at these sites, the details of their mechanism of formation, and their likely functional significance. Finally, we consider the key uncertainties regarding these elusive structures in cells. This article is categorized under: Translation Translation > Mechanisms RNA Export and Localization > RNA Localization Translation > Regulation.

Live cell imaging of β-tubulin mRNA reveals spatiotemporal expression dynamics in the filamentous fungus Aspergillus oryzae

In filamentous fungi, microtubules are important for polar growth and morphological maintenance and serve as rails for intracellular trafficking. The molecular mechanisms associated with microtubules have been analyzed. However, little is known about when and where tubulin, a component of microtubules, is biosynthesized in multinuclear and multicellular filamentous fungi. In this study, we visualized microtubules based on the enhanced green fluorescence protein (EGFP)-labeled α-tubulin and β-tubulin mRNA tagged by the EGFP-mediated MS2 system in living yellow Koji mold Aspergillus oryzae cells in order to understand the spatiotemporal production mechanism of tubulin. We found that mRNA of btuA, encoding for β-tubulin, localized at dot-like structures through the apical, middle and basal regions of the hyphal cells. In addition, some btuA mRNA dots showed microtubule-dependent motor protein-like dynamics in the cells. Furthermore, it was found that btuA mRNA dots were decreased in the cytoplasm just before mitosis but increased immediately after mitosis, followed by a gradual decrease. In summary, the localization and abundance of β-tubulin mRNA is spatiotemporally regulated in living A. oryzae hyphal cells.

Inferring Stochastic Rates from Heterogeneous Snapshots of Particle Positions

Many imaging techniques for biological systems-like fixation of cells coupled with fluorescence microscopy-provide sharp spatial resolution in reporting locations of individuals at a single moment in time but also destroy the dynamics they intend to capture. These snapshot observations contain no information about individual trajectories, but still encode information about movement and demographic dynamics, especially when combined with a well-motivated biophysical model. The relationship between spatially evolving populations and single-moment representations of their collective locations is well-established with partial differential equations (PDEs) and their inverse problems. However, experimental data is commonly a set of locations whose number is insufficient to approximate a continuous-in-space PDE solution. Here, motivated by popular subcellular imaging data of gene expression, we embrace the stochastic nature of the data and investigate the mathematical foundations of parametrically inferring demographic rates from snapshots of particles undergoing birth, diffusion, and death in a nuclear or cellular domain. Toward inference, we rigorously derive a connection between individual particle paths and their presentation as a Poisson spatial process. Using this framework, we investigate the properties of the resulting inverse problem and study factors that affect quality of inference. One pervasive feature of this experimental regime is the presence of cell-to-cell heterogeneity. Rather than being a hindrance, we show that cell-to-cell geometric heterogeneity can increase the quality of inference on dynamics for certain parameter regimes. Altogether, the results serve as a basis for more detailed investigations of subcellular spatial patterns of RNA molecules and other stochastically evolving populations that can only be observed for single instants in their time evolution.

Polarity-dependent expression and localization of secretory glucoamylase mRNA in filamentous fungal cells

In multinuclear and multicellular filamentous fungi little is known about how mRNAs encoding secreted enzymes are transcribed and localized spatiotemporally. To better understand this process we analyzed mRNA encoding GlaA, a glucoamylase secreted in large amounts by the industrial filamentous fungus Aspergillus oryzae, by the MS2 system, in which mRNA can be visualized in living cells. We found that glaA mRNA was significantly transcribed and localized near the hyphal tip and septum, which are the sites of protein secretion, in polarity-dependent expression and localization manners. We also revealed that glaA mRNA exhibits long-range dynamics in the vicinity of the endoplasmic reticulum (ER) in a manner that is dependent on the microtubule motor proteins kinesin-1 and kinesin-3, but independent of early endosomes. Moreover, we elucidated that although glaA mRNA localized to stress granules (SGs) and processing bodies (PBs) under high temperature, glaA mRNA was not seen under ER stress, suggesting that there are different regulatory mechanisms of glaA mRNA by SG and PB under high temperature and ER stress. Collectively, this study uncovers a dynamic regulatory mechanism of mRNA encoding a secretory enzyme in filamentous fungi.

Tropomyosin 1-I/C coordinates kinesin-1 and dynein motors during oskar mRNA transport

Dynein and kinesin motors mediate long-range intracellular transport, translocating towards microtubule minus and plus ends, respectively. Cargoes often undergo bidirectional transport by binding to both motors simultaneously. However, it is not known how motor activities are coordinated in such circumstances. In the Drosophila female germline, sequential activities of the dynein-dynactin-BicD-Egalitarian (DDBE) complex and of kinesin-1 deliver oskar messenger RNA from nurse cells to the oocyte, and within the oocyte to the posterior pole. We show through in vitro reconstitution that Tm1-I/C, a tropomyosin-1 isoform, links kinesin-1 in a strongly inhibited state to DDBE-associated oskar mRNA. Nuclear magnetic resonance spectroscopy, small-angle X-ray scattering and structural modeling indicate that Tm1-I/C suppresses kinesin-1 activity by stabilizing its autoinhibited conformation, thus preventing competition with dynein until kinesin-1 is activated in the oocyte. Our work reveals a new strategy for ensuring sequential activity of microtubule motors.

Inferring stochastic rates from heterogeneous snapshots of particle positions

Many imaging techniques for biological systems - like fixation of cells coupled with fluorescence microscopy - provide sharp spatial resolution in reporting locations of individuals at a single moment in time but also destroy the dynamics they intend to capture. These snapshot observations contain no information about individual trajectories, but still encode information about movement and demographic dynamics, especially when combined with a well-motivated biophysical model. The relationship between spatially evolving populations and single-moment representations of their collective locations is well-established with partial differential equations (PDEs) and their inverse problems. However, experimental data is commonly a set of locations whose number is insufficient to approximate a continuous-in-space PDE solution. Here, motivated by popular subcellular imaging data of gene expression, we embrace the stochastic nature of the data and investigate the mathematical foundations of parametrically inferring demographic rates from snapshots of particles undergoing birth, diffusion, and death in a nuclear or cellular domain. Toward inference, we rigorously derive a connection between individual particle paths and their presentation as a Poisson spatial process. Using this framework, we investigate the properties of the resulting inverse problem and study factors that affect quality of inference. One pervasive feature of this experimental regime is the presence of cell-to-cell heterogeneity. Rather than being a hindrance, we show that cell-to-cell geometric heterogeneity can increase the quality of inference on dynamics for certain parameter regimes. Altogether, the results serve as a basis for more detailed investigations of subcellular spatial patterns of RNA molecules and other stochastically evolving populations that can only be observed for single instants in their time evolution.

Dominant negative variants in KIF5B cause osteogenesis imperfecta via down regulation of mTOR signaling

BACKGROUND

Kinesin motor proteins transport intracellular cargo, including mRNA, proteins, and organelles. Pathogenic variants in kinesin-related genes have been implicated in neurodevelopmental disorders and skeletal dysplasias. We identified de novo, heterozygous variants in KIF5B, encoding a kinesin-1 subunit, in four individuals with osteogenesis imperfecta. The variants cluster within the highly conserved kinesin motor domain and are predicted to interfere with nucleotide binding, although the mechanistic consequences on cell signaling and function are unknown.

METHODS

To understand the in vivo genetic mechanism of KIF5B variants, we modeled the p.Thr87Ile variant that was found in two patients in the C. elegans ortholog, unc-116, at the corresponding position (Thr90Ile) by CRISPR/Cas9 editing and performed functional analysis. Next, we studied the cellular and molecular consequences of the recurrent p.Thr87Ile variant by microscopy, RNA and protein analysis in NIH3T3 cells, primary human fibroblasts and bone biopsy.

RESULTS

C. elegans heterozygous for the unc-116 Thr90Ile variant displayed abnormal body length and motility phenotypes that were suppressed by additional copies of the wild type allele, consistent with a dominant negative mechanism. Time-lapse imaging of GFP-tagged mitochondria showed defective mitochondria transport in unc-116 Thr90Ile neurons providing strong evidence for disrupted kinesin motor function. Microscopy studies in human cells showed dilated endoplasmic reticulum, multiple intracellular vacuoles, and abnormal distribution of the Golgi complex, supporting an intracellular trafficking defect. RNA sequencing, proteomic analysis, and bone immunohistochemistry demonstrated down regulation of the mTOR signaling pathway that was partially rescued with leucine supplementation in patient cells.

CONCLUSION

We report dominant negative variants in the KIF5B kinesin motor domain in individuals with osteogenesis imperfecta. This study expands the spectrum of kinesin-related disorders and identifies dysregulated signaling targets for KIF5B in skeletal development.

Spatially resolved transcriptomics: a comprehensive review of their technological advances, applications, and challenges

The ability to explore life kingdoms is largely driven by innovations and breakthroughs in technology, from the invention of the microscope 350 years ago to the recent emergence of single-cell sequencing, by which the scientific community has been able to visualize life at an unprecedented resolution. Most recently, the Spatially Resolved Transcriptomics (SRT) technologies have filled the gap in probing the spatial or even three-dimensional organization of the molecular foundation behind the molecular mysteries of life, including the origin of different cellular populations developed from totipotent cells and human diseases. In this review, we introduce recent progresses and challenges on SRT from the perspectives of technologies and bioinformatic tools, as well as the representative SRT applications. With the currently fast-moving progress of the SRT technologies and promising results from early adopted research projects, we can foresee the bright future of such new tools in understanding life at the most profound analytical level.

Control of mRNA fate by its encoded nascent polypeptide

Cells tightly regulate mRNA processing, localization, and stability to ensure accurate gene expression in diverse cellular states and conditions. Most of these regulatory steps have traditionally been thought to occur before translation by the action of RNA-binding proteins. Several recent discoveries highlight multiple co-translational mechanisms that modulate mRNA translation, localization, processing, and stability. These mechanisms operate by recognition of the nascent protein, which is necessarily coupled to its encoding mRNA during translation. Hence, the distinctive sequence or structure of a particular nascent chain can recruit recognition factors with privileged access to the corresponding mRNA in an otherwise crowded cellular environment. Here, we draw on both well-established and recent examples to provide a conceptual framework for how cells exploit nascent protein recognition to direct mRNA fate. These mechanisms allow cells to dynamically and specifically regulate their transcriptomes in response to changes in cellular states to maintain protein homeostasis.

Mechanism of ribosome-associated mRNA degradation during tubulin autoregulation

Microtubules play crucial roles in cellular architecture, intracellular transport, and mitosis. The availability of free tubulin subunits affects polymerization dynamics and microtubule function. When cells sense excess free tubulin, they trigger degradation of the encoding mRNAs, which requires recognition of the nascent polypeptide by the tubulin-specific ribosome-binding factor TTC5. How TTC5 initiates the decay of tubulin mRNAs is unknown. Here, our biochemical and structural analysis reveals that TTC5 recruits the poorly studied protein SCAPER to the ribosome. SCAPER, in turn, engages the CCR4-NOT deadenylase complex through its CNOT11 subunit to trigger tubulin mRNA decay. SCAPER mutants that cause intellectual disability and retinitis pigmentosa in humans are impaired in CCR4-NOT recruitment, tubulin mRNA degradation, and microtubule-dependent chromosome segregation. Our findings demonstrate how recognition of a nascent polypeptide on the ribosome is physically linked to mRNA decay factors via a relay of protein-protein interactions, providing a paradigm for specificity in cytoplasmic gene regulation.

Analysis of RNA processing directly from spatial transcriptomics data reveals previously unknown regulation

Technical advances have led to an explosion in the amount of biological data available in recent years, especially in the field of RNA sequencing. Specifically, spatial transcriptomics (ST) datasets, which allow each RNA molecule to be mapped to the 2D location it originated from within a tissue, have become readily available. Due to computational challenges, ST data has rarely been used to study RNA processing such as splicing or differential UTR usage. We apply the ReadZS and the SpliZ, methods developed to analyze RNA process in scRNA-seq data, to analyze spatial localization of RNA processing directly from ST data for the first time. Using Moran's I metric for spatial autocorrelation, we identify genes with spatially regulated RNA processing in the mouse brain and kidney, re-discovering known spatial regulation in Myl6 and identifying previously-unknown spatial regulation in genes such as Rps24, Gng13, Slc8a1, Gpm6a, Gpx3, ActB, Rps8, and S100A9. The rich set of discoveries made here from commonly used reference datasets provides a small taste of what can be learned by applying this technique more broadly to the large quantity of Visium data currently being created.

RNA localization mechanisms transcend cell morphology

RNA molecules are localized to specific subcellular regions through interactions between RNA regulatory elements and RNA binding proteins (RBPs). Generally, our knowledge of the mechanistic details behind the localization of a given RNA is restricted to a particular cell type. Here, we show that RNA/RBP interactions that regulate RNA localization in one cell type predictably regulate localization in other cell types with vastly different morphologies. To determine transcriptome-wide RNA spatial distributions across the apicobasal axis of human intestinal epithelial cells, we used our recently developed RNA proximity labeling technique, Halo-seq. We found that mRNAs encoding ribosomal proteins (RP mRNAs) were strongly localized to the basal pole of these cells. Using reporter transcripts and single-molecule RNA FISH, we found that pyrimidine-rich motifs in the 5' UTRs of RP mRNAs were sufficient to drive basal RNA localization. Interestingly, the same motifs were also sufficient to drive RNA localization to the neurites of mouse neuronal cells. In both cell types, the regulatory activity of this motif was dependent on it being in the 5' UTR of the transcript, was abolished upon perturbation of the RNA-binding protein LARP1, and was reduced upon inhibition of kinesin-1. To extend these findings, we compared subcellular RNAseq data from neuronal and epithelial cells. We found that the basal compartment of epithelial cells and the projections of neuronal cells were enriched for highly similar sets of RNAs, indicating that broadly similar mechanisms may be transporting RNAs to these morphologically distinct locations. These findings identify the first RNA element known to regulate RNA localization across the apicobasal axis of epithelial cells, establish LARP1 as an RNA localization regulator, and demonstrate that RNA localization mechanisms cut across cell morphologies.