Radiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Morton, D. W.
Right arrow Articles by Winn, H. R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Morton, D. W.
Right arrow Articles by Winn, H. R.
(Radiology. 2001;218:598-601.)
© RSNA, 2001


Technical Developments

Clinically Relevant Rat Model for Testing BOLD Functional MR Imaging Techniques by Using Single-Shot Echo-planar Imaging at 1.5 T1

Douglas W. Morton, MD, PhD, Kenneth R. Maravilla, MD, Joseph R. Meno, PhD and H. Richard Winn, MD

1 From the Departments of Radiology (D.W.M., K.R.M.) and Neurological Surgery (K.R.M., J.R.M., H.R.W.), University of Washington, 1959 NE Pacific St, Box 357115, Seattle, WA 98195. From the 1999 RSNA scientific assembly. Received December 30, 1999; revision requested February 28, 2000; revision received May 15; accepted June 1. D.W.M. supported by an RSNA Research Resident Grant from the RSNA Research and Education Foundation. Work also supported by National Institutes of Health grant NS21076. Address correspondence to D.W.M.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 Materials and Methods
 Results
 Discussion
 REFERENCES
 
By using a 1.5-T whole-body magnetic resonance (MR) imager, a high-spatial-resolution single-shot echo-planar technique was developed to perform blood oxygen level dependent functional MR imaging of rat sensory cortex during forepaw stimulation. This technique produced cubic 1-mm3 voxels. Signal-to-noise ratio was 140–160 (43–44 dB). Optimal effective echo time was 50 msec. This system should prove useful for developing new functional MR imaging techniques with rapid adaptation to human use.

Index terms: Animals, 10.121412, 10.121415, 10.121416 • Brain, MR, 10.121412, 10.121415, 10.121416 • Experimental study, 10.121412, 10.121415, 10.121416 • Magnetic resonance (MR), echo planar, 10.121416 • Magnetic resonance (MR), functional imaging, 10.121412, 10.121415, 10.121416


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 Materials and Methods
 Results
 Discussion
 REFERENCES
 
Several blood oxygen level dependent (BOLD) functional magnetic resonance (MR) imaging rat models have been described (14); they generally use high-field-strength (>=2.0 T) imagers and slower gradient-echo pulse sequences. The limitations of these models include the fact that they are not using human subject echo-planar imaging techniques and that they require specialized small-bore research magnets. Human functional MR imaging is typically performed by using a 1.5-T whole-body imager with single-shot gradient-echo echo-planar pulse sequences.

To develop and test new human functional MR imaging methods, it will become increasingly important to have a small-animal model of functional MR imaging that uses techniques similar to those used in humans. Using a previously described rat model of electric stimulation of the forepaw to activate the contralateral sensory cortex (5,1), we determined parameters for BOLD functional MR imaging in a 1.5-T whole-body MR imager by using a single-shot echo-planar imaging sequence to see if a reliable rat functional MR imaging model could be developed with techniques similar to those used in human studies (6).


    Materials and Methods
 TOP
 ABSTRACT
 INTRODUCTION
 Materials and Methods
 Results
 Discussion
 REFERENCES
 
Animal Preparation
Our institutional animal care committee approved all animal procedures. We used a well-characterized rat preparation (5) and adopted the experimental paradigm for forepaw stimulation described by Hyder et al (1). Eight adult male Sprague-Dawley rats (350–450 g; B & K, Kent, Wash) were initially anesthetized with 4% halothane (Halocarbon Laboratories, River Edge, NJ) and then maintained with 1%–2% halothane and 100% oxygen. Surgical incision sites were infiltrated with subcutaneous injections of 2% mepivicaine hydrochloride (Polocaine-MPF; Astra, Westborough, Mass). The right femoral artery and vein were surgically exposed and then catheterized by using 24- and 22-gauge catheters (AngioCatheter; BD Instruments, Franklin Lakes, NJ), respectively. The femoral arterial catheter was used for continuous monitoring of arterial blood pressure and heart rate (BPA-2000 analyzer; Micro-Med, Louisville, Ky) and for measuring arterial blood gas. The femoral venous catheter was used for intravenous drug administration.

A 22-gauge intraperitoneal catheter was then placed through a small, anterior midline abdominal incision. After this, tracheotomy was performed, and the trachea was cannulated by using a 14-gauge catheter. Mechanical ventilation was administered by using a rodent respirator. The rats were then paralyzed with intravenous administration of 1 mg of D-tubocurarine chloride (Abbott Laboratories, North Chicago, Ill) per kilogram of body weight. The inspired gas mixture of oxygen and room air was adjusted to produce a fraction of inspired oxygen of 60%. The halothane anesthesia was gradually discontinued over the course of 15 minutes and replaced with an intraperitoneally administered mixture of {alpha}-chloralose (Part C8091; Sigma-Aldrich, St Louis, Mo) and urethane (Part U2500; Sigma-Aldrich) in a 1:10 ({alpha}-chloralose:urethane, by weight) ratio. The initial dose of this mixture was 30/300 mg/kg ({alpha}-chloralose/urethane).

Subcutaneous needle electrodes made from 30-gauge copper transformer wire were placed in the radial and ulnar aspects of the left forepaw. The animal was then placed prone in an MR-compatible stereotactic head holder constructed from polycarbonate (Dupont, Wilmington, Del) and acetyl resin (Delrin; Dupont) to minimize added susceptibility artifacts. Functional imaging did not begin until at least 1 hour after discontinuation of halothane anesthesia to ensure adequate anesthetic washout.

On the basis of arterial blood gas analysis, changes to ventilator settings and fraction of inspired oxygen were made to maintain pH, PaCO2, and PaO2 within the physiologic range. Body temperature was maintained by using a circulating water warming blanket (Gaymar Industries, Orchard Park, NY). Supplemental intravenous doses of D-tubocurarine chloride (1 mg/kg) were administered every 30 minutes for the duration of the experiment. Supplemental doses of an {alpha}-chloralose/urethane mixture (15/150–20/200 mg/kg) were administered intraperitoneally approximately every 60 minutes for the duration of the experiment. At the conclusion of each experiment, each animal was killed while anesthetized, according to the protocol approved by the institutional animal care committee.

Functional MR Imaging
All experiments were performed with a 1.5-T whole-body MR imager (Horizon Echo Speed; GE Medical Systems, Milwaukee, Wis) in a radio-frequency-shielded room. A 2.5-cm-diameter transmit-receive surface coil designed in our laboratory (Hayes C, unpublished data, 1997) was used for all measurements.

Functional imaging was performed by using a single-shot multisection two-dimensional standard gradient-echo echo-planar pulse sequence with the following parameters: 2,000/40, 50, or 60 (repetition time msec/echo time msec [effective]); number of signals acquired, one; flip angle, 90°; field of view, 7 x 4; matrix, 70 x 40; echo train length, 40; linear phase encode ordering; and 1-mm section thickness. The animal was placed prone in the imager with the 7-cm field of view axis and the frequency acquisition axis oriented in the left-right direction.

Both fat saturation and ramp sampling were enabled. Fat saturation forces the imager to use a narrow-band, sync-type radio-frequency pulse that produces a sharper section profile than the default spatial-spectral radio-frequency pulse, which allows submillimeter section thickness with the standard two-dimensional gradient-echo echo-planar pulse sequence. Fat saturation also eliminates a checkerboardlike artifact seen with the spatial-spectral pulse, which is thought to be related to off-axis fat excitation. With ramp sampling, an adaptive wideband sampling technique is used to allow shorter effective echo times while still collecting full echoes. With ramp sampling enabled, the maximal receiver bandwidth was 65 kHz. These parameters resulted in cubic 1-mm3 voxels.

Ten different section locations in the rat brain were imaged during each functional MR imaging measurement. One set of 10 sections was acquired during each repetition time, with individual sections acquired sequentially in an interleaved fashion—odd sections first, then even. By using this strategy, each section was stimulated once during each repetition time, while the MR imager generated a radio-frequency pulse and received an echo every repetition time divided by 10 (ie, every 200 msec).

Each functional MR imaging measurement consisted of 51 phases (ie, 51 sets of 10 sections), lasted 102 seconds, and generated 510 images. There was a rest period of at least 10 minutes between functional MR imaging measurements in animals undergoing multiple functional MR imaging measurements.

Electric forepaw stimulation was accomplished with 5-V, 0.3-msec, 3-Hz pulses administered by means of a pulse generator (SD9 Stimulator; Grass Instruments, Quincy, Mass). The maximal current delivered during each pulse was 0.75–1.00 mA. A boxcar pattern of forepaw stimulation was used as follows: off, on, off, on, and off for 11, 10, 10, 10, and 10 phases, respectively, for a total of 51 phases. The pulse generator was controlled by using a personal computer (PowerMac 8100/100; Apple Computer, Cupertino, Calif) with a commercially available software package (PSYSCOPE; New Micros, Dallas, Tex) (7) via a hardware signal interface (PSYSCOPE Button Box; New Micros). By using a 40-msec monostable multivibrator (74123; Texas Instruments, Dallas, Tex), the scope trigger output from the integrated pulse generator, or IPG, module of the MR imager was shaped into a pulse compatible with transistor-transistor logic, or TTL, which was then buffered with a 50-{Omega} line driver (74S140, Texas Instruments) and passed to an optically isolated input of the hardware signal interface. This allowed the software to count each image in the functional MR imaging measurement as it was recorded and use this information to enable and disable the pulse generator at precise predetermined times during each functional MR imaging acquisition.

Magnet Shim
We found that manually shimming the magnet, particularly along the z axis, was critical for BOLD activations to be observed. After initial imaging sequences to achieve three-axis localization within the magnet bore, an automated shim was performed with the shim volume limited to the sensory cortex. Regions near the paranasal sinuses and mastoid air cells were specifically excluded from this volume. The results of the automated shim were then used as a starting point for a manual gradient shim. We found it possible to achieve an adequate shim by adjusting only the first-order x, y, and z gradients. Once shimming had been completed, the gradient shim and the animal’s position within the magnet remained unchanged for the duration of the experiment.

Data Processing
Rapid preliminary image analysis (FUNCTOOLS; GE Medical Systems) was performed with a UNIX-based workstation during the experiment (D.W.M.) to verify the presence of BOLD activations. After completion of the experiment, images were extracted into digital imaging and communications in medicine, or DICOM, format, and a detailed analysis was performed (MEDx; Sensor Systems, Sterling, Va) with a LINUX-based computer. z scores were generated for each voxel by using an unpaired Student t test referenced to a boxcar waveform that matched the pattern of forepaw stimulation with a 4-second hemodynamic delay. Activation maps were generated by performing thresholding of these z maps, with activation defined as present in any voxel with a z score greater than or equal to 3. A cluster analysis was used to identify clusters of activated voxels. Clusters were analyzed for the number of active voxels and the peak BOLD response observed in a single activated voxel.

These parameters were entered into a spreadsheet program for final analysis (performed by D.W.M. and reviewed by K.R.M.). Activation maps were overlaid on high-spatial-resolution anatomic images by using the preliminary image analysis software package and then exported (PHOTOSHOP; Adobe Software, Seattle, Wash) for final formatting and printing.

Volume of Activated Sensory Cortex during Forepaw Stimulation
In eight different animals, we performed a single BOLD measurement to estimate the volume of sensory cortex responding to forepaw stimulation. Effective echo time was 50 msec in all cases. The number of activated voxels for each animal was used as an estimate of the activated cortex volume.

Determination of Optimum Echo Time
In three of these eight animals, we also performed BOLD measurements while varying effective echo time from 40 to 60 msec in 10-msec increments. A single BOLD measurement was made for each effective echo time in each animal. The number of activated voxels and the BOLD response amplitude in the maximally responding activated voxel were used for analysis.

Signal-to-Noise Ratio Calculations
Signal-to-noise ratio (SNR) was defined as follows: SNR = 20 · log10[RMS(signal)/RMS(noise)], where RMS is the root-mean-square of the indicated time-varying value. The root-mean-square signal intensity was estimated by calculating the root-mean-square of the signal intensity within a region of interest containing only the sensory cortex over the course of a BOLD measurement. Root-mean-square noise intensity was estimated by calculating the root-mean-square signal intensity within a region of interest containing only air over the course of a BOLD measurement.


    Results
 TOP
 ABSTRACT
 INTRODUCTION
 Materials and Methods
 Results
 Discussion
 REFERENCES
 
We consistently observed signal-to-noise ratios of 140–160 (43–44 dB) within the brain in the echo-planar BOLD images. Our functional MR imaging measurements demonstrated a highly localized area of activation centered 1–2 mm anterior to the bregma and 4 mm lateral to the midline, contralateral to the stimulated forepaw (Fig 1a). At an effective echo time of 50 msec, the activated volume was 11 mm3 ± 4 (mean ± SD; n = 8). No other areas of activation were observed in the imaged portions of the brain. The response from a single activated voxel is shown in Figure 1b. The optimal effective echo time for our BOLD activation model was 50 msec, which produced the maximal BOLD response amplitude and maximal number of activated voxels (Fig 2).



View larger version (76K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 1a. BOLD functional MR imaging response to left forepaw electric stimulation. Data from a single 1-mm section are shown. (a) Activation map shows a cluster of activated voxels in the expected location of the forepaw sensory cortex. The arrow indicates the most significantly activated voxel (highest score). (b) Graph shows time course of the BOLD response in the most significantly activated voxel in a. Horizontal bars at the top of b indicate 20-second periods of left forepaw electric stimulation.

 


View larger version (20K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 1b. BOLD functional MR imaging response to left forepaw electric stimulation. Data from a single 1-mm section are shown. (a) Activation map shows a cluster of activated voxels in the expected location of the forepaw sensory cortex. The arrow indicates the most significantly activated voxel (highest score). (b) Graph shows time course of the BOLD response in the most significantly activated voxel in a. Horizontal bars at the top of b indicate 20-second periods of left forepaw electric stimulation.

 


View larger version (30K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 2a. BOLD response as a function of effective echo time (TE; n = 3). (a) Graph shows the number of activated voxels is maximal for an effective echo time near 50 msec (P < .05, Student t test). Error bars denote the SD. (b) Graph shows the BOLD response amplitude of the most significantly activated voxel is maximal for an effective echo time near 50 msec (P < .05, Student t test). Error bars denote the SD.

 


View larger version (32K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 2b. BOLD response as a function of effective echo time (TE; n = 3). (a) Graph shows the number of activated voxels is maximal for an effective echo time near 50 msec (P < .05, Student t test). Error bars denote the SD. (b) Graph shows the BOLD response amplitude of the most significantly activated voxel is maximal for an effective echo time near 50 msec (P < .05, Student t test). Error bars denote the SD.

 
The Table documents all physiologic parameters recorded for all animals immediately prior to the first BOLD measurement. Mean arterial blood pressure during each BOLD measurement never changed by more than 2 mm Hg from its immediate premeasurement value.


View this table:
[in this window]
[in a new window]

 
Physiologic Parameters for the Eight Animals in This Study
 

    Discussion
 TOP
 ABSTRACT
 INTRODUCTION
 Materials and Methods
 Results
 Discussion
 REFERENCES
 
In this study, we describe a rat model for BOLD functional MR imaging by using techniques similar to those used with human subjects but with spatial resolution and section thickness optimized for small-animal studies. BOLD functional MR imaging measurements in which echo-planar imaging techniques with a 2-T magnet were used have been demonstrated recently in rats (2). However, to our knowledge, this is the first report of small-animal functional MR imaging performed with a conventional 1.5-T whole-body MR imager by using echo-planar imaging techniques.

We observed a localized focus of activation, centered 1–2 mm anterior and 4 mm lateral to the bregma, in response to forepaw stimulation, which is consistent with the known location of the forepaw somatosensory cortex (8,9). We determined the activation volume to be 11 mm3 ± 4, which is similar to the findings of Hyder et al (1) and Gyngell et al (3) but somewhat larger than that observed by Hsu et al (4). Differences between the activation volumes reported by these groups are likely to represent differences in a combination of factors, including forepaw electrode placement, stimulation paradigm, imaging pulse sequence, animal weight and age, and statistical threshold for activation.

We observed maximal BOLD activation with a 4-second hemodynamic delay incorporated into our boxcar reference waveform. This 4-second delay is consistent with the temporal course of deoxyhemoglobin levels in response to visual stimulation reported by Malonek et al (10). Given the more rapid time course of sensory evoked potentials (approximately 50 msec) and sensory cortex arteriolar vasodilatation (approximately 1.5 seconds) in response to peripheral nerve stimulation (5,11), our results are consistent with the hypothesis that the BOLD response is primarily related to changes in deoxyhemoglobin levels associated with neural activation.

The primary advantage of imaging rats with a whole-body 1.5-T MR imager and a standard human echo-planar imaging sequence is the potential to use this system as a tool for testing and developing human functional MR imaging techniques. New pulse sequence parameters or unique pulse sequence designs could first be worked out in this small-animal system and then rapidly applied to human studies.

An additional advantage of the described technique is that it permits small-animal functional MR imaging research studies to be conducted on a clinical 1.5-T MR platform with only minor hardware and software modifications. These experiments were previously performed only with small-bore high-field-strength systems.

Signal-to-noise ratio is the primary limitation of small-animal functional MR imaging at 1.5 T. With careful coil design and pulse sequence parameter selection, we were able to achieve a high signal-to-noise ratio that permitted resolution of cubic 1-mm3 voxels within the sensory cortex. This resolution should be sufficient to observe activations in many rat brain structures, such as the visual cortex, forepaw and hind paw sensory cortex, or olfactory cortex.

Manual shimming is essential to achieve good results. Whole-body imagers are optimized to image large structures by using a large field of view. Rat functional MR imaging measurements require a small field of view and a small-geometry radio-frequency coil. However, the small field of view and small-geometry radio-frequency coil provide insufficient signal for a traditional automated shim. We found it sufficient to adjust only the first-order shims while optimizing the height and width of the water resonance peak.

In conclusion, it is possible to develop a robust small-animal model of BOLD functional MR imaging by using single-shot echo-planar imaging techniques with a 1.5-T whole-body imager. Since this model uses techniques commonly used for human functional MR imaging studies, it should prove useful for testing and developing new human functional MR imaging techniques, as well as for enabling basic research into the neurovascular physiology underlying the BOLD response by using a 1.5-T clinical MR system.


    FOOTNOTES
 
Abbreviation: BOLD = blood oxygen level dependent

Author contributions: Guarantor of integrity of entire study, D.W.M.; study concepts, all authors; study design, D.W.M., K.R.M.; definition of intellectual content, D.W.M., K.R.M.; literature research, D.W.M., K.R.M.; experimental studies, D.W.M.; data acquisition and analysis, D.W.M.; statistical analysis, D.W.M.; manuscript preparation and editing, D.W.M.; manuscript review, all authors; manuscript final version approval, D.W.M., K.R.M.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 Materials and Methods
 Results
 Discussion
 REFERENCES
 

  1. Hyder F, Behar KL, Martin MA, Blamire AM, Shulman RG. Dynamic magnetic resonance imaging of the rat brain during forepaw stimulation. J Cereb Blood Flow Metab 1994; 14:649-655.[Medline]
  2. Mandeville JB, Marota JJ, Kosofsky BE, et al. Dynamic functional imaging of relative cerebral blood volume during rat forepaw stimulation. Magn Reson Med 1998; 39:615-624.[Medline]
  3. Gyngell ML, Bock C, Schmitz B, Hoehn-Berlage M, Hossmann KA. Variation of functional MRI signal in response to frequency of somatosensory stimulation in alpha-chloralose anesthetized rats. Magn Reson Med 1996; 36:13-15.[Medline]
  4. Hsu EW, Hedlund LW, MacFall JR. Functional MRI of the rat somatosensory cortex: effects of hyperventilation. Magn Reson Med 1998; 40:421-426.[Medline]
  5. Morii S, Ngai AC, Winn HR. Reactivity of rat pial arterioles and venules to adenosine and carbon dioxide: with detailed description of the closed cranial window technique in rats. J Cereb Blood Flow Metab 1986; 6:34-41.[Medline]
  6. Morton DW, Maravilla KR, Meno JR, Winn HR. A clinically relevant rat model for testing BOLD fMRI techniques using single-shot EPI at 1.5T (abstr). Radiology 1999; 213(P):395.[Abstract/Free Full Text]
  7. Cohen JD, MacWhinney B, Flatt M, Provost J. PsyScope: a new graphic interactive environment for designing psychology experiments. Behav Res Meth Instruments Comput 1993; 25:257-271.
  8. Chapin JK, Lin CS. Mapping the body representation in the SI cortex of anesthetized and awake rats. J Comp Neurol 1984; 229:199-213.[Medline]
  9. Riddle DR, Purves D. Individual variation and lateral asymmetry of the rat primary somatosensory cortex. J Neurosci 1995; 15:4184-4195.[Abstract]
  10. Malonek D, Dirnagl U, Lindauer U, Yamada K, Kanno I, Grinvald A. Vascular imprints of neuronal activity: relationships between the dynamics of cortical blood flow, oxygenation, and volume changes following sensory stimulation. Proc Natl Acad Sci U S A 1997; 94:14826-14831.[Abstract/Free Full Text]
  11. Ko KR, Ngai AC, Winn HR. Role of adenosine in regulation of regional cerebral blood flow in sensory cortex. Am J Physiol 1990; 259(6 pt 2):H1703-H1708.[Abstract/Free Full Text]



This article has been cited by other articles:


Home page
Am. J. Neuroradiol.Home page
D.W. Morton, B. Keogh, K. Lim, and K.R. Maravilla
Functional Brain Imaging Using a Long Intravenous Half-Life Gadolinium-Based Contrast Agent
AJNR Am. J. Neuroradiol., August 1, 2006; 27(7): 1467 - 1471.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Neuroradiol.Home page
D. W. Morton, K. R. Maravilla, J. R. Meno, and H. R. Winn
Systemic Theophylline Augments the Blood Oxygen Level--Dependent Response to Forepaw Stimulation in Rats
AJNR Am. J. Neuroradiol., April 1, 2002; 23(4): 588 - 593.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Morton, D. W.
Right arrow Articles by Winn, H. R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Morton, D. W.
Right arrow Articles by Winn, H. R.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
RADIOLOGY RADIOGRAPHICS RSNA JOURNALS ONLINE