Published: January 21, 2020

Introduction {#sec1}
============

In the neocortex, sensory integration is accomplished through the coordinated activity of neuronal networks across different cortical layers and columns (reviewed in [@bib2], [@bib22], [@bib23]). These hardwired anatomical connectivity patterns between several neuron types define specific pathways and enable salient flow of information across and within different cortical layers.

Functional cortical networks result from direct contact between neurons and indirect feedforward and feedback connections from intercalated neurons, whose recruitment strength and excitability contribute to the formation and dissolution of neuronal ensembles ([@bib13]). These interposed neurons are mostly (but not only) inhibitory. Importantly, fast synaptic inhibition plays a fundamental role in shaping and sculpting virtually all forms of cortical activity ([@bib5], [@bib15], [@bib17], [@bib39], [@bib64], [@bib68], [@bib69], [@bib71], [@bib77]).

Neocortical inhibition is provided by a rich diversity of GABAergic interneurons ([@bib3], [@bib36], [@bib51]). In particular, parvalbumin (PV)-positive, perisomatic targeting interneurons account for nearly half of all cortical inhibitory cells ([@bib8], [@bib26]). Perisomatic inhibition from PV cells exerts a plethora of functions within neocortical circuits, including control of pyramidal neuron (PN) output activity and gain during sensory perception ([@bib5], [@bib16], [@bib39], [@bib67], [@bib77]) and the ability to drive fast cortical oscillations underlying several aspects of cognition ([@bib7], [@bib14], [@bib17], [@bib33], [@bib68], [@bib74]), thus likely acting as a synchronizing mechanism and likely promoting the generation of cell assemblies ([@bib13]).

PV cell-mediated perisomatic inhibition equalizes excitation proportionally, resulting in a tight balanced excitation-to-inhibition (E/I) ratio in single PNs ([@bib4], [@bib78]). Notably, perturbation of PN activity *in vivo* has been shown to modify PV-mediated perisomatic GABAergic inhibition ([@bib78]). We have found that postsynaptic depolarization or bursts of action potentials in layer (L) 5 PNs of the mouse barrel cortex induces long-term potentiation of inhibition (LTPi), which is selective for PV cell transmission and sharpens the time window of synaptic integration ([@bib42]).

In the last four decades, long-term plasticity of synaptic transmission has been studied extensively at glutamatergic synapses because it is considered the cellular correlate of learning and memory ([@bib45]). Importantly, inhibitory synapses are also highly plastic ([@bib18], [@bib19], [@bib20], [@bib28], [@bib29], [@bib38], [@bib51], [@bib58], [@bib72]), but the functional role of GABAergic plasticity is unknown (but see [@bib52], [@bib73]).

Here we set out to investigate how the plasticity of PV cell-mediated perisomatic GABAergic synapses modulates several computations performed by single L5 PNs of the mouse barrel cortex (S1). Using *in utero* electroporation, we expressed light-sensitive opsins in L2/3 PNs of the mouse S1. We demonstrate that activation of L2/3 induces robust feedforward inhibition (FFI) on L5 PNs, mostly mediated by PV basket cells ([@bib37], [@bib50]). FFI could be strongly potentiated by cell-autonomous postsynaptic paradigms. LTPi-FFI modified the input/output relationship of L5 PNs and strongly modulated information flow across cortical layers. Moreover, LTPi-inducing bursts affected the temporal association of PN spiking with γ-oscillations both *in vitro* and *in vivo*. These results were captured by a computational model, indicating that plasticity of PV cell-dependent perisomatic inhibition can strongly modulate single PNs at the single-cell and network levels.

Results {#sec2}
=======

Burst Firing of L5 PNs Selectively Potentiates Feedforward GABAergic Inputs {#sec2.1}
---------------------------------------------------------------------------

L2/3 PNs play a prominent role in activating deep cortical output layers. We have demonstrated previously that bursts of action potentials (APs) in L5 PNs induce LTPi selectively at synapses from PV cells ([@bib42]). We therefore asked whether L2/3 activation triggers PV cell-mediated FFI onto L5 PNs and whether this is plastic.

Using *in utero* electroporation, we expressed the light-sensitive opsin channelrhodospin 2 (ChR2) in a large fraction of L2/3 PNs ([Figure 1](#fig1){ref-type="fig"}A). We then performed whole-cell patch-clamp recordings from large L5 PNs in acute slices of the barrel cortex (S1, barrel field) from mice that had been electroporated previously *in utero* ([Figures 1](#fig1){ref-type="fig"}A and 1B). Brief (1-ms) stimulation with blue light (λ = 470 nm) of L2/3 ChR2+ PNs induced a composite postsynaptic potential (PSP) in L5 PNs recorded in current clamp mode. This composite PSP was made of an early excitatory PSP (EPSP) triggered by L2/3 PNs ([Figure S1](#mmc1){ref-type="supplementary-material"}), which was followed by a GABAergic inhibitory PSP (IPSP; [Figures 1](#fig1){ref-type="fig"}B and 1C, top panel; [Figure S1](#mmc1){ref-type="supplementary-material"}). This inhibitory component had typical disynaptic latencies ([Figure S1](#mmc1){ref-type="supplementary-material"}) and was likely triggered by perisomatically targeting PV cells. Indeed, activation of L2/3 PNs recruited L5 PV interneurons efficiently ([Figure S1](#mmc1){ref-type="supplementary-material"}; [@bib37], [@bib50], [@bib21]), and late GABAergic responses recorded in voltage-clamp had rise time values compatible with light-evoked synaptic events mediated by PV cells ([Figures S1](#mmc1){ref-type="supplementary-material"}G and S1H; [@bib42]).Figure 1Burst Firing of L5 PNs Selectively Potentiates Feedforward GABAergic Input(A) *In utero* electroporation of ChR2 and red fluorescent protein (RFP) in L2/3 PNs of the mouse S1.(B) Scheme of the recording configuration.(C) Average (10 sweeps) current-clamp traces of the EPSP-IPSP composite response recorded in L5 PNs upon photostimulation of L2/3 PNs before (black, top) and after (bottom, red) inducing LTPi.(D) LTPi of light-IPSPs (top graph) of the cell shown in (C). The bottom graphs indicate light-EPSP slope, input resistance (Rin), and resting membrane potential (Vm) of the same cell.(E) Population graph of LTPi in L5 PNs.(F) Plot illustrating the relative change of light-IPSPs in response to burst firing of individual PNs (after 20 min). Dark circles, LTPi-expressing PNs; light gray circles, PNs not expressing LTPi; Bsl, baseline.(G) Graphs showing average depolarizing slopes, areas, and the EPSP/IPSP ratio of composite PSPs in Bsl and after postsynaptic bursts. In some cases, the error bars are too small to be visible. n.s., not significant. ^∗^p \< 0.05, ^∗∗^p \< 0.01, with paired t test.Population data are illustrated as mean ± SEM.

Importantly, in response to postsynaptic AP bursts of L5 PNs (5 APs at 100 Hz, repeated 15 times every 10 s), we observed an increase in amplitude of the GABAergic component of the light-evoked PSP (light-IPSP) that persisted for more than 30 min ([Figures 1](#fig1){ref-type="fig"}C--1F), which we termed long-term potentiation of feedforward inhibition (LTPi-FFI). LTPi-FFI was observed in ∼65% of the recorded PNs ([Figure S1](#mmc1){ref-type="supplementary-material"}I, includes all cells tested for LTPi-FFI in [Figures 1](#fig1){ref-type="fig"}, [2](#fig2){ref-type="fig"}, and [4](#fig4){ref-type="fig"}), and it induced an increase in light-IPSPs amplitude of at least 50% of the baseline amplitude (0.73 ± 0.2 versus 2.38 ± 0.58 mV, light-IPSP baseline versus 20 min after AP bursts, respectively; n = 10, p = 0.0098, Wilcoxon matched-pairs signed-rank test; [Figure 1](#fig1){ref-type="fig"}F; [Figure S3](#mmc1){ref-type="supplementary-material"}A). We reported previously that this form of cell-autonomous plasticity is sensitive to retrograde nitric oxide (NO). Moreover, NO-mediated plasticity was expressed at PV-PN but not somatostatin (SST) interneuron-PN synapses ([@bib42]). Likewise, LTPi-FFI was also sensitive to pharmacological inhibition of the canonical NO receptor guanylylcyclase (GC) with 1H-{1,2,4}oxadiazolo{4,3-a}quinoxalin-&-dione (ODQ, 10 μM; 0.73 ± 0.16 versus 0.89 ± 0.31 mV, light-IPSP baseline versus 20 min after AP bursts, respectively, in the continuous presence of ODQ; n = 8, p = 0.74, Wilcoxon matched-pairs signed-rank test; [Figures S1](#mmc1){ref-type="supplementary-material"}J and S1K). Importantly, LTPi-FFI did not affect the slope of the depolarizing EPSP component (baseline, 0.904 ± 0.19 mV/ms; LTPi, 1.024 ± 0.0.19 mV, n = 10, p = 0.16, Wilcoxon matched-pairs signed-rank test; [Figures 1](#fig1){ref-type="fig"}D and 1G1). However, because of the potentiation of the hyperpolarizing (IPSP) component, the area of the light-evoked EPSP decreased significantly (baseline, 42.46 ± 7.47 mV/ms; LTPi, 15.73 ± 4.8 mV/ms; n = 10, p = 0.0036, paired t test; [Figure 1](#fig1){ref-type="fig"}G2). Consequently, LTPi-FFI strongly reduced the E/I ratio, measured as the EPSP area divided by the total composite PSP area (0.64 ± 0.07 versus 0.28 ± 0.11, baseline versus LTPi; n = 10, p = 0.0023, paired t test; [Figure 1](#fig1){ref-type="fig"}G3). In the presence of the NO receptor inhibitor ODQ, postsynaptic AP bursts failed to induce changes in the E/I ratio (0.79 ± 0.28 versus 0.636 ± 0.1 mV, baseline versus 20 min after AP bursts, in the continuous presence of ODQ; n = 8, p = 0.12, Wilcoxon matched-pairs signed-rank test; [Figure S1](#mmc1){ref-type="supplementary-material"}L).Figure 2LTPi-FFI Induces a Multiplicative Gain Modulation of PN Output Spiking(A) Schematic of the recording configuration.(B) Averaged (10 sweeps) current-clamp traces of EPSP-IPSP responses in L5 PNs upon photostimulation of ChETA+ L2/3 PNs before (black) and after (red) LTPi induction.(C) Vm (left panel) and Rin (right panel) throughout the experiment.(D) Spiking output of L5 PNs upon photostimulation of L2/3 PNs at 15-Hz frequency during Bsl and after LTPi-FFI (red trace).(E) Averaged output firing rate of L5 PNs upon photostimulation of L2/3 PNs at different frequencies before (Bsl, light gray) and 15 min after LTPi induction (LTPi-FFI, dark gray); ^∗∗∗∗^p \< 0.0001, Friedman test. Dashed lines are fits to a Hill function.(F) Offset (light gray) and gain (dark gray) changes in the I/O function after LTPi-FFI from the fit in (D) (^∗∗^p \< 0.01, Wilcoxon matched-pairs signed-rank test, theoretical median = 0).Population data are illustrated as mean ± SEM.

In sum, FFI onto L5 PNs was triggered by L2/3 PNs and could be strongly potentiated by postsynaptic firing activity alone. This form of plasticity was mediated by NO signaling and is thus consistent with the involvement of PV cells ([@bib42]). Plasticity of FFI resulted in prolonged changes of the E/I balance in L5 PNs.

LTPi-FFI Induces a Divisive Gain Modulation of PN Output Spiking {#sec2.2}
----------------------------------------------------------------

Inhibitory transmission alters the computations performed by PNs, strongly affecting their input/output (I/O) relationship ([@bib16], [@bib67]), but how plastic changes of inhibitory synaptic strength could modulate I/O is unknown. To address this question, we measured the I/O relationship before and after inducing LTPi-FFI in L5 PNs by postsynaptic AP bursts ([Figures 2](#fig2){ref-type="fig"}A--2C). To reliably evoke presynaptic spike trains, we expressed the fast light-sensitive opsin ChETA in L2/3 PNs via *in utero* electroporation. We then tested its effectiveness by recording spike activity in ChETA-expressing L2/3 PNs in response to trains of brief blue light pulses over a range of different frequencies ([Figures S4](#mmc1){ref-type="supplementary-material"}A--S4C). *In utero* electroporated L2/3 PNs linearly followed ChETA stimulation at frequencies up to 30 Hz ([Figure S2](#mmc1){ref-type="supplementary-material"}). The effect of LTPi-FFI on synaptic integration was then investigated by stimulating L2/3 ChETA(+) neurons with 10 pulses of light at different frequencies and measuring the mean output firing rate of L5 PNs ([Figures 2](#fig2){ref-type="fig"}D and 2E) before and after inducing LTPi-FFI. We observed an overall difference between baseline and LTPi-FFI with a major effect on higher frequencies between 15 and 30 Hz (F~(11,165)~ = 13.76, p \< 0.0001; p \< 0.05 for 15 Hz, p \< 0.0001 for 20 Hz, and p \< 0.05 for 30 Hz; Friedman test followed by Dunn's multiple comparison test; [Figure 2](#fig2){ref-type="fig"}E). The I/O relationship was quantified by fitting the data to Hill-like equations, which allowed us to extract slope and offset of I/O curves ([@bib53], [@bib63]). We observed a shift in the input frequency required to achieve half maximum of the output frequency, although this was not significantly different (Δ~Offset~ = 1.56 ± 1.3 Hz, n = 14, p = 0.8, Wilcoxon matched-pairs signed-rank test, theoretical median = 0; [Figure 2](#fig2){ref-type="fig"}F, light gray bar). Despite an absence of significant changes in the subtractive shift behavior of the I/O curve, we observed a significant change in the slope (Δ~Gain~ = −36.29% ± 11.79%, n = 14, p = 0.01, Wilcoxon matched-pairs signed-rank test, theoretical median = 0; [Figure 2](#fig2){ref-type="fig"}F, dark gray bar). Importantly, in the minority of L5 PNs that failed to express LTPi-FFI ([Figure S1](#mmc1){ref-type="supplementary-material"}I), we did not observe any change in gain modulation ([Figure S2](#mmc1){ref-type="supplementary-material"}D). Thus, these results indicate that LTPi-FFI has as an almost purely divisive effect on the PN I/O relationship ([Figures 2](#fig2){ref-type="fig"}E and F).

LTPi-FFI Affects the Flow of Information across Cortical Layers {#sec2.3}
---------------------------------------------------------------

Sustained firing of L2/3 PNs induces lateral suppression of superficial cortical layers and simultaneously activates L5 PNs within the same column. This latter phenomenon is known as feedforward facilitation (FFF) and is due to differences in E/I ratios in different cortical layers ([@bib1]). Given the abovementioned effects of LTPi-FFI on the E/I ratio and gain modulation of L5 PN output spiking, we hypothesized that L2/3-induced facilitation of deep cortical layer PNs could be affected by potentiation of local perisomatic inhibitory transmission. FFF can be induced in cortical slices by depolarizing ChR2+ L2/3 PNs with a 1-s-long ramp of blue light while simultaneously depolarizing L5 PNs ([@bib1]). L5 PNs were depolarized with 1-s-long current injections to trigger AP firing (average rate, 5.4 ± 0.29 Hz, range, 2--8 Hz; n = 23; [Figures 3](#fig3){ref-type="fig"}A and 3C, baseline \[Bsl\], black trace, and gray bar), which was significantly increased by simultaneous photostimulation of L2/3 (average rate, 8.26 ± 0.35 Hz; range, 5--12 Hz; n = 23, [Figures 3](#fig3){ref-type="fig"}A and 3C, blue trace and bar, Bsl_light; F(3,66) = 35.83, p \< 0.0001; p \< 0.0001 for Bsl versus Bsl_light, Friedman test followed by Dunn's multiple comparison test). LTPi-FFI did not change the Bsl spike rate of L5 PNs (average rate, 5.873 ± 0.35 Hz; range, 3--8 Hz; n = 23; [Figures 3](#fig3){ref-type="fig"}B and 3C, black trace and dotted bar, LTPi), ruling out the possibility that postsynaptic burst firing alone altered L5 PN excitability. Remarkably, when the same experiment was repeated after inducing LTPi-FFI, facilitation of L5 PN excitability, triggered by L2/3 activation, was largely decreased (average rate, 6.22 ± 0.41 Hz; range, 2--10 Hz; n = 23; [Figures 3](#fig3){ref-type="fig"}B and 3C, red trace and red bar, LTPi_light; F(3,66) = 35.83, p \< 0.0001; p = 0.7218 for LTPi versus LTPi_light, Friedman test followed by Dunn's multiple comparisons test).Figure 3LTPi-FFI Affects the Flow Information across Cortical Layers(A) Response of a L5 PN to a depolarizing current injection (+70 pA) under Bsl conditions without (left trace, Bsl) and with light stimulation (right trace, Bsl_light), inducing feedforward facilitation (FFF).(B) The same cell as in (A) upon LTPi-inducing AP bursts. APs were clipped for display purposes.(C) Average population data of FFF during Bsl conditions in the absence and presence of light (circles) and upon LTPI-FFI (diamonds, ^∗∗∗∗^p \< 0.0001) and Bsl_light and LTPi-FFI light conditions (^∗∗^p \< 0.001, Friedman test).(D) Summary scheme of the effect of LTPi in preventing FFF across cortical layers. In the absence of LTPi, L2/3 activity efficiently spreads to L5 (left). LTPi prevents this prominent information transfer from L2/3 to L5 PNs.Population data are illustrated as mean ± SEM.

Altogether, these results indicate that the cell-autonomous long-term strengthening of perisomatic inhibition can reduce or completely abolish coordination of cortical activity across cortical layers ([Figure 3](#fig3){ref-type="fig"}D).

LTPi-FFI Induces a Shift in Temporal Association of L5 PN Firing during Photo-Induced Rhythmic Activity {#sec2.4}
-------------------------------------------------------------------------------------------------------

Perisomatic inhibition, particularly from PV basket cells, drives large populations of PNs, entraining them to network oscillations in the β-γ-frequency range ([@bib15], [@bib17], [@bib68]). Therefore, we hypothesized that changes in perisomatic inhibitory strength might affect the temporal association of PN firing with ongoing γ-activity. A well-characterized protocol to induce γ-oscillations in cortical slices is optogenetic stimulation with ramps of blue light, resulting in prolonged firing of ChR2-expressing L2/3 PNs ([@bib1], [@bib21]). This protocol evoked robust γ-oscillations that depend on both GABAergic and glutamatergic synaptic transmission ([@bib1], [@bib61], [@bib65]). Importantly, photo-induced γ-activity is faithfully transmitted vertically to L5 PNs and PV cells belonging to the same cortical column ([@bib1], [@bib59], [@bib21]).

We induced oscillatory activity by photostimulation of ChR2+ L2/3 PNs while simultaneously recording inhibitory postsynaptic currents (IPSCs) and excitatory postsynaptic currents (EPSCs) in voltage clamp in ChR2-negative PNs in L2/3 and L5, respectively (average amplitude of IPSCs, 703.9 ± 82.45 pA; average amplitude of EPSCs, 266.8 ± 56.87 pA; n = 15; [Figures 4](#fig4){ref-type="fig"}A and 4B). Photostimulation induced robust oscillations in the γ-frequency range, involving both IPSCs and EPSCs, in L2/3 and L5 PNs, respectively (frequency IPSCs, 26.04 ± 0.78 Hz; frequency EPSCs, 27.26 ± 1.37 Hz; n = 15, p = 0.45, Wilcoxon matched-pairs signed-rank test; [Figure 4](#fig4){ref-type="fig"}C). We then recorded L5 PNs in current-clamp mode to analyze how their spikes were temporally associated with ongoing IPSC rhythmic activity, recorded in voltage-clamp in L2/3 PNS ([Figures 4](#fig4){ref-type="fig"}D and 4E). During baseline (Bsl), L5 PNs were tuned to γ-oscillations and discharged, on average, ∼2 ms before the peak of the γ-cycle (spike time to oscillation peak during Bsl, −2.44 ± 0.44 ms, n = 9; [Figures 4](#fig4){ref-type="fig"}E and 4F). Remarkably, after inducing LTPi-FFI with bursts of APs, L5 PN firing was significantly anticipated (spike time to oscillations after AP bursts, −4.22 ± 0.52 ms, n = 9, p = 0.0313, Wilcoxon matched-pairs signed-rank test; [Figures 4](#fig4){ref-type="fig"}E and 4F). We restricted our analysis to a pool of PNs exhibiting sufficient spiking activity (\>1 Hz) following LTPi. In these cells, the observed change in spike timing was not accompanied by changes in firing rate or resting membrane potential ([Figures 4](#fig4){ref-type="fig"}G and 4H, firing rate: 6.3 ± 0.61 Hz versus 4.8 ± 0.58 Hz, Bsl versus LTPi, respectively; paired t test, p = 0.0804; membrane potential: p = 0.49 for Bsl versus LTPi, paired t test). Importantly, however, when the firing frequency was analyzed for all PNs exhibiting LTPi, the firing frequency was significantly reduced in response to LTPi-inducing bursts (5.643 ± 0.7 versus 3.88 ± 0.67 Hz, Bsl versus LTPi, respectively; p = 0.0266, paired t test; [Figure 4](#fig4){ref-type="fig"}G). Interestingly, in cells in which LTPi-FFI could not be induced, we did not observe a change in timing association (Bsl time to oscillation peak, −5.33 ± 1.12 ms; after bursts, −5.67 ± 1.2 ms, n = 6; p = 0.36, paired t test; [Figure 4](#fig4){ref-type="fig"}F). Intriguingly, we found that PNs that did not express LTPi displayed a significantly longer time to oscillation peak (p = 0.0269, Kruskal-Wallis test followed by Dunn's multiple comparisons test), which was similar to that of PNs exhibiting plasticity (p = 1.0).Figure 4LTPi-FFI Induces a Shift in Temporal Association of L5 PN Firing during Photo-Induced Rhythmic Activity(A) Scheme of the recording configuration.(B) Representative traces during a light ramp (blue, 1-s duration) protocol inducing rhythmic activity of L2/3 IPSCs (blue trace) and L5 EPSCs (red trace).(C) Power spectra of the EPSCs and IPSCs of the cells shown in (B).(D) Representative firing activity of a L5 PN during γ-activity in L2/3 before (Bsl, black trace) and after induction of LTPi-FFI (red trace).(E) Spike probability during a cycle (blue trace) of the same L5 PN before (black) and after LTPi-FFI (red).(F) Left: representative traces illustrating the association of a spike of a L5 PN (red trace) with the peak of oscillating IPSCs from a L2/3 PN (black trace). Right: population data of the time to oscillation peak in cells that underwent LTPI-FFI (circles) and cells in which the postsynaptic bursts had no effect (diamonds). ^∗^p \< 0.05, Wilcoxon matched-pairs signed-rank test (Bsl versus LTPi) and ^∗^p \< 0.05, Kruskal-Wallis test followed by Dunn's multiple comparisons test (Bsl versus Bsl).(G) Firing rates of LTPI-FFI (circles) and no LTPi cells (diamonds). Black and red dots (left), and black and gray diamonds (right) refer to PNs used for temporal association analysis. Cells indicated with blue and/or pink symbols had to be removed from the temporal association analysis because of spike sparsity. ^∗^p \< 0.05, paired t test (Bsl versus LTPi).(H) Average population data of the membrane resting potentials of LTPI-FFI (circles) and no LTPi cells (diamonds).Population data are illustrated as mean ± SEM.

These results indicate that potentiation of feedforward perisomatic inhibition alters the temporal association of L5 PNs during γ-oscillations.

Bursts of APs Decrease the Firing Rate in a Subset of L5 PNs *In Vivo* {#sec2.5}
----------------------------------------------------------------------

Do LTPi-inducing bursts also decrease L5 PN spikes *in vivo*? We performed whole-cell recordings from L5 PNs in the barrel cortex of anesthetized mice ([STAR Methods](#sec4){ref-type="sec"}) during spontaneous activity and in the presence of sensory stimulation induced by air puffs to the contralateral whisker pads. Location and PN identity were confirmed by anatomy in some experiments. Notably, 100% of recordings were in L5, as confirmed by labeling L4 in scnn1a-cre crossed with a td-Tomato reporter line ([Figure S3](#mmc1){ref-type="supplementary-material"}). PNs could be separated into two populations according to their change in firing rate following AP bursts (5 APs at 100 Hz, repeated 15 times every 10 s). Type A cells displayed an increased firing rate upon LTPi-inducing bursts during spontaneous activity and in the presence of sensory stimulation (1.37 ± 0.28 versus 2.43 ± 0.52 Hz, pre-bursts versus after bursts, respectively; n = 10, p = 0.002, Wilcoxon matched-pairs signed-rank test; [Figures 5](#fig5){ref-type="fig"}A and 5B; [Figure S4](#mmc1){ref-type="supplementary-material"}). Conversely, type B cells exhibited a marked decrease in firing rate (2.4 ± 0.56 versus 1.04 ± 0.26 Hz, pre-bursts versus after bursts, respectively; n = 16, p \< 0.0001, Wilcoxon matched-pairs signed-rank test; [Figures 5](#fig5){ref-type="fig"}C, 5D and [4](#fig4){ref-type="fig"}). Importantly, these changes were not associated with significant variations in membrane potential (type A: −57.9 ± 2.8 versus −55.4 ± 2.9 mV, pre-bursts versus after bursts; n = 10, p = 0.08, Wilcoxon matched-pairs signed-rank test; type B: −58.09 ± 2.9 versus −56.31 ± 2.0 mV, pre-bursts versus after bursts; n = 16, p = 0.25, Wilcoxon matched-pairs signed-rank test). The depth of patched neurons *in vivo* also did not differ between type A and B neurons (666.7 ± 21.08 μm depth versus 714.4 ± 21.35 μm depth, type A versus type B, respectively; p = 0.2440, Mann-Whitney test). Moreover, the Bsl spike frequency of type A and B neurons was not significantly different (p = 0.3767, Kruskal-Wallis test followed by Dunn's multiple comparisons test). Type A and B accounted for 38.5% and 61.5% of recorded PNs.Figure 5Bursts of APs Decrease the Firing Rate in a Subset of L5 PNs *In Vivo*(A) Left: scheme of the recording configuration. Right: representative experiment of a type A PN (top, raster plot; bottom, representative current-clamp trace). Air puff stimulation is indicated as a gray rectangle during a pre- and post-burst period.(B) Plot of the mean firing rate from individual PNs before (x axis) versus 5--10 min after postsynaptic bursts (y axis). Type A L5 PNs displayed a small but significant increase in firing rate after bursts. A dotted line indicates unitary values (no change). Blue-filled symbols refer to PNs receiving air puff stimulation; open circles refer to PNs, which did not receive air puff stimulation. Note that these are different cells from different recordings. ^∗∗^p \< 0.01, Wilcoxon matched-pairs signed-rank test.(C and D) The same as in (A) and (B) but for type B cells (C and D, respectively), in which postsynaptic bursts induced a decrease in firing rate under both conditions. ^∗∗∗^p \< 0.0001, Wilcoxon matched-pairs signed-rank test.Population data are illustrated as mean ± SEM.

These results indicate that, *in vivo*, a prominent fraction of PNs respond to repetitive AP burst firing with a decreased firing rate, possibly through potentiation of perisomatic inhibition, similarly to our slice results.

A Reduced Firing Rate Is Associated with Increased Tuning of PN Spiking Activity with γ-Oscillations *In Vivo* {#sec2.6}
--------------------------------------------------------------------------------------------------------------

We found that LTPi is associated with altered temporal association of PN firing to photo-induced γ-oscillations in acute cortical slices ([Figure 4](#fig4){ref-type="fig"}). Photo-activated γ-oscillations are a very helpful strategy to study network activity in slice. However, their activation is somewhat artificial, and optogenetic hyper-synchrony can affect oscillation temporal dynamics. We therefore examined whether postsynaptic bursting activity changes the temporal association of L5 PN firing with naturally occurring rhythmic activity *in vivo*. Previous evidence indicates that cortical PNs are poorly coupled to spontaneous and sensory evoked γ-activity ([@bib55]). Could increases in perisomatic inhibition facilitate tuning of PN spiking with ongoing network oscillations? *In vivo* whole-cell recordings of L5 PNs were coupled to local field potential (LFP) recordings obtained with a separate electrode ([Figure S5](#mmc1){ref-type="supplementary-material"}; [STAR Methods](#sec4){ref-type="sec"}). We analyzed the relationship of spike probability of L5 PNs to γ-activity embedded in the LFP before and after bursting, in the presence and absence of sensory stimulation. As reported previously in the visual cortex ([@bib55]), PN spiking activity was poorly tuned in the presence and absence of whisker stimulation ([Figure 6](#fig6){ref-type="fig"}; [Figure S6](#mmc1){ref-type="supplementary-material"}). On average, spike distributions did not reveal a significant phase preference under both conditions ([Figure 6](#fig6){ref-type="fig"}; [Figure S6](#mmc1){ref-type="supplementary-material"}). However, in 8 of 11 cells recorded during spontaneous activity (n = 10 mice) and in 11 of 15 cells recorded during air puff stimulation (n = 13 mice), we observed a slight but significant increase in spike tuning with γ-activity. This was indicated by nonrandom spike-phase distributions crossing the Monte Carlo simulation threshold ([STAR Methods](#sec4){ref-type="sec"}; [Figure 6](#fig6){ref-type="fig"}; [Figure S6](#mmc1){ref-type="supplementary-material"}).Figure 6A Reduced Firing Rate Is Associated with Increased Tuning of PN Spiking Activity with γ-Oscillations *In Vivo*(A) Top: current-clamp traces showing firing activity in a type A L5 PN (top trace) and filtered LFP (between 20 and 100 Hz, bottom trace). The red bar indicates air puff stimulation.(B) Polar plots displaying the circular distribution of APs with γ-phases in the same type A neuron before (left, pre-bursts, black trace) and after bursts (right, post-bursts, blue trace). Gray lines refer to confidence intervals.(C) Population spike distribution across the phase before (black trace) and after bursts (blue trace). Note the point of maximum (max) spike probability (red arrow).(D) Summary plot of max spike probability extracted from each cell before (light blue circles) and after bursts (dark blue circles).(E--H) The same as in (A)--(D) but for type B cells (E, F, G, and H, respectively). In these PNs, postsynaptic bursts induced a marked increase in max spike probability. ^∗∗^p \< 0.01, Wilcoxon matched-pairs signed-rank test.Population data are illustrated as mean ± SEM.

Interestingly, in sensory evoked conditions, although type A cells did not exhibit any significant change in maximum (max) spike probability ([Figures 6](#fig6){ref-type="fig"}A--6C), type B displayed a significant increase in max spike probability after bursts ([Figures 6](#fig6){ref-type="fig"}D--6F; type A: 0.11 ± 0.016 versus 0.092 ± 0.011, Bsl versus after bursts, n = 6, p = 0.22, Wilcoxon matched-pairs signed-rank test; type B: 0.084 ± 0.007 versus 0.113 ± 0.007, Bsl versus after bursts, n = 9, p = 0.0039, Wilcoxon matched-pairs signed-rank test). Accordingly, after AP bursts, the fraction of type B neurons tuned to γ-activity increased (5 of 9 versus 8 of 9, Bsl versus after bursts, respectively). This was also reflected by an increased delta pairwise phase consistency (PPC), which is a measure of phase locking, independent of spike rate ([@bib55], [@bib71]; [Figure S6](#mmc1){ref-type="supplementary-material"}H).

These results indicate that type B L5 PNs increase their tuning to γ-oscillations both during spontaneous and sensory evoked activity. This suggests that potentiation of inhibitory transmission likely promotes orchestrated activity of single PNs *in vivo*.

A Computational Model Reveals Two Separate Effects of LTPI-FFI {#sec2.7}
--------------------------------------------------------------

Our *in vitro* and *in vivo* results ([Figures 4](#fig4){ref-type="fig"} and [6](#fig6){ref-type="fig"}) suggest that LTPi alters the temporal association of PN firing during γ-oscillations. This prompts the question whether this is due to the actual strength of PN perisomatic inhibition, which sets the actual E/I level. To address this question, we examined a computational model of the FFI circuit between L2/3 and L5 ([Figure 7](#fig7){ref-type="fig"}A). In this model, a L5 PN was represented by an integrate-and-fire neuron, which was driven by L2/3 oscillatory activity (30 Hz) through a composite postsynaptic current (PSC, [Figure 7](#fig7){ref-type="fig"}B). As determined experimentally in response to brief L2/3 optical stimulation ([Figures 1](#fig1){ref-type="fig"}A--1C) or synaptic recordings ([@bib41]), the composite PSC consisted of an early excitatory component, followed, after a short delay, by an inhibitory component corresponding to FFI ([STAR Methods](#sec4){ref-type="sec"}). The amplitude of the inhibitory component was systematically varied to investigate the effects of LTPi-FFI in the model. Additional inputs to the L5 PN were modeled as background noise.Figure 7A Computational Model Reveals Two Separate Effects of LTPI-FFI(A) Schematic of the model.(B) Composite postsynaptic current received by the model L5 PN in response to a brief L2/3 activation for four increasing values of inhibitory strength.(C) Activity in the model L5 PN in response to oscillatory L2/3 inputs for four increasing values of inhibitory strength (left to right). From top to bottom: L2/3 oscillation, membrane potential traces in the model L5 PN, rastergram of model L5 PN APs over 300 repeats of the input, and histogram of the average activity computed over 3,000 repeats of the input.(D and E) Phase and amplitude of the average L5PN activity as function of the I/E ratio. Dots, results of cosine fits to simulation data (illustrated in the bottom panels in C); continuous line, theoretical prediction based on the shape of the composite PSC ([STAR Methods](#sec4){ref-type="sec"}). In (D), the phase is determined with respect to L2/3 input. In (E), the amplitude was normalized to 1 for an I/E ratio of 0.05. The theoretical prediction is based on a linear approximation. This approximation fails at large values of inhibition, as seen in the deviations between predicted and measured amplitude (E). Non-linear effects lead to an increase in peri-oscillation time histogram (POTH) amplitude and, therefore, firing precision, with respect to the linear prediction.

L2/3 oscillatory inputs filtered through the composite PSC led to temporal modulation of the activity in the model L5 PN ([Figure 7](#fig7){ref-type="fig"}C). To determine the temporal relationship between L5 and L2/3 activity, we computed the histograms of temporal activity in L5 and extracted its phase and amplitude by fitting cosine functions ([Figure 7](#fig7){ref-type="fig"}C). The phase of the histogram quantifies the average timing of L5 activity with respect to L2/3 oscillations, whereas the amplitude of the modulation represents the precision of L5 spikes. A mathematical analysis of the circuit dynamics predicted that the phase and amplitude of L5 PN activity are determined by the phase and amplitude of the Fourier transform of the composite PSC, evaluated at the frequency of 30 Hz, which corresponds to L2/3 oscillatory inputs ([Figure 7](#fig7){ref-type="fig"}B; [STAR Methods](#sec4){ref-type="sec"}). A comparison with numerical simulations confirmed this prediction ([Figures 7](#fig7){ref-type="fig"}D and 7E).

The circuit model revealed that increasing the amplitude of the inhibitory component of the composite PSC led to two separate effects. The first effect is that increasing inhibition at moderate strengths shifts the average timing of L5 activity with respect to L2/3 oscillations but leaves the precision of the activity essentially unchanged. This timing shift is consistent with the experimental measurements reported in [Figure 4](#fig4){ref-type="fig"} (note that, in [Figure 4](#fig4){ref-type="fig"}, the absolute timing was quantified with respect to L2/3 IPSCs, whereas in the model, the timing is quantified with respect to L2/3 average firing activity). This shift can be understood by considering the histograms of activity in [Figure 7](#fig7){ref-type="fig"}C. If the composite PSC is dominantly excitatory ([Figure 7](#fig7){ref-type="fig"}C, leftmost column), then a peak in L2/3 activity ([Figure 7](#fig7){ref-type="fig"}C, top row) is followed by a peak in the L5 activity histogram ([Figure 7](#fig7){ref-type="fig"}C, bottom row). If the compound PSC is dominantly inhibitory ([Figure 7](#fig7){ref-type="fig"}C, rightmost column), a peak in L2/3 activity ([Figure 7](#fig7){ref-type="fig"}C, top row) is instead followed by a trough in the L5 activity histogram ([Figure 7](#fig7){ref-type="fig"}C, bottom row). In between these two extremes, increasing inhibition progressively shifts the trough in L5 activity earlier with respect to the peak of L2/3 inputs ([Figure 7](#fig7){ref-type="fig"}C, bottom row, dashed lines show the peak of L2/3 inputs).

The model shows that the range of possible timing shift is limited because the phase shift saturates as inhibition is increased. Interestingly, however, large increases in inhibition lead to a second, separate effect in which the timing of L5 activity does not change anymore, but its precision increases with stronger inhibition, as shown previously ([@bib42]). Accordingly, this effect is consistent with increased tuning *in vivo*, shown in [Figure 6](#fig6){ref-type="fig"}. The difference in the overall amplitude of tuning to γ-oscillations between slices and *in vivo* conditions could be ascribed to different background noise levels, typical of slices versus *in vivo* conditions ([Figure S7](#mmc1){ref-type="supplementary-material"}). However, both phase shift and increased tuning were present over a very broad range of background noise level ([Figure S7](#mmc1){ref-type="supplementary-material"}).

Discussion {#sec3}
==========

In this study, we provide evidence that bursting activity of L5 PNs can potentiate incoming FFI recruited by a descending excitatory pathway, which is one of the most prominent in the neocortex. This plasticity likely originated from PV basket cells, strongly modulated PN gain, affected coordinated activity across cortical layers, and altered the temporal association of PNs with γ-oscillations both *in vitro* and *in vivo*.

In principle, FFI can result from long-range recruitment of different interneuron subclasses ([@bib59]). However, potentiation of FFI likely involves perisomatic inhibition originating almost exclusively from PV basket cells. Indeed, here we provide evidence that optogenetic activation of L2/3 PNs efficiently recruits L5 PV basket cells ([Figure S2](#mmc1){ref-type="supplementary-material"}). Moreover, LTPi-FFI was completely abolished by the canonical NO receptor inhibitor ODQ, which we previously showed to be selective for PV cell-mediated LTPi, and did not involve dendrite-targeting SST interneurons ([@bib42]). Moreover, SST interneurons are recruited more efficiently in response to prolonged trains because of strongly facilitating glutamatergic synapses ([@bib35], [@bib66], [@bib75]). Although we cannot exclude that other interneuron types might be involved in LTPi-FFI, we believe that this plasticity relies on the strength of PV cell synapses. Accordingly, PV basket cells have been demonstrated to be the most prominent, if not exclusive, providers of perisomatic inhibition onto large deep-layer PNs of the somatosensory cortex ([@bib8]) and to be the main source of evoked inhibition onto L5 PNs ([@bib37], [@bib50]).

Our data indicate that postsynaptic firing activity alone can effectively alter the E/I ratio of a prominent input pathway, impinging L5 PNs. A tight balance between excitation and inhibition is believed to guarantee proper functioning of neural circuits ([@bib33], [@bib49], [@bib32], [@bib54], [@bib27], [@bib76], [@bib24], [@bib31]). However, loosening of the E/I lock is necessary for sensory processing and refinement of sensory maps ([@bib27]). Moreover, it has been shown that the E/I ratio is different across individual cortical principal cells, depending on the specific layer ([@bib1]) and the intrinsic activity state of each PN ([@bib78]), demonstrating that the tight lock of the E/I ratio can be disrupted by perturbing pyramidal cell activity. Therefore, plasticity of GABAergic perisomatic FFI can be a relevant mechanism to tune PN sensitivity to process sensory information within the cortical column.

LTPi-FFI results in a strong divisive output modulation of the I/O curve because it significantly changed its slope without altering the rheobase. From a computational perspective, this means that, for rate-coded neuronal signaling, the modulation operated by LTPi-FFI results in a change of PN gain that affects the dynamic range of its response, likely preventing saturation of firing ([@bib67]). However, whether inhibition plays an additive or multiplicative role has been debated, depending on the location (perisomatic versus dendritic) and strength of the inhibitory response as well as on the morphological complexity of the postsynaptic neuron ([@bib43], [@bib67]). In particular, perisomatic PV cell-mediated modulation of orientation tuning of visual cortical PNs has been shown to have both an additive and multiplicative effect with strong modulation of PN gain during sensory processing ([@bib5], [@bib39], [@bib77]). Here we show that cell-autonomous potentiation of perisomatic inhibition has a strong effect on PN gain. This was not due to changes in intrinsic excitability (e.g., changes in membrane potential or resistance) but due to a selective increase in FFI. Therefore, by inducing retrograde potentiation of perisomatic inhibition, burst firing of PNs altered their functional network connectivity and, thus, the dynamic pattern of their activation through divisive modulation of their gain.

Cortical PNs are known to respond to sensory stimuli using a sparse population coding ([@bib56]). In this regime, neurons act as coincidence detectors of temporally correlated input; changes in the gain of the I/O curve of PNs might shape the time window in which inputs can be integrated to generate a spike ([@bib42], [@bib60]). It is therefore likely that burst firing-induced potentiation of inhibition controls the temporal properties of signals propagating through the network.

The dynamic modulation of PN firing, induced by plasticity of FFI had a powerful effect on columnar integration of activity across cortical layers. Indeed, activation of L2/3 PNs has been shown to increase the firing of PNs in L5 of the same column ([@bib1]; but see [@bib59]). LTP of FFI strongly diminished (in fact, on average, it abolished) L2/3-dependent activation of L5 PNs, affecting the flow of neuronal communication across cortical layers.

Here we show that induction of long-term plasticity of feedforward perisomatic inhibition shifted the temporal association of L5 PN spikes with photo-induced γ-oscillations. This effect was due to increased inhibition because glutamatergic neurotransmission and passive properties were not altered by postsynaptic burst firing.

Photo-activated γ-oscillations are a very useful tool to dissect the cellular and synaptic mechanisms underlying network synchronization ([@bib1], [@bib62]). However, this approach suffers from some limitations; namely, hyper-synchronous activity induced by simultaneous prolonged activation of a large number of neurons in a reduced preparation. We therefore investigated whether LTPi-inducing trains could affect spike timing relative to endogenous γ-activity in a more intact preparation. We found that our *in vitro* and *in vivo* results converged: burst firing induced a decrease in firing rate in a prominent fraction of L5 PNs and improved tuning during spontaneous and sensory evoked γ-activity. The modulation of spike association to network oscillations both *in vitro* and *in vivo* we show here is in line with the known function of perisomatic inhibition (particularly from PV basket cells) to entrain PNs during γ-activity ([@bib7], [@bib74], [@bib14], [@bib17], [@bib33], [@bib68]; but see [@bib71]).

Moreover, both *in vitro* and *in vivo* we found a similar low percentage of neurons that did not display LTPi and did not exhibit decreased firing rates, respectively. We termed these neurons type A in our *in vivo* recordings, and we speculate that they correspond to the fraction of PNs not expressing LTPi *in vitro*. In these PNs, coupling with network activity was not affected by burst firing. It is currently unclear whether these PNs belong to different cell types. Indeed, it has been demonstrated that thick-tufted, large, cortico-fugal PNs of the prefrontal cortex are preferentially innervated by PV cells compared with more slender, thin-tufted, cortico-cortical PNs ([@bib2], [@bib40]). Alternatively, in this minority of PNs, the lack of LTPi expression could be due to already potentiated (and, hence, saturated) inhibition. This could explain the consistent negative shift of their spike times *in vitro* compared with PNs that did undergo LTPi. Moreover, bi-directional, activity-dependent changes in firing were reported in L5 PNs of the rat neocortex ([@bib44]). Future experiments will be necessary to unequivocally determine the detailed morpho-functional properties of the small percentage of PNs lacking expression of GABAergic plasticity.

Importantly, our computational model revealed that increasing inhibition in the feedforward circuit induced two separate effects for different levels of the E/I ratio. If excitation dominates, then increasing inhibition modulates the timing of L5 spikes but not their precision with respect to L2/3 oscillations. This is indeed what we found during photo-activated γ-oscillations in neocortical slices. If, however, inhibition becomes stronger than excitation, then further increases in inhibition improve the precision of L5 spikes (as found in [@bib42]) but do not additionally affect their timing. This could explain the improvement in tuning we detect *in vivo*, where assessing the actual E/I ratio is technically challenging and prone to strong biases. Modifying inhibition strength within different ranges may therefore lead to different functional consequences. The strength of inhibitory connections is known to control coherence in recurrent networks ([@bib6], [@bib9], [@bib11]). Here, in contrast, we demonstrated its specific role in temporal organization of activity in FFI circuits independent of background noise level.

Therefore, plasticity of feedforward perisomatic inhibition by governing spike-time association of single PNs with ongoing γ-activity might be responsible for shifting their participation to distinct cell assemblies, reconfiguring the local network ([@bib52]).

In conclusion, LTPi of feedforward inhibition can be a simple mechanism modulating the functional connectivity of single PNs, largely influencing cortical networks and subnetworks. LTPi of PV cell-mediated FFI in sensory cortices could be a fundamental dynamic property of cortical networks, providing the basis of diverse cognitive functions such as sensory perception and attention.

STAR★Methods {#sec4}
============

Key Resources Table {#sec4.1}
-------------------

REAGENT or RESOURCESOURCEIDENTIFIER**Chemicals, Peptides, and Recombinant Proteins**ODQR&D Systems Europe0880SR 95531 hydrobromideR&D Systems Europe15753BiocytinSigmaB4261Bupivacaine hydrochlorideSigmaB1160000Buprenorphine (Buprecare)CentravetBUP002AtropineICM Animal House PharmacyN/ADexametosone (DEXAFORT)CentravetDEX211Fast GreenSigmaF7252**Critical Commercial Assays**DAB substrate kit for peroxidaseVECTASTAIN EliteSK-4100ABC kitVECTASTAIN ElitePK-6100**Experimental Models: Organisms/Strains**Mice: C57BL6/6JJanvierSC-C57J-FMice: Ai14JacksonRRID:IMSR_JAX:007914Mice: Scnn1a-Tg3-CreJacksonRRID:IMSR_JAX:009613Mice: SwissJanvierSN-SWISS-F**Recombinant DNA**p-Lenti-CaMKIIa-ChETA-EYFPAddgeneRRID:Addgene_26967pCAG-mRFPAddgeneRRID:Addgene_28311pCAG-ChR2-VenusAddgeneRRID:Addgene_15753p-ChETA-EYFPiVector (ICM)Subcloned from p-Lenti-CaMKIIa-ChETA-EYFP**Software and Algorithms**pClamp 10.3Molecular Devices<https://www.moleculardevices.com/products/axon-patch-clamp-system/acquisition-and-analysis-software/pclamp-software-suite>IGOR PRO 5.0Wavemetrics<https://www.wavemetrics.com/>PrismPrism-GraphPad<https://www.graphpad.com/scientific-software/prism/>MATLABMathworks<https://www.mathworks.com/>Detector -- custom written softwareJ. R. Huguenard, Stanford UniversityJ. R. Huguenard, Stanford University<https://hlab.stanford.edu/wdetecta.php>**Deposited Data**DatasetsMendeley Data<https://data.mendeley.com/datasets/75t97ncb24/2>Theoretical ModelGitHub<https://github.com/MathildeBigot>

Lead Contact and Materials Availability {#sec4.2}
---------------------------------------

Further information and requests for resources and reagents should be directed to and will be fulfilled by the Lead Contact, Alberto Bacci (<alberto.bacci@icm-institute.org>). All unique reagents generated in this study are available from the Lead Contact without restriction.

Experimental Model and Subject Details {#sec4.3}
--------------------------------------

Experimental procedures followed national and European (2010/63/EU) guidelines, and have been approved by the authors' institutional review boards and national authorities. All efforts were made to minimize suffering and reduce the number of animals. Experiments were performed on 15- to 28-day-old C57BL/6 wild-type mice (Janvier Labs, France) or Scnn1a-Tg3-Cre mice expressing Cre recombinase under the control of the *Scnn1a* (sodium channel, nonvoltage-gated 1 alpha; Jackson) crossed with Ai9 mice, which have a *loxP*-flanked STOP cassette preventing transcription of a CAG promoter-driven red fluorescent protein variant (tdTomato) - all inserted into the *Gt(ROSA)26Sor* locus (Jackson). Mice used in this study were of both sexes.

Method Details {#sec4.4}
--------------

### In utero electroporation {#sec4.4.1}

Timed-pregnant C57BL/6 wild-type female mice (15.5 days postcoitum) were anaesthetized with 1%--2% isoflurane. The abdomen was cleaned with 70% ethanol and swabbed with betadine. Buprenorphine (0.05 mg/kg) was administered subcutaneously for preoperative analgesia and local anesthetic bupivacaine (2.5mg/kg) was injected between the skin and the abdomen 5 min before incision. A midline ventral laparotomy (∼2 cm) was performed, and the uterus gently exposed and moistened with PBS, pre-warmed at 37°C. Using glass beveled capillaries, DNA plasmids mixed in saline (PBS) solution and 0.025% Fast Green (Sigma) were injected through the uterine wall into the lateral ventricle of each embryo. Embryos were injected with pCAG-mRFP (0.8 μg/μl) ([@bib46]) (Addgene \#28311) plasmid DNA mixed with either pCAG-ChR2-Venus ([@bib57]) (Addgene \#15753) or pCAG-ChETA-EYFP (1.5 μg/μL). ChETA-EYFP was subcloned from p-Lenti-CaMKIIa-ChETA-EYFP ([@bib30]) (Addgene \#26967) into the backbone of pCAG-ChR2. After each injection, the embryos were moistened with PBS. DNA was electroporated via 5 square electrical pulses of 40 V amplitude and 50 ms duration through forceps-type circular electrodes positioned at 0° angle with respect to the rostral-caudal axis of the head of the embryos. After electroporation, the uterus was placed back into the peritoneal cavity and moistened with PBS. The abdomen and skin were then sutured and the latter cleaned with betadine. The procedure typically lasted maximum 40 min starting from anesthesia induction. Pups were born by natural birth and placed with a Swiss foster mother (timed-pregnant at the same time as C57BL/6 females) after they were screened for location and strength of transfection by trans-cranial epifluorescence under a fluorescence stereoscope.

### *In Vitro* Slice Preparation and Electrophysiology {#sec4.4.2}

Coronal slices (400-μm-thick) from somatosensory cortex were obtained from 15- to 28-day-old C57BL/6 mice. Animals were deeply anesthetized with isofluorane and decapitated. Brains were quickly removed and immerse in "cutting" solution (4°C) containing the following (in mM): 87 NaCl, 25 NaHCO~3~, 2.5 KCl, 1.25 NaH~2~PO~4~, 7.5 MgSO~4~, 0.5 CaCl~2~, 2 pyruvic acid, 3 myo-inositol, 0.4 ascorbic acid, 25 glucose and 70 sucrose (equilibrated with 95% O~2~ / 5% CO~2~). Slices were cut with a vibratome (Leica) in cutting solution and then incubated in oxygenated artificial cerebrospinal fluid (ASCF) containing the following (in mM): 125 NaCl, 3 KCl, 2.5 CaCl~2~, 1.3 MgSO~4~, 1.25 mM NaH~2~PO~4~, 26 mM NaHCO~3~, 2 pyruvic acid, 3 myo-inositol, 0.4 ascorbic acid and 16 mM glucose (pH 7.4), initially at 34°C for 30 min, and subsequently at room temperature, before being transfer to the recording chamber. Recordings were obtained at 30°C. Synaptic events were recorded in whole-cell, voltage- or current-clamp mode from L2/3 and deep L5 PNs of mouse primary barrel somatosensory cortex visually identified using infrared video microscopy ([@bib42]). For voltage clamp experiments of L2/3 PNs, electrodes (with a tip resistance of 2-4 MΩ) were filled with a cesium-based internal solution (in mM): 120 CsMeSO~4~, 8 CsCl, 10 HEPES, 10 BAPTA, 4 NaCl, 2 CaCl~2~, 4 Mg-ATP, 0.3 Na-GTP, 4 phosphocreatine di(tris), 0.5 QX-314-Cl; pH adjusted to 7.2 with CsOH; 280-300 mOsm. Under these recording conditions, activation of GABA~A~ receptors resulted in outward currents at a holding potential (*V*~h~) of +10 mV. In current clamp experiments electrodes were filled with a potassium-based intracellular solution containing (in mM): 130 K-gluconate, 3 KCl, 10 HEPES, 0.2 EGTA, 2 MgCl~2,~ 4 Mg-ATP, 0.3 Na-GTP, 4 phosphocreatine di(tris); pH adjusted to 7.2 with KOH; 280-300 mOsm. E~Cl~ was approximately −77 mV based on the Nernst equation, without correction for gluconate-generated liquid junction potential. In current-clamp mode cells were recorded at their resting membrane potential unless for [Figure 4](#fig4){ref-type="fig"} where occasionally depolarization (max. current injection 30pA) was required in order to induce firing of L5 PN upon L2/3 light activation. In voltage-clamp experiments, access resistance was on average \< 20 MΩ and monitored throughout the experiment. Recordings were discarded from analysis if the resistance changed by \> 20% over the course of the experiment. In current-clamp experiments, input resistance was monitored with small current steps (−30 pA for 600 ms) and cells were excluded if it changed by \> 25%. ODQ was obtain from R&D Systems Europe.

### *In vitro* data analysis {#sec4.4.3}

Signals were amplified, using a Multiclamp 700B patch-clamp amplifier (Axon Instruments, Foster City, California, United States), sampled at 50 kHz and filtered at 4 or 10 kHz for voltage and current-clamp mode, respectively. Data were analyzed using pClamp (Axon Instruments), IGOR PRO 5.0, (Wavemetrics), MATLAB (MathWorks) and GraphPad Prism software.

#### IPSCs during rhythmic activity *in vitro* {#sec4.4.3.1}

Custom written software (Detector, courtesy J. R. Huguenard, Stanford University) was used for analyzing GABAergic events, as previously described ([@bib47], [@bib70]). Briefly, individual events were detected with a threshold-triggered process from a differentiated copy of the real trace. For each cell, the detection criteria (threshold and duration of trigger for detection) were adjusted to ignore slow membrane fluctuations and electric noise while allowing maximal discrimination of IPSCs. Detection frames were regularly inspected visually to ensure that the detector was working properly.

#### Spike probability of L5 PNs during rhythmic activity *in vitro* {#sec4.4.3.2}

Analysis of the temporal relationship between L5 PNs spikes and ongoing rhythmic activity was analyzed using custom written scripts in MATLAB (The MathWorks, Inc., Natick, Massachusetts, United States). Briefly, spikes were extracted using a threshold of −10 mV on the membrane potential trace, and the times of the action potential peaks were extracted after cubic spline data interpolation of the waveforms around the spikes. Next, IPSC peak positions from the detector (see above) were adjusted with a 50-samples smoothing of the current waveform, and delays between action potential timing and IPSC occurrence were computed. Finally, histograms between spike times and IPSC occurrences were generated using 2 ms binning and converted into probability distributions. Due to strong reduction of spike rates following LTPi induction, analysis of spike timing was restricted to PNs exhibiting a sufficient number of spikes after LTPi-inducing bursts.

### Photostimulation {#sec4.4.4}

ChR2 or ChETA activation was induced by light flashes on cortical slices, using a 20 mW LED (λ = 470 nm, Cairn research, UK) collimated and coupled to the epifluorescence path of a Zeiss AxioExaminer microscope, using a 40X water immersion (N.A. 1) lens. In order to trigger robust oscillations ([Figure 4](#fig4){ref-type="fig"}), light ramps had a duration of 1--2 s, started at zero intensity and reached a final intensity of 9 mW/mm^2^. The illuminated area was measured by bleaching coverslipped fluorescein. The light power was calculated with a power-meter at the tip of the 40X objective and was not corrected for the actual intensity at the tissue. The stimulus intensity was adapted to each slice, depending on the opsin expression and it was repeated with a frequency of 0.025 Hz. ChR2 activation was also obtained by brief square light pulses (ranging between 0.5 and 1 ms) evoking postsynaptic potentials in L5 PNs ([Figure 1](#fig1){ref-type="fig"}).

### Immunofluorescence {#sec4.4.5}

In order to check proper electroporation of both plasmids in the somatosensory cortex, in some cases ([Figure 1](#fig1){ref-type="fig"}A), slices used for electrophysiology experiments were fixed overnight in 4% paraformaldehyde in phosphate buffer (PB, pH 7.4) at 4°C. Slices were then rinsed three times at room temperature (10 min each time) in PB and were then rinsed three times in PB (10 min each) at room temperature and coverslipped in mounting medium. Immunofluorescence was then observed with an ApoTome.2 microscope (Zeiss) and images were acquired using a 10x objective.

### Preparation for *in vivo* electrophysiology {#sec4.4.6}

Two- to three-week-old naive C57BL/6 or Scnn1a x tdTomato mice of both sexes were anesthetized via intraperitoneal (i.p.) injection with 18% urethane (1.8 g/kg in sodium lactate ringer solution) and placed on a stereotaxic apparatus. The body temperature was constantly monitored and kept at 37°C with a heating blanket. Eye ointment was applied to prevent dehydration. To ensure a deep and constant level of anesthesia, vibrissae movement, eyelid reflex, response to tail, and toe pinching were visually controlled before and during the surgery. Subcutaneous injections of atropine (0.07 mg/kg) and dexamethasone (0.2 mg/kg) were used to maintain clear airways and prevent edema, respectively. A mix of local lidocaine and bupivacaine injection was performed over the cranial area of interest and, after a few minutes, a longitudinal incision was performed to expose the skull. A stainless steel head post was sealed on to the mouse skull using dental acrylic cement. A small craniotomy (≥1 mm diameter) was made on the right hemisphere to target the primary somatosensory cortex according to stereotaxic coordinates at −1 from bregma: 3 lateral). A second small craniotomy (approximately 200 μm far from the previous) was drilled for the separate entry of the LFP pipette. Dura was not removed. The exposed cortical surface was superfused with warm HEPES-buffered extracellular solution (in mM: 125 NaCl, 5 KCl, 10 glucose, 10 HEPES, 1.8 CaCl2 and 1 MgSO4 (pH 7.2)) in order to maintain ionic balance and prevent desiccation.

### *In Vivo* LFP and whole-cell patch recording {#sec4.4.7}

Patch pipettes (5-7 MΩ) of 1.5 mm external diameter borosilicate glass (WPI) were pulled on a Narishige P100 Vertica Puller and filled with (in mM): 135 K-gluconate, 6 KCl, 10 HEPES, 1 EGTA, 4 MgATP, Na2ATP and 8 phosphocreatin, pH adjusted to 7.2 with KOH, 290-295 mOsm. Pipette capacitance was neutralized before break-in. Whole-cell patch-clamp recordings of L5 pyramidal neurons were performed following the standard techniques for blind patching ([@bib48]). High positive pressure was applied to the pipette to prevent tip occlusion. After breaking the meninges, the positive pressure was immediately reduced to prevent cortical damage. Once reached L5 depth (∼700 μm depth), the pipette was then advanced in 2-μm steps, and pipette resistance was monitored in the conventional voltage clamp configuration. When the pipette resistance suddenly increased, positive pressure was relieved and small negative pressure pulses were applied to achieve GΩ seal formation. Seal resistances were always \> 1 GΩ. Recordings were made in current-clamp mode, and no holding current was applied. Typical recording durations were ∼15-20 min (which usually allowed 5 minutes recording of baseline and 5-10 minutes after burst trains). Air puff stimulation of the whisker pad was achieved by 1 s-long pulses of compressed air delivered by a picospritzer unit via a 1 mm diameter plastic tube placed at ∼20 mm from the mouse snout. To record local field potential (LFP) patch pipette (1-2MΩ) were filled with HEPES-buffered extracellular solution and inserted in the cortex. Data were acquired at 50 kHz using a Multiclamp 700B Amplifier (Molecular Devices).

### Biocytin filling {#sec4.4.8}

To certify that deep L5 PNs of S1 were being target, experiments were performed in Scnn1a x tdTomato mice allowing labeling of layer 4. Biocytin (Sigma) was added to the intracellular solution at a high concentration (0.5 g / 100 ml). Neurons were injected with large depolarizing currents in current clamp mode for fifteen times (100 ms, 1-2 nA, 1 Hz). Briefly, after *in vivo* experiments, mice were perfused and brain slices of 200 μm were performed as described above for *in vitro* electrophysiology. Slices were then fixed with 4% paraformaldehyde in phosphate buffer saline (PBS, Sigma) for at least 48 h. Following fixation, slices were incubated with the avidin-biotin complex (VECTASTAIN Elite) and a high concentration of detergent (Triton X-100, 5%) for at least two days before staining with 3,3′Diaminobenzidine (DAB, VECTASTAIN Elite) following ([@bib34]). Slices were then rinsed in PBS and coverslipped in mounting medium. Immunofluorescence and DAB staining were observed in a Micro Zeiss Routine microscope and images were acquired using a 10-40x objectives.

### Analysis of the phase modulation of AP-LFP coupling {#sec4.4.9}

The analysis of the coupling between APs and local field potentials (LFP) was done in MATLAB (Mathworks, Natick MA, US), employing custom-written scripts. Briefly, AP occurrence times were extracted from the intracellular membrane potential, recorded and sampled at 50 kHz. This was based on a supervised peak-detection algorithm, using the findpeak() MATLAB function with a −15 mV detection threshold and a 4 ms dead-time. LFP were first down sampled to 1 KHz and then digitally filtered, between 20 and 100 Hz, by a 16^th^ order bandpass Butterworth filter. The waveform was then Hilbert-transformed and its instantaneous phase component unrolled in the range \[0 ;2π \]. For each AP, the corresponding LFP phase was extracted as the corresponding value of the instantaneous phase component. A histogram was used as polar-plot or equivalently as a Cartesian-plot to display the distribution of AP-LFP phases. The confidence interval was estimated by Monte Carlo methods from 100 independent surrogate dataset, obtained upon randomly jittering each AP times (i.e., Gaussian distributed, zero mean, 10 ms standard deviation) and again estimating the distribution of the corresponding LFP phases from the Hilbert instantaneous phase waveform.

### Computational model {#sec4.4.10}

L5 PNs were modeled as integrate-and-fire neurons, with membrane potential dynamics given by$$\mathbf{c}_{\mathbf{m}}\frac{\mathbf{dV}}{\mathbf{dt}} = - \mathbf{g}_{\mathbf{m}}\mathbf{V} + \mathbf{I}\left( \mathbf{t} \right)$$where the membrane potential $\mathbf{V}$ is determined with respect to the resting potential of the cell. The membrane capacitance was $\mathbf{c}_{\mathbf{m}} =$100 pF, and the membrane conductance $\mathbf{g} =$10 nS. An action potential was emitted when the membrane potential reached a threshold value V~T~ = 40 mV. The membrane potential was subsequently reset to a value V~R~ = 0 mV. The total input to the neuron was given by$$\mathbf{I}\left( \mathbf{t} \right) = \mathbf{I}_{\mathbf{osc}}\left( \mathbf{t} \right) + \mathbf{I}_{\mathbf{bkg}}\left( \mathbf{t} \right)$$Here I~osc~ represents the oscillatory current received from L2/3 via feed-forward excitation and inhibition. It is therefore given by L2/3 oscillatory activity r~2/3~ (t),$$\mathbf{r}_{2/3}\left( \mathbf{t} \right) = \mathbf{r}_{\mathbf{f}}\mathbf{\cos}\ \left( {2\mathbf{\pi}\mathbf{f}} \right)$$filtered by a composite post-synaptic current K(t), modeled as a difference of excitatory and inhibitory components$$\mathbf{K}\left( \mathbf{t} \right) = \mathbf{K}_{\mathbf{E}}\mathbf{\exp}\left( {- \frac{\mathbf{t} - \mathbf{\delta}_{5}}{\mathbf{\tau}_{\mathbf{E}}}} \right)\mathbf{\Theta}\left( {\mathbf{t} - \mathbf{\delta}_{5}} \right) + \mathbf{K}_{\mathbf{I}}\mathbf{\exp}\left( {- \frac{\mathbf{t} - \mathbf{\delta} - \mathbf{\delta}_{5}}{\mathbf{\tau}_{\mathbf{I}}}} \right)\mathbf{\Theta}\left( {\mathbf{t} - \mathbf{\delta} - \mathbf{\delta}_{5}} \right)$$where $\mathbf{\tau}_{\mathbf{E}} =$ 2 ms and $\mathbf{\tau}_{\mathbf{I}} =$ 7 ms are the timescales of excitation and inhibition, $\mathbf{\delta}$ = 2 ms is the delay between excitatory and inhibitory inputs, $\mathbf{\delta}_{5}$ is an overall delay between L2/3 activity and L5 inputs, $\mathbf{K}_{\mathbf{E}}$ and $\mathbf{K}_{\mathbf{I}}$ represent the strengths of excitation and inhibition, and $\mathbf{\Theta}\left( \mathbf{t} \right)$ is the Heaviside step function (0 for t \< 0, 1 for t \> 1).

The oscillatory input was therefore given by the convolution between K(t) and $\mathbf{r}_{2/3}\left( \mathbf{t} \right)$**:**$$\mathbf{I}_{\mathbf{osc}}\left( \mathbf{t} \right) = {\int\limits_{- \infty}^{+ \infty}{\mathbf{d}\mathbf{\tau}\ \mathbf{r}_{2/3}}}\left( {\mathbf{t} - \mathbf{\tau}} \right)\mathbf{K}\left( \mathbf{\tau} \right).$$In the simulations, the amplitude of L2/3 oscillations was set to $\mathbf{r}_{\mathbf{f}}$ = 3 spk/s, the excitatory strength of the PSC was set to $\mathbf{K}_{\mathbf{E}} =$20 pA and $\mathbf{K}_{\mathbf{I}}$ was varied in the range from 0 to 32 pA, so that oscillations led to composite postsynaptic currents of about 100-350 pA.

Other inputs to the L5 PN were modeled as background noise$$\mathbf{I}_{\mathbf{bkg}}\left( \mathbf{t} \right) = \mathbf{\sigma}\sqrt{\mathbf{c}_{\mathbf{m}}\mathbf{g}_{\mathbf{m}}}\mathbf{\xi}\left( \mathbf{t} \right)$$where $\mathbf{\xi}\left( \mathbf{t} \right)$ is Gaussian white noise of zero mean and unit variance, I~0~ is the mean input and $\mathbf{\sigma}$ is the background noise amplitude. In the simulations in [Figure 7](#fig7){ref-type="fig"} we used $\sigma = 22\ mV$ , while in [Figure S7](#mmc1){ref-type="supplementary-material"} we show the results of simulations for half and double noise amplitude (adjusting the mean input to keep the mean firing rate identical).

The oscillatory input from L2/3 entrains the activity of the L5 PN, which therefore acquires an oscillatory temporal structure. To determine the temporal relationship between L2/3 inputs and the output from L5PNs, we computed the trial-averaged firing rate $\mathbf{r}_{5}\left( \mathbf{t} \right)$ (histograms in [Figure 7](#fig7){ref-type="fig"}C bottom). If the oscillatory inputs are relatively weak, L5 PN activity can be approximated as$$\mathbf{r}_{5}\left( \mathbf{t} \right) = \mathbf{r}_{5}^{(0)} + \mathbf{r}_{5}^{(1)}\mathbf{\cos}\ \left( {2\mathbf{\pi} - \mathbf{\varphi}} \right)$$where the first term represents constant activity, and the second term describes an oscillation around that baseline. The phase $\mathbf{\varphi}$ quantifies the timing of the output with respect to the L2/3 oscillatory input, while the amplitude $\mathbf{r}_{5}^{(1)}$ quantifies the precision of this timing -- the larger $\mathbf{r}_{5}^{(1)}$, the more L5 PN action potentials are concentrated at the specific phase of the input given by $\mathbf{\varphi}$.

Our aim was to understand how the shape of the composite post-synaptic current K(t) influences the timing and precision of the L5 PN output. Previous theoretical works (for a review, see [@bib10]) showed that for weak inputs the phase and amplitude of the output can be decomposed into a synaptic and a neuronal contribution:$$\mathbf{\varphi} = \mathbf{\varphi}_{\mathbf{syn}} + \mathbf{\varphi}_{\mathbf{n}}$$$$\mathbf{r}_{5}^{(1)} = \mathbf{A}_{\mathbf{syn}}\mathbf{A}_{\mathbf{n}.}$$Here the neuronal contributions $\mathbf{\varphi}_{\mathbf{n}}\ $and $\mathbf{A}_{\mathbf{n}}$ to the phase and amplitude depend on the background noise input and the specific neural model ([@bib12], [@bib25]). As we did not vary the background noise parameters, we treated them as arbitrary constants.

The synaptic phase lag $\mathbf{\varphi}_{\mathbf{syn}}$ and amplitude $\mathbf{A}_{\mathbf{syn}}$ are given by the phase and amplitude of the Fourier transform $\widetilde{\mathbf{K}}\left( \mathbf{f} \right)$ of the Fourier transform of the composite post-synaptic current K(t):$$\widetilde{\mathbf{K}}\left( \mathbf{f} \right) = \frac{\mathbf{e}^{- 2\mathbf{\pi}\mathbf{if}\mathbf{\delta}_{5}}}{\sqrt{2\mathbf{\pi}}}\left( {\frac{\mathbf{K}_{\mathbf{E}}\mathbf{\tau}_{\mathbf{E}}}{2\mathbf{\pi}\mathbf{if}\mathbf{\tau}_{\mathbf{E}} + 1} - \frac{\mathbf{K}_{\mathbf{I}}\mathbf{\tau}_{\mathbf{I}}}{2\mathbf{\pi}\mathbf{if}\mathbf{\tau}_{\mathbf{I}} + 1}\mathbf{e}^{- 2\mathbf{\pi}\mathbf{if}\mathbf{\delta}}} \right).$$The phase and amplitude of $\widetilde{\mathbf{K}}\left( \mathbf{f} \right)$ depend on the frequency of the oscillation (as well as on synaptic parameters). Since we were looking at the response to L2/3 inputs oscillating at 30Hz, $\widetilde{\mathbf{K}}\left( \mathbf{f} \right)$ was evaluated at f = 30 Hz, and a range of $\mathbf{K}_{\mathbf{I}}/\mathbf{K}_{\mathbf{E}}$ values was used (see [Figure 7](#fig7){ref-type="fig"}). The phase and the amplitude of the obtained values were then compared with direct fits of a cosine function to simulation results ([Figure 7](#fig7){ref-type="fig"}). Note that the above theoretical prediction for the phase and amplitude of L5 PN output rely on a linear approximation, expected to be accurate if the amplitude of the oscillations is not too strong. The comparison between the prediction and simulations indeed shows deviations from the predicted amplitude when inhibition becomes very strong ([Figure 7](#fig7){ref-type="fig"}E).

Quantification and Statistical Analysis {#sec4.5}
---------------------------------------

Analysis of LTP was performed by comparing the mean amplitude of light evoked inhibitory postsynaptic currents in the 5 last minutes of the plasticity to the baseline period ([Figure 1](#fig1){ref-type="fig"}). Unless otherwise indicated, statistical comparisons were done between raw values.

Normality of the data was assessed (D'Agostino & Pearson omnibus normality test). Normal distributions were statistically compared using paired t test two-tailed and Wilcoxon matched-pairs signed rank test was used as a non-parametric test. Group data was analyzed using the either Friedman test followed by Dunn's multiple comparison test when comparing repeated-measurements (e.g., [Figure 2](#fig2){ref-type="fig"}D for different frequencies and [Figure 3](#fig3){ref-type="fig"}D), or Kruskal-Wallis Test (e.g., [Figure 4](#fig4){ref-type="fig"}F). Differences were considered significant if p \< 0.05. Values are presented as mean ± SEM of *n* experiments.

Gain modulation ([Figure 2](#fig2){ref-type="fig"}F) was calculated from the average slope (F′) of the fits between 5% and 75% of its maximum value ([@bib63]). Only cells where the fit was possible were included in the calculation. Changes in gain (ΔGain) were computed as follows: $$\Delta Gain = \frac{\left( {\text{F}^{\text{'}}\text{LTPi} - \text{~F}^{\text{'}}\text{bsl}} \right)}{\text{F'~bsl}} \times 100.$$Additive offset shifts (ΔOffset) were defined as the difference between the half-maximum frequencies of the fits for the two conditions LTPi and baseline.

Data and Code Availability {#sec4.6}
--------------------------

The code used to generate the model in [Figures 7](#fig7){ref-type="fig"} and [S7](#mmc1){ref-type="supplementary-material"} is available at <https://github.com/MathildeBigot>. All the other scripts supporting this study have not been deposited in a public repository, because they are standard analysis routines. However, they are available from the corresponding author upon request. The source datasets is available in Mendeley Data (<https://data.mendeley.com/datasets/75t97ncb24/2>).

Supplemental Information {#app2}
========================

Document S1. Figures S1--S7Document S2. Article plus Supplemental Information

We thank Caroline Mailhes and Geeske M. van Woerden for initial help with *in utero* electroporation. We thank Michael Graupner and Christoph Schmidt-Hieber for advice regarding establishing the *in vivo* whole-cell recording technique. This work was supported by the European Research Council (ERC) under the European Community 7th Framework Program (FP7/2007-2013/ERC grant agreement 200808), "Investissements d'avenir" ANR-10-IAIHU-06, Agence Nationale de la Recherche (ANR-13-BSV4-0015-01, ANR-16-CE16-0007-02, ANR-17-CE16-0026-01, and ANR-18-CE16-0001-01), Fondation Recherche Médicale (Equipe FRM DEQ20150331684), a NARSAD independent investigator grant, and a grant from the Institut du Cerveau et de la Moelle Épinière (Paris) (to A.B.). S.O. was supported by the program Emergences of the City of Paris and the program "Investissements d'Avenir" launched by the French government and implemented by the ANR, with references ANR-10-LABX-0087 IEC and ANR- 11-IDEX-0001-02 PSL Research University.

Author Contributions {#sec5}
====================

Conceptualization, J.L. and A.B.; Methodology, J.L., A.M.D.S., A.P., C.D., A.A., M.G., M.B., and S.O.; Investigation, J.L., A.M.D.S., and M.B.; Formal Analysis, J.L., A.P., A.M.D.S., M.B., S.O., M.G., and A.B.; Writing -- Original Draft, J.L. and A.B.; Writing -- Review & Editing, J.L., S.O., and A.B.; Resources, S.O. and A.B.; Visualization, J.L., S.O., and A.B.; Supervision, J.L., S.O., and A.B.; Project Administration, J.L. and A.B.; Funding Acquisition, S.O. and A.B.

Declaration of Interests {#sec6}
========================

The authors declare no competing interests.

Supplemental Information can be found online at <https://doi.org/10.1016/j.celrep.2019.12.052>.

[^1]: Present address: Istituto Superiore di Sanità, 00161 Rome, Italy

[^2]: Lead Contact
