Open-ST: High-resolution spatial transcriptomics in 3D.

Spatial transcriptomics (ST) methods unlock molecular mechanisms underlying tissue development, homeostasis, or disease. However, there is a need for easy-to-use, high-resolution, cost-efficient, and 3D-scalable methods. Here, we report Open-ST, a sequencing-based, open-source experimental and computational resource to address these challenges and to study the molecular organization of tissues in 2D and 3D. In mouse brain, Open-ST captured transcripts at subcellular resolution and reconstructed cell types. In primary head-and-neck tumors and patient-matched healthy/metastatic lymph nodes, Open-ST captured the diversity of immune, stromal, and tumor populations in space, validated by imaging-based ST. Distinct cell states were organized around cell-cell communication hotspots in the tumor but not the metastasis. Strikingly, the 3D reconstruction and multimodal analysis of the metastatic lymph node revealed spatially contiguous structures not visible in 2D and potential biomarkers precisely at the 3D tumor/lymph node boundary. All protocols and software are available at https://rajewsky-lab.github.io/openst.

Sites of transcription initiation drive mRNA isoform selection.

The generation of distinct messenger RNA isoforms through alternative RNA processing modulates the expression and function of genes, often in a cell-type-specific manner. Here, we assess the regulatory relationships between transcription initiation, alternative splicing, and 3' end site selection. Applying long-read sequencing to accurately represent even the longest transcripts from end to end, we quantify mRNA isoforms in Drosophila tissues, including the transcriptionally complex nervous system. We find that in Drosophila heads, as well as in human cerebral organoids, 3' end site choice is globally influenced by the site of transcription initiation (TSS). "Dominant promoters," characterized by specific epigenetic signatures including p300/CBP binding, impose a transcriptional constraint to define splice and polyadenylation variants. In vivo deletion or overexpression of dominant promoters as well as p300/CBP loss disrupted the 3' end expression landscape. Our study demonstrates the crucial impact of TSS choice on the regulation of transcript diversity and tissue identity.

A variety of dicer substrates in human and C. elegans.

The endoribonuclease Dicer is known for its central role in the biogenesis of eukaryotic small RNAs/microRNAs. Despite its importance, Dicer target transcripts have not been directly mapped. Here, we apply biochemical methods to human cells and C. elegans and identify thousands of Dicer-binding sites. We find known and hundreds of additional miRNAs with high sensitivity and specificity. We also report structural RNAs, promoter RNAs, and mitochondrial transcripts as Dicer targets. Interestingly, most Dicer-binding sites reside on mRNAs/lncRNAs and are not significantly processed into small RNAs. These passive sites typically harbor small, Dicer-bound hairpins within intact transcripts and generally stabilize target expression. We show that passive sites can sequester Dicer and reduce microRNA expression. mRNAs with passive sites were in human and worm significantly associated with processing-body/granule function. Together, we provide the first transcriptome-wide map of Dicer targets and suggest conserved binding modes and functions outside of the miRNA pathway.

Paternal diet defines offspring chromatin state and intergenerational obesity.

The global rise in obesity has revitalized a search for genetic and epigenetic factors underlying the disease. We present a Drosophila model of paternal-diet-induced intergenerational metabolic reprogramming (IGMR) and identify genes required for its encoding in offspring. Intriguingly, we find that as little as 2 days of dietary intervention in fathers elicits obesity in offspring. Paternal sugar acts as a physiological suppressor of variegation, desilencing chromatin-state-defined domains in both mature sperm and in offspring embryos. We identify requirements for H3K9/K27me3-dependent reprogramming of metabolic genes in two distinct germline and zygotic windows. Critically, we find evidence that a similar system may regulate obesity susceptibility and phenotype variation in mice and humans. The findings provide insight into the mechanisms underlying intergenerational metabolic reprogramming and carry profound implications for our understanding of phenotypic variation and evolution.

Context-specific regulation of cell survival by a miRNA-controlled BIM rheostat.

Knockout of the ubiquitously expressed miRNA-17∼92 cluster in mice produces a lethal developmental lung defect, skeletal abnormalities, and blocked B lymphopoiesis. A shared target of miR-17∼92 miRNAs is the pro-apoptotic protein BIM, central to life-death decisions in mammalian cells. To clarify the contribution of miR-17∼92:Bim interactions to the complex miR-17∼92 knockout phenotype, we used a system of conditional mutagenesis of the nine Bim 3' UTR miR-17∼92 seed matches. Blocking miR-17∼92:Bim interactions early in development phenocopied the lethal lung phenotype of miR-17∼92 ablation and generated a skeletal kinky tail. In the hematopoietic system, instead of causing the predicted B cell developmental block, it produced a selective inability of B cells to resist cellular stress; and prevented B and T cell hyperplasia caused by Bim haploinsufficiency. Thus, the interaction of miR-17∼92 with a single target is essential for life, and BIM regulation by miRNAs serves as a rheostat controlling cell survival in specific physiological contexts.

Spatiotemporal, optogenetic control of gene expression in organoids.

Organoids derived from stem cells have become an increasingly important tool for studying human development and modeling disease. However, methods are still needed to control and study spatiotemporal patterns of gene expression in organoids. Here we combined optogenetics and gene perturbation technologies to activate or knock-down RNA of target genes in programmable spatiotemporal patterns. To illustrate the usefulness of our approach, we locally activated Sonic Hedgehog (SHH) signaling in an organoid model for human neurodevelopment. Spatial and single-cell transcriptomic analyses showed that this local induction was sufficient to generate stereotypically patterned organoids and revealed new insights into SHH's contribution to gene regulation in neurodevelopment. With this study, we propose optogenetic perturbations in combination with spatial transcriptomics as a powerful technology to reprogram and study cell fates and tissue patterning in organoids.

miR-7 controls glutamatergic transmission and neuronal connectivity in a Cdr1as-dependent manner.

The circular RNA (circRNA) Cdr1as is conserved across mammals and highly expressed in neurons, where it directly interacts with microRNA miR-7. However, the biological function of this interaction is unknown. Here, using primary cortical murine neurons, we demonstrate that stimulating neurons by sustained depolarization rapidly induces two-fold transcriptional upregulation of Cdr1as and strong post-transcriptional stabilization of miR-7. Cdr1as loss causes doubling of glutamate release from stimulated synapses and increased frequency and duration of local neuronal bursts. Moreover, the periodicity of neuronal networks increases, and synchronicity is impaired. Strikingly, these effects are reverted by sustained expression of miR-7, which also clears Cdr1as molecules from neuronal projections. Consistently, without Cdr1as, transcriptomic changes caused by miR-7 overexpression are stronger (including miR-7-targets downregulation) and enriched in secretion/synaptic plasticity pathways. Altogether, our results suggest that in cortical neurons Cdr1as buffers miR-7 activity to control glutamatergic excitatory transmission and neuronal connectivity important for long-lasting synaptic adaptations.

LifeTime and improving European healthcare through cell-based interceptive medicine.

Here we describe the LifeTime Initiative, which aims to track, understand and target human cells during the onset and progression of complex diseases, and to analyse their response to therapy at single-cell resolution. This mission will be implemented through the development, integration and application of single-cell multi-omics and imaging, artificial intelligence and patient-derived experimental disease models during the progression from health to disease. The analysis of large molecular and clinical datasets will identify molecular mechanisms, create predictive computational models of disease progression, and reveal new drug targets and therapies. The timely detection and interception of disease embedded in an ethical and patient-centred vision will be achieved through interactions across academia, hospitals, patient associations, health data management systems and industry. The application of this strategy to key medical challenges in cancer, neurological and neuropsychiatric disorders, and infectious, chronic inflammatory and cardiovascular diseases at the single-cell level will usher in cell-based interceptive medicine in Europe over the next decade.

Best practice standards for circular RNA research.

Circular RNAs (circRNAs) are formed in all domains of life and via different mechanisms. There has been an explosion in the number of circRNA papers in recent years; however, as a relatively young field, circRNA biology has an urgent need for common experimental standards for isolating, analyzing, expressing and depleting circRNAs. Here we propose a set of guidelines for circRNA studies based on the authors' experience. This Perspective will specifically address the major class of circRNAs in Eukarya that are generated by a spliceosome-catalyzed back-splicing event. We hope that the implementation of best practice principles for circRNA research will help move the field forward and allow a better functional understanding of this fascinating group of RNAs.

Epigenomic Profiling of Human CD4+ T Cells Supports a Linear Differentiation Model and Highlights Molecular Regulators of Memory Development.

The impact of epigenetics on the differentiation of memory T (Tmem) cells is poorly defined. We generated deep epigenomes comprising genome-wide profiles of DNA methylation, histone modifications, DNA accessibility, and coding and non-coding RNA expression in naive, central-, effector-, and terminally differentiated CD45RA+ CD4+ Tmem cells from blood and CD69+ Tmem cells from bone marrow (BM-Tmem). We observed a progressive and proliferation-associated global loss of DNA methylation in heterochromatic parts of the genome during Tmem cell differentiation. Furthermore, distinct gradually changing signatures in the epigenome and the transcriptome supported a linear model of memory development in circulating T cells, while tissue-resident BM-Tmem branched off with a unique epigenetic profile. Integrative analyses identified candidate master regulators of Tmem cell differentiation, including the transcription factor FOXP1. This study highlights the importance of epigenomic changes for Tmem cell biology and demonstrates the value of epigenetic data for the identification of lineage regulators.

Gene expression cartography.

Multiplexed RNA sequencing in individual cells is transforming basic and clinical life sciences1-4. Often, however, tissues must first be dissociated, and crucial information about spatial relationships and communication between cells is thus lost. Existing approaches to reconstruct tissues assign spatial positions to each cell, independently of other cells, by using spatial patterns of expression of marker genes5,6-which often do not exist. Here we reconstruct spatial positions with little or no prior knowledge, by searching for spatial arrangements of sequenced cells in which nearby cells have transcriptional profiles that are often (but not always) more similar than cells that are farther apart. We formulate this task as a generalized optimal-transport problem for probabilistic embedding and derive an efficient iterative algorithm to solve it. We reconstruct the spatial expression of genes in mammalian liver and intestinal epithelium, fly and zebrafish embryos, sections from the mammalian cerebellum and whole kidney, and use the reconstructed tissues to identify genes that are spatially informative. Thus, we identify an organization principle for the spatial expression of genes in animal tissues, which can be exploited to infer meaningful probabilities of spatial position for individual cells. Our framework ('novoSpaRc') can incorporate prior spatial information and is compatible with any single-cell technology. Additional principles that underlie the cartography of gene expression can be tested using our approach.

Translation of CircRNAs.

Circular RNAs (circRNAs) are abundant and evolutionarily conserved RNAs of largely unknown function. Here, we show that a subset of circRNAs is translated in vivo. By performing ribosome footprinting from fly heads, we demonstrate that a group of circRNAs is associated with translating ribosomes. Many of these ribo-circRNAs use the start codon of the hosting mRNA, are bound by membrane-associated ribosomes, and have evolutionarily conserved termination codons. In addition, we found that a circRNA generated from the muscleblind locus encodes a protein, which we detected in fly head extracts by mass spectrometry. Next, by performing in vivo and in vitro translation assays, we show that UTRs of ribo-circRNAs (cUTRs) allow cap-independent translation. Moreover, we found that starvation and FOXO likely regulate the translation of a circMbl isoform. Altogether, our study provides strong evidence for translation of circRNAs, revealing the existence of an unexplored layer of gene activity.

Circ-ZNF609 Is a Circular RNA that Can Be Translated and Functions in Myogenesis.

Circular RNAs (circRNAs) constitute a family of transcripts with unique structures and still largely unknown functions. Their biogenesis, which proceeds via a back-splicing reaction, is fairly well characterized, whereas their role in the modulation of physiologically relevant processes is still unclear. Here we performed expression profiling of circRNAs during in vitro differentiation of murine and human myoblasts, and we identified conserved species regulated in myogenesis and altered in Duchenne muscular dystrophy. A high-content functional genomic screen allowed the study of their functional role in muscle differentiation. One of them, circ-ZNF609, resulted in specifically controlling myoblast proliferation. Circ-ZNF609 contains an open reading frame spanning from the start codon, in common with the linear transcript, and terminating at an in-frame STOP codon, created upon circularization. Circ-ZNF609 is associated with heavy polysomes, and it is translated into a protein in a splicing-dependent and cap-independent manner, providing an example of a protein-coding circRNA in eukaryotes.

SLAM-Drop-seq reveals mRNA kinetic rates throughout the cell cycle.

RNA abundance is tightly regulated in eukaryotic cells by modulating the kinetic rates of RNA production, processing, and degradation. To date, little is known about time­dependent kinetic rates during dynamic processes. Here, we present SLAM­Drop­seq, a method that combines RNA metabolic labeling and alkylation of modified nucleotides in methanol­fixed cells with droplet­based sequencing to detect newly synthesized and preexisting mRNAs in single cells. As a first application, we sequenced 7280 HEK293 cells and calculated gene­specific kinetic rates during the cell cycle using the novel package Eskrate. Of the 377 robust­cycling genes that we identified, only a minor fraction is regulated solely by either dynamic transcription or degradation (6 and 4%, respectively). By contrast, the vast majority (89%) exhibit dynamically regulated transcription and degradation rates during the cell cycle. Our study thus shows that temporally regulated mRNA degradation is fundamental for the correct expression of a majority of cycling genes. SLAM­Drop­seq, combined with Eskrate, is a powerful approach to understanding the underlying mRNA kinetics of single­cell gene expression dynamics in continuous biological processes.

The Lupus Autoantigen La Prevents Mis-channeling of tRNA Fragments into the Human MicroRNA Pathway.

The Lupus autoantigen La is an RNA-binding protein that stabilizes RNA polymerase III (Pol III) transcripts and supports RNA folding and has in addition been implicated in the mammalian microRNA (miRNA) pathway. Here, we have analyzed effects of La depletion on Argonaute (Ago)-bound small RNAs in human cells. We find that in the absence of La, distinct tRNA fragments are loaded into Ago proteins. Thus, La functions as gatekeeper ensuring correct tRNA maturation and protecting the miRNA pathway from potentially functional tRNA fragments. However, one specific isoleucin pre-tRNA produces both a functional tRNA and a miRNA even when La is present. We demonstrate that the fully complementary 5' leader and 3' trailer of the pre-tRNA-Ile form a double-stranded RNA molecule that has low affinity to La. Instead, Exportin-5 (Xpo5) recognizes it as miRNA precursor and transports it into the cytoplasm for Dicer processing and Ago loading.

Dicer ablation affects antibody diversity and cell survival in the B lymphocyte lineage.

To explore the role of Dicer-dependent control mechanisms in B lymphocyte development, we ablated this enzyme in early B cell progenitors. This resulted in a developmental block at the pro- to pre-B cell transition. Gene-expression profiling revealed a miR-17 approximately 92 signature in the 3'UTRs of genes upregulated in Dicer-deficient pro-B cells; a top miR-17 approximately 92 target, the proapoptotic molecule Bim, was highly upregulated. Accordingly, B cell development could be partially rescued by ablation of Bim or transgenic expression of the prosurvival protein Bcl-2. This allowed us to assess the impact of Dicer deficiency on the V(D)J recombination program in developing B cells. We found intact Ig gene rearrangements in immunoglobulin heavy (IgH) and kappa chain loci, but increased sterile transcription and usage of D(H) elements of the DSP family in IgH, and increased N sequence addition in Igkappa due to deregulated transcription of the terminal deoxynucleotidyl transferase gene.

Circular RNAs in the Mammalian Brain Are Highly Abundant, Conserved, and Dynamically Expressed.

Circular RNAs (circRNAs) are an endogenous class of animal RNAs. Despite their abundance, their function and expression in the nervous system are unknown. Therefore, we sequenced RNA from different brain regions, primary neurons, isolated synapses, as well as during neuronal differentiation. Using these and other available data, we discovered and analyzed thousands of neuronal human and mouse circRNAs. circRNAs were extraordinarily enriched in the mammalian brain, well conserved in sequence, often expressed as circRNAs in both human and mouse, and sometimes even detected in Drosophila brains. circRNAs were overall upregulated during neuronal differentiation, highly enriched in synapses, and often differentially expressed compared to their mRNA isoforms. circRNA expression correlated negatively with expression of the RNA-editing enzyme ADAR1. Knockdown of ADAR1 induced elevated circRNA expression. Together, we provide a circRNA brain expression atlas and evidence for important circRNA functions and values as biomarkers.

Long-term imaging reveals behavioral plasticity during C. elegans dauer exit.

BACKGROUND: During their lifetime, animals must adapt their behavior to survive in changing environments. This ability requires the nervous system to undergo adjustments at distinct temporal scales, from short-term dynamic changes in expression of neurotransmitters and receptors to longer-term growth, spatial and connectivity reorganization, while integrating external stimuli. The nematode Caenorhabditis elegans provides a model of nervous system plasticity, in particular its dauer exit decision. Under unfavorable conditions, larvae will enter the non-feeding and non-reproductive stress-resistant dauer stage and adapt their behavior to cope with the harsh new environment, with active reversal under improved conditions leading to resumption of reproductive development. However, how different environmental stimuli regulate the exit decision mechanism and thereby drive the larva's behavioral change is unknown. To fill this gap and provide insights on behavioral changes over extended periods of time, we developed a new open hardware method for long-term imaging (12h) of C. elegans larvae. RESULTS: Our WormObserver platform comprises open hardware and software components for video acquisition, automated processing of large image data (> 80k images/experiment) and data analysis. We identified dauer-specific behavioral motifs and characterized the behavioral trajectory of dauer exit in different environments and genetic backgrounds to identify key decision points and stimuli promoting dauer exit. Combining long-term behavioral imaging with transcriptomics data, we find that bacterial ingestion triggers a change in neuropeptide gene expression to establish post-dauer behavior. CONCLUSIONS: Taken together, we show how a developing nervous system can robustly integrate environmental changes activate a developmental switch and adapt the organism's behavior to a new environment. WormObserver is generally applicable to other research questions within and beyond the C. elegans field, having a modular and customizable character and allowing assessment of behavioral plasticity over longer periods.

Urinary single-cell sequencing captures kidney injury and repair processes in human acute kidney injury.

Acute kidney injury (AKI) is a major health issue, the outcome of which depends primarily on damage and reparative processes of tubular epithelial cells. Mechanisms underlying AKI remain incompletely understood, specific therapies are lacking and monitoring the course of AKI in clinical routine is confined to measuring urine output and plasma levels of filtration markers. Here we demonstrate feasibility and potential of a novel approach to assess the cellular and molecular dynamics of AKI by establishing a robust urine-to-single cell RNA sequencing (scRNAseq) pipeline for excreted kidney cells via flow cytometry sorting. We analyzed 42,608 single cell transcriptomes of 40 urine samples from 32 patients with AKI and compared our data with reference material from human AKI post-mortem biopsies and published mouse data. We demonstrate that tubular epithelial cells transcriptomes mirror kidney pathology and reflect distinct injury and repair processes, including oxidative stress, inflammation, and tissue rearrangement. We also describe an AKI-specific abundant urinary excretion of adaptive progenitor-like cells. Thus, single cell transcriptomics of kidney cells excreted in urine provides noninvasive, unprecedented insight into cellular processes underlying AKI, thereby opening novel opportunities for target identification, AKI sub-categorization, and monitoring of natural disease course and interventions.