Biological magnetic resonance imaging using laser-polarized 129Xe

What can computed tomography and magnetic resonance imaging tell us about ventilation?

This review provides a summary of pulmonary functional imaging approaches for determining pulmonary ventilation, with a specific focus on multi-detector x-ray computed tomography and magnetic resonance imaging (MRI). We provide the important functional definitions of pulmonary ventilation typically used in medicine and physiology and discuss the fact that some of the imaging literature describes gas distribution abnormalities in pulmonary disease that may or may not be related to the physiological definition or clinical interpretation of ventilation. We also review the current state-of-the-field in terms of the key physiological questions yet unanswered related to ventilation and gas distribution in lung disease. Current and emerging imaging research methods are described, including their strengths and the challenges that remain to translate these methods to more wide-spread research and clinical use. We also examine how computed tomography and MRI might be used in the future to gain more insight into gas distribution and ventilation abnormalities in pulmonary disease.

Use of weekly 4DCT-based ventilation maps to quantify changes in lung function for patients undergoing radiation therapy

PURPOSE

A method has been proposed to calculate ventilation maps from four-dimensional computed tomography (4DCT) images. Weekly 4DCT data were acquired throughout the course of radiation therapy for patients with lung cancer. The purpose of our work was to use ventilation maps calculated from weekly 4DCT data to study how ventilation changed throughout radiation therapy.

METHODS

Quantitative maps representing ventilation were generated for six patients. Deformable registration was used to link corresponding lung volume elements between the inhale and exhale phases of the 4DCT dataset. Following spatial registration, corresponding Hounsfield units were input into a density-change-based model for quantifying the local ventilation. The ventilation data for all weeks were registered to the pretreatment ventilation image set. We quantitatively analyzed the data by defining regions of interest (ROIs) according to dose (V(20)) and lung lobe and by tracking the weekly ventilation of each ROI throughout treatment. The slope of the linear fit to the weekly ventilation data was used to evaluate the change in ventilation throughout treatment. A positive slope indicated an increase in ventilation, a negative slope indicated a decrease in ventilation, and a slope of 0 indicated no change. The dose ROI ventilation and slope data were used to study how ventilation changed throughout treatment as a function of dose. The lung lobe ROI ventilation data were used to study the impact of the presence of tumor on pretreatment ventilation. In addition, the lobe ROI data were used to study the impact of tumor reduction on ventilation change throughout treatment.

RESULTS

Using the dose ROI data, we found that three patients had an increase in weekly ventilation as a function of dose (slopes of 1.1, 1.4, and 1.5) and three patients had no change or a slight decrease in ventilation as a function of dose (slopes of 0.3, -0.6, -0.5). Visually, pretreatment ventilation appeared to be lower in the lobes that contained tumor. Pretreatment ventilation was 39% for lobes that contained tumor and 54% for lobes that did not contain tumor. The difference in ventilation between the two groups was statistically significant (p = 0.017). When the weekly lobe ventilation data were qualitatively observed, two distinct patterns emerged. When the tumor volume in a lobe was reduced, ventilation increased in the lobe. When the tumor volume was not reduced, the ventilation distribution did not change. The average slope of the group of lobes that contained tumors that shrank was 1.18, while the average slope of the group that did not contain tumors (or contained tumors that did not shrink) was -0.32. The slopes for the two groups were significantly different (p = 0.014).

CONCLUSIONS

We did not find a consistent pattern of ventilation change as a function of radiation dose. Pretreatment ventilation was significantly lower for lobes that contained tumor, due to occlusion of the central airway. The weekly lobe ventilation data indicated that when tumor volume shrinks, ventilation increases, and when the thoracic anatomy is not visibly changed, ventilation is likely to remain unchanged.

Synchrotron-based dynamic computed tomography of tissue motion for regional lung function measurement

During breathing, lung inflation is a dynamic process involving a balance of mechanical factors, including trans-pulmonary pressure gradients, tissue compliance and airway resistance. Current techniques lack the capacity for dynamic measurement of ventilation in vivo at sufficient spatial and temporal resolution to allow the spatio-temporal patterns of ventilation to be precisely defined. As a result, little is known of the regional dynamics of lung inflation, in either health or disease. Using fast synchrotron-based imaging (up to 60 frames s(-1)), we have combined dynamic computed tomography (CT) with cross-correlation velocimetry to measure regional time constants and expansion within the mammalian lung in vivo. Additionally, our new technique provides estimation of the airflow distribution throughout the bronchial tree during the ventilation cycle. Measurements of lung expansion and airflow in mice and rabbit pups are shown to agree with independent measures. The ability to measure lung function at a regional level will provide invaluable information for studies into normal and pathological lung dynamics, and may provide new pathways for diagnosis of regional lung diseases. Although proof-of-concept data were acquired on a synchrotron, the methodology developed potentially lends itself to clinical CT scanning and therefore offers translational research opportunities.

Multiple-exchange-time xenon polarization transfer contrast (MXTC) MRI: initial results in animals and healthy volunteers

Hyperpolarized xenon-129 is a noninvasive contrast agent for lung MRI, which upon inhalation dissolves in parenchymal structures, thus mirroring the gas-exchange process for oxygen in the lung. Multiple-exchange-time xenon polarization transfer contrast (MXTC) MRI is an implementation of the XTC MRI technique in four dimensions (three spatial dimensions plus exchange time). The aim of this study was to evaluate the sensitivity of MXTC MRI for the detection of microstructural deformations of the healthy lung in response to gravity-induced tissue compression and the degree of lung inflation. MXTC MRI was performed in four rabbits and in three healthy human volunteers. Two lung function parameters, one related to tissue- to alveolar-volume ratio and the other to average septal-wall thickness, were determined regionally. A significant gradient in MXTC MRI parameters, consistent with gravity-induced lung tissue deformation in the supine imaging position, was found at low lung volumes. At high lung volumes, parameters were generally lower and the gradient in parameter values was less pronounced. Results show that MXTC MRI permits the quantification of subtle changes in healthy lung microstructure. Further, only structures participating in gas exchange are represented in MXTC MRI data, which potentially makes the technique especially sensitive to pathological changes in lung microstructure affecting gas exchange.

Pulmonary magnetic resonance imaging for airway diseases

Pulmonary magnetic resonance (MR) imaging has been put forward as a new research and diagnostic tool mainly to overcome the limitations of computed tomography and nuclear medicine studies. However, pulmonary MR imaging has been difficult to use because of inherently low proton density, a multitude of air-tissue interfaces, which create significant magnetic field distortions and are commonly referred to as susceptibility artifacts; diminishing signal in the lung; and respiratory and/or cardiac motion artifacts. To overcome these drawbacks of pulmonary MR imaging, technical advances made during the last decade in sequencing, scanner and coil, adaptation of parallel imaging techniques, and utilization of contrast media have been reported as being useful for functional and morphologic assessment of various pulmonary diseases including airway diseases. This review article covers (1) pulmonary MR techniques for morphologic and functional assessment of airway diseases, and (2) pulmonary MR imaging for cystic fibrosis, asthma, and chronic obstructive pulmonary disease. Pulmonary MR imaging provides not only morphology-related but also pulmonary function-related information. It has the potential to replace nuclear medicine studies for the identification of regional pulmonary function and may perform a complementary role in airway disease assessment instead of nuclear medicine study. We believe that the findings of further basic studies as well as clinical applications of this new technique will validate the real significance of pulmonary MR imaging for the future of airway disease assessment and its usefulness for diagnostic radiology and pulmonary medicine.

Strongly hyperpolarized gas from parahydrogen by rational design of ligand-capped nanoparticles

The production of hyperpolarized fluids in continuous mode would broaden substantially the range of applications in chemistry, materials science, and biomedicine. Here we show that the rational design of a heterogeneous catalyst based on a judicious choice of metal type, nanoparticle size and surface decoration with appropriate ligands leads to highly efficient pairwise addition of dihydrogen across an unsaturated bond. This is demonstrated in a parahydrogen-induced polarization (PHIP) experiment by a 508-fold enhancement (±78) of a CH₃ proton signal and a corresponding 1219-fold enhancement (±187) of a CH₂ proton signal using nuclear magnetic resonance (1H-NMR). In contrast, bulk metal catalyst does not show this effect due to randomization of reacting dihydrogen. Our approach results in the largest gas-phase NMR signal enhancement by PHIP known to date. Sensitivity-enhanced NMR with this technique could be used to image microfluidic reactions in-situ, to probe nonequilibrium thermodynamics or for the study of metabolic reactions.

Principles of magnetic resonance imaging

This article aims to provide an educational document of magnetic resonance imaging principles for applied biomedical users of the technology. Basic principles are illustrated using simple experimental models on a preclinical imaging system.

Chronic obstructive pulmonary disease: safety and tolerability of hyperpolarized 129Xe MR imaging in healthy volunteers and patients

PURPOSE

To evaluate the safety and tolerability of inhaling multiple 1-L volumes of undiluted hyperpolarized xenon 129 ((129)Xe) followed by up to a 16-second breath hold and magnetic resonance (MR) imaging.

MATERIALS AND METHODS

This study was approved by the institutional review board and was HIPAA compliant. Written informed consent was obtained. Forty-four subjects (19 men, 25 women; mean age, 46.1 years ± 18.8 [standard deviation]) were enrolled, consisting of 24 healthy volunteers, 10 patients with chronic obstructive pulmonary disease (COPD), and 10 age-matched control subjects. All subjects received three or four 1-L volumes of undiluted hyperpolarized (129)Xe, followed by breath-hold MR imaging. Oxygen saturation, heart rate and rhythm, and blood pressure were continuously monitored. These parameters, along with respiratory rate and subjective symptoms, were assessed after each dose. Subjects' serum biochemistry and hematology were recorded at screening and at 24-hour follow-up. A 12-lead electrocardiogram (ECG) was obtained at these times and also within 2 hours prior to and 1 hour after (129)Xe MR imaging. Xenon-related symptoms were evaluated for relationship to subject group by using a χ(2) test and to subject age by using logistic regression. Changes in vital signs were tested for significance across subject group and time by using a repeated-measures multivariate analysis of variance test.

RESULTS

The 44 subjects tolerated all xenon inhalations, no subjects withdrew, and no serious adverse events occurred. No significant changes in vital signs (P > .27) were observed, and no subjects exhibited changes in laboratory test or ECG results at follow-up that were deemed clinically important or required intervention. Most subjects (91%) did experience transient xenon-related symptoms, most commonly dizziness (59%), paresthesia (34%), euphoria (30%), and hypoesthesia (30%). All symptoms resolved without clinical intervention in 1.6 minutes ± 0.9.

CONCLUSION

Inhalation of hyperpolarized (129)Xe is well tolerated in healthy subjects and in those with mild or moderate COPD. Subjects do experience mild, transient, xenon-related symptoms, consistent with its known anesthetic properties.

The future of ultra-high field MRI and fMRI for study of the human brain

MRI and fMRI have been used for about three and two decades respectively and much has changed over this time period, both in the quality of the data and in the range of applications for studying the brain. Apart from resolution improvements from around 4mm in the early days to below 0.5mm with modern technology, novel uses of contrast have led to the ability to sensitize images to some of the brain's structural properties at the cellular scale as well as study the localization and organization of brain function at the level of cortical columns. These developments have in part been facilitated by a continuing drive to increase the magnetic field strength. Will the next few decades see similar improvements? Here we will discuss current state of high field MRI, expected further increases in field strength, and improvements expected with these increases.

Increasing the sensitivity of magnetic resonance spectroscopy and imaging

Nuclear magnetic resonance and magnetic resonance imaging are two of the most important techniques in analytical chemistry and noninvasive medical imaging, respectively. They share a common physical basis, one aspect of which is a low intrinsic sensitivity relative to complementary techniques. Encouragingly, recent advances in physics, chemistry, engineering, and data processing have enabled significant increases in sensitivity, as measured by both increased signal-to-noise and reduced data acquisition times, allowing previously unattainable data to be acquired and also new types of experiments to be designed.

Hyperpolarized Xenon-129 gas-exchange imaging of lung microstructure: first case studies in subjects with obstructive lung disease

PURPOSE

To develop and test a method to noninvasively assess the functional lung microstructure.

MATERIALS AND METHODS

The Multiple exchange time Xenon polarization Transfer Contrast technique (MXTC) encodes xenon gas-exchange contrast at multiple delay times permitting two lung-function parameters to be derived: (i) MXTC-F, the long exchange-time depolarization value, which is proportional to the tissue to alveolar-volume ratio and (ii) MXTC-S, the square root of the xenon exchange-time constant, which characterizes thickness and composition of alveolar septa. Three healthy volunteers, one asthmatic, and two chronic obstructive pulmonary disease (COPD) (GOLD stage I and II) subjects were imaged with MXTC MRI. In a subset of subjects, hyperpolarized xenon-129 ADC MRI and CT imaging were also performed.

RESULTS

The MXTC-S parameter was found to be elevated in subjects with lung disease (P-value = 0.018). In the MXTC-F parameter map it was feasible to identify regional loss of functional tissue in a COPD patient. MXTC-F maps showed excellent regional correlation with CT and ADC (P ≥ 0.90) in one COPD subject.

CONCLUSION

The functional tissue-density parameter MXTC-F showed regional agreement with other imaging techniques. The newly developed parameter MXTC-S, which characterizes the functional thickness of alveolar septa, has potential as a novel biomarker for regional parenchymal inflammation or thickening.

A simple and low-cost device for generating hyperpolarized contrast agents using parahydrogen

A detailed description of the construction and use of a device for hyperpolarization of select contrast agents is presented. The device is based on molecular incorporation of the spin-order inherent to parahydrogen, followed by order transfer to a metastable heteronuclear alignment. Design considerations and experimental results relating to catalyst/solvent choice and handling, solvent heating, efficient gas entrainment and spin-order transfer are described. The resulting degree of hyperpolarization is shown to be substantial, ranging from a few to over 50%, depending on the choice of target molecule. Finally, the use of the hyperpolarized agent is demonstrated in a series of in vivo images.

Could 13C MRI assist clinical decision-making for patients with heart disease?

Even at this early stage of development, it is clear that the imaging of hyperpolarized (13)C-enriched molecules and their metabolic products offers a new approach to the study of the physiology and disease of the heart. The technology is practical in humans and, for this reason, we consider whether a role in clinical decision-making should motivate further development. The range of interventions available to treat coronary and valvular heart disease is already extensive, and new options are imminent. Yet the appropriate management of patients with left ventricular dysfunction can be challenging because the mechanism of reduced function may be unclear and the ability of the ventricle to respond to therapy may be difficult to predict. Pyruvate is a promising early target for development as a diagnostic agent because it lies at a critical branch point in cardiac biochemistry. The rate of metabolism of hyperpolarized pyruvate to CO(2) relative to lactate may prove to be a useful indicator of preserved mitochondrial function, and therefore provide a specific signal of viable myocardium. Other species including physiological substrates and nonphysiological molecules may provide additional information. Once suitable technology becomes available, it is likely that clinical research will progress quickly. The ability to monitor directly specific metabolic pathways may lead to an improvement in the selection of patients who will benefit from interventions, pharmacologic or otherwise.

Using hyperpolarized 3He MRI to evaluate treatment efficacy in cystic fibrosis patients

PURPOSE

To use hyperpolarized (HP) (3)He MR imaging to assess functional lung ventilation in subjects with cystic fibrosis (CF) before and after treatment.

MATERIALS AND METHODS

We performed HP (3)He static ventilation MRI scans on three subjects, using a Philips 3.0 Tesla (T) Achieva MRI scanner, before and after 11 days of in-patient treatment with combined intravenous and inhaled therapies for pulmonary exacerbations of CF. We also collected spirometry data. We quantified pulmonary ventilation volume measured with HP (3)He MRI using an advanced semi-automated analysis technique.

RESULTS

Following 11 days of treatment with intravenous antibiotics, hypertonic saline, and rhDNase, HP (3)He MR images in one subject displayed a 25% increase in total ventilation volume. Total ventilation volume in the other two subjects slightly decreased. All three subjects showed increases in FEV(1) and FVC following treatment.

CONCLUSION

In all subjects, the HP (3)He MR images provided detailed information on precisely where in the lungs gas was reaching. These data provide additional support for the conclusion that HP noble gas MRI can be a powerful tool for evaluating lung ventilation in patients with cystic fibrosis, but also raise important questions about the correlation between spirometry and HP gas MRI measurements.

Theoretical predictions of the spectroscopic parameters in noble-gas molecules: HXeOH and its complex with water

We employ state-of-the-art methods and basis sets to study the effect of inserting the Xe atom into the water molecule and the water dimer on their NMR parameters. Our aim is to obtain predictions for the future experimental investigation of novel xenon complexes by NMR spectroscopy. Properties such as molecular structure and energetics have been studied by supermolecular approaches using HF, MP2, CCSD, CCSD(T) and MP4 methods. The bonding in HXeOH···H(2)O complexes has been analyzed by Symmetry-Adapted Perturbation Theory to provide the intricate insight into the nature of the interaction. We focus on vibrational spectra, NMR shielding and spin-spin coupling constants-experimental signals that reflect the electronic structures of the compounds. The parameters have been calculated at electron-correlated and Dirac-Hartree-Fock relativistic levels. This study has elucidated that the insertion of the Xe atom greatly modifies the NMR properties, including both the electron correlation and relativistic effects, the (129)Xe shielding constants decrease in HXeOH and HXeOH···H(2)O in comparison to Xe atom; the (17)O, as a neighbour of Xe, is deshielded too. The HXeOH···H(2)O complex in its most stable form is stabilized mainly by induction and dispersion energies.

Hyperpolarized 129Xe MR imaging with balanced steady-state free precession in spontaneously breathing mouse lungs

PURPOSE

We investigated the characteristics of hyperpolarized (HP) (129)Xe magnetic resonance (MR) imaging obtained from balanced steady-state free precession (SSFP) measurement of mouse lungs, especially under spontaneous breathing, and compared the results with those obtained using traditional spoiled gradient echo (SPGR) method, focusing on improved signal-to-noise ratio (SNR) and reduced total acquisition time.

METHODS

We calculated magnetization response of the HP (129)Xe gas for the balanced SSFP sequence under spontaneous breathing to derive optimal conditions for the imaging experiment. We then placed an anesthetized mouse in the magnet (9.4T) supplied with oxygen gas and a mixture of HP (129)Xe gas supplied from a continuous-flow hyperpolarizing system. We obtained an axial plane image of the lung through balanced SSFP and SPGR sequences, changing the various magnetic resonance (MR) imaging parameters, and measured the SNR of these images.

RESULTS

We demonstrated the clear dependence of image intensity on flip angle and number of shots. The SNR was higher in balanced SSFP than in SPGR and 2.3-fold higher compared at each maximum. In contrast, total acquisition time in balanced SSFP was shortened to about one-eighth that of SPGR using a one-shot acquisition mode.

CONCLUSION

In HP (129)Xe MR imaging of the lung of a spontaneously breathing mouse, balanced SSFP sequence with multi-shot and centric order acquisition provides higher SNR in a shorter acquisition time than SPGR.

Distribution of hyperpolarized xenon in the brain following sensory stimulation: preliminary MRI findings

In hyperpolarized xenon magnetic resonance imaging (HP (129)Xe MRI), the inhaled spin-1/2 isotope of xenon gas is used to generate the MR signal. Because hyperpolarized xenon is an MR signal source with properties very different from those generated from water-protons, HP (129)Xe MRI may yield structural and functional information not detectable by conventional proton-based MRI methods. Here we demonstrate the differential distribution of HP (129)Xe in the cerebral cortex of the rat following a pain stimulus evoked in the animal's forepaw. Areas of higher HP (129)Xe signal corresponded to those areas previously demonstrated by conventional functional MRI (fMRI) methods as being activated by a forepaw pain stimulus. The percent increase in HP (129)Xe signal over baseline was 13-28%, and was detectable with a single set of pre and post stimulus images. Recent innovations in the production of highly polarized (129)Xe should make feasible the emergence of HP (129)Xe MRI as a viable adjunct method to conventional MRI for the study of brain function and disease.

Anisotropic hyperfine interactions limit the efficiency of spin-exchange optical pumping of ³He nuclei

We use accurate ab initio and quantum scattering calculations to demonstrate that the maximum ³He spin polarization that can be achieved in spin-exchange collisions with potassium (³⁹K) and silver (¹⁰⁷Ag) atoms is limited by the anisotropic hyperfine interaction. We find that spin exchange in Ag-He collisions occurs much faster than in K-He collisions over a wide range of temperatures (10-600 K). Our analysis indicates that measurements of trap loss rates of ²S atoms in the presence of cold ³He gas may be used to probe anisotropic spin-dependent interactions in atom-He collisions.

Analysis of cancer metabolism by imaging hyperpolarized nuclei: prospects for translation to clinical research

A major challenge in cancer biology is to monitor and understand cancer metabolism in vivo with the goal of improved diagnosis and perhaps therapy. Because of the complexity of biochemical pathways, tracer methods are required for detecting specific enzyme-catalyzed reactions. Stable isotopes such as (13)C or (15)N with detection by nuclear magnetic resonance provide the necessary information about tissue biochemistry, but the crucial metabolites are present in low concentration and therefore are beyond the detection threshold of traditional magnetic resonance methods. A solution is to improve sensitivity by a factor of 10,000 or more by temporarily redistributing the populations of nuclear spins in a magnetic field, a process termed hyperpolarization. Although this effect is short-lived, hyperpolarized molecules can be generated in an aqueous solution and infused in vivo where metabolism generates products that can be imaged. This discovery lifts the primary constraint on magnetic resonance imaging for monitoring metabolism-poor sensitivity-while preserving the advantage of biochemical information. The purpose of this report was to briefly summarize the known abnormalities in cancer metabolism, the value and limitations of current imaging methods for metabolism, and the principles of hyperpolarization. Recent preclinical applications are described. Hyperpolarization technology is still in its infancy, and current polarizer equipment and methods are suboptimal. Nevertheless, there are no fundamental barriers to rapid translation of this exciting technology to clinical research and perhaps clinical care.

Lung morphometry with hyperpolarized 129Xe: theoretical background

The (3) He lung morphometry technique, based on MRI measurements of hyperpolarized (3) He gas diffusion in lung airspaces, provides unique information on the lung microstructure at the alveolar level. In vivo 3D tomographic images of standard morphological parameters (airspace chord length, lung parenchyma surface-to-volume ratio, and number of alveoli per unit volume) can be generated from a rather short (several seconds) MRI scan. The technique is based on a theory of gas diffusion in lung acinar airways and experimental measurements of diffusion-attenuated MRI signal. The present work aims at developing the theoretical background of a similar technique based on hyperpolarized (129) Xe gas. As the diffusion coefficient and gyromagnetic ratio of (129) Xe gas are substantially different from those of (3) He gas, the specific details of the theory and experimental measurements with (129) Xe should be amended. We establish phenomenological relationships between acinar airway geometrical parameters and the diffusion-attenuated MR signal for human and small animal lungs, both normal lungs and lungs with mild emphysema. Optimal diffusion times are shown to be about 5 ms for human and 1.3 ms for small animals. The expected uncertainties in measuring main morphometrical parameters of the lungs are estimated in the framework of Bayesian probability theory.

Effects of pulmonary inhalation on hyperpolarized krypton-83 magnetic resonance T1 relaxation

The (83)Kr magnetic resonance (MR) relaxation time T(1) of krypton gas in contact with model surfaces was previously found to be highly sensitive to surface composition, surface-to-volume ratio, and surface temperature. The work presented here explored aspects of pulmonary (83)Kr T(1) relaxation measurements in excised lungs from healthy rats using hyperpolarized (hp) (83)Kr with approximately 4.4% spin polarization. MR spectroscopy without spatial resolution was applied to the ex vivo lungs that actively inhale hp (83)Kr through a custom designed ventilation system. Various inhalation schemes were devised to study the influence of anatomical dead space upon the measured (83)Kr T(1) relaxation times. The longitudinal (83)Kr relaxation times in the distal airways and the respiratory zones were independent of the lung inhalation volume, with T(1) = 1.3 s and T(1) = 1.0 s, depending only on the applied inhalation scheme. The obtained data were highly reproducible between different specimens. Further, the (83)Kr T(1) relaxation times in excised lungs were unaffected by the presence of up to 40% oxygen in the hp gas mixture. The results support the possible importance of (83)Kr as a biomarker for evaluating lung function.

Signal compensation and compressed sensing for magnetization-prepared MR angiography

Magnetization-prepared acquisitions offer a trade-off between image contrast and scan efficiency for magnetic resonance imaging. Because the prepared signals gradually decay, the contrast can be improved by frequently repeating the preparation, which in turn significantly increases the scan time. A common solution is to perform the data collection progressing from low- to high-spatial-frequency samples following each preparation. Unfortunately, this leads to loss of spatial resolution, and thereby image blurring. In this work, a new technique is proposed that first corrects the signal decay in high-frequency data to mitigate the resolution loss and improve the image contrast without reducing the scan efficiency. The proposed technique then employs a sparsity-based nonlinear reconstruction to further improve the image quality. In addition to reducing the amplified high-frequency noise, this reconstruction extrapolates missing k-space samples in the case of undersampled compressed-sensing acquisitions. The technique is successfully demonstrated for noncontrast-enhanced flow-independent angiography of the lower extremities, an application that substantially benefits from both the signal compensation and the nonlinear reconstruction.

Measurement of lung airways in three dimensions using hyperpolarized helium-3 MRI

Large airway measurement is clinically important in cases of airway disease and trauma. The gold standard is computed tomography (CT), which allows for airway measurement. However, the ionizing radiation dose associated with CT is a major limitation in longitudinal studies and trauma. To avoid ionizing radiation from CT, we present a method for measuring the large airway diameter in humans using hyperpolarized helium-3 (HPHe) MRI in conjunction with a dynamic 3D radial acquisition. An algorithm is introduced which utilizes the significant airway contrast for semi-automated segmentation and skeletonization which is used to derive the airway lumen diameter. The HPHe MRI method was validated with quantitative CT in an excised and desiccated porcine lung (linear regression R(2) = 0.974 and slope = 0.966 over 32 airway segments). The airway lumen diameters were then compared in 24 human subjects (22 asthmatics and 2 normals; linear regression R(2) value of 0.799 and slope = 0.768 over 309 airway segments). The feasibility for airway path analysis to areas of ventilation defect is also demonstrated.

Image-based modeling of lung structure and function

The current state-of-the-art in image-based modeling allows derivation of patient-specific models of the lung, lobes, airways, and pulmonary vascular trees. The application of traditional engineering analyses of fluid and structural mechanics to image-based subject-specific models has the potential to provide new insight into structure-function relationships in the individual via functional interpretation that complements imaging and experimental studies. Three major issues that are encountered in studies of airflow through the bronchial airways are the representation of airway geometry, the imposition of physiological boundary conditions, and the treatment of turbulence. Here we review some efforts to resolve each of these issues, with particular focus on image-based models that have been developed to simulate airflow from the mouth to the terminal bronchiole, and subjected to physiologically meaningful boundary conditions via image registration and soft-tissue mechanics models.

Imaging of lung function using hyperpolarized helium-3 magnetic resonance imaging: Review of current and emerging translational methods and applications

During the past several years there has been extensive development and application of hyperpolarized helium-3 (HP (3)He) magnetic resonance imaging (MRI) in clinical respiratory indications such as asthma, chronic obstructive pulmonary disease, cystic fibrosis, radiation-induced lung injury, and transplantation. This review focuses on the state-of-the-art of HP (3)He MRI and its application to clinical pulmonary research. This is not an overview of the physics of the method, as this topic has been covered previously. We focus here on the potential of this imaging method and its challenges in demonstrating new types of information that has the potential to influence clinical research and decision making in pulmonary medicine. Particular attention is given to functional imaging approaches related to ventilation and diffusion-weighted imaging with applications in chronic obstructive pulmonary disease, cystic fibrosis, asthma, and radiation-induced lung injury. The strengths and challenges of the application of (3)He MRI in these indications are discussed along with a comparison to established and emerging imaging techniques.

Lung perfusion measured using magnetic resonance imaging: New tools for physiological insights into the pulmonary circulation

Since the lung receives the entire cardiac output, sophisticated imaging techniques are not required in order to measure total organ perfusion. However, for many years studying lung function has required physiologists to consider the lung as a single entity: in imaging terms as a single voxel. Since imaging, and in particular functional imaging, allows the acquisition of spatial information important for studying lung function, these techniques provide considerable promise and are of great interest for pulmonary physiologists. In particular, despite the challenges of low proton density and short T2* in the lung, noncontrast MRI techniques to measure pulmonary perfusion have several advantages including high reliability and the ability to make repeated measurements under a number of physiologic conditions. This brief review focuses on the application of a particular arterial spin labeling (ASL) technique, ASL-FAIRER (flow sensitive inversion recovery with an extra radiofrequency pulse), to answer physiologic questions related to pulmonary function in health and disease. The associated measurement of regional proton density to correct for gravitational-based lung deformation (the "Slinky" effect (Slinky is a registered trademark of Pauf-Slinky incorporated)) and issues related to absolute quantification are also discussed.

Investigation of four-dimensional computed tomography-based pulmonary ventilation imaging in patients with emphysematous lung regions

A pulmonary ventilation imaging technique based on four-dimensional (4D) computed tomography (CT) has advantages over existing techniques. However, physiologically accurate 4D-CT ventilation imaging has not been achieved in patients. The purpose of this study was to evaluate 4D-CT ventilation imaging by correlating ventilation with emphysema. Emphysematous lung regions are less ventilated and can be used as surrogates for low ventilation. We tested the hypothesis: 4D-CT ventilation in emphysematous lung regions is significantly lower than in non-emphysematous regions. Four-dimensional CT ventilation images were created for 12 patients with emphysematous lung regions as observed on CT, using a total of four combinations of two deformable image registration (DIR) algorithms: surface-based (DIR(sur)) and volumetric (DIR(vol)), and two metrics: Hounsfield unit (HU) change (V(HU)) and Jacobian determinant of deformation (V(Jac)), yielding four ventilation image sets per patient. Emphysematous lung regions were detected by density masking. We tested our hypothesis using the one-tailed t-test. Visually, different DIR algorithms and metrics yielded spatially variant 4D-CT ventilation images. The mean ventilation values in emphysematous lung regions were consistently lower than in non-emphysematous regions for all the combinations of DIR algorithms and metrics. V(HU) resulted in statistically significant differences for both DIR(sur) (0.14 ± 0.14 versus 0.29 ± 0.16, p = 0.01) and DIR(vol) (0.13 ± 0.13 versus 0.27 ± 0.15, p < 0.01). However, V(Jac) resulted in non-significant differences for both DIR(sur) (0.15 ± 0.07 versus 0.17 ± 0.08, p = 0.20) and DIR(vol) (0.17 ± 0.08 versus 0.19 ± 0.09, p = 0.30). This study demonstrated the strong correlation between the HU-based 4D-CT ventilation and emphysema, which indicates the potential for HU-based 4D-CT ventilation imaging to achieve high physiologic accuracy. A further study is needed to confirm these results.

Measurement of alveolar oxygen partial pressure in the rat lung using Carr-Purcell-Meiboom-Gill spin-spin relaxation times of hyperpolarized 3He and 129Xe at 74 mT

Regional measurement of alveolar oxygen partial pressure can be obtained from the relaxation rates of hyperpolarized noble gases, (3) He and (129) Xe, in the lungs. Recently, it has been demonstrated that measurements of alveolar oxygen partial pressure can be obtained using the spin-spin relaxation rate (R(2) ) of (3) He at low magnetic field strengths (<0.1 T) in vivo. R(2) measurements can be achieved efficiently using the Carr-Purcell-Meiboom-Gill pulse sequence. In this work, alveolar oxygen partial pressure measurements based on Carr-Purcell-Meiboom-Gill R(2) values of hyperpolarized (3) He and (129) Xe in vitro and in vivo in the rat lung at low magnetic field strength (74 mT) are presented. In vitro spin-spin relaxivity constants for (3) He and (129) Xe were determined to be (5.2 ± 0.6) × 10(-6) Pa(-1) sec(-1) and (7.3 ± 0.4) × 10(-6) Pa(-1) s(-1) compared with spin-lattice relaxivity constants of (4.0 ± 0.4) × 10(-6) Pa(-1) s(-1) and (4.3 ± 1.3) × 10(-6) Pa(-1) s(-1), respectively. In vivo experimental measurements of alveolar oxygen partial pressure using (3) He in whole rat lung show good agreement (r(2) = 0.973) with predictions based on lung volumes and ventilation parameters. For (129) Xe, multicomponent relaxation was observed with one component exhibiting an increase in R(2) with decreasing alveolar oxygen partial pressure.

Interdependence of in-cell xenon density and temperature during Rb/129Xe spin-exchange optical pumping using VHG-narrowed laser diode arrays

The (129)Xe nuclear spin polarization (P(Xe)) that can be achieved via spin-exchange optical pumping (SEOP) is typically limited at high in-cell xenon densities ([Xe](cell)), due primarily to corresponding reductions in the alkali metal electron spin polarization (e.g. P(Rb)) caused by increased non-spin-conserving Rb-Xe collisions. While demonstrating the utility of volume holographic grating (VHG)-narrowed lasers for Rb/(129)Xe SEOP, we recently reported [P. Nikolaou et al., JMR 197 (2009) 249] an anomalous dependence of the observed P(Xe) on the in-cell xenon partial pressure (p(Xe)), wherein P(Xe) values were abnormally low at decreased p(Xe), peaked at moderate p(Xe) (~300 torr), and remained surprisingly elevated at relatively high p(Xe) values (>1000 torr). Using in situ low-field (129)Xe NMR, it is shown that the above effects result from an unexpected, inverse relationship between the xenon partial pressure and the optimal cell temperature (T(OPT)) for Rb/(129)Xe SEOP. This interdependence appears to result directly from changes in the efficiency of one or more components of the Rb/(129)Xe SEOP process, and can be exploited to achieve improved P(Xe) with relatively high xenon densities measured at high field (including averaged P(Xe) values of ~52%, ~31%, ~22%, and ~11% at 50, 300, 500, and 2000 torr, respectively).

Synchronous acquisition of hyperpolarised 3He and 1H MR images of the lungs - maximising mutual anatomical and functional information

The development of hybrid medical imaging scanners has allowed imaging with different detection modalities at the same time, providing different anatomical and functional information within the same physiological time course with the patient in the same position. Until now, the acquisition of proton MRI of lung anatomy and hyperpolarised gas MRI of lung function required separate breath-hold examinations, meaning that the images were not spatially registered or temporally synchronised. We demonstrate the spatially registered concurrent acquisition of lung images from two different nuclei in vivo. The temporal and spatial registration of these images is demonstrated by a high degree of mutual consistency that is impossible to achieve in separate scans and breath holds.

MRI of stroke using hyperpolarized 129Xe

Because there is no background signal from xenon in biological tissue, and because inhaled xenon is delivered to the brain by blood flow, we would expect a perfusion deficit, such as is seen in stroke, to reduce the xenon concentration in the region of the deficit. Thermal polarization yields negligible xenon signal relative to hyperpolarized xenon; therefore, hyperpolarized xenon can be used as a tracer of cerebral blood flow. Using a rat permanent right middle cerebral artery occlusion model, we demonstrated that hyperpolarized (129)Xe MRI is able to detect, in vivo, the hypoperfused area of focal cerebral ischemia, that is the ischemic core area of stroke. To the best of our knowledge, this is the first time that hyperpolarized (129)Xe MRI has been used to explore normal and abnormal cerebral perfusion. Our study shows a novel application of hyperpolarized (129)Xe MRI for imaging stroke, and further demonstrates its capacity to serve as a complementary tool to proton MRI for the study of the pathophysiology during brain hypoperfusion.

MRI of the lung: non-invasive protocols and applications to small animal models of lung disease

Magnetic resonance imaging (MRI) can be used in pre-clinical studies as a non-invasive imaging tool for assessing the morphological and functional impact of lung diseases and for evaluating the efficacy of potential treatments for airways diseases. Hyperpolarized gases ((3)He or (129)Xe) MRI provides insight into the lung ventilation function. Lung proton MRI provides information on lung diseases associated with inflammatory activity or with changes in lung tissue density. These imaging techniques can be implemented with non-invasive protocols appropriate for longitudinal investigations in small animal models of lung diseases. This chapter will detail two (3)He and proton lung MR imaging protocols applied on two models of lung pathology in rodents.

Hyperpolarized molecules in solution

Hyperpolarization is a technique to enhance the nuclear polarization and thereby increase the available signal in magnetic resonance (MR). This chapter provides an introduction to the concept of hyperpolarization as well as an overview of dynamic nuclear polarization (DNP) and para-hydrogen induced polarization (PHIP), two methods used to generate hyperpolarized molecules in aqueous solution.

Hyperpolarized noble gases as contrast agents

Hyperpolarized noble gases ((3)He and (129)Xe) can provide NMR signal enhancements of 10,000 to 100,000 times that of thermally polarized gases and have shown great potential for applications in lung magnetic resonance imaging (MRI) by greatly enhancing the sensitivity and contrast. These gases obtain a highly polarized state by employing a spin exchange optical pumping technique. In this chapter, the underlying physics of spin exchange optical pumping for production of hyperpolarized noble gases is explained and the basic components and procedures for building a polarizer are described. The storage and delivery strategies of hyperpolarized gases for in vivo imaging are discussed. Many of the problems that are likely to be encountered in practical experiments and the corresponding detailed approaches to overcome them are also discussed.

Hyperpolarized (131)Xe NMR spectroscopy

Hyperpolarized (hp) (131)Xe with up to 2.2% spin polarization (i.e., 5000-fold signal enhancement at 9.4 T) was obtained after separation from the rubidium vapor of the spin-exchange optical pumping (SEOP) process. The SEOP was applied for several minutes in a stopped-flow mode, and the fast, quadrupolar-driven T(1) relaxation of this spin I = 3/2 noble gas isotope required a rapid subsequent rubidium removal and swift transfer into the high magnetic field region for NMR detection. Because of the xenon density dependent (131)Xe quadrupolar relaxation in the gas phase, the SEOP polarization build-up exhibits an even more pronounced dependence on xenon partial pressure than that observed in (129)Xe SEOP. (131)Xe is the only stable noble gas isotope with a positive gyromagnetic ratio and shows therefore a different relative phase between hp signal and thermal signal compared to all other noble gases. The gas phase (131)Xe NMR spectrum displays a surface and magnetic field dependent quadrupolar splitting that was found to have additional gas pressure and gas composition dependence. The splitting was reduced by the presence of water vapor that presumably influences xenon-surface interactions. The hp (131)Xe spectrum shows differential line broadening, suggesting the presence of strong adsorption sites. Beyond hp (131)Xe NMR spectroscopy studies, a general equation for the high temperature, thermal spin polarization, P, for spin I ≥ 1/2 nuclei is presented.

The lung physiome: merging imaging-based measures with predictive computational models

Global measurements of the lung provided by standard pulmonary function tests do not give insight into the regional basis of lung function and lung disease. Advances in imaging methodologies, computer technologies, and subject-specific simulations are creating new opportunities to study structure-function relationships in the lung through multidisciplinary research. The digital Human Lung Atlas is an image-based resource compiled from male and female subjects spanning several decades of age. The Atlas comprises both structural and functional measures, and includes computational models derived to match individual subjects for personalized prediction of function. The computational models in the Atlas form part of the Lung Physiome project, which is an international effort to develop integrative models of lung function at all levels of biological organization. The computational models provide mechanistic interpretation of imaging measures; the Atlas provides structural data on which to base model geometry, and functional data against which to test hypotheses. The example of simulating airflow on a subject-specific basis is considered. Methods for deriving multiscale models of the airway geometry for individual subjects in the Atlas are outlined, and methods for modeling turbulent flows in the airway are reviewed.

Kinetics from indirectly detected hyperpolarized NMR spectroscopy by using spatially selective coherence transfers

An important recent development in NMR spectroscopy is the advent of ex situ dynamic nuclear polarization (DNP) approaches, which are capable of yielding liquid-state sensitivities that exceed considerably those afforded by the highest-field spectrometers. This increase in sensitivity has triggered new research avenues, particularly concerning the in vivo monitoring of metabolism and disease by NMR spectroscopy. So far such gains have mainly materialized for experiments that focus on nonprotonated, low-γ nuclei; targets favored by relatively long relaxation times T(1), which enable them to withstand the transfer from the cryogenic hyperpolarizer to the reacting centers of interest. Recent studies have also shown that transferring this hyperpolarization to protons by indirectly detected methods could successfully give rise to (1)H NMR spectra of hyperpolarized compounds with a high sensitivity. The present study demonstrates that, when merged with spatially encoded methods, indirectly detected (1)H NMR spectroscopy can also be exploited as time-resolved hyperpolarized spectroscopy. A methodology is thus introduced that can successfully deliver a series of hyperpolarized (1)H NMR spectra over a minutes-long timescale. The principles and opportunities presented by this approach are exemplified by following the in vitro phosphorylation of choline by choline kinase, a potential metabolic marker of cancer; and by tracking acetylcholine's hydrolysis by acetylcholine esterase, an important enzyme partaking in synaptic transmission and neuronal degradation.

Vertical distribution of specific ventilation in normal supine humans measured by oxygen-enhanced proton MRI

Specific ventilation (SV) is the ratio of fresh gas entering a lung region divided by its end-expiratory volume. To quantify the vertical (gravitationally dependent) gradient of SV in eight healthy supine subjects, we implemented a novel proton magnetic resonance imaging (MRI) method. Oxygen is used as a contrast agent, which in solution changes the longitudinal relaxation time (T1) in lung tissue. Thus alterations in the MR signal resulting from the regional rise in O(2) concentration following a sudden change in inspired O(2) reflect SV-lung units with higher SV reach a new equilibrium faster than those with lower SV. We acquired T1-weighted inversion recovery images of a sagittal slice of the supine right lung with a 1.5-T MRI system. Images were voluntarily respiratory gated at functional residual capacity; 20 images were acquired with the subject breathing air and 20 breathing 100% O(2), and this cycle was repeated five times. Expired tidal volume was measured simultaneously. The SV maps presented an average spatial fractal dimension of 1.13 ± 0.03. There was a vertical gradient in SV of 0.029 ± 0.012 cm(-1), with SV being highest in the dependent lung. Dividing the lung vertically into thirds showed a statistically significant difference in SV, with SV of 0.42 ± 0.14 (mean ± SD), 0.29 ± 0.10, and 0.24 ± 0.08 in the dependent, intermediate, and nondependent regions, respectively (all differences, P < 0.05). This vertical gradient in SV is consistent with the known gravitationally induced deformation of the lung resulting in greater lung expansion in the dependent lung with inspiration. This SV imaging technique can be used to quantify regional SV in the lung with proton MRI.

Xenon ventilation imaging using dual-energy computed tomography in asthmatics: initial experience

PURPOSE

To assess the feasibility of xenon ventilation computed tomography (CT) for evaluating ventilation abnormality in asthmatics.

MATERIALS AND METHODS

Twenty-two, stable asthmatics (M:F = 10:12; mean age, 57.6) were included. Single-phase, whole-thorax, dual-energy CT was performed using dual-source CT (Somatom Definition, Siemens) after subjects had inhaled 30% xenon for 90 seconds. Parameters include 512 x 512 matrix; 14 x 1.2 mm collimation; 40/140 eff. mAs at 140/80 kV; 0.45 pitch; and 0.33 seconds rotation time. On the color-coded xenon map, the extent of the ventilation defect was visually assessed using a 5-point scale in each lobe (0, absent defect; 1, 0%-25%; 2, 25%-50%; 3, 50%-75%; and 4, 75%-100%), which was defined as defect score. On the weighted average image, airway wall dimensions were measured at 4 segmental bronchi in both upper and lower lobes. Patients were classified into a defect group and a defect-free group based on the presence of defects shown on the xenon map. Pulmonary function test parameters and airway wall dimensions were compared in those 2 groups. Correlation analyses between the defect score, pulmonary function test, and airway wall dimensions were performed.

RESULTS

Sixteen asthmatics showed peripheral ventilation defects on the xenon map (defect score, 6.6 +/- 4.2). The defect group had a significantly lower forced expiratory volume in 1 second (FEV1) and thicker airway wall than that of the defect-free group (P = 0.04 and 0.02, respectively). The defect score correlated negatively with a ratio of FEV1 and forced vital capacity (FEV1/FVC) (r = -0.44, P = 0.04) and corrected diffusing capacity (r = -0.76, P = 0.04) and correlated positively with total lung capacity, functional residual volume, and residual volume (r = 0.90, P = 0.005; r= 0.99, P < 0.001; r = 0.88, P = 0.008, respectively).

CONCLUSIONS

The ventilation defects appeared on xenon ventilation CT in asthmatics with more severe airflow limitation and airway wall thickening. The extent of the ventilation defects showed correlations with parameters of pulmonary function test.

Quantitative assessment of lung microstructure in healthy mice using an MR-based 3He lung morphometry technique

The recently developed technique of lung morphometry using hyperpolarized (3)He diffusion magnetic resonance (MR) (Yablonskiy DA, Sukstanskii AL, Woods JC, Gierada DS, Quirk JD, Hogg JC, Cooper JD, Conradi MS. J Appl Physiol 107: 1258-1265, 2009) permits in vivo study of lung microstructure at the alveolar level. Originally proposed for human lungs, it also has the potential to study small animals. The technique relies on theoretical developments in the area of gas diffusion in lungs linking the diffusion attenuated MR signal to the lung microstructure. To adapt this technique to small animals, certain modifications in MR protocol and data analysis are required, reflecting the smaller size of mouse alveoli and acinar airways. This is the subject of the present paper. Herein, we established empirical relationships relating diffusion measurements to geometrical parameters of lung acinar airways with dimensions typical for mice and rats by using simulations of diffusion in the airways. We have also adjusted the MR protocol to acquire data with much shorter diffusion times compared with humans to accommodate the substantially smaller acinar airway length. We apply this technique to study mouse lungs ex vivo. Our MR-based measurements yield mean values of lung surface-to-volume ratio of 670 cm(-1), alveolar density of 3,200 per mm(3), alveolar depth of 55 μm, and mean chord length of 62 μm, all consistent with published data obtained histologically in mice by unbiased methods. The proposed technique can be used for in vivo experiments, opening a door for longitudinal studies of lung morphometry in mice and other small animals.

Ventilation from four-dimensional computed tomography: density versus Jacobian methods

Two calculation methods to produce ventilation images from four-dimensional computed tomography (4DCT) acquired without added contrast have been reported. We reported a method to obtain ventilation images using deformable image registration (DIR) and the underlying CT density information. A second method performs the ventilation image calculation from the DIR result alone, using the Jacobian determinant of the deformation field to estimate the local volume changes resulting from ventilation. For each of these two approaches, there are variations on their implementation. In this study, two implementations of the Jacobian-based methodology are evaluated, as well as a single density change-based model for calculating the physiologic specific ventilation from 4DCT. In clinical practice, (99m)Tc-labeled aerosol single photon emission computed tomography (SPECT) is the standard method used to obtain ventilation images in patients. In this study, the distributions of ventilation obtained from the CT-based ventilation image calculation methods are compared with those obtained from the clinical standard SPECT ventilation imaging. Seven patients with 4DCT imaging and standard (99m)Tc-labeled aerosol SPECT/CT ventilation imaging obtained on the same day as part of a prospective validation study were selected. The results of this work demonstrate the equivalence of the Jacobian-based methodologies for quantifying the specific ventilation on a voxel scale. Additionally, we found that both Jacobian- and density-change-based methods correlate well with global measurements of the resting tidal volume. Finally, correlation with the clinical SPECT was assessed using the Dice similarity coefficient, which showed statistically higher (p-value < 10(-4)) correlation between density-change-based specific ventilation and the clinical reference than did either Jacobian-based implementation.

Quantum rotor induced hyperpolarization

Despite its broad applicability NMR has always been limited by its inherently low sensitivity. Hyperpolarization methods have the potential to overcome this limitation and, in the case of ex situ dynamic nuclear polarization (DNP), large enhancement factors have been achieved. Although many other polarization methods have been described in the past, including chemically and parahydrogen-induced polarization and optical pumping, DNP has recently been the most popular. Here we present an additional polarization mechanism arising from quantum rotor effects in methyl groups, which generates polarizations at temperatures < 1.5 K and interferes with DNP at such temperatures. The polarization generated by this mechanism is efficiently transferred via carbon bound protons. Although quantum rotor polarizations have been studied for a small range of molecules in great detail, we observe such effects for a much broader range of substances with very different polarization rates at temperatures < 1.5 K. Moreover, we report transfer of quantum rotor polarization across a chain of protons. The observed effect not only influences the polarization in low-temperature DNP experiments but also opens a new independent avenue to generate enhanced sensitivity for NMR.

Magnetized boxes for housing polarized spins in homogeneous fields

We present novel types of permanently magnetized as well as current powered boxes built from soft-ferromagnetic materials. They provide shielded magnetic fields which are homogeneous within a large fraction of the enclosed volume, thus minimizing size, weight, and costs. For the permanently magnetized solutions, homogenization is achieved either by an optimized distribution of the permanent field sources or by jacketing the field with a soft-ferromagnetic cylindrical shell which is magnetized in parallel to the enclosed field. The latter principle may be applied up to fields of about 0.1T. With fields of about 1mT, such boxes are being used for shipping spin-polarized (3)He worldwide for MRI purposes. For current powered boxes, we present concepts and realizations of uniaxial and tri-axial shielded magnetic fields which are homogeneous on the level of 10(-4) within the entire shielded volume. This is achieved by inserting tightly fitting solenoids into a box from soft-magnetic material. The flexible tri-axial solution suits in particular laboratory applications, e.g. for establishing a spin quantization axis.