the presence and modulation of neuronal responses within specific time and frequency bands has been implicated in the implementation of perception and cognition within the brain .
there is also increasing evidence that the most plausible mechanism for the generation of temporally organized network activity is in reciprocally connected neuronal networks containing mixtures of interconnected glutamatergic ( excitatory ) pyramidal and stellate cells , and gabaergic ( inhibitory ) interneurons ( bartos et al , 2007 ; traub et al , 1996 ) .
however , the neuropharmacological bases of the spatially summated , population - level , neurophysiological responses ( both evoked and induced ) that are recorded non - invasively in humans with magnetoencephalographic / electroencephalographic ( meg / eeg ) are largely unknown . in a recent observational study we found that , across individuals , the frequency of stimulus - induced network gamma oscillations in primary visual cortex was positively correlated with the concentration of gaba measured with edited magnetic resonance spectroscopy ( mrs ) ( muthukumaraswamy et al , 2009 ) .
however , mrs is an indirect measure of synaptic gaba function and could be influenced by a number of anatomical , biochemical / metabolic , or even genetic factors . in this experiment
, we wished to test the hypothesis that modulation of endogenous gaba levels would modify either the amplitude or frequency of stimulus - related gamma oscillations in human visual cortex . to do this
, we manipulated endogenous gaba levels using the drug tiagabine in a set of healthy control participants in a single - blind crossover design .
tiagabine , used in the treatment of epilepsy , binds with high affinity and selectivity to gaba transporter 1 ( gat-1 ; borden et al , 1994 ) , the primary gaba reuptake transporter in the human cerebral cortex ( conti et al , 2004 ) .
gat-1 transporters are located on both neurons and glia ( minelli et al , 1995 ) and their presence at fast - inhibitory gabaa synapses terminates gaba activity and shapes ipsp activity .
the effect of tiagabine , measurable by microdialysis , is to elevate the extracellular / synaptic concentration of gaba ( dalby , 2000 ; fink - jensen et al , 1992 ) and therefore to enhance the actions of endogenous gaba .
the electrophysiological effect of tiagabine is to enhance evoked inhibitory postsynaptic current ( ipsc ) duration without modifying ipsc amplitude ( roepstorff and lambert , 1994 ; thompson and gahwiler , 1992 ) . to directly measure neuronal activity in humans , we combined meg recordings with source estimations using an experimental paradigm that we have previously demonstrated to be highly robust to potentially confounding factors such as day , session , or repetition effects ( muthukumaraswamy et al , 2010 )
. meg is an ideal technique to capture higher frequency oscillations as it is more robust to high - frequency artefact sources such as saccades and neck muscle artefacts than eeg ( fries et al , 2008 ; whitham et al , 2007 , 2008 ; yuval - greenberg et al , 2008 , 2009 ) .
the principal hypothesis we tested was that modulation of endogenous gaba levels would modify either the amplitude or frequency of stimulus - related gamma oscillations .
subsidiary to that hypothesis , we tested for gabaergic modulations of several neurophysiological parameters , including stimulus - induced alpha amplitude decreases and evoked response components , which are also known to occur with these stimuli ( muthukumaraswamy et al , 2010 ) .
eighteen right - handed volunteers ( 14 men and 4 women ) participated in the experiment after giving informed consent , with all procedures approved by the uk national research ethics service ( south east wales ) .
volunteers were screened and excluded for personal histories of neurological and psychiatric disease , current recreational or prescription drug use , and impaired liver function ( by standard liver function tests ) .
early in the data collection , three women volunteers were not able to complete the experiment after becoming too heavily sedated . from this point in the data collection
, women were excluded from recruitment , leaving completed data sets from 14 men and 1 woman to analyze ( mean age 25.5 , range 2032 ) .
all results and conclusions are similar whether this one woman was included / excluded from the analysis .
participants were scanned on 2 separate days , separated by a minimum 7-day washout period , at approximately the same time - of - day . on each day , an initial  pre ' meg recording was obtained then participants orally ingested a capsule containing either placebo or 15  mg of tiagabine ( gabitril ) .
this oral dose is similar to the 16  mg used in a pet study by frankle et al ( 2009 ) and identical to several ongoing imaging studies by our group .
participants were blinded to the contents of the capsule and placebo / control session order was counterbalanced across participants . following drug / placebo administration , meg recordings were obtained from participants at 1 , 3 , and 5-h time - points . at the conclusion of every meg recording
, participants completed several psychological questionnaires including , the biphasic alcohol effects scale ( martin et al , 1993 ) , the subjective high assessment scale ( schuckit , 1980 ) and asked to qualitatively report their psychological state . during each meg recording , participants sat in a magnetically shielded room and were presented with a visual stimulus consisting of a , vertical , stationary , maximum contrast , three cycles per degree , square - wave grating presented on a mean luminance background .
the stimulus was presented in the lower left visual field and subtended 8 both horizontally and vertically .
a small red fixation square was located at the top right - hand edge of the stimulus and remained on for the entire stimulation protocol ( muthukumaraswamy , 2010 ; swettenham et al , 2009 ) .
the stimulus was presented on a mitsubishi diamond pro 2070 monitor controlled by the psychophysics toolbox ( brainard , 1997 ; pelli , 1997 ) .
the screen size was 1024 by 768 pixels and the monitor frame rate was 100  hz . the monitor was outside the magnetically shielded room and viewed directly from within , at 2.15  m , through a cut - away portal in the shield .
the duration of each stimulus was 11.5  s followed by 1.5  s of fixation square only .
participants were instructed to fixate for the entire experiment and , in order to maintain attention , were instructed to 
responses slower than 750  ms from the termination of the stimulus triggered a visual warning for participants . in each recording session , 120 stimuli were presented and participants responded with their right hand .
whole - head meg recordings were made using a ctf 275-channel radial gradiometer system sampled at 1200  hz ( 0300  hz bandpass ) .
an additional 29 reference channels were recorded for noise cancellation purposes and the primary sensors were analyzed as synthetic third - order gradiometers ( vrba and robinson , 2001 ) .
three of the 275 channels were turned off because of excessive sensor noise . at the onset of each stimulus presentation
participants were fitted with three electromagnetic head coils ( nasion and pre - auriculars ) , which they wore for the entire recording day , and were localized relative to the meg system immediately before and after the recording session .
each participant had a 1  mm isotropic mri scan available for source localization analysis . to achieve mri / meg co - registration ,
the fiduciary markers were placed at fixed distances from anatomical landmarks identifiable in participants ' anatomical mris ( tragus , eye centre ) .
offline , data were first epoched from 1 to 1  s around stimulus onset and each trial visually inspected for data quality .
data with gross artefacts , such as , head movements and muscle clenching were excluded from further analysis .
reaction times were calculated for each trial relative to the ( variable ) offset of each grating stimulus .
any trial where there was no behavioural responses or a response > 750  ms were classified as incorrect responses and excluded from further analysis . for each participant , we ensured that equal number of trials , from each recording was used in the analysis . to achieve this
, trials were removed from the end of each recording run so that each participant had an equal number of trials submitted for each condition in the analysis . following these trial removal procedures , the mean number of trials remaining to analyze was 105.5 per participant ( range 82117 ) .
three source localizations were performed on each data set using synthetic aperture magnetometry ( robinson and vrba , 1999 ) , one for gamma ( 3080  hz ) , one for alpha ( 815  hz ) , and one for evoked responses ( samerf ; robinson , 2004 ) .
correspondingly , three global covariance matrices were calculated for each data set , one for gamma , one for alpha , and one for samerf ( 0100  hz ) .
based on these covariance matrices , using the beamformer algorithm ( robinson and vrba , 1999 ) , three sets of beamformer weights were computed for the entire brain at 4  mm isotropic voxel resolution . a multiple local - spheres
( huang et al , 1999 ) volume conductor model was derived by fitting spheres to the brain surface extracted by fsl 's brain extraction tool ( smith , 2002 ) .
readers unfamiliar with sam are referred to several methodological reviews ( hillebrand et al , 2005 ; vrba and robinson , 2001 ) . for gamma - band sam image reconstruction virtual sensors
were constructed for each beamformer voxel and student 's t - images of source power changes computed using a baseline period of 1 to 0  s and an active period of 0 to 1  s , while for alpha - band sam image reconstruction a baseline period of 1 to 0.3  s and an active period of 0.3 to 1  s was used .
this latter time window was used to avoid confounding alpha responses with the time - locked evoked response . after generating the source reconstruction images ,
the voxel with the strongest power increase ( in the contralateral occipital lobe ) was located for gamma and the voxel with the strongest source power decrease located for alpha .
the source analysis and peak voxel identification procedures were carried out separately for each recording as participants were moved in and out of the meg between recordings , which would alter their head position in the meg dewar . to reveal time - frequency responses at these locations , the virtual sensors at these locations
were repeatedly band - passed filtered between 1 and 160  hz at 0.5  hz frequency step intervals , filtering with an 8  hz bandpass , third - order butterworth filter ( le van quyen et al , 2001 ; muthukumaraswamy et al , 2010 ) .
the hilbert transform was used to obtain the amplitude envelope and time - frequency spectra were computed as a percentage changed from the pre - stimulus baselines for each frequency band . from these spectra ,
the peak frequency and amplitude at that frequency for late gamma amplitude increases ( 3080  hz , 0.31 
s ) were quantified by collapsing across the time dimension and finding the maximal ( or minimal for alpha ) frequency .
we also checked for confounding effects of using relative change baselines by computing the average absolute values of the source amplitude estimates within the respective time - frequency bins .
for samerf , the computed evoked response was passed through the 0100  hz beamformer weights to generate samerf images ( robinson , 2004 ) at 10  ms intervals from 50150  ms .
the image window ( usually 7080  ms or 8090  ms ) with the maximal response in visual cortex was identified and the maximal voxel selected as the peak location for further analysis . for this , the evoked field was computed for this virtual sensor ( 0.2 to 0-s baseline ) and the peak amplitude and latency of the m80 and m150 responses quantified . for the production of grand - average sam maps , individual sam images were first spatially normalized onto the mni ( t1 ) average brain using fmrib 's linear affine registration tool ( jenkinson and smith , 2001 ) .
this was done by first obtaining a set of warping parameters by registering the participant 's anatomical mri with the average brain and then applying these parameters to the sam source power maps .
the grand - average images presented in figure 2 are the grand - average across both participants and conditions .
for these analyses , there were main effects of  day ' ,  session ' , and an interaction term  drug ' , presented in that order . here
, the critically important term is the interaction term , indicative of a tiagabine drug - effect .
it should be noted that for f(1 , 14 ) the critical value is 4.6 and for f(3 , 42 ) the critical value is 2.82 ( alpha=0.05 ) .
eighteen right - handed volunteers ( 14 men and 4 women ) participated in the experiment after giving informed consent , with all procedures approved by the uk national research ethics service ( south east wales ) .
volunteers were screened and excluded for personal histories of neurological and psychiatric disease , current recreational or prescription drug use , and impaired liver function ( by standard liver function tests ) .
early in the data collection , three women volunteers were not able to complete the experiment after becoming too heavily sedated . from this point in the data collection
, women were excluded from recruitment , leaving completed data sets from 14 men and 1 woman to analyze ( mean age 25.5 , range 2032 ) .
all results and conclusions are similar whether this one woman was included / excluded from the analysis .
participants were scanned on 2 separate days , separated by a minimum 7-day washout period , at approximately the same time - of - day . on each day , an initial  pre ' meg recording was obtained then participants orally ingested a capsule containing either placebo or 15  mg of tiagabine ( gabitril ) .
this oral dose is similar to the 16  mg used in a pet study by frankle et al ( 2009 ) and identical to several ongoing imaging studies by our group .
participants were blinded to the contents of the capsule and placebo / control session order was counterbalanced across participants . following drug / placebo administration , meg recordings were obtained from participants at 1 , 3 , and 5-h time - points . at the conclusion of every meg recording
, participants completed several psychological questionnaires including , the biphasic alcohol effects scale ( martin et al , 1993 ) , the subjective high assessment scale ( schuckit , 1980 ) and asked to qualitatively report their psychological state . during each meg recording , participants sat in a magnetically shielded room and were presented with a visual stimulus consisting of a , vertical , stationary , maximum contrast , three cycles per degree , square - wave grating presented on a mean luminance background .
the stimulus was presented in the lower left visual field and subtended 8 both horizontally and vertically .
a small red fixation square was located at the top right - hand edge of the stimulus and remained on for the entire stimulation protocol ( muthukumaraswamy , 2010 ; swettenham et al , 2009 ) .
the stimulus was presented on a mitsubishi diamond pro 2070 monitor controlled by the psychophysics toolbox ( brainard , 1997 ; pelli , 1997 ) .
the screen size was 1024 by 768 pixels and the monitor frame rate was 100  hz . the monitor was outside the magnetically shielded room and viewed directly from within , at 2.15  m , through a cut - away portal in the shield .
the duration of each stimulus was 11.5  s followed by 1.5  s of fixation square only .
participants were instructed to fixate for the entire experiment and , in order to maintain attention , were instructed to 
responses slower than 750  ms from the termination of the stimulus triggered a visual warning for participants . in each recording session , 120 stimuli were presented and participants responded with their right hand .
whole - head meg recordings were made using a ctf 275-channel radial gradiometer system sampled at 1200  hz ( 0300  hz bandpass ) .
an additional 29 reference channels were recorded for noise cancellation purposes and the primary sensors were analyzed as synthetic third - order gradiometers ( vrba and robinson , 2001 ) .
three of the 275 channels were turned off because of excessive sensor noise . at the onset of each stimulus presentation
participants were fitted with three electromagnetic head coils ( nasion and pre - auriculars ) , which they wore for the entire recording day , and were localized relative to the meg system immediately before and after the recording session .
each participant had a 1  mm isotropic mri scan available for source localization analysis . to achieve mri / meg co - registration ,
the fiduciary markers were placed at fixed distances from anatomical landmarks identifiable in participants ' anatomical mris ( tragus , eye centre ) .
offline , data were first epoched from 1 to 1  s around stimulus onset and each trial visually inspected for data quality .
data with gross artefacts , such as , head movements and muscle clenching were excluded from further analysis .
reaction times were calculated for each trial relative to the ( variable ) offset of each grating stimulus .
any trial where there was no behavioural responses or a response > 750  ms were classified as incorrect responses and excluded from further analysis . for each participant
, we ensured that equal number of trials , from each recording was used in the analysis . to achieve this
, trials were removed from the end of each recording run so that each participant had an equal number of trials submitted for each condition in the analysis . following these trial removal procedures , the mean number of trials remaining to analyze was 105.5 per participant ( range 82117 ) .
three source localizations were performed on each data set using synthetic aperture magnetometry ( robinson and vrba , 1999 ) , one for gamma ( 3080  hz ) , one for alpha ( 815  hz ) , and one for evoked responses ( samerf ; robinson , 2004 ) .
correspondingly , three global covariance matrices were calculated for each data set , one for gamma , one for alpha , and one for samerf ( 0100  hz ) .
based on these covariance matrices , using the beamformer algorithm ( robinson and vrba , 1999 ) , three sets of beamformer weights were computed for the entire brain at 4  mm isotropic voxel resolution .
a multiple local - spheres ( huang et al , 1999 ) volume conductor model was derived by fitting spheres to the brain surface extracted by fsl 's brain extraction tool ( smith , 2002 ) .
readers unfamiliar with sam are referred to several methodological reviews ( hillebrand et al , 2005 ; vrba and robinson , 2001 ) . for gamma - band sam image reconstruction virtual sensors
were constructed for each beamformer voxel and student 's t - images of source power changes computed using a baseline period of 1 to 0  s and an active period of 0 to 1  s , while for alpha - band sam image reconstruction a baseline period of 1 to 0.3  s and an active period of 0.3 to 1  s was used .
this latter time window was used to avoid confounding alpha responses with the time - locked evoked response . after generating the source reconstruction images ,
the voxel with the strongest power increase ( in the contralateral occipital lobe ) was located for gamma and the voxel with the strongest source power decrease located for alpha .
the source analysis and peak voxel identification procedures were carried out separately for each recording as participants were moved in and out of the meg between recordings , which would alter their head position in the meg dewar . to reveal time - frequency responses at these locations ,
the virtual sensors at these locations were repeatedly band - passed filtered between 1 and 160  hz at 0.5  hz frequency step intervals , filtering with an 8  hz bandpass , third - order butterworth filter ( le van quyen et al , 2001 ; muthukumaraswamy et al , 2010 ) .
the hilbert transform was used to obtain the amplitude envelope and time - frequency spectra were computed as a percentage changed from the pre - stimulus baselines for each frequency band . from these spectra ,
the peak frequency and amplitude at that frequency for late gamma amplitude increases ( 3080  hz , 0.31 
s ) were quantified by collapsing across the time dimension and finding the maximal ( or minimal for alpha ) frequency .
we also checked for confounding effects of using relative change baselines by computing the average absolute values of the source amplitude estimates within the respective time - frequency bins .
for samerf , the computed evoked response was passed through the 0100  hz beamformer weights to generate samerf images ( robinson , 2004 ) at 10  ms intervals from 50150  ms .
the image window ( usually 7080  ms or 8090  ms ) with the maximal response in visual cortex was identified and the maximal voxel selected as the peak location for further analysis . for this , the evoked field was computed for this virtual sensor ( 0.2 to 0-s baseline ) and the peak amplitude and latency of the m80 and m150 responses quantified . for the production of grand - average sam maps , individual sam images were first spatially normalized onto the mni ( t1 ) average brain using fmrib 's linear affine registration tool ( jenkinson and smith , 2001 ) .
this was done by first obtaining a set of warping parameters by registering the participant 's anatomical mri with the average brain and then applying these parameters to the sam source power maps .
the grand - average images presented in figure 2 are the grand - average across both participants and conditions .
for these analyses , there were main effects of  day ' ,  session ' , and an interaction term  drug ' , presented in that order . here
, the critically important term is the interaction term , indicative of a tiagabine drug - effect .
it should be noted that for f(1 , 14 ) the critical value is 4.6 and for f(3 , 42 ) the critical value is 2.82 ( alpha=0.05 ) .
figure 1 shows summed scores from the subjective high assessment scale ( schuckit , 1980 ) that participants completed at the end of each recording session .
individual items that showed the largest increase relative to placebo were  feeling of intoxication ' ,  dizziness ' , and 
task response accuracy data showed a significant interaction effect ( f(1 , 14)=3.71 , f(3 , 42)=2.34 , f(3 , 42)=2.99 , p<0.041 ) with post hoc testing showing a weak effect at 3  h for the tiagabine session to have more errors ( tiagabine 94.9(2.1 ) % , placebo 98.5(0.6 ) % , t(14)=2.99 , p=0.06 ) .
reaction time data showed a main effect of day and an interaction effect ( f(1 , 14)=12.7 , p<0.003 , f(3 , 42)=0.73 , f(3 , 42)=3.11 , p<0.036 ) .
post hoc tests showed reaction times were significantly slowed at 3  h ( tiagabine 271(17 )  ms , placebo 215(13 )  ms , t(14)=4.719 , p<0.001 ) and 5  h ( tiagabine 251(20 )  ms , placebo 208(14 )  ms , t(14)=3.733 , p<0.002 ) after tiagabine .
figures 2a  c show grand - averaged source reconstructions for gamma - band ( 3080  hz ) , alpha - band ( 815  hz ) and evoked field ( 8090  ms ) responses to presentation of the grating stimulus and as previously reported ( brookes et al , 2005 ) these are generated in the medial visual cortex in the hemisphere opposite to the side of visual stimulation and closely colocalized .
the spatial separation of the peak locations from each other was approximately 1  cm .
time - frequency analyses for the gamma - band locations are presented in figure 2d , grand - averaged for each day and recording session , and show the typical component morphology following this type of visual stimulus .
there is an initial transient broadband ( 50100  ms ) amplitude increase in the gamma - frequency ( > 30  hz ) band (  early gamma ' ) , followed by a longer - lasting elevation of gamma - frequency amplitude in a narrower frequency range (  late gamma ' ) . in the lower frequencies
, there is an initial transient onset response , which is a characteristic signature of the evoked response ( clapp et al , 2006 ; muthukumaraswamy et al , 2010 ; swettenham et al , 2009 ) , followed by a sustained alpha amplitude decrease .
response components were parameterized ( using the optimal virtual sensor locations for each component ) for each participant and subjected to statistical analysis ( figure 3 ) .
neither early gamma frequency ( f(1 , 14)=0.093 , f(3 , 42)=2.17 , f(3 , 42)=0.73 ) or late gamma frequency ( f(1 , 14)=1.42 , f(3 , 42)=1.16 , f(3 , 42)=1.70 ) showed any significant effects . similarly , neither early gamma amplitude ( f(1 , 14)=0.23 , f(3 , 42)=0.25 , f(3 , 42)=1.57 ) or late gamma amplitude ( f(1 , 14)=0.02 , f(3 , 42)=2.24 , f(3 , 42)=1.55 ) demonstrated any effects .
alpha activity showed some non - drug specific effects including , alpha frequency showing a main effect of session but no interaction effect ( f(1 , 14)=2.32 , f(3 , 42)=3.58 , p<0.03 , f(3 , 42)=0.53 ) , while alpha amplitude showed a main effect of day effect but no interaction effect ( f(1 , 14)=7.78 , p<0.02 , f(3 , 42)=0.46 , f(3 , 42)=0.93 ) .
in addition , we performed similar statistical analyses on the individual time - frequency  pixels ' to look for subtle time - frequency effects ( data not shown ) and these demonstrated similar null effects .
finally , we checked to see if there were any confounding effects of using relative change baselines by extracting absolute amplitude values ( no baselining performed ) for the various time - frequency windows . for gamma ,
no effects were seen for the early ( f(1 , 14)=3.57 , f(3 , 42)=0.87 , f(3 , 42)=0.14 ) , late ( f(1 , 14)=3.30 , f(3 , 42)=0.15 , f(3 , 42)=0.26 ) , or baseline periods ( f(1 , 14)=1.05 , f(3  42)=0.08 , f(3 , 42)=0.21 ) . for alpha ,
no effects were seen in the active ( f(1 ,  14)=3.98 , f(3 , 42)=1.41 , f(3 , 42)=0.67 ) or the baseline period ( f(1 , 14)=0.15 , f(3 , 42)=0.08 , f(3 , 42)=0.14 ) . in figure 4
these demonstrate a clear reduction in the amplitude of the first evoked response component that occurs just after 80  ms ( labeled  m80 ' ) following tiagabine and an enhanced second component ( labeled  m150 ' ) .
statistical analysis of the m80 amplitudes data showed a main effect of day , session , and critically an interaction effect ( f(1 , 14)=24.7 , p<0.001 , f(3 , 42)=30.9 , p<0.001 , f(3 , 42)=13.1 , p<0.001 ) .
m150 amplitudes showed significant effects of day , session , and an interaction effect ( f(1 , 14)=7.26 , p<0.001 , f(3 , 42)=3.60 , p<0.03 , f(3 , 42)=4.00 , p<0.02 ) .
the m80 latency showed a day effect but no other effects ( f(1 , 14)=8.06 , p<0.02 , f(3 , 42)=0.49 , f(3 , 42)=0.14 ) and m150 latency showed no significant effects ( f(1 , 14)=0.14 , f(3 , 42)=2.34 , f(3 , 42)=0.35 ) .
no correlations were observed between the shas measures or reaction times and the evoked field magnitudes and their respective change from baseline scores .
finally , we conducted an analysis to check whether there were any systematic differences in virtual sensor locations that might potentially have confounded the results .
the virtual sensor locations for gamma , alpha , and evoked response locations were tested using anova in x , z , and z directions .
no significant interaction effects were observed , indicating no systematic changes in response locations with tiagabine .
in the current experiment , we elevated endogenous gaba in human participants by administering the gat-1 blocker tiagabine , while measuring evoked and induced visual responses with meg .
the data showed a clear modification of visual evoked responses but , contrary to our hypotheses , we found no alteration in stimulus - induced oscillatory responses in either the alpha or gamma - frequency bands .
we believe our approach did significantly increase gaba based on previous pet studies using a similar dose of tiagabine that show gaba - related changes in tracer binding ( frankle et al , 2009 , 2012 ; myers et al , 2012 ) .
one particularly striking feature of the dissociation between evoked and induced responses is that this occurred in the same temporal segments of meg data .
that is , there was a clear reduction in the m80 evoked response but no alteration in the temporally congruent early gamma spike ( 50100  ms ) .
the early gamma response may reflect early ascending input from the lateral geniculate nucleus ( castelo - branco et al , 1998 ) and is manifested across all cortical layers , although particularly in layers ii and iii ( xing et al , 2012 ) , whereas early evoked responses are mostly generated in layer iv ( kraut et al , 1985 ) .
the pharmacological basis of this selective modulation presumably lies in the complex distribution of gabaa and gabab receptor subunits across the lamina of v1 ( hendry et al , 1994 ; zilles et al , 2004 ) .
in particular , the inhibition provided by gabaa activation is thought to control response gain to stimuli in primary visual cortex ( katzner et al , 2011 ) and to control the layer - specific spread of excitatory activity caused by visual stimulation ( olivas et al , 2012 ) .
reduction of this horizontal spread of excitation ( decreased spatial summation ) during stimulation onset may be the mechanism underlying the diminished evoked responses we have observed .
this selective modulation also suggests other more - simplistic interpretations of the data are incorrect .
for example , one might argue that the reduction in m80 amplitude is due to reduced task vigilance , attention , or even fixation control as participants changed psychological state after tiagabine ingestion .
however , each of these explanations would also be expected to modulate oscillatory gamma amplitude ( hoogenboom et al , 2010 ; kahlbrock et al , 2012 ) .
similarly , as these dissociable parameters are all based on the same data , they all have the same relative signal - to - noise ratio .
the 45% peak reduction we observed in m80 amplitude is significantly larger than other studies examining visual evoked potentials ( veps ) with benzodiazepenes that affect gabaergic function . in all , 20  mg diazepam administered intravenously ( loughnan et al , 1987 ) or 10  mg administered orally ( bartel et al , 1988 ) has no effect on the vep .
however , 6 and 12  mg of the benzodiazepine bromazepam has been shown to reduce erp components beginning with the n1 ( bauer et al , 1997 ) and 10  g / kg of intravenous midazolam caused a relatively small reduction in n75-p100 amplitude ( van leeuwen et al , 1992 ) .
in fact , it has been shown that tiagabine , which increases gaba rather than acting via the benzodiazepine receptor , actually increases the amplitude of early somatosensory evoked potential components ( restuccia et al , 2002 ) .
the lack of correlation between the subjective measures in the shas and evoked field changes suggest that the changes in evoked field are not simply attributable to the arousal state of participants .
another possibility is that the decreased evoked responses we observed might be due to decreased fixation control during the tiagabine - affected conditions .
however , loss of fixation control might be expected to alter the response location of both retinotopic evoked responses ( di russo et al , 2002 ) and also gamma - band responses as these are retinotopically organized ( perry et al , 2011 ) .
our analysis of source peak activity locations indicated no shifts in spatial location for either of these response components .
further , the amplitude of gamma - band activity is very sensitive to the retinotopic position of the stimulus and falls off more rapidly with eccentricity ( swettenham et al , 2009 ) than evoked components .
nevertheless , measurement of fixation position via either eye - tracking or electrooculography would be a useful addition to future experiments .
direct manipulation of gaba levels in this experiment does not appear to have altered gamma frequency , and potential reasons for this are worth considering .
it is possible that the lack of effect in the oscillatory domain is due to a lack of statistical power , or perhaps because an inadequate dose of tiagabine was given .
however , we do not believe that it would be ethically or feasibly possible to give a much higher dose to control participants and still have them complete the required tasks given the profound behavioural and subjective effects on participants .
we also note that overall data quality was very high , with gamma activity identifiable in every subject .
further , although there is no guarantee that raising gaba levels with this dose of tiagabine has a direct effect on the synaptic ( eg , time constants ) properties of neurons in visual cortex , this is unlikely since similar doses have been shown to change 11c - flumazenil and 11c - ro15 4513 binding potentials ( frankle et al , 2009 ; myers et al , 2012 ) .
in addition , tiagabine - induced increases in gaba could potentially modulate gabaergic synapses on all neuronal subtypes , including pyramidal cells and inhibitory interneurons , and receptor subtypes ( gabaa ( six subtypes ) and gabab ) with unpredictable effects on oscillatory dynamics until the role of all these gaba receptors is elucidated . in a previous study we found , across individuals , that the frequency of stimulus - induced network gamma oscillations in primary visual cortex was positively correlated with the concentration of gaba measured with mrs .
mrs measures the total bulk of gaba in a given area of the brain ( usually from a rather large voxel , eg , 27  cm ) .
consistent with invasive evidence , mrs returns concentrations of gaba in the millimolar range ( puts and edden , 2012 ) . in the synaptic cleft
, gaba concentration can reach the millimolar range ( mody et al , 1994 ) but this rapidly declines ( < 1  ms ) .
although estimates for ambient levels gaba in the extracellular space range from tens of nanomolar to micromolar ( farrant and nusser , 2005 ) these are still orders of magnitude less than the total bulk of gaba
. therefore , mrs measures of gaba concentration probably depend on interneuron numbers / inhibitory synapse numbers , and/or the rate constants of glutamic acid decarboxylase , gaba transaminase , and gat-1 .
in particular , interneuron density / synapse numbers , which are inherently structural parameters ( schwarzkopf et al , 2012 ) , would not be susceptible to short - term pharmacological manipulations , as we have attempted here .
although our specific hypothesis was not confirmed , the data here are promising for the relatively unexplored field of pharmaco - meg ( hall et al , 2010 ) .
we note another very recent study using the cholinergic agonist phystostigmine found a selective modulation of alpha oscillation amplitude ( but not gamma ) in response to visual stimuli in humans with meg ( bauer et al , 2012 ) .
these studies therefore demonstrate that meg can noninvasively characterize the effect of pharmacologically targeted and physiologically selective agents using alterations in quantitative neuronal biomarkers .
future studies systematically investigating the pharmacological sensitivity of meg responses will allow the development of data - driven pharmacological models of these responses . given the wide range of safe and selective interventional agents that can be safely administered to interact with neurotransmitter systems , this bodes well for our future understanding of the neuropharmacological basis of meg / eeg .
 the electroencephalographic / magnetoencephalographic ( eeg / meg ) signal is generated primarily by the summation of the postsynaptic currents of cortical principal cells . at a microcircuit level , these glutamatergic principal cells are reciprocally connected to gabaergic interneurons . 
 here we investigated the relative sensitivity of visual evoked and induced responses to altered levels of endogenous gabaergic inhibition . to do this , we pharmacologically manipulated the gaba system using tiagabine , which blocks the synaptic gaba transporter 1 , and so increases endogenous gaba levels . in a single - blinded and placebo - controlled crossover study of 15 healthy participants , we administered either 15  mg of tiagabine or a placebo . 
 we recorded whole - head meg , while participants viewed a visual grating stimulus , before , 1 , 3 and 5  h post tiagabine ingestion . using beamformer source localization , we reconstructed responses from early visual cortices . 
 our results showed no change in either stimulus - induced gamma - band amplitude increases or stimulus - induced alpha amplitude decreases . 
 however , the same data showed a 45% reduction in the evoked response component at 80  ms . 
 these data demonstrate that , in early visual cortex the evoked response shows a greater sensitivity compared with induced oscillations to pharmacologically increased endogenous gaba levels . 
 we suggest that previous studies correlating gaba concentrations as measured by magnetic resonance spectroscopy to gamma oscillation frequency may reflect underlying variations such as interneuron / inhibitory synapse density rather than functional synaptic gaba concentrations . 
INTRODUCTION
PARTICIPANTS AND METHODS
Participants, Design, and Stimuli
MEG Acquisition and Analysis
RESULTS
DISCUSSION