First experiences of whole-head OP-MEG recordings from a patient with Epilepsy

Optically Pumped Magnetometer based Magnetoencephalography (OP-MEG) has significant potential for clinical use in pre-surgical planning in epilepsy. Unlike current clinical MEG, the sensors do not require cryogenic cooling and so can be placed directly on the patient's head. This allows the patient to move during the recording and means that the sensor positions can be chosen to suit the patient's head shape and suspected epileptogenic focus. However, OP-MEG is a new technology, and more work is needed to demonstrate this potential. We present OP-MEG recordings from a patient (male in their 30s) with a radiologically identifiable focal cortical dysplasia (FCD) in the right superior frontal sulcus approximately 1.9 cm3 in volume. Previous scalp EEG studies and prolonged video-EEG telemetry did not identify any interictal epileptiform abnormalities. We recorded 30 minutes of OP-MEG with 62-channel, whole-head sensor coverage. During the experiment, the patient's head was unconstrained. We localised interictal epileptiform discharges (IEDs) from this recording with a beamformer and by fitting a dipole to the averaged IED data. Both the beamformer peak and dipole fit locations were within 2.3 cm of the MRI lesion boundary. This single subject, proof-of-concept recording provides further evidence that OP-MEG can be a useful and minimally invasive tool in the clinical evaluation of epilepsy.

and predict seizure freedom after surgery (Englot et al., 2015). In a recent retrospective study of 1000 35 patients (Rampp et al., 2019), it was shown that complete resection of MEG localisations was 36 associated with significantly higher probability of Engel 1 outcome (free from disabling seizures) 37 over the 10 years proceeding surgery. 38 However, MEG is often not clinically available. A survey of epilepsy surgery centres found only 1/3 39 of the centres had access to MEG (Mouthaan et al., 2016). Optically Pumped magnetometer-based 40 MEG (OP-MEG) has many advantages which could make it more accessible in a clinical setting than 41 the currently more standard superconducting quantum interference device (SQUID) based MEG. The 42 sensors do not need to be cryogenically cooled and can be worn directly on the head, meaning that 43 there is scope for the patient to move. This is particularly important for less compliant patient groups 44 such as children (Wehner et al., 2008;Larson and Taulu, 2017) or patients requiring long-term 45 monitoring for the purpose of capturing ictal events. Additionally, the sensors can be placed directly 46 on the scalp, rather than the one-size-fits-all array which is typical for SQUID-MEG. This means that 47 it is possible to capture maximal signal for any head-size. This can increase the SNR of observed 48 interictal epileptiform discharges (Feys et al., 2021) and it has been shown with high-TC on-scalp 49 SQUID-MEG, that higher proximity of the sensors to the scalp can increase the number of interictal 50 discharges that are observed (Westin et al., 2020). This improvement would be particularly important 51 for children, due to their smaller head size (Boto et al., 2016;Hill et al., 2019) and is pertinent for 52 epilepsy since there is an incidence peak in childhood (Kotsopoulos et al., 2002;Fiest et al., 2017 printed to match their scalp surface, as measured from a previous structural 3T MRI (Meyer et al.,74 2017). The sensors were positioned approximately evenly around the head, with greater coverage 75 over the right, frontal cortex, as shown in Figure 1. The magnetic field perturbations along two axes 76 of the OPMsone radial to the scalp and one tangentialwere recorded, meaning that there were 74 77 recording channels in total. 4 OPMs were later removed from the analysis due to abnormally high 78 noise floors and 2 were removed due to faults in the OPM heating hardware creating erroneous 79 signals. This left 62 recording channels. 80 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. at the centre of the UCL room is approximately 2 nT and the field gradient is approximately 1 nT per 84 metre (Mellor et al., 2021). No additional active shielding was used. The recording was split into 85 three runs of 10 minutes each. In the first and second runs, the patient was reading a book. In the 86 second run, for the first and last minute of the 10-minute recording, we asked the participant to move 87 their head. In the third and final run, the participant sat still and at rest with their eyes open. No overt 88 task was performed during the recording. 89 The OPM data was recorded using a custom acquisition software, built in Labview. Each OPM 90 channel was measured as a voltage (±5 V, 500 Hz antialiasing hardware filter) and sampled at 6 kHz 91 by a National Instruments (NI) analogue-to-digital converter (NI-9205, 16-bit, ± 10 V input range) 92 using QuSpin's adapter (https://quspin.com/products-qzfm/ni-9205-data-acquisition-unit/). This 93 voltage was then multiplied by a calibration factor to calculate the recorded magnetic field. 94 In addition to the OPM data, the patient's movements were recorded with an array of 6 OptiTrack 95 Flex 13 cameras. The cameras were positioned in three corners of the room, with one low and one 96 close to the ceiling in each corner. All cameras were positioned so that the participant's head was 97 within their field of view. The cameras were calibrated according to the manufacturer's instructions, 98 using the OptiTrack CW-500 calibration wand. To track the participant's head, four retroreflective 99 markers were placed onto the OPM scanner-cast. These are shown in blue in Figure 1. A rigid body 100 was created from the four markers in Motive, the acquisition programme associated with the 101 OptiTrack cameras. The rigid body position and orientation was recorded by the Flex 13 motion 102 tracking cameras and then the corresponding positions and orientations of the OPM channels were 103 calculated at each time point. 104 The marker positions were recorded at 120 Hz. A 5 V pulse was used to synchronize the motion 105 tracking with the MEG recordings. Where markers were occluded, usually due to the wires from the 106 OPMs, the missing data was filled in using Motive. Firstly, for any gap where only one marker was 107 missing, a pattern-based method was used to fill in the missing marker data. In this case, the other 108 three markers are used to estimate the position of the missing one. Then for the remaining gaps where 109 multiple markers were occluded, cubic spline interpolation was used to fill in the missing data. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted September 30, 2021. ; https://doi.org/10.1101/2021.09.28.21264047 doi: medRxiv preprint Before IED detection, we performed the following pre-processing steps. Firstly, the OPM data was 112 downsampled to 240 Hz using SPM12. The start and end were then trimmed to match the motion 113 tracking recording and the motion data was upsampled to 240 Hz using linear interpolation in order 114 to match the OPM data. A 10 Hz, 6th order Butterworth low-pass filter was applied to the rigid body 115 position and orientation using FieldTrip to remove high frequency noise caused by marker vibrations. 116 Three 5th order Butterworth notch filters were applied to the MEG data using FieldTrip. One was at 117 50 Hz to remove the mains noise; the other two were at 37 Hz and 83 Hz to remove noise caused by 118 the motion tracking cameras. Environmental noise was then reduced by modelling the background 119 magnetic field as a homogeneous field at each timepoint (Tierney et al., 2021). A 5th order 120 Butterworth high-pass filter was then applied to the MEG data at 1 Hz. 121 IEDs were identified by visual inspection of the pre-processed data. IEDs were identified and 122 manually marked on the OP-MEG traces for each session by consensus between a researcher (SM) 123 and an experienced clinician (UV). The reviewers were blinded to the magnetometer positions. These 124 selections were then used to create 2 s trials. The centre of the epileptiform trials was determined by 125 maximising the cross-correlation with the fourth identified IED trial. For each IED trial, the central 126 timepoint was shifted by ±0.5 s in steps of 1/240 s and the mean (across channels) correlation 127 between the shifted trial and fourth IED trial was calculated. The shift which corresponded to the 128 maximum correlation was chosen to be the centre of the IED trial. The remaining data was split into 129 2 s baseline trials. To ensure that the identified epileptiform activity was not caused by movement 130 artefact, we set a threshold for the amount of variance in the OPM data in the trial which could be 131 explained by the patient's position and orientation, above which a trial was rejected. We set this 132 threshold based on runs 1 and 3, since there was more movement in run 2 which made it 133 unrepresentative of the other two recordings. The threshold was set so that 20% of trials in runs 1 and 134 3 were rejected. This equated to 13.8% of the variance in the trial and led to rejection of 43% of the 135 data, leaving 11 good epileptiform trials and 500 baseline trials. The implications of this conservative 136 rejection criteria are considered further in the discussion. 137 Source localisation was performed using an F contrast on the source power output of an LCMV 138 beamformer from DAiSS (https://github.com/spm/DAiSS). The epileptiform activity was compared 139 with other sections of the data where no abnormal or epileptiform activity was identified. A 140 confidence volume was obtained by bootstrapping which trials were used in the source localisation 141 and then looking at the global peak location for 100 different trial selections. We also performed a 142 dipole fit of the average spike peak (between 0.9 s and 1.1 s in the averaged trial) in FieldTrip. In all 143 cases, we used the Nolte single shell head model (Nolte, 2003). 144

3
Results 145 In total, the data was split into 893, 2 s trials, 382 of which were rejected due to the degree of the data 146 which was explained by movement in the trial. The percentage of variance in the OPM data which is 147 explained by the rotation and position of the patient's head in each trial is shown in Figure 2A. 148 Encouragingly, for most trials, less than 20% of the data can be explained by movement. The  149 distance travelled by the patient in each accepted 2 s trial is shown in Figure 2B. This suggests that 150 these trials were accepted because there is little movement in them. 151 Figure 3 shows the sensor level results of the accepted spike trials. In Figure 3A, a single spike trial 152 is shown on all channels. The selected trial is shaded in blue. Each individual spike trial is shown in 153 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted September 30, 2021. ; https://doi.org/10.1101/2021.09.28.21264047 doi: medRxiv preprint Figure 3B for a single channel. The average of all 11 trials is also shown. The topography of the 154 average is shown in Figure 3C for the radially oriented OPM channels. 155 The source localisation from this epileptiform activity is shown in Figure 4. For both localisation 156 methods, the dysplasia has been circled where visible. The results in Figure 4A are from a 157 beamformer F contrast between epileptiform and non-epileptiform trials. A 2 mm volumetric grid 158 was used for the possible source locations. The peak F-value is shown by the crosshairs. The We then took the beamformer results in two directions. Firstly, we constructed virtual electrodes 169 across a grid to look at the time series of this beamformed activity. This is shown in Figure 5. The 170 time series in Figure 5 corresponds to the time around the first identified IED, as shown at the OPM 171 sensor level in Figure 3A. The spike, as identified in the OP-MEG data, is clearly visible. 172 Constructing virtual electrodes in this way, either inside the head or on the scalp, could be one way to 173 reduce the impact of interference in OP-MEG (Seymour et al., 2021). 174 We also bootstrapped the beamformer, in order to give a confidence volume for the localisation. We

A) B)
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review) preprint
The copyright holder for this this version posted September 30, 2021. ; https://doi.org/10.1101/2021.09.28.21264047 doi: medRxiv preprint from the boundary. The 95% confidence volume was less than 8 mm 3 , the grid spacing of the 180 possible source locations. 181

Discussion 182
We report whole-head OP-MEG recordings of epileptiform activity. This patient had an anatomically 183 identifiable FCD and our source localisation is consistent with epileptiform activity arising from a 184 region bordering the FCD. The patient was seated and head-movement was unconstrained. This is 185 positive for future studies involving recordings over long time periods or with cohorts who are less 186 able to remain still. 187 We have argued here that MRI-lesion positive is a sensible patient population to evaluate OP-MEG, 188 since the lesion location is known. The lesion is a clear anatomical marker and it is highly likely that 189 epileptiform activity arises from either within or around this highly epileptogenic lesion (Jackson,190 Kuzniecky and Berkovic, 2005). Findings from other imaging modalities have shown, however, that 191  . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.  . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. scenarios, the timing of the trials is set by an external stimulus. Here however, we required a 206 researcher or clinician to manually select the interesting events in the sensor level data; sudden head-207

A) B)
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted September 30, 2021. ; https://doi.org/10.1101/2021.09.28.21264047 doi: medRxiv preprint movements or movements of cables had a spike-like appearance which made this selection difficult. 208 We therefore used a conservative criterion to reject epochs showing movement related activity. 209 Consequently, whilst we did not constrain head movement in this study, we have not analyzed data 210 during the movement itself. 211 There are hardware solutions to reduce these issues. In this study, the most pernicious noise was due 212 to the OPM cables moving relative to one another. from the output of a beamformer, as we showed in Figure 5, could foreseeably be a sensible way to 234 present OP-MEG data. This would reduce interference and improve the probability of detecting more 235 subtle epileptiform activity (Van Klink, Hillebrand and Zijlmans, 2016). This is already used in 236 epilepsy in synthetic aperture magnetometry and excess kurtosis mapping (SAM(g2)), in which a 237 beamformer is used to create virtual electrodes in a grid across the cortex. The excess kurtosis in each 238 virtual electrode is calculated and mapped (Robinson et al., 2004). For pediatric epilepsy surgery in 239 particular, this pipeline has been shown to be predictive of seizure outcome (Gofshteyn et al., 2019). 240 As such, there is still much scope for research into how best to process and present OP-MEG data for 241 epilepsy. 242

Conclusion 243
Here we have presented whole-head OP-MEG recordings from a patient with epilepsy. The 244 localisation of the observed interictal activity was within 2.3 cm of the boundary of the FCD lesion, 245 as found in MRI, and appropriately located for the patient's epilepsy for both beamformer and dipole 246 localisations. The patient was seated and their head was unconstrained, which is encouraging for 247 future studies with cohorts who are less able to remain still and for studies over long time periods. 248

6
Conflict of Interest 249 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted September 30, 2021. ; https://doi.org/10.1101/2021.09.28.21264047 doi: medRxiv preprint