Publications

Outcomes for metastatic Ewing sarcoma and osteosarcoma are dismal and have not changed for decades. Oxidative stress attenuates melanoma metastasis, and melanoma cells must reduce oxidative stress to metastasize. We explored this in sarcomas by screening for oxidative stress sensitizers, which identified the class I HDAC inhibitor MS-275 as enhancing vulnerability to reactive oxygen species (ROS) in sarcoma cells. Mechanistically, MS-275 inhibits YB-1 deacetylation, decreasing its binding to 5'-UTRs of NFE2L2 encoding the antioxidant factor NRF2, thereby reducing NFE2L2 translation and synthesis of NRF2 to increase cellular ROS. By global acetylomics, MS-275 promotes rapid acetylation of the YB-1 RNA-binding protein at lysine-81, blocking binding and translational activation of NFE2L2, as well as known YB-1 mRNA targets, HIF1A, and the stress granule nucleator, G3BP1. MS-275 dramatically reduces sarcoma metastasis in vivo, but an MS-275-resistant YB-1K81-to-alanine mutant restores metastatic capacity and NRF2, HIF1α, and G3BP1 synthesis in MS-275-treated mice. These studies describe a novel function for MS-275 through enhanced YB-1 acetylation, thus inhibiting YB-1 translational control of key cytoprotective factors and its pro-metastatic activity.

Inhibition of the androgen receptor (AR) is the mainstay treatment for advanced prostate cancer. Ralaniten (formally EPI-002) prevents AR transcriptional activity by binding to its N-terminal domain (NTD) which is essential for transcriptional activity. Ralaniten acetate (EPI-506) the triacetate pro-drug of ralaniten, remains the only AR-NTD inhibitor to have entered clinical trials (NCT02606123). While well tolerated, the trial was ultimately terminated due to poor pharmacokinetic properties and resulting pill burden. Here we discovered that ralaniten was glucuronidated which resulted in decreased potency. Long-term treatment of prostate cancer cells with ralaniten results in upregulation of UGT2B enzymes with concomitant loss of potency. This has proven to be a useful model with which to facilitate the development of more potent second-generation AR-NTD inhibitors. Glucuronidated metabolites of ralaniten were also detected in the serum of patients in Phase 1 clinical trials. Therefore, we tested an analogue of ralaniten (EPI-045) which was resistant to glucuronidation and demonstrated superiority to ralaniten in our resistant model. These data support that analogues of ralaniten designed to mitigate glucuronidation may optimize clinical responses to AR-NTD inhibitors.

Despite advances in diagnostic tools and therapeutic options, treatment resistance remains a challenge for many cancer patients. Recent studies have found evidence that autophagy, a cellular pathway that delivers cytoplasmic components to lysosomes for degradation and recycling, contributes to treatment resistance in different cancer types. A role for autophagy in resistance to chemotherapies and targeted therapies has been described based largely on associations with various signaling pathways, including MAPK and PI3K/AKT signaling. However, our current understanding of the molecular mechanisms underlying the role of autophagy in facilitating treatment resistance remains limited. Here we provide a comprehensive summary of the evidence linking autophagy to major signaling pathways in the context of treatment resistance and tumor progression, and then highlight recently emerged molecular mechanisms underlying autophagy and the p62/KEAP1/NRF2 and FOXO3A/PUMA axes in chemoresistance.

Extra-cranial malignant rhabdoid tumors (MRTs) and cranial atypical teratoid RTs (ATRTs) are heterogeneous pediatric cancers driven primarily by SMARCB1 loss. To understand the genome-wide molecular relationships between MRTs and ATRTs, we analyze multi-omics data from 140 MRTs and 161 ATRTs. We detect similarities between the MYC subgroup of ATRTs (ATRT-MYC) and extra-cranial MRTs, including global DNA hypomethylation and overexpression of HOX genes and genes involved in mesenchymal development, distinguishing them from other ATRT subgroups that express neural-like features. We identify five DNA methylation subgroups associated with anatomical sites and SMARCB1 mutation patterns. Groups 1, 3, and 4 exhibit cytotoxic T cell infiltration and expression of immune checkpoint regulators, consistent with a potential role for immunotherapy in rhabdoid tumor patients.

In the modern genomics era, genome sequence assemblies are routine practice. However, depending on the methodology, resulting drafts may contain considerable base errors. Although utilities exist for genome base polishing, they work best with high read coverage and do not scale well. We developed ntEdit, a Bloom filter-based genome sequence editing utility that scales to large mammalian and conifer genomes.

Despite being the most common childhood bone tumor, the genomic characterization of osteosarcoma remains incomplete. In particular, very few osteosarcoma metastases have been sequenced to date, critical to better understand mechanisms of progression and evolution in this tumor. We performed an integrated whole genome and exome sequencing analysis of paired primary and metastatic pediatric osteosarcoma specimens to identify recurrent genomic alterations. Sequencing of 13 osteosarcoma patients including 13 primary, 10 metastatic, and 3 locally recurring tumors revealed a highly heterogeneous mutational landscape, including cases of hypermutation and microsatellite instability positivity, but with virtually no recurrent alterations except for mutations involving the tumor suppressor genes RB1 and TP53. At the germline level, we detected alterations in multiple cancer related genes in the majority of the cohort, including those potentially disrupting DNA damage response pathways. Metastases retained only a minimal number of short variants from their corresponding primary tumors, while copy number alterations showed higher conservation. One recurrently amplified gene, KDR, was highly expressed in advanced cases and associated with poor prognosis. © 2019 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Accurate measurement of clonal genotypes, mutational processes, and replication states from individual tumor-cell genomes will facilitate improved understanding of tumor evolution. We have developed DLP+, a scalable single-cell whole-genome sequencing platform implemented using commodity instruments, image-based object recognition, and open source computational methods. Using DLP+, we have generated a resource of 51,926 single-cell genomes and matched cell images from diverse cell types including cell lines, xenografts, and diagnostic samples with limited material. From this resource we have defined variation in mitotic mis-segregation rates across tissue types and genotypes. Analysis of matched genomic and image measurements revealed correlations between cellular morphology and genome ploidy states. Aggregation of cells sharing copy number profiles allowed for calculation of single-nucleotide resolution clonal genotypes and inference of clonal phylogenies and avoided the limitations of bulk deconvolution. Finally, joint analysis over the above features defined clone-specific chromosomal aneuploidy in polyclonal populations.

The ORCA bioinformatics environment is a Docker image that contains hundreds of bioinformatics tools and their dependencies. The ORCA image and accompanying server infrastructure provide a comprehensive bioinformatics environment for education and research. The ORCA environment on a server is implemented using Docker containers, but without requiring users to interact directly with Docker, suitable for novices who may not yet have familiarity with managing containers. ORCA has been used successfully to provide a private bioinformatics environment to external collaborators at a large genome institute, for teaching an undergraduate class on bioinformatics targeted at biologists, and to provide a ready-to-go bioinformatics suite for a hackathon. Using ORCA eliminates time that would be spent debugging software installation issues, so that time may be better spent on education and research.

Cytotoxic CD8{{sup}}+{{/sup}} T cells recognize and eliminate infected or malignant cells that present peptide epitopes derived from intracellularly processed antigens on their surface. However, comprehensive profiling of specific major histocompatibility complex (MHC)-bound peptide epitopes that are naturally processed and capable of eliciting a functional T cell response has been challenging. Here, we report a method for deep and unbiased T cell epitope profiling, using in vitro co-culture of CD8{{sup}}+{{/sup}} T cells together with target cells transduced with high-complexity, epitope-encoding minigene libraries. Target cells that are subject to cytotoxic attack from T cells in co-culture are isolated prior to apoptosis by fluorescence-activated cell sorting, and characterized by sequencing the encoded minigenes. We then validate this highly parallelized method using known murine T cell receptor/peptide-MHC pairs and diverse minigene-encoded epitope libraries. Our data thus suggest that this epitope profiling method allows unambiguous and sensitive identification of naturally processed and MHC-presented peptide epitopes.

High-grade B-cell lymphomas with MYC and BCL2 and/or BCL6 rearrangements (HGBL-DH/THs) include a group of diffuse large B-cell lymphomas (DLBCLs) with inferior outcomes after standard chemoimmunotherapy. We recently described a gene expression signature that identifies 27% of germinal center B-cell DLBCLs (GCB-DLBCLs) as having a double-hit-like expression pattern (DHITsig) and inferior outcomes; however, only half of these cases have both MYC and BCL2 translocations identifiable using standard breakapart fluorescence in situ hybridization (FISH). Here, 20 DHITsig+ GCB-DLBCLs apparently lacking MYC and/or BCL2 rearrangements underwent whole-genome sequencing. This revealed 6 tumors with MYC or BCL2 rearrangements that were cryptic to breakapart FISH. Copy-number analysis identified 3 tumors with MYC and 6 tumors with MIR17HG gains or amplifications, both of which may contribute to dysregulation of MYC and its downstream pathways. Focal deletions of the PVT1 promoter were observed exclusively among DHITsig+ tumors lacking MYC translocations; this may also contribute to MYC overexpression. These results highlight that FISH fails to identify all HGBL-DH/THs, while revealing a range of other genetic mechanisms potentially underlying MYC dysregulation in DHITsig+ DLBCL, suggesting that gene expression profiling is more sensitive for identifying the biology underlying poor outcomes in GCB-DLBCL.