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.

Selective depletion of HBV-infected hepatocytes by class A capsid assembly modulators requires high levels of intrahepatic HBV core protein.

Capsid assembly mediated by hepatitis B virus (HBV) core protein (HBc) is an essential part of the HBV replication cycle, which is the target for different classes of capsid assembly modulators (CAMs). While both CAM-A ("aberrant") and CAM-E ("empty") disrupt nucleocapsid assembly and reduce extracellular HBV DNA, CAM-As can also reduce extracellular HBV surface antigen (HBsAg) by triggering apoptosis of HBV-infected cells in preclinical mouse models. However, there have not been substantial HBsAg declines in chronic hepatitis B (CHB) patients treated with CAM-As to date. To investigate this disconnect, we characterized the antiviral activity of tool CAM compounds in HBV-infected primary human hepatocytes (PHHs), as well as in HBV-infected human liver chimeric mice and mice transduced with adeno-associated virus-HBV. Mechanistic studies in HBV-infected PHH revealed that CAM-A, but not CAM-E, induced a dose-dependent aggregation of HBc in the nucleus which is negatively regulated by the ubiquitin-binding protein p62. We confirmed that CAM-A, but not CAM-E, induced HBc-positive cell death in both mouse models via induction of apoptotic and inflammatory pathways and demonstrated that the degree of HBV-positive cell loss was positively correlated with intrahepatic HBc levels. Importantly, we determined that there is a significantly lower level of HBc per hepatocyte in CHB patient liver biopsies than in either of the HBV mouse models. Taken together, these data confirm that CAM-As have a unique secondary mechanism with the potential to kill HBc-positive hepatocytes. However, this secondary mechanism appears to require higher intrahepatic HBc levels than is typically observed in CHB patients, thereby limiting the therapeutic potential.

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.

Targeted Long-Read Sequencing Reveals Comprehensive Architecture, Burden, and Transcriptional Signatures from Hepatitis B Virus-Associated Integrations and Translocations in Hepatocellular Carcinoma Cell Lines.

Hepatitis B virus (HBV) can integrate into the chromosomes of infected hepatocytes, creating potentially oncogenic lesions that can lead to hepatocellular carcinoma (HCC). However, our current understanding of integrated HBV DNA architecture, burden, and transcriptional activity is incomplete due to technical limitations. A combination of genomics approaches was used to describe HBV integrations and corresponding transcriptional signatures in three HCC cell lines: huH-1, PLC/PRF/5, and Hep3B. To generate high-coverage, long-read sequencing data, a custom panel of HBV-targeting biotinylated oligonucleotide probes was designed. Targeted long-read DNA sequencing captured entire HBV integration events within individual reads, revealing that integrations may include deletions and inversions of viral sequences. Surprisingly, all three HCC cell lines contain integrations that are associated with host chromosomal translocations. In addition, targeted long-read RNA sequencing allowed for the assignment of transcriptional activity to specific integrations and resolved the contribution of overlapping HBV transcripts. HBV transcripts chimeric with host sequences were resolved in their entirety and often included >1,000 bp of host sequence. This study provides the first comprehensive description of HBV integrations and associated transcriptional activity in three commonly utilized HCC-derived cell lines. The application of novel methods sheds new light on the complexity of these integrations, including HBV bidirectional transcription, nested transcripts, silent integrations, and host genomic rearrangements. The observation of multiple HBV-associated chromosomal translocations gives rise to the hypothesis that HBV is a driver of genetic instability and provides a potential new mechanism for HCC development. IMPORTANCE HCC-derived cell lines have served as practical models to study HBV biology for decades. These cell lines harbor multiple HBV integrations and express only HBV surface antigen (HBsAg). To date, an accurate description of the integration burden, architecture, and transcriptional profile of these cell lines has been limited due to technical constraints. We have developed a targeted long-read sequencing assay that reveals the entire architecture of integrations in these cell lines. In addition, we identified five chromosomal translocations with integrated HBV DNA at the interchromosomal junctions. Incorporation of long-read transcriptome sequencing (RNA-Seq) data indicated that many integrations and translocations were transcriptionally silent. The observation of multiple HBV-associated translocations has strong implications regarding the potential mechanisms for the development of HBV-associated HCC.

Prp8 impacts cryptic but not alternative splicing frequency.

Pre-mRNA splicing must occur with extremely high fidelity. Spliceosomes assemble onto pre-mRNA guided by specific sequences (5' splice site, 3' splice site, and branchpoint). When splice sites are mutated, as in many hereditary diseases, the spliceosome can aberrantly select nearby pseudo- or "cryptic" splice sites, often resulting in nonfunctional protein. How the spliceosome distinguishes authentic splice sites from cryptic splice sites is poorly understood. We performed a Caenorhabditis elegans genetic screen to find cellular factors that affect the frequency with which the spliceosome uses cryptic splice sites and identified two alleles in core spliceosome component Prp8 that alter cryptic splicing frequency. Subsequent complementary genetic and structural analyses in yeast implicate these alleles in the stability of the spliceosome's catalytic core. However, despite a clear effect on cryptic splicing, high-throughput mRNA sequencing of these prp-8 mutant C. elegans reveals that overall alternative splicing patterns are relatively unchanged. Our data suggest the spliceosome evolved intrinsic mechanisms to reduce the occurrence of cryptic splicing and that these mechanisms are distinct from those that impact alternative splicing.

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.

Targeted long-read sequencing reveals clonally expanded HBV-associated chromosomal translocations in patients with chronic hepatitis B.

Background & Aims: HBV infects over 257 million people worldwide and is associated with the development of hepatocellular carcinoma (HCC). Integration of HBV DNA into the host genome is likely a key driver of HCC oncogenesis. Here, we utilise targeted long-read sequencing to determine the structure of HBV DNA integrations as well as full isoform information of HBV mRNA with more accurate quantification than traditional next generation sequencing platforms. Methods: DNA and RNA were isolated from fresh frozen liver biopsies collected within the GS-US-174-0149 clinical trial. A pan-genotypic panel of biotinylated oligos was developed to enrich for HBV sequences from sheared genomic DNA (∼7 kb) and full-length cDNA libraries from poly-adenylated RNA. Samples were sequenced on the PacBio long-read platform and analysed using a custom bioinformatic pipeline. Results: HBV-targeted long-read DNA sequencing generated high coverage data spanning entire integrations. Strikingly, in 13 of 42 samples (31%) we were able to detect HBV sequences flanked by 2 different chromosomes, indicating a chromosomal translocation associated with HBV integration. Chromosomal translocations were unique to each biopsy sample, suggesting that each originated randomly, and in some cases had evidence of clonal expansion. Using targeted long-read RNA sequencing, we determined that upwards of 95% of all HBV transcripts in patients who are HBeAg-positive originate from cccDNA. In contrast, patients who are HBeAg-negative expressed mostly HBsAg from integrations. Conclusions: Targeted lso-Seq allowed for accurate quantitation of the HBV transcriptome and assignment of transcripts to either cccDNA or integration origins. The existence of multiple unique HBV-associated inter-chromosomal translocations in non-HCC CHB patient liver biopsies suggests a novel mechanism with mutagenic potential that may contribute to progression to HCC. Lay summary: Fresh frozen liver biopsies from patients infected with HBV were subjected to targeted long-read RNA and DNA sequencing. Long-read RNA sequencing captures entire HBV transcripts in a single read, allowing for resolution of overlapping transcripts from the HBV genome. This resolution allowed us to quantify the burden of transcription from integrations vs. cccDNA origin in individual patients. Patients who were HBeAg-positive had a significantly larger fraction of the HBV transcriptome originating from cccDNA compared with those who were HBeAg-negative. Long-read DNA sequencing captured entire integrated HBV sequences including multiple kilobases of flanking host sequence within single reads. This resolution allowed us to describe integration events flanked by 2 different host chromosomes, indicating that integrated HBV DNA are associated with inter-chromosomal translocations. This may lead to significant transcriptional dysregulation and drive progression to HCC.

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.