Eduardo Eyras

Control of alternative splicing through siRNA-mediated transcriptional gene silencing.

When targeting promoter regions, small interfering RNAs (siRNAs) trigger a previously proposed pathway known as transcriptional gene silencing by promoting heterochromatin formation. Here we show that siRNAs targeting intronic or exonic sequences close to an alternative exon regulate the splicing of that exon. The effect occurred in hepatoma and HeLa cells with siRNA antisense strands designed to enter the silencing pathway, suggesting hybridization with nascent pre-mRNA. Unexpectedly, in HeLa cells the sense strands were also effective, suggesting that an endogenous antisense transcript, detectable in HeLa but not in hepatoma cells, acts as a target. The effect depends on Argonaute-1 and is counterbalanced by factors favoring chromatin opening or transcriptional elongation. The increase in heterochromatin marks (dimethylation at Lys9 and trimethylation at Lys27 of histone H3) at the target site, the need for the heterochromatin-associated protein HP1α and the reduction in RNA polymerase II processivity suggest a mechanism involving the kinetic coupling of transcription and alternative splicing.

DGCR8 HITS-CLIP reveals novel functions for the Microprocessor

The Drosha–DGCR8 complex (Microprocessor) is required for microRNA (miRNA) biogenesis. DGCR8 recognizes the RNA substrate, whereas Drosha functions as the endonuclease. Using high-throughput sequencing and cross-linking immunoprecipitation (HITS-CLIP) we identified RNA targets of DGCR8 in human cells. Unexpectedly, miRNAs were not the most abundant targets. DGCR8-bound RNAs also comprised several hundred mRNAs as well as small nucleolar RNAs (snoRNAs) and long noncoding RNAs. We found that the Microprocessor controlled the abundance of several mRNAs as well as of MALAT1. By contrast, DGCR8-mediated cleavage of snoRNAs was independent of Drosha, suggesting the involvement of DGCR8 in cellular complexes with other endonucleases. Binding of DGCR8 to cassette exons is a new mechanism for regulation of the relative abundance of alternatively spliced isoforms. These data provide insights in the complex role of DGCR8 in controlling the fate of several classes of RNAs.

RBM5, 6, and 10 Differentially Regulate NUMB Alternative Splicing to Control Cancer Cell Proliferation

RBM5, a regulator of alternative splicing of apoptotic genes, and its highly homologous RBM6 and RBM10 are RNA-binding proteins frequently deleted or mutated in lung cancer. We report that RBM5/6 and RBM10 antagonistically regulate the proliferative capacity of cancer cells and display distinct positional effects in alternative splicing regulation. We identify the Notch pathway regulator NUMB as a key target of these factors in the control of cell proliferation. NUMB alternative splicing, which is frequently altered in lung cancer, can regulate colony and xenograft tumor formation, and its modulation recapitulates or antagonizes the effects of RBM5, 6, and 10 in cell colony formation. RBM10 mutations identified in lung cancer cells disrupt NUMB splicing regulation to promote cell growth. Our results reveal a key genetic circuit in the control of cancer cell proliferation.

ASTD: The Alternative Splicing and Transcript Diversity database

The Alternative Splicing and Transcript Diversity database (ASTD) gives access to a vast collection of alternative transcripts that integrate transcription initiation, polyadenylation and splicing variant data. Alternative transcripts are derived from the mapping of transcribed sequences to the complete human, mouse and rat genomes using an extension of the computational pipeline developed for the ASD (Alternative Splicing Database) and ATD (Alternative Transcript Diversity) databases, which are now superseded by ASTD. For the human genome, ASTD identifies splicing variants, transcription initiation variants and polyadenylation variants in 68%, 68% and 62% of the gene set, respectively, consistent with current estimates for transcription variation. Users can access ASTD through a variety of browsing and query tools, including expression state-based queries for the identification of tissue-specific isoforms. Participating laboratories have experimentally validated a subset of ASTD-predicted alternative splice forms and alternative polyadenylation forms that were not previously reported. The ASTD database can be accessed at http://www.ebi.ac.uk/astd.

Nucleosome-Driven Transcription Factor Binding and Gene Regulation

Elucidating the global function of a transcription factor implies the identification of its target genes and genomic binding sites. The role of chromatin in this context is unclear, but the dominant view is that factors bind preferentially to nucleosome-depleted regions identified as DNaseI-hypersensitive sites (DHS). Here we show by ChIP, MNase, and DNaseI assays followed by deep sequencing that the progesterone receptor (PR) requires nucleosomes for optimal binding and function. In breast cancer cells treated with progestins, we identified 25,000 PR binding sites (PRbs). The majority of these sites encompassed several copies of the hexanucleotide TGTYCY, which is highly abundant in the genome. We found that functional PRbs accumulate around progesterone-induced genes, mainly in enhancers. Most of these sites overlap with DHS but exhibit high nucleosome occupancy. Progestin stimulation results in remodeling of these nucleosomes with displacement of histones H1 and H2A/H2B dimers. Our results strongly suggest that nucleosomes are crucial for PR binding and hormonal gene regulation.

Key contribution of CPEB4-mediated translational control to cancer progression

Malignant transformation, invasion and angiogenesis rely on the coordinated reprogramming of gene expression in the cells from which the tumor originated. Although deregulated gene expression has been extensively studied at genomic and epigenetic scales, the contribution of the regulation of mRNA-specific translation to this reprogramming is not well understood. Here we show that cytoplasmic polyadenylation element binding protein 4 (CPEB4), an RNA binding protein that mediates meiotic mRNA cytoplasmic polyadenylation and translation, is overexpressed in pancreatic ductal adenocarcinomas and glioblastomas, where it supports tumor growth, vascularization and invasion. We also show that, in pancreatic tumors, the pro-oncogenic functions of CPEB4 originate in the translational activation of mRNAs that are silenced in normal tissue, including the mRNA of tissue plasminogen activator, a key contributor to pancreatic ductal adenocarcinoma malignancy. Taken together, our results document a key role for post-transcriptional gene regulation in tumor development and describe a detailed mechanism for gene expression reprogramming underlying malignant tumor progression.

Genome-Wide Association between Branch Point Properties and Alternative Splicing

The branch point (BP) is one of the three obligatory signals required for pre-mRNA splicing. In mammals, the degeneracy of the motif combined with the lack of a large set of experimentally verified BPs complicates the task of modeling it in silico, and therefore of predicting the location of natural BPs. Consequently, BPs have been disregarded in a considerable fraction of the genome-wide studies on the regulation of splicing in mammals. We present a new computational approach for mammalian BP prediction. Using sequence conservation and positional bias we obtained a set of motifs with good agreement with U2 snRNA binding stability. Using a Support Vector Machine algorithm, we created a model complemented with polypyrimidine tract features, which considerably improves the prediction accuracy over previously published methods. Applying our algorithm to human introns, we show that BP position is highly dependent on the presence of AG dinucleotides in the 3′ end of introns, with distance to the 3′ splice site and BP strength strongly correlating with alternative splicing. Furthermore, experimental BP mapping for five exons preceded by long AG-dinucleotide exclusion zones revealed that, for a given intron, more than one BP can be chosen throughout the course of splicing. Finally, the comparison between exons of different evolutionary ages and pseudo exons suggests a key role of the BP in the pathway of exon creation in human. Our computational and experimental analyses suggest that BP recognition is more flexible than previously assumed, and it appears highly dependent on the presence of downstream polypyrimidine tracts. The reported association between BP features and the splicing outcome suggests that this, so far disregarded but yet crucial, element buries information that can complement current acceptor site models.

The Microprocessor controls the activity of mammalian retrotransposons

More than half of the human genome is made of transposable elements whose ongoing mobilization is a driving force in genetic diversity; however, little is known about how the host regulates their activity. Here, we show that the Microprocessor (Drosha-DGCR8), which is required for microRNA biogenesis, also recognizes and binds RNAs derived from human long interspersed element 1 (LINE-1), Alu and SVA retrotransposons. Expression analyses demonstrate that cells lacking a functional Microprocessor accumulate LINE-1 mRNA and encoded proteins. Furthermore, we show that structured regions of the LINE-1 mRNA can be cleaved in vitro by Drosha. Additionally, we used a cell culture–based assay to show that the Microprocessor negatively regulates LINE-1 and Alu retrotransposition in vivo. Altogether, these data reveal a new role for the Microprocessor as a post-transcriptional repressor of mammalian retrotransposons and a defender of human genome integrity.

Detection of recurrent alternative splicing switches in tumor samples reveals novel signatures of cancer

The determination of the alternative splicing isoforms expressed in cancer is fundamental for the development of tumor-specific molecular targets for prognosis and therapy, but it is hindered by the heterogeneity of tumors and the variability across patients. We developed a new computational method, robust to biological and technical variability, which identifies significant transcript isoform changes across multiple samples. We applied this method to more than 4000 samples from the The Cancer Genome Atlas project to obtain novel splicing signatures that are predictive for nine different cancer types, and find a specific signature for basal-like breast tumors involving the tumor-driver CTNND1. Additionally, our method identifies 244 isoform switches, for which the change occurs in the most abundant transcript. Some of these switches occur in known tumor drivers, including PPARG, CCND3, RALGDS, MITF, PRDM1, ABI1 and MYH11, for which the switch implies a change in the protein product. Moreover, some of the switches cannot be described with simple splicing events. Surprisingly, isoform switches are independent of somatic mutations, except for the tumor-suppressor FBLN2 and the oncogene MYH11. Our method reveals novel signatures of cancer in terms of transcript isoforms specifically expressed in tumors, providing novel potential molecular targets for prognosis and therapy. Data and software are available at: http://dx.doi.org/10.6084/m9.figshare.1061917 and https://bitbucket.org/regulatorygenomicsupf/iso-ktsp.

Pyicos: a versatile toolkit for the analysis of high-throughput sequencing data

Motivation: High-throughput sequencing (HTS) has revolutionized gene regulation studies and is now fundamental for the detection of protein–DNA and protein–RNA binding, as well as for measuring RNA expression. With increasing variety and sequencing depth of HTS datasets, the need for more flexible and memory-efficient tools to analyse them is growing. Results: We describe Pyicos, a powerful toolkit for the analysis of mapped reads from diverse HTS experiments: ChIP-Seq, either punctuated or broad signals, CLIP-Seq and RNA-Seq. We prove the effectiveness of Pyicos to select for significant signals and show that its accuracy is comparable and sometimes superior to that of methods specifically designed for each particular type of experiment. Pyicos facilitates the analysis of a variety of HTS datatypes through its flexibility and memory efficiency, providing a useful framework for data integration into models of regulatory genomics. Availability: Open-source software, with tutorials and protocol files, is available at http://regulatorygenomics.upf.edu/pyicos or as a Galaxy server at http://regulatorygenomics.upf.edu/galaxy Contact: eduardo.eyras@upf.edu Supplementary Information: Supplementary data are available at Bioinformatics online.