Abstract
Paucity of physiologically relevant cardiac models has limited the widespread application of human induced pluripotent stem cell (hiPSC)-derived cardiomyocytes in drug development. Here, we performed comprehensive characterization of hiPSC-derived cardiomyocyte subtypes from 2D and 3D cultures and established a novel 3D model to study impulse initiation and propagation. Directed differentiation approaches were used to generate sinoatrial nodal (SANCM), atrial (ACM) and ventricular cardiomyocytes (VCM). Single cell RNA sequencing established that the protocols yield distinct cell populations in line with expected identities, which was also confirmed by electrophysiological characterization. In 3D EHT cultures of all subtypes, we observed prominent expression of stretch-responsive genes such as NPPA. Response to rate modulating drugs noradrenaline, carbachol and ivabradine were comparable in single cells and EHTs. Differences in the speed of impulse propagation between the subtypes were more pronounced in EHTs compared with 2D monolayers owing to a progressive increase in conduction velocities in atrial and ventricular cardiomyocytes, in line with a more mature phenotype. In a novel binary EHT model of pacemaker-atrial interface, the SANCM end of the tissue consistently paced the EHTs under baseline conditions, which was inhibited by ivabradine. Taken together, our data provide comprehensive insights into molecular and electrophysiological properties of hiPSC-derived cardiomyocyte subtypes, facilitating the creation of next generation composite cardiac models for drug discovery, disease modeling and cell-based regenerative therapies.