ABSTRACT: Some of the most clinically consequential aspects of focal epilepsy, e.g. loss of consciousness, arise from the generalization or propagation of seizures through local and large-scale neocortical networks. Yet, the dynamics of such neocortical propagation remain poorly understood. Here, we studied the microdynamics of focal seizure propagation in neocortical patches (4×4 mm) recorded via high-density microelectrode arrays (MEAs) implanted in people with pharmacologically resistant epilepsy. Our main findings are threefold: (1) a newly developed stage segmentation method, applied to local field potentials (LFPs) and multiunit activity (MUA), revealed a succession of discrete seizure stages, each lasting several seconds. These different stages showed characteristic evolutions in overall activity and spatial patterns, which were relatively consistent across seizures within each of the 5 patients studied. Interestingly, segmented seizure stages based on LFPs or MUA showed a dissociation of their spatiotemporal dynamics, likely reflecting different contributions of non-local synaptic inputs and local network activity. (2) As previously reported, some of the seizures showed a peak in MUA that happened several seconds after local seizure onset and slowly propagated across the MEA. However, other seizures had a more complex structure characterized by, for example, several MUA peaks, more consistent with the succession of discrete stages than the slow propagation of a simple wavefront of increased MUA. In both cases, nevertheless, seizures characterized by spike-wave discharges (SWDs, ~2-3 Hz) eventually evolved into patterns of phase-locked MUA and LFPs. (3) Individual SWDs or gamma oscillation cycles (25-60 Hz), characteristic of two different types of recorded seizures, tended to propagate with varying degrees of directionality, directions of propagation and speeds, depending on the identified seizure stage. However, no clear relationship was observed between the MUA peak onset time (in seizures where such peak onset occurred) and changes in MUA or LFP propagation patterns. Overall, our findings indicate that the recruitment of neocortical territories into ictal activity undergoes complex spatiotemporal dynamics evolving in slow discrete states, which are consistent across seizures within each patient. Furthermore, ictal states at finer spatiotemporal scales (individual SWDs or gamma oscillations) are organized by slower time scale network dynamics evolving through these discrete stages.
