Abstract:Discovery of molecular mechanisms responsible for trait associations as discovered by genome-wide association studies (GWAS) is hampered by difficulty in identifying causal genetic variants due to linkage disequilibrium. Typical assays of genetic function are low throughput or evaluate sequences in heterologous ectopic settings. Genome editing enables perturbation of trait-associated genetic sequences within relevant genomic, chromatin and cellular context. Here we perform comprehensive analysis of genetic variants associated with red blood cell traits by pooled CRISPR screening. We performed a genome-wide Cas9 gene knockout screen in immortalized erythroid precursors (HUDEP-2 cells) during erythroid maturation to define functional erythroid genes required for cell growth or differentiation. We evaluated 952 loci associated with nine red blood cell traits (Astle et al, Cell 2016) comprising 24,843 SNPs. We linked 7,187 (28.9%) of these SNPs to genes by at least one of four routes: sharing topological associated domain, physical proximity (<20 kb), long-range chromatin interaction (promoter HiC), or eQTL with a functional erythroid gene. We designed ~5 guide RNAs per SNP requiring cleavage position within at least 50 bp and exceeding an off-target score threshold, resulting in 32,710 sgRNAs testing 5,592 SNPs at 481 loci. We utilized four editors: Cas9 nuclease to produce indels, dCas9-VP64 for gene activation, dCas9-KRAB for gene repression, and dCas9 as a DNA targeting control. By pooled lentiviral transduction, erythroid differentiation culture, and guide RNA library deep sequencing, we found reproducible results across biological replicates, with guide count clustered by Cas9 protein type. We performed fine-mapping of association results by Bayesian inference to calculate posterior probability of inclusion (PPI). We found a strong correlation between PPI and CRISPR significance score indicating agreement between CRISPR screening and genetic fine-mapping, despite examples of validated CRISPR signals with low PPI scores. We identified numerous CRISPR-implicated functional SNPs at regulatory elements including promoters and enhancers but also at noncoding sequences lacking chromatin marks. The editing of CRISPR -implicated functional SNPs correlated well with expression of linked genes. Editing of CRISPR-implicated functional SNPs caused altered proliferation and/or differentiation of both HUDEP-2 cells and CD34+ HSPC-derived primary erythroid precursors. We validated a functional SNP at the BCL2L1 enhancer, where both Cas9 disruption and dCas9-KRAB inhibition resulted in decreased expression of BCL2L1,and reduced cell survival during erythropoiesis. We validated several CRISPR-implicated functional SNPs at GFI1B that controlled erythroid cell differentiation, including one at a distal enhancer element and another that impacted exon splicing. These results demonstrate the potential and challenges of comprehensive integrated genomic perturbation to complement genetic, biochemical, and statistical approaches to uncover molecular underpinnings of human blood cell traits.

Variant-to-function dissection of rare non-coding GWAS loci with high impact on blood traits

From GWAS Variant to Function: A Study of ∼148,000 Variants for Blood Cell Traits

Systematic Functional Dissection of Common Genetic Variation Affecting Red Blood Cell Traits

Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations

Discovery of target genes and pathways of blood trait loci using pooled CRISPR screens and single cell RNA sequencing

Transcriptome-wide Association Study in UK Biobank Europeans Identifies Associations with Blood Cell Traits

Whole-genome Sequencing Association Analysis of Quantitative Red Blood Cell Phenotypes: the NHLBI TOPMed Program.

The Polygenic and Monogenic Basis of Blood Traits and Diseases

Comprehensive Integrated Genomic Perturbations Reveal Molecular Mechanisms of Red Blood Cell Trait Associations

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

A large-scale transcriptome-wide association study (TWAS) of ten blood cell phenotypes reveals complexities of TWAS fine-mapping

The functional impact of rare variation across the regulatory cascade

Principles for the Post-Gwas Functional Characterisation of Risk Loci

Mapping structural variants to rare disease genes using long-read whole genome sequencing and trait-relevant polygenic scores

Genetic Analyses of Blood Cell Structure for Biological and Pharmacological Inference

DeepWAS: Multivariate Genotype-Phenotype Associations by Directly Integrating Regulatory Information Using Deep Learning

The Allelic Landscape of Human Blood Cell Trait Variation and Links to Common Complex Disease

Whole genome sequencing identifies structural variants contributing to hematologic traits in the NHLBI TOPMed program

Whole-genome Sequencing in Diverse Subjects Identifies Genetic Correlates of Leukocyte Traits: the NHLBI TOPMed Program

Genome-wide association testing beyond SNPs

Unraveling the polygenic architecture of complex traits using blood eQTL metaanalysis