Accuracy assessment of fusion transcript detection via read-mapping and de novo fusion transcript assembly-based methods

Mesenchymal tumor organoid models recapitulate rhabdomyosarcoma subtypes

Rhabdomyosarcomas (RMS) are mesenchyme-derived tumors and the most common childhood soft tissue sarcomas. Treatment is intense, with a nevertheless poor prognosis for high-risk patients. Discovery of new therapies would benefit from additional preclinical models. Here, we describe the generation of a collection of 19 pediatric RMS tumor organoid (tumoroid) models (success rate of 41%) comprising all major subtypes. For aggressive tumors, tumoroid models can often be established within 4-8 weeks, indicating the feasibility of personalized drug screening. Molecular, genetic, and histological characterization show that the models closely resemble the original tumors, with genetic stability over extended culture periods of up to 6 months. Importantly, drug screening reflects established sensitivities and the models can be modified by CRISPR/Cas9 with TP53 knockout in an embryonal RMS model resulting in replicative stress drug sensitivity. Tumors of mesenchymal origin can therefore be used to generate organoid models, relevant for a variety of preclinical and clinical research questions.

Experimentally Deduced Criteria for Detection of Clinically Relevant Fusion 3′ Oncogenes from FFPE Bulk RNA Sequencing Data

Drugs targeting receptor tyrosine kinase (RTK) oncogenic fusion proteins demonstrate impressive anti-cancer activities. The fusion presence in the cancer is the respective drug prescription biomarker, but their identification is challenging as both the breakpoint and the exact fusion partners are unknown. RNAseq offers the advantage of finding both fusion parts by screening sequencing reads. Paraffin (FFPE) tissue blocks are the most common way of storing cancer biomaterials in biobanks. However, finding RTK fusions in FFPE samples is challenging as RNA fragments are short and their artifact ligation may appear in sequencing libraries. Here, we annotated RNAseq reads of 764 experimental FFPE solid cancer samples, 96 leukemia samples, and 2 cell lines, and identified 36 putative clinically relevant RTK fusions with junctions corresponding to exon borders of the fusion partners. Where possible, putative fusions were validated by RT-PCR (confirmed for 10/25 fusions tested). For the confirmed 3′RTK fusions, we observed the following distinguishing features. Both moieties were in-frame, and the tyrosine kinase domain was preserved. RTK exon coverage by RNAseq reads upstream of the junction site were lower than downstream. Finally, most of the true fusions were present by more than one RNAseq read. This provides the basis for automatic annotation of 3′RTK fusions using FFPE RNAseq profiles.

Accurate detection of tumor-specific gene fusions reveals strongly immunogenic personal neo-antigens

Cancer-associated gene fusions are a potential source for highly immunogenic neoantigens, but the lack of computational tools for accurate, sensitive identification of personal gene fusions has limited their targeting in personalized cancer immunotherapy. Here we present EasyFuse, a machine learning computational pipeline for detecting cancer-specific gene fusions in transcriptome data obtained from human cancer samples. EasyFuse predicts personal gene fusions with high precision and sensitivity, outperforming previously described tools. By testing immunogenicity with autologous blood lymphocytes from patients with cancer, we detected pre-established CD4⁺ and CD8⁺ T cell responses for 10 of 21 (48%) and for 1 of 30 (3%) identified gene fusions, respectively. The high frequency of T cell responses detected in patients with cancer supports the relevance of individual gene fusions as neoantigens that might be targeted in personalized immunotherapies, especially for tumors with low mutation burden.

Molecular profiling identifies targeted therapy opportunities in pediatric solid cancer

To evaluate the clinical impact of molecular tumor profiling (MTP) with targeted sequencing panel tests, pediatric patients with extracranial solid tumors were enrolled in a prospective observational cohort study at 12 institutions. In the 345-patient analytical population, median age at diagnosis was 12 years (range 0-27.5); 298 patients (86%) had 1 or more alterations with potential for impact on care. Genomic alterations with diagnostic, prognostic or therapeutic significance were present in 61, 16 and 65% of patients, respectively. After return of the results, impact on care included 17 patients with a clarified diagnostic classification and 240 patients with an MTP result that could be used to select molecularly targeted therapy matched to identified alterations (MTT). Of the 29 patients who received MTT, 24% had an objective response or experienced durable clinical benefit; all but 1 of these patients received targeted therapy matched to a gene fusion. Of the diagnostic variants identified in 209 patients, 77% were gene fusions. MTP with targeted panel tests that includes fusion detection has a substantial clinical impact for young patients with solid tumors.

Molecular Characterization of the Tumor Microenvironment in Renal Medullary Carcinoma

Renal medullary carcinoma (RMC) is a highly aggressive disease associated with sickle hemoglobinopathies and universal loss of the tumor suppressor gene SMARCB1. RMC has a relatively low rate of incidence compared with other renal cell carcinomas (RCCs) that has hitherto made molecular profiling difficult. To probe this rare disease in detail we performed an in-depth characterization of the RMC tumor microenvironment using a combination of genomic, metabolic and single-cell RNA-sequencing experiments on tissue from a representative untreated RMC patient, complemented by retrospective analyses of archival tissue and existing published data. Our study of the tumor identifies a heterogenous population of malignant cell states originating from the thick ascending limb of the Loop of Henle within the renal medulla. Transformed RMC cells displayed the hallmarks of increased resistance to cell death by ferroptosis and proteotoxic stress driven by MYC-induced proliferative signals. Specifically, genomic characterization of RMC tumors provides substantiating evidence for the recently proposed dependence of SMARCB1-difficient cancers on proteostasis modulated by an intact CDKN2A-p53 pathway. We also provide evidence that increased cystine-mTORC-GPX4 signaling plays a role in protecting transformed RMC cells against ferroptosis. We further propose that RMC has an immune landscape comparable to that of untreated RCCs, including heterogenous expression of the immune ligand CD70 within a sub-population of tumor cells. The latter could provide an immune-modulatory role that serves as a viable candidate for therapeutic targeting.

A Comprehensive Understanding of the Genomic Bone Tumor Landscape: A Multicenter Prospective Study

Complexity and heterogeneity increases the difficulty of diagnosis and treatment of bone tumors. We aimed to identify the mutational characterization and potential biomarkers of bone tumors. In this study, a total of 357 bone tumor patients were recruited and the next generation sequencing (NGS)-based YuanSu450 panel, that includes both DNA and RNA sequencing, was performed for genomic alteration identification. The most common mutated genes in bone tumors included TP53, NCOR1, VEGFA, RB1, CCND3, CDKN2A, GID4, CCNE1, TERT, and MAP2K4. The amplification of genes such as NCOR1, VEGFA, and CCND3 mainly occurred in osteosarcoma. Germline mutation analysis reveal a high frequency of HRD related mutations (46.4%, 13/28) in this cohort. With the assistance of RNA sequencing, 16.8% (19/113) gene fusions were independently detected in 20% (16/79) of patients. Nearly 34.2% of patients harbored actionable targeted mutations, of which the most common mutation is CDKN2A deletion. The different mutational characterizations between juvenile patients and adult patients indicated the potential effect of age in bone tumor treatment. According to the genomic alterations, the diagnosis of 26 (7.28%) bone tumors were corrected. The most easily misdiagnosed bone tumor included malignant giant cell tumors of bone (2.8%, 10/357) and fibrous dysplasia of bone (1.7%, 6/357). Meanwhile, we found that the mutations of MUC16 may be a potential biomarker for the diagnosis of mesenchymal chondrosarcomas. Our results indicated that RNA sequencing effectively complements DNA sequencing and increased the detection rate of gene fusions, supporting that NGS technology can effectively assist the diagnosis of bone tumors.

Circular and Fusion RNAs in Medulloblastoma Development

Simple Summary

Expression of circular RNAs is known to be deregulated in cancer. Here the most comprehensive set of differentially expressed RNA circles in medulloblastoma compared to cerebellum is provided. Additionally, fusion RNAs are also identified in both cancerous and normal cerebellar tissue. Some of the fusions detected in medulloblastoma are generated by genomic rearrangements that link different genes. However, fusion RNAs are also detected in normal cerebellum. In fact, there are cases where the same fusion RNA is also found in medulloblastoma. This observation underscores that the formation of fusion transcripts may not be limited to chromosomal events but could also result from mechanisms that act at the RNA level. These include read-through transcription of neighboring genes and intermolecular splicing of pre-mRNAs from different genes Importantly, these RNA “recombination” events occur not only in normal but also in cancerous tissue.

Abstract

Background. The cerebellar cancer medulloblastoma is the most common childhood cancer in the brain. Methods. RNA sequencing of 81 human biospecimens of medulloblastoma using pipelines to detect circular and fusion RNAs. Validation via PCR and Sanger sequencing. Results. 27, 56, 28 and 11 RNA circles were found to be uniquely up-regulated, while 149, 7, 20 and 15 uniquely down-regulated in the SHH, WNT, Group 3, and Group 4 medulloblastoma subtypes, respectively. Moreover, linear and circular fusion RNAs containing exons from distinct genes joined at canonical splice sites were also identified. These were generally expressed less than the circular RNAs, however the expression of both the linear and the circular fusions was comparable. Importantly, the expression of the fusions in medulloblastoma was also comparable to that of cerebellum. Conclusions. A significant number of fusions in tumor may be generated by mechanisms similar to the ones generating fusions in normal tissue. Some fusions could be rationalized by read-through transcription of two neighboring genes. However, for other fusions, e.g., a linear fusion with an exon from a downstream gene joined 5′ to 3′ with an exon from an upstream gene, more complicated splicing mechanisms, e.g., trans-splicing, have to be postulated.

Recurrent chromosomal translocations in sarcomas create a megacomplex that mislocalizes NuA4/TIP60 to Polycomb target loci

Chromosomal translocations frequently promote carcinogenesis by producing gain-of-function fusion proteins. Recent studies have identified highly recurrent chromosomal translocations in patients with endometrial stromal sarcomas (ESSs) and ossifying fibromyxoid tumors (OFMTs), leading to an in-frame fusion of PHF1 (PCL1) to six different subunits of the NuA4/TIP60 complex. While NuA4/TIP60 is a coactivator that acetylates chromatin and loads the H2A.Z histone variant, PHF1 is part of the Polycomb repressive complex 2 (PRC2) linked to transcriptional repression of key developmental genes through methylation of histone H3 on lysine 27. In this study, we characterize the fusion protein produced by the EPC1-PHF1 translocation. The chimeric protein assembles a megacomplex harboring both NuA4/TIP60 and PRC2 activities and leads to mislocalization of chromatin marks in the genome, in particular over an entire topologically associating domain including part of the HOXD cluster. This is linked to aberrant gene expression-most notably increased expression of PRC2 target genes. Furthermore, we show that JAZF1-implicated with a PRC2 component in the most frequent translocation in ESSs, JAZF1-SUZ12-is a potent transcription activator that physically associates with NuA4/TIP60, its fusion creating outcomes similar to those of EPC1-PHF1 Importantly, the specific increased expression of PRC2 targets/HOX genes was also confirmed with ESS patient samples. Altogether, these results indicate that most chromosomal translocations linked to these sarcomas use the same molecular oncogenic mechanism through a physical merge of NuA4/TIP60 and PRC2 complexes, leading to mislocalization of histone marks and aberrant Polycomb target gene expression.

The DNA methylation landscape of five pediatric-tumor types

Fewer DNA mutations have been identified in pediatric tumors than in adult tumors, suggesting that alternative tumorigenic mechanisms, including aberrant DNA methylation, may play a prominent role. In one epigenetic process of regulating gene expression, methyl groups are attached at the 5-carbon of the cytosine ring, leading to 5-methylcytosine (5mC). In somatic cells, 5mC occurs mostly in CpG islands, which are often within promoter regions. In Wilms tumors and acute myeloid leukemias, increased levels of epigenetic silencing have been associated with worse patient outcomes. However, to date, researchers have studied methylation primarily in adult tumors and for specific genes-but not on a pan-pediatric cancer scale. We addressed these gaps first by aggregating methylation data from 309 noncancerous samples, establishing baseline expectations for each probe and gene. Even though these samples represent diverse, noncancerous tissue types and population ancestral groups, methylation levels were consistent for most genes. Second, we compared tumor methylation levels against the baseline values for 489 pediatric tumors representing five cancer types: Wilms tumors, clear cell sarcomas of the kidney, rhabdoid tumors, neuroblastomas, and osteosarcomas. Tumor hypomethylation was more common than hypermethylation, and as many as 41.7% of genes were hypomethylated in a given tumor, compared to a maximum of 34.2% for hypermethylated genes. However, in known oncogenes, hypermethylation was more than twice as common as in other genes. We identified 139 probes (31 genes) that were differentially methylated between at least one tumor type and baseline levels, and 32 genes that were differentially methylated across the pediatric tumor types. We evaluated whether genomic events and aberrant methylation were mutually exclusive but did not find evidence of this phenomenon.

Identification of a novel PHIP::BRAF gene fusion in infantile fibrosarcoma

INTRODUCTION

The ETV6::NTRK3 fusion is the most common gene alteration in infantile fibrosarcoma, a soft tissue tumor affecting patients under two years of age. Less frequently, these tumors harbor fusions of genes encoding other kinases, such as BRAF, which activates MEK in the mitogen-activated protein kinase pathway. The identification and characterization of these oncogenes is crucial to facilitate diagnosis, validate new treatments and better understand the pathophysiology of these neoplasms.

METHODS

Herein, we analyzed an ETV6::NTRK3-negative infantile fibrosarcoma from a 5-day-old patient by RNA-sequencing to identify new fusion transcripts. Functional exploration of the fusion of interest was performed by in vitro assays to study its activity, oncogenicity and sensitivity to the MEK inhibitor trametinib.

RESULTS

We identified a novel fusion involving the PHIP and BRAF genes. The corresponding fusion protein constitutively activated the mitogen-activated protein kinase pathway, resulting in fibroblast transformation. Treatment of transfected cells with trametinib effectively inhibited signaling by PHIP::BRAF.

CONCLUSION

PHIP::BRAF is a novel fusion oncogene that can be targeted by trametinib in infantile fibrosarcoma. This article is protected by copyright. All rights reserved.

Comprehensive genomic and tumour immune profiling reveals potential therapeutic targets in malignant pleural mesothelioma

Background

Malignant pleural mesothelioma (MPM) has a poor overall survival with few treatment options. Whole genome sequencing (WGS) combined with the immune features of MPM offers the prospect of identifying changes that could inform future clinical trials.

Methods

We analysed somatic mutations from 229 MPM samples, including previously published data and 58 samples that had undergone WGS within this study. This was combined with RNA-seq analysis to characterize the tumour immune environment.

Results

The comprehensive genome analysis identified 12 driver genes, including new candidate genes. Whole genome doubling was a frequent event that correlated with shorter survival. Mutational signature analysis revealed SBS5/40 were dominant in 93% of samples, and defects in homologous recombination repair were infrequent in our cohort. The tumour immune environment contained high M2 macrophage infiltrate linked with MMP2, MMP14, TGFB1 and CCL2 expression, representing an immune suppressive environment. The expression of TGFB1 was associated with overall survival. A small subset of samples (less than 10%) had a higher proportion of CD8 T cells and a high cytolytic score, suggesting a ‘hot’ immune environment independent of the somatic mutations.

Conclusions

We propose accounting for genomic and immune microenvironment status may influence therapeutic planning in the future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01060-8.

Detecting anaplastic lymphoma kinase (ALK) gene rearrangements with next-generation sequencing remains a reliable approach in patients with non-small-cell lung cancer

Next-generation sequencing (NGS) is rapidly becoming routine in clinical oncology practice to identify therapeutic biomarkers, including gene rearrangements in anaplastic lymphoma kinase (ALK). Our study investigated the concordance of ALK positivity evaluated by DNA-based NGS with orthogonal ALK testing methods such as fluorescence in situ hybridization (FISH), immunohistochemistry (IHC), and RNA-based NGS (RNA-NGS). Thirty-eight patients with lung adenocarcinoma who were detected with ALK rearrangements using DNA-NGS and also had adequate tissue samples submitted for FISH, IHC, and RNA-NGS, were included in this study. Of the 38 patients, RNA samples from 3 patients failed quality control for RNA-NGS. The concordance of ALK positivity was calculated relative to DNA-NGS results. The concordance rates were 97.1% (34/35) for RNA-NGS, 94.7% (36/38) for IHC, and 97.4% (37/38) for FISH. DNA-NGS detected single ALK rearrangements in 14 (35.0%) patients and complex ALK rearrangements in 26 (65.0%). RNA-NGS detected only single transcripts of the primary ALK fusions. A novel LANCL1-ALK (L7:A20) detected using DNA-NGS was detected as EML4-ALK (E13:A20) transcripts using RNA-NGS. Interestingly, patients with single ALK rearrangements were more likely to be detected with atypical isolated red signals (p < 0.001), while patients with complex ALK rearrangements were more likely to be detected with atypical split red and green signals less than 2 signal diameters apart (p < 0.001). Our study highlights the reliability of NGS in the accurate detection of specific ALK fusion variants and concomitant mutations that are crucial for individualized treatment decisions in patients with lung cancer.

Feasibility of whole genome and transcriptome profiling in pediatric and young adult cancers

The utility of cancer whole genome and transcriptome sequencing (cWGTS) in oncology is increasingly recognized. However, implementation of cWGTS is challenged by the need to deliver results within clinically relevant timeframes, concerns about assay sensitivity, reporting and prioritization of findings. In a prospective research study we develop a workflow that reports comprehensive cWGTS results in 9 days. Comparison of cWGTS to diagnostic panel assays demonstrates the potential of cWGTS to capture all clinically reported mutations with comparable sensitivity in a single workflow. Benchmarking identifies a minimum of 80× as optimal depth for clinical WGS sequencing. Integration of germline, somatic DNA and RNA-seq data enable data-driven variant prioritization and reporting, with oncogenic findings reported in 54% more patients than standard of care. These results establish key technical considerations for the implementation of cWGTS as an integrated test in clinical oncology.

Cancer whole-genome and transcriptome sequencing (cWGTS) has been challenging to implement in clinical settings. Here, the authors develop a workflow to deliver robust cWGTS analyses and reports within clinically-relevant timeframes for paediatric, adolescent and young adult solid tumour patients.

Phase II clinical trial of Eribulin-gemcitabine combination therapy in previously treated patients with advanced liposarcoma or leiomyosarcoma

PURPOSE

Monotherapy with eribulin or gemcitabine has been found to be moderately effective in treating soft-tissue sarcomas (STS). In this study, we evaluated the efficacy and safety of eribulin-gemcitabine combination therapy for the two most common histologic types of STS, liposarcoma and leiomyosarcoma.

PATIENTS AND METHODS

In this non-randomized, multicenter, phase II study, we included patients with progressive disease who had received one or two courses of chemotherapy that included doxorubicin. Patients were administered 1.4 mg/m2 eribulin and 1,000 mg/m2 gemcitabine on days 1 and 8 every 3 weeks. The primary endpoint was progression-free survival rate at 12 weeks (PFSR12wks), with null and alternative hypotheses of PFSR12wks {less than or equal to}20.0% and {greater than or equal to}40.0%, respectively. Exploratory biomarker analyses with next-generation sequencing (NGS) were performed on pretreatment tumor samples.

RESULTS

Among the 37 patients included, the overall PFSR12wks was 73.0%, achieving the primary endpoint. The objective response rate, disease control rate, median progression-free survival, and median overall survival were 16.2%, 78.4%, 5.6 months, and 31.9 months, respectively, without differences according to histologic type. New safety signals and treatment-related deaths were not documented. NGS-based transcriptome analysis revealed that functional enrichment in the transforming growth factor-ß pathway was mostly associated with a poor outcome, whereas single genetic alterations largely failed to predict treatment outcome.

CONCLUSIONS

Eribulin-gemcitabine combination therapy showed promising activity and an acceptable safety profile in patients with liposarcoma or leiomyosarcoma. Gene expression profiling with pathway enrichment analysis would have possibilities to have predictive value for survival outcome, necessitating further investigation to confirm.

Classification and clustering of RNA crosslink-ligation data reveal complex structures and homodimers

The recent development and application of methods based on the general principle of "crosslinking and proximity ligation" (crosslink-ligation) are revolutionizing RNA structure studies in living cells. However, extracting structure information from such data presents unique challenges. Here, we introduce a set of computational tools for the systematic analysis of data from a wide variety of crosslink-ligation methods, specifically focusing on read mapping, alignment classification, and clustering. We design a new strategy to map short reads with irregular gaps at high sensitivity and specificity. Analysis of previously published data reveals distinct properties and bias caused by the crosslinking reactions. We perform rigorous and exhaustive classification of alignments and discover eight types of arrangements that provide distinct information on RNA structures and interactions. To deconvolve the dense and intertwined gapped alignments, we develop a network/graph-based tool Crosslinked RNA Secondary Structure Analysis using Network Techniques (CRSSANT), which enables clustering of gapped alignments and discovery of new alternative and dynamic conformations. We discover that multiple crosslinking and ligation events can occur on the same RNA, generating multisegment alignments to report complex high-level RNA structures and multi-RNA interactions. We find that alignments with overlapped segments are produced from potential homodimers and develop a new method for their de novo identification. Analysis of overlapping alignments revealed potential new homodimers in cellular noncoding RNAs and RNA virus genomes in the Picornaviridae family. Together, this suite of computational tools enables rapid and efficient analysis of RNA structure and interaction data in living cells.

Comprehensive assessment of anaplastic lymphoma kinase in localized and metastatic prostate cancer reveals targetable alterations

Anaplastic lymphoma kinase (ALK) is a tyrosine kinase with genomic and expression changes in many solid tumors. ALK inhibition is first line therapy for lung cancers with ALK alterations, and an effective therapy in other tumor types, but has not been well-studied in prostate cancer. Here, we aim to delineate the role of ALK genomic and expression changes in primary and metastatic prostate cancer. We determined ALK expression by immunohistochemistry and RNA-Seq, and genomic alterations by NGS. We assessed functional consequences of ALK overexpression and pharmacological ALK inhibition by cell proliferation and cell viability assays. Among 372 primary prostate cancer cases we identified one case with uniformly high ALK protein expression. Genomic analysis revealed a SLC45A3-ALK fusion which promoted oncogenesis in in vitro assays. We observed ALK protein expression in 5/52 (9%) of metastatic prostate cancer cases, of which 4 of 5 had neuroendocrine features. ALK-expressing neuroendocrine prostate cancer had a distinct transcriptional program, and earlier disease progression. An ALK-expressing neuroendocrine prostate cancer model was sensitive to pharmacological ALK inhibition. In summary, we found that ALK overexpression is rare in primary prostate cancer, but more frequent in metastatic prostate cancers with neuroendocrine differentiation. Further, ALK fusions similar to lung cancer are an occasional driver in prostate cancer. Our data suggest that ALK-directed therapies could be an option in selected patients with advanced prostate cancer.

Computational analysis of cancer genome sequencing data

Distilling biologically meaningful information from cancer genome sequencing data requires comprehensive identification of somatic alterations using rigorous computational methods. As the amount and complexity of sequencing data have increased, so has the number of tools for analysing them. Here, we describe the main steps involved in the bioinformatic analysis of cancer genomes, review key algorithmic developments and highlight popular tools and emerging technologies. These tools include those that identify point mutations, copy number alterations, structural variations and mutational signatures in cancer genomes. We also discuss issues in experimental design, the strengths and limitations of sequencing modalities and methodological challenges for the future.

Transcriptomic Portraits and Molecular Pathway Activation Features of Adult Spinal Intramedullary Astrocytomas

In this study, we report 31 spinal intramedullary astrocytoma (SIA) RNA sequencing (RNA-seq) profiles for 25 adult patients with documented clinical annotations. To our knowledge, this is the first clinically annotated RNA-seq dataset of spinal astrocytomas derived from the intradural intramedullary compartment. We compared these tumor profiles with the previous healthy central nervous system (CNS) RNA-seq data for spinal cord and brain and identified SIA-specific gene sets and molecular pathways. Our findings suggest a trend for SIA-upregulated pathways governing interactions with the immune cells and downregulated pathways for the neuronal functioning in the context of normal CNS activity. In two patient tumor biosamples, we identified diagnostic KIAA1549-BRAF fusion oncogenes, and we also found 16 new SIA-associated fusion transcripts. In addition, we bioinformatically simulated activities of targeted cancer drugs in SIA samples and predicted that several tyrosine kinase inhibitory drugs and thalidomide analogs could be potentially effective as second-line treatment agents to aid in the prevention of SIA recurrence and progression.

Identifying the oncogenic potential of gene fusions exploiting miRNAs

It is estimated that oncogenic gene fusions cause about 20% of human cancer morbidity. Identifying potentially oncogenic gene fusions may improve affected patients' diagnosis and treatment. Previous approaches to this issue included exploiting specific gene-related information, such as gene function and regulation. Here we propose a model that profits from the previous findings and includes the microRNAs in the oncogenic assessment. We present ChimerDriver, a tool to classify gene fusions as oncogenic or not oncogenic. ChimerDriver is based on a specifically designed neural network and trained on genetic and post-transcriptional information to obtain a reliable classification. The designed neural network integrates information related to transcription factors, gene ontologies, microRNAs and other detailed information related to the functions of the genes involved in the fusion and the gene fusion structure. As a result, the performances on the test set reached 0.83 f1-score and 96% recall. The comparison with state-of-the-art tools returned comparable or higher results. Moreover, ChimerDriver performed well in a real-world case where 21 out of 24 validated gene fusion samples were detected by the gene fusion detection tool Starfusion. ChimerDriver integrates transcriptional and post-transcriptional information in an ad-hoc designed neural network to effectively discriminate oncogenic gene fusions from passenger ones. ChimerDriver source code is freely available at https://github.com/martalovino/ChimerDriver.

Unscrambling cancer genomes via integrated analysis of structural variation and copy number

Complex somatic genomic rearrangements and copy number alterations are hallmarks of nearly all cancers. We have developed an algorithm, LINX, to aid interpretation of structural variant and copy number data derived from short-read, whole-genome sequencing. LINX classifies raw structural variant calls into distinct events and predicts their effect on the local structure of the derivative chromosome and the functional impact on affected genes. Visualizations facilitate further investigation of complex rearrangements. LINX allows insights into a diverse range of structural variation events and can reliably detect pathogenic rearrangements, including gene fusions, immunoglobulin enhancer rearrangements, intragenic deletions, and duplications. Uniquely, LINX also predicts chained fusions that we demonstrate account for 13% of clinically relevant oncogenic fusions. LINX also reports a class of inactivation events that we term homozygous disruptions that may be a driver mutation in up to 9% of tumors and may frequently affect PTEN, TP53, and RB1.

Cancer Neoantigens: Challenges and Future Directions for Prediction, Prioritization, and Validation

Prioritization of immunogenic neoantigens is key to enhancing cancer immunotherapy through the development of personalized vaccines, adoptive T cell therapy, and the prediction of response to immune checkpoint inhibition. Neoantigens are tumor-specific proteins that allow the immune system to recognize and destroy a tumor. Cancer immunotherapies, such as personalized cancer vaccines, adoptive T cell therapy, and immune checkpoint inhibition, rely on an understanding of the patient-specific neoantigen profile in order to guide personalized therapeutic strategies. Genomic approaches to predicting and prioritizing immunogenic neoantigens are rapidly expanding, raising new opportunities to advance these tools and enhance their clinical relevance. Predicting neoantigens requires acquisition of high-quality samples and sequencing data, followed by variant calling and variant annotation. Subsequently, prioritizing which of these neoantigens may elicit a tumor-specific immune response requires application and integration of tools to predict the expression, processing, binding, and recognition potentials of the neoantigen. Finally, improvement of the computational tools is held in constant tension with the availability of datasets with validated immunogenic neoantigens. The goal of this review article is to summarize the current knowledge and limitations in neoantigen prediction, prioritization, and validation and propose future directions that will improve personalized cancer treatment.

Identification of alternative splicing-derived cancer neoantigens for mRNA vaccine development

Messenger RNA (mRNA) vaccines have shown great potential for anti-tumor therapy due to the advantages in safety, efficacy and industrial production. However, it remains a challenge to identify suitable cancer neoantigens that can be targeted for mRNA vaccines. Abnormal alternative splicing occurs in a variety of tumors, which may result in the translation of abnormal transcripts into tumor-specific proteins. High-throughput technologies make it possible for systematic characterization of alternative splicing as a source of suitable target neoantigens for mRNA vaccine development. Here, we summarized difficulties and challenges for identifying alternative splicing-derived cancer neoantigens from RNA-seq data and proposed a conceptual framework for designing personalized mRNA vaccines based on alternative splicing-derived cancer neoantigens. In addition, several points were presented to spark further discussion toward improving the identification of alternative splicing-derived cancer neoantigens.

Transcriptomic analysis of papillary thyroid cancer focused on immune-subtyping, oncogenic fusion, and recurrence

Objectives

Thyroid cancer is the most common endocrine tumor, with rapidly increasing incidence worldwide. However, its transcriptomic characteristics associated with immunological signatures, driver fusions, and recurrence markers remain unclear. We aimed to investigate the transcriptomic characteristics of advanced papillary thyroid cancer.

Methods

This study included 282 papillary thyroid cancer tumor samples and 155 normal samples from Chungnam National University Hospital and Seoul National University Hospital. Transcriptomic quantification was determined by high-throughput RNA sequencing. We investigated the associations of clinical parameters and molecular signatures using RNA sequencing. We validated predictive biomarkers using the Cancer Genome Atlas database.

Results

Through a comparison of differentially expressed genes, gene sets, and pathways in papillary thyroid cancer compared to normal tumor-adjacent tissue, we found increased immune signaling associated with cytokines or T cells and decreased thyroid hormone synthetic pathways. In addition, patients with recurrence presented increased CD8+ T-cell and Th1-cell signatures. Interestingly, we found differentially overexpressed genes related to immune-escape signaling such as CTLA4, IDO1, LAG3, and PDCD1 in advanced papillary thyroid cancer with a low thyroid differentiation score. Fusion analysis showed that the PI3K and mitogen-activated protein kinase (MAPK) signaling pathways were regulated differently according to the RET fusion partner genes (CCDC6 or NCOA4). Finally, we identified HOXD9 as a novel molecular biomarker that predicts the recurrence of thyroid cancer in addition to known risk factors (tumor size, lymph node metastasis, and extrathyroidal extension).

Conclusion

We identified a high association with immune-escape signaling in the immune-hot group with aggressive clinical characteristics among Korean thyroid cancer patients. Moreover, RET fusion differentially regulated PI3K and MAPK signaling depending on the partner gene of RET, and HOXD9 was found to be a recurrence marker for advanced papillary thyroid cancer.

FusionAI, a DNA-sequence-based deep learning protocol reduces the false positives of human fusion gene prediction

Even though there were many tool developments of fusion gene prediction from NGS data, too many false positives are still an issue. Wise use of the genomic features around the fusion gene breakpoints will be helpful to identify reliable fusion genes efficiently. For this aim, we developed FusionAI, a deep learning pipeline predicting human fusion gene breakpoints from DNA sequence. FusionAI is freely available via https://compbio.uth.edu/FusionGDB2/FusionAI.

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

•

FusionAI can predict the fusion breakpoints from the given DNA sequence

•

FusionAI can reduce the false positives of the predicted fusion genes by other tools

•

FusionAI can identify the genomic features related to the genomic breakage

•

FusionAI creates a landscape image of 44 human genomic features around the breakpoints

Papillary Renal Neoplasm With Reverse Polarity Is Often Cystic: Report of 7 Cases and Review of 93 Cases in the Literature

Papillary renal neoplasm with reverse polarity (PRNRP) is a newly proposed entity with distinct histology and frequent KRAS mutations. To date, 93 cases of PRNRPs have been reported. In this study, we present 7 new cases of PRNRP and review the literature. Most of the pathologic features in our 7 cases are similar to those previously reported cases. However, all 7 of our cases showed at least partial cystic changes, which was not stressed in prior studies. Single-nucleotide polymorphism-microarray based chromosomal analysis demonstrated no trisomy or other alteration of chromosomes 7 or 17; and no loss or other alteration of chromosome Y was detected in all 7 cases. Next-generation sequencing detected KRAS missense mutations in 4 of 7 cases. No fusion genes were detected. In summary, PRNRP is a small, well-circumscribed often encapsulated and cystic neoplasm with loose papillary formations. Cuboidal tumor cells always have eosinophilic cytoplasm and nuclei located at the pole opposite the basement membrane with a low World Health Organization (WHO)/International Society of Urologic Pathologists (ISUP) nuclear grade. The fibrovascular cores can be hyalinized or edematous. Macrophage aggregates and intracellular hemosiderin are uncommon, and no psammoma bodies or necrosis should be seen. Immunophenotypically, this tumor is always positive for CK7 and GATA3, and negative for CD117 and vimentin. CD10 and AMACR can be positive, but often weakly and focally. PRNRP often has KRAS mutations, however, only 32% of cases have chromosomal abnormalities in chromosomes 7, 17, and Y. No recurrences, metastases, or tumor-related deaths have been reported following complete resection.

Shedding light on the hidden human proteome expands immunopeptidome in cancer

Unrestrained cellular growth and immune escape of a tumor are associated with the incidental errors of the genome and transcriptome. Advances in next-generation sequencing have identified thousands of genomic and transcriptomic aberrations that generate variant peptides that assemble the hidden proteome, further expanding the immunopeptidome. Emerging next-generation sequencing technologies and a number of computational methods estimated the abundance of immune infiltration from bulk transcriptome have advanced our understanding of tumor microenvironments. Here, we will characterize several major types of tumor-specific antigens arising from single-nucleotide variants, insertions and deletions, gene fusion, alternative splicing, RNA editing and non-coding RNAs. Finally, we summarize the current state-of-the-art computational and experimental approaches or resources and provide an integrative pipeline for the identification of candidate tumor antigens. Together, the systematic investigation of the hidden proteome in cancer will help facilitate the development of effective and durable immunotherapy targets for cancer.

Fusion Gene Detection Using Whole-Exome Sequencing Data in Cancer Patients

Several fusion genes are directly involved in the initiation and progression of cancers. Numerous bioinformatics tools have been developed to detect fusion events, but they are mainly based on RNA-seq data. The whole-exome sequencing (WES) represents a powerful technology that is widely used for disease-related DNA variant detection. In this study, we build a novel analysis pipeline called Fuseq-WES to detect fusion genes at DNA level based on the WES data. The same method applies also for targeted panel sequencing data. We assess the method to real datasets of acute myeloid leukemia (AML) and prostate cancer patients. The result shows that two of the main AML fusion genes discovered in RNA-seq data, PML-RARA and CBFB-MYH11, are detected in the WES data in 36 and 63% of the available samples, respectively. For the targeted deep-sequencing of prostate cancer patients, detection of the TMPRSS2-ERG fusion, which is the most frequent chimeric alteration in prostate cancer, is 91% concordant with a manually curated procedure based on four other methods. In summary, the overall results indicate that it is challenging to detect fusion genes in WES data with a standard coverage of ∼ 15–30x, where fusion candidates discovered in the RNA-seq data are often not detected in the WES data and vice versa. A subsampling study of the prostate data suggests that a coverage of at least 75x is necessary to achieve high accuracy.

Genion, an accurate tool to detect gene fusion from long transcriptomics reads

BACKGROUND

The advent of next-generation sequencing technologies empowered a wide variety of transcriptomics studies. A widely studied topic is gene fusion which is observed in many cancer types and suspected of having oncogenic properties. Gene fusions are the result of structural genomic events that bring two genes closely located and result in a fused transcript. This is different from fusion transcripts created during or after the transcription process. These chimeric transcripts are also known as read-through and trans-splicing transcripts. Gene fusion discovery with short reads is a well-studied problem, and many methods have been developed. But the sensitivity of these methods is limited by the technology, especially the short read length. Advances in long-read sequencing technologies allow the generation of long transcriptomics reads at a low cost. Transcriptomic long-read sequencing presents unique opportunities to overcome the shortcomings of short-read technologies for gene fusion detection while introducing new challenges.

RESULTS

We present Genion, a sensitive and fast gene fusion detection method that can also detect read-through events. We compare Genion against a recently introduced long-read gene fusion discovery method, LongGF, both on simulated and real datasets. On simulated data, Genion accurately identifies the gene fusions and its clustering accuracy for detecting fusion reads is better than LongGF. Furthermore, our results on the breast cancer cell line MCF-7 show that Genion correctly identifies all the experimentally validated gene fusions.

CONCLUSIONS

Genion is an accurate gene fusion caller. Genion is implemented in C++ and is available at https://github.com/vpc-ccg/genion .

Splice variants denote differences between a cancer stem cell side population of EWSR1‑ERG‑based Ewing sarcoma cells, its main population and EWSR1‑FLI‑based cells

Ewing sarcoma is a challenging cancer entity, which, besides the characteristic presence of a fusion gene, is driven by multiple alternative splicing events. So far, splice variants in Ewing sarcoma cells were mainly analyzed for EWSR1‑FLI1. The present study provided a comprehensive alternative splicing study on CADO‑ES1, an Ewing model cell line for an EWSR1‑ERG fusion gene. Based on a well‑-characterized RNA‑sequencing dataset with extensive control mechanisms across all levels of analysis, the differential spliced genes in Ewing cancer stem cells were ATP13A3 and EPB41, while the main population was defined by ACADVL, NOP58 and TSPAN3. All alternatively spliced genes were further characterized by their Gene Ontology (GO) terms and by their membership in known protein complexes. These results confirm and extend previous studies towards a systematic whole‑transcriptome analysis. A highlight is the striking segregation of GO terms associated with five basic splice events. This mechanistic insight, together with a coherent integration of all observations with prior knowledge, indicates that EWSR1‑ERG is truly a close twin to EWSR1‑FLI1, but still exhibits certain individuality. Thus, the present study provided a measure of variability in Ewing sarcoma, whose understanding is essential both for clinical procedures and basic mechanistic insight.

Wholistic approach - transcriptomic analysis and beyond using archival material for molecular diagnosis

Many neoplasms remain unclassified after histopathological examination, which requires further molecular analysis. To this regard, mesenchymal neoplasms are particularly challenging due to the combination of their rarity and the large number of subtypes, and many entities still lack robust diagnostic hallmarks. RNA transcriptomic profiles have proven to be a reliable basis for the classification of previously unclassified tumors and notably for mesenchymal neoplasms. Using exome-based RNA capture sequencing on more than 5000 samples of archival material (FFPE), the combination of expression profiles analyzes (including several clustering methods), fusion genes, and small nucleotide variations has been developed at the Centre Léon Bérard (CLB) in Lyon for the molecular diagnosis of challenging neoplasms and the discovery of new entities. The molecular basis of the technique, the protocol, and the bioinformatics algorithms used are described herein, as well as its advantages and limitations.

Assessment of ALK Fusions in Uncommon Inflammatory Myofibroblastic Tumors With ALK IHC Positivity but FISH-Equivocal Findings by Targeted RNA Sequencing

CONTEXT.—

Fluorescence in situ hybridization (FISH) and immunohistochemistry (IHC) are common methods to detect ALK status in inflammatory myofibroblastic tumors (IMTs). However, equivocal ALK FISH signals and inconsistency between FISH and IHC are occasionally observed.

OBJECTIVE.—

To study the inconsistency between FISH and IHC, and clarify ALK status in IMT by targeted RNA sequencing (RNAseq).

DESIGN.—

In this study, 12 consultation cases preliminarily diagnosed as uncommon IMTs with ALK IHC positivity but FISH negativity, plus 3 ALK-positive and 3 ALK-negative IMTs, were re-analyzed by IHC, FISH, and RNAseq.

RESULTS.—

As a result, 1 case with FUS-TFCP2 fusion was detected by RNAseq, which was previously misdiagnosed as IMT. In the other 11 uncommon IMTs, 90.9% (10 of 11) showed equivocal ALK FISH signals, and all were confirmed to harbor ALK fusion by RNAseq, except for 1 failure, suggesting that a low threshold for ALK FISH might be proposed in IMT. Furthermore, RNAseq also identified IGFBP5-ALK in 1 case with ALK IHC positivity but typical FISH negativity, suggesting the possibility of false negatives for ALK FISH. For the typical IMTs, ALK fusion was identified by RNAseq in all 3 ALK-positive IMTs as expected, and additionally FN1-ROS1 fusions were identified in 2 of 3 ALK-negative IMTs.

CONCLUSIONS.—

These findings indicated that RNAseq can simultaneously detect multiple gene fusions and provide fusion forms and breakpoints, which is of great value for differential diagnosis, especially for those uncommon IMTs with equivocal FISH findings, or inconsistency between IHC and FISH.

Mutations in JAK/STAT and NOTCH1 Genes Are Enriched in Post-Transplant Lymphoproliferative Disorders

Post-transplant lymphoproliferative disorders (PTLD) are diseases occurring in immunocompromised patients after hematopoietic stem cell transplantation (HCT) or solid organ transplantation (SOT). Although PTLD occurs rarely, it may be associated with poor outcomes. In most cases, PTLD is driven by Epstein-Barr virus (EBV) infection. Few studies have investigated the mutational landscape and gene expression profile of PTLD. In our study, we performed targeted deep sequencing and RNA-sequencing (RNA-Seq) on 16 cases of florid follicular hyperplasia (FFH) type PTLD and 15 cases of other PTLD types that include: ten monomorphic (M-PTLD), three polymorphic (P-PTLD), and two classic Hodgkin lymphoma type PTLDs (CHL-PTLD). Our study identified recurrent mutations in JAK3 in five of 15 PTLD cases and one of 16 FFH-PTLD cases, as well as 16 other genes that were mutated in M-PTLD, P-PTLD, CHL-PTLD and FFH-PTLD. Digital image analysis demonstrated significant differences in single cell area, major axis, and diameter when comparing cases of M-PTLD and P-PTLD to FFH-PTLD. No morphometric relationship was identified with regards to a specific genetic mutation. Our findings suggest that immune regulatory pathways play an essential role in PTLD, with the JAK/STAT pathway affected in many PTLDs.

RNA-Seq for the detection of gene fusions in solid tumors: development and validation of the JAX FusionSeq™ 2.0 assay

Whole transcriptome sequencing (RNA-Seq) has gained prominence for the detection of fusions in solid tumors. Here, we describe the development and validation of an in-house RNA-Seq-based test system (FusionSeq™ 2.0) for the detection of clinically actionable gene fusions, in formalin-fixed paraffin-embedded (FFPE) specimens, using seventy tumor samples with varying fusion status. Conditions were optimized for RNA input of 50 ng, shown to be adequate to call known fusions at as low as 20% neoplastic content. Evaluation of assay performance between FFPE and fresh-frozen (FF) tissues exhibited little to no difference in fusion calling capability. Performance analysis of the assay validation data determined 100% accuracy, sensitivity, specificity, and reproducibility. This clinically developed and validated RNA-Seq-based approach for fusion detection in FPPE samples was shown to be on par if not superior to off-the-shelf commercially offered assays. With gene fusions implicated in a variety of cancer types, offering high-quality, low-cost molecular testing services for FFPE specimens will serve to best benefit the patient and the advancement of precision medicine in molecular oncology. KEY MESSAGES: A custom RNA-Seq-based test system (FusionSeq™ 2.0) for the detection of clinically actionable gene fusions, Evaluation of assay performance between FFPE and fresh-frozen (FF) tissues exhibited little to no difference in fusion calling capability. The assay can be performed with low RNA input and neoplastic content. Performance characteristics of the assay validation data determined 100% accuracy, sensitivity, specificity, and reproducibility.

Proteogenomic characterization identifies clinically relevant subgroups of intrahepatic cholangiocarcinoma

We performed proteogenomic characterization of intrahepatic cholangiocarcinoma (iCCA) using paired tumor and adjacent liver tissues from 262 patients. Integrated proteogenomic analyses prioritized genetic aberrations and revealed hallmarks of iCCA pathogenesis. Aflatoxin signature was associated with tumor initiation, proliferation, and immune suppression. Mutation-associated signaling profiles revealed that TP53 and KRAS co-mutations may contribute to iCCA metastasis via the integrin-FAK-SRC pathway. FGFR2 fusions activated the Rho GTPase pathway and could be a potential source of neoantigens. Proteomic profiling identified four patient subgroups (S1-S4) with subgroup-specific biomarkers. These proteomic subgroups had distinct features in prognosis, genetic alterations, microenvironment dysregulation, tumor microbiota composition, and potential therapeutics. SLC16A3 and HKDC1 were further identified as potential prognostic biomarkers associated with metabolic reprogramming of iCCA cells. This study provides a valuable resource for researchers and clinicians to further identify molecular pathogenesis and therapeutic opportunities in iCCA.

Single-cell analysis of human primary prostate cancer reveals the heterogeneity of tumor-associated epithelial cell states

Prostate cancer is the second most common malignancy in men worldwide and consists of a mixture of tumor and non-tumor cell types. To characterize the prostate cancer tumor microenvironment, we perform single-cell RNA-sequencing on prostate biopsies, prostatectomy specimens, and patient-derived organoids from localized prostate cancer patients. We uncover heterogeneous cellular states in prostate epithelial cells marked by high androgen signaling states that are enriched in prostate cancer and identify a population of tumor-associated club cells that may be associated with prostate carcinogenesis. ERG-negative tumor cells, compared to ERG-positive cells, demonstrate shared heterogeneity with surrounding luminal epithelial cells and appear to give rise to common tumor microenvironment responses. Finally, we show that prostate epithelial organoids harbor tumor-associated epithelial cell states and are enriched with distinct cell types and states from their parent tissues. Our results provide diagnostically relevant insights and advance our understanding of the cellular states associated with prostate carcinogenesis.

JAFFAL: detecting fusion genes with long-read transcriptome sequencing

In cancer, fusions are important diagnostic markers and targets for therapy. Long-read transcriptome sequencing allows the discovery of fusions with their full-length isoform structure. However, due to higher sequencing error rates, fusion finding algorithms designed for short reads do not work. Here we present JAFFAL, to identify fusions from long-read transcriptome sequencing. We validate JAFFAL using simulations, cell lines, and patient data from Nanopore and PacBio. We apply JAFFAL to single-cell data and find fusions spanning three genes demonstrating transcripts detected from complex rearrangements. JAFFAL is available at https://github.com/Oshlack/JAFFA/wiki .

Frankly Invasive Carcinoma Ex-intraductal Carcinoma: Expanding on an Emerging and Perplexing Concept in Salivary Gland Tumor Pathology

Intraductal carcinoma (IDC) of the salivary glands is an uncommon and enigmatic tumor, our understanding of which is rapidly evolving. Recent studies have demonstrated multiple IDC subtypes and consistent gene fusions, most frequently involving RET. Because IDC is a ductal proliferation surrounded by flattened myoepithelial cells, it was previously presumed to be analogous to breast ductal carcinoma in situ, but recent evidence has shown that the myoepithelial cells of fusion-positive IDC harbor the same genetic alterations of the ductal cells and are therefore neoplastic. In addition, there are rare reports of fusion-positive IDC with overt areas of irregular invasion lacking myoepithelial cells, but this phenomenon is not well documented or understood. This study aims to better characterize these frankly invasive carcinoma ex-IDC. All cases of frankly invasive carcinoma ex-IDC were obtained from the authors' files. Inclusion criteria included a component of concurrent or antecedent IDC and/or a fusion known to be associated with IDC. Immunohistochemistry (S100, SOX10, mammaglobin, androgen receptor, p63, p40) and molecular analysis (targeted RNA sequencing or large panel DNA next generation sequencing) was performed. Clinical follow-up was obtained from medical records. Ten cases of frankly invasive carcinoma ex-IDC were identified. The tumors occurred in 8 men and 2 women ranging from 33 to 82 years (mean, 66.3). All but one case arose in the parotid gland. In 4 cases, the IDC component was intercalated duct type. It was mixed apocrine/intercalated duct in two, and in the remaining 4 cases, no residual IDC was identified. The frankly invasive carcinomas were remarkably heterogeneous, ranging from minimally to widely invasive beyond the confines of the IDC, low-grade to high-grade, with morphologies that varied from duct-forming to those having clear cell or sarcomatoid features, to frankly apocrine. The original diagnoses for these cases were (adeno) carcinoma, not otherwise specified (n = 6), salivary duct carcinoma (n = 3), and secretory carcinoma (n = 1). All cases harbored fusions: NCOA4::RET (n = 6), TRIM33::RET (n = 2), TRIM27::RET (n = 1), and STRN::ALK (n = 1). Clinically, one tumor recurred locally, cervical lymph node metastases occurred in five patients, and distant metastasis later developed in four of these patients. Our findings highlight striking diversity in frankly invasive carcinomas that arise from fusion-positive IDC, a tumor which may serve as a precursor neoplasm like pleomorphic adenoma. These carcinomas vary in their extent of invasion, grade, histologic appearances, and clinical behavior. Importantly, in contrast to pure IDC, which is believed to be indolent, many frankly invasive cases were aggressive. Because RET and ALK fusions are targetable, it is important to recognize the broad spectrum of frankly invasive carcinomas that can arise from IDC, particularly because some cases are completely overrun or recur without any recognizable IDC component. These results suggest fusion analysis may be of clinical benefit on any salivary gland (adeno) carcinoma, not otherwise specified or salivary duct carcinoma.

Fusion gene detection by RNA-sequencing complements diagnostics of acute myeloid leukemia and identifies recurring NRIP1-MIR99AHG rearrangements

Identification of fusion genes in clinical routine is mostly based on cytogenetics and targeted molecular genetics, such as metaphase karyotyping, fluorescence in situ hybridization and reverse-transcriptase polymerase chain reaction. However, sequencing technologies are becoming more important in clinical routine as processing time and costs per sample decrease. To evaluate the performance of fusion gene detection by RNAsequencing compared to standard diagnostic techniques, we analyzed 806 RNA-sequencing samples from patients with acute myeloid leukemia using two state-of-the-art software tools, namely Arriba and FusionCatcher. RNA-sequencing detected 90% of fusion events that were reported by routine with high evidence, while samples in which RNA-sequencing failed to detect fusion genes had overall lower and inhomogeneous sequence coverage. Based on properties of known and unknown fusion events, we developed a workflow with integrated filtering strategies for the identification of robust fusion gene candidates by RNA-sequencing. Thereby, we detected known recurrent fusion events in 26 cases that were not reported by routine and found discrepancies in evidence for known fusion events between routine and RNA-sequencing in three cases. Moreover, we identified 157 fusion genes as novel robust candidates and comparison to entries from ChimerDB or Mitelman Database showed novel recurrence of fusion genes in 14 cases. Finally, we detected the novel recurrent fusion gene NRIP1- MIR99AHG resulting from inv(21)(q11.2;q21.1) in nine patients (1.1%) and LTN1-MX1 resulting from inv(21)(q21.3;q22.3) in two patients (0.25%). We demonstrated that NRIP1-MIR99AHG results in overexpression of the 3' region of MIR99AHG and the disruption of the tricistronic miRNA cluster miR-99a/let-7c/miR-125b-2. Interestingly, upregulation of MIR99AHG and deregulation of the miRNA cluster, residing in the MIR99AHG locus, are known mechanisms of leukemogenesis in acute megakaryoblastic leukemia. Our findings demonstrate that RNA-sequencing has a strong potential to improve the systematic detection of fusion genes in clinical applications and provides a valuable tool for fusion discovery.

Targeted RNA-sequencing analysis for fusion transcripts detection in tumor diagnostics: assessment of bioinformatic tools reliability in FFPE samples

Aim:

Diagnostic laboratories are progressively introducing next-generation sequencing (NGS) technologies in the routine workflow to meet the increasing clinical need for comprehensive molecular characterization in cancer patients for diagnosis and precision medicine, including fusion-transcripts detection. Nevertheless, the low quality of messenger RNA (mRNA) extracted from formalin-fixed paraffin-embedded (FFPE) samples may affect the transition from traditional single-gene testing approaches [like fluorescence in situ hybridization (FISH), immunohistochemistry (IHC), or polymerase chain reaction (PCR)] to NGS. The present study is aimed at assessing the overall accuracy of RNA fusion transcripts detection by NGS analysis in FFPE samples in real-world diagnostics.

Methods:

Herein, NGS data from 190 soft tissue tumors (STTs) and carcinoma cases, discussed in the context of the institutional Molecular Tumor Board, are reported and analyzed by FusionPlex^© Solid tumor kit through the manufacturer’s pipeline and by two well-known fast and accurate open-source tools [Arriba (ARR) and spliced transcripts alignment to reference (STAR)-fusion (SFU)].

Results:

The combination of FusionPlex^© Solid tumor with ArcherDX^® Analysis suite (ADx) analysis package has been proven to be sensitive and specific in STT samples, while partial loss of sensitivity has been found in carcinoma specimens.

Conclusions:

Albeit ARR and SFU showed lower sensitivity, the use of additional fusion-detection tools can contribute to reinforcing or extending the output obtained by ADx, particularly in the case of low-quality input data. Overall, our results sustain the clinical use of NGS for the detection of fusion transcripts in FFPE material.

ChiTaH: a fast and accurate tool for identifying known human chimeric sequences from high-throughput sequencing data

Fusion genes or chimeras typically comprise sequences from two different genes. The chimeric RNAs of such joined sequences often serve as cancer drivers. Identifying such driver fusions in a given cancer or complex disease is important for diagnosis and treatment. The advent of next-generation sequencing technologies, such as DNA-Seq or RNA-Seq, together with the development of suitable computational tools, has made the global identification of chimeras in tumors possible. However, the testing of over 20 computational methods showed these to be limited in terms of chimera prediction sensitivity, specificity, and accurate quantification of junction reads. These shortcomings motivated us to develop the first ‘reference-based’ approach termed ChiTaH (Chimeric Transcripts from High–throughput sequencing data). ChiTaH uses 43,466 non–redundant known human chimeras as a reference database to map sequencing reads and to accurately identify chimeric reads. We benchmarked ChiTaH and four other methods to identify human chimeras, leveraging both simulated and real sequencing datasets. ChiTaH was found to be the most accurate and fastest method for identifying known human chimeras from simulated and sequencing datasets. Moreover, especially ChiTaH uncovered heterogeneity of the BCR-ABL1 chimera in both bulk and single-cells of the K-562 cell line, which was confirmed experimentally.

Discovery of clinically relevant fusions in pediatric cancer

Background

Pediatric cancers typically have a distinct genomic landscape when compared to adult cancers and frequently carry somatic gene fusion events that alter gene expression and drive tumorigenesis. Sensitive and specific detection of gene fusions through the analysis of next-generation-based RNA sequencing (RNA-Seq) data is computationally challenging and may be confounded by low tumor cellularity or underlying genomic complexity. Furthermore, numerous computational tools are available to identify fusions from supporting RNA-Seq reads, yet each algorithm demonstrates unique variability in sensitivity and precision, and no clearly superior approach currently exists. To overcome these challenges, we have developed an ensemble fusion calling approach to increase the accuracy of identifying fusions.

Results

Our Ensemble Fusion (EnFusion) approach utilizes seven fusion calling algorithms: Arriba, CICERO, FusionMap, FusionCatcher, JAFFA, MapSplice, and STAR-Fusion, which are packaged as a fully automated pipeline using Docker and Amazon Web Services (AWS) serverless technology. This method uses paired end RNA-Seq sequence reads as input, and the output from each algorithm is examined to identify fusions detected by a consensus of at least three algorithms. These consensus fusion results are filtered by comparison to an internal database to remove likely artifactual fusions occurring at high frequencies in our internal cohort, while a “known fusion list” prevents failure to report known pathogenic events. We have employed the EnFusion pipeline on RNA-Seq data from 229 patients with pediatric cancer or blood disorders studied under an IRB-approved protocol. The samples consist of 138 central nervous system tumors, 73 solid tumors, and 18 hematologic malignancies or disorders. The combination of an ensemble fusion-calling pipeline and a knowledge-based filtering strategy identified 67 clinically relevant fusions among our cohort (diagnostic yield of 29.3%), including RBPMS-MET, BCAN-NTRK1, and TRIM22-BRAF fusions. Following clinical confirmation and reporting in the patient’s medical record, both known and novel fusions provided medically meaningful information.

Conclusions

The EnFusion pipeline offers a streamlined approach to discover fusions in cancer, at higher levels of sensitivity and accuracy than single algorithm methods. Furthermore, this method accurately identifies driver fusions in pediatric cancer, providing clinical impact by contributing evidence to diagnosis and, when appropriate, indicating targeted therapies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08094-z.

Computational cancer neoantigen prediction: current status and recent advances

Over the last few decades, immunotherapy has shown significant therapeutic efficacy in a broad range of cancer types. Antitumor immune responses are contingent on the recognition of tumor-specific antigens, which are termed neoantigens. Tumor neoantigens are ideal targets for immunotherapy since they can be recognized as non-self antigens by the host immune system and thus are able to elicit an antitumor T-cell response. There are an increasing number of studies that highlight the importance of tumor neoantigens in immunoediting and in the sensitivity to immune checkpoint blockade. Therefore, one of the most fundamental tasks in the field of immuno-oncology research is the identification of patient-specific neoantigens. To this end, a plethora of computational approaches have been developed in order to predict tumor-specific aberrant peptides and quantify their likelihood of binding to patients' human leukocyte antigen molecules in order to be recognized by T cells. In this review, we systematically summarize and present the most recent advances in computational neoantigen prediction, and discuss the challenges and novel methods that are being developed to resolve them.

•

Tumors have the ability to acquire immune escape mechanisms.

•

Tumor-specific aberrant peptides (neoantigens) can elicit an immune response by the host immune system.

•

The identification of neoantigens is one of the most fundamental tasks in the field of immuno-oncology research.

•

A plethora of computational approaches have been developed in order to predict patient-specificneoantigens.

Microsecretory Adenocarcinoma of Salivary Glands: An Expanded Series of 24 Cases

Microsecretory adenocarcinoma (MSA) is a recently described salivary gland tumor with a characteristic histologic and immunophenotypic profile and recurrent MEF2C-SS18 fusions. Because only six cases of MSA have been published, its complete clinicopathologic spectrum is unclear, and its biologic behavior has not been documented. Here, we present an updated and expanded experience of 24 MSA cases. All cases of MSA were obtained from the authors' files. Immunohistochemistry for S100, SOX10, p63, p40, SMA, calponin, and mammaglobin was performed. Molecular analysis was performed by targeted RNA sequencing, SS18 break apart fluorescence in situ hybridization, and/or reverse transcriptase polymerase chain reaction for MEF2C-SS18 fusion. Clinical follow-up was obtained from medical records. A total of 24 MSA cases were collected, from 13 women and 11 men, ranging from 17 to 83 years (mean 49.5 years). The vast majority (23 of 24) arose in the oral cavity, with the palate (n = 14) and buccal mucosa (n = 6) as the most frequent subsites. Tumors showed consistent histologic features including: (1) microcystic tubules, (2) flattened intercalated duct-like cells, (3) monotonous oval hyperchromatic nuclei, (4) abundant basophilic luminal secretions, (5) fibromyxoid stroma, and (6) circumscribed borders with subtle infiltration. The tumors were very consistently positive for S100 (24 of 24), p63 (24 of 24), and SOX10 (14 of 14) and negative for p40 (0 of 21), calponin (0 of 12) and mammaglobin (0 of 16), while SMA (4 of 20) was variable. MEF2C-SS18 fusion was demonstrated in 21 of 24 cases; in the remaining 3 cases with insufficient RNA, SS18 break apart FISH was positive. Treatment information was available in 17 cases, all of which were managed with surgery only. In 14 cases with follow-up (1-216 months, mean 30), no cases recurred or metastasized. MSA is a distinct salivary gland neoplasm with remarkably consistent clinical, histologic, immunophenotypic, and genetic features that generally behaves in an indolent manner following surgery alone. These observations solidify MSA as a unique, low-grade salivary gland carcinoma that warrants inclusion in the next version of the WHO classification of head and neck tumors.

ecDNA hubs drive cooperative intermolecular oncogene expression

Extrachromosomal DNA (ecDNA) is prevalent in human cancers and mediates high expression of oncogenes through gene amplification and altered gene regulation1. Gene induction typically involves cis-regulatory elements that contact and activate genes on the same chromosome2,3. Here we show that ecDNA hubs-clusters of around 10-100 ecDNAs within the nucleus-enable intermolecular enhancer-gene interactions to promote oncogene overexpression. ecDNAs that encode multiple distinct oncogenes form hubs in diverse cancer cell types and primary tumours. Each ecDNA is more likely to transcribe the oncogene when spatially clustered with additional ecDNAs. ecDNA hubs are tethered by the bromodomain and extraterminal domain (BET) protein BRD4 in a MYC-amplified colorectal cancer cell line. The BET inhibitor JQ1 disperses ecDNA hubs and preferentially inhibits ecDNA-derived-oncogene transcription. The BRD4-bound PVT1 promoter is ectopically fused to MYC and duplicated in ecDNA, receiving promiscuous enhancer input to drive potent expression of MYC. Furthermore, the PVT1 promoter on an exogenous episome suffices to mediate gene activation in trans by ecDNA hubs in a JQ1-sensitive manner. Systematic silencing of ecDNA enhancers by CRISPR interference reveals intermolecular enhancer-gene activation among multiple oncogene loci that are amplified on distinct ecDNAs. Thus, protein-tethered ecDNA hubs enable intermolecular transcriptional regulation and may serve as units of oncogene function and cooperative evolution and as potential targets for cancer therapy.

The CLIP1-LTK fusion is an oncogenic driver in non-small-cell lung cancer

Lung cancer is one of the most aggressive tumour types. Targeted therapies stratified by oncogenic drivers have substantially improved therapeutic outcomes in patients with non-small-cell lung cancer (NSCLC)¹. However, such oncogenic drivers are not found in 25-40% of cases of lung adenocarcinoma, the most common histological subtype of NSCLC². Here we identify a novel fusion transcript of CLIP1 and LTK using whole-transcriptome sequencing in a multi-institutional genome screening platform (LC-SCRUM-Asia, UMIN000036871). The CLIP1-LTK fusion was present in 0.4% of NSCLCs and was mutually exclusive with other known oncogenic drivers. We show that kinase activity of the CLIP1-LTK fusion protein is constitutively activated and has transformation potential. Treatment of Ba/F3 cells expressing CLIP1-LTK with lorlatinib, an ALK inhibitor, inhibited CLIP1-LTK kinase activity, suppressed proliferation and induced apoptosis. One patient with NSCLC harbouring the CLIP1-LTK fusion showed a good clinical response to lorlatinib treatment. To our knowledge, this is the first description of LTK alterations with oncogenic activity in cancers. These results identify the CLIP1-LTK fusion as a target in NSCLC that could be treated with lorlatinib.

CLINICOPATHOLOGIC AND MOLECULAR ANALYSIS OF A SINGLE BCOR-CCNB3+ UNDIFFERENTIATED SARCOMA OF THE KIDNEY CONFERS SIGNIFICANT EPIGENETIC ALTERATIONS

Undifferentiated soft tissue sarcomas (UDSTSs) are a group of mesenchymal tumors that remain a diagnostic challenge because of their morphologic heterogeneity and unclear histologic origin (Peters et al., Mod Pathol28: 575 [2015]). In this case report, we present the first multiomics molecular signature for a BCOR–CCNB3 sarcoma (BCS) that includes mutation analysis, gene expression, DNA methylation, and micro RNA (miRNA) expression. We identify a paucity of additional mutations in this tumor and detail that there is significant dysregulation of gene expression of epigenetic remodeling agents including key members of the PRC, Sin3A/3b, NuRD, and NcoR/SMRT complexes and the DNA methyltransferases DNMT1, DNMT3a, and DNMT3b. This is accompanied by significant DNA methylation changes and dysregulation of multiple miRNAs with known links to tumorigenesis. This study significantly increases our understanding of the BCOR effects on fusion-positive undifferentiated sarcomas at both the genomic and epigenomic level and suggests that as better-tailored and more refined treatment algorithms continue to evolve, epigenetic modifying agents should be further evaluated for their efficacy against these tumors.

Comprehensive Analysis of Large-Scale Transcriptomes from Multiple Cancer Types

Various abnormalities of transcriptional regulation revealed by RNA sequencing (RNA-seq) have been reported in cancers. However, strategies to integrate multi-modal information from RNA-seq, which would help uncover more disease mechanisms, are still limited. Here, we present PipeOne, a cross-platform one-stop analysis workflow for large-scale transcriptome data. It was developed based on Nextflow, a reproducible workflow management system. PipeOne is composed of three modules, data processing and feature matrices construction, disease feature prioritization, and disease subtyping. It first integrates eight different tools to extract different information from RNA-seq data, and then used random forest algorithm to study and stratify patients according to evidences from multiple-modal information. Its application in five cancers (colon, liver, kidney, stomach, or thyroid; total samples n = 2024) identified various dysregulated key features (such as PVT1 expression and ABI3BP alternative splicing) and pathways (especially liver and kidney dysfunction) shared by multiple cancers. Furthermore, we demonstrated clinically-relevant patient subtypes in four of five cancers, with most subtypes characterized by distinct driver somatic mutations, such as TP53, TTN, BRAF, HRAS, MET, KMT2D, and KMT2C mutations. Importantly, these subtyping results were frequently contributed by dysregulated biological processes, such as ribosome biogenesis, RNA binding, and mitochondria functions. PipeOne is efficient and accurate in studying different cancer types to reveal the specificity and cross-cancer contributing factors of each cancer.It could be easily applied to other diseases and is available at GitHub.