Patent Publication Number: US-9841480-B2

Title: Localization of magnetic particles by means of swift-MRI

Description:
CLAIM OF PRIORITY 
     This patent application is a nationalization under 35 U.S.C. 371, and claims the benefit of priority, of PCT/US 2010/000979, filed Apr. 1, 2010 and published as WO 2010/114609 A1, on Oct. 7, 2010, which claims the benefit of priority, under 35 U.S.C. Section 119(e), to Curtis A. Corum et al., U.S. Provisional Patent Application Ser. 61/166,108, entitled “DETECTION AND QUANTITATION OF MAGNETIC PARTICLES WITH MAGNETIC RESONANCE IMAGING,” filed on Apr. 2, 2009, which applications and publication are hereby incorporated herein by reference in their entirety. 
    
    
     STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT 
     This invention was made with government support under award number BTRR P 41 RR008079 from the National Institutes of Health (NIH). The government has certain rights in this invention. 
    
    
     BACKGROUND 
     Under certain circumstances, images generated using magnetic resonance can exhibit undesirable artifacts related to the subject under examination. For example, tissue including magnetically susceptible elements can yield an image with negative contrast. Present systems are inadequate for imaging tissue including magnetically susceptible elements. 
     OVERVIEW 
     The present subject matter enables improved detection and quantitation of a magnetic agent or a magnetic particle using MRI. In one example, a phase flattening algorithm is applied (for B1-) in order to obtain a short dead time complex MRI image, such as from Sweep Imaging with Fourier Transform (SWIFT), SWIFT with Limited TE (SWIFT-LiTE), or ultrashort TE (UTE). In one example, post-processing is used to observe the imaginary (Im) component. The Im component may include information about magnetic particles. In various examples, a region of interest (ROI) can be drawn (or otherwise specified by way of a user-interface) around the region of enhancement by a user or the ROI can be automatically determined. A fitting algorithm can be used to determine the field strength quantifying the mass of iron if the particle type is known. In addition signal enhancement (from T 1 ) can be quantified giving a separate path to quantization. 
     An example of the present subject matter can be used to detect and quantify constituent elements in a region of interest. For instance, cancer diagnosis, stem cell tracking and liver MRI can be conducted using examples of the present subject matter. 
     This overview is intended to provide an overview of subject matter of the present patent application. It is not intended to provide an exclusive or exhaustive explanation of the invention. The detailed description is included to provide further information about the present patent application. 
    
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
       In the drawings, which are not necessarily drawn to scale, like numerals may describe similar components in different views. Like numerals having different letter suffixes may represent different instances of similar components. The drawings illustrate generally, by way of example, but not by way of limitation, various embodiments discussed in the present document. 
         FIG. 1  includes a magnetic resonance system according to one example. 
         FIGS. 2A, 2B, and 2C  include diagrams for a pulse sequence for SWIFT according to one example. 
         FIGS. 3A, 3B  include images of a heart. 
         FIGS. 4A and 4B  includes images for a slice based on  FIG. 3B . 
         FIGS. 5A, 5B, and 5C  illustrate views of a heart. 
         FIG. 6  includes a method according to one example. 
     
    
    
     DETAILED DESCRIPTION 
     Part 1 includes a description of an example of a magnetic resonance system. Part 2 describes an example of the SWIFT imaging protocol. Part 3 describes an example of detecting an iron-containing agent in the myocardium. Part 4 describes iterative reconstruction for off-resonance effects in SWIFT imaging. 
     Part 1 
       FIG. 1  includes a block diagram of magnetic resonance system  100 . Magnetic resonance system  100 , in one example, depicts an imaging system  100  having magnet  105 . In one example, system  100  includes an electron paramagnetic resonance system. Magnet  105  can provide a biasing magnetic field. Coil  115  and subject  110  are positioned within the field of magnet  105 . Subject  110  can include a human body, an animal, a phantom, or other specimen. Coil  115 , sometimes referred to as an antenna, can include a transmit coil, a receive coil, a separate transmit coil and receive coil, or a transceiver coil. Coil  115  is in communication with transmitter/receiver unit  120  and with processor  130 . In various examples, coil  115  both transmits and receives radio frequency (RF) signals relative to subject  110 . Transmitter/receiver unit  120  can include a transmit/receive switch, an analog-to-digital converter (ADC), a digital-to-analog converter (DAC), an amplifier, a filter, or other modules configured to excite coil  115  and to receive a signal from coil  115 . Transmitter/receiver unit  120  is coupled to processor  330 . 
     Processor  130  can include a digital signal processor, a microprocessor, a controller, or other module. Processor  130 , in one example, is configured to generate an excitation signal (for example, a pulse sequence) for coil  115 . Processor  130 , in one example, is configured to perform a post-processing operation on the signal received from coil  115 . Processor  130  is also coupled to storage  125 , display  135  and output unit  140 . 
     Storage  125  can include a memory for storing data. The data can include image data as well as results of processing performed by processor  130 . In one example, storage  125  provides storage for executable instructions for use by processor  130 . The instructions can be configured to generate and deliver a particular pulse sequence or to implement a particular algorithm. 
     Display  135  can include a screen, a monitor, or other device to render a visible image corresponding to subject  110 . For example, display  135  can be configured to display a radial projection, a Cartesian coordinate projection, or other view corresponding to subject  110 . Output unit  140  can include a printer, a storage device, a network interface or other device configured to receive processed data. 
     The present subject matter can be implemented using system  100  or other system suitable for nuclear magnetic resonance, magnetic resonance imaging, electron spin resonance (ESR), electron paramagnetic resonance (EPR), ferromagnetic resonance (FMR), or ferromagnetic spin wave resonance (FSWR). 
     Part 2 
     In nuclear magnetic resonance (NMR, also abbreviated as magnetic resonance, MR), RF excitation can be described as sequential, simultaneous, and random. Three different corresponding NMR techniques are used, including continuous wave (CW), pulsed, and stochastic. 
     Pulsed FT spectroscopy can be used with high resolution NMR. MRI has additional technical requirements over high resolution NMR. Because the objects of interest are much larger than a test tube, inevitably the static and RF fields used in MRI are more inhomogeneous than those used in high resolution NMR. 
     As in CW, the SWIFT method uses RF sweep excitation and uses a sweep rate that exceeds the sweep rate of the CW method by more than a few orders of magnitude. Unlike the CW method in which the signal is acquired in the frequency domain, in SWIFT, the signal is considered as a time function, as in the pulsed FT method. In addition, SWIFT uses the correlation method similar to stochastic NMR in order to extract proper spectral information from the spin system response. 
     The rapid-scan FT technique and SWIFT technique have some common properties but are different in point of view to system response on excitation. Rapid-scan FT considers the system response in frequency domain and SWIFT considers the system response in the time domain. As a result, the spectra obtained using SWIFT is insensitive to the linearity of the sweep rate. This permits use of a broad class of frequency modulated pulses having more uniform excitation profiles than the chirp excitation required in rapid-scan FT. SWIFT also provides virtually simultaneous excitation and acquisition of signal. Accordingly, SWIFT has a “zero echo time”, and so is well-suited for studying objects having very fast spin-spin relaxation (or very short T 2 ). SWIFT can be used for MRI of quadrupolar nuclei, such as sodium-23, potassium-39, and boron-11. 
     SWIFT Methodology 
     SWIFT can be modeled by the method presented in  FIG. 2A . SWIFT employs a sequence of frequency-modulated pulses with short repetition time TR that exceeds the pulse length T p  by at least the amount of time needed for setting a new value (or orientation) of a magnetic field gradient used to encode spatial information. The images are processed using 3D back-projection reconstruction. In one example, frequency-modulated pulses from the hyperbolic secant family (HSn pulses) are used. In  FIG. 2B , one shaped pulse is represented which includes N different sub-pulse elements with time-dependent amplitudes and phases. During the FM pulse, an isochromat follows the effective RF field vector until the instant resonance is attained. At resonance, the isochromat is released from the RF pulse&#39;s “hug” and thereafter almost freely precesses with a small decaying modulation, yielding spectral contamination. Thus, to extract spectral information from such a spin system response, processing is performed using a cross-correlation method similar to the method of recovering phase information in stochastic NMR. The theoretically achievable signal-to-noise ratio (SNR) per unit time for SWIFT for TR&lt;&lt;T 1  is the same as that for pulsed FT. During SWIFT acquisition, the applied imaging gradients usually exceed all intrinsic gradients due to susceptibility or inhomogeneity. For this condition the images obtained are fully independent of transverse relaxation and signal intensity depends only on T 1  and spin density. The maximum T 1  contrast depends on effective flip angle and the best compromise between sensitivity and contrast will have flip angles exceeding two times the Ernst angle. If flip angles are very small, T 1  contrast is negligible, and contrast comes entirely from spin density. Other kinds of contrast can be reached by an appropriate preparation sequence prior to or interleaved with the image acquisition. 
     SWIFT provides novel and beneficial properties for MRI, including the following: 
     (a) fast: SWIFT eliminates the delays associated with refocusing pulses or gradient inversion, and also time for an excitation pulse, which is integrated with the acquisition period. As in other fast imaging sequences, SWIFT is limited by existing imaging system hardware and chosen compromise between acquisition speed, spatial resolution and SNR. 
     (b) sensitive to short T 2 : SWIFT is sensitive to excited spins having T 2 &gt;1/SW (SW=spectral width). To be specifically resolved, T 2 &gt;N/SW must be satisfied, which is theoretically feasible even for solid objects by increasing SW. 
     (c) reduced motion artifacts: Because SWIFT has no “echo time” it is less sensitive to motion artifacts. It loses less signal due to either diffusion in the presence of a gradient or uncompensated motion than other fast sequences. 
     (d) reduced dynamic range requirement: Because the different frequencies are excited sequentially the resulting signal is distributed in time with decreased amplitude of the acquired signal. This allows more effective utilization of the dynamic range of the digitizer. 
     (e) quiet: SWIFT uses a small step when changing gradients between projections, and thus, fast gradient switching that creates loud noise can be avoided. SWIFT can also be operated in rapid updated mode to reach high temporal resolution in dynamic studies. This pseudo-temporal resolution is possible because projection reconstruction, unlike Fourier imaging, samples the center of k-space with every acquisition. 
     Part 3 
     SWIFT can be used to detect superparamagnetic iron-oxide (SPIO) labeled stem cells grafted in the myocardium. 
     In addition to elements labeled using SPIO, other magnetic or paramagnetic materials can also be used as a label or contrast agent. For example, ultrasmall superparamagentic iron oxide particles USPIO, or monocrystalline iron oxide nanoparticles (MION) can also be used. The term magnetic includes the term paramagnetic. A magnetic material can be ferromagnetic (a material that produces its own persistent magnetic field), paramagnetic (a material that is attracted to a magnetic field), or diamagnetic (a material that is repulsed by a magnetic field). 
     For example, superparamagnetic iron-oxide labeling of stem cells allows their detection by MRI with high sensitivity via their T 2 * effect. Consequently, T 2 /T 2 * weighted gradient focused echo (GRE) detection leads to negative contrast (signal dropout) with low SNR and poor specificity because a similar signal void would occur at tissue boundaries or in hemorrhagic region. 
     A method can be used to convert a negative contrast to positive contrast and such methods can be applied to track SPIO labels in cardiovascular system. Since the on-resonance water signal is suppressed (or not excited), a separate anatomical image is usually necessary to register with the image containing SPIO signals. Such requirement might pose a limitation in a particular cardiovascular system, where tissue motion is likely to induce errors for co-registration. 
     Sweep imaging with Fourier transformation (SWIFT) technique can be used to image spins having extremely fast transverse relaxation rate. 
     In addition, SWIFT can be used to visualize SPIO labeled cells. 
     The present subject matter concerns SWIFT technique that can visualize the susceptibility artifacts leading to enhanced detection of magnetic particles. 
     SNR. 
     Example Methods 
     Murine embryonic stem cells can be labeled with magnetic particles and such labeled cells can be grafted directly into the myocardium of a rat. In vivo detection can be achieved using a fast GRE sequence with cardiac and respiratory gating. 
     Short and long axis cine images of the heart can be acquired at 4.7 T. Two-dimensional GRE and 3D SWIFT imaging of the heart can be performed using a 4.7 T and 9.4 T magnet system respectively. The figures illustrate in-plane resolution of 100 2  μm 2  with 500 μm slice thickness using GRE and an isotropic resolution of 98 3  μm 3  using SWIFT. For the SWIFT acquisition, bandwidth is 125 kHz, flip angle is 20°, TR is 2 ms and TE is 6 μs, and the 3D data can be acquired in 15 min with 4 averages. 
     Example Results 
       FIG. 3A  illustrates in vivo GRE detection of SPIO labeled stem cells (arrows) and  FIG. 3B  illustrates ex vivo GRE detection of SPIO labeled stem cells (arrows). The GRE images depict a short axis view of the same heart. 
       FIG. 4A  includes the magnitude image and  FIG. 4B  includes the imaginary SWIFT image corresponding approximately to the slice shown in  FIG. 3B  (note that slice thickness of ex vivo GRE is 5-fold of SWIFT). While the GRE image of  FIG. 3B  shows the typical “blooming effect” associated with susceptibility, such artifacts are suppressed by SWIFT: the hypointense region in the magnitude images is well defined with enhanced boundaries ( FIG. 4A ), which are shown on the phase images ( FIG. 4B ). Imperfections in refocusing the signal in certain boundary region (arrows) can be corrected by post-processing the data. The magnitude images provide T 1 W signals of the myocardial wall. Therefore, it is not necessary to acquire any reference image for co-registration. 
     The data from the magnetic resonance system can include complex data having a first component and a second component. For example, the complex data can include a magnitude component and a phase component or the complex data can include a real component and an imaginary component. 
     SWIFT can positive contrast imaging and is suitable for cardiovascular applications. 
     Part 4 
     The following part describes iterative reconstruction for off-resonance effects in SWIFT imaging. 
     Strong local field inhomogeneities can be challenging for MRI. The challenges are three-fold: (a) excite the larger bandwidth, (b) acquire the signal from regions with short T 2 *, and (c) determine the correct and undistorted spatial frequency of the signal. The excitation and acquisition of fast relaxing signals are well captured by the SWIFT sequence. 
     Even with high bandwidth, as used for SWIFT, the acquired signal can still be spatially distorted. Traditional methods, such a two acquisitions with different TE&#39;s, for estimating field maps, are not always possible options. This is a concern for applications involving long acquisitions that can not be repeated, in-vivo imaging of moving organs (heart, lungs etc.), or when imaging short T 2  signals, where the signal is dephased even at the shortest possible TE. For correction of field inhomogeneities when mapping is not possible, auto-focusing methods can be used for smoothly varying field inhomogenieties. 
     For 3D radial imaging the feasibility of quickly reconstructing the acquired signal to a different off-resonance frequency by applying radial phase shifts in the gridded k-space makes it tractable to use iterative methods. 
     To demonstrate that signal recovery is feasible using SWIFT signals, an example of the present subject matter can implement an iterative frequency shift method for correction of the unknown local magnetic fields which are not necessarily slowly varying. 
     Methods 
     As described elsewhere in this document, murine embryonic stem cells can be labeled. 
     After reconstruction (through gridding) to a Cartesian grid of the standard SWIFT signal, a 3D region of interest (ROI) covering the insertion area of the SPIO&#39;s can be selected for auto-focusing. In an example described herein, the volume of interest (VOI) selected encompassed 100×60×50 voxels. A series of images with different off-resonance shift can be used as a basis for estimating the correct reallocation of signal. A minimization function for minimizing the square root of the magnitude of the imaginary part of the signal can be used. For sequential refinement, decreasingly smaller VOI&#39;s of 13 3 , 9 3 , 7 3 , 5 3 , 3 3  voxels can be used for averaging to prevent noise amplification. 
     Results 
     Imaging with SWIFT confirms a locally and moderately strong T 2 *. In  FIG. 5A , the location of the ROI in a slice is identified by the rectangular box.  FIG. 5A  includes a cropped FOV of representative short axis images of in-vitro rat heart using SWIFT for 6 sequential slices, with some ballooning due to T 2 *. The box indicates the size of the VOI used for off-resonance correction. 
     The row of images in  FIG. 5B  illustrates hyperintense signal due to T 2 * at the border of the SPIO injection site for sequential slices. After correction, more homogenous signal distribution of the voids observed in the SWIFT acquisition is achieved as shown in the row of images at  FIG. 5C .  FIG. 5C  includes reconstruction of signal after iterative off-resonance correction. 
     The residual bright signal is not inconsistent with bright contrast from SPIO&#39;s when imaged with SWIFT. 
     Additional Notes 
       FIG. 6  illustrates method  600  according to one example. Method  600  includes; at  610 , labeling cells of a specimen with a magnetic particle. The magnetic particle can include a paramagnetic tag. At  620 , the method includes generating magnetic resonance (MR) image data for the specimen using swept frequency excitation and substantially simultaneous signal acquisition in a time-shared mode, the MR image data including a first component and a second component. At  630 , the method includes detecting the labeled cells using an image generated using the second component and using an anatomical image generated using the first component. 
     In one example, detecting the labeled cells using the image generated using the second component includes using at least one of a phase component or an imaginary component. One example further includes generating the anatomical image using at least one of a magnitude component or a real component of the MR image data. 
     One example further includes using a local frequency shift to correct for an off-resonance artifact during signal acquisition. Generating MR image data can include performing at least one of SWIFT, SWIFT-LiTE, or UTE. According to one example, detecting the labeled cells includes identifying a position of the labeled cells. 
     The sweep imaging with Fourier transformation (SWIFT) method can be used to visualize SPIO labeled cells grafted to tissue. SWIFT uses a frequency-swept excitation and virtually simultaneous signal acquisition in a time-shared mode. SWIFT is suited for imaging objects with fast spin-spin relaxation rates and produces data with reduced “ballooning” susceptibility artifacts. The real and imaginary parts of the SWIFT image can be used in identifying the SPIO-labeled cells. 
     The magnitude and imaginary images generated from 3D-SWIFT for heart tissue is shown in the figures along with GRE counterparts. The GRE images show “blooming effect” associated with magnetic field inhomogeneity induced by the presence of SPIO. The magnetic susceptibility artifacts are suppressed using SWIFT. The SWIFT imaginary images present a positive contrast for spins dephased by the presence of SPIO labels and the majority of the moderate T 2 * effects are not visible, thus simplifying detection of SPIO labeled cells. 
     Imperfections in refocusing signal in certain boundary region and a method utilizing iterative frequency shift can be implemented to further correct the image. The SWIFT magnitude images present anatomical details and an enhanced boundary surrounding the SPIO-associated hypointense region. The magnitude image can serve as the reference image on which the SPIO-associated signal is fused without the need to co-register the two. 
     A method according to one example can detect magnetic content in tissue including intrinsic iron. The intrinsic iron content can include ferretin, hemoglobin, or a cytochrome. 
     The above detailed description includes references to the accompanying drawings, which form a part of the detailed description. The drawings show, by way of illustration, specific embodiments in which the invention can be practiced. These embodiments are also referred to herein as “examples.” Such examples can include elements in addition to those shown and described. However, the present inventors also contemplate examples in which only those elements shown and described are provided. 
     All publications, patents, and patent documents referred to in this document are incorporated by reference herein in their entirety, as though individually incorporated by reference. In the event of inconsistent usages between this document and those documents so incorporated by reference, the usage in the incorporated reference(s) should be considered supplementary to that of this document; for irreconcilable inconsistencies, the usage in this document controls. 
     In this document, the terms “a” or “an” are used, as is common in patent documents, to include one or more than one, independent of any other instances or usages of “at least one” or “one or more.” In this document, the term “or” is used to refer to a nonexclusive or, such that “A or B” includes “A but not B,” “B but not A,” and “A and B,” unless otherwise indicated. In the appended claims, the terms “including” and “in which” are used as the plain-English equivalents of the respective terms “comprising” and “wherein.” Also, in the following claims, the terms “including” and “comprising” are open-ended, that is, a system, device, article, or process that includes elements in addition to those listed after such a term in a claim are still deemed to fall within the scope of that claim. Moreover, in the following claims, the terms “first,” “second,” and “third,” etc. are used merely as labels, and are not intended to impose numerical requirements on their objects. 
     Method examples described herein can be machine or computer-implemented at least in part. Some examples can include a computer-readable medium or machine-readable medium encoded with instructions operable to configure an electronic device to perform methods as described in the above examples. An implementation of such methods can include code, such as microcode, assembly language code, a higher-level language code, or the like. Such code can include computer readable instructions for performing various methods. The code may form portions of computer program products. Further, the code may be tangibly stored on one or more volatile or non-volatile computer-readable media during execution or at other times. These computer-readable media may include, but are not limited to, hard disks, removable magnetic disks, removable optical disks (e.g., compact disks and digital video disks), magnetic cassettes, memory cards or sticks, random access memories (RAMs), read only memories (ROMs), and the like. 
     The above description is intended to be illustrative, and not restrictive. For example, the above-described examples (or one or more aspects thereof) may be used in combination with each other. Other embodiments can be used, such as by one of ordinary skill in the art upon reviewing the above description. The Abstract is provided to comply with 37 C.F.R. §1.72(b), to allow the reader to quickly ascertain the nature of the technical disclosure. It is submitted with the understanding that it will not be used to interpret or limit the scope or meaning of the claims. Also, in the above Detailed Description, various features may be grouped together to streamline the disclosure. This should not be interpreted as intending that an unclaimed disclosed feature is essential to any claim. Rather, inventive subject matter may lie in less than all features of a particular disclosed embodiment. Thus, the following claims are hereby incorporated into the Detailed Description, with each claim standing on its own as a separate embodiment. The scope of the invention should be determined with reference to the appended claims, along with the full scope of equivalents to which such claims are entitled.