Fine-mapping causal tissues and genes at disease-associated loci

Heritable diseases often manifest in a highly tissue-specific manner, with different disease loci mediated by genes in distinct tissues or cell types. We propose Tissue-Gene Fine-Mapping (TGFM), a fine-mapping method that infers the posterior probability (PIP) for each gene-tissue pair to mediate a disease locus by analyzing GWAS summary statistics (and in-sample LD) and leveraging eQTL data from diverse tissues to build cis-predicted expression models; TGFM also assigns PIPs to causal variants that are not mediated by gene expression in assayed genes and tissues. TGFM accounts for both co-regulation across genes and tissues and LD between SNPs (generalizing existing fine-mapping methods), and incorporates genome-wide estimates of each tissue's contribution to disease as tissue-level priors. TGFM was well-calibrated and moderately well-powered in simulations; unlike previous methods, TGFM was able to attain correct calibration by modeling uncertainty in cis-predicted expression models. We applied TGFM to 45 UK Biobank diseases/traits (average N = 316 K ) using eQTL data from 38 GTEx tissues. TGFM identified an average of 147 PIP > 0.5 causal genetic elements per disease/trait, of which 11% were gene-tissue pairs. Implicated gene-tissue pairs were concentrated in known disease-critical tissues, and causal genes were strongly enriched in disease-relevant gene sets. Causal gene-tissue pairs identified by TGFM recapitulated known biology (e.g., TPO-thyroid for Hypothyroidism), but also included biologically plausible novel findings (e.g., SLC20A2-artery aorta for Diastolic blood pressure). Further application of TGFM to single-cell eQTL data from 9 cell types in peripheral blood mononuclear cells (PBMC), analyzed jointly with GTEx tissues, identified 30 additional causal gene-PBMC cell type pairs at PIP > 0.5-primarily for autoimmune disease and blood cell traits, including the well-established role of CTLA4 in CD8+ T cells for All autoimmune disease. In conclusion, TGFM is a robust and powerful method for fine-mapping causal tissues and genes at disease-associated loci.

Copper induces cell death by targeting lipoylated TCA cycle proteins

Copper is an essential cofactor for all organisms, and yet it becomes toxic if concentrations exceed a threshold maintained by evolutionarily conserved homeostatic mechanisms. How excess copper induces cell death, however, is unknown. Here, we show in human cells that copper-dependent, regulated cell death is distinct from known death mechanisms and is dependent on mitochondrial respiration. We show that copper-dependent death occurs by means of direct binding of copper to lipoylated components of the tricarboxylic acid (TCA) cycle. This results in lipoylated protein aggregation and subsequent iron-sulfur cluster protein loss, which leads to proteotoxic stress and ultimately cell death. These findings may explain the need for ancient copper homeostatic mechanisms.

Abstract 3400: Adenosine receptor antagonists exhibit potent and selective off-target killing of FOXA1-high cancers

Drugs targeting adenosine receptors were originally developed for the treatment of Parkinson's disease and are now being tested in immuno-oncology clinical trials in combination with checkpoint inhibitors. We recently reported the killing activity of 4,518 drugs against 578 diverse cancer cell lines determined using the PRISM molecular barcoding approach. Surprisingly, three established adenosine receptor antagonists (CGS-15943, MRS-1220, and SCH-58261) showed potent and selective killing of FOXA1-high cancer cell lines without the need for immune cells. FOXA1 is a lineage-restricted transcription factor in luminal breast cancer, hepatocellular carcinoma, and prostate cancer without known small molecule inhibitors. We find that cytotoxic activity is limited to adenosine antagonists with a three-member aromatic core bound to a furan group, thus indicating a potential off-target mechanism of action. To identify genomic modulators of drug response, we performed genome-wide CRISPR/Cas9 knockout modifier screens. Killing by CGS-15943 and MRS-1220 was rescued by knockout of the aryl hydrocarbon receptor (AHR) and its nuclear partner ARNT. In confirmatory studies, knockout of AHR completely rescued killing by CGS-15943 in multiple cell types. Co-treatment with an AHR small molecule antagonist also rescued cell viability. Knockout of adenosine receptors did not alter drug response. Given that AHR is a known transcriptional regulator, we performed global mRNA sequencing to assess transcriptional changes induced by CGS-15943. The top two genes induced were the p450 enzymes CYP1A1 and CYP1B1. To determine sufficiency, we overexpressed CYP1A1 in a resistant cell line. Ectopic CYP1A1 expression sensitized to CGS-15943-mediated killing. Mass spectrometry revealed covalent trapping of a reactive metabolite by glutathione and potassium cyanide following in vitro incubation with liver microsomes. In addition, treatment of breast cancer cells with CGS-15943 for 24 hours resulted in increased γ-H2AX phosphorylation by western blot, indicative of DNA double stranded breaks. In summary, we identified off-target anti-cancer activity of multiple established adenosine receptor antagonists mediated by activation of AHR. Future studies will evaluate the functional contribution of FOXA1 and activity in vivo. Starting from a phenotypic screening hit, we leverage functional genomics to unlock the underlying mechanism of action. This project will pave the way for developing more effective therapies for biomarker-selected cancers, with potential to improve the care of patients with liver, breast, and prostate cancer.

Citation Format: Steven M. Corsello, Ryan D. Spangler, Ranad Humeidi, Caitlin N. Harrington, Rohith T. Nagari, Ritu Singh, Vickie Wang, Mustafa Kocak, Jordan Rossen, Amael Madec, Nancy Dumont, Todd R. Golub. Adenosine receptor antagonists exhibit potent and selective off-target killing of FOXA1-high cancers [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 3400.

Discovering the anti-cancer potential of non-oncology drugs by systematic viability profiling

Anti-cancer uses of non-oncology drugs have occasionally been found, but such discoveries have been serendipitous. We sought to create a public resource containing the growth inhibitory activity of 4,518 drugs tested across 578 human cancer cell lines. We used PRISM, a molecular barcoding method, to screen drugs against cell lines in pools. An unexpectedly large number of non-oncology drugs selectively inhibited subsets of cancer cell lines in a manner predictable from the cell lines' molecular features. Our findings include compounds that killed by inducing PDE3A-SLFN12 complex formation; vanadium-containing compounds whose killing depended on the sulfate transporter SLC26A2; the alcohol dependence drug disulfiram, which killed cells with low expression of metallothioneins; and the anti-inflammatory drug tepoxalin, which killed via the multi-drug resistance protein ABCB1. The PRISM drug repurposing resource (https://depmap.org/repurposing) is a starting point to develop new oncology therapeutics, and more rarely, for potential direct clinical translation.

Discovering the anti-cancer potential of non-oncology drugs by systematic viability profiling

Anti-cancer uses of non-oncology drugs have occasionally been found, but such discoveries have been serendipitous. We sought to create a public resource containing the growth inhibitory activity of 4,518 drugs tested across 578 human cancer cell lines. We used PRISM, a molecular barcoding method, to screen drugs against cell lines in pools. An unexpectedly large number of non-oncology drugs selectively inhibited subsets of cancer cell lines in a manner predictable from the cell lines' molecular features. Our findings include compounds that killed by inducing PDE3A-SLFN12 complex formation; vanadium-containing compounds whose killing depended on the sulfate transporter SLC26A2; the alcohol dependence drug disulfiram, which killed cells with low expression of metallothioneins; and the anti-inflammatory drug tepoxalin, which killed via the multi-drug resistance protein ABCB1. The PRISM drug repurposing resource (https://depmap.org/repurposing) is a starting point to develop new oncology therapeutics, and more rarely, for potential direct clinical translation.

Abstract B121: An oncogene-linked prodrug strategy in lung cancer

KEAP1-mutant non-small cell lung cancer is a high prevalence indication that responds poorly to conventional chemotherapy owing to constitutive activation of NRF2 and its associated drug metabolism target genes. Using an approach to identify compounds that selectively kill cancer cell lines in a biomarker-driven manner, we identify BRD-K19050021 (K1905), a compound displaying strong toxicity in cells expressing CYP4F11 - a cytochrome p450 family member and NRF2 transcriptional target. Using genome-wide pooled CRISPR screens, we show that CYP4F11 activity is necessary and sufficient for response to K1905, and in acquired resistance models developed from multiple cell lines, we show a minor subpopulation of cells mechanistically converge on suppressed CYP4F11 expression to bypass K1905 response. CYP4F11 converts K1905 from a prodrug into a covalent active metabolite that alkylates several cellular targets, triggering all three canonical arms of the unfolded protein response pathway and culminating in cell death. We propose that CYP4F11 and similar metabolic enzyme activities promoted by oncogenic drivers represent a unique opportunity to restrict prodrug activation within tumors, provided that normal tissue expression of such prodrug-converting enzymes do not diminish therapeutic index.

Citation Format: Aaron Boudreau, Jennifer Roth, James Rong, Niclas Olsson, Chris Mader, Jordan Byran, Jordan Rossen, Li Wang, Kevin Larpenteur, Amy Goodale, Colin Trepicchio, Samantha Bender, Aviad Tsherniak, Aravind Subramanian, Mustafa Kocak, Federica Piccioni, Josh Bittker, Cong Zhu, Frank Li, Nick Eriksson, Daphne Koller, Fiona McAllister, Todd Golub, Jeff Settleman, Ari Firestone, David Stokoe. An oncogene-linked prodrug strategy in lung cancer [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference on Molecular Targets and Cancer Therapeutics; 2019 Oct 26-30; Boston, MA. Philadelphia (PA): AACR; Mol Cancer Ther 2019;18(12 Suppl):Abstract nr B121. doi:10.1158/1535-7163.TARG-19-B121

Abstract 2690: Massively parallel multiplexed methods to screen hundreds of barcoded cancer cell line models with small molecules or genetic perturbations using next-generation sequencing

Phenotypic screening is a valuable tool to identify compounds to treat cancer, but is limited as it is time and resource intensive to screen hundreds of cancer cell lines. In order to increase the throughput of phenotypic screening, we set out to expand and further develop the PRISM method (Profiling Relative Inhibition Simultaneously in Mixtures) [Yu et al., 2016]. PRISM is a method to barcode cell lines with a unique 24 nucleotide barcode and mix them together to screen simultaneously in a pool, decreasing the time and cost of screening. The previously described PRISM method was limited in that it used only 100 adherent cell lines of one cancer model context, it required the use of a highly specialized Luminex bead-based detection system, and it was only applicable to small molecule screening.

Here we report on the expansion of our barcoded cell line collection to 800 cancer cell lines of over 45 lineages, an improved method called PRISMseq using next-generation sequencing, and a novel method to perform genetic perturbations in 500 cancer cell lines simultaneously called PRISPR (PRISM/CRISPR). PRISMseq allows for assaying compound cell line sensitivity profiles in a large pool of hundreds of cell lines and detecting the DNA barcodes using next-generation sequencing (NGS). We recapitulate the known biology for established oncology drugs like the BRAF, EGFR, BCR-ABL and MDM2 inhibitor classes using this method with our cell line panel. We also further extended this method to be applicable to genetic perturbation using CRISPR/Cas9 knockout of individual genes concurrently in hundreds of cancer cell lines. We are able to recover expected biomarkers or genetic dependencies for our CRISPR/Cas9 validation reagents.

This expanded PRISM profiling approach increases statistical power due to the addition of cancer contexts in our cell line collection, and improves versatility by enabling the screening of both small molecules and genetic perturbations. We believe that this method has overcome the limitations of the original PRISM method. It has also overcome the limitations of phenotypic screening, as the resources required to screen hundreds of cell lines has been decreased by orders of magnitude.

Citation Format: Samantha Bender, Cong Zhu, Li Wang, Michael Rothberg, Joshua Dempster, Brenton Paolella, Mustafa Kocak, Massami Laird, Jordan Rossen, Karolina Stumbraite, Jordan Bryan, Vickie Wang, John Doench, Francisca Vazquez, Aviad Tsherniak, Todd Golub, Jennifer Roth. Massively parallel multiplexed methods to screen hundreds of barcoded cancer cell line models with small molecules or genetic perturbations using next-generation sequencing [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 2690.

Abstract 2948: Novel cell line barcoding method reveals tepoxalin as a selective drug against MDR1-high tumor cells

Drug repurposing for cancer therapy has seen notable success, such as the reintroduction of thalidomide. To discover new indications, all existing drugs known to be safe should be systematically evaluated for anti-cancer properties. However, given cost and throughput considerations, phenotypic screens are typically limited along the cell line (e.g., NCI-60 cell lines) or compound dimension (e.g., testing only several hundred cancer drugs). To identify new repurposing opportunities at scale, we tested more than 4,600 existing drugs for cytotoxicity against 578 diverse cancer cell lines using PRISM, a recently developed multiplexed cell viability assay. High gene expression of the ABCB1 transporter (MDR1/p-glycoprotein) was strongly correlated with resistance to numerous approved cancer therapeutics including taxanes, anthracyclines, vinca alkaloids, and proteasome inhibitors. Surprisingly, we identified a single drug, tepoxalin, with the opposite profile: selective activity against MDR1-high cancer cell lines.

Tepoxalin, an oral cyclooxygenase inhibitor, is FDA-approved for treatment of osteoarthritis in dogs and was previously found to be safe in human trials. To evaluate for alternative mechanism of action and genomic mediators of tepoxalin resistance, we conducted pooled genome-wide CRISPR-Cas9 modifier screens in the LS1034 colorectal cancer cell line. The CRISPR knockout screen revealed ABCB1 as the top enriched gene mediating tepoxalin resistance, while the CRISPR activation screen revealed ABCB1 as the top depleted gene. Overexpression of ABCB1 in ovarian cancer cell lines also sensitized to tepoxalin killing. Small molecule combination testing demonstrated rescue of tepoxalin killing by known ABCB1 inhibitors. Future studies will include biophysical characterization of the interaction between tepoxalin and the ABCB1 protein as well as in vivo efficacy testing. Our results demonstrate that tepoxalin or related derivatives may be useful in the treatment of chemotherapy-resistant colorectal cancer, ovarian cancer, lymphoma, and other malignancies.

Citation Format: Steven M. Corsello, Ryan D. Spangler, Rohith T. Nagari, Mustafa Kocak, Jordan Rossen, Patrick O'Hearn, Jennifer Roth, Alfredo Gonzalez, Nancy Dumont, John Doench, Jesse S. Boehm, Francisca Vazquez, Aviad Tsherniak, Todd R. Golub. Novel cell line barcoding method reveals tepoxalin as a selective drug against MDR1-high tumor cells [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 2948.

Author Correction: Mitochondrial metabolism promotes adaptation to proteotoxic stress

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Mitochondrial metabolism promotes adaptation to proteotoxic stress

The mechanisms by which cells adapt to proteotoxic stress are largely unknown, but key to understanding how tumor cells, particularly in vivo, are largely resistant to proteasome inhibitors. Analysis of cancer cell lines, mouse xenografts and patient-derived tumor samples all showed an association between mitochondrial metabolism and proteasome inhibitor sensitivity. When cells were forced to use oxidative phosphorylation rather than glycolysis, they became proteasome inhibitor-resistant. This mitochondrial state, however, creates a unique vulnerability: sensitivity to the small-molecule compound elesclomol. Genome-wide CRISPR/Cas9 screening showed that a single gene, encoding the mitochondrial reductase FDX1, could rescue elesclomol-induced cell death. Enzymatic function and NMR-based analyses further showed that FDX1 is the direct target of elesclomol, which promotes a unique form of copper-dependent cell death. These studies elucidate a fundamental mechanism by which cells adapt to proteotoxic stress and suggests strategies to mitigate proteasome inhibitor-resistance.

Reduction of the nosocomial meticillin-resistant Staphylococcus aureus incidence density by a region-wide search and follow-strategy in forty German hospitals of the EUREGIO, 2009 to 2011

Meticillin-resistant Staphylococcus aureus (MRSA) disseminates between hospitals serving one patient catchment area. Successful prevention and control requires concerted efforts and regional surveillance. Forty hospitals located in the German EUREGIO have established a network for combating MRSA. In 2007 they agreed upon a synchronised strategy for screening of risk patients and a standard for transmissionbased precautions (search and follow). The same year, the hospitals started synchronised MRSA prevention and annually reporting MRSA-data to the public health authorities. The median rate of screening cultures per 100 patients admitted increased from 4.38 in 2007 to 34.4 in 2011 (p<0.0001). Between 2007 and 2011, the overall incidence density of MRSA (0.87 MRSA cases/1,000 patient days vs 1.54; p<0.0001) increased significantly. In contrast, both the incidence density of nosocomial MRSA cases (0.13 nosocomial MRSA cases/1,000 patient days in 2009 vs 0.08 in 2011; p=0.0084) and the MRSA-days-associated nosocomial MRSA rate (5.51 nosocomial MRSA cases/1,000 MRSA days in 2009 vs 3.80 in 2011; p=0.0437) decreased significantly after the second year of the project. We documented adherence to the regional screening strategy resulting in improved detection of MRSA carriers at admission. Subsequently, after two years the nosocomial MRSA-incidence density was reduced. Regional surveillance data, annually provided as benchmarking to the regional hospitals and public health authorities, indicated successful prevention.

Monte Verde: seaweed, food, medicine, and the peopling of South America

The identification of human artifacts at the early archaeological site of Monte Verde in southern Chile has raised questions of when and how people reached the tip of South America without leaving much other evidence in the New World. Remains of nine species of marine algae were recovered from hearths and other features at Monte Verde II, an upper occupational layer, and were directly dated between 14,220 and 13,980 calendar years before the present ( approximately 12,310 and 12,290 carbon-14 years ago). These findings support the archaeological interpretation of the site and indicate that the site's inhabitants used seaweed from distant beaches and estuarine environments for food and medicine. These data are consistent with the ideas that an early settlement of South America was along the Pacific coast and that seaweeds were important to the diet and health of early humans in the Americas.

Interaction of mouse hepatitis virus (MHV) spike glycoprotein with receptor glycoprotein MHVR is required for infection with an MHV strain that expresses the hemagglutinin-esterase glycoprotein

In addition to the spike (S) glycoprotein that binds to carcinoembryonic antigen-related receptors on the host cell membrane, some strains of mouse coronavirus (mouse hepatitis virus [MHV]) express a hemagglutinin esterase (HE) glycoprotein with hemagglutinating and acetylesterase activity. Virions of strains that do not express HE, such as MHV-A59, can infect mouse fibroblasts in vitro, showing that the HE glycoprotein is not required for infection of these cells. The present work was done to study whether interaction of the HE glycoprotein with carbohydrate moieties could lead to virus entry and infection in the absence of interaction of the S glycoprotein with its receptor glycoprotein, MHVR. The DVIM strain of MHV expresses large amounts of HE glycoprotein, as shown by hemadsorption, acetylesterase activity, and immunoreactivity with antibodies directed against the HE glycoprotein of bovine coronavirus. A monoclonal anti-MHVR antibody, MAb-CC1, blocks binding of virus S glycoprotein to MHVR and blocks infection of MHV strains that do not express HE. MAb-CC1 also prevented MHV-DVIM infection of mouse DBT cells and primary mouse glial cell cultures. Although MDCK-I cells express O-acetylated sialic acid residues on their plasma membranes, these canine cells were resistant to infection with MHV-A59 and MHV-DVIM. Transfection of MDCK-I cells with MHVR cDNA made them susceptible to infection with MHV-A59 and MHV-DVIM. Thus, the HE glycoprotein of an MHV strain did not lead to infection of cultured murine neural cells or of nonmurine cells that express the carbohydrate ligand of the HE glycoprotein. Therefore, interaction of the spike glycoprotein of MHV with its carcinoembryonic antigen-related receptor glycoprotein is required for infectivity of MHV strains whether or not they express the HE glycoprotein.