- 1Department of Neurology, University of Texas Health Science Center at Houston, Houston, TX, United States
- 2Epilepsy Center, University Hospitals Cleveland Medical Center, Cleveland, OH, United States
- 3Department of Neurology, Case Western Reserve University, Cleveland, OH, United States
- 4Department of Electrical Engineering and Computer Science, Case Western Reserve University, Cleveland, OH, United States
- 5Department of Neurology, University of Iowa School of Medicine, Iowa City, IA, United States
- 6NYU Langone School of Medicine, New York, NY, United States
- 7Sidney Kimmel Medical College, Thomas Jefferson University, Philadelphia, PA, United States
- 8Division of Pulmonary, Critical Care and Sleep Medicine, University Hospitals Medical Center, Cleveland, OH, United States
- 9Institute of Neurology, University College London, London, United Kingdom
- 10Feinberg School of Medicine, Northwestern University, Chicago, IL, United States
- 11Department of Neurobiology and the Brain Research Institute, University of California, Los Angeles, Los Angeles, CA, United States
- 12Department of Neurology, Columbia University, New York, NY, United States
- 13Mercy Health St. Mary's Campus, Grand Rapids, MI, United States
Introduction: Peri-ictal breathing dysfunction was proposed as a potential mechanism for SUDEP. We examined the incidence and risk factors for both ictal (ICA) and post-convulsive central apnea (PCCA) and their relationship with potential seizure severity biomarkers (i. e., post-ictal generalized EEG suppression (PGES) and recurrence.
Methods: Prospective, multi-center seizure monitoring study of autonomic, and breathing biomarkers of SUDEP in adults with intractable epilepsy and monitored seizures. Video EEG, thoraco-abdominal excursions, capillary oxygen saturation, and electrocardiography were analyzed. A subgroup analysis determined the incidences of recurrent ICA and PCCA in patients with ≥2 recorded seizures. We excluded status epilepticus and obscured/unavailable video. Central apnea (absence of thoracic-abdominal breathing movements) was defined as ≥1 missed breath, and ≥5 s. ICA referred to apnea preceding or occurring along with non-convulsive seizures (NCS) or apnea before generalized convulsive seizures (GCS).
Results: We analyzed 558 seizures in 218 patients (130 female); 321 seizures were NCS and 237 were GCS. ICA occurred in 180/487 (36.9%) seizures in 83/192 (43.2%) patients, all with focal epilepsy. Sleep state was related to presence of ICA [RR 1.33, CI 95% (1.08–1.64), p = 0.008] whereas extratemporal epilepsy was related to lower incidence of ICA [RR 0.58, CI 95% (0.37–0.90), p = 0.015]. ICA recurred in 45/60 (75%) patients. PCCA occurred in 41/228 (18%) of GCS in 30/134 (22.4%) patients, regardless of epilepsy type. Female sex [RR 11.30, CI 95% (4.50–28.34), p < 0.001] and ICA duration [RR 1.14 CI 95% (1.05–1.25), p = 0.001] were related to PCCA presence, whereas absence of PGES was related to absence of PCCA [0.27, CI 95% (0.16–0.47), p < 0.001]. PCCA duration was longer in males [HR 1.84, CI 95% (1.06–3.19), p = 0.003]. In 9/17 (52.9%) patients, PCCA was recurrent.
Conclusion: ICA incidence is almost twice the incidence of PCCA and is only seen in focal epilepsies, as opposed to PCCA, suggesting different pathophysiologies. ICA is likely to be a recurrent semiological phenomenon of cortical seizure discharge, whereas PCCA may be a reflection of brainstem dysfunction after GCS. Prolonged ICA or PCCA may, respectively, contribute to SUDEP, as evidenced by two cases we report. Further prospective cohort studies are needed to validate these hypotheses.
Introduction
Sudden Unexpected Death in Epilepsy (SUDEP) is the leading cause of premature mortality in patients with intractable epilepsy (1). The main SUDEP phenotype comprises frequent generalized convulsive seizures in patients with early onset, longstanding epilepsy (2, 3). Both cardiac and respiratory mechanisms likely contribute to SUDEP pathophysiology (3, 4), although video electroencephalogram monitored (VEEG) deaths suggest that terminal cardiac arrest is almost always preceded by central apnea (5). Central, obstructive and mixed apneas have all been proposed as SUDEP mechanisms, and may occur during or after seizures (6–9). Whereas, ictal central apnea (ICA) is common, prolonged ICA with profound oxygen desaturation may pose SUDEP risk, as may post-convulsive central apnea (PCCA) (10–14). The latter, when combined with bradycardia/asystole, comprised the majority of observed deaths in the MORTEMUS series, and two near-SUDEP instances in one observational SUDEP risk study (14). Thus, it is evident that breathing dysfunction plays a major role in SUDEP, although the exact characteristics of respiratory compromise that contribute to death are unknown (15). Since central apnea (prolonged ictal or post-convulsive) seems a viable, agonal mechanism, we set out to determine incidence, recurrence, and characteristics of peri-ictal central apnea. We assessed its influence on potential seizure severity biomarkers, such as hypoxemia extent and post-ictal generalized electroencephalographic (EEG) suppression (PGES) in a prospective study of SUDEP risk biomarkers. Further, we describe two additional cases of near-SUDEP due to prolonged, exaggerated peri-ictal central apnea.
Material And Methods
Patient Selection
All patients were prospectively recruited participants in the NINDS Center for SUDEP Research's Autonomic and Imaging Biomarkers of SUDEP multi-center project (U01-NS090407), and its preliminary phase, the Prevention and Risk Identification of SUDEP Mortality (PRISM) Project (P20NS076965). This study was carried out in accordance with the recommendations of University Hospitals Case Medical Center Institutions Review Boards (UHIRB) and University of Iowa, School of Medicine, Iowa City (IA) Institutions Review Boards. The protocol was approved by UHIRB and University of Iowa, School of Medicine, Iowa City (IA) Institutions Review Boards. All subjects gave written informed consent in accordance with the Declaration of Helsinki. Informed written consent was obtained for publication of two clinical cases. Patients with intractable epilepsy (failure of adequate trials of two or more antiepileptic medications) (16) aged ≥18 years underwent video-electroencephalography (VEEG) evaluation in the adult epilepsy monitoring units (EMU) of participating centers from September 2011 until April 2018. We included patients with recorded electroencephalographic seizures with or without clinical correlate, who had complete polygraphic physiological recordings sufficient for analysis. Exclusion criteria were status epilepticus (SE), obscured or unavailable video, and insufficient multimodal physiological recording quality.
Demographic and clinical data were collected. Semiology was classified into two types (17). (1) Generalized convulsive seizures (GCS) which included generalized tonic-clonic seizures, focal to bilateral tonic-clonic seizures, and focal onset motor bilateral clonic seizures. (2) Non-convulsive seizures (NCS), which included seizures with focal onset without evolution to bilateral tonic-clonic seizures, myoclonic seizures, absence seizures, and electroencephalographic seizures without clinical correlate. We determined state (awake or asleep) at seizure onset. We defined the putative epileptogenic zone based on clinical history, seizure semiology, neuroimaging, and VEEG.
Data Collection
All patients underwent prolonged surface VEEG monitoring using the 10–20 International Electrode System. EEG and electrocardiogram (EKG) were acquired using Nihon Kohden (Tokyo, Japan), Micromed (Modigliani Veneto, Italy), and Xltek (Natus) acquisition platforms. Peripheral capillary oxygen saturation (SpO2) was monitored using pulse oximetry (Nellcor OxiMax N-600x [Convidien], Masimo Radical-7 [Irvine], and SenTec Digital Monitoring System [Therwil BL]) and plethysmography (Ambu [Ballerup, Denmark] Sleepmate and Perfect Fit 2 [Dymedix]). Chest wall and abdominal excursions were recorded using inductance plethysmography (Ambu, Ballerup, Denmark and Sleepmate or Perfect Fit 2, Dymedix, St Paul, MN, USA).
Breathing analysis for apnea utilized careful composite analysis of inductance plethysmography, EEG breathing artifact and visually inspected thoraco-abdominal excursions 2 min before seizure onset (clinical or electrographic, whichever that occurred first). Central apnea (cessation of thoraco-abdominal breathing movements) was defined as ≥1 missed breath without any other explanation (i.e., speech, movement, or intervention), with a minimal duration of 5 s. ICA referred to apnea in NCS, or apnea occurring in the pre-convulsive phase of GCS. PCCA referred to apnea after GCS; we preferred this term to post-ictal central apnea since it could occur after convulsions but before EEG seizure end. Incidences of ICA and PCCA were determined, as well as their durations. Apnea could not be, and was not assessed during tonic or clonic movements, because of invariable artifact presence in breathing channels. A subgroup analysis identified recurrences of ICA and PCCA (in ≥2 seizures).
Hypoxemia was defined as SpO2 < 95% and where baseline SpO2 was already <95% a >1% drop was considered significant. Oxygen desaturations were classified as mild (SpO2 90–94%), moderate (75–89%), and severe (<75%). For SpO2 evaluation, several time points were considered. Firstly, SpO2 was determined at baseline, 2 min pre-ictally as mean SpO2 in a 15 s, artifact free epoch. For GCS and NCS, the overall desaturation nadir referred to the lowest SpO2 value registered during and up to 3 min after the seizure. To evaluate respiration in the pre-convulsive phase of GCS, an additional desaturation nadir was considered during this phase in patients with ICA.
Presence and duration of PGES (18) after GCS was determined by a validated automated EEG suppression detection tool (19), and supplemented with visual analysis by two epilepsy neurophysiologists when the tool gave no solution. Presence and duration of post-ictal EEG burst suppression were also determined. Combined PGES and burst suppression comprised EEG “recovery” duration.
Statistical Analysis
Statistical analyses were conducted using SPSS (version 24; IBM Corp, Armol, NY, USA) and SAS for Windows 9.4 (SAS Institute Inc., NC, USA). Summary statistics were reported as mean ± standard deviation (SD; median, range). Relative risk (RR) for the primary outcome of ICA and PCCA at a seizure level was assessed by Generalized Estimating Equation (GEE) with same subject exchangeable correlation. All variables with a p < 0.20 in a univariate analysis were included in a multivariate Poisson GEE regression (20). Variables related to ICA and PCCA durations were determined using Cox Regression with robust sandwich covariance estimation (21). Lastly, recurrence of ICA and PCCA for each patient with at least two seizures in this data were categorized as binary outcomes and patient-level covariates were included in respective logistic regressions. Corresponding 95% CIs of risk and hazard ratios were generated from these models. The significance level for hypothesis testing was set at p < 0.05.
Results
Demographic and Clinical Characteristics
Among 218 patients (130 female), 558 seizures met inclusion criteria. Four hundred and twenty-six seizures were previously reported in two different studies (13, 14).
Mean age at study was 40.2 years ± 14.7 (39; 18–77), mean epilepsy duration was 16.6 years ± 13.8 years (1 month−58 years) and mean age at epilepsy onset was 23.5 years ± 17.2 (20; 1 month-69 years). There were 321 NCS (in 128 patients) and 237 GCS (in 137 patients). State was sleep in 239 seizures and wakefulness in 318 seizures. One seizure arose from post-ictal coma in a patient with a seizure cluster.
There were 182 patients (493 seizures) with focal epilepsy and 33 patients (60 seizures) with generalized epilepsy. One patient had both focal and generalized epilepsy (2 seizures). Epilepsy type was unknown in 2 patients (3 seizures) (Table 1).
PGES could be determined in all but one GCS, where electrode artifact prevented accurate interpretation. PGES occurred in 165/236 (69.9%) GCS in 106/136 (77.9%) patients. Mean PGES duration was 38.9 s ± 21.2 (37; 1–169) and mean EEG recovery duration was 85 s ± 107.9 (54; 1–1,091).
Ictal Central Apnea (ICA) Incidence, Duration, and Recurrence
ICA could not be confidently ascertained in 71 seizures, all GCS, due to plethysmographic signal acquisition artifact. ICA occurred in 180/487 (36.9%) seizures in 83/192 (43.2%) patients: 65/166 (39.2%) in GCS and 115/321 (35.8%) in NCS (p = 0.960).
ICA preceded EEG seizure onset in 85/180 (47.2%) seizures by 9 s ± 9.5 (7; 1–58). ICA occurred after EEG seizure onset in 74/180 (41.1%) seizures, with a delay of 25.8 s ± 62 (8; 1–436). ICA coincided with EEG seizure onset in the remaining 21/180 (11.7%) seizures (Figure 1A).
ICA was the sole manifestation in 14/180 (7.8%) seizures. It coincided with clinical onset in 24/166 (14.5%) seizures, started before clinical onset in 102/166 (61.4%), preceding it by 16.3 ± 12.7 s (12; 1–66 s), and started after clinical onset in 40/180 (24.1%) seizures, with a difference of 26.4 s ± 36 (14.5; 1–195) (Figures 1B, 2A,B).
Figure 2. Example of ictal central apnea (ICA). Sensitivity 7 μV, High Frequency Filter: 70 Hz, Time constant: 0.1 s. (A) Twelve seconds after the electrographic onset, ICA is noted without any other clinical signs. (B) Continuation of ICA, with a total duration of 14 s, followed by breathing resumption. ABD, abdominal; EEG sz onset, electrographic seizure onset; EKG, electrocardiogram; ICA, Ictal central apnea; s, seconds; THOR, thoracic.
Information regarding nadir SpO2 in NCS and the pre-convulsive phase of GCS was available in 141/180 (78.3%) seizures with ICA, with a mean value of 87.7% ± 9.4 (91; 46–99).
All 180 seizures with ICA were seen in focal epilepsies, and none in patients with generalized epilepsy. In patients with focal epilepsies ICA was more frequent in temporal lobe epilepsies than extratemporal (p = 0.002) but there was no association with laterality (p = 0.215). ICA incidence did not show any differences regarding sex (p = 0.171) or epilepsy duration (p = 0.077) but it was related to older age at study (p = 0.004) and older age at epilepsy onset (p < 0.001). ICA was more frequent in seizures arising from sleep than during wakefulness (p = 0.013). ICA was not related to PGES or PCCA and did not affect PGES, EEG recovery, hypoxemia or PCCA durations or SpO2 nadir in GCS (p > 0.050) (Table 2).
After multivariate analysis, sleep state was related to presence of ICA [RR 1.326, CI95% (1.075–1.637), p = 0.008] whereas extratemporal epilepsy was related to lower incidence of ICA [RR 0.579, CI 95% (0.373–0.900), p = 0.015].
Mean ICA duration was 20.9 s ± 17.5 (14; 5–97 s) and was longer in patients with NCS without subsequent GCS than those with subsequent GCS [HR 2.276; CI 95% (1.565–3.311), p < 0.001] and in temporal lobe epilepsy compared to extratemporal [HR 1.753, CI 95% (1.065–2.885), p = 0.027]. ICA duration did not correlate with awake/sleep state at seizure onset (p > 0.050). Longer ICA duration was associated with lower SpO2 nadir during the NCS phase of seizures [HR 1.098, CI 95% (1.064–1.133), p < 0.001] and longer EEG recovery duration [HR 1.002, CI95% (1.001–1.003), p = 0.003]. ICA duration did not significantly correlate with hypoxemia duration, PGES duration, or SpO2 nadir in GCS.
Mean number of seizures per patient was 2.6 ± 1.7 (2; 1–8). In the total sample, 92/218 (42.2%) patients either had only one seizure, or only one analyzable seizure for ICA. The remaining 126/218 (57.8%) patients had recurrent seizures and comment on ICA could be made. ICA occurred in 60/126 (47.6%) patients, and recurred in 45/60 (75%) of the patients. No clinical characteristics (age, age at epilepsy onset, epilepsy duration, sex, and epileptogenic zone) were related to ICA recurrence (p > 0.05).
Clinical Case 1-Prolonged ICA and near-SUDEP
A 36 year-old right handed man with intractable right temporal lobe epilepsy of unknown etiology since the age of nine was enrolled into the study. His seizure semiology consisted of psychic aura followed by auditory aura with impaired awareness, and rare secondarily generalization. The last generalized convulsion had occurred 4 years before the admission. He had co-morbid depression. Previous antiepileptic drugs (AEDs) were carbamazepine, phenytoin, valproic acid and zonisamide. At admission for presurgical evaluation he was on oxcarbazepine 1,800 mg/day. Physical and neurological examinations were normal. Brain MRI was normal and the interictal PET scan showed bilateral mesial temporal hypometabolism, more pronounced on the right. Interictal recordings showed right temporal sharp waves (maximum at T8>F8). Retrospective review of older (non-study) VEEG records revealed a near-SUDEP incident (not included in the above analysis). The patient had 4 seizures during that admission. The first one, was a brief (<10 s duration auditory aura). The second and third seizures, were brief auras with rapid secondarily generalization, one arising from wakefulness and the other one arising from sleep. No comment about presence of ICA could be made on those seizures due to lack of plethysmography and rapid secondary generalization. No PCCA was noted in any of the GCS and regular breathing resumed immediately after clinical seizure end. These GCS occurred within 12 h of the fourth and last seizure. This was an apneic seizure with impaired awareness, and respiratory arrest lasting for 285 s, as evidenced by video analysis and oxygen desaturation. After a period of several shallow breaths, breathing finally resumed normally 311 s after seizure onset. Ictal EEG showed rhythmic alpha activity arising over the right antero-mesial temporal lobe with bilateral spread. No alteration in heart rhythm was noted apart from tachycardia. The patient was repositioned, oxygen administered, and ventilated with a face mask. He later underwent invasive evaluation, had further seizures without apnea, and a right temporal lobectomy in 2016 which resulted in seizure freedom (Engel Class 1, [>2 years]) Video 1.
Post-convulsive Central Apnea (PCCA) Incidence, Duration, and Recurrence
Presence of PCCA could not be confidently ascertained in 9/237 (3.8%) GCS in 3/137 (2.2%) patients due to movement artifact. PCCA was present in 41/228 (18%) of GCS in 30/134 (22.4%) patients.
In 24/41 (58.65%) seizures (in 19 patients), PCCA was observed without EEG seizure. In 14/41 (34.1%) seizures (in 12 patients), PCCA occurred with ongoing EEG seizure activity. In 3 seizures (in 3 patients) PCCA recurred in the same seizure, occurring initially with EEG seizure discharges and then after 1–2 breaths, recurring without accompanying seizure discharge(Figures 3A,B). In 13 seizures (in 11 patients), PCCA immediately followed clinical seizure end. In 25 seizures (in 20 patients), “delayed” PCCA occurred several breaths after clinical seizure end.
Figure 3. Example of post-convulsive central apnea (PCCA). Sensitivity 20 μV, High Frequency Filter: 70 Hz, Time constant: 0.1 s. (A). After the end of convulsive phase, the patient is apneic for 7 s. After the electrographic end, there are 2 noticeable breaths which are followed by another brief apnea. Lastly, regular breathing resumes. (B) Continuation from (A), regular breathing continues. ABD, abdominal; Clinical sz end, end of clinical seizure. EEG sz end, electrographic seizure end; EKG, electrocardiogram; PCCA, post-convulsive central apnea; s, seconds; THOR, thoracic.
PCCA was more frequent in women than in men (p = 0.004) and occurred more often in generalized than focal epilepsies (p = 0.016). In focal epilepsy, it was more frequently seen in extratemporal than temporal (p = 0.020) patients, but there was no relationship with lateralization (p = 0.148). PCCA was unrelated to age at study or epilepsy onset, epilepsy duration and awake or sleep states (p > 0.050). Whereas, ICA presence was not related to PCCA occurrence, PCCA was significantly associated with longer ICA duration (p = 0.001).
Presence of PCCA was not related to PGES duration, and was not associated with EEG recovery duration and total hypoxemia duration (p > 0.050). However, PCCA was associated with longer SpO2 recovery times to mild hypoxemia (>90%). [RR 1.01, CI95% (1.003–1.017), p = 0.003] (Table 3).
After multivariate regression analysis, female sex [RR 11.297, CI 95% (4.50–28.34), p < 0.001] and ICA duration [RR 1.149 CI 95% (1.053–1.254), p = 0.001] were related to PCCA, whereas absence of PGES was related to absence of PCCA [RR = 0.274, CI 95% (0.159–0.471), p < 0.001].
Mean PCCA duration was 8.9 s ± 4.9 (5–32). PCCA duration was longer in males [HR 1.844, CI 95% (1.06–3.19), p = 0.003]. Epilepsy type, awake/sleep state did not influence PCCA occurrence. PCCA duration did not correlate with age, epilepsy duration, PGES duration, EEG recovery, hypoxemia duration, or time to recovery to mild hypoxemia (p > 0.050).
Mean number of GCS per patient was 1.7 ± 1 (1; 1–5). In patients with GCS, 77/137 (56.2%) had only one GCS and the remaining 60/137 (47.8%) had two or more GCS. In the group of patients with repeated GCS, comment on PCCA could be made on 57/60 (95%) patients. PCCA was seen in 17/57 (29.8%). In 9/17 (52.9%) patients, PCCA was recurrent.
Clinical Case 2-Prolonged PCCA and near-SUDEP
A 15 year-old right handed girl with epilepsy since age 5 years was admitted for evaluation. She was not an enrolled study patient. Seizures occurred once or twice a month and lasted up to two with whole body sensory aura (tingling) followed by oral automatisms with impaired awareness. This was rarely followed by secondary generalized convulsions lasting 1–2 min. On several occasions, paramedics were summoned as an emergency because of cyanosis and unresponsiveness after generalized convulsions. On admission she was on lamotrigine 200 mg/day and levetiracetam 3,000 mg/day, having previously failed multiple other AEDS. She had no epilepsy risk factors and no family history of epilepsy. Her physical and neurological examinations were normal. Epilepsy protocol MRI brain scans were normal on two occasions. Inter-ictal brain FDG-PET showed focal hypometabolism in the anterior left temporal lobe tip. Non–invasive VEEG monitoring showed left temporal sharp waves, maximum at F7/T7/FT9. Four habitual clinical seizures were recorded without secondarily generalization. EEG onsets were left hemispheric but not further localizable.
She underwent invasive EEG monitoring for better localization of the epileptogenic zone. A left subdural grid (8 × 6) was implanted along with strips covering the left orbitofrontal, superior temporal, inferior temporal regions, as well as left anterior-anterior, anterior-middle and anterior-posterior temporal, left middle temporal, left middle-middle, and middle-posterior. A left anterior temporal seizure was recorded, with typical automatisms and impaired awareness, right face clonic movements, and a secondary generalized tonic clonic seizure. After clinical seizure end, the patient was immediately apneic (as evidenced by video analysis, cyanosis, and severe O2 desaturation) for 126 s, followed by an isolated breath. A second period of apnea/hypopnea was then seen until regular breathing pattern resumed a total of 187 s after clinical seizure end. EEG seizure discharges were seen up to 25 s after clinical seizure end. Thirty nine seconds after clinical seizure end, there was concurrent progressive bradycardia followed by 10 s of asystole. Cardiac rhythm resumed, with a pattern of bradycardia and normal sinus rhythm, for 75 s, after which EKG signal was lost, but pulse artifact was evident on EEG. EEG suppression duration (all invasive electrodes), was ~254 s. During the episode, there was repeated tactile nursing intervention. Further, her head was re-positioned and O2 administered. No active resuscitation measures were performed. Due to continuing seizures, the patient underwent responsive neurostimulation (RNS® System) and was temporarily seizure free for 3 years, until recent recurrence of focal seizures at last follow up Video 2.
Discussion
In this study we summarize incidence and risk factors for both ICA and PCCA. Additionally, we describe two near-SUDEP instances of prolonged ICA and PCCA, respectively.
ICA incidence in our study (43.2%), on a larger number of patients, was similar to those previously reported (10, 13). Consistent with our previous reports, ICA was not observed in patients with generalized epilepsy, and was more frequent in patients with temporal rather than extratemporal epilepsy (13, 14). ICA preceded other clinical signs in the vast majority of seizures and in almost half of them, also preceded EEG seizure onset. The observation of ICA being an exclusive feature for focal epilepsies, and especially in those from the temporal lobe, is consistent with previous human stimulation studies pointing out the amygdala, hippocampus, and mesial temporal pole, regardless of lateralization, as the symptomatogenic zone for ICA (11, 13, 22). The absence of ICA in the GCS of generalized epilepsy is in large part due to immediate or rapid generalization where breathing compromise may be partly or wholly due to generalized muscle tonicity that includes respiratory musculature, rather than solely due to unequivocal central apnea. However, we cannot be sure that these patients truly do not have ICA.
Unlike previous publications (13), where no differences between awake/sleep states at seizure onset were found, in this more robustly powered study, ICA occurred more frequently in seizures arising from sleep. One possible explanation is that ICA is easier to detect in the sleep state, where acquisition artifact is less prevalent. However, awake recordings were not disproportionately excluded because of artifact and such disparities can be explained by physiological differences in breathing control during sleep and wakefulness (23). Breathing is under automatic control through multiple pontomedullary nuclei, the pre-Bötzinger complex (pre-BötC) comprising the main rhythm generator (24). Cortical and subcortical structures, such as thalamus, hypothalamus, amygdalo-hippocampal complex, cerebellum, and mesencephalic nuclei relay to pontomedullary respiratory centers and along with peripheral sensory feedback, modulate breathing output (25). There is increased evidence that serotoninergic neurons lying in the midline raphe nuclei play an active role in both arousal and chemoreception (26, 27). These neurons tonically excitate multiple components of the brainstem respiratory network, with interconnections with the pre-BötC, and also act as central chemosensors, detecting changes in tissue CO2/H+ modulating the aforementioned tonic excitatory drive to adjust ventilation accordingly (28). Mice with genetically deleted medullary serotoninergic neurons lack any arousal response to inhalation of CO2 but have normal arousal responses to other stimuli such as hypoxia, sound and air puffs (26). Moreover, the activity of medullary raphe serotoninergic neurons is highest during wakefulness and absent during REM sleep (28). Most SUDEP cases occur at night, and in the MORTEMUS study this was true in the majority of monitored cases (5) animal studies reveal brainstem serotoninergic dysfunction during and after seizures, with decreased firing of the medullary raphe neurons during the ictal and post-ictal periods (29). Therefore, in the setting of a potentially dysfunctional serotoninergic network in epileptic patients, sleep would constitute a vulnerable period for breathing disturbances (30).
The near-SUDEP case with prolonged ICA, in whom breathing resumed after seizure end raises a number of interesting issues. First, it seems possible that unobserved and in the absence of active intervention, the outcome could have been fatal. This supports the contention that prolonged ICA is dangerous and potentially lethal. Second, the near-fatal seizure episode may be more akin to focal status epilepticus, and supports the view that the latter cannot be excluded in at-home, unobserved deaths that are usually labeled SUDEP (31). Revisiting SUDEP definition (which mandates exclusion of status epilepticus) may be necessary. The episode mimics a sheep model of status epilepticus and apneic death (32). Third, seizure termination and resumption of breathing with seizure end, suggest that ICA is driven by seizure discharge (likely in the mesial temporal structures), rather than other mechanisms (13). Finally, why the majority of ICA instance are self-terminating and some become prolonged, is unresolved, but may reflect the consequences of damage caused by early onset, long-standing epilepsy, and frequent GCS to key breathing control sites (amygdala, hippocampus, dorsal thalamus, anterior cingulate, ventrolateral medulla etc.) (33–35), rendering greater susceptibility to exaggerated apneic responses.
PCCA incidence (22.4%) was almost half of ICA incidence, more frequently observed in female subjects (although more likely to be longer in male patients), and commoner in those with longer ICA duration. In contrast to ICA, PCCA was observed in both patients with focal, and generalized epilepsy. These differences suggest differing pathophysiologies (14). Whereas, ICA appears to be a semiological phenomenon most often resulting from seizure activity in the amygdalo-hippocampal complex, PCCA most likely results from seizure spread to the brainstem during GCS, regardless of epileptogenic zone (36). Breathing cessation may be derived either from active depolarization and activation of crucial breathing centers that generate apneic responses, such as the periaqueductal gray (37, 38), or disruption of the normal functioning of rhythm-generating neurons and its intricate network, leading to breathing cessation (39). This is consistent with animal models of SUDEP showing post-ictal depolarization in dorsal medulla (40), in which apnea and PGES precede cardiac arrest, and resemble the clinical phenotype of monitored SUDEP patients in the MORTEMUS study (5). Human neuroimaging and neuropathological studies have shown damage in key brainstem structures that modulate breathing, such as the medullary raphe and ventrolateral medulla, in SUDEP and high SUDEP-risk patients (33, 41, 42).
In our study, ICA duration conferred higher risk for PCCA and was related to lower SpO2 nadirs in the pre-convulsive phase. This may indicate greater hypoxemia induced brainstem compromise, leading to PCCA. Seizure induced focal brainstem hypoxia, due to vasospasm, has been described in animal models and posited as a potential mechanism for SUDEP (43–45). Moreover, we found that the presence of PCCA was related to longer hypoxemia recovery times (SpO2 > 90%). Although causality could not be established, longer hypoxemia recovery times may be a consequence of PCCA. However, PCCA duration itself, was not related to hypoxemia severity or duration. PCCA durations were typically short, and hypoxemia severity may be more related to GCS severity, although PCCA impact on hypoxemia may become independently important in instances of prolonged PCCA.
Although PCCA occurred preferentially in women, its duration was longer in men, consistent with the SUDEP phenotype (2). Sex differences in breathing function and the protective role of estrogens in respiratory diseases may explain these findings (46–48), since differences in epilepsy phenotype do not explain these differences. As with ICA, duration rather than presence, may primarily influence SUDEP risk in PCCA.
PGES is a frequent finding in GCS, particularly in those arising from sleep and is related to the symmetric tonic phase, postictal immobility, lack of early oxygen administration, duration of oxygen desaturation and lower SpO2 nadir values (49–51). PGES has been postulated as a SUDEP biomarker, especially if prolonged (>50 s) (18). Its relationship with ICA or PCCA has not been established, except indirectly through O2 desaturation findings (50). In our study, PCCA was proportionally seen more frequently in seizures with PGES than in seizures without. The pathogenesis of PGES is not well-determined. Cortical neuronal exhaustion or a disruption of ascending inputs after a GCS, or a combination of both, are viable hypotheses. Disruption of ascending pathways such as from the reticular activating system may conceivably prolong the comatose post-ictal state, as well as modulate cortical neuronal activity, and thus, impair the protective behavioral effect of arousal to overcome PCCA (52).
ICA recurred in 75% of the cases, whereas PCCA recurred in 52.9%. This may further reinforce that ICA is a semiological, and therefore recurrent, phenomenon. However, PCCA seems only slightly less frequent, and may also be semiological rather than probabilistic, although our two case reports suggest that prolonged instances of either, are what potentially determine mortality risk. PCCA instances combined with bradycardia/asystole may be particularly dangerous. Our second clinical case of prolonged PCCA accompanied by asystole, resembles the clinical phenotype described in the MORTEMUS study and in a recent analysis of this cohort in a smaller number of patients (14). Invasive monitoring did not show ongoing seizure activity that was concurrent with apnea, reinforcing once again, the different pathophysiologies of ICA, and PCCA, with higher likelihood of involvement of subcortical structures, such as the brainstem, in PCCA (53, 54).
Our study has several limitations. First, it is an observational study in a select group of patients (i.e., primarily treatment resistant epilepsy) and does not necessarily reflect seizure phenomenology or SUDEP risk in a treatment responsive population. Detection of ICA was heavily dependent on extent of acquisition artifact, and hence we may have underestimated incidence. Alternatively, postictal immobility after GCS allowed PCCA identification in the majority of cases. Breathing analysis through polygraphic study was limited to thoraco-abdominal movement and pulse oximetry. Thus, additional information on the presence of mixed central/obstructive apneas, is unavailable. Our apnea definition differs from previous literature, which is based on the 10 s sleep-study criterion, and therefore, ICA and PCCA incidence may be overestimated in our study. However, based on our brain stimulation experiments on breathing modulation, where brief stimulation periods result in immediate and brief (<10 s) apneas, we believe our definition is both accurate and sensitive (12, 22). Our conclusions are based on a relatively small group of number seizures in the primary generalized epilepsy group compared to patients with focal epilepsy. Lastly, our study was based on surface EEG and persistence of intracranial seizure in deep, apnea causing structures (12) cannot be completely excluded in patients with PCCA. However, Case 2 above, along with a previous case reported in literature, suggest that apnea in epileptic patients can occur in the absence of electrographic seizure (53).
Conclusions
Peri-ictal central apnea takes two main forms, as ICA or PCCA. ICA incidence is almost twice PCCA incidence and is only seen in focal epilepsies, suggesting different pathophysiologies. Both ICA and PCCA may be recurrent, but prolonged instances leading to SUDEP and near-SUDEP may be probabilistic instances. Prolonged ICA is related to presence of PCCA, possibly due to greater effect of ICA-induced hypoxemia on brainstem function. Absence of PCCA is associated with absence of PGES, suggesting that PCCA presence directly correlates with GCS severity. Alternatively, brainstem structures responsible for arousal and breathing may obviate PGES occurrence. Apnea preceding both EEG as well as clinical seizure onset in a substantial number of patients suggests that plethysmographic respiratory monitoring in regular clinical practice may have seizure detection value. Moreover, such monitoring may facilitate detection of prolonged ICA and PCCA, thus allowing SUDEP risk quantification although further evidence is required to confirm this. Both prolonged ICA and PCCA may contribute to SUDEP. Further prospective cohort studies are needed to validate these hypotheses.
Data Availability
The datasets analyzed in this study are available from the corresponding author on request.
Author Contributions
LV had a major role in the acquisition and analysis of data, interpreted the data, drafted the manuscript for intellectual content. NL designed and conceptualized study, interpreted the data, revised the manuscript for intellectual content. JH analysed the data, performed statistical analysis, images and video editing. MR acted as a recruiter, revised the manuscript for intellectual content. KL analysed the data. RS, DF, MN, BG, SS, JO, RH, BD, LB, OD, and GR revised the manuscript for intellectual content. KS interpreted the data, revised the manuscript for intellectual content. LA, CS, BZ, NH, NS, XZ, VR-M, and AT performed data acquisition. AZ led and coordinated communication among sites; revised the manuscript for intellectual content. CT analysed the data and performed statistical analysis. SL designed and conceptualized study, performed analysis and interpretation of data, revised the manuscript for intellectual content.
Funding
LB, DF, MN, RH, BD, LV, MR, BG, NH, AZ, RS, JO, LA, and OD are funded by the Center for SUDEP Research: NIH/NINDS U01-NS090407. SS is on the speaker's bureau for Sunovion and Eisai. He receives grant support from NINDS (RFA-NS-14-004), NIH/NINDS U01-NS090407, and the Danny Did Foundation. SL is funded by the Center for SUDEP Research: NIH/NINDS U01-NS090405 and NIH/NINDS U01-NS090407. GR is funded by the Center for SUDEP Research: NIH/NINDS U01-NS090414 and NIH/NINDS U01-NS090407.
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
We would like to thank patients and personnel from the different monitoring units for their invaluable contribution to understanding epilepsy and SUDEP pathophysiology.
Supplementary Material
The Supplementary Material for this article can be found online at: https://doi.org/10.6084/m9.figshare.7188998
References
1. Tomson T, Nashef L, Ryvlin P. Sudden unexpected death in epilepsy: current knowledge and future directions. Lancet Neurol. (2008) 7:1021–31. doi: 10.1016/S1474-4422(08)70202-3
2. Hesdorffer DC, Tomson T, Benn E, Sander JW, Nilsson L, Langan Y, et al. Combined analysis of risk factors for SUDEP. Epilepsia. (2011) 52:1150–9. doi: 10.1111/j.1528-1167.2010.02952.x
3. Devinsky O, Hesdorffer DC, Thurman DJ, Lhatoo S, Richerson G. Sudden unexpected death in epilepsy: epidemiology, mechanisms, and prevention. Lancet Neurol. (2016) 15:1075–88. doi: 10.1016/S1474-4422(16)30158-2
4. Nascimento FA, Tseng ZH, Palmiere C, Maleszewski JJ, Shiomi T, McCrillis A, et al. Pulmonary and cardiac pathology in sudden unexpected death in epilepsy (SUDEP). Epilepsy Behav. (2017) 73:119–25. doi: 10.1016/j.yebeh.2017.05.013
5. Ryvlin P, Nashef L, Lhatoo SD, Bateman LM, Bird J, Bleasel A, et al. Incidence and mechanisms of cardiorespiratory arrests in epilepsy monitoring units (MORTEMUS): a retrospective study. Lancet Neurol. (2013) 12:966–77. doi: 10.1016/S1474-4422(13)70214-X
6. Tavee J, Morris H III. Severe postictal laryngospasm as a potential mechanism for sudden unexpected death in epilepsy: a near-miss in an EMU. Epilepsia. (2008) 49:2113–7. doi: 10.1111/j.1528-1167.2008.01781.x
7. Nakase K, Kollmar R, Lazar J, Arjomandi H, Sundaram K, Silverman J, et al. Laryngospasm, central and obstructive apnea during seizures: defining pathophysiology for sudden death in a rat model. Epilepsy Res. (2016) 128:126–39. doi: 10.1016/j.eplepsyres.2016.08.004
8. Amir J, Ashkenazi S, Schonfeld T, Weitz R, Nitzan M. Laryngospasm as a single manifestation of epilepsy. Arch Dis Child. (1983) 58:151–3. doi: 10.1136/adc.58.2.151
9. Lacuey N, Vilella L, Hampson JP, Sahadevan J, Lhatoo SD. Ictal laryngospasm monitored by video-EEG and polygraphy: a potential SUDEP mechanism. Epileptic Disord. (2018) 20:146–50. doi: 10.1684/epd.2018.096
10. Bateman LM, Li CS, Seyal M. Ictal hypoxemia in localization-related epilepsy: analysis of incidence, severity and risk factors. Brain. (2008) 131:3239–45. doi: 10.1093/brain/awn277
11. Dlouhy BJ, Gehlbach BK, Kreple CJ, Kawasaki H, Oya H, Buzza C, et al. Breathing Inhibited When Seizures Spread to the Amygdala and upon Amygdala Stimulation. J Neurosci. (2015) 35:10281–9. doi: 10.1523/JNEUROSCI.0888-15.2015
12. Lacuey N, Zonjy B, Londono L, Lhatoo SD. Amygdala and hippocampus are symptomatogenic zones for central apneic seizures. Neurology. (2017) 88:701–5. doi: 10.1212/WNL.0000000000003613
13. Lacuey N, Zonjy B, Hampson JP, Rani MRS, Zaremba A, Sainju RK, et al. The incidence and significance of periictal apnea in epileptic seizures. Epilepsia. (2018) 59:573–82. doi: 10.1111/epi.14006
14. Vilella L, Lacuey N, Hampson JP, et al. Post-convulsive central apnea as a biomarker for sudden unexpected death in epilepsy (SUDEP). Neurology. (2019) 92:e171–82. doi: 10.1212/WNL.0000000000006785
15. Bateman LM, Spitz M, Seyal M. Ictal hypoventilation contributes to cardiac arrhythmia and SUDEP: report on two deaths in video-EEG-monitored patients. Epilepsia. (2010) 51:916–20. doi: 10.1111/j.1528-1167.2009.02513.x
16. Kwan P, Arzimanoglou A, Berg AT, Brodie MJ, Allen Hauser W, Mathern G, et al. Definition of drug resistant epilepsy: consensus proposal by the ad hoc task force of the ILAE commission on therapeutic strategies. Epilepsia. (2010) 51:1069–77. doi: 10.1111/j.1528-1167.2009.02397.x
17. Fisher RS, Cross JH, French JA, Higurashi N, Hirsch E, Jansen FE, et al. Operational classification of seizure types by the international league against epilepsy: position paper of the ILAE commission for classification and terminology. Epilepsia. (2017) 58:522–30. doi: 10.1111/epi.13670
18. Lhatoo SD, Faulkner HJ, Dembny K, Trippick K, Johnson C, Bird JM. An electroclinical case-control study of sudden unexpected death in epilepsy. Ann Neurol. (2010) 68:787–96. doi: 10.1002/ana.22101
19. Theeranaew W, McDonald J, Zonjy B, Kaffashi F, Moseley BD, Friedman D, et al. Automated detection of postictal generalized EEG suppression. IEEE Trans Biomed Eng. (2018) 65:371–7. doi: 10.1109/TBME.2017.2771468
20. Zou G. A modified poisson regression approach to prospective studies with binary data. Am J Epidemiol. (2004) 159:702–6. doi: 10.1093/aje/kwh090
21. Gharibvand L, Liu L. Analysis of Survival data with Clustered Events [online]. Available online at: https://support.sas.com/resources/papers/proceedings09/237-2009.pdf.
22. Lacuey N, Hampson JP, Harper RM, Miller JP, Lhatoo SD. Limbic and paralimbic structures driving ictal central apnea. Neurology. (2019) 11:10.1212/WNL.0000000000006920. doi: 10.1212/WNL.0000000000006920
23. Zaidi S, Gandhi J, Vatsia S, Smith NL, Khan SA. Congenital central hypoventilation syndrome: an overview of etiopathogenesis, associated pathologies, clinical presentation, and management. Auton Neurosci. (2018) 210:1–9. doi: 10.1016/j.autneu.2017.11.003
24. Feldman JL, Del Negro CA. Looking for inspiration: new perspectives on respiratory rhythm. Nat Rev Neurosci. (2006) 7:232–42. doi: 10.1038/nrn1871
25. Hilaire G, Pasaro R. Genesis and control of the respiratory rhythm in adult mammals. News Physiol Sci. (2003) 18:23–8. doi: 10.1152/nips.01406.2002
26. Buchanan GF, Richerson GB. Central serotonin neurons are required for arousal to CO2. Proc Natl Acad Sci USA. (2010) 107:16354–9. doi: 10.1073/pnas.1004587107
27. Richerson GB. Serotonergic neurons as carbon dioxide sensors that maintain pH homeostasis. Nat Rev Neurosci. (2004) 5:449–61. doi: 10.1038/nrn1409
28. Benarroch EE. Medullary serotonergic system: organization, effects, and clinical correlations. Neurology. (2014) 83:1104–11. doi: 10.1212/WNL.0000000000000806
29. Zhan Q, Buchanan GF, Motelow JE, Andrews J, Vitkovskiy P, Chen WC, et al. Impaired serotonergic brainstem function during and after seizures. J Neurosci. (2016) 36:2711–22. doi: 10.1523/JNEUROSCI.4331-15.2016
30. Buchanan GF, Murray NM, Hajek MA, Richerson GB. Serotonin neurones have anti-convulsant effects and reduce seizure-induced mortality. J Physiol. (2014) 592:4395–410. doi: 10.1113/jphysiol.2014.277574
31. Devinsky O, Spruill T, Thurman D, Friedman D. Recognizing and preventing epilepsy-related mortality: a call for action. Neurology. (2016) 86:779–86. doi: 10.1212/WNL.0000000000002253
32. Johnston SC, Siedenberg R, Min JK, Jerome EH, Laxer KD. Central apnea and acute cardiac ischemia in a sheep model of epileptic sudden death. Ann Neurol. (1997) 42:588–94. doi: 10.1002/ana.410420409
33. Patodia S, Somani A, O'Hare M, Venkateswaran R, Liu J, Michalak Z, et al. The ventrolateral medulla and medullary raphe in sudden unexpected death in epilepsy. Brain. (2018) 141:1719–33. doi: 10.1093/brain/awy078
34. Allen LA, Harper RM, Kumar R, Guye M, Ogren JA, Lhatoo SD, et al. Dysfunctional brain networking among autonomic regulatory structures in temporal lobe epilepsy patients at high risk of sudden unexpected death in epilepsy. Front Neurol. (2017) 8:544. doi: 10.3389/fneur.2017.00544
35. Ogren JA, Tripathi R, Macey PM, Kumar R, Stern JM, Eliashiv DS, et al. Regional cortical thickness changes accompanying generalized tonic-clonic seizures. Neuroimage Clin. (2018) 20:205–15. doi: 10.1016/j.nicl.2018.07.015
36. Salam MT, Montandon G, Genov R, Devinsky O, Del Campo M, Carlen PL. Mortality with brainstem seizures from focal 4-aminopyridine-induced recurrent hippocampal seizures. Epilepsia. (2017) 58:1637–44. doi: 10.1111/epi.13846
37. Subramanian HH. Descending control of the respiratory neuronal network by the midbrain periaqueductal grey in the rat in vivo. J Physiol. (2013) 591:109–22. doi: 10.1113/jphysiol.2012.245217
38. Subramanian HH, Balnave RJ, Holstege G. The midbrain periaqueductal gray control of respiration. J Neurosci. (2008) 28:12274–83. doi: 10.1523/JNEUROSCI.4168-08.2008
39. Faingold CL. Locomotor behaviors in generalized convulsions are hierarchically driven from specific brain-stem nuclei in the network subserving audiogenic seizure. Ann N Y Acad Sci. (1998) 860:566–9.
40. Aiba I, Noebels JL. Spreading depolarization in the brainstem mediates sudden cardiorespiratory arrest in mouse SUDEP models. Sci Transl Med. (2015) 7:282ra246. doi: 10.1126/scitranslmed.aaa4050
41. Wandschneider B, Koepp M, Scott C, Micallef C, Balestrini S, Sisodiya SM, et al. Structural imaging biomarkers of sudden unexpected death in epilepsy. Brain. (2015) 138:2907–19. doi: 10.1093/brain/awv233
42. Mueller SG, Nei M, Bateman LM. Brainstem network disruption: a pathway to sudden unexplained death in epilepsy? Hum Brain Mapp. (2018) 39:4820–30. doi: 10.1002/hbm.24325
43. Farrell JS, Gaxiola-Valdez I, Wolff MD, David LS, Dika HI, Geeraert BL, et al. Postictal behavioural impairments are due to a severe prolonged hypoperfusion/hypoxia event that is COX-2 dependent. Elife. (2016) 5:e19352. doi: 10.7554/eLife.19352
44. Farrell JS, Colangeli R, Wolff MD, Wall AK, Phillips TJ, George A, et al. Postictal hypoperfusion/hypoxia provides the foundation for a unified theory of seizure-induced brain abnormalities and behavioral dysfunction. Epilepsia. (2017) 58:1493–501. doi: 10.1111/epi.13827
45. Wall AK. Seizure-Induced Brainstem Hypoxia as a Possible Mechanism of Sudden Unexpected Death in Epilepsy. Alberta: University of Calgary (2017).
46. LoMauro A, Aliverti A. Sex differences in respiratory function. Breathe. (2018) 14:131–40. doi: 10.1183/20734735.000318
47. Townsend EA, Miller VM, Prakash YS. Sex differences and sex steroids in lung health and disease. Endocr Rev. (2012) 33:1–47. doi: 10.1210/er.2010-0031
48. Macey PM, Prasad JP, Ogren JA, Moiyadi AS, Aysola RS, Kumar R, et al. Sex-specific hippocampus volume changes in obstructive sleep apnea. Neuroimage Clin. (2018) 20:305–17. doi: 10.1016/j.nicl.2018.07.027
49. Alexandre V, Mercedes B, Valton L, Maillard L, Bartolomei F, Szurhaj W, et al. Risk factors of postictal generalized EEG suppression in generalized convulsive seizures. Neurology. (2015) 85:1598–603. doi: 10.1212/WNL.0000000000001949
50. Peng W, Danison JL, Seyal M. Postictal generalized EEG suppression and respiratory dysfunction following generalized tonic-clonic seizures in sleep and wakefulness. Epilepsia. (2017) 58:1409–14. doi: 10.1111/epi.13805
51. Kuo J, Zhao W, Li CS, Kennedy JD, Seyal M. Postictal immobility and generalized EEG suppression are associated with the severity of respiratory dysfunction. Epilepsia. (2016) 57:412–7. doi: 10.1111/epi.13312
52. Dlouhy BJ, Gehlbach BK, Richerson GB. Sudden unexpected death in epilepsy: basic mechanisms and clinical implications for prevention. J Neurol Neurosurg Psychiatry. (2016) 87:402–13. doi: 10.1136/jnnp-2013-307442
53. Ba-Armah DM, Donner EJ, Ochi A. “Saved by the Bell”: near SUDEP during intracranial EEG monitoring. Epilepsia Open. (2018) 3:98–102. doi: 10.1002/epi4.12093
Keywords: apnea, breathing, epilepsy, ictal central apnea (ICA), seizures, sudden unexpected death in epilepsy (SUDEP), post-convulsive central apnea (PCCA)
Citation: Vilella L, Lacuey N, Hampson JP, Rani MRS, Loparo K, Sainju RK, Friedman D, Nei M, Strohl K, Allen L, Scott C, Gehlbach BK, Zonjy B, Hupp NJ, Zaremba A, Shafiabadi N, Zhao X, Reick-Mitrisin V, Schuele S, Ogren J, Harper RM, Diehl B, Bateman LM, Devinsky O, Richerson GB, Tanner A, Tatsuoka C and Lhatoo SD (2019) Incidence, Recurrence, and Risk Factors for Peri-ictal Central Apnea and Sudden Unexpected Death in Epilepsy. Front. Neurol. 10:166. doi: 10.3389/fneur.2019.00166
Received: 02 October 2018; Accepted: 08 February 2019;
Published: 01 March 2019.
Edited by:
Andrea Romigi, Mediterranean Neurological Institute (IRCCS), ItalyReviewed by:
Jose F. Tellez-Zenteno, University of Saskatchewan, CanadaMarino M. Bianchin, Universidade Federal do Rio Grande do Sul (UFRGS), Brazil
Copyright © 2019 Vilella, Lacuey, Hampson, Rani, Loparo, Sainju, Friedman, Nei, Strohl, Allen, Scott, Gehlbach, Zonjy, Hupp, Zaremba, Shafiabadi, Zhao, Reick-Mitrisin, Schuele, Ogren, Harper, Diehl, Bateman, Devinsky, Richerson, Tanner, Tatsuoka and Lhatoo. 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) and the copyright owner(s) 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: Laura Vilella, lvilellabertran@gmail.com