Differentiated evolutionary rates in alternative exons and the implications for splicing regulation

BACKGROUND: Alternatively spliced exons play an important role in the diversification of gene function in most metazoans and are highly regulated by conserved motifs in exons and introns. Two contradicting properties have been associated to evolutionary conserved alternative exons: higher sequence conservation and higher rate of non-synonymous substitutions, relative to constitutive exons. In order to clarify this issue, we have performed an analysis of the evolution of alternative and constitutive exons, using a large set of protein coding exons conserved between human and mouse and taking into account the conservation of the transcript exonic structure. Further, we have also defined a measure of the variation of the arrangement of exonic splicing enhancers (ESE-conservation score) to study the evolution of splicing regulatory sequences. We have used this measure to correlate the changes in the arrangement of ESEs with the divergence of exon and intron sequences. RESULTS: We find evidence for a relation between the lack of conservation of the exonic structure and the weakening of the sequence evolutionary constraints in alternative and constitutive exons. Exons in transcripts with non-conserved exonic structures have higher synonymous (dS) and non-synonymous (dN) substitution rates than exons in conserved structures. Moreover, alternative exons in transcripts with non-conserved exonic structure are the least constrained in sequence evolution, and at high EST-inclusion levels they are found to be very similar to constitutive exons, whereas alternative exons in transcripts with conserved exonic structure have a dS significantly lower than average at all EST-inclusion levels. We also find higher conservation in the arrangement of ESEs in constitutive exons compared to alternative ones. Additionally, the sequence conservation at flanking introns remains constant for constitutive exons at all ESE-conservation values, but increases for alternative exons at high ESE-conservation values. CONCLUSION: We conclude that most of the differences in dN observed between alternative and constitutive exons can be explained by the conservation of the transcript exonic structure. Low dS values are more characteristic of alternative exons with conserved exonic structure, but not of those with non-conserved exonic structure. Additionally, constitutive exons are characterized by a higher conservation in the arrangement of ESEs, and alternative exons with an ESE-conservation similar to that of constitutive exons are characterized by a conservation of the flanking intron sequences higher than average, indicating the presence of more intronic regulatory signals

RNA Dynamics in Alzheimer’s Disease

Alzheimer's disease (AD) is the most common age-related neurodegenerative disorder that heavily burdens healthcare systems worldwide. There is a significant requirement to understand the still unknown molecular mechanisms underlying AD. Current evidence shows that two of the major features of AD are transcriptome dysregulation and altered function of RNA binding proteins (RBPs), both of which lead to changes in the expression of different RNA species, including microRNAs (miRNAs), circular RNAs (circRNAs), long non-coding RNAs (lncRNAs), and messenger RNAs (mRNAs). In this review, we will conduct a comprehensive overview of how RNA dynamics are altered in AD and how this leads to the differential expression of both short and long RNA species. We will describe how RBP expression and function are altered in AD and how this impacts the expression of different RNA species. Furthermore, we will also show how changes in the abundance of specific RNA species are linked to the pathology of AD

Methanol fixation is the method of choice for droplet-based single-cell transcriptomics of neural cells

The main critical step in single-cell transcriptomics is sample preparation. Several methods have been developed to preserve cells after dissociation to uncouple sample handling from library preparation. Yet, the suitability of these methods depends on the cell types to be processed. In this project, we perform a systematic comparison of preservation methods for droplet-based single-cell RNA-seq on neural and glial cells derived from induced pluripotent stem cells. Our results show that while DMSO provides the highest cell quality in terms of RNA molecules and genes detected per cell, it strongly affects the cellular composition and induces the expression of stress and apoptosis genes. In contrast, methanol fixed samples display a cellular composition similar to fresh samples and provide a good cell quality and little expression biases. Taken together, our results show that methanol fixation is the method of choice for performing droplet-based single-cell transcriptomics experiments on neural cell populations. Methanol fixation appears to be the method of choice for droplet-based scRNA-seq on neural cell populations, based on a broad analysis of how various fixation or preservation methods impact the single-cell transcriptomes of hiPSC-derived neural cells

<i>Drosophila</i> Imp iCLIP identifies an RNA assemblage coordinating F-actin formation

BACKGROUND: Post-transcriptional RNA regulons ensure coordinated expression of monocistronic mRNAs encoding functionally related proteins. In this study, we employ a combination of RIP-seq and short- and long-wave individual-nucleotide resolution crosslinking and immunoprecipitation (iCLIP) technologies in Drosophila cells to identify transcripts associated with cytoplasmic ribonucleoproteins (RNPs) containing the RNA-binding protein Imp. RESULTS: We find extensive binding of Imp to 3′ UTRs of transcripts that are involved in F-actin formation. A common denominator of the RNA–protein interface is the presence of multiple motifs with a central UA-rich element flanked by CA-rich elements. Experiments in single cells and intact flies reveal compromised actin cytoskeletal dynamics associated with low Imp levels. The former shows reduced F-actin formation and the latter exhibits abnormal neuronal patterning. This demonstrates a physiological significance of the defined RNA regulon. CONCLUSIONS: Our data imply that Drosophila Imp RNPs may function as cytoplasmic mRNA assemblages that encode proteins which participate in actin cytoskeletal remodeling. Thus, they may facilitate coordinated protein expression in sub-cytoplasmic locations such as growth cones. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s13059-015-0687-0) contains supplementary material, which is available to authorized users

PAGA: graph abstraction reconciles clustering with trajectory inference through a topology preserving map of single cells.

Single-cell RNA-seq quantifies biological heterogeneity across both discrete cell types and continuous cell transitions. Partition-based graph abstraction (PAGA) provides an interpretable graph-like map of the arising data manifold, based on estimating connectivity of manifold partitions ( https://github.com/theislab/paga ). PAGA maps preserve the global topology of data, allow analyzing data at different resolutions, and result in much higher computational efficiency of the typical exploratory data analysis workflow. We demonstrate the method by inferring structure-rich cell maps with consistent topology across four hematopoietic datasets, adult planaria and the zebrafish embryo and benchmark computational performance on one million neurons.Wellcome Trust, MRC, CRUK, Bloodwise, Swedish Research Council, Helmholtz Association, German Center for Cardiovascular Research, German Research Foundatio

Cell type atlas and lineage tree reconstruction of whole adult animals by single-cell transcriptomics

Flatworms of the species Schmidtea mediterranea are immortal –adult animals contain a large pool of pluripotent stem cells that continuously differentiate to all adult cell types. Therefore, single-cell transcriptome profiling of adult animals should reveal mature and progenitor cells. Here, by combining perturbation experiments, gene expression analysis, a computational method that predicts future cell states from the transcriptional changes, and a novel lineage reconstruction method, we placed all major cell types onto a single lineage tree that connects all cells to a single stem cell compartment. We characterize gene expression changes during differentiation and discover cell types important for regeneration. Our results demonstrate the importance of single-cell transcriptome analysis for mapping and reconstructing fundamental processes of developmental and regenerative biology at unprecedented resolution

Transplanting rejuvenated blood stem cells extends lifespan of aged immunocompromised mice

One goal of regenerative medicine is to rejuvenate tissues and extend lifespan by restoring the function of endogenous aged stem cells. However, evidence that somatic stem cells can be targeted in vivo to extend lifespan is still lacking. Here, we demonstrate that after a short systemic treatment with a specific inhibitor of the small RhoGTPase Cdc42 (CASIN), transplanting aged hematopoietic stem cells (HSCs) from treated mice is sufficient to extend the healthspan and lifespan of aged immunocompromised mice without additional treatment. In detail, we show that systemic CASIN treatment improves strength and endurance of aged mice by increasing the myogenic regenerative potential of aged skeletal muscle stem cells. Further, we show that CASIN modifies niche localization and H4K16ac polarity of HSCs in vivo. Single-cell profiling reveals changes in HSC transcriptome, which underlie enhanced lymphoid and regenerative capacity in serial transplantation assays. Overall, we provide proof-of-concept evidence that a short systemic treatment to decrease Cdc42 activity improves the regenerative capacity of different endogenous aged stem cells in vivo, and that rejuvenated HSCs exert a broad systemic effect sufficient to extend murine health- and lifespan

Heterogeneity and Cancer-Related Features in Lymphangioleiomyomatosis Cells and Tissue

Lymphangioleiomyomatosis (LAM) is a rare, low-grade metastasizing disease characterized by cystic lung destruction. LAM can exhibit extensive heterogeneity at the molecular, cellular, and tissue levels. However, the molecular similarities and differences among LAM cells and tissue, and their connection to cancer features are not fully understood. By integrating complementary gene and protein LAM signatures, and single-cell and bulk tissue transcriptome profiles, we show sources of disease heterogeneity, and how they correspond to cancer molecular portraits. Subsets of LAM diseased cells differ with respect to gene expression profiles related to hormones, metabolism, proliferation, and stemness. Phenotypic diseased cell differences are identified by evaluating lumican (LUM) proteoglycan and YB1 transcription factor expression in LAM lung lesions. The RUNX1 and IRF1 transcription factors are predicted to regulate LAM cell signatures, and both regulators are expressed in LAM lung lesions, with differences between spindle-like and epithelioid LAM cells. The cancer single-cell transcriptome profiles most similar to those of LAM cells include a breast cancer mesenchymal cell model and lines derived from pleural mesotheliomas. Heterogeneity is also found in LAM lung tissue, where it is mainly determined by immune system factors. Variable expression of the multifunctional innate immunity protein LCN2 is linked to disease heterogeneity. This protein is found to be more abundant in blood plasma from LAM patients than from healthy women.This research was partially supported by AELAM (ICO-IDIBELL agreement, to M.A. Pujana), The LAM Foundation Seed Grant 2019, to M.A. Pujana, Carlos III Institute of Health grant PI18/01029, to M.A. Pujana and ICI19/00047 to M. Molina-Molina [co-funded by European Regional Development Fund (ERDF), a way to build Europe], Generalitat de Catalunya SGR grant 2017-449, to M.A. Pujana, the CERCA Program for IDIBELL institutional support, and ZonMW-TopZorg grant 842002003, to C.H.M. van Moorsel. M. Plass was supported by a “Ramón y Cajal” contract of the Spanish Ministry of Science and Innovation (RYC2018-024564-I) and J. Moss was supported by the Intramural Research Program of NIH/NHLBI

DGCR8 HITS-CLIP reveals novel functions for the Microprocessor

The Drosha-DGCR8 complex (Microprocessor) is required for microRNA (miRNA) biogenesis. DGCR8 recognizes the RNA substrate, whereas Drosha functions as the endonuclease. High-throughput sequencing and crosslinking immunoprecipitation (HITS-CLIP) was used to identify RNA targets of DGCR8 in human cells. Unexpectedly, miRNAs were not the most abundant targets. DGCR8-bound RNAs also comprised several hundred mRNAs as well as snoRNAs and long non-coding RNAs. We found that the Microprocessor controls the abundance of several mRNAs as well as of MALAT-1. By contrast, DGCR8-mediated cleavage of snoRNAs is independent of Drosha, suggesting the involvement of DGCR8 in cellular complexes with other endonucleases. Interestingly, binding of DGCR8 to cassette exons, acts as a novel mechanism to regulate the relative abundance of alternatively spliced isoforms. Collectively, these data provide new insights in the complex role of DGCR8 in controlling the fate of several classes of RNAs

DGCR8 acts as an adaptor for the exosome complex to degrade double-stranded structured RNAs

The Microprocessor complex (DGCR8/Drosha) is required for microRNA (miRNA) biogenesis but also binds and regulates the stability of several types of cellular RNAs. Of particular interest, DGCR8 controls the stability of mature small nucleolar RNA (snoRNA) transcripts independently of Drosha, suggesting the existence of alternative DGCR8 complex(es) with other nucleases to process a variety of cellular RNAs. Here, we found that DGCR8 copurifies with subunits of the nuclear exosome, preferentially associating with its hRRP6-containing nucleolar form. Importantly, we demonstrate that DGCR8 is essential for the recruitment of the exosome to snoRNAs and to human telomerase RNA. In addition, we show that the DGCR8/exosome complex controls the stability of the human telomerase RNA component (hTR/TERC). Altogether, these data suggest that DGCR8 acts as an adaptor to recruit the exosome complex to structured RNAs and induce their degradation