Loss of heterozygosity of essential genes represents a widespread class of potential cancer vulnerabilities.

Alterations in non-driver genes represent an emerging class of potential therapeutic targets in cancer. Hundreds to thousands of non-driver genes undergo loss of heterozygosity (LOH) events per tumor, generating discrete differences between tumor and normal cells. Here we interrogate LOH of polymorphisms in essential genes as a novel class of therapeutic targets. We hypothesized that monoallelic inactivation of the allele retained in tumors can selectively kill cancer cells but not somatic cells, which retain both alleles. We identified 5664 variants in 1278 essential genes that undergo LOH in cancer and evaluated the potential for each to be targeted using allele-specific gene-editing, RNAi, or small-molecule approaches. We further show that allele-specific inactivation of either of two essential genes (PRIM1 and EXOSC8) reduces growth of cells harboring that allele, while cells harboring the non-targeted allele remain intact. We conclude that LOH of essential genes represents a rich class of non-driver cancer vulnerabilities.

Introduction

Despite progress in precision cancer drug discovery, few highly selective therapies exist in the clinic. A current paradigm focuses on drugging driver alterations in cancer; however, many driver genes have proven difficult to target therapeutically[1,2], and in many cancers no easily targeted drivers exist. Alterations in non-driver genes represent an alternative target class that merits further investigation.

Loss of heterozygosity (LOH) may generate cancer-specific vulnerabilities by eliminating genetic redundancy in cancer cells. LOH occurs when a cancer cell that is originally heterozygous at a locus loses one of its two alleles at that locus, either by simple deletion of one allele (copy-loss LOH), or by deletion of one allele accompanied by duplication of the remaining allele (copy-neutral LOH). In either case, the cancer cell then relies on the gene products encoded by a single allele, in contrast to normal cells, which retain both alleles. When a cancer cell undergoes LOH of an essential gene, further loss or inhibition specifically of the allele retained in the tumor should not be tolerated, whereas normal cells will be able to survive relying solely on the remaining allele[3] (Fig. 1a). We term this target class GEMINI vulnerabilities, after the twins from Greek mythology Castor and Pollux, one of which was mortal and the other immortal.

Fig. 1

Genomic rates of LOH and allelic variation in normal and cancer genomes.

a Schematic indicating how loss of heterozygosity (LOH) of essential genes represents a potentially targetable difference between cancer and normal cells. b Violin plot of minor allele frequency of polymorphisms in essential versus non-essential genes in the ExAC cohort. Intersecting lines represent median values: essential = 0.141, non-essential = 0.146; one-tailed Student’s t-test, **p = 0.005. c (Left) Overlap between genes with common polymorphisms in the ExAC database (pink circle) and essential genes (blue circle). (Right) Fraction of essential genes with common polymorphisms. d Percent of genome affected by LOH across 9686 cancers from TCGA. e Stacked histogram representing the number of genes with copy-loss (yellow) or copy-neutral LOH (purple) across 9686 cancers from TCGA. f Dot plot of the number of essential genes affected by LOH across 33 TCGA tumor types. Tumor types are indicated by TCGA abbreviations (see https://gdc.cancer.gov/resources-tcga-users/tcga-code-tables/tcga-study-abbreviations). Each dot represents an individual sample. Lines indicate median values.

While previous reports have described individual GEMINI vulnerabilities[4,5], these studies have not systematically evaluated the landscape of potential targets, taking into account genome-scale assessments of gene essentiality, variation in human genomes, and rates of LOH across cancers. Open questions include which essential genes exhibit widespread variation in human populations and frequent LOH in cancers, providing potential GEMINI vulnerabilities, and at what rates these vulnerabilities occur. Moreover, different GEMINI vulnerabilities may require different therapeutic approaches to exploit them, due to the location of the variant within each gene and its effects on the amino acid composition of the protein. These differences have not been explored. Furthermore, GEMINI vulnerabilities have never been validated in isogenic systems to confirm specificity.

To address these questions, we integrated genome-scale copy number, germline allelic variation, and gene essentiality data to identify a list of polymorphisms in cell essential genes that undergo LOH in cancer, serving as a compendium of potential GEMINI targets. We also performed proof-of-principle validation of GEMINI vulnerabilities for two candidate genes in this list, PRIM1 and EXOSC8, using allele-specific CRISPR in both patient-derived and isogenic models. These results rigorously validate the GEMINI class of vulnerabilities and define its potential scope.

Results

Genome-wide identification of GEMINI vulnerabilities

To identify potential targets for our approach, we first characterized the landscape of cell-essential genes. We integrated genome-wide gene essentiality data from loss-of-function genetic screens and CCLE cell lines to conservatively estimate 1481 genes that are essential across lineages (Supplementary Data 1; Methods). This list is enriched for genes involved in essential cellular processes including rRNA processing, mRNA splicing, and translation (functional enrichment analysis performed with DAVID[6,7], version 6.8; Supplementary Data 2).

We then assessed germline heterozygosity resulting from normal human genetic variation in coding regions and 5′ and 3′ untranslated regions (UTRs) using allele frequencies across 60,706 individuals in the Exome Aggregation Consortium database[8]. Variants at 90,409 loci were observed to be present among at least 1% of alleles. As expected, polymorphisms in essential genes are slightly less common than those in non-essential genes (median minor allele frequency: essential = 0.141, non-essential = 0.146; p = 0.005, one-tailed Student’s t-test; Fig. 1b). However, essential genes still contain an abundance of genetic variation: 86% (1278/1481) harbor at least one common germline variant (Fig. 1c), with 49% (730/1481) harboring at least one missense variant. The median essential gene contains 3 germline polymorphisms. The median polymorphism in an essential gene is heterozygous in 13.9% of individuals (Supplementary Fig. 1a).

We were interested in how much of this heterozygosity in essential genes is lost in cancer. Loss of heterozygosity (LOH) in cancer frequently results from copy number alterations (CNAs) that can alter dozens to thousands of genes in cancer genomes[9,10]. Most LOH is due to strict copy-loss (copy-loss LOH), where allelic loss occurs in the context of a decrease in gene copy number. However, copy-neutral LOH is also frequently observed, whereby an allele is lost but the number of gene copies remains the same or in some cases even increases due to a duplication event. LOH has been frequently described[10,11], but to our knowledge there has not yet been a systematic analysis of the frequency of LOH events across cancer types.

We therefore analyzed copy number and LOH calls from 9686 patient samples across 33 TCGA tumor types (Methods)[12]. On average and across all cancers, 16% of genes undergo LOH (Fig. 1d). Genome-wide LOH rates vary widely by tumor type, ranging from a median of 45% in adenoid cystic carcinoma to 0.01% in thyroid carcinoma (Supplementary Fig. 1b). Approximately 28% of genes undergoing LOH undergo copy-neutral LOH (Fig. 1e), and on average across all cancers, 4.4% of all genes undergo copy-neutral LOH.

Rates of LOH are no lower for cell-essential genes relative to the rest of the genome (essential: 16.4%, non-essential: 15.6%; p = 1, one-tailed Student’s t-test; Supplementary Fig. 1c), suggesting that LOH of essential genes does not impose negative selection pressure. As a result, tumors harbored an average of 189 essential genes with LOH (Fig. 1f).

We hypothesized that the widespread nature of LOH of essential genes could represent a new opportunity to target essential genes that are heterozygous in normal tissue but undergo LOH in cancer. Among individuals with heterozygous SNPs within an essential gene, cancer cells with LOH of that gene would rely solely on the gene product encoded by one allele, in contrast to somatic cells, which would retain both alleles. We therefore hypothesized that allele-specific inactivation of the allele that had been retained in the cancer would selectively kill the cancer cells (Fig. 1a).

Our analysis identified 5664 polymorphisms in 1278 cell-essential genes, representing a compendium of potential GEMINI vulnerabilities (Supplementary Data 3). These GEMINI genes are enriched for similar pathways as the wider set of essential genes, including rRNA processing, mRNA splicing, and translation (functional enrichment analysis performed with DAVID[6,7], version 6.8; Supplementary Data 4). Among the 5664 GEMINI variants, 1688 lead to missense changes in amino acid composition of an essential protein, raising the possibility that they could be distinguished by molecules that interact with the protein directly. We focused on two of these missense SNPs for further functional analysis.

Validation of PRIM1rs2277339 as a GEMINI vulnerability

Variants residing in putative CRISPR protospacer adjacent motif (PAM) sites have previously been shown to enable allele-specific gene disruption[13-15]. For nuclease activity, S. pyogenes Cas9 requires a PAM site of the canonical motif 5′-NGG-3′ downstream of the 20-nucleotide target site; deviations from this motif abrogate Cas9-mediated target cleavage[16,17]. Therefore, we hypothesized that in the case in which one allele of a SNP generates a novel PAM site, Cas9 would be able to disrupt the “CRISPR-sensitive” (S), G allele that maintains the PAM sequence while leaving the other, “CRISPR-resistant” (R) allele intact (Fig. 2a).

Fig. 2

Validation of PRIM1rs2277339 as a GEMINI vulnerability.

a Schematic indicating allele-specific CRISPR approach. “Preexisting genome” represents individuals heterozygous for a germline SNP in a S. pyogenes Cas9 protospacer adjacent motif (PAM) site. A “G” allele (blue) in the PAM retains Cas9 activity at the target site, making this allele CRISPR-sensitive (S). An allele other than “G,” represented by “X” (red) abrogates Cas9 activity at the target site, making this allele CRISPR-resistant (R). Expression of an allele-specific (AS) CRISPR sgRNA targeting the polymorphic PAM site leads to specific inactivation of the S allele. b Schematic of PRIM1 SNP rs2277339 locus showing target sites for positive control, non-allele specific (NA) sgRNA and experimental, allele-specific (AS) sgRNA. Alleles appear in bold. c Crystal structure of PRIM1 gene product[88] shows the amino acid encoded by rs2277339 (teal) lies on the surface of the primase catalytic subunit (gray) near a potentially small-molecule accessible location. d Immunoblot of PRIM1 protein levels in indicated patient-derived cell lines expressing LacZ, PRIM1 NA, or PRIM1 AS sgRNA (n = 1 biological replicate). e Representative growth curves of indicated patient-derived cell lines expressing LacZ (black), PRIM1 NA (red), or PRIM1 AS (blue) sgRNA, as measured by CellTiter-Glo luminescence, relative to day of assay plating. n = 5 technical replicates. Data are presented as mean values ± s.d. See Supplementary Fig. 2 for additional biological replicates. f Representative growth curves of indicated isogenic cell lines expressing LacZ (black), PRIM1 NA (red), or PRIM1 AS (blue) sgRNA, as measured by CellTiter-Glo luminescence, relative to day of assay plating. n = 5 technical replicates. Data are presented as mean values ± s.d. See Supplementary Fig. 2 for additional biological replicates. g Disruption of PRIM1 in isogenic hemizygous PRIM1 resistant (PRIM1R) or PRIM1 sensitive (PRIM1S) cells expressing PRIM1 NA or AS sgRNA. Unaltered alleles (black), alleles with in-frame insertions or deletions (gray), and alleles with frameshift alterations (yellow) were assessed by deep sequencing of PRIM1 four days post-infection with sgRNA. Source data for Fig. 2d–g are provided as a Source Data file.

We identified such a SNP in the essential gene PRIM1 as a promising candidate for proof-of-principle validation. PRIM1 encodes the catalytic subunit of DNA primase and has previously been determined to be an essential gene[18-20]. It contains two common SNPs, of which one (rs2277339) leads to a change in the amino acid sequence: a T to G substitution resulting in conversion of an aspartate on the protein surface to an alanine (Fig. 2b–c, Supplementary Fig. 2a). The minor allele is common (minor allele frequency = 0.177), leading to heterozygosity at this locus in 29% of individuals represented in the ExAC database. This locus also undergoes frequent LOH. Across the 33 cancer types profiled, LOH was observed at the rs2277339 locus in 9% of cancers, including 21% of lung adenocarcinomas, 18% of ovarian cancers, and 17% of pancreatic cancers (Supplementary Fig. 2b).

PRIM1rs2277339 lies in a polymorphic PAM site—the “CRISPR-sensitive,” G allele generates a canonical S. pyogenes Cas9 PAM site, while the “CRISPR-resistant,” T allele disrupts the NGG PAM motif. We tested allele-specific PRIM1 disruption using an allele-specific (AS) CRISPR single guide RNA (sgRNA) designed to target only the G allele at rs2277339, encoding the alanine version of the protein (Fig. 2b). In the context of CRISPR experiments, because the G allele should be sensitive to allele-specific disruption, we use an “S” to designate cells with this allele and an “R” to designate cells with the other, resistant allele: for example, PRIM1S/– and PRIM1R/S genotypes reflect cells with one copy of the sensitive G allele and cells with one copy of each allele, respectively.

We transduced four patient-derived cancer cell lines that naturally exhibit either rs2277339 allele with AS sgRNA and verified that AS sgRNA disrupts PRIM1 in PRIM1S genetic contexts (Fig. 2d). PRIM1S/– and PRIM1S/S cells expressing AS sgRNA show decreased proliferation relative to LacZ-targeting control, whereas cells retaining the resistant allele (PRIM1R/–, PRIM1R/R, or PRIM1R/S) show no such defects (Fig. 2e, Supplementary Fig. 2c–f).

The specificity of the AS sgRNA for PRIM1S cell lines was not due to a lack of Cas9 activity or PRIM1 essentiality in the PRIM1R cell lines. We confirmed this finding by transducing four cell lines with a non-allele specific (NA) PRIM1-targeting sgRNA. We successfully ablated PRIM1 expression in all contexts (Fig. 2d), and cells expressing PRIM1-targeting sgRNA showed dramatic and significant decreases in proliferation relative to LacZ-targeting control even in cases where expression of the AS sgRNA did not significantly limit growth (p < 0.01 in all cases, one-tailed Student’s t-test; Fig. 2e, Supplementary Fig. 2c–f).

We further tested isogenic cell lines harboring either allele. Using SNU-175 and SNU-C4 cells, which are heterozygous for PRIM1rs2277339, as a base, we transiently transfected a vector expressing Cas9 and two sgRNAs that flank the PRIM1 gene. We then screened single cell clones for PRIM1 deletion by PCR. Among deletion-positive clones, we identified heterozygous (PRIM1R/S), hemizygous sensitive (PRIM1S), and hemizygous resistant (PRIM1R) lines (Supplementary Fig. 3a–e). Using these isogenic cells, we confirmed PRIM1S/– cells expressing AS sgRNA show decreased proliferation relative to LacZ-targeting control, whereas cells retaining the resistant allele (PRIM1R/– or PRIM1R/S) show no such defects (Fig. 2f, Supplementary Fig. 2g–k).

Within these isogenic lines, we also confirmed that AS CRISPR disrupts PRIM1 in a PRIM1rs2277339-dependent manner using deep sequencing. In order to assay PRIM1 knockout efficiency before cells started to die due to loss of the essential PRIM1 protein, we performed deep sequencing on DNA collected from cells at the early time point of four days post-infection. At this early time point, isogenic PRIM1S and PRIM1R cells infected with the NA sgRNA showed comparable fractions of disrupted alleles (Fig. 2g), suggesting both lines exhibit similar levels of Cas9 activity. However, while PRIM1S cells expressing AS sgRNA showed approximately 40% disrupted alleles, resistant cells under the same condition showed 0 disrupted alleles (p < 0.0001, Chi-square with Yates correction; Fig. 2g). This result confirms that AS PRIM1 sgRNA targets PRIM1 in a SNP-specific manner. We also verified that allele-specific inactivation of essential genes is possible in a heterozygous context (Supplementary Fig. 2l, Supplementary Note 1).

Validation of EXOSC8rs117135638 as a GEMINI vulnerability

We also performed proof-of-principle validation for another candidate SNP in the essential gene EXOSC8. EXOSC8 codes for Rrp43, a component of the RNA exosome. The RNA exosome is an essential multi-protein complex involved in RNA degradation and processing, including processing of pre-rRNA[21-23]. Two common SNPs have been described within EXOSC8, one of which (rs117135638) represents a C to A change in DNA sequence; this SNP leads to a proline to histidine substitution on the interface between Rrp43 and exosome complex member Mtr3 (Fig. 3a, b, Supplementary Fig. 4a). This candidate SNP is heterozygous in 2% of individuals and undergoes LOH in 29% of cancers, including 72% of lung squamous cell carcinomas, 62% of ovarian cancers, 46% of lung adenocarcinomas, and 40% of breast cancers (Supplementary Fig. 4b).

Fig. 3

Validation of EXOSC8rs117135638 as a GEMINI vulnerability.

a Schematic of EXOSC8 SNP rs117135638 locus showing target sites for positive control, non-allele specific (NA) sgRNA and experimental, allele-specific (AS) sgRNA. Alleles appear in bold. b Crystal structure of EXOSC8 gene product, Rrp43[89] (gray) shows the amino acid encoded by rs117135638 (teal) lies on the surface of the Rrp43 protein near the interface with exosome complex subunit Mtr3 (orange). c Disruption of EXOSC8 in patient-derived EXOSC8 resistant (EXOSC8R) or EXOSC8 sensitive (EXOSC8S) cells expressing EXOSC8 non-allele specific (NA) positive control sgRNA or allele-specific (AS) experimental sgRNA. Unaltered alleles (black), alleles with in-frame insertions or deletions (gray), and alleles with frameshift alterations (yellow) were assessed by deep sequencing of EXOSC8 four days post-infection with sgRNA. d Immunoblot of EXOSC8 protein levels in indicated patient-derived and isogenic cell lines expressing LacZ, EXOSC8 NA, or EXOSC8 AS sgRNA (n = 2 technical replicates of 1 biological sample). e Representative growth curves of indicated patient-derived and isogenic cell lines expressing LacZ (black), EXOSC8 NA (red), or EXOSC8 AS (blue) sgRNA, as measured by CellTiter-Glo luminescence, relative to day of assay plating. n = 5 technical replicates. Data are presented as mean values ± s.d. See Supplementary Fig. 4 for additional biological replicates. f Immunoblot of EXOSC8 protein levels in indicated isogenic cell lines expressing LacZ, EXOSC8 NA, or EXOSC8 AS sgRNA (n = 2 technical replicates of 1 biological sample). g Representative growth curves of indicated isogenic cell lines expressing LacZ (black), EXOSC8 NA (red), or EXOSC8 AS (blue) sgRNA, as measured by CellTiter-Glo luminescence, relative to day of assay plating. n = 5 technical replicates. Data are presented as mean values ± s.d. See Supplementary Fig. 4 for additional biological replicates. Source data for Fig. 3c–g are provided as a Source Data file.

We first tested allele-specific (AS) EXOSC8 disruption using an sgRNA designed to target only the C allele at rs117135638, encoding the proline version of the protein. We designated cells as EXOSC8S (for “sensitive”) if they contained this allele, and as EXOSC8R (for “resistant”) if they contained the A allele. Using patient-derived cell lines, we verified the allele-specific effects of EXOSC8 AS sgRNA on both the DNA and the protein level (Fig. 3c, d, Supplementary Note 2). Consistent with these observations, EXOSC8S/– and EXOSC8S/S cells expressing AS sgRNA showed decreased proliferation relative to LacZ-targeting control, whereas cells retaining the resistant allele (EXOSC8R/– or EXOSC8R/R) showed no such defects (Fig. 3e, Supplementary Fig. 4c–f, Supplementary Note 3).

We next determined that both copy-loss and copy-neutral LOH of EXOSC8 represents a vulnerability in an isogenic context. We generated diploid and single-copy knockout isogenic cells representing EXOSC8S and EXOSC8R genotypes by Cas9-mediated homology-directed repair (HDR) editing (Methods; Supplementary Fig. 3f–i), and then infected these isogenic Cas9-stable lines with constructs expressing EXOSC8 NA or AS sgRNA. As expected, EXOSC8 NA sgRNA ablated EXOSC8 expression in all contexts, while AS sgRNA ablated EXOSC8 expression only in EXOSC8S cells (Fig. 3f). EXOSC8S/S and EXOSC8S/Δ cells expressing AS sgRNA showed decreased proliferation relative to LacZ-targeting control, whereas cells retaining the resistant allele (EXOSC8R/R and EXOSC8R/Δ) showed no such defects (Fig. 3g, Supplementary Fig. 4g–j).

Potential approaches to targeting GEMINI vulnerabilities

We were interested in understanding the potential scope of patients that could benefit from therapeutic approaches targeting GEMINI vulnerabilities. For each GEMINI variant, we calculated the number of new patients in the US per year that exhibit LOH of the hypothetical “targetable” allele (Methods). Across the 33 tumor types we profiled, the median GEMINI vulnerability could be targetable in 17,747 US patients per year (Fig. 4a). PRIM1 rs2277339 and EXOSC8rs117135638 could be targetable in a theoretical 22,470 and 5,307 patients per year in the US, respectively (Supplementary Figs. 2a and 4a).

Fig. 4

Potential therapeutic approaches to targeting GEMINI vulnerabilities.

a Number of GEMINI variants (vertical axis) plotted against the number of patients per year in the US whose tumors might respond to therapeutics targeting those variants (i.e., have lost the resistant allele from a heterozygous germline; horizontal axis). Bin width = 1000 patients. b Growth of heterozygous (red circles) and hemizygous cells (pink circles) expressing positive control, non-allele specific PRIM1-targeting shRNAs versus PRIM1 mRNA expression. Cell growth measured by CellTiter-Glo luminescence relative to day 2 post-infection and shGFP (n = 5 technical replicates). PRIM1 mRNA expression assessed by qRT-PCR (n = 3 technical replicates). Data are presented as mean values ± s.d. Dashed gray line indicates PRIM1 expression threshold below which substantial decreases in cell viability are observed. c Summary table representing challenges to developing allele-specific small molecules that target GEMINI vulnerabilities and associated analyses to prioritize targets. Source data for Fig. 4b is provided as a Source Data file.

The major challenge to exploiting GEMINI vulnerabilities is identifying means to target them in humans. Three approaches that may be contemplated are DNA-targeting CRISPR effectors (e.g., Cas9), RNA-targeting approaches (e.g., RNAi), and allele-specific small molecules. We characterized each GEMINI vulnerability according to criteria that would indicate its amenability to targeting by each of these approaches and, when practical, performed proof-of-concept in vitro experiments testing each approach.

We first focused on CRISPR effectors as a potential therapeutic modality for targeting GEMINI vulnerabilities. Because we had already demonstrated that S. pyogenes Cas9 can disrupt PRIM1 and EXOSC8 in an allele-specific manner, we next analyzed the list of GEMINI vulnerabilities to identify the full range of polymorphisms whose targeting on the DNA level may enable allele-selective gene disruption by a CRISPR-based approach. For this analysis, we included both the canonical S. pyogenes PAM, NGG, as well as the weaker, non-canonical PAM, NAG[17,24]. Of the 4648 GEMINI vulnerabilities in open reading frames, 23% (1088/4648, or 19% of all GEMIMI vulnerabilities) generate a PAM site in one allele but not the other, suggesting the potential for allele-specific knockout (Supplementary Data 3).

In theory, every GEMINI variant could be the target of allele-specific RNAi reagents. However, it is possible that, for some GEMINI variants, RNAi reagents would be unable to suppress expression sufficiently to reduce cell viability, or that sufficient allelic specificity might not be achieved. For example, we tested the hypothesis that PRIM1rs2277339 may be targetable in an allele-specific manner using RNAi. For these experiments, we refer to the PRIM1 alleles by their identifying nucleotide; for example, heterozygous cells are referred to as PRIM1T/G, and cells hemizygous for the minor allele are referred to as PRIM1G/–. We first sought to determine the level of PRIM1 knockdown required to substantially decrease cell proliferation. Accordingly, we infected hemizygous and heterozygous isogenic cells with non-allele specific PRIM1-targeting shRNAs and assessed PRIM1 expression and cell growth. We observed that substantial decreases in cell growth were possible under conditions of robust PRIM1 knockdown (>80%) (Fig. 4b).

We then asked whether allele-specific shRNAs targeting the PRIM1rs2277339 locus could decrease growth in cells representing the fully matched genotype. PRIM1T/– and PRIM1G/– cells were infected with constructs encoding fully complementary shRNAs tiling across the SNP and assessed for growth. Only one shRNA, shG7 (targeting the minor, G allele at position 7) significantly reduced cell growth relative to GFP-targeting control (Supplementary Fig. 5a, b). We then selected the four PRIM1 SNP-targeting shRNAs that yielded the lowest average cell growth relative to GFP-targeting control and assessed their ability to decrease cell growth in an allele-specific manner. Heterozygous cells (PRIM1T/G) and hemizygous cells of the targeted genotype (PRIM1T/– or PRIM1G/–) were infected with constructs encoding the appropriate shRNA. No putative allele-specific shRNAs were found to significantly decrease cell growth in hemizygous cells of the targeted genotype relative to heterozygous cells (Supplementary Fig. 5c, d). We conclude that PRIM1rs2277339 may not represent an optimal candidate for allele-specific shRNA-mediated inhibition.

Given the large number of additional GEMINI variants that may be suitable for RNAi-mediated targeting, we sought to prioritize GEMINI genes that may be amenable to allele-specific inhibition using mRNA-targeting approaches. RNAi-mediated knockdown of some essential genes may be more effective at inducing cell death than others, based on differential expression thresholds required for cell survival[25-27]. We hypothesized that GEMINI genes representing strong dependencies in shRNA screens would be most amenable to potential targeting using an RNAi-based therapeutic. We therefore analyzed shRNA data representing 17,212 genes in 712 cell lines[28], including 1183 GEMINI genes, and looked for genes whose suppression led to at least a moderately strong response in most of the cell lines (median DEMETER2 score < −0.5; Methods). Approximately 35% of GEMINI genes (413/1183), including PRIM1 (median DEMETER2 score = −0.52), fit this category, representing 35% (1804/5196) of GEMINI vulnerabilities (Supplementary Data 3). In comparison, only 3.6% of all genes profiled (623/17,212) passed this dependency threshold, indicating a significant enrichment for GEMINI genes (p < 0.0001, binomial proportion test). However, this level of dependency was not observed for all GEMINI or essential genes. For example, the median DEMETER2 score for EXOSC8 was only −0.14, despite our and others’ extensive data showing its essentiality in multiple cell types[21,22,23]. These results raise the possibility that RNAi-based approaches may not be able to exploit many GEMINI vulnerabilities.

Both CRISPR- and RNAi-based therapeutic approaches suffer from difficulties in effectively delivering reagents to all cancer cells in an animal. Small molecule–based approaches can overcome such delivery issues, but substantial obstacles exist to developing allele-specific small molecules that target GEMINI vulnerabilities. These challenges include identifying GEMINI genes that are amenable to small-molecule inhibition, determining which GEMINI variants lie near potentially druggable pockets, and predicting which GEMINI variants are most likely to facilitate allele-specific drug binding. In comparison to the CRISPR- and RNAi-based approaches explored above, allele-specific small molecules are more tractable for delivery/efficacy. However, no clear allele-specific small molecules exist for any of the 1749 protein-altering GEMINI variants identified in our analysis. Therefore, we focused on in silico analyses to identify and prioritize these GEMINI vulnerabilities (missense, insertion, and deletion variants) for potential allele-specific drug development (Supplementary Data 5).

To identify GEMINI genes that may be amenable to small-molecule inhibition, we annotated those containing protein-altering alleles using the canSAR Protein Annotation Tool[29,30] (cansar.icr.ac.uk; Methods). This tool uses publicly available structural and chemical information to generate structure- and ligand-based druggability scores. While these scores do not necessarily reflect potential for allele-specific small molecule inhibition based on the GEMINI variant of interest, they may nonetheless allow prioritization of targets based on general druggability. This analysis found that of the 1734 protein-altering variants in genes assessed by canSAR, 12% (212) reside in proteins with a small molecule ligand–bound structure (Fig. 4c). Additional assessments of potential small-molecule binding sites on structures with and without existing ligands indicated that 39% of protein-altering variants (674) lie in proteins with molecular structures that are predicted to be druggable (drug-like compound modulates activity in vivo) or tractable (tool compound modulates activity in vitro) (Fig. 4c). Furthermore, 160 GEMINI variants reside in proteins in the top 90th percentile of ligand-based druggability as assessed by the physiochemical properties of small molecules tested against the protein or its homologs (Fig. 4c). We also found that 25% of protein-altering GEMINI variants (441/1734) reside in enzymes as defined by their annotation with an Enzyme Commission (EC) number[31] (Fig. 4c).

To assess which variants may reside in protein regions amenable to small molecule binding, we performed a p-blast of the 1749 protein-altering variants against protein sequences for molecular structures found in the Protein Data Bank[32] (rcsb.org; Methods). This analysis identified 153 variants characterized in a homologous structure. We then visually scored 81 missense and indel variants in X-ray crystal structures for their proximity to solvent-exposed pockets or known small-molecule binding sites using a scale of 0 to 4 (Methods). Of the variants analyzed, 15 were near a potential binding pocket on the surface of the protein (score = 3), with two of these pockets containing a small molecule ligand (score = 4) (Fig. 4c).

We also assessed protein-altering GEMINI variants to prioritize those that may be most amenable for allele-specific small-molecule inhibition. For this analysis, we scored variants that altered the number or sign of residue charges. For example, a variant that changes the charge of a residue from neutral to negative or that adds an additional negative charge through an inserted residue would qualify as a charge-altering variant. Of the 1749 protein-altering variants, 584 induced a change in residue charge (Fig. 4c). We further hypothesized that variants introducing a cysteine residue could provide additional allele selectivity by enabling the potential development of a covalent inhibitor. Among the missense and indel GEMINI variants, 95 generate a cysteine in one allele.

We then integrated each of these analyses to characterize the potential druggability landscape of these protein-altering GEMINI vulnerabilities. Every variant was given a score from 0 to 7 based on the number of analyses in which it scored among the top candidates. One variant, TGS1rs7823773, earned a score of 6, including in the visual scoring and cysteine categories. Nine additional variants earned a score of 5 (Supplementary Data 5). These may be among the highest-priority candidates for further exploration.

Discussion

Leveraging synthetic lethal interactions in cancer cells represents a promising avenue to targeting genomic differences between tumor and normal tissue. Synthetic lethality between genes occurs when singly inactivating one gene or the other maintains viability, but inactivating both genes simultaneously causes lethality[33]. Over the past 20 years, many efforts have been directed toward discovering synthetic lethal interactions with genetic driver alterations of oncogenes and tumor suppressor genes[34,35]. However, the number of genetically activated oncogenes and inactivated tumor suppressor genes in any given tumor is limited and, in many cancer types, is vastly outnumbered by genetically altered non-driver genes (e.g., due to passenger events). Therefore, identifying synthetic lethalities with genetic alterations affecting non-driver genes (also termed “collateral lethalities”[36]) would increase the scope of potential therapeutic approaches. While individual GEMINI genes have been described previously[3], our work integrated genome-wide assessments of gene essentiality, genetic variation, and LOH to generate the first systematic analysis of this target class.

GEMINI vulnerabilities represent one of four classes of collateral lethalities. In addition to GEMINI vulnerabilities, deletion of paralogs can result in dependency on the remaining paralog; loss or gain of function of a non-driver pathway can lead to dependencies on alternative non-driver pathways;[36] and hemizygous loss of essential genes can result in dependency on the remaining copy (CYCLOPS)[25,26].

Prior analyses have indicated CYCLOPS genes to be the most frequent class of these synthetic lethal interactions[26,27], but we find that GEMINI vulnerabilities provide similar numbers of potential targets. In comparison, fewer paralog dependencies have been described[27,37-40]. Larger numbers of paralog vulnerabilities have been predicted[41], but it is unclear whether these predictions represent viable candidates[36]. The 1278 GEMINI genes that we identified also exceed the 299 known driver genes[42], many of which are proposed therapeutic targets. Expanding the search for GEMINI vulnerabilities beyond pan-essential genes to include variants in lineage-specific essential genes may also increase the number of potential GEMINI targets.

In comparison to CYCLOPS, targeting GEMINI vulnerabilities has two distinct advantages. First, whereas CYCLOPS genes must lie in regions of copy loss, GEMINI genes encompass genes that undergo both copy-loss and copy-neutral LOH. Second, while CYCLOPS vulnerabilities rely on relative differences between tumor and normal cells (differential expression of target genes), GEMINI vulnerabilities exploit absolute differences (the presence or absence of the allele that has undergone LOH). Thus, the possibility of allele-specific targeting presented by GEMINI genes may widen prospective therapeutic windows. In 269 cases, GEMINI vulnerabilities we detected reside in CYCLOPS genes[26] (Supplementary Data 3). If GEMINI and CYCLOPS vulnerabilities are additive, targeting these genes might offer an even wider therapeutic window in cancers where CYCLOPS-GEMINI genes suffer LOH due to copy loss. However, individual GEMINI alterations may be less common among patients than individual CYCLOPs alterations due to the requirement that the germline genome be heterozygous at the GEMINI locus.

Like any target class, we expect resistance mechanisms to arise in response to targeting GEMINI vulnerabilities. Base pair substitutions that replace the targeted allele with an alternative are likely to occur in one in every 108–109 cells, given observed mutation rates per cell division[43]. Additional alterations affecting nearby nucleic or amino acids could interfere with genetic (e.g., CRISPR and RNAi) or protein-targeting approaches. It is also possible that alternative pathways exist for some GEMINI genes whereby alterations of other genes compensate for inhibition of a GEMINI gene. However, our list of GEMINI genes is highly enriched for components of universal cellular processes, such as DNA, RNA, and protein biogenesis, for which no alternative pathways exist to compensate their loss[44].

Biomarkers for detection of patients who may benefit from a GEMINI approach are relatively straightforward: one would select patients who are heterozygous for the targeted allele and for whom the tumor is found to have lost the alternative allele. One consideration is tumor heterogeneity; if the LOH event is present in only part of the tumor, resistance would be expected to arise quickly. However, in prior analyses[43,45-47], a majority of somatic copy-number alterations, including LOH events, appeared to be clonal, although the fraction of clonal events can be lower in some loci in some tissues[48]. One approach to minimize clonal variation in LOH is to prioritize GEMINI genes that lie on chromosomes or chromosome arms that are characteristically lost early in oncogenesis (e.g., 3p in renal clear cell carcinoma)[37].

While we show that cells heterozygous for PRIM1rs2277339 exhibit no substantial proliferation defects upon ablation of the targeted allele (Fig. 2e, f), systemic knockout of one allele of an essential gene across all cells in a patient is not likely to be a tractable therapeutic strategy. Thus, potential allele-specific gene editing approaches to leverage GEMINI vulnerabilities in the clinic would rely on a highly cancer cell–specific delivery system to avoid knockout of the targeted allele in normal tissue. While much work remains to achieve the necessary targeting specificity, advances in nanoparticle delivery systems present the possibility of targeting Cas9 DNA, mRNA, or protein in a tumor-specific manner[49-52]. Additionally, S. pyogenes Cas9 enzymes with altered[53] or expanded[54,55] PAM specificities or CRISPR effectors from other species[53,56-58] have broadened the total number of targetable loci in the genome and, thereby, the number of targetable variants.

GEMINI vulnerabilities could be also targetable by reversible genetic approaches. Such reversible genetic inhibitors have seen recent success in other disease indications (e.g., the RNAi-based patisiran for treatment of hereditary transthyretin amyloidosis;[59] the antisense oligonucleotide [ASO]–based IONIS-HTTRX for treatment of Huntington’s Disease[60]). Notably, early studies of the GEMINI genes POLR2A and RPA1 achieved allele-specific growth suppression of cancer cells using ASO[4,61,62] and RNAi[63] reagents. Peptide nucleic acids, or PNAs, which can suppress both transcription and translation of target genes, represent another potential allele-specific genetic approach[64-67]. Finally, the RNA-targeting CRISPR effector Cas13 has shown the ability to knock down target genes[68] and decrease proliferation of cancer cells[69] in an allele-specific manner. Unlike Cas9, the Cas13 enzyme from L. wadei previously used for allele-specific RNA cleavage does not require a downstream PAM-like motif[68], potentially expanding the number of targetable sites beyond those tractable with DNA-targeting CRISPR effectors. The use of genetic targeting approaches would dramatically increase the number of potentially targetable GEMINI vulnerabilities by including silent as well as protein-altering variants. However, unlike the promising therapies for genetic disorders mentioned above, a genetic inhibitor of a GEMINI vulnerability would need to be delivered to all cells in a tumor. Thus, like Cas9-based modalities, the use of Cas13 and other reversible genetic approaches to exploit GEMINI vulnerabilities would require the development of novel delivery systems.

Allele-specific small molecule inhibitors present another attractive possibility for drugging GEMINI vulnerabilities. Allele-specific therapeutics in clinical use include rationally designed drugs (e.g., mutant EGFR inhibitors[70]) as well as those whose genotype-specific effects were identified by pharmacogenomic studies (e.g., warfarin and VKORC1[71]). However, GEMINI vulnerabilities present an additional challenge for allele-specific inhibitor development because most variants in cell-essential genes do not reside in or near an active site (Fig. 4c) or other functionally critical protein region. This challenge may be addressed through alternative small-molecule approaches, such as proteolysis-targeting chimera (PROTAC)-mediated degradation[72]. SNPs for which one allele is a cysteine could be prioritized for this approach because of the possibility of engineering a covalent inhibitor[73].

While we have rigorously validated PRIM1 and EXOSC8 as genetic dependencies in cancer, further work is necessary to explore potential therapeutic modalities for targeting them (Supplementary Discussion). The design of a tractable therapeutic that targets any single GEMINI gene in an allele-specific manner is a substantial challenge. However, the sheer number of potential candidates suggests that some of these GEMINI vulnerabilities may represent viable targets.

Methods

Variant lists

A list of 228,440 potentially targetable variants was downloaded from the Exome Aggregation Consortium (ExAC) database (exac.broadinstitute.org)[8]. Potentially targetable variants were defined as those in the following classes: 3_prime_UTR_variant, 5_prime_UTR_variant, frameshift_variant, inframe_deletion, inframe_insertion, initiator_codon_variant, missense_variant, splice_acceptor_variant, splice_donor_variant, splice_region_variant, stop_gained, stop_lost, stop_retained_variant, synonymous_variant. These variants were filtered to include only PASSing, common variants (global minor allele frequency ≥ 0.01) in genes for which copy number calls were available through the NCI Genomic Data Commons (see below for further details of copy number analyses).

All variant classes were included in the analysis of potential target SNPs for reversible genetic therapeutic approaches. All variant classes except 3_prime_UTR_variant and 5_prime_UTR_variant were included in the determination of variants generating or disrupting an S. pyogenes PAM site.

Genomic analyses of copy number and LOH from TCGA

Patient-derived genome-wide copy number and LOH data were downloaded from the TCGA Pan-Can project via the NCI Genomic Data Commons (https://gdc.cancer.gov/about-data/publications/pancan-aneuploidy) first reported in[12]. For copy number, gene-level log2 relative data were calculated by GISTIC 2.0, referenced in the output file “all_data_by_genes_whitelisted.tsv”. Copy-loss was defined as log2 relative values ≤−0.1 and copy-neutral was defined as >−0.1.

For LOH calls, we used TCGA analyses[12]. Briefly, SNP array and exome sequencing data from both tumor samples and paired normal DNA were used as inputs to ABSOLUTE[74], which calculated absolute allelic copy numbers genome-wide for each tumor. These absolute allelic copy numbers took into account the purity and ploidy of each tumor, as determined by ABSOLUTE. Autosomal regions for which the absolute copy number of one allele was zero were considered to have undergone LOH. These ABSOLUTE calls are in the file “TCGA_mastercalls.abs_segtabs.fixed.txt” (https://gdc.cancer.gov/about-data/publications/pancan-aneuploidy). These calls were transformed into per-gene calls for all subsequent analyses.

Essential gene list

Candidate essential genes were nominated using data from three genome-scale loss-of-function screens of haploid human cell lines (KBM7 with CRISPR-Cas9 gene inactivation or mutagenized with gene trapping[75], and pluripotent stem cells with CRISPR-Cas9 gene inactivation[76]). Briefly, all genes that passed a threshold of <10% FDR for a given cell line were included as a candidate essential gene. FDR-corrected p-values from the original publications were used for both CRISPR screens; FDR q-values for the KBM7 gene trap scores were calculated using a binomial model (representing equal probability of gene trap inserting in a sense versus anti-sense orientation) and correction for multiple hypotheses using Benjamini and Hochberg. This initial candidate list contained 3431 genes, with 633 scoring as essential in all three screens.

These candidate essential genes were then filtered using CCLE gene copy-number and RNA expression data to determine if loss-of-function genetic alterations were observed in human cell lines. Genes that met any of the following criteria were excluded: homozygously deleted in >2 cell lines (log2 copy-number <−5); very low RNA expression (<0.5 RPKM) in >5 cell lines; or homozygously deleted in 1 cell line that also has low RNA expression (<1.0 RPKM). This analysis reduced the list to 2566 candidate essential genes. Genes were then filtered based on mean CERES score from CRISPR knockout screens of 517 cell lines (https://depmap.org/portal/download; derived from the file “gene_effects.csv”)[77]. Genes with CERES scores >−0.4 were excluded, yielding a list of 1499 essential genes. To account for instances in which CCLE copy-number and/or RNA expression data were not available for a particular gene, genes were rescued from the CCLE filter if they scored as essential in two of the three haploid cell line screens and had mean a CERES score <−0.4. This rescue yielded 17 genes, bringing the total number of candidate essential genes to 1516. This list was further filtered to remove genes classified as Tier 1 tumor suppressor genes in the COSMIC Cancer Gene Census (https://cancer.sanger.ac.uk/census/)[78], yielding a final list of 1482 essential genes. (One essential gene in this list, AK6, was not characterized in TCGA copy-number and LOH data and so was excluded from further analyses.)

Cell line identification and cell culture

Human cancer cell lines of the appropriate genotypes for PRIM1 and EXOSC8 were identified using whole exome sequencing and absolute gene copy number data from the Cancer Cell Line Encyclopedia (https://portals.broadinstitute.org/ccle)[79]. All lines were genotyped for the SNP of interest using Sanger sequencing. Cell lines were maintained in RPMI-1640 supplemented with 10% fetal bovine serum and 1% penicillin. Lines were not assessed for contamination with mycoplasma. No commonly misidentified cell lines defined by the International Cell Line Authentication Committee have been used in these studies.

Plasmids

lentiCas9-Blast (Addgene plasmid # 52962) and lentiGuide-Puro (Addgene plasmid # 52963) were gifts from Feng Zhang[80]. A Cas9 construct co-expressing GFP and two sgRNAs was a gift from Peter Choi[26]. pLKO.1–TRC cloning vector was a gift from David Root (Addgene plasmid # 10878)[81].

CRISPR sgRNAs

To identify target sites for CRISPR-Cas9–mediated knockout, the genetic sequences of PRIM1 and EXOSC8 were obtained from the UCSC genome browser (http://genome.ucsc.edu) using the human assembly GRC38/Hg38 (December 2013). The 20 nucleotides upstream of the polymorphic PAM site containing the SNP for each gene constitutes the AS sgRNA for that gene. All other sgRNAs were designed using the CRISPR sgRNA design tool from the Zhang lab (http://crispr.mit.edu). sgRNAs were cloned into the appropriate vector as described previously[80,82]. Briefly, plasmids were cut and dephosphorylated with BsmBI (New England Biolabs) and FastAP (Fermentas) at 37 °C for 2 h. Oligonucleotides for each sgRNA guide sequence (Integrated DNA Technologies) were phosphorylated using T4 polynucleotide kinase (New England Biolabs) at 37 °C for 30 min and then annealed by heating to 95 °C for 5 min and cooling to 25 °C at 1.5 °C/min. Using Quick Ligase (New England Biolabs), annealed oligos were ligated into gel purified vectors (Qiagen) at 25 °C for 5 min. Cloned plasmids were amplified using a endotoxin-free maxi-prep kit (Qiagen).

The sgRNA sequences were as follows:

LacZ: GTTCGCATTATCCGAACCAT

PRIM1 AS: CAGCTCGGGCAGCTCGGTGG

PRIM1 NA: CGCTGGCTCAACTACGGTGG

EXOSC8 AS: CGGAATCTCGATGAACACAG

EXOSC8 NA: ACCGGAATCTCGATGAACAC

Cell growth assays

Cells were plated in opaque 96-well plates (Corning) at 500, 1000, or 2500 cells per well on the indicated day post–lentiviral infection. Cell number was inferred by ATP-dependent luminescence using CellTiter-Glo reagent (Promega) and normalized to the relative luminescence on the day of plating.

Generation of PRIM1-loss and EXOSC8-loss cells

A Cas9 construct co-expressing GFP and two sgRNAs with target sites flanking PRIM1 was used to delete a 20.6 kb region encoding PRIM1. Cell lines heterozygous for PRIM1rs2277339 (SNU-C4 and SNU-175) were transfected with this construct using LipoD293 transfection reagent (SignaGen), and single GFP+ cells were sorted by FACS and plated at low density for single-cell cloning or single-cell sorted into 96-well tissue culture plates containing a 50:50 mix of conditioned and fresh RPMI-1640 media, 20% serum, 1% penicillin-streptomycin, and 10 µM ROCK inhibitor Y-27632. Clones were expanded and validated by PCR to harbor the 20.6 kb deletion encoding PRIM1, and the retained allele was genotyped by Sanger sequencing. These clones were designated SNU-175PRIM1 S/–, SNU-175PRIM1 R/–, SNU-C4PRIM1 S/– for subsequent experiments. Other clones were determined by PCR and Sanger sequencing to retain both PRIM1 alleles and not to harbor this deletion and were designated as control cell lines for subsequent experiments (SNU-175PRIM1 R/S and SNU-C4PRIM1 R/S). The same procedure was employed using a cell line diploid for the EXOSC8R SNP (SW-579) to generate EXOSC8R/– cell lines harboring a 7.1 kb deletion and EXOSC8R/R control lines.

The sgRNA sequences were as follows:

PRIM1 upstream: GCGCGGAACTCGCCACGGTA

PRIM1 downstream: CAGAGCTCCTCAAACCATTG

EXOSC8 upstream: GGTTTCTCGGCCGAGCGCCG

EXOSC8 downstream: TGTACCCATCTACTTAAGTT

Primers used to verify gene deletion by PCR were as follows:

PRIM1 deletion genotyping F: ACTGTATGCACCACCACACC

PRIM1 deletion genotyping R: AGTTCACGTGGAGCATCCTT

EXOSC8 deletion genotyping F: TTTGGGGCATACTCATGCTT

EXOSC8 deletion genotyping R: TCCACCTCCAATTATTTGTTCC

Generation of EXOSC8 isogenic cell lines

Cas9 RNPs and a ssODN repair template were used to edit the EXOSC8S SNP to the EXOSC8R SNP. S. pyogenes Cas9-NLS (Synthego) and an sgRNA (sequence: ACCGGAATCTCGATGAACAC) targeting the EXOSC8 SNP region (Synthego) were complexed as described previously[83]. Briefly, 100 pmol of Cas9-NLS was diluted to a final volume of 5 µL with Cas9 buffer (20 mM HEPES [pH 7.5], 150 mM KCl, 1 mM MgCl2, 10% glycerol, and 1 mM TCEP) and mixed with 5 µL of Cas9 buffer containing 120 pmol of sgRNA. This mixture was incubated for 10 min at RT to allow RNP formation. DV-90 cells (EXOSC8S/S) were nucleofected with resulting RNPs, a 50:50 mix of EXOSC8S and EXOSC8R ssODN (IDT), and a GFP-expressing plasmid (pMAX-GFP) (Lonza). ssODN repair templates contained a synonymous mutation introducing a novel Mnl1 restriction site for downstream genotyping as well as a silent blocking mutation to prevent repeated Cas9 cleavage. Single GFP+ cells were single-cell sorted by FACS into 96-well tissue culture plates containing a 50:50 mix of conditioned and fresh RPMI-1640 media, 20% serum, 1% penicillin-streptomycin, and 10 µM ROCK inhibitor Y-27632. Clones were expanded and evaluated for HDR-mediated editing by PCR and restriction digest, and positive clones were genotyped by next-generation sequencing (NGS; MGH DNA Core).

CRISPR variant sequencing

Cellular pellets were collected from Cas9-stable cells 4 or 18 days post-infection with lentiGuide-Puro virus encoding the indicated sgRNA. Genomic DNA was isolated using a DNAMini kit (Qiagen), and the target region for each gene was amplified by PCR (EMD Millipore). Amplicons were submitted to NGS CRISPR sequencing by the MGH DNA Core. Non-altered alleles as well as those containing in-frame or frameshift indels were determined manually using the CRISPR variant output file. PCR primer sequences were as follows:

PRIM1 MGH F: GCACAGAAGGCGCTTCATA

PRIM1 MGH R: CGCCAATTCCTGTGGTAATC

EXOSC8 MGH F: AGCTGCAGAGTGTTTCTTTCA

EXOSC8 MGH R: AGAGCAAAGTAAATGAAAAGCCCAA

Western blotting

Cells were washed in ice-cold PBS and lysed in 1x RIPA buffer (10 mM Tris-Cl pH 8.0, 1 mM EDTA, 1% Triton X-100, 0.1% SDS, and 140 mM NaCl) supplemented with 1x protease and phosphatase inhibitor cocktail (PI-290, Boston Bioproducts). Lysates were sonicated in a bioruptor (Diagenode) for 5 min (medium intensity) and cleared by centrifugation at 15,000g for 15 min at 4°C. Proteins were electrophoresed on polyacrylamide gradient gels (Life Technologies) and detected by chemiluminescence (Bio-rad).

Antibodies used were as follows:

EXOSC8: Proteintech #11979-1-AP

PRIM1: Cell Signaling Technology #4725

Vinculin: Sigma #V9131

See Source Data file for unprocessed scans of the most important blots.

shRNA sequences

pLKO.1 GFP shRNA (target sequence: GCAAGCTGACCCTGAAGTTCAT) was a gift from David Sabatini (Addgene plasmid # 30323)[84]. Lentiviral expression constructs for non-allele specific shRNA-mediated suppression of PRIM1 were obtained through the Broad Institute of MIT and Harvard Genomic Perturbation Platform (https://portals.broadinstitute.org/gpp/public/). The names, clone IDs, and target sequences used in our studies are as follows:

shPRIM1 (TRCN0000275194): AGCATCGTCTCTGGGTATATT

TRCN0000151860: CCGAGCTGCTTAAACTTTATT

TRCN0000275194: AGCATCGTCTCTGGGTATATT

TRCN0000275195: GATTGATATAGGCGCAGTATA

TRCN0000275196: CCGAGCTGCTTAAACTTTATT

Allele-specific shRNA sequences were cloned into the vector pLKO.1 as described previously[81]. Briefly, the pLKO.1 plasmid was cut with AgeI and EcoRI (New England Biolabs) at 37 °C for 2 h. Oligonucleotides for each shRNA sequence (Integrated DNA Technologies) were annealed by heating to 95 °C for 5 min and cooling to 25 °C at 1.5 °C/min. Using T4 ligase (New England Biolabs), annealed oligos were ligated into gel-purified vectors (Qiagen) at 25 °C for 30 min. Cloned plasmids were amplified using a endotoxin-free maxi-prep kit (Qiagen). The shRNA sequences were as follows:

PRIM1rs2277339 major-allele (T) targeting:

sh3T: TCAATGGAGACGTTTGACC

sh4T: CAATGGAGACGTTTGACCC

sh5T: AATGGAGACGTTTGACCCC

sh6T: ATGGAGACGTTTGACCCCA

sh7T: TGGAGACGTTTGACCCCAC

sh8T: GGAGACGTTTGACCCCACC

sh9T: GAGACGTTTGACCCCACCG

sh10T: AGACGTTTGACCCCACCGA

sh11T: GACGTTTGACCCCACCGAG

sh16T: TTGACCCCACCGAGCTGCC

PRIM1rs2277339 minor-allele (G) targeting:

sh3G: TCAATGGAGACGTTTGCCC

sh4G: CAATGGAGACGTTTGCCCC

sh5G: AATGGAGACGTTTGCCCCC

sh6G: ATGGAGACGTTTGCCCCCA

sh7G: TGGAGACGTTTGCCCCCAC

sh8G: GGAGACGTTTGCCCCCACC

sh9G: GAGACGTTTGCCCCCACCG

sh10G: AGACGTTTGCCCCCACCGA

sh11G: GACGTTTGCCCCCACCGAG

sh16G: TTGCCCCCACCGAGCTGCC

Quantitative and reverse transcription PCR

RNA was extracted using the RNeasy Mini kit (Qiagen) and subjected to on-column DNase treatment. cDNA was synthesized with the Superscript II Reverse Transcriptase kit (Life Technologies) with no–reverse-transcriptase samples serving as negative controls. Gene expression was quantified by Power Sybr Green Master Mix (Applied Biosystems). PRIM1 expression values were normalized to vinculin (VCL) and the fold change calculated by the DDCt method. Primers used in our studies are as follows:

PRIM1-F: GCTCAACTACGGTGGAGTGAT

PRIM1-R: GGTTGTTGAAGGATTGGTAGCG

VCL-F: CGCTGAGGTGGGTATAGGTG

VCL-R: TTGGATGGCATTAACAGCAG.

Calculations of theoretical patient numbers

To determine number of patients in the US that could benefit from a therapeutic approach targeting each GEMINI vulnerability, we used the following formula:

in which

κ = # new pan-cancer cases per year in the US (1,735,350)

χ = rate of heterozygosity of GEMINI variant

λ = pan-cancer rate of LOH of GEMINI gene

0.5 = fraction of patients with LOH that undergo loss of theoretical “targetable”

allele, assuming that the allele lost during LOH is random

Π = theoretical number of new patients per year in the US.

Estimate of new US pan-cancer cases per year derived from SEER Cancer Statistics Review[85]. Rate of heterozygosity estimated using 2pq from Hardy-Weinberg equation[86,87].

DEMETER2 analyses

For a detailed description of the screening and analysis methodology used to generate DEMETER2 scores, please see[28]. Briefly, DEMETER2 generates an absolute dependency score for each gene suppressed in each cell line. A score of 0 signifies no dependency and a score of 1 signifies a strong dependency as estimated by scaling the effect to a panel of known pan-essential genes. DEMETER2 scores were obtained from the Cancer Dependency Map Portal (https://depmap.org/portal/download/) using the file “D2_combined_gene_dep_scores.csv”. We classified genes that exhibited a median DEMETER2 score of ≤−0.5 across all cell lines as moderately strong dependencies.

canSAR protein annotation

The canSAR protein annotation tool (cansar.icr.ac.uk) was run on a list of 741 unique genes containing 1749 insertion, deletion, and missense variants. Structures with >90% sequence homology were included in structural druggability and chemical matter analyses.

Determination of structures corresponding to variants

To determine which variants were present in PDB, DNA sequences (30mer) encapsulating 1749 insertion, deletion, and missense variants were translated in all 6 frames using the Bio.Seq Python module. Output was blasted using the Bio.Blast Python module against the PDB database with E-value thresholds of 0.001 or less, resulting in hits for 267 variants. We manually curated these structures to verify the presence of the variant within the PDB file and eliminated structures for which correspondence between the PDB protein sequence, ExAC amino acid prediction, and UCSC Genome Browser amino acid sequence was inconclusive. This curation yielded 153 protein-altering variants in proteins with homologous molecular structures.

Visual scoring was performed on 81 protein-altering variants that lie in X-ray crystal structures. Variants were scored using the following scale: 0 = no clear pockets on the protein surface, 1 = SNP far from pocket on protein surface, 2 = SNP near pocket on protein surface, 3 = SNP in pocket on protein surface, 4 = SNP near pocket containing small molecule.