A distinct class of pan-cancer susceptibility genes revealed by an alternative polyadenylation transcriptome-wide association study.

Alternative polyadenylation plays an important role in cancer initiation and progression; however, current transcriptome-wide association studies mostly ignore alternative polyadenylation when identifying putative cancer susceptibility genes. Here, we perform a pan-cancer 3' untranslated region alternative polyadenylation transcriptome-wide association analysis by integrating 55 well-powered (n > 50,000) genome-wide association studies datasets across 22 major cancer types with alternative polyadenylation quantification from 23,955 RNA sequencing samples across 7,574 individuals. We find that genetic variants associated with alternative polyadenylation are co-localized with 28.57% of cancer loci and contribute a significant portion of cancer heritability. We further identify 642 significant cancer susceptibility genes predicted to modulate cancer risk via alternative polyadenylation, 62.46% of which have been overlooked by traditional expression- and splicing- studies. As proof of principle validation, we show that alternative alleles facilitate 3' untranslated region lengthening of CRLS1 gene leading to increased protein abundance and promoted proliferation of breast cancer cells. Together, our study highlights the significant role of alternative polyadenylation in discovering new cancer susceptibility genes and provides a strong foundational framework for enhancing our understanding of the etiology underlying human cancers.

scQTLbase: an integrated human single-cell eQTL database.

Genome-wide association studies (GWAS) have identified numerous genetic variants associated with diseases and traits. However, the functional interpretation of these variants remains challenging. Expression quantitative trait loci (eQTLs) have been widely used to identify mutations linked to disease, yet they explain only 20-50% of disease-related variants. Single-cell eQTLs (sc-eQTLs) studies provide an immense opportunity to identify new disease risk genes with expanded eQTL scales and transcriptional regulation at a much finer resolution. However, there is no comprehensive database dedicated to single-cell eQTLs that users can use to search, analyse and visualize them. Therefore, we developed the scQTLbase (http://bioinfo.szbl.ac.cn/scQTLbase), the first integrated human sc-eQTLs portal, featuring 304 datasets spanning 57 cell types and 95 cell states. It contains ∼16 million SNPs significantly associated with cell-type/state gene expression and ∼0.69 million disease-associated sc-eQTLs from 3 333 traits/diseases. In addition, scQTLbase offers sc-eQTL search, gene expression visualization in UMAP plots, a genome browser, and colocalization visualization based on the GWAS dataset of interest. scQTLbase provides a one-stop portal for sc-eQTLs that will significantly advance the discovery of disease susceptibility genes.

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.

Neural cell diversity in the light of single-cell transcriptomics.

The development of highly parallel and affordable high-throughput single-cell transcriptomics technologies has revolutionized our understanding of brain complexity. These methods have been used to build cellular maps of the brain, its different regions, and catalog the diversity of cells in each of them during development, aging and even in disease. Now we know that cellular diversity is way beyond what was previously thought. Single-cell transcriptomics analyses have revealed that cell types previously considered homogeneous based on imaging techniques differ depending on several factors including sex, age and location within the brain. The expression profiles of these cells have also been exploited to understand which are the regulatory programs behind cellular diversity and decipher the transcriptional pathways driving them. In this review, we summarize how single-cell transcriptomics have changed our view on the cellular diversity in the human brain, and how it could impact the way we study neurodegenerative diseases. Moreover, we describe the new computational approaches that can be used to study cellular differentiation and gain insight into the functions of individual cell populations under different conditions and their alterations in disease.

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.

The X-linked splicing regulator MBNL3 has been co-opted to restrict placental growth in eutherians.

Understanding the regulatory interactions that control gene expression during the development of novel tissues is a key goal of evolutionary developmental biology. Here, we show that Mbnl3 has undergone a striking process of evolutionary specialization in eutherian mammals resulting in the emergence of a novel placental function for the gene. Mbnl3 belongs to a family of RNA-binding proteins whose members regulate multiple aspects of RNA metabolism. We find that, in eutherians, while both Mbnl3 and its paralog Mbnl2 are strongly expressed in placenta, Mbnl3 expression has been lost from nonplacental tissues in association with the evolution of a novel promoter. Moreover, Mbnl3 has undergone accelerated protein sequence evolution leading to changes in its RNA-binding specificities and cellular localization. While Mbnl2 and Mbnl3 share partially redundant roles in regulating alternative splicing, polyadenylation site usage and, in turn, placenta maturation, Mbnl3 has also acquired novel biological functions. Specifically, Mbnl3 knockout (M3KO) alone results in increased placental growth associated with higher Myc expression. Furthermore, Mbnl3 loss increases fetal resource allocation during limiting conditions, suggesting that location of Mbnl3 on the X chromosome has led to its role in limiting placental growth, favoring the maternal side of the parental genetic conflict.

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.

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. IMPLICATIONS: This study identifies LAM molecular and cellular features, master regulators, cancer similarities, and potential causes of disease heterogeneity.

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.

Cell type atlas and lineage tree of a whole complex animal 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 into all adult cell types. Therefore, single-cell transcriptome profiling of adult animals should reveal mature and progenitor cells. By combining perturbation experiments, gene expression analysis, a computational method that predicts future cell states from transcriptional changes, and a 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 characterized gene expression changes during differentiation and discovered 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 high resolution.

Loss of a mammalian circular RNA locus causes miRNA deregulation and affects brain function.

Hundreds of circular RNAs (circRNAs) are highly abundant in the mammalian brain, often with conserved expression. Here we show that the circRNA Cdr1as is massively bound by the microRNAs (miRNAs) miR-7 and miR-671 in human and mouse brains. When the Cdr1as locus was removed from the mouse genome, knockout animals displayed impaired sensorimotor gating-a deficit in the ability to filter out unnecessary information-which is associated with neuropsychiatric disorders. Electrophysiological recordings revealed dysfunctional synaptic transmission. Expression of miR-7 and miR-671 was specifically and posttranscriptionally misregulated in all brain regions analyzed. Expression of immediate early genes such as Fos, a direct miR-7 target, was enhanced in Cdr1as-deficient brains, providing a possible molecular link to the behavioral phenotype. Our data indicate an in vivo loss-of-function circRNA phenotype and suggest that interactions between Cdr1as and miRNAs are important for normal brain function.

Highly accessible AU-rich regions in 3' untranslated regions are hotspots for binding of regulatory factors.

Post-transcriptional regulation is regarded as one of the major processes involved in the regulation of gene expression. It is mainly performed by RNA binding proteins and microRNAs, which target RNAs and typically affect their stability. Recent efforts from the scientific community have aimed at understanding post-transcriptional regulation at a global scale by using high-throughput sequencing techniques such as cross-linking and immunoprecipitation (CLIP), which facilitates identification of binding sites of these regulatory factors. However, the diversity in the experimental procedures and bioinformatics analyses has hindered the integration of multiple datasets and thus limited the development of an integrated view of post-transcriptional regulation. In this work, we have performed a comprehensive analysis of 107 CLIP datasets from 49 different RBPs in HEK293 cells to shed light on the complex interactions that govern post-transcriptional regulation. By developing a more stringent CLIP analysis pipeline we have discovered the existence of conserved regulatory AU-rich regions in the 3'UTRs where miRNAs and RBPs that regulate several processes such as polyadenylation or mRNA stability bind. Analogous to promoters, many factors have binding sites overlapping or in close proximity in these hotspots and hence the regulation of the mRNA may depend on their relative concentrations. This hypothesis is supported by RBP knockdown experiments that alter the relative concentration of RBPs in the cell. Upon AGO2 knockdown (KD), transcripts containing "free" target sites show increased expression levels compared to those containing target sites in hotspots, which suggests that target sites within hotspots are less available for miRNAs to bind. Interestingly, these hotspots appear enriched in genes with regulatory functions such as DNA binding and RNA binding. Taken together, our results suggest that hotspots are functional regulatory elements that define an extra layer of regulation of post-transcriptional regulatory networks.

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.

Drosophila 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.

Alternative polyadenylation of tumor suppressor genes in small intestinal neuroendocrine tumors.

The tumorigenesis of small intestinal neuroendocrine tumors (SI-NETs) is poorly understood. Recent studies have associated alternative polyadenylation (APA) with proliferation, cell transformation, and cancer. Polyadenylation is the process in which the pre-messenger RNA is cleaved at a polyA site and a polyA tail is added. Genes with two or more polyA sites can undergo APA. This produces two or more distinct mRNA isoforms with different 3' untranslated regions. Additionally, APA can also produce mRNAs containing different 3'-terminal coding regions. Therefore, APA alters both the repertoire and the expression level of proteins. Here, we used high-throughput sequencing data to map polyA sites and characterize polyadenylation genome-wide in three SI-NETs and a reference sample. In the tumors, 16 genes showed significant changes of APA pattern, which lead to either the 3' truncation of mRNA coding regions or 3' untranslated regions. Among these, 11 genes had been previously associated with cancer, with 4 genes being known tumor suppressors: DCC, PDZD2, MAGI1, and DACT2. We validated the APA in three out of three cases with quantitative real-time-PCR. Our findings suggest that changes of APA pattern in these 16 genes could be involved in the tumorigenesis of SI-NETs. Furthermore, they also point to APA as a new target for both diagnostic and treatment of SI-NETs. The identified genes with APA specific to the SI-NETs could be further tested as diagnostic markers and drug targets for disease prevention and treatment.

Approaches to link RNA secondary structures with splicing regulation.

In higher eukaryotes, alternative splicing is usually regulated by protein factors, which bind to the pre-mRNA and affect the recognition of splicing signals. There is recent evidence that the secondary structure of the pre-mRNA may also play an important role in this process, either by facilitating or hindering the interaction with factors and small nuclear ribonucleoproteins (snRNPs) that regulate splicing. Moreover, the secondary structure could play a fundamental role in the splicing of yeast species, which lack many of the regulatory splicing factors present in metazoans. This chapter describes the steps in the analysis of the secondary structure of the pre-mRNA and its possible relation to splicing. As a working example, we use the case of yeast and the problem of the recognition of the 3' splice site (3'ss).

Drosha regulates gene expression independently of RNA cleavage function.

Drosha is the main RNase III-like enzyme involved in the process of microRNA (miRNA) biogenesis in the nucleus. Using whole-genome ChIP-on-chip analysis, we demonstrate that, in addition to miRNA sequences, Drosha specifically binds promoter-proximal regions of many human genes in a transcription-dependent manner. This binding is not associated with miRNA production or RNA cleavage. Drosha knockdown in HeLa cells downregulated nascent gene transcription, resulting in a reduction of polyadenylated mRNA produced from these gene regions. Furthermore, we show that this function of Drosha is dependent on its N-terminal protein-interaction domain, which associates with the RNA-binding protein CBP80 and RNA Polymerase II. Consequently, we uncover a previously unsuspected RNA cleavage-independent function of Drosha in the regulation of human gene expression.

The Microprocessor controls the activity of mammalian retrotransposons.

More than half of the human genome is made of transposable elements whose ongoing mobilization is a driving force in genetic diversity; however, little is known about how the host regulates their activity. Here, we show that the Microprocessor (Drosha-DGCR8), which is required for microRNA biogenesis, also recognizes and binds RNAs derived from human long interspersed element 1 (LINE-1), Alu and SVA retrotransposons. Expression analyses demonstrate that cells lacking a functional Microprocessor accumulate LINE-1 mRNA and encoded proteins. Furthermore, we show that structured regions of the LINE-1 mRNA can be cleaved in vitro by Drosha. Additionally, we used a cell culture-based assay to show that the Microprocessor negatively regulates LINE-1 and Alu retrotransposition in vivo. Altogether, these data reveal a new role for the Microprocessor as a post-transcriptional repressor of mammalian retrotransposons and a defender of human genome integrity.