Four-dimensional flow cardiovascular magnetic resonance in tetralogy of Fallot: a systematic review

BACKGROUND

Patients with repaired Tetralogy of Fallot (rTOF) often develop cardiovascular dysfunction and require regular imaging to evaluate deterioration and time interventions such as pulmonary valve replacement. Four-dimensional flow cardiovascular magnetic resonance (4D flow CMR) enables detailed assessment of flow characteristics in all chambers and great vessels. We performed a systematic review of intra-cardiac 4D flow applications in rTOF patients, to examine clinical utility and highlight optimal methods for evaluating rTOF patients.

METHODS

A comprehensive literature search was undertaken in March 2020 on Google Scholar and Scopus. A modified version of the Critical Appraisal Skills Programme (CASP) tool was used to assess and score the applicability of each study. Important clinical outcomes were assessed including similarities and differences.

RESULTS

Of the 635 articles identified, 26 studies met eligibility for systematic review. None of these were below 59% applicability on the modified CASP score. Studies could be broadly classified into four groups: (i) pilot studies, (ii) development of new acquisition methods, (iii) validation and (vi) identification of novel flow features. Quantitative comparison with other modalities included 2D phase contrast CMR (13 studies) and echocardiography (4 studies). The 4D flow study applications included stroke volume (18/26;69%), regurgitant fraction (16/26;62%), relative branch pulmonary artery flow(4/26;15%), systolic peak velocity (9/26;35%), systemic/pulmonary total flow ratio (6/26;23%), end diastolic and end systolic volume (5/26;19%), kinetic energy (5/26;19%) and vorticity (2/26;8%).

CONCLUSIONS

4D flow CMR shows potential in rTOF assessment, particularly in retrospective valve tracking for flow evaluation, velocity profiling, intra-cardiac kinetic energy quantification, and vortex visualization. Protocols should be targeted to pathology. Prospective, randomized, multi-centered studies are required to validate these new characteristics and establish their clinical use.

Statistical shape modeling of the left ventricle: myocardial infarct classification challenge

Statistical shape modeling is a powerful tool for visualizing and quantifying geometric and functional patterns of the heart. After myocardial infarction (MI), the left ventricle typically remodels in response to physiological challenges. Several methods have been proposed in the literature to describe statistical shape changes. Which method best characterizes left ventricular remodeling after MI is an open research question. A better descriptor of remodeling is expected to provide a more accurate evaluation of disease status in MI patients. We therefore designed a challenge to test shape characterization in MI given a set of three-dimensional left ventricular surface points. The training set comprised 100 MI patients, and 100 asymptomatic volunteers (AV). The challenge was initiated in 2015 at the Statistical Atlases and Computational Models of the Heart workshop, in conjunction with the MICCAI conference. The training set with labels was provided to participants, who were asked to submit the likelihood of MI from a different (validation) set of 200 cases (100 AV and 100 MI). Sensitivity, specificity, accuracy and area under the receiver operating characteristic curve were used as the outcome measures. The goals of this challenge were to (1) establish a common dataset for evaluating statistical shape modeling algorithms in MI, and (2) test whether statistical shape modeling provides additional information characterizing MI patients over standard clinical measures. Eleven groups with a wide variety of classification and feature extraction approaches participated in this challenge. All methods achieved excellent classification results with accuracy ranges from 0.83 to 0.98. The areas under the receiver operating characteristic curves were all above 0.90. Four methods showed significantly higher performance than standard clinical measures. The dataset and software for evaluation are available from the Cardiac Atlas Project website1.

An Open Benchmark Challenge for Motion Correction of Myocardial Perfusion MRI

Cardiac magnetic resonance perfusion examinations enable noninvasive quantification of myocardial blood flow. However, motion between frames due to breathing must be corrected for quantitative analysis. Although several methods have been proposed, there is a lack of widely available benchmarks to compare different algorithms. We sought to compare many algorithms from several groups in an open benchmark challenge. Nine clinical studies from two different centers comprising normal and diseased myocardium at both rest and stress were made available for this study. The primary validation measure was regional myocardial blood flow based on the transfer coefficient (Ktrans), which was computed using a compartment model and the myocardial perfusion reserve (MPR) index. The ground truth was calculated using contours drawn manually on all frames by a single observer, and visually inspected by a second observer. Six groups participated and 19 different motion correction algorithms were compared. Each method used one of three different motion models: rigid, global affine, or local deformation. The similarity metric also varied with methods employing either sum-of-squared differences, mutual information, or cross correlation. There were no significant differences in Ktrans or MPR compared across different motion models or similarity metrics. Compared with the ground truth, only Ktrans for the sum-of-squared differences metric, and for local deformation motion models, had significant bias. In conclusion, the open benchmark enabled evaluation of clinical perfusion indices over a wide range of methods. In particular, there was no benefit of nonrigid registration techniques over the other methods evaluated in this study. The benchmark data and results are available from the Cardiac Atlas Project ( www.cardiacatlas.org).

Left ventricular function and regional strain with subtly-tagged steady-state free precession feature tracking

PURPOSE

To provide regional strain and ventricular volume from a single acquisition, using subtly tagged steady-state free precession (SubTag SSFP) feature tracking.

MATERIALS AND METHODS

The effects on regional strain of tag strength in gradient recalled echo (GRE) tagging, flip angle in untagged balanced SSFP, and both in SubTag SSFP were examined in the mid left ventricle of 15 healthy volunteers at 3T. Optimal parameters were determined from varying both tag strength and SSFP flip angle using full tag saturation GRE as the reference standard. SubTag SSFP was acquired in 15 additional healthy volunteers for whole-heart volume and strain assessment using the optimized parameters. Values measured by two image analysts were compared to clinical reference standards from untagged SSFP (volumes) and GRE tagging (strains).

RESULTS

Regional strain accuracy was maintained with decreasing total tagging flip angle (β); less than 3% differences for β ≥ 26°. For untagged SSFP flip angle (α), whole-wall strain differences became statistically significant when α < 40°. A SubTag SSFP acquisition with α = 40° and β = 46° showed the best combination of tagging strength, blood-myocardial contrast, and tag persistence at end-systole for regional strain estimation. SubTag SSFP also showed excellent agreement with untagged SSFP for volumetrics (percent difference: end-diastolic volume = 0.6%, end-systolic volume = 0.4%, stroke volume = 1.2%, ejection fraction = 0.6%, mass = 1.1%).

CONCLUSION

Feature tracking for regional myocardial strain assessment is dependent on image features, mainly the tag strength, persistence, and image contrast. SubTag SSFP balances these criteria to provide accurate regional strain and volumetric assessment in a single acquisition.

LEVEL OF EVIDENCE

1 Technical Efficacy Stage: 2 J. Magn. Reson. Imaging 2018;47:787-797.

Orthogonal decomposition of left ventricular remodeling in myocardial infarction

Left ventricular size and shape are important for quantifying cardiac remodeling in response to cardiovascular disease. Geometric remodeling indices have been shown to have prognostic value in predicting adverse events in the clinical literature, but these often describe interrelated shape changes. We developed a novel method for deriving orthogonal remodeling components directly from any (moderately independent) set of clinical remodeling indices. Results: Six clinical remodeling indices (end-diastolic volume index, sphericity, relative wall thickness, ejection fraction, apical conicity, and longitudinal shortening) were evaluated using cardiac magnetic resonance images of 300 patients with myocardial infarction, and 1991 asymptomatic subjects, obtained from the Cardiac Atlas Project. Partial least squares (PLS) regression of left ventricular shape models resulted in remodeling components that were optimally associated with each remodeling index. A Gram–Schmidt orthogonalization process, by which remodeling components were successively removed from the shape space in the order of shape variance explained, resulted in a set of orthonormal remodeling components. Remodeling scores could then be calculated that quantify the amount of each remodeling component present in each case. A one-factor PLS regression led to more decoupling between scores from the different remodeling components across the entire cohort, and zero correlation between clinical indices and subsequent scores. Conclusions: The PLS orthogonal remodeling components had similar power to describe differences between myocardial infarction patients and asymptomatic subjects as principal component analysis, but were better associated with well-understood clinical indices of cardiac remodeling. The data and analyses are available from www.cardiacatlas.org.

Three-Dimensional Volumetric Assessment of Diastolic Function by Cardiac Magnetic Resonance Imaging: The Multi-Ethnic Study of Atherosclerosis (MESA)

Background:

Cardiac Magnetic Resonance is in need of a simple and robust method for diastolic function assessment that can be done with routine protocol sequences.

Objective:

To develop and validate a three-dimensional (3D) model-based volumetric assessment of diastolic function using cardiac magnetic resonance (CMR) imaging and compare the results obtained with the model with those obtained by echocardiography.

Methods:

The study participants provided written informed consent and were included if having undergone both echocardiography and cine steady-state free precession (SSFP) CMR on the same day. Guide points at the septal and lateral mitral annulus were used to define the early longitudinal relaxation rate (E'), while a time-volume curve from the 3D model was used to assess diastolic filling parameters. We determined the correlation between 3D CMR and echocardiography and the accuracy of CMR in classifying the diastolic function grade.

Results:

The study included 102 subjects. The E/A ratio by CMR was positively associated with the E/A ratio by echocardiography (r = 0.71, p < 0.0001). The early diastolic relaxation velocity by tissue Doppler and longitudinal relaxation rate for the lateral mitral annulus displacement were positively associated (p = 0.007), as were the ratio between Doppler E/e' and CMR E/E' (p = 0.01). CMR-determined normalized peak E (NE) and deceleration time (DT) were able to predict diastolic dysfunction (areas under the curve [AUCs] = 0.70 and 0.72, respectively). In addition, the lateral E/E' ratio showed good utility in identifying diastolic dysfunction (AUC = 0.80). Overall, echocardiography and CMR interobserver and intraobserver agreements were excellent (intraclass correlation coefficient range 0.72 - 0.97).

Conclusion:

3D modeling of standard cine CMR images was able to identify study subjects with reduced diastolic function and showed good reproducibility, suggesting a potential for a routine diastolic function assessment by CMR.

Fundamento:

A ressonância magnética cardíaca necessita de um método simples e robusto para a avaliação da função diastólica que pode ser feito com sequências protocolares de rotina.

Objetivo:

Desenvolver e validar a avaliação volumétrica da função diastólica através de um modelo tridimensional (3D) com utilização de imagens de ressonância magnética cardíaca (RMC) e comparar os resultados obtidos com este modelo com os obtidos por ecocardiografia.

Métodos:

Os participantes do estudo assinaram um termo de consentimento e foram incluídos se tivessem sido submetidos no mesmo dia tanto à ecocardiografia quanto à cine RMC com precessão livre no estado estacionário (steady-state free precession, SSFP). Pontos-guia foram utilizados no anel mitral septal e lateral para definir a velocidade de estiramento no início da diástole (E'), enquanto curvas de volume-tempo do modelo 3D foram utilizadas para avaliar os parâmetros de enchimento diastólico. Foram determinadas a correlação entre a RMC 3D e a ecocardiografia, além da acurácia da RMC em classificar o grau de função diastólica.

Resultados:

Ao todo, 102 sujeitos foram incluídos no estudo. A razão E/A pela RMC esteve positivamente associada com a razão E/A obtida pela ecocardiografia (r = 0,71, p < 0,0001). Estiveram positivamente associadas a velocidade de relaxamento diastólico inicial ao Doppler tecidual e a velocidade de relaxamento longitudinal de deslocamento do anel mitral lateral (p = 0,007), bem como a razão entre E/e' por Doppler e E/E' pela RMC (p = 0,01). A velocidade normalizada de pico de enchimento (EM) determinada pela RMC e o tempo de desaceleração (TD) foram capazes de predizer a disfunção diastólica (áreas sob a curva [AUCs] = 0,70 e 0,72, respectivamente). Além disso, a razão E/E' lateral mostrou boa utilidade para a identificação da disfunção diastólica (AUC = 0,80). No geral, a ecocardiografia e a RMC apresentaram excelente concordância interobservador e intraobservador (coeficiente de correlação intraclasse 0,72 - 0,97).

Conclusão:

Uma modelagem 3D de imagens padrões de cine RMC foi capaz de identificar os indivíduos do estudo com função diastólica reduzida e mostrou uma boa reprodutibilidade, sugerindo ter potencial na avaliação rotineira da função diastólica por RMC.

Long-term cardiovascular outcome following fetal anaemia and intrauterine transfusion: a cohort study

OBJECTIVE

To compare long-term cardiovascular outcomes in survivors of fetal anaemia and intrauterine transfusion with those of non-anaemic siblings.

DESIGN

Retrospective cohort study.

SETTING

Auckland, New Zealand.

PARTICIPANTS

Adults who received intrauterine transfusion for anaemia due to rhesus disease (exposed) and their unexposed sibling(s).

EXPOSURE

Fetal anaemia requiring intrauterine transfusion.

MAIN OUTCOME MEASURES

Anthropometry, blood pressure, lipids, heart rate variability and cardiac MRI, including myocardial perfusion.

RESULTS

Exposed participants (n=95) were younger than unexposed (n=92, mean±SD 33.7±9.3 vs 40.1±10.9 years) and born at earlier gestation (34.3±1.7 vs 39.5±2.1 weeks). Exposed participants had smaller left ventricular volumes (end-diastolic volume/body surface area, difference between adjusted means -6.1, 95% CI -9.7 to -2.4 mL/m²), increased relative left ventricular wall thickness (difference between adjusted means 0.007, 95% CI 0.001 to 0.012 mm.m²/mL) and decreased myocardial perfusion at rest (ratio of geometric means 0.86, 95% CI 0.80 to 0.94). Exposed participants also had increased low frequency-to-high frequency ratio on assessment of heart rate variability (ratio of geometric means 1.53, 95% CI 1.04 to 2.25) and reduced high-density lipoprotein concentration (difference between adjusted means -0.12, 95% CI -0.24 to 0.00 mmol/L).

CONCLUSIONS

This study provides the first evidence in humans that cardiovascular development is altered following exposure to fetal anaemia and intrauterine transfusion, with persistence of these changes into adulthood potentially indicating increased risk of cardiovascular disease. These findings are relevant to the long-term health of intrauterine transfusion recipients, and may potentially also have implications for adults born preterm who were exposed to anaemia at a similar postconceptual age.

Image-Based Investigation of Human in Vivo Myofibre Strain

Cardiac myofibre deformation is an important determinant of the mechanical function of the heart. Quantification of myofibre strain relies on 3D measurements of ventricular wall motion interpreted with respect to the tissue microstructure. In this study, we estimated in vivo myofibre strain using 3D structural and functional atlases of the human heart. A finite element modelling framework was developed to incorporate myofibre orientations of the left ventricle (LV) extracted from 7 explanted normal human hearts imaged ex vivo with diffusion tensor magnetic resonance imaging (DTMRI) and kinematic measurements from 7 normal volunteers imaged in vivo with tagged MRI. Myofibre strain was extracted from the DTMRI and 3D strain from the tagged MRI. We investigated: i) the spatio-temporal variation of myofibre strain throughout the cardiac cycle; ii) the sensitivity of myofibre strain estimates to the variation in myofibre angle between individuals; and iii) the sensitivity of myofibre strain estimates to variations in wall motion between individuals. Our analysis results indicate that end systolic (ES) myofibre strain is approximately homogeneous throughout the entire LV, irrespective of the inter-individual variation in myofibre orientation. Additionally, inter-subject variability in myofibre orientations has greater effect on the variabilities in myofibre strain estimates than the ventricular wall motions. This study provided the first quantitative evidence of homogeneity of ES myofibre strain using minimally-invasive medical images of the human heart and demonstrated that image-based modelling framework can provide detailed insight to the mechanical behaviour of the myofibres, which may be used as a biomarker for cardiac diseases that affect cardiac mechanics.

Estimation of myocardial strain from non-rigid registration and highly accelerated cine CMR

Sparsely sampled cardiac cine accelerated acquisitions show promise for faster evaluation of left-ventricular function. Myocardial strain estimation using image feature tracking methods is also becoming widespread. However, it is not known whether highly accelerated acquisitions also provide reliable feature tracking strain estimates. Twenty patients and twenty healthy volunteers were imaged with conventional 14-beat/slice cine acquisition (STD), 4× accelerated 4-beat/slice acquisition with iterative reconstruction (R4), and a 9.2× accelerated 2-beat/slice real-time acquisition with sparse sampling and iterative reconstruction (R9.2). Radial and circumferential strains were calculated using non-rigid registration in the mid-ventricle short-axis slice and inter-observer errors were evaluated. Consistency was assessed using intra-class correlation coefficients (ICC) and bias with Bland-Altman analysis. Peak circumferential strain magnitude was highly consistent between STD and R4 and R9.2 (ICC = 0.876 and 0.884, respectively). Average bias was -1.7 ± 2.0 %, p < 0.001, for R4 and -2.7 ± 1.9 %, p < 0.001 for R9.2. Peak radial strain was also highly consistent (ICC = 0.829 and 0.785, respectively), with average bias -11.2 ± 18.4 %, p < 0.001, for R4 and -15.0 ± 21.2 %, p < 0.001 for R9.2. STD circumferential strain could be predicted by linear regression from R9.2 with an R² of 0.82 and a root mean squared error of 1.8 %. Similarly, radial strain could be predicted with an R² of 0.67 and a root mean squared error of 21.3 %. Inter-observer errors were not significantly different between methods, except for peak circumferential strain R9.2 (1.1 ± 1.9 %) versus STD (0.3 ± 1.0 %), p = 0.011. Although small systematic differences were observed in strain, these were highly consistent with standard acquisitions, suggesting that accelerated myocardial strain is feasible and reliable in patients who require short acquisition durations.

Dynamically scaled phantom phase contrast MRI compared to true-scale computational modeling of coronary artery flow

PURPOSE

To examine the feasibility of combining computational fluid dynamics (CFD) and dynamically scaled phantom phase-contrast magnetic resonance imaging (PC-MRI) for coronary flow assessment.

MATERIALS AND METHODS

Left main coronary bifurcations segmented from computed tomography with bifurcation angles of 33°, 68°, and 117° were scaled-up ∼7× and 3D printed. Steady coronary flow was reproduced in these phantoms using the principle of dynamic similarity to preserve the true-scale Reynolds number, using blood analog fluid and a pump circuit in a 3T MRI scanner. After PC-MRI acquisition, the data were segmented and coregistered to CFD simulations of identical, but true-scale geometries. Velocities at the inlet region were extracted from the PC-MRI to define the CFD inlet boundary condition.

RESULTS

The PC-MRI and CFD flow data agreed well, and comparison showed: 1) small velocity magnitude discrepancies (2-8%); 2) with a Spearman's rank correlation ≥0.72; and 3) a velocity vector correlation (including direction) of r(2) ≥ 0.82. The highest agreement was achieved for high velocity regions with discrepancies being located in slow or recirculating zones with low MRI signal-to-noise ratio (SNRv ) in tortuous segments and large bifurcating vessels.

CONCLUSION

Characterization of coronary flow using a dynamically scaled PC-MRI phantom flow is feasible and provides higher resolution than current in vivo or true-scale in vitro methods, and may be used to provide boundary conditions for true-scale CFD simulations. J. MAGN. RESON. IMAGING 2016;44:983-992.

Creating shape templates for patient specific biventricular modeling in congenital heart disease

Survival rates for infants with congenital heart disease (CHD) are improving, resulting in a growing population of adults with CHD. However, the analysis of left and right ventricular function is very time-consuming owing to the variety of congenital morphologies. Efficient customization of patient geometry and function depends on high quality shape templates specifically designed for the application. In this paper, we combine a method for creating finite element shape templates with an interactive template customization to patient MRI examinations. This enables different templates to be chosen depending on patient morphology. To demonstrate this pipeline, a new biventricular template with 162 elements was created and tested in place of an existing 82-element template. The method was able to provide fast interactive biventricular analysis with 0.31 sec per edit response time. The new template was customized to 13 CHD patients with similar biventricular topology, showing improved performance over the previous template and good agreement with clinical indices.

An interactive tool for rapid biventricular analysis of congenital heart disease

Cardiac malformations are the most common birth defect. Better interventions in early life have improved mortality for children with congenital heart disease, but heart failure is a significant problem in adulthood. These patients require regular imaging and analysis of biventricular (left and right ventricular) function. In this study, we describe a rapid method to analyse left and right ventricular shape and function from cardiac MRI examinations. A 4D (3D+time) finite element model template is interactively customized to the anatomy and motion of the biventricular unit. The method was validated in 17 patients and 10 ex-vivo hearts. Interactive model updates were achieved through preconditioned conjugate gradient optimization on a multithread system, and by precomputing points predicted from a coarse mesh optimization.

Information maximizing component analysis of left ventricular remodeling due to myocardial infarction

BACKGROUND

Although adverse left ventricular shape changes (remodeling) after myocardial infarction (MI) are predictive of morbidity and mortality, current clinical assessment is limited to simple mass and volume measures, or dimension ratios such as length to width ratio. We hypothesized that information maximizing component analysis (IMCA), a supervised feature extraction method, can provide more efficient and sensitive indices of overall remodeling.

METHODS

IMCA was compared to linear discriminant analysis (LDA), both supervised methods, to extract the most discriminatory global shape changes associated with remodeling after MI. Finite element shape models from 300 patients with myocardial infarction from the DETERMINE study (age 31-86, mean age 63, 20 % women) were compared with 1991 asymptomatic cases from the MESA study (age 44-84, mean age 62, 52 % women) available from the Cardiac Atlas Project. IMCA and LDA were each used to identify a single mode of global remodeling best discriminating the two groups. Logistic regression was employed to determine the association between the remodeling index and MI. Goodness-of-fit results were compared against a baseline logistic model comprising standard clinical indices.

RESULTS

A single IMCA mode simultaneously describing end-diastolic and end-systolic shapes achieved best results (lowest Deviance, Akaike information criterion and Bayesian information criterion, and the largest area under the receiver-operating-characteristic curve). This mode provided a continuous scale where remodeling can be quantified and visualized, showing that MI patients tend to present larger size and more spherical shape, more bulging of the apex, and thinner wall thickness.

CONCLUSIONS

IMCA enables better characterization of global remodeling than LDA, and can be used to quantify progression of disease and the effect of treatment. These data and results are available from the Cardiac Atlas Project ( http://www.cardiacatlas.org ).

Real-time aortic pulse wave velocity measurement during exercise stress testing

BACKGROUND

Pulse wave velocity (PWV), a measure of arterial stiffness, has been demonstrated to be an independent predictor of adverse cardiovascular outcomes. This can be derived non-invasively using cardiovascular magnetic resonance (CMR). Changes in PWV during exercise may reveal further information on vascular pathology. However, most known CMR methods for quantifying PWV are currently unsuitable for exercise stress testing.

METHODS

A velocity-sensitive real-time acquisition and evaluation (RACE) pulse sequence was adapted to provide interleaved acquisition of two locations in the descending aorta (at the level of the pulmonary artery bifurcation and above the renal arteries) at 7.8 ms temporal resolution. An automated method was used to calculate the foot-to-foot transit time of the velocity pulse wave. The RACE method was validated against a standard gated phase contrast (STD) method in flexible tube phantoms using a pulsatile flow pump. The method was applied in 50 healthy volunteers (28 males) aged 22-75 years using a MR-compatible cycle ergometer to achieve moderate work rate (38 ± 22 W, with a 31 ± 12 bpm increase in heart rate) in the supine position. Central pulse pressures were estimated using a MR-compatible brachial device. Scan-rescan reproducibility was evaluated in nine volunteers.

RESULTS

Phantom PWV was 22 m/s (STD) vs. 26 ± 5 m/s (RACE) for a butyl rubber tube, and 5.5 vs. 6.1 ± 0.3 m/s for a latex rubber tube. In healthy volunteers PWV increased with age at both rest (R(2) = 0.31 p < 0.001) and exercise (R(2) = 0.40, p < 0.001). PWV was significantly increased at exercise relative to rest (0.71 ± 2.2 m/s, p = 0.04). Scan-rescan reproducibility at rest was -0.21 ± 0.68 m/s (n = 9).

CONCLUSIONS

This study demonstrates the validity of CMR in the evaluation of PWV during exercise in healthy subjects. The results support the feasibility of using this method in evaluating of patients with systemic aortic disease.

Challenges of cardiac image analysis in large-scale population-based studies

Large-scale population-based imaging studies of preclinical and clinical heart disease are becoming possible due to the advent of standardized robust non-invasive imaging methods and infrastructure for big data analysis. This gives an exciting opportunity to gain new information about the development and progression of heart disease across population groups. However, the large amount of image data and prohibitive time required for image analysis present challenges for obtaining useful derived data from the images. Automated analysis tools for cardiac image analysis are only now becoming available. This paper reviews the challenges and possible solutions to the analysis of big imaging data in population studies. We also highlight the potential of recent large epidemiological studies using cardiac imaging to discover new knowledge on heart health and well-being.

Quantification of LV function and mass by cardiovascular magnetic resonance: multi-center variability and consensus contours

BACKGROUND

High reproducibility of LV mass and volume measurement from cine cardiovascular magnetic resonance (CMR) has been shown within single centers. However, the extent to which contours may vary from center to center, due to different training protocols, is unknown. We aimed to quantify sources of variation between many centers, and provide a multi-center consensus ground truth dataset for benchmarking automated processing tools and facilitating training for new readers in CMR analysis.

METHODS

Seven independent expert readers, representing seven experienced CMR core laboratories, analyzed fifteen cine CMR data sets in accordance with their standard operating protocols and SCMR guidelines. Consensus contours were generated for each image according to a statistical optimization scheme that maximized contour placement agreement between readers.

RESULTS

Reader-consensus agreement was better than inter-reader agreement (end-diastolic volume 14.7 ml vs 15.2-28.4 ml; end-systolic volume 13.2 ml vs 14.0-21.5 ml; LV mass 17.5 g vs 20.2-34.5 g; ejection fraction 4.2 % vs 4.6-7.5 %). Compared with consensus contours, readers were very consistent (small variability across cases within each reader), but bias varied between readers due to differences in contouring protocols at each center. Although larger contour differences were found at the apex and base, the main effect on volume was due to small but consistent differences in the position of the contours in all regions of the LV.

CONCLUSIONS

A multi-center consensus dataset was established for the purposes of benchmarking and training. Achieving consensus on contour drawing protocol between centers before analysis, or bias correction after analysis, is required when collating multi-center results.

Big heart data: advancing health informatics through data sharing in cardiovascular imaging

The burden of heart disease is rapidly worsening due to the increasing prevalence of obesity and diabetes. Data sharing and open database resources for heart health informatics are important for advancing our understanding of cardiovascular function, disease progression and therapeutics. Data sharing enables valuable information, often obtained at considerable expense and effort, to be reused beyond the specific objectives of the original study. Many government funding agencies and journal publishers are requiring data reuse, and are providing mechanisms for data curation and archival. Tools and infrastructure are available to archive anonymous data from a wide range of studies, from descriptive epidemiological data to gigabytes of imaging data. Meta-analyses can be performed to combine raw data from disparate studies to obtain unique comparisons or to enhance statistical power. Open benchmark datasets are invaluable for validating data analysis algorithms and objectively comparing results. This review provides a rationale for increased data sharing and surveys recent progress in the cardiovascular domain. We also highlight the potential of recent large cardiovascular epidemiological studies enabling collaborative efforts to facilitate data sharing, algorithms benchmarking, disease modeling and statistical atlases.

Atlas-based quantification of cardiac remodeling due to myocardial infarction

Myocardial infarction leads to changes in the geometry (remodeling) of the left ventricle (LV) of the heart. The degree and type of remodeling provides important diagnostic information for the therapeutic management of ischemic heart disease. In this paper, we present a novel analysis framework for characterizing remodeling after myocardial infarction, using LV shape descriptors derived from atlas-based shape models. Cardiac magnetic resonance images from 300 patients with myocardial infarction and 1991 asymptomatic volunteers were obtained from the Cardiac Atlas Project. Finite element models were customized to the spatio-temporal shape and function of each case using guide-point modeling. Principal component analysis was applied to the shape models to derive modes of shape variation across all cases. A logistic regression analysis was performed to determine the modes of shape variation most associated with myocardial infarction. Goodness of fit results obtained from end-diastolic and end-systolic shapes were compared against the traditional clinical indices of remodeling: end-diastolic volume, end-systolic volume and LV mass. The combination of end-diastolic and end-systolic shape parameter analysis achieved the lowest deviance, Akaike information criterion and Bayesian information criterion, and the highest area under the receiver operating characteristic curve. Therefore, our framework quantitatively characterized remodeling features associated with myocardial infarction, better than current measures. These features enable quantification of the amount of remodeling, the progression of disease over time, and the effect of treatments designed to reverse remodeling effects.

Atlas-based Anatomical Modeling and Analysis of Heart Disease

Heart shape and function are major determinants of disease severity and predictors of future morbidity and mortality. Many studies now rely on non-invasive cardiac imaging techniques to quantify structural and functional changes. Statistical anatomical modeling of heart shape and motion provides a new tool for the quantification and evaluation of heart disease. This review surveys recent progress in the evaluation of statistical shape measures across populations and sub-cohorts, and highlights collaborative efforts to facilitate data sharing and atlas-based shape analysis.

Left ventricular shape variation in asymptomatic populations: the Multi-Ethnic Study of Atherosclerosis

BACKGROUND

Although left ventricular cardiac geometric indices such as size and sphericity characterize adverse remodeling and have prognostic value in symptomatic patients, little is known of shape distributions in subclinical populations. We sought to quantify shape variation across a large number of asymptomatic volunteers, and examine differences among sub-cohorts.

METHODS

An atlas was constructed comprising 1,991 cardiovascular magnetic resonance (CMR) cases contributed from the Multi-Ethnic Study of Atherosclerosis baseline examination. A mathematical model describing regional wall motion and shape was used to establish a coordinate map registered to the cardiac anatomy. The model was automatically customized to left ventricular contours and anatomical landmarks, corrected for breath-hold mis-registration between image slices. Mathematical techniques were used to characterize global shape distributions, after removal of translations, rotations, and scale due to height. Differences were quantified among ethnicity, sex, smoking, hypertension and diabetes sub-cohorts.

RESULTS

The atlas construction process yielded accurate representations of global shape (errors between manual and automatic surface points in 244 validation cases were less than the image pixel size). After correction for height, the dominant shape component was associated with heart size, explaining 32% of the total shape variance at end-diastole and 29% at end-systole. After size, the second dominant shape component was sphericity at end-diastole (13%), and concentricity at end-systole (10%). The resulting shape components distinguished differences due to ethnicity and risk factors with greater statistical power than traditional mass and volume indices.

CONCLUSIONS

We have quantified the dominant components of global shape variation in the adult asymptomatic population. The data and results are available at cardiacatlas.org. Shape distributions were principally explained by size, sphericity and concentricity, which are known correlates of adverse outcomes. Atlas-based global shape analysis provides a powerful method for quantifying left ventricular shape differences in asymptomatic populations.

TRIAL REGISTRATION

ClinicalTrials.gov NCT00005487.

A collaborative resource to build consensus for automated left ventricular segmentation of cardiac MR images

A collaborative framework was initiated to establish a community resource of ground truth segmentations from cardiac MRI. Multi-site, multi-vendor cardiac MRI datasets comprising 95 patients (73 men, 22 women; mean age 62.73±11.24years) with coronary artery disease and prior myocardial infarction, were randomly selected from data made available by the Cardiac Atlas Project (Fonseca et al., 2011). Three semi- and two fully-automated raters segmented the left ventricular myocardium from short-axis cardiac MR images as part of a challenge introduced at the STACOM 2011 MICCAI workshop (Suinesiaputra et al., 2012). Consensus myocardium images were generated based on the Expectation-Maximization principle implemented by the STAPLE algorithm (Warfield et al., 2004). The mean sensitivity, specificity, positive predictive and negative predictive values ranged between 0.63 and 0.85, 0.60 and 0.98, 0.56 and 0.94, and 0.83 and 0.92, respectively, against the STAPLE consensus. Spatial and temporal agreement varied in different amounts for each rater. STAPLE produced high quality consensus images if the region of interest was limited to the area of discrepancy between raters. To maintain the quality of the consensus, an objective measure based on the candidate automated rater performance distribution is proposed. The consensus segmentation based on a combination of manual and automated raters were more consistent than any particular rater, even those with manual input. The consensus is expected to improve with the addition of new automated contributions. This resource is open for future contributions, and is available as a test bed for the evaluation of new segmentation algorithms, through the Cardiac Atlas Project (www.cardiacatlas.org).

Atlas-based analysis of cardiac shape and function: correction of regional shape bias due to imaging protocol for population studies

BACKGROUND

Cardiovascular imaging studies generate a wealth of data which is typically used only for individual study endpoints. By pooling data from multiple sources, quantitative comparisons can be made of regional wall motion abnormalities between different cohorts, enabling reuse of valuable data. Atlas-based analysis provides precise quantification of shape and motion differences between disease groups and normal subjects. However, subtle shape differences may arise due to differences in imaging protocol between studies.

METHODS

A mathematical model describing regional wall motion and shape was used to establish a coordinate system registered to the cardiac anatomy. The atlas was applied to data contributed to the Cardiac Atlas Project from two independent studies which used different imaging protocols: steady state free precession (SSFP) and gradient recalled echo (GRE) cardiovascular magnetic resonance (CMR). Shape bias due to imaging protocol was corrected using an atlas-based transformation which was generated from a set of 46 volunteers who were imaged with both protocols.

RESULTS

Shape bias between GRE and SSFP was regionally variable, and was effectively removed using the atlas-based transformation. Global mass and volume bias was also corrected by this method. Regional shape differences between cohorts were more statistically significant after removing regional artifacts due to imaging protocol bias.

CONCLUSIONS

Bias arising from imaging protocol can be both global and regional in nature, and is effectively corrected using an atlas-based transformation, enabling direct comparison of regional wall motion abnormalities between cohorts acquired in separate studies.

Treatment with a copper-selective chelator causes substantive improvement in cardiac function of diabetic rats with left-ventricular impairment

BACKGROUND

Defective copper regulation is implicated as a causative mechanism of organ damage in diabetes. Treatment with trientine, a divalent-copper-selective chelator, improves arterial and renal structure/function in diabetes, wherein it also ameliorates left-ventricular (LV) hypertrophy. However, direct in vivo evidence that trientine can improve cardiac function in heart failure has hitherto been lacking.

METHODS

To determine whether trientine treatment could improve in vivo outcome, we measured cardiac function in groups of trientine-treated diabetic (TETA-DIA), non-drug-treated diabetic (DIA) and sham-treated control (SHAM) rats, by using in vivo high-field cardiac magnetic-resonance imaging (cMRI) and an ex vivo isolated-perfused working heart method. Forty age-matched animals underwent a cMRI scan after which 12 were randomized to the SHAM group and 28 underwent streptozotocin-injection; of these, 25 developed stable diabetes, and 12 were then randomized to receive no treatment for 16 weeks (DIA) and the other 13 to undergo 8-weeks' untreated diabetes followed by 8-weeks' drug treatment (TETA-DIA). Animals were studied again by cMRI at 8 and 16 weeks following disease induction, and finally by measurement of ex vivo cardiac function.

RESULTS

After eight weeks diabetes, rats (DIA/TETA-DIA) had developed significant impairment of LV function, as judged by impairment of ejection fraction (LVEF), cardiac output (CO), and LV mass (LVM)/body-mass (all P < 0.001), as well as other functional indexes. LVEF, CO (both P < 0.001) and the other indexes deteriorated further at 16 weeks in DIA, whereas trientine (TETA-DIA) improved cardiac function by elevating LVEF and CO (both P < 0.001), and also partially reversed the increase in LVM/body-mass (P < 0.05). In ex vivo hearts from DIA, the CO response to increasing preload pressure was deficient compared with SHAM (P < 0.001) whereas the preload-CO relationship was significantly improved in TETA-DIA animals (P < 0.001).

CONCLUSIONS

Trientine treatment significantly improved cardiac function in diabetic rats with substantive LV impairment. These results implicate impaired copper regulation in the pathogenesis of impaired cardiac function caused by diabetic cardiomyopathy, and support ongoing studies of trientine treatment in patients with heart failure.

Evaluation of left ventricular torsion by cardiovascular magnetic resonance

Recently there has been considerable interest in LV torsion and its relationship with symptomatic and pre-symptomatic disease processes. Torsion gives useful additional information about myocardial tissue performance in both systolic and diastolic function. CMR assessment of LV torsion is simply and efficiently performed. However, there is currently a wide variation in the reporting of torsional motion and the procedures used for its calculation. For example, torsion has been presented as twist (degrees), twist per length (degrees/mm), shear angle (degrees), and shear strain (dimensionless). This paper reviews current clinical applications and shows how torsion can give insights into LV mechanics and the influence of LV geometry and myocyte fiber architecture on cardiac function. Finally, it provides recommendations for CMR measurement protocols, attempts to stimulate standardization of torsion calculation, and suggests areas of useful future research.

Relationship between QRS duration and left ventricular mass and volume in patients at high cardiovascular risk

OBJECTIVE

Longer QRS duration on the ECG is associated with increased cardiovascular (CV) mortality. To evaluate potential mechanisms, we examined in this study the relationship between QRS duration and left ventricular (LV) mass and LV end systolic and end diastolic volume in patients with known CV disease or high-risk diabetes.

METHODS

In a substudy of the ONTARGET/TRANSCEND (Ongoing Telmisartan Alone and in Combination with Ramipril Global Endpoint Trial/Telmisartan Randomised Assessment Study in ACE Intolerant Subjects with Cardiovascular Disease) clinical trials, 368 patients had a cardiac magnetic resonance scan to measure LV mass, LV end systolic volume, LV end diastolic volume and LV ejection fraction at baseline and after 2 years of follow-up. Relationships between QRS duration on the 12-lead ECG and LV mass and volumes were evaluated at both assessments.

RESULTS

Each 10-ms increase in QRS duration both within and above the normal reference range was associated with an 8.3% (95% CI 6.7% to 9.9%) increase in LV mass, a 9.2% (95% CI 7.4% to 10%) increase in LV end diastolic volume and a 7.8% (95% CI 6.4% to 9.3%) increase in LV end systolic volume. QRS duration increased with body size, but associations with LV mass and volumes remained strong after indexing measurements to height(2.7) (p<0.001 for all) and were similar for subjects with an otherwise normal and abnormal ECG.

CONCLUSION

A longer QRS duration both within and above the normal reference range is associated with a greater LV mass and larger LV end systolic and end diastolic volumes. This may explain the known association of longer QRS duration with increased CV mortality.

The Cardiac Atlas Project--an imaging database for computational modeling and statistical atlases of the heart

Motivation: Integrative mathematical and statistical models of cardiac anatomy and physiology can play a vital role in understanding cardiac disease phenotype and planning therapeutic strategies. However, the accuracy and predictive power of such models is dependent upon the breadth and depth of noninvasive imaging datasets. The Cardiac Atlas Project (CAP) has established a large-scale database of cardiac imaging examinations and associated clinical data in order to develop a shareable, web-accessible, structural and functional atlas of the normal and pathological heart for clinical, research and educational purposes. A goal of CAP is to facilitate collaborative statistical analysis of regional heart shape and wall motion and characterize cardiac function among and within population groups.

Results: Three main open-source software components were developed: (i) a database with web-interface; (ii) a modeling client for 3D + time visualization and parametric description of shape and motion; and (iii) open data formats for semantic characterization of models and annotations. The database was implemented using a three-tier architecture utilizing MySQL, JBoss and Dcm4chee, in compliance with the DICOM standard to provide compatibility with existing clinical networks and devices. Parts of Dcm4chee were extended to access image specific attributes as search parameters. To date, approximately 3000 de-identified cardiac imaging examinations are available in the database. All software components developed by the CAP are open source and are freely available under the Mozilla Public License Version 1.1 (http://www.mozilla.org/MPL/MPL-1.1.txt).

Availability:http://www.cardiacatlas.org

Contact:a.young@auckland.ac.nz

Supplementary information:Supplementary data are available at Bioinformatics online.

Changes in mitral annular geometry and dynamics with ß-blockade in patients with degenerative mitral valve disease

BACKGROUND

remodeling of the mitral annulus contributes to progression of mitral regurgitation (MR). In patients with moderate-to-severe MR, short-term treatment with β-blockers has been shown to increase left ventricular (LV) end-diastolic and end-systolic volume, and this could deleteriously increase mitral valve annular dimensions. The objective of this study was to quantify the effects of a short duration of β-blocker treatment on mitral annular dimensions and dynamics in patients with MR due to primary degenerative valve disease.

METHODS AND RESULTS

twenty-five patients with moderate-to-severe degenerative MR and normal LV systolic function were studied in a double-blind crossover experiment using a β1-selective adrenergic blocker and placebo administered for 14±3 days. Cardiac MRI images were acquired after each treatment period to quantify mitral annular dimensions. At end diastole, there was no change in annular area (1659±331 versus 1632±299 mm(2); P<0.19), annular perimeter (154.3±16.4 versus 152±13.9 mm; P<0.13), septal-lateral (SL) dimension (38.0±5 versus 39.0±4.5 mm; P<0.15), or annular height (9.8±3.8 versus 9.5±2.5 mm; P<0.53). β-blockade resulted in significant end-diastole decreases in commissure-commissure dimension (48.9±4.6 versus 47.2±4.0 mm; P<0.01) and eccentricity (1.3±0.2 versus 1.2±0.1; P<0.01). At end systole (ES), β-blockade conferred a small, but significant decrease in annular perimeter (161.0±19.3 versus 156.8±16.9 mm; P<0.04) and eccentricity (1.2±0.1 versus 1.1±0.1; P<0.02), and the SL dimension significantly increased (41.5±5.7 versus 43.0±5.3 mm; P<0.03). Commissure-commissure dimension, annular area, and annular height at ES were not significantly different.

CONCLUSIONS

despite significant increases in LV end-diastolic and end-systolic volume, short-term β-blocker treatment of patients with moderate-to-severe MR reduced or preserved all mitral annular dimensions except SL at ES.

Left ventricular mass and volume with telmisartan, ramipril, or combination in patients with previous atherosclerotic events or with diabetes mellitus (from the ONgoing Telmisartan Alone and in Combination With Ramipril Global Endpoint Trial [ONTARGET])

The ONgoing Telmisartan Alone and in combination with Ramipril Global Endpoint Trial (ONTARGET) showed that the angiotensin receptor blocker telmisartan 80 mg was not inferior to the angiotensin-converting enzyme inhibitor ramipril 10 mg, and the combination no more effective than ramipril alone, in decreasing morbidity and mortality in patients with cardiovascular disease or high-risk diabetes. Although therapy targeting angiotensin II is known to decrease left ventricular (LV) mass and volume, the relative influence of angiotensin-converting enzyme inhibitor inhibitors and angiotensin receptor blocker, and their combination, on the heart remains unclear in this population. Magnetic resonance imaging was performed in 287 patients enrolled in ONTARGET, across 8 centers in 6 countries, at randomization and after 2-year treatment (90, 100, and 97 patients in the ramipril, telmisartan, and combination therapy groups, respectively). Baseline patient characteristics showed higher frequencies of coronary artery disease, Asian ethnicity, and use of statins and beta blockers than the main ONTARGET trial. LV mass decreased in all groups (p <0.0001 for each), but there were no significant differences in change in LV mass or volume among groups, except that LV mass index decreased more on combination versus telmisartan (p = 0.04). Key determinants of LV mass decrease were a history of hypertension (p = 0.03), baseline mass (p <0.0001), and decrease in systolic blood pressure (p <0.0001). The best magnetic resonance imaging predictor of composite events was end-systolic volume (p <0.0001). In conclusion, telmisartan and ramipril had similar effects on LV mass and volume, and combination therapy was not more effective, in high-risk patients with cardiovascular disease. These results are consistent with the major outcome findings of the main ONTARGET study.

Aortic valve stenotic area calculation from phase contrast cardiovascular magnetic resonance: the importance of short echo time

BACKGROUND

Cardiovascular magnetic resonance (CMR) can potentially quantify aortic valve area (AVA) in aortic stenosis (AS) using a single-slice phase contrast (PC) acquisition at valve level: AVA = aortic flow/aortic velocity-time integral (VTI). However, CMR has been shown to underestimate aortic flow in turbulent high velocity jets, due to intra-voxel dephasing. This study investigated the effect of decreasing intra-voxel dephasing by reducing the echo time (TE) on AVA estimates in patients with AS.

METHOD

15 patients with moderate or severe AS, were studied with three different TEs (2.8 ms/2.0 ms/1.5 ms), in the main pulmonary artery (MPA), left ventricular outflow tract (LVOT) and 0 cm/1 cm/2.5 cm above the aortic valve (AoV). PC estimates of stroke volume (SV) were compared with CMR left ventricular SV measurements and PC peak velocity, VTI and AVA were compared with Doppler echocardiography. CMR estimates of AVA obtained by direct planimetry from cine acquisitions were also compared with the echoAVA.

RESULTS

With a TE of 2.8 ms, the mean PC SV was similar to the ventricular SV at the MPA, LVOT and AoV0 cm (by Bland-Altman analysis bias +/- 1.96 SD, 1.3 +/- 20.2 mL/-6.8 +/- 21.9 mL/6.5 +/- 50.7 mL respectively), but was significantly lower at AoV1 and AoV2.5 (-29.3 +/- 31.2 mL/-21.1 +/- 35.7 mL). PC peak velocity and VTI underestimated Doppler echo estimates by approximately 10% with only moderate agreement. Shortening the TE from 2.8 to 1.5 msec improved the agreement between ventricular SV and PC SV at AoV0 cm (6.5 +/- 50.7 mL vs 1.5 +/- 37.9 mL respectively) but did not satisfactorily improve the PC SV estimate at AoV1 cm and AoV2.5 cm. Agreement of CMR AVA with echoAVA was improved at TE 1.5 ms (0.00 +/- 0.39 cm2) versus TE 2.8 (0.11 +/- 0.81 cm2). The CMR method which agreed best with echoAVA was direct planimetry (-0.03 cm2 +/- 0.24 cm2).

CONCLUSION

Agreement of CMR AVA at the aortic valve level with echo AVA improves with a reduced TE of 1.5 ms. However, flow measurements in the aorta (AoV 1 and 2.5) are underestimated and 95% limits of agreement remain large. Further improvements or novel, more robust techniques are needed in the CMR PC technique in the assessment of AS severity in patients with moderate to severe aortic stenosis.

Modelling passive diastolic mechanics with quantitative MRI of cardiac structure and function

The majority of patients with clinically diagnosed heart failure have normal systolic pump function and are commonly categorized as suffering from diastolic heart failure. The left ventricle (LV) remodels its structure and function to adapt to pathophysiological changes in geometry and loading conditions, which in turn can alter the passive ventricular mechanics. In order to better understand passive ventricular mechanics, a LV finite element (FE) model was customized to geometric data segmented from in vivo tagged magnetic resonance images (MRI) data and myofibre orientation derived from ex vivo diffusion tensor MRI (DTMRI) of a canine heart using nonlinear finite element fitting techniques. MRI tissue tagging enables quantitative evaluation of cardiac mechanical function with high spatial and temporal resolution, whilst the direction of maximum water diffusion in each voxel of a DTMRI directly corresponds to the local myocardial fibre orientation. Due to differences in myocardial geometry between in vivo and ex vivo imaging, myofibre orientations were mapped into the geometric FE model using host mesh fitting (a free form deformation technique). Pressure recordings, temporally synchronized to the tagging data, were used as the loading constraints to simulate the LV deformation during diastole. Simulation of diastolic LV mechanics allowed us to estimate the stiffness of the passive LV myocardium based on kinematic data obtained from tagged MRI. Integrated physiological modelling of this kind will allow more insight into mechanics of the LV on an individualized basis, thereby improving our understanding of the underlying structural basis of mechanical dysfunction under pathological conditions.

A copper(II)-selective chelator ameliorates left-ventricular hypertrophy in type 2 diabetic patients: a randomised placebo-controlled study

AIMS/HYPOTHESIS

Cu(II)-selective chelation with trientine ameliorates cardiovascular and renal disease in a model of diabetes in rats. Here, we tested the hypothesis that Cu(II)-selective chelation might improve left ventricular hypertrophy (LVH) in type 2 diabetic patients.

METHODS

We performed a 12 month randomised placebo-controlled study of the effects of treatment with the Cu(II)-selective chelator trientine (triethylenetetramine dihydrochloride, 600 mg given orally twice daily) on LVH in diabetic patients (n = 15/group at baseline) in an outpatient setting wherein participants, caregivers and those assessing outcomes were blinded to group assignment. Using MRI, we measured left ventricular variables at baseline, and at months 6 and 12. The change from baseline in left ventricular mass indexed to body surface area (LVM(bsa)) was the primary endpoint variable.

RESULTS

Diabetic patients had LVH with preserved ejection fraction at baseline. Trientine treatment decreased LVM(bsa) by 5.0 +/- 7.2 g/m(2) (mean +/- SD) at month 6 (when 14 trientine-treated and 14 placebo-treated participants were analysed; p = 0.0056 compared with placebo) and by 10.6 +/- 7.6 g/m(2) at month 12 (when nine trientine-treated and 13 placebo-treated participants were analysed; p = 0.0088), whereas LVM(bsa) was unchanged by placebo treatment. In a multiple-regression model that explained ~75% of variation (R (2) = 0.748, p = 0.001), cumulative urinary Cu excretion over 12 months was positively associated with trientine-evoked decreases in LVM(bsa).

CONCLUSIONS/INTERPRETATION

Cu(II)-selective chelation merits further exploration as a potential pharmacotherapy for diabetic heart disease.

TRIAL REGISTRATION

Australian New Zealand Clinical Trials Registry ACTRN 12609000053224 FUNDING: The Endocore Research Trust; Lottery Health New Zealand; the Maurice and Phyllis Paykel Trust; the Foundation of Research, Science and Technology (New Zealand); the Health Research Council of New Zealand; the Ministry of Education (New Zealand) through the Maurice Wilkins Centre for Molecular Biodiscovery; and the Protemix Corporation.