Why SIT works: normal function despite typical myofiber pattern in Situs Inversus Totalis (SIT) hearts derived by shear-induced myofiber reorientation

Establishment of Cardiac Laterality

Formation of the vertebrate heart with its complex arterial and venous connections is critically dependent on patterning of the left-right axis during early embryonic development. Abnormalities in left-right patterning can lead to a variety of complex life-threatening congenital heart defects. A highly conserved pathway responsible for left-right axis specification has been uncovered. This pathway involves initial asymmetric activation of a nodal signaling cascade at the embryonic node, followed by its propagation to the left lateral plate mesoderm and activation of left-sided expression of the Pitx2 transcription factor specifying visceral organ asymmetry. Intriguingly, recent work suggests that cardiac laterality is encoded by intrinsic cell and tissue chirality independent of Nodal signaling. Thus, Nodal signaling may be superimposed on this intrinsic chirality, providing additional instructive cues to pattern cardiac situs. The impact of intrinsic chirality and the perturbation of left-right patterning on myofiber organization and cardiac function warrants further investigation. We summarize recent insights gained from studies in animal models and also some human clinical studies in a brief overview of the complex processes regulating cardiac asymmetry and their impact on cardiac function and the pathogenesis of congenital heart defects.

Using Systolic Local Mechanical Load to Predict Fiber Orientation in Ventricles

A simple rule adopted for myofiber reorientation in the ventricles is pursued by taking the microscopic branching network of myocytes into account. The macroscopic active tension generated on the microscopic branching structure is modeled by a multidirectional active stress tensor, which is defined as a function of the strains in the branching directions. In our reorientation algorithm, the principal direction of the branching network is updated so that it turns in the direction of greater active tension in the isovolumetric systole. Updates are performed step-by-step after the mechanical equilibrium has been attained with the current fiber structure. Starting from a nearly flat distribution of the principal fiber orientation along the circumferential direction, the reoriented fiber helix angles range from 70 to 40° at epicardium and from 60 to 80° at endocardium, in agreement with experimental observations. The helical ventricular myocardial band of Torrent-Guasp’s model and the apical spiral structure of Rushmer’s model are also reconstructed by our algorithm. Applying our algorithm to the infarcted ventricle model, the fiber structure near the infarcted site is remodeled so that the helix angle becomes steeper with respect to the circumferential direction near the epicardial surface. Based on our numerical analysis, we draw the following conclusions. (i) The multidirectional active tension based on the microscopic branching network is potentially used to seek tighter connection with neighboring aggregates. (ii) The thickening and thinning transitions in response to active tension in each myocyte allow the macroscopic principal fiber orientation of the microscopic branching network to move toward the direction of greater active tension. (iii) The force–velocity relationship is the key factor in transferring the fiber shortening strain to the magnitude of active tensions used in the myofiber reorientation. (iv) The algorithm naturally leads to homogeneity in the macroscopic active tension and the fiber shortening strain, and results in near-optimal pumping performance. (v) However, the reorientation mechanism may degrade the pumping performance if there is severely inhomogeneous contractility resulting from infarction. Our goal is to provide a tool to predict the fiber architecture of various heart disease patients for numerical simulations of their treatment plans.

Optical coherence tomography imaging of cardiac substrates

Cardiovascular disease is the leading cause of morbidity and mortality in the United States. Knowledge of a patient's heart structure will help to plan procedures, potentially identifying arrhythmia substrates, critical structures to avoid, detect transplant rejection, and reduce ambiguity when interpreting electrograms and functional measurements. Similarly, basic research of numerous cardiac diseases would greatly benefit from structural imaging at cellular scale. For both applications imaging on the scale of a myocyte is needed, which is approximately 100 µm × 10 µm. The use of optical coherence tomography (OCT) as a tool for characterizing cardiac tissue structure and function has been growing in the past two decades. We briefly review OCT principles and highlight important considerations when imaging cardiac muscle. In particular, image penetration, tissue birefringence, and light absorption by blood during in vivo imaging are important factors when imaging the heart with OCT. Within the article, we highlight applications of cardiac OCT imaging including imaging heart tissue structure in small animal models, quantification of myofiber organization, monitoring of radiofrequency ablation (RFA) lesion formation, structure-function analysis enabled by functional extensions of OCT and multimodal analysis and characterizing important substrates within the human heart. The review concludes with a summary and future outlook of OCT imaging the heart, which is promising with progress in optical catheter development, functional extensions of OCT, and real time image processing to enable dynamic imaging and real time tracking during therapeutic procedures.

Application of an Adaptive Polynomial Chaos Expansion on Computationally Expensive Three-Dimensional Cardiovascular Models for Uncertainty Quantification and Sensitivity Analysis

When applying models to patient-specific situations, the impact of model input uncertainty on the model output uncertainty has to be assessed. Proper uncertainty quantification (UQ) and sensitivity analysis (SA) techniques are indispensable for this purpose. An efficient approach for UQ and SA is the generalized polynomial chaos expansion (gPCE) method, where model response is expanded into a finite series of polynomials that depend on the model input (i.e., a meta-model). However, because of the intrinsic high computational cost of three-dimensional (3D) cardiovascular models, performing the number of model evaluations required for the gPCE is often computationally prohibitively expensive. Recently, Blatman and Sudret (2010, “An Adaptive Algorithm to Build Up Sparse Polynomial Chaos Expansions for Stochastic Finite Element Analysis,” Probab. Eng. Mech., 25(2), pp. 183–197) introduced the adaptive sparse gPCE (agPCE) in the field of structural engineering. This approach reduces the computational cost with respect to the gPCE, by only including polynomials that significantly increase the meta-model’s quality. In this study, we demonstrate the agPCE by applying it to a 3D abdominal aortic aneurysm (AAA) wall mechanics model and a 3D model of flow through an arteriovenous fistula (AVF). The agPCE method was indeed able to perform UQ and SA at a significantly lower computational cost than the gPCE, while still retaining accurate results. Cost reductions ranged between 70–80% and 50–90% for the AAA and AVF model, respectively.

Assessment and comparison of left ventricular shear in normal and situs inversus totalis hearts by means of magnetic resonance tagging

Situs inversus totalis (SIT) is characterized by complete mirroring of gross cardiac anatomy and position combined with an incompletely mirrored myofiber arrangement, being normal at the apex but inverted at the base of the left ventricle (LV). This study relates myocardial structure to mechanical function by analyzing and comparing myocardial deformation patterns of normal and SIT subjects, focusing especially on circumferential-radial shear. In nine control and nine SIT normotensive human subjects, myocardial deformation was assessed from magnetic resonance tagging (MRT) image sequences of five LV short-axis slices. During ejection, no significant difference in either circumferential shortening (εcc) or its axial gradient (Δεcc) is found between corresponding LV levels in control and SIT hearts. Circumferential-radial shear (εcr) has a clear linear trend from apex-to-base in controls, while in SIT it hovers close to zero at all levels. Torsion as well as axial change in εcr (Δεcr) is as in controls in apical sections of SIT hearts but deviates significantly towards the base, changing sign close to the LV equator. Interindividual variability in torsion and Δεcr values is higher in SIT than in controls. Apex-to-base trends of torsion and Δεcr in SIT, changing sign near the LV equator, further substantiate a structural transition in myofiber arrangement close to the LV equator itself. Invariance of εcc and Δεcc patterns between controls and SIT subjects shows that normal LV pump function is achieved in SIT despite partial mirroring of myocardial structure leading to torsional and shear patterns that are far from normality.

Effects of activation pattern and active stress development on myocardial shear in a model with adaptive myofiber reorientation

It has been hypothesized that myofiber orientation adapts to achieve a preferred mechanical loading state in the myocardial tissue. Earlier studies tested this hypothesis in a combined model of left ventricular (LV) mechanics and remodeling of myofiber orientation in response to fiber cross-fiber shear, assuming synchronous timing of activation and uniaxial active stress development. Differences between computed and measured patterns of circumferential-radial shear strain Ecr were assumed to be caused by limitations in either the LV mechanics model or the myofiber reorientation model. Therefore, we extended the LV mechanics model with a physiological transmural and longitudinal gradient in activation pattern and with triaxial active stress development. We investigated the effects on myofiber reorientation, LV function, and deformation. The effect on the developed pattern of the transverse fiber angle αt,0 and the effect on global pump function were minor. Triaxial active stress development decreased amplitudes of Ecr towards values within the experimental range and resulted in a similar base-to-apex gradient during ejection in model computed and measured Ecr. The physiological pattern of mechanical activation resulted in better agreement between computed and measured strain in myofiber direction, especially during isovolumic contraction phase and first half of ejection. In addition, remodeling was favorable for LV pump and myofiber function. In conclusion, the outcome of the combined model of LV mechanics and remodeling of myofiber orientation is found to become more physiologic by extending the mechanics model with triaxial active stress development and physiological activation pattern.

Extracting three-dimensional orientation and tractography of myofibers using optical coherence tomography

Abnormal changes in orientation of myofibers are associated with various cardiac diseases such as arrhythmia, irregular contraction, and cardiomyopathy. To extract fiber information, we present a method of quantifying fiber orientation and reconstructing three-dimensional tractography of myofibers using optical coherence tomography (OCT). A gradient based algorithm was developed to quantify fiber orientation in three dimensions and particle filtering technique was employed to track myofibers. Prior to image processing, three-dimensional image data set were acquired from all cardiac chambers and ventricular septum of swine hearts using OCT system without optical clearing. The algorithm was validated through rotation test and comparison with manual measurements. The experimental results demonstrate that we are able to visualize three-dimensional fiber tractography in myocardium tissues.