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

Genomic and transcriptomic profiling of inflammatory breast cancer reveals distinct molecular characteristics to non-inflammatory breast cancers

PURPOSE

Inflammatory breast cancer (IBC), a rare and highly aggressive form of breast cancer, accounts for 10% of breast cancer-related deaths. Previous omics studies of IBC have focused solely on one of genomics or transcriptomics and did not discover common differences that could distinguish IBC from non-IBC.

METHODS

Seventeen IBC patients and five non-IBC patients as well as additional thirty-three Asian breast cancer samples from TCGA-BRCA were included for the study. We performed whole-exon sequencing (WES) to investigate different somatic genomic alterations, copy number variants, and large structural variants between IBC and non-IBC. Bulk RNA sequencing (RNA-seq) was performed to examine the differentially expressed genes, pathway enrichment, and gene fusions. WES and RNA-seq data were further investigated in combination to discover genes that were dysregulated in both genomics and transcriptomics.

RESULTS

Copy number variation analysis identified 10 cytobands that showed higher frequency in IBC. Structural variation analysis showed more frequent deletions in IBC. Pathway enrichment and immune infiltration analysis indicated increased immune activation in IBC samples. Gene fusions including CTSC-RAB38 were found to be more common in IBC. We demonstrated more commonly dysregulated RAS pathway in IBC according to both WES and RNA-seq. Inhibitors targeting RAS signaling and its downstream pathways were predicted to possess promising effects in IBC treatment.

CONCLUSION

We discovered differences unique in Asian women that could potentially explain IBC etiology and presented RAS signaling pathway as a potential therapeutic target in IBC treatment.

Characterization and integrated analysis of extrachromosomal DNA amplification in hematological malignancies

The study of extrachromosomal DNA (ecDNA), an element existing beyond classical chromosomes, contributes to creating a more comprehensive map of the cancer genome. In hematological malignancies, research on ecDNA has lacked comprehensive investigation into its frequency, structure, function, and mechanisms of formation. We re-analyzed WGS data from 208 hematological cancer samples across 11 types, focusing on ecDNA characteristics. Amplification of ecDNA was observed in 7 of these cancer types, with no instances found in normal blood cells. Patients with leukemia carrying ecDNA showed a low induction therapy remission rate (<30 %), a high relapse rate (75 %) among those who achieved complete remission, and a significantly lower survival rate compared to the general leukemia population, even those with complex chromosomal karyotypes. Among the 55 identified ecDNA amplicons, 268 genes were detected, of which 38 are known cancer-related genes exhibiting significantly increased copy numbers. By integrating RNA-Seq data, we discovered that the increased copy number, resulting in a higher amount of available DNA templates, indeed leads to the elevated expression of genes encoded on ecDNA. Additionally, through the integration of H3K4me3/H3K27ac chromatin immunoprecipitation sequencing, assay for transposase-accessible chromatin with sequencing, and high-throughput chromosome conformation capture data, we identified that ecDNA amplifications can also facilitate efficient, copy number-independent amplification of oncogenes. This process is linked to active histone modifications, improved chromatin accessibility, and enhancer hijacking, all of which are effects of ecDNA amplification. Mechanistically, chromothripsis and dysfunction of the DNA repair pathway can, to some extent, explain the origin of ecDNA.

Epigenetic alterations affecting hematopoietic regulatory networks as drivers of mixed myeloid/lymphoid leukemia

Leukemias with ambiguous lineage comprise several loosely defined entities, often without a clear mechanistic basis. Here, we extensively profile the epigenome and transcriptome of a subgroup of such leukemias with CpG Island Methylator Phenotype. These leukemias exhibit comparable hybrid myeloid/lymphoid epigenetic landscapes, yet heterogeneous genetic alterations, suggesting they are defined by their shared epigenetic profile rather than common genetic lesions. Gene expression enrichment reveals similarity with early T-cell precursor acute lymphoblastic leukemia and a lymphoid progenitor cell of origin. In line with this, integration of differential DNA methylation and gene expression shows widespread silencing of myeloid transcription factors. Moreover, binding sites for hematopoietic transcription factors, including CEBPA, SPI1 and LEF1, are uniquely inaccessible in these leukemias. Hypermethylation also results in loss of CTCF binding, accompanied by changes in chromatin interactions involving key transcription factors. In conclusion, epigenetic dysregulation, and not genetic lesions, explains the mixed phenotype of this group of leukemias with ambiguous lineage. The data collected here constitute a useful and comprehensive epigenomic reference for subsequent studies of acute myeloid leukemias, T-cell acute lymphoblastic leukemias and mixed-phenotype leukemias.

Real-World Clinical Utility of Targeted RNA Sequencing in Leukemia Diagnosis and Management

Gene fusions are key drivers in acute leukemia, impacting diagnosis and treatment decisions. We analyzed 264 leukemia patients using targeted RNA sequencing with conventional karyotyping and reverse transcription polymerase chain reaction (RT-PCR). Leukemic fusions were detected in 127 patients (48.1%). The new guidelines introduced additional diagnostic criteria, expanding the spectrum of gene fusions. We discovered three novel fusions (RUNX1::DOPEY2, RUNX1::MACROD2, and ZCCHC7::LRP1B). We analyzed recurrent breakpoints for the KMT2A and NUP98 rearrangements. Targeted RNA sequencing showed consistent results with RT-PCR in all tested samples. However, when compared to conventional karyotyping, we observed an 83.3% concordance rate, with 29 cases found only in targeted RNA sequencing, 7 cases with discordant results, and 5 cases found only in conventional karyotyping. For the five cases where known leukemic gene rearrangements were suspected only in conventional karyotyping, we conducted additional messenger RNA sequencing in four cases and proved no pathogenic gene rearrangements. Targeted RNA sequencing proved advantageous for the rapid and accurate interpretation of gene rearrangements. The concurrent use of multiple methods was essential for a comprehensive evaluation. Comprehensive molecular analysis enhances our understanding of leukemia's genetic basis, aiding diagnosis and classification. Advanced molecular techniques improve clinical decision-making, offering potential benefits.

Comprehensive molecular characterization of long-term glioblastoma survivors

Fewer than 5 % glioblastoma (GBM) patients survive over five years and are termed long-term survivors (LTS), yet their molecular background is unclear. The present cohort included 72 isocitrate dehydrogenase (IDH)-wildtype GBM patients, consisting of 35 LTS and 37 short-term survivors (STS), and we employed whole exome sequencing, RNA-seq and DNA methylation array to delineate this largest LTS cohort to date. Although LTS and STS demonstrated analogous clinical characters and classical GBM biomarkers, CASC5 (P = 0.002) and SPEN (P = 0.013) mutations were enriched in LTS, whereas gene-to-gene fusions were concentrated in STS (P = 0.007). Importantly, LTS exhibited higher tumor mutation burden (P < 0.001) and copy number (CN) increase (P = 0.013), but lower mutant-allele tumor heterogeneity score (P < 0.001) and CN decrease (P = 0.026). Additionally, LTS demonstrated hypermethylated genome (P < 0.001) relative to STS. Differentially expressed and methylated genes both enriched in olfactory transduction. Further, analysis of the tumor microenvironment revealed higher infiltration of M1 macrophages (P = 0.043), B cells (P = 0.016), class-switched memory B cells (P = 0.002), central memory CD4+ T cells (P = 0.031) and CD4+ Th1 cells (P = 0.005) in LTS. We also separately analyzed a subset of patients who were methylation class-defined GBM, contributing 70.8 % of the entire cohort, and obtained similar results relative to prior analyses. Finally, we demonstrated that LTS and STS could be distinguished using a subset of molecular features. Taken together, the present study delineated unique molecular attributes of LTS GBM.

Passive Smoking-Induced Mutagenesis as a Promoter of Lung Carcinogenesis

INTRODUCTION

The International Agency for Research on Cancer has classified passive smoking (PS) or secondhand smoke exposure as a group 1 carcinogen linked to lung cancer. However, in contrast to active smoking, the mutagenic properties of PS remain unclear.

METHODS

A consecutive cohort of 564 lung adenocarcinoma samples from female never-smokers, who provided detailed information about their exposure to PS during adolescence and in their thirties through a questionnaire, was prepared. Of these, all 291 cases for whom frozen tumor tissues were available were subjected to whole exome sequencing to estimate tumor mutational burden, and the top 84 cases who were exposed daily, or not, to PS during adolescence, in their thirties or in both periods, were further subjected to whole genome sequencing.

RESULTS

A modest yet statistically significant increase in tumor mutational burden was observed in the group exposed to PS compared with the group not exposed to PS (median values = 1.44 versus 1.29 per megabase, respectively; p = 0.020). Instead of inducing driver oncogene mutations, PS-induced substantial subclonal mutations exhibiting APOBEC-type signatures, including SMAD4 and ADGRG6 hotspot mutations. A polymorphic APOBEC3A/3B allele-specific to the Asian population that leads to up-regulated expression of APOBEC3A accentuated the mutational load in individuals exposed daily to PS during adolescence.

CONCLUSIONS

This study reveals that PS-induced mutagenesis can promote lung carcinogenesis. The APOBEC3A/3B polymorphism may serve as a biomarker for identifying passive nonsmoking individuals at high risk of developing lung cancer.

Brd4::Nutm1 fusion gene initiates NUT carcinoma in vivo

NUT carcinoma (NC) is an aggressive cancer with no effective treatment. About 70% of NUT carcinoma is associated with chromosome translocation events that lead to the formation of a BRD4::NUTM1 fusion gene. Because the BRD4::NUTM1 gene is unequivocally cytotoxic when ectopically expressed in cell lines, questions remain on whether the fusion gene can initiate NC. Here, we report the first genetically engineered mouse model for NUT carcinoma that recapitulates the human t(15;19) chromosome translocation in mice. We demonstrated that the mouse t(2;17) syntenic chromosome translocation, forming the Brd4::Nutm1 fusion gene, could induce aggressive carcinomas in mice. The tumors present histopathological and molecular features similar to human NC, with enrichment of undifferentiated cells. Similar to the reports of human NC incidence, Brd4::Nutm1 can induce NC from a broad range of tissues with a strong phenotypical variability. The consistent induction of poorly differentiated carcinoma demonstrated a strong reprogramming activity of BRD4::NUTM1. The new mouse model provided a critical preclinical model for NC that will lead to better understanding and therapy development for NC.

The Variable Genomic Landscape During Osteosarcoma Progression: Insights From a Longitudinal WGS Analysis

Osteosarcoma is a primary bone tumor that exhibits a complex genomic landscape characterized by gross chromosomal abnormalities. Osteosarcoma patients often develop metastatic disease, resulting in limited therapeutic options and poor survival rates. To gain knowledge on the mechanisms underlying osteosarcoma heterogeneity and metastatic process, it is important to obtain a detailed profile of the genomic alterations that accompany osteosarcoma progression. We performed WGS on multiple tissue samples from six patients with osteosarcoma, including the treatment naïve biopsy of the primary tumor, resection of the primary tumor after neoadjuvant chemotherapy, local recurrence, and distant metastases. A comprehensive analysis of single-nucleotide variants (SNVs), structural variants, copy number alterations (CNAs), and chromothripsis events revealed the genomic heterogeneity during osteosarcoma progression. SNVs and structural variants were found to accumulate over time, contributing to an increased complexity of the genome of osteosarcoma during disease progression. Phylogenetic trees based on SNVs and structural variants reveal distinct evolutionary patterns between patients, including linear, neutral, and branched patterns. The majority of osteosarcomas showed variable copy number profiles or gained whole-genome doubling in later occurrences. Large proportions of the genome were affected by loss of heterozygosity (LOH), although these regions remain stable during progression. Additionally, chromothripsis is not confined to a single early event, as multiple other chromothripsis events may appear in later occurrences. Together, we provide a detailed analysis of the complex genome of osteosarcomas and show that five of six osteosarcoma genomes are highly dynamic and variable during progression.

Novel gene fusions in human oropharyngeal carcinoma

Few reports have analyzed the fusion genes involved in carcinogenesis in the oropharynx, where the incidence of human papillomavirus-associated tumors is relatively low. The aim of this study was to identify novel driver fusion genes in patients with oropharyngeal cancer. The study enrolled fifty-seven patients who were diagnosed with oropharyngeal carcinoma. RNA sequencing data from fresh-frozen specimens were used to identify candidate fusion genes via the JAFFA, arriba, and STAR-Fusion pipelines. Candidate fusion genes were confirmed by direct sequencing. The expression level of a candidate fusion gene was compared to that of tumors without fusion genes. Finally, filtering was performed for driver genes using the annoFuse pipeline. In addition, the VIRTUS pipeline was used to analyze the presence of human papillomavirus in the tumors. We identified 5 (8.8 %) novel potential driver in-frame fusion genes, MKNK2::MOB3A, ICMT::RPS6KA3, ATP1B3::GRK7, CSNK2A1::KIF16B, and FGFR3::MAEA, and 1 (1.8 %) known in-frame fusion gene, FGFR3::TACC3, in 57 patients with pharyngeal carcinoma. Our results suggest that sporadic fusion genes may contribute to tumorigenesis in oropharyngeal carcinomas.

Sarcoma_CellminerCDB: A tool to interrogate the genomic and functional characteristics of a comprehensive collection of sarcoma cell lines

Sarcomas are a diverse group of rare malignancies composed of multiple different clinical and molecular subtypes. Due to their rarity and heterogeneity, basic, translational, and clinical research in sarcoma has trailed behind that of other cancers. Outcomes for patients remain generally poor due to an incomplete understanding of disease biology and a lack of novel therapies. To address some of the limitations impeding preclinical sarcoma research, we have developed Sarcoma_CellMinerCDB, a publicly available interactive tool that merges publicly available sarcoma cell line data and newly generated omics data to create a comprehensive database of genomic, transcriptomic, methylomic, proteomic, metabolic, and pharmacologic data on 133 annotated sarcoma cell lines. The reproducibility, functionality, biological relevance, and therapeutic applications of Sarcoma_CellMinerCDB described herein are powerful tools to address and generate biological questions and test hypotheses for translational research. Sarcoma_CellMinerCDB (https://discover.nci.nih.gov/SarcomaCellMinerCDB) aims to contribute to advancing the preclinical study of sarcoma.

Oncogene-Driven Non-Small Cell Lung Cancers in Patients with a History of Smoking Lack Smoking-Induced Mutations

Non-small cell lung cancers (NSCLC) in nonsmokers are mostly driven by mutations in the oncogenes EGFR, ERBB2, and MET and fusions involving ALK and RET. In addition to occurring in nonsmokers, alterations in these "nonsmoking-related oncogenes" (NSRO) also occur in smokers. To better understand the clonal architecture and genomic landscape of NSRO-driven tumors in smokers compared with typical-smoking NSCLCs, we investigated genomic and transcriptomic alterations in 173 tumor sectors from 48 NSCLC patients. NSRO-driven NSCLCs in smokers and nonsmokers had similar genomic landscapes. Surprisingly, even in patients with prominent smoking histories, the mutational signature caused by tobacco smoking was essentially absent in NSRO-driven NSCLCs, which was confirmed in two large NSCLC data sets from other geographic regions. However, NSRO-driven NSCLCs in smokers had higher transcriptomic activities related to the regulation of the cell cycle. These findings suggest that, whereas the genomic landscape is similar between NSRO-driven NSCLC in smokers and nonsmokers, smoking still affects the tumor phenotype independently of genomic alterations.

SIGNIFICANCE

Non-small cell lung cancers driven by nonsmoking-related oncogenes do not harbor genomic scars caused by smoking regardless of smoking history, indicating that the impact of smoking on these tumors is mainly nongenomic.

Myeloid sarcoma with plasmacytoid dendritic cell-like proliferation associated with IKZF1, ETV6 and DNMT3A mutations

The classification of clonal plasmacytoid dendritic cell (pDC) proliferation associated with myeloid neoplasms remains a topic of ongoing debate. Although the fifth edition of the World Health Organization classification classifies clonal pDC proliferation into two categories, it is unclear whether this classification adequately captures the complexities of clonal pDC pathogenesis. We present a clinical case featuring myeloid sarcoma with pDC-like cells in cervical lymph nodes and bone marrow (BM). Analysis of biopsy specimens and BM aspirate revealed two distinct cellular populations expressing myeloid and pDC markers. One population exhibited myeloid leukemia and monocyte markers, including MPO, CD13, CD33, CD11b, and CD14, while the other manifested an immunophenotype reminiscent of pDCs, characterized by expression of CD56 and CD123. Additionally, whole exome sequencing and RNA sequencing of BM mononuclear cells were conducted to explore the pathophysiology of this rare malignancy, and unveiled pDC-like cell proliferation driven by IKZF1 and ETV6 mutations originating from clonal hematopoiesis initiated by a DNMT3A mutation. Notably, venetoclax-based therapy exhibited efficacy for achieving and sustaining complete remission. This case provides pivotal insights into the mechanistic aspects of pDC/pDC-like cell proliferation in myeloid sarcoma, offering valuable perspectives on therapeutic strategies.

A regionally based precision medicine implementation initiative in North Africa:The PerMediNA consortium

Precision Medicine is being increasingly used in the developed world to improve health care. While several Precision Medicine (PM) initiatives have been launched worldwide, their implementations have proven to be more challenging particularly in low- and middle-income countries. To address this issue, the "Personalized Medicine in North Africa" initiative (PerMediNA) was launched in three North African countries namely Tunisia, Algeria and Morocco. PerMediNA is coordinated by Institut Pasteur de Tunis together with the French Ministry for Europe and Foreign Affairs, with the support of Institut Pasteur in France. The project is carried out along with Institut Pasteur d'Algérie and Institut Pasteur du Maroc in collaboration with national and international leading institutions in the field of PM including Institut Gustave Roussy in Paris. PerMediNA aims to assess the readiness level of PM implementation in North Africa, to strengthen PM infrastructure, to provide workforce training, to generate genomic data on North African populations, to implement cost effective, affordable and sustainable genetic testing for cancer patients and to inform policy makers on how to translate research knowledge into health products and services. Gender equity and involvement of young scientists in this implementation process are other key goals of the PerMediNA project. In this paper, we are describing PerMediNA as the first PM implementation initiative in North Africa. Such initiatives contribute significantly in shortening existing health disparities and inequities between developed and developing countries and accelerate access to innovative treatments for global health.

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.

Aggressive systemic mastocytosis with the co-occurrence of PRKG2::PDGFRB, KAT6A::NCOA2, and RXRA::NOTCH1 fusion transcripts and a heterozygous RUNX1 frameshift mutation

Systemic mastocytosis (SM) is a myeloproliferative neoplasm displaying abnormal mast cell proliferation. It is subdivided into different forms, including aggressive systemic mastocytosis (ASM) and systemic mastocytosis with an associated hematologic neoplasm (SM-AHN). Oncogenic genetic alterations include point mutations, mainly the KIT D816V, conferring poor prognosis and therapy resistance, and fusion genes, with those involving PDGFRA/PDGFRB as the most recurrent events. We here describe an ASM case negative to the KIT D816V and JAK2 V617F alterations but showing a RUNX1 frameshift heterozygous mutation and the co-occurrence of three fusion transcripts. The first one, PRKG2::PDGFRB, was generated by a balanced t(4;5)(q24;q32) translocation as the sole abnormality. Other two novel chimeras, KAT6A::NCOA2 and RXRA::NOTCH1, originated from cryptic intra-chromosomal abnormalities. The patient rapidly evolved towards SM-AHN, characterized by the persistence of the PRKG2::PDGFRB chimera, due to the presence of an extra copy of the der(5)t(4;5)(q24;q34) chromosome and an increase in the RUNX1 mutation allelic frequency. The results indicated that the transcriptional landscape and the mutational profile of SM deserve attention to predict the evolution and prognosis of this complex disease, whose classification criteria are still a matter of debate.

Prediction of tumor-specific splicing from somatic mutations as a source of neoantigen candidates

Motivation

Neoantigens are promising targets for cancer immunotherapies and might arise from alternative splicing. However, detecting tumor-specific splicing is challenging because many non-canonical splice junctions identified in tumors also appear in healthy tissues. To increase tumor-specificity, we focused on splicing caused by somatic mutations as a source for neoantigen candidates in individual patients.

Results

We developed the tool splice2neo with multiple functionalities to integrate predicted splice effects from somatic mutations with splice junctions detected in tumor RNA-seq and to annotate the resulting transcript and peptide sequences. Additionally, we provide the tool EasyQuant for targeted RNA-seq read mapping to candidate splice junctions. Using a stringent detection rule, we predicted 1.7 splice junctions per patient as splice targets with a false discovery rate below 5% in a melanoma cohort. We confirmed tumor-specificity using independent, healthy tissue samples. Furthermore, using tumor-derived RNA, we confirmed individual exon-skipping events experimentally. Most target splice junctions encoded neoepitope candidates with predicted major histocompatibility complex (MHC)-I or MHC-II binding. Compared to neoepitope candidates from non-synonymous point mutations, the splicing-derived MHC-I neoepitope candidates had lower self-similarity to corresponding wild-type peptides. In conclusion, we demonstrate that identifying mutation-derived, tumor-specific splice junctions can lead to additional neoantigen candidates to expand the target repertoire for cancer immunotherapies.

Availability and implementation

The R package splice2neo and the python package EasyQuant are available at https://github.com/TRON-Bioinformatics/splice2neo and https://github.com/TRON-Bioinformatics/easyquant, respectively.

Synovial Sarcoma Chromatin Dynamics Reveal a Continuum in SS18:SSX Reprograming

Synovial sarcoma (SyS) is an aggressive soft-tissue malignancy characterized by a pathognomonic chromosomal translocation leading to the formation of the SS18::SSX fusion oncoprotein. SS18::SSX associates with mammalian BAF complexes suggesting deregulation of chromatin architecture as the oncogenic driver in this tumour type. To examine the epigenomic state of SyS we performed comprehensive multi-omics analysis on 52 primary pre-treatment human SyS tumours. Our analysis revealed a continuum of epigenomic states across the cohort at fusion target genes independent of rare somatic genetic lesions. We identify cell-of-origin signatures defined by enhancer states and reveal unexpected relationships between H2AK119Ub1 and active marks. The number of bivalent promoters, dually marked by the repressive H3K27me3 and activating H3K4me3 marks, has strong prognostic value and outperforms tumor grade in predicting patient outcome. Finally, we identify SyS defining epigenomic features including H3K4me3 expansion associated with striking promoter DNA hypomethylation in which SyS displays the lowest mean methylation level of any sarcoma subtype. We explore these distinctive features as potential vulnerabilities in SyS and identify H3K4me3 inhibition as a promising therapeutic strategy.

Proteogenomic Characterization Reveals Estrogen Signaling as a Target for Never-Smoker Lung Adenocarcinoma Patients without EGFR or ALK Alterations

Never-smoker lung adenocarcinoma (NSLA) is prevalent in Asian populations, particularly in women. EGFR mutations and anaplastic lymphoma kinase (ALK) fusions are major genetic alterations observed in NSLA, and NSLA with these alterations have been well studied and can be treated with targeted therapies. To provide insights into the molecular profile of NSLA without EGFR and ALK alterations (NENA), we selected 141 NSLA tissues and performed proteogenomic characterization, including whole genome sequencing (WGS), transcriptomic, methylation EPIC array, total proteomic, and phosphoproteomic analyses. Forty patients with NSLA harboring EGFR and ALK alterations and seven patients with NENA with microsatellite instability were excluded. Genome analysis revealed that TP53 (25%), KRAS (22%), and SETD2 (11%) mutations and ROS1 fusions (14%) were the most frequent genetic alterations in NENA patients. Proteogenomic impact analysis revealed that STK11 and ERBB2 somatic mutations had broad effects on cancer-associated genes in NENA. DNA copy number alteration analysis identified 22 prognostic proteins that influenced transcriptomic and proteomic changes. Gene set enrichment analysis revealed estrogen signaling as the key pathway activated in NENA. Increased estrogen signaling was associated with proteogenomic alterations, such as copy number deletions in chromosomes 14 and 21, STK11 mutation, and DNA hypomethylation of LLGL2 and ST14. Finally, saracatinib, an Src inhibitor, was identified as a potential drug for targeting activated estrogen signaling in NENA and was experimentally validated in vitro. Collectively, this study enhanced our understanding of NENA NSLA by elucidating the proteogenomic landscape and proposed saracatinib as a potential treatment for this patient population that lacks effective targeted therapies.

SIGNIFICANCE

The proteogenomic landscape in never-smoker lung cancer without known driver mutations reveals prognostic proteins and enhanced estrogen signaling that can be targeted as a potential therapeutic strategy to improve patient outcomes.

RT-PCR assay to detect FGFR3::TACC3 fusions in formalin-fixed, paraffin-embedded glioblastoma samples

Background

One targeted treatment option for isocitrate dehydrogenase (IDH)-wild-type glioblastoma focuses on tumors with fibroblast growth factor receptor 3::transforming acidic coiled-coil-containing protein 3 (FGFR3::TACC3) fusions. FGFR3::TACC3 fusion detection can be challenging, as targeted RNA next-generation sequencing (NGS) is not routinely performed, and immunohistochemistry is an imperfect surrogate marker. Fusion status can be determined using reverse transcription polymerase chain reaction (RT-PCR) on fresh frozen (FF) material, but sometimes only formalin-fixed, paraffin-embedded (FFPE) tissue is available.

Aim

To develop an RT-PCR assay to determine FGFR3::TACC3 status in FFPE glioblastoma samples.

Methods

Twelve tissue microarrays with 353 historical glioblastoma samples were immunohistochemically stained for FGFR3. Samples with overexpression of FGFR3 (n = 13) were subjected to FGFR3::TACC3 RT-PCR on FFPE, using 5 primer sets for the detection of 5 common fusion variants. Fusion-negative samples were additionally analyzed with NGS (n = 6), FGFR3 Fluorescence In Situ Hybridization (n = 6), and RNA sequencing (n = 5).

Results

Using RT-PCR on FFPE material of the 13 samples with FGFR3 overexpression, we detected an FGFR3::TACC3 fusion in 7 samples, covering 3 different fusion variants. For 5 of these FF was available, and the presence of the fusion was confirmed through RT-PCR on FF. With RNA sequencing, 1 additional sample was found to harbor an FGFR3::TACC3 fusion (variant not covered by current RT-PCR for FFPE). The frequency of FGFR3::TACC3 fusion in this cohort was 9/353 (2.5%).

Conclusions

RT-PCR for FGFR3::TACC3 fusions can successfully be performed on FFPE material, with a specificity of 100% and (due to limited primer sets) a sensitivity of 83.3%. This assay allows for the identification of potential targeted treatment options when only formalin-fixed tissue is available.

Central neurocytoma exhibits radial glial cell signatures with FGFR3 hypomethylation and overexpression

We explored the genomic events underlying central neurocytoma (CN), a rare neoplasm of the central nervous system, via multiomics approaches, including whole-exome sequencing, bulk and single-nuclei RNA sequencing, and methylation sequencing. We identified FGFR3 hypomethylation leading to FGFR3 overexpression as a major event in the ontogeny of CN that affects crucial downstream events, such as aberrant PI3K-AKT activity and neuronal development pathways. Furthermore, we found similarities between CN and radial glial cells based on analyses of gene markers and CN tumor cells and postulate that CN tumorigenesis is due to dysregulation of radial glial cell differentiation into neurons. Our data demonstrate the potential role of FGFR3 as one of the leading drivers of tumorigenesis in CN.

Integrated multiomic profiling of breast cancer in the Chinese population reveals patient stratification and therapeutic vulnerabilities

Molecular profiling guides precision treatment of breast cancer; however, Asian patients are underrepresented in publicly available large-scale studies. We established a comprehensive multiomics cohort of 773 Chinese patients with breast cancer and systematically analyzed their genomic, transcriptomic, proteomic, metabolomic, radiomic and digital pathology characteristics. Here we show that compared to breast cancers in white individuals, Asian individuals had more targetable AKT1 mutations. Integrated analysis revealed a higher proportion of HER2-enriched subtype and correspondingly more frequent ERBB2 amplification and higher HER2 protein abundance in the Chinese HR+HER2+ cohort, stressing anti-HER2 therapy for these individuals. Furthermore, comprehensive metabolomic and proteomic analyses revealed ferroptosis as a potential therapeutic target for basal-like tumors. The integration of clinical, transcriptomic, metabolomic, radiomic and pathological features allowed for efficient stratification of patients into groups with varying recurrence risks. Our study provides a public resource and new insights into the biology and ancestry specificity of breast cancer in the Asian population, offering potential for further precision treatment approaches.

Feasibility of functional precision medicine for guiding treatment of relapsed or refractory pediatric cancers

Children with rare, relapsed or refractory cancers often face limited treatment options, and few predictive biomarkers are available that can enable personalized treatment recommendations. The implementation of functional precision medicine (FPM), which combines genomic profiling with drug sensitivity testing (DST) of patient-derived tumor cells, has potential to identify treatment options when standard-of-care is exhausted. The goal of this prospective observational study was to generate FPM data for pediatric patients with relapsed or refractory cancer. The primary objective was to determine the feasibility of returning FPM-based treatment recommendations in real time to the FPM tumor board (FPMTB) within a clinically actionable timeframe (<4 weeks). The secondary objective was to assess clinical outcomes from patients enrolled in the study. Twenty-five patients with relapsed or refractory solid and hematological cancers were enrolled; 21 patients underwent DST and 20 also completed genomic profiling. Median turnaround times for DST and genomics were within 10 days and 27 days, respectively. Treatment recommendations were made for 19 patients (76%), of whom 14 received therapeutic interventions. Six patients received subsequent FPM-guided treatments. Among these patients, five (83%) experienced a greater than 1.3-fold improvement in progression-free survival associated with their FPM-guided therapy relative to their previous therapy, and demonstrated a significant increase in progression-free survival and objective response rate compared to those of eight non-guided patients. The findings from our proof-of-principle study illustrate the potential for FPM to positively impact clinical care for pediatric and adolescent patients with relapsed or refractory cancers and warrant further validation in large prospective studies. ClinicalTrials.gov registration: NCT03860376 .

Targeted next-generation sequencing of 491 lung cancers in clinical practice: Implications for future detection strategy and targeted therapy

Although lung cancer remains the most common cause of global cancer-related mortality, the identification of oncogenic driver alterations and the development of targeted drugs has dramatically altered the therapeutic landscape. In this retrospective study, we found that 97.7% samples carried at least one mutation in the 25 genes tested in our cohort. 53.6% samples were positive for EGFR mutations, followed by TP53 (41.1%), KRAS (11.8%), ERBB2 (4.3%). EGFR mutations were mainly found in female adenocarcinomas, while TP53 was mainly found in male non-adenocarcinomas. Significant differences can be found in the mutation rate of EGFR (60.9% vs 11.9%), KRAS (12.2% vs 25.0%), STK11 (1.5% vs 11.9%), FGFR3 (2.4% vs 0.0%) and ERBB4 (1.2% vs 6.1%) between adenocarcinoma in our cohort and TCGA-LUAD data (all p < 0.001). What's more, we found that the mutation of EGFR increased significantly from adenocarcinomas in situ (AIS, 21.4%) to microinvasive adenocarcinomas (MIA, 52.4%) and invasive adenocarcinomas (IA, 61.1%), while the mutation of ERBB2 dropped markedly from AIS (21.4%) to MIA (9.5%) and IA (4.1%). At last, comparations between targeted NGS and ARMS-based single gene test in the detection of EGFR showed a 94.6% consistence. In conclusion, targeted NGS can provide a comprehensive mutational profile of lung cancer. Considering the high mutation rate of EGFR in NSCLC of Asian populations, a specialized detection strategy should be conducted.

Anchored-fusion enables targeted fusion search in bulk and single-cell RNA sequencing data

Here, we present Anchored-fusion, a highly sensitive fusion gene detection tool. It anchors a gene of interest, which often involves driver fusion events, and recovers non-unique matches of short-read sequences that are typically filtered out by conventional algorithms. In addition, Anchored-fusion contains a module based on a deep learning hierarchical structure that incorporates self-distillation learning (hierarchical view learning and distillation [HVLD]), which effectively filters out false positive chimeric fragments generated during sequencing while maintaining true fusion genes. Anchored-fusion enables highly sensitive detection of fusion genes, thus allowing for application in cases with low sequencing depths. We benchmark Anchored-fusion under various conditions and found it outperformed other tools in detecting fusion events in simulated data, bulk RNA sequencing (bRNA-seq) data, and single-cell RNA sequencing (scRNA-seq) data. Our results demonstrate that Anchored-fusion can be a useful tool for fusion detection tasks in clinically relevant RNA-seq data and can be applied to investigate intratumor heterogeneity in scRNA-seq data.

Repurposing CRISPR-Cas13 systems for robust mRNA trans-splicing

Type VI CRISPR enzymes have been developed as programmable RNA-guided Cas proteins for eukaryotic RNA editing. Notably, Cas13 has been utilized for site-targeted single base edits, demethylation, RNA cleavage or knockdown and alternative splicing. However, the ability to edit large stretches of mRNA transcripts remains a significant challenge. Here, we demonstrate that CRISPR-Cas13 systems can be repurposed to assist trans-splicing of exogenous RNA fragments into an endogenous pre-mRNA transcript, a method termed CRISPR Assisted mRNA Fragment Trans-splicing (CRAFT). Using split reporter-based assays, we evaluate orthogonal Cas13 systems, optimize guide RNA length and screen for optimal trans-splicing site(s) across a range of intronic targets. We achieve markedly improved editing of large 5' and 3' segments in different endogenous mRNAs across various mammalian cell types compared to other spliceosome-mediated trans-splicing methods. CRAFT can serve as a versatile platform for attachment of protein tags, studying the impact of multiple mutations/single nucleotide polymorphisms, modification of untranslated regions (UTRs) or replacing large segments of mRNA transcripts.

Integrating Optical Genome Mapping and Whole Genome Sequencing in Somatic Structural Variant Detection

Structural variants drive tumorigenesis by disrupting normal gene function through insertions, inversions, translocations, and copy number changes, including deletions and duplications. Detecting structural variants is crucial for revealing their roles in tumor development, clinical outcomes, and personalized therapy. Presently, most studies rely on short-read data from next-generation sequencing that aligns back to a reference genome to determine if and, if so, where a structural variant occurs. However, structural variant discovery by short-read sequencing is challenging, primarily because of the difficulty in mapping regions of repetitive sequences. Optical genome mapping (OGM) is a recent technology used for imaging and assembling long DNA strands to detect structural variations. To capture the structural variant landscape more thoroughly in the human genome, we developed an integrated pipeline that combines Bionano OGM and Illumina whole-genome sequencing and applied it to samples from 29 pediatric B-ALL patients. The addition of OGM allowed us to identify 511 deletions, 506 insertions, 93 duplications/gains, and 145 translocations that were otherwise missed in the short-read data. Moreover, we identified several novel gene fusions, the expression of which was confirmed by RNA sequencing. Our results highlight the benefit of integrating OGM and short-read detection methods to obtain a comprehensive analysis of genetic variation that can aid in clinical diagnosis, provide new therapeutic targets, and improve personalized medicine in cancers driven by structural variation.

Deep whole-genome analysis of 494 hepatocellular carcinomas

Over half of hepatocellular carcinoma (HCC) cases diagnosed worldwide are in China1-3. However, whole-genome analysis of hepatitis B virus (HBV)-associated HCC in Chinese individuals is limited4-8, with current analyses of HCC mainly from non-HBV-enriched populations9,10. Here we initiated the Chinese Liver Cancer Atlas (CLCA) project and performed deep whole-genome sequencing (average depth, 120×) of 494 HCC tumours. We identified 6 coding and 28 non-coding previously undescribed driver candidates. Five previously undescribed mutational signatures were found, including aristolochic-acid-associated indel and doublet base signatures, and a single-base-substitution signature that we termed SBS_H8. Pentanucleotide context analysis and experimental validation confirmed that SBS_H8 was distinct to the aristolochic-acid-associated SBS22. Notably, HBV integrations could take the form of extrachromosomal circular DNA, resulting in elevated copy numbers and gene expression. Our high-depth data also enabled us to characterize subclonal clustered alterations, including chromothripsis, chromoplexy and kataegis, suggesting that these catastrophic events could also occur in late stages of hepatocarcinogenesis. Pathway analysis of all classes of alterations further linked non-coding mutations to dysregulation of liver metabolism. Finally, we performed in vitro and in vivo assays to show that fibrinogen alpha chain (FGA), determined as both a candidate coding and non-coding driver, regulates HCC progression and metastasis. Our CLCA study depicts a detailed genomic landscape and evolutionary history of HCC in Chinese individuals, providing important clinical implications.

Mosaic BRAF Fusions Are a Recurrent Cause of Congenital Melanocytic Nevi Targetable by MAPK Pathway Inhibition

Among children with multiple congenital melanocytic nevi, 25% have no established genetic cause, of whom many develop a hyperproliferative and severely pruritic phenotype resistant to treatment. Gene fusions have been reported in individual cases of congenital melanocytic nevi. We studied 169 patients with congenital melanocytic nevi in this study, 38 of whom were double wild type for pathogenic NRAS/BRAF variants. Nineteen of these 38 patients had sufficient tissue to undergo RNA sequencing, which revealed mosaic BRAF fusions in 11 of 19 patients and mosaic RAF1 fusions in 1 of 19. Recurrently, fusions involved the loss of the 5´ regulatory domain of BRAF or RAF1 but preserved the kinase domain. We validated all cases and detected the fusions in two separate nevi in 5 of 12 patients, confirming clonality. The absence of the fusion in blood in 8 of 12 patients indicated mosaicism. Primary culture of BRAF-fusion nevus cells from 3 of 12 patients demonstrated highly increased MAPK activation, despite only mildly increased BRAF expression, suggesting additional mechanisms of kinase activation. Trametinib quenched MAPK hyperactivation in vitro, and treatment of two patients caused rapid improvement in bulk tissue, improving bodily movement and reducing inflammation and severe pruritus. These findings offer a genetic diagnosis to an additional group of patients and trametinib as a treatment option for the severe associated phenotypes.

Recent advances in the investigation of fusion RNAs and their role in molecular pathology of cancer

Gene fusions have had a significant role in the development of various types of cancer, oftentimes involved in oncogenic activities through dysregulation of gene expression or signalling pathways. Some cancer-associated chromosomal translocations can undergo backsplicing, resulting in fusion-circular RNAs, a more stable isoform immune to RNase degradation. This stability makes fusion circular RNAs a promising diagnostic biomarker for cancer. While the detection of linear fusion RNAs and their function in certain cancers have been described in literature, fusion circular RNAs lag behind due to their low abundance in cancer cells. This review highlights current literature on the role of linear and circular fusion transcripts in cancer, tools currently available for detecting of these chimeric RNAs and their function and how they play a role in tumorigenesis.

Spindle cell neoplasms with novel LTK fusion - Expanding the spectrum of kinase fusion-positive soft tissue tumors

AIMS

Kinase fusion-positive soft tissue tumors represent an emerging, molecularly defined group of mesenchymal tumors with a wide morphologic spectrum and diverse activating kinases. Here, we present two cases of soft tissue tumors with novel LTK fusions.

METHODS AND RESULTS

Both cases presented as acral skin nodules (big toe and middle finger) in pediatric patients (17-year-old girl and 2-year-old boy). The tumors measured 2 and 3 cm in greatest dimension. Histologically, both cases exhibited bland-looking spindle cells infiltrating adipose tissue and accompanied by collagenous stroma. One case additionally displayed perivascular hyalinization and band-like stromal collagen. Both cases exhibited focal S100 staining, and one case had patchy coexpression of CD34. Targeted RNA-seq revealed the presence of novel in-frame MYH9::LTK and MYH10::LTK fusions, resulting in upregulation of LTK expression. Of interest, DNA methylation-based unsupervised clustering analysis in one case showed that the tumor clustered with dermatofibrosarcoma protuberans (DFSP). One tumor was excised with amputation with no local recurrence or distant metastasis at 18-month follow-up. The other case was initially marginally excised with local recurrence after one year, followed by wide local excision, with no evidence of disease at 10 years of follow-up.

CONCLUSIONS

This is the first reported case series of soft tissue tumors harboring LTK fusion, expanding the molecular landscape of soft tissue tumors driven by activating kinase fusions. Furthermore, studies involving a larger number of cases and integrated genomic analyses will be warranted to fully elucidate the pathogenesis and classification of these tumors.

CTAT-LR-fusion: accurate fusion transcript identification from long and short read isoform sequencing at bulk or single cell resolution

Gene fusions are found as cancer drivers in diverse adult and pediatric cancers. Accurate detection of fusion transcripts is essential in cancer clinical diagnostics, prognostics, and for guiding therapeutic development. Most currently available methods for fusion transcript detection are compatible with Illumina RNA-seq involving highly accurate short read sequences. Recent advances in long read isoform sequencing enable the detection of fusion transcripts at unprecedented resolution in bulk and single cell samples. Here we developed a new computational tool CTAT-LR-fusion to detect fusion transcripts from long read RNA-seq with or without companion short reads, with applications to bulk or single cell transcriptomes. We demonstrate that CTAT-LR-fusion exceeds fusion detection accuracy of alternative methods as benchmarked with simulated and real long read RNA-seq. Using short and long read RNA-seq, we further apply CTAT-LR-fusion to bulk transcriptomes of nine tumor cell lines, and to tumor single cells derived from a melanoma sample and three metastatic high grade serous ovarian carcinoma samples. In both bulk and in single cell RNA-seq, long isoform reads yielded higher sensitivity for fusion detection than short reads with notable exceptions. By combining short and long reads in CTAT-LR-fusion, we are able to further maximize detection of fusion splicing isoforms and fusion-expressing tumor cells. CTAT-LR-fusion is available at https://github.com/TrinityCTAT/CTAT-LR-fusion/wiki.

Neurotrophic-tyrosine receptor kinase gene fusion in papillary thyroid cancer: A clinicogenomic biobank and record linkage study from Finland

Selective tropomyosin receptor kinase (TRK) inhibitors are approved targeted therapies for patients with solid tumors harboring a neurotrophic tyrosine receptor kinase (NTRK) gene fusion. Country-specific estimates of NTRK gene fusion frequency, and knowledge on the characteristics of affected patients, are limited. We identified patients with histologically-confirmed papillary thyroid cancer (PTC) from Finland's Auria Biobank. TRK protein expression was determined by pan-TRK immunohistochemistry. Immuno-stained tumor samples were scored by a certified pathologist. Gene fusions and other co-occurring gene alterations were identified by next generation sequencing. Patient characteristics and vital status were determined from linked hospital electronic health records (EHRs). Patients were followed from 1 year before PTC diagnosis until death. 6/389 (1.5%) PTC patients had an NTRK gene fusion (all NTRK3); mean age 43.8 years (and none had comorbidities) at PTC diagnosis. Gene fusion partners were EML4 (n = 3), ETV6 (n = 2), and RBPMS (n = 1). Of 3/6 patients with complete EHRs, all received radioactive iodine ablation only and were alive at end of follow-up (median observation, 9.12 years). In conclusion, NTRK gene fusion is infrequent in patients with PTC. Linkage of biobank samples to EHRs is feasible in describing the characteristics and outcomes of patients with PTC and potentially other cancer types.

Disruption of the TP53 locus in osteosarcoma leads to TP53 promoter gene fusions and restoration of parts of the TP53 signalling pathway

TP53 is the most frequently mutated gene in human cancer. This gene shows not only loss-of-function mutations but also recurrent missense mutations with gain-of-function activity. We have studied the primary bone malignancy osteosarcoma, which harbours one of the most rearranged genomes of all cancers. This is odd since it primarily affects children and adolescents who have not lived the long life thought necessary to accumulate massive numbers of mutations. In osteosarcoma, TP53 is often disrupted by structural variants. Here, we show through combined whole-genome and transcriptome analyses of 148 osteosarcomas that TP53 structural variants commonly result in loss of coding parts of the gene while simultaneously preserving and relocating the promoter region. The transferred TP53 promoter region is fused to genes previously implicated in cancer development. Paradoxically, these erroneously upregulated genes are significantly associated with the TP53 signalling pathway itself. This suggests that while the classical tumour suppressor activities of TP53 are lost, certain parts of the TP53 signalling pathway that are necessary for cancer cell survival and proliferation are retained. In line with this, our data suggest that transposition of the TP53 promoter is an early event that allows for a new normal state of genome-wide rearrangements in osteosarcoma. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Whole genome sequencing in clinical practice

Whole genome sequencing (WGS) is becoming the preferred method for molecular genetic diagnosis of rare and unknown diseases and for identification of actionable cancer drivers. Compared to other molecular genetic methods, WGS captures most genomic variation and eliminates the need for sequential genetic testing. Whereas, the laboratory requirements are similar to conventional molecular genetics, the amount of data is large and WGS requires a comprehensive computational and storage infrastructure in order to facilitate data processing within a clinically relevant timeframe. The output of a single WGS analyses is roughly 5 MIO variants and data interpretation involves specialized staff collaborating with the clinical specialists in order to provide standard of care reports. Although the field is continuously refining the standards for variant classification, there are still unresolved issues associated with the clinical application. The review provides an overview of WGS in clinical practice - describing the technology and current applications as well as challenges connected with data processing, interpretation and clinical reporting.

Linking Gene Fusions to Bone Marrow Failure and Malignant Transformation in Dyskeratosis Congenita

Dyskeratosis Congenita (DC) is a multisystem disorder intrinsically associated with telomere dysfunction, leading to bone marrow failure (BMF). Although the pathology of DC is largely driven by mutations in telomere-associated genes, the implications of gene fusions, which emerge due to telomere-induced genomic instability, remain unexplored. We meticulously analyzed gene fusions in RNA-Seq data from DC patients to provide deeper insights into DC's progression. The most significant DC-specific gene fusions were subsequently put through in silico assessments to ascertain biophysical and structural attributes, including charge patterning, inherent disorder, and propensity for self-association. Selected candidates were then analyzed using deep learning-powered structural predictions and molecular dynamics simulations to gauge their potential for forming higher-order oligomers. Our exploration revealed that genes participating in fusion events play crucial roles in upholding genomic stability, facilitating hematopoiesis, and suppressing tumors. Notably, our analysis spotlighted a particularly disordered polyampholyte fusion protein that exhibits robust higher-order oligomerization dynamics. To conclude, this research underscores the potential significance of several high-confidence gene fusions in the progression of BMF in DC, particularly through the dysregulation of genomic stability, hematopoiesis, and tumor suppression. Additionally, we propose that these fusion proteins might hold a detrimental role, specifically in inducing proteotoxicity-driven hematopoietic disruptions.

Rare Case of a Turner Syndrome Patient with Metastatic Dysgerminoma and No Y-Chromosomal Material with Pathogenic Variants Found in KIT and MTOR

INTRODUCTION

The presence of Y-chromosomal material in females with Turner syndrome (TS) is a well-established risk factor for developing gonadoblastoma and malignant transformations thereof. However, these events are rarely seen in TS patients with no Y-chromosomal material. Thus, it is the current understanding that parts of the Y-chromosome are essential for the malignant transformation of gonadoblastoma in the dysgenetic gonad.

METHODS

We report a case of a TS female with an apparent 46,X,idic(Xq) karyotype, who was diagnosed with a metastatic dysgerminoma. Whole exome sequencing of the tumor and blood, along with RNA sequencing of the tumor, was performed to comprehensively search for cryptic Y-chromosomal material and pathogenic variants.

RESULTS

No Y-chromosomal material was detected in either tumor or blood. Whole exome-sequencing of DNA and RNA revealed a pathogenic somatic gain-of-function mutation in KIT and a pathogenic missense mutation in MTOR. The patient underwent total hysterectomy with bilateral salpingo-oophorectomy, followed by adjuvant chemotherapy. Unfortunately, she died due to chemotherapy-induced pneumonitis 7 months after the initial diagnosis.

CONCLUSION

Females with TS can develop metastatic dysgerminoma even in the absence of Y-chromosomal material. This questions the current understanding of Y-chromosomal material being essential for the malignant transformation of a gonadoblastoma in the dysgenetic gonad.

FusionNeoAntigen: a resource of fusion gene-specific neoantigens

Among the diverse sources of neoantigens (i.e. single-nucleotide variants (SNVs), insertions or deletions (Indels) and fusion genes), fusion gene-derived neoantigens are generally more immunogenic, have multiple targets per mutation and are more widely distributed across various cancer types. Therefore, fusion gene-derived neoantigens are a potential source of highly immunogenic neoantigens and hold great promise for cancer immunotherapy. However, the lack of fusion protein sequence resources and knowledge prevents this application. We introduce 'FusionNeoAntigen', a dedicated resource for fusion-specific neoantigens, accessible at https://compbio.uth.edu/FusionNeoAntigen. In this resource, we provide fusion gene breakpoint crossing neoantigens focused on ∼43K fusion proteins of ∼16K in-frame fusion genes from FusionGDB2.0. FusionNeoAntigen provides fusion gene information, corresponding fusion protein sequences, fusion breakpoint peptide sequences, fusion gene-derived neoantigen prediction, virtual screening between fusion breakpoint peptides having potential fusion neoantigens and human leucocyte antigens (HLAs), fusion breakpoint RNA/protein sequences for developing vaccines, information on samples with fusion-specific neoantigen, potential CAR-T targetable cell-surface fusion proteins and literature curation. FusionNeoAntigen will help to develop fusion gene-based immunotherapies. We will report all potential fusion-specific neoantigens from all possible open reading frames of ∼120K human fusion genes in future versions.

Gene Rearrangement and Expression of PRKACA and PRKACB Govern Morphobiology of Pancreatobiliary Oncocytic Neoplasms

Intraductal oncocytic papillary neoplasms (IOPNs) are distinct from intraductal papillary mucinous neoplasms based on characteristic morphologic and genetic features represented by fusion genes involving PRKACA or PRKACB (PRKACA/B). However, pancreatic and biliary tumors with partial oncocytic features are often encountered clinically, and their molecular features are yet to be clarified. This study included 80 intraductal papillary neoplasms: 32 tumors with mature IOPN morphology (typical), 28 with partial or subclonal oncocytic features (atypical), and 20 without oncocytic features (control). We analyzed PRKACA/B fusion genes, including ATP1B1::PRKACA, DNAJB1::PRKACA, and ATP1B1::PRKACB, by reverse-transcription PCR; mRNA expression of fusion genes and nonrearranged PRKACA/B genes by quantitative reverse-transcription PCR; mutations in KRAS, BRAF, and GNAS by targeted sequencing or droplet digital PCR; and the expression of cyclic adenosine monophosphate (cAMP)-dependent protein kinase catalytic subunits α (PRKACA) and β (PRKACB), phosphorylated cAMP response element-binding protein, and aberrations of p16, p53, SMAD4, STK11, and β-catenin by immunohistochemistry. PRKACA/B fusion genes were detected in 100% (32/32) of typical, 46% (13/28) of atypical, and 0% (0/20) of control (P < .05). Expression of PRKACA, PRKACB, and phosphorylated cAMP response element-binding protein was upregulated in neoplasms with PRKACA/B fusion genes (P < .05). mRNA expression of the PRKACA/B fusion genes and protein expression of PRKACA or PRKACB tended to be higher in typical than in atypical cases (mRNA, P = .002; protein expression, P = .054). In some atypical neoplasms with mixed subtypes, PRKACA/B fusion genes were superimposed exclusively on oncocytic components. Typical IOPNs harbored fewer KRAS and GNAS mutations than control samples and fewer alterations in p53 and STK11 than atypical samples (P < .05). In conclusion, PRKACA/B fusion genes not only are the characteristic drivers of IOPNs but also play a crucial role in the development of subclonal oncocytic neoplasms. Moreover, oncocytic morphology is strongly associated with upregulation of PRKACA/B, which may provide clues for potential therapeutic options.

RNASeq Analysis for Accurate Identification of Fusion Partners in Tumor Specific Translocations Detected by Standard FISH Probes in Hematologic Malignancies

Background

Fluorescence labeled DNA probes and in situ hybridization methods had shorter turn round time for results revolutionized their clinical application. Signals obtained from these probes are highly specific, yet they can produce fusion signals not necessarily representing fusion of actual genes due to other genes included in the probe design. In this study we evaluated discordance between cytogenetic, FISH and RNAseq results in 3 different patients with hematologic malignancies and illustrated the need to perform next generation sequencing (NGS) or RNASeq to accurately interpret FISH results.

Methods

Bone marrow or peripheral blood karyotypes and FISH were performed to detect recurring translocations associated with hematologic malignancies in clinical samples routinely referred to our clinical cytogenetics laboratory. When required, NGS was performed on DNA and RNA libraries to detect somatic alterations and gene fusions in some of these specimens. Discordance in results between these methods is further evaluated.

Results

For a patient with plasma cell leukemia standard FGFR3 / IGH dual fusion FISH assay detected fusion that was interpreted as FGFR3-positive leukemia, whereas NGS/RNASeq detected NSD2::IGH. For a pediatric acute lymphoblastic leukemia patient, a genetic diagnosis of PDGFRB-positive ALL was rendered because the PDGFRB break-apart probe detected clonal rearrangement, whereas NGS detected MEF2D::CSF1R. A MYC-positive B-prolymphocytic leukemia was rendered for another patient with a cytogenetically identified t(8;14) and MYC::IGH by FISH, whereas NGS detected a novel PVT1::RCOR1 not previously reported.

Conclusions

These are 3 cases in a series of several other concordant results, nevertheless, elucidate limitations when interpreting FISH results in clinical applications, particularly when other genes are included in probe design. In addition, when the observed FISH signals are atypical, this study illustrates the necessity to perform complementary laboratory assays, such as NGS and/or RNASeq, to accurately identify fusion genes in tumorigenic translocations.

Molecular heterogeneity in prostate cancer and the role of targeted therapy

Data collected from large-scale studies has shown that the incidence of prostate cancer globally is on the rise, which could be attributed to an overall increase in lifespan. So, the question is how has modern science with all its new technologies and clinical breakthroughs mitigated or managed this disease? The answer is not a simple one as prostate cancer exhibits various subtypes, each with its unique characteristics or signatures which creates challenges in treatment. To understand the complexity of prostate cancer these signatures must be deciphered. Molecular studies of prostate cancer samples have identified certain genetic and epigenetic alterations, which are instrumental in tumorigenesis. Some of these candidates include the androgen receptor (AR), various oncogenes, tumor suppressor genes, and the tumor microenvironment, which serve as major drivers that lead to cancer progression. These aberrant genes and their products can give an insight into prostate cancer development and progression by acting as potent markers to guide future therapeutic approaches. Thus, understanding the complexity of prostate cancer is crucial for targeting specific markers and tailoring treatments accordingly.

Osteosarcomas With Few Chromosomal Alterations or Adult Onset Are Genetically Heterogeneous

Osteosarcoma is the most common primary bone malignancy, often detected in children and adolescents and commonly associated with TP53 alterations along with a high number of chromosomal rearrangements. However, osteosarcoma can affect patients of any age, and some tumors display less genetic complexity. Besides TP53 variants, data on key driving mutations are lacking for many osteosarcomas, particularly those affecting adults. To detect osteosarcoma-specific alterations, we screened transcriptomic and genomic sequencing and copy number data from 150 bone tumors originally diagnosed as osteosarcomas. To increase the precision in gene fusion detection, we developed a bioinformatic tool denoted as NAFuse, which extracts gene fusions that are verified at both the genomic and transcriptomic levels. Apart from the already reported genetic subgroups of osteosarcoma with TP53 structural variants, or MDM2 and/or CDK4 amplification, we did not identify any recurrent genetic driver that signifies the remaining cases. Among the plethora of mutations identified, we found genetic alterations characteristic of, or similar to, those of other bone and soft tissue tumors in 8 cases. These mutations were found in tumors with relatively few other genetic alterations or in adults. Due to the lack of clinical context and available tissue, we can question the diagnosis only on a genetic basis. However, our findings support the notion that osteosarcomas with few chromosomal alterations or adult onset seem genetically distinct from conventional osteosarcomas of children and adolescents.

Improved detection of clinically relevant fusion transcripts in cancer by machine learning classification

BACKGROUND

Genomic rearrangements in cancer cells can create fusion genes that encode chimeric proteins or alter the expression of coding and non-coding RNAs. In some cancer types, fusions involving specific kinases are used as targets for therapy. Fusion genes can be detected by whole genome sequencing (WGS) and targeted fusion panels, but RNA sequencing (RNA-Seq) has the advantageous capability of broadly detecting expressed fusion transcripts.

RESULTS

We developed a pipeline for validation of fusion transcripts identified in RNA-Seq data using matched WGS data from The Cancer Genome Atlas (TCGA) and applied it to 910 tumors from 11 different cancer types. This resulted in 4237 validated gene fusions, 3049 of them with at least one identified genomic breakpoint. Utilizing validated fusions as true positive events, we trained a machine learning classifier to predict true and false positive fusion transcripts from RNA-Seq data. The final precision and recall metrics of the classifier were 0.74 and 0.71, respectively, in an independent dataset of 249 breast tumors. Application of this classifier to all samples with RNA-Seq data from these cancer types vastly extended the number of likely true positive fusion transcripts and identified many potentially targetable kinase fusions. Further analysis of the validated gene fusions suggested that many are created by intrachromosomal amplification events with microhomology-mediated non-homologous end-joining.

CONCLUSIONS

A classifier trained on validated fusion events increased the accuracy of fusion transcript identification in samples without WGS data. This allowed the analysis to be extended to all samples with RNA-Seq data, facilitating studies of tumor biology and increasing the number of detected kinase fusions. Machine learning could thus be used in identification of clinically relevant fusion events for targeted therapy. The large dataset of validated gene fusions generated here presents a useful resource for development and evaluation of fusion transcript detection algorithms.

Germline gene fusions across species reveal the chromosomal instability regions and cancer susceptibility

The canine transmissible venereal tumor (CTVT) is a clonal cell-mediated cancer with a long evolutionary history and extensive karyotype rearrangements in its genome. However, little is known about its genetic similarity to human tumors. Here, using multi-omics data we identified 11 germline gene fusions (GGFs) in CTVT, which showed higher genetic susceptibility than others. Additionally, we illustrate a mechanism of a complex gene fusion of three gene segments (HSD17B4-DMXL1-TNFAIP8) that we refer to "greedy fusion". Our findings also provided evidence that expressions of GGFs are downregulated during the tumor regressive phase, which is associated with DNA methylation level. This study presents a comprehensive landscape of gene fusions (GFs) in CTVT, which offers a valuable genetic resource for exploring potential genetic mechanisms underlying the development of cancers in both dogs and humans.

A comprehensive assessment of exome capture methods for RNA sequencing of formalin-fixed and paraffin-embedded samples

RNA-Seq analysis of Formalin-Fixed and Paraffin-Embedded (FFPE) samples has emerged as a highly effective approach and is increasingly being used in clinical research and drug development. However, the processing and storage of FFPE samples are known to cause extensive degradation of RNAs, which limits the discovery of gene expression or gene fusion-based biomarkers using RNA sequencing, particularly methods reliant on Poly(A) enrichment. Recently, researchers have developed an exome targeted RNA-Seq methodology that utilizes biotinylated oligonucleotide probes to enrich RNA transcripts of interest, which could overcome these limitations. Nevertheless, the standardization of this experimental framework, including probe designs, sample multiplexing, sequencing read length, and bioinformatic pipelines, remains an essential requirement. In this study, we conducted a comprehensive comparison of three main commercially available exome capture kits and evaluated key experimental parameters, to provide the overview of the advantages and limitations associated with the selection of library preparation protocols and sequencing platforms. The results provide valuable insights into the best practices for obtaining high-quality data from FFPE samples.

Utilizing immunogenomic approaches to prioritize targetable neoantigens for personalized cancer immunotherapy

Advancements in sequencing technologies and bioinformatics algorithms have expanded our ability to identify tumor-specific somatic mutation-derived antigens (neoantigens). While recent studies have shown neoantigens to be compelling targets for cancer immunotherapy due to their foreign nature and high immunogenicity, the need for increasingly accurate and cost-effective approaches to rapidly identify neoantigens remains a challenging task, but essential for successful cancer immunotherapy. Currently, gene expression analysis and algorithms for variant calling can be used to generate lists of mutational profiles across patients, but more care is needed to curate these lists and prioritize the candidate neoantigens most capable of inducing an immune response. A growing amount of evidence suggests that only a handful of somatic mutations predicted by mutational profiling approaches act as immunogenic neoantigens. Hence, unbiased screening of all candidate neoantigens predicted by Whole Genome Sequencing/Whole Exome Sequencing may be necessary to more comprehensively access the full spectrum of immunogenic neoepitopes. Once putative cancer neoantigens are identified, one of the largest bottlenecks in translating these neoantigens into actionable targets for cell-based therapies is identifying the cognate T cell receptors (TCRs) capable of recognizing these neoantigens. While many TCR-directed screening and validation assays have utilized bulk samples in the past, there has been a recent surge in the number of single-cell assays that provide a more granular understanding of the factors governing TCR-pMHC interactions. The goal of this review is to provide an overview of existing strategies to identify candidate neoantigens using genomics-based approaches and methods for assessing neoantigen immunogenicity. Additionally, applications, prospects, and limitations of some of the current single-cell technologies will be discussed. Finally, we will briefly summarize some of the recent models that have been used to predict TCR antigen specificity and analyze the TCR receptor repertoire.

A novel chimeric RNA originating from BmCPV S4 and Bombyx mori HDAC11 transcripts regulates virus proliferation

Polymerases encoded by segmented negative-strand RNA viruses cleave 5'-m7G-capped host transcripts to prime viral mRNA synthesis ("cap-snatching") to generate chimeric RNA, and trans-splicing occurs between viral and cellular transcripts. Bombyx mori cytoplasmic polyhedrosis virus (BmCPV), an RNA virus belonging to Reoviridae, is a major pathogen of silkworm (B. mori). The genome of BmCPV consists of 10 segmented double-stranded RNAs (S1-S10) from which viral RNAs encoding a protein are transcribed. In this study, chimeric silkworm-BmCPV RNAs, in which the sequence derived from the silkworm transcript could fuse with both the 5' end and the 3' end of viral RNA, were identified in the midgut of BmCPV-infected silkworms by RNA_seq and further confirmed by RT-PCR and Sanger sequencing. A novel chimeric RNA, HDAC11-S4 RNA 4, derived from silkworm histone deacetylase 11 (HDAC11) and the BmCPV S4 transcript encoding viral structural protein 4 (VP4), was selected for validation by in situ hybridization and Northern blotting. Interestingly, our results indicated that HDAC11-S4 RNA 4 was generated in a BmCPV RNA-dependent RNA polymerase (RdRp)-independent manner and could be translated into a truncated BmCPV VP4 with a silkworm HDAC11-derived N-terminal extension. Moreover, it was confirmed that HDAC11-S4 RNA 4 inhibited BmCPV proliferation, decreased the level of H3K9me3 and increased the level of H3K9ac. These results indicated that during infection with BmCPV, a novel mechanism, different from that described in previous reports, allows the genesis of chimeric silkworm-BmCPV RNAs with biological functions.

Somatic mutation burden in relation to aging and functional life span: implications for cellular reprogramming and rejuvenation

The accrual of somatic mutations has been implicated as causal factors in aging since the 1950s. However, the quantitative analysis of somatic mutations has posed a major challenge due to the random nature of de novo mutations in normal tissues, which has limited analysis to tumors and other clonal lineages. Advances in single-cell and single-molecule next-generation sequencing now allow to obtain, for the first time, detailed insights into the landscape of somatic mutations in different human tissues and cell types as a function of age under various conditions. Here, we will briefly recapitulate progress in somatic mutation analysis and discuss the possible relationship between somatic mutation burden with functional life span, with a focus on differences between germ cells, stem cells, and differentiated cells.

The subgenome Saccharum spontaneum contributes to sugar accumulation in sugarcane as revealed by full-length transcriptomic analysis

INTRODUCTION

Modern sugarcane cultivars (Saccharum spp. hybrids) derived from crosses between S. officinarum and S. spontaneum, with high-sugar traits and excellent stress tolerance inherited respectively. However, the contribution of the S. spontaneum subgenome to sucrose accumulation is still unclear.

OBJECTIVE

To compensate for the absence of a high-quality reference genome, a transcriptome analysis method is needed to analyze the molecular basis of differential sucrose accumulation in sugarcane hybrids and to find clues to the contribution of the S. spontaneum subgenome to sucrose accumulation.

METHODS

PacBio full-length sequencing was used to complement genome annotation, followed by the identification of differential genes between the high and low sugar groups using differential alternative splicing analysis and differential expression analysis. At the subgenomic level, the factors responsible for differential sucrose accumulation were investigated from the perspective of transcriptional and post-transcriptional regulation.

RESULTS

A full-length transcriptome annotated at the subgenomic level was provided, complemented by 263,378 allele-defined transcript isoforms and 139,405 alternative splicing (AS) events. Differential alternative splicing (DA) analysis and differential expression (DE) analysis identified differential genes between high and low sugar groups and explained differential sucrose accumulation factors by the KEGG pathways. In some gene models, different or even opposite expression patterns of alleles from the same gene were observed, reflecting the potential evolution of these alleles toward novel functions in polyploid sugarcane. Among DA and DE genes in the sucrose source-sink complex pathway, we found some alleles encoding sucrose accumulation-related enzymes derived from the S. spontaneum subgenome were differentially expressed or had DA events between the two contrasting sugarcane hybrids.

CONCLUSION

Full-length transcriptomes annotated at the subgenomic level could better characterize sugarcane hybrids, and the S. spontaneum subgenome was found to contribute to sucrose accumulation.

Genomic profiling of subcutaneous patient-derived xenografts reveals immune constraints on tumor evolution in childhood solid cancer

Subcutaneous patient-derived xenografts (PDXs) are an important tool for childhood cancer research. Here, we describe a resource of 68 early passage PDXs established from 65 pediatric solid tumor patients. Through genomic profiling of paired PDXs and patient tumors (PTs), we observe low mutational similarity in about 30% of the PT/PDX pairs. Clonal analysis in these pairs show an aggressive PT minor subclone seeds the major clone in the PDX. We show evidence that this subclone is more immunogenic and is likely suppressed by immune responses in the PT. These results suggest interplay between intratumoral heterogeneity and antitumor immunity may underlie the genetic disparity between PTs and PDXs. We further show that PDXs generally recapitulate PTs in copy number and transcriptomic profiles. Finally, we report a gene fusion LRPAP1-PDGFRA. In summary, we report a childhood cancer PDX resource and our study highlights the role of immune constraints on tumor evolution.