Subjects were recruited from the community by newspaper ads, flyers, and the Internet. Screening interviews and written informed consent were obtained at Boston University School of Medicine General Clinical Research Unit. Eligible subjects were randomized in permuted blocks (
n
=

4) to a 12-week intervention of either Iyengar yoga or walking for three 60-minute sessions per week, with a maximum of 36 sessions. All subjects had three MRS scans: Scan 1 at baseline; Scan 2 after their 12-week intervention; immediately after Scan 2, all subjects completed a 60-minute yoga or walking intervention, depending on group assignment, which was immediately followed by Scan 3.
Participants were 18–45 years old with no current Axis I diagnosis. Nonpsychoactive medications were allowed if the subject had been on a stable dose for at least 1 month with no anticipated changes during the study. The following items were exclusionary: any yoga practice in the previous 3 months, or a lifetime history of one yoga session/week for ≥4 weeks; current participation in psychotherapy, prayer groups, or any mind–body disciplines; a neurological disorder or medical condition that would compromise subject safety or scan data; treatment within the previous 3 months with medications that might affect the GABA system; use of tobacco products (known to affect GABA levels)
12; alcohol consumption >4 drinks/day; and contraindication to magnetic resonance evaluation.
The following instruments were used for screening: the Structured Clinical Interview for Diagnostic and Statistical Manual of Mental Disorders IV to identify Axis I Disorders
13 and the Time Line Follow Back to assess alcohol consumption.
14,15 Two (2) reliable and valid psychologic scales were selected to monitor the effects of the interventions on mood and anxiety over time. Mood was assessed with the Exercise-Induced Feeling Inventory (EIFI), which has four subscales: Positive Engagement, Revitalization, Tranquility, and Physical Exhaustion.
16 Anxiety was assessed with the State scale of the Spielberger State–Trait Anxiety Inventory (STAI).
17 The EIFI and STAI-State were given before each scan, prior to the first intervention session (week 0), and after completion of sessions at weeks 4, 8, and 12.
Metabolic equivalents (METs) are used to rate and compare the physical demands of various activities.
18 The American College of Sports Medicine list of metabolic equivalents was consulted to match the 60-minute Iyengar yoga intervention with a 60-minute walking intervention at 2.5 miles per hour (mph) on a flat surface rated at 3.0

METs. During the intervention, the Physical Activity Recall (PAR), a valid and widely used instrument, was used to convert each subject's weekly physical activity outside of the intervention into a METs score.
19,20 For each group, the mean weekly PAR METs scores was computed.
Certified Iyengar yoga instructors taught the yoga interventions, which were monitored by the Principal Investigator to ensure consistency in presentation of weekly posture sequences. Written lists of the weekly sequences and pictures of the postures were given to the subjects.
21 After 4 weeks of instruction, subjects were encouraged to practice at home. The intrascan yoga sequence was taught in class and monitored by research staff during Imaging Session II. The structure of the walking intervention was designed to be similar to that of the yoga intervention, with weekly group sessions in which subjects walked around the gym perimeter at 2.5

mph for 60 minutes. The intrascan walking session was done on a treadmill set to 2.5

mph with 0 incline. This design controlled for group effects and interaction time with research staff.
Imaging
Subjects were scanned on a 4-Tesla full-body MR scanner (Varian/UnityInova, Varian Inc., Palo Alto, CA) at Mclean Hospital in Belmont, MA. Scout images confirmed optimal positioning. After global shimming on unsuppressed water, T1-weighted anatomical images were taken in sagittal and axial planes [echo time (TE)/repetition time (TR)

=

6.2 seconds/11.4 milliseconds, field-of-view

=

22

×

22

×

8

cm (sagittal) and 22

×

22

×

16

cm (axial), readout duration

=

4

ms, receive bandwidth

=

±32

kHz, in-plane matrix size

= 128

×

256

×

16 (sagittal) and 256

×

256

×

64 (axial), in-plane resolution

=

0.94

×

1.9

mm (sagittal) and 0.94

×

0.94

mm (axial), readout points

=

512, slice-thickness

=

2.5

mm, flip-angle

=

11°].
In our previous study, a
post-hoc regional analysis that used multivoxel spectroscopic imaging showed that the greatest increase in GABA levels after the yoga intervention was in the thalamus.
8 The selection of the left thalamus was based on evidence that the left side has greater parasympathetic innervations and that GABA levels are lower in the left thalamus in post-traumatic stress disorder subjects.
22,23 For this study, an algorithm was developed to position a 2

×

2

×

3-cm voxel over the left thalamus. Proton spectroscopy implemented a MEGAPRESS [MEscher-GArwood Point-Resolved Echo Spectroscopy Sequence] difference-editing sequence specifically tuned for GABA.
24 Manual voxel shimming yielded global water-line widths ranging from 8 to 15

Hz. The MEGAPRESS sequence collected 68-millisecond echo-time spectra in an interleaved fashion where the GABA editing pulse was applied on every second transient. Additional MEGAPRESS acquisition parameters were: TR

=

2 seconds, spectral-bandwidth

=

2

kHz, readout-duration

=

512 milliseconds, Number of Excitations (NEX)

=

384, and total scan duration

=

13 minutes.
In order to quantify GABA, the difference-edited spectra were processed and then fitted with LCModel using basis sets acquired at 4

T. A separate LCModel template was used to fit the unedited 68-milisecond subspectrum to obtain creatine (Cr). All fitted metabolite areas were normalized to the fitted Cr resonance from the 68-millisecond subspectrum. One (1) spectrum from the MEGAPRESS acquisition in the thalamus was excluded from analysis due to low signal-to-noise. GABA/Cr ratios are referred to as GABA levels. In order to ascertain the gray and white matter contribution to each voxel, the axial T
1-weighted images were segmented into gray matter, white matter, and cerebrospinal fluid compartments using the commercial software package FSL 4.1 (FMRIB Software Library; Analysis Group, FMRIB; Oxford, UK).
Statistical analysis
The primary outcome variables were mood scores, anxiety scores, and thalamic GABA levels. Continuous measures were summarized by means

±

standard deviations; within-group comparisons were performed using paired
t-tests, while between-group comparisons were performed using two-sample
t-tests. Discrete measures were summarized by raw counts for numerators and denominators, as well as the associated percentages, and were compared by Fisher's exact test due to the limited sample size. Linear regression analysis was used to quantify the association between the primary outcome variables and potential predictor variables. In order to take into account within-subject correlations arising from repeated longitudinal measurements, generalized estimated equations (GEEs) were used to analyze within-group trends in mood and anxiety scores, as well as to perform between-group analyses.
25,26 All hypothesis tests were two-tailed and conducted at the α

=

0.05 significance level. Confidence intervals were two-sided and were constructed with 95% confidence. Stata 10.0 (College Station, TX) was used for analysis.