Multiplexed assays of variant effects contribute to a growing genotype-phenotype atlas

Saturation genome editing-based clinical classification of BRCA2 variants

Sequencing-based genetic tests have uncovered a vast array of BRCA2 sequence variants1. Owing to limited clinical, familial and epidemiological data, thousands of variants are considered to be variants of uncertain significance2-4 (VUS). Here we have utilized CRISPR-Cas9-based saturation genome editing in a humanized mouse embryonic stem cell line to determine the functional effect of VUS. We have categorized nearly all possible single nucleotide variants (SNVs) in the region that encodes the carboxylate-terminal DNA-binding domain of BRCA2. We have generated function scores for 6,551 SNVs, covering 96.4% of possible SNVs in exons 15-26 spanning BRCA2 residues 2479-3216. These variants include 1,282 SNVs that are categorized as missense VUS in the clinical variant database ClinVar, with 77.2% of these classified as benign and 20.4% classified as pathogenic using our functional score. Our assay provides evidence that 3,384 of the SNVs in the region are benign and 776 are pathogenic. Our classification aligns closely with pathogenicity data from ClinVar, orthogonal functional assays and computational meta predictors. We have integrated our embryonic stem cell-based BRCA2-saturation genome editing dataset with other available evidence and utilized the American College of Medical Genetics and Genomics/Association for Molecular Pathology guidelines for clinical classification of all possible SNVs. This classification is available as a sequence-function map and serves as a valuable resource for interpreting unidentified variants in the population and for physicians and genetic counsellors to assess BRCA2 VUS in patients.

Functional evaluation and clinical classification of BRCA2 variants

Germline BRCA2 loss-of function variants, which can be identified through clinical genetic testing, predispose to several cancers1-5. However, variants of uncertain significance limit the clinical utility of test results. Thus, there is a need for functional characterization and clinical classification of all BRCA2 variants to facilitate the clinical management of individuals with these variants. Here we analysed all possible single-nucleotide variants from exons 15 to 26 that encode the BRCA2 DNA-binding domain hotspot for pathogenic missense variants. To enable this, we used saturation genome editing CRISPR-Cas9-based knock-in endogenous targeting of human haploid HAP1 cells6. The assay was calibrated relative to nonsense and silent variants and was validated using pathogenic and benign standards from ClinVar and results from a homology-directed repair functional assay7. Variants (6,959 out of 6,960 evaluated) were assigned to seven categories of pathogenicity based on a VarCall Bayesian model8. Single-nucleotide variants that encode loss-of-function missense variants were associated with increased risks of breast cancer and ovarian cancer. The functional assay results were integrated into models from ClinGen, the American College of Medical Genetics and Genomics, and the Association for Molecular Pathology9 for clinical classification of BRCA2 variants. Using this approach, 91% were classified as pathogenic or likely pathogenic or as benign or likely benign. These classified variants can be used to improve clinical management of individuals with a BRCA2 variant.

A missense variant effect map for the human tumor-suppressor protein CHK2

The tumor suppressor CHEK2 encodes the serine/threonine protein kinase CHK2 which, upon DNA damage, is important for pausing the cell cycle, initiating DNA repair, and inducing apoptosis. CHK2 phosphorylation of the tumor suppressor BRCA1 is also important for mitotic spindle assembly and chromosomal stability. Consistent with its cell-cycle checkpoint role, both germline and somatic variants in CHEK2 have been linked to breast and other cancers. Over 90% of clinical germline CHEK2 missense variants are classified as variants of uncertain significance, complicating diagnosis of CHK2-dependent cancer. We therefore sought to test the functional impact of all possible missense variants in CHK2. Using a scalable multiplexed assay based on the ability of human CHK2 to complement DNA sensitivity of Saccharomyces cerevisiae cells lacking the CHEK2 ortholog, RAD53, we generated a systematic "missense variant effect map" for CHEK2 missense variation. The map reflects known biochemical features of CHK2 while offering new biological insights. It also provides strong evidence toward pathogenicity for some clinical missense variants and supporting evidence toward benignity for others. Overall, this comprehensive missense variant effect map contributes to understanding of both known and yet-to-be-observed CHK2 variants.

Multiplexed Assays of Variant Effect and Automated Patch Clamping Improve KCNH2-LQTS Variant Classification and Cardiac Event Risk Stratification

BACKGROUND

Long QT syndrome is a lethal arrhythmia syndrome, frequently caused by rare loss-of-function variants in the potassium channel encoded by KCNH2. Variant classification is difficult, often because of lack of functional data. Moreover, variant-based risk stratification is also complicated by heterogenous clinical data and incomplete penetrance. Here we sought to test whether variant-specific information, primarily from high-throughput functional assays, could improve both classification and cardiac event risk stratification in a large, harmonized cohort of KCNH2 missense variant heterozygotes.

METHODS

We quantified cell-surface trafficking of 18 796 variants in KCNH2 using a multiplexed assay of variant effect (MAVE). We recorded KCNH2 current density for 533 variants by automated patch clamping. We calibrated the strength of evidence of MAVE data according to ClinGen guidelines. We deeply phenotyped 1458 patients with KCNH2 missense variants, including QTc, cardiac event history, and mortality. We correlated variant functional data and Bayesian long QT syndrome penetrance estimates with cohort phenotypes and assessed hazard ratios for cardiac events.

RESULTS

Variant MAVE trafficking scores and automated patch clamping peak tail currents were highly correlated (Spearman rank-order ρ=0.69; n=433). The MAVE data were found to provide up to pathogenic very strong evidence for severe loss-of-function variants. In the cohort, both functional assays and Bayesian long QT syndrome penetrance estimates were significantly predictive of cardiac events when independently modeled with patient sex and corrected QT interval (QTc); however, MAVE data became nonsignificant when peak tail current and penetrance estimates were also available. The area under the receiver operator characteristic curve for 20-year event outcomes based on patient-specific sex and QTc (area under the curve, 0.80 [0.76-0.83]) was improved with prospectively available penetrance scores conditioned on MAVE (area under the curve, 0.86 [0.83-0.89]) or attainable automated patch clamping peak tail current data (area under the curve, 0.84 [0.81-0.88]).

CONCLUSIONS

High-throughput KCNH2 variant MAVE data meaningfully contribute to variant classification at scale, whereas long QT syndrome penetrance estimates and automated patch clamping peak tail current measurements meaningfully contribute to risk stratification of cardiac events in patients with heterozygous KCNH2 missense variants.

There will always be variants of uncertain significance. Analysis of VUSs

The ACMG/AMP guidelines include five categories of which variants of uncertain significance (VUSs) have received increasing attention. Recently, Fowler and Rehm claimed that all or most VUSs could be reclassified as pathogenic or benign within few years. To test this claim, we collected validated benign, pathogenic, VUS and conflicting variants from ClinVar and LOVD and investigated differences at gene, protein, structure, and variant levels. The gene and protein features included inheritance patterns, actionability, functional categories for housekeeping, essential, complete knockout, lethality and haploinsufficient proteins, Gene Ontology annotations, and protein network properties. Structural properties included the location at secondary structural elements, intrinsically disordered regions, transmembrane regions, repeats, conservation, and accessibility. Gene features were distributions of nucleotides, their groupings, codons, and location to CpG islands. The distributions of amino acids and their groups were investigated. VUSs did not markedly differ from other variants. The only major differences were the accessibility and conservation of pathogenic variants, and reduced ratio of repeat-locating variants in VUSs. Thus, all VUSs cannot be distinguished from other types of variants. They display one form of natural biological heterogeneity. Instead of concentrating on eradicating VUSs, the community would benefit from investigating and understanding factors that contribute to phenotypic heterogeneity.

ReNU syndrome - a newly discovered prevalent neurodevelopmental disorder

Two recent papers have identified genetic variants in the noncoding gene RNU4-2 to cause a frequent neurodevelopmental disorder. This work will have a substantial impact on the rare disease community, leading to thousands of diagnoses worldwide. These studies also highlight the untapped diagnostic potential of noncoding regions.

Decoding biology with massively parallel reporter assays and machine learning

Massively parallel reporter assays (MPRAs) are powerful tools for quantifying the impacts of sequence variation on gene expression. Reading out molecular phenotypes with sequencing enables interrogating the impact of sequence variation beyond genome scale. Machine learning models integrate and codify information learned from MPRAs and enable generalization by predicting sequences outside the training data set. Models can provide a quantitative understanding of cis-regulatory codes controlling gene expression, enable variant stratification, and guide the design of synthetic regulatory elements for applications from synthetic biology to mRNA and gene therapy. This review focuses on cis-regulatory MPRAs, particularly those that interrogate cotranscriptional and post-transcriptional processes: alternative splicing, cleavage and polyadenylation, translation, and mRNA decay.

Leaving no patient behind! Expert recommendation in the use of innovative technologies for diagnosing rare diseases

Genetic diagnosis plays a crucial role in rare diseases, particularly with the increasing availability of emerging and accessible treatments. The International Rare Diseases Research Consortium (IRDiRC) has set its primary goal as: "Ensuring that all patients who present with a suspected rare disease receive a diagnosis within one year if their disorder is documented in the medical literature". Despite significant advances in genomic sequencing technologies, more than half of the patients with suspected Mendelian disorders remain undiagnosed. In response, IRDiRC proposes the establishment of "a globally coordinated diagnostic and research pipeline". To help facilitate this, IRDiRC formed the Task Force on Integrating New Technologies for Rare Disease Diagnosis. This multi-stakeholder Task Force aims to provide an overview of the current state of innovative diagnostic technologies for clinicians and researchers, focusing on the patient's diagnostic journey. Herein, we provide an overview of a broad spectrum of emerging diagnostic technologies involving genomics, epigenomics and multi-omics, functional testing and model systems, data sharing, bioinformatics, and Artificial Intelligence (AI), highlighting their advantages, limitations, and the current state of clinical adaption. We provide expert recommendations outlining the stepwise application of these innovative technologies in the diagnostic pathways while considering global differences in accessibility. The importance of FAIR (Findability, Accessibility, Interoperability, and Reusability) and CARE (Collective benefit, Authority to control, Responsibility, and Ethics) data management is emphasized, along with the need for enhanced and continuing education in medical genomics. We provide a perspective on future technological developments in genome diagnostics and their integration into clinical practice. Lastly, we summarize the challenges related to genomic diversity and accessibility, highlighting the significance of innovative diagnostic technologies, global collaboration, and equitable access to diagnosis and treatment for people living with rare disease.

Assessing the role of evolutionary information for enhancing protein language model embeddings

Embeddings from protein Language Models (pLMs) are replacing evolutionary information from multiple sequence alignments (MSAs) as the most successful input for protein prediction. Is this because embeddings capture evolutionary information? We tested various approaches to explicitly incorporate evolutionary information into embeddings on various protein prediction tasks. While older pLMs (SeqVec, ProtBert) significantly improved through MSAs, the more recent pLM ProtT5 did not benefit. For most tasks, pLM-based outperformed MSA-based methods, and the combination of both even decreased performance for some (intrinsic disorder). We highlight the effectiveness of pLM-based methods and find limited benefits from integrating MSAs.

Machine Learning of Three-Dimensional Protein Structures to Predict the Functional Impacts of Genome Variation

Research in the human genome sciences generates a substantial amount of genetic data for hundreds of thousands of individuals, which concomitantly increases the number of variants of unknown significance (VUS). Bioinformatic analyses can successfully reveal rare variants and variants with clear associations with disease-related phenotypes. These studies have had a significant impact on how clinical genetic screens are interpreted and how patients are stratified for treatment. There are few, if any, computational methods for variants comparable to biological activity predictions. To address this gap, we developed a machine learning method that uses protein three-dimensional structures from AlphaFold to predict how a variant will influence changes to a gene's downstream biological pathways. We trained state-of-the-art machine learning classifiers to predict which protein regions will most likely impact transcriptional activities of two proto-oncogenes, nuclear factor erythroid 2 (NFE2L2)-related factor 2 (NRF2) and c-Myc. We have identified classifiers that attain accuracies higher than 80%, which have allowed us to identify a set of key protein regions that lead to significant perturbations in c-Myc or NRF2 transcriptional pathway activities.

Benchmarking computational variant effect predictors by their ability to infer human traits

BACKGROUND

Computational variant effect predictors offer a scalable and increasingly reliable means of interpreting human genetic variation, but concerns of circularity and bias have limited previous methods for evaluating and comparing predictors. Population-level cohorts of genotyped and phenotyped participants that have not been used in predictor training can facilitate an unbiased benchmarking of available methods. Using a curated set of human gene-trait associations with a reported rare-variant burden association, we evaluate the correlations of 24 computational variant effect predictors with associated human traits in the UK Biobank and All of Us cohorts.

RESULTS

AlphaMissense outperformed all other predictors in inferring human traits based on rare missense variants in UK Biobank and All of Us participants. The overall rankings of computational variant effect predictors in these two cohorts showed a significant positive correlation.

CONCLUSION

We describe a method to assess computational variant effect predictors that sidesteps the limitations of previous evaluations. This approach is generalizable to future predictors and could continue to inform predictor choice for personal and clinical genetics.

Sampling Strategies for Experimentally Mapping Molecular Fitness Landscapes Using High-Throughput Methods

Empirical studies of genotype-phenotype-fitness maps of proteins are fundamental to understanding the evolutionary process, in elucidating the space of possible genotypes accessible through mutations in a landscape of phenotypes and fitness effects. Yet, comprehensively mapping molecular fitness landscapes remains challenging since all possible combinations of amino acid substitutions for even a few protein sites are encoded by an enormous genotype space. High-throughput mapping of genotype space can be achieved using large-scale screening experiments known as multiplexed assays of variant effect (MAVEs). However, to accommodate such multi-mutational studies, the size of MAVEs has grown to the point where a priori determination of sampling requirements is needed. To address this problem, we propose calculations and simulation methods to approximate minimum sampling requirements for multi-mutational MAVEs, which we combine with a new library construction protocol to experimentally validate our approximation approaches. Analysis of our simulated data reveals how sampling trajectories differ between simulations of nucleotide versus amino acid variants and among mutagenesis schemes. For this, we show quantitatively that marginal gains in sampling efficiency demand increasingly greater sampling effort when sampling for nucleotide sequences over their encoded amino acid equivalents. We present a new library construction protocol that efficiently maximizes sequence variation, and demonstrate using ultradeep sequencing that the library encodes virtually all possible combinations of mutations within the experimental design. Insights learned from our analyses together with the methodological advances reported herein are immediately applicable toward pooled experimental screens of arbitrary design, enabling further assay upscaling and expanded testing of genotype space.

High-throughput assays to assess variant effects on disease

Interpreting the wealth of rare genetic variants discovered in population-scale sequencing efforts and deciphering their associations with human health and disease present a critical challenge due to the lack of sufficient clinical case reports. One promising avenue to overcome this problem is deep mutational scanning (DMS), a method of introducing and evaluating large-scale genetic variants in model cell lines. DMS allows unbiased investigation of variants, including those that are not found in clinical reports, thus improving rare disease diagnostics. Currently, the main obstacle limiting the full potential of DMS is the availability of functional assays that are specific to disease mechanisms. Thus, we explore high-throughput functional methodologies suitable to examine broad disease mechanisms. We specifically focus on methods that do not require robotics or automation but instead use well-designed molecular tools to transform biological mechanisms into easily detectable signals, such as cell survival rate, fluorescence or drug resistance. Here, we aim to bridge the gap between disease-relevant assays and their integration into the DMS framework.

High-throughput functional mapping of variants in an arrhythmia gene, KCNE1, reveals novel biology

BACKGROUND

KCNE1 encodes a 129-residue cardiac potassium channel (IKs) 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.

METHODS

In this study, we leveraged 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.

RESULTS

We comprehensively assayed KCNE1 variant cell surface expression (2554/2709 possible single-amino-acid variants) and function (2534 variants). Our study identified 470 loss- or partial loss-of-surface expression and 574 loss- or partial loss-of-function variants. Of the 574 loss- or partial loss-of-function variants, 152 (26.5%) had reduced cell surface expression, indicating that most functionally deleterious variants affect channel gating. Nonsense variants at residues 56-104 generally had WT-like trafficking scores but decreased functional scores, indicating that 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 (with > 70% loss-of-function variants) were in predicted close contact with binding partners KCNQ1 or calmodulin. Our functional assay data were consistent with gold standard electrophysiological data (ρ =  - 0.64), population and patient cohorts (32/38 presumed benign or pathogenic variants with consistent scores), and computational predictors (ρ =  - 0.62). Our data provide moderate-strength evidence for the American College of Medical Genetics/Association of Molecular Pathology 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.

PRKN-linked familial Parkinson's disease: cellular and molecular mechanisms of disease-linked variants

Parkinson's disease (PD) is a common and incurable neurodegenerative disorder that arises from the loss of dopaminergic neurons in the substantia nigra and is mainly characterized by progressive loss of motor function. Monogenic familial PD is associated with highly penetrant variants in specific genes, notably the PRKN gene, where homozygous or compound heterozygous loss-of-function variants predominate. PRKN encodes Parkin, an E3 ubiquitin-protein ligase important for protein ubiquitination and mitophagy of damaged mitochondria. Accordingly, Parkin plays a central role in mitochondrial quality control but is itself also subject to a strict protein quality control system that rapidly eliminates certain disease-linked Parkin variants. Here, we summarize the cellular and molecular functions of Parkin, highlighting the various mechanisms by which PRKN gene variants result in loss-of-function. We emphasize the importance of high-throughput assays and computational tools for the clinical classification of PRKN gene variants and how detailed insights into the pathogenic mechanisms of PRKN gene variants may impact the development of personalized therapeutics.

Deep mutational scanning reveals a correlation between degradation and toxicity of thousands of aspartoacylase variants

Unstable proteins are prone to form non-native interactions with other proteins and thereby may become toxic. To mitigate this, destabilized proteins are targeted by the protein quality control network. Here we present systematic studies of the cytosolic aspartoacylase, ASPA, where variants are linked to Canavan disease, a lethal neurological disorder. We determine the abundance of 6152 of the 6260 ( ~ 98%) possible single amino acid substitutions and nonsense ASPA variants in human cells. Most low abundance variants are degraded through the ubiquitin-proteasome pathway and become toxic upon prolonged expression. The data correlates with predicted changes in thermodynamic stability, evolutionary conservation, and separate disease-linked variants from benign variants. Mapping of degradation signals (degrons) shows that these are often buried and the C-terminal region functions as a degron. The data can be used to interpret Canavan disease variants and provide insight into the relationship between protein stability, degradation and cell fitness.

Analysis of Tumor-Associated AXIN1 Missense Mutations Identifies Variants That Activate β-Catenin Signaling

AXIN1 is a major component of the β-catenin destruction complex and is frequently mutated in various cancer types, particularly liver cancers. Truncating AXIN1 mutations are recognized to encode a defective protein that leads to β-catenin stabilization, but the functional consequences of missense mutations are not well characterized. Here, we first identified the GSK3β, β-catenin, and RGS/APC interaction domains of AXIN1 that are the most critical for proper β-catenin regulation. Analysis of 80 tumor-associated variants in these domains identified 18 that significantly affected β-catenin signaling. Coimmunoprecipitation experiments revealed that most of them lost binding to the binding partner corresponding to the mutated domain. A comprehensive protein structure analysis predicted the consequences of these mutations, which largely overlapped with the observed effects on β-catenin signaling in functional experiments. The structure analysis also predicted that loss-of-function mutations within the RGS/APC interaction domain either directly affected the interface for APC binding or were located within the hydrophobic core and destabilized the entire structure. In addition, truncated AXIN1 length inversely correlated with the β-catenin regulatory function, with longer proteins retaining more functionality. These analyses suggest that all AXIN1-truncating mutations at least partially affect β-catenin regulation, whereas this is only the case for a subset of missense mutations. Consistently, most colorectal and liver cancers carrying missense variants acquire mutations in other β-catenin regulatory genes such as APC and CTNNB1. These results will aid the functional annotation of AXIN1 mutations identified in large-scale sequencing efforts or in individual patients.

SIGNIFICANCE

Characterization of 80 tumor-associated missense variants of AXIN1 reveals a subset of 18 mutations that disrupt its β-catenin regulatory function, whereas the majority are passenger mutations.

FiTMuSiC: leveraging structural and (co)evolutionary data for protein fitness prediction

Systematically predicting the effects of mutations on protein fitness is essential for the understanding of genetic diseases. Indeed, predictions complement experimental efforts in analyzing how variants lead to dysfunctional proteins that in turn can cause diseases. Here we present our new fitness predictor, FiTMuSiC, which leverages structural, evolutionary and coevolutionary information. We show that FiTMuSiC predicts fitness with high accuracy despite the simplicity of its underlying model: it was among the top predictors on the hydroxymethylbilane synthase (HMBS) target of the sixth round of the Critical Assessment of Genome Interpretation challenge (CAGI6) and performs as well as much more complex deep learning models such as AlphaMissense. To further demonstrate FiTMuSiC's robustness, we compared its predictions with in vitro activity data on HMBS, variant fitness data on human glucokinase (GCK), and variant deleteriousness data on HMBS and GCK. These analyses further confirm FiTMuSiC's qualities and accuracy, which compare favorably with those of other predictors. Additionally, FiTMuSiC returns two scores that separately describe the functional and structural effects of the variant, thus providing mechanistic insight into why the variant leads to fitness loss or gain. We also provide an easy-to-use webserver at https://babylone.ulb.ac.be/FiTMuSiC , which is freely available for academic use and does not require any bioinformatics expertise, which simplifies the accessibility of our tool for the entire scientific community.

Detection of autism spectrum disorder-related pathogenic trio variants by a novel structure-based approach

BACKGROUND

Glutamatergic synapse dysfunction is believed to underlie the development of Autism Spectrum Disorder (ASD) and Intellectual Disability (ID) in many individuals. However, identification of genetic markers that contribute to synaptic dysfunction in these individuals is notoriously difficult. Based on genomic analysis, structural modeling, and functional data, we recently established the involvement of the TRIO-RAC1 pathway in ASD and ID. Furthermore, we identified a pathological de novo missense mutation hotspot in TRIO's GEF1 domain. ASD/ID-related missense mutations within this domain compromise glutamatergic synapse function and likely contribute to the development of ASD/ID. The number of ASD/ID cases with mutations identified within TRIO's GEF1 domain is increasing. However, tools for accurately predicting whether such mutations are detrimental to protein function are lacking.

METHODS

Here we deployed advanced protein structural modeling techniques to predict potential de novo pathogenic and benign mutations within TRIO's GEF1 domain. Mutant TRIO-9 constructs were generated and expressed in CA1 pyramidal neurons of organotypic cultured hippocampal slices. AMPA receptor-mediated postsynaptic currents were examined in these neurons using dual whole-cell patch clamp electrophysiology. We also validated these findings using orthogonal co-immunoprecipitation and fluorescence lifetime imaging (FLIM-FRET) experiments to assay TRIO mutant overexpression effects on TRIO-RAC1 binding and on RAC1 activity in HEK293/T cells.

RESULTS

Missense mutations in TRIO's GEF1 domain that were predicted to disrupt TRIO-RAC1 binding or stability were tested experimentally and found to greatly impair TRIO-9's influence on glutamatergic synapse function. In contrast, missense mutations in TRIO's GEF1 domain that were predicted to have minimal effect on TRIO-RAC1 binding or stability did not impair TRIO-9's influence on glutamatergic synapse function in our experimental assays. In orthogonal assays, we find most of the mutations predicted to disrupt binding display loss of function but mutants predicted to disrupt stability do not reflect our results from neuronal electrophysiological data.

LIMITATIONS

We present a method to predict missense mutations in TRIO's GEF1 domain that may compromise TRIO function and test for effects in a limited number of assays. Possible limitations arising from the model systems employed here can be addressed in future studies. Our method does not provide evidence for whether these mutations confer ASD/ID risk or the likelihood that such mutations will result in the development of ASD/ID.

CONCLUSIONS

Here we show that a combination of structure-based computational predictions and experimental validation can be employed to reliably predict whether missense mutations in the human TRIO gene impede TRIO protein function and compromise TRIO's role in glutamatergic synapse regulation. With the growing accessibility of genome sequencing, the use of such tools in the accurate identification of pathological mutations will be instrumental in diagnostics of ASD/ID.

Scaling up Functional Analyses of the G Protein-Coupled Receptor Rhodopsin

Eukaryotic cells use G protein-coupled receptors (GPCRs) to convert external stimuli into internal signals to elicit cellular responses. However, how mutations in GPCR-coding genes affect GPCR activation and downstream signaling pathways remain poorly understood. Approaches such as deep mutational scanning show promise in investigations of GPCRs, but a high-throughput method to measure rhodopsin activation has yet to be achieved. Here, we scale up a fluorescent reporter assay in budding yeast that we engineered to study rhodopsin's light-activated signal transduction. Using this approach, we measured the mutational effects of over 1200 individual human rhodopsin mutants, generated by low-frequency random mutagenesis of the GPCR rhodopsin (RHO) gene. Analysis of the data in the context of rhodopsin's three-dimensional structure reveals that transmembrane helices are generally less tolerant to mutations compared to flanking helices that face the lipid bilayer, which suggest that mutational tolerance is contingent on both the local environment surrounding specific residues and the specific position of these residues in the protein structure. Comparison of functional scores from our screen to clinically identified rhodopsin disease variants found many pathogenic mutants to be loss of function. Lastly, functional scores from our assay were consistent with a complex counterion mechanism involved in ligand-binding and rhodopsin activation. Our results demonstrate that deep mutational scanning is possible for rhodopsin activation and can be an effective method for revealing properties of mutational tolerance that may be generalizable to other transmembrane proteins.

Genetic Testing of Movements Disorders: A Review of Clinical Utility

Currently, pathogenic variants in more than 500 different genes are known to cause various movement disorders. The increasing accessibility and reducing cost of genetic testing has resulted in increasing clinical use of genetic testing for the diagnosis of movement disorders. However, the optimal use case(s) for genetic testing at a patient level remain ill-defined. Here, we review the utility of genetic testing in patients with movement disorders and also highlight current challenges and limitations that need to be considered when making decisions about genetic testing in clinical practice.

Highlights

The utility of genetic testing extends across multiple clinical and non-clinical domains. Here we review different aspects of the utility of genetic testing for movement disorders and the numerous associated challenges and limitations. These factors should be weighed on a case-by-case basis when requesting genetic tests in clinical practice.

COV2Var, a function annotation database of SARS-CoV-2 genetic variation

The COVID-19 pandemic, caused by the coronavirus SARS-CoV-2, has resulted in the loss of millions of lives and severe global economic consequences. Every time SARS-CoV-2 replicates, the viruses acquire new mutations in their genomes. Mutations in SARS-CoV-2 genomes led to increased transmissibility, severe disease outcomes, evasion of the immune response, changes in clinical manifestations and reducing the efficacy of vaccines or treatments. To date, the multiple resources provide lists of detected mutations without key functional annotations. There is a lack of research examining the relationship between mutations and various factors such as disease severity, pathogenicity, patient age, patient gender, cross-species transmission, viral immune escape, immune response level, viral transmission capability, viral evolution, host adaptability, viral protein structure, viral protein function, viral protein stability and concurrent mutations. Deep understanding the relationship between mutation sites and these factors is crucial for advancing our knowledge of SARS-CoV-2 and for developing effective responses. To fill this gap, we built COV2Var, a function annotation database of SARS-CoV-2 genetic variation, available at http://biomedbdc.wchscu.cn/COV2Var/. COV2Var aims to identify common mutations in SARS-CoV-2 variants and assess their effects, providing a valuable resource for intensive functional annotations of common mutations among SARS-CoV-2 variants.

Will variants of uncertain significance still exist in 2030?

In 2020, the National Human Genome Research Institute (NHGRI) made ten "bold predictions," including that "the clinical relevance of all encountered genomic variants will be readily predictable, rendering the diagnostic designation 'variant of uncertain significance (VUS)' obsolete." We discuss the prospects for this prediction, arguing that many, if not most, VUS in coding regions will be resolved by 2030. We outline a confluence of recent changes making this possible, especially advances in the standards for variant classification that better leverage diverse types of evidence, improvements in computational variant effect predictor performance, scalable multiplexed assays of variant effect capable of saturating the genome, and data-sharing efforts that will maximize the information gained from each new individual sequenced and variant interpreted. We suggest that clinicians and researchers can realize a future where VUSs have largely been eliminated, in line with the NHGRI's bold prediction. The length of time taken to reach this future, and thus whether we are able to achieve the goal of largely eliminating VUSs by 2030, is largely a consequence of the choices made now and in the next few years. We believe that investing in eliminating VUSs is worthwhile, since their predominance remains one of the biggest challenges to precision genomic medicine.

Saturation genome editing-based functional evaluation and clinical classification of BRCA2 single nucleotide variants

Germline BRCA2 loss-of function (LOF) variants identified by clinical genetic testing predispose to breast, ovarian, prostate and pancreatic cancer. However, variants of uncertain significance (VUS) (n>4000) limit the clinical use of testing results. Thus, there is an urgent need for functional characterization and clinical classification of all BRCA2 variants. Here we report on comprehensive saturation genome editing-based functional characterization of 97% of all possible single nucleotide variants (SNVs) in the BRCA2 DNA Binding Domain hotspot for pathogenic missense variants that is encoded by exons 15 to 26. The assay was based on deep sequence analysis of surviving endogenously targeted haploid cells. A total of 7013 SNVs were characterized as functionally abnormal (n=955), intermediate/uncertain, or functionally normal (n=5224) based on 95% agreement with ClinVar known pathogenic and benign standards. Results were validated relative to batches of nonsense and synonymous variants and variants evaluated using a homology directed repair (HDR) functional assay. Breast cancer case-control association studies showed that pooled SNVs encoding functionally abnormal missense variants were associated with increased risk of breast cancer (odds ratio (OR) 3.89, 95%CI: 2.77-5.51). In addition, 86% of tumors associated with abnormal missense SNVs displayed loss of heterozygosity (LOH), whereas 26% of tumors with normal variants had LOH. The functional data were added to other sources of information in a ClinGen/ACMG/AMP-like model and 700 functionally abnormal SNVs, including 220 missense SNVs, were classified as pathogenic or likely pathogenic, while 4862 functionally normal SNVs, including 3084 missense SNVs, were classified as benign or likely benign. These classified variants can now be used for risk assessment and clinical care of variant carriers and the remaining functional scores can be used directly for clinical classification and interpretation of many additional variants.

Summary

Germline BRCA2 loss-of function (LOF) variants identified by clinical genetic testing predispose to several types of cancer. However, variants of uncertain significance (VUS) limit the clinical use of testing results. Thus, there is an urgent need for functional characterization and clinical classification of all BRCA2 variants to facilitate current and future clinical management of individuals with these variants. Here we show the results from a saturation genome editing (SGE) and functional analysis of all possible single nucleotide variants (SNVs) from exons 15 to 26 that encode the BRCA2 DNA Binding Domain hotspot for pathogenic missense variants. The assay was based on deep sequence analysis of surviving endogenously targeted human haploid HAP1 cells. The assay was calibrated relative to ClinVar known pathogenic and benign missense standards and 95% prevalence thresholds for functionally abnormal and normal variants were identified. Thresholds were validated based on nonsense and synonymous variants. SNVs encoding functionally abnormal missense variants were associated with increased risks of breast and ovarian cancer. The functional assay results were integrated into a ClinGen/ACMG/AMP-like model for clinical classification of the majority of BRCA2 SNVs as pathogenic/likely pathogenic or benign/likely benign. The classified variants can be used for improved clinical management of variant carriers.

Protocol for the saturation and multiplexing of genetic variants using CRISPR-Cas9

Here, we present a multiplexed assay for variant effect protocol to assess the functional impact of all possible genetic variations within a particular genomic region. We describe steps for saturation genome editing by designing and cloning of single-guide RNA (sgRNA). We then detail steps for nucleofection of sgRNA, testing drug response on variants, and amplification of genomic DNA for next-generation sequencing. For complete details on the use and execution of this protocol, please refer to Sahu et al.1.

Constructing and interpreting a large-scale variant effect map for an ultrarare disease gene: Comprehensive prediction of the functional impact of PSAT1 genotypes

Reduced activity of the enzymes encoded by PHGDH, PSAT1, and PSPH causes a set of ultrarare, autosomal recessive diseases known as serine biosynthesis defects. These diseases present in a broad phenotypic spectrum: at the severe end is Neu-Laxova syndrome, in the intermediate range are infantile serine biosynthesis defects with severe neurological manifestations and growth deficiency, and at the mild end is childhood disease with intellectual disability. However, L-serine supplementation, especially if started early, can ameliorate and in some cases even prevent symptoms. Therefore, knowledge of pathogenic variants can improve clinical outcomes. Here, we use a yeast-based assay to individually measure the functional impact of 1,914 SNV-accessible amino acid substitutions in PSAT. Results of our assay agree well with clinical interpretations and protein structure-function relationships, supporting the inclusion of our data as functional evidence as part of the ACMG variant interpretation guidelines. We use existing ClinVar variants, disease alleles reported in the literature and variants present as homozygotes in the primAD database to define assay ranges that could aid clinical variant interpretation for up to 98% of the tested variants. In addition to measuring the functional impact of individual variants in yeast haploid cells, we also assay pairwise combinations of PSAT1 alleles that recapitulate human genotypes, including compound heterozygotes, in yeast diploids. Results from our diploid assay successfully distinguish the genotypes of affected individuals from those of healthy carriers and agree well with disease severity. Finally, we present a linear model that uses individual allele measurements to predict the biallelic function of ~1.8 million allele combinations corresponding to potential human genotypes. Taken together, our work provides an example of how large-scale functional assays in model systems can be powerfully applied to the study of ultrarare diseases.

Classification of PTEN missense VUS through exascale simulations

Phosphatase and tensin homolog (PTEN), a tumor suppressor with dual phosphatase properties, is a key factor in PI3K/AKT signaling pathway. Pathogenic germline variation in PTEN can abrogate its ability to dephosphorylate, causing high cancer risk. Lack of functional evidence lets numerous PTEN variants be classified as variants of uncertain significance (VUS). Utilizing Molecular Dynamics (MD) simulations, we performed a thorough evaluation for 147 PTEN missense VUS, sorting them into 66 deleterious and 81 tolerated variants. Utilizing replica exchange molecular dynamic (REMD) simulations, we further assessed the variants situated in the catalytic core of PTEN's phosphatase domain and uncovered conformational alterations influencing the structural stability of the phosphatase domain. There was a high degree of agreement between our results and the variants classified by Variant Abundance by Massively Parallel Sequencing, saturation mutagenesis, multiplexed functional data and experimental assays. Our extensive analysis of PTEN missense VUS should benefit their clinical applications in PTEN-related cancer.

SIGNIFICANCE STATEMENT

Classification of PTEN variants affecting its lipid phosphatase activity is important for understanding the roles of PTEN variation in the pathogenesis of hereditary and sporadic malignancies. Of the 3000 variants identified in PTEN, 1296 (43%) were assigned as VUS. Here, we applied MD and REMD simulations to investigate the effects of PTEN missense VUS on the structural integrity of the PTEN phosphatase domain consisting the WPD, P and TI active sites. We classified a total of 147 missense VUS into 66 deleterious and 81 tolerated variants by referring to the control group comprising 54 pathogenic and 12 benign variants. The classification was largely in concordance with these classified by experimental approaches.

An Atlas of Variant Effects to understand the genome at nucleotide resolution

Sequencing has revealed hundreds of millions of human genetic variants, and continued efforts will only add to this variant avalanche. Insufficient information exists to interpret the effects of most variants, limiting opportunities for precision medicine and comprehension of genome function. A solution lies in experimental assessment of the functional effect of variants, which can reveal their biological and clinical impact. However, variant effect assays have generally been undertaken reactively for individual variants only after and, in most cases long after, their first observation. Now, multiplexed assays of variant effect can characterise massive numbers of variants simultaneously, yielding variant effect maps that reveal the function of every possible single nucleotide change in a gene or regulatory element. Generating maps for every protein encoding gene and regulatory element in the human genome would create an 'Atlas' of variant effect maps and transform our understanding of genetics and usher in a new era of nucleotide-resolution functional knowledge of the genome. An Atlas would reveal the fundamental biology of the human genome, inform human evolution, empower the development and use of therapeutics and maximize the utility of genomics for diagnosing and treating disease. The Atlas of Variant Effects Alliance is an international collaborative group comprising hundreds of researchers, technologists and clinicians dedicated to realising an Atlas of Variant Effects to help deliver on the promise of genomics.

Translating Genetic Discovery into a Mechanistic Understanding of Pediatric Movement Disorders: Lessons from Genetic Dystonias and Related Disorders

The era of next-generation sequencing has increased the pace of gene discovery in the field of pediatric movement disorders. Following the identification of novel disease-causing genes, several studies have aimed to link the molecular and clinical aspects of these disorders. This perspective presents the developing stories of several childhood-onset movement disorders, including paroxysmal kinesigenic dyskinesia, myoclonus-dystonia syndrome, and other monogenic dystonias. These stories illustrate how gene discovery helps focus the research efforts of scientists trying to understand the mechanisms of disease. The genetic diagnosis of these clinical syndromes also helps clarify the associated phenotypic spectra and aids the search for additional disease-causing genes. Collectively, the findings of previous studies have led to increased recognition of the role of the cerebellum in the physiology and pathophysiology of motor control-a common theme in many pediatric movement disorders. To fully exploit the genetic information garnered in the clinical and research arenas, it is crucial that corresponding multi-omics analyses and functional studies also be performed at scale. Hopefully, these integrated efforts will provide us with a more comprehensive understanding of the genetic and neurobiological bases of movement disorders in childhood.

The functional impact of 1,570 individual amino acid substitutions in human OTC

Deleterious mutations in the X-linked gene encoding ornithine transcarbamylase (OTC) cause the most common urea cycle disorder, OTC deficiency. This rare but highly actionable disease can present with severe neonatal onset in males or with later onset in either sex. Individuals with neonatal onset appear normal at birth but rapidly develop hyperammonemia, which can progress to cerebral edema, coma, and death, outcomes ameliorated by rapid diagnosis and treatment. Here, we develop a high-throughput functional assay for human OTC and individually measure the impact of 1,570 variants, 84% of all SNV-accessible missense mutations. Comparison to existing clinical significance calls, demonstrated that our assay distinguishes known benign from pathogenic variants and variants with neonatal onset from late-onset disease presentation. This functional stratification allowed us to identify score ranges corresponding to clinically relevant levels of impairment of OTC activity. Examining the results of our assay in the context of protein structure further allowed us to identify a 13 amino acid domain, the SMG loop, whose function appears to be required in human cells but not in yeast. Finally, inclusion of our data as PS3 evidence under the current ACMG guidelines, in a pilot reclassification of 34 variants with complete loss of activity, would change the classification of 22 from variants of unknown significance to clinically actionable likely pathogenic variants. These results illustrate how large-scale functional assays are especially powerful when applied to rare genetic diseases.

High-throughput functional mapping of variants in an arrhythmia gene, KCNE1, reveals novel biology

Background

KCNE1 encodes a 129-residue cardiac potassium channel (IKs) 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 IKs channel biology and help reclassify variants of uncertain significance.

satmut_utils: a simulation and variant calling package for multiplexed assays of variant effect

The impact of millions of individual genetic variants on molecular phenotypes in coding sequences remains unknown. Multiplexed assays of variant effect (MAVEs) are scalable methods to annotate relevant variants, but existing software lacks standardization, requires cumbersome configuration, and does not scale to large targets. We present satmut_utils as a flexible solution for simulation and variant quantification. We then benchmark MAVE software using simulated and real MAVE data. We finally determine mRNA abundance for thousands of cystathionine beta-synthase variants using two experimental methods. The satmut_utils package enables high-performance analysis of MAVEs and reveals the capability of variants to alter mRNA abundance.

Accurate Prediction of Transcriptional Activity of Single Missense Variants in HIV Tat with Deep Learning

Tat is an essential gene for increasing the transcription of all HIV genes, and affects HIV replication, HIV exit from latency, and AIDS progression. The Tat gene frequently mutates in vivo and produces variants with diverse activities, contributing to HIV viral heterogeneity as well as drug-resistant clones. Thus, identifying the transcriptional activities of Tat variants will help to better understand AIDS pathology and treatment. We recently reported the missense mutation landscape of all single amino acid Tat variants. In these experiments, a fraction of double missense alleles exhibited intragenic epistasis. However, it is too time-consuming and costly to determine the effect of the variants for all double mutant alleles through experiments. Therefore, we propose a combined GigaAssay/deep learning approach. As a first step to determine activity landscapes for complex variants, we evaluated a deep learning framework using previously reported GigaAssay experiments to predict how transcription activity is affected by Tat variants with single missense substitutions. Our approach achieved a 0.94 Pearson correlation coefficient when comparing the predicted to experimental activities. This hybrid approach can be extensible to more complex Tat alleles for a better understanding of the genetic control of HIV genome transcription.

HAP1, a new revolutionary cell model for gene editing using CRISPR-Cas9

The use of next-generation sequencing (NGS) technologies has been instrumental in the characterization of the mutational landscape of complex human diseases like cancer. But despite the enormous rise in the identification of disease candidate genetic variants, their functionality is yet to be fully elucidated in order to have a clear implication in patient care. Haploid human cell models have become the tool of choice for functional gene studies, since they only contain one copy of the genome and can therefore show the unmasked phenotype of genetic variants. Over the past few years, the human near-haploid cell line HAP1 has widely been consolidated as one of the favorite cell line models for functional genetic studies. Its rapid turnover coupled with the fact that only one allele needs to be modified in order to express the subsequent desired phenotype has made this human cell line a valuable tool for gene editing by CRISPR-Cas9 technologies. This review examines the recent uses of the HAP1 cell line model in functional genetic studies and high-throughput genetic screens using the CRISPR-Cas9 system. It covers its use in an attempt to develop new and relevant disease models to further elucidate gene function, and create new ways to understand the genetic basis of human diseases. We will cover the advantages and potential of the use of CRISPR-Cas9 technology on HAP1 to easily and efficiently study the functional interpretation of gene function and human single-nucleotide genetic variants of unknown significance identified through NGS technologies, and its implications for changes in clinical practice and patient care.

Opportunities and Challenges with Artificial Intelligence in Genomics

The development of artificial intelligence and machine learning algorithms may allow for advances in patient care. There are existing and potential applications in cancer diagnosis and monitoring, identification of at-risk groups of individuals, classification of genetic variants, and even prediction of patient ancestry. This article provides an overview of some current and future applications of artificial intelligence in genomic medicine, in addition to discussing challenges and considerations when bringing these tools into clinical practice.

Predicting the functional effect of compound heterozygous genotypes from large scale variant effect maps

Background

Pathogenic variants in PHGDH, PSAT1 , and PSPH cause a set of rare, autosomal recessive diseases known as serine biosynthesis defects. Serine biosynthesis defects present in a broad phenotypic spectrum that includes, at the severe end, Neu-Laxova syndrome, a lethal multiple congenital anomaly disease, intermediately in the form of infantile serine biosynthesis defects with severe neurological manifestations and growth deficiency, and at the mild end, as childhood disease with intellectual disability. However, because L-serine supplementation, especially if started early, can ameliorate and in some cases even prevent symptoms, knowledge of pathogenic variants is highly actionable.

Methods

Recently, our laboratory established a yeast-based assay for human PSAT1 function. We have now applied it at scale to assay the functional impact of 1,914 SNV-accessible amino acid substitutions. In addition to assaying the functional impact of individual variants in yeast haploid cells, we can assay pairwise combinations of PSAT1 alleles that recapitulate human genotypes, including compound heterozygotes, in yeast diploids.

Results

Results of our assays of individual variants (in haploid yeast cells) agree well with clinical interpretations and protein structure-function relationships, supporting the use of our data as functional evidence under the ACMG interpretation guidelines. Results from our diploid assay successfully distinguish patient genotypes from those of healthy carriers and agree well with disease severity. Finally, we present a linear model that uses individual allele measurements (in haploid yeast cells) to accurately predict the biallelic function (in diploid yeast cells) of ~ 1.8 million allele combinations corresponding to potential human genotypes.

Conclusions

Taken together, our work provides an example of how large-scale functional assays in model systems can be powerfully applied to the study of a rare disease.

Deep mutational scanning: A versatile tool in systematically mapping genotypes to phenotypes

Unveiling how genetic variations lead to phenotypic variations is one of the key questions in evolutionary biology, genetics, and biomedical research. Deep mutational scanning (DMS) technology has allowed the mapping of tens of thousands of genetic variations to phenotypic variations efficiently and economically. Since its first systematic introduction about a decade ago, we have witnessed the use of deep mutational scanning in many research areas leading to scientific breakthroughs. Also, the methods in each step of deep mutational scanning have become much more versatile thanks to the oligo-synthesizing technology, high-throughput phenotyping methods and deep sequencing technology. However, each specific possible step of deep mutational scanning has its pros and cons, and some limitations still await further technological development. Here, we discuss recent scientific accomplishments achieved through the deep mutational scanning and describe widely used methods in each step of deep mutational scanning. We also compare these different methods and analyze their advantages and disadvantages, providing insight into how to design a deep mutational scanning study that best suits the aims of the readers' projects.

Shared Cancer Dataset Analysis Identifies and Predicts the Quantitative Effects of Pan-Cancer Somatic Driver Variants

Driver mutations endow tumors with selective advantages and produce an array of pathogenic effects. Determining the function of somatic variants is important for understanding cancer biology and identifying optimal therapies. Here, we compiled a shared dataset from several cancer genomic databases. Two measures were applied to 535 cancer genes based on observed and expected frequencies of driver variants as derived from cancer-specific rates of somatic mutagenesis. The first measure comprised a binary classifier based on a binomial test; the second was tumor variant amplitude (TVA), a continuous measure representing the selective advantage of individual variants. TVA outperformed all other computational tools in terms of its correlation with experimentally derived functional scores of cancer mutations. TVA also highly correlated with drug response, overall survival, and other clinical implications in relevant cancer genes. This study demonstrates how a selective advantage measure based on a large cancer dataset significantly impacts our understanding of the spectral effect of driver variants in cancer. The impact of this information will increase as cancer treatment becomes more precise and personalized to tumor-specific mutations.

SIGNIFICANCE

A new selective advantage estimation assists in oncogenic driver identification and relative effect measurements, enabling better prognostication, therapy selection, and prioritization.

RosettaDDGPrediction for high-throughput mutational scans: From stability to binding

Reliable prediction of free energy changes upon amino acid substitutions (ΔΔGs) is crucial to investigate their impact on protein stability and protein-protein interaction. Advances in experimental mutational scans allow high-throughput studies thanks to multiplex techniques. On the other hand, genomics initiatives provide a large amount of data on disease-related variants that can benefit from analyses with structure-based methods. Therefore, the computational field should keep the same pace and provide new tools for fast and accurate high-throughput ΔΔG calculations. In this context, the Rosetta modeling suite implements effective approaches to predict folding/unfolding ΔΔGs in a protein monomer upon amino acid substitutions and calculate the changes in binding free energy in protein complexes. However, their application can be challenging to users without extensive experience with Rosetta. Furthermore, Rosetta protocols for ΔΔG prediction are designed considering one variant at a time, making the setup of high-throughput screenings cumbersome. For these reasons, we devised RosettaDDGPrediction, a customizable Python wrapper designed to run free energy calculations on a set of amino acid substitutions using Rosetta protocols with little intervention from the user. Moreover, RosettaDDGPrediction assists with checking completed runs and aggregates raw data for multiple variants, as well as generates publication-ready graphics. We showed the potential of the tool in four case studies, including variants of uncertain significance in childhood cancer, proteins with known experimental unfolding ΔΔGs values, interactions between target proteins and disordered motifs, and phosphomimetics. RosettaDDGPrediction is available, free of charge and under GNU General Public License v3.0, at https://github.com/ELELAB/RosettaDDGPrediction.

Integrating deep mutational scanning and low-throughput mutagenesis data to predict the impact of amino acid variants

BACKGROUND

Evaluating the impact of amino acid variants has been a critical challenge for studying protein function and interpreting genomic data. High-throughput experimental methods like deep mutational scanning (DMS) can measure the effect of large numbers of variants in a target protein, but because DMS studies have not been performed on all proteins, researchers also model DMS data computationally to estimate variant impacts by predictors.

RESULTS

In this study, we extended a linear regression-based predictor to explore whether incorporating data from alanine scanning (AS), a widely used low-throughput mutagenesis method, would improve prediction results. To evaluate our model, we collected 146 AS datasets, mapping to 54 DMS datasets across 22 distinct proteins.

CONCLUSIONS

We show that improved model performance depends on the compatibility of the DMS and AS assays, and the scale of improvement is closely related to the correlation between DMS and AS results.

Saturation-scale functional evidence supports clinical variant interpretation in Lynch syndrome

BACKGROUND

Lynch syndrome (LS) is a cancer predisposition syndrome affecting more than 1 in every 300 individuals worldwide. Clinical genetic testing for LS can be life-saving but is complicated by the heavy burden of variants of uncertain significance (VUS), especially missense changes.

RESULT

To address this challenge, we leverage a multiplexed analysis of variant effect (MAVE) map covering >94% of the 17,746 possible missense variants in the key LS gene MSH2. To establish this map's utility in large-scale variant reclassification, we overlay it on clinical databases of >15,000 individuals with LS gene variants uncovered during clinical genetic testing. We validate these functional measurements in a cohort of individuals with paired tumor-normal test results and find that MAVE-based function scores agree with the clinical interpretation for every one of the MSH2 missense variants with an available classification. We use these scores to attempt reclassification for 682 unique missense VUS, among which 34 scored as deleterious by our function map, in line with previously published rates for other cancer predisposition genes. Combining functional data and other evidence, ten missense VUS are reclassified as pathogenic/likely pathogenic, and another 497 could be moved to benign/likely benign. Finally, we apply these functional scores to paired tumor-normal genetic tests and identify a subset of patients with biallelic somatic loss of function, reflecting a sporadic Lynch-like Syndrome with distinct implications for treatment and relatives' risk.

CONCLUSION

This study demonstrates how high-throughput functional assays can empower scalable VUS resolution and prospectively generate strong evidence for variant classification.

Reducing uncertainty in genetic testing with Saturation Genome Editing

Accurate interpretation of human genetic data is critical for optimizing outcomes in the era of genomic medicine. Powerful methods for testing genetic variants for functional effects are allowing researchers to characterize thousands of variants across disease genes. Here, we review experimental tools enabling highly scalable assays of variants, focusing specifically on Saturation Genome Editing (SGE). We discuss examples of how this technique is being implemented for variant testing at scale and describe how SGE data for BRCA1 have been clinically validated and used to aid variant interpretation. The initial success at predicting variant pathogenicity with SGE has spurred efforts to expand this and related techniques to many more genes.

Scalable Functional Assays for the Interpretation of Human Genetic Variation

Scalable sequence-function studies have enabled the systematic analysis and cataloging of hundreds of thousands of coding and noncoding genetic variants in the human genome. This has improved clinical variant interpretation and provided insights into the molecular, biophysical, and cellular effects of genetic variants at an astonishing scale and resolution across the spectrum of allele frequencies. In this review, we explore current applications and prospects for the field and outline the principles underlying scalable functional assay design, with a focus on the study of single-nucleotide coding and noncoding variants.

Pooled image-base screening of mitochondria with microraft isolation distinguishes pathogenic mitofusin 2 mutations

Most human genetic variation is classified as variants of uncertain significance. While advances in genome editing have allowed innovation in pooled screening platforms, many screens deal with relatively simple readouts (viability, fluorescence) and cannot identify the complex cellular phenotypes that underlie most human diseases. In this paper, we present a generalizable functional genomics platform that combines high-content imaging, machine learning, and microraft isolation in a method termed "Raft-Seq". We highlight the efficacy of our platform by showing its ability to distinguish pathogenic point mutations of the mitochondrial regulator Mitofusin 2, even when the cellular phenotype is subtle. We also show that our platform achieves its efficacy using multiple cellular features, which can be configured on-the-fly. Raft-Seq enables a way to perform pooled screening on sets of mutations in biologically relevant cells, with the ability to physically capture any cell with a perturbed phenotype and expand it clonally, directly from the primary screen.

Intragenic compensation through the lens of deep mutational scanning

A significant fraction of mutations in proteins are deleterious and result in adverse consequences for protein function, stability, or interaction with other molecules. Intragenic compensation is a specific case of positive epistasis when a neutral missense mutation cancels effect of a deleterious mutation in the same protein. Permissive compensatory mutations facilitate protein evolution, since without them all sequences would be extremely conserved. Understanding compensatory mechanisms is an important scientific challenge at the intersection of protein biophysics and evolution. In human genetics, intragenic compensatory interactions are important since they may result in variable penetrance of pathogenic mutations or fixation of pathogenic human alleles in orthologous proteins from related species. The latter phenomenon complicates computational and clinical inference of an allele's pathogenicity. Deep mutational scanning is a relatively new technique that enables experimental studies of functional effects of thousands of mutations in proteins. We review the important aspects of the field and discuss existing limitations of current datasets. We reviewed ten published DMS datasets with quantified functional effects of single and double mutations and described rates and patterns of intragenic compensation in eight of them.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-022-01005-w.

High-throughput approaches to functional characterization of genetic variation in yeast

Expansion of sequencing efforts to include thousands of genomes is providing a fundamental resource for determining the genetic diversity that exists in a population. Now, high-throughput approaches are necessary to begin to understand the role these genotypic changes play in affecting phenotypic variation. Saccharomyces cerevisiae maintains its position as an excellent model system to determine the function of unknown variants with its exceptional genetic diversity, phenotypic diversity, and reliable genetic manipulation tools. Here, we review strategies and techniques developed in yeast that scale classic approaches of assessing variant function. These approaches improve our ability to better map quantitative trait loci at a higher resolution, even for rare variants, and are already providing greater insight into the role that different types of mutations play in phenotypic variation and evolution not just in yeast but across taxa.

Embeddings from protein language models predict conservation and variant effects

The emergence of SARS-CoV-2 variants stressed the demand for tools allowing to interpret the effect of single amino acid variants (SAVs) on protein function. While Deep Mutational Scanning (DMS) sets continue to expand our understanding of the mutational landscape of single proteins, the results continue to challenge analyses. Protein Language Models (pLMs) use the latest deep learning (DL) algorithms to leverage growing databases of protein sequences. These methods learn to predict missing or masked amino acids from the context of entire sequence regions. Here, we used pLM representations (embeddings) to predict sequence conservation and SAV effects without multiple sequence alignments (MSAs). Embeddings alone predicted residue conservation almost as accurately from single sequences as ConSeq using MSAs (two-state Matthews Correlation Coefficient-MCC-for ProtT5 embeddings of 0.596 ± 0.006 vs. 0.608 ± 0.006 for ConSeq). Inputting the conservation prediction along with BLOSUM62 substitution scores and pLM mask reconstruction probabilities into a simplistic logistic regression (LR) ensemble for Variant Effect Score Prediction without Alignments (VESPA) predicted SAV effect magnitude without any optimization on DMS data. Comparing predictions for a standard set of 39 DMS experiments to other methods (incl. ESM-1v, DeepSequence, and GEMME) revealed our approach as competitive with the state-of-the-art (SOTA) methods using MSA input. No method outperformed all others, neither consistently nor statistically significantly, independently of the performance measure applied (Spearman and Pearson correlation). Finally, we investigated binary effect predictions on DMS experiments for four human proteins. Overall, embedding-based methods have become competitive with methods relying on MSAs for SAV effect prediction at a fraction of the costs in computing/energy. Our method predicted SAV effects for the entire human proteome (~ 20 k proteins) within 40 min on one Nvidia Quadro RTX 8000. All methods and data sets are freely available for local and online execution through bioembeddings.com, https://github.com/Rostlab/VESPA , and PredictProtein.

Mutational scanning identified amino acids of the CLOCK exon 19-domain essential for circadian rhythms

AIM

In the mammalian circadian clock, the CLOCK/BMAL1 heterodimer binds to E-box enhancer elements in the promoters of its target genes to activate transcription. The classical Clock mice, the first circadian mouse mutant discovered, are behaviourally arrhythmic. In this mutant, CLOCK lacks a 51 amino acid domain corresponding to exon 19 (CLOCKΔ19), which is required for normal transactivation. While the importance of this CLOCK domain for circadian rhythms is well established, the exact molecular mechanism is still unclear.

METHODS

Using CRISPR/Cas9 technology, we created a CLOCK knockout - CLOCK rescue system in human circadian reporter cells and performed systematic mutational scanning to assess the functionality of individual amino acids within the CLOCK exon 19-domain.

RESULTS

CLOCK knockout cells were arrhythmic, and circadian rhythms could be rescued by introducing wild-type CLOCK, but not CLOCKΔ19. In addition, we identified several residues, whose mutation failed to rescue rhythms in CLOCK knockout cells. Many of these are part of the hydrophobic binding interface of the predicted dimer of the CLOCK exon 19-domain.

CONCLUSION

Our data not only indicate that CLOCK/BMAL1 oligomerization mediated by the exon 19-domain is important for circadian dynamics but also suggest that the exon 19-domain provides a platform for binding coactivators and repressors, which in turn is required for normal circadian rhythms.

Machine Learning Models for Accurate Prioritization of Variants of Uncertain Significance

The growing use of next generation sequencing technologies on genetic diagnosis has produced an exponential increase in the number of Variants of Uncertain Significance (VUS). In this manuscript we compare three machine learning methods to classify VUS as Pathogenic or No pathogenic, implementing a Random Forest (RF), a Support Vector Machine (SVM), and a Multilayer Perceptron (MLP). To train the models, we extracted high quality variants from ClinVar that were previously classified as VUS. For each variant, we retrieved 9 conservation scores, the loss of function tool and allele frequencies. For the RF and SVM models, hyperparameters were tuned using cross validation with a grid search. The three models were tested on a non-overlapping set of variants that had been classified as VUS any time along the last three years but had been reclassified in august 2020. The three models yielded superior accuracy on this set compared to the benchmarked tools. The RF based model yielded the best performance across different variant types and was used to create VusPrize, an open source software tool for prioritization of variants of uncertain significance. We believe that our model can improve the process of genetic diagnosis in research and clinical settings. This article is protected by copyright. All rights reserved.

Variant Library Annotation Tool (VaLiAnT): an oligonucleotide library design and annotation tool for saturation genome editing and other deep mutational scanning experiments

MOTIVATION

CRISPR/Cas9-based technology allows for the functional analysis of genetic variants at single nucleotide resolution whilst maintaining genomic context. This approach, known as saturation genome editing (SGE), a form of deep mutational scanning, systematically alters each position in a target region to explore its function. SGE experiments require the design and synthesis of oligonucleotide variant libraries which are introduced into the genome. This technology is applicable to diverse fields such as disease variant identification, drug development, structure-function studies, synthetic biology, evolutionary genetics and host-pathogen interactions. Here, we present the Variant Library Annotation Tool (VaLiAnT) which can be used to generate variant libraries from user-defined genomic coordinates and standard input files. The software can accommodate user-specified species, reference sequences and transcript annotations.

RESULTS

Coordinates for a genomic range are provided by the user to retrieve a corresponding oligonucleotide reference sequence. A user-specified range within this sequence is then subject to systematic, nucleotide and/or amino acid saturating mutator functions. VaLiAnT provides a novel way to retrieve, mutate and annotate genomic sequences for oligonucleotide library generation. Specific features for SGE library generation can be employed. In addition, VaLiAnT is configurable, allowing for cDNA and prime editing saturation library generation, with other diverse applications possible.

AVAILABILITY AND IMPLEMENTATION

VaLiAnT is a command line tool written in Python. Source code, testing data, example input and output files and executables are available (https://github.com/cancerit/VaLiAnT) in addition to a detailed user manual (https://github.com/cancerit/VaLiAnT/wiki). VaLiAnT is licensed under AGPLv3.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.