Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2013; 8(8): e72795.
Published online 2013 August 26. doi:  10.1371/journal.pone.0072795
PMCID: PMC3753231

Diffusion MRI Tractography of the Developing Human Fetal Heart

Wolfgang Rudolf Bauer, Editor



Human myocardium has a complex and anisotropic 3D fiber pattern. It remains unknown, however, when in fetal life this anisotropic pattern develops and whether the human heart is structurally fully mature at birth. We aimed here to use diffusion tensor MRI (DTI) tractography to characterize the evolution of fiber architecture in the developing human fetal heart.


Human fetal hearts (n = 5) between 10–19 weeks of gestation were studied. The heart from a 6-day old neonate and an adult human heart served as controls. The degree of myocardial anisotropy was measured by calculating the fractional anisotropy (FA) index. In addition, fiber tracts were created by numerically integrating the primary eigenvector field in the heart into coherent streamlines.


At 10–14 weeks the fetal hearts were highly isotropic and few tracts could be resolved. Between 14–19 weeks the anisotropy seen in the adult heart began to develop. Coherent fiber tracts were well resolved by 19 weeks. The 19-week myocardium, however, remained weakly anisotropic with a low FA and no discernable sheet structure.


The human fetal heart remains highly isotropic until 14–19 weeks, at which time cardiomyocytes self-align into coherent tracts. This process lags 2–3 months behind the onset of cardiac contraction, which may be a prerequisite for cardiomyocyte maturation and alignment. No evidence of a connective tissue scaffold guiding this process could be identified by DTI. Maturation of the heart’s sheet structure occurs late in gestation and evolves further after birth.


The myocardium is an anisotropic tissue in which the constituent cardiomyocytes (CMs) are arranged in a complex 3D network [1]. Groups of CMs in this network form laminar sheets, which allow the myocardium to shear and thicken during systole [2]. While the myocardium is a continuum, it is instructive to group aligned CMs into myofiber tracts. This approach has been used by anatomists, histologists, physicists and cardiologists to study myocardial microstructure in both health and disease [3][10].

In the adult myocardium, the alignment of the myocytes is highly anisotropic, creating a transmural gradient in fiber orientation [1], [7]. Fibers in the subendocardium have a positive helix angle (HA) and those in the subepicardium a negative HA, while fibers in the midmyocardium are circumferential [1], [7]. It remains unknown, however, when in fetal life and by what mechanisms this architectural pattern develops. The aim of the current study was, therefore, to use diffusion tensor magnetic resonance imaging (DTI) tractography to characterize the development of myofiber architecture in the human fetal heart. We hypothesize that this knowledge has the potential to elucidate the pathogenesis of a broad range of congenital heart defects and also to provide a valuable template for the regeneration of functional myocardium in the adult.

Materials and Methods

Ethics Statement: The Cardiac Registry collection at Boston Children’s Hospital was started in the 1940s and includes a large number of postmortem hearts, including fetal hearts, for teaching and research. Written consent has been obtained from the parents/guardians at the time of donation to the Registry (according to the hospital policy in force at the time) stating that the hearts can be used for research purposes. However, prior to the last two decades hearts were sometimes donated anonymously to the registry for teaching and research with no personal information or identifiers. All the hearts in this study were completely anonymized prior to imaging. A research protocol to perform this study was specifically approved by the Institutional Review Board of the Massachusetts General Hospital (MGH). The hearts were imaged at the MGH in a completely anonymized fashion. Only the gestational age of the hearts was known.

Six normal human hearts were obtained from the cardiac registry collection at Children’s Hospital in Boston. The gestational ages of the five fetal hearts studied were: 10, 13, 14, 14 and 19 weeks. The sixth heart, a normal neonatal heart, was obtained 6 days post birth (P6). All hearts were fixed in formaldehyde and were structurally intact. DTI of the fetal hearts was performed at 9.4 Tesla (Biospec, Bruker, Billerica MA) with a 3D fat-suppressed single-shot spin echo EPI (echoplanar) sequence with the following parameters: spatial resolution 130×130×130 µm3, TR = 1000 ms, TE = 18 ms, b-value of 1400–1600 s/mm2, 24 diffusion-encoding directions plus one b = 0 s/mm2 image. DTI of the P6 (day 6 neonatal heart) was performed using the identical sequence at 4.7 Tesla (Biospec, Bruker) with the following parameters: spatial resolution 500×500×500 µm3, TR = 1200 ms, TE = 41 ms, b-value of 1560 s/mm2, 24 diffusion-encoding directions plus one b = 0 s/mm2 image.

The fetal and P6 hearts were compared with each other and an adult human heart, which was imaged as previously described [10]. In brief, a 2D fat-suppressed single-shot spin echo EPI sequence was used on a 3 Tesla clinical imaging system (TRIO, Siemens, Erlangen Germany). Parameters included: voxel size = 2×2×2 mm3, TR/TE = 8430/96 ms, 50–70 short axis slices (without gaps) to cover the entire heart, b-value of 2000 s/mm2, 6 diffusion-encoding directions plus one b = 0 s/mm2 image, 24 averages.

In all hearts, diagonalization of the diffusion tensor in each voxel was performed to derive the primary, secondary and tertiary eigenvectors (ê1, ê2, and ê3) and their associated eigenvalues (λ1, λ2, λ3). The image acquisition parameters above were chosen to ensure that the transmural thickness of the myocardium in each heart would contain more than 8 voxels. The structural maturation of the fetal myocardium was evaluated by measuring its mean diffusivity (MD) and the degree of tissue anisotropy, estimated by the fractional anisotropy (FA) index [11], [12]. Further characterization of the diffusion tensor field within the myocardium was performed using an ellipsoid-based representation. In addition, the primary eigenvector (ê1) field in the lateral wall of the left ventricle (LV) was numerically integrated into streamlines using a 4th order Runge-Kutta approach to create myofiber tracts [10].

Following MR image acquisition, the fetal hearts were sectioned in a plane tangential to the epicardium in the lateral wall. Phase-contrast microscopy was performed on H&E stained histological sections of the hearts at 20×magnification. Values in the manuscript are reported as the mean ± standard deviation of all voxels in a given heart. Comparison of the MD and FA values in the fetal, P6 and adult hearts was performed with ANOVA and a Tukey’s test.


Fiber architecture in the adult LV is shown in Figure 1. The subendocardial and subepicardial fibers cross over each other at an angle of approximately 120o. The ellipsoidal glyphs form a coherent and ordered pattern throughout the LV and have a prolate shape representing a high degree of structural anisotropy.

Figure 1
Fiber architecture in the adult human heart.

DTI of the human fetal hearts showed that myocardial anisotropy develops fairly late in fetal life. At 10 weeks of gestation, the ellipsoids in the myocardium were spherical (all eigenvalues similar in magnitude), consistent with a largely isotropic tissue (Figure 2A–B). At 14 weeks, the ellipsoids became oblate and oriented in a transmural pattern resembling those seen in the adult heart (Figure 2C). By 19 weeks of gestation, the ellipsoidal glyphs began to assume a more prolate shape as anisotropy in the myocardium began to increase (Figure 2D).

Figure 2
Transition from tissue isotropy to anisotropy in the developing human fetal heart.

FA in the hearts was first computed within a region-of-interest in the lateral wall and then in the entire left ventricle. FA (mean ± SD) in the lateral wall at 10 weeks of gestation was 0.09±0.01, remained at 0.09±0.01 at 14 weeks, and increased to 0.13±0.03 (p<0.01) at 19 weeks. By P6, FA in the lateral wall had increased to 0.27±0.05 (p<0.01). A corresponding decrease in MD from 1.34±0.09 at 19 weeks to 0.9±0.09 (p<0.01) at P6 was seen. The transition from isotropy to anisotropy thus started between 14–19 weeks of gestation but was a slow process, and was still incomplete at birth. FA in the lateral wall remained significantly lower in the P6 heart (0.27±0.05) than in the adult heart (0.35±0.04, p<0.01).

FA values derived from the entire left ventricle are summarized in Table 1. A steady increase was seen in FA from 14 weeks of gestation (0.07±0.02) to 19 weeks (0.1±0.04), P6 (0.2±0.08) and adulthood (0.31±0.06). Individual eigenvalues and their ratios are listed in Table 1. Little change was seen in λ2 and λ3 between 10–19 weeks. However, a 46% decrease in λ2 and a 51% decrease in λ3 were seen between the 19-week and adult hearts. The ratios of λ12 and λ13 remained low at 19 weeks (1.12±0.07 and 1.21±0.09, respectively) indicating that myocardial sheet structure had not yet matured. These ratios were significantly higher in the adult heart (1.53±0.17 and 1.83±0.23, respectively, p<0.01). A similar pattern of structural development was thus observed in the lateral wall and the entire left ventricle. Of note, significant differences in FA, the eigenvalues and their ratios were seen between the P6 and adult heart (Table 1). The λ12 and λ13 ratios in the P6 heart (1.23±0.13 and 1.53±0.5, respectively) were significantly lower than the adult heart (p<0.01), indicating that myocardial sheet architecture is not fully mature at birth.

Table 1
Fractional anisotropy and eigenvalues in the entire left ventricle (mean ± SD).

FA maps and histograms of a 10-week, 14-week, 19-week, P6 and adult heart are shown in Figure 3. The maps and histograms confirm that myofiber organization begins to develop between 14–19 weeks of gestation and continues to evolve after birth. Myofiber tractography in the lateral LV wall revealed a similar pattern of development (Figure 4). At 10 weeks of gestation few fiber tracts could be reconstructed, consistent with tissue isotropy (Figure 4A). At 14 weeks more tracts could be reconstructed but their length and density remained well below those observed in adult myocardium (Figure 4B). By 19 weeks of gestation the characteristic crossing helical pattern seen in adult hearts could be well resolved (Figure 4C, 4E). The 19-week myocardium, however, remained a weakly anisotropic tissue. While sufficient anisotropy had developed to form coherent tracts, the density of the tracts was low and the myocardium was still characterized by a very low FA and a high MD.

Figure 3
Fractional anisotropy (FA) in the developing heart.
Figure 4
Tractography of fiber tracts (color-coded by HA) passing through a region-of-interest in the lateral LV wall.

Histology of the lateral wall of the LV (Figure 5) confirmed the DTI findings. At 10 weeks, the fetal myocardium was highly isotropic with no discernable architectural pattern (Figure 5A–C). At 14 weeks, coherent patterns of myocyte alignment could be seen in some parts of the myocardium (Figure 5D–F). The characteristic transmural evolution of CM alignment was consistently seen at 19 weeks (Figure 5G–I).

Figure 5
Histological sections of the lateral LV wall.


The evaluation of myocardial architecture in the human fetal heart has been performed using histology and microscopy [13], [14]. These techniques, however, involve the destruction of the heart, are subject to tissue shearing/tearing and are two-dimensional. Here, for the first time we evaluate the microstructure of the human fetal heart with DTI. This allowed the hearts to be imaged in a completely intact state, 3D tracts of myofiber architecture to be resolved and measures of diffusion anisotropy to be quantified.

Our data suggest that fiber architecture in the human fetal heart develops in two distinct phases. At 10 weeks of gestation the myocardium is isotropic with no coordinated alignment of its CMs. Between 14 to 19 weeks of gestation the CMs align themselves into resolvable tracts that resemble those seen in the adult heart. However, in addition to FA, the λ12 and λ13 ratios remain low indicating that water diffusion has not yet become restricted by densely packed myofiber sheets. Sheet architecture, which is crucial to the function of a fully mature adult heart, thus develops relatively late in fetal life.

The evolution of fiber architecture in the fetal heart remains poorly understood and characterized. Studies in mice have shown that fetal CMs are spherical, and thus isotropic, until Carnegie Stage 18 of gestation [15]. At this time they begin to elongate, hypertrophy and form end-end connections supporting the generation of myofiber tracts. This process, however, is slow and is not yet complete at birth [15]. The results of our study suggest that a similarly slow process occurs in the human fetal heart. Myofiber anisotropy, and tracts resembling those seen in the adult heart, developed between 14–19 weeks of gestation, corresponding to the initiation of CM elongation in the human fetal heart. However, FA at 19 weeks remained extremely low, indicating that the CMs had only undergone a mild degree of elongation [11], [12].

Animal studies have shown that the cytoarchitecture of the heart is not mature at birth [15]. Further CM elongation and hypertrophy continue for at least the first month of post-natal life. Likewise, the fiber architecture of the neonatal pig heart has recently been shown to be highly plastic [16]. These findings are consistent with our observation that FA and the ratios of λ12 and λ13 in the P6 human heart were significantly lower than in the adult heart. The human heart thus likely also undergoes further structural maturation within the first months of life. This may have important implications for the timing of therapy in congenital heart disease.

Prior to the elongation of its constituent CMs the fetal heart undergoes four important processes: looping, trabeculation, septation and myocardial compaction [13]. These processes, however, are completed well before 10 weeks of gestation and are therefore unlikely to directly affect myocardial anisotropy [13]. Rather, it seems most likely that the process of CM elongation and polarization, which is most prominent in the latter stages of gestation [15], produces the 3D pattern of CMs seen in the adult heart. The signals responsible for this have not been fully elucidated, but likely involve the noncanonical Wnt planar cell polarity (PCP) pathway [17], [18], and possibly the ERK5 pathway [15].

The results of this study provide several intriguing insights into the structure-function relationship in the human fetal heart. Most intriguingly, the development of coherent myofiber tracts lagged behind the onset of cardiac contraction by 2–3 months [13]. It should be noted, however, that the extremely low impedance of the placenta provides a very low load against which the fetal heart contracts. Altering this load experimentally impacts CM orientation in the fetal heart [19], suggesting that mechanical forces may also play a role in CM alignment.

No evidence of a preexisting connective tissue scaffold with the ability to restrict water diffusion, and accordingly be resolvable with DTI, was seen in our study. This is consistent with prior observations that the deposition of the laminin network in the heart coincides with, and is a consequence of, CM elongation in the latter part of gestation [15]. This suggests that complex tissue scaffolds may not be a prerequisite for successful stem cell therapy in the heart. However, the developmental timeline of fiber tracts in the fetal heart suggests that a prolonged period (2–3 months) of preexisting cyclic contraction may be required to condition primitive CMs to mature, elongate and align.

The procurement of completely intact (non-dissected) human fetal hearts is extremely challenging and limited the number of hearts imaged in this study. Our findings, however, are in complete agreement with prior observations made by histology and microscopy [14]. This study is, to the best of our knowledge, the first application of diffusion tensor imaging and tractography in the human fetal heart, and provided valuable new insights into the transition of the heart from an isotropic tissue into a highly anisotropic organ. Moreover, while major technical hurdles remain, the performance of tractography of the human fetal heart in-utero is theoretically possible [20]. This could have profound implications for the diagnosis, treatment and prevention of congenital heart disease.


We thank Howard Chen for his assistance with histology.

Funding Statement

This work was funded in part by the following grants from the National Institutes of Health: R01 HL093038 and P41RR14075. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. No additional external funding received for this study.


1. Streeter DD Jr, Spotnitz HM, Patel DP, Ross J Jr, Sonnenblick EH (1969) Fiber orientation in the canine left ventricle during diastole and systole. Circ Res 24: 339–347 [PubMed]
2. LeGrice IJ, Smaill BH, Chai LZ, Edgar SG, Gavin JB, et al. (1995) Laminar structure of the heart: ventricular myocyte arrangement and connective tissue architecture in the dog. Am J Physiol 269: H571–582 [PubMed]
3. Chen J, Song SK, Liu W, McLean M, Allen JS, et al. (2003) Remodeling of cardiac fiber structure after infarction in rats quantified with diffusion tensor MRI. Am J Physiol Heart Circ Physiol 285: H946–954 [PubMed]
4. Helm PA, Younes L, Beg MF, Ennis DB, Leclercq C, et al. (2006) Evidence of structural remodeling in the dyssynchronous failing heart. Circ Res 98: 125–132 [PubMed]
5. Wu EX, Wu Y, Nicholls JM, Wang J, Liao S, et al. (2007) MR diffusion tensor imaging study of postinfarct myocardium structural remodeling in a porcine model. Magn Reson Med 58: 687–695 [PubMed]
6. Wu MT, Su MY, Huang YL, Chiou KR, Yang P, et al. (2009) Sequential changes of myocardial microstructure in patients postmyocardial infarction by diffusion-tensor cardiac MR: correlation with left ventricular structure and function. Circ Cardiovasc Imaging 2: 32–40 [PubMed]
7. Scollan DF, Holmes A, Winslow R, Forder J (1998) Histological validation of myocardial microstructure obtained from diffusion tensor magnetic resonance imaging. Am J Physiol 275: H2308–2318 [PubMed]
8. Smerup M, Nielsen E, Agger P, Frandsen J, Vestergaard-Poulsen P, et al. (2009) The three-dimensional arrangement of the myocytes aggregated together within the mammalian ventricular myocardium. Anat Rec (Hoboken) 292: 1–11 [PubMed]
9. Sosnovik DE, Wang R, Dai G, Wang T, Aikawa E, et al. (2009) Diffusion spectrum MRI tractography reveals the presence of a complex network of residual myofibers in infarcted myocardium. Circ Cardiovasc Imaging 2: 206–212 [PMC free article] [PubMed]
10. Mekkaoui C, Huang S, Chen HH, Dai G, Reese TG, et al. (2012) Fiber architecture in remodeled myocardium revealed with a quantitative diffusion CMR tractography framework and histological validation. J Cardiovasc Magn Reson 14: 70. [PMC free article] [PubMed]
11. Basser PJ (1995) Inferring microstructural features and the physiological state of tissues from diffusion-weighted images. NMR Biomed 8: 333–344 [PubMed]
12. Pierpaoli C, Basser PJ (1996) Toward a quantitative assessment of diffusion anisotropy. Magn Reson Med 36: 893–906 [PubMed]
13. Henderson DJ, Anderson RH (2009) The development and structure of the ventricles in the human heart. Pediatr Cardiol 30: 588–596 [PubMed]
14. Jouk PS, Usson Y, Michalowicz G, Grossi L (2000) Three-dimensional cartography of the pattern of the myofibres in the second trimester fetal human heart. Anat Embryol (Berl) 202: 103–118 [PubMed]
15. Hirschy A, Schatzmann F, Ehler E, Perriard JC (2006) Establishment of cardiac cytoarchitecture in the developing mouse heart. Dev Biol 289: 430–441 [PubMed]
16. Zhang L, Allen J, Hu L, Caruthers SD, Wickline SA, et al. (2013) Cardiomyocyte architectural plasticity in fetal, neonatal, and adult pig hearts delineated with diffusion tensor MRI. Am J Physiol Heart Circ Physiol 304: H246–252 [PubMed]
17. Phillips HM, Rhee HJ, Murdoch JN, Hildreth V, Peat JD, et al. (2007) Disruption of planar cell polarity signaling results in congenital heart defects and cardiomyopathy attributable to early cardiomyocyte disorganization. Circ Res 101: 137–145 [PubMed]
18. Henderson DJ, Chaudhry B (2011) Getting to the heart of planar cell polarity signaling. Birth Defects Res A Clin Mol Teratol 91: 460–467 [PubMed]
19. Tobita K, Garrison JB, Liu LJ, Tinney JP, Keller BB (2005) Three-dimensional myofiber architecture of the embryonic left ventricle during normal development and altered mechanical loads. Anat Rec A Discov Mol Cell Evol Biol 283: 193–201 [PubMed]
20. Roy CW, Seed M, van Amerom JF, Al Nafisi B, Grosse-Wortmann L, et al. .. (2013) Dynamic imaging of the fetal heart using metric optimized gating. Magn Reson Med. 10.1002/mrm.24614. [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science