Loss of PTEN, the negative regulator of PI3K activity, is frequent in glioblastomas (GBMs). However, the role of the two major PI3K isoforms, p110α and p110β, in PTEN-deficient gliomagenesis remains unknown. We show that PTEN-deficient GBM largely depends on p110α for proliferation and p110β for migration. Genetic ablation of either isoform delays tumor progression in mice, but only ablating both isoforms completely blocks GBM driven by the concurrent ablation of Pten and p53. BKM120 (buparlisib) treatment only modestly prolongs survival in mice bearing intracranial Pten/p53 null tumors due to partial pathway inhibition. BKM120 extends the survival of mice bearing intracranial tumors in which p110β, but not p110α, has been genetically ablated in the Pten/p53 null glioma, indicating that BKM120 fails to inhibit p110β effectively. Our study suggests that the failure of PI3K inhibitors in GBM may be due to insufficient inhibition of p110β and indicates a need to develop brain-penetrant p110α/β inhibitors.

Alternative splicing (AS) is an important source of proteome diversity in eukaryotes. However, how this affects protein repertoires at a single-cell level remains an open question. Here, we show that many 3'-terminal exons are persistently co-expressed with their alternatives in mammalian neurons. In an important example of this scenario, cell polarity gene Cdc42, a combination of polypyrimidine tract-binding, protein-dependent, and constitutive splicing mechanisms ensures a halfway switch from the general (E7) to the neuron-specific (E6) alternative 3'-terminal exon during neuronal differentiation. Perturbing the nearly equimolar E6/E7 ratio in neurons results in defects in both axonal and dendritic compartments and suggests that Cdc42E7 is involved in axonogenesis, whereas Cdc42E6 is required for normal development of dendritic spines. Thus, co-expression of a precise blend of functionally distinct splice isoforms rather than a complete switch from one isoform to another underlies proper structural and functional polarization of neurons.

The protein-level translational status and function of many alternative splicing events remain poorly understood. We use an RNA sequencing (RNA-seq)-guided proteomics method to identify protein alternative splicing isoforms in the human proteome by constructing tissue-specific protein databases that prioritize transcript splice junction pairs with high translational potential. Using the custom databases to reanalyze ∼80 million mass spectra in public proteomics datasets, we identify more than 1,500 noncanonical protein isoforms across 12 human tissues, including ∼400 sequences undocumented on TrEMBL and RefSeq databases. We apply the method to original quantitative mass spectrometry experiments and observe widespread isoform regulation during human induced pluripotent stem cell cardiomyocyte differentiation. On a proteome scale, alternative isoform regions overlap frequently with disordered sequences and post-translational modification sites, suggesting that alternative splicing may regulate protein function through modulating intrinsically disordered regions. The described approach may help elucidate functional consequences of alternative splicing and expand the scope of proteomics investigations in various systems.

The epithelial-specific splicing regulators Esrp1 and Esrp2 are required for mammalian development, including establishment of epidermal barrier functions. However, the mechanisms by which Esrp ablation causes defects in epithelial barriers remain undefined. We determined that the ablation of Esrp1 and Esrp2 impairs epithelial tight junction (TJ) integrity through loss of the epithelial isoform of Rho GTP exchange factor Arhgef11. Arhgef11 is required for the maintenance of TJs via RhoA activation and myosin light chain (MLC) phosphorylation. Ablation or depletion of Esrp1/2 or Arhgef11 inhibits MLC phosphorylation and only the epithelial Arhgef11 isoform rescues MLC phosphorylation in Arhgef11 KO epithelial cells. Mesenchymal Arhgef11 transcripts contain a C-terminal exon that binds to PAK4 and inhibits RhoA activation byArhgef11. Deletion of the mesenchymal-specific Arhgef11 exon in Esrp1/2 KO epithelial cells using CRISPR/Cas9 restored TJ function, illustrating how splicing alterations can be mechanistically linked to disease phenotypes that result from impaired functions of splicing regulators.

The regulation of hepatic gene expression has been extensively studied at the transcriptional level; however, the control of metabolism through posttranscriptional gene regulation by RNA-binding proteins in physiological and disease states is less understood. Here, we report a major role for the hormone-sensitive RNA-binding protein (RBP) APOBEC1 complementation factor (A1CF) in the generation of hepatocyte-specific and alternatively spliced transcripts. Among these transcripts are isoforms for the dominant and high-affinity fructose-metabolizing ketohexokinase C and glycerol kinase, two key metabolic enzymes that are linked to hepatic gluconeogenesis and found to be markedly reduced upon hepatic ablation of A1cf. Consequently, mice lacking A1CF exhibit improved glucose tolerance and are protected from fructose-induced hyperglycemia, hepatic steatosis, and development of obesity. Our results identify a previously unreported function of A1CF as a regulator of alternative splicing of a subset of genes influencing hepatic glucose production through fructose and glycerol metabolism.

Adenosine-to-inosine RNA editing, catalyzed by adenosine deaminase acting on RNA (ADAR) enzymes, alters RNA sequences from those encoded by DNA. These editing events are dynamically regulated, but few trans regulators of ADARs are known in vivo. Here, we screen RNA-binding proteins for roles in editing regulation with knockdown experiments in the Drosophila brain. We identify zinc-finger protein at 72D (Zn72D) as a regulator of editing levels at a majority of editing sites in the brain. Zn72D both regulates ADAR protein levels and interacts with ADAR in an RNA-dependent fashion, and similar to ADAR, Zn72D is necessary to maintain proper neuromuscular junction architecture and fly mobility. Furthermore, Zn72D's regulatory role in RNA editing is conserved because the mammalian homolog of Zn72D, Zfr, regulates editing in mouse primary neurons. The broad and conserved regulation of ADAR editing by Zn72D in neurons sustains critically important editing events.

Methods for single-cell RNA sequencing (scRNA-seq) have greatly advanced in recent years. While droplet- and well-based methods have increased the capture frequency of cells for scRNA-seq, these technologies readily produce technical artifacts, such as doublet cell captures. Doublets occurring between distinct cell types can appear as hybrid scRNA-seq profiles, but do not have distinct transcriptomes from individual cell states. We introduce DoubletDecon, an approach that detects doublets with a combination of deconvolution analyses and the identification of unique cell-state gene expression. We demonstrate the ability of DoubletDecon to identify synthetic, mixed-species, genetic, and cell-hashing cell doublets from scRNA-seq datasets of varying cellular complexity with a high sensitivity relative to alternative approaches. Importantly, this algorithm prevents the prediction of valid mixed-lineage and transitional cell states as doublets by considering their unique gene expression. DoubletDecon has an easy-to-use graphical user interface and is compatible with diverse species and unsupervised population detection algorithms.

Emerging evidence suggests that long non-coding RNAs (lncRNAs) play diverse and critical roles in neural development, function, and disease. Here, we examine neuronal lncRNAs in a model system that offers enormous advantages for deciphering their functions: the Drosophila olfactory system. This system is numerically simple, its neurons are exquisitely well defined, and it drives multiple complex behaviors. We undertake a comprehensive survey of linear and circular lncRNAs in the Drosophila antenna and identify a wealth of lncRNAs enriched in it. We generate an unprecedented lncRNA-to-neuron map, which reveals that olfactory receptor neurons are defined not only by their receptors but also by the combination of lncRNAs they express. We identify species-specific lncRNAs, including many that are expressed primarily in pheromone-sensing neurons and that may act in modulation of pheromonal responses or in speciation. This resource opens many new opportunities for investigating the roles of lncRNAs in the nervous system.

The RNA N6-methyladenosine (m6A) methylation is installed by the METTL3-METTL14 methyltransferase complex. This modification has critical regulatory roles in various biological processes. Here, we report that deletion of Mettl14 dramatically reduces mRNA m6A methylation in developing B cells and severely blocks B cell development in mice. Deletion of Mettl14 impairs interleukin-7 (IL-7)-induced pro-B cell proliferation and the large-pre-B-to-small-pre-B transition and causes dramatic abnormalities in gene expression programs important for B cell development. Suppression of a group of transcripts by cytoplasmic m6A reader YTHDF2 is critical to the IL-7-induced pro-B cell proliferation. In contrast, the block in the large-pre-B-to-small-pre-B transition is independent of YTHDF1 or YTHDF2 but is associated with a failure to properly upregulate key transcription factors regulating this transition. Our data highlight the important regulatory roles of the RNA m6A methylation and its reader proteins in early B cell development.

Degradation of transcripts in human nuclei is primarily facilitated by the RNA exosome. To obtain substrate specificity, the exosome is aided by adaptors; in the nucleoplasm, those adaptors are the nuclear exosome-targeting (NEXT) complex and the poly(A) (pA) exosome-targeting (PAXT) connection. How these adaptors guide exosome targeting remains enigmatic. Employing high-resolution 3' end sequencing, we demonstrate that NEXT substrates arise from heterogenous and predominantly pA- 3' ends often covering kilobase-wide genomic regions. In contrast, PAXT targets harbor well-defined pA+ 3' ends defined by canonical pA site use. Irrespective of this clear division, NEXT and PAXT act redundantly in two ways: (1) regional redundancy, where the majority of exosome-targeted transcription units produce NEXT- and PAXT-sensitive RNA isoforms, and (2) isoform redundancy, where the PAXT connection ensures fail-safe decay of post-transcriptionally polyadenylated NEXT targets. In conjunction, this provides a two-layered targeting mechanism for efficient nuclear sorting of the human transcriptome.

RNA editing is a process by which post-transcriptional changes of mRNA nucleotides alter protein function through modification of the amino acid content. The 5HT2C serotonin receptor, which undergoes 32 distinct RNA-editing events leading to 24 protein isoforms, is a notable example of this process. These 5HT2C isoforms display differences in constitutive activity, agonist/inverse agonist potencies, and efficacies. To elucidate the molecular mechanisms responsible for these effects of RNA editing, we present four active-state 5HT2C-transducer-coupled structures of three representative isoforms (INI, VGV, and VSV) with the selective drug lorcaserin (Belviq) and the classic psychedelic psilocin. We also provide a comprehensive analysis of agonist activation and constitutive activity across all 24 protein isoforms. Collectively, these findings reveal a unique hydrogen-bonding network located on intracellular loop 2 that is subject to RNA editing, which differentially affects GPCR constitutive and agonist signaling activities.

Neurexins are key synaptic organizers that are expressed in thousands of alternatively spliced isoforms. Because transsynaptic neurexin interactions with different postsynaptic molecules are largely isoform dependent, a cell type-level census of different neurexin isoforms could predict molecular interactions relating to synapse identity and function. Using single-cell transcriptomics to study the origin of neurexin diversity in multiple murine mature and embryonic cell types, we have discovered shared neurexin expression patterns in developmentally related cells. By comparing neurexin profiles in immature embryonic neurons, we show that neurexin profiles are specified during early development and remain unchanged throughout neuronal maturation. Thus, our findings reveal ontogenetic stability and provide a cell type-level census of neurexin isoform expression in the cortex.

Organoids allow the recapitulation of intestinal homeostasis and cancerogenesis in vitro; however, RNA sequencing (RNA-seq)-based methods for drug screens are missing. We develop targeted organoid sequencing (TORNADO-seq), a high-throughput, high-content drug discovery platform that uses targeted RNA-seq to monitor the expression of large gene signatures for the detailed evaluation of cellular phenotypes in organoids. TORNADO-seq is a fast, highly reproducible time- and cost-effective ($5 per sample) method that can probe cell mixtures and their differentiation state in the intestinal system. We apply this method to isolate drugs that enrich for differentiated cell phenotypes and show that these drugs are highly efficacious against cancer compared to wild-type organoids. Furthermore, TORNADO-seq facilitates in-depth insight into the mode of action of these drugs. Our technology can easily be adapted to many other systems and will allow for more systematic, large-scale, and quantitative approaches to study the biology of complex cellular systems.

Although adult neurogenesis recapitulates processes that occur during embryonic development, it exhibits distinct characteristics from the embryonic counterpart. However, the intrinsic mechanism underlying the differential regulation of neurogenesis between these two stages remains unclear. Herein, we show that the ablation of RNA-binding protein HuR in NSCs impairs adult but not embryonic neurogenesis. HuR is predominantly expressed in the cytoplasm of embryonic NSCs but translocates into the nucleus of adult NSCs. Transcriptomic analysis of HuR-deficient adult NSCs revealed that HuR primarily regulates alternative splicing of numerous premRNA transcripts, including focal adhesion kinase (FAK). HuR-deficient adult NSCs generate increased FAK mRNA isoforms with shorter 5'-UTRs, leading to enhanced FAK mRNA translation and hyperactivated FAK signaling, and inhibition of FAK ameliorates defective adult neurogenesis and impaired hippocampus-dependent learning in HuR-deficient mice. These findings provide mechanistic insights into the differential regulation of embryonic and adult neurogenesis through developmental cytoplasmic-to-nuclear translocation of HuR.

Drosophila melanogaster Down syndrome cell adhesion molecule (Dscam1) can generate 38,016 different isoforms through largely stochastic, yet highly biased, alternative splicing. These isoforms are required for nervous functions. However, the functional significance of splicing bias remains unknown. Here, we provide evidence that Dscam1 splicing bias is required for mushroom body (MB) axonal wiring. We generate mutant flies with normal overall protein levels and an identical number but global changes in exon 4 and 9 isoform bias (DscamΔ4D-/- and DscamΔ9D-/-), respectively. In contrast to DscamΔ4D-/-, DscamΔ9D-/- exhibits remarkable MB defects, suggesting a variable domain-specific requirement for isoform bias. Importantly, changes in isoform bias cause axonal defects but do not influence the self-avoidance of axonal branches. We conclude that, in contrast to the isoform number that provides the molecular basis for neurite self-avoidance, isoform bias may play a role in MB axonal wiring by influencing non-repulsive signaling.

Linker histones are highly abundant chromatin-associated proteins with well-established structural roles in chromatin and as general transcriptional repressors. In addition, it has been long proposed that histone H1 exerts context-specific effects on gene expression. Here, we identify a function of histone H1 in chromatin structure and transcription using a range of genomic approaches. In the absence of histone H1, there is an increase in the transcription of non-coding RNAs, together with reduced levels of m6A modification leading to their accumulation on chromatin and causing replication-transcription conflicts. This strongly suggests that histone H1 prevents non-coding RNA transcription and regulates non-coding transcript turnover on chromatin. Accordingly, altering the m6A RNA methylation pathway rescues the replicative phenotype of H1 loss. This work unveils unexpected regulatory roles of histone H1 on non-coding RNA turnover and m6A deposition, highlighting the intimate relationship between chromatin conformation, RNA metabolism, and DNA replication to maintain genome performance.

Ependymoma (EPN) is a brain tumor commonly presenting in childhood that remains fatal in most children. Intra-tumoral cellular heterogeneity in bulk-tumor samples significantly confounds our understanding of EPN biology, impeding development of effective therapy. We, therefore, use single-cell RNA sequencing, histology, and deconvolution to catalog cellular heterogeneity of the major childhood EPN subgroups. Analysis of PFA subgroup EPN reveals evidence of an undifferentiated progenitor subpopulation that either differentiates into subpopulations with ependymal cell characteristics or transitions into a mesenchymal subpopulation. Histological analysis reveals that progenitor and mesenchymal subpopulations co-localize in peri-necrotic zones. In conflict with current classification paradigms, relative PFA subpopulation proportions are shown to determine bulk-tumor-assigned subgroups. We provide an interactive online resource that facilitates exploration of the EPN single-cell dataset. This atlas of EPN cellular heterogeneity increases understanding of EPN biology.

The nuclear exosome targeting (NEXT) complex is responsible for specific nuclear RNA degradation in mammalian cells. However, its function in development remains unknown. Here, we find that the depletion of a central factor of the NEXT complex, Zcchc8, in mouse results in developmental defects, a shortened lifespan, and infertility. We find that Zcchc8-deficient embryonic stem cells (ESCs) exhibit proliferation abnormalities and reduced developmental potencies. Importantly, the transcripts of retrotransposon element LINE1 are found to be targeted by Zcchc8 and degraded by a Zcchc8-mediated mechanism. We further find that sustained expression of higher levels of LINE1 RNA is detected in maternal Zcchc8-depleted oocytes and embryos. Zcchc8-depleted oocytes show higher chromatin accessibility and developmental defects in both meiotic maturation and embryogenesis after fertilization. Collectively, our study defines Zcchc8-mediated RNA degradation as an important post-transcription regulation of LINE1 transcripts in early embryos and ESCs, which play vital roles in the pluripotency and early development.

Silencing of transposable elements (TEs) drives the evolution of numerous redundant mechanisms of transcriptional regulation. Arabidopsis MBD5, MBD6, and SILENZIO act as TE repressors downstream of DNA methylation. Here, we show, via single-nucleus RNA-seq of developing male gametophytes, that these repressors are critical for TE silencing in the pollen vegetative cell, a companion cell important for fertilization that undergoes chromatin decompaction. Instead, other silencing mutants (met1, ddm1, mom1, morc) show loss of silencing in all pollen nucleus types and somatic cells. We show that TEs repressed by MBD5/6 gain chromatin accessibility in wild-type vegetative nuclei despite remaining silent, suggesting that loss of DNA compaction makes them sensitive to loss of MBD5/6. Consistently, crossing mbd5/6 to histone 1 mutants, which have decondensed chromatin in leaves, reveals derepression of MBD5/6-dependent TEs in leaves. MBD5/6 and SILENZIO thus act as a silencing system particularly important when chromatin compaction is compromised.

Longitudinals lacking (lola) is one of the most complex genes in Drosophila melanogaster, encoding up to 20 protein isoforms that include key transcription factors involved in axonal pathfinding and neural reprogramming. Most previous studies have employed loss-of-function alleles that disrupt lola common exons, making it difficult to delineate isoform-specific functions. To overcome this issue, we have generated isoform-specific mutants for all isoforms using CRISPR/Cas9. This enabled us to study specific isoforms with respect to previously characterized roles for Lola and to demonstrate a specific function for one variant in axon guidance via activation of the microtubule-associated factor Futsch. Importantly, we also reveal a role for a second variant in preventing neurodegeneration via the positive regulation of a key enzyme of the octopaminergic pathway. Thus, our comprehensive study expands the functional repertoire of Lola functions, and it adds insights into the regulatory control of neurotransmitter expression in vivo.