Weimin Zhang1, Ran Brosh1, Laura H McCulloch1, Yinan Zhu1, Hannah Ashe1, Gwen Ellis1, Brendan R Camellato1, Sang Yong Kim2, Matthew T Maurano1,2, Jef D Boeke1,3,4. 1. Institute for Systems Genetics, NYU Langone Health, New York, NY 10016, USA. 2. Department of Pathology, NYU Langone Health, New York, NY 10016, USA. 3. Department of Biochemistry and Molecular Pharmacology, NYU Langone Health, New York, NY 10016, USA. 4. Department of Biomedical Engineering, NYU Tandon School of Engineering, Brooklyn, NY 11201, USA.
Abstract
Overwriting counterselectable markers is an efficient strategy for removing wild-type DNA or replacing it with payload DNA of interest. Currently, one bottleneck of efficient genome engineering in mammals is the shortage of counterselectable (negative selection) markers that work robustly without affecting organismal developmental potential. Here, we report a conditional Piga knockout strategy that enables efficient proaerolysin-based counterselection in mouse embryonic stem cells. The conditional Piga knockout cells show similar proaerolysin resistance as full (non-conditional) Piga deletion cells, which enables the use of a PIGA transgene as a counterselectable marker for genome engineering purposes. Native Piga function is readily restored in conditional Piga knockout cells to facilitate subsequent mouse development. We also demonstrate the generality of our strategy by engineering a conditional knockout of endogenous Hprt. Taken together, our work provides a new tool for advanced mouse genome writing and mouse model establishment.
Overwriting counterselectable markers is an efficient strategy for removing wild-type DNA or replacing it with payload DNA of interest. Currently, one bottleneck of efficient genome engineering in mammals is the shortage of counterselectable (negative selection) markers that work robustly without affecting organismal developmental potential. Here, we report a conditional Piga knockout strategy that enables efficient proaerolysin-based counterselection in mouse embryonic stem cells. The conditional Piga knockout cells show similar proaerolysin resistance as full (non-conditional) Piga deletion cells, which enables the use of a PIGA transgene as a counterselectable marker for genome engineering purposes. Native Piga function is readily restored in conditional Piga knockout cells to facilitate subsequent mouse development. We also demonstrate the generality of our strategy by engineering a conditional knockout of endogenous Hprt. Taken together, our work provides a new tool for advanced mouse genome writing and mouse model establishment.
Mouse models have long served as key tools for studying disease mechanisms. Advances in stem cell manipulation and genome engineering have facilitated the development of many useful genetically engineered mouse models (GEMMs). In recent decades, as biologists advanced the ability to synthesize and integrate megabase-sized DNA sequences of interest across systems from Mycoplasma subspecies to Saccharomyces cerevisiae (Gibson et al., 2008, 2010; Dymond et al., 2011; Annaluru et al., 2014; Richardson et al., 2017; Shen et al., 2017; Xie et al., 2017; Mitchell et al., 2017; Wu et al., 2017; Zhang et al., 2017, 2020), increased attention has turned to how such technologies could be applied to mouse genome writing. Large-scale mouse genome engineering on the order of hundreds of kilobases (kb) at a time makes it possible to design and integrate complex sequences that would be difficult to complete using one-edit-at-a-time approaches. Additionally, fully overwriting a few mouse genes with their human counterparts allows for better recapitulation of human gene functions in mice (Lee et al., 2014; Mitchell et al., 2021; Wallace et al., 2007). At present, challenges in mouse genome writing have mainly focused on delivering large pieces of DNA into mouse embryonic stem cells (mESCs) in a site-specific manner, which enables facile analysis of the resultant animals.Existing big DNA delivery technologies in non-mammalian organisms provide useful insights for mouse genome writing. For example, the Synthetic Yeast Genome Project (Sc2.0) uses the Switch Auxotrophies Progressively for Integration (SwAP-In) method to integrate large sequences of interest in a stepwise manner. In more detail, SwAP-In alternates auxotrophic markers during each engineering step to select for insertion of synthetic DNA and against the preexisting wild-type (WT) sequences (Dymond and Boeke, 2012; Richardson et al., 2017). A similar approach relying on both selection and counterselection was deployed for the assembly of the 61-codon synthetic Escherichia coli genome (Wang et al., 2016; Fredens et al., 2019). In mESCs, many positive selection markers, including the puromycin-N-acetyltransferase gene (pac), blasticidin S-resistance gene (bsr), and hygromycin phosphotransferase gene (hph), are readily available. However, effective counterselectable markers are generally in short supply, limiting the development of efficient DNA delivery methods for mESCs.Recombinase-mediated cassette exchange has proven promising for targeted DNA delivery to cells, and has been adapted for delivery of payloads larger than 100 kb (Iacovino et al., 2011; Lee et al., 2014; Wallace et al., 2007). Recently, the Big-IN method, which relies on overwriting an intermediate landing pad harboring a counterselectable marker, was shown to enable site-specific >100 kb payload delivery, leaving only recombinase binding sites as genomic scars (Brosh et al., 2021). Additionally, if iteratively delivering DNA to mESCs, counterselectable markers can be embedded near the end of the payload DNA and used for selecting against residual unwanted DNA in the next engineering step. However, many common counterselectable markers used in mESCs have significant limitations. One popular counterselectable marker, the herpes simplex virus type 1–thymidine kinase (HSV1-TK, or TK), produces a bystander effect by which toxic ganciclovir-triphosphate can diffuse from TK positive cells to TK negative cells via gap junctions, resulting in unwanted TK negative cell death following ganciclovir treatment and limiting TK’s counterselection accuracy (Brosh et al., 2021; Elshami et al., 1996; Mesnil et al., 1996). Diphtheria toxin binds to human heparin-binding EGF-like growth factor (HB-EGF) and abolishes host protein synthesis. Introducing the diphtheria toxin receptor (DTR) in mice allows targeted cell ablation in transgenic mice (Saito et al., 2001), but it cannot be used as a counterselectable marker in human cells. Other counterselectable markers rely on knocking out important endogenous genes, which may require additional engineering steps to reconstitute the knocked-out gene using a minigene before phenotyping the engineered cells or converting mESCs to mice. Neither of these strategies for reintroducing genes is ideal; indeed, minigenes lack the natural regulatory sequences associated with their native loci of interest. Moreover, their reintroduction can be technically challenging and can lead to additional off-target genome modifications.One such option for counterselection in mESCs is Piga (phosphatidylinositol glycan anchor biosynthesis class A), an essential gene in the glycosylphosphatidylinositol (GPI)-anchor biosynthesis pathway (Miyata et al., 1993). Piga-expressing cells are sensitive to proaerolysin, a bacterial prototoxin that binds cell-surface GPI anchors and perforates the cell membrane, causing rapid cell death (Abrami et al., 1998). Knocking out the single copy of the X-linked Piga gene in male genomes readily generates proaerolysin-resistant cells. This enables the use of a Piga minigene as a counterselectable marker on intermediate landing pads and consequently, efficient isolation of cells in Piga-overwritten cells following genome writing with no bystander effect (Brosh et al., 2021; Li et al., 2021). However, as PIGA is essential for murine development (Rosti et al., 1997), PIGA function must be restored before mESCs can give rise to mice.Here, we demonstrate conditionally knocking out and restoring Piga using a floxed STOP cassette inserted into the second intron of Piga. Using this strategy, we induce a temporal endogenous Piga knockout state, enabling the use of a Piga transgene as a counterselectable marker for genome writing. We demonstrate highly efficient excision of the STOP cassette using Cre recombinase, which can be applied once mESC engineering is complete to enable mouse development. This easily deployable “counterselection switch” preserves the native regulatory structure (and thus native levels and splice isoforms of the gene in multiple tissues) of Piga. This approach provides an efficient counterselectable marker for mESCs engineering, while preserving the potential of developing GEMMs from engineered mESCs that contain large, targeted DNA modifications of interest.
Results
Construction of a reversible Piga knockout mouse embryonic stem cell lines
To conditionally knock out Piga in a male mESC line, we inserted an STOP cassette consisting of a β-geo (β-galactosidase and neomycin fusion) marker and three tandemly repeated SV40 polyadenylation signals (pA) into the second intron of the Piga gene to prematurely terminate its transcription (Figure 1A). The STOP cassette is flanked by two forward-facing loxP sites to enable Cre-mediated excision. Integration of the STOP cassette was achieved by co-transfecting mESCs with a plasmid co-expressing Cas9 and a guide RNA (gRNA) targeting the second intron of Piga (chrX:164,425,912, mm10) along with a plasmid carrying the STOP cassette flanked by homology arm (HA) sequences (Figure S1A). Successful STOP cassette insertion (henceforth referred to as Piga-STOP) was initially verified by identifying the novel genomic junctions by PCR in G418-resistant clones (Figure 1B). Two candidate clones that passed the junction PCR screening were subjected to targeted capture sequencing for further verification of engineering precision (Brosh et al., 2021). We used a modular mapping approach, in which the sequencing reads were separately mapped to three references: the STOP cassette, the plasmid backbone of the STOP cassette, and the mouse genome (mm10). We found that both Piga-STOP clones had reads mapping to the STOP cassette (Figure 1C), indicating that the STOP cassette was successfully inserted into the genome. We then applied bamintersect integration site analysis (Brosh et al., 2021), which identifies reads mapping to two different reference genomes, enabling unbiased detection of the STOP cassette’s integration site. This analysis detected no off-target junction for the two sequenced clones (Table S1). We also found both Piga-STOP clones had reads mapping to the plasmid backbone of the STOP cassette (Figure 1D). Also, a 2-fold increase in sequencing coverage in the Piga homology arm regions relative to the surrounding genome coverage was observed (Figure S1B), indicating that two copies of the STOP cassette and a copy of the plasmid backbone might have integrated. Genotyping PCR spanning the backbone and STOP cassette confirmed this duplicated configuration (Figure S1C). Nonetheless, the duplicated STOP cassette should still fully retain the conditional knockout’s reversibility as the entire insertion is still flanked by two loxP sites.
Figure 1
Piga transcriptional knockout by a floxed STOP cassette
(A) Schematic for Piga conditional knockout. Dark blue lines, Piga introns, dark blue boxes, coding exons or untranslated region (UTR). L-HA and R-HA, left and right homology arms (∼900 bp), respectively. pCAG, synthetic promoter consists of CMV enhancer, chicken beta-Actin promoter and rabbit beta-Globin splice acceptor site. Arrows indicate the position of genotyping primers. Scissors indicate the Cas9-gRNA cutting site.
(B) Genotyping PCR of STOP cassette insertion clones. L, left junction (primers P1 and P2); R, right junction (primers P3 and P4). Ladder, 1 kb plus DNA ladder (NEB).
(C) Capture sequencing analysis of Piga-STOP mESC clones. Coverage plots of reads mapping to STOP cassette.
(D) Capture sequencing analysis of Piga-STOP mESC clones. Coverage plots of reads mapping to the STOP cassette backbone.
Piga transcriptional knockout by a floxed STOP cassette(A) Schematic for Piga conditional knockout. Dark blue lines, Piga introns, dark blue boxes, coding exons or untranslated region (UTR). L-HA and R-HA, left and right homology arms (∼900 bp), respectively. pCAG, synthetic promoter consists of CMV enhancer, chicken beta-Actin promoter and rabbit beta-Globin splice acceptor site. Arrows indicate the position of genotyping primers. Scissors indicate the Cas9-gRNA cutting site.(B) Genotyping PCR of STOP cassette insertion clones. L, left junction (primers P1 and P2); R, right junction (primers P3 and P4). Ladder, 1 kb plus DNA ladder (NEB).(C) Capture sequencing analysis of Piga-STOP mESC clones. Coverage plots of reads mapping to STOP cassette.(D) Capture sequencing analysis of Piga-STOP mESC clones. Coverage plots of reads mapping to the STOP cassette backbone.
Evaluation of Piga-STOP mESCs
To evaluate whether transcription of Piga can be terminated by the STOP cassette, we examined the mRNA level of Piga using primers targeting exons downstream of the STOP cassette integration site. As expected, Piga transcription was significantly decreased in Piga-STOP cells compared to wild-type cells (Figure 2A). Consistent with this, PIGA protein was significantly depleted in Piga-STOP cells as determined by western blotting (Figure 2B).
Figure 2
Evaluation of proaerolysin resistance in Piga-STOP mESCs
(A) RT-qPCR analysis of Piga mRNA expression in WT and Piga-STOP mESC lines. Values are normalized to Actb mRNA. Bars represent mean ± SD of four technical replicates.
(B) Western blot analysis of PIGA in wild-type and Piga-STOP mESCs. α-tubulin serves as a loading control.
(C) Illustration of proaerolysin-mediated cell death; images were generated with bioRender.
(D) Bright-field images of mESCs after 24 h of 5 nM proaerolysin treatment.
(E) Cell viability comparison between proaerolysin-treated and untreated mESCs (5 nM proaerolysin for 24 h). Cell viability was quantified using PrestoBlue (see STAR Methods). Bars represent mean ± SD of two replicates.
Evaluation of proaerolysin resistance in Piga-STOP mESCs(A) RT-qPCR analysis of Piga mRNA expression in WT and Piga-STOP mESC lines. Values are normalized to Actb mRNA. Bars represent mean ± SD of four technical replicates.(B) Western blot analysis of PIGA in wild-type and Piga-STOP mESCs. α-tubulin serves as a loading control.(C) Illustration of proaerolysin-mediated cell death; images were generated with bioRender.(D) Bright-field images of mESCs after 24 h of 5 nM proaerolysin treatment.(E) Cell viability comparison between proaerolysin-treated and untreated mESCs (5 nM proaerolysin for 24 h). Cell viability was quantified using PrestoBlue (see STAR Methods). Bars represent mean ± SD of two replicates.Proaerolysin-mediated cell death is initiated when the active form of the proaerolysin monomer binds to the glycosylphosphatidylinositol (GPI) anchors on the cell membrane. The activation of proaerolysin leads to a conformational change followed by the formation of a heptameric oligomer. The heptamer inserts into cell membrane, forming a membrane pore that leads to rapid cell death (Iacovache et al., 2011, 2016) (Figure 2C). Blocking Piga transcriptionally should abolish GPI-anchor biosynthesis, thus preventing proaerolysin from binding GPI-anchored proteins and rendering Piga-STOP mESCs resistant to proaerolysin. As a control, we used a complete Piga deletion (ΔPiga) mESC line (Brosh et al., 2021). A proaerolysin kill curve assay was first performed to determine the lowest effective working concentration. We found that 0.5 nM and higher concentrations of proaerolysin eliminated almost all the sensitive mESCs (Figure S2A). To make sure no sensitive cells escaped proaerolysin-mediated cell death, we withdrew proaerolysin after 24 h and continued to culture any residual cells for three days. Consistent with our initial observations, concentrations of 0.5 nM proaerolysin and above produced no escapees (Figure S2B). Next, we tested the resistance of Piga-STOP mESCs to proaerolysin. Severe cell shrinkage was observed in wild-type mESCs during proaerolysin treatment, suggesting loss of cell content due to cell membrane perforation. By contrast, both ΔPiga and Piga-STOP mESCs showed no obvious morphological changes following proaerolysin treatment (Figure 2D). We then evaluated cell viability in proaerolysin-treated and untreated cells. The viability ratio between proaerolysin-treated and untreated cells reflects the degree of resistance to proaerolysin (with zero indicating full proaerolysin sensitivity and one indicating complete resistance). The Piga-STOP lines displayed nearly complete resistance to proaerolysin in comparison to the ΔPiga mESCs (Figure 2E), consistent with functional inactivation of PIGA by the STOP cassette.
The human PIGA minigene serves as a counterselectable marker in Piga-STOP mESCs
Next, we asked whether Piga-STOP mESCs can serve as a clean genetic background for proaerolysin-based counterselection. We first restored proaerolysin sensitivity by integrating a human PIGA transgene (hPIGA) -containing marker cassette into the mouse Hprt locus of Piga-STOP mESCs (Figure 3A). Cells were selected for the presence of the marker cassette with puromycin and for the loss of HPRT activity using 6-thioguanine (6-TG). The integration was validated by the presence of novel junctions between the integrated hPIGA marker cassette and Hprt flanking regions (Figure S2C, left). The presence of the STOP cassette in Piga was confirmed in these clones (Figure S2C, right). Capture sequencing of two Piga-STOP-hPIGA clones confirmed on-target cassette integration and no plasmid backbone integration (Figures 3B, 3C, S2D, and Table S1). To test whether we can select for hPIGA loss in the Piga-STOP background, we subsequently disrupted hPIGA by introducing two Cas9-gRNAs into Piga-STOP-hPIGA mESCs. We found that following “no DNA control” transfections, both Piga-STOP-hPIGA lines produced zero colonies after proaerolysin selection, suggesting the hPIGA transgene renders mESCs fully sensitive to proaerolysin. In contrast, following the targeting of the hPIGA marker cassette by Cas9-gRNAs, hundreds of colonies were resistant to proaerolysin (Figure 3D). Furthermore, Sanger sequencing revealed that at least one of the two gRNA sites had indels (Figure 3E).
Figure 3
hPIGA integration and deletion
(A) Schematic for integration of the hPIGA marker cassette into the Hprt locus.
(B) Capture sequencing of two Piga-STOP-hPIGA clones. Reads were mapped to mm10 (top) or the marker cassette (bottom).
(C) Bamintersect integration site analysis of junctions spanning hPIGA marker cassette and mm10.
(D) Crystal violet staining of Piga-STOP-hPIGA mESC colonies with or without hPIGA targeting Cas9-gRNAs.
(E) Sanger sequencing of hPIGA gRNA sites in Cas9-gRNA-transfected clones.
hPIGA integration and deletion(A) Schematic for integration of the hPIGA marker cassette into the Hprt locus.(B) Capture sequencing of two Piga-STOP-hPIGA clones. Reads were mapped to mm10 (top) or the marker cassette (bottom).(C) Bamintersect integration site analysis of junctions spanning hPIGA marker cassette and mm10.(D) Crystal violet staining of Piga-STOP-hPIGA mESC colonies with or without hPIGA targeting Cas9-gRNAs.(E) Sanger sequencing of hPIGA gRNA sites in Cas9-gRNA-transfected clones.
Piga restoration by efficient removal of STOP cassette
Given the essential role of PIGA in murine development, Piga expression necessitates restoration before generating viable mice from Piga-STOP mESCs. To facilitate this, we transiently expressed Cre recombinase in Piga-STOP mESCs to excise the STOP cassette (Figure 4A). mESCs were sparsely plated after transfecting Cre-expressing plasmid, allowing single cell-derived colony formation. We randomly selected 94 mESC colonies for PCR genotyping, and 58.5% of the colonies showed deletion of the STOP cassette (Figures 4B, S3A, and S3B). These results demonstrate that the removal of the STOP cassette is highly efficient, despite the extra integrated copy of the STOP cassette and its plasmid backbone (Figure 1D). Capture sequencing of Piga-restored clones revealed that the duplicated homology arm regions were excised together with the STOP cassette, as expected (Figure S3C), and Sanger sequencing of the excision-spanning PCR amplicons confirmed the expected 76 bp scar (Figure S3D). Robust Piga transcription was detected in STOP cassette removal lines by RT-qPCR (Figure 4C). Importantly, we found reacquisition of proaerolysin sensitivity in STOP cassette-excised mESCs (Figures 4D and 4E).
Figure 4
Piga restoration by Cre-mediated STOP cassette excision
(A) Schematic diagram showing Piga gene before and after STOP cassette excision.
(B) PCR genotyping of WT, Piga-STOP, and Piga-restoration clones. L, left junction of the STOP cassette; R, right junction of the STOP cassette; F, STOP excision. The expected size of the F product in pre-excision cells is 6.4 kb, which is challenging to amplify in this PCR assay. Ladder, 1 kb plus DNA ladder (NEB).
(C) RT-qPCR analysis of Piga expression in Piga-restored clones. Fold change was calculated by normalizing to WT. Bars represent mean ± SD of four technical replicates.
(D) Bright field view of mESCs 24 h after proaerolysin (5 nM) treatment.
(E) Cell viability comparison (measured with PrestoBlue) between proaerolysin-treated and untreated mESCs. Bars represent mean ± SD of three replicates.
Piga restoration by Cre-mediated STOP cassette excision(A) Schematic diagram showing Piga gene before and after STOP cassette excision.(B) PCR genotyping of WT, Piga-STOP, and Piga-restoration clones. L, left junction of the STOP cassette; R, right junction of the STOP cassette; F, STOP excision. The expected size of the F product in pre-excision cells is 6.4 kb, which is challenging to amplify in this PCR assay. Ladder, 1 kb plus DNA ladder (NEB).(C) RT-qPCR analysis of Piga expression in Piga-restored clones. Fold change was calculated by normalizing to WT. Bars represent mean ± SD of four technical replicates.(D) Bright field view of mESCs 24 h after proaerolysin (5 nM) treatment.(E) Cell viability comparison (measured with PrestoBlue) between proaerolysin-treated and untreated mESCs. Bars represent mean ± SD of three replicates.
Applying the STOP cassette to Hprt
To demonstrate the general applicability of the conditional knockout strategy to other counterselectable markers in mESCs, we next sought to introduce the STOP cassette into the mouse hypoxanthine phosphoribosyltransferase (Hprt) gene. The deletion of Hprt renders cells resistant to 6-thioguanine (6-TG), a purine analog that is converted by HPRT to the toxin 6-thioguanosine monophosphate. We determined the 6-TG concentration that eliminates all escapees (Figure S4A). Similar to the approach we used for Piga-STOP engineering, we precisely inserted the STOP cassette into the first intron of Hprt (chrX:52,992,027, mm10) (Figures 5A, S4B, and Table S1). A control 6-TG-resistant mESC line (ΔHprt) was generated by a complete Hprt knockout using paired CRISPR-Cas9 deletion (Figure S4C). We then treated WT, Hprt-STOP, and ΔHprt mESCs with 6-TG, and found the Hprt-STOP mESCs displayed the same level of 6-TG resistance as ΔHprt cells (Figure 5B), suggesting that the STOP cassette-mediated gene inactivation effectively eliminated Hprt expression.
Figure 5
Engineering counterselectable Hprt knockout by STOP cassette
(A) Schematic of Hprt-STOP engineering.
(B) Cell viability comparison (measured with PrestoBlue) between 6-TG treated and untreated mESCs. Bars represent mean ± SD of three replicates.
Engineering counterselectable Hprt knockout by STOP cassette(A) Schematic of Hprt-STOP engineering.(B) Cell viability comparison (measured with PrestoBlue) between 6-TG treated and untreated mESCs. Bars represent mean ± SD of three replicates.
Discussion
Advances in genome writing demand the development of new tools to bridge the gap between DNA delivery to mESCs and eventual mouse model generation. Here, we have created a conditional Piga knockout in mESCs to enable the utilization of Piga as a counterselectable marker during genetic engineering. Piga-STOP mESCs demonstrate abrogated Piga mRNA levels and acquired proaerolysin resistance, similar to full Piga deletion lines. Following transient expression of Cre recombinase, we observed restoration of mRNA levels and proaerolysin sensitivity, consistent with normal PIGA function. Together, our findings suggest that this conditional counterselectable marker can be readily introduced to and removed from mESC genomes during the course of mESC engineering.This approach represents a proof of concept for increasing the flexibility of mESC engineering. While our system does leave behind a modest-sized intronic scar (76 bp), theory and prior work suggest that small modifications to introns are unlikely to affect gene function (Lee and Rio, 2015; Yoshimatsu and Nagawa, 1989). Existing delivery technologies, such as inducible cassette exchange (ICE), routinely leave scars following DNA integration, and the resulting cell lines are able to generate mice (Iacovino et al., 2011; Wallace et al., 2007). Thus, we expect the restored gene to behave normally and support organismal development.More broadly, we hope to adapt this approach to work with a variety of DNA delivery approaches. LoxP sites could be replaced with FRT sites or other loxP variants to improve compatibility with the Big-IN delivery system (Brosh et al., 2021), or with other DNA integration strategies that already use loxP sites for Cre-mediated recombination. We could further increase system flexibility by adding an inducible Cre segment within the STOP cassette, such as Cre-ERT2 (Indra et al., 1999). Subsequent tamoxifen treatment would trigger cassette removal and restoration of gene function, without an additional Cre recombinase introduction step. This approach is also feasible for studying noncoding RNA genes to achieve temporal control.Conditional counterselectable markers may also prove valuable for genome engineering efforts in other species, including rats, flies, and other commonly used model organisms (Murata et al., 2012). As in mice, such tools may help strike a balance between ease of genome editing and downstream organismal viability. Additionally, conditional counterselectable markers could be used in other cell lines, such as human-induced pluripotent stem cells (iPSCs). While such cell lines may not give rise to full organisms, there may be benefits to restoring the natural regulatory environment surrounding counterselectable markers following cellular engineering, particularly when studying genes or pathways related to the counterselectable marker of interest.As mammalian genome writing efforts expand, demand for tools to accelerate the pipeline from design and DNA synthesis to downstream applications will likely continue to grow. Ultimately, conditional counterselection systems like the one piloted here may become key components of the genome writing toolbox, accelerating our ability to study basic biology and model diseases.
Limitations of study
Here, we present a conditional knockout of two X-linked endogenous genes, Piga and Hprt, to enable proaerolysin and 6-TG-based counterselections in a male mESC line. We have discussed some limitations of this study, including the compatibility to recombinase-based big DNA delivery methods. The insertion and excision of the STOP cassette is simplified due to the single copy nature of X chromosome in male mESC. Biallelic insertion of the STOP cassette is required when implementing this strategy to autosomal genes, thus complicating the engineering process. Further optimization of the STOP cassette by introduction of a positive selection marker other from NeoR (for example PuroR), co-delivering both STOP cassettes with distinct selection markers will facilitate the selection of biallelic insertions.
STAR★Methods
Key resources table
Resource availability
Lead contact
Further information should be directed to and will be fulfilled by the lead contact Jef D. Boeke (jef.boeke@nyulangone.org).
Materials availability
All unique/stable reagents generated in this study are available from the Lead contact with a completed Materials Transfer Agreement.
Experimental models and subject details
mESCs
C57BL/6J male mESCs (MK6) were provided by NYU Langone Health, Rodent Genetic Engineering Laboratory. mESCs were cultured in 80/20 medium [80% of 2i medium (2i basal medium supplement with 3 μM CHIR99021 and 1 μM PD0325901) mixed 20% of ES medium [Knockout DMEM, 15% fetal bovine serum, 1% MEM Non-Essential Amino Acids, 103 U/mL Leukemia Inhibitory Factor, 0.1 mM 2-mercaptoethanol, 1% Pen-Strep]) on plates coated with 0.1% gelatin (EMD Milipore, ES-006-B). mESCs were grown in a humidified tissue culture incubator at 37°C, 5% CO2. Medium was exchanged daily.
Method details
Plasmid construction
The STOP cassette plasmid is derived from pWA203 used in a previous study (Han et al., 2008). The STOP cassette fragment was released from pWA203 by AscI and KpnI digestion. Piga and Hprt homology arms were amplified from mouse BACs RP23-32H22 and RP23-412J16 using primers shown in Table S2. Gibson assembly was used to construct the STOP cassette with flanking homology arms into the multiple cloning site of pRS413 vector. Puro-hPIGA-mScarlet plasmid was constructed by stitching all the components together using Gibson assembly enzyme mix. Guide RNAs (gRNAs, listed in Table S3) were chosen using the CRISPOR (crispor.tefor.net/) online tool, and subsequently cloned into pX330 (Addgene, 42230) using Golden Gate assembly (Ran et al., 2013). The pCAG-iCre plasmid was purchased from Addgene (89573).
Nucleofection
The mouse embryonic stem cell nucleofector kit (Lonza, VPH-1001) was used for the nucleofection. mESCs were harvested by trypsinization followed by a PBS wash and were counted using a hemocytometer. Approximately 2 × 106 cells were used for each reaction. A total of 10 μg DNA was mixed with mESCs prior to the nucleofection. The nucleofection was conducted using the Lonza 2b Nucleofector, program A-023. mESCs were resuspended in culture medium after the nucleofection and plated immediately onto gelatin-coated 10 cm plates.
Cell viability assay
Approximately 1 × 105 mESCs were seeded in a 24 well plate. After growing for 24 h, mESCs were treated either with 5 nM of proaerolysin for 24 h or 5 μM of 6-thioguanine for 3 days. Both treated and untreated cells were harvested by trypsinization. mESCs were resuspended in 1 mL of 80/20 medium, and 180 μL cell suspension was mixed with 20 μL of PrestoBlue reagent (Invitrogen, A13261) in a 96 well plate for 30 min at 37°C, 5% CO2. Absorbance was read using a plate reader following manufacturer’s instructions (Biotek, synergy H1).
RT-qPCR
mESCs were harvested by trypsinization. Total RNA was isolated using the RNeasy kit (QIAGEN, 74106) following the manufacturer’s instructions. 1 μg of total RNA was used for reverse transcription (Invitrogen, 18090200). 1 μL of 10-fold diluted cDNA was used as template in a 10 μL SYBR (Roche, 04887352001) qPCR reaction on a LightCycler 480 instrument. Relative expression was calculated using the ΔΔCT method.
Immunoblotting
mESCs were harvested by trypsinization and pellets were resuspended in 200 μL RIPA buffer supplemented with protease inhibitor cocktail (Roche, 11873580001). Cell suspensions were then incubated on ice for 30 min. After the on-ice incubation, 1 μL Benzonase (Sigma-Aldrich, E1014) was added to the cell lysate for a 10 min 37°C digestion of genomic DNA. Protein lysates were mixed with LDS loading buffer (Invitrogen, NP0007), and boiled at 70°C for 10 min. Proteins were separated using a NuPAGE 4–12% Bis-Tris gel (Life Technologies, WG1402BOX) and transferred to PVDF membranes (Millipore, IPFL00010), which were blotted with PIGA (Abcam, ab69768) and Tubulin (Sigma-Aldrich, T5168) primary antibodies, followed by goat anti-mouse (LI-COR Biosciences, 926-32210) and goat anti-rabbit (LI-COR Biosciences, 926-68071) secondary antibodies. Blots were visualized using the LI-COR system.
Targeted capture sequencing
Targeted capture sequencing was performed as previously described (Brosh et al., 2021). In brief, mESC genomic DNA was isolated with QIAamp gDNA extraction kit (QIAGEN, 51306) following the manufacturer’s instructions. 1 μg of genomic DNA was used as input for the large fragment size (>550 bp) DNA library preparation (NEBNext Ultra II FS). Mouse BACs RP23-32H22 and RP23-412J16 (BACPAC Resources Center), STOP cassette plasmid, Puro-hPIGA-mScarlet and Cas9-gRNA plasmids were used for nick translation to generate the biotinylated bait set. Captured DNA library was sequenced using a NextSeq 500 75 cycles high output kit (Illumina, 20024906).
Quantification and statistical analysis
For RT-qPCR assay, data were collected from four technical replicates. ΔΔCT method was used for computing the relative expression level. For cell viability assays, two (Figure 2E) or three (Figures 4E and 5B) wells of cells were used as replicates. All graphs were created in GraphPad Prism 9.
REAGENT or RESOURCE
SOURCE
IDENTIFIER
Antibodies
Anti-PIGA antibody
Abcam
Cat# ab69768, RRID: AB_1566593
Mouse monoclonal anti-alpha-Tubulin antibody
Sigma-Aldrich
Cat# T5168, RRID: AB_477579
Chemicals, peptides, and recombinant proteins
Crystal violet solution
Sigma-Aldrich
V5265-500ML
LightCycler 480 SYBR green I master
Roche
04887352001
Geneticin (G418 sulfate)
Gibco
10131027
Critical commercial assays
PrestoBlue cell viability assay
Invitrogen
A13261
Mouse ES cell nucleofection kit
Lonza
VPH-1001
NEBNext Ulta II FS DNA library prep kit
NEB
E7645S
SuperScript IV reverse transcriptase
Invitrogen
18090200
NextSeq High-output 75-cycle V2.5 Kit
Illumina
20024906
Experimental models: Cell lines
C57BL/6J (MK6) mouse embryonic stem cells
NYU Langone Health Rodent Genetic Engineering Laboratory
Authors: Shane Neph; M Scott Kuehn; Alex P Reynolds; Eric Haugen; Robert E Thurman; Audra K Johnson; Eric Rynes; Matthew T Maurano; Jeff Vierstra; Sean Thomas; Richard Sandstrom; Richard Humbert; John A Stamatoyannopoulos Journal: Bioinformatics Date: 2012-05-09 Impact factor: 6.937
Authors: M Saito; T Iwawaki; C Taya; H Yonekawa; M Noda; Y Inui; E Mekada; Y Kimata; A Tsuru; K Kohno Journal: Nat Biotechnol Date: 2001-08 Impact factor: 54.908
Authors: F Ann Ran; Patrick D Hsu; Jason Wright; Vineeta Agarwala; David A Scott; Feng Zhang Journal: Nat Protoc Date: 2013-10-24 Impact factor: 13.491
Authors: Daniel G Gibson; Gwynedd A Benders; Cynthia Andrews-Pfannkoch; Evgeniya A Denisova; Holly Baden-Tillson; Jayshree Zaveri; Timothy B Stockwell; Anushka Brownley; David W Thomas; Mikkel A Algire; Chuck Merryman; Lei Young; Vladimir N Noskov; John I Glass; J Craig Venter; Clyde A Hutchison; Hamilton O Smith Journal: Science Date: 2008-01-24 Impact factor: 47.728
Authors: Sarah M Richardson; Leslie A Mitchell; Giovanni Stracquadanio; Kun Yang; Jessica S Dymond; James E DiCarlo; Dongwon Lee; Cheng Lai Victor Huang; Srinivasan Chandrasegaran; Yizhi Cai; Jef D Boeke; Joel S Bader Journal: Science Date: 2017-03-10 Impact factor: 47.728
Authors: Leslie A Mitchell; Ann Wang; Giovanni Stracquadanio; Zheng Kuang; Xuya Wang; Kun Yang; Sarah Richardson; J Andrew Martin; Yu Zhao; Roy Walker; Yisha Luo; Hongjiu Dai; Kang Dong; Zuojian Tang; Yanling Yang; Yizhi Cai; Adriana Heguy; Beatrix Ueberheide; David Fenyö; Junbiao Dai; Joel S Bader; Jef D Boeke Journal: Science Date: 2017-03-10 Impact factor: 47.728
Authors: Donghui Li; Xiaoji Sun; Fangzhou Yu; Mary Ann Perle; David Araten; Jef D Boeke Journal: Nucleic Acids Res Date: 2021-03-18 Impact factor: 16.971