Skip to main content

ORIGINAL RESEARCH article

Front. Cardiovasc. Med., 11 November 2020
Sec. Cardiovascular Imaging

Neural-Network-Based Diagnosis Using 3-Dimensional Myocardial Architecture and Deformation: Demonstration for the Differentiation of Hypertrophic Cardiomyopathy

\nAlessandro Satriano
Alessandro Satriano1*Yarmaghan AfzalYarmaghan Afzal1Muhammad Sarim AfzalMuhammad Sarim Afzal1Ali Fatehi HassanabadAli Fatehi Hassanabad2Cody WuCody Wu1Steven DykstraSteven Dykstra2Jacqueline Flewitt,,Jacqueline Flewitt1,2,3Patricia FeuchterPatricia Feuchter1Rosa SandonatoRosa Sandonato1Bobak HeydariBobak Heydari1Naeem Merchant,,Naeem Merchant1,3,4Andrew G. Howarth,,Andrew G. Howarth1,2,3Carmen P. Lydell,,Carmen P. Lydell1,3,4Aneal KhanAneal Khan5Nowell M. FineNowell M. Fine3Russell Greiner,Russell Greiner6,7James A. White,,James A. White1,2,3
  • 1Stephenson Cardiac Imaging Center, Calgary, AB, Canada
  • 2Division of Cardiology, School of Medicine, University of Calgary, Calgary, AB, Canada
  • 3Libin Cardiovascular Institute of Alberta, Calgary, AB, Canada
  • 4Department of Diagnostic Imaging, University of Calgary, Calgary, AB, Canada
  • 5Department of Medical Genetics, University of Calgary, Calgary, AB, Canada
  • 6Department of Computing Science, University of Alberta, Edmonton, AB, Canada
  • 7Alberta Machine Learning Institute, Edmonton, AB, Canada

The diagnosis of cardiomyopathy states may benefit from machine-learning (ML) based approaches, particularly to distinguish those states with similar phenotypic characteristics. Three-dimensional myocardial deformation analysis (3D-MDA) has been validated to provide standardized descriptors of myocardial architecture and deformation, and may therefore offer appropriate features for the training of ML-based diagnostic tools. We aimed to assess the feasibility of automated disease diagnosis using a neural network trained using 3D-MDA to discriminate hypertrophic cardiomyopathy (HCM) from its mimic states: cardiac amyloidosis (CA), Anderson–Fabry disease (AFD), and hypertensive cardiomyopathy (HTNcm). 3D-MDA data from 163 patients (mean age 53.1 ± 14.8 years; 68 females) with left ventricular hypertrophy (LVH) of known etiology was provided. Source imaging data was from cardiac magnetic resonance (CMR). Clinical diagnoses were as follows: 85 HCM, 30 HTNcm, 30 AFD, and 18 CA. A fully-connected-layer feed-forward neural was trained to distinguish HCM vs. other mimic states. Diagnostic performance was compared to threshold-based assessments of volumetric and strain-based CMR markers, in addition to baseline clinical patient characteristics. Threshold-based measures provided modest performance, the greatest area under the curve (AUC) being 0.70. Global strain parameters exhibited reduced performance, with AUC under 0.64. A neural network trained exclusively from 3D-MDA data achieved an AUC of 0.94 (sensitivity 0.92, specificity 0.90) when performing the same task. This study demonstrates that ML-based diagnosis of cardiomyopathy states performed exclusively from 3D-MDA is feasible and can distinguish HCM from mimic disease states. These findings suggest strong potential for computer-assisted diagnosis in clinical practice.

Introduction

The application of machine-learning (ML) techniques for computer-assisted diagnosis is an anticipated milestone for cardiac imaging. However, robust approaches to the classification of cardiomyopathy states requires concurrent consideration of architectural features as well as regional patterns of deformation. This set of requirements poses challenges for the development of ML-assisted diagnosis in cardiovascular diagnostics.

Three-dimensional myocardial deformation analysis (3D-MDA) is a recently validated (19) image post-processing technique for the transformation of dynamic imaging data into a standardized, time-resolved 3D model of cardiac phenotype. Using foundational techniques transferrable across cardiac magnetic resonance (CMR) (1, 4, 68), computed tomography (CT) (10) and 3D echocardiography (3), this offers a common approach to convert source images into data required to describe chamber phenotypes using spatially-resolved measures of wall thickness and principal strain (PS), the latter an objective measure of tissue deformation in the direction resulting from locally-engaged myocardial fibers (1, 9). We hypothesized that a ML-based model could use this standardized phenotypic data to classify cardiomyopathy state. As a sentinel demonstration of this novel paradigm, we chose to assess diagnostic accuracy of a neural network trained exclusively from 3D-MDA data (wall architecture and deformation) to distinguish hypertrophic cardiomyopathy (HCM) from its commonly encountered mimic states: cardiac amyloidosis (CA), Anderson–Fabry disease (AFD), and hypertensive cardiomyopathy (HTNcm).

Materials and Methods

Study Population

Patients were identified from the Cardiovascular Imaging Registry of Canada (CIROC), a prospective initiative of the Libin Cardiovascular Institute at the University of Calgary (Calgary, Alberta) (NCT04367220). Data management is executed by commercial software (cardioDI™, Cohesic Inc., Calgary). This data repository was interrogated to identify patients with a confirmed etiology of left ventricular hypertrophy (LVH), defined as an indexed LV mass index ≥2 SD above age and sex-based reference values (11) and a maximal LV wall thickness ≥13 mm by linear measurement on at least 1 short-axis (SAX) cine view. Etiology of LVH was established by review of medical records and all available diagnostic testing. HCM diagnosis required either genotype confirmation of a pathologic sarcomere protein mutation, a 1st degree relative with the same, or a maximal wall thickness ≥15 mm, plus typical mid-wall patchy fibrosis by late gadolinium enhancement (LGE) imaging and no other identifiable cause. CA diagnosis required a typical diffuse pattern of LGE with confirmatory testing by either fat pad biopsy (for light chain amyloidosis) or 99m-technetium-pyrophosphate scintigraphy imaging (for transthyretin variant amyloidosis). AFD diagnosis required positive genotype for a pathologic mutation of the alpha galactosidase encoding gene. Finally, HTNcm diagnosis required a history of hypertension for ≥10 years, no family history of HCM or AFD, and no LGE-based features suggestive of an alternate myocardial disease state.

All CMR studies were de-identified and transferred to a study server for blinded analysis. The study design was approved by the Conjoint Health Research Ethics Board at the University of Calgary and all subjects provided written informed consent. All research activities were performed in accordance with the Declaration of Helsinki.

MRI Image Acquisition

All CMR studies were performed using 3T scanners (Prisma or Skyra, Siemens Healthineers, Erlangen, Germany). A standardized imaging protocol was performed in all patients inclusive of both cine and LGE imaging. Cine imaging was performed in sequential SAX and standard long-axis (LAX) imaging planes at end-expiration using a steady-state free precession pulse sequence. Typical imaging parameters were: repetition time 3.1 ms, 6 lines per segment, echo time 1.3 ms, flip angle 45°, field of view 276 × 360 mm2, matrix 156 × 192, slice thickness 6 mm, gap 2 mm, parallel imaging factor of 2, reconstructed to 30 cardiac phases. LGE imaging was performed 7-10 min following intravenous administration of 0.1 mmol/kg gadolinium chelate (Gadovist, Bayer Inc., Canada) in matched slice orientations using a standard inversion-recovery gradient pulse sequence with phase sensitive image reconstruction (12).

Conventional Volumetric Chamber Analysis

All analyses were performed by trained core laboratory personnel blinded to clinical data. Each type of image analysis was performed independently and in random order. LV volumes and mass were calculated using commercially available software (cvi42, Circle Cardiovascular Inc., Calgary, Canada) from sequential SAX cine images using semi-automated tracing of the endocardial and epicardial contours to derive LV end-diastolic volume (LVEDV), LV end-systolic volume (LVESV), LV ejection fraction (LVEF), and mass. Papillary muscles were considered part of the LV mass. Chamber volumes and mass were indexed to body surface area (BSA) according to the Mostellar formula.

3D Myocardial Deformation Analysis

3D-MDA of the LV was performed using validated in-house software, leveraging previously established technologies for image registration, temporal tracking, mesh generation, and computation of deformation (14, 6, 13, 14). Briefly, a 4D displacement field is generated using optical-flow-based feature tracking (7, 8) of all pixels from co-registered multi-planar 2D cines (1, 4, 6, 7), automatically adjusting for breath-hold motion (15). This field is used to deform an end-diastolic 3D mesh model of the LV (1, 3, 4, 6), defined by ~500 nodal elements present on each of the endocardial and epicardial layers. This results in a cross-subject (and cross-modality) anatomically coherent model throughout 30 phases of the cardiac cycle (Figure 1) (1, 4, 6). Using this standardized dynamic mesh, peak-systolic strain amplitude, time to peak-systolic strain, maximum systolic and diastolic strain rates, and end-diastolic thickness are estimated from pre-defined hexahedral structured mesh elements (16, 17). For this study, corresponding analyses were expressed as (i) global mean values and (ii) segmental mean values according to the American Heart Association (AHA) standard (18). All variables were calculated for conventional (longitudinal, circumferential, and radial) and principal directions (minimum and maximum) of deformation, with regional wall thickness calculated in the radial direction (1, 2, 9). Data were reported for subendocardial, subepicardial, and transmural layers, as appropriate.

FIGURE 1
www.frontiersin.org

Figure 1. Illustration of 3D-MDA workflow. Top left: Endocardial and epicardial contour tracing performed at end-diastole on long-axis cine images with automated registration to short-axis images and 3D mesh generation. Bottom left: A 4D displacement field is generated (1, 4, 68, 15), and used to deform the end-diastolic phase 3D LV mesh (1, 3, 4, 6) throughout the cardiac cycle (1, 4, 6). Right: From this dynamic mesh model strain, conventional and principal strain quantification is performed (16, 17). In the example, principal strain amplitude and direction lines are shown for a patient with confirmed hypertrophic cardiomyopathy (HCM).

Figure 2 provides atlas-based representations of 3D-MDA features that were used to predict disease etiology, simultaneously presenting architectural and deformation-based features used in the classifier. Atlases have been rendered at peak-systole frame with color-coding of peak strain amplitude (minimum principal strain) along directional lines of deformation. The latter aims to provide a visual representation of the intrinsic (local) direction of dominant tissue deformation as described by minimum principal strain (1), this being a resulting effect of myofibril engagement between the endocardial and epicardial layers (9).

FIGURE 2
www.frontiersin.org

Figure 2. Average minimum principal strain and direction average atlases for patients with hypertrophic cardiomyopathy (HCM: N = 85), cardiac amyloidosis (CA: N = 18), Anderson–Fabry Disease (AFD: N = 30), and Hypertensive cardiomyopathy (HTNcm: N = 30) (9).

Regional-Pattern-Based Classifier

A feed-forward neural network model was used as classifier. We aimed to capture composite features indicative of HCM phenotype vs. any other disease phenocopy state (CA, AFD or HTNcm) using 3D-MDA data. Our goal was a ML-based model that could use 917 segmental architectural and deformation-based features. These features were provided as AHA segmental, layer-specific (i.e., subepicardial, subendocardial, and transmural) systolic amplitude, time to systolic peak, peak-systolic rate, peak-diastolic strain rate, and wall thickness. Each AHA segmental feature was obtained by averaging values across all nodes within the corresponding AHA segment. To develop the ML-based model we used a neural network taking all 917 features as inputs. Two fully-connected hidden layers (30 and 5 neurons, respectively, each using a hyperbolic tangent sigmoid transformation) were used. A final read-out single-perceptron layer, followed by a log-sigmoid activation function, provided a real-valued output between zero and one to predict the presence of HCM (assigned to HCM when >0.5).

The choice of architecture (i.e., number of neurons in each layer, and number of hidden layers) was based on using, for each layer, the square root of the number of features from the previous layer, in order to establish a proof of concept model. Following this criterion to increasingly reduce the number of features throughout the network, only two hidden layers were deemed necessary. Our learner used a cross-entropy loss function as a reference point for the backward pass, and trained for 100 epochs.

The full pipeline from contouring to meshing, tracking, strain analysis, and neural network building can be replicated following the previously described 3D-MDA methodology (1), and building a fully-connected network according to the above specified hyperparameters for size of the layers, network depth, and used activation functions. Image acquisition parameters are provided within the “MRI Image acquisition” section.

While preferable to use one set of instances for training and a holdout subset to estimate the accuracy of the learned model, given our small sample size, we estimated diagnostic performance using a 5-fold cross-validation procedure (19). This involved using 4/5-ths of the patients to develop a model whose performance was evaluated on the remaining 1/5-th of patients, which were not previously seen by the training process. This procedure was repeated five times, rotating the 4/5-th training portion and 1/5-th testing portion. It is worth noting that, on each fold, feature selection and parameter estimation (training) were performed from scratch on the first 4/5-ths, thus relegating the last 1/5-th to an external holdout at each iteration, and reducing the chance of overfitting.

Performance was evaluated as area under the curve (AUC) of the receiver operating characteristic (ROC) curve as well as respective measures of sensitivity, specificity, positive and negative predictive values, and accuracy. These values were reported as mean with standard deviation range derived from the 5-fold cross-validation procedure (19). Training and validation procedures were implemented via Matlab Statistics and Machine Learning Toolbox (R2019b, The Mathworks, Inc., Natick, MA, USA).

Statistical Analysis

Categorical variables were presented as counts with percentages, while continuous variables were expressed as mean with standard deviation. Categorical variables were compared using the Fisher's exact test, and comparisons for continuous data were performed with the 2-sample t-test or Mann-Whitney U test, where appropriate. A two-sided p-value of <0.05 was considered statistically significant. Multiple comparisons between patients stratified by phenocopy disease state were performed by ANOVA. All statistical analysis was performed using the Matlab Statistics and Machine Learning Toolbox (R2019b, The Mathworks, Inc., Natick, MA, USA).

Results

Study Population

A total of 163 patients with LVH of known etiology were studied (85 HCM, 30 AFD, 30 HTNcm, 18 CA). Patient characteristics are provided in Table 1. The mean age was 52.7 ± 14.7 years, with 67 females (41%). No significant differences in baseline clinical characteristics were observed between HCM and non-HCM patients, with the exception of age and BSA. However, neither of these two variables provided significant predictive value for the classification of HCM (AUC = 0.63 for each). The following count was found for the HCM sub-phenotypes: (i) isolated basal-septal: 57 (67%), (ii) reverse septal curvature: 12 (14%), (iii) apical: 11 (13%), and (iv) other sub-phenotypes: 5 (6%).

TABLE 1
www.frontiersin.org

Table 1. Clinical characteristics of all patients with left ventricular hypertrophy (LVH), patients with LVH and confirmed etiology of hypertrophic cardiomyopathy (LVH+ HCM+) and patients with LVH and a mimicking cardiomyopathy state (LVH+ HCM–).

Disease Classification Using Conventional Volumetric Cine Imaging Markers

Of all conventionally reported CMR characteristics (Table 2), significant differences between HCM and non-HCM etiology were identified for LVEDVI (81.0 ± 15.1 mL/m2 vs. 75.5 ± 20.4 mL/m2, respectively— p < 0.01), LVEF (70.0 ± 7.5% vs. 65.4 ± 10.7%, respectively— p < 0.05), LVMI (84.2 ± 26.3 g/m2 vs. 75.9 ± 21.8 g/m2, respectively— p < 0.05), RVESVI (28.1 ± 8.8 mL/m2 vs. 35.6 ± 17.4 mL/m2, respectively— p < 0.01), RVEF (62.8 ± 7.5 ± vs. 57.7% ± 9.7%, respectively— p < 0.001), and LAEDVI (103.7 ± 38.5 mL/m2 vs. 78.0 ± 26.3 mL/m2, respectively —p < 0.0001). The highest AUC achieved of all such variables was for LAEDVI, reaching 0.70 (0.54— 0.90).

TABLE 2
www.frontiersin.org

Table 2. CMR-based volumetric measurements of all patients with left ventricular hypertrophy (LVH), patients with LVH and confirmed etiology of hypertrophic cardiomyopathy (LVH+ HCM+) and patients with LVH and a mimicking cardiomyopathy state (LVH+ HCM–).

Disease Classification Using Global Strain Amplitude and Time to Peak Systolic Strain

To establish a baseline, we studied the capacity of conventional strain-based measures to classify disease. Global peak-systolic strain amplitude values of the study population are presented in Table 3. Of conventional, axis-dependent measures of strain, global circumferential strain amplitude measured at the subendocardial layer provided the highest AUC for the classification of HCM vs. non-HCM etiology with an AUC of 0.62 (0.54—0.71). Global minimum principal strain measured at the subendocardial layer showed similar performance with an AUC of 0.64 (0.52—0.76), whereas transmurally the AUC was documented at 0.60 (0.49—0.71) in the radial direction.

TABLE 3
www.frontiersin.org

Table 3. Global peak systolic strain amplitude and end-diastolic wall thickness of all patients with left ventricular hypertrophy (LVH), patients with LVH and confirmed etiology of hypertrophic cardiomyopathy (LVH+HCM+) and patients with LVH and a mimicking cardiomyopathy state (LVH+ HCM–).

Time to peak-systolic strain amplitude data (indexed by duration of the cardiac cycle) is reported in Table 4. In the circumferential direction, HCM and non-HCM etiology were significantly different in strain timing, both when measured subendocardially (41.4 ± 5.1% vs. 43.2 ± 5.9%, respectively— p < 0.05) and transmurally (40.6 ± 5.0 ± 42.7 ± 5.6, respectively— p < 0.01). The latter achieved the highest timing-related AUC of 0.62 (0.53—0.72).

TABLE 4
www.frontiersin.org

Table 4. Global time to peak systolic strain (% of the cardiac cycle length) of all patients with left ventricular hypertrophy (LVH), patients with LVH and confirmed etiology of hypertrophic cardiomyopathy (LVH+HCM+) and patients with LVH and a mimicking cardiomyopathy state (LVH+HCM–).

Disease Classification Using Regional Wall Thickness

The sole use of architectural features to classify disease was explored, this considering segmental end-diastolic mean wall thickness derived from each 3D-MDA mesh-based model. This led to an AUC of 0.75 (0.71–0.80) for the discrimination of HCM vs. non-HCM etiology (Table 3).

Neural-Net-Based Disease Classification Using 3D-MDA-Derived Data

Finally, neural network-based classification performance of 3D-MDA data was assessed using all available architectural and deformation features, inclusive of regional measures of wall thickness, strain amplitude and time to peak-systolic strain. A neural network-based model, developed from all AHA segmental architectural and deformation-based features, demonstrated diagnostic performance with a mean AUC (ROC) of 0.94 (0.89–0.99) for the classification of HCM etiology, as shown in Figure 3. Mean sensitivity and specificity for HCM etiology was 0.92 (0.85–0.98) and 0.90 (0.83–0.97), respectively, with the mean accuracy being 0.91 (0.87–0.95). Mean positive and negative predictive values were 0.91 (0.85–0.97) and 0.91 (0.85−0.98), respectively.

FIGURE 3
www.frontiersin.org

Figure 3. Receiver operating characteristic (ROC) curves for the classification of patients with left ventricular hypertrophy (LVH) due to hypertrophic cardiomyopathy (LVH+HCM+) obtained during the 5-fold cross validation. The heavy black line indicates the median ROC curve.

Discussion

This study demonstrates capacity for ML-based methodologies to execute disease classification tasks, when provided a highly standardized set of myocardial architectural and deformation features, in this case through 3D-MDA of routine non-contrast cine CMR.

To our knowledge, this is the first study in the literature demonstrating feasibility of neural networks to classify the diagnosis of a hypertrophic cardiomyopathy state using 3D measures of myocardial architecture and deformation. Our neural network was trained using over 900 phenotypic features that were generated by 3D-MDA, this providing intrinsic temporo-spatial data registration across individual subjects. Using this structured data, we were able to deliver good performance using a two-layer neural network, achieving high discriminatory performance (AUC 0.94) for classifying HCM from its mimic cardiomyopathy states.

A single prior study has evaluated ML-based classification of HCM using 2D-based measures of deformation, these derived from speckle-tracking echocardiography (5). In this study, reported by Narula et al. (5), a classifier was trained to differentiate HCM from the presumptive diagnosis of athletic heart, achieving an AUC of 0.795. In contrast to this study, ours examined 3D features of architecture and deformation among clinically confirmed pathologic causes of hypertrophy. Furthermore, we compared model performance to the performance based on other clinically available markers, inclusive of clinical referral data, non-strain based volumetric CMR data, global strain-based markers, and 3D-MDA based measures modeled without ML techniques. We observed modest performance for conventional global strain-based measures to differentiate HCM from common phenocopy states. This finding is similar to a prior report by Neisius et al. (20), who attempted to distinguish HCM from HTNcm subjects and identified similar performance. By contrast, we demonstrated substantial gains in diagnostic accuracy for ML-based algorithms trained from 3D-MDA.

Dawes et al. (21) has previously described disease classification from 3D modeled descriptors of deformation for the right ventricle (RV). In this study, they elegantly demonstrated an atlas-based approach to perform RV segmentation from SAX cine images followed by modeling of chamber deformation (21). This allowed for construction of a dynamic 3D mesh model to derive RV deformation features in conventional axis-dependent directions of deformation. Following principal component analysis for feature reduction, this framing a regression-based approach, they trained models to identify patients at elevated risk of clinical outcomes, demonstrating superior performance to conventional CMR-based markers. This study highlighted a potential role for 3D myocardial deformation analyses to provide relevant information for ML-based tool development. The same research group reinforced this concept in predicting survival among a cohort of 302 patients with pulmonary hypertension (22).

While deformation is considered central to the diagnosis of cardiomyopathies, Swift et al. (23) explored an alternative approach to perform feature extraction directly from static (single phase) 2D cine CMR images. In this work a tensor-based machine-learned model was developed to evaluate pixel-based features from SAX and 4-chamber images, demonstrating the ability to classify presence vs. absence of pulmonary hypertension with good accuracy. However, this model did not perform significantly better than conventional CMR-based markers. How such methods will facilitate the differentiation of similar diseases remains uncertain at this time.

Our observation that regional deformation features provide incremental value to global architectural and global functional features in classifying disease etiology is consistent with emerging evidence that such measures are intimately coupled to underlying tissue characteristics (24). Strain-based markers have shown strong regional correlation with underlying markers of fibrosis (i.e., interstitial expansion) among patients with HCM (2, 2528) and CA (13, 29), as well as in non-hypertrophic states such as ischemic (3032) and dilated cardiomyopathy (33). Accordingly, capacity exists for regionally encoded markers of deformation to provide unique insights into underlying myocardial tissue health, and in-turn provide relevant information for the classification of disease etiology.

The use of routine strain-based markers to assist in the discrimination of hypertrophic disease states has been explored using both echocardiography (3436), and CMR-based (37) techniques. Williams et al. demonstrated that 2D strain analysis of cine CMR data shows greater relative preservation of apical global longitudinal strain in CA vs. HCM and AFD (37), and this pattern of apical sparing has similarly been described across echocardiography-based studies (3840). Reduction in global longitudinal strain have also been shown in biopsy proven HCM relative to hypertensive hypertrophy (36). While these studies have aimed to deliver simple and practical “rules of thumb” for identifying patients at greater likelihood of a specific disease, incremental value is inherently provided through ML-based modeling of raw deformation data, allowing for non-linear relationships to be identified across a larger number of features.

Study Limitations

Our study has several limitations. Our training cohort was modest in size and based on single-center data. As a consequence, we decided to contain model complexity by constructing it with only fully-connected hidden layers. Accordingly, our model may benefit from training dataset expansion and requires external validation in a multi-center clinical setting. Our evaluation was based on 5-fold cross-validation, as a separate hold-out validation cohort was not available. To address this, for each of the 5 cross-validations, training was restricted to the first 4/5-ths of the population followed by external validation on a 1/5-th hold-out population, this being an effective strategy to reduce the chance of overfitting. However, generalizability of our model to other patient populations (i.e., from other centers) has not been assessed. We chose to describe each patient using only segmentally-coded data from 3D-MDA in our neural network, despite having access to features derived for each hexahedral element of the mesh model. This decision was based on the size of our available population and practical application of the AHA segmental model as a feature reduction strategy. Consistent with clinical practice, we did not mandate genotype positivity for establishing a clinical diagnosis of HCM, rather also accepted phenotype positive patients if they had a first-degree family member with positive genotype or had a typical LGE-based pattern of fibrosis with no alternative explanatory disease. While potentially excluding mild HCM phenotypes without genotypic or LGE criteria, we aimed to maximize robustness of diagnostic classification for all LVH etiologies. The size of our neural-network was small and we used a simple fully-connected architecture to drastically reduce the number of features from each layer to the next one. This approach was taken to address the modest size of the available population and limit the number of parameters to be trained within the network. With exposure to larger populations, deeper networks with less aggressive feature reduction between layers are possibly applicable, and this is a priority for future studies. However, despite such limitations, strong discriminatory performance was achieved to distinguish HCM from non-HCM disease states in this study. We acknowledge that a significant difference in age and sex was observed between the respective disease cohorts, and that the influence of this on biomechanical disease profiles remains unknown. How training cohort demographics influence generalizability of ML-based diagnostic tools is critically important to be address in future studies. Finally, in anticipation of reductions in image quality, we a-priori excluded patients in atrial fibrillation. Accordingly, we cannot confirm performance of the described ML-based model for patients not in sinus rhythm.

Conclusions

We have demonstrated feasibility and acceptable performance of a neural network-based approach for the automated discrimination of HCM vs. its known phenocopy states. This unique approach, leveraging standardized phenotypic data from 3D-MDA, offers expanded potential for ML-assisted diagnostics and justifies broader investigation for other disease phenotypes.

Data Availability Statement

The datasets presented in this article are not readily available because the data underlying this article will be shared on reasonable request. Requests to access the datasets should be directed to James A. White, jawhit@ucalgary.ca.

Ethics Statement

The studies involving human participants were reviewed and approved by Conjoint Health Research Ethics Board at the University of Calgary. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.

Author Contributions

AS and JW contributed to conception and design of the study. AS contributed software and model development. SD and RG assisted with model development. YA, MS, and CW made measurements. PF assisted with image acquisition. RS contributed with patient recruitment. JF and JW acquired ethics approval. AS performed the statistical analysis, model assessment, and drafted the manuscript. AF, BH, NM, AH, CL, SD, AK, RG, and JW revised the manuscript critically for important intellectual content. All authors provided approval for publication of the manuscript.

Conflict of Interest

JF, AH, and JW are shareholders of Cohesic, Inc. AK is a member of the Canadian Fabry Disease Initiative (CFDI). AK and the CFDI have received sponsorship from any combination of Sanofi Genzyme, Takeda, and Amicus as research grants. The data contributed to this study were not collected as part of a clinical trial, under any form of sponsorship nor part of government or commercially funded enterprise, and utilized for academic purposes only.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors thank Dr. Gianni Pedrizzetti for his assistance with the graphical representation of principal directions.

References

1. Satriano A, Heydari B, Narous M, Exner DV, Mikami Y, Attwood MM, et al. Clinical feasibility and validation of 3D principal strain analysis from cine MRI: comparison to 2D strain by MRI and 3D speckle tracking echocardiography. Int J Cardiovasc Imaging. (2017) 33:1979–92. doi: 10.1007/s10554-017-1199-7

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Satriano A, Heydari B, Guron N, Fenwick K, Cheung M, Mikami Y, et al. 3-Dimensional regional and global strain abnormalities in hypertrophic cardiomyopathy. Int J Cardiovasc Imaging. (2019) doi: 10.1007/s10554-019-01631-8

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Satriano A, Pournazari P, Hirani N, Helmersen D, Thakrar M, Weatherald J, et al. Characterization of right ventricular deformation in pulmonary arterial hypertension using three-dimensional principal strain analysis. J Am Soc Echocardiogr. (2019) 32:385–93. doi: 10.1016/j.echo.2018.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Haber I, Metaxas DN, Axel L. Three-dimensional motion reconstruction and analysis of the right ventricle using tagged MRI. Med Image Anal. (2000) 4:335–55. doi: 10.1016/S1361-8415(00)00028-1

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Narula S, Shameer K, Salem Omar AM, Dudley JT, Sengupta PP. Machine-learning algorithms to automate morphological and functional assessments in 2D echocardiography. J Am Coll Cardiol. (2016) 68:2287–95. doi: 10.1016/j.jacc.2016.08.062

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Young A, Axel L. Non-rigid heart wall motion using MR tagging. In: Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition, Champaign, IL (1992). p. 399–404. doi: 10.1109/CVPR.1992.223158

CrossRef Full Text | Google Scholar

7. Oregan DP, Shi W, Ariff B, Baksi AJ, Durighel G, Rueckert D, et al. Remodeling after acute myocardial infarction: mapping ventricular dilatation using three dimensional CMR image registration. J Cardiovasc Magn Reson. (2012) 14:1–9. doi: 10.1186/1532-429X-14-41

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Xu C, Pilla J, Isaac G, Gorman J, Blom A, Gorman R, et al. Deformation analysis of 3D tagged cardiac images using an optical flow method. J Cardiovasc Magn Reson. (2010) 12:1–14. doi: 10.1186/1532-429X-12-19

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Pedrizzetti G, Sengupta S, Caracciolo G, Park CS, Amaki M, Goliasch G, et al. Three-dimensional principal strain analysis for characterizing subclinical changes in left ventricular function. J Am Soc Echocardiogr. (2014) 27:1041–50.e1. doi: 10.1016/j.echo.2014.05.014

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Satriano A, Guenther Z, White JA, Merchant N, Di Martino ES, Al-Qoofi F, et al. Three-dimensional thoracic aorta principal strain analysis from routine ECG-gated computerized tomography: feasibility in patients undergoing transcatheter aortic valve replacement. BMC Cardiovasc Disord. (2018) 18:1–11. doi: 10.1186/s12872-018-0818-0

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Kawel-Boehm N, Maceira A, Valsangiacomo-Buechel ER, Vogel-Claussen J, Turkbey EB, Williams R, et al. Normal values for cardiovascular magnetic resonance in adults and children. J Cardiovasc Magn Reson. (2015) 17:1–33. doi: 10.1186/s12968-015-0111-7

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Schieda N, Blaichman JI, Costa AF, Glikstein R, Hurrell C, James M, et al. Gadolinium-based contrast agents in kidney disease: a comprehensive review and clinical practice guideline issued by the canadian association of radiologists. Can J Kidney Heal Dis. (2018) 5:205435811877857. doi: 10.1177/2054358118778573

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Avitzur N, Satriano A, Afzal M, Narous M, Mikami Y, Hansen R, et al. 3D myocardial deformation analysis from cine MRI as a marker of amyloid protein burden in cardiac amyloidosis: validation vs. T1 mapping. Int J Cardiovasc Imaging. (2018) 34:1937–1946. doi: 10.1007/s10554-018-1410-5

CrossRef Full Text | Google Scholar

14. Todd A, Satriano A, Fenwick K, Merchant N, Lydell CP, Howarth AG, et al. Intra-thoracic adiposity is associated with impaired contractile function in patients with coronary artery disease: a cardiovascular magnetic resonance imaging study. Int J Cardiovasc Imaging. (2019) 35:121–31. doi: 10.1007/s10554-018-1430-1

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Carminati MC, Maffessanti F, Caiani EG. Nearly automated motion artifacts correction between multi breath-hold short-axis and long-axis cine CMR images. Comput Biol Med. (2014) 46:42–50. doi: 10.1016/j.compbiomed.2013.12.013

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Hallquist JO. LS-Dyna Theory Manual. Livermore, CA: Livermore Software Technology Corporation (LSTC) (2006). Available online at: http://www.lstc.com/pdf/ls-dyna_theory_manual_2006.pdf

17. Bonet J, Wood RD. Nonlinear Continuum Mechanics for Finite Element Analysis. 2nd ed. Cambridge: Cambridge University Press (2008). doi: 10.1017/CBO9780511755446

CrossRef Full Text | Google Scholar

18. Cerqueira MD, Weissman NJ, Dilsizian V, Jacobs AK, Kaul S, Laskey WK, et al. Standardized myocardial segmentation and nomenclature for tomographic imaging of the heart. A statement for healthcare professionals from the Cardiac Imaging Committee of the Council on Clinical Cardiology of the American Heart Association. Int J Cardiovasc Imaging. (2002) 18:539–42. doi: 10.1067/mje.2002.123374

PubMed Abstract | CrossRef Full Text | Google Scholar

19. James G, Witten D, Hastie T, Tibshirani R. An Introduction to Statistical Learning. New York, NY: Springer New York (2013). doi: 10.1007/978-1-4614-7138-7

CrossRef Full Text | Google Scholar

20. Neisius U, Myerson L, Fahmy AS, Nakamori S, El-Rewaidy H, Joshi G, et al. Cardiovascular magnetic resonance feature tracking strain analysis for discrimination between hypertensive heart disease and hypertrophic cardiomyopathy. PLoS ONE. (2019) 14:e0221061. doi: 10.1371/journal.pone.0221061

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Dawes TJW, De Marvao A, Shi W, Fletcher T, Watson GMJ, Wharton J, et al. Machine learning of threedimensional right ventricular motion enables outcome prediction in pulmonary hypertension: a cardiac MR imaging study. Radiology. (2017) 283:381–90. doi: 10.1148/radiol.2016161315

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Bello GA, Dawes TJW, Duan J, Biffi C, de Marvao A, Howard LSGE, et al. Deep-learning cardiac motion analysis for human survival prediction. Nat Mach Intell. (2019) 1:95–104. doi: 10.1038/s42256-019-0019-2

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Swift AJ, Lu H, Uthoff J, Garg P, Cogliano M, Taylor J, et al. A machine learning cardiac magnetic resonance approach to extract disease features and automate pulmonary arterial hypertension diagnosis. Eur Heart J Cardiovasc Imaging. (2020) 44:1–10. doi: 10.1093/ehjci/jeaa001

CrossRef Full Text | Google Scholar

24. Moharram MA, Lamberts RR, Whalley G, Williams MJA, Coffey S. Myocardial tissue characterisation using echocardiographic deformation imaging. Cardiovasc Ultrasound. (2019) 17:1–11. doi: 10.1186/s12947-019-0176-9

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Haland TF, Almaas VM, Hasselberg NE, Saberniak J, Leren IS, Hopp E. Strain echocardiography is related to fibrosis and ventricular arrhythmias in hypertrophic cardiomyopathy. Eur Heart J Cardiovasc Imaging. (2016) 17:613–21. doi: 10.1093/ehjci/jew005

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Popović ZB, Kwon DH, Mishra M, Buakhamsri A, Greenberg NL, Thamilarasan M, et al. Association between regional ventricular function and myocardial fibrosis in hypertrophic cardiomyopathy assessed by speckle tracking echocardiography and delayed hyperenhancement magnetic resonance imaging. J Am Soc Echocardiogr. (2008) 21:1299–305. doi: 10.1016/j.echo.2008.09.011

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Delgado V, Marsan NA. Global and regional longitudinal strain assessment in hypertrophic cardiomyopathy. Circ Cardiovasc Imaging. (2019) 12:1–3. doi: 10.1093/ehjci/jeu274

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Pagourelias ED, Mirea O, Duchenne J, Van Cleemput J, Delforge M, Bogaert J, et al. Echo parameters for differential diagnosis in cardiac amyloidosis: a head-to-head comparison of deformation and nondeformation parameters. Circ Cardiovasc Imaging. (2017) 10:1–11. doi: 10.1161/CIRCIMAGING.116.005588

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Gotschy A, Von Deuster C, Van Gorkum RJH, Gastl M, Vintschger E, Schwotzer R, et al. Characterizing cardiac involvement in amyloidosis using cardiovascular magnetic resonance diffusion tensor imaging. J Cardiovasc Magn Reson. (2019) 21:1–9. doi: 10.1186/s12968-019-0563-2

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Di Franco A, Kim J, Rodriguez-Diego S, Khalique O, Siden JY, Goldburg SR, et al. Multiplanar strain quantification for assessment of right ventricular dysfunction and non-ischemic fibrosis among patients with ischemic mitral regurgitation. PLoS ONE. (2017) 12:e0185657. doi: 10.1371/journal.pone.0185657

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Erley J, Genovese D, Tapaskar N, Alvi N, Rashedi N, Besser SA, et al. Echocardiography and cardiovascular magnetic resonance based evaluation of myocardial strain and relationship with late gadolinium enhancement. J Cardiovasc Magn Reson. (2019) 21:1–11. doi: 10.1186/s12968-019-0559-y

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Stathogiannis K, Mor-Avi V, Rashedi N, Lang RM, Patel AR. Regional myocardial strain by cardiac magnetic resonance feature tracking for detection of scar in ischemic heart disease. Magn Reson Imaging. (2020) 68:190–196. doi: 10.1016/j.mri.2020.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Taylor RJ, Umar F, Lin ELS, Ahmed A, Moody WE, Mazur W, et al. Mechanical effects of left ventricular midwall fibrosis in non-ischemic cardiomyopathy. J Cardiovasc Magn Reson. (2016) 18:1–8. doi: 10.1186/s12968-015-0221-2

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Sun JP, Stewart WJ, Yang XS, Donnell RO, Leon AR, Felner JM, et al. Differentiation of hypertrophic cardiomyopathy and cardiac amyloidosis from other causes of ventricular wall thickening by two-dimensional strain imaging echocardiography. Am J Cardiol. (2009) 103:411–415. doi: 10.1016/j.amjcard.2008.09.102

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Williams LK, Frenneaux MP, Steeds RP. Echocardiography in hypertrophic cardiomyopathy diagnosis, prognosis, and role in management. Eur J Echocardiogr. (2009) 10:9–14. doi: 10.1093/ejechocard/jep157

CrossRef Full Text | Google Scholar

36. Kato TS, Noda A, Izawa H, Yamada A, Obata K, Nagata K, et al. Discrimination of nonobstructive hypertrophic cardiomyopathy from hypertensive left ventricular hypertrophy on the basis of strain rate imaging by tissue Doppler ultrasonography. Circulation. (2004) 110:3808–3814. doi: 10.1161/01.CIR.0000150334.69355.00

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Williams LK, Forero JF, Popovic ZB, Phelan D, Delgado D, Rakowski H, et al. Patterns of CMR measured longitudinal strain and its association with late gadolinium enhancement in patients with cardiac amyloidosis and its mimics. J Cardiovasc Magn Reson. (2017) 19:1–10. doi: 10.1186/s12968-017-0376-0

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Nardozza M, Chiodi E, Mele D. Left ventricle relative apical sparing in cardiac amyloidosis. J Cardiovasc Echogr. (2017) 27:141–142. doi: 10.4103/jcecho.jcecho_22_17

CrossRef Full Text | Google Scholar

39. Phelan D, Collier P, Thavendiranathan P, Popović ZB, Hanna M, Plana JC, et al. Relative apical sparing of longitudinal strain using two-dimensional speckle-tracking echocardiography is both sensitive and specific for the diagnosis of cardiac amyloidosis. Heart. (2012) 98:1442–1448. doi: 10.1136/heartjnl-2012-302353

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Lagies R, Beck BB, Hoppe B, Sreeram N, Ten Cate FEAU. Apical sparing of longitudinal strain, left ventricular rotational abnormalities, and short-axis dysfunction in primary hyperoxaluria type 1. Circ Hear Fail. (2013) 6:45–47. doi: 10.1161/CIRCHEARTFAILURE.113.000432

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: machine learning, neural network, strain analysis, magnetic resonance, cardiomyopathy, hypertrophic

Citation: Satriano A, Afzal Y, Sarim Afzal M, Fatehi Hassanabad A, Wu C, Dykstra S, Flewitt J, Feuchter P, Sandonato R, Heydari B, Merchant N, Howarth AG, Lydell CP, Khan A, Fine NM, Greiner R and White JA (2020) Neural-Network-Based Diagnosis Using 3-Dimensional Myocardial Architecture and Deformation: Demonstration for the Differentiation of Hypertrophic Cardiomyopathy. Front. Cardiovasc. Med. 7:584727. doi: 10.3389/fcvm.2020.584727

Received: 17 July 2020; Accepted: 09 October 2020;
Published: 11 November 2020.

Edited by:

Matteo Cameli, University of Siena, Italy

Reviewed by:

Giulia Elena Mandoli, University of Siena, Italy
Alicia M. Maceira, Hospital Arnau de Vilanova, Spain

Copyright © 2020 Satriano, Afzal, Sarim Afzal, Fatehi Hassanabad, Wu, Dykstra, Flewitt, Feuchter, Sandonato, Heydari, Merchant, Howarth, Lydell, Khan, Fine, Greiner and White. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Alessandro Satriano, asatrian@ucalgary.ca

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.