- Research
- Open access
- Published:
Cardiac structure discontinuities revealed by ex-vivo microstructural characterization. A focus on the basal inferoseptal left ventricle region
Journal of Cardiovascular Magnetic Resonance volume 25, Article number: 78 (2023)
Abstract
Background
While the microstructure of the left ventricle (LV) has been largely described, only a few studies investigated the right ventricular insertion point (RVIP). It was accepted that the aggregate cardiomyocytes organization was much more complex due to the intersection of the ventricular cavities but a precise structural characterization in the human heart was lacking even if clinical phenotypes related to right ventricular wall stress or arrhythmia were observed in this region.
Methods
MRI-derived anatomical imaging (150 µm3) and diffusion tensor imaging (600 µm3) were performed in large mammalian whole hearts (human: N = 5, sheep: N = 5). Fractional anisotropy, aggregate cardiomyocytes orientations and tractography were compared within both species. Aggregate cardiomyocytes orientation on one ex-vivo sheep whole heart was then computed using structure tensor imaging (STI) from 21 µm isotropic acquisition acquired with micro computed tomography (MicroCT) imaging. Macroscopic and histological examination were performed. Lastly, experimental cardiomyocytes orientation distribution was then compared to the usual rule-based model using electrophysiological (EP) modeling. Electrical activity was modeled with the monodomain formulation.
Results
The RVIP at the level of the inferior ventricular septum presented a unique arrangement of aggregate cardiomyocytes. An abrupt, mid-myocardial change in cardiomyocytes orientation was observed, delimiting a triangle-shaped region, present in both sheep and human hearts. FA’s histogram distribution (mean ± std: 0.29 ± 0.06) of the identified region as well as the main dimension (22.2 mm ± 5.6 mm) was found homogeneous across samples and species. Averaged volume is 0.34 cm3 ± 0.15 cm3. Both local activation time (LAT) and morphology of pseudo-ECGs were strongly impacted with delayed LAT and change in peak-to-peak amplitude in the simulated wedge model.
Conclusion
The study was the first to describe the 3D cardiomyocytes architecture of the basal inferoseptal left ventricle region in human hearts and identify the presence of a well-organized aggregate cardiomyocytes arrangement and cardiac structural discontinuities. The results might offer a better appreciation of clinical phenotypes like RVIP-late gadolinium enhancement or uncommon idiopathic ventricular arrhythmias (VA) originating from this region.
Introduction
The growing number of imaging technologies that have emerged over the last decades has greatly extended our knowledge of cardiac microstructure, and given detailed information about 3D cardiomyocyte architecture [1]. Nevertheless, there is still a gap in our understanding of both the course of events of pathophysiological processes like cardiac remodeling [2], and the link with clinical phenotypes [3]. Recently, some concerns have been raised regarding the interpretation of the presence of late gadolinium enhancement (LGE) at the right ventricle attachment or insertion point (RVIP) [4,5,6,7,8,9].
The origin of the LGE pattern has not been extensively described in the literature. The existing studies [10,11,12] examined it using histopathology of human postmortem hearts in the context of hypertrophic cardiomyopathy (HCM) or pulmonary arterial hypertension (PAH). More recently, Friedberg et al. [13] investigated the role of regional myocardial remodeling in rat or rabbit in the context of PAH using echocardiography and histopathology. All these studies reported the presence of myocardial disarray, expanded extracellular space, disorderly arranged myocytes, or an increase of regional fibrosis. The only published 3D descriptions were two studies that quantified either myoarchitectural disarray in HCM rats using high-resolution episcopic microscopy [14] or infarcted myocardium using tractrography processing in patients [15]. In contrast to the above-mentioned studies, we recently reported an abrupt mid-myocardial change in cardiomyocytes orientation at the RVIP in ex-vivo large mammalian samples [16] and did not report the presence of myoarchitectural disarray but the study was limited to a few sheep samples.
The purpose of the present work is to characterize in ex-vivo large mammalian (human and sheep) hearts, the basal inferoseptal (BIS) cardiomyocyte orientation and assess its functional impact on electrophysiology. In particular, it aims to detail the 3D structure organization by identifying the cardiomyocyte organization with high-resolution diffusion MRI at 600 µm3, and X-ray imaging at 21 µm3 and 2D histology. Experimental characterization of the microstructure organization is also largely used in realistic or patient-specific models for biomechanical, or electrophysiological models. Recently, Doste et al. [17] underlined the discrepancy between cardiomyocyte orientation distributions computed by rule-based models with histological or MRI data in the IVS due to the dual-layer organization [18], and in the RVIP. Thus, to expand the scope of the result, based on the observed human experimental data, we model the effects of cardiomyocyte orientation on cardiac electrophysiology (EP), computing local activation times and pseudo-ECGs. Lastly, we discuss the hypothetical structural determinism relating aggregate cardiomyocyte orientation with the presence of focal LGE in this area.
Methods
Sample preparation
All hearts were fixed with a solution containing formalin (10%) and a gadolinium-based contrast agent (Dotarem, Guerbet, Paris, France) at 0.2% of total volume of perfusion, by retrograde perfusion from the aorta. Sample preparation was performed as previously described [16, 18,19,20].
Human samples
Human samples (N = 5) were derived from the human donor program (providing access to heart samples from patients under cerebral death for scientific research purposes) approved by the ‘‘Agence Française de la Biomedecine’’ and with a written informed consent of the patient’s family. The experiment was conducted in accordance with the declaration of Helsinki and the institutional ethics committee. Five hearts (1 male and 4 females, 71 ± 15.5 years old) were used in the study. A brief summary of patient characteristics is available in Table 1. Hearts were reported with cardiac pathology and not eligible for cardiac transplantation. Donors of hearts #1 and #2 had a history of prior myocardial infarction, donor #3 had a history of mitral valve prolapse and donor #4 had a history of hypertensive cardiomyopathy and atrial fibrillation. For one heart (#5), no significant cardiac pathology was revealed during the clinical evaluation performed before explanting the heart [20]. This heart was thus labeled “control”.
Sheep samples for MRI and MicroCT
Sheep samples (N = 6; 5 for MRI, 1 for MicroCT) and protocols used in this study were approved by the Animal Research Ethics Committee (CEEA—050 Comité d’éthique pour l’expérimentation animale Bordeaux) in accordance with the European rules for animal experimentation. Hearts were explanted via sternal thoracotomy under general anesthesia using the following protocol. In brief, anesthesia was induced with an intravenous bolus of ketamine (15 mg/kg) and midazolam (1.5 mg/kg). Animals were then intubated and ventilated and received an injection of heparin (2.5 mg/kg). Anesthesia was maintained with ketamine and midazolam (40 mg kg−1 h−1 and 2 mg kg−1 h−1, respectively). The thorax was opened, and blood from each animal was collected via the introduction of an 8-Fr sheath into the right jugular vein. Heparin (15 mg/l) was added in the reservoir to avoid coagulation. Sheep were euthanized by intraperitoneal injection with sodium pentobarbital (35 mg/kg) and the hearts were quickly excised. The aorta was cannulated and perfused with cold cardioplegic (4 °C) solution supplemented with heparin (5 U/ml). The hearts were mostly arrested in the systolic phase (#S1-S4) with the exception of sample #S5 which was arrested in a cardiac phase closer to diastole. The myocardium is also relaxed in comparison to in-vivo state due to this first cardioplegic flushing before formalin fixation macroscopic examination was performed to verify the absence of overt cardiac disease.
Additional protocol for microCT
A novel tissue preparation technique to alleviate background photon attenuation was performed using a tissue air-drying approach [21]. In this case, left and right coronary ostia were cannulated individually for optimal tissue perfusion. The #S6 heart (N = 1) was perfused with phosphate-buffered saline (PBS, pH = 7.4) followed by ethylenediaminetetraacetic acid (EDTA) [10 mM] solution with Papaverine [30 μM]. The heart was fixed 2 h with 4% formaldehyde, rinsed three times for 20 min in PBS and then dehydrated using ethanol and perfused with hexamethyldisilazane (HMDS) to prevent tissue deformation. Finally, the heart was hung to air-dry under the fume hood and inside small containers to reduce airflow for 7 days.
Data acquisition
MRI acquisition
All experiments were performed at 9.4T with an inner bore size of 30 cm (BioSpin MRI; Bruker, Ettlingen, Germany) equipped with 300 mT/m gradient insert of 200 mm inner diameter adapted for large samples. Images were acquired using a dedicated radiofrequency volume array coil with seven elements in transmit and receive. All MRI scans and experiment acquisitions were similar to those described previously in [16, 18,19,20].
Diffusion-weighted (DW) imaging
A 3D DW spin-echo sequence was used to acquire DW images (TE/TR: 22/500 ms, acquisition time: 24 h). The DW dataset consisted of six noncollinear diffusion encoding directions acquired with a b value of 1000 s/mm2 with an isotropic spatial resolution of 600 µm3. The same parameters were set for each acquisition except for the field of view which was adjusted for each sample.
A 2D DW spin-echo sequence was used to acquire DW images with six noncollinear diffusion encoding directions to test the influence of b-value. The slice was positioned in SA view in mid-ventricular region on sheep sample #S2. 5 acquisitions (TE/TR = 22/500 ms, 5 averages, acquisition time: 24′13″) were done at 0.6 × 0.6 × 1 mm3 voxel resolution with b-values ranging at 100, 350, 500, 750 and 1000 s/mm2.
Anatomical imaging
3D anatomical image was acquired with a gradient echo sequence (TE/TR: 9/30 ms, acquisition time: 18 h) at an isotropic resolution of 150 µm3. The resulting volume is referred to as the anatomical image.
MicroCT acquisitions
Heart #S6 was imaged using microCT (SkyScan 1276, Bruker, Belgium). X-ray transmission images were acquired using X-ray source energies of 55 kV and 150 μA, and an aluminum plate of thickness 0.5 mm. An average of three acquisitions were captured at rotation steps incrementing by 0.18° over 360°. The total acquisition time was 8 h. Short-axis (SA) images were obtained by 3D tomographic reconstruction of the raw images using the cone-beam FDK11 algorithm using the NRecon software (Bruker). Images were reconstructed with an isotropic voxel resolution of 21 μm. The matrix size was 1801 × 1344 × 3872 for a corresponding field of view of 38 × 28 × 81 mm3.
Image analysis
MRI diffusion tensor estimation
Each DW image was up-sampled by a factor of 2 using trilinear interpolation to reach a voxel size of 0.3 × 0.3 × 0.3 mm3. Diffusion tensor calculations were performed prior to any registration to avoid difficult transformations of the diffusion encoding matrix. Diffusion tensor maps [eigenvalues: λ1, λ2, λ3, apparent diffusion coefficient, fractional anisotropy (FA), and color-coded FA (cFA), also known as red–green–blue colormap] were obtained with MRtrix3 software (https://www.mrtrix.org) [22]. The first diffusion tensor eigenvector v1 was associated with the main aggregate cardiomyocytes orientation. To obtain a similar alignment among hearts and ease the visualization of the data, the long axis (LA) of the LV were aligned to the z-axis of each volume using rigid transformation. The LV was subdivided into regions using a polar coordinate system and the helix angles (HA) were calculated for each sample.
Retrospective sampling to the standard in-vivo resolution
The diffusion tensor images were retrospectively downsampled to two standard in-vivo resolutions of 2.0 × 2.0 × 5.0 mm3 and 2.5 × 2.5 × 8 mm3. The downsampling was performed in the log-Euclidean space. Diffusion tensor maps were computed as described in the previous section in the previous paragraph.
Estimation of the whole heart label
A radio-frequency (RF) inhomogeneity field correction was applied to non-diffusion weighted images (b0) with ITK-N4 [23] to reduce the intensity bias and signal drop-off in the apical and basal regions associated with the RF coil. A binary mask was created and applied on both anatomical and DW images to segment myocardial tissue from regions with hyperintense signal corresponding to residual formaldehyde. Low and high cutoff thresholds were applied on the FA, trace, and DW images to define the binary mask as described previously in [19].
Tractography processing
Tractography was used for visualization of global cardiomyocytes arrangement. Streamlines were generated with MRtrix3 [22] using the principal eigenvector of the diffusion tensor. Unless specified, the tractography parameters were defined as an FA-stopping threshold of 0.1, a maximum angle between steps of 45, a step size of 0.05 mm, a maximum length of 20 mm, and a minimum length value of 1 mm. Thus, streamlines terminated if a streamline made a sharp turn (angles larger than 45°), crossed a voxel with an FA of less than 0.1, extended to 20 mm, or extended outside of the masked region.
Delineation of the aggregate cardiomyocytes bundle
Tractography processing (summarized in Fig. 1) was adapted to delimit and further quantify the size of aggregate cardiomyocytes bundles located at the RVIP, also referred to as the singularity.
i) Tractograms generation with preferential orientation.
As the aggregate cardiomyocytes bundle was visualized with a preferential orientation to the LA, tractograms were generated with a preferential orientation in the LA direction to isolate the aggregate cardiomyocytes bundle from the surrounding cardiomyocytes having a circumferential orientation (Fig. 1A). 16 tractrograms were generated by varying FA-stopping threshold and maximum angle as described in Fig. 1B. One set of parameters was picked arbitrarily (FA-stopping threshold of 0.2 and a maximum angle of 30°). The same values were used for all samples.
ii) Conversion of tractogram into structure delineation.
Track Density Imaging (TDI) was used to interpolate streamline information in a Cartesian frame [24]. Thus, for each element of a grid, the total number of streamlines was counted and converted into a scalar volume image. Then, intensity-based thresholding was applied on the scalar volume image to obtain a mask, delineating the aggregate cardiomyocytes in the heart with a preferential orientation to LA as shown in Fig. 1C. Lastly, refinement of the mask was manually performed using 3DSlicer (https://www.slicer.org/) to remove voxels outside the region of interest.
iii) Characterization of the aggregate cardiomyocytes bundles.
The volumes of the aggregate cardiomyocytes bundle and myocardium of the whole heart (without fat) were calculated by counting the corresponding number of voxels. The distribution of FA in the aggregate cardiomyocytes bundle and the whole heart were compared.
To quantify the size of the aggregate cardiomyocytes bundle compared to the whole heart, different distances have been used. These distances were measured manually using the software 3DSlicer and are described in. Figure 1D. The first (dBA) was the distance between the lowest point of the heart (Apex) and the highest point between the mitral valve, the same height of the top of the septum (Base). dCR and dCL were two distances between the center or the singularity and the closest point of the RV and LV respectively. The wall thickness (dWT) was calculated with the two points of the LV and RV chosen for dCR and dCL. The dSLA was the distance of the singularity on the LA view.
STI estimation and tractography on microCT
3D STI was used on microCT images to estimate the aggregate cardiomyocytes orientation similarly to [14, 25]. First, the outer product of the intensity gradient vectors was computed to estimate the derivatives. The derivatives were then convolved with a Gaussian kernel with standard deviation σ = 6 and assembled to generate a symmetric second-moment matrix or structure tensor per voxel. The principal directions of the structure tensor in each voxel were extracted using eigen analysis. The third eigenvector (smallest magnitude eigenvalue) corresponded to cardiomyocytes orientation and was used to generate the streamlines. The tractography parameters were adapted to the signal-to-noise ratio and imaging resolution of the microCT images with an FA-stopping threshold of 0.1, a maximum angle between steps of 15, a step size of 0.05 mm, a maximum length of 20 mm, and a minimum length value of 0.5 mm.
Histology
The sheep heart #1 was dissected to assess the presence of cardiomyocytes in the region depicting a triangular shape. The sample was transected from base to apex in large slices of 1 cm thick. The anterior part and the RV wall were removed to isolate the posterior wall of the heart. Sample was first rinsed in phosphate-buffered saline, then agitated at 4 °C for 24 h. After samples have been dehydrated, it was embedded in paraffin and sectioned at 6 µm in the transmural direction (or SA view). Tissue sections were stained with Masson’s trichrome for structural identification: collagen fibers in tissues were green; nuclei were black; and myocytes were red/pink. Slices were examined at 10× magnification on a tissue slide scanner (Axio Scan Z1, Carl Zeiss SAS, Marly le Roi, France). Images were visualized using Zen lite software (ZEN Blue 2.6, Carl Zeiss Microscopy, Thornwood, New York, USA).
Electrophysiology simulation
Details about the simulations of the electrical activation propagation are provided in Fig. 2 and in Methods in the Additional file 1.
Results
Visualization of the aggregate cardiomyocytes bundle in the BIS wall in sheep hearts
Figure 3 displays a whole sheep heart #4 (A) and the location (blue cube) of the BIS region using 3D rendering (B, C). Anatomical data shows a relatively homogeneous tissue (D, E). The FA maps computed from the diffusion tensor images (F) show a more heterogeneous distribution. A triangular pattern presenting lower FA values on the edge is somehow visible (yellow arrows). The triangular pattern has a clear image signature on cFA or tractography map (indicated by pink arrows) with the aggregate cardiomyocytes oriented in base to apex direction starting at the inferobasal crux and ending at the middle of the inferoseptal wall of the LV. The surrounding aggregate cardiomyocytes have, in the endocardial part of the LV, a standard circumferential orientation. On the opposite side, at the RVIP, a corridor of aggregate cardiomyocytes, indicated by a purple arrow in the SA view, goes from the septum to connect with the tissues of the RV wall.
Additional file 1: Fig. S1 plots the profile of the FA map with the corresponding SA view at the location of the identified pattern. The mean value of FA along the line is: 0.20 ± 0.08 and a decreased value is noticeable in the edge of the pattern indicating potential aggregate cardiomyocytes discontinuity or crossing in this region.
Figure 4 investigates the reproducibility across different samples of the same species in five sheep hearts. In anatomical images, the intensity is relatively homogeneous except for a few areas with reduced signal in the epicardial fat. Nevertheless, small variation in gray contrast in samples #S1–S3 are indicated by black arrows. #S1–S4 (fixed in a pseudo-systolic phase) show a large wall thickness, the last sample #S5 (fixed in a pseudo-diastolic phase) presents a small wall thickness. On hearts #S1–S4, the same decrease in the FA value is visible (indicated by the yellow arrow) on the edge of the triangular pattern, while it is not visible on sample #S5. Smooth transition of aggregate cardiomyocytes orientation is visible (gradient of green to purple to red) indicating a radial (in-plane) aggregate cardiomyocytes orientation in the endocardial part of the LV while aggregate cardiomyocytes orientation changes abruptly between adjacent voxels in the interventricular septum (#S1–S3, #S5) and in the middle of the RVIP for all samples with a blue color that depicts an apex-base orientation (#S1–S5). The contrast of the color, slightly different in sample #S5, indicates a difference in orientation, but similar discontinuities in aggregate cardiomyocytes orientation remain. The LA view is available in Additional file 1: Fig. S2 that displays the extension of the aggregate cardiomyocytes bundles that vanished close to the mid-ventricular location. Additional file 1: Fig. S3 displays both SA and LA view with HA maps.
Figure 5 displays the cardiomyocytes architecture using microCT and STI. The left panel displays the whole sheep heart #S6 (A) using 3D rendering and the location of the BIS area using a blue cube (B). On the right panel, the anatomical contrast offers a direct visualization of the laminar structure. In the SA view, circumferential cardiomyocytes are visible (green arrow) while cardiomyocytes crossing the view are located in the middle of the inferior RVIP (yellow arrows). The longitudinal views also show the presence of cardiomyocytes bundles with specific base to apex orientation (yellow arrows). Visually, the cardiomyocytes’ density appears more pronounced in comparison to the surrounding cleavage planes showing alternating cardiomyocytes and air layers. The reconstructed cFA maps and the streamlines display a triangular pattern with cardiomyocytes bundles in base to apex orientation (white arrows).
Macroscopic examination and histological examination were performed on sheep heart #S1. Additional file 1: Fig. S4 shows a photograph of the sample after short-axis transection, confirming the macroscopic size of the identified region. Visual agreement with the anatomical MRI image is evident. Figure 6 investigates the presence and orientation of cardiomyocytes in the region. The histology slice was first compared to MRI-derived data using landmarks (yellow asterisks) at the vessel locations (Panels 6A–F). Tissue staining revealed the presence of cardiomyocytes in the entire triangular regions and a clear image signature about myocyte orientation was found. Inside the region (panel 6G), myocytes are in transverse orientation depicting a base-apex direction. The surrounding myocytes outside the region (panel 6H and 6I), have a standard longitudinal orientation. Lastly, visualization of abrupt mid-myocardial change in myocytes orientation is clearly visible (Panels 6J–L) at the junction confirming the cardiac structure discontinuities revealed by DTI measurements.
Visualization of the aggregate cardiomyocytes bundles in the BIS area in human hearts
Figure 7 investigates the presence and the reproducibility across different human samples. Samples #H1–2 presented myocardial infarction scars. Samples #H3 and #H4 came from patients with minimal cardiac disease and #H5 came from a healthy, albeit aged, subject. Gross morphology as well as myocardial viability was highly heterogeneous. A large infarct scar is visible on sample #H2 with reduced wall thickness on the inferior and inferolateral walls. The presence of epicardial fat is noticeable in all samples. Anatomical data shows a relatively homogeneous tissue with the exception of sample #H2. On the FA map, the edge of the triangular pattern is well visible for #H1, #H2, and #H4 samples (yellow arrows) in the SA view, while it is not clearly visible on sample #H3 and #H5. The cFA maps and the streamline rendering display the typical pattern in each sample with aggregate cardiomyocytes bundles in base to apex orientation indicated by the white arrow. Although it is up to the reader's appreciation, the geometry remains triangular for samples #H1–4 but with a larger heterogeneity. The sample #H5 presents relatively small aggregate cardiomyocytes bundles at the junction and the pattern does not look triangular. The LA view is available in Fig. 8, and shows the singularity dives from base to apex at a depth of ~ 2 cm. Additional file 1: Fig. S5 displays both SA and LA view with helix angle maps.
Characterization of the fractional anisotropy in the singularity
Figure 9 shows the distribution of FA in the singularity (using roughly 15,000 voxels) and the whole heart (using roughly 5.106 voxels) for all samples. In average, there is a small shift between the FA distribution in the whole heart and in the singularity for both sheep (except #S4) and human samples (except #H3). The average mean and standard deviation values (summarized in Additional file 1: Table S1) in the whole heart are (mean ± std: #S: 0.26 ± 0.09, #H: 0.24 ± 0.10) and in the singularity (mean ± std: #S 0.29 ± 0.06, #H: 0.27 ± 0.05) for sheep and human samples respectively. The histogram distribution of the FA in the singularity is relatively homogeneous across samples and species, and a narrow distribution of FA was found in comparison to the whole heart distribution.
Sizing of the singularity and potential of in vivo imaging
Additional file 1: Fig. S6 and Table S2 summarize the size and volume estimation of the aggregate cardiomyocytes bundles after structure delineation using TDI. Base-Apex distances (dBA) of the ventricle ranged from 70 to 85 mm for both species. Wall thickness (dWT) was 15 mm on average and showed higher variability for human hearts. The long axis of the singularity (dSLA) was estimated to 23.0 ± 5.6 mm (both for human and sheep hearts) with the exception of the sample #H5. The singularity is closer to the RV cavity than the LV cavity as the wall thickness of the LV is higher. In all cases, the singularity is mostly inside the LV assuming that the RV/LV limit is a B-slice interpolation of the LV epicardium. The volume of the singularity ranges from 0.1 to 0.6 cm3 with a mean value of 0.35 ± 0.15 cm3. Additional file 1: Fig. S7 shows the cFA maps of all human samples downsampled to 2.0 × 2.0 × 5 mm3 and 2.5 × 2.5 × 8 mm3 and display a few remaining voxels of the singularity. Additional file 1: Fig. S8 shows the RVIP of #S2 with different b-values. The first row shows the average of DWI volumes and the second shows the cFA as a function of b-values. The cardiomyocyte orientation starts to be visible from b = 350 s/mm2 and the aggregate of myocytes in base apex orientation seems more visible from b = 500 s/mm2 and display a similar pattern like in Fig. 4.
Cardiac electrophysiology simulation
Spatial distribution of local activation times is shown in Fig. 10 using the usual rule-based "fiber" orientation (1), experimentally measuremented “fibers” (2), and with a non-conducting interface at the region border (3). For all activation patterns, incorporating experimental “fibers” slows wavefront propagation, leading to significant delays ranging from 10 to 40 ms. Such differences are reflected in LAT maps. The average LAT differences between cases (1) and (2) are 21.4 ± 8.0 ms, 19.9 ± 6.6 ms and 17.2 ± 4.8 ms for the top to bottom, left to right and front to back propagations, respectively. LAT differences between cases (1) and (3) are 22.0 ± 9.8 ms, 21.3 ± 9.7 ms and 19.0 ± 9.0 ms as before.
“Fiber” orientation influence is reflected in the pseudo-ECGs (Additional file 1: Fig. S4). Only electrode pairs perpendicular to the propagation wavefronts are shown. Incorporating experimental “fibers” changed the morphology and amplitudes of all pseudo-ECGs. For top-to-bottom stimulation, the highest amplitude of the B1-A1 voltage changes from 4.96 mV to 4.8 mV to 5.68 mV, for rule-based, experimental, and experimental + interface, respectively. Analogously, for front-to-back activation, T1-T2 voltage amplitude changed from 7.69 to 4.05 and 3.45 mV, for rule-based, experimental, and experimental + interface, respectively. Finally, for left-to-right activation, the E1-E2 electrode voltage amplitude changed from 1.97 mV to 4.85 and 4.58 mV, as previously.
The non-conducting interface is more apparent for bipolar points crossing it. The amplitude discrepancies between test cases (2) and (3) are more pronounced for top to bottom activation and right to left activation.
The RVIP bundle interface effect is clearly seen in electrograms traversing it. Current flow is affected by two factors: (1) the interface is non-conducting with holes which slow down propagation, since the bundle is excited at discrete points, rather than a plane wave, and (2) fiber direction within the RVIP bundle is distinct from adjacent myocardium, producing different conduction velocities within the bundle relative to the electrodes. Also, experimentally measured fibers result in faster trans-slab conduction to change the time of the second spike in the activation complex, corresponding to far surface excitation. This peak was the same for experimental fibers regardless of interface, but different for rule-based fibers. The interface is responsible for the pseudo-ECG morphology. When electrical coupling across an interface depolarizes transmembrane voltage, but insufficiently to initiate excitation, ST elevation results [26]. When electrode pairs do not cross the interface (F1-F2 & E2-E3), pseudo-ECGs are almost equal.
Discussion
In this study, the microstructure organization of two large mammalian species was investigated in the BIS myocardium using three distinct imaging methods. First, we highlight the existence of abrupt mid-myocardial cardiomyocyte orientation changes using DTI at an isotropic spatial resolution of 600 µm, delimiting a triangle-shaped region, in both sheep and human samples. The discontinuity can be identified either using the FA maps that quantify the anisotropy of water molecule motion in a voxel, or using directional information of the diffusion tensor via cFA maps or streamlines representation. FA distribution was found systematically equal or higher in both species in the identified aggregate cardiomyocytes bundles than in the whole hearts. The mean FA found in the whole hearts was in agreement with previous ex-vivo studies [27]. The proposed description was then confirmed by microCT acquisition that offers direct visualization of the laminar structure in this region.
In sheep hearts, the triangular pattern shows comparable configurations across samples corresponding to the systolic phase. The sample fixed in the diastolic phase reports an overall shift in aggregate cardiomyocytes orientation but the discontinuity remains clearly visible. In human hearts, a more heterogeneous pattern is noticeable in particular for the sample #H2 which has an important scar along the free wall (Fig. 7) but is visible in all samples. We quantified the volume and size of the aggregate cardiomyocyte bundles using state of art tractographic processing. The segmentation of the singularity was performed using TDI, which allows the delineation of specific regions by converting streamline statistics in scalar images [24]. The mean length of dSLA and volume of the aggregate cardiomyocytes were 23.7 ± 4.1 mm and 0.42 ± 0.10 cm3 for sheep, and 20.8 ± 7.0 mm and 0.26 ± 0.15 cm3 for human hearts, respectively. To assess the presence of cardiomyocytes in the region depicting a triangular shape, macroscopic examination and histological examination were performed on sheep heart #S1. The histological results (Fig. 6) confirm the presence of cardiomyocytes in the region and demonstrated an excellent match between DTI and histological technique to follow the main cardiomyocytes orientation. Lastly, the existence of cardiac structure discontinuities (or abrupt change in cardiomyocytes orientation) revealed by DTI measurements was confirmed at the junction of the identified region. Our study is proof of existence, but generalization of our observations is difficult regarding the low number of samples (N = 11) and several questions remain open. The findings are at odds with literature that reports either highly organized myocyte architecture or disorganized myocyte architecture at the RVIP. Indeed, a recent meta-analysis on ventricular cardiomyocytes [28] refutes the presence of internal tissue structures that could be interpreted as boundaries.
The physiological hypothesis
The presence of these aggregate cardiomyocytes bundles in healthy and relatively young animals could suggest that this aggregate cardiomyocyte bundles organization is physiological. A naive and speculative hypothesis attributes the presence of this structure for stiffening the region located at the intersection of the ventricular and atrial chambers. A similar role has been reported for the papillary muscles [29]. While different configurations and sizes are reported in the study, the inclusion of ex-vivo human samples could not confirm this hypothesis as most donors had a history of heart disease. In such a context, the human sample #H2 that presents a large scar in the inferior and inferolateral wall is an interesting case. This heart had undergone a significant reduction of the myocardial wall thickness in the LV, and a massive replacement of the working myocardium by fibrosis, visible on the anatomical images (Fig. 6). However, we can distinguish on the FA map, an isolated region with a higher value of FA where the aggregate cardiomyocyte bundles are present in base to apex direction. This could be interpreted that the region is stiffer and less sensitive to remodeling. Additionally, rare examples are also present in the literature such as Fig. 2 from [30], that showed a similar pattern in a fetal heart after 20 weeks of gestation or in canine heart [31]. Prior work by Mekkaoui et al. [15] used tractographic techniques for myocardial structural characterization and showed the presence of discontinuity at the insertion point. The propagation angle, a tractography-based metric quantifying the curvature of streamlines, computed in this study was consistently found higher at the RVIP. While under debate in the literature, we can also underline that gross pathology description (limited to the content of 3.2 paragraph: “The intervetricular sulci”) by Kocica et al. [32] closely match our current and previous [16] experimental findings.
The cardiac remodeling hypothesis
The cardiomyocyte orientation at the RVIP follows the rule-based helix angle orientation in the literature, either in microstructural or modeling studies [33]. Nevertheless, it is often accepted that the cardiomyocyte organization is much more complex due to the intersection of the ventricular cavities. Stretch-induced adaptive reorientation of cardiomyocytes has been reported in the literature [34] and is supposed to play a large role in the pathological process. Thus, cardiac remodeling could be considered sufficient to substantially alter the orientation in the LV of circumferential cardiomyocytes toward base to apex direction. The myocardium at the junction could be subjected to a complex mechanical pressure that might lead to heterogeneous shear force or stress. Nevertheless, such a hypothesis would imply an asymmetrical and localized response to the mechanical pressure. First, no evidence of similar aggregate cardiomyocyte bundles was found in the superior wall. Secondly, aggregate cardiomyocyte orientation changes rapidly and not gradually across the inferior RVIP. Computational mechanical models could be envisioned to explore different scenarios of remodeling but are out of the scope of the paper.
Fibrosis implications
RVIP-LGE was reported in multiple patient populations including hypertrophic cardiomyopathy [6, 35], dilated cardiomyopathy (DCM) [5], non-ischemic cardiomyopathy (NICM) [7], in heart failure patients with preserved ejection fraction [8], in highly trained endurance athletes [9, 36], and in all patients with pulmonary arterial hypertension (PAH) [3, 13]. The clinical relevance of RVIP-LGE and prognostic differ in the literature [2, 37, 38]. The LGE pattern was found to be a nonspecific clinical marker in patients without cardiac damage [4] but also a marker of PH for more advanced disease [37]. The described aggregate cardiomyocyte bundles probably play a role in the aforementioned clinical studies. The presence of focal LGE could be a chronic mechanical stress marker.
Two mechanisms could be considered within and at the border of the aggregate cardiomyocytes bundles: i) in this region, the forces related to ventricular contraction could be expressed differently at the cell level of these two aggregate cardiomyocytes populations (circumferentially and apico-basal oriented) by producing a longitudinal and perpendicular stretch respectively which may result in cell remodeling, and ii) the presence of an abrupt change in orientation might also be a potential substrate for interstitial fibrosis. Experimentally, a decrease of FA was noticeable at the interface of the aggregate cardiomyocyte bundles. This could indicate either cross orientation of cardiomyocytes within a voxel, or also the presence of interstitial fibrosis [26].
Electrophysiology implication
Adequate simulation of arrhythmic phenomena requires accurate patient-specific data. Although cardiomyocyte orientation has been shown to have major effects on wavefront propagation in electrophysiology modeling [39], most studies rely on rule-based information derived from ex-vivo canine samples [40] that do not fully reflect the heterogeneity of the human myoarchitecture. In this study a realistic cardiomyocyte orientation distribution was compared to the usual rule-based model in a block of tissue located in the BIS myocardium. Both spatial distribution of LAT and pseudo-ECGs were strongly impacted with delayed LATs and changes in amplitude and morphology.
The impact of fibrosis is usually driven via the in-silico integration of heterogeneous tissue composition at the cellular level [41]. Here we suggested the presence of cardiomyocyte orientation discontinuity as a potential substrate for fibrosis and leading to conduction slowing. To investigate such a scenario, a simulation was performed with the presence of a non-conducting interface at the border of the identified region. The morphology of pseudo-ECGs was changed, specifically by adding the interface, a notch was observed in the bipolar signals in each of the activations. Reentry or electrical wavefront discontinuity leading to arrhythmia was not observed but were unlikely due to the small dimension of the block of tissue and the length of the action potential duration. Nevertheless, uncommon idiopathic ventricular arrhythmias (VA) [42,43,44] originating from the basal inferoseptal left ventricular myocardium have been reported and could originate from the abovementioned aggregate cardiomyocyte bundles.
Limitations
The extension of these findings is, however, limited by a large number of factors with notably the small number of samples and the contractile state of the ex vivo samples. In vivo imaging of the identified structure could be envisioned using new DTI-CMR acquisition [45] but was out scope of the paper. To investigate such a possibility, the diffusion tensor of all human samples were downsampled up to 2.5 × 2.5 × 8 mm3. The new voxel size (50 mm3) is 231 times larger than the previous one (0.216 mm3) and increases the partial volume effect. Nevertheless, a few remaining voxels (< 10 per sample) displayed in Additional file 1: Fig. S7 still shows characteristic cardiomyocyte orientation in the RVIP in samples 1 to 4. However, the transition to in-vivo acquisition is extremely challenging due to the combined cardiac and respiratory motion. We also investigated the influence of b-values and found that a b-value b > = 500 s/mm2 might be necessary to visualize the identified structure. It is reasonable to underline that the two main factors to image the region will be a sufficient spatial resolution in plane and in slice thickness to avoid too large a partial volume and a drop in FA, and of course adequate signal-to-noise ratio. The usual limitations regarding the interpretation of tractography and streamlines apply due to the large difference between the myocyte size and maximum streamline length.
In-vivo cardiac DTI could serve as a biomarker of structural remodeling and may be used in the close future for risk stratification. Relationships between FA, fibrosis and ventricular arrhythmia have recently been investigated in this way [46]. Such clinical perspectives are encouraging for improved diagnosis of cardiac pathologies but also underline the need for a precise description of the baseline cardiomyocyte organization at multiple scales and imaging modalities. To some extent, the concept of myocardial disarray must be used with caution, as the more complex aggregate cardiomyocyte orientation in the inferior RVIP may not be considered as pathological myoarchitectural disarray but might be a substrate for pathological mechanisms.
The development of high-resolution accelerated diffusion MRI acquisition coupled with standardized sample preparation and post-processing pipeline could help in designing a modern atlas of the cardiac aggregate cardiomyocyte architecture [15, 47, 48]. To support our findings and promote the reproducibility of the results, the raw data and diffusion tensor images used in the article and minimal scripts to reproduce the figure of the article were released (see data availability section).
Conclusion
The study was the first to describe the 3D cardiomyocyte architecture of the BIS ventricular region in large mammalian hearts. A triangular region depicting a arrangement of cardiomyocytes in the basal–apical direction starting at the inferobasal crux and ending at the middle of the inferoseptal wall of the LV was found in both species using diffusion tensor imaging methods. The results were confirmed by microCT imaging and histology. EP computational models were then used to investigate whether the presence of cardiomyocytes orientation discontinuity, being a potential substrate for fibrosis, could promote a conduction slowing. Lastly, this peculiar arrangement may provide an explanation as to why this region may be damaged by increased right ventricular wall stress. It is hoped that these basic research findings could help in providing a better understanding of the biological processes underlying cardiac remodeling or arrhythmia vulnerability and will pave the way for linking microstructural observations to distinct clinical manifestations.
Availability of data and materials
Link to the data is indicated here (https://github.com/valeryozenne/Cardiac-Structure-Database/tree/master/Article-4). Results are visible with the mrview viewer using the panels overlay/tensor/tractography. Command scripts that compute in native space (diffusion tensor, first eigenvectors, tractography) and command-line instructions for generating figures are also available on the same webpage.
References
Sosnovik DE. Magnetic resonance-based characterization of myocardial architecture. Heart Fail Clin. 2021;17(1):85–101.
Andersen S, et al. Right ventricular fibrosis. Circulation. 2019;139(2):269–85.
Glasenapp A, et al. Molecular imaging of inflammation and fibrosis in pressure overload heart failure. Circ Res. 2021;129(3):369–82.
Grigoratos C, et al. Clinical importance of late gadolinium enhancement at right ventricular insertion points in otherwise normal hearts. Int J Cardiovasc Imaging. 2020;36(5):913–20.
Mikami Y, et al. Right ventricular insertion site fibrosis in a dilated cardiomyopathy referral population: phenotypic associations and value for the prediction of heart failure admission or death. J Cardiovasc Magn Reson. 2021;23(1):79.
Chan RH, et al. Prognostic value of quantitative contrast-enhanced cardiovascular magnetic resonance for the evaluation of sudden death risk in patients with hypertrophic cardiomyopathy. Circulation. 2014;130(6):484–95.
Yi JE, et al. Prognostic implications of late gadolinium enhancement at the right ventricular insertion point in patients with non-ischemic dilated cardiomyopathy: a multicenter retrospective cohort study. PLoS ONE. 2018;13(11): e0208100.
Nitsche C, et al. Native T1 time of right ventricular insertion points by cardiac magnetic resonance: relation with invasive haemodynamics and outcome in heart failure with preserved ejection fraction. Eur Heart J Cardiovasc Imaging. 2020;21(6):683–91.
Domenech-Ximenos B, et al. Prevalence and pattern of cardiovascular magnetic resonance late gadolinium enhancement in highly trained endurance athletes. J Cardiovasc Magn Reson. 2020;22(1):62.
Kuribayashi T, Roberts WC. Myocardial disarray at junction of ventricular septum and left and right ventricular free walls in hypertrophic cardiomyopathy. Am J Cardiol. 1992;70(15):1333–40.
Bulkley BH, D’Amico B, Taylor AL. Extensive myocardial fiber disarray in aortic and pulmonary atresia. Relevance to hypertrophic cardiomyopathy. Circulation. 1983;67(1):191–8.
Bradlow WM, Assomull R, Kilner PJ, Gibbs JSR, Sheppard MN, Mohiaddin RH. Understanding late gadolinium enhancement in pulmonary hypertension. Circ Cardiovasc Imaging. 2010;3(4):501–3.
Friedberg MK, Redington AN. Right versus left ventricular failure: differences, similarities, and interactions. Circulation. 2014;129(9):1033–44.
Garcia-Canadilla P, et al. Myoarchitectural disarray of hypertrophic cardiomyopathy begins pre-birth. J Anat. 2019;235(5):962–76.
Mekkaoui C, Jackowski MP, Kostis WJ, et al. Myocardial scar delineation using diffusion tensor magnetic resonance tractography. J Am Heart Assoc. 2018;7(3): e007834. https://doi.org/10.1161/JAHA.117.007834.
Magat J, et al. A groupwise registration and tractography framework for cardiac myofiber architecture description by diffusion MRI: an application to the ventricular junctions. PLoS ONE. 2022;17(7): e0271279.
Doste R, et al. A rule-based method to model myocardial fiber orientation in cardiac biventricular geometries with outflow tracts. Int J Numer Method Biomed Eng. 2019;35(4): e3185.
Rodriguez-Padilla J, et al. Impact of intraventricular septal fiber orientation on cardiac electromechanical function. Am J Physiol Heart Circ Physiol. 2022;322(6):H936–52.
Magat J, et al. 3D MRI of explanted sheep hearts with submillimeter isotropic spatial resolution: comparison between diffusion tensor and structure tensor imaging. MAGMA. 2021;34:741.
Haliot K, et al. A 3D high resolution MRI method for the visualization of cardiac fibro-fatty infiltrations. Sci Rep. 2021;11(1):9266.
Pallares-Lupon, N et al. Optimizing large organ scale micro computed tomography imaging in pig and human hearts using a novel air-drying technique. bioRxiv, 2021: p. 2021.07.29.454121.
Tournier JD, Calamante F, Connelly A. Robust determination of the fibre orientation distribution in diffusion MRI: non-negativity constrained super-resolved spherical deconvolution. Neuroimage. 2007;35(4):1459–72.
Tustison NJ, et al. N4ITK: improved N3 bias correction. IEEE Trans Med Imaging. 2010;29(6):1310–20.
Calamante F, et al. Track-density imaging (TDI): super-resolution white matter imaging using whole-brain track-density mapping. Neuroimage. 2010;53(4):1233–43.
Gilbert S et al. Measurement of myocardial structure: 3D structure tensor analysis of high resolution MRI quantitatively compared to DT-MRI. 2013. 7746: p. 207–214.
Vigmond EJ, Efimov IR, Rentschler SL, Coronel R, Boukens BJ. Fractionated electrograms with ST-segment elevation recorded from the human right ventricular outflow tract. HeartRhythm Case Rep. 2017;3(11):546–50.
Abdullah OM, et al. Characterization of diffuse fibrosis in the failing human heart via diffusion tensor imaging and quantitative histological validation. NMR Biomed. 2014;27(11):1378–86.
MacIver DH, et al. The end of the unique myocardial band: Part I. Anatomical considerations. Eur J Cardio Thorac Surg. 2018;53(1):112–9.
Axel L. Papillary muscles do not attach directly to the solid heart wall. Circulation. 2004;109(25):3145–8.
Nishitani S, et al. Development of helical myofiber tracts in the human fetal heart: analysis of myocardial fiber formation in the left ventricle from the late human embryonic period using diffusion tensor magnetic resonance imaging. J Am Heart Assoc. 2020;9(19): e016422.
Dierckx, H., et al. Intravoxel fibre structure of the left ventricular free wall and posterior left-right ventricular insertion site in canine myocardium using Q-ball imaging. In functional imaging and modeling of the heart. 2009. Berlin, Heidelberg: Springer Berlin Heidelberg.
Kocica MJ, Corno AF, Carreras-Costa F, Ballester-Rodes M, Moghbel MC, Cueva CN, Lackovic V, Kanjuh VI, Torrent-Guasp F. The helical ventricular myocardial band: global, three-dimensional, functional architecture of the ventricular myocardium. Eur J Cardio-Thorac Surg. 2006;29(Supplement_1):S21–40.
Vigmond EJ, et al. Solvers for the cardiac bidomain equations. Prog Biophys Mol Biol. 2008;96(1–3):3–18.
Salameh A, et al. Cyclic mechanical stretch induces cardiomyocyte orientation and polarization of the gap junction protein connexin43. Circ Res. 2010;106(10):1592–660.
Hughes RK, et al. Myocardial perfusion defects in hypertrophic cardiomyopathy mutation carriers. J Am Heart Assoc. 2021;10(15): e020227.
Zorzi A, et al. Nonischemic left ventricular scar as a substrate of life-threatening ventricular arrhythmias and sudden cardiac death in competitive athletes. Circ Arrhythm Electrophysiol. 2016;9(7): e004229.
Swift AJ, et al. LGE patterns in pulmonary hypertension do not impact overall mortality. JACC Cardiovasc Imaging. 2014;7(12):1209–17.
Freed BH, et al. Late gadolinium enhancement cardiovascular magnetic resonance predicts clinical worsening in patients with pulmonary hypertension. J Cardiovasc Magn Reson. 2012;14:11.
Roberts DE, Hersh LT, Scher AM. Influence of cardiac fiber orientation on wavefront voltage, conduction velocity, and tissue resistivity in the dog. Circ Res. 1979;44(5):701–12.
Vadakkumpadan F, Arevalo H, Trayanova NA. Patient-specific modeling of the heart: estimation of ventricular fiber orientations. JoVE. 2013. https://doi.org/10.3791/50125-v.
Boukens BJ, Potse M, Coronel R. Fibrosis and conduction abnormalities as basis for overlap of brugada syndrome and early repolarization syndrome. Int J Mol Sci. 2021;22(4):1570.
Kawamura M, et al. Idiopathic ventricular arrhythmia originating from the cardiac crux or inferior septum: epicardial idiopathic ventricular arrhythmia. Circ Arrhythm Electrophysiol. 2014;7(6):1152–8.
Li A, et al. Endocardial ablation of ventricular ectopic beats arising from the basal inferoseptal process of the left ventricle. Heart Rhythm. 2018;15(9):1356–62.
Liang JJ, et al. Electrocardiographic and electrophysiologic characteristics of idiopathic ventricular arrhythmias originating from the basal inferoseptal left ventricle. JACC Clin Electrophysiol. 2019;5(7):833–42.
Moulin K, Verzhbinsky IA, Maforo NG, Perotti LE, Ennis DB. Probing cardiomyocyte mobility with multi-phase cardiac diffusion tensor MRI. PLoS ONE. 2020;15(11): e0241996.
Ariga R, et al. Identification of myocardial disarray in patients with hypertrophic cardiomyopathy and ventricular arrhythmias. J Am Coll Cardiol. 2019;73(20):2493–502.
Lombaert H, Peyrat JM, Croisille P, Rapacchi S, Fanton L, Cheriet F, Clarysse P, Magnin I, Delingette H, Ayache N. Human atlas of the cardiac fiber architecture: study on a healthy population. IEEE Trans Med Imaging. 2012;31(7):1436–47.
Peyrat JM, Sermesant M, Pennec X, Delingette H, Xu C, McVeigh E, Ayache N. Towards a statistical atlas of cardiac fiber structure. In International Conference on Medical Image Computing and Computer-Assisted Intervention (pp. 297–304). Springer, Berlin, Heidelberg; 2006
Acknowledgements
The authors gratefully thank Pr. Frédéric Schnell for helpful discussions and suggestions on the topic. All clinical staff involved in the CADENCE project are gratefully acknowledged for their valuable contributions.
Funding
This study received financial support from the French government in the framework of the University of Bordeaux's France 2030 program / RRI "IMPACT" and the ANR (Agence Nationale de la Recherche) ANR-10-IAHU04-LIRYC.
Author information
Authors and Affiliations
Contributions
JM collected the hearts, prepared them for fixation, implement and optimize the sequences, carried out MR experiments and image reconstruction. MH, JR. L.B. and OB provided the human healthy heart thought the CADENCE project. FV, MC, DB, LP, VD, BQ collected the hearts at the hospital and prepared them for fixation. JM carried out MR experiments. PC, VO, YBP and JM conducted the MRI-study, developed the method (registration, DW processing, tractography pipeline) and processed the data. PP and PJ provided the sheep heart thought the UNMASC project. NPL, RW collected the heart and carried out microCT experiment. GR, VO and RW conducted the microCT-study, developed the method (structure tensor processing, tractography pipeline) and processed the data. MC conducted the histology-study. JRP conducted the EP-study. PC, VO, JRP, JD and EV drafted the manuscript. VO supervised the study, and participated in its design and coordination. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Human samples were derived from the Human donor program (providing access to heart samples from patients under cerebral death for scientific research purposes) approved by the ‘‘Agence Française de la Biomedecine’’ and with a written informed consent of the patient’s family. The experiment was conducted in accordance with the declaration of Helsinki and the institutional ethics committee.
Consent for publication
The authors declare that they have no competing of 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.
Supplementary figures and tables.
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 http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Cabanis, P., Magat, J., Rodriguez-Padilla, J. et al. Cardiac structure discontinuities revealed by ex-vivo microstructural characterization. A focus on the basal inferoseptal left ventricle region. J Cardiovasc Magn Reson 25, 78 (2023). https://doi.org/10.1186/s12968-023-00989-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12968-023-00989-y