Mapping iron in human heart tissue with synchrotron x-ray fluorescence microscopy and cardiovascular magnetic resonance
© House et al.; licensee BioMed Central Ltd. 2014
Received: 25 June 2014
Accepted: 10 September 2014
Published: 27 September 2014
MRI assessment of cardiac iron is particularly important for assessing transfusion-dependent anaemia patients. However, comparing the iron distribution from histology or bulk samples to MRI is not ideal. Non-destructive, high-resolution imaging of post-mortem samples offers the ability to examine iron distributions across large samples at resolutions closer to those used in MRI. The aim of this ex vivo case study was to compare synchrotron X-ray fluorescence microscopy (XFM) elemental iron maps with magnetic resonance transverse relaxation rate maps of cardiac tissue samples from an iron-loaded patient.
Two 5 mm thick slices of formalin fixed cardiac tissue from a Diamond Blackfan anaemia patient were imaged in a 1.5 T MR scanner. R2 and R2* transverse relaxation rate maps were generated for both slices using RF pulse recalled spin echo and gradient echo acquisition sequences. The tissue samples were then imaged at the Australian Synchrotron on the X-ray Fluorescence Microscopy beamline using a focussed incident X-ray beam of 18.74 keV and the Maia 384 detector. The event data were analyzed to produce elemental iron maps (uncalibrated) at 25 to 60 microns image resolution.
The R2 and R2* maps and profiles for both samples showed very similar macro-scale spatial patterns compared to the XFM iron distribution. Iron appeared to preferentially load into the lateral epicardium wall and there was a strong gradient of decreasing iron, R2 and R2* from the epicardium to the endocardium in the lateral wall of the left ventricle and to a lesser extent in the septum. On co-registered images XFM iron was more strongly correlated to R2* (r = 0.86) than R2 (r = 0.79). There was a strong linear relationship between R2* and R2 (r = 0.87).
The close qualitative and quantitative agreement between the synchrotron XFM iron maps and MR relaxometry maps indicates that iron is a significant determinant of R2 and R2* in these ex vivo samples. The R2 and R2* maps of human heart tissue give information on the spatial distribution of tissue iron deposits.
Quantitative cardiovascular magnetic resonance imaging (CMR) techniques have been developed to non-invasively assess iron loading in the heart. CMR assessment of cardiac iron is particularly important for assessing transfusion-dependent anaemia patients, who can accumulate dangerous levels of iron in the heart that can eventually lead to heart failure. Chelation therapy has been successful in treating these patients, but directly monitoring heart iron loading is difficult.
Historically iron distribution has been measured and mapped using histological staining, which can provide some detailed information on iron distribution across small areas , but these stains are specific to the chemical speciation of iron and only semi-quantitative. Biochemical analysis of post-mortem heart samples suggests that there is a gradient in iron concentration from the epicardial to the endocardial layer ,. However, these bulk-sample techniques are relatively crude, destructive, and lack detailed spatial resolution. Biopsies are invasive and potentially unreliable because of the uneven distribution of iron in the heart.
To overcome the limitations inherent in biopsied material, quantitative CMR techniques have been developed to assess iron loading in the heart. These non-invasive techniques rely on measuring the proton transverse relaxation time or T2*. Typically a T2* image of a cross-section through the heart is generated and a region of interest in the septum selected to calculate an average T2* value. A low value of T2* (less than 20 ms) indicates iron overload in the heart tissue  and for patients with T2* less than 10 ms the risk of heart failure is significantly increased .
While the concentration of iron has a strong effect on T2*, it is not the sole determinant as the distribution and size of the iron deposits and susceptibility artefacts (e.g. proximity of blood vessels) can affect the relaxation time. The potential susceptibility artefacts in vivo limit the reliable region of the heart that can be assessed using T2* to the interventricular septum, a bridge of tissue separating the left and right ventricles ,.
Determining how closely CMR measurements represent the actual iron distribution in the human heart has been hindered by technical issues and artefacts associated with measuring the beating heart in vivo and a lack of quantitative elemental data at sufficient spatial resolution. CMR images are comprised of voxels, each voxel representing a sample volume of about 4 mm3 (ex-vivo), in comparison to the sample volumes of bulk biochemical methods, which are typically 150 times larger. Hence, relating CMR measurements to elemental analysis results, for example, is likely to be affected by partial volume sampling bias. In principle these errors can be minimised by making ex vivo CMR measurements on heart samples ,, but until recently obtaining elemental maps at a comparable resolution to CMR has been difficult.
Elemental mapping techniques based on synchrotron XFM have successfully imaged iron and other metals in human and rat brain tissue at resolutions approaching one micron -. We hypothesised that a similar elemental mapping technique applied to the heart, combined with ex vivo quantitative CMR measurements, would allow a much closer comparison and better understanding of the relationship between R2*, R2 and iron in the heart. We now present a single case history examination comparing iron maps generated by synchrotron XFM and quantitative CMR measurements.
The heart tissue for this case study came from a larger study of 12 patients investigating the amount and distribution of iron in the heart . The 22 year old male patient had Diamond Blackfan anaemia and died of cardiac failure and pneumonia. He was transfusion dependent and received chelation therapy from the age of 10, but compliance was poor. His mean left ventricular iron concentration was 3.91 mg/g dry weight . The study protocol was approved by all local research ethics committees, and local consent was obtained.
The heart was fixed in 10% neutral buffered formalin. Two slices of 5 mm thickness were cut in the ventricular short axis. The smaller slice was taken from a location closer to the heart apex and was labelled slice 7968. The larger slice was taken from a location closer to the mid-section and was labelled slice 7970.
CMR and analysis protocols
All CMR measurements were made on Siemens 1.5 T Avanto scanner (Siemens Medical Systems, Erlangen, Germany) using a phased-array body coil. R2 measurements were made with the tissue samples in air. A single spin echo sequence was used to acquire R2 data with seven echo-times between 6 and 30 ms, a TR of 750 ms, a flip angle of 90 degrees, and two signal averages. The field of view was 90 by 120 mm and the matrix was 192 by 256 pixels. For R2* measurements the heart slices were placed between Perspex plates and immersed in ultrapure water in such a way as to minimise the presence of bubbles. The temperature of the water containing the samples was 21.5°C. A multi-echo gradient echo sequence was used to acquire R2* data with eight echo-times (TE) between 2.47 and 18.95 ms, a repetition time (TR) of 22 ms, a flip angle of 35 degrees, and four signal averages. The field of view was 130 by 130 mm interpolated onto a 256 by 256 matrix. Coronal slices of 4 mm thickness were acquired parallel to the top and bottom surfaces of the tissue. Using high-resolution localiser images the slices were positioned such that the majority of the slice was within the margin of the tissue.
to the data, where S(0) is the initial amplitude. Voxels were excluded where the goodness of fit parameter (R2) was less than 0.9. Maps of R2 were filtered to reduce noise by using a median filter with a radius of two pixels.
Synchrotron sample preparation and imaging
XFM, R2 and R2* images for sample 7970 were co-registered at the same resolution (0.5 × 0.5 mm) using ImageJ (NIH) following resampling of the high resolution XFM image to match the coarser resolution of the CMR images. To quantitatively compare the XFM and relaxometry data, 114 circular regions of interest were overlayed on the left ventricle of the co-registered XFM, R2 and R2* maps. The R2 and R2* values were plotted against the XFM intensity and assessed by the Pearson correlation coefficient.
The wet weight iron concentrations for the two heart slices scanned in this study were extracted from ICP-AES data acquired in a previous study on adjoining slices from the same heart . The XFM iron maps were digitally overlayed onto the photographs of each sample using Photoshop CS4 (Adobe) and the approximate segments used in Carpenter et al. study  were digitally marked on images and annotated with the measured iron concentrations.
XFM iron maps
In the interventricular septum the iron is distributed differently, particularly for slice 7968 where two discontinuous bands of higher iron are apparent (Figure 2B) within the myocardium layer. In heart slice 7970 there appears to be one dominant band of elevated iron signal, also within the myocardium (Figure 3B). The XFM iron signal is, in general, more uniform through the septum. The right ventricle of both samples shows a discontinuous band of high iron signal within the epicardium layer.
In heart slice 7968 (Figure 2B, black arrows), two oval regions of low iron signal are evident in the inferior myocardium and the anteroseptal epicardium. On the photograph of the sample (Figure 2A, black arrows), the anomaly in the inferior myocardium corresponds to a lighter coloured region of tissue about 3–4 mm across. In heart slice 7970, there is distinct linear structure that transects the inferior wall of the heart through the myocardium and endocardium layers that corresponds to a low iron signal (Figure 3B, black arrow) and a light coloured feature in the photograph of the sample (Figure 3A, black arrow). The papillary muscles also appeared to contain moderate to high iron concentrations.
Comparison between XFM Iron Maps and CMR R2 and R2* Maps
XFM versus elemental iron concentrations
To our knowledge this is the first demonstration of high-resolution synchrotron XFM elemental imaging of an entire cross section of human heart tissue and the first study to compare the distribution of XFM iron maps to quantitative CMR R2 and R2* images. The fast scanning approach coupled with the Maia detector allowed the acquisition of XFM iron maps that displayed the internal iron distribution of the heart in great detail. Both the R2 and R2* maps displayed a high level of correspondence to the macro-scale distribution of iron in the XFM images and profiles. On co-registered images, R2 and R2* also showed strong linear correlations with the XFM signal.
The gradient of decreasing iron from the epicardium to endocardium we observed in XFM and transverse relaxation rates is consistent with the observations of Carpenter et al.  and Ghugre et al. . Compilation of the iron concentrations from all the samples from this heart (Additional file 1: Figure S1) also confirm that mean iron concentrations in the epicardium are significantly higher compared to the mean iron concentrations in the myocardium and endocardium. It is not clear why iron appears to preferentially load into the lateral epicardium of this patient. Consistent with our case study, an in-vivo study on Thalassemia Major and Sickle Cell Disease patients by Ghugre et al.  observed elevated R2* values in the lateral segment of the heart for patients with T2* greater than 20 ms compared to patients with T2* less than 20 ms. In contrast, the case-study of a 24 year old thalassaemic heart by Ghugre et al.  measured similar iron concentrations to our study, but found the outer septal layer had the highest mean iron concentration, while the lateral epicardium had the lowest iron concentration in the left ventricle epicardium layer, as measured by elemental analysis. Interestingly, the lateral epicardium in that study had the highest mean R2* and R2, consistent with our case study. Other studies have reported segmental variation in R2* for TM patients  consistent with the variations we observed in inferior and anterior segments. Notwithstanding the results of these workers, observations from this single case study on the distribution of iron in heart tissue may not necessarily be generalizable to other patients with cardiac iron overload.
Although a gradient in XFM iron exists in the septum of both samples, it is not as strong or consistent as in the lateral parts of the left ventricle. The images and profiles suggest that, at least for this case study, the distribution of iron, R2 and R2* in the septal wall is more uniform overall with less extremes in the range of iron concentrations, particularly in the outer layer. The lack of a strong iron gradient in the septum potentially relates to the absence of an epicardium layer in the septum, which instead comprises a left ventricle endocardium, a mesocardium and a right ventricle endocardium . However, there is still discernable structure evident in the septal region of both hearts (Additional file 2: Figure S2).
Qualitatively and quantitatively, the R2* data agrees more closely with the XFM iron data than the R2 data. Apart from the greater sensitivity of R2* to iron concentrations, as suggested by the larger regression slope between R2* and XFM iron, there may be other reasons for the better agreement of the R2* data. The shorter data acquisition time of the gradient echo sequence allowed 4 signals to be averaged for each gradient echo, compared to 2 averages for the spin echo data. Hence, the R2* data may have benefited from higher signal to noise compared to the R2 data. In addition, R2 could potentially to be more affected by non-iron mediated relaxation mechanisms (e.g. hydration) at the iron concentrations observed in this case study.
He et al.  investigated in-vivo cardiac R2 and R2* and observed that the relationship between R2 and R2* is curvilinear. Our in-vitro data on fixed cardiac tissue, and that of Ghugre et al.  from fresh cardiac tissue, both show strong linear relationships between R2 and R2* over a comparable range in R2* values to He et al. . While differences may arise comparing in-vivo to in-vitro measurements, another possible explanation for the discrepancy is the difference in R2 acquisition sequences. While the multi-echo gradient echo sequences between all three studies are similar, the R2 measurements in our study and Ghugre et al.  used a single spin echo acquisition where as He et al.  were restricted to a fast multi-echo spin echo acquisition to enable a short breath-hold acquisition in vivo. A multi-echo spin echo acquisition will potentially reduce the measured R2 relative to a single spin echo measurement as the multi-echo refocusing pulses suppress the effects of diffusion. Such a reduction in R2 may explain the curvilinear relationship observed by He et al. . Alternatively, the simple fitting of all echo times we used in this study could, in the presence of noise, underestimate high R2* values potentially masking a more curvilinear relationship. In regard to the curvilinear relationship between R2 and R2* reported in the liver , we note that at relatively low iron concentrations, such as in this case study (maximum iron concentration 6 mg/g dry weight) and in Ghugre et al. , the relationship between R2 and R2* is essentially linear.
It is difficult to extrapolate the CMR observations from this in-vitro study to the in-vivo measurement of cardiac T2*. Our study was performed on formalin fixed tissue and we reduced the magnetic susceptibility artefacts from interfaces found between living tissue and air or blood by acquiring our R2* data with the sample immersed in water. In-vitro imaging also does not have to contend with the motion of the beating heart. The septal region of the heart is typically used in cardiac T2* imaging for iron assessment , as it is bounded by blood on both sides and hence susceptibility variations are potentially minimised . Assuming the in-vivo and in-vitro septum are likely to be the most similar in terms of their R2* response, our results are consistent with R2* over the whole septum being a reasonable proxy for heart iron concentration in this region.
The ability to map elements in large tissue specimens using a synchrotron source can highlight new insights into fundamental biological systems and importantly provides a non-destructive technique for comparison and validation of other imaging modalities. The close qualitative and quantitative agreement between the synchrotron XFM iron maps and CMR relaxometry maps indicates that iron is a significant contributor to R2 and R2* in these ex vivo samples and that R2 and R2* maps of human heart tissue give information on the spatial distribution of tissue iron deposits.
The study protocol was approved by all local research ethics committees, and local consent was obtained.
MH, TSP: experimental conception and design, acquisition, analysis and interpretation of data, drafting of manuscript; AF, MDJ, DP, DH: experimental design, acquisition, analysis and interpretation of data, critical and intellectual revision of manuscript; JPC, DP: sample acquisition, critical and intellectual revision of manuscript. All authors read and approved the final manuscript.
We are indebted to the generosity of the patients and families without whose help we would have been unable to perform this important study. We would also like to thank SKG Radiology for their radiological support. This research was undertaken on the X-ray Fluorescence Microscopy beamline at the Australian Synchrotron, Victoria, Australia.
- Barosi G, Arbustini E, Gavazzi A, Grasso M, Pucci A: Myocardial iron grading by endomyocardial biopsy. A clinico-pathologic study on iron overloaded patients. Eur J Haematol. 1989, 42: 382-388. 10.1111/j.1600-0609.1989.tb01229.x.View ArticlePubMedGoogle Scholar
- Buja LM, Roberts WC: Iron in the heart: Etiology and clinical significance. Am J Med. 1971, 51: 209-221. 10.1016/0002-9343(71)90240-3.View ArticlePubMedGoogle Scholar
- Ghugre NR, Enriquez CM, Gonzalez I, Nelson MD, Coates TD, Wood JC: MRI detects myocardial iron in the human heart. Magn Reson Med. 2006, 56: 681-686. 10.1002/mrm.20981.PubMed CentralView ArticlePubMedGoogle Scholar
- Pennell DJ: T2* Magnetic Resonance and Myocardial Iron in Thalassemia. Ann N Y Acad Sci. 2005, 1054: 373-378. 10.1196/annals.1345.045.View ArticlePubMedGoogle Scholar
- Kirk P, Roughton M, Porter JB, Walker JM, Tanner MA, Patel J, Wu D, Taylor J, Westwood MA, Anderson LJ, Pennell DJ: Cardiac T2* Magnetic Resonance for Prediction of Cardiac Complications in Thalassemia Major. Circulation. 2009, 120: 1961-1968. 10.1161/CIRCULATIONAHA.109.874487.PubMed CentralView ArticlePubMedGoogle Scholar
- Anderson LJ, Holden S, Davis B, Prescott E, Charrier CC, Bunce NH, Firmin DN, Wonke B, Porter J, Walker JM, Pennell DJ: Cardiovascular T2-star (T2*) magnetic resonance for the early diagnosis of myocardial iron overload. Eur Heart J. 2001, 22: 2171-2179. 10.1053/euhj.2001.2822.View ArticlePubMedGoogle Scholar
- Wood JC, Tyszka JM, Carson S, Nelson MD, Coates TD: Myocardial iron loading in transfusion-dependent thalassemia and sickle cell disease. Blood. 2004, 103: 1934-1936. 10.1182/blood-2003-06-1919.View ArticlePubMedGoogle Scholar
- Carpenter J-P, He T, Kirk P, Roughton M, Anderson LJ, de Noronha SV, Sheppard MN, Porter JB, Walker JM, Wood JC, Galanello R, Forni G, Catani G, Matta G, Fucharoen S, Fleming A, House MJ, Black G, Firmin DN, St. Pierre TG, Pennell DJ: On T2* Magnetic Resonance and Cardiac Iron. Circulation. 2011, 123: 1519-1528. 10.1161/CIRCULATIONAHA.110.007641.PubMed CentralView ArticlePubMedGoogle Scholar
- Popescu BFG, George MJ, Bergmann U, Garachtchenko AV, Kelly ME, McCrea RPE, Luning K, Devon RM, George GN, Hanson AD, Harder SM, Chapman LD, Pickering IJ, Nichol H: Mapping metals in Parkinson’s and normal brain using rapid-scanning x-ray fluorescence. Phys Med Biol. 2009, 54: 651-663. 10.1088/0031-9155/54/3/012.View ArticlePubMedGoogle Scholar
- McCrea RPE, Harder SL, Martin M, Buist R, Nichol H: A comparison of rapid-scanning X-ray fluorescence mapping and magnetic resonance imaging to localize brain iron distribution. Eur J Radiol. 2008, 68: S109-S113. 10.1016/j.ejrad.2008.04.048.View ArticlePubMedGoogle Scholar
- Hopp K, Popescu BFG, McCrea RPE, Harder SL, Robinson CA, Haacke ME, Rajput AH, Rajput A, Nichol H: Brain iron detected by SWI high pass filtered phase calibrated with synchrotron X-ray fluorescence. J Magn Reson Imaging. 2010, 31: 1346-1354. 10.1002/jmri.22201.View ArticlePubMedGoogle Scholar
- Zheng W, Nichol H, Liu S, Cheng Y-CN, Haacke EM: Measuring iron in the brain using quantitative susceptibility mapping and X-ray fluorescence imaging. NeuroImage. 2013, 78: 68-74. 10.1016/j.neuroimage.2013.04.022.View ArticlePubMedGoogle Scholar
- Zheng W, Haacke EM, Webb SM, Nichol H: Imaging of stroke: a comparison between X-ray fluorescence and magnetic resonance imaging methods. Magn Reson Imaging. 2012, 30: 1416-1423. 10.1016/j.mri.2012.04.011.View ArticlePubMedGoogle Scholar
- James S, Myers D, Jonge M, Vogt S, Ryan C, Sexton B, Hoobin P, Paterson D, Howard D, Mayo S, Altissimo M, Moorhead G, Wilkins S: Quantitative comparison of preparation methodologies for x-ray fluorescence microscopy of brain tissue. Anal Bioanal Chem. 2011, 401: 853-864. 10.1007/s00216-011-4978-3.View ArticlePubMedGoogle Scholar
- Paterson D, de Jonge MD, Howard DL, Lewis W, McKinlay J, Starritt A, Kusel M, Ryan CG, Kirkham R, Moorhead G, Siddons DP: The X-ray Fluorescence Microscopy Beamline at the Australian Synchrotron. AIP Conf Proc. 2011, 1365: 219-222. 10.1063/1.3625343.View ArticleGoogle Scholar
- Kirkham R, Dunn PA, Kuczewski AJ, Siddons DP, Dodanwela R, Moorhead GF, Ryan CG, De Geronimo G, Beuttenmuller R, Pinelli D, Pfeffer M, Davey P, Jensen M, Paterson DJ, de Jonge MD, Howard DL, Küsel M, McKinlay J: The Maia Spectroscopy Detector System: Engineering for Integrated Pulse Capture, Low-Latency Scanning and Real-Time Processing. AIP Conf Proc. 2010, 1234: 240-243. 10.1063/1.3463181.View ArticleGoogle Scholar
- Ryan CG, Etschmann BE, Vogt S, Maser J, Harland CL, van Achterbergh E, Legnini D: Nuclear microprobe – synchrotron synergy: Towards integrated quantitative real-time elemental imaging using PIXE and SXRF. Nucl Instrum Methods Physics Res Sec B: Beam Interact Mater Atoms. 2005, 231: 183-188. 10.1016/j.nimb.2005.01.054.View ArticleGoogle Scholar
- Ghugre NR, Enriquez CM, Coates TD, Nelson MD, Wood JC: Improved R2* measurements in myocardial iron overload. J Magn Reson Imaging. 2006, 23: 9-16. 10.1002/jmri.20467.PubMed CentralView ArticlePubMedGoogle Scholar
- Positano V, Pepe A, Santarelli MF, Ramazzotti A, Meloni A, De Marchi D, Favilli B, Cracolici E, Midiri M, Spasiano A, Lombardi M, Landini L: Multislice multiecho T2* cardiac magnetic resonance for the detection of heterogeneous myocardial iron distribution in thalassaemia patients. NMR Biomed. 2009, 22: 707-715. 10.1002/nbm.1382.View ArticlePubMedGoogle Scholar
- Greenbaum RA, Ho SY, Gibson DG, Becker AE, Anderson RH: Left ventricular fibre architecture in man. Br Heart J. 1981, 45: 248-263. 10.1136/hrt.45.3.248.PubMed CentralView ArticlePubMedGoogle Scholar
- He T, Smith GC, Gatehouse PD, Mohiaddin RH, Firmin DN, Pennell DJ: On using T2 to assess extrinsic magnetic field inhomogeneity effects on T2* measurements in myocardial siderosis in thalassemia. Magn Reson Med. 2009, 61: 501-506. 10.1002/mrm.21874.View ArticlePubMedGoogle Scholar
- Wood JC, Enriquez C, Ghugre N, Tyzka JM, Carson S, Nelson MD, Coates TD: MRI R2 and R2* mapping accurately estimates hepatic iron concentration in transfusion-dependent thalassemia and sickle cell disease patients. Blood. 2005, 106: 1460-1465. 10.1182/blood-2004-10-3982.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.