A robust computational framework for estimating 3D Bi-Atrial chamber wall thickness

Improved algorithm for generating evenly-spaced streamlines from an orientation field on a triangulated surface

BACKGROUND

Vector fields such as cardiac fiber orientation can be visualized on a surface using streamlines. The application of evenly-spaced streamline generation to the construction of interconnected cable structure for cardiac propagation models has more stringent requirements imperfectly fulfilled by current algorithms.

METHOD

We developed an open-source C++/python package for the placement of evenly-spaced streamlines on a triangulated surface. The new algorithm improves upon previous works by more accurately handling streamline extremities, U-turns and limit cycles, by providing stronger geometrical guarantees on inter-streamline minimal distance, particularly when a high density of streamlines (up to 10μm spacing) is desired, and by making a more efficient parallel implementation available. The approach requires finding intersections between geometrical capsules and triangles to update an occupancy mask defined on the triangles. This enables fast streamline integration from thousands of seed points to identify optimal streamline placement.

RESULTS

The algorithm was assessed qualitatively on different left atrial models of fiber orientation with varying mesh resolutions (up to 375k triangles) and quantitatively by measuring streamline lengths and distribution of inter-streamline minimal distance. The complexity and the computational performance of the algorithm were studied as a function of streamline spacing in relation to triangular mesh resolution.

CONCLUSION

More accurate geometrical computations, attention to details and fine-tuning led to an algorithm more amenable to applications that require precise positioning of streamlines.

Computing thickness of irregularly-shaped thin walls using a locally semi-implicit scheme with extrapolation to solve the Laplace equation: Application to the right ventricle

Cardiac Magnetic Resonance (CMR) Imaging is currently considered the gold standard imaging modality in cardiology. However, it is accompanied by a tradeoff between spatial resolution and acquisition time. Providing accurate measures of thin walls relative to the image resolution may prove challenging. One such anatomical structure is the cardiac right ventricle. Methods for measuring thickness of wall-like anatomical structures often rely on the Laplace equation to provide point-to-point correspondences between both boundaries. This work presents limex, a novel method to solve the Laplace equation using ghost nodes and providing extrapolated values, which is tested on three different datasets: a mathematical phantom, a set of biventricular segmentations from CMR images of ten pigs and the database used at the RV Segmentation Challenge held at MICCAI'12. Thickness measurements using the proposed methodology are more accurate than state-of-the-art methods, especially with the coarsest image resolutions, yielding mean L1 norms of the error between 43.28% and 86.52% lower than the second-best methods on the different test datasets. It is also computationally affordable. Limex has outperformed other state-of-the-art methods in classifying RV myocardial segments by their thickness.

Chamber-specific wall thickness features in human atrial fibrillation

Persistent atrial fibrillation (AF) is not effectively treated due to a lack of adequate tools for identifying patient-specific AF substrates. Recent studies revealed that in 30-50% of patients, persistent AF is maintained by localized drivers not only in the left atrium (LA) but also in the right atrium (RA). The chamber-specific atrial wall thickness (AWT) features underlying AF remain elusive, though the important role of AWT in AF is widely acknowledged. We aimed to provide direct evidence of the existence of distinguished RA and LA AWT features underlying AF drivers by analysing functionally and structurally mapped human hearts ex vivo. Coronary-perfused intact human atria (n = 7, 47 ± 14 y.o.; two female) were mapped using panoramic near-infrared optical mapping during pacing-induced AF. Then the hearts were imaged at approximately 170 µm3 resolution by 9.4 T gadolinium-enhanced MRI. The heart was segmented, and 3D AWT throughout atrial chambers was estimated and analysed. Optical mapping identified six localized RA re-entrant drivers in four hearts and four LA drivers in three hearts. All RA AF drivers were anchored to the pectinate muscle junctions with the crista terminalis or atrial walls. The four LA AF drivers were in the posterior LA. RA (n = 4) with AF drivers were thicker with greater AWT variation than RA (n = 3) without drivers (5.4 ± 2.6 mm versus 5.0 ± 2.4 mm, T-test p < 0.05; F-test p < 0.05). Furthermore, AWT in RA driver regions was thicker and varied more than in RA non-driver regions (5.1 ± 2.5 mm versus 4.4 ± 2.2 mm, T-test p < 0.05; F-test p < 0.05). On the other hand, LA (n = 3) with drivers was thinner than the LA (n = 4) without drivers. In particular, LA driver regions were thinner than the rest of LA regions (3.4 ± 1.0 mm versus 4.2 ± 1.0 mm, T-test p < 0.05). This study demonstrates chamber-specific AWT features of AF drivers. In RA, driver regions are thicker and have more variable AWT than non-driver regions. By contrast, LA drivers are thinner than non-drivers. Robust evaluation of patient-specific AWT features should be considered for chamber-specific targeted ablation.

A novel method to estimate adult age from the lumbar vertebral body using 3D PMCT images in Japanese

This study evaluated the age-related changes in the vertebral body using 3D Postmortem CT (PMCT) images and proposed an alternative age estimation formula. The PMCT images of 200 deceased individuals aged 25 to 99 years (126 males, 74 females) were retrospectively reviewed and included in the study. Using the open-source software ITK-SNAP and MeshLab, a 3D surface mesh of the fourth lumbar vertebral body (L4) and its convex hull models were created from the PMCT data. Using their inbuilt tools, volumes (in mm3) of the L4 surface mesh and convex hull models were subsequently computed. We derived VD, defined as the difference in volumes between the convex hull and L4 surface mesh normalized by L4 mesh volume, and VR, defined as the ratio of L4 mesh volume to convex hull volume based on individual L4. Correlation and regression analyses were performed between VD, VR, and chronological age. A statistically significant positive correlation (P < 0.001) between chronological age and VD, (r_s = 0.764, males; r_s = 0.725, females), and a significant negative correlation between chronological age and VR (r_s = -0.764, males; r_s = -0.725, females) was obtained in both sexes. The lowest standard error of the estimate was demonstrated by the VR at 11.9 years and 12.5 years for males and females, respectively. As such, their regression models to estimate adult age were Age = 248.9-2.5VR years, males; Age = 258.1-2.5VR years, females. These regression equations may be useful for estimating age in Japanese adults in forensic settings.

Medical image analysis on left atrial LGE MRI for atrial fibrillation studies: A review

Late gadolinium enhancement magnetic resonance imaging (LGE MRI) is commonly used to visualize and quantify left atrial (LA) scars. The position and extent of LA scars provide important information on the pathophysiology and progression of atrial fibrillation (AF). Hence, LA LGE MRI computing and analysis are essential for computer-assisted diagnosis and treatment stratification of AF patients. Since manual delineations can be time-consuming and subject to intra- and inter-expert variability, automating this computing is highly desired, which nevertheless is still challenging and under-researched. This paper aims to provide a systematic review on computing methods for LA cavity, wall, scar, and ablation gap segmentation and quantification from LGE MRI, and the related literature for AF studies. Specifically, we first summarize AF-related imaging techniques, particularly LGE MRI. Then, we review the methodologies of the four computing tasks in detail and summarize the validation strategies applied in each task as well as state-of-the-art results on public datasets. Finally, the possible future developments are outlined, with a brief survey on the potential clinical applications of the aforementioned methods. The review indicates that the research into this topic is still in the early stages. Although several methods have been proposed, especially for the LA cavity segmentation, there is still a large scope for further algorithmic developments due to performance issues related to the high variability of enhancement appearance and differences in image acquisition.

A Transmural Path Model Improves the Definition of the Orthotropic Tissue Structure in Heart Simulations

In the past decades, the structure of the heart, human as well as other species, has been explored in a detailed way, e.g., via histological studies or diffusion tensor magnetic resonance imaging. Nevertheless, the assignment of the characteristic orthotropic structure in a patient-specific finite element model remains a challenging task. Various types of rule-based models, which define the local fiber and sheet orientation depending on the transmural depth, have been developed. However, the correct assessment of the transmural depth is not trivial. Its accuracy has a substantial influence on the overall mechanical and electrical properties in rule-based models. The main purpose of this study is the development of a finite element-based approach to accurately determine the transmural depth on a general unstructured grid. Instead of directly using the solution of the Laplace problem as the transmural depth, we make use of a well-established model for the assessment of the transmural thickness. It is based on two hyperbolic first-order partial differential equations for the definition of a transmural path, whereby the transmural thickness is defined as the arc length of this path. Subsequently, the transmural depth is determined based on the position on the transmural path. Originally, the partial differential equations were solved via finite differences on structured grids. In order to circumvent the need of two grids and mapping between the structured (to determine the transmural depth) and unstructured (electromechanical heart simulation) grids, we solve the equations directly on the same unstructured tetrahedral mesh. We propose a finite-element-based discontinuous Galerkin approach. Based on the accurate transmural depth, we assign the local material orientation of the orthotropic tissue structure in a usual fashion. We show that this approach leads to a more accurate definition of the transmural depth. Furthermore, for the left ventricle, we propose functions for the transmural fiber and sheet orientation by fitting them to literature-based diffusion tensor magnetic resonance imaging data. The proposed functions provide a distinct improvement compared to existing rules from the literature.

Understanding PITX2-Dependent Atrial Fibrillation Mechanisms through Computational Models

Atrial fibrillation (AF) is a common arrhythmia. Better prevention and treatment of AF are needed to reduce AF-associated morbidity and mortality. Several major mechanisms cause AF in patients, including genetic predispositions to AF development. Genome-wide association studies have identified a number of genetic variants in association with AF populations, with the strongest hits clustering on chromosome 4q25, close to the gene for the homeobox transcription PITX2. Because of the inherent complexity of the human heart, experimental and basic research is insufficient for understanding the functional impacts of PITX2 variants on AF. Linking PITX2 properties to ion channels, cells, tissues, atriums and the whole heart, computational models provide a supplementary tool for achieving a quantitative understanding of the functional role of PITX2 in remodelling atrial structure and function to predispose to AF. It is hoped that computational approaches incorporating all we know about PITX2-related structural and electrical remodelling would provide better understanding into its proarrhythmic effects leading to development of improved anti-AF therapies. In the present review, we discuss advances in atrial modelling and focus on the mechanistic links between PITX2 and AF. Challenges in applying models for improving patient health are described, as well as a summary of future perspectives.

Comprehensive evaluation of electrophysiological and 3D structural features of human atrial myocardium with insights on atrial fibrillation maintenance mechanisms

Atrial fibrillation (AF) occurrence and maintenance is associated with progressive remodeling of electrophysiological (repolarization and conduction) and 3D structural (fibrosis, fiber orientations, and wall thickness) features of the human atria. Significant diversity in AF etiology leads to heterogeneous arrhythmogenic electrophysiological and structural substrates within the 3D structure of the human atria. Since current clinical methods have yet to fully resolve the patient-specific arrhythmogenic substrates, mechanism-based AF treatments remain underdeveloped. Here, we review current knowledge from in-vivo, ex-vivo, and in-vitro human heart studies, and discuss how these studies may provide new insights on the synergy of atrial electrophysiological and 3D structural features in AF maintenance. In-vitro studies on surgically acquired human atrial samples provide a great opportunity to study a wide spectrum of AF pathology, including functional changes in single-cell action potentials, ion channels, and gene/protein expression. However, limited size of the samples prevents evaluation of heterogeneous AF substrates and reentrant mechanisms. In contrast, coronary-perfused ex-vivo human hearts can be studied with state-of-the-art functional and structural technologies, such as high-resolution near-infrared optical mapping and contrast-enhanced MRI. These imaging modalities can resolve atrial arrhythmogenic substrates and their role in reentrant mechanisms maintaining AF and validate clinical approaches. Nonetheless, longitudinal studies are not feasible in explanted human hearts. As no approach is perfect, we suggest that combining the strengths of direct human atrial studies with high fidelity approaches available in the laboratory and in realistic patient-specific computer models would elucidate deeper knowledge of AF mechanisms. We propose that a comprehensive translational pipeline from ex-vivo human heart studies to longitudinal clinically relevant AF animal studies and finally to clinical trials is necessary to identify patient-specific arrhythmogenic substrates and develop novel AF treatments.

Artificial Intelligence and Machine Learning in Arrhythmias and Cardiac Electrophysiology

Artificial intelligence (AI) and machine learning (ML) in medicine are currently areas of intense exploration, showing potential to automate human tasks and even perform tasks beyond human capabilities. Literacy and understanding of AI/ML methods are becoming increasingly important to researchers and clinicians. The first objective of this review is to provide the novice reader with literacy of AI/ML methods and provide a foundation for how one might conduct an ML study. We provide a technical overview of some of the most commonly used terms, techniques, and challenges in AI/ML studies, with reference to recent studies in cardiac electrophysiology to illustrate key points. The second objective of this review is to use examples from recent literature to discuss how AI and ML are changing clinical practice and research in cardiac electrophysiology, with emphasis on disease detection and diagnosis, prediction of patient outcomes, and novel characterization of disease. The final objective is to highlight important considerations and challenges for appropriate validation, adoption, and deployment of AI technologies into clinical practice.

Mini Review: Deep Learning for Atrial Segmentation From Late Gadolinium-Enhanced MRIs

Segmentation and 3D reconstruction of the human atria is of crucial importance for precise diagnosis and treatment of atrial fibrillation, the most common cardiac arrhythmia. However, the current manual segmentation of the atria from medical images is a time-consuming, labor-intensive, and error-prone process. The recent emergence of artificial intelligence, particularly deep learning, provides an alternative solution to the traditional methods that fail to accurately segment atrial structures from clinical images. This has been illustrated during the recent 2018 Atrial Segmentation Challenge for which most of the challengers developed deep learning approaches for atrial segmentation, reaching high accuracy (>90% Dice score). However, as significant discrepancies exist between the approaches developed, many important questions remain unanswered, such as which deep learning architectures and methods to ensure reliability while achieving the best performance. In this paper, we conduct an in-depth review of the current state-of-the-art of deep learning approaches for atrial segmentation from late gadolinium-enhanced MRIs, and provide critical insights for overcoming the main hindrances faced in this task.

Unified mechanism of local drivers in a percolation model of atrial fibrillation

The mechanisms of atrial fibrillation (AF) are poorly understood, resulting in disappointing success rates of ablative treatment. Different mechanisms defined largely by different atrial activation patterns have been proposed and, arguably, this dispute has slowed the progress of AF research. Recent clinical evidence suggests a unifying mechanism of local drivers based on sustained reentrant circuits in the complex atrial architecture. Here, we present a percolation inspired computational model showing spontaneous emergence of AF that strongly supports, and gives a theoretical explanation for, the clinically observed diversity of activation. We show that the difference in surface activation patterns is a direct consequence of the thickness of the discrete network of heart muscle cells through which electrical signals percolate to reach the imaged surface. The model naturally follows the clinical spectrum of AF spanning sinus rhythm, paroxysmal AF, and persistent AF as the decoupling of myocardial cells results in the lattice approaching the percolation threshold. This allows the model to make the prediction that, for paroxysmal AF, reentrant circuits emerge near the endocardium, but in persistent AF they emerge deeper in the bulk of the atrial wall. If experimentally verified, this may go towards explaining the lowering ablation success rate as AF becomes more persistent.

Identifying Potential Re-Entrant Circuit Locations From Atrial Fibre Maps

Re-entrant circuits have been identified as potential drivers of atrial fibrillation (AF). In this paper, we develop a novel computational framework for finding the locations of re-entrant circuits from high resolution fibre orientation data. The technique follows a statistical approach whereby we generate continuous fibre tracts across the tissue and couple adjacent fibres stochastically if they are within a given distance of each other. By varying the connection distance, we identify which regions are most susceptible to forming re-entrant circuits if muscle fibres are uncoupled, through the action of fibrosis or otherwise. Our results highlight the sleeves of the pulmonary veins, the posterior left atrium and the left atrial appendage as the regions most susceptible to re-entrant circuit formation. This is consistent with known risk locations in clinical AF. If the model can be personalised for individual patients undergoing ablation, future versions may be able to suggest suitable ablation targets.