Functional arrays of human pluripotent stem cell-derived cardiac microtissues

To accelerate the cardiac drug discovery pipeline, we set out to develop a platform that would be capable of quantifying tissue-level functions such as contractile force and be amenable to standard multiwell-plate manipulations. We report a 96-well-based array of 3D human pluripotent stem cell (hPSC)-derived cardiac microtissues - termed Cardiac MicroRings (CaMiRi) - in custom 3D-print-molded multiwell plates capable of contractile force measurement. Within each well, two elastomeric microcantilevers are situated above a circumferential ramp. The wells are seeded with cell-laden collagen, which, in response to the gradual slope of the circumferential ramp, self-organizes around tip-gated microcantilevers to form contracting CaMiRi. The contractile force exerted by the CaMiRi is measured and calculated using the deflection of the cantilevers. Platform responses were robust and comparable across wells, and we used it to determine an optimal tissue formulation. We validated the contractile force response of CaMiRi using selected cardiotropic compounds with known effects. Additionally, we developed automated protocols for CaMiRi seeding, image acquisition, and analysis to enable the measurement of contractile force with increased throughput. The unique tissue fabrication properties of the platform, and the consequent effects on tissue function, were demonstrated upon adding hPSC-derived epicardial cells to the system. This platform represents an open-source contractile force screening system useful for drug screening and tissue engineering applications.

introduction of cells into these devices, including a tissue transfer step to another device after tissue formation for long-term culture 8,10,[13][14][15] , or a transfer step to a measurement device such as a force transducer 7,9,16,17 to measure contractile force. These steps not only complicate and extend the overall process but may also introduce unwanted variability in baseline tissue conditions, a significant hurdle in designing quality screening applications. Additionally, the ability to measure contractile force, a key functional metric in heart health, is not always easily integrated into these engineered heart tissue screening platforms, limiting the functional assessment and predictive ability of these platforms.
We previously determined key design criteria for the formulation of human pluripotent stem cell (hPSC)-derived cardiac microtissues including an optimal ratio of input cells for tissue formation and induction of intratissue cell alignment 3 . Here we outline the incorporation of these criteria in the design of a 3D printing-based 96-well plate screening platform that supports both cardiac microtissue self-organization and contractile force measurement. Our new platform is simple to manufacture and is seeded with cells that form Cardiac MicroRings (CaMiRi) -cardiac microtissue that is configured for contractile force measurements, in addition to the standard calcium handling and conductance parameters conventionally measured in engineered cardiac tissue 3 . Notably, to maximize accuracy and reproducibility, we have designed a simple process for measuring contractile force in situ with automated data collection using only a bright field camera. Briefly, CaMiRi are formed by seeding cell-laden collagen into a reservoir in each well. The cell-laden collagen forms an organoid around two elastomeric microcantilevers located at the base of each well. The contractions of the CaMiRi cause the microcantilevers to deflect towards the center of the well. The displacement of this deflection can be imaged via a digital camera in each well and used to calculate a total magnitude contractile force, using both the structural and material properties of the elastomer microcantilevers. We optimized this platform by conducting a Central Composite Design (CCD) experiment to test the effects of several input variables on the contractile force exerted by the CaMiRi. We then validated our platform with a cardiotoxicity study using compounds with known effects on cardiomyocyte contractile force response and automated key bottlenecks in the experimental process to enable high throughput measurement of contractile force. Finally, we added fractions of hPSC-derived epicardial cells to our optimized cardiac tissue formulation to test effects on contractile force in CaMiRi. This platform provides a widely accessible solution for functional screening and validation of cardiac-associated drugs.

Development of 3D printing-based method to produce 96-well plates for Cardiac MicroRing (CaMiRi) formation.
We designed a platform that would allow us to culture and measure the contractile forces exerted by cardiac microtissues. The platform consists of wells containing dual polydimethylsiloxane (PDMS) cantilevers around which CaMiRi form due to the compaction of a collagen 1-based matrix by the cardiomyocytes and cardiac fibroblasts (cFB) contained within (Fig. 1A). Our elastomeric microcantilevers provide resistance during auxotonic contraction 18 and serve as sensors for measuring contractile force 1,6,19 . To build this construct, we designed and rendered a 3D Computer Aided Design (3D CAD) of the 96-well tissue culture plate model with outer dimensions and inter-well distances based on Costar ™ 96 well plates to be compatible with conventional multi-well cell culture and imaging equipment, such as multi-well pipettors and automated imaging microscope systems (Fig. 1A,B). This design was printed using an Objet 30 Pro 3D printer to generate a positive mold. Next, we produced a negative PDMS mold that was used to produce a complementary positive PDMS mold. The positive PDMS mold was then used to create a polyurethane-based negative master mold that can generate multiple PDMS plates. This process can be used to create 96-well plates with any desired geometry. The open-source designs for our plates will be made available for researchers to further explore and build customized designs. The final device and micro-cantilevers (Fig. 1C,D) were tested with a force transducer (Microsquisher) to determine a force-displacement standard curve for force measurements. The theoretical curve ( Supplementary  Fig. 1D) for this relationship was determined and confirmed empirically (Fig. 1E).
We designed the well reservoir volumes to precisely match the cell seeding volumes to facilitate the formation of a fully closed ring in the reservoir which is critical for proper tissue formation ( Supplementary Figs 1 and 2). To ensure that any residual debris from the cell seeding step would not interfere with the formation and maturation of the CaMiRi, and to maximize tissue surface area exposed to culture medium (which is important for nutrient diffusion), we designed our wells to promote passive movement of the self-organized CaMiRi away and upward from the point of initial tissue formation and into a toroidal geometry ( Supplementary Fig. 3). To this end, we added a ramp in the middle of the cell seeding reservoir to promote passive tissue movement away from the reservoir during the tissue compaction phase and toward the vertical tip of the cantilevers. To prevent the CaMiRi from slipping off of the cantilevers and to ensure that the fully formed CaMiRi would be positioned at the same vertical location on the cantilevers during contractile force measurement, we incorporated an anchor by angling the cantilever tips away from the well center (Fig. 1C). The consistent positioning of the CaMiRi is critical as the absolute magnitude force measured via cantilever deflection is a function of the vertical positioning of the contracting CaMiRi. The anchor length and angle selected had to be sufficient to trap the tissue yet not too long or acute so as not to survive the demolding process. We screened three angles (30 degrees, 60 degrees, and 90 degrees to the horizontal) and 3 anchor lengths (0.5 mm, 0.75 mm, and 1.0 mm) (Fig. 1D) in four demolding tests and found only the 0.5 mm anchors with either 60 or 45 degree angles remained fully intact. The primary criteria of the resulting mold to remain intact was used to choose the anchor geometry. We chose the 0.5 mm length anchors with the 45-degree (more acute) angle to achieve better tissue anchorage. The complete set of key steps from fabrication to cell seeding is summarized in Supplementary Fig. 4.
Formation of CaMiRi in 96 well plate platforms suitable for contractile force measurements. Next, we developed a seeding protocol for the formation of CaMiRi amenable to high-throughput robotic pipetting systems. First, chilled cell-laden collagen 1-based extracellular matrix is pipetted into the wells www.nature.com/scientificreports www.nature.com/scientificreports/ of the plate ( Supplementary Fig. 3). Within 24 hours, the collagen polymerizes and a CaMiRi forms around the two elastic micro-cantilevers located at the base of each well (Fig. 1F). As the CaMiRi contracts around the micro-cantilevers, the contractile force deflects the cantilevers (Supplementary Video 1). Custom edge-detection software can be used to analyze and calculate a total contractile force based on the deflection measurements, using both the structural and material properties of the micro-cantilevers ( Supplementary Fig. 5).
We tested the effect of varying input cell numbers on CaMiRi formation to test the capability of forming a range of tissue sizes. A preliminary set of cell input numbers ranging from 25,000 to over 150,000 cells per tissue were tested (Fig. 1G) as a proof of concept. The cells in the CaMiRi maintained viability over the course of remodeling conducted after 48 hours (Fig. 1H).
To maximize successful CaMiRi formation and tissue force generation, we next optimized our seeding protocol with respect to cardiomyocyte composition, collagen concentration and seeding density. These factors were assessed using a three-level face centred central composite design (CCD) which allowed us to create a linear regression model to obtain a range of response values for each tested factor. The range of values selected in the CCD for these three factors were based upon previous studies in the formulation and function of cardiac microtissues 1,3,6,20 . Specifically, the CCD explored input cardiomyocyte compositions of 50%, 70%, and 90% (with the remaining tissue composition consisting of cFB acquired from Promocell); collagen concentrations of 1.8 mg/mL, 2.0 mg/mL, and 2.2 mg/mL; and total cells per CaMiRi of 40,000; 70,000; and 100,000 ( Supplementary Fig. 6A,B). A quadratic model was selected to fit the data and an ANOVA was performed to assess the significance of each term as a means of maximizing successful CaMiRi formation and force generation.
Contraction frequency data were used to generate response surface models to determine the optimal formulation that would lead to maximal contractile force and dynamic measurement range ( Supplementary Fig. 7). Fitting a quadratic model to the data demonstrated that the percentage of cardiomyocytes seeded (p < 0.05) and the square of collagen concentration (p < 0.01) were statistically significant. From this model, the condition leading to maximum micro-cantilever deflection was predicted to be 90% cardiomyocytes, 2.03 mg/mL collagen, and 75,400 input cells per well. For practicality, a collagen concentration of 2.0 mg/mL and a cell concentration of 75,000 per well were subsequently used. Whole mount immunostaining of least optimal condition (7), optimal www.nature.com/scientificreports www.nature.com/scientificreports/ condition (10), and predicted optimal condition CaMiRi was performed to observe cardiac troponin T alignment, which was most strongly displayed in the optimal CaMiRi condition ( Fig. 3).
Assessing the CaMiRi platform to measure contractile force response to cardiotropic drugs with known effects. To test the consistency of our platform for measuring contractile force, we studied CaMiRi response to cardiotropic drugs. We first generated CaMiRi based using the optimized seeding conditions determined from our CCD experiment ( Supplementary Fig. 6D). The tissues were cultured for 2 weeks before drug screening assays were conducted. We chose 3 cardiotropic drugs with known effects: Blebbistatin and Nifedipine are negative inotropes which decrease contractile force, and Isoproterenol is a positive inotrope which increases contractile force. We tested 5 doses based on previous studies using these drugs 11,21 . The three drugs were administered, and tissues were incubated for 20 minutes and then assayed relative to a vehicle control containing 0.1% DMSO.
Blebbistatin, a potent inhibitor of myosin, has been previously shown to substantially reduce the force of contraction of explanted rabbit myocardial tissues, with little to no effect on heart pacing 22 . Increasing concentrations of Blebbistatin resulted in diminishing contractile force responses in the CaMiRi (Fig. 4A). Several studies have demonstrated that the administration of Isoproterenol, an adrenoreceptor agonist, causes an increase in the force of contraction of myocardial tissues ex-vivo 23,24 . In contrast to Blebbistatin, administration of Isoproterenol led to a gradual increase in observed contractile force, which peaked at 1 µM (Fig. 4B). Finally, Nifedipine, an L-type Ca2+ channel blocker, caused a decrease in relative contractile force (Fig. 4C) with increasing dose, consistent with results reported in ex vivo myocardial tissues 25 and clinical trials 26 . Overall, our results are consistent with previously reported electrophysiology analysis using the xCELLigence and microelectrode arrays (MEAs) systems on iCell cardiomyocytes 21 . We also used the CaMiRi system to measure the effect of 7 growth factors with reported effects on cardiomyocytes on two electrophysiological parameters of cardiac tissue: maximum capture rate (a measure of tissue response to stimulation) and excitation threshold (a measure of tissue excitability). www.nature.com/scientificreports www.nature.com/scientificreports/ Although there were no significant effects relative to control on maximum capture rate, we observed that both Insulin Growth Factor (IGF)-1 and Histidine-glyco-protein (HRG) elicited significant decreases in excitation threshold relative to control (p < 0.05). These growth factors could thus be the basis for a drug to treat conduction block in patients with problematic electrophysiology. There was no noticeable change in either response to the remaining 5 growth factors ( Supplementary Fig. 8A,B).
Assay automation for high throughput measurement of contractile force. To address experimental bottlenecks, we developed automation protocols to accelerate and standardize CaMiRi seeding, image acquisition, and image analysis. Cell-laden collagen was dispensed in discrete droplets by an Agilent Bravo robotic liquid handler from a source plate to a destination plate ( Supplementary Fig. 9A,B), both maintained at 4 °C to avoid gelation. Droplets were then driven to well bottoms via cooled centrifugation, again at 4 °C ( Supplementary  Fig. 9C). Centrifugation was optimized at 200 g for 5 minutes, both qualitatively by visual inspection of ring formation and functionally by maximizing post deflection ( Supplementary Fig. 9D). Collagen concentration was also increased to 2.2 mg/mL to promote uniform CaMiRi generation based on visual inspection of intact formation of ring structures. Uniform CaMiRi subsequently developed in 96 well format within 7 days ( Supplementary  Fig. 9E).
Automated imaging of CaMiRi in the 96 well PDMS plate presented significant challenges due to slight spacing between wells and minor variability in absolute post heights relative to the plate bottom, both due to the variation involved in plate curation. To overcome these issues, we employed a targeted imaging approach on the microscope in which wells were first imaged at low magnification to identify post coordinates. An image-based autofocus routine was then performed at high magnification, and videos of post edges were acquired at a high frame rate which were used to produce kymographs (Fig. 5A). This system, equipped with environmental control and liquid handling capabilities, also enabled automation of compound addition studies, reducing experiment workflow from five to one hour of hands-on time. Automated image analysis was also employed to accelerate analysis of post deflection videos. Kymographs were generated from video streams and analyzed in a custom MATLAB (Mathworks) program to export beats per minute, contraction amplitude, contraction velocity, and relaxation velocity for each CaMiRi (Fig. 5B). Image analysis software was designed to achieve a high degree of correlation between automated and manual measurements over a wide range of post deflection amplitudes (Fig. 5C). Finally, we validated our automated workflow by confirming positive and negative inotropic responses of isoproterenol and blebbistatin, respectively (Fig. 5D), which were in line with the expected range of effect, along with relative beating frequency (Fig. 5E). Using the CaMiRi system to rapidly explore cell mixing studies on tissue formation. We next sought to engineer a CaMiRi composed of hPSC-derived epicardial and hPSC-CM layers to represent a more anatomically-relevant cardiac tissue 27 . In the mammalian heart, the epicardium is located on the outer region of the heart, facing inward to the pericardium and the pericardial cavity (Fig. 6A). Epicardial cells are present in both the adult, where they give rise to cell types such as the cFB that are critical during both homeostasis and repair, and the developing heart 27 . Epicardial cells secrete factors that promote cardiomyocyte hyperplasia and hypertrophy in the myocardium, which in turn leads to an increase in contractile force. To determine the effects of epicardium on cardiomyocyte contractility in our CaMiRi, we received hPSC-derived epicardial cells using a published 31-day differentiation protocol 27 .
We first labeled hPSC-CM with cell tracker green CMFDA and epicardial cells with cell tracker red CMTPX. We then dosed in epicardial cell fractions (5%, 9%, 20%, and 33%) into our optimized base CaMiRi formulation to determine its functional effects. We also conducted a control set in parallel where additional cFB fractions (5%, 9%, 20%, and 33%) were integrated into the optimized base CaMiRi formulation. The mixed cell formulations were seeded (Fig. 6B upper panel) and after 2 weeks (Fig. 6B bottom panel) CaMiRi formation was observed. We measured contractile force for each condition and observed maximum contractile force in CaMiRi composed of 9% epicardial cell fraction (Fig. 6C). The force generated in CaMiRi composed of 9% epicardial cell fraction was ~2-fold higher than the associated controls. Epicardial cells do not contribute directly to the active contractile force in cardiac tissue 27 , suggesting their interaction with hPSC-CM leads to larger force output of the hPSC-CM or they affect the transmission of force via structural changes in the tissue. As the fraction of both epicardial cells and cFB increased, specifically in the 20% and 33% fractions, contractile force lowered, and instances of tissue failure increased. In the 33% condition, all tissues failed due to tension-induced tears. Optimization of cardiac tissue formulation using epicardial cells and eventually including other stromal cell types will be critical in creating physiologically relevant and highly predictive models. Our platform has been shown to facilitate these studies effectively.
As a final step, we attempted to create a process to build a more complex tissue with the goal of understanding the limits of creating sophisticated models to explore complex biology in our platform. We did this by taking advantage of the well reservoir that becomes vacated by the tissue as it forms and moves along the micro-cantilever to integrate around the primary CaMiRi with a secondary cell type (Supplementary Fig. 10). As a proof-of-principle experiment, we first seeded our platform with Green Fluorescent Protein (GFP)-expressing fibroblasts, which formed a ring around the dual micro-cantilevers, and subsequently seeded Red Fluorescent Protein (RFP)-expressing fibroblasts into the reservoir. The RFP+ cells formed a secondary tissue layer, also in ring form, around the initial GFP+ ring to create a tissue with distinct cell type separation and minimal cell type mixing (Supplementary Fig. 10). This sequential seeding method not only demonstrates that our platform can be used to create more complex tissue structures compared to simultaneous seeding methods, but also provides a model to study tissue-tissue interactions.

Discussion
In this work, we developed an engineered cardiac microtissue platform that enables high-content contractile force measurements in a 96-well plate configuration. Our CCD study identified that the optimal conditions to form CaMiRi involved seeding 75,000 cells composed of 90% CM and 10% cFB in a matrix of 2.0 mg/mL of Collagen 1. It is important to note that the CM:cFB ratio identified here varies from our previous studies 3 . This discrepancy is perhaps due to the nuanced functional differences of the input cell types resulting from the differentiation protocol used to generate the hPSC-CM or primary tissue from which the input cells are derived.
Interestingly, the adult mammalian heart is composed of 20-30% cardiomyocytes and 70-80% non-myocytes 28 . However, in newborn hearts, the percentage of cardiomyocytes are much higher at nearly 70% 29 . The hPSC-CM we use represent an embryonic stage of development, providing a possible explanation for the higher percentages of CM needed for optimal contractile function in our CaMiRi. Perhaps once there are protocols to efficiently mature hPSC-CM to approach the morphology and phenotype that more closely resembles mature heart cells, the optimal ratios in our engineered tissues will also more closely resemble ratios found in native heart tissue. www.nature.com/scientificreports www.nature.com/scientificreports/ Although our CaMiRi exhibited the expected responses to the cardiotropic drugs tested, greater changes in contractile force may be detected if the microtissues we generated possessed a mature electrical pacing apparatus. Interestingly, in a study of five patients, Ross et al. demonstrated that the infusion of isoproterenol, a positive inotrope, causes an increase in heart rate but with inconsistent effects on stroke volume 30 . However, when the heart rate was held constant through electrical stimulation, the infusion of isoproterenol successfully provoked an increase in stroke volume in all patients. Thus, the lack of an even larger increase in measured contractile force in our platform due to isoproterenol treatment may result from a lack of mature pacemaker cells or external electrical stimulation during our assay.
It is important to note that observed effects of drug exposure are also dependent on the assay technology employed and its readout. For example, Micro Electrode Arrays (MEA) have been used to assess changes in electrical field potentials of electrically coupled populations of cardiac cells as a surrogate measure for the changes in the force of contraction. Likewise, the xCELLigence system uses changes in impedance measurements associated with electrically coupled cardiac cells as a surrogate for changes in force of contraction. Using these two assay platforms to assess the effect of isoproterenol on hPSC-CM monolayers, Guo et al. did not observe an increase in calcium flux electrical activity with the MEA system but did observe an increase in impedance on the xCELLigence system 21 . In contrast to both technologies, the CaMiRi platform enables a direct measurement of contractile force. Furthermore, we have developed an automation pipeline to enable high throughput measurement of contractile force in a 96-well format.
Our studies with the epicardial and myocardial mixed tissues show the ability of our system to study cell-cell interactions within a tissue setting. Although our tissues lack spatial organization between the myocardial and epicardial fractions, we showed the possibility of using sequential seeding to study such conditions in a sequential tissue layering CaMiRi model. We envision that this concept will help with tissue studies in many systems www.nature.com/scientificreports www.nature.com/scientificreports/ that involve 2 or more functional cell types that interact either through cell-cell contact, secreted factors, or even through cell migration via chemokines. Additionally, this model provides the capability to study tissue structure in both healthy and disease states using electrophysiological metrics previously determined ex vivo, for example, in the context of ventricular fibrillation 31 .
We are interested in further investigating the effect of the epicardial fraction on contractile force in our CaMiRi platform. The epicardium may be mediating this effect via three potential mechanisms: either by providing structural support for contractile force transmission within the tissue; through the secretion of soluble factors that directly affect the CM fraction within the microtissues; or by providing additional fibroblasts following epithelial-mesenchymal transition. With respect to the second potential mechanism, many factors have been shown to be secreted by epicardial cells including IGF 27 . We observed the positive effect of IGF on cardiac electrophysiological function, specifically on the excitation threshold of cardiac microtissues which is a measure of tissue excitability 32 , which may play a role in this case as well, to increase contractility.
In conclusion, we have validated the advantages and demonstrated the potential of our simple-to-use and rapid-to-manufacture platform for studying engineered heart tissue. We hope this new method sets a new standard for cardiac tissue engineering and contractility testing, along with further possibilities for engineered custom geometry tissue culture plates applied to other tissue types and applications.

Fabrication of cardiac microring (CaMiRi) 96-well plates.
A modified 96-well plate with embedded microcantilevers was designed using Solid Edge Three-Dimensional Computer Aided Design (3D CAD) software (Siemens, Student Edition) and printed with an Objet 30 Pro 3D printer (Stratasys). A polydimethylsiloxane (PDMS) negative was molded from the 3D printed plate, and then a positive PDMS plate molded following it. Next, SmoothCast310 polyurethane (Smooth-On) was used to fabricate a master mold. CaMiRi PDMS devices were fabricated by pouring Sylgard 184 (Dow Chemicals) into the polyurethane mold, followed by degassing to remove bubbles and curing for 48 hours at 37 °C. The resulting PDMS CaMiRi device was then carefully removed from the mould and visually inspected for defects with a Leica MZ6 stereomicroscope.

Cardiac differentiation of human pluripotent stem cells. Cardiac differentiation of HES2 (ES Cell
International) human embryonic stem cells (hESC) was carried out as previously reported 3 . The hESC were maintained and expanded as described previously. Briefly, HES2 cells were cultured on mouse embryonic feeders (MEF) for 6 days in HES2 maintenance media (80% DMEM/F12, 20% KOSR, 20 ng/mL bFGF, 0.5% P/S, 1% NEAA, 1% BME) and media was changed daily. Cells are maintained in normoxia at 37 °C in a 5.0% CO 2 atmosphere. The cells were then trypsinized along with MEFs and plated onto Matrigel (diluted at 1:30) coated plates at a split ratio of 1:3 for MEF depletion. After two days of MEF depletion, HES2 cells were again trypsinized and seeded into AggreWells manufactured in-house to form hEB. The hEB were generated using 400 µm microwell PDMS inserts cast from a silicon master mould. The inserts were cut and glued into 24-well tissue culture plates and then sterilized using ethanol. The microwells were then coated with 5% Pluronic Acid for at least an hour and washed with PBS before cell seeding. A single cell suspension of aggregation media containing base media and T0 cytokines supplemented with ROCK inhibitor Y-27632 was then seeded into the wells and allowed to aggregate overnight after centrifuging at 200 × g. Cells are maintained in hypoxia at 37 °C in a 5.0% CO 2 and 5.0% O 2 atmosphere. After 24 hours, hEB were formed and aggregation media was exchanged for T1 media. On day 4, hEB were removed from AggreWells and placed in Low cluster 6-well plates (NUNC). Corresponding media for T4, T8, T12 was freshly made and exchanged. On T12, cells were returned to normoxia at 37 °C in a 5.0% CO 2 atmosphere. Media was replaced every 8 days onward.
For select experiments, cardiomyocytes were differentiated from human induced pluripotent stem cells (iPSC; ThermoFisher) using modifications to the monolayer-based differentiation protocol carried out by Lian et al. and licensed from the Wisconsin Alumni Research Foundation 33 . Briefly, iPSCs were seeded on PLO-laminin (Sigma; ThermoFisher) coated tissue culture plates at ~40,000 cells/cm 2 in mTeSR1 (StemCell Technologies) until confluency reached 85-95%. Cells were then treated with 12 µM CHIR 99021 (Tocris) in RPMI-B27 basal medium (ThermoFisher) for 22-26 hours, before CHIR was replaced by RPMI/B27 basal medium. After 48 hours in basal medium, cells were treated with 5 µM IWP2 (Tocris) and a potent and proprietary selective tankyrase inhibitor (IC 50 < 10 nM in DLD-1 epithelial cells) for 48 hours in RPMI/B27 basal medium. IWP2 and tankyrase inhibitor were then removed and cells were grown in RPMI/B27 basal medium with insulin with medium change every other day until harvest. Cells were harvested on day 11 to day 15 post-differentiation using TrypLE (ThermoFisher). The percentage of cTnT-positive cells were determined using flow cytometry as reported previously 33 .
IPSC-CM were thawed 7 days prior to CaMiRi seeding. T75 flasks (Nalgene) were coated with 25 µg/mL Synthemax II-SC (Corning) overnight at 37 °C. Synthemax solution was then aspirated immediately prior to cell seeding. Cryopreserved iPSC-CM were thawed in at 37 °C water bath for 4 minutes then transferred to a 50 mL conical tube containing 5 mL of growth media. Cells were centrifuged at 220 × g for 3 minutes and cell pellet was resuspended in 10 mL growth media per T75 flask. After 48 hours, a half media exchange was performed; thereafter, media was fully replaced every 48-72 hours.
Seeding and cultivation of CaMiRi. hPSC-CM were suspended in a collagen mastermix and seeded into cardiac microtissue wells at desired density (optimal condition determined was 90% cardiomyocytes, a collagen concentration of 2.0 mg/mL, and a total cell concentration of 75,000 per microtissue). Microwell substrates were prepared by sterilizing with ethanol, washing and coating with 5% Pluronic Acid for at least an hour each. During coating, hESC-or iPSC-CM and hCF were prepared. Aggregates from hESC-CM differentiation were incubated in Collagenase for 3 hours at 37 °C, 5% CO 2 . Aggregates were then immersed in wash solution (10% FBS and 90% DMEM/F12) and triturated ~10 times. Once aggregates were dissociated into single cells, samples were counted. iPSC-derived CM and hCF were dissociated from T75 flasks by first washing with phosphate-buffered saline without calcium nor magnesium (PBS −/− ) followed by a brief rinse with TrypLE (ThermoFisher); cells were then incubated for 10-20 minutes at 37 °C until dissociated and counted. The collagen mastermix was prepared by combining the following: 10X M199 (GIBCO), Glutamax (GIBCO), Collagen 1 (3.66 mg/mL) (BD), Glucose (0.3 g/mL) (GIBCO), NaOH (SIGMA), NaHCO3 (0.075 g/mL) (SIGMA), Hepes (GIBCO), GFR Matrigel (BD), ddH20 at appropriate ratios for desired collagen concentrations. Collagen mastermix was maintained on ice under 4 °C to prevent premature crosslinking. Finally, 12 uL of mastermix was pipetted into each well (of 96-well plate), after which the entire plate was placed into a normoxic incubator for 15 minutes. After 15 minutes, 300 µL of cell culture media was slowly added to each well so as to not disrupt the polymerized collagen microtissues. Media was exchanged every 2 (iPSC-CM) or 4 (hESC-CM) days. CaMiRi remodeled between 1-3 days depending on input cell composition. Imaging of CaMiRi was performed in-situ. Samples were fixed, permeablized, and stained inside the microwells and imaged using a fluorescence microscope.
Automated seeding of CaMiRi. A Bravo automated liquid handling platform (Agilent) equipped with a 9-pod deck was utilized for automated seeding of CaMiRi. 2 pods were outfitted with Thermoshake (Inheco) units maintained at 4 °C in order to keep source and destination plates cooled to prevent premature gelation of mastermix. Mastermix, containing cells, was evenly distributed to wells E9:H12 of a 96 well U-bottom plate (Falcon). A custom protocol was then employed to distribute 4×2.5 µL drops of mastermix from the source plate to the north, east, south, and west positions of 16 destination wells at a time; this was repeated 6 times to seed the entire plate and helped minimize dead volume. The destination plate was then centrifuged at 400 x g for 5 minutes at 4 °C to drive the mastermix down into a continuous ring at the bottom of each well. Plate was then incubated for 30 minutes 37 °C, after which 300 µL of maintenance media, comprised of DMEM/F12 (Life Technologies), 2% FBS (Sigma), 10 µg/mL insulin (Life Technologies), and 1x antibiotic/antimycotic (ThermoFisher), was slowly added to all wells. Plate was returned to 37 °C, 5% CO 2 incubator.
Video acquisition and analysis. Throughout the 14-day culture period, images of CaMiRi were acquired every 2 to 3 days using a 4× objective on an Olympus CKX41 microscope equipped with an OptixCam Summit OCS-3.0 digital microscope camera. On day 14, 20 s long videos were acquired using 4× and 10× objectives for each CaMiRi well. Cantilever deflection and contraction frequency were determined by direct measurement of cantilever deflection within video frames via the ImageJ digital ruler function or segregating the 10× videos into images and importing them into a custom-designed cantilever tracking software developed in collaboration with CellScale (Waterloo, Canada).

Automated video acquisition and analysis. Automated videos of CaMiRi-induced post deflections in
96 well plates were acquired at 20x on an Image Xpress Micro 4 (Molecular Devices) by employing a custom journaling solution in Metamorph (Molecular Devices). CaMiRi were maintained at 37 °C and 5% CO 2 during experiments. Briefly, wells were initially imaged at 2x and image analysis was used to identify the x-and y-coordinates of the left post in each well (previous experiments supported the need to only image a single post due to like measurements from both posts; data not shown); magnification was then increased to 20× and a contrast-maximizing autofocus routine was used to focus on the top edge of each post. 20 s videos were then acquired at 100 frames per second using a limited region of interest that captured post displacement. For experiments involving compound addition, compounds were prepared at 5x in a separate U-bottom 96 well plate (Falcon) and maintained on the system's incubated stage at 37 °C. Following acquisition of baseline post deflection, 50 µL of compound was added from source well (one well at a time) to destination well containing 200 µL of medium. Compounds were incubated with CaMiRi for 2 minutes prior to acquisition of an additional 20 s, post-compound addition video.
To extract contractility measurements, videos were converted to 2-dimensional kymographs within Metamorph by projecting each video's center row of pixels over time in the y-direction. Kymographs were then analysed with a custom MATLAB (Mathworks) script to extract beats per minute, contraction amplitude, contraction velocity, and relaxation velocity.
Immunohistochemistry and confocal imaging. CaMiRi were rinsed with phosphate buffered saline (PBS) and fixed with 2% paraformaldehyde (PFA) overnight at 4 °C. They were then rinsed again with PBS and permeabilized in 0.1% Triton-X in blocking solution (2% BSA in PBS) for 5 min. Microtissues were then incubated in blocking solution for 30 min. Primary antibody was diluted in a 1% BSA in PBS solution and was added and left for 3 days at 4 °C. Mouse IgG1 anti-cardiac troponin (Thermofisher MS-295-P) was used at a dilution of 1:100, while Rabbit IgG anti-Vimentin (Abcam AB16700) was used at a dilution of 1:500. After three Scientific RepoRtS | (2020) 10:6919 | https://doi.org/10.1038/s41598-020-62955-3 www.nature.com/scientificreports www.nature.com/scientificreports/ consecutive washes in PBS of 1 hour each, secondary antibodies were applied for 1 day at 4 °C. Alexa Fluor ® 488 goat anti-mouse IgG1 and Alexa Fluor ® 555 goat anti-rabbit (Life Technologies) antibodies were used at a dilution of 1:1000 for cardiac troponin and vimentin respectively. After the removal of the secondary antibodies, nuclear staining was performed using a DAPI solution (1 µg/mL in PBS) for 15 min. The microtissues were then washed in PBS (3 × 1 hour) before mounting on glass slides using fluoro-safe mounting media (Dako, S3023). Mounted samples were allowed to cure for 1 day and then imaged using a confocal microscope (Zeiss LSM700 Confocal Microscope).
Statistical analysis and data representation. Statistical analysis for the Central Composite Design (CCD) experiment was performed with Design Expert v.9 software (Stat-Ease Inc.). The CCD was used to determine the optimal tissue formulation for CaMiRi based upon 3 parameters of collagen concentration, total cells per tissue, and percentage of cardiomyocytes and cFB in the tissue. The two factors used as metrics were contraction frequency and contractile force (as measured by cantilever deflection). The design consisted of 64 experimental runs with four replicates at the factorial and axial points, and eight replicates at the centre point ( Supplementary Fig. 7A). The experimental design and data analysis were performed using Design Expert v.9 (Stat-Ease Inc.). Statistical analyses between groups were computed by an ANOVA using PSPP v.0.8.3 software, with pairwise comparisons performed by Tukey's post-hoc test. A quadratic model was selected to fit the data and an ANOVA was performed to assess the significance of each term as a means of maximizing successful wire formation and cantilever deflection. Statistical analysis for the drug screen experiment was conducted with pairwise comparisons performed by the Mann-Whitney U test. All error bars represent the standard deviation of three or more biological replicates. Asterisks (*) indicate statistical significance between conditions of P < 0.05. All data analyses, including graphical representations, were performed using Excel (Microsoft, Redmond, WA).