A left ventricle model to predict post-revascularization ejection fraction based on cine magnetic resonance images

The left ventricular ejection fraction (LVEF) is an important clinical indicator of the cardiac function and long-term outcome for patients with coronary artery disease. A biomechanical model of the left ventricle was developed to quantitatively predict post-revascularization LVEF based on noninvasive magnetic resonance imaging. The myocardium was categorized into normal, hibernating, and infarcted regions from the ventricular short-axis images. Assuming that hibernating tissue would potentially regain contractility after revascularization, the expected maximum post-revascularization LVEF was calculated for four patients with chronic left ventricular dysfunction. The predictions were within three ejection fraction points of the follow-up LVEFs. This model may be useful to estimate the outcome and efficacy of revascularization plans.     

Changes in magnetic resonance images in human skeletal muscle after eccentric exercise

To investigate the time-course of changes in transverse relaxation time (T₂) and cross-sectional area (CSA) of the quadriceps muscle after a single session of eccentric exercise, magnetic resonance imaging was performed on six healthy male volunteers before and at 0, 7, 15, 20, 30 and 60 min and 12, 24, 36, 48, 72 and 168 h after exercise. Although there was almost no muscle soreness immediately after exercise, it started to increase 1 day after, peaking 1–2 days after the exercise (P<0.01). Immediately after exercise, T₂ increased significantly in the rectus femoris, vastus lateralis and intermedius muscles (P<0.05) and decreased quickly continuing until 60 min after exercise. At and after the 12th h, a significant increase was perceived again in the T₂ values of the vastus lateralis and intermedius muscles (P<0.01) [maximum 9.3 (SEM 2.8)% and 10.9 (SEM 2.2)%, respectively]. The maximal values were exhibited at 24–36 h after exercise. In contrast, the rectus femoris muscle showed no delayed-stage increase. Also, in CSA, an increase after 12 h was observed in addition to the one immediately after exercise in the vastus lateralis, intermedius and medialis and quadriceps muscles as a whole (P < 0.01), reaching the maximal values at 12–24 h after exercise. The plasma creative kinase activity remained unchanged up to 24 h after and then increased significantly 48 h after exercise (P < 0.05). Beginning 12 h after exercise, the subjects whose T₂ and CSA increased less than the others displayed a faster decrease in muscle soreness. These results suggested that T₂ and CSA displayed bimodal responses after eccentric exercise and the time-courses of changes in them were similar to those in muscle soreness.     

Measurement of strain in physical models of brain injury: a method based on HARP analysis of tagged magnetic resonance images (MRI)

Two-dimensional (2-D) strain fields were estimated non-invasively in two simple experimental models of closed-head brain injury. In the first experimental model, shear deformation of a gel was induced by angular acceleration of its spherical container In the second model the brain of a euthanized rat pup was deformed by indentation of its skull. Tagged magnetic resonance images (MRI) were obtained by gated image acquisition during repeated motion. Harmonic phase (HARP) images corresponding to the spectral peaks of the original tagged MRI were obtained, following procedures proposed by Osman, McVeigh and Prince. Two methods of HARP strain analysis were applied, one based on the displacement of tag line intersections, and the other based on the gradient of harmonic phase. Strain analysis procedures were also validated on simulated images of deformed grids. Results show that it is possible to visualize deformation and to quantify strain efficiently in animal models of closed head injury.     

An experimental analysis of the flow field in a three-dimensional model of the human carotid artery bifurcation

Steady flow measurements were carried out in a rigid three-dimensional model of the human carotid artery bifurcation at a Reynolds number of 640 and a flow division ratio of 50/50. Both axial and secondary velocities were measured with a laser-Doppler anemometer. In the bulb opposite to the flow divider a zone with negative axial velocities was found with a maximal diameter of about 60% of the local diameter of the branch and a cross-sectional extent of about 25% of the local cross-sectional area. In the bulb the maximum axial velocity shifted towards the divider wall and at the end of the bulb an axial velocity plateau arose near the non-divider wall. Halfway through the bulb, secondary flow showed a vortex through which fluid flowed towards the divider wall near the bifurcation plane and back towards the non-divider wall near the upper walls.     

Three-dimensional geometry of the human carotid artery

Accurate characterization of carotid artery geometry is vital to our understanding of the pathogenesis of atherosclerosis. Three-dimensional computer reconstructions based on medical imaging are now ubiquitous; however, mean carotid artery geometry has not yet been comprehensively characterized. The goal of this work was to build and study such geometry based on data from 16 male patients with severe carotid artery disease. Results of computerized tomography angiography were used to analyze the cross-sectional images implementing a semiautomated segmentation algorithm. Extracted data were used to reconstruct the mean three-dimensional geometry and to determine average values and variability of bifurcation and planarity angles, diameters and cross-sectional areas. Contrary to simplified carotid geometry typically depicted and used, our mean artery was tortuous exhibiting nonplanarity and complex curvature and torsion variations. The bifurcation angle was 36 deg?±?11 deg if measured using arterial centerlines and 15 deg?±?14 deg if measured between the walls of the carotid bifurcation branches. The average planarity angle was 11 deg?±?10 deg. Both bifurcation and planarity angles were substantially smaller than values reported in most studies. Cross sections were elliptical, with an average ratio of semimajor to semiminor axes of 1.2. The cross-sectional area increased twofold in the bulb compared to the proximal common, but then decreased 1.5-fold for the combined area of distal internal and external carotid artery. Inter-patient variability was substantial, especially in the bulb region; however, some common geometrical features were observed in most patients. Obtained quantitative data on the mean carotid artery geometry and its variability among patients with severe carotid artery disease can be used by biomedical engineers and biomechanics vascular modelers in their studies of carotid pathophysiology, and by endovascular device and materials manufacturers interested in the mean geometrical features of the artery to target the broad patient population.     

In vivo detection of c-MET expression in a rat hepatocarcinogenesis model using molecularly targeted magnetic resonance imaging

The multifunctional growth factor scatter factor/hepatocyte growth factor and its tyrosine kinase receptor, c-MET, have been implicated in the genesis and malignant progression of numerous human malignancies, including hepatocellular carcinomas. The incidence of hepatocellular carcinomas in the United States has increased noticeably over the past two decades and is listed as the fifth major cancer in men worldwide. In this study, we used a choline-deficient l-amino acid (CDAA)-defined rat hepatocarcinogenesis model to visualize increased in vivo expression of the c-MET antigen in neoplastic lesion formation with the use of a super paramagnetic iron oxide (SPIO)-anti-c-MET molecularly targeted magnetic resonance imaging (MRI) contrast agent. SPIO-anti-c-MET was used for the first time to detect overexpression of c-MET in neoplastic nodules and tumors within the livers of CDAA-treated rats, as determined by a decrease in MRI signal intensity and a decrease in regional T(2) values. Specificity for the binding of the molecularly targeted anti-c-MET contrast agent was determined using rat hepatoma (H4-II-E-C3) cell cultures and immunofluorescence microscopic imaging of the targeting agents within neoplastic liver tissue 1 to 2 hours following intravenous administration of SPIO-anti-c-MET and MRI investigation. This method has the ability to visualize in vivo the overexpression of c-MET at early developmental stages of tumor formation.     

In vivo tracking of the human patella using cine phase contrast magnetic resonance imaging

Improper patellar tracking is often considered to be the cause of patellar-femoral pain. Unfortunately, our knowledge of patellar-femoral-tibial (knee) joint kinematics is severely limited due to a lack of three-dimensional, noninvasive, in vivo measurement techniques. This study presents the first large-scale, dynamic, three-dimensional, noninvasive, in vivo study of nonimpaired knee joint kinematics during volitional leg extensions. Cine-phase contrast magnetic resonance imaging was used to measure the velocity profiles of the patella, femur, and tibia in 18 unimpaired knees during leg extensions, resisted by a 34 N weight. Bone displacements were calculated through integration and then converted into three-dimensional orientation angles. We found that the patella displaced laterally, superiorly, and anteriorly as the knee extended. Further, patellar flexion lagged knee flexion, patellar tilt was variable, and patellar rotation was fairly constant throughout extension.     

A numerical and experimental analysis of the flow field in a two-dimensional model of the human carotid artery bifurcation

Axial velocities were measured in an enlarged, two-dimensional, rigid model of the carotid artery bifurcation by means of a laser-Doppler anemometer, under both steady and unsteady flow conditions. Also a numerical model was developed, based on the finite element approximation of the Navier-Stokes and continuity equations. From this study it appeared that the numerically predicted velocities agree well with the experimentally obtained values. Besides, the bifurcation hardly influenced the upstream flow in the main branch (common carotid artery), high velocity gradients were observed at the divider walls of the daughter branches (internal and external carotid arteries) and large zones with reversed flow were present near the nondivider walls of these branches. For steady flow the maximal diameter of this zone at the entrance of the internal carotid artery (carotid sinus) was about 25% of the local diameter of this branch. For unsteady flow this zone was absent during the initial phase of flow acceleration and maximal at the end of flow deceleration with a maximal diameter of about 50% of the local diameter of the carotid sinus.     

Data distributions in magnetic resonance images: A review

Many image processing methods applied to magnetic resonance (MR) images directly or indirectly rely on prior knowledge of the statistical data distribution that characterizes the MR data. Also, data distributions are key in many parameter estimation problems and strongly relate to the accuracy and precision with which parameters can be estimated. This review paper provides an overview of the various distributions that occur when dealing with MR data, considering both single-coil and multiple-coil acquisition systems. The paper also summarizes how knowledge of the MR data distributions can be used to construct optimal parameter estimators and answers the question as to what precision may be achieved ultimately from a particular MR image.     

Comparison of model and nuclear magnetic resonance structures for the human inflammatory protein C5a

The model structure previously proposed for human C5a, based upon the crystal structure of the homologous protein human C3a, is compared to the solution structure of human C5a recently determined by nuclear magnetic resonance (NMR) methods in our laboratory. The general folding and helix topography of the C5a protein were modeled very well. The N-terminus, which is disordered in the C3a crystal, was correctly predicted in the C5a model both as to its being a helix and as to its docking site on the rest of the molecule. On the other hand, the NMR data show that the biologically important C-terminal residues are disordered in solution, unlike the model and the C3a crystal structure where this region was helical.     

Variability of a three-dimensional finite element model constructed using magnetic resonance images of a knee for joint contact stress analysis

Magnetic resonance (MR) imaging has been widely used to evaluate the thickness and volume of articular cartilage both in vivo and in vitro. While morphological information on the cartilage can be obtained using MR images, image processing for extracting geometric boundaries of the cartilage may introduce variations in the thickness of the cartilage. To evaluate the variability of using MR images to construct finite element (FE) knee cartilage models, five investigators independently digitized the same set of MR images of a human knee. The topology of cartilage thickness was determined using a minimal distance algorithm. Less than 8 percent variation in cartilage thickness was observed from the digitized data. The effect of changes in cartilage thickness on contact stress analysis was then investigated using five FE models of the knee. One FE model (average FE model) was constructed using the mean values of the digitized contours of the cartilage, and the other four were constructed by varying the thickness of the average FE model by +/- 5 percent and +/- 10 percent, respectively. The results demonstrated that under axial tibial compressive loading (up to 1,400 N), variations of cartilage thickness caused by digitization of MR images may result in a difference of approximately 10 percent in peak contact stresses (surface pressure, von Mises stress, and hydrostatic pressure) in the cartilage. A reduction of cartilage thickness caused increases of contact stresses, while an increase of cartilage thickness reduced contact stresses. Furthermore, the effect of variation of material properties of the cartilage on contact stress analysis was investigated. The peak contact stress increased almost linearly with the Young's modulus of the cartilage. The peak von Mises stress was dramatically reduced when the Poisson,s ratio was increased from 0.05 to 0.49 under an axial compressive load of 1,400 N, while peak hydrostatic pressure was dramatically increased. Peak surface pressure was also increased with the Poisson's ratio, but with a lower magnitude compared to von Mises stress and hydrostatic pressure. In conclusion, the imaging process may cause 10 percent variations in peak contact stress, and the predicted stress distribution is sensitive to the accuracy of the material properties of the cartilage model, especially to the variation of Poisson's ratio.     

Automated analysis of craniofacial morphology using magnetic resonance images

Quantitative analysis of craniofacial morphology is of interest to scholars working in a wide variety of disciplines, such as anthropology, developmental biology, and medicine. T1-weighted (anatomical) magnetic resonance images (MRI) provide excellent contrast between soft tissues. Given its three-dimensional nature, MRI represents an ideal imaging modality for the analysis of craniofacial structure in living individuals. Here we describe how T1-weighted MR images, acquired to examine brain anatomy, can also be used to analyze facial features. Using a sample of typically developing adolescents from the Saguenay Youth Study (N = 597; 292 male, 305 female, ages: 12 to 18 years), we quantified inter-individual variations in craniofacial structure in two ways. First, we adapted existing nonlinear registration-based morphological techniques to generate iteratively a group-wise population average of craniofacial features. The nonlinear transformations were used to map the craniofacial structure of each individual to the population average. Using voxel-wise measures of expansion and contraction, we then examined the effects of sex and age on inter-individual variations in facial features. Second, we employed a landmark-based approach to quantify variations in face surfaces. This approach involves: (a) placing 56 landmarks (forehead, nose, lips, jaw-line, cheekbones, and eyes) on a surface representation of the MRI-based group average; (b) warping the landmarks to the individual faces using the inverse nonlinear transformation estimated for each person; and (3) using a principal components analysis (PCA) of the warped landmarks to identify facial features (i.e. clusters of landmarks) that vary in our sample in a correlated fashion. As with the voxel-wise analysis of the deformation fields, we examined the effects of sex and age on the PCA-derived spatial relationships between facial features. Both methods demonstrated significant sexual dimorphism in craniofacial structure in areas such as the chin, mandible, lips, and nose.     

39K nuclear magnetic resonance and a mathematical model of K+ transport in human erythrocytes

³⁹K nuclear magnetic resonance was used to measure the efflux of K⁺ from suspensions of human erythrocytes [red blood cells (RBCs)], that occurred in response to the calcium ionophore, A23187 and calcium ions; the latter activate the Gárdos channel. Signals from the intra- and extracellular populations of ³⁹K⁺ were selected on the basis of their longitudinal relaxation times, T ₁, by using an inversion- recovery pulse sequence with the mixing time, τ₁, chosen to null one or other of the signals. Changes in RBC volume consequent upon efflux of the ions also changed the T ₁ values so a new theory was implemented to obviate a potential artefact in the data analysis. The velocity of the K⁺ efflux mediated by the Gárdos channel was 1.19±0.40 mmol (L RBC)⁻¹ min⁻¹ at 37°C.     

In vivo magnetic resonance imaging of the reproductive organs in a passerine bird species

For researchers working on the reproductive anatomy of birds, the standard tool has been the scalpel. We demonstrate that magnetic resonance imaging (MRI) allows us to study non-invasively the testes and ovaries of small birds. Thus, repeated measurements of the same individual can be obtained at different life-history stages or before and after experimental treatment with a minimum of harm done to the birds. We present cross-sectional images of male and female Garden Warblers Sylvia borin at different stages of their gonadal maturation cycle to illustrate the high resolution that can be achieved through MRI. Volumetric analyses enabled us to determine the total testicular volume with high accuracy.     

In vivo magnetic resonance spectroscopy in the evaluation of mitochondrial disorders

Magnetic resonance spectroscopy (MRS) has been utilized to study several metabolic pathways in vivo and in live tissues in vitro non-invasively. Despite its inherited lack of sensitivity, its application has extended all the way to in situ human tissues and organs since proper technical advancements were devised. Examples of its application described here demonstrate the value of in vivo MRS as a technique that determines parameters of mitochondrial dysfunction directly and indirectly which could be of value for the diagnosis, prognosis, and follow-up of mitochondrial disorders.     

In vivo magnetic resonance imaging tracking of SPIO-labeled human umbilical cord mesenchymal stem cells

Human umbilical cord mesenchymal stem cells (hUC-MSCs) can be efficiently labeled by superparamagnetic iron oxide (SPIO) nanoparticles, which produces low signal intensity on magnetic resonance imaging (MRI) in vitro. This study was to evaluate the feasibility of in vivo tracking for hUC-MSCs labeled by SPIO with noninvasive MRI. SPIO was added to cultures at concentrations equivalent to 0, 7, 14, 28, and 56 μg Fe/ml (diluted with DMEM/F12) and incubated for 16 h. Prussian Blue staining was used to determinate the labeling efficiency. Rats were randomly divided into three groups, control group, hUC-MSCs group, and SPIO-labeled hUC-MSCs group. All groups were subjected to spinal cord injury (SCI) by weight drop device. Rats were examined for neurological function. In vivo MRI was used to track SPIO-labeled hUC-MSCs transplanted in rats spinal cord. Survival and migration of hUC-MSCs were also explored using immunofluorescence. Significant improvements in locomotion were observed in the hUC-MSCs groups. There was statistical significance compared with control group. In vivo MRI 1 and 3 weeks after injection showed a large reduction in signal intensity in the region transplanted with SPIO-labeled hUC-MSCs. The images from unlabeled hUC-MSCs showed a smaller reduction in signal intensity. Transplanted hUC-MSCs engrafted within the injured rats spinal cord and survived for at least 8 weeks. In conclusion, hUC-MSCs can survive and migrate in the host spinal cord after transplantation, which promote functional recovery after SCI. Noninvasive imaging of transplanted SPIO-labeled hUC-MSCs is feasible.