MR imaging of left-ventricular function : novel image acquisition and analysis techniques.

Many cardiac diseases, such as myocardial ischemia, secondary to coronary artery disease, may be identified and localized through the analysis of cardiac deformations. Early efforts for quantifying ventricular wall motion used surgical implantation and tracking of radiopaque markers with X-ray imaging in canine hearts [1]. Such techniques are invasive and affect the regional motion pattern of the ventricular wall during the marker tracking process and, clearly are not feasible clinically. Noninvasive imaging techniques are vital and have been widely applied to the clinic. MRI is a noninvasive imaging technique with the capability to monitor and assess the progression of cardiovascular diseases (CVD) so that effective procedures for the care and treatment of patients can be developed by physicians and researchers. It is capable of providing 3D analysis of global and regional cardiac function with great accuracy and reproducibility. In the past few years, numerous efforts have been devoted to cardiac motion recovery and deformation analysis from MR imaging sequences. In order to assess cardiac function, there are two categories of indices that are used: global and regional indices. Global indices include ejection fraction, cavity volume, and myocardial mass [2]. They are important indices for cardiac disease diagnosis. However, these global indices are not specific for regional analysis. A quantitative assessment of regional parameters may prove beneficial for the diagnosis of disease and evaluation of severity and the quantification of treatment [3]. Local measures, such as wall deformation and strain in all regions of the heart, can provide objective regional quantification of ventricular wall function and relate to the location and extent of ischemic injury. This dissertation is concerned with the development of novel MR imaging techniques and image postprocessing algorithms to analyze left ventricular deformations. A novel pulse sequence, termed Orthogonal CSPAMM (OCSPAMM), has been proposed which results in the same acquisition time as SPAMM for 2D deformation estimation while keeping the main advantages of CSPAMM [4,5]: i.e., maintaining tag contrast through-out the ECG cycle. Different from CSPAMM, in OCSPAMM the second tagging pulse orientation is rotated 90 degrees relative to the first one so that motion information can be obtained simultaneously in two directions. This reduces the acquisition time by a factor of two as compared to the traditional CSPAMM, in which two separate imaging sequences are applied per acquisition. With the application of OCSPAMM, the effect of tag fading encountered in SPAMM tagging due to Tl relaxation is mitigated and tag deformations can be visualized for the entire cardiac cycle, including diastolic phases. A multilevel B-spline fitting method (MBS) has been proposed which incorporates phase-based displacement information for accurate calculation of 2D motion and strain from tagged MRI [6, 7]. The proposed method combines the advantages of continuity and smoothness of MBS, and makes use of phase information derived from tagged MR images. Compared to previous 2D B-spline-based deformation analysis methods, MBS has the following advantages: 1) It can simultaneously achieve a smooth deformation while accurately approximating the given data set; 2) Computationally, it is very fast; and 3) It can produce more accurate deformation results. Since the tag intersections (intersections between two tag lines) can be extracted accurately and are more or less distributed evenly over the myocardium, MBS has proven effective for 2D cardiac motion tracking. To derive phase-based displacements, 2D HARP and SinMod analysis techniques [8,9] were employed. By producing virtual tags from HARP /SinMod and calculating intersections of virtual tag lines, more data points are obtained. In the reference frame, virtual tag lines are the isoparametric curves of an undeformed 2D B-spline model. In subsequent frames, the locations of intersections of virtual tag lines over the myocardium are updated with phase-based displacement. The advantage of the technique is that in acquiring denser myocardial displacements, it uses both real and virtual tag line intersections. It is fast and more accurate than 2D HARP and SinMod tracking. A novel 3D sine wave modeling (3D SinMod) approach for automatic analysis of 3D cardiac deformations has been proposed [10]. An accelerated 3D complementary spatial modulation of magnetization (CSPAMM) tagging technique [11] was used to acquire complete 3D+t tagged MR data sets of the whole heart (3 dynamic CSPAMM tagged MRI volume with tags in different orientations), in-vivo, in 54 heart beats and within 3 breath-holds. In 3D SinMod, the intensity distribution around each pixel is modeled as a cosine wave front. The principle behind 3D SinMod tracking is that both phase and frequency for each voxel are determined directly from the frequency analysis and the displacement is calculated from the quotient of phase difference and local frequency. The deformation fields clearly demonstrate longitudinal shortening during systole. The contraction of the LV base towards the apex as well as the torsional motion between basal and apical slices is clearly observable from the displacements. 3D SinMod can automatically process the image data to derive measures of motion, deformations, and strains between consecutive pair of tagged volumes in 17 seconds. Therefore, comprehensive 4D imaging and postprocessing for determination of ventricular function is now possible in under 10 minutes. For validation of 3D SinMod, 7 3D+t CSPAMM data sets of healthy subjects have been processed. Comparison of mid-wall contour deformations and circumferential shortening results by 3D SinMod showed good agreement with those by 3D HARP. Tag lines tracked by the proposed technique were also compared with manually delineated ones. The average errors calculated for the systolic phase of the cardiac cycles were in the sub-pixel range

Subendocardial contractile impairment in chronic ischemic myocardium: assessment by strain analysis of 3T tagged CMR

<p>Abstract</p> <p>Background</p> <p>The purpose of this study was to quantify myocardial strain on the subendocardial and epicardial layers of the left ventricle (LV) using tagged cardiovascular magnetic resonance (CMR) and to investigate the transmural degree of contractile impairment in the chronic ischemic myocardium.</p> <p>Methods</p> <p>3T tagged CMR was performed at rest in 12 patients with severe coronary artery disease who had been scheduled for coronary artery bypass grafting. Circumferential strain (C-strain) at end-systole on subendocardial and epicardial layers was measured using the short-axis tagged images of the LV and available software (Intag; Osirix). The myocardial segment was divided into stenotic and non-stenotic segments by invasive coronary angiography, and ischemic and non-ischemic segments by stress myocardial perfusion scintigraphy. The difference in C-strain between the two groups was analyzed using the Mann-Whitney U-test. The diagnostic capability of C-strain was analyzed using receiver operating characteristics analysis.</p> <p>Results</p> <p>The absolute subendocardial C-strain was significantly lower for stenotic (-7.5 ± 12.6%) than non-stenotic segment (-18.8 ± 10.2%, p < 0.0001). There was no difference in epicardial C-strain between the two groups. Use of cutoff thresholds for subendocardial C-strain differentiated stenotic segments from non-stenotic segments with a sensitivity of 77%, a specificity of 70%, and areas under the curve (AUC) of 0.76. The absolute subendocardial C-strain was significantly lower for ischemic (-6.7 ± 13.1%) than non-ischemic segments (-21.6 ± 7.0%, p < 0.0001). The absolute epicardial C-strain was also significantly lower for ischemic (-5.1 ± 7.8%) than non-ischemic segments (-9.6 ± 9.1%, p < 0.05). Use of cutoff thresholds for subendocardial C-strain differentiated ischemic segments from non-ischemic segments with sensitivities of 86%, specificities of 84%, and AUC of 0.86.</p> <p>Conclusions</p> <p>Analysis of tagged CMR can non-invasively demonstrate predominant impairment of subendocardial strain in the chronic ischemic myocardium at rest.</p

Myocardial tagging by Cardiovascular Magnetic Resonance: evolution of techniques--pulse sequences, analysis algorithms, and applications

Cardiovascular magnetic resonance (CMR) tagging has been established as an essential technique for measuring regional myocardial function. It allows quantification of local intramyocardial motion measures, e.g. strain and strain rate. The invention of CMR tagging came in the late eighties, where the technique allowed for the first time for visualizing transmural myocardial movement without having to implant physical markers. This new idea opened the door for a series of developments and improvements that continue up to the present time. Different tagging techniques are currently available that are more extensive, improved, and sophisticated than they were twenty years ago. Each of these techniques has different versions for improved resolution, signal-to-noise ratio (SNR), scan time, anatomical coverage, three-dimensional capability, and image quality. The tagging techniques covered in this article can be broadly divided into two main categories: 1) Basic techniques, which include magnetization saturation, spatial modulation of magnetization (SPAMM), delay alternating with nutations for tailored excitation (DANTE), and complementary SPAMM (CSPAMM); and 2) Advanced techniques, which include harmonic phase (HARP), displacement encoding with stimulated echoes (DENSE), and strain encoding (SENC). Although most of these techniques were developed by separate groups and evolved from different backgrounds, they are in fact closely related to each other, and they can be interpreted from more than one perspective. Some of these techniques even followed parallel paths of developments, as illustrated in the article. As each technique has its own advantages, some efforts have been made to combine different techniques together for improved image quality or composite information acquisition. In this review, different developments in pulse sequences and related image processing techniques are described along with the necessities that led to their invention, which makes this article easy to read and the covered techniques easy to follow. Major studies that applied CMR tagging for studying myocardial mechanics are also summarized. Finally, the current article includes a plethora of ideas and techniques with over 300 references that motivate the reader to think about the future of CMR tagging

Self-induced topological transition in phononic crystals by nonlinearity management

A new design paradigm of topology has recently emerged to manipulate the flow of phonons. At its heart lies a topological transition to a nontrivial state with exotic properties. This framework has been limited to linear lattice dynamics so far. Here we show a topological transition in a nonlinear regime and its implication in emerging nonlinear solutions. We employ nonlinearity management such that the system consists of masses connected with two types of nonlinear springs, "stiffening" and "softening" types, alternating along the length. We show, analytically and numerically, that the lattice makes a topological transition simply by changing the excitation amplitude and invoking nonlinear dynamics. Consequently, we witness the emergence of a new family of finite-frequency edge modes, not observed in linear phononic systems. We also report the existence of kink solitons at the topological transition point. These correspond to heteroclinic orbits that form a closed curve in the phase portrait separating the two topologically-distinct regimes. These findings suggest that nonlinearity can be used as a strategic tuning knob to alter topological characteristics of phononic crystals. These also provide fresh perspectives towards understanding a new family of nonlinear solutions in light of topology.Comment: 14 pages, 8 figure

A multiscale model for collagen alignment in wound healing

It is thought that collagen alignment plays a significant part in scar tissue formation during dermal wound healing. We present a multiscale model for collagen deposition and alignment during this process. We consider fibroblasts as discrete units moving within an extracellular matrix of collagen and fibrin modelled as continua. Our model includes flux induced alignment of collagen by fibroblasts, and contact guidance of fibroblasts by collagen fibres. We can use the model to predict the effects of certain manipulations, such as varying fibroblast speed, or placing an aligned piece of tissue in the wound. We also simulate experiments which alter the TGF-β concentrations in a healing dermal wound and use the model to offer an explanation of the observed influence of this growth factor on scarring

Cardiovascular magnetic resonance tagging imaging correlates with myocardial dysfunction and T2 mapping in idiopathic dilated cardiomyopathy

To evaluate the details of myocardial dysfunction in dilated cardiomyopathy (DCM) patients using tagging images and the correlation of tagging imaging with tissue characteristics. C

Cancer modelling: Getting to the heart of the problem

Paradoxically, improvements in healthcare that have enhanced the life expectancy of humans in the Western world have, indirectly, increased the prevalence of certain types of cancer such as prostate and breast. It remains unclear whether this phenomenon should be attributed to the ageing process itself or the cumulative effect of prolonged exposure to harmful environmental stimuli such as ultraviolet light, radiation and carcinogens (Franks and Teich, 1988). Equally, there is also compelling evidence that certain genetic abnormalities can predispose individuals to specific cancers (Ilyas et al., 1999). The variety of factors that have been implicated in the development of solid tumours stems, to a large extent, from the fact that ‘cancer’ is a generic term, often used to characterize a series of disorders that share common features. At this generic level of description, cancer may be viewed as a cellular disease in which controls that usually regulate growth and maintain homeostasis are disrupted. Cancer is typically initiated by genetic mutations that lead to enhanced mitosis of a cell lineage and the formation of an avascular tumour. Since it receives nutrients by diffusion from the surrounding tissue, the size of an avascular tumour is limited to several millimeters in diameter. Further growth relies on the tumour acquiring the ability to stimulate the ingrowth of a new, circulating blood supply from the host vasculature via a process termed angiogenesis (Folkman, 1974). Once vascularised, the tumour has access to a vast nutrient source and rapid growth ensues. Further, tumour fragments that break away from the primary tumour, on entering the vasculature, may be transported to other organs in which they may establish secondary tumours or metastases that further compromise the host. Invasion is another key feature of solid tumours whereby contact with the tissue stimulates the production of enzymes that digest the tissue, liberating space into which the tumour cells migrate. Thus, cancer is a complex, multiscale process. The spatial scales of interest range from the subcellular level, to the cellular and macroscopic (or tissue) levels while the timescales may vary from seconds (or less) for signal transduction pathways to months for tumour doubling times The variety of phenomena involved, the range of spatial and temporal scales over which they act and the complex way in which they are inter-related mean that the development of realistic theoretical models of solid tumour growth is extremely challenging. While there is now a large literature focused on modelling solid tumour growth (for a review, see, for example, Preziosi, 2003), existing models typically focus on a single spatial scale and, as a result, are unable to address the fundamental problem of how phenomena at different scales are coupled or to combine, in a systematic manner, data from the various scales. In this article, a theoretical framework will be presented that is capable of integrating a hierarchy of processes occurring at different scales into a detailed model of solid tumour growth (Alarcon et al., 2004). The model is formulated as a hybrid cellular automaton and contains interlinked elements that describe processes at each spatial scale: progress through the cell cycle and the production of proteins that stimulate angiogenesis are accounted for at the subcellular level; cell-cell interactions are treated at the cellular level; and, at the tissue scale, attention focuses on the vascular network whose structure adapts in response to blood flow and angiogenic factors produced at the subcellular level. Further coupling between the different spatial scales arises from the transport of blood-borne oxygen into the tissue and its uptake at the cellular level. Model simulations will be presented to illustrate the effect that spatial heterogeneity induced by blood flow through the vascular network has on the tumour’s growth dynamics and explain how the model may be used to compare the efficacy of different anti-cancer treatment protocols

On motion in dynamic magnetic resonance imaging: Applications in cardiac function and abdominal diffusion

La imagen por resonancia magnética (MRI), hoy en día, representa una potente herramienta para el diagnóstico clínico debido a su flexibilidad y sensibilidad a un amplio rango de propiedades del tejido. Sus principales ventajas son su sobresaliente versatilidad y su capacidad para proporcionar alto contraste entre tejidos blandos. Gracias a esa versatilidad, la MRI se puede emplear para observar diferentes fenómenos físicos dentro del cuerpo humano combinando distintos tipos de pulsos dentro de la secuencia. Esto ha permitido crear distintas modalidades con múltiples aplicaciones tanto biológicas como clínicas. La adquisición de MR es, sin embargo, un proceso lento, lo que conlleva una solución de compromiso entre resolución y tiempo de adquisición (Lima da Cruz, 2016; Royuela-del Val, 2017). Debido a esto, la presencia de movimiento fisiológico durante la adquisición puede conllevar una grave degradación de la calidad de imagen, así como un incremento del tiempo de adquisición, aumentando así tambien la incomodidad del paciente. Esta limitación práctica representa un gran obstáculo para la viabilidad clínica de la MRI. En esta Tesis Doctoral se abordan dos problemas de interés en el campo de la MRI en los que el movimiento fisiológico tiene un papel protagonista. Éstos son, por un lado, la estimación robusta de parámetros de rotación y esfuerzo miocárdico a partir de imágenes de MR-Tagging dinámica para el diagnóstico y clasificación de cardiomiopatías y, por otro, la reconstrucción de mapas del coeficiente de difusión aparente (ADC) a alta resolución y con alta relación señal a ruido (SNR) a partir de adquisiciones de imagen ponderada en difusión (DWI) multiparamétrica en el hígado.Departamento de Teoría de la Señal y Comunicaciones e Ingeniería TelemáticaDoctorado en Tecnologías de la Información y las Telecomunicacione