Changyu SUN1, Lihui WANG2, Feng YANG3, and Yuemin ZHU1
1CREATIS; CNRS UMR 5220; Inserm U 1044; INSA of Lyon; University of Lyon;, Lyon, France, 2School of computer science and technology, Guizhou University, Guiyang, China, People's Republic of, 3School of Computer and Information Technology, Beijing JiaoTong University, Beijing, China, People's Republic of
Synopsis
Fiber crossing in the myocardium of the heart is controversial despite some reports
on the heart of animals. The purpose of this study is to investigate the fiber
configuration complexity in a voxel, such as fiber crossing, in the human
heart. We applied two different reconstruction techniques commonly used in high
angular resolution diffusion imaging (HARDI) to conventional diffusion tensor
imaging (DTI) data of the human heart acquired on clinical MRI system. The
results consistently demonstrated that the human myocardium exhibits some complex
fiber configurations such as fiber crossing in a voxel.Introduction
Unlike
in the brain, fiber crossing in the myocardium of the heart is controversial
despite some reports on the heart of animals [1] [2], and even more so when there is no
consensus on the definition itself of fiber in the heart. Nevertheless, with a
spatial resolution of the order of millimeter of clinical MRI systems, we can
by the abuse of language use the term fiber to designate a group of myocytes
having a privileged orientation in the heart. The purpose of this study is to investigate the fiber configuration complexity in a voxel,
such as fiber crossing, in the human heart.
Materials and
Methods
The
real data concerns an ex vivo human heart. It was acquired using a conventional
diffusion tensor imaging (DTI) protocol on a Siemens 1.5T Magnetom Sonata with
the following settings: TE=73ms, TR=6400ms, FOV=256×256mm, slice thickness=2
mm, slice spacing=2 mm, number of slices=54, image size=128×128, diffusion
sensitivity b-value=1000s/mm
2, gradient directions=64 and a single
T2-weighted (b-value=0) image. The heart was located in a plastic container and
fixed by hydrophilic gel to maintain a diastolic shape. This setup has a low
dielectric effect and eliminates unwanted susceptibility artifacts near the
boundaries of the heart. To reconstruct the maps of fiber configurations, we employed
two different techniques in high angular resolution diffusion imaging (HARDI): the
constrained spherical deconvolution (CSD) [3] and the analytical q-ball imaging
(AQBI) [4][5], which are commonly used for brain, but until now never applied
to heart. Strictly speaking, the CSD generates the fiber orientation
distribution (FOD) map, and the AQBI the orientation distribution function
(ODF) map. For simplicity, we call both of them the ODF.
Results
The
ODF maps reconstructed using CSD and AQBI are shown in Fig. 1 (21
th
slice). The zones marked by the red rectangle contain the voxels with complex
fiber configurations. To show more clearly the complex fiber configurations revealed
by CSD and AQBI, in Fig. 2(a) to (c) are visualized their zoomed versions. The
ODFs obtained from both methods show more or less the existence of intravoxel
complex fiber configurations in the middle of the red rectangular region. The
fiber crossings are more pronounced by CSD, but with AQBI, we can still observe
that it does not concern simple single fiber configurations. As an illustration,
a clearly crossing fiber (reconstructed by CSD) in the voxel circled by the black
box is shown in Fig. 2(c). Other fiber crossing voxels were also found near
this voxel. By exploring the neighboring slices around the 21th
slice, at the same region, we also observed the existence of complex fiber
configurations, as shown in Fig. 3 (ODFs obtained using CSD).
Discussion
On
the same DTI datasets, two different reconstruction methods led to the same
observation that the myocardium of the human heart presents complex fiber
configurations. This result is not due to the partial volume since the voxel is
in the middle of the myocardium. It is not due to noise either since the fiber
crossing occurs also at exactly the same position in the neighboring slices
around the current slice. However, the two methods did not yield the same
degree of resolving the fiber crossing. It would then be necessary to validate
such findings either via simulation or via some physical measurements like
polarized light imaging.
Conclusion
This
work has investigated for the first time the fiber configuration complexity of
the human heart. The results consistently demonstrated that the human
myocardium exhibits clearly some complex fiber configurations in a voxel such
as fiber crossing when using both CSD and AQBI reconstruction techniques, which
is particularly interesting since the results are obtained on conventional DTI
data (so shorter acquisition time). Since the two different reconstruction
methods led to slightly different complex fiber configuration results, further
study both on data acquisition schemes such as HARDI and on the reconstruction
of fiber configurations would be interesting for confirming these preliminary
findings or searching fiber configuration information other than fiber
orientation.
Acknowledgements
The authors would like to thank P. Croisille, S.
Rapacchi and M. Viallon for their help in acquiring the DTI data. The work is
supported by French ANR MOSIFAH.References
[1] Sosnovik D et al, J.
Cardiovasc. Magn. Reson, 47: 1-15, 2009; [2] Shi Y et al, In Proc Intl Soc Mag Reson
Med, 1470, 2007; [3] Tournier JD et al,
Neuroimage, 42:
617-625, 2008; [4] Maxime D et al, Magn
Reson Med, 3:
497-510, 2007; [5] SUN CY et al, Phys.
Med. Biol., 60:8417-8436,2015.