Next Article in Journal
Research and Development Journey and Future Trends of Hollow Fiber Membranes for Purification Applications (1970–2020): A Bibliometric Analysis
Next Article in Special Issue
Protein Kinase C Regulates Expression and Function of the Cav3.2 T-Type Ca2+ Channel during Maturation of Neonatal Rat Cardiomyocyte
Previous Article in Journal
Human Mesenchymal Stem Cells Seeded on the Natural Membrane to Neurospheres for Cholinergic-like Neurons
Previous Article in Special Issue
Regulation of GABAA Receptors Induced by the Activation of L-Type Voltage-Gated Calcium Channels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

L-Type Calcium Channel: Predicting Pathogenic/Likely Pathogenic Status for Variants of Uncertain Clinical Significance

by
Svetlana I. Tarnovskaya
1,2,*,
Anna A. Kostareva
1,3 and
Boris S. Zhorov
1,4,5,*
1
Almazov National Medical Research Centre, 197341 St. Petersburg, Russia
2
St. Petersburg Academic University, Russian Academy of Sciences, 194021 St. Petersburg, Russia
3
Department of Women’s and Children’s Health and Center for Molecular Medicine, Karolinska Institute, 17177 Stockholm, Sweden
4
Sechenov Institute of Evolutionary Physiology & Biochemistry, Russian Academy of Sciences, 194223 St. Petersburg, Russia
5
Department of Biochemistry and Biomedical Sciences, McMaster University, Hamilton, ON L8S 4K1, Canada
*
Authors to whom correspondence should be addressed.
Membranes 2021, 11(8), 599; https://doi.org/10.3390/membranes11080599
Submission received: 5 July 2021 / Revised: 1 August 2021 / Accepted: 4 August 2021 / Published: 7 August 2021
(This article belongs to the Special Issue Physiology, Pathophysiology and Pharmacology of Calcium Channels)

Abstract

:
(1) Background: Defects in gene CACNA1C, which encodes the pore-forming subunit of the human Cav1.2 channel (hCav1.2), are associated with cardiac disorders such as atrial fibrillation, long QT syndrome, conduction disorders, cardiomyopathies, and congenital heart defects. Clinical manifestations are known only for 12% of CACNA1C missense variants, which are listed in public databases. Bioinformatics approaches can be used to predict the pathogenic/likely pathogenic status for variants of uncertain clinical significance. Choosing a bioinformatics tool and pathogenicity threshold that are optimal for specific protein families increases the reliability of such predictions. (2) Methods and Results: We used databases ClinVar, Humsavar, gnomAD, and Ensembl to compose a dataset of pathogenic/likely pathogenic and benign variants of hCav1.2 and its 20 paralogues: voltage-gated sodium and calcium channels. We further tested the performance of sixteen in silico tools in predicting pathogenic variants. ClinPred demonstrated the best performance, followed by REVEL and MCap. In the subset of 309 uncharacterized variants of hCav1.2, ClinPred predicted the pathogenicity for 188 variants. Among these, 36 variants were also categorized as pathogenic/likely pathogenic in at least one paralogue of hCav1.2. (3) Conclusions: The bioinformatics tool ClinPred and the paralogue annotation method consensually predicted the pathogenic/likely pathogenic status for 36 uncharacterized variants of hCav1.2. An analogous approach can be used to classify missense variants of other calcium channels and novel variants of hCav1.2.

1. Introduction

L-type Cav1.2 calcium channels are expressed in various excitable cells including cardiomyocytes [1]. Defects in gene CACNA1C, which encodes the pore-forming α1 subunit of the hCav1.2 channel, underlie cardiac disorders such as atrial fibrillation, long QT syndrome, conduction disorders, cardiomyopathies, and congenital heart defects [2]. With the advent of whole-exome sequencing data, public databases are rapidly replenished with new gene variants. Over 350 CACNA1C missense variants are listed in public databases but clinical significance is known for only 12% of the variants. The guideline of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) recommends the employment of computational tools to predict damaging variants [3]. Numerous computational tools, which are based on different principles, have been developed to predict the pathogenicity and tolerance of genetic variants [4]. The success rate of these tools varies from 60% to 80% [5]. The ACMG/AMP guideline recommends the employment of multiple software programs for variants’ interpretation because individual programs and underlying algorithms have their own strengths and weaknesses. Therefore, the choice of bioinformatics tools is critical for reliable variant interpretation.
The performance of in silico tools is highly variable and depends on the disease phenotype [6]. For instance, tools MCap, MetaSVM, and MetaLR demonstrated the top performance in predicting the pathogenicity for variants associated with abnormalities in the cardiovascular system [6]. However, some methods yielded many false-positive and false-negative predictions of pathogenicity in individual protein families [7,8,9,10]. For example, MetaSVM, which reportedly has a high accuracy in general [11] and in particular for variants associated with cardiovascular diseases, predicted a deleterious effect for 75% of benign variants of the cardiac sodium channel Nav1.5 [7]. Thus, choosing a tool with a high rate of correct predictions for specific protein families allows to improve predictions [7,12].
Many statistical and machine learning tools, which use only protein sequences, have been proposed to predict variant pathogenicity. An alternative approach is the paralogous annotation method [7,13] that is based on analysis of variants in the multiple sequence alignment of functionally and structurally related proteins. The method assumes that if one protein is known to have a damaging mutation in a certain position, analogous mutation in the sequentially matching position of another protein is likely a damaging one. Consensus predictions of pathogenicity by both sequence-based methods and the paralogous annotation within the family of structurally related voltage-gated sodium and calcium channels is expected to provide more reliable pathogenicity prediction for hCav1.2 than individual approaches alone.
A comprehensive predictor assessment requires a benchmark with both positive (pathogenic/likely pathogenic) and negative (benign) variants. Here, we composed a dataset to test the performance of various predictors in identifying known pathogenic/likely pathogenic (P/LP) and benign variants in the families of voltage-gated sodium and calcium channels. We collected common (benign) missense variants from the gnomAD database and P/LP missense variants from three databases that are referred to in the next section. We evaluated the performance of 16 popular prediction tools and identified top-performing tools for the hCav1.2 channel and its paralogs. The best-performing bioinformatics tool, ClinPred, and the paralogue annotation method consensually predicted that 36 hCav1.2 variants of uncertain clinical significance are putative pathogenic/likely pathogenic variants.

2. Methods

2.1. Data Collection and Preprocessing

Paralogues of the hCav1.2 channel were previously identified [7]. Sequences of hCav1.2 and its paralogues were obtained from the UniProt database [14]. Pathogenic/likely pathogenic variants of hCav1.2 and its paralogues were collected from three databases: Humsavar (https://www.uniprot.org/docs/humsavar, last visited 26 January 2021), Ensembl Variation [15], and ClinVar [16]. Only ‘disease’ variants were extracted from the databases Ensembl and Humsavar. From the ClinVar database, we selected variants that are characterized as ‘pathogenic’ or ‘likely pathogenic’ and are associated with specific clinical conditions. Benign (neutral) variants along with their minor allele frequencies (AF) were obtained from the population database gnomAD (data release 2.1.1, October 2018) [17]. Variants with AF > 0.0001 were considered as benign [17,18]. Variants of uncertain clinical significance were extracted from ClinVar and Ensembl. All variants were combined into one broad dataset (Table S1).

2.2. hCav1.2 Topology

Region borders of the hCav1.2 channel were determined according to Uniprot entry Q13936. The pore-forming α1 subunit of the Cav1.2 channel folds from a single polypeptide chain of four homologous repeat domains (DI–DIV), which are connected by intracellular linkers. Each domain has six transmembrane helices (S1–S6) and a large extracellular membrane-reentering P-loop (Figure 1). In each repeat, helices S1–S4 constitute a voltage sensing domain, whereas helices S5, S6, and the P-loop contribute a quarter to the pore domain.

2.3. Multiple Sequence Alignment and Paralogue Annotation

The paralogue annotation method specifically identifies P/LP missense variants by transferring annotations across families of related proteins [13,19]. Recently a modified method of paralogue annotation [13] was used to predict the pathogenicity of variants in the cardiac sodium channel hNav1.5 [7]. This approach is applied here to select potential P/LP variants from the large set of VUS in the Cav1.2 channel.
We have chosen 20 paralogues of hCav1.2: ten voltage-gated sodium channels, nine voltage-gated calcium channels, and a voltage-independent non-selective sodium leak channel (Table 1). For each paralogue, P/LP variants were collected (Section 2.1). Amino acid sequences of Cav1.2 and its paralogues were aligned by Clustal Omega [20]. Each paralogue mutation was mapped on the Cav1.2 sequence according to the alignment (Table S2). Position-specific conservation scores (Cs) were calculated using the Zvelebil method [21] as implemented in the Amino Acid Conservation Calculation Service [22]. Variants in positions with conservation scores > 0.3 were considered as P/LP variants according to References [7,19].

2.4. Annotation of Missense Variants

We used the dbNSFPv4 database [23] to evaluate the performance of popular tools in predicting P/LP and benign variants. The database contains pre-computed scores for all potential amino acid substitutions that are taken from the following tools: SIFT, Polyphen HVAR, Polyphen HDIV, Mutation Taster, Mutation Assessor, PROVEAN, FATHMM, FATHMM_XF, MetaSVM, MetaLR, CADD, ClinPred, REVEL, PrimateAI, Eigen, and MCap [4]. The performance of algorithms was evaluated using the broad dataset, which includes P/LP and benign variants from hCav1.2 and its 20 paralogues. Variant prioritization tools do not always provide scores for every variant reported in dbNSFP. We selected only those tools in cases where scores are missing for less than 30% of the variants in our dataset. Altogether, 16 different scores were considered. We calculated the area under the ROC (receiver operating characteristic) curve (AUC) using the library pROC in R [24]. The higher the AUC score, the better is the tool performance in the dataset. ROC curves were obtained by plotting sensitivity against (1–specificity) at each threshold for each tool.
To increase the reliability of the analysis, we determined the optimal pathogenicity threshold using the R package pROC (Table 2). All variants were divided into two categories (P/LP or benign) according to the calculated thresholds. To evaluate the performance of these tools, we calculated the following standard characteristics:
Sensitivity = T P T P + F N
Specificity = T N T N + F P
Accuracy = T P + T N T P + F P + T N + F N
where TP (true positive) is the number of P/LP variants correctly predicted as pathogenic; FN (false negative) is the number of P/LP variants incorrectly predicted as benign; TN (true negative) is the number of benign variants correctly predicted as benign; and FP (false positive) is the number of benign variants incorrectly predicted as pathogenic.
Custom threshold is the custom pathogenicity threshold that divides variants in two categories: pathogenic or benign. The larger or smaller the score than the threshold, the more likely the variant is damaging. Sensitivity characterizes the number of P/LP variants, which were predicted as P/LP by the tool, while specificity characterizes the number of benign variants, which were predicted as benign by the tool. Accuracy indicates the predictive accuracy of the tool [25].

3. Results

3.1. Composing a Broad Dataset of Missense Variants for Channel hCav1.2 and Its Paralogues

For the 21 channels listed in Table 1, we collected a total of 7164 missense variants from the databases gnomAD, ClinVar, Uniprot, and Humsavar (Table S1). These include 1549 P/LP variants, 763 benign variants (with AF > 0.0001), and 4852 uncharacterized variants or VUS. We further refer to this dataset as the “broad dataset”. The largest numbers of P/LP variants were found for channels hNav1.1, hNav1.5, and hNav1.2 (605, 350, and 166, respectively). No P/LP variants were found for channels hCav3.3 and hNav2.1 (Table 1). For channel hCav1.2, we found 22 P/LP variants, 21 benign variants (with AF > 0.0001), and 309 VUS (Table 1 and Table S1).

3.2. Distribution of Missense Variants in Topological Regions of hCav1.2

To identify hCav1.2 regions with pathogenic or benign variants, we explored the occurrence of mutation in the regions. About 70% of hCav1.2 missense variants are localized in cytoplasmic linkers DI-DII and DII-DIII, and in N and C-terminal parts. Most of the P/LP variants appear in linkers DI-DII and DII-DII. Most of the benign variants are localized in the C-terminal part of the channel protein (Figure 1). Highly conserved transmembrane domains contain relatively few P/LP variants.
Most of the P/LP variants of hCav1.2 are associated with long QT syndrome (56%). Other variants are associated with Timothy syndrome, Brugada syndrome, inborn genetic diseases, and other CACNA1C-related disorders (Table S1). Variants causing long QT syndrome and Timothy syndrome most often are localized in interdomain linkers DI-DII and DII-DIII. Some variants are associated with two or more diseases.

3.3. Amino Acid Substitutions in P/LP and Benign Variants

We analyzed statistics of amino acid substitutions in the 21 channels. Arg and Leu have high mutation rates in the P/LP variants, whereas Arg and Ala are highly mutable in benign variants. The most frequent event, irrespective of disease relevance, is the substitution of a hydrophobic residue by another hydrophobic residue. Hydrophobic substitutions of Leu by Pro are frequently associated with diseases, whereas in neutral variants, substitutions of Pro by Leu are most frequent (Table 3).

3.4. Paralogue Annotation of Variants Identified in hCav1.2

The human genome has 20 paralogues of the Cav1.2 channel: ten voltage-gated sodium channels, nine voltage-gated calcium channels, and a voltage-independent non-selective sodium leak channel (Table 1). Using the multiple sequence alignment of hCav1.2 and its paralogues (Section 2.3), we mapped each paralogue protein residue with a known P/LP variant onto the sequentially matching amino acid of hCav1.2 (Table S2).
A total of 146 known P/LP variants in paralogues were mapped to 89 variants in Cav1.2 (Table S2). Among these, 11 variants correspond to known P/LP variants of Cav1.2, two correspond to benign variants, and 76 correspond to VUS. Most of the paralogue variants were mapped to linkers DI-DII and DII-DIII, N and C-terminal parts, and S5/S6 loops of repeat domains DI and DIII (Table S2).

3.5. Comparing Performance of the Computational Tools

We compared the performance of 16 in-silico tools that predict the pathogenicity of variants (Table 2) [23]. Pre-computed algorithm scores were retrieved from the dbNSFPv4 database. ROC curves and AUC are shown in Figure 2. For each tool, we determined the optimal pathogenicity threshold and calculated sensitivity, specificity, and accuracy (Table 2).
We used AUC as the main measure of performance. ClinPred demonstrated the best predictive performance (AUC = 0.97) in the broad dataset (Figure 2), followed by REVEL (AUC = 0.93) and MCap (AUC = 0.91). Moreover, ClinPred is the only tool that performed with accuracy as (AUC) > 0.90. It correctly classified 95% of the P/LP variants as pathogenic variants and 94% of the benign variants as benign variants (Table 2). MutationAssessor, Polyphen, FATHMM, and MutationTaster performed with an accuracy of < 0.80. The accuracy of other algorithms ranged from 0.80 to 0.88. The lowest accuracy across all methods was found for MutationTaster (accuracy = 0.62, AUC = 0.66). For the hCav1.2 channel, ClinPred, CADD, MetaSVM, Polyphen HDIV, and MutationTaster correctly classified all 22 P/LP variants as pathogenic. However, these methods also assigned pathogenic status for some benign variants. ClinPred predicted 4/21 (19%) benign variants as pathogenic. Polyphen HDIV and MutationTaster classified less than 60% of benign variants as benign.
The results indicate that ClinPred is the best-performing pathogenicity predictor for variants in the family of voltage-gated sodium and calcium channels.

3.6. Reclassifying VUS Variants of hCav1.2 with ClinPred and Paralogue Annotation

Most of the Cav1.2 variants in our broad dataset are currently classified as VUS. We used the best-performing tool, ClinPred, to predict the pathogenicity of VUS’. ClinPred predicted 188 VUS variants as pathogenic/likely pathogenic (pathogenicity threshold > 0.66). Among the 188 VUS, we further selected only those variants that are annotated as pathogenic or likely pathogenic in at least one of 20 paralogs of Cav1.2 (conservation score across paralogues Cs > 0.3). Both methods consensually predicted 36 of 309 VUS as P/LP variants. We reclassified these as putative P/LP variants (Table 4). The variants were localized mainly in the extracellular loops DI-S5/S6 and DIII-S5/S6 (Figure 1).

4. Discussion

Numerous bioinformatics methods are used to predict the deleteriousness of missense variants. These methods rely on supervised machine learning models that are trained on a collection of manually annotated variants to predict the probable pathogenicity for each amino acid substitution. Our recent analysis revealed that in the cardiac sodium channel hNav1.5, the MetaSVM method, which reportedly has a high accuracy in general [11] and in particular for variants associated with cardiovascular diseases [6], predicted a deleterious effect for 75% of variants that are annotated as benign [7]. It is recommended to find a specific predictor with the optimal threshold value for each family of proteins [12].
In the present study, we have shown that various popular bioinformatics tools yield different predictions of pathogenicity for known P/LP and benign variants in the family of the hCav1.2 channel and its paralogues. In our broad dataset, ClinPred performed best with an area of 0.97 under the receiver operating characteristic curve. ClinPred is a meta-predictor that combines commonly used and recently developed individual prediction tool scores, as well as the allele frequency of variants in different populations from the gnomAD database [26]. Its high AUC and prediction accuracy values in our dataset suggest that ClinPred was trained on large newer datasets, which overlap with the dataset used in this study.

5. Conclusions

Here, we created a broad dataset that includes all known missense variants in the hCav1.2 channel and 20 paralogues voltage-gated calcium and sodium channels. Most of the known pathogenic/likely pathogenic variants are found in intracellular linkers DI-DII and DII-DIII. Among sixteen pathogenicity prediction tools, which were tested using our broad dataset, ClinPred demonstrated the best performance in distinguishing pathogenic/likely pathogenic variants from benign ones. The best-performing tool is expected to improve in silico assessment of clinically relevant variants of hCav1.2 and its paralogues. ClinPred and the paralogue annotation method consensually predicted that 36 variants of hCav1.2, which are currently classified as variants of unknown significance, are pathogenic/likely pathogenic variants. Most of these variants are located in the extracellular loops DI-S5/S6 and DIII-S5/S6. The reclassified variants can be used for diagnostics of cardiac diseases. They are also promising targets for further experimental and theoretical studies.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/membranes11080599/s1. Table S1: Broad dataset of Cav1.2 variants and its paralogues, and Table S2: Mapping of P/LP paralogue variants on the human Cav1.2 sequence.

Author Contributions

Conceptualization, B.S.Z.; bioinformatics research, S.I.T.; data curation, B.S.Z. and A.A.K.; writing—original draft preparation, S.I.T.; writing—review and editing, A.A.K. and B.S.Z.; project administration, A.A.K.; funding acquisition, B.S.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by grants to BSZ from the Russian Science Foundation (17-15-01292-П) and Natural Sciences and Engineering Research Council of Canada (RGPIN-2020-07100).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AF, allele frequencies; AUC, area under the ROC curve; DI-DIV, repeat domains in eukaryotic sodium and calcium channels; P/LP, pathogenic/likely pathogenic variants; ROC, receiver operating characteristic; S1-S6, transmembrane segments; S1/S2, S3/S4, and S5/S6, extracellular loops between transmembrane segments; S2/S3, cytoplasmic loops between transmembrane segments; S4/S5, linker helices between voltage-sensing domains and the pore domain; and VUS, variant of uncertain significance.

References

  1. Zamponi, G.W.; Striessnig, J.; Koschak, A.; Dolphin, A.C. The physiology, pathology, and pharmacology of voltage-gated calcium channels and their future therapeutic potential. Pharmacol. Rev. 2015, 67, 821–870. [Google Scholar] [CrossRef] [Green Version]
  2. Benitah, J.-P.; Alvarez, J.L.; Gómez, A.M. L Type Ca2+ current in ventricular cardiomyocytes. J. Mol. Cell. Cardiol. 2010, 48, 26–36. [Google Scholar] [CrossRef]
  3. Richards, S.; Aziz, N.; Bale, S.; Bick, D.; Das, S.; Gastier-Foster, J.; Grody, W.W.; Hegde, M.; Lyon, E.; Spector, E.; et al. Standards and guidelines for the interpretation of sequence variants: A joint consensus recommendation of the American college of medical genetics and genomics and the association for molecular pathology. Genet. Med. 2015, 17, 405–423. [Google Scholar] [CrossRef]
  4. Dong, C.; Wei, P.; Jian, X.; Gibbs, R.; Boerwinkle, E.; Wang, K.; Liu, X. Comparison and integration of deleteriousness prediction methods for nonsynonymous SNVs in whole exome sequencing studies. Hum. Mol. Genet. 2015, 24, 2125–2137. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Ghosh, R.; Oak, N.; Plon, S.E. Evaluation of in silico algorithms for use with ACMG/AMP clinical variant interpretation guidelines. Genome Biol. 2017, 18, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Anderson, D.; Lassmann, T. A phenotype centric benchmark of variant prioritisation tools. NPJ Genom. Med. 2018, 3, 5. [Google Scholar] [CrossRef] [Green Version]
  7. Tarnovskaya, S.I.; Korkosh, V.S.; Zhorov, B.S.; Frishman, D. Predicting novel disease mutations in the cardiac sodium channel. Biochem. Biophys. Res. Commun. 2020, 521, 603–611. [Google Scholar] [CrossRef] [PubMed]
  8. Sun, H.; Yu, G. New insights into the pathogenicity of non-synonymous variants through multi-level analysis. Sci. Rep. 2019, 9, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Niroula, A.; Vihinen, M. How good are pathogenicity predictors in detecting benign variants? PLoS Comput. Biol. 2019, 15, 1–17. [Google Scholar] [CrossRef] [Green Version]
  10. Miosge, L.A.; Field, M.A.; Sontani, Y.; Cho, V.; Johnson, S.; Palkova, A.; Balakishnan, B.; Liang, R.; Zhang, Y.; Lyon, S.; et al. Comparison of predicted and actual consequences of missense mutations. Proc. Natl. Acad. Sci. USA 2015, 112, E5189–E5198. [Google Scholar] [CrossRef] [Green Version]
  11. Kim, S.; Jhong, J.-H.; Lee, J.; Koo, J.-Y. Meta-analytic support vector machine for integrating multiple omics data. BioData Min. 2017, 10, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Zaucha, J.; Heinzinger, M.; Tarnovskaya, S.; Rost, B.; Frishman, D. Family-specific analysis of variant pathogenicity prediction tools. NAR Genom. Bioinform. 2020, 2, 1–8. [Google Scholar] [CrossRef] [Green Version]
  13. Ware, J.S.; Walsh, R.; Cunningham, F.; Birney, E.; Cook, S.A. Paralogous annotation of disease-causing variants in long qt syndrome genes. Hum. Mutat. 2012, 33, 1188–1191. [Google Scholar] [CrossRef] [Green Version]
  14. The UniProt Consortium. UniProt: A hub for protein information. Nucleic Acids Res. 2014, 43, D204–D212. [Google Scholar] [CrossRef]
  15. Chen, Y.; Cunningham, F.; Rios, D.; McLaren, W.M.; Smith, J.; Pritchard, B.; Spudich, G.M.; Brent, S.; Kulesha, E.; Marin-Garcia, P.; et al. Ensembl variation resources. BMC Genom. 2010, 11, 293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Landrum, M.J.; Lee, J.M.; Benson, M.; Brown, G.; Chao, C.; Chitipiralla, S.; Gu, B.; Hart, J.; Hoffman, D.; Hoover, J.; et al. ClinVar: Public archive of interpretations of clinically relevant variants. Nucleic Acids Res. 2016, 44, D862–D868. [Google Scholar] [CrossRef] [Green Version]
  17. Kaltman, J.R.; Evans, F.; Fu, Y.-P. Re-evaluating pathogenicity of variants associated with the long QT syndrome. J. Cardiovasc. Electrophysiol. 2018, 29, 98–104. [Google Scholar] [CrossRef] [PubMed]
  18. Walsh, R.; Thomson, K.L.; Ware, J.S.; Funke, B.H.; Woodley, J.; McGuire, K.J.; Mazzarotto, F.; Blair, E.; Seller, A.; Taylor, J.C.; et al. Reassessment of mendelian gene pathogenicity using 7,855 cardiomyopathy cases and 60,706 reference samples. Genet. Med. 2017, 19, 192–203. [Google Scholar] [CrossRef] [Green Version]
  19. Walsh, R.; Peters, N.S.; Cook, S.A.; Ware, J.S. Paralogue annotation identifies novel pathogenic variants in patients with brugada syndrome and catecholaminergic polymorphic ventricular tachycardia. J. Med. Genet. 2014, 51, 35–44. [Google Scholar] [CrossRef] [Green Version]
  20. Sievers, F.; Wilm, A.; Dineen, D.; Gibson, T.J.; Karplus, K.; Li, W.; Lopez, R.; McWilliam, H.; Remmert, M.; Söding, J.; et al. Fast, scalable generation of high-quality protein multiple sequence alignments using clustal omega. Mol. Syst. Biol. 2011, 7, 539. [Google Scholar] [CrossRef]
  21. Zvelebil, M.J.; Barton, G.J.; Taylor, W.R.; Sternberg, M.J.E. Prediction of protein secondary structure and active sites using the alignment of homologous sequences. J. Mol. Biol. 1987, 195, 957–961. [Google Scholar] [CrossRef]
  22. Manning, J.R.; Jefferson, E.R.; Barton, G.J. The Contrasting Properties of Conservation and Correlated Phylogeny in Protein Functional Residue Prediction. BMC Bioinforma. 2008, 9, 51. [Google Scholar] [CrossRef] [Green Version]
  23. Liu, X.; Li, C.; Mou, C.; Dong, Y.; Tu, Y. DbNSFP v4: A comprehensive database of transcript-specific functional predictions and annotations for human nonsynonymous and splice-site SNVs. Genome Med. 2020, 12, 1–8. [Google Scholar] [CrossRef]
  24. Robin, X.; Turck, N.; Hainard, A.; Tiberti, N.; Lisacek, F.; Sanchez, J.-C.; Müller, M. PROC: An open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinform. 2011, 12, 77. [Google Scholar] [CrossRef]
  25. Niroula, A.; Vihinen, M. Variation interpretation predictors: Principles, types, performance, and choice. Hum. Mutat. 2016, 37, 579–597. [Google Scholar] [CrossRef] [PubMed]
  26. Alirezaie, N.; Kernohan, K.D.; Hartley, T.; Majewski, J.; Hocking, T.D. ClinPred: Prediction tool to identify disease-relevant nonsynonymous single-nucleotide variants. Am. J. Hum. Genet. 2018, 103, 474–483. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Topology of the hCav1.2 channel. Four repeat domains (DI–DIV) are connected by intracellular linkers. Each domain contains six transmembrane segments (S1–S6). Variants of uncertain clinical significance, which are reclassified here as P/LP variants, are marked by green triangles; known P/LP variants are marked by red triangles; and benign variants are marked by blue triangles.
Figure 1. Topology of the hCav1.2 channel. Four repeat domains (DI–DIV) are connected by intracellular linkers. Each domain contains six transmembrane segments (S1–S6). Variants of uncertain clinical significance, which are reclassified here as P/LP variants, are marked by green triangles; known P/LP variants are marked by red triangles; and benign variants are marked by blue triangles.
Membranes 11 00599 g001
Figure 2. The ROC curves of the 16 tools tested on the broad dataset. The higher AUC score (shown in brackets), the better the performance.
Figure 2. The ROC curves of the 16 tools tested on the broad dataset. The higher AUC score (shown in brackets), the better the performance.
Membranes 11 00599 g002
Table 1. Number of missense variants in the human Cav1.2 and its paralogs.
Table 1. Number of missense variants in the human Cav1.2 and its paralogs.
GeneChannelBenignP/LPVUS
CACNA1AhCav2.1560346
CACNA1BhCav2.248010
CACNA1ChCav1.22122309
CACNA1DhCav1.34746
CACNA1EhCav2.3361718
CACNA1FhCav1.4392839
CACNA1GhCav3.142545
CACNA1HhCav3.211621454
CACNA1IhCav3.36101
CACNA1ShCav1.15711290
NALCNhNavi2.1183117
SCN1AhNav1.118605482
SCN2AhNav1.215166301
SCN3AhNav1.3168202
SCN4AhNav1.44778312
SCN5AhNav1.543350705
SCN7AhNav2.16500
SCN8AhNav1.6693249
SCN9AhNav1.71132505
SCN10AhNav1.8663291
SCN11AhNav1.92912230
Table 2. Performance of variant interpretation tools.
Table 2. Performance of variant interpretation tools.
ToolCustom ThresholdSensitivitySpecificityAccuracy
ClinPred>0.660.950.940.95
REVEL>0.730.890.870.88
MetaLR>0.910.860.840.85
Eigen>0.550.810.870.84
MCap>0.590.820.870.84
PrimateAI>0.690.880.790.84
MetaSVM>0.920.880.850.87
FATHMM_XF>0.850.820.840.83
SIFT<0.00250.830.830.83
PROVEAN<−3.420.800.800.80
CADD>240.880.720.80
MutationAssessor>2.510.770.810.79
Polyphen HVAR>0.60.800.700.75
FATHMM<−4.3050.760.730.74
Polyphen HDIV>0.960.740.710.73
MutationTaster>0.990.940.310.62
Table 3. Occurrence of top five residue types in P/LP and benign variants.
Table 3. Occurrence of top five residue types in P/LP and benign variants.
Type of DataWT ResidueMutant ResiduePreferred Substitutions
All variants
(n = 2312)
R, A, V, L, GV, T, S, R, IR > Q, L > P, R > H, A > T, R > C
P/LP variants
(n = 1549)
R, L, V, G, AM, R, C, D, HL > P, G > R, R > Q, R > H, R > C
Benign variants
(n = 763)
R, A, P, V, GV, S, T, A, RA > T, R > Q, P > L, R > H, V > I
Table 4. hCav1.2 variants whose status was changed from “uncertain significance” to “pathogenic/likely pathogenic”.
Table 4. hCav1.2 variants whose status was changed from “uncertain significance” to “pathogenic/likely pathogenic”.
#VariantLocationParalogue
1A180TDI-S2/S3SCN5A-A178G, SCN1A-A175V, SCN1A-A175T
2T330MDI-S5/S6SCN1A-Y349C
3G342SDI-S5/S6SCN5A-G351V, SCN5A-G351D
4T344IDI-S5/S6SCN1A-T363R, SCN1A-T363P, SCN5A-T353I
5T366MDI-S5/S6SCN1A-E385Q
6G377ADI-S5/S6SCN5A-G386R, SCN5A-G386E
7G402RDI-DIICACNA1E-G348R, SCN1A-A420V, CACNA1D-G403D, CACNA1F-G369D
8W528GDII-S1SCN1A-L772P
9A562TDII-S2CACNA1H-V831M
10A565TDII-S2SCN5A-G758E
11I630VDII-S4SCN1A-L869S, SCN1A-L869F
12L644WDII-S4/S5SCN9A-I859T, SCN3A-I875T, SCN4A-I693T, SCN8A-I868T
13I651NDII-S4/S5SCN5A-L839P, CACNA1E-I603L, SCN8A-L875Q, CACNA1A-I614M, SCN1A-L890P
14I651VDII-S4/S5SCN5A-L839P, CACNA1E-I603L, SCN8A-L875Q, CACNA1A-I614M, SCN1A-L890P
15L658PDII-S5NALCN-T513N, SCN1A-L897S, SCN1A-L897F
16F692IDII-S5/S6SCN2A-F928C
17G705RDII-S5/S6SCN1A-G950E, SCN1A-G950R
18I751VDII-DIIICACNA1D-I750F, CACNA1A-I712V, CACNA1E-I701V, CACNA1F-I756T, CACNA1D-I750M
19A757PDII-DIIISCN2A-S987I
20E760KDII-DIIISCN1A-D998G
21L912IDIII-S1SCN1A-L1230F
22I932TDIII-S2SCN1A-M1251R, SCN5A-L1238P
23F936CDIII-S2SCN1A-A1255D, SCN1A-A1255P
24K1074RDIII-S5/S6SCN5A-K1359N
25D1119NDIII-S5/S6SCN1A-D1416G
26A1124VDIII-S5/S6SCN1A-G1421R, SCN1A-G1421E, SCN5A-G1408R
27V1131IDIII-S5/S6CACNA1A-V1456L, SCN1A-V1428A, SCN1A-V1428F
28E1135GDIII-S5/S6SCN1A-K1432I, SCN5A-K1419E, SCN2A-K1422E
29E1135KDIII-S5/S6SCN1A-K1432I, SCN5A-K1419E, SCN2A-K1422E
30G1136ADIII-S5/S6SCN1A-G1433V, SCN5A-G1420R, SCN5A-G1420V, SCN1A-G1433R, SCN1A-G1433E
31G1153DDIII-S5/S6SCN1A-Q1450K, SCN1A-Q1450R
32A1174SDIII-S6SCN5A-S1458Y, SCN1A-S1471F
33F1246LDIV-S1SCN4A-M1360V, SCN2A-M1538I
34I1373VDIV-S4SCN1A-P1632S
35I1421VDIV-S5SCN1A-I1683F, SCN4A-I1495F, SCN1A-I1683T
36A1617TC-termSCN5A-V1861I
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tarnovskaya, S.I.; Kostareva, A.A.; Zhorov, B.S. L-Type Calcium Channel: Predicting Pathogenic/Likely Pathogenic Status for Variants of Uncertain Clinical Significance. Membranes 2021, 11, 599. https://doi.org/10.3390/membranes11080599

AMA Style

Tarnovskaya SI, Kostareva AA, Zhorov BS. L-Type Calcium Channel: Predicting Pathogenic/Likely Pathogenic Status for Variants of Uncertain Clinical Significance. Membranes. 2021; 11(8):599. https://doi.org/10.3390/membranes11080599

Chicago/Turabian Style

Tarnovskaya, Svetlana I., Anna A. Kostareva, and Boris S. Zhorov. 2021. "L-Type Calcium Channel: Predicting Pathogenic/Likely Pathogenic Status for Variants of Uncertain Clinical Significance" Membranes 11, no. 8: 599. https://doi.org/10.3390/membranes11080599

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop