Single cell‐ and spatial 'omics revolutionize physiology

Single Cell multi- 'Omics and Spatial Transcriptomics are prominent technological highlights of recent years, and both fields still witness a ceaseless firework of novel approaches for high resolution profiling and for additional omics layers. As all life processes in organs and organisms are based on the functions of their fundamental building blocks, the individual cells and their interactions, these methods are of utmost worth for the study of physiology in health and disease. Recent discoveries on embryonic development, tumor immunology, detailed cellular composition and function of complex tissues like for example the kidney or the brain, different roles of the same cell type in different organs, the oncogenic program of individual tumor entities, or the architecture of immunopathology in infected tissue are based on single cell and spatial transcriptomics experiments. In this review, we will give a broad overview of technological concepts for single cell and spatial analysis, showing both advantages and limitations, and illustrate their impact with some particularly impressive case studies

Discovering a junction tree behind a Markov network by a greedy algorithm

In an earlier paper we introduced a special kind of k-width junction tree, called k-th order t-cherry junction tree in order to approximate a joint probability distribution. The approximation is the best if the Kullback-Leibler divergence between the true joint probability distribution and the approximating one is minimal. Finding the best approximating k-width junction tree is NP-complete if k>2. In our earlier paper we also proved that the best approximating k-width junction tree can be embedded into a k-th order t-cherry junction tree. We introduce a greedy algorithm resulting very good approximations in reasonable computing time. In this paper we prove that if the Markov network underlying fullfills some requirements then our greedy algorithm is able to find the true probability distribution or its best approximation in the family of the k-th order t-cherry tree probability distributions. Our algorithm uses just the k-th order marginal probability distributions as input. We compare the results of the greedy algorithm proposed in this paper with the greedy algorithm proposed by Malvestuto in 1991.Comment: The paper was presented at VOCAL 2010 in Veszprem, Hungar

A protocol for laser microdissection (LMD) followed by transcriptome analysis of plant reproductive tissue in phylogenetically distant angiosperms

BACKGROUND: Plant development is controlled by the action of many, often connected gene regulatory networks. Differential gene expression controlled by internal and external cues is a major driver of growth and time specific differentiation in plants. Transcriptome analysis is the state-of-the-art method to detect spatio-temporal changes in gene expression during development. Monitoring changes in gene expression at early stages or in small plant organs and tissues requires an accurate technique of tissue isolation, which subsequently results in RNA of sufficient quality and quantity. Laser-microdissection enables such accurate dissection and collection of desired tissue from sectioned material at a microscopic level for RNA extraction and subsequent downstream analyses, such as transcriptome, proteome, genome or miRNA. RESULTS: A protocol for laser-microdissection, RNA extraction and RNA-seq was optimized and verified for three distant angiosperm species: Arabidopsis thaliana (Brassicaceae), Oryza sativa (Poaceae) and Eschscholzia californica (Papaveraceae). Previously published protocols were improved in processing speed by reducing the vacuum intensity and incubation time during tissue fixation and incubation time and cryoprotection and by applying adhesive tape. The sample preparation and sectioning of complex and heterogenous flowers produced adequate histological quality and subsequent RNA extraction from micro-dissected gynoecia reliably generated samples of sufficient quality and quantity on all species for RNA-seq. Expression analysis of growth stage specific A. thaliana and O. sativa transcriptomes showed distinct patterns of expression of chromatin remodelers on different time points of gynoecium morphogenesis from the initiation of development to post-meiotic stages. CONCLUSION: Here we describe a protocol for plant tissue preparation, cryoprotection, cryo-sectioning, laser microdissection and RNA sample preparation for Illumina sequencing of complex plant organs from three phyletically distant plant species. We are confident that this approach is widely applicable to other plant species to enable transcriptome analysis with high spatial resolution in non-model plant species. The protocol is rapid, produces high quality sections of complex organs and results in RNA of adequate quality well suited for RNA-seq approaches. We provide detailed description of each stage of sample preparation with the quality and quantity measurements as well as an analysis of generated transcriptomes

RNA modification mapping with JACUSA2

Several high-throughput antibody-free methods for RNA modification detection from sequencing data have been developed. We present JACUSA2 as a versatile software solution and comprehensive analysis framework for RNA modification detection assays that are based on either the Illumina or Nanopore platform. Importantly, JACUSA2 can integrate information from multiple experiments, such as replicates and different conditions, and different library types, such as first- or second-strand cDNA libraries. We demonstrate its utility, showing analysis workflows for N6-methyladenosine (m6A) and pseudouridine (Ψ) detection on Illumina and Nanopore sequencing data sets. Our software and its R helper package are available as open source solutions

UPF3A and UPF3B are redundant and modular activators of nonsense-mediated mRNA decay in human cells

The paralogous human proteins UPF3A and UPF3B are involved in recognizing mRNAs targeted by nonsense-mediated mRNA decay (NMD). While UPF3B has been demonstrated to support NMD, contradicting reports describe UPF3A either as an NMD activator or inhibitor. Here, we present a comprehensive functional analysis of UPF3A and UPF3B in human cells using combinatory experimental approaches. Overexpression or knockout of UPF3A as well as knockout of UPF3B did not detectably change global NMD activity. In contrast, the co-depletion of UPF3A and UPF3B resulted in a marked NMD inhibition and a transcriptome-wide upregulation of NMD substrates, demonstrating a functional redundancy between both NMD factors. Although current models assume that UPF3 bridges NMD-activating exon-junction complexes (EJC) to the NMD factor UPF2, UPF3B exhibited normal NMD activity in rescue experiments when UPF2 or EJC binding was impaired. Further rescue experiments revealed partially redundant functions of UPF3B domains in supporting NMD, involving both UPF2 and EJC interaction sites and the central region of UPF3. Collectively, UPF3A and UPF3B serve as fault-tolerant NMD activators in human cells

Exon junction complex-associated multi-adapter RNPS1 nucleates splicing regulatory complexes to maintain transcriptome surveillance

The exon junction complex (EJC) is an RNA-binding multi-protein complex with critical functions in post-transcriptional gene regulation. It is deposited on the mRNA during splicing and regulates diverse processes including pre-mRNA splicing and nonsense-mediated mRNA decay (NMD) via various interacting proteins. The peripheral EJC-binding protein RNPS1 was reported to serve two insufficiently characterized functions: suppressing mis-splicing of cryptic splice sites and activating NMD in the cytoplasm. The analysis of transcriptome-wide effects of EJC and RNPS1 knockdowns in different human cell lines supports the conclusion that RNPS1 can moderately influence NMD activity, but is not a globally essential NMD factor. However, numerous aberrant splicing events strongly suggest that the main function of RNPS1 is splicing regulation. Rescue analyses revealed that the RRM and C-terminal domain of RNPS1 both contribute partially to regulate RNPS1-dependent splicing events. We defined the RNPS1 core interactome using complementary immunoprecipitations and proximity labeling, which identified interactions with splicing-regulatory factors that are dependent on the C-terminus or the RRM domain of RNPS1. Thus, RNPS1 emerges as a multifunctional splicing regulator that promotes correct and efficient splicing of different vulnerable splicing events via the formation of diverse splicing-promoting complexes

Human UPF3A and UPF3B enable fault-tolerant activation of nonsense-mediated mRNA decay

The paralogous human proteins UPF3A and UPF3B are involved in recognizing mRNAs targeted by nonsense-mediated mRNA decay (NMD). UPF3B has been demonstrated to support NMD, presumably by bridging an exon junction complex (EJC) to the NMD factor UPF2. The role of UPF3A has been described either as a weak NMD activator or an NMD inhibitor. Here, we present a comprehensive functional analysis of UPF3A and UPF3B in human cells using combinatory experimental approaches. Overexpression or knockout of UPF3A as well as knockout of UPF3B did not substantially change global NMD activity. In contrast, the co-depletion of UPF3A and UPF3B resulted in a marked NMD inhibition and a transcriptome-wide upregulation of NMD substrates, demonstrating a functional redundancy between both NMD factors. In rescue experiments, UPF2 or EJC binding-deficient UPF3B largely retained NMD activity. However, combinations of different mutants, including deletion of the middle domain, showed additive or synergistic effects and therefore failed to maintain NMD. Collectively, UPF3A and UPF3B emerge as fault-tolerant, functionally redundant NMD activators in human cells

Linear mitochondrial DNA is rapidly degraded by components of the replication machinery.

Emerging gene therapy approaches that aim to eliminate pathogenic mutations of mitochondrial DNA (mtDNA) rely on efficient degradation of linearized mtDNA, but the enzymatic machinery performing this task is presently unknown. Here, we show that, in cellular models of restriction endonuclease-induced mtDNA double-strand breaks, linear mtDNA is eliminated within hours by exonucleolytic activities. Inactivation of the mitochondrial 5'-3'exonuclease MGME1, elimination of the 3'-5'exonuclease activity of the mitochondrial DNA polymerase POLG by introducing the p.D274A mutation, or knockdown of the mitochondrial DNA helicase TWNK leads to severe impediment of mtDNA degradation. We do not observe similar effects when inactivating other known mitochondrial nucleases (EXOG, APEX2, ENDOG, FEN1, DNA2, MRE11, or RBBP8). Our data suggest that rapid degradation of linearized mtDNA is performed by the same machinery that is responsible for mtDNA replication, thus proposing novel roles for the participating enzymes POLG, TWNK, and MGME1

Full-length spatial transcriptomics reveals the unexplored isoform diversity of the myocardium post-MI

We introduce Single-cell Nanopore Spatial Transcriptomics (scNaST), a software suite to facilitate the analysis of spatial gene expression from second- and third-generation sequencing, allowing to generate a full-length near-single-cell transcriptional landscape of the tissue microenvironment. Taking advantage of the Visium Spatial platform, we adapted a strategy recently developed to assign barcodes to long-read single-cell sequencing data for spatial capture technology. Here, we demonstrate our workflow using four short axis sections of the mouse heart following myocardial infarction. We constructed a de novo transcriptome using long-read data, and successfully assigned 19,794 transcript isoforms in total, including clinically-relevant, but yet uncharacterized modes of transcription, such as intron retention or antisense overlapping transcription. We showed a higher transcriptome complexity in the healthy regions, and identified intron retention as a mode of transcription associated with the infarct area. Our data revealed a clear regional isoform switching among differentially used transcripts for genes involved in cardiac muscle contraction and tissue morphogenesis. Molecular signatures involved in cardiac remodeling integrated with morphological context may support the development of new therapeutics towards the treatment of heart failure and the reduction of cardiac complications