Introduction
ClinVar is a
freely available, public archive of human genetic variants that provides
clinical classifications for whether a variant is likely benign or
pathogenic. The AlphaMissense publication
uses the ClinVar data to evaluate and calibrate the predictions
generated by their model. A table containing ClinVar information for
82872 variants across 7951 proteins was derived from the supplemental
data of the AlphaMissense paper, and is made available through this
package for benchmarking and visualization purposes.
Access ClinVar classifications with AlphaMissense predictions
The ClinVar table can be accessed using clinvar_data()
from the database.
clinvar_data()
#> * [16:30:25][info] downloading or finding local file
#> * [16:30:25][info] creating database table 'clinvar'
#> * [16:30:25][info] disconnecting all registered connections
#> # Source: table<clinvar> [?? x 5]
#> # Database: DuckDB v1.1.1 [biocbuild@Linux 6.8.0-47-generic:R 4.5.0//home/biocbuild/.cache/R/BiocFileCache/2f7e82265d846_2f7e82265d846]
#> variant_id transcript_id protein_variant AlphaMissense label
#> <chr> <chr> <chr> <dbl> <fct>
#> 1 chr1_925969_C_T_hg38 ENST00000342066.8 Q96NU1:P10S 0.967 benign
#> 2 chr1_930165_G_A_hg38 ENST00000342066.8 Q96NU1:R28Q 0.663 benign
#> 3 chr1_930204_G_A_hg38 ENST00000342066.8 Q96NU1:R41Q 0.0866 benign
#> 4 chr1_930245_G_A_hg38 ENST00000342066.8 Q96NU1:D55N 0.134 benign
#> 5 chr1_930248_G_A_hg38 ENST00000342066.8 Q96NU1:G56S 0.100 benign
#> 6 chr1_930282_G_A_hg38 ENST00000342066.8 Q96NU1:R67Q 0.0635 benign
#> 7 chr1_930285_G_A_hg38 ENST00000342066.8 Q96NU1:R68Q 0.0629 benign
#> 8 chr1_930314_C_T_hg38 ENST00000342066.8 Q96NU1:H78Y 0.110 benign
#> 9 chr1_930320_C_T_hg38 ENST00000342066.8 Q96NU1:R80C 0.0918 benign
#> 10 chr1_931058_G_A_hg38 ENST00000342066.8 Q96NU1:V92M 0.196 benign
#> # ℹ more rows
The ClinVar table is now available for exploration or parsing.
Compare ClinVar and AlphaMissense
This section uses the clinvar_plot() function to
generate a scatterplot for benchmarking and comparing ClinVar
classification with AlphaMissense predictions. By default, the function
takes one UniProt accession identifier, derives AlphaMissense scores
from am_data("aa_substitution"), and pulls ClinVar
classifications from the data.frame previously obtained. Alternatively,
it is possible to pass a custom AlphaMissense or ClinVar table to the
function. See function details for more information.
clinvar_plot(uniprotId = "P37023")
#> * [16:30:26][info] 'alphamissense_table' not provided, using default 'am_data("aa_substitution")' table accessed through the AlphaMissenseR package
#> * [16:30:28][info] 'clinvar_table' not provided, using default ClinVar dataset in AlphaMissenseR package
We returned a ggplot object which overlays ClinVar
classifications onto AlphaMissense predicted scores. Blue, gray, and red
colors represent pathogenicity classifications for “likely benign”,
“ambiguous”, or “likely pathogenic”, respectively. Large, bolded points
are ClinVar variants colored according to their clinical classification,
while smaller points in the background are AlphaMissense
predictions.
We can note a discrepancy between the clinically-validated
annotations and the AlphaMissense predictions around position 50.
AlphaMissense seems to predict several variants in that region as likely
benign, while ClinVar identifies them as pathogenic.
Because the ClinVar dataset is not exhaustive (not all proteins have
been clinically-assessed), there may be proteins where information is
not available. In this case, the function will provide an error.
Remember to disconnect from the database.
db_disconnect_all()
#> * [16:30:31][info] disconnecting all registered connections