Analysis of alternative cleavage and polyadenylation in mature and differentiating neurons using RNA-seq data

Background: Most eukaryotic protein-coding genes exhibit alternative cleavage and polyadenylation (APA), resulting in mRNA isoforms with different 3′ untranslated regions (3′ UTRs). Studies have shown that brain cells tend to express long 3′ UTR isoforms using distal cleavage and polyadenylation sites (PASs). Methods: Using our recently developed, comprehensive PAS database PolyA_DB, we developed an efficient method to examine APA, named Significance Analysis of Alternative Polyadenylation using RNA-seq (SAAP-RS). We applied this method to study APA in brain cells and neurogenesis. Results: We found that neurons globally express longer 3′ UTRs than other cell types in brain, and microglia and endothelial cells express substantially shorter 3′ UTRs. We show that the 3′ UTR diversity across brain cells can be corroborated with single cell sequencing data. Further analysis of APA regulation of 3′ UTRs during differentiation of embryonic stem cells into neurons indicates that a large fraction of the APA events regulated in neurogenesis are similarly modulated in myogenesis, but to a much greater extent. Conclusion: Together, our data delineate APA profiles in different brain cells and indicate that APA regulation in neurogenesis is largely an augmented process taking place in other types of cell differentiation.     

CRISPR-iPAS: a novel dCAS13-based method for alternative polyadenylation interference

Alternative polyadenylation (APA) plays an important role in gene regulation. With the recent application of novel sequencing technology in APA profiling, an ever-increasing number of APA genes/sites have been identified. However, the phenotypic relevance of most of these APA isoforms remains elusive, which is largely due to the lack of a convenient genetics tool for APA interference. To address this issue, herein, an efficient method is developed based on the CRISPR-dCas13 system, termed as CRISPR-iPAS. Out of eight different dCas13 proteins, Porphyromonas gulae (Pgu) dCas13b, is identified as the most effective one in blocking the usage of the polyadenylation site (PAS). With guide RNAs targeting at core regulatory elements, dPguCas13b enabled APA regulation of endogenous genes with different APA types, including tandem 3′UTR, alternative terminal exon, as well as intronic PAS. Finally, we demonstrated that the proposed APA perturbation tool could be used to investigate the functional relevance of APA isoforms.     

选择性多聚腺苷酸化的生物学效应及其调控机制研究进展

真核生物基因的前体mRNA(pre-mRNA)及一些lncRNA在成熟过程中其3'端会发生剪切和多聚腺苷酸化反应(cleavage and polyadenylation, C/P),C/P的发生需要多聚腺苷酸化信号(polyadenylation signal, PAS)的存在。选择性多聚腺苷酸化(alternative cleavage and polyadenylation, APA)是指具有多个PAS的基因,在其mRNA3'端成熟过程中,由于选择不同的PAS,导致产生出多个3'UTR长度和序列组成不同的转录异构体。3'UTR长度和序列的不同会影响mRNA的稳定性、翻译效率、运输和细胞定位等,因此APA是真核生物的一个重要转录后调控方式。近年来,对大量动物、植物及酵母的基因组测序分析发现,APA在真核生物广泛存在,针对APA的生物学效应和调控机制开展了一系列研究。目前已鉴定出许多APA调控的顺式调控元件和反式作用因子。本文重点介绍了APA生物学效应和调控机制的最新研究进展,并探讨了未来APA调控的研究方向。     

Use of minigene systems to dissect alternative splicing elements

Pre-mRNA splicing is an essential step for gene expression in higher eukaryotes. The splicing efficiency of individual exons is determined by multiple features involving gene architecture, a variety of cis-acting elements within the exons and flanking introns, and interactions with components of the basal splicing machinery (called the spliceosome) and auxiliary regulatory factors which transiently co-assemble with the spliceosome. Both alternative and constitutive exons are recognized by multiple weak protein:RNA interactions and different exons differ in the interactions which are determinative for exon usage. Alternative exons are often regulated according to cell-specific patterns and regulation is mediated by specific sets of cis-acting elements and trans-acting factors. Transient expression of minigenes is a commonly used in vivo assay to identify the intrinsic features of a gene that control exon usage, identify specific cis-acting elements that control usage of constitutive and alternative exons, identify cis-acting elements that control cell-specific usage of alternative exons, and once regulatory elements have been identified, to identify the trans-acting factors that bind to these elements and modulate splicing. This chapter describes approaches and strategies for using minigenes to define the cis-acting elements that determine splice site usage and to identify and characterize the trans-acting factors that bind to these elements and regulate alternative splicing.     

Predicting genetic interactions with random walks on biological networks

Background  

Several studies have demonstrated that synthetic lethal genetic interactions between gene mutations provide an indication of functional redundancy between molecular complexes and pathways. These observations help explain the finding that organisms are able to tolerate single gene deletions for a large majority of genes. For example, system-wide gene knockout/knockdown studies in S. cerevisiae and C. elegans revealed non-viable phenotypes for a mere 18% and 10% of the genome, respectively. It has been postulated that the low percentage of essential genes reflects the extensive amount of genetic buffering that occurs within genomes. Consistent with this hypothesis, systematic double-knockout screens in S. cerevisiae and C. elegans show that, on average, 0.5% of tested gene pairs are synthetic sick or synthetic lethal. While knowledge of synthetic lethal interactions provides valuable insight into molecular functionality, testing all combinations of gene pairs represents a daunting task for molecular biologists, as the combinatorial nature of these relationships imposes a large experimental burden. Still, the task of mapping pairwise interactions between genes is essential to discovering functional relationships between molecular complexes and pathways, as they form the basis of genetic robustness. Towards the goal of alleviating the experimental workload, computational techniques that accurately predict genetic interactions can potentially aid in targeting the most likely candidate interactions. Building on previous studies that analyzed properties of network topology to predict genetic interactions, we apply random walks on biological networks to accurately predict pairwise genetic interactions. Furthermore, we incorporate all published non-interactions into our algorithm for measuring the topological relatedness between two genes. We apply our method to S. cerevisiae and C. elegans datasets and, using a decision tree classifier, integrate diverse biological networks and show that our method outperforms established methods.     

Identification of alternative 5'/3' splice sites based on the mechanism of splice site competition

Alternative splicing plays an important role in regulating gene expression. Currently, most efficient methods use expressed sequence tags or microarray analysis for large-scale detection of alternative splicing. However, it is difficult to detect all alternative splice events with them because of their inherent limitations. Previous computational methods for alternative splicing prediction could only predict particular kinds of alternative splice events. Thus, it would be highly desirable to predict alternative 5'/3' splice sites with various splicing levels using genomic sequences alone. Here, we introduce the competition mechanism of splice sites selection into alternative splice site prediction. This approach allows us to predict not only rarely used but also frequently used alternative splice sites. On a dataset extracted from the AltSplice database, our method correctly classified approximately 70% of the splice sites into alternative and constitutive, as well as approximately 80% of the locations of real competitors for alternative splice sites. It outperforms a method which only considers features extracted from the splice sites themselves. Furthermore, this approach can also predict the changes in activation level arising from mutations in flanking cryptic splice sites of a given splice site. Our approach might be useful for studying alternative splicing in both computational and molecular biology.     

CRISPRpas: programmable regulation of alternative polyadenylation by dCas9

Most human protein-coding genes produce alternative polyadenylation (APA) isoforms that differ in 3′ UTR size or, when coupled with splicing, have variable coding sequences. APA is an important layer of gene expression program critical for defining cell identity. Here, by using a catalytically dead Cas9 and coupling its target site with polyadenylation site (PAS), we develop a method, named CRISPRpas, to alter APA isoform abundance. CRISPRpas functions by enhancing proximal PAS usage, whose efficiency is influenced by several factors, including targeting strand of DNA, distance between PAS and target sequence and strength of the PAS. For intronic polyadenylation (IPA), splicing features, such as strengths of 5′ splice site and 3′ splice site, also affect CRISPRpas efficiency. We show modulation of APA of multiple endogenous genes, including IPA of PCF11, a master regulator of APA and gene expression. In sum, CRISPRpas offers a programmable tool for APA regulation that impacts gene expression.     

Species-specific regulation of alternative splicing in the C-terminal region of the p53 tumor suppressor gene

Alternative splicing occurs in the C-terminal region of the p53 tumor suppressor gene between two alternative 3′ splice sites in intron 10. This alternative splicing event has been detected in murine cells, but not in rat or human tissues. In this paper, we have characterized the pattern of p53 alternative splicing in cell lines from five different species. Our results confirm that p53 alternative splicing is species-specific, being detected only in cell lines of rodent origin. Using transient transfection assays, we have established that the rat p53 gene undergoes efficient alternative splicing in both mouse and rat cell lines, thus demonstrating that it has all the necessary cis-acting sequences to be alternatively spliced. In contrast, we were unable to detect any usage of the human alternative 3′ splice site under the same experimental conditions. Thus, the low levels or absence of alternatively spliced p53 mRNA in rat and human cell lines seems to be the result of different mechanisms. Our results support the hypothesis that there are species-specific mechanisms implicated in the regulation of p53 activity.     

Performance analysis of MPI collective operations

Previous studies of application usage show that the performance of collective communications are critical for high-performance computing. Despite active research in the field, both general and feasible solution to the optimization of collective communication problem is still missing. In this paper, we analyze and attempt to improve intra-cluster collective communication in the context of the widely deployed MPI programming paradigm by extending accepted models of point-to-point communication, such as Hockney, LogP/LogGP, and PLogP, to collective operations. We compare the predictions from models against the experimentally gathered data and using these results, construct optimal decision function for broadcast collective. We quantitatively compare the quality of the model-based decision functions to the experimentally-optimal one. Additionally, in this work, we also introduce a new form of an optimized tree-based broadcast algorithm, splitted-binary. Our results show that all of the models can provide useful insights into various aspects of the different algorithms as well as their relative performance. Still, based on our findings, we believe that the complete reliance on models would not yield optimal results. In addition, our experimental results have identified the gap parameter as being the most critical for accurate modeling of both the classical point-to-point-based pipeline and our extensions to fan-out topologies.