Hyperactivated PTP1B phosphatase in parvalbumin neurons alters anterior cingulate inhibitory circuits and induces autism-like behaviors

Individuals with autism spectrum disorder (ASD) have social interaction deficits and difficulty filtering information. Inhibitory interneurons filter information at pyramidal neurons of the anterior cingulate cortex (ACC), an integration hub for higher-order thalamic inputs important for social interaction. Humans with deletions including LMO4, an endogenous inhibitor of PTP1B, display intellectual disabilities and occasionally autism. PV-Lmo4KO mice ablate Lmo4 in PV interneurons and display ASD-like repetitive behaviors and social interaction deficits. Surprisingly, increased PV neuron-mediated peri-somatic feedforward inhibition to the pyramidal neurons causes a compensatory reduction in (somatostatin neuron-mediated) dendritic inhibition. These homeostatic changes increase filtering of mediodorsal-thalamocortical inputs but reduce filtering of cortico-cortical inputs and narrow the range of stimuli ACC pyramidal neurons can distinguish. Simultaneous ablation of PTP1B in PV-Lmo4KO neurons prevents these deficits, indicating that PTP1B activation in PV interneurons contributes to ASD-like characteristics and homeostatic maladaptation of inhibitory circuits may contribute to deficient information filtering in ASD.

I ndividuals with autism spectrum disorder (ASD) have impaired social communication, inflexible and repetitive stereotyped behaviors, and cognitive deficits tied to deficient filtering of task-irrelevant information 1 . Inhibitory GABAergic interneurons are key to filtering information in the cortex. While mutations in >100 genes that affect synapse/network formation have been tied to ASD 2,3 , targeting ablation of ASD-risk genes only to fast-spiking parvalbumin (PV) interneurons is sufficient to produce similar ASD-like phenotypes as in mice with global ablation 4 . Mice lacking PV interneurons display ASD-like behaviors 5 . In the prefrontal cortex (PFC) of ASD subjects postmortem 6 and in some ASD mouse models 7 , there is a selective loss of PV, but not other GABAergic interneurons. Prenatal exposure to the anti-epileptic valproic acid (VPA) that increases the risk of ASD~6-10 fold 8 selectively affects PV interneurons 9 . Thus, altered function of PV interneurons affecting the balance of excitatory and inhibitory (E/I) transmission has emerged as an important mechanism underlying ASD.
The medial prefrontal cortex (mPFC) performs higher cognitive processes, including decision making, attention, working memory, emotional control, and social interaction 10 . The mPFC (particularly, the dorsal anterior cingulate cortex, dACC) makes reciprocal connections with the mediodorsal (MD) thalamus 11 . Imaging studies reveal altered function at the anterior cingulate cortex (ACC) and aberrant connectivity between the MD thalamus and ACC in ASD subjects during task performance 12,13 . Hypofunction or lesions at these two regions reduce social interaction in animals [14][15][16] . An elegant recent study confirmed the relevance of the ACC to social interaction behaviors; selective ablation of the autism-risk gene Shank3 in ACC pyramidal neurons caused social interaction deficits while restoration of Shank3 only in ACC neurons rescued social deficits in global Shank3 mutant mice 17 .
Unlike other regions of the thalamus that receive and relay sensory inputs to the cortex, the MD thalamus is a "higher-order" thalamic nucleus that integrates "already processed" information from the layer 5 cortex and relays it back to the layer 2/3 (L2/3) ACC 11 . MD thalamocortical projections excite L2/3 glutamatergic pyramidal neurons and simultaneously activate GABAergic PV interneurons that synapse onto the peri-somata of L2/3 pyramidal neurons and provide feedforward inhibition (FFI) 18 . The fastspiking PV-mediated FFI provides a temporal filter to limit the "window of opportunity" during which pyramidal neurons integrate excitatory inputs 19 . In addition to peri-somatic FFI, pyramidal neurons also receive dendritic inhibition from somatostatin (SST) interneurons 20 . Together, these inhibitory inputs enable the recruitment of a population of pyramidal neurons in a progressive manner over a wide dynamic range of afferent input strengths; without them, cortical pyramidal neurons would be recruited in an all-or-none fashion that would limit their ability to respond to afferent inputs of different intensities 21,22 . These mechanisms are crucial for information processing and the execution of complex tasks, including social interaction. However, it remains uncertain to what extent deficits in peri-somatic and/or dendritic inhibition contribute to ASD.
Rare single allele deletions containing LMO4 in humans are related to several cases of intellectual disability and one case of autism (Decipher database, https://decipher.sanger.ac.uk/). LMO4 expression is reduced in lymphoblasts of patients with autism 23 and in human cells that carry a mutation in MeCP2 that causes Rett syndrome 24 , another neurodevelopmental disorder with ASD-like behavior deficits. Of note, LMO4 is an endogenous inhibitor of the tyrosine phosphatase PTP1B 25,26 that was recently implicated in Rett syndrome 27 . Inhibition of PTP1B eliminated repetitive behaviors and improved motor function in MeCP2-deficient mice 27 , suggesting that a severe neurodevelopmental disorder could be ameliorated by postnatal pharmacological intervention targeting PTP1B.
To determine how unopposed PTP1B function in PV interneurons affects cortical function, we examined local-and longrange circuits of the ACC in mice with PV interneuron-specific deletion of the endogenous PTP1B inhibitor LMO4 25,26 . We found that PV-Lmo4KO mice display ASD-like behaviors and have markedly increased PV interneuron excitability that heightens thalamocortical FFI onto the peri-somata of L2/3 ACC pyramidal neurons. Surprisingly, dendritic inhibitory inputs to these pyramidal neurons from non-PV interneurons was reciprocally reduced. Importantly, simultaneous ablation of PTP1B and Lmo4 in PV interneurons was sufficient to restore these inhibitory synaptic circuits and to prevent ASD-like behaviors.

Results
PV-Lmo4KO mice display ASD-like behaviors. Lmo4 is expressed in the medial and lateral ganglionic eminences that give rise to cortical interneurons, including PV neurons 28 . To avoid the embryonic lethality of global Lmo4 deletion 29 , we selectively knocked out Lmo4 in PV neurons (PV-Lmo4KO). PV-Lmo4KO mice display reduced social interaction (Fig. 1a) and repetitive stereotyped behaviors (enhanced grooming, enhanced digging) (Fig. 1b, c). Mutant mice are not anxious and display normal olfactory responses to both social cues and non-social cues (Supplementary Fig. 1), suggesting that impaired social interaction probably did not involve olfactory dysfunction. While social interaction increased the number of c-Fos+ neurons (an index of neuronal activation) in layer 2/3 of the dACC nearly eightfold compared to naive wild-type mice, fewer c-Fos+ neurons were counted in PV-Lmo4KO mice (Fig. 1d). Despite these deficits, PV-Lmo4KO mice show no loss of PV expression or PV neuron numbers ( Supplementary Fig. 2), suggesting that loss of Lmo4 in PV neurons affects other neuronal properties leading to ASD-like behaviors.
Increased excitability of PV interneurons. To determine the functional consequence of Lmo4 ablation in PV interneurons at the dACC, we measured their electrophysiological properties. AAV9 vectors expressing Cre-dependent mCherry were injected to the dACC of PV-Cre/Lmo4 WT/WT (WT) or PV-Cre/Lmo4 flox/flox (PV-Lmo4KO) mice to label PV interneurons (Fig. 2a). PV interneurons in L2/3 of the dACC showed a relatively depolarized resting membrane potential (RMP) (Fig. 2b) with increased membrane excitability (Fig. 2c) in PV-Lmo4KO mice. Importantly, the latency to the first action potential (AP) was markedly shortened (~1.45 ms) (Fig. 2d-f); many features could account for this change, including reduced membrane conductance and/or reduced I A currents (e.g., Kv1.4 (Kcna4), Kv4.2 (Kcnd2) or Kv4.3 (Kcnd3)). However, further examination revealed no deficit in I A currents in PV neurons of KO mice ( Supplementary Fig. 3), suggesting that the main deficit lies in the reduced resting membrane conductance which is mainly contributed by the leak conductance. Reduced latency to PV interneuron firing would narrow the temporal window to integrate excitatory inputs at pyramidal neurons. Consistent with reduced resting membrane conductance (Fig. 2g), I-V curves of PV neurons (Fig. 2h, i) revealed loss of two conduction components: a voltageindependent component and a voltage-dependent delayed outward rectifier. From these altered kinetics, two types of conductances are likely altered in PV-Lmo4KO mice; the voltageindependent two-pore leak (e.g., Kcnk3) and the voltagedependent delayed outward rectifying potassium conductances (e.g., Kv1.2 (Kcna2), Kv3.1 (Kcnc1) or Kv3.3 (Kcnc2)). This was reflected in the elevated resting membrane potential (Fig. 2b) and increased width of the AP, with slower decay slope (Supplementary Fig. 3), respectively. Fast after-hyperpolarization (fAHP) was also reduced in mutant PV interneurons ( Supplementary  Fig. 3), perhaps due to altered function of BK channels. Direct photoactivation of PV interneurons produced less short-term depression of PV-mediated IPSC ( Supplementary Fig. 4b, c) implicating a more prolonged and sustained inhibition of L2/3 pyramidal neurons in PV-Lmo4KO mice upon repetitive stimulation.
In addition to these altered PV interneuron properties, we also observed increased spontaneous inhibitory but not excitatory synaptic inputs onto L2/3 ACC pyramidal neurons and the E/I ratio was drastically reduced (Fig. 3j, k), with a divisive change in the sIPSC-sEPSC relationship (Fig. 3i). Lack of correlation between sIPSC decay time and rise slope indicates these sIPSCs occur at or near the soma of L2/3 pyramidal neurons, where PV interneurons synapse (Fig. 3b, similar to Supplementary Fig. 4d upon direct photoactivation of PV interneurons).
Heightened MD thalamocortical FFI. Given the increased excitability of PV interneurons, we next examined how thalamocortical FFI is affected in PV-Lmo4KO mice. MD thalamocortical projections to the dACC were photo-activated after expressing the photosensitive channelrhodopsin ChR2 and mCherry fluorescent marker in these projections with a viral vector (Fig. 4a). Light stimulation activates a modest monosynaptic EPSC and a strong di-synaptic IPSC (from PV interneurons) onto L2/3 pyramidal neurons of the dACC, as described 18 (Fig. 4b). When the light stimulus was adjusted to produce a monosynaptic EPSC in WT and PV-Lmo4KO, we found that FFI was 4× stronger in PV-Lmo4KO mice than in controls ( Fig. 4b-d), resulting in a reduced E/I ratio (Fig. 4e). The eIPSC-eEPSC amplitude relationship for each individual L2/3 pyramidal neuron showed a divisive change in PV-Lmo4KO mice relative to WT, indicating that a larger feedforward inhibitory current occurs with any given excitatory input in PV-Lmo4KO mice (Fig. 4c). Of note, the divisive change in the I-E relationship here was similar to what we observed above with the sIPSC-sEPSC amplitude (Fig. 3i). The latency of the di-synaptic IPSC (occurring after the EPSC) recorded from L2/3 pyramidal neurons was much shorter (Fig. 4f), consistent with the faster firing of PV neurons from PV-Lmo4KO mice (Fig. 2d). Inhibitory synaptic release probability (Fig. 4g, h) and short-term depression (Fig. 4i, j) were reduced in PV-Lmo4KO mice upon repetitive activation of thalamocortical projections; similar synaptic properties were observed with direct photoactivation of PV interneurons ( Supplementary Fig. 4b, c), supporting that these PV interneurons mediate MD thalamocortical FFI 18 . In contrast, no difference in excitatory synaptic responses of L2/3 pyramidal neurons to MD thalamocortical activation was observed (Supplementary Fig. 5), suggesting that increased FFI is due to increased PV neuron function rather than altered MD thalamocortical excitatory inputs.
Altered FFI will affect the input-output relationship of pyramidal neurons 21 . Consistent with heightened PV-mediated FFI, L2/3 pyramidal neurons were less responsive to MD thalamocortical inputs (Fig. 4k-p). When 10 sample traces recorded 10 s apart were overlaid for WT and KO upon optogenetic stimulation of MD thalamic inputs, we observed that spike timing was less precise in the KO (Fig. 4k). Spike probability of each pyramidal neuron was measured by the spike number per each ten photo-stimuli at various strengths until saturated, and fitted to a sigmoid function 21 . With the same EPSC input (measured at Vh = −60 mV, Fig. 4l, inset), spike probabilities of the PV-Lmo4KO cells were suppressed (Fig. 4l, m), by 30% overall (Fig. 4n). The threshold required to induce AP was increased ( Fig. 4o) and output gain was much reduced (Fig. 4p). It is intriguing that some WT L2/3 pyramidal neurons responded to MD thalamic activation with a double spike at threshold (Fig. 4m), but this never occurred in PV-Lmo4KO mice, perhaps reflecting the shortened latency of PV-mediated FFI (Fig. 2d).
Unexpected reduction in dendritic inhibitory inputs to ACC. L2/3 pyramidal neurons receive at least two types of inhibitory inputs: 1 >peri-somatic inhibitory inputs from L2/3 PV neurons, and 2 >L1 dendritic inhibitory inputs from SST neurons residing in layer 2/3 30 (Fig. 5a). Dendritic inhibition versus proximal (peri-somatic) inhibition affects pyramidal neuron activation differently: dendritic inhibition vetoes dendritic electrogenesis and raises the threshold of pyramidal neuron activation whereas proximal peri-somatic inhibition is important to set the gain (the slope of the output responses relative to input stimuli) and affects how pyramidal neurons integrate incoming signals as well as the timing and synchrony of output (spiking) 22,31,32 . ACC L2/3 pyramidal neurons output onto L5 pyramidal neurons that project to MD thalamus as well as to the upper layers of the cortex; the latter cross-laminar projections are thought to be important for inter-laminar feedback integration. As part of this feedback loop, L5 pyramidal neurons activate SST inhibitory interneurons in L2/3 that synapse onto dendrites of L2/ 3 pyramidal neurons [33][34][35][36] .
To examine this feedback circuit, electrical stimulation was applied to layer 5. Remarkably, in PV-Lmo4KO mice, this produced a smaller FFI onto L2/3 pyramidal neurons with no change in the monosynaptic EPSC ( Fig. 5b-d), resulting in a markedly increased E/I ratio (Fig. 5e). There was a subtractive Rise up slope (ms/pA) Rise up slope (ms/pA)  36 . We observed a positive correlation between the signal rise (ms/pA) and decay (tau, ms), proportional to the distance traveled, providing evidence that these are (SST-mediated) dendritic inputs, rather than PV-mediated peri-somatic inputs ( Supplementary Fig. 6b). Upon layer 5 stimulation, excitatory inputs to L2/3 showed an increase in paired-pulse ratio and reduced short-term depression for PV-Lmo4KO mice ( Supplementary Fig. 7), which was different from what we observed for excitatory inputs from the MDT ( Supplementary Fig. 5) and may reflect changes in the plasticity of selected local excitatory circuits.
Beyond the feedback loop from L5, L2/3 pyramidal neurons also integrate cortico-cortical and other higher-order thalamic excitatory inputs onto their dendrites in L1. At the same time, these distal dendrites receive inhibitory inputs mediated largely by SST interneurons (and some other local interneurons), but not from PV interneurons 39,40 , to provide context-dependent modulation of cortico-cortical excitatory inputs 20 . To determine whether FFI from the cortico-cortical inputs shows a similar homeostatic change as observed in the feedback loop from L5 stimulation (Fig. 5), we carried out electrical stimulation in L1.
Interestingly, in PV-Lmo4KO mice, L1 stimulation elicited a reduced FFI response similar to L5 stimulation (compare Supplementary Fig. 8e, f to Fig. 5c, e), suggesting that these (SST) interneurons are also affected by a homeostatic response to Lmo4 ablation in PV neurons. The reduced dendritic inhibition observed here was surprising since Lmo4 is ablated in PV but not SST interneurons of PV-Lmo4KO mice and suggests a homeostatic compensation in the inhibitory circuits.
Reduced dendritic inhibition, response threshold, and range. As a consequence of reduced dendritic FFI, the threshold to activate L2/3 pyramidal neurons in response to layer 1 stimulation (Fig. 6a) was also markedly reduced ( Fig. 6b-e). The gain (i.e., the output response relative to input stimuli) of individual neurons was not changed (Fig. 6f), but the population dynamic range (i.e., the range of stimuli that the population of neurons responds to and can distinguish) was much reduced in L2/3 pyramidal neurons of PV-Lmo4KO mice (Fig. 6g). At 3.25 V, a stimulus that recruited 92 % of KO cells (22/24) only recruited 30% of WT cells (7/24). The dynamic range to recruit 92% of WT neurons (black dashed line) was twice that of KO neurons.
As reported by Pouille et al. 21 , the thresholds and output gains showed no correlation with resting cell membrane potential or conductance of L2/3 pyramidal neurons ( Supplementary Fig. 9), indicating that integration of dendritic inputs rather than the pyramidal neurons' intrinsic properties controls their response threshold and gain.
The leftward shift in the input threshold to recruit pyramidal neurons (without altering their gain) in PV-Lmo4KO mice could be reproduced in wild-type mice using low-dose bicuculline (2 μM) to partially inhibit GABAa receptors and was reversed upon washout ( Supplementary Fig. 10b-d). Interestingly, if GABA receptors were completely blocked with 10 μM bicuculline, both the threshold and the gain were changed ( Supplementary Fig. 11e, f), an effect mediated by APV-sensitive NMDA receptors ( Supplementary Fig. 11a, b, e); this is consistent with a previous report that NMDA receptors at distal dendrites increase the output gain 22 . In contrast, upon photo-activation of MD thalamocortical inputs, we did not detect an NMDA component either by voltage clamp (Supplementary Fig. 12a) or current clamp ( Supplementary Fig. 12b), indicating that MD thalamocortical inputs target the apical trunk of proximal dendrites 33 that lack NMDAR 41 , rather than distal dendrites.

Ablation of PTP1B in PV neurons prevents ASD-like deficits.
Since LMO4 is a dual function protein affecting transcription factors 42 as well as an endogenous inhibitor of tyrosine phosphatase PTP1B 25,26 , it was important to determine to what extent PTP1B hyperactivity might contribute to the functional deficits and phenotype that we observed in PV-Lmo4KO mice. To this end, we generated PV-DKO (double knock-out) mice that ablate PTP1B together with Lmo4 in the PV neurons. The social interaction deficits (Fig. 7a) and repetitive behaviors (Fig. 7b) observed in PV-Lmo4KO mice were prevented in PV-DKO mice (also see Supplementary Figs. 13, 14). Consistent with rescued behaviors, the number of c-Fos+ neurons after social interaction were not different between WT and PV-DKO mice (Fig. 7c).
PV interneurons of PV-DKO mice showed normalized resting membrane potential (Fig. 7d), excitability (Fig. 7e) and latency to first action potential (Fig. 7f). Membrane conductance (Fig. 7g) was restored, with a normal leak conductance, yet with an exaggerated delayed rectifying current (Fig. 7h, i). The action potential width and fast after hyperpolarization (fAHP) were also normalized ( Supplementary Fig. 15). Of note, no difference in basic membrane properties of L2/3 pyramidal neurons was detected (resting membrane potential: In line with the normalization of most electrophysiological properties of PV interneurons in DKO mice, optogenetic Fig. 4 Increased feedforward inhibition in the MD thalamus-dACC circuit in PV-Lmo4KO mice. a Synaptic inputs recorded at dACC layer 2/3 pyramidal neurons (PyNs) after optogenetic activation of thalamocortical projections (left, middle). Right, diagram of MD thalamus-dACC FFI circuit. b Sample traces of evoked inhibitory (eIPSC) and excitatory (eEPSC) currents recorded at layer 2/3 pyramidal neurons of wild-type (WT, black) and PV-Lmo4KO (KO, red) mice. Light stimulus was adjusted to trigger monosynaptic EPSC in WT and PV-Lmo4KO. Traces were enlarged to reveal short latency of monosynaptic EPSC versus longer latency of disynaptic IPSC. c E-I relationship fitted with linear regression from paired eIPSC and eEPSC amplitudes from each individual cell. d Increased eIPSC amplitude but no difference in eEPSC amplitude, resulted in a reduced (E/I) ratio (e). f Reduced latency between the rising phase of inhibition and excitatory synaptic responses (I-E) in KO. g Sample traces of the eIPSC at a given paired-pulse interval. h Paired-pulse ratio. i Sample traces of eIPSC in response to repetitive stimuli at 20 Hz reveal short-term depression j. Inset, the mean of 5th/1st eIPSC ratio, which represents the degree of short-term depression at various stimulation frequencies. n = 21 cells/7 WT, 23 cells/8 KO mice for (c-j). k-p Under current clamp, L2/3 pyramidal neurons were less responsive to thalamocortical inputs in PV-Lmo4KO mice. k Sample traces of 10 sweeps, 10 s apart were overlaid for WT and KO. l The spike probability from a representative WT or KO L2/3 neuron, plotted against stimulation (light intensity*duration). Insets, representative traces (l) and histogram (m) show EPSC were of similar amplitudes at holding potential −6 mV with 2 ms light stimulation. m Spike probabilities of all WT (15 cells/6 mice) and KO (8 cells/4 mice) L2/3 neurons recorded were sigmoid fitted to obtain (n) maximum spike probability, o stimulation threshold to elicit 50% peak spike and (p) output gain. n = 15 cells/6 WT, 8 cells/4 KO mice for (l-p). *, **, ***, p < 0.05, 0.01, 0.001, respectively. NATURE COMMUNICATIONS | https://doi.org/10.1038/s41467-020-14813-z ARTICLE NATURE COMMUNICATIONS | (2020) 11:1017 | https://doi.org/10.1038/s41467-020-14813-z | www.nature.com/naturecommunications activation of the MD thalamic projections also showed a normalized PV-mediated FFI (Fig. 8a), with restored E/I ratio (Fig. 8b) and latency of inhibition (Fig. 8c). The short-term depression of PV-mediated FFI was also normalized ( Supplementary Fig. 16). Furthermore, the stimulation threshold (Fig. 8d, e, g), the maximum spike probability (Fig. 8f), and the output gain (Fig. 8h) of L2/3 neurons to MD thalamocortical stimulation were all restored to WT levels in DKO mice.   Most importantly, SST-mediated FFI of L2/3 pyramidal neurons upon L1 stimulation was also normalized in DKO mice (Fig. 8i-m). Together, these results show that the ASD-like behaviors and electrophysiological changes observed in PV-Lmo4KO mice are almost entirely due to disinhibition of PTP1B through loss of Lmo4.

Discussion
Our study shows that selective ablation of Lmo4 in PV interneurons causes ASD-like behaviors. At the ACC, L2/3 pyramidal neurons receive faster and stronger MD thalamocortical FFI due to increased excitability of PV interneurons and their shortened latency to fire. At the same time, dendritic inhibition of these pyramidal neurons is reduced, lowering their activation threshold without altering their output gain. This compensatory reduction in dendritic inhibition occurs at the expense of the range of stimuli to which L2/3 pyramidal neurons can respond and distinguish. These changes were all dependent on the selective activation of the tyrosine phosphatase PTP1B, since they were resolved by simultaneous ablation of PTP1B and Lmo4 in the same PV interneurons in PV-DKO mice.
The most parsimonious explanation is that ablation of the endogenous PTP1B inhibitor Lmo4 25,26 activates PTP1B and leads to a myriad of post-translational tyrosine dephosphorylation events. Our previous study reported that ablation of Lmo4 in glutamatergic neurons (in Camk2αCre/Lmo4flox mice) hyperactivates PTP1B that dephosphorylates mGluR5 and impairs mGluR5 function 26 . mGluR5 function may be similarly affected in PV-Lmo4KO mice. Altered mGluR5 function has been associated with autism-like behaviors in mice lacking Gprasp2 43 . Moreover, PV neuron-specific ablation of mGluR5 induces social novelty deficits and compulsive repetitive behaviors associated with loss of PV neurons and reduced inhibitory currents (mIPSC) at the hippocampus 44 . In contrast, our PV-Lmo4KO mice display social deficits and repetitive behaviors without detectable loss of PV neurons. Moreover, both spontaneous and evoked PV-mediated inhibitory currents are increased in L2/3 ACC pyramidal neurons of PV-Lmo4KO mice. Although the present study focused on PV neuron function in the ACC, PV neurons in other brain regions including cerebellum 4 and hippocampus 45 , could also contribute to ASD-like behavior deficits in PV-Lmo4KO mice. Whether compensatory changes in SST-mediated inhibitory inputs are specific to the ACC or also occur in other brain regions is an important question for future studies. Our previous work indicated that PTP1B activity is normally suppressed by LMO4 in glutamatergic neurons, since knockdown or pharmacological inhibition of PTP1B had no observable effect on behaviors of wild-type mice 26 . Similarly, here we saw no overt effect of PTP1B ablation in PV neurons on social interaction or repetitive behaviors compared to wild-type mice, suggesting that PTP1B function is also normally quiescent in PV neurons. Here, we observed that loss of Lmo4 in PV interneurons led to a depolarized resting membrane potential and a broadened duration of the action potential likely through reduced leak conductance (Kcnk3) and delayed rectifying potassium (Kv1.2 or Kv3) conductance, respectively. Both functions are regulated by tyrosine phosphorylation [46][47][48] . These changes in PV-Lmo4KO mice were rescued by PTP1B ablation, suggesting that tyrosine phosphorylation-dependent functional deficits in potassium conductances may account for the altered PV function in these mice.
Since PV neuron-specific genetic ablation of PTP1B was sufficient to restore electrophysiological properties of Lmo4-deficient PV neurons and to prevent the behavioral and inhibitory circuit deficits in PV-Lmo4KO mice, this raises important questions: Might PTP1B hyperactivity in PV interneurons contribute to ASD in humans? Given that the inhibitory circuit is plastic and adaptive changes in adults can occur within days 49 , would pharmacological PTP1B inhibition normalize PV neuron properties and restore SST-mediated inhibition?
The di-synaptic nature of fast-spiking PV interneuronmediated FFI ensures that inhibitory inputs arrive with a short delay and set the window of opportunity during which L2/3 pyramidal neurons can integrate excitatory dendritic inputs and generate action potentials. However, in PV-Lmo4KO mice, this window of opportunity is markedly narrowed due to shortened latency of PV interneuron firing coupled with a fourfold heightened PV-mediated inhibition (Fig. 9a). These changes predict NATURE COMMUNICATIONS | https://doi.org/10.1038/s41467-020-14813-z ARTICLE that MD thalamocortical inputs would be less effective in activating L2/3 pyramidal neurons (Fig. 9b, c). Indeed, we observed a 30% reduction in spike probability (Fig. 4n).
On the other hand, we found an overall reduction in inhibitory inputs from L1 and L5 to distal dendrites of L2/3 pyramidal cells in PV-Lmo4KO mice, a feature that should greatly increase pyramidal neuron firing in response to cortico-cortical or other thalamocortical excitatory inputs (Fig. 9d, e). Reduced SSTmediated inhibition would lower the threshold to activate L2/3 pyramidal neurons and background noise would fail to be filtered in PV-Lmo4KO mice. Moreover, this would limit the dynamic range of the L2/3 pyramidal neuron population; a stimulus intensity that recruits fewer than one third of L2/3 pyramidal neurons in wild-type mice would have recruited nearly all of the L2/3 pyramidal neurons in PV-Lmo4KO mice (Fig. 6g). The deficit in dendritic inhibition of L2/3 pyramidal neurons was surprising, since Lmo4 is ablated in PV (but not SST) interneurons of PV-Lmo4KO mice. Also surprising was the fact that both cortico-cortical FFI in layer 1 and intralaminar FFI from layer 5 to layer 1 were similarly reduced. This could reflect a compensatory change of the local cortical circuit 49 in response to PV interneuron hyperactivity in order to maintain homeostatic inhibition of pyramidal neurons 3 . Alternatively, this could be a hard-wired direct consequence of increased PV-mediated inhibition of SST neurons 50 . This could set up a vicious cycle, since SST neurons tonically inhibit PV neurons 20 and suppressing SST neurons enhances PV-mediated FFI to regulate thalamocortical inputs at the ACC 18 . Whether increased PV-mediated inhibition also causes homeostatic changes in other inhibitory interneurons (e.g., VIP) remains to be determined.
Increased thalamocortical FFI was also reported in a mouse model of ASD induced by fetal exposure to valproic acid 51 , where reduced function of Kcnc1 (Kv3.1) in PV interneurons at the mPFC has been suggested 52 . The notion that altered E/I balance underlies ASD-related deficits may be over-simplified, as both reduced 2 and increased 3 E/I balance have been measured in mouse models of ASD. A recent elegant study of L2/3 pyramidal neurons of the somatosensory cortex in three ASD mouse models by Feldman's group reports that increased E/I ratios are a homeostatic mechanism stabilizing synaptic drive, rather than predicting network hyperexcitability in autism 3 .
Our finding of a homeostatic regulation of the inhibitory circuits at the ACC of PV-Lmo4KO mice, compartmentalized to peri-somatic versus dendritic inputs, adds a further layer of complexity: a reduced E/I ratio from the MD thalamocortical circuit was accompanied by an increased E/I ratio from the cortico-cortical and other higher-order thalamic inputs. It will be interesting to examine whether compartmentalized changes in the E/I ratio also occur in inhibitory circuits of other brain regions of PV-Lmo4KO mice, e.g., in the somatosensory cortex. Another pressing question is whether compartmentalized reciprocal changes in inhibitory circuits occur in other ASD models.
The strength of dendritic and peri-somatic inhibition of cortical pyramidal neurons is dynamically regulated. The compartmentalized switch between dendritic and peri-somatic inhibition offers context-dependent modulation of excitatory activity of pyramidal neurons, increases their flexibility and experiencedependent learning in processing different information inflows 20,22 that are important for social interaction. Activity state-dependent and cortical layer-specific changes in dendritic versus peri-somatic inhibition have also been reported 30,53 . During exploratory behaviors (whisking), L2/3 pyramidal neurons of the somatosensory cortex receive reduced dendritic inhibition from L2/3 SST interneurons but increased peri-somatic inhibition from L2/3 PV interneurons 53 . On the other hand, L4 SST inhibition that targets L4 PV interneurons is increased during whisking 30 , and thereby decreases thalamocortical FFI. Unlike the somatosensory cortex, the ACC lacks layer 4 pyramidal neurons and layer 2/3 pyramidal neurons receive direct thalamic inputs. Whether similar changes in inhibition also function at the ACC during exploration and social interaction remain to be determined.
Our c-Fos data revealed that social interaction increased activation of L2/3 pyramidal neurons in wild-type mice but this response was reduced in PV-Lmo4KO mice with social interaction deficits. Reduced dendritic (SST-mediated) inhibition onto ACC L2/3 pyramidal neurons observed in PV-Lmo4KO mice hints that mutant mice might already be overwhelmed with excessive cortico-cortical and other higher-ordered thalamic . c Thus, fewer MDT excitatory inputs would succeed in eliciting action potentials at L2/3 pyramidal neurons (PyN); not only background noise (thin black arrows), some valid signals (thick black arrows) would be also filtered away. d Adaptive reduction of SST-mediated dendritic inhibition lowers the activation threshold of L2/3 pyramidal neurons, as reflected by a subtractive change between the E-I. This would compromise the ability of the L2/3 pyramidal neuron population to respond/distinguish the range of stimulus intensities from cortical-cortical inputs (CC) or projections from other higher-order thalamic nuclei (TC) (Fig. 6e). e In contrast to wild-type mice (PVWT), background noise (thin black arrows) would fail to be filtered in PV-LMO4KO (PVKO) mice.
inflows at rest and may be incapable of further increasing attention to (or to distinguish) relevant information during exploration or social interactions. Impaired capacity of neurons to regulate state-dependent synaptic homeostasis may be a major mechanism underlying ASD.

Methods
Mice. PV-Cre/Lmo4 flox/flox (PV-Lmo4KO) mice were obtained by mating PV-Cre mice (gift from Dr. Sylvia Arber 1 , Jackson Lab # 017320) to Lmo4 flox/flox mice (gift from Dr. Jane Visvader 2 ). PV-DKO mice were generated by mating PV-Cre/ Lmo4 flox/flox mice to PTP1B flox/flox mice (gift from Dr. Benjamin Neel 3 , Jackson Lab #012679). PV-Cre/Lmo4 flox/flox mice were bred with wild-type mice to obtain PV-Cre/Lmo4 flox/WT mice for optogenetic studies using viral vectors. All mice were bred on the C57BL6 background for over 12 generations. Mice were fed with regular chow, and randomly assigned to experimental groups. All procedures for animal use were approved by the University of Ottawa Animal Care and Veterinary Service, and were performed according to institutional guidelines and in accordance with those of the Canadian Council on Animal Care. All studies were conducted using equal numbers of male and female mice.
Behavioral tests. Mice~2-4 months of age were subjected to the following behavior tests modified from ref. 4 with at least 2 days rest in between each test. Investigators were blinded to genotype. Mice were placed in an open field box (45 cm × 45 cm × 45 cm; MED Associates, St. Albans,VT). Red lighting was used for this test to minimize anxiety form being in a novel environment. The social target mouse (novel or littermate) is contained within a 5.5 × 9.6 × 30 cm 3 (W × L × H) wire mesh cage placed on the opposite side of the square box as described in Fig. 1a. Mice were first habituated for 5 min in the arena with empty mesh cages. Following a 5-min rest in their home cages, they were then videotaped and monitored for 5 min with or without a social target mouse from a mounted camera connected to a computer using Ethovision video tracking software (Noldus). The amount of time they spent at the interactive zone (8 cm away from each side of the mesh cage) was recorded. The discriminative index for sociability was calculated by the formula Each mouse was placed in a new cage and videotaped and monitored for 10 min. The time spent licking itself was recorded as self-grooming time.
Mice were placed individually in a new cage with fresh bedding and videotaped. The time spent digging the bedding over 5 min was recorded.
Mice aged 8-12 weeks were habituated in the testing room 30 min prior to testing in their normal home cage. Olfaction was tested as described 4 , by measuring the duration of the sniffing response to non-social, then social odors. Odors were presented on cotton-tipped applicators in three consecutive trials per odorant stimulus (2 min per trial) in the following order: water, almond extract, banana extract, social odor 1, social odor 2. Social odors were swipes from cages containing age-matched unfamiliar same sex and opposite sex mice.
After >40 min recovery time at room temperature, slices were transferred to a chamber perfused with ACSF 2 ml/min at room temperature. For experiments to test MD thalamic-ACC inputs, slices containing MD thalamus and dACC from each mouse were cut at a 300-μm thickness and imaged to examine the location and extent of ChR2 expression in the MD. Mice were excluded if the extent of infection was too large and had leaked into surrounding brain regions. Rodent MD thalamus lacks interneurons; therefore, all ChR2-infected neurons are expected to be relay projection neurons 7 . For PV-IN, recorded neurons were visually identified using immunofluorescence and infrared differential interference contrast with video microscopy.
In ACC, layers were defined by distance from the pial surface: L2/3: 200-350 μm; L5: 400-550 μm; we recorded our L2/3 pyramidal neurons between 250 and 350 μm. For L5 stimulation, electrode was placed in the middle of L5 which is 450-500 μm away from pia 8 . For L1 stimulation, the electrode was placed on the outer layer of pia. For both L1 or L5 stimulation, the pulse was limited to 0.1 ms with various strengths.
During voltage clamp to obtain spontaneous EPSC/IPSC or evoke EPSC and IPSC to calculate the E/I ratio, recording electrodes were filled with a solution of (in mM): 130 CsMeSO 3 , 10 KCl, 10 HEPES, 1 EGTA, 2 MgCl 2 , 10 Na +phosphocreatine, 4 Na + -ATP, 0.3 mM GTP (pH 7.3 adjusted with CsOH, 295-300 mOsM); junction potential: 14 mV. In our preparation, the reversal potential for eIPSC was~−45 mV measured under NBQX plus D-APV, so we clamped cells at −45 mV for sEPSC/eEPSC. The reversal potential of EPSC was measured at +20 mV in the presence of bicuculline or picrotoxin, so the cells were held at +20 mV to record sIPSC/eIPSC. During current clamp for spike probability, the 130 mM CsMeSO 3 was replaced with 137 mM K + -gluconate and 10 mM KCl was reduced to 3 mM in the recording solution (the Cl − reversal potential was~−74.4 mV); junction potential: 15 mV. Unless specified, all potentials mentioned were not adjusted for junction potentials. To ensure slice viability and input continuity between measurements of individual neurons, voltage clamp experiments were carried out in the same cells after current clamp protocols. This approach was also used to obtain basic properties of PV interneurons and L2/3 pyramidal neurons.
Series resistance was compensated at~50% and monitored during the process, cells with series resistances larger than 25 MΩ were excluded. To characterize PV interneuron and pyramidal neuron membrane properties, we measured resting membrane potential immediately after the rupture of the neuronal membrane. Resting membrane conductance was measured by linear fit I-V curve under voltage clamp between −70 and −40 mV with holding at −50 mV.
To compare fast-inactivating A-type (I A ) potassium conductance, TTX (1 μM), Cd 2+ (100 μM), TEA 20 mM, 4-AP (80 μM) and Ba 2+ (100 μM) were included in the recording bath solution and currents were evoked in response to 10 mV membrane potential steps from −80 to +40 mV with holding membrane potentials at −80 and −40 mV respectively. The difference between the two currents at these two holding potentials revealed the fast-inactivating A-type (I A ) potassium conductance. The delayed rectifier conductance (K DR ) component from each single cell was isolated by subtracting the voltage-independent component (extrapolated from linear fit between −70 and −40 mV) from the outwardly rectifying currents at each step potential after −40 mV.
Cell membrane excitability was measured under current clamp beginning from −60 mV and the number of action potentials evoked by a series of current injections was determined for all cells. Action potentials elicited by a 140 pA 500 ms square pulse were pooled obtain parameters such as half width, decay slope and fAHP for distribution comparisons and the first ten spikes for real value comparisons ( Supplementary Figs. 2, 15). For spontaneous events, 150-200 events recorded from each individual cell were analyzed for amplitude, instant frequency and decay tau analysis; the average values of the 150-200 events from each individual cell were pooled together and compared between groups. The duration of light stimulation for feed-forward inhibition experiments was limited to no more than 3 ms at various intensities.
Spike probability of each pyramidal neuron was measured by the spike number per each ten photo-stimuli or electrical inputs at various strengths until saturated. Data for each single cell or group of cells (membrane excitability and spike probability measurements) were plotted to a sigmoid function to obtain slope or gain and threshold for 50% maximum APs probability, as described 9 . In our experiments, the threshold refers to the electrical stimulation or photoactivation input strength to reach 50% of the maximum response.
To evoke synaptic transmission by activating ChR2, blue light pulses were delivered from a single wavelength LED system (λ = 465 nm; Plexon) triggered by a TTL (transistor-transistor logic) signal from the Clampex software (Molecular Devices). A Multiclamp 700B amplifier connected to a Digidata 1440 interface (Molecular Devices) was used for whole cell recordings and analyzed using pClampfit 10.1. Recordings were filtered at 3 kHz, digitized at 10 kHz. To compare spike probabilities in response to MD thalamocortical inputs and to test how FFI affects input-output relationships in WT, KO and DKO mice, every cell was held at −60 mV to measure EPSC (close to the Clreversal potential,− 74.4 mV, after adjusting for the 15 mV junction potential) with a 2 ms maximum blue light strength (325 mA) before obtaining spike probability. To ensure cells were given similar input strengths, only cells responding to a photo-stimulus with a monosynaptic EPSC between −300 and −400 pA were used for analysis.
Where spike probability was measured following L1 or optogenetic MD fiber stimulation under current clamp, cells were subsequently subjected to voltage clamp and held at −40 mV to measure an EPSC -IPSC sequence and ensure FFI was triggered by the stimulation or input range (Supplementary Fig. 8). The relative amplitudes of "EPSC" and "IPSC" were analyzed.
Statistical analysis. All results are presented as means ± SEM. For between two group comparisons, a t-test or one-way ANOVA was used to compare normally distributed data. Data that were not normally distributed as determined by the Shapiro-Wilk normality test were compared with Mann-Whitney Rank Sum test; Z-tests were used to compare parameters obtained from sigmoid and exponential decay functions. For three groups or more, one-way ANOVA or ANOVA on ranks (if normality not passed) was applied, and Bonferroni correction was applied for multiple pairwise testing using two-tailed Student's t test for post hoc analysis; for paired-pulse ratio (PPR), short-term depression (STD) and action potential parameters comparisons, two-factor ANOVA compared the main effects of genotype and two-factor interactions; repeat ANOVA was also used for STD analysis. All tests were carried out with SPSS statistical software and differences in means were considered significant at p < 0.05. See Supplementary Tables 1, 2 for summaries of statistical analyses.
Reporting summary. Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability
The source data underlying all data