Abstract
Cardiac microtissues provide a promising platform for disease modeling and developmental studies, which require the close monitoring of the multimodal excitation-contraction dynamics. However, no existing assessing tool can track these multimodal dynamics across the live tissue. We develop a tissue-like mesh bioelectronic system to track these multimodal dynamics. The mesh system has tissue-level softness and cell-level dimensions to enable stable embedment in the tissue. It is integrated with an array of graphene sensors, which uniquely converges both bioelectrical and biomechanical sensing functionalities in one device. The system achieves stable tracking of the excitation-contraction dynamics across the tissue and throughout the developmental process, offering comprehensive assessments for tissue maturation, drug effects, and disease modeling. It holds the promise to provide more accurate quantification of the functional, developmental, and pathophysiological states in cardiac tissues, creating an instrumental tool for improving tissue engineering and studies.
Similar content being viewed by others
Introduction
Cardiac diseases are the leading cause of human morbidity and mortality1. Cardiotoxicity is also the main side effect preventing one-in-three drugs from being used clinically2. The lack of a close animal model makes in vitro cardiac tissues the important alternative3. Cardiac tissue built from human stem cell-derived cardiomyocytes (SC-CMs) is most frequently used3,4. They are also capable of retaining patient information for studying genetic diseases and personalized drug screening3,5. Clinical transplantation of these tissues may offer long-term hope for heart disease treatments6. Compared to planar tissue culture, three-dimensional (3D) cardiac microtissues (CMTs) are preferred tissue models because they can recapitulate cell phenotypes, microenvironment, and cell-cell interactions closer to the living organ7.
The assessment of developmental state, drug effect, and disease mechanism often requires the close monitoring of physiological responses3,4,5. The electrical and mechanical activities are the most relevant parameters8. Importantly, as these activities are inherently connected through the excitation-contraction (EC) coupling9, the simultaneous measurements of the correlated dynamics are highly desirable8. For example, a single-modal measurement of the electrical or mechanical response only is inherently inadequate to differentiate drug effects that work on the opposite part or both. In other cases, the dysfunction of myocytes in chronic myocardial infarction is linked to an impaired EC coupling10, and many arrhythmias are caused by a weakened EC coupling8. Thus, spatial mapping by tracking the electrical or mechanical response only is limited in identifying the weakened EC link. For developmental studies, tissue maturity or aging is often better characterized with correlated EC dynamics (e.g., evolutional details in the contractile phases)11,12 than the typical information of electrical or mechanical amplitude only.
Compared to planar tissue culture, CMTs introduce additional challenges to the effective assessment. Specifically, traditional optical means to track electrical activity (e.g., by voltage-sensitive dye)13 and mechanical motion (e.g., by micro post array)14 are now inherently limited to the surface region due to the opacity in the tissue. Electronic sensors integrated on a planar or flexible substrate, which work well with planar tissue culture15,16,17,18,19,20, now fall into the same scenario. The recent development of flexible mesh electronics provides the opportunity to embed addressable sensors in deep tissue21,22,23,24,25,26,27,28,29,30. In particular, the confined ribbon feature brings the mechanical properties in the mesh close to the tissue, which reduces the invasiveness and prolongs the interfacing stability28,29. However, existing platforms can only probe the single modality of either electrical or mechanical activity only21,22,23,24,25,26,27,28,29,30,31. This may be adequate for some (e.g., neural) tissues28 but falls short of characterizing cardiac systems, in which the correlated mechanical and electrical responses are of paramount importance as discussed before8.
The limitation is not without cause. Almost all existing electronic sensors are unifunctional and can probe the single property of either electrical or mechanical response only21,22,23,24,25,26,27,28,29,30,31. The conventional way of combining two complementary sensor types, which may work with the planar substrate, is now inherently limited for a confined ribbon-substrate integration. Among the limited number of demonstrations of simultaneous EC recordings in planar tissue culture, none of them can be translated to mesh integration, because at least one of the sensor technologies is either inaccessible32 to deep tissue or too big in device size33.
Fundamentally, the strategy of combining two sensor types inevitably leads to heterogeneity in integration and/or signal acquisition, which increases the overall footprint in interconnects and devices. This may be accommodated by a surface interfacing but is particularly unfavored by the deep-tissue interfacing that requires minimized substrate size. Moreover, spatial separation between the two sensors inevitably introduces inaccuracy in signal correlation, which is unfavored as the EC coupling originates from a very localized cell level.
Recently, we developed a 2-in-1 sensor concept of converging electrical and mechanical sensing in one device, by exploiting the field effect and piezoresistive effect in a bottom-up semiconducting silicon nanowire34. The nanowire was specifically assembled into 3D suspended structure on a planar substrate to attain geometric freedom for enabling multifunctionality34,35. Both electrical action potential and mechanical contraction were simultaneously detected in planar cardiac tissue culture34. However, the device integration in CMT has not yet been achieved, given that attaining scalable 3D nanowire assembly from the bottom up on a freestanding substrate can be nontrivial. Overall, no existing sensing platform has realized the tracking of the multimodal EC dynamics in live CMTs.
Results
Basic device design and characterizations
We developed a sensing platform to overcome the above limitations. The overall theme was to still exploit converging the multifunctional sensing in one nanoelectronic sensor, integrate the sensor array in an ultra-flexible mesh electronics system, and innervate the system in the CMT (Fig. 1a). To facilitate scalable integration, the sensor was made from a graphene transistor (Fig. 1b–i) designed to uniquely achieve the converged sensing of the multimodal EC dynamics. Specifically, on the one hand, the electrical cellular activity (e.g., action potential) can be detected by the transistor through the field effect36 (Fig. 1b–ii). On the other hand, the strain induced by the mechanical cellular motion can be transmitted by the ultra-flexible ribbon substrate and electrically detected by the transistor through the piezoresistive effect37 (Fig. 1b–iii). The ultra-flexible ribbon substrate is also important for avoiding the transfer of global bending to enable strict local motion detection. Together, one device can detect both electrical and mechanical activities simultaneously. The spatiotemporal convergence of both signals in the same device can enable the precise tracking of their correlated dynamics. The convergence can also lead to an equivalent reduction in overall device size and addressing budget, which is desirable for mesh integration.
The mesh electronics system was composed of a serpentine ribbon network to enable stretchability and accommodate tissue contraction. The ribbon was made of multiple layers (Fig. 1c; Figs. S1, S2). The bottom SU-8 layer (~400 nm thick) was lithographically defined for support and insulation. Monolayer graphene synthesized by chemical vapor deposition38 (Fig. S3) was transferred onto the SU-8 layer, patterned into array, and electrically addressed by an interconnect layer (e.g., Pd/Au/Pd, 15/40/5 nm) to form the sensor array. The interconnects were passivated by a top dielectric layer (e.g., 70 nm Si3N4) for insulation34. All the materials used were studied to be biocompatible and stable in physiological environment28,29,39. The mesh was fabricated on a sacrificial layer (e.g., 120 nm Ge) to be dissolved for release, except for the input/output (I/O) contacts for connecting to the recording system.
The fabricated mesh covered an area of 5.8 × 3.8 mm2 (Fig. 1d), with a ribbon width ~20 µm (e.g., smaller than cell) and a spacing between neighboring ribbons ~200 µm (e.g., larger than cell) (Fig. 1e). The total thickness was ~530 nm (Fig. S4). These features yielded a highly porous system with a softness comparable to tissue23, which is important for establishing intimate tissue interfacing and effective signal transduction. The size of the transistor (e.g., 20 × 20 µm2) was comparable to cell size to enable localized detection of cellular action potential and mechanical contraction (Fig. 1f). Raman spectrum confirmed the presence of the graphene transistor (Fig. 1g), showing the typical38 G peak at 1594 cm−1 and a 2D peak at 2650 cm−1 against the peak-less SU-8 substrate (gray curve).
Electrical characterizations were performed in the integrated transistors to reveal the potential for recording electrical and mechanical cellular responses. First, the electrical sensitivity was obtained by measuring the water-gate response in cell culture medium. The devices exhibited typical ambipolar behavior with the maximal transconductance ~2.2 ± 0.4 mS/V (Fig. 1h; Fig. S5), consistent with typical values obtained on rigid substrates36. This enhanced transconductance (e.g., compared to values in silicon nanowire transistors34) can improve detection. Note that the sensor detection resolution is still affected by the noise level in individual devices. Our characterizations revealed the typical noise level ~20–40 nS (Fig. S5), equivalent to an estimated detection resolution ~30–60 μV (e.g., with a signal-to-noise ratio >3) that was sufficient for detecting extracellular action potentials15. The transconductance maintained stable value during repeated cycles of mesh folding (e.g., up to 180o) and stretching (e.g., 30% biaxial strain covering cardiac contractile range) (Fig. S6). Second, the mechanical sensitivity was obtained by measuring the piezoresistive response, which showed a linear decrease in conductance with the increase of bending stain (Figs. 1i, S7). The gauge factor (i.e., the ratio of fractional resistance change to the fractional strain change) was calculated to be ~82 ± 24, consistent with the typical range observed previously37. With the relative noise level ~10−5 (Fig. S5), the obtained gauge factor suggested that devices could detect a strain level <10−6, sufficient for capturing cellular motion that typically induces local strain >10−5 34. These properties suggested that the integrated graphene transistors were capable of detecting cellular electrical and mechanical activities.
Tissue integration and recording
The integrated mesh system was released in a chamber and fixed by a Matrigel layer, with the I/O contacts connected to the recording system (Figs. 2a, S8). Human embryonic stem cell-derived cardiomyocytes (hESC-CMs) were seeded on the mesh. During the continuous development, the tissue sheet grew, folded, and engulfed the mesh by the increasing cell-cell interactions to form a mesh-innervated CMT (Figs. 2b, S9). The CMT continued to grow and maintained the integrity and contractile function across the entire developmental time of one month (Fig. 2c; Supplementary movie 1), showing the initial evidence of a stable mesh-tissue interface.
Multi-channel recordings revealed synchronized and periodic signals charactering the beating (~0.6 Hz) in the CMT (Fig. 2d). The recording yield (e.g., 9/14 ~ 64%) was reasonable by taking consideration of yields in device fabrication (e.g., ~90–95%), functional cardiomyocytes (e.g., 80–98%) during differentiation40, and potential inner necrotic part in the CMT resulted from supply limit of nutrient and oxygen by diffusion41. In each period, the broad peak was preceded by a sharp spike (dashed box, Fig. 2e). These sharp spikes showed a biphasic feature with the duration ~27 ± 6 ms (Fig. 2f) and calibrated amplitude of ~70–200 µV (Fig. S10), consistent with the properties of the extracellular action potential in hESC-CMs24,27. The noticeable variation in the shape features of the recorded action potentials can be attributed from the mixed cardiomyocyte subtypes differentiated from hESCs and cell states affected by diffusion gradient in oxygen and nutrient41,42. The time discrepancy between these action-potential signals (Fig. 2f) showed that the sensor array had a high temporal resolution to capture the propagating dynamics across different locations.
Additional investigations were carried out to correspond the broad peak to the mechanical response in the CMT. First, graphene transistors integrated on an unreleased mesh on the rigid substrate only detected the action-potential spikes but not the broad peak (Fig. S11), suggesting that the mechanical contraction was the source of the broad peak. Second, mesh containing only interconnects (e.g., without transistors) innervated in the CMT did not detect any signal (Fig. S12), confirming that the broad peaks were not mechanical artifacts from interconnects but active piezoresistive signals from the transistors. Third, shifting the gate reference (Vg) to move the graphene transistor from a p-type transport region to a n-type one reversed the signs in the broad peaks (Fig. S13). This is consistent with the general observation that the piezoresistive coefficient reverses the sign between n- and p-type conduction43, which was further confirmed by the Vg-dependent strain response in the graphene transistors (Fig. S14). Finally, the broad peaks showed close evolution to that of the Ca2+ dynamics governing the contractile dynamics in cardiomyocytes (Fig. S15)44. These results confirmed that the integrated graphene transistors also captured the detailed mechanical responses in the CMT.
Tracking maturation process
Improving maturation in CMT is key to recapitulating close organ function for precision studies27,45. An effective assessing means is indispensable to tissue engineering and studies for accelerating maturation8. Mature cardiomyocytes are generally characterized by concerted electrical and mechanical activities (e.g., fast depolarization, strong contraction, and close EC coupling)7. Previous unifunctional sensors only tracked one of these physiological components21,22,23,24,25,26,27,28,29,31, falling short of a comprehensive assessment. The multifunctional graphene sensor could track all these components and provide an enriched set of parameters (Fig. 3a). Therefore, we intended to investigate if the integrated mesh system can provide a much more detailed tracking of the CMT development for the comprehensive assessment of maturation.
Continuous recordings across the entire maturing stage (e.g., > 1 month)27,46 in the CMT were carried out. The transistor sensors maintained stable performance throughout the recording period (Fig. S16), forming the basis for capturing state evolution in the tissue. The key parameters charactering the electrical, mechanical activities, and their correlation (Fig. 3a) were extracted for the analysis of the CMT state. The recordings started from the 16th day of cell differentiation (corresponding to the 5th day of cell seeding on the mesh, Fig. S17). Representative recording showed a gradual increase in both the beating frequency f (e.g., from 0.5 to 1.4 Hz) and contractile amplitude AC (e.g., ~500% increase) during the process (Fig. 3b). An overlay of the mechanical signals also revealed a noticeable decrease in the contractile duration t2, mainly contributed from the relaxation phase t5 (Fig. 3c). A zoom-in overlay of the action potentials revealed an increase in the amplitude AE but a decrease in the duration t1 (Fig. 3d)
Statistical data from other devices featured the consistent trends of increasing beating frequency (f), contractile amplitude (AC), action-potential amplitude (AE), but decreasing action-potential duration (t1) and contractile duration (t2) mainly contributed from the decreasing relaxation phase (t5) (Fig. 3e). These trends were consistent with the correlated cellular mechanisms in a maturing cardiomyocyte. Specifically, maturing cells are featured with an increasing rate of activation in the Na+ channels (e.g., upstroke velocity)47, which was reflected from the shortening of the action-potential duration t1. Maturing cells also recruit more Na+ channels45 to yield an enhancing action-potential amplitude AE as observed. The increasing contractile amplitude AC suggested the improving contractile force contributed from the strengthening myofibrils, which was confirmed by fluorescence imaging showing increasing intensity of cardiac troponin T (cTnT)48 and improving sarcomere alignment (Fig. 3f, g). The multi-channel recordings also captured the temporal signal delay between different sites (Fig. S18), revealing an increasing conduction velocity (v) consistent with the enhancing electrical signaling in maturing tissue.
The simultaneous measurements of the action-potential and mechanical signals provided a unique window to peek into the Ca2+-handling dynamics underlying the EC coupling. A faster upstroke velocity in the action potential leads to the faster triggering of the L-type Ca2+ channel for extracellular Ca2+ recruitment44, which further induces the Ca2+-induced Ca2+ release (CICR) in the sarcoplasmic reticulum (SR) for contraction. This improving EC coupling was uniquely captured by the integrated sensors, which showed a decreasing EC temporal delay (t3). The convergence of both electrical and mechanical sensing at the same location ensured the precision in the EC correlation, which was not possible with previous sensor technologies21,22,23,24,25,26,27,28,29,31. Meanwhile, the increasing rate of Ca2+ uptake by the SR leads to a shortening relaxation phase in maturing cells49, which was clearly reflected by the shortening t5 (also contributing to the overall decrease in the entire contractile duration t2). This unique capability of tracking high-resolution contractile details may be employed in studying cardiac aging, which can be subtly indicative of a prolonged relaxation in the contractile process without altering the force amplitude11,12. All these features were consistent with improving Ca2+ handling and EC coupling during maturation. While post-dissection imaging revealed morphological evolution at sporadic stages to confirm the maturing process (Fig. 3f, g), the innervated mesh system enabled real-time tracking of the live CMT across the entire development and provided enriched physiological parameters beyond previous platforms.
Drug tests
One of the main purposes of developing cardiac tissue models is to study drug effects because the heart is very susceptible to side effects from drugs2. Cardiotoxicity accounts for one-third adverse drug reactions and many are associated with ion-channel disturbance and result in arrhythmia50. The real-time tracking of the live CMT with enriched physiological parameters provides unique opportunity for the efficient detecting and differentiating drug effects. Drugs that work on different cell mechanisms were used to demonstrate the potential. Dimethyl sulfoxide (DMSO), the common carrier of drugs, showed negligible effect on tissue state (Fig. S19), establishing the baseline for following tests. First, the innervated CMT was treated with Blebbistatin (10 µM), an inhibitor of the ATPase activity associated with the actin-myosin complex, to suppress contraction51. Electrical recordings showed a prominent suppression of the mechanical signals after the treatment (Figs. 4a, S20). Analysis of the action-potential signals showed negligible change in the amplitude AE and duration t1 (Fig. 4b, c), suggesting a non-ionotropic mechanism consistent with the drug effect. Verapamil, another inhibitor of cell contraction, relies on the blockage of L-type Ca2+ channels52, which in turn also reduces the Na+-channel activity53. As a result, the suppression of the mechanical contraction was coupled with weakened action potentials, which were reliably captured by electrical recordings (Fig. S21). These results show that the enriched physiological parameters provided by the innervated mesh system can discern mechanistic differences in close drug effects. The mesh-tissue intimacy also enabled the tracking of the evolution in drug effects. In one case, the mesh system provided a detailed comparison of the tissue state before and after recovery (e.g., washing out Blebbistatin) (Fig. S22). In another case, the mesh system revealed the dose-dependent drug effect of Verapamil, showing negative inotropic and chronotropic responses31 with the increase of dosages (Fig. S23).
Conversely, we also constructed an electrically dysfunctional model by applying quinidine (1 μM), a Na+-channel blocker54, to the innervated CMT and recorded the signals (Fig. 4d). Examination of the overlayed signals revealed a decrease in action-potential amplitude AE but a prolongation in its duration t1 (Fig. 4e, f). The trend was consistent with the drug mechanism in reducing and slowing Na+ influx, which is the key component contributing to the recorded action-potential signal. The weakened Na+-channel activity also affected the efficiency in Ca2+ channels55, leading to a decreased amplitude AC (Fig. 4d) and a prolongation t2 in contraction (Fig. 4e, f). Dose-dependent studies enabled by the mesh further revealed that a dosage beyond 10 µM led to the full suppression of action potential and subsequent suppression of contraction (Fig. S24).
The above drug tests demonstrated that the innervated mesh system can track details to differentiate drug mechanisms that otherwise yield similar contractile impairment. The fine tracking may also help to capture subtle functional deviation. We therefore used the innervated sensing system to detect potential cardiotoxicity from doxorubicin, an effective anticancer agent commonly used for chemotherapy56. Although under normal dosage (e.g., <1 µM) the tissue maintained the basic electrical and mechanical functions (Fig. 4g, S25), a close examination of the recorded signals revealed changes in the EC dynamics (Fig. 4h). Specifically, although the decrease in the contractile amplitude AC was small (e.g., ~5%), a prominent prolongation in the contractile duration t2 (e.g., >20%) was observed. The prolongation indicated the slowing in Ca2+-handling dynamics, which was further suggested by the increased (~20%) EC delay t3. However, unlike previous cases where the contractile prolongation mainly resided in the relaxation phase t5 (Fig. 3c, 4e), here the contraction phase t4 contributed proportionally (Fig. 4h, i). Meanwhile, the action potential only showed a slight decrease in the amplitude AE, suggesting that an immediate ionotropic mechanism was not the direct cause. Although the mechanism is not fully known, the above observations are consistent with some existing proposed mechanisms, in which the doxorubicin produces a toxic metabolite that inhibits the sodium-calcium exchanger; this will further lead to a slowdown in the Ca2+ uptake by sarcoplasmic reticulum and also a decrease in the L-type calcium channel activity57, and hence, the prolongation in both relaxation and contraction phases. The innervated mesh not only readily captured the contractile deviation for identifying possible cardiotoxicity from doxorubicin, but also provided details for examining proposed molecular mechanisms.
Disease modeling
We continued to demonstrate the potential of the innervated CMT for disease modeling. Myocardial infarction, which is usually caused by the occlusion of blood flow or deficient oxygen supply in local regions to result in cell malfunction or death, is the most common cause of heart failure58. The border zone between the infarcted and distal healthy tissue can undergo significant remodeling in the electrical and mechanical properties59, which can eventually lead to heart failure. Tracking the EC-dynamics evolution in the border zone is crucial for the mechanistic understanding and eventual development of effective intervention and treatment60. While CMTs provide a convenient modeling platform, the typical use of an enclosed space to control the oxygen content creates a huddle for the real-time probing of the modeling result.
The innervated CMT is free of the huddle and can be readily placed in an enclosed microfluidic system capable of oxygen-flow control for the modeling study (Fig. 5a–c, S26). The CMT was subject to a hypoxic environment (10% oxygen) for 10 h before returning to a normoxia (21% oxygen) (Fig. S27) to emulate the conditions of ischemia and reperfusion, respectively61. The electrical recordings showed the development of tachycardia featuring a much faster (e.g., ~65% enhancement) beating rate after the hypoxic stress for 2 h (Fig. 5d–i, ii), which was attributed to the enhanced Na+ current during hypoxia to result in Ca2+ overload62. Within the next 4 h, the beating rate gradually reduced to value below the initial level (Fig. S28). Arrhythmia that features irregular beating was observed after 8-h hypoxic treatment (Fig. 5d–iii), which was a characteristic symptom of ischemia63. Reestablishing a normoxic environment led to the gradual recovery of the beating rate (at 12 h, Fig. 5d–iv) close to the initial level (at 0 h). However, the recovery in the electrical conduction velocity, suggested by the temporal delay between the spatially distributed sensors (right panels, Fig. 5d), was slower after normoxia (Fig. S29). This suggested a decrease in the gap-junction Connexin-43 expression and an increase in tissue impedance consistent with previous studies64. The mechanical response recorded by channel-3 (blue curve) was completely suppressed after 8-h hypoxic treatment and did not recover by reperfusion (at 12 h), suggesting the irreversible loss of cell activity in the local region. As the locations of the corresponding sensors could be determined by later-stage imaging (Fig. S30), region-specific physiological activities could be captured to identify the ischemia-reperfusion injury65. Single-modal measurement such as electrophysiological recording alone is incapable of capturing such spatially differentiated EC decoupling effect for detailed tracking.
Analysis in the other parameters revealed additional information. A consistent trend of decreasing amplitude (AE) and increasing duration (t1) in the action potentials was observed after 6 h of hypoxia (Fig. 5e–i, ii), which was consistent with results from in vivo models of isolated guinea-pig ventricular myocytes and rabbit atrioventricular node. The result was also consistent with the mechanism of hypoxia-induced inactivation of the fast Na+ channels66, which can result in a decrease in both the amplitude and speed of the rapid depolarization phase in the action potential.
The mechanical contraction showed a slight decrease in the amplitude (AC) but more prominent prolongation (t2) that was mainly contributed from the contractile phase (t5) after 6 h of hypoxia (Fig. 5e–iii–v; Fig. S31). Examination of the EC delay t3 revealed a similar prolongation (Fig. 5e-vi; Fig. S32), suggesting a slowed triggering of the CICR process. The prolongations did not fully recover after the reperfusion (at 12 h). These features were expected from slowed Ca2+-handling dynamics associated with the slowed upstroke velocity in the action potential (i.e., t1 prolongation)67. Together, the lack of a full recovery in both the electrical and mechanical responses after the re-establishment of normal oxygen levels suggested that the cardiac tissue suffered irreversible damage during the process. The innervated mesh provided a window to observe the temporal and regional susceptibility of the cardiac tissue in response to hypoxic conditions. Capturing the differential details in response was essential for comprehending the complexities of cardiac function and developing targeted therapeutic strategies for ischemia-reperfusion injury.
Discussion
In conclusion, our study has demonstrated a mesh sensing system with converged multifunctionality to track the detailed EC dynamics in cardiac tissues beyond previous platforms. The convergence of the multifunctional sensing in one device circumvents the challenges in signal synchronization, device addressing, and integration heterogeneity involved in the typical strategy of combining multiple sensor types. These budget reductions are very important for mesh integration and tissue interfacing, in which minimized invasiveness is always preferred. The mesh platform provides tissue-level softness and cell-level feature size to enable intimate interfacing with cardiac tissue. Together, the innervated sensing platform enables stable tracking of the enriched EC dynamics across the tissue and throughout the developmental process. These performance advances in the system enable comprehensive assessment for tracking cardiac tissue maturation, differentiating drug effects, and modeling disease as demonstrated in the study, providing enriched data sets far beyond what was possible with previous single-modal sensor technologies21,22,23,24,25,26,27,28,29,30,31. The sensing platform holds the promise to provide more accurate quantification of the functional, developmental, and pathophysiological states in cardiac tissues, creating an instrumental tool for improving cardiac tissue engineering and studies. Previous studies have shown that graphene can maintain functional stability68,69 and the mesh electronic scaffold can establish long-term interfacing intimacy28,70,71,72 in in vivo tissue environments. These properties suggest that the graphene-integrated mesh platform, combined with efficient delivery and addressing strategy73, can also be used in in vivo cardiac and neural systems for chronic recordings and studies. The platform also provides a unique tool to converge cell mechanics and electrophysiology studies—two very important but traditionally less overlapped fields.
Methods
Graphene synthesis
For monolayer-graphene synthesis, we used single-crystal Cu foil obtained by annealing polycrystalline Cu foil as the growth substrate. Graphene was then synthesized on the Cu foil in a low-pressure chemical vapor deposition (LPCVD) system using a modified two-step growth recipe38 (Fig. S3). In the first step, the temperature was ramped from 25 °C to 1050 °C within 30 min. Then, 40 sccm of methane (CH4) and 0.5 sccm of hydrogen (H2) were introduced to the LPCVD system for the graphene synthesis. This resulted in the growth of multilayer graphene on the front side of the Cu foil. The system was maintained at 1050 °C for 40 min during this step. In the second step, the ratio of CH4 to H2 was changed to 1.5:40 (sccm), and the temperature was maintained at 1050 °C for 120 min. After this, the temperature gradually cooled down to 25 °C in 40 min. Adlayer-free, single-crystal graphene was synthesized on the backside of the Cu foil during this second step.
After the synthesis of graphene on Cu foil, a layer of PMMA 950 A4 (Kayaku Adv. Mater., Inc.) was spin-coated on the backside of the Cu foil and baked for 5 min at 80 °C. The frontside graphene was etched by oxygen plasma (50 W, 2 min, 50 sccm O2). The sample was then floated on the Cu etchant (Transene, Inc.) for 30 min to remove the Cu foil. The released PMMA/graphene film was rinsed with deionized (DI) water and ready for the subsequent transfer process.
Fabrication of the integrated mesh system
Detailed procedure is illustrated in Figs. S1, S2. The key steps include: (1) a patterned sacrificial layer (120-nm Ge) was deposited on the substrate (Si/SiO2 or glass wafer). (2) A bottom SU-8 layer (2000.5; Kayaku Adv. Mater., Inc.) was spin-coated on the substrate and hard baked (For fluorescence imaging, 0.0008 wt% of Rhodamine 6 G power (Sigma-Aldrich) was mixed into the SU-8 solution). (3) The PMMA/Graphene film was then transferred onto the SU-8 layer, with the PMMA removed in acetone. (4) The graphene array was defined by lithographically patterning a protective resist array (PMGI SF6/S1805; Kayaku Adv. Mater., Inc.), oxygen-plasma etching (50 W, O2 50 sccm for 1 min) of unprotected graphene, and removal of the resist array. (5) A similar step was used to define a protective resist layer for defining the mesh pattern, etch the unprotected SU-8 film by reactive ion etching (RIE), and removal of the protective resist layer. (6) Standard photolithography and metal deposition were used to define the metal contacts and interconnects (Cr/Pd/Au/Pd, 3/15/40/5 nm). The contacts/interconnects were further passivated by a layer of Si3N4 (~70 nm). (7) The fabricated mesh system was released from the substrate by etching the Ge layer in 1% H2O2 solution for 30 min and rinsed in DI water before cell culture.
Raman spectrum
The Raman spectrum of the fabricated graphene devices was performed using a DXR Raman microscope (Thermo Fisher Scientific) with 600-nm excitation laser. The spectrum was recorded with an acquisition time of 40 s at a power of 4 mW.
Electrical measurements
The water-gate effect and mechanical response in the graphene transistors were conducted in DPBS (Gibco) by using a semiconductor analyzer. The gate voltage was applied to a gold wire immersed in the DPBS solution. In the water-gate characterization, the conductance in the devices was measured under a fixed DC bias (10 mV) with a gate voltage (sweeping from 0 to 0.3 V) applied to the gold-wire gate electrode. The piezoresistive responses in the devices were performed by bending the substrate. Details can be found in Fig. S7. The in vitro electrical recordings from the cardiac tissues were performed with a DC bias (5–10 mV) applied to the drain. The drain current was amplified with a home-built amplifier system, and the output was acquired at 20 kHz by an A/D converter (Digidata 1440 A, Molecular Devices) interfaced with computerized software (pClamp 10.7, Molecular Devices). Data analysis was carried out using OriginPro (version 2022, OriginLab).
Cell culture
hESCs (WAe009-A, H9, WiCell, NIH registration no. 0062, female. The cell line is commercially available and has been tested or authenticated by the manufacturers by standard techniques including morphology check, isoenzyme analysis, and mycoplasma detection as shown on their webpages: https://www.wicell.org/home/stem-cells/catalog-of-stem-cell-lines/wa09.cmsx) were used for the differentiation of cardiomyocytes by the following methods (Fig. S9)40. hESCs were maintained in the 60-mm tissue culture dishes coated with Matrigel (10 μg/ml, Corning) in DMEM-F12 (Gibco) using Essential 8 medium (Gibco) and subpassaged every 3 to 4 days. Cardiomyocyte differentiation from hESCs was conducted following the Wnt signaling modulation protocol24,40. During differentiation, hESCs were seeded in a 12-well plate for 2 to 3 days until 100% confluency. On day 0 of differentiation, the medium was changed to RPMI 1640 medium (Gibco) and 1% B27-insulin (Gibco) and 10 μM CHIR99021 (Tocris Bioscience). After 24 h (day 1), the medium was replaced to RPMI 1640 plus 1% B27-insulin. On day 3, the medium was changed to RPMI 1640 plus 1% B27-insulin, supplemented with 5 μM IWR-1-endo (Cayman Chemical). On days 5 and 7, the medium was changed to RPMI 1640 plus 1% B27-insulin, and RPMI 1640 plus 1% B27 (Gibco), respectively. The cells were maintained with RPMI 1640 plus 1% B27 every other day. The contraction of cells was usually observed on day 8. During days 10 to 12, the cardiomyocytes were ready for transfer and tissue culture.
For planar tissue culture (for device control test), the cells were rinsed with 1× DPBS to remove Ca2+ and inhibit contraction. Cells were dissociated by treating 0.5 ml of 0.5 mM trypsin-EDTA (Ethlyenediaminetetraacetic acid, Gibco) per well for 4–5 min in a 37 °C incubator. The trypsin-EDTA was then aspirated, and cells were dissociated by gently pipetting with 2 ml RPMI 1640 plus 1% B27 using a 1-ml pipette tip. The cells were transferred to a 15-ml conical tube and centrifuged at 250 g for 3 min and then resuspended with 2 ml RPMI 1640 plus 1% B27 with 10 μM Y-27632 ROCK inhibitor (Tocris Bioscience). The substrate having the unreleased mesh was sterilized with 70% ethanol solution (1 h) and ultraviolet (1 h) and coated with Matrigel (20 μg/ml) in RPMI 1640 for 1 h at 37 °C. Cells were seeded on the device substrate at the density of 3–5 × 105/cm2. The cells were maintained using RPMI 1640 plus 1% B27 with daily change. Electrical recordings were typically performed starting from day 7 after the cell seeding.
The innervation of the integrated mesh system into the CMT involved following steps: (1) Releasing the mesh system from substrate by etching the Ge sacrificial layer in 1% H2O2 solution for 30 min. (2) The released mesh was rinsed by DI water and sterilized by 70% ethanol for 1 h. (3) The mesh was washed with DPBS three times followed by the overnight treatment of Poly-D-lysine (0.1 mg/ml, Gibco) in incubator. (4) The Poly-D-lysine solution was aspirated, and the mesh was embedded in 100 μL Matrigel (10 mg/ml and incubated for 30 min to solidify), which also provided the initial fixation of the mesh on the substrate. (5) The resuspended cells in RPMI 1640 plus 1% B27 with 10 μM Y-27632 ROCK inhibitor were transferred onto the cured Matrigel with a density (3–4 × 106 cells/cm2) and incubated for 24 h. Then the cells were maintained using RPMI 1640 plus1% B27 by with the daily change of the medium. The electrical recordings typically started at day 5 of seeding (day 15–17 of the differentiation) and continually for more than 30 days.
Calcium imaging
The mesh-innervated CMT in the chamber was rinsed with 1 ml of Hepes-buffered Tyrode’s solution (Thermo Fisher Scientific) and stained with Fluo-8 AM (5 mM, AAT Bioquest). The CMT was then incubated for 1 h at 37 °C. Subsequently, the system was rinsed with Tyrode’s solution three times. The Ca2+ transient signal was recorded with a confocal microscope (Nikon A1R). Electrical recordings were performed simultaneously. The recorded signals were analyzed by ImageJ.
CMT cryosection, immunostaining and imaging
CMTs at different culture days were rinsed with DPBS three times and fixed with 4% paraformaldehyde (PFA)/DPBS at 4 °C overnight. The CMTs were then transferred to 30% sucrose/DPBS solution for 24–48 h. The samples were embedded in optimal cutting temperature compound (Fisher Scientific) and cryosectioned into 30 μm-thick slides. For immunostaining, the CMT slices were treated with washing solution three times (0.2% Triton, Thermo Fisher Scientific) for 15 min each time. Slices were then washed in blocking solution (10% normal goat serum, Thermo Fisher Scientific) for 1–2 h in a humidified chamber at room temperature and subsequently replaced with primary antibody (Cardiac Troponin T Monoclonal Antibody Ctnt (clone number 13–11), diluted in the blocking solution with a ratio of 1:250, Invitrogen, Catalog # MA5-12960, LOT 2504985) diluted in blocking solution at 4 °C overnight. The secondary antibodies (Alexa Fluor™ 647, Goat anti-Mouse IgG (H + L) Cross-Adsorbed Secondary Antibody, diluted in the blocking solution with a ratio of 1:500, Invitrogen, Catalog # A-21235, LOT 1987304.), 4’,6-diamidino-2-phenylindole (DAPI, diluted in the blocking solution with a ratio of 1:1000, Invitrogen), and Phalloidin (Alexa Fluor™ 488, diluted in the blocking solution with a ratio of 1:50, Invitrogen) were co-stained for 1–2 h at room temperature. The slices were washed with DPBS three times before imaging. All samples were imaged by Nikon A1R confocal microscopy. Laser wavelengths of 405 nm, 488 nm, 561 nm, and 640 nm were applied to excite components labeled with DAPI, Phalloidin, Rhodamine 6 G, and cTNT, respectively. Images and fluorescence intensities were analyzed by NIS-Elements version 5.3.
Whole-mount CMT immunostaining
The CMT was fixed with 4% PFA/DPBS solution at 4 °C overnight and washed with DPBS three times. The whole CMT was treated with 1 ml washing solution and 1 ml blocking solution with primary antibody (diluted in the blocking solution with a ratio of 1:250) at 4 °C overnight, respectively. The CMT was washed in blocking solution two times for 2 h and subsequently incubated with secondary antibody (diluted in blocking solution with a ratio of 1:500), DAPI (diluted in the blocking solution with a ratio of 1:1000), and Phalloidin (diluted in the blocking solution with a ratio of 1:50) overnight at 4 °C. Then the CMT was washed in DPBS two times for 0.5 h and mounted with Fluoromount-G™ Mounting Medium (Invitrogen). The imaging of the whole CMT was performed by Nikon A1R confocal microscopy with a step of 5 μm and the 3D view of fluorescence images of CMT was reconstructed by NIS-Elements version 5.3.
Drug assays
Drug assays were performed when the CMTs showed stable electrical and mechanical signals after culture for several days. Before drug testing, the medium was changed and the drug assays were conducted after 4 h of recovery. Drugs were dissolved in dimethyl sulfoxide (DMSO, Thermo Fisher Scientific) with a final concentration of less than 1%. Cumulative dosing of drugs was performed in the CMT until the recording signals presented stability. Drugs were purchased from Cayman Chemical.
Integration with microfluidic system
The assembly of the microfluidic chamber and environment control system are illustrated in Figs. S26, S27 in detail. Briefly, the microfluidic chamber template/mold was laser cut from a poly(methyl methacrylate) plate (radius: 4.5 mm, height: 7 mm, volume: ~445 μL). The microfluidic chamber was fabricated by casting polydimethylsiloxane (PDMS, base: cure agent=10:1, Sylgard 184), and the released PDMS chamber was attached to the substrate with a thin layer of PDMS glue. The template/mold for the microfluidic chamber cap with micropillar arrays of different radii (0.5, 0.4, and 0.3 mm) was also laser cut from a poly(methyl methacrylate) plate and PDMS was casted to prepare the cap. The micropillar arrays were designed to eliminate the bubbles in the chamber74. The cap and chamber were bonded with double-coated spacer tape (ARseal 94119, Adhesives Research). Fluid was delivered to the channel using polyethylene tubes (BTPE-50, Instech), with the speed programmably controlled by a syringe pump (Harvard Apparatus, Standard Infuse/Withdraw Pump 11 Pico Plus Elite).
Hypoxic medium was prepared by bubbling ultrapure nitrogen and oxygen (9:1, Airgas) saturated with water to the culture medium (RPMI-1640 with HEPES (Gibco), pulsing 1% B27). The oxygen concentration was calibrated by an oxygen sensor (PCE Instruments). After 10-h hypoxia, the cell medium was replaced by normoxic medium. During the medium switching process, the flow rate was increased from 250 μL/h to 250 μL/min for ~2 min to flush the channel with required medium.
Statistics and reproducibility
Microsoft Excel (version 2021) was used to analyze all data in this work. No statistical method was used to predetermine the sample size. No data were excluded from the analyses. The experiments were not randomized. The investigators were blinded to group allocation during data collection and analysis.
Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Data availability
The data supporting the findings of this study are available within the article and its supplementary files. Any additional requests for information can be directed to, and will be fulfilled by, the corresponding author. Source data are provided in this paper.
References
Nabel, E. G. Cardiovascular disease. N. Engl. J. Med. 349, 60–72 (2003).
Laverty, H. G. et al. How can we improve our understanding of cardiovascular safety liabilities to develop safer medicines? Br. J. Pharm. 163, 675–693 (2011).
Yang, C. B. et al. Concise review: cardiac disease modeling using induced pluripotent stem cells. Stem Cells 33, 2643–2651 (2015).
Mathur, A., Ma, Z., Loskill, P., Jeeawoody, S. & Healy, K. E. In vitro cardiac tissue models: current status and future prospects. Adv. Drug Deliv. Rev. 96, 203–213 (2016).
Ashley, E. A. The precision medicine initiative a new national effort. JAMA J. Am. Med Assoc. 313, 2119–2120 (2015).
Laflamme, M. A. & Murry, C. E. Heart regeneration. Nature 473, 326–335 (2011).
Arslan, U., Orlova, V. V. & Mummery, C. L. Perspectives for future use of cardiac microtissues from human pluripotent stem cells. ACS Biomater. Sci. Eng. 8, 4605–4609 (2022).
van Meer, B. J., Tertoolen, L. G. J. & Mummery, C. L. Concise review: measuring physiological responses of human pluripotent stem cell derived cardiomyocytes to drugs and disease. Stem Cells 34, 2008–2015 (2016).
Bers, D. M. Cardiac excitation-contraction coupling. Nature 415, 198–205 (2002).
Kim, Y. K. et al. Altered excitation-contraction coupling in myocytes from remodeled myocardium after chronic myocardial infarction. J. Mol. Cell Cardiol. 34, 63–73 (2002).
Fares, E. & Howlett, S. E. Effect of age on cardiac excitation-contraction coupling. Clin. Exp. Pharm. P 37, 1–7 (2010).
Strait, J. B. & Lakatta, E. G. Aging-associated cardiovascular changes and their relationship to heart failure. Heart Fail Clin. 8, 143–146 (2012).
Herron, T. J., Lee, P. & Jalife, J. Optical imaging of voltage and calcium in cardiac cells & tissues. Circ. Res. 110, 609–623 (2012).
Unal, M. et al. Micro and nano-scale technologies for cell mechanics. Nanobiomedicine 1, 1–5 (2014).
Spira, M. E. & Hai, A. Multi-electrode array technologies for neuroscience and cardiology. Nat. Nanotechnol. 8, 83–94 (2013).
Fendyur, A. & Spira, M. E. Toward on-chip, in-cell recordings from cultured cardiomyocytes by arrays of gold mushroom-shaped microelectrodes. Front. Neuroeng. 5, 21 (2012).
Xie, C., Lin, Z. L., Hanson, L., Cui, Y. & Cui, B. X. Intracellular recording of action potentials by nanopillar electroporation. Nat. Nanotechnol. 7, 185–190 (2012).
Jahed, Z. et al. Nanocrown electrodes for parallel and robust intracellular recording of cardiomyocytes. Nat. Commun. 13, 2253 (2022).
Gu, Y. et al. Three-dimensional transistor arrays for intra-and inter-cellular recording. Nat. Nanotechnol. 17, 292–300 (2022).
Abbott, J., Ye, T. Y., Ham, D. & Park, H. Optimizing nanoelectrode arrays for scalable intracellular electrophysiology. Acc. Chem. Res. 51, 600–608 (2018).
Tian, B. et al. Macroporous nanowire nanoelectronic scaffolds for synthetic tissues. Nat. Mater. 11, 986–994 (2012).
Feiner, R. et al. Engineered hybrid cardiac patches with multifunctional electronics for online monitoring and regulation of tissue function. Nat. Mater. 15, 679–685 (2016).
Dai, X. C., Zhou, W., Gao, T., Liu, J. & Lieber, C. M. Three-dimensional mapping and regulation of action potential propagation in nanoelectronics-innervated tissues. Nat. Nanotechnol. 11, 776–782 (2016).
Li, Q. et al. Cyborg organoids: implantation of nanoelectronics via organogenesis for tissue-wide electrophysiology. Nano Lett. 19, 5781–5789 (2019).
Kalmykov, A. et al. Organ-on-e-chip: three-dimensional self-rolled biosensor array for electrical interrogations of human electrogenic spheroids. Sci. Adv. 5, eaax0729 (2019).
Le Floch, P. et al. Stretchable mesh nanoelectronics for 3D single-cell chronic electrophysiology from developing brain organoids. Adv. Mater. 34, e2106829 (2022).
Lin, Z. W. et al. Tissue-embedded stretchable nanoelectronics reveal endothelial cell-mediated electrical maturation of human 3D cardiac microtissues. Sci. Adv. 9, eade8513 (2023).
Hong, G. S. & Lieber, C. M. Novel electrode technologies for neural recordings. Nat. Rev. Neurosci. 20, 330–345 (2019).
Liu, R., Zhao, S. Y. & Liu, J. From lithographically patternable to genetically patternable electronic materials for miniaturized, scalable, and soft implantable bioelectronics to interface with nervous and cardiac systems. ACS Appl. Electron. Mater. 3, 101–118 (2021).
Tang, X., He, Y. & Liu, J. Soft bioelectronics for cardiac interfaces. Biophys. Rev. 3, 011301 (2022).
Lind, J. U. et al. Instrumented cardiac microphysiological devices via multimaterial three-dimensional printing. Nat. Mater. 16, 303–308 (2017).
Desbiolles, B. X. E. et al. Volcano-shaped scanning probe microscopy probe for combined force-electrogram recordings from excitable cells. Nano Lett. 20, 4520–4529 (2020).
Qian, F. et al. Simultaneous electrical recording of cardiac electrophysiology and contraction on chip. Lab. Chip 17, 1732–1739 (2017).
Gao, H. Y. et al. Bioinspired two-in-one nanotransistor sensor for the simultaneous measurements of electrical and mechanical cellular responses. Sci. Adv. 8, eabn2485 (2022).
Gao, H. Y. et al. Deterministic assembly of three-dimensional suspended nanowire structures. Nano Lett. 19, 5647–5652 (2019).
Cohen-Karni, T., Qing, Q., Li, Q., Fang, Y. & Lieber, C. M. Graphene and nanowire transistors for cellular interfaces and electrical recording. Nano Lett. 10, 1098–1102 (2010).
Irani, F. S. et al. Graphene as a Piezoresistive material in strain sensing applications. Micromachines 13, 119 (2022).
Wang, H. Z. et al. Low-temperature copper bonding strategy with graphene interlayer. ACS Nano 12, 2395–2402 (2018).
Garcia-Cortadella, R. et al. Graphene active sensor arrays for long-term and wireless mapping of wide frequency band epicortical brain activity. Nat. Commun. 12, 211 (2021).
Lian, X. J. et al. Directed cardiomyocyte differentiation from human pluripotent stem cells by modulating Wnt/beta-catenin signaling under fully defined conditions. Nat. Protoc. 8, 162–175 (2013).
Zuppinger, C. 3D cardiac cell culture: a critical review of current technologies and applications. Front. Cardiovasc. Med. 6, 87 (2019).
Lyra-Leite, D. M. et al. A review of protocols for human iPSC culture, cardiac differentiation, subtype-specification, maturation, and direct reprogramming. STAR Protoc. 3, 101560 (2022).
Wang, L. et al. Global strain-induced scalar potential in graphene devices. Commun. Phys. 4, 147 (2021).
Guo, Y. X. & Pu, W. L. T. Cardiomyocyte maturation new phase in development. Circ. Res. 126, 1086–1106 (2020).
Nunes, S. S. et al. Biowire: a platform for maturation of human pluripotent stem cell-derived cardiomyocytes. Nat. Methods 10, 781–787 (2013).
Ergir, E. et al. Generation and maturation of human iPSC-derived 3D organotypic cardiac microtissues in long-term culture. Sci. Rep. 12, 17409 (2022).
Buchanan, J. W., Saito, T. & Gettes, L. S. The effects of antiarrhythmic drugs, stimulation frequency, and potassium-induced resting membrane-potential changes on conduction-velocity and Dv/Dtmax in Guinea-Pig Myocardium. Circ. Res 56, 696–703 (1985).
Davis, J. P. & Tikunova, S. B. Ca2+ exchange with troponin C and cardiac muscle dynamics. Cardiovasc. Res. 77, 619–626 (2008).
Brini, M. & Carafoli, E. The plasma membrane Ca2+ATPase and the plasma membrane sodium calcium exchanger cooperate in the regulation of cell calcium. Cold Spring Harb. Perspect. Biol. 3, a004168 (2011).
Stevens, J. L. & Baker, T. K. The future of drug safety testing: expanding the view and narrowing the focus. Drug Discov. Today 14, 162–167 (2009).
Kovács, M., Tóth, J., Hetényi, C., Málnási-Csizmadia, A. & Sellers, J. R. Mechanism of blebbistatin inhibition of myosin II. J. Biol. Chem. 279, 35557–35563 (2004).
Lee, K. S. & Tsien, R. W. Mechanism of Calcium-Channel blockade by Verapamil, D600, diltiazem and nitrendipine in single dialyzed heart-cells. Nature 302, 790–794 (1983).
Lee, H. A., Hyun, S. A., Park, S. G., Kim, K. S. & Kim, S. J. Comparison of electrophysiological effects of calcium channel blockers on cardiac repolarization. Korean J. Physiol. Pha 20, 119–127 (2016).
Roden, D. M. Pharmacology and toxicology of Nav1.5-class 1 anti-arrhythmic drugs. Card. Electrophysiol. Clin. 6, 695–704 (2014).
Michel, D., Wegener, J. W. & Nawrath, H. Effects of quinine and quinidine on the transient outward and on the L-type Ca2+ current in rat ventricular cardiomyocytes. Pharmacology 65, 187–192 (2002).
Zhang, S. et al. Identification of the molecular basis of doxorubicin-induced cardiotoxicity. Nat. Med. 18, 1639–1642 (2012).
Rawat, P. S., Jaiswal, A., Khurana, A., Bhatti, J. S. & Navik, U. Doxorubicin-induced cardiotoxicity: an update on the molecular mechanism and novel therapeutic strategies for effective management. Biomed. Pharmacother. 139, 111708 (2021).
Prabhu, S. D. & Frangogiannis, N. G. The biological basis for cardiac repair after myocardial infarction from inflammation to fibrosis. Circ. Res. 119, 91–112 (2016).
Rumsey, W. L., Pawlowski, M., Lejavardi, N. & Wilson, D. F. Oxygen-pressure distribution in the heart in-vivo and evaluation of the ischemic border zone. Am. J. Physiol. 266, H1676–H1680 (1994).
Rexius-Hall, M. L. et al. A myocardial infarct border-zone-on-a-chip demonstrates distinct regulation of cardiac tissue function by an oxygen gradient. Sci. Adv. 8, eabn7097 (2022).
Richards, D. J. et al. Human cardiac organoids for the modelling of myocardial infarction and drug cardiotoxicity. Nat. Biomed. Eng. 4, 446–462 (2020).
Saint, D. A. The role of the persistent Na+ current during cardiac ischemia and hypoxia. J. Cardiovasc. Electrophysiol. 17, S96–S103 (2006).
Liu, H. T. et al. Heart-on-a-chip model with integrated extra- and intracellular bioelectronics for monitoring cardiac electrophysiology under acute hypoxia. Nano Lett. 20, 2585–2593 (2020).
Cascio, W. E., Yang, H., Muller-Borer, B. J. & Johnson, T. A. Ischemia-induced arrhythmia: the role of connexins, gap junctions, and attendant changes in impulse propagation. J. Electrocardiol. 38, 55–59 (2005).
Yeung, C. K. et al. To establish a pharmacological experimental platform for the study of cardiac hypoxia using the microelectrode array. J. Pharm. Toxicol. Methods 59, 146–152 (2009).
Klabunde, R. E. Cardiac electrophysiology: normal and ischemic ionic currents and the ECG. Adv. Physiol. Educ. 41, 29–37 (2017).
Sah, R. et al. Regulation of cardiac excitation-contraction coupling by action potential repolarization: role of the transient outward potassium current (Ito). J. Physiol. 546, 5–18 (2003).
Nimbalkar, S. et al. Graphene on glassy carbon microelectrodes demonstrate long-term structural and functional stability in neurophysiological recording and stimulation. J. Neural Eng. 18, 056035 (2021).
Bonaccini Calia, A. et al. Full-bandwidth electrophysiology of seizures and epileptiform activity enabled by flexible graphene microtransistor depth neural probes. Nat. Nanotechnol. 17, 301–309 (2022).
Fu, T. M. et al. Stable long-term chronic brain mapping at the single-neuron level. Nat. Methods 13, 875–882 (2016).
Zhao, S. et al. Tracking neural activity from the same cells during the entire adult life of mice. Nat. Neurosci. 26, 696–710 (2023).
Yang, X. et al. Bioinspired neuron-like electronics. Nat. Mater. 18, 510–517 (2019).
Schuhmann, T. G., Yao, J., Hong, G., Fu, T. M. & Lieber, C. M. Syringe-injectable electronics with a plug-and-play Input/Output interface. Nano Lett. 17, 5836–5842 (2017).
Zhang, Y. S. et al. Multisensor-integrated organs-on-chips platform for automated and continual in situ monitoring of organoid behaviors. Proc. Natl Acad. Sci. USA 114, E2293–E2302 (2017).
Acknowledgements
J.Y. acknowledges support from the Army Research Office W911NF2210027 and National Science Foundation (NSF) CBET-1844904. J.Y. also acknowledges supports from National Institutes of Health (NIH) 5R21 EB030216, NSF ECCS-1917630, and NSF DMR-2027102. Z.W. and J.K. acknowledges support from the Semiconductor Research Corporation Center 7 in JUMP 2.0 (award no. 145105-21913). Y.S. acknowledges support from NSF CMMI 1846866 and NIH R01DK129990. X.L. acknowledges support from a Link Foundation Energy Fellowship. Part of the device fabrication work was conducted in the clean room of the Center for Hierarchical Manufacturing (CHM), an NSF Nanoscale Science and Engineering Center (NSEC) located at the University of Massachusetts, Amherst.
Author information
Authors and Affiliations
Contributions
H.G. and J.Y. conceived the project. H.G. performed experiments in device fabrication, characterizations, cell culture, imaging, and in vitro electrical recordings. Z.W. and J.K. performed graphene synthesis. F.Y. and Y.S. contributed to cell supply and helped in cell culture. S.W. and X.W. assisted in device fabrication and cell culture. Q.Z. and X.L. assisted in microfluidic chamber fabrication. H.G. and J.Y. wrote the paper. All authors discussed the results and implications and commented on the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks Jong Seob Choi, Guosong Hong, Lukas Bongartz and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Gao, H., Wang, Z., Yang, F. et al. Graphene-integrated mesh electronics with converged multifunctionality for tracking multimodal excitation-contraction dynamics in cardiac microtissues. Nat Commun 15, 2321 (2024). https://doi.org/10.1038/s41467-024-46636-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-024-46636-7
This article is cited by
-
Atomically thin bioelectronics
Nature Reviews Materials (2024)