Inferring global levels of alternative splicing isoforms using a generative model of microarray data

MOTIVATION: Alternative splicing (AS) is a frequent step in metozoan gene expression whereby the exons of genes are spliced in different combinations to generate multiple isoforms of mature mRNA. AS functions to enrich an organism's proteomic complexity and regulates gene expression. Despite its importance, the mechanisms underlying AS and its regulation are not well understood, especially in the context of global gene expression patterns. We present here an algorithm referred to as the Generative model for the Alternative Splicing Array Platform (GenASAP) that can predict the levels of AS for thousands of exon skipping events using data generated from custom microarrays. GenASAP uses Bayesian learning in an unsupervised probability model to accurately predict AS levels from the microarray data. GenASAP is capable of learning the hybridization profiles of microarray data, while modeling noise processes and missing or aberrant data. GenASAP has been successfully applied to the global discovery and analysis of AS in mammalian cells and tissues. RESULTS: GenASAP was applied to data obtained from a custom microarray designed for the monitoring of 3126 AS events in mouse cells and tissues. The microarray design included probes specific for exon body and junction sequences formed by the splicing of exons. Our results show that GenASAP provides accurate predictions for over one-third of the total events, as verified by independent RT-PCR assays. SUPPLEMENTARY INFORMATION: http://www.psi.toronto.edu/GenASAP.     

Discovery and analysis of evolutionarily conserved intronic splicing regulatory elements

Knowledge of the functional cis-regulatory elements that regulate constitutive and alternative pre-mRNA splicing is fundamental for biology and medicine. Here we undertook a genome-wide comparative genomics approach using available mammalian genomes to identify conserved intronic splicing regulatory elements (ISREs). Our approach yielded 314 ISREs, and insertions of ~70 ISREs between competing splice sites demonstrated that 84% of ISREs altered 5′ and 94% altered 3′ splice site choice in human cells. Consistent with our experiments, comparisons of ISREs to known splicing regulatory elements revealed that 40%–45% of ISREs might have dual roles as exonic splicing silencers. Supporting a role for ISREs in alternative splicing, we found that 30%–50% of ISREs were enriched near alternatively spliced (AS) exons, and included almost all known binding sites of tissue-specific alternative splicing factors. Further, we observed that genes harboring ISRE-proximal exons have biases for tissue expression and molecular functions that are ISRE-specific. Finally, we discovered that for Nova1, neuronal PTB, hnRNP C, and FOX1, the most frequently occurring ISRE proximal to an alternative conserved exon in the splicing factor strongly resembled its own known RNA binding site, suggesting a novel application of ISRE density and the propensity for splicing factors to auto-regulate to associate RNA binding sites to splicing factors. Our results demonstrate that ISREs are crucial building blocks in understanding general and tissue-specific AS regulation and the biological pathways and functions regulated by these AS events.     

Global, comparative analysis of tissue-specific promoter CpG methylation

Understanding cell-type-specific epigenetic codes on a global level is a major challenge after the sequencing of the human genome has been completed. Here we applied methyl-CpG immunoprecipitation (MCIp) to obtain comparative methylation profiles of coding and noncoding genes in three human tissues, testis, brain, and monocytes. Forty-four mainly testis-specific promoters were independently validated using bisulfite sequencing or single-gene MCIp, confirming the results obtained by the MCIp microarray approach. We demonstrate the previously unknown somatic hypermethylation at many CpG-rich, testis-specific gene promoters, in particular in ampliconic areas of the Y chromosome. We also identify a number of miRNA genes showing tissue-specific methylation patterns. The comparison of the obtained tissue methylation profiles with corresponding gene expression data indicates a significant association between tissue-specific promoter methylation and gene expression, not only in CpG-rich promoters. In summary, our study highlights the exceptional epigenetic status of germ-line cells in testis and provides a global insight into tissue-specific DNA methylation patterns.     

Functional coordination of alternative splicing in the mammalian central nervous system

Background

Alternative splicing (AS) functions to expand proteomic complexity and plays numerous important roles in gene regulation. However, the extent to which AS coordinates functions in a cell and tissue type specific manner is not known. Moreover, the sequence code that underlies cell and tissue type specific regulation of AS is poorly understood.

Results

Using quantitative AS microarray profiling, we have identified a large number of widely expressed mouse genes that contain single or coordinated pairs of alternative exons that are spliced in a tissue regulated fashion. The majority of these AS events display differential regulation in central nervous system (CNS) tissues. Approximately half of the corresponding genes have neural specific functions and operate in common processes and interconnected pathways. Differential regulation of AS in the CNS tissues correlates strongly with a set of mostly new motifs that are predominantly located in the intron and constitutive exon sequences neighboring CNS-regulated alternative exons. Different subsets of these motifs are correlated with either increased inclusion or increased exclusion of alternative exons in CNS tissues, relative to the other profiled tissues.

Conclusion

Our findings provide new evidence that specific cellular processes in the mammalian CNS are coordinated at the level of AS, and that a complex splicing code underlies CNS specific AS regulation. This code appears to comprise many new motifs, some of which are located in the constitutive exons neighboring regulated alternative exons. These data provide a basis for understanding the molecular mechanisms by which the tissue specific functions of widely expressed genes are coordinated at the level of AS.     

Tissue-specific splicing of disordered segments that embed binding motifs rewires protein interaction networks

Alternative inclusion of exons increases the functional diversity of proteins. Among alternatively spliced exons, tissue-specific exons play a critical role in maintaining tissue identity. This raises the question of how tissue-specific protein-coding exons influence protein function. Here we investigate the structural, functional, interaction, and evolutionary properties of constitutive, tissue-specific, and other alternative exons in human. We find that tissue-specific protein segments often contain disordered regions, are enriched in posttranslational modification sites, and frequently embed conserved binding motifs. Furthermore, genes containing tissue-specific exons tend to occupy central positions in interaction networks and display distinct interaction partners in the respective tissues, and are enriched in signaling, development, and disease genes. Based on these findings, we propose that tissue-specific inclusion of disordered segments that contain binding motifs rewires interaction networks and signaling pathways. In this way, tissue-specific splicing may contribute to functional versatility of proteins and increases the diversity of interaction networks across tissues.     

Design of recombinant antibody microarrays for membrane protein profiling of cell lysates and tissue extracts

Generating global protein expression profiles, including also membrane proteins, will be crucial for our understanding of biological processes in health and disease. In this study, we have expanded our antibody microarray technology platform and designed the first human recombinant antibody microarray for membrane proteins targeting crude cell lysates and tissue extracts. We have optimized all key technological parameters and successfully developed a setup for extracting, labeling and analyzing non-fractionated membrane proteomes under non-denaturing conditions. Finally, the platform was also extended and shown to be compatible with simultaneous profiling of both membrane proteins and water-soluble proteins.     

Computational comparative analyses of alternative splicing regulation using full-length cDNA of various eukaryotes

We previously reported a computational approach to infer alternative splicing patterns from Mus musculus full-length cDNA clones and microarray data. Although we predicted a large number of unreported splice variants, the general mechanisms regulating alternative splicing were yet unknown. In the present study, we compared alternative exons and constitutive exons in terms of splice-site strength and frequency of potential regulatory sequences. These regulatory features were further compared among five different species: Homo sapiens, M. musculus, Arabidopsis thaliana, Oryza sativa, and Drosophila melanogaster. Solid statistical validations of our comparative analyses indicated that alternative exons have (1) weaker splice sites and (2) more potential regulatory sequences than constitutive exons. Based on our observations, we propose a combinatorial model of alternative splicing mechanisms, which suggests that alternative exons contain weak splice sites regulated alternatively by potential regulatory sequences on the exons.     

RASSF1A对黑色素瘤A375细胞基因网络的影响

RASSF1A(Ras association domain family 1 isoform A)是定位于染色体3p21.3区域的抑瘤基因,编码一个由340个氨基酸残基构成的微管相关蛋白.该基因在包括恶性黑色素瘤在内的多种肿瘤中因启动子高甲基化而表达沉默.本研究建立了RASSF1A稳定表达的恶性黑色素瘤A375细胞系,通过全基因组表达谱基因芯片分析RASSF1A过表达对A375细胞基因表达谱的影响,发现RASSF1A引起184个基因表达上调,26个基因表达下调.通过Realtime RT-PCR对部分差异表达基因进行验证,结果表明与芯片筛选结果一致.RASSF1A影响的差异表达基因功能上归属于细胞生长与增殖、细胞周期、细胞凋亡、细胞间黏附、信号传导等生物过程.采用STRING软件构建了RASSF1A影响的差异表达基因调控网络,结果表明RASSF1A调控的差异表达基因构成一个高连接度的基因网络.其中,炎症细胞因子、转录因子位于网络中央.RASSF1A通过影响炎症细胞因子与转录因子之间的表达,影响A375细胞基因网络,调节黑色素瘤恶性生物学行为.