Chaolin Zhang

An RNA-Sequencing Transcriptome and Splicing Database of Glia, Neurons, and Vascular Cells of the Cerebral Cortex

The major cell classes of the brain differ in their developmental processes, metabolism, signaling, and function. To better understand the functions and interactions of the cell types that comprise these classes, we acutely purified representative populations of neurons, astrocytes, oligodendrocyte precursor cells, newly formed oligodendrocytes, myelinating oligodendrocytes, microglia, endothelial cells, and pericytes from mouse cerebral cortex. We generated a transcriptome database for these eight cell types by RNA sequencing and used a sensitive algorithm to detect alternative splicing events in each cell type. Bioinformatic analyses identified thousands of new cell type-enriched genes and splicing isoforms that will provide novel markers for cell identification, tools for genetic manipulation, and insights into the biology of the brain. For example, our data provide clues as to how neurons and astrocytes differ in their ability to dynamically regulate glycolytic flux and lactate generation attributable to unique splicing of PKM2, the gene encoding the glycolytic enzyme pyruvate kinase. This dataset will provide a powerful new resource for understanding the development and function of the brain. To ensure the widespread distribution of these datasets, we have created a user-friendly website (http://web.stanford.edu/group/barres_lab/brain_rnaseq.html) that provides a platform for analyzing and comparing transciption and alternative splicing profiles for various cell classes in the brain.

FMRP Stalls Ribosomal Translocation on mRNAs Linked to Synaptic Function and Autism

FMRP loss of function causes Fragile X syndrome (FXS) and autistic features. FMRP is a polyribosome-associated neuronal RNA-binding protein, suggesting that it plays a key role in regulating neuronal translation, but there has been little consensus regarding either its RNA targets or mechanism of action. Here, we use high-throughput sequencing of RNAs isolated by crosslinking immunoprecipitation (HITS-CLIP) to identify FMRP interactions with mouse brain polyribosomal mRNAs. FMRP interacts with the coding region of transcripts encoding pre- and postsynaptic proteins and transcripts implicated in autism spectrum disorders (ASD). We developed a brain polyribosome-programmed translation system, revealing that FMRP reversibly stalls ribosomes specifically on its target mRNAs. Our results suggest that loss of a translational brake on the synthesis of a subset of synaptic proteins contributes to FXS. In addition, they provide insight into the molecular basis of the cognitive and allied defects in FXS and ASD and suggest multiple targets for clinical intervention.

Mapping in vivo protein-RNA interactions at single-nucleotide resolution from HITS-CLIP data

Mammalian RNA complexity is regulated through interactions of RNA-binding proteins (RBPs) with their target transcripts. High-throughput sequencing together with UV-crosslinking and immunoprecipitation (HITS-CLIP) is able to globally map RBP-binding footprint regions at a resolution of ∼30-60 nucleotides. Here we describe a systematic way to analyze HITS-CLIP data to identify exact crosslink sites, and thereby determine protein-RNA interactions at single-nucleotide resolution. We found that reverse transcriptase used in CLIP frequently skips the crosslinked amino-acid-RNA adduct, resulting in a nucleotide deletion. Genome-wide analysis of these crosslinking-induced mutation sites (CIMS) in HITS-CLIP data for Nova and Argonaute (Ago) proteins in mouse brain tissue revealed deletions in ∼8–20% of mRNA tags, which mapped to Nova and Ago binding sites on mRNA or miRNA. CIMS analysis provides a general and more precise means of mapping protein-RNA interactions than currently available methods and insight into the biochemical properties of such interactions in living tissues.

Defining the regulatory network of the tissue-specific splicing factors Fox-1 and Fox-2.

The precise regulation of many alternative splicing (AS) events by specific splicing factors is essential to determine tissue types and developmental stages. However, the molecular basis of tissue-specific AS regulation and the properties of splicing regulatory networks (SRNs) are poorly understood. Here we comprehensively predict the targets of the brain- and muscle-specific splicing factor Fox-1 (A2BP1) and its paralog Fox-2 (RBM9) and systematically define the corresponding SRNs genome-wide. Fox-1/2 are conserved from worm to human, and specifically recognize the RNA element UGCAUG. We integrate Fox-1/2-binding specificity with phylogenetic conservation, splicing microarray data, and additional computational and experimental characterization. We predict thousands of Fox-1/2 targets with conserved binding sites, at a false discovery rate (FDR) of approximately 24%, including many validated experimentally, suggesting a surprisingly extensive SRN. The preferred position of the binding sites differs according to AS pattern, and determines either activation or repression of exon recognition by Fox-1/2. Many predicted targets are important for neuromuscular functions, and have been implicated in several genetic diseases. We also identified instances of binding site creation or loss in different vertebrate lineages and human populations, which likely reflect fine-tuning of gene expression regulation during evolution.

Dissecting the unique role of the retinoblastoma tumor suppressor during cellular senescence

The RB protein family (RB, p107, and p130) has overlapping and compensatory functions in cell-cycle control. However, cancer-associated mutations are almost exclusively found in RB, implying that RB has a nonredundant role in tumor suppression. We demonstrate that RB preferentially associates with E2F target genes involved in DNA replication and is uniquely required to repress these genes during senescence but not other growth states. Consequently, RB loss leads to inappropriate DNA synthesis following a senescence trigger and, together with disruption of a p21-mediated cell-cycle checkpoint, enables extensive proliferation and rampant genomic instability. Our results identify a nonredundant RB effector function that may contribute to tumor suppression and reveal how loss of RB and p53 cooperate to bypass senescence.

Integrative modeling defines the Nova splicing-regulatory network and its combinatorial controls.

Making the Final Cut RNA splicing, which involves selectively cutting and pasting messenger RNA to generate different proteins, is critical in regulating human physiology and diseases. However, our knowledge of the underlying rules governing splicing regulation remains incomplete. The recent emergence of next-generation sequencing and other high-throughput technologies has provided an opportunity to transform our understanding of RNA regulation. Zhang et al. (p. 439, published online 17 June) combined multiple data sets to generate a robust and relatively complete picture of splicing regulation by the mammalian neuronal splicing factor Nova in the brain. About 700 splicing events were identified, including many novel target exons, some likely to be involved in neurological disease. Combining genomic studies with computational biology also yielded insight into the regulation of alternative splicing. Network modeling reveals an RNA splicing network in the brain and uncovers complex mechanisms of regulation. The control of RNA alternative splicing is critical for generating biological diversity. Despite emerging genome-wide technologies to study RNA complexity, reliable and comprehensive RNA-regulatory networks have not been defined. Here, we used Bayesian networks to probabilistically model diverse data sets and predict the target networks of specific regulators. We applied this strategy to identify ~700 alternative splicing events directly regulated by the neuron-specific factor Nova in the mouse brain, integrating RNA-binding data, splicing microarray data, Nova-binding motifs, and evolutionary signatures. The resulting integrative network revealed combinatorial regulation by Nova and the neuronal splicing factor Fox, interplay between phosphorylation and splicing, and potential links to neurologic disease. Thus, we have developed a general approach to understanding mammalian RNA regulation at the systems level.

Muscleblind-like 2-Mediated Alternative Splicing in the Developing Brain and Dysregulation in Myotonic Dystrophy

The RNA-mediated disease model for myotonic dystrophy (DM) proposes that microsatellite C(C)TG expansions express toxic RNAs that disrupt splicing regulation by altering MBNL1 and CELF1 activities. While this model explains DM manifestations in muscle, less is known about the effects of C(C)UG expression on the brain. Here, we report that Mbnl2 knockout mice develop several DM-associated central nervous system (CNS) features including abnormal REM sleep propensity and deficits in spatial memory. Mbnl2 is prominently expressed in the hippocampus and Mbnl2 knockouts show a decrease in NMDA receptor (NMDAR) synaptic transmission and impaired hippocampal synaptic plasticity. While Mbnl2 loss did not significantly alter target transcript levels in the hippocampus, misregulated splicing of hundreds of exons was detected using splicing microarrays, RNA-seq, and HITS-CLIP. Importantly, the majority of the Mbnl2-regulated exons examined were similarly misregulated in DM. We propose that major pathological features of the DM brain result from disruption of the MBNL2-mediated developmental splicing program.

Ptbp2 represses adult-specific splicing to regulate the generation of neuronal precursors in the embryonic brain

Two polypyrimidine tract RNA-binding proteins (PTBs), one near-ubiquitously expressed (Ptbp1) and another highly tissue-restricted (Ptbp2), regulate RNA in interrelated but incompletely understood ways. Ptbp1, a splicing regulator, is replaced in the brain and differentiated neuronal cell lines by Ptbp2. To define the roles of Ptbp2 in the nervous system, we generated two independent Ptbp2-null strains, unexpectedly revealing that Ptbp2 is expressed in neuronal progenitors and is essential for postnatal survival. A HITS-CLIP (high-throughput sequencing cross-linking immunoprecipitation)-generated map of reproducible Ptbp2-RNA interactions in the developing mouse neocortex, combined with results from splicing-sensitive microarrays, demonstrated that the major action of Ptbp2 is to inhibit adult-specific alternative exons by binding pyrimidine-rich sequences upstream of and/or within them. These regulated exons are present in mRNAs encoding proteins associated with control of cell fate, proliferation, and the actin cytoskeleton, suggesting a role for Ptbp2 in neurogenesis. Indeed, neuronal progenitors in the Ptbp2-null brain exhibited an aberrant polarity and were associated with regions of premature neurogenesis and reduced progenitor pools. Thus, Ptbp2 inhibition of a discrete set of adult neuronal exons underlies early brain development prior to neuronal differentiation and is essential for postnatal survival.

Mapping Argonaute and conventional RNA-binding protein interactions with RNA at single-nucleotide resolution using HITS-CLIP and CIMS analysis

The identification of sites where RNA-binding proteins (RNABPs) interact with target RNAs opens the door to understanding the vast complexity of RNA regulation. UV cross-linking and immunoprecipitation (CLIP) is a transformative technology in which RNAs purified from in vivo cross-linked RNA-protein complexes are sequenced to reveal footprints of RNABP:RNA contacts. CLIP combined with high-throughput sequencing (HITS-CLIP) is a generalizable strategy to produce transcriptome-wide maps of RNA binding with higher accuracy and resolution than standard RNA immunoprecipitation (RIP) profiling or purely computational approaches. The application of CLIP to Argonaute proteins has expanded the utility of this approach to mapping binding sites for microRNAs and other small regulatory RNAs. Finally, recent advances in data analysis take advantage of cross-link–induced mutation sites (CIMS) to refine RNA-binding maps to single-nucleotide resolution. Once IP conditions are established, HITS-CLIP takes ∼8 d to prepare RNA for sequencing. Established pipelines for data analysis, including those for CIMS, take 3–4 d.

Profiling alternatively spliced mRNA isoforms for prostate cancer classification

BackgroundProstate cancer is one of the leading causes of cancer illness and death among men in the United States and world wide. There is an urgent need to discover good biomarkers for early clinical diagnosis and treatment. Previously, we developed an exon-junction microarray-based assay and profiled 1532 mRNA splice isoforms from 364 potential prostate cancer related genes in 38 prostate tissues. Here, we investigate the advantage of using splice isoforms, which couple transcriptional and splicing regulation, for cancer classification.ResultsAs many as 464 splice isoforms from more than 200 genes are differentially regulated in tumors at a false discovery rate (FDR) of 0.05. Remarkably, about 30% of genes have isoforms that are called significant but do not exhibit differential expression at the overall mRNA level. A support vector machine (SVM) classifier trained on 128 signature isoforms can correctly predict 92% of the cases, which outperforms the classifier using overall mRNA abundance by about 5%. It is also observed that the classification performance can be improved using multivariate variable selection methods, which take correlation among variables into account.ConclusionThese results demonstrate that profiling of splice isoforms is able to provide unique and important information which cannot be detected by conventional microarrays.