Decellularized cardiac extracellular matrix scaffolds with preserved composition and architecture can be used in tissue engineering to reproduce the complex cardiac extracellular matrix. However, evaluating the extent of cardiomyocyte repopulation of decellularized cardiac extracellular matrix scaffolds after recellularization attempts is challenging. Here, we describe a unique combination of biochemical, biomechanical, histological, and physiological parameters for quantifying recellularization efficiency of tissue-engineered cardiac patches compared with native cardiac tissue. Human embryonic stem cell-derived cardiomyocytes were seeded into rat heart atrial and ventricular decellularized cardiac extracellular matrix patches. Confocal and atomic force microscopy showed cell integration within the extracellular matrix basement membrane that was accompanied by restoration of native cardiac tissue passive mechanical properties. Multi-electrode array and immunostaining (connexin 43) were used to determine synchronous field potentials with electrical coupling. Myoglobin content (~60%) and sarcomere length measurement (>45% vs 2D culture) were used to evaluate cardiomyocyte maturation of integrated cells. The combination of these techniques allowed us to demonstrate that as cellularization efficiency improves, cardiomyocytes mature and synchronize electrical activity, and tissue mechanical/biochemical properties improve toward those of native tissue.

Human pluripotent stem cells (PSC) have been used for disease modelling, after differentiation into the desired cell type. Electrophysiologic properties of cardiomyocytes derived from pluripotent stem cells are extensively used to model cardiac arrhythmias, in cardiomyopathies and channelopathies. This requires strict control of the multiple variables that can influence the electrical properties of these cells. In this article, we report the action potential variability of 780 cardiomyocytes derived from pluripotent stem cells obtained from six healthy donors. We analyze the overall distribution of action potential (AP) data, the distribution of action potential data per cell line, per differentiation protocol and batch. This analysis indicates that even using the same cell line and differentiation protocol, the differentiation batch still affects the results. This variability has important implications in modeling arrhythmias and imputing pathogenicity to variants encountered in patients with arrhythmic diseases. We conclude that even when using isogenic cell lines to ascertain pathogenicity to variants associated to arrythmias one should use cardiomyocytes derived from pluripotent stem cells using the same differentiation protocol and batch and pace the cells or use only cells that have very similar spontaneous beat rates. Otherwise, one may find phenotypic variability that is not attributable to pathogenic variants.