The vertebral venous system in healthy and scoliotic adolescent spines - a 3D MRI investigation

Introduction. The venous drainage system within vertebral bodies (VBs) has been well documented previously in cadaveric specimens. Advances in 3D imaging and image processing now allow for in vivo quantification of larger venous vessels, such as the basivertebral vein. Differences between healthy and scoliotic VB veins can therefore be investigated. Methods. 20 healthy adolescent controls and 21 AIS patients were recruited (with ethics approval) to undergo 3D MRI, using a 3 Tesla, T1-weighted 3D gradient echo sequence, resulting in 512 slices across the thoraco-lumbar spine, with a voxel size of 0.5x0.5x0.5mm. Using Amira Filament Editor, five transverse slices through the VB were examined simultaneously and the resulting observable vascular network traced. Each VB was assessed, and a vascular network recorded when observable. A local coordinate system was created in the centre of each VB and the vascular networks aligned to this. The length of the vascular network on the left and right sides (with a small central region) of the VB was calculated, and the spatial patterning of the networks assessed level-by-level within each subject. Results. An average of 6 (range 4-10) vascular networks, consistent with descriptions of the basivertebral vein, were identifiable within each subject, most commonly between T10-L1. Differences were seen in the left/right distribution of vessels in the control and AIS subjects. Healthy controls saw a percentage distribution of 29:18:53 across the left:centre:right regions respectively, whereas the AIS subjects had a slightly shifted distribution of 33:25:42. The control group showed consistent spatial patterning of the vascular networks across most levels, but this was not seen in the AIS group. Conclusion. Observation and quantification of the basivertebral vein in vivo is possible using 3D MRI. The AIS group lacked the spatial pattern repetition seen in the control group and minor differences were seen in the left/right distribution of vessels.

Sequential MRI of growing scoliotic spines reveals individual level deformities that are clinically masked by the cobb angle

INTRODUCTION. Clinically, the Cobb angle method measures the overall scoliotic curve in the coronal plane but does not measure individual vertebra and disc wedging. The contributions of the vertebrae and discs in the growing scoliotic spine were measured to investigate coronal plane deformity progression with growth. METHODS. A 0.49mm isotropic 3D MRI technique was developed to investigate the level-by-level changes that occur in the growing spine of a group of Adolescent Idiopathic Scoliosis (AIS) patients, who received two to four sequential scans (spaced 3-12 months apart). The coronal plane wedge angles of each vertebra and disc in the major curve were measured to capture any changes that occurred during their adolescent growth phase. RESULTS. Seventeen patients had at least two scans. Mean patient age was 12.9 years (SD 1.5 years). Sixteen were classified as right-sided major thoracic Lenke Type 1 (one left sided). Mean standing Cobb angle at initial presentation was 31° (SD 12°). Six received two scans, nine three scans and two four scans, with 65% showing a Cobb angle progression of 5° or more between scans. Overall, there was no clear pattern of deformity progression of individual vertebrae and discs, nor between patients who progressed and those who didn’t. There were measurable changes in the wedging of the vertebrae and discs in all patients. In sequential scans, change in direction of wedging was also seen. In several patients there was reverse wedging in the discs that counteracted increased wedging of the vertebrae such that no change in overall Cobb angle was seen. CONCLUSION. Sequential MRI data showed complex patterns of deformity progression. Changes to the wedging of individual vertebrae and discs may occur in patients who have no increase in Cobb angle measure; the Cobb method alone may be insufficient to capture the complex mechanisms of deformity progression.

Sequential MRI reveals individual level deformities in the growing scoliotic spine that are not seen clinically by the cobb angle

Clinically, the Cobb angle method measures the overall scoliotic curve in the coronal plane but does not measure individual vertebra and disc wedging. The contributions of the vertebrae and discs in the growing scoliotic spine were measured to investigate coronal plane deformity progression with growth. Sequential MRI data in this project showed complex patterns of deformity progression. Changes to the wedging of individual vertebrae and discs may occur in patients who have no increase in Cobb angle measure; the Cobb method alone may be insufficient to capture the complex mechanisms of deformity progression.

Imaging tumour distribution of a polymeric drug delivery platformin vivoby PET-MRI

Background Over the past decade, molecular imaging has played a key role in the progression of drug delivery platforms from concept to commercialisation. Of the molecular imaging techniques commonly utilised, positron emission tomography (PET) can yield a breadth of information not easily accessible by other methodologies and when combined with other complementary imaging modalities, is a powerful tool for pre- and clinical development of therapeutics. However, very little research has focussed on the information available from complimentary imaging modalities. This paper reports on the data-rich methodologies of contrast enhanced PET/CT and PET/MRI for probing efficacy of polymer drug delivery platforms. Results The information available from an ExiTron nano 6000 contrast enhanced PET/CT and a gadolinium (Gd) enhanced PET/MRI image of a 64Cu labeled HBP in the same mouse was qualitatively compared. Conclusions Gd contrast enhanced PET/MRI offers a powerful methodology for investigating the distribution of polymer drug delivery platforms in vivo and throughout a tumour volume. Furthermore, information about depth of penetration away from primary blood vessels can be gleaned, potentially leading to development of more efficacious delivery vehicles for clinical use.

MRI-based strain and strain rate analysis of left ventricle: A modified hierarchical transformation model

Background Different from other indicators of cardiac function, such as ejection fraction and transmitral early diastolic velocity, myocardial strain is promising to capture subtle alterations that result from early diseases of the myocardium. In order to extract the left ventricle (LV) myocardial strain and strain rate from cardiac cine-MRI, a modified hierarchical transformation model was proposed. Methods A hierarchical transformation model including the global and local LV deformations was employed to analyze the strain and strain rate of the left ventricle by cine-MRI image registration. The endocardial and epicardial contour information was introduced to enhance the registration accuracy by combining the original hierarchical algorithm with an Iterative Closest Points using Invariant Features algorithm. The hierarchical model was validated by a normal volunteer first and then applied to two clinical cases (i.e., the normal volunteer and a diabetic patient) to evaluate their respective function. Results Based on the two clinical cases, by comparing the displacement fields of two selected landmarks in the normal volunteer, the proposed method showed a better performance than the original or unmodified model. Meanwhile, the comparison of the radial strain between the volunteer and patient demonstrated their apparent functional difference. Conclusions The present method could be used to estimate the LV myocardial strain and strain rate during a cardiac cycle and thus to quantify the analysis of the LV motion function.

A follow up MRI-based geometry and computational fluid dynamics study of carotid bifurcation

Cardiovascular disease is the leading causes of death in the developed world. Wall shear stress (WSS) is associated with the initiation and progression of atherogenesis. This study combined the recent advances in MR imaging and computational fluid dynamics (CFD) and evaluated the patient-specific carotid bifurcation. The patient was followed up for 3 years. The geometry changes (tortuosity, curvature, ICA/CCA area ratios, central to the cross-sectional curvature, maximum stenosis) and the CFD factors (Velocity distribute, Wall Shear Stress (WSS) and Oscillatory Shear Index (OSI)) were compared at different time points.The carotid stenosis was a slight increase in the central to the cross-sectional curvature, and it was minor and variable curvature changes for carotid centerline. The OSI distribution presents ahigh-values in the same region where carotid stenosis and normal border, indicating complex flow and recirculation.The significant geometric changes observed during the follow-up may also cause significant changes in bifurcation hemodynamics.

A new method for the in vivo identification of mechanical properties in arteries from cine MRI images: Theoretical framework and validation

Quantifying the stiffness properties of soft tissues is essential for the diagnosis of many cardiovascular diseases such as atherosclerosis. In these pathologies it is widely agreed that the arterial wall stiffness is an indicator of vulnerability. The present paper focuses on the carotid artery and proposes a new inversion methodology for deriving the stiffness properties of the wall from cine-MRI (magnetic resonance imaging) data. We address this problem by setting-up a cost function defined as the distance between the modeled pixel signals and the measured ones. Minimizing this cost function yields the unknown stiffness properties of both the arterial wall and the surrounding tissues. The sensitivity of the identified properties to various sources of uncertainty is studied. Validation of the method is performed on a rubber phantom. The elastic modulus identified using the developed methodology lies within a mean error of 9.6%. It is then applied to two young healthy subjects as a proof of practical feasibility, with identified values of 625 kPa and 587 kPa for one of the carotid of each subject.

MRI-based fatigue analysis predicts plaque vulnerability: A study comparing symptomatic and asymptomatic individuals

BACKGROUND: Rupture of atheromatous plaque in the carotid artery often leads to thrombosis and subsequent stroke. The mechanism of plaque rupture is not entirely clear but is thought to be a multi-factorial process involving thinning and weakening of the fibrous cap and biomechanical stress as the trigger leading to plaque rupture. As the cardiovascular system is a classic fatigue environment, the weakening of plaque leading to rupture may be a fatigue process, which is a symptomatically quiescent but potentially progressive failure process. In this study, we used a fatigue analysis based on in vivo magnetic resonance imaging (MRI) to investigate the rupture initiation location, crack propagation path and fatigue life within plaques of asymptomatic and symptomatic individuals. METHODS: Forty non-consecutive subjects (20 symptomatic and 20 asymptomatic) underwent high-resolution multi-sequence in vivo MRI of the carotid bifurcation. Fatigue analysis was performed based on the plaque geometry derived from in vivo MRI of the carotid artery at the point of maximum stenosis. Paris’ Law in fracture mechanics is adopted to determine the fatigue crack growth rate. Incremental crack propagation was dynamically simulated based on stress distributions. Plaque initiation location, crack propagation path and fatigue cycle of symptomatic and asymptomatic individuals were compared. RESULTS: Cracks were often found to begin at the lumen wall at areas of stress concentration. The preferred rupture direction was radial from the lumen center. The crack initially advanced slowly but accelerated as it developed, depending on plaque morphology. The fatigue cycles of symptomatic plaques were significantly less than those in the asymptomatic group (2.3 ± 0.9 vs 3.1 ± 0.7 (x106); p = 0.003). CONCLUSIONS: The number of cycles to rupture in symptomatic patients was higher than those predicted in asymptomatic patients by fatigue analysis, suggesting the possibility that plaques with a less fatigue life may be more prone to be symptomatic and rupture. If further validated by large-scale longitudinal studies, fatigue analysis based on high resolution in vivo MRI could potentially act as a useful tool for risk assessment of carotid atheroma.

In vivo velocity vector imaging and time-resolved strain rate measurements in the wall of blood vessels using MRI

In this paper, we present a new approach for velocity vector imaging and time-resolved measurements of strain rates in the wall of human arteries using MRI and we prove its feasibility on two examples: in vitro on a phantom and in vivo on the carotid artery of a human subject. Results point out the promising potential of this approach for investigating the mechanics of arterial tissues in vivo.

Normalized wall index specific and MRI-based stress analysis of atherosclerotic carotid plaques: A study comparing acutely symptomatic and asymptomatic patients

Background: Biomechanical stresses play an important role in determining plaque stability. Quantification of these simulated stresses can be potentially used to assess plaque vulnerability and differentiate different patient groups. Methods and Results: 54 asymptomatic and 45 acutely symptomatic patients underwent in vivo multicontrast magnetic resonance imaging (MRI) of the carotid arteries. Plaque geometry used for finite element analysis was derived from in vivo MRI at the sites of maximum and minimum plaque burden. In total, 198 slices were used for the computational simulations. A pre-shrink technique was used to refine the simulation. Maximum principle stress at the vulnerable plaque sites (ie, critical stress) was extracted for the selected slices and a comparison was performed between the 2 groups. Critical stress in the slice with maximum plaque burden is significantly higher in acutely symptomatic patients as compared to asymptomatic patients (median, inter quartile range: 198.0 kPa (119.8-359.0 kPa) vs 138.4 kPa (83.8-242.6 kPa), P=0.04). No significant difference was found in the slice with minimum plaque burden between the 2 groups (196.7 kPa (133.3-282.7 kPa) vs 182.4 kPa (117.2-310.6 kPa), P=0.82). Conclusions: Acutely symptomatic carotid plaques have significantly high biomechanical stresses than asymptomatic plaques. This might be potentially useful for establishing a biomechanical risk stratification criteria based on plaque burden in future studies.

Arterial luminal curvature and fibrous-cap thickness affect critical stress conditions within atherosclerotic plaque: An in vivo MRI-based 2D finite-element study

High mechanical stress in atherosclerotic plaques at vulnerable sites, called critical stress, contributes to plaque rupture. The site of minimum fibrous cap (FC) thickness (FCMIN) and plaque shoulder are well-documented vulnerable sites. The inherent weakness of the FC material at the thinnest point increases the stress, making it vulnerable, and it is the big curvature of the lumen contour over FC which may result in increased plaque stress. We aimed to assess critical stresses at FCMIN and the maximum lumen curvature over FC (LCMAX) and quantify the difference to see which vulnerable site had the highest critical stress and was, therefore, at highest risk of rupture. One hundred patients underwent high resolution carotid magnetic resonance (MR) imaging. We used 352 MR slices with delineated atherosclerotic components for the simulation study. Stresses at all the integral nodes along the lumen surface were calculated using the finite-element method. FCMIN and LCMAX were identified, and critical stresses at these sites were assessed and compared. Critical stress at FC MIN was significantly lower than that at LCMAX (median: 121.55 kPa; inter quartile range (IQR) = [60.70-180.32] kPa vs. 150.80 kPa; IQR = [91.39-235.75] kPa, p < 0.0001). If critical stress at FCMIN was only used, then the stress condition of 238 of 352 MR slices would be underestimated, while if the critical stress at LCMAX only was used, then 112 out of 352 would be underestimated. Stress analysis at FCMIN and LCMAX should be used for a refined mechanical risk assessment of atherosclerotic plaques, since material failure at either site may result in rupture.

Stress analysis of carotid plaque based on in vivo MRI of acute symptomatic and asymptomatic patients

Stress analysis within carotid plaques based on in vivo MR imaging has shown to be useful for the identification of vulnerable atheroma. This study is to investigate whether magnetic resonance imaging (MRI) based-biomechanical stress analysis of carotid plaques can differentiate acute symptomatic and asymptomatic patients. 54 asymptomatic and 45 acute symptomatic patients underwent in vivo multi-contrast MRI of the carotid arteries. Plaque geometry used for finite element analysis was derived from in vivo MR images at the site of maximum and minimum plaque burden. In total 198 slices were used for the computational simulations. A pre shrink technique was used to refine the simulation. Maximum principle stress at the vulnerable plaque sites (i.e. critical stress) was extracted for the selected slices and a comparison was performed between the two groups. Critical stress at the site of maximum plaque burden is significantly higher in acute symptomatic patients as compared to asymptomatic patients [median: 198.0kPa (inter quartile range (IQR) = (119.8 - 359.0) vs. 138.4kPa (83.8, 242.6), p=0.04]. No significant difference was found at the minimum plaque burden site between the two groups [196.7kPa (133.3- 282.7) vs. 182.4kPa (117.2 - 310. 6), p=0.82). Stress analysis at the site of maximal plaque burden can be effectively used for differentiating acute symptomatic carotid plaques from asymptomatic plaques. This maybe potentially used for development of biomechanical risk stratification criteria based on plaque burden in future studies.

In vivo mri-based estimation of time-dependent elastic modulus in healthy arteries

The elastic properties of the arterial wall have been the subject of physiological, clinical and biomedical research for many years. There is convincing evidence that the elastic properties of the large arteries are seriously impaired in the presence of cardiovascular disease (CVD), due to alterations in the intrinsic structural and functional characteristics of vessels [1]. Early detection of changes in the elastic modulus of arteries would provide a powerful tool for both monitoring patients at high cardiovascular risk and testing the effects of pharmaceuticals aimed at stabilizing existing plaques by stiffening them or lowering the lipids.

Temporal dependence of in vivo USPIO-enhanced MRI signal changes in human carotid atheromatous plaques

Introduction: Ultrasmall superparamagnetic iron oxide (USPIO)-enhanced MRI has been shown to be a useful modality to image activated macrophages in vivo, which are principally responsible for plaque inflammation. This study determined the optimum imaging time-window to detect maximal signal change post-USPIO infusion using T1-weighted (T1w), T2*- weighted (T2*w) and quantitative T2*(qT 2*) imaging. Methods: Six patients with an asymptomatic carotid stenosis underwent high resolution T1w, T2*w and qT2*MR imaging of their carotid arteries at 1.5 T. Imaging was performed before and at 24, 36, 48, 72 and 96 h after USPIO (Sinerem™, Guerbet, France) infusion. Each slice showing atherosclerotic plaque was manually segmented into quadrants and signal changes in each quadrant were fitted to an exponential power function to model the optimum time for post-infusion imaging. Results: The power function determining the mean time to convergence for all patients was 46, 41 and 39 h for the T1w, T 2*w and qT2*sequences, respectively. When modelling each patient individually, 90% of the maximum signal intensity change was observed at 36 h for three, four and six patients on T1w, T 2*w and qT2*, respectively. The rates of signal change decrease after this period but signal change was still evident up to 96 h. Conclusion: This study showed that a suitable imaging window for T 1w, T2*w and qT2*signal changes post-USPIO infusion was between 36 and 48 h. Logistically, this would be convenient in bringing patients back for one post-contrast MRI, but validation is required in a larger cohort of patients.

MRI-based measurement of carotid mechanical properties in asymptomatic patients

Arterial mechanical property may be a potential variable for risk stratification. Large artery and central arterial compliance have been shown not only to correlate well with overall cardiovascular outcome in large epidemiological studies [1, 2] but also to correlate with coronary atherosclerotic burden as measured by conventional angiography [3]. Until recently, real-time B-mode ultrasound combined with simultaneous blood pressure measurements have been used to assess large artery compliance [4]. These techniques have an excellent temporal resolution but are unable to provide adequate spatial resolution to determine changes in vessel area as opposed to diameter and make the assumption that the vessel is perfectly round. Attempts to use MR imaging to measure large artery compliance have been published previously [5]. However, they have not utilised simultaneous blood pressure measurements during sequence acquisition. We report a technique using regular and simultaneous blood pressure measurement during 2 dimensional phase contrast magnetic resonance imaging 2DPC-MRI to determine local carotid compliance.