CAR-T cell expansion platforms yield distinct T cell differentiation states.

With investigators looking to expand engineered T cell therapies such as CAR-T to new tumor targets and patient populations, a variety of cell manufacturing platforms have been developed to scale manufacturing capacity using closed and/or automated systems. Such platforms are particularly useful for solid tumor targets, which typically require higher CAR-T cell doses. Although T cell phenotype and function are key attributes that often correlate with therapeutic efficacy, how manufacturing platforms influence the final CAR-T cell product is currently unknown. We compared 4 commonly used T cell manufacturing platforms (CliniMACS Prodigy, Xuri W25 rocking platform, G-Rex gas-permeable bioreactor, static bag culture) using identical media, stimulation, culture length, and donor starting material. Selected CD4+CD8+ cells were transduced with lentiviral vector incorporating a CAR targeting FGFR4, a promising target for pediatric sarcoma. We observed significant differences in overall expansion over the 14-day culture; bag cultures had the highest capacity for expansion while the Prodigy had the lowest (481-fold versus 84-fold, respectively). Strikingly, we also observed considerable differences in the phenotype of the final product, with the Prodigy significantly enriched for CCR7+CD45RA+ naïve/stem central memory (Tn/scm)-like cells at 46% compared to bag and G-Rex with 16% and 13%, respectively. Gene expression analysis also showed that Prodigy CAR-Ts are more naïve, less cytotoxic and less exhausted than bag, G-Rex, and Xuri CAR-Ts, and pointed to differences in cell metabolism that were confirmed via metabolic assays. We hypothesized that dissolved oxygen level, which decreased substantially during the final 3 days of the Prodigy culture, may contribute to the observed differences in T cell phenotype. By culturing bag and G-Rex cultures in 1% O2 from day 5 onward, we could generate >60% Tn/scm-like cells, with longer time in hypoxia correlating with a higher percentage of Tn/scm-like cells. Intriguingly, our results suggest that oxygenation is responsible, at least in part, for observed differences in T cell phenotype among bioreactors and suggest hypoxic culture as a potential strategy prevent T cell differentiation during expansion. Ultimately, our study demonstrates that selection of bioreactor system may have profound effects not only on the capacity for expansion, but also on the differentiation state of the resulting CAR-T cells.

Lineage specific transcription factor waves reprogram neuroblastoma from self-renewal to differentiation.

Temporal regulation of super-enhancer (SE) driven transcription factors (TFs) underlies normal developmental programs. Neuroblastoma (NB) arises from an inability of sympathoadrenal progenitors to exit a self-renewal program and terminally differentiate. To identify SEs driving TF regulators, we use all-trans retinoic acid (ATRA) to induce NB growth arrest and differentiation. Time-course H3K27ac ChIP-seq and RNA-seq reveal ATRA coordinated SE waves. SEs that decrease with ATRA link to stem cell development (MYCN, GATA3, SOX11). CRISPR-Cas9 and siRNA verify SOX11 dependency, in vitro and in vivo. Silencing the SOX11 SE using dCAS9-KRAB decreases SOX11 mRNA and inhibits cell growth. Other TFs activate in sequential waves at 2, 4 and 8 days of ATRA treatment that regulate neural development (GATA2 and SOX4). Silencing the gained SOX4 SE using dCAS9-KRAB decreases SOX4 expression and attenuates ATRA-induced differentiation genes. Our study identifies oncogenic lineage drivers of NB self-renewal and TFs critical for implementing a differentiation program.

KDM3B inhibitors disrupt the oncogenic activity of PAX3-FOXO1 in fusion-positive rhabdomyosarcoma.

Fusion-positive rhabdomyosarcoma (FP-RMS) is an aggressive pediatric sarcoma driven primarily by the PAX3-FOXO1 fusion oncogene, for which therapies targeting PAX3-FOXO1 are lacking. Here, we screen 62,643 compounds using an engineered cell line that monitors PAX3-FOXO1 transcriptional activity identifying a hitherto uncharacterized compound, P3FI-63. RNA-seq, ATAC-seq, and docking analyses implicate histone lysine demethylases (KDMs) as its targets. Enzymatic assays confirm the inhibition of multiple KDMs with the highest selectivity for KDM3B. Structural similarity search of P3FI-63 identifies P3FI-90 with improved solubility and potency. Biophysical binding of P3FI-90 to KDM3B is demonstrated using NMR and SPR. P3FI-90 suppresses the growth of FP-RMS in vitro and in vivo through downregulating PAX3-FOXO1 activity, and combined knockdown of KDM3B and KDM1A phenocopies P3FI-90 effects. Thus, we report KDM inhibitors P3FI-63 and P3FI-90 with the highest specificity for KDM3B. Their potent suppression of PAX3-FOXO1 activity indicates a possible therapeutic approach for FP-RMS and other transcriptionally addicted cancers.

A SNAI2/CTCF Interaction is Required for NOTCH1 Expression in Rhabdomyosarcoma.

Rhabdomyosarcoma (RMS) is a pediatric malignancy of the muscle with characteristics of cells blocked in differentiation. NOTCH1 is an oncogene that promotes self-renewal and blocks differentiation in the fusion negative-RMS sub-type. However, how NOTCH1 expression is transcriptionally maintained in tumors is unknown. Analyses of SNAI2 and CTCF chromatin binding and HiC analyses revealed a conserved SNAI2/CTCF overlapping peak downstream of the NOTCH1 locus marking a sub-topologically associating domain (TAD) boundary. Deletion of the SNAI2-CTCF peak showed that it is essential for NOTCH1 expression and viability of FN-RMS cells. Reintroducing constitutively activated NOTCH1-ΔE in cells with the SNAI2-CTCF peak deleted restored cell-viability. Ablation of SNAI2 using CRISPR/Cas9 reagents resulted in the loss of majority of RD and SMS-CTR FN-RMS cells. However, the few surviving clones that repopulate cultures have recovered NOTCH1. Cells that re-establish NOTCH1 expression after SNAI2 ablation are unable to differentiate robustly as SNAI2 shRNA knockdown cells; yet, SNAI2-ablated cells continued to be exquisitely sensitive to ionizing radiation. Thus, we have uncovered a novel mechanism by which SNAI2 and CTCF maintenance of a sub-TAD boundary promotes rather than represses NOTCH1 expression. Further, we demonstrate that SNAI2 suppression of apoptosis post-radiation is independent of SNAI2/NOTCH1 effects on self-renewal and differentiation.

Piperacetazine Directly Binds to the PAX3::FOXO1 Fusion Protein and Inhibits Its Transcriptional Activity.

The tumor-specific chromosomal translocation product, PAX3::FOXO1, is an aberrant fusion protein that plays a key role for oncogenesis in the alveolar subtype of rhabdomyosarcoma (RMS). PAX3::FOXO1 represents a validated molecular target for alveolar RMS and successful inhibition of its oncogenic activity is likely to have significant clinical applications. Even though several PAX3::FOXO1 function-based screening studies have been successfully completed, a directly binding small-molecule inhibitor of PAX3::FOXO1 has not been reported. Therefore, we screened small-molecule libraries to identify compounds that were capable of directly binding to PAX3::FOXO1 protein using surface plasmon resonance technology. Compounds that directly bound to PAX3::FOXO1 were further evaluated in secondary transcriptional activation assays. We discovered that piperacetazine can directly bind to PAX3::FOXO1 protein and inhibit fusion protein-derived transcription in multiple alveolar RMS cell lines. Piperacetazine inhibited anchorage-independent growth of fusion-positive alveolar RMS cells but not embryonal RMS cells. On the basis of our findings, piperacetazine is a molecular scaffold upon which derivatives could be developed as specific inhibitors of PAX3::FOXO1. These novel inhibitors could potentially be evaluated in future clinical trials for recurrent or metastatic alveolar RMS as novel targeted therapy options. SIGNIFICANCE: RMS is a malignant soft-tissue tumor mainly affecting the pediatric population. A subgroup of RMS with worse prognosis harbors a unique chromosomal translocation creating an oncogenic fusion protein, PAX3::FOXO1. We identified piperacetazine as a direct inhibitor of PAX3::FOXO1, which may provide a scaffold for designing RMS-specific targeted therapy.

Preclinical development of a chimeric antigen receptor T cell therapy targeting FGFR4 in rhabdomyosarcoma.

Pediatric patients with relapsed or refractory rhabdomyosarcoma (RMS) have dismal cure rates, and effective therapy is urgently needed. The oncogenic receptor tyrosine kinase fibroblast growth factor receptor 4 (FGFR4) is highly expressed in RMS and lowly expressed in healthy tissues. Here, we describe a second-generation FGFR4-targeting chimeric antigen receptor (CAR), based on an anti-human FGFR4-specific murine monoclonal antibody 3A11, as an adoptive T cell treatment for RMS. The 3A11 CAR T cells induced robust cytokine production and cytotoxicity against RMS cell lines in vitro. In contrast, a panel of healthy human primary cells failed to activate 3A11 CAR T cells, confirming the selectivity of 3A11 CAR T cells against tumors with high FGFR4 expression. Finally, we demonstrate that 3A11 CAR T cells are persistent in vivo and can effectively eliminate RMS tumors in two metastatic and two orthotopic models. Therefore, our study credentials CAR T cell therapy targeting FGFR4 to treat patients with RMS.

A stem cell epigenome is associated with primary nonresponse to CD19 CAR T cells in pediatric acute lymphoblastic leukemia.

CD19 chimeric antigen receptor T-cell therapy (CD19-CAR) has changed the treatment landscape and outcomes for patients with pre-B-cell acute lymphoblastic leukemia (B-ALL). Unfortunately, primary nonresponse (PNR), sustained CD19+ disease, and concurrent expansion of CD19-CAR occur in 20% of the patients and is associated with adverse outcomes. Although some failures may be attributable to CD19 loss, mechanisms of CD19-independent, leukemia-intrinsic resistance to CD19-CAR remain poorly understood. We hypothesize that PNR leukemias are distinct compared with primary sensitive (PS) leukemias and that these differences are present before treatment. We used a multiomic approach to investigate this in 14 patients (7 with PNR and 7 with PS) enrolled in the PLAT-02 trial at Seattle Children's Hospital. Long-read PacBio sequencing helped identify 1 PNR in which 47% of CD19 transcripts had exon 2 skipping, but other samples lacked CD19 transcript abnormalities. Epigenetic profiling discovered DNA hypermethylation at genes targeted by polycomb repressive complex 2 (PRC2) in embryonic stem cells. Similarly, assays of transposase-accessible chromatin-sequencing revealed reduced accessibility at these PRC2 target genes, with a gain in accessibility of regions characteristic of hematopoietic stem cells and multilineage progenitors in PNR. Single-cell RNA sequencing and cytometry by time of flight analyses identified leukemic subpopulations expressing multilineage markers and decreased antigen presentation in PNR. We thus describe the association of a stem cell epigenome with primary resistance to CD19-CAR therapy. Future trials incorporating these biomarkers, with the addition of multispecific CAR T cells targeting against leukemic stem cell or myeloid antigens, and/or combined epigenetic therapy to disrupt this distinct stem cell epigenome may improve outcomes of patients with B-ALL.

Predicting Molecular Subtype and Survival of Rhabdomyosarcoma Patients Using Deep Learning of H&E Images: A Report from the Children's Oncology Group.

PURPOSE: Rhabdomyosarcoma (RMS) is an aggressive soft-tissue sarcoma, which primarily occurs in children and young adults. We previously reported specific genomic alterations in RMS, which strongly correlated with survival; however, predicting these mutations or high-risk disease at diagnosis remains a significant challenge. In this study, we utilized convolutional neural networks (CNN) to learn histologic features associated with driver mutations and outcome using hematoxylin and eosin (H&E) images of RMS. EXPERIMENTAL DESIGN: Digital whole slide H&E images were collected from clinically annotated diagnostic tumor samples from 321 patients with RMS enrolled in Children's Oncology Group (COG) trials (1998-2017). Patches were extracted and fed into deep learning CNNs to learn features associated with mutations and relative event-free survival risk. The performance of the trained models was evaluated against independent test sample data (n = 136) or holdout test data. RESULTS: The trained CNN could accurately classify alveolar RMS, a high-risk subtype associated with PAX3/7-FOXO1 fusion genes, with an ROC of 0.85 on an independent test dataset. CNN models trained on mutationally-annotated samples identified tumors with RAS pathway with a ROC of 0.67, and high-risk mutations in MYOD1 or TP53 with a ROC of 0.97 and 0.63, respectively. Remarkably, CNN models were superior in predicting event-free and overall survival compared with current molecular-clinical risk stratification. CONCLUSIONS: This study demonstrates that high-risk features, including those associated with certain mutations, can be readily identified at diagnosis using deep learning. CNNs are a powerful tool for diagnostic and prognostic prediction of rhabdomyosarcoma, which will be tested in prospective COG clinical trials.

An optimized bicistronic chimeric antigen receptor against GPC2 or CD276 overcomes heterogeneous expression in neuroblastoma.

Chimeric antigen receptor (CAR) T cell therapies targeting single antigens have performed poorly in clinical trials for solid tumors due to heterogenous expression of tumor-associated antigens (TAAs), limited T cell persistence, and T cell exhaustion. Here, we aimed to identify optimal CARs against glypican 2 (GPC2) or CD276 (B7-H3), which were highly but heterogeneously expressed in neuroblastoma (NB), a lethal extracranial solid tumor of childhood. First, we examined CAR T cell expansion in the presence of targets by digital droplet PCR. Next, using pooled competitive optimization of CAR by cellular indexing of transcriptomes and epitopes by sequencing (CITE-Seq), termed P-COCC, we simultaneously analyzed protein and transcriptome expression of CAR T cells to identify high-activity CARs. Finally, we performed cytotoxicity assays to identify the most effective CAR against each target and combined the CARs into a bicistronic "OR" CAR (BiCisCAR). BiCisCAR T cells effectively eliminated tumor cells expressing GPC2 or CD276. Furthermore, the BiCisCAR T cells demonstrated prolonged persistence and resistance to exhaustion when compared with CARs targeting a single antigen. This study illustrated that targeting multiple TAAs with BiCisCAR may overcome heterogenous expression of target antigens in solid tumors and identified a potent, clinically relevant CAR against NB. Moreover, our multimodal approach integrating competitive expansion, P-COCC, and cytotoxicity assays is an effective strategy to identify potent CARs among a pool of candidates.

Dentithecamides A-H, Diacylated Zoanthoxanthin Derivatives with PAX3-FOXO1 Inhibitory Activity from the Hydroid Dentitheca habereri.

Chemical investigation of the marine hydroid Dentitheca habereri led to the identification of eight new diacylated zoanthoxanthin alkaloids, named dentithecamides A-H (1-8), along with three previously reported analogues, zoamides B-D (9-11). The structures of compounds 1-11 were elucidated by spectroscopic and spectrometric analyses, including IR, HRESIMS, and NMR experiments, and by comparison with literature data. Compounds 1-11 are the first zoanthoxanthin alkaloids to be reported from a hydroid. Dentithecamides A (1) and B (2) along with zoamides B-D (9-11), which all share a conformationally mobile cycloheptadiene core, inhibited PAX3-FOXO1 regulated transcriptional activity and thus provided a structural framework for the potential development of more potent PAX3-FOXO1 inhibitors.

Tumor Mutation Burden, Expressed Neoantigens and the Immune Microenvironment in Diffuse Gliomas.

BACKGROUND: A consistent correlation between tumor mutation burden (TMB) and tumor immune microenvironment has not been observed in gliomas as in other cancers. METHODS: Driver germline and somatic mutations, TMB, neoantigen, and immune cell signatures were analyzed using whole exome sequencing (WES) and transcriptome sequencing of tumor and WES of matched germline DNA in a cohort of 66 glioma samples (44 IDH-mutant and 22 IDH-wildtype). RESULTS: Fourteen samples revealed a hypermutator phenotype (HMP). Eight pathogenic (P) or likely pathogenic (LP) germline variants were detected in 9 (19%) patients. Six of these 8 genes were DNA damage repair genes. P/LP germline variants were found in 22% of IDH-mutant gliomas and 12.5% of IDH-wildtype gliomas (p = 0.7). TMB was correlated with expressed neoantigen but showed an inverse correlation with immune score (R = -0.46, p = 0.03) in IDH-wildtype tumors and no correlation in IDH-mutant tumors. The Antigen Processing and Presentation (APP) score correlated with immune score and was surprisingly higher in NHMP versus HMP samples in IDH-wildtype gliomas, but higher in HMP versus NHMP in IDH-mutant gliomas. CONCLUSION: TMB was inversely correlated with immune score in IDH-wildtype gliomas and showed no correlation in IDH-mutant tumors. APP was correlated with immune score and may be further investigated as a biomarker for response to immunotherapy in gliomas. Studies of germline variants in a larger glioma cohort are warranted.

BAF complexes drive proliferation and block myogenic differentiation in fusion-positive rhabdomyosarcoma.

Rhabdomyosarcoma (RMS) is a pediatric malignancy of skeletal muscle lineage. The aggressive alveolar subtype is characterized by t(2;13) or t(1;13) translocations encoding for PAX3- or PAX7-FOXO1 chimeric transcription factors, respectively, and are referred to as fusion positive RMS (FP-RMS). The fusion gene alters the myogenic program and maintains the proliferative state while blocking terminal differentiation. Here, we investigated the contributions of chromatin regulatory complexes to FP-RMS tumor maintenance. We define the mSWI/SNF functional repertoire in FP-RMS. We find that SMARCA4 (encoding BRG1) is overexpressed in this malignancy compared to skeletal muscle and is essential for cell proliferation. Proteomic studies suggest proximity between PAX3-FOXO1 and BAF complexes, which is further supported by genome-wide binding profiles revealing enhancer colocalization of BAF with core regulatory transcription factors. Further, mSWI/SNF complexes localize to sites of de novo histone acetylation. Phenotypically, interference with mSWI/SNF complex function induces transcriptional activation of the skeletal muscle differentiation program associated with MYCN enhancer invasion at myogenic target genes, which is recapitulated by BRG1 targeting compounds. We conclude that inhibition of BRG1 overcomes the differentiation blockade of FP-RMS cells and may provide a therapeutic strategy for this lethal childhood tumor.

Immuno-transcriptomic profiling of extracranial pediatric solid malignancies.

We perform an immunogenomics analysis utilizing whole-transcriptome sequencing of 657 pediatric extracranial solid cancer samples representing 14 diagnoses, and additionally utilize transcriptomes of 131 pediatric cancer cell lines and 147 normal tissue samples for comparison. We describe patterns of infiltrating immune cells, T cell receptor (TCR) clonal expansion, and translationally relevant immune checkpoints. We find that tumor-infiltrating lymphocytes and TCR counts vary widely across cancer types and within each diagnosis, and notably are significantly predictive of survival in osteosarcoma patients. We identify potential cancer-specific immunotherapeutic targets for adoptive cell therapies including cell-surface proteins, tumor germline antigens, and lineage-specific transcription factors. Using an orthogonal immunopeptidomics approach, we find several potential immunotherapeutic targets in osteosarcoma and Ewing sarcoma and validated PRAME as a bona fide multi-pediatric cancer target. Importantly, this work provides a critical framework for immune targeting of extracranial solid tumors using parallel immuno-transcriptomic and -peptidomic approaches.

Genomic and Transcriptomic Analysis of Relapsed and Refractory Childhood Solid Tumors Reveals a Diverse Molecular Landscape and Mechanisms of Immune Evasion.

Children with treatment-refractory or relapsed (R/R) tumors face poor prognoses. As the genomic underpinnings driving R/R disease are not well defined, we describe here the genomic and transcriptomic landscapes of R/R solid tumors from 202 patients enrolled in Beat Childhood Cancer Consortium clinical trials. Tumor mutational burden (TMB) was elevated relative to untreated tumors at diagnosis, with one-third of tumors classified as having a pediatric high TMB. Prior chemotherapy exposure influenced the mutational landscape of these R/R tumors, with more than 40% of tumors demonstrating mutational signatures associated with platinum or temozolomide chemotherapy and two tumors showing treatment-associated hypermutation. Immunogenomic profiling found a heterogenous pattern of neoantigen and MHC class I expression and a general absence of immune infiltration. Transcriptional analysis and functional gene set enrichment analysis identified cross-pathology clusters associated with development, immune signaling, and cellular signaling pathways. While the landscapes of these R/R tumors reflected those of their corresponding untreated tumors at diagnosis, important exceptions were observed, suggestive of tumor evolution, treatment resistance mechanisms, and mutagenic etiologies of treatment. SIGNIFICANCE: Tumor heterogeneity, chemotherapy exposure, and tumor evolution contribute to the molecular profiles and increased mutational burden that occur in treatment-refractory and relapsed childhood solid tumors.

Genomic Classification and Clinical Outcome in Rhabdomyosarcoma: A Report From an International Consortium.

PURPOSE: Rhabdomyosarcoma is the most common soft tissue sarcoma of childhood. Despite aggressive therapy, the 5-year survival rate for patients with metastatic or recurrent disease remains poor, and beyond PAX-FOXO1 fusion status, no genomic markers are available for risk stratification. We present an international consortium study designed to determine the incidence of driver mutations and their association with clinical outcome. PATIENTS AND METHODS: Tumor samples collected from patients enrolled on Children's Oncology Group trials (1998-2017) and UK patients enrolled on malignant mesenchymal tumor and RMS2005 (1995-2016) trials were subjected to custom-capture sequencing. Mutations, indels, gene deletions, and amplifications were identified, and survival analysis was performed. RESULTS: DNA from 641 patients was suitable for analyses. A median of one mutation was found per tumor. In FOXO1 fusion-negative cases, mutation of any RAS pathway member was found in > 50% of cases, and 21% had no putative driver mutation identified. BCOR (15%), NF1 (15%), and TP53 (13%) mutations were found at a higher incidence than previously reported and TP53 mutations were associated with worse outcomes in both fusion-negative and FOXO1 fusion-positive cases. Interestingly, mutations in RAS isoforms predominated in infants < 1 year (64% of cases). Mutation of MYOD1 was associated with histologic patterns beyond those previously described, older age, head and neck primary site, and a dismal survival. Finally, we provide a searchable companion database (ClinOmics), containing all genomic variants, and clinical annotation including survival data. CONCLUSION: This is the largest genomic characterization of clinically annotated rhabdomyosarcoma tumors to date and provides prognostic genetic features that refine risk stratification and will be incorporated into prospective trials.

Interaction between SNAI2 and MYOD enhances oncogenesis and suppresses differentiation in Fusion Negative Rhabdomyosarcoma.

Rhabdomyosarcoma (RMS) is an aggressive pediatric malignancy of the muscle, that includes Fusion Positive (FP)-RMS harboring PAX3/7-FOXO1 and Fusion Negative (FN)-RMS commonly with RAS pathway mutations. RMS express myogenic master transcription factors MYOD and MYOG yet are unable to terminally differentiate. Here, we report that SNAI2 is highly expressed in FN-RMS, is oncogenic, blocks myogenic differentiation, and promotes growth. MYOD activates SNAI2 transcription via super enhancers with striped 3D contact architecture. Genome wide chromatin binding analysis demonstrates that SNAI2 preferentially binds enhancer elements and competes with MYOD at a subset of myogenic enhancers required for terminal differentiation. SNAI2 also suppresses expression of a muscle differentiation program modulated by MYOG, MEF2, and CDKN1A. Further, RAS/MEK-signaling modulates SNAI2 levels and binding to chromatin, suggesting that the differentiation blockade by oncogenic RAS is mediated in part by SNAI2. Thus, an interplay between SNAI2, MYOD, and RAS prevents myogenic differentiation and promotes tumorigenesis.

Rhabdoid Tumors Are Sensitive to the Protein-Translation Inhibitor Homoharringtonine.

PURPOSE: Rhabdoid tumors are devastating pediatric cancers in need of improved therapies. We sought to identify small molecules that exhibit in vitro and in vivo efficacy against preclinical models of rhabdoid tumor. EXPERIMENTAL DESIGN: We screened eight rhabdoid tumor cell lines with 481 small molecules and compared their sensitivity with that of 879 other cancer cell lines. Genome-scale CRISPR-Cas9 inactivation screens in rhabdoid tumors were analyzed to confirm target vulnerabilities. Gene expression and CRISPR-Cas9 data were queried across cell lines and primary rhabdoid tumors to discover biomarkers of small-molecule sensitivity. Molecular correlates were validated by manipulating gene expression. Subcutaneous rhabdoid tumor xenografts were treated with the most effective drug to confirm in vitro results. RESULTS: Small-molecule screening identified the protein-translation inhibitor homoharringtonine (HHT), an FDA-approved treatment for chronic myelogenous leukemia (CML), as the sole drug to which all rhabdoid tumor cell lines were selectively sensitive. Validation studies confirmed the sensitivity of rhabdoid tumor to HHT was comparable with that of CML cell lines. Low expression of the antiapoptotic gene BCL2L1, which encodes Bcl-XL, was the strongest predictor of HHT sensitivity, and HHT treatment consistently depleted Mcl-1, the synthetic-lethal antiapoptotic partner of Bcl-XL. Rhabdoid tumor cell lines and primary-tumor samples expressed low BCL2L1, and overexpression of BCL2L1 induced resistance to HHT in rhabdoid tumor cells. Furthermore, HHT treatment inhibited rhabdoid tumor cell line and patient-derived xenograft growth in vivo. CONCLUSIONS: Rhabdoid tumor cell lines and xenografts are highly sensitive to HHT, at least partially due to their low expression of BCL2L1. HHT may have therapeutic potential against rhabdoid tumors.

Histone hyperacetylation disrupts core gene regulatory architecture in rhabdomyosarcoma.

Core regulatory transcription factors (CR TFs) orchestrate the placement of super-enhancers (SEs) to activate transcription of cell-identity specifying gene networks, and are critical in promoting cancer. Here, we define the core regulatory circuitry of rhabdomyosarcoma and identify critical CR TF dependencies. These CR TFs build SEs that have the highest levels of histone acetylation, yet paradoxically the same SEs also harbor the greatest amounts of histone deacetylases. We find that hyperacetylation selectively halts CR TF transcription. To investigate the architectural determinants of this phenotype, we used absolute quantification of architecture (AQuA) HiChIP, which revealed erosion of native SE contacts, and aberrant spreading of contacts that involved histone acetylation. Hyperacetylation removes RNA polymerase II (RNA Pol II) from core regulatory genetic elements, and eliminates RNA Pol II but not BRD4 phase condensates. This study identifies an SE-specific requirement for balancing histone modification states to maintain SE architecture and CR TF transcription.

SRAssembler: Selective Recursive local Assembly of homologous genomic regions.

BACKGROUND: The falling cost of next-generation sequencing technology has allowed deep sequencing across related species and of individuals within species. Whole genome assemblies from these data remain high time- and resource-consuming computational tasks, particularly if best solutions are sought using different assembly strategies and parameter sets. However, in many cases, the underlying research questions are not genome-wide but rather target specific genes or sets of genes. We describe a novel assembly tool, SRAssembler, that efficiently assembles only contigs containing potential homologs of a gene or protein query, thus enabling gene-specific genome studies over large numbers of short read samples. RESULTS: We demonstrate the functionality of SRAssembler with examples largely drawn from plant genomics. The workflow implements a recursive strategy by which relevant reads are successively pulled from the input sets based on overlapping significant matches, resulting in virtual chromosome walking. The typical workflow behavior is illustrated with assembly of simulated reads. Applications to real data show that SRAssembler produces homologous contigs of equivalent quality to whole genome assemblies. Settings can be chosen to not only assemble presumed orthologs but also paralogous gene loci in distinct contigs. A key application is assembly of the same locus in many individuals from population genome data, which provides assessment of structural variation beyond what can be inferred from read mapping to a reference genome alone. SRAssembler can be used on modest computing resources or used in parallel on high performance computing clusters (most easily by invoking a dedicated Singularity image). CONCLUSIONS: SRAssembler offers an efficient tool to complement whole genome assembly software. It can be used to solve gene-specific research questions based on large genomic read samples from multiple sources and would be an expedient choice when whole genome assembly from the reads is either not feasible, too costly, or unnecessary. The program can also aid decision making on the depth of sequencing in an ongoing novel genome sequencing project or with respect to ultimate whole genome assembly strategies.