L-type voltage-gated calcium channel regulation of in vitro human cortical neuronal networks

The combination of in vitro multi-electrode arrays (MEAs) and the neuronal differentiation of stem cells offers the capability to study human neuronal networks from patient or engineered human cell lines. Here, we use MEA-based assays to probe synaptic function and network interactions of hiPSC-derived neurons. Neuronal network behaviour first emerges at approximately 30 days of culture and is driven by glutamate neurotransmission. Over a further 30 days, inhibitory GABAergic signalling shapes network behaviour into a synchronous regular pattern of burst firing activity and low activity periods. Gene mutations in L-type voltage gated calcium channel subunit genes are strongly implicated as genetic risk factors for the development of schizophrenia and bipolar disorder. We find that, although basal neuronal firing rate is unaffected, there is a dose-dependent effect of L-type voltage gated calcium channel inhibitors on synchronous firing patterns of our hiPSC-derived neural networks. This demonstrates that MEA assays have sufficient sensitivity to detect changes in patterns of neuronal interaction that may arise from hypo-function of psychiatric risk genes. Our study highlights the utility of in vitro MEA based platforms for the study of hiPSC neural network activity and their potential use in novel compound screening.

Understanding how neurons interact both within local networks or as components of neural circuits is a key objective for the investigation of brain function throughout development and in neuropsychiatric disorders. Study of neural interactions in animals in vivo, in intact brain slices or in isolated primary neurons in culture has been a mainstay of neuroscience for more than a century. This has been augmented by human studies investigating functional (EEG and MEG) and structural (MRI) connectivity. However, the study of the development and function of human neural networks at the cellular level has been more difficult, largely due to the difficulty of combining cell culture methods and longitudinal electrophysiological experimentation. This is now changing with the opportunity to develop neurons derived from human induced pluripotent stem cells (hiPSC) and the advent of user-friendly multi-electrode array (MEA) systems. These offer the potential to generate small-scale human neural networks within a familiar context and investigate their development and function in vitro 1 .
Dissociated primary rodent neurons in vitro self-organise to form functional networks [2][3][4] , which spontaneously exhibit coordinated action potential bursts (synchronised bursts; SBs) that can be detected across spatially separated electrodes of an MEA 5,6 . This culture-wide, synchronised bursting behaviour correlates with an increased basal firing rate and emerges after an approximate two-week period of asynchronous firing. Recently, it has been shown that this synchronous behaviour may develop further to form a temporal pattern of regular low and high activity periods lasting for tens of seconds 7 , behaviour which is sensitive to pharmacological manipulation and in particular AMPA receptor activity.
MEA methods to study neuronal networks derived from hiPSC are less established. However, a number of MEA studies have shown that human stem cell-derived neurons develop spontaneous activity 8 , which increases during an extended period of development [9][10][11][12] . With prolonged incubation, the emergence of short SBs across multiple electrodes of the array has been observed 10 . These studies have established the use of MEAs for pharmacological profiling and toxicology testing 8,10,13 . However, it is currently not known whether such methods have the discriminatory sensitivity to investigate the effects of genetic risk in patient hiPSCs derived cultures, thereby offering a platform for drug development.
In this study, we demonstrate the use of MEA assays to probe synaptic function and functional neural networks in hiPSC-derived neuronal cultures. We first investigate the development of network behaviour in in vitro neuronal cultures on MEAs and use pharmacology to demonstrate a need for synaptic activity. We then use our hiPSC-based MEAs to probe the roles of L-type voltage-gated calcium channel activity in network behaviour. These channels have been implicated in several aspects of network function in vivo, including the synchronisation of calcium events 14,15 and the regulation of hippocampal oscillatory activity 16,17 . Significantly, genetic variants in CACNA1C and CACNA1D (coding for the alpha-1 subunits Ca v 1.2 and Ca v 1.3 respectively) are some of the most robust genetic risk alleles associated with the development of mental health disorders, including schizophrenia and bipolar disorder. By use of selective L-type channel inhibitors we demonstrate a specific modulation of SB pattern in hiPSC-derived neural networks without changes in basal neuronal activity. These results demonstrate the potential of hiPSC-based MEA assays to detect functional changes arising from genetic risk loci for psychiatric disorders.

Results
creation of hipSc-derived neuronal networks. Neurons were differentiated from hiPSC by a dual-SMAD monolayer protocol and re-plated on MEAs at 30DIV. MEA cultures were maintained in astrocyte conditioned medium (ACM) and in a 2% O 2 atmosphere to enhance neuronal maturation 18,19 . Parallel cell cultures to those on MEAs were examined by immunocytochemistry at 50 days post re-plating (50 DPP; corresponding to 50 days after cells were re-plated on MEAs or 80 days of total differentiation; Fig. 1). All imaged neurons expressed the neuron-specific microtubule associated protein MAP2, of which 64 ± 9.0% co-expressed the vesicular glutamate transporter vGLUT1 and 3.6 ± 1.2% of cells expressed the GABA synthesising enzyme GAD67 (Fig. 1A,D), indicating a mixed culture of predominately glutamatergic neurons with a small population of GABAergic cells. We also found that 73.5 ± 9.7% of nuclei expressed CTIP2 with 37.7 ± 8.4% also expressing SATB2 (Fig. 1B,D), suggesting a deeplayer cortical identity of most neurons primarily corresponding to layers 6 and 5a. Higher magnification images revealed that neurons also expressed the key excitatory synaptic scaffolding protein PSD95 and the obligatory glycine-binding NMDA receptor subunit, GluN1, suggesting the presence of excitatory synapses (Fig. 1C).
To confirm functional maturation of our hiPSC-derived neurons, we performed single-cell patch clamp experiments ( Fig. 2A) to examine the basal physiological state of the neurons. Resting membrane potential (V m ) and input resistance of neurons significantly reduced throughout differentiation (−30.6 ± 10.2 mV to −44.8 ± 10.0 mV and 1.26 ± 0.66 GΩ to 0.92 ± 0.41 GΩ respectively), indicating that neurons developed over time (Fig. 2B), although no change of the membrane time constant (tau) was observed. Patched neurons were classified according to their ability to fire induced and spontaneous action potentials (iAP and sAPs respectively; Fig. 2A,C). Cells were held at −70 mV via constant current injection and iAPs were elicited by injecting additional depolarising current steps. Neurons were subsequently coded based upon their voltage response ( Fig. S1 and methods). At 50 DPP, 83.9% of neurons were classed as active, which included neurons showing single iAPs, attempts at trains or full iAP trains ( Fig. 2A). The AP threshold and average event amplitude of these iAPs was unchanged over development, however neurons did exhibit AP with faster rise times and shorter half-widths at 50 DPP compared to 30 DPP (Fig. 2B), highlighting the maturation of the cells over these 20 days. Finally, at 50 DPP, 22.9% of patched neurons exhibited sAPs, as determined by continuous recording of cells at their resting membrane potential (I = 0; Fig. 2C).
Our single cell electrophysiology and imaging experiments indicated that our protocols yielded mature and functionally active neurons. To examine the development of networks in the cultures, differentiated hiPSC-derived neurons were plated as drop cultures onto MEAs possessing an 8 × 8 grid of substrate embedded electrodes, spaced 200 μm apart. Spike detection from filtered raw voltage recording made from each electrode utilised the threshold-based method of QuianQuiroga and co-workers 20 and following quality control spikes were timestamped for subsequent analysis. Spontaneous activity of cultures changed markedly over the 60 days of recorded development (Fig. 3A). Initial activity was very low with a mean firing rate of 0.065 ± 0.021 Hz and no detectable spike bursts but increased to a firing rate of 0.75 ± 0.21 Hz by 20DPP with 232 ± 103 bursts. The most substantial increase in activity occurred between 20 to 30 DPP with a 2-3-fold increase in spike rate and burst number, while after 30 DPP basal activity levels remain largely stable. These results are consistent with an increase in excitability previously seen throughout hiPSC-derived culture development [8][9][10][11] .
The emergence of synchronised bursting (SB) coincided with the switch to high neuronal activity levels seen at 30 DPP (Fig. 3A). This change in behaviour can be seen more clearly in raster and array-wide spike detection rate (ASDR) plots (Fig. 3B,C). At 30 DPP, SBs were short in length (1.25 ± 0.37 s) with a firing rate of 107.60 ± 87.06 Hz (Fig. 3D). Over the next 30 days, the number of SBs decreased, with a concomitant increase in both the length of the bursts and the interval between them (17.07 ± 8.34 s and 78.33 ± 33.00 s respectively at 60 DPP). Although SB firing can be detected from 30 DPP, there was a clear shift in the pattern of array-wide activity from 50 DPP into a regular pattern of more active periods (MAPs) and less active periods (LAPs) each lasting >30 seconds. Once networks have entered this state, the MAP/LAP pattern evolves further by increasing in MAP (SB) and LAP (SB interval) lengths (Fig. 3D). Finally, it should be noted that the changes in synchronised activity patterns occur with little change of overall basal neuronal activity (Fig. 3A).
Pharmacological profiling of hiPSC derived neuronal networks. To establish whether the extended SB firing patterns arise via synaptic activity, we probed our MEA cultures with pharmacological agents that alter glutamate/GABA signalling. At 20 DPP, when no SB activity is seen, neuronal cultures are unresponsive to AMPA, NMDA and GABA A receptor inhibition, but sensitive to 10 μM GABA, which strongly attenuated spontaneous activity (Fig. S2). At the later 50 DPP stage when sustained SB firing has fully matured, inhibition of glutamate signalling by application of either the AMPA receptor antagonist CNQX or the NMDA receptor antagonist APV completely suppressed synchronised bursting (Fig. 4). In both cases, SB patterns rapidly re-establisherd after drug washout. Importantly, neither drug had a substantial effect on the basal rate of neuronal firing (Spike Rate) (Fig. 4Ca), indicating that both AMPA and NMDA receptor mediated neurotransmission is required for network synchrony, not single cell firing. Although the hiPSC differentiation protocol used here generates primarily cortical glutamatergic neurons, the resulting cultures also contain a small proportion on GABAergic neurons, which, in the nervous system, would be expected to exert an inhibitory effect on neural circuits and networks. To investigate the potential role of this inhibitory signalling on SB patterns, GABA or GABA A receptor antagonists were applied acutely to the MEA cultures. Exposure to 1 μM GABA caused an almost complete, but reversible, loss of both basal and hence synchronised neuronal activity (Fig. 4), demonstrating that our hiPSC-derived glutamatergic neurons are GABA sensitive. In contrast, lowering the effect of endogenous GABAergic signalling within the culture by application of bicuculline, a GABA A receptor blocker, significantly increased the number of SBs recorded (t 10 = 2.81; p < 0.05). This was accompanied by a concomitant decrease in SB firing rate (t 10 = 1.91; p < 0.05) and in SB interval length ( Fig. 4C; t 10 = 2.14; p < 0.05). As a consequence, the general spike rate was not significantly altered, so that in effect reducing the effect of GABA increased the SB frequency without decreasing neuronal activity. These changes were reproduced when cultures were exposed to the non-competitive GABA A antagonist, picrotoxin ( Fig. S3), strongly suggesting that this regulation is driven by the attenuation of GABAergic function.

Synchronised network activity in MeA cultures is regulated by L-type calcium channels.
Calcium signalling via L-type voltage gated calcium (VGCCs) has been shown to be involved in the regulation of oscillatory behaviour in certain brain regions 21,22 . NMDA-mediated network signalling is thought to be regulated via calcium homeostasis mechanisms, partly through L-type VGCCs 23,24 . Moreover, recent genetic studies have strongly implicated L-type VGCCs in increased risk for a range of neuropsychiatric disorders, including schizophrenia, autism, depression and epilepsies; all pathologies with strong evidence for aberrant network signalling 25,26 . To investigate further the physiology of the neuronal cultures and to study the potential role of calcium in the underlying mechanism contributing to the network activity observed here, cells were exposed to the L-type VCCC blocker, Diltiazem.
Acute exposure of increasing doses of Diltiazem in 60 DPP day cultures had little effect on the basal firing properties of the neurons, specifically spike rate and single-unit bursts (Fig. 5Ba,Bb). However, application of Diltiazem induced a dose -dependent increase in the number of synchronised bursts (SB) and corresponding reduction in the interval between SBs, decreasing from 62.6 ± 37.3 s, to 39.5 ± 11.4 s (1 μM), to 22.9 ± 12.7 s (2 μM) and finally to 20.5 ± 9.4 s (5 μM; Fig. 5Bc,Bd). Dunnett's multiple comparison tests following one-way ANOVA (F = 6.004, R 2 = 0.62, p < 0.01) revealed that the reduction in SB interval reached significance at 2 and 5 μM (both p < 0.05). Average SB intervals returned to baseline following washes. The attenuation in SB intervals was replicated when cultures were exposed to the dihydropyridine-type L-type VGCC blocker Nifedipine www.nature.com/scientificreports www.nature.com/scientificreports/ ( Fig. 6A,B), strongly indicating that this observation was due to a reduction in L-type VGCC function. Both CACNA1C and CACNA1D genes, encoding the alpha 1 subunit of Cav1.2 and Cav1.3 respectively are expressed in neurons of these 60 DPP MEA cultures, however as Nifedipine has a 3-4 higher affinity to Cav1.2 than Cav1.3, inhibition of the former may be sufficient to alter network behaviour on MEAs.
To determine whether the effect on SB pattern was specific to the action of L-type VGCCs or Ca v channel function more generally, the Ca v P/Q type blocker ω-agotoxin TK and the Ca v T-type inhibitor ML218 were also applied acutely to cultures. Exposure of the cultures to both inhibitors had no effect on the basal firing rate and importantly on the SB interval (Fig. 6C,D). We did observe an apparent change in the shape profile of SBs after inhibiting T and P/Q type channels, however this did not translate to a significant reduction in SB firing rates (Fig. 6Db) at these concentrations. These results, therefore, strongly suggest that the regulation of synchronised network activity observed in these hiPSC derived neuronal cultures is due to the action of L-type VGCCs.

Discussion
We have established conditions to monitor the development of human neuronal networks derived from hiPSCs as they self-assemble in vitro and transition from uncoordinated, spontaneously active cultures to complex oscillatory networks with close similarities to those previously reported with dissociated rodent primary neurons. We further show that this complex network behaviour can be targeted specifically by pharmacological agents without change to the underlying basal neuronal activity. To relate this process to those changes in network behaviour associated with mental health disorders, we demonstrated the specific effects of a pharmacological model of L-type voltage gated calcium channel hypo-function, a class of ion channels which have been strongly implicated as risk factors for psychiatric disorders. www.nature.com/scientificreports www.nature.com/scientificreports/ Although the presence of synchronised array-wide activity has been reported previously in hiPSC-derived neuronal cultures 10,27 , this previously reported behaviour did not progress further than the short SB firing that we observed at 40 DPP and did not evolve to the bi-stable period of oscillatory firing we report here. Indeed, this extended behaviour actually compares favourably to what has previously been described by several groups using similar MEA systems to study networks in dissociated rodent primary neurons 3,5,28,29 . In particular, the extended MAP/LAP oscillations we observed in hiPSC-derived neurons at 50 DPP matches to those reported in rodent studies at ~20 DPP 7,29 .
Pharmacological profiling of our cultures showed that glutamate/GABA signalling is the primary driver of the complex network behaviour observed. Whereas at initial stages spontaneous neuronal activity was insensitive to AMPA or NMDA receptor inhibition, we observed that synchronised network activity in these cultures was completely eliminated by inhibition of NMDA and AMPA receptor signalling. It has previously been suggested that blockade of NMDA but not AMPA receptors abolishes similar slow coordinated activity both in vivo and in in vitro rodent slices [30][31][32][33] . However, work with dissociated rodent neurons has also shown that the coordinated burst firing is indeed blocked by both AMPA and NMDA receptor inhibition 2,34,35 . It is therefore possible that the difference seen in dissociated neurons compared to 'intact' models is a function of the more random and heterogeneous nature of such cultures, which lack the intrinsic complexity, highly regulated developmental structure and region specific networking seen in slices and in vivo. Nonetheless, this adds to the utility of hiPSC-derived neuronal networks, as they can monitor both NMDA and AMPA receptor-based glutamate signalling.
Inhibition of GABA A receptor activity modifies the SB pattern in our mature networks, decreasing the interval between SB periods, a finding consistent with reports from both rodent dissociated primary cultures 3,7,35 and hiPSC-derived neurons 9,10 . This indicates that although glutamate drives the complex network behaviour, GABA shapes it. The regulation of coordinated network firing by GABAergic interneurons is well documented, ranging from coordinating oscillatory activity across brain regions 34,36,37 , synchronised networks within structures (e.g. the hippocampus 38,39 ) to regulation of small, localised networks 40,41 . Furthermore, GABAergic activity is thought to be highly important throughout development, where it has been implicated in the correct formation of networks and the regulation of synaptic plasticity [42][43][44] . Importantly, inappropriate GABAergic function leads to a non-physiological shift in the inhibitory/excitatory balance which is thought to underlie several neuropsychiatric disorders including schizophrenia, ASD and epilepsies.
By capturing the properties of glutamate/GABA interaction of neuronal networks, the MEA neuronal networks we report here offer a cell assay platform to study the effects of risk alleles associated with neuropsychiatric disorders. To pursue this further, we examined the effects of inhibition of the L-type VGCCs Ca v 1.2 and Ca v 1.3, the genes encoding the core subunits of which, CACNA1C and CACNA1D, are strongly associated with neuropsychiatric disorders. The L-type VGCC blockers of Diltiazem and Nifedipine caused a dose-dependent reduction of the SB interval with little change to the basal excitatory profile. This effect is specific to L-type channels and was not observed when T or P/Q channels were similarly blocked. www.nature.com/scientificreports www.nature.com/scientificreports/ Ca v 1.2 and Ca v 1.3 alpha 1 subunits are expressed throughout the mammalian brain and are involved in wide array of calcium regulatory mechanisms 45,46 . In particular, they are strongly implicated in NMDA-dependant LTP/LTD, via regulation of local intracellular calcium concentrations [47][48][49] . This regulation of synaptic plasticity is thought to underlie the alterations to hippocampal -dependent learning seen with decreased or absent expression of L-type channels 47,50,51 . Furthermore, similar rodent KO studies have implicated L-type VGCCs in hippocampal-independent mechanisms of fear learning 52-54 , regulation of axon growth 55 , trafficking of AMPA receptor subunits 56 and gene expression 57,58 .
A role for L-type VGCCs has also been described in a range of oscillatory activities in rodent neurons. Synchronised calcium transients have been shown to be mediated and controlled by L-type VGCCs in dissociated primary neurons 22 , while further studies have shown that these channels can regulate calcium oscillations in intact systems 14,15,23 . Importantly, it has also been shown that L-type VGCC currents can modulate the network response of neurons during physiological oscillatory behaviour in the hippocampus 16,17,59 and during modelled epileptiform activity 60,61 . Involvement of L-type VGCCs in these systems is likely to be down to their role in the regulation of the post-burst after-hyperpolarisation (AHP), a period of hyper-polarisation which terminates high frequency firing 62 , where they have been shown to be a key mediator of the size and duration of the AHP 62-64 .
Our study presents for the first time the regulation of network activity in hiPSC derived neurons by modulating the function of L-type VGCCs. Importantly, this has implications for disease-modelling as mutations in the genes coding for L-type VGCCs represent some of the most strongly associated genetic risk factors for the development of mental-health disorders. Specifically, mutations in CACNA1C are amongst the most consistently detected genetic risk factors in both bipolar disorder and schizophrenia. Moreover, a gain-of function mutation in the pore forming region of Ca v 1.2 is responsible for development of Timothy syndrome, a neurodevelopmental disorder characterised by cardiac arrhythmia, heart malformations and ASD 65 .
In our study, the action of Diltiazem and Nifedipine induced a similar effect on synchronised network behaviour (SB intervals) as that caused by inhibition of GABA A receptors. Is it possible that the two pharmacological pathways share an underlying mechanism? A potential scenario is that Diltiazem is acting via presynaptic L-type VGCCs on GABAergic interneurons nerve terminals. Indeed, L-type VGCCs have been shown be expressed in Figure 6. Effects on periodicity of synchronised burst firing in hiPSC derived neuronal cultures is specific to L-type and not T-type or P/Q type voltage-gated calcium channels. To determine whether the attenuation of intervals between synchronised bursts (SBs) was due to a generic reduction of general voltage gated calcium channel activity, neuronal cultures were exposed to the alternative L-type channel inhibitor Nifedipine, the T-type channel inhibitor ML218 and the P/Q-type blocker ω-agatoxin TK. www.nature.com/scientificreports www.nature.com/scientificreports/ presynaptic terminals in the hippocampus 66 and in interneurons 67 , while they have also been shown to interact with aspects of the exocytotic pathways and facilitate neurotransmitter release in certain neuronal populations 68,69 . Furthermore, calcium currents via L-type VGCCs have been shown to act through the MAPK ERK1/2 signalling pathway, which itself has been implicated in the regulation of vesicular exocytosis 70,71 . Importantly, this link has been observed directly as inhibition of the ERK1/2 pathway has been shown to increase neurotransmitter release via increased calcium influx via L-type VGCCs 72 . This therefore provides a potential mechanism by which L-type VGCCs may be acting to regulate the activity of GABAergic interneurons in these iPS cell derived cultures: synchronised network activity develops in cultures after maturation of functional synapses and is dependent on both AMPA and NMDA signalling; this activity is regulated by GABAergic innervation, inhibition of which attenuates the period between high frequency culture-wide firing; blocking L-type VGCCs mimics the effect of GABA A antagonism by reducing vesicular release at the interneuron-projection neuron synapse.
Our findings have shown that hiPSC derived neurons self-assemble in vitro to form complex, extended coordinated network firing. This has only previously been described in primary rodent cultures. This behaviour can be pharmacologically manipulated owing to the role of both intrinsic glutamatergic and GABAergic innervation and as such, serves as a physiologically relevant in vitro model of human neural development. Importantly, networks are sensitive to the specific attenuation of L-type VGCC function, thereby demonstrating the use of such models to monitor the function of disease relevant genetic hypo-function. The scalability of the developed MEA assays, together with increased knowledge about the genetics of mental health disorders and the increased availability of patient or engineered hiPSC lines should provide exciting opportunities for probing the function of human neural networks in disease relevant models and the potential development of novel therapeutics.

Methods
Cells culture and neuronal differentiation. All experiments were conducted using the IBJ4 cell line, a control iPS cell line derived from the BJ fibroblast cell line (ATCC; CRL-2522). iPS cells were differentiated into forebrain projection neurons using an adapted dual SMAD protocol 73 . Briefly, iPS cells were maintained on matrigel (Corning) and in mTeSR1 medium (StemCell Technologies). Prior to the start of differentiation, cells were passaged onto 12 well plated coated with growth factor reduced matrigel (Corning). Once cells were >80% confluent (D0), medium was switched to N2B27 (2/3 DMEM/F12; 1/3 Neurobasal; B27-RA; N2; 1xPSG; 0.1 mM β-mecaptoethanol) +100 nM SB and 100 nM LDN 193189. After 8 days SB and LDN were removed and cells were passaged onto fibronectin coated 12 well plates at a ratio of 2:3. After a further 8 days, cells were passaged onto PDL/laminin coated coverslips as single cells at a density of 400 cells/mm 2 . Media was switched to N2B27 with B27+ retinoic acid, and 10 μM DAPT for 7 days, after which medium was changed to astrocyte conditioned BrainPhys (Stem Cell Technologies 74 ), supplemented with 1x B27, 1x PSG, 20 μM ascorbic acid and 20 ng/mL BDNF, in which cells remained for the remainder of the cultures. For astrocyte conditioning of BrainPhys, 50 mL of medium was added to 60-80% confluent flasks of <P5 Primary normal human astrocytes (Lonza), which were maintained previously in astrocyte medium (DMEM/F12 + 10% FBS and 1x PSG). Medium was conditioned for 72 hours at 37 °C in 20% O 2 , after which conditioned medium was removed from flasks, sterilised using a 0.22 μm syringe PVDF filter and stored at at 4 °C for use with 7 days or −20 °C for longer storage. Conditioned medium was used 1:1 with fresh BrainPhys.
Single cell electrophysiology. Patch clamping was performed on neurons using an electrophysiology rig comprising an Olympus BX51 WI upright microscope with 10x and 40x water immersion objectives; a Q-imaging Rolera Bolt CMOS camera; Luigs and Neumann Junior manipulators; Axon Instruments HS-2 Unity gain headstages; Axon Instruments Multiclamp 700B amplifier and Digidata 1550B; and a PC running Multiclamp commander. Cells were patched using borosilicate filament glass forged into pipettes using a Flaming/Brown puller (Sutter Instruments), resulting in resistances of 4-6 MΩ. All cells were patched at room temperature (around 22 °C) in a basic physiological extracellular fluid (ECF) comprising 142 mM NaCl, 2.5 mM KCL, 2 mM CaCl 2 , 1 mM MgCl 2 , 10 mM HEPES buffer and 30 mM D-glucose; pH adjusted to 7.4 with 4 M NaOH. The intracellular solution for pipettes consisted of 142 mM potassium gluconate, 1 mM CaCl 2 , 2 mM MgCl 2 , 10 mM HEPES and 10 mM EGTA; adjusted to pH 7.4 with KOH. The osmolality of internal solutions was adjusted to 290 mOsm, using a Vapour pressure osmometer (ELITech). After braking in to cells, resting membrane potentials were determined by taking the mean voltage from 2 minutes of gap free recording (I = 0). Occurrence of induced action potentials (iAPs) was determined by holding cells at around −70 mV and injecting positive current steps (Δ5-10 pA; 1 second) until APs were seen. If no APs were seen when the membrane potential reached −10 mV, it was assumed none would be seen at all. iAP categorisation of cells was based on visual assessment of traces. To be deemed a full AP, the overshoot had to be greater than 0 mV otherwise they were deemed as 'attempted' . AP trains were determined as two or more full APs seen within the stimulus period (Fig. S1). Spontaneous action potentials were observed with gap free recording with cells at resting membrane potential (I = 0). Compensation for fast and slow capacitive transients was applied after seal formation and series resistance (bridge balance) values were corrected after break in and were monitored throughout experiments. Recordings were taken using Clampex 10 (Molecular Devices, USA) and were subsequently analysed using Clampfit 10 (Molecular Devices, USA).
Multi electrode array culturing and recording. 60MEA200/30iR-Ti-gr MEAs were used throughout the study (Multi channel systems; MCS); planar MEAs with 60 titanium nitrate electrodes (59 + 1 internal reference) embedded within a silicon nitrate substrate. Each electrode had a diameter of 30 μm, with 200 μm spacing between electrodes. MEAs were recorded using a MEA2100-HS2x60 headstage amplifier, attached into MSC-IFB-3.0 analogue/digital interface board. Cultures were recorded using the MCRack data acquisition software running on a high-performance PC. The recording head stage was isolated within a custom-made faraday cage.
www.nature.com/scientificreports www.nature.com/scientificreports/ Coverslips of cells earmarked for arrays were re-plated onto MEAs at D40-45. Clean MEAs were stored at 4 °C in the dark, with the cultured area submerged in distilled water. Before cell plating, MEAs were first pre-treated with 1 ml FBS for 1 hour at 4 °C in the dark, followed by washes with distilled water. Culture surfaces were then treated with 0.01% PEI (Sigma) and incubated for 1 hour at 37 °C, were then washed with distilled water and left to dry completely in a sterile cell culture hood 1 hour prior to re-plating. Neurons were dissociated from coverslips with Accutase (Thermo Fisher) and plated as high density drop cultures (50,000 cells/20 μL, which equated to around 1800 cells/mm 2 ) containing 10 μg/mL laminin. After 1 hour, 1 mL of conditioned medium was then carefully flooded into each MEA and after a further 24 hours incubation, 1 ml of fresh medium was added to each array. MEA cultures were maintained in 1:1 fresh/ACM BrainPhys and in an incubator with a 2% O 2 environment. A low oxygen incubator environment has been shown to increase the survival of neural precursors 75 as well as promoting neuronal maturation 18 . All recordings were made at 37 °C.
All drugs used in this study were from Tocris (Bio-Techne Ltd.). For MEA drug experiments, stock concentrations were diluted as required in 1 mL 1:1 Brainphys:Brainphys ACM. Medium containing diluted drugs were added to MEAs and subsequently incubated for 10 mins at 37 °C in a standard 5% CO 2 incubator environment. After recordings, medium was removed from MEAs and cultures were washed 3 times with PBS. 1 ml of fresh medium was then added and cultures were incubated for 5 minutes as before. This medium was then removed and replaced with fresh medium for 'Wash' recordings.
Raw MEA data was recorded at a sample rate of 25000 Hz and filtered online with a 200 Hz high pass and a 5000 Hz low pass filter (both 2 nd order butterworth). Recorded data was converted to ascii files for offline analysis using the MC Data tool (MCS). Offline analysis was achieved with custom scripts written in Matlab. Briefly, spikes were detected from filtered data using an automatic threshold-based method set at −5.5 × σ , where σ is an estimate of the noise of each electrode based upon the median absolute deviation (MAD 20 ). Spike timestamps were analysed to provide statistics on the general excitability of cultures. Network activity was analysed by creating array-wide spike detection rate (ASDR) plots with a bin width of 200 ms. Synchronised bursts (SBs) were detected from ASDR plots by a 4-step process: 1; A 'start threshold' was determined as ~10% of the maximum ASDR -this was varied to set the threshold just above baseline activity; 2. Data passing the start threshold was then required to pass a 'confirmation threshold' within 1 second -this was set as 70% of the maximum ASDR; 3. The end of SBs was determined as there being 3 seconds of <start threshold activity; 4. SBs were finalised by confirming that they lasted at least 2 seconds. Average results for every measure were calculated as medians for each culture and data from single electrodes was excluded if it contained less than 5% of the spikes detected in the most active electrode (see Supplementary Information for full description of the analysis parameters).
Statistical analysis. All descriptive and comparative statistics were completed with Prism 6 (Graphpad. To determine the route of analysis (parametric or non-parametric), all data was first processed using histograms, q-q plots and normality tests (kurtosis and skew tests) in R. Where groups of data from the same experiment, presented with contrasting distributions, parametric tests were used as they are, in general, better equipped to cope with non-gaussian distributed data. Unless otherwise stated, all comparisons between pre-and post-treated cultures were done with paired t-tests and tests statistics, degrees of freedom and two-tailed p values are all reported. All summary plots of data show means + standard deviation. The numbers of arrays used for each experiment vary and are reported in individual figures. All results presented here were collected from at least three independent differentiations.