Discovery of target genes and pathways at GWAS loci by pooled single-cell CRISPR screens.

Most variants associated with complex traits and diseases identified by genome-wide association studies (GWAS) map to noncoding regions of the genome with unknown effects. Using ancestrally diverse, biobank-scale GWAS data, massively parallel CRISPR screens, and single-cell transcriptomic and proteomic sequencing, we discovered 124 cis-target genes of 91 noncoding blood trait GWAS loci. Using precise variant insertion through base editing, we connected specific variants with gene expression changes. We also identified trans-effect networks of noncoding loci when cis target genes encoded transcription factors or microRNAs. Networks were themselves enriched for GWAS variants and demonstrated polygenic contributions to complex traits. This platform enables massively parallel characterization of the target genes and mechanisms of human noncoding variants in both cis and trans.

Expanding the Scope of Bacterial CRISPR Activation with PAM-Flexible dCas9 Variants.

CRISPR-Cas transcriptional tools have been widely applied for programmable regulation of complex biological networks. In comparison to eukaryotic systems, bacterial CRISPR activation (CRISPRa) has stringent target site requirements for effective gene activation. While genes may not always have an NGG protospacer adjacent motif (PAM) at the appropriate position, PAM-flexible dCas9 variants can expand the range of targetable sites. Here we systematically evaluate a panel of PAM-flexible dCas9 variants for their ability to activate bacterial genes. We observe that dxCas9-NG provides a high dynamic range of gene activation for sites with NGN PAMs while dSpRY permits modest activity across almost any PAM. Similar trends were observed for heterologous and endogenous promoters. For all variants tested, improved PAM-flexibility comes with the trade-off that CRISPRi-mediated gene repression becomes less effective. Weaker CRISPR interference (CRISPRi) gene repression can be partially rescued by expressing multiple sgRNAs to target many sites in the gene of interest. Our work provides a framework to choose the most effective dCas9 variant for a given set of gene targets, which will further expand the utility of CRISPRa/i gene regulation in bacterial systems.

A genome-scale screen for synthetic drivers of T cell proliferation.

The engineering of autologous patient T cells for adoptive cell therapies has revolutionized the treatment of several types of cancer1. However, further improvements are needed to increase response and cure rates. CRISPR-based loss-of-function screens have been limited to negative regulators of T cell functions2-4 and raise safety concerns owing to the permanent modification of the genome. Here we identify positive regulators of T cell functions through overexpression of around 12,000 barcoded human open reading frames (ORFs). The top-ranked genes increased the proliferation and activation of primary human CD4+ and CD8+ T cells and their secretion of key cytokines such as interleukin-2 and interferon-γ. In addition, we developed the single-cell genomics method OverCITE-seq for high-throughput quantification of the transcriptome and surface antigens in ORF-engineered T cells. The top-ranked ORF-lymphotoxin-ß receptor (LTBR)-is typically expressed in myeloid cells but absent in lymphocytes. When overexpressed in T cells, LTBR induced profound transcriptional and epigenomic remodelling, leading to increased T cell effector functions and resistance to exhaustion in chronic stimulation settings through constitutive activation of the canonical NF-κB pathway. LTBR and other highly ranked genes improved the antigen-specific responses of chimeric antigen receptor T cells and γδ T cells, highlighting their potential for future cancer-agnostic therapies5. Our results provide several strategies for improving next-generation T cell therapies by the induction of synthetic cell programmes.

Integrative approach identifies SLC6A20 and CXCR6 as putative causal genes for the COVID-19 GWAS signal in the 3p21.31 locus.

To date, the locus with the most robust human genetic association to COVID-19 severity 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 modulate COVID-19 risk and highlight the usefulness of this integrative approach to bridge the divide between correlational and causal studies of human biology.

An autoimmune disease risk variant: A trans master regulatory effect mediated by IRF1 under immune stimulation?

Functional mechanisms remain unknown for most genetic loci associated to complex human traits and diseases. In this study, we first mapped trans-eQTLs in a data set of primary monocytes stimulated with LPS, and discovered that a risk variant for autoimmune disease, rs17622517 in an intron of C5ORF56, affects the expression of the transcription factor IRF1 20 kb away. The cis-regulatory effect specific to IRF1 is active under early immune stimulus, with a large number of trans-eQTL effects across the genome under late LPS response. Using CRISPRi silencing, we showed that perturbation of the SNP locus downregulates IRF1 and causes widespread transcriptional effects. Genome editing by CRISPR had suggestive recapitulation of the LPS-specific trans-eQTL signal and lent support for the rs17622517 site being functional. Our results suggest that this common genetic variant affects inter-individual response to immune stimuli via regulation of IRF1. For this autoimmune GWAS locus, our work provides evidence of the functional variant, demonstrates a condition-specific enhancer effect, identifies IRF1 as the likely causal gene in cis, and indicates that overactivation of the downstream immune-related pathway may be the cellular mechanism increasing disease risk. This work not only provides rare experimental validation of a master-regulatory trans-eQTL, but also demonstrates the power of eQTL mapping to build mechanistic hypotheses amenable for experimental follow-up using the CRISPR toolkit.

Integrative approach identifies SLC6A20 and CXCR6 as putative causal genes for the COVID-19 GWAS signal in the 3p21.31 locus.

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.

The Spike D614G mutation increases SARS-CoV-2 infection of multiple human cell types.

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 from 1.5- to 7.7-fold in A549ACE2 and Huh7.5ACE2 overexpressing ACE2. Furthermore, trans-complementation of SARS-CoV-2 virus with Spike D614G showed an increased infectivity in human cells. Although there is minimal difference in ACE2 receptor binding between the D614 and G614 Spike variants, the G614 variant is more resistant to proteolytic cleavage, suggesting a possible mechanism for the increased transduction.

Identification of Required Host Factors for SARS-CoV-2 Infection in Human Cells.

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.

Massively parallel Cas13 screens reveal principles for guide RNA design.

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 Spike D614G mutation increases SARS-CoV-2 infection of multiple human cell types.

A novel isolate 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 to introduce the D614G variant and show that in multiple cell lines, including human lung epithelial cells, that the D614G mutation is up to 8-fold more effective at transducing cells than wild-type. We demonstrate increased infection using both Spike-pseudotyped lentivirus and intact SARS-CoV-2 virus. Although there is minimal difference in ACE2 receptor binding between the Spike variants, we show that the G614 variant is more resistant to proteolytic cleavage in vitro and in human cells, suggesting a possible mechanism for the increased transduction. This result has important implications for the efficacy of Spike-based vaccines currently under development in protecting against this recent and highly-prevalent SARS-CoV-2 isolate.

High-Throughput Screens of PAM-Flexible Cas9 Variants for Gene Knockout and Transcriptional Modulation.

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.

Resolution of human ribosomal DNA occurs in anaphase, dependent on tankyrase 1, condensin II, and topoisomerase IIα.

Formation of individualized sister chromatids is essential for their accurate segregation. In budding yeast, while most of the genome segregates at the metaphase to anaphase transition, resolution of the ribosomal DNA (rDNA) repeats is delayed. The timing and mechanism in human cells is unknown. Here we show that resolution of human rDNA occurs in anaphase after the bulk of the genome, dependent on tankyrase 1, condensin II, and topoisomerase IIα. Defective resolution leads to rDNA bridges, rDNA damage, and aneuploidy of an rDNA-containing acrocentric chromosome. Thus, temporal regulation of rDNA segregation is conserved between yeast and man and is essential for genome integrity.

Loss of Tumor Suppressor STAG2 Promotes Telomere Recombination and Extends the Replicative Lifespan of Normal Human Cells.

Sister chromatids are held together by cohesin, a tripartite ring with a peripheral SA1/2 subunit, where SA1 is required for telomere cohesion and SA2 for centromere cohesion. The STAG2 gene encoding SA2 is often inactivated in human cancer, but not in in a manner associated with aneuploidy. Thus, how these tumors maintain chromosomal cohesion and how STAG2 loss contributes to tumorigenesis remain open questions. Here we show that, despite a loss in centromere cohesion, sister chromatids in STAG2 mutant tumor cells maintain cohesion in mitosis at chromosome arms and telomeres. Telomere maintenance in STAG2 mutant tumor cells occurred by either telomere recombination or telomerase activation mechanisms. Notably, these cells were refractory to telomerase inhibitors, indicating recombination can provide an alternative means of telomere maintenance. STAG2 silencing in normal human cells that lack telomerase led to increased recombination at telomeres, delayed telomere shortening, and postponed senescence onset. Insofar as telomere shortening and replicative senescence prevent genomic instability and cancer by limiting the number of cell divisions, our findings suggest that extending the lifespan of normal human cells due to inactivation of STAG2 could promote tumorigenesis by extending the period during which tumor-driving mutations occur. Cancer Res; 77(20); 5530-42. ©2017 AACR.

SA1 binds directly to DNA through its unique AT-hook to promote sister chromatid cohesion at telomeres.

Sister chromatid cohesion relies on cohesin, a complex comprising a tri-partite ring and a peripheral subunit Scc3, which is found as two related isoforms SA1 and SA2 in vertebrates. There is a division of labor between the vertebrate cohesin complexes; SA1-cohesin is required at telomeres and SA2-cohesin at centromeres. Depletion of SA1 has dramatic consequences for telomere function and genome integrity, but the mechanism by which SA1-cohesin mediates cohesion at telomeres is not well understood. Here we dissect the individual contribution of SA1 and the ring subunits to telomere cohesion and show that telomeres rely heavily on SA1 and to a lesser extent on the ring for cohesion. Using chromatin immunoprecipitation we show that SA1 is highly enriched at telomeres, is decreased at mitosis when cohesion is resolved, and is increased when cohesion persists. Overexpression of SA1 alone was sufficient to induce cohesion at telomeres, independent of the cohesin ring and dependent on its unique (not found in SA2) N-terminal domain, which we show binds to telomeric DNA through an AT-hook motif. We suggest that a specialized cohesion mechanism may be required to accommodate the high level of DNA replication-associated repair at telomeres.