A Microwell Cell Capture Device Reveals Variable Response to Dobutamine in Isolated Cardiomyocytes
J. Alexander Clark,1 Jonathan D. Weiss,1 and Stuart G. Campbell1,*
ABSTRACT Isolated ventricular cardiomyocytes exhibit substantial cell-to-cell variability, even when obtained from the same small volume of myocardium. In this study, we investigated the possibility that cardiomyocyte responses to b-adrenergic stim- ulus are also highly heterogeneous. To achieve the throughput and measurement duration desired for these experiments, we designed and validated a novel microwell system that immobilizes and uniformly orients isolated adult cardiomyocytes. In this configuration, detailed drug responses of dozens of cells can be followed for intervals exceeding 1 h. At the conclusion of an experiment, specific cells can also be harvested via a precision aspirator for single-cell gene expression profiling. Using this system, we followed changes in Ca2+ signaling and contractility of individual cells under sustained application of either do- butamine or omecamtiv mecarbil. Both compounds increased average cardiomyocyte contractility over the course of an hour, but responses of individual cells to dobutamine were significantly more variable. Surprisingly, some dobutamine-treated cardi- omyocytes augmented Ca2+ release without increasing contractility. Other cells responded with increased contractility despite unchanged Ca2+ release. Single-cell gene expression analysis revealed significant co-expression of b-adrenergic pathway genes PKA regulatory subunit type I, PKA regulatory subunit type II, and Ca2+/calmodulin-dependent protein kinase II across cardiomyocytes. Other data supported a connection between the effects of dobutamine on relaxation rate and the expression of protein phosphatase 2. These findings suggest that variable drug responses among cells are not merely experimental artifacts. By enabling direct comparison of the functional behavior of an individual cell and the genes it expresses, this new system con- stitutes a unique tool for interrogating cardiomyocyte drug responses and discovering the genes that modulate them.
SIGNIFICANCE We have created a microwell capture device that allows drug responses of dozens of isolated adult cardiomyocytes to be monitored for extended intervals. Using the device, we observed striking diversity in the Ca2+ handling and contractility responses of each cell to a b-adrenergic agonist. This included cells that responded to dobutamine by doubling the amplitude of Ca2+ release while decreasing contractility and vice versa. Our data suggest that these diverse responses are related to gene expression differences between cells. This work demonstrates the feasibility of linking the drug responses of individual cells with their gene expression. It opens the possibility of exploiting cell-to-cell variation to discover new genes that participate in and modulate regulatory cascades.
INTRODUCTION
Individual isolated cardiomyocytes display an astonishing de- gree of functional heterogeneity. We recently reported that cells isolated from the same myocardial region vary in their intrinsic relaxation rate by more than 25% of the population mean (1). At least part of this variability can be explained bytroponin I (TnI) phosphorylation, which was found to be elevated in the most rapidly relaxing cells. Because TnI phos- phorylation is a known end point of the cardiac b-adrenergic signaling pathway, our observations provide implicit evidence of cell-to-cell variation in this regulatory cascade.Seeking direct evidence to support this hypothesis, we studied single-cell responses to dobutamine, a b-adrenergic agonist and potent positive inotrope often used in clinical settings (2). Omecamtiv mecarbil (OM), another positive inotrope (3), was selected as an adrenoceptor-independent control. Unlike dobutamine, which ultimately affectsmany intracellular targets, OM is a direct small-molecule activator of cardiac myosin with high molecular specificity. We therefore expected cardiomyocyte contraction to respond more uniformly to OM than to dobutamine.To facilitate these studies, we designed and validated a novel device that traps cardiomyocytes into a grid of micropat- terned wells. These wells provide a uniform orientation and consistent positioning of isolated cardiomyocytes, even during fluid flow. The device enables the rapid and robust measure- ment of dozens of cells in the same experiment, with the ability to track cells during 1 h or more of drug exposure.By coupling this microwell array with a precision aspi- rator (1), it was also possible to collect specific cells after an experiment and perform single-cell gene expression assays. To our knowledge, the experiments reported here constitute the first successful attempts to relate a single car- diomyocyte drug response to gene expression in the same cell. The resulting data expose surprising cell-to-cell varia- tion in dobutamine responses and that these variations are due to stochastic gene expression.
MATERIALS AND METHODS
All animal procedures were approved by the Yale University Institutional Animal Care and Use Committee (approval number 2015-11528), compliant with the regulations of the Animal Welfare Act, Public Health Service, and the United States Department of Agriculture. Studies were performed using seven female Sprague Dawley rats, comprising retired breeders 4–6 months of age purchased from Charles River Laboratories (Wilmington, MA). Animals were housed under a standard 12:12 h light/dark cycle and fed ad libitum in accordance with an approved Yale University protocol. Animals were exposed to 15 min of 500 mL min—1 isoflurane inhalation and then were subjected to a terminal bilateral thoracotomy for removal of the heart. All experiments were conducted in accordance with the relevant guidelines established standards.
Cardiomyocyte isolation
Left ventricular adult rat cardiomyocytes were isolated via Langendorff perfusion of the heart with a collagenase-containing solution. Adult female rats were anesthetized with isoflurane for 15 min, injected with 0.5 mL of —1 digested through mechanical agitation on a shaking incubator set at 115 rpm for 10 min at 37◦C and then gently triturated to liberate individual cells from the tissue. Undigested tissue chunks were then transferred to fresh digestion buffer in a separate microcentrifuge tube, and the process was repeated as many as three additional times or until all tissue had been di- gested. Cells were removed from digestion buffer by centrifugation and resus- pended in a sequence of perfusion buffers supplemented with fetal bovine serum and with gradually increasing calcium concentrations (0.05– 1.1 mM). Functional experiments were performed after allowing cells to rest for at least 1 h. Unless otherwise stated, each experiment utilized cells isolated from the same small volume of myocardium within a single rat.
Microwell device and transcriptomics overview
To measure sarcomere-length changes and Ca2+ transient records from Fura-2-loaded cardiomyocytes consistently over multiple time points, we designed a custom microwell polydimethylsiloxane (PDMS) substrate within a three-dimensional printed cell bath that kept the cells uniformly oriented (Fig. 1 A). An array of rectangular posts allowed cells to be exposed to superfused physiological buffer and electrical stimulus from platinum pacing electrodes while preventing them from drifting away from the previous measurement window (Fig. 1, B and C). After taking functional records, some of the cardiomyocytes were individually aspirated and deposited via an automated single-cell aspirator into microcentrifuge wells filled with lysis buffer for transcriptomic analysis.
PDMS microwell manufacturing
SYLGARD 184 Silicone Elastomer Base and Curing Agent (Dow Chemi- cal, Midland, MI) were mixed thoroughly in a 8:1 mass ratio and vacuum desiccated for 20 min until all air bubbles were removed. The resulting so- lution was spin coated (Laurell Technologies, North Wales, PA) on a SU8 silicon master wafer (soft lithography negative photoresist) containing the microwell pattern at 400 rpm for 60 s, resulting in a uniform ~0.3 mm PDMS thickness. The wafer was baked in a 70◦C oven for at least 4 h. Wells were cut using a custom razor blade stencil on a hard steel plate, and both sides were spray washed with ethanol and then air dried. The device had dimensions of 42 × 34 × 5 mm, and a 20 × 6.8 × 0.3 mm PDMS microwell array was attached as the floor of the main bath using liquid-tight silicon adhesive and a glass coverslip. Two channels for 0.5-mm diameter platinum wire electrodes (95% platinum, 5% ruthenium; Hoover & Strong, North Chesterfield, VA) ran the length of the bath and were used to field stimulate cells at 0.5 or 1 Hz. A temperature probe and four 3U power resistors along both sides of the bath were used to maintain an internal bath temperature of 36 5 1◦C. Inflow and outflow channels with barbed fittings allowed for superfusion with Tyrode’s solution driven by external syringe pumps at 0.4 mL min—1. At this rate, the bath volume was exchanged every 2.3 min. of the heart. Excised hearts were cannulated via the aorta and mounted onto a Langendorff perfusion apparatus within 3 min of excising the heart. Hearts were immediately perfused for 15 min with 37◦C calcium-free perfu- sion buffer containing (in mM) 118 NaCl, 4.8 KCl, 1.25 KH2PO4, 1.25 MgSO4, 10 2,3-butanedione monoxime, 25 Hepes, and 11 glucose, supple- mented with 10 mM (5)-propranolol (hydrochloride) (Cayman Chemical, Ann Arbor, MI) (pH 7.3). A digestion buffer was next perfused through the heart for 20 min. The digestion buffer consisted of perfusion buffer supple- mented with (in mM) 2.5 carnitine, 5 taurine, 2.5 glutamic acid, 0.025 CaCl2, 120 U mL—1 collagenase type II (Worthington Biochemical, Lake- wood, NJ), and 0.96 U mL—1 protease (Sigma-Aldrich, St. Louis, MO) (pH 7.3). The heart was then cut down from the apparatus, and the left ven- tricular free wall was removed. The subepicardial portion of this left ventric- ular excision was cut free, and the innermost third of that subepicardial layer was isolated to avoid edge effects from the base and apex. The resulting strip was segmented into three equal regions of ~1 × 1 × 1 mm (length × width × thickness) in size per region.
Device chassis preparation
The chassis of the device was printed using a resin-cure lithographic printer (Form 2; Formlabs, Somerville, MA) and washed in an isopropanol bath for 20 min. Compressed air was used to clear all channels to prevent debris and resin buildup, and the device was allowed to dry for at least 1 h. After the coverslip with PDMS microwells was secured onto the bottom of the de- vice, the device was plasma treated with oxygen for 90 s under a 200 mtorr vacuum in preparation for open-faced cell seeding.
Functional testing of cardiomyocytes
Cardiomyocytes were imaged on a custom microscope apparatus in Tyrode’s solution (containing in mM: 140 NaCl, 5.4 KCl, 1.8 CaCl2, 1 MgCl2, 25 Hepes, and 10 glucose; pH 7.3). Before study, cells were incubated in
Single-cardiomyocyte transcriptomics
FIGURE 1 Overview of experimental apparatus and PDMS microwell device for cell capture and transcriptomics. (A) Schematic of device with optical data collection system is shown below. Arrows indicate direction of fluid flow. Device contains platinum pacing electrodes, heating elements, and a temper- ature sensor. (B) Enlarged diagram of PDMS microwell pattern is shown. Each filled rectangle is a pillar between which cells are captured. (C) Bright-field micrograph (left) of a captured cardiomyocyte is outlined in blue, under 30× magnification. Cells can be collected with an automated single-cell aspirator (right) for single cardiomyocyte transcriptomic analysis. To see this figure in color, go online darkness for 10 min with Tyrode’s solution supplemented with 2.2 mM Fura- 2AM (MilliporeSigma, Burlington, MA) and with Pluronic F-127 (0.022%; Sigma-Aldrich) for Ca2+ fluorescence imaging (except for dye-free cardio- myocytes in Fig. 3). After loading, cells were resuspended in fresh Tyrode’s solution and allowed to settle naturally until imaging. The cells were vor- texed and resuspended to a concentration of 25 cells/mL. 30 mL of cell-laden solution was added to the device (already containing 100 mL Tyrode’s solu- tion, perfusion halted). After 30 s, perfusion was introduced to wash out cells that did not settle into the wells. Cardiomyocyte Ca2+ transients and un- loaded shortening contractions were measured using an inverted microscope (Eclipse Ti-U; Nikon, Tokyo, Japan) equipped with either the traditional or PDMS microwell perfusion bath.
The baths were temperature controlled (Cell MicroControls, Norfolk, VA) and continuously perfused with 36 5 1◦C Tyrode’s solution. Cells were field stimulated at a constant 35 mA at 1 Hz for comparison between traditional and PDMS microwells and at Redmond, WA) and processed with custom software written in MATLAB (The MathWorks, Natick, MA). Only rod-shaped cardiomyocytes with well-defined sarcomere striations were measured. Cardiomyocytes that did not respond to pacing were excluded regardless of appearance. From each recording, peak sarcomere-length shortening, time to peak shortening (TTP), time from peak shortening to 50% relengthening (RT50), diastolic sarcomere length, magnitude of Ca2+ transient, and Ca2+ decay constant were calculated.
Multiple cells were studied during the same experiment by translating the computer-controlled microscope stage to marked cell locations, visiting each cell in turn. This route was then traversed repeatedly during the exper- iment to monitor cell responses over time. The throughput of the device, in terms of the number of cells that can be followed in a single experiment (ncells), is determined by the following relationship: t 0.5 Hz for all other experiments to reduce effects from experimental rundown. Cells were preconditioned under electrical stimulus and flow for >12 min before collecting functional records.Contractile events were imaged in real time using a sarcomere-length camera system at an acquisition rate of 2 kHz (HVSL; Aurora Scientific, Aurora, ON, Canada). Ca2+ transient measurements were recorded by using a low-pass filtered illumination with oscillating excitation wavelengths of 340 and 380 nm controlled by a RatioMaster fluorescence system (Horiba, Kyoto, Japan) with wavelengths switched every 10 ms. Fura-2 photobleach- ing over many repeat excitation cycles was compensated for using an empirically derived linear correction.
This involved fitting data from un- treated, Fura-2-loaded cells with a linear correction coefficient based on time and repetition of measurements. The correction coefficient was used to offset the bleaching of the fluorescent signal. Data signals were re- corded with a DAP5216a data acquisition system (Microstar Laboratories, where tperiod is the desired period between repeat measurements of the same cell, tacquisition is the time required to translate the stage to the next cell and perform a single acquisition (15 s under the stated conditions), and d is the dropout rate or probability that a cell will become unsuitable for measure- ment during the experiment (~0.25 in our studies). The relation described by Eq. 1 makes evident the trade-off between the total number of cells fol- lowed and how frequently the behavior of each cell is sampled. For instance, if temporal resolution of one measurement every minute is desired, the number of cells that can be followed is limited to three. In most experiments performed in this study, we considered 5 min to be suf- ficient temporal resolution (tperiod = 300 s), meaning that roughly 20 cells could be followed in parallel throughout a long experiment.
Drug treatments
Cells in the drug comparison experiments were preconditioned for >12 min before perfusion was switched to either a drug-containing or control solu-
tion. The ‘‘Control’’ group of cells was perfused with Tyrode’s solution only. The ‘‘DMSO’’ group of cells was treated with Tyrode’s solution sup- plemented with 0.1% dimethyl sulfoxide (DMSO), the same concentration as the vehicle in the drug-containing solutions. The ‘‘Dobutamine’’ group of cells was treated with Tyrode’s solution supplemented with 40 nM of dobut- amine (hydrochloride) (Cayman Chemical) carried in 0.1% DMSO. The ‘‘OM’’ group of cells was treated with Tyrode’s solution supplemented with 150 nM CK-1827452 (Cayman Chemical) carried in 0.1% DMSO. Each group of cells was constantly superfused with each respective solution after preconditioning until the completion of the experiment.
Single-cell RT-qPCR
After completing functional testing, individual cardiomyocytes were imme- diately aspirated and deposited into separate microcentrifuge tubes contain- ing SingleShot Cell Lysis buffer (Bio-Rad Laboratories, Hercules, CA) supplemented with DNase and proteinase K. Aliquots were frozen at
—80◦C overnight to aid in cell lysis and then heated to 37◦C for 5 min to digest genomic DNA, followed by an increase in temperature to 75◦C for 5 min to inactivate DNase. Complementary DNA (cDNA) was synthesized using the iScript Advanced cDNA Synthesis Kit (Bio-Rad Laboratories) and then preamplified using SsoAdvanced PreAmp Supermix (Bio-Rad Laboratories) to amplify target cDNAyield. A nested primer design was im- plemented to reduce nonspecific binding. qPCR was performed using SsoAdvanced Universal SYBR Green Supermix on a CFX Connect Detec- tion System (Bio-Rad Laboratories). Genes of interest were Camk2d, Gapdh, Ppp1ca, Ppp2cb, protein kinase A gene, Prkar1a, and Prkar2a. Transcript abundance was normalized to Gapdh and then to the respective transcript abundance in control cells. Gene expression FCs were hierar- chically clustered with functional properties to generate a heatmap of —log2 (FC) values. All primer sequences are listed in Table S2.
Statistical analysis
Data were analyzed using Prism (GraphPad Software, San Diego, CA), Minitab (Minitab, State College, PA), and ClustVis (BIIT Research Group, University of Tartu, Tartu, Estonia (4)). All data were presented as mean 5 95% confidence interval. Comparisons were conducted via one-way anal- ysis of variance test followed by multiple comparisons (Kurskal-Wallis post hoc analysis) or via an unpaired two-tailed Student’s t-test. A linear regression was applied to assess the correlation between variables, using an F-test to determine whether the slope of the fit line differed from 0. Sig-
nificance was defined by p < 0.05. For comparison of variance, Levene’s test (assuming non-Gaussian distributions) was used.
RESULTS
To assess whether the PDMS microwell array had any sig- nificant effect on the integrity of the cardiomyocytes, sarco- mere shortening dynamics, or Ca2+ transient measurements, we compared the PDMS microwell device (n = 57 cells) to a more traditional setup in which cardiomyocytes are depos- ited directly onto a coverglass that forms the bottom of the superfusion bath (n = 58 cells). After allowing cell ali- quots to precondition for at least 12 min under 1 Hz pacing, one set of sarcomere-length shortening and Ca2+ transient records were collected from each cell within the time frame of 45 min after the start of pacing. Bright-field micrographs taken before and after 5 min of perfusion highlight cell drift in the traditional setup, which is one of the key limitations for capturing multiple functional records of cardiomyocytes over time using a traditional setup (Fig. 2 A). By contrast, cell loss or drift was not observed in the microwell device. Representative Ca2+ transients and sarcomere shortening records show the overall similarity between measurements made on the respective devices (Fig. 2 B). No significant dif- ferences between the mean Ca2+ release magnitude, Ca2+ decay constant, peak sarcomere shortening, cell length, or cell width were detected (Fig. 2 C). Both TTP (p < 0.001) and RT50 (p < 0.05) were reduced in the cells measured in PDMS microwells. Resting sarcomere length was shown to be slightly longer in cells tested in the new device (p < 0.05). Although statistically significant, the observed mean differences between the two methods were relatively small (~9%).
If anything, the larger mean resting sarcomere length of cells in the PDMS microwells may indicate bene- ficial effects of shielding cells from direct fluid shear, result- ing in slightly better cardiomyocyte viability. Fluid shear could cause incidental membrane damage, leading to Ca2+ leaking into the cell, slightly higher resting Ca2+, and slightly shorter resting sarcomere length. This could also result in slower relaxation, which is also seen in cells from the traditional setup. We next sought to establish the baseline variability of car- diomyocyte contractile behavior during experiments lasting many tens of minutes. This inherent variability over time is the backdrop against which drug responses occur. The behavior of cells was followed over the course of 1 h with repeated measurements made roughly every 2 min. To investigate the influence of the Fura-2 Ca2+ indicator on long-term cell dynamics, we tracked 12 cells loaded with in- dicator and 12 cells without. The cells were stimulated at 0.5 Hz to limit the effect of repeated pacing events over a long period of superfusion. Ca2+ transients and sarcomere shortening records of a Fura-2 loaded cell over five time points after preconditioning show no obvious rundown or degradation of cell integrity (Fig. 3 A).
Properties such as peak Ca2+ release and peak sarcomere shortening were calculated from measurements at each time point. Comparing these values to those that manifested in the same cell at the initiation of the experiment allows the fold change (FC) of the various functional metrics to be rep- resented over time (only eight representative cells from each group are shown for clarity; Fig. 3 B). Each cell’s FC values over time for a given parameter were fit with a Hill regres- sion to summarize time-based trends in the data as the cells reached a dynamic steady state (Fig. 3 C). The deviation of FC values from the Hill equation (the residuals) of each cell indicate that both Fura-2-loaded and dye-free cells fluctuate continuously over the entire course of the experiment (Fig. 3 D). These fits had an average root mean-square error (RMSE) of 0.062 FC for peak Ca2+ transient amplitude and 0.073 FC for peak sarcomere shortening in the FIGURE 2 Comparison of traditional and PDMS microwell methods. (A) Shown are bright-field micrographs of cardiomyocytes at 30× magnification in traditional bath (left) and PDMS microwells (right) before and after a 5 min interval to demonstrate cell drifting. Colored arrows indicate the same cells at different time points. (B) Shown are representative normalized Ca2+ transients and sarcomere shortening using traditional bath (left) and PDMS microwells (right). (C) Comparison of cell properties between traditional and PDMS microwell methods is shown. Statistical significance of means indicated by *p < 0.05 and **p < 0.001. Error bars indicate 95% confidence interval. To see this figure in color, go online.
Fura-2-loaded cells, with an average RMSE of 0.0528 FC for peak sarcomere shortening in dye-free cells. Fura-2 loading did not appear to affect the variation of peak sarco- mere shortening, diastolic sarcomere length, TTP, or RT50 over time based on a comparison of the RMSE means from all 12 cells in each group (see Fig. S1, A and B). Records of TTP, RT50, diastolic sarcomere length, and Ca2+ decay constant FCs from the representative cells are included in the Fig. S1 C. Having established the viability of cells in the microwell device and characterized their stability during prolonged ex- periments, we next used the device to quantify the responses of several individual cells to dobutamine and OM. We found that ensuring steady-state drug responses across all cells required more than 40 min of exposure. This behavior is evident in responses of four dobutamine-treated cells that were measured with high temporal resolution (once each minute; see Fig. S2). Cardiomyocytes were superfused with either drug (40 nM dobutamine, n = 41 cells; 150 nM OM, n = 34 cells) for up to 1 h after 12 min of pre- conditioning. For control groups, cells were also superfused with Tyrode’s solution alone (Control, n = 19 cells) as well as a vehicle-only control group superfused with 0.1%
DMSO in Tyrode’s (DMSO, n = 17 cells). Cells used in this experiment came from two different rat hearts, repre- sented in roughly equal proportion within each drug treat- ment group. Analysis of the data using a linear mixed- effect model confirmed that animal-to-animal variation was not significant. Representative Ca2+ transients and sarcomere shortening records illustrate the mean behavior of cardiomyocytes under the various test conditions (Fig. 4 A). Gray insets within the sarcomere shortening re- cords show the unadjusted transients to emphasize differ- ences in diastolic sarcomere length. Eight representative time-course records of the changes in Ca2+ and contractile properties during drug treatment are shown in Fig. 4 B. These responses show the positive inotropic effects of both dobutamine and OM. They also suggest substantial het- erogeneity in the responses of individual cells. Representa- tive time-course records of the Ca2+ decay constant and RT50 can be found in Fig. S3. As before, the time courses of cellular properties were fitted with a Hill regression to describe their dynamic steady-state responses (Fig. 4 B, gray lines). FC residuals from Hill fits were similar for the most part across all treat- ment groups, indicating that time variation in cellular behavior was not affected by the drugs studied (comparisons of RMSE between groups can be found in Table S1). From these fits, steady-state FC values for Ca2+ and contractile properties were aggregated for each treatment and compared via a one-way analysis of variance (Fig. 5). On average, dobutamine caused increases in Ca2+ release
(p < 0.0001) and peak shortening (p < 0.05) with little change in kinetics. OM also increased the magnitude of sarcomere shortening (p < 0.01) but with a concomitant lengthening of TTP (p < 0.0001) and RT50 (p < 0.0001) and a large decrease in the diastolic sarcomere length (p < 0.0001). The variance of steady-state response proper- ties was also assessed and compared between dobutamine- and OM-treated groups. Peak Ca2+ release (p < 0.0001), Ca2+ decay constant (p < 0.05), and peak sarcomere short- ening (p < 0.05) responses were all significantly more var- iable in dobutamine-treated cardiomyocytes compared to OM (Levene’s Test).
Dobutamine-treated cardiomyocytes also exhibited a sur- prising dissociation between increases in peak Ca2+ release and increases in peak sarcomere shortening. This is clearly seen in the records of two dobutamine-treated cells (Fig. 6 A). In the first case, the drug causes a robust increase in Ca2+ release, whereas sarcomere shortening actually de- creases. In another cell, the Ca2+ transient is unchanged even as sarcomere shortening doubles. These stark contrasts in response are still clearer when visualized as changes in each cell’s Ca2+ sarcomere shortening hysteresis loop (Fig. 6 B). Examining FC in steady-state Ca2+ release and peak sarcomere shortening for each of the dobutamine- treated cells shows that some form of mismatch is present in the majority of cells (Fig. 6 C). The same magnitude of mismatch is not seen in the control, DMSO, or OM groups shown in Fig. S4. Furthermore, 16 of 41 cells showed either no change or a decrease in contractility after dobutamine treatment.
FIGURE 4 Effects of dobutamine and OM on functional properties. (A) Shown are representative Ca2+ transients and sarcomere shortening records from each treatment group after preconditioning (gray) and after drug exposure (colored). Both are scaled to maintain the same peak sarcomere shortening values between the pretreatment transients. Sarcomere shortening records are aligned at the same diastolic sarcomere length. Gray subplots show both shortening records unaligned. Scale bar on the y axis indicates normalized baseline record magnitude. (B) Shown are eight representative time courses of functional parameters for control (n = 19), 0.1% DMSO control (n = 17), dobutamine (n = 41), and OM (n = 34) groups. Hill-fit regressions are shown in gray behind the time course data. Columns are labeled as in (A). To see this figure in color, go online.
This finding reinforced our original hypothesis concern- ing the variability of cardiomyocyte responses to b-adren- ergic stimulus. It also underscored the likelihood that the abundance of gene transcripts associated with this pathway would differ significantly among cardiomyocytes. Accord- ingly, we resolved to examine the transcriptomic activity of individual cells and attempt to correlate this information with their particular responses to dobutamine. A new batch of cardiomyocytes was superfused with 40 nM dobutamine for 1 h while monitoring cellular re- sponses (n = 31). The magnitudes of Ca2+ and shortening responses were used to rank order cells as before (Fig. 7 A). After functional measurements were taken, each cell was individually aspirated and deposited into its own microcentrifuge tube containing cell lysis buffer. Cell lysates underwent reverse transcriptase followed by a preamplification step to amplify our genes of interest: Ca2+/calmodulin-dependent protein kinase II (CaMKII) (Camk2d), protein phosphatase 1, protein phosphatase 2 (Ppp2ca), protein kinase A (PKA) gene (Prkaca), PKA reg- ulatory subunit type I (Prkar1a), PKA regulatory subunit type II (Prkar2a), and Gapdh (housekeeping). PKA phos- phorylates regulatory Ca2+ and sarcomeric proteins as a result of b-adrenergic stimulation. Protein phosphatase 1 and 2 dephosphorylate those same proteins. CaMKII phos- phorylates L-type Ca2+ channels as well as many sarco- meric proteins. Real-time PCR was used to quantify the abundance of each transcript relative to Gapdh and then compared to the transcript abundance in 12 DMSO control cells. Hierarchical clustering of gene expression FC with functional properties was used to generate a heatmap of —log2 (FC) values (Fig. 7 B). We subsequently probed six
pairs of cell properties that were related by clustering (Ppp2ca/RT50, peak sarcomere shortening/Prkaca, Peak Ca2+/Prkaca, protein phosphate 1/Camk2d, Camk2d/ Prkar1a, and Prkar1a/Prkar2a) to see if they were signifi- cantly correlated. This analysis yielded a significant positive correlation between the gene expression FC of Prkar1a and Prkar2a (p < 0.0001) as well as Prkar1a and Camk2d (p < 0.0001; Fig. 7 C). A weak correlation was also found between FC of Ppp2ca gene expression and RT50 (p = 0.066; Fig. 7 D).
DISCUSSION
This study comes at a time of increased interest in character- izing the functional behavior and gene expression of indi- vidual cardiomyocytes. For instance, micropatterned surfaces have been employed to generate large arrays of sin- gle induced pluripotent stem cell-derived cardiomyocytes (5,6). These arrays allow a variety of functional measure- ments to be performed on many individual cells. Others have succeeded in relating functional measurements of indi- vidual induced pluripotent stem cell-derived cardiomyo- cytes with cell-specific myosin heavy chain isoform expression (7). The microwell capture device described here is a technical advance that provides new and useful ways to link isolated adult cardiomyocyte function with gene expression. We have also demonstrated the ability of the device to quantify cell-to-cell variation in sustained drug responses. Our initial objective was to test the hypothesis that an inotrope operating via multiple intracellular intermediaries (dobutamine) would produce more variable responses than one targeted to a specific molecule (OM). This was clearly supported by the data (Fig. 5) and suggests the presence of substantial noise in the expression of genes that partici- pate in the b-adrenergic signaling pathway. Cardiomyocyte responses to OM were consistent with other recent reports FIGURE 6 Relative changes in peak Ca2+ and peak sarcomere shortening in dobutamine-treated cardiomyocytes. (A) Shown are representative Ca2+ transients and sarcomere shortening records of a car- diomyocyte that displayed a positive change in Ca2+ release with a negative change in peak sarcomere shortening (left) and a cell that displayed no change in Ca2+ release yet a positive change in peak sarco- mere shortening (right). Gray transients indicate re- cords taken immediately after preconditioning, whereas colored transients were taken after sustained dobutamine treatment. Both are scaled to maintain the same peak sarcomere shortening values between the pretreatment transients. Scale bar on the y axis in- dicates normalized baseline record magnitude. (B)
Ca2+ sarcomere shortening hysteresis loops of each FC in phenotype parameter of the select transients from panel A emphasize the uniaxial effects of dobutamine on each cell. (C) Each data pair correspond to a single cell’s FC response for peak Ca2+ and peak sarcomere short- ening, rank ordered by the difference in FC between the two functional properties. Cells from (A and B) are highlighted with their respective colors. To see this figure in color, go online(3,8), including the observation that OM increases contrac- tility at the expense of diastolic function. Perhaps the most striking finding in this study is the frequent occurrence of discordant Ca2+ contraction coupling responses in cardiomyocytes after dobutamine administra- tion (Fig. 6). Because intracellular Ca2+ is the primary signal regulating cardiomyocyte contraction, these two properties are expected to be tightly correlated. However, we observed several examples of dobutamine-treated cardiomyocytes that experienced decreased contractility despite elevated Ca2+ release, and others that showed large increases in contrac- tility with no corresponding increase in Ca2+. These seemingly paradoxical responses are less surprising when the complexity of the b-adrenergic signaling pathway and the diversity of its subcellular targets are considered. At the level of the intact myocardium, b-adrenergic stimulus causes both increased contractility (inotropy) and enhanced relaxation rate (lusitropy) (9). This highly desirable combi- nation of increased strength and speed is accomplished through a multitude of simultaneous subcellular changes. On one hand, Ca2+ release is augmented via PKA phosphor- ylation of L-type Ca2+ channels and other targets (10). PKA also phosphorylates myosin-binding protein C, which pri- marily increases the rate of cross-bridge binding (11). All of these events tend toward increases in the magnitude and duration of contraction. On the other hand, PKA phosphory- lates phospholamban and TnI. These events promote relaxa- tion and tend to reduce force production (12–14). When enacted together, these apparently antagonistic effects combine to produce increases in contractility without pro- longing the duration of the twitch (15).
The presence of antagonistic subcellular processes that are all activated by the b-adrenergic signaling cascade raises the possibility of extreme cell behavior if the opposing effects are not perfectly balanced. We believe that this explains the di- versity of dobutamine responses revealed in our experiments (Fig. 6). For instance, if a cell were randomly deficient in reg- ulatory proteins that anchored PKA specifically to myosin- binding protein C, such as myomegalin (16), it might display the expected increase in Ca2+ release but fail to show increased contractility. PKA could still be localized to the myofilament by other anchoring proteins such as troponin T (17), but the Ca2+-desensitizing effects of TnI phosphory- lation (13) would operate with less opposition. Interestingly, such scenarios were contemplated by Negroni et al., who constructed an integrative cardiomyocyte model represent- ing the various effects of b-adrenergic signaling on electrical, Ca2+ handing, and contractile systems of the cell (14). To illustrate the importance of each phosphorylation target on achieving the combined inotropic and lusitropic effects, they simulated cardiomyocyte Ca2+ transients and unloaded shortening while eliminating each target in turn. Their simu- lations show dramatic shifts away from the expected behavior when one or more components of the response are suppressed. In fact, the diversity of functional outputs they were able to predict is strikingly similar to the modes of discordant Ca2+ contraction responses we observed among measured single-cell responses.
Despite a high degree of individual variation, on average, the cardiomyocytes exposed to dobutamine behaved in the expected manner, namely an increase in cell shortening without prolonging the twitch duration (Fig. 5). It is possible In experiments that link aspects of the dobutamine response to gene expression (Fig. 7), we established the technical feasibility of this approach and illustrated the types of insights that can be gained. Hierarchical clustering of single-cell data indicated an association between RT50 (describing the lusitropic effects of dobutamine) and the expression of Ppp2ca (Fig. 7 D). This protein is responsible for reversing PKA-mediated phosphorylation of TnI. Its ac- tions would therefore tend to increase Ca2+ sensitivity of the myofilaments and slow relaxation. These data suggest that excessive Ppp2ca expression can abolish the lusitropic effects of b-adrenergic pathway activity. The most striking feature of the single-cell gene expres- sion data is the co-regulation of PKA regulatory subunits (Prkar1a and Prkar2a) and the gene encoding CaMKII (Camk2d) (Fig. 7 C). In one of the few studies examining sustained b-adrenergic pathway activation in cardiomyo- cytes, Wang et al. (18) demonstrated that levels of cAMP quickly surge and then gradually taper off. Meanwhile, CaMKII activity is slow to develop but then remains active indefinitely. For this reason, we had hypothesized that the expression of Prkaca (representing the acute response) and Camk2d (representing the chronic response) might be negatively correlated. In fact, no correlation was observed between the two genes (data not shown). The observation that expression of Camk2d and genes encoding PKA-associ- ated regulatory subunits were strongly correlated is there- fore surprising and underscores the need for additional experimentation. Expanded studies of larger populations of cardiomyocytes using an unbiased method for assessing gene expression (RNA sequencing) will allow for a more detailed accounting of the factors that are impacting vari- able cellular responses to dobutamine.
It must be acknowledged that phenomena identified in FIGURE 7 Gene expression analysis of dobutamine-treated cells corre- lated with FC in functional properties. (A) Each data pair correspond to a single cell’s FC response for peak Ca2+ and peak sarcomere shortening, rank ordered by the difference in FC between the two functional properties (n = 31). (B) Shown is an hierarchically clustered heatmap of functional properties and gene expression —log2 (FC) values. (C) Prkar1a gene expression FC correlated with both Prkar2a and Camk2d FC. Significant correlation represented by ****p < 0.0001. (D) Ppp2ca gene expression FC correlated with RT50 FC. To see this figure in color, go online that electrical and mechanical coupling between cardiomyo- cytes in the intact myocardium tempers cell-to-cell differ- ences. Nevertheless, the existence of this heterogeneity is meaningful in at least two ways. First, it provides eviden- tiary support for the premise advanced by Negroni et al., namely that each of the disparate subcellular functions tar- geted by b-adrenergic signaling is critical to its physiolog- ical role of simultaneously enhancing ventricular inotropy and lusitropy. Second, when coupled with single-cell biochemical assays, cell-to-cell differences can be leveraged to discover new features of the b-adrenergic pathway that could be exploited for cardiac drug development isolated rat cardiomyocytes are not guaranteed to apply unambiguously to physiological situations. The process of isolation disconnects cardiomyocytes from the extracellular matrix and from neighboring cells, with which they share mechanical and electrical connections (19). Undoubtedly, this perturbs cellular properties and may interfere Omecamtiv mecarbil with signaling microdomains critical to the b-adrenergic response (20,21). For these reasons, it will ultimately be necessary to confirm findings obtained with this system in multicellular preparations.
ACKNOWLEDGMENTS
The authors thank Adriel S. Sumathipala and Burak Dura for contributions toward early prototypes of the PDMS microwells.
This work was supported by a National Science Foundation grant (CMMI- 1562587) to S.G.C.