Skip to main content

ORIGINAL RESEARCH article

Front. Comput. Neurosci., 20 February 2017
Volume 11 - 2017 | https://doi.org/10.3389/fncom.2017.00008

Contributions of EEG-fMRI to Assessing the Epileptogenicity of Focal Cortical Dysplasia

  • 1Department of Neurology and Neurosurgery, Montreal Neurological Institute and Hospital, McGill University, Québec, QC, Canada
  • 2Neurology Department, Geneva University Hospitals, Geneva, Switzerland

Purpose: To examine the ability of the BOLD response to EEG spikes to assess the epileptogenicity of the lesion in patients with focal cortical dysplasia (FCD).

Method: Patients with focal epilepsy and FCD who underwent 3T EEG-fMRI from 2006 to 2010 were included. Diagnosis of FCD was based on neuroradiology (MRI+), or histopathology in MRI-negative cases (MRI−). Patients underwent 120 min EEG-fMRI recording session. Spikes similar to those recorded outside the scanner were marked in the filtered EEG. The lesion (in MRI+) or the removed cortex (in MRI−) was marked on the anatomical T1 sequence, blindly to the BOLD response, after reviewing the FLAIR images. For each BOLD response we assessed the concordance with the spike field and with the lesion in MRI+ or the removed cortex in MRI−. BOLD responses were considered “concordant” if the maximal t-value was inside the marking. Follow-up after resection was used as gold-standard.

Results: Twenty patients were included (13 MRI+, 7 MRI−), but in seven the EEG was not active or there were artifacts during acquisition. In all 13 studied patients, at least one BOLD response was concordant with the spike field; in 9/13 (69%) at least one BOLD response was concordant with the lesion: in 6/7 (86%) MRI+ and in 3/6 (50%) MRI− patients.

Conclusions: Our study shows a high level of concordance between FCD and BOLD response. This data could provide useful information especially for MRI negative patients. Moreover, it shows in almost all FCD patients, a metabolic involvement of remote cortical or subcortical structures, corroborating the concept of epileptic network.

Introduction

Focal cortical dysplasia (FCD), a frequent cause of drug resistant focal epilepsy, is a neuronal migration disorder characterized by abnormalities of the laminar structure of the cortex, variably associated with cytopathological features including giant (or cytomegalic) neurons, dysmorphic neurons, and balloon cells (Palmini et al., 2004). The classification of FCD was reviewed by an ILAE task force (Blumcke et al., 2011), which proposes three categories: FCD type I characterized by radial (FCD type Ia) or tangential (FCD type Ib) dyslamination of the neocortex; FCD type II refers to a cortical dyslamination and dysmorphic neurons without (type IIa) or with balloon cells (type IIb); and FCD type III occurs in combination with other types of lesion (hippocampal sclerosis in type IIIa, tumors in type IIIb, vascular malformations in type IIIc, and epileptogenic lesions acquired in early life in type IIId). Structural MRI can detect dysplastic changes in a variable percentage of cases and more easily in patients exhibiting balloon cells (type IIb) (Tassi et al., 2012); 3T MRI fails to localize the dysplastic lesion in up to 87% of patients with FCD type I and 33% of those with FCD type II (Bernasconi et al., 2011). Surgical results are variable (Palmini et al., 1991; Kloss et al., 2002; Tassi et al., 2002; Chassoux et al., 2012; Guerrini et al., 2015) and depend on the ability to delineate and fully excise the entire region of dysplastic cortex (Paolicchi et al., 2000; Wagner et al., 2011). The lesion may be small and difficult to detect on routine MRI (De Ciantis et al., 2016), parts of the dysplastic tissue may be undetectable (Bernasconi et al., 2011), or additional dysplastic cortical areas remote from the lesion exist (Fauser et al., 2009), all reasons that explain the variability of the surgical outcomes (Fauser et al., 2015).

Recording simultaneously functional MRI (fMRI) and EEG is a noninvasive method detecting cerebral hemodynamic changes related to interictal epileptic discharges (IEDs) on scalp EEG. Several studies demonstrated the ability of EEG-fMRI to characterize various forms of focal and generalized epilepsy (Gotman et al., 2005; Laufs and Duncan, 2007; Gotman, 2008) and the clinical utility of this approach (Zijlmans et al., 2007; Moeller et al., 2009; Pittau et al., 2012; An et al., 2013). Previous small series of FCD patients studied with EEG-fMRI (Federico et al., 2005; Tyvaert et al., 2008) showed that focal IEDs are associated with metabolic changes (IEDs-related BOLD activations or deactivations) in the lesion itself and the overlying cortex but also in areas and neural networks extending beyond the FCD. A larger study of patients with type II FCD (Thornton et al., 2011) comparing EEG-fMRI with intracranial EEG findings, confirmed that EEG-fMRI provides useful information on the extension and epileptogenicity of FCDs and helps to identify patients with good surgical outcome (those with discreet and solitary lesion and well localized epileptogenic area) from patients who would be less likely to benefit from resective surgery (those with remote epileptogenic areas).

The aim of our study is to further characterize EEG-fMRI changes in a population of patients with focal cortical dysplasia demonstrated by MRI or histopathology. More specifically, we wish to assess if EEG-fMRI can define the localization and extension of the FCD, and if it can contribute to determine the localization and extension of the epileptogenic area(s) explained by the primary dysplastic lesion.

Methods

Subjects

We included consecutive patients with a diagnosis of focal epilepsy and FCD who underwent EEG-fMRI from October 2006 to July 2010. The diagnosis of FCD was based on the neuroradiological findings (MRI+) and when applicable confirmed by histopathology or, in case of negative MRI (MRI−), on histopathology reports only. Some patients had surgery and resections were guided by the integration of clinical and laboratory findings excluding the EEG-fMRI results. The study was approved by the institutional research ethics board. Each subject gave written informed consent in accordance with the Research Ethics Committee of the Montreal Neurological Institute and Hospital.

EEG-fMRI Acquisition

EEG was continuously recorded inside a 3T MRI scanner (Siemens, Trio, Germany). No sedation was given. The EEG acquisition was performed with 25 MR compatible electrodes (Ag/AgCl) placed on the scalp using the 10-20 (reference at FCz) and the 10-10 (F9, T9, P9, F10, T10, and P10) placement systems. Two additional electrodes were placed to record the electrocardiogram. The head of the patient was immobilized with a pillow filled with foam microspheres (Siemens, Germany) to minimize movement artifacts and for patient's comfort. Data were transmitted from a Brain Amp amplifier (Brain Products, Munich, Germany, 5 kHz sampling rate) to the EEG monitor located outside the scanner room via an optic fiber cable.

A T1-weighted anatomical acquisition was first done (1 mm slices, 256 × 256 matrix, TE = 7.4 ms, TR = 23 ms, flip angle 30°) and used to superimpose functional images. The functional data were acquired in runs of 6 min each with the patient in the resting state using a T2*-weighted EPI sequence (64 × 64 matrix; either 25 slices, 5 × 5 × 5 mm, TE = 30 ms, TR = 1.7 s, or 33 slices, 3.7 × 3.7 × 3.7 mm, TE = 25 ms, TR = 1.9 s; flip angle 90°). The data are available upon request to the authors.

EEG-fMRI Processing

EEG

Brain Vision Analyser software (Brain Products, Munich, Germany) was used for off-line correction of the gradient artifact (Allen et al., 2000). A 50-Hz low-pass filter was also applied to remove the remaining artifact. The ballistocardiogram artifact was removed by independent component analysis (Bénar et al., 2003). A neurologist reviewed the EEG recording and marked IEDs, according to those observed during clinical monitoring (outside the scanner).

fMRI

The EPI images were motion corrected and smoothed (6 mm full width at half maximum) using the software package from the Brain Imaging Center of the Montreal Neurological Institute (http://www.bic.mni.mcgill.ca/software/). Data were then analyzed as an event-related design using fMRIstat software (Worsley et al., 2002). The EPI frames were realigned using a linear 6-parameter rigid-body transformation (3 translations and 3 rotations) to correct for movement effects. To account for residual movement artifacts, the six parameters used for the realignment were also integrated in the analysis as confound regressors in the general linear model. A regressor for each type of IED was built using the timing and duration of each event and convolved with four hemodynamic response functions (HFRs) with peaks at 3, 5, 7, and 9 s (Bagshaw et al., 2004). All these regressors were included in the same general linear model. A statistic t map was obtained for each regressor using the other regressors as confounds (a study was performed for each type of interictal event) in the fMRI analysis (fMRIstat) (Worsley et al., 2002). At each voxel, the maximum t-value was taken from the four individual t maps created with the four HRFs. To be significant, a response needed to have a spatial extent threshold of five contiguous voxels having a t > 3.1 corresponding to p < 0.05, corrected for multiple comparisons (Family Wise Error rate = FWE) resulting from the number of voxels in the brain and the use of four HRFs. The t map results were represented using red-yellow scale corresponding to positive BOLD changes (activation) and blue-white scale to negative BOLD changes (deactivation). BOLD responses related to each IED type were reviewed by three experts; responses outside the brain parenchyma were ignored.

EEG-fMRI Analysis

For each patient, the analysis proceeded as follows:

EEG and BOLD

The spike field (the region thought to generate the spike) was estimated at the sub-lobar level by visual inspection of the scalp EEG. For example, spikes seen at electrodes F7-T3-T5 corresponded to the anterior and mid aspect of left temporal lobe, or spikes generated at Fp1-F3-F7 corresponded to the anterior aspect of left frontal lobe. Each BOLD response (positive and negative) for each type of IED was considered “concordant” if the max t-value corresponded to the localization of the spike field determined by EEG. If a patient had both activation and deactivation, we selected the response with the highest t-value. We considered “not concordant” a BOLD response with a max t-value outside the spike field. We reported also the BOLD responses in subcortical areas and in areas far from the spike field, to allow the description of the correlated epileptic network.

BOLD and Lesion (MRI Positive or Resected Cortex in MRI Negative)

To allow the co-registration between the BOLD response and the lesion or the resected cortex, firstly we had to mark the lesion or the resected cortex in the anatomical T1 sequence acquired during the EEG-fMRI session. These markings were performed using MNI-Display software (http://www.bic.mni.mcgill.ca/software/Display/Display.html); they were made by a trained physician blind to the EEG-fMRI results, after revising the FLAIR images acquired separately from the EEG-fMRI session. The concordance between the lesion and the marked area were judged independently and blindly with respect to the BOLD results by two neurologists. Then, the co-registration between the marking and the EEG-fMRI results was made with Anatomist (http://brainvisa.info/download.html). Each BOLD response (positive and negative) for each type of IED was considered “concordant” if the maximal t-value overlapped with the marking. For MRI—we considered concordant also a distance up to 1 cm away from the border of the marking, because of the possible post-surgical brain displacement.

To estimate the extent of the BOLD response compared to the anatomical lesion, we computed, for the MRI positive patients, the volume of the lesion and the volume of the BOLD cluster (with t > 3.1) concordant with the lesion.

Results

Subjects

Twenty patients (10 males) with FCD and refractory focal epilepsy underwent EEG-fMRI studies. Seven were excluded from the analysis because of the lack of IEDs (n = 5) or artifacts (n = 2) during scanning. Therefore, we studied 13 subjects, all with definite FCD. Mean age at study was 30 years (range, 11–50). Mean age at seizure onset was 8.2 years (range, 1–17). Clinical details of these patients are described in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Clinical characteristics of the 13 patients showing IEDs during EEG-fMRI.

In seven patients FCD was detected by MRI: in six the MRI was characteristic for a transmantle FCD and in one (#4) a thickness and blurring of the cortex was associated with a suspicious sulcus appearance (Bernasconi et al., 2011). Two MRI+ patients were operated and the histopathological analysis showed the presence of FCD type IIb. In six MRI− patients, FCD was detected by histopathological analysis of the surgical sample. This analysis showed the presence of FCD type IIb in 3 cases, IIa in 2 cases, and Ia in 1 case.

Six patients underwent invasive EEG investigations, one MRI+ and five MRI− (Table 1). The dysplastic lesion was frontal or fronto-central in 11 cases and parietal in two.

Eight patients were operated and type 1 was found in one (MRI−), type IIa in two (both MRI −) and type IIb in five (2 MRI+ and 3 MRI−).

Concerning post-surgical follow-up, two patients (1 MRI+ and 1 MRI−) were in Engel 1; 5 were in Engel 3 or 4 (1 MRI+, 4 MRI−). Follow up was not available for patient 11.

EEG-fMRI Data

Of the 13 patients with active EEG during acquisition, 10 patients had a single type of IED and three patients had two types of IED, so 16 analyses were computed. In 14 analyses IEDs were unilateral, in one they were bilateral and in one they were bilateral and diffuse. The number of events recorded during the fMRI session ranged from 4 to 459 (average 90, median 29).

A BOLD response was observed for each analysis, and it was characterized by both activation and deactivation. In all 13 patients at least one study was concordant with the spike field (Table 2). In 9/13 (69%) patients at least one BOLD response was concordant with the lesion: in 6/7 (86%) MRI+ patients (Figures 15) and in 3/6 (50%) MRI− patients (Figure 6). In eight patients the response concordant with the lesion was an activation (Figures 1, 2, 4, 5) and in one a deactivation (Figure 3, Table 2).

TABLE 2
www.frontiersin.org

Table 2. EEG-fMRI results in the 13 patients showing IEDs during EEG-fMRI; concordance between BOLD, spike field and lesion.

FIGURE 1
www.frontiersin.org

Figure 1. Thirty two year old patient (case 1) with FCD and macrogyria involving the lips of the middle and inferior left frontal gyri (indicated with a green arrow in the FLAIR sequences acquired outside the EEG-fMRI, and marked in white in the anatomical sequence acquired during EEG-fMRI). The marked events were spikes and slow waves maximum at C3 (average reference). The BOLD response showed a focal activation with a maximum t-value in the lesion, a cluster in the homologous contralateral region, and deactivation clusters in the precuneus, globi pallidi, and left parahippocampal gyrus.

FIGURE 2
www.frontiersin.org

Figure 2. Thirty four year old patient (case 2) with FCD (marked in white) in the paramedial portion of the left precuneus (the green arrow indicates the lesion in the coronal FLAIR sequence). The marked events were polyspikes max at CzPzP3 (average reference). The BOLD response showed a focal activation with a maximum t-value in the lesion.

FIGURE 3
www.frontiersin.org

Figure 3. Twenty six year old patient (case 3) with FCD (marked in white in the anatomical sequence acquired during EEG-fMRI and indicated with a green arrow in the coronal FLAIR sequence) in the L second frontal gyrus. The marked events were low amplitude spikes max at F3 (referential montage). The BOLD response showed a focal deactivation with a maximum t-value in the lesion.

FIGURE 4
www.frontiersin.org

Figure 4. Twenty three year old patient (case 4). (A) Structural MRI shows a thickening of the cortical gray matter and blurring of the gray-white matter interface in the right frontal operculum (in white). (B) The marked events were spikes with maximum amplitude at Fp2 F8 F10 (reference Fz). (C) BOLD response was characterized by different clusters of activation. That with max t value is in the right frontal operculum, corresponding to the lesion. (D) The marked events were burst of spike and wave complexes bilateral and diffuse (bipolar montage), but more evident on the right fronto temporal region. (E) BOLD response was characterized by different clusters of activation: those with the highest t-values were in the right frontal operculum, in the right mid cingulate cortex and cerebellum; clusters of deactivation were present in the areas of the DMN.

FIGURE 5
www.frontiersin.org

Figure 5. Two patients (cases 6, 7) whose structural MRI showed a FCD lesion (marked in white on the anatomical acquired during EEG-fMRI after comparison with the FLAIR sequences). The BOLD response showed an activation with max t-value concordant with the marked lesion.

FIGURE 6
www.frontiersin.org

Figure 6. Three patients (cases 9, 10, and 11) whose structural MRI was negative. The removed cortex (marked in white on the anatomical acquired during EEG-fMRI after comparison with the post-surgical MRI) showed the presence of FCD. The BOLD response showed an activation with a max t-value concordant with the removed cortex.

In 10/16 (62%) fMRI studies, concordance was found between highest BOLD response and lesion (in 9 the highest t-value was inside the lesion or removed cortex and in 1 it was within 1 cm from the border of the resected cortex). In six studies the BOLD response with max t-value was not concordant with the lesion/resection: patient 5 (MRI+, 2 studies) with the MRI lesion located in the right fronto-polar region and the two analyses showing the highest activation in the pons and gyrus rectus bilaterally (Figure 7); patient 8 (MRI−, one study) had the BOLD max t-value in the left post-central cortex, but the resection was just anterior to the motor cortex (Figure 8); patient 12 (MRI−, 2 studies) had one study with the max t-value localized in the right fusiform gyrus (Figure 9) and the other with the highest t-value in the right superior frontal gyrus, and in this patient the right supplementary motor area was resected; and, finally, in patient 13 (MRI−, one study), the max t-value was in the left anterior frontal lobe but more anterior than the area resected (Figure 9).

FIGURE 7
www.frontiersin.org

Figure 7. Thirty year old patient (case 5) with right fronto polar FCD (marked in white). The marked event was a prolonged burst of rhythmic delta wave with equipotential at F10-T10 (bipolar montage). The BOLD response was characterized by a cluster with max t-value on both giry recti and an activation with a lower t-value in the area subsequently removed. We considered this response as not concordant.

FIGURE 8
www.frontiersin.org

Figure 8. Twenty year old patient (case 8) with seizures characterized by a tingling sensation starting in the right hand digits and then propagating to the forearm, followed by right arm stiffness and jerks and sometimes secondarily generalization. The structural MRI did not show any lesion. The marked events were spikes with phase reversal at C3 (bipolar montage). The BOLD showed different clusters of activation: that with the highest t-value is in the postrolandic mid-sensory cortex. The intracerebral chronic electrocorticogram confirmed this finding. During surgery the pre-central gyrus was removed (white marking) and the histopatology revealed the presence of dysplastic cells. This case was considered as not concordant.

FIGURE 9
www.frontiersin.org

Figure 9. Above: case 12. The marked events were spikes max at P10 O2 (the electrode T6 had an artifact during the acquisition). The BOLD response showed an activation with max t value in the right fusiform gyrus. The removed cortex cointaining the FCD was the right supplementary motor area (marked in white). Below: case 13. The marked events were T3 spikes. This patient had three operations and the surgical bed is evident on the structural MRI. The cortex removed during the last operation (marked in white) showed the presence of a FCD. The BOLD response showed a max deactivation more anterior than the removed cortex. These 2 cases are not concordant.

In five of the seven MRI+ patients, the volume of the BOLD cluster concordant with the lesion was larger than the lesion itself (Table 2). In six of seven patients this cluster had the highest t-value of the analysis (six concordant MRI+ patients).

In 15 studies, concomitant cortical activation clusters with a lower t-value were seen remotely from the lesion: in two the activation was in the contralateral homologous area, and in the others they were either in the same or in different lobes. Nine studies had activation clusters in subcortical structures: five in the cerebellum, one in the pons, one in thalami, one in putamen, and both caudate nuclei, one in hypothalamus, and cerebellum.

Deactivation cortical clusters were also observed in all 16 studies, remotely from the lesion: in six cases deactivation responses were present in brain regions belonging to the default mode network, in the others they were in different lobes. Six studies had deactivation clusters in subcortical structures: two in caudate nucleus, one in mesencephalon, and hippocampus, one in hippocampus, one in globus pallidus, and one in nucleus accumbens.

Discussion

This study describes the EEG-fMRI findings of patients with focal epilepsy and FCD. The first finding is that all patients (13/13) showing an active EEG during the acquisition had a BOLD response (activation and deactivation). This high ratio is concordant with a previous 3T study (Thornton et al., 2011) where 92% of the patients (11/12) also showed significant spike-related hemodynamic changes. Although we cannot exclude that these results are explained by a higher sensitivity of 3T EEG-fMRI compared to 1.5T, as expected from higher magnet strength (Kruger et al., 2001; Gholipour et al., 2011; Grouiller et al., 2016), this high rate of responses could be due to the high degree of epileptogenicity of FCD (Palmini et al., 1991; Gambardella et al., 1996). This intrinsic high epileptogenicity has been related to the malfunction of GABA-mediated inhibition in the dysplastic neocortex (Spreafico et al., 1998; Avoli et al., 1999; Sakakibara et al., 2012; Medici et al., 2016) and in peri-lesional tissue (Hodozuka et al., 2006), as demonstrated by electrophysiological and histochemical studies. MR spectroscopy (MRS)/BOLD studies have also shown an inverted relationship between the resting state GABA concentration, and amplitude of BOLD responses (Muthukumaraswamy et al., 2012). Hence, the high number of BOLD responses found in patients with FCD may be related to the decreases or malfunction of GABA-mediated inhibition described in the brain of the patients with this type of lesion.

This study showed overall a good concordance between localization of the highest BOLD response and localization of the lesion: this was found in 69% of the patients that had an active EEG during acquisition. In patients with a lesion detectable by structural MRI, concordance was excellent and found in 86%, but EEG-fMRI was also able to localize the lesion in 50% of patients with negative structural MRI. In such MRI− cases, the location of BOLD responses could be used to help the detection of occult lesions by targeted revision of MRI studies or to guide intracranial electrode implantation in patients with focal epilepsy likely explained by an invisible FCD, after integration with information from the clinical context and other diagnostic techniques (Guerrini et al., 2015).

In our series, just 8 of 13 included patients underwent surgical treatment; indeed often the epileptogenic zone was located in eloquent areas. Although with very small numbers, we can observe that the operated cases (one MRI+ and one MRI−) who were seizure free at >12 months of follow up, had a BOLD response concordant with the lesion. Concerning patients who kept having seizures (one MRI+ and 4 MRI−), the BOLD was not concordant to the resected area in 4/5 cases, although it was concordant with the spike field in all but one. This finding, although on small number, strengthens the concept of negative predictive value, as expressed by An et al. (2013) and Coan et al. (2016). The design of our study allowed us indeed to detect the presence of FCD in the areas that were surgically removed, but does not allow us to assess if the areas where we observed a BOLD responses contain FCD. This can be assumed at least in some cases, particularly because all the BOLD responses were concordant with the spike field. For example in patient 8, seizure semiology and EEG data obtained from a cortical grid were highly suggestive of a post-central focus, and the EEG-fMRI study from spikes with phase reversal at C3-P3 showed the highest BOLD response in that region. However, during surgery the pre-central gyrus appeared macroscopically abnormal and was removed, and the histopathology revealed FCD type IIb. The patient was not seizure-free after the resection, suggesting a more widespread lesion and a postcentral gyrus possibly also containing dysplastic tissue. Another example is patient 12, in whom intracerebral EEG revealed multifocal interictal epileptic activity arising from the right supplementary motor area (SMA) and right posterior hippocampus. The semiology and ictal EEG activity supported an early SMA involvement, and this region was resected showing a type Ia FCD. The BOLD response with the highest t-value was however obtained from the posterior quadrant spikes, with activation in the right fusiform gyrus. This patient is not seizure-free, and the EEG-fMRI findings may have again suggested a diffuse or multifocal cortical dysplasia. Previous studies showed that FCD may be associated with multiple areas of epileptogenicity (Chassoux et al., 2000; Fauser et al., 2009), some appearing structurally normal. An EEG-fMRI study correlating spike-related hemodynamic changes, intracranial EEG findings and surgical outcome (Thornton et al., 2011) supports this evidence.

Our results also suggest that the BOLD response, and hence the metabolic demand, is not dependent on the histological FCD subtype. Eight patients were operated and type I was found in one (MRI−), type IIa in two (also MRI−) and type IIb in five (2 MRI+ and 3 MRI−).

The BOLD responses with the highest t value were activations in 14 analyses, and deactivations in two. In one patient (patient 3), the negative BOLD response was concordant with the spike field and with the lesion, and in another (patient 13), although the BOLD was concordant with the spike field, it was not concordant with the lesion. This is in agreement with previous EEG-fMRI studies in FCD patients that found deactivations in the epileptic focus in only a small proportion of cases (Federico et al., 2005; Thornton et al., 2011). The localizing value of the negative BOLD is still debated and its significance is not completely clarified (Kobayashi et al., 2006a). Indeed localizing deactivation has been associated in a few cases with an earlier activation (Rathakrishnan et al., 2010); in other cases, it can occur in the region of IED generation, mostly if such IEDs are accompanied by high amplitude slow waves (Pittau et al., 2013). In our study, we distinguished the deactivations belonging to the DMN from those with a possible localizing value in order to divide responses coming from a network inhibition from those useful to the localization of the focus. A deactivation in the areas of DMN was present in 38% of the analyses; as demonstrated by other studies, this proportion is less consistent than the DMN deactivation present in other type of common pathology, as mesial temporal lobe epilepsy (Laufs and Duncan, 2007; Coan et al., 2016), most probably because HRF peak amplitude in deactivation clusters is larger in the mesial temporal lobe epilepsy group than in the FCD, when the deactivation occurred in DMN regions (Watanabe et al., 2014).

In this study, we found that BOLD responses are generally much larger than the visible lesion. We observed that in only one study (patient 2) the hemodynamic response was very focal and coinciding with the FCD; even in this case the volume of the BOLD cluster concordant with the lesion was bigger than the visible anatomical anomaly. This discrepancy between the volumes of the anatomical lesion and the BOLD cluster containing the lesion is present in almost all cases, and it could depend on the rate of IEDs or on the statistical models used for the analysis. However, the same phenomenon is common in another technique used to localize interictal focus, i.e., Positron Emission Tomography (PET). The focal interictal hypometabolism of radio-labeled fluoro-deoxy-glucose-PET is usually larger than the epileptogenic cortex, reflecting probably the altered function not only of the interictal focus, but also of the areas involved by the first ictal spread (Nelissen et al., 2006; Pittau et al., 2014a). In all the other studies, remote cortical or subcortical structures were activated suggesting that focal IEDs recorded from scalp EEG may represent only a fraction of broader metabolic events that involve widespread epileptic network (for review: Pittau et al., 2014b). These findings are in agreement with previous studies on focal epilepsy from malformations of cortical development (Kobayashi et al., 2006b) or from heterogeneous etiology (Fahoum et al., 2012). These findings are strengthened by other localizing technology, as stereo-EEG (SEEG): Varotto et al. (2012) in a series of patients with type II FCD studied with SEEG demonstrated a specific connectivity pattern (mainly in the gamma band) between the epileptogenic zone and other distant cortical regions during the ictal, inter- and pre-ictal periods. They suggested that the lesional nodes act as the center of the epileptic network where seizures originate and are sustained, and that the cortical regions beyond the dysplasia, also involved in the ictal activity, essentially act as “secondary” generators of synchronous activity.

A limitation of our study is that concordance between EEG and BOLD has been established by visual analysis. Electric source analysis (ESI) is an alternative and precise method to localize interictal spikes, when obtained from an electrode setup of at least 64 electrodes (Lantz et al., 2003). Indeed it has been shown (Brodbeck et al., 2011) that ESI has a high sensitivity (84.1%) and specificity (87.5%) on localization precision when obtained with high-resolution electric source imaging (>128-electrodes) and individual MRI as the head model. However, in the same study, sensitivity and specificity decreased to 66 and 54%, respectively when EEG was recorded with a low number of electrodes (19–29 channels). Studies combining high density EEG and fMRI will solve this limitation (Vulliemoz et al., 2009).

It is known that electrodes placed inside a FCD usually show almost continuous spiking activity (Palmini et al., 1995), but on scalp EEG much less frequent spikes are detected. How is it possible that these scalp spikes also demonstrate a source that appears to be in the FCD? One can speculate that occasionally the dysplastic focus generates a widespread discharge that reaches the scalp, and that at the time of this discharge the metabolic activity is even more intense in the dysplasia than during “background” more localized spiking, thus generating a maximal BOLD response in the dysplasia.

To conclude, our study shows a high level of concordance between FCD and BOLD response. In some MRI-negative cases the BOLD response did not correspond to the resection, but we provided evidence that the region of BOLD response could nevertheless correspond to dysplastic and epileptogenic tissue. This could be useful especially for MRI negative patients, where the location of BOLD responses can be used to detect occult lesions by targeted revision of MRI studies, or for guiding intracranial electrode implantation. However, in almost all FCD patients, a metabolic involvement of remote cortical or subcortical structures was also found, corroborating the idea that focal IEDs recorded from scalp EEG may represent only a fraction of broad metabolic events that involve a widespread epileptic network in spite of their focal appearance.

Ethics Statement

This study was carried out in accordance with the recommendations of Neurosciences Panel of the MUHC Research Ethics Board with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the Neurosciences Panel of the MUHC Research Ethics Board.

Author Contributions

FP: Substantial contributions to the conception or design of the work. Recording and analysis EEG-fMRI data. Writing the manuscript. Revision and interpretation data. LF: Substantial contributions to the conception or design of the work. Analysis EEG-fMRI and making figure. FF: Recording and analysis EEG-fMRI data. Revising it critically for important intellectual content; FD: recruiting patient. Revision and interpretation data. Revising it critically for important intellectual content; JG: idea of the project. Revision and interpretation data. All authors: Agreement to be accountable for all of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. Final approval of the version to be published.

Funding

This work was supported by the Canadian Institutes of Health Research (CIHR) grant MOP-38079. FP was supported by the Savoy Foundation for epilepsy.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors thank Natalja Zazubovits for helping to collect and analyze the data and Dr. D. K. Nguyen for contributing with patient data.

References

Allen, P. J., Josephs, O., and Turner, R. (2000). A method for removing imaging artifact from continuous EEG recorded during functional MRI. Neuroimage 12, 230–239. doi: 10.1006/nimg.2000.0599

PubMed Abstract | CrossRef Full Text | Google Scholar

An, D., Fahoum, F., Hall, J., Olivier, A., Gotman, J., and Dubeau, F. (2013). Electroencephalography/functional magnetic resonance imaging responses help predict surgical outcome in focal epilepsy. Epilepsia 54, 2184–2194. doi: 10.1111/epi.12434

PubMed Abstract | CrossRef Full Text | Google Scholar

Avoli, M., Bernasconi, A., Mattia, D., Olivier, A., and Hwa, G. G. (1999). Epileptiform discharges in the human dysplastic neocortex: in vitro physiology and pharmacology. Ann. Neurol. 46, 816–826. doi: 10.1002/1531-8249(199912)46:6<816::AID-ANA3>3.0.CO;2-O

PubMed Abstract | CrossRef Full Text | Google Scholar

Bagshaw, A. P., Aghakhani, Y., Benar, C. G., Kobayashi, E., Hawco, C., Dubeau, F., et al. (2004). EEG-fMRI of focal epileptic spikes: analysis with multiple haemodynamic functions and comparison with gadolinium-enhanced MR angiograms. Hum. Brain Mapp. 22, 179–192. doi: 10.1002/hbm.20024

PubMed Abstract | CrossRef Full Text | Google Scholar

Bénar, C., Aghakhani, Y., Wang, Y., Izenberg, A., Al-Asmi, A., Dubeau, F., et al. (2003). Quality of EEG in simultaneous EEG-fMRI for epilepsy. Clin. Neurophysiol. 114, 569–580. doi: 10.1016/S1388-2457(02)00383-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernasconi, A., Bernasconi, N., Bernhardt, B. C., and Schrader, D. (2011). Advances in MRI for 'cryptogenic' epilepsies. Nat. Rev. Neurol. 7, 99–108. doi: 10.1038/nrneurol.2010.199

PubMed Abstract | CrossRef Full Text | Google Scholar

Blumcke, I., Thom, M., Aronica, E., Armstrong, D. D., Vinters, H. V., Palmini, A., et al. (2011). The clinicopathologic spectrum of focal cortical dysplasias: a consensus classification proposed by an ad hoc task force of the ILAE diagnostic methods commission. Epilepsia 52, 158–174. doi: 10.1111/j.1528-1167.2010.02777.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Brodbeck, V., Spinelli, L., Lascano, A. M., Wissmeier, M., Vargas, M. I., Vulliemoz, S., et al. (2011). Electroencephalographic source imaging: a prospective study of 152 operated epileptic patients. Brain 134(Pt 10), 2887–2897. doi: 10.1093/brain/awr243

PubMed Abstract | CrossRef Full Text | Google Scholar

Chassoux, F., Devaux, B., Landre, E., Turak, B., Nataf, F., Varlet, P., et al. (2000). Stereoelectroencephalography in focal cortical dysplasia: a 3D approach to delineating the dysplastic cortex. Brain 123(Pt 8), 1733–1751. doi: 10.1093/brain/123.8.1733

PubMed Abstract | CrossRef Full Text | Google Scholar

Chassoux, F., Landre, E., Mellerio, C., Turak, B., Mann, M. W., Daumas-Duport, C., et al. (2012). Type II focal cortical dysplasia: electroclinical phenotype and surgical outcome related to imaging. Epilepsia 53, 349–358. doi: 10.1111/j.1528-1167.2011.03363.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Coan, A. C., Chaudhary, U. J., Grouiller, F., Campos, B. M., Perani, S., De Ciantis, A., et al. (2016). EEG-fMRI in the presurgical evaluation of temporal lobe epilepsy. J. Neurol. Neurosurg. Psychiatry 87, 642–649. doi: 10.1136/jnnp-2015-310401

PubMed Abstract | CrossRef Full Text | Google Scholar

De Ciantis, A., Barba, C., Tassi, L., Cosottini, M., Tosetti, M., Costagli, M., et al. (2016). 7T MRI in focal epilepsy with unrevealing conventional field strength imaging. Epilepsia 57, 445–454. doi: 10.1111/epi.13313

PubMed Abstract | CrossRef Full Text | Google Scholar

Fahoum, F., Lopes, R., Pittau, F., Dubeau, F., and Gotman, J. (2012). Widespread epileptic networks in focal epilepsies: EEG-fMRI study. Epilepsia 53, 1618–1627. doi: 10.1111/j.1528-1167.2012.03533.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Fauser, S., Essang, C., Altenmüller, D. M., Staack, A. M., Steinhoff, B. J., Strobl, K., et al. (2015). Long-term seizure outcome in 211 patients with focal cortical dysplasia. Epilepsia 56, 66–76. doi: 10.1111/epi.12876

CrossRef Full Text | Google Scholar

Fauser, S., Sisodiya, S. M., Martinian, L., Thom, M., Gumbinger, C., Huppertz, H. J., et al. (2009). Multi-focal occurrence of cortical dysplasia in epilepsy patients. Brain, 132(Pt 8), 2079–2090. doi: 10.1093/brain/awp145

PubMed Abstract | CrossRef Full Text | Google Scholar

Federico, P., Archer, J. S., Abbott, D. F., and Jackson, G. D. (2005). Cortical/subcortical BOLD changes associated with epileptic discharges: an EEG-fMRI study at 3 T. Neurology 64, 1125–1130. doi: 10.1212/01.WNL.0000156358.72670.AD

PubMed Abstract | CrossRef Full Text | Google Scholar

Gambardella, A., Palmini, A., Andermann, F., Dubeau, F., Da Costa, J. C., Quesney, L. F., et al. (1996). Usefulness of focal rhythmic discharges on scalp EEG of patients with focal cortical dysplasia and intractable epilepsy. Electroencephalogr. Clin. Neurophysiol. 98, 243–249. doi: 10.1016/0013-4694(95)00266-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Gholipour, T., Moeller, F., Pittau, F., Dubeau, F., and Gotman, J. (2011). Reproducibility of interictal EEG-fMRI results in patients with epilepsy. Epilepsia 52, 433–442. doi: 10.1111/j.1528-1167.2010.02768.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Gotman, J. (2008). Epileptic networks studied with EEG-fMRI. Epilepsia 49(Suppl. 3), 42–51. doi: 10.1111/j.1528-1167.2008.01509.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Gotman, J., Grova, C., Bagshaw, A., Kobayashi, E., Aghakhani, Y., and Dubeau, F. (2005). Generalized epileptic discharges show thalamocortical activation and suspension of the default state of the brain. Proc. Natl. Acad. Sci. U.S.A. 102, 15236–15240. doi: 10.1073/pnas.0504935102

PubMed Abstract | CrossRef Full Text | Google Scholar

Grouiller, F., Jorge, J., Pittau, F., van der Zwaag, W., Iannotti, G. R., Michel, C. M., et al. (2016). Presurgical brain mapping in epilepsy using simultaneous EEG and functional MRI at ultra-high field: feasibility and first results. MAGMA 29, 605–616. doi: 10.1007/s10334-016-0536-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Guerrini, R., Duchowny, M., Jayakar, P., Krsek, P., Kahane, P., Tassi, L., et al. (2015). Diagnostic methods and treatment options for focal cortical dysplasia. Epilepsia 56, 1669–1686. doi: 10.1111/epi.13200

PubMed Abstract | CrossRef Full Text | Google Scholar

Hodozuka, A., Tsuda, H., Hashizume, K., and Tanaka, T. (2006). Focal cortical dysplasia: pathophysiological approach. Childs Nerv. Syst. 22, 827–833. doi: 10.1007/s00381-006-0136-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kloss, S., Pieper, T., Pannek, H., Holthausen, H., and Tuxhorn, I. (2002). Epilepsy surgery in children with focal cortical dysplasia (FCD): results of long-term seizure outcome. Neuropediatrics 33, 21–26. doi: 10.1055/s-2002-23595

PubMed Abstract | CrossRef Full Text | Google Scholar

Kobayashi, E., Bagshaw, A. P., Grova, C., Dubeau, F., and Gotman, J. (2006a). Negative BOLD responses to epileptic spikes. Hum. Brain Mapp. 27, 488–497. doi: 10.1002/hbm.20193

PubMed Abstract | CrossRef Full Text | Google Scholar

Kobayashi, E., Bagshaw, A. P., Grova, C., Gotman, J., and Dubeau, F. (2006b). Grey matter heterotopia: what EEG-fMRI can tell us about epileptogenicity of neuronal migration disorders. Brain 129(Pt 2), 366–374. doi: 10.1093/brain/awh710

PubMed Abstract | CrossRef Full Text | Google Scholar

Kruger, G., Kastrup, A., and Glover, G. H. (2001). Neuroimaging at 1.5 T and 3.0 T: comparison of oxygenation-sensitive magnetic resonance imaging. Magn. Reson. Med. 45, 595–604. doi: 10.1002/mrm.1081

PubMed Abstract | CrossRef Full Text | Google Scholar

Lantz, G., Grave de Peralta, R., Spinelli, L., Seeck, M., and Michel, C. M. (2003). Epileptic source localization with high density EEG: how many electrodes are needed? Clin. Neurophysiol. 114, 63–69. doi: 10.1016/S1388-2457(02)00337-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Laufs, H., and Duncan, J. S. (2007). Electroencephalography/functional MRI in human epilepsy: what it currently can and cannot do. Curr. Opin. Neurol. 20, 417–423. doi: 10.1097/WCO.0b013e3282202b92

PubMed Abstract | CrossRef Full Text | Google Scholar

Medici, V., Rossini, L., Deleo, F., Tringali, G., Tassi, L., Cardinale, F., et al. (2016). Different parvalbumin and GABA expression in human epileptogenic focal cortical dysplasia. Epilepsia 57, 1109–1119. doi: 10.1111/epi.13405

PubMed Abstract | CrossRef Full Text | Google Scholar

Moeller, F., Tyvaert, L., Nguyen, D. K., LeVan, P., Bouthillier, A., Kobayashi, E., et al. (2009). EEG-fMRI: adding to standard evaluations of patients with nonlesional frontal lobe epilepsy. Neurology 73, 2023–2030. doi: 10.1212/WNL.0b013e3181c55d17

PubMed Abstract | CrossRef Full Text | Google Scholar

Muthukumaraswamy, S. D., Evans, C. J., Edden, R. A., Wise, R. G., and Singh, K. D. (2012). Individual variability in the shape and amplitude of the BOLD-HRF correlates with endogenous GABAergic inhibition. Hum. Brain Mapp. 33, 455–465. doi: 10.1002/hbm.21223

PubMed Abstract | CrossRef Full Text | Google Scholar

Nelissen, N., Van Paesschen, W., Baete, K., Van Laere, K., Palmini, A., Van Billoen, H., et al. (2006). Correlations of interictal FDG-PET metabolism and ictal SPECT perfusion changes in human temporal lobe epilepsy with hippocampal sclerosis. NeuroImage 32, 684–695. doi: 10.1016/j.neuroimage.2006.04.185

PubMed Abstract | CrossRef Full Text | Google Scholar

Palmini, A., Andermann, F., Aicardi, J., Dulac, O., Chaves, F., Ponsot, G., et al. (1991). Diffuse cortical dysplasia, or the 'double cortex' syndrome: the clinical and epileptic spectrum in 10 patients. Neurology 41, 1656–1662. doi: 10.1212/WNL.41.10.1656

PubMed Abstract | CrossRef Full Text | Google Scholar

Palmini, A., Gambardella, A., Andermann, F., Dubeau, F., da Costa, J. C., Olivier, A., et al. (1995). Intrinsic epileptogenicity of human dysplastic cortex as suggested by corticography and surgical results. Ann. Neurol. 37, 476–487. doi: 10.1002/ana.410370410

PubMed Abstract | CrossRef Full Text | Google Scholar

Palmini, A., Najm, I., Avanzini, G., Babb, T., Guerrini, R., Foldvary-Schaefer, N., et al. (2004). Terminology and classification of the cortical dysplasias. Neurology 62(6 Suppl. 3), S2–S8. doi: 10.1212/01.wnl.0000114507.30388.7e

PubMed Abstract | CrossRef Full Text | Google Scholar

Paolicchi, J. M., Jayakar, P., Dean, P., Yaylali, I., Morrison, G., Prats, A., et al. (2000). Predictors of outcome in pediatric epilepsy surgery. Neurology 54, 642–647. doi: 10.1212/WNL.54.3.642

PubMed Abstract | CrossRef Full Text | Google Scholar

Pittau, F., Dubeau, F., and Gotman, J. (2012). Contribution of EEG/fMRI to the definition of the epileptic focus. Neurology 78, 1479–1487. doi: 10.1212/WNL.0b013e3182553bf7

PubMed Abstract | CrossRef Full Text | Google Scholar

Pittau, F., Fahoum, F., Zelmann, R., Dubeau, F., and Gotman, J. (2013). Negative BOLD response to interictal epileptic discharges in focal epilepsy. Brain Topogr. 26, 627–640. doi: 10.1007/s10548-013-0302-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Pittau, F., Grouiller, F., Spinelli, L., Seeck, M., Michel, C. M., and Vulliemoz, S. (2014a). The role of functional neuroimaging in pre-surgical epilepsy evaluation. Front. Neurol. 5:31. doi: 10.3389/fneur.2014.00031

PubMed Abstract | CrossRef Full Text | Google Scholar

Pittau, F., Megevand, P., Sheybani, L., Abela, E., Grouiller, F., Spinelli, L., et al. (2014b). Mapping epileptic activity: sources or networks for the clinicians? Front. Neurol. 5:218. doi: 10.3389/fneur.2014.00218

CrossRef Full Text | Google Scholar

Rathakrishnan, R., Moeller, F., Levan, P., Dubeau, F., and Gotman, J. (2010). BOLD signal changes preceding negative responses in EEG-fMRI in patients with focal epilepsy. Epilepsia 51, 1837–1845. doi: 10.1111/j.1528-1167.2010.02643.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakakibara, T., Sukigara, S., Otsuki, T., Takahashi, A., Kaneko, Y., Kaido, T., et al. (2012). Imbalance of interneuron distribution between neocortex and basal ganglia: consideration of epileptogenesis of focal cortical dysplasia. J. Neurol. Sci. 323, 128–133. doi: 10.1016/j.jns.2012.08.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Spreafico, R., Battaglia, G., Arcelli, P., Andermann, F., Dubeau, F., Palmini, A., et al. (1998). Cortical dysplasia: an immunocytochemical study of three patients. Neurology 50, 27–36. doi: 10.1212/WNL.50.1.27

PubMed Abstract | CrossRef Full Text | Google Scholar

Tassi, L., Colombo, N., Garbelli, R., Francione, S., Lo Russo, G., Mai, R., et al. (2002). Focal cortical dysplasia: neuropathological subtypes, EEG, neuroimaging and surgical outcome. Brain 125(Pt 8), 1719–1732. doi: 10.1093/brain/awf175

PubMed Abstract | CrossRef Full Text | Google Scholar

Tassi, L., Garbelli, R., Colombo, N., Bramerio, M., Russo, G. L., Mai, R., et al. (2012). Electroclinical, MRI and surgical outcomes in 100 epileptic patients with type II FCD. Epileptic Disord. 14, 257–266. doi: 10.1684/epd.2012.0525

PubMed Abstract | CrossRef Full Text | Google Scholar

Thornton, R., Vulliemoz, S., Rodionov, R., Carmichael, D. W., Chaudhary, U. J., Diehl, B., et al. (2011). Epileptic networks in focal cortical dysplasia revealed using electroencephalography-functional magnetic resonance imaging. Ann. Neurol. 70, 822–837. doi: 10.1002/ana.22535

PubMed Abstract | CrossRef Full Text | Google Scholar

Tyvaert, L., Hawco, C., Kobayashi, E., LeVan, P., Dubeau, F., and Gotman, J. (2008). Different structures involved during ictal and interictal epileptic activity in malformations of cortical development: an EEG-fMRI study. Brain 131(Pt 8), 2042–2060. doi: 10.1093/brain/awn145

PubMed Abstract | CrossRef Full Text | Google Scholar

Varotto, G., Tassi, L., Franceschetti, S., Spreafico, R., and Panzica, F. (2012). Epileptogenic networks of type II focal cortical dysplasia: a stereo-EEG study. Neuroimage 61, 591–598. doi: 10.1016/j.neuroimage.2012.03.090

PubMed Abstract | CrossRef Full Text | Google Scholar

Vulliemoz, S., Thornton, R., Rodionov, R., Carmichael, D. W., Guye, M., Lhatoo, S., et al. (2009). The spatio-temporal mapping of epileptic networks: combination of EEG-fMRI and EEG source imaging. Neuroimage 46, 834–843. doi: 10.1016/j.neuroimage.2009.01.070

PubMed Abstract | CrossRef Full Text | Google Scholar

Wagner, J., Urbach, H., Niehusmann, P., von Lehe, M., Elger, C. E., and Wellmer, J. (2011). Focal cortical dysplasia type IIb: completeness of cortical, not subcortical, resection is necessary for seizure freedom. Epilepsia 52, 1418–1424. doi: 10.1111/j.1528-1167.2011.03158.x

CrossRef Full Text | Google Scholar

Watanabe, S., An, D., Safi-Harb, M., Dubeau, F., and Gotman, J. (2014). Hemodynamic response function (HRF) in epilepsy patients with hippocampal sclerosis and focal cortical dysplasia. Brain Topogr. 27, 613–619. doi: 10.1007/s10548-014-0362-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Worsley, K. J., Liao, C. H., Aston, J., Petre, V., Duncan, G. H., Morales, F., et al. (2002). A general statistical analysis for fMRI data. Neuroimage 15, 1–15. doi: 10.1006/nimg.2001.0933

PubMed Abstract | CrossRef Full Text | Google Scholar

Zijlmans, M., Huiskamp, G., Hersevoort, M., Seppenwoolde, J. H., van Huffelen, A. C., and Leijten, F. S. (2007). EEG-fMRI in the preoperative work-up for epilepsy surgery. Brain 130(Pt 9), 2343–2353. doi: 10.1093/brain/awm141

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: EEG-fMRI, FCD, neuronal migration disorder, epileptic network, localization

Citation: Pittau F, Ferri L, Fahoum F, Dubeau F and Gotman J (2017) Contributions of EEG-fMRI to Assessing the Epileptogenicity of Focal Cortical Dysplasia. Front. Comput. Neurosci. 11:8. doi: 10.3389/fncom.2017.00008

Received: 01 December 2016; Accepted: 02 February 2017;
Published: 20 February 2017.

Edited by:

Daqing Guo, University of Electronic Science and Technology of China, China

Reviewed by:

Paolo Federico, University of Calgary, Canada
Xu Lei, Southwest University, China

Copyright © 2017 Pittau, Ferri, Fahoum, Dubeau and Gotman. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Francesca Pittau, francesca.pittau@gmail.com

Download