Global transcriptional profiles of beating clusters derived from human induced pluripotent stem cells and embryonic stem cells are highly similar

Background Functional and molecular integrity of cardiomyocytes (CMs) derived from induced pluripotent stem (iPS) cells is essential for their use in tissue repair, disease modelling and drug screening. In this study we compared global transcriptomes of beating clusters (BCs) microdissected from differentiating human iPS cells and embryonic stem (ES) cells. Results Hierarchical clustering and principal component analysis revealed that iPS-BCs and ES-BCs cluster together, are similarly enriched for cardiospecific genes and differ in expression of only 1.9% of present transcripts. Similarly, sarcomeric organization, electrophysiological properties and calcium handling of iPS-CMs were indistinguishable from those of ES-CMs. Gene ontology analysis revealed that among 204 genes that were upregulated in iPS-BCs vs ES-BCs the processes related to extracellular matrix, cell adhesion and tissue development were overrepresented. Interestingly, 47 of 106 genes that were upregulated in undifferentiated iPS vs ES cells remained enriched in iPS-BCs vs ES-BCs. Most of these genes were found to be highly expressed in fibroblasts used for reprogramming and 34% overlapped with the recently reported iPS cell-enriched genes. Conclusions These data suggest that iPS-BCs are transcriptionally highly similar to ES-BCs. However, iPS-BCs appear to share some somatic cell signature with undifferentiated iPS cells. Thus, iPS-BCs may not be perfectly identical to ES-BCs. These minor differences in the expression profiles may occur due to differential cellular composition of iPS-BCs and ES-BCs, due to retention of some genetic profile of somatic cells in differentiated iPS cell-derivatives, or both.

Although comparisons between iPS cells and conventional ES cells revealed that they have highly similar growth characteristics, gene expression profiles, epigenetic status and developmental potential [3,5,6,27,28], recent comprehensive analyses of various undifferentiated ES and iPS cell lines showed that iPS cells may not be perfectly identical to conventional ES cells at the molecular level [29][30][31]. These studies demonstrated that iPS cells are characterized by a unique gene and miRNA expression signature as well as a CpG methylation pattern, which distinguish them from ES cells. However, the comparison of global transcriptomes of mature cells differentiated from ES and iPS cells has not yet been performed, and it is not clear whether molecular differences between iPS and ES cells are retained upon their differentiation into mature cells.
Full transcriptional profiles of human ES-CMs have been reported by several groups. Synnergren and coworkers differentiated the human ES cell line SA002 to cardiac lineage in an embryoid body (EB) system and CMs were enriched by mechanical dissection of spontaneously beating clusters (BCs) [32]. The results of this study indicated that human ES-BCs, despite being composed of different cell types, are highly enriched for CM-specific transcripts and display high similarities to human fetal heart tissue. Cao and coworkers reported the transcriptional profile of Percoll density gradientenriched human H9 ES-CMs that were isolated by centrifugation to a purity of about 45% [33], and most recently, microarray analyses were carried out with highly purified CMs (> 99%) that were generated by drug selection from the transgenic human ES cell line HES3 [34] and H9 [35]. These studies provide the first comprehensive characterization of the molecular signature of purified human ES-CMs revealing high similarities between expression profiles of ES-CMs and human native CMs.
The purpose of this study was to determine the degree of molecular similarity between CMs and other differentiated cell types present in microdissected BCs derived from human iPS cells and conventional ES cells. Our data indicates that iPS-BCs are transcriptionally highly similar to ES-BCs. However, iPS-BCs appear to share some somatic cell signature with undifferentiated iPS cells and express higher levels of transcripts encoding for some extracellular matrix components and cell adhesion molecules than ES-BCs. This expression profile may reflect either the differential cellular composition or activity in iPS-and ES-BCs, partial retention of genetic signature of somatic cells in differentiated iPS cell-derivatives, or both. These findings emphasize the necessity for a careful assessment of the differentiation capacity of iPS cells of various somatic cell origins and molecular profiles of specific mature iPS and ES cell-derivatives, so as to determine whether any alterations may exist that could affect their use for regenerative medicine and research.

Cardiac differentiation of human iPS and ES cells
The human iPS cell line derived from foreskin fibroblasts, clone 1 (C1), [6] and the human ES cell line HES-2 were differentiated to CMs using the END2 coculture system [15] (Additional file 1, Figure S1). Upon inducing differentiation, spontaneously contracting clusters were first observed at day 11 in cultures of both cell lines. At day 15 of differentiation, the fraction of clusters that exhibited spontaneous contractions in iPS and ES cell cultures was 12.8 ± 3.5% and 16.0 ± 4.8%, respectively (n = 9, p = 0.128), (Additional file 1, Figure  S2A). With further cultivation no additional beating areas appeared. In BCs derived from iPS and ES cells the CMs represented, respectively, 16.7 ± 5.1% and 15.0 ± 2.6% of all cells (Additional file 1, Figure S2B).
The contraction rates of these BCs were similar in iPS and ES cell cultures (46 ± 6.8 beats/min for iPS cells and 43 ± 4.5 beats/min for ES cells), (Additional file 1, Figure S2C). These data indicate that the iPS cell line Foreskin C1 differentiates to CMs at rates comparable to the ES cell line HES-2 and possesses similar beating characteristics.

Expression of cardiospecific genes and sarcomeric proteins in BCs
To confirm the presence of CMs in BCs isolated from differentiating iPS and ES cells we examined the expression of cardiospecific markers using qRT-PCR (primer sequences are given in Additional file 1, Table S1) and immunocytochemistry. The data on expression of the endogenous and viral pluripotency marker OCT4 in undifferentiated and differentiated iPS and ES cells is discussed in Additional file 1 ( Figure S3 and Supplementary results). The expression of mRNAs for transcription factors NKX2.5 and GATA4, and the cardiac structural proteins myosin light chain 2 ventricular isoform (MLC2v), α-actinin (ACTN2) and α-myosin heavy chain (MYH6) was undetectable or very low in undifferentiated cells ( Figure 1A). However, their expression was strongly upregulated in BCs at day 18 of differentiation ( Figure 1A). The specific marker of skeletal muscle, myosin heavy chain 2 (MYH2), was not detectable in these aggregates ( Figure 1A). Similar levels of the cardiac genes in microdissected BCs of both cell types indicate that they contain comparable amounts of CMs, which is in agreement with the CM-counting data in the Additional file 1, Figure S2B. Immunocytochemical analyses revealed that CMs in day 18 iPS-and ES-BCs stain positively for cardiac proteins α-actinin and troponin T (cTnT) and display typical pattern of cross-striations indicative of sarcomeric organization in CMs ( Figure 1B). The highly ordered striated pattern is highlighted in high resolution images and clearly recapitulates the normal architecture of the contractile apparatus in functional CMs.

Intact adrenergic and muscarinergic signaling in iPS-CMs
MEA measurements were performed to verify the functional integrity of CMs in microdissected BCs. In the Figure 1 Expression of cardiac transcripts and proteins in cardiomyocytes derived from human iPS and ES cells. A. qRT-PCR assessment of expression of cardiac specific genes NKX2.5, MLC2v, ACTN2, MYH6, GATA4 and a skeletal muscle gene MYH2 in undifferentiated iPS and ES cells and their corresponding beating clusters (BC) microdissected at day 18 of differentiation. The data are means of triplicate analyses ± SEM. B. Anti-sarcomeric cardiac actinin and cardiac troponin T stainings in HES2-(upper panels) and iPS(Forskin) C1-derived (lower panels) cardiomyocytes. Spontaneously beating outgrowths were microdissected at day 15 of differentiation and replated on fibronectin-coated plates. Adherent beating clusters were used for staining 3-5 days later. Scale bar: 20 μm. absence of pharmacological stimulation, iPS-CMs and ES-CMs exhibited similar spontaneous electrical activity with the field potential (FP) frequencies of 1.04 ± 0.39 Hz (n = 11) and 1.15 ± 0.66 Hz (n = 9), respectively ( Figure 2 and Additional file 1, Table S2). Application of β-adrenergic receptor agonist isoproterenol (ISO) significantly increased the basal FP frequency in both cell populations to 1.75 ± 0.44 Hz (iPS cells, n = 5) and 1  its reduction in iPS-BCs and ES-BCs to 0.23 ± 0.21 Hz (n = 6) and 0.35 ± 0.42 Hz (n = 5), respectively ( Figure  2, and Additional file 1, Table S2). These results suggest that iPS-CMs express functional β-adrenergic and muscarinic receptors and their associated signaling pathways in a manner indistinguishable from ES-CMs.

Spontaneous Ca 2+ transients in iPS-CMs
To further characterize the functional properties of iPS-CMs, we have compared spontaneous Ca 2+ release in BCs of human ES and iPS origin using fluorometric Ca 2+ imaging with Fura-2 AM ( Figure 3A). Both cell types showed rhythmic Ca 2+ transients with identical baseline and amplitude, maximum upstroke-and maximum decay velocity, time to peak, total transient time and time to half peak relaxation ( Figure 3B, 3C and Additional file 1,    Table S6). Principal component analysis (PCA) with this data set showed that samples cluster in four groups corresponding to a) undifferentiated ES and iPS cells, b) microdissected ESand iPS-BCs, c) fetal heart and d) adult heart ( Figure  4A). Hierarchical clustering also showed that undifferentiated iPS cells cluster together with undifferentiated ES cells and iPS-BCs together with ES-BCs ( Figure  4B). The transcriptional profiles of FH and AH were more related to those of iPS-BCs and ES-BCs than to those of their undifferentiated counterparts. Furthermore, scatter-plot analysis of every probe set on the array further emphasized a close correlation of gene expression levels between iPS and ES cells (r 2 = 0.96 ± 0.01, n = 9) as well as between iPS-BCs and ES-BCs (r 2 = 0.93 ± 0.04, n = 9) ( Figure 4C). In contrast, correlation between iPS cells and iPS-BCs as well as between undifferentiated ES cells and ES-BCs was significantly lower (r 2 = 0.87 ± 0.02, n = 9). A closer look at differentially expressed genes (p-value < 0.05; > 2fold change; difference between mean intensity signals > 100) revealed that only 198 probe sets (1.19% of present transcripts or 4.65% of variable genes) significantly differ in their expression levels between ES and iPS cells and only 330 probe sets (1.91% of present transcripts and 7.76% of variable genes) were differentially expressed between ES-BCs and iPS-BCs ( Figure 4D). Complete lists of up-and down-regulated genes in these comparisons are provided in Additional file 3 (Tables S7 and S8) and Additional file 4 (Tables S9  and S10).

Genes commonly upregulated in iPS cells and iPS-BCs
In order to determine whether any of the 106 unique transcripts found to be overexpressed in iPS cells in iPS-versus ES-cell comparison also remain overexpressed in iPS-BCs compared to ES-BCs we searched for overlapping genes in these two lists of overexpressed genes and found that among 106 unique iPS cellenriched transcripts 47 (44.3%) were also present on the list of 204 iPS-BC-enriched transcripts ( Figure 4E and Table 1). Using Fisher's exact test we have determined that this overlap was highly statistically significant (p = 6,7e-28). Moreover, 38 (80.9%) of these overlapping transcripts remained expressed at similar levels in undifferentiated iPS cells and differentiated iPS-BCs ( Table  1). The gene ontology (GO) analysis of iPS cell-upregulated transcripts did not reveal significant enrichment for any GO-term in this group of genes (data not shown). Verification of these microarray data with quantitative PCR revealed that out of 11 genes selected from this list 10 genes (COMT, DYNLT3, NME4, OXCT1, MGMT, PTGR1, MGC3207, CKLF, ZNF167, ZNF626, and RPL39L) were expressed at very low levels or not at all in ES cells and ES-BCs, but were highly upregulated in iPS cells and iPS-BCs ( Figure 5, for comparison see microarray expression data in Table 1). Interestingly, all of these 10 iPS/iPS-BC-enriched genes were also strongly expressed in fibroblasts used for reprogramming, their expression levels in fibroblasts being similar to those in iPS cells and iPS-BCs ( Figure 5). These data suggest that this group of genes may belong to fibroblast specific-genes that could not be repressed in the process of reprogramming and still remain upregulated in iPS cells. Indeed, two recent studies demonstrated that iPS cells are characterized by a unique gene expression signature, which distinguishes them from ES cells [29,31]. By comparing these data with our list of 47 transcripts upregulated in iPS cells and iPS-BCs we found that 16 of these transcripts (34%) were also detected in these previous studies to be specifically expressed in multiple human iPS cell lines in comparison with various ES cell lines (Table 1). Thus, iPS cell-signature genes appear to be specifically expressed in iPS-BCs but not ES-BCs.

Gene ontology analysis of genes upregulated in iPS-BCs vs ES-BCs
In order to determine which biological terms are enriched in iPS-BCs compared to ES-BCs, we have performed a functional annotation clustering analysis of 204 unique transcripts that were found to be upregulated in iPS-BCs (see Additional file 4, Table S7). Considering only the significant GO-terms with the p-value < 0.001 this analysis revealed the overrepresentation in  iPS-BCs of functional categories related to collagen, extracellular matrix, cell adhesion and tissue development and morphogenesis (Table 2). Genes in these categories included various types of collagens (COL1A1, COL3A1, COL4A1, COL5A1, COL5A2, COL6A3) and other extracellular matrix components (LUM, LAMB1, CTGF), cell surface proteins involved in cell-cell interactions (CD36, CD44, CD47) and a variety of transcription factors involved in tissue development and morphogenesis (HOXA5, HOXB5, SNAI2, NR2F2). These transcripts were increased in iPS-BCs vs ES-BCs in average by 3,3 ± 2,1 -fold (Additional file 5, Table S11). Since these transcripts are expressed by specific cell types (e.g. collagens by fibroblasts) or involved in specific biological processes (e.g. HOXB5 in lung and gut development), these findings indicate that transcriptional differences between iPS-BCs and ES-BCs may partially occur due to higher abundance (or higher transcriptional activity) of specific non-cardiac cell types in these clusters.

Microarray analysis of human ES cell-and CM-specific and enriched genes
We have further compared the expression levels of selected transcripts known to be specifically expressed or enriched in undifferentiated ES cells or CMs. 93.4% of human ES cell-specific or -enriched genes were expressed at similar levels in iPS and ES cells ( Figure 6 and Additional file 6, Tables S12 and S13). Furthermore, 89.2% of 83 cardiac specific or enriched genes were expressed at similar levels in iPS-BCs and ES-BCs and out of 34 BC-or CM-associated genes identified in previous studies [32][33][34][35][36] 28 (82.4%) were enriched to a similar extent in both iPS-BCs and ES-BCs in our study ( Figure 6 and Additional file 6, Table S14). Further elaboration of these microarray data is given in the Additional file 1, Supplementary results. Microarray data were validated by comparative qPCR analysis of undifferentiated iPS and ES cells with Human embryonic stem cell 96 StellARray™ qPCR Array (Lonza). This analysis revealed that 82 genes (91.1%) were detected with both methods as being expressed at undistinguishable levels in ES and iPS cells (p > 0.05 and fold change < 2) (Additional file 7, Table S15).

Comparison of transcriptional profiles between undifferentiated and differentiated iPS and ES cells
While transcriptomes of BCs derived from iPS and ES cells were highly similar, as many as 2394 (13.73% of present transcripts) and 2106 probe sets (12.75% of present transcripts) were differentially expressed between undifferentiated iPS cells and iPS-BCs and between undifferentiated ES cells and ES-BCs, respectively (Figure 4D). Among the most upregulated genes in iPS-BCs and ES-BCs were known cardiac genes and the most downregulated genes in both types of BCs belonged to pluripotency markers. More extensive discussion on this subject is provided in the Supplementary results (Additional file 1) and the lists of differentially expressed genes can be found in Additional file 8, Tables S16-S19. Functional annotation analysis was performed in DAVID to identify the most important biological processes that are significantly enriched in iPS-and ES-BCs compared to their undifferentiated counterparts. Taking into account only the GO-terms that had more than 3 genes and p-value < 0,001 this analysis revealed that in categories biological process, molecular function, cellular component and KEGG pathway the overlap between terms enriched in iPS-BCs and ES-BCs was highly significant as determined by Fisher's exact test (Figure 7 and Additional file 9, Table S20). Among the most enriched GO-terms related to cardiomyocytes were heart development, blood vessel and muscle tissue development, myofibril assembly, muscle contraction, ECM structural constituent, calcium ion binding, actin binding, contractile fiber, myofibril, sarcomere, actin cytoskeleton, collagen, Z disc, I band and A band, focal   adhesion, dilated, hyperthrophic and arrhythmogenic cardiomyopathy, and ALK pathway in cardiac myocytes ( Figure 7 and Additional file 9, Table S20). However, a number of GO-terms were differentially enriched between iPS-BCs and ES-BCs. For example, in the category biological process the iPS-BCs but not ES-BCs were enriched for GO-terms such as kidney development, mesenchyme development, angiogenesis, neural crest cell development, central nervous system development and others, while ES-BCs were enriched, among others, for GO-terms sarcomere organization, muscle thin filament assembly, heart septum morphogenesis and skeletal myofibril assembly (Additional file 9, Table  S20). Although these differences may reflect a certain degree of heterogeneity in cellular composition of iPSand ES-BCs, most categories were enriched in both iPS-BCs and ES-BCs to the same extent, suggesting that they are very similar.

Discussion
In this study we show that global gene expression profiles of functionally intact CMs and other differentiated cell types enriched in microdissected iPS-BCs and ES-BCs are highly similar. Cardiac specific genes identified in previous reports using enriched [32,33,36] or highly purified human ES-CM-preparations [34,35] were similarly expressed in both iPS-BCs and ES-BCs and only a few percent of transcripts classified as present among all samples were differentially expressed between these two cell populations. The small difference in transcriptomes between iPS-BCs and ES-BCs was not significantly greater than the difference between undifferentiated iPS and ES cells or between different ES cell or different iPS cell lines as reported by others [37][38][39][40].
Of interest is the observation that 44.3% of transcripts (47 out of 106) that were upregulated in iPS cells compared to ES cells also remained upregulated in iPS-BCs compared to ES-BCs. Statistical analysis revealed that this overlap is unlikely to occur by chance. Many of these genes were expressed at very high levels in iPS cells and iPS-BCs but were absent from both undifferentiated ES cells and ES-BCs. We could not identify these genes on the list of 530 genes reported to be upregulated in microdissected BCs derived from human ES cell line SA002 by Synnergren and coworkers [32]. However, the significant fraction of these genes (16 out of 47) has been found by other groups to be specifically expressed in various human iPS cell lines compared to ES cell lines [29,31]. In one another comparison, one third of our 47 genes was also found to be overexpressed in at least five different undifferentiated iPS cell lines when compared to several independently derived human ES cell lines (Natalia Polouliakh, Sony Computer Science Laboratories, Tokyo, Japan, personal communication). Thus, these genes may belong to a unique gene expression signature of iPS cells [29,31]. Interestingly, many of these 47 genes have been found to be expressed in foreskin fibroblasts [6] and this was also confirmed by qRT-PCR in our study for 10 randomly selected genes. Thus, this subset of genes may represent fibroblast-associated genes that are expressed in iPS cells because they are refractory to epigenetic silencing during reprogramming, as suggested previously [29,31]. Recent meta-analysis of multiple transcriptional data sets derived from human fibroblasts, ES and iPS cells also revealed that iPS cells share four-times more genes in common with fibroblasts from which they were derived than ES cells and fibroblasts [40].
While the retention of epigenetic memory of somatic cells in undifferentiated iPS cells may explain a subset of differentially expressed genes in iPS vs ES cell comparison, it is more difficult to discern the mechanism behind the highly significant overlap between genes upregulated in iPS vs ES cells and those upregulated in   iPS-BC vs ES-BCs. It is unlikely that differential expression of these genes resulted from differential amounts of contaminating murine feeder cells in BCs derived from two types of stem cells, because iPS and ES cells were maintained and differentiated under the same conditions for all three respective biological replicates. Differential contamination of BCs with pluripotent stem cells is also unlikely, because qRT-PCR analyses demonstrated that the expression of endogenous OCT4 decreased to a similar extent upon differentiation of both ES and iPS cells, and SOX2 could not even be detected in both types of BCs. We have also excluded the possibility that some of the common iPS cell-and iPS-BC-enriched genes are upregulated as a consequence of ectopic reactivation of lentiviral OCT4 gene in iPS-BCs, because these genes had not been reported to be the Oct4 targets in previous studies [41][42][43]. Instead, our analyses indicate that transcriptional differences between iPS-BCs and ES-BCs may, at least partially, be caused by differential cellular composition in these clusters. As revealed by gene ontology analysis, biological terms related to collagen, extracellular matrix, cell adhesion and tissue development were enriched in iPS-BCs compared to ES-BCs by about 3-fold. It remains to be determined whether higher content of fibroblasts in iPS-BCs (or other cell types expressing the overrepresented transcripts) is due to specific properties of the iPS cell line used in this study or due to inherently better capability of different iPS cell lines to differentiate towards cells of their somatic origin. Just recently, two studies using murine iPS cell lines derived from different somatic cell types demonstrated that iPS cells retain some DNA methylation signature characteristic of their somatic tissue of origin and that their cellular origin influences the in vitro differentiation propensity of iPS cells [44,45]. The disadvantage of using microdissected beating clusters in our study, which contain CMs and multiple other cell types, is that we are not able to provide a definitive answer to the question whether epigenetic memory of somatic cells may still be retained in some mature iPS cell derivatives. To answer this question, transcriptional profiles of undifferentiated iPS cells and highly purified differentiated cells derived from several iPS cell lines should be compared. We have compared transcriptional profiles of highly purified cardiomyocytes generated by drug selection from a transgenic murine iPS and ES cell line, but could not find any evidence for retention of somatic cell signature in murine iPS-CMs (Fatima A., manuscript in preparation). The same observation has been also recently reported for murine cardiac Nkx2.5-positive progenitor cells isolated by drug selection from several murine iPS and ES cell lines [46]. In this study, the variability in gene expression profiles that was seen between undifferentiated iPS and ES cell lines was even reduced in highly purified lineage-specific cardiac progenitors. Thus, it is unlikely that CMs present in iPS-BCs express any of upregulated genes listed in our Table 1. However, we can not exclude the possibility that other cell types in iPS-BCs still ectopically express some of the genes in this list as a consequence of retention of epigenetic memory of somatic cells. Clearly, additional studies with different cell lines and purified differentiated cell populations will be required to discern between these various possibilities.
A significant progress made in the generation of genetically intact iPS cells [47,48] renders the investigation of consequences of transgene reactivation on the physiology of differentiated cells less relevant. The iPS cell line used in our present report contain sequences of reprogramming transgenes stably integrated into their genome [6] and this may have caused some of the differences in gene expression profiles detected between ES and iPS cells. This notion is supported by observation that iPS cells generated by non-integrating reprogramming methods are transcriptionally more similar to ES cells than those generated with stably integrating viral vectors [40]. Therefore, genetically intact iPS cell lines should be utilized in future studies. However, it is worth mentioning that the genetic integration of viral sequences and incomplete silencing of the lentiviral OCT4 transgene in iPS-BCs as observed in our study, and corroborated by Zhang and coworkers [22], apparently did not adversely affect the functional as well as molecular properties of iPS-CMs. This is in agreement with the observation that ectopic Oct4 expression in the heart of adult mice did not lead to any detectable alterations of cell phenotype [49]. Similarly, 303 transcripts differentially expressed between iPS-BC and ES-BCs did not appear to exert a measurable affect on functional properties of CMs contained in iPS-BCs in vitro.

Conclusions
Transcriptional profiles of undifferentiated iPS and ES cells have been extensively analyzed in the past. In this study we provide molecular evidence that not only undifferentiated iPS and ES cells but also their differentiated derivatives contained in spontaneously beating areas have highly similar global transcriptional profiles. Minor differences in gene expression that were identified between iPS and ES cells in undifferentiated and differentiated state did not appear to significantly affect the structural or functional properties of iPS cell-derived cardiomyocytes. However, iPS-BCs appear to retain some of the iPS cell genetic signature partially due to enrichment of specific cell types (such as those expressing collagen and ECM-components) and thus may not be perfectly identical to ES-BCs. These findings emphasize the necessity for detailed functional and molecular analyses of mature iPS cell-derivatives so as to determine the consequences of their expression for the physiology and safety of iPS cell-derivatives.

Culture of undifferentiated human ES and iPS cells
The human iPS cell line derived from foreskin fibroblasts, clone 1 (C1), was kindly provided by James Thomson (University of Wisconsin, Madison, WI, USA) [6]. The informed consent of the patient donating the tissue probe for iPS cell generation was obtained by the original group (see the reference 6). For comparison we used the human ES cell line HES-2 that was generated by the ES Cell International (Singapore, http://www. escellinternational.com/about_esi/index.html) and obtained from the repository at the WiCell Research Institute (Madison, WI, USA, http://www.wicell.org/). The iPS and ES cells were maintained on irradiated murine embryonic fibroblasts (MEFs) in DMEM/F12 medium supplemented with Glutamax, 20% knockout serum replacer, 1% nonessential amino acids (NAA), 0.1 mmol/L β-mercaptoethanol (βME) (Invitrogen, Carlsbad, CA, http://www.invitrogen.com) and with 4 or 100 ng/ml basic fibroblast growth factor (Peprotech, Rocky Hill, NJ, USA, http://www.peprotech.com) for ES cells or iPS cells, respectively. Culture media were changed daily and the cells were passaged by manual dissection of cell clusters every 5-6 days. Work with human ES cells has been approved by the regulatory authorities at the Robert Koch Institute, Berlin, Germany (permission number 1710-79-1-4-2-A10).

Cardiac differentiation
Cardiac differentiation of human iPS and ES cells was carried out on the murine visceral endoderm-like cell line END2 [15,50]. END2 cell cultures were mitotically inactivated for 3 hours with 10 μg/ml mitomycin C (Sigma-Aldrich) and 1.2 × 10 6 cells were plated on 60 mm tissue culture dishes one day before starting the coculture. To initiate co-cultures, iPS and ES cell colonies were dissociated into clumps by either using collagenase IV (1 mg/ml in DMEM/F-12 at 37°C for 5-10 minutes) or by manual cutting. The differentiation was carried out in Knockout-DMEM having 1 mM Lglutamine, 1% NAA, 0.1 mmol/L (βME) and Penicillin/Streptomycin (100 U/ml and 100 μg/ml, respectively) devoid of serum and serum replacement (all reagents were from Invitrogen). Medium change was performed at day 5, 9, 12 and 15 after initiating the co-cultures.

RT-PCR and quantitative RT-PCR
Total RNA was isolated using TRIzol Reagent (Invitrogen) from iPS or ES cells and from 40-60 iPS cellderived BCs (iPS-BCs) or ES cell-derived BCs (ES-BCs) microdissected at day 18 of differentiation. DNase Itreated total RNA (500 ng) was reverse-transcribed using Superscript II RTase (Invitrogen) and random hexamers. cDNA was diluted 1:4 with sterile tri-destilled water and 5 μl were amplified using JumpStart™RedTaq ReadyMix™PCR Reaction Mix (Sigma). Negative controls were generated in RT reactions in which all reaction components were included except RTase. Reactions were terminated at the exponential phase of amplification and products were analyzed by agarose gel electrophoresis. For quantitative RT-PCR the cDNA probes were diluted 1:40 and 2 μl was amplified using SYBR Green PCR Master Mix (Qiagen, Hilden, Germany, http://www.qiagen.com) in triplicate for each sample and each gene. Real-time PCRs were performed in a 7500 Fast System Real Time Cycler (Applied Biosystems, Foster City, CA, USA, http://www.appliedbiosystems. com) and analyzed with SDSShell 1.4 software (Applied Biosystems). GAPDH was used for normalization of expression levels of individual genes. Primers used are listed in Additional file 1, Table S1. Microarray results obtained for undifferentiated iPS and ES cells were validated using Human Embryonic Stem Cell 96 StellAR-ray™ qPCR array (Lonza, Cologne, Germany, http:// www.lonza.com) (see Supplementary methods in Additional file 1).
Immuncytochemistry iPS-BCs and ES-BCs were microdissected on day 15 of differentiation and plated on fibronectin-coated (2 μg/ ml) μ-dishes 35 mm, low (Ibidi GmbH, Munich, Germany, http://www.ibidi.de). Three to five days after plating, the cells were fixed with ice-cold methanol at -20°C. After washing with phosphate buffered saline (PBS, pH 7.4), the cells were blocked for one hour with 5% bovine serum albumin (BSA) in PBS.

Measurement of intracellular Ca 2+
Ca 2+ imaging experiments were performed as described previously [53]. Briefly, BCs were isolated on day 18 of differentiation and loaded with 2.5 μM Fura-2 AM (Molecular Probes, Eugene, OR, USA, http://www.molecularprobes.com) and 0.02% Pluronic-127 (Sigma-Aldrich) in modified Tyrode solution containing (in mmol/l): NaCl 140, KCl 5.4, MgCl 2 1, sodium pyruvate 2, CaCl 2 2, HEPES 10, glucose 10 (pH 7.4) for 30 minutes at room temperature [54]. Following the loading procedure the BCs were washed and incubated for 30 minutes at 37°C for deesterification before measurements. After alternating monochromatic excitation with 340 nm and 380 nm, emitted light was detected using a charge-coupled device cooled camera (TILL IMAGO CCD, TILL Photonics, Planegg, Germany, http://www. till-photonics.com) with a sampling rate of 40 Hz. The cell aggregates were perfused with Tyrode solution at 1 ml/min. Two perfusion systems were used: chamber perfusion with approx. 1 cm distance to the clusters and 'puff'-perfusion placed directly in front of the cells. The latter was used for caffeine (10 mmol/L) application in order to estimate the sarcoplasmic reticulum (SR) Ca 2+ load. Ratio F340/F380 data were background subtracted (TILLvisION, TILL Photonics) and transient parameters were analyzed with Chart software v5 (ADInstruments, Castle Hill, Australia, http://www.adinstruments.com). Results are expressed as mean ± SEM. Statistical analysis was performed by Student's t-test or paired t-test (p-value < 0.05 was considered significant).

Microarray procedure and analysis
Total RNA was extracted from undifferentiated ES and iPS cells and their corresponding BCs microdissected at day 18 of differentiation using the TRIzol reagent. Human RNA from fetal (FH) and adult hearts (AH) was purchased from Clontech (Saint-Germain-en-Laye, France, http://www.clontech-europe.com/). The characteristics of these samples are summarized in Additional file 1, Table S5. Biotinlabeled cRNA preparation for the Illumina platform was performed using the Ambion® Illumina RNA amplification kit (Ambion Europe, Huntington, Cambridgeshire, UK). Afterwards, biotin-labeled cRNA (750 ng) was hybridized to Sentrix® whole genome bead chips comprising 48803 probe sets (Human HT-12_V3, Illumina, San Diego, CA, USA, http://www. illumina.com/) and scanned on the Illumina® BeadStation 500×. Across all samples, 16837 ± 1311 probe sets (34.5 ± 2.7%) received present calls as defined by the detection p-value of < 0.05 (n = 3). Raw data extraction of mRNA microarrays was performed with Beadstudio 3.1.1.0 software using the Beadstudio Gene Expression Analysis Module 3.1.8. All further analysis was performed in R (http://www.r-project.org, version 2.8.0) using Bioconductor packages. For further analysis we used quantile normalization implemented in the affy package. Variable genes were defined by a coefficient of variation (SD/mean) between 0.5-10. Determination of present calls was based on the detection p-value assessed by Beadstudio software; a gene was called present if the detection p-value was < 0.05. Otherwise the mRNA transcript was called absent. Differentially expressed genes were selected using a fold change/pvalue filter with the following criteria: only p-values smaller than 0.05 and an expression change higher than 2 fold and a difference between mean intensity signals greater than 100 were considered statistically significant for further analysis. The Benjamini-Hochberg method was used to adjust the raw p-values to control the false discovery rate. The fold-change was calculated by dividing the mean intensity of the genes in one group by that in the other group. If this number was less than one, the negative reciprocal was used. Hierarchical cluster analysis was performed using the hcluster method in R. Before clustering, the data were log2 transformed. Distances of the samples were calculated using Pearson correlation and clusters were formed by taking the average of each cluster. Principal component (PC) analysis was performed using the pcurve package in R. When visualizing PC analysis results, the first 3 principal components (coordinates) were z-transformed (mean = 0, standard deviation = 1) and subsequently plotted in 3 D. The first PC axis accounted for 63.8% of the variance in the data set of variable transcripts, while the second and third PC axes accounted for 20.9% and 0.07% of the variance, respectively All heat maps were visualized using MAYDAY [55]. The enrichment of specific biological processes, molecular functions, cellular components and pathways among differentially expressed genes was analyzed using the Database for Annotation, Visualization and Integrated Discovery (DAVID, http://david.abcc. ncifcrf.gov/) bioinformatics resource [56,57] and functional categories and biological pathways annotated by the Gene ontology tool [58] and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database [59]. Annotations were considered significantly overrepresented when the p-value of the Fisher's exact test as used by DAVID was < 0.001 (enrichment score in functional annotation clustering > 3). Fisher's exact test was performed using an online calculator at http://www. langsrud.com/fisher.htm.

Additional material
Additional file 1: Collection of Figures S1-S3, Tables S1-S5, Supplemental results, methods and references. This file is a PDF file containing Figure S1 (characterization of pluripotency markers in undifferentiated human iPS and ES cells), Figure S2 (assessment of the cardiogenic potential of human iPS and ES cells), Figure S3 (expression of pluripotency genes in beating clusters derived from human iPS and ES cells), Table S1 (primers used for RT-PCR analyses), Additional file 2: List of genes that are differentially expressed across all samples. This file is an Excel file containing Table S6 with the list of 4254 variable genes and their expression levels in all samples. Variable genes were defined by a coefficient of variation (SD/mean) between 0.5-10.
Additional file 3: List of genes differentially expressed between undifferentiated iPS and ES cells. This Excel file contains Tables S7 and S8 with the list of genes that are, respectively, upregulated and downregulated in undifferentiated iPS cells compared to ES cells.
Additional file 4: List of genes differentially expressed between iPS cell-and ES cell-derived beating clusters (BCs). This Excel file contains Tables S9 and S10 with the list of genes that are, respectively, upregulated and downregulated in iPS-BCs compared to ES-BCs.
Additional file 5: Expression levels of selected genes included in GO-terms that were found to be overrepresented in iPS-BCs compared to ES-BCs. This PDF file contains the Table S11 displaying the expression levels of selected fibroblastassociated genes found to be overexpressed in iPS-BCs compared to ES-BCs.
Additional file 6: Expression levels of selected human ES cell-and cardiacspecific or enriched genes across all samples. This is a PDF file in a landscape format containing Table S12 (list of 27 ES cell-specific genes and their expression levels across all samples), Table S13 (list of 49 human ES cell-enriched genes with their expression levels in all samples) and Table S14 (list of 83 cardiac specific or enriched genes and their expression levels across all samples). Additional file 7: qRT-PCR validation of microarray data with Human ES cell StellArray qPCR Array (Lonza). This file is a PDF document containing Table S15 showing the results of qRT-PCR analysis of expression of selected pluripotency and differentiations genes in undifferentiated human ES and iPS cells in comparison with microarray data obtained in this study.
Additional file 8: List of differentially expressed genes in iPS-BC versus iPS cell comparison and ES-BC versus ES cell comparison. This file is an Excel file containing Table S16 (list of 1235 genes upregulated in iPS-BC compared to iPS cells), Table S17 (list of 1014 genes downregulated in iPS-BC compared to iPS cells), Table S18 (list of 875 genes upregulated in ES-BC compared to ES cells), and Table S19 (list of 1103 genes downregulated in ES-BC compared to ES cells).
Additional file 9: Gene ontology analysis of genes upregulated in iPS-BCs and ES-BCs as compared to their respective undifferentiated counterparts. This file is a PDF document containing Table S20 with the list of all overrepresented GO terms in iPS-BCs compared to iPS cells and in ES-BCs compared to ES cells.