Transcriptome variation in human tissues revealed by long-read sequencing.

Regulation of transcript structure generates transcript diversity and plays an important role in human disease1-7. The advent of long-read sequencing technologies offers the opportunity to study the role of genetic variation in transcript structure8-16. In this Article, we present a large human long-read RNA-seq dataset using the Oxford Nanopore Technologies platform from 88 samples from Genotype-Tissue Expression (GTEx) tissues and cell lines, complementing the GTEx resource. We identified just over 70,000 novel transcripts for annotated genes, and validated the protein expression of 10% of novel transcripts. We developed a new computational package, LORALS, to analyse the genetic effects of rare and common variants on the transcriptome by allele-specific analysis of long reads. We characterized allele-specific expression and transcript structure events, providing new insights into the specific transcript alterations caused by common and rare genetic variants and highlighting the resolution gained from long-read data. We were able to perturb the transcript structure upon knockdown of PTBP1, an RNA binding protein that mediates splicing, thereby finding genetic regulatory effects that are modified by the cellular environment. Finally, we used this dataset to enhance variant interpretation and study rare variants leading to aberrant splicing patterns.

Reversing a model of Parkinson's disease with in situ converted nigral neurons.

Parkinson's disease is characterized by loss of dopamine neurons in the substantia nigra1. Similar to other major neurodegenerative disorders, there are no disease-modifying treatments for Parkinson's disease. While most treatment strategies aim to prevent neuronal loss or protect vulnerable neuronal circuits, a potential alternative is to replace lost neurons to reconstruct disrupted circuits2. Here we report an efficient one-step conversion of isolated mouse and human astrocytes to functional neurons by depleting the RNA-binding protein PTB (also known as PTBP1). Applying this approach to the mouse brain, we demonstrate progressive conversion of astrocytes to new neurons that innervate into and repopulate endogenous neural circuits. Astrocytes from different brain regions are converted to different neuronal subtypes. Using a chemically induced model of Parkinson's disease in mouse, we show conversion of midbrain astrocytes to dopaminergic neurons, which provide axons to reconstruct the nigrostriatal circuit. Notably, re-innervation of striatum is accompanied by restoration of dopamine levels and rescue of motor deficits. A similar reversal of disease phenotype is also accomplished by converting astrocytes to neurons using antisense oligonucleotides to transiently suppress PTB. These findings identify a potentially powerful and clinically feasible approach to treating neurodegeneration by replacing lost neurons.

RNA-binding protein Ptbp1 is essential for BCR-mediated antibody production.

The RNA-binding protein polypyrimidine tract-binding protein-1 (Ptbp1) binds to the pyrimidine-rich sequence of target RNA and controls gene expression via post-transcriptional regulation such as alternative splicing. Although Ptbp1 is highly expressed in B lymphocytes, its role to date is largely unknown. To clarify the role of Ptbp1 in B-cell development and function, we generated B-cell-specific Ptbp1-deficient (P1BKO) mice. B-cell development in the bone marrow, spleen and peritoneal cavity of the P1BKO mice was nearly normal. However, the P1BKO mice had significantly lower levels of natural antibodies in serum compared with those of the control mice. To investigate the effect of Ptbp1 deficiency on the immune response in vivo, we immunized the P1BKO mice with T-cell-independent type-2 (TI-2) antigen NP-Ficoll and T-cell-dependent (TD) antigen NP-CGG. We found that B-cell-specific Ptbp1 deficiency causes an immunodeficiency phenotype due to defective production of antibody against both TI-2 and TD antigen. This immunodeficiency was accompanied by impaired B-cell receptor (BCR)-mediated B-cell activation and plasmablast generation. These findings demonstrate that Ptbp1 is essential for the humoral immune response.

PTBP1 contributes to spermatogenesis through regulation of proliferation in spermatogonia.

Polypyrimidine tract-binding protein 1 (PTBP1) is a highly conserved RNA-binding protein that is a well-known regulator of alternative splicing. Testicular tissue is one of the richest tissues with respect to the number of alternative splicing mRNA isoforms, but the molecular role(s) of PTBP1 in the regulation of these isoforms during spermatogenesis is still unclear. Here, we developed a germ cell-specific Ptbp1 conditional knockout (cKO) mouse model by using the Cre-loxP system to investigate the role of PTBP1 in spermatogenesis. Testis weight in Ptbp1 cKO mice was comparable to that in age-matched controls until 3 weeks of age; at ≥ 2 months old, testis weight was significantly lighter in cKO mice than in age-matched controls. Sperm count in Ptbp1 cKO mice at 2 months old was comparable to that in controls, whereas sperm count significantly decreased at 6 months old. Seminiferous tubules that exhibited degeneration in spermatogenic function were more evident in the 2-month-old Ptbp1 cKO mice than in controls. In addition, the early neonatal proliferation of spermatogonia, during postnatal days 1-5, was significantly retarded in Ptbp1 cKO mice compared with that in controls. An in vitro spermatogonia culture model (germline stem cells) revealed that hydroxytamoxifen-induced deletion of PTBP1 from germline stem cells caused severe proliferation arrest accompanied by an increase of apoptotic cell death. These data suggest that PTBP1 contributes to spermatogenesis through regulation of spermatogonia proliferation.

Single-cell transcriptomics reconstructs fate conversion from fibroblast to cardiomyocyte.

Direct lineage conversion offers a new strategy for tissue regeneration and disease modelling. Despite recent success in directly reprogramming fibroblasts into various cell types, the precise changes that occur as fibroblasts progressively convert to the target cell fates remain unclear. The inherent heterogeneity and asynchronous nature of the reprogramming process renders it difficult to study this process using bulk genomic techniques. Here we used single-cell RNA sequencing to overcome this limitation and analysed global transcriptome changes at early stages during the reprogramming of mouse fibroblasts into induced cardiomyocytes (iCMs). Using unsupervised dimensionality reduction and clustering algorithms, we identified molecularly distinct subpopulations of cells during reprogramming. We also constructed routes of iCM formation, and delineated the relationship between cell proliferation and iCM induction. Further analysis of global gene expression changes during reprogramming revealed unexpected downregulation of factors involved in mRNA processing and splicing. Detailed functional analysis of the top candidate splicing factor, Ptbp1, revealed that it is a critical barrier for the acquisition of cardiomyocyte-specific splicing patterns in fibroblasts. Concomitantly, Ptbp1 depletion promoted cardiac transcriptome acquisition and increased iCM reprogramming efficiency. Additional quantitative analysis of our dataset revealed a strong correlation between the expression of each reprogramming factor and the progress of individual cells through the reprogramming process, and led to the discovery of new surface markers for the enrichment of iCMs. In summary, our single-cell transcriptomics approaches enabled us to reconstruct the reprogramming trajectory and to uncover intermediate cell populations, gene pathways and regulators involved in iCM induction.

Concordance of deregulated mechanisms unveiled in underpowered experiments: PTBP1 knockdown case study.

BACKGROUND: Genome-wide transcriptome profiling generated by microarray and RNA-Seq often provides deregulated genes or pathways applicable only to larger cohort. On the other hand, individualized interpretation of transcriptomes is increasely pursued to improve diagnosis, prognosis, and patient treatment processes. Yet, robust and accurate methods based on a single paired-sample remain an unmet challenge. METHODS: "N-of-1-pathways" translates gene expression data profiles into mechanism-level profiles on single pairs of samples (one p-value per geneset). It relies on three principles: i) statistical universe is a single paired sample, which serves as its own control; ii) statistics can be derived from multiple gene expression measures that share common biological mechanisms assimilated to genesets; iii) semantic similarity metric takes into account inter-mechanisms' relationships to better assess commonality and differences, within and cross study-samples (e.g. patients, cell-lines, tissues, etc.), which helps the interpretation of the underpinning biology. RESULTS: In the context of underpowered experiments, N-of-1-pathways predictions perform better or comparable to those of GSEA and Differentially Expressed Genes enrichment (DEG enrichment), within-and cross-datasets. N-of-1-pathways uncovered concordant PTBP1-dependent mechanisms across datasets (Odds-Ratios≥13, p-values≤1 × 10-5), such as RNA splicing and cell cycle. In addition, it unveils tissue-specific mechanisms of alternatively transcribed PTBP1-dependent genesets. Furthermore, we demonstrate that GSEA and DEG Enrichment preclude accurate analysis on single paired samples. CONCLUSIONS: N-of-1-pathways enables robust and biologically relevant mechanism-level classifiers with small cohorts and one single paired samples that surpasses conventional methods. Further, it identifies unique sample/ patient mechanisms, a requirement for precision medicine.

Drosophila polypyrimidine tract-binding protein (DmPTB) regulates dorso-ventral patterning genes in embryos.

The Drosophila polypyrimidine tract-binding protein (dmPTB or hephaestus) plays an important role during embryogenesis. A loss of function mutation, heph(03429), results in varied defects in embryonic developmental processes, leading to embryonic lethality. However, the suite of molecular functions that are disrupted in the mutant remains unknown. We have used an unbiased high throughput sequencing approach to identify transcripts that are misregulated in this mutant. Misregulated transcripts show evidence of significantly altered patterns of splicing (exon skipping, 5' and 3' splice site switching), alternative 5' ends, and mRNA level changes (up and down regulation). These findings are independently supported by reverse-transcription-polymerase chain reaction (RT-PCR) analysis and in situ hybridization. We show that a group of genes, such as Zerknüllt, z600 and screw are among the most upregulated in the mutant and have been functionally linked to dorso-ventral patterning and/or dorsal closure processes. Thus, loss of dmPTB function results in specific misregulated transcripts, including those that provide the missing link between the loss of dmPTB function and observed developmental defects in embryogenesis. This study provides the first comprehensive repertoire of genes affected in vivo in the heph mutant in Drosophila and offers insight into the role of dmPTB during embryonic development.

The polypyrimidine tract binding protein regulates desaturase alternative splicing and PUFA composition.

The Δ6 desaturase, encoded by FADS2, plays a crucial role in omega-3 and omega-6 fatty acid synthesis. These fatty acids are essential components of the central nervous system, and they act as precursors for eicosanoid signaling molecules and as direct modulators of gene expression. The polypyrimidine tract binding protein (PTB or hnRNP I) is a splicing factor that regulates alternative pre-mRNA splicing. Here, PTB is shown to bind an exonic splicing silencer element and repress alternative splicing of FADS2 into FADS2 AT1. PTB and FADS2AT1 were inversely correlated in neonatal baboon tissues, implicating PTB as a major regulator of tissue-specific FADS2 splicing. In HepG2 cells, PTB knockdown modulated alternative splicing of FADS2, as well as FADS3, a putative desaturase of unknown function. Omega-3 fatty acids decreased by nearly one half relative to omega-6 fatty acids in PTB knockdown cells compared with controls, with a particularly strong decrease in eicosapentaenoic acid (EPA) concentration and its ratio to arachidonic acid (ARA). This is a rare demonstration of a mechanism specifically altering the cellular omega-3 to omega-6 fatty acid ratio without any change in diet/media. These findings reveal a novel role for PTB, regulating availability of membrane components and eicosanoid precursors for cell signaling.

Polypyrimidine tract-binding protein 1 regulates the alternative splicing of dopamine receptor D2.

Dopamine receptor D(2) (DRD2) has two splicing isoforms, a long form (D2L) and short form (D2S), which have distinct functions in the dopaminergic system. However, the regulatory mechanism of the alternative splicing of DRD2 is unknown. In this study, we examined which splicing factors regulate the expression of D2L and D2S by over-expressing several RNA-binding proteins in HEK293 cells. In a cellular splicing assay, the over-expression of polypyrimidine tract-binding protein 1 (PTBP1) reduced the expression of D2S, whereas the knockdown of PTBP1 increased the expression of D2S. We also identified the regions of DRD2 that are responsive to PTBP1 using heterologous minigenes and deletion mutants. Our results indicate that PTBP1 regulates the alternative splicing of DRD2. Considering that DRD2 inhibits cAMP-dependent protein kinase A, which modulates the intracellular localization of PTBP1, PTBP1 may contribute to the autoregulation of DRD2 by regulating the expression of its isoforms.

RNAi-mediated PTB depletion leads to enhanced exon definition.

Mutually exclusive use of exons IIIb or IIIc in FGF-R2 transcripts requires the silencing of exon IIIb. This repression is mediated by silencer elements upstream and downstream of the exon. Both silencers bind the polypyrimidine tract binding protein (PTB) and PTB binding sites within these elements are required for efficient silencing of exon IIIb. Recruitment of MS2-PTB fusion proteins upstream or downstream of exon IIIb causes repression of this exon. Depletion of endogenous PTB using RNAi increases exon IIIb inclusion in transcripts derived from minigenes and from the endogenous FGF-R2 gene. These data demonstrate that PTB is a negative regulator of exon definition in vivo.