Search

Multiplexed assays of variant effect are powerful tools for assessing the impact of protein sequence variation, but are limited to measuring a single protein property and often rely on indirect readouts of intracellular protein function. Here, we developed LAbeling with Barcodes and Enrichment for biochemicaL analysis by sequencing (LABEL-seq), a platform for the multimodal profiling of thousands of protein variants in cultured human cells. Multimodal measurement of ~20,000 variant effects for ~1,600 BRaf variants using LABEL-seq revealed that variation at positions that are frequently mutated in cancer had minimal effects on folding and intracellular abundance but could dramatically alter activity, protein-protein interactions, and druggability. Integrative analysis of our multimodal measurements identified networks of positions with similar roles in regulating BRaf's signaling properties and enabled predictive modeling of variant effects on complex processes such as cell proliferation and small molecule-promoted degradation. LABEL-seq provides a scalable approach for the direct measurement of multiple biochemical effects of protein variants in their native cellular context, yielding insight into protein function, disease mechanisms, and druggability.

Long read sequencing technologies, an attractive solution for many applications, often suffer from higher error rates. Alignment of multiple reads can improve base-calling accuracy, but some applications, e.g. sequencing mutagenized libraries where multiple distinct clones differ by one or few variants, require the use of barcodes or unique molecular identifiers. Unfortunately, sequencing errors can interfere with correct barcode identification, and a given barcode sequence may be linked to multiple independent clones within a given library. Here we focus on the target application of sequencing mutagenized libraries in the context of multiplexed assays of variant effects (MAVEs). MAVEs are increasingly used to create comprehensive genotype-phenotype maps that can aid clinical variant interpretation. Many MAVE methods use long-read sequencing of barcoded mutant libraries for accurate association of barcode with genotype. Existing long-read sequencing pipelines do not account for inaccurate sequencing or non-unique barcodes. Here, we describe Pacybara, which handles these issues by clustering long reads based on the similarities of (error-prone) barcodes while also detecting barcodes that have been associated with multiple genotypes. Pacybara also detects recombinant (chimeric) clones and reduces false positive indel calls. In three example applications, we show that Pacybara identifies and correctly resolves these issues.

Background: KCNE1 encodes a 129-residue cardiac potassium channel (I Ks ) subunit. KCNE1 variants are associated with long QT syndrome and atrial fibrillation. However, most variants have insufficient evidence of clinical consequences and thus limited clinical utility., Results: Here, we demonstrate the power of variant effect mapping, which couples saturation mutagenesis with high-throughput sequencing, to ascertain the function of thousands of protein coding KCNE1 variants. We comprehensively assayed KCNE1 variant cell surface expression (2,554/2,709 possible single amino acid variants) and function (2,539 variants). We identified 470 loss-of-surface expression and 588 loss-of-function variants. Out of the 588 loss-of-function variants, only 155 had low cell surface expression. The latter half of the protein is dispensable for protein trafficking but essential for channel function. 22 of the 30 KCNE1 residues (73%) highly intolerant of variation were in predicted close contact with binding partners KCNQ1 or calmodulin. Our data were highly concordant with gold standard electrophysiological data (ρ = -0.65), population and patient cohorts (32/38 concordant variants), and computational metrics (ρ = -0.55). Our data provide moderate-strength evidence for the ACMG/AMP functional criteria for benign and pathogenic variants., Conclusions: Comprehensive variant effect maps of KCNE1 can both provide insight into I Ks channel biology and help reclassify variants of uncertain significance., Competing Interests: Competing interests The authors declare that they have no competing interests.

Standard transgenic cell line generation requires screening 100-1000s of colonies to isolate correctly edited cells. We describe CR ISPR a On- T arget E diting R etrieval (CRaTER) which enriches for cells with on-target knock-in of a cDNA-fluorescent reporter transgene by transient activation of the targeted locus followed by flow sorting to recover edited cells. We show CRaTER recovers rare cells with heterozygous, biallelic-editing of the transcriptionally-inactive MYH7 locus in human induced pluripotent stem cells (hiPSCs), enriching on average 25-fold compared to standard antibiotic selection. We leveraged CRaTER to enrich for heterozygous knock-in of a library of single nucleotide variants (SNVs) in MYH7 , a gene in which missense mutations cause cardiomyopathies, and recovered hiPSCs with 113 different MYH7 SNVs. We differentiated these hiPSCs to cardiomyocytes and show MYH7 fusion proteins can localize as expected. Thus, CRaTER substantially reduces screening required for isolation of gene-edited cells, enabling generation of transgenic cell lines at unprecedented scale.