Single-nucleotide evolutionary constraint scores highlight disease-causing mutations

By Gregory M. Cooper, David L. Goode, Sarah B. H. Ng, Arend Sidow, Michael Bamshad, Jay Shendure, Deborah A. Nickerson. 2010

208 citations recorded by the source corpus.

Identifying disease-causing genetic variants in individual human genomes is a major challenge, even in protein-coding exons (the `exome'). Analysis of nucleotide-level sequence conservation may help address this challenge, on the assumption that purifying selection `constrains' evolutionary divergence at phenotypically important nucleotides. In contrast to functional classifiers (for example, non-synonymous mutations), constraint scores are quantitative and applicable to any genomic position1. However, it remains unclear if constraint scores can facilitate causal variant discovery, as statistical power is estimated to be marginal at the single-nucleotide level given current genome alignments1,2. We therefore applied and assessed a nucleotide-level evolutionary metric to prioritize causal variants in genomes of 16 individuals. We analyzed exomes from four individuals with Freeman-Sheldon syndrome (FSS; Online Mendelian Inheritance in Man (OMIM) database identifier 193700), a dominant disease caused by mutations in MYH33; four individuals with Miller syndrome (OMIM identifier 263750), a recessive disease caused by mutations in DHODH4; and eight HapMap samples3. We generated constraint scores by Genomic Evolutionary Rate Profiling (GERP)1 on the mammalian subset of the 44-way MULTIZ/TBA alignments (http://genome.ucsc.edu; for details see2). For each aligned site, GERP defines a `rejected substitution' (RS) score by estimating the actual number of substitutions at that site and subtracting it from the number expected assuming neutrality (~5.82 substitutions per site). Selectively constrained sites tolerate fewer substitutions than neutral sites and have positive RS scores1,2. We first defined the consensus nucleotide from the chimpanzee, gorilla, orangutan, and macaque genomes as ancestral and determined the derived allele frequency (DAF) for each variant in the eight HapMap exomes. We found a significant inverse correlation between DAF and RS score (Fig. 1; p < 0.0001; R2 = 1.4%, slope estimated as −1% DAF per RS). No correlation existed between DAF and the RS score for the nucleotide adjacent to the variant (Supplementary Fig. 1). While the DAF-RS correlation resulted partly from enrichment for singletons at sites with high RS scores, it was significant even within common variants (p < 0.0001; Supplementary Fig. 2). We also found that segregating sites, regardless of DAF, were enriched at sites with low RS scores and progressively depleted as the RS scores increased (Supplementary Fig. 3). Consistent with previous data2, these results suggest that RS scores enrich site-specifically for deleterious variants and non-variant positions at which new mutations would be deleterious. Figure 1 RS scores inversely correlate with DAF of single-nucleotide variants (n=48,750) in 8 HapMap exomes. The average DAF (Y-axis) is plotted for all variants at a site within a given RS bin (X-axis). Error bars show 1 standard error unit. Next, we tested whether constraint scores could enrich for FSS or Miller syndrome causal variants. We identified candidate disease genes as those in which the affected individuals had variants not seen in the HapMap exomes that affected a nucleotide with a high constraint score. For a comparison, we used functional definitions of deleteriousness, namely non-synonymous, splice-site, or insertion-deletion (indel)3,4. We first used a threshold of RS > 0 (fewer substitutions than expected), and found that this narrowed candidate gene lists nearly as effectively as functional annotations. For example, there were 21 genes in which all FSS samples had a rare, functionally annotated variant3,4 versus 24 genes in which all FSS samples had a rare variant with RS > 0 (Fig. 2a). Increasing the RS threshold, which cannot readily be done with functional annotations, reduced candidate gene lists. At a threshold of RS > 4, for example, MYH3 was one of only five FSS candidate genes, while DHODH was the only Miller syndrome candidate (Fig. 2b). Figure 2 Constraint scores enrich for disease-causing genes. (a). Number of genes (Y-axis, log-scale) in which at least the given number of FSS individuals (X-axis) has a rare variant that is functionally defined (white), or with increasing RS scores (light gray ... We note that protein-based approaches could similarly be used to reduce candidate gene lists. For example, there were only seven genes in which all FSS individuals harbored a rare variant annotated by PolyPhen5 as `possibly' or `probably' damaging. However, PolyPhen (and related approaches) is restricted to non-synonymous variants, does not facilitate ranking of candidates (see below), and excluded DHODH as a Miller syndrome candidate4. Finally, we exploited the quantitative nature of constraint scores and ranked genes by the average score of all rare and deleterious (RS > 0) variants in the affected individuals. MYH3 and DHODH ranked highly for their associated diseases, even under models allowing for the possibility of multiple causal genes. For example, requiring only that at least two individuals shared the same causal gene, MYH3 ranked 9th among 666 genes. If we assumed FSS and Miller syndrome to be monogenic, MYH3 and DHODH ranked as the top candidates, respectively (Fig. 2). RS scores for known or user-defined variants can be obtained from the Genome Variation Server (http://gvs.gs.washington.edu/GVS/) or SeattleSeq annotation pipeline (http://gvs.gs.washington.edu/SeattleSeqAnnotation/). Constraint scores facilitate threshold flexibility and candidate ranking and do not require functional annotations. Even in exomes, this allows for the possibility that synonymous variants contribute to disease6. More importantly, this independence offers exciting potential for the discovery of causal variation in arbitrary genomic segments (for example, linkage peaks) and ultimately resequenced genomes.

Subfields: Genetics, Molecular Biology

How researchers see this paper

The researchers whose trusted networks rank it highest, strongest advocate first. Every researcher gets a different order.

In Genetics

In Molecular Biology

Citations in the corpus

CADD: predicting the deleteriousness of variants throughout the human genome (2018, positive): “…prioritize variants based on functional or evolutionary properties such as sequence conservation, genic effects and regulatory element annotations can serve to improve power and ultimately the success of disease studies, for both Mendelian phenotypes (3) as well as common traits and diseases (4).”
A general framework for estimating the relative pathogenicity of human genetic variants (2014, background): “Second, it is well established that annotations like PolyPhen and conservation are valuable in the sequencing-based identification of disease-causal genes by virtue of their ability to highly rank pathogenic variants 1,2,38 .”
dbNSFP: A Lightweight Database of Human Non-synonymous SNPs and Their Functional Predictions (2011, background): “More importantly, functional prediction has become an indispensable step in identifying genes causing rare Mendelian diseases using an exome-sequencing approach [Cooper et al., 2010; Ng et al., 2010a,b,c].”
Evaluation of in silico algorithms for use with ACMG/AMP clinical variant interpretation guidelines (2017, methodology): “All the largely evolutionary conservation algorithms such as phyloP[25], phastCons[26], GERP[27] and Siphy[28] belong to different clusters from the metapredictors REVEL, MetaSVM and MetaLR (Fig.”
Exome sequencing identifies MLL2 mutations as a cause of Kabuki syndrome (2010, unknown): “We have previously shown that the quantitative ranking of candidate genes by the rejected substitution scores of their variants can facilitate the exome-based analysis of Mendelian disorder”
Needles in stacks of needles: finding disease-causal variants in a wealth of genomic data (2011, unknown): “Recently, we showed that nucleotide constraint scores defined by GERP 67 provide enrichments for deleterious causal mutations that are similar to those provided by protein-based annotations for two Mendelian diseases 72 .”
Exome sequencing in sporadic autism spectrum disorders identifies severe de novo mutations (2011, unknown): “…were enriched in the aggregate for disruptive events by considering two independent quantitative measures: the nature of the amino-acid replacement (Grantham matrix score 17 ) and the degree of nucleotide-level evolutionary conservation (Genomic Evolutionary Rate Profiling (GERP) 18,19 ) (Fig.”
Comparison and integration of deleteriousness prediction methods for nonsynonymous SNVs in whole exome sequencing studies. (2015, methodology): “Multiple algorithms were developed for predicting such deleteriousness based on different information of the variant, such as its sequence homology (2), protein structure (3,4) and evolutionary conservation (5).”
Exome sequencing as a tool for Mendelian disease gene discovery (2011, unknown): “Implementations of this strategy include phastCONS 39 , phyloP and Genetic Evolutionary Rate Profiling (GERP) 40 .”
dbNSFP v2.0: A Database of Human Non-synonymous SNVs and Their Functional Predictions and Annotations (2013, methodology): “Although GERP RS scores were typically used to measure the conservation of a nucleotide site in Mendelian disease studies (e.g. Cooper et al. 2010), an alternative measure might be a scaled RS score with the corresponding neutral rate (NR) of the site (i.e. RS/NR ratio).”