The Genome Aggregation Database (gnomAD) - Predicted Constraint Metrics track set contains
metrics of pathogenicity per-gene as predicted for gnomAD v2.1.1 and identifies genes subject to
strong selection against various classes of mutation.
This track includes several subtracks of constraint metrics calculated at gene (canonical
transcript), transcript and transcript-region level. For more information see the following
The metrics include:
- Observed and expected variant counts per transcript/gene
- Observed/Expected ratio (O/E)
- Z-scores of the observed counts compared to expected
- Probability of loss of function intolerance (pLI), for predicted loss-of-function (pLoF) variation only
- Chi-Squared difference of observed to expected counts, for the regional missense constraint track only
Display Conventions and Configuration
There are two "groups" of tracks in this set:
All tracks follow the general configuration settings for bigBed tracks. Mouseover on the
Gene/Transcript Constraint tracks shows the pLI score and the loss-of-function
observed/expected upper bound fraction (LOEUF), while mouseover on the Regional
Constraint track shows only the missense O/E ratio. Clicking on items in any track brings
up a table of constraint metrics.
- Gene/Transcript Constraint tracks: Predicted constraint metrics at the whole gene
level or whole transcript level for three different types of variation: missense, synonymous,
and predicted loss-of-function. The Gene Constraint track displays metrics for a canonical
transcript per gene defined as the longest isoform. The Transcript Constraint track displays
metrics for all transcript isoforms. Items on both tracks are shaded according to the pLI score,
with outlier items shaded in grey.
- Regional Constraint track: Missense-variation constrained regions at the sub-genic
level. This track displays metrics for transcripts that have two or more regions with
significantly different levels of missense constraint. All items are colored black.
Clicking the grey box to the left of the track, or right-clicking and choosing the Configure option,
brings up the interface for filtering items based on their pLI score, or labeling the items
based on their Ensembl identifier and/or Gene Name.
Please see the gnomAD browser help page and FAQ for further explanation of the topics below.
Observed and Expected Variant Counts
Observed count: The number of unique single-nucleotide variants in each transcript/gene
with 123 or fewer alternative alleles (MAF < 0.1%).
Expected count: A depth corrected probability prediction model that takes into account
sequence context, coverage, and methylation, was used to predict expected
variant counts. For more information please see Lek et al., 2016.
Variants found in exons with a median depth < 1 were removed from both counts.
The O/E constraint score is the ratio of the observed/expected variants in that gene. Each item in
this track shows the O/E ratio for three different types of variation: missense, synonymous, and
loss-of-function. The O/E ratio is a continuous measurement of how tolerant a gene or
transcript is to a certain class of variation. When a gene has a low O/E value, it is under stronger
selection for that class of variation than a gene with a higher O/E value. Because Counts depend on
gene size and sample size, the precision of the values varies a lot from one gene to the next.
Therefore, the 90% confidence interval (CI) is also displayed along with the O/E ratio to better
assist interpretation of the scores.
When evaluating how constrained a gene is, it is essential to consider the CI when using O/E. In
research and clinical interpretation of Mendelian cases, pLI > 0.9 has widely been used for
filtering. Accordingly, the Gnomad team suggests using the upper bound of the O/E confidence interval
LOEUF < 0.35 as a threshold if needed.
Please see the Methods section below for more information about how the scores were calculated.
pLI and Z-scores
The pLI and Z-scores of the deviation of observed variant counts relative to the expected number
are intended to measure how constrained or intolerant a gene or transcript is to a specific type of
variation. Genes or transcripts that are particularly depleted of a specific class of variation
(as observed in the gnomAD data set) are considered intolerant of that specific type of variation.
Z-scores are available for the missense and synonynmous categories and pLI scores are available for
the loss-of-function variation.
NOTE: The Regional Constraint track data reflects regions within transcripts that are
intolerant of missense variation within the ExAc dataset and was calculated with the method
described by Samocha et al., 2017.
Missense and Synonymous: Positive Z-scores indicate more constraint (fewer observed
variants than expected), and negative scores indicate less constraint (more observed variants than
expected). A greater Z-score indicates more intolerance to the class of variation. Z-scores
were generated by a sequence-context-based mutational model that predicted the number of expected
rare (< 1% MAF) variants per transcript. The square root of the chi-squared value of the
deviation of observed counts from expected counts was multiplied by -1 if the observed count was
greater than the expected and vice versa. For the synonymous score, each Z-score was corrected by
dividing by the standard deviation of all synonymous Z-scores between -5 and 5. For the missense
scores, a mirrored distribution of all Z-scores between -5 and 0 was created, and then all missense
Z-scores were corrected by dividing by the standard deviation of the Z-score of the mirror
Loss-of-function: pLI closer to 1 indicates that the gene or transcript cannot tolerate
protein truncating variation (nonsense, splice acceptor and splice donor variation). The gnomAD
team recommends transcripts with a pLI >= 0.9 for the set of transcripts extremely intolerant
to truncating variants. pLI is based on the idea that transcripts can be classified into three
An expectation-maximization algorithm was then used to assign a probability of belonging in each
class to each gene or transcript. pLI is the probability of belonging in the haploinsufficient class.
- null: heterozygous or homozygous protein truncating variation is completely tolerated
- recessive: heterozygous variants are tolerated but homozygous variants are not
- haploinsufficient: heterozygous variants are not tolerated
Please see Samocha et al., 2014 and Lek et al., 2016 for further discussion of these metrics.
Transcripts from GENCODE v19 were filtered according to the following criteria:
After filtering the transcript set, 18225 transcripts were left over.
- Must have methionine at start of coding sequence
- Must have stop codon at end of coding sequence
- Must be divisible by 3
- Must have at least one observed variant when removing exons with median depth < 1
- Must have reasonable number of missense and synonymous variants as determined by a Z-score cutoff
UCSC Track Methods
Gene and Transcript Constraint tracks
Per gene and per transcript data were downloaded from the gnomAD Google Storage bucket:
These data were then joined to the Gencode v19 set of genes/transcripts available at the UCSC
Genome Browser and then transformed into a bigBed 12+5. For the full list of commands used to
make this track please see the "gnomAD 2 pLI and other loss-of-function metrics" section
Regional Constraint track
Supplementary Table 4 from the
associated publication was downloaded and joined to the Gencode v19 set of transcripts
available at UCSC and then transformed into a bigBed 12+6. For the full list of commands
used to make this track please the "gnomAD Missense Constraint Scores" section of the
The raw data can be explored interactively with the Table Browser, or
the Data Integrator. For automated access, this track, like all
others, is available via our API. However, for bulk
processing, it is recommended to download the dataset. The genome annotation is stored in a bigBed
file that can be downloaded from the
download server. The exact
filenames can be found in the track configuration file. Annotations can be converted to ASCII text
by our tool
bigBedToBed which can be compiled from the source code or downloaded as
a precompiled binary for your system. Instructions for downloading source code and binaries can be
found here. The tool
can also be used to obtain only features within a given range, for example:
bigBedToBed http://hgdownload.soe.ucsc.edu/gbdb/hg19/gnomAD/pLI/pliByTranscript.bb -chrom=chr6 -start=0 -end=1000000 stdout
Please refer to our
mailing list archives
for questions and example queries, or our
Data Access FAQ
for more information.
Thanks to the gnomAD team, especially Anne O'Donnell Luria and Emily O'Heir for their help with
making these data and recommendations on how to display the data, and to Christopher Lee, Lou
Nassar and Anna Benet-Pagès of the Browser team for implementation.
Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alfoldi J, Wang Q, Collins RL, Laricchia KM, Gann A, Birnbaum DP et al.
Variation across 141,456 human exomes and genomes reveals the spectrum of loss-of-function intolerance across human protein-coding genes.
bioRxiv. 2019 August 13.
Lek M, Karczewski KJ, Minikel EV, Samocha KE, Banks E, Fennell T, O'Donnell-Luria AH, Ware JS, Hill
AJ, Cummings BB et al.
Analysis of protein-coding genetic variation in 60,706 humans.
Nature. 2016 Aug 18;536(7616):285-91.
PMID: 27535533; PMC: PMC5018207
Samocha KE, Robinson EB, Sanders SJ, Stevens C, Sabo A, McGrath LM, Kosmicki JA, Rehnström K,
Mallick S, Kirby A et al.
A framework for the interpretation of de novo mutation in human disease.
Nat Genet. 2014 Sep;46(9):944-50.
PMID: 25086666; PMC: PMC4222185
Samocha KE, Kosmicki JA, Karczewski KJ, O'Donnell-Luria AH, Pierce-Hoffman E, MacArthur DG, Neale BM, Daly MJ.
Regional missense constraint improves variant deleteriousness prediction.
bioRxiv. 2017 June 12.