Chimeric Transcripts Research Articles

Replication Stress is an Actionable Genetic Vulnerability in Desmoplastic Small Round Cell Tumors.

Desmoplastic small round cell tumor (DSRCT) is an aggressive sarcoma subtype that is driven by the EWS-WT1 chimeric transcription factor. The prognosis for DSRCT is poor, and major advances in treating DSCRT have not occurred for over two decades. To identify effective therapeutic approaches to target DSRCT, we conducted a high-throughput drug sensitivity screen in a DSRCT cell line assessing chemosensitivity profiles for 79 small-molecule inhibitors. DSRCT cells were sensitive to PARP and ATR inhibitors (PARPi, ATRi), as monotherapies and in combination. These effects were recapitulated using multiple clinical PARPi and ATRi in three biologically distinct, clinically-relevant models of DSRCT, including cell lines, a patient-derived xenograft (PDX)-derived organoid model, and a cell line-derived xenograft mouse model. Mechanistically, exposure to a combination of PARPi and ATRi caused increased DNA damage, G2/M checkpoint activation, micronuclei accumulation, replication stress, and R-loop formation. EWS-WT1 silencing abrogated these phenotypes and was epistatic with exogenous expression of the R-loop resolution enzyme RNase H1 in reversing the sensitivity to PARPi and ATRi monotherapies. The combination of PARPi and ATRi also induced EWS-WT1-dependent cell-autonomous activation of the cGAS/STING innate immune pathway and cell surface expression of PD-L1. Taken together, these findings point towards a role for EWS-WT1 in generating R-loop-dependent replication stress that leads to a targetable vulnerability, providing a rationale for the clinical assessment of PARPi and ATRi in DSRCT.

Diversification of Antibodies: From V(D)J Recombination to Somatic Exon Shuffling.

Antibodies that gain specificity by a large insert encoding for an extra domain were described for the first time in 2016. In malaria-exposed individuals, an exon deriving from the leukocyte-associated immunoglobulin-like 1 (LAIR1) gene integrated via a copy-and-paste insertion into the immunoglobulin heavy chain encoding region. A few years later, a second example was identified, namely a dual exon integration from the leukocyte immunoglobulin-like receptor B1 (LILRB1) gene that is located in close proximity to LAIR1. A dedicated high-throughput characterization of chimeric immunoglobulin heavy chain transcripts unraveled, that insertions from distant genomic regions (including mitochondrial DNA) can contribute to human antibody diversity. This review describes the modalities of insert-containing antibodies. The role of known DNA mobility aspects, such as genomic translocation, gene conversion, and DNA fragility, is discussed in the context of insert-antibody generation. Finally, the review covers why insert antibodies were omitted from the past repertoire analyses and how insert antibodies can contribute to protective immunity or an autoreactive response.

Investigation of chimeric transcripts derived from LINE-1 and Alu retrotransposons in cerebellar tissues of individuals with autism spectrum disorder (ASD)

LINE-1 and Alu retrotransposons are components of the human genome and have been implicated in many human diseases. These elements can influence human transcriptome plasticity in various mechanisms. Chimeric transcripts derived from LINE-1 and Alu can also impact the human transcriptome, such as exonization and post-transcriptional modification. However, its specific role in ASD neuropathology remains unclear, particularly in the cerebellum tissues. We performed RNA-sequencing of post-mortem cerebellum tissues from ASD and unaffected individuals for transposable elements profiling and chimeric transcript identification. The majority of free transcripts of transposable elements were not changed in the cerebellum tissues of ASD compared with unaffected individuals. Nevertheless, we observed that chimeric transcripts derived from LINE-1 and Alu were embedded in the transcripts of differentially expressed genes in the cerebellum of ASD, and these genes were related to developments and abnormalities of the cerebellum. In addition, the expression levels of these genes were correlated with the significantly decreased thickness of the molecular layer in the cerebellum of ASD. We also found that global methylation and expression of LINE-1 and Alu elements were not changed in ASD, but observed in the ASD sub-phenotypes. Our findings showed associations between transposable elements and cerebellar abnormalities in ASD, particularly in distinct phenotypic subgroups. Further investigations using appropriate models are warranted to elucidate the structural and functional implications of LINE-1 and Alu elements in ASD neuropathology.

Transformation into acute myeloid leukemia with t(8;21)(q22;q22.1); RUNX1::RUNX1T1 from JAK2-mutated essential thrombocythemia: a case report

BackgroundBlast transformation is a rare but well-recognized event in Philadelphia-negative myeloproliferative neoplasms associated with a poor prognosis. Secondary acute myeloid leukemias evolving from myeloproliferative neoplasms are characterized by a unique set of cytogenetic and molecular features distinct from de novo disease. t(8;21) (q22;q22.1); RUNX1::RUNX1T1, one of the most frequent cytogenetic abnormalities in de novo acute myeloid leukemia, is rarely observed in post-myeloproliferative neoplasm acute myeloid leukemia. Here we report a case of secondary acute myeloid leukemia with t(8;21) evolving from JAK2-mutated essential thrombocythemia.Case presentationThe patient was a 74-year-old Japanese woman who was referred because of thrombocytosis (platelets 1046 × 109/L). Bone marrow was hypercellular with increase of megakaryocytes. Chromosomal analysis presented normal karyotype and genetic test revealed JAK2 V617F mutation. She was diagnosed with essential thrombocythemia. Thrombocytosis had been well controlled by oral administration of hydroxyurea; 2 years after the initial diagnosis with ET, she presented with leukocytosis (white blood cells 14.0 × 109/L with 82% of blasts), anemia (hemoglobin 91 g/L), and thrombocytopenia (platelets 24 × 109/L). Bone marrow was hypercellular and filled with 80% of myeloperoxidase-positive blasts bearing Auer rods. Chromosomal analysis revealed t(8;21) (q22;q22.1) and flow cytometry presented positivity of CD 13, 19, 34, and 56. Molecular analysis showed the coexistence of RUNX1::RUNX1T1 chimeric transcript and heterozygous JAK2 V617F mutation in leukemic blasts. She was diagnosed with secondary acute myeloid leukemia with t(8;21)(q22;q22.1); RUNX1::RUNX1T1 evolving from essential thrombocythemia. She was treated with combination chemotherapy with venetoclax and azacytidine. After the first cycle of the therapy, blasts disappeared from peripheral blood and decreased to 1.4% in bone marrow. After the chemotherapy, RUNX1::RUNX1T1 chimeric transcript disappeared, whereas mutation of JAK2 V617F was still present in peripheral leukocytes.ConclusionsTo our best knowledge, the present case is the first one with JAK2 mutation preceding the acquisition of t(8;21). Our result suggests that t(8;21); RUNX1::RUNX1T1 can be generated as a late event in the progression of JAK2-mutated myeloproliferative neoplasms. The case presented typical morphological and immunophenotypic features associated with t(8;21) acute myeloid leukemia.

Bioinformatic approaches for detection of fusion genes and <i>trans</i>-splicing products

Chimeric genes and transcripts can be biological markers as well as the reasons for tumor progression and development. Modern algorithms and high-throughput sequencing are the complementary clues to the question of the tumor origin and cancer detection as well as to the fundamental question of chimeric genes origin and their influence on molecular processes of the cell. A wide-range of algorithms for chimeric genes detection was developed, with various differences in computing speed, sensitivity, specificity, and focus on the experimental design. There exist three main types of bioinformatic approaches, which act according to the sequencing read length. Algorithms, which focus on short-read high-throughput sequencing (about 50–300 bр of read length) or long-read sequencing (about 5000–100000 bр of read length) exclusively or algorithms, which combine the results of both short and long-read sequencing. These algorithms are further subdivided into: 1) mapping-first approaches (STAR-Fusion, Arriba), which map reads to the genome or transcriptome directly and search the reads supporting the fused gene or transcript; 2) assembly-first approaches (Fusion-Bloom), which assemble the genome or transcriptome from the overlapping reads, and then compare the results to the reference transcriptome or genome to find transcripts or genes not present in the reference and therefore raising questions; 3) pseudoalignment approaches, which do not make local alignment, but just search for the closest transcript subsequence to the reads seed, following the precomputed index for all reference transcripts and provides the results. This article describes the main classes of available software tools for chimeric gene detection, provides the characteristics of these programs, their advantages and disadvantages. To date the most resource intensive and slowest are still assembly-first algorithms. Mapping-first approaches are quite fast and rather accurate at fusion detection, still the fastest and resource-saving are the pseudoalignment algorithms, but, worth noting, that the quick search is carried out at the expense of chimeras search quality decrease.

Expanding salivary biomarker detection by creating a synthetic neuraminic acid sensor via chimeragenesis.

Accurate and timely diagnosis of oral squamous cell carcinoma (OSCC) is crucial in preventing its progression to advanced stages with a poor prognosis. As such, the construction of sensors capable of detecting previously established disease biomarkers for the early and non-invasive diagnosis of this and many other conditions has enormous therapeutic potential. In this work, we apply synthetic biology techniques for the development of a whole-cell biosensor (WCB) that leverages the physiology of engineered bacteria in vivo to promote the expression of an observable effector upon detection of a soluble molecule. To this end, we have constructed a bacterial strain expressing a novel chimeric transcription factor (Sphnx) for the detection of N-acetylneuraminic acid (Neu5Ac), a salivary biomolecule correlated with the onset of OSCC. This WCB serves as the proof-of-concept of a platform that can eventually be applied to clinical screening panels for a multitude of oral and systemic medical conditions whose biomarkers are present in saliva.

Cancer fusion transcripts with human non-coding RNAs.

Cancer chimeric, or fusion, transcripts are thought to most frequently appear due to chromosomal aberrations that combine moieties of unrelated normal genes. When being expressed, this results in chimeric RNAs having upstream and downstream parts relatively to the breakpoint position for the 5'- and 3'-fusion components, respectively. As many other types of cancer mutations, fusion genes can be of either driver or passenger type. The driver fusions may have pivotal roles in malignisation by regulating survival, growth, and proliferation of tumor cells, whereas the passenger fusions most likely have no specific function in cancer. The majority of research on fusion gene formation events is concentrated on identifying fusion proteins through chimeric transcripts. However, contemporary studies evidence that fusion events involving non-coding RNA (ncRNA) genes may also have strong oncogenic potential. In this review we highlight most frequent classes of ncRNAs fusions and summarize current understanding of their functional roles. In many cases, cancer ncRNA fusion can result in altered concentration of the non-coding RNA itself, or it can promote protein expression from the protein-coding fusion moiety. Differential splicing, in turn, can enrich the repertoire of cancer chimeric transcripts, e.g. as observed for the fusions of circular RNAs and long non-coding RNAs. These and other ncRNA fusions are being increasingly recognized as cancer biomarkers and even potential therapeutic targets. Finally, we discuss the use of ncRNA fusion genes in the context of cancer detection and therapy.

Single-cell multiomics profiling reveals heterogeneous transcriptional programs and microenvironment in DSRCTs

Desmoplastic small round cell tumor (DSRCT) is a rare, aggressive sarcoma driven by the EWSR1::WT1 chimeric transcription factor. Despite this unique oncogenic driver, DSRCT displays a polyphenotypic differentiation of unknown causality. Using single-cell multi-omics on 12 samples from five patients, we find that DSRCT tumor cells cluster into consistent subpopulations with partially overlapping lineage- and metabolism-related transcriptional programs. Invitro modeling shows that high EWSR1::WT1 DNA-binding activity associates with most lineage-related states, in contrast to glycolytic and profibrotic states. Single-cell chromatin accessibility analysis suggests that EWSR1::WT1 binding site variability may drive distinct lineage-related transcriptional programs, supporting some level of cell-intrinsic plasticity. Spatial transcriptomics reveals that glycolytic and profibrotic states specifically localize within hypoxic niches at the periphery of tumor cell islets, suggesting an additional role of tumor cell-extrinsic microenvironmental cues. We finally identify a single-cell transcriptomics-derived epithelial signature associated with improved patient survival, highlighting the clinical relevance of our findings.

Transposable elements alter gene expression and may impact response to cisplatin therapy in ovarian cancer.

Cisplatin is widely employed for cancer treatment; therefore, understanding resistance to this drug is critical for therapeutic practice. While studies have delved into differential gene expression in the context of cisplatin resistance, findings remain somewhat scant. We performed a comprehensive investigation of transposable elements (TEs) expression and their impact in host genes in two cisplatin-treated ovarian cancer cell lines. RNA-seq, ATAC-seq, and in-depth bioinformatics analysis were used to compare cisplatin-sensitive and -resistant ovarian cancer cell lines. Our results reveal that cisplatin therapy alters not only the expression of protein-coding genes, but also key TEs, including LINE1, Alu, and endogenous retroviruses, in both cisplatin-sensitive and -resistant cell lines. By co-expressing with downstream genes or by creating chimeric transcripts with host genes at their insertion sites, these TEs seem to control the expression of protein-coding genes, including tumor-related genes. Our model uncovers TEs influencing the expression of cancer genes and cancer pathways. Collectively, our findings indicate that TE alterations associated with cisplatin treatment occur in critical cancer genes and cellular pathways synergically. This research highlights the importance of considering the entire spectrum of transcribed elements in the genome, especially TE expression, for a complete understanding of complex models like cancer response to treatment.

Fusion of Platelet Derived Growth Factor Receptor Alpha (PDGFRA) With Ubiquitin Specific Peptidase 8 (USP8) in a Calcified Chondroid Mesenchymal Neoplasm Harboring t(4;15)(q12;q21) as a Sole Aberration.

The term "calcified chondroid mesenchymal neoplasm" was introduced in 2021 to describe a group of tumors characterized by various morphological features, including the formation of cartilage or chondroid matrix. These tumors frequently carry chimeric genes where the 5'-end partner gene is fibronectin 1 and the 3'-end partner gene codes for receptor tyrosine kinase. Our study explores fusion of the genes platelet-derived growth factor receptor alpha (PDGFRA) and ubiquitin-specific peptidase 8 (USP8) in calcified chondroid mesenchymal neoplasm. Genetic investigations were conducted on a tumor located in the leg of a 71-year-old woman. G-banding analysis of short-term cultured tumor cells revealed the karyotype 46,XX,t(4;15)(q12;q21)[6]/46,XX[4]. RNA sequencing detected in-frame PDGFRA::USP8 and USP8::PDGFRA chimeric transcripts, which were validated by RT-PCR/Sanger sequencing. The PDGFRA::USP8 chimeric protein is predicted to have cell membrane location and functions as a chimeric ubiquitinyl hydrolase. The USP8::PDGFRA protein was predicted to be nuclear and function as a positive regulator of cellular metabolic process. We report, for the first time, a calcified chondroid mesenchymal neoplasm carrying a balanced t(4;15)(q12;q21) chromosomal translocation, resulting in the generation of both PDGFRA::USP8 and USP8::PDGFRA chimeras. The PDGFRA::USP8 protein is located on the cell membrane and functions as a chimeric ubiquitinyl hydrolase, activated by PDGFs. Conversely, USP8::PDGFRA is a nuclear protein regulating metabolic processes.

SAFB restricts contact domain boundaries associated with L1 chimeric transcription

Long interspersed element-1 (LINE-1 or L1) comprises 17% of the human genome, continuously generates genetic variations, and causes disease in certain cases. However, the regulation and function of L1 remain poorly understood. Here, we uncover that L1 can enrich RNA polymerase IIs (RNA Pol IIs), express L1 chimeric transcripts, and create contact domain boundaries in human cells. This impact of L1 is restricted by a nuclear matrix protein scaffold attachment factor B (SAFB) that recognizes transcriptionally active L1s by binding L1 transcripts to inhibit RNA Pol II enrichment. Acute inhibition of RNA Pol II transcription abolishes the domain boundaries associated with L1 chimeric transcripts, indicating a transcription-dependent mechanism. Deleting L1 impairs domain boundary formation, and L1 insertions during evolution have introduced species-specific domain boundaries. Our data show that L1 can create RNA Pol II-enriched regions that alter genome organization and that SAFB regulates L1 and RNA Pol II activity to preserve gene regulation.

Discovery of a polymorphic gene fusion via bottom-up chimeric RNA prediction.

Gene fusions and their chimeric products are commonly linked with cancer. However, recent studies have found chimeric transcripts in non-cancer tissues and cell lines. Large-scale efforts to annotate structural variations have identified gene fusions capable of generating chimeric transcripts even in normal tissues. In this study, we present a bottom-up approach targeting population-specific chimeric RNAs, identifying 58 such instances in the GTEx cohort, including notable cases such as SUZ12P1-CRLF3, TFG-ADGRG7and TRPM4-PPFIA3, which possess distinct patterns across different ancestry groups. We provide direct evidence for an additional 29 polymorphic chimeric RNAs with associated structural variants, revealing 13 novel rare structural variants. Additionally, we utilize the All of Us dataset and a large cohort of clinical samples to characterize the association of the SUZ12P1-CRLF3-causing variant with patient phenotypes. Our study showcases SUZ12P1-CRLF3 as a representative example, illustrating the identification of elusive structural variants by focusing on those producing population-specific fusion transcripts.

Immunophenotyping of blood and bone marrow cells as a way to search for differentiation syndrome risk factors in acute promyelocytic leukemia

Background. Differentiation syndrome (DS) is a potentially fatal complication of therapy for acute promyelocytic leukemia (APL) with an incidence of up to 48 %. To date, no reliable DS risk factors have been found, with the exception of leukocytosis at the APL onset.Aim. To determine the risk factors associated with DS in patients with APL during induction therapy with arsenic trioxide (ATO) and tretinoin (ATRA).Materials and methods. The study included 39 patients with APL, 29 (74.4 %) of them were classified as low-risk according to ELN (European Leukemia Net), 10 (25.6 %) were classified as high-risk. At the disease onset, cytological and molecular (chimeric transcript PML::RARα, FLT3-ITD mutation) bone marrow studies were performed, the expression of 28 differentiation antigens by blood and bone marrow blast cells was determined (markers of early precursors, myeloid and lymphoid differentiation, cell adhesion molecules, chemokine receptors, integrins, selectin), body mass index (BMI) and the leukocytes number dynamics during induction course were assessed. All patients received ATRA and ATO therapy. Patients from the high-risk group at the onset received 1–3 injections of idarubicin (12 mg/m2) and dexamethasone (8–10 mg/m2 2 times a day) to prevent DS until leukocytosis reduced. In cases of DS, dexamethasone was prescribed at a dose of 10 mg/m2 2 times a day; in cases of severe DS, the induction course was interrupted.Results. Of the 39 patients, 12 (30.8 %) were diagnosed with DS: 20 % of high-risk patients (2/10) and 34.5 % of low-risk patients (10/29). There was no statistically significant association of leukocytosis more than 10 × 109 /L at onset, microgranular morphology of blast cells, bcr3-variant PML::RARα, FLT3-ITD mutation with DS. In multivariate analysis, the probability of DS was associated with BMI ≥30 kg/m2 and mean fluorescence intensity of CD38 antigen by blast cells, regardless of risk group. based on the results of the ROC-analysis, the threshold value of mean CD38 fluorescence intensity was set at 25,000 cu, if exceeded, DS is highly likely to develop.Conclusion. the high incidence of DS among low-risk patients is probably due to the lack of prophylactic glucocorticosteroids administration for the development of leukocytosis during ATRA and ATO therapy. BMI ≥30 kg/m2 and mean CD38 fluorescence intensity more than 25,000 cu were identified as statistically significant DS risk factors.

Endogenous retroviral solo-LTRs in human genome.

Human endogenous retroviruses (HERVs) are derived from the infection and integration of exogenetic retroviruses. HERVs account for 8% of human genome, and the majority of HERVs are solitary LTRs (solo-LTRs) due to homologous recombination. Multiple findings have showed that solo-LTRs could provide an enormous reservoir of transcriptional regulatory sequences involved in diverse biological processes, especially carcinogenesis and cancer development. The link between solo-LTRs and human diseases still remains poorly understood. This review focuses on the regulatory modules of solo-LTRs, which contribute greatly to the diversification and evolution of human genes. More importantly, although inactivating mutations, insertions and deletions have been identified in solo-LTRs, the inherited regulatory elements of solo-LTRs initiate the expression of chimeric lncRNA transcripts, which have been reported to play crucial roles in human health and disease. These findings provide valuable insights into the evolutionary and functional mechanisms underlying the presence of HERVs in human genome. Taken together, in this review, we will present evidences showing the regulatory and encoding capacity of solo-LTRs as well as the significant impact on various aspects of human biology.

Abstract 143: Molecular profiling of tumor-initiation events in rhabdomyosarcoma uncovers mitochondrial vulnerabilities

Abstract Background: Alveolar Rhabdomyosarcoma (ARMS) is a very aggressive soft tissue tumor affecting children and adolescents. It is found that pediatric patients with RMS harboring chimeric PAX fusion transcription factors (TFs) exhibit a greater incidence of tumor relapse, metastasis, and poor survival outcome, thereby underscoring the urgent need to develop effective therapies to treat this subtype of childhood cancer. Methods: We generated a novel ARMS model by stable transduction of myogenic progenitor cells with a cassette encoding the exact PAX3/7-FOXO1 fusion protein found in ARMS tumors. We analyzed genome-wide expression levels with RNA-seq and utilized Cut&Run to map the binding sites of PAX3/7-FOXO1, active histone marks, enhancers, and super-enhancers in our novel model. Hi-C and promoter capture HiC (pCHiC) were also performed to generate 3D chromatin conformation maps in cells expressing the oncogenic TFs. We compared these data with analogous datasets from orthotopic patient-derived xenografts (O-PDXs) of PAX3/7-FOXO1 fusion-positive ARMS tumors from St. Jude and established ARMS cell lines to validate whether our model system recapitulates ARMS tumors. CRISPRi was performed to silence key targets in two representative ARMS cell lines and proliferation was scored after silencing. Next, we found that ARMS tumors displayed anomalies in mitochondrial mass and numbers, metabolic parameters, and membrane potential by performing TEM, Seahorse assays, and FACS-based TMRM and mitotracker Red dye uptake in a panel of ARMS cells and PDX samples. We performed cell viability and colony forming assays in ARMS cells treated with single drugs and combinations of drugs targeting FGFR4 and mitochondrial metabolism. Results: Our data provide evidence to support muscle progenitors as potential cells-of-origin for ARMS. We showed that rewiring of chromatin architecture mediated by PAX3/7-FOXO1 TFs is evident across established cell lines, primary tumors as well as O-PDXs. By integrating PAX3/7-FOXO1 binding, RNA-seq, and pCHiC data, we identified key promoter-enhancer (P-En) interactions that are erased or established de novo at multiple loci of master TFs like MYOD, MYOG and MYCN. We also identified key regulators of proliferation, including FGFR4, and mitochondrial function that are essential for tumor growth and invasiveness. Our data indicated aberrant and enhanced mitochondrial metabolism in ARMS as well as selective sensitivity to mitochondrial inhibitors in RMS cells. Administration of these inhibitors preferentially inhibited cell viability and clonogenicity in ARMS, and prompted a pre-clinical phase II study to test these inhibitors at St. Jude Research Hospital. Conclusion: We present a novel treatment strategy leveraging the mitochondrial dependencies of fusion positive PAX-driven tumors which can be exploited for precision therapy in high-risk ARMS. Citation Format: Bhargab Kalita, Gerard Martinez-Cebrian, Justina McEvoy, Brittney Gordon, Kaley Blankenship, Asa Karlstorm, Elizabeth Stewart, Michael Dyer, Brian Dynlacht. Molecular profiling of tumor-initiation events in rhabdomyosarcoma uncovers mitochondrial vulnerabilities [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 143.

Abstract 2861: Rhabdomyosarcoma fusion-oncogene alters the chromatin landscape to initially drive a neural signature in vivo

Abstract Fusion-positive rhabdomyosarcoma is an aggressive pediatric cancer lacking curative therapies, and outcomes for children with this disease have not improved in decades. The chimeric transcription factor PAX3::FOXO1 is the most common and lethal driver of fusion-positive rhabdomyosarcoma, and it consists of the DNA binding domains of PAX3 fused to the transactivation domain of FOXO1. Despite its well-established essentiality for tumorigenesis, how PAX3::FOXO1 initially generates a tumorigenic cell state has been challenging to study due to the lack of in vivo models expressing the fusion-oncogene at accessible early developmental time points. Here, we developed a novel zebrafish mRNA injection model that ubiquitously expresses PAX3::FOXO1 during gastrulation representing all three germ layers. Using this high-throughput model, we used ChIP-seq to identify initial in vivo PAX3::FOXO1 binding sites and targets. PAX3::FOXO1 has pioneering activity, where it can utilize partial-motif recognition to bind to inaccessible chromatin and alter chromatin structure. With biochemical and 2D chromatin sequencing approaches in our model, we find evidence for in vivo PAX3::FOXO1 pioneering activity. PAX3::FOXO1 consists of two DNA binding domains, a paired and homeobox domain, and we demonstrate a new mode of pioneering activity driven by its homeobox domain given a strong enrichment of homeobox-related motifs at PAX3::FOXO1 binding sites. PAX3::FOXO1 activity in developing zebrafish embryos results in an arrested development phenotype and transcriptional signatures enriched for pathways such as an inhibition of segmentation and myogenic development. Critically, we find that PAX3::FOXO1 directly activates neural-related gene targets by increasing chromatin accessibility and re-distribution of the active histone mark H3K27Ac to these loci. These discoveries are particularly striking because rhabdomyosarcoma is traditionally associated with skeletal muscle characteristics. With our findings, we hypothesize neural gene activation is a critical mechanism for fusion-positive rhabdomyosarcoma tumor initiation and could be a critical pathway across various stages of rhabdomyosarcoma tumorigenesis. This hypothesis agrees with our previous discovery of a novel PAX3::FOXO1 cooperating neural transcription factor, HES3, which we find is directly activated by PAX3::FOXO1 in our mRNA injection model. Our long-term goal is to utilize this versatile model to functionally evaluate cooperation between PAX3::FOXO1 and both canonical rhabdomyosarcoma and neural transcription factors. Our approach will employ a cross-species comparative analysis across this injection model, our zebrafish tumor models, other model systems, and patient data. This work will identify conserved targets and mechanisms of PAX3::FOXO1 tumorigenesis that could become new therapeutic vulnerabilities. Citation Format: Jack Kucinski, Cenny Taslim, Alexi Tallan, Matthew Cannon, Katherine Silvius, Benjamin Stanton, Genevieve Kendall. Rhabdomyosarcoma fusion-oncogene alters the chromatin landscape to initially drive a neural signature in vivo [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 2861.

Abstract 2851: Characterization of neogenes in desmoplastic small round cell tumors

Abstract Desmoplastic Small Round Cell Tumor (DSRCT) is a rare, soft tissue sarcoma caused by the EWS::WT1 fusion protein (FP). The FP acts as an aberrant transcription factor with oncogenic properties. Recently, oncogenic chimeric transcription factors (OCTF), such as EWS::WT1, were discovered to induce the expression of novel transcripts in otherwise silent genomic regions. These neogenes (NGs) are highly specific to each OCTF and not detected in normal tissue. 37 NGs were found to be associated with the EWS::WT1 FP.We explored the NGs to characterize their expression and regulation in a cohort of MD Anderson Cancer Center patients. We annotated and counted the NGs from the bulk RNA-seq data of 37 specimens. After fusion detection, only 22/37 specimens harbored the EWS::WT1 FP. We found that only the same 22 specimens robustly expressed DSRCT-specific NGs. Interestingly the NGs were heterogeneously expressed, suggesting that there may be patient-specific differences. Single-nucleus RNA-sequencing (snRNA-seq) of Ewing sarcoma (ES), CIC::DUX4 sarcoma (CDS), DSRCT, and osteosarcoma revealed that only the three sarcomas driven by an OCTF robustly expressed NGs. Notably, the DSRCT-specific NGs were only expressed in DSRCT. Similarly, ES and CDS only expressed their subtype-specific NGs. Given the robust expression of the DSRCT-specific NGs in DSRCT, we asked if the EWS::WT1 FP regulated the expression of the NGs. Knockdown (KD) of the EWS::WT1 FP in four DSRCT cell lines showed a consistent reduction in DSRCT-specific NGs. While there was variability in NG expression and EWS::WT1 FP KD sensitivity, we identified seven NGs decreased by at least 75% in all four cell lines and an additional eight NGs with 75% reduction across three cell lines. BOD-DSRCT and SK-DSRCT2 had the largest number of NGs reduced after KD, with 25/37 and 20/37 NGs with 75% expression reduction, respectively. EWS::WT1 FP has two known isoforms generated by alternative splicing in the WT1 domain that includes or excludes three amino acids, lysine, threonine, and serine (KTS) between zinc fingers 3 and 4, producing E+KTS and E-KTS isoforms. By transducing one or both isoforms into the mesothelial cell line LP9, we found 3/37 NGs were induced by the E+KTS isoform and 17/37 were induced by the E-KTS isoform, suggesting that the E-KTS isoform is primarily responsible for NG expression. Current work is underway investigating the relationship between EWS::WT1 FP binding and NG expression and the potential for NG to produce novel peptides that could serve as DSRCT immunotherapy targets. Citation Format: Danh Truong, Justin Magrath, Sean Lee, Joseph Ludwig. Characterization of neogenes in desmoplastic small round cell tumors [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 2851.

Abstract 4739: Targeting resistance to BET inhibitors in Fusion-Positive Rhabdomyosarcoma

Abstract Rhabdomyosarcoma (RMS) is the most frequent pediatric soft tissue sarcoma. Fusion-positive (FP-)RMS expressing the oncogenic chimeric transcription factor (TF) PAX3-FOXO1 (P3F) is at high risk of recurrence. FP-RMS cells show survival dependency on P3F. Despite P3F, as a TF, is considered undruggable, we and others demonstrated that P3F levels/functions can be epigenetically modulated. The BET protein and epigenetic reader BRD4 binds super-enhancers (SEs) to foster oncogenic transcription in cancer. We showed that P3F rewires the RMS enhancer landscape by recruiting BRD4 on oncogenic SEs. Accordingly, JQ1, a BET inhibitor (BETi) shuts down P3F functions halting tumor growth in vitro and in vivo. Therefore, BRD4 inhibition in FP-RMS results in a tumor subtype-specific vulnerability. Nonetheless, as for other targeted therapies, cancer cells can acquire resistance to BETi, thus, we investigated whether FP-RMS cells could develop BETi resistance. Therefore, we established resistant cells from two sensitive FP-RMS cell lines chronically exposed to escalating doses of JQ1, as a model of naturally acquired resistance. Our data show that the two established FP-RMS cell lines, acquired resistance by increasing JQ1 IC50 over the sensitive cells both in 2D and 3D settings. Compared to the sensitive ones, both the resistant cell lines show cross-resistance to other monovalent and bivalent BETi and degraders. Resistant cells also maintain the ability to form colonies in a clonogenic assay and to invade in a 3D Matrigel assay in the presence of JQ1. Moreover, conversely to sensitive cells, JQ1 treatment is unable to down-regulate P3F, MYCN and BCL2 (P3F target) levels in resistant cells. Furthermore, integration of transcriptomic and proteomic data reveals differentially modulated subset of genes/proteins among which MYCN and P3F targets, suggesting a restored activity of the TFs circuitry in resistant cells. Accordingly, MYCN silencing completely halts growth of resistant cells. In line, MYCN protein levels are reduced by JQ1 in a dose-dependent manner in sensitive but not in resistant cells and its half-life is 3 folds longer in resistant than in sensitive cells. Being MYCN protein stability affected by specific post-translational modifications among which the phosphorylation at Ser62 by ERKs, we evaluated the expression of phosphorylated MYCN and the activated form of ERK (pERK) and found both of them increased in resistant cells. In agreement, resistant cells are more sensitive to MEK-ERKs pathway inhibition, showing a lower IC50 of Trametinib (MEKi) compared to sensitive cells. Moreover, Trametinib and JQ1 co-treatment downregulates MYCN protein levels in resistant but not in sensitive cells, suggesting that ERKs overactivation could participate to BETi acquired resistance by stabilizing MYCN, thus representing an acquired vulnerability. The study has received funding from Ministero della Salute to SP. Citation Format: Matteo Cassandri, Erika Ferraro, Lucrezia D'Archivio, Francesca Aiello, Hsien-Chao Chou, Young K. Song, Berkley E. Gryder, Jun S. Wei, Simone Sidoli, Robert G. Hawley, Rossella Rota, Franco Locatelli, Javed Khan, Silvia Pomella. Targeting resistance to BET inhibitors in Fusion-Positive Rhabdomyosarcoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 4739.

Abstract 138: UCHL5 and PAX3-FOXO1 resolve replication stress in rhabdomyosarcoma

Abstract PAX3-FOXO1 positive rhabdomyosarcoma (RMS) is an aggressive pediatric sarcoma with a terrible prognosis. The resultant fusion acts as a chimeric transcription factor, in part remodeling chromatin to inhibit terminal differentiation and maintain cancer cell proliferation. As this oncoprotein is not directly druggable, we have used genetic screens to identify dependencies downstream of PAX3-FOXO1 to improve our understanding of pathogenesis and develop new therapeutic targets. We performed pooled CRISPRi screening in an isogenic system of a patient-derived RMS cell line with or without knockdown of PAX3-FOXO1. Loss of the ubiquitin hydroxylase UCHL5 led to selective depletion of PAX3-FOXO1 positive Rh30 cells compared to PAX3-FOXO1 knockdown cells, suggesting that the fusion creates UCHL5 dependence. We validated this finding in competition growth assays in additional fusion positive (FP) and fusion negative (FN) RMS cell lines, and observed UCHL5 dependency only in the former. By contrast, UCHL5 was dispensable for cell growth and survival in FN RMS cell lines. These data were recapitulated in the large DepMap dataset, which also identified the top genetic co-dependencies of UCHL5 as INO80-complex (INO80c) members. We therefore investigated whether the fusion selective effects of UCHL5 were mediated by the INO80c, of which UCHL5 is an N-terminal member. We observed that knockdown of another N-terminal component of INO80c, NFRKB, phenocopies the fusion-selective effects of UCHL5 loss. Conversely, the catalytic ATPase INO80 was pan-essential, suggesting that only N-terminal subunits of INO80c regulate PAX3-FOXO1 specific functions. INO80c plays an essential role in mitigating DNA replication stress. Indeed, loss of UCHL5 led to DNA damage as measured by γH2AX and slowing of replication forks as measured by fiber assays in FP, but not FN RMS. Unexpectedly, acute degradation of PAX3-FOXO1 also slowed DNA replication fork speed, suggesting a genetic model in which UCHL5 and PAX3-FOXO1 cooperate to resolve replication stress in FP RMS. Finally, in vivo experiments using cell line xenografts confirms that loss of UCHL5 slows the growth of FP RMS, but not FN RMS. These findings suggest that loss of UCHL5 induces replication stress in a fusion-specific manner, and offers both an opportunity to study the mechanisms by which the PAX3-FOXO1 oncoprotein mitigates replication stress, and to develop UCHL5 inhibitors as precision therapies for FP RMS patients. Citation Format: Amit J. Sabnis, Pushpendra K. Sahu, Yue Pan. UCHL5 and PAX3-FOXO1 resolve replication stress in rhabdomyosarcoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 138.

Abstract 7418: LongFuse: Detecting gene fusion transcripts from high throughput long-read single cell RNA sequencing data

Abstract Gene fusion resulted from genomic translocation plays important oncogenic roles including but not limited to gene activation and changing protein coding sequences. Historically, fusion gene analysis has been limited to bulk levels resulting in mixed oncogene activation models. Single cell Nanopore RNA sequencing (scNanoRNAseq) has the unique strength in detecting full-length chimeric transcripts in thousands of single cells in parallel. However, robust computational tools for detecting fusion genes in single cells from scNanoRNAseq data are missing. Here, we developed LongFuse, a fast and accurate tool, to detect fusion genes and their splicing isoforms in single cells from scNanoRNAseq data. LongFuse implements an XOR logic gate algorithm to speed up fusion candidate searching and stringent criteria to eliminate false discoveries due to trans-splicing events and randomly ligated chimeras. We applied LongFUSE onto scNanoRNAseq data of both human prostate cancer cell line, VCaP, and tumor tissues with known TMPRSS2-ERG fusion. Our results showed that both cancer cell line and prostate tumors harbored a mixture of fusion positive and fusion negative tumor cells. Gene expression analysis revealed transcriptional differences between fusion positive and negative cells including genes involved in cancer cell proliferation and invasion programs such as ELL2, ZBTB16, IGF1R, ANK3, MYPOP and PLPP1 genes. LongFUSE also calculates the splicing isoforms of fusion gene transcripts, allowing prediction of protein sequences of fusion gene transcripts. Taken together, our data demonstrated that LongFUSE is a robust computational tool for detecting gene fusions and their splicing isoforms at single cell level from scNanoRNAseq data. Citation Format: Cheng-Kai Shiau, Yueying He, Hsiao-Yun Lin, Lina Lu, Xiao-Dong Lu, Jindan Yu, Ruli Gao. LongFuse: Detecting gene fusion transcripts from high throughput long-read single cell RNA sequencing data [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 7418.