Faculty Bibliography

To date the locus with the most robust human genetic association to COVID-19 susceptibility is 3p21.31. Here, we integrate genome-scale CRISPR loss-of-function screens and eQTLs in diverse cell types and tissues to pinpoint genes underlying COVID-19 risk. Our findings identify SLC6A20 and CXCR6 as putative causal genes that mediate COVID-19 risk and highlight the usefulness of this integrative approach to bridge the divide between correlational and causal studies of human biology.

A novel variant of the SARS-CoV-2 virus carrying a point mutation in the Spike protein (D614G) has recently emerged and rapidly surpassed others in prevalence. This mutation is in linkage disequilibrium with an ORF1b protein variant (P314L), making it difficult to discern the functional significance of the Spike D614G mutation from population genetics alone. Here, we perform site-directed mutagenesis on wild-type human codon optimized Spike to introduce the D614G variant. Using multiple human cell lines, including human lung epithelial cells, we found that the lentiviral particles pseudotyped with Spike D614G are more effective at transducing cells than ones pseudotyped with wild-type Spike. The increased transduction with Spike D614G ranged from 1.3 to 2.4-fold in Caco-2 and Calu-3 cells expressing endogenous ACE2, and 1.5 to 7.7-fold in A549^ACE2 and Huh7.5^ACE2 overexpressing ACE2. Furthermore, trans-complementation of SARS-CoV-2 virus with Spike D614G showed an increased infectivity of human cells. Although there is minimal difference in ACE2 receptor binding between the D614 and G614 Spike variants, we show that the G614 variant is more resistant to proteolytic cleavage in human cells, suggesting a possible mechanism for the increased transduction.

To better understand host-virus genetic dependencies and find potential therapeutic targets for COVID-19, we performed a genome-scale CRISPR loss-of-function screen to identify host factors required for SARS-CoV-2 viral infection of human alveolar epithelial cells. Top-ranked genes cluster into distinct pathways, including the vacuolar ATPase proton pump, Retromer, and Commander complexes. We validate these gene targets using several orthogonal methods such as CRISPR knockout, RNA interference knockdown, and small-molecule inhibitors. Using single-cell RNA-sequencing, we identify shared transcriptional changes in cholesterol biosynthesis upon loss of top-ranked genes. In addition, given the key role of the ACE2 receptor in the early stages of viral entry, we show that loss of RAB7A reduces viral entry by sequestering the ACE2 receptor inside cells. Overall, this work provides a genome-scale, quantitative resource of the impact of the loss of each host gene on fitness/response to viral infection.

Recently, a novel isolate of the SARS-CoV-2 virus carrying a point mutation in the Spike protein (D614G) has emerged and rapidly surpassed others in prevalence, including the original SARS-CoV-2 isolate from Wuhan, China. This Spike variant is a defining feature of the most prevalent clade (A2a) of SARS-CoV-2 genomes worldwide. Using phylogenomic data, several groups have proposed that the D614G variant may confer increased transmissibility leading to positive selection, while others have claimed that currently available evidence does not support positive selection. Furthermore, in the A2a clade, this mutation is in linkage disequilibrium with a ORF1b protein variant (P314L), making it difficult to discern the functional significance of the Spike D614G mutation from population genetics alone. Here, we perform site-directed mutagenesis on a human codon-optimized spike protein to introduce the D614G variant and produce SARS-CoV-2-pseudotyped lentiviral particles (S-Virus) with this variant and with D614 Spike. We show that in multiple cell lines, including human lung epithelial cells, that S-Virus carrying the D614G mutation is up to 8-fold more effective at transducing cells than wild-type S-Virus. This provides functional evidence that the D614G mutation in the Spike protein increases transduction of human cells. Further we show that the G614 variant is more resistant to cleavage in vitro and in human cells, which may suggest a possible mechanism for the increased transduction. Given that several vaccines in development and in clinical trials are based on the initial (D614) Spike sequence, this result has important implications for the efficacy of these vaccines in protecting against this recent and highly-prevalent SARS-CoV-2 isolate.

Type VI CRISPR enzymes are RNA-targeting proteins with nuclease activity that enable specific and robust target gene knockdown without altering the genome. To define rules for the design of Cas13d guide RNAs (gRNAs), we conducted massively parallel screens targeting messenger RNAs (mRNAs) of a green fluorescent protein transgene, and CD46, CD55 and CD71 cell-surface proteins in human cells. In total, we measured the activity of 24,460 gRNAs with and without mismatches relative to the target sequences. Knockdown efficacy is driven by gRNA-specific features and target site context. Single mismatches generally reduce knockdown to a modest degree, but spacer nucleotides 15-21 are largely intolerant of target site mismatches. We developed a computational model to identify optimal gRNAs and confirm their generalizability, testing 3,979 guides targeting mRNAs of 48 endogenous genes. We show that Cas13 can be used in forward transcriptomic pooled screens and, using our model, predict optimized Cas13 gRNAs for all protein-coding transcripts in the human genome.

The emergence of CRISPR-Cas9 gene-editing technologies and genome-wide CRISPR-Cas9 libraries enables efficient unbiased genetic screening that can accelerate the process of therapeutic discovery for genetic disorders. Here, we demonstrate the utility of a genome-wide CRISPR-Cas9 loss-of-function library to identify therapeutic targets for facioscapulohumeral muscular dystrophy (FSHD), a genetically complex type of muscular dystrophy for which there is currently no treatment. In FSHD, both genetic and epigenetic changes lead to misexpression of DUX4, the FSHD causal gene that encodes the highly cytotoxic DUX4 protein. We performed a genome-wide CRISPR-Cas9 screen to identify genes whose loss-of-function conferred survival when DUX4 was expressed in muscle cells. Genes emerging from our screen illuminated a pathogenic link to the cellular hypoxia response, which was revealed to be the main driver of DUX4-induced cell death. Application of hypoxia signaling inhibitors resulted in increased DUX4 protein turnover and subsequent reduction of the cellular hypoxia response and cell death. In addition, these compounds proved successful in reducing FSHD disease biomarkers in patient myogenic lines, as well as improving structural and functional properties in two zebrafish models of FSHD. Our genome-wide perturbation of pathways affecting DUX4 expression has provided insight into key drivers of DUX4-induced pathogenesis and has identified existing compounds with potential therapeutic benefit for FSHD. Our experimental approach presents an accelerated paradigm toward mechanistic understanding and therapeutic discovery of a complex genetic disease, which may be translatable to other diseases with well-established phenotypic selection assays.

A key limitation of the widely used CRISPR enzyme S. pyogenes Cas9 is the strict requirement of an NGG protospacer-adjacent motif (PAM) at the target site. This constraint can be limiting for genome editing applications that require precise Cas9 positioning. Recently, two Cas9 variants with a relaxed PAM requirement (NG) have been developed (xCas9 and Cas9-NG), but their activity has been measured at only a small number of endogenous sites. Here, we devise a high-throughput Cas9 pooled competition screen to compare the performance of Cas9 variants at thousands of genomic loci for gene knockout, transcriptional activation, and inhibition. We show that PAM flexibility comes at a substantial cost of decreased DNA targeting and cleavage. Of the PAM-flexible variants, we find that Cas9-NG outperforms xCas9 regardless of genome engineering modality or PAM. Finally, we combine xCas9 mutations with those of Cas9-NG, creating a stronger transcriptional modulator than existing PAM-flexible Cas9 variants.

Multimodal single-cell assays provide high-resolution snapshots of complex cell populations, but are mostly limited to transcriptome plus an additional modality. Here, we describe expanded CRISPR-compatible cellular indexing of transcriptomes and epitopes by sequencing (ECCITE-seq) for the high-throughput characterization of at least five modalities of information from each single cell. We demonstrate application of ECCITE-seq to multimodal CRISPR screens with robust direct single-guide RNA capture and to clonotype-aware multimodal phenotyping of cancer samples.

Genome-scale forward genetic screens elucidate the genetic basis of therapeutic resistance, tumor evolution, and metastasis in diverse human cancers.

Alexander disease (AxD) is a leukodystrophy that primarily affects astrocytes and is caused by mutations in the astrocytic filament gene GFAP. While astrocytes are thought to have important roles in controlling myelination, AxD animal models do not recapitulate critical myelination phenotypes and it is therefore not clear how AxD astrocytes contribute to leukodystrophy. Here, we show that AxD patient iPSC-derived astrocytes recapitulate key features of AxD pathology such as GFAP aggregation. Moreover, AxD astrocytes inhibit proliferation of human iPSC-derived oligodendrocyte progenitor cells (OPCs) in co-culture and reduce their myelination potential. CRISPR/Cas9-based correction of GFAP mutations reversed these phenotypes. Transcriptomic analyses of AxD astrocytes and postmortem brains identified CHI3L1 as a key mediator of AxD astrocyte-induced inhibition of OPC activity. Thus, this iPSC-based model of AxD not only recapitulates patient phenotypes not observed in animal models, but also reveals mechanisms underlying disease pathology and provides a platform for assessing therapeutic interventions.