Heart Size Difference Drives Sex-Specific Response to Cardiac Resynchronization Therapy: A Post Hoc Analysis of the MORE-MPP CRT Trial.

BACKGROUND: Studies have reported that female sex predicts superior cardiac resynchronization therapy (CRT) response. One theory is that this association is related to smaller female heart size, thus increased relative dyssynchrony at a given QRS duration (QRSd). Our objective was to investigate the mechanisms of sex-specific CRT response relating to heart size, relative dyssynchrony, cardiomyopathy type, QRS morphology, and other patient characteristics. METHODS AND RESULTS: This is a post hoc analysis of the MORE-CRT MPP (More Response on Cardiac Resynchronization Therapy with Multipoint Pacing)  trial (n=3739, 28% women), with a subgroup analysis of patients with nonischemic cardiomyopathy and left bundle-branch block (n=1308, 41% women) to control for confounding characteristics. A multivariable analysis examined predictors of response to 6 months of conventional CRT, including sex and relative dyssynchrony, measured by QRSd/left ventricular end-diastolic volume (LVEDV). Women had a higher CRT response rate than men (70.1% versus 56.8%, P<0.0001). In subgroup analysis, regression analysis of the nonischemic cardiomyopathy left bundle-branch block subgroup identified QRSd/LVEDV, but not sex, as a modifier of CRT response (P<0.0039). QRSd/LVEDV was significantly higher in women (0.919) versus men (0.708, P<0.001). CRT response was 78% for female patients with QRSd/LVEDV greater than the median value, compared with 68% with QRSd/LVEDV less than the median value (P=0.012). The association between CRT response and QRSd/LVEDV was strongest at QRSd <150 ms. CONCLUSIONS: In the nonischemic cardiomyopathy left bundle-branch block population, increased relative dyssynchrony in women, who have smaller heart sizes than their male counterparts, is a driver of sex-specific CRT response, particularly at QRSd <150 ms. Women may benefit from CRT at a QRSd <130 ms, opening the debate on whether sex-specific QRSd cutoffs or QRS/LVEDV measurement should be incorporated into clinical guidelines.

Observable Atrial and Ventricular Fibrillation Episode Durations Are Conformant With a Power Law Based on System Size and Spatial Synchronization.

BACKGROUND: Atrial fibrillation (AF) and ventricular fibrillation (VF) episodes exhibit varying durations, with some spontaneously ending quickly while others persist. A quantitative framework to explain episode durations remains elusive. We hypothesized that observable self-terminating AF and VF episode lengths, whereby durations are known, would conform with a power law based on the ratio of system size and correlation length ([Formula: see text]. METHODS: Using data from computer simulations (2-dimensional sheet and 3-dimensional left-atrial), human ischemic VF recordings (256-electrode sock, n=12 patients), and human AF recordings (64-electrode basket-catheter, n=9 patients; 16-electrode HD-grid catheter, n=42 patients), conformance with a power law was assessed using the Akaike information criterion, Bayesian information criterion, coefficient of determination (R2, significance=P<0.05) and maximum likelihood estimation. We analyzed fibrillatory episode durations and [Formula: see text], computed by taking the ratio between system size ([Formula: see text], chamber/simulation size) and correlation length ([Formula: see text], measured from pairwise correlation coefficients over electrode/node distance). RESULTS: In all computer models, the relationship between episode durations and [Formula: see text] was conformant with a power law (Aliev-Panfilov R2: 0.90, P<0.001; Courtemanche R2: 0.91, P<0.001; Luo-Rudy R2: 0.61, P<0.001). Observable clinical AF/VF durations were also conformant with a power law relationship (VF R2: 0.86, P<0.001; AF basket R2: 0.91, P<0.001; AF grid R2: 0.92, P<0.001). [Formula: see text] also differentiated between self-terminating and sustained episodes of AF and VF (P<0.001; all systems), as well as paroxysmal versus persistent AF (P<0.001). In comparison, other electrogram metrics showed no statistically significant differences (dominant frequency, Shannon Entropy, mean voltage, peak-peak voltage; P>0.05). CONCLUSIONS: Observable fibrillation episode durations are conformant with a power law based on system size and correlation length.

Evaluation of Deep Learning Estimation of Whole Heart Anatomy from Automated Cardiovascular Magnetic Resonance Short- and Long-Axis Analyses in UK Biobank.

AIMS: Standard methods of heart chamber volume estimation in cardiovascular magnetic resonance (CMR) typically utilize simple geometric formulae based on a limited number of slices. We aimed to evaluate whether an automated deep learning neural network prediction of 3D anatomy of all four chambers would show stronger associations with cardiovascular risk factors and disease than standard volume estimation methods in the UK Biobank. METHODS: A deep learning network was adapted to predict 3D segmentations of left and right ventricles (LV, RV) and atria (LA, RA) at ∼1mm isotropic resolution from CMR short and long axis 2D segmentations obtained from a fully automated machine learning pipeline in 4723 individuals with cardiovascular disease (CVD) and 5733 without in the UK Biobank. Relationships between volumes at end-diastole (ED) and end-systole (ES) and risk/disease factors were quantified using univariate, multivariate and logistic regression analyses. Strength of association between deep learning volumes and standard volumes was compared using the area under the receiving operator characteristic curve (AUC). RESULTS: Univariate and multivariate associations between deep learning volumes and most risk and disease factors were stronger than for standard volumes (higher R2 and more significant P values), particularly for sex, age, and body mass index. AUC for all logistic regressions were higher for deep learning volumes than standard volumes (p<0.001 for all four chambers at ED and ES). CONCLUSIONS: Neural network reconstructions of whole heart volumes had significantly stronger associations with cardiovascular disease and risk factors than standard volume estimation methods in an automatic processing pipeline.

Cardiac length-dependent activation driven by force-dependent thick-filament dynamics.

The length-dependent activation (LDA) of maximum force and calcium sensitivity are established features of cardiac muscle contraction but the dominant underlying mechanisms remain to be fully clarified. Alongside the well-documented regulation of contraction via the thin filaments, experiments have identified an additional force-dependent thick-filament activation, whereby myosin heads parked in a so-called off state become available to generate force. This process produces a feedback effect that may potentially drive LDA. Using biomechanical modeling of a human left-ventricular myocyte, this study investigates the extent to which the off-state dynamics could, by itself, plausibly account for LDA, depending on the specific mathematical formulation of the feedback. We hypothesized four different models of the off-state regulatory feedback based on (A) total force, (B) active force, (C) sarcomere strain, and (D) passive force. We tested if these models could reproduce the isometric steady-state and dynamic LDA features predicted by an earlier published model of a human left-ventricle myocyte featuring purely phenomenological length dependences. The results suggest that only total-force feedback (A) is capable of reproducing the expected behaviors, but that passive tension could provide a length-dependent signal on which to initiate the feedback. Furthermore, by attributing LDA to off-state dynamics, our proposed model also qualitatively reproduces experimentally observed effects of the off-state-stabilizing drug mavacamten. Taken together, these results support off-state dynamics as a plausible primary mechanism underlying LDA.

The Virtual Child.

SUMMARY: We are building the world's first Virtual Child-a computer model of normal and cancerous human development at the level of each individual cell. The Virtual Child will "develop cancer" that we will subject to unlimited virtual clinical trials that pinpoint, predict, and prioritize potential new treatments, bringing forward the day when no child dies of cancer, giving each one the opportunity to lead a full and healthy life.

Results of the 2023 ISBI challenge to reduce GABA-edited MRS acquisition time.

PURPOSE: Use a conference challenge format to compare machine learning-based gamma-aminobutyric acid (GABA)-edited magnetic resonance spectroscopy (MRS) reconstruction models using one-quarter of the transients typically acquired during a complete scan. METHODS: There were three tracks: Track 1: simulated data, Track 2: identical acquisition parameters with in vivo data, and Track 3: different acquisition parameters with in vivo data. The mean squared error, signal-to-noise ratio, linewidth, and a proposed shape score metric were used to quantify model performance. Challenge organizers provided open access to a baseline model, simulated noise-free data, guides for adding synthetic noise, and in vivo data. RESULTS: Three submissions were compared. A covariance matrix convolutional neural network model was most successful for Track 1. A vision transformer model operating on a spectrogram data representation was most successful for Tracks 2 and 3. Deep learning (DL) reconstructions with 80 transients achieved equivalent or better SNR, linewidth and fit error compared to conventional 320 transient reconstructions. However, some DL models optimized linewidth and SNR without actually improving overall spectral quality, indicating a need for more robust metrics. CONCLUSION: DL-based reconstruction pipelines have the promise to reduce the number of transients required for GABA-edited MRS.

Neural network emulation of the human ventricular cardiomyocyte action potential for more efficient computations in pharmacological studies.

Computer models of the human ventricular cardiomyocyte action potential (AP) have reached a level of detail and maturity that has led to an increasing number of applications in the pharmaceutical sector. However, interfacing the models with experimental data can become a significant computational burden. To mitigate the computational burden, the present study introduces a neural network (NN) that emulates the AP for given maximum conductances of selected ion channels, pumps, and exchangers. Its applicability in pharmacological studies was tested on synthetic and experimental data. The NN emulator potentially enables massive speed-ups compared to regular simulations and the forward problem (find drugged AP for pharmacological parameters defined as scaling factors of control maximum conductances) on synthetic data could be solved with average root-mean-square errors (RMSE) of 0.47 mV in normal APs and of 14.5 mV in abnormal APs exhibiting early afterdepolarizations (72.5% of the emulated APs were alining with the abnormality, and the substantial majority of the remaining APs demonstrated pronounced proximity). This demonstrates not only very fast and mostly very accurate AP emulations but also the capability of accounting for discontinuities, a major advantage over existing emulation strategies. Furthermore, the inverse problem (find pharmacological parameters for control and drugged APs through optimization) on synthetic data could be solved with high accuracy shown by a maximum RMSE of 0.22 in the estimated pharmacological parameters. However, notable mismatches were observed between pharmacological parameters estimated from experimental data and distributions obtained from the Comprehensive in vitro Proarrhythmia Assay initiative. This reveals larger inaccuracies which can be attributed particularly to the fact that small tissue preparations were studied while the emulator was trained on single cardiomyocyte data. Overall, our study highlights the potential of NN emulators as powerful tool for an increased efficiency in future quantitative systems pharmacology studies.

Whole-heart electromechanical simulations using Latent Neural Ordinary Differential Equations.

Cardiac digital twins provide a physics and physiology informed framework to deliver personalized medicine. However, high-fidelity multi-scale cardiac models remain a barrier to adoption due to their extensive computational costs. Artificial Intelligence-based methods can make the creation of fast and accurate whole-heart digital twins feasible. We use Latent Neural Ordinary Differential Equations (LNODEs) to learn the pressure-volume dynamics of a heart failure patient. Our surrogate model is trained from 400 simulations while accounting for 43 parameters describing cell-to-organ cardiac electromechanics and cardiovascular hemodynamics. LNODEs provide a compact representation of the 3D-0D model in a latent space by means of an Artificial Neural Network that retains only 3 hidden layers with 13 neurons per layer and allows for numerical simulations of cardiac function on a single processor. We employ LNODEs to perform global sensitivity analysis and parameter estimation with uncertainty quantification in 3 hours of computations, still on a single processor.

A three-dimensional left atrial motion estimation from retrospective gated computed tomography: application in heart failure patients with atrial fibrillation.

Background: A reduced left atrial (LA) strain correlates with the presence of atrial fibrillation (AF). Conventional atrial strain analysis uses two-dimensional (2D) imaging, which is, however, limited by atrial foreshortening and an underestimation of through-plane motion. Retrospective gated computed tomography (RGCT) produces high-fidelity three-dimensional (3D) images of the cardiac anatomy throughout the cardiac cycle that can be used for estimating 3D mechanics. Its feasibility for LA strain measurement, however, is understudied. Aim: The aim of this study is to develop and apply a novel workflow to estimate 3D LA motion and calculate the strain from RGCT imaging. The utility of global and regional strains to separate heart failure in patients with reduced ejection fraction (HFrEF) with and without AF is investigated. Methods: A cohort of 30 HFrEF patients with (n = 9) and without (n = 21) AF underwent RGCT prior to cardiac resynchronisation therapy. The temporal sparse free form deformation image registration method was optimised for LA feature tracking in RGCT images and used to estimate 3D LA endocardial motion. The area and fibre reservoir strains were calculated over the LA body. Universal atrial coordinates and a human atrial fibre atlas enabled the regional strain calculation and the fibre strain calculation along the local myofibre orientation, respectively. Results: It was found that global reservoir strains were significantly reduced in the HFrEF + AF group patients compared with the HFrEF-only group patients (area strain: 11.2 ± 4.8% vs. 25.3 ± 12.6%, P = 0.001; fibre strain: 4.5 ± 2.0% vs. 15.2 ± 8.8%, P = 0.001), with HFrEF + AF patients having a greater regional reservoir strain dyssynchrony. All regional reservoir strains were reduced in the HFrEF + AF patient group, in whom the inferior wall strains exhibited the most significant differences. The global reservoir fibre strain and LA volume + posterior wall reservoir fibre strain exceeded LA volume alone and 2D global longitudinal strain (GLS) for AF classification (area-under-the-curve: global reservoir fibre strain: 0.94 ± 0.02, LA volume + posterior wall reservoir fibre strain: 0.95 ± 0.02, LA volume: 0.89 ± 0.03, 2D GLS: 0.90 ± 0.03). Conclusion: RGCT enables 3D LA motion estimation and strain calculation that outperforms 2D strain metrics and LA enlargement for AF classification. Differences in regional LA strain could reflect regional myocardial properties such as atrial fibrosis burden.

Performance of Image-navigated and Diaphragm-navigated 3D Late Gadolinium-enhanced Cardiac MRI for the Assessment of Atrial Fibrosis.

Purpose To perform a qualitative and quantitative evaluation of the novel image-navigated (iNAV) 3D late gadolinium enhancement (LGE) cardiac MRI imaging strategy in comparison with the conventional diaphragm-navigated (dNAV) 3D LGE cardiac MRI strategy for the assessment of left atrial fibrosis in atrial fibrillation (AF). Materials and Methods In this prospective study conducted between April and September 2022, 26 consecutive participants with AF (mean age, 61 ± 11 years; 19 male) underwent both iNAV and dNAV 3D LGE cardiac MRI, with equivalent spatial resolution and timing in the cardiac cycle. Participants were randomized in the acquisition order of iNAV and dNAV. Both, iNAV-LGE and dNAV-LGE images were analyzed qualitatively using a 5-point Likert scale and quantitatively (percentage of atrial fibrosis using image intensity ratio threshold 1.2), including testing for overlap in atrial fibrosis areas by calculating Dice score. Results Acquisition time of iNAV was significantly lower compared with dNAV (4.9 ± 1.1 minutes versus 12 ± 4 minutes, P < .001, respectively). There was no evidence of a difference in image quality for all prespecified criteria between iNAV and dNAV, although dNAV was the preferred image strategy in two-thirds of cases (17/26, 65%). Quantitative assessment demonstrated that mean fibrosis scores were lower for iNAV compared with dNAV (12 ± 8% versus 20 ± 12%, P < .001). Spatial correspondence between the atrial fibrosis maps was modest (Dice similarity coefficient, 0.43 ± 0.15). Conclusion iNAV-LGE acquisition in individuals with AF was more than twice as fast as dNAV acquisition but resulted in a lower atrial fibrosis score. The differences between these two strategies might impact clinical interpretation. ©RSNA, 2024.

Enhancing transvenous lead extraction risk prediction: Integrating imaging biomarkers into machine learning models.

BACKGROUND: Machine learning (ML) models have been proposed to predict risk related to transvenous lead extraction (TLE). OBJECTIVE: The purpose of this study was to test whether integrating imaging data into an existing ML model increases its ability to predict major adverse events (MAEs; procedure-related major complications and procedure-related deaths) and lengthy procedures (≥100 minutes). METHODS: We hypothesized certain features-(1) lead angulation, (2) coil percentage inside the superior vena cava (SVC), and (3) number of overlapping leads in the SVC-detected from a pre-TLE plain anteroposterior chest radiograph (CXR) would improve prediction of MAE and long procedural times. A deep-learning convolutional neural network was developed to automatically detect these CXR features. RESULTS: A total of 1050 cases were included, with 24 MAEs (2.3%) . The neural network was able to detect (1) heart border with 100% accuracy; (2) coils with 98% accuracy; and (3) acute angle in the right ventricle and SVC with 91% and 70% accuracy, respectively. The following features significantly improved MAE prediction: (1) ≥50% coil within the SVC; (2) ≥2 overlapping leads in the SVC; and (3) acute lead angulation. Balanced accuracy (0.74-0.87), sensitivity (68%-83%), specificity (72%-91%), and area under the curve (AUC) (0.767-0.962) all improved with imaging biomarkers. Prediction of lengthy procedures also improved: balanced accuracy (0.76-0.86), sensitivity (75%-85%), specificity (63%-87%), and AUC (0.684-0.913). CONCLUSION: Risk prediction tools integrating imaging biomarkers significantly increases the ability of ML models to predict risk of MAE and long procedural time related to TLE.

High Prevalence of New Clinically Significant Findings in Patients With Embolic Stroke of Unknown Source Evaluated by Cardiac Magnetic Resonance Imaging.

BACKGROUND: Embolic stroke of unknown source (ESUS) accounts for 1 in 6 ischemic strokes. Current guidelines do not recommend routine cardiac magnetic resonance (CMR) imaging in ESUS, and beyond the identification of cardioembolic sources, there are no data assessing new clinical findings from CMR in ESUS. This study aimed to assess the prevalence of new cardiac and noncardiac findings and to determine their impact on clinical care in patients with ESUS. METHODS AND RESULTS: In this prospective, multicenter, observational study, CMR imaging was performed within 3 months of ESUS. All scans were reported according to standard clinical practice. A new clinical finding was defined as one not previously identified through prior clinical evaluation. A clinically significant finding was defined as one resulting in further investigation, follow-up, or treatment. A change in patient care was defined as initiation of medical, interventional, surgical, or palliative care. From 102 patients recruited, 96 underwent CMR imaging. One or more new clinical findings were observed in 59 patients (61%). New findings were clinically significant in 48 (81%) of these patients. Of 40 patients with a new clinically significant cardiac finding, 21 (53%) experienced a change in care (medical therapy, n=15; interventional/surgical procedure, n=6). In 12 patients with a new clinically significant extracardiac finding, 6 (50%) experienced a change in care (medical therapy, n=4; palliative care, n=2). CONCLUSIONS: CMR imaging identifies new clinically significant cardiac and noncardiac findings in half of patients with recent ESUS. Advanced cardiovascular screening should be considered in patients with ESUS. REGISTRATION: URL: https://www.clinicaltrials.gov; Unique identifier: NCT04555538.

Atrial fibrosis by cardiac MRI is a correlate for atrial stiffness in patients with atrial fibrillation.

A relationship between left atrial strain and pressure has been demonstrated in many studies, but not in an atrial fibrillation (AF) cohort. In this work, we hypothesized that elevated left atrial (LA) tissue fibrosis might mediate and confound the LA strain vs. pressure relationship, resulting instead in a relationship between LA fibrosis and stiffness index (mean LA pressure/LA reservoir strain). Sixty-seven patients with AF underwent a standard cardiac MR exam including long-axis cine views (2 and 4-ch) and a free-breathing high resolution three-dimensional late gadolinium enhancement (LGE) of the atrium (N = 41), within 30 days prior to AF ablation, at which procedure invasive mean left atrial pressure (LAP) was measured. LV and LA Volumes, EF, and comprehensive analysis of LA strains (strain and strain rates and strain timings during the atrial reservoir, conduit and active, i.e. active atrial contraction, phases) were measured and LA fibrosis content (LGE (ml)) was assessed from 3D LGE volumes. LA LGE was well correlated to atrial stiffness index overall (R = 0.59, p < 0.001), and among patient subgroups. Pressure was only correlated to maximal LA volume (R = 0.32) and the time to peak reservoir strain rate (R = 0.32) (both p < 0.01), among all functional measurements. LA reservoir strain was strongly correlated with LAEF (R = 0.95, p < 0.001) and LA minimum volume (r = 0.82, p < 0.001). In our AF cohort, pressure is correlated to maximum LA volume and time to peak reservoir strain. LA pressure/ LA reservoir strain, a metric of stiffness, correlates with LA fibrosis (LA LGE), reflecting Hook's Law.

Developing Cardiac Digital Twins at Scale: Insights from Personalised Myocardial Conduction Velocity.

Large-cohort studies using cardiovascular imaging and diagnostic datasets have assessed cardiac anatomy, function, and outcomes, but typically do not reveal underlying biological mechanisms. Cardiac digital twins (CDTs) provide personalized physics- and physiology-constrained in-silico representations, enabling inference of multi-scale properties tied to these mechanisms. We constructed 3464 anatomically-accurate CDTs using cardiac magnetic resonance images from UK biobank and personalised their myocardial conduction velocities (CVs) from electrocardiograms (ECG), through an automated framework. We found well-known sex-specific differences in QRS duration were fully explained by myocardial anatomy, as CV remained consistent across sexes. Conversely, significant associations of CV with ageing and increased BMI suggest myocardial tissue remodelling. Novel associations were observed with left ventricular ejection fraction and mental-health phenotypes, through a phenome-wide association study, and CV was also linked with adverse clinical outcomes. Our study highlights the utility of population-based CDTs in assessing intersubject variability and uncovering strong links with mental health.

Two-Stage Self-Supervised Contrastive Learning Aided Transformer for Real-Time Medical Image Segmentation.

The availability of large, high-quality annotated datasets in the medical domain poses a substantial challenge in segmentation tasks. To mitigate the reliance on annotated training data, self-supervised pre-training strategies have emerged, particularly employing contrastive learning methods on dense pixel-level representations. In this work, we proposed to capitalize on intrinsic anatomical similarities within medical image data and develop a semantic segmentation framework through a self-supervised fusion network, where the availability of annotated volumes is limited. In a unified training phase, we combine segmentation loss with contrastive loss, enhancing the distinction between significant anatomical regions that adhere to the available annotations. To further improve the segmentation performance, we introduce an efficient parallel transformer module that leverages Multiview multiscale feature fusion and depth-wise features. The proposed transformer architecture, based on multiple encoders, is trained in a self-supervised manner using contrastive loss. Initially, the transformer is trained using an unlabeled dataset. We then fine-tune one encoder using data from the first stage and another encoder using a small set of annotated segmentation masks. These encoder features are subsequently concatenated for the purpose of brain tumor segmentation. The multiencoder-based transformer model yields significantly better outcomes across three medical image segmentation tasks. We validated our proposed solution by fusing images across diverse medical image segmentation challenge datasets, demonstrating its efficacy by outperforming state-of-the-art methodologies.

Using Gaussian process for velocity reconstruction after coronary stenosis applicable in positron emission particle tracking: An in-silico study.

Accurate velocity reconstruction is essential for assessing coronary artery disease. We propose a Gaussian process method to reconstruct the velocity profile using the sparse data of the positron emission particle tracking (PEPT) in a biological environment, which allows the measurement of tracer particle velocity to infer fluid velocity fields. We investigated the influence of tracer particle quantity and detection time interval on flow reconstruction accuracy. Three models were used to represent different levels of stenosis and anatomical complexity: a narrowed straight tube, an idealized coronary bifurcation with stenosis, and patient-specific coronary arteries with a stenotic left circumflex artery. Computational fluid dynamics (CFD), particle tracking, and the Gaussian process of kriging were employed to simulate and reconstruct the pulsatile flow field. The study examined the error and uncertainty in velocity profile reconstruction after stenosis by comparing particle-derived flow velocity with the CFD solution. Using 600 particles (15 batches of 40 particles) released in the main coronary artery, the time-averaged error in velocity reconstruction ranged from 13.4% (no occlusion) to 161% (70% occlusion) in patient-specific anatomy. The error in maximum cross-sectional velocity at peak flow was consistently below 10% in all cases. PEPT and kriging tended to overestimate area-averaged velocity in higher occlusion cases but accurately predicted maximum cross-sectional velocity, particularly at peak flow. Kriging was shown to be useful to estimate the maximum velocity after the stenosis in the absence of negative near-wall velocity.

Constructing bilayer and volumetric atrial models at scale.

To enable large in silico trials and personalized model predictions on clinical timescales, it is imperative that models can be constructed quickly and reproducibly. First, we aimed to overcome the challenges of constructing cardiac models at scale through developing a robust, open-source pipeline for bilayer and volumetric atrial models. Second, we aimed to investigate the effects of fibres, fibrosis and model representation on fibrillatory dynamics. To construct bilayer and volumetric models, we extended our previously developed coordinate system to incorporate transmurality, atrial regions and fibres (rule-based or data driven diffusion tensor magnetic resonance imaging (MRI)). We created a cohort of 1000 biatrial bilayer and volumetric models derived from computed tomography (CT) data, as well as models from MRI, and electroanatomical mapping. Fibrillatory dynamics diverged between bilayer and volumetric simulations across the CT cohort (correlation coefficient for phase singularity maps: left atrial (LA) 0.27 ± 0.19, right atrial (RA) 0.41 ± 0.14). Adding fibrotic remodelling stabilized re-entries and reduced the impact of model type (LA: 0.52 ± 0.20, RA: 0.36 ± 0.18). The choice of fibre field has a small effect on paced activation data (less than 12 ms), but a larger effect on fibrillatory dynamics. Overall, we developed an open-source user-friendly pipeline for generating atrial models from imaging or electroanatomical mapping data enabling in silico clinical trials at scale (https://github.com/pcmlab/atrialmtk).

Heart size disparity drives sex-specific response to cardiac resynchronization therapy: a post-hoc analysis of the MORE-MPP CRT trial.

Background: Studies have reported that female sex predicts superior cardiac resynchronization therapy (CRT) response. One theory is that this association is related to smaller female heart size, thus increased "relative dyssynchrony" at given QRS durations (QRSd). Objective: To investigate the mechanisms of sex-specific CRT response relating to heart size, relative dyssynchrony, cardiomyopathy type, QRS morphology, and other patient characteristics. Methods: A post-hoc analysis of the MORE-CRT MPP trial (n=3739, 28% female), with a sub-group analysis of patients with non-ischaemic cardiomyopathy (NICM) and left bundle branch block (LBBB) (n=1308, 41% female) to control for confounding characteristics. A multivariable analysis examined predictors of response to 6 months of conventional CRT, including sex and relative dyssynchrony, measured by QRSd/LVEDV (left ventricular end-diastolic volume). Results: Females had a higher CRT response rate than males (70.1% vs. 56.8%, p<0.0001). Subgroup analysis: Regression analysis of the NICM LBBB subgroup identified QRSd/LVEDV, but not sex, as a modifier of CRT response (p<0.0039). QRSd/LVEDV was significantly higher in females (0.919) versus males (0.708, p<0.001). CRT response was 78% for female patients with QRSd/LVEDV>median value, compared to 68% < median value (p=0.012). Association between CRT response and QRSd/LVEDV was strongest at QRSd<150ms. Conclusions: In the NICM LBBB population, increased relative dyssynchrony in females, who have smaller heart sizes than their male counterparts, is a driver of sex-specific CRT response, particularly at QRSd <150ms. Females may benefit from CRT at a QRSd <130ms, opening the debate on whether sex-specific QRSd cut-offs or QRS/LVEDV measurement should be incorporated into clinical guidelines.