Isogenic human pluripotent stem cell pairs reveal the role of a KCNH2 mutation in long-QT syndrome (original) (raw)

Abstract

Patient-specific induced pluripotent stem cells (iPSCs) will assist research on genetic cardiac maladies if the disease phenotype is recapitulated in vitro. However, genetic background variations may confound disease traits, especially for disorders with incomplete penetrance, such as long-QT syndromes (LQTS). To study the LQT2-associated c.A2987T (N996I) KCNH2 mutation under genetically defined conditions, we derived iPSCs from a patient carrying this mutation and corrected it. Furthermore, we introduced the same point mutation in human embryonic stem cells (hESCs), generating two genetically distinct isogenic pairs of LQTS and control lines. Correction of the mutation normalized the current (_I_Kr) conducted by the HERG channel and the action potential (AP) duration in iPSC-derived cardiomyocytes (CMs). Introduction of the same mutation reduced _I_Kr and prolonged the AP duration in hESC-derived CMs. Further characterization of N996I-HERG pathogenesis revealed a trafficking defect. Our results demonstrated that the c.A2987T KCNH2 mutation is the primary cause of the LQTS phenotype. Precise genetic modification of pluripotent stem cells provided a physiologically and functionally relevant human cellular context to reveal the pathogenic mechanism underlying this specific disease phenotype.

Keywords: gene targeting, HERG, human embryonic stem cells, induced pluripotent stem cells, long-QT syndrome

Introduction

Transgenic animal models of human cardiac diseases and cultured cells have been crucial to cardiovascular research, contributing to both our basic understanding of normal development and disease mechanisms. However, these approaches often only partially recapitulate the molecular and cellular phenotypes observed in patients because human heart physiology, as well as the cardiomycyte’s gene and protein expression profile, differs from almost all other experimentally accessible species (Davis et al, 2011).

Human induced pluripotent stem cells (hiPSCs) (Takahashi et al, 2007b; Yu et al, 2007) are creating exciting new opportunities for biomedical research by providing platforms to study mechanisms of genetic disease pathogenesis that could lead to new therapies or reveal drug sensitivities (Braam et al, 2013; Liang et al, 2013). Patient-specific hiPSCs that carry all disease-relevant genetic alterations are important not only for understanding monogenic and complex disease mechanisms in the affected cell types in vitro and identifying potential treatments (Gold-von Simson et al, 2009; Lee et al, 2009; Marchetto et al, 2010; Axelrod et al, 2011; Pini et al, 2012), but also for providing insights into factors that predispose individuals to develop overt symptoms. However, the application of hiPSCs in disease modelling has been limited by the paucity of experimental tools for distinguishing subtle but disease-relevant phenotypic changes from background-related variations.

Inherited long-QT syndrome (LQTS) is a life-threatening, often autosomal dominant, disorder characterized by prolonged ventricular repolarization, a propensity to polymorphic ventricular tachycardia and sudden cardiac death in young patients (Crotti et al, 2008). Among genotyped individuals, ∼90% of affected patients have mutations in the genes encoding repolarizing K+ channels of the delayed rectifier currents, _I_Ks (KCNQ1 in LQT1) and _I_Kr (KCNH2 in LQT2) (Morita et al, 2008). LQTS presents clinically with a broad range of phenotypes even among family members with identical mutations (Giudicessi and Ackerman, 2013). This heterogeneity is thought to result from the interplay of complex networks of direct and indirect causal factors, including epistatic effects of genetic background, resulting in incomplete penetrance and variable severity. LQTS was one of the first cardiac diseases recapitulated using hiPSCs as a model (Moretti et al, 2010; Itzhaki et al, 2011; Matsa et al, 2011; Yazawa et al, 2011; Davis et al, 2012; Egashira et al, 2012; Lahti et al, 2012). However, there is a growing realization that comparing patient cells to unrelated healthy control cells, no matter how they are matched, does not control for any effects of genetic modifier loci (Bellin et al, 2012; Cherry and Daley, 2013; Sinnecker et al, 2013). Furthermore, additional key factors may contribute to variability among individual hiPSC clones, including the number and location of viral integrations following reprogramming (Gonzalez et al, 2011), genetic heterogeneity in the starting somatic cell population (Gore et al, 2011), copy number variations (Panopoulos et al, 2011; Abyzov et al, 2012), and their epigenetic status (Nazor et al, 2012). For these reasons, it is preferable to compare cell lines in which the disease-causing genetic lesion of interest is the only modified variable. Gene correction in hiPSCs from patients with known genetic mutations is a powerful tool for overcoming the limitations imposed by individual variability among independent controls, and with recent technological improvements genome editing is now feasible in both human embryonic stem cells (hESCs) and hiPSCs (Costa et al, 2007; Zou et al, 2009; Hockemeyer et al, 2011; Asuri et al, 2012). The generation of genetically matched controls that differ exclusively in well-validated susceptibility variants for maladies is a generally applicable solution to the key problem of distinguishing pathologically relevant phenotypic changes from other background-related variations (Soldner et al, 2011; Reinhardt et al, 2013). However, to date, no isogenic cardiac disease iPSC lines have been reported, and neither has the disease phenotype been investigated in parallel in both hiPSCs and hESCs harbouring the same mutation.

In this study, we sought to generate and characterize a panel of control and LQT2-related human pluripotent stem cell (hPSC) lines. In order to investigate the specific effects of the c.A2987T (N996I) KCNH2 mutation without the confounding elements that could result from individual genetic background variability, we studied this mutation under two genetically defined conditions. First, by genetically correcting the KCNH2 mutation in the LQT2-hiPSCs derived from a patient, and second, by introducing the same c.A2987T (N996I) mutation into an NKX2.5 eGFP/w hESC reporter line that allows the cardiac cells to be selected from a mixed differentiated population (Elliott et al, 2011). Cardiomyocytes (CMs) differentiated from these two pairs of isogenic hPSCs demonstrated that the N996I KCNH2 mutation causes an ∼30–40% _I_Kr reduction with consequential action potential (AP) duration (APD) prolongation, and indicated a trafficking defect as the associated pathological mechanism responsible for the disease electrophysiological phenotype.

Results

Generation of hiPSCs with c.A2987T (N996I) KCNH2 mutation

Dermal fibroblasts were obtained from a 38-year-old woman with a diagnosis of familial type-2 LQTS (Supplementary Figure S1A and B). Genetic screening indicated a heterozygous c.A2987T mutation in exon 13 of the KCNH2 gene (Figure 1A), resulting in an asparagine to isoleucine substitution at position 996 (N996I) within the intracellular C-terminal region of the encoded HERG potassium channel (Figure 1B). We generated hiPSCs (LQT2-hiPSCsN996I) from the primary skin fibroblasts by retroviral transduction of four reprogramming genes (OCT4, SOX2, KLF4, and MYC) (Takahashi et al, 2007a) (Figure 1C), and three of the resulting clones underwent further characterization. These clones expressed the hESC markers SSEA4 and NANOG, as detected by immunofluorescence staining (Figure 1D). Additionally, qRT–PCR indicated that the hiPSC lines had reactivated the endogeneous pluripotency genes OCT4, SOX2, LEFTYA, NANOG, REX1, and TDGF1, and silenced the retroviral transgenes (Supplementary Figure S2A and B). Finally, differentiation of the hiPSCs as embryoid bodies demonstrated the capacity of all lines to generate derivatives of the three embryonic germ layers in vitro (Supplementary Figure S3). On the basis of these findings, one hiPSC clone was selected for targeted genetic correction.

Figure 1.

Figure 1

Generation of hiPSCs from a patient with type-2 long-QT syndrome. (A) Genetic screening in the patient revealed the heterozygous single-nucleotide mutation A→T in exon 13 of the KCNH2 gene, in position 2987 of the coding sequence (CDS) (c.A2987T, NM_000238.3), resulting in the substitution of an asparagine with an isoleucine at position 996 of the protein (N996I, NP_000229.1). (B) The N996I mutation (red dot) is located in the C-terminal of the HERG protein, which is made of six trans-membrane domains (S1–S6), an amino (NH2) domain, a carboxyl (COOH) domain, and a pore (P) region. (C) Example of a hiPSC colony harbouring the c.A2987T (N996I) KCNH2 mutation (LQT2-hiPSCsN996I). Scale bar: 400 μm. (D) Immunofluorescence analysis of pluripotency markers SSEA4 (green) and NANOG (red) in a representative LQT2-hiPSCN996I clone, with nuclear staining (DNA, blue). The image on the right is a magnification of the area framed in the left image. Scale bars: 100 μm (left image); 50 μm (right image).

Targeted gene correction in LQT2-hiPSCs N996I and targeted gene mutation in NKX2.5 eGFP/w hESCs

The study design is illustrated in Figure 2A. A conventional homologous recombination strategy was performed for either targeted gene correction in LQT2-hiPSCsN996I or targeted gene mutation in the NKX2.5 eGFP/w hESCs. We constructed two targeting vectors that included either a wild-type A or a mutated T nucleotide at KCNH2 c.2987 (KCNH2-A-loxP-pGK-Neo-loxP and KCNH2-T-loxP-pGK-Neo-loxP, respectively) (Figure 2B). Targeting of either the mutant allele (hiPSCs) or one of the two wild-type alleles (hESCs) is expected to result in a single base pair change that either repairs (LQT2-hiPSCscorr) or introduces (NKX2.5 eGFP/w hESCsN996I) the N996I HERG mutation. The coding regions and the exon–intron boundaries of KCNH2 in the hiPSCs and hESCs were sequenced to exclude the presence of any other variants.

Figure 2.

Gene targeting by homologous recombination in LQT2-hiPSCsN996I and in NKX2.5 eGFP/w hESCs. (A) Schematic showing the project rationale. Patient-specific LQT2-hiPSCs harbouring the N996I mutation were corrected and NKX2.5 eGFP/w hESCs were mutated by gene targeting. Parental and genetically modified hPSC lines were differentiated into CMs and their electrophysiological phenotypes were analysed and compared. (B) The strategy for precise genomic modification of KCNH2. Top line, structure of the KCNH2 locus. Numbered black boxes indicate exons 6–15. Exon 13 is mutated (red T) in LQT2-hiPSCsN996I and wild type (black A) in NKX2.5 eGFP/w hESCs. The gene targeting vector for correcting the mutation in LQT2-hiPSCs has the wild-type adenine nucleotide, whereas the gene targeting vector for introducing the mutation in NKX2.5 eGFP/w hESCs has the mutated thymine nucleotide. NeoR, the PGK-Neo cassette encoding G418 resistance flanked by loxP sequences (grey triangles), was inserted in the reverse direction. PCR primers (a, b) and (c, d) were used to identify the targeted clone. (C) PCR analysis using these primers generated specific bands of 4.8 kb (5′ homology arm) and 7.3 kb (3′ homology arm) from targeted clones (LQT2-hiPSCsNeoR and NKX2.5 eGFP/w hESCsNeoR). (D) Sequence analysis of PCR-amplified genomic DNA showing correction of the c.A2987T mutation in the LQT2-hiPSCscorr line and mutation in the NKX2.5 eGFP/w hESCsN996I line. The wild-type reference sequence (Ref) is shown in the top line.

Source data for this figure is available on the online supplementary information page.

Figure 2

The linearized targeting constructs were electroporated into the hiPSCs and hESCs and G418-resistant colonies isolated. Targeted clones were identified using a PCR-based screen that amplified the novel junction fragment between the genomic DNA and the integrated targeting vector generated following homologous recombination (Figure 2B and C). These PCR products were sequenced to confirm homologous recombination and to determine whether the mutation was corrected or introduced. We obtained one corrected LQT2-hiPSC clone (LQT2-hiPSCscorr) and one mutated NKX2.5 eGFP/w hESC clone (LQT2-hESCsN996I) (Figure 2D). The G418-resistance cassette was excised and the cell line cloned. Similar to the parental cell lines, the LQT2-hiPSCscorr and NKX2.5 eGFP/w hESCsN996I expressed markers of undifferentiated hPSCs (Supplementary Figure S4A and B), could differentiate into the cellular derivatives of the three embryonic germ layers in vitro (Supplementary Figure S5), passed the PluriTest with a high ‘pluripotency score’ and a low ‘novelty score’, indicating that they resemble normal hPSCs (Supplementary Figure S6A), and retained a normal karyotype (Supplementary Figure S6B and C).

The c.A2987T (N996I) KCNH2 mutation causes a long-QT phenotype in hPSC-derived CMs

The LQT2-hiPSCscorr and NKX2.5 eGFP/w hESCsN996I, together with their parental cell lines, were differentiated to CMs, with spontaneously contracting foci first observed at 10–12 days of differentiation. These were mechanically micro-dissected and allowed to mature further (Moretti et al, 2010). Dissociated cells displayed characteristic sarcomeric structures that were positive for troponin I (TNNI) and alpha-actinin (Figure 3A), and continued to spontaneously contract. Quantitative RT–PCR in cardiac cells, corresponding to either micro-dissected beating areas (hiPSCs) or NKX2.5 eGFP+ cell populations (hESCs), revealed upregulation of both the KCNH2-1a and -1b transcript variants, as well as other ion channels essential for AP generation, compared to their undifferentiated counterparts (Figure 3B).

Figure 3.

Figure 3

Differentiation of hPSCs into the cardiac lineage. (A) Confocal immunofluorescence images of cardiac sarcomeric proteins TNNI (green) and α-actinin (red) in human CMs generated from LQT2-hiPSCsN996I, LQT2-hiPSCscorr, NKX2.5 eGFP/w hESCs, and NKX2.5 eGFP/w hESCsN996I. Nuclei are stained in blue. Bottom panels are a magnification of the area framed in the upper images. Top panels, scale bar: 25 μm; bottom panels, scale bar: 10 μm. (B) Transcriptional profile of human CMs generated from mutated and corrected hiPSCs (LQT2-hiPSCsN996I and LQT2-hiPSCscorr, respectively) and from wild-type and mutated hESCs (NKX2.5 eGFP/w hESCs and NKX2.5 eGFP/w hESCsN996I, respectively). Undifferentiated cells from each hPSC line are also shown (Undiff.). Quantitative RT–PCR analysis was performed on the cardiac troponin gene (TNNT2), to show enrichment for the cardiomyocyte population, on the HERG channel gene (KCNH2), and on other key genes encoding for ion channels involved in the generation of the action potential in cardiac cells. All values are normalized to GAPDH and are relative to undifferentiated NKX2.5 eGFP/w hESCs. Raw minimum (min) and raw maximum (max) values were taken as a reference for heatmap representation.

To confirm that the loxP sequence that remained in the KCNH2 locus following excision of the selection cassette did not interfere with HERG expression, we took advantage of another LQT2-hiPSC clone that was targeted and contained the residual loxP sequence, but was not corrected (LQT2-hiPSCsN996I/loxP-control). Analysing the response of the CMs to the selective _I_Kr channel blocker E-4031 using multielectrode arrays revealed that, consistent with earlier reports (Matsa et al, 2011), increasing concentrations of E-4031 caused prolongation of field potential duration (FPD) in all analysed LQT2-hiPSC lines (original LQT2-hiPSCsN996I, LQT2-hiPSCsN996I/loxP-control, and LQT2-hiPSCscorr) (Supplementary Figure S7A). Importantly, the field potential response did not differ between LQT2-hiPSCN996I and LQT2-hiPSCN996I/loxP-control CMs and was prolonged when compared to LQT2-hiPSCcorr CMs, indicating not only that the N996I mutation conferred increased sensitivity to _I_Kr blockade but also that the remaining loxP sequence did not appear to interfere with their electrophysiological phenotype. We therefore used the LQT2-hiPSCsN996I/loxP-control line for further electrophysiological characterization. Moreover, similar analysis in CMs from NKX2.5 eGFP/w hESCs and NKX2.5 eGFP/w hESCsN996I showed that introduction of the N996I mutation in a genetically distinct background resulted in a significant, E-4031-dependent FPD prolongation that was comparable to that measured in mutated LQT2-hiPSC CMs, suggesting an exclusive effect of this mutation on the observed phenotype (Supplementary Figure S7A–D).

To further investigate whether the correction of the c.A2987T (N996I) mutation resulted in an increase in the repolarizing potassium current specifically transduced by the HERG channel, we recorded _I_Kr in individual hiPSC-derived CMs. Typical examples of current traces recorded before and after addition of 1 μM of E-4031, in mutated and corrected LQT2-hiPSC-derived CMs are shown in Figure 4A and B, where _I_Kr is defined as E-4031-sensitive current. Average _I_Kr density in LQT2-hiPSCN996I CMs was significantly smaller than in the LQT2-hiPSCcorr CMs (Figure 4C and D). At a membrane potential of −20 mV, both _I_Kr density (measured at the end of the test pulse) and peak-tail _I_Kr density were decreased by 33%. To confirm whether the c.A2987T (N996I) mutation is the cause of the LQT2 phenotype observed in hPSC-derived myocytes, we measured _I_Kr in CMs derived from NKX2.5 eGFP/w hESCs and NKX2.5 eGFP/w hESCsN996I that represent a different genetic background. Representative examples of current traces recorded in wild-type and mutated hESC-CMs are shown in Figure 5A and B. Average _I_Kr density was significantly reduced in the mutated hESC-CMs compared with the wild-type cells (Figure 5C and D). At a membrane potential of −20 mV, _I_Kr density (measured at the end of the test pulse) and peak-tail _I_Kr density were decreased by 43 and 45%, respectively, corroborating the results obtained in the patient-specific hiPSC model. Importantly, while _I_Kr densities were affected by the mutation, voltage dependence of activation parameters (half-maximal voltage (_V_1/2) and slope factor (k)) were similar, both in mutated and corrected LQT2-hiPSC-derived CMs and in wild-type and mutated hESC-derived CMs (Figure 6A and B). Furthermore, the time constant of _I_Kr activation (τ) (Figure 6C and D), as well as the fast (τf) and the slow (τs) deactivation time constants were all unchanged when comparing wild-type and mutated cells (Figure 6E and F). Taken together, these findings suggest that the N996I HERG mutation is the only relevant mutation that influences _I_Kr current in the heterozygous LQT2 CMs.

Figure 4.

Figure 4

_I_Kr densities in mutated and corrected LQT2-hiPSC-derived CMs. (A) Representative current traces elicited upon 4 s depolarizing voltage steps to −30, −20, and −10 mV from a holding potential of −40 mV, before and after the application of 1 μM E-4031. Inset: voltage protocol. (B) Typical tail currents measured after the depolarizing step to −10 mV in corrected (black) and mutated (red) LQT2-hiPSC-CMs, showing a bi-exponential decay. (C) Average current–voltage (I_–_V) relationships for _I_Kr, measured at the end of the test pulses, in mutated (red) and corrected (black) LQT2-hiPSC-derived CMs. Inset: voltage protocol. * indicates statistical significance (_P_=0.046, two-way rmANOVA; Holm–Sidak test post hoc analysis: −20 mV: _P_=0.015, −10 mV: _P_=0.007, 0 mV: _P_=0.018, 10 mV: _P_=0.011, 20 mV: _P_=0.013). (D) Average I_–_V relationships for peak tail currents in mutated (red) and corrected (black) LQT2-hiPSC-derived CMs. * indicates statistical significance (_P_=0.039, two-way rmANOVA; Holm–Sidak test post-hoc analysis: −20 mV: _P_=0.008, −10 mV: _P_=0.004, 0 mV: _P_=0.005, 10 mV: _P_=0.004, 20 mV: _P_=0.006). Inset: voltage protocol.

Figure 5.

Figure 5

_I_Kr densities in wild-type and mutated hESC-derived CMs. (A) Representative current traces elicited upon 4 s depolarizing voltage steps to −30, −20, and −10 mV from a holding potential of −40 mV, before and after the application of 1 μM E-4031. Inset: voltage protocol. (B) Typical tail currents measured after the depolarizing step to −10 mV in wild-type (black) and mutated (red) hESC-CMs, showing a bi-exponential decay. (C) Average I_–_V relationships for _I_Kr, measured at the end of the test pulses, in wild-type (black) and mutated (red) hESC-derived CMs. Inset: voltage protocol. * indicates statistical significance (_P_=0.026, two-way rmANOVA; Holm–SIdak test post hoc analysis: −30 mV: _P_=0.004; −20 mV: _P_=0.000, −10 mV: _P_=0.002, 0 mV: _P_=0.009). (D) Average I_–_V relationships for peak tail currents in wild-type (black) and mutated (red) hESC-derived CMs. * indicates statistical significance (_P_=0.014, two-way rmANOVA; Holm–Sidak test post-hoc analysis: −30 mV: _P_=0.022, −20 mV: _P_=0.001, −10 mV: _P_=0.002, 0 mV: _P_=0.001, 10 mV: _P_=0.001, 20 mV: _P_=0.003). Inset: voltage protocol.

Figure 6.

Figure 6

_I_Kr activation and deactivation properties in mutated and corrected LQT2-hiPSC-derived CMs and in wild-type and mutated hESC-derived CMs. (A, B) Average peak tail current normalized to the maximal current following repolarization to −40 mV in mutated (red) and corrected (black) LQT2-hiPSC-derived CMs (A) and in wild-type (black) and mutated (red) hESC-derived CMs (B). Inset: voltage protocol. (C, D) Average time constants of _I_Kr activation (τ) in mutated (red) and corrected (black) LQT2-hiPSC-derived CMs (C) and in wild-type (black) and mutated (red) hESC-derived CMs (D). (E, F) Average slow and fast time constants of _I_Kr deactivation (τs and τf, respectively) in mutated (red) and corrected (black) LQT2-hiPSC-derived CMs (E) and in wild-type (black) and mutated (red) hESC-derived CMs (F).

To evaluate the impact of the _I_Kr change on the APD, APs were recorded in single spontaneously contracting cells (Figure 7). Representative APs (paced at 1 Hz) from mutated and corrected hiPSC-derived CMs, and from wild-type and mutated hESC-derived CMs are shown in Figure 7A and C, respectively. The APD at 50 and 90% of repolarization (APD50 and APD90, respectively) was significantly reduced in the LQT2-hiPSCcorr CMs compared with their mutated counterparts (Figure 7B). In the hESC model, only the APD50 of NKX2.5 eGFP/w hESCN996I CMs was significantly prolonged compared with the NKX2.5 eGFP/w hESC CMs (Figure 7D). In contrast, in both the hiPSC and hESC paradigms, no differences in upstroke velocity (_V_max), AP amplitude (APA), and maximal diastolic potential (MDP) were detected between the mutated and corrected patient lines, and between the wild-type and mutated hESC groups (Figure 7B and D, respectively).

Figure 7.

Figure 7

AP characteristics in mutated and corrected LQT2-hiPSC-derived CMs and in wild-type and mutated hESC-derived CMs. (A, B) Representative examples of AP measured at 1 Hz (A) and average APD50, APD90, _V_max, APA, and MDP (B) in mutated (red) and corrected (black) LQT2-hiPSC-derived CMs. * indicates statistical significance (APD50: _P_=0.045, APD90: _P_=0.026; _t_-test). (C, D) Representative examples of AP measured at 1 Hz (C) and average APD50, APD90, _V_max, APA, and MDP (D) in wild-type (black) and mutated (red) hESC-derived CMs. * indicates statistical significance (APD50: _P_=0.049, _t_-test).

The c.A2987T (N996I) KCNH2 mutation causes a trafficking defect in the HERG channel

To investigate the functional consequences of the N996I-HERG mutation further, we examined the cellular distribution of the HERG channel in hiPSC- and hESC-derived CMs. Immunocytochemistry in both hPSC-CMs revealed a distribution of the HERG channel over the cell surface (Figure 8A), as well as in the endoplasmic reticulum (ER) and in the Golgi apparatus (Supplementary Figures S8 and S9). However, it was not possible to determine accurately by immunofluorescence whether CMs generated from the LQT2-hPSCsN996I models presented any differences in subcellular localization of HERG channels when compared to their wild-type or corrected counterparts. Therefore, we exploited the fact that we had introduced the c.A2987T (N996I) KCNH2 mutation into an NKX2.5 eGFP/w hESC line in which GFP reports the expression of the cardiac associated transcription factor NKX2.5 (Elliott et al, 2011). This allowed us to purify a population enriched for CMs by flow cytometry (Figure 8B), and quantify HERG channel levels within this fraction by western blot analysis. Two HERG protein bands were identified consistent with core glycosylated (135 kDa, immature) and more complex glycosylated (155 kDa, mature) isoforms (Figure 8C), with the latter present predominantly in the _NKX2.5_-eGFP+ fraction. In the _NKX2.5-_eGFP+ hESCN996I cells, the 155-kDa protein band, representing the form transported to the cell membrane through the Golgi, was reduced by ∼2-fold compared to the wild-type cells, while the 135-kDa band, which corresponds to the protein located in the ER, was unaffected (Figure 8D). This resulted in a significant decrease in HERG trafficking efficiency, calculated as a ratio of the fully glycosylated 155 kDa band over total HERG protein (155 kDa band+135 kDa band), of ∼40% (Supplementary Figure S10A). Similar analysis in LQT2-hiPSCcorr and LQT2-hiPSCN996I CMs confirmed a specific reduction in the complex glycosylated HERG protein in the mutated cells compared to the corrected counterparts (Figure 8E and F; Supplementary Figure S10B). Thus, the N996I mutation appears to interfere with the maturation and trafficking of the channel to the membrane.

Figure 8.

Trafficking defect in CMs harbouring the c.A2987T (N996I) KCNH2 mutation. (A) lmmunofluorescence images of HERG channel (green) and actin (red) in representative human CMs derived from mutated and corrected LQT2-hiPSCs, and from wild-type and mutated hESCs. Nuclei are stained in blue. Bottom panels are a magnification of the area framed in the upper corresponding images. Top panels, scale bar: 25 μm; bottom panels, scale bar: 5 μm. (B) Flow cytometry purification of NKX2.5 eGFP+ and NKX2.5 eGFP− hESC population from embryoid bodies differentiated from wild-type (left) and mutated (right) hESCs. The pictures show an individual eGFP-expressing (green) embryoid body; BF: bright field; scale bars: 400 μm. The representative dot plots show flow cytometric isolation of eGFP+ (green) and eGFP− (orange) cell populations. (C) Representative western blot analysis of HERG protein in eGFP+ and eGFP− cell populations purified from differentiated wild-type and mutated hESCs. Core- and complex-glycosylated HERG (135 and 155 kDa, respectively) are indicated. Actin is shown as a loading control. (D) Densitometric quantification of the 155-kDa and 135-kDa bands corresponding to the complex- and core-glycosylated HERG channel, respectively; ADU: arbitrary densitometric units; values are presented as mean±s.e.m., _n_=4. (E) Representative western blot analysis of HERG protein in LQT2-hiPSCscorr- and LQT2-hiPSCsN996I-derived CMs. Core- and complex-glycosylated HERG (135 and 155 kDa, respectively) are indicated. Actin is shown as a loading control. (F) Densitometric quantification of the 155-kDa and 135-kDa bands corresponding to the complex- and core-glycosylated HERG channel, respectively; ADU: arbitrary densitometric units; values are presented as mean±s.e.m., _n_=2.

Source data for this figure is available on the online supplementary information page.

Figure 8

To further investigate the pathological mechanism of the N996I mutation, we tested whether this mutation resulted in a misfolded HERG protein by analysing the unfolded protein response (UPR) in CMs purified from NKX2.5 eGFP/w hESCs and NKX2.5 eGFP/w hESCsN996I. UPR is an ER stress pathway that increases the synthesis of chaperones, such as Calnexin and Calreticulin, which in turn target misfolded proteins for proteasome degradation (Ellgaard and Helenius, 2003; Hetz, 2012). Western blot for the Activating Transcription Factor 6 (ATF6), a key regulator of transcriptional control in the mammalian UPR (Haze et al, 1999; Chen et al, 2002), revealed the presence of several forms of ATF6, including both the unprocessed protein embedded in the ER (90 kDa band) and some smaller forms resulting from the cleavage of the full-length ATF6 during ER stress, with no detectable differences between wild-type and mutated _NKX2.5_-eGFP+ cells (Figure 9A). Moreover, there was no evidence of an increase in the expression of the chaperones Calnexin and Calreticulin in NKX2.5 eGFP/w hESCN996I CMs (Figure 9B and C), indicating that the N996I HERG mutation does not appear to stimulate the UPR pathway. Nevertheless, to assess the involvement of proteasome-mediated and/or lysosome-mediated degradation of HERG channel further, we determined the effects of the proteasome inhibitor lactacystin and the lysosome inhibitor leupeptin on the steady-state HERG protein levels in wild-type and mutated hESC-CMs. Following purification of the _NKX2.5_-eGFP+ cell population by flow cytometry, cells were treated for 24 h with either lactacystin (20 μM) or leupeptin (100 μM) and cell lysates were then immunoblotted with anti-HERG antibody (Figure 9D). Treatment with the proteasome inhibitor lactacystin augmented the level of the mature, highly glycosylated form of HERG in NKX2.5 eGFP/w hESCN996I eGFP+ cells, resulting in an ∼1.4 fold increase in channel trafficking efficiency. No increase was observed in the wild-type counterparts or after treatment with the lysosome inhibitor leupeptin (Figure 9D and E), suggesting a role of proteasomes in the degradation of N996I HERG protein.

Figure 9.

UPR pathway analysis in wild-type and mutated hESC-CMs. (A–C) Western blot analysis of ATF6 (A), Calnexin (B), and Calreticulin (C) in eGFP+ and eGFP− cell populations purified from differentiated wild-type and mutated hESCs. Actin is shown as a loading control. (D, E) Western blot analysis of HERG channel, calnexin, and calreticulin (D) and quantification of trafficking efficiency (E) under basal conditions (CTR) and upon proteasome (LACT) or lysosome (LEUP) inhibition in eGFP+ cell populations purified from differentiated wild-type and mutated hESCs. Actin is shown as a loading control. Trafficking efficiency=fg/(fg+cg), where fg=fully-glycosylated 155 kDa band and cg=core-glycosylated 135 kDa band.

Source data for this figure is available on the online supplementary information page.

Figure 9

Discussion

The ability to reprogramme somatic cells from any individual of choice has created unprecedented opportunities to study disease development in vitro and create novel platforms for drug screening and discovery. However, one crucial limitation for disease modelling in hiPSCs has been that, unlike experimental animal models, humans are genetically heterogeneous (Goldstein, 2009; McKernan et al, 2009), and genetically matched controls are necessary to distinguish disease-relevant phenotypic changes from normal background-related variations due to polymorphisms in other genes. Here, we generated a panel of isogenic control and cardiac disease lines and investigated the disease phenotype in parallel in both hiPSCs and hESCs harbouring the same N996I KCNH2 mutation. Although this mutation has been described as a novel LQT2 variant and included in studies exploring possible genotype–phenotype correlations (Khositseth et al, 2004; Tester et al, 2005; Tan et al, 2006), its mechanism of action has remained unclear. Here, we demonstrated that this specific mutation is the cause of the LQT2 phenotype observed in the hPSC-derived CMs. Furthermore, we investigated the biophysical disease mechanism in functionally relevant human cardiac cells. Our findings demonstrated that the N996I KCNH2 mutation caused a similar (∼30–40%) _I_Kr reduction in both mutated hiPSC- and hESC-derived CMs. In particular, because a decrease of ≤50% in the cardiac repolarizing potassium current indicates haploinsufficiency (Moss et al, 2007; Vandenberg et al, 2012), our results suggested this as the mechanism underlying _I_Kr reduction, rather than a dominant-negative effect of mutant subunits on the wild-type subunits. Furthermore, the reduction in the fully glycosylated mature 155 kDa HERG protein in the mutated hPSC-CMs and its partial rescue after treatment with the proteasome inhibitor lactacystin suggest that the mutation might induce, to some extent, protein clearance through the proteasomes, thus resulting in defective channel trafficking (Vandenberg et al, 2012). Further analysis of HERG protein stability and internalization is required to unravel the additional mechanisms through which mutated N996I HERG channels may be processed. It is noteworthy that, in contrast to previous reports examining trafficking-deficient HERG mutants in heterologous systems (Smith et al, 2011; Guo et al, 2012; Wang et al, 2012), we did not observe any accumulation of the core-glycosylated form of HERG in the ER nor an activation of the UPR in the mutated hPSC-derived CMs. This might reflect either the specific effect of this particular mutation or the different cell systems used to investigate the pathophysiological mechanism. In overexpression systems, the levels of both wild-type and mutated HERG proteins are far from physiological; this could alter molecular stoichiometry and cause a more pronounced stimulation of ER stress and clearance pathways than in the native conditions of a CM.

It is widely recognized that hPSC-derived CMs are immature, with electrical properties, gene expression and contraction forces only equivalent to those of fetal or neonatal human myocytes (de Boer et al, 2010; Hoekstra et al, 2012). Nevertheless, our results show that these cells are capable of capturing specific traits of an electrical disease of the heart, such as LQT2, corroborating earlier findings (Itzhaki et al, 2011; Matsa et al, 2011). Of note, AP parameters (especially upstroke velocity) suggested that hESC-derived CMs were electrophysiologically less mature than those derived from the hiPSCs used here. While this did not preclude detection of _I_Kr reduction as a specific effect of the N996I-HERG mutation, the electrophysiological immaturity of the hESC-CMs might have influenced the impact of _I_Kr reduction on APD, where we observed a significant prolongation of APD50 but not of both APD50 and APD90 as observed in hiPSC-CMs. In this regard, the longer plateau phase of the AP in hiPSC-CMs compared to hESC-CMs would allow a higher number of HERG channels to inactivate, resulting in a different contribution of _I_Kr to AP repolarization (Rudy, 2008). These differences could also be due to inherent variability among individual hiPSC and hESC lines (Narsinh et al, 2011; Mummery et al, 2012) but do not detract from the conclusion since isogenic pairs were always compared. In addition, while the _I_Kr traces recorded in our lines are in agreement with _I_Kr previously reported in hPSC-derived CMs (Moretti et al, 2010; Itzhaki et al, 2011; Ma et al, 2011; Lahti et al, 2012), the current density measured in our hESC-CMs was higher (3.5 pA/pF) than those described in the literature for both hiPSC-derived (0.55–1.9 p/pF) and native human (0.25–0.6 pA/pF) CMs (Hoekstra et al, 2012). The molecular basis of this discrepancy requires further independent investigation, but a larger _I_Kr could impact on the AP characteristics, for example, contributing to a faster repolarization resulting in a shorter AP. Indeed, while the APD values measured in the LQT2-hiPSCcorr CMs fall into the range reported for hPSC-derived CMs (Hoekstra et al, 2012), in the NKX2.5 eGFP/w hESC CMs the APD is rather small.

Another aspect of our work is noteworthy: we measured very modest APD prolongation in both LQT2-hiPSCN996I and NKX2.5 eGFP/w hESCN996I CMs when compared to their wild-type counterparts, with values shorter than those reported in other LQT2 hiPSC models (Itzhaki et al, 2011; Matsa et al, 2011; Lahti et al, 2012), and did not observe any early-after depolarization in the mutated cells. These differences are likely due to the specific N996I mutation, which results in a less severe LQT2 phenotype than those caused by other trafficking-deficient mutations (Moss et al, 2002). However, it should be noted that in all previous studies using LQT2 hiPSC models the control group consisted of wild-type hiPSCs derived from a healthy individual unrelated to the LQT2 patient, therefore over-estimation of differences cannot be excluded under these circumstances. The advantage of using isogenic lines is indeed highlighted by our results, where even subtle differences could be appreciated.

Most importantly, our data provide evidence that the pathogenesis of the N996I-HERG mutation can be modelled exclusively in CMs generated from hPSCs, without the need for heterologous and/or overexpression systems. In particular, we detected and analysed physiological levels of HERG protein (by immunostaining and western blot), as well as molecular chaperones and transcriptional regulators of the UPR, in functionally relevant human cardiac cells. Moreover, the use of an NKX2.5 eGFP/w hESC reporter line allowed selection of the cardiac cell population from other cell types, such as neurons, endocrine cells, and undifferentiated proliferating cells, that also express KCNH2 (Rosati et al, 2000; Pillozzi et al, 2002; Huffaker et al, 2009; Vandenberg et al, 2012). This highlights the particular value of combining patient-specific iPSC technology and genetic modification tools.

Since our study compared genetically matched lines differing uniquely in one LQT2 point mutation, a clear genotype–phenotype relationship could be ascertained. This is of particular value in avoiding over- or under-estimation of the impact of disease-related changes as a result of inadvertent bias when choosing a random control (even if age and gender matched with the same reprogramming and differentiation method). Therefore, this approach is a powerful means of identifying links between specific genotypes and cardiac disease predisposition, as also recently demonstrated for neurodegenerative disease (Reinhardt et al, 2013).

Finally, targeted gene correction of KCNH2 in hiPSCs rescued _I_Kr density in hiPSC-CMs and normalized AP values. Although precise gene modification of human CMs in situ is not an immediate option for regenerative medicine or for treating LQTS, gene targeting as carried out here provides essential tools for deciphering molecular mechanisms, studying genetic–phenotypic interaction and perhaps, in the future, analysing the potential disease-causing SNPs identified by genome-wide association studies in humans. Other rescue approaches, such as RNA interference, may only result in partial knock down of the mutated mRNA (Matsa et al, 2013), while genetic correction should restore the normal genetic phenotype in monogenic diseases. This can then demonstrate whether the specific genetic lesion is pathogenic or whether there are likely to be other modifier loci involved. Considering the increasing awareness of the broad influence of genetic background (Montgomery et al, 2010) and the biological differences between hPSC lines (Boulting et al, 2011), our experimental system may overcome some of the shortcomings of conventional hiPSC approaches in identifying disease-related phenotypes. As more LQTS-hiPSC lines are investigated using this approach, the molecular and genetic reasons for variability of expression and for incomplete penetrance seen in this, and other channelopathies, may be elucidated.

In summary, we have demonstrated that genetic correction of the N996I KCNH2 mutation associated with LQT2 restores _I_Kr density and normalizes APD in patient-specific LQT2-hiPSC-CMs. The same mutation introduced into a different genetic background (hESC) results in a very similar electrophysiological phenotype. Furthermore, the reduction in the mature complex-glycosylated 155 kDa protein band in the mutated hPSC-CMs suggests that, due to the mutation, a smaller number of channels are present at the cell membrane. This could be partially due to proteasome targeting for degradation, but might not be the only mechanism for this N996I-specific mutation. Further studies are required to assess more in detail the exact contribution of proteasome clearance pathway, protein channel stability, and lysosome degradation. Our study demonstrates that isogenic pairs of hPSCs can be used (1) to prove the authenticity of the genotype–phenotype correlation and (2) to unravel the pathophysiological mechanism seen in a genetically inherited cardiac disease. As such, this approach represents a robust strategy to study not only cardiac disease mechanisms but also other genetic disorders.

Materials and methods

Clinical history and genetic phenotype

A 38-year-old Caucasian woman presented with a prolonged QT interval (QT interval corrected for heart rate (QTc) 617 ms) at the ECG and was diagnosed with type-2 LQTS. Genetic screening indicated a heterozygous mutation c.A2987T in exon 13 of the KCNH2 gene, resulting in the N996I missense mutation. The patient has thus far been asymptomatic. She is now treated with β-blockers.

Generation of patient-specific hiPSCs, differentiation into CMs, and drug treatment

We recruited a 38-year-old Caucasian female LQT2 patient for dermal biopsy after obtaining written informed consent. Reprogramming of primary skin fibroblasts was performed as described previously (Moretti et al, 2010; Jung et al, 2012). Briefly, fibroblasts were infected with retroviruses encoding OCT4, SOX2, KLF4, and MYC and cultured on mouse embryonic feeder (MEF) cells until hiPSC colonies could be picked. Both hiPSCs and hESCs were maintained in culture using standard procedures (Takahashi et al, 2007a) and were differentiated using standard cardiogenic protocols (Mummery et al, 2003; Elliott et al, 2011). Spontaneously contracting areas were manually dissected and cultured further until days 20–60 of differentiation. Cells for physiological experiments were enzymatically dissociated into single cells (Moretti et al, 2010), plated on fibronectin-coated glass coverslips, and analysed within 7–15 days.

For drug treatment, _NKX2.5_-eGFP+ cell populations were purified by flow cytometry and maintained in low attachment plates while treated for 24 h with 20 μM lactacystin (Sigma) or 100 μM leupeptin (Sigma); cells were then lysated and subjected to immunoblotting.

Genomic sequencing

Genomic DNA was isolated from cultured cells using the Gentra PureGene Cell Kit (Qiagen). To confirm the presence or absence of the c.A2987T KCNH2 mutation in the hPSCs, and to exclude the presence of other variants in the KCNH2 locus, the coding region and the exon–intron boundaries were PCR amplified from genomic DNA (Supplementary Table S1). The PCR products were purified using the QIAquick PCR Purification kit (Qiagen) and sequenced.

Generation and identification of targeted hiPSCs and hESCs

A BAC (RP11-10L20) containing the ∼33-kb human KCNH2 locus was modified stepwise by recombineering to generate the final targeting vector (Fu et al, 2010). First, a _loxP_-flanked positive selection cassette (loxP-PGK-neo-loxP) comprising a mammalian promoter (pGK), a bacterial promoter (gb3), and a G418/kanamycin-resistance gene (neo) was amplified from an R6K plasmid and inserted 0.9 kb downstream of the mutation site. A 15.6-kb fragment, including a 6.6-kb 5′-homology arm, the inserted G418-resistance cassette, and a 4.3-kb 3′-homology arm, was then sub-cloned into a minimal vector (Gene Bridges), generating the wild-type KCNH2 targeting vector (KCNH2-A-loxP-pGK-Neo-loxP) used to correct the mutation in the LQT2-hiPSCsN996I. For hESC targeting, the c.A2987T mutation was introduced into this targeting vector using the QuickChange XL Site-Directed Mutagenesis Kit (Stratagene) (KCNH2-T-loxP-pGK-Neo-loxP). The two vectors, wild type and mutated, were linearized with the restriction enzyme _Not_I before electroporation (Costa et al, 2007). Targeted clones were identified using a PCR-based screening strategy. Homologous integration of both the 5′ and 3′ homology arms, as well as correction or introduction of the c.A2987T KCNH2 mutation in the resulting clones were confirmed by sequencing. The _loxP_-flanked G418-resistance cassette was excised using Cre recombinase from both a corrected- and non-corrected hiPSC clone, and a mutated hESC clone, and the resulting lines subcloned as described previously (Davis et al, 2008). Clonal lines were screened for the loss of the neomycin-resistance cassette and for the absence of the Cre recombinase expression plasmid by PCR (Supplementary Table S1). These PCR products were also sequenced to confirm excision of the neomycin-resistance cassette. Karyotype analysis was performed using COBRA-FISH as described elsewhere (Szuhai and Tanke, 2006). Twenty metaphase spreads for each sample were analysed.

Immunofluorescence and western blot analysis

Cells were fixed in 4% paraformaldehyde, permeabilised with phosphate buffer saline (PBS)/0.1% Triton X-100 (Sigma-Aldrich), and blocked with 10% FCS (Life Technologies). Samples were incubated overnight at 4°C with primary antibodies specific for the following: NANOG (goat polyclonal, R&D Systems), SSEA4 and TNNI (mouse monoclonal and rabbit polyclonal, respectively, both from Santa Cruz Biotechnology), α-ACTININ and HERG (mouse monoclonal and rabbit polyclonal, respectively, both from Sigma-Aldrich), PDI (mouse monoclonal, Abcam), and Golgin-97 (mouse monoclonal, Molecular Probes). Primary antibodies were detected with either Cy3- or Alexa-Fluor 488-conjugated antibodies. Nuclei were visualized with DAPI (Invitrogen) and F-actin with Phalloidin-Alexa-Fluor-594-conjugate (Invitrogen). Images were captured using either a Leica DMI6000-AF6000 fluorescence microscope or a Leica SP5 confocal laser-scanning microscope (both from Leica Microsystems).

Western blotting on whole-cell lysate of either eGFP+ and eGFP− cells purified from differentiated NKX2.5 eGFP/w hESCs and NKX2.5 eGFP/w hESCsN996I or of manually microdissected beating areas from differentiated LQT2-hiPSCsN996I and LQT2-hiPSCscorr was performed with standard protocols (Kurien and Scofield, 2006) using 40–80 μg proteins and the following primary antibodies: HERG (mouse monoclonal, Enzo Life Sciences), ACTIN (mouse monoclonal, Millipore), ATF6 (mouse monoclonal, ActiveMotif), Calnexin and Calreticulin (both rabbit polyclonal, Abcam). Western blots were quantified using the ImageJ software ( http://rsb.info.nih.gov/ij/). Whole film images are provided in Supplementary data.

Gene expression analysis

For qRT–PCR, total RNA was purified using the RNeasy Mini Kit (Qiagen) according to the manufacturer’s protocol. RNA samples were reverse transcribed using the SuperScript II First-Strand Synthesis kit (Invitrogen). Gene expression was assessed by qRT–PCR as described previously (Moretti et al, 2010). Gene expression levels were normalized to GAPDH and represented in a heatmap using the GENE-E software ( http://www.broadinstitute.org/cancer/software/GENE-E/index.html). Primer sequences are provided in Supplementary Table S1.

For whole-genome microarray, total RNA was isolated from undifferentiated hPSCs using the NucleoSpin miRNA kit (Macherey-Nagel) and hybridized on Illumina HT12v4 microarrays. The raw microarray data were analysed with the PluriTest algorithm ( http://www.pluritest.org; Muller et al, 2011).

Electrophysiological characterization

Data acquisition and analysis

Recordings were performed on single CMs 7–15 days after cell dissociation. For _I_Kr measurements, CMs were identified on the basis of their typical morphology, while for AP measurements, spontaneously contracting cells that could be paced at the stimulation frequency of 1 Hz were selected. Data were collected from at least three independent differentiations per line. _I_Kr and APs were recorded with the ruptured and the perforated patch-clamp technique, respectively, using an Axopatch 200B amplifier (Molecular Devices). Voltage control and data acquisition of _I_Kr and APs were performed with pClamp10.1 (Axon Instruments). Membrane currents were analysed with Clampfit 10.1 (Axon Instruments), while analysis of APs was performed with the custom-made software. Patch pipettes had a tip resistance of ∼2 MΩ and series resistance (Rs) was compensated by 80%. _I_Kr and APs were filtered (2 and 5 kHz, respectively), and digitized (2 and 40 kHz, respectively). Potentials were corrected for liquid junction potentials (Barry and Lynch, 1991), which were calculated using Clampex 10.1 (Axon Instruments) and were equal to 15 and 17 mV for AP and _I_Kr measurements, respectively.

Voltage-clamp experiments

_I_Kr was recorded at 37°C using a pipette solution containing (mM): K-gluconate 125, KCl 20, K2-ATP 5, HEPES 10, EGTA 10; pH 7.2 (KOH). The bath solution was Tyrode’s solution containing (mM): NaCl 140, KCl 5.4, CaCl2 1.8, MgCl2 1.0, glucose 5.5, HEPES 5; pH 7.4 (NaOH). Nifedipine (5 μM) and JNJ-303 (1 μM) were added to the external solution to block the L-type calcium current (_I_CaL) and the slow component of the rectifier potassium current (_I_Ks), respectively. Currents were elicited using 4 s hyper- and depolarizing pulses from a holding potential of −40 mV, as indicated in the voltage protocols depicted in Figures 4 and 5. The cycle length of the voltage clamp steps was 10 s. _I_Kr was measured as a 1-μM E-4031-sensitive current, by subtraction of the current recorded before and after E-4031 application. _I_Kr densities were calculated by dividing current amplitude (pA), measured either at the end of the test pulses (Figures 4C and 5C) or at the peak of the tail current (Figures 4D and 5D), by cell membrane capacitance (pF). Cell membrane capacitance was measured by dividing the decay time constant of the capacitive transient in response to 5 mV depolarizing steps from −40 mV, by the Rs. Average cell capacitance for both hiPSC- and hESC-CM groups are shown in Supplementary Figure S11.

Activation curves, determined from peak tail current normalized to the maximal value and plotted against the test pulse voltage, were fitted with Boltzmann equation (_y_=[1+exp{(_V_−_V_1/2)/k}]−1), where _V_1/2 is the half-maximal voltage of activation, and k is the slope factor.

Time constants of activation (τ) were obtained by fitting _I_Kr during depolarizing pulses with a single-exponential function (_y_=_A_0+_A_1exp[−t/τ]). Finally, the slow (τs) and fast (τf) time constants of deactivation were determined by fitting _I_Kr tail currents, at voltage steps positive to −10 mV, with a bi-exponential function (_y_=_A_0+_A_fexp[−t/τf]+_A_sexp[−t/τs]).

Current-clamp experiments

APs were recorded at 37°C in Tyrode’s solution. The pipette solution contained (mM): K-gluconate 125, KCl 20, NaCl 5, amphotericin-B 0.22, HEPES 10; pH 7.2 (KOH). APs were elicited at the stimulation frequency of 1 Hz by 3 ms, 1.2 × threshold current pulses through the patch pipette. MDP, _V_max, APA, and APD50 and APD90 were analysed. Data from 10 consecutive APs were averaged.

MEA electrophysiology

MEA standard measurements were performed in DMEM culture medium supplemented with 2% FCS as previously described (Braam et al, 2013). E-4031 was dissolved in DMSO at 10 mM and serial dilutions were made in culture medium.

Statistical analysis

The SigmaStat 3.5 software was used for statistical analysis. Results are expressed as mean±s.e.m. The software performed tests for normality and equality of variance on the data sets prior to application of a statistical test. Group comparisons were made using unpaired Student’s _t_-test when both conditions were fulfilled or Mann–Whitney test when the conditions of normality and equal variance were not met. For repetitive measurements, two-way rmANOVA followed by a Holm–Sidak test for post-hoc analysis was used. P<0.05 was considered as statistically significant.

Ethics statement

The protocols for research involving human subjects and for stem-cell research were approved by the institutional review board and the committee charged with oversight of embryonic stem-cell research at the Technical University of Munich. Written informed consent was received from the participant prior to inclusion in the study. The studies on hESCs were performed exclusively in the Netherlands where the ethical committee for research at the LUMC provided approval for their use.

Supplementary Material

Supplementary Information

Review Process File

Acknowledgments

We thank C Höhnke (Technical University of Munich) for obtaining the skin biopsy; K Takashi and S Yamanaka (Kyoto University) for providing viral vectors through Addgene; F Stewart (Technical University of Dresden) for providing the R6K plasmid; MJM van der Burg, K Szuhai, and H Tanke (Leiden University Medical Center) for karyotyping analysis; AO Verkerk (University of Amsterdam), M Birket (Leiden University Medical Center), D Sinnecker (Technical University of Munich) and R Dirschinger (Technical University of Munich) for comments on the manuscript; S Böhringer (Leiden University Medical Center) for statistical assistance. We are especially thankful to the type-2 long-QT syndrome patient for tissue donation. This work was supported by grants from the European Research Council, ERC 261053 (K-LL) and ERC 323182 (CLM); the German Research Foundation, Research Unit 923, Mo 2217/1-1 (AM), La 1238 3-1/4-1 (K-LL), and WE 1639/3-1 (AW); German Centre for Cardiovascular Research (K-LL and AM); ZoNMW Animal Alternatives, 114000101 (CLM, SC, and LT); Netherlands Proteomics Consortium, 050-040-250 (CLM); Netherlands Institute of Regenerative Medicine (CLM and RD); EU Marie Curie FP7-people-2011-IEF programme, HPSCLQT 29999 (MB) and a short-term EMBO fellowship, ASTF 387.00-2011 (CD’A).

Author contributions: MB, AM, CLM, and K-LL conceived the project and designed experimental details. MB performed hiPSC reprogramming together with CBJ. SC performed the single-cell electrophysiology measurements together with AW. LGJT did the MEA experiments. The targeting strategy was designed by RPD and MB, and experimental targeting work was done by MB and DW. DE provided the _NKX2.5_-hESC line. CD’A performed the FACS sort and western blot experiments. JH did the HERG immunofluorescence staining. The manuscript was written by MB, CLM, and AM.

Footnotes

Christine Mummery is co-founder and advisor of Pluriomics. The remaining authors declare that they have no conflict of interest.

References

  1. Abyzov A, Mariani J, Palejev D, Zhang Y, Haney MS, Tomasini L, Ferrandino AF, Rosenberg Belmaker LA, Szekely A, Wilson M, Kocabas A, Calixto NE, Grigorenko EL, Huttner A, Chawarska K, Weissman S, Urban AE, Gerstein M, Vaccarino FM (2012) Somatic copy number mosaicism in human skin revealed by induced pluripotent stem cells. Nature 492: 438–442 [DOI] [PMC free article] [PubMed] [Google Scholar]
  2. Asuri P, Bartel MA, Vazin T, Jang JH, Wong TB, Schaffer DV (2012) Directed evolution of adeno-associated virus for enhanced gene delivery and gene targeting in human pluripotent stem cells. Mol Ther 20: 329–338 [DOI] [PMC free article] [PubMed] [Google Scholar]
  3. Axelrod FB, Liebes L, Gold-Von Simson G, Mendoza S, Mull J, Leyne M, Norcliffe-Kaufmann L, Kaufmann H, Slaugenhaupt SA (2011) Kinetin improves IKBKAP mRNA splicing in patients with familial dysautonomia. Pediatr Res 70: 480–483 [DOI] [PMC free article] [PubMed] [Google Scholar]
  4. Barry PH, Lynch JW (1991) Liquid junction potentials and small cell effects in patch-clamp analysis. J Membr Biol 121: 101–117 [DOI] [PubMed] [Google Scholar]
  5. Bellin M, Marchetto MC, Gage FH, Mummery CL (2012) Induced pluripotent stem cells: the new patient? Nat Rev Mol Cell Biol 13: 713–726 [DOI] [PubMed] [Google Scholar]
  6. Boulting GL, Kiskinis E, Croft GF, Amoroso MW, Oakley DH, Wainger BJ, Williams DJ, Kahler DJ, Yamaki M, Davidow L, Rodolfa CT, Dimos JT, Mikkilineni S, MacDermott AB, Woolf CJ, Henderson CE, Wichterle H, Eggan K (2011) A functionally characterized test set of human induced pluripotent stem cells. Nat Biotechnol 29: 279–286 [DOI] [PMC free article] [PubMed] [Google Scholar]
  7. Braam SR, Tertoolen L, Casini S, Matsa E, Lu HR, Teisman A, Passier R, Denning C, Gallacher DJ, Towart R, Mummery CL (2013) Repolarization reserve determines drug responses in human pluripotent stem cell derived cardiomyocytes. Stem Cell Res 10: 48–56 [DOI] [PubMed] [Google Scholar]
  8. Chen X, Shen J, Prywes R (2002) The luminal domain of ATF6 senses endoplasmic reticulum (ER) stress and causes translocation of ATF6 from the ER to the Golgi. J Biol Chem 277: 13045–13052 [DOI] [PubMed] [Google Scholar]
  9. Cherry AB, Daley GQ (2013) Reprogrammed cells for disease modeling and regenerative medicine. Annu Rev Med 64: 277–290 [DOI] [PMC free article] [PubMed] [Google Scholar]
  10. Costa M, Dottori M, Sourris K, Jamshidi P, Hatzistavrou T, Davis R, Azzola L, Jackson S, Lim SM, Pera M, Elefanty AG, Stanley EG (2007) A method for genetic modification of human embryonic stem cells using electroporation. Nat Protoc 2: 792–796 [DOI] [PubMed] [Google Scholar]
  11. Crotti L, Celano G, Dagradi F, Schwartz PJ (2008) Congenital long QT syndrome. Orphanet J Rare Dis 3: 18. [DOI] [PMC free article] [PubMed] [Google Scholar]
  12. Davis RP, Casini S, van den Berg CW, Hoekstra M, Remme CA, Dambrot C, Salvatori D, Oostwaard DW, Wilde AA, Bezzina CR, Verkerk AO, Freund C, Mummery CL (2012) Cardiomyocytes derived from pluripotent stem cells recapitulate electrophysiological characteristics of an overlap syndrome of cardiac sodium channel disease. Circulation 125: 3079–3091 [DOI] [PubMed] [Google Scholar]
  13. Davis RP, Costa M, Grandela C, Holland AM, Hatzistavrou T, Micallef SJ, Li X, Goulburn AL, Azzola L, Elefanty AG, Stanley EG (2008) A protocol for removal of antibiotic resistance cassettes from human embryonic stem cells genetically modified by homologous recombination or transgenesis. Nat Protoc 3: 1550–1558 [DOI] [PubMed] [Google Scholar]
  14. Davis RP, van den Berg CW, Casini S, Braam SR, Mummery CL (2011) Pluripotent stem cell models of cardiac disease and their implication for drug discovery and development. Trends Mol Med 17: 475–484 [DOI] [PubMed] [Google Scholar]
  15. de Boer TP, van Veen TA, Jonsson MK, Kok BG, Metz CH, Sluijter JP, Doevendans PA, de Bakker JM, Goumans MJ, van der Heyden MA (2010) Human cardiomyocyte progenitor cell-derived cardiomyocytes display a maturated electrical phenotype. J Mol Cell Cardiol 48: 254–260 [DOI] [PubMed] [Google Scholar]
  16. Egashira T, Yuasa S, Suzuki T, Aizawa Y, Yamakawa H, Matsuhashi T, Ohno Y, Tohyama S, Okata S, Seki T, Kuroda Y, Yae K, Hashimoto H, Tanaka T, Hattori F, Sato T, Miyoshi S, Takatsuki S, Murata M, Kurokawa J et al. (2012) Disease characterization using LQTS-specific induced pluripotent stem cells. Cardiovasc Res 95: 419–429 [DOI] [PubMed] [Google Scholar]
  17. Ellgaard L, Helenius A (2003) Quality control in the endoplasmic reticulum. Nat Rev Mol Cell Biol 4: 181–191 [DOI] [PubMed] [Google Scholar]
  18. Elliott DA, Braam SR, Koutsis K, Ng ES, Jenny R, Lagerqvist EL, Biben C, Hatzistavrou T, Hirst CE, Yu QC, Skelton RJ, Ward-van Oostwaard D, Lim SM, Khammy O, Li X, Hawes SM, Davis RP, Goulburn AL, Passier R, Prall OW et al. (2011) NKX2-5(eGFP/w) hESCs for isolation of human cardiac progenitors and cardiomyocytes. Nat Methods 8: 1037–1040 [DOI] [PubMed] [Google Scholar]
  19. Fu J, Teucher M, Anastassiadis K, Skarnes W, Stewart AF (2010) A recombineering pipeline to make conditional targeting constructs. Methods Enzymol 477: 125–144 [DOI] [PubMed] [Google Scholar]
  20. Giudicessi JR, Ackerman MJ (2013) Determinants of incomplete penetrance and variable expressivity in heritable cardiac arrhythmia syndromes. Transl Res 161: 1–14 [DOI] [PMC free article] [PubMed] [Google Scholar]
  21. Gold-von Simson G, Goldberg JD, Rolnitzky LM, Mull J, Leyne M, Voustianiouk A, Slaugenhaupt SA, Axelrod FB (2009) Kinetin in familial dysautonomia carriers: implications for a new therapeutic strategy targeting mRNA splicing. Pediatr Res 65: 341–346 [DOI] [PubMed] [Google Scholar]
  22. Goldstein DB (2009) Common genetic variation and human traits. N Engl J Med 360: 1696–1698 [DOI] [PubMed] [Google Scholar]
  23. Gonzalez F, Boue S, Izpisua Belmonte JC (2011) Methods for making induced pluripotent stem cells: reprogramming a la carte. Nat Rev Genet 12: 231–242 [DOI] [PubMed] [Google Scholar]
  24. Gore A, Li Z, Fung HL, Young JE, Agarwal S, Antosiewicz-Bourget J, Canto I, Giorgetti A, Israel MA, Kiskinis E, Lee JH, Loh YH, Manos PD, Montserrat N, Panopoulos AD, Ruiz S, Wilbert ML, Yu J, Kirkness EF, Izpisua Belmonte JC et al. (2011) Somatic coding mutations in human induced pluripotent stem cells. Nature 471: 63–67 [DOI] [PMC free article] [PubMed] [Google Scholar]
  25. Guo J, Zhang X, Hu Z, Zhuang Z, Zhu Z, Chen Z, Chen W, Zhao Z, Zhang C, Zhang Z (2012) A422T mutation in HERG potassium channel retained in ER is rescurable by pharmacologic or molecular chaperones. Biochem Biophys Res Commun 422: 305–310 [DOI] [PubMed] [Google Scholar]
  26. Haze K, Yoshida H, Yanagi H, Yura T, Mori K (1999) Mammalian transcription factor ATF6 is synthesized as a transmembrane protein and activated by proteolysis in response to endoplasmic reticulum stress. Mol Biol Cell 10: 3787–3799 [DOI] [PMC free article] [PubMed] [Google Scholar]
  27. Hetz C (2012) The unfolded protein response: controlling cell fate decisions under ER stress and beyond. Nat Rev Mol Cell Biol 13: 89–102 [DOI] [PubMed] [Google Scholar]
  28. Hockemeyer D, Wang H, Kiani S, Lai CS, Gao Q, Cassady JP, Cost GJ, Zhang L, Santiago Y, Miller JC, Zeitler B, Cherone JM, Meng X, Hinkley SJ, Rebar EJ, Gregory PD, Urnov FD, Jaenisch R (2011) Genetic engineering of human pluripotent cells using TALE nucleases. Nat Biotechnol 29: 731–734 [DOI] [PMC free article] [PubMed] [Google Scholar]
  29. Hoekstra M, Mummery CL, Wilde AA, Bezzina CR, Verkerk AO (2012) Induced pluripotent stem cell derived cardiomyocytes as models for cardiac arrhythmias. Front Physiol 3: 346. [DOI] [PMC free article] [PubMed] [Google Scholar]
  30. Huffaker SJ, Chen J, Nicodemus KK, Sambataro F, Yang F, Mattay V, Lipska BK, Hyde TM, Song J, Rujescu D, Giegling I, Mayilyan K, Proust MJ, Soghoyan A, Caforio G, Callicott JH, Bertolino A, Meyer-Lindenberg A, Chang J, Ji Y et al. (2009) A primate-specific, brain isoform of KCNH2 affects cortical physiology, cognition, neuronal repolarization and risk of schizophrenia. Nat Med 15: 509–518 [DOI] [PMC free article] [PubMed] [Google Scholar]
  31. Itzhaki I, Maizels L, Huber I, Zwi-Dantsis L, Caspi O, Winterstern A, Feldman O, Gepstein A, Arbel G, Hammerman H, Boulos M, Gepstein L (2011) Modelling the long QT syndrome with induced pluripotent stem cells. Nature 471: 225–229 [DOI] [PubMed] [Google Scholar]
  32. Jung CB, Moretti A, Mederos y Schnitzler M, Iop L, Storch U, Bellin M, Dorn T, Ruppenthal S, Pfeiffer S, Goedel A, Dirschinger RJ, Seyfarth M, Lam JT, Sinnecker D, Gudermann T, Lipp P, Laugwitz KL (2012) Dantrolene rescues arrhythmogenic RYR2 defect in a patient-specific stem cell model of catecholaminergic polymorphic ventricular tachycardia. EMBO Mol Med 4: 180–191 [DOI] [PMC free article] [PubMed] [Google Scholar]
  33. Khositseth A, Tester DJ, Will ML, Bell CM, Ackerman MJ (2004) Identification of a common genetic substrate underlying postpartum cardiac events in congenital long QT syndrome. Heart Rhythm 1: 60–64 [DOI] [PubMed] [Google Scholar]
  34. Kurien BT, Scofield RH (2006) Western blotting. Methods 38: 283–293 [DOI] [PubMed] [Google Scholar]
  35. Lahti AL, Kujala VJ, Chapman H, Koivisto AP, Pekkanen-Mattila M, Kerkela E, Hyttinen J, Kontula K, Swan H, Conklin BR, Yamanaka S, Silvennoinen O, Aalto-Setala K (2012) Model for long QT syndrome type 2 using human iPS cells demonstrates arrhythmogenic characteristics in cell culture. Dis Model Mech 5: 220–230 [DOI] [PMC free article] [PubMed] [Google Scholar]
  36. Lee G, Papapetrou EP, Kim H, Chambers SM, Tomishima MJ, Fasano CA, Ganat YM, Menon J, Shimizu F, Viale A, Tabar V, Sadelain M, Studer L (2009) Modelling pathogenesis and treatment of familial dysautonomia using patient-specific iPSCs. Nature 461: 402–406 [DOI] [PMC free article] [PubMed] [Google Scholar]
  37. Liang P, Lan F, Lee AS, Gong T, Sanchez-Freire V, Wang Y, Diecke S, Sallam K, Knowles JW, Wang PJ, Nguyen PK, Bers DM, Robbins RC, Wu JC (2013) Drug screening using a library of human induced pluripotent stem cell-derived cardiomyocytes reveals disease-specific patterns of cardiotoxicity. Circulation 127: 1677–1691 [DOI] [PMC free article] [PubMed] [Google Scholar]
  38. Ma J, Guo L, Fiene SJ, Anson BD, Thomson JA, Kamp TJ, Kolaja KL, Swanson BJ, January CT (2011) High purity human-induced pluripotent stem cell-derived cardiomyocytes: electrophysiological properties of action potentials and ionic currents. Am J Physiol Heart Circ Physiol 301: H2006–H2017 [DOI] [PMC free article] [PubMed] [Google Scholar]
  39. Marchetto MC, Carromeu C, Acab A, Yu D, Yeo GW, Mu Y, Chen G, Gage FH, Muotri AR (2010) A model for neural development and treatment of Rett syndrome using human induced pluripotent stem cells. Cell 143: 527–539 [DOI] [PMC free article] [PubMed] [Google Scholar]
  40. Matsa E, Dixon JE, Medway C, Georgiou O, Patel MJ, Morgan K, Kemp PJ, Staniforth A, Mellor I, Denning C (2013) Allele-specific RNA interference rescues the long-QT syndrome phenotype in human-induced pluripotency stem cell cardiomyocytes. Eur Heart J (advance online publication, 6 March 2013; doi:; DOI: 10.1093/eurheartj/eht067) [DOI] [PMC free article] [PubMed] [Google Scholar]
  41. Matsa E, Rajamohan D, Dick E, Young L, Mellor I, Staniforth A, Denning C (2011) Drug evaluation in cardiomyocytes derived from human induced pluripotent stem cells carrying a long QT syndrome type 2 mutation. Eur Heart J 32: 952–962 [DOI] [PMC free article] [PubMed] [Google Scholar]
  42. McKernan KJ, Peckham HE, Costa GL, McLaughlin SF, Fu Y, Tsung EF, Clouser CR, Duncan C, Ichikawa JK, Lee CC, Zhang Z, Ranade SS, Dimalanta ET, Hyland FC, Sokolsky TD, Zhang L, Sheridan A, Fu H, Hendrickson CL, Li B et al. (2009) Sequence and structural variation in a human genome uncovered by short-read, massively parallel ligation sequencing using two-base encoding. Genome Res 19: 1527–1541 [DOI] [PMC free article] [PubMed] [Google Scholar]
  43. Montgomery SB, Sammeth M, Gutierrez-Arcelus M, Lach RP, Ingle C, Nisbett J, Guigo R, Dermitzakis ET (2010) Transcriptome genetics using second generation sequencing in a Caucasian population. Nature 464: 773–777 [DOI] [PMC free article] [PubMed] [Google Scholar]
  44. Moretti A, Bellin M, Welling A, Jung CB, Lam JT, Bott-Flugel L, Dorn T, Goedel A, Hohnke C, Hofmann F, Seyfarth M, Sinnecker D, Schomig A, Laugwitz KL (2010) Patient-specific induced pluripotent stem-cell models for long-QT syndrome. N Engl J Med 363: 1397–1409 [DOI] [PubMed] [Google Scholar]
  45. Morita H, Wu J, Zipes DP (2008) The QT syndromes: long and short. Lancet 372: 750–763 [DOI] [PubMed] [Google Scholar]
  46. Moss AJ, Shimizu W, Wilde AA, Towbin JA, Zareba W, Robinson JL, Qi M, Vincent GM, Ackerman MJ, Kaufman ES, Hofman N, Seth R, Kamakura S, Miyamoto Y, Goldenberg I, Andrews ML, McNitt S (2007) Clinical aspects of type-1 long-QT syndrome by location, coding type, and biophysical function of mutations involving the KCNQ1 gene. Circulation 115: 2481–2489 [DOI] [PMC free article] [PubMed] [Google Scholar]
  47. 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 (2002) 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 105: 794–799 [DOI] [PubMed] [Google Scholar]
  48. Muller FJ, Schuldt BM, Williams R, Mason D, Altun G, Papapetrou EP, Danner S, Goldmann JE, Herbst A, Schmidt NO, Aldenhoff JB, Laurent LC, Loring JF (2011) A bioinformatic assay for pluripotency in human cells. Nat Methods 8: 315–317 [DOI] [PMC free article] [PubMed] [Google Scholar]
  49. Mummery C, Ward-van Oostwaard D, Doevendans P, Spijker R, van den Brink S, Hassink R, van der Heyden M, Opthof T, Pera M, de la Riviere AB, Passier R, Tertoolen L (2003) Differentiation of human embryonic stem cells to cardiomyocytes: role of coculture with visceral endoderm-like cells. Circulation 107: 2733–2740 [DOI] [PubMed] [Google Scholar]
  50. Mummery CL, Zhang J, Ng ES, Elliott DA, Elefanty AG, Kamp TJ (2012) Differentiation of human embryonic stem cells and induced pluripotent stem cells to cardiomyocytes: a methods overview. Circ Res 111: 344–358 [DOI] [PMC free article] [PubMed] [Google Scholar]
  51. Narsinh KH, Plews J, Wu JC (2011) Comparison of human induced pluripotent and embryonic stem cells: fraternal or identical twins? Mol Ther 19: 635–638 [DOI] [PMC free article] [PubMed] [Google Scholar]
  52. Nazor KL, Altun G, Lynch C, Tran H, Harness JV, Slavin I, Garitaonandia I, Muller FJ, Wang YC, Boscolo FS, Fakunle E, Dumevska B, Lee S, Park HS, Olee T, D’Lima DD, Semechkin R, Parast MM, Galat V, Laslett AL et al. (2012) Recurrent variations in DNA methylation in human pluripotent stem cells and their differentiated derivatives. Cell Stem Cell 10: 620–634 [DOI] [PMC free article] [PubMed] [Google Scholar]
  53. Panopoulos AD, Ruiz S, Izpisua Belmonte JC (2011) iPSCs: induced back to controversy. Cell Stem Cell 8: 347–348 [DOI] [PubMed] [Google Scholar]
  54. Pillozzi S, Brizzi MF, Balzi M, Crociani O, Cherubini A, Guasti L, Bartolozzi B, Becchetti A, Wanke E, Bernabei PA, Olivotto M, Pegoraro L, Arcangeli A (2002) HERG potassium channels are constitutively expressed in primary human acute myeloid leukemias and regulate cell proliferation of normal and leukemic hemopoietic progenitors. Leukemia 16: 1791–1798 [DOI] [PubMed] [Google Scholar]
  55. Pini G, Scusa MF, Congiu L, Benincasa A, Morescalchi P, Bottiglioni I, Di Marco P, Borelli P, Bonuccelli U, Della-Chiesa A, Prina-Mello A, Tropea D (2012) IGF1 as a potential treatment for Rett syndrome: safety assessment in six Rett patients. Autism Res Treat 2012: 679801. [DOI] [PMC free article] [PubMed] [Google Scholar]
  56. Reinhardt P, Schmid B, Burbulla LF, Schondorf DC, Wagner L, Glatza M, Hoing S, Hargus G, Heck SA, Dhingra A, Wu G, Muller S, Brockmann K, Kluba T, Maisel M, Kruger R, Berg D, Tsytsyura Y, Thiel CS, Psathaki OE et al. (2013) Genetic correction of a LRRK2 mutation in human iPSCs links Parkinsonian neurodegeneration to ERK-dependent changes in gene expression. Cell Stem Cell 12: 354–367 [DOI] [PubMed] [Google Scholar]
  57. Rosati B, Marchetti P, Crociani O, Lecchi M, Lupi R, Arcangeli A, Olivotto M, Wanke E (2000) Glucose- and arginine-induced insulin secretion by human pancreatic beta-cells: the role of HERG K(+) channels in firing and release. FASEB J 14: 2601–2610 [DOI] [PubMed] [Google Scholar]
  58. Rudy Y (2008) Molecular basis of cardiac action potential repolarization. Ann NY Acad Sci 1123: 113–118 [DOI] [PubMed] [Google Scholar]
  59. Sinnecker D, Goedel A, Laugwitz KL, Moretti A (2013) Induced pluripotent stem cell-derived cardiomyocytes: a versatile tool for arrhythmia research. Circ Res 112: 961–968 [DOI] [PubMed] [Google Scholar]
  60. Smith JL, McBride CM, Nataraj PS, Bartos DC, January CT, Delisle BP (2011) Trafficking-deficient hERG K(+) channels linked to long QT syndrome are regulated by a microtubule-dependent quality control compartment in the ER. Am J Physiol Cell Physiol 301: C75–C85 [DOI] [PMC free article] [PubMed] [Google Scholar]
  61. Soldner F, Laganiere J, Cheng AW, Hockemeyer D, Gao Q, Alagappan R, Khurana V, Golbe LI, Myers RH, Lindquist S, Zhang L, Guschin D, Fong LK, Vu BJ, Meng X, Urnov FD, Rebar EJ, Gregory PD, Zhang HS, Jaenisch R (2011) Generation of isogenic pluripotent stem cells differing exclusively at two early onset Parkinson point mutations. Cell 146: 318–331 [DOI] [PMC free article] [PubMed] [Google Scholar]
  62. Szuhai K, Tanke HJ (2006) COBRA: combined binary ratio labeling of nucleic-acid probes for multi-color fluorescence in situ hybridization karyotyping. Nat Protoc 1: 264–275 [DOI] [PubMed] [Google Scholar]
  63. Takahashi K, Okita K, Nakagawa M, Yamanaka S (2007a) Induction of pluripotent stem cells from fibroblast cultures. Nat Protoc 2: 3081–3089 [DOI] [PubMed] [Google Scholar]
  64. Takahashi K, Tanabe K, Ohnuki M, Narita M, Ichisaka T, Tomoda K, Yamanaka S (2007b) Induction of pluripotent stem cells from adult human fibroblasts by defined factors. Cell 131: 861–872 [DOI] [PubMed] [Google Scholar]
  65. Tan HL, Bardai A, Shimizu W, Moss AJ, Schulze-Bahr E, Noda T, Wilde AA (2006) Genotype-specific onset of arrhythmias in congenital long-QT syndrome: possible therapy implications. Circulation 114: 2096–2103 [DOI] [PubMed] [Google Scholar]
  66. Tester DJ, Will ML, Haglund CM, Ackerman MJ (2005) Compendium of cardiac channel mutations in 541 consecutive unrelated patients referred for long QT syndrome genetic testing. Heart Rhythm 2: 507–517 [DOI] [PubMed] [Google Scholar]
  67. Vandenberg JI, Perry MD, Perrin MJ, Mann SA, Ke Y, Hill AP (2012) hERG K(+) channels: structure, function, and clinical significance. Physiol Rev 92: 1393–1478 [DOI] [PubMed] [Google Scholar]
  68. Wang Y, Huang X, Zhou J, Yang X, Li D, Mao H, Sun HH, Liu N, Lian J (2012) Trafficking-deficient G572R-hERG and E637K-hERG activate stress and clearance pathways in endoplasmic reticulum. PLoS One 7: e29885. [DOI] [PMC free article] [PubMed] [Google Scholar]
  69. Yazawa M, Hsueh B, Jia X, Pasca AM, Bernstein JA, Hallmayer J, Dolmetsch RE (2011) Using induced pluripotent stem cells to investigate cardiac phenotypes in Timothy syndrome. Nature 471: 230–234 [DOI] [PMC free article] [PubMed] [Google Scholar]
  70. Yu J, Vodyanik MA, Smuga-Otto K, Antosiewicz-Bourget J, Frane JL, Tian S, Nie J, Jonsdottir GA, Ruotti V, Stewart R, Slukvin II, Thomson JA (2007) Induced pluripotent stem cell lines derived from human somatic cells. Science 318: 1917–1920 [DOI] [PubMed] [Google Scholar]
  71. Zou J, Maeder ML, Mali P, Pruett-Miller SM, Thibodeau-Beganny S, Chou BK, Chen G, Ye Z, Park IH, Daley GQ, Porteus MH, Joung JK, Cheng L (2009) Gene targeting of a disease-related gene in human induced pluripotent stem and embryonic stem cells. Cell Stem Cell 5: 97–110 [DOI] [PMC free article] [PubMed] [Google Scholar]

Associated Data

This section collects any data citations, data availability statements, or supplementary materials included in this article.

Supplementary Materials

Supplementary Information

Review Process File