Gene fusions associated with recurrent amplicons represent a class of passenger aberrations in breast cancer

Transcriptome-wide gene expression outlier analysis pinpoints therapeutic vulnerabilities in colorectal cancer

Multiple strategies are continuously being explored to expand the drug target repertoire in solid tumors. We devised a novel computational workflow for transcriptome-wide gene expression outlier analysis that allows the systematic identification of both overexpression and underexpression events in cancer cells. Here, it was applied to expression values obtained through RNA sequencing in 226 colorectal cancer (CRC) cell lines that were also characterized by whole-exome sequencing and microarray-based DNA methylation profiling. We found cell models displaying an abnormally high or low expression level for 3533 and 965 genes, respectively. Gene expression abnormalities that have been previously associated with clinically relevant features of CRC cell lines were confirmed. Moreover, by integrating multi-omics data, we identified both genetic and epigenetic alternations underlying outlier expression values. Importantly, our atlas of CRC gene expression outliers can guide the discovery of novel drug targets and biomarkers. As a proof of concept, we found that CRC cell lines lacking expression of the MTAP gene are sensitive to treatment with a PRMT5-MTA inhibitor (MRTX1719). Finally, other tumor types may also benefit from this approach.

Breaking paradigms: Long non-coding RNAs forming gene fusions with potential implications in cancer

Long non-coding RNAs (lncRNAs) are non-coding RNAs longer than 200 nucleotides with dynamic regulatory functions. They interact with a wide range of molecules such as DNA, RNA, and proteins to modulate diverse cellular functions through several mechanisms and, if deregulated, they can lead to cancer development and progression. Recently, it has been described that lncRNAs are susceptible to form gene fusions with mRNAs or other lncRNAs, breaking the paradigm of gene fusions consisting mainly of protein-coding genes. However, their biological significance in the tumor phenotype is still uncertain. Therefore, their recent identification opens a new line of research to study their biological role in tumorigenesis, and their potential as biomarkers with clinical relevance or as therapeutic targets. The present study aimed to review the lncRNA fusions identified so far and to know which of them have been associated with a potential function. We address the current challenges to deepen their study as well as the reasons why they represent a future therapeutic window in cancer.

Frequency and functional characterization of fusion genes in squamous cell carcinoma of the lung

INTRODUCTION

In contrast to lung adenocarcinoma (LUAD), targetable genetic alterations are less frequently detected in squamous cell carcinoma of the lung (LUSC). Over the last years, gene fusions have become promising targets in many solid cancers. Here, we analysed a cohort of LUSC, identified recurrent fusion genes and functionally characterised these tumour genomes.

METHODS

A subset of 1608 squamous cell carcinomas of the lung was analysed by means of the FusionPlex® Lung Panel to identify potentially targetable gene fusions using targeted next-generation sequencing. Cases harbouring recurrent gene fusions were further analysed using FISH, Cytoscan HD arrays and cell culture experiments.

RESULTS

We found both, known and novel gene fusions in about 3 % of the cases. Known fusions occurring in lung cancer included ALK::EML4, EGFRvIII, EZR::ROS1 and FGFR3::TACC. We further identified recurrent gene fusions of currently unknown biological function, involving EGFR::VSTM2A and NSD3::FGFR1 and showed that the occurrence of the EGFR::VSTM2A fusion is accompanied by high-level amplification of EGFR. Our analyses further revealed that the genomes of these LUSC patients are chromosomally unstable, which leads us to believe that such non-actionable genomic rearrangements may be a result of "chromosomal chaos" most probably not representing exclusive cancer-driving genes in this cancer entity.

CONCLUSIONS

We emphasise that caution should be taken when novel fusions are found and that the appearance of new gene fusions should always be interpreted in the molecular context of the respective disease.

Systematic discovery of gene fusions in pediatric cancer by integrating RNA-seq and WGS

BACKGROUND

Gene fusions are important cancer drivers in pediatric cancer and their accurate detection is essential for diagnosis and treatment. Clinical decision-making requires high confidence and precision of detection. Recent developments show RNA sequencing (RNA-seq) is promising for genome-wide detection of fusion products but hindered by many false positives that require extensive manual curation and impede discovery of pathogenic fusions.

METHODS

We developed Fusion-sq to overcome existing disadvantages of detecting gene fusions. Fusion-sq integrates and "fuses" evidence from RNA-seq and whole genome sequencing (WGS) using intron-exon gene structure to identify tumor-specific protein coding gene fusions. Fusion-sq was then applied to the data generated from a pediatric pan-cancer cohort of 128 patients by WGS and RNA sequencing.

RESULTS

In a pediatric pan-cancer cohort of 128 patients, we identified 155 high confidence tumor-specific gene fusions and their underlying structural variants (SVs). This includes all clinically relevant fusions known to be present in this cohort (30 patients). Fusion-sq distinguishes healthy-occurring from tumor-specific fusions and resolves fusions in amplified regions and copy number unstable genomes. A high gene fusion burden is associated with copy number instability. We identified 27 potentially pathogenic fusions involving oncogenes or tumor-suppressor genes characterized by underlying SVs, in some cases leading to expression changes indicative of activating or disruptive effects.

CONCLUSIONS

Our results indicate how clinically relevant and potentially pathogenic gene fusions can be identified and their functional effects investigated by combining WGS and RNA-seq. Integrating RNA fusion predictions with underlying SVs advances fusion detection beyond extensive manual filtering. Taken together, we developed a method for identifying candidate gene fusions that is suitable for precision oncology applications. Our method provides multi-omics evidence for assessing the pathogenicity of tumor-specific gene fusions for future clinical decision making.

Advancements in Focal Amplification Detection in Tumor/Liquid Biopsies and Emerging Clinical Applications

Focal amplifications (FAs) are crucial in cancer research due to their significant diagnostic, prognostic, and therapeutic implications. FAs manifest in various forms, such as episomes, double minute chromosomes, and homogeneously staining regions, arising through different mechanisms and mainly contributing to cancer cell heterogeneity, the leading cause of drug resistance in therapy. Numerous wet-lab, mainly FISH, PCR-based assays, next-generation sequencing, and bioinformatics approaches have been set up to detect FAs, unravel the internal structure of amplicons, assess their chromatin compaction status, and investigate the transcriptional landscape associated with their occurrence in cancer cells. Most of them are tailored for tumor samples, even at the single-cell level. Conversely, very limited approaches have been set up to detect FAs in liquid biopsies. This evidence suggests the need to improve these non-invasive investigations for early tumor detection, monitoring disease progression, and evaluating treatment response. Despite the potential therapeutic implications of FAs, such as, for example, the use of HER2-specific compounds for patients with ERBB2 amplification, challenges remain, including developing selective and effective FA-targeting agents and understanding the molecular mechanisms underlying FA maintenance and replication. This review details a state-of-the-art of FA investigation, with a particular focus on liquid biopsies and single-cell approaches in tumor samples, emphasizing their potential to revolutionize the future diagnosis, prognosis, and treatment of cancer patients.

Clinically relevant fusion oncogenes: detection and practical implications

Mechanistically, chimeric genes result from DNA rearrangements and include parts of preexisting normal genes combined at the genomic junction site. Some rearranged genes encode pathological proteins with altered molecular functions. Those which can aberrantly promote carcinogenesis are called fusion oncogenes. Their formation is not a rare event in human cancers, and many of them were documented in numerous study reports and in specific databases. They may have various molecular peculiarities like increased stability of an oncogenic part, self-activation of tyrosine kinase receptor moiety, and altered transcriptional regulation activities. Currently, tens of low molecular mass inhibitors are approved in cancers as the drugs targeting receptor tyrosine kinase (RTK) oncogenic fusion proteins, that is, including ALK, ABL, EGFR, FGFR1-3, NTRK1-3, MET, RET, ROS1 moieties. Therein, the presence of the respective RTK fusion in the cancer genome is the diagnostic biomarker for drug prescription. However, identification of such fusion oncogenes is challenging as the breakpoint may arise in multiple sites within the gene, and the exact fusion partner is generally unknown. There is no gold standard method for RTK fusion detection, and many alternative experimental techniques are employed nowadays to solve this issue. Among them, RNA-seq-based methods offer an advantage of unbiased high-throughput analysis of only transcribed RTK fusion genes, and of simultaneous finding both fusion partners in a single RNA-seq read. Here we focus on current knowledge of biology and clinical aspects of RTK fusion genes, related databases, and laboratory detection methods.

circPVT1 and PVT1/AKT3 show a role in cell proliferation, apoptosis, and tumor subtype-definition in small cell lung cancer

Small cell lung cancer (SCLC) is treated as a homogeneous disease, although the expression of NEUROD1, ASCL1, POU2F3, and YAP1 identifies distinct molecular subtypes. The MYC oncogene, amplified in SCLC, was recently shown to act as a lineage-specific factor to associate subtypes with histological classes. Indeed, MYC-driven SCLCs show a distinct metabolic profile and drug sensitivity. To disentangle their molecular features, we focused on the co-amplified PVT1, frequently overexpressed and originating circular (circRNA) and chimeric RNAs. We analyzed hsa_circ_0001821 (circPVT1) and PVT1/AKT3 (chimPVT1) as examples of such transcripts, respectively, to unveil their tumorigenic contribution to SCLC. In detail, circPVT1 activated a pro-proliferative and anti-apoptotic program when over-expressed in lung cells, and knockdown of chimPVT1 induced a decrease in cell growth and an increase of apoptosis in SCLC in vitro. Moreover, the investigated PVT1 transcripts underlined a functional connection between MYC and YAP1/POU2F3, suggesting that they contribute to the transcriptional landscape associated with MYC amplification. In conclusion, we have uncovered a functional role of circular and chimeric PVT1 transcripts in SCLC; these entities may prove useful as novel biomarkers in MYC-amplified tumors.

Gene regulation on extrachromosomal DNA

Oncogene amplification on extrachromosomal DNA (ecDNA) is prevalent in human cancer and is associated with poor outcomes. Clonal, megabase-sized circular ecDNAs in cancer are distinct from nonclonal, small sub-kilobase-sized DNAs that may arise during normal tissue homeostasis. ecDNAs enable profound changes in gene regulation beyond copy-number gains. An emerging principle of ecDNA regulation is the formation of ecDNA hubs: micrometer-sized nuclear structures of numerous copies of ecDNAs tethered by proteins in spatial proximity. ecDNA hubs enable cooperative and intermolecular sharing of DNA regulatory elements for potent and combinatorial gene activation. The 3D context of ecDNA shapes its gene expression potential, selection for clonal heterogeneity among ecDNAs, distribution through cell division, and reintegration into chromosomes. Technologies for studying gene regulation and structure of ecDNA are starting to answer long-held questions on the distinct rules that govern cancer genes beyond chromosomes.

Mutations in RPS19 may affect ribosome function and biogenesis in Diamond Blackfan Anemia

Ribosomes, the cellular organelles translating the genetic code to proteins, are assemblies of RNA chains and many proteins (RPs) arranged in precise fine-tuned interwoven structures. Mutated ribosomal genes cause ribosomopathies, including Diamond Blackfan Anemia (DBA, a rare heterogeneous red-cell aplasia connected to ribosome malfunction) or failed biogenesis. Combined bioinformatical, structural, and predictive analyses of potential consequences of possibly expressed mutations in eS19, the protein product of the highly mutated RPS19, suggests that mutations in its exposed surface could alter its positioning during assembly and consequently prevent biogenesis, implying a natural selective strategy to avoid malfunctions in ribosome assembly. A search for RPS19 pseudogenes indicated >90% sequence identity with the wild type, hinting at its expression in cases of absent or truncated gene products.

Utilisation of semiconductor sequencing for the detection of predictive biomarkers in glioblastoma

The standard treatment for glioblastoma involves a combination of surgery, radiation and chemotherapy but have limited impact on survival. The exponential increase in targeted agents directed at pivotal oncogenic pathways now provide new therapeutic opportunities for this tumour type. However, lack of comprehensive precision oncology testing at diagnosis means such therapeutic opportunities are potentially overlooked. To investigate the role of semiconductor sequencing for detection of predictive biomarkers in routine glioblastoma samples we have undertaken analysis of test trending data generated by a clinically validated next generation sequencing platform designed to capture actionable genomic variants distributed across 505 genes. Analysis was performed across a cohort of 55 glioblastoma patients. Analysis of trending data has revealed a complex and rich actionable mutational landscape in which 166 actionable mutations were detected across 36 genes linked to 17 off label targeted therapy protocols and 111 clinical trials. The majority of patients harboured three or more actionable mutations affecting key cancer related regulatory networks including the PI3K/AKT/MTOR and RAS/RAF/MEK/MAPK signalling pathways, DNA-damage repair pathways and cell cycle checkpoints. Linkage with immunotherapy and PARP inhibitors was identified in 44% of glioblastoma patients as a consequence of alterations in DNA-damage repair genes. Taken together our data indicates that precision oncology testing utilising semiconductor sequencing can be used to identify a broad therapeutic armamentarium of targeted therapies and immunotherapies that can be potentially employed for the improved clinical management of glioblastoma patients.

Regulatory networks and 5' partner usage of miRNA host gene fusions in breast cancer

Genomic rearrangements in cancer cells can create gene fusions where the juxtaposition of two different genes leads to the production of chimeric proteins or altered gene expression through promoter‐swapping. We have previously shown that fusion transcripts involving microRNA (miRNA) host genes contribute to deregulation of miRNA expression regardless of the protein‐coding potential of these transcripts. Many different genes can also be used as 5′ partners by a miRNA host gene in what we named recurrent miRNA‐convergent fusions. Here, we have explored the properties of 5′ partners in fusion transcripts that involve miRNA hosts in breast tumours from The Cancer Genome Atlas (TCGA). We hypothesised that firstly, 5′ partner genes should belong to pathways and transcriptional programmes that reflect the tumour phenotype and secondly, there should be a selection for fusion events that shape miRNA expression to benefit the tumour cell through the known hallmarks of cancer. We found that the set of 5′ partners in miRNA host fusions is non‐random, with overrepresentation of highly expressed genes in pathways active in cancer including epithelial‐to‐mesenchymal transition, translational regulation and oestrogen signalling. Furthermore, many miRNAs were upregulated in samples with host gene fusions, including established oncogenic miRNAs such as mir‐21 and the mir‐106b~mir‐93~mir‐25 cluster. To the list of mechanisms for deregulation of miRNA expression, we have added fusion transcripts that change the promoter region. We propose that this adds material for genetic selection and tumour evolution in cancer cells and that miRNA host fusions can act as tumour ‘drivers’.

Profiling Receptor Tyrosine Kinase Fusions in Chinese Breast Cancers

Background

Receptor tyrosine kinases (RTKs) are a class of tyrosine kinases that regulate cell-to-cell communication and control a variety of complex biological functions. Dysregulation of RTK signaling partly due to chromosomal rearrangements leads to novel tyrosine kinase fusion oncoproteins that are possibly driver alterations to cancers. Targeting some RTK fusions with specific tyrosine kinases inhibitors (TKIs) is an effective therapeutic strategy across a spectrum of RTK fusion-related cancers. However, there is still a paucity of extensive RTK fusion investigations in breast cancer. This study aims to characterize RTK fusions in Chinese breast cancer patients.

Methods

An in-house DNA sequencing database of 1440 Chinese breast cancer patients with a capture-based panel (520 gene or 108 gene-panel) was thoroughly reviewed. A total of 2,229 samples including 1,045 tissues and 1,184 plasmas were analyzed. RTK fusion was defined as an in-frame fusion with the tyrosine kinase domain of the RTK completely retained. Concomitant mutations were also analyzed and tumor mutational burden (TMB) was calculated. Patients' clinical characteristics were retrieved from case records.

Results

A total of 30 RTK fusion events were identified from 27 breast cancer patients with a prevalence of 1.875%%. FGFR2 fusions were seen the most commonly (n=7), followed by RET (n=5), ROS1 (n=3), NTRK3 (n=3), BRAF (n=2), and NTRK1 (n=2). Other RTK fusions including ALK, EGFR, FGFR1, FGFR3, MET, and NTRK2 were identified in one patient each. A total of 27 unique resultant fusion proteins (22 with a novel partner) were discovered including 19 intrachromosomal rearrangements and 8 interchromosomal ones. Twenty-one fusions had the tyrosine kinase domain in-frame fused with a partner gene and six were juxtaposed with an intergenic space. Among the 27 fusions, FGFR2-WDR11 (E17: intergenic) (n=3) and ETV6-NTRK3 (E5:E15) (n=2) occurred recurrently. Of note, the normalized abundance of RTK fusion (fusion AF/max AF) correlated negatively with TMB (r=-0.48, P=0.017). Patients with TMB < 8 (Mutations/Mb) displayed a higher fusion abundance than those with TMB ≥ 8 (Mutations/Mb) (P=0.025). Moreover, CREBBP mutation only co-occurred with FGFR2 fusion (P=0.012), while NTRK3 fusion and TP53 mutation were mutually exclusive (P=0.019).

Conclusion

This is the first study comprehensively delineating the prevalence and spectrum of RTK fusions in Chinese breast cancers. Further study is ongoing to identify the enriched subpopulation who may benefit from RTK fusion inhibitors.

Toward best practice in cancer mutation detection with whole-genome and whole-exome sequencing

Clinical applications of precision oncology require accurate tests that can distinguish true cancer-specific mutations from errors introduced at each step of next-generation sequencing (NGS). To date, no bulk sequencing study has addressed the effects of cross-site reproducibility, nor the biological, technical and computational factors that influence variant identification. Here we report a systematic interrogation of somatic mutations in paired tumor-normal cell lines to identify factors affecting detection reproducibility and accuracy at six different centers. Using whole-genome sequencing (WGS) and whole-exome sequencing (WES), we evaluated the reproducibility of different sample types with varying input amount and tumor purity, and multiple library construction protocols, followed by processing with nine bioinformatics pipelines. We found that read coverage and callers affected both WGS and WES reproducibility, but WES performance was influenced by insert fragment size, genomic copy content and the global imbalance score (GIV; G > T/C > A). Finally, taking into account library preparation protocol, tumor content, read coverage and bioinformatics processes concomitantly, we recommend actionable practices to improve the reproducibility and accuracy of NGS experiments for cancer mutation detection.

PVT1: A long non-coding RNA recurrently involved in neoplasia-associated fusion transcripts

NGS technologies and bioinformatics tools allow the rapid identification of chimeric transcripts in cancer. More than 40,000 fusions are so far reported in the literature; however, for most of them, the role in oncogenesis is still not fully understood. This is the case for fusions involving the long non-coding RNA (lncRNA) Plasmacytoma variant translocation 1 (PVT1) (8q24.21). This lncRNA displays oncogenic functions in several cancer types interacting with microRNAs and proteins, but the role of PVT1 fusion transcripts is more obscure. These chimeras have been identified in both hematological malignancies and solid tumors, mainly arising from rearrangements and/or amplification of the 8q24 chromosomal region. In this review, we detail the full spectrum of PVT1 fusions in cancer, summarizing current knowledge about their genesis, function, and role as biomarkers.

Differential Expression Analysis of RNA-seq Reads: Overview, Taxonomy, and Tools

Analysis of RNA-sequence (RNA-seq) data is widely used in transcriptomic studies and it has many applications. We review RNA-seq data analysis from RNA-seq reads to the results of differential expression analysis. In addition, we perform a descriptive comparison of tools used in each step of RNA-seq data analysis along with a discussion of important characteristics of these tools. A taxonomy of tools is also provided. A discussion of issues in quality control and visualization of RNA-seq data is also included along with useful tools. Finally, we provide some guidelines for the RNA-seq data analyst, along with research issues and challenges which should be addressed.

Case Study: Systematic Detection and Prioritization of Gene Fusions in Cancer by RNA-Seq: A DIY Toolkit

RNA-seq provides an efficient and sensitive methodology to identify fusion transcripts in cancer tissues. Chimeric reads mapping across two different genes represent potential gene fusions. Various methodologies have been implemented in the detection of gene fusions by RNA-seq. Here we describe a general methodology used in processing and filtering of RNA-seq data, followed by filtering of multiple varieties of artifacts to nominate potentially relevant gene fusions. Functional relevance of gene fusions is assessed based on the predicted domain architecture of the putative fusion proteins.

JAK2/PD-L1/PD-L2 (9p24.1) amplifications in renal cell carcinomas with sarcomatoid transformation: implications for clinical management

Amplifications of JAK2, PD-L1, and PD-L2 at 9p24.1 lead to constitutive expression of PD-L1. This, coupled with JAK2-activation dependent upregulation of PD-L1 and adaptive/induced expression leads to higher tumor PD-L1 expression and immune evasion. Renal tumors were therefore evaluated for 9p24.1 amplifications. A combination of next generation sequencing-based copy number analysis, fluorescence in situ hybridization for JAK2/INSL6 and PD-L1/PD-L2 and immunohistochemistry for phospho-STAT3 (downstream target of JAK2), PD-L1, PD-L2, and PD-1 was performed. In this study we interrogated a "Discovery" cohort of 593 renal tumors, a "Validation" cohort of 398 high-grade renal tumors, The Cancer Genome Atlas (879 cases) and other public datasets (846 cases). 9p24.1 amplifications were significantly enriched in renal tumors with sarcomatoid transformation (5.95%, 15/252) when compared to all histologic subtypes in the combined "Discovery", "Validation" and public datasets (16/2636, 0.6%, p < 0.00001). Specifically, 9p24.1 amplifications amongst sarcomatoid tumors in public datasets, the "Discovery" and "Validation" cohorts were 7.7% (6/92), 15.1% (5/33), and 3.1% (4/127), respectively. Herein, we describe 13 cases and amplification status for these was characterized using next generation sequencing (n = 9) and/or fluorescence in situ hybridization (n = 10). Correlation with PD-L1 immunohistochemistry (n = 10) revealed constitutive expression (mean H-score: 222/300, n = 10). Analysis of outcomes based on PD-L1 expression in tumor cells performed on 282 cases ("Validation" cohort) did not reveal a significant prognostic effect and was likely reflective of advanced disease. A high incidence of constitutive PD-L1 expression in tumor cells in the "Validation" cohort (H-Score ≥250/300) was noted amongst 83 rhabdoid (6%) and 127 sarcomatoid renal tumors (7.1%). This suggests additional mechanisms of constitutive expression other than amplification events. Importantly, two patients with 9p24.1-amplified sarcomatoid renal tumors showed significant response to immunotherapy. In summary, a subset of renal tumors with sarcomatoid transformation exhibits constitutive PD-L1 overexpression and these patients should be evaluated for enhanced response to immunotherapy.

High expression of the vacuole membrane protein 1 (VMP1) is a potential marker of poor prognosis in HER2 positive breast cancer

Background

Fusion genes result from genomic structural changes, which can lead to alterations in gene expression that supports tumor development. The aim of the study was to use fusion genes as a tool to identify new breast cancer (BC) genes with a role in BC progression.

Methods

Fusion genes from breast tumors and BC cell lines were collected from publications. RNA-Seq data from tumors and cell lines were retrieved from databanks and analyzed for fusions with SOAPfuse or the analysis was purchased. Fusion genes identified in both tumors (n = 1724) and cell lines (n = 45) were confirmed by qRT-PCR and sequencing. Their individual genes were ranked by selection criteria that included correlation of their mRNA level with copy number. The expression of the top ranked gene was measured by qRT-PCR in normal tissue and in breast tumors from an exploratory cohort (n = 141) and a validation cohort (n = 277). Expression levels were correlated with clinical and pathological factors as well as the patients’ survival. The results were followed up in BC cohorts from TCGA (n = 818) and METABRIC (n = 2509).

Results

Vacuole membrane protein 1 (VMP1) was the most promising candidate based on specific selection criteria. Its expression was higher in breast tumor tissue than normal tissue (p = 1x10^-4), and its expression was significantly higher in HER2 positive than HER2 negative breast tumors in all four cohorts analyzed. High expression of VMP1 associated with breast cancer specific survival (BCSS) in cohort 1 (hazard ratio (HR) = 2.31, CI 1.27–4.18) and METABRIC (HR = 1.26, CI 1.02–1.57), and also after adjusting for HER2 expression in cohort 1 (HR = 2.03, CI 1.10–3.72). BCSS was not significant in cohort 2 or TCGA cohort, which may be due to differences in treatment regimens.

Conclusions

The results suggest that high VMP1 expression is a potential marker of poor prognosis in HER2 positive BC. Further studies are needed to elucidate how VMP1 could affect pathways supportive of tumorigenesis.

FusionPro, a Versatile Proteogenomic Tool for Identification of Novel Fusion Transcripts and Their Potential Translation Products in Cancer Cells

Fusion proteoforms are translation products derived from gene fusion. Although very rare, the fusion proteoforms play important roles in biomedical science. For example, fusion proteoforms influence the development of tumors by serving as cancer markers or cell cycle regulators. Although numerous studies have reported bioinformatics tools that can predict fusion transcripts, few proteogenomic tools are available that can predict and identify proteoforms. In this study, we develop a versatile proteogenomic tool "FusionPro," which facilitates the identification of fusion transcripts and their potential translatable peptides. FusionPro provides an independent gene fusion prediction module and can build sequence databases for annotated fusion proteoforms. FusionPro shows greater sensitivity than the available fusion finders when analyzing simulated or real RNA sequencing data sets. We use FusionPro to identify 18 fusion junction peptides and three potential fusion-derived peptides by MS/MS-based analysis of leukemia cell lines (Jurkat and K562) and ovarian cancer tissues from the Clinical Proteomic Tumor Analysis Consortium. Among the identified fusion proteins, we molecularly validate two fusion junction isoforms and a translation product of FAM133B:CDK6. Moreover, sequence analysis suggests that the fusion protein participates in the cell cycle progression. In addition, our prediction results indicate that fusion transcripts often have multiple fusion junctions and that these fusion junctions tend to be distributed in a nonrandom pattern at both the chromosome and gene levels. Thus, FusionPro allows users to detect various types of fusion translation products using a transcriptome-informed approach and to gain a comprehensive understanding of the formation and biological roles of fusion proteoforms.

Tyrosine Threonine Kinase Inhibition Eliminates Lung Cancers by Augmenting Apoptosis and Polyploidy

The spindle assembly checkpoint maintains genomic integrity. A key component is tyrosine threonine kinase (TTK, also known as Mps1). TTK antagonism is hypothesized to cause genomic instability and cell death. Interrogating The Cancer Genome Atlas revealed high TTK expression in lung adenocarcinomas and squamous cell cancers versus the normal lung (P < 0.001). This correlated with an unfavorable prognosis in examined lung adenocarcinoma cases (P = 0.007). TTK expression profiles in lung tumors were independently assessed by RNA in situ hybridization. CFI-402257 is a highly selective TTK inhibitor. Its potent antineoplastic effects are reported here against a panel of well-characterized murine and human lung cancer cell lines. Significant antitumorigenic activity followed independent treatments of athymic mice bearing human lung cancer xenografts (6.5 mg/kg, P < 0.05; 8.5 mg/kg, P < 0.01) and immunocompetent mice with syngeneic lung cancers (P < 0.001). CFI-402257 antineoplastic mechanisms were explored. CFI-402257 triggered aneuploidy and apoptotic death of lung cancer cells without changing centrosome number. Reverse phase protein arrays (RPPA) of vehicle versus CFI-402257-treated lung cancers were examined using more than 300 critical growth-regulatory proteins. RPPA bioinformatic analyses discovered CFI-402257 enhanced MAPK signaling, implicating MAPK antagonism in augmenting TTK inhibitory effects. This was independently confirmed using genetic and pharmacologic repression of MAPK that promoted CFI-402257 anticancer actions. TTK antagonism exerted marked antineoplastic effects against lung cancers and MAPK inhibition cooperated. Future work should determine whether CFI-402257 treatment alone or with a MAPK inhibitor is active in the lung cancer clinic.

Molecular Signatures of Fusion Proteins in Cancer

Although gene fusions are recognized as driver mutations in a wide variety of cancers, the general molecular mechanisms underlying oncogenic fusion proteins are insufficiently understood. Here, we employ large-scale data integration and machine learning and (1) identify three functionally distinct subgroups of gene fusions and their molecular signatures; (2) characterize the cellular pathways rewired by fusion events across different cancers; and (3) analyze the relative importance of over 100 structural, functional, and regulatory features of ∼2200 gene fusions. We report subgroups of fusions that likely act as driver mutations and find that gene fusions disproportionately affect pathways regulating cellular shape and movement. Although fusion proteins are similar across different cancer types, they affect cancer type-specific pathways. Key indicators of fusion-forming proteins include high and nontissue specific expression, numerous splice sites, and higher centrality in protein-interaction networks. Together, these findings provide unifying and cancer type-specific trends across diverse oncogenic fusion proteins.

Most gene fusions in cancer are stochastic events

Cancer-associated gene fusions resulting in chimeric proteins or aberrant expression of one or both partner genes are pathogenetically and clinically important in several hematologic malignancies and solid tumors. Since the advent of different types of massively parallel sequencing (MPS), the number of identified gene fusions has increased dramatically, prompting the question whether they all have a biologic impact. By ascertaining the chromosomal locations of 8934 genes involved in 10 861 gene fusions reported in the literature, we here show that there is a highly significant association between gene content of chromosomes and chromosome bands and number of genes involved in fusions. This strongly suggests that a clear majority of gene fusions detected by MPS are stochastic events associated with the number of genes available to participate in fusions and that most reported gene fusions are passengers without any pathogenetic importance.

Discovery of New Fusion Transcripts in a Cohort of Pediatric Solid Cancers at Relapse and Relevance for Personalized Medicine

We hypothetized that pediatric cancers would more likely harbor fusion transcripts. To dissect the complexity of the fusions landscape in recurrent solid pediatric cancers, we conducted a study on 48 patients with different relapsing or resistant malignancies. By analyzing RNA sequencing data with a new in-house pipeline for fusions detection named ChimComp, followed by verification by real-time PCR, we identified and classified the most confident fusion transcripts (FTs) according to their potential biological function and druggability. The majority of FTs were predicted to affect key cancer pathways and described to be involved in oncogenesis. Contrary to previous descriptions, we found no significant correlation between the number of fusions and mutations, emphasizing the particularity to study pre-treated pediatric patients. A considerable proportion of FTs containing tumor suppressor genes was detected, reflecting their importance in pediatric cancers. FTs containing non-receptor tyrosine kinases occurred at low incidence and predominantly in brain tumors. Remarkably, more than 30% of patients presented a potentially druggable high-confidence fusion. In conclusion, we detected new oncogenic FTs in relapsing pediatric cancer patients by establishing a robust pipeline that can be applied to other malignancies, to detect and prioritize experimental validation studies leading to the development of new therapeutic options.

Mesenchymal Chondrosarcoma: a Review with Emphasis on its Fusion-Driven Biology

Mesenchymal chondrosarcoma is a rare but deadly form of chondrosarcoma that typically affects adolescents and young adults. While curative intent is possible for patients with localized disease, few options exist for patients in the unresectable/metastatic setting. Thus, it is imperative to understand the fusion-driven biology of this rare malignant neoplasm so as to lead to the future development of better therapeutics for this disease. This manuscript will briefly review the clinical and pathologic features of mesenchymal chondrosarcoma followed by an appraisal of existing data linked to the fusions, HEY1-NCOA2 and IRF2BP2-CDX1, and the associated downstream pathways.

Genomic and transcriptomic characterization of skull base chordoma

Skull base chordoma is a primary rare malignant bone-origin tumor showing relatively slow growth pattern and locally destructive lesions, which can only be characterized by histologic components. There is no available prognostic or therapeutic biomarker to predict clinical outcome or treatment response and the molecular mechanisms underlying chordoma development still remain unexplored. Therefore, we sought out to identify novel somatic variations that are associated with chordoma progression and potentially employed as therapeutic targets. Thirteen skull base chordomas were subjected for whole-exome and/or whole-transcriptome sequencing. In process, we have identified chromosomal aberration in 1p, 7, 10, 13 and 17q, high frequency of functional germline SNP of the T gene, rs2305089 (P = 0.0038) and several recurrent alterations including MUC4, NBPF1, NPIPB15 mutations and novel gene fusion of SAMD5-SASH1 for the first time in skull base chordoma.

Genomic hotspots but few recurrent fusion genes in breast cancer

The advent of next generation sequencing technologies has boosted the interest in exploring the role of fusion genes in the development and progression of solid tumors. In breast cancer, most of the detected gene fusions seem to be "passenger" events while the presence of recurrent and driver fusions is still under study. We performed RNA sequencing in 55 well-characterized breast cancer samples and 10 adjacent normal breast tissues, complemented by an analysis of SNP array data. We explored the presence of fusion genes and defined their association with breast cancer subtypes, clinical-pathologic characteristics and copy number aberrations. Overall, 370 fusions were detected across the majority of the samples. HER2+ samples had significantly more fusions than triple negative and luminal subtypes. The number of fusions was correlated with histological grade, Ki67 and tumor size. Clusters of fusion genes were observed across the genome and a significant correlation of fusions with copy number aberrations and more specifically amplifications was also revealed. Despite the large number of fusion events, only a few were recurrent, while recurrent individual genes forming fusions with different partners were also detected including the estrogen receptor 1 gene in the previously detected ESR1-CCDC170 fusion. Overall we detected novel gene fusion events while we confirmed previously reported fusions. Genomic hotspots of fusion genes, differences between subtypes and small number of recurrent fusions are the most relevant characteristics of these events in breast cancer. Further investigation is necessary to comprehend the biological significance of these fusions.

Cancer transcriptome profiling at the juncture of clinical translation

Methodological breakthroughs over the past four decades have repeatedly revolutionized transcriptome profiling. Using RNA sequencing (RNA-seq), it has now become possible to sequence and quantify the transcriptional outputs of individual cells or thousands of samples. These transcriptomes provide a link between cellular phenotypes and their molecular underpinnings, such as mutations. In the context of cancer, this link represents an opportunity to dissect the complexity and heterogeneity of tumours and to discover new biomarkers or therapeutic strategies. Here, we review the rationale, methodology and translational impact of transcriptome profiling in cancer.

Frequent miRNA-convergent fusion gene events in breast cancer

Studies of fusion genes have mainly focused on the formation of fusions that result in the production of hybrid proteins or, alternatively, on promoter-switching events that put a gene under the control of aberrant signals. However, gene fusions may also disrupt the transcriptional control of genes that are encoded in introns downstream of the breakpoint. By ignoring structural constraints of the transcribed fusions, we highlight the importance of a largely unexplored function of fusion genes. Here, we show, using breast cancer as an example, that miRNA host genes are specifically enriched in fusion genes and that many different, low-frequency, 5ʹ partners may deregulate the same miRNA irrespective of the coding potential of the fusion transcript. These results indicate that the concept of recurrence, defined by the rate of functionally important aberrations, needs to be revised to encompass convergent fusions that affect a miRNA independently of transcript structure and protein-coding potential.

Fusion gene research traditionally focuses on fusions that result in hybrid proteins or promoter switching events. Here, the authors demonstrate enrichment of fusions in miRNA host genes in breast cancer, highlighting that disparate fusions could have convergent impact on miRNA.

Renal Cell Carcinoma With Chromosome 6p Amplification Including the TFEB Gene: A Novel Mechanism of Tumor Pathogenesis?

Amplification of chromosome 6p has been implicated in aggressive behavior in several cancers, but has not been characterized in renal cell carcinoma (RCC). We identified 9 renal tumors with amplification of chromosome 6p including the TFEB gene, 3 by fluorescence in situ hybridization, and 6 from the Cancer Genome Atlas (TCGA) databases. Patients' ages were 28 to 78 years (median, 61 y). Most tumors were high stage (7/9 pT3a, 2/9 pN1). Using immunohistochemistry, 2/4 were positive for melanocytic markers and cathepsin K. Novel TFEB fusions were reported by TCGA in 2; however, due to a small composition of fusion transcripts compared with full-length transcripts (0.5/174 and 3.3/132 FPKM), we hypothesize that these represent secondary fusions due to amplification. Five specimens (4 TCGA, 1 fluorescence in situ hybridization) had concurrent chromosome 3p copy number loss or VHL deletion. However, these did not resemble clear cell RCC, had negative carbonic anhydrase IX labeling, lacked VHL mutation, and had papillary or unclassified histology (2/4 had gain of chromosome 7 or 17). One tumor each had somatic FH mutation and SMARCB1 mutation. Chromosome 6p amplification including TFEB is a previously unrecognized cytogenetic alteration in RCC, associated with heterogenous tubulopapillary eosinophilic and clear cell histology. The combined constellation of features does not fit cleanly into an existing tumor category (unclassified), most closely resembling papillary or translocation RCC. The tendency for high tumor stage, varied tubulopapillary morphology, and a subset with melanocytic marker positivity suggests the possibility of a unique tumor type, despite some variation in appearance and genetics.

TFEB-VEGFA (6p21.1) co-amplified renal cell carcinoma: a distinct entity with potential implications for clinical management

A subset of renal cell carcinomas shows TFEB overexpression secondary to MALAT1-TFEB gene fusion. As alternate mechanisms of TFEB overexpression are likely to have the same effect, we sought to determine the frequency of amplification of TFEB and the adjacent VEGFA gene at 6p21.1. As patients with metastatic renal cell carcinomas are managed with anti-VEGF therapies, we retrospectively assessed therapeutic response in patients with amplified tumors. Amplification status was analyzed for 875 renal cell carcinomas from our institution, a consultative case and 794 cases from The Cancer Genome Atlas. Cases were classified as having low level (5-10 copies), and high-level amplification (>10 copies), and were further analyzed for adjacent oncogene copy number status (n=6; 3 single-nucleotide polymorphism genomic microarray, 3 The Cancer Genome Atlas) and structural rearrangements (n=1; mate-pair sequencing). These were then reviewed for histopathology, immunophenotype, and response to VEGF-targeted therapy on follow-up. In all, 10/875 (1.1%) institutional cases, 1 consultative case, and 3/794 (0.4%) of The Cancer Genome Atlas cases showed TFEB high-level amplification, while 14/875 (1.6%) cases showed TFEB low-level amplification. All cases had associated VEGFA amplification. This was confirmed with evaluation for copy number changes (n=6). The 6p21.1 high and low-level amplified tumors occurred in adults (mean age: 66), with over half being ≥pT3 (13/25, 52%), and most showed oncocytic, tubulopapillary features and high grade (≥grade 3: 20/22, 91%). These were aggressive tumors with metastasis and death from renal cell carcinoma in 11 (of 24, 46%) cases. Four patients received targeted therapy and had a mean survival of 31 months (range: 17-50) post nephrectomy. In summary, a group of aggressive renal cell carcinomas show genomic amplification of the 6p21.1 region including TFEB and VEGFA genes and share morphologic features. Additional studies are warranted to determine whether these patients respond to anti-VEGF therapy.

Relevance of Fusion Genes in Pediatric Cancers: Toward Precision Medicine

Pediatric cancers differ from adult tumors, especially by their very low mutational rate. Therefore, their etiology could be explained in part by other oncogenic mechanisms such as chromosomal rearrangements, supporting the possible implication of fusion genes in the development of pediatric cancers. Fusion genes result from chromosomal rearrangements leading to the juxtaposition of two genes. Consequently, an abnormal activation of one or both genes is observed. The detection of fusion genes has generated great interest in basic cancer research and in the clinical setting, since these genes can lead to better comprehension of the biological mechanisms of tumorigenesis and they can also be used as therapeutic targets and diagnostic or prognostic biomarkers. In this review, we discuss the molecular mechanisms of fusion genes and their particularities in pediatric cancers, as well as their relevance in murine models and in the clinical setting. We also point out the difficulties encountered in the discovery of fusion genes. Finally, we discuss future perspectives and priorities for finding new innovative therapies in childhood cancer.

Functions of bromodomain-containing proteins and their roles in homeostasis and cancer

Bromodomains (BRDs) are evolutionarily conserved protein-protein interaction modules that are found in a wide range of proteins with diverse catalytic and scaffolding functions and are present in most tissues. BRDs selectively recognize and bind to acetylated Lys residues - particularly in histones - and thereby have important roles in the regulation of gene expression. BRD-containing proteins are frequently dysregulated in cancer, they participate in gene fusions that generate diverse, frequently oncogenic proteins, and many cancer-causing mutations have been mapped to the BRDs themselves. Importantly, BRDs can be targeted by small-molecule inhibitors, which has stimulated many translational research projects that seek to attenuate the aberrant functions of BRD-containing proteins in disease.

Recurrent and pathological gene fusions in breast cancer: current advances in genomic discovery and clinical implications

Gene fusions have long been considered principally as the oncogenic events of hematologic malignancies, but have recently gained wide attention in solid tumors due to several milestone discoveries and the advancement of deep sequencing technologies. With the progress in deep sequencing studies of breast cancer transcriptomes and genomes, the discovery of recurrent and pathological gene fusions in breast cancer is on the focus. Recently, driven by new deep sequencing studies, several recurrent or pathological gene fusions have been identified in breast cancer, including ESR1-CCDC170, SEC16A-NOTCH1, SEC22B-NOTCH2, and ESR1-YAP1 etc. More important, most of these gene fusions are preferentially identified in the more aggressive breast cancers, such as luminal B, basal-like, or endocrine-resistant breast cancer, suggesting recurrent gene fusions as additional key driver events in these tumors other than the known drivers such as the estrogen receptor. In this paper, we have comprehensively summarized the newly identified recurrent or pathological gene fusion events in breast cancer, reviewed the contributions of new genomic and deep sequencing technologies to new fusion discovery and the integrative bioinformatics tools to analyze these data, highlighted the biological relevance and clinical implications of these fusion discoveries, and discussed future directions of gene fusion research in breast cancer.

Landscape of somatic mutations in 560 breast cancer whole-genome sequences

We analysed whole-genome sequences of 560 breast cancers to advance understanding of the driver mutations conferring clonal advantage and the mutational processes generating somatic mutations. We found that 93 protein-coding cancer genes carried probable driver mutations. Some non-coding regions exhibited high mutation frequencies, but most have distinctive structural features probably causing elevated mutation rates and do not contain driver mutations. Mutational signature analysis was extended to genome rearrangements and revealed twelve base substitution and six rearrangement signatures. Three rearrangement signatures, characterized by tandem duplications or deletions, appear associated with defective homologous-recombination-based DNA repair: one with deficient BRCA1 function, another with deficient BRCA1 or BRCA2 function, the cause of the third is unknown. This analysis of all classes of somatic mutation across exons, introns and intergenic regions highlights the repertoire of cancer genes and mutational processes operating, and progresses towards a comprehensive account of the somatic genetic basis of breast cancer.

Discovering and understanding oncogenic gene fusions through data intensive computational approaches

Although gene fusions have been recognized as important drivers of cancer for decades, our understanding of the prevalence and function of gene fusions has been revolutionized by the rise of next-generation sequencing, advances in bioinformatics theory and an increasing capacity for large-scale computational biology. The computational work on gene fusions has been vastly diverse, and the present state of the literature is fragmented. It will be fruitful to merge three camps of gene fusion bioinformatics that appear to rarely cross over: (i) data-intensive computational work characterizing the molecular biology of gene fusions; (ii) development research on fusion detection tools, candidate fusion prioritization algorithms and dedicated fusion databases and (iii) clinical research that seeks to either therapeutically target fusion transcripts and proteins or leverages advances in detection tools to perform large-scale surveys of gene fusion landscapes in specific cancer types. In this review, we unify these different-yet highly complementary and symbiotic-approaches with the view that increased synergy will catalyze advancements in gene fusion identification, characterization and significance evaluation.

Emerging Role of Genomic Rearrangements in Breast Cancer: Applying Knowledge from Other Cancers

Significant advances in our knowledge of cancer genomes are rapidly changing the way we think about tumor biology and the heterogeneity of cancer. Recent successes in genomically-guided treatment approaches accompanied by more sophisticated sequencing techniques have paved the way for deeper investigation into the landscape of genomic rearrangements in cancer. While considerable research on solid tumors has focused on point mutations that directly alter the coding sequence of key genes, far less is known about the role of somatic rearrangements. With many recurring alterations observed across tumor types, there is an obvious need for functional characterization of these genomic biomarkers in order to understand their relevance to tumor biology, therapy, and prognosis. As personalized therapy approaches are turning toward genomic alterations for answers, these biomarkers will become increasingly relevant to the practice of precision medicine. This review discusses the emerging role of genomic rearrangements in breast cancer, with a particular focus on fusion genes. In addition, it raises several key questions on the therapeutic value of such rearrangements and provides a framework to evaluate their significance as predictive and prognostic biomarkers.

Novel kinase fusion transcripts found in endometrial cancer

Recent advances in RNA-sequencing technology have enabled the discovery of gene fusion transcripts in the transcriptome of cancer cells. However, it remains difficult to differentiate the therapeutically targetable fusions from passenger events. We have analyzed RNA-sequencing data and DNA copy number data from 25 endometrial cancer cell lines to identify potential therapeutically targetable fusion transcripts, and have identified 124 high-confidence fusion transcripts, of which 69% are associated with gene amplifications. As targetable fusion candidates, we focused on three in-frame kinase fusion transcripts that retain a kinase domain (CPQ-PRKDC, CAPZA2-MET, and VGLL4-PRKG1). We detected only CPQ-PRKDC fusion transcript in three of 122 primary endometrial cancer tissues. Cell proliferation of the fusion-positive cell line was inhibited by knocking down the expression of wild-type PRKDC but not by blocking the CPQ-PRKDC fusion transcript expression. Quantitative real-time RT-PCR demonstrated that the expression of the CPQ-PRKDC fusion transcript was significantly lower than that of wild-type PRKDC, corresponding to a low transcript allele fraction of this fusion, based on RNA-sequencing read counts. In endometrial cancers, the CPQ-PRKDC fusion transcript may be a passenger aberration related to gene amplification. Our findings suggest that transcript allele fraction is a useful predictor to find bona-fide therapeutic-targetable fusion transcripts.

A novel gammaretroviral shuttle vector insertional mutagenesis screen identifies SHARPIN as a breast cancer metastasis gene and prognostic biomarker

Breast cancer (BC) is the second leading cause of malignancy among U.S. women. Metastasis results in a poor prognosis and increased mortality, but the molecular mechanisms by which metastatic tumors occur are not well understood. Identifying the genes that drive the metastatic process could provide targets for improved therapy and biomarkers to improve BC patient outcomes. Using a forward mutagenesis screen, BC cells mutagenized with a replication-incompetent gammaretroviral vector (γRV) were xenotransplanted into the mammary fat pad of immunodeficient mice. In this approach the vector provirus dysregulates nearby genes, providing a selective advantage to transduced cells to form metastases. Metastatic tumors were analyzed for proviral integration sites to identify nearby candidate metastasis genes. The γRV has a transgene cassette that allows for rescue in bacteria and rapid identification of vector integration sites. Using this approach, we identified the previously described metastasis gene WWTR1 (TAZ), and three other novel candidate metastasis genes including SHARPIN. SHARPIN was independently validated in vivo as a BC metastasis gene. Analysis of patient data showed that SHARPIN expression predicts metastasis-free survival after adjuvant therapy. Our approach has broad potential to identify genes involved in oncogenic processes for BC and other cancers. We show here it can identify both known (WWTR1) and novel (SHARPIN) BC metastasis genes.

INTEGRATE: gene fusion discovery using whole genome and transcriptome data

While next-generation sequencing (NGS) has become the primary technology for discovering gene fusions, we are still faced with the challenge of ensuring that causative mutations are not missed while minimizing false positives. Currently, there are many computational tools that predict structural variations (SV) and gene fusions using whole genome (WGS) and transcriptome sequencing (RNA-seq) data separately. However, as both WGS and RNA-seq have their limitations when used independently, we hypothesize that the orthogonal validation from integrating both data could generate a sensitive and specific approach for detecting high-confidence gene fusion predictions. Fortunately, decreasing NGS costs have resulted in a growing quantity of patients with both data available. Therefore, we developed a gene fusion discovery tool, INTEGRATE, that leverages both RNA-seq and WGS data to reconstruct gene fusion junctions and genomic breakpoints by split-read mapping. To evaluate INTEGRATE, we compared it with eight additional gene fusion discovery tools using the well-characterized breast cell line HCC1395 and peripheral blood lymphocytes derived from the same patient (HCC1395BL). The predictions subsequently underwent a targeted validation leading to the discovery of 131 novel fusions in addition to the seven previously reported fusions. Overall, INTEGRATE only missed six out of the 138 validated fusions and had the highest accuracy of the nine tools evaluated. Additionally, we applied INTEGRATE to 62 breast cancer patients from The Cancer Genome Atlas (TCGA) and found multiple recurrent gene fusions including a subset involving estrogen receptor. Taken together, INTEGRATE is a highly sensitive and accurate tool that is freely available for academic use.

Revealing the Complexity of Breast Cancer by Next Generation Sequencing

Over the last few years the increasing usage of "-omic" platforms, supported by next-generation sequencing, in the analysis of breast cancer samples has tremendously advanced our understanding of the disease. New driver and passenger mutations, rare chromosomal rearrangements and other genomic aberrations identified by whole genome and exome sequencing are providing missing pieces of the genomic architecture of breast cancer. High resolution maps of breast cancer methylomes and sequencing of the miRNA microworld are beginning to paint the epigenomic landscape of the disease. Transcriptomic profiling is giving us a glimpse into the gene regulatory networks that govern the fate of the breast cancer cell. At the same time, integrative analysis of sequencing data confirms an extensive intertumor and intratumor heterogeneity and plasticity in breast cancer arguing for a new approach to the problem. In this review, we report on the latest findings on the molecular characterization of breast cancer using NGS technologies, and we discuss their potential implications for the improvement of existing therapies.

Aberrant RNA splicing in cancer; expression changes and driver mutations of splicing factor genes

Alternative splicing is a widespread process contributing to structural transcript variation and proteome diversity. In cancer, the splicing process is commonly disrupted, resulting in both functional and non-functional end-products. Cancer-specific splicing events are known to contribute to disease progression; however, the dysregulated splicing patterns found on a genome-wide scale have until recently been less well-studied. In this review, we provide an overview of aberrant RNA splicing and its regulation in cancer. We then focus on the executors of the splicing process. Based on a comprehensive catalog of splicing factor encoding genes and analyses of available gene expression and somatic mutation data, we identify cancer-associated patterns of dysregulation. Splicing factor genes are shown to be significantly differentially expressed between cancer and corresponding normal samples, and to have reduced inter-individual expression variation in cancer. Furthermore, we identify enrichment of predicted cancer-critical genes among the splicing factors. In addition to previously described oncogenic splicing factor genes, we propose 24 novel cancer-critical splicing factors predicted from somatic mutations.

Characterization of fusion genes and the significantly expressed fusion isoforms in breast cancer by hybrid sequencing

We developed an innovative hybrid sequencing approach, IDP-fusion, to detect fusion genes, determine fusion sites and identify and quantify fusion isoforms. IDP-fusion is the first method to study gene fusion events by integrating Third Generation Sequencing long reads and Second Generation Sequencing short reads. We applied IDP-fusion to PacBio data and Illumina data from the MCF-7 breast cancer cells. Compared with the existing tools, IDP-fusion detects fusion genes at higher precision and a very low false positive rate. The results show that IDP-fusion will be useful for unraveling the complexity of multiple fusion splices and fusion isoforms within tumorigenesis-relevant fusion genes.

The emerging complexity of gene fusions in cancer

Structural chromosome rearrangements may result in the exchange of coding or regulatory DNA sequences between genes. Many such gene fusions are strong driver mutations in neoplasia and have provided fundamental insights into the disease mechanisms that are involved in tumorigenesis. The close association between the type of gene fusion and the tumour phenotype makes gene fusions ideal for diagnostic purposes, enabling the subclassification of otherwise seemingly identical disease entities. In addition, many gene fusions add important information for risk stratification, and increasing numbers of chimeric proteins encoded by the gene fusions serve as specific targets for treatment, resulting in dramatically improved patient outcomes. In this Timeline article, we describe the spectrum of gene fusions in cancer and how the methods to identify them have evolved, and also discuss conceptual implications of current, sequencing-based approaches for detection.

The landscape and therapeutic relevance of cancer-associated transcript fusions

Transcript fusions as a result of chromosomal rearrangements have been a focus of attention in cancer as they provide attractive therapeutic targets. To identify novel fusion transcripts with the potential to be exploited therapeutically, we analyzed RNA sequencing, DNA copy number and gene mutation data from 4,366 primary tumor samples. To avoid false positives, we implemented stringent quality criteria that included filtering of fusions detected in RNAseq data from 364 normal tissue samples. Our analysis identified 7,887 high confidence fusion transcripts across 13 tumor types. Our fusion prediction was validated by evidence of a genomic rearrangement for 78 of 79 fusions in 48 glioma samples where whole genome sequencing data was available. Cancers with higher levels of genomic instability showed a corresponding increase in fusion transcript frequency, whereas tumor samples harboring fusions contained statistically significantly fewer driver gene mutations, suggesting an important role for tumorigenesis. We identified at least one in-frame protein kinase fusion in 324 of 4,366 samples (7.4%). Potentially druggable kinase fusions involving ALK, ROS, RET, NTRK, and FGFR gene families were detected in bladder carcinoma (3.3%), glioblastoma (4.4%), head and neck cancer (1.0%), low grade glioma (1.5%), lung adenocarcinoma (1.6%), lung squamous cell carcinoma (2.3%), and thyroid carcinoma (8.7%), suggesting a potential for application of kinase inhibitors across tumor types. In-frame fusion transcripts involving histone methyltransferase or histone demethylase genes were detected in 111 samples (2.5%) and may additionally be considered as therapeutic targets. In summary, we described the landscape of transcript fusions detected across a large number of tumor samples and revealed fusion events with clinical relevance that have not been previously recognized. Our results support the concept of basket clinical trials where patients are matched with experimental therapies based on their genomic profile rather than the tissue where the tumor originated.

The landscape of kinase fusions in cancer

Human cancer genomes harbour a variety of alterations leading to the deregulation of key pathways in tumour cells. The genomic characterization of tumours has uncovered numerous genes recurrently mutated, deleted or amplified, but gene fusions have not been characterized as extensively. Here we develop heuristics for reliably detecting gene fusion events in RNA-seq data and apply them to nearly 7,000 samples from The Cancer Genome Atlas. We thereby are able to discover several novel and recurrent fusions involving kinases. These findings have immediate clinical implications and expand the therapeutic options for cancer patients, as approved or exploratory drugs exist for many of these kinases.

Kinases activated by gene fusions represent potentially important targets for the development of cancer drugs. Here, the authors develop a method for detecting gene fusion events in RNA sequencing data from The Cancer Genome Atlas and identify several novel recurrent fusions involving kinases.

Cancer subclonal genetic architecture as a key to personalized medicine

The future of personalized oncological therapy will likely rely on evidence-based medicine to integrate all of the available evidence to delineate the most efficacious treatment option for the patient. To undertake evidence-based medicine through use of targeted therapy regimens, identification of the specific underlying causative mutation(s) driving growth and progression of a patient's tumor is imperative. Although molecular subtyping is important for planning and treatment, intraclonal genetic diversity has been recently highlighted as having significant implications for biopsy-based prognosis. Overall, delineation of the clonal architecture of a patient's cancer and how this will impact on the selection of the most efficacious therapy remain a topic of intense interest.

Characterization of the genomic features and expressed fusion genes in micropapillary carcinomas of the breast

Micropapillary carcinoma (MPC) is a rare histological special type of breast cancer, characterized by an aggressive clinical behaviour and a pattern of copy number aberrations (CNAs) distinct from that of grade- and oestrogen receptor (ER)-matched invasive carcinomas of no special type (IC-NSTs). The aims of this study were to determine whether MPCs are underpinned by a recurrent fusion gene(s) or mutations in 273 genes recurrently mutated in breast cancer. Sixteen MPCs were subjected to microarray-based comparative genomic hybridization (aCGH) analysis and Sequenom OncoCarta mutation analysis. Eight and five MPCs were subjected to targeted capture and RNA sequencing, respectively. aCGH analysis confirmed our previous observations about the repertoire of CNAs of MPCs. Sequencing analysis revealed a spectrum of mutations similar to those of luminal B IC-NSTs, and recurrent mutations affecting mitogen-activated protein kinase family genes and NBPF10. RNA-sequencing analysis identified 17 high-confidence fusion genes, eight of which were validated and two of which were in-frame. No recurrent fusions were identified in an independent series of MPCs and IC-NSTs. Forced expression of in-frame fusion genes (SLC2A1-FAF1 and BCAS4-AURKA) resulted in increased viability of breast cancer cells. In addition, genomic disruption of CDK12 caused by out-of-frame rearrangements was found in one MPC and in 13% of HER2-positive breast cancers, identified through a re-analysis of publicly available massively parallel sequencing data. In vitro analyses revealed that CDK12 gene disruption results in sensitivity to PARP inhibition, and forced expression of wild-type CDK12 in a CDK12-null cell line model resulted in relative resistance to PARP inhibition. Our findings demonstrate that MPCs are neither defined by highly recurrent mutations in the 273 genes tested, nor underpinned by a recurrent fusion gene. Although seemingly private genetic events, some of the fusion transcripts found in MPCs may play a role in maintenance of a malignant phenotype and potentially offer therapeutic opportunities.