Exome sequencing helped the fine diagnosis of two siblings afflicted with atypical Timothy syndrome (TS2)
- Sebastian Fröhler†1,
- Moritz Kieslich†2,
- Claudia Langnick1,
- Mirjam Feldkamp1,
- Bernd Opgen-Rhein2,
- Felix Berger2,
- Joachim C Will†2Email author and
- Wei Chen†1Email author
© Fröhler et al.; licensee BioMed Central Ltd. 2014
Received: 12 June 2013
Accepted: 25 April 2014
Published: 29 April 2014
Long-QT syndrome (LQTS) causes a prolongation of the QT-interval in the ECG leading to life threatening tachyarrhythmia and ventricular fibrillation. One atypical form of LQTS, Timothy syndrome (TS), is associated with syndactyly, immune deficiency, cognitive and neurological abnormalities as well as distinct cranio-facial abnormalities.
On a family with both children diagnosed with clinical LQTS, we performed whole exome sequencing to comprehensively screen for causative mutations after a targeted candidate gene panel screen for Long-QT syndrome target genes failed to identify any underlying genetic defect. Using exome sequencing, we identified in both affected children, a p.402G > S mutation in exon 8 of the CACNA1C gene, a voltage-dependent Ca2+ channel. The mutation was inherited from their father, a mosaic mutation carrier. Based on this molecular finding and further more careful clinical examination, we refined the diagnosis to be Timothy syndrome (TS2) and thereby were able to present new therapeutic approaches.
Our study highlights the difficulties in accurate diagnosis of patients with rare diseases, especially those with atypical clinical manifestation. Such challenge could be addressed with the help of comprehensive and unbiased mutation screening, such as exome sequencing.
KeywordsTimothy syndrome Exome sequencing Mosaic mutation LQTS
Long-QT syndrome [1, 2] (LQTS) is a cardiac chanellopathy characterized by QT prolongation and T-wave abnormalities in the ECG causing syncope, life threatening arrhythmia and sudden cardiac death due to torsade de pointes (TdP) and ventricular fibrillation.
On a molecular level, most LQTS patients have mutations in a set of ten genes (KCNQ1, KCNH2, SCN5A, KCNE1, KCNE2, CAV3, SCN4B, AKAP9, SNTA1 and KCNJ5) with ~97% of all mutations within this gene panel falling into the first three genes . The remaining seven genes represent some rare cases of LQTS. However, in about 25% of all studied families showing LQTS symptoms, no mutations within the set of 10 candidate genes can be found , indicating either a limitation in the testing methods employed on the identified candidate genes or the presence of additional genes related to LQTS.
A distinct form of LQTS is the syndactyly related TS (also called LQT8) a rare autosome dominant disorder. The first clinical cases were described as the ‘heart-hand syndrome’ in the 1990’s [4–6] with a prolonged QT-interval and syndactyly being present in 100% of described cases. Additional cardinal features of TS include: bradycardia, AV-Block, congenital heart disease, immune deficiency, hypoglycemia, neurological and cognitive abnormalities such as intellectual disability, autism and seizures as well as distinctive cranio-facial abnormalities. Also, an association of TS with hypertrophic cardiomyopathy has been described .
In 2004, a missense mutation (p.406G > R) located in exon 8A of the gene CACNA1C, encoding the transmembrane segment S6 of domain I of CaV1.2 (Voltage-dependent L-type calcium channel subunit alpha-1C), has been identified to cause TS . The mutation leads to almost complete loss of voltage dependent channel-inactivation, thereby resulting in intracellular calcium overload and subsequently delayed repolarization of the respective cell. The vast majority of described cases result from de-novo missense mutations. Nevertheless, in one family with two affected children, mosaicism for the p.406G > R mutation was discovered in the mother . In 2005, Splawski et al. described two unrelated individuals with severe prolongation of the QTc-interval without syndactyly, in which two de novo mutations (p.406G > R and p.402G > S) in Exon 8 of CACNA1C have been identified to cause a related but distinct phenotype, classified as Timothy syndrome 2 (TS2) . In contrast to the splice variant containing Exon 8A, Exon 8 is widely expressed in heart and brain, consistent with the phenotype of the described patients. Whereas patient one, described by Splawski, presented with severe central nervous system (CNS) and cardiac affection, patient two presented with severe QTc prolongation but milder CNS affection which was attributed to the somatic mosaicism for the mutation.
In this report, we describe a family with two affected children initially diagnosed to have LQTS. After a candidate gene panel test for LQTS failed to identify any causative genetic defect, we referred to whole-exome sequencing and identified a causative p.402G > S (exon 8) mutation in the CACNA1C gene in both affected children. The mutation was inherited from the father, who is mosaic for the specific mutation. As a result, we refined the diagnosis to be Timothy syndrome type 2 (TS2) although both patients do not manifest the full clinical spectrum of classical TS2. On one hand, our study highlights the difficulties in accurate diagnosis of patients with rare diseases, especially those with atypical clinical manifestation. On the other hand, we demonstrated such challenge could be addressed with the help of comprehensive and unbiased mutation screening, such as exome sequencing.
The index patient is the 8 year old daughter (Patient 1) who presented VF/ TdP leading to cardiac arrest and was successful resuscitated at the age of 2 months. Antiarrhythmic therapy was started with beta blockers, nicorandil and mexiletine, controlling TdP best with high doses of propranolol and mexiletine. At the age of 4, her parents agreed to have an internal cardiac defibrillator (ICD) implanted, who delivered 2 adequate shock since then. The patient developed generalized epilepsy which is currently well controlled with levetiracetam. Today, at eight years of age, the patient is unable to speak and learned walking only at the age of seven. She additionally shows evidence of a vast precocious puberty which is not known to be associated with TS. The neurological findings (epilepsy, developmental delay) are rather attributed to hypoxic brain damage caused by cardiac arrest and resuscitation than being part of TS. Especially since the early neonatal development of the patient has been entirely unremarkable. Moreover, patient 2 unlike his sister reached the milestones of psychomotor development at the expected times and clinical workup failed to show any neurological abnormalities. He is the second child of the family and is now six years old. He developed VF/ TdP at the age of 3 years and was successfully resuscitated using his sister’s automatic external defibrillator and received an ICD two months later. As his sister, he was put on propranolol and mexiletine. Since then, the patient has had frequent episodes of VF and non-sustained ventricular tachycardia and has received overall more than 40 adequate and appropriate shocks. His antiarrhythmic therapy has been modified several time, he was administered several antiarrhythmic drugs (AAD): high-dose betablockers, potassium-, magnesium- and calcium-antagonists as well as Vaughn-Williams class Ia, Ib and Ic agents. None of these AAD combinations shortened QTc-intervals or reduced VF episodes. Left- and right-cardiac sympathetic denervation could also not reduce the arrhythmic episodes.
Today, quality of life of both patients is still significantly affected by the continuing relapses of VF, especially in the boy, who seems to be especially insusceptible to standard therapeutic approaches.
Written informed consent was obtained from each family member for the subsequent studies according to the German Genetic Diagnostic Law (GenDG).
Candidate gene panel screening for LQTS
Given the diagnosis of LQTS for both children, LQTS candidate gene panels I, II and III, including genes: KCNQ1, KCNH2, SCN5A, KCNE1 and KCNE2 were PCR amplified and Sanger-sequenced to identify mutations. Neither this screen, nor a multiplex ligation probe amplification (MLPA) of 40 regions of the target genes identified any mutations or insertions/ deletions in the regions screened (both analysis were performed at the Center for Human Genetics and Laboratory Medicine, Martinsried, Germany).
Exome sequencing and bioinformatics
All four family members were subjected to exome sequencing. Genomic DNA was isolated from blood samples using standard methods. Five micrograms of genomic DNA were enriched using the Agilent Human All Exon V4 kit (Agilent Technologies, Santa Clara, CA, USA) following the manufacturer's protocol. A list of screened genes and a list of assayed regions for these probes is available on the Agilent website (https://earray.chem.agilent.com/suredesign/).
Whole-exome libraries were sequenced on an Illumina HiSeq 2000 system for 1x101 cycles following the manufacturer's instructions (Illumina, San Diego, CA, USA). All raw sequencing reads were mapped onto UCSC hg19  using BWA 0.5.9-r16 , mappings were converted into BAM file format using samtools 0.1.18 . Initial mappings were post-processed using GATK 1.3-21  following their ‘best practices V3’. In brief, reads were realigned around sites of known insertions and deletions (INDELs). Then, likely PCR duplicates were detected using Picard 1.48 . Finally, raw base quality scores were empirically recalibrated. Single Nucleotide Polymorphisms (SNPs) and short INDELs were identified using the UnifiedGenotyper from GATK and in parallel using the variant caller included in samtools. Variants were classified as novel or known variants deposited in dbSNP 135 . Functional consequences of each variant were annotated using snpEff 2.0.5d  for UCSC hg19 RefSeq genes and ENSEMBL 65 human gene models [10, 16]. The potential deleterious effect was evaluated using PolyPhen 2 , SIFT , PhyloP , MutationTaster , GERP++ , LRT  and OMIM  if available. Candidate variants from GATK and from samtools were compared to increase both sensitivity and specificity.
To confirm the presence and identity of the p.402G > S variant in CACNA1C exon 8, Sanger sequencing was performed on PCR amplicons from genomic DNA covering that variant position. A restriction assay on the same amplicons using the enzyme DdeI that specifically cleaves amplicons carrying the mutation was also performed to validate the mosaicism in the father.
Summary statistics for exome sequencing
# mapped reads
# uniquely mapped reads
# non-redundant uniquely mapped reads
# non-redundant reads uniquely mapped in targeted regions
Average sequencing depth in target regions
Using GATK and samtools, we detected 26,501 – 30,024 SNPs and INDELs in the exome of each family member, of which 94.3 – 97.1% are known variants deposited in dbSNP 135 (see Additional file 2). Given the pedigree, we first searched for non-silent variants present in both affected children but absent in both apparently healthy parents. Whereas no such variants were found by GATK, samtools outputs one such variant in exon 8 of the gene CACNA1C (c.1204G > A, p.402G > S).
We then checked for the reason leading to the inconsistent results from samtools and GATK. While both tools called the same genotypes in affected siblings as well as their mother, the father was identified to carry heterozygous or homozygous reference allele by GATK (PHRED-scaled genotype quality 99) and samtools (PHRED-scaled genotype quality 50), respectively. After retrieving the original mapping data from the father, we found 15 high-quality reads (9.9% of all 151 reads aligned at this position) carrying the variant allele. This relatively low non-reference allele frequency, indicating a mosaic carrier, allowed GATK to call a heterozygous genotype, whereas it does not pass the default samtools’ variant calling threshold.
Sanger sequencing on the blood samples of all family members not only confirmed the heterozygous p.402G > S variant in both affected children but also the mosaic genotype for this variant in the father (Figure 1b). The father's mosaic genotype was also observed in his oral mucosa swap sample at a similar allele frequency (see Additional file 3). Both heterozygous and mosaic genotypes could also be validated by a restriction assay using DdeI, specifically cleaving amplicons containing the variant (Figure 1c).
Refined diagnosis and treatment
The p.G402S variant in exon 8 of the gene CACNA1C has been previously associated with Timothy Syndrome TS2 (MIM: 601005). In literature, treatment with multipotent ion-channel-blocker ranozaline in combination with verapamil has been reported to be beneficial in a patient with this distinct mutation . This therapeutical approach is currently applied in our patients with short follow up time.
Cell-culture-based experiments with Cav1.2. and Timothy channels show that calcineurin inhibitors (e.g. FK 506, cyclosporine A) stabilize the action potential of the respective cells by increasing voltage dependent channel-inactivation . Theoretically this could be another therapeutical approach in case of TS refractory to conventional medical therapy.
Using whole-exome sequencing we identified a p.402G > S mutation in exon 8 of the gene CACNA1C in a family with two affected children, both presenting atypical form of Timothy syndrome (TS2). This variant has been previously described to alter the transmembrane segment S6 of domain I in the CACNA1C gene, thereby almost completely disrupting voltage-dependence of this Ca2+-channel, leading to a massively increased QT-interval in an electrocardiogram . The CACNA1C isoform containing exon 8 is reported to be the predominant isoform expressed in heart and brain, consistent with the phenotype of both patients.
A previous LQTS candidate gene panel screen returned no diagnostic result since, at the time of evaluation, both children only unequivocally manifest the LQTS part of the Timothy Syndrome phenotype, together with an only partial syndactyly that they even share with their father - a feature quite easy to oversee by a cardiologist focusing on their severe heart arrhythmia. Exome sequencing, as an unbiased diagnostic aid however, clearly identified the causative mutation after a careful analysis.
The p.402G > S mutation was inherited from the mosaic father. However, a mosaic genotype, manifested with a relatively low non-reference allele frequency in the exome sequencing results is not susceptible to most variant calling algorithms. Interestingly, although samtools lead to the discovery of this mutation, a (somehow) wrong genotype call was made by samtools, whereas the correct genotype call of GATK resulted in filtering out this variant from the list of autosome-dominant events, given the father is not afflicted with Timothy syndrome. Luckily, mosaicism is a rare event and most studies will identify the same causative mutation using different variant callers. However, care should be taken in case of pedigrees with an atypical inheritance pattern where a combination of variant callers with different genotyping models/thresholds can provide a valuable list of discordant genotype calls, thereby enriching for candidate mosaicisms. Those variants might just be filtered out when relying on a single variant caller.
Diagnosing patients with rare diseases is a complex mission that can easily take several years of clinical examinations. Many different gene defects can give rise to clinically indistinguishable phenotypes. Recent data show that up to 40% of all initial diagnosis are wrong [26, 27]. Current whole-exome/whole-genome sequencing technology, applied in clinical settings, can aid in correctly diagnosing those patients, as demonstrated in this study. With the ever decreasing sequencing costs, a global, truly unbiased genomic screen, at a cost even below that of targeted gene panel screens today, will be available in the clinic in only a few years. Combined with a genotype-phenotype-treatment resource, constantly expanding during such kind of studies in the clinic, new therapeutic approaches could be offered to the patient as an immediate benefit. However, in such clinical context, controlled vocabularies for phenotyping patients together with rigid data protection strategies must be warranted at all times. Additional legislative effort might be required for recording, storing and accessing this kind of data. Furthermore, a clear consent must be obtained regarding how to proceed with incidental findings. Whether those should be conveyed to the patient if requested, whether to only convey results on ‘treatable diseases’ or even only to convey findings immediately relevant for the disease studied must be agreed upon before initiating such kind of studies. Finally, diagnosing a ‘syndrome’ might stigmatize the patient, so proper genetic counseling by specialized physicians will be required for all such molecular genetic diagnostics.
We have described a family with two affected children initially diagnosed to have LQTS. After a candidate gene panel screen for LQTS failed, we referred to whole-exome sequencing of the whole family as a global, unbiased method. As a result, we have identified a causative p.402G > S mutation in the CACNA1C gene in both affected children and showed that this mutation was inherited from their mosaic father. Thereby, we could refine the diagnosis for both children to be Timothy Syndrome (TS2), a very rare atypical form of LQTS. This refined diagnosis immediately opened new therapeutic options for the patient. Our study strongly supports the benefit of exome sequencing in clinical diagnosis of complex cases and soon might become a widely used tool in clinical settings.
S.F., C.L., M.F. and W.C are supported by the Helmholtz Association, the German Ministry of Education and Research (BMBF) and the Senate of Berlin by funds to the Berlin Institute for Medical Systems Biology (BIMSB).
- Romano C, Gemme G, Pongiglione R: Rare cardiac arrythmias of the pediatric age. II. syncopal attacks due to paroxysmal ventricular fibrillation. (Presentation of 1st case in Italian pediatric literature). Clin Pediatr (Bologna). 1963, 45: 656-683.Google Scholar
- Ward OC: A new familial cardiac syndrome in children. J Ir Med Assoc. 1964, 54: 103-106.PubMedGoogle Scholar
- NCBI Gene Reviews: Long-QT Syndrome. http://www.ncbi.nlm.nih.gov/books/NBK1129/,
- Reichenbach H, Meister EM, Theile H: The heart-hand syndrome. A new variant of disorders of heart conduction and syndactylia including osseous changes in hands and feet. Kinderarztl Prax. 1992, 60 (2): 54-56.PubMedGoogle Scholar
- Marks ML, Whisler SL, Clericuzio C, Keating M: A new form of long QT syndrome associated with syndactyly. J Am Coll Cardiol. 1995, 25 (1): 59-64. 10.1016/0735-1097(94)00318-K.View ArticlePubMedGoogle Scholar
- Marks ML, Trippel DL, Keating MT: Long QT syndrome associated with syndactyly identified in females. Am J Cardiol. 1995, 76 (10): 744-745. 10.1016/S0002-9149(99)80216-1.View ArticlePubMedGoogle Scholar
- Lo-A-Njoe SM, Wilde AA, van Erven L, Blom NA: Syndactyly and long QT syndrome (CaV1.2 missense mutation G406R) is associated with hypertrophic cardiomyopathy. Heart Rhythm. 2005, 2 (12): 1365-1368. 10.1016/j.hrthm.2005.08.024.View ArticlePubMedGoogle Scholar
- Splawski I, Timothy KW, Sharpe LM, Decher N, Kumar P, Bloise R, Napolitano C, Schwartz PJ, Joseph RM, Condouris K, Tager-Flusberg H, Priori SG, Sanguinetti MC, Keating MT: Ca(V)1.2 calcium channel dysfunction causes a multisystem disorder including arrhythmia and autism. Cell. 2004, 119 (1): 19-31. 10.1016/j.cell.2004.09.011.View ArticlePubMedGoogle Scholar
- Splawski I, Timothy KW, Decher N, Kumar P, Sachse FB, Beggs AH, Sanguinetti MC, Keating MT: Severe arrhythmia disorder caused by cardiac L-type calcium channel mutations. Proc Natl Acad Sci U S A. 2005, 102 (23): 8089-8096. 10.1073/pnas.0502506102. discussion 8086–8088View ArticlePubMedPubMed CentralGoogle Scholar
- Dreszer TR, Karolchik D, Zweig AS, Hinrichs AS, Raney BJ, Kuhn RM, Meyer LR, Wong M, Sloan CA, Rosenbloom KR, Roe G, Rhead B, Pohl A, Malladi VS, Li CH, Learned K, Kirkup V, Hsu F, Harte RA, Guruvadoo L, Goldman M, Giardine BM, Fujita PA, Diekhans M, Cline MS, Clawson H, Barber GP, Haussler D, James Kent W: The UCSC Genome Browser database: extensions and updates 2011. Nucleic Acids Res. 2012, 40 (Database issue): D918-D923.View ArticlePubMedGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25 (14): 1754-1760. 10.1093/bioinformatics/btp324.View ArticlePubMedPubMed CentralGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Subgroup GPDP: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-2079. 10.1093/bioinformatics/btp352.View ArticlePubMedPubMed CentralGoogle Scholar
- DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, Philippakis AA, del Angel G, Rivas MA, Hanna M, McKenna A, Fennell TJ, Kernytsky AM, Sivachenko AY, Cibulskis K, Gabriel SB, Altshuler D, Daly MJ: A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011, 43 (5): 491-498. 10.1038/ng.806.View ArticlePubMedPubMed CentralGoogle Scholar
- Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, Sirotkin K: dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001, 29 (1): 308-311. 10.1093/nar/29.1.308.View ArticlePubMedPubMed CentralGoogle Scholar
- Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, Land SJ, Lu X, Ruden DM: A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin). 2012, 6 (2): 80-92. 10.4161/fly.19695.View ArticleGoogle Scholar
- Flicek P, Amode MR, Barrell D, Beal K, Brent S, Carvalho-Silva D, Clapham P, Coates G, Fairley S, Fitzgerald S, Gil L, García-Girón C, Gordon L, Hourlier T, Hunt S, Juettemann T, Kähäri AK, Keenan S, Komorowska M, Kulesha E, Longden I, Maurel T, McLaren WM, Muffato M, Nag R, Overduin B, Pignatelli M, Pritchard B, Pritchard E, Riat HS, et al: Ensembl 2012. Nucleic Acids Res. 2012, 40 (Database issue): D84-D90.View ArticlePubMedGoogle Scholar
- Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, Kondrashov AS, Sunyaev SR: A method and server for predicting damaging missense mutations. Nat Methods. 2010, 7 (4): 248-249. 10.1038/nmeth0410-248.View ArticlePubMedPubMed CentralGoogle Scholar
- Kumar P, Henikoff S, Ng PC: Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc. 2009, 4 (7): 1073-1081.View ArticlePubMedGoogle Scholar
- Pollard KS, Hubisz MJ, Rosenbloom KR, Siepel A: Detection of nonneutral substitution rates on mammalian phylogenies. Genome Res. 2010, 20 (1): 110-121. 10.1101/gr.097857.109.View ArticlePubMedPubMed CentralGoogle Scholar
- Schwarz JM, Rödelsperger C, Schuelke M, Seelow D: MutationTaster evaluates disease-causing potential of sequence alterations. Nat Methods. 2010, 7 (8): 575-576. 10.1038/nmeth0810-575.View ArticlePubMedGoogle Scholar
- Davydov EV, Goode DL, Sirota M, Cooper GM, Sidow A, Batzoglou S: Identifying a high fraction of the human genome to be under selective constraint using GERP++. PLoS Comput Biol. 2010, 6 (12): e1001025-10.1371/journal.pcbi.1001025.View ArticlePubMedPubMed CentralGoogle Scholar
- Chun S, Fay JC: Identification of deleterious mutations within three human genomes. Genome Res. 2009, 19 (9): 1553-1561. 10.1101/gr.092619.109.View ArticlePubMedPubMed CentralGoogle Scholar
- Online Mendelian Inheritance in Man. http://omim.org,
- Shah DP, Baez-Escudero JL, Weisberg IL, Beshai JF, Burke MC: Ranolazine safely decreases ventricular and atrial fibrillation in Timothy syndrome (LQT8). Pacing Clin Electrophysiol. 2012, 35 (3): e62-64. 10.1111/j.1540-8159.2010.02913.x.View ArticlePubMedGoogle Scholar
- Cohen-Kutner M, Yahalom Y, Trus M, Atlas D: Calcineurin Controls Voltage-Dependent-Inactivation (VDI) of the Normal and Timothy Cardiac Channels. Sci Rep. 2012, 2: 366-View ArticlePubMedPubMed CentralGoogle Scholar
- Ropers HH: On the future of genetic risk assessment. J Community Genet. 2012, 3 (3): 229-236. 10.1007/s12687-012-0092-2.View ArticlePubMedPubMed CentralGoogle Scholar
- EURORDIS: Rare Disease Europe. http://www.eurordis.org,
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/15/48/prepub
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 credited.