Personalized Computational Modeling Identifies Embolic Stroke of Undetermined Source Patients with Potential Arrhythmic Substrate

Patient Characteristics

Ninety patient-derived models were included in our analysis: 45 post-stroke ESUS and 45 pre-ablation AFib patients. Demographic information about both patient groups is provided in Table 1. There was a significant difference in LA surface area; however, the trend towards higher body mass index (BMI) in AFib patients and the lack of difference in LA volume index between the two groups suggests the latter trend (i.e., higher surface area in AFib) would be lessened if the values were normalized to body surface area. LA fibrosis burden was not significantly different between ESUS (13.6 ± 6.2%) and AFib patients (14.2 ± 4.5%) (P=0.91), consistent with previous findings (4).

Patient-derived model in silico AFib induction and fibrosis quantification

Personalized LA bilayer models were generated for all ESUS and AFib patients. Examples of physiological detail incorporated in models can be seen in Fig. 1, including patient-specific patterns of fibrotic remodeling, realistic atrial fiber orientations, and locations of electric pacing sites; further detail can be found in the Methods section. Rapid electric stimulation caused RD-sustained arrhythmia in 28 of 45 AFib models (62.2%) and 23 of 45 ESUS models (51.1%). Thus, the capability of the fibrotic substrate to sustain RDs was not significantly different between the two groups (P=0.39, Fig. 2A). ESUS and AFib patients were then sorted by amount of global LA fibrosis and arranged into quartiles. For five patients (21.7%) in the first quartile (fibrosis<9.75%), six patients (29.0%) in the second quartile (fibrosis<12.6%), nineteen (79.2%) patients in the third quartile (fibrosis<17.45%), and twenty-one (95.5%) patients in the top quartile, simulations in the corresponding LA model revealed at least one pacing site for which stimulation produced an episode of RD-sustained arrhythmia (Fig. 2B).

• Download figure
• Open in new tab

Figure 1

Model generation (A) Reconstruction of LA geometry with anatomical features labelled (RIPV/RSPV/LIPV/LSPV, right/left inferior/superior pulmonary veins; LAA, LA appendage). The LA is modeled as a bilayer comprising nested endocardial and epicardial shells linked in both fibrotic and non-fibrotic regions by 1D linear elements. (B) LA fiber orientations for the endocardium and epicardium, mapped from human atlas geometry as described in Methods. (C) AFib trigger sites as pacing sites (posterior/anterior LIPV, LSPV, RSPV, RIPV, LAA base, mitral valve annulus and posterior wall). (D) Regions of the LA generated as described in methods: (1) atrial floor, (2) anterior wall and LAA, (3) posterior wall, (4) left PVs, (5) right PVs.

• Download figure
• Open in new tab

Figure 2

Summary of patient-derived model fibrosis with respect to RD inducibility. (A) Histogram of AFib (28/45) and ESUS (23/45) inducible patients. Inducibility was not significantly different by χ² test. (B) Patients with ESUS and AFib arranged by percentage of LA fibrosis. Dotted lines indicate the quartiles of fibrosis observed for all 90 patients. Circles are indicative of stable reentry observed in the model from at least one pacing site after in silico pacing protocol. Triangles indicate no RDs after pacing from all 15 pacing sites independently. Cases that lacked RDs despite high fibrosis (inset I) or were inducible despite low fibrosis (inset II) are highlighted.

To explore potential pro-arrhythmic substrate properties in ESUS and AFib, we analyzed fibrosis burden in the sub-groups of each cohort in which RD-sustained arrhythmias were inducible and non-inducible (Fig. 3A). Fibrosis burden was not significantly different between ESUS and AFib with (P=0.068, confidence interval; CI: [-0.2, 6.1], Fig. 3A) or without induced reentry (P=0.58, CI: [-2.4, 1.5], Fig. 3A). However, when fibrosis burdens for inducible and non-inducible models were aggregated across ESUS and AFib groups, a significant difference was evident (Fig. 3A; 16.8 ± 5.04% vs. 10.19 ± 3.14%; P<0.0001, CI: [4.4, 8.3]).

• Download figure
• Open in new tab

Figure 3

Summary of RD characteristics between ESUS and AFib models. (A) Boxplot of fibrosis percentage in ESUS and AFib models where reentry was induced (ESUS: N=23, IQR=8.1; AFib: N=28, IQR=5.95) and where reentry was not induced (ESUS: N=22, IQR=4.98; AFib N=17, IQR=3.5). Across ESUS and AFib models, fibrosis burden for RD inducible and RD non-inducible models was significantly different (P<0.0001). (B) Boxplot of number of unique reentrant morphologies elicited by all 15 pacing sites (P=0.79). (C) Correlation plot of fibrosis vs number of RDs (R=0.69, P<0.0001) (D) Boxplot of the number of pacing sites which induced reentry (P=0.98) (E) Correlation plot of fibrosis vs number of pacing sites that induced reentry (R=0.67, P<0.0001).

We also investigated potential differences in the number of unique RDs and the number of AFib triggers that induced reentry in each patient-derived model. In all 51 inducible models, the median number of unique RDs was 1.5 for AFib and 2 for ESUS models (range for both groups: 1-5). There was no significant difference in the number of unique reentrant morphologies per model between the two groups (Fig. 3B, P=0.79), CI: [-1.6×10⁻⁵, 1]. The number of unique RDs was positively correlated with atrial fibrosis burden (Fig. 3C, R=0.69, P<0.0001). For all RD-inducible cases, the median number of stimulation sites from which rapid pacing led to RD formation was 2.5 for AFib and 3 for ESUS models (AFib range: 1-11; ESUS range: 1-7). No significant difference was found in the number of pacing sites that induced reentry per model between AFib and ESUS (Fig. 3D, P=0.98, CI [-1,1]). The number of RD inducing pacing sites was also significantly correlated with fibrosis burden (Fig. 3E, R=0.67, P<0.0001).

Arrhythmia Dynamics

Analysis of simulated reentry episodes revealed no qualitative differences in arrhythmia dynamics between AFib and ESUS models. Figure 4 shows examples of RD-perpetuated in silico arrhythmia and instances where stimulation failed to induce reentry for both groups. Of note, this figure highlights two inducible low-fibrosis ESUS models (Fig. 4A: 10.8% fibrosis, RD near the LIPV; Fig. 4B: 10.0% fibrosis, RD in anterior wall) and a non-inducible high-fibrosis ESUS model (Fig. 4C: 19.4% fibrosis). In the latter case, dense fibrosis on the posterior wall resulted in conduction block as indicated. Similarly, Figs. 4D and 4E present examples of RD-driven arrhythmia in AFib models (9.0% and 24.1% fibrosis, respectively), and Fig. 4F shows an AFib model (12.3% fibrosis) in which reentry was not induced due to wavefront collision in the posterior wall region. Overall, this analysis shows that both ESUS and AFib models exhibited activation patterns consistent with previous definitions of RD-driven arrhythmia; examples of inducible low-fibrosis and non-inducible high-fibrosis models emphasize that fibrosis burden alone is an insufficient predictor for a potential arrhythmic substrate.

• Download figure
• Open in new tab

Figure 4.

Maps of fibrotic tissue distribution (left) and activation time (right) for ESUS and AFib models in which pacing succeeded (rows 1-2) or failed (row 3) to induce RD-driven arrhythmia. Black arrows indicate directions of wavefront propagation in RDs. Double lines indicate sites of conduction block. Black-shaded regions in activation maps indicate locations where activation did not occur during the analysis interval. (A) ESUS model with 10.8% fibrosis and reentry inferior to LIPV. (B) ESUS model with 10.0% fibrosis and reentry on anterior wall. (C) ESUS model with 19.4% fibrosis with wavefront termination through fibrosis on posterior wall. (D) AFib model with 9.0% fibrosis and reentry observed adjacent to RIPV on posterior wall. (E) AFib model with 24.1% fibrosis and reentry observed inferior to the LIPV and on posterior wall. (F) AFib model with 12.3% fibrosis with wavefront collision on posterior wall.

For 50/51 inducible models, AFib was driven by persistent singular RDs. In one AFib model we observed unique arrhythmia dynamics previously undocumented in atrial simulations, best described as dual RD feedback and characterized by interaction between two distinct wavefronts wherein one follows in the wake of another (Fig. 5 and Supplemental Video 1). Reentry appeared to occur as a result of a large central refractory region arising from wavefront collision adjacent to the RIPV on the septal face. This reentry morphology was stable for the duration of the simulation.

• Download figure
• Open in new tab

Figure 5.

Sequential activation maps showing one RD cycle of dual RD feedback viewed from posterior (top) and anterior (bottom) perspectives. Black-shaded region indicate locations where activation did not occur during the analysis interval. Activation maps at 5450ms progress to 5900ms in increments of 150ms. Two distinct wavefronts (labeled 1 and 2 in the panels above) exist in this form of reentry and “chase” after one another for duration of the simulation. See also Supplemental Video 1.

Properties of RD Localization

As described in Methods, each model was automatically subdivided into five anatomical regions (see schematic illustrations in Fig. 1D and Supplementary Fig. 1) and region-wise inducibility score analysis (IdS, as described in Methods) was used to gauge likelihood of RD induction in response to rapid electrical stimulation from different locations in the LA. While ESUS and AFib models had a statistically similar pattern of inducibility rates (P=0.73, by χ² test), stimulation from the posterior wall was ∼3 times more likely to induce RDs in AFib models than ESUS models (Fig. 6A, IdS = 8 vs IdS = 3). In other words, with all other factors held equal, our simulations suggest that triggered activity in the posterior wall may be more likely to initiate reentrant arrhythmia in AFib patients compared to ESUS patients. The same IdS values plotted in Fig. 6A were mapped onto representative LA models to facilitate visual comparison of regional sensitivity to rapid pacing (Fig. 6B). Next, we considered the number of unique RD localization sites in each LA region across the different model groups (i.e., AFib vs. ESUS). The LPV region was most likely to harbor RDs in the AFib cohort (Fig. 6C, N=17). In the ESUS cohort, the atrial floor was the most likely region to contain an RD (N=12), whereas in the AFib cohort this was one of the least likely locations to harbor RDs (Fig. 6C, N=6). Of special note for the ESUS cohort, the atrial floor was both the most sensitive to rapid pacing in terms of RD initiation and the most likely area of RD localization (Fig. 6D). The association between the type of model (ESUS vs AFib) and RD localization was not significant (P=0.29, χ² test).

• Download figure
• Open in new tab

Figure 6.

Summary of IdS and RD localization characteristics. (A) Region-wise IdS for both ESUS and AFib LA models (B) Heat map of the regions in which triggers are most likely to induce arrhythmias depicted as representative ESUS and AFib models. (C) Histogram of RDs across all AFib and ESUS models binned by localization to specific LA regions. (D) Heat map of regions in which RDs are most likely to localize depicted as representative ESUS and AFib models.

Inducibility of reentry and fibrosis quantification in patient-derived atrial models

Experimental findings have shown that atrial fibrosis results in changes that promote reentry (8, 9), but the exact mechanism of this connection is not fully understood. Previous modeling studies have linked RD localization to specific spatial patterns of fibrotic remodeling in AFib (5, 10). Recent clinical findings indicate that atrial fibrosis burden does not differ significantly between AFib and ESUS patients, and yet (by definition) ESUS patients do not demonstrate AFib at the time of stroke or during ambulatory monitoring (4). Given the findings summarized above, a potential explanation is that, notwithstanding the fact that ESUS patients have substantial fibrosis, the particular spatial distribution of fibrotic remodeling in their atria is not conducive to arrhythmia perpetuation. Our findings suggest that this is likely not the case. From the standpoint of computational models derived from patient LGE-MRI scans, fibrotic substrate in individuals with ESUS is indistinguishable from that in patients with AFib in terms of the fundamental capacity to sustain RDs. Reentrant arrhythmias were induced in 28/45 AFib models and 23/45 ESUS models (Fig. 2A). Across the combined set of 51 inducible models from both cohorts, we found that high fibrosis models were more likely to exhibit RDs irrespective of whether they corresponded to ESUS or AFib patients. This correlation has been characterized in prior AFib modeling studies (5, 6) and is consistent with established clinical thinking regarding the relationship between fibrosis and AFib outcomes (11). Notably, we did observe several cases in which models defied inducibility expectations based on fibrosis alone. Such models exhibited RDs despite low fibrosis or were non-inducible despite high fibrosis. This observation confirms that, as observed previously in analogous AFib modeling studies (5, 6), assessment of raw fibrosis burden LGE-MRI scans alone is insufficient to fully characterize arrhythmogenic capacity of potentially pro-arrhythmic substrate in ESUS.

Consistent with the goal of this research to understand the contribution of fibrotic substrate to potential RD formation in ESUS, we purposefully excluded arrhythmias perpetuated by other mechanisms from our study design (e.g., self-sustaining activity driven solely by focal sources). Potential contributions from the right atrium (RA) were also excluded, since only the LA was segmented from LGE-MRI as part of the clinical workflow. Either of these factors may explain the absence of simulated arrhythmia in 17 of 45 AFib models, many of which had very little LA fibrosis (i.e., AFib in these individuals might have been predominantly focal in nature or sustained by RDs in the RA). This rate of inducibility is consistent with previous studies (e.g.,13 out of 20 in computational models reconstructed from LGE-MRI scans of persistent AFib patients) (5), supporting the notion that LA fibrosis is associated with increased arrhythmia inducibility but fails to tell the whole story. Importantly, neither of these model constraints repudiates the central finding of our study, which suggests there is no difference between ESUS and AFib patients in terms of the fundamental capacity of the fibrotic substrate to potentially harbor RDs.

Reentrant driver localization dynamics and morphology in patient-derived atrial models

As discussed above, our qualitative findings suggest that ESUS patients’ fibrotic substrate is no different than that of AFib patients in terms the capacity to sustain RDs per se. We performed additional analysis to assess whether specific consequences of fibrotic remodeling influenced any characteristics of RD inducibility in different ways for simulations in models corresponding to ESUS vs. AFib patients. First, we found that there was no significant difference in global fibrosis burden between the 28 inducible AFib models and the 23 inducible ESUS models (Fig. 3A). This finding suggests that intrinsic pro-arrhythmic characteristics of the fibrotic substrate in ESUS and AFib patients are indistinguishable, which is a key finding of this study. Given this result, one potential alternative explanation for the lack of arrhythmia in ESUS patients is a lack of suitable triggers, despite an abundance of fibrotic substrate on par with that observed in AFib. The plausibility of this explanation is strengthened by the fact that the pacing sites from which episodes of reentry were induced in our study were based on common AFib trigger sites as identified in a recent clinical study (12); this is in contrast to previous modeling studies, which simulated triggered activity from evenly-distributed atrial sites (6). However,, the premise of this aspect of our work was not to assess an absence of triggers, but instead ask the following: if the atria of these ESUS patients were subjected to the same type of triggered activity known to occur in typical AFib patients, is it possible the result be would sustained arrhythmia? In more than half of the cohort (23/45), our analysis suggests that the answer is yes.

A potential implication of our study is that in AFib patients, pathological changes in geometry, specifically increased LA surface area (134±40 cm² for AFib vs. 109±26 cm² for ESUS; P=0.0007; see Table 1), and fibrotic remodeling may interact synergistically to promote RD formation. Other measurements (BMI, LA volume index) suggest this difference might be abolished by normalizing these values to body surface area. However, we opted to present absolute LA surface area values since a relevant metric for arrhythmia perpetuation is effective heart size (i.e., ratio between absolute heart size and cardiac wavelength) (13). In effect, larger effective heart size implies that reentrant wavefronts have more “elbow room” to establish and anchor themselves, increasing the likelihood of sustained RD occurrence. Thus, this may explain the trend towards higher burden of fibrosis in inducible ESUS models vs. inducible AFib models (Fig. 3A, P=0.068). Since ESUS patients lack the anatomical LA remodeling observed in AFib patients, a higher threshold of fibrotic tissue may be necessary in order for RD-driven arrhythmia to be even theoretically possible. More work is needed to concretely validate this finding.

Further analysis was performed to probe potential differences between ESUS and AFib patients that went beyond consideration of RD inducibility as a binary variable. Specifically, we found no significant difference (P=0.79) in the number of unique model-predicted RDs between ESUS and AFib (range for both groups: 1-5; Fig. 3B). Instead, the number of unique RDs was highly correlated (Fig. 3C, R=0.69, P<0.0001) with LA fibrosis burden, which is consistent with the concept that high fibrosis models are more likely to exhibit RDs. Secondly, analysis of the number of pacing sites that induced RDs revealed high variability (range: 0-11 out of 15), but again the ESUS and AFib groups were similar in this regard (Fig. 3D, P=0.98). This variable was also correlated to LA fibrosis burden (R=0.67, P<0.0001), suggesting greater fibrotic remodeling leads to susceptibility to pacing-induced reentry by triggered activity from a larger number of locations. This finding further substantiates our principal claim that no significant differences exist between the detected fibrotic substrate in ESUS and AFib, in that it holds true for the substrate’s capacity to sustain reentry and its susceptibility to triggered activity. The general implication is that in the presence of simulated triggered activity, both of these characteristics are closely linked to global fibrosis burden. RDs identified by non-invasive electrocardiographic imaging (ECGI) and in silico phase singularity identification have been shown to co-localize with fibrosis boundary zones identified by LGE-MRI (5, 14-16). RD localization dynamics in this study were consistent with these findings, as illustrated by representative LA fibrotic tissue distributions and corresponding RD activation patterns in Fig. 4. RD morphology in this study, for both AFib and ESUS models, largely corroborated previous findings – arrhythmia episodes were perpetuated by one RD at a time, with activity in the periphery including conduction block, transient reentry, and wavefront collision (5, 6). There was one noteworthy exception, which was the dual reciprocating RD interaction documented in Fig. 5, in which a region of conduction block formed dynamically in the myocardial region between the two circulating wavefronts. Similar excitation sequences have been documented previously in embryonic chick cardiac cells (17) and in patients with atypical flutter (18). It is unclear why we observed this behavior for the first time in this study despite the fact that over 100 episodes of RD-driven arrhythmia have been simulated in dozens of similarly constructed models (5, 6, 19). A potential explanation is the use of increased tissue conductivity in the present study, resulting in faster conduction velocity.

Insights from analysis of RD inducibility and localization by LA region

Our findings thus far demonstrate that the properties of simulated RDs are indistinguishable between ESUS and AFib. We then performed region-wise analysis to assess any systematic differences in the spatial distribution of pro-arrhythmic fibrosis that could render certain regions more susceptible to either triggered activity or RD-localization. Even when analyzing the atria by distinct regions, we found that the spatial properties of the fibrotic substrate between AFib and ESUS patients are not intrinsically different. ESUS and AFib models had no significant association to either reentrant driver localization or region-wise inducibility (P=0.73 and P=0.29, by χ² test). However, our region-wise analysis (Fig. 6) did yield several noteworthy findings that suggest subtle distinctions in the spatial distribution of pro-arrhythmic fibrosis may exist between the two groups. While non-PV LA triggers have been identified (posterior wall, appendage, mitral valve) (12), the PVs are considered the main source of focal activity responsible for initiating AFib episodes (20, 21). In light of this fact, it is notable that our modeling suggests triggered activity in the LPV and RPV regions is ∼50% more likely to induce reentrant arrhythmia in AFib vs. ESUS models (i.e., IdS ≈ 7.5 vs. 5). An even more striking differential exists for the posterior wall region, where the IdS score was ∼3x higher in AFib compared to ESUS models. This observation suggests that even in the cases when common potential LA triggers, such as (i.e., PV or posterior wall ectopy) do occur in ESUS patients, they are less likely to engage the fibrotic substrate and initiate sustained reentry compared to the same activity in AFib patients. This possibility does not contradict our hypothesis that the lack of arrhythmia in ESUS is due to a dearth of triggers; rather, it is a complementary corollary that can be put to the test in future clinical and computational analysis.

While understanding of RD localization dynamics in AFib remains limited, evidence from ECGI mapping indicates that reentrant activity occurs most frequently in the PV and posterior wall regions (15, 20, 21). Our findings are consistent with these data, with those same three areas harboring a majority (∼70%) of all observed RDs in models corresponding to AFib patients. In contrast, in the ESUS population the atrial floor was the RD localization hotspot. The importance of this finding is unclear as tendencies toward reentrant activity in particular LA areas has not been meaningfully correlated to clinical arrhythmia properties. However, it provides a path for future validation studies: if incident AFib in patients who previously presented with ESUS can be characterized by intracardiac mapping, the hypothesis that RDs localize preferentially to the atrial floor can be tested.

Finally, across both ESUS and AFib cohorts we observed that LA regions with the highest IdS generally corresponded to areas most likely to harbor RDs. The implication is that despite some interesting and noteworthy trends in region-wise substrate susceptibility between ESUS and AFib, the overarching conclusion of our analysis remains unchanged: we hypothesize that if the ESUS substrate were subjected to suitable triggered activity, it would be fully capable of sustaining RDs indistinguishable from those that perpetuate AFib.

Limitations

In this study, atrial tissue is modeled as a bilayer to drastically reduce computational load. Previous studies have used this modeling framework (14, 22) to represent human atria effectively, but the framework remains a simplification compared to volumetric 3D models. Moreover, clinical-grade MRI resolution limits our ability to detect fine details in anatomical structure and spatial distribution of potentially arrhythmogenic substrate, for instance slow-conducting tracks of fibrotic atrial tissue that could underlie microreentrant circuits (23). While these models are patient-specific in terms of LA anatomy and each individual’s unique pattern of fibrotic remodeling, they do not incorporate inter-patient variability in CV and electrophysiological properties such as ion channel expression. Nevertheless, our previous analysis indicates that this representation of atrial architecture with generic “average AFib” electrophysiology is appropriate for use in patient-derived modeling (10, 24).

As in previous studies (25, 26), our models do not differentiate between cell- or tissue-scale properties of atrial electrophysiology between patients with paroxysmal and persistent forms of AFib. Likewise, our approach to characterizing potential arrhythmia propensity in ESUS patients assumes cell- and tissue-scale remodeling based on experimental and clinical data from the AFib milieu. Although this is relevant as a limitation and must be considered when interpreting our results, this aspect of our approach is also one of the major advantages of the modeling and simulation methodology. Specifically, it allows us to assess whether there are any relevant differences in the spatial pattern of fibrotic remodeling between ESUS and AFib patients in the absence of other potentially confounding variables.

Finally, the mechanism of stroke in ESUS patients may be independent of the presence of AFib and atrial fibrosis through decreased atrial function might be a contributor to thrombus formation in the absence of AFib. Currently, secondary stroke prophylaxis is dependent on detecting AFib and predicting, through computational modeling, which atria are more prone to manifest AFib may be of clinical value. Another future research direction that could prove highly fruitful in the near future would be to create multi-scale, multi-physics image-based models of the fibrotic atria to assess each individual’s risk of clot formation in a patient-specific manner.

Patient Population

Patients were recruited to undergo cardiac LGE-MRI from the University of Washington (Seattle, WA) and Klinikum Coburg (Coburg, Germany) between July 2016 and June 2019. This study was approved by the Institutional Review Board (IRB) of the University of Washington (UW) and the Ethikkommission der Bayerischen Ländesärztekammer München, Bayern, Deutschland; all participants provided written informed consent. Patients with ESUS met published diagnostic criteria (27). Patients with paroxysmal (27/45, 60%) or persistent AFib and without stroke were recruited from the UW Cardiac Arrhythmia Data Repository, an IRB-approved database for arrhythmia patients. Exclusion criteria for AFib patients included those who had undergone LA catheter ablation before MRI and those with only atrial flutter. Patients with cardiac implantable electronic devices, severe claustrophobia, renal dysfunction, and other contraindications to MRI or gadolinium-based contrast were excluded.

MRI acquisition

Cardiac LGE-MRI was obtained on all participants to quantify the extent of LA fibrosis using previously described protocols (11). Scans were performed on Philips Ingenia and Siemens Avanto clinical scanners, 15 to 25 minutes after contrast injection, using a 3-dimensional inversion-recovery, respiration-navigated, ECG-gated, gradient echo pulse sequence. Acquisition parameters included transverse imaging volume with a voxel size of 1.25 x 1.25 x 2.5 mm (reconstructed to 0.625 x 0.625 x 1.25 mm). Scan time was 5 to 10 minutes dependent on respiration and heart rate. Fat saturation sequences were used to suppress signal from fatty tissue.

Reconstruction of 3D patient-derived atrial models from LGE-MRI

Geometric models were reconstructed from LGE-MRI and the relative extent of fibrosis in the LA was quantified via an adaptive histogram thresholding algorithm (28). Clinical-grade meshes (i.e., coarse discretization) produced by Merisight Inc. (Salt Lake City, UT) were resampled with a target resolution of 200 μm using an automated process based on gmsh (29). Each LA model was represented as a bilayer comprising of nested endocardial and epicardial shells, linked at every point by linear connections (σ = 0.8 S m⁻¹) (Fig. 1A). In each patient-derived model, realistic myocardial fiber orientations were mapped from an atlas geometry (22) using the universal atrial coordinates (UAC) approach. Briefly, this process assigned epicardial and endocardial fibers from a previously published bilayer model to the target atrial geometry (Fig. 1B) (30, 31). In all finite element LA meshes, the average element edge length was ∼188 μm and the number of nodes ranged from ∼600,000 to ∼1.4 million, depending on LA size. This mesh resolution is consistent with previously established benchmarks for minimizing numerical error due to spatial discretization in simulations of cardiac wavefront propagation (32).

Modeling of atrial electrophysiology in fibrotic and non-fibrotic regions

Our methodology for computational modeling at the cell and tissue scale of the fibrotic and non-fibrotic atrial electrophysiology can be found in previously published papers (5, 15, 19). Briefly, in non-fibrotic regions, a human atrial action potential model (33) was used to represent membrane kinetics, including parameter modifications to fit clinical monophasic action potential recordings from AFib patients (I_Kur, I_to, I_CaL decreased by 50%, 50%, 70%, respectively) (5, 34). At the tissue scale, non-fibrotic atrial tissue, conductivity tensor values (longitudinal: σ_L = 0.409 S m⁻¹; transverse: σ_T = 0.08195 S m⁻¹) were calibrated to obtain effective conduction velocity (CV) values of 71.49 cm s⁻¹ and 37.14 cm s⁻¹ (longitudinal and transverse). These conductivities were chosen to match, CV values measured in AFib patients (61±6 cm/s) (35). In fibrotic regions, modifications to the AFib-like action potential model (ICaL, I_Na, and I_K1 decreased by 50%, 40%, 50%, respectively) were implemented as in prior studies (5, 36), resulting in a 15.4% increase in action potential duration and a 49.6% decrease in upstroke velocity. These changes represented the effect of elevated transforming growth factor-β1, a key component of the fibrogenic signaling pathway. As in previous studies (5, 15, 19), tissue-scale effects of interstitial fibrosis and gap junction remodeling were represented by reducing overall conductivity and exaggerating the anisotropy ratio (σ_L:σ_T) from 5:1 to 8:1 (σ_L = 0.17708 S m⁻¹; σ_T = 0.022135 S m⁻¹).

Simulation of electrical activity and numerical aspects

Electrical propagation in LA models was simulated by solving the monodomain equation using the finite element method. This system was coupled with ordinary differential and algebraic equations representing myocyte membrane dynamics at each node in the mesh, as described in the prior section. All simulations were executed on the Hyak supercomputer system at the University of Washington using the CARP software package (37, 38). A free version of CARP is available for academic use (see: https://openCarp.org). The compute time required to complete each unique simulation ranged from 1-10 hours. The total CPU time for all simulations conducted in all models was 13.4 years.

Induction and analysis of reentrant atrial arrhythmias

Simulations were performed to assess the pro-arrhythmic propensity of the fibrotic substrate in each patient-derived model. Arrhythmia induction via rapid pacing was attempted from 15 pacing sites derived from AFib trigger sites (Fig. 1C, see caption for detailed anatomical site descriptions) (12). Clinically relevant AFib trigger sites were chosen over a random pacing schematic to specifically capture RDs that arise from locations demonstrated to induce AFib. Similar to previous publications, a clinically relevant pacing sequence of 12 electrical stimuli was delivered at each of the 15 locations (5, 36). Individual cell-scale ionic models were paced to steady-state at a rate basic cycle length of 500 ms. The electrical stimulus consisted of two initial pulses with a coupling interval of 300ms, followed by pulses ramping down to 200ms in 20ms intervals. After the delivery of the final stimulus, simulations were monitored for self-sustaining electrical wavefront propagation. For all cases in which activity persisted for at least 5,000ms post-pacing, we applied further analysis to determine whether the cause was an induced RD or macroscopic re-entry (i.e. continuous repetitive, self-sustaining activation propagating around a non-conductive obstacle such as the mitral valve or pulmonary vein(s)), which we consider flutter-like reentry. Instances of macroscopic re-entry were excluded from further analysis.

For each AFib-inducible simulation, we documented whether each pacing site induced reentry and analyzed patterns of RD localization. Unique RD morphologies in each patient-derived model were classified as being located in one of five anatomical regions (Fig. 1D). These regions were delineated automatically in a two-step process summarized in Supplemental Fig. 1. First, the LA was subdivided into three broad anatomical areas (region 1: LA floor, 2: posterior wall; 3: anterior wall including LAA) using standardized cutoff values in the UAC space (30). Second, the left and right PV areas (regions 4 and 5, respectively) were established using a region-growing approach such that each accounted for 15% of the total LA surface area. We then defined region-wise inducibility scores (IdS) across all models in a particular group (ESUS or AFib) as the proportion of pacing sites within a given region from which rapid pacing resulted in initiation of an RD. For example, since the LPV region contains four pacing sites (anterior/posterior LSPV/LIPV), the corresponding ESUS IdS value would be derived by summing the number of instances in which pacing from those locations induced RD across inducible ESUS models then dividing by four. This ensured our ability to assess spatial heterogeneity of sensitivity to triggered activity in a manner that was unbiased to the relative abundance of pacing sites in some LA regions.

Statistical Analysis

LA models for ESUS and AFib patients were divided into quartiles based on the extent of fibrotic remodeling as measured by LGE-MRI. Continuous variables were compared pairwise between groups using Wilcoxon rank-sum tests and were reported as mean ± standard deviation. Confidence intervals were calculated as the interval for the true difference in mean with 95% certainty. Categorical variables were compared using a χ² test. After classifying unique RDs and number of pacing sites that induced reentry, correlation with fibrosis was assessed with logistic regression. Statistical significance was established at two-tailed P≤0.05. All statistical analysis was performed using R (39).