We estimate the penetrance of LQTS for SCN5A R1626H around
20% and the Brugada syndrome penetrance around
5%. SCN5A R1626H was found in a total of
14 carriers in 4 papers and/or in gnomAD:
0 had Brugada syndrome, 2 had LQTS.
R1626H is present in 12 alleles in gnomAD.
R1626H has been functionally characterized in 4 papers.
This residue is located in a Non_Hotspot region for Brugada syndrome and a Hotspot
region for LQTS.
In silico predictions, functional data (if available), and location in structure are equivalent to phenotyping
10 individuals for Brugada syndrome (1 diagnosed with Brugada syndrome) and 5 individuals for LQTS (1
with LQTS). These data combined with observations of carriers lead us to estimate the LQTS penetrance for SCN5A
R1626H around 20% (3/24)
and the Brugada syndrome penetrance around 5%
(1/24).
In Silico Data
PROVEAN
PolyPhen-2
BLAST-PSSM
REVEL
Penetrance Density BrS (%)
Penetrance Density LQT3 (%)
-4.66
0.998
-3.78
0.986
5
49
PROVEAN scores less than -2 are considered deleterious. REVEL scores higher than 0.5 or 0.75 are considered
likely pathogenic (higher sensitivity with the former cutoff, higher specificity with the latter cutoff).
A PolyPhen-2 score of 0.85 or greater is considered likely pathogenic.
BLAST-PSSM reflects the evolutionary conservation of residue substitutions, more negative numbers indicate
fewer observations of the specific substitution than is expected. Penetrance Density is our previously published method
to calculate the average BrS/LQTS probability density in a shell of residues surrounding a residue of interest
(Kroncke et al. 2019).
Summary totals might not agree with the literature table because of duplicate patients, which were excluded from the
total counts. We do not distinguish here between multiple missense codons. Missense variants are combined across
degenerate codon substitutions since codon-level data were not consistently available for curation.
Functional Data
Peak and late/persistent current are relative to wildtype (100% being no different from wildtype).
V0.5 act/inact are the voltages at which half of the maximal current is reached during an activation and
inactivation protocol, each is in units of mV and relative to wildtype.
R1626H has 71 previously observed neighbors within 15 angstroms
A residue within a folded protein on average has nearest neighbors that fall roughly into two shells: a "nearest"
neighbor around 5-6 angstroms and a second shell around 11 angstroms. All variants shown in the rightmost
column have been observed in at least one individual in the literature or gnomAD.