Next Article in Journal
Predicting Successful Weaning through Sonographic Measurement of the Rapid Shallow Breathing Index
Previous Article in Journal
HLA-A29 Negative Birdshot-like Chorioretinopathy Associated with Vitiligo—Case Report
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computer-Assisted Algorithm for Quantification of Fibrosis by Native Cardiac CT: A Pilot Study

by
Diana Gonciar
1,
Alexandru-George Berciu
2,*,
Eva-Henrietta Dulf
2,3,
Rares Ilie Orzan
1,
Teodora Mocan
4,5,
Alex Ede Danku
2,
Noemi Lorenzovici
2 and
Lucia Agoston-Coldea
1
1
2nd Department of Internal Medicine, Iuliu Hațieganu University of Medicine and Pharmacy, 400012 Cluj-Napoca, Romania
2
Automation Department, Faculty of Automation and Computer Science, Energy Transition Research Center, Technical University of Cluj-Napoca, 400114 Cluj-Napoca, Romania
3
Physiological Controls Research Center, University Research and Innovation Center, Obuda University, 1034 Budapest, Hungary
4
Physiology Department, “Iuliu Hațieganu” University of Medicine and Pharmacy, 400012 Cluj-Napoca, Romania
5
Department of Nanomedicine, Regional Institute of Gastroenterology and Hepatology, 400158 Cluj-Napoca, Romania
*
Author to whom correspondence should be addressed.
J. Clin. Med. 2024, 13(16), 4807; https://doi.org/10.3390/jcm13164807
Submission received: 4 August 2024 / Revised: 8 August 2024 / Accepted: 9 August 2024 / Published: 15 August 2024
(This article belongs to the Section Cardiology)

Abstract

:
Background/Objectives: Recent advances in artificial intelligence, particularly in cardiac imaging, can potentially enhance patients’ diagnosis and prognosis and identify novel imaging markers. We propose an automated, computer-aided algorithm utilizing native cardiac computed tomography (CT) imaging to identify myocardial fibrosis. This study aims to evaluate its performance compared to CMR markers of fibrosis in a cohort of patients diagnosed with breast cancer. Methods: The study included patients diagnosed with early HER2+ breast cancer, who presented LV dysfunction (LVEF < 50%) and myocardial fibrosis detected on CMR at the time of diagnosis. The patients were also evaluated by cardiac CT, and the extracted images were processed for the implementation of the automatic, computer-assisted algorithm, which marked as fibrosis every pixel that fell within the range of 60–90 HU. The percentage of pixels with fibrosis was subsequently compared with CMR parameters. Results: A total of eight patients (n = 8) were included in the study. High positive correlations between the algorithm’s result and the ECV fraction (r = 0.59, p = 0.126) and native T1 (r = 0.6, p = 0.112) were observed, and a very high positive correlation with LGE of the LV(g) and the LV-LGE/LV mass percentage (r = 0.77, p = 0.025; r = 0.81, p = 0.015). A very high negative correlation was found with GLS (r = −0.77, p = 0.026). The algorithm presented an intraclass correlation coefficient of 1 (95% CI 0.99–1), p < 0.001. Conclusions: The present pilot study proposes a novel promising imaging marker for myocardial fibrosis, generated by an automatic algorithm based on native cardiac CT images.

1. Introduction

Myocardial fibrosis represents the excessive accumulation of collagen in the extracellular space and is classified into two main types: replacement and diffuse. Due to the heart’s limited regenerative capacity, extensive cardiomyocyte necrosis triggers a significant tissue repair process, resulting in the formation of a fibrous scar (replacement fibrosis), which is essential to prevent the rupture of the ischemic myocardium [1]. In contrast, diffuse fibrosis arises from various physiological factors (such as aging) and pathological conditions (including volume overload, pressure overload, inflammatory factors, and cardiomyopathies) [2]. It affects the interstitium or the perivascular area, leading to the onset or aggravation of heart failure (HF) due to impaired microanatomy and function [3]. Identifying and quantifying diffuse fibrosis is critical, not only because it holds prognostic value but also because it is potentially reversible with appropriate management [3].
Cardiotoxicity among patients with breast cancer and human epidermal growth factor receptor 2 (HER2) gene amplification, primarily manifested as asymptomatic left ventricular (LV) dysfunction, is a well-known complication of the standard treatment. This condition is induced by both anthracyclines [4] and targeted therapy with monoclonal antibodies (trastuzumab) [5]. The cumulative incidence of HF and cardiomyopathy over three years increased from 26.7% to 28.2% with the simultaneous administration of these two drugs [6]. Although echocardiography is used for the initial assessment, cardiac magnetic resonance (CMR) imaging provides superior sensitivity and reproducibility, allowing for the identification of myocardial fibrosis that precedes clinical cardiac dysfunction [7].
CMR is central to the multimodal imaging workup, providing anatomical and functional evaluation and tissue characterization. Late gadolinium enhancement (LGE) is the standard non-invasive tool for evaluating replacement fibrosis [8], resulting from extracellular space expansion or insufficient drainage [9]. However, LGE is not suitable for evaluating diffuse fibrosis due to poor spatial resolution and the potential absence of unaffected myocardium [3]. Native CMR T1 mapping effectively highlights interstitial fibrosis through prolonged T1 relaxation times, correlating well with histology assessment [10]. This technique is particularly useful in patients with chronic renal disease, where gadolinium-based contrast agents cannot be used [8], and it offers reduced time and costs [10].
Myocardial strain represents the percentage change in the size of a myocardial segment compared to its initial size [11]. Global longitudinal strain (GLS) is a marker of subclinical myocardial dysfunction, observed to change in patients undergoing anti-HER2 therapy [12] and those receiving anthracyclines [13]. In vivo studies have shown that GLS correlates with the collagen volume fraction [14] and has a predictive role in the mortality of patients with cardiomyopathy [15].
Despite its advantages, the clinical use of CMR is limited by its reduced availability, high costs, and contraindications [16]. Consequently, computed tomography (CT) is gaining attention as an alternative for evaluating myocardial function and tissue characterization, in addition to its established role in assessing coronary anatomy. Several studies have reported a good correlation between the extracellular volume (ECV) fraction determined by CT and CMR [17,18]. However, a significant limitation is the increased dose of contrast medium required [16], and to our knowledge, no methods have been reported that emphasize myocardial fibrosis on native cardiac CT (CCT).
In recent years, automated diagnostic methods based on machine learning have gained popularity and have been successfully used in areas such as detection of colorectal polyps [19], celiac disease [20], or even cancer [21].
The present paper proposes a novel computer-aided diagnostic method for identifying and quantifying cardiac fibrosis on native CCT images. The study aimed to test its performance on a retrospective cohort of patients diagnosed with HER2-positive breast cancer.

2. Materials and Methods

2.1. Study Population

An observational study was conducted involving patients with HER2-positive early breast cancer and a diagnosis of LV dysfunction associated with chemotherapy. The study was conducted at the 2nd Medical Clinic of the County Emergency Hospital Cluj-Napoca and the Hiperdia Diagnostic Imaging Centre Cluj-Napoca between October 2021 and October 2022, in accordance with the principles of the Declaration of Helsinki. The study protocol was approved by the Ethics Committee of the “Iuliu Haţieganu” University of Medicine and Pharmacy Cluj-Napoca (approval number 257/30.06.2021). Eligible patients were informed about the study protocol and provided written informed consent.
Inclusion criteria: Patients were included if they had LV dysfunction (LV ejection fraction (LVEF) ≤ 50%) and the presence of myocardial fibrosis on CMR, were over 18 years old, and were recently diagnosed with early-stage breast cancer (stage I to IIIA) HER2 positive, and treated with anthracycline-based chemotherapy and trastuzumab.
Exclusion criteria: Patients were excluded if they: (1) refused to participate in the study after being informed; (2) had poor acoustic windows; (3) had contraindications for CMR examination (e.g., cardiac stimulators, implantable cardioverter defibrillators, cardiac resynchronization therapy devices, metal implants or impacted metal objects); (4) had a life expectancy of less than one year due to the progression of myocardial disease or other causes; (5) had a history of significant cardiovascular diseases (e.g., cardiac malformations, myocardial infarction, valvulopathies, myocarditis, cardiomyopathies; atrial fibrillation); (6) had other conditions that could contribute to myocardial fibrosis and influence the results (e.g., history of cardiac surgery); (7) had severe renal insufficiency or severe allergy to contrast agents; (8) had a history of radio- and/or chemotherapy for a prior malignancy; or (9) were pregnant or breastfeeding.
Demographic data collected for all patients included age, sex, height, weight, medical history, cardiovascular symptoms (chest pain, dyspnea, syncope, palpitations), and current medication. Additionally, a 12-lead electrocardiogram (ECG) was recorded. Twenty-four-hour ECG Holter monitoring and transthoracic echocardiography were also performed. CCT was conducted within the first 10 days after CMR.

2.2. CMR Protocol and Image Analysis

CMR was performed using a 1.5 Tesla imaging system (Magnetom Altea, Siemens Medical Systems Solutions, Erlanger, Germany), ECG-gated, and conducted during complete apnea, following a standard protocol in accordance with international guidelines [22].
The function and the mass of the LV were evaluated in longitudinal and short axes, from the base to the apex, using steady-state free precession (SSFP) sequences. Standardized acquisition parameters ensured adequate blood-cavitary wall discrimination: repetition time (RT) = 3.6 ms, eco-time (ET) = 1.8 ms, flip angle = 60°, section thickness = 6 mm, field of view = 360 mm, image matrix = 192 × 192 pixels, voxel size = 1.9 × 1.9 × 6 mm, and temporal resolution = 25–40 ms reconstructed to 25 phases of the cardiac cycle [23].
The endocardial and epicardial regions were demarcated in telesystole and telediastole using Syngo.via software (version VB20A_HF08, Argus, Siemens Medical Solutions, Oakville, ON, Canada).
A modified inversion recovery Look-Locker sequence was used for T1 mapping, focusing on one mid-ventricular short-axis slice. T1 relaxation time was determined by establishing a conservative region of interest in the myocardial septum [24]. The limit for native T1 times was considered to be 950 ± 21 ms [25].
Images for T1 measurements were obtained before and after an intravenous infusion of 0.2 mmol/kg gadobutrol (Gadovist, Bayer Schering Pharma AG, Berlin, Germany) was administered. The ECV was calculated using the formula E C V = 1 h a e m a t o c r i t × [ ( 1 p T 1 m 1 n T 1 m ) / ( 1 p T 1 b 1 n T 1 b ) ] , where nT1 represents native T1, pT1 is post-contrast T1, m refers to myocardium, and b refers to blood pool, with reference values of 25 ± 4% [25].

2.3. CCT Protocol and Image Analysis

CCT images were acquired using a second-generation single-source CT scanner (Siemens SOMATOM Definition Edge, Siemens Healthcare, Erlangen, Germany) with a 0.6 mm slice thickness and 1.5 mm reconstruction interval [26]. The X-ray tube voltage ranged from 70 to 140 kV, the tube current was between 500 and 650 mAs per rotation, and the section thickness was 3 mm.
Pre-contrast calcium scoring was assessed semi-automatically using Syngo Calcium Scoring software (CT VC28, Siemens AG, Erlangen, Germany). The Agatston score was calculated, defining a lesion as at least four contiguous pixels, an area greater than 1 mm2, and a density of ≥130 Hounsfield Units (HU) [27].
Coronary images were acquired covering the entire cardiac cycle following an intravenous injection of 1.5 mL/kg of contrast medium, divided into two boluses: 80–100 mL of non-ionic contrast medium (Omnipaque 350 mgI/mL, GE Healthcare, Princeton, NJ, USA) administered with a dual-head power injector (SCT 210; Medrad, Warrendale, PA, USA), followed by 50–80 mL of saline solution. Data acquisition started 5 s after reaching the signal attenuation threshold. For patients with heart rates below 65 bpm, a prospective ECG-gated protocol was implemented, scanning 70–80% of the R-R interval. A postcontrast cardiac CT scan was performed 20 min later with parameters identical to those used for the pre-contrast calcium score scan.
All CCT and CMR interpretations were performed by two Level III-trained experts who were blinded to the clinical data.

2.4. Algorithm Protocol

A novel algorithm was developed by the researchers from the ADAPTED Research Group, part of the Energy Transition Research Center within the Technical University of Cluj-Napoca, Romania. The initial phase involves enumerating the patients within the database and assigning unique identification numbers to each, ensuring a comprehensive index essential for subsequent image analysis and processing.
The number of CT images available for each patient is identified, cataloging the volume of images per patient to provide a clear understanding of the data scope for each individual. The original images are read in DICOM format, a standardized format in medical imaging that contains the necessary information for further processing and analysis.
The application of the fibrosis detection algorithm to the DICOM images involves several sub-steps. The DICOM images, formatted as a 3D numerical matrix, are converted into a matrix that includes HU data, which will be utilized in the fibrosis detection process. Heart contour detection is performed using the numerical matrix. This begins with determining a customized set of starting parameters for each image through a series of manual tests, defining the area of interest where the Snake algorithm will search for the desired contour. The Snake algorithm, implemented using Python version 3.11.7 and the scikit-image library version 0.24.0 (specifically the active_contour function), is configured with the following parameters: a. Gaussian filter with a sigma of 3 and preserve_range set to False; b. Initial contour from the previous step; c. Alpha: 0.15, Beta: 5, Gamma: 0.001; d. 10 iterations. The initial contour is represented as a red dashed line, and the detected contour is depicted as a solid blue line. The region of interest (ROI) is defined as the area between the minimum and maximum coordinates where the contour detected by the Snake algorithm evolves, thereby extracting the heart region as a rectangle from the original image.
Using the detected ROI, the fibrosis detection algorithm is applied at the pixel level. A copy of the image showing the heart contour is created, and the ROI is traversed pixel by pixel. Fibrosis is identified if a pixel’s value falls within the range of 60–90 HU, an interval established through repeated comparative measurements between normal and pathological myocardium. Pixels identified as having fibrosis are marked with a different color in the new image created from the initially identified ROI.

2.5. Statistics

Data distribution was assessed using the Kolmogorov–Smirnov test. Results are presented as mean ± standard deviation (SD) or 95% confidence interval (CI) if normally distributed, as median (25–75% interquartile range [IQR]) if non-normally distributed, or as percentages. Correlation coefficients (r) were interpreted as follows: >0.7 as very high correlation, 0.5–0.7 as high correlation, 0.3–0.5 as medium correlation, 0.1–0.3 as low correlation, and <0.1 no correlation. Differences between CMR and CT groups were assessed using the paired Student t-test.
To evaluate the reliability of the proposed algorithm, the intraclass correlation coefficient (ICC) was determined. ICC values were interpreted as follows: <0.50 indicated poor reliability, 0.50–0.75 suggest moderate reliability, 0.75–0.90 indicated good reliability, and >0.90 indicated excellent reliability [28]. Repeatability was assessed using the coefficient of variance, determined by the root mean square method.
A p-value < 0.05 was considered statistically significant. All statistical analyses were performed using DATAtab: Online Statistics Calculator (DATAtab e.U. Graz, Austria. URL https://datatab.net) and MedCalc Statistical Software version 22.032 (MedCalc Software Ltd., Ostend, Belgium; https://www.medcalc.org; 2020).

3. Results

A total of eight female patients were included in the study (Figure 1). The mean age of the patients was 54.88 ± 12.05, with an age range of 36–70 years.
Dyspnea was the most prevalent symptom, observed in 75% of the study population. The clinical characteristics of the population are summarized in Table 1.
The parameters determined by CMR and CCT are presented in Table 2. At the moment of diagnosis, the mean LVEF was 46.3% ± 12.05 by CMR and 48.8% ± 1.1 by CCT, with 87.2% of the patients presenting LGE, predominantly in the septal area (3/8).
The results obtained from applying the developed algorithm to the CT images of patient 1 are presented in Figure 2. The left column illustrates the regions automatically identified by the algorithm as containing fibrotic tissue, indicated by the pixel value, measured in HU, within the range of 60–90 HU. The second column displays the heart region detection using the Snake algorithm applied to the entire CT image. The starting region is marked in red, while the final region, determined after completing the specific iterations, is outlined in blue. Notably, the algorithm successfully identifies both the heart region and areas of fibrosis across different images of the same patient, demonstrating its consistency and reliability when applied to native CT images.
In this situation, the ROI, and also the returned percentages of fibrosis, will be slightly different (Table 3). The coefficient of variance in this case ranged between 1.93% and 15.56%.
To evaluate the robustness of the algorithm, the effect was analyzed of increasing the number of iterations by 10 across three successive runs of the Snake algorithm for heart region detection and ROI identification. The initial value was determined through manual testing, and for iterations 2 and 3, the increment was held constant. The results obtained from this analysis are presented in Figure 3, focusing on the first patient and their first CT.
After running the algorithm three times, on the same selected ROI, by the same observer (Table 4), an intraclass correlation coefficient was found of 1 (95% CI 0.99–1), p < 0.001. The coefficient of variance ranged between 0.08% and 5.2%.
For further evaluation, the mean of three fibrosis percentage determinations was calculated for each patient across three different CT images (Table 3). Pearson correlation analysis indicated a high positive correlation between the fibrosis percentage determined by the algorithm and ECV fraction (r = 0.59, p = 0.126), as well as native T1 (r = 0.6, p = 0.112). A very high positive correlation was reported between the algorithm’s results and the LGE of the LV(g) and the LV-LGE/LV mass percentage (r = 0.77, p = 0.025; r = 0.81, p = 0.015) (Figure 4). Additionally, the analysis revealed a very high negative correlation between the GLS percentage and the fibrosis percentage determined by the algorithm (r = −0.77, p = 0.026) (Figure 5), which was also observed between GLS percentage and ECV fraction (r = −0.82, p = 0.014). The correlation between NT-proBNP levels and the fibrosis percentage, ECV fraction, native T1, LGE(g), or LV-LGE/LV mass percentage was not statistically significant (p = 0.697, p = 0.716, p = 0.911, p = 0.744, p = 0.842).

4. Discussions

In the present study, patients with LV dysfunction, defined by a reduction in LVEF, were included. The observed LVEF values ranged from 39.2% to 50%, consistent with a HF with mid-interval LVEF [29]. Data from the literature indicate that LVEF can change as early as two weeks after the initiation of chemotherapy [30], with possible mechanisms including increased vascular resistance or aortic stiffening secondary to the antitumor treatment [31].
In the context of chemotherapy-induced cardiotoxicity, myocardial strain emerges as the earliest altered parameter, often preceding changes in LVEF [32]. Longitudinal strain is represented by negative values, indicating longitudinal shortening of the heart from base to apex [33]. In our study, the GLS was −12.25 ± 3.87%, which is above the pathological threshold of −17%. GLS is particularly sensitive to lesions in the subepicardial and subendocardial regions and has been correlated with the extent of LGE, suggesting it may serve as an indirect parameter for fibrosis [33]. Consistent with our findings, a high statistically significant negative correlation was observed between GLS and LGE, as well as GLS and LGE/LV mass (Pearson correlation coefficients r = −0.79, p = 0.02; r = −0.83, p = 0.01).
LGE patterns can be broadly classified into ischemic and non-ischemic types [34]. However, in the context of anticancer therapy, it has been noted that the pattern and location of LGE do not conform to a single model; any pattern can manifest in any location [35]. Our study similarly did not reveal a clear predilection for a specific LGE model.
In patients treated with anthracyclines, T1 relaxation time can change relatively early, even before reductions in LVEF are observed [4,36], and increased ECV fraction is directly correlated with anthracycline administration [37]. Changes in these parameters among chemotherapy-treated patients occur independently of tumor presence or other cardiovascular risk factors [38]. In our study, both elevated T1 (1114 ± 51 ms) and ECV (31.2 ± 3.1%) were observed, indicating the presence of diffuse fibrosis, although we did not achieve a significant correlation with decreased LVEF (p = 0.527). The ECV fraction showed a statistically significant correlation with GLS, supporting previous reports of an association between myocardial strain and diffuse fibrosis in patients with lower LVEF [39].
The automatic algorithm proposed in this study demonstrated reliability, with excellent correlation in repeated measurements performed by the same observer. The coefficient of variation analysis revealed an average percentage of 1.52 ± 1.65%, indicative of high precision. This algorithm, utilizing native CT images, effectively highlighted myocardial fibrosis in HER2-positive breast cancer patients, exhibiting a very high and statistically significant correlation with the amount of LGE, along with a strong correlation with ESV and T1 relaxation time. Notably, the percentage of fibrosis generated by the algorithm correlates with both focal and diffuse fibrosis, and the statistically significant correlation with GLS may be attributed to underlying fibrosis influencing myocardial strain.
In recent years, determining ECV via CT has gained considerable attention, yet it requires a specialized protocol involving delayed-phase acquisitions, which are limited by the use of contrast agents and associated risks, as well as increased radiation exposure [40]. The implementation of a deep-learning algorithm has the potential to partially improve these shortcomings, being able to highlight fibrosis by evaluating early-enhanced CCT sequences [41]. Overall, artificial intelligence solutions have increasingly supported clinicians and patients, enabling automatic segmentation, lesion identification, image reconstruction, and assessment of coronary stenosis in CT images [42].
There are several limitations to our study. Firstly, the study included only eight female patients, which limits the generalizability of the findings and may affect the statistical power; however, it is the first to propose an automatic method that can be performed retrospectively on CT without requiring contrast agents or a special acquisition sequence. Further research with larger cohorts is essential to validate this method comprehensively and to ascertain the specific type of fibrosis (focal vs. diffuse) that this parameter may represent in the future. The study included only female patients with HER2-positive breast cancer, which may limit the applicability of the findings to other populations, including males or patients with different cancer types. Another important limitation of the study is the absence of correlation with histopathological findings, which is why it cannot be known if the percentage calculated in the 60–90 HU range is exclusively caused by fibrosis and not by another associated process, such as microvascular ischemia. The study relies on specific imaging modalities (CT and CMR) for the detection of myocardial fibrosis. Variability in imaging quality or protocols may impact the results. The identification of LGE plays an important prognostic role, both in ischemic and non-ischemic cardiomyopathies, correlating with an increased odds ratio for sudden cardiac death [34]. Given the observed correlation between the fibrosis percentage determined by the algorithm and LGE, it would be valuable to assess in future studies the extent to which fibrosis identification on native CT sequences can provide similar prognostic information. However, the study aimed to offer preliminary data needed for prospective sample size calculations. For an alpha error of 0.05, using a beta error threshold of 0.20, and based on fibrosis percentages obtained for the patients in this study, we calculated that a minimum of 62 patients would need to be enrolled for a final validation of our algorithm. Given the observed correlation between the fibrosis percentage determined by the algorithm and LGE, we strongly feel such a final validation would be a valuable assessment by offering fibrosis identification on native CT sequences with similar prognostic information.

5. Conclusions

The present pilot study demonstrates encouraging results for the newly developed, automatic, computer-assisted method that proposes a novel imaging marker for identifying and quantifying myocardial fibrosis. A significant advantage of this innovative approach is that it extracts data from native CT images, making it a faster and more cost-effective alternative to CMR. The method was tested on a small cohort of HER2-positive breast cancer patients, evaluated using both CT and CMR, and exhibited optimal reliability along with high to very high correlations with CMR-derived fibrosis parameters.

Author Contributions

Conceptualization, E.-H.D. and L.A.-C.; Methodology, D.G., A.-G.B. and T.M.; Software, A.-G.B., E.-H.D., A.E.D. and N.L.; Validation, D.G., R.I.O., T.M. and L.A.-C.; Data curation, A.E.D.; Writing—original draft, D.G., A.-G.B., A.E.D. and N.L.; Writing—review and editing, E.-H.D., R.I.O., T.M. and L.A.-C.; Supervision, L.A.-C. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by “Iuliu Hatieganu” University of Medicine and Pharmacy, Cluj-Napoca, Romania, contract no.2461/29/17.01.2020. Alexandru G. Berciu acknowledge project 38 PFE in the frame of the programme PDI-PFE-CDI 2021.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Ethics Committee of “Iuliu Hatieganu” University of Medicine and Pharmacy (257/30.06.2021).

Informed Consent Statement

Written informed consent has been obtained from the patients to publish this paper.

Data Availability Statement

The data presented in this study are available on request from the corresponding author due to privacy constraints.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Talman, V.; Ruskoaho, H. Cardiac Fibrosis in Myocardial Infarction—From Repair and Remodeling to Regeneration. Cell Tissue Res. 2016, 365, 563–581. [Google Scholar] [CrossRef] [PubMed]
  2. Bengel, F.M.; Diekmann, J.; Hess, A.; Jerosch-Herold, M. Myocardial Fibrosis: Emerging Target for Cardiac Molecular Imaging and Opportunity for Image-Guided Therapy. J. Nucl. Med. 2023, 64, 49S–58S. [Google Scholar] [CrossRef] [PubMed]
  3. López, B.; Ravassa, S.; Moreno, M.U.; José, G.S.; Beaumont, J.; González, A.; Díez, J. Diffuse Myocardial Fibrosis: Mechanisms, Diagnosis and Therapeutic Approaches. Nat. Rev. Cardiol. 2021, 18, 479–498. [Google Scholar] [CrossRef]
  4. Cardinale, D.; Colombo, A.; Bacchiani, G.; Tedeschi, I.; Meroni, C.A.; Veglia, F.; Civelli, M.; Lamantia, G.; Colombo, N.; Curigliano, G.; et al. Early Detection of Anthracycline Cardiotoxicity and Improvement With Heart Failure Therapy. Circulation 2015, 131, 1981–1988. [Google Scholar] [CrossRef]
  5. Payne, D.L.; Nohria, A. Prevention of Chemotherapy Induced Cardiomyopathy. Curr. Heart Fail. Rep. 2017, 14, 398–403. [Google Scholar] [CrossRef] [PubMed]
  6. Chen, J.; Long, J.B.; Hurria, A.; Owusu, C.; Steingart, R.M.; Gross, C.P. Incidence of Heart Failure or Cardiomyopathy After Adjuvant Trastuzumab Therapy for Breast Cancer. J. Am. Coll. Cardiol. 2012, 60, 2504–2512. [Google Scholar] [CrossRef]
  7. Ferreira de Souza, T.; Silva, T.Q.A.C.; Neilan, T.G.; Coelho-Filho, O.R. Assessment of Cardiotoxicity of Cancer Chemotherapy: The Value of Cardiac MRI. Magn. Reson. Imaging Clin. N. Am. 2019, 27, 533–544. [Google Scholar] [CrossRef]
  8. Haaf, P.; Garg, P.; Messroghli, D.R.; Broadbent, D.A.; Greenwood, J.P.; Plein, S. Cardiac T1 Mapping and Extracellular Volume (ECV) in Clinical Practice: A Comprehensive Review. J. Cardiovasc. Magn. Reson. 2016, 18, 89. [Google Scholar] [CrossRef] [PubMed]
  9. Aquaro, G.D.; De Gori, C.; Faggioni, L.; Parisella, M.L.; Cioni, D.; Lencioni, R.; Neri, E. Diagnostic and Prognostic Role of Late Gadolinium Enhancement in Cardiomyopathies. Eur. Heart J. Suppl. J. Eur. Soc. Cardiol. 2023, 25, C130–C136. [Google Scholar] [CrossRef]
  10. Bull, S.; White, S.K.; Piechnik, S.K.; Flett, A.S.; Ferreira, V.M.; Loudon, M.; Francis, J.M.; Karamitsos, T.D.; Prendergast, B.D.; Robson, M.D.; et al. Human Non-Contrast T1 Values and Correlation with Histology in Diffuse Fibrosis. Heart 2013, 99, 932–937. [Google Scholar] [CrossRef]
  11. Kwong, R.Y.; Kramer, C.M.; Chandrashekhar, Y. CMR Global Longitudinal Strain: A Better Tool for Unraveling the Links to Heart Failure Mortality. JACC Cardiovasc. Imaging 2018, 11, 1554–1555. [Google Scholar] [CrossRef] [PubMed]
  12. Cheng, S.; Wang, J.; Wang, Y.; Qi, L.; Li, F.; Liu, J.; Chen, J.; Fan, Y.; Xie, L. Longitudinal Assessment of Cardiac Parameters through MRI in Breast Cancer Patients Treated with Anti-HER2 Therapy. Eur. Radiol. Exp. 2023, 7, 22. [Google Scholar] [CrossRef] [PubMed]
  13. Houbois, C.P.; Nolan, M.; Somerset, E.; Shalmon, T.; Esmaeilzadeh, M.; Lamacie, M.M.; Amir, E.; Brezden-Masley, C.; Koch, C.A.; Thevakumaran, Y.; et al. Serial Cardiovascular Magnetic Resonance Strain Measurements to Identify Cardiotoxicity in Breast Cancer: Comparison With Echocardiography. JACC Cardiovasc. Imaging 2021, 14, 962–974. [Google Scholar] [CrossRef] [PubMed]
  14. Wang, R.; Zhou, Z.; Schoepf, U.J.; Varga-Szemes, A.; Strigenz, A.; Wang, H.; Liu, Y.; Xu, L. Monitoring of Anthracycline-Induced Myocardial Injury Using Serial Cardiac Magnetic Resonance: An Animal Study. Int. J. Cardiol. 2021, 328, 111–116. [Google Scholar] [CrossRef] [PubMed]
  15. Romano, S.; Judd, R.M.; Kim, R.J.; Kim, H.W.; Klem, I.; Heitner, J.F.; Shah, D.J.; Jue, J.; White, B.E.; Indorkar, R.; et al. Feature-Tracking Global Longitudinal Strain Predicts Death in a Multicenter Population of Patients With Ischemic and Nonischemic Dilated Cardiomyopathy Incremental to Ejection Fraction and Late Gadolinium Enhancement. JACC Cardiovasc. Imaging 2018, 11, 1419–1429. [Google Scholar] [CrossRef] [PubMed]
  16. Conte, E.; Mushtaq, S.; Muscogiuri, G.; Formenti, A.; Annoni, A.; Mancini, E.; Ricci, F.; Melotti, E.; Gigante, C.; Lorenza, Z.; et al. The Potential Role of Cardiac CT in the Evaluation of Patients With Known or Suspected Cardiomyopathy: From Traditional Indications to Novel Clinical Applications. Front. Cardiovasc. Med. 2021, 8, 709124. [Google Scholar] [CrossRef] [PubMed]
  17. Kim, N.Y.; Im, D.J.; Youn, J.-C.; Hong, Y.J.; Choi, B.W.; Kang, S.-M.; Lee, H.-J. Synthetic Extracellular Volume Fraction Derived Using Virtual Unenhanced Attenuation of Blood on Contrast-Enhanced Cardiac Dual-Energy CT in Nonischemic Cardiomyopathy. AJR Am. J. Roentgenol. 2022, 218, 454–461. [Google Scholar] [CrossRef] [PubMed]
  18. Nacif, M.S.; Kawel, N.; Lee, J.J.; Chen, X.; Yao, J.; Zavodni, A.; Sibley, C.T.; Lima, J.A.C.; Liu, S.; Bluemke, D.A. Interstitial Myocardial Fibrosis Assessed as Extracellular Volume Fraction with Low-Radiation-Dose Cardiac CT. Radiology 2012, 264, 876–883. [Google Scholar] [CrossRef] [PubMed]
  19. Dulf, E.-H.; Bledea, M.; Mocan, T.; Mocan, L. Automatic Detection of Colorectal Polyps Using Transfer Learning. Sensors 2021, 21, 5704. [Google Scholar] [CrossRef]
  20. Stoleru, C.-A.; Dulf, E.H.; Ciobanu, L. Automated Detection of Celiac Disease Using Machine Learning Algorithms. Sci. Rep. 2022, 12, 4071. [Google Scholar] [CrossRef]
  21. Danku, A.E.; Dulf, E.H.; Banut, R.P.; Silaghi, H.; Silaghi, C.A. Cancer Diagnosis With the Aid of Artificial Intelligence Modeling Tools. IEEE Access 2022, 10, 20816–20831. [Google Scholar] [CrossRef]
  22. Kramer, C.M.; Barkhausen, J.; Flamm, S.D.; Kim, R.J.; Nagel, E. Society for Cardiovascular Magnetic Resonance Board of Trustees Task Force on Standardized Protocols Standardized Cardiovascular Magnetic Resonance (CMR) Protocols 2013 Update. J. Cardiovasc. Magn. Reson. Off. J. Soc. Cardiovasc. Magn. Reson. 2013, 15, 91. [Google Scholar] [CrossRef]
  23. Popa, A.; Cionca, C.; Agoston, R.; Rusu, F.; Tarcau, B.M.; Negru, A.; Orzan, R.I.; Agoston-Coldea, L. The Role of Magnetic Resonance Imaging in Risk Stratification of Patients with Acute Myocarditis. Diagnostics 2024, 14, 1426. [Google Scholar] [CrossRef]
  24. Rogers, T.; Dabir, D.; Mahmoud, I.; Voigt, T.; Schaeffter, T.; Nagel, E.; Puntmann, V.O. Standardization of T1 Measurements with MOLLI in Differentiation between Health and Disease—The ConSept Study. J. Cardiovasc. Magn. Reson. 2013, 15, 78. [Google Scholar] [CrossRef] [PubMed]
  25. Dabir, D.; Child, N.; Kalra, A.; Rogers, T.; Gebker, R.; Jabbour, A.; Plein, S.; Yu, C.-Y.; Otton, J.; Kidambi, A.; et al. Reference Values for Healthy Human Myocardium Using a T1 Mapping Methodology: Results from the International T1 Multicenter Cardiovascular Magnetic Resonance Study. J. Cardiovasc. Magn. Reson. 2014, 16, 69. [Google Scholar] [CrossRef] [PubMed]
  26. Orzan, R.I.; Gligor, R.I.; Agoston, R.; Cionca, C.; Zlibut, A.; Pais, R.; Seicean, A.; Agoston-Coldea, L. JCDD|Free Full-Text|Metabolic Dysfunction Associated Liver Disease in Patients Undergoing Coronary Computed Tomography Angiography. J. Cardiovasc. Dev. Dis. 2024, 11, 77. [Google Scholar]
  27. Van Der Bijl, N.; Joemai, R.M.S.; Geleijns, J.; Bax, J.J.; Schuijf, J.D.; De Roos, A.; Kroft, L.J.M. Assessment of Agatston Coronary Artery Calcium Score Using Contrast-Enhanced CT Coronary Angiography. Am. J. Roentgenol. 2010, 195, 1299–1305. [Google Scholar] [CrossRef]
  28. Han, X. On Statistical Measures for Data Quality Evaluation. J. Geogr. Inf. Syst. 2020, 12, 178–187. Available online: https://www.scirp.org/journal/paperinformation?paperid=100819 (accessed on 1 August 2024). [CrossRef]
  29. Lopatin, Y. Heart Failure with Mid-Range Ejection Fraction and How to Treat It. Card. Fail. Rev. 2018, 4, 9–13. [Google Scholar] [CrossRef]
  30. Cove-Smith, L.; Woodhouse, N.; Hargreaves, A.; Kirk, J.; Smith, S.; Price, S.A.; Galvin, M.; Betts, C.J.; Brocklehurst, S.; Backen, A.; et al. An Integrated Characterization of Serological, Pathological, and Functional Events in Doxorubicin-Induced Cardiotoxicity. Toxicol. Sci. 2014, 140, 3–15. [Google Scholar] [CrossRef]
  31. Li, W.; Croce, K.; Steensma, D.P.; McDermott, D.F.; Ben-Yehuda, O.; Moslehi, J. Vascular and Metabolic Implications of Novel Targeted Cancer Therapies: Focus on Kinase Inhibitors. J. Am. Coll. Cardiol. 2015, 66, 1160–1178. [Google Scholar] [CrossRef] [PubMed]
  32. Jordan, J.H.; Todd, R.M.; Vasu, S.; Hundley, W.G. Cardiovascular Magnetic Resonance in the Oncology Patient. JACC Cardiovasc. Imaging 2018, 11, 1150–1172. [Google Scholar] [CrossRef] [PubMed]
  33. Scatteia, A.; Baritussio, A.; Bucciarelli-Ducci, C. Strain Imaging Using Cardiac Magnetic Resonance. Heart Fail. Rev. 2017, 22, 465–476. [Google Scholar] [CrossRef] [PubMed]
  34. Kramer, C.M. Role of Cardiac MR Imaging in Cardiomyopathies. J. Nucl. Med. 2015, 56, 39S–45S. [Google Scholar] [CrossRef] [PubMed]
  35. Modi, K.; Joppa, S.; Chen, K.-H.A.; Athwal, P.S.S.; Okasha, O.; Velangi, P.S.; Hooks, M.; Nijjar, P.S.; Blaes, A.H.; Shenoy, C. Myocardial Damage Assessed by Late Gadolinium Enhancement on Cardiovascular Magnetic Resonance Imaging in Cancer Patients Treated with Anthracyclines and/or Trastuzumab. Eur. Heart J. Cardiovasc. Imaging 2020, 22, 427–434. [Google Scholar] [CrossRef] [PubMed]
  36. Hong, Y.J.; Park, H.S.; Park, J.K.; Han, K.; Park, C.H.; Kim, T.K.; Yoo, S.J.; Lee, J.Y.; Kim, P.K.; Hur, J.; et al. Early Detection and Serial Monitoring of Anthracycline-Induced Cardiotoxicity Using T1-Mapping Cardiac Magnetic Resonance Imaging: An Animal Study. Sci. Rep. 2017, 7, 2663. [Google Scholar] [CrossRef] [PubMed]
  37. Tham, E.B.; Haykowsky, M.J.; Chow, K.; Spavor, M.; Kaneko, S.; Khoo, N.S.; Pagano, J.J.; Mackie, A.S.; Thompson, R.B. Diffuse Myocardial Fibrosis by T1-Mapping in Children with Subclinical Anthracycline Cardiotoxicity: Relationship to Exercise Capacity, Cumulative Dose and Remodeling. J. Cardiovasc. Magn. Reson. Off. J. Soc. Cardiovasc. Magn. Reson. 2013, 15, 48. [Google Scholar] [CrossRef] [PubMed]
  38. Jordan, J.H.; Vasu, S.; Morgan, T.M.; D’Agostino, R.B.; Meléndez, G.C.; Hamilton, C.A.; Arai, A.E.; Liu, S.; Liu, C.-Y.; Lima, J.A.C.; et al. Anthracycline-Associated T1 Mapping Characteristics Are Elevated Independent of the Presence of Cardiovascular Comorbidities in Cancer Survivors. Circ. Cardiovasc. Imaging 2016, 9, e004325. [Google Scholar] [CrossRef] [PubMed]
  39. Zhang, Y.; Zhang, X.; Wang, Y.; Hu, X.; Wang, B.; Yang, J.; Zhao, X.; Zhang, L. Relationship between Diffuse Fibrosis Assessed by CMR and Depressed Myocardial Strain in Different Stages of Heart Failure. Eur. J. Radiol. 2023, 164, 110848. [Google Scholar] [CrossRef]
  40. Cundari, G.; Galea, N.; Mergen, V.; Alkadhi, H.; Eberhard, M. Myocardial Extracellular Volume Quantification with Computed Tomography—Current Status and Future Outlook. Insights Imaging 2023, 14, 156. [Google Scholar] [CrossRef]
  41. Penso, M.; Babbaro, M.; Moccia, S.; Baggiano, A.; Carerj, M.L.; Guglielmo, M.; Fusini, L.; Mushtaq, S.; Andreini, D.; Pepi, M.; et al. A Deep-Learning Approach for Myocardial Fibrosis Detection in Early Contrast-Enhanced Cardiac CT Images. Front. Cardiovasc. Med. 2023, 10, 1151705. [Google Scholar] [CrossRef] [PubMed]
  42. Tatsugami, F.; Nakaura, T.; Yanagawa, M.; Fujita, S.; Kamagata, K.; Ito, R.; Kawamura, M.; Fushimi, Y.; Ueda, D.; Matsui, Y.; et al. Recent Advances in Artificial Intelligence for Cardiac CT: Enhancing Diagnosis and Prognosis Prediction. Diagn. Interv. Imaging 2023, 104, 521–528. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Study population was determined after applying inclusion and exclusion criteria (BC: breast cancer; HER2: human epidermal growth factor receptor 2; LVEF: left ventricle ejection fraction; CMR: cardiac magnetic resonance; CCT: cardiac computer tomography).
Figure 1. Study population was determined after applying inclusion and exclusion criteria (BC: breast cancer; HER2: human epidermal growth factor receptor 2; LVEF: left ventricle ejection fraction; CMR: cardiac magnetic resonance; CCT: cardiac computer tomography).
Jcm 13 04807 g001
Figure 2. Results generated by the algorithm when using different CT images of the same patient.
Figure 2. Results generated by the algorithm when using different CT images of the same patient.
Jcm 13 04807 g002
Figure 3. Results generated by the algorithm when using same CT image of patient 1 with different values for the number of iterations of the Snake algorithm.
Figure 3. Results generated by the algorithm when using same CT image of patient 1 with different values for the number of iterations of the Snake algorithm.
Jcm 13 04807 g003
Figure 4. Scatter plot of fibrosis mean (%) compared to LGE/LV mass (%) and LGE (g).
Figure 4. Scatter plot of fibrosis mean (%) compared to LGE/LV mass (%) and LGE (g).
Jcm 13 04807 g004
Figure 5. Scatter plot of fibrosis mean (%) compared to GLS (%).
Figure 5. Scatter plot of fibrosis mean (%) compared to GLS (%).
Jcm 13 04807 g005
Table 1. Summary of the clinical characteristics.
Table 1. Summary of the clinical characteristics.
Clinical CharacteristicsPatients (n = 8)
Age54.88 ± 12.05
Weight (kg)87.38 ± 15.14
Height (cm)164.38 ± 7.98
BMI (kg/m2)32.21 ± 4.3
Pulse rate (bpm)72.88 ± 10.93
SBP (mmHg)133.75 ± 13.02
DBP (mmHg)78.75 ± 13.82
Total cholesterol (mg/dL)203 ± 46.91
HDL-C (mg/dL)45.13 ± 8.43
LDL-C (mg/dl)125.38 ± 37.97
Tg (mg/dL)161.75 ± 71.21
Blood glucose (mmol/L)113.38 ± 11.87
Creatinine (mg/dL)0.76 ± 0.18
GFR (ml/min/1.73 m2)87.88 ± 22.15
Hypertension (%)50
Diabetes (%)12.5
Smoker status
Active (%)12.5
Non-smoker (%)75
Ex-smoker (%)12.5
Symptoms
Dyspnea (%)75
Chest pain (%) 25
Palpitations (%)37.5
ECG findings
Non-specific ST-T changes (%)62.5
Medication
Beta blocker (%)87.5
ACE inhibitors (%)75
Diuretics (%)62.5
Calcium blocker (%)0
Statins (%)50
Antiplatelets drug (%)25
NT-proBNP, (IQR) pg/mL785 (290–960)
BMI body mass index; SBP systolic blood pressure; DBP diastolic blood pressure; PP: pulse pressure; Tg triglycerides; GFR glomerular filtration rate; ECG electrocardiogram; ACE angiotensin-converting enzyme; The results are reported as mean ± SD, if not otherwise specified.
Table 2. CMR and CCT characteristics of the enrolled patients.
Table 2. CMR and CCT characteristics of the enrolled patients.
CMR
(n = 8)
CT
(n = 8)
p
LVEDVI, mean (SD), mL/m265.0 (11.6)60.4 (13.1)0.07
LVESVI, mean (SD), mL/m235.0 (6.4)30.9 (8.4)0.076
LVEF, mean (SD), %46.3 (3.7)48.8 (1.1)0.579
LV mass, g/m2 (SD) 71.8 (11.2)64 (15.3)0.027
GLS, mean (SD), %−12.25 (3.87)-NA
LGE+, n (%) 7 (87.2)-NA
Myocardial localization—septal/lateral/anterior/inferior/circumferential 3/1/2/2/0-NA
Myocardial pattern—subepicardial/nodal/midwall 2/5/1-NA
LV-LGE, mean (SD), g24.8 (12.3)-NA
LV-LGE/LV mass, %17.9 (8.9)-NA
Native T1 mapping, mean (SD), ms 1114 +/−51-NA
ECV, mean (SD), %31.2 (3.1)-NA
Pericardial effusion+, n (%) 3 (37.5)-NA
LVEDVI left ventricle end-diastolic volume index; LVESVI left ventricle end-systolic volume index; LVEF left ventricle ejection fraction; LV left ventricle; GLS global longitudinal strain; LGE late gadolinium enhancement; ECV extracellular volume.
Table 3. Percentages of fibrosis generated by the algorithm when different CT images of the same patient are chosen.
Table 3. Percentages of fibrosis generated by the algorithm when different CT images of the same patient are chosen.
Fibrosis Percentage as Returned by the Algorithm (%) with 3 Different Images of the Same Patient
Patient#1#2#3MeanSDCV (%)
18.99.28.798.960.211.93
25.777.576.516.620.9011.16
37.519.759.598.951.2511.40
413.8711.0810.5211.821.7912.39
56.539.179.48.371.5915.56
611.5310.411.0410.990.574.21
75.684.865.735.420.497.35
89.719.948.979.540.514.34
CV coefficient of variance; # number of image.
Table 4. Fibrosis percentage, as reported by the algorithm, after running it 3 times by the same observer, for the region of interest identified using a different number of iterations for the Snake algorithm.
Table 4. Fibrosis percentage, as reported by the algorithm, after running it 3 times by the same observer, for the region of interest identified using a different number of iterations for the Snake algorithm.
Fibrosis Percentage as Returned by the Algorithm (%)
Patient#1#2#3MeanCV (%)
18.99.098.828.941.27
25.775.815.815.800.33
37.517.797.47.572.17
413.8713.8513.6413.790.75
56.536.786.536.611.78
611.5311.5511.5311.540.08
75.685.115.755.515.20
89.719.839.719.750.58
CV coefficient of variance; # number of measurement.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gonciar, D.; Berciu, A.-G.; Dulf, E.-H.; Orzan, R.I.; Mocan, T.; Danku, A.E.; Lorenzovici, N.; Agoston-Coldea, L. Computer-Assisted Algorithm for Quantification of Fibrosis by Native Cardiac CT: A Pilot Study. J. Clin. Med. 2024, 13, 4807. https://doi.org/10.3390/jcm13164807

AMA Style

Gonciar D, Berciu A-G, Dulf E-H, Orzan RI, Mocan T, Danku AE, Lorenzovici N, Agoston-Coldea L. Computer-Assisted Algorithm for Quantification of Fibrosis by Native Cardiac CT: A Pilot Study. Journal of Clinical Medicine. 2024; 13(16):4807. https://doi.org/10.3390/jcm13164807

Chicago/Turabian Style

Gonciar, Diana, Alexandru-George Berciu, Eva-Henrietta Dulf, Rares Ilie Orzan, Teodora Mocan, Alex Ede Danku, Noemi Lorenzovici, and Lucia Agoston-Coldea. 2024. "Computer-Assisted Algorithm for Quantification of Fibrosis by Native Cardiac CT: A Pilot Study" Journal of Clinical Medicine 13, no. 16: 4807. https://doi.org/10.3390/jcm13164807

APA Style

Gonciar, D., Berciu, A. -G., Dulf, E. -H., Orzan, R. I., Mocan, T., Danku, A. E., Lorenzovici, N., & Agoston-Coldea, L. (2024). Computer-Assisted Algorithm for Quantification of Fibrosis by Native Cardiac CT: A Pilot Study. Journal of Clinical Medicine, 13(16), 4807. https://doi.org/10.3390/jcm13164807

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop