Precise therapeutic gene correction by a simple nuclease-induced double-stranded break.

Current programmable nuclease-based methods (for example, CRISPR-Cas9) for the precise correction of a disease-causing genetic mutation harness the homology-directed repair pathway. However, this repair process requires the co-delivery of an exogenous DNA donor to recode the sequence and can be inefficient in many cell types. Here we show that disease-causing frameshift mutations that result from microduplications can be efficiently reverted to the wild-type sequence simply by generating a DNA double-stranded break near the centre of the duplication. We demonstrate this in patient-derived cell lines for two diseases: limb-girdle muscular dystrophy type 2G (LGMD2G)1 and Hermansky-Pudlak syndrome type 1 (HPS1)2. Clonal analysis of inducible pluripotent stem (iPS) cells from the LGMD2G cell line, which contains a mutation in TCAP, treated with the Streptococcus pyogenes Cas9 (SpCas9) nuclease revealed that about 80% contained at least one wild-type TCAP allele; this correction also restored TCAP expression in LGMD2G iPS cell-derived myotubes. SpCas9 also efficiently corrected the genotype of an HPS1 patient-derived B-lymphoblastoid cell line. Inhibition of polyADP-ribose polymerase 1 (PARP-1) suppressed the nuclease-mediated collapse of the microduplication to the wild-type sequence, confirming that precise correction is mediated by the microhomology-mediated end joining (MMEJ) pathway. Analysis of editing by SpCas9 and Lachnospiraceae bacterium ND2006 Cas12a (LbCas12a) at non-pathogenic 4-36-base-pair microduplications within the genome indicates that the correction strategy is broadly applicable to a wide range of microduplication lengths and can be initiated by a variety of nucleases. The simplicity, reliability and efficacy of this MMEJ-based therapeutic strategy should permit the development of nuclease-based gene correction therapies for a variety of diseases that are associated with microduplications.

Enhanced Cas12a editing in mammalian cells and zebrafish.

Type V CRISPR-Cas12a systems provide an alternate nuclease platform to Cas9, with potential advantages for specific genome editing applications. Here we describe improvements to the Cas12a system that facilitate efficient targeted mutagenesis in mammalian cells and zebrafish embryos. We show that engineered variants of Cas12a with two different nuclear localization sequences (NLS) on the C terminus provide increased editing efficiency in mammalian cells. Additionally, we find that pre-crRNAs comprising a full-length direct repeat (full-DR-crRNA) sequence with specific stem-loop G-C base substitutions exhibit increased editing efficiencies compared with the standard mature crRNA framework. Finally, we demonstrate in zebrafish embryos that the improved LbCas12a and FnoCas12a nucleases in combination with these modified crRNAs display high mutagenesis efficiencies and low toxicity when delivered as ribonucleoprotein complexes at high concentration. Together, these results define a set of enhanced Cas12a components with broad utility in vertebrate systems.

Efficient Homology-Directed Repair with Circular Single-Stranded DNA Donors.

While genome editing has been revolutionized by the advent of CRISPR-based nucleases, difficulties in achieving efficient, nuclease-mediated, homology-directed repair (HDR) still limit many applications. Commonly used DNA donors such as plasmids suffer from low HDR efficiencies in many cell types, as well as integration at unintended sites. In contrast, single-stranded DNA (ssDNA) donors can produce efficient HDR with minimal off-target integration. In this study, we describe the use of ssDNA phage to efficiently and inexpensively produce long circular ssDNA (cssDNA) donors. These cssDNA donors serve as efficient HDR templates when used with Cas9 or Cas12a, with integration frequencies superior to linear ssDNA (lssDNA) donors. To evaluate the relative efficiencies of imprecise and precise repair for a suite of different Cas9 or Cas12a nucleases, we have developed a modified traffic light reporter (TLR) system (TLR-multi-Cas variant 1 [MCV1]) that permits side-by-side comparisons of different nuclease systems. We used this system to assess editing and HDR efficiencies of different nuclease platforms with distinct DNA donor types. We then extended the analysis of DNA donor types to evaluate efficiencies of fluorescent tag knockins at endogenous sites in HEK293T and K562 cells. Our results show that cssDNA templates produce efficient and robust insertion of reporter tags. Targeting efficiency is high, allowing production of biallelic integrants using cssDNA donors. cssDNA donors also outcompete lssDNA donors in template-driven repair at the target site. These data demonstrate that circular donors provide an efficient, cost-effective method to achieve knockins in mammalian cell lines.

Allelic variation in class I HLA determines CD8+ T cell repertoire shape and cross-reactive memory responses to SARS-CoV-2.

Effective presentation of antigens by human leukocyte antigen (HLA) class I molecules to CD8+ T cells is required for viral elimination and generation of long-term immunological memory. In this study, we applied a single-cell, multiomic technology to generate a unified ex vivo characterization of the CD8+ T cell response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) across four major HLA class I alleles. We found that HLA genotype conditions key features of epitope specificity, TCRα/ß sequence diversity, and the utilization of pre-existing SARS-CoV-2-reactive memory T cell pools. Single-cell transcriptomics revealed functionally diverse T cell phenotypes of SARS-CoV-2-reactive T cells, associated with both disease stage and epitope specificity. Our results show that HLA variations notably influence the CD8+ T cell repertoire shape and utilization of immune recall upon SARS-CoV-2 infection.

Valves Are a Conserved Feature of the Zebrafish Lymphatic System.

The lymphatic system comprises blind-ended tubes that collect interstitial fluid and return it to the circulatory system. In mammals, unidirectional lymphatic flow is driven by muscle contraction working in conjunction with valves. Accordingly, defective lymphatic valve morphogenesis results in backflow leading to edema. In fish species, studies dating to the 18th century failed to identify lymphatic valves, a precedent that currently persists, raising the question of whether the zebrafish could be used to study the development of these structures. Here, we provide functional and morphological evidence of valves in the zebrafish lymphatic system. Electron microscopy revealed valve ultrastructure similar to mammals, while live imaging using transgenic lines identified the developmental origins of lymphatic valve progenitors. Zebrafish embryos bearing mutations in genes required for mammalian valve morphogenesis show defective lymphatic valve formation and edema. Together, our observations provide a foundation from which to further investigate lymphatic valve formation in zebrafish.

Highly efficient therapeutic gene editing of human hematopoietic stem cells.

Re-expression of the paralogous γ-globin genes (HBG1/2) could be a universal strategy to ameliorate the severe ß-globin disorders sickle cell disease (SCD) and ß-thalassemia by induction of fetal hemoglobin (HbF, α2γ2)1. Previously, we and others have shown that core sequences at the BCL11A erythroid enhancer are required for repression of HbF in adult-stage erythroid cells but are dispensable in non-erythroid cells2-6. CRISPR-Cas9-mediated gene modification has demonstrated variable efficiency, specificity, and persistence in hematopoietic stem cells (HSCs). Here, we demonstrate that Cas9:sgRNA ribonucleoprotein (RNP)-mediated cleavage within a GATA1 binding site at the +58 BCL11A erythroid enhancer results in highly penetrant disruption of this motif, reduction of BCL11A expression, and induction of fetal γ-globin. We optimize conditions for selection-free on-target editing in patient-derived HSCs as a nearly complete reaction lacking detectable genotoxicity or deleterious impact on stem cell function. HSCs preferentially undergo non-homologous compared with microhomology-mediated end joining repair. Erythroid progeny of edited engrafting SCD HSCs express therapeutic levels of HbF and resist sickling, while those from patients with ß-thalassemia show restored globin chain balance. Non-homologous end joining repair-based BCL11A enhancer editing approaching complete allelic disruption in HSCs is a practicable therapeutic strategy to produce durable HbF induction.

Determination of ubiquitin fitness landscapes under different chemical stresses in a classroom setting.

Ubiquitin is essential for eukaryotic life and varies in only 3 amino acid positions between yeast and humans. However, recent deep sequencing studies indicate that ubiquitin is highly tolerant to single mutations. We hypothesized that this tolerance would be reduced by chemically induced physiologic perturbations. To test this hypothesis, a class of first year UCSF graduate students employed deep mutational scanning to determine the fitness landscape of all possible single residue mutations in the presence of five different small molecule perturbations. These perturbations uncover 'shared sensitized positions' localized to areas around the hydrophobic patch and the C-terminus. In addition, we identified perturbation specific effects such as a sensitization of His68 in HU and a tolerance to mutation at Lys63 in DTT. Our data show how chemical stresses can reduce buffering effects in the ubiquitin proteasome system. Finally, this study demonstrates the potential of lab-based interdisciplinary graduate curriculum.

Systematic exploration of ubiquitin sequence, E1 activation efficiency, and experimental fitness in yeast.

The complexity of biological interaction networks poses a challenge to understanding the function of individual connections in the overall network. To address this challenge, we developed a high-throughput reverse engineering strategy to analyze how thousands of specific perturbations (encompassing all point mutations in a central gene) impact both a specific edge (interaction to a directly connected node) and an overall network function. We analyzed the effects of ubiquitin mutations on activation by the E1 enzyme and compared these to effects on yeast growth rate. Using this approach, we delineated ubiquitin mutations that selectively impacted the ubiquitin-E1 edge. We find that the elasticity function relating the efficiency of ubiquitin-E1 interaction to growth rate is non-linear and that a greater than 50-fold decrease in E1 activation efficiency is required to reduce growth rate by 2-fold. Despite the robustness of fitness to decreases in E1 activation efficiency, the effects of most ubiquitin mutations on E1 activation paralleled the effects on growth rate. Our observations indicate that most ubiquitin mutations that disrupt E1 activation also disrupt other functions. The structurally characterized ubiquitin-E1 interface encompasses the interfaces of ubiquitin with most other known binding partners, and we propose that this enables E1 in wild-type cells to selectively activate ubiquitin protein molecules capable of binding to other partners from the cytoplasmic pool of ubiquitin protein that will include molecules with chemical damage and/or errors from transcription and translation.

Resistance to vemurafenib resulting from a novel mutation in the BRAFV600E kinase domain.

Resistance to the BRAF inhibitor vemurafenib poses a significant problem for the treatment of BRAFV600E-positive melanomas. It is therefore critical to prospectively identify all vemurafenib resistance mechanisms prior to their emergence in the clinic. The vemurafenib resistance mechanisms described to date do not result from secondary mutations within BRAFV600E. To search for possible mutations within BRAFV600E that can confer drug resistance, we developed a systematic experimental approach involving targeted saturation mutagenesis, selection of drug-resistant variants, and deep sequencing. We identified a single nucleotide substitution (T1514A, encoding L505H) that greatly increased drug resistance in cultured cells and mouse xenografts. The kinase activity of BRAFV600E/L505H was higher than that of BRAFV600E, resulting in cross-resistance to a MEK inhibitor. However, BRAFV600E/L505H was less resistant to several other BRAF inhibitors whose binding sites were further from L505 than that of PLX4720. Our results identify a novel vemurafenib-resistant mutant and provide insights into the treatment for melanomas bearing this mutation.

Viewing protein fitness landscapes through a next-gen lens.

High-throughput sequencing has enabled many powerful approaches in biological research. Here, we review sequencing approaches to measure frequency changes within engineered mutational libraries subject to selection. These analyses can provide direct estimates of biochemical and fitness effects for all individual mutations across entire genes (and likely compact genomes in the near future) in genetically tractable systems such as microbes, viruses, and mammalian cells. The effects of mutations on experimental fitness can be assessed using sequencing to monitor time-dependent changes in mutant frequency during bulk competitions. The impact of mutations on biochemical functions can be determined using reporters or other means of separating variants based on individual activities (e.g., binding affinity for a partner molecule can be interrogated using surface display of libraries of mutant proteins and isolation of bound and unbound populations). The comprehensive investigation of mutant effects on both biochemical function and experimental fitness provide promising new avenues to investigate the connections between biochemistry, cell physiology, and evolution. We summarize recent findings from systematic mutational analyses; describe how they relate to a field rich in both theory and experimentation; and highlight how they may contribute to ongoing and future research into protein structure-function relationships, systems-level descriptions of cell physiology, and population-genetic inferences on the relative contributions of selection and drift.

Analyses of the effects of all ubiquitin point mutants on yeast growth rate.

The amino acid sequence of a protein governs its function. We used bulk competition and focused deep sequencing to investigate the effects of all ubiquitin point mutants on yeast growth rate. Many aspects of ubiquitin function have been carefully studied, which enabled interpretation of our growth analyses in light of a rich structural, biophysical and biochemical knowledge base. In one highly sensitive cluster on the surface of ubiquitin, almost every amino acid substitution caused growth defects. In contrast, the opposite face tolerated virtually all possible substitutions. Surface locations between these two faces exhibited intermediate mutational tolerance. The sensitive face corresponds to the known interface for many binding partners. Across all surface positions, we observe a strong correlation between burial at structurally characterized interfaces and the number of amino acid substitutions compatible with robust growth. This result indicates that binding is a dominant determinant of ubiquitin function. In the solvent-inaccessible core of ubiquitin, all positions tolerated a limited number of substitutions, with hydrophobic amino acids especially interchangeable. Some mutations null for yeast growth were previously shown to populate folded conformations indicating that, for these mutants, subtle changes to conformation caused functional defects. The most sensitive region to mutation within the core was located near the C-terminus that is a focal binding site for many critical binding partners. These results indicate that core mutations may frequently cause functional defects through subtle disturbances to structure or dynamics.

Fitness analyses of all possible point mutations for regions of genes in yeast.

Deep sequencing can accurately measure the relative abundance of hundreds of mutations in a single bulk competition experiment, which can give a direct readout of the fitness of each mutant. Here we describe a protocol that we previously developed and optimized to measure the fitness effects of all possible individual codon substitutions for 10-aa regions of essential genes in yeast. Starting with a conditional strain (i.e., a temperature-sensitive strain), we describe how to efficiently generate plasmid libraries of point mutants that can then be transformed to generate libraries of yeast. The yeast libraries are competed under conditions that select for mutant function. Deep-sequencing analyses are used to determine the relative fitness of all mutants. This approach is faster and cheaper per mutant compared with analyzing individually isolated mutants. The protocol can be performed in ∼4 weeks and many 10-aa regions can be analyzed in parallel.