Skip to main content

Comparative study of left atrium epicardial fat tissue pattern using persistent homology approach

Abstract

Objective

Atrial Fibrillation (A-fib) is an abnormal heartbeat condition in which the heart races and beats in an uncontrollable way. It is observed that the presence of increased epicardial fat/fatty tissue in the atrium can lead to A-fib. Persistent homology using topological features can be used to recapitulate enormous amounts of spatially complicated medical data into a visual code to identify a specific pattern of epicardial fat tissue with non-fat tissue. Our aim is to evaluate the topological pattern of left atrium epicardial fat tissue with non-fat tissue.

Results

A topological data analysis approach was acquired to study the imaging pattern between the left atrium epicardial fat tissue and non-fat tissue patches. The patches of eight patients from CT images of the left atrium heart were used and categorized into “left atrium epicardial fat tissue” and “non-fat tissue” groups. The features that distinguish the “epicardial fat tissue” and “non-fat tissue” groups are extracted using persistent homology (PH). Our result reveals that our proposed research can discriminate between left atrium epicardial fat tissue and non-fat tissue. Specifically, the range of Betti numbers in the epicardial tissue is smaller (0–30) than the non-fat tissue (0–100), indicating that non-fat tissue has good topology.

Introduction

Atrial fibrillation (A-fib) is a heart disorder characterized by an abnormal heart beat in which the heart races and beats in an unpredictable way. It is seen that A-fib results in the development of blood clots in the upper chamber of the heart which later flow to different organs and confine the blood supply to tissues and elicit heart attacks. It is the most prevalent type of arrhythmia, and causes a higher mortality rate worldwide. It is assumed the mortality and morbidity related to A-fib can significantly increase up to three-fold within the next 30 years [1,2,3]. The presence of elevated epicardial fat exists in patients with permanent A-fib [4]. The epicardial fat is found on the top of the left atrium close to the left atrial appendage and lateral to the mitral isthmus and is more in the superior half of the left atrium than in the inferior half [4, 5]. The peri-atrial epicardial fat may contribute to the formation of A-fib by obstructing the flow of electrical signals to all regions of the heart [4]. Here in this study, our main concern will be on the epicardial fat present in the left atrium. In past few decades, medical imaging data has been piling up in abundance. Medical imaging provides tremendously meaningful information which can be used to extract patterns and can be saved for future reference. With the recent up gradation of computational methods, researchers have offered a variety of concepts and tools that may help to recognize a specific disease pattern through topological features. Topological Data Analysis (TDA) gives a fresh perspective in this field to analyze the datasets using topological approaches. The fundamental goal of TDA is to build tools for studying qualitative characteristics of data using ideas and results from geometry and topology [6]. PH is one of the TDA approaches for detecting topological features of data. This method is based on algebraic computation, which provides an effective theoretical framework for understand qualitative characteristics of data with complex structures [6]. Here, we introduce a method that allows the extraction of topological features that differentiate “left atrium epicardial fat tissue” and “non-fat tissue” group patients.

Main text

Methods

Data information and preparation

This work included eight male patients with A-fib in 52–61 age groups. Our work was based on delayed enhancement cardiac CT image datasets obtained by ten sets of time frames at the same position using different contrast mediums, each containing one complete cardiac cycle. Philips computerized tomography instrument was utilized to obtain the image dataset of. Each patient's data set included 392 images with 512 × 512 dimensions that covered the entire thoracic region.

All the 392 images were stacked up to get a 3D view of the thoracic region. Then we extracted only those images which included only the heart. These images were 204 in number. Since the thoracic region contains the heart, ribcage and spinal cord, which can make difficult the visualization of the heart, that is the rationale for our clearing out the surrounding structures by removing the pixels representing these structures. The 3D view was visualized through the Volume Viewer App of MATLAB R2018a. After that, a masking filter was applied to these images to highlight the required area by changing the color of the original images. An image display threshold setting of –190 to –30 Hounsfield Units (HU) was identified for epicardial fat on grayscale; other research also suggested a threshold of -200 to -50 HU [5, 7]. The images were sent to mask the pixel values which lie between -190 to -30 HU. The fat tissues on the atrium of the heart were identified using the pixel value range masked over the entire dataset. The above-mentioned methodology has already been described in our previously published study [8]. We extracted the desired dataset of 36 images from each patient’s data which contains the fat region of the atrium. We cropped these fat regions from all the images and made 36 patches of 32 × 32 dimensions. These 36 patches were concatenated to form one large patch. Additionally another 36 patches are extracted from the atrium containing non-fat tissues and preceded the same way as mentioned (Fig. 1).

Persistent homology filtration

Key objects in topological data analysis are filtered simplicial complexes, called filtrations. Homology is a method for in an accurate measuring the shape of a geometric object by counting holes or features of various dimensions. Its output is a vector space for each non-negative integer k whose dimension (called the k-th Betti number) is the feature count for holes of dimension k. This notion extends to filtered simplicial complex by producing a “persistence vector space”, the analogue of vector spaces, as well as the analogue of dimension, which is called a persistence barcode or persistence diagram. These PH diagrams can be thought of as finite collections (unordered) of intervals, and one uses various algebraic combinations of the lengths of the intervals as well as their midpoints [9].

Filtration is the first step of PH which produces a series of simplicial complexes for a scale proximity parameters (ε). We observe the filtration procedure for a point cloud data set where each point is surrounded by a sphere of radius ε. We draw an edge between two spots at each intersection of two spheres. Filtering data creates a simplicial complex space from which PH quantifies the presence of n-dimensional holes, which include 0-dimensional holes, 1-dimensional holes, circles/loops/tunnels, and 2-dimensional holes. Since the best value for the scale ε cannot be determined, the primary principle of PH is to move through all possible values (0) to see how the homology of these components changes [10]. We evaluate the times of birth (ε emerges) and death (ε vanishes) for each n-dimensional structure (Fig. 2).

Fig. 1
figure 1

(i) Patch representation of epicardial fat tissue. (ii) Patch representation of Non-fat tissue

Fig. 2
figure 2

Topological data analysis workflow with imaging signal input

Fig. 3
figure 3

(i) PH of epicardial fat tissue. (ii) PH of Non-fat tissue.

Results and discussion

In left atrium epicardial fat tissue, we see the range of Betti numbers varies less (0–30), while in the non-fat tissue, the range of the Betti numbers is large (0–100). After Betti number, we plotted PH. Figure 3 shows the PH diagrams of left atrium epicardial fat and non-fat tissue. In the diagram of epicardial fat tissue, the points are closer towards the diagonal line, which means a topology with many small holes, but in non-fat tissue, the points are more scattered and concentrated away from diagonal that means topology is smoother on a small scale. This is the other way to calculate the pattern between two groups of tissues.

Fat tissues are the adipose tissue that helps in storing energy, and in distress, they provide support to the system for the proper execution of body functions. They insulate the body and act as an endocrine organ. The adipose tissue between the visceral pericardium and the myocardium is known as epicardial fat [11]. The increased epicardial fat on the left atrium wall hampers electrical conduction. Several studies have been reported that provide evidence of epicardial fat in the atrium causing A-fib [12, 13]. To measure cardiac fats, a couple of imaging modalities such as magnetic resonance imaging (MRI), echocardiography, and computed tomography (CT) are currently accessible in the market. However, detection of the fat can need either manual strategies which can be tedious or using some available strategies [14,15,16,17]. TDA combines algebraic topology and statistical learning techniques to provide a mathematical foundation for studying the shape of data. TDA also offers dimensionality reduction and noise stability. Here, we applied PH, an algebraic method of TDA that discover the topology of data to find new and distinctive features. The basic idea behind PH is to replace data points with a parametrized family of simplicial complexes, which can generally be depicted as a union of points, edges, triangles, tetrahedrons, and higher-dimensional polytopes, and encode the change of the simplicial complexes' topological features (such as the number of connected components, holes, and voids) across various parameters for data analysis [18,19,20]. Our PH diagrams represent the distinct comparison between fat tissue and non-fat tissues in our data. Our results have shown that this method may help in identifying the fat tissue from non-fat tissue for better stratification. In conclusion, by analyzing CT patches using a topological data analysis approach known as persistent homology, we have identified patterns in barcodes and persistence diagrams that discriminate cardiac patients whoexperience epicardial fat tissue versus those who do not have epicardial fat tissue. While promising, our results will need to be validated on a larger cohort.

Limitations of the study

We have used very small number of the patients with A-fib. The data was validated under the supervision of cardiologists.

Availability of data and materials

The dataset of the current study is available from the first author on reasonable request.

Abbreviations

A-fib:

Atrial Fibrillation

PH:

Persistent homology

TDA:

Topological Data Analysis

HU:

Hounsfield Units

References

  1. Heeringa J, van der Kuip DA, Hofman A, Kors JA, van Herpen G, Stricker BH, Stijnen T, Lip GY, Witteman JC. Prevalence, incidence and lifetime risk of atrial fibrillation: the Rotterdam study. Eur Heart J. 2006;27:949–53.

    Article  Google Scholar 

  2. Schnabel RB, Yin X, Gona P, Larson MG, Beiser AS, McManus DD, NewtonCheh C, Lubitz SA, Magnani JW, Ellinor PT, Seshadri S, Wolf PA, Vasan RS, Benjamin EJ, Levy D. 50 year trends in atrial fibrillation prevalence, incidence, risk factors, and mortality in the Framingham Heart Study: a cohort study. Lancet. 2015;386:154–62.

    Article  Google Scholar 

  3. Willems S, Meyer C, de Bono J, Brandes A, Eckardt L, Elvan A, et al. Cabins, castles, and constant hearts: rhythm control therapy in patients with atrial fibrillation. Eur Heart J. 2019;40(46):3793–3799c.

    Article  Google Scholar 

  4. Al Chekakie MO, Akar JG. Epicardial fat and atrial fibrillation: a review. J Atr Fibrillation 2012; 4(6).

  5. Hsuan-Ming T, et al. Quantitative analysis of quantity and distribution of epicardial adipose tissue surrounding the left atrium in patients with atrial fibrillation and effect of recurrence after ablation. Am J Cardiol. 2011;107(10):1498–503.

    Article  Google Scholar 

  6. Otter N, Porter MA, Tillmann U, Grindrod P, Harrington HA. A roadmap for the computation of persistent homology. EPJ Data Sci. 2017;6:1–38.

    Article  Google Scholar 

  7. Al Chekakie MO, et al. Pericardial fat is independently associated with human atrial fibrillation. J Am Coll Cardiol. 2010;56(10):784–8.

    Article  Google Scholar 

  8. Deepa D, Singh Y, Wang MC, Hu W. An automated method for detecting atrial fat using convolutional neural network. Proc Inst Mech Eng Part H J Eng Med. 2021;235(11):1329–34.

    Article  Google Scholar 

  9. Wasserman L. Topological data analysis. Ann Rev Stat Appl. 2018;5:501–32.

    Article  Google Scholar 

  10. Zomorodian A, Carlsson G. Computing persistent homology. Discret Comput Geom. 2005;33(2):249–74.

    Article  Google Scholar 

  11. Sacher F, Roberts-Thomson K, Maury P, Tedrow U, Nault I, Steven D, et al. Epicardial ventricular tachycardia ablation: a multicenter safety study. J Am Coll Cardiol. 2010;55(21):2366–72.

    Article  Google Scholar 

  12. Yang F, et al. Towards a mechanistic understanding and treatment of a progressive disease: atrial fibrillation. J Atr Fibrillation. 2017. https://doi.org/10.4022/jafib.1627.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Singh Y, Shakyawar D, Weichih Hu. Non-ischemic endocardial scar geometric remodeling toward topological machine learning. Proc Inst Mech Eng [H]. 2020;234(9):1029–35.

    Article  Google Scholar 

  14. Singh Y, et al. Geometrical evaluation of the Scar in Left ventricle using TDA. In: 2021 IEEE 11th annual computing and communication workshop and conference (CCWC). IEEE; 2021.

  15. Polonsky TS, et al. Coronary artery calcium score and risk classification for coronary heart disease prediction. JAMA. 2010;303(16):1610–6.

    CAS  Article  Google Scholar 

  16. Rodrigues ÉO, et al. A novel approach for the automated segmentation and volume quantification of cardiac fats on computed tomography. Comput Methods Programs Biomed. 2016;123:109–28.

    CAS  Article  Google Scholar 

  17. Singh Y, Jons W, Conte GM, Jagtap J, Zhang K, Sobek JD, Rouzrokh P, Eaton JE, Erickson BJ. Persistent homology approach distinguishes potential pattern between “Early” and “Not Early” hepatic decompensation groups using MRI modalities. Curr Dir Biomed Eng. 2021;7(2):488–91.

    Article  Google Scholar 

  18. Ghrist R. Barcodes: the persistent topology of data. Bull Am Math Soc. 2008;45(1):61–75.

    Article  Google Scholar 

  19. Aktas ME, Akbas E, El Fatmaoui A. Persistence homology of networks: methods and applications. Appl Netw Sci. 2019;4(1):1–28.

    Article  Google Scholar 

  20. Adams H, Johnathan B, Joshua M. Topological data analysis. Data science for mathematicians. Chapman and Hall/CRC; 2020. p. 441–474.

Download references

Acknowledgements

The authors are thankful to Hospital Yang-Ming University Taiwan, CYCU, Biomedical department and Dr. Carlsson.

Funding

The study was funded by Biomedical engineering, Chung Yuan Christian University, Taiwan and the University of Dubai, Dubai.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization, DD, YS, and WH; methodology, YS, GC, RP; investigation, WH,YS, GC; writing—original draft preparation, D,WM, YS; writing—review and editing YS, GC, RP; project administration, WH. All authors have read and approved the final manuscript.

Corresponding authors

Correspondence to Yashbir Singh or Wathiq Mansoor.

Ethics declarations

Ethics approval and consent to participate

The data was provided by the National Institute of Hospital Yang-Ming University Taiwan, and the Institutional Review Board approved the study. Patients signed consent was also obtained.

Consent for publication

Not applicable.

Competing interests

The authors have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Deepa, D., Singh, Y., Mansoor, W. et al. Comparative study of left atrium epicardial fat tissue pattern using persistent homology approach. BMC Res Notes 15, 299 (2022). https://doi.org/10.1186/s13104-022-06173-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13104-022-06173-2

Keywords

  • Atrial fibrillation
  • Barcode
  • Persistent homology
  • Pixel value masking
  • Topology