Xuesong Li1, Nigel J O'Neil1, Noushin Moshgabadi1, Philip Hieter2. 1. Michael Smith Laboratories, University of British Columbia, Vancouver, British Columbia V6T 1Z4, Canada. 2. Michael Smith Laboratories, University of British Columbia, Vancouver, British Columbia V6T 1Z4, Canada hieter@msl.ubc.ca.
Abstract
Many tumors contain mutations that confer defects in the DNA-damage response and genome stability. DNA-damaging agents are powerful therapeutic tools that can differentially kill cells with an impaired DNA-damage response. The response to DNA damage is complex and composed of a network of coordinated pathways, often with a degree of redundancy. Tumor-specific somatic mutations in DNA-damage response genes could be exploited by inhibiting the function of a second gene product to increase the sensitivity of tumor cells to a sublethal concentration of a DNA-damaging therapeutic agent, resulting in a class of conditional synthetic lethality we call synthetic cytotoxicity. We used the Saccharomyces cerevisiae nonessential gene-deletion collection to screen for synthetic cytotoxic interactions with camptothecin, a topoisomerase I inhibitor, and a null mutation in TEL1, the S. cerevisiae ortholog of the mammalian tumor-suppressor gene, ATM. We found and validated 14 synthetic cytotoxic interactions that define at least five epistasis groups. One class of synthetic cytotoxic interaction was due to telomere defects. We also found that at least one synthetic cytotoxic interaction was conserved in Caenorhabditis elegans. We have demonstrated that synthetic cytotoxicity could be a useful strategy for expanding the sensitivity of certain tumors to DNA-damaging therapeutics.
Many tumors contain mutations that confer defects in the DNA-damage response and genome stability. DNA-damaging agents are powerful therapeutic tools that can differentially kill cells with an impaired DNA-damage response. The response to DNA damage is complex and composed of a network of coordinated pathways, often with a degree of redundancy. Tumor-specific somatic mutations in DNA-damage response genes could be exploited by inhibiting the function of a second gene product to increase the sensitivity of tumor cells to a sublethal concentration of a DNA-damaging therapeutic agent, resulting in a class of conditional synthetic lethality we call synthetic cytotoxicity. We used the Saccharomyces cerevisiae nonessential gene-deletion collection to screen for synthetic cytotoxic interactions with camptothecin, a topoisomerase I inhibitor, and a null mutation in TEL1, the S. cerevisiae ortholog of the mammaliantumor-suppressor gene, ATM. We found and validated 14 synthetic cytotoxic interactions that define at least five epistasis groups. One class of synthetic cytotoxic interaction was due to telomere defects. We also found that at least one synthetic cytotoxic interaction was conserved in Caenorhabditis elegans. We have demonstrated that synthetic cytotoxicity could be a useful strategy for expanding the sensitivity of certain tumors to DNA-damaging therapeutics.
CHANGES in genome structure and sequence underlie oncogenesis (Charames and Bapat 2003; Schvartzman ). These tumor-specific somatic mutations can be exploited to target tumor cells for killing relative to normal cells (Hartwell ). Genome stability genes, such as TP53, ATM, BRCA1, BRCA2, and MSH2, are frequently mutated in tumors and many current antitumor therapeutics exploit genome stability defects through the use of DNA-damaging genotoxic agents (Helleday ). The resultant DNA-damage results in cell-cycle arrest and cell death. However, DNA-damaging therapeutic agents often have a low therapeutic index relative to toxicity and can cause deleterious side effects and generate new mutations that can result in therapeutic resistance or secondary cancers.Over a decade ago, two important concepts were proposed to facilitate the development of new anticancer drugs. First, that the somatic mutations in cancers could selectively sensitize tumor cells to therapies that inhibit the function of a second gene product resulting in synthetic lethality (SL). Second, that SL interactions, which are defined between two genes when disruption of function of either gene product is viable but disruption of function of both simultaneously results in death, could be screened in genetically amenable model organisms to identify those that may be relevant to the treatment of tumors (Hartwell ).Not all tumors contain mutations that can be exploited by treatment with DNA-damaging agents (because, for example, of redundancy in the DNA-damage response), or by synthetic lethal approaches (because, for example, SL partners do not exist or the SL partners are not “druggable” targets). It is possible that these tumors could be sensitive to a combination of the two approaches. In a manner similar to a SL interaction, a somatic tumor-specific mutation together with inhibition of a second gene product could increase the sensitivity of tumor cells to a low, sublethal concentration of DNA-damaging agent resulting in a conditional synthetic lethality that we are calling synthetic cytotoxicity (SC) (Figure 1). For instance, cells with mutations affecting a DNA-damage repair pathway may rely on parallel or adapted DNA-repair pathways when treated with DNA-damaging agents (Bandyopadhyay ; Guenole ). Hence, modulation of DNA-damage responses genetically with mutations or chemically with small molecule inhibitors of DNA-repair enzymes could selectively enhance the sensitivity of cancer cells to DNA-damaging therapies resulting in SC.
Figure 1
Schematic of cytotoxic therapy, synthetic lethality, and synthetic cytotoxicity. Selective killing of tumor cells using DNA-damaging therapeutic agents and DNA-repair enzyme inhibitors. Cytotoxic therapy: A mutation in the chromosome stability gene CIN1 sensitizes tumor cells to DNA-damaging agents. Unrepaired induced DNA damage leads to cell death. Synthetic lethality: A mutation in CIN2 is synthetic lethal with inhibition of DDR1, a DNA-damage response (DDR) protein. Endogenous DNA damage cannot be repaired in the absence of both CIN2 and DDR1 and leads to cell death. This outcome is analogous to the synthetic lethality observed when cells with BRCA1 or BRCA2 mutations are treated with PARP inhibitor. Synthetic cytotoxicity: Inhibition of DDR1 in the CIN3 mutant background does not result in synthetic lethality but loss of function of both proteins sensitizes the tumor cell to a low sublethal dose of DNA-damaging agent, enhancing the differential killing of tumor cells relative to a CIN3 wild-type background.
Schematic of cytotoxic therapy, synthetic lethality, and synthetic cytotoxicity. Selective killing of tumor cells using DNA-damaging therapeutic agents and DNA-repair enzyme inhibitors. Cytotoxic therapy: A mutation in the chromosome stability gene CIN1 sensitizes tumor cells to DNA-damaging agents. Unrepaired induced DNA damage leads to cell death. Synthetic lethality: A mutation in CIN2 is synthetic lethal with inhibition of DDR1, a DNA-damage response (DDR) protein. Endogenous DNA damage cannot be repaired in the absence of both CIN2 and DDR1 and leads to cell death. This outcome is analogous to the synthetic lethality observed when cells with BRCA1 or BRCA2 mutations are treated with PARP inhibitor. Synthetic cytotoxicity: Inhibition of DDR1 in the CIN3 mutant background does not result in synthetic lethality but loss of function of both proteins sensitizes the tumor cell to a low sublethal dose of DNA-damaging agent, enhancing the differential killing of tumor cells relative to a CIN3 wild-type background.Mapping the large number of genetic interactions needed to identify SC in human cells is feasible but techniques are not as robust as those currently available in budding yeast. Synthetic genetic arrays (SGA) in yeast facilitate the collection and analysis of genetic interaction data (Tong ; Collins ; Costanzo ). The use of the model organisms, yeast, and Caenorhabditis elegans, to identify conserved genetic interactions with potential cancer therapeutic value has proven effective (McLellan ; McManus ; van Pel ).As a proof of principle study of SC, we screened the Saccharomyces cerevisiae collection of nonessential gene deletions for SC interactions with the topoisomerase I poison camptothecin (CPT) and a deletion of the yeastATM (ataxia telangiectasia mutated) ortholog, was chosen in our study to represent a frequently mutated tumor-associated genome stability gene. The ATM gene is mutated in the hereditary cancer-prone human disease ataxia telangiectasia (AT) (Savitsky ). Cells from AT patients exhibit radiosensitivity, altered DNA repair, cell-cycle checkpoint defects, short telomeres, and chromosome instability (Lavin. 2008). Although AT is rare, with a frequency of ∼1/40,000, up to 1% of the population could be heterozygous carriers of ATM mutations, and these carrier individuals are estimated to have an increased risk of pancreatic cancer and breast cancer (Savitsky ; Thompson ; Renwick ; Roberts ). In addition to hereditary cancer, ATM mutations are frequently detected in diverse sporadic cancers (Stankovic ; Gumy-Pause ; Hall. 2005; Kang ). Importantly for this study, the function of ATM in genome stability appears to be well conserved between in yeast, worm, and human (Fritz ; Garcia-Muse and Boulton 2005; Jones ).We performed digenic SGA screens with as a query mutation in the presence or absence of CPT to identify common processes required for resistance to camptothecin when Tel1 is absent. We identified and validated 14 gene mutations that resulted in SC with and CPT. We also demonstrated that at least one SC interaction is conserved in C. elegans. This study demonstrates the utility of model organisms in screening for SC to DNA-damaging chemotherapeutics and raises the potential for detecting candidate combination therapies in simple model organisms.
Materials and Methods
Media and growth conditions
Yeast was grown in rich media at 30°. Plasmid-bearing strains were grown in synthetic complete media lacking the appropriate nutrient. All strains are BY4743 background (Brachmann ). Single-gene knockout alleles were generated through tetrad dissection using the heterozygous gene deletion collection as starting resource (kindly provided by Jef Boeke) (Pan ). All double-mutant heterozygous strains were then constructed by mating each of the respective single mutants. All haploid double-mutant strains were generated by sporulation of diploid heterozygous double-mutant strains followed by tetrad dissection. The genotypes of all the strains used in the experiment were checked by PCR for confirmation of the gene deletions (Supporting Information, Table S1). was deleted to maintain the viability of cells with (Zhao ).
Synthetic genetic array construction and analysis
SGA analyses were performed using a Singer RoToR as described in Tong . The MATa yeast deletion mutant set (ykoΔ::kanMX) was arrayed at a density of 1536 colonies/plate by robotic pining. The media for SGA were the same as described in Tong . MATa deletion mutant array (DMA) for SGA was obtained from Research Genetics Company (http://www.sequence.stanford.edu/group/yeast_deletion_project/deletions3.html). Statistical analysis was performed in R (J. Stoepel, J. Bryan, K. Ushey, B. P. Young, C. J. Loewen, and P. Hieter, unpublished results) using a script that compares the colony size measurements for each gene pair. The program output (EC value) is a measure of the difference in colony size, a proxy for strain fitness, of corresponding colonies on the single- and double-selection plates, and a measure of statistical significance (P-value).
Plate assays and growth curve analysis
To identify SC interactions we utilized qualitative plate assays as well as quantitative growth-curve analyses comparing the agent sensitivity of double-mutant strains with the two single-mutant parental strains. Different camptothecin concentrations were used to maximize growth differences for single-, double-, and triple-mutant strains and allow for identification of additive and suppressing effects. For plate assays, cells were grown in YPD media until log phase before spotting on YPD plates and growth assessed by visual inspection. For growth-curve analysis, cells were grown in YPD media until log phase before addition to YPD liquid media containing sublethal doses of the DNA-damaging agents. Cell density was normalized by OD600 readings. For each agent, the concentration was optimized to noticeably inhibit at least one single parental mutant strain growth and further diluted 4–8× to determine the sensitivity of the double mutants compared to wild type. For plate spot assays, an identical optical density (OD600) of cells was serially diluted fivefold and spotted on the indicated plate at the indicated temperature for 72 hr. For growth-curve analysis, logarithmic phase cultures were diluted to an OD600 value of 0.15 in 96-well plates and grown for 24 hr in a TECAN M200 plate reader at 30°. Each strain was tested in three replicates. Strain fitness F was based on the area under the curve (AUC) of mutants relative to wild type as previously described (McLellan ). Genetic interactions were quantified through the comparison of fitness of double mutants and the expected phenotype based on the two single mutants. Mutations in independent genes (two genes with a neutral interaction) often combine to generate fitness (growth relative to WT) in a multiplicative manner. The expected fitness of the resulting double mutant is assumed to be Fexpected = F × F. The interaction can be quantitatively measured by comparing double mutant Fab against Fabexpected (Mani ).
C. elegans brood and camptothecin assays
C. elegans strains were cultured at 20° under standard conditions. Bristol N2 was used as the wild-type strain. The following mutations were used in this study: , , , and . Hatching rates and male production were assessed by individually plating 10–20 L4 larvae of each genotype and transferring the animals every 24 hr until the onset of sterility. The number of embryos was counted immediately after transferring and the number of adult worms was counted 48 hr later. The percentage of males among adults was determined. To assess CPT sensitivity, early adults were treated with 0, 5, 10, or 20 nM CPT in M9 buffer containing E. coliOP50 for 19 hr. Fifty animals were plated 10 per plate on an OP50 seeded NGM plate and cultured to produce progeny. After 2 hr, treated P0 animals were removed. The number of inviable embryos and adults was scored 24 and 72 hr later, respectively. Percentage embryonic survival was calculated by dividing the total number of adults by the total progeny (adults + inviable embryos).
C. elegans DAPI staining
Young adult hermaphrodites were picked to a watch glass containing 10–20 μl of M9 buffer. A total of 500 μl 96% ethanol containing 0.1μg/ml DAPI was added to the watch glass and incubated at room temperature for 30 min in the dark. Fixed and stained worms were washed with 2 ml M9 buffer 3× (30 min per wash) before mounting on 1% agarose-coated microscope slide. The condensed meiotic chromosomes were visualized using a 100× objective on a Zeiss Axioplan 2 microscope.
Results
tel1Δ rad27Δ is synthetic cytotoxic to camptothecin
Although Tel1 plays a role in the DNA-damage response, few synthetic lethal interactions have been identified with mutants (Stark ). We tested whether loss of a second DNA-repair gene sensitized the mutant to sublethal concentrations of DNA-damaging agents resulting in SC. The function of Tel1 partially overlaps with the related kinase Mec1 (Shiloh. 2003; Budd ; Chakhparonian ); therefore, we tested both and deletions for SC to DNA-damaging agents in combination with deletions in three DNA-repair enzyme genes: (tyrosyl-DNA-phosphodiesterase 1), (three prime phosphatase 1), and (radiation sensitive 27). Viable double mutants were tested for SC by comparing the growth of double mutant strains to the two single mutant parental strains using a quantitative liquid growth curve assay in the presence or absence of four different DNA-damaging agents: bleomycin (BLEO), 5-fluorouracil (5FU), CPT, and hydroxyurea (HU). Δ Δ was synthetic lethal, which is consistent with previous reports (Pan ). Δ Δ was SC to sublethal concentrations of CPT but not to other DNA-damaging agents (Figure 2A). Quantitative liquid growth-curve analysis at a range of CPT concentrations confirmed that the double-mutant growth was worse than what was predicted from a multiplicative model based on the effect on each single mutant (Figure 2, B and C).
Figure 2
Synthetic cytotoxicity of rad27Δ tel1Δ to CPT. (A) Sensitivity to CPT of single and double mutants in spot assays. rad52Δ (positive control) and rad27Δ tel1Δ mutants are hypersensitive to CPT. tel1Δ is mildly sensitive compared to wild type, while rad27Δ shows a mild growth defect. (B) Fitness data from growth-curve assay demonstrating the synthetic cytotoxicity of rad27Δ tel1Δ to CPT. Quantitative measurement of genetic interaction between rad27Δ and tel1Δ. (C) Plot showing the analysis of growth-curve values derived from area under the curve (AUC). The expected fitness values of double mutants are obtained by multiplying two single mutant fitness values together. The actual fitness of rad27Δ tel1Δ in the presence of 1 μg/ml CPT is ∼7% lower than expected, indicating that rad27Δ tel1Δ is SC to CPT. Error bars represent standard deviation.
Synthetic cytotoxicity of rad27Δ tel1Δ to CPT. (A) Sensitivity to CPT of single and double mutants in spot assays. rad52Δ (positive control) and rad27Δ tel1Δ mutants are hypersensitive to CPT. tel1Δ is mildly sensitive compared to wild type, while rad27Δ shows a mild growth defect. (B) Fitness data from growth-curve assay demonstrating the synthetic cytotoxicity of rad27Δ tel1Δ to CPT. Quantitative measurement of genetic interaction between rad27Δ and tel1Δ. (C) Plot showing the analysis of growth-curve values derived from area under the curve (AUC). The expected fitness values of double mutants are obtained by multiplying two single mutant fitness values together. The actual fitness of rad27Δ tel1Δ in the presence of 1 μg/ml CPT is ∼7% lower than expected, indicating that rad27Δ tel1Δ is SC to CPT. Error bars represent standard deviation.
SC interactions with tel1Δ and CPT identified by SGA analysis
The SC of Δ Δ to CPT led us to screen for SC interactions with CPT and Δ on a genome-wide scale. We constructed double mutants by mating a Δ query strain to an ordered array of ∼4800 nonessential gene deletion mutants that represent ∼80% of all yeast genes (Tong and Boone 2006). The single- and double-deletion haploid strains were replicated onto media containing 4 µg/ml CPT (Figure S1). SC was identified when a double mutant exhibited a growth defect in the presence of CPT that was greater than that predicted using a multiplicative model based on the growth of single mutants (Figure S2).We identified 173 double mutants that showed significant growth defects compared to single mutants (P-value <0.05), which included all potential SC/SL interaction genes as well as single mutants that are sensitive to CPT. After removing single mutants that were sensitive to CPT and five false-positive interactions that were due to linkage to the query gene , the SGA screen resulted in 22 candidates (EC < −0.4 P-value < 0.05) as potential SL/SC interactions partner genes with Δ (Table 1).
Table 1
SL/SC candidate validation
Function
Gene
EC value
P-value
Random spore
Growth curve
Interaction
Shu complex
CSM2
−0.4399
0.004
+
+
SC
PSY3
−0.4109
<0.001
+
SC
SHU1
−0.2937
0.136
+
SC
SHU2
−0.2430
0.054
+
SC
Ku complex
YKU70
−0.4369
0.013
+
SC
YKU80
−0.2227
0.020
+
SC
DNA helicase
RRM3
−0.5544
<0.001
+
+
SC
PIF1
0.32135
0.659
—
NI
Spindle assembly checkpoint (SAC)
LTE1
−0.4332
0.029
+
+
SC
BFA1
−0.3709
0.060
+
+
SC
BUB2
−0.3067
0.019
+
+
SC
BUB3
−0.4356
0.346
+
+
SC
Casein kinase 2
CKB2
−0.4374
0.002
+
+
SC
CKA2
−0.3128
0.117
—
NI
CKB1
−0.2738
0.052
+
+
SC
CKA1
−0.0636
0.779
—
NI
Other function
ERG5
−0.4210
0.014
+
+
SC
YPR1
−0.5277
0.047
—
NI
CKI1
−0.8473
<0.001
—
NI
MVB12
−0.5303
0.006
—
NI
RUB1
−0.4708
0.006
—
NI
FPK1
−0.4598
0.005
—
NI
TRE1
−0.4370
0.003
—
NI
AGX1
−0.4214
<0.001
—
NI
RAS1
−1.128
<0.001
—
NI
HBT1
−1.1812
<0.001
—
NI
Unknown function
YJR128W
−0.7980
0.031
—
NI
YBR134W
−0.4113
0.033
—
NI
YCL046W
−0.4728
0.008
—
NI
YPR063C
−0.5446
0.014
—
NI
YOR223W
−0.4290
0.044
—
NI
RRT6
−0.5046
0.037
—
NI
The underlined text indicates interactions that did not satisfy the initial P-value <0.05 or interaction magnitude (EC) value <−0.4 cutoffs for SGA analysis but were chosen for further analysis based on the fact that they are part of multisubunit complexes for which at least one component was identified as SC with tel1Δ and CPT. For random spore and tetrad dissection assays, + and − indicate confirmed true and false positives, respectively. NI, no genetic interaction. The absence of a symbol indicates that the interaction was not tested.
The underlined text indicates interactions that did not satisfy the initial P-value <0.05 or interaction magnitude (EC) value <−0.4 cutoffs for SGA analysis but were chosen for further analysis based on the fact that they are part of multisubunit complexes for which at least one component was identified as SC with tel1Δ and CPT. For random spore and tetrad dissection assays, + and − indicate confirmed true and false positives, respectively. NI, no genetic interaction. The absence of a symbol indicates that the interaction was not tested.
Validation and expansion to related genes
The 22 candidate SC genes identified by SGA were tested by random spore assay or tetrad dissection. Newly constructed haploid strains containing gene-deletion alleles were derived from a heterozygous diploid collection (a gift from Jef Boeke) and used to construct the double-deletion strains for validation. Using alleles from a different source than those used in the SGA increased confidence in the interaction. Newly derived double mutants were tested by spot assays and quantitative liquid growth curve analysis. Of the 22 candidate gene mutations tested, seven showed significant SC interactions with in these validation assays (Table 1).Five of the seven SC gene mutations were components of multisubunit complexes or processes: and encode components of the SHU complex, which contains two other subunits, Shu1 and Shu2; encodes a component of the spindle assembly checkpoint; encodes a protein in the nonhomologous end joining (NHEJ) DNA-repair pathway; and encodes one of the β-subunits of the tetrameric casein kinase 2 (CK2). Other components of these pathways or complexes were not identified as SC according to our cutoff criteria in the SGA screen (Table 1). To test if these genes were false negatives, we constructed double mutants with and assayed for SC to CPT by spot assays and growth-curve analysis. We found that deletions of the SHU complex genes and , the NHEJ gene , the casein kinase 2 subunit but not or , and the spindle assembly checkpoint genes , , and resulted in SC to CPT (Figure 3 and Table 1).
Figure 3
Validation of SC interactions identified by SGA. Serial dilutions of cells were spotted on YPD plates containing 16 µg/ml CPT to test SC.
Validation of SC interactions identified by SGA. Serial dilutions of cells were spotted on YPD plates containing 16 µg/ml CPT to test SC.Many DNA-repair gene mutants are very sensitive to CPT. These genes were filtered out during analysis because the single mutant was inviable when exposed to 4 µg/ml CPT. To determine if interactions were missed due to sensitivity of single mutants, we used a SGA mini-array containing ∼200 DNA-repair gene mutants and scored for SC with at two concentrations of CPT (1 and 4 µg/ml). At 4 µg/ml CPT, the mini-array identified all of the SC interactions that were identified in the whole nonessential genome SGA and three subunits of the SWI/SNF nucleosome-remodeling complex, , , , that were not detected in the whole nonessential genome screen (Table 2). At 1 µg/ml CPT, the mini-array identified SC interactions with four subunits of the PCNA-like DNA-damage checkpoint complex, , , , , that were not identified at 4 µg/ml because the single mutants are exquisitely sensitive to CPT. The SC interactions with PCNA-like checkpoint genes are consistent with a previous study (Guenole )
Table 2
Additional SC interactions identified using a mini-array and screening with 1 µg/ml and 4 µg/ml camptothecin
Function
Gene
EC value
P-value
1 µg/ml CPT
4 µg/ml CPT
SWI–SNF complexes
SNF2
−0.7816
0.007
—
SC
SNF5
−0.2278
0.008
—
SC
SNF6
−0.3235
0.304
—
SC
RRM3
−0.3059
0.042
SC
RRM3
−0.7822
0.003
SC
BUB3
−0.3898
0.005
—
SC
SHU complex
SHU1
−0.3732
0.002
—
SC
CSM2
−0.3269
0.004
—
SC
PSY3
−0.2323
0.012
—
SC
CHL1
−0.2994
0.002
—
SC
RecQ–Top3 complex
RMI1
−0.2913
0.030
—
SC
SGS1
−0.2653
0.041
—
SC
RDH54
−0.2362
<0.001
—
SC
YKU70
−0.2300
0.008
—
SC
PCNA-like checkpoint
RAD17
−0.6640
<0.001
SC
—
RAD24
−0.6585
<0.001
SC
—
DDC1
−0.5103
0.002
SC
—
MEC3
−0.4901
<0.001
SC
—
HDA1
−0.2562
0.001
SC
—
SRS2
−0.2122
0.014
SC
—
SSN8
−0.2008
0.007
SC
—
Due to differences in array composition, a larger percentage of slow-growing or DNA-damage-sensitive mutant strains, an EC value <−0.2 was used as a cutoff for SC for the mini-array SGA analysis. Underlined text indicates that the gene was also identified in the larger scale SGA screen.
Due to differences in array composition, a larger percentage of slow-growing or DNA-damage-sensitive mutant strains, an EC value <−0.2 was used as a cutoff for SC for the mini-array SGA analysis. Underlined text indicates that the gene was also identified in the larger scale SGA screen.
SC is specific to DNA damage and CPT
To test if the SC interactions were specific to DNA damage or were common to other stresses, we tested the fitness of the double mutants in the presence of cycloheximide (CHX), which inhibits protein synthesis. None of the double mutants were SC to CHX, indicating that the identified SC interactions were not due to a general response to stress. Only and were sensitive to treatment with CHX (Figure 4). Cells deleted for genes encoding ergosterol biosynthetic enzymes show sensitivity to many different drugs (Mukhopadhyay ; Parsons ) suggesting that the SC to CPT observed in was the result of increased intracellular levels of CPT coupled with the loss of Tel1, which is required for resistance to CPT.
Figure 4
SC with tel1Δ is dependent on the specific DNA-damaging agent. The sensitivity of SC interactions to camptothecin (CPT), hydroxyurea (HU), and cycloheximide (CHX) tested in spot assays.
SC with tel1Δ is dependent on the specific DNA-damaging agent. The sensitivity of SC interactions to camptothecin (CPT), hydroxyurea (HU), and cycloheximide (CHX) tested in spot assays.We next tested whether the SC was specific to CPT by assaying the double mutants for sensitivity to the replication inhibitor HU. Both HU and CPT increase DNA replication stress during S phase leading to DNA damage. However, HU acts by a different mechanism than CPT does. HU is an inhibitor of ribonucleotide reductase, slows replication fork progression, and induces replication fork stalling by reducing dNTP pools (Lopes ; Alvino ). We found that , and were SC to both CPT and HU, suggesting that these double mutants are sensitive to replication stress in general. In contrast, the SC interactions between and the CKB subunit mutations, SHU complex mutations, and the Ku complex mutations were specific to CPT (Figure 4).
Epistasis of SC groups
To determine if the SC to CPT was due to a common mechanism, we constructed triple mutants with and pairwise combinations of the interacting genes to determine if the SC to CPT was additive or epistatic. Based on sensitivity to CPT we were able to define five genetic epistasis groups: SHU complex genes , , , ; Ku genes , ; the spindle checkpoint genes , , , ; casein kinase 2 β-subunit genes , ; and the DNA helicase gene . Triple mutants displayed additive sensitivity to CPT in the background when knocking out two genes, one gene from each group (Figure S3).
The SC of tel1Δ yku70Δ and tel1Δ yku80Δ is not due to loss of NHEJ
DSBs can be repaired by either homologous recombination (HR), which is the predominant form of repair during late S/G2 phase, or NHEJ, which is the predominant form of repair during G1/early S phase (Mao ). As core proteins in NHEJ repair machinery, the Ku70–Ku80 heterodimer is recruited rapidly to DSBs and forms a ring-like structure that binds to DNA ends and initiates NHEJ (Daley ). Both and when combined with were SC to CPT and as expected triple mutants did not exhibit more severe SC to CPT (Figure S3). Although both and were SC to CPT, a deletion of the NHEJ-associated DNA ligase, , did not result in SC to CPT when combined with (Figure 5A). As Dnl4 is essential for NHEJ, the SC to CPT of and or was not due to the loss of NHEJ suggesting an NHEJ-independent mechanism. Recent work has demonstrated that Ku proteins have roles beyond NHEJ repair. Ku proteins protect dsDNA break ends from resection and potentially restrict HR activity (Langerak ). However, deletion of or , which are the nucleases required for end resection, did not alleviate the SC of suggesting that SC was not the result of increased DNA end resection (Figure S4).
Figure 5
Investigating the mechanism of SC. (A) Comparison of growth of tel1Δ yku70Δ and tel1Δ dnl4Δ mutants on CPT in spot assays. (B) Growth of tel1Δ yku70Δ, tel1Δ yku80Δ, and tel1Δ yku70Δ yku80Δ mutants containing pRS415–CDC13–EST1 or vector alone with and without CPT. (C) Comparison of sensitivity to CPT of tel1Δ psy3Δ, tel1Δ shu1Δ double mutants in the presence and absence of Srs2.
Investigating the mechanism of SC. (A) Comparison of growth of tel1Δ yku70Δ and tel1Δ dnl4Δ mutants on CPT in spot assays. (B) Growth of tel1Δ yku70Δ, tel1Δ yku80Δ, and tel1Δ yku70Δ yku80Δ mutants containing pRS415–CDC13–EST1 or vector alone with and without CPT. (C) Comparison of sensitivity to CPT of tel1Δ psy3Δ, tel1Δ shu1Δ double mutants in the presence and absence of Srs2.
SC to CPT of tel1Δ yku70Δ and tel1Δ yku80Δ is suppressed by telomere lengthening
Mutations in , , , or result in the shortening of telomeres (Boulton and Jackson 1998; Ivessa ). CPT treatment induces telomeric DNA damage in human cells (Kang ). CPT may also affect telomeres in yeast and exacerbate the telomere-shortening defects observed in , , and mutants. To examine whether telomere elongation can suppress the SC to CPT, we expressed a Cdc13–Est1 fusion protein, which promotes telomere elongation (Evans and Lundblad 1999), in and . Increasing telomerase activity by expression of Cdc13–Est1 protein suppressed SC to CPT in , , and but did not suppress SC to CPT in (Figure 5B). These data suggest that the SC to CPT in and is the result of telomere shortening.
Analysis of tel1Δshu1Δ SC
The yeast Shu complex, consisting of Shu1/Shu2/Psy3/Csm2, facilitates efficient HR repair of specific replication-induced lesions (Shor ). Shu1 and Psy3 share homologies with the humanRAD51 paralogs XRCC2 and RAD51D, respectively, which act as recombination mediators (Martin ). It was proposed that the Shu complex stabilizes Rad51 filaments to promote HR by inhibiting the disassembly reaction of the anti-recombinase Srs2 (Bernstein ). To determine if the SC was due to increased anti-recombinase activity of Srs2, we deleted in and double mutants. In the triple mutant, deletion of did not rescue the SC to CPT of or suggesting that the SC interactions are Srs2 independent (Figure 5C).
SC interactions are conserved in Caenorhabditis elegans
To determine if the SC interactions observed in yeast were conserved in a more complex organism, we tested two classes of the SC interactions found in yeast in the model metazoan, C. elegans. We first tested whether the SC to CPT of was conserved. C. elegans has a well-characterized Tel1 ortholog, ATM-1 (Garcia-Muse and Boulton. 2005; Bailly ; Jones ), and a single RAD-51 paralog, RFS-1, which when mutated shares many phenotypic similarities with yeast Shu complex mutants (Ward ; Yanowitz 2008; Ward ). We constructed ; double mutants and measured the sensitivity of the single and double mutants to a range of CPT concentrations. Even in the absence of CPT, ; double mutants exhibited a reduced brood size and increased frequency of arrested embryos and male progeny, which are phenotypes consistent with increased chromosome instability (Figure 6, A and B, and Figure 7A). The ; double mutant was significantly more sensitive to CPT than either single mutant, demonstrating that the SC interaction was conserved between yeast and C. elegans (Figure 6C).
Figure 6
Conservation of SC interactions in C. elegans. (A) The percentage of embryonic arrest observed in the progeny of single and double mutants. Error bars represent SEM of at least eight broods. (B) The total brood size (arrested embryos and viable progeny) of self-fertilizing single and double mutants. Error bars represent SEM of at least eight broods. (C) Comparison of sensitivity to CPT of atm-1, rfs-1, and atm-1; rfs-1 mutants. Error bars represent SEM of at least three experiments. (*) P < 1E-2, (**) P < 1E-3, (***) P < 1E-4, (****) P < 1E-5 by Student’s t-test analysis. (D) Comparison of sensitivity to CPT of atm-1, cku-80, and atm-1; cku-80 mutants. Error bars represent SEM of at least three experiments.
Figure 7
Conservation of the telomere stability defect in atm-1; cku-80. (A) Frequency of males observed in progeny of single and double mutants. (*) P < 0.05, (****) P < 5E-5. Kolmogorov–Smirnov (KS) test was used instead of Student’s t-test to compare atm-1; cku-80 to the other strains because the frequency of males in different broods did not exhibit a normal distribution. (B) atm-1; cku-80 strains, although viable, exhibit reduced chromosome numbers in diakinesis oocytes, which is indicative of telomere to telomere fusions. (C) Two representative images of diakinesis nuclei of atm-1; cku-80 animals exhibiting a strong Him phenotype.
Conservation of SC interactions in C. elegans. (A) The percentage of embryonic arrest observed in the progeny of single and double mutants. Error bars represent SEM of at least eight broods. (B) The total brood size (arrested embryos and viable progeny) of self-fertilizing single and double mutants. Error bars represent SEM of at least eight broods. (C) Comparison of sensitivity to CPT of atm-1, rfs-1, and atm-1; rfs-1 mutants. Error bars represent SEM of at least three experiments. (*) P < 1E-2, (**) P < 1E-3, (***) P < 1E-4, (****) P < 1E-5 by Student’s t-test analysis. (D) Comparison of sensitivity to CPT of atm-1, cku-80, and atm-1; cku-80 mutants. Error bars represent SEM of at least three experiments.Conservation of the telomere stability defect in atm-1; cku-80. (A) Frequency of males observed in progeny of single and double mutants. (*) P < 0.05, (****) P < 5E-5. Kolmogorov–Smirnov (KS) test was used instead of Student’s t-test to compare atm-1; cku-80 to the other strains because the frequency of males in different broods did not exhibit a normal distribution. (B) atm-1; cku-80 strains, although viable, exhibit reduced chromosome numbers in diakinesis oocytes, which is indicative of telomere to telomere fusions. (C) Two representative images of diakinesis nuclei of atm-1; cku-80 animals exhibiting a strong Him phenotype.We next tested whether the SC interaction was conserved in C. elegans. There are C. elegans orthologs of Yku70, Yku80, and Dnl4, CKU-70, CKU-80, and LIG-4, respectively (Boulton ; Clejan ). We constructed ; and ; double mutants and analyzed the phenotypes. In the absence of CPT, ; and ; double mutants showed evidence of increased chromosome instability. However, their phenotypes were markedly different. ; double mutants had reduced brood sizes and increased frequency of arrested embryos but no significant increase in the frequency of males when compared to . In contrast, ; mutants exhibited similar reductions in brood size and frequency of arrested embryos compared to ; but unlike ; , ; mutants exhibited a high frequency of males, indicating increased X chromosome loss (Figure 7A). The frequency of males in ; broods varied greatly between parental animals; some animals gave large broods with few males and arrested embryos while others had severely reduced broods and very high frequencies of males and arrested embryos. As the mutant has increased frequency of apparent telomere to telomere chromosome fusions, which can result in a variable high frequency of males (Jones ), we examined whether the loss of CKU-80 was exacerbating this phenotype. DAPI staining of the ; double mutants that had very high frequencies of males revealed a reduced number of diakinetic chromosomes, suggesting that chromosome fusions had occurred. This phenotype was specific to ; and was not observed in ; (Figure 7, B and C). Although the overlapping roles of Ku proteins and ATM-1/Tel1 in maintaining telomere stability appears to be conserved between yeast and C. elegans, treatment with CPT did not further exacerbate the interaction between and in C. elegans. It is possible that the effect of CPT on C. elegans telomeres is not as great as that on yeast telomeres or because the CPT sensitivity assay we employed in C. elegans was not designed to score the effects on telomeres.
Discussion
We screened for genetic interactions with a null mutation of the yeastATM ortholog, , that increase the cytotoxicity of the topoisomerase I poison camptothecin. We identified and validated 14 synthetic cytotoxic interactions, demonstrating that applying conditions, such as genotoxic stress, increase the number of genetic interactions that can be potentially exploited for therapy.
Synthetic cytotoxic interactions with DNA-damaging agents and tel1Δ
We found no SL or strong negative genetic interactions with the null mutant, which is consistent with previous reports (Chatr-Aryamontri ). In contrast, mutations affecting the related kinase have far more SL interactions (Chatr-Aryamontri ). This observation is likely the result of Mec1 having a more prominent role than Tel1 in DDR checkpoint signaling (Morrow ; Sanchez ; Usui ). Application of DNA damage uncovers genetic interactions that are not obvious in the absence of exogenous damage. We, and others, have uncovered genetic interactions with in response to DNA-damaging agents. Two previous studies demonstrated that in combination with mutations affecting 9-1-1 DDR checkpoint genes are more sensitive to CPT, ionizing radiation, and the alkylating agent MMS (Guenole ; Piening ). Piening and co-workers also identified 10 interactions with that were specific to treatment with MMS. We observed the 9-1-1 checkpoint interactions and 2 of the 10 interactions with MMS, in response to CPT. However, several interactions were specific to CPT. These CPT-specific interactions may be due to the fact that CPT generates DNA–TopI adducts and Tel1 is hyperactivated by DNA ends that are covalently bound to proteins, such as the Top1–DNA adduct (Fukunaga ).
Classes of gene mutations that result in increased sensitivity to CPT in tel1 mutants
CK2 is composed of catalytic and regulatory subunits and is ubiquitous in eukaryotic organisms (Pinna 1990). CK2 is involved in a myriad of cellular processes, including cell growth and proliferation and is responsible for phosphorylating >300 substrates, including topoisomerase I and II (Meggio and Pinna 2003). It has been shown that CK2 protein levels and activity are increased in many cancers (Toczyski ; Ahmad ). The free catalytic subunits of CK2 also display high catalytic activity in the absence of the regulatory subunits, which appear to operate as a docking platform for binding substrates and/or substrate-directed effectors, thus modulating CK2 substrate specificity rather than its activity. Interestingly, we found that CK2 β-subunits, but not α-subunits, exhibit SC to CPT with . CK2β appears to have functions apart from the regulation of CK2 catalytic subunits (Bibby and Litchfield. 2005). In S. cerevisiae, only CK2β mutants show adaptation defects in response to DNA damage and this defect is independent of the catalytic CK2 subunits (Toczyski ). In mammalian cells, CK2β is involved in CK2-independent interactions with other proteins, including the DNA-damage checkpoint Chk1 and the HR-associated BRCA1 (O’Brien ; Guerra ). Given the apparent upregulation of CK2 in tumors, the fact that many CK2 substrates act in the DNA-damage response, and our finding that the loss of the CK2 regulatory subunit results in SC to CPT in , CK2 is a potential target for inhibition in ATM-deficient tumor cells.The spindle assembly checkpoint is made up of the evolutionarily conserved proteins Bub1, Bub3, Mad1, Mad2, and Mad3. The major trigger of the SAC is unattached kinetochores, which leads to a “wait” signal that pauses anaphase. For example, nocodazole, an inhibitor of microtubule polymerization, destabilizes the spindle and triggers the SAC (Musacchio and Salmon 2007; Lara-Gonzalez and Taylor 2012). In budding yeast, stalled replication forks also activate the SAC to block anaphase. Mutations in partially relieve replication defect-induced arrest. The loss of the replication checkpoint kinase further relieves the arrest in mutants, demonstrating the additive nature of the SAC and the replication checkpoint in response to HU-induced replication stalling (Garber and Rine 2002). In fission yeast, a similar coordination of the SAC and DNA replication checkpoint is observed in response to CPT (Collura ; Kim and Burke 2008; Chila ). The interconnectedness of the SAC and the replication checkpoint is further illustrated by the observations that ATM/Tel1 and ATR/Mec1 can directly affect spindle assembly (Kim and Burke 2008; Smith ). Our finding that loss of SAC components results in SC to CPT and HU in is consistent with the role for the SAC in responding to replication stress and the significant interplay between the SAC and the role of ATM/Tel1 in the replication checkpoint.Rrm3 is a 5′ to 3′ DNA helicase that migrates with the DNA replication fork and relieves replication fork pauses. Loss of Rrm3 results in delayed DNA replication at ∼1000 genomic loci, including tRNA genes, inactive replication origins, centromeres, rDNA repeat, and telomeres (Azvolinsky ). Although cells are DNA-repair proficient, the associated replication defects lead to Mec1/Rad53 checkpoint activation (Taylor ; Bochman ). Loss of Rrm3 also increases the probability of replication fork stalling and leads to Tel1- and Mec1-dependent phosphorylation of H2A (Szilard ). Given the role of Rrm3 in promoting replication progression and the role of Tel1 in the replication checkpoint, it is not surprising that loss of both proteins result in SC when replication is further perturbed by treatment with either CPT or HU.We found that the SC to CPT of was not due to loss of NHEJ activity but rather due to telomere defects. mutants exhibit abnormally short telomeres (Lustig and Petes 1986) and Tel1 has been shown to associate with short telomeres to promote the recruitment of telomerase and telomere elongation (Hector ; Sabourin ). Furthermore, Tel1 and Mec1 phosphorylate Cdc13 to promote telomerase activity (Tseng ). Mutations affecting the yeast orthologs of the Ku genes are also associated with shortened telomeres (Boulton and Jackson 1996; Porter ) and and double mutants have even shorter telomeres than the yku or single mutants, suggesting that the proteins act independently to promote telomere elongation (Porter ; Grandin ; Piening ). CPT treatment could further affect telomeres, as CPT has been shown to cause telomeric damage in human cells (Kang ). CPT-induced telomeric damage coupled with loss of Tel1 and Yku80/Yku70, which function independently to maintain telomere length, could result in SC to CPT through telomere damage and the inability to protect and elongate telomeres. In contrast to the SC to CPT, which was rescued by expression of a Cdc13–EST1 fusion protein, the SC of to MMS was not rescued by Cdc13–EST2 fusion protein (Piening ). The differences may be explained by the effects on the telomere of the different DNA-damage lesions caused by MMS and CPT.Although we did not observe SC to CPT in ; , we did observe a negative genetic interaction that was associated with telomere maintenance. Based on the high incidence of males and the presence of apparent chromosome fusions, it appears that loss of CKU-80 exacerbates the telomere fusions that occur in the mutant (Jones ). This was unexpected as loss of YKU-80 alone does not result in obvious telomeric defects in C. elegans (Lowden ). Our data suggest that CKU-80 may protect telomeres in the absence of ATM-1 and that the role of CKU-80 is apparent only when telomere length is compromised. It also supports the hypothesis that telomere–telomere chromosome fusions are not the product of nonhomologous end joining as ; animals have frequent apparent chromosome fusions. While we did not see an obvious SC to CPT in ; based on decreased embryonic viability, it is possible that treatment with CPT exacerbates the telomere defect in the double mutant, which would become apparent only through analysis of the F2 progeny of CPT-treated animals.We observed SC to CPT when is combined with mutations in components of the Shu complex. The Shu complex has been proposed to mediate the assembly of Rad51 filaments, which are required for HR (Ball ; Bernstein ; Sasanuma ). Recent analysis of the Shu complex demonstrated that there are two Shu subcomplexes, Psy3–Csm2, which constitutes a core subcomplex with DNA-binding activity, and Shu1–Shu2, which results in milder phenotypes than Psy3 or Csm2 when mutated (Tao ). Consistent with these data, we found that confers a more severe SC to CPT interaction with than either or . In mammalian cells, the Psy3 homolog, Rad51D, has a role in telomere maintenance (Tarsounas ) and in C. elegans, mutations in the RAD51 paralog gene, , result in telomeric repeat instability (Yanowitz 2008). However, unlike , elongation of telomeres by expression of the Cdc13–Est1 fusion protein did not suppress the SC of (data not shown). We also showed that increased Srs2 antirecombinase activity, which occurs in the absence of the Shu complex (Bernstein ), was not responsible for the SC. The fact that the SC to CPT of and Shu complex mutations was specific to CPT and not HU or the radiomimetic bleomycin suggests that the interaction could be specifically in response to aberrant DNA structures that form at CPT-stalled replication forks. Similar to Shu, which is not essential for HR-mediated DSB repair, the C. elegansRAD51 paralog, RFS-1, is required for RAD-51 foci formation after treatment with genotoxic agents that block replication fork progression such as CPT and interstrand cross-linking agents but not other replication perturbing agents, such as ionizing radiation or hydroxyurea (Ward ). Treatment of cells with CPT results in lesions that present a physical barrier to replication and cause replication fork slowing and reversal (Ray Chaudhuri ). Based on the SC to CPT when Shu mutations are combined with in yeast and when is combined with in the worm, it is possible that the Shu complex and Tel1/ATM-1 define redundant pathways that mediate the loading of Rad51 at sites of CPT-induced replication fork stalls to promote HR-mediated repair or bypass.
Synthetic cytotoxicity as a therapeutic approach
Synthetic lethal-based therapies are emerging as viable antitumor treatments (Chan and Giaccia 2011). The number of tumors susceptible to SL approaches could be increased by identifying genetic interactions that result in synthetic lethality under endogenous conditions that can affect tumors such as hypoxia, aneuploidy, or replication stress, or exogenous (applied) conditions such as increased DNA damage from DNA-damaging therapeutic agents, replication stress from replication inhibitors, or spindle destabilization by microtubule poisons. Identifying these conditional SL interactions can expand the number of tumor genotypes that can be treated by SL approaches. In this study we have defined a subclass of conditional synthetic lethality that we have called synthetic cytotoxicity. By exploiting genetic interactions, synthetic cytotoxicity has the potential to expand the spectrum of tumor genotypes that can be targeted for cytotoxic therapy and can lead to lower doses of cytotoxic therapeutic agents, which could reduce off-target cytotoxicity.
Authors: K Savitsky; A Bar-Shira; S Gilad; G Rotman; Y Ziv; L Vanagaite; D A Tagle; S Smith; T Uziel; S Sfez; M Ashkenazi; I Pecker; M Frydman; R Harnik; S R Patanjali; A Simmons; G A Clines; A Sartiel; R A Gatti; L Chessa; O Sanal; M F Lavin; N G Jaspers; A M Taylor; C F Arlett; T Miki; S M Weissman; M Lovett; F S Collins; Y Shiloh Journal: Science Date: 1995-06-23 Impact factor: 47.728
Authors: Jessica X Chong; Kati J Buckingham; Shalini N Jhangiani; Corinne Boehm; Nara Sobreira; Joshua D Smith; Tanya M Harrell; Margaret J McMillin; Wojciech Wiszniewski; Tomasz Gambin; Zeynep H Coban Akdemir; Kimberly Doheny; Alan F Scott; Dimitri Avramopoulos; Aravinda Chakravarti; Julie Hoover-Fong; Debra Mathews; P Dane Witmer; Hua Ling; Kurt Hetrick; Lee Watkins; Karynne E Patterson; Frederic Reinier; Elizabeth Blue; Donna Muzny; Martin Kircher; Kaya Bilguvar; Francesc López-Giráldez; V Reid Sutton; Holly K Tabor; Suzanne M Leal; Murat Gunel; Shrikant Mane; Richard A Gibbs; Eric Boerwinkle; Ada Hamosh; Jay Shendure; James R Lupski; Richard P Lifton; David Valle; Deborah A Nickerson; Michael J Bamshad Journal: Am J Hum Genet Date: 2015-07-09 Impact factor: 11.025
Authors: Romulo Segovia; Yaoqing Shen; Scott A Lujan; Steven J M Jones; Peter C Stirling Journal: Proc Natl Acad Sci U S A Date: 2017-02-21 Impact factor: 11.205
Authors: Souhrid Mukherjee; Joy D Cogan; John H Newman; John A Phillips; Rizwan Hamid; Jens Meiler; John A Capra Journal: Am J Hum Genet Date: 2021-09-15 Impact factor: 11.025
Authors: Carolyn A Turcotte; Solomon A Sloat; Julia A Rigothi; Erika Rosenkranse; Alexandra L Northrup; Nicolas P Andrews; Paula M Checchi Journal: Genetics Date: 2018-01-16 Impact factor: 4.562