Non-invasive cardiovascular magnetic resonance assessment of pressure recovery distance after aortic valve stenosis.

BACKGROUND: Decisions in the management of aortic stenosis are based on the peak pressure drop, captured by Doppler echocardiography, whereas gold standard catheterization measurements assess the net pressure drop but are limited by associated risks. The relationship between these two measurements, peak and net pressure drop, is dictated by the pressure recovery along the ascending aorta which is mainly caused by turbulence energy dissipation. Currently, pressure recovery is considered to occur within the first 40-50 mm distally from the aortic valve, albeit there is inconsistency across interventionist centers on where/how to position the catheter to capture the net pressure drop. METHODS: We developed a non-invasive method to assess the pressure recovery distance based on blood flow momentum via 4D Flow cardiovascular magnetic resonance (CMR). Multi-center acquisitions included physical flow phantoms with different stenotic valve configurations to validate this method, first against reference measurements and then against turbulent energy dissipation (respectively n = 8 and n = 28 acquisitions) and to investigate the relationship between peak and net pressure drops. Finally, we explored the potential errors of cardiac catheterisation pressure recordings as a result of neglecting the pressure recovery distance in a clinical bicuspid aortic valve (BAV) cohort of n = 32 patients. RESULTS: In-vitro assessment of pressure recovery distance based on flow momentum achieved an average error of 1.8 ± 8.4 mm when compared to reference pressure sensors in the first phantom workbench. The momentum pressure recovery distance and the turbulent energy dissipation distance showed no statistical difference (mean difference of 2.8 ± 5.4 mm, R2 = 0.93) in the second phantom workbench. A linear correlation was observed between peak and net pressure drops, however, with strong dependences on the valvular morphology. Finally, in the BAV cohort the pressure recovery distance was 78.8 ± 34.3 mm from vena contracta, which is significantly longer than currently accepted in clinical practise (40-50 mm), and 37.5% of patients displayed a pressure recovery distance beyond the end of the ascending aorta. CONCLUSION: The non-invasive assessment of the distance to pressure recovery is possible by tracking momentum via 4D Flow CMR. Recovery is not always complete at the ascending aorta, and catheterised recordings will overestimate the net pressure drop in those situations. There is a need to re-evaluate the methods that characterise the haemodynamic burden caused by aortic stenosis as currently clinically accepted pressure recovery distance is an underestimation.

CRIMSON: An open-source software framework for cardiovascular integrated modelling and simulation.

In this work, we describe the CRIMSON (CardiovasculaR Integrated Modelling and SimulatiON) software environment. CRIMSON provides a powerful, customizable and user-friendly system for performing three-dimensional and reduced-order computational haemodynamics studies via a pipeline which involves: 1) segmenting vascular structures from medical images; 2) constructing analytic arterial and venous geometric models; 3) performing finite element mesh generation; 4) designing, and 5) applying boundary conditions; 6) running incompressible Navier-Stokes simulations of blood flow with fluid-structure interaction capabilities; and 7) post-processing and visualizing the results, including velocity, pressure and wall shear stress fields. A key aim of CRIMSON is to create a software environment that makes powerful computational haemodynamics tools accessible to a wide audience, including clinicians and students, both within our research laboratories and throughout the community. The overall philosophy is to leverage best-in-class open source standards for medical image processing, parallel flow computation, geometric solid modelling, data assimilation, and mesh generation. It is actively used by researchers in Europe, North and South America, Asia, and Australia. It has been applied to numerous clinical problems; we illustrate applications of CRIMSON to real-world problems using examples ranging from pre-operative surgical planning to medical device design optimization.

Noninvasive quantification of cerebrovascular pressure changes using 4D Flow MRI.

PURPOSE: Hemodynamic alterations are indicative of cerebrovascular disease. However, the narrow and tortuous cerebrovasculature complicates image-based assessment, especially when quantifying relative pressure. Here, we present a systematic evaluation of image-based cerebrovascular relative pressure mapping, investigating the accuracy of the routinely used reduced Bernoulli (RB), the extended unsteady Bernoulli (UB), and the full-field virtual work-energy relative pressure ( ν WERP) method. METHODS: Patient-specific in silico models were used to generate synthetic cerebrovascular 4D Flow MRI, with RB, UB, and ν WERP performance quantified as a function of spatiotemporal sampling and image noise. Cerebrovascular relative pressures were also derived in 4D Flow MRI from healthy volunteers ( n=8 ), acquired at two spatial resolutions (dx = 1.1 and 0.8 mm). RESULTS: The in silico analysis indicate that accurate relative pressure estimations are inherently coupled to spatial sampling: at dx = 1.0 mm high errors are reported for all methods; at dx = 0.5 mm ν WERP recovers relative pressures at a mean error of 0.02 ± 0.25 mm Hg, while errors remain higher for RB and UB (mean error of -2.18 ± 1.91 and -2.18 ± 1.87 mm Hg, respectively). The dependence on spatial sampling is also indicated in vivo, albeit with higher correlative dependence between resolutions using ν WERP (k = 0.64, R2 = 0.81 for dx = 1.1 vs. 0.8 mm) than with RB or UB (k = 0.04, R2 = 0.03, and k = 0.07, R2 = 0.07, respectively). CONCLUSION: Image-based full-field methods such as ν WERP enable cerebrovascular relative pressure mapping; however, accuracy is directly dependent on utilized spatial resolution.

False lumen pressure estimation in type B aortic dissection using 4D flow cardiovascular magnetic resonance: comparisons with aortic growth.

BACKGROUND: Chronic type B aortic dissection (TBAD) is associated with poor long-term outcome, and accurate risk stratification tools remain lacking. Pressurization of the false lumen (FL) has been recognized as central in promoting aortic growth. Several surrogate imaging-based metrics have been proposed to assess FL hemodynamics; however, their relationship to enlarging aortic dimensions remains unclear. We investigated the association between aortic growth and three cardiovascular magnetic resonance (CMR)-derived metrics of FL pressurization: false lumen ejection fraction (FLEF), maximum systolic deceleration rate (MSDR), and FL relative pressure (FL ΔPmax). METHODS: CMR/CMR angiography was performed in 12 patients with chronic dissection of the descending thoracoabdominal aorta, including contrast-enhanced CMR angiography and time-resolved three-dimensional phase-contrast CMR (4D Flow CMR). Aortic growth rate was calculated as the change in maximal aortic diameter between baseline and follow-up imaging studies over the time interval, with patients categorized as having either 'stable' (< 3 mm/year) or 'enlarging' (≥ 3 mm/year) growth. Three metrics relating to FL pressurization were defined as: (1) FLEF: the ratio between retrograde and antegrade flow at the TBAD entry tear, (2) MSDR: the absolute difference between maximum and minimum systolic acceleration in the proximal FL, and (3) FL ΔPmax: the difference in absolute pressure between aortic root and distal FL. RESULTS: FLEF was higher in enlarging TBAD (49.0 ± 17.9% vs. 10.0 ± 11.9%, p = 0.002), whereas FL ΔPmax was lower (32.2 ± 10.8 vs. 57.2 ± 12.5 mmHg/m, p = 0.017). MSDR and conventional anatomic variables did not differ significantly between groups. FLEF showed positive (r = 0.78, p = 0.003) correlation with aortic growth rate whereas FL ΔPmax showed negative correlation (r = - 0.64, p = 0.026). FLEF and FL ΔPmax remained as independent predictors of aortic growth rate after adjusting for baseline aortic diameter. CONCLUSION: Comparative analysis of three 4D flow CMR metrics of TBAD FL pressurization demonstrated that those that focusing on retrograde flow (FLEF) and relative pressure (FL ΔPmax) independently correlated with growth and differentiated patients with enlarging and stable descending aortic dissections. These results emphasize the highly variable nature of aortic hemodynamics in TBAD patients, and suggest that 4D Flow CMR derived metrics of FL pressurization may be useful to separate patients at highest and lowest risk for progressive aortic growth and complications.

A partition of unity approach to fluid mechanics and fluid-structure interaction.

For problems involving large deformations of thin structures, simulating fluid-structure interaction (FSI) remains a computationally expensive endeavour which continues to drive interest in the development of novel approaches. Overlapping domain techniques have been introduced as a way to combine the fluid-solid mesh conformity, seen in moving-mesh methods, without the need for mesh smoothing or re-meshing, which is a core characteristic of fixed mesh approaches. In this work, we introduce a novel overlapping domain method based on a partition of unity approach. Unified function spaces are defined as a weighted sum of fields given on two overlapping meshes. The method is shown to achieve optimal convergence rates and to be stable for steady-state Stokes, Navier-Stokes, and ALE Navier-Stokes problems. Finally, we present results for FSI in the case of 2D flow past an elastic beam simulation. These initial results point to the potential applicability of the method to a wide range of FSI applications, enabling boundary layer refinement and large deformations without the need for re-meshing or user-defined stabilization.

A class of analytic solutions for verification and convergence analysis of linear and nonlinear fluid-structure interaction algorithms.

Fluid-structure interaction (FSI) problems are pervasive in the computational engineering community. The need to address challenging FSI problems has led to the development of a broad range of numerical methods addressing a variety of applicationspecific demands. While a range of numerical and experimental benchmarks are present in the literature, few solutions are available that enable both verification and spatiotemporal convergence analysis. In this paper, we introduce a class of analytic solutions to FSI problems involving shear in channels and pipes. Comprised of 16 separate analytic solutions, our approach is permuted to enable progressive verification and analysis of FSI methods and implementations, in two and three dimensions, for static and transient scenarios as well as for linear and hyperelastic solid materials. Results are shown for a range of analytic models exhibiting progressively complex behavior. The utility of these solutions for analysis of convergence behavior is further demonstrated using a previously published monolithic FSI technique. The resulting class of analytic solutions addresses a core challenge in the development of novel FSI algorithms and implementations, providing a progressive testbed for verification and detailed convergence analysis.

An efficient and accurate method for modeling nonlinear fractional viscoelastic biomaterials.

Computational biomechanics plays an important role in biomedical engineering: using modeling to understand pathophysiology, treatment and device design. While experimental evidence indicates that the mechanical response of most tissues is viscoelastic, current biomechanical models in the computational community often assume hyperelastic material models. Fractional viscoelastic constitutive models have been successfully used in literature to capture viscoelastic material response; however, the translation of these models into computational platforms remains limited. Many experimentally derived viscoelastic constitutive models are not suitable for three-dimensional simulations. Furthermore, the use of fractional derivatives can be computationally prohibitive, with a number of current numerical approximations having a computational cost that is 𝒪 ( N T 2 ) and a storage cost that is 𝒪(NT ) (NT denotes the number of time steps). In this paper, we present a novel numerical approximation to the Caputo derivative which exploits a recurrence relation similar to those used to discretize classic temporal derivatives, giving a computational cost that is 𝒪(NT ) and a storage cost that is fixed over time. The approximation is optimized for numerical applications, and an error estimate is presented to demonstrate the efficacy of the method. The method, integrated into a finite element solid mechanics framework, is shown to be unconditionally stable in the linear viscoelastic case. It was then integrated into a computational biomechanical framework, with several numerical examples verifying the accuracy and computational efficiency of the method, including in an analytic test, in an analytic fractional differential equation, as well as in a computational biomechanical model problem.

Stiffness reconstruction methods for MR elastography.

Assessment of tissue stiffness is desirable for clinicians and researchers, as it is well established that pathophysiological mechanisms often alter the structural properties of tissue. Magnetic resonance elastography (MRE) provides an avenue for measuring tissue stiffness and has a long history of clinical application, including staging liver fibrosis and stratifying breast cancer malignancy. A vital component of MRE consists of the reconstruction algorithms used to derive stiffness from wave-motion images by solving inverse problems. A large range of reconstruction methods have been presented in the literature, with differing computational expense, required user input, underlying physical assumptions, and techniques for numerical evaluation. These differences, in turn, have led to varying accuracy, robustness, and ease of use. While most reconstruction techniques have been validated against in silico or in vitro phantoms, performance with real data is often more challenging, stressing the robustness and assumptions of these algorithms. This article reviews many current MRE reconstruction methods and discusses the aforementioned differences. The material assumptions underlying the methods are developed and various approaches for noise reduction, regularization, and numerical discretization are discussed. Reconstruction methods are categorized by inversion type, underlying assumptions, and their use in human and animal studies. Future directions, such as alternative material assumptions, are also discussed.

Bridging Three Orders of Magnitude: Multiple Scattered Waves Sense Fractal Microscopic Structures via Dispersion.

Wave scattering provides profound insight into the structure of matter. Typically, the ability to sense microstructure is determined by the ratio of scatterer size to probing wavelength. Here, we address the question of whether macroscopic waves can report back the presence and distribution of microscopic scatterers despite several orders of magnitude difference in scale between wavelength and scatterer size. In our analysis, monosized hard scatterers 5 µm in radius are immersed in lossless gelatin phantoms to investigate the effect of multiple reflections on the propagation of shear waves with millimeter wavelength. Steady-state monochromatic waves are imaged in situ via magnetic resonance imaging, enabling quantification of the phase velocity at a voxel size big enough to contain thousands of individual scatterers, but small enough to resolve the wavelength. We show in theory, experiments, and simulations that the resulting coherent superposition of multiple reflections gives rise to power-law dispersion at the macroscopic scale if the scatterer distribution exhibits apparent fractality over an effective length scale that is comparable to the probing wavelength. Since apparent fractality is naturally present in any random medium, microstructure can thereby leave its fingerprint on the macroscopically quantifiable power-law exponent. Our results are generic to wave phenomena and carry great potential for sensing microstructure that exhibits intrinsic fractality, such as, for instance, vasculature.

Aortic relative pressure components derived from four-dimensional flow cardiovascular magnetic resonance.

PURPOSE: To describe the assessment of the spatiotemporal distribution of relative aortic pressure quantifying the magnitude of its three major components. METHODS: Nine healthy volunteers and three patients with aortic disease (bicuspid aortic valve, dissection, and Marfan syndrome) underwent 4D-flow CMR. Spatiotemporal pressure maps were computed from the CMR flow fields solving the pressure Poisson equation. The individual components of pressure were separated into time-varying inertial ("transient"), spatially varying inertial ("convective"), and viscous components. RESULTS: Relative aortic pressure is primarily caused by transient effects followed by the convective and small viscous contributions (64.5, 13.6, and 0.3 mmHg/m, respectively, in healthy subjects), although regional analysis revealed prevalent convective effects in specific contexts, e.g., Sinus of Valsalva and aortic arch at instants of peak velocity. Patients showed differences in peak transient values and duration, and localized abrupt convective changes explained by abnormalities in aortic geometry, including the presence of an aneurysm, a pseudo-coarctation, the inlet of a dissection, or by complex flow patterns. CONCLUSION: The evaluation of the three components of relative pressure enables the quantification of mechanistic information for understanding and stratifying aortic disease, with potential future implications for guiding therapy.

A displacement-based finite element formulation for incompressible and nearly-incompressible cardiac mechanics.

The Lagrange Multiplier (LM) and penalty methods are commonly used to enforce incompressibility and compressibility in models of cardiac mechanics. In this paper we show how both formulations may be equivalently thought of as a weakly penalized system derived from the statically condensed Perturbed Lagrangian formulation, which may be directly discretized maintaining the simplicity of penalty formulations with the convergence characteristics of LM techniques. A modified Shamanskii-Newton-Raphson scheme is introduced to enhance the nonlinear convergence of the weakly penalized system and, exploiting its equivalence, modifications are developed for the penalty form. Focusing on accuracy, we proceed to study the convergence behavior of these approaches using different interpolation schemes for both a simple test problem and more complex models of cardiac mechanics. Our results illustrate the well-known influence of locking phenomena on the penalty approach (particularly for lower order schemes) and its effect on accuracy for whole-cycle mechanics. Additionally, we verify that direct discretization of the weakly penalized form produces similar convergence behavior to mixed formulations while avoiding the use of an additional variable. Combining a simple structure which allows the solution of computationally challenging problems with good convergence characteristics, the weakly penalized form provides an accurate and efficient alternative to incompressibility and compressibility in cardiac mechanics.

Generalized super-resolution 4D Flow MRI - using ensemble learning to extend across the cardiovascular system.

4D Flow Magnetic Resonance Imaging (4D Flow MRI) is a non-invasive measurement technique capable of quantifying blood flow across the cardiovascular system. While practical use is limited by spatial resolution and image noise, incorporation of trained super-resolution (SR) networks has potential to enhance image quality post-scan. However, these efforts have predominantly been restricted to narrowly defined cardiovascular domains, with limited exploration of how SR performance extends across the cardiovascular system; a task aggravated by contrasting hemodynamic conditions apparent across the cardiovasculature. The aim of our study was therefore to explore the generalizability of SR 4D Flow MRI using a combination of existing super-resolution base models, novel heterogeneous training sets, and dedicated ensemble learning techniques; the latter-most being effectively used for improved domain adaption in other domains or modalities, however, with no previous exploration in the setting of 4D Flow MRI. With synthetic training data generated across three disparate domains (cardiac, aortic, cerebrovascular), varying convolutional base and ensemble learners were evaluated as a function of domain and architecture, quantifying performance on both in-silico and acquired in-vivo data from the same three domains. Results show that both bagging and stacking ensembling enhance SR performance across domains, accurately predicting high-resolution velocities from low-resolution input data in-silico. Likewise, optimized networks successfully recover native resolution velocities from downsampled in-vivo data, as well as show qualitative potential in generating denoised SR-images from clinicallevel input data. In conclusion, our work presents a viable approach for generalized SR 4D Flow MRI, with the novel use of ensemble learning in the setting of advanced fullfield flow imaging extending utility across various clinical areas of interest.

MicroBundleCompute: Automated segmentation, tracking, and analysis of subdomain deformation in cardiac microbundles.

Advancing human induced pluripotent stem cell derived cardiomyocyte (hiPSC-CM) technology will lead to significant progress ranging from disease modeling, to drug discovery, to regenerative tissue engineering. Yet, alongside these potential opportunities comes a critical challenge: attaining mature hiPSC-CM tissues. At present, there are multiple techniques to promote maturity of hiPSC-CMs including physical platforms and cell culture protocols. However, when it comes to making quantitative comparisons of functional behavior, there are limited options for reliably and reproducibly computing functional metrics that are suitable for direct cross-system comparison. In addition, the current standard functional metrics obtained from time-lapse images of cardiac microbundle contraction reported in the field (i.e., post forces, average tissue stress) do not take full advantage of the available information present in these data (i.e., full-field tissue displacements and strains). Thus, we present "MicroBundleCompute," a computational framework for automatic quantification of morphology-based mechanical metrics from movies of cardiac microbundles. Briefly, this computational framework offers tools for automatic tissue segmentation, tracking, and analysis of brightfield and phase contrast movies of beating cardiac microbundles. It is straightforward to implement, runs without user intervention, requires minimal input parameter setting selection, and is computationally inexpensive. In this paper, we describe the methods underlying this computational framework, show the results of our extensive validation studies, and demonstrate the utility of exploring heterogeneous tissue deformations and strains as functional metrics. With this manuscript, we disseminate "MicroBundleCompute" as an open-source computational tool with the aim of making automated quantitative analysis of beating cardiac microbundles more accessible to the community.

In silico modeling of shear-stress-induced nitric oxide production in endothelial cells through systems biology.

Nitric oxide (NO) produced by vascular endothelial cells is a potent vasodilator and an antiinflammatory mediator. Regulating production of endothelial-derived NO is a complex undertaking, involving multiple signaling and genetic pathways that are activated by diverse humoral and biomechanical stimuli. To gain a thorough understanding of the rich diversity of responses observed experimentally, it is necessary to account for an ensemble of these pathways acting simultaneously. In this article, we have assembled four quantitative molecular pathways previously proposed for shear-stress-induced NO production. In these pathways, endothelial NO synthase is activated 1), via calcium release, 2), via phosphorylation reactions, and 3), via enhanced protein expression. To these activation pathways, we have added a fourth, a pathway describing actual NO production from endothelial NO synthase and its various protein partners. These pathways were combined and simulated using CytoSolve, a computational environment for combining independent pathway calculations. The integrated model is able to describe the experimentally observed change in NO production with time after the application of fluid shear stress. This model can also be used to predict the specific effects on the system after interventional pharmacological or genetic changes. Importantly, this model reflects the up-to-date understanding of the NO system, providing a platform upon which information can be aggregated in an additive way.

A viscoelastic constitutive model for human femoropopliteal arteries.

High failure rates present challenges for surgical and interventional therapies for peripheral artery disease of the femoropopliteal artery (FPA). The FPA's demanding biomechanical environment necessitates complex interactions with repair devices and materials. While a comprehensive understanding of the FPA's mechanical characteristics could improve medical treatments, the viscoelastic properties of these muscular arteries remain poorly understood, and the constitutive model describing their time-dependent behavior is absent. We introduce a new viscoelastic constitutive model for the human FPA grounded in its microstructural composition. The model is capable of detailing the contributions of each intramural component to the overall viscoelastic response. Our model was developed utilizing fractional viscoelasticity and tested using biaxial experimental data with hysteresis and relaxation collected from 10 healthy human subjects aged 57 to 65 and further optimized for high throughput and automation. The model accurately described the experimental data, capturing significant nonlinearity and hysteresis that were particularly pronounced circumferentially, and tracked the contribution of passive smooth muscle cells to viscoelasticity that was twice that of the collagen fibers. The high-throughput parameter estimation procedure we developed included a specialized objective function and modifications to enhance convergence for the common exponential-type fiber laws, facilitating computational implementation. Our new model delineates the time-dependent behavior of human FPAs, which will improve the fidelity of computational simulations investigating device-artery interactions and contribute to their greater physical accuracy. Moreover, it serves as a useful tool to investigate the contribution of arterial constituents to overall tissue viscoelasticity, thereby expanding our knowledge of arterial mechanophysiology. STATEMENT OF SIGNIFICANCE: The demanding biomechanical environment of the femoropopliteal artery (FPA) necessitates complex interactions with repair devices and materials, but the viscoelastic properties of these muscular arteries remain poorly understood with the constitutive model describing their time-dependent behavior being absent. We hereby introduce the first viscoelastic constitutive model for the human FPA grounded in its microstructures. This model was tested using biaxial mechanical data collected from 10 healthy human subjects between the ages of 57 to 65. It can detail the contributions of each intramural component to the overall viscoelastic response, showing that the contribution of passive smooth muscle cells to viscoelasticity is twice that of collagen fibers. The usefulness of this model as tool to better understand arterial mechanophysiology was demonstrated.

Evaluation of aortic stenosis: From Bernoulli and Doppler to Navier-Stokes.

Uni-dimensional Doppler echocardiography data provide the mainstay of quantative assessment of aortic stenosis, with the transvalvular pressure drop a key indicator of haemodynamic burden. Sophisticated methods of obtaining velocity data, combined with improved computational analysis, are facilitating increasingly robust and reproducible measurement. Imaging modalities which permit acquisition of three-dimensional blood velocity vector fields enable angle-independent valve interrogation and calculation of enhanced measures of the transvalvular pressure drop. This manuscript clarifies the fundamental principles of physics that underpin the evaluation of aortic stenosis and explores modern techniques that may provide more accurate means to grade aortic stenosis and inform appropriate management.

A data-driven computational model for engineered cardiac microtissues.

Engineered heart tissues (EHTs) present a potential solution to some of the current challenges in the treatment of heart disease; however, the development of mature, adult-like cardiac tissues remains elusive. Mechanical stimuli have been observed to improve whole-tissue function and cardiomyocyte (CM) maturation, although our ability to fully utilize these mechanisms is hampered, in part, by our incomplete understanding of the mechanobiology of EHTs. In this work, we leverage experimental data, produced by a mechanically tunable experimental setup, to introduce a tissue-specific computational modeling pipeline of EHTs. Our new modeling pipeline generates simulated, image-based EHTs, capturing ECM and myofibrillar structure as well as functional parameters estimated directly from experimental data. This approach enables the unique estimation of EHT function by data-based estimation of CM active stresses. We use this experimental and modeling pipeline to study different mechanical environments, where we contrast the force output of the tissue with the computed active stress of CMs. We show that the significant differences in measured experimental forces can largely be explained by the levels of myofibril formation achieved by the CMs in the distinct mechanical environments, with active stress showing more muted variations across conditions. The presented model also enables us to dissect the relative contributions of myofibrils and extracellular matrix to tissue force output, a task difficult to address experimentally. These results highlight the importance of tissue-specific modeling to augment EHT experiments, providing deeper insights into the mechanobiology driving EHT function. STATEMENT OF SIGNIFICANCE: Engineered heart tissues (EHTs) have the potential to revolutionize the way heart disease is treated. However, developing mature cardiomyocytes (CM) in these tissues remains a challenge due, in part, to our incomplete understanding of the fundamental biomechanical mechanisms that drive EHT development. This work integrates the experimental data of an EHT platform developed to study the influence of mechanics in CM maturation with computational biomechanical models. This approach is used to augment conclusions obtained in-vitro - by measuring quantities such as cell stress and strain - and to dissect the relevance of each component in the whole tissue performance. Our results show how a combination of specialized in-silico and in-vitro approaches can help us better understand the mechanobiology of EHTs.

lGeneralized super-resolution 4D Flow MRI-using ensemble learning to extend across the cardiovascular system.

4D Flow Magnetic Resonance Imaging (4D Flow MRI) is a non-invasive measurement technique capable of quantifying blood flow across the cardiovascular system. While practical use is limited by spatial resolution and image noise, incorporation of trained super-resolution (SR) networks has potential to enhance image quality post-scan. However, these efforts have predominantly been restricted to narrowly defined cardiovascular domains, with limited exploration of how SR performance extends across the cardiovascular system; a task aggravated by contrasting hemodynamic conditions apparent across the cardiovasculature. The aim of our study was to explore the generalizability of SR 4D Flow MRI using a combination of heterogeneous training sets and dedicated ensemble learning. With synthetic training data generated across three disparate domains (cardiac, aortic, cerebrovascular), varying convolutional base and ensemble learners were evaluated as a function of domain and architecture, quantifying performance on both in-silico and acquired in-vivo data from the same three domains. Results show that both bagging and stacking ensembling enhance SR performance across domains, accurately predicting high-resolution velocities from low-resolution input data in-silico. Likewise, optimized networks successfully recover native resolution velocities from downsampled in-vivo data, as well as show qualitative potential in generating denoised SR-images from clinicallevel input data. In conclusion, our work presents a viable approach for generalized SR 4D Flow MRI, with ensemble learning extending utility across various clinical areas of interest.

Matrix architecture and mechanics regulate myofibril organization, costamere assembly, and contractility of engineered myocardial microtissues.

The mechanical function of the myocardium is defined by cardiomyocyte contractility and the biomechanics of the extracellular matrix (ECM). Understanding this relationship remains an important unmet challenge due to limitations in existing approaches for engineering myocardial tissue. Here, we established arrays of cardiac microtissues with tunable mechanics and architecture by integrating ECM-mimetic synthetic, fiber matrices and induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs), enabling real-time contractility readouts, in-depth structural assessment, and tissue-specific computational modeling. We find that the stiffness and alignment of matrix fibers distinctly affect the structural development and contractile function of pure iPSC-CM tissues. Further examination into the impact of fibrous matrix stiffness enabled by computational models and quantitative immunofluorescence implicates cell-ECM interactions in myofibril assembly and notably costamere assembly, which correlates with improved contractile function of tissues. These results highlight how iPSC-CM tissue models with controllable architecture and mechanics can inform the design of translatable regenerative cardiac therapies.

Personalization of biomechanical simulations of the left ventricle by in-vivo cardiac DTI data: Impact of fiber interpolation methods.

Simulations of cardiac electrophysiology and mechanics have been reported to be sensitive to the microstructural anisotropy of the myocardium. Consequently, a personalized representation of cardiac microstructure is a crucial component of accurate, personalized cardiac biomechanical models. In-vivo cardiac Diffusion Tensor Imaging (cDTI) is a non-invasive magnetic resonance imaging technique capable of probing the heart's microstructure. Being a rather novel technique, issues such as low resolution, signal-to noise ratio, and spatial coverage are currently limiting factors. We outline four interpolation techniques with varying degrees of data fidelity, different amounts of smoothing strength, and varying representation error to bridge the gap between the sparse in-vivo data and the model, requiring a 3D representation of microstructure across the myocardium. We provide a workflow to incorporate in-vivo myofiber orientation into a left ventricular model and demonstrate that personalized modelling based on fiber orientations from in-vivo cDTI data is feasible. The interpolation error is correlated with a trend in personalized parameters and simulated physiological parameters, strains, and ventricular twist. This trend in simulation results is consistent across material parameter settings and therefore corresponds to a bias introduced by the interpolation method. This study suggests that using a tensor interpolation approach to personalize microstructure with in-vivo cDTI data, reduces the fiber uncertainty and thereby the bias in the simulation results.