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.

Effects of tenofovir disoproxil fumarate on intrahepatic viral burden and liver immune microenvironment in patients with chronic hepatitis B.

BACKGROUND: The impact of nucleos(t)ide analogues on intrahepatic viral burden and immune microenvironment in patients with chronic hepatitis B (CHB) is not clear. OBJECTIVE: We aimed to characterise the effects of tenofovir disoproxil fumarate (TDF) on intrahepatic viral burden and the liver immune microenvironment in patients with CHB. DESIGN: Core liver biopsies were collected at baseline and year 3 from patients with CHB with minimally raised serum alanine aminotransferase in a double-blind placebo-controlled trial (NCT01522625). Paired biopsies were analysed by RNA-sequencing (n=119 pairs), a custom multiplex immunofluorescence assay (n=30 pairs), and HBV-targeted long-read DNA sequencing (n=49 pairs). RESULTS: Both non-integrated and integrated HBV DNA were present in all patients at baseline, with >65% having interchromosomal translocations. Treatment significantly reduced the frequency of HBV core+ hepatocytes and intrahepatic (integrated and non-integrated) HBV DNA, but had no effect on HBsAg+ hepatocytes. Clonally expanded integrations were enriched for HBsAg coding regions and showed dysregulation of nearby genes. At baseline, there was significant enrichment of intrahepatic CD8+ T cell proximity to HBV core+ hepatocytes, but not to HBsAg+ cells. The densities of T cells and B cells were significantly reduced by TDF. Transcriptomic analyses found TDF induced widespread downregulation of immune-related genes including inhibitory and regulatory genes. CONCLUSION: TDF significantly reduced intrahepatic integrated and non-integrated HBV DNA, exerting disparate effects on HBV core+ and HBsAg+ cells and on different immune cell subsets. Our data suggest there may be differential cytotoxic T cell-mediated killing of HBV core+ versus HBsAg+ hepatocytes, providing insights for HBV cure strategies.

Author Correction: Nanopore native RNA sequencing of a human poly(A) transcriptome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Nanopore native RNA sequencing of a human poly(A) transcriptome.

High-throughput complementary DNA sequencing technologies have advanced our understanding of transcriptome complexity and regulation. However, these methods lose information contained in biological RNA because the copied reads are often short and modifications are not retained. We address these limitations using a native poly(A) RNA sequencing strategy developed by Oxford Nanopore Technologies. Our study generated 9.9 million aligned sequence reads for the human cell line GM12878, using thirty MinION flow cells at six institutions. These native RNA reads had a median length of 771 bases, and a maximum aligned length of over 21,000 bases. Mitochondrial poly(A) reads provided an internal measure of read-length quality. We combined these long nanopore reads with higher accuracy short-reads and annotated GM12878 promoter regions to identify 33,984 plausible RNA isoforms. We describe strategies for assessing 3' poly(A) tail length, base modifications and transcript haplotypes.

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.

U2AF1 S34F enhances tumorigenic potential of lung cells by exhibiting synergy with KRAS mutation and altering response to environmental stress.

Although U2AF1 S34F is a recurrent splicing factor mutation in lung adenocarcinoma (ADC), U2AF1 S34F alone is insufficient for producing tumors in previous models. Because lung ADCs with U2AF1 S34F frequently have co-occurring KRAS mutations and smoking histories, we hypothesized that tumor-forming potential arises from U2AF1 S34F interacting with oncogenic KRAS and environmental stress. To elucidate the effect of U2AF1 S34F co-occurring with a second mutation, we generated human bronchial epithelial cells (HBEC3kt) with co-occurring U2AF1 S34F and KRAS G12V . Transcriptome analysis revealed that co-occurring U2AF1 S34F and KRAS G12V differentially impacts inflammatory, cell cycle, and KRAS pathways. Subsequent phenotyping found associated suppressed cytokine production, increased proliferation, anchorage-independent growth, and tumors in mouse xenografts. Interestingly, HBEC3kts harboring only U2AF1 S34F display increased splicing in stress granule protein genes and viability in cigarette smoke concentrate. Our results suggest that U2AF1 S34F may potentiate transformation by granting precancerous cells survival advantage in environmental stress, permitting accumulation of additional mutations like KRAS G12V , which synergize with U2AF1 S34F to transform the cell.

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.

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.

On the Function of Trans-Splicing: No Evidence for Widespread Proteome Diversification in Trypanosomes.

A long-standing mystery of genomic/transcriptomic structure involves spliced leader trans-splicing (SLTS), in which short RNA "tags" transcribed from a distinct genomic locus is added near the 5' end of RNA transcripts by the spliceosome. SLTS has been observed in diverse eukaryotes in a phylogenetic pattern implying recurrent independent evolution. This striking convergence suggests important functions for SLTS, however no general novel function is known. Recent findings of frequent alternative SLTS (ALT-TS) suggest that ALT-TS could impart widespread functionality. Here, we tested the hypothesis that ALT-TS diversifies proteomes by comparing splicing patterns in orthologous genes between two deeply diverged trypanosome parasites. We also tested proteome diversification functions of ALT-TS by utilizing ribosome profiling sequence data. Finally, we investigated ALT-TS as a mechanism to regulate the expression of unproductive transcripts. Although our results indicate the functional importance of some cases of trans-splicing, we find no evidence for the hypothesis that proteome diversification is a general function of trans-splicing.