Skip to main content

Ex vivo cardiovascular magnetic resonance diffusion weighted imaging in congenital heart disease, an insight into the microstructures of tetralogy of Fallot, biventricular and univentricular systemic right ventricle



Common types of congenital heart disease exhibit a variety of structural and functional variations which may be accompanied by changes in the myocardial microstructure. We aimed to compare myocardial architecture from magnetic resonance diffusion tensor imaging (DTI) in preserved pathology specimens.

Materials and methods

Pathology specimens (n = 24) formalin-fixed for 40.8 ± 7.9 years comprised tetralogy of Fallot (TOF, n = 10), dextro-transposition of great arteries (D-TGA, n = 8) five with ventricular septal defect (VSD), systemic right ventricle (n = 4), situs inversus totalis (SIT, n = 1) and levo-TGA (L-TGA, n = 1). Specimens were imaged using a custom spin-echo sequence and segmented automatically according to tissue volume fraction. In each specimen T1, T2, fractional anisotropy, mean diffusivity, helix angle (HA) and sheet angle (E2A) were quantified. Pathologies were compared according to their HA gradient, HA asymmetry and E2A mean value in each myocardial segment (anterior, posterior, septal and lateral walls).


TOF and D-TGA with VSD had decreased helix angle gradient by − 0.34°/% and remained symmetric in the septum in comparison to D-TGA without VSD. Helix angle range was decreased by 45°. It was associated with a decreased HA gradient in the right ventricular (RV) wall, i.e. predominant circumferential myocytes. The sheet angle in the septum of TOF was opposing those of the left ventricular (LV) free wall. Univentricular systemic RV had the lowest HA gradient (− 0.43°/%) and the highest HA asymmetry (75%). HA in SIT was linear, asymmetric, and reversed with a sign change at about 70% of the depth at mid-ventricle. In L-TGA with VSD, HA was asymmetric (90%) and its gradients were decreased in the septum, anterior and lateral wall.


The organization of the myocytes as determined by DTI differs between TOF, D-TGA, L-TGA, systemic RV and SIT specimens. These differences in cardiac structure may further enlighten our understanding of cardiac function in these diverse congenital heart diseases.


Congenital heart disease (CHD) is the most common malformation arising during fetal development, with increasing prevalence worldwide [1]. Advances in paediatric cardiac surgery and intensive care medicine have increased survival in the younger population and there are now more adults with CHD than children [2]. These adults have a high risk of myocardial dysfunction and heart failure. Many surgical interventions are more palliative than curative, with many patients requiring multiple surgeries through their lifetime [3, 4]. Although reparative or palliative surgery will provide adequate anatomical correction, it is likely that abnormalities persist within the myocardial structure [5]. Moreover, surgeries on newborns may alter the structural growth pattern, resulting in intrinsic tissue problems and regional remodeling [6, 7], leading to differences in blood flow patterns through the heart and early heart failure.

Diffusion tensor imaging (DTI) provides quantitative information on tissue microstructure. In myocardium, myocytes are packed with collagen into sheetlets to facilitate wall thickening [8, 9]. DTI exploits the anisotripic diffusion of water in the myocardium, with the degree of signal attenuation related to the diffusion along the direction of diffusion encoding gradients [10] [11]. Determining the diffusion-related signal attenuation over a number of different diffusion encoding gradient directions allows a diffusion tensor to be constructed. The eigenvalues of the diffusion tensor characterize the mean diffusivity and the fractional anisotropy [10], with the first eigenvector aligned with the main direction of the myocytes [12,13,14] and the second eigenvector with the sheetlet orientation. The angle of the myocytes relative to the cardiac short axis plane is known as helix angle (HA); the angle known as E2A defines the sheetlet angle in relation to the radial direction.

Although ex vivo and in vivo myocardial structure have been investigated with DTI in a number of studies [15,16,17,18,19], little work has been done to date in assessing the myocardial architecture of CHD hearts. We hypothesise that the typical helical architecture is altered in CHD with a variety of pathologies. Many tertiary institutions have legacy pathology specimens used for teaching and research, which comprise a valuable resource for non-destructive imaging studies. The information obtained from CHD pathological specimens can provide insight into the characteristics of CHD in particular patient groups. We developed custom tools for imaging historical pathology specimens, many of which have short T2, and automatic evaluation of the images. We used these methods to examine regional differences in common and rare types of CHD.

Materials and methods


DTI was performed on 24 historical formalin-fixed heart specimens, consisting of tetralogy of Fallot (TOF, n = 10), dextro transposition of the great arteries (D-TGA, n = 8), systemic right ventricle (RV) pathologies (n = 4, 2 single ventricle and two TGA with hypoplastic LV), situs inversus totalis (SIT, n = 1) and levo-TGA (L-TGA, n = 1). Of these 24 specimens, 17 had a ventricular septal defect (VSD), including all the TOF, 5 D-TGA specimens and 2 systemic RVs.

Image acquisition

One of the main challenges in scanning historical ex vivo specimens is the low signal-to-noise ratio (SNR) due to the reduction of T2 from the formalin fixation process. In this study, the specimens were scanned on a 3 T (Skyra, Siemens Healthineers, Erlangen Germany) with a custom spin-echo sequence with monopolar diffusion encoding gradients, eliminating susceptibility artefacts and minimising echo time (TE) as compared to typical echo planar imaging (EPI)-based techniques. To preserve the integrity of the tissues, all specimens were scanned in the same formalin-filled container that was used for their long-term storage. To comply with our ethical permissions, each specimen was scanned in its jar filled with formalin and placed at isocentre and scanned with a body 18 3 T Tim coil and the 32 channel spine coil. The diffusion sequence design, the data analysis for each specimen and congenital heart diseases groups are provided on

Imaging parameters used were: TR = 1 s, TE = 56.88 ms, δ = 20.67ms,  = 27.09 ms, slices = 3 (at mid-ventricle, or under the VSD when occurring), slice thickness = 4 mm, matrix = 100 × 100, field of view = 200 × 200 mm2, acquisition voxel = 2x2x4 mm3, 2 averages; 32 diffusion-encoding gradients on a single half-shell calculated by dmritool [20]. Total acquisition time for each specimen was 1 h. Non-diffusion weighted images (b = 0 s/mm2) were acquired every ten directions to minimize the load on the gradients. The optimal highest b value required without loss of diffusion information was calculated by b = 1.11/ADC [21, 22]. Since the mean diffusivity across specimens was 1.30 ± 0.27 × 10−3mm2/s, a b-value of b = 853 s/mm2 was used.

Definition of the radial vector

The radial vector was defined as the direction of minimal distance between epicardium and endocardium in both the left ventricle (LV) and RV. The location of each pixel between the endocardium and epicardium along the radial vector was calculated and expressed as a percentage, with 0% corresponding to the endocardium and 100% to the epicardium. The circumferential direction was defined to be perpendicular to the radial direction and the longitudinal direction was defined in the apex-base direction (out of the plane of the image).

Helix and sheet angles

Figure 1 details the calculation of the HA and E2A, respectively, according to the eigenvectors of the diffusion tensor with \( \overrightarrow{E_1} \) the main eigenvector which aligns with the orientation of the myocytes [9, 12, 23, 24]. The HA [°] was defined as the angle between the projection of the first eigenvector of the diffusion tensor onto the circumferential-longitudinal plane, and the circumferential vector. The E2A [°] was defined by the angle between the projection of the second eigenvector onto the radial-longitudinal plane and the radial vector. E2A was averaged over whole segment.

Fig. 1
figure 1

Helix angle (HA in red) and sheet angle (E2A in yellow) calculated from the projection of the eigenvectors of the diffusion tensor onto the circumferential and radial planes, respectively

The myocardium was manually segmented with papillary muscles removed from the myocardial mass. The myocardium was divided into five regions with four LV regions (anterior wall, lateral wall, posterior wall, and septum) and one region representing the RV free wall (Fig. 2).

Fig. 2
figure 2

A specimen with Levo transposition of the great arteries (L-TGA). a spin echo diffusion weighted image at b = 0 s/mm2; b myocardial segment (Anterior, Septum, Posterior, Lateral walls and RV wall); c helix angle; d sheet angle. LV, left ventricle, RV, right ventricle.

Helix angle gradient and asymmetry

HA values were averaged for each percentage of depth over the three acquired short axis slices acquired beneath the VSD or at the base of the ventricles in the absence of VSD. We observed that most segments showed a linear variation of HA with transmural depth over the range 20–80% depth (Fig. 3). An average HA line was fitted for each group and each segment using all points to calculate the linear gradient of the HA per degree of depth [°/%]. Endocardial HA were predominantly positive (right-handed) and epicardial angles were negative (left handed), leading to a negative gradient. Differences in slope between groups and between segments were tested using a t-test statistic incorporating error variances from each slope and weighing each of them by their degrees of freedom [25].

Fig. 3
figure 3

a Fitted line showing helix angle (HA) gradient calculated between 20 and 80% depth in biventricular dextro-tranposition of the great arteries (D-TGA) with and without ventricular septal defect (VSD). b Boxplot for D-TGA without VSD (n = 3). c Boxplot for D-TGA with VSD (n = 5). Middle line corresponds to the mean HA while outer dash lines are the 95% confident interval of the mean HA. a (b, n = 3) and with (c, n = 5). The Boxplot represents all HA values at each percentage of depth with the HA median as the red line, the first quartiles are the boxes and the dashed lines are the whiskers. The maximum whisker length is the interquartile range. HA values beyond the whiskers are displayed using “+”

Oblique myocytes may not necessarily realign circumferentially at mid-wall, influencing the amount of right or left handed myocytes across the wall. Asymmetry can be computed from the intersection of the fitted line and the x axis, i.e. the depth at which the HA is zero. Hence, HA0 will be less than 50% depth for predominantly left handed myocytes (predominant epicardial fibres) and HA0 will be greater than 50% depth for predominantly right handed myocytes (predominant endocardial fibres).


The collection of specimens

The median age at death was 8 years (quartile = [3:16], range = [1:46]) (Table 1). The median time spent in formalin was 41 years (quartile = [36:42.5], range = [29:57]). The mean T2 value over all specimens was 31.9 ± 6.5 ms and mean T1 was 146 ± 44 ms. The fractional anisotropy in the specimens was 0.24 ± 0.01.

Table 1 Information on the CHD specimens category, time in formalin, age of the individual at death, mean and standard deviation of T1 and T2, female(F) or male(M), the presence of VSD (Y:yes), biventricular (bi.) or univentricular (uni.) systemic right ventricle (RV), pulmonary valves stenosis (PVS) or atresia (PVA) is also specified. D-TGA = dextro-transposition of the great arteries; L-TGA – levo-transposition of the great arteries; TOF = tetralogy of Fallot;VSD = ventricular septal defect

Effect of a ventricular septal defect in D-TGA

Figure 3 shows how the helix angle changes with depth through the myocardium in biventricular D-TGA specimens with and without a VSD. The HA gradient was lower when a VSD was present compared to specimens without VSD in the septum and the RV free wall. The HA range decreased from [+ 55°:-42°] without VSD to [+ 28°:-33°] with VSD. Mean E2A in the septum changed from + 17° ± 18° in non VSD to − 20° ± 49° in VSD.

Tetralogy of Fallot and Dextro-transposition of the great arteries

Figure 4 shows the differences in HA gradient between the TOF group and each of the D-TGA groups (with and without VSD). In the septum, HA gradient was lower in TOF relative to D-TGA without VSD, and higher than D-TGA with VSD. In the lateral wall, HA gradient was higher in TOF than both D-TGA with VSD and D-TGA without VSD. In the RV free wall and LV anterior, the TOF group was similar to the D-TGA with VSD group but significantly lower than the D-TGA without VSD group (p < 0.007).

Fig. 4
figure 4

Helix angle gradient in tetralogy of Fallot (TOF), dextro-transposition of the great arteries (D-TGA) with and without ventricular septal defect (VSD). Significant differences between HA gradient (p < 0.01) are highlighted in green

The HA asymmetry in the septum of the TOF group was different from that in a bi-ventricular D-TGA with VSD (p < 0.001) but not from D-TGA without VSD (Fig. 5). HA asymmetry was different in RV wall of TOF and specimens with D-TGA (with or without VSD). Regardless of the pathology the HA asymmetry in the lateral wall was significantly different (p < 0.001).

Fig. 5
figure 5

HA asymmetry (intercept of HA) with the 95% confidence interval ([lower limit: upper limit]) in specimen of tetralogy of Fallot and biventricular transposition of the great arteries with or without ventricular septal defect (VSD). The p value measures the significant difference between two HA intercepts from two congenital heart disease groups in each segment

The mean E2A were opposing in the septum of TOF and biventricular D-TGA with VSD (p < 0.001) although not in biventricular D-TGA without VSD (Fig. 6). When a VSD was present, the E2A in the septum were opposite from those of the LV free wall (p < 0.001). D-TGA with and without VSD had opposing E2A in the LV. RV E2A were not different across CHD pathologies.

Fig. 6
figure 6

Sheet angle (E2A) in tetralogy of Fallot (TOF, n = 10), dextro-transposition of the great arteries (D-TGA) with and without ventricular septal defect (VSD). Significant differences (p < 0.01) are highlighted in green

Systemic right ventricle

Among all CHD groups, specimens presenting univentricular systemic RV had the most decreased HA gradient (Table 2) and the highest HA asymmetry (Table 3). These asymmetries characterised dominant right-handed angles in the single ventricle. In contrast, looking at the biventricular systemic RV group, HA gradient was steeper and HA symmetry was less left handed than the single ventricles.

Table 2 Helix angle (HA) gradient according to the myocardial regions. Specimens number #13 and #17 make the group of univentricular systemic RV while specimens #15 and #18 make the group of biventricular systemic RV. LV = left ventricle; RV = right ventricle
Table 3 Helix angle (HA) asymmetry calculated from 20 to 80% depth of the myocardial wall

The group with biventricular systemic RV, in comparison to the univentricular group, had similar sheet angles with the exception of the lateral wall (p = 0.008) (Table 4). Within this biventricular group, the sheet angle seen in the RV wall was similar to those in the LV, while the sheet angle in the septum could be differentiated from the free wall (p < 0.001). The sheet angle in biventricular systemic RV was positive, which was significantly different from the negative sheet angle seen in the TOF group (p < 0.001).

Table 4 Mean sheet angle (E2A) and standard deviation calculated from 20 to 80% depth of the myocardial wall

Situs inversus totalis

The HA gradient in the single SIT specimen was reversed (Table 5) relative to those seen in all other cases. However, HA was not linear in the region 20–80% and described a transition zone at about 70% of the depth at mid-ventricle which varied longitudinally from apex to base and according to the segment. Therefore, HA gradient was calculated within 20 to 70% depth. The HA asymmetries were segment dependent (Table 5). A scan along the longitudinal axis showed the transition from reversed helix (left-handed orientation) at the base to right-handed helix at the apex. The complete transition to normal HA appeared at about 1.2 cm from the apex and finished before segment 17.

Table 5 Helix angle (HA) gradient and asymmetry, as well as mean sheet angle (E2A), calculated from 20 to 70% depth of the myocardial wall in Situs Inversus Totalis

Levo-transposition of the great arteries

L-TGA with VSD had a more decreased HA gradient in the septum and anterior walls compared with the D-TGA with VSD group (Fig. 2, Table 6). In these segments, HA asymmetry was over 87% while the posterior and lateral wall remained with symmetric HA. HA asymmetry showed a rightward shift relative to D-TGA with VSD in the LV but not in the RV. The mean E2A could be statistically differentiated between dilated LV and hypertrophied RV and between the septum and the LV free wall (p < 0.001). E2A in the septum were more similar to D-TGA with VSD than D-TGA without VSD.

Table 6 Helix angle (HA) gradient and asymmetry, as well as mean sheet angle (E2A), calculated from 20 to 80% depth of the myocardial wall in L-TGA


Common CHD subtypes such as TOF and TGA are associated with increased prevalence of heart failure and adverse events [26,27,28,29], which may be partly associated with changes in myoarchitecture. Pathology specimens represent a valuable resource for studying the characteristics of CHD lesions. We found marked differences in cardiac microstructure associated with various CHD lesions in pathology specimens. In particular, we found that HA gradient was significantly reduced in the presence of a VSD, whether in TOF or D-TGA.

Ventricular septal defect in TOF and D-TGA

Sanchez-Quintana [30] observed in VSD of TOF from a visual inspection that outer angles of the septum in the RV and LV are more longitudinal until the VSD where they bend circumferentially to anchor on the sides of the VSD. The measured decreased HA in our results is consistent with this qualitative description.

In healthy hearts, the HA of the septum is sigmoidal (high HA gradient) because of the dominant longitudinal myocytes across the wall [31]. From biomechanical models, an increased HA gradient is associated with more longitudinal fibres contributing to longitudinal shortening and sigmoidal HA overcomes the linear HA configurations in all measured metrics (stroke volume, ejection fraction, base and apex thickening, shortening ratio, degree of torsion) [25]. The lack of connection of longitudinal layers in the septum was also observed in pathological human fetal heart [32].

The right ventricle in TOF and D-TGA

The presence of a VSD was also associated with a decreased HA gradient in the RV wall of TOF and D-TGA with VSD. The decreased HA gradient in the RV was not observed in the D-TGA specimens without VSD. The abnormal predominance of circumferential myocytes at mid ventricle in the RV may facilitate apical and basal dilatation, acting like a girdle on a balloon [33].

In the healthy heart, the epicardial myocytes of the RV wall wrap the LV free wall at the epicardium [32, 34] to insert in the endocardium by a crossover [35, 36]. This crossover contributes to the change of asymmetry observed in the anterior and lateral LV segment. In contrast, the endocardial myocytes of the RV merge into the septum to form the opposing oriented myocytes with the LV endocardial myocytes [31]. A change in the RV myocardial structure, such as seen in TOF and dTGA with VSD, might therefore be responsible of the changes observed in the epicardial structure of the LV at the septum, anterior and lateral walls.

Sheet angle in TOF and D-TGA

The E2A is independent from the HA organisation [37]. In the septum of TOF, the E2A was significantly different from the rest of the LV wall (p < 0.001). The septum bulges toward the LV [33, 38] which may explain a change in the sheet angle. We observed that septal sheets are opposing adjacent sheets in the anterior and posterior wall. The opposite sign of the E2A between TOF and D-TGA may come from the difference of E2A in the septal bulging between TOF and D-TGA without VSD caused by the difference of ventricular pressure [39]. The difference of the E2A between myocardial segments may result from the hemodynamic change due to case-specific anatomic, vascular and post-surgical characteristics.

Systemic right ventricle

The RV, regardless of the CHD pathology, had HA asymmetry beyond 70% (dominant endocardial orientation) with dominant circumferential myocytes across the myocardial wall (low helix angle gradient). The dominant endocardial orientation facilitates the contraction of the inner chamber at the expense of reduced epicardial lever arm force made of oblique and longitudinal orientations. Ventricular twist in RV is thus decreased and the torsion is lower, possibly increasing endocardial stress [40]. Consequently, a univentricular systemic RV has a structural disadvantage over univentricular systemic LV. The difference of HA asymmetry between univentricular and biventricular systemic RV would suggest that the absence of LV myocytes translates into the absence of LV shared parietal myocytes at the epicardium and its crossover toward the endocardium. Preserving the non-systemic ventricle during CHD surgery may maintain the orientation of the opposing myocytes necessary for efficient ventricular torsion.

In this study, univentricular systemic RV had the lowest HA gradient among all the investigated CHD groups, which is consistent with previous observations [41]. The biventricular systemic RV has previously been investigated with cardiac diffusion weighting in one patient [42]. The entire LV had an helicoidal loss with a predominance of circumferential myocytes (HA between 0 ° and 20 °) [42]. These preliminary findings are also in agreement with our measurements.


The changes observed in myocardial microstructure of these CHD specimens may affect their electrophysiology. Electrical activation encounters less resistivity along the myocytes’ length [43,44,45]. There is a strong interdependence between epicardial fibre direction, conduction velocity, resistivity of the myocardium, and the potential field surrounding and generated by a wave of depolarization [45]. As a result, any changes in the helical fibre angle can affect its electrical propagation properties in charge of the cardiac function. The helicoidal fibers minimize diffusion bias of the electrical wave propagation [46]. The changes we measured in myocardial microstructure in TOF may therefore be arrhythmogenic and explain the commonly observed nonuniform and delayed polarization that is known to be associated with malignant ventricular arrhythmia and sudden death in these patients [47, 48].

Situs Inversus Totalis and Levo transposition of the great arteries

The reversed HA in this case of SIT was in agreement with previous findings [49,50,51,52,53] and simulation by finite element modelling of the transition zone [54,55,56,57]. Our specimen shows the existence of the HA transition zone previously suggested in simulation to explain the reversed HA (positive HA gradient) at the base and normal HA at the apex (negative HA gradient). The transition pattern influences the sign and amplitude of systolic torsion - from positive to negative torsion [55]. The longitudinal location of these transitions can affect the systolic function [58] which may explain the inter-individual differences of torsion and gradient of circumferential radial shear in SIT. In fact, the transmural location of these HA transition pattern, specified as HA asymmetry, were observed in the helical arrangement of the basal part and vary across SIT patients [50,51,52], thus directly impacting the torsion pattern in the apical part.

Although a L-TGA has ventricles and valves reversed, the myocytes orientations are not reversed [52]. As previously observed in other VSD specimens, sheetlets in the septum were significantly different from the LV free wall. HA asymmetry of the LV can be a factor influencing RV output as it determines the percentage of oriented epicardial myocytes throughout the wall. The epicardial myocytes were more circumferential than longitudinal which might affect the resultant RV shortening by the LV.

The analysis of single case specimen either in L-TGA or SIT limits the conclusion about the myocytes’ orientation for this CHD group. The difference may either be due to specific defects from the individual rather than from the group or by way of noise from a single data point measurement.


The specimens investigated in this study comprised numerous specific and individualised diagnoses (Table S1). Further, it is possible that the histological assessment did not reference all the present pathology which may bias the interpretation of the sub-categories of defects. As a result, specimens were only categorised in broad terms according to the type of CHD, and further specific analysis of sub-categories was not considered.

The collection of pathological specimens spans many decades of changing treatment and outcome. In many of cases, the conditions were untreated, or the treatment was unsuccessful. For example, some of the cases with transposition had an atrial switch procedure but did not survive long after operation. Others succumbed from RV failure and others from pulmonary vascular obstructive disease and left heart failure. Thus, we might observe a combination of abnormal fiber orientations associated with the congenital condition and then modified by the evolving pathology. Further, the results of congenital cardiac surgery have improved significantly improved during the last five decades [3, 59]. As a result, those cases in a single category (D-TGA for example) could represent a heterogeneous group with different surgical repair techniques affecting the microstructure.

Owing to the nature of the tissue collection used in this study, preserving the integrity of the specimens was paramount. Other studies assessing myocardial fiber structure have used specific tissue preparation techniques prior to scanning, or performed histological analysis on tissues [17]. In our study, manipulating the specimens or performing any destructive procedures was not possible. As such we were not able to validate findings with histology.

The phase of the cardiac cycle affects the HA and E2A [17]. Given that we average the HA or the E2A over several specimens of the same category we expect to be more robust by getting an average estimate of the E2A regardless of the cardiac phase at fixation. This average and their respective standard deviation represent the micro-structure features specific to the CHD rather than to the cardiac phase cycle. However, with smaller number of specimens we become more subject to this cardiac fixation stage and the specific defect or noise measurements from the specimen(s).

The CHD collection investigated in our study has a disparate heart maturation which can affect the interpretation of fibre orientations as a consequence of aging rather than the pathology itself. Sanchez-Quintana [38] found that neonate hearts have circumferential epicardial helix while few weeks before birth they become oblique. Before 15 years old, the endocardial angles in the RV and LV run longitudinally and become slightly oblique to arch around the tricuspid valves. Specimens without VSD or individuals at the age of 15 years old may present more longitudinal endocardial angles than others (high HA gradient).


Myocardial architecture can be non-destructively examined in specimens of CHD with CMR DTI. The differences in microstructure observed between different types of CHD specimens may explain some of the anatomical and functional observations made in CHD patients.

Availability of data and materials

the datasets supporting the conclusions of this article are available in the github repository, It contains the design of the diffusion sequences, the data analysis for each specimen and congenital heart diseases groups.



Apparent diffusion coefficient


Congenital heart disease


Diffusion tensor imaging


Dextro-transposition of the great arteries


Sheet angle (second eigenvector angle)


Echo planar imaging


Fractional anisotropy


Helix angle


Levo-transposition of the great arteries


Left ventricle/left ventricular


Right ventricle/right ventricular


Situs inversus totalis


Signal-to-noise ratio


Transverse relaxation time


Spin-lattice relaxation time


Echo time


Transposition of the great arteries


Repetition time


Tetralogy of Fallot


Ventricular septal defect


  1. Marelli AJ, Ionescu-Ittu R, Mackie AS, et al. Lifetime prevalence of congenital heart disease in the general population from 2000 to 2010. Circulation. 2014;130:749–56.

    Article  PubMed  Google Scholar 

  2. Ou P, Iserin L, Raisky O, et al. Post-operative cardiac lesions after cardiac surgery in childhood. Pediatr Radiol. 2010;40:885–94.

    Article  PubMed  Google Scholar 

  3. Popelová JR, Gebauer R, Černý Š, et al. Operations of adults with congenital heart disease – single center experience with 10 years results. Cor Vasa. 2016;58:e317–27.

    Article  Google Scholar 

  4. Putman LM, van Gameren M, Meijboom FJ, et al. Seventeen years of adult congenital heart surgery: a single Centre experience. Eur J Cardiothorac Surg. 2009;36:96–104; discussion 104.

    Article  PubMed  Google Scholar 

  5. Bolger AP, Coats AJS, Gatzoulis MA. Congenital heart disease: the original heart failure syndrome. Eur Heart J. 2003;24:970–6.

    Article  PubMed  Google Scholar 

  6. Hein S, Arnon E, Kostin S, et al. Progression from compensated hypertrophy to failure in the pressure-overloaded human: heart structural deterioration and compensatory mechanisms. Circulation. 2003;107:984–91.

    Article  PubMed  Google Scholar 

  7. Reynolds HR, Tunick PA, Grossi EA, et al. Paradoxical septal motion after cardiac surgery: a review of 3,292 cases. Clin Cardiol. 2007;30:621–3.

    Article  PubMed  PubMed Central  Google Scholar 

  8. LeGrice IJ, Smaill BH, Chai LZ, et al. Laminar structure of the heart: ventricular myocyte arrangement and connective tissue architecture in the dog. Am J Phys. 1995;269:H571–82.

    CAS  Google Scholar 

  9. Scollan DF, Holmes A, Zhang J, Winslow RL. Reconstruction of cardiac ventricular geometry and Fiber orientation using magnetic resonance imaging. Ann Biomed Eng. 2000;28:934–44.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. Basser PJ, Mattiello J, LeBihan D. MR diffusion tensor spectroscopy and imaging. Biophys J. 1994;66:259–67.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. Le Bihan D, Mangin JF, Poupon C, et al. Diffusion tensor imaging: concepts and applications. J Magn Reson Imaging. 2001;13:534–46.

    Article  PubMed  Google Scholar 

  12. Reese TG, Weisskoff RM, Smith RN, et al. Imaging myocardial fiber architecture in vivo with magnetic resonance. Magn Reson Med. 1995;34:786–91.

    Article  PubMed  CAS  Google Scholar 

  13. Scollan DF, Holmes A, Winslow R, Forder J. Histological validation of myocardial microstructure obtained from diffusion tensor magnetic resonance imaging. Am J Phys. 1998;275:H2308–18.

    CAS  Google Scholar 

  14. Hsu EW, Muzikant AL, Matulevicius SA, et al. Magnetic resonance myocardial fiber-orientation mapping with direct histological correlation. Am J Phys. 1998;274:H1627–34.

    CAS  Google Scholar 

  15. Lombaert H. Statistical Atlas of Human Cardiac Fibers Comparison with Abnormal Hearts; 2012.

    Book  Google Scholar 

  16. Lombaert H, Peyrat JM, Fanton L, Cheriet F (2012) Variability of the Human Cardiac Laminar Structure. Accessed 26 May 2016.

  17. Nielles-Vallespin S, Khalique Z, Ferreira PF, et al. Assessment of myocardial microstructural dynamics by in vivo diffusion tensor cardiac magnetic resonance. J Am Coll Cardiol. 2017;69:661–76.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Ferreira P, Kilner PJ, McGill L-A, et al. Aberrant myocardial sheetlet mobility in hypertrophic cardiomyopathy detected using in vivo cardiovascular magnetic resonance diffusion tensor imaging. J Cardiovasc Magn Reson. 2014;16:P338.

    Article  PubMed Central  Google Scholar 

  19. Khalique Z, Ferreira PF, Scott AD, et al. Diffusion tensor cardiovascular magnetic resonance of microstructural recovery in dilated cardiomyopathy. JACC Cardiovasc Imaging. 2018.

  20. Cheng J, Shen D, Yap P-T, Basser PJ. Novel single and multiple Shell uniform sampling schemes for diffusion MRI using spherical codes. Cham: Springer; 2015. p. 28–36.

    Google Scholar 

  21. Conturo TE, McKinstry RC, Aronovitz JA, Neil JJ. Diffusion MRI: precision, accuracy and flow effects. NMR Biomed. 1995;8:307–32.

    Article  PubMed  CAS  Google Scholar 

  22. Jones DK, Horsfield MA, Simmons A. Optimal strategies for measuring diffusion in anisotropic systems by magnetic resonance imaging. Magn Reson Med. 1999;42:515–25.<515::AID-MRM14>3.0.CO;2-Q.

    Article  PubMed  CAS  Google Scholar 

  23. Garrido L, Wedeen VJ, Kwong KK, et al. Anisotropy of water diffusion in the myocardium of the rat. Circ Res. 1994;74:789–93.

    Article  PubMed  CAS  Google Scholar 

  24. Holmes AA, Scollan DF, Winslow RL. Direct histological validation of diffusion tensor MRI in formaldehyde- fixed myocardium. Magn Reson Med. 2000;44:157–61.<157::AID-MRM22>3.0.CO;2-F.

    Article  PubMed  CAS  Google Scholar 

  25. Colombo NM. Numerical modelling of ventricular mechanics : role of the myofibre architecture; 2014.

    Google Scholar 

  26. Zomer AC, Uiterwaal CSPM, Der Velde ET Van, et al (2011) Mortality in adult congenital heart disease: are national registries reliable for cause of death? Int J Cardiol 152:212–217.

  27. Engelings CC, Helm PC, Abdul-Khaliq H, et al. Cause of death in adults with congenital heart disease - an analysis of the German National Register for congenital heart defects. Int J Cardiol. 2016;211:31–6.

    Article  PubMed  Google Scholar 

  28. Brickner ME, Hillis LD, Lange RA. Congenital heart disease in adults. N Engl J Med. 2000;342:256–63.

    Article  PubMed  CAS  Google Scholar 

  29. Hoffman JI, Kaplan S, Liberthson RR. Prevalence of congenital heart disease. Am Heart J. 2004;147:425–39.

    Article  PubMed  Google Scholar 

  30. Sanchez-Quintana. Ventricular myoarchitecture in tetralogy of Fallot. Heart. 1996;76:280–6.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. Varray F, Mirea I, Langer M, et al. Extraction of the 3D local orientation of myocytes in human cardiac tissue using X-ray phase-contrast micro-tomography and multi-scale analysis. Med Image Anal. 2017;38:117–32.

    Article  PubMed  Google Scholar 

  32. Jouk P-S, Usson Y, Michalowicz G, Grossi L. Three-dimensional cartography of the pattern of the myofibres in the second trimester fetal human heart. Anat Embryol (Berl). 2000;202:103–18.

    Article  CAS  Google Scholar 

  33. Sheehan FH, Ge S, Vick GW III, et al. Three-dimensional shape analysis of right ventricular remodeling in repaired tetralogy of Fallot. Am J Cardiol. 2008;101:107–13.

    Article  PubMed  Google Scholar 

  34. Smerup M, Nielsen E, Agger P, et al. The three-dimensional arrangement of the myocytes aggregated together within the mammalian ventricular myocardium. Anat Rec. 2009;292:1–11.

    Article  Google Scholar 

  35. Bovendeerd PHM, Huyghe JM, Arts T, et al. Influence of endocardial-epicardial crossover of muscle fibers on left ventricular wall mechanics. J Biomech. 1994;27:941–51.

    Article  PubMed  CAS  Google Scholar 

  36. Greenbaum RA, Ho SY, Gibson DG, et al. Left ventricular fibre architecture in man. Heart. 1981;45:248–63.

    Article  CAS  Google Scholar 

  37. Helm PA, Younes L, Beg MF, et al. Evidence of structural remodeling in the dyssynchronous failing heart. Circ Res. 2006;98:125–32.

    Article  PubMed  CAS  Google Scholar 

  38. Sanchez-Quintana D, Garcia-Martinez V, Climent V, et al. Morphological changes in the normal pattern of ventricular myoarchitecture in the developing human heart. Anat Rec. 1995;243:483–95.

    Article  PubMed  CAS  Google Scholar 

  39. Trzebiatowska-Krzynska A, Swahn E, Wallby L, et al. Afterload dependence of right ventricular myocardial deformation: a comparison between tetralogy of Fallot and atrially corrected transposition of the great arteries in adult patients. PLoS One. 2018;13:e0204435.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Young AA, Cowan BR. Evaluation of left ventricular torsion by cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2012;14:49.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Corno AF, Kocica MJ. Potential implications of the helical heart in congenital heart defects. Semin Thorac Cardiovasc Surg Pediatr Card Surg Annu. 2007;10:61–7.

    Article  Google Scholar 

  42. ​Harmer J, Toussaint N, Pushparajah K, Stoeck CT, Chan RW, Razavi R, Kozerke S. In-vivo Diffusion Tensor Imaging of the Systemic Right Ventricle at 3T. Salt Lake City. Proc Intl Soc Mag Reson Med. 2013;21(3098).

  43. Chang M-C, Wu M-T, Weng K-P, et al. Left ventricular regional myocardial motion and twist function in repaired tetralogy of Fallot evaluated by magnetic resonance tissue phase mapping. Eur Radiol. 2018;28:104–14.

    Article  PubMed  Google Scholar 

  44. Colli Franzone P, Guerri L, Pennacchio M, Taccardi B. Spread of excitation in 3-D models of the anisotropic cardiac tissue. II. Effects of fiber architecture and ventricular geometry. Math Biosci. 1998;147:131–71.

    Article  Google Scholar 

  45. Roberts D, Hersh L, Scher A. Influence of cardiac Fiber orientation on Wavefront voltage, conduction velocity, and tissue resistivity in the dog. Circ Res. 1979;44:701–12.

    Article  PubMed  CAS  Google Scholar 

  46. Aumentado-Armstrong T, Kadivar A, Savadjiev P, et al. Conduction in the Heart Wall: Helicoidal fibers minimize diffusion Bias. Sci Rep. 2018;8:7165.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. Gatzoulis MA, Balaji S, Webber SA, et al. Risk factors for arrhythmia and sudden cardiac death late after repair of tetralogy of Fallot: a multicentre study. Lancet. 2000;356:975–81.

    Article  PubMed  CAS  Google Scholar 

  48. Uebing A, Gibson DG, Babu-Narayan SV, et al. Right ventricular mechanics and QRS duration in patients with repaired tetralogy of Fallot implications of Infundibular disease. Circulation. 2007;116:1532–9.

    Article  PubMed  Google Scholar 

  49. Baillie M (1789) An account of a remarkable transposition of the viscera in the human body. London Med J 10:178–197.

  50. Taussig HB. The anatomy of the heart in two cases of situs transversus. Bull Johns Hopkins Hosp. 1926;39:199–202. .

  51. Matsumura H, Aizawa Y, Kumaki K. Myocardial architecture in situs inversus viscerum totalis. In: Developmental Cardiology: Morphogenesis and Function, edited by Clark EB and Takao A. Mount Kisco: Futura; 1990. p. 605–24.

  52. Asami I, Koizumi K. The vortex cordis is never reversely directed, even in situs inversus and L-loop anomaly. Kaibogaku Zasshi. 1989;64:36–45.

    PubMed  CAS  Google Scholar 

  53. Khalique Z, Ferreira PF, Scott AD, et al. Deranged Myocyte microstructure in Situs Inversus Totalis demonstrated by diffusion tensor cardiac magnetic resonance. JACC Cardiovasc Imaging. 2018.

  54. Delhaas T, Decaluwe W, Rubbens M, et al. Cardiac fiber orientation and the left-right asymmetry determining mechanism. Ann N Y Acad Sci. 2004;1015:190–201.

    Article  PubMed  Google Scholar 

  55. Delhaas T, Kroon W, Decaluwe W, et al. Structure and torsion of the normal and situs inversus totalis cardiac left ventricle. I. Experimental data in humans. Am J Physiol Circ Physiol. 2008;295:H197–201.

    Article  CAS  Google Scholar 

  56. Kroon W, Delhaas T, Arts T, Bovendeerd P. Constitutive modeling of cardiac tissue growth. In: Functional imaging and modeling of the heart. Berlin Heidelberg, Berlin, Heidelberg: Springer; 2007. p. 340–9.

    Chapter  Google Scholar 

  57. Rossi AC, Pluijmert M, Bovendeerd PHM, et al. Assessment and comparison of left ventricular shear in normal and situs inversus totalis hearts by means of magnetic resonance tagging. Am J Physiol Circ Physiol. 2015;308:H416–23.

    Article  CAS  Google Scholar 

  58. Pluijmert M, Kroon W, Delhaas T, Bovendeerd PHM. Adaptive reorientation of cardiac myofibers: the long-term effect of initial and boundary conditions. Mech Res Commun. 2012;42:60–7.

    Article  Google Scholar 

  59. Vaujois L, Gorincour G, Alison M, et al. Imaging of postoperative tetralogy of Fallot repair. Diagn Interv Imaging. 2016;97:549–60.

    Article  PubMed  CAS  Google Scholar 

Download references


The authors would like to acknowledge the Auckland District Health Board, as well as the Heart Registry Governance Group for their part in facilitating the use of the heart specimens for this project, Rhonda Holloway for her valuable assistance in accessing and transporting the specimens, and Dame Naida Glavish for her guidance with obtaining ethical approval for this study.


This study was funded by the National Heart Foundation of New Zealand with support from The Lady Alport Barker Trust and The T M Hosking Charitable Trust.

Author information

Authors and Affiliations



C.T developed the diffusion sequence, analysed and interpreted the data. A.Y, T.G. and B.P were involved with the ethical approval and interpreted the data. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Beau P. Pontré.

Ethics declarations

Ethics approval and consent to participate

Ethical approval for the use of the specimens was obtained from the national health and disability ethics committee. The Auckland District Health Board Heart Registry provides access to human heart tissue for education, training and study purposes with strictly adhered to policies and procedures and under the guidance of a Heart Registry Custodian. The Registry is overseen by the Heart Registry Governance Group which consented to the use of tissue for this study. This article does not contain any studies with human participants or animals performed by any of the authors.

Consent for publication

Not applicable.

Competing interests

C.T., A.Y, T.G. and B.P declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Complete diagnosis of the specimens according to their corresponding scan order (CHD number #).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tous, C., Gentles, T.L., Young, A.A. et al. Ex vivo cardiovascular magnetic resonance diffusion weighted imaging in congenital heart disease, an insight into the microstructures of tetralogy of Fallot, biventricular and univentricular systemic right ventricle. J Cardiovasc Magn Reson 22, 69 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: