Peptide-mediated delivery of CRISPR enzymes for the efficient editing of primary human lymphocytes.

CRISPR-mediated genome editing of primary human lymphocytes is typically carried out via electroporation, which can be cytotoxic, cumbersome and costly. Here we show that the yields of edited primary human lymphocytes can be increased substantially by delivering a CRISPR ribonucleoprotein mixed with an amphiphilic peptide identified through screening. We evaluated the performance of this simple delivery method by knocking out genes in T cells, B cells and natural killer cells via the delivery of Cas9 or Cas12a ribonucleoproteins or an adenine base editor. We also show that peptide-mediated ribonucleoprotein delivery paired with an adeno-associated-virus-mediated homology-directed repair template can introduce a chimaeric antigen receptor gene at the T-cell receptor α constant locus, and that the engineered cells display antitumour potency in mice. The method is minimally perturbative, does not require dedicated hardware, and is compatible with multiplexed editing via sequential delivery, which minimizes the risk of genotoxicity. The peptide-mediated intracellular delivery of ribonucleoproteins may facilitate the manufacturing of engineered T cells.

Severe Acute Respiratory Syndrome Coronavirus 2 Delta Variant Genomic Variation Associated With Breakthrough Infection in Northern California: A Retrospective Cohort Study.

BACKGROUND: The association between severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomic variation and breakthrough infection is not well defined among persons with Delta variant SARS-CoV-2 infection. METHODS: In a retrospective cohort, we assessed whether individual nonlineage defining mutations and overall genomic variation (including low-frequency alleles) were associated with breakthrough infection, defined as SARS-CoV-2 infection after coronavirus disease 2019 primary vaccine series. We identified all nonsynonymous single-nucleotide polymorphisms, insertions, and deletions in SARS-CoV-2 genomes with ≥5% allelic frequency and population frequency of ≥5% and ≤95%. Using Poisson regression, we assessed the association with breakthrough infection for each individual mutation and a viral genomic risk score. RESULTS: Thirty-six mutations met our inclusion criteria. Among 12 744 persons infected with Delta variant SARS-CoV-2, 5949 (47%) were vaccinated and 6795 (53%) were unvaccinated. Viruses with a viral genomic risk score in the highest quintile were 9% more likely to be associated with breakthrough infection than viruses in the lowest quintile, but including the risk score improved overall predictive model performance (measured by C statistic) by only +0.0006. CONCLUSIONS: Genomic variation within SARS-CoV-2 Delta variant was weakly associated with breakthrough infection, but several potential nonlineage defining mutations were identified that might contribute to immune evasion by SARS-CoV-2.

High-level correction of the sickle mutation is amplified in vivo during erythroid differentiation.

Background: A point mutation in sickle cell disease (SCD) alters one amino acid in the ß-globin subunit of hemoglobin, with resultant anemia and multiorgan damage that typically shortens lifespan by decades. Because SCD is caused by a single mutation, and hematopoietic stem cells (HSCs) can be harvested, manipulated, and returned to an individual, it is an attractive target for gene correction. Results: An optimized Cas9 ribonucleoprotein (RNP) with an ssDNA oligonucleotide donor together generated correction of at least one ß-globin allele in more than 30% of long-term engrafting human HSCs. After adopting a high-fidelity Cas9 variant, efficient correction with minimal off-target events also was observed. In vivo erythroid differentiation markedly enriches for corrected ß-globin alleles, indicating that erythroblasts carrying one or more corrected alleles have a survival advantage. Significance: These findings indicate that the sickle mutation can be corrected in autologous HSCs with an optimized protocol suitable for clinical translation.

Identification of novel HPFH-like mutations by CRISPR base editing that elevate the expression of fetal hemoglobin.

Naturally occurring point mutations in the HBG promoter switch hemoglobin synthesis from defective adult beta-globin to fetal gamma-globin in sickle cell patients with hereditary persistence of fetal hemoglobin (HPFH) and ameliorate the clinical severity. Inspired by this natural phenomenon, we tiled the highly homologous HBG proximal promoters using adenine and cytosine base editors that avoid the generation of large deletions and identified novel regulatory regions including a cluster at the -123 region. Base editing at -123 and -124 bp of HBG promoter induced fetal hemoglobin (HbF) to a higher level than disruption of well-known BCL11A binding site in erythroblasts derived from human CD34+ hematopoietic stem and progenitor cells (HSPC). We further demonstrated in vitro that the introduction of -123T > C and -124T > C HPFH-like mutations drives gamma-globin expression by creating a de novo binding site for KLF1. Overall, our findings shed light on so far unknown regulatory elements within the HBG promoter and identified additional targets for therapeutic upregulation of fetal hemoglobin.

CRISPR-Cas9-mediated knockout of CYP79D1 and CYP79D2 in cassava attenuates toxic cyanogen production.

Cassava (Manihot esculenta) is a starchy root crop that supports over a billion people in tropical and subtropical regions of the world. This staple, however, produces the neurotoxin cyanide and requires processing for safe consumption. Excessive consumption of insufficiently processed cassava, in combination with protein-poor diets, can have neurodegenerative impacts. This problem is further exacerbated by drought conditions which increase this toxin in the plant. To reduce cyanide levels in cassava, we used CRISPR-mediated mutagenesis to disrupt the cytochrome P450 genes CYP79D1 and CYP79D2 whose protein products catalyze the first step in cyanogenic glucoside biosynthesis. Knockout of both genes eliminated cyanide in leaves and storage roots of cassava accession 60444; the West African, farmer-preferred cultivar TME 419; and the improved variety TMS 91/02324. Although knockout of CYP79D2 alone resulted in significant reduction of cyanide, mutagenesis of CYP79D1 did not, indicating these paralogs have diverged in their function. The congruence of results across accessions indicates that our approach could readily be extended to other preferred or improved cultivars. This work demonstrates cassava genome editing for enhanced food safety and reduced processing burden, against the backdrop of a changing climate.

Transmission, infectivity, and neutralization of a spike L452R SARS-CoV-2 variant.

We identified an emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant by viral whole-genome sequencing of 2,172 nasal/nasopharyngeal swab samples from 44 counties in California, a state in the western United States. Named B.1.427/B.1.429 to denote its two lineages, the variant emerged in May 2020 and increased from 0% to >50% of sequenced cases from September 2020 to January 2021, showing 18.6%-24% increased transmissibility relative to wild-type circulating strains. The variant carries three mutations in the spike protein, including an L452R substitution. We found 2-fold increased B.1.427/B.1.429 viral shedding in vivo and increased L452R pseudovirus infection of cell cultures and lung organoids, albeit decreased relative to pseudoviruses carrying the N501Y mutation common to variants B.1.1.7, B.1.351, and P.1. Antibody neutralization assays revealed 4.0- to 6.7-fold and 2.0-fold decreases in neutralizing titers from convalescent patients and vaccine recipients, respectively. The increased prevalence of a more transmissible variant in California exhibiting decreased antibody neutralization warrants further investigation.

Transmission, infectivity, and antibody neutralization of an emerging SARS-CoV-2 variant in California carrying a L452R spike protein mutation.

We identified a novel SARS-CoV-2 variant by viral whole-genome sequencing of 2,172 nasal/nasopharyngeal swab samples from 44 counties in California. Named B.1.427/B.1.429 to denote its 2 lineages, the variant emerged around May 2020 and increased from 0% to >50% of sequenced cases from September 1, 2020 to January 29, 2021, exhibiting an 18.6-24% increase in transmissibility relative to wild-type circulating strains. The variant carries 3 mutations in the spike protein, including an L452R substitution. Our analyses revealed 2-fold increased B.1.427/B.1.429 viral shedding in vivo and increased L452R pseudovirus infection of cell cultures and lung organoids, albeit decreased relative to pseudoviruses carrying the N501Y mutation found in the B.1.1.7, B.1.351, and P.1 variants. Antibody neutralization assays showed 4.0 to 6.7-fold and 2.0-fold decreases in neutralizing titers from convalescent patients and vaccine recipients, respectively. The increased prevalence of a more transmissible variant in California associated with decreased antibody neutralization warrants further investigation.

Genome editing to model and reverse a prevalent mutation associated with myeloproliferative neoplasms.

Myeloproliferative neoplasms (MPNs) cause the over-production of blood cells such as erythrocytes (polycythemia vera) or platelets (essential thrombocytosis). JAK2 V617F is the most prevalent somatic mutation in many MPNs, but previous modeling of this mutation in mice relied on transgenic overexpression and resulted in diverse phenotypes that were in some cases attributed to expression level. CRISPR-Cas9 engineering offers new possibilities to model and potentially cure genetically encoded disorders via precise modification of the endogenous locus in primary cells. Here we develop "scarless" Cas9-based reagents to create and reverse the JAK2 V617F mutation in an immortalized human erythroid progenitor cell line (HUDEP-2), CD34+ adult human hematopoietic stem and progenitor cells (HSPCs), and immunophenotypic long-term hematopoietic stem cells (LT-HSCs). We find no overt in vitro increase in proliferation associated with an endogenous JAK2 V617F allele, but co-culture with wild type cells unmasks a competitive growth advantage provided by the mutation. Acquisition of the V617F allele also promotes terminal differentiation of erythroid progenitors, even in the absence of hematopoietic cytokine signaling. Taken together, these data are consistent with the gradually progressive manifestation of MPNs and reveals that endogenously acquired JAK2 V617F mutations may yield more subtle phenotypes as compared to transgenic overexpression models.

Controlled Cycling and Quiescence Enables Efficient HDR in Engraftment-Enriched Adult Hematopoietic Stem and Progenitor Cells.

Genome editing often takes the form of either error-prone sequence disruption by non-homologous end joining (NHEJ) or sequence replacement by homology-directed repair (HDR). Although NHEJ is generally effective, HDR is often difficult in primary cells. Here, we use a combination of immunophenotyping, next-generation sequencing, and single-cell RNA sequencing to investigate and reprogram genome editing outcomes in subpopulations of adult hematopoietic stem and progenitor cells. We find that although quiescent stem-enriched cells mostly use NHEJ, non-quiescent cells with the same immunophenotype use both NHEJ and HDR. Inducing quiescence before editing results in a loss of HDR in all cell subtypes. We develop a strategy of controlled cycling and quiescence that yields a 6-fold increase in the HDR/NHEJ ratio in quiescent stem cells ex vivo and in vivo. Our results highlight the tension between editing and cellular physiology and suggest strategies to manipulate quiescent cells for research and therapeutic genome editing.

ATF4 Regulates MYB to Increase γ-Globin in Response to Loss of ß-Globin.

ß-Hemoglobinopathies can trigger rapid production of red blood cells in a process known as stress erythropoiesis. Cellular stress prompts differentiating erythroid precursors to express high levels of fetal γ-globin. However, the mechanisms underlying γ-globin production during cellular stress are still poorly defined. Here, we use CRISPR-Cas genome editing to model the stress caused by reduced levels of adult ß-globin. We find that decreased ß-globin is sufficient to induce robust re-expression of γ-globin, and RNA sequencing (RNA-seq) of differentiating isogenic erythroid precursors implicates ATF4 as a causal regulator of this response. ATF4 binds within the HBS1L-MYB intergenic enhancer and regulates expression of MYB, a known γ-globin regulator. Overall, the reduction of ATF4 upon ß-globin knockout decreases the levels of MYB and BCL11A. Identification of ATF4 as a key regulator of globin compensation adds mechanistic insight to the poorly understood phenomenon of stress-induced globin compensation and could inform strategies to treat hemoglobinopathies.

The Histone Chaperone FACT Induces Cas9 Multi-turnover Behavior and Modifies Genome Manipulation in Human Cells.

Cas9 is a prokaryotic RNA-guided DNA endonuclease that binds substrates tightly in vitro but turns over rapidly when used to manipulate genomes in eukaryotic cells. Little is known about the factors responsible for dislodging Cas9 or how they influence genome engineering. Unbiased detection through proximity labeling of transient protein interactions in cell-free Xenopus laevis egg extract identified the dimeric histone chaperone facilitates chromatin transcription (FACT) as an interactor of substrate-bound Cas9. FACT is both necessary and sufficient to displace dCas9, and FACT immunodepletion converts Cas9's activity from multi-turnover to single turnover. In human cells, FACT depletion extends dCas9 residence times, delays genome editing, and alters the balance between indel formation and homology-directed repair. FACT knockdown also increases epigenetic marking by dCas9-based transcriptional effectors with a concomitant enhancement of transcriptional modulation. FACT thus shapes the intrinsic cellular response to Cas9-based genome manipulation most likely by determining Cas9 residence times.

Timed inhibition of CDC7 increases CRISPR-Cas9 mediated templated repair.

Repair of double strand DNA breaks (DSBs) can result in gene disruption or gene modification via homology directed repair (HDR) from donor DNA. Altering cellular responses to DSBs may rebalance editing outcomes towards HDR and away from other repair outcomes. Here, we utilize a pooled CRISPR screen to define host cell involvement in HDR between a Cas9 DSB and a plasmid double stranded donor DNA (dsDonor). We find that the Fanconi Anemia (FA) pathway is required for dsDonor HDR and that other genes act to repress HDR. Small molecule inhibition of one of these repressors, CDC7, by XL413 and other inhibitors increases the efficiency of HDR by up to 3.5 fold in many contexts, including primary T cells. XL413 stimulates HDR during a reversible slowing of S-phase that is unexplored for Cas9-induced HDR. We anticipate that XL413 and other such rationally developed inhibitors will be useful tools for gene modification.

CRISPR off-target detection with DISCOVER-seq.

DISCOVER-seq (discovery of in situ Cas off-targets and verification by sequencing) is a broadly applicable approach for unbiased CRISPR-Cas off-target identification in cells and tissues. It leverages the recruitment of DNA repair factors to double-strand breaks (DSBs) after genome editing with CRISPR nucleases. Here, we describe a detailed experimental protocol and analysis pipeline with which to perform DISCOVER-seq. The principle of this method is to track the precise recruitment of MRE11 to DSBs by chromatin immunoprecipitation followed by next-generation sequencing. A customized open-source bioinformatics pipeline, BLENDER (blunt end finder), then identifies off-target sequences genome wide. DISCOVER-seq is capable of finding and measuring off-targets in primary cells and in situ. The two main advantages of DISCOVER-seq are (i) low false-positive rates because DNA repair enzyme binding is required for genome edits to occur and (ii) its applicability to a wide variety of systems, including patient-derived cells and animal models. The whole protocol, including the analysis, can be completed within 2 weeks.

Author Correction: SLC19A1 transports immunoreactive cyclic dinucleotides.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

SLC19A1 transports immunoreactive cyclic dinucleotides.

The accumulation of DNA in the cytosol serves as a key immunostimulatory signal associated with infections, cancer and genomic damage1,2. Cytosolic DNA triggers immune responses by activating the cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) pathway3. The binding of DNA to cGAS activates its enzymatic activity, leading to the synthesis of a second messenger, cyclic guanosine monophosphate-adenosine monophosphate (2'3'-cGAMP)4-7. This cyclic dinucleotide (CDN) activates STING8, which in turn activates the transcription factors interferon regulatory factor 3 (IRF3) and nuclear factor κ-light-chain-enhancer of activated B cells (NF-κB), promoting the transcription of genes encoding type I interferons and other cytokines and mediators that stimulate a broader immune response. Exogenous 2'3'-cGAMP produced by malignant cells9 and other CDNs, including those produced by bacteria10-12 and synthetic CDNs used in cancer immunotherapy13,14, must traverse the cell membrane to activate STING in target cells. How these charged CDNs pass through the lipid bilayer is unknown. Here we used a genome-wide CRISPR-interference screen to identify the reduced folate carrier SLC19A1, a folate-organic phosphate antiporter, as the major transporter of CDNs. Depleting SLC19A1 in human cells inhibits CDN uptake and functional responses, and overexpressing SLC19A1 increases both uptake and functional responses. In human cell lines and primary cells ex vivo, CDN uptake is inhibited by folates as well as two medications approved for treatment of inflammatory diseases, sulfasalazine and the antifolate methotrexate. The identification of SLC19A1 as the major transporter of CDNs into cells has implications for the immunotherapeutic treatment of cancer13, host responsiveness to CDN-producing pathogenic microorganisms11 and-potentially-for some inflammatory diseases.

Oligogenic inheritance of a human heart disease involving a genetic modifier.

Complex genetic mechanisms are thought to underlie many human diseases, yet experimental proof of this model has been elusive. Here, we show that a human cardiac anomaly can be caused by a combination of rare, inherited heterozygous mutations. Whole-exome sequencing of a nuclear family revealed that three offspring with childhood-onset cardiomyopathy had inherited three missense single-nucleotide variants in the MKL2, MYH7, and NKX2-5 genes. The MYH7 and MKL2 variants were inherited from the affected, asymptomatic father and the rare NKX2-5 variant (minor allele frequency, 0.0012) from the unaffected mother. We used CRISPR-Cas9 to generate mice encoding the orthologous variants and found that compound heterozygosity for all three variants recapitulated the human disease phenotype. Analysis of murine hearts and human induced pluripotent stem cell-derived cardiomyocytes provided histologic and molecular evidence for the NKX2-5 variant's contribution as a genetic modifier.

Unbiased detection of CRISPR off-targets in vivo using DISCOVER-Seq.

CRISPR-Cas genome editing induces targeted DNA damage but can also affect off-target sites. Current off-target discovery methods work using purified DNA or specific cellular models but are incapable of direct detection in vivo. We developed DISCOVER-Seq (discovery of in situ Cas off-targets and verification by sequencing), a universally applicable approach for unbiased off-target identification that leverages the recruitment of DNA repair factors in cells and organisms. Tracking the precise recruitment of MRE11 uncovers the molecular nature of Cas activity in cells with single-base resolution. DISCOVER-Seq works with multiple guide RNA formats and types of Cas enzymes, allowing characterization of new editing tools. Off-targets can be identified in cell lines and patient-derived induced pluripotent stem cells and during adenoviral editing of mice, paving the way for in situ off-target discovery within individual patient genotypes during therapeutic genome editing.

A most wanted list of conserved microbial protein families with no known domains.

The number and proportion of genes with no known function are growing rapidly. To quantify this phenomenon and provide criteria for prioritizing genes for functional characterization, we developed a bioinformatics pipeline that identifies robustly defined protein families with no annotated domains, ranks these with respect to phylogenetic breadth, and identifies them in metagenomics data. We applied this approach to 271 965 protein families from the SFams database and discovered many with no functional annotation, including >118 000 families lacking any known protein domain. From these, we prioritized 6 668 conserved protein families with at least three sequences from organisms in at least two distinct classes. These Function Unknown Families (FUnkFams) are present in Tara Oceans Expedition and Human Microbiome Project metagenomes, with distributions associated with sampling environment. Our findings highlight the extent of functional novelty in sequence databases and establish an approach for creating a "most wanted" list of genes to prioritize for further characterization.

Automated and Accurate Estimation of Gene Family Abundance from Shotgun Metagenomes.

Shotgun metagenomic DNA sequencing is a widely applicable tool for characterizing the functions that are encoded by microbial communities. Several bioinformatic tools can be used to functionally annotate metagenomes, allowing researchers to draw inferences about the functional potential of the community and to identify putative functional biomarkers. However, little is known about how decisions made during annotation affect the reliability of the results. Here, we use statistical simulations to rigorously assess how to optimize annotation accuracy and speed, given parameters of the input data like read length and library size. We identify best practices in metagenome annotation and use them to guide the development of the Shotgun Metagenome Annotation Pipeline (ShotMAP). ShotMAP is an analytically flexible, end-to-end annotation pipeline that can be implemented either on a local computer or a cloud compute cluster. We use ShotMAP to assess how different annotation databases impact the interpretation of how marine metagenome and metatranscriptome functional capacity changes across seasons. We also apply ShotMAP to data obtained from a clinical microbiome investigation of inflammatory bowel disease. This analysis finds that gut microbiota collected from Crohn's disease patients are functionally distinct from gut microbiota collected from either ulcerative colitis patients or healthy controls, with differential abundance of metabolic pathways related to host-microbiome interactions that may serve as putative biomarkers of disease.

Development of DNA confirmatory and high-risk diagnostic testing for newborns using targeted next-generation DNA sequencing.

PURPOSE: Genetic testing is routinely used for second-tier confirmation of newborn sequencing results to rule out false positives and to confirm diagnoses in newborns undergoing inpatient and outpatient care. We developed a targeted next-generation sequencing panel coupled with a variant processing pipeline and demonstrated utility and performance benchmarks across multiple newborn disease presentations in a retrospective clinical study. METHODS: The test utilizes an in silico gene filter that focuses directly on 126 genes related to newborn screening diseases and is applied to the exome or a next-generation sequencing panel called NBDx. NBDx targets the 126 genes and additional newborn-specific disorders. It integrates DNA isolation from minimally invasive biological specimens, targeted next-generation screening, and rapid characterization of genetic variation. RESULTS: We report a rapid parallel processing of 8 to 20 cases within 105 hours with high coverage on our NBDx panel. Analytical sensitivity of 99.8% was observed across known mutation hotspots. Concordance calls with or without clinical summaries were 94% and 75%, respectively. CONCLUSION: Rapid, automated targeted next-generation sequencing and analysis are practical in newborns for second-tier confirmation and neonatal intensive care unit diagnoses, laying a foundation for future primary DNA-based molecular screening of additional disorders and improving existing molecular testing options for newborns.