High-throughput screening of the Saccharomyces cerevisiae genome for 2-amino-3-methylimidazo [4,5-f] quinoline resistance identifies colon cancer-associated genes

Abstract Heterocyclic aromatic amines (HAAs) are potent carcinogenic agents found in charred meats and cigarette smoke. However, few eukaryotic resistance genes have been identified. We used Saccharomyces cerevisiae (budding yeast) to identify genes that confer resistance to 2-amino-3-methylimidazo[4,5-f] quinoline (IQ). CYP1A2 and NAT2 activate IQ to become a mutagenic nitrenium compound. Deletion libraries expressing human CYP1A2 and NAT2 or no human genes were exposed to either 400 or 800 µM IQ for 5 or 10 generations. DNA barcodes were sequenced using the Illumina HiSeq 2500 platform and statistical significance was determined for exactly matched barcodes. We identified 424 ORFs, including 337 genes of known function, in duplicate screens of the “humanized” collection for IQ resistance; resistance was further validated for a select group of 51 genes by growth curves, competitive growth, or trypan blue assays. Screens of the library not expressing human genes identified 143 ORFs conferring resistance to IQ per se. Ribosomal protein and protein modification genes were identified as IQ resistance genes in both the original and “humanized” libraries, while nitrogen metabolism, DNA repair, and growth control genes were also prominent in the “humanized” library. Protein complexes identified included the casein kinase 2 (CK2) and histone chaperone (HIR) complex. Among DNA Repair and checkpoint genes, we identified those that function in postreplication repair (RAD18, UBC13, REV7), base excision repair (NTG1), and checkpoint signaling (CHK1, PSY2). These studies underscore the role of ribosomal protein genes in conferring IQ resistance, and illuminate DNA repair pathways for conferring resistance to activated IQ.


Introduction
The heterocyclic aromatic amine (HAA) 2-amino-3methylimidazo[4,5-f ] quinoline (IQ) is a mutagen (Kasai et al. 1980) and a class 2A carcinogen (IARC 1993;Olalekan Adeyeye and Ashaolu 2021), present in charred meat and cigarette smoke.IQ and other heterocyclic aromatic amines (HAAs) are generated by Maillard reactions involving creatine, sugar, and an amino acid (Jägerstad et al. 1991); these reactions occur when charring meat and may suggest why red meat consumption increases colon cancer risk (Sinha et al. 1999;Sinha et al. 2001;Linseisen et al. 2002;Aykan 2015;Gongora et al. 2019).While charred meats contain a mixture of HAAs, IQ exposure has been directly correlated to inducing cancer.For example, IQ exposure induces liver, kidney, and colon cancers in rats (Turesky and Markovic 1995).Although the exact mechanism by which IQ induces cancer is unknown, IQ exposure is associated with oxidative stress (Carvalho et al. 2015) and higher levels of DNA damage (de Carvalho et al. 2016).Genotoxic endpoints of IQ exposure include higher frequencies of both mutations (Aeschbacher and Turesky 1991) and sister chromatid exchange (Sawada et al. 1994).These studies suggest that IQ acts as a genotoxin contributing to the carcinogenicity associated with red meat or cigarette smoke.
IQ pre-incubated with rat liver S9 is more mutagenic in the Ames assay than IQ alone (Kasai et al. 1980), indicating that IQ's genotoxicity requires bioactivation.In humans, bioactivation occurs by a 2-step mechanism (Fig. 1; Kim and Guengerich 2005).First, cytochrome P450 enzymes (CYPs), such as CYP1A2, hydroxylate IQ to form N-hydroxy-IQ, and second, N-acetyltransferases (NATs), such as NAT1 and NAT2, acetylate the N-hydroxy-IQ to form N-acetoxy-IQ (Fig. 1).This latter compound is unstable and undergoes spontaneous heterolytic cleavage to generate a highly reactive nitrenium ion (Kato 1986;Turesky and Vouros 2004;Kim and Guengerich 2005), which forms C 8 and N 2 Gua adducts.Interestingly, the N 2 Gua adduct appears refractory to nucleotide excision repair (NER) and is bypassed by error-prone polymerase (Stavros et al. 2014), while the C 8 adduct is more efficiently removed (Turesky et al. 1996).NAT2 expression increases the mutagenicity of IQ in both CYP1A2-expressing strains of Salmonella typhimurium and in Chinese hamster ovary cells (CHO) (Grant et al. 1992;Josephy et al. 1995) and the recombinogenicity of IQ in CYP1A2-expressing strains of Saccharomyces cerevisiae (budding yeast, Paladino et al 1999).These studies indicate that expression of both cytochrome P450 and NAT enzymes are required for high levels of IQ genotoxicity.
Genetic factors that increase the mutagenicity of HAAs have been correlated to higher incidence of smoking and diet-associated cancers.For example, the human NAT2*4 rapid acetylator allele, whose expression in CHO cells increases HAA-associated mutations more than that of NAT2 slow acetylator alleles (Probst-Hensch et al. 1995;Grant et al. 1997;Fretland et al. 2001), is a colon cancer risk factor among Caucasians (Butler et al. 2008) and smokers (Slattery et al. 1998;Nöthlings et al. 2009).Individuals who frequently consume red meat and exhibit rapid acetylator NAT phenotypes, carrying NAT2*4 or NAT1*10 alleles, are at a higher risk for developing colon cancer (Lilla et al. 2006;Metry et al. 2010;Wang et al. 2011;Le Marchand 2021).Although HAAs are found in varying concentrations in the diet and cigarette smoke, the epidemiological studies suggest that HAA genotoxicity aggravates cancer risk.
Besides DNA adducts, IQ exposure also generates DNA strand breaks and oxidative stress, suggesting that DNA repair genes may function in suppressing IQ-associated genotoxicity (Pfau et al. 1999;Pezdirc et al. 2013).The importance of this question is underscored by observations that high-risk factors for colon cancer include biallelic mutations of 2 base excision repair genes, NTHL1 (Weren et al. 2015;Te Paske et al. 2020), and hMYH (Al-Tassan et al 2002; Wooden et al. 2004).Low penetrant colon cancer risk factors include polymorphisms in the DNA repair genes RAD18, ERCC5, XPC, PARP, and APE1 (Pan et al. 2012;Aggarwal et al. 2017;Matejcic et al. 2021).Whether DNA repair gene polymorphisms found specifically in colon cancer confer higher levels of HAA genotoxicity is unknown.
Budding yeast is a model organism for profiling genes that confer resistance to genotoxins.Over 30% of the yeast genome is orthologous to humans (O'Brien et al. 2005), and human orthologs can functionally substitute for 469 essential yeast genes (Kachroo et al. 2015).Each nonessential gene has been individually knocked out and tagged with both upstream and downstream barcodes (uptags and downtags, Giaever et al. 2002); this deletion library has already been shown to be a powerful reagent in identifying drug targets and toxicant resistance genes (Jo et al. 2009;Giaever and Nislow 2014;De La Rosa et al. 2017).While budding yeast does not bioactivate HAAs, expression of human CYP1A2 and NAT2 increases the HAA-associated translocation frequencies in budding yeast (Paladino et al. 1999).A "humanized" yeast deletion library expressing CYP1A2 and NAT2 has already been used to profile the yeast genome for aflatoxin B 1 (AFB 1 ) resistance (St. John et al. 2020).In this study, we profiled the same "humanized" deletion collection for resistance to IQ.We chose IQ because activated IQ is the most recombinogenic HAA in yeast (Paladino et al. 1999) and is highly mutagenic in strains of S. typhimurium (Knasmüller et al. 1999;Nohmi and Watanabe 2021).We identified budding yeast genes that confer IQ resistance; mammalian orthologues of a subset of the IQ resistance genes are associated risk factors for colon cancer.

Media and chemicals
Standard media were used for the culture of yeast and bacterial strains (Burke et al. 2000).The bacterial DH1 strains containing the vectors pCS316 and pCYP1A2_NAT2 were cultured in LB-Amp (Luria Broth media with 100 µg/mL ampicillin).Media used for the culture of yeast cells included YPD (yeast extract, peptone, dextrose), SC (synthetic complete, dextrose), and SC-URA (SC lacking uracil).A mouse monoclonal antibody to CYP1A2 (ab22717) and a mouse polyclonal antibody to NAT2 (ab88443) were purchased from Abcam.Purified recombinant NAT2 protein was purchased from OriGene (TP761755).A goat antimouse IgG-HRP antibody was purchased from Santa Cruz Biotechnology (sc-2005).2-amino-3-methylimidazo [4,5-f ] quinoline (IQ) was purchased from Toronto Research Chemicals.Dimethyl sulfoxide (DMSO), methanol, and glacial acetic acid were purchased from Sigma.AFB1 was dissolved in DMSO, while IQ was dissolved in methanol containing 0.1% acetic acid, and is simply referred to as MeOH.

Strains and plasmids
Yeast strains used in this study were derived from either BY4741, BY4743 (Brachmann et al. 1998), YA289, or YB198 (Fasullo et al. 2001) and are of the S288C background.Strain genotypes are listed in Supplementary Table 1.The yeast deletion collection is derived from BY4743, whose genotype is MATa/α his3Δ1/his3Δ1 leu2Δ0/ leu2Δ0 LYS2/lys2Δ0 met15Δ0/MET15 ura3Δ0/ura3 Δ0.The diploid and haploid homozygous deletion libraries were purchased from Open Biosystems and are now available from Horizon Discovery (https://horizondiscovery.com/en/non-mammalian-research-tools/ products/yeast-knockout).The pooled diploid homozygous deletion library (N = 4607) was a gift of C. Vulpe (University of Florida).The CYP1A2 expression plasmid, pCS316 (Sengstag et al. 1996), was obtained by CsCl centrifugation (Ausubel 2003).The plasmid used in this study, pCYP1A2_NAT2 (St. John et al. 2020), was constructed by removing the human oxidoreductase sequence from pCS316 and replacing it with the Not1 fragment containing human NAT2 from pRS424-CYP1A2_NAT2.The cDNA of NAT2 was previously positioned in the pGAC.24(Lesser and Guthrie 1993) expression plasmid under the control of the glycerol phosphate dehydrogenase (GPD) promoter.

Transformation protocol
The modified diploid homozygous deletion libraries were prepared as previously described (St. John et al. 2020).In brief, 5-fluoroorotic acid resistant isolates were taken from each strain in the diploid homozygous deletion collection and inoculated in 100 μL of YPD in 96-well plates.After overnight incubation at 30°C, plates were centrifuged, washed, and resuspended in 1-step buffer (50% PEG, MW 3300, 0.2 N lithium acetate, 100 mM dithiothreitol, 500 mg/mL denatured salmon sperm DNA).Each well received 1 µg of pCYP1A2_NAT2 and plates were incubated for 30 mins at 30°C, before pipetting 10 μL onto duplicate SC-URA plates.Two Ura + transformants were chosen from each strain and frozen in SC-URA with 7.5% DMSO in 96-well plates.This strategy was successful in transforming approximately 90% of the diploid homozygous deletion library.
To determine the retention of pCYP1A2_NAT2 in wild type and diploid mutants, Ura + transformants were selected in SC-URA and then diluted, approximately 10 5 cells were exposed to no treatment, 2% MeOH, 400 µM IQ or 800 µM IQ, and incubated for 20 hours.An appropriate dilution was then plated on YPD to obtain single colonies.URA + colonies were counted after replica plating onto SC-URA plates and a 2-day incubation.

Functional profiling of the yeast genome
We profiled the yeast genome for resistance to the HAA IQ using 2 different pooled collections of the diploid nonessential deletion strains, derived from BY4743.The first pool was a yeast deletion collection library of 4607 knock-out strains, each strain containing a single deletion in a nonessential gene (Jo et al. 2009).The "humanized" pool expressed CYP1A2 and NAT2 in approximately 4,900 strains derived from the 5,417 diploid strains (St. John et al. 2020).
Tubes of the pooled library (10 plates per tube) were thawed on ice and 250 µL of each tube's content were added to a 15 mL tube and allowed to recover by incubation at 30°C for 2 hours in a shaking incubator.From the recovered pooled libraries, 100 μL was inoculated into culture tubes containing 2 mL of either YPD for the original yeast library or SC-URA for the yeast deletion collection containing the pCYP1A2_NAT2 plasmid.Quadruplicate samples from each library were then treated with 2% MeOH, 400 μM IQ, or 800 μM IQ.Samples were incubated on a rotary incubator at 30°C for either 15 hours (for the YPD cultures) or 20 hours (for the SC-URA cultures), approximately 10 generations.The cells were then centrifuged and washed thrice with sterile water.
To amplify molecular barcodes, DNA from each treatment tube was isolated using the "smash and grab" method (Hoffman and Winston 1987) and resuspended in TE buffer (10 mM Tris-HCl, 1 mM EDTA, pH 8.0) at concentrations of 2-4 µg/µL.Polymerase chain reaction (PCR) was performed on each sample using customized oligos (Supplementary Table 2) that hybridize to the regions just upstream and downstream of the barcode "uptag" sequences; the oligos contained a multiplex index and sequencing primers for use with the Illumina Hi-Seq 2500 (Smith et al. 2009;Robinson et al. 2014).PCR protocols for barcode amplification were previously described (Robinson, et al. 2014).The ∼150 bp PCR products were resolved on a 10% polyacrylamide nondenaturing gel, extracted using the "crush and soak" method (Green and Sambrook 2019), and quantified.A pooled sample containing 50 ng of DNA from each of the 32 samples was sent to the Bauer Sequencing Core Facility at Harvard University for quality assurance and sequencing using the Illumina HiSeq 2500 platform.For experiments comparing 5-generation vs 10-generation IQ exposures, samples were sent to the University at Buffalo Genomics and Bioinformatics Core Facility for quality assurance and sequencing on the Illumina HiSeq 2000 platform.

Barcode analysis
The sequencing data obtained in the fast all with quality (FASTQ) files was demultiplexed to separate reads from different exposures.The individual files were then subjected to a FastQC report, and then trimmed to include only the barcode region with 2-3 extra bases on either side.Barcodes were counted with an in-house program and associated with strains.Replicate treatments were grouped and compared to solvent treatment.Comparisons of treatment were performed using Bioconductor's "TCC" package in RStudio (Sun et al. 2013), processed in batches using RStudio version 1.4.1717 and R version 4.1.0using the code in Supplementary File S1. To

Western blots
Expression of CYP1A2 and NAT2 in yeast cells was determined by Western blot.Cells were inoculated in YPD (Original) or SC-URA medium (Containing pCYP1A2_NAT2) and grown to log phase (A 600 between 0.5 and 1) and then concentrated.Protein extracts were prepared as described previously by Foiani et al. (1994).

CYP1A2 and NAT2 enzyme activity assays
To determine whether the expressed protein was functional, the methoxyresorufin O-demethylase (MROD) activity of CYP1A2 was measured, using a method previously described (Fasullo et al. 2014).Microsomes were prepared using a modified protocol originally described by Pompon et al. (1996) and St. John et al. (2020).Microsomal protein concentrations were 15-20 mg/ml.The reaction mix consisted of adding 1 μL microsomes to 5 μM 7-methoxyresorufin in 10 mM Tris, pH 7.4 at 37°C.Reactions were started with the addition of 500 mM NADPH and fluorescence was recorded every minute for 1 hour.Fluorescence of resorufin was measured on a Tecan Infinite M200 plate reader with an excitation wavelength of 535 nm and emission wavelength of 580 nm.Microsomes derived from BY4743 containing no CYP1A2 were used for the negative control.For a positive control, rat liver S9 fractions were used.
To measure the NAT activity of NAT2, we used cytosolic extracts from yeast expressing CYP1A2 and NAT2, prepared as described by Paladino et al. (1999).Briefly, a 100 mL culture, grown to midlog phase, was pelleted and cells were washed twice before being resuspended in 2 mL ice-cold disruption buffer (100 mM Tris, pH 7.4, 100 mM DTT, 1 mM EDTA, 10% (v/v) glycerol, and 1 × yeast proteinase inhibitor cocktail IV (Thermo Scientific).Cells were lysed in the bullet blender after addition of 200 mg 0.45-0.55mm glass beads.After 10 mins on the bullet blender, the lysate was removed to a new tube and the beads were washed twice in an equal volume of ice-cold disruption buffer.The 3 lysate fractions were pooled in a new tube.Cell debris was pelleted by centrifugation at 21,000 × g for 20 mins at 4°C and the supernatant was further centrifuged for 1 hour at 100,000 × g, before being aliquoted for assays.Extract (50 μL) was then added to a new tube containing 27.5 μL reaction buffer (225 mM Tris, pH 7.5, 4.5 mM DTT, 4.5 mM EDTA, 1 × proteinase inhibitor cocktail IV), 4.5 μL of 2 mM sulfamethazine (SMZ), and reactions were initiated with the addition of 18 μL of 2.5 mM acetyl-CoA.Negative controls included water in place of acetyl-CoA.Reactions were terminated after 30 minutes by the addition of 50 μL 25% trichloroacetic acid.To each reaction tube, 760 μL of 5% 4-dimethylamino-benzaldehyde in acetonitrile was added and mixed.Tubes were spun for 1 minute at 2,000 × g to remove precipitate.Aliquots of 100 μL were then added to a 96-well plate and absorbances at 450 nm were measured.The presence of acetyl-SMZ was determined by loss of signal at 450 nm.

Growth curves
In brief, individual saturated cultures were prepared for each yeast strain.Cell density was adjusted to 0.8 × 10 7 cells/mL for all cultures.For single mutants, strains containing pCYP1A2_NAT2 were maintained in selective medium (SC-URA), and exposed in duplicate MeOH, 400 μM IQ, or 800 μM IQ.The plate was covered with optical adhesive tape and incubated at 30°C in a TECAN Infinite M200 plate reader and absorbance at 600 nm was recorded every 10 minutes for 24 hours, as previously described (Fasullo et al. 2010;Fasullo et al. 2014;St. John et al. 2020).Absorbance was plotted against time and the areas under the curve (AUCs) were calculated for the time interval of 0 to 20 hours using freely available software (https://www.padowan.dk/download/).We calculated percent growth in the presence of toxin using the formula (AUC IQ /AUC MeOH )×100% (Toussaint et al. 2006).Statistical significance of differences between growth percentages for diploid strains and BY4743 was determined by Student's t-test, assuming constant variance between samples.

Competitive growth assays
We measured competitive fitness of a strain using a modified protocol (North et al. 2012;De la Rosa et al. 2017).In brief, the strain was grown to a stationary phase in either YPD or SC-URA (for strains containing pCYP1A2_NAT2).In addition to the strain of interest, a wild-type strain (YB676), constitutively expressing GFP containing the pCYP1A2_NAT2 plasmid, was grown to saturation.Cells from the strain of interest and the GFP-expressing strain were combined at a set ratio and inoculated into 2 mL of culture at a concentration of 3× 10 5 cells/mL.These cultures were then treated in duplicate with either solvent, 400 μM IQ, or 800 μM IQ and incubated for 20 hours on a rotary shaker at 30°C.Cells were then washed thrice with water and resuspended in 1 mL of water.Cells (2 × 10 4 ) were then counted on an Amnis ImageStream X imaging flow cytometer (N = 2).Cells were imaged with brightfield and a 488 nm laser at 50 mW power.The percentage of GFP-expressing cells to a total number of cells was compared for solvent-exposed and IQ-exposed cells.Data were analyzed with the IDEAS v6.2 software.In-focus cells were gated, followed by single cells, which included budded cells.GFP fluorescence intensity was quantified for the focused, single-cell population which resulted in 2 distinct peaks, 1 for fluorescent cells and 1 for nonfluorescent cells.Differences in GFP percent between IQ and methanol treatments of the strains of interest were compared to the difference in GFP percent of the IQ and methanol treatments with BY4743.The statistical significance of these differences was determined by Student's t-test, assuming constant variance between samples.

Ames test
To ensure the mutagenicity of IQ, we performed the Ames test with the preincubation protocol (Maron and Ames 1983).In brief, cultures of Salmonella typhimurium TA98 and TA100 were grown overnight at 37°C to a concentration of 1-2 × 10 8 cells/mL.The cultures were placed on ice until preincubation.An S9 master mix was made according to the protocol (33 mM KCl, 8 mM MgCl 2 , 5 mM glucose-6-phosphate, 4 mM NADP, 100 mM phosphate buffer pH 7.4), and to it was added rat liver S9 (4%) or yeast-derived microsomes, and 495 μL of this master mix were added to culture tubes, in triplicate, along with 100 μL of the chilled bacteria.IQ was added at either 4 or 40 µM concentration before incubating the tubes in a heat block at 37°C for 20 minutes.After incubation, 2 mL of molten 0.6% top agar, heated to 48°C, was added and the mix was pipetted onto minimal glucose agar plates.Plates were cooled for 5 minutes, then inverted and incubated for 48 hours at 37°C before His + colonies were counted.

Trypan blue staining
To measure cell viability after IQ exposure, we performed a trypan blue exclusion assay.Selected strains expressing CYP1A2 were inoculated in SC-URA until cultures reached an A 600 of 0.1-0.5, and then exposed to either 800 μM IQ or 2% solvent alone.After incubating for 3 hours, cells were washed twice in sterile phosphatebuffered saline (PBS) and stained with trypan blue at a final concentration of 2 mg/ml (Liesche et al. 2015).Both the total number and blue-stained (dead) cells were counted on a light microscope using a hemocytometer.At least 200 cells were counted for each experiment, and percentages of live and dead cells were calculated, N = 3.The differences between percent dead cells after IQ treatment and methanol treatment were calculated.To determine statistical significance, the differences in percent dead for each strain of interest were compared to the difference in percent dead in wild type (WT, YB679, N = 10) using Student's 1-tailed t-test.

Expression of CYP1A2 and NAT2 enhances IQ toxicity in yeast rad + and rad4 rad51 mutants
To confirm that the expression of CYP1A2 and NAT2 bioactivate IQ, the expression plasmid pCYP1A2_NAT2 was introduced in wild type yeast (BY4743) and a rad4 rad51 strain by selecting Ura + transformants.The rad4 rad51 strain is defective in both nucleotide excision and recombinational repair and is sensitive to DNA damaging agents that form bulky DNA adducts as well as DNA strand breaks.Expression of CYP1A2 and NAT2 proteins was then confirmed by Western blots (Fig. 2a and 2b, respectively).We performed MROD assays to quantify CYP1A2 activity (Fig. 2c, Fasullo et al. 2014) and NAT assays to detect NAT2 activity (Fig. 2d).Yeast extracts from these assays also activated IQ in the Ames assay; compared to extracts from the strain expressing no human genes, the mutagenicity of IQ was highest (P < 0.05, Supplementary Table 3) when IQ was pre-incubated using microsomal and cytosolic extracts obtained from strains expressing CYP1A2 and NAT2, but still detectable using extract from strains expressing CYP1A2 alone (Supplementary Table 3).These data indicate that CYP1A2 and NAT2 were expressed in yeast and were sufficient to increase the mutagenicity of IQ.
By calculating AUCs, we determined IQ concentrations that reduced growth in both the wild-type Rad + diploid (BY4743) and the rad4 rad51 haploid mutant, which either did ("humanized") or did not (original) express CYP1A2 and NAT2.Exposure to 400 μM IQ and 800 μM IQ, reduced growth of the "humanized" Rad + diploid by 6 and 17%, respectively (Fig. 2e), and in the "humanized" rad4 rad51 haploid by 40 and 59%, respectively (Fig. 2f).On the other hand, exposure to 400 μM IQ and 800 μM IQ, reduced the growth of the original wild type by 1 and 7%, respectively (Fig. 2g), and in the original rad4 rad51 haploid by 6 and 13% (Fig. 2h).In comparison to the Rad + diploid (BY4743), the higher level of IQ toxicity in the "humanized" Rad + diploid is significant only for the 800 µM IQ exposure (N = 8, P = 0.035).In comparison to rad4 rad51, the higher level of IQ toxicity in the "humanized" rad4 rad51 is significantly different both at 400 µM IQ (N = 3, P = 0.019) and 800 µM IQ (P = 0.007) exposure.These data indicate the "humanized" strains exhibit a higher level of IQ genotoxicity.

Identification of resistance genes by barcode analysis
We identified genes that confer resistance to either IQ per se, or activated IQ in 6 independent screens (Fig. 3).We exposed the original pooled diploid collection to 400 μM and 800 μM IQ for 10 generations (10G).We also exposed the pooled "humanized" yeast diploid collection to 400 μM and 800 μM IQ for 10G.In addition, we exposed the "humanized" collection to 800 μM IQ exposure for 5 generations (5G) and 10G.The control for all screens was exposed to 2% methanol with 0.1% acetic acid (MeOH) and 4 biological replicates were performed for each exposure and 2 technical replicates were sequenced.Approximately 98% (4513/4607) of strains from the original library and 89% (4343/4900) of strains from the "humanized" library were detected in the pooled cultures.Genes conferring IQ resistance were identified by differential abundance calculations of barcode counts where m-values <0 and q-values <0.1.
Across 4 screens of the "humanized" library, 1159 ORFs from the "humanized" library were identified that conferred IQ resistance at any condition (see Supplementary Table 4), of which, 424 ORFs were identified in at least 2 independent screens.In the time-dependent study, 12/21 (57%) of the ORFs identified at 5G were present in the 10G (113 ORFs) dataset (Fig. 3, Panel A).In the concentration-dependent study, 412/671 (62%) ORFs from the 400 µM screen were also present in the 800 µM dataset (792 ORFs) (Fig. 3, Panels b and c).These data thus show considerable overlap between screens of the "humanized" library.
A random selection of the concentration-dependent screen was used to estimate false positives and false negatives.To estimate the number of false positives among the overlapping set of Panel E is a growth curve of the Rad + diploid strain containing pCYP1A2_NAT2 (YB679), N = 8.Percent growth with 400 µM IQ and 800 µM IQ were 94 and 83%, respectively.Panel F is a growth curve of the rad4 rad51 haploid yeast strain containing pCYP1A2_NAT2 (YB400), N = 3. Percent growth with 400 µM IQ and 800 µM IQ were 60 and 41%, respectively.Panel G is a growth curve of the Rad + diploid strain (BY4743), N = 8.Percent growth with 400 µM IQ and 800 µM IQ were 99 and 93%, respectively.Panel H is a growth curve the rad4 rad51 (YB226) strains, N = 3. Percent growth with 400 µM IQ and 800 µM IQ were 94 and 87%, respectively.study, ORFs were alphabetized and 40 were randomly selected.Of these 40, 39 (98%) strains exhibited sensitivity to IQ in the trypan blue assay (Supplementary Table 5).Strains not identified as significant (2,049 strains) were also alphabetized and 10 strains were randomly selected; of these, 9 were not significantly IQ sensitive (Supplementary Table 5).We thus estimate that our false positive percentage was less than 5%, while our failure to detect IQ-sensitive strains was approximately 10%.

ORFs conferring resistance in the "humanized" library
Based on our estimate that over 95% of duplicated ORFs confer IQ resistance in the "humanized" library, we identified key interactions that mediate IQ resistance from our list of 424 ORFs (Fig. 4a), using ELIXIR's Core Data Resource STRING (https:// string-db.org/,Szklarczyk et al. 2019 ).A network of proteins from ORFs conferring IQ resistance in at least 1 screen is presented in Supplementary Fig. 1.In both sets, the largest and most heavily connected cluster was ribosomal proteins (Fig. 4a, Supplementary Fig. 1).Additional clusters in both sets are encoded by DNA damage tolerance genes, such as REV7, UBC13, and RAD18, and nitrogen metabolism genes, such ALT1, ADE6, and LYS20.These data thus illustrate several key gene interaction networks involved in protein metabolism and DNA damage repair.
Genes selected for validation included those from the major interactome clusters, genes that were missing from the library but whose presence was implied based on related significant genes, as well as genes with ambiguity in whether they conferred resistance or sensitivity based on the outcomes of multiple screens.We validated IQ resistance genes using growth curves, competitive growth assays, and trypan blue staining; while growth curves and competitive growth assays measure toxin resistance after chronic exposures, trypan blue was used to determine cell viability after an acute 3 h IQ exposure.The interactome of proteins from all validated ORFs (Fig. 4b) contains a larger cluster of DNA repair genes than the interactome of proteins from all ORFs appearing at least twice (Fig. 4b).
In total, we tested 51 selected genes by both competitive growth and trypan blue assays, of which, 18 were also tested by growth curves (Table 1).Twenty-one genes were validated by only 1 method, 30 (59%) were validated by at least 2 methods, and 7 were validated by all 3. We previously observed that after 10 generations, 10-40% of total cells, depending on the strain, still lose the pCYP1A2_NAT2 plasmid, even in selective media, rendering it difficult to accurately determine toxin sensitivity by trypan blue viability or growth curves (St. John et al. 2020).On the other hand, CYP1A2 expression increases the number of nonfluorescent cells in the "wild-type" strain containing the GFP construct ([6 ± 4] %, N = 28) thus undercounting the true number of "wild-type" cells after competitive growth.This may account for why over half the validated genes (8 out of 15) that have m-values greater than −0.7 did not validate by competitive growth assays.Nonetheless, despite these limitations, most genes are validated by multiple independent methods for monitoring toxin sensitivity.
As an example, mutants representative of the ribosome (rpl33b), nitrogen metabolism (alt1), and DNA repair (sae2) clusters, exhibited particularly high IQ sensitivity either by growth curves or competitive growth assays, as well as validating with trypan blue (Fig. 5).ALT1 functions in utilizing alanine as a nitrogen source, and alt1 mutants are sensitive to oxidative stress.RPL33B encodes the ribosomal 60S genes and is essential when its paralog RPL33A is defective.SAE2 is an endonuclease involved in several DNA repair and maintenance functions.Only significant ORFs conferring resistance in at least 1 screen and which were validated by at least 1 method are included in further analysis (82 ORFs).
We identified protein complexes using the Princeton gene ontology (GO) Term finder (https://go.princeton.edu/cgi-bin/LAGO, Boyle et al. 2004).Protein complexes in global cell regulation and housekeeping were identified among the validated ORFs using the Princeton GO Term finder (P < 0.05), and redundancy was reduced using Revigo (Supek et al. 2011).The resulting 5 components GO groups included the nonmembrane bound organelle, ribosome, ribosomal subunit, the protein kinase CK2 complex, and the HIR complex (Table 2).The CK2 kinase complex regulates cell cycle progression, while the HIR complex, also referred to as the replication-independent HIR complex, regulates chromatin structure and affects gene expression.These complexes underscore the importance of regulating cell cycle progression and protein translation and gene expression in IQ resistance.For all GO Terms from the complete set of overlapping ORFs, see Supplementary Table 6.
In addition to the protein complexes, we identified process and function GO groups represented by the 82 ORFs (P < 0.01) using the Princeton LAGO tool and reduced redundancy using Revigo, resulting in 4 "functions" and 16 "processes" (Table 3).Key constituents that were expected included the structural constituent of the ribosome, DNA repair, DNA damage response, and nitrogen cycle metabolic process GO terms.Other terms included alpha-amino acid metabolic process, regulation of nucleobase-containing compound metabolic process, DNA binding, regulation of ergosterol biosynthetic process, and regulation of transcription by RNA polymerase III.These process terms further emphasize the role of nitrogen metabolism in conferring IQ resistance, as well as the role of nucleic acid binding.
Fourteen IQ-validated resistance genes that function in either DNA replication checkpoints, DNA repair, or DNA damageassociated gene expression (Cherry et al. 2012) are listed in Table 4.While NER genes were not prominent, we identified the dual function base excision repair gene, NTG1.Several genes are directly associated with DNA damage tolerance, such as UBC13, REV7, and RAD18, while others are involved in cell cycle progression in the presence of DNA damage, including PSY2, CKA1, and CKB1.ADE6, which functions in de novo purine synthesis, suppresses chromosome loss, as revealed in a MATa-like faker screen (ALF, Yuen et al. 2007).Recombinational repair genes include SAE2.Although RAD52, MRE11, and XRS2 were missing in all the screens of the "humanized" library, the rad52, mre11, xrs2 strains were subsequently observed to be IQ sensitive by trypan blue assay.MPH1 functions in blocking recombination and facilitating replication (Prakash et al. 2009).CHK1 is a checkpoint signaling genes downstream of RAD9 (For review, see Chen and Sanchez 2004); while RAD9 was not initially identified in the screen, the rad9 mutant is IQ sensitive by both trypan blue and competitive growth assays.These genes emphasize the role of postreplication repair in conferring IQ resistance.

ORFs conferring resistance in the original yeast library
To distinguish between ORFs that confer resistance to activated IQ or IQ per se, we screened the original yeast deletion library, lacking Italicized values did not meet the threshold of q < 0.1 in 800 µM study.m-values and q-values from 10G screen, from which these genes were identified.human CYP1A2 and NAT2, for IQ resistance.We identified 835 ORFs conferring resistance to either 400 μM or 800 μM IQ, of which 143 ORFs were identified in both screens (Fig. 3, Supplementary Table 7).In the interactome of overlapping ORFs in the original library, the prominent clusters are related to asparagine-linked glycosylation, and the Rpd3L histone deacetylase (Supplementary Fig. 2a).The over-represented GO terms in the screens of the original library were negative regulation of nucleobase-containing compound metabolic process (GO:0045934) and cellular protein modification process (GO:0006464).The protein interactome of all 835 ORFs further expands the clusters of protein modification process and glycosylation, and clearly denotes ribosomal proteins, vacuolar ATPase, and DNA damage repair clusters (Supplementary Fig. 2b).These ORFs indicate protein synthesis and modification are prominent in conferring IQ resistance.
There were 254 ORFs in at least 1 screen of the "humanized" library which overlapped with at least 1 screen of the original library (Supplementary Tables 4 and 7).Of these, notable GO terms were negative regulation of DNA recombination, nucleotide-excision repair, DNA damage recognition, lipid binding, and vacuole (Supplementary Table 8).The m-values of 211 ORFs (83%) were lower from the "humanized" library than the original library (Supplementary File S2).The percentage was slightly greater (88%, 14/16) when comparing DNA damage repair genes.These data suggest there is a subset of genes required for conferring resistance to both IQ per se and "metabolically activated" IQ.
When comparing the interactomes of the unique ORFs identified in the "humanized" (Supplementary Fig. 1) and original (Supplementary Fig. 2b) libraries, the ribosome cluster is prominent in both; 15 genes were identified in the original library, while 38 were identified "humanized" library.Four ribosomal genes, RPL22A, RPL26A, RPL37B, and RPL6B, were identified in profiling both the "humanized" yeast library and the original library.These data suggest that activated IQ may target different subunits of the ribosome than IQ exposure per se.

ORFs conferring sensitivity to activated IQ
We also identified ORFs conferring sensitivity to IQ (m > 1, q < 0.1, Supplementary Tables 4 and 7).While 54 were identified in duplicate screens of the original library, 502 were identified in duplicate screens of the "humanized" library.Genes conferring IQ sensitivity included those that function in nitrogen compounds and organic substance transport (19%); we speculate that transporting IQ into cellular compartments may enhance sensitivity.Genes that confer IQ sensitivity may also be required for or enhancing expression of CYP1A2 or NAT2, such as those encoding cytochrome B reductase (CBR1) and cytochrome C oxidase (COX14).Fifty-one ORFs appeared to either confer both IQ sensitivity (m > 1, q < 0.1) or resistance (m < 1, q < 0.1).Among these were the DNA repair genes, NTG1 and RAD18.Both trypan blue and competitive growth assays, however, indicate that these are validated resistance genes (Table 1).Both haploid rad18 and ntg1 strains also exhibit enhanced IQ sensitivity (Fig. 6), compared to the wild-type.Thus, NTG1 and RAD18 confer IQ resistance in both haploid and diploid strains.

Enhanced sensitivity of haploid double mutants defective in DNA repair
We hypothesize that 1 reason we did not identify more DNA repair genes is that there are multiple repair functions for IQ-associated DNA damage.To test this hypothesis, we made rad4 rad51, rad4 rad18, rad4 ntg1, and rad18 rad51 haploid double mutants.We compared the double mutants and the corresponding single mutant haploids using growth curves, competitive growth, and Genes were significant in at least 2 screens and were independently validated as conferring IQ resistance.Bolded genes were significant in at least 1 screen, and were independently validated as conferring IQ resistance.Genes were significant in at least 2 screens and were independently validated as conferring IQ resistance.Bolded genes were significant in at least 1 screen and were independently validated as conferring IQ resistance.viability assays.Competitive growth data for these double mutants, along with the haploid single mutants, expressing CYP1A2 and NAT2 are seen in Fig. 6.IQ-associated cell death data are shown in Supplementary Table 5.Growth curves for double and single mutants are shown in Supplementary Fig. 3. Compared to respective single mutants, the IQ sensitivity was most enhanced for the rad4 ntg1 and rad4 rad51 double mutants.Collectively these data suggest that NTG1 and RAD4, and RAD4 and RAD51 function in multiple DNA repair pathways for IQ-associated DNA damage.

Human orthologs of yeast IQ resistance genes associated with cancer
The genes appearing in at least 2 screens were checked for human orthologs using the Saccharomyces Genome Database (SGD, https:// www.yeastgenome.org).Human orthologs were then checked for disease association using Alliance of Genome Resources (https:// www.alliancegenome.org) to determine if any were implicated in human cancer (Sternberg et al. 2022).Of the human orthologs of IQ resistance genes, there were 11 genes associated with colon cancer and 10 genes associated with other or multiple cancers (Table 5).This is of particular interest because these compounds are often found in cooked red meats, which is a probable human carcinogen, and a risk factor for CRC.Biallelic mutations in the human ortholog of NTG1, NTHL1, are a high-risk factor for familial adenomatous polyposis (Groden et al. 1991).The human homolog for RAD18 has been noted to participate in multiple cancer progression stages and over-expression confers resistance to chemotherapeutic agents (Yang et al. 2018, Li et al. 2020).Alleles in other colon cancer risk-associated genes include SET2 and SHH3 (Habano et al. 2003, Choi et al. 2014).In addition, other genes were associated with multiple cancers, including breast and colon cancer.These include ADE6, which encodes phosphoribosylformylglycinamidine (FGAM) synthase, and TEP1, which is the yeast ortholog for the tumor suppressor gene PTEN.

Discussion
IQ is a class 2A and rodent carcinogen, present in charred meats in cigarette smoke.In this study, we profiled the yeast genome for IQ resistance by screening a "humanized" yeast diploid deletion library that expresses both CYP1A2 and NAT2.We identified 424 ORFs that appeared in at least 2 of 4 independent screens of the "humanized" library, and validated 82 using either growth curves, competitive growth assays, or viability assays.Major ontology groups included DNA binding, phosphatase, and ribosomal protein genes, while protein complexes included HIR and CK2.Human orthologs of 11 yeast genes are associated with colon cancer.We speculate that the colon cancer risk associated with these genetic factors could be increased by exposure to HAAs.
Several conclusions can be drawn from our studies screening the "humanized" and original yeast library.First, both screens identified ribosomal protein genes as resistance genes, thus suggesting that IQ targets ribosomes, and by extension, protein synthesis.Second, prominent among DNA repair genes were postreplication repair genes and NTG1, while NER genes were less prominent.Third, a comparison of the toxicity of IQ with other toxins underscores the role of several important genes and complexes involved in DNA damage and genotoxin tolerance, including RAD18 and the CK2 complex.
These conclusions were based on profiling the yeast genome for IQ resistance using a "humanized" library that expresses both CYP1A2 and NAT2.As expected, more ORFs were identified in the screen using a higher concentration of IQ in the "humanized" library and most genes identified at the lower concentration were also identified at the higher concentration.Nonetheless, there were several limitations to this study.First, many DNA repairdeficient strains were not detected even after exposure for 5 generations; a similar limitation was observed in a screen of a "humanized" library for AFB 1 resistance and was partially attributed to CYP1A2 expression, which confers an enhanced growth disadvantage (St John et al. 2020).For example, mre11, xrs2, and rad52 strains were not detected at all in the pooled samples, even without treatment.Second, genome profiling necessitated relatively high IQ concentrations that likely conferred toxicity by additional nongenotoxic mechanisms.While IQ sensitivity is enhanced in strains expressing CYP1A2 and NAT2, we do not know the exact proportion or amounts of activated IQ metabolites and all their associated targets.Among the ribosomal protein genes required for resistance to both "activated" IQ and IQ, the oxidative stress response gene, RPL22A (Chan et al. 2012), is present in both sets, while genes that function in rRNA export from the nucleus were only present among "activated" IQ resistance genes.Both IQ and activated IQ resistance genes encode rRNA proteins located in the large and small ribosomal subunits; however, there are more "activated" IQ resistance genes that function in a wider set of RNA translation functions and ribosomal protein transport (Supplementary Table 4 and Fig. 1).This set included genes that facilitate RNA translation, such as DEG1, modify tRNAs, such as PUS7, MOD5, and DUS1, and import ribosomal proteins into the nucleus, such as SYO1.Taken together, these data suggest that in comparison to resistance to IQ, resistance to activated IQ requires a larger set of genes involved in protein translation and rRNA transport.
Many ribosomal protein genes that confer IQ resistance also confer resistance to other toxins.Typical cellular reactions to multiple stresses involve downregulation of ribosomal protein gene transcription (Warner 1999;Gasch et al. 2000), which may render ribosomal protein mutants sensitive to mycotoxins and antifungal agents.Ribosome quality control is a central feature of resistance to multiple toxins, including trichothecene mycotoxins and cycloheximide (Kugler et al. 2016).Nonetheless, major cycloheximide resistance genes, RPL41 and RPL29 (Käufer et al. 1983;Dehoux et al. 1993;Bae et al. 2018), did not appear in our screen for IQ resistance but instead appeared to confer IQ sensitivity (Supplementary Tables 4 and 7).
The notion that activated IQ targets proteins is supported by observations that human serum albumin (SA)Cys 34 readily reacts with HAAs (Turesky and Le Marchand 2011;Bellamri et al. 2021) and that IQ reacts with heme-containing proteins (Turesky et al. 1987), suggesting that similar reactions could occur with yeast proteins.IQ-associated protein modification may explain why ammonia, amino acid metabolism, and protein synthesis were heavily represented among HAA resistance genes (Supplementary Table 4).Since protein modification genes (ALG3, ALG5, ALG6, ALG8, ALG12) were also identified in the IQ screen of the original deletion library (Supplementary Table 7 and Fig. 2), we presume that some cellular damage could also occur by IQ exposure per se or by IQ metabolites generated by endogenous yeast enzymes.
Our data supports the notion that replication bypass functions in IQ resistance, as was suggested by in vitro experiments using oligonucleotides containing IQ-N 2 Gua adducts (Stavros 2014).DNA repair genes that conferred IQ resistance were largely confined to those participating in postreplication repair.These included RAD18, REV7, and which are involved in DNA damage tolerance pathways that bypass DNA damage by errorprone polymerases (for review, see Xu et al. 2015).We also identified MPH1, a gene whose helicase function resolves stalled DNA structures, and which functions in error-free repair (Schürer et al. 2004).
RAD18 and REV7, which confer IQ resistance, are also common resistance genes for multiple cross-linking agents and agents that cause DNA adducts.These include trichloroethylene (De la Rosa et al. 2017) and AFB 1 (St. John et al. 2020), as well as chemotherapeutic drug derivatives, such as cisplatin, oxaliplatin, and mitomycin C (Wu et al. 2004).RAD18 may also function in error-free mechanisms that efficiently bypass N 2 -Gua adducts (Washington et al. 2004).The contribution of both error-free and error-prone mechanisms in bypassing IQ-associated DNA adducts will require additional studies.
Other protein kinases and phosphatases may have secondary roles in conferring tolerance to IQ-associated damage.Both the CK2 complex genes and PSY2 confer tolerance to multiple DNA-damaging agents (St. John et al. 2020).PSY2 encodes a phosphatase that deactivates Rad53 after exposure to a subset of DNA-damaging agents (O'Neill et al. 2007), while the CK2 complex is required for toleration of irreparable double-strand breaks (Toczyski et al. 1997).
In the repair of IQ-N 2 Gua or oxidized DNA bases, NER may not serve as a primary DNA repair mechanism, as suggested by Stavros et al. (2014).This may partially explain why rad4 mutants exhibited modest IQ sensitivity, while the double rad4 rad51 and rad4 ntg1 mutants exhibit enhanced IQ sensitivities.A synergistic relationship between NTG1 and RAD4 has been previously observed for oxidative stress resistance (Gellon et al. 2001).Further studies are necessary to determine whether the IQ-N 2 Gua adduct can be recognized by NER enzymes.
Eleven genes of the 424 ORFs that confer IQ resistance are associated with colon cancer (Table 4).Alleles mapping in human orthologs RAD18 and NTHL1, increase the risk of colon cancer.RAD18 Arg302Gln increases the colon cancer risk in Japanese patients (Kanzaki et al. 2007), while biallelic mutations in NTHL1, result in NTHL1-associated polyposis (Weren et al. 2015).Both NTG1 and NTHL1 encode bifunctional glycosylases, containing both glycosylase and lyase activity.NTG1 is noted to function in a redundant pathway to reduce G:C to T:A transversions, while in humans, formamido-pyrimidine residues, but not 8-oxyG, are recognized by NTHL1 (Das et al. 2020).Decreased activity is correlated with a mutation signature 30, where C is converted to T, while increased activity correlates with sequestration of the NER gene XPG, the mammalian RAD2 ortholog.These studies further underscore the primary role of base excision repair and a secondary role of NER in repairing oxidative DNA damage.
In summary, we profiled the yeast genome for IQ resistance using 2 different yeast libraries; 1 which was capable of bioactivation and the other which was not.Similarities between the IQ resistance genes indicated that protein was an IQ target, while in the "humanized" library, we identified DNA damage tolerance genes and NTG1.Considering that human orthologs NTHL1 and RAD18 are risk factors for colon cancer, our results provoke further studies suggesting that polymorphic NTHL1 alleles may confer HAA sensitivities.

Fig. 1 .
Fig. 1.Bioactivation of the HAA IQ by human CYP and NAT proteins.The chemical structures of IQ, N-hydroxy IQ, and N-acetoxy-IQ are shown in the figure.IQ is first hydroxylated by CYPs to form N-hydroxy-IQ.The hydroxylated IQ is then acetylated by N-acetyl transferases (NATs) to form N-acetyl-IQ.This compound is unstable and generates the nitrenium ion shown in the figure (for review see, Turesky and Le Marchand 2011).

Fig. 2 .
Fig.2.IQ sensitivity and expression of CYP1A2 and NAT2 in rad + and DNA repair deficient strains.Panel A is a Western blot that detects the human CYP1A2 protein present in yeast protein extracts.The black arrow indicates the position of CYP1A2, which is ∼58 kDa.The first and third lanes from the left contain the molecular weight standards.The second lane contains the protein extract from wild type (BY4743) and the last lane contains the extract from wild-type strain containing pCYP1A2_NAT2 (YB679).Panel B is a Western blot that detects the human NAT2 protein in yeast protein extracts.The red arrow indicates the position of the NAT2 protein, which is ∼37 kDa.The first lane from the left contains the molecular weight standards.The second contains a purified NAT2 protein.The third and fourth lanes, respectively, contain the wild-type (BY4743) and the wild type expressing pCYP1A2_NAT2 (YB679).Panel C shows the 7-MROD activity of CYP1A2 ±1 standard deviation, N = 3. Panel D shows the acetyltransferase activity of Nat2.Acetyl-CoA is a cofactor, without which the reaction does not occur.Error bars represent 1 standard deviation, N = 3. Panels E-H are growth curves of Rad + or rad4 rad51 yeast strains exposed to 2% MeOH (black), 400 µM IQ (Light red), or 800 µM IQ (Dark red), plotted as absorbance (A 600 ) against time (Hours).Error bars represent 1 standard deviation at 1-hour intervals.Panel E is a growth curve of the Rad + diploid strain containing pCYP1A2_NAT2 (YB679), N = 8.Percent growth with 400 µM IQ and 800 µM IQ were 94 and 83%, respectively.Panel F is a growth curve of the rad4 rad51 haploid yeast strain containing pCYP1A2_NAT2 (YB400), N = 3. Percent growth with 400 µM IQ and 800 µM IQ were 60 and 41%, respectively.Panel G is a growth curve of the Rad + diploid strain (BY4743), N = 8.Percent growth with 400 µM IQ and 800 µM IQ were 99 and 93%, respectively.Panel H is a growth curve the rad4 rad51 (YB226) strains, N = 3. Percent growth with 400 µM IQ and 800 µM IQ were 94 and 87%, respectively.

Fig. 3 .
Fig. 3. Venn diagrams of IQ resistance genes identified in screens that vary either exposure time or IQ concentrations.All ORFs were identified from the yeast diploid deletion library expressing CYP1A2 and NAT2.ORFs include genes of known function, genes of unknown function, and uncharacterized open Reading frames, which confer IQ resistance.Panel A is the overlap of the ORFs identified after 800 µM IQ exposure for 5 and 10 generations.Twenty-one ORFs were identified after 5 generations and 113 ORFs were identified after 10 generations.Panel B is the overlap of the IQ resistance genes identified after 400 µM and 800 µM IQ exposures for 10 generations.Individually, 671 ORFs were identified after 400 µM IQ exposure and 792 ORFs were identified after 800 µM IQ exposures.The original library treated with either IQ concentration identified 835 unique ORFs across both 400 µM and 800 µM IQ treatments.Panel C expands the combined original library from B into ORFs from the individual 400 µM (565 ORFs) and 800 µM IQ exposures (432 ORFs), displaying the degree of overlap.

Fig. 4 .
Fig.4.Protein interactome identified from screens of the diploid deletion library containing pCYP1A2_NAT2.The interactome was curated using String V11 (https://string-db.org,Szklarczyk et al. 2019).Nodes represent proteins, while edges represent interactions between 2 proteins.More evidence of interaction between proteins is denoted by thicker lines.Colored nodes indicate proteins involved in translation (green), DNA damage tolerance (blue), and amino acid metabolism (orange).Genes with human orthologs implicated in colon cancer are shown as diamond-shaped nodes.Disconnected nodes have been removed to improve readability.Panel A represents the functional interactome of proteins encoded by the 424 ORFs that were identified in both screens of the "humanized" library.Disconnected nodes in the bottom left were included because they have human orthologs implicated in colorectal cancer (CRC).Protein interactions with "high" confidence (0.7) of interaction are displayed.Panel B represents the functional interactome of each protein encoded by the 93 ORFs from the "humanized" yeast deletion library, which have been individually validated.Protein interactions with a "medium" confidence (0.4) of interaction are displayed.

c
Difference = percent of GFP-expressing cells when treated with IQ minus percent GFP-expressing cells when treated with MeOH; Bolded values indicate agreement that a gene confers resistance with P < 0.05.d Significance determined by comparing difference of tested strain with difference in the wild-type.e Difference = percent dead cells when treated with IQ minus percent dead cells when treated with MeOH; Bolded values indicate agreement that a gene confers resistance with P < 0.05.f Significance determined by comparing differential death of tested strains with the wild-type. g

Fig. 5 .
Fig.5.Growth curves and competition assays for selected IQ resistance genes identified in high-throughput screens.Panels A, B, C, and D show the wild-type strain, alt1, rpl33B, and sae2, respectively.The latter 3 represent major gene clusters highlighted in Fig.4, including nitrogen metabolism (ALT1), ribosome protein (RPL33B) and DNA repair (SAE2).The top row contains growth curves for selected strains containing pCYP1A2_NAT2.A 600 is plotted against time (Hours) with the background absorbance was subtracted, N ≥ 3. Error bars represent 1 standard deviation for 1-hour measurements.The middle row shows competitive growth data for the indicated strains containing pCYP1A2_NAT2.The y-axis shows the percent of GFP-containing wild type expressing CYP1A2 and NAT2 (YB676).Competitive cultures were inoculated into SC-Ura at approximately 10% YB676 and 90% the indicated strain. 2 × 10 4 cells were acquired per sample.Error bars represent 1 standard deviation.For wild type (YB679), N ≥ 5; for all others, N = 2. Significance was determined by Dunnett's test and asterisks indicate P < 0.05 (*), P < 0.01 (**) and P < 0.001(***).The bottom row shows competitive growth data for the indicated strains not expressing CYP1A2 or NAT2.The y-axis shows the % of GFP-containing WT (YB675).Competitive cultures were inoculated into SC at approximately 10% YB675 and 90% the indicated strain and grown for 20 hours. 2 10 4 cells were acquired per sample.For wild type (BY4743), N = 12; for all others, N = 2. Error bars represent 1 standard deviation.Significance was determined by Dunnett's test and asterisks indicate P < 0.05 (*), P < 0.01 (**) and P < 0.001(***).
a Bolded genes were showed statistically significant differences in AUCs in growth curves.b

Table 2 .
Component GO terms for validated IQ resistance genes.Genes from the 82 validated ORFs out of the total number of yeast ORFs belonging to each GO group, based on reference list of 7166 ORFs.
b c

Table 3 .
Function and process GO terms for validated IQ resistance genes.Genes from the 82 validated ORFs out of the total number of yeast ORFs belonging to each GO group, based on reference list of 7166 ORFs.
b c

Table 4 .
Validated checkpoint, DNA damage repair, and DNA damage-inducible genes.Asterisk indicates genes not directly involved in DNA repair, but whose expression is induced in response to DNA damage. a

Table 5 .
Yeast genes conferring resistance to activated IQ with human orthologs associated with cancer.
GPT, Glutamic-Pyruvic Transaminase.a Asterisk (*) indicates genes that have not been independently validated.b Bolded genes are associated with colorectal cancer.