- Research article
- Open Access
Evolutionary analyses of KCNQ1 and HERG voltage-gated potassium channel sequences reveal location-specific susceptibility and augmented chemical severities of arrhythmogenic mutations
© Jackson and Accili; licensee BioMed Central Ltd. 2008
- Received: 15 April 2008
- Accepted: 30 June 2008
- Published: 30 June 2008
Mutations in HERG and KCNQ1 potassium channels have been associated with Long QT syndrome and atrial fibrillation, and more recently with sudden infant death syndrome and sudden unexplained death. In other proteins, disease-associated amino acid mutations have been analyzed according to the chemical severity of the changes and the locations of the altered amino acids according to their conservation over metazoan evolution. Here, we present the first such analysis of arrhythmia-associated mutations (AAMs) in the HERG and KCNQ1 potassium channels.
Using evolutionary analyses, AAMs in HERG and KCNQ1 were preferentially found at evolutionarily conserved sites and unevenly distributed among functionally conserved domains. Non-synonymous single nucleotide polymorphisms (nsSNPs) are under-represented at evolutionarily conserved sites in HERG, but distribute randomly in KCNQ1. AAMs are chemically more severe, according to Grantham's Scale, than changes observed in evolution and their severity correlates with the expected chemical severity of the involved codon. Expected chemical severity of a given amino acid also correlates with its relative contribution to arrhythmias. At evolutionarily variable sites, the chemical severity of the changes is also correlated with the expected chemical severity of the involved codon.
Unlike nsSNPs, AAMs preferentially locate to evolutionarily conserved, and functionally important, sites and regions within HERG and KCNQ1, and are chemically more severe than changes which occur in evolution. Expected chemical severity may contribute to the overrepresentation of certain residues in AAMs, as well as to evolutionary change.
- Amino Acid Change
- Pore Region
- Cystic Fibrosis Transmembrane Regulator
- Disease Mutation
- Voltage Sense Domain
Two voltage-gated potassium ion channel genes, KCNQ1 and KCNH2 (HERG), encode for channels that underlie the slowly- and rapidly-activated delayed rectifier potassium currents (IKs and IKr), respectively [1–4]. Efflux of potassium ions through these channels is critical for repolarization of the cardiac action potential. Mutations that disrupt normal biosynthesis and function of KCNQ1 and HERG have been associated with three cardiac arrhythmias: Short QT syndrome (SQTS) [5–7], atrial fibrillation [8, 9] and Long QT syndrome (LQTS) [10, 11]. To date, there are close to 200 reported arrhythmia-associated mutations (AAMs) in each channel, with more than 95% of which are linked to LQTS. This arrhythmia, which affects an estimated 1 in 5000–10000 people worldwide, is characterized by a prolongation of the QT interval on an electrocardiogram and can lead to Torsade de Pointes (TdP), ventricular fibrillation and sudden cardiac death [12, 13]. More recently, the progression from LQTS into TdP has been proposed as a cause of sudden infant death syndrome (SIDS) [14–16] and sudden unexplained death syndrome (SUDS) . Ninety percent of all known LQTS-associated mutations occur in HERG and KCNQ1.
In studies of larger groups of proteins, or individual proteins other than HERG and KCNQ1, disease-associated amino acid mutations (DAMs) have been analyzed according to the chemical severity of the change, as determined from Grantham's Scale , and the location and/or context of the altered amino acid . An up to two-fold increase in clinically observable disease occurs in parallel with increases in the amount of chemical change . In many proteins, DAMs are chemically more severe than changes over the course of evolution (called interspecific changes) and polymorphic changes [21, 22]. In rhodopsin, the chemical severity of DAMs also correlates with the expected chemical severity of a given codon , as determined by comparison with the normally occurring human codon.
The importance of an amino acid to protein function can be inferred from its conservation over the course of metazoan evolution. In many proteins, DAMs are overabundant at evolutionarily conserved and slowly evolving sites [21–24], presumably because sites that have experienced little interspecific variation are critical for function. Mutations at these sites are likely deleterious and would be removed from the population by natural selection, if given enough time. In some proteins, DAMs are unevenly distributed among functionally conserved domains , even after accounting for the length and number of evolutionarily conserved amino acids. These data imply that functionally conserved domains, like conserved sites, are less tolerant to mutations because of their greater importance to the overall protein function. However, according to the Neutral Theory of Molecular Evolution, most nucleotide substitutions are phenotypically neutral and avoid natural selection . Most nucleotide changes, except for disease-associated changes, might then be expected to distribute randomly throughout a protein. In a study of a large number of proteins, synonymous single nucleotide polymorphisms (sSNPs) were shown to distribute randomly consistent with a neutral phenotype . Initial studies of non-synonymous (ns) SNPs in the cystic fibrosis transmembrane regulator and tuberous sclerosis complex 2 gene, were also shown to distribute randomly . However, a more recent, study, using many proteins, showed that nsSNPs preferentially locate at variable sites and sites with high evolutionary rates, and are underrepresented at sites that are evolutionarily conserved and have low evolutionary rates . This discrepancy may be due to differences in the nature of the nsSNPs available e.g. ns SNPs are less disruptive when located at conserved sites in some proteins versus others, or to a difference in the nature of the conserved sites e.g. conserved sites in some proteins are less tolerant to amino acid changes.
In this study, we take advantage of the > 200 mutations reported per channel to quantitatively analyze, for the first time, the distribution of AAMs and non-synonymous SNPs within HERG and KCNQ1 and to determine the chemical severity of these changes, as compared to interspecific changes.
Human protein and mRNA sequences for HERG and KCNQ1 were obtained from NCBI. Other vertebrate sequences were collected from a BLASTP  search at ENSEMBL (default parameters) and the NCBI non-redundant database. To eliminate the inclusion of sequences from closely related paralogs, a manual reciprocal best-hit analysis was used. Sequences were aligned by ClustalX  using default settings. Less conserved regions in the N and/or C-termini were removed, leaving sequences corresponding to residues 1–906 and 83–599 of human HERG and human KCNQ1, respectively.
Collection of arrhythmia-associated mutations (AAMs) and non-synonymous polymorphisms (nsSNPs)
Breakdown of clinical phenotype of disease mutations included in analyses
Phylogenetic analysis and determination of interspecific variability
Aligned protein sequences, with gaps removed, were used to construct a phylogenetic tree using the maximum likelihood analysis in PHYLIP . SEQBOOT, PROML, and CONSENSE programs were used with a JTT model of evolution. Ancestral sequences were determined using the maximum likelihood method from PAMLv3.15  under a Poisson model of amino acid evolution. Discrete values for interspecific variability (0 through n) were determined for each residue in the protein from differences between the ancestral and descendent sequences throughout the provided tree.
Association between AAMs or nsSNPs and evolutionarily conserved sites
To determine whether AAMs and nsSNPs are preferentially associated with evolutionarily conserved sites, we compared the number of observed mutations in HERG and KCNQ1 to the number of mutations expected at each site in both proteins based on neutral substitution . Sites were binned according to counts of interspecific variability (0 through i) which were determined using PAML. The expected number of mutations was determined using the following equationDiexpected = Dtotal × Ni/N (1)
where Diexpected is the expected number of mutations at sites that have undergone i substitutions, Dtotal is the total number of disease mutations for each channel given by ΣDiobserved, Ni is the observed number of sites in the alignment that have undergone i substitutions and N is the total number of residues in the gene being examined. Therefore, Ni/N is the fraction of sites in the gene that belong to a particular variability class (i), and if disease mutations distribute randomly throughout the protein, then Diexpected will be proportional to the fraction of the total sites and the total number of disease mutations observed in each gene.
To determine whether differences in the distribution pattern between observed values and those expected from neutral theory were significant, the X2 statistic was calculated  and compared to a critical value for the given degrees of freedom (i-1) using the following equation,X2 = Σ0i (Dobserved-Dexpected)2/Dexpected(2)
Determination of codon evolutionary rate of change
The evolutionary rates of change for codons were estimated using the maximum likelihood method implemented in the CODEML program of PAML [34, 23], using a discrete gamma model (eight categories). The shape parameter was either fixed or free to vary and a likelihood ratio test was performed to evaluate model fitting. Evolutionary rates based on a Poisson model of evolution were established for every site and normalized to the maximum rate observed for each protein. Values between 0 and 1 were binned into eight different categories and used to represent eight different levels of evolutionary change. The analysis was performed on nucleotide sequence alignments of core regions, with gaps removed, of human and three closely related vertebrate orthologs, guided by protein sequence alignments. The expected number of mutations at the codons belonging to each rate category was calculated using a modification of Eq. 1 where Ni is the number of codons belonging to i category, N is the total number of amino acid positions, and Dtotal is the total number of disease mutations for each protein used in the analysis.
Distribution of AAMs among functionally important regions of the channels
To quantify over- or underrepresentation of AAMs in functionally conserved regions of the channels, we compared their distribution to that expected from a uniform or evolutionary hypothesis . First, we tested whether AAMs were distributed uniformly across the protein. The expected number of mutations in a given region was determined using the following equationDjexpected = (Rj/R) × Dtotal(3)
where Djexpected is the expected number of disease mutations in a particular region, j, Rj is the number of residues found in region j, R is the total number of residues used in the analysis (ΣRj) and Dtotal = ΣDjobserved or the total number of disease mutations used in the analysis.
The X2 statistic was calculated and compared to a critical value for the given degrees of freedom equal to j-1, where j is the number of different regions in the channel being analyzed, using the following equation,X2 = Σj (Dobserved-Dexpected)2/Dexpected(4)
Second, we tested whether the distribution of AAMs in different regions was related to the distribution of evolutionarily conserved sites. If AAMs are overrepresented at conserved sites, the number of AAMs for a given region will be proportional to the number of conserved sites found within that region. The expected number of mutations per region was determined using the following equation,Djexpected = Σ0i ((aij/ai) × Diobserved)(5)
where ai is the total number of sites in a protein belonging to variability class i, aij is the number of sites of variability class i found within the region j, and Diobserved is the total number of disease mutations found at variability sites, i, across the entire protein.
Chemical severity of amino acid changes
The interspecific chemical severity of a given site was determined by the average severity (according to Grantham's Scale ) of all ancestor-descendent amino acid differences at that site throughout the tree, as reported by PAML. Only those interspecific changes that result from a single point mutation were included and each type of amino acid change at a given site was counted once to account for common ancestry . The expected chemical severity at each site in HERG and KCNQ1 was determined by computing the average severity of all non-synonymous changes produced by a single point mutation from the human reference codon .
Weighted average for amino acid expected chemical severity
To examine the involvement of specific amino acid residues in disease, the proportion of AAMs at a particular amino acid was calculated as percentage of the total AAMs in both channels. To determine whether these findings were due to an overrepresentation of certain amino acids in the proteins, data was normalized for the total number of codons for a particular residue. Finally, the weighted average for amino acid expected chemical severity was calculated by the sum of the average of each of the individual expected codon chemical severities of the residue multiplied by its contribution to the total number of codons for the residue in the two channels combined.
Channel Structure and Mutation Mapping
Disease mutation distribution by channel region
nsSNPs distribute differently from AAMs in HERG and KCNQ1 (Figure 1). In HERG, nsSNPs are more commonly found in the cytosolic regions (15/16) compared to the transmembrane regions (1/16). In KCNQ1, 8/12 nsSNPs are found in the cytosolic regions and 4/12 nsSNPs are found in the transmembrane region. Most of these occur at sites that do not have associated disease mutations, suggesting that location is an important determinant of disease. In both channels, however, two sites situated in cytosolic regions harbor both disease and polymorphic mutations, suggesting that amino acid identity plays a role in channel dysfunction at certain sites.
AAMs occur preferentially at sites conserved throughout vertebrate evolution and at those with lower evolutionary rates of change
The majority of AAMs mapped to sites completely conserved throughout the evolution of the respective channel: 146/172 mutations (84.9%) in HERG, mapping to 114/906 (12.6%) of the sites used in the analysis, and 153/174 (87.9%) in KCNQ1, mapping to 110/517 (21.3%) of the sites used in the analysis. Conversely, 22/172 AAMs in HERG and 20/174 AAMs in KCNQ1, map to sites with interspecific variation (Figure 1). Only two sites in HERG and one in KCNQ1 involve an amino acid change that is observed in both the AAM and in interspecific change. Because the amino acid sequences of HERG and KCNQ1 are highly conserved (over 65% complete identity between fish and human sequences in both channels for the regions used), a neutral mutational process would still produce a large number of mutations at conserved sites.
To ascertain directly whether AAMs associate preferentially with sites that experience low rates of evolutionary change, we utilized an approach similar to the site analysis used above together with codon evolutionary rate obtained from CODEML in PAML . For both channels, AAMs are found at sites with lower evolutionary rates (Figure 3c). In both channels, certain AAMs were located at codons belonging to the 1st and 2nd evolutionary rate categories as well as to variability classes 0, 1 and 2. This implies that, when AAMs are found at variable sites, they may preferentially occur at those with low evolutionary rates.
Disease mutations are not equally distributed among functional regions of the channels
To examine whether the overall conservation of a region (average variability/site in domain) exerts a non-additive influence on the disease susceptibility , the average number of AAMs per site in a given region of HERG and KCNQ1 was plotted against its average variability per site. These values were correlated for KCNQ1, but not for HERG (Figure 4b). The slope of this relationship in KCNQ1 was greater than those based on the uniform or evolutionary hypothesis. This indicates that AAMs are overabundant in conserved regions and less than expected in variable regions, which suggests that the regional sequence conservation may play a role in KCNQ1 disease susceptibility.
The chemical severities of AAMs are different than changes observed throughout evolution
Average chemical differences of amino acid changes
n = 172
n = 16
67.37 (1.97) *
n = 296
n = 174
52.25 (8.66) †
n = 12
58.32 (2.74) *
n = 168
The chemical severity of amino acid changes may depend on the expected chemical severity of mutations that can arise from a single nucleotide substitution in the codon involved . For KCNQ1 and HERG, the chemical severity of AAMs and expected chemical severity were correlated (Figure 5b). Interspecific and expected chemical severities, however, were not correlated for KCNQ1, and had a very small correlation coefficient for HERG (Figure 5c) (which becomes non-significant when three outliers are removed; not shown). These findings are consistent with previous studies, which suggest that interspecific chemical severity is not influenced by the expected chemical severity of the codon involved but rather by the process of natural selection . Nonetheless, the chemical severities of changes tolerated at variable sites may be influenced by the codon's expected chemical severity. When completely conserved sites were removed (Figure 5d), the expected chemical severity of the involved codon was correlated with interspecific chemical severity, for both HERG and KCNQ1, but to a lesser extent than with disease chemical severity based on slope and correlation co-efficient.
We next compared the average chemical severities of AAMs and nsSNPs (Table 3). We found a significant difference between AAMs and nsSNPs in KCNQ1, but not in HERG, suggesting that other factors, such as location, may play a larger role in causing channel dysfunction for AAMs in HERG compared to KCNQ1. There were no significant differences between nsSNP and interspecific severities in either channel.
Involvement of specific amino acids in arrhythmogenic disease
The high involvement of some residues in disease may be due to a high mutation rate that occurs at CpG dinucleotides  which is a result of a cytosine to thymine transition. This transition is possible in triplets coding for only five of the twenty amino acids, which are: arginine (4/6 codons), serine (1/6 codons), proline (1/4 codons), threonine (1/4 codons) and alanine (1/4 codons). Of the total number of AAMs in HERG and KCNQ1, only 13.3% are due to a C/T transition at a CpG dinucleotide. Of all arginine mutations resulting in disease, 47% are due to this specific nucleotide transition. These numbers suggest that CpG dinucleotides may contribute to the disease process but that this is not the only factor responsible for the high numbers of mutations in this sub-set of amino acids. A role for factors other that CpG dinucleotide hypermutability is also supported by the high involvement of glycines and tryptophans in disease (Figure 6b), which do not possess CpG dinucleotides.
Our findings that the overall chemical severity of AAMs is greater than those of interspecific variation and SNPs, and that the chemical severity of AAMs correlates with expected chemical severity at those individual codons, suggest that chemical severity is on a continuum and that a threshold severity exists which, once crossed, results in disease. We might then expect that whether mutations at specific amino acids cause disease at all also depends on the expected chemical severity of the involved codon. To examine this, the proportion of total AAMs at a particular residue was plotted against the calculated weighted average of expected chemical severity (Figure 6c). A significant, positive correlation was found between these suggesting that expected chemical severity of a site contributes to the probability of obtaining a disease mutation, as well as to the severity of that mutation.
In this study, we show that that AAMs are overabundant at evolutionarily conserved and slowly evolving sites, which are likely critical for channel function and thus intolerant of changes in amino acid sequence. Because KCNQ1 and HERG are highly conserved, an underlying neutral mutational process could produce large numbers of mutations at conserved sites. Our data provide the quantitative backing to support an over representation of AAMs at evolutionarily conserved positions in these channels. A smaller than expected, but still substantial, numbers of AAMs were found at sites that show interspecific variation. However, only two sites in HERG and one in KCNQ1 are converted to residues that are both an AAM and an interspecific change. Thus, the identity of the residue, rather than its location, is most responsible for producing the disease at sites that have undergone evolutionary change.
HERG and KCNQ1 possess structurally defined domains with specific functions that have been strongly preserved throughout the course of evolution [36, 38]. We found that AAMs are found preferentially in some domains, even after accounting for their size and evolutionary conservation. In both channels, the pore region possesses an overabundance of mutations, while both N- and C-termini possess an underabundance of mutations. The extracellular linker region had an overabundance of mutations in HERG, whereas the intracellular linker had an overabundance of mutations in KCNQ1. This implies that the extracellular linker may be more important to overall function in HERG, whereas the intracellular linker is more critical for function in KCNQ1. The overabundance of mutations in the PAS domain, not present in KCNQ1, highlights its functional importance and suggests that that the addition of a functionally important domain in a protein can increase susceptibility to disease. Overall regional conservation, which takes into account the average variability per site in a domain, contributes to the uneven regional distribution of disease-causing mutations in some proteins , but we found this only for KCNQ1 (Figure 4b). Therefore, factors other than domain size, site conservation and regional conservation must influence the domain specific distribution of AAMs in HERG, and possibly in KCNQ1 as well.
NsSNPs are underabundant at conserved sites in HERG, but distribute randomly in KCNQ1. The latter finding may be due to a small sample size, although the numbers of nsSNPs analyzed in HERG were similarly small. This difference between the channels may be because nsSNPs are less disruptive when located at conserved sites in KCNQ1, or because conserved sites in HERG are less tolerant to amino acid changes. Nevertheless, the distribution of analyzed nsSNPs in both channels suggests that they are phenotypically neutral. The different patterns of nsSNP distribution between the KCNQ1 and HERG underscore the need to identify and quantitatively analyze the distribution of more nsSNPs on each protein, and to ascertain their impact on channel function and arrhythmia susceptibility. In HERG, it is known that some, but not all, polymorphisms may alter channel function , and also contribute to an increased QT interval duration .
In both HERG and KCNQ1, AAMs are chemically more severe than interspecific changes tolerated throughout evolution or polymorphisms that are not associated with disease. Codons with a higher expected chemical severity are associated with disease mutations with a high chemical severity. These data are in keeping with those found previously in rhodopsin  and suggest that the intrinsic potential of the involved codon contributes to disease chemical severity.
We also provide novel evidence that the expected chemical severity of a codon contributes to the overrepresentation of certain amino acids in HERG and KCNQ1 AAMs, especially arginine, tryptophan and glycine. In a recent study of 437 proteins, these three amino acids were also highly overrepresented in DAMs . Therefore, we predict that expected chemical severity plays an important role in determining the propensity of a given codon to cause disease in many proteins, in addition to other factors such as the residue's roles in biosynthesis, function and stability of the channels. Finally, the chemical severities of interspecific changes in KCNQ1 and HERG also correlate with those expected for the codons, when only evolutionarily variable sites are considered. These novel data argue that, in addition to a predicted role of natural selection, the expected chemical severity of the codon contributes to variation observed over the course of evolution in these channels. Despite the presumed predominance of natural selection, the genetic code has been shown to influence the mutational process when evolutionary divergence is low . These data are significant given the uncertainty as to the role of natural selection versus non-adaptive forces in shaping genotypic and phenotypic variation [42, 43].
Our analyses may be influenced by the fact that some mutations may, in a systematic way, never be detected. For example, AAMs that lead to death before natural birth, or to SUDs, may never be identified unless the fetus or victim is subsequently screened for arrhythmogenic mutations. This could, in turn, reduce the number of observed mutations unique to certain sites or functional domains. The evaluation of AAMs has broadened, and identification of mutations in KCNQ1 and HERG, and in other candidate genes, associated with SIDs and SUDs has been carried out [14–17]. Identification of AAMs in a more broad population may reveal a different sub-set of mutations that localize to unique regions within the channels, or more strongly support the susceptibility of sites and functional domains identified in this study to arrhythmogenic mutations.
Our study represents the first quantitative evolutionary and chemical severity analysis of AAMs in the HERG and KCNQ1 potassium channel genes. Unlike nsSNPs, AAMs preferentially locate to evolutionarily conserved, and functionally important, sites and regions within HERG and KCNQ1, and are chemically more severe than changes which occur in evolution. Expected chemical severity may contribute to the overrepresentation of certain residues in AAMs, as well as to changes observed throughout evolution.
Our findings, together with those from other studies, suggest that novel DAMs and AAMs may be recognized quickly by surveying naturally occurring variation among species . If a SNP identified in an individual does not appear in other species at that position, then it is likely to be disease-causing. The location of AAMs (to conserved or variable regions and/or residues) may correlate with clinical severity or other characteristics of the diseases. In the case of Long QT syndrome, genotype and specific mutations have been shown to contribute to phenotype [44–47] and the underlying genetic defects contribute to risk stratification, prevention and therapy [12, 48]. Unfortunately, there is still considerable variation in Long QT phenotype and age of disease onset . Therefore, continued discovery and mapping of mutations, as done in this study, along with parallel studies on disease phenotype will ultimately lead to a better understanding of the genotype-phenotype relationship, help to better predict the outcome of novel disease mutations and aid in development of mutation-specific therapies.
We would like to acknowledge Christian Peters for helpful feedback on the manuscript as well as James Ho, Chris Fung, Thomas Li, and Navinder Dhaliwat for preliminary work and discussion on the project. Furthermore, we would like to acknowledge three anonymous reviewers that provided helpful feedback on the manuscript.
Supported by grants from the Natural Sciences and Engineering Research Council of Canada (HAJ and EAA). E.A.A. is also the recipients of a Tier II Canada Research Chair.
- Sanguinetti MC, Curran ME, Zou A, Shen J, Spector PS, Atkinson DL, Keating MT: Coassembly of K(V)LQT1 and minK (IsK) proteins to form cardiac I(Ks) potassium channel. Nature. 1996, 384 (6604): 80-83.View ArticlePubMedGoogle Scholar
- Sanguinetti MC, Jiang C, Curran ME, Keating MT: A mechanistic link between an inherited and an acquired cardiac arrhythmia: HERG encodes the IKr potassium channel. Cell. 1995, 81 (2): 299-307.View ArticlePubMedGoogle Scholar
- Trudeau MC, Warmke JW, Ganetzky B, Robertson GA: HERG, a human inward rectifier in the voltage-gated potassium channel family. Science. 1995, 269 (5220): 92-95.View ArticlePubMedGoogle Scholar
- Barhanin J, Lesage F, Guillemare E, Fink M, Lazdunski M, Romey G: K(V)LQT1 and lsK (minK) proteins associate to form the I(Ks) cardiac potassium current. Nature. 1996, 384 (6604): 78-80.View ArticlePubMedGoogle Scholar
- Brugada R, Hong K, Dumaine R, Cordeiro J, Gaita F, Borggrefe M, Menendez TM, Brugada J, Pollevick GD, Wolpert C, Burashnikov E, Matsuo K, Wu YS, Guerchicoff A, Bianchi F, Giustetto C, Schimpf R, Brugada P, Antzelevitch C: Sudden death associated with short-QT syndrome linked to mutations in HERG. Circulation. 2004, 109 (1): 30-35.View ArticlePubMedGoogle Scholar
- Gaita F, Giustetto C, Bianchi F, Wolpert C, Schimpf R, Riccardi R, Grossi S, Richiardi E, Borggrefe M: Short QT Syndrome: a familial cause of sudden death. Circulation. 2003, 108 (8): 965-970.View ArticlePubMedGoogle Scholar
- Bellocq C, van Ginneken AC, Bezzina CR, Alders M, Escande D, Mannens MM, Baro I, Wilde AA: Mutation in the KCNQ1 gene leading to the short QT-interval syndrome. Circulation. 2004, 109 (20): 2394-2397.View ArticlePubMedGoogle Scholar
- Chen YH, Xu SJ, Bendahhou S, Wang XL, Wang Y, Xu WY, Jin HW, Sun H, Su XY, Zhuang QN, Yang YQ, Li YB, Liu Y, Xu HJ, Li XF, Ma N, Mou CP, Chen Z, Barhanin J, Huang W: KCNQ1 gain-of-function mutation in familial atrial fibrillation. Science. 2003, 299 (5604): 251-254.View ArticlePubMedGoogle Scholar
- Hong K, Bjerregaard P, Gussak I, Brugada R: Short QT syndrome and atrial fibrillation caused by mutation in KCNH2. J Cardiovasc Electrophysiol. 2005, 16 (4): 394-396.View ArticlePubMedGoogle Scholar
- Wang Q, Curran ME, Splawski I, Burn TC, Millholland JM, VanRaay TJ, Shen J, Timothy KW, Vincent GM, de Jager T, Schwartz PJ, Toubin JA, Moss AJ, Atkinson DL, Landes GM, Connors TD, Keating MT: Positional cloning of a novel potassium channel gene: KVLQT1 mutations cause cardiac arrhythmias. Nat Genet. 1996, 12 (1): 17-23.View ArticlePubMedGoogle Scholar
- Curran ME, Splawski I, Timothy KW, Vincent GM, Green ED, Keating MT: A molecular basis for cardiac arrhythmia: HERG mutations cause long QT syndrome. Cell. 1995, 80 (5): 795-803.View ArticlePubMedGoogle Scholar
- Schwartz PJ: Management of long QT syndrome. Nat Clin Pract Cardiovasc Med. 2005, 2 (7): 346-351.View ArticlePubMedGoogle Scholar
- Roden DM: Clinical practice. Long-QT syndrome. N Engl J Med. 2008, 358 (2): 169-176.View ArticlePubMedGoogle Scholar
- Schwartz PJ, Priori SG, Bloise R, Napolitano C, Ronchetti E, Piccinini A, Goj C, Breithardt G, Schulze-Bahr E, Wedekind H, Nastoli J: Molecular diagnosis in a child with sudden infant death syndrome. Lancet. 2001, 358 (9290): 1342-1343.View ArticlePubMedGoogle Scholar
- Arnestad M, Crotti L, Rognum TO, Insolia R, Pedrazzini M, Ferrandi C, Vege A, Wang DW, Rhodes TE, George AL, Schwartz PJ: Prevalence of long-QT syndrome gene variants in sudden infant death syndrome. Circulation. 2007, 115 (3): 361-367.View ArticlePubMedGoogle Scholar
- Wang DW, Desai RR, Crotti L, Arnestad M, Insolia R, Pedrazzini M, Ferrandi C, Vege A, Rognum T, Schwartz PJ, George AL: Cardiac sodium channel dysfunction in sudden infant death syndrome. Circulation. 2007, 115 (3): 368-376.View ArticlePubMedGoogle Scholar
- Tester DJ, Ackerman MJ: Postmortem long QT syndrome genetic testing for sudden unexplained death in the young. J Am Coll Cardiol. 2007, 49 (2): 240-246.View ArticlePubMedGoogle Scholar
- Grantham R: Amino acid difference formula to help explain protein evolution. Science. 1974, 185 (4154): 862-864.View ArticlePubMedGoogle Scholar
- Botstein D, Risch N: Discovering genotypes underlying human phenotypes: past successes for mendelian disease, future approaches for complex disease. Nat Genet. 2003, 33 Suppl: 228-237.View ArticlePubMedGoogle Scholar
- Krawczak M, Ball EV, Cooper DN: Neighboring-nucleotide effects on the rates of germ-line single-base-pair substitution in human genes. Am J Hum Genet. 1998, 63 (2): 474-488.PubMed CentralView ArticlePubMedGoogle Scholar
- Briscoe AD, Gaur C, Kumar S: The spectrum of human rhodopsin disease mutations through the lens of interspecific variation. Gene. 2004, 332: 107-118.View ArticlePubMedGoogle Scholar
- Miller MP, Kumar S: Understanding human disease mutations through the use of interspecific genetic variation. Hum Mol Genet. 2001, 10 (21): 2319-2328.View ArticlePubMedGoogle Scholar
- Subramanian S, Kumar S: Evolutionary anatomies of positions and types of disease-associated and neutral amino acid mutations in the human genome. BMC Genomics. 2006, 7: 306-PubMed CentralView ArticlePubMedGoogle Scholar
- Mooney SD, Klein TE: The functional importance of disease-associated mutation. BMC Bioinformatics. 2002, 3: 24-PubMed CentralView ArticlePubMedGoogle Scholar
- Miller MP, Parker JD, Rissing SW, Kumar S: Quantifying the intragenic distribution of human disease mutations. Ann Hum Genet. 2003, 67 (Pt 6): 567-579.View ArticlePubMedGoogle Scholar
- Nei M: Selectionism and neutralism in molecular evolution. Mol Biol Evol. 2005, 22 (12): 2318-2342.PubMed CentralView ArticlePubMedGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.View ArticlePubMedGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25 (24): 4876-4882.PubMed CentralView ArticlePubMedGoogle Scholar
- Inherited Arrhythmias Database. [http://www.fsm.it/cardmoc/]
- Mank-Seymour AR, Richmond JL, Wood LS, Reynolds JM, Fan YT, Warnes GR, Milos PM, Thompson JF: Association of torsades de pointes with novel and known single nucleotide polymorphisms in long QT syndrome genes. Am Heart J. 2006, 152 (6): 1116-1122.View ArticlePubMedGoogle Scholar
- Aydin A, Bahring S, Dahm S, Guenther UP, Uhlmann R, Busjahn A, Luft FC: Single nucleotide polymorphism map of five long-QT genes. J Mol Med. 2005, 83 (2): 159-165.View ArticlePubMedGoogle Scholar
- Millat G, Chevalier P, Restier-Miron L, Da Costa A, Bouvagnet P, Kugener B, Fayol L, Gonzalez Armengod C, Oddou B, Chanavat V, Froidefond E, Perraudin R, Rousson R, Rodriguez-Lafrasse C: Spectrum of pathogenic mutations and associated polymorphisms in a cohort of 44 unrelated patients with long QT syndrome. Clin Genet. 2006, 70 (3): 214-227.View ArticlePubMedGoogle Scholar
- Felsenstein J: Inferring phylogenies from protein sequences by parsimony, distance, and likelihood methods. Methods Enzymol. 1996, 266: 418-427.View ArticlePubMedGoogle Scholar
- Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci. 1997, 13 (5): 555-556.PubMedGoogle Scholar
- Splawski I, Shen J, Timothy KW, Lehmann MH, Priori S, Robinson JL, Moss AJ, Schwartz PJ, Towbin JA, Vincent GM, Keating MT: Spectrum of mutations in long-QT syndrome genes. KVLQT1, HERG, SCN5A, KCNE1, and KCNE2. Circulation. 2000, 102 (10): 1178-1185.View ArticlePubMedGoogle Scholar
- Sanguinetti MC, Tristani-Firouzi M: hERG potassium channels and cardiac arrhythmia. Nature. 2006, 440 (7083): 463-469.View ArticlePubMedGoogle Scholar
- Vitkup D, Sander C, Church GM: The amino-acid mutational spectrum of human genetic disease. Genome Biol. 2003, 4 (11): R72-PubMed CentralView ArticlePubMedGoogle Scholar
- Jespersen T, Grunnet M, Olesen SP: The KCNQ1 potassium channel: from gene to physiological function. Physiology (Bethesda). 2005, 20: 408-416.View ArticleGoogle Scholar
- Anson BD, Ackerman MJ, Tester DJ, Will ML, Delisle BP, Anderson CL, January CT: Molecular and functional characterization of common polymorphisms in HERG (KCNH2) potassium channels. Am J Physiol Heart Circ Physiol. 2004, 286 (6): H2434-41.View ArticlePubMedGoogle Scholar
- Newton-Cheh C, Guo CY, Larson MG, Musone SL, Surti A, Camargo AL, Drake JA, Benjamin EJ, Levy D, D'Agostino RB, Hirschhorn JN, O'Donnell C J: Common genetic variation in KCNH2 is associated with QT interval duration: the Framingham Heart Study. Circulation. 2007, 116 (10): 1128-1136.View ArticlePubMedGoogle Scholar
- Benner SA, Cohen MA, Gonnet GH: Amino acid substitution during functionally constrained divergent evolution of protein sequences. Protein Eng. 1994, 7 (11): 1323-1332.View ArticlePubMedGoogle Scholar
- Yi SV: Non-adaptive evolution of genome complexity. Bioessays. 2006, 28 (10): 979-982.View ArticlePubMedGoogle Scholar
- Lynch M: The evolution of genetic networks by non-adaptive processes. Nat Rev Genet. 2007, 8 (10): 803-813.View ArticlePubMedGoogle Scholar
- Brink PA, Crotti L, Corfield V, Goosen A, Durrheim G, Hedley P, Heradien M, Geldenhuys G, Vanoli E, Bacchini S, Spazzolini C, Lundquist AL, Roden DM, George AL, Schwartz PJ: Phenotypic variability and unusual clinical severity of congenital long-QT syndrome in a founder population. Circulation. 2005, 112 (17): 2602-2610.View ArticlePubMedGoogle Scholar
- Crotti L, Spazzolini C, Schwartz PJ, Shimizu W, Denjoy I, Schulze-Bahr E, Zaklyazminskaya EV, Swan H, Ackerman MJ, Moss AJ, Wilde AA, Horie M, Brink PA, Insolia R, De Ferrari GM, Crimi G: The common long-QT syndrome mutation KCNQ1/A341V causes unusually severe clinical manifestations in patients with different ethnic backgrounds: toward a mutation-specific risk stratification. Circulation. 2007, 116 (21): 2366-2375.View ArticlePubMedGoogle Scholar
- Moss AJ, Zareba W, Kaufman ES, Gartman E, Peterson DR, Benhorin J, Towbin JA, Keating MT, Priori SG, Schwartz PJ, Vincent GM, Robinson JL, Andrews ML, Feng C, Hall WJ, Medina A, Zhang L, Wang Z: Increased risk of arrhythmic events in long-QT syndrome with mutations in the pore region of the human ether-a-go-go-related gene potassium channel. Circulation. 2002, 105 (7): 794-799.View ArticlePubMedGoogle Scholar
- Tan HL, Bardai A, Shimizu W, Moss AJ, Schulze-Bahr E, Noda T, Wilde AA: Genotype-specific onset of arrhythmias in congenital long-QT syndrome: possible therapy implications. Circulation. 2006, 114 (20): 2096-2103.View ArticlePubMedGoogle Scholar
- Priori SG, Schwartz PJ, Napolitano C, Bloise R, Ronchetti E, Grillo M, Vicentini A, Spazzolini C, Nastoli J, Bottelli G, Folli R, Cappelletti D: Risk stratification in the long-QT syndrome. N Engl J Med. 2003, 348 (19): 1866-1874.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.