Skip to main content

Qualitative transcriptional signatures for evaluating the maturity degree of pluripotent stem cell-derived cardiomyocytes



Pluripotent stem cell-derived cardiomyocytes (PSC-CMs) are widely used models for regenerative medicine and disease research. However, PSC-CMs are usually immature in morphology and functionality and the maturity of PSC-CMs could not be determined accurately. In order to reasonably interpret the experimental results obtained by PSC-CMs, it is necessary to evaluate the maturity of PSC-CMs and find the key genes related to maturation.


Using the gene expression profiles of normal adult cardiac tissue and embryonic stem cell (ESC) samples, we identified gene pairs with identically relative expression orderings (REOs) within adult cardiac tissue but reversely identical in ESCs. Then, for a PSC-CM model, we calculated the maturity score as the percentage of these gene pairs that exhibit the same REOs in adult cardiac tissue. Lastly, the CellComp method was used to identify the maturation-related genes.


The maturity score increased gradually from 0.8401 for 18-week fetal cardiac tissue to 0.9997 for adult cardiac tissue. For four human PSC-CM models, the mature scores increased with prolonged culture time but were all below 0.8. The genes involved in energy metabolism, angiogenesis, immunity, and proliferation were dysregulated in the 1-year PSC-CMs compared with adult cardiac tissue.


We proposed a qualitative transcriptional signature to score the maturity degree of PSC-CMs. This score can reasonably track the maturity of PSC-CMs and be used to compare different PSC-CM culture methods.


Pluripotent stem cells, including embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSCs), have the potential to differentiate into any cell type including cardiomyocytes (CMs) [1]. Thus, pluripotent stem cell-derived cardiomyocytes (PSC-CMs) are widely used as models for drug screening in vitro [2] or for potential therapeutic applications in vivo [3]. However, the applications of PSC-CMs are limited due to the uncertainty of the maturation status [4,5,6,7]. Therefore, it is necessary to assess whether PSC-CMs are sufficiently mature in culture to make them suitable for their intended purpose.

To date, the maturity of CM is commonly evaluated by morphological and electrophysiological characteristics and cardiac-specific gene markers [8,9,10]. However, the morphological and electrophysiological characteristics are easily affected by subjective factors and thus difficult to quantify as strict maturity markers [11]. Also, there are many problems about cardiac-specific gene markers, such as MYH6, RYR2, KCNJ2, and cTnT [8, 12]. First, these markers are not exclusively expressed in cardiac myocytes. For example, cTnT is also expressed in non-cardiac cells such as smooth muscle cells [13,14,15]. Second, the cardiac-specific genes expressed in the early stage of heart development can also be reactive during a variety of pathophysiologic conditions, including hypoxia, ischemia, hypertrophy, atrophy, diabetes, and hypothyroidism [16,17,18,19]. Third, these markers cannot be used to evaluate the maturity degree of CM.

Given the abovementioned problems, several signatures based on the quantitative measurements of gene expression for evaluating CM maturity have been proposed [11, 20]. However, they usually could not be applied directly to other independent data because of experimental batch effects and data normalization biases [21,22,23]. In contrast, the within-sample relative expression orderings (REOs) of genes, which are qualitative transcriptional characteristics, are insensitive to experimental batch effects. Our previous researches have reported that REOs of gene pairs are highly stable in a type of tissues [24] or cells [25,26,27] in a particular status but tend to be widely disrupted in another status. Based on these unique advantages of REOs, Ao et al. used the hepatocellular carcinoma (HCC)-specific REO patterns of gene pairs to evaluate the similarity between HCC cell lines and HCC tissues [25]. Similarly, it should be possible to integrate heart-related data produced by different laboratories and select the adult cardiac tissue-specific and PSC-specific REO patterns of gene pairs to evaluate the similarity between PSC-CMs and adult cardiac tissue.

In this study, we identified a set of gene pairs with highly stable REOs within adult cardiac tissue but reversely stable REOs in ESCs. Then, for a PSC-CM model, we calculated the maturity score as the percentage of the selected gene pairs that exhibit the stable REOs in adult cardiac tissue. Using this scoring model to evaluate four human PSC-CMs, we found that the maturity scores of PSC-CMs were on the rise with the extension of culture time but were all below 0.8. Even when the culture time was extended to 1 year, the maturity of PSC-CMs reached 0.8274 but the rising degree was not obvious. Finally, using the REO-based CellComp method [27], we identified some key genes and functional pathways which could be responsible for the CM immaturity in vitro.


Data and preprocessing

All gene expression data analyzed in this study were downloaded from the GEO ( [28] and ArrayExpress ( [29] repositories. We identified and screened the gene expression datasets by searching the following keyword: cardiac, heart, cardiomyocyte, embryonic stem cell, pluripotent stem cell, PSC-CMs, and ESC-CMs. Accordingly, a diverse set of cardiac-related gene expression data was collected comprising 17 datasets for human normal adult cardiac tissue, 41 datasets for human ESCs, 2 datasets for fetal cardiac tissue, and 4 datasets for PSC-CMs. The datasets for human normal adult cardiac tissue and ESCs samples were presented in Additional file 1: Table S1. The datasets for human PSC-CMs were described in Table 1. In Additional file 1: Table S2, we have briefly described the cell culture and isolation from tissue for each of the datasets according to the corresponding descriptions in the GEO or ArrayExpress repositories. The human cardiac samples took all the cells from the ventricle, atrium, or septum tissue and were approved by the ethics committee of hospital, local, or school.

Table 1 The datasets for PSC-CM from different culture times

For the raw mRNA expression data (.CEL files) measured by the Affymetrix platform, the Robust Multi-array Average algorithm [30] was used to do background adjustment. For the data measured by the Illumina platform, we directly downloaded the processed data. Each probeset ID was mapped to Entrez gene ID according to the platform file. If a probeset was mapped to multiple or zero gene, then the data of this probeset was deleted. If multiple probesets were mapped to the same gene, the expression value for the gene was defined as the arithmetic mean of the value of multiple probesets. For GSE57338, we divided the 136 normal adult cardiac tissue samples into two parts according to the GSM series numbers of samples: the first 68 samples for establishing the highly stable gene pairs and the remaining 68 samples for evaluating the maturity.

Identifying the gene pairs with stably reversal REOs between adult cardiac tissue and ESCs

For a particular tissue or cell line, pairwise comparisons were performed for all genes to identify gene pairs with stable expression ordering in accumulated samples from different data sources. For each gene pair (Gi, Gj), being viewed as an event with only two possible outcomes (Gi > Gj or Gi < Gj), the gene pairs that the expression level of Gi was higher (or lower) than that of Gj in all accumulated adult cardiac tissue or ESCs samples were defined as highly stable gene pair.

A stable reversal gene pair (Gi and Gj) was selected when its REO, Gi > Gj or Gi < Gj in expression level, was identical in all adult cardiac tissue samples and was reversed (Gi < Gj or Gi > Gj) in all ESCs samples.

Maturity score

The REOs of stably reversal gene pairs in adult cardiac tissue samples were treated as the golden standard. The maturity score in every PSC-CMs sample was calculated as k/n, where n was the number of the stably reversal gene pairs and k was the number of reversal gene pairs with the consistent REOs in the PSC-CMs and adult cardiac tissue. If the REOs of the stable reversal gene pairs in PSC-CMs are closer to that of adult cardiac tissue, the maturity score is closer to 1. For PSC-CMs with two or three technical replicates, we calculated the average maturity score of these PSC-CMs as the final maturity score.

CellComp algorithm

For small-scale cell line data commonly with only two or three technical replicates, the traditional methods such as SAM, RP, and FC lack statistical power or statistical control for differential expression analysis [27]. Thus, we used the REO-based CellComp algorithm which has higher statistical power for small-scale data sets to identify differentially expressed genes (DEGs).

In our study, we adjusted the reference for the highly stable gene pairs from accumulated normal cardiac tissue samples. Then, based on the Fisher exact test, DEGs were identified in PSC-CMs through finding genes with large expression changes that may lead to the stably reversal REOs in two or three PSC-CM replicates compared to the reference, as described in detail in our previous work [27]. The Benjamini–Hochberg method was used to control false discovery rate (FDR) in the multiple tests. The CellComp algorithms are freely available online at

Functional enrichment analysis

Functional enrichment analyses were performed based on Gene Ontology (GO) [31]. First, we used the CellComp algorithm to identify DEGs. Then, the Go-function algorithm [32] was adopted to identify the GO annotations by the following categories: biological processes, molecular functions, and cellular components. The Benjamini–Hochberg method was adopted to estimate the FDR.


The REO-based maturity score

We collected 195 human normal adult cardiac tissue and 116 human ESCs (hESCs) samples from 57 datasets measured by ten platforms (see Additional file 1: Table S1). Then, we identified 16,576,526 and 30,068,382 gene pairs with identical REOs in all adult cardiac tissue and all hESCs samples, respectively. Among these gene pairs, 223,004 gene pairs had reversal REOs between adult cardiac tissue and hESCs. Based on the REOs of 223,004 gene pairs in adult cardiac tissue, we calculated the maturity score of every PSC-CM sample as the percentage of theses gene pairs that showed the same REOs in adult cardiac tissue (see the “Methods” section).

Firstly, we evaluated the maturity degree of adult cardiac tissue and hESCs. The average maturity score of 68 normal adult cardiac tissue samples from GSE57338 was 0.9997. Conversely, the average maturity score of 72 hESCs samples from GSE38662 was 0.0004 (see Additional file 1: Table S3). Similarly, we evaluated the fetal cardiac tissue samples from two datasets and found the maturity score gradually increased from 0.8401 to 0.9665 with the prolonged gestation from 18 to 39 weeks, as showed in detail in Table 2. These results demonstrated that this score can reasonably evaluate the maturity of the cardiac tissue.

Table 2 The maturity scores of fetal cardiac tissue in different gestational ages

Then, we used four human PSC-CM datasets generated by four experimental methods [33,34,35,36] to evaluate the maturity degree of PSC-CMs at different culture time points (Table 1). The result showed that, for each of the four experiments, the maturity score was on the rise with prolonged culture time (Fig. 1). The maturity scores of the four PSC-CMs, which were cultured for 30, 31, 56, and 120 days, reached 0.5371, 0.5943, 0.7529, and 0.7638, respectively. Notably, the maturity scores of PSC-CMs from different sources were all below 0.8. Recent studies reported that long-term cultures may enhance CM maturation [8, 37, 38]. However, for human PSC-CMs cultured for 1 year (GSE62913), the maturity score was only 0.8274 and the rising degree of maturity was not obvious.

Fig 1

The maturity score of four human PSC-CMs from different culture time points

The above results indicated that current differentiation protocols in vitro could not induce mature CMs. Thus, there must be some key molecules that are aberrant in the regulation of CM maturation.

The key genes and pathways influencing the maturation of PSC-CMs

Taking 195 human adult cardiac tissue samples as the reference, we applied the REO-based CellComp method [27] to identify DEGs in PSC-CMs compared to adult cardiac tissue (FDR < 0.05). The CellComp is an efficient tool for analyzing small-scale cell line data with only two or three technical replicates, and it can identify DEGs whose change can disrupt the rank or correlation structure of the transcriptome.

In the 1-year PSC-CMs with three technical replicates (GSE62913), 1323 DEGs (886 up- and 437 downregulated) were detected by the CellComp method. Functional enrichment analysis (FDR < 0.05, Additional file 1: Table S4) showed that the downregulated genes were significantly enriched in the circulatory system, angiogenesis, metabolism, and immune system-related pathways, whereas the upregulated genes were enriched in cell proliferation-related pathways such as nuclear division and sister chromatid segregation.

The top five DEGs, participating in the largest numbers of gene pairs with reversal REOs between the ESCs and the adult cardiac tissue, were showed in Table 3. The five genes were all downregulated in 1-year PSC-CMs relative to adult CMs. It implied that the aberrant expression of these genes could induce great variation of transcriptome structure which can retard the transformation of PSC-CMs from an immature state to mature state. CASQ2, participating in 4987 reversal gene pairs, can facilitate the maturation of Ca2+ handling by modulating RyR activities to coordinate the rates of sarcoplasmic reticulum Ca2+ release and load, which is crucial for CM maturation [39]. CKM is related with mitochondrial functions, including intracellular energy transport and ATP generation, which can promote the oxidative metabolism [40, 41]. It has been reported that, as the CMs mature, mitochondrial oxidative metabolism increases with fatty acid oxidation, providing 90% of the heart’s energy demands [12]. CD36, a cell adhesion molecule, can regulate fatty acid transport [42]. Cell adhesion molecules are related with intercellular communication. Studies pointed out that the cell–cell interaction increasing would promote CM maturation [43, 44]. CMYA5, a Z-disc-related protein, involved in the construction of myocardial fiber [45, 46]. The overexpression of CMYA5 is associated with left ventricular hypertrophy [47, 48]. DCN, a structural constituent of the extracellular matrix, plays roles in proliferation, autophagy, angiogenesis, and fibrin organization [49, 50], and its overexpression can facilitate the formation of blood vessels, which is important for the formation of highly vascularized adult cardiac tissue [51, 52].

Table 3 The top five DEGs participating in the largest numbers of gene pairs with reversal REOs in ESCs compared with the mature adult CMs

The above results indicated that these genes and pathways may be important factors for the CM maturation in the late stage of myocardial differentiation.


In this study, we established a REOs-based score method to assess the maturity degree of PSC-CMs. It is highly robust against large measurement variations introduced by experimental batch effects and platform differences. Essentially, the maturity scores are 1 for adult cardiac tissues and 0 for ESCs, and the scores for the human PSC-CMs increase with the extension of culture time. These results suggested that this score can reasonably evaluate the maturity degree of human PSC-CMs.

Our analysis showed that the maturity of PSC-CMs in vitro can hardly exceed 0.8. Even when PSC-CMs were cultured for 1 year, the maturity score was just 0.8274. Thus, we screened the aberrantly expressed genes in 1-year PSC-CMs compared with adult cardiac tissue samples. The genes involved in energy metabolism, proliferation, angiogenesis, and immunity were dysregulated. We identified some DEGs such as CASQ2, CKM, CD36, CMYA5, and DCN that participated in largest numbers of reversal gene pairs, which could induce the great variation of transcriptome correlation structure and retard the transformation of PSC-CMs from an immature state to mature state. These DEGs might be of special biological significance for CM maturation, because the functionally related genes tend to express coordinately in a stable state but widely altered in another state [53]. Further studies are necessary to determine their roles for PSC-CM maturation.

Many approaches have been reported to improve the maturation of PSC-CMs [43, 54,55,56]. However, there are no definitive markers available to determine which method can better promote the CM maturation. The maturity score can be used to compare different culture methods. For example, the maturity scores for three human PSC-CMs cultured for 14 days from protocol 1 (GSE35671), protocol 2 (GSE64189), and protocol 3 (GSE84815) were 0.3811, 0.5087, and 0.4210, respectively. Moreover, the maturity score of 56-day PSC-CMs from protocol 2 (0.7529) was higher than that of 90-day PSC-CMs from protocol 1 (0.7339). The maturity score of 21-day PSC-CMs from protocol 2 (0.5909) was close to that of 31-day PSC-CMs from protocol 3 (0.5943). Therefore, we can infer that protocol 2 may be more effective to promote CM maturation than protocols 1 and 3.


In summary, we identified a qualitative transcriptional signature which can reasonably track the maturity of PSC-CMs. As a score, it can be used to compare different culture methods for PSC-CM maturation. By REOs-based method, we further identified some genes and pathways which may be important factors for the CM maturation in the late stage of myocardial differentiation.





Differentially expressed genes


Embryonic stem cells


False discovery rate


Hepatocellular carcinoma


Induced pluripotent stem cells


Pluripotent stem cell-derived cardiomyocytes


Relative expression orderings


  1. 1.

    Takahashi K, Tanabe K, Ohnuki M, Narita M, Ichisaka T, Tomoda K, et al. Induction of pluripotent stem cells from adult human fibroblasts by defined factors. Cell. 2007;131:861–72.

    CAS  Article  Google Scholar 

  2. 2.

    Sun N, Yazawa M, Liu J, Han L, Sanchez-Freire V, Abilez OJ, et al. Patient-specific induced pluripotent stem cells as a model for familial dilated cardiomyopathy. Sci Transl Med. 2012;4:130ra147.

    Article  Google Scholar 

  3. 3.

    Menasche P, Vanneaux V, Hagege A, Bel A, Cholley B, Cacciapuoti I, et al. Human embryonic stem cell-derived cardiac progenitors for severe heart failure treatment: first clinical case report. Eur Heart J. 2015;36:2011–7.

    Article  Google Scholar 

  4. 4.

    Tan SH, Ye L. Maturation of pluripotent stem cell-derived cardiomyocytes: a critical step for drug development and cell therapy. J Cardiovasc Transl Res. 2018.

  5. 5.

    Chen HS, Kim C, Mercola M. Electrophysiological challenges of cell-based myocardial repair. Circulation. 2009;120:2496–508.

    Article  Google Scholar 

  6. 6.

    Veerman CC, Kosmidis G, Mummery CL, Casini S, Verkerk AO, Bellin M. Immaturity of human stem-cell-derived cardiomyocytes in culture: fatal flaw or soluble problem? Stem Cells Dev. 2015;24:1035–52.

    CAS  Article  Google Scholar 

  7. 7.

    Bhute VJ, Bao X, Dunn KK, Knutson KR, McCurry EC, Jin G, et al. Metabolomics identifies metabolic markers of maturation in human pluripotent stem cell-derived cardiomyocytes. Theranostics. 2017;7:2078–91.

    CAS  Article  Google Scholar 

  8. 8.

    Lundy SD, Zhu WZ, Regnier M, Laflamme MA. Structural and functional maturation of cardiomyocytes derived from human pluripotent stem cells. Stem Cells Dev. 2013;22:1991–2002.

    CAS  Article  Google Scholar 

  9. 9.

    Feric NT, Radisic M. Maturing human pluripotent stem cell-derived cardiomyocytes in human engineered cardiac tissues. Adv Drug Deliv Rev. 2016;96:110–34.

    CAS  Article  Google Scholar 

  10. 10.

    Shen N, Knopf A, Westendorf C, Kraushaar U, Riedl J, Bauer H, et al. Steps toward maturation of embryonic stem cell-derived cardiomyocytes by defined physical signals. Stem Cell Reports. 2017;9:122–35.

    CAS  Article  Google Scholar 

  11. 11.

    Bedada FB, Chan SS, Metzger SK, Zhang L, Zhang J, Garry DJ, et al. Acquisition of a quantitative, stoichiometrically conserved ratiometric marker of maturation status in stem cell-derived cardiac myocytes. Stem Cell Reports. 2014;3:594–605.

    CAS  Article  Google Scholar 

  12. 12.

    Kuppusamy KT, Jones DC, Sperber H, Madan A, Fischer KA, Rodriguez ML, et al. Let-7 family of microRNA is required for maturation and adult-like metabolism in stem cell-derived cardiomyocytes. Proc Natl Acad Sci U S A. 2015;112:E2785–94.

    CAS  Article  Google Scholar 

  13. 13.

    Lu TY, Lin B, Kim J, Sullivan M, Tobita K, Salama G, et al. Repopulation of decellularized mouse heart with human induced pluripotent stem cell-derived cardiovascular progenitor cells. Nat Commun. 2013;4:2307.

    Article  Google Scholar 

  14. 14.

    Mahmoud AI, Kocabas F, Muralidhar SA, Kimura W, Koura AS, Thet S, et al. Meis1 regulates postnatal cardiomyocyte cell cycle arrest. Nature. 2013;497:249–53.

    CAS  Article  Google Scholar 

  15. 15.

    Xin M, Kim Y, Sutherland LB, Murakami M, Qi X, McAnally J, et al. Hippo pathway effector Yap promotes cardiac regeneration. Proc Natl Acad Sci U S A. 2013;110:13839–44.

    CAS  Article  Google Scholar 

  16. 16.

    Nakao K, Minobe W, Roden R, Bristow MR, Leinwand LA. Myosin heavy chain gene expression in human heart failure. J Clin Invest. 1997;100:2362–70.

    CAS  Article  Google Scholar 

  17. 17.

    Reiser PJ, Portman MA, Ning XH, Schomisch Moravec C. Human cardiac myosin heavy chain isoforms in fetal and failing adult atria and ventricles. Am J Physiol Heart Circ Physiol. 2001;280:H1814–20.

    CAS  Article  Google Scholar 

  18. 18.

    Dhar K, Moulton AM, Rome E, Qiu F, Kittrell J, Raichlin E, et al. Targeted myocardial gene expression in failing hearts by RNA sequencing. J Transl Med. 2016;14:327.

    Article  Google Scholar 

  19. 19.

    Tanaka M, Hiroe M, Ito H, Nishikawa T, Adachi S, Aonuma K, et al. Differential localization of atrial natriuretic peptide and skeletal alpha-actin messenger RNAs in left ventricular myocytes of patients with dilated cardiomyopathy. J Am Coll Cardiol. 1995;26:85–92.

    CAS  Article  Google Scholar 

  20. 20.

    Uosaki H, Cahan P, Lee DI, Wang S, Miyamoto M, Fernandez L, et al. Transcriptional landscape of cardiomyocyte maturation. Cell Rep. 2015;13:1705–16.

    CAS  Article  Google Scholar 

  21. 21.

    Guan Q, Chen R, Yan H, Cai H, Guo Y, Li M, et al. Differential expression analysis for individual cancer samples based on robust within-sample relative gene expression orderings across multiple profiling platforms. Oncotarget. 2016.

  22. 22.

    Nygaard V, Rodland EA, Hovig E. Methods that remove batch effects while retaining group differences may lead to exaggerated confidence in downstream analyses. Biostatistics. 2016;17:29–39.

    PubMed  Google Scholar 

  23. 23.

    Eddy JA, Sung J, Geman D, Price ND. Relative expression analysis for molecular cancer diagnosis and prognosis. Technol Cancer Res Treat. 2010;9:149–59.

    CAS  Article  Google Scholar 

  24. 24.

    Wang H, Sun Q, Zhao W, Qi L, Gu Y, Li P, et al. Individual-level analysis of differential expression of genes and pathways for personalized medicine. Bioinformatics. 2015;31:62–8.

  25. 25.

    Ao L, Guo Y, Song X, Guan Q, Zheng W, Zhang J, et al. Evaluating hepatocellular carcinoma cell lines for tumour samples using within-sample relative expression orderings of genes. Liver Int. 2017;37:1688–96.

    CAS  Article  Google Scholar 

  26. 26.

    He J, Yan H, Cai H, Li X, Guan Q, Zheng W, et al. Statistically controlled identification of differentially expressed genes in one-to-one cell line comparisons of the CMAP database for drug repositioning. J Transl Med. 2017;15:198.

    Article  Google Scholar 

  27. 27.

    Li X, Cai H, Wang X, Ao L, Guo Y, He J, et al. A rank-based algorithm of differential expression analysis for small cell line data with statistical control. Brief Bioinform. 2017. doi:10.1093

  28. 28.

    Clough E, Barrett T. The Gene Expression Omnibus database. Methods Mol Biol. 2016;1418:93–110.

    Article  Google Scholar 

  29. 29.

    Parkinson H, Kapushesky M, Shojatalab M, Abeygunawardena N, Coulson R, Farne A, et al. ArrayExpress--a public database of microarray experiments and gene expression profiles. Nucleic Acids Res. 2007;35:D747–50.

    CAS  Article  Google Scholar 

  30. 30.

    Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003;4:249–64.

    Article  Google Scholar 

  31. 31.

    Harris MA, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, et al. The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res. 2004;32:D258–61.

    CAS  Article  Google Scholar 

  32. 32.

    Wang J, Zhou X, Zhu J, Gu Y, Zhao W, Zou J, et al. GO-function: deriving biologically relevant functions from statistically significant functions. Brief Bioinform. 2012;13:216–27.

    Article  Google Scholar 

  33. 33.

    Babiarz JE, Ravon M, Sridhar S, Ravindran P, Swanson B, Bitter H, et al. Determination of the human cardiomyocyte mRNA and miRNA differentiation network by fine-scale profiling. Stem Cells Dev. 2012;21:1956–65.

    CAS  Article  Google Scholar 

  34. 34.

    Zhang M, Schulte JS, Heinick A, Piccini I, Rao J, Quaranta R, et al. Universal cardiac induction of human pluripotent stem cells in two and three-dimensional formats: implications for in vitro maturation. Stem Cells. 2015;33:1456–69.

    CAS  Article  Google Scholar 

  35. 35.

    Tompkins JD, Jung M, Chen CY, Lin Z, Ye J, Godatha S, et al. Mapping human pluripotent-to-cardiomyocyte differentiation: Methylomes, transcriptomes, and exon DNA methylation “memories”. EBioMedicine. 2016;4:74–85.

    Article  Google Scholar 

  36. 36.

    Nakano H, Minami I, Braas D, Pappoe H, Wu X, Sagadevan A, et al. Glucose inhibits cardiac muscle maturation through nucleotide biosynthesis. Elife. 2017;6.

  37. 37.

    Kamakura T, Makiyama T, Sasaki K, Yoshida Y, Wuriyanghai Y, Chen J, et al. Ultrastructural maturation of human-induced pluripotent stem cell-derived cardiomyocytes in a long-term culture. Circ J Soc. 2013;77:1307–14.

    CAS  Article  Google Scholar 

  38. 38.

    Yang X, Pabon L, Murry CE. Engineering adolescence: maturation of human pluripotent stem cell-derived cardiomyocytes. Circ Res. 2014;114:511–23.

    CAS  Article  Google Scholar 

  39. 39.

    Liu J, Lieu DK, Siu CW, Fu JD, Tse HF, Li RA. Facilitated maturation of Ca2+ handling properties of human embryonic stem cell-derived cardiomyocytes by calsequestrin expression. Am J Physiol Cell Physiol. 2009;297:C152–9.

    CAS  Article  Google Scholar 

  40. 40.

    Gupta A, Akki A, Wang Y, Leppo MK, Chacko VP, Foster DB, et al. Creatine kinase-mediated improvement of function in failing mouse hearts provides causal evidence the failing heart is energy starved. J Clin Invest. 2012;122:291–302.

    CAS  Article  Google Scholar 

  41. 41.

    Chung S, Dzeja PP, Faustino RS, Terzic A. Developmental restructuring of the creatine kinase system integrates mitochondrial energetics with stem cell cardiogenesis. Ann N Y Acad Sci. 2008;1147:254–63.

    CAS  Article  Google Scholar 

  42. 42.

    Yang X, Okamura DM, Lu X, Chen Y, Moorhead J, Varghese Z, et al. CD36 in chronic kidney disease: novel insights and therapeutic opportunities. Nat Rev Nephrol. 2017;13:769–81.

    CAS  Article  Google Scholar 

  43. 43.

    Herron TJ, Rocha AM, Campbell KF, Ponce-Balbuena D, Willis BC, Guerrero-Serna G, et al. Extracellular matrix-mediated maturation of human pluripotent stem cell-derived cardiac monolayer structure and electrophysiological function. Circ Arrhythm Electrophysiol. 2016;9:e003638.

    CAS  Article  Google Scholar 

  44. 44.

    Lee S, Serpooshan V, Tong X, Venkatraman S, Lee M, Lee J, et al. Contractile force generation by 3D hiPSC-derived cardiac tissues is enhanced by rapid establishment of cellular interconnection in matrix with muscle-mimicking stiffness. Biomaterials. 2017;131:111–20.

    CAS  Article  Google Scholar 

  45. 45.

    Poon E, Yan B, Zhang S, Rushing S, Keung W, Ren L, et al. Transcriptome-guided functional analyses reveal novel biological properties and regulatory hierarchy of human embryonic stem cell-derived ventricular cardiomyocytes crucial for maturation. PLoS One. 2013;8:e77784.

    CAS  Article  Google Scholar 

  46. 46.

    Benson MA, Tinsley CL, Blake DJ. Myospryn is a novel binding partner for dysbindin in muscle. J Biol Chem. 2004;279:10450–8.

    CAS  Article  Google Scholar 

  47. 47.

    Kumar S, Mishra A, Srivastava A, Bhatt M, Garg N, Agarwal SK, et al. Role of common sarcomeric gene polymorphisms in genetic susceptibility to left ventricular dysfunction. J Genet. 2016;95:263–72.

    CAS  Article  Google Scholar 

  48. 48.

    Sarparanta J. Biology of myospryn: what’s known? J Muscle Res Cell Motil. 2008;29:177–80.

    CAS  Article  Google Scholar 

  49. 49.

    Jarvinen TA, Prince S. Decorin: a growth factor antagonist for tumor growth inhibition. Biomed Res Int. 2015;2015:654765.

    Article  Google Scholar 

  50. 50.

    Dugan TA, Yang VW, McQuillan DJ, Hook M. Decorin modulates fibrin assembly and structure. J Biol Chem. 2006;281:38208–16.

    CAS  Article  Google Scholar 

  51. 51.

    Li L, Okada H, Takemura G, Kosai K, Kanamori H, Esaki M, et al. Postinfarction gene therapy with adenoviral vector expressing decorin mitigates cardiac remodeling and dysfunction. Am J Physiol Heart Circ Physiol. 2009;297:H1504–13.

    CAS  Article  Google Scholar 

  52. 52.

    Jarvelainen H, Sainio A, Wight TN. Pivotal role for decorin in angiogenesis. Matrix Biol. 2015;43:15–26.

    CAS  Article  Google Scholar 

  53. 53.

    Wang D, Cheng L, Zhang Y, Wu R, Wang M, Gu Y, et al. Extensive up-regulation of gene expression in cancer: the normalised use of microarray data. Mol BioSyst. 2012;8:818–27.

    CAS  Article  Google Scholar 

  54. 54.

    Ulmer BM, Stoehr A, Schulze ML, Patel S, Gucek M, Mannhardt I, et al. Contractile work contributes to maturation of energy metabolism in hiPSC-derived cardiomyocytes. Stem Cell Reports. 2018;10:834–47.

    CAS  Article  Google Scholar 

  55. 55.

    Shadrin IY, Allen BW, Qian Y, Jackman CP, Carlson AL, Juhas ME, et al. Cardiopatch platform enables maturation and scale-up of human pluripotent stem cell-derived engineered heart tissues. Nat Commun. 2017;8:1825.

    Article  Google Scholar 

  56. 56.

    Ronaldson-Bouchard K, Ma SP, Yeager K, Chen T, Song L, Sirabella D, et al. Advanced maturation of human cardiac tissue grown from pluripotent stem cells. Nature. 2018;556:239–43.

    CAS  Article  Google Scholar 

Download references


Not applicable


This work was supported by grants from the National Natural Science Foundation of China (81870242, 81872396, 81572935, 81501203), the Shanghai Committee of Science and Technology (18JC1414300, 16140904200), the Fund of the Key Laboratory of Regenerative Biology of Chinese Academy of Science (KLRB201609) and the Fundamental Research Fund for the Central Universities to Li Li.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its supplementary information files.

Author information




ZG and RC contributed to the conception and design. LP, LL, and JZ contributed financial support; QL and JZ contributed to the collection and assembly of data. JH and YG contributed to the data analysis and interpretation; DW contributed to the production of figure and table. RC and YW did the manuscript writing. ZG and LL approved the final version of the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Zheng Guo or Li Li.

Ethics declarations

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Table S1. The datasets for human normal adult cardiac tissue and ESCs samples. Table S2. The brief description of the cell culture and isolation from cardiac tissue for the datasets. Table S3. The maturity scores for 68 normal adult cardiac tissue and 72 ESC samples. Table S4. The enriched pathways of DEGs for 1-year PSC-CMs. (XLSX 26 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Chen, R., He, J., Wang, Y. et al. Qualitative transcriptional signatures for evaluating the maturity degree of pluripotent stem cell-derived cardiomyocytes. Stem Cell Res Ther 10, 113 (2019).

Download citation


  • Pluripotent stem cell
  • Cardiomyocytes
  • Maturity score
  • Relative expression orderings