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

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.

Identification of tumor rejection antigens and the immunologic landscape of medulloblastoma

BACKGROUND

The current standard of care treatments for medulloblastoma are insufficient as these do not take tumor heterogeneity into account. Newer, safer, patient-specific treatment approaches are required to treat high-risk medulloblastoma patients who are not cured by the standard therapies. Immunotherapy is a promising treatment modality that could be key to improving survival and avoiding morbidity. For an effective immune response, appropriate tumor antigens must be targeted. While medulloblastoma patients with subgroup-specific genetic substitutions have been previously reported, the immunogenicity of these genetic alterations remains unknown. The aim of this study is to identify potential tumor rejection antigens for the development of antigen-directed cellular therapies for medulloblastoma.

METHODS

We developed a cancer immunogenomics pipeline and performed a comprehensive analysis of medulloblastoma subgroup-specific transcription profiles (n = 170, 18 WNT, 46 SHH, 41 Group 3, and 65 Group 4 patient tumors) available through International Cancer Genome Consortium (ICGC) and European Genome-Phenome Archive (EGA). We performed in silico antigen prediction across a broad array of antigen classes including neoantigens, tumor-associated antigens (TAAs), and fusion proteins. Furthermore, we evaluated the antigen processing and presentation pathway in tumor cells and the immune infiltrating cell landscape using the latest computational deconvolution methods.

RESULTS

Medulloblastoma patients were found to express multiple private and shared immunogenic antigens. The proportion of predicted TAAs was higher than neoantigens and gene fusions for all molecular subgroups, except for sonic hedgehog (SHH), which had a higher neoantigen burden. Importantly, cancer-testis antigens, as well as previously unappreciated neurodevelopmental antigens, were found to be expressed by most patients across all medulloblastoma subgroups. Despite being immunologically cold, medulloblastoma subgroups were found to have distinct immune cell gene signatures.

CONCLUSIONS

Using a custom antigen prediction pipeline, we identified potential tumor rejection antigens with important implications for the development of immunotherapy for medulloblastoma.

Medulloblastomas Initiated by Homologous Recombination Defects in Mice

Germline mutations of homologous-recombination (HR) genes are among the top contributors to medulloblastomas. A significant portion of human medulloblastomas exhibited genomic signatures of HR defects. We asked whether ablation of Brca2, Palb2, and their related Brca1 and Bccip genes, in the mouse brain can differentially initiate medulloblastomas. We established conditional knockout mouse models of these HR-genes, and a conditional knockdown of Bccip (shBccip-KD). Deletion of any of these genes led to microcephaly and neurological defects, with Brca1- and Bccip- producing the worst. Trp53 co-deletion significantly rescued the microcephaly with Brca1, Palb2, and Brca2 deficiency, but it exhibited limited impact on Bccip- mice. For the first time, inactivation of either Brca1 or Palb2 with Trp53 was found to induce medulloblastomas. Bccip/Trp53 deletions failed, despite that shBccip-CKD was highly penetrative, to induce medulloblastomas. The tumors displayed diverse immunohistochemical features and chromosome copy number variation. While there were widespread upregulations of cell proliferative pathways, most of the tumors expressed biomarkers of the Sonic Hedgehog subgroup. The MBs developed from Brca1-, Palb2-, and Brca2- mice were highly sensitive to a PARP inhibitor, but not the ones from shBccip-CKD mice. Our models recapitulate the spontaneous medulloblastoma development with high penetrance and narrow time-window, providing ideal platforms to test therapeutic agents with the ability to differentiate HR defective and proficient tumors.

Uterine tumors mimicking ovarian sex cord tumors with rhabdoid differentiation: a clinicopathologic study of 4 cases: A case series analysis

RATIONALE

Uterine tumors resembling ovarian sex cord tumors (UTROSCT) with rhabdoid features are uncommon mesenchymal neoplasms exhibiting diverse histological patterns, including significant rhabdoid morphology. A thorough comprehension of their clinicopathologic features is crucial for precise diagnosis and effective management.

PATIENT CONCERNS

This study presents 4 cases of UTROSCT with rhabdoid features, diagnosed in patients aged 31 to 58. Varied recurrence patterns were observed, including similar recurrent lesions to the primary tumors with subsequent mortality, initial invasion and lymph node metastasis, and presence of only primary tumor.

DIAGNOSES

Histopathological examination revealed diverse morphological patterns, prominently featuring rhabdoid differentiation. Immunohistochemical analysis showed expression of hormone receptors, sex cord, smooth muscle, and epithelial markers, notably WT1, CD56, and CD99. Molecular analysis identified ESR1-NCOA2 fusions and ESR1 and NCOA2/3 rearrangements, indicating a potential association between these genetic alterations and extensive rhabdoid differentiation.

INTERVENTIONS

Various treatments were administered post-recurrence, including chemotherapy and targeted therapies. However, poor clinical outcomes were observed in all cases.

OUTCOMES

Despite aggressive treatments, including chemotherapy and targeted therapies, poor clinical outcomes were observed, highlighting the aggressive nature of UTROSCT with significant rhabdoid differentiation.

LESSONS

This case series emphasizes the importance of detailed pathological reporting, comprehensive molecular testing, and thorough tumor staging in UTROSCT cases with rhabdoid features. Enhanced understanding of the clinicopathologic characteristics of UTROSCT with rhabdoid differentiation is crucial for accurate diagnosis, prognostication, and management strategies.

Misexpression of inactive genes in whole blood is associated with nearby rare structural variants

Gene misexpression is the aberrant transcription of a gene in a context where it is usually inactive. Despite its known pathological consequences in specific rare diseases, we have a limited understanding of its wider prevalence and mechanisms in humans. To address this, we analyzed gene misexpression in 4,568 whole-blood bulk RNA sequencing samples from INTERVAL study blood donors. We found that while individual misexpression events occur rarely, in aggregate they were found in almost all samples and a third of inactive protein-coding genes. Using 2,821 paired whole-genome and RNA sequencing samples, we identified that misexpression events are enriched in cis for rare structural variants. We established putative mechanisms through which a subset of SVs lead to gene misexpression, including transcriptional readthrough, transcript fusions, and gene inversion. Overall, we develop misexpression as a type of transcriptomic outlier analysis and extend our understanding of the variety of mechanisms by which genetic variants can influence gene expression.

Prognostic genome and transcriptome signatures in colorectal cancers

Colorectal cancer is caused by a sequence of somatic genomic alterations affecting driver genes in core cancer pathways1. Here, to understand the functional and prognostic impact of cancer-causing somatic mutations, we analysed the whole genomes and transcriptomes of 1,063 primary colorectal cancers in a population-based cohort with long-term follow-up. From the 96 mutated driver genes, 9 were not previously implicated in colorectal cancer and 24 had not been linked to any cancer. Two distinct patterns of pathway co-mutations were observed, timing analyses identified nine early and three late driver gene mutations, and several signatures of colorectal-cancer-specific mutational processes were identified. Mutations in WNT, EGFR and TGFβ pathway genes, the mitochondrial CYB gene and 3 regulatory elements along with 21 copy-number variations and the COSMIC SBS44 signature correlated with survival. Gene expression classification yielded five prognostic subtypes with distinct molecular features, in part explained by underlying genomic alterations. Microsatellite-instable tumours divided into two classes with different levels of hypoxia and infiltration of immune and stromal cells. To our knowledge, this study constitutes the largest integrated genome and transcriptome analysis of colorectal cancer, and interlinks mutations, gene expression and patient outcomes. The identification of prognostic mutations and expression subtypes can guide future efforts to individualize colorectal cancer therapy.

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.

Prognostic value of structural variants in early breast cancer patients

Genomic analysis of structural variants(SVs) in breast cancer (BC) patients has been conducted, but the relationship between genomic alterations and BC prognosis remains unclear. We performed RNA sequencing of 297 early BC fresh-frozen tissues. We identified SVs using three tools (STAR.Arriba, STAR.fusion, and STAR.SEQR) with the COSMIC and Mitelman databases as guide references. We found a median of five to eight fusions per sample. In BC intrinsic subtypes, normal subtype had the fewest fusions (median: 1, interquartile range [IQR]: 0, 3) followed by luminal A (median: 5.5, IQR: 2.75, 10.25), luminal B (median: 9, IQR: 6, 16.5), HER2-enriched (median: 9, IQR: 6, 16.5) and basal (median 10, IQR: 6, 15.5) subtypes (p < 0.05). Intrachromosomal fusion was more frequent observed rather than interchromosomal fusion. In location, chromosome 17 had the most fusions followed by chromosome 1 and 11. When samples were divided into high and low fusion groups based on a cut-off value of 11 fusions, five-year event-free survival (5Y-EFS) was 68.1% in the high fusion group (n = 72) and 80.1% in the low fusion group (n = 125) (p = 0.024) while 75.6% among all patients (95% confidence interval: 0.699, 0.819). Among BC subtype, TNBCs with more fusions had shorter EFS compared to those with fewer fusions (5Y-EFS rate: 65.1% vs. 85.7%; p = 0.013) but no EFS differences were observed in other BC subtypes. ESTIMATE ImmuneScore was also associated with the number of fusions in TNBC (p < 0.005) and TNBCs with high ImmuneScore had better 5Y-EFS compared to those with low ImmuneScore (p = 0.041). In conclusion, diverse fusions were observed by BC subtype, and the number of fusions was associated with BC survival outcome and immune status in TNBC.

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.

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.

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.

A Protocol for the Detection of Fusion Transcripts Using RNA-Sequencing Data

Fusion transcripts are formed when two genes or their mRNAs fuse to produce a novel gene or chimeric transcript. Fusion genes are well-known cancer biomarkers used for cancer diagnosis and as therapeutic targets. Gene fusions are also found in normal physiology and lead to the evolution of novel genes that contribute to better survival and adaptation for an organism. Various in vitro approaches, such as FISH, PCR, RT-PCR, and chromosome banding techniques, have been used to detect gene fusion. However, all these approaches have low resolution and throughput. Due to the development of high-throughput next-generation sequencing technologies, the detection of fusion transcript becomes feasible using whole genome sequencing, RNA-Seq data, and bioinformatics tools. This chapter will overview the general computational protocol for fusion transcript detection from RNA-sequencing datasets.

Proteogenomic Workflow for Characterization of Microphthalmia Transcription Factor (MiT) Family Translocation Renal Cell Carcinoma

Microphthalmia transcription factor (MiT) family translocation renal cell carcinoma (tRCC) is a rare, aggressive, and heterogeneous subtype of kidney cancer, which is not well characterized. Since genetic alterations are always associated with carcinogenesis, and proteins are the major executors of biological features, multi-omics studies can reveal the systematic tRCC biological process comprehensively. Here, we describe the proteogenomic workflow for characterization of tRCC in detail to provide the knowledge foundation for integrated proteogenomic analysis of tRCC and other malignant tumors in the future.

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.