Digital twinning of cardiac electrophysiology for congenital heart disease.

In recent years, blending mechanistic knowledge with machine learning has had a major impact in digital healthcare. In this work, we introduce a computational pipeline to build certified digital replicas of cardiac electrophysiology in paediatric patients with congenital heart disease. We construct the patient-specific geometry by means of semi-automatic segmentation and meshing tools. We generate a dataset of electrophysiology simulations covering cell-to-organ level model parameters and using rigorous mathematical models based on differential equations. We previously proposed Branched Latent Neural Maps (BLNMs) as an accurate and efficient means to recapitulate complex physical processes in a neural network. Here, we employ BLNMs to encode the parametrized temporal dynamics of in silico 12-lead electrocardiograms (ECGs). BLNMs act as a geometry-specific surrogate model of cardiac function for fast and robust parameter estimation to match clinical ECGs in paediatric patients. Identifiability and trustworthiness of calibrated model parameters are assessed by sensitivity analysis and uncertainty quantification.

Uncertainty quantification of the wall thickness and stiffness in an idealized dissected aorta.

Personalized treatment informed by computational models has the potential to markedly improve the outcome for patients with a type B aortic dissection. However, existing computational models of dissected walls significantly simplify the characteristic false lumen, tears and/or material behavior. Moreover, the patient-specific wall thickness and stiffness cannot be accurately captured non-invasively in clinical practice, which inevitably leads to assumptions in these wall models. It is important to evaluate the impact of the corresponding uncertainty on the predicted wall deformations and stress, which are both key outcome indicators for treatment optimization. Therefore, a physiology-inspired finite element framework was proposed to model the wall deformation and stress of a type B aortic dissection at diastolic and systolic pressure. Based on this framework, 300 finite element analyses, sampled with a Latin hypercube, were performed to assess the global uncertainty, introduced by 4 uncertain wall thickness and stiffness input parameters, on 4 displacement and stress output parameters. The specific impact of each input parameter was estimated using Gaussian process regression, as surrogate model of the finite element framework, and a δ moment-independent analysis. The global uncertainty analysis indicated minor differences between the uncertainty at diastolic and systolic pressure. For all output parameters, the 4th quartile contained the major fraction of the uncertainty. The parameter-specific uncertainty analysis elucidated that the material stiffness and relative thickness of the dissected membrane were the respective main determinants of the wall deformation and stress. The uncertainty analysis provides insight into the effect of uncertain wall thickness and stiffness parameters on the predicted deformation and stress. Moreover, it emphasizes the need for probabilistic rather than deterministic predictions for clinical decision making in aortic dissections.

Digital twinning of cardiac electrophysiology for congenital heart disease.

In recent years, blending mechanistic knowledge with machine learning has had a major impact in digital healthcare. In this work, we introduce a computational pipeline to build certified digital replicas of cardiac electrophysiology in pediatric patients with congenital heart disease. We construct the patient-specific geometry by means of semi-automatic segmentation and meshing tools. We generate a dataset of electrophysiology simulations covering cell-to-organ level model parameters and utilizing rigorous mathematical models based on differential equations. We previously proposed Branched Latent Neural Maps (BLNMs) as an accurate and efficient means to recapitulate complex physical processes in a neural network. Here, we employ BLNMs to encode the parametrized temporal dynamics of in silico 12-lead electrocardiograms (ECGs). BLNMs act as a geometry-specific surrogate model of cardiac function for fast and robust parameter estimation to match clinical ECGs in pediatric patients. Identifiability and trustworthiness of calibrated model parameters are assessed by sensitivity analysis and uncertainty quantification.

Balancing the ventricular outputs of pulsatile total artificial hearts.

BACKGROUND: Maintaining balanced left and right cardiac outputs in a total artificial heart (TAH) is challenging due to the need for continuous adaptation to changing hemodynamic conditions. Proper balance in ventricular outputs of the left and right ventricles requires a preload-sensitive response and mechanisms to address the higher volumetric efficiency of the right ventricle. METHODS: This review provides a comprehensive overview of various methods used to balance left and right ventricular outputs in pulsatile total artificial hearts, categorized based on their actuation mechanism. RESULTS: Reported strategies include incorporating compliant materials and/or air cushions inside the ventricles, employing active control mechanisms to regulate ventricular filling state, and utilizing various shunts (such as hydraulic or intra-atrial shunts). Furthermore, reducing right ventricular stroke volume compared to the left often serves to balance the ventricular outputs. Individually controlled actuation of both ventricles in a pulsatile TAH seems to be the simplest and most effective way to achieve proper preload sensitivity and left-right output balance. Pneumatically actuated TAHs have the advantage to respond passively to preload changes. CONCLUSION: Therefore, a pneumatic TAH that comprises two individually actuated ventricles appears to be a more desirable option-both in terms of simplicity and efficacy-to respond to changing hemodynamic conditions.

Validating MRI-Derived Myocardial Stiffness Estimates Using In Vitro Synthetic Heart Models.

Impaired cardiac filling in response to increased passive myocardial stiffness contributes to the pathophysiology of heart failure. By leveraging cardiac MRI data and ventricular pressure measurements, we can estimate in vivo passive myocardial stiffness using personalized inverse finite element models. While it is well-known that this approach is subject to uncertainties, only few studies quantify the accuracy of these stiffness estimates. This lack of validation is, at least in part, due to the absence of ground truth in vivo passive myocardial stiffness values. Here, using 3D printing, we created soft, homogenous, isotropic, hyperelastic heart phantoms of varying geometry and stiffness and simulate diastolic filling by incorporating the phantoms into an MRI-compatible left ventricular inflation system. We estimate phantom stiffness from MRI and pressure data using inverse finite element analyses based on a Neo-Hookean model. We demonstrate that our identified softest and stiffest values of 215.7 and 512.3 kPa agree well with the ground truth of 226.2 and 526.4 kPa. Overall, our estimated stiffnesses revealed a good agreement with the ground truth ([Formula: see text] error) across all models. Our results suggest that MRI-driven computational constitutive modeling can accurately estimate synthetic heart material stiffnesses in the range of 200-500 kPa.

Sex Matters: A Comprehensive Comparison of Female and Male Hearts.

Cardiovascular disease in women remains under-diagnosed and under-treated. Recent studies suggest that this is caused, at least in part, by the lack of sex-specific diagnostic criteria. While it is widely recognized that the female heart is smaller than the male heart, it has long been ignored that it also has a different microstructural architecture. This has severe implications on a multitude of cardiac parameters. Here, we systematically review and compare geometric, functional, and structural parameters of female and male hearts, both in the healthy population and in athletes. Our study finds that, compared to the male heart, the female heart has a larger ejection fraction and beats at a faster rate but generates a smaller cardiac output. It has a lower blood pressure but produces universally larger contractile strains. Critically, allometric scaling, e.g., by lean body mass, reduces but does not completely eliminate the sex differences between female and male hearts. Our results suggest that the sex differences in cardiac form and function are too complex to be ignored: the female heart is not just a small version of the male heart. When using similar diagnostic criteria for female and male hearts, cardiac disease in women is frequently overlooked by routine exams, and it is diagnosed later and with more severe symptoms than in men. Clearly, there is an urgent need to better understand the female heart and design sex-specific diagnostic criteria that will allow us to diagnose cardiac disease in women equally as early, robustly, and reliably as in men. Systematic Review Registration: https://livingmatter.stanford.edu/.

Growth and remodeling in the pulmonary autograft: Computational evaluation using kinematic growth models and constrained mixture theory.

Computational investigations of how soft tissues grow and remodel are gaining more and more interest and several growth and remodeling theories have been developed. Roughly, two main groups of theories for soft tissues can be distinguished: kinematic-based growth theory and theories based on constrained mixture theory. Our goal was to apply these two theories on the same experimental data. Within the experiment, a pulmonary artery was exposed to systemic conditions. The change in diameter was followed-up over time. A mechanical and microstructural analysis of native pulmonary artery and pulmonary autograft was conducted. Whereas the kinematic-based growth theory is able to accurately capture the growth of the tissue, it does not account for the mechanobiological processes causing this growth. The constrained mixture theory takes into account the mechanobiological processes including removal, deposition and adaptation of all structural constituents, allowing us to simulate a changing microstructure and mechanical behavior.

Sex Differences in Drug-Induced Arrhythmogenesis.

The electrical activity in the heart varies significantly between men and women and results in a sex-specific response to drugs. Recent evidence suggests that women are more than twice as likely as men to develop drug-induced arrhythmia with potentially fatal consequences. Yet, the sex-specific differences in drug-induced arrhythmogenesis remain poorly understood. Here we integrate multiscale modeling and machine learning to gain mechanistic insight into the sex-specific origin of drug-induced cardiac arrhythmia at differing drug concentrations. To quantify critical drug concentrations in male and female hearts, we identify the most important ion channels that trigger male and female arrhythmogenesis, and create and train a sex-specific multi-fidelity arrhythmogenic risk classifier. Our study reveals that sex differences in ion channel activity, tissue conductivity, and heart dimensions trigger longer QT-intervals in women than in men. We quantify the critical drug concentration for dofetilide, a high risk drug, to be seven times lower for women than for men. Our results emphasize the importance of including sex as an independent biological variable in risk assessment during drug development. Acknowledging and understanding sex differences in drug safety evaluation is critical when developing novel therapeutic treatments on a personalized basis. The general trends of this study have significant implications on the development of safe and efficacious new drugs and the prescription of existing drugs in combination with other drugs.

Bayesian Physics-Based Modeling of Tau Propagation in Alzheimer's Disease.

Amyloid-ß and hyperphosphorylated tau protein are known drivers of neuropathology in Alzheimer's disease. Tau in particular spreads in the brains of patients following a spatiotemporal pattern that is highly sterotypical and correlated with subsequent neurodegeneration. Novel medical imaging techniques can now visualize the distribution of tau in the brain in vivo, allowing for new insights to the dynamics of this biomarker. Here we personalize a network diffusion model with global spreading and local production terms to longitudinal tau positron emission tomography data of 76 subjects from the Alzheimer's Disease Neuroimaging Initiative. We use Bayesian inference with a hierarchical prior structure to infer means and credible intervals for our model parameters on group and subject levels. Our results show that the group average protein production rate for amyloid positive subjects is significantly higher with 0.019±0.27/yr, than that for amyloid negative subjects with -0.143±0.21/yr (p = 0.0075). These results support the hypothesis that amyloid pathology drives tau pathology. The calibrated model could serve as a valuable clinical tool to identify optimal time points for follow-up scans and predict the timeline of disease progression.

Effects of B.1.1.7 and B.1.351 on COVID-19 Dynamics: A Campus Reopening Study.

The timing and sequence of safe campus reopening has remained the most controversial topic in higher education since the outbreak of the COVID-19 pandemic. By the end of March 2020, almost all colleges and universities in the United States had transitioned to an all online education and many institutions have not yet fully reopened to date. For a residential campus like Stanford University, the major challenge of reopening is to estimate the number of incoming infectious students at the first day of class. Here we learn the number of incoming infectious students using Bayesian inference and perform a series of retrospective and projective simulations to quantify the risk of campus reopening. We create a physics-based probabilistic model to infer the local reproduction dynamics for each state and adopt a network SEIR model to simulate the return of all undergraduates, broken down by their year of enrollment and state of origin. From these returning student populations, we predict the outbreak dynamics throughout the spring, summer, fall, and winter quarters using the inferred reproduction dynamics of Santa Clara County. We compare three different scenarios: the true outbreak dynamics under the wild-type SARS-CoV-2, and the hypothetical outbreak dynamics under the new COVID-19 variants B.1.1.7 and B.1.351 with 56% and 50% increased transmissibility. Our study reveals that even small changes in transmissibility can have an enormous impact on the overall case numbers. With no additional countermeasures, during the most affected quarter, the fall of 2020, there would have been 203 cases under baseline reproduction, compared to 4727 and 4256 cases for the B.1.1.7 and B.1.351 variants. Our results suggest that population mixing presents an increased risk for local outbreaks, especially with new and more infectious variants emerging across the globe. Tight outbreak control through mandatory quarantine and test-trace-isolate strategies will be critical in successfully managing these local outbreak dynamics.

Visualizing the invisible: The effect of asymptomatic transmission on the outbreak dynamics of COVID-19.

Understanding the outbreak dynamics of the COVID-19 pandemic has important implications for successful containment and mitigation strategies. Recent studies suggest that the population prevalence of SARS-CoV-2 antibodies, a proxy for the number of asymptomatic cases, could be an order of magnitude larger than expected from the number of reported symptomatic cases. Knowing the precise prevalence and contagiousness of asymptomatic transmission is critical to estimate the overall dimension and pandemic potential of COVID-19. However, at this stage, the effect of the asymptomatic population, its size, and its outbreak dynamics remain largely unknown. Here we use reported symptomatic case data in conjunction with antibody seroprevalence studies, a mathematical epidemiology model, and a Bayesian framework to infer the epidemiological characteristics of COVID-19. Our model computes, in real time, the time-varying contact rate of the outbreak, and projects the temporal evolution and credible intervals of the effective reproduction number and the symptomatic, asymptomatic, and recovered populations. Our study quantifies the sensitivity of the outbreak dynamics of COVID-19 to three parameters: the effective reproduction number, the ratio between the symptomatic and asymptomatic populations, and the infectious periods of both groups For nine distinct locations, our model estimates the fraction of the population that has been infected and recovered by Jun 15, 2020 to 24.15% (95% CI: 20.48%-28.14%) for Heinsberg (NRW, Germany), 2.40% (95% CI: 2.09%-2.76%) for Ada County (ID, USA), 46.19% (95% CI: 45.81%-46.60%) for New York City (NY, USA), 11.26% (95% CI: 7.21%-16.03%) for Santa Clara County (CA, USA), 3.09% (95% CI: 2.27%-4.03%) for Denmark, 12.35% (95% CI: 10.03%-15.18%) for Geneva Canton (Switzerland), 5.24% (95% CI: 4.84%-5.70%) for the Netherlands, 1.53% (95% CI: 0.76%-2.62%) for Rio Grande do Sul (Brazil), and 5.32% (95% CI: 4.77%-5.93%) for Belgium. Our method traces the initial outbreak date in Santa Clara County back to January 20, 2020 (95% CI: December 29, 2019 - February 13, 2020). Our results could significantly change our understanding and management of the COVID-19 pandemic: A large asymptomatic population will make isolation, containment, and tracing of individual cases challenging. Instead, managing community transmission through increasing population awareness, promoting physical distancing, and encouraging behavioral changes could become more relevant.

The reproduction number of COVID-19 and its correlation with public health interventions.

Throughout the past six months, no number has dominated the public media more persistently than the reproduction number of COVID-19. This powerful but simple concept is widely used by the public media, scientists, and political decision makers to explain and justify political strategies to control the COVID-19 pandemic. Here we explore the effectiveness of political interventions using the reproduction number of COVID-19 across Europe. We propose a dynamic SEIR epidemiology model with a time-varying reproduction number, which we identify using machine learning. During the early outbreak, the basic reproduction number was 4.22 ± 1.69, with maximum values of 6.33 and 5.88 in Germany and the Netherlands. By May 10, 2020, it dropped to 0.67 ± 0.18, with minimum values of 0.37 and 0.28 in Hungary and Slovakia. We found a strong correlation between passenger air travel, driving, walking, and transit mobility and the effective reproduction number with a time delay of 17.24 ± 2.00 days. Our new dynamic SEIR model provides the flexibility to simulate various outbreak control and exit strategies to inform political decision making and identify safe solutions in the benefit of global health.

The reproduction number of COVID-19 and its correlation with public health interventions.

Throughout the past six months, no number has dominated the public media more persistently than the reproduction number of COVID-19. This powerful but simple concept is widely used by the public media, scientists, and political decision makers to explain and justify political strategies to control the COVID-19 pandemic. Here we explore the effectiveness of political interventions using the reproduction number of COVID-19 across Europe. We propose a dynamic SEIR epidemiology model with a time-varying reproduction number, which we identify using machine learning. During the early outbreak, the basic reproduction number was 4.22+/-1.69, with maximum values of 6.33 and 5.88 in Germany and the Netherlands. By May 10, 2020, it dropped to 0.67+/-0.18, with minimum values of 0.37 and 0.28 in Hungary and Slovakia. We found a strong correlation between passenger air travel, driving, walking, and transit mobility and the effective reproduction number with a time delay of 17.24+/-2.00 days. Our new dynamic SEIR model provides the flexibility to simulate various outbreak control and exit strategies to inform political decision making and identify safe solutions in the benefit of global health.

Outbreak dynamics of COVID-19 in Europe and the effect of travel restrictions.

For the first time in history, on March 17, 2020, the European Union closed all its external borders in an attempt to contain the spreading of the coronavirus 2019, COVID-19. Throughout two past months, governments around the world have implemented massive travel restrictions and border control to mitigate the outbreak of this global pandemic. However, the precise effects of travel restrictions on the outbreak dynamics of COVID-19 remain unknown. Here we combine a global network mobility model with a local epidemiology model to simulate and predict the outbreak dynamics and outbreak control of COVID-19 across Europe. We correlate our mobility model to passenger air travel statistics and calibrate our epidemiology model using the number of reported COVID-19 cases for each country. Our simulations show that mobility networks of air travel can predict the emerging global diffusion pattern of a pandemic at the early stages of the outbreak. Our results suggest that an unconstrained mobility would have significantly accelerated the spreading of COVID-19, especially in Central Europe, Spain, and France. Ultimately, our network epidemiology model can inform political decision making and help identify exit strategies from current travel restrictions and total lockdown.

Outbreak dynamics of COVID-19 in China and the United States.

On March 11, 2020, the World Health Organization declared the coronavirus disease 2019, COVID-19, a global pandemic. In an unprecedented collective effort, massive amounts of data are now being collected worldwide to estimate the immediate and long-term impact of this pandemic on the health system and the global economy. However, the precise timeline of the disease, its transmissibility, and the effect of mitigation strategies remain incompletely understood. Here we integrate a global network model with a local epidemic SEIR model to quantify the outbreak dynamics of COVID-19 in China and the United States. For the outbreak in China, in [Formula: see text] provinces, we found a latent period of 2.56 ± 0.72 days, a contact period of 1.47 ± 0.32 days, and an infectious period of 17.82 ± 2.95 days. We postulate that the latent and infectious periods are disease-specific, whereas the contact period is behavior-specific and can vary between different provinces, states, or countries. For the early stages of the outbreak in the United States, in [Formula: see text] states, we adopted the disease-specific values from China and found a contact period of 3.38 ± 0.69 days. Our network model predicts that-without the massive political mitigation strategies that are in place today-the United States would have faced a basic reproduction number of 5.30 ± 0.95 and a nationwide peak of the outbreak on May 10, 2020 with 3 million infections. Our results demonstrate how mathematical modeling can help estimate outbreak dynamics and provide decision guidelines for successful outbreak control. We anticipate that our model will become a valuable tool to estimate the potential of vaccination and quantify the effect of relaxing political measures including total lockdown, shelter in place, and travel restrictions for low-risk subgroups of the population or for the population as a whole.

Visualizing the invisible: The effect of asymptomatic transmission on the outbreak dynamics of COVID-19.

Understanding the outbreak dynamics of the COVID-19 pandemic has important implications for successful containment and mitigation strategies. Recent studies suggest that the population prevalence of SARS-CoV-2 antibodies, a proxy for the number of asymptomatic cases, could be an order of magnitude larger than expected from the number of reported symptomatic cases. Knowing the precise prevalence and contagiousness of asymptomatic transmission is critical to estimate the overall dimension and pandemic potential of COVID-19. However, at this stage, the effect of the asymptomatic population, its size, and its outbreak dynamics remain largely unknown. Here we use reported symptomatic case data in conjunction with antibody seroprevalence studies, a mathematical epidemiology model, and a Bayesian framework to infer the epidemiological characteristics of COVID-19. Our model computes, in real time, the time-varying contact rate of the outbreak, and projects the temporal evolution and credible intervals of the effective reproduction number and the symptomatic, asymptomatic, and recovered populations. Our study quantifies the sensitivity of the outbreak dynamics of COVID-19 to three parameters: the effective reproduction number, the ratio between the symptomatic and asymptomatic populations, and the infectious periods of both groups. For nine distinct locations, our model estimates the fraction of the population that has been infected and recovered by Jun 15, 2020 to 24.15% (95% CI: 20.48%-28.14%) for Heinsberg (NRW, Germany), 2.40% (95% CI: 2.09%-2.76%) for Ada County (ID, USA), 46.19% (95% CI: 45.81%-46.60%) for New York City (NY, USA), 11.26% (95% CI: 7.21%-16.03%) for Santa Clara County (CA, USA), 3.09% (95% CI: 2.27%-4.03%) for Denmark, 12.35% (95% CI: 10.03%-15.18%) for Geneva Canton (Switzerland), 5.24% (95% CI: 4.84%-5.70%) for the Netherlands, 1.53% (95% CI: 0.76%-2.62%) for Rio Grande do Sul (Brazil), and 5.32% (95% CI: 4.77%-5.93%) for Belgium. Our method traces the initial outbreak date in Santa Clara County back to January 20, 2020 (95% CI: December 29, 2019-February 13, 2020). Our results could significantly change our understanding and management of the COVID-19 pandemic: A large asymptomatic population will make isolation, containment, and tracing of individual cases challenging. Instead, managing community transmission through increasing population awareness, promoting physical distancing, and encouraging behavioral changes could become more relevant.

Kinematic boundary conditions substantially impact in silico ventricular function.

Computational cardiac mechanical models, individualized to the patient, have the potential to elucidate the fundamentals of cardiac (patho-)physiology, enable non-invasive quantification of clinically significant metrics (eg, stiffness, active contraction, work), and anticipate the potential efficacy of therapeutic cardiovascular intervention. In a clinical setting, however, the available imaging resolution is often limited, which limits cardiac models to focus on the ventricles, without including the atria, valves, and proximal arteries and veins. In such models, the absence of surrounding structures needs to be accounted for by imposing realistic kinematic boundary conditions, which, for prognostic purposes, are preferably generic and thus non-image derived. Unfortunately, the literature on cardiac models shows no consistent approach to kinematically constrain the myocardium. The impact of different approaches (eg, fully constrained base, constrained epi-ring) on the predictive capacity of cardiac mechanical models has not been thoroughly studied. For that reason, this study first gives an overview of current approaches to kinematically constrain (bi) ventricular models. Next, we developed a patient-specific in silico biventricular model that compares well with literature and in vivo recorded strains. Alternative constraints were introduced to assess the influence of commonly used mechanical boundary conditions on both the predicted global functional behavior of the in-silico heart (cavity volumes, stroke volume, ejection fraction) and local strain distributions. Meaningful differences in global functioning were found between different kinematic anchoring strategies, which brought forward the importance of selecting appropriate boundary conditions for biventricular models that, in the near future, may inform clinical intervention. However, whilst statistically significant differences were also found in local strain distributions, these differences were minor and mostly confined to the region close to the applied boundary conditions.

A modular inverse elastostatics approach to resolve the pressure-induced stress state for in vivo imaging based cardiovascular modeling.

Patient-specific biomechanical modeling of the cardiovascular system is complicated by the presence of a physiological pressure load given that the imaged tissue is in a pre-stressed and -strained state. Neglect of this prestressed state into solid tissue mechanics models leads to erroneous metrics (e.g. wall deformation, peak stress, wall shear stress) which in their turn are used for device design choices, risk assessment (e.g. procedure, rupture) and surgery planning. It is thus of utmost importance to incorporate this deformed and loaded tissue state into the computational models, which implies solving an inverse problem (calculating an undeformed geometry given the load and the deformed geometry). Methodologies to solve this inverse problem can be categorized into iterative and direct methodologies, both having their inherent advantages and disadvantages. Direct methodologies are typically based on the inverse elastostatics (IE) approach and offer a computationally efficient single shot methodology to compute the in vivo stress state. However, cumbersome and problem-specific derivations of the formulations and non-trivial access to the finite element analysis (FEA) code, especially for commercial products, refrain a broad implementation of these methodologies. For that reason, we developed a novel, modular IE approach and implemented this methodology in a commercial FEA solver with minor user subroutine interventions. The accuracy of this methodology was demonstrated in an arterial tube and porcine biventricular myocardium model. The computational power and efficiency of the methodology was shown by computing the in vivo stress and strain state, and the corresponding unloaded geometry, for two models containing multiple interacting incompressible, anisotropic (fiber-embedded) and hyperelastic material behaviors: a patient-specific abdominal aortic aneurysm and a full 4-chamber heart model.

An in silico framework to analyze the anisotropic shear wave mechanics in cardiac shear wave elastography.

Shear wave elastography (SWE) is a potential tool to non-invasively assess cardiac muscle stiffness. This study focused on the effect of the orthotropic material properties and mechanical loading on the performance of cardiac SWE, as it is known that these factors contribute to complex 3D anisotropic shear wave propagation. To investigate the specific impact of these complexities, we constructed a finite element model with an orthotropic material law subjected to different uniaxial stretches to simulate SWE in the stressed cardiac wall. Group and phase speed were analyzed in function of tissue thickness and virtual probe rotation angle. Tissue stretching increased the group and phase speed of the simulated shear wave, especially in the direction of the muscle fiber. As the model provided access to the true fiber orientation and material properties, we assessed the accuracy of two fiber orientation extraction methods based on SWE. We found a higher accuracy (but lower robustness) when extracting fiber orientations based on the location of maximal shear wave speed instead of the angle of the major axis of the ellipsoidal group speed surface. Both methods had a comparable performance for the center region of the cardiac wall, and performed less well towards the edges. Lastly, we also assessed the (theoretical) impact of pathology on shear wave physics and characterization in the model. It was found that SWE was able to detect changes in fiber orientation and material characteristics, potentially associated with cardiac pathologies such as myocardial fibrosis. Furthermore, the model showed clearly altered shear wave patterns for the fibrotic myocardium compared to the healthy myocardium, which forms an initial but promising outcome of this modeling study.

Erratum to: An in silico biomechanical analysis of the stent-esophagus interaction.

In the original publication of the article, Tables 2 and 3 were published with error. The correct tables are provided below (Tables 2, 3). The original version of the article has also been corrected.