Selection of chromosomal DNA libraries using a multiplex CRISPR system in Saccharomyces cerevisiae

The directed evolution of biomolecules to improve or change their activity is central to many engineering and synthetic biology efforts. However, selecting improved variants from gene libraries in living cells requires plasmid expression systems that suffer from variable copy number effects, or the use of complex marker-­‐dependent chromosomal integration strategies. We developed quantitative gene assembly and DNA library insertion into the Saccharomyces cerevisiae genome by optimizing an efficient single-­‐step and marker-­‐free genome editing system using CRISPR-­‐Cas9. With this Multiplex CRISPR (CRISPRm) system, we selected an improved cellobiose utilization pathway in diploid yeast in a single round of mutagenesis and selection, which increased cellobiose fermentation rates by over ten-­‐fold. Mutations recovered in the best cellodextrin transporters reveal synergy between substrate binding and transporter dynamics, and demonstrate the power of CRISPRm to accelerate selection experiments and discoveries of the molecular determinants that enhance biomolecule function.


Introduction
Directed evolution using living systems allows selections for improved biomolecule function to be directly coupled to phenotype. However, present directed evolution systems require the use of extrachromosomal gene libraries encoded in bacteriophage or plasmids, which suffer from high levels of copy number variation (Yokobayashi et al. 2002;Turner 2009;Esvelt et al. 2011). For example, in the yeast Saccharomyces cerevisiae, common yeast centromere containing (CEN) plasmids vary widely in copy number per cell resulting in highly variable expression levels from cell to cell (Figure 1). Ideally, the selective pressure applied to evolve improved or new biomolecule function should be limited to the DNA sequence level and not to gene copy number in a given cell (Zhou et al. 2012). To overcome copy number variation in yeast, the standard genetic tool of integrating linear DNA into the genome by homologous recombination is too inefficient to insert DNA libraries without multiple steps that rely on dominant selectable markers for chromosomal integrations (Wingler and Cornish 2011) or specialized strain backgrounds .
Bacterial type II CRISPR--Cas9 genome editing has been used successfully in several eukaryotic organisms (Cong et al. 2013;DiCarlo et al. 2013;Mali et al. 2013;Mali et al. 2013;Shalem et al. 2014) but has not been adapted for selection experiments. Genome editing CRISPR--Cas systems require a Cas9 endonuclease that is targeted to specific DNA sequences by a non--coding single guide RNA (sgRNA) (Jinek, et al. 2012;Jinek et al. 2014). The Cas9--sgRNA ribonucleoprotein complex precisely generates double--strand breaks (DSBs) in eukaryotic genomes at sites specified by a twenty--nucleotide guide sequence at the 5' end of the sgRNA that base pairs with the protospacer DNA sequence preceding a genomic Protospacer Adjacent Motif (PAM) . The presence of the Cas9--produced DSB in genomic DNA can increase the rate of homology--directed repair (HDR) with linear DNA at the DSB locus by several thousand--fold (Storici et al. 2003;) potentially enabling high--throughput experimental methods. CRISPR--Cas9 could therefore be useful in probing industrial Saccharomyces cerevisiae strains, which are more robust compared to laboratory strains (Kerr and Service 2005;Farrell et al. 2006;Rubin 2008), but for which few genetic tools are available, due to the fact that these strains are often polyploid with low--efficiency mating and sporulation.

Results
To optimize the efficiency of HDR in S. cerevisiae mediated by Cas9, we first constructed a plasmid expressing Streptococcus pyogenes Cas9 fused to a nuclear localization signal (Cas9--NLS) that efficiently localizes to the nucleus of yeast cells (Figure 2A, B) (Jinek et al. 2012). We used an intermediate strength promoter, PRNR2, to express Cas9 for genome targeting experiments, because Cas9 (Cas9--NLS--His8x) expression from the PRNR2 promoter resulted in yeast strains with near wild type fitness whereas Cas9 expressed using strong yeast promoters such as PTDH3 (Lee et al. 2013) reduced yeast fitness relative to wild type cells ( Figure 2C).
Second, we designed a new sgRNA architecture in order to improve its expression All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; and function. Cellular levels of sgRNAs correlate with the efficiency of Cas9-mediated genome targeting in mammalian cells . Controlled expression of nuclear--localized RNAs and the correlation between sgRNA levels and Cas9--editing in yeast has not been explored. To better control cellular levels of correctly folded sgRNA, we developed a modular design by fusing the sgRNA  to the 3' end of the self--cleaving hepatitis delta virus (HDV) ribozyme (Ke et al. 2007;Webb et al. 2009) (Figure 2D). We reasoned that the structured ribozyme would protect the 5' end of the sgRNA from 5' exonucleases (Houseley and Tollervey 2009). Additionally, the HDV ribozyme would cleave the RNA immediately 5' of the ribozyme, removing extraneous RNA sequences and allowing flexibility in which promoters can be used, such as tRNAs whose DNA sequences also serve as RNA polymerase III promoters (Marck et al. 2006). Ribozyme removal of the tRNA would then separate transcription initiation from within the tRNA sequence (Marck et al. 2006) and tRNA nuclear export (Köhler and Hurt 2007) from the process of forming functional sgRNA complexes with Cas9.
To quantify the efficiency and specificity of HDR with the new Cas9/sgRNA format, phenotype screens were performed by the co--transformation of a single plasmid that expresses the Cas9 protein and one or more HDV ribozyme--sgRNAs, containing a 20--nucleotide guide sequence that matches a specific site in the yeast genome, along with a linear double--stranded DNA ( Figure 2E). The linear dsDNA, which contains a unique 20--mer barcode (Giaever et al. 2002) flanked by common primer sequences and fifty base pairs of DNA homologous to the regions proximal to All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; the PAM motif ( Figure 2F), acts as a template for DNA repair by HDR, resulting in a unique markerless and barcoded insertion allele. With a single sgRNA, we found tRNA sequences used as RNA polymerase III promoters resulted in nearly 100% efficient barcode insertion in diploid yeast at the URA3 locus, resulting in 5-fluororotic acid resistance (Boeke et al. 1984) whereas the non--tRNA promoters were mostly ineffective ( Figure 3A and Supplementary file 1A). We also assessed target sequence bias by targeting 11 unlinked yeast genes in diploid S288C yeast cells using tRNA Tyr as the RNA polymerase III promoter (Supplementary file 2A).
We found that targeting was highly efficient for every genomic target except LEU2, which at first was weak but was restored to 100% by selecting a different guide sequence (LEU2-2) as the sgRNA ( Figure 3B). We also verified that our CRISPR--Cas9/sgRNA system did not result in non--specific genome targeting by sequencing 9 of the targeted strains (Supplementary file 1B). We also tested efficacy of our approach in the polyploid industrial Saccharomyces cerevisiae strain ATCC4124 that has superior tolerance and productivity phenotypes (Ness et al. 1993). The non--tRNA promoter previously used in haploid yeast , PSNR52, functioned in diploid S288C yeast but failed to result in targeting in the polyploid industrial yeast ATCC4124; in contrast, the tRNA Pro promoter enabled highly-efficient barcode insertion into all copies of the URA3 locus in ATCC4124, which we term cis--multiplexing, i.e. introducing a double--stranded DNA break at a single genomic locus across all chromosomes ( Figure 3C and Supplementary file 1A).
Notably, tRNA Tyr was inefficient at targeting in ATCC4124 though it was effective in S288C diploid (Supplementary file 1A), indicating that the choice of tRNA used for All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; the expression of the HDV--sgRNA impacts multiplexing in a strain--dependent manner. Since the efficiency of cis--multiplexing with our CRISPR--Cas9/sgRNA system is nearly 100% using one sgRNA, we next tested its ability for trans--multiplexed CRISPR--Cas9 genome editing, i.e. introducing double--stranded DNA breaks at multiple genomic loci simultaneously, using tRNAs to drive the expression of the HDV--sgRNAs. We simultaneously targeted two and three unlinked genomic loci in diploid cells for loss--of--function mutagenesis by introducing a single plasmid expressing Cas9 and two or three sgRNAs with two or three gene--specific barcoded DNA molecules, respectively. The efficiency of duplex and triplex targeting (requiring 4 and 6 chromosome cuts, respectively) with two sgRNAs (URA3 and LYP1) or three (URA3, LYP1 and COX10) was 43% and 19%, respectively ( Figure 3D and Supplementary file 1A). This multiplexed efficiency was dependent on the presence of the HDV ribozyme 5' of the guide sequence as duplex targeting dropped to 3.5% in cells with sgRNA lacking the 5' HDV ribozyme ( Figure 3D and Supplementary file 1A). We measured the relative cellular abundance of sgRNAs expressed by PSNR52, with and without the 5' HDV ribozyme, by reverse transcription quantitative polymerase chain reaction (RT--qPCR) and found that the presence of the ribozyme increases the intracellular abundance of the sgRNAs by six--fold ( Figure 3E), consistent with the model that sgRNA abundance is rate limiting for CRISPR/Cas targeting . Haploid trans--multiplexing showed a minimal loss in activity scaling from one to three loci (up to three chromosome cuts) All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; ( Figure 3F and Supplementary file 1A). These experiments demonstrate that our multiplexed CRISPR/Cas9 system, which we term CRISPRm, is powerful enough to generate multiple marker--free targeted mutations in the yeast genome in a single experiment and that the 5' tRNA--HDV ribozyme sequence is required for higher order multiplex targeting.
Having established CRISPRm for cis--and trans--multiplex targeting, we next tested its capabilities in engineering genes and pathways integrated in the yeast genome. We first tested whether CRISPRm could be used to assemble foreign genes as chromosomal integrations. We tested the efficiency of in vivo assembly of a functional nourseothricin--resistance (Nat R ) gene from three overlapping PCR products encoding a transcription promoter, protein open reading frame (ORF) and transcription terminator at a single locus using cis--multiplexing in the diploid S288C strain. The efficiency of Cas9--mediated integration and assembly of these three DNA fragments to the correct (URA3) locus was measured by a combination of 5-fluoroorotic acid resistance (5--FOA R ) and Nat R . We found 85% efficiency of targeting and assembly in both copies of the URA3 locus in diploid yeast S288C cells, and 70% in polyploid ATCC4124 by using the tRNA Pro sequence as the sgRNA promoter ( Figure 3C and Supplementary file 1A), whereas PSNR52 --dependent targeting was weaker in S288C cells and non--existent in ATCC4124 cells ( Figure 3C and Supplementary file 1A). Thus, CRISPRm enables the one--step markerless assembly of functional genes from multiple fragments into the S. cerevisiae genome, resulting in homozygous insertions into both diploid and polyploid yeast strains. All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; Finally, we tested whether CRISPRm could be used for in vivo selections of improved protein function. We targeted a cellodextrin utilization pathway from the cellulolytic fungus Neurospora crassa, comprised of the cellodextrin transporter (cdt-1) and an intracellular β--glucosidase (gh1-1), to assemble in vivo in diploid S. cerevisiae. In combination, these two genes result in yeast cells capable of consuming the disaccharide cellobiose, which could be used to improve biofuel production (Galazka, et al. 2010). To select for improved cellobiose utilizing strains we used error--prone PCR to amplify the cdt-1 gene. We co--transformed the library of mutated cdt-1 alleles, along with overlapping dsDNA fragments for the promoter and terminator, for assembly and integration into the URA3 locus (ura3::PPGK1-cdt-1-TADH1 library) of a diploid strain with a previously integrated β--glucosidase gh1-1 gene (lyp1::PTDH3-gh1-1-TCYC1) using CRISPRm. We then grew the transformants in medium containing cellobiose as the sole carbon source to select for functional cdt-1 alleles and identified one strain with 2.6--fold (homozygous cdt-1 N209S/I354N/S360P/T406S/W531L allele) and one strain with 1.7--fold (homozygous cdt-1 Q45H/F262Y/F533Y allele) increased cellobiose utilization capacity over wild type cdt-1 ( Figure 4A). The mutations were mapped in CDT--1 onto structures of the homologous major facilitator superfamily (MFS) hexose--like transporters XylE from Escherichia coli (Sun et al. 2012) and GlcP from Staphylococcus epidermidis (Iancu et al. 2013) (Figure 4B). Notably, one site maps to a site that likely interacts directly with the sugar substrate (N209S) (Figure 4B).
Whereas two mutations in the 2 nd allele (Q45H, F533Y) map to non--conserved All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; regions of the transporters, one (F262Y) maps to a site immediately adjacent to the nearly universally--conserved PESPR motif in hexose transporters that is involved in transporter dynamics (Sun et al. 2012). Furthermore, the mutation from phenylalanine to tyrosine results in an expanded motif (PESPRY) that is present in all of the major hexose transporters in Saccharomyces cerevisiae ( Figure 4B). We then used CRISPRm to introduce the N209S mutation (nucleotide G626A) in the cdt-1 Q45H/F262Y/F533Y diploid yeast background. The quadruple mutant cdt-1 allele resulted in synergistic cellobiose utilization with 3.9--fold increased growth over strains expressing wild--type cdt-1 ( Figure 4A). We next tested the selected cdt-1 alleles in anaerobic fermentations using cellobiose as the sole carbon source. We constructed a cdt-1 N209S/F262Y double--mutant in diploid yeast using CRISPRm, and compared it to wild--type cdt-1 and the cdt-1 Q45H/N209S/F262Y/F533Y quadruple mutant. We found that wild--type cdt-1 integrated into the chromosome was barely able to support cellobiose fermentation ( Figure 4C). By contrast, the cdt-1 Q45H/N209S/F262Y/F533Y quadruple mutant or the cdt-1 N209S/F262Y double--mutant expressed from the chromosome resulted in complete fermentation of cellobiose, comparable to wild-type cdt-1 expressed from a high--copy 2µ plasmid ( Figure 4C) (Ha et al. 2011).

Discussion
We have demonstrated that CRISPRm can be used to quickly generate markerless loss--of--function alleles, heterologous insertion alleles, allele swaps and engineer proteins in yeast by in vivo selections for greatly improved metabolic activity. Although we selected for improved protein function, as suggested by the All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; location of the mutations in CDT--1 (Figure 4B), it is also possible that selection enriched for sequences that confer improved mRNA stability or translation efficiency, due to the fact that the selection was carried out in vivo. Future experiments will be needed to distinguish between these possibilities.
CRISPRm could be used to integrate DNA libraries to interrogate cis-regulatory elements such as transcription promoters, 5'--and 3'--untranslated regions (5'--UTRs and 3'--UTRs) in mRNAs, or riboswitches. CRISPRm enables these types of experiments because it allows for near quantitative integration of DNA libraries, which can be a single fragment or complex assemblies of multiple fragments. This could enable integration and selection of multi--gene pathways, or specific mutation and selection of subsets of a protein or multiple DNA--encoded functions.
The ability to use chromosomally--integrated DNA libraries removes a key limitation in in vivo selection experiments by eliminating biomolecule expression variation due to variable plasmid copy numbers. Although we used a relatively small library size in the present experiments, even in this small library we observed impressive gains in transporter function. We observed 2--3 mutations per kilobase in the selected transporter genes, resulting in 3--5 amino acid substitutions per transporter. We were then able to combine the mutations in a single transporter with superior fermentation performance. Library sizes introduced by CRISPRm should be scalable to the limits of transformation efficiency, to 10 7 or larger (Kawai All rights reserved. No reuse allowed without permission. was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; et al. 2010). CRISPRm should also enable multiple rounds of selection for directed evolution experiments. We envision that each cycle would involve error--prone PCR of the gene of interest, either from a starting sequence or an enriched and improved cell population, followed by CRISPRm--mediated integration of the amplified gene and selection. This cycle would compare favorably with that required for plasmid-based selections, yet would avoid the problem of plasmid copy number variation.
CRISPRm should also be adaptable to screens that rely on sorting large populations of cells, i.e. by fluorescence--activated cell sorting. We envision that CRISPRm could improve the discovery of human therapeutics by yeast display (Krugel et al. 1988), and could be used to map the importance of protein--protein interactions (Ovchinnikov et al. 2014) via multiplexed library insertion.
To achieve high efficiency multiplexing with CRISPRm, we found it was necessary to modify the sgRNA design by including a 5' HDV ribozyme, and to assay multiple RNA polymerase III promoters for a given strain background. The ease of discovering tRNAs in any genome and the universality of the sgRNA construct means that, with few modifications, CRISPRm can be used directly in non--model fungal hosts, such as pathogens or organisms used in the biotechnology industry.
The genetics of many of these organisms have not been studied in any depth due to the technological limitations of available genetic manipulation techniques. For example, industrial Saccharomyces cerevisiae strains are more stress tolerant and produce much higher yields of desired biofuel or renewable chemical end products than laboratory strains (Kerr and Service 2005;Farrell et al. 2006;Rubin 2008). All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; However, linking genotypes of industrial yeasts to their phenotypes remains difficult. CRISPRm should serve as a rapid and high throughput means for connecting the genotypes of these organisms to their phenotypes, i.e. for generating marker--free barcoded alleles for large--scale pooled fitness studies of loss--of-function mutants in these organisms (Giaever et al. 2002). In the future, CRISPRm may also be applicable to interrogating mammalian cells, if the levels of HDR are sufficiently high. Finally, the ability to use CRISPRm for multiplexed targeting paves the way for applying directed evolution to cellular pathways and genetic circuits for higher order synthetic biology applications in any host strain.
All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.

Cloning the pCAS plasmid backbone
Gibson Assembly Mastermix (E2611L)(New England Biolabs, Ipswich, MA) (Gibson et al. 2009) was used to fuse the KANMX (Yeastdeletionpages.com) cassette to the pUC bacterial origin of replication from pESC--URA (Agilent Technologies, Santa Clara, CA). Restriction -free (RF)(van den Ent and Lowe 2006) cloning was used to add a yeast 2µ origin of replication from pESC--URA to the pCAS backbone.
The resulting pCAS backbone plasmid was propagated in yeast to confirm functionality.

Cas9 expression constructs
The Cas9 gene from Streptococcus pyogenes was amplified from clone MJ824 (Jinek et al. 2012) and cloned into the pCAS backbone plasmid by RF cloning. A yeast nuclear localization signal (NLS) sequence, codon optimized using IDT software (Integrated DNA Technologies, Coralville, IO), was then cloned into the plasmid by RF cloning. Additional elements fused by RF cloning to the Cas9--NLS sequence included the GFP gene, the CYC1 terminator from S. cerevisiae strain S288c (Yeastgenome.org) and the promoters (800 base pairs upstream) from the genes TDH3, TEF1, RNR2 and REV1, also taken from strain S288c (Lee et al. 2013). For genome editing experiments, the GFP sequence was removed from the Cas9 gene and replaced with a C--terminal His8 affinity tag, by RF cloning.

Engineering of sgRNA constructs
All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; Synthetic DNA (Integrated DNA Technologies, Coralville, IO) for the sgRNA and for a catalytically active form of the Hepatitis Delta Virus (HDV) Ribozyme was sequentially cloned by RF cloning into the Cas9 containing vector (pCAS). The terminator (200bp) of SNR52 (Yeastgenome.org) was cloned 3' of the ribozyme--sgRNA sequence by RF cloning. A series of RNA polymerase III (Pol III) promoters were PCR amplified from S288c genomic DNA and cloned 5' of the ribozyme--sgRNA sequence by RF cloning. The tRNA promoters included the full--length tRNA plus 100 base pairs upstream of the tRNA gene (Supplementary file 1C). The sgRNAs used for multiplex targeting were PCR amplified using primers containing 5' and 3' restriction sequences and sub--cloned into pCAS by ligation dependent cloning into SalI, SpeI and SacII unique restriction sites.
Yeast strains used in this study S. cerevisiae strain S288c (ATCC 204508) was used as the haploid, and then mated to form the homozygous diploid for the diploid experiments. Yeast strain ATCC4124 is an industrial polyploid yeast isolated from a molasses distillery, and was obtained from ATCC.

Cas9-GFP Localization and Expression
Expression and localization of Cas9--GFP was verified by imaging haploid S.
cerevisiae S288c cells transformed with pCas9--GFP::KAN using fluorescence microscopy (Leica Epifluorescence, Leica Microsystems, Buffalo Grove, IL). Cells were grown overnight and nuclear localization visualized at 100X magnification. All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; were re--suspended in 250 µL YPD at 30 °C for two hours and then the entire All rights reserved. No reuse allowed without permission.

Fitness analysis of Cas9 expressed by different promoters
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
Cells were grown for 48 hours at 37 °C, and colonies imaged using the Biorad ChemiDoc Imager (Biorad, Hercules, CA) before replica plating onto phenotype-- CAN1 mutants were selected on 2.0 g/L Yeast nitrogen base without amino acids or ammonium sulfate, 5.0 g/L Ammonium sulfate, 1.0 g/L CSM--arginine (MP Biosciences 4510--112), 20 g/L Glucose, 20 g/L Agar + canavanine sulfate (50mg/L) (Sigma C9758); the remaining auxotrophic mutants were selected on 2.0g/L Yeast nitrogen base without amino acids or ammonium sulfate, 5.0 g/L Ammonium sulfate, 1.0 g/L CSM , 20 g/L Glucose, 20 g/L Agar; and aerobic respiration deficient mutants (petites) were selected on 20 g/L Peptone, 10 g/L Yeast Extract, 20 g/L Agar, 0.15 g/L Adenine hemisulfate, 20 g/L Glycerol (Sigma G5516). Colonies from the YPD+G418 plates were picked and grown overnight in 0.8mL of YPD. Genomic DNA was extracted from these cultures using the MasterPure Yeast DNA Extraction Kit (Epicentre MPY80200). PCR confirmation of the 60--mer integration allele was All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; performed using primers flanking the target site. PCR products were purified by Exo--SAP--IT (Affymetrix 78201) and Sanger sequenced at the UC Berkeley, Sequencing Facility (Berkeley, CA) to confirm barcode sequence in the amplicon.

RT-qPCR of sgRNAs
Cells containing the pCAS plasmid with sgRNA inserts were grown in 900 µL of YPD+G418 medium for 24 hours at 30 C and 750 rpm. Total RNA was extracted from exponentially growing yeast cells using Ambion RNA RiboPure TM Yeast Kit

Multiplex genome targeting by Cas9
All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; Multiplex targeting was performed as described above using pCAS plasmids containing more than one sgRNA expression construct cloned into one of the restriction sites by ligation dependent cloning. Single versus double mutant efficiency was scored relative to the number of colonies present on the YPD+G418 plate. Genomic DNA isolation and PCR of the integration site was performed as described above.

Multiplex in vivo assembly of DNA using Cas9
Drug resistance cassettes were assembled in vivo from three linear double-stranded DNA fragments PCR amplified from the Ashbya gossypii TEF1 promoter (AgPTEF1), the nourseothricin open reading frame (Nat R ) and Ashbya gossypii TEF1 (AgTTEF1) terminator in separate reactions. The primers used to amplify the promoter and terminators contained 50 bp of homology to the nourseothricin ORF and 50 bp of homology to the genomic target. Colonies exhibiting drug resistance (nourseothricin 100 mg/L) (Goldbio N--500--1) following replica plating were compared to the number of colonies on the YPD+G418 to determine efficiency of multiplex assembly.

Selection of improved variants of cellodextrin transporter CDT1
To generate CDT1 mutant allele libraries, the GeneMorph II Random was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; ScPPGK1 promoter and ScTADH1 terminator into a diploid S288c yeast strain containing a previously--integrated gh1-1 gene. The gh1-1 gene included the ScPTDH3 promoter, the N. crassa gh1-1 open reading frame and ScTCYC1 terminator. The primers used to amplify the promoters and terminators contained 50 bp of homology to either the cdt-1 or gh1-1 ORFs and 50 bp of homology to the respective the genomic targets. Five µg of each of the three PCR products (promoter, open reading frame, terminator) were co--transformed with the pCAS plasmid containing the URA3 guide sequences and screened for G418 resistance as described. Approximately 1600 G418--resistant colonies were pooled and resuspended in minimal cellobiose medium (SCel) (2.0g/L Yeast nitrogen base without amino acids or ammonium sulfate, 5.0 g/L Ammonium sulfate, 1.0 g/L CSM , 20 g/L Cellobiose). Resuspended cells were immediately spread evenly on SCel plates for initial analysis prior to cellobiose selection (t=0 samples). Ten microliters of the resuspended cells were inoculated in 50 mL of SC medium in biological triplicate. Cells were harvested after three days and spread onto SCel plates. Cells were grown at 30 C for four days, and the largest colonies were chosen for further analysis.

Tecan growth analyzer and fitness calculation of cdt-1 alleles
Cells were grown overnight in 0.5 mL of Synthetic Dextrose (2%) (SD) in 96 well plates. Cultures were diluted 1:500 in SCel media (4% cellobiose) and 150 µL were grown using the Tecan Sunrise (Tecan Systems Inc., San Jose, CA) in biological triplicate for three days at 30 C. Average and standard deviation was calculated for each biological sample. Relative fitness was calculated by measuring area between All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.

Generating the double mutant and quadruple cdt-1 alleles
Integration of the repair oligonucleotide into the chromosomal wild--type cdt-1 gene or cdt-1 Q45H/F262Y/F533Y was performed as described for integrating barcoded DNA. The sgRNA sequence used to target the wild type allele of cdt-1 was cloned into the pCAS vector with the sequence 5'--TGCACTGGCTTCTACAACTG--3'. The repair oligonucleotide used to make the A626G (i.e. N209S) mutation was 5'--CGGCCGCTGCACTGGCTTCTACAGCTGCGGTTGGTTCGGAGGTTCATTCCTGCCGCCTG --3' with fifty more base pairs of homology to cdt-1, on both sides of the above oligonucleotide. Genomic DNA was extracted as described and Sanger sequencing confirmed the incorporation of the A626G mutation. Similarly, the T785A mutation (F262Y) was incorporated into the cdt-1 N209S allele using an sgRNA with guide sequence 5'--CCTCGCTTCCTATTTGCCAA--3', and a repair oligonucleotide with sequence 5'--AGAATCCCCTCGCTACCTATTTGCCAACGGCCGCGACGCTGAGGCTGTTGCCTTTCTTGT--3'. Genomic DNA was extracted as described and Sanger sequencing confirmed the incorporation of the T785A mutation.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
Synthetic defined media lacking leucine and uracil (SD--LU) was made by adding 6.7 g/L Difco Yeast Nitrogen Base without amino acids; 2 g/L Drop--out Mix Synthetic Minus Leucine, Uracil without Yeast Nitrogen Base (US Biological); and 20 g/L Dextrose to distilled water.
Colonies were picked into SD--LU media and grown to exponential phase at 30 °C. Cultures were concentrated by centrifugation, spotted onto plain glass slides, All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; and examined on a Zeiss Observer D1 microscope using a 100x DIC objective. Images were captured using a Hamamatsu Orca--flash 4.0 (C11440) camera using auto--exposure. Fluorescence images were taken using an X--Cite Series 120 lamp, Zeiss filter set 45 (excitation at 560/40 nm and emission at 630/75 nm) for mRuby2, and Zeiss filter set 46 (excitation at 500/20 nm and emission at 535/30 nm) for Venus. Images were analyzed and composites were created using Fiji (http://fiji.sc).

Illumina sequencing to identify potential off-target mutations
To rule out that our CRISPR--Cas9/sgRNA system resulted in non--specific genome targeting, we performed whole genome sequencing for URA3 and LYP1 was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. corresponding to the 20 nucleotides immediately 5' of the "NGG" PAM site), and obtain the genome coordinates ranging from 30 nucleotides upstream and downstream of the PAM site. Cas9 target sequences were added to VCF files as custom annotations using snpEff (v3.3h), and SnpSift (v3.3h) was used to extract desired fields into tables for analysis with a custom R script (Supplementary file   4). Needleman-Wunsch global alignments between our guide sequences and Cas9 target sequences were performed using the pairwiseAlignment function (Biostrings package, Bioconductor) in R, with a substitution matrix of --1 for mismatches and 2 All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; for matches, produced with the nucleotideSubstitutionMatrix function (Biostrings package, Bioconductor). The probability of there being a better match for the guide sequence to a given Cas9 target sequence was calculated as the frequency of Cas9 target sequences with better alignments to the same guide sequence, amongst 10,000 randomly selected Cas9 target sequences. To compile counts of all variants and various subclasses, a GATKReport was generated from the VCF files with GATK's VariantEval walker, read into R using the "gsalib" library, and the desired categories were extracted with a custom R script (Supplementary file 4).
Without the repair DNA template required for HDR, the majority of mutations caused by Cas9 are expected to be INDELs, SNPs or MNPs that initiate within the protospacer sequence, which is the 20 nucleotide (nt) sequence 5' of the PAM (Fu et al. 2013). Therefore, for completeness, we searched for all PAM sites within 30 nt upstream and downstream for each of the 21 variants. We then compared the 20 nt URA3 and LYP1 guide sequences to the putative protospacer sequences within the regions flanking each variant. An end--to--end alignment identified 10 or fewer nucleotide matches between the URA3 or LYP1 guide sequences and the variant sequences (Supplementary file 1B). Cas9 requires at least 12 perfect base pair matches within the guide--target sequence (Hsu, Scott et al. 2013) so it is highly unlikely that the URA3 and LYP1 guide sequences directed Cas9 to any of these potential protospacers that lie within variant sites. These mutations likely arose as natural variants during the course of the experiment. All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; As a second method to evaluate the likelihood of off--target mutations, we performed local alignments of our guide sequences to all Cas9 target sequences whose PAM site was within 30 nucleotides upstream and downstream of a detected variant, as well as to 10,000 randomly selected Cas9 target sequences from the genome. Since our guide sequences are expected to have a better match to 13% or more of all Cas9 target sequences (~126,000 or more sites) than to the best matching Cas9 target sequence with a nearby variant, and the number of nucleotide matches in end--to--end alignments is at most 10, we argue that the variants was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; mL serum flasks. After inoculation, the flasks were sealed with rubber stoppers and clamped with an aluminum seal. To achieve anaerobic conditions, the headspaces of the sealed flasks were purged with nitrogen gas for 30 minutes. These were then cultivated at 30 °C and 220 rpm. was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; and (D) integrated in tandem at a single locus (URA3) in the genome. Each image was auto--exposed for both red and green channels, with yellow showing co--All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; expression of mRuby2 and Venus. The variability in the relative expression of the two fluorescent proteins is reduced by moving from two plasmids to one, and the variability in total expression is reduced by moving to the genome. The difference between integrating into one or two loci in the genome is minimal.
All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
(E) Schematic of yeast Cas9 targeting. Cas9 and sgRNA are expressed from a single plasmid, form a complex, and cleave targeted genomic DNA, which is repaired using a barcoded oligonucleotide. (F) The linear barcoded repair DNA molecule. Each repair DNA contains a 5' TAA stop codon (gold), a forward primer sequence (green), a unique molecular barcode (grey), and a 3' reverse primer (green) (Giaever et al. 2002). Barcodes are amplified using a forward primer that contains 50bp of homology (blue) to the genome targeting site and a reverse primer that contains 50bp of homology to the genome targeting site. The 50bp of genomic targeting sequence are each 10bp proximal to the PAM motif, resulting in a 20nt deletion and barcode oligonucleotide integration.
All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; Single fragment barcode integration (blue) and three--fragment Nat R cassette All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
For each experiment the promoter and strain are indicated as promoter--strain (e.g. All rights reserved. No reuse allowed without permission.

SNR52--S288C). (D) Efficiency of multiplex insertion of barcoded
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014; was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which . http://dx.doi.org/10.1101/007310 doi: bioRxiv preprint first posted online Jul. 21, 2014 in the inward--facing configuration (lower right, PDB entry 4LDS). Amino acids in parentheses are the sequences in XylE or GlcP. (C) Fermentation of cellobiose by wild--type cdt-1 expressed from a chromosomally--integrated copy in diploid S288C yeast. Cellobiose is indicated as G2, ethanol as EtOH. The rate of cellobiose consumption was 0.13 g L --1 hr --1 . (D) Fermentation of cellobiose by mutant versions of cdt-1 expressed from chromosomally--integrated copies. Glucose consumption (grey) was identical in both strains. The rate of cellobiose consumption for both strains was over 2.0 g L --1 hr --1 . In panels C and D, values and error bars represent the means and standard deviations of three independent biological replicates.
All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
The copyright holder for this preprint (which Supplementary file 4. Custom R script for yeast genome analysis.
All rights reserved. No reuse allowed without permission.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.
was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.