Introduction {#Sec1}
============

Left atrial and ventricular volumes and ejection fraction are important diagnostic and prognostic parameters, widely used in daily practice. Indication for cardiac surgery in patients with valvular abnormalities such as bicuspid aortic valve disease, relies heavily on accurate left ventricular (LV) volume and function assessment. For many years, two-dimensional echocardiography (2DE) has been the most widely used modality for LV volumetric assessments; however, it is based on geometric assumptions which cause inaccuracy and the reproducibility remains suboptimal. Three-dimensional echocardiography (3DE) has largely overcome these drawbacks as it has the ability to visualize cardiac structures from any perspective, entailing an accurate quantitative and more reproducible evaluation of cardiac chambers. However, the use of 3DE in daily clinical practice has been hampered, because there is a learning curve for data acquisition, and 3D data analysis can be a time-consuming process, moreover there is need for experienced observers.

A fully automated and user-interference free algorithm may improve the feasibility of 3DE in daily clinical practice. Such a method has now been proposed in the new 'Heart Model' Software (LHM), which promises a rapid and accurate automated quantification of left atrial (LAV) and LV volumes and ejection fraction (LVEF). The feasibility and accuracy of the new "Heart Model" was recently compared to cardiac magnetic resonance imaging by Tsang et al. who concluded that this technique is strongly correlated with CMR, with a high reproducibility and short analysis time \[[@CR1], [@CR2]\]. However, it has not been reported whether this also applies to patients with valvular heart disease where high reproducibility and feasibility is very important.

Therefore, this study aims to assess feasibility and reproducibility of the 'Heart model' software in a prospective cohort study by comparing the results between echo and CT using conventional 2D, 2D-xPlane, 3D transthoracic echocardiography TTE (3DE) techniques in patients with a bicuspid aortic valve (BAV) with moderate to severe aortic valve stenosis or regurgitation.

Methods {#Sec2}
=======

Patients selection {#Sec3}
------------------

Patients with a BAV who visited the outpatient clinic for regular follow-up between October 2014 and March 2016 were consecutively included into this prospective study. All patients underwent the full study protocol on the same day. The study protocol consisted of physical examination, electrocardiography (ECG), 2D and 3D echocardiography and a cardiac CT scan. The study was approved by the medical ethical committee of the Erasmus medical center and informed consent was given by all patients who participated in the study. Two protocols were used reference methods were used to compare the LFM data: in the first protocol 3DE was used as a reference standard, in the second protocol the LHM was compared to CT.

Echocardiography {#Sec4}
----------------

### Image acquisition {#Sec5}

Two experienced sonographers (J.S.M., W.B.V.) performed a standard two-dimensional transthoracic echocardiogram (2DE). All studies were acquired in the left lateral decubitus position, in harmonic imaging using an EPIQ7 ultrasound system (Philips Medical Systems, Best, the Netherlands) equipped with a x5--1 matrix-array transducer (composed of 3040 elements with 1--5 MHz) A non-foreshortened apical four-chamber view (A4C) and two-chamber view (A2C) were recorded with manual or electronic rotation (iRotate) followed by a focused LV, A4C and A2C image. From the focused LV-A4C a true perpendicular image view (A2C) was acquired with xPlane mode, in order to retrieve both views from the same heart-beat. This was repeated with the focus on the true long axis of the LA \[[@CR3]\]. Real-time 3D-TTE was performed immediately after the 2D-TTE with the same ultrasound unit and transducer. A four-or six-beat full volume dataset of the LV and LA was acquired from the apical window during a single breath hold. Two extra datasets were acquired from the A4C view in the dedicated 'Heart Model' acquisition mode.

### Image analysis {#Sec6}

Analysis was performed by A.T.H. and J.S.M. All measurements were blinded to patient specific information. Before measurements were performed, the quality of each dataset was evaluated by both observers. Patients were excluded from further analysis in cases of poor image quality (e.g. poor endocardial visualization).

#### 2D echocardiography {#Sec7}

LV end diastolic volume (LVEDV), LV end systolic volume (LVESV) and LVEF were calculate using the Simpsons bi-plane method of disk summation, as stated in the guidelines, from the standard A4C and A2C and apical xPlane images \[[@CR4]\]. LAV was calculated using the area length method.

#### 3D echocardiography {#Sec8}

Manual: LV volumes and LVEF were measured using commercially available software (QLab-3DQ, Philips medical systems). The user aligned the multiplanar view to obtain the true long axis of the LV in the A4C and A2C view. Landmarks were placed on the mitral annulus and apex. The endocardial border was traced automatically and adjusted manually where needed. For the LA volume, the true long-axis was aligned using the multiplanar mode in the end-systolic frame and the contour traced as mentioned above. The 3D dataset was scored feasible when the entire cardiac contour could be traced.

Automatic: Offline fully automatic analysis of the datasets was performed using the Q Lab advanced 'Heart Model' analysis software. This software has previously described in detail before by Tsang et al. \[[@CR1], [@CR2]\]. In brief, this software detects the endocardial surfaces by using an adaptive algorithm. This identifies a global end-diastolic shape which it uses in combination with motion detection to determine an end-systolic cavity \[[@CR1], [@CR2]\]. The program combines information from a database of 1000 3D TTE studies and its endocardial surface detection to model the LA and LV. Afterwards it matches features from the known datasets to the current patient for which it needs, much like manual measurements, a minimum of approximately 14 or 15 LV segments. The final model (Fig. [1](#Fig1){ref-type="fig"}) is displayed with the possibility to manually edit the contours if the user deems this necessary. To better estimate the added value over existing 3D techniques we chose not to manually edit the contours.

Fig. 1Example of a model created by the LHM. The model as generated by the LHM of left atrial (LA and left ventricular (LV) volumes

Cardiac computed tomography {#Sec9}
---------------------------

### Data acquisition {#Sec10}

Acquisition was performed on a 3rd generation dual-source CT (Somatom Force, Siemens Healthcare, Forchheim, Germany). Retrospective ECG gated spiral acquisition was used, with a mean dose length product (DLP) of (estimated effective dose 5 mSv, using k = 0.017) 362 mGy-cm kVp was modulated to patient size. Besides patient size and the selected kV the mA was modulated to the heartrate to provide a high mA pulse during 1--40% of the RR interval. The pitch was adapted to increase proportionally with higher heartrates. Reconstructions were made with a medium smooth kernel, with a slice thickness of 1.5 mm at an increment of 0.4 mm at every 5% of the RR interval.

### Image processing {#Sec11}

Images were analyzed semi-automatically using Syngo Via software (vb 10., Siemens, Forcheim Germany). Image analysis was performed by AC, with 1 year of cardiac CT experience. All cardiac phases were analyzed the software automatically detects the ED and ES phase, this was manually changed if needed. The endo and epicardial contours are automatically placed by the software and manually adjusted were need. The papillary muscles and if present trabeculations were included into the LV lumen. The basal plane was selected perpendicular to the short axis at the level of the mitral valve. Care was used to make sure the basal plane was on the same level in the ED and ES phase.

Statistics {#Sec12}
----------

The IBM SPSS^®^ statistics 21.0 software was used to analyze the data. Continuous variables were presented as mean ± standard deviation (SD) or as median with an interquartile range. Categorical variables were presented as frequencies and percentages. We tested for normality by calculating Z-values of skewness and kurtosis, using the Shapiro--Wilk test and by visually assessing the data. For comparison of normally distributed continuous variables between two groups the student's *t* test was used. To quantify correlations the Pearson or Spearman correlation test was applied. Intra-observer and inter-observer agreement between two investigators (A.T.H, J.S.M.) were assessed by repeated analysis of the same images in a third of the dataset a month after initial analysis at the same images and blinded to the initial results. The limits of agreement between two measurements were determined as the mean of the differences (bias) ± 1.96 SD and presented in a Bland--Altman plot \[[@CR5]\]. Additionally, the coefficient of variation (COV) was provided to compare the dispersion of two variables. The COV was defined as either the SD of the differences of two measurements divided by the mean of their means. The statistical tests were two sided and a p-value below 0.05 was considered significant.

Results {#Sec13}
=======

Fifty-three patients with a BAV were included; baseline characteristics are presented in Table [1](#Tab1){ref-type="table"}. Eight patients also had Turner syndrome. Four patients previously underwent aortic coarctation resection and three patients underwent a balloon dilatation of their stenotic aortic valve. One patient underwent closure of a type II atrial septal defect.

Table 1Baseline characteristics of the study population (n = 37)ParameterMedian (IQR)Baseline Men, n (%)25 (68) Age (years)35.2 (23) Height (cm)178 (26) Weight (kg)72 (24) BMI (kg/m^2^)23.9 (3.0) SBP (mmHg)122 (20) DPB (mmHg)80 (19) 3DE LVEF \< 50%, n (%)16 (43)Mitral valve E-wave (m/s)0.70 (0.2) A-wave (m/s)0.50 (0.2) E/A-ratio1.2 (0.9) DT (ms)209 (67) E' septal, cm/s7.8 (3.1) Ee'-ratio9.0 (3.6)Aortic valve BAV (n = 32)  No AoI, n= (%)5 (16)  Mild AoI, n= (%)21 (66)  Moderate AoI, n= (%)5 (16)  Severe AoI, n= (%)1 (3)  Peak velocity (m/s)2.65 (1.6)  VT (cm)52.8 (45)  Gradient (m/s)28 (33) TS\* (n = 5)  No AoI, n= (%)4 (80)  Mild AoI, n= (%)1 (20)  Moderate, AoI n= (%)0 (0)  Severe AoI, n= (%)0 (0)  Peak velocity (m/s)1.4 (0.6)  VTI (cm)28.4 (12)  Gradient (m/s)8 (6.5)Data are expressed as median and inter quartile range (IQR) or as 'n=, (%)' for the variables 'gender', 3DE LVEF\<50% and 'aortic insufficiency'*AoI* aortic insufficiency, *BMI* body mass index, *SBP* systolic blood pressure, *DBP* diastolic blood pressure, *DT* deceleration time\*Subgroup of patients with Turner syndrome (TS) and a bicuspid aortic valve (BAV)

In all 53 patients, 2DE measurements could be performed and functional echo parameters were measured, as presented in Table [1](#Tab1){ref-type="table"}. In four patients, no 3D dataset could be acquired. Additionally, 12 patients were excluded based on poor imaging quality.

Left ventricle {#Sec14}
--------------

### Protocol 1: LHM versus conventional echocardiography {#Sec15}

There was a good correlation between the LHM and the manual bi-and xPlane measurements for LVEDV, LVESV, and LVEF. In Table [2](#Tab2){ref-type="table"} the measurements of LV volumes and function and LA maximal volume are presented for all different methods. Figure [2](#Fig2){ref-type="fig"}a, b show the Bland--Altman plots for EF by 2DE bi-plane and xPlane compared to the LHM. The results of the agreement analysis between the LHM and the measurements based on the 2DE are shown in Table [2](#Tab2){ref-type="table"}. For clarity purposes, we did not include mutual correlations between the different echo modalities; however, Bi-plane and xPlane had a high correlation for LVEDV (r = 0.977), LVESV (r = 0.978) and LVEF (r = 0.702). Additionally, both methods correlated strongly with conventional 3D as expected.

Table 2Correlations between LHM and four different methods of volumetric chamber quantificationMethodPhaseMean, SDPearson's r ^†^BiasLower LOAUpper LOALHM (n = 37)EDV (ml)146 ± 48ESV (ml)77 ± 29EF (%)47 ± 5LA (ml)61 ± 19LHM 2nd dataset (n = 12)EDV (ml)143 ± 610.99\*\*−0.8−97ESV (ml)78 ± 340.99\*\*1−79EF (%)45 ± 40.58\*−1.5−74LA (ml)57 ± 190.98\*\*−1−972DE Bi-plane (n = 37)EDV (ml)145 ± 540.93\*\*0.06−4141ESV (ml)77 ± 310.88\*\*0.2−29.129.4EF (%)47 ± 80.63\*\*0.5−11.812.8LA (ml)43 ± 160.61\*\*17−13.647.42DE xPlane (n = 37)EDV (ml)143 ± 520.94\*\*2−3236.4ESV (ml)77 ± 300.88\*\*0.8−27.128.8EF (%)47 ± 80.43\*\*0.6−1314.2LA (ml)42 ± 160.69\*\*18−8.845.23D (n = 37)EDV (ml)143 ± 500.97\*\*2−2327.6ESV (ml)71 ± 260.91\*\*6−16.929.8EF (%)50 ± 70.51\*\*−3−15.59.8LA (ml)53 ± 190.81\*\*9−13.132.1CT (n = 37)EDV (ml)185 ± 630.88\*\*−42−102.318.9ESV (ml)67 ± 240.81\*\*10−24.143.6EF (%)64 ± 50.2416−28.3−4.7Data are presented as mean and SD*EDV* end-diastolic volume,*ESV* end-systolic volume,*EF* ejection fraction,*LA*left atrium,*LOA* limit of agreement,*COV*coefficient of variation^†^Compared with the LHM, \*\*p\<0.01, \*p\<0.05. A negative mean implies a smaller value was given by the LHM

Fig. 2Bland--Altman plots demonstrating inter-modality agreement of EF and LAV in 2DE biplane (*panel* **a** and **e**), 2DE xPlane (*panel* **b** and **f**), 3DE (*panel* **c** and **g**) compared to the LHM. Ejection fraction measured by bi-plane 2DE (**a**), by xPlane 2DE (**b**) and by 3DE (**c**). Left atrial volume measured by bi-plane 2DE (**e**), by xPlane 2DE (**f**) and by 3DE (**g**). The*solid lines* depict the mean difference of the two measurements; the*dashed lines* depict the limits of agreement.*COV* coefficient of variation

The LHM correlated strongly with the 3D LV measurements as shown in Table [2](#Tab2){ref-type="table"}. The LHM seems to estimate slightly larger LVEDV and LVESV and a smaller LVEF compared to manual 3D (Table [2](#Tab2){ref-type="table"}). Figure [2](#Fig2){ref-type="fig"}c shows the Bland--Altman plot for EF measured by 3DE. 16 patients (30%) had to be excluded due to poor imaging quality, which is considered acceptable in a routine clinical setting for 3DE.

### Protocol 2: LHM versus CT {#Sec16}

In the second protocol, where CT was used as a golden standard method, agreement was found for EDV and ESV (r = 0.88 and r = 0.81); mean differences were 42 ml (23%) and 10 ml (15%) respectively (p \< 0.001 and p = 0.002). However, compared to CT the LHM seems to estimate smaller LVEDV (Table [2](#Tab2){ref-type="table"}) and consequently also produces a smaller LVEF (mean difference 16%, p \< 0.001) as shown by the positive difference in EF between CT and the LHM in Fig. [3](#Fig3){ref-type="fig"}. The correlations for Bi-plane, xPlane and conventional 3D with CT for LVEDV, LVESV and LVEF (Fig. [4](#Fig4){ref-type="fig"}) were comparable; a relatively small LVEDV and therefore a lower LVEF compared to CT. Three patients did not undergo a CT scan for personal reasons, those three patients are excluded from analysis.

Fig. 3Bland--Altman plots demonstrating inter-modality agreement between CT and 3DE and LHM. Agreement in end-diastolic (EDV) and end-systolic volume (ESV) and ejection fraction (EF) comparing CT versus conventional 3D echocardiography (*top row*) and CT versus the LHM (*bottom row*)

Fig. 4Bland--Altman plots demonstrating inter-modality agreement between CT and 3DE and LHM. Agreement in end-diastolic (EDV) and end-systolic volume (ESV) and ejection fraction (EF) comparing CT versus 2D bi-plane echocardiography (*top row*) and CT versus xPlane echocardiography (*bottom row*)

### Inter- and intra-observer variability {#Sec17}

The LHM has no intra-, or inter-observer variation as the model produced the exact same measurements when the same dataset is used (mean difference: 0 ± 0, p \< 0.001). In addition, when the measurements were repeated on a second acquisition of the same patients (n = 12) there was some variation but generally there was good agreement reflected by a small bias and narrow limits of agreement (Table [2](#Tab2){ref-type="table"}). The inter-observer variability for the LVEDV and LVESV on conventional 3DE was calculated (r = 0.909, p \< 0.001 and r = 0.862, p \< 0.001) in 14 (38%) patients.

Left atrial volume {#Sec18}
------------------

The LAV was estimated consistently significantly larger by the LHM (Table [2](#Tab2){ref-type="table"}) compared with on both 2DE and 3DE as is shown in Fig. [2](#Fig2){ref-type="fig"}e--g. This 'overestimation' was the most evident comparing to the 2D methods and less pronounced compared with manual 3D measurements. The correlation between 2D Bi-plane and 2D xPlane measurements of the LA volume was very high (r = 0.946, p \< 0.001). When comparing the LAV measured on 2D with the LA volume measured on 3D, xPlane outperforms Bi-plane (r = 0.632, p \< 0.001 and r = 0.551, p = 0.002 respectively).

Discussion {#Sec19}
==========

The main findings of this study can be summarized as follows:

Automated chamber quantification is feasible in patients with bicuspid aortic valve disease in a routine clinical setting.The 'Heart Model' provides accurate automatic measurements of LVEDV, LVESV and LVEF and LAV compared to 2D and 3D echocardiography.There is no inter or intra-observer variability and very little 'inter-dataset' variability.

Left ventricular assessment {#Sec20}
---------------------------

In daily clinical practice LV function is routinely estimated by bi-plane Simpson method of disk-summation; however, this requires sufficient experience and has limitations in accuracy and reproducibility, especially due to the geometric assumptions of the shape of the ventricular or atrial cavity inherent to 2DE \[[@CR4]\]. Moreover, the lack of a third dimension is generally considered to result in high inter-measurement variability and limits endocardial visualization, predominantly of apical lateral segments. Foreshortening of the LV often performed in an attempt to alleviate this shortcoming causes reduced accuracy and reproducibility \[[@CR6]\]. Volumetric quantification from 3D data sets allows frame-by-frame detection of endocardial surface and does not require manual image plane positioning or geometric assumption \[[@CR7]\] and has furthermore been shown to have a higher reproducibility than 2DE \[[@CR8]--[@CR10]\].

The bias of CT and 3D echocardiography is approximately as low as the bias of MRI when estimating LVEF \[[@CR11]\]. Additionally, previous studies showed that for LVEF CT has the best correlation with MRI. In this light it is remarkable that in our study the correlation with CT is weak, as the EDV seems to be systematically estimated to be larger by CT than by the LHM or other echo modalities \[[@CR11], [@CR12]\]. This systematic underestimation may in part reflect an inter observer variability, as more of the trabeculations were included in the LV cavity for CT, entailing a larger EDV. Also, the trabeculations in this population may be more pronounced than in healthy subjects, adding to the observed difference.

Left heart chamber quantification by 3DE is hampered by several factors, mainly the ease with which these techniques can be used in daily clinical practice as a degree of experience is required. However, especially in the growing and aging population of patients with valvular and congenital heart disease, where patients regularly undergo extensive echocardiographic evaluation, user-independent and non-invasive follow-up imaging is much needed.

Added value of LHM {#Sec21}
------------------

Accurate and reproducible measurements of left chamber volumes are very important in clinical practice, as they correlate with prognosis and determine treatment strategies. Moreover, in order to test new therapies, changes in these parameters must be accurate to demonstrate the efficacy of medical therapy or intervention. Our results indicate that intra-observer, inter-observer, and test--retest reproducibility of this method are very high and even exceptionally for intra-observer variability (0%) since this algorithm has no human interaction (i.e., phase selection and contouring are automatic).

We had to exclude 16 patients (30%) in total due to imaging quality limitations which is comparable to conventional 3DE in a routine clinical setting.

The main concern with 3D echocardiography is the accessibility in terms of time and skill required to produce accurate and reproducible results. The algorithm described in this study promises to improve on these points. This study demonstrates a high feasibility and accuracy in a population of patients with BAV disease, a population where fast, reliable, user-independent and non-invasive follow-up imaging by echocardiography is imperative.

Left atrial volume {#Sec22}
------------------

Another remarkable finding is the relatively large LA volume estimated by the LHM compared to 2DE measurements. First, we suggested that this could be explained by the inclusion of the pulmonary vein orifice into the left atrial cavity by the algorithm. However, when carefully re-evaluating the LA contours, most seemed to adequately follow the left atrial walls. Another explanation could be that it is actually an underestimation of the 2DE measurements. This can partly be explained by the use of the length area formula which assumes the LA to be ellipsoid, which is evidently not always the case \[[@CR4]\]. The correlation of LA volume measurements with 3D TTE was better. Still, inherent flaws of 2D imaging in combination with LA anatomy could contribute to this discrepancy. In the standard 2DE the LA measurements are performed on the A4C and A2C views which are focused on the true long axis of LV. The true long axis LA may not be in the same plane as the LV and therefore may appear foreshortened on the apical four chamber view. In xPlane a special focus view on the true LA long axis was acquired which better correlated with the LHM.

Limitations {#Sec23}
-----------

The golden standard for quantitative volumetric heart chamber assessment is cardiac magnetic resonance imaging (MRI). Unfortunately, no cine cardiac images were available in this study. Therefore, no comparison with MRI could be made; consequently, we used volumetric data from CT as an additional source of validation. Additionally, we feel that EDV by 2DE has been slightly underestimated and therefore 3DE is best used as a reference modality. We were strict when it came to image quality therefore we had to exclude 12 patients from analysis. However, a sufficiently large cohort could be analyzed.

Conclusion {#Sec24}
==========

Automated chamber quantification is feasible and accurate in patients with bicuspid aortic valve disease in a routine clinical setting. The 'Heart Model' provides accurate automatic measurements of LVEDV, LVESV and LVEF and LAV and has a high reproducibility between dataset and no inter or intra-observer variability. It does however seem to underestimate end-diastolic volume compared to CT.

Funding {#FPar1}
=======

This study was funded by a grant of the Dutch Heart Foundation (The Hague, The Netherlands, Grant No. 2013T093). The Dutch Heart Foundation had no influence on study design; in the collection, analysis and interpretation of data; in the writing of the report; or the decision to submit the article for publication.

Conflict of interest {#FPar2}
====================

All authors declare that they have no conflict of interest.
