We estimate the penetrance of LQTS for SCN5A G1631D around
51% and the Brugada syndrome penetrance around
19%. SCN5A G1631D was found in a total of
2 carriers in 3 papers and/or in gnomAD:
0 had Brugada syndrome, 2 had LQTS.
G1631D is not present in gnomAD.
G1631D has been functionally characterized in 4 papers.
This residue is located in a Mild_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 (2 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
G1631D around 51% (3/12)
and the Brugada syndrome penetrance around 19%
(2/12).
In Silico Data
PROVEAN
PolyPhen-2
BLAST-PSSM
REVEL
Penetrance Density BrS (%)
Penetrance Density LQT3 (%)
-6.47
0.954
-4.02
0.982
29
40
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.
G1631D has 58 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.