Abstract 摘要
Mapping the entire frequency bandwidth of brain electrophysiological signals is of paramount importance for understanding physiological and pathological states. The ability to record simultaneously DC-shifts, infraslow oscillations (<0.1 Hz), typical local field potentials (0.1–80 Hz) and higher frequencies (80–600 Hz) using the same recording site would particularly benefit preclinical epilepsy research and could provide clinical biomarkers for improved seizure onset zone delineation. However, commonly used metal microelectrode technology suffers from instabilities that hamper the high fidelity of DC-coupled recordings, which are needed to access signals of very low frequency. In this study we used flexible graphene depth neural probes (gDNPs), consisting of a linear array of graphene microtransistors, to concurrently record DC-shifts and high-frequency neuronal activity in awake rodents. We show here that gDNPs can reliably record and map with high spatial resolution seizures, pre-ictal DC-shifts and seizure-associated spreading depolarizations together with higher frequencies through the cortical laminae to the hippocampus in a mouse model of chemically induced seizures. Moreover, we demonstrate the functionality of chronically implanted devices over 10 weeks by recording with high fidelity spontaneous spike-wave discharges and associated infraslow oscillations in a rat model of absence epilepsy. Altogether, our work highlights the suitability of this technology for in vivo electrophysiology research, and in particular epilepsy research, by allowing stable and chronic DC-coupled recordings.
绘制整个脑电生理信号频率带宽对于理解生理和病理状态至关重要。使用同一记录位点同时记录直流偏移、极慢振荡(<0.1 Hz)、典型的局部场电位(0.1-80 Hz)和高频(80-600 Hz)的能力,将特别有利于临床前癫痫研究,并可能为改善癫痫发作起始区划分提供临床生物标志物。然而,常用的金属微电极技术在直流耦合记录的高保真度方面存在不稳定性,这阻碍了对非常低频信号的访问。在本研究中,我们使用由线性阵列的石墨烯微晶体管组成的柔性石墨烯深度神经探针(gDNPs),在清醒的啮齿动物中同时记录直流偏移和高频神经元活动。我们在这里表明,gDNPs 可以可靠地记录和以高空间分辨率绘制癫痫发作、发作前直流偏移和与癫痫发作相关的扩散去极化,以及通过皮质层到海马的高频,在化学诱导癫痫发作的小鼠模型中。 此外,我们通过在大鼠失神癫痫模型中记录高保真的自发性尖波放电和相关慢波振荡,展示了植入式设备在 10 周以上的慢性功能。总的来说,我们的工作突出了这项技术在活体电生理学研究,特别是癫痫研究中的适用性,因为它允许稳定的慢性直流耦合记录。
Similar content being viewed by others
其他人在查看的类似内容
Main 主要
Pathological electrophysiological signals resulting from epilepsy span a wide frequency range, from fast activity, at hundreds of hertz or higher, to less explored signals below 0.1 Hz, including DC-shifts (sustained over seconds to minutes) and infraslow oscillations1,2,3. In contrast to the widely used AC-coupled recordings, DC-coupled recordings allow more comprehensive assessment of the pathophysiological processes associated with either an increase in neuronal activity during seizures or a loss of neuronal activity during a spreading depolarization (SD)4. Clinically, electrophysiological biomarkers have been explored to guide surgical resection of the seizure onset zone5,6; however, improvement and refinements are still needed7,8. Indeed, intracranial clinical studies using large metallic electrodes suggest that the efficacy of current electrophysiological biomarkers could be improved by the inclusion of DC-shifts9. However, broad adoption of DC-shifts as a novel biomarker is limited by the poor capabilities of microelectrode technology to reliably record this type of electrophysiological signal in vivo with high spatiotemporal fidelity. This is due to polarization-induced drift and signal attenuation causing distortion of the measured signal by metallic electrodes10,11,12. Currently, the gold standard for high-fidelity DC-coupled recordings is the glass micropipette; however, its use is limited to few-point measurements, hampering the spatial resolution required for seizure localization, and is not suitable for clinical applications. An alternative to commonly used passive electrodes are field-effect transistors (FETs), which are active transducers offering important advantages in electrophysiology13. Among the few FET technologies that have been validated for in vivo electrophysiology, graphene-based FETs are particularly attractive because of the combination of material properties, including chemical and electrochemical inertness, high electrical mobility, biocompatibility, as well as facile integration into flexible and ultrathin substrates14. Recent reports demonstrate the potential of graphene solution-gated field-effect transistors (gSGFETs) for neural interfacing15,16. A proof-of-concept of high-fidelity DC-coupled recordings enabled by gSGFETs was demonstrated by mapping chemically induced cortical SDs in anaesthetized rats17. To explore its usefulness for epilepsy research, we have developed implantable, flexible graphene depth neural probes (gDNPs) capable of recording localized full-bandwidth neuronal activity through cortical columns and subcortical structures. Here, we present a wafer-scale microtechnology process to fabricate ultrathin gDNPs consisting of a linear array of graphene microtransistors imbedded in a polymeric flexible substrate. To penetrate through the mouse cortex and reach the hippocampus without buckling, we adapted an insertion protocol that uses silk fibroin (SF)18,19 to temporarily stiffen flexible gDNPs. The ability to detect electrophysiological biomarkers of epileptiform activity in preclinical rodent models of induced seizures and chronic epilepsy has been validated experimentally20. These included seizure events such as spike-and-wave discharges (SWDs)21 and high-frequency oscillations (HFOs)22,23. Additionally, we highlight the suitability of graphene transistor technology to record concurrently fast and slow brain electrophysiological signals by benchmarking against a solution-filled glass micropipette recording, demonstrating the high fidelity of the DC-coupled recordings obtained by the gDNPs. Here, we demonstrate the capability of gDNPs to record both active and passive ictal baseline shifts24,25,26 as well as SD associated with seizure activity. Altogether, our results illustrate the potential of graphene microtransistor technology to reveal mechanistic insights into the involvement of infraslow oscillations and DC-shifts associated with seizures in vivo in awake brain.
癫痫引起的病理生理电信号跨越广泛的频率范围,从数百赫兹或更高的快速活动,到尚未充分探索的低于 0.1 赫兹的信号,包括直流偏移(持续数秒至数分钟)和极慢振荡1,2,3。与广泛使用的交流耦合记录相比,直流耦合记录允许更全面地评估与癫痫发作期间神经元活动增加或去极化扩散(SD)期间神经元活动丧失相关的病理生理过程4。在临床上,已经探索了电生理生物标志物来指导癫痫发作起始区的手术切除5,6;然而,仍需改进和细化7,8。实际上,使用大型金属电极的颅内临床研究表明,通过包括直流偏移,可以改善当前电生理生物标志物的有效性9。 然而,微电极技术在体内以高时空保真度可靠记录这种电生理信号的能力有限,这限制了 DC 偏移作为一种新型生物标志物的广泛应用。这是由于金属电极引起的极化诱导漂移和信号衰减,导致测量的信号失真10,11,12。目前,高保真度直流耦合记录的金标准是玻璃微管;然而,其使用仅限于少量点测量,阻碍了癫痫定位所需的空間分辨率,并且不适合临床应用。一种替代常用被动电极的方法是场效应晶体管(FETs),它们是主动传感器,在电生理学中提供了重要的优势13。 在为数不多的已被验证用于体内电生理学的 FET 技术中,基于石墨烯的 FETs 因其材料特性的组合而特别吸引人,这些特性包括化学和电化学惰性、高电迁移率、生物相容性,以及易于集成到灵活和超薄基板中14。最近的报告表明,石墨烯溶液栅场效应晶体管(gSGFETs)在神经接口方面具有潜力15,16。通过映射麻醉大鼠中化学诱导的皮质 SDs,展示了 gSGFETs 实现高保真度直流耦合记录的原理证明17。为了探索其在癫痫研究中的实用性,我们开发了可植入的、灵活的石墨烯深度神经探针(gDNPs),能够通过皮质柱和皮层下结构记录局部全带宽神经元活动。在这里,我们提出了一种晶圆级微技术工艺,用于制造由嵌入在聚合物柔性基板中的线性阵列石墨烯微晶体管组成的超薄 gDNPs。 为了穿透小鼠皮层并到达海马体而不弯曲,我们采用了一种插入协议,该协议使用丝素蛋白(SF)18,19来暂时使柔性的 gDNPs 变硬。在诱导性癫痫和慢性癫痫的啮齿动物模型中检测癫痫样活动的电生理生物标志物的能力已经通过实验得到验证20。这些包括癫痫事件,如尖波和波放电(SWDs)21和高频振荡(HFOs)22,23。此外,我们通过将 gDNPs 与溶液填充的玻璃微管记录进行基准测试,突出了石墨烯晶体管技术同时记录快速和慢速脑电生理信号的适用性,证明了 gDNPs 获得的直流耦合记录的高保真度。在这里,我们展示了 gDNPs 记录活动性和被动性癫痫发作基线偏移的能力24,25,26以及与癫痫活动相关的 SD。 总的来说,我们的结果表明,石墨烯微晶体管技术在揭示清醒大脑中与癫痫发作相关的亚慢振荡和直流偏移的机制方面具有潜力。
Microfabrication, characterization and stiffening of the gDNPs
微加工、gDNPs 的表征和强化
A graphene-based SGFET is a three-terminal device in which single-layer graphene prepared by chemical vapour deposition (CVD) is used as the channel material in contact with the drain and source terminals. Graphene is the sensing part of the device directly exposed to the neural tissue. The current in the graphene channel can be modulated or pinned by a third terminal given by a reference electrode (gate), which is also in contact with the conducting neural tissue (Fig. 1a). Thus, variations in the electrical potential in the tissue can be transduced into variations in the channel current; this transduction mechanism has been shown to offer a very versatile sensing platform for electrophysiology27,28. The flexible gDNP prepared in this study consists of a linear array of 14 recording transistors, each with an active area of 60 × 60 μm2, and a pitch of 100 μm. The probe’s tip design consists of a polyimide shank of 200 μm width and 1.6 mm length (Fig. 1b). A gDNP design with smaller SGFETs (20 × 20 μm2) and narrower shank (150 μm width) has also been fabricated (Supplementary Fig. 1). The gDNPs were fabricated on a 10-μm-thick flexible polyimide (PI) substrate using a wafer-scale microfabrication process previously reported16 (Methods). A two-level metallization strategy, with the metal levels interconnected by via holes (Fig. 1b), reduces track resistance and improves sensor performance. To characterize gDNPs in saline solution we simultaneously measure the drain–source current (IDS) versus the applied gate–source voltage (VGS) for all the transistors on the shank with a fixed drain–source voltage (VDS) using customized electronics (Methods). gSGFETs exhibit very high transconductance (gm) values due to direct contact with the tissue, possible as a result of its biocompatibility and electrochemical inertness, and to the very high charge carrier mobility of graphene29. Figure 1c shows the transfer curves as well as the normalized transconductance (|gm| = gm/VDS) for all 14 gSGFETs of an exemplary gDNP device. The small dispersion of the charge neutrality point (CNP), defined as the value of VGS where the IDS reaches a minimum, attests to the homogeneity of the gDNP. Furthermore, gm shows a very stable response over a wide applied gate frequency range (up to 5 kHz), both in the hole regime, VGS < VCNP, and in the electron regime, VGS > VCNP (Fig. 1d). Such a constant frequency response is important for a proper calibration of the recorded signals30. The detection limit of the sensors was evaluated by means of the effective gate noise (that is, root-mean-square gate voltage noise, VRMS) integrated between 1 Hz and 2 kHz, with averaged values of 25–30 µV for all fabricated gDNPs (Supplementary Fig. 2).
基于石墨烯的 SGFET 是一种三端器件,其中使用化学气相沉积(CVD)制备的单层石墨烯作为与漏极和源极端子接触的通道材料。石墨烯是直接暴露于神经组织的器件的传感部分。石墨烯通道中的电流可以通过由参考电极(栅极)提供的第三个端子进行调制或固定,该参考电极也接触导电的神经组织(图1a)。因此,组织中的电势变化可以被转换为通道电流的变化;这种转换机制已被证明为电生理学提供了一个非常通用的传感平台27,28。本研究中制备的柔性 gDNP 由一个线性阵列的 14 个记录晶体管组成,每个晶体管的活性面积为 60×60μm2,间距为 100μm。探针的尖端设计由 200μm 宽和 1.6mm 长的聚酰亚胺支架组成(图1b)。还制造了具有更小 SGFET(20×20μm2)和更窄支架(150μm 宽)的 gDNP 设计(补充图。 1) gDNPs 是在 10 微米厚的柔性聚酰亚胺(PI)基板上使用之前报道的晶圆级微加工工艺制造的16 (方法)。一种两级金属化策略,金属层通过通孔(图1b)互连,降低了轨道电阻并提高了传感器性能。为了表征盐水溶液中的 gDNPs,我们同时测量了所有在支架上的晶体管的漏源电流(IDS)相对于施加的栅源电压(VGS),使用固定的漏源电压(VDS)和定制的电子设备(方法)。gSGFETs 表现出非常高的跨导(gm)值,这是由于与组织的直接接触,可能是由于其生物相容性和电化学惰性,以及石墨烯非常高的电荷载体迁移率29。 图 1c 展示了所有 14 个 gSGFET 的传输曲线以及归一化跨导(|gm| = gm/VDS)的示例 gDNP 器件。电荷中性点(CNP)的小分散性,定义为 VGS 的值,其中 IDS 达到最小值,证明了 gDNP 的一致性。此外,gm 在宽范围的施加栅极频率(高达 5 kHz)下表现出非常稳定的响应,无论是在空穴区,VGS < VCNP,还是在电子区,VGS > VCNP(图 1d)。这种恒定的频率响应对于正确校准记录的信号30非常重要。 传感器的检测限通过有效门噪声(即均方根门电压噪声,VRMS)在 1 Hz 和 2 kHz 之间的积分来评估,所有制造的 gDNPs 的平均值为 25–30 µV(补充图 2)。
gDNPs are highly flexible, compared with traditional rigid depth electrodes, and although flexibility is advantageous once inserted into tissue, this provides a challenge during insertion. To insert these probes we temporarily stiffened the gDNP using SF18,19. The stiffening technique (Methods) consisted of a moulding process in which the gDNP is back-coated with SF, allowing the preparation of a rigid shank with a defined shape and thickness. We tuned the thickness of the SF by controlling the mould’s trench depth, achieving two typical thicknesses of 80 ± 10 and 150 ± 12 μm (Fig. 1e). The mechanical properties of the SF-coated gDNPs were assessed using a buckling test, in which perpendicularly positioned probes were driven against a flat, hard surface (Fig. 1f). An initial linear increase in force was observed for both coating thicknesses tested, the probes remaining straight before buckling. Continued application of force resulted in buckling and bending, characterized by a peak in the force–displacement curve. The obtained peak forces, 101 ± 21 mN for the 150-µm-thick SF and 29 ± 13 mN for the 80-µm-thick SF, are in good agreement with the previously reported values of peak forces for similar SF-coated neural probes19,31. We electrically characterized the gDNPs before and after the SF stiffening process, as well as before and after insertion into and removal from an agarose gel brain model. Figure 1g shows the distribution of the normalized gm as well as VRMS for all 14 transistors in a gDNP, confirming that neither the stiffening process nor insertion in an agarose brain model impairs gDNP performance in terms of transconductance or noise (Supplementary Fig. 3).
gDNPs 具有高度灵活性,与传统刚性深度电极相比,尽管在插入组织后灵活性具有优势,但在插入过程中这带来了挑战。为了插入这些探针,我们使用 SF18,19暂时使 gDNP 变硬。硬化技术(方法)包括一个成型过程,其中 gDNP 背面涂覆 SF,允许制备具有定义形状和厚度的刚性柄。我们通过控制模具的沟槽深度来调整 SF 的厚度,实现了两种典型的厚度,分别为 80±10 和 150±12μm(图1e)。使用屈曲测试评估了 SF 涂覆的 gDNPs 的机械性能,其中垂直放置的探针被驱动到平坦、硬的表面上(图1f)。观察到两种测试的涂层厚度在初始阶段力线性增加,探针在屈曲前保持直线。继续施加力导致屈曲和弯曲,在力-位移曲线上表现为峰值。 获得的峰值力,150-µm 厚的 SF 为 101±21 mN,80-µm 厚的 SF 为 29±13 mN,与之前报道的类似 SF 涂层神经探针的峰值力值19,31非常吻合。我们在 SF 硬化过程前后,以及插入和从琼脂糖凝胶脑模型中取出前后,对 gDNPs 进行了电学表征。图1g显示了所有 14 个 gDNP 中归一化的gm以及VRMS的分布,证实了硬化过程或插入琼脂糖脑模型中都不会损害 gDNP 在跨导或噪声方面的性能(补充图3)。
Figure 1h depicts video frames of a SF-coated gDNP after insertion (insertion speed 400 μm s–1) in an agarose brain model; the images help to visualize the bending-free insertion of the gDNP shank and the fast water absorption (<10 s) by SF. Moreover, it can be seen that the SF was completely delaminated from the polymeric shank after removal from the agar gel (180 s).
图 1h 展示了在琼脂脑模型中插入(插入速度 400 μm s–1)后的 SF 涂层 gDNP 的视频帧;这些图像有助于可视化 gDNP 腿部的无弯曲插入和 SF 的快速吸水(<10 s)。此外,还可以看到,在从琼脂凝胶中取出后(180 s),SF 完全从聚合物腿部剥离。
These observations confirm the suitability of the SF back-coating stiffening procedure to enable the flexible gDNPs to reach the target position in a brain-like tissue without bending; the fast dissolution of the shuttle once inserted enables the probe to quickly regain flexibility. This insertion strategy is only suitable for single-time use; for further insertions, the flexible gDNP needs to be recoated. We have verified that the gDNP technology can withstand multiple coating and insertion procedures (Supplementary Fig. 4).
这些观察结果证实了 SF 背涂加固程序适用于使柔性 gDNPs 在类似大脑组织中到达目标位置而不弯曲;一旦插入,穿梭的快速溶解使探针能够迅速恢复灵活性。这种插入策略仅适用于单次使用;对于进一步的插入,柔性 gDNP 需要重新涂覆。我们已经验证了 gDNP 技术可以承受多次涂覆和插入程序(补充图4)。
Awake in vivo full-bandwidth recording with gDNPs
使用 gDNPs 的清醒状态活体全带宽记录
We assessed full-bandwidth recording capability by implanting a gDNP into awake, head-fixed mice. The electrophysiological signal measured by the graphene transistors was acquired with a customized electronic system that enables simultaneous recording in two frequency bands with different gains preventing amplifier saturation (Fig. 2a and Methods). gDNPs were implanted in the right hemisphere visual cortex (V1) and lowered until the tip reached hippocampal tissue. Baseline activity was recorded for 10–20 min. To induce network discharges and synchronicity of neuronal bursting, 200 nl of 4-aminopyridine (4-AP, 50 mM), a selective blocker of Kv1 potassium channels32,33, was focally injected into the cortex adjacent to the gDNP (Fig. 2b).
我们通过在清醒、头部固定的老鼠中植入 gDNP 来评估全带宽记录能力。由石墨烯晶体管测量的电生理信号是通过一个定制的电子系统获得的,该系统允许在两个具有不同增益的频率带中同时记录,以防止放大器饱和(图2a和方法)。gDNPs 被植入右侧半球视觉皮层(V1)并降低,直到尖端到达海马组织。记录了 10-20 分钟的基线活动。为了诱导网络放电和神经元爆发的同步性,将 200 nl 的 4-氨基吡啶(4-AP,50 mM),一种选择性阻断 Kv1 钾通道的阻断剂32,33,聚焦注射到 gDNP 旁边的皮层中(图2b)。
Full-bandwidth recordings
Figure 2c displays an uninterrupted electrophysiological recording over 2 h (only the uppermost and the lowest channels of the implanted gDNP are displayed); the complete dataset is shown in Supplementary Fig. 5. The ability of the graphene transistors to provide long and stable full-bandwidth recordings without the need for electronic off-set readjustments contrasts the limitation of DC-coupled passive electrodes34. The recording shows the full-bandwidth signal (>0.001 Hz) and the signal high-pass (HP)-filtered above 0.5 Hz (which is the expected signal recorded by AC-coupled electrodes)35. The recording is broken down into different experimental conditions during the recording: baseline, needle-induced pin-prick SD36, interictal activity5 induced by chemoconvulsant drugs and seizures. During baseline recording (Fig. 2c), the lower channels exhibit theta activity, correlated with animal movement, indicating that the gDNP reached the hippocampus, confirmed post hoc by histological analysis of fixed brain sections (Supplementary Fig. 6). After the injection of 4-AP, epileptiform spiking evolved and five seizures (over 60 min) were elicited in this example, one of which was followed by a SD. A second cortical injection of 4-AP induced two additional seizures, both followed by SDs that were detected first in the hippocampus. In five different mice injected with 4-AP, an average number of 7 ± 3 seizures were recorded in 60 min after drug injection. In this chemoconvulsant model, SDs could be observed initially either in superficial cortical layers or in the hippocampus (Fig. 2c and Supplementary Fig. 7).
The fidelity of the DC-coupled recordings was validated by simultaneous recordings using a solution-filled glass micropipette. Figure 2d shows the full-bandwidth recording obtained with the gDNP and micropipette after the injection of 4-AP. Both recordings reveal DC-shifts preceding seizure, and SD occurring after seizure. The DC-shifts measured by the gDNP have a similar shape, magnitude and temporal duration to the signal recorded by the glass micropipette. A cross-correlation analysis (signal filtered <5 Hz) of the signals recorded by the glass micropipette and the 14 gDNP transistors demonstrated a very high correlation (above 90%) for the channels located at the same cortical depth as the micropipette. The same high correlation was observed during hour-long recordings and for small, sudden DC-shifts (Supplementary Fig. 23).
Assessment of the detection limits of gDNP
Post-mortem recordings were acquired to characterize the electrical noise level of the gDNP in the activity-free brain state and, consequently, to quantify the detection limit of the gDNP. Figure 2e shows the averaged power spectral density (PSD) calculated using the recordings of all channels in a gDNP, obtained from different brain states (baseline, after injection of 4-AP and post-mortem). Compared with the baseline PSD, the large amplitude of the PSD at low frequencies (<1 Hz) after 4-AP injection provides an indication of the interictal and DC-shifts in the brain. We calculated the signal-to-noise ratio (SNR) in three different bands, that is, 1–70, 80–200 and 200–4,000 Hz (Fig. 2f). This analysis shows that the gDNPs are able to record electrophysiological signals in typical bandwidths with a SNR ratio higher than 1 dB.
Electrophysiological recording of characteristic epilepsy biomarkers
Epilepsy alterations of physiological brain activity include interictal spikes, fast activity at hundreds of hertz (such as HFOs) and pre-seizure active DC-shifts. We used two different pro-convulsive drugs (4-AP or picrotoxin (PTX)) to induce and evaluate epileptic activity. Figure 3a shows the interictal activity and associated HFOs (>80 Hz)22,23 recorded by three of the transistors of a gDNP, each located at a different depth in the mouse brain. Filtering between 80 and 600 Hz revealed layer-specific bursting of HFOs and sharp wave ripples during interictal spikes with characteristic oscillations of 200–300 and 400–600 Hz in the cortical and hippocampal channels, respectively37 (Supplementary Fig. 8). Entrained interictal epileptiform activity was found in all channels before each seizure. Figure 3b illustrates characteristic examples of sharp wave ripples and HFOs induced by 4-AP and PTX recorded by the lowest channel of the gDNP (hippocampal CA1 region). The HFO and ripple traces shown in Fig. 3b exhibit high-frequency tones up to 600 Hz. The filtered traces (>200 Hz) are compared with the original traces (full bandwidth) for verification of the ripples. The advantage of the gDNPs for monitoring concurrent DC-shifts and high-frequency activity is illustrated in Fig. 3c, which shows an epileptic spike-triggered SD arising from the hippocampus. The layer-dependent silencing of the neural activity by the hippocampal SD is represented in Fig. 3d in terms of the percentage of activity variation. The layer-dependent amplitude of the SD and the following hyperpolarization, also shown in Fig. 3d (more details of this analysis are provided in the Methods and Supplementary Fig. 9), reveal that the silencing of the neural activity in the hippocampus is correlated with the amplitude and subsequent hyperpolarization wave of the SD. The silencing of neural activity in the hippocampus by the SD is visualized with more clarity in Fig. 3e, where the spectrograms for the uppermost and lowest channels are compared.
DC-shifts are typically associated with seizures24 and can be differentiated into active shifts that precede the seizure, and passive shifts that follow the intense neuronal firing at seizure onset38. For these chemically induced seizures events, we observed that the DC-shift preceded the seizure onset (Ts) by 2–3 s (n = 4 seizures; Fig. 3f and Supplementary Fig. 10). Moreover, we observed a layer-dependent amplitude of the DC-shift, with negative amplitudes in the cortical layers (close to the 4-AP injection site) and positive amplitudes recorded in the channels close to the hippocampus (Fig. 3f, bottom left). The spatiotemporal dynamics of the DC-coupled extracellular potential for the pre-ictal to ictal transition is depicted in Fig. 3f (bottom right).
The inversion of DC-shifts can be used to identify current sources and sinks through the cortical laminae. Current source density (CSD) analysis is a technique to identify source activation in a variety of focal neurological disorders, including epilepsy39,40. CSD analysis (Methods) of the data in Fig. 3c revealed a large net ionic sink in the hippocampal extracellular space after the seizure, followed by a large source at the beginning of the hyperpolarization wave (Fig. 3g). Enlarging the seizure onset region, four sink and source regions can be identified from the laminae profile. CSD analysis computed with a HP-filtered signal (0.5 Hz < f < 70 Hz) failed to report the ionic sinks preceding and during the seizure in the bottom layers (Fig. 3g, AC-coupled), illustrating the importance of using full-bandwidth recordings for CSD analysis to avoid misinterpretation of the CSD sinks and sources (Supplementary Fig. 11 illustrates additional examples of CSD analyses).
Chronic functional validation and biocompatibility assessment
We implanted gDNPs in the right-hemisphere somatosensory cortex of WAG/Rij rats (n = 4), a rodent model of absence epilepsy20, and obtained chronic full-bandwidth recordings over a 10-week period (Fig. 4a). The WAG/Rij rats exhibited frequent spontaneous SWDs, a characteristic thalamocortical oscillation of 8–10 Hz (ref. 20; see Methods). Implanted animals were connected 1–2 times per week for tethered recordings (using a commutator to enable free movement of the rats; Fig. 4a and Supplementary Fig. 12).
Transistor curves were measured in each recording session to assess device stability, changes at the device–tissue interface and importantly, to permit selection of an optimal VGS to maximize the SNR, a feature possible with active sensor devices28. Figure 4b shows the averaged transfer curves of a gDNP measured over 10 weeks (Supplementary Fig. 13). The stability of the transistors’ performance is illustrated in Fig. 4c, which depicts the position of the CNP and the maximum value of gm over the implantation period. The averaged gm value remains approximately constant over the whole study. Because gm is directly related to the tissue–graphene interfacial capacitance and to the carrier mobility in graphene, the stability of gm strongly suggests little or no variation of these two parameters. The CNP (Fig. 4c) shows a notable shift (200 mV) during the few first days after implantation, and then remains stable for the rest of the experiment. We tentatively attribute the initial shift to the adsorption of negatively charged species, which reduce the intrinsic p-type doping of the graphene transistors29.
The long-term functionality of the gDNP was assessed by evaluating the quality of the recorded signals over the implantation period using two parameters, normalized SWD power amplitude and the transistor noise (VRMS). For further details, see Methods and Supplementary Figs. 13 and 14. The variation of these two parameters over time is shown in Fig. 4d, demonstrating the ability of the implanted devices to monitor seizure activity with high fidelity during the whole implantation period. The slight increase in the normalized SWD activity could result from a strengthened coupling between neural tissue and the gDNP or from an increase in the seizure power and duration as the animal ages41.
Correlation between SWDs and infraslow oscillations
The WAG/Rij rat chronic model of absence epilepsy offers the possibility to investigate correlations between infraslow oscillations (ISOs) and SWD events42. Owing to the full-bandwidth capability of the gDNP, we were able to observe these correlations as well as a dipole across the neocortex upper and lower layers. Figure 4e shows the uppermost and lowest channels having opposite phases during each SWD (Supplementary Fig. 15). The information contained in the DC-shift is lost after application of a HP filter (>0.5 Hz), as typically used with conventional AC-coupled microelectrodes (Fig. 4e, right). This observation persists over the implantation period, as shown in Fig. 4f. To further evaluate the relation of ISOs and SWDs, we calculated the phase of the signals (limited to the 0.01–0.1 Hz band) and the power of the neural activity associated with the SWDs (5–9 Hz). Figure 4g displays the full band and the ISO component (0.001–0.1 Hz) of the signal recorded from one of the upper channels, together with the RMS of the SWD signal (5–9 Hz) and the ISO phase. Prior to the SWD events, whose onset is marked by the abrupt increase of the SWD RMS, the ISO phase is consistently observed to change (Supplementary Figs. 15 and 16). The correlation between the ISO phase and the SWD can be better illustrated by the joint probability distribution plot (Fig. 4g, right) obtained from a long recording (1,600 s, >35 SWDs; see Methods). Figure 4h depicts the correlation between the ISO phase and the RMS of the SWDs for each channel, revealing an inversion in the lower layer of the neocortex (Supplementary Fig. 16). To further confirm the capability to observe these ISOs preceding the seizure events, we calculated the depth profile of the DC-coupled signal (averaged for n = 9 SWD events) in a timeframe around the start of the SWDs. Preceding the SWD (with its characteristic positive DC-shift) by a couple of seconds, a small negative DC-shift was observed at superficial to medial layers (Supplementary Fig. 17).
This correlation between ISO and SWD was also observed in the other implanted WAG/Rij rats (n = 4; Supplementary Fig. 18).
Chronic biocompatibility of implanted gDNP
We also conducted an extensive chronic biocompatibility study to assess any potential neuroinflammation caused by the invasive nature of the penetrating neural probes, the presence of CVD graphene or by the release of SF following implantation. Neural probes with or without graphene at the recording sites were implanted in adult, male Sprague Dawley rats (n = 20). Histological and immunohistochemical studies were conducted at 2, 6 and 12 weeks post-implantation and compared with the contralateral hemisphere without device implantation (Supplementary Figs. 19 and 20). An additional small cohort of animals (n = 3) were implanted with commercially available rigid depth neural probes with similar dimensions (NeuroNexus, Supplementary Fig. 19). Figure 4i shows a fluorescence image of glial fibrillary acidic protein (GFAP) immunostaining (positive marker for astrocyte cells) in the area of insertion 12 weeks post-implantation (brain sections at 800 μm from the pia). There was no significant difference in the number and morphology phenotype of astrocyte cells (typically associated with inflammation) in the area directly surrounding the implanted probe site, or on comparison between probes with or without graphene (Fig. 4i, right). No significant difference was observed at 2, 6 and 12 weeks post-implantation, and values are in the same range as in the contralateral hemisphere where no probe was inserted (Supplementary Fig. 19). The activation state of microglia, assessed by morphological analysis of cells stained for ionized calcium binding adaptor molecule 1, showed no significant increase in the abundance of activated microglia present in the area surrounding the implantation site (Fig. 4j). Additional immunohistochemical analysis showed no sign of an inflammatory response over the 12-week period for either device used, comparable to commercially available rigid probes (Supplementary Fig. 19). Recently published biocompatibility data from some of our laboratories using flexible gFET devices designed with large graphene surface areas implanted epicortically also showed no significant adverse tissue response chronically43. Altogether, the chronic biocompatibility study reported here indicates that gDNPs are suitable for deep brain structure chronic implantations, inducing no significant tissue damage nor neuroinflammatory responses.
Outlook
We have demonstrated in this study the capability of gDNPs based on linear arrays of gSGFETs to record electrophysiological brain activity in a wide temporal range, with fidelity similar to solution-filled glass micropipettes but with superior spatial sampling. The gDNP is a novel technology that allows concurrent recordings of DC-shifts, ISOs and HFOs in preclinical epileptic models through the cortex laminae to upper hippocampal layers.
Together with the validation of the chronic functionality of implanted gDNPs and their biocompatibility, our work underlines the very distinct advantages of this technology for in vivo epilepsy research. In particular, the gDNPs have allowed us to reveal in rodent models of drug-induced seizures and chronic epilepsy20 associations between DC-shifts (including ictal baseline shifts and SDs24,25,26,44,45) and higher-frequency activity (including SWDs21 and HFOs22,23).
The reliability, stability and high fidelity of the recordings observed in this work suggest that this technology could be eventually used to record clinically relevant electrophysiological biomarkers. In this respect, we envision that clinical development of this technology could be directed towards neural probes capable of simultaneous high-quality wide-band (DC to HFO) recordings from multiple brain regions during pre-surgical monitoring. Moreover, this technology could be applied combining depth and subdural gSGFET arrays capable of stable DC-coupled recordings with high spatial resolution, which will offer advantages over current technologies for epilepsy and neurocritical care11,46. In particular, it is expected that this technology can help improve the identification of seizure onset zones and ultimately surgical outcome5. Currently, long-term seizure freedom after surgery is only around ~50% (ref. 47); thus, improvements in a technology that more precisely defines surgical margins is warranted. The few studies that have investigated DC-shifts report that these signals co-localize with seizure onset zones and are more spatially restricted than conventional electroencephalogram (EEG) recordings48. Therefore, routine inclusion of accurate DC-coupled recordings in surgical monitoring, which can be enabled by the technology described in this work, could result in less extensive, yet more effective surgical resections9. Beyond its use for epilepsy research, we envision that this technology can be applied to advance our understanding of other neurological diseases associated with SD, including traumatic brain injury, stroke and migraine49.
Methods
Microfabrication of flexible gDNPs
Flexible neural probes were fabricated using standard microelectronic fabrication technology on a rigid 10-cm sacrificial Si/SiO2 wafer. A 10-µm-thick PI (PI-2611, HD MicroSystems) layer was spin-coated and cured at 350 °C in a N2 atmosphere. To reduce the shank width of the depth neural probes, we used a two-metal-level strategy in which the metal tracks, separated by PI, are interconnected by via holes. After evaporation and definition (by lift off) of a first metal layer of Ti/Au (10/100 nm), a 3-µm-thick PI layer was spin-coated and cured. A protective mask of Al (200 nm) was used to etch the second PI layer using oxygen plasma and form the via holes. On top of the via holes a second metal layer of Ti/Au (10/100 nm) was applied to interconnect the two metal layers. CVD graphene was transferred onto the patterned wafer as described in the Supplementary Information Methods. The graphene active areas were defined by oxygen-based reactive ion etching. A sandwich-like contact strategy was used to improve the contact at the drain and source terminals; the top metal structure consisted of Ni/Au (20/200 nm). For passivation, a 2-μm-thick chemically resistant polymer was deposited (SU8-2005 MicroChem) with open windows in the channel region. Finally, the gDNP structure was defined in a deep reactive ion etching process using a thick AZ9260 positive photoresist (Clariant) as etching mask. The PI probes were then released from the SiO2 wafer and placed in a zero insertion force connector to interface our custom electronics. The transconductance of the flexible gSGEFTs could be improved by ultraviolet ozone cleaning of the contact region at the graphene–metal interface, as recently demonstrated for flexible graphene neural sensors50. For CVD graphene growth, transfer and its characterization we refer to the Supplementary Information.
Characterization of gDNPs in saline
The gSGFETs on the neural probes were characterized in PBS solution (150 mM). Drain–source currents (IDS) were measured varying the gate–source voltage (VGS) versus a Ag/AgCl reference electrode, which was set to ground. Steady state was ensured by acquiring the current only after its time derivative was below a threshold (10−7 A s−1). The detection limit of the gSGFET was assessed by measuring the PSD of the DC current at each polarization point VGS. Integrating the PSD over the frequencies of interest (1 Hz–2 kHz) and using the transconductance allowed us to calculate VRMS. The noise measurement was performed in a Faraday cage, with DC-batteries powering the amplifiers, to reduce any 50 Hz coupling or pick-up noise. Additionally, the frequency response of the device’s transconductance was measured by applying a sum of sinusoidal signals to the electrolyte solution through the reference electrode and by measuring the modulation of the drain current. The acquired signals were split into two bands: low frequencies (~0–10 Hz), at which the drain–source current was simultaneously acquired for all transistors in an array, and higher frequencies (10 Hz–30 kHz), at which each transistor was recorded individually.
Back-coating of gDNP with silk fibroin
Compared with other natural biopolymers, SF offers excellent mechanical properties, extremely good biocompatibility, biodegradability and the versatility of structural readjustments51,52,53. Further, the byproducts of SF degradation by enzymes (for example, proteases) have low antigenicity and non-inflammatory characteristics54,55. Temporary stiffening of the gDNPs with SF was achieved using a microstructured polydimethylsiloxane (PDMS) mould with the shape of the neural probes. To fabricate the moulds, PDMS was cast on a standard 4-inch silicon wafer with prepatterned 100- and 200-μm-thick SU8 (SU8-2050) epoxy resin. The back-coating procedure was as follows: first, the probe was placed in the mould trench filled previously with water, with the transistor side facing down. Through surface tension effects, the probe self-aligned in the mould. After evaporation of the water, SF (Silk, Fibroin Solution 50 mg ml–1, Sigma Aldrich) was applied through a syringe to the mould’s trench. We double-coated the shank in drying intervals of 20 min and then slowly increased the temperature on a hotplate to 80 °C, leaving the SF to cure for 90 min. To have a delayed dissolution time compared with the SF curing time at room temperature, we increased the duration of the water annealing step. After curing, the coated probe can be easily removed from the PDMS mould (Supplementary Fig. 3). In all the in vivo experiments presented here, we implanted the flexible gDNPs with a 150-µm-thick SF back-coating.
Assessment of mechanical properties of the stiffened gDNPs
Standard compression tests against a hard Si substrate were performed to assess the mechanical properties of our SF-coated probes. Buckling experiments were carried out in a UMIS nanoindenter (Fischer-Cripps Laboratories). A custom clamp was fabricated to fix the probes at the end of the indenter shaft, which, in turn, was connected to an actuator and load cell. Buckling tests were carried out at a loading rate of 8.8 mN s–1. Once the indenter detected noticeable buckling the test was automatically stopped. The maximum applied load that the indenter can apply is 500 mN. Applied force versus displacement was measured until the probe started to buckle and eventually broke down. We additionally measured the Young’s modulus of the SF cured at 80 °C by nanoindentation tests. SF was dropcast on a 2 × 2 cm2 Si chip and cured. Indentation was measured using an NHT2 Nanoindentation Tester (Anton-Paar) equipped with a Berkovich pyramidal-shaped diamond tip. A maximum applied load of 5 mN was applied with a loading segment of 30 s followed by a load holding segment of 10 s and an unloading segment of 30 s. The hardness and reduced Young’s modulus are reported as average values of at least 20 indentations, performed on top of each sample (in the central region). Young’s modulus values in the range of 10 GPa were measured for SF cured at 80 °C (Supplementary Fig. 22).
Electronics for in vivo recordings with gDNPs
The experimental set-up used to perform the in vivo recordings provided VS and VD bias control and current-to-voltage conversion for up to 16 channels (g.RAPHENE, g.tec medical engineering). The instrumentation split the recorded signals into two bands with different gains: low-pass filtered (LPF; <0.16 Hz, 104 gain) and band-pass filtered (BPF; 0.16 Hz < f < 160 kHz, 106 gain). Two custom Simulink models were used (1) to generate the transfer curve of the microtransistors once inserted and at the end of the experiment, and (2) to set the VS and VD bias and acquire the recorded signals. Signals were sampled at 9.6 and 19 kHz, depending on the type of experiment (Supplementary Table 1).
Ethical approval and animal handling for acute and chronic experiments
Animal experiments were conducted in accordance with the United Kingdom Animal (Scientific Procedures) Act 1986, and approved by the Animals in Scientific Research Unit at the Home Office (licence PPL70-13691). C57BL/ mice were bred (2–4-month-old males), whereas WAG rats were imported (Charles River, 6–9 months of age). Animals were housed in a 12 h/12 h dark/light cycle, and food and water were given ad libitum. The animals were group-housed prior to headbar surgery, but after this, animals were individually housed.
Acute preparation surgeries for headbar attachment and craniotomy
For both surgeries, aseptic techniques were used with mice anaesthetized using isoflurane (988–3245, Henry Schein) and placed in a stereotaxic frame (David Kopf Instruments). Viscotears were applied (Bausch + Lomb) and pain relief, which consisted of subcutaneous bupenorphine (0.5 mg per kg mouse (hereafter denoted mg kg–1 mouse); Ceva) and metacam (15 mg kg–1 mouse; Boehringer Ingelheim), was injected. Saline was administered just before recovery or every 45 min depending on the length of surgery. To apply the headbars for the Neurotar system, the skin on the top of the head was cut to expose the skull. The skull was cleaned and dried, which enabled drilling (RA1 008, Henry Schein) of a small hole in the left-hand visual cortex for a metal support screw (00-96 × 3-32, Plastics One). The headplate (Model 9, Neurotar) was firmly attached using vetbond (1469SB, 3M) and strengthened using dental cement before the exposed skull was covered with Kwik-Cast (WPI). Mice were checked daily to ensure recovery. After at least 5 days of recovery, habituation was performed by placing the mouse in the Neurotar frame for increasing periods of time (15–60 min) over several days. On the day of recording, a craniotomy was performed. Under isoflurane anaesthesia, with administration of pain medication and intramuscular dexamethasone (1 mg per kg of mouse, intramuscular; 7247104, MSD Animal Health), two areas were exposed. A large (2 × 2 mm2) craniotomy was performed over the somatosensory and visual cortex on the right-hand side and a small drill hole over the motor cortex on the left-hand side. Cold Cortex-buffered saline was continually applied to the craniotomies. After completion, the exposed dura was covered with Cortex-buffered saline, sterilized Sylgard (~200-μm thick) and a Kwik-Cast layer. After ~2 h, the animal was moved to the Neurotar frame and the craniotomies were exposed by removal of the Kwik-Cast and Sylguard layers. The gDNP was carefully connected to a printed circuit board (PCB) and lowered using a micromanipulator to just above the dura over the visual cortex. The dura was gently pierced using either microdissection scissors or a 26-gauge needle and the gDNP lowered ~2 mm into the brain. A reference wire (Ag/AgCl2) was placed in the ipsilateral motor cortex and g.tec hardware (see the section Electronics for in vivo recordings with gDNPs) was used to generate the transfer curve to determine the optimal VGS and initiate recordings (Supplementary Fig. 4). Chemoconvulsant was injected into the brain using a Nanofil injection system (WPI). At the end of the experiment, sodium pentobarbital was administered intraperitoneally.
Recording with solution-filled glass micropipette
Borosilicate capillary tubes (outer diameter 1.50 mm, inner diameter 0.86 mm; Warner Instruments) were pulled using a horizontal puller (Sutter Instruments P-97, resistance of 3–5 MΩ), filled with artificial cerebral spinal fluid (119 mM NaCl, 2.5 mM KCl, 2.5 mM CaCl2, 1.3 mM MgSO4, 1.25 mM NaH2PO4, 25 mM NaHCO3 and 10 mM glucose) and attached to an CV-7B headstage (Axon Instruments). A micromanipulator (MM3301, WPI) was used to position the pipette above the cortical surface before insertion approximately 400 µm into the cortex. The headstage was provided with the same reference as the gDNP, a chlorinated silver wire touching the ipsilateral motor cortex. The headstage was connected to a Multiclamp 700B amplifier (Axon Instruments) operating in current clamp mode. Analogue–digital conversion and transistor–transistor logic pulse delivery for temporal synchronization was achieved using a Micro1401 MkII unit (CED). Data were acquired by WINEDR sampling at 20 kHz with a 4 kHz Bessel filter.
Chronic preparation surgery and recording
First, the gDNP was coated with fibroin, as described above, to aid insertion. The rat was anaesthetized to a surgical depth using Isoflurane. After placement in a stereotaxic frame, Viscotears were applied and pain medication, which consisted of bupenorphine (0.15 mg kg–1 mouse, subcutaneous; Ceva) and metacam (4.5 mg kg–1 mouse, subcutaneous; Boehringer Ingelheim), was injected. The skull was cleaned and dried. Small burr holes (~1 mm) were drilled at four positions: (1) the somatosensory cortex for gDNP (as the perioral somatosensory cortex is the focal area for SWDs56), (2) the contralateral cerebellum for a reference Ag/AgCl wire held in place by a nylon screw, (3) the motor cortex, ipsilateral, for a support screw and (4) the visual cortex, ipsilateral, for a support screw. The metal screws were inserted and provided structural support for the dental cement. Next, the gDNP and reference wire were inserted, and a DC characterization curve confirmed that the transistors were performing optimally. Dental cement, mixed with vetbond, was applied around the PCB for support. Animals were weighed daily and their physiology was monitored to ensure a full recovery. For recording, animals were anaesthesia-free and moving, with the PCB interface on the head connected to an Omnetics cable (A79635, Omnetics) that interfaced with the g.tec recording hardware as described above. After allowing ~5 min for settling, a DC characterization curve was recorded to allow accurate calibration of the gSGFETs. A script was used to calculate the optimal VGS based on the transfer curves. Recordings were performed for ~10–60 min twice a week for 10 weeks. After recording, the Omnetics wire was disconnected and a protective cap was applied.
WAG/Rij strain of rats
WAG/Rij rats were used as a model of absence epilepsy. Spike-wave activity in the cortical EEG of adult WAG/Rij rats has a frequency of 7–10 Hz and can last from just a few, to up to 40 s in duration. From an age of 6 months, WAG/Rij rats display several hundred SWDs a day20. In addition to electrophysiological signs, mild behavioural phenotypes concurrent with SWD similar to human absence epilepsy are displayed. These include facial myoclonic jerks, twitching, accelerated breathing, head tilting and eye blinking. Drugs prescribed for human absence epilepsy suppress SWDs in WAG/Rij rats, and these rats have been used to predict the action of antiepileptic medication and other experimental treatments57.
Data analysis
All electrophysiological data were analysed using Python v.3.7 packages (Matplotlib, Numpy, Neo and Elephant) and the custom library PhyREC (https://github.com/aguimera/PhyREC). The conversion of the recorded current signals (LPF and BPF) to a voltage signal was performed by summation of the two signals and inverse interpolation in the in vivo/chronic measured transfer curve of the corresponding gSGFET17. Even for large signals (tens of millivolts), this calibration method yields distortion-free signals, in contrast to using a fixed transconductance value30. Transfer curves were always measured at the beginning and end of every recording to ensure that no significant variations were present and to detect any malfunctioning transistor. Moreover, all recordings presented in the manuscript have been calibrated with the nearest-recorded transfer curve to ensure high fidelity in the voltage-converted signals.
To assess the detection limit in the conventional frequency bands (>1 Hz), baseline recordings were directly compared with post-mortem recordings (Fig. 2f). For instance, applying a digital filter (>200 Hz) and comparing post-mortem with baseline validates the ability of gDNP to record spontaneous high-frequency activity (>200 Hz) arising from groups of neurons in a non-pathological brain state (upper plot in Fig. 2f). Beyond this qualitative comparison, the SNR shown in Fig. 2f was evaluated from the ratio of the RMS mean value over 25 s of recording of baseline (spontaneous activity) and post-mortem (no activity) signals. The signals were BPF in three different bands, corresponding to local field potential activity (1–70 Hz), high-frequency activity (80–200 Hz) and very high-frequency activity (200–4,000 Hz). RMS values were calculated with a sliding window of 500 ms for the 1–70 Hz band and with a sliding window of 10 ms for the other two bands. The whisker plots were created with a Python library (seaborn.boxplot) to show the data distribution for working transistors in the gDNP (gDNP1 = 10, gDNP2 = 12, gDNP3 = 12 and gDNP4 = 14). In the whisker plots, the line in the box corresponds to the median value, the box shows the quartiles of the dataset and the whiskers extend to show the maximum and minimum values. The SNRs for the different bands were evaluated from a total of four in vivo experiments with four different gDNPs (Fig. 2f). The SNR is expressed in decibels (20 × ln(RMS(S)/RMS(N)), where S is the ‘signal’ during the baseline recording and N is the ‘noise’ measured during the post-mortem recording. The silencing of neuronal activity shown in Fig. 3d was extracted from the AC-coupled recording (HP-filtered, >0.5 Hz). Then, the RMS values of the pre-ictal phase (calculated with a sliding window of 1 s) were averaged over 50 s. A similar analysis was performed during the SD (15 s, shaded areas in Fig. 3d). The ratio of the two averaged RMS values corresponds to the neuronal activity variation (before and during the SD). The amplitudes of the hippocampal SD and the hyperpolarization wave in Fig. 3d were evaluated using the recording LPF in the infraslow regime (<0.5 Hz) and resampled at 3 Hz (instead of at 9.6 kHz, used in the original recording). The zero of the voltage was set using the mean value of the signal 50 s before the pre-ictal phase (Supplementary Fig. 9), and the minimum and maximum values for each channel were extracted (corresponding to the SD and hyperpolarization amplitude, respectively).
For the analysis of the DC-shifts pre-seizure (Fig. 3f and Supplementary Fig. 17), 20 min recordings starting from 1 min before the first seizure following 4-AP injection were selected. Seizures with concurrent high-amplitude depolarizations such as SDs were excluded from analysis. The 20 min recordings were filtered in two bands: DC-coupled (LPF, <0.1 Hz) and seizure band (bandpass second-order filter, 4–12 Hz). Signals were downsampled to 96 Hz. The RMS value of the seizure band was calculated with a 0.5 s sliding window and used to determine seizure onset and end. Epochs around −5 to 2 s relative to seizure start of the DC-coupled signal were extracted for analysing the DC-shifts. Zero voltage was set to the mean of the first 2 s of the epoch, and the maximum amplitude of the period from −5 to 0 s was then extracted and used as the amplitude of the pre-seizure DC-shift.
CSD analysis applied to the low-frequency part of the potential (local field potential) was calculated with the Python open source Elephant library (Elephant Electrophysiology Analysis Toolkit) using the class ‘Current Source Density analysis (CSD)’, and the method 1D – StandardCSD was chosen for the linear gDNP array. A homogeneous conductivity of the neural tissue of σ = 0.3 S m–1 across the different layers was used for the calculations.
The concurrency of ISO with the SWD shown in Fig. 4h was evaluated by performing the Hilbert transform to extract the phase of the DC-coupled signal (0.01–0.1 Hz) and associating it with the RMS in the range 5–9 Hz evaluated with a sliding window of 0.5 s (typical bandwidth for the SWD). Supplementary Fig. 16 shows in more detail the dependency of the ISO phase and SWD amplitude. The data shown in Fig. 4h illustrate the joint density probability, indicating the higher probability of a given RMS amplitude to ISO phase pair.
Chronic biocompatibility study
Details of the surgery, tissue collection, analysis of histology and data statistics used for the chronic biocompatibility study are provided in Supplementary Figs. 19 and 20.
Reporting Summary
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
Data availability
All relevant data obtained to evaluate the main findings of the paper are openly available in Zenodo at https://doi.org/10.5281/zenodo.5655535. All other raw data are available from the corresponding author upon reasonable request.
Code availability
Python v.3.7 packages (Matplotlib v.3.2.0 and Numpy v.1.17.4) and the following Python library were used for electric characterization of the gSGFET arrays: https://github.com/aguimera/PyGFET. A custom Simulink model was used for graphene microtransistor electrophysiological data acquisition; contact g.tec medical engineering for code access. Electrophysiological data were analysed using Python v.3.7 packages (Matplotlib v.3.2.0, Numpy v.1.17.4, Pandas v.0.25.3, seaborn v.0.9.0, Neo v.0.8.0 and Elephant) and the custom library PhyREC (https://github.com/aguimera/PhyREC/tree/PhyREC4). Custom scripts can be found at Zenodo (https://doi.org/10.5281/zenodo.5655535). Immunohistochemical data analysis was performed using Python v.3.7 script (https://github.com/kebarr/biocompatibility_study).
References
Perucca, P., Dubeau, F. & Gotman, J. Intracranial electroencephalographic seizure-onset patterns: effect of underlying pathology. Brain 137, 183–196 (2014).
Modur, P. N. High frequency oscillations and infraslow activity in epilepsy. Ann. Indian Acad. Neurol. 17, S99–S106 (2014).
Revankar, G. S. et al. in Seizures in Critical Care: A Guide to Diagnosis and Therapeutics (eds Varelas, P. N. and Claassen, J.) 77–90 (Humana Press, 2017); https://doi.org/10.1007/978-3-319-49557-6_5
Dreier, J. P. The role of spreading depression, spreading depolarization and spreading ischemia in neurological disease. Nat. Med. 17, 439–447 (2011).
Staba, R. J., Stead, M. & Worrell, G. A. Electrophysiological biomarkers of epilepsy. Neurotherapeutics 11, 334–346 (2014).
Dell, K. L., Cook, M. J. & Maturana, M. I. Deep brain stimulation for epilepsy: biomarkers for optimization. Curr. Treat. Options Neurol. 21, 47 (2019).
Kuhlmann, L., Lehnertz, K., Richardson, M. P., Schelter, B. & Zaveri, H. P. Seizure prediction—ready for a new era. Nat. Rev. Neurol. 14, 618–630 (2018).
Chari, A., Thornton, R. C., Tisdall, M. M. & Scott, R. C. Microelectrode recordings in human epilepsy: a case for clinical translation. Brain Commun. 2, fcaa082 (2020).
Lee, S. et al. DC shifts, high frequency oscillations, ripples and fast ripples in relation to the seizure onset zone. Seizure 77, 52–58 (2020).
Li, C. et al. Evaluation of microelectrode materials for direct-current electrocorticography. J. Neural Eng. 13, 16008 (2015).
Hartings, J. A. How slow can you go? Nat. Mater. 18, 194–196 (2019).
Major, S., Gajovic-Eichelmann, N., Woitzik, J. & Dreier, J. P. Oxygen-induced and pH-induced direct current artifacts on invasive platinum/iridium electrodes for electrocorticography. Neurocrit. Care 35, 146–159 (2021).
Khodagholy, D. et al. In vivo recordings of brain activity using organic transistors. Nat. Commun. 4, 1575 (2013).
Kostarelos, K., Vincent, M., Hebert, C. & Garrido, J. A. Graphene in the design and engineering of next-generation neural interfaces. Adv. Mater. 29, 1700909 (2017).
Blaschke, B. M. et al. Mapping brain activity with flexible graphene micro-transistors. 2D Mater. 4, 025040 (2017).
Hébert, C. et al. Flexible graphene solution-gated field-effect transistors: efficient transducers for micro-electrocorticography. Adv. Funct. Mater. 28, 1703976 (2017).
Masvidal-Codina, E. et al. High-resolution mapping of infraslow cortical brain activity enabled by graphene microtransistors. Nat. Mater. 18, 280–288 (2019).
Weltman, A., Yoo, J. & Meng, E. Flexible, penetrating brain probes enabled by advances in polymer microfabrication. Micromachines 7, 180 (2016).
Tien, L. W. et al. Silk as a multifunctional biomaterial substrate for reduced glial scarring around brain-penetrating electrodes. Adv. Funct. Mater. 23, 3185–3193 (2013).
Coenen, A. M. L. & Van Luijtelaar, E. L. J. M. Genetic animal models for absence epilepsy: a review of the WAG/Rij strain of rats. Behav. Genet. 33, 635–655 (2003).
Terlau, J. et al. Spike-wave discharges in absence epilepsy: segregation of electrographic components reveals distinct pathways of seizure activity. J. Physiol. https://doi.org/10.1113/JP279483 (2020).
Zijlmans, M. et al. High-frequency oscillations as a new biomarker in epilepsy. Ann. Neurol. 71, 169–178 (2012).
Jacobs, J. et al. High-frequency oscillations (HFOs) in clinical epilepsy. Prog. Neurobiol. 98, 302–315 (2012).
Ikeda, A. et al. Focal ictal direct current shifts in human epilepsy as studied by subdural and scalp recording. Brain 122, 827–838 (1999).
Wu, S. et al. Role of ictal baseline shifts and ictal high-frequency oscillations in stereo-electroencephalography analysis of mesial temporal lobe seizures. Epilepsia 55, 690–698 (2014).
Vanhatalo, S. et al. Very slow EEG responses lateralize temporal lobe seizures: an evaluation of non-invasive DC-EEG. Neurology 60, 1098–1104 (2003).
Duan, X. et al. Intracellular recordings of action potentials by an extracellular nanoscale field-effect transistor. Nat. Nanotechnol. 7, 174–179 (2012).
Hess, L. H. Graphene transistors for biosensing and bioelectronics. Proc. IEEE 101, 1780–1792 (2013).
Hess, L. H. et al. High-transconductance graphene solution-gated field effect transistors. Appl. Phys. Lett. 99, 033503 (2011).
Garcia-Cortadella, R. et al. Distortion-free sensing of neural activity using graphene transistors. Small 16, 1906640 (2020).
Lecomte, A. et al. Silk and PEG as means to stiffen a parylene probe for insertion in the brain: toward a double time-scale tool for local drug delivery. J. Micromech. Microeng. 25, 125003 (2015).
Fueta, Y. & Avoli, M. Effects of antiepileptic drugs on 4-aminopyridine-induced epileptiform activity in young and adult rat hippocampus. Epilepsy Res. 12, 207–215 (1992).
Padmanabhan, K. & Urban, N. N. Disrupting information coding via block of 4-AP-sensitive potassium channels. J. Neurophysiol. 112, 1054–1066 (2014).
Zakharov, A., Chernova, K., Burkhanova, G., Holmes, G. L. & Khazipov, R. Segregation of seizures and spreading depolarization across cortical layers. Epilepsia 60, 2386–2397 (2019).
Hartings, J. A. et al. Direct current electrocorticography for clinical neuromonitoring of spreading depolarizations. J. Cereb. Blood Flow Metab. 37, 1857–1870 (2017).
Harriott, A. M., Takizawa, T., Chung, D. Y. & Chen, S. P. Spreading depression as a preclinical model of migraine. J. Headache Pain 20, 45 (2019).
Buzsáki, G. & Lopes da Silva, F. L High frequency oscillations in the intact brain. Prog. Neurobiol. 98, 241–249 (2012).
Ikeda, A. et al. Active direct current (DC) shifts an “Red slow”: two new concepts for seizure mechanisms and identification of the epileptogenic zone. Neurosci. Res. 156, 95–101 (2020).
Kamarajan, C., Pandey, A. K., Chorlian, D. B. & Porjesz, B. The use of current source density as electrophysiological correlates in neuropsychiatric disorders: a review of human studies. Int. J. Psychophysiol. 97, 310–322 (2015).
Flynn, S. P., Barrier, S., Scott, R. C., Lenck-Santini, P. P. & Holmes, G. L. Status epilepticus induced spontaneous dentate gyrus spikes: in vivo current source density analysis. PLoS ONE 10, e0132630 (2015).
Coenen, A. M. L. & Van Luijtelaar, E. L. J. M. The WAG/Rij rat model for absence epilepsy: age and sex factors. Epilepsy Res. 1, 297–301 (1987).
Orlowska-Feuer, P. et al. Infra-slow modulation of fast beta/gamma oscillations in the mouse visual system. J. Physiol. 599, 1631–1650 (2021).
Garcia-Cortadella, R. et al. Graphene active sensor arrays for long-term and wireless mapping of wide frequency band epicortical brain activity. Nat. Commun. 12, 211 (2021).
Bahari, F. et al. Seizure-associated spreading depression is a major feature of ictal events in two animal models of chronic epilepsy. Preprint at bioRxiv https://doi.org/10.1101/455519 (2020).
Dreier, J. P. et al. Spreading convulsions, spreading depolarization and epileptogenesis in human cerebral cortex. Brain 135, 259–275 (2012).
Dreier, J. P. et al. Recording, analysis, and interpretation of spreading depolarizations in neurointensive care: review and recommendations of the COSBID research group. J. Cereb. Blood Flow Metab. 37, 1595–1625 (2017).
De Tisi, J. et al. The long-term outcome of adult epilepsy surgery, patterns of seizure remission, and relapse: a cohort study. Lancet 378, 1388–1395 (2011).
Kanazawa, K. et al. Intracranially recorded ictal direct current shifts may precede high frequency oscillations in human epilepsy. Clin. Neurophysiol. 126, 47–59 (2015).
Lauritzen, M. et al. Clinical relevance of cortical spreading depression in neurological disorders: migraine, malignant stroke, subarachnoid and intracranial hemorrhage, and traumatic brain injury. J. Cereb. Blood Flow Metab. 31, 17–35 (2011).
Schaefer, N. et al. Improved metal-graphene contacts for low-noise, high-density microtransistor arrays for neural sensing. Carbon 161, 647–655 (2020).
Jin, H. J. et al. Water-stable silk films with reduced β-sheet content. Adv. Funct. Mater. 15, 1241–1247 (2005).
Rockwood, D. N. et al. Materials fabrication from Bombyx mori silk fibroin. Nat. Protoc. 6, 1612–1631 (2011).
Vepari, C. & Kaplan, D. L. Silk as a biomaterial. Prog. Polym. Sci. 32, 991–1007 (2007).
Cao, Y. & Wang, B. Biodegradation of silk biomaterials. Int. J. Mol. Sci. 10, 1514–1524 (2009).
Gobin, A. S., Froude, V. E. & Mathur, A. B. Structural and mechanical characteristics of silk fibroin and chitosan blend scaffolds for tissue regeneration. J. Biomed. Mater. Res. A 74, 465–473 (2005).
Russo, E. et al. Upholding WAG/Rij rats as a model of absence epileptogenesis: hidden mechanisms and a new theory on seizure development. Neurosci. Biobehav. Rev. 71, 388–408 (2016).
van Luijtelaar, G. & van Oijen, G. Establishing drug effects on electrocorticographic activity in a genetic absence epilepsy model: advances and pitfalls. Front. Pharmacol. 11, 395 (2020).
Acknowledgements
This work has received funding from the European Union’s Horizon 2020 research and innovation programme under Grant Agreement No. 881603 (GrapheneCore3). ICN2 is supported by the Severo Ochoa Centres of Excellence programme, funded by the Spanish Research Agency (AEI, grant no. SEV-2017–0706), and by the CERCA Programme/Generalitat de Catalunya. A.B.C. is supported by the International PhD Programme La Caixa-Severo Ochoa (Programa Internacional de Becas ‘la Caixa’-Severo Ochoa). This work has made use of the Spanish ICTS Network MICRONANOFABS, partially supported by MICINN and the ICTS ‘NANBIOSIS’, more specifically by the Micro-NanoTechnology Unit of the CIBER in Bioengineering, Biomaterials and Nanomedicine (CIBER-BBN) at the IMB-CNM. We also acknowledge funding from the Generalitat de Catalunya (2017 SGR 1426), and the 2DTecBio project (FIS2017-85787-R) funded by the Ministerio de Ciencia, Innovación y Universidades of Spain, the Agencia Estatal de Investigación (AEI) and the Fondo Europeo de Desarrollo Regional (FEDER/UE). Part of this work was co-funded by the European Regional Development Funds (ERDF) allocated to the Programa operatiu FEDER de Catalunya 2014–2020, with the support of the Secretaria d’Universitats i Recerca of the Departament d’Empresa i Coneixement of the Generalitat de Catalunya for emerging technology clusters devoted to the valorization and transfer of research results (GraphCAT 001-P-001702). A.B.C. acknowledges that this work has been carried out within the framework of a PhD in Electrical and Telecommunication Engineering at the Universitat Autònoma de Barcelona. R.C.W. is funded by a Senior Research Fellowship awarded by the Worshipful Company of Pewterers. D.R. is a Biotechnology and Biological Sciences Research Council (BBSRC) LIDo sponsored PhD student. We thank M. Walker and L. Lemieux (UCL Queen Square Institute of Neurology) for their comments on the manuscript.
Author information
Authors and Affiliations
Contributions
A.B.C. carried out most of the fabrication and characterization of the gDNPs, contributed to the design and performance of the in vivo experiments, analysed the data and wrote the manuscript. E.M.-C. contributed to the design and planning of the in vivo experiments, to the data analysis and particularly to the DC-shift and SD analysis of the in vivo data. R.C.W. and T.M.S. performed the in vivo experiments. D.R. contributed to the in vivo experiments and DC-coupled recordings with the glass micropipette. N.S., E.R.-L., X.I. and J.M.D.l.C. contributed to the fabrication and characterization of the gDNPs. E.D.C., J.B. and C.H. contributed to the growth, transfer and characterization of the CVD graphene used in the gDNPs. E.P.-A., A.H. and E.R.-L. contributed to the optimization of the SF stiffening protocol of the gDNPs. J.M.-A. contributed to the fabrication of the custom electronic instrumentation and development of a Python-based user interface. D.V. contributed to the Python scripts and technical discussions. J.R.S. reviewed the manuscript. J.F. and J.S. contributed to the mechanical assessment of the SF and the SF back-coated gDNPs. M.D. performed all surgeries for the biocompatibility study. A.D. and K.B. contributed to the capture of histological images and image processing and analysis. S.S. and K.B. contributed to the preparation and review of the manuscript. A.G.-B. contributed to the design and fabrication of the custom electronic instrumentation, the development of a custom gSGFET Python library and analysis of the data. R.V., K.K., R.C.W., A.G.-B. and J.A.G. participated in the design of all experiments and thoroughly reviewed the manuscript. All authors read and reviewed the manuscript.
Corresponding authors
Ethics declarations
Competing interests
C.G. is the owner of g.tec medical engineering and Guger Technologies. J.A.G, K.K and A.G.-B declare financial interest in INBRAIN Neuroelectronics. All other authors have no competing interests.
Additional information
Peer review information Nature Nanotechnology thanks Matthew Nelson, Bozhi Tian and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Supplementary Information
Supplementary Methods, Figs. 1–23 and Table 1.
Rights and permissions
About this article
Cite this article
Bonaccini Calia, A., Masvidal-Codina, E., Smith, T.M. et al. Full-bandwidth electrophysiology of seizures and epileptiform activity enabled by flexible graphene microtransistor depth neural probes. Nat. Nanotechnol. 17, 301–309 (2022). https://doi.org/10.1038/s41565-021-01041-9
Subjects
This article is cited by
-
Clinical translation of graphene-based medical technology
Nature Reviews Electrical Engineering (2024)
-
High-density transparent graphene arrays for predicting cellular calcium activity at depth from surface potential recordings
Nature Nanotechnology (2024)
-
Control of polymers’ amorphous-crystalline transition enables miniaturization and multifunctional integration for hydrogel bioelectronics
Nature Communications (2024)
-
Graphene-integrated mesh electronics with converged multifunctionality for tracking multimodal excitation-contraction dynamics in cardiac microtissues
Nature Communications (2024)
-
The ultra-thin, minimally invasive surface electrode array NeuroWeb for probing neural activity
Nature Communications (2023)