Chromatin interactions connect distal regulatory elements to target gene promoters guiding stimulus- and lineage-specific transcription. Few factors securing chromatin interactions have so far been identified. Here, by integrating chromatin interaction maps with the large collection of transcription factor-binding profiles provided by the ENCODE project, we demonstrate that the zinc-finger protein ZNF143 preferentially occupies anchors of chromatin interactions connecting promoters with distal regulatory elements. It binds directly to promoters and associates with lineage-specific chromatin interactions and gene expression. Silencing ZNF143 or modulating its DNA-binding affinity using single-nucleotide polymorphisms (SNPs) as a surrogate of site-directed mutagenesis reveals the sequence dependency of chromatin interactions at gene promoters. We also find that chromatin interactions alone do not regulate gene expression. Together, our results identify ZNF143 as a novel chromatin-looping factor that contributes to the architectural foundation of the genome by providing sequence specificity at promoters connected with distal regulatory elements.

Patient-derived tumor organoids (PDOs) are a highly promising preclinical model that recapitulates the histology, gene expression, and drug response of the donor patient tumor. Currently, PDO culture relies on basement-membrane extract (BME), which suffers from batch-to-batch variability, the presence of xenogeneic compounds and residual growth factors, and poor control of mechanical properties. Additionally, for the development of new organoid lines from patient-derived xenografts, contamination of murine host cells poses a problem. We propose a nanofibrillar hydrogel (EKGel) for the initiation and growth of breast cancer PDOs. PDOs grown in EKGel have histopathologic features, gene expression, and drug response that are similar to those of their parental tumors and PDOs in BME. In addition, EKGel offers reduced batch-to-batch variability, a range of mechanical properties, and suppressed contamination from murine cells. These results show that EKGel is an improved alternative to BME matrices for the initiation, growth, and maintenance of breast cancer PDOs.

Serial circulating tumor DNA (ctDNA) monitoring is emerging as a non-invasive strategy to predict and monitor immune checkpoint blockade (ICB) therapeutic efficacy across cancer types. Yet, limited data exist to show the relationship between ctDNA dynamics and tumor genome and immune microenvironment in patients receiving ICB. Here, we present an in-depth analysis of clinical, whole-exome, transcriptome, and ctDNA profiles of 73 patients with advanced solid tumors, across 30 cancer types, from a phase II basket clinical trial of pembrolizumab (NCT02644369) and report changes in genomic and immune landscapes (primary outcomes). Patients stratified by ctDNA and tumor burden dynamics correspond with survival and clinical benefit. High mutation burden, high expression of immune signatures, and mutations in BRCA2 are associated with pembrolizumab molecular sensitivity, while abundant copy-number alterations and B2M loss-of-heterozygosity corresponded with resistance. Upon treatment, induction of genes expressed by T cell, B cell, and myeloid cell populations are consistent with sensitivity and resistance. We identified the upregulated expression of PLA2G2D, an immune-regulating phospholipase, as a potential biomarker of adaptive resistance to ICB. Together, these findings provide insights into the diversity of immunogenomic mechanisms that underpin pembrolizumab outcomes.

Statins, a family of FDA-approved cholesterol-lowering drugs that inhibit the rate-limiting enzyme of the mevalonate metabolic pathway, have demonstrated anticancer activity. Evidence shows that dipyridamole potentiates statin-induced cancer cell death by blocking a restorative feedback loop triggered by statin treatment. Leveraging this knowledge, we develop an integrative pharmacogenomics pipeline to identify compounds similar to dipyridamole at the level of drug structure, cell sensitivity and molecular perturbation. To overcome the complex polypharmacology of dipyridamole, we focus our pharmacogenomics pipeline on mevalonate pathway genes, which we name mevalonate drug-network fusion (MVA-DNF). We validate top-ranked compounds, nelfinavir and honokiol, and identify that low expression of the canonical epithelial cell marker, E-cadherin, is associated with statin-compound synergy. Analysis of remaining prioritized hits led to the validation of additional compounds, clotrimazole and vemurafenib. Thus, our computational pharmacogenomic approach identifies actionable compounds with pathway-specific activities.

Next-generation sequencing technologies have recently been used in pharmacogenomic studies to characterize large panels of cancer cell lines at the genomic and transcriptomic levels. Among these technologies, RNA-sequencing enable profiling of alternatively spliced transcripts. Given the high frequency of mRNA splicing in cancers, linking this feature to drug response will open new avenues of research in biomarker discovery. To identify robust transcriptomic biomarkers for drug response across studies, we develop a meta-analytical framework combining the pharmacological data from two large-scale drug screening datasets. We use an independent pan-cancer pharmacogenomic dataset to test the robustness of our candidate biomarkers across multiple cancer types. We further analyze two independent breast cancer datasets and find that specific isoforms of IGF2BP2, NECTIN4, ITGB6, and KLHDC9 are significantly associated with AZD6244, lapatinib, erlotinib, and paclitaxel, respectively. Our results support isoform expressions as a rich resource for biomarkers predictive of drug response.

Glioblastoma (GBM) is a deadly cancer in which cancer stem cells (CSCs) sustain tumor growth and contribute to therapeutic resistance. Protein arginine methyltransferase 5 (PRMT5) has recently emerged as a promising target in GBM. Using two orthogonal-acting inhibitors of PRMT5 (GSK591 or LLY-283), we show that pharmacological inhibition of PRMT5 suppresses the growth of a cohort of 46 patient-derived GBM stem cell cultures, with the proneural subtype showing greater sensitivity. We show that PRMT5 inhibition causes widespread disruption of splicing across the transcriptome, particularly affecting cell cycle gene products. We identify a GBM splicing signature that correlates with the degree of response to PRMT5 inhibition. Importantly, we demonstrate that LLY-283 is brain-penetrant and significantly prolongs the survival of mice with orthotopic patient-derived xenografts. Collectively, our findings provide a rationale for the clinical development of brain penetrant PRMT5 inhibitors as treatment for GBM.

In acute myeloid leukemia (AML), molecular heterogeneity across patients constitutes a major challenge for prognosis and therapy. AML with NPM1 mutation is a distinct genetic entity in the revised World Health Organization classification. However, differing patterns of co-mutation and response to therapy within this group necessitate further stratification. Here we report two distinct subtypes within NPM1 mutated AML patients, which we label as primitive and committed based on the respective presence or absence of a stem cell signature. Using gene expression (RNA-seq), epigenomic (ATAC-seq) and immunophenotyping (CyToF) analysis, we associate each subtype with specific molecular characteristics, disease differentiation state and patient survival. Using ex vivo drug sensitivity profiling, we show a differential drug response of the subtypes to specific kinase inhibitors, irrespective of the FLT3-ITD status. Differential drug responses of the primitive and committed subtype are validated in an independent AML cohort. Our results highlight heterogeneity among NPM1 mutated AML patient samples based on stemness and suggest that the addition of kinase inhibitors to the treatment of cases with the primitive signature, lacking FLT3-ITD, could have therapeutic benefit.

Pineoblastoma is a rare pediatric cancer induced by germline mutations in the tumor suppressors RB1 or DICER1. Presence of leptomeningeal metastases is indicative of poor prognosis. Here we report that inactivation of Rb plus p53 via a WAP-Cre transgene, commonly used to target the mammary gland during pregnancy, induces metastatic pineoblastoma resembling the human disease with 100% penetrance. A stabilizing mutation rather than deletion of p53 accelerates metastatic dissemination. Deletion of Dicer1 plus p53 via WAP-Cre also predisposes to pineoblastoma, albeit with lower penetrance. In silico analysis predicts tricyclic antidepressants such as nortriptyline as potential therapeutics for both pineoblastoma models. Nortriptyline disrupts the lysosome, leading to accumulation of non-functional autophagosome, cathepsin B release and pineoblastoma cell death. Nortriptyline further synergizes with the antineoplastic drug gemcitabine to effectively suppress pineoblastoma in our preclinical models, offering new modality for this lethal childhood malignancy.

The human genome is partitioned into a collection of genomic features, inclusive of genes, transposable elements, lamina interacting regions, early replicating control elements and cis-regulatory elements, such as promoters, enhancers, and anchors of chromatin interactions. Uneven distribution of these features within chromosomes gives rise to clusters, such as topologically associating domains (TADs), lamina-associated domains, clusters of cis-regulatory elements or large organized chromatin lysine (K) domains (LOCKs). Here we show that LOCKs from diverse histone modifications discriminate primitive from differentiated cell types. Active LOCKs (H3K4me1, H3K4me3 and H3K27ac) cover a higher fraction of the genome in primitive compared to differentiated cell types while repressive LOCKs (H3K9me3, H3K27me3 and H3K36me3) do not. Active LOCKs in differentiated cells lie proximal to highly expressed genes while active LOCKs in primitive cells tend to be bivalent. Genes proximal to bivalent LOCKs are minimally expressed in primitive cells. Furthermore, bivalent LOCKs populate TAD boundaries and are preferentially bound by regulators of chromatin interactions, including CTCF, RAD21 and ZNF143. Together, our results argue that LOCKs discriminate primitive from differentiated cell populations.

Human cancers exhibit strong phenotypic differences that can be visualized noninvasively by medical imaging. Radiomics refers to the comprehensive quantification of tumour phenotypes by applying a large number of quantitative image features. Here we present a radiomic analysis of 440 features quantifying tumour image intensity, shape and texture, which are extracted from computed tomography data of 1,019 patients with lung or head-and-neck cancer. We find that a large number of radiomic features have prognostic power in independent data sets of lung and head-and-neck cancer patients, many of which were not identified as significant before. Radiogenomics analysis reveals that a prognostic radiomic signature, capturing intratumour heterogeneity, is associated with underlying gene-expression patterns. These data suggest that radiomics identifies a general prognostic phenotype existing in both lung and head-and-neck cancer. This may have a clinical impact as imaging is routinely used in clinical practice, providing an unprecedented opportunity to improve decision-support in cancer treatment at low cost.

Neuroendocrine prostate cancer (NEPC), a lethal form of the disease, is characterized by loss of androgen receptor (AR) signaling during neuroendocrine transdifferentiation, which results in resistance to AR-targeted therapy. Clinically, genomically and epigenetically, NEPC resembles other types of poorly differentiated neuroendocrine tumors (NETs). Through pan-NET analyses, we identified ONECUT2 as a candidate master transcriptional regulator of poorly differentiated NETs. ONECUT2 ectopic expression in prostate adenocarcinoma synergizes with hypoxia to suppress androgen signaling and induce neuroendocrine plasticity. ONEUCT2 drives tumor aggressiveness in NEPC, partially through regulating hypoxia signaling and tumor hypoxia. Specifically, ONECUT2 activates SMAD3, which regulates hypoxia signaling through modulating HIF1α chromatin-binding, leading NEPC to exhibit higher degrees of hypoxia compared to prostate adenocarcinomas. Treatment with hypoxia-activated prodrug TH-302 potently reduces NEPC tumor growth. Collectively, these results highlight the synergy between ONECUT2 and hypoxia in driving NEPC, and emphasize the potential of hypoxia-directed therapy for NEPC patients.

Reproducibility is essential to open science, as there is limited relevance for findings that can not be reproduced by independent research groups, regardless of its validity. It is therefore crucial for scientists to describe their experiments in sufficient detail so they can be reproduced, scrutinized, challenged, and built upon. However, the intrinsic complexity and continuous growth of biomedical data makes it increasingly difficult to process, analyze, and share with the community in a FAIR (findable, accessible, interoperable, and reusable) manner. To overcome these issues, we created a cloud-based platform called ORCESTRA ( orcestra.ca ), which provides a flexible framework for the reproducible processing of multimodal biomedical data. It enables processing of clinical, genomic and perturbation profiles of cancer samples through automated processing pipelines that are user-customizable. ORCESTRA creates integrated and fully documented data objects with persistent identifiers (DOI) and manages multiple dataset versions, which can be shared for future studies.

Triple negative breast cancer (TNBC) is a deadly form of breast cancer due to the development of resistance to chemotherapy affecting over 30% of patients. New therapeutics and companion biomarkers are urgently needed. Recognizing the elevated expression of glucose transporter 1 (GLUT1, encoded by SLC2A1) and associated metabolic dependencies in TNBC, we investigated the vulnerability of TNBC cell lines and patient-derived samples to GLUT1 inhibition. We report that genetic or pharmacological inhibition of GLUT1 with BAY-876 impairs the growth of a subset of TNBC cells displaying high glycolytic and lower oxidative phosphorylation (OXPHOS) rates. Pathway enrichment analysis of gene expression data suggests that the functionality of the E2F pathway may reflect to some extent OXPHOS activity. Furthermore, the protein levels of retinoblastoma tumor suppressor (RB1) strongly correlate with the degree of sensitivity to GLUT1 inhibition in TNBC, where RB1-negative cells are insensitive to GLUT1 inhibition. Collectively, our results highlight a strong and targetable RB1-GLUT1 metabolic axis in TNBC and warrant clinical evaluation of GLUT1 inhibition in TNBC patients stratified according to RB1 protein expression levels.