Cellular and synaptic phenotypes lead to disrupted information processing in Fmr1-KO mouse layer 4 barrel cortex

Domanski, Aleksander P. F.; Booker, Sam A.; Wyllie, David J. A.; Isaac, John T. R.; Kind, Peter C.

doi:10.1038/s41467-019-12736-y

Download PDF

Article
Open access
Published: 23 October 2019

Cellular and synaptic phenotypes lead to disrupted information processing in Fmr1-KO mouse layer 4 barrel cortex

Aleksander P. F. Domanski ORCID: orcid.org/0000-0001-9582-5192^1,2,3,4,
Sam A. Booker^2,3,5,
David J. A. Wyllie ORCID: orcid.org/0000-0002-4957-6049^2,3,5,6,
John T. R. Isaac⁴^nAff7 &
…
Peter C. Kind^2,3,5,6

Nature Communications volume 10, Article number: 4814 (2019) Cite this article

5330 Accesses
24 Citations
3 Altmetric
Metrics details

Subjects

Abstract

Sensory hypersensitivity is a common and debilitating feature of neurodevelopmental disorders such as Fragile X Syndrome (FXS). How developmental changes in neuronal function culminate in network dysfunction that underlies sensory hypersensitivities is unknown. By systematically studying cellular and synaptic properties of layer 4 neurons combined with cellular and network simulations, we explored how the array of phenotypes in Fmr1-knockout (KO) mice produce circuit pathology during development. We show that many of the cellular and synaptic pathologies in Fmr1-KO mice are antagonistic, mitigating circuit dysfunction, and hence may be compensatory to the primary pathology. Overall, the layer 4 network in the Fmr1-KO exhibits significant alterations in spike output in response to thalamocortical input and distorted sensory encoding. This developmental loss of layer 4 sensory encoding precision would contribute to subsequent developmental alterations in layer 4-to-layer 2/3 connectivity and plasticity observed in Fmr1-KO mice, and circuit dysfunction underlying sensory hypersensitivity.

Formation of memory assemblies through the DNA-sensing TLR9 pathway

Article Open access 27 March 2024

Vladimir Jovasevic, Elizabeth M. Wood, … Jelena Radulovic

Simultaneous single-cell three-dimensional genome and gene expression profiling uncovers dynamic enhancer connectivity underlying olfactory receptor choice

Article Open access 15 April 2024

Honggui Wu, Jiankun Zhang, … X. Sunney Xie

Neuronal activity rapidly reprograms dendritic translation via eIF4G2:uORF binding

Article Open access 08 April 2024

Ezgi Hacisuleyman, Caryn R. Hale, … Robert B. Darnell

Introduction

Individuals affected by many types of Autism Spectrum Disorder (ASD) and Intellectual Disabilities (ID) commonly exhibit sensory perceptual disturbances and tactile reactivity that span multiple modalities^1,2,3. Fragile X syndrome (FXS) is a leading heritable cause of ASD/ID⁴ with symptoms including seizures, tactile hypersensitivity and abnormal behaviours that affect early sensory and cognitive development. FXS is caused by loss of FMRP protein following transcriptional silencing of the Fmr1 gene. Like FXS, the Fmr1-KO mouse model⁵ lacks FMRP and exhibits sensory, behavioural and cognitive deficits^6,7. Sensory dysfunction in FXS and related ASDs have been proposed to underlie a range of behavioural and cognitive symptoms. In support of this hypothesis, a recent study has indicated a causal link between sensory dysfunction and social and repetitive behaviours in a mouse model of autism⁸. Hence a detailed understanding the sensory function in FXS may be critical to developing novel therapies.

Rodent models demonstrate that the sensory hypersensitivities associated with Fmr1 deletion are mirrored by an increase in circuit excitability^{9,10,11,12,13}. However, numerous cellular processes contribute to circuit hyperexcitability in Fmr1-KO mice¹⁴ and the potential number of mechanisms is even greater. FMRP has the potential to regulate the translation of diverse classes of neuronal mRNA^15,16,17 including ion channels, neurotransmitter receptor subunits, and intracellular signalling molecules. Furthermore, protein–protein interactions with voltage gated ion channels directly link FMRP to maintenance of intrinsic neuronal properties^11,18,19,20. Finally alterations of the excitatory/inhibitory (E/I) balance, proposed to form a key component of circuit dysfunction in neurodevelopmental disorders^21,22 have been reported in Fmr1-KO mice^{9,23,24,25,26,27} although the precise cellular mechanism underlying this E/I imbalance is unknown. Hence, a detailed dissection of the contribution of individual cellular phenotypes underlying the emergent circuit pathophysiology is required to understand the sensory processing deficits associated with FXS¹⁰.

In the somatosensory cortex of Fmr1-KO mice, circuit dysfunction arises very early in development correlating with the peak of FMRP expression during the second postnatal week²⁸, a developmental stage marked by both the end of critical period refinement for thalamocortical (TC) synaptic input and the coordinated maturation of cortical layer 4 cell-intrinsic properties and recurrent network circuitry^{29,30,31,32,33}. Loss of FMRP delays the onset and termination of the critical period for synaptic plasticity at TC synapses²⁸. The consequences of this delay, both in terms of cellular and circuit function, are unknown, however, it is notable that active whisking begins soon after³⁴. Hence, changes in the cellular and circuit physiology in layer 4 at this age could dramatically alter the nature of the sensory information being transmitted to layer 2/3 that drives further experience-dependent development. Importantly, at later ages Fmr1-KO mice display disrupted functional connectivity in layer 4 and an altered synaptic E/I balance^24,35,36,37. However, whether these differences arise as a direct result of the loss of FMRP or are compensatory changes resulting from earlier developmental alterations in cellular physiology is not known. Furthermore, it is unclear if the cellular abnormalities that result from a delay in the sensitive period for synaptic plasticity at TC synapses underlies the altered E/I balance.

To address these questions, we examined the cellular and circuit properties of excitatory and feed-forward inhibition providing Fast-spiking GABAergic neurons (FS) in layer 4 of barrel cortex at postnatal days 10–11 (P10–11)^29,38,39, immediately after the termination of the delayed sensitive period for synaptic plasticity in Fmr1-KO mice and immediately prior to the onset of whisking behaviour. By combining brain slice electrophysiology and computational modelling, we show that an array of cellular-level pathologies is observed in Fmr1-KO layer 4: in connectivity, in cellular intrinsic properties and in synaptic function. The net result of these pathologies is a circuit with a lower threshold for action potential generation in response to TC input, but less well-timed firing relative to input stimuli. Our modelling shows that the cellular-level pathologies observed in the Fmr1-KO are often antagonistic in terms of circuit function suggesting that some ‘pathologies’ rather may be compensatory adaptations. Despite the compensation, the layer 4 circuit is dysfunctional in Fmr1-KO with a reduced ability to encode information manifesting as a reduction in pattern classification accuracy as relayed to layer 2/3.

Results

Altered membrane properties and excitability in Fmr1-KO mice

Previous work has identified a number of cellular, synaptic and circuit changes in Fmr1-KO mouse barrel cortex^9,24,40. However, these experiments were conducted 5 days after the end of the delayed critical period for long-term potentiation (LTP) in Fmr1-KO mice. Therefore, it is unclear whether these changes are pathological or compensatory to network function or indeed how these cellular and synaptic mechanisms interact to produce circuit level deficits. To address this, we systematically analysed cellular, synaptic and network changes in thalamocortical brain slices^41,42 from younger Fmr1-KO mice immediately following the period of LTP at thalamocortical synapses and evaluated the effects of these mechanisms at a circuit level using computer simulations of layer 4 barrel cortex.

We first investigated whether there were changes in passive membrane properties in layer 4 barrel cortex neurons in acute slices from P10/11 Fmr1-KO mice compared with wild-type (WT) littermates, using whole-cell patch-clamp recordings. The principal cell type in layer 4 barrel cortex is the stellate cell (SC), which are recurrently connected glutamatergic neurons that project to layer 2/3^{43,44,45,46,47}. Assessing intrinsic excitability, compared with WTs, Fmr1-KO SCs required less injected current to fire an action potential (lower rheobase) and exhibited an increase in both input resistance and membrane time constant (Fig. 1a), but no change in membrane capacitance or resting membrane potential (Fig. S1). SCs in Fmr1-KOs also had enhanced excitability with an increase in the number of action potentials elicited by 500 ms depolarising current steps (Fig. 1b); however, this increase in action potential number was associated with a decrease in action potential frequency during the early part of the depolarisation resulting in a reduction in action potential frequency adaptation, as assessed over a duration necessary to fire a standardised number of spikes (Fig. 1c). Moreover, in agreement with a previous finding in Fmr1-KO mouse hippocampus²⁰ action potential kinetics were slowed in Fmr1-KOs, with an increase in width and a decrease in amplitude (Fig. 1d).

A similar analysis was performed on the intrinsic excitability of fast-spiking interneurons. FS provide strong feed-forward inhibition (FFI) onto SCs and play a critical role in determining the integration of TC input and action potential output of SCs^{30,31,38,48,49}. Similar to SCs, we observed an increase in input resistance and membrane time constant and reduced rheobase (Fig. 1e), and no change in whole-cell capacitance or resting membrane potential (Fig. S1b). However, in contrast to SCs, FS exhibited a reduction in action potential number during depolarisation (Fig. 1f) accompanied by a decrease in both action potential frequency and frequency accommodation (Fig. 1g) and a slowing of action potential kinetics (Fig. 1h, i).

Thus, in layer 4 of P10/11 Fmr1-KO mice both SCs and FS have altered passive membrane properties such that they produce action potentials in response to less depolarisation, compared with wild type. During sustained depolarisation, Fmr1-KO SCs continue to exhibit increased action potential output albeit with alterations in frequency and kinetics. In contrast, FS show decreased action potential output in response to sustained depolarisation as well as changes in frequency and kinetics.

Increased response to depolarisation in SCs in Fmr1-KO mice

SCs are the output neurons from layer 4 and as such play an important role to integrate ascending TC input and provide output to layer 2/3, encoding stimulus features such as the presence of an object and its surface quality in spike frequency and times^38,39. The capacity of neurons to faithfully follow rhythmic synaptic input, i.e. their impedance, is fundamentally constrained by their intrinsic properties⁵⁰. We therefore next investigated whether the changes in intrinsic properties of SCs in Fmr1-KO mice alter their ability to transform inputs at ethologically relevant frequencies associated with somatosensation^{51,52,53,54,55}. We assessed SC impedance properties^56,57, using a sinusoidal current of progressively increasing frequency and measured the resulting voltage response. SCs from Fmr1-KO mice exhibited an increased impedance between 0.5 and 4 Hz (Fig. S2a), consistent with predictions from the change in passive membrane properties (Fig. S3). This frequency-dependent response suggests that SCs in Fmr1-KO mice should exhibit alterations in action potential generation in response to membrane depolarisations of the same frequency range. To test this, we applied suprathreshold sinusoidal depolarisations at specific frequencies and found an increase in the number of action potentials elicited at low frequencies in Fmr1-KO compared with wild-type mice (0.5–4 Hz; Fig. S2B). Moreover, a significant phase shift in the timing of action potential generation was also observed, for example at 10 Hz (Fig. S2C).

Thus, the alteration in the passive membrane properties of SCs in Fmr1-KO mice results in a selective increase in membrane depolarisation in response to low frequency depolarising currents. This produces both an increase in action potential output and a shift in action potential timing in response to low frequency depolarising currents.

Altered feed-forward inhibition on SCs in Fmr1-KO mice

FFI is a critical mechanism that governs the integration of TC input by SCs^{30,31,32,48,58}. In layer 4, barrel cortex FFI is mediated by FS neurons which provide the TC-evoked inhibition onto SCs. This form of evoked E/I balance sets an integration window for SC action potential output between the onset of excitatory and the strong but delayed inhibitory inputs⁴⁸. To determine if there are alterations in FFI in Fmr1-KO mice, we first investigated connectivity between FS and SCs using simultaneous whole-cell patch-clamp recordings. Similar to findings by Gibson et al. at P14²⁴, at P10/11, we found reductions in connectivity from SC to FS, but additionally from FS to SC and reciprocally between SC and FS, in Fmr1-KO mice compared with wild types (Fig. 2a). In contrast, however, the strength of the connections between FS and SCs was unchanged in Fmr1-KO mice compared with wild types (Fig. 2b) and the kinetics of the IPSCs onto SCs and of EPSCs into FS were also unchanged (Fig. 2c). Furthermore, there was no difference in SC–SC connectivity in Fmr1-KO mice compared with wild type, nor in synaptic strength or kinetics (Fig. 2d–f).

We next directly measured FFI in SCs, by stimulating TC inputs and comparing the size of the monosynaptic TC-evoked EPSC (at −70 mV holding potential) and the size of the di-synaptic feed-forward IPSC in the same SCs (at 0 mV, following methods in ref. ³²). Here, we report the E/I balance as the strength of evoked inhibition normalised to the driving excitatory input strength: In SCs from Fmr1-KO mice, we found an increase in the peak amplitude of the feed-forward IPSC relative to the monosynaptic EPSC (‘GABA/AMPA ratio’ [G/A]; Fig. 3a). Notably however, a subset of Fmr1-KO neurons lacked any FFI, whereas FFI was observed in all tested WT neurons. The feed-forward IPSC also exhibited slower kinetics and a longer onset lag time (compared with the monosynaptic EPSC) in SCs from Fmr1-KO mice (Fig. 3b).

To investigate the consequences of these changes in the relative magnitude and timing of FFI we measured the resulting postsynaptic potential (PSP, Fig. 3c). As we and others have previously shown^32,48,59, the postsynaptic potential (PSP) half-width in SCs is strongly influenced by FFI, which truncates the PSP. We found that PSP full width at half-maximum amplitude (half-width) is increased in the Fmr1-KO mice (Fig. 3d) suggesting a decrease in functional FFI, despite the increased feed-forward inhibitory synaptic input onto SCs. This reduction in functional FFI is likely due to the changes in passive membrane properties of SCs; the mechanism is further explored as described later. One additional mechanism that could contribute to the lack of increase in functional FFI in SCs is a change in chloride reversal potential in the Fmr1-KO mice. To test this, we used perforated patch-clamp recordings and found no difference in chloride reversal potential in SCs from the two genotypes (Fig. S4).

These findings, therefore, show that there is diminished functional FFI onto SCs in Fmr1-KO mice despite an increase in G/A ratio to a single stimulus. Considering the critical role feed-forward inhibition plays in determining action potential generation and timing in SCs^30,48,60, this deficit is likely to have important consequences for layer 4 network function at this crucial stage of activity-dependent development.

Enhanced short-term synaptic depression in Fmr1-KO mice

Short-term synaptic plasticity is an important mechanism for information processing in cortical networks^61,62 including determining the efficacy of FFI^30,31,32,48. We therefore compared short-term depression of the TC excitatory input and of the feed-forward TC-evoked IPSC onto SCs. Both excitatory and feed-forward inhibitory transmission in the Fmr1-KO mice showed an increase in short-term depression during trains of five stimuli over a frequency range from 5 to 50 Hz compared with WT (Fig. 4a). However, the relative increase in short-term depression was greater for the IPSC than the EPSC leading to a stronger increase in E/I ratio during trains of activity in Fmr1-KO compared with wild-type mice (Fig. 4b, c). We also investigated short-term plasticity at SC-FS and FS-SC synapses in connected simultaneous recordings. In contrast to the normal short-term plasticity observed by Gibson et al.²⁴ in Fmr1-KO recordings at 2 and 4 weeks old, we found increased short-term depression at P10/11, both for the excitatory input onto FS from SC and for the inhibitory input from FS to SC (Fig. 4d, e). In SC to SC connections, although there was a trend for a decrease in short-term depression during the train, there was no robust change observed (Fig. 4f, g).

Increased postsynaptic summation and spiking in Fmr1-KO mice

SCs are the output neurons for layer 4 barrel cortex and serve to integrate TC input and transform it to action potential output, modulated by local feed-forward inhibition^46,49,63. Our results thus far reveal a number of cellular and synaptic phenotypes in layer 4 barrel cortex neurons in the Fmr1-KO mice that likely will effect SC integration of TC input. For FFI, although we show that the G/A current ratios are typically elevated in SCs from P10–11 Fmr1-KO mice compared with WT, the activity-dependent depression of FF inhibitory input is faster and greater in KO SCs compared with the depression of thalamic EPSCs (group data in Fig. 4c). Coupled with the increased intrinsic excitability of the KO neurons, this effect is predicted to produce an increase in voltage summation and contribute to abnormal AP generation in response to thalamic input.

Therefore, to investigate the net effect of these phenotypes on the layer 4 circuit we studied the membrane potential response in SCs elicited by TC stimulation. Using whole-cell patch-clamp recordings in current clamp, we found that trains of five stimuli at frequencies from 5 to 50 Hz produced greater postsynaptic summation in SCs from Fmr1-KO mice (Fig. 5a) and a larger fraction of stimuli eliciting action potentials (Fig. 5b, c). We used cell-attached patch recordings to make less invasive measure of spike output. Using this approach, we confirmed that SCs from Fmr1-KO mice produce spike output at lower TC stimulus frequencies than wild type. We further found that the timing precision of action potential firing during stimulation (50 Hz), a physiologically relevant frequency, was strongly impaired, exhibiting a lower instantaneous frequency and greater variability in rate and timing (Figs. 6a–d, S5). These compound effects on trial-to-trial fidelity and lower firing rate by Fmr1-KO neurons reduced their population average spike density function (Fig. 6b, c individual neurons shown in S5C,D). No significant genotype-dependent changes were observed in the number of spikes fired per trial. Thus, the net effect of the cellular and synaptic changes in the Fmr1-KO mice is that SCs that produce action potentials more readily, but with less precise timing.

Antagonistic cellular and synaptic changes in Fmr1-KO mice

Our experimental data show that loss of Fmr1 produces a variety of cellular and synaptic alterations leading to an overall change in the transformation of TC input to action potential output by SCs. The various mechanisms we describe likely interact in a complex manner to produce the overall circuit phenotype. To better understand this, we first developed a single-cell model of thalamocortical integration to systematically explore the interaction of our findings in Fmr1-KO SCs. We used a leaky integrate-and-fire model neuron with a single dendritic segment, receiving bulk glutamatergic synaptic input connected to a soma receiving somatic GABAergic input, representing thalamocortical and FFI inputs, respectively. We adjusted the intrinsic and synaptic kinetics and short-term plasticity behaviour of the model to match WT and Fmr1-KO mean values and altered the strength and timing of the GABAergic input (relative to that of the glutamatergic) to explore the parameter space described by our experimental data (Fig. 7a). We asked which frequencies of train stimulation between 5 and 50 Hz led to spike firing by the model (Fig. 7b). In the absence of FFI (G/A = 0, GABAergic conductance silenced) both WT and Fmr1-KO models fired spike(s) at all tested frequencies. Gradual addition of FFI progressively restricted the ability of trains to elicit spikes to those at higher frequencies, but this effect was less pronounced in Fmr1-KO model compared with WT. Notably, the range of G/A over which this effect was observed are similar to the G/A ranges we measured experimentally. The Fmr1-KO model also reproduced the experimentally observed increases in first spike delay, spike timing jitter and spike number compared with WT (Fig. 7b insert, c).

We used our single-cell model to explore the role of groups of parameters in regulating spike output generated by thalamocortical input over a range of G/A values (Fig. 8). To simplify the model parameter landscape to those involving similar physiological processes, we assigned parameters into four groups. These were namely: short-term plasticity (coefficients constraining synaptic depression behaviour of TC and FFI inputs during repetitive stimulation), excitatory-inhibitory synaptic input delay (lag between simulated TC and FFI input onset times), intrinsic excitability (leak conductances and capacitance of model), and synaptic kinetics (rise and decay kinetics of synaptic conductance). We then systematically replaced Fmr1-KO values with WT values for the various parameter groups to simulate different rescue scenarios and measured spike output over a range of G/A values for all of these combinations (four parameter groups to the power of two potential genotypes = 16 simulated scenarios) in response to 5 stimuli at 5–50 Hz. For example, when we ‘rescued’ excitatory-inhibitory synaptic input delay on the Fmr1-KO background there was an 86% increase in spike output by the model across the parameter space explored, relative to the number of spikes fired by the WT model. Conversely, rescuing intrinsic excitability on the Fmr1-KO background produced a 26% reduction in spike output. Thus, manipulation of different parameter groups can have opposing effects on spike output. Overall, by rescuing various combinations of parameter groups we observed that some have antagonistic effects on spike output (e.g. simultaneous rescue of intrinsic excitability and excitatory-inhibitory synaptic input delay) while others drive spike output in the same direction (e.g. excitatory-inhibitory synaptic input delay and short-term plasticity). This analysis indicates that the elevated intrinsic excitability observed in Fmr1-KO SCs is a dominating feature driving the changes in spike output. However, changes to other parameters limit the physiological manifestation of this cellular phenotype at the output level of spike generation, suggesting that some of the parameter changes observed in the Fmr1-KO thalamocortical phenotype reflect a compensation in response to the underlying pathology that limits circuit dysfunction.

Network dysfunction in layer 4 in Fmr1-KO mice

Having explored the mechanisms underpinning the altered transformation of the thalamocortical input in layer 4 Fmr1-KO, we next asked how the cellular and synaptic deficits result in altered network activity and what the consequences are for information processing by layer 4. We developed a spiking network model of a cortical layer 4 barrel, based on our electrophysiological recordings (Fig. 9a and S6; see also Methods and Supplemental Table 1). The model consists of a network of interconnected SCs and FS, stimulated by TC input. The great advantage of layer 4 barrel cortex in this regard is the numbers of SCs and FS in a barrel are known as is their probabilistic connectivity⁴³. The model was randomly connected, constrained by the experimentally determined probabilities (as detailed above) for WT and Fmr1-KO phenotypes. Cells were represented as single-compartment leaky integrate-and-fire (LIF) neurons connected by inhibitory (GABA_A) and excitatory (AMPA and NMDA) synapses displaying short-term plasticity. In addition, model parameter distributions matched recorded biological diversity in intrinsic neuronal properties, synaptic strength and kinetics, transduction delays and short-term plasticity of synapses (Fig. S6). Simulated membrane depolarisation and action potentials were measured in response to the same patterns of train stimulation of TC axons as we used in the slice experiments.

Ab initio, the random model reproduced well the main features of our experimental data in both WT and Fmr1-KOs, both at the single cell and at the network level. The Fmr1-KO displayed increased subthreshold summation and increased spiking activity to TC input at low train frequencies at the single-cell level (Fig. 9b, c). Similar to our experimental findings and in agreement with our single-neuron thalamocortical model, the timing of first spikes fired in the Fmr1-KO model was delayed and showed decreased inter-trial fidelity (Fig. S7a). Similarly, fluctuations in trial-to-trial mean firing rate were enhanced in the Fmr1-KO model (Fig. S7b). Notably, approximately half of the per-trial-averaged numbers of spikes fired by neurons in the recurrently connected model could be accounted for by those fired by our single-cell model lacking recurrent connectivity (Fig. 7c cf. Fig. S6a) demonstrating a strong thalamocortical contribution to the generation of spiking in layer 4.

At the network level, high-frequency TC stimulation evoked a transient reverberating population response, characterised by high rates of simultaneous model neuron firing (highlighted by peaks in the grey firing histograms overlaid in Fig. 9d). In the WT model these synchronous population events were initially tightly locked to latter stimuli during the train before transitioning to self-sustaining network activity persisting after the last stimulus. Fmr1-KO model network activation required fewer stimuli and was recruited immediately following the first TC-locked spike volley but did not produce the self-sustaining network activity observed in WT. We explored this observation further by examining the stimulus-evoked population firing rates of excitatory and inhibitory neurons in the model network (Fig. 9e). The first TC stimulus led to spiking in a greater fraction of Fmr1-KO inhibitory neurons compared with WT. Strong reductions in the fraction of firing of inhibitory neurons were observed for successive stimuli for both genotypes, but this was much greater for the Fmr1-KO model, consistent with experimental data. In the WT model, firing rate increases in the excitatory and inhibitory populations were strong and temporally matched, leading to a slowly decaying limit cycle (Fig. 9f) for recurrent network firing. In contrast, the rates of synchronous firing for excitatory and inhibitory populations were disorganised in the Fmr1-KO model, with lower overall synchrony, producing unstable limit cycles.

Overall, the network model captures many of the important features of our experimental data and suggests a mechanism by which single-cell physiological defects in the Fmr1-KO layer 4 network (particularly weaker Ex–In connectivity and cellular hyperexcitability) generate impaired network activity.

Impaired pattern discrimination in a Fmr1-KO layer 4 model

Thus far we have demonstrated that a simple data-driven network model of Fmr1-KO layer 4 reproduces the main experimentally obtained differences in responses to patterned TC input. Since a major role of barrel cortex layer 4 is the representation of sensory details in the rate and timing of excitatory neurons action potentials⁶⁵, we next examined how individual neurons and small cell ensembles in our WT and Fmr1-KO layer 4 models performed when challenged to detect subtle temporal features of simulated sensory input (Fig. 10). This was represented as an extra stimulus inserted into the input stimulus train and can be conceptualised as the perturbation of TC neurons’ firing patterns in response to a small surface bump encountered during active sensation.

The majority of WT Ex neurons detected the oddball stimulus as evidenced by a change in firing rate (indicated by the off-diagonal clustering of points in Fig. 10a). However, the majority of simulated Fmr1-KO neurons were less likely to detect the oddball stimulus. In addition, WT neurons represented the presence of an oddball stimulus in the timing of the first spike fired (Fig. 10b) such that, despite an interdependence, WT neurons could encode the presence of an oddball as a change in spike rate and/or a shift in first spike timing. Conversely, Fmr1-KO neurons showed a reduced change in first spike timing in response to the oddball stimulus. To quantify both of these spike train effects, we used a distance metric that encompasses both spike train and rate^66,67 to demonstrate the distribution of single-cell oddball stimulus representations (Fig. 10c). The distribution of regular-oddball spike train dissimilarities was uniformly reduced in the KO, implying that the net effect of spike train rate and timing leads to a global reduction across the network in the ability to discriminate a fine temporal feature of simulated sensory input.

The principal output from the layer 4 circuit is a sparse population code provided by the ascending synaptic connection to layer 2/3 neurons^{51,52,68,69,70,71}, the activity of which plays an instructive role in the development of the layer 2/3 network^44,72. Therefore, it is important to consider how the defects in layer 4 stimulus representation might relay at the subsequent processing level. We therefore used a linear population decoder⁷³ to classify input stimulus patterns based on the layer 4 output firing pattern (Fig. 10d). For population codes formed from ensembles of between 10 and 500 neurons (i.e. representing convergence of approximately 1/80 to 2/3 of model layer 4 neurons to a layer 2/3 ‘readout’ neuron), more neurons were required for the Fmr1-KO model to reach the same degree of classification accuracy as the WT model (Fig. 10e). Even drawing 500 neurons (~2/3 of the total excitatory network), population coding of input pattern was impaired in the Fmr1-KO model.

Taken together, this result suggests that the developing synaptic pathway from layer 4 to layer 2/3, as schematised here as the layer 4 network’s input pattern classification accuracy would be mis-informative about the quality of sensory input in the Fmr1-KO model, and that information gleaned from a larger pool of layer 4 neurons would be necessary for an informative representation of the sensory input.

Discussion

Altered sensory reactivities are a common and debilitating feature of FXS and of ASDs more generally (reviewed in refs. ^1,74). Sensory symptoms are often distressing for affected individuals as they lead to anxiety and social withdrawal, as well as self-injurious behaviour. Unfortunately, very little is known about when and how altered sensory function arises during development, nor how altered neuronal properties affect experience-dependent development.

Here, we dissect the cellular mechanisms by which FMRP deletion results in a paradoxical increase in both G/A ratio and feed-forward summation of TC input induced layer 4 SC EPSCs (see Supplementary Discussion for comparison to existing literature on altered FFI in Fmr1-KOs). Despite this increase in local excitability in the layer 4 feed-forward circuit, stimulation of TCs at behaviourally-relevant frequencies leads to decreased and disorganised spiking of L4 SCs during recurrent activity. Modelling this range of phenotypes predicts an altered temporal coincidence detection of thalamocortical inputs at a range of physiologically relevant stimulation frequencies, a prediction that is borne out by experimental data. Together these data indicate that the sensory information being relayed to layer 2/3 is dramatically altered at the commencement of whisking (P12), which would have a deleterious effect on experience-dependent development and hence sensory processing at later ages. These findings also highlight the power of integrating experimental and mathematical approaches to understanding and forming testable hypothesis about the nature of the developmental mechanisms underlying neuronal circuit dysfunction in ASD.

An altered E/I balance has been causally linked to a range of psychiatric disorders from ASD to schizophrenia^{21,75,76,77,78}. In primary somatosensory cortex, the E/I balance is essential for regulating both the frequency (spike rate) and quality (spike timing) of action potentials, and hence for accurately processing sensory input⁶⁵. In particular, the precise temporal tuning of excitation and inhibition is critical to signal processing of cortical networks⁷⁹, as reported elsewhere⁸⁰. Here, we find a reduction in the E/I balance in a subset of layer 4 SCs that causes an overall decrease in the E/I at the population level despite the observation of a subset of Fmr1-KO SCs that completely lack functional FFI. The reduction in the E/I balance is accompanied by a gain in overall excitability because of elevated SC intrinsic excitability and an increase in the delay between TC excitation and FFI onto spiny stellate neurons (Ex–In lag). This delay results in a broadening of the TC EPSP and a large increase in the summation of EPSCs following TC stimulation that is especially pronounced at higher stimulation frequencies. Compounding this increase in excitatory summation, we find a more dramatic loss of inhibition relative to excitation during physiologically relevant repetitive stimulations. The rapid decrease of inhibition results in a further increase in summation of the EPSPs. These data highlight the vagaries of simply examining E/I imbalance as a key characteristic of circuit dysfunction. It is important to remember that the E/I balance is dynamically regulated^{51,79,81,82,83} and is only one factor that regulates neuronal processing within a complex circuit and its role in mediating sensory dysfunction^9,10. These nuanced and dynamic circuit effects further highlight the need for reductive circuit modelling in disentangling the processes underlying circuit dysfunction in animal models of psychiatric disease^80,84.

Gibson et al.²⁴ have previously reported normal thalamocortical FFI in the Fmr1-KO mouse at P14 and P21, implicating impaired excitatory drive onto layer 4 FS neurons as the mechanism underlying layer 4 network hyperexcitability. The disparity between these findings and the present study could result from the age of the animals being examined and the static nature of single-stimulus responses to thalamic stimulation employed by Gibson et al., which lacks the information of the dynamic circuit behaviour invoked by our stimulation paradigms.

Our finding of an increase in input resistance of layer 4 SCs is in good agreement with previous findings for this cell type²⁴ and is perhaps the most consistent physiological feature across neuronal cell types in the Fmr1-KO mouse^24,85,86,87. It is important to note that several Fmr1-KO studies spanning prefrontal and entorhinal cortices as well as hippocampus have demonstrated altered spike firing and excitability in the absence of altered input resistance^88,89,90, indicating that alteration in conductances contributing to cellular input resistance are only one (potentially co-occurring) mechanism underlying distorted cell-intrinsic response to input. A further pertinent observation supporting this claim comes from the MeCP2 mouse model of Rett syndrome^91,92, in which levels of cortical neuronal firing are dampened by reduced E/I tone in the absence of any change to intrinsic cellular excitability. In the present study, elevation in input resistance and resultant slowing in membrane time constant raises the impedance to sinusoidal stimuli. By increasing the intrinsic sensitivity of Fmr1-KO SC neurons to low frequency stimulation this has the compound effect of many more action potentials being generated in response to low frequency stimuli, and a dampening of responsivity to high-frequency inputs. Importantly, these input resistance-driven changes in intrinsic properties also result in a phase shift in the voltage response to input currents, for example producing a phase shift in action potential generate. Such a shift in phase can have dramatic effects on spike-time dependent plasticity (STDP⁹³) which may contribute to the shift in the timing of the critical period for LTP at thalamocortical synapses²⁸. Crucially, the maturation of the layer 4 excitatory network in the brief developmental window beyond P10 is sensitive to sensory experience³³ and relies on spike timing-dependent plasticity mechanisms⁹³. The cellular and circuit disruptions to spike-time representation reported here would be expected to dramatically alter experience-dependent development of primary somatosensory cortex in Fmr1-KO mice and could contribute to the circuit hyperexcitability reported at older ages. Consistent with this idea, impaired STDP is reported in the Fmr1-KO mouse prefrontal cortex⁹⁴.

Our findings indicate that the circuit deficits in layer 4 of somatosensory cortex of Fmr1-KO mice arise from a range of cellular deficits, however, their relative contribution to the overall layer 4 circuit function and hence the information propagated to layer 2/3 is not known, prompting our use of circuit modelling.

The distributed nature of the physiological changes that result from the loss of FMRP suggest that compensatory changes may be activated even at these early stages of development. However, it is not clear whether the phenotypes reported here, or indeed many of those reported for other cell types (for review see ref. ¹⁴), arise as a direct result of the loss of FMRP or homeostatic compensation in an attempt to restore circuit function. In this context, the recent work of Antoine et al.⁸⁴ is particularly relevant. They demonstrated that altered E/I balance in layer 3 of primary somatosensory cortex was a common feature of multiple ASD mouse models. Furthermore, modelling indicated that the changes in E/I balance was tuned to stabilise synaptic drive and spike output for cells near threshold strongly indicating that altered E/I balance was a homeostatic response of the neurons to normalise spiking. For the phenotypes described in this study, it is not possible to determine which cellular features are directly related to the loss of FMRP and which reflect homeostatic changes. However, what is clear is at this early stage of development, E/I balance is not tuned sufficiently to normalise spike output of layer 4 SCs. Future studies examining later stages of development combined with their electrical and/or pharmacological manipulation will be needed to determine the role of each of these physiological parameters in determining the overall circuit dysfunction and the overall hyperexcitability phenotype. It should be noted that, irrespective of causality, each of these phenotypes could provide a potential for therapeutic intervention for FXS but the efficacy of treatment is likely to be specific to particular developmental time-points.

Recent studies have attempted to pharmacologically normalise abnormal E/I balance in mouse models of diverse forms of intellectual disability^{92,95,96,97,98,99}. Our findings indicate that E/I balance may be a poor proxy for underlying cellular disturbances and subsequent circuit dysfunction and hence should be approached with caution when devising and determining the effectiveness of potential treatment strategies.

The approach taken in this study was to build a mechanistic understanding of a cellular and circuit dysfunction through combined and bidirectional physiology and predictive modelling. By studying seemingly small pathophysiological changes in the component parts of a well-described dynamic circuit, our models were able to predict large physiological increases in the sensitivity to ethological stimuli. The local circuit model lacks the recurrently connected architecture of L4, thus it is useful to disentangle the contributions of direct vs. feedback processing on the firing patterns of the SCs. The single-cell FFI model specifically predicted the altered summation properties at a range of stimulation frequencies that were verified through experimentation. Furthermore, by substituting each parameter or combination of parameters from Fmr1-KO animals with that from their WT counterparts, we were able to begin examining the relative contribution of each phenotype to the overall local circuit dysfunction.

This model also provides a powerful framework to study compensation within neocortical neurons that arise from particular genetic mutations. It could also be used to examine the effects of targeted treatments, for example by simulating the effect of a single parameter repair or combination repair strategy, as we do in the current study.

To encompass the recurrent circuit activity initiated by thalamocortical activity, we generated a network model to investigate the nature of the output of the of the layer 4 circuit. The primary goal of this model was to the understand the quality and fidelity of the information leaving layer 4 in Fmr1-KO compared with WT which is key to driving experience-dependent development of layer 2/3 neurons⁴⁴. Using data-tuned cell models as a starting point, we incorporated experimentally derived connection probabilities with known layer SC and FS cell numbers and GABA, NMDA and AMPA currents. The model predicts a decrease in spike firing and loss of trial-to-trial fidelity in the Fmr1-KO animal, accurately reflecting our data obtained using patch-clamp recordings from slices. To simulate a subtle temporal sensory disturbance during active whisking, we examined the effect of an oddball stimulus introduced during the stimulus train. The WT model was able to detect the presence of an oddball stimulus through a change in the firing rate and timing of the spikes fired by Excitatory neurons. On the contrary, the Fmr1-KO model was less likely to alter its firing rate and the timing of the first spike in response to the oddball stimulus, indicating a robust decrease in the circuit’s ability to discriminate changes in the sensory environment. These findings strongly suggest that the ability of the layer 4 circuit to relay details about patterned sensory information, a key feature of experience-dependent plasticity, is dramatically reduced in Fmr1-KO mice and make some strong and testable predictions for human psychophysical experiments and their animal experimental in vivo counterparts. The consequence of reduced sensory driven plasticity would be a retardation of circuit development. Importantly, Bureau and colleagues¹⁰⁰ found a delay in the developmental increase in connection probability between layer 4 and layer 2/3 cells in somatosensory cortex of Fmr1-KO mice that was mirrored by a delay in the pruning of SC axons projecting to layer 2/3. This delay could be accounted for the shift in the sensitive period for inducing LTP at TC inputs in layer 4 in somatosensory cortex which is delayed in Fmr1-KO animals²⁸. Importantly while some key features of somatosensory cortex development are delayed, this is not a global phenomenon. For example, layer 4 SCs show typical developmental trajectories for dendritic elaboration, spinogenesis, and synaptogenesis^28,101. They also show the typical activity-dependent restriction of their dendrites to the regions of thalamocortical axon terminals¹⁰¹. These findings provide support for the idea that a mismatch in the maturation of physiological and/or anatomical properties that are necessary for the subsequent development of circuit function that may underlie pathology in FXS¹⁰².

The mismatch in developmental trajectories demonstrated here may subsequently induce altered development in downstream areas. At P10, we find a dramatic reduction in the information being propagated to layer 2/3 at this age. Hence, for a layer 2/3 neuron, the loss of FMRP essentially mimics sensory deprivation^31,32,33,93. Intriguingly, Gainey et al.¹⁰³ found a compensatory change in E/I balance following experience-induced plasticity to whisker trimming similar to that found by Antoine et al.⁸⁴ in layer 2/3 of older Fmr1-KO mice. Hence the changes in E/I balance seen in layer 2/3 at P21 may result from the altered circuit activity we observe here. Several studies have examined changes in somatosensory circuit function in Fmr1-KO mice at more advanced stages of development^9,24,26,84 and it will be important to determine whether these alteration could similarly arise from the altered circuit activity at younger stages of development and the resulting changes in experience-dependent development.

Despite these differences in age between studies, we have confirmed key findings from previous work. Notably, although we report normal synaptic connection probability between pairs of L4 SCs at P10–11, but this was reduced by P14 in Fmr1-KO mice (reported in legend for Fig. 2d) in agreement with Gibson et al.²⁴. We hypothesize that an initial developmental wave of synaptic connectivity emerges unperturbed between SCs in Fmr1-KO L4 despite the delayed critical period for thalamocortical plasticity²⁸ as well as local circuit and cell-intrinsic effects reported in the present study.

Methods

Mice

Fmr1-KO mice (B6.129-Fmr1^tm1Cgr, RRID:MGI_MGI:3815018) were obtained from the Contactor laboratory (Northwestern University, Chicago IL, USA). Mice were maintained on a C57Bl6/J background. Hemizygous male Fmr1-KO (Fmr1^−/Y) and wild-type littermate (WT, Fmr1^+/Y) pups were used at postnatal days 10–11 (P10–11). All procedures were carried out according to UK Home Office and NIH IACUC guidelines for animal welfare. N’s of neurons/animals for each experiment are indicated in the individual figure legends.

Brain slice preparation

Mice were decapitated and brains were rapidly removed and placed in ice-cold carbonated (95% O₂/5% CO₂) cutting solution. Thalamocortical brain slices (400-μm thick) were prepared according to reference⁴¹ using a vibrating microslicer (Leica VT1200). Slices were transferred to artificial cerebro-spinal fluid (aCSF) and stored at room temperature. aCSF contained (in mM): NaCl 119, KCl 2.5, NaH₂PO₄ 1, NaHCO₃ 26.5, Glucose 11, MgSO₄ 1.3, CaCl₂ 2.5. Cutting solution was identical except for substitution with 7 mM MgSO₄ and 0.5 mM CaCl₂.

Electrophysiology

Patch-clamp recordings were performed using internal solutions containing (in millimolar): CsMeSO₄ 130, NaCl 8.5, HEPES 10, EGTA 0.5, Mg-ATP 4, Na-GTP 0.3, QX-314 (Tocris) 5 (for voltage-clamp recordings) or KMeSO₄ 130, NaCl 8.5, HEPES 5, EGTA 0.5, Mg-ATP 4, Na-GTP 0.3 (for current-clamp recordings) pH 7.4 at 32 ± 0.5 °C, adjusted to 285 mOsm. Junction potential was uncorrected. Slices were superfused with aCSF saturated with 95% O₂/5% CO₂ at 32 ± 0.5 °C at 8 ml/min, containing in (mM): NaCl 119, KCl 2.5, NaH₂PO₄ 1, NaHCO₃ 26.2, Glucose 11, MgSO₄ 1.3, CaCl₂ 2.5, adjusted to 305 mOsm, pH 7.4 at 32 ± 0.5 °C. Patch recordings were performed with Multiclamp 700B amplifiers (Molecular Devices), using pipettes (4–7 MΩ open bath resistance) fabricated from thin-walled borosilicate glass (Warner Instruments). Signals were low-pass Bessel filtered at 10 KHz and digitized at 20 KHz using either a Molecular Devices Digidata 1322A board and pClamp 10.2 for recordings performed in the USA, or a National Instruments PCI-6110 board, using acquisition software custom-written in C++/MATLAB incorporating modules from Ephus (http://www.ephus.org)¹⁰⁴ for recordings performed in the UK. Custom stimulating electrodes were manufactured using twisted Ni:Cr wire. For TC fibre stimulation, stimulation electrodes were inserted into VB thalamus. Stimulation in VB thalamus provides monosynaptic input to layer 4, and with single stimuli evokes direct excitatory and FFI input to layer 4 SCs. Biphasic, 100 μs constant voltage pulses were delivered by an optically isolated stimulus generator (AMPI systems) at 0.03 Hz inter-pulse frequency. Layer 4 SCs were recorded in barrels that showed >100 μV TC-evoked field potentials. Cells were selected for recording using DIC optics under infrared illumination based upon somatic morphology and laminar position. Resting potential was measured in bridge balance configuration immediately after breaking in. Recording quality was monitored on-line incorporating the following criteria: resting potentials more hyperpolarized than −50 mV, stable (<25 MΩ, <20% drift) access resistance and holding currents <100 pA. Cell-attached recordings were performed according to reference¹⁰⁵. Briefly, stable greater than gigaohm seals were obtained between electrode and neuron with minimal suction to minimise mechanical membrane stress. Recordings were performed in voltage-clamp configuration by holding cells at pipette potentials that were empirically determined in current clamp to minimise holding current. Recordings were accepted if all action currents (spikes) fired were <−200 pA in peak amplitude, seal resistance and holding current were stable for the duration of the experiment. In our hands, the statistics of network-evoked firing responses were stable for over an hour (Figs. 1, 5, 6 and Supplementary Fig. 4 and data not shown), in contrast to a previous study^106,107. For recordings from putative FS interneurons, we targeted neurons with large, cigar-shaped somata in the barrel hollow close to the bottom of L4. Since multiple electrophysiologically defined classes of FS interneurons have been reported in barrel cortex¹⁰⁸, careful inspection of firing patterns was performed before further analysis. Only neurons firing early in current step, with a non-stuttering firing profile were included for analysis in this study (Fig. S1a). Simultaneous whole-cell recordings from two neurons (‘paired recordings’) co-located within the same barrel were performed in the same manner. Synaptic connectivity was tested by injecting a suprathreshold current step into the presynaptic neuron whilst monitoring synaptic currents in the postsynaptic neuron in current clamp at holding potentials of either −70 mV (for Glutamatergic synapses) or 0 mV (GABAergic synapses). A minimum of 20 trials were used to ascertain connectivity.

Impedance profiling

Sinusoidal current waveforms (ZAPs), with $I_{{\mathrm{inject}}}\left( t \right) = A\sin \left( {2\pi ft} \right)$ were injected into current-clamped neurons, with amplitude (A) of ±40 pA and frequency (f) increasing between 0.5 and 50 Hz over 25 s. Fourier-transformed voltage and current waveforms were used to derive the neurons’ complex impedance ${\bar Z\left( f \right)}$ according to ${\bar Z\left( f \right)} = \frac{{{\bar V(f)}}}{{{\bar I}_{{\mathrm{Inject}}}(f)}}$^{109,110,111,112}, where superscripts denote the Fourier operator⁵⁰. The impedance magnitude, $|{\bar Z\left( f \right)}|$ is thus a real valued number between R_membrane for sustained DC current (i.e. Ohmic resistance, ${\bar Z}(f = 0) = R_{{\mathrm{membrane}}}$) and zero, decreasing to zero in the limit of f → ∞. The phase shift ϕ_(f) between the input current and voltage output was taken as the inverse tangent of the ratio between the real and imaginary components of the complex impedance:

$$\phi _{({\mathrm{f}})} = \angle {\bar Z}_{({\mathrm{f}})} = {\mathrm{arctan}}\left( {\frac{{{\bar Z}_{\left( {\mathrm{f}} \right){\mathrm{imaginary}}}}}{{{\bar Z}_{\left( {\mathrm{f}} \right){\mathrm{real}}}}}} \right)$$

(1)

For Bode filter analysis of impedance profiles, the frequency-dependent system gain is expressed as log-power ratios of impedance at higher frequencies relative to that at the lowest tested (f = 0.5 Hz). Little additional output attenuation compared with DC is expected for either genotype at 0.5 Hz, i.e. where $R_{{\mathrm{membrane}}} \approx Z_{{(f} = 0.5{\mathrm{Hz}})}$, hence this is a fair normalization strategy¹¹³. Model first-order RC low-pass filter responses were assumed as:

$$Z = \frac{1}{{\sqrt {1/R_{{\mathrm{membrane}}} + {\mathrm{i}}\omega C_{{\mathrm{membrane}}}} }},{\mathrm{where}}\,{\mathrm{\omega }} = 2{\mathrm{\pi}} f$$

(2)

The bandwidth, or cut-off frequency f_cut-off, of an ideal first-order low-pass filter, above which the output power relative to DC is attenuated by greater than −3 dB (i.e. approximately halved) is given as:

$$f_{{\mathrm{cut}} - {\mathrm{off}}} = \frac{1}{{2{\mathrm{\pi \tau }}_{{\mathrm{membrane}}}}} = \frac{1}{{2{\mathrm{\pi}} RC}}$$

(3)

where here R and C are R_membrane and C_membrane, respectively. Above f_cut-off the voltage response decays with a frequency-dependent roll-off of −20 dB/decade. It can therefore be expected that an increase in either R_membrane and C_membrane will, by increasing the membrane time constant, lower the cut-off frequency of the filter effect.

Oscillatory spike locking

To find firing threshold, current-clamped neurons were first depolarized by bias current injection until repetitive action potential firing was observed. Holding potentials were then corrected to 5 mV hyperpolarized from this value at the start of each trial (V_Baseline¹⁰⁶). Sinusoidal current stimulation was repeated five times for each frequency. To avoid affecting AP initiation by fluctuating voltage level, cells in which V_Baseline changed by >2 mV were discarded. Phase-locking of AP firing was assessed by registering each AP’s peak time to the phase of the injected current waveform for each cycle, calculating using the Hilbert transform in MATLAB.

Data analysis

All data were analysed using custom-written routines in MATLAB, except for additional toolboxes detailed below. Where beneficial to improve analysis speed, C++ code was compiled and called from within MATLAB as MEX files.

Statistical analysis

Data is shown summarised as mean ± SEM unless otherwise stated. Individual parameter distributions were tested for normality (Kolmogorov–Smirnov tests, p > 0.05) before comparisons were made. Non-parametric tests were used under conditions of deviation from normality, as stated in the individual figure legends.

Spike train statistics

For spike probability analysis, spike times from cell-attached recordings were down-sampled to 1 ms resolution and converted to spike probability functions. In this approach each spike time is discretized at 1 ms and convolved with a Gaussian kernel with a standard deviation of 5 ms. Each spike’s resulting total integral was adjusted to 1 across a range of 3 standard deviations. Correspondingly, spike times are smoothed such that a spike occurs at time t contributes an extra ~0.1 spike probability to that of the next spike if it arises 30 ms later. The resulting spike density convolution is a continuous function rather than a discrete train; more amenable to comparison between trials and recordings, and less sensitive to binning artefacts than with a raw spike count histogram approach.

Thalamocortical integration simulations

NEURON code for simulations that repeat the main findings of these results and a graphical user interface for exploring interacting parameters in a feed-forward inhibitory circuit can be found at https://github.com/apfdomanski/Thalamocortical-Synaptic-Integration-in-Fmr1-KO-cortex. A single-compartment model neuron consisting of a soma with one dendrite was instantiated in NEURON¹¹⁴. Dendrite length and global leak conductance were tuned to match genotype mean values for input resistance and whole-cell capacitance of Fmr1-KO and wild-type recordings. To recreate TC and FFI inputs, two Exp2Syn model synapses were placed at the soma, with reversal potentials of 0 mV and −71 mV, respectively. Synaptic kinetics (rise and fall times) matched voltage-clamp data. The model TC input was activated with a peak conductance of 1 nS, yielding an ~60 pA peak inward current from a leak reversal potential of −60 mV. The model FFI input was then activated with an onset latency matching voltage-clamp recordings. To recreate the effects of FFI under conditions of different G/A ratios, peak conductance that was varied on successive trials relative to the peak TC conductance between 0 and 10×, in increments of 0.5.

Short-term plasticity model

Short-term depression dynamics of TC-EPSC and FF-IPSC current components were quantified using bi-exponential decay fits to trial-averaged voltage-clamp recordings, normalized to initial steady-state current amplitudes. Curves in Fig. 5c are best fits to within-cell measurements of G/A ratios, normalized to steady-state G/A ratio for each cell. Least-squares optimal fits to G/A ratio depression were obtained by using the fminsearch algorithm in MATLAB to obtain terms for underlying TC-EPSC and FF-IPSC components. For implementation in NEURON (v7.3), Depressing current inputs were re-fit with a phenomenological model of STP¹¹⁵. Briefly, starting from steady-state amplitude A₀, the stimulus-evoked response amplitude A of each current was multiplied by two dynamic depression factors, D₁ and D₂ (constrained <1):

$$A = A_0D_1D_2.$$

(4)

After each stimulus, D₁ and D₂ were multiplied by constants d₁ and d₂, representing the amount of depression per presynaptic action potential (i.e. D₁ → D₁d₁, and D₂ → D₂d₂). Between stimuli, D variables recovered exponentially back towards 1 with first-order kinetics governed by recovery time constants τ_D1 and τ_D2:

$${\mathrm{\tau }}_{{\mathrm{D}}1}\frac{{{\mathrm{d}}D_1}}{{{\mathrm{d}}t}} = 1 - D_1,\,{\mathrm{and}}\,{\mathrm{\tau }}_{{\mathrm{D}}2}\frac{{{\mathrm{d}}D_2}}{{{\mathrm{d}}t}} = 1 - D_2.$$

(5)

Accumulation of synaptic depression was therefore observed when the inter-stimulus interval was shorter than time required for recovery. Synaptic rise and decay time constants, as well as synaptic plasticity constants d₁, d₂, τ_D1 and τ_D2 were obtained from the genotype mean synaptic kinetics and short-term depression of EPSCs and FF-IPSCs individually. Fitting of short-term depression terms was performed on normalized synaptic depression curves for 5, 10, 20 and 50 Hz TC stimulation, equally weighted and optimised using a Levenberg-Marquardt search implemented in MATLAB.

Layer 4 network simulations

MATLAB code to reproduce the main findings of this paper is hosted at https://github.com/apfdomanski/Fmr1-KO-cortical-layer-4-spiking-network. A recurrently connected spiking network of single-compartment conductance-based leaky integrate-and-fire (LIF) neurons was instantiated in MATLAB. The subthreshold potential of the ith neuron evolved according to

$${\mathrm{\tau }}_{\mathrm{m}}\frac{{{\mathrm{d}}v_{\mathrm{i}}}}{{{\mathrm{d}}t}} = - \left[ {v_{\mathrm{i}}\left( t \right) - v_{{\mathrm{rest}}}} \right] + R_{\mathrm{m}}I_{{\mathrm{tot}}}^i(t)$$

(6)

where $I_{{\mathrm{tot}}}^i(t)$ is the summed synaptic current from recurrent and external inputs at time t. Spikes were emitted at $t^{\mathrm{f}} = \left\{ {t|v\left( t \right) = v_{{\mathrm{threhsold}}}} \right\}$ when v > v_threshold leading to the reset condition $v\left( {t + {\mathrm{d}}t} \right) \to v_{{\mathrm{reset}}} \, < \, v_{{\mathrm{threshold}}}$ for a finite refractory period of 1.5 ms.

$I_{{\mathrm{tot}}}^i(t)$ was calculated at each time step as:

$$I_{{\mathrm{tot}}}^i(t)= \mathop {\sum }\limits_{N_{\left( {i,{\mathrm{AMPArec}}} \right)}} I_{{\mathrm{AMPArec}}}^i\left( t \right) + \mathop {\sum }\limits_{N_{\left( {i,{\mathrm{NMDArec}}} \right)}} I_{{\mathrm{NMDArec}}}^i\left( t \right) \\ + \mathop {\sum }\limits_{N_{\left( {i,{\mathrm{GABArec}}} \right)}} I_{{\mathrm{GABArec}}}^i\left( t \right) + I_{{\mathrm{AMPArec}}}^i\left( t \right) + I_{{\mathrm{NMDArec}}}^i\left( t \right),$$

(7)

with N_{(i, AMPArec)}, N_{(i, NMDArec)} and N_{(i, GSBArec)} being the set of N recurrent synapses projecting onto the ith neuron, and $I_{{\mathrm{AMPArec}}}^{\mathrm{i}}\left( t \right)$, $I_{{\mathrm{NMDArec}}}^i\left( t \right)$, $I_{{\mathrm{GABArec}}}^i\left( t \right)$, $I_{{\mathrm{AMPAext}}}^i\left( t \right)$, $I_{{\mathrm{NMDAext}}}^i\left( t \right)$ the synaptic inputs from recurrent AMPA, NMDA and GABA and external AMPA and NMDA synapses, respectively.

Synaptic input current depended also on the postsynaptic membrane potential:

$$I_{{\mathrm{syn}}}^i\left( t \right) = A^i(t) \cdot \overline {g_{{\mathrm{syn}}}^i} \cdot s_{{\mathrm{syn}}}\left( t \right)\left( {v\left( t \right) - E_{{\mathrm{syn}}}} \right),$$

(8)

where Aⁱ(t) is the short-term depression state of synapse i described by Eq. (4), ${\bar g}_{{\mathrm{syn}}}$ and E_syn are, respectively, the maximal synaptic conductance and reversal potential of the synapse. s_syn(t) is a delayed-sum-of-exponentials (alpha) synapse model (after ref. ¹¹⁶):

$$S_{{\mathrm{syn}}}(t) = f\left( {{\mathrm{e}}^{ - \frac{{t_{{\mathrm{conduct}}} - \Delta t}}{{{\mathrm{\tau }}_{{\mathrm{decay}}}}}} - {\mathrm{e}}^{\frac{{t_{{\mathrm{conduct}}} - \Delta t}}{{{\mathrm{\tau }}_{{\mathrm{rise}}}}}}} \right){\cal{H}}(\Delta t),$$

(9)

in which t_conduct is a finite synaptic conduction delay, Δt > −∞ is the elapsed time t − t^f since the time of the last presynaptic spike. The Heaviside step function ${\cal{H}}$ ensures causality in time. The function f(x) is an amplitude normalisation factor:

$$f(x) = \frac{1}{{ - {\mathrm{e}}^{ - (t_{{\mathrm{peak}}} - t_{{\mathrm{conduct}}} - \Delta t)/{\mathrm{\tau }}_{{\mathrm{rise}}}} + {\mathrm{e}}^{ - (t_{{\mathrm{peak}}} - t_{{\mathrm{conduct}}} - \Delta t)/{\mathrm{\tau }}_{{\mathrm{rise}}}}}},$$

(10)

where the amplitude of the conductance is maximal at time t_peak:

$$t_{{\mathrm{peak}}} = \Delta t + t_{{\mathrm{conduct}}} + \frac{{{\mathrm{\tau }}_{{\mathrm{rise}}}{\mathrm{\tau }}_{{\mathrm{decay}}}}}{{{\mathrm{\tau }}_{{\mathrm{decay}}} - {\mathrm{\tau }}_{{\mathrm{rise}}}}}{\mathrm{ln}}\left( {\frac{{{\mathrm{\tau }}_{{\mathrm{decay}}}}}{{{\mathrm{\tau }}_{{\mathrm{rise}}}}}} \right),$$

(11)

NMDA conductances displayed additional voltage dependence:

$$I_{{\mathrm{NMDA}}}^{\mathrm{i}}\left( t \right) = {A^i\left( t \right)} \cdot G\left( {v\left( t \right)} \right) \cdot \overline {g^i} _{{\mathrm{NMDA}}} \cdot s_{{\mathrm{syn}}}\left( t \right)\left( {v\left( t \right) - E_{{\mathrm{NMDA}}/{\mathrm{AMPA}}}} \right),$$

(12)

In which G(v(t)) describes the voltage-dependent Mg²⁺-sensitive blockade of NMDA receptors:

$$G\left( {v\left( t \right)} \right) = \frac{1}{{1 + \frac{{{\mathrm{e}}^{ - av(t)}\left[ {{\mathrm{Mg}}^{2 + }} \right]_{{\mathrm{out}}}}}{b}}},$$

(13)

with a = 0.062 mV⁻¹ and b = 3.57 mM^117,118 and using our experimental extracellular Magnesium concentration [Mg²⁺]_out of 1.3 mM.

Network architecture

800 Excitatory and 150 Inhibitory neurons were randomly connected with probabilities obtained experimentally in Fig. 2, excluding autapses. All classes of synaptic connection were randomly assigned fractional weights up to ${\bar g}_{{\mathrm{syn}}}$, with the exception of Ex–Ex connections, which were log-normally distributed in line with our experimental findings (Fig. 3e) and previous reports^33,43,119. Simulated external thalamocortical inputs were provided to 80% of Ex and In neurons. ${\bar g}_{{\mathrm{syn}}}$ of external NMDA and AMPA were jittered between trials and thalamorecipient neurons were shuffled to introduce variability between trials. Short-term plasticity at TC synapses was modelled using the same fits to data as described above. Simulations (1000 ms) were carried out using forward Euler integration at 0.5 ms resolution. Spike times were discretized to 1 ms resolution. Modelling results show responses of 800 Ex neurons across N = 10 network simulation trials for each stimulation condition each from five random network seeds. For ensemble decoding analysis shown inf Fig. 10, results are calculated for 10 random permutations of each ensemble size, four oddball stimulus positions. Simulations were performed on five random network seeds, ten trials for each oddball position and frequency combination.

Parameter distributions

Cell-intrinsic, synaptic and connectivity parameters used are detailed in Supplemental Table 1. Where two numbers are shown, they represent mean and standard deviation of jittered variables. All parameters were assumed to vary independently, excluding input resistance and membrane capacitance, which displayed both private and shared variance.

Spike train metrics and ensemble decoding

Spike train distance metrics were computed with the method described by van Rossum⁶⁶ with a decay time constant of 50 ms. We used an optimised calculation of the van Rossum metric implemented in C++ by Houghton and Kreuz¹²⁰. Ensemble decoding was performed using a multivariate linear decoder using ensemble spike trains as covariates. To prevent over-fitting, pooled covariance matrices for regular and oddball conditions were regularised by a factor of 0.05 on the identity matrix. Decoder accuracy was evaluated with tenfold leave-one-out cross-validation, with decoders trained on nine trials asked to predict the stimulus condition that lead to the outcome of the remaining withheld trial. Training/test sets were split 70/30% and independent evaluation sets were used.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability

All raw data generated/analysed in this study are deposited as Supplementary Material on the Nature Publishing Group website.

Code availability

NEURON code for the single-cell model and MATLAB/C++ code for the Layer 4 network model are publicly available on Github via the above repository links. MATLAB/C++ code for analysis of data and simulations will be available from the authors upon reasonable request.

References

Marco, E. J., Hinkley, L. B. N., Hill, S. S. & Nagarajan, S. S. Sensory processing in autism: a review of neurophysiologic findings. Pediatr. Res. 69, 48R–54R (2011).
Article PubMed PubMed Central Google Scholar
Crane, L., Goddard, L. & Pring, L. Sensory processing in adults with autism spectrum disorders. Autism 13, 215–228 (2009).
Article PubMed Google Scholar
Rogers, S. J., Hepburn, S. & Wehner, E. Parent reports of sensory symptoms in toddlers with autism and those with other developmental disorders. J. Autism Dev. Disord. 33, 631–642 (2003).
Article PubMed Google Scholar
Hagerman, R. J. et al. Advances in the treatment of fragile X syndrome. Pediatrics 123, 378–390 (2009).
Article PubMed Google Scholar
Consortium, T. D. F. X. Fmr1 knockout mice: a model to study Fragile X mental retardation. Cell 78, 23–33 (1994).
Google Scholar
Bernardet, M. & Crusio, W. W. E. Fmr1 KO mice as a possible model of autistic features. Sci. World J. 6, 1164–1176 (2006).
Article CAS Google Scholar
Mineur, Y., Huynh, L. & Crusio, W. Social behavior deficits in the Fmr 1 mutant mouse. Behav. Brain. Res. 168, 172–175 (2006).
Article CAS PubMed Google Scholar
Shin Yim, Y. et al. Reversing behavioural abnormalities in mice exposed to maternal inflammation. Nature 549, 482–487 (2017).
Article ADS CAS PubMed Google Scholar
Gonçalves, J. T., Anstey, J. E., Golshani, P. & Portera-Cailliau, C. Circuit level defects in the developing neocortex of Fragile X mice. Nat. Neurosci. https://doi.org/10.1038/nn.3415 (2013).
Article PubMed PubMed Central Google Scholar
O’Donnell, C., Gonçalves, J. T., Sejnowski, T. J. & Portera-cailliau, C. Beyond excitation/inhibition imbalance in multidimensional models of neural circuit changes in brain disorders. eLife 6, 1–28 (2017).
Google Scholar
Zhang, Y. et al. Regulation of neuronal excitability by interaction of Fragile X mental retardation protein with slack potassium channels. J. Neurosci. 32, 15318–15327 (2012).
Article CAS PubMed PubMed Central Google Scholar
He, Q. et al. Critical period inhibition of NKCC1 rectifies synapse plasticity in the somatosensory cortex and restores adult tactile response maps in fragile X mice. Mol. Psychiatry 1, https://doi.org/10.1038/s41380-018-0048-y (2018).
He, C. X. et al. Tactile defensiveness and impaired adaptation of neuronal activity in the Fmr1 knock-out mouse model of autism. J. Neurosci. 37, 6475–6487 (2017).
Article CAS PubMed PubMed Central Google Scholar
Contractor, A., Klyachko, V. A. & Portera-Cailliau, C. Altered neuronal and circuit excitability in Fragile X syndrome. Neuron 87, 699–715 (2015).
Article CAS PubMed PubMed Central Google Scholar
Brown, V. et al. Microarray identification of FMRP-associated brain mRNAs and altered mRNA translational profiles in fragile X syndrome. Cell 107, 477–487 (2001).
Article CAS PubMed Google Scholar
Darnell, J. C. et al. Fragile X mental retardation protein targets G quartet mRNAs important for neuronal function. Cell 107, 489–499 (2001).
Article CAS PubMed Google Scholar
Darnell, J. C. et al. FMRP stalls ribosomal translocation on mRNAs linked to synaptic function and autism. Cell 146, 247–261 (2011).
Article CAS PubMed PubMed Central Google Scholar
Strumbos, J. G., Brown, M. R., Kronengold, J., Polley, D. B. & Kaczmarek, L. K. Fragile X mental retardation protein is required for rapid experience-dependent regulation of the potassium channel Kv3.1b. J. Neurosci. 30, 10263–10271 (2010).
Article CAS PubMed PubMed Central Google Scholar
Brown, M. R. et al. Fragile X mental retardation protein controls gating of the sodium-activated potassium channel Slack. Nat. Neurosci. 5–7, https://doi.org/10.1038/nn.2563 (2010).
Deng, P.-Y. et al. FMRP regulates neurotransmitter release and synaptic information transmission by modulating action potential duration via BK channels. Neuron 77, 696–711 (2013).
Article CAS PubMed PubMed Central Google Scholar
Rubenstein, J. L. R. & Merzenich, M. M. Model of autism: increased ratio of excitation/inhibition in key neural systems. Genes, Brain Dev. 255–267, https://doi.org/10.1046/j.1601-183x.2003.00037.x (2003).
Nelson, S. B. & Valakh, V. Excitatory/inhibitory balance and circuit homeostasis in autism spectrum disorders. Neuron 87, 684–698 (2015).
Article CAS PubMed PubMed Central Google Scholar
He, Q., Nomura, T., Xu, J. & Contractor, A. The developmental switch in GABA polarity is delayed in fragile X mice. J. Neurosci. 34, 446–450 (2014).
Article CAS PubMed PubMed Central Google Scholar
Gibson, J. R., Bartley, A. F. A., Hays, S. A. & Huber, K. M. Imbalance of neocortical excitation and inhibition and altered UP states reflect network hyperexcitability in the mouse model of fragile X syndrome. J. Neurophysiol. 100, 2615–2626 (2008).
Article PubMed PubMed Central Google Scholar
Paluszkiewicz, S. M., Olmos-Serrano, J. L., Corbin, J. G. & Huntsman, M. M. Impaired inhibitory control of cortical synchronization in fragile x syndrome. J. Neurophysiol. https://doi.org/10.1152/jn.00421.2011 (2011).
Article PubMed PubMed Central Google Scholar
Vislay, R. L. et al. Homeostatic responses fail to correct defective amygdala inhibitory circuit maturation in fragile X syndrome. J. Neurosci. 33, 7548–7558 (2013).
Article CAS PubMed PubMed Central Google Scholar
Cea-Del Rio, Ca & Huntsman, M. M. The contribution of inhibitory interneurons to circuit dysfunction in Fragile X Syndrome. Front. Cell. Neurosci. 8, 1–7 (2014).
Article CAS Google Scholar
Harlow, E. et al. Critical period plasticity is disrupted in the barrel cortex of Fmr1 knockout mice. Neuron 65, 385–398 (2010).
Article CAS PubMed PubMed Central Google Scholar
Daw, M. I., Scott, H. L. & Isaac, J. T. R. Developmental synaptic plasticity at the thalamocortical input to barrel cortex: mechanisms and roles. Mol. Cell. Neurosci. 34, 493–502 (2007).
Article CAS PubMed PubMed Central Google Scholar
Daw, M. I., Bannister, N. V. & Isaac, J. T. R. Rapid, activity-dependent plasticity in timing precision in neonatal barrel cortex. J. Neurosci. 26, 4178–4187 (2006).
Article CAS PubMed PubMed Central Google Scholar
Daw, M. I., Ashby, M. C. & Isaac, J. T. R. Coordinated developmental recruitment of latent fast spiking interneurons in layer IV barrel cortex. Nat. Neurosci. 10, 453–461 (2007).
Article CAS PubMed Google Scholar
Chittajallu, R. & Isaac, J. T. R. Emergence of cortical inhibition by coordinated sensory-driven plasticity at distinct synaptic loci. Nat. Neurosci. 13, 1240–1248 (2010).
Article CAS PubMed PubMed Central Google Scholar
Ashby, M. C. & Isaac, J. T. R. Maturation of a recurrent excitatory neocortical circuit by experience-dependent unsilencing of newly formed dendritic spines. Neuron 70, 510–521 (2011).
Article CAS PubMed PubMed Central Google Scholar
Landers, M. & Philip Zeigler, H. Development of rodent whisking: Trigeminal input and central pattern generation. Somatosens. Mot. Res. 23, 1–10 (2006).
Article PubMed Google Scholar
Paluszkiewicz, S. M., Martin, B. S. & Huntsman, M. M. Fragile X syndrome: the GABAergic system and circuit dysfunction. Dev. Neurosci. 20010, 1–16 (2011).
Google Scholar
Baudouin, S. J. et al. Shared synaptic pathophysiology in syndromic and nonsyndromic rodent models of autism. Science 338, 128–132 (2012).
Cellot, G. & Cherubini, E. Reduced inhibitory gate in the barrel cortex of Neuroligin3R451C knock-in mice, an animal model of autism spectrum disorders. Physiol. Rep. 2, 1–12 (2014).
Article CAS Google Scholar
Miller, K. D., Pinto, D. J. & Simons, D. J. Processing in layer 4 of the neocortical circuit: new insights from visual and somatosensory cortex. Curr. Opin. Neurobiol. 11, 488–497 (2001).
Article CAS PubMed Google Scholar
Favorov, O. V. & Kursun, O. Neocortical layer 4 as a pluripotent function linearizer. J. Neurophysiol. 105, 1342–1360 (2011).
Article PubMed Google Scholar
Hays, S. A., Huber, K. M. & Gibson, J. R. Altered neocortical rhythmic activity states in Fmr1 KO mice are due to enhanced mGluR5 signaling and involve changes in excitatory circuitry. J. Neurosci. 31, 14223–14234 (2011).
Article CAS PubMed PubMed Central Google Scholar
Agmon, A. & Connors, B. W. Thalamocortical responses of mouse somatosensory (barrel) cortex in vitro. Neuroscience 41, 365–379 (1991).
Article CAS PubMed Google Scholar
Feldman, D. E., Nicoll, Ra & Malenka, R. C. Synaptic plasticity at thalamocortical synapses in developing rat somatosensory cortex: LTP, LTD, and silent synapses. J. Neurobiol. 41, 92–101 (1999).
Article CAS PubMed Google Scholar
Lefort, S., Tomm, C., Floyd Sarria, J.-C., Petersen, C. C. H. & Sarria, J. F. The excitatory neuronal network of the C2 barrel column in mouse primary somatosensory cortex. Neuron 61, 301–316 (2009).
Article CAS PubMed Google Scholar
Feldman, D. E. Timing-based LTP and LTD at vertical inputs to layer II/III pyramidal cells in rat barrel cortex. Neuron 27, 45–56 (2000).
Article CAS PubMed Google Scholar
Feldmeyer, D. & Sakmann, B. Synaptic connections between layer 4 spiny neurone-layer 2/3 pyramidal cell pairs in juvenile rat barrel cortex: physiology and anatomy of interlaminar signalling within a cortical column. Cell 803–822, https://doi.org/10.1013/jphysiol.2001.012959 (2002).
Feldmeyer, D. et al. Barrel cortex function. Prog. Neurobiol. 1–25, https://doi.org/10.1016/j.pneurobio.2012.11.002 (2012).
Helmstaedter, M., Staiger, J. F., Sakmann, B. & Feldmeyer, D. Efficient recruitment of layer 2/3 interneurons by layer 4 input in single columns of rat somatosensory cortex. J. Neurosci. 28, 8273–8284 (2008).
Article CAS PubMed PubMed Central Google Scholar
Gabernet, L., Jadhav, S. P., Feldman, D. E., Carandini, M. & Scanziani, M. Somatosensory integration controlled by dynamic thalamocortical feed-forward inhibition. Neuron 48, 315–327 (2005).
Article CAS PubMed Google Scholar
Swadlow, H. A. Thalamocortical control of feed-forward inhibition in awake somatosensory’barrel’cortex. Philos. Trans. R. Soc. B Biol. Sci. 357, 1717 (2002).
Article Google Scholar
Koch, C. P. Biophysics of Computation. (Oxford University Press, 1999).
Crochet, S., Poulet, J. F. A., Kremer, Y. & Petersen, C. C. H. Synaptic mechanisms underlying sparse coding of active touch. Neuron 69, 1160–1175 (2011).
Article CAS PubMed Google Scholar
Jadhav, S. P., Wolfe, J. & Feldman, D. E. Sparse temporal coding of elementary tactile features during active whisker sensation. Nat. Neurosci. 12, 792–800 (2009).
Article CAS PubMed Google Scholar
Ferezou, I. et al. Spatiotemporal dynamics of cortical sensorimotor integration in behaving mice. Neuron 56, 907–923 (2007).
Article CAS PubMed Google Scholar
O’Connor, D. H., Peron, S. P., Huber, D. & Svoboda, K. Neural activity in barrel cortex underlying vibrissa-based object localization in mice. Neuron 67, 1048–1061 (2010).
Article CAS PubMed Google Scholar
Xu, N. et al. Nonlinear dendritic integration of sensory and motor input during an active sensing task. Nature. https://doi.org/10.1038/nature11601 (2012).
Article PubMed PubMed Central Google Scholar
Lawrence, J. J., Statland, J. M., Grinspan, Z. M. & McBain, C. J. Cell type-specific dependence of muscarinic signalling in mouse hippocampal stratum oriens interneurones. J. Physiol. 570, 595–610 (2006).
Article CAS PubMed Google Scholar
Erchova, I., Kreck, G., Heinemann, U., Herz, A. V. M. & Johannes, M. Dynamics of rat entorhinal cortex layer II and III cells: characteristics of membrane potential resonance at rest predict oscillation properties near threshold. Society 1, 89–110 (2004).
Google Scholar
Cruikshank, S. J., Urabe, H., Nurmikko, A. V & Connors, B. W. Pathway-specific feedforward circuits between thalamus and neocortex revealed by selective optical stimulation of axons. Neuron 230–245, https://doi.org/10.1016/j.neuron.2009.12.025 (2010).
Pouille, F. et al. Enforcement of temporal fidelity in pyramidal cells by somatic feed-forward inhibition. Science 293, 1159–1163 (2001).
Article CAS PubMed Google Scholar
Swadlow, H. A. Sharp, local synchrony among putative feed-forward inhibitory interneurons of rabbit somatosensory cortex. J. Neurophysiol. 79, 567–582 (1998).
Abbott, L. F. & Regehr, W. G. Synaptic computation. Nature 431, 796–803 (2004).
Article ADS CAS PubMed Google Scholar
Silver, R. A. Neuronal arithmetic. Nat. Rev. Neurosci. 11, 474–489 (2010).
Article CAS PubMed PubMed Central Google Scholar
Petersen, C. C. H. & Loos, V. Der. The functional organization of the barrel cortex. October 339–355, https://doi.org/10.1016/j.neuron.2007.09.017 (2007).
Booker, S. A. et al. Altered dendritic spine function and integration in a mouse model of fragile X Syndrome. bioRxiv, 396986 (2018).
Panzeri, S., Petersen, R. S., Schultz, S. R., Lebedev, M. & Diamond, M. E. The role of spike timing in the coding of stimulus location in rat somatosensory cortex. Neuron 29, 769–777 (2001).
Article CAS PubMed MATH Google Scholar
van Rossum, M. C. W. A novel spike distance. Neural Comput. 13, 751–763 (2001).
Article PubMed MATH Google Scholar
Victor, J. D. Spike train metrics. Curr. Opin. Neurobiol. 15, 585–592 (2005).
Article CAS PubMed PubMed Central Google Scholar
Petersen, R. S., Panzeri, S. & Diamond, M. E. Population coding of stimulus location in rat somatosensory cortex. Neuron 32, 503–514 (2001).
Article CAS PubMed Google Scholar
Olshausen, Ba & Field, D. J. Sparse coding of sensory inputs. Curr. Opin. Neurobiol. 14, 481–487 (2004).
Article CAS PubMed Google Scholar
Golshani, P. et al. Internally mediated developmental desynchronization of neocortical network activity. J. Neurosci. 29, 10890–10899 (2009).
Article CAS PubMed PubMed Central Google Scholar
Olshausen, B. A. & Field, D. J. Sparse coding with an overcomplete basis set: a strategy employed by V1? Vis. Res. 37, 3311–3325 (1997).
Article CAS PubMed Google Scholar
Celikel, T., Szostak, V. A. & Feldman, D. E. Modulation of spike timing by sensory deprivation during induction of cortical map plasticity. Nat. Neurosci. 7, 534–541 (2004).
Article CAS PubMed PubMed Central Google Scholar
Bishop, C. M. Pattern recognition and Machine Learning. (Springer, 2011).
Cascio, C. J. Somatosensory processing in neurodevelopmental disorders. J. Neurodev. Disord. 2, 62–69 (2010).
Article PubMed PubMed Central Google Scholar
Chao, H.-T. et al. Dysfunction in GABA signalling mediates autism-like stereotypies and Rett syndrome phenotypes. Nature 468, 263–269 (2010).
Article ADS CAS PubMed PubMed Central Google Scholar
Yizhar, O. et al. Neocortical excitation/inhibition balance in information processing and social dysfunction. Nature 1–8, https://doi.org/10.1038/nature10360 (2011).
Uhlhaas, P. J. Dysconnectivity, large-scale networks and neuronal dynamics in schizophrenia. Curr. Opin. Neurobiol. 23, 283–290 (2013).
Article CAS PubMed Google Scholar
Lisman, J. Excitation, inhibition, local oscillations, or large-scale loops: what causes the symptoms of schizophrenia? Curr. Opin. Neurobiol. 22, 537–544 (2012).
Article CAS PubMed Google Scholar
Haider, B. & Mccormick, D. A. Rapid neocortical dynamics: cellular and network mechanisms. Neuron 62, 171–189 (2009).
Article CAS PubMed PubMed Central Google Scholar
Wahlstrom-Helgren, S. & Klyachko, V. A. Dynamic balance of excitation and inhibition rapidly modulates spike probability and precision in feed-forward hippocampal circuits. J. Neurophysiol. 116, 2564–2575 (2016).
Article PubMed PubMed Central Google Scholar
Gentet, L. J., Avermann, M., Matyas, F., Staiger, J. F. & Petersen, C. C. H. Membrane potential dynamics of GABAergic neurons in the barrel cortex of behaving mice. Neuron 65, 422–435 (2010).
Article CAS PubMed Google Scholar
Haider, B., Häusser, M. & Carandini, M. Inhibition dominates sensory responses in the awake cortex. Nature 2–12, https://doi.org/10.1038/nature11665 (2012).
Hasenstaub, A. R. et al. Inhibitory postsynaptic potentials carry synchronized frequency information in active cortical networks. Neuron 47, 423–435 (2005).
Article CAS PubMed Google Scholar
Antoine, M. W., Langberg, T., Schnepel, P. & Feldman, D. E. Increased excitation-inhibition ratio stabilizes synapse and circuit excitability in four autism mouse models. Neuron 101, 648–661.e4 (2019).
Article CAS PubMed PubMed Central Google Scholar
Meredith, R. M. & Mansvelder, H. D. STDP and mental retardation: dysregulation of dendritic excitability in fragile X syndrome. Front. Synaptic Neurosci. 2, 10 (2010).
PubMed PubMed Central Google Scholar
Brager, D. H., Akhavan, A. R. & Johnston, D. Impaired dendritic expression and plasticity of h-channels in the fmr1(-/y) mouse model of fragile X syndrome. Cell Rep. 1, 225–233 (2012).
Article CAS PubMed PubMed Central Google Scholar
Olmos-Serrano, J. L. et al. Defective GABAergic neurotransmission and pharmacological rescue of neuronal hyperexcitability in the amygdala in a mouse model of fragile X syndrome. J. Neurosci. 30, 9929–9938 (2010).
Article CAS PubMed PubMed Central Google Scholar
Deng, P.-Y. & Klyachko, V. A. Increased persistent sodium current causes neuronal hyperexcitability in the entorhinal cortex of Fmr1 knockout mice. Cell Rep. 16, 3157–3166 (2016).
Article CAS PubMed PubMed Central Google Scholar
Deng, P.-Y. et al. Voltage-independent SK-channel dysfunction causes neuronal hyperexcitability in the hippocampus of Fmr1 knock-out mice. J. Neurosci. 39, 28–43 (2019).
Article CAS PubMed PubMed Central Google Scholar
Routh, B. N., Johnston, D. & Brager, D. H. Loss of functional A-type potassium channels in the dendrites of CA1 pyramidal neurons from a mouse model of fragile X syndrome. J. Neurosci. 33, 19442–19450 (2013).
Article CAS PubMed PubMed Central Google Scholar
Dani, V. S. et al. Reduced cortical activity due to a shift in the balance between excitation and inhibition in a mouse model of Rett syndrome. Proc. Natl Acad. Sci. USA 102, 12560–12565 (2005).
Article ADS CAS PubMed PubMed Central Google Scholar
Lo, F.-S., Blue, M. E. & Erzurumlu, R. S. Enhancement of postsynaptic GABA_A and extrasynaptic NMDA receptor-mediated responses in the barrel cortex of Mecp2-null mice. J. Neurophysiol. 115, 1298–1306 (2016).
Article PubMed Google Scholar
Egger, V., Feldmeyer, D. & Sakmann, B. Coincidence detection and changes of synaptic efficacy in spiny stellate neurons in rat barrel cortex. Nat. Neurosci. 2, 1098–1105 (1999).
Article CAS PubMed Google Scholar
Meredith, R. M., Holmgren, C. & Weidum, M. Increased threshold for spike-timing-dependent plasticity is caused by unreliable calcium signaling in mice lacking fragile X gene Fmr1. Neuron 627–638, https://doi.org/10.1016/j.neuron.2007.04.028 (2007).
Pacey, L. K. K., Heximer, S. P. & Hampson, D. R. Increased GABAB receptor-mediated signaling reduces the susceptibility of fragile X knockout mice to audiogenic seizures. Mol. Pharmacol. 76, 18–24 (2009).
Article CAS PubMed Google Scholar
Henderson, C. et al. Reversal of disease-related pathologies in the fragile X mouse model by selective activation of GABAB receptors with arbaclofen. Sci. Transl. Med. 4, 152ra128–152ra128 (2012).
Article CAS PubMed PubMed Central Google Scholar
Sinclair, D. et al. GABA-B agonist baclofen normalizes auditory-evoked neural oscillations and behavioral deficits in the Fmr1 knockout mouse model of fragile X syndrome. Eneuro 4, ENEURO.0380-16.2017 (2017).
Lo, F.-S. & Erzurumlu, R. S. Insulin receptor sensitization restores neocortical excitation/inhibition balance in a mouse model of autism. Mol. Autism 9, 13 (2018).
Article CAS PubMed PubMed Central Google Scholar
Lo, F.-S., Erzurumlu, R. S. & Powell, E. M. Insulin-independent GABA_A receptor-mediated response in the barrel cortex of mice with impaired met activity. J. Neurosci. 36, 3691–3697 (2016).
Article CAS PubMed PubMed Central Google Scholar
Bureau, I., Shepherd, G. M. G. & Svoboda, K. Circuit and plasticity defects in the developing somatosensory cortex of FMR1 knock-out mice. J. Neurosci. 28, 5178–5188 (2008).
Article CAS PubMed PubMed Central Google Scholar
Till, S. M. et al. Altered maturation of the primary somatosensory cortex in a mouse model of fragile X syndrome. Hum. Mol. Genet. 1–14, https://doi.org/10.1093/hmg/dds030 (2012).
Meredith, R. M., Dawitz, J. & Kramvis, I. Sensitive time-windows for susceptibility in neurodevelopmental disorders. Trends Neurosci. 35, 335–344 (2012).
Article CAS PubMed Google Scholar
Gainey, M. A., Aman, J. W. & Feldman, D. E. Rapid disinhibition by adjustment of PV intrinsic excitability during whisker map plasticity in mouse S1. J. Neurosci. 38, 4749–4761 (2018).
Article CAS PubMed PubMed Central Google Scholar
Suter, B. Ephus: multipurpose data acquisition software for neuroscience experiments. Front. Neural Circuits 4, 1–12 (2010).
Article Google Scholar
Perkins, K. L. Cell-attached voltage-clamp and current-clamp recording and stimulation techniques in brain slices. J. Neurosci. Methods 154, 1–18 (2006).
Article CAS PubMed PubMed Central Google Scholar
Lawrence, J. J., Grinspan, Z. M., Statland, J. M. & McBain, C. J. Muscarinic receptor activation tunes mouse stratum oriens interneurones to amplify spike reliability. J. Physiol. 571, 555–562 (2006).
Article CAS PubMed PubMed Central Google Scholar
Alcami, P., Franconville, R., Llano, I. & Marty, A. Measuring the firing rate of high-resistance neurons with cell-attached recording. J. Neurosci. 32, 3118–3130 (2012).
Article CAS PubMed PubMed Central Google Scholar
Li, P. & Huntsman, M. M. Two functional inhibitory circuits are comprised of a heterogeneous population of fast-spiking cortical interneurons. Neuroscience 265, 60–71 (2014).
Article CAS PubMed Google Scholar
Cole, K. S. & Curtis, H. J. Electrical impedance of the squid giant axon during activity. Cold Spring Harb. Symp. Quant. Biol. 22, 649–670 (1936).
Cole, K. S. & Curtis, H. J. Electrical impedance of nerve during activity. Nature 142, 209 (1938).
Carandini, M., Mechler, F., Leonard, C. S. & Movshon, Ja Spike train encoding by regular-spiking cells of the visual cortex. J. Neurophysiol. 76, 3425–3441 (1996).
Article CAS PubMed Google Scholar
Puil, E., Gimbarzevsky, B. & Miura, R. M. Quantification of membrane properties of trigeminal root ganglion neurons in guinea pigs. J. Neurophysiol. 55, 995–1016 (1986).
Article CAS PubMed Google Scholar
Zemankovics, R., Szabolcs, K., Paulsen, O., Freund, F. & Norbert, H. Differences in subthreshold resonance of hippocampal pyramidal cells and interneurons: the role of h-current and passive membrane characteristics. J. Physiol. 12, 2109–2132 (2011).
Google Scholar
Hines, M. L. et al. The NEURON simulation environment. Neuron 23, 349–398 (2002).
Google Scholar
Varela, J. A. et al. A quantitative description of short-term plasticity at excitatory synapses in layer 2/3 of rat primary visual cortex. J. Neurosci. 17, 7926–7940 (1997).
Article CAS PubMed PubMed Central Google Scholar
Dayan, P. & Abbott, L. F. Theoretical Neuroscience: Computational and Mathematical Modeling of Neural Systems (Computational Neuroscience Series). (2000) MIT Press.
CE, J. & Stevens, C. F. Voltage dependence of NMDA-activated macroscopic conductances predicted by single-channel kinetics. J. Neurosci. 10, 3178–3182 (1990).
Article Google Scholar
CE, J. & Stevens, C. F. A quantitative description of NMDA receptor-channel kinetic behavior. J. Neurosci. 10, 1830–1837 (1990).
Article Google Scholar
Petersen, C. C. H. & Sakmann, B. The excitatory neuronal network of rat layer 4 barrel cortex. J. Neurosci. 20, 7579–7586 (2000).
Article CAS PubMed PubMed Central Google Scholar
Houghton, C. & Kreuz, T. On the efficient calculation of van Rossum distances. Netw. Comput. Neural Syst. 23, 48–58 (2012).
Article Google Scholar

Download references

Acknowledgements

The authors wish to thank Ramesh Chittajallu, Chris McBain, Michael Ashby, Michael Daw, Timothy o’Leary, Mark van Rossum and past/current members of the Kind, McBain, Isaac and Wyllie laboratories for invaluable scientific input, and Owen Dando and Zrinko Kozic for statistical discussions. Funders: Simons Foundation Autism Research Initiative (529085), The Patrick Wild Centre, Medical Research Council UK (MR/P006213/1), The Shirley Foundation and the RS Macdonald Charitable Trust, Wellcome-National Institutes of Health Collaborative Scholarship (APFD).

Author information

John T. R. Isaac
Present address: Janssen Neuroscience, J&J London Innovation Centre, J&J London Innovation Centre, One Chapel Place, London, W1G 0B, UK

Authors and Affiliations

School of Physiology, Pharmacology & Neuroscience, University of Bristol, Bristol, UK
Aleksander P. F. Domanski
Centre for Discovery Brain Sciences, University of Edinburgh, Hugh Robson Building, George Square, Edinburgh, EH8 9XD, UK
Aleksander P. F. Domanski, Sam A. Booker, David J. A. Wyllie & Peter C. Kind
Patrick Wild Centre, University of Edinburgh, Hugh Robson Building, George Square, Edinburgh, EH8 9XD, UK
Aleksander P. F. Domanski, Sam A. Booker, David J. A. Wyllie & Peter C. Kind
Developmental Synaptic Plasticity Section, NINDS, NIH, Bethesda, MD, 20892, USA
Aleksander P. F. Domanski & John T. R. Isaac
Simons Initiative for the Developing Brain, University of Edinburgh, Hugh Robson Building, George Square, Edinburgh, EH8 9XD, UK
Sam A. Booker, David J. A. Wyllie & Peter C. Kind
Centre for Brain Development and Repair, NCBS, GKVK Campus, Bangalore, 560065, India
David J. A. Wyllie & Peter C. Kind

Authors

Aleksander P. F. Domanski
View author publications
You can also search for this author in PubMed Google Scholar
Sam A. Booker
View author publications
You can also search for this author in PubMed Google Scholar
David J. A. Wyllie
View author publications
You can also search for this author in PubMed Google Scholar
John T. R. Isaac
View author publications
You can also search for this author in PubMed Google Scholar
Peter C. Kind
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

A.P.F.D.: designed and interpreted, performed experiments, developed and deployed simulations, analysed data and wrote the paper; S.A.B.: designed and performed experiments, analysed/interpreted data; D.J.A.W.: designed experiments, analysed/interpreted data; J.T.R.I.: designed experiments, analysed/interpreted data, obtained funding and wrote the paper P.C.K.: designed experiments, analysed/interpreted data, obtained funding and wrote the paper.

Corresponding authors

Correspondence to Aleksander P. F. Domanski, John T. R. Isaac or Peter C. Kind.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information: Nature Communications thanks Reha Erzurumlu, Molly Huntsman and Vitaly Klyachko for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Supplementary Information

Peer Review File

Reporting Summary

Source data

Source Data

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Domanski, A.P.F., Booker, S.A., Wyllie, D.J.A. et al. Cellular and synaptic phenotypes lead to disrupted information processing in Fmr1-KO mouse layer 4 barrel cortex. Nat Commun 10, 4814 (2019). https://doi.org/10.1038/s41467-019-12736-y

Download citation

Received: 15 November 2018
Accepted: 23 September 2019
Published: 23 October 2019
DOI: https://doi.org/10.1038/s41467-019-12736-y

This article is cited by

Identifying foetal forebrain interneurons as a target for monogenic autism risk factors and the polygenic 16p11.2 microdeletion
- Yifei Yang
- Sam A. Booker
- Thomas Pratt
BMC Neuroscience (2023)
Fmr1-KO mice failure to detect object novelty associates with a post-test decrease of structural and synaptic plasticity upstream of the hippocampus
- Antonella Borreca
- Mariassunta De Luca
- Martine Ammassari-Teule
Scientific Reports (2023)
Elevated levels of FMRP-target MAP1B impair human and mouse neuronal development and mouse social behaviors via autophagy pathway
- Yu Guo
- Minjie Shen
- Xinyu Zhao
Nature Communications (2023)
Experience-dependent changes in hippocampal spatial activity and hippocampal circuit function are disrupted in a rat model of Fragile X Syndrome
- Antonis Asiminas
- Sam A. Booker
- Emma R. Wood
Molecular Autism (2022)
Reduced expression of the psychiatric risk gene DLG2 (PSD93) impairs hippocampal synaptic integration and plasticity
- Simonas Griesius
- Cian O’Donnell
- Jack R. Mellor
Neuropsychopharmacology (2022)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.