Transformer-based spatial-temporal detection of apoptotic cell death in live-cell imaging.

Intravital microscopy has revolutionized live-cell imaging by allowing the study of spatial-temporal cell dynamics in living animals. However, the complexity of the data generated by this technology has limited the development of effective computational tools to identify and quantify cell processes. Amongst them, apoptosis is a crucial form of regulated cell death involved in tissue homeostasis and host defense. Live-cell imaging enabled the study of apoptosis at the cellular level, enhancing our understanding of its spatial-temporal regulation. However, at present, no computational method can deliver robust detection of apoptosis in microscopy timelapses. To overcome this limitation, we developed ADeS, a deep learning-based apoptosis detection system that employs the principle of activity recognition. We trained ADeS on extensive datasets containing more than 10,000 apoptotic instances collected both in vitro and in vivo, achieving a classification accuracy above 98% and outperforming state-of-the-art solutions. ADeS is the first method capable of detecting the location and duration of multiple apoptotic events in full microscopy timelapses, surpassing human performance in the same task. We demonstrated the effectiveness and robustness of ADeS across various imaging modalities, cell types, and staining techniques. Finally, we employed ADeS to quantify cell survival in vitro and tissue damage in mice, demonstrating its potential application in toxicity assays, treatment evaluation, and inflammatory dynamics. Our findings suggest that ADeS is a valuable tool for the accurate detection and quantification of apoptosis in live-cell imaging and, in particular, intravital microscopy data, providing insights into the complex spatial-temporal regulation of this process.

Development and Validation of a Concise Objectifiable Risk Evaluation Score for Non-Relapse Mortality after Allogeneic Hematopoietic Stem Cell Transplantation.

We aimed to develop a concise objectifiable risk evaluation (CORE) tool for predicting non-relapse mortality (NRM) and overall survival (OS) after allogeneic hematopoietic stem cell transplantation (allo-HCT). A total of 1120 adult patients who had undergone allo-HCT at our center between 2013 and 2020 were divided into training, first, and second validation cohorts. Objectifiable, patient-related factors impacting NRM in univariate and multivariate analyses were: serum albumin, serum creatinine, serum C-reactive protein (CRP), heart function (LVEF), lung function (VC, FEV1), and patient age. Hazard ratios were assigned points (0-3) based on their impact on NRM and summed to the individual CORE HCT score. The CORE HCT score stratified patients into three distinct low-, intermediate-, and high-risk groups with two-year NRM rates of 9%, 22%, and 46%, respectively, and OS rates of 73%, 55%, and 35%, respectively (p < 0.001). These findings were confirmed in a first and a second recently treated validation cohort. Importantly, the CORE HCT score remained informative across various conditioning intensities, disease-specific subgroups, and donor types, but did not impact relapse incidence. A comparison of CORE HCT vs. HCT Comorbidity Index (HCT-CI) in the second validation cohort revealed better performance of the CORE HCT score with c-statistics for NRM and OS of 0.666 (SE 0.05, p = 0.001) and 0.675 (SE 0.039, p < 0.001) vs. 0.431 (SE 0.057, p = 0.223) and 0.535 (SE 0.042, p = 0.411), respectively. The CORE HCT score is a concise and objectifiable risk evaluation tool for adult patients undergoing allo-HCT for malignant disease. External multicenter validation is underway.

Fast deep learning reconstruction techniques for preclinical magnetic resonance fingerprinting.

We propose a deep learning (DL) model and a hyperparameter optimization strategy to reconstruct T1 and T2 maps acquired with the magnetic resonance fingerprinting (MRF) methodology. We applied two different MRF sequence routines to acquire images of ex vivo rat brain phantoms using a 7-T preclinical scanner. Subsequently, the DL model was trained using experimental data, completely excluding the use of any theoretical MRI signal simulator. The best combination of the DL parameters was implemented by an automatic hyperparameter optimization strategy, whose key aspect is to include all the parameters to the fit, allowing the simultaneous optimization of the neural network architecture, the structure of the DL model, and the supervised learning algorithm. By comparing the reconstruction performances of the DL technique with those achieved from the traditional dictionary-based method on an independent dataset, the DL approach was shown to reduce the mean percentage relative error by a factor of 3 for T1 and by a factor of 2 for T2 , and to improve the computational time by at least a factor of 37. Furthermore, the proposed DL method enables maintaining comparable reconstruction performance, even with a lower number of MRF images and a reduced k-space sampling percentage, with respect to the dictionary-based method. Our results suggest that the proposed DL methodology may offer an improvement in reconstruction accuracy, as well as speeding up MRF for preclinical, and in prospective clinical, investigations.

A dynamic clamping approach using in silico IK1 current for discrimination of chamber-specific hiPSC-derived cardiomyocytes.

Human induced pluripotent stem cell (hiPSC)-derived cardiomyocytes (CM) constitute a mixed population of ventricular-, atrial-, nodal-like cells, limiting the reliability for studying chamber-specific disease mechanisms. Previous studies characterised CM phenotype based on action potential (AP) morphology, but the classification criteria were still undefined. Our aim was to use in silico models to develop an automated approach for discriminating the electrophysiological differences between hiPSC-CM. We propose the dynamic clamp (DC) technique with the injection of a specific IK1 current as a tool for deriving nine electrical biomarkers and blindly classifying differentiated CM. An unsupervised learning algorithm was applied to discriminate CM phenotypes and principal component analysis was used to visualise cell clustering. Pharmacological validation was performed by specific ion channel blocker and receptor agonist. The proposed approach improves the translational relevance of the hiPSC-CM model for studying mechanisms underlying inherited or acquired atrial arrhythmias in human CM, and for screening anti-arrhythmic agents.

Tracking unlabeled cancer cells imaged with low resolution in wide migration chambers via U-NET class-1 probability (pseudofluorescence).

Cell migration is a pivotal biological process, whose dysregulation is found in many diseases including inflammation and cancer. Advances in microscopy technologies allow now to study cell migration in vitro, within engineered microenvironments that resemble in vivo conditions. However, to capture an entire 3D migration chamber for extended periods of time and with high temporal resolution, images are generally acquired with low resolution, which poses a challenge for data analysis. Indeed, cell detection and tracking are hampered due to the large pixel size (i.e., cell diameter down to 2 pixels), the possible low signal-to-noise ratio, and distortions in the cell shape due to changes in the z-axis position. Although fluorescent staining can be used to facilitate cell detection, it may alter cell behavior and it may suffer from fluorescence loss over time (photobleaching).Here we describe a protocol that employs an established deep learning method (U-NET), to specifically convert transmitted light (TL) signal from unlabeled cells imaged with low resolution to a fluorescent-like signal (class 1 probability). We demonstrate its application to study cancer cell migration, obtaining a significant improvement in tracking accuracy, while not suffering from photobleaching. This is reflected in the possibility of tracking cells for three-fold longer periods of time. To facilitate the application of the protocol we provide WID-U, an open-source plugin for FIJI and Imaris imaging software, the training dataset used in this paper, and the code to train the network for custom experimental settings.

Electromechanical factors associated with favourable outcome in cardiac resynchronization therapy.

AIMS: Electromechanical coupling in patients receiving cardiac resynchronization therapy (CRT) is not fully understood. Our aim was to determine the best combination of electrical and mechanical substrates associated with effective CRT. METHODS AND RESULTS: Sixty-two patients were prospectively enrolled from two centres. Patients underwent 12-lead electrocardiogram (ECG), cardiovascular magnetic resonance (CMR), echocardiography, and anatomo-electromechanical mapping (AEMM). Remodelling was measured as the end-systolic volume (ΔESV) decrease at 6 months. CRT was defined effective with ΔESV ≤ -15%. QRS duration (QRSd) was measured from ECG. Area strain was obtained from AEMM and used to derive systolic stretch index (SSI) and total left-ventricular mechanical time. Total left-ventricular activation time (TLVAT) and transeptal time (TST) were derived from AEMM and ECG. Scar was measured from CMR. Significant correlations were observed between ΔESV and TST [rho = 0.42; responder: 50 (20-58) vs. non-responder: 33 (8-44) ms], TLVAT [-0.68; 81 (73-97) vs. 112 (96-127) ms], scar [-0.27; 0.0 (0.0-1.2) vs. 8.7 (0.0-19.1)%], and SSI [0.41; 10.7 (7.1-16.8) vs. 4.2 (2.9-5.5)], but not QRSd [-0.13; 155 (140-176) vs. 167 (155-177) ms]. TLVAT and SSI were highly accurate in identifying CRT response [area under the curve (AUC) > 0.80], followed by scar (AUC > 0.70). Total left-ventricular activation time (odds ratio = 0.91), scar (0.94), and SSI (1.29) were independent factors associated with effective CRT. Subjects with SSI >7.9% and TLVAT <91 ms all responded to CRT with a median ΔESV ≈ -50%, while low SSI and prolonged TLVAT were more common in non-responders (ΔESV ≈ -5%). CONCLUSION: Electromechanical measurements are better associated with CRT response than conventional ECG variables. The absence of scar combined with high SSI and low TLVAT ensures effectiveness of CRT.

Fast Characterization of Inducible Regions of Atrial Fibrillation Models With Multi-Fidelity Gaussian Process Classification.

Computational models of atrial fibrillation have successfully been used to predict optimal ablation sites. A critical step to assess the effect of an ablation pattern is to pace the model from different, potentially random, locations to determine whether arrhythmias can be induced in the atria. In this work, we propose to use multi-fidelity Gaussian process classification on Riemannian manifolds to efficiently determine the regions in the atria where arrhythmias are inducible. We build a probabilistic classifier that operates directly on the atrial surface. We take advantage of lower resolution models to explore the atrial surface and combine seamlessly with high-resolution models to identify regions of inducibility. We test our methodology in 9 different cases, with different levels of fibrosis and ablation treatments, totalling 1,800 high resolution and 900 low resolution simulations of atrial fibrillation. When trained with 40 samples, our multi-fidelity classifier that combines low and high resolution models, shows a balanced accuracy that is, on average, 5.7% higher than a nearest neighbor classifier. We hope that this new technique will allow faster and more precise clinical applications of computational models for atrial fibrillation. All data and code accompanying this manuscript will be made publicly available at: https://github.com/fsahli/AtrialMFclass.

CANCOL, a Computer-Assisted Annotation Tool to Facilitate Colocalization and Tracking of Immune Cells in Intravital Microscopy.

Two-photon intravital microscopy (2P-IVM) has become a widely used technique to study cell-to-cell interactions in living organisms. Four-dimensional imaging data obtained via 2P-IVM are classically analyzed by performing automated cell tracking, a procedure that computes the trajectories followed by each cell. However, technical artifacts, such as brightness shifts, the presence of autofluorescent objects, and channel crosstalking, affect the specificity of imaging channels for the cells of interest, thus hampering cell detection. Recently, machine learning has been applied to overcome a variety of obstacles in biomedical imaging. However, existing methods are not tailored for the specific problems of intravital imaging of immune cells. Moreover, results are highly dependent on the quality of the annotations provided by the user. In this study, we developed CANCOL, a tool that facilitates the application of machine learning for automated tracking of immune cells in 2P-IVM. CANCOL guides the user during the annotation of specific objects that are problematic for cell tracking when not properly annotated. Then, it computes a virtual colocalization channel that is specific for the cells of interest. We validated the use of CANCOL on challenging 2P-IVM videos from murine organs, obtaining a significant improvement in the accuracy of automated tracking while reducing the time required for manual track curation.

All-at-once multigrid approaches for one-dimensional space-fractional diffusion equations.

We focus on a time-dependent one-dimensional space-fractional diffusion equation with constant diffusion coefficients. An all-at-once rephrasing of the discretized problem, obtained by considering the time as an additional dimension, yields a large block linear system and paves the way for parallelization. In particular, in case of uniform space-time meshes, the coefficient matrix shows a two-level Toeplitz structure, and such structure can be leveraged to build ad-hoc iterative solvers that aim at ensuring an overall computational cost independent of time. In this direction, we study the behavior of certain multigrid strategies with both semi- and full-coarsening that properly take into account the sources of anisotropy of the problem caused by the grid choice and the diffusion coefficients. The performances of the aforementioned multigrid methods reveal sensitive to the choice of the time discretization scheme. Many tests show that Crank-Nicolson prevents the multigrid to yield good convergence results, while second-order backward-difference scheme is shown to be unconditionally stable and that it allows good convergence under certain conditions on the grid and the diffusion coefficients. The effectiveness of our proposal is numerically confirmed in the case of variable coefficients too and a two-dimensional example is given.

Space-time shape uncertainties in the forward and inverse problem of electrocardiography.

In electrocardiography, the "classic" inverse problem is the reconstruction of electric potentials at a surface enclosing the heart from remote recordings at the body surface and an accurate description of the anatomy. The latter being affected by noise and obtained with limited resolution due to clinical constraints, a possibly large uncertainty may be perpetuated in the inverse reconstruction. The purpose of this work is to study the effect of shape uncertainty on the forward and the inverse problem of electrocardiography. To this aim, the problem is first recast into a boundary integral formulation and then discretised with a collocation method to achieve high convergence rates and a fast time to solution. The shape uncertainty of the domain is represented by a random deformation field defined on a reference configuration. We propose a periodic-in-time covariance kernel for the random field and approximate the Karhunen-Loève expansion using low-rank techniques for fast sampling. The space-time uncertainty in the expected potential and its variance is evaluated with an anisotropic sparse quadrature approach and validated by a quasi-Monte Carlo method. We present several numerical experiments on a simplified but physiologically grounded two-dimensional geometry to illustrate the validity of the approach. The tested parametric dimension ranged from 100 up to 600. For the forward problem, the sparse quadrature is very effective. In the inverse problem, the sparse quadrature and the quasi-Monte Carlo method perform as expected, except for the total variation regularisation, where convergence is limited by lack of regularity. We finally investigate an H1/2 regularisation, which naturally stems from the boundary integral formulation, and compare it to more classical approaches.

GEASI: Geodesic-based earliest activation sites identification in cardiac models.

The identification of the initial ventricular activation sequence is a critical step for the correct personalization of patient-specific cardiac models. In healthy conditions, the Purkinje network is the main source of the electrical activation, but under pathological conditions the so-called earliest activation sites (EASs) are possibly sparser and more localized. Yet, their number, location and timing may not be easily inferred from remote recordings, such as the epicardial activation or the 12-lead electrocardiogram (ECG), due to the underlying complexity of the model. In this work, we introduce GEASI (Geodesic-based Earliest Activation Sites Identification) as a novel approach to simultaneously identify all EASs. To this end, we start from the anisotropic eikonal equation modeling cardiac electrical activation and exploit its Hamilton-Jacobi formulation to minimize a given objective function, for example, the quadratic mismatch to given activation measurements. This versatile approach can be extended to estimate the number of activation sites by means of the topological gradient, or fitting a given ECG. We conducted various experiments in 2D and 3D for in-silico models and an in-vivo intracardiac recording collected from a patient undergoing cardiac resynchronization therapy. The results demonstrate the clinical applicability of GEASI for potential future personalized models and clinical intervention.

Automatic reconstruction of the left atrium activation from sparse intracardiac contact recordings by inverse estimate of fibre structure and anisotropic conduction in a patient-specific model.

AIMS: Electric conduction in the atria is direction-dependent, being faster in fibre direction, and possibly heterogeneous due to structural remodelling. Intracardiac recordings of atrial activation may convey such information, but only with high-quality data. The aim of this study was to apply a patient-specific approach to enable such assessment even when data are scarce, noisy, and incomplete. METHODS AND RESULTS: Contact intracardiac recordings in the left atrium from nine patients who underwent ablation therapy were collected before pulmonary veins isolation and retrospectively included in the study. The Personalized Inverse Eikonal Model from cardiac Electro-Anatomical Maps (PIEMAP), previously developed, has been used to reconstruct the conductivity tensor from sparse recordings of the activation. Regional fibre direction and conduction velocity were estimated from the fitted conductivity tensor and extensively cross-validated by clustered and sparse data removal. Electrical conductivity was successfully reconstructed in all patients. Cross-validation with respect to the measurements was excellent in seven patients (Pearson correlation r > 0.93) and modest in two patients (r = 0.62 and r = 0.74). Bland-Altman analysis showed a neglectable bias with respect to the measurements and the limit-of-agreement at -22.2 and 23.0 ms. Conduction velocity in the fibre direction was 82 ± 25 cm/s, whereas cross-fibre velocity was 46 ± 7 cm/s. Anisotropic ratio was 1.91±0.16. No significant inter-patient variability was observed. Personalized Inverse Eikonal model from cardiac Electro-Anatomical Maps correctly predicted activation times in late regions in all patients (r = 0.88) and was robust to a sparser dataset (r = 0.95). CONCLUSION: Personalized Inverse Eikonal model from cardiac Electro-Anatomical Maps offers a novel approach to extrapolate the activation in unmapped regions and to assess conduction properties of the atria. It could be seamlessly integrated into existing electro-anatomic mapping systems. Personalized Inverse Eikonal model from cardiac Electro-Anatomical Maps also enables personalization of cardiac electrophysiology models.

Efficient identification of scars using heterogeneous model hierarchies.

AIMS: Detection and quantification of myocardial scars are helpful for diagnosis of heart diseases and for personalized simulation models. Scar tissue is generally characterized by a different conduction of excitation. We aim at estimating conductivity-related parameters from endocardial mapping data. Solving this inverse problem requires computationally expensive monodomain simulations on fine discretizations. We aim at accelerating the estimation by combining electrophysiology models of different complexity. METHODS AND RESULTS: Distributed parameter estimation is performed by minimizing the misfit between simulated and measured electrical activity on the endocardial surface, subject to the monodomain model and regularization. We formulate this optimization problem, including the modelling of scar tissue and different regularizations, and design an efficient solver. We consider grid hierarchies and monodomain-eikonal model hierarchies in a recursive multilevel trust-region method. With numerical examples, efficiency and estimation quality, depending on the data, are investigated. The multilevel solver is significantly faster than a comparable single level solver. Endocardial mapping data of realistic density appears to be sufficient to provide quantitatively reasonable estimates of location, size, and shape of scars close to the endocardial surface. CONCLUSION: In several situations, scar reconstruction based on eikonal and monodomain models differ significantly, suggesting the use of the more involved monodomain model for this purpose. Eikonal models can accelerate the computations considerably, enabling the use of complex electrophysiology models for estimating myocardial scars from endocardial mapping data.

Fully coupled dynamic simulations of bioprosthetic aortic valves based on an embedded strategy for fluid-structure interaction with contact.

AIMS: This work aims at presenting a fully coupled approach for the numerical solution of contact problems between multiple elastic structures immersed in a fluid flow. The key features of the computational model are (i) a fully coupled fluid-structure interaction with contact, (ii) the use of a fibre-reinforced material for the leaflets, (iii) a stent, and (iv) a compliant aortic root. METHODS AND RESULTS: The computational model takes inspiration from the immersed boundary techniques and allows the numerical simulation of the blood-tissue interaction of bioprosthetic heart valves (BHVs) as well as the contact among the leaflets. First, we present pure mechanical simulations, where blood is neglected, to assess the performance of different material properties and valve designs. Secondly, fully coupled fluid-structure interaction simulations are employed to analyse the combination of haemodynamic and mechanical characteristics. The isotropic leaflet tissue experiences high-stress values compared to the fibre-reinforced material model. Moreover, elongated leaflets show a stress concentration close to the base of the stent. We observe a fully developed flow at the systolic stage of the heartbeat. On the other hand, flow recirculation appears along the aortic wall during diastole. CONCLUSION: The presented FSI approach can be used for analysing the mechanical and haemodynamic performance of a BHV. Our study suggests that stresses concentrate in the regions where leaflets are attached to the stent and in the portion of the aortic root where the BHV is placed. The results from this study may inspire new BHV designs that can provide a better stress distribution.

Synergistic antiarrhythmic effect of inward rectifier current inhibition and pulmonary vein isolation in a 3D computer model for atrial fibrillation.

AIMS: Recent clinical studies showed that antiarrhythmic drug (AAD) treatment and pulmonary vein isolation (PVI) synergistically reduce atrial fibrillation (AF) recurrences after initially successful ablation. Among newly developed atrial-selective AADs, inhibitors of the G-protein-gated acetylcholine-activated inward rectifier current (IKACh) were shown to effectively suppress AF in an experimental model but have not yet been evaluated clinically. We tested in silico whether inhibition of inward rectifier current or its combination with PVI reduces AF inducibility. METHODS AND RESULTS: We simulated the effect of inward rectifier current blockade (IK blockade), PVI, and their combination on AF inducibility in a detailed three-dimensional model of the human atria with different degrees of fibrosis. IK blockade was simulated with a 30% reduction of its conductivity. Atrial fibrillation was initiated using incremental pacing applied at 20 different locations, in both atria. IK blockade effectively prevented AF induction in simulations without fibrosis as did PVI in simulations without fibrosis and with moderate fibrosis. Both interventions lost their efficacy in severe fibrosis. The combination of IK blockade and PVI prevented AF in simulations without fibrosis, with moderate fibrosis, and even with severe fibrosis. The combined therapy strongly decreased the number of fibrillation waves, due to a synergistic reduction of wavefront generation rate while the wavefront lifespan remained unchanged. CONCLUSION: Newly developed blockers of atrial-specific inward rectifier currents, such as IKAch, might prevent AF occurrences and when combined with PVI effectively supress AF recurrences in human.

Learning atrial fiber orientations and conductivity tensors from intracardiac maps using physics-informed neural networks.

Electroanatomical maps are a key tool in the diagnosis and treatment of atrial fibrillation. Current approaches focus on the activation times recorded. However, more information can be extracted from the available data. The fibers in cardiac tissue conduct the electrical wave faster, and their direction could be inferred from activation times. In this work, we employ a recently developed approach, called physics informed neural networks, to learn the fiber orientations from electroanatomical maps, taking into account the physics of the electrical wave propagation. In particular, we train the neural network to weakly satisfy the anisotropic eikonal equation and to predict the measured activation times. We use a local basis for the anisotropic conductivity tensor, which encodes the fiber orientation. The methodology is tested both in a synthetic example and for patient data. Our approach shows good agreement in both cases and it outperforms a state of the art method in the patient data. The results show a first step towards learning the fiber orientations from electroanatomical maps with physics-informed neural networks.

In Vivo Motility Patterns Displayed by Immune Cells Under Inflammatory Conditions.

The migration of immune cells plays a key role in inflammation. This is evident in the fact that inflammatory stimuli elicit a broad range of migration patterns in immune cells. Since these patterns are pivotal for initiating the immune response, their dysregulation is associated with life-threatening conditions including organ failure, chronic inflammation, autoimmunity, and cancer, amongst others. Over the last two decades, thanks to advancements in the intravital microscopy technology, it has become possible to visualize cell migration in living organisms with unprecedented resolution, helping to deconstruct hitherto unexplored aspects of the immune response associated with the dynamism of cells. However, a comprehensive classification of the main motility patterns of immune cells observed in vivo, along with their relevance to the inflammatory process, is still lacking. In this review we defined cell actions as motility patterns displayed by immune cells, which are associated with a specific role during the immune response. In this regard, we summarize the main actions performed by immune cells during intravital microscopy studies. For each of these actions, we provide a consensus name, a definition based on morphodynamic properties, and the biological contexts in which it was reported. Moreover, we provide an overview of the computational methods that were employed for the quantification, fostering an interdisciplinary approach to study the immune system from imaging data.

Large scale simulation of pressure induced phase-field fracture propagation using Utopia.

Non-linear phase field models are increasingly used for the simulation of fracture propagation problems. The numerical simulation of fracture networks of realistic size requires the efficient parallel solution of large coupled non-linear systems. Although in principle efficient iterative multi-level methods for these types of problems are available, they are not widely used in practice due to the complexity of their parallel implementation. Here, we present Utopia, which is an open-source C++ library for parallel non-linear multilevel solution strategies. Utopia provides the advantages of high-level programming interfaces while at the same time a framework to access low-level data-structures without breaking code encapsulation. Complex numerical procedures can be expressed with few lines of code, and evaluated by different implementations, libraries, or computing hardware. In this paper, we investigate the parallel performance of our implementation of the recursive multilevel trust-region (RMTR) method based on the Utopia library. RMTR is a globally convergent multilevel solution strategy designed to solve non-convex constrained minimization problems. In particular, we solve pressure-induced phase-field fracture propagation in large and complex fracture networks. Solving such problems is deemed challenging even for a few fractures, however, here we are considering networks of realistic size with up to 1000 fractures.

Reconstruction of three-dimensional biventricular activation based on the 12-lead electrocardiogram via patient-specific modelling.

AIMS: Non-invasive imaging of electrical activation requires high-density body surface potential mapping. The nine electrodes of the 12-lead electrocardiogram (ECG) are insufficient for a reliable reconstruction with standard inverse methods. Patient-specific modelling may offer an alternative route to physiologically constraint the reconstruction. The aim of the study was to assess the feasibility of reconstructing the fully 3D electrical activation map of the ventricles from the 12-lead ECG and cardiovascular magnetic resonance (CMR). METHODS AND RESULTS: Ventricular activation was estimated by iteratively optimizing the parameters (conduction velocity and sites of earliest activation) of a patient-specific model to fit the simulated to the recorded ECG. Chest and cardiac anatomy of 11 patients (QRS duration 126-180 ms, documented scar in two) were segmented from CMR images. Scar presence was assessed by magnetic resonance (MR) contrast enhancement. Activation sequences were modelled with a physiologically based propagation model and ECGs with lead field theory. Validation was performed by comparing reconstructed activation maps with those acquired by invasive electroanatomical mapping of coronary sinus/veins (CS) and right ventricular (RV) and left ventricular (LV) endocardium. The QRS complex was correctly reproduced by the model (Pearson's correlation r = 0.923). Reconstructions accurately located the earliest and latest activated LV regions (median barycentre distance 8.2 mm, IQR 8.8 mm). Correlation of simulated with recorded activation time was very good at LV endocardium (r = 0.83) and good at CS (r = 0.68) and RV endocardium (r = 0.58). CONCLUSION: Non-invasive assessment of biventricular 3D activation using the 12-lead ECG and MR imaging is feasible. Potential applications include patient-specific modelling and pre-/per-procedural evaluation of ventricular activation.