Indirect Correlative Light and Electron Microscopy (iCLEM): A Novel Pipeline for Multiscale Quantification of Structure From Molecules to Organs.

Correlative light and electron microscopy (CLEM) methods are powerful methods that combine molecular organization (from light microscopy) with ultrastructure (from electron microscopy). However, CLEM methods pose high cost/difficulty barriers to entry and have very low experimental throughput. Therefore, we have developed an indirect correlative light and electron microscopy (iCLEM) pipeline to sidestep the rate-limiting steps of CLEM (i.e., preparing and imaging the same samples on multiple microscopes) and correlate multiscale structural data gleaned from separate samples imaged using different modalities by exploiting biological structures identifiable by both light and electron microscopy as intrinsic fiducials. We demonstrate here an application of iCLEM, where we utilized gap junctions and mechanical junctions between muscle cells in the heart as intrinsic fiducials to correlate ultrastructural measurements from transmission electron microscopy (TEM), and focused ion beam scanning electron microscopy (FIB-SEM) with molecular organization from confocal microscopy and single molecule localization microscopy (SMLM). We further demonstrate how iCLEM can be integrated with computational modeling to discover structure-function relationships. Thus, we present iCLEM as a novel approach that complements existing CLEM methods and provides a generalizable framework that can be applied to any set of imaging modalities, provided suitable intrinsic fiducials can be identified.

Emergent activity, heterogeneity, and robustness in a calcium feedback model of the sinoatrial node.

The sinoatrial node (SAN) is the primary pacemaker of the heart. SAN activity emerges at an early point in life and maintains a steady rhythm for the lifetime of the organism. The ion channel composition and currents of SAN cells can be influenced by a variety of factors. Therefore, the emergent activity and long-term stability imply some form of dynamical feedback control of SAN activity. We adapt a recent feedback model-previously utilized to describe control of ion conductances in neurons-to a model of SAN cells and tissue. The model describes a minimal regulatory mechanism of ion channel conductances via feedback between intracellular calcium and an intrinsic target calcium level. By coupling a SAN cell to the calcium feedback model, we show that spontaneous electrical activity emerges from quiescence and is maintained at steady state. In a 2D SAN tissue model, spatial variability in intracellular calcium targets lead to significant, self-organized heterogeneous ion channel expression and calcium transients throughout the tissue. Furthermore, multiple pacemaking regions appear, which interact and lead to time-varying cycle length, demonstrating that variability in heart rate is an emergent property of the feedback model. Finally, we demonstrate that the SAN tissue is robust to the silencing of leading cells or ion channel knockouts. Thus, the calcium feedback model can reproduce and explain many fundamental emergent properties of activity in the SAN that have been observed experimentally based on a minimal description of intracellular calcium and ion channel regulatory networks.

Modeling incomplete penetrance in long QT syndrome type 3 through ion channel heterogeneity: an in silico population study.

Many cardiac diseases are characterized by an increased late sodium current, including heart failure, hypertrophic cardiomyopathy, and inherited long QT syndrome type 3 (LQT3). The late sodium current in LQT3 is caused by a gain-of-function mutation in the voltage-gated sodium channel Nav1.5. Despite a well-defined genetic cause of LQT3, treatment remains inconsistent because of incomplete penetrance of the mutation and variability of antiarrhythmic efficacy. Here, we investigate the relationship between LQT3-associated mutation incomplete penetrance and variability in ion channel expression, simulating a population of 1,000 individuals with the O'Hara-Rudy model of the human ventricular myocyte. We first simulate healthy electrical activity (i.e., in the absence of a mutation) and then incorporate heterozygous expression for three LQT3-associated mutations (Y1795C, I1768V, and ΔKPQ), to directly compare the effects of each mutation on individuals across a diverse population. For all mutations, we find that susceptibility, defined by either the presence of an early afterdepolarization (EAD) or prolonged action potential duration (APD), primarily depends on the balance between the conductance of IKr and INa, for which individuals with a higher IKr-to-INa ratio are less susceptible. Furthermore, we find distinct differences across the population, observing individuals susceptible to zero, one, two, or all three mutations. Individuals tend to be less susceptible with an appropriate balance of repolarizing currents, typically via increased IKs or IK1. Interestingly, the more critical repolarizing current is mutation specific. We conclude that the balance between key currents plays a significant role in mutant-specific presentation of the disease phenotype in LQT3.NEW & NOTEWORTHY An in silico population approach investigates the relationship between variability in ion channel expression and gain-of-function mutations in the voltage-gated sodium channel associated with the congenital disorder long QT syndrome type 3 (LQT3). We find that ion channel variability can contribute to incomplete penetrance of the mutation, with mutant-specific differences in ion channel conductances leading to susceptibility to proarrhythmic action potential duration prolongation or early afterdepolarizations.

Long-term changes in heart rate and electrical remodeling contribute to alternans formation in heart failure: a patient-specific in silico study.

Individuals with chronic heart failure (CHF) have an increased risk of ventricular arrhythmias, which has been linked to pathological cellular remodeling and may also be mediated by changes in heart rate. Heart rate typically fluctuates on a timescale ranging from seconds to hours, termed heart rate variability (HRV). This variability is reduced in CHF, and this HRV reduction is associated with a greater risk for arrhythmias. Furthermore, variations in heart rate influence the formation of proarrhythmic alternans, a beat-to-beat alternation in the action potential duration (APD), or intracellular calcium (Ca). In this study, we investigate how long-term changes in heart rate and electrical remodeling associated with CHF influence alternans formation. We measure key statistical properties of the RR-interval sequences from ECGs of individuals with normal sinus rhythm (NSR) and CHF. Patient-specific RR-interval sequences and synthetic sequences (randomly generated to mimicking these statistical properties) are used as the pacing protocol for a discrete time-coupled map model that governs APD and intracellular Ca handling of a single cardiac myocyte, modified to account for pathological electrical remodeling in CHF. Patient-specific simulations show that beat-to-beat differences in APD vary temporally in both populations, with alternans formation more prevalent in CHF. Parameter studies using synthetic sequences demonstrate that increasing the autocorrelation time or mean RR-interval reduces APD alternations, whereas increasing the RR-interval standard deviation leads to higher alternans magnitudes. Importantly, we find that although both the CHF-associated changes in heart rate and electrical remodeling influence alternans formation, variations in heart rate may be more influential.NEW & NOTEWORTHY Using patient-specific data, we show that both the changes in heart rate and electrical remodeling associated with chronic heart failure influence the formation of proarrhythmic alternans in the heart.

Ephaptic Coupling as a Resolution to the Paradox of Action Potential Wave Speed and Discordant Alternans Spatial Scales in the Heart.

Previous computer simulations have suggested that existing models of action potential wave propagation in the heart are not consistent with observed wave propagation behavior. Specifically, computer models cannot simultaneously reproduce the rapid wave speeds and small spatial scales of discordant alternans patterns measured experimentally in the same simulation. The discrepancy is important, because discordant alternans can be a key precursor to the development of abnormal and dangerous rapid rhythms in the heart. In this Letter, we show that this paradox can be resolved by allowing so-called ephaptic coupling to play a primary role in wave front propagation in place of conventional gap-junction coupling. With this modification, physiological wave speeds and small discordant alternans spatial scales both occur with gap-junction resistance values that are more in line with those observed in experiments. Our theory thus also provides support to the hypothesis that ephaptic coupling plays an important role in normal wave propagation.

A data-assimilation approach to predict population dynamics during epithelial-mesenchymal transition.

Epithelial-mesenchymal transition (EMT) is a biological process that plays a central role in embryonic development, tissue regeneration, and cancer metastasis. Transforming growth factor-ß (TGFß) is a potent inducer of this cellular transition, comprising transitions from an epithelial state to partial or hybrid EMT state(s), to a mesenchymal state. Recent experimental studies have shown that, within a population of epithelial cells, heterogeneous phenotypical profiles arise in response to different time- and TGFß dose-dependent stimuli. This offers a challenge for computational models, as most model parameters are generally obtained to represent typical cell responses, not necessarily specific responses nor to capture population variability. In this study, we applied a data-assimilation approach that combines limited noisy observations with predictions from a computational model, paired with parameter estimation. Synthetic experiments mimic the biological heterogeneity in cell states that is observed in epithelial cell populations by generating a large population of model parameter sets. Analysis of the parameters for virtual epithelial cells with biologically significant characteristics (e.g., EMT prone or resistant) illustrates that these sub-populations have identifiable critical model parameters. We perform a series of in silico experiments in which a forecasting system reconstructs the EMT dynamics of each virtual cell within a heterogeneous population exposed to time-dependent exogenous TGFß dose and either an EMT-suppressing or EMT-promoting perturbation. We find that estimating population-specific critical parameters significantly improved the prediction accuracy of cell responses. Thus, with appropriate protocol design, we demonstrate that a data-assimilation approach successfully reconstructs and predicts the dynamics of a heterogeneous virtual epithelial cell population in the presence of physiological model error and parameter uncertainty.

Computational modeling of aberrant electrical activity following remuscularization with intramyocardially injected pluripotent stem cell-derived cardiomyocytes.

Acute engraftment arrhythmias (EAs) remain a serious complication of remuscularization therapy. Preliminary evidence suggests that a focal source underlies these EAs stemming from the automaticity of immature pluripotent stem cell-derived cardiomyocytes (PSC-CMs) in nascent myocardial grafts. How these EAs arise though during early engraftment remains unclear. In a series of in silico experiments, we probed the origin of EAs-exploring aspects of altered impulse formation and altered impulse propagation within nascent PSC-CM grafts and at the host-graft interface. To account for poor gap junctional coupling during early PSC-CM engraftment, the voltage dependence of gap junctions and the possibility of ephaptic coupling were incorporated. Inspired by cardiac development, we also studied the contributions of another feature of immature PSC-CMs, circumferential sodium channel (NaCh) distribution in PSC-CMs. Ectopic propagations emerged from nascent grafts of immature PSC-CMs at a rate of <96 bpm. Source-sink effects dictated this rate and contributed to intermittent capture between host and graft. Moreover, ectopic beats emerged from dynamically changing sites along the host-graft interface. The latter arose in part because circumferential NaCh distribution in PSC-CMs contributed to preferential conduction slowing and block of electrical impulses from host to graft myocardium. We conclude that additional mechanisms, in addition to focal ones, contribute to EAs and recognize that their relative contributions are dynamic across the engraftment process.

Automaticity in ventricular myocyte cell pairs with ephaptic and gap junction coupling.

Spontaneous electrical activity, or automaticity, in the heart is required for normal physiological function. However, irregular automaticity, in particular, originating from the ventricles, can trigger life-threatening cardiac arrhythmias. Thus, understanding mechanisms of automaticity and synchronization is critical. Recent work has proposed that excitable cells coupled via a shared narrow extracellular cleft can mediate coupling, i.e., ephaptic coupling, that promotes automaticity in cell pairs. However, the dynamics of these coupled cells incorporating both ephaptic and gap junction coupling has not been explored. Here, we show that automaticity and synchronization robustly emerges via a Hopf bifurcation from either (i) increasing the fraction of inward rectifying potassium channels (carrying the IK1 current) at the junctional membrane or (ii) by decreasing the cleft volume. Furthermore, we explore how heterogeneity in the fraction of potassium channels between coupled cells can produce automaticity of both cells or neither cell, or more rarely in only one cell (i.e., automaticity without synchronization). Interestingly, gap junction coupling generally has minor effects, with only slight changes in regions of parameter space of automaticity. This work provides insight into potentially new mechanisms that promote spontaneous activity and, thus, triggers for arrhythmias in ventricular tissue.

Initiation and entrainment of multicellular automaticity via diffusion limited extracellular domains.

Electrically excitable cells often spontaneously and synchronously depolarize in vitro and in vivo preparations. It remains unclear how cells entrain and autorhythmically activate above the intrinsic mean activation frequency of isolated cells with or without pacemaking mechanisms. Recent studies suggest that cyclic ion accumulation and depletion in diffusion-limited extracellular volumes modulate electrophysiology by ephaptic mechanisms (nongap junction or synaptic coupling). This report explores how potassium accumulation and depletion in a restricted extracellular domain induces spontaneous action potentials in two different computational models of excitable cells without gap junctional coupling: Hodgkin-Huxley and Luo-Rudy. Importantly, neither model will spontaneously activate on its own without external stimuli. Simulations demonstrate that cells sharing a diffusion-limited extracellular compartment can become autorhythmic and entrained despite intercellular electrical heterogeneity. Autorhythmic frequency is modulated by the cleft volume and potassium fluxes through the cleft. Additionally, inexcitable cells can suppress or induce autorhythmic activity in an excitable cell via a shared cleft. Diffusion-limited shared clefts can also entrain repolarization. Critically, this model predicts a mechanism by which diffusion-limited shared clefts can initiate, entrain, and modulate multicellular automaticity in the absence of gap junctions.

Mechanisms underlying age-associated manifestation of cardiac sodium channel gain-of-function.

Cardiac action potentials are initiated by sodium ion (Na+) influx through voltage-gated Na+ channels. Na+ channel gain-of-function (GOF) can arise in inherited conditions due to mutations in the gene encoding the cardiac Na+ channel, such as Long QT syndrome type 3 (LQT3). LQT3 can be a "concealed" disease, as patients with LQT3-associated mutations can remain asymptomatic until later in life; however, arrhythmias can also arise early in life in LQT3 patients, demonstrating a complex age-associated manifestation. We and others recently demonstrated that cardiac Na+ channels preferentially localize at the intercalated disc (ID) in adult cardiac tissue, which facilitates ephaptic coupling and formation of intercellular Na+ nanodomains that regulate pro-arrhythmic early afterdepolarization (EAD) formation in tissue with Na+ channel GOF. Several properties related to ephaptic coupling vary with age, such as cell size and Na+ channel and gap junction (GJ) expression and distribution: neonatal cells have immature IDs, with Na+ channels and GJs primarily diffusively distributed, while adult myocytes have mature IDs with preferentially localized Na+ channels and GJs. Here, we perform an in silico study varying critical age-dependent parameters to investigate mechanisms underlying age-associated manifestation of Na+ channel GOF in a model of guinea pig cardiac tissue. Simulations predict that total Na+ current conductance is a critical factor in action potential duration (APD) prolongation. We find a complex cell size/ Na+ channel expression relationship: increases in cell size (without concurrent increases in Na+ channel expression) suppress EAD formation, while increases in Na+ channel expression (without concurrent increases in cell size) promotes EAD formation. Finally, simulations with neonatal and early age-associated parameters predict normal APD with minimal dependence on intercellular cleft width; however, variability in cellular properties can lead to EADs presenting in early developmental stages. In contrast, for adult-associated parameters, EAD formation is highly dependent on cleft width, consistent with a mechanism underlying the age-associated manifestation of the Na+ channel GOF.

Hypernatremia and intercalated disc edema synergistically exacerbate long-QT syndrome type 3 phenotype.

Cardiac voltage-gated sodium channel gain-of-function prolongs repolarization in the long-QT syndrome type 3 (LQT3). Previous studies suggest that narrowing the perinexus within the intercalated disc, leading to rapid sodium depletion, attenuates LQT3-associated action potential duration (APD) prolongation. However, it remains unknown whether extracellular sodium concentration modulates APD prolongation during sodium channel gain-of-function. We hypothesized that elevated extracellular sodium concentration and widened perinexus synergistically prolong APD in LQT3. LQT3 was induced with sea anemone toxin (ATXII) in Langendorff-perfused guinea pig hearts (n = 34). Sodium concentration was increased from 145 to 160 mM. Perinexal expansion was induced with mannitol or the sodium channel ß1-subunit adhesion domain antagonist (ßadp1). Epicardial ventricular action potentials were optically mapped. Individual and combined effects of varying clefts and sodium concentrations were simulated in a computational model. With ATXII, both mannitol and ßadp1 significantly widened the perinexus and prolonged APD, respectively. The elevated sodium concentration alone significantly prolonged APD as well. Importantly, the combination of elevated sodium concentration and perinexal widening synergistically prolonged APD. Computational modeling results were consistent with animal experiments. Concurrently elevating extracellular sodium and increasing intercalated disc edema prolongs repolarization more than the individual interventions alone in LQT3. This synergistic effect suggests an important clinical implication that hypernatremia in the presence of cardiac edema can markedly increase LQT3-associated APD prolongation. Therefore, to our knowledge, this is the first study to provide evidence of a tractable and effective strategy to mitigate LQT3 phenotype by means of managing sodium levels and preventing cardiac edema in patients.NEW & NOTEWORTHY This is the first study to demonstrate that the long-QT syndrome type 3 (LQT3) phenotype can be exacerbated or concealed by regulating extracellular sodium concentrations and/or the intercalated disc separation. The animal experiments and computational modeling in the current study reveal a critically important clinical implication: sodium dysregulation in the presence of edema within the intercalated disc can markedly increase the risk of arrhythmia in LQT3. These findings strongly suggest that maintaining extracellular sodium within normal physiological limits may be an effective and inexpensive therapeutic option for patients with congenital or acquired sodium channel gain-of-function diseases.

Distributed synthesis of sarcolemmal and sarcoplasmic reticulum membrane proteins in cardiac myocytes.

It is widely assumed that synthesis of membrane proteins, particularly in the heart, follows the classical secretory pathway with mRNA translation occurring in perinuclear regions followed by protein trafficking to sites of deployment. However, this view is based on studies conducted in less-specialized cells, and has not been experimentally addressed in cardiac myocytes. Therefore, we undertook direct experimental investigation of protein synthesis in cardiac tissue and isolated myocytes using single-molecule visualization techniques and a novel proximity-ligated in situ hybridization approach for visualizing ribosome-associated mRNA molecules for a specific protein species, indicative of translation sites. We identify here, for the first time, that the molecular machinery for membrane protein synthesis occurs throughout the cardiac myocyte, and enables distributed synthesis of membrane proteins within sub-cellular niches where the synthesized protein functions using local mRNA pools trafficked, in part, by microtubules. We also observed cell-wide distribution of membrane protein mRNA in myocardial tissue from both non-failing and hypertrophied (failing) human hearts, demonstrating an evolutionarily conserved distributed mechanism from mouse to human. Our results identify previously unanticipated aspects of local control of cardiac myocyte biology and highlight local protein synthesis in cardiac myocytes as an important potential determinant of the heart's biology in health and disease.

Immunofluorescence Image Feature Analysis and Phenotype Scoring Pipeline for Distinguishing Epithelial-Mesenchymal Transition.

Epithelial­mesenchymal transition (EMT) is an essential biological process, also implicated in pathological settings such as cancer metastasis, in which epithelial cells transdifferentiate into mesenchymal cells. We devised an image analysis pipeline to distinguish between tissues comprised of epithelial and mesenchymal cells, based on extracted features from immunofluorescence images of differing biochemical markers. Mammary epithelial cells were cultured with 0 (control), 2, 4, or 10 ng/mL TGF-ß1, a well-established EMT-inducer. Cells were fixed, stained, and imaged for E-cadherin, actin, fibronectin, and nuclei via immunofluorescence microscopy. Feature selection was performed on different combinations of individual cell markers using a Bag-of-Features extraction. Control and high-dose images comprised the training data set, and the intermediate dose images comprised the testing data set. A feature distance analysis was performed to quantify differences between the treatment groups. The pipeline was successful in distinguishing between control (epithelial) and the high-dose (mesenchymal) groups, as well as demonstrating progress along the EMT process in the intermediate dose groups. Validation using quantitative PCR (qPCR) demonstrated that biomarker expression measurements were well-correlated with the feature distance analysis. Overall, we identified image pipeline characteristics for feature extraction and quantification of immunofluorescence images to distinguish progression of EMT.

Cell Fate Forecasting: A Data-Assimilation Approach to Predict Epithelial-Mesenchymal Transition.

Epithelial-mesenchymal transition (EMT) is a fundamental biological process that plays a central role in embryonic development, tissue regeneration, and cancer metastasis. Transforming growth factor-ß (TGFß) is a potent inducer of this cellular transition, which is composed of transitions from an epithelial state to intermediate or partial EMT state(s) to a mesenchymal state. Using computational models to predict cell state transitions in a specific experiment is inherently difficult for reasons including model parameter uncertainty and error associated with experimental observations. In this study, we demonstrate that a data-assimilation approach using an ensemble Kalman filter, which combines limited noisy observations with predictions from a computational model of TGFß-induced EMT, can reconstruct the cell state and predict the timing of state transitions. We used our approach in proof-of-concept "synthetic" in silico experiments, in which experimental observations were produced from a known computational model with the addition of noise. We mimic parameter uncertainty in in vitro experiments by incorporating model error that shifts the TGFß doses associated with the state transitions and reproduces experimentally observed variability in cell state by either shifting a single parameter or generating "populations" of model parameters. We performed synthetic experiments for a wide range of TGFß doses, investigating different cell steady-state conditions, and conducted parameter studies varying properties of the data-assimilation approach including the time interval between observations and incorporating multiplicative inflation, a technique to compensate for underestimation of the model uncertainty and mitigate the influence of model error. We find that cell state can be successfully reconstructed and the future cell state predicted in synthetic experiments, even in the setting of model error, when experimental observations are performed at a sufficiently short time interval and incorporate multiplicative inflation. Our study demonstrates the feasibility and utility of a data-assimilation approach to forecasting the fate of cells undergoing EMT.

Intercellular Sodium Regulates Repolarization in Cardiac Tissue with Sodium Channel Gain of Function.

In cardiac myocytes, action potentials are initiated by an influx of sodium (Na+) ions via voltage-gated Na+ channels. Na+ channel gain of function (GOF), arising in both inherited conditions associated with mutation in the gene encoding the Na+ channel and acquired conditions associated with heart failure, ischemia, and atrial fibrillation, enhance Na+ influx, generating a late Na+ current that prolongs action potential duration (APD) and triggering proarrhythmic early afterdepolarizations (EADs). Recent studies have shown that Na+ channels are highly clustered at the myocyte intercalated disk, facilitating formation of Na+ nanodomains in the intercellular cleft between cells. Simulations from our group have recently predicted that narrowing the width of the intercellular cleft can suppress APD prolongation and EADs in the presence of Na+ channel mutations because of increased intercellular cleft Na+ ion depletion. In this study, we investigate the effects of modulating multiple extracellular spaces, specifically the intercellular cleft and bulk interstitial space, in a novel computational model and experimentally via osmotic agents albumin, dextran 70, and mannitol. We perform optical mapping and transmission electron microscopy in a drug-induced (sea anemone toxin, ATXII) Na+ channel GOF isolated heart model and modulate extracellular spaces via osmotic agents. Single-cell patch-clamp experiments confirmed that the osmotic agents individually do not enhance late Na+ current. Both experiments and simulations are consistent with the conclusion that intercellular cleft narrowing or expansion regulates APD prolongation; in contrast, modulating the bulk interstitial space has negligible effects on repolarization. Thus, we predict that intercellular cleft Na+ nanodomain formation and collapse critically regulates cardiac repolarization in the setting of Na+ channel GOF.

Dual regulation by subcellular calcium heterogeneity and heart rate variability on cardiac electromechanical dynamics.

Heart rate constantly varies under physiological conditions, termed heart rate variability (HRV), and in clinical studies, low HRV is associated with a greater risk of cardiac arrhythmias. Prior work has shown that HRV influences the temporal patterns of electrical activity, specifically the formation of pro-arrhythmic alternans, a beat-to-beat alternation in the action potential duration (APD), or intracellular calcium (Ca) levels. We previously showed that HRV may be anti-arrhythmic by disrupting APD and Ca alternations in a homogeneous cardiac myocyte. Here, we expand on our previous work, incorporating variation in subcellular Ca handling (also known to influence alternans) into a nonlinear map model of a cardiac myocyte composed of diffusively coupled Ca release units (CRUs). Ca-related parameters and initial conditions of each CRU are varied to mimic subcellular Ca heterogeneity, and a stochastic pacing sequence reproduces HRV. We find that subcellular Ca heterogeneity promotes the formation of spatially discordant subcellular alternans patterns, which decreases whole cell Ca and APD alternation for low and moderate HRV, while high subcellular Ca heterogeneity and HRV both promote electromechanical desynchronization. Finally, we find that for low and moderate HRV, both the specific subcellular Ca-related parameters and the pacing sequences influence measures of electromechanical dynamics, while for high HRV, these measures depend predominantly on the pacing sequence. Our results suggest that pro-arrhythmic subcellular discordant alternans tend to form for low levels of HRV, while high HRV may be anti-arrhythmic due to mitigated influence from subcellular Ca heterogeneity and desynchronization of APD from Ca instabilities.

Death Certification in Northern Alberta: Error Occurrence Rate and Educational Intervention.

Errors in death certification can directly affect the decedent's survivors and the public register. We assessed the effectiveness of an educational seminar targeting frequent and important errors identified by local death certificate (DC) evaluation. Retrospective review of 1500 DCs categorized errors and physician specialty. A 60-minute didactic/case-based seminar was subsequently designed for family medicine physician (FAM) participants, with administration of presurvey, immediate post, and 2-month postsurveys. Most DCs were completed by FAM (73%), followed by internists (18%) and surgeons (3%). Error occurrence (EO) rate ranged between 32 and 75% across all specialities. Family medicine physician experienced in palliative care had the lowest EO rate (32%), significantly lower (P < 0.001) than FAM without interest in palliative care (62%), internal medicine (62%), and surgery (75%). Common errors were use of abbreviations (26%), mechanism as underlying cause of death (23%), and no underlying cause of death recorded (22%). Presurvey participants (n = 72) had an overall EO rate of 72% (64% excluding formatting errors). Immediate postsurvey (n = 75) and 2-month postsurvey (n = 24) participants demonstrated significantly lower overall EO (34% and 24%, respectively), compared with the Pre-S (P < 0.05). A 60-minute seminar on death certification reduced EO rate with perceived long-term effects.

Mechanochemical Coupling and Junctional Forces during Collective Cell Migration.

Cell migration, a fundamental physiological process in which cells sense and move through their surrounding physical environment, plays a critical role in development and tissue formation, as well as pathological processes, such as cancer metastasis and wound healing. During cell migration, dynamics are governed by the bidirectional interplay between cell-generated mechanical forces and the activity of Rho GTPases, a family of small GTP-binding proteins that regulate actin cytoskeleton assembly and cellular contractility. These interactions are inherently more complex during the collective migration of mechanically coupled cells because of the additional regulation of cell-cell junctional forces. In this study, we adapted a recent minimal modeling framework to simulate the interactions between mechanochemical signaling in individual cells and interactions with cell-cell junctional forces during collective cell migration. We find that migration of individual cells depends on the feedback between mechanical tension and Rho GTPase activity in a biphasic manner. During collective cell migration, waves of Rho GTPase activity mediate mechanical contraction/extension and thus synchronization throughout the tissue. Further, cell-cell junctional forces exhibit distinct spatial patterns during collective cell migration, with larger forces near the leading edge. Larger junctional force magnitudes are associated with faster collective cell migration and larger tissue size. Simulations of heterogeneous tissue migration exhibit a complex dependence on the properties of both leading and trailing cells. Computational predictions demonstrate that collective cell migration depends on both the emergent dynamics and interactions between cellular-level Rho GTPase activity and contractility and multicellular-level junctional forces.