The Role of Parvalbumin-positive Interneurons in Auditory Steady-State Response Deficits in Schizophrenia

Despite an increasing body of evidence demonstrating subcellular alterations in parvalbumin-positive (PV+) interneurons in schizophrenia, their functional consequences remain elusive. Since PV+ interneurons are involved in the generation of fast cortical rhythms, these changes have been hypothesized to contribute to well-established alterations of beta and gamma range oscillations in patients suffering from schizophrenia. However, the precise role of these alterations and the role of different subtypes of PV+ interneurons is still unclear. Here we used a computational model of auditory steady-state response (ASSR) deficits in schizophrenia. We investigated the differential effects of decelerated synaptic dynamics, caused by subcellular alterations at two subtypes of PV+ interneurons: basket cells and chandelier cells. Our simulations suggest that subcellular alterations at basket cell synapses rather than chandelier cell synapses are the main contributor to these deficits. Particularly, basket cells might serve as target for innovative therapeutic interventions aiming at reversing the oscillatory deficits.

Fast-spiking (FS), PV + , γ-amino-butyric acidergic (GABAergic) interneurons seem to be a major contributor to gamma oscillations 18,19 . However, this class of inhibitory interneurons can be further divided into at least two subgroups: basket cells (BCs) and chandelier cells (ChCs) 20,21 . The main difference between basket and chandelier cells is their axonal target. While basket cells mainly target the perisomatic region of excitatory pyramidal cells and other inhibitory interneurons, chandelier cells exclusively target the axon hillock of pyramidal cells (PCs) 20 . Furthermore, although their electrophysiological properties are similar subtle differences, such as different afterhyperpolarization amplitudes and firing frequencies, between the two subclasses exist 22 . Lastly, while BC synapses display fast and strong inhibition, ChC synapses have been shown to have a depolarizing effect in vitro [23][24][25][26] . However, whether this depolarization also exists in vivo and whether it has an excitatory or rather a predominantly shunting effect remains elusive. Moreover, it is not clear so far whether these differences also imply a different involvement in the generation of gamma oscillations.
Several alterations of PV + neurons have been found in schizophrenia. While the density of PV + neurons does not seem to be altered in subjects with schizophrenia 27,28 , the levels of mRNA for the calcium-binding protein PV are markedly reduced 29 . Furthermore, 50% of PV + neurons in dorsolateral prefrontal cortex lack detectable levels of the 67 kDa isoform of glutamate decarboxylase (GAD 67 ), an enzyme that catalyzes the decarboxylation of glutamate to GABA 29 . It has been suggested that the reduced expression of GAD67 mRNA, the gene encoding the GAD 67 enzyme, likely implies a reduction in GABA synthesis in cortical GABA neurons, which in turn would lead to smaller amplitude IPSCs at the postsynaptic site 30 . Furthermore, a reduction in the plasma membrane GABA transporter GAT1, which is a major contributor to the specificity of synapses by preventing spillover to neighbouring synapses 31,32 , has been found in PV + interneurons in patients with schizophrenia 27 . A reduction in GAT1 leads to a prolongation of IPSC durations 31,32 . Interestingly, in patients with schizophrenia, immunoreactivity for the GABA-A α 2 subunit is increased 33 in the axon initial segment of pyramidal cells, the exclusive target of chandelier cell cartridges. GABA receptors with an α 2 subunit, compared to receptors with an α 1 subunit, also show an increase in IPSC duration in tranformed human embryonic cells 34 . In summary, the reduced expression of GAD67, the redcution of GAT1 and increased α 2 subunit expression in SCZ patients have two main effects: a reduction of inhibitory strength of PV + neurons (both BCs and ChCs) and an increased IPSC duration at PV + synapses (also for both BCs and ChCs, however, considerably longer for ChCs). These cellular differences in PV + neurons are likely to have an influence on oscillatory activity. Since gamma oscillation generation and maintenance crucially depends on fast-spiking PV + interneurons, changes in peak amplitude and duration of IPSCs might particularly affect oscillations in this frequency band.
In addition to evidence suggesting an immediately compromised GABAergic system, the finding that N-methyl-D-aspartate receptor (NMDAR) antagonists (such as phencyclidine and ketamine) produce symptoms very similar to key clinical features of schizophrenia 35 , sparked the idea that the underlying cause of the disturbance of the GABAergic cicuitry might be a secondary effect of an NMDAR hypofunction at GABAergic interneurons. While studies exploring the idea that the hypofunction was due to less receptors have produced mixed results 36 , convergent lines of evidence underpin the hypothesis that NMDARs are dysfunctional in schizophrenia patients 37 . Direct evidence in favour of this hypothesis are for example, changes in NMDAR-associated protein levels at the postsynaptic site 38 , a reduction in NMDAR-mediated signalling following neuregulin 1 activation of ErbB4 receptors 39 , lower levels of glutathione (a modulator at the redox-sensitive site of NMDARs) 40 and a reduction of kynurenine 3-monooxygenase which might increase kynurenic acid (an NMDAR antagonist) levels 41 .
Earlier modelling efforts by Vierling-Claassen et al. 8,42 , proposed that increased IPSC decay times at chandelier cell synapses might be sufficient to explain the gamma and beta range ASSR deficits in schizophrenia patients. However, while demonstrating that, in principle, increased IPSC decay times at chandelier cell synapses could result in experimentally observed ASSR deficits for patients with SCZ, Vierling-Claassen et al. did not consider the relatively small number of chandelier cells compared to basket cells 43 . Furthermore, they did not take into account potentially increased IPSC decay times at basket cell synapses, changes in excitability of interneurons due to NMDAR hypofunction and the putative excitatory effect of ChC-mediated neurotransmission. Therefore, we extended their computational model to investigate the differential role of prolonged IPSC decay time at basket and chandelier cell synapses, at realistic ratios of chandelier to basket cells, in gamma and beta range ASSRs. We found that, while prolonged IPSC decay times at chandelier cells could explain experimental results when chandelier cells form 50% of inhibitory neurons in the network, they failed to do so at lower, more realistic percentages such as 10% or even 5%. In contrast, for physiologically plausible chandelier cell percentages of 10% or 5%, prolonged IPSC decay times at basket cell synapses were able to reproduce experimental findings. Increased IPSC decay times at synapses of both interneuron subtypes resulted in qualitatively very similar ASSRs as for changes at basket cell synapses only, suggesting that chandelier cell alterations might not be crucially involved in gamma and beta range ASSR deficits in schizophrenia. We then wondered whether a potentially excitatory effect of chandelier cells onto pyramidal cells could explain the deficits but found that, at low chandelier cell percentages, it had a negligible effect on the overall network in general and on gamma and beta range ASSRs specifically. Lastly, decreased excitability of chandelier cells at low chandelier cell percentages, as a consequence of reduced, tonic NMDAergic excitation, did not produce ASSR deficits consistent with experimental observations, in fact, it hardly changed the network behaviour at all.

Results
In order to explore the roles chandelier and basket cells play in the emergence of oscillatory entrainment deficits in schizophrenia patients, we simulated the effects of cellular and molecular abnormalities of both cell types found in schizophrenia, focusing on GABAergic synapses. We started with an investigation of abnormalities of chandelier cells in isolation, followed by abnormalities of basket cells in isolation and an exploration of interactions of abnormalities in both cell types, followed by an exploration of the potentially excitatory ChC-PC coupling. We concluded with an exploration of NMDAR hypofunction at ChC interneurons and its interplay whith changes at GABAergic synapses in determining the oscillatory behaviour of the network.
We simulated the effect of GABAergic abnormalities for four different conditions with a BC/ChC ratio of 50%, 25%, 10% and 5%, respectively. As results at high ChC percentages (50 and 25%) were similar and results at low ChC percentages (10 and 5%) were also similar, we only present results for 50% and 10%. Based on the investigation of the GABAergic alterations, we focused on low ChC percentages of 10% when exploring the excitatory ChC-PC coupling and the NMDAR hypofunction.
Alterations of chandelier cells. As explained earlier, the plasma membrane GABA transporter GAT1 is reduced in chandelier cell synapses onto the axon initial segment of cortical pyramidal neurons, resulting in a prolongation of IPSC decay times. Additionally, the α subunit composition of the GABA A receptors at the axon initial segment is skewed towards the α 2 subunit in comparison to receptors at the dendrites, which leads to an even more prominent prolongation of decay times. This increased time course of inhibition at chandelier cell synapses has been proposed to account for the characteristic oscillatory entrainment changes in the gamma and beta range by Vierling-Claassen et al. 8,42 . However, in their studies two important aspects of the local microcircuitry have not been accounted for: 1) the relatively small number of chandelier cells (compared to basket cells) 43 and 2) the specific connectivity scheme of chandelier cells, i.e. the fact that chandelier cells do not form connections with themselves 44 . We explored whether prolonged IPSC decay times at chandelier cell synapses (τ = 28 ms ChC , as in Vierling-Claassen et al. 8 ) could still explain schizophrenia-like entrainment deficits when these two aspects were included in the model.
At a ChC percentage of 50, the 40/40 measure was drastically reduced in the 'schizophrenia-like' network, reflecting a disturbance in the generation of gamma oscillations, which is in line with experimental evidence (see Table 1). However, the other entrainment measures did not match with experimental findings. The 30/30 measure was also strongly reduced, meaning that the disturbance in the generation of oscillations was not restricted to the gamma range at around 40 Hz. Furthermore, the 20/20 measure did not differ strongly between the 'schizophrenia-like' and the control network, while a difference was found in experimental studies (see Thune et al. 17 ). Lastly, the emergence of a large 20/40 measure, which has been observed in the experimental study of Vierling-Claassen et al. as well as in their modelling work 8 , could not be seen. As can be seen in Fig. 1, the 20/40 measure, although higher than in the control network, remained small (note the scale in the 20/40 plot of Fig. 1(b)).
At lower, more realistic, values of ChC percentage hardly any difference of the MEG signals and the five entrainment measures between control and 'schizophrenia-like' networks could be seen (see Fig. 1). Note that this did not change with further increases of the IPSC decay times (we tested τ = 35 ms chc and τ = 40 ms chc ; data not shown). This suggests that prolonged IPSC decay times at ChC synapses are not sufficient to induce schizophrenia-like alterations in oscillatory entrainment.
Alterations of basket cells. The plasma membrane GABA transporter GAT1 is also reduced in basket cell synapses onto dendrites in the perisomatic region of cortical pyramidal neurons, again resulting in a prolongation of IPSC decay times at these synapses 27 . However, opposed to chandelier cells, GABA A receptors at basket cell synapses predominantly contain the α 1 subunit 33 . This means that the IPSC decay time increases at BC-PC synapses probably also exist in schizophrenia patients but the increase is likely to be smaller than that at the ChC-PC synapses we examined before. Therefore, we initially chose an increased IPSC decay time of τ = 20 ms BC . At high ChC percentages, this alterations had hardly any effect on the oscillatory behaviour of the network in response to any of the three driving frequencies tested (see Fig. 2(a) top row and Fig. 2(b)).
At low, more realistic ChC percentages this alteration led to changes of oscillatory entrainment more in line with experimental evidence (see Fig. 2). Specifically, the gamma range entrainment was signifcantly disturbed, reflected by the strong decrease of the 40/40 measure. Furthermore, a characteristic increase of the 20/20 measure was also clearly visible. The changes were restricted to the gamma and beta band leaving the respective border intact, i.e. the 30/30 measure was not significantly altered. However, the emergence of a strong 20/40 measure, seen by Vierling-Claassen et al. 8 , both in experiment and model, was not visible. Figure 2(b) (middle panel) does show an increase of the 20/40 measure compared to control, however it remains small. Therefore, we increased the IPSC decay time at BC-PC synapses further to τ = 28 ms BC . For high ChC percentages, the further increase of τ BC led to a strong decrease of the 40/40 measure but only had minor effects on all other entrainment measures (lower row of (a) and the third panel of (b) of Fig. 3).
As shown in Fig. 3 (lower row of (a) and the third panel of (b)), in the case of low ChC percentages, this led to a prominent increase of the 20/40 measure, similar to what was found in experiments 8 . Figure 4(a) shows the evolution of this 20 Hz peak for 40 Hz drive over three different values of τ = 20 ms BC (black), 24 ms (red) and 28 ms (green). While the PCs are still able to recover from inhibition within every gamma cycle for τ = 20 ms BC , for larger values of τ BC , more and more PCs are suppressed during every other gamma cycle, giving rise to a prominent 20 Hz peak. Interestingly, however, the longer IPSC decay time of τ = 28 ms BC did not only effect the 20/40 measure but also led to a strong decrease of the 30/30 measure, which then became markedly different from the 30/30 measure of the control network (see Figs. 3 and 4(b)). This is in disagreement with experimental findings.

Alterations of chandelier and basket cells. The results presented in the previous sections suggest that
ChCs (at realistic numbers) are unable to produce 'schizophrenia-like' gamma and beta entrainment deficits in isolation and that BCs are the major contributor to these. However, alterations at ChCs might nevertheless modulate the deficits produced by BC alterations. Therefore, we next looked at the interactions of simultaneous alterations at BCs and ChCs.
First, we simulated the combined effect of a moderate increase of IPSC decay times at both ChC-PC and BC-PC synapses, with a stronger increase at ChC-PC synapses (i.e. τ = 28 ms  Table 1. Summary of ASSR deficits in patients with schizophrenia in the three studies considered here (↓: significantly lower in patients, ↑: significantly higher in patients, -: no significant difference between controls and patients). www.nature.com/scientificreports www.nature.com/scientificreports/ Comparison of the 5 entrainment measures between control (black) and 'schizophrenia-like' (red) for high and low ChC percentages. Note that 'schizophrenia-like' here means that IPSC decay times at basket cell synapses were increased from 8 ms to 20 ms. All traces here depict the mean signal averaged over 20 trials with different realisations of the Poisson process underlying the generation of noise. Importantly, averaging was performed in the time domain prior to a transformation into the frequency domain via Fourier transform, which is also the standard procedure in experiments. www.nature.com/scientificreports www.nature.com/scientificreports/ condition, in disagreement with experimental findings 7,8 . For biologically realistic low ChC percentages we also found a strong reduction of the 40/40 measure (Fig. 5(a)), and a small reduction of the 40/20 measure (Fig. 5(b)). Interestingly, in this case the 20/20 measure was increased and, although still low, the 20/40 measure was also increased, while the 30/30 measure was only marginally different. Overall, the changes in spectral composition were consistent with experimental findings suggesting a general gamma deficit and beta increase in schizophrenia 8 , similar to the case of τ = 20 ms BC . Lastly, we increased the IPSC times at BC-PC synapses further to 28 ms and, simultaneously, also increased IPSC times at ChC-PC synapses to 40 ms. In the low ChC percentage condition, this further reduced the 40/40 amplitude while increasing the 20/20 measure. Additionally, a 20/40 peak emerged. However, there was also a strong reduction of the 30/30 measure (Fig. 6). Overall, these results here are very similar to the results from the simulations where only IPSC times at BC-PC synapses were modified, suggesting a minor role of ChCs in the generation of ASSR deficits in schizophrenia.

Excitatory effects of chandelier cells. Several in vitro studies have demonstrated an excitatory effect
of ChC synapses onto PCs, most likely due to differences in the GABAergic reversal potential at the axon initial segment of PCs compared to somatic or dendritic segments. Therefore, we next wondered whether an excitatory effect of ChC-PC neurotransmission together with the increased IPSC decay times at these synapses might explain the experimentally observed changes in gamma and beta range ASSRs. We simulated a newly configured network model, where the ChC-PC synapses had a positive instead of a negative weight factor α (see Methods for details on the synapse models in the network). First, we confirmed that this new configuration qualitatively reproduced the ASSR behaviour of the control network (see Fig. 7(a)). Next, we increased the IPSC decay times at ChC-PC synapses from 8 ms to 28 ms, as before. As can be seen in Fig. 7(b), excitatory ChC-PC coupling together www.nature.com/scientificreports www.nature.com/scientificreports/ www.nature.com/scientificreports www.nature.com/scientificreports/ www.nature.com/scientificreports www.nature.com/scientificreports/ with prolonged IPSC times did not lead to substantial changes in any of the 5 entrainment measures. Last, we tested whether excitatory ChC-PC synapses would modulate the ASSR changes caused by prolonged IPSC decay times at BC-PC synapses we had seen before. Therefore, we increased the IPSC decay times at BC-PC synapses from 8 ms to 28 ms in the newly configured network model. Figure 7(c) shows that the excitatory ChC-PC coupling also did not substantially alter the ASSR changes caused by BC-PC synaptic changes. NMDAR hypofunction. Lastly, we explored the consequences of an NMDAR hypofunction at ChC interneurons by first investigating the effect of NMDAR hypoactivation in isolation and subsequently, together with increased IPSC at ChC-PC and BC-PC synapses, respectively. Since our network does not explicitly model www.nature.com/scientificreports www.nature.com/scientificreports/ NMDA receptors, we modelled NMDA hypoactivation as a reduction of ChC excitability (see Methods for implementation details) similar to previous modelling efforts 45 . We found that NMDAR hypofunction at ChC interneurons did not substantially alter the networks response to the oscillatory drive. Figure 8 shows that none of the 5 entrainment measures showed a large difference between the control and NMDAR hypofunction condition despite the decrease in NMDAR-mediated excitation modelled being relatively large. Furthermore, NMDAR hypofucntion at ChCs combined with an increased IPSC decay time did not change the effect increased IPSC times showed in isolation, neither for IPSC increases at ChCs (compare Figs. 9(a) and 1(b)) nor for IPSC changes at BCs (compare Figs. 9(b) and 3(b)).

Discussion
Here, we explored the contribution of slower IPSC dynamics at ChC-PC synapses to gamma and beta range ASSR deficits in patients with schizophrenia using a model cortical network. Generally, at realistic, low ratios of ChCs to BCs, increased IPSC decay times at ChC synapses alone hardly had any influence on ASSRs. At low ChC/BC ratios, increased decay times at BC-PC synapses, on the other hand, strongly altered the spectral composition of ASSRs, replicating several characteristics found experimentally. For a moderate increase (i.e. τ = 20 ms BC ), we found a strong decrease of the 40/40 measure, an increase of the 20/20 measure and an increase of the 20/40 Figure 8. Comparison of the 5 entrainment measures between control (black) and 'schizophrenia-like' (red) for low ChC percentages. Note that 'schizophrenia-like' here means that NMDAR hypofunction was modelled as a strong reduction of excitability of chandelier cells. All traces here depict the mean signal averaged over 20 trials with different realisations of the Poisson process underlying the generation of noise. Importantly, averaging was performed in the time domain prior to a transformation into the frequency domain via Fourier transform, which is also the standard procedure in experiments. Figure 9. Comparison of the 5 entrainment measures between control (black) and 'schizophrenia-like' (red) for low ChC percentages. Note that 'schizophrenia-like' here means that NMDAR hypofunction was modelled as a strong reduction of excitability of chandelier cells and additionally the IPSC decay times at (a) chandelier cell synapses and at (b) basket cell synapses were increased from 8 ms to 28 ms, respectively. All traces here depict the mean signal averaged over 20 trials with different realisations of the Poisson process underlying the generation of noise. Importantly, averaging was performed in the time domain prior to a transformation into the frequency domain via Fourier transform, which is also the standard procedure in experiments. measure, although it remained low (Fig. 2). This suggests that this change shifts the preferred frequency of the network from the gamma range more towards the beta range and gives an explanation of the reduced gamma seen across many experiments. Furthermore, the 30/30 measure remained virtually unaffected. However, it did not explain the strong 20/40 component seen by Vierling-Claassen et al. 8 . For a stronger increase (i.e. τ = 28 ms BC ), we found an even stronger decrease of the 40/40 measure, a further increase of the 20/20 measure and a much stronger increase of the 20/40 measure (see Fig. 3). However, we also saw an increase of the 40/20 measure, contradictory to the gamma to beta shift hypothesis. Furthermore, the 30/30 measure becomes greatly reduced ( Fig. 3(b)), which is not seen experimentally. Looking at the power spectra for 40 Hz and 30 Hz drive for increasing τ BC , we see that with longer decay time the 20/40 measure gets stronger but at the same time the change to the network increasingly affects the 30/30 measure (Fig. 4). Lastly, when looking at the effect of increased decay times of ChC-PC synapses in a network where also BC-PC synapse dynamics are slowly decaying, we see that the additional contribution of ChC-mediated effects is very weak for different combinations of decay times τ BC and τ ChC (see Figs. 5 and 6). One interesting effect, however, was the smaller 20/40 measure here compared to the condition where only BC-PC synapses had longer decay times (compare Figs. 3 and 6). Additional to the more classical view that GABAergic chandelier cells provide inhibitory input to pyramidal cells, we have also explored the potential effect of an excitatory ChC-PC coupling, which has been proposed recently [23][24][25][26] . However, at low ChC percentages, we did not see a substantial effect of this excitatory coupling, neither did it change the response to oscillatory drive in the control network, i.e. without alterations of the dynamics of GABAergic synapses, nor did it modulate the effect of changed GABAergic dynamics (see Fig. 7). Since changes to the GABAergic system in schizophrenia have been hypothesized to be a secondary effect due to NMDAR hypofunction at inhibitory interneurons, we also explored the potential contribution of reduced slow NMDA currents at ChC cells. Again, we saw that strongly reduced NMDAR function did hardly affect the ASSRs of the network at low ChC percentages.
Overall, our findings suggest that at realistic, low ratios of ChCs to BCs in a microcircuit, gamma and beta range changes to the spectral composition of ASSRs as seen in SCZ patients, can be well explained by an increase of GABAergic decay times at BC-PC synapses. While postmortem studies indicate that the reduction in GAT1 in SCZ patients preferentially affects ChCs 46 , this might simply reflect easier detection due to the specifc morphology of ChC synaptic cartridges compared to BC synapses and not a general difference in GAT1 reduction between the two subtypes. Nevertheless, prolonged IPSC decay times as a result of reduced GAT1 levels requires a high synapse density 31 and, thus, it is not immediately clear if a GAT1 reduction would be sufficient to cause a prolonged IPSC time course at BC synapses. Interestingly, however, repeated activation together with a reduction in GAT1 also resulted in an increase in IPSC decay times in the study of Overstreet et al. 31 . Therefore, increased IPSC times might also be present at BC synapses, especially in a gamma and beta entrainment setting where the network is driven with a relatively high frequency.
Our results demonstrate that prolonged IPSC decay times at BC-PC synapses robustly lead to a weakening of gamma synchronization, especially during gamma drive but also during beta drive. This is consistent with a large body of experimental evidence reporting reduced gamma band power during ASSR paradigms 7-17 but also during various sensory (e.g. a visual Gestalt task 5 ), and cognitive tasks (e.g. a working memory task 47 ). We also found an increase of beta synchonization during beta drive, and, although in most conditions still very small, increased beta synchronization during gamma drive. This reproduces the findings of Vierling-Claassen et al. 8 that slow dynamics at FS PV + interneurons might not only desynchronize gamma oscillations but also shift the preferred frequency of cortical networks from gamma to the beta band. However, our results further demonstrate that, most likely, BCs are the main interneuron subtype causing this shift and not ChCs, as hypothesized by Vierling-Claassen et al. 8 . Furthermore, we found that while longer decay times at BC-PC synapses led to a stronger and more robust peak of the 20/40 measure, they also introduced a stronger reduction of the 30/30 measure. Given that studies have consistently reported that the 30/30 measure remains unaffected in SCZ patients 7,8,17 , together with the fact that the strong 20/40 component has not been reported by experimental studies except for Vierling-Claassen et al. 8 , our results suggest that a strong 20/40 component might not be a core ASSR deficit in patients with schizophrenia. This is further underpinned by the finding that alterations at ChCs together with BC alterations, led to a reduction of the 20/40 measure compared to the condition where only BC-PC had prolonged decay times.
While the results presented here suggest a modulatory influence of ChC alterations on gamma and beta oscillations at most, we have to note that the employed model has substantial limitations. Our model, for example, does not explore the possibility that, due to the location directly at the axon initial segment, ChCs might exert significantly stronger control over PC firing, which in turn might counterbalance their small numbers. While we have investigated the potential consequences of NMDAR hypofunction affecting ChCs, we haven't explored NMDAR hypofunction in BCs because we were primarily interested to understand the contribution of ChCs within this model of prolonged IPSC times. Experimental and computational studies have linked NMDAR hypofucntion to gamma deficits 48,49 and further research in this area is definitely warranted. However, we must admit that our model is not very well suited for this task because NMDARs are not explicitly modelled. Finally, the theta neuron model is a simple and abstract neuron model and does not capture electrophysiological differences between BCs and ChCs found experimentally 22 . Interestingly, genetic variants of risk genes for schizophrenia include a variety of genes coding for ion channels or ion transporters, which have been shown to alter cell excitability 50,51 and network oscillations 52 . Such variants might potentially amplify electrophysiological differences between BCs and ChCs and thereby change their relative roles in the generation of gamma/beta oscillations.
In summary, our findings suggest that ChCs are likely not strongly involved in gamma/beta range ASSR deficits in schizophrenia, and, therefore future modelling efforts should focus on alterations at BC-PC synapses. Furthermore, ChC-PC synapses might not represent a valuable target for therapeutic interventions aiming at restoring the gamma ASSR in patients with schizophrenia. (2019) 9:18525 | https://doi.org/10.1038/s41598-019-53682-5 www.nature.com/scientificreports www.nature.com/scientificreports/

Methods
The model. The model proposed here is based on a recent reimplementation 53 of the simple model presented by Vierling-Claassen et al. 8 . However, we extended the model to consist of two distinct inhibitory populations (instead of only one), one representing basket cells and one representing chandelier cells (see Fig. 10).
Single cell model. Individual cells are modeled as theta neurons (for a detailed discussion of this neuron model see e.g. Boergers et al. 54 ).
The kth neuron in a network is described by a single variable θ k , which can be regarded as the neuron state, subject to the following dynamics where b is an externally applied current, S is the total synaptic input to the cell and N t ( ) is a time-varying noise input. Total synaptic input to a cell in a network is calculated as where n is the number of presynaptic neurons, α j controls excitation and inhibition, i.e. is +1 for excitatory synapses and −1 for inhibitory ones, g jk is the synaptic strength from cell j to cell k and s jk is the synaptic gating variable from cell j to cell k. Synaptic gating variables are subject to the following dynamics where τ j is the synaptic decay time, τ R the synaptic rise time and η is a scaling parameter. The network receives excitatory drive input at click train frequency from a single pacemaker cell. Additionally, Poissonian noise input is also given to all cells in the network. A noise spike at time t n elicits the following excitatory postsynaptic potential ('EPSP') where ⋅ A g gmax is the strength of the noise, τ exc is the synaptic decay time, τ R the synaptic rise time, and H the Heaviside function. An overview of all the model parameters used in this study can be found in Table 2.
Network. We combined 80 excitatory, pyramidal cells together with 40 inhibitory cells of the two different inhibitory subtypes (basket and chandelier cells) into a network model, following the earlier work of Vierling-Claassen et al. 8 . The percentages of the two different subtypes varied.
The connectivity of the network is summarized in Fig. 10. Pyramidal cells make connections to other pyramidal cells, as well as to both types of inhibitory interneurons. Basket cells inhibit pyramidal cells, themselves and chandelier cells. Chandelier cells only inhibit pyramidal cells (see e.g. Tremblay et al. 21 ). The connectivity between www.nature.com/scientificreports www.nature.com/scientificreports/ any two populations is all-to-all. All populations also have two sources of input, the oscillatory drive input and a background noise input. The drive input periodically sends spikes with a given frequency to all three populations. In order to generate the spikes a drive cell is implemented (also a theta neuron), which receives an applied current so that its spike frequency matches the driving frequency. This drive cell is then connected to all cells in the network. The background noise input sends noise spikes at times drawn from a Poisson distribution.
Entrainment measures. In order to evaluate the oscillatory entrainment we record the spiking of the three populations as well as a simulated 'EEG/MEG' signal, calculated as the sum of all excitatory synaptic variables over all pyramidal cells 8 Table 1, a decrease in the 40/40 measure is a robust finding in all three considered studies. Furthermore, the 30/30 measure is not affected in all three studies. Beyond the 40/40 measure, only the study by Vierling-Claassen et al. 8 found significant differences for other measures, namely, a decrease in the 40/20 measure, underpinning the deficit in gamma generation, and an increase in the 20/20 and the 20/40 measures, suggesting a transfer of power from the gamma to the beta band in patients with schizophrenia. It is important to note here, that both the Kwon et al. 7 and the Krishnan et al. 9 studies, used EEG whereas the Vierling-Claassen et al. 8 study used MEG and that it has been argued that MEG might have a higher sensitivity, which could explain the absence of findings other than for the 40/40 measure in the first two studies. Overall, given the experimental evidence discussed above, we think that the five measures chosen provide a good coverage of the frequency bands that are important to understand gamma and beta entrainment deficits in schizophrenia. Implementation. The model was implemented using Python 2.7.9 and numpy 1.9.3. Analysis and visualization of the model output was also done in Python using the numpy and matplotlib packages (matplotlib 1.4.3).
All differential equations were solved using a simple forward Euler scheme. A single simulation simulated a 500 ms trial and the time step was chosen such that this results in = 2 8192 13 data points. However, the main results were unaffected by using a smaller time step.
All code (for simulations, analysis and visualization) has been made publicly available in the following github repository: https://github.com/ChristophMetzner/Chandelier-Basket-Model.
Exploration of circuit abnormalities. We explored the effects of natural differences between BCs and ChCs together with their differential schizophrenia-associated abnormalities on oscillatory entrainment in the gamma and beta frequency range. Specifically, we investigated the ratio of BCs versus ChCs, a prolonged IPSC decay times of inhibitory synapses of both interneuron subtypes (because of a reduced expression of GAT1 for both BCs and ChCs and an increased expression of α 2 subunits for ChCs), the putative excitatory effect of ChC-PC synapses and reduced levels of ChC excitability due to NMDAR hypoactivation. Table 3 details how these changes were modelled.  www.nature.com/scientificreports www.nature.com/scientificreports/ Because of the Poissonian background noise input, simulation results varied from trial to trial. To capture the robustness of findings we always performed 20 simulation trials, each with a different realisation of the underlying Poissonian noise process. When analysing the effect of circuit abnormalities, we averaged the 20 trials for each condition in time to get an average simulated MEG signal and calculated a power spectral density of this signal.