Neural network for autonomous segmentation and volumetric assessment of clot and edema in acute and subacute intracerebral hemorrhages

INTRODUCTION

Minimally-invasive surgical techniques for intracerebral hemorrhage (ICH) evacuation use imaging to guide the suction, lysing and/or drainage from the hemorrhage site via various designs. A previous international surgical study has shown that reduction of hematoma volume below 15 ml is indicative of improved long term patient outcomes. The study noted a need for tools to periodically visualize remaining clot during intervention to increase the likelihood of evacuating sufficient clot volumes without endangering rebleeds. Robust segmentation of MRI could guide surgeons and radiologists regarding remaining regions and approaches for prudent evacuation. We thus propose a Convolutional Neural Network (CNN) to identify and autonomously segment clot and peripheral edema in MR images of the brain and generate an estimate of the remaining clot volume.

MATERIALS AND METHODS

We used a retrospective, locally-acquired dataset of ICH patient scans taken on 3 T MRI scanners. Three sets of ground truth manual segmentations were independently generated by two imaging scientists and one radiology fellow. Evaluation of clot age was determined based on relative contrast of hemorrhage components and reviewed by a neurosurgeon. Model accuracy was determined by pixel-wise Dice coefficient (DC) calculations between each ground truth manual segmentation and the machine-derived autonomous segmentations.

RESULTS

The model produced autonomous segmentations of clot core with an average DC of 0.75 ± 0.21 relative to manual segmentations of the same scans. For edema, it produced segmentations with an average DC of 0.68 ± 0.16 relative to manual. From these pixel-wise segmentations, clot volume can be calculated. Model-produced segmentations underestimated clot volumes by an average of 17% relative to ground-truth.

CONCLUSION

The machine learning models were able to identify and segment volumes of ICH components swiftly and accurately.

Real-time trajectory guide tracking for intraoperative MRI-guided neurosurgery

PURPOSE

In current intraoperative MRI (IMRI) methods, an iterative approach is used to aim trajectory guides at intracerebral targets: image MR-visible features, determine current aim by fitting model to image, manipulate device, repeat. Infrequent updates are produced by such methods, compared to rapid optically tracked stereotaxy used in the operating room. Our goal was to develop a real-time interactive IMRI method for aiming.

METHODS

The current trajectory was computed from two points along the guide's central axis, rather than by imaging the entire device. These points were determined by correlating one-dimensional spokes from a radial sequence with the known cross-sectional projection of the guide. The real-time platform RTHawk was utilized to control MR sequences and data acquisition. On-screen updates were viewed by the operator while simultaneously manipulating the guide to align it with the planned trajectory. Accuracy was quantitated in a phantom, and in vivo validation was demonstrated in nonhuman primates undergoing preclinical gene (n = 5 $$ n=5 $$ ) and cell (n = 4 $$ n=4 $$ ) delivery surgeries.

RESULTS

Updates were produced at 5 Hz In 10 phantom experiments at a depth of 48 mm, the cannula tip was placed with radial error of (min, mean, max) = (0.16, 0.29, 0.68) mm. Successful in vivo delivery of payloads to all 14 targets was demonstrated across nine surgeries with depths of (min, mean, max) = (33.3, 37.9, 42.5) mm.

CONCLUSION

A real-time interactive update rate was achieved, reducing operator fatigue without compromising accuracy. Qualitative interpretation of images during aiming was rendered unnecessary by objectively computing device alignment.

Multiscale Label-Free Imaging of Fibrillar Collagen in the Tumor Microenvironment

With recent advances in cancer therapeutics, there is a great need for improved imaging methods for characterizing cancer onset and progression in a quantitative and actionable way. Collagen, the most abundant extracellular matrix protein in the tumor microenvironment (and the body in general), plays a multifaceted role, both hindering and promoting cancer invasion and progression. Collagen deposition can defend the tumor with immunosuppressive effects, while aligned collagen fiber structures can enable tumor cell migration, aiding invasion and metastasis. Given the complex role of collagen fiber organization and topology, imaging has been a tool of choice to characterize these changes on multiple spatial scales, from the organ and tumor scale to cellular and subcellular level. Macroscale density already aids in the detection and diagnosis of solid cancers, but progress is being made to integrate finer microscale features into the process. Here we review imaging modalities ranging from optical methods of second harmonic generation (SHG), polarized light microscopy (PLM), and optical coherence tomography (OCT) to the medical imaging approaches of ultrasound and magnetic resonance imaging (MRI). These methods have enabled scientists and clinicians to better understand the impact collagen structure has on the tumor environment, at both the bulk scale (density) and microscale (fibrillar structure) levels. We focus on imaging methods with the potential to both examine the collagen structure in as natural a state as possible and still be clinically amenable, with an emphasis on label-free strategies, exploiting intrinsic optical properties of collagen fibers.

[18F]FEPPA PET imaging for monitoring CD68-positive microglia/macrophage neuroinflammation in nonhuman primates

PURPOSE

The aim of this study was to examine whether the translocator protein 18-kDa (TSPO) PET ligand [18F]FEPPA has the sensitivity for detecting changes in CD68-positive microglial/macrophage activation in hemiparkinsonian rhesus macaques treated with allogeneic grafts of induced pluripotent stem cell-derived midbrain dopaminergic neurons (iPSC-mDA).

METHODS

In vivo positron emission tomography (PET) imaging with [18F]FEPPA was used in conjunction with postmortem CD68 immunostaining to evaluate neuroinflammation in the brains of hemiparkinsonian rhesus macaques (n = 6) that received allogeneic iPSC-mDA grafts in the putamen ipsilateral to MPTP administration.

RESULTS

Based on assessment of radiotracer uptake and confirmed by visual inspection of the imaging data, nonhuman primates with allogeneic grafts showed increased [18F]FEPPA binding at the graft sites relative to the contralateral putamen. From PET asymmetry analysis of the images, the mean asymmetry index of the monkeys was AI = - 0.085 ± 0.018. Evaluation and scoring of CD68 immunoreactivity by an investigator blind to the treatment identified significantly more neuroinflammation in the grafted areas of the putamen compared to the contralateral putamen (p = 0.0004). [18F]FEPPA PET AI showed a positive correlation with CD68 immunoreactivity AI ratings in the monkeys (Spearman's ρ = 0.94; p = 0.005).

CONCLUSION

These findings reveal that [18F]FEPPA PET is an effective marker for detecting increased CD68-positive microglial/macrophage activation and demonstrates sufficient sensitivity to detect changes in neuroinflammation in vivo following allogeneic cell engraftment.

Dorsal Amygdala Neurotrophin-3 Decreases Anxious Temperament in Primates

BACKGROUND

An early-life anxious temperament (AT) is a risk factor for the development of anxiety, depression, and comorbid substance abuse. We validated a nonhuman primate model of early-life AT and identified the dorsal amygdala as a core component of AT's neural circuit. Here, we combine RNA sequencing, viral-vector gene manipulation, functional brain imaging, and behavioral phenotyping to uncover AT's molecular substrates.

METHODS

In response to potential threat, AT and brain metabolism were assessed in 46 young rhesus monkeys. We identified AT-related transcripts using RNA-sequencing data from dorsal amygdala tissue (including central nucleus of the amygdala [Ce] and dorsal regions of the basal nucleus). Based on the results, we overexpressed the neurotrophin-3 gene, NTF3, in the dorsal amygdala using intraoperative magnetic resonance imaging-guided surgery (n = 5 per group).

RESULTS

This discovery-based approach identified AT-related alterations in the expression of well-established and novel genes, including an inverse association between NTRK3 expression and AT. NTRK3 is an interesting target because it is a relatively unexplored neurotrophic factor that modulates intracellular neuroplasticity pathways. Overexpression of the transcript for NTRK3's endogenous ligand, NTF3, in the dorsal amygdala resulted in reduced AT and altered function in AT's neural circuit.

CONCLUSIONS

Together, these data implicate neurotrophin-3/NTRK3 signaling in the dorsal amygdala in mediating primate anxiety. More generally, this approach provides an important step toward understanding the molecular underpinnings of early-life AT and will be useful in guiding the development of treatments to prevent the development of stress-related psychopathology.

Use of a platelet-rich plasma-collagen scaffold as a bioenhanced repair treatment for management of partial cruciate ligament rupture in dogs

Dogs are commonly affected with cruciate ligament rupture (CR) and associated osteoarthritis (OA), and frequently develop a second contralateral CR. Platelet rich plasma (PRP) is a component of whole blood that contains numerous growth factors, which in combination with a collagen scaffold may act to promote bioenhanced primary repair of ligament. This study tested the hypothesis that treatment of partial stable CR stifles with an intra-articular collagen scaffold and PRP would decrease the disease progression, synovitis and risk of complete CR over a 12-month study period. We conducted a prospective cohort study of 29 client-owned dogs with an unstable stifle due to complete CR and stable contralateral stifle with partial CR. All dogs were treated with tibial plateau leveling osteotomy (TPLO) on the unstable stifle and a single intra-articular application of PRP-collagen in the stable partial CR stifle. Dogs were evaluated at the time of diagnosis, and at 10-weeks and 12-months after treatment. We evaluated correlation between both development of complete CR and time to complete CR with diagnostic tests including bilateral stifle radiographs, 3.0 Tesla magnetic resonance (MR) imaging, and bilateral stifle arthroscopy. Additionally, histologic evaluation of synovial biopsies, C-reactive protein (CRP) concentrations in serum and synovial fluid, and synovial total nucleated cell count, were determined. Results indicated that a single application of PRP-collagen in partial CR stifles of client owned dogs is not an effective disease-modifying therapy for the prevention of progression to complete CR. Radiographic effusion, arthroscopic evaluation of cranial cruciate ligament (CrCL) damage, and MR assessment of ligament fiber tearing in partial CR stifles correlated with progression to complete CR over the 12-month follow-up period. We determined that the best predictive model for development of complete CR in PRP-collagen treated partial CR stifles included variables from multiple diagnostic modalities.

Feasibility of high spatiotemporal resolution for an abbreviated 3D radial breast MRI protocol

PURPOSE

To develop a volumetric imaging technique with 0.8-mm isotropic resolution and 10-s/volume rate to detect and analyze breast lesions in a bilateral, dynamic, contrast-enhanced MRI exam.

METHODS

A local low-rank temporal reconstruction approach that also uses parallel imaging and spatial compressed sensing was designed to create rapid volumetric frame rates during a contrast-enhanced breast exam (vastly undersampled isotropic projection [VIPR] spatial compressed sensing with temporal local low-rank [STELLR]). The dynamic-enhanced data are subtracted in k-space from static mask data to increase sparsity for the local low-rank approach to maximize temporal resolution. A T₁ -weighted 3D radial trajectory (VIPR iterative decomposition with echo asymmetry and least squares estimation [IDEAL]) was modified to meet the data acquisition requirements of the STELLR approach. Additionally, the unsubtracted enhanced data are reconstructed using compressed sensing and IDEAL to provide high-resolution fat/water separation. The feasibility of the approach and the dual reconstruction methodology is demonstrated using a 16-channel breast coil and a 3T MR scanner in 6 patients.

RESULTS

The STELLR temporal performance of subtracted data matched the expected temporal perfusion enhancement pattern in small and large vascular structures. Differential enhancement within heterogeneous lesions is demonstrated with corroboration from a basic reconstruction using a strict 10-second temporal footprint. Rapid acquisition, reliable fat suppression, and high spatiotemporal resolution are presented, despite significant data undersampling.

CONCLUSION

The STELLR reconstruction approach of 3D radial sampling with mask subtraction provides a high-performance imaging technique for characterizing enhancing structures within the breast. It is capable of maintaining temporal fidelity, while visualizing breast lesions with high detail over a large FOV to include both breasts.

Review of quantitative multiscale imaging of breast cancer

Breast cancer is the most common cancer among women worldwide and ranks second in terms of overall cancer deaths. One of the difficulties associated with treating breast cancer is that it is a heterogeneous disease with variations in benign and pathologic tissue composition, which contributes to disease development, progression, and treatment response. Many of these phenotypes are uncharacterized and their presence is difficult to detect, in part due to the sparsity of methods to correlate information between the cellular microscale and the whole-breast macroscale. Quantitative multiscale imaging of the breast is an emerging field concerned with the development of imaging technology that can characterize anatomic, functional, and molecular information across different resolutions and fields of view. It involves a diverse collection of imaging modalities, which touch large sections of the breast imaging research community. Prospective studies have shown promising results, but there are several challenges, ranging from basic physics and engineering to data processing and quantification, that must be met to bring the field to maturity. This paper presents some of the challenges that investigators face, reviews currently used multiscale imaging methods for preclinical imaging, and discusses the potential of these methods for clinical breast imaging.

Radiographic and magnetic resonance imaging predicts severity of cruciate ligament fiber damage and synovitis in dogs with cranial cruciate ligament rupture

Cruciate ligament rupture (CR) and associated osteoarthritis (OA) is a common condition in dogs. Dogs frequently develop a second contralateral CR. This study tested the hypothesis that the degree of stifle synovitis and cranial cruciate ligament (CrCL) matrix damage in dogs with CR is correlated with non-invasive diagnostic tests, including magnetic resonance (MR) imaging. We conducted a prospective cohort study of 29 client-owned dogs with an unstable stifle due to complete CR and stable contralateral stifle with partial CR. We evaluated correlation of stifle synovitis and CrCL fiber damage with diagnostic tests including bilateral stifle radiographs, 3.0 Tesla MR imaging, and bilateral stifle arthroscopy. Histologic grading and immunohistochemical staining for CD3⁺ T lymphocytes, TRAP⁺ activated macrophages and Factor VIII⁺ blood vessels in bilateral stifle synovial biopsies were also performed. Serum and synovial fluid concentrations of C-reactive protein (CRP) and carboxy-terminal telopeptide of type I collagen (ICTP), and synovial total nucleated cell count were determined. Synovitis was increased in complete CR stifles relative to partial CR stifles (P<0.0001), although total nucleated cell count in synovial fluid was increased in partial CR stifles (P<0.01). In partial CR stifles, we found that 3D Fast Spin Echo Cube CrCL signal intensity was correlated with histologic synovitis (S_R = 0.50, P<0.01) and that radiographic OA was correlated with CrCL fiber damage assessed arthroscopically (S_R = 0.61, P<0.001). Taken together, results of this study show that clinical diagnostic tests predict severity of stifle synovitis and cruciate ligament matrix damage in stable partial CR stifles. These data support use of client-owned dogs with unilateral complete CR and contralateral partial CR as a clinical trial model for investigation of disease-modifying therapy for partial CR.

Fast Realistic MRI Simulations Based on Generalized Multi-Pool Exchange Tissue Model

We present MRiLab, a new comprehensive simulator for large-scale realistic MRI simulations on a regular PC equipped with a modern graphical processing unit (GPU). MRiLab combines realistic tissue modeling with numerical virtualization of an MRI system and scanning experiment to enable assessment of a broad range of MRI approaches including advanced quantitative MRI methods inferring microstructure on a sub-voxel level. A flexible representation of tissue microstructure is achieved in MRiLab by employing the generalized tissue model with multiple exchanging water and macromolecular proton pools rather than a system of independent proton isochromats typically used in previous simulators. The computational power needed for simulation of the biologically relevant tissue models in large 3D objects is gained using parallelized execution on GPU. Three simulated and one actual MRI experiments were performed to demonstrate the ability of the new simulator to accommodate a wide variety of voxel composition scenarios and demonstrate detrimental effects of simplified treatment of tissue micro-organization adapted in previous simulators. GPU execution allowed  ∼ 200× improvement in computational speed over standard CPU. As a cross-platform, open-source, extensible environment for customizing virtual MRI experiments, MRiLab streamlines the development of new MRI methods, especially those aiming to infer quantitatively tissue composition and microstructure.

Real-Time Intraoperative MRI Intracerebral Delivery of Induced Pluripotent Stem Cell-Derived Neurons

Induced pluripotent stem cell (iPSC)-derived neurons represent an opportunity for cell replacement strategies for neurodegenerative disorders such as Parkinson's disease (PD). Improvement in cell graft targeting, distribution, and density can be key for disease modification. We have previously developed a trajectory guide system for real-time intraoperative magnetic resonance imaging (RT-IMRI) delivery of infusates, such as viral vector suspensions for gene therapy strategies. Intracerebral delivery of iPSC-derived neurons presents different challenges than viral vectors, including limited cell survival if cells are kept at room temperature for prolonged periods of time, precipitation and aggregation of cells in the cannula, and obstruction during injection, which must be solved for successful application of this delivery approach. To develop procedures suitable for RT-IMRI cell delivery, we first performed in vitro studies to tailor the delivery hardware (e.g., cannula) and defined a range of parameters to be applied (e.g., maximal time span allowable between cell loading in the system and intracerebral injection) to ensure cell survival. Then we performed an in vivo study to evaluate the feasibility of applying the system to nonhuman primates. Our results demonstrate that the RT-IMRI delivery system provides valuable guidance, monitoring, and visualization during intracerebral cell delivery that are compatible with cell survival.

Overexpressing Corticotropin-Releasing Factor in the Primate Amygdala Increases Anxious Temperament and Alters Its Neural Circuit

BACKGROUND

Nonhuman primate models are critical for understanding mechanisms underlying human psychopathology. We established a nonhuman primate model of anxious temperament (AT) for studying the early-life risk to develop anxiety and depression. Studies have identified the central nucleus of the amygdala (Ce) as an essential component of AT's neural substrates. Corticotropin-releasing factor (CRF) is expressed in the Ce, has a role in stress, and is linked to psychopathology. Here, in young rhesus monkeys, we combined viral vector technology with assessments of anxiety and multimodal neuroimaging to understand the consequences of chronically increased CRF in the Ce region.

METHODS

Using real-time intraoperative magnetic resonance imaging-guided convection-enhanced delivery, five monkeys received bilateral dorsal amygdala Ce-region infusions of adeno-associated virus serotype 2 containing the CRF construct. Their cagemates served as unoperated control subjects. AT, regional brain metabolism, resting functional magnetic resonance imaging, and diffusion tensor imaging were assessed before and 2 months after viral infusions.

RESULTS

Dorsal amygdala CRF overexpression significantly increased AT and metabolism within the dorsal amygdala. Additionally, we observed changes in metabolism in other AT-related regions, as well as in measures of functional and structural connectivity.

CONCLUSIONS

This study provides a translational roadmap that is important for understanding human psychopathology by combining molecular manipulations used in rodents with behavioral phenotyping and multimodal neuroimaging measures used in humans. The results indicate that chronic CRF overexpression in primates not only increases AT but also affects metabolism and connectivity within components of AT's neural circuitry.

Multicomponent T2 analysis of articular cartilage with synovial fluid partial volume correction

PURPOSE

To investigate the use of a three-pool model to account for the confounding effects of synovial fluid on multicomponent T2 analysis of articular cartilage using Multicomponent Driven Equilibrium Single Shot Observation of T1 and T2 (mcDESPOT).

MATERIALS AND METHODS

mcDESPOT was performed on the knee of eight asymptomatic volunteers and eight patients with osteoarthritis at 3.0T with multicomponent T2 maps created using the two-pool model and a three-pool model containing a nonexchanging synovial fluid water pool. The fraction of the fast-relaxing water component (FF ) and the T2 relaxation times for the fast-relaxing (T2F ) and slow-relaxing (T2S ) water components were measured in the superficial and deep layers of patellar cartilage using the two-pool and three-pool models in asymptomatic volunteers and patients with osteoarthritis and were compared using Wilcoxon signed rank tests.

RESULTS

Within the superficial layer of patellar cartilage, FF was 22.5% and 25.6% for asymptomatic volunteers and 21.3% and 22.8% for patients with osteoarthritis when using the two-pool and three-pool models, respectively, while T2S was 73.9 msec and 62.0 msec for asymptomatic volunteers and 72.0 msec and 63.1 msec for patients with osteoarthritis when using the two-pool and three-pool models, respectively. For both asymptomatic volunteers and patients with osteoarthritis, the two-pool model provided significantly (P < 0.05) lower FF and higher T2S than the three-pool model, likely due to the effects of synovial fluid partial volume averaging.

CONCLUSION

The effects of partial volume averaging between superficial cartilage and synovial fluid may result in biased multicomponent T2 measurements that can be corrected using an mcDESPOT three-pool model containing a nonexchanging synovial fluid water pool.

Rapid multicomponent relaxometry in steady state with correction of magnetization transfer effects

PURPOSE

To study the effects of magnetization transfer (MT) on multicomponent T2 parameters obtained using mcDESPOT in macromolecule-rich tissues and to propose a new method called mcRISE to correct MT-induced biases.

METHODS

The two-pool mcDESPOT model was modified by the addition of an exchanging macromolecule proton pool to model the MT effect in cartilage. The mcRISE acquisition scheme was developed to provide sensitivity to all pools. An incremental fitting was applied to estimate MT and relaxometry parameters with minimized coupling. The interaction between MT and relaxometry parameters, efficacy of MT correction, and feasibility of mcRISE in vivo were investigated in simulations and in healthy volunteers.

RESULTS

The MT effect caused significant errors in multicomponent T1/T2 values and in fast-relaxing water fraction fF , which is consistent with previous experimental observations. fF increased significantly with macromolecule content if MT was ignored. mcRISE resulted in a multifold reduction of MT biases and yielded decoupled multicomponent T1/T2 relaxometry and quantitative MT parameters.

CONCLUSION

mcRISE is an efficient approach for correcting MT biases in multicomponent relaxometry based on steady state sequences. Improved specificity of mcRISE may help to elucidate the sources of the previously described high sensitivity of noncorrected mcDESPOT parameters to disease-related changes in cartilage and the brain.

3D FSE Cube and VIPR-aTR 3.0 Tesla magnetic resonance imaging predicts canine cranial cruciate ligament structural properties

Estimation of cranial cruciate ligament (CrCL) structural properties in client-owned dogs with incipient cruciate rupture would be advantageous. The objective of this study was to determine whether magnetic resonance imaging (MRI) measurement of normal CrCL volume in an ex-vivo canine model predicts structural properties. Stifles from eight dogs underwent 3.0 Tesla 3D MRI. CrCL volume and normalized median grayscale values were determined using 3D Fast Spin Echo (FSE) Cube and Vastly under-sampled Isotropic PRojection (VIPR)-alternative repetition time (aTR) sequences. Stifles were then mechanically tested. After joint laxity testing, CrCL structural properties were determined, including displacement at yield, yield load, load to failure, and stiffness. Yield load and load to failure (R(2)=0.56, P <0.01) were correlated with CrCL volume determined by VIPR-aTR. Yield load was also correlated with CrCL volume determined by 3D FSE Cube (R(2)=0.32, P <0.05). Structural properties were not related to median grayscale values. Joint laxity and CrCL stiffness were not related to MRI parameters, but displacement at yield load was related to CrCL volume for both sequences during testing (R(2)>0.57, P <0.005). In conclusion, 3D MRI offers a predictive method for estimating canine CrCL structural properties. 3D MRI may be useful for monitoring CrCL properties in clinical trials.

High-spatial and high-temporal resolution dynamic contrast-enhanced perfusion imaging of the liver with time-resolved three-dimensional radial MRI

PURPOSE

Detection, characterization, and monitoring the treatment of hepatocellular carcinomas (HCC) in patients with cirrhosis is challenging because of their variable and rapid arterial enhancement. Multiphase dynamic contrast-enhanced MRI is used clinically for HCC assessment; however, the method suffers from limited temporal resolution and difficulty in coordinating imaging and breath-hold timing within a narrow temporal window of interest. In this article, a volumetric, high-spatial resolution, and high-temporal resolution dynamic contrast-enhanced liver imaging method for improved detection and characterization of HCC is demonstrated.

METHODS

A time-resolved three-dimensional radial acquisition with iterative sensitivity-encoding reconstruction images the entire abdomen and thorax with high spatial and temporal resolution, using real-time three-dimensional fluoroscopy to match the breath hold to contrast arrival. The sequence was tested on 17 subjects, including eight patients with HCC or other hypervascular focal lesions.

RESULTS

This technique was successful in acquiring volumetric imaging of the entire liver with 2.1-mm isotropic spatial and true 4-s temporal resolution.

CONCLUSION

This technique may be suitable for detecting, characterizing, and monitoring the treatment of HCC. It also holds significant potential for perfusion modeling, which may provide a noninvasive means to rapidly determine the efficacy of chemotherapeutic agents in these tumors over the entire liver volume.

Articular Cartilage of the Human Knee Joint: In Vivo Multicomponent T2 Analysis at 3.0 T

PURPOSE

To compare multicomponent T2 parameters of the articular cartilage of the knee joint measured by using multicomponent driven equilibrium single-shot observation of T1 and T2 (mcDESPOT) in asymptomatic volunteers and patients with osteoarthritis.

MATERIALS AND METHODS

This prospective study was performed with institutional review board approval and with written informed consent from all subjects. The mcDESPOT sequence was performed in the knee joint of 13 asymptomatic volunteers and 14 patients with osteoarthritis of the knee. Single-component T2 (T2(Single)), T2 of the fast-relaxing water component (T2F) and of the slow-relaxing water component (T2S), and the fraction of the fast-relaxing water component (F(F)) of cartilage were measured. Wilcoxon rank-sum tests and multivariate linear regression models were used to compare mcDESPOT parameters between volunteers and patients with osteoarthritis. Receiver operating characteristic analysis was used to assess diagnostic performance with mcDESPOT parameters for distinguishing morphologically normal cartilage from morphologically degenerative cartilage identified at magnetic resonance imaging in eight cartilage subsections of the knee joint.

RESULTS

Higher cartilage T2(Single) (P < .001), lower cartilage F(F) (P < .001), and similar cartilage T2F (P = .079) and T2S (P = .124) values were seen in patients with osteoarthritis compared with those in asymptomatic volunteers. Differences in T2(Single) and F(F) remained significant (P < .05) after consideration of age differences between groups of subjects. Diagnostic performance was higher with F(F) than with T2(Single) for distinguishing between normal and degenerative cartilage (P < .05), with greater areas under the curve at receiver operating characteristic analysis.

CONCLUSION

Patients with osteoarthritis of the knee had significantly higher cartilage T2(Single) and significantly lower cartilage F(F) than did asymptomatic volunteers, and receiver operating characteristic analysis results suggested that F(F) may allow greater diagnostic performance than that with T2(Single) for distinguishing between normal and degenerative cartilage.

Utilization of a balanced steady state free precession signal model for improved fat/water decomposition

PURPOSE

Chemical shift based fat/water decomposition methods such as IDEAL are frequently used in challenging imaging environments with large B0 inhomogeneity. However, they do not account for the signal modulations introduced by a balanced steady state free precession (bSSFP) acquisition. Here we demonstrate improved performance when the bSSFP frequency response is properly incorporated into the multipeak spectral fat model used in the decomposition process.

THEORY AND METHODS

Balanced SSFP allows for rapid imaging but also introduces a characteristic frequency response featuring periodic nulls and pass bands. Fat spectral components in adjacent pass bands will experience bulk phase offsets and magnitude modulations that change the expected constructive and destructive interference between the fat spectral components. A bSSFP signal model was incorporated into the fat/water decomposition process and used to generate images of a fat phantom, and bilateral breast and knee images in four normal volunteers at 1.5 Tesla.

RESULTS

Incorporation of the bSSFP signal model into the decomposition process improved the performance of the fat/water decomposition.

CONCLUSION

Incorporation of this model allows rapid bSSFP imaging sequences to use robust fat/water decomposition methods such as IDEAL. While only one set of imaging parameters were presented, the method is compatible with any field strength or repetition time.

Rapid in vivo multicomponent T2 mapping of human knee menisci

PURPOSE

To compare multicomponent T2 parameters of menisci measured using Multicomponent Driven Equilibrium Single Pulse Observation of T1 and T2 (mcDESPOT) in asymptomatic volunteers and osteoarthritis (OA) patients with intact and torn menisci.

MATERIALS AND METHODS

The prospective study was performed with Institutional Review Board approval and with all subjects signing written informed consent. mcDESPOT was performed on the knee joint of 12 asymptomatic volunteers and 14 patients with knee OA. Single-component T2 relaxation time (T2Single ), T2 relaxation time of the fast relaxing water component (T2F ), and the slow relaxing water component (T2S ), and fraction of the fast relaxing water component (FF ) of the medial and lateral menisci were measured. Multivariate linear regression models were used to compare mcDESPOT parameters between normal menisci in asymptomatic volunteers, intact menisci in OA patients, and torn menisci in OA patients with adjustment for differences in age between subjects.

RESULTS

The mean mcDESPOT parameters for normal menisci in asymptomatic volunteers, intact menisci in OA patients, and torn menisci in OA patients were respectively 16.1 msec, 18.8 msec, and 22.7 msec for T2Single ; 9.0 msec, 10.0 msec, and 11.1 msec for T2F ; 24.4 msec, 27.7 msec, and 31.4 msec for T2S ; and 34%, 32%, 27% for FF . There were significant differences (P < 0.05) in T2Single , T2F , T2S , and FF between the three groups of menisci.

CONCLUSION

The menisci of OA patients had significantly higher T2Single , T2F , and T2S and significantly lower FF than normal menisci in asymptomatic volunteers with greater changes in multicomponent T2 parameters noted in torn than intact menisci in OA patients.

Articular cartilage zonal differentiation via 3D Second-Harmonic Generation imaging microscopy

PURPOSE

The collagen structure throughout the patella has not been thoroughly investigated by 3D imaging, where the majority of the existing data come from histological cross sections. It is important to have a better understanding of the architecture in normal tissues, where this could then be applied to imaging of diseased states.

METHODS

To address this shortcoming, we investigated the combined use of collagen-specific Second-Harmonic Generation (SHG) imaging and measurement of bulk optical properties to characterize collagen fiber orientations of the histologically defined zones of bovine articular cartilage. Forward and backward SHG intensities of sections from superficial, middle and deep zones were collected as a function of depth and analyzed by Monte Carlo simulations to extract the SHG creation direction, which is related to the fibrillar assembly.

RESULTS

Our results revealed differences in SHG forward-backward response between the three zones, where these are consistent with a previously developed model of SHG emission. Some of the findings are consistent with that from other modalities; however, SHG analysis showed the middle zone had the most organized fibril assembly. While not distinct, we also report bulk optical property values for these different zones within the patella.

CONCLUSIONS

Collectively, these results provide quantitative measurements of structural changes at both the fiber and fibril assembly of the different cartilage zones and reveals structural information not possible by other microscope modalities. This can provide quantitative insight to the collagen fiber network in normal cartilage, which may ultimately be developed as a biomarker for osteoarthritis.

Cross-relaxation imaging of human patellar cartilage in vivo at 3.0T

OBJECTIVE

To compare quantitative magnetization transfer (qMT) parameters of patellar cartilage measured using cross-relaxation imaging (CRI) in asymptomatic volunteers and patients with osteoarthritis.

DESIGN

The study was performed with Institutional Review Board approval and with all subjects signing informed consent. CRI of the knee joint was performed at 3.0T on 20 asymptomatic volunteers and 11 patients with osteoarthritis. The fraction of macromolecular bound protons (f), the exchange rate constant between macromolecular bound protons and free water protons (k), and the T2 relaxation time of macromolecular bound protons (T2(B)) of patellar cartilage were measured. Mann-Whitney-Wilcoxon rank-sum tests were used to compare qMT parameters between asymptomatic volunteers and patients with osteoarthritis.

RESULTS

Average f, k, and T2(B) of patellar cartilage was 12.46%, 7.22 s(-1), and 6.49 μs respectively for asymptomatic volunteers and 12.80%, 6.13 s(-1), and 6.80 μs respectively for patients with osteoarthritis. There were statistically significant differences between groups of subjects for k (P < 0.01) and T2(B) (P < 0.0001) but not f (P = 0.38) of patellar cartilage.

CONCLUSION

Patients with osteoarthritis had significantly lower k and significantly higher T2(B) of patellar cartilage than asymptomatic volunteers which suggests that qMT parameters can detect changes in the macromolecular matrix of degenerative cartilage.

Quantitative hepatic perfusion modeling using DCE-MRI with sequential breathholds

PURPOSE

To develop and demonstrate the feasibility of a new formulation for quantitative perfusion modeling in the liver using interrupted DCE-MRI data acquired during multiple sequential breathholds.

MATERIALS AND METHODS

A new mathematical formulation to estimate quantitative perfusion parameters using interrupted data was developed. Using this method, we investigated whether a second degree-of-freedom in the tissue residue function (TRF) improves quality-of-fit criteria when applied to a dual-input single-compartment perfusion model. We subsequently estimated hepatic perfusion parameters using DCE-MRI data from 12 healthy volunteers and 9 cirrhotic patients with a history of hepatocellular carcinoma (HCC); and examined the utility of these estimates in differentiating between healthy liver, cirrhotic liver, and HCC.

RESULTS

Quality-of-fit criteria in all groups were improved using a Weibull TRF (2 degrees-of-freedom) versus an exponential TRF (1 degree-of-freedom), indicating nearer concordance of source DCE-MRI data with the Weibull model. Using the Weibull TRF, arterial fraction was greater in cirrhotic versus normal liver (39 ± 23% versus 15 ± 14%, P = 0.07). Mean transit time (20.6 ± 4.1 s versus 9.8 ± 3.5 s, P = 0.01) and arterial fraction (39 ± 23% versus 73 ± 14%, P = 0.04) were both significantly different between cirrhotic liver and HCC, while differences in total perfusion approached significance.

CONCLUSION

This work demonstrates the feasibility of estimating hepatic perfusion parameters using interrupted data acquired during sequential breathholds.

Rapid isotropic resolution cartilage assessment using radial alternating repetition time balanced steady-state free-precession imaging

PURPOSE

To compare a balanced steady-state free-precession sequence with a radial k-space trajectory and alternating repetition time fat suppression (Radial-ATR) with other currently used fat-suppressed 3D sequences for evaluating the articular cartilage of the knee joint at 3.0T.

MATERIALS AND METHODS

Radial-ATR, fast spin-echo (FSE-Cube), gradient recall-echo acquired in the steady-state (GRASS), and spoiled gradient recall-echo (SPGR) sequences with similar voxel volumes and identical scan times were performed at 3.0T on both knee joints of five volunteers. Signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) measurements were performed for all sequences using a double acquisition method and compared using Mann-Whitney Wilcoxon tests. Radial-ATR sequences with 0.3 mm and 0.4 mm isotropic resolution were also performed on the knee joints of seven volunteers and three patients with osteoarthritis.

RESULTS

Average SNR values for cartilage, synovial fluid, and bone marrow were 54.7, 153.3, and 12.9, respectively, for Radial ATR, 30.8, 44.1, and 1.9, respectively, for FSE-Cube, 13.3, 46.9, and 3.3, respectively, for GRASS, and 19.1, 8.1, and 2.1, respectively, for SPGR. Average CNR values between cartilage and synovial fluid and between cartilage and bone marrow were 98.6 and 41.8, respectively, for VIPR-ATR, 13.4 and 28.8, respectively, for FSE-Cube, 33.6 and 10.0, respectively, for GRASS, and 11.0 and 16.9, respectively, for SPGR. Radial-ATR had significantly higher (P < 0.001) cartilage, synovial fluid, and bone marrow SNR and significantly higher (P < 0.01) CNR between cartilage and synovial fluid and between cartilage and bone marrow than FSE-Cube, GRASS, and SPGR. Radial-ATR provided excellent visualization of articular cartilage at high isotropic resolution with no image degradation due to off-resonance banding artifacts.

CONCLUSION

Radial-ATR had superior SNR efficiency to other fat-suppressed 3D cartilage imaging sequences and produced high isotropic resolution images of the knee joint which could be used for evaluating articular cartilage at 3.0T.

Non-contrast enhanced 3D SSFP MRA of the renal allograft vasculature: a comparison between radial linear combination and Cartesian inflow-weighted acquisitions

Renal transplant patients often require imaging to ensure appropriate graft placement, to assess integrity of transplant vessel anastomosis and to evaluate for stenosis that can be a cause of graft failure. Because there is risk for nephrogenic systemic fibrosis in the setting of renal insufficiency, the use of non-contrast MRA in these patients is helpful. In this study, the ability of two non-contrast MRA methods - 3D radial linear combination balanced SSFP (VIPR-SSFP) and inflow-weighted Cartesian SSFP (IFIR) - to visualize the transplant renal vessels is compared. Twenty-one renal transplant patients were scanned using the VIPR-SSFP and IFIR sequences. Diagnostic efficacy of the sequences was scored using a four point Likert scale according to the following criteria: overall image quality, fat suppression, and arterial/venous visualization quality. Average scores for each criterion were compared using the Wilcoxon signed-rank test. In addition to significantly improved venous visualization, the VIPR-SSFP sequence provided significantly improved fat suppression quality (p<0.03) compared to IFIR. VIPR-SSFP also identified several pathologies such as renal arterial pseudoaneurysm that were not visible on the IFIR images. However, IFIR afforded superior quality of arterial visualization (p<0.005). These two methods of non-contrast MR imaging each have significant strengths and are complementary to each other in evaluating the vasculature of renal allografts.

Noncontrast dynamic 3D intracranial MR angiography using pseudo-continuous arterial spin labeling (PCASL) and accelerated 3D radial acquisition

PURPOSE

To develop a novel dynamic 3D noncontrast magnetic resonance angiography (MRA) technique that combines dynamic pseudo-continuous arterial spin labeling (dynamic PCASL), accelerated 3D radial sampling (VIPR), and time-of-arrival (TOA) mapping to provide quantitative assessment of arterial flow.

MATERIALS AND METHODS

Digital simulations were performed to investigate the effects of acquisition scheme and sequence parameters on image quality and TOA mapping fidelity. Five patients with vascular malformations (arteriovenous malformation [AVM] = 3, dural arteriovenous fistula [DAVF] = 2) were scanned and the images were compared to digital subtraction angiography (DSA) for the ability to identify the arterial supply, AVM location, nidus size, and venous drainage.

RESULTS

Digital simulations demonstrated reduced image artifacts and improved TOA accuracy using radial acquisition over Cartesian. TOA mapping accuracy is more sensitive to sampling window length than time spacing. Dynamic PCASL MRA depicted seven of eight arterial pedicles, and accurately measured the AVM nidus size when the nidus was compact. The venous drainage in the AVM patients was not consistently visualized.

CONCLUSION

Dynamic 3D PCASL-VIPR with TOA mapping is able to acquire both high temporal and spatial resolution inflow dynamics that could improve diagnosis of high-flow intracranial vascular diseases.

Rapid multicomponent T2 analysis of the articular cartilage of the human knee joint at 3.0T

PURPOSE

To determine the feasibility of using multicomponent-driven equilibrium single-shot observation of T1 and T2 (mcDESPOT) for evaluating the human knee joint at 3.0T and to investigate depth-dependent and regional-dependent variations in multicomponent T2 parameters within articular cartilage.

MATERIALS AND METHODS

mcDESPOT was performed on the knee joint of 10 asymptomatic volunteers at 3.0T. Single-component T2 relaxation time (T2single ), multicomponent T2 relaxation time for water tightly bound to proteoglycan (T2PG ) and bulk water loosely bound to the macromolecular matrix (T2BW ), and fraction of water tightly bound to proteoglycan (FPG ) were measured in eight cartilage subsections and within the superficial and deep layers of patellar cartilage. Statistical analysis was used to investigate depth-dependent and regional-dependent variations in parameters.

RESULTS

There was lower (P = 0.001) T2single and T2PG and higher (P < 0.001) FPG in the deep than superficial layer of patellar cartilage. There was higher (P < 0.001) FPG on the weight-bearing surfaces than nonweight-bearing surfaces. There was higher (P < 0.001) T2single , T2PG , and T2BW on the trochlea and posterior medial and lateral femoral condyles than the patella, central medial and lateral femoral condyles, and medial and lateral tibia plateaus.

CONCLUSION

Multicomponent T2 parameters of the articular cartilage of the human knee joint can be measured at 3.0T using mcDESPOT and show depth-dependent and regional-dependent variations.

High-resolution 3D radial bSSFP with IDEAL

Radial trajectories facilitate high-resolution balanced steady state free precession (bSSFP) because the efficient gradients provide more time to extend the trajectory in k-space. A number of radial bSSFP methods that support fat-water separation have been developed; however, most of these methods require an environment with limited B0 inhomogeneity. In this work, high-resolution bSSFP with fat-water separation is achieved in more challenging B0 environments by combining a 3D radial trajectory with the IDEAL chemical species separation method. A method to maintain very high resolution within the timing constraints of bSSFP and IDEAL is described using a dual-pass pulse sequence. The sampling of a unique set of radial lines at each echo time is investigated as a means to circumvent the longer scan time that IDEAL incurs as a multiecho acquisition. The manifestation of undersampling artifacts in this trajectory and their effect on chemical species separation are investigated in comparison to the case in which each echo samples the same set of radial lines. This new bSSFP method achieves 0.63 mm isotropic resolution in a 5-min scan and is demonstrated in difficult in vivo imaging environments, including the breast and a knee with ACL reconstruction hardware at 1.5 T.

Noncontrast-enhanced three-dimensional (3D) intracranial MR angiography using pseudocontinuous arterial spin labeling and accelerated 3D radial acquisition

Pseudocontinuous arterial spin labeling (PCASL) can be used to generate noncontrast magnetic resonance angiograms of the cerebrovascular structures. Previously described PCASL-based angiography techniques were limited to two-dimensional projection images or relatively low-resolution three-dimensional (3D) imaging due to long acquisition time. This work proposes a new PCASL-based 3D magnetic resonance angiography method that uses an accelerated 3D radial acquisition technique (VIPR, spoiled gradient echo) as the readout. Benefiting from the sparsity provided by PCASL and noise-like artifacts of VIPR, this new method is able to obtain submillimeter 3D isotropic resolution and whole head coverage with a 8-min scan. Intracranial angiography feasibility studies in healthy (N = 5) and diseased (N = 5) subjects show reduced saturation artifacts in PCASL-VIPR compared with a standard time-of-flight protocol. These initial results show great promise for PCASL-VIPR for static, dynamic, and vessel selective 3D intracranial angiography.

Rapid measurement and correction of phase errors from B0 eddy currents: impact on image quality for non-Cartesian imaging

Non-Cartesian imaging sequences and navigational methods can be more sensitive to scanner imperfections that have little impact on conventional clinical sequences, an issue which has repeatedly complicated the commercialization of these techniques by frustrating transitions to multicenter evaluations. One such imperfection is phase errors caused by resonant frequency shifts from eddy currents induced in the cryostat by time-varying gradients, a phenomenon known as B(0) eddy currents. These phase errors can have a substantial impact on sequences that use ramp sampling, bipolar gradients, and readouts at varying azimuthal angles. We present a method for measuring and correcting phase errors from B(0) eddy currents and examine the results on two different scanner models. This technique yields significant improvements in image quality for high-resolution joint imaging on certain scanners. This result suggests that correcting short-time B(0) eddy currents that do not affect conventional clinical sequences may simplify the adoption of non-Cartesian methods.

Improved least squares MR image reconstruction using estimates of k-space data consistency

This study describes a new approach to reconstruct data that has been corrupted by unfavorable magnetization evolution. In this new framework, images are reconstructed in a weighted least squares fashion using all available data and a measure of consistency determined from the data itself. The reconstruction scheme optimally balances uncertainties from noise error with those from data inconsistency, is compatible with methods that model signal corruption, and may be advantageous for more accurate and precise reconstruction with many least squares-based image estimation techniques including parallel imaging and constrained reconstruction/compressed sensing applications. Performance of the several variants of the algorithm tailored for fast spin echo and self-gated respiratory gating applications was evaluated in simulations, phantom experiments, and in vivo scans. The data consistency weighting technique substantially improved image quality and reduced noise as compared to traditional reconstruction approaches.

Self-calibrated multiple-echo acquisition with radial trajectories using the conjugate gradient method (SMART-CG)

PURPOSE

To remove phase inconsistencies between multiple echoes, an algorithm using a radial acquisition to provide inherent phase and magnitude information for self correction was developed. The information also allows simultaneous support for parallel imaging for multiple coil acquisitions.

MATERIALS AND METHODS

Without a separate field map acquisition, a phase estimate from each echo in multiple echo train was generated. When using a multiple channel coil, magnitude and phase estimates from each echo provide in vivo coil sensitivities. An algorithm based on the conjugate gradient method uses these estimates to simultaneously remove phase inconsistencies between echoes, and in the case of multiple coil acquisition, simultaneously provides parallel imaging benefits. The algorithm is demonstrated on single channel, multiple channel, and undersampled data.

RESULTS

Substantial image quality improvements were demonstrated. Signal dropouts were completely removed and undersampling artifacts were well suppressed.

CONCLUSION

The suggested algorithm is able to remove phase cancellation and undersampling artifacts simultaneously and to improve image quality of multiecho radial imaging, the important technique for fast three-dimensional MRI data acquisition.

Frequency response of multipoint chemical shift-based spectral decomposition

PURPOSE

To provide a framework for characterizing the frequency response of multipoint chemical shift based species separation techniques.

MATERIALS AND METHODS

Multipoint chemical shift based species separation techniques acquire complex images at multiple echo times and perform maximum likelihood estimation to decompose signal from different species into separate images. In general, after a nonlinear process of estimating and demodulating the field map, these decomposition methods are linear transforms from the echo-time domain to the chemical-shift-frequency domain, analogous to the discrete Fourier transform (DFT). In this work we describe a technique for finding the magnitude and phase of chemical shift decomposition for input signals over a range of frequencies using numerical and experimental modeling and examine several important cases of species separation.

RESULTS

Simple expressions can be derived to describe the response to a wide variety of input signals. Agreement between numerical modeling and experimental results is very good.

CONCLUSION

Chemical shift-based species separation is linear, and therefore can be fully described by the magnitude and phase curves of the frequency response. The periodic nature of the frequency response has important implications for the robustness of various techniques for resolving ambiguities in field inhomogeneity.

Single breathhold cardiac CINE imaging with multi-echo three-dimensional hybrid radial SSFP acquisition

PURPOSE

To achieve single breathhold whole heart cardiac CINE imaging with improved spatial resolution and temporal resolution by using a multi-echo three-dimensional (3D) hybrid radial SSFP acquisition.

MATERIALS AND METHODS

Multi-echo 3D hybrid radial SSFP acquisitions were used to acquire cardiac CINE imaging within a single breathhold. An optimized interleaving scheme was developed for view ordering throughout the cardiac cycle.

RESULTS

Whole heart short axis views were acquired with a spatial resolution of 1.3 x 1.3 x 8.0 mm(3) and temporal resolution of 45 ms, within a single 17 s breathhold. The technique was validated on eight healthy volunteers by measuring the left ventricular volume throughout the cardiac cycle and comparing with the conventional 2D multiple breathhold technique. The left ventricle functional measurement bias of our proposed 3D technique from the conventional 2D technique: end diastolic volume -3.3 mL +/- 13.7 mL, end systolic volume 1.4 mL +/- 6.1 mL, and ejection fraction -1.7% +/- 4.3%, with high correlations 0.94, 0.97, and 0.91, accordingly.

CONCLUSION

A multi-echo 3D hybrid radial SSFP acquisition was developed to allow for a whole heart cardiac CINE exam in a single breathhold. Cardiac function measurements in volunteers compared favorably with the standard multiple breathhold exams.

Cartilage morphology at 3.0T: assessment of three-dimensional magnetic resonance imaging techniques

PURPOSE

To compare six new three-dimensional (3D) magnetic resonance (MR) methods for evaluating knee cartilage at 3.0T.

MATERIALS AND METHODS

We compared: fast-spin-echo cube (FSE-Cube), vastly undersampled isotropic projection reconstruction balanced steady-state free precession (VIPR-bSSFP), iterative decomposition of water and fat with echo asymmetry and least-squares estimation combined with spoiled gradient echo (IDEAL-SPGR) and gradient echo (IDEAL-GRASS), multiecho in steady-state acquisition (MENSA), and coherent oscillatory state acquisition for manipulation of image contrast (COSMIC). Five-minute sequences were performed twice on 10 healthy volunteers and once on five osteoarthritis (OA) patients. Signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) were measured from the volunteers. Images of the five volunteers and the five OA patients were ranked on tissue contrast, articular surface clarity, reformat quality, and lesion conspicuity. FSE-Cube and VIPR-bSSFP were compared to IDEAL-SPGR for cartilage volume measurements.

RESULTS

FSE-Cube had top rankings for lesion conspicuity, overall SNR, and CNR (P < 0.02). VIPR-bSSFP had top rankings in tissue contrast and articular surface clarity. VIPR and FSE-Cube tied for best in reformatting ability. FSE-Cube and VIPR-bSSFP compared favorably to IDEAL-SPGR in accuracy and precision of cartilage volume measurements.

CONCLUSION

FSE-Cube and VIPR-bSSFP produce high image quality with accurate volume measurement of knee cartilage.

POCS-enhanced correction of motion artifacts in parallel MRI

A new method for correction of MRI motion artifacts induced by corrupted k-space data, acquired by multiple receiver coils such as phased arrays, is presented. In our approach, a projections onto convex sets (POCS)-based method for reconstruction of sensitivity encoded MRI data (POCSENSE) is employed to identify corrupted k-space samples. After the erroneous data are discarded from the dataset, the artifact-free images are restored from the remaining data using coil sensitivity profiles. The error detection and data restoration are based on informational redundancy of phased-array data and may be applied to full and reduced datasets. An important advantage of the new POCS-based method is that, in addition to multicoil data redundancy, it can use a priori known properties about the imaged object for improved MR image artifact correction. The use of such information was shown to improve significantly k-space error detection and image artifact correction. The method was validated on data corrupted by simulated and real motion such as head motion and pulsatile flow.

Dual half-echo phase correction for implementation of 3D radial SSFP at 3.0 T

Fat/water separation methods such as fluctuating equilibrium magnetic resonance and linear combination steady-state free precession have not yet been successfully implemented at 3.0 T due to extreme limitations on the time available for spatial encoding with the increase in magnetic field strength. We present a method to utilize a three-dimensional radial sequence combined with linear combination steady-state free precession at 3.0 T to take advantage of the increased signal levels over 1.5 T and demonstrate high spatial resolution compared to Cartesian techniques. We exploit information from the two half-echoes within each pulse repetition time to correct the accumulated phase on a point-by-point basis, thereby fully aligning the phase of both half-echoes. The correction provides reduced sensitivity to static field (B(0)) inhomogeneity and robust fat/water separation. Resultant images in the knee joint demonstrate the necessity of such a correction, as well as the increased isotropic spatial resolution attainable at 3.0 T. Results of a clinical study comparing this sequence to conventional joint imaging sequences are included.

Characterizing and correcting gradient errors in non-cartesian imaging: Are gradient errors linear time-invariant (LTI)?

Non-Cartesian and rapid imaging sequences are more sensitive to scanner imperfections such as gradient delays and eddy currents. These imperfections vary between scanners and over time and can be a significant impediment to successful implementation and eventual adoption of non-Cartesian techniques by scanner manufacturers. Differences between the k-space trajectory desired and the trajectory actually acquired lead to misregistration and reduction in image quality. While early calibration methods required considerable scan time, more recent methods can work more quickly by making certain approximations. We examine a rapid gradient calibration procedure applied to multiecho three-dimensional projection reconstruction (3DPR) acquisitions in which the calibration runs as part of every scan. After measuring the trajectories traversed for excitations on each of the orthogonal gradient axes, trajectories for the oblique projections actually acquired during the scan are synthesized as linear combinations of these measurements. The ability to do rapid calibration depends on the assumption that gradient errors are linear and time-invariant (LTI). This work examines the validity of these assumptions and shows that the assumption of linearity is reasonable, but that gradient errors can vary over short time periods (due to changes in gradient coil temperature) and thus it is important to use calibration data matched to the scan data.

Pilot study of improved lesion characterization in breast MRI using a 3D radial balanced SSFP technique with isotropic resolution and efficient fat-water separation

PURPOSE

To assess a 3D radial balanced steady-state free precession (SSFP) technique that provides submillimeter isotropic resolution and inherently registered fat and water image volumes in comparison to conventional T2-weighted RARE imaging for lesion characterization in breast magnetic resonance imaging (MRI).

MATERIALS AND METHODS

3D projection SSFP (3DPR-SSFP) combines a dual half-echo radial k-space trajectory with a linear combination fat/water separation technique (linear combination SSFP). A pilot study was performed in 20 patients to assess fat suppression and depiction of lesion morphology using 3DPR-SSFP. For all patients fat suppression was measured for the 3DPR-SSFP image volumes and depiction of lesion morphology was compared against corresponding T2-weighted fast spin echo (FSE) datasets for 15 lesions in 11 patients.

RESULTS

The isotropic 0.63 mm resolution of the 3DPR-SSFP sequence demonstrated improved depiction of lesion morphology in comparison to FSE. The 3DPR-SSFP fat and water datasets were available in a 5-minute scan time while average fat suppression with 3DPR-SSFP was 71% across all 20 patients.

CONCLUSION

3DPR-SSFP has the potential to improve the lesion characterization information available in breast MRI, particularly in comparison to conventional FSE. A larger study is warranted to quantify the effect of 3DPR-SSFP on specificity.

3D diffusion tensor MRI with isotropic resolution using a steady-state radial acquisition

PURPOSE

To obtain diffusion tensor images (DTI) over a large image volume rapidly with 3D isotropic spatial resolution, minimal spatial distortions, and reduced motion artifacts, a diffusion-weighted steady-state 3D projection (SS 3DPR) pulse sequence was developed.

MATERIALS AND METHODS

A diffusion gradient was inserted in a SS 3DPR pulse sequence. The acquisition was synchronized to the cardiac cycle, linear phase errors were corrected along the readout direction, and each projection was weighted by measures of consistency with other data. A new iterative parallel imaging reconstruction method was also implemented for removing off-resonance and undersampling artifacts simultaneously.

RESULTS

The contrast and appearance of both the fractional anisotropy and eigenvector color maps were substantially improved after all correction techniques were applied. True 3D DTI datasets were obtained in vivo over the whole brain (240 mm field of view in all directions) with 1.87 mm isotropic spatial resolution, six diffusion encoding directions in under 19 minutes.

CONCLUSION

A true 3D DTI pulse sequence with high isotropic spatial resolution was developed for whole brain imaging in under 20 minutes. To minimize the effects of brain motion, a cardiac synchronized, multiecho, DW-SSFP pulse sequence was implemented. Motion artifacts were further reduced by a combination of linear phase correction, corrupt projection detection and rejection, sampling density reweighting, and parallel imaging reconstruction. The combination of these methods greatly improved the quality of 3D DTI in the brain.

Improved 3D phase contrast MRI with off-resonance corrected dual echo VIPR

Phase contrast (PC) magnetic resonance imaging with a three-dimensional, radially undersampled acquisition allows for the acquisition of high resolution angiograms and velocimetry in dramatically reduced scan times. However, such an acquisition is sensitive to blurring and artifacts from off-resonance and trajectory errors. A dual-echo trajectory is proposed with a novel trajectory calibration from prescan data coupled with a multi-frequency reconstruction to correct for these errors. Comparisons of phantom data and in vivo results from volunteer, and patients with arteriovenous malformations patients are presented with and without these corrections and show significant improvement of image quality when both corrections are applied. The results demonstrate significantly improved visualization of vessels, allowing for highly accelerated PC acquisitions without sacrifice in image quality.

Prolonged mechanical and physiological changes in the upper extremity following short-term simulated power hand tool use

This study investigated in-vivo changes in upper limb dynamic mechanical properties and magnetic resonance imaging (MRI) parameters following short-term power hand tool operation. Previous studies have found reduction in mechanical properties following short-term power tool usage at long build-up times. This study advances that work by having participants operate a simulated pistol grip power hand tool and evaluating changes in mechanical properties, strength, discomfort level and MRI prior to tool operation and daily for 3 d after tool operation. Twenty-four participants were randomly assigned to operate a simulated power hand tool for either a high peak reaction force of 123 N (peak torque=8 Nm, build-up time=250 ms) or at a low peak reaction force of 5 N (peak torque=2 Nm, build-up time=50 ms). Subjects operated the tool for 60 min at the rate of six times per min. A reduction in stiffness (27%, p<0.05) was observed 24 h after tool operation for the high force group and this change persisted (26%, p<0.05) up to 72 h after tool operation. Similar changes were not observed for the low force group. No changes were observed in mass moment of inertia, damping, isometric strength and damping for either group (p>0.05). There was a signal intensity increase (12%, CI 19%, 5.06%) in the supinator muscle MRI for both groups 24 h after tool operation but only the high force group remained elevated (10%, CI 13.7%, 0.06%) 72 h after tool operation. Persistent short-term changes in mechanical and MRI parameters at high force levels could indicate increased strain on the upper limb and may negatively affect ability to react during rapid forceful loading of the upper limb. This research can ultimately lead to better ergonomic interventions through quantitative power hand tool design guidelines and work practices based on understanding the damaging effects of exposure to specific levels of reaction force, build-up time and repetition, as well as providing new outcome measures for epidemiological studies.

Mechanical and magnetic resonance imaging changes following eccentric or concentric exertions

BACKGROUND

Prior work has shown that changes in mechanical parameters and magnetic resonance imaging parameters occur following submaximal eccentric activity but it is unclear whether similar changes occur following submaximal concentric activity. This study compared mechanical response parameters and MRI relaxation parameters following submaximal concentric or eccentric exertions.

METHODS

This single site, randomized study investigated in vivo changes in human upper limb dynamic mechanical properties following exposure to short term repetitive submaximal eccentric or concentric exertions. Eighteen subjects were assigned to either an eccentric or concentric group and exercised for 30 min at 50% of isometric forearm maximum voluntary contraction. Changes in strength, symptom intensity, magnetic resonance imaging T2 relaxation measurements, which are indicative of edema, and dynamic mechanical parameters (stiffness, effective mass, and damping) were ascertained prior to exercise, 1h after, and 24h later.

FINDINGS

Strength decreased following exercise (P<0.01), however only the eccentric exercise group exhibited a reduction in mechanical stiffness (55%, P<0.01) and damping (31%, P<0.05), and an increase (17%, P<0.05) in magnetic resonance imaging T2 relaxation time.

INTERPRETATION

The changes in mechanical parameters and magnetic resonance imaging findings following repetitive submaximal eccentric activity could negatively impact the ability of the arm to react to rapid forceful loading during repetitive industrial work activities and may result in increased strain on the upper limb. Similar changes were not observed following concentric exercise.

3D hyperpolarized He-3 MRI of ventilation using a multi-echo projection acquisition

A method is presented for high-resolution 3D imaging of the whole lung using inhaled hyperpolarized (HP) He-3 MR with multiple half-echo radial trajectories that can accelerate imaging through undersampling. A multiple half-echo radial trajectory can be used to reduce the level of artifact for undersampled 3D projection reconstruction (PR) imaging by increasing the amount of data acquired per unit time for HP He-3 lung imaging. The point spread functions (PSFs) for breath-held He-3 MRI using multiple half-echo trajectories were evaluated using simulations to predict the effects of T(2)* and gas diffusion on image quality. Results from PSF simulations were consistent with imaging results in volunteer studies showing improved image quality with increasing number of echoes using up to 8 half-echoes. The 8-half-echo acquisition is shown to accommodate lost breath-holds as short as 6 sec using a retrospective reconstruction at reduced resolution and also to allow reduced breath-hold time compared with an equivalent Cartesian trajectory. Furthermore, preliminary results from a 3D dynamic inhalation-exhalation maneuver are demonstrated using the 8-half-echo trajectory. Results demonstrate the first high-resolution 3D PR imaging of ventilation and respiratory dynamics in humans using HP He-3 MR.

Self-calibrated GRAPPA method for 2D and 3D radial data

A fast parallel MRI (pMRI) reconstruction method is presented for 2D and 3D radial trajectories. A limitation of the radial generalized autocalibrating partially parallel acquisitions (GRAPPA) method is the need to acquire training data prior to the actual scan. This can be eliminated by the use of self-calibration when synthesizing the missing data for each coil reconstruction. The training data for each coil are estimated by multiplying the conventionally reconstructed composite image, which contains the aliasing artifacts from undersampling, with each coil's spatial sensitivity profile. An estimate of the individual receiver spatial sensitivity profiles is obtained from the k-space data that fulfill the Nyquist sampling criterion. The frequency domain representation of the training data is then calculated at the acquired k-space sample points and at the unacquired locations at which we desire to synthesize k-space data. Fitting the acquired k-space samples to the unacquired points creates reconstruction weights that are used to synthesize unacquired radial lines. The in vivo feasibility of the method for 2D radial trajectories is illustrated with an example of 2D abdominal imaging. Preliminary results obtained after applying the method on a 3D radial steady-state free precession (SSFP) data set are also demonstrated.

Balanced SSFP imaging of the musculoskeletal system

Magnetic resonance imaging (MRI), with its unique ability to image and characterize soft tissue noninvasively, has emerged as one of the most accurate imaging methods available to diagnose bone and joint pathology. Currently, most evaluation of musculoskeletal pathology is done with two-dimensional acquisition techniques such as fast spin echo (FSE) imaging. The development of three-dimensional fast imaging methods based on balanced steady-state free precession (SSFP) shows great promise to improve MRI of the musculoskeletal system. These methods may allow acquisition of fluid sensitive isotropic data that can be reformatted into arbitrary planes for improved detection and visualization of pathology. Sensitivity to fluid and fat suppression are important issues in these techniques to improve delineation of cartilage contours, for detection of marrow edema and derangement of other joint structures.

Evaluation of ganglion cysts using vastly undersampled isotropic projection reconstruction (VIPR)

For some atypical para-articular ganglia, the presence of a joint connection is highly controversial. The proper preoperative diagnosis and identification of this joint connection for ganglion cysts is important for patient treatment and outcome. MRI is the imaging modality of choice when evaluating such lesions, but the detection of subtle joint connections remains difficult with conventional MR protocols. We investigated the utility of a steady-state free-precession acquisition with isotropic high resolution using the vastly undersampled isotropic projection reconstruction (VIPR) pulse sequence to determine if joint connections for ganglion cysts could be seen more effectively, using the knee region as a model. We evaluated four patients: two with peroneal intraneural ganglion cysts, one with adventitial cystic disease of the popliteal artery, and one patient with a more typical extraneural (intramuscular) cyst. Both conventional MR and VIPR techniques were used. In our clinical experience, we found VIPR to be superior to conventional MR techniques in detecting and depicting joint connections in typical and atypical ganglion cysts around the knee.

3D fluoroscopy with real-time 3D non-cartesian phased-array contrast-enhanced MRA

For optimized CE-MRA of the chest and abdomen, the scan time and breath-hold must be coordinated with the arrival of contrast. A 3D fluoroscopy system is demonstrated that performs real-time 3D projection reconstruction acquisition, reconstruction, and visualization using only the standard scanner hardware and operator console workstation. Unlike 2D fluorotriggering techniques, no specification of a monitoring slab or careful placement of the imaging volume is required. 3DPR data are acquired continuously throughout the examination using an eight-channel receiver and 1 s interleaved subframes. The data are reconstructed using 1 s segments for real-time monitoring with 0.8-cm isotropic spatial resolution over the entire torso, allowing full-volume axial, coronal, and sagittal MIPs to be displayed simultaneously with minimal latency. The system later uses the same scan data to generate high-spatial-resolution time-resolved sequences of the breath-hold interval. The 3D fluoroscopy system was validated on phantoms and human volunteers.