Probing chromatin accessibility with small molecule DNA intercalation and nanopore sequencing

Genome-wide identification of chromatin organization and structure has been generally probed by measuring accessibility of the underlying DNA to nucleases or methyltransferases. These methods either only observe the positioning of a single nucleosome or rely on large enzymes to modify or cleave the DNA. We developed adduct sequencing (Add-seq), a method to probe chromatin accessibility by treating chromatin with the small molecule angelicin, which preferentially intercalates into DNA not bound to core nucleosomes. We show that Nanopore sequencing of the angelicin-modified DNA is possible and allows visualization and analysis of long single molecules with distinct chromatin structure. The angelicin modification can be detected from the Nanopore current signal data using a neural network model trained on unmodified and modified chromatin-free DNA. Applying Add-seq to Saccharomyces cerevisiae nuclei, we identified expected patterns of accessibility around annotated gene loci in yeast. We also identify individual clusters of single molecule reads displaying different chromatin structure at specific yeast loci, which demonstrates heterogeneity in the chromatin structure of the yeast population. Thus, using Add-seq, we are able to profile DNA accessibility in the yeast genome across long molecules.

Full-length transcript alterations in human bronchial epithelial cells with U2AF1 S34F mutations

U2AF1 is one of the most recurrently mutated splicing factors in lung adenocarcinoma and has been shown to cause transcriptome-wide pre-mRNA splicing alterations; however, the full-length altered mRNA isoforms associated with the mutation are largely unknown. To better understand the impact U2AF1 has on full-length isoform fate and function, we conducted high-throughput long-read cDNA sequencing from isogenic human bronchial epithelial cells with and without a U2AF1 S34F mutation. We identified 49,366 multi-exon transcript isoforms, more than half of which did not match GENCODE or short-read-assembled isoforms. We found 198 transcript isoforms with significant expression and usage changes relative to WT, only 68% of which were assembled by short reads. Expression of isoforms from immune-related genes is largely down-regulated in mutant cells and without observed splicing changes. Finally, we reveal that isoforms likely targeted by nonsense-mediated decay are down-regulated in U2AF1 S34F cells, suggesting that isoform changes may alter the translational output of those affected genes. Altogether, our work provides a resource of full-length isoforms associated with U2AF1 S34F in lung cells.

MET exon 14 skipping is overexpressed in an allele-specific manner in lung adenocarcinoma primary samples

MET exon 14 skipping ( METΔ14 ) is a well-characterized oncogene in the Ras-MAPK pathway driving lung adenocarcinoma (LUAD). Previous studies on METΔ14 revealed this aberrantly spliced oncogene is expressed in LUAD primary samples and is associated with heterozygous somatic mutations and deletions near exon 14 splice sites. Upon further examination of DNA and RNA sequencing data from primary samples, we highlight that METΔ14 is overexpressed in an allele-specific manner. These data suggest that dose-dependence of METΔ14 may be critical to oncogenesis.

Systematic assessment of long-read RNA-seq methods for transcript identification and quantification

The Long-read RNA-Seq Genome Annotation Assessment Project (LRGASP) Consortium was formed to evaluate the effectiveness of long-read approaches for transcriptome analysis. The consortium generated over 427 million long-read sequences from cDNA and direct RNA datasets, encompassing human, mouse, and manatee species, using different protocols and sequencing platforms. These data were utilized by developers to address challenges in transcript isoform detection and quantification, as well as de novo transcript isoform identification. The study revealed that libraries with longer, more accurate sequences produce more accurate transcripts than those with increased read depth, whereas greater read depth improved quantification accuracy. In well-annotated genomes, tools based on reference sequences demonstrated the best performance. When aiming to detect rare and novel transcripts or when using reference-free approaches, incorporating additional orthogonal data and replicate samples are advised. This collaborative study offers a benchmark for current practices and provides direction for future method development in transcriptome analysis.

Detecting haplotype-specific transcript variation in long reads with FLAIR2

Background

RNA-Seq has brought forth significant discoveries regarding aberrations in RNA processing, implicating these RNA variants in a variety of diseases. Aberrant splicing and single nucleotide variants in RNA have been demonstrated to alter transcript stability, localization, and function. In particular, the upregulation of ADAR, an enzyme which mediates adenosine-to-inosine editing, has been previously linked to an increase in the invasiveness of lung ADC cells and associated with splicing regulation. Despite the functional importance of studying splicing and SNVs, short read RNA-Seq has limited the community's ability to interrogate both forms of RNA variation simultaneously.

Results

We employed long-read technology to obtain full-length transcript sequences, elucidating cis-effects of variants on splicing changes at a single molecule level. We have developed a computational workflow that augments FLAIR, a tool that calls isoform models expressed in long-read data, to integrate RNA variant calls with the associated isoforms that bear them. We generated nanopore data with high sequence accuracy of H1975 lung adenocarcinoma cells with and without knockdown of ADAR. We applied our workflow to identify key inosine-isoform associations to help clarify the prominence of ADAR in tumorigenesis.

Conclusions

Ultimately, we find that a long-read approach provides valuable insight toward characterizing the relationship between RNA variants and splicing patterns.

An expression-based variant impact phenotyping protocol to predict the impact of gene variants in cell lines

We describe a bioinformatics protocol for eVIP2 (expression-based variant impact phenotyping). eVIP2 can predict a gene variant’s functional impact by comparing gene expression signatures induced by introduction of wild-type versus mutant cDNAs in cell lines. The predicted functional outcomes of the variants include gain-of-function, loss-of-function, change-of-function, or neutral. eVIP2 improves upon eVIP by being applicable to RNA-seq data and providing pathway-level functional predictions for each mutation. Here, we detail how to run eVIP2 on RNA-seq data from two RNF43 variants.

For complete details on the use and execution of this protocol, please refer to Thornton et al. (2021).

•

A user-friendly computational tool to predict variant impact

•

Profiles many variants across multiple genes in a single run

•

Identification of specific signaling pathways affected by each variant

•

Static and interactive visualization of results

Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Cell Painting predicts impact of lung cancer variants

Most variants in most genes across most organisms have an unknown impact on the function of the corresponding gene. This gap in knowledge is especially acute in cancer, where clinical sequencing of tumors now routinely reveals patient-specific variants whose functional impact on the corresponding genes is unknown, impeding clinical utility. Transcriptional profiling was able to systematically distinguish these variants of unknown significance as impactful vs. neutral in an approach called expression-based variant-impact phenotyping. We profiled a set of lung adenocarcinoma-associated somatic variants using Cell Painting, a morphological profiling assay that captures features of cells based on microscopy using six stains of cell and organelle components. Using deep-learning-extracted features from each cell's image, we found that cell morphological profiling (cmVIP) can predict variants' functional impact and, particularly at the single-cell level, reveals biological insights into variants that can be explored at our public online portal. Given its low cost, convenient implementation, and single-cell resolution, cmVIP profiling therefore seems promising as an avenue for using non-gene specific assays to systematically assess the impact of variants, including disease-associated alleles, on gene function.

SRSF1 governs progenitor-specific alternative splicing to maintain adult epithelial tissue homeostasis and renewal

Alternative splicing generates distinct mRNA variants and is essential for development, homeostasis, and renewal. Proteins of the serine/arginine (SR)-rich splicing factor family are major splicing regulators that are broadly required for organ development as well as cell and organism viability. However, how these proteins support adult organ function remains largely unknown. Here, we used the continuously growing mouse incisor as a model to dissect the functions of the prototypical SR family protein SRSF1 during tissue homeostasis and renewal. We identified an SRSF1-governed alternative splicing network that is specifically required for dental proliferation and survival of progenitors but dispensable for the viability of differentiated cells. We also observed a similar progenitor-specific role of SRSF1 in the small intestinal epithelium, indicating a conserved function of SRSF1 across adult epithelial tissues. Thus, our findings define a regulatory mechanism by which SRSF1 specifically controls progenitor-specific alternative splicing events to support adult tissue homeostasis and renewal.

eVIP2: Expression-based variant impact phenotyping to predict the function of gene variants

While advancements in genome sequencing have identified millions of somatic mutations in cancer, their functional impact is poorly understood. We previously developed the expression-based variant impact phenotyping (eVIP) method to use gene expression data to characterize the function of gene variants. The eVIP method uses a decision tree-based algorithm to predict the functional impact of somatic variants by comparing gene expression signatures induced by introduction of wild-type (WT) versus mutant cDNAs in cell lines. The method distinguishes between variants that are gain-of-function, loss-of-function, change-of-function, or neutral. We present eVIP2, software that allows for pathway analysis (eVIP Pathways) and usage with RNA-seq data. To demonstrate the eVIP2 software and approach, we characterized two recurrent frameshift variants in RNF43, a negative regulator of Wnt signaling, frequently mutated in colorectal, gastric, and endometrial cancer. RNF43 WT, RNF43 R117fs, RNF43 G659fs, or GFP control cDNA were overexpressed in HEK293T cells. Analysis with eVIP2 predicted that the frameshift at position 117 was a loss-of-function mutation, as expected. The second frameshift at position 659 has been previously described as a passenger mutation that maintains the RNF43 WT function as a negative regulator of Wnt. Surprisingly, eVIP2 predicted G659fs to be a change-of-function mutation. Additional eVIP Pathways analysis of RNF43 G659fs predicted 10 pathways to be significantly altered, including TNF-α via NFκB signaling, KRAS signaling, and hypoxia, highlighting the benefit of a more comprehensive approach when determining the impact of gene variant function. To validate these predictions, we performed reporter assays and found that each pathway activated by expression of RNF43 G659fs, but not expression of RNF43 WT, was identified as impacted by eVIP2, supporting that RNF43 G659fs is a change-of-function mutation and its effect on the identified pathways. Pathway activation was further validated by Western blot analysis. Lastly, we show primary colon adenocarcinoma patient samples with R117fs and G659fs variants have transcriptional profiles similar to BRAF missense mutations with activated RAS/MAPK signaling, consistent with KRAS signaling pathways being GOF in both variants. The eVIP2 method is an important step towards overcoming the current challenge of variant interpretation in the implementation of precision medicine. eVIP2 is available at https://github.com/BrooksLabUCSC/eVIP2.

Inflammation drives alternative first exon usage to regulate immune genes including a novel iron-regulated isoform of Aim2

Determining the layers of gene regulation within the innate immune response is critical to our understanding of the cellular responses to infection and dysregulation in disease. We identified a conserved mechanism of gene regulation in human and mouse via changes in alternative first exon (AFE) usage following inflammation, resulting in changes to the isoforms produced. Of these AFE events, we identified 95 unannotated transcription start sites in mice using a de novo transcriptome generated by long-read native RNA-sequencing, one of which is in the cytosolic receptor for dsDNA and known inflammatory inducible gene, Aim2. We show that this unannotated AFE isoform of Aim2 is the predominant isoform expressed during inflammation and contains an iron-responsive element in its 5′UTR enabling mRNA translation to be regulated by iron levels. This work highlights the importance of examining alternative isoform changes and translational regulation in the innate immune response and uncovers novel regulatory mechanisms of Aim2.

Reintroduction of the archaic variant of NOVA1 in cortical organoids alters neurodevelopment

The evolutionarily conserved splicing regulator neuro-oncological ventral antigen 1 (NOVA1) plays a key role in neural development and function. NOVA1 also includes a protein-coding difference between the modern human genome and Neanderthal and Denisovan genomes. To investigate the functional importance of an amino acid change in humans, we reintroduced the archaic allele into human induced pluripotent cells using genome editing and then followed their neural development through cortical organoids. This modification promoted slower development and higher surface complexity in cortical organoids with the archaic version of NOVA1 Moreover, levels of synaptic markers and synaptic protein coassociations correlated with altered electrophysiological properties in organoids expressing the archaic variant. Our results suggest that the human-specific substitution in NOVA1, which is exclusive to modern humans since divergence from Neanderthals, may have had functional consequences for our species' evolution.

A Pan-cancer Transcriptome Analysis Reveals Pervasive Regulation through Alternative Promoters

Most human protein-coding genes are regulated by multiple, distinct promoters, suggesting that the choice of promoter is as important as its level of transcriptional activity. However, while a global change in transcription is recognized as a defining feature of cancer, the contribution of alternative promoters still remains largely unexplored. Here, we infer active promoters using RNA-seq data from 18,468 cancer and normal samples, demonstrating that alternative promoters are a major contributor to context-specific regulation of transcription. We find that promoters are deregulated across tissues, cancer types, and patients, affecting known cancer genes and novel candidates. For genes with independently regulated promoters, we demonstrate that promoter activity provides a more accurate predictor of patient survival than gene expression. Our study suggests that a dynamic landscape of active promoters shapes the cancer transcriptome, opening new diagnostic avenues and opportunities to further explore the interplay of regulatory mechanisms with transcriptional aberrations in cancer.

Proteasome inhibitor-induced modulation reveals the spliceosome as a specific therapeutic vulnerability in multiple myeloma

Enhancing the efficacy of proteasome inhibitors (PI) is a central goal in myeloma therapy. We proposed that signaling-level responses after PI may reveal new mechanisms of action that can be therapeutically exploited. Unbiased phosphoproteomics after treatment with the PI carfilzomib surprisingly demonstrates the most prominent phosphorylation changes on splicing related proteins. Spliceosome modulation is invisible to RNA or protein abundance alone. Transcriptome analysis after PI demonstrates broad-scale intron retention, suggestive of spliceosome interference, as well as specific alternative splicing of protein homeostasis machinery components. These findings lead us to evaluate direct spliceosome inhibition in myeloma, which synergizes with carfilzomib and shows potent anti-tumor activity. Functional genomics and exome sequencing further support the spliceosome as a specific vulnerability in myeloma. Our results propose splicing interference as an unrecognized modality of PI mechanism, reveal additional modes of spliceosome modulation, and suggest spliceosome targeting as a promising therapeutic strategy in myeloma.

Full-length transcript characterization of SF3B1 mutation in chronic lymphocytic leukemia reveals downregulation of retained introns

While splicing changes caused by somatic mutations in SF3B1 are known, identifying full-length isoform changes may better elucidate the functional consequences of these mutations. We report nanopore sequencing of full-length cDNA from CLL samples with and without SF3B1 mutation, as well as normal B cell samples, giving a total of 149 million pass reads. We present FLAIR (Full-Length Alternative Isoform analysis of RNA), a computational workflow to identify high-confidence transcripts, perform differential splicing event analysis, and differential isoform analysis. Using nanopore reads, we demonstrate differential 3' splice site changes associated with SF3B1 mutation, agreeing with previous studies. We also observe a strong downregulation of intron retention events associated with SF3B1 mutation. Full-length transcript analysis links multiple alternative splicing events together and allows for better estimates of the abundance of productive versus unproductive isoforms. Our work demonstrates the potential utility of nanopore sequencing for cancer and splicing research.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

Genomic basis for RNA alterations in cancer

Transcript alterations often result from somatic changes in cancer genomes1. Various forms of RNA alterations have been described in cancer, including overexpression2, altered splicing3 and gene fusions4; however, it is difficult to attribute these to underlying genomic changes owing to heterogeneity among patients and tumour types, and the relatively small cohorts of patients for whom samples have been analysed by both transcriptome and whole-genome sequencing. Here we present, to our knowledge, the most comprehensive catalogue of cancer-associated gene alterations to date, obtained by characterizing tumour transcriptomes from 1,188 donors of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)5. Using matched whole-genome sequencing data, we associated several categories of RNA alterations with germline and somatic DNA alterations, and identified probable genetic mechanisms. Somatic copy-number alterations were the major drivers of variations in total gene and allele-specific expression. We identified 649 associations of somatic single-nucleotide variants with gene expression in cis, of which 68.4% involved associations with flanking non-coding regions of the gene. We found 1,900 splicing alterations associated with somatic mutations, including the formation of exons within introns in proximity to Alu elements. In addition, 82% of gene fusions were associated with structural variants, including 75 of a new class, termed 'bridged' fusions, in which a third genomic location bridges two genes. We observed transcriptomic alteration signatures that differ between cancer types and have associations with variations in DNA mutational signatures. This compendium of RNA alterations in the genomic context provides a rich resource for identifying genes and mechanisms that are functionally implicated in cancer.

A Murine Model of Chronic Lymphocytic Leukemia Based on B Cell-Restricted Expression of Sf3b1 Mutation and Atm Deletion

SF3B1 is recurrently mutated in chronic lymphocytic leukemia (CLL), but its role in the pathogenesis of CLL remains elusive. Here, we show that conditional expression of Sf3b1-K700E mutation in mouse B cells disrupts pre-mRNA splicing, alters cell development, and induces a state of cellular senescence. Combination with Atm deletion leads to the overcoming of cellular senescence and the development of CLL-like disease in elderly mice. These CLL-like cells show genome instability and dysregulation of multiple CLL-associated cellular processes, including deregulated B cell receptor signaling, which we also identified in human CLL cases. Notably, human CLLs harboring SF3B1 mutations exhibit altered response to BTK inhibition. Our murine model of CLL thus provides insights into human CLL disease mechanisms and treatment.

Splicing modulation sensitizes chronic lymphocytic leukemia cells to venetoclax by remodeling mitochondrial apoptotic dependencies

The identification of targetable vulnerabilities in the context of therapeutic resistance is a key challenge in cancer treatment. We detected pervasive aberrant splicing as a characteristic feature of chronic lymphocytic leukemia (CLL), irrespective of splicing factor mutation status, which was associated with sensitivity to the spliceosome modulator, E7107. Splicing modulation affected CLL survival pathways, including members of the B cell lymphoma-2 (BCL2) family of proteins, remodeling antiapoptotic dependencies of human and murine CLL cells. E7107 treatment decreased myeloid cell leukemia-1 (MCL1) dependence and increased BCL2 dependence, sensitizing primary human CLL cells and venetoclax-resistant CLL-like cells from an Eμ-TCL1-based adoptive transfer murine model to treatment with the BCL2 inhibitor venetoclax. Our data provide preclinical rationale to support the combination of venetoclax with splicing modulators to reprogram apoptotic dependencies in CLL for treating venetoclax-resistant CLL cases.

An alternative splicing switch in FLNB promotes the mesenchymal cell state in human breast cancer

Alternative splicing of mRNA precursors represents a key gene expression regulatory step and permits the generation of distinct protein products with diverse functions. In a genome-scale expression screen for inducers of the epithelial-to-mesenchymal transition (EMT), we found a striking enrichment of RNA-binding proteins. We validated that QKI and RBFOX1 were necessary and sufficient to induce an intermediate mesenchymal cell state and increased tumorigenicity. Using RNA-seq and eCLIP analysis, we found that QKI and RBFOX1 coordinately regulated the splicing and function of the actin-binding protein FLNB, which plays a causal role in the regulation of EMT. Specifically, the skipping of FLNB exon 30 induced EMT by releasing the FOXC1 transcription factor. Moreover, skipping of FLNB exon 30 is strongly associated with EMT gene signatures in basal-like breast cancer patient samples. These observations identify a specific dysregulation of splicing, which regulates tumor cell plasticity and is frequently observed in human cancer.

As the human body develops, countless cells change from one state into another. Two important cell states are known as epithelial and mesenchymal. Cells in the epithelial state tend to be tightly connected and form barriers, like skin cells. Mesenchymal state cells are loosely organized, move around more and make up connective tissues. Some cells alternate between these states via an epithelial-to-mesenchymal transition (EMT for short) and back again. Without this transition, certain organs would not develop and wounds would not heal. Yet, cancer cells also use this transition to spread to distant sites of the body. Such cancers are often the most aggressive, and therefore the most deadly.

The epithelial-to-mesenchymal transition is dynamically regulated in a reversible manner. For example, the genes for some proteins might only be active in the epithelial state and further reinforce this state by turning on other ‘epithelial genes’. Alternatively, there might be differences in the processing of mRNA molecules – the intermediate molecules between DNA and protein – that result in the production of different proteins in epithelial and mesenchymal cells.

Li, Choi et al. wanted to know which of the thousands of human genes can endow epithelial state cells with mesenchymal characteristics. A better understanding of the switch could help to prevent cancers undergoing an epithelial-to-mesenchymal transition.

From a large-scale experiment in human breast cancer cells, Li, Choi et al. found that a group of proteins that bind and modify mRNA molecules are important for the epithelial-to-mesenchymal transition. Two proteins in particular promoted the transition, most likely by binding to the mRNA of a third protein called FLNB and removing a small piece of it. FLNB normally works to prevent the epithelial-to-mesenchymal transition, but the smaller protein encoded by the shorter mRNA promoted the transition by turning on ‘mesenchymal genes’.

This switching between different FLNB proteins happens in some of the more aggressive breast cancers, which also contain mesenchymal cells. Finding out which FLNB protein is made in a given cancer may provide an indication of its aggressiveness. Also, looking for drugs that can target the mRNA-binding proteins or FLNB may one day lead to new treatments for some of the most aggressive breast cancers.

Abstract 2466: Identifying confidently measured genes in single pediatric cancer patient samples using RNA sequencing

In the UC Santa Cruz Treehouse Childhood Cancer Initiative (treehousegenomics.soe.ucsc.edu), we are exploring the utility of using RNA-Seq analysis of tumor samples from children to identify potential novel therapeutic options for each individual. Within a single RNA-Seq data set, the gene expression measurements are not equally accurate. The identification of activated, druggable pathways requires accurate gene-level expression measurements.

We receive samples from a variety of clinical and research settings, and the quantity and complexity of the available input material and the depth of sequencing differ. These factors inspired us to develop a tool that will allow us to identify accurate measurements in most RNA-Seq samples we receive.

First, we characterized the relationship between depth of sequencing and the accuracy of the gene expression measurement. We analyzed subsets of reads in samples with more than 50 million Uniquely Mapped, Exonic, Non-duplicate (UMEND) reads. UMEND reads typically constitute over 80% of the reads in a high quality experiment with sufficient starting material. We compared gene expression across the subsets of reads to calculate how many UMEND reads are required to produce consistent measurements. We found that, on average, genes expressed at 1-5 TPM in our data require 30 million reads to be accurately measured. For this calculation, we define accuracy as the condition in which 75% of genes are measured to within 25% of the true value.

Secondly, we use these known relationships to identify genes that have been accurately measured in our tumor RNA-Seq samples. For a sample with 15 million UMEND reads, we find that genes expressed above 5 TPM can be accurately measured and are retained. In the first twelve samples analyzed, samples with more than 10 million UMEND reads retained at least 46% of the genes expressed above zero. We exclude as references those samples with fewer than 10 million UMEND reads due to the marked gene loss after thresholding for this group.

Using accurately measured genes allows us to more confidently assess similarity to other samples, identify enriched pathways, and confirm the expression of drug targets and related molecules under consideration. For example, we reconsidered the CDK4 inhibitor Palbociclib in one patient because the expression of RB1, downstream effector required for Palbociclib-mediated tumor cell death, was under our accuracy threshold. Accuracy thresholds can also be used in experiment planning.

Accuracy thresholding allows us to better assess the value of an RNA-Seq data set and, if necessary, identify the subset of genes whose expression can be confidently considered in a clinical setting. Our experience points to the importance of careful quality control in this process.

Citation Format: Holly Beale, Du Linh Lam, John Vivian, Yulia Newton, Avanthi Tayi Shah, Isabel Bjork, Ted Goldstein, Angela N. Brooks, Josh Stuart, Sofie Salama, E. Alejandro Sweet-Cordero, David Haussler1, Olena Morozova. Identifying confidently measured genes in single pediatric cancer patient samples using RNA sequencing [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 2466. doi:10.1158/1538-7445.AM2017-2466

Abstract 389: Integrating diverse transcriptomic alterations to identify cancer-relevant genes

Introduction:

We present a novel method to identify cancer driver genes that jointly examines any number of diverse transcriptomic alterations with the goal to uncover highly recurrent and heterogeneous patterns in 1190 samples across 26 cancer types as part of the PanCancer Analysis of Whole Genomes (PCAWG) of the International Cancer Genome Consortium (ICGC).

Motivation:

Previous pan-cancer genomic studies have focused on the analysis of somatic mutations as the driver of phenotypic changes. Here, we propose a method to integrate a wide variety of RNA and DNA changes to redefine the concept of driver events and account for the transcriptome’s role in tumorigenesis. PTK2 provides a motivating example, since it has many RNA alterations that correlate with patient survival, such as overexpression, exon-skips, and alternative promoter usage.

In our analysis, we integrate an unprecedented amount of various alterations including gene fusions, RNA editing, alternative splicing, expression outliers, alternative promoters, allele specific expression, and somatic mutations. This enables us to also identify mutually exclusive (MutE) and co-occurring (CoO) patterns between different types of alterations within a gene.

Methods:

Our method has 3 main strengths: flexibility to handle any number or type of alteration, sensitivity to different frequencies of alterations so rare events are not lost in the recurrence analysis, and diversity of ranking such that genes with multiple alterations are prioritized. Our method is summarized in two steps:

1) Identify genes that are both recurrently and heterogeneously altered across many samples by calculating a rank-based score for each gene.

2) Identify MutE and CoO patterns between alteration types for the genes identified in the previous step.

To ensure that alterations were comparable, we applied a thresholding model to binarize all alterations for gene-sample pairs, allowing us to account for the properties of the different modalities involved.

Step 1 of our method calculates a score for each gene that takes into account: 1) the number of alterations to a gene across all samples, 2) the rarity of each alteration, and 3) how many types of alterations are observed per gene. The score is then used to rank the genes and top genes are considered for MutE and CoO analyses.

Results:

Our top 100 ranked genes were highly enriched for cancer census genes (adjusted p-value: 2.06e-9), indicating that we identify cancer relevant genes. Our top five ranked cancer census genes were IGF2, ERBB2, RARA, CREBBP, and ARID1A; all of which had at least 4 of 7 possible alterations, showing our scoring method prioritizes genes with diverse alterations. We also found that alternative promoter usage and alternative splicing were highly co-occurring alterations, with PTK2 having the highest co-occurrence between them. In summary, we propose a new method to analyze various RNA disruptions and show it can yield new insights beyond genomic variation.

Citation Format: Natalie R. Davidson, PanCancer Analysis of Whole Genomes 3 (PCAWG-3) for ICGC, Alvis Brazma, Angela N. Brooks, Claudia Calabrese, Nuno A. Fonseca, Jonathan Goke, Yao He, Xueda Hu, Andre Kahles, Kjong-Van Lehmann, Fenglin Liu, Gunnar Rätsch, Siliang Li, Roland F. Schwarz, Mingyu Yang, Zemin Zhang, Fan Zhang, Liangtao Zheng. Integrating diverse transcriptomic alterations to identify cancer-relevant genes [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 389. doi:10.1158/1538-7445.AM2017-389

RBM25 is a global splicing factor promoting inclusion of alternatively spliced exons and is itself regulated by lysine mono-methylation

In eukaryotes, precursor mRNA (pre-mRNA) splicing removes non-coding intron sequences to produce mature mRNA. This removal is controlled in part by RNA-binding proteins that regulate alternative splicing decisions through interactions with the splicing machinery. RNA binding motif protein 25 (RBM25) is a putative splicing factor strongly conserved across eukaryotic lineages. However, the role of RBM25 in global splicing regulation and its cellular functions are unknown. Here we show that RBM25 is required for the viability of multiple human cell lines, suggesting that it could play a key role in pre-mRNA splicing. Indeed, transcriptome-wide analysis of splicing events demonstrated that RBM25 regulates a large fraction of alternatively spliced exons throughout the human genome. Moreover, proteomic analysis indicated that RBM25 interacts with components of the early spliceosome and regulators of alternative splicing. Previously, we identified an RBM25 species that is mono-methylated at lysine 77 (RBM25K77me1), and here we used quantitative mass spectrometry to show that RBM25K77me1 is abundant in multiple human cell lines. We also identified a region of RBM25 spanning Lys-77 that binds with high affinity to serine- and arginine-rich splicing factor 2 (SRSF2), a crucial protein in exon definition, but only when Lys-77 is unmethylated. Together, our findings uncover a pivotal role for RBM25 as an essential regulator of alternative splicing and reveal a new potential mechanism for regulation of pre-mRNA splicing by lysine methylation of a splicing factor.

Transcriptomic Characterization of SF3B1 Mutation Reveals Its Pleiotropic Effects in Chronic Lymphocytic Leukemia

Mutations in SF3B1, which encodes a spliceosome component, are associated with poor outcome in chronic lymphocytic leukemia (CLL), but how these contribute to CLL progression remains poorly understood. We undertook a transcriptomic characterization of primary human CLL cells to identify transcripts and pathways affected by SF3B1 mutation. Splicing alterations, identified in the analysis of bulk cells, were confirmed in single SF3B1-mutated CLL cells and also found in cell lines ectopically expressing mutant SF3B1. SF3B1 mutation was found to dysregulate multiple cellular functions including DNA damage response, telomere maintenance, and Notch signaling (mediated through KLF8 upregulation, increased TERC and TERT expression, or altered splicing of DVL2 transcript, respectively). SF3B1 mutation leads to diverse changes in CLL-related pathways.

High-throughput Phenotyping of Lung Cancer Somatic Mutations

Recent genome sequencing efforts have identified millions of somatic mutations in cancer. However, the functional impact of most variants is poorly understood. Here we characterize 194 somatic mutations identified in primary lung adenocarcinomas. We present an expression-based variant-impact phenotyping (eVIP) method that uses gene expression changes to distinguish impactful from neutral somatic mutations. eVIP identified 69% of mutations analyzed as impactful and 31% as functionally neutral. A subset of the impactful mutations induces xenograft tumor formation in mice and/or confers resistance to cellular EGFR inhibition. Among these impactful variants are rare somatic, clinically actionable variants including EGFR S645C, ARAF S214C and S214F, ERBB2 S418T, and multiple BRAF variants, demonstrating that rare mutations can be functionally important in cancer.

Abstract 4368: High-throughput phenotyping of lung cancer somatic mutations

Recent cancer genome sequencing and analysis has identified millions of somatic mutations in cancer. However, the functional impact of most variants is poorly understood, limiting the use of this genetic knowledge for clinical decision-making. Here we describe a new high-throughput approach, expression-based variant impact phenotyping (eVIP), which uses gene expression changes to infer somatic mutation impact. We generated a lentiviral expression library representing 53 genes and 194 somatic mutations identified in primary lung adenocarcinomas. Next, we introduced this library into A549 lung adenocarcinoma cells and 96 hours later performed gene expression profiling using Luminex-based L1000 profiling. We built a computational pipeline, eVIP, to compare mutant and wild-type expression signatures to infer whether variants were gain-of-function, change-of-function, loss-of-function, or neutral. Overall, eVIP identified 69% of mutations as impactful whereas 31% appeared functionally neutral. A very high rate, 92%, of missense mutations in the KEAP1 and STK11 tumor suppressor genes were found to inactivate or diminish protein function. As a complementary approach, we assessed which mutations are epistatic to EGFR or capable of initiating xenograft tumor formation in vivo. A subset of the impactful mutations identified by eVIP could induce xenograft tumor formation in mice and/or confer resistance to cellular EGFR inhibition. Among these mutations were 20 rare or non-canonical somatic variants in clinically-actionable or -relevant oncogenes including EGFR S645C, ARAF S214C and S214F, ERBB2 S418T, and PIK3CA E600K. eVIP can, in principle, characterize any genetic variant, independent of prior knowledge of gene function. Further application of eVIP should significantly advance the pace of functional characterization of mutations identified from genome sequencing.

Citation Format: Alice H. Berger, Angela N. Brooks, Xiaoyun Wu, Yashaswi Shrestha, Candace Chouinard, Federica Piccioni, Mukta Bagul, Atanas Kamburov, Marcin Imielinski, Larson Hogstrom, Cong Zhu, Xiaoping Yang, Sasha Pantel, Ryo Sakai, Nathan Kaplan, David Root, Rajiv Narayan, Ted Natoli, David Lahr, Itay Tirosh, Pablo Tamayo, Gad Getz, Bang Wong, John Doench, Aravind Subramanian, Todd R. Golub, Matthew Meyerson, Jesse S. Boehm. High-throughput phenotyping of lung cancer somatic mutations. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 4368.

Abstract 669: Compound heterozygous Sf3b1-K700E mutation and Atm deletion in B cells leads to CLL in mice

Unbiased next-generation sequencing using primary samples has identified SF3B1 as among the most frequently mutated genes in chronic lymphocytic leukemia (CLL). These mutations localize to a restricted gene region, with more than 50% at the K700E site. The presence of mutation is associated with poor clinical outcomes. These lines of evidence together emphasize the high priority for gaining an understanding of role of SF3B1 in CLL. However, the mechanistic basis for how mutated SF3B1 impacts CLL remains unknown. Prior transcriptomic studies using primary CLL samples have led to the appreciation of altered RNA splicing induced by this mutated gene, but studies of the functional impact of mutated SF3B1 in primary human samples have been complicated by its variable mutant allele frequency across samples as well as its common co-occurrence with other diverse gene mutations. We therefore generated a mouse line that conditionally expressed heterozygous Sf3b1-K700E mutation in B lineage cells.

We sought to characterize the effects of Sf3b1-K700E on RNA splicing and B cell function in this in vivo model. By RNA-sequencing of splenic B cells from wild-type and mutant mice (n = 3, each), we detected, classified, and quantified 54 differentially spliced transcripts (p10%) using the tool JuncBASE. Consistent with prior findings in human CLL samples, the splice variants in our mouse model were highly enriched with altered selection of 3’ splice sites (49 of 54 events, p<0.001), suggesting that Sf3b1-K700E mutation in the mouse lines cause altered RNA splicing in vivo and functions in a manner faithful to the human leukemia. B cell characterization of the mutant mice revealed that expression of Sf3b1-K700E caused subtle impairments in B cell development and function, and induced a state of cellular senescence based on gene expression and secretome characterization. Mice with heterozygous Sf3b1-K700E did not develop evidence of CLL-associated immunopathologic changes, even with prolonged aging to 24 months. Nevertheless, expression of Sf3b1-K700E with heterozygous deletion of Atm in B cells led to the overcoming of cellular senescence and resulted in the accumulation of clonal B220+CD5+ cells in peripheral blood, marrow, and spleen of elderly mice. We thus dissect the impact of mutated SF3B1 individually and in combination with Atm deletion on paving the path towards B cell leukemogenesis. This novel murine model faithfully recapitulates features of human CLL and has the potential to shed light on mechanisms of cooperation between Atm deletion and SF3B1 mutation in the pathogenesis of CLL.

Citation Format: Lili Wang, Rutendo Gambe, Jean Fan, Angela N. Brooks, Jing Sun, Donna Neuberg, Peter Kharchenko, Matthew Meyerson, Mark D. Fleming, Benjamin L. Ebert, Ruben Carrasco, Catherine J. Wu. Compound heterozygous Sf3b1-K700E mutation and Atm deletion in B cells leads to CLL in mice. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 669.

Comprehensive Molecular Characterization of Papillary Renal-Cell Carcinoma

BACKGROUND

Papillary renal-cell carcinoma, which accounts for 15 to 20% of renal-cell carcinomas, is a heterogeneous disease that consists of various types of renal cancer, including tumors with indolent, multifocal presentation and solitary tumors with an aggressive, highly lethal phenotype. Little is known about the genetic basis of sporadic papillary renal-cell carcinoma, and no effective forms of therapy for advanced disease exist.

METHODS

We performed comprehensive molecular characterization of 161 primary papillary renal-cell carcinomas, using whole-exome sequencing, copy-number analysis, messenger RNA and microRNA sequencing, DNA-methylation analysis, and proteomic analysis.

RESULTS

Type 1 and type 2 papillary renal-cell carcinomas were shown to be different types of renal cancer characterized by specific genetic alterations, with type 2 further classified into three individual subgroups on the basis of molecular differences associated with patient survival. Type 1 tumors were associated with MET alterations, whereas type 2 tumors were characterized by CDKN2A silencing, SETD2 mutations, TFE3 fusions, and increased expression of the NRF2-antioxidant response element (ARE) pathway. A CpG island methylator phenotype (CIMP) was observed in a distinct subgroup of type 2 papillary renal-cell carcinomas that was characterized by poor survival and mutation of the gene encoding fumarate hydratase (FH).

CONCLUSIONS

Type 1 and type 2 papillary renal-cell carcinomas were shown to be clinically and biologically distinct. Alterations in the MET pathway were associated with type 1, and activation of the NRF2-ARE pathway was associated with type 2; CDKN2A loss and CIMP in type 2 conveyed a poor prognosis. Furthermore, type 2 papillary renal-cell carcinoma consisted of at least three subtypes based on molecular and phenotypic features. (Funded by the National Institutes of Health.).

Abstract B2-21: Identification of somatic RNA splicing alterations in human cancers

Recent whole-exome sequencing studies have found that somatic mutations frequently occur in splicing factors across multiple cancer types, supporting the need to systematically and globally characterize splicing alterations across human cancers. Somatic alterations that affect RNA processing of cancer genes may be caused by multiple mechanisms, including mutations in splicing factors, mutations in cis-acting splice sites and other regulatory sequences, alterations in methylation, or alterations in chromatin structure. For a global look at somatic RNA splicing alterations, we are applying the Cancer Outlier Profile Analysis (COPA) approach on ~7,000 cancer transcriptomes from 12 tumor types in The Cancer Genome Atlas (TCGA) and investigating the underlying somatic mutations that cause these splicing alterations using matched whole-exome, whole-genome, and methylation data from these samples.

To identify and quantify alternative splicing in RNA-Seq data, including unannotated splicing events, we have further developed a computational pipeline called JuncBASE. To distinguish between cancer-specific splicing alterations and normal transcriptome variation, we are utilizing ~700 RNA-Seq libraries from healthy individuals from the Genotype-Tissue Expression (GTEx) project.

COPA analysis of outlier splicing events in lung adenocarcinomas and glioblastomas successfully identified known altered splicing events in MET and EGFR, respectively, that are caused by DNA level somatic mutations. The genomic mechanisms for novel somatic splicing events in known cancer genes such as FGFR3, as well as genes previously uncharacterized in cancer genomes, are currently being investigated.

To identify splicing events that may be novel somatic driver alterations, these events have been profiled using RNA-Seq data from the Cancer Cell Line Encyclopedia and are being used as biomarkers to identify genetic vulnerabilities in high-throughput shRNA screens. This work will have a significant impact on the identification of somatic splicing events that contribute to cancer pathogenesis.

Citation Format: Angela N. Brooks, Yawei Ge, Kevin Chau, Samuel S. Freeman, Gordon Saksena, Chandra Sekhar Pedamallu, Matthew Meyerson. Identification of somatic RNA splicing alterations in human cancers. [abstract]. In: Proceedings of the AACR Special Conference on Computational and Systems Biology of Cancer; Feb 8-11 2015; San Francisco, CA. Philadelphia (PA): AACR; Cancer Res 2015;75(22 Suppl 2):Abstract nr B2-21.

Regulation of alternative splicing in Drosophila by 56 RNA binding proteins

Alternative splicing is regulated by RNA binding proteins (RBPs) that recognize pre-mRNA sequence elements and activate or repress adjacent exons. Here, we used RNA interference and RNA-seq to identify splicing events regulated by 56 Drosophila proteins, some previously unknown to regulate splicing. Nearly all proteins affected alternative first exons, suggesting that RBPs play important roles in first exon choice. Half of the splicing events were regulated by multiple proteins, demonstrating extensive combinatorial regulation. We observed that SR and hnRNP proteins tend to act coordinately with each other, not antagonistically. We also identified a cross-regulatory network where splicing regulators affected the splicing of pre-mRNAs encoding other splicing regulators. This large-scale study substantially enhances our understanding of recent models of splicing regulation and provides a resource of thousands of exons that are regulated by 56 diverse RBPs.

Abstract PR08: NF1, MET, and RIT1 mutations are RAS-pathway driver events in lung adenocarcinoma

Lung adenocarcinoma is comprised of distinct mutational subtypes characterized by mutually exclusive oncogenic alterations such as EGFR and KRAS mutations. Identification of these oncogenic events has transformed treatment strategies and drug development for lung adenocarcinoma. However, such mutations have been identified in only 50-60% of lung adenocarcinoma cases in the United States.

We analyzed SNP array, exome, and RNA sequencing data from 230 primary lung adenocarcinomas to identify known and novel mutually exclusive driver events. MET exon 14 skipping was present in 10/230 tumors (4.3%) and only RNA sequencing data allowed robust detection of this event. To identify novel driver events, we performed a focused statistical analysis of the 87/230 tumors (38%) that lacked a mutation in MET or other known oncogenes. Copy number analysis revealed rare, focal amplifications of MET and ERBB2 that were enriched in this subset. Exome analysis revealed significant enrichment of TP53, KEAP1, NF1, and RIT1 mutations in the previously oncogene-negative subset. Of these mutated genes, only NF1 and RIT1 displayed mutual exclusivity with the other known and candidate driver events.

Whereas TP53, KEAP1, and NF1 are well-documented cancer genes, little is known about the role of RIT1 in human cancer. RIT1 is a RAS-subfamily small GTPase with significant homology to HRAS, NRAS, and KRAS. Mutations in RIT1 were present in ∼2% of cases and clustered in the switch II domain. Ectopic expression of mutated RIT1 in murine fibroblasts and human immortalized lung epithelial cells conferred anchorage-independent growth and tumor forming capability. cDNA sequencing identified a p.M90I missense mutation in the NCI-H2110 cell line. Knockdown of RIT1 in NCI-H2110 cells revealed the role of RIT1 in regulation of PI3K and MEK signaling. Consistently, combination PI3K/mTOR/MEK inhibition could inhibit cellular transformation induced by ectopic expression of mutated RIT1, and PI3K/mTOR inhibition suppressed in vivo tumor growth of NCI-H2110 cells.

Recently, RIT1 mutations were described in myeloid malignancies and Noonan syndrome. Our data, together with these reported observations, indicate that RIT1 is a RAS-pathway oncogene in lung adenocarcinoma. In addition, PI3K/mTOR inhibition with or without MEK inhibition should be explored as a therapeutic strategy for RIT1-mutated lung or other RIT1-mutated cancers. This work expands the number of U.S. lung adenocarcinoma tumors with identified driver events from 62% to 75% and nominates these additional alterations as candidate therapeutic targets.

This abstract is also presented as Poster B26.

Citation Format: Alice H. Berger, Angela N. Brooks, Marcin Imielinski, Andrew Cherniack, Fujiko Duke, Nathan Kaplan, Jeremiah Wala, Matthew Meyerson. NF1, MET, and RIT1 mutations are RAS-pathway driver events in lung adenocarcinoma. [abstract]. In: Proceedings of the AACR Special Conference on RAS Oncogenes: From Biology to Therapy; Feb 24-27, 2014; Lake Buena Vista, FL. Philadelphia (PA): AACR; Mol Cancer Res 2014;12(12 Suppl):Abstract nr PR08. doi: 10.1158/1557-3125.RASONC14-PR08

Abstract 1391: SOX2-p63 interaction and genomic co-localization in squamous cell carcinomas

SOX2 is a transcription factor essential for pluripotent stem cells and the development and maintenance of squamous epithelium. We previously reported SOX2 to be an oncogene subject to highly recurrent genomic amplification in squamous cell carcinomas (SCCs). Here, we aim to further characterize SOX2's function in SCC. We demonstrate that SOX2 binds to distinct genomic loci in SCCs than in embryonic stem (ES) cells through comparative ChIP-seq analysis. Through mass-spectrometric analysis following tandem-affinity immunopurification in SCCs we identify SOX2 interacts with another master squamous transcription factor, p63, instead of its partner in ES cells, OCT4. We find that genomic occupancy of SOX2 in SCC overlaps with that of p63 at a large number of loci and this SOX2-p63 coordinate binding is absent in ES cells. We further demonstrate that SOX2 and p63 jointly regulate gene expression, including the oncogene ETV4, which we find essential for SOX2-amplified SCC cell survival. Together, these findings demonstrate that SOX2's actions in SCC differ substantially from its role in pluripotency. The identification of novel SOX2-p63 interaction enables deeper characterization of SOX2's function in SCC and normal squamous epithelial physiology.

Citation Format: Hideo Watanabe, Qiuping Ma, Shouyong Peng, Guillaume Adelmant, Danielle Swain, Wenyu Song, Cameron Fox, Joshua M. Francis, Chandra Sekhar Pedamallu, David S. Deluca, Angela N. Brooks, Jianwen Que, Anil K. Rustgi, Kwok-kin Wong, Keith L. Ligon, X. Shirley Liu, Jarrod A. Marto, Matthew Meyerson, Adam J. Bass. SOX2-p63 interaction and genomic co-localization in squamous cell carcinomas. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 1391. doi:10.1158/1538-7445.AM2014-1391

SOX2 and p63 colocalize at genetic loci in squamous cell carcinomas

The transcription factor SOX2 is an essential regulator of pluripotent stem cells and promotes development and maintenance of squamous epithelia. We previously reported that SOX2 is an oncogene and subject to highly recurrent genomic amplification in squamous cell carcinomas (SCCs). Here, we have further characterized the function of SOX2 in SCC. Using ChIP-seq analysis, we compared SOX2-regulated gene profiles in multiple SCC cell lines to ES cell profiles and determined that SOX2 binds to distinct genomic loci in SCCs. In SCCs, SOX2 preferentially interacts with the transcription factor p63, as opposed to the transcription factor OCT4, which is the preferred SOX2 binding partner in ES cells. SOX2 and p63 exhibited overlapping genomic occupancy at a large number of loci in SCCs; however, coordinate binding of SOX2 and p63 was absent in ES cells. We further demonstrated that SOX2 and p63 jointly regulate gene expression, including the oncogene ETV4, which was essential for SOX2-amplified SCC cell survival. Together, these findings demonstrate that the action of SOX2 in SCC differs substantially from its role in pluripotency. The identification of the SCC-associated interaction between SOX2 and p63 will enable deeper characterization the downstream targets of this interaction in SCC and normal squamous epithelial physiology.

A pan-cancer analysis of transcriptome changes associated with somatic mutations in U2AF1 reveals commonly altered splicing events

Although recurrent somatic mutations in the splicing factor U2AF1 (also known as U2AF35) have been identified in multiple cancer types, the effects of these mutations on the cancer transcriptome have yet to be fully elucidated. Here, we identified splicing alterations associated with U2AF1 mutations across distinct cancers using DNA and RNA sequencing data from The Cancer Genome Atlas (TCGA). Using RNA-Seq data from 182 lung adenocarcinomas and 167 acute myeloid leukemias (AML), in which U2AF1 is somatically mutated in 3-4% of cases, we identified 131 and 369 splicing alterations, respectively, that were significantly associated with U2AF1 mutation. Of these, 30 splicing alterations were statistically significant in both lung adenocarcinoma and AML, including three genes in the Cancer Gene Census, CTNNB1, CHCHD7, and PICALM. Cell line experiments expressing U2AF1 S34F in HeLa cells and in 293T cells provide further support that these altered splicing events are caused by U2AF1 mutation. Consistent with the function of U2AF1 in 3' splice site recognition, we found that S34F/Y mutations cause preferences for CAG over UAG 3' splice site sequences. This report demonstrates consistent effects of U2AF1 mutation on splicing in distinct cancer cell types.

Abstract 3150: Characterizing the effects of somatic mutations in splice factors on the transcriptome

Recent studies have identified recurrent somatic mutations in splice factors, including SF3B1 and U2AF1, in myelodysplastic syndromes, chronic lymphocytic leukemias, breast cancer, and lung adenocarcinoma. In one study, 91 tumor and matched normal whole-exomes were sequenced from patients with CLL and the second most frequently mutated gene was SF3B1. Analysis of 187 whole-exomes in lung adenocarcinoma and matched normals revealed recurrent somatic mutations in the splice factor U2AF1 in 3% of cases. In addition to U2AF1 mutations, recurrent truncating mutations were observed in a putative splice factor, RBM10, in 4% of lung adenocarcinomas. The recurrence of somatic mutations in splice factors in multiple tumor types suggests that alterations in RNA splicing may play an important role in tumorigenesis.

Although these mutations have been identified in cancer genomes, the effect on the cancer transcriptome has yet to be fully explored; therefore, we do not understand how mutations in splice factors can contribute to tumorigenesis. We hypothesize that somatic mutations in splice factors cause alterations in mRNA splicing of important cancer genes.

We are using a computational tool called JuncBASE to identify alternative splicing in tumor transcriptomes from RNA-Seq data. We are then associating somatic mutations in splice factors, found from whole-exome sequencing, with splicing signatures identified through RNA-Seq.

Preliminary results on the identification of splice events associated with SF3B1 mutations in CLL reveal novel isoforms of multiple genes specifically expressed in SF3B1 mutant samples. Additionally, we examined the effects of SF3B1, U2AF1, and RBM10 mutations on tumor transcriptomes using RNA-Seq data from The Cancer Genome Atlas.

Citation Format: Angela N. Brooks, Youzhong Wan, Peter Choi, Rui Jing, David S. DeLuca, Carrie Sougnez, Juliann Chmielecki, Marcin Imielinski, Gad Getz, Catherine J. Wu, Matthew Meyerson. Characterizing the effects of somatic mutations in splice factors on the transcriptome. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 3150. doi:10.1158/1538-7445.AM2013-3150

Identification and experimental validation of splicing regulatory elements in Drosophila melanogaster reveals functionally conserved splicing enhancers in metazoans

RNA sequence elements involved in the regulation of pre-mRNA splicing have previously been identified in vertebrate genomes by computational methods. Here, we apply such approaches to predict splicing regulatory elements in Drosophila melanogaster and compare them with elements previously found in the human, mouse, and pufferfish genomes. We identified 99 putative exonic splicing enhancers (ESEs) and 231 putative intronic splicing enhancers (ISEs) enriched near weak 5' and 3' splice sites of constitutively spliced introns, distinguishing between those found near short and long introns. We found that a significant proportion (58%) of fly enhancer sequences were previously reported in at least one of the vertebrates. Furthermore, 20% of putative fly ESEs were previously identified as ESEs in human, mouse, and pufferfish; while only two fly ISEs, CTCTCT and TTATAA, were identified as ISEs in all three vertebrate species. Several putative enhancer sequences are similar to characterized binding-site motifs for Drosophila and mammalian splicing regulators. To provide additional evidence for the function of putative ISEs, we separately identified 298 intronic hexamers significantly enriched within sequences phylogenetically conserved among 15 insect species. We found that 73 putative ISEs were among those enriched in conserved regions of the D. melanogaster genome. The functions of nine enhancer sequences were verified in a heterologous splicing reporter, demonstrating that these sequences are sufficient to enhance splicing in vivo. Taken together, these data identify a set of predicted positive-acting splicing regulatory motifs in the Drosophila genome and reveal regulatory sequences that are present in distant metazoan genomes.

Identification of functional elements and regulatory circuits by Drosophila modENCODE

To gain insight into how genomic information is translated into cellular and developmental programs, the Drosophila model organism Encyclopedia of DNA Elements (modENCODE) project is comprehensively mapping transcripts, histone modifications, chromosomal proteins, transcription factors, replication proteins and intermediates, and nucleosome properties across a developmental time course and in multiple cell lines. We have generated more than 700 data sets and discovered protein-coding, noncoding, RNA regulatory, replication, and chromatin elements, more than tripling the annotated portion of the Drosophila genome. Correlated activity patterns of these elements reveal a functional regulatory network, which predicts putative new functions for genes, reveals stage- and tissue-specific regulators, and enables gene-expression prediction. Our results provide a foundation for directed experimental and computational studies in Drosophila and related species and also a model for systematic data integration toward comprehensive genomic and functional annotation.

Conservation of an RNA regulatory map between Drosophila and mammals

Alternative splicing is generally controlled by proteins that bind directly to regulatory sequence elements and either activate or repress splicing of adjacent splice sites in a target pre-mRNA. Here, we have combined RNAi and mRNA-seq to identify exons that are regulated by Pasilla (PS), the Drosophila melanogaster ortholog of mammalian NOVA1 and NOVA2. We identified 405 splicing events in 323 genes that are significantly affected upon depletion of ps, many of which were annotated as being constitutively spliced. The sequence regions upstream and within PS-repressed exons and downstream from PS-activated exons are enriched for YCAY repeats, and these are consistent with the location of these motifs near NOVA-regulated exons in mammals. Thus, the RNA regulatory map of PS and NOVA1/2 is highly conserved between insects and mammals despite the fact that the target gene orthologs regulated by PS and NOVA1/2 are almost entirely nonoverlapping. This observation suggests that the regulatory codes of individual RNA binding proteins may be nearly immutable, yet the regulatory modules controlled by these proteins are highly evolvable.

Genome-wide analysis of alternative pre-mRNA splicing and RNA-binding specificities of the Drosophila hnRNP A/B family members

Heterogeneous nuclear ribonucleoproteins (hnRNPs) have been traditionally seen as proteins packaging RNA nonspecifically into ribonucleoprotein particles (RNPs), but evidence suggests specific cellular functions on discrete target pre-mRNAs. Here we report genome-wide analysis of alternative splicing patterns regulated by four Drosophila homologs of the mammalian hnRNP A/B family (hrp36, hrp38, hrp40, and hrp48). Analysis of the global RNA-binding distributions of each protein revealed both small and extensively bound regions on target transcripts. A significant subset of RNAs were bound and regulated by more than one hnRNP protein, revealing a combinatorial network of interactions. In vitro RNA-binding site selection experiments (SELEX) identified distinct binding motif specificities for each protein, which were overrepresented in their respective regulated and bound transcripts. These results indicate that individual heterogeneous ribonucleoproteins have specific affinities for overlapping, but distinct, populations of target pre-mRNAs controlling their patterns of RNA processing.

Maternal nutrient restriction in early pregnancy programs hepatic mRNA expression of growth-related genes and liver size in adult male sheep

The liver is a major metabolic and endocrine organ of critical importance in the regulation of growth and metabolism. Its function is determined by a complex interaction of nutritionally regulated counter-regulatory hormones. The extent to which hepatic endocrine sensitivity can be programed in utero and whether the resultant adaptations persist into adulthood is unknown and was therefore the subject of this study. Young adult male sheep born to mothers that were fed either a control diet (i.e.100% of total live weight-maintenance requirements) throughout gestation or 50% of that intake (i.e. nutrient restricted (NR)) from 0 to 95 days gestation and thereafter 100% of requirements (taking into account increasing fetal mass) were entered into the study. All mothers gave birth normally at term, the singleton offspring were weaned at 16 weeks, and then reared at pasture until 3 years of age when their livers were sampled. NR offspring were of similar birth and body weights at 3 years of age when they had disproportionately smaller livers than controls. The abundance of mRNA for GH, prolactin, and IGF-II receptors, plus hepatocyte growth factor and suppressor of cytokine signaling-3 were all lower in livers of NR offspring. In contrast, the abundance of the mitochondrial protein voltage-dependent anion channel and the pro-apoptotic factor Bax were up regulated relative to controls. In conclusion, maternal nutrient restriction in early gestation results in adult offspring with smaller livers. This may be mediated by alterations in both hepatic mitogenic and apoptotic factors.

Structural basis for double-stranded RNA processing by Dicer

The specialized ribonuclease Dicer initiates RNA interference by cleaving double-stranded RNA (dsRNA) substrates into small fragments about 25 nucleotides in length. In the crystal structure of an intact Dicer enzyme, the PAZ domain, a module that binds the end of dsRNA, is separated from the two catalytic ribonuclease III (RNase III) domains by a flat, positively charged surface. The 65 angstrom distance between the PAZ and RNase III domains matches the length spanned by 25 base pairs of RNA. Thus, Dicer itself is a molecular ruler that recognizes dsRNA and cleaves a specified distance from the helical end.

The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs

Transfer RNAs (tRNAs) and small nucleolar RNAs (snoRNAs) are two of the largest classes of non-protein-coding RNAs. Conventional gene finders that detect protein-coding genes do not find tRNA and snoRNA genes because they lack the codon structure and statistical signatures of protein-coding genes. Previously, we developed tRNAscan-SE, snoscan and snoGPS for the detection of tRNAs, methylation-guide snoRNAs and pseudouridylation-guide snoRNAs, respectively. tRNAscan-SE is routinely applied to completed genomes, resulting in the identification of thousands of tRNA genes. Snoscan has successfully detected methylation-guide snoRNAs in a variety of eukaryotes and archaea, and snoGPS has identified novel pseudouridylation-guide snoRNAs in yeast and mammals. Although these programs have been quite successful at RNA gene detection, their use has been limited by the need to install and configure the software packages on UNIX workstations. Here, we describe online implementations of these RNA detection tools that make these programs accessible to a wider range of research biologists. The tRNAscan-SE, snoscan and snoGPS servers are available at http://lowelab.ucsc.edu/tRNAscan-SE/, http://lowelab.ucsc.edu/snoscan/ and http://lowelab.ucsc.edu/snoGPS/, respectively.

Impact of maternal undernutrition and fetal number on glucocorticoid, growth hormone and insulin-like growth factor receptor mRNA abundance in the ovine fetal kidney

Epidemiological and animal studies strongly indicate that the environment experienced in utero determines, in part, an individual's likelihood of developing cardiovascular disease in later life. This risk has been further linked to impaired kidney function, as a result of compromised development during fetal life. The present study therefore examined the influence of maternal nutrient restriction (NR), targeted at specific periods of kidney development during early to mid gestation, on the mRNA abundance of receptors for glucocorticoid (GCR), growth hormone (GHR) and insulin-like growth factors-I (IGF-IR) and -II (IGF-IIR), and the IGF-I and -II ligands. This was undertaken in both singleton and twin fetuses. At conception ewes were randomly allocated to either an adequately fed control group or one of four nutrient-restricted groups that were fed half the control amount from 0 to 30, 31 to 65, 66 to 110 or 0 to 110 days gestation. At 110 days gestation all ewes were humanely euthanased and fetal kidneys and surrounding adipose tissue sampled. There was no effect of NR or fetal number on kidney weight, shape or nephron number, but the surrounding fat mass was increased in singleton fetuses exposed to NR for 110 days. An increase in kidney mRNA abundance with NR only occurred in singleton fetuses where IGF-IR mRNA was enhanced with NR from 66-110 days gestation. In twin fetuses, NR had no effect on mRNA abundance. However, for all genes examined mRNA expression was lower in the kidneys of twin compared with singleton fetuses following NR, and the magnitude of the effect was dependent on the timing of NR. In conclusion, the abundance of mRNA for receptors which regulate fetal kidney development are lower in twin animals compared with singletons following periods of nutrient deficiency. This may impact on later kidney development and function.

Programming of adult cardiovascular function after early maternal undernutrition in sheep

The prenatal nutritional environment influences the subsequent risk of hypertension in adulthood. Animal studies have used, generally, the rat as a model species to illustrate the association between maternal nutrient intake and blood pressure in the resulting adult offspring. No study to date has shown programming of adult cardiovascular function in the sheep through maternal dietary intervention. We therefore fed pregnant sheep to either 100% recommended intake from day 0 of gestation to term [ approximately 147 days gestational age (dGA); controls n = 8] or to 50% recommended intake from day 0 to 95 dGA and thereafter to 100% intake (NR; n = 9). Sheep lambed naturally, offspring were weaned at 16 wk, and the male offspring were reared on pasture until 3 yr of age. At this time, cardiovascular catheters were inserted under halothane anesthesia and sheep were allowed 2-4 days recovery. Basal cardiovascular status and pressor responses to infusion of norepinephrine, angiotensin II, and captopril were then assessed alongside basal plasma concentrations of glucose, cortisol, and leptin. NR sheep were of similar birth weight to controls but at 3 yr of age had higher blood pressure before, but not after, feeding. Peripheral sensitivity to vasoconstrictor infusion was similar between dietary groups, although a reflex bradycardia was not apparent in NR sheep during norepinephrine infusion. Circulating leptin correlated well with fat mass and increased more after vasoconstrictor infusion in NR sheep. In conclusion, early NR has been shown to program aspects of cardiovascular control and adipocyte function in adult sheep.

Environmental influences on the fetus and neonate--timing, mechanisms of action and effects on subsequent adult function

Environmental influences on fetal and neonatal development can affect neural, reproductive, immune and cardiovascular function in adult humans and animals. The effects can be exerted at many different stages of development from before conception to after birth. Effects may even be exerted during a preceding generation. Some known and some possible mechanisms are reviewed. Systems likely to be affected include the brain, hypothalamus, pituitary and adrenal glands and the gonads. The effects may be exerted through altered gene expression at any stage of development or through changes in organ structure or physiology.

The effects of undernutrition, in utero, on reproductive function in adult male and female sheep

The aim of this study was to determine the effects of maternal undernutrition during pregnancy on adult reproductive function in male and female offspring. Groups of ewes were fed rations providing either 100% (High, H) or 50% (Low, L) of estimated metabolisable energy (ME) requirements for pregnancy, from mating until day 95 of gestation, and thereafter were conventionally managed. At 20 months of age, LH and FSH profiles, and LH responses to exogenous GnRH were measured in male and female offspring and, in males, testicular responses to exogenous LH (as measured by testosterone concentrations) were also measured. Undernutrition had no effect on the mean birth weights of lambs of either sex, or on testicular size in male animals at either 6 weeks or 20 months of age. L males exhibited significantly higher FSH concentrations than H males (P < 0.05) but there were no differences with treatment in FSH profiles in females, basal LH profiles or gonadotrophin responses to GnRH in offspring of either sex, and no difference in basal testosterone concentrations or in the testosterone response to exogenous LH administration in males. Semen quality at 20 months of age was unaffected by pre-natal undernutrition but ovulation rate was significantly reduced in L compared to H female offspring (P < 0.05). It is concluded that pre-natal undernutrition had no effect on male reproductive development and adult function, but reduced ovulation rate in female progeny. This effect was not associated with a change in gonadotrophin profiles or pituitary responsiveness.

Effect of maternal undernutrition on fetal testicular steroidogenesis during the CNS androgen-responsive period in male sheep fetuses

The aim of this study was to determine the effects of maternal undernutrition, applied during physiologically relevant stages of development of the reproductive system, on reproductive development in male sheep fetuses. Groups of ewes (n = 11-19) were fed rations providing either 100% (high; H) or 50% (low; L) of metabolizable energy requirements for live weight maintenance during selected 'windows', bounded by days 0, 30, 50, 65 and 110 after mating. Ewes of control groups (HH (Expts 1 and 2) and HHH (Expt 3)) were fed the H ration from mating until they were killed at day 50 (Expt 1), day 65 (Expt 2) or day 110 (Expt 3) of gestation, whereas ewes of other groups were fed the L ration for the periods days 0-30 of gestation (LH and LHH), days 31-50 or days 31-65 of gestation (HL and HLH), days 65-110 of gestation (HHL), or day 0 to day 50, day 65 or day 110 of gestation (LL and LLL) when the animals were killed. At day 50 of gestation, there was no effect of nutritional treatment on mean fetal mass or fetal testicular mass, but there was increased expression of mRNA for steroidogenic acute regulatory protein (StAR) in the testes of LL animals (P < 0.05) compared with HH controls. Compared with HH animals, the mean plasma testosterone concentrations of LL fetuses tended to be higher, but this result did not reach significance. At day 65 of gestation there were no significant differences between treatments in mean fetal masses, testicular masses, mean plasma testosterone concentrations or StAR mRNA content. At day 110 of gestation, fetal masses in the LLL group were lower (P < 0.01) than those of control fetuses, although no differences in testicular size or fetal plasma testosterone concentrations were recorded. It is concluded that the effects of undernutrition on reproductive development of male sheep fetuses are dependent on the timing of the period of undernutrition.