Three-dimensional histology: tools and application to quantitative assessment of cell-type distribution in rabbit heart

Whole-Heart Histological and Electroanatomic Assessment of Postinfarction Cardiac Magnetic Resonance Imaging Scar and Conducting Channels

BACKGROUND

Cardiac magnetic resonance imaging (CMR)-defined ventricular scar and anatomic conduction channels (CMR-CCs) offer promise in delineating ventricular tachycardia substrate. No studies have validated channels with coregistered histology, nor have they ascertained the histological characteristics of deceleration zones (DZs) within these channels. We aimed to validate CMR scar and CMR-CCs with whole-heart histology and electroanatomic mapping in a postinfarction model.

METHODS

Five sheep underwent anteroseptal infarction. CMR (116±20 days post infarct) was postprocessed using ADAS-3D, varying pixel intensity thresholds (5545, 6040, 6535, and 7030). DZs were identified by electroanatomic mapping (129±12 days post infarct). Explanted hearts were sectioned and stained with Picrosirius red, and whole-heart histopathologic shells were generated. Scar topography as well as percentage fibrosis, adiposity, and remaining viable myocardium within 3 mm histological biopsies and within CMR-CCs were determined.

RESULTS

Using the standard 6040 thresholding, CMR had 83.8% accuracy for identifying histological scar in the endocardium (κ, 0.666) and 61.4% in the epicardium (κ, 0.276). Thirty-seven CMR-CCs were identified by varying thresholding; 23 (62%) were unique. DZs colocalized to 19 of 23 (83%) CMR-CCs. Twenty (87%) CMR-CCs were histologically confirmed. Within-channel histological fibrosis did not differ by the presence of DZs (P=0.242). Within-channel histological adiposity was significantly higher at sites with versus without DZs (24.1% versus 8.3%; P<0.001).

CONCLUSIONS

Postprocessed CMR-derived scars and channels were validated by histology and electroanatomic mapping. Regions of CMR-CCs at sites of DZs had higher adiposity but similar fibrosis than regions without DZs, suggesting that lipomatous metaplasia may contribute to arrhythmogenicity of postinfarction scar.

A groupwise registration and tractography framework for cardiac myofiber architecture description by diffusion MRI: An application to the ventricular junctions

Background

Knowledge of the normal myocardial–myocyte orientation could theoretically allow the definition of relevant quantitative biomarkers in clinical routine to diagnose heart pathologies. A whole heart diffusion tensor template representative of the global myofiber organization over species is therefore crucial for comparisons across populations. In this study, we developed a groupwise registration and tractography framework to resolve the global myofiber arrangement of large mammalian sheep hearts. To demonstrate the potential application of the proposed method, a novel description of sub-regions in the intraventricular septum is presented.

Methods

Three explanted sheep (ovine) hearts (size ~12×8×6 cm3, heart weight ~ 150 g) were perfused with contrast agent and fixative and imaged in a 9.4T magnet. A group-wise registration of high-resolution anatomical and diffusion-weighted images were performed to generate anatomical and diffusion tensor templates. Diffusion tensor metrics (eigenvalues, eigenvectors, fractional anisotropy …) were computed to provide a quantitative and spatially-resolved analysis of cardiac microstructure. Then tractography was performed using deterministic and probabilistic algorithms and used for different purposes: i) Visualization of myofiber architecture, ii) Segmentation of sub-area depicting the same fiber organization, iii) Seeding and Tract Editing. Finally, dissection was performed to confirm the existence of macroscopic structures identified in the diffusion tensor template.

Results

The template creation takes advantage of high-resolution anatomical and diffusion-weighted images obtained at an isotropic resolution of 150 μm and 600 μm respectively, covering ventricles and atria and providing information on the normal myocardial architecture. The diffusion metric distributions from the template were found close to the one of the individual samples validating the registration procedure. Small new sub-regions exhibiting spatially sharp variations in fiber orientation close to the junctions of the septum and ventricles were identified. Each substructure was defined and represented using streamlines. The existence of a fiber-bundles in the posterior junction was validated by anatomical dissection. A complex structural organization of the anterior junction in comparison to the posterior junction was evidenced by the high-resolution acquisition.

Conclusions

A new framework combining cardiac template generation and tractography was applied on the whole sheep heart. The framework can be used for anatomical investigation, characterization of microstructure and visualization of myofiber orientation across samples. Finally, a novel description of the ventricular junction in large mammalian sheep hearts was proposed.

Influence of respiration and tissue contact on ventricular substrate identification during high density mapping: results from an ovine infarct model

INTRODUCTION

Multi-electrode mapping (MEM) and automated point collection are important enhancements to substrate mapping in ventricular tachycardia ablation. The effects of tissue contact and respiration on electrogram voltage with differing depolarisation wavefronts with MEM catheters are unclear.

METHODS

Bipolar and unipolar voltages were collected from control (n=5) and infarcted (n=7) animals with a multi-spline MEM catheter. Electro-anatomic maps were created in sinus rhythm, and right and left ventricular pacing. Analysis was performed across three collection settings: standard (SS), respiratory-phase gating (RG), and electrode-tissue proximity (TP). Comparison was made to scar detected by cardiac MRI (cMRI).

RESULTS

Compared to SS and RG acquisition, median bipolar and unipolar voltages were higher using TP, regardless of the depolarization wavefront. In infarct animals, bipolar voltages were 30.7-50.5% higher for bipolar and 8.7-13.8% higher on unipolar voltages with TP, compared to SS. The effect of RG on bipolar and unipolar voltages was minimal. Percentage of local abnormal ventricular activities was not impacted by acquisition settings or wavefront direction in infarct animals. Compared with cMRI defined scar, all three acquisition settings overestimated scar area using standard voltage-based cutoffs. RG improved the low voltage area concordance with MRI by 1.6-5.1% whereas TP improved by 5.9-8.4%.

CONCLUSIONS

High density voltage mapping with a MEM catheter is influenced by point collection settings. Tissue contact filters reduced low voltage areas and improved agreement with cMRI fibrosis in infarcted ovine hearts. These findings have critical implications for optimising filter settings for high density substrate mapping in the left ventricle. This article is protected by copyright. All rights reserved.

Impact of Intraventricular Septal Fiber Orientation on Cardiac Electromechanical Function

Cardiac fiber direction is an important factor determining the propagation of electrical activity, as well as the development of mechanical force. In this article, we imaged the ventricles of several species with special attention to the intraventricular septum to determine the functional consequences of septal fiber organization. First, we identified a dual-layer organization of the fiber orientation in the intraventricular septum of ex vivo sheep hearts using diffusion tensor imaging at high field MRI. To expand the scope of the results, we investigated the presence of a similar fiber organization in five mammalian species (rat, canine, pig, sheep, and human) and highlighted the continuity of the layer with the moderator band in large mammalian species. We implemented the measured septal fiber fields in three-dimensional electromechanical computer models to assess the impact of the fiber orientation. The downward fibers produced a diamond activation pattern superficially in the right ventricle. Electromechanically, there was very little change in pressure volume loops although the stress distribution was altered. In conclusion, we clarified that the right ventricular septum has a downwardly directed superficial layer in larger mammalian species, which can have modest effects on stress distribution.

NEW & NOTEWORTHY A dual-layer organization of the fiber orientation in the intraventricular septum was identified in ex vivo hearts of large mammals. The RV septum has a downwardly directed superficial layer that is continuous with the moderator band. Electrically, it produced a diamond activation pattern. Electromechanically, little change in pressure volume loops were noticed but stress distribution was altered. Fiber distribution derived from diffusion tensor imaging should be considered for an accurate strain and stress analysis.

Quantitative collagen assessment in right ventricular myectomies from patients with tetralogy of Fallot

Aims

Patients with tetralogy of Fallot (TOF) are often affected by right ventricular fibrosis, which has been associated with arrhythmias. This study aimed to assess fibrosis distribution in right ventricular outflow tract (RVOT) myocardium of TOF patients to evaluate the utility of single histology-section analyses, and to explore the possibility of fibrosis quantification in unlabelled tissue by second harmonic generation imaging (SHGI) as an alternative to conventional histology-based assays.

Methods and results

We quantified fibrosis in 11 TOF RVOT samples, using a tailor-made automated image analysis method on Picrosirius red-stained sections. In a subset of samples, histology- and SHGI-based fibrosis quantification approaches were compared. Fibrosis distribution was highly heterogeneous, with significant and comparable variability between and within samples. We found that, on average, 67.8 mm² of 10 µm thick, histologically processed tissue per patient had to be analysed for accurate fibrosis quantification. SHGI provided data faster and on live tissue, additionally enabling quantification of collagen anisotropy.

Conclusion

Given the high intra-individual heterogeneity, fibrosis quantification should not be conducted on single sections of TOF RVOT myectomies. We provide an analysis algorithm for fibrosis quantification in histological images, which enables the required extended volume analyses in these patients.

3D Digital Pathology for a Chemical-Functional Analysis of Glomeruli in Health and Pathology

Determining the filtration function and biochemical status of kidney at the single glomerulus level remains hardly accessible, even from biopsies. Here, we provide evidence that IR spectro-microscopy is a suitable method to account for the filtration capacity of individual glomeruli along with related physio-pathological condition. A ∼4 μm voxel resolution 3D IR image reconstruction is built from consecutive tissue sections, thus, providing a 3D IR spectrum matrix of an individual glomerulus. The filtration capacity of glomeruli was quantitatively determined after BaSO₄ perfusion, and additional chemical data could be used to determined oxidative stress effects and fibrosis, thus, combining functional and biochemical information from the same 3D IR spectrum matrix. This analytical approach was applied on mice with unilateral ureteral obstruction (UUO) inducing chronic kidney disease. Compared to the healthy condition, UUO induced a significant drop in glomeruli filtration capacity (-17 ± 8% at day 4 and -48 ± 14% at day 14) and volume (36 ± 10% at day 4 and 67 ± 13% at day 14), along a significant increase of oxidative stress (+61 ± 19% at day 4 and +84 ± 17% at day 14) and a change in the lipid-to-protein ratio (-8.2 ± 3.6% at day 4 and -18.1 ± 5.9% at day 14). Therefore, IR spectro-microscopy might be developed as a new 3D pathology resource for analyzing functional and biochemical parameters of glomeruli.

Transformation diffusion reconstruction of three-dimensional histology volumes from two-dimensional image stacks

•

A method for 3D reconstruction of serial 2D histology image stacks is proposed.

•

Pre-alignment to an external pre-cut reference (blockface) prevents shape artifacts.

•

Formulated as diffusion of transformations from each slice to its neighbors.

•

Registrations replaced by much faster transformation operations.

Traditional histology is the gold standard for tissue studies, but it is intrinsically reliant on two-dimensional (2D) images. Study of volumetric tissue samples such as whole hearts produces a stack of misaligned and distorted 2D images that need to be reconstructed to recover a congruent volume with the original sample's shape. In this paper, we develop a mathematical framework called Transformation Diffusion (TD) for stack alignment refinement as a solution to the heat diffusion equation. This general framework does not require contour segmentation, is independent of the registration method used, and is trivially parallelizable. After the first stack sweep, we also replace registration operations by operations in the space of transformations, several orders of magnitude faster and less memory-consuming. Implementing TD with operations in the space of transformations produces our Transformation Diffusion Reconstruction (TDR) algorithm, applicable to general transformations that are closed under inversion and composition. In particular, we provide formulas for translation and affine transformations. We also propose an Approximated TDR (ATDR) algorithm that extends the same principles to tensor-product B-spline transformations. Using TDR and ATDR, we reconstruct a full mouse heart at pixel size 0.92 µm × 0.92 µm, cut 10 µm thick, spaced 20 µm (84G). Our algorithms employ only local information from transformations between neighboring slices, but the TD framework allows theoretical analysis of the refinement as applying a global Gaussian low-pass filter to the unknown stack misalignments. We also show that reconstruction without an external reference produces large shape artifacts in a cardiac specimen while still optimizing slice-to-slice alignment. To overcome this problem, we use a pre-cutting blockface imaging process previously developed by our group that takes advantage of Brewster's angle and a polarizer to capture the outline of only the topmost layer of wax in the block containing embedded tissue for histological sectioning.

Human-based approaches to pharmacology and cardiology: an interdisciplinary and intersectorial workshop

Both biomedical research and clinical practice rely on complex datasets for the physiological and genetic characterization of human hearts in health and disease. Given the complexity and variety of approaches and recordings, there is now growing recognition of the need to embed computational methods in cardiovascular medicine and science for analysis, integration and prediction. This paper describes a Workshop on Computational Cardiovascular Science that created an international, interdisciplinary and inter-sectorial forum to define the next steps for a human-based approach to disease supported by computational methodologies. The main ideas highlighted were (i) a shift towards human-based methodologies, spurred by advances in new in silico, in vivo, in vitro, and ex vivo techniques and the increasing acknowledgement of the limitations of animal models. (ii) Computational approaches complement, expand, bridge, and integrate in vitro, in vivo, and ex vivo experimental and clinical data and methods, and as such they are an integral part of human-based methodologies in pharmacology and medicine. (iii) The effective implementation of multi- and interdisciplinary approaches, teams, and training combining and integrating computational methods with experimental and clinical approaches across academia, industry, and healthcare settings is a priority. (iv) The human-based cross-disciplinary approach requires experts in specific methodologies and domains, who also have the capacity to communicate and collaborate across disciplines and cross-sector environments. (v) This new translational domain for human-based cardiology and pharmacology requires new partnerships supported financially and institutionally across sectors. Institutional, organizational, and social barriers must be identified, understood and overcome in each specific setting.

Resolving Fine Cardiac Structures in Rats with High-Resolution Diffusion Tensor Imaging

Cardiac architecture is fundamental to cardiac function and can be assessed non-invasively with diffusion tensor imaging (DTI). Here, we aimed to overcome technical challenges in ex vivo DTI in order to extract fine anatomical details and to provide novel insights in the 3D structure of the heart. An integrated set of methods was implemented in ex vivo rat hearts, including dynamic receiver gain adjustment, gradient system scaling calibration, prospective adjustment of diffusion gradients, and interleaving of diffusion-weighted and non-diffusion-weighted scans. Together, these methods enhanced SNR and spatial resolution, minimised orientation bias in diffusion-weighting, and reduced temperature variation, enabling detection of tissue structures such as cell alignment in atria, valves and vessels at an unprecedented level of detail. Improved confidence in eigenvector reproducibility enabled tracking of myolaminar structures as a basis for segmentation of functional groups of cardiomyocytes. Ex vivo DTI facilitates acquisition of high quality structural data that complements readily available in vivo cardiac functional and anatomical MRI. The improvements presented here will facilitate next generation virtual models integrating micro-structural and electro-mechanical properties of the heart.

Mapping cardiac microstructure of rabbit heart in different mechanical states by high resolution diffusion tensor imaging: A proof-of-principle study

Myocardial microstructure and its macroscopic materialisation are fundamental to the function of the heart. Despite this importance, characterisation of cellular features at the organ level remains challenging, and a unifying description of the structure of the heart is still outstanding. Here, we optimised diffusion tensor imaging data to acquire high quality data in ex vivo rabbit hearts in slack and contractured states, approximating diastolic and systolic conditions. The data were analysed with a suite of methods that focused on different aspects of the myocardium. In the slack heart, we observed a similar transmural gradient in helix angle of the primary eigenvector of up to 23.6°/mm in the left ventricle and 24.2°/mm in the right ventricle. In the contractured heart, the same transmural gradient remained largely linear, but was offset by up to +49.9° in the left ventricle. In the right ventricle, there was an increase in the transmural gradient to 31.2°/mm and an offset of up to +39.0°. The application of tractography based on each eigenvector enabled visualisation of streamlines that depict cardiomyocyte and sheetlet organisation over large distances. We observed multiple V- and N-shaped sheetlet arrangements throughout the myocardium, and insertion of sheetlets at the intersection of the left and right ventricle. This study integrates several complementary techniques to visualise and quantify the heart's microstructure, projecting parameter representations across different length scales. This represents a step towards a more comprehensive characterisation of myocardial microstructure at the whole organ level.

Three-dimensional visualisation of soft biological structures by X-ray computed micro-tomography

Whereas the two-dimensional (2D) visualisation of biological samples is routine, three-dimensional (3D) imaging remains a time-consuming and relatively specialised pursuit. Current commonly adopted techniques for characterising the 3D structure of non-calcified tissues and biomaterials include optical and electron microscopy of serial sections and sectioned block faces, and the visualisation of intact samples by confocal microscopy or electron tomography. As an alternative to these approaches, X-ray computed micro-tomography (microCT) can both rapidly image the internal 3D structure of macroscopic volumes at sub-micron resolutions and visualise dynamic changes in living tissues at a microsecond scale. In this Commentary, we discuss the history and current capabilities of microCT. To that end, we present four case studies to illustrate the ability of microCT to visualise and quantify: (1) pressure-induced changes in the internal structure of unstained rat arteries, (2) the differential morphology of stained collagen fascicles in tendon and ligament, (3) the development of Vanessa cardui chrysalises, and (4) the distribution of cells within a tissue-engineering construct. Future developments in detector design and the use of synchrotron X-ray sources might enable real-time 3D imaging of dynamically remodelling biological samples.

How to improve the overall quality of cardiac morphometric data

By the mid-1990s, experts realized that drugs leading to improved ventricular remodeling were doing something remarkable in cardiac patients. The "age of cardiac remodeling" had begun. This created an experimental need for high-quality assessment of changes in cardiac tissue composition, including myocyte shape, myocardial fibrosis/collagen, and vascular remodeling. Many working in the field today have little or no training related to recognition of fixation artifacts or common errors associated with quantitative morphology. Unfortunately, such skills had become somewhat of a lost art during the ages of cardiac physiology in the mid-20th century and molecular biology, gaining prominence by the mid-1970s. Consequently, cardiac remodeling studies today are often seriously flawed to the point where data are not reproducible and subsequent researchers may be chasing the molecular basis of a nonexistent or erroneous phenotype. The current unacceptably high incidence of irreproducible data is a serious waste of time and resources as recently noted in comments by the National Institutes of Health director. The goal of this "how to" article is to share some lessons I have learned during nearly 40 years of assessing morphological changes in the heart. It is possible for any laboratory to routinely publish highly reproducible morphological data that stand the test of time and contribute to our fundamental knowledge of cardiac remodeling and the molecular mechanisms that drive it.

Three-dimensional cardiac computational modelling: methods, features and applications

The combination of computational models and biophysical simulations can help to interpret an array of experimental data and contribute to the understanding, diagnosis and treatment of complex diseases such as cardiac arrhythmias. For this reason, three-dimensional (3D) cardiac computational modelling is currently a rising field of research. The advance of medical imaging technology over the last decades has allowed the evolution from generic to patient-specific 3D cardiac models that faithfully represent the anatomy and different cardiac features of a given alive subject. Here we analyse sixty representative 3D cardiac computational models developed and published during the last fifty years, describing their information sources, features, development methods and online availability. This paper also reviews the necessary components to build a 3D computational model of the heart aimed at biophysical simulation, paying especial attention to cardiac electrophysiology (EP), and the existing approaches to incorporate those components. We assess the challenges associated to the different steps of the building process, from the processing of raw clinical or biological data to the final application, including image segmentation, inclusion of substructures and meshing among others. We briefly outline the personalisation approaches that are currently available in 3D cardiac computational modelling. Finally, we present examples of several specific applications, mainly related to cardiac EP simulation and model-based image analysis, showing the potential usefulness of 3D cardiac computational modelling into clinical environments as a tool to aid in the prevention, diagnosis and treatment of cardiac diseases.

Interrogation of living myocardium in multiple static deformation states with diffusion tensor and diffusion spectrum imaging

Diffusion tensor magnetic resonance imaging (MRI) reveals valuable insights into tissue histo-anatomy and microstructure, and has steadily gained traction in the cardiac community. Its wider use in small animal cardiac imaging in vivo has been constrained by its extreme sensitivity to motion, exaggerated by the high heart rates usually seen in rodents. Imaging of the isolated heart eliminates respiratory motion and, if conducted on arrested hearts, cardiac pulsation. This serves as an important intermediate step for basic and translational studies. However, investigating the micro-structural basis of cardiac deformation in the same heart requires observations in different deformation states.

Here, we illustrate the imaging of isolated rat hearts in three mechanical states mimicking diastole (cardioplegic arrest), left-ventricular (LV) volume overload (cardioplegic arrest plus LV balloon inflation), and peak systole (lithium-induced contracture). An optimised MRI-compatible Langendorff perfusion setup with the radio-frequency (RF) coil integrated into the wet chamber was developed for use in a 9.4T horizontal bore scanner. Signal-to-noise ratio improved significantly, by 75% compared to a previous design with external RF coil, and stability tests showed no significant changes in mean T₁, T₂ or LV wall thickness over a 170 min period. In contracture, we observed a significant reduction in mean fractional anisotropy from 0.32 ± 0.02 to 0.28 ± 0.02, as well as a significant rightward shift in helix angles with a decrease in the proportion of left-handed fibres, as referring to the locally prevailing cell orientation in the heart, from 24.9% to 23.3%, and an increase in the proportion of right-handed fibres from 25.5% to 28.4%. LV overload, in contrast, gave rise to a decrease in the proportion of left-handed fibres from 24.9% to 21.4% and an increase in the proportion of right-handed fibres from 25.5% to 26.0%. The modified perfusion and coil setup offers better performance and control over cardiac contraction states.

We subsequently performed high-resolution diffusion spectrum imaging (DSI) and 3D whole heart fibre tracking in fixed ex vivo rat hearts in slack state and contracture. As a model-free method, DSI augmented the measurements of water diffusion by also informing on multiple intra-voxel diffusion orientations and non-Gaussian diffusion. This enabled us to identify the transition from right- to left-handed fibres from the subendocardium to the subepicardium, as well as voxels in apical regions that were traversed by multiple fibres. We observed that both the mean generalised fractional anisotropy and mean kurtosis were lower in hearts in contracture compared to the slack state, by 23% and 9.3%, respectively. While its heavy acquisition burden currently limits the application of DSI in vivo, ongoing work in acceleration techniques may enable its use in live animals and patients. This would provide access to the as yet unexplored dimension of non-Gaussian diffusion that could serve as a highly sensitive marker of cardiac micro-structural integrity.