idiopathic generalized epilepsy ( ige ) is characterized by various combinations of generalized tonic  clonic seizures , absence seizures , myoclonus ( blumenfeld , 2005 ; zhang et al . , 2011 ) , and generalized spike - and - wave discharges ( gswds ) observed during interictal periods on electroencephalography ( eeg ) recordings ( hamandi et al . , 2006 ; moeller et al . , 2011 ; zhang et al . , 2011 ) .
focal spike  wave complexes occasionally are observed interictally in patients with ige , although ictal eeg recordings show only generalized - onset seizures and no focal - onset seizures in ige ( drury and henry , 1993 ; seneviratne et al . , 2012 ) .
unlike focal or partial epilepsy , which has a confined range of influence , ige affects the whole or a larger portion of the brain often without obvious , known cause ( engel , 2001 ) . among drug - resistant epilepsy ,
patients with focal epilepsy may receive surgical resection to become seizure free . on the other hand , patients with generalized epilepsy do not have such a treatment option .
recently deep brain stimulation ( johnson et al . , 2013 ) has been hypothesized as a way to treat epilepsy patients ( fisher et al . , 2010 ) .
therefore , it is important to distinguish the driver ( or source ) versus recipient ( or sink ) to understand how the epileptic activities propagate to the entire brain . in the recent decades
, it has been generally agreed that the highly interconnected circuitry of the cortex and thalamus plays a crucial role for generalized epilepsy ( blumenfeld , 2005 ) .
there is general agreement that both cortex and thalamus participate in the generation of typical spike  wave seizures , but their relative importance is still unclear .
previous works using eeg  functional magnetic resonance imaging ( fmri ) and anatomical mr - based study ( bernhardt et al . , 2009 ) indicated the involvements of thalamus , default mode network ( fox et al .
2001 ) , cerebellum , caudate nuclei and corticocortical networks in the generations of gswds ( bernhardt et al .
however , the exact interplay between the cortical and sub - cortical structures remains to be further explored .
the goal of the present study is to use noninvasive , multimodal imaging techniques to elucidate the underlying mechanisms that generate gswds in ige patients . specifically , we aim to map the cortical networks associated with gswds and investigate the causality between cortex and thalamus during gswds .
we first performed eeg - informed fmri analysis and identified regions of interest ( roi ) .
we then tested the specific connectivity patterns by seed - based connectivity analysis in fmri data .
the seeds include both regions determined by the eeg - informed fmri analysis and additional ones identified by eeg waveforms .
rois that exhibited network properties , i.e. the ones that share temporal profile with remote regions , were further subjected to the granger causality analysis to identify sources and sinks within the networks .
eeg and fmri are two noninvasive neuroimaging tools used in epilepsy research and clinical applications .
eeg has the benefit of having high temporal resolution but often suffers from limited spatial resolution . with the advancement of source imaging techniques ( he et al . ,
2011 ; he and ding , 2013 ; holmes et al . , 2010 ; koessler et al . , 2010 ; lantz et al . ,
2003 ) , successfully localizing epileptic activity in focal epilepsy is possible ( ding et al . , 2007 ; he and ding , 2013 ; lai et al .
; lu et al . , 2012a ; michel et al . , 2004 ; sohrabpour et al . , 2015 ; yang et al . , 2011 ) . however , very few studies have looked into using dense - array eeg to study the temporal dynamics of the sources in ige ( jung et al . , 2005 ) .
with complementary high spatial and temporal resolution , simultaneous fmri and eeg ( he and liu , 2008 ; liu and he , 2008 ) has been shown to provided valuable information in diagnosis of epilepsy ( pittau et al . , 2012 ; zhang et al . , 2015 ) .
seed - based connectivity analysis using resting state fmri is another common technique in the studying of both healthy and diseased neurological networks ( greicius et al .
, 2009 ; moeller et al . , 2011 ; o'muircheartaigh et al . , 2012 ) . using this technique , remote areas that share the same temporal characteristics can be identified .
compared with eeg - informed fmri , seed - based connectivity analysis does not rely on precise knowledge of the eeg event timing , which can be difficult to obtain given the noisy signal collected concurrently with fmri .
however , it requires prior knowledge in determining the rois as seeds . in our study , we used regions identified by the eeg - informed fmri result as seed in conjunction with other areas that may potentially be involved in the network activities in generating gswds . by using seed - based analysis ,
we can obtain a more specific network level activity than using eeg - informed fmri alone .
once we establish a network that is involved in generating gswds , the next question is which node in the network is driving the others .
directed connectivity measures based on the concept of granger causality ( granger , 1969 ) has been proposed ( babiloni et al .
, 2005 ; goebel et al . , 2003 ; kaminski and blinowska , 1991 ) to discern the causal relationship among different temporal series .
the direction can be estimated with the following rationale : the driver is earlier than the recipient implying that the driver contains information about the future of the recipient not contained in the past of the recipient while the reverse is not the case . the directed transfer function ( dtf ) has been used to quantify the directionality and strength of the connectivity profile among different brain regions ( kaminski and blinowska , 1991 ) .
dtf has been successfully applied in the field of epilepsy research , to identify sources ( active or efferent sources ) and sinks ( passive or afferent sources ) that may play important roles in generating seizures and interictal activities ( ding et al .
built upon the dtf method , which computes the overall connectivity strength in a given time window , an adaptive dtf ( adtf ) method was developed to study the time - variant propagation of interictal spikes ( wilke et al .
the adtf method may be able to capture the temporal dynamics of the propagation and shed light in the inter - play among the networks in the genesis and propagation of gswds .
the overall aim of this work is to use noninvasive , multimodal imaging techniques to map the brain networks associated with gswds and investigate the functional connectivity between cortex and thalamus during gswds .
the novelty of the current work lies in using multimodal neuroimaging approaches to identify the critical roles of thalamocortical networks in generalized spike and wave discharges , including electrophysiological recordings ( eeg ) , hemodynamic recordings ( bold fmri ) , eeg - informed fmri and functional connectivity mapping .
ten patients ( mean age 33    14 , three females ) with idiopathic generalized epilepsy syndromes were recruited from the department of neurology of the university of minnesota , usa . the patients included in this study were selected based on the criteria that there were visible interictal gswds recorded from clinical eeg ( with rare focal spikes only if gswds were frequent interictally ) , normal brain mri ( or normal brain ct in subjects 5 and 10 ) , and a clinical diagnosis of ige .
the study has been approved by the institutional review board of the university of minnesota .
ten healthy volunteers ( ages 2131  years , 26    2.4  years , 6 males ) participated in this study .
all subjects had written consent according to a protocol approved by the institutional review board of the university of minnesota .
a total of two functional mri scans , each lasting for 6  min , were recorded from each subject while lying in the mr scanner quietly .
one electrode was placed on patients ' back to record cardiac activity for noise removal purpose .
electrode impedances were brought below 20  k. the eeg was amplified using mr - compatible amplifiers ( brainamp mr 64 plus , brainproducts , germany ) and recorded at 1000  hz .
two sessions of eeg were recorded both inside and outside of the scanner experiment . during outside scanner recording ,
we recorded simultaneous eeg and fmri for at least a total of 40  min for each patient .
the mr gradient artifact was removed using a principal component analysis ( pca)-based optimal basis set ( obs ) algorithm ( niazy et al . , 2005 ) . when detecting and removing the cardioballistic artifact ( cba ) , ecg signal from a single electrode on the subject 's back ,
the timing of each heartbeat artifact in this channel was determined using an r - peak detection algorithm adapted from liu et al .
the final artifact correction procedure is based on a combination of ica , obs , and an information - theoretic rejection criterion ( liu et al . , 2012 ) .
briefly , the signal is decomposed into independent components , which are rejected if the mutual information between the component 's time course and the cba artifact is sufficiently high .
the remaining components are then divided into epochs around each heartbeat and an optimal basis set is obtained across all epochs to fit and remove the artifacts .
detection of bad electrodes and data epochs was performed before cba detection , and again after cba correction .
electrodes were first re - referenced to a common average of electrodes connected to the same amplifier , and then to the combined average . together with eeg data obtained from outside of mr scanner , the eeg signal was filtered and down - sampled to 256  hz .
we used 3  t siemens magnetom trio and skyra mr scanners ( germany ) with 16 channel head coil .
the echo planar imaging ( epi ) volumes underwent several preprocessing steps including three - dimensional ( 3-d ) motion correction , slice scan time correction and linear trend removal .
all fmri data were pre - processed for slice scan time correction , 3-d motion correction and temporal filtering .
matlab based toolbox spm8 ( ashburner et al . , 2010 ) was used for eeg - informed fmri analysis .
brainvoyager qx software ( brain innovation , maastricht , netherlands ) was used for the seed - based connectivity analysis .
, there were 176 contiguous sagittal slices with 1  mm slice thickness ( matrix size : 256    256 ; fov : 256  mm    256 mm ; tr / te  =  20  ms/3.3  ms ) .
whole - brain functional images with bold contrast were acquired using gradient echo planar imaging sequence ( 32 axial 3-mm thick interleaved slices with 0.3-mm gap ; tr / te  =  2000  ms/30  ms ; flip angle  =  90 ; matrix size : 64    64 ; fov : 192  mm    192 mm ) .
structural mris were normalized via alignment to the anterior  posterior commissural line and then transformation into talairach space .
fmri data were spatially coregistered to the anatomical mri . independent component analysis ( ica ) is a widely used data - driven technique to separate spatio - temporal signals into spatial components that are independent from each other through the selected time segment .
infomax ica algorithm ( bell and sejnowski , 1995 ; delorme and makeig , 2004 ) was used to decompose the spatio - temporal electrophysiological data into multiple independent components ( ics ) using a time - by - space formulation .
ica was performed on eeg obtained both in and out of the scanner to identify gswd related components for the subsequent eeg - informed fmri analysis . in the eeg - informed fmri analysis ,
the important issue is the identification of gswd timing based on eeg collected in the scanner .
although the artifact removal algorithm used was adequate in removing the majority of noise , it is still possible that some gswds were distorted by the residual noise and were rendered difficult to identify using visual inspection . compared to baseline activity , gswds are characterized by synchronized large amplitude discharges that are present in multiple channels .
such large changes in activities are visible in both raw eeg and independent components related to the gswd activities .
therefore , the temporal correlation between the two signals can be used in this study for the selection of gswd components .
this method is similar to what was previously described in seizure imaging ( yang et al . , 2011 ) .
once an ic of interest was identified from out - scanner eeg , it can be used as a benchmark for the ic selection from in - scanner eeg .
the detailed steps are as follows . for eeg obtained outside of the scanner , timing of each gswd
each 10  s time window containing one or multiple gswds was selected and concatenated to form a gswd - dense eeg . the cross correlation between the time course of each ic and the averaged time course of all of the eeg channels
the ic with the maximum absolute cross correlation valued was selected as most representative of the gswd activity .
the correlation between the spatial map of the selected ic from in - scanner eeg and that from out - scanner eeg was computed to ensure accuracy .
since both spatial maps from in and out of the scanner represent the same gswd activities , the two should share similar spatial pattern .
the timings of gswd were then identified basing on the time course of the selected ic .
the regressor of the general linear model ( glm ) was constructed using identified time points convolving with the canonical hemodynamic response function ( hrf ) ( jann et al . , 2008 ; marques et al . , 2009 ) .
the final design matrix was composed of the regressor that represents gswd activity and the 6 movement parameters as previously described ( marques et al . , 2009 ) .
the group level analysis was performed using spm8 ( ashburner et al . , 2010 ) .
individual t - statistic images were averaged using random - filed theory to correct for multiple comparison errors .
fmri data were filtered using a band - pass filter ( 0.0090.15  hz ) to reduce low frequency drift and high - frequency noise ( lowe et al . ,
the rois used in this study were informed by the eeg - informed fmri results , with additional seeds added based on eeg waveform in bilateral superior frontal regions where there are strong gswd discharges observed on eeg .
such areas include the left and right superior gyri and the middle frontal superior gyrus .
a seed in the anterior nucleus of the thalamus was also included , as it was the stimulation target of the sante trial ( fisher et al . , 2010 ) .
seeds were selected in the brainvoyager software by referencing talairach client 's ( lancaster et al .
, 2000 , 1997 ) archive of talairach labels and selecting a central coordinate for each seed .
the time courses of both seed coordinates was regressed against all brain voxel time courses to create two brain maps of r - values for each fmri scan .
a p - value threshold less than 0.05 with correction via bonferroni multiple comparisons was used to identify which voxels were significantly correlated with the seed location .
all images were smoothed using a 2.0  mm full width at half maximum ( fwhm ) gaussian kernel within brainvoyager .
the resulting voxels were clustered and counted to record a total volume of significantly correlated connectivity for each fmri scan . at the group level , a second - level , random - effects analysis was performed .
connectivity maps were created with the same threshold levels and smoothing parameters described above . only voxels with correlation less than the p - value of 0.05 corrected using the bonferroni method ,
are reported as significant across subjects . in each of the selected well - formed gswds ,
an epoch of approximately 400  ms before and 600  ms after the peak of the spike was extracted for the subsequent continuous source localization as described above .
the distributed current density of the underlying neuronal activity was estimated to obtain the source waveform at each voxel .
time series of the source waveforms corresponding to three rois , i.e. left mediodorsal nucleus of the thalamus , right mediodorsal nucleus of the thalamus and the medial frontal cortex , were selected .
these anatomical locations were chosen based on analysis previously described in the seed - based connectivity analysis .
the three source waveforms were subjected to the dtf computation , similar to the procedure previously described ( ding et al . , 2007 ; lu et al . , 2012b ) .
nonparametric permutation tests were conducted to test the significance of the obtained directional dtf values .
the threshold was set at p  <  0.01 to consider a dtf value as significant . as a step toward computing the overall dtf value between two time series , the contribution of each frequency point , at 1  hz increment , was computed automatically . since different frequencies might carry information differently , dtf output at each frequency point was averaged across all spikes spanning from 1 to 125  hz to show the contributions of different frequency bands .
additionally , adaptive dtf ( wilke et al . , 2011 , 2008 ) was performed on individual spikes to study the time varying feature of the information flow at different time points of the spike  slow wave complex .
this measure will aid in delineating the temporal changes of the connectivity strength and determining the dynamics in initiation and propagation of the gswds .
1c shows the histogram of the temporal cross - correlation coefficients between each ic and the global field potential .
the cross correlation between the spatial weight of this ic and that of out - scanner eeg is 1 .
the cross correlation between the time course of the selected ic and the averaged time course of all the eeg is 0.73 . on a group level ,
the correlation between the spatial weight of the selected ic and that of out - scanner eeg is 0.99    0.01 .
the cross correlation between the time course of the selected ic and the mean global field potential is 0.71    0.05 .
the highest activations were observed in anterior cingulate cortex , bilateral mediodorsal nuclei , left caudate nucleus , bilateral insula and bilateral sensorimotor areas .
seed - based connectivity analysis with seed in left , right superior gyrus , and the insula showed these regions are connected bilaterally ( fig .
3a , b and e ) . but seeds in middle frontal superior gyrus , caudate nucleus and sensorimotor cortex are only temporally correlated to the close vicinities ( fig .
group average of seed - based analysis in acc for both patients with ige and healthy controls are shown in fig .
the total voxel counts among all activities that are correlated to the acc were significantly different between patients with ige and healthy controls ( p  <  0.05 by nonparametric wilcoxon test , fig .
4b ) . there was a linear positive trend between the degree of connectivity between acc and the thalami , reflected by the number of voxels in the thalami that are correlated with the acc time course , and the frequency of gswds ( fig .
4c ) . activities in the thalami that are correlated with acc are located in the bilateral mediodorsal nuclei .
group comparison of seed - based analysis in mediodorsal nuclei of the thalamus between patients with ige and healthy controls are shown in fig .
the total voxel counts among all activities that are correlated to the mediodorsal nuclei were significantly different between patients with ige and healthy controls ( p  <  0.05 , fig .
5b ) . there was a similar linear trend between the degree of connectivity and the frequency of gswds ( fig .
the blue and two red dots represent the acc area in the medial frontal cortex and the two mediodorsal nuclei respectively .
the arrows between the dots indicate the directionality of the flow , where the arrows are pointing from one area , the  source  toward another , the  sink  .
the strengths of the information flowing from the medial frontal cortex is only about half in strength as the reverse direction .
the difference in strength between the two directions is statistically significant ( p  <  0.05 ) . the time varying feature of the causality as results of the adtf analysis is shown in fig .
the blue and red traces show the group averaged temporal changes of the connectivity during gswd .
the most significant exchange seems to occur as early as 50  ms before the peak of the spike , initiated by the thalamus .
the averaged information flow spanning the entire frequency range ( 1125  hz ) is plotted in fig .
there is a considerable variability but alpha and low gamma bands appear to have the most contribution .
seed - based connectivity analysis was performed with a focus on two specific structures : anterior cingulate cortex in the medial frontal lobe and the mediodorsal nuclei in the thalamus .
these two structures have been previously mentioned in works of ige but not been specifically studied .
however , it was also thought that the sources can be much more distributed only with center of gravity located near to the midline instead of having the actual focus in the medial frontal cortex .
using the seed - based connectivity analysis , we were able to specifically target this region and find its remote connections in the thalamus .
surprisingly , activity of just the mediodorsal nuclei , not the entire thalamus , is significantly correlated to the medial frontal lobe .
activation in mediodorsal nuclei has been previously seen anecdotally in a patient with history of generalized tonic - clonic seizure ( aghakhani , 2004 ) . in agreement with this finding
clonic seizure . another previous study by moeller and colleague found that in 6 children patients with ige ( moeller et al . , 2008 ) .
they saw symmetrical medial thalamic activation , which is also in agreement with our observation .
when we put the seed in the mediodorsal nuclei , its connectivity with the medial frontal cortex was replicated as expected in the ige patients and again absent in healthy controls .
in fact , the specific cortical projections of mediodorsal nucleus to the frontal cortex has long been reported in rodents and monkeys ( krettek and price , 1977 ; leonard , 1969 ; price and drevets , 2010 ) .
however , the different degree of connection between ige patient group and the control group is unexpected . it may be explained that a strengthened connection between the two structures in each patient can make the spread of gswds and seizure activity progress quickly .
furthermore , the degrees of connectivity between these two structures also seem to be positively correlated , albeit not statistically significant , with the discharge rate of gswds .
on ige patients showed slight difference in connectivity for seed in the left mediodorsal nucleus of the thalamus compared with control .
the reason they did not see bithalamic correlation , or the connection between the thalamus and the cingulate cortex could be due to the fact that moeller and colleagues selected the gswd - free periods for their analysis , where we did not make such selection . therefore , the effect they observed is weaker compared with ours .
however , bilateral superior gyri and the middle frontal superior gyrus , where strong gswds are observed from raw eeg , do not seem to be involved in a thalamocortical level network activity .
bilateral insulae showed such correlation as well . despite the preprocessing steps taken in reducing the effect of noise , the periphery of the cortex is still prone to movement artifacts in the fmri recording
therefore , seed - based analysis in close vicinity of the periphery may contain erroneous connection caused by noise rather than neurophysiological activities , as seen in fig .
( 1995 ) reported that the anterior executive region formed by acc around the rostrum of the corpus callosum has numerous projections into motor systems , which can be linked to the motor response such as uncontrolled jerking movement in tonic 
there is also evidence indicating correlation between the neural activity in the acc and the degree of consciousness in patients with disorders of consciousness ( doc ) ( qin et al . , 2010 ) .
it was shown that slow delta wave activity was generated in the frontal area accompanying loss of consciousness post secondarily generalized partial seizures and complex partial seizures ( blumenfeld et al . , 2009 ; yang et al . , 2012
this may partially explain another significant symptom of generalized epilepsy , such as the momentary loss of awareness during or post seizures .
the central role of mediodorsal nuclei has been shown in the interconnected medial frontal cortico - striato - pallido - thalamic and amygdalo - striato - pallido - thalamic networks in multiple animal models ranging from rats to monkeys ( price and drevets , 2010 ; ray and price , 1993 , 1992 ; russchen et al . , 1987 ) .
the projections from mediodorsal nuclei to the amygdala and hippocampus form important circuitries that regulate emotions and memory ( price and drevets , 2010 ) .
not surprisingly , depression and memory loss are some of the prominent comorbidities in patients with epilepsy ( hesdorffer et al . , 2000 ;
however , we did not see explicit change in the seed - based functional connectivity between mediodorsal nuclei and amygdala or hippocampus
we initially hypothesized that the degree of global connectivity to the thalamus , which is defined as the total voxel counts within the brain that are correlated to the seed in the thalamus , may be indicative of the rate of occurrence of gswds .
however , it appeared not to be the case . instead , the number of voxels in the thalamus that are connected with the acc alone appears to be positively correlated to the rate of gswds .
since acc is probably related to the gswd signals , the degree of connection between thalamus and acc may be able to serve as a predictor of the level of activity observed on scalp .
however , because of the limited recording period ( up to 2.5  h total ) of eeg in each patient , the actual gswd rate may deviate from the value we obtained .
this limitation could influence the strength of the correlation between the connectivity and discharge rate . in this study , we used canonical hrf to convolve with spike timing to construct regressor in the eeg - informed fmri step .
canonical hrf is the most widely used function to represent the link between electrophysiological activity and the corresponding hemodynamic activity .
however , there have been several recent studies pointing out that a canonical hrf may not be the best in presenting the actual function , as it may vary from person to person or may even change from location to location ( bai et al . , 2010 ; levan et al . ,
, we opted for canonical hrf basis functions with time and dispersion derivatives which can model small differences in the latency and the duration of the peak response ( ashburner et al . , 2010 ) .
dtf and adtf of eeg activity have been previously applied in finding driving sources in epileptic activities ( ding et al .
, 2007 ; he et al . , 2011 ; lin et al . , 2009 ;
our dtf findings showed a reciprocal causal relationship between the frontal cortex and the thalamus , where thalamus serves more as a driver .
specifically , mediodorsal nuclei of the thalamus have strong projection to the medial frontal cortex in the acc area .
the reverse projection is much weaker in comparison , approximately at half of the strength .
the reciprocal directionality of granger causality is generally accepted ( ding and he , 2013 ) .
this particular thalamocortical reciprocal relationship is also in agreement with the understanding of the interconnected thalamocortical circuitry in generating spontaneous spike  waves in ige ( blumenfeld , 2005 ) .
the crucial role of thalamus as part of the cortico - thalamo - cortical network , in sustaining seizures has been shown recently by paz et al . in a rat model with optogenetics ( paz et al . ,
6c ) , which shows the temporal evolution of the connectivity strength , thalamus seems to play an important role from the initiation to the propagation of the gswds .
the main changes in the connectivity strength from the thalamus to the cortex occur as early as 50  ms before the start of the spike .
fmri , information exchange among the thalamocortical precedes the peak of the spike gswd , though the changes we observed occur only 50  ms prior to the gswd events , as opposed to the seconds level window observed by moeller et al .
first , due to volume conduction effect , the time courses extracted using weighted minimum norm type of inverse source estimation algorithms might be smeared .
other source localization approaches , such as nulling , where nulling constraints are used to cancel signals from specific cortical locations beamforming ( cheung et al .
2010 ) , might be able to decrease such effects . however , these methods are not without their own limitations .
for example , nulling beamformer requires that we know the locations and extent of the sources to be canceled . such information is often not available in practical applications .
another potential issue with our model is that we only studied the relationships between the medial frontal cortex and the mediodorsal nuclei of the thalamus .
this is because these structures were implicated by our eeg - informed fmri and fmri connectivity analysis .
other anatomical entities , such as other nuclei in the thalamus that were not included in the present study , or hippocampus , may potentially play a role in the network of generating gswds as well .
although seed - based connectivity analysis was performed on anterior nucleus and caudate nucleus of the thalamus ( supplementary fig .
clusters that are temporally correlated to the seeds were confined to the close vicinity of each seed itself .
it is possible , but unlikely , that the network activity was not in the form of temporal correlation and thus was not detected by our analysis method .
lastly , other causality measures , such as phase slope index ( psi ) ( marzetti et al . , 2008 ; nolte et al . , 2008 ) , effective connectivity methods ( dcm ) ( friston , 2009 ; murta et al . , 2012 ) and structural equation models ( sem ) may provide additional insight .
since such methods all require prior anatomically motivated assumptions , this acc - mediodorsal nuclei network may serve as a model framework .
a few other approaches using different recording systems may also help circumvent the volume conduction issue .
it was disputed whether fmri can be a viable tool to study causal relationships , especially in the context of epilepsy research ( david et al .
it has been shown recently that by using a faster sampling rate , at the order of 250500  ms repetition time ( tr ) it was possible to detect a multivariate network using granger causality in several simulation studies ( deshpande et al . , 2010a , b ) .
unfortunately , the tr used in our study ranged from 2000 to 2500  ms .
it may be too slow to delineate meaningful causal information at a scale of milliseconds or lower .
with advancement in mri acquisition techniques using multiband approaches , tr can be shortened to 400  ms ( feinberg et al .
, 2010 ; uurbil et al . , 2013 ) . at this rate , we may be able to extract causal information using fmri time courses in future studies .
invasive recording using intracranial electrodes planted in the thalamus and cortex may provide the most direct measures of the neurophysiological activity and dynamic changes of gswds .
similar approach using granger causality and electrocorticalgraphy ( ecog ) have been previously applied in seizure imaging with success ( wilke et al . , 2011 , 2010 , 2008 ) .
the source of gswd is located in the medial frontal lobe , more specifically , at the anterior cingulate cortex ( acc ) .
( 1994 ) used regional dipoles to localize source of generalized epilepsy activities , and the location of equivalent dipole was around the midline of baso - frontal area . this location suggested by eeg source localization appeared to be plausible ,
given the general predominance of eeg in fronto - central region ( aghakhani , 2004 ; montalenti et al . , 2001 ) .
however a few recently published works suggested that potential spurious results may be yielded when applying source estimation methods to wide spread spike and wave discharges ( daunizeau et al .
, 2010 ; kobayashi et al . , 2005 ; wennberg and cheyne , 2013 ) .
wennberg and cheyne ( 2013 ) reported that despite intracranial evidence of cortical origins , the scalp eeg during k - complex was localized to deep brain regions using either dipole localization or distributed current density source imaging .
while this finding was based on a low - density scalp electrode configuration ( 27 electrodes ) in patients in whom intracranial eeg was available , wennberg and cheyne 's results suggested the possibility of mislocalization of widespread bi - hemispheric activities to mid - deep brain ( these authors also examined 87-ch scalp eeg in one healthy human volunteer without ieeg recordings ) . in our study , since there were no intracranial recordings in the patients studied , we can not preclude the possibility of mislocalization of eeg sources due to technical limitations in solving the eeg inverse problem
. however , our high - density eeg recording- ( 64-channel ) based source analysis did return source locations that are in agreement with the eeg - informed fmri results , a technique ( gotman et al . , 2005 ) that is not based on inverse solutions .
the middle frontal region was implicated by both approaches in the group of patients we studied .
another approach is to integrate bold fmri to improve the eeg source localization accuracy ( dale and sereno , 1993 ; daunizeau et al . , 2010 , 2009 ;
by combining the complementary strengths of eeg and fmri , we showed consistent results concerning the originating and propagation of gswds .
eeg - informed fmri revealed multiple brain regions that may be involved in gswds . by means of seed - based fmri
, we tested the specific network level activity and found temporal correlation between cortical and bithalamic bold activities . according to the granger causality analysis the mediodorsal nuclei of the thalamus appears to serve as the main driver from the initiation and throughout the propagation of the gswds .
once validated , this work may provide insight in understanding the enigmatic etiology of generalized epilepsy and offer guidance in treatments .
b. time course of an example of independent component that has the highest correlation with the gswd activity . c. source localization using the spatial map of the independent component in b. d. source localization at different time points of the gswds .
b. time course of an example of independent component that has the highest correlation with the gswd activity .
c. source localization using the spatial map of the independent component in b. d. source localization at different time points of the gswds .
 unlike focal or partial epilepsy , which has a confined range of influence , idiopathic generalized epilepsy ( ige ) often affects the whole or a larger portion of the brain without obvious , known cause . 
 it is important to understand the underlying network which generates epileptic activity and through which epileptic activity propagates . 
 the aim of the present study was to investigate the thalamocortical relationship using non - invasive imaging modalities in a group of ige patients . 
 we specifically investigated the roles of the mediodorsal nuclei in the thalami and the medial frontal cortex in generating and spreading ige activities . 
 we hypothesized that the connectivity between these two structures is key in understanding the generation and propagation of epileptic activity in brains affected by ige . using three imaging techniques of eeg , fmri and eeg - informed fmri 
 , we identified important players in generation and propagation of generalized spike - and - wave discharges ( gswds ) . 
 eeg - informed fmri suggested multiple regions including the medial frontal area near to the anterior cingulate cortex , mediodorsal nuclei of the thalamus , caudate nucleus among others that related to the gswds . 
 the subsequent seed - based fmri analysis revealed a reciprocal cortical and bi - thalamic functional connection . through eeg - based granger causality analysis using ( dtf ) and adaptive dtf , within the reciprocal thalamocortical circuitry 
 , thalamus seems to serve as a stronger source in driving cortical activity from initiation to the propagation of a gswd . 
 such connectivity change starts before the gswds and continues till the end of the slow wave discharge . 
 thalamus , especially the mediodorsal nuclei , may serve as potential targets for deep brain stimulation to provide more effective treatment options for patients with drug - resistant generalized epilepsy . 