Mechanics Meets Echocardiographic Imaging: Investigation on the Principal Strain Lines in Human Left Ventricle

. From a mechanical point of view, the most intense work is performed along the pattern from point 2 to point 3, that is, along the systolic phase, when pressures are high and muscle contraction too. Typically, critical behaviors of the ventricular function are evidenced in this phase, and mechanics can suggest which are good indicators of cardiac function. A relevant condition about these indicators is the possibility to catch them through noninvasive analyses.




A329170_1_En_3_Fig1_HTML.gif


Fig. 1
Cartoon sketching the phases of the cardiac cycle of a normal human subject. 1 Mitral valve closes; isovolumetric contraction. 2 Aortic valve opens; ejection. 3 Aortic valve closes; relaxation. 4 Mitral valve opens, filling. The green area represents the stroke work

A well-known example is ventricular torsion [23]. The role played by the LV torsional rotation with respect to LV ejection and filling was only recently recognized by application of speckle tacking echocardiography, whose output includes, among the other things, the pattern of ventricular torsion along the cardiac cycle [4, 5, 1012, 16]. As ventricular torsion is altered when a few pathologies are present (see [3, 13, 18, 22, 24, 28]), it can be used as an indicator of cardiac function which can be noninvasively investigated through 3-dimensional speckle tracking echocardiography (3DSTE).

Detection of principal strain lines in LV may emerge as another possible noninvasive tool to discriminate among different LVs as well as a tool that can help clinicians to identify cardiac diseases at the early stages. On the other hand, and differently from ventricular torsion, PSLs are not delivered as output by 3DSTE devices, and a post-processing analysis of 3DSTE data is needed to identify them, based on concepts borrowed from continuum mechanics. In [17, 20], it was initially proposed to look at PSLs to identify muscle fiber architecture on the endocardial surface. Therein, the echocardiographic analysis was limited to the endocardial surface, and it was noted as due to the high contractions suffered by muscle fibres along the systolic phase, PSLs may well determine just muscle fiber directions. Successively, in [8] an accurate protocol of measurement of PSLs was proposed, tested, and succesfully verified through a computational model. The conclusions of this last work were partially in contrast with the ones in [20]. It was demonstrated, firstly, that on the endocardial surface of healthy LVs, primary strain lines identify circumferential material directions; secondly, that on the epicardial surface primary strain lines are similar to muscle fiber directions. In [6, 8, 9], a comparison between a real human LV and a corresponding model was implemented by the same Authors; the conclusions of [8] were confirmed, and made precise through a statistical analysis involving real and computational data.

What is emerging, even if further investigations are needed, is that endocardial PSLs coincide with circumferential material lines, due to the relevant stiffening effect of the circumferential material lines when high pressures are involved, as it occurs along the systolic phase, and to the capacity of the same material fibers to contrast the LV dilation. It would mean that these visible functional strain lines are related to the capacity of elastic response of the cardiac tissue to the high systolic pressure, and that it might be important to follow this pattern when, due to pathological conditions, this capacity is missing.



2 Continuum Cardio-Mechanics


Typically, when mechanics is applied to biology, it is named biomechanics; we use cardio-mechanics to mean that specific branch of mechanics which has been successfully applied to the analysis and investigation of the cardiovascular system, whose center is the heart pump. Discuss the contribution of cardio-mechanics to clinics is beyond our aims (refer to [7] for extended references); we only aim to shortly discuss two different approaches of mechanics to clinics, with specific reference to the heart pump.

The first approach is based on the continuum analysis of the heart, aimed to model the cardiac activity. It starts from the construction of anatomical models of the heart and of constitutive model describing the passive and active material response of cardiac tissues. Electromechanical interaction are sometimes accounted for, if one has interest to also investigate cardiac electrophysiology. Typically, these models are implemented within a finite element code. Critical points are the constitutive models of the tissue, being cardiac tissue highly non homogeneous and contracting, due to an electrophysiological stimulus, and anatomical data about muscle architecture, which strongly influences the mechanical performances of the heart [14, 15, 27]. Once the model is complete, specific cardiac diseases can be included within the modeling, and the consequences on the heart activity studied. Typical examples are given by the investigations on the heart remodeling due to left anterior descending artery occlusion, to ventricular hypertrophy, to aortic stenosis, well discussed in literature [10, 18, 24]. Likewise, the model can be improved to account consequences of an infarct in different places of the left ventricular walls, with the aim of studying the ability of the pump to get on or not its work [13].

Here, we are more interested in the second approach, which starts from an analysis of real data extracted from the heart through appropriate tools such as Magnetic Risonance Image (MRI) and 3DSTE. Of course, such data are only concerning heart kinematics, and say nothing about stresses within heart walls. However, as it is recently shown, an accurate and careful analysis of real data allows to deeply investigate on heart. Let us note that the interest in in-vivo myocardial deformation dates back to 
$'80$
; in [26], the normal in-vivo three-dimensional finite strains were studied in dogs, through the application of appropriate markers whose coordinated were followed along the cardiac cycle through high-speed biplane cineradiography. Of course, the analysis was highly invasive. The recent techniques of visualization realized by 3DSTE make possible follow the coordinates of natural markers, automatically identified by the device, appropriately supported by an operator.

A typical example of the outcomes of a deep analysis on 3DSTE real data comes from PSLs, which are the core of this contribution. In mechanics, it is well-known that stresses and strains within a body are limited above and below by their principal counterparts; this allows for the discussion and verification of the mechanical state of that body. Moreover, the principal stress and strain lines (which are the same only when special symmetry conditions are verified) determine the directions where the largest strains and/or stresses are to be expected. Due to these characteristics, the mechanics of fiber-reinforced bodies are often based on the detection of the principal strain lines and, wherever needed, fiber architecture is conceived in order to make the fiber lines coincide with the PSLs. Fibers make a tissue highly anisotropic; hence, principal strain and stress lines may be distinct. Whereas principal strains can be measured starting with the analysis of tissue motion, being only dependent on the three-dimensional strain state of the tissue, principal stresses can only be inferred. Thus, the PSL have a predominant role where the analysis of the mechanics of a body is concerned, and can reveal which are the lines where largest strains are expected, and how they change when diseases occur.

Key point is the evaluation at any place within the body of the nonlinear strain tensor 
${\mathbf{C}}$
, whose eigenpairs (eigenvalues-eigenvectors) deliver principal strains and PSLs, respectively. Fixed a body, identified with the region 
${\mathcal{B}}$
of the three-dimensional Euclidean space 
${\mathcal{E}}$
it occupies at a time t o denoted as reference configuration of the body, we are interested in following the motion of the body at any time 
$t\in{\mathcal{I}}\subset{\mathcal{R}}$
, with the time interval 
${\mathcal{I}}$
identifying the duration of a human cardiac cycle (hence, different from subject to subject, as it is discussed later). The displacement field 
${\mathbf{u}}$
, that is a map from 
${\mathcal{B}}\times{\mathcal{T}}\to{\mathcal{V}}=T{\mathcal{E}}$
, delivers at any time and for any point 
$y\in{\mathcal{B}}$
the position 
$p(y,t)$
of that point at that time: 
$p(y,t)=y+{\mathbf{u}}(y,t)$
. Strains are related to displacement gradients within the body; precisely, it can be shown as, introduced the deformation gradient 
${\mathbf{F}}=\nabla\,p={\mathbf{I}}+\nabla{\mathbf{u}}$
, the nonlinear Cauchy-Green strain tensor is



$${\mathbf{C}}={\mathbf{F}}^T{\mathbf{F}} = {\mathbf{I}} + \nabla{\mathbf{u}}+\nabla{\mathbf{u}}^T+\nabla{\mathbf{u}}^T\nabla{\mathbf{u}},$$

(1)
being 
${\mathbf{I}}$
the identity tensor in 
${\mathcal{V}}$
. In general, 
${\mathbf{C}}$
is a three-dimensional tensor, describing the strain state at any point y and time t of the body. If there is within the body a distinguished surface 
${\mathcal{S}}$
, whose unit normal field is described by the unit vector field 
${\mathbf{n}}$
, the corresponding surface strain tensor 
$\hat{\mathbf{C}}$
can be obtained through a preliminary projection of 
${\mathbf{C}}$
onto that surface. The projector 
${\mathbf{P}}={\mathbf{I}}-{\mathbf{n}}\otimes{\mathbf{n}}$
, leads to the following definition:



$$\hat{\mathbf{C}}={\mathbf{P}}{\mathbf{C}}{\mathbf{P}}.$$

(2)
It is expected that 
$\hat{\mathbf{C}}$
will represent a plane strain state, hence, that it will have a zero eigenvalue corresponding to the eigenvector 
${\mathbf{n}}$
. The primary strain lines on the surface will be the streamlines of the eigenvector 
${\mathbf{c}}_2$
, which lies on the surface and corresponds to the smallest non-zero eigenvalues; the secondary strain lines are the streamlines corresponding to the eigenvector 
${\mathbf{c}}_3$
. Of course, when the strain tensor 
${\mathbf{C}}$
is ab initio evaluated from surface deformation gradients 
$\check{\mathbf{F}}={\mathbf{P}}{\mathbf{F}}{\mathbf{P}}$
, it naturally arises as a plane tensor.


3 Speckle Tracking Echocardiography


Speckle tracking echocardiography (STE) is an application of pattern-matching technology to ultrasound cine data and is based on the tracking of the ‘speckles’ in a 2D plane or in a 3D volume (2DSTE and 3DSTE, respectively). Speckles are disturbances in ultrasounds caused by reflections in the ultrasound beam: each structure in the body has a unique speckle pattern that moves with tissue (Fig. 2, left panel). A square or cubic template image is created using a local myocardial region in the starting frame of the image data. The size of the template image is around 1 cm2 in 2D or 1 cm3 in 3D. In the successive frame, the algorithm identifies the local speckle pattern that most closely matches the template (see [29] for further details). A displacement vector is created using the location of the template and the matching image in the subsequent frame. Multiple templates can be used to observe displacements of the entire myocardium. By using hundreds of these samples in a single image, it is possible to provide regional information on the displacement of the LV walls, and thus, other parameters such as strain, rotation, twist and torsion can be derived.



A329170_1_En_3_Fig2_HTML.gif


Fig. 2
Speckles moving with tissue as viewed through STE (left); the apical four chamber view (A); the second apical view orthogonal to plane A (B); three short-axis planes (C), in the apical region (C1), in the mid-ventricle (C2), and at the basal portion of the LV (C3) (right) (unmodified from the original ARTIDA image)

Echocardiographic examinations were performed with an Aplio-Artida ultrasound system (Toshiba Medical Systems Co, Tochigi, Japan). Full-volume ECG-gated 3D data sets were acquired from apical positions using a 1–4 MHz 3D matrix array transducer to visualize the entire LV in a volumetric image. To obtain these 3D data sets, four or six sectors were scanned from consecutive cardiac cycles and combined to provide a larger pyramidal volume covering the entire LV. The final LV geometry was reconstructed by starting from a set of 6 homologous landmarks (see Fig. 2), manually detected by the operator for all subjects under study. The manual detection for a given set of landmarks is crucial because it allows recording spatial coordinates in perfectly comparable anatomical structures of different subjects (following a homology principle). The results of our 3DSTE system is a time-sequence of shapes, each constituted by 1297 landmarks, assumed to be homologous, for both the epicardial and endocardial surfaces, positioned along 36 horizontal circles, each comprised of 36 landmarks, plus the apex (see Fig. 3).

Only gold members can continue reading. Log In or Register to continue

Stay updated, free articles. Join our Telegram channel

Jun 14, 2017 | Posted by in GENERAL SURGERY | Comments Off on Mechanics Meets Echocardiographic Imaging: Investigation on the Principal Strain Lines in Human Left Ventricle

Full access? Get Clinical Tree

Get Clinical Tree app for offline access