{
"nbformat": 4,
"nbformat_minor": 4,
"metadata": {
"title": "Precise excitation-inhibition balance controls gain and timing in the hippocampus",
"authors": [
{
"givenNames": [
"Aanchal"
],
"familyNames": [
"Bhatia"
],
"type": "Person",
"affiliations": [
{
"name": "National Centre for Biological Sciences",
"parentOrganization": {
"name": "Tata Institute of Fundamental Research",
"type": "Organization"
},
"address": {
"addressCountry": "India",
"addressLocality": "Bangalore",
"type": "PostalAddress"
},
"type": "Organization"
}
],
"name": "Aanchal Bhatia"
},
{
"givenNames": [
"Sahil"
],
"familyNames": [
"Moza"
],
"type": "Person",
"affiliations": [
{
"name": "National Centre for Biological Sciences",
"parentOrganization": {
"name": "Tata Institute of Fundamental Research",
"type": "Organization"
},
"address": {
"addressCountry": "India",
"addressLocality": "Bangalore",
"type": "PostalAddress"
},
"type": "Organization"
}
],
"name": "Sahil Moza"
},
{
"givenNames": [
"Upinder",
"Singh"
],
"familyNames": [
"Bhalla"
],
"type": "Person",
"emails": [
"bhalla@ncbs.res.in"
],
"affiliations": [
{
"name": "National Centre for Biological Sciences",
"parentOrganization": {
"name": "Tata Institute of Fundamental Research",
"type": "Organization"
},
"address": {
"addressCountry": "India",
"addressLocality": "Bangalore",
"type": "PostalAddress"
},
"type": "Organization"
}
],
"name": "Upinder Singh Bhalla"
}
],
"editors": [
{
"givenNames": [
"Ronald",
"L"
],
"familyNames": [
"Calabrese"
],
"type": "Person",
"affiliations": [
{
"name": "Emory University",
"address": {
"addressCountry": "United States",
"type": "PostalAddress"
},
"type": "Organization"
}
]
}
],
"datePublished": {
"value": "2019-04-25",
"type": "Date"
},
"dateReceived": {
"value": "2018-11-06",
"type": "Date"
},
"dateAccepted": {
"value": "2019-04-10",
"type": "Date"
},
"description": "Excitation-inhibition (EI) balance controls excitability, dynamic range, and input gating in many brain circuits. Subsets of synaptic input can be selected or 'gated' by precise modulation of finely tuned EI balance, but assessing the granularity of EI balance requires combinatorial analysis of excitatory and inhibitory inputs. Using patterned optogenetic stimulation of mouse hippocampal CA3 neurons, we show that hundreds of unique CA3 input combinations recruit excitation and inhibition with a nearly identical ratio, demonstrating precise EI balance at the hippocampus. Crucially, the delay between excitation and inhibition decreases as excitatory input increases from a few synapses to tens of synapses. This creates a dynamic millisecond-range window for postsynaptic excitation, controlling membrane depolarization amplitude and timing via subthreshold divisive normalization. We suggest that this combination of precise EI balance and dynamic EI delays forms a general mechanism for millisecond-range input gating and subthreshold gain control in feedforward networks.",
"isPartOf": {
"volumeNumber": 8,
"isPartOf": {
"title": "eLife",
"issns": [
"2050-084X"
],
"identifiers": [
{
"name": "nlm-ta",
"propertyID": "https://registry.identifiers.org/registry/nlm-ta",
"value": "elife",
"type": "PropertyValue"
},
{
"name": "publisher-id",
"propertyID": "https://registry.identifiers.org/registry/publisher-id",
"value": "eLife",
"type": "PropertyValue"
}
],
"publisher": {
"name": "eLife Sciences Publications, Ltd",
"type": "Organization"
},
"type": "Periodical"
},
"type": "PublicationVolume"
},
"licenses": [
{
"url": "http://creativecommons.org/licenses/by/4.0/",
"content": [
{
"content": [
"This article is distributed under the terms of the ",
{
"content": [
"Creative Commons Attribution License"
],
"target": "http://creativecommons.org/licenses/by/4.0/",
"type": "Link"
},
", which permits unrestricted use and redistribution provided that the original author and source are credited."
],
"type": "Paragraph"
}
],
"type": "CreativeWork"
}
],
"keywords": [
"hippocampus",
"EI balance",
"subthreshold",
"channelrhodopsin",
"neural coding",
"inhibition",
"Mouse"
],
"identifiers": [
{
"name": "publisher-id",
"propertyID": "https://registry.identifiers.org/registry/publisher-id",
"value": 43415,
"type": "PropertyValue"
},
{
"name": "doi",
"propertyID": "https://registry.identifiers.org/registry/doi",
"value": "10.7554/eLife.43415",
"type": "PropertyValue"
},
{
"name": "elocation-id",
"propertyID": "https://registry.identifiers.org/registry/elocation-id",
"value": "e43415",
"type": "PropertyValue"
}
],
"fundedBy": [
{
"identifiers": [
{
"value": "UGC/ISF No. F 6-18/2014 (IC)",
"type": "PropertyValue"
}
],
"funders": [
{
"name": "University Grants Commission",
"type": "Organization"
}
],
"type": "MonetaryGrant"
},
{
"identifiers": [
{
"value": "UGC/ISF No. F 6-18/2014 (IC)",
"type": "PropertyValue"
}
],
"funders": [
{
"name": "Israel Science Foundation",
"type": "Organization"
}
],
"type": "MonetaryGrant"
},
{
"identifiers": [
{
"value": "Senior Research Fellowship",
"type": "PropertyValue"
}
],
"funders": [
{
"name": "Council of Scientific and Industrial Research",
"type": "Organization"
}
],
"type": "MonetaryGrant"
},
{
"identifiers": [
{
"value": "Graduate Student Fellowship",
"type": "PropertyValue"
}
],
"funders": [
{
"name": "National Centre for Biological Sciences",
"type": "Organization"
}
],
"type": "MonetaryGrant"
}
],
"about": [
{
"name": "Neuroscience",
"type": "DefinedTerm"
}
],
"genre": [
"Research Article"
],
"references": [
{
"id": "bib1",
"authors": [
{
"givenNames": [
"OJ"
],
"familyNames": [
"Ahmed"
],
"type": "Person"
},
{
"givenNames": [
"MR"
],
"familyNames": [
"Mehta"
],
"type": "Person"
}
],
"title": "The hippocampal rate code: anatomy, physiology and theory",
"datePublished": {
"value": "2009",
"type": "Date"
},
"pageStart": 329,
"pageEnd": 338,
"isPartOf": {
"volumeNumber": 32,
"isPartOf": {
"name": "Trends in Neurosciences",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib2",
"authors": [
{
"givenNames": [
"JS"
],
"familyNames": [
"Anderson"
],
"type": "Person"
},
{
"givenNames": [
"M"
],
"familyNames": [
"Carandini"
],
"type": "Person"
},
{
"givenNames": [
"D"
],
"familyNames": [
"Ferster"
],
"type": "Person"
}
],
"title": "Orientation tuning of input conductance, excitation, and inhibition in cat primary visual cortex",
"datePublished": {
"value": "2000",
"type": "Date"
},
"pageStart": 909,
"pageEnd": 926,
"isPartOf": {
"volumeNumber": 84,
"isPartOf": {
"name": "Journal of Neurophysiology",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib3",
"authors": [
{
"givenNames": [
"BV"
],
"familyNames": [
"Atallah"
],
"type": "Person"
},
{
"givenNames": [
"M"
],
"familyNames": [
"Scanziani"
],
"type": "Person"
}
],
"title": "Instantaneous modulation of gamma oscillation frequency by balancing excitation with inhibition",
"datePublished": {
"value": "2009",
"type": "Date"
},
"pageStart": 566,
"pageEnd": 577,
"isPartOf": {
"volumeNumber": 62,
"isPartOf": {
"name": "Neuron",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib4",
"authors": [
{
"givenNames": [
"HC"
],
"familyNames": [
"Barron"
],
"type": "Person"
},
{
"givenNames": [
"TP"
],
"familyNames": [
"Vogels"
],
"type": "Person"
},
{
"givenNames": [
"TE"
],
"familyNames": [
"Behrens"
],
"type": "Person"
},
{
"givenNames": [
"M"
],
"familyNames": [
"Ramaswami"
],
"type": "Person"
}
],
"title": "Inhibitory engrams in perception and memory",
"datePublished": {
"value": "2017",
"type": "Date"
},
"isPartOf": {
"volumeNumber": 13,
"isPartOf": {
"name": "PNAS",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib5",
"authors": [
{
"givenNames": [
"AF"
],
"familyNames": [
"Bartley"
],
"type": "Person"
},
{
"givenNames": [
"LE"
],
"familyNames": [
"Dobrunz"
],
"type": "Person"
}
],
"title": "Short-term plasticity regulates the excitation/inhibition ratio and the temporal window for spike integration in CA1 pyramidal cells",
"datePublished": {
"value": "2015",
"type": "Date"
},
"pageStart": 1402,
"pageEnd": 1415,
"isPartOf": {
"volumeNumber": 41,
"isPartOf": {
"name": "European Journal of Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib6",
"authors": [
{
"givenNames": [
"J"
],
"familyNames": [
"Basu"
],
"type": "Person"
},
{
"givenNames": [
"JD"
],
"familyNames": [
"Zaremba"
],
"type": "Person"
},
{
"givenNames": [
"SK"
],
"familyNames": [
"Cheung"
],
"type": "Person"
},
{
"givenNames": [
"FL"
],
"familyNames": [
"Hitti"
],
"type": "Person"
},
{
"givenNames": [
"BV"
],
"familyNames": [
"Zemelman"
],
"type": "Person"
},
{
"givenNames": [
"A"
],
"familyNames": [
"Losonczy"
],
"type": "Person"
},
{
"givenNames": [
"SA"
],
"familyNames": [
"Siegelbaum"
],
"type": "Person"
}
],
"title": "Gating of hippocampal activity, plasticity, and memory by entorhinal cortex long-range inhibition",
"datePublished": {
"value": "2016",
"type": "Date"
},
"isPartOf": {
"volumeNumber": 351,
"isPartOf": {
"name": "Science",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib7",
"authors": [
{
"givenNames": [
"MJ"
],
"familyNames": [
"Bezaire"
],
"type": "Person"
},
{
"givenNames": [
"I"
],
"familyNames": [
"Soltesz"
],
"type": "Person"
}
],
"title": "Quantitative assessment of CA1 local circuits: knowledge base for interneuron-pyramidal cell connectivity",
"datePublished": {
"value": "2013",
"type": "Date"
},
"pageStart": 751,
"pageEnd": 785,
"isPartOf": {
"volumeNumber": 23,
"isPartOf": {
"name": "Hippocampus",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib8",
"authors": [
{
"givenNames": [
"KC"
],
"familyNames": [
"Bittner"
],
"type": "Person"
},
{
"givenNames": [
"C"
],
"familyNames": [
"Grienberger"
],
"type": "Person"
},
{
"givenNames": [
"SP"
],
"familyNames": [
"Vaidya"
],
"type": "Person"
},
{
"givenNames": [
"AD"
],
"familyNames": [
"Milstein"
],
"type": "Person"
},
{
"givenNames": [
"JJ"
],
"familyNames": [
"Macklin"
],
"type": "Person"
},
{
"givenNames": [
"J"
],
"familyNames": [
"Suh"
],
"type": "Person"
},
{
"givenNames": [
"S"
],
"familyNames": [
"Tonegawa"
],
"type": "Person"
},
{
"givenNames": [
"JC"
],
"familyNames": [
"Magee"
],
"type": "Person"
}
],
"title": "Conjunctive input processing drives feature selectivity in hippocampal CA1 neurons",
"datePublished": {
"value": "2015",
"type": "Date"
},
"pageStart": 1133,
"pageEnd": 1142,
"isPartOf": {
"volumeNumber": 18,
"isPartOf": {
"name": "Nature Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib9",
"authors": [
{
"givenNames": [
"KC"
],
"familyNames": [
"Bittner"
],
"type": "Person"
},
{
"givenNames": [
"AD"
],
"familyNames": [
"Milstein"
],
"type": "Person"
},
{
"givenNames": [
"C"
],
"familyNames": [
"Grienberger"
],
"type": "Person"
},
{
"givenNames": [
"S"
],
"familyNames": [
"Romani"
],
"type": "Person"
},
{
"givenNames": [
"JC"
],
"familyNames": [
"Magee"
],
"type": "Person"
}
],
"title": "Behavioral time scale synaptic plasticity underlies CA1 place fields",
"datePublished": {
"value": "2017",
"type": "Date"
},
"pageStart": 1033,
"pageEnd": 1036,
"isPartOf": {
"volumeNumber": 357,
"isPartOf": {
"name": "Science",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib10",
"authors": [
{
"givenNames": [
"RM"
],
"familyNames": [
"Bruno"
],
"type": "Person"
}
],
"title": "Synchrony in sensation",
"datePublished": {
"value": "2011",
"type": "Date"
},
"pageStart": 701,
"pageEnd": 708,
"isPartOf": {
"volumeNumber": 21,
"isPartOf": {
"name": "Current Opinion in Neurobiology",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib11",
"authors": [
{
"givenNames": [
"M"
],
"familyNames": [
"Carandini"
],
"type": "Person"
},
{
"givenNames": [
"DJ"
],
"familyNames": [
"Heeger"
],
"type": "Person"
}
],
"title": "Normalization as a canonical neural computation",
"datePublished": {
"value": "2011",
"type": "Date"
},
"pageStart": 51,
"pageEnd": 62,
"isPartOf": {
"volumeNumber": 13,
"isPartOf": {
"name": "Nature Reviews Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib12",
"authors": [
{
"givenNames": [
"S"
],
"familyNames": [
"Cash"
],
"type": "Person"
},
{
"givenNames": [
"R"
],
"familyNames": [
"Yuste"
],
"type": "Person"
}
],
"title": "Linear summation of excitatory inputs by CA1 pyramidal neurons",
"datePublished": {
"value": "1999",
"type": "Date"
},
"pageStart": 383,
"pageEnd": 394,
"isPartOf": {
"volumeNumber": 22,
"isPartOf": {
"name": "Neuron",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib13",
"authors": [
{
"givenNames": [
"JA"
],
"familyNames": [
"D'amour"
],
"type": "Person"
},
{
"givenNames": [
"RC"
],
"familyNames": [
"Froemke"
],
"type": "Person"
}
],
"title": "Inhibitory and excitatory spike-timing-dependent plasticity in the auditory cortex",
"datePublished": {
"value": "2015",
"type": "Date"
},
"pageStart": 514,
"pageEnd": 528,
"isPartOf": {
"volumeNumber": 86,
"isPartOf": {
"name": "Neuron",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib14",
"authors": [
{
"givenNames": [
"S"
],
"familyNames": [
"Denève"
],
"type": "Person"
},
{
"givenNames": [
"CK"
],
"familyNames": [
"Machens"
],
"type": "Person"
}
],
"title": "Efficient codes and balanced networks",
"datePublished": {
"value": "2016",
"type": "Date"
},
"pageStart": 375,
"pageEnd": 382,
"isPartOf": {
"volumeNumber": 19,
"isPartOf": {
"name": "Nature Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib15",
"authors": [
{
"givenNames": [
"JS"
],
"familyNames": [
"Dittman"
],
"type": "Person"
},
{
"givenNames": [
"AC"
],
"familyNames": [
"Kreitzer"
],
"type": "Person"
},
{
"givenNames": [
"WG"
],
"familyNames": [
"Regehr"
],
"type": "Person"
}
],
"title": "Interplay between facilitation, depression, and residual calcium at three presynaptic terminals",
"datePublished": {
"value": "2000",
"type": "Date"
},
"pageStart": 1374,
"pageEnd": 1385,
"isPartOf": {
"volumeNumber": 20,
"isPartOf": {
"name": "The Journal of Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib16",
"authors": [
{
"givenNames": [
"R"
],
"familyNames": [
"Enoki"
],
"type": "Person"
},
{
"givenNames": [
"M"
],
"familyNames": [
"Inoue"
],
"type": "Person"
},
{
"givenNames": [
"Y"
],
"familyNames": [
"Hashimoto"
],
"type": "Person"
},
{
"givenNames": [
"Y"
],
"familyNames": [
"Kudo"
],
"type": "Person"
},
{
"givenNames": [
"H"
],
"familyNames": [
"Miyakawa"
],
"type": "Person"
}
],
"title": "GABAergic control of synaptic summation in hippocampal CA1 pyramidal neurons",
"datePublished": {
"value": "2001",
"type": "Date"
},
"pageStart": 683,
"pageEnd": 689,
"isPartOf": {
"volumeNumber": 11,
"isPartOf": {
"name": "Hippocampus",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib17",
"authors": [
{
"givenNames": [
"RC"
],
"familyNames": [
"Froemke"
],
"type": "Person"
},
{
"givenNames": [
"MM"
],
"familyNames": [
"Merzenich"
],
"type": "Person"
},
{
"givenNames": [
"CE"
],
"familyNames": [
"Schreiner"
],
"type": "Person"
}
],
"title": "A synaptic memory trace for cortical receptive field plasticity",
"datePublished": {
"value": "2007",
"type": "Date"
},
"pageStart": 425,
"pageEnd": 429,
"isPartOf": {
"volumeNumber": 450,
"isPartOf": {
"name": "Nature",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib18",
"authors": [
{
"givenNames": [
"RC"
],
"familyNames": [
"Froemke"
],
"type": "Person"
}
],
"title": "Plasticity of cortical excitatory-inhibitory balance",
"datePublished": {
"value": "2015",
"type": "Date"
},
"pageStart": 195,
"pageEnd": 219,
"isPartOf": {
"volumeNumber": 38,
"isPartOf": {
"name": "Annual Review of Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib19",
"authors": [
{
"givenNames": [
"L"
],
"familyNames": [
"Gabernet"
],
"type": "Person"
},
{
"givenNames": [
"SP"
],
"familyNames": [
"Jadhav"
],
"type": "Person"
},
{
"givenNames": [
"DE"
],
"familyNames": [
"Feldman"
],
"type": "Person"
},
{
"givenNames": [
"M"
],
"familyNames": [
"Carandini"
],
"type": "Person"
},
{
"givenNames": [
"M"
],
"familyNames": [
"Scanziani"
],
"type": "Person"
}
],
"title": "Somatosensory integration controlled by dynamic thalamocortical feed-forward inhibition",
"datePublished": {
"value": "2005",
"type": "Date"
},
"pageStart": 315,
"pageEnd": 327,
"isPartOf": {
"volumeNumber": 48,
"isPartOf": {
"name": "Neuron",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib20",
"authors": [
{
"givenNames": [
"C"
],
"familyNames": [
"Grienberger"
],
"type": "Person"
},
{
"givenNames": [
"AD"
],
"familyNames": [
"Milstein"
],
"type": "Person"
},
{
"givenNames": [
"KC"
],
"familyNames": [
"Bittner"
],
"type": "Person"
},
{
"givenNames": [
"S"
],
"familyNames": [
"Romani"
],
"type": "Person"
},
{
"givenNames": [
"JC"
],
"familyNames": [
"Magee"
],
"type": "Person"
}
],
"title": "Inhibitory suppression of heterogeneously tuned excitation enhances spatial coding in CA1 place cells",
"datePublished": {
"value": "2017",
"type": "Date"
},
"pageStart": 417,
"pageEnd": 426,
"isPartOf": {
"volumeNumber": 20,
"isPartOf": {
"name": "Nature Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib21",
"authors": [
{
"givenNames": [
"B"
],
"familyNames": [
"Haider"
],
"type": "Person"
},
{
"givenNames": [
"A"
],
"familyNames": [
"Duque"
],
"type": "Person"
},
{
"givenNames": [
"AR"
],
"familyNames": [
"Hasenstaub"
],
"type": "Person"
},
{
"givenNames": [
"DA"
],
"familyNames": [
"McCormick"
],
"type": "Person"
}
],
"title": "Neocortical network activity in vivo is generated through a dynamic balance of excitation and inhibition",
"datePublished": {
"value": "2006",
"type": "Date"
},
"pageStart": 4535,
"pageEnd": 4545,
"isPartOf": {
"volumeNumber": 26,
"isPartOf": {
"name": "Journal of Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib22",
"authors": [
{
"givenNames": [
"JE"
],
"familyNames": [
"Heiss"
],
"type": "Person"
},
{
"givenNames": [
"Y"
],
"familyNames": [
"Katz"
],
"type": "Person"
},
{
"givenNames": [
"E"
],
"familyNames": [
"Ganmor"
],
"type": "Person"
},
{
"givenNames": [
"I"
],
"familyNames": [
"Lampl"
],
"type": "Person"
}
],
"title": "Shift in the balance between excitation and inhibition during sensory adaptation of S1 neurons",
"datePublished": {
"value": "2008",
"type": "Date"
},
"pageStart": 13320,
"pageEnd": 13330,
"isPartOf": {
"volumeNumber": 28,
"isPartOf": {
"name": "Journal of Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib23",
"authors": [
{
"givenNames": [
"G"
],
"familyNames": [
"Hennequin"
],
"type": "Person"
},
{
"givenNames": [
"EJ"
],
"familyNames": [
"Agnes"
],
"type": "Person"
},
{
"givenNames": [
"TP"
],
"familyNames": [
"Vogels"
],
"type": "Person"
}
],
"title": "Inhibitory plasticity: balance, control, and codependence",
"datePublished": {
"value": "2017",
"type": "Date"
},
"pageStart": 557,
"pageEnd": 579,
"isPartOf": {
"volumeNumber": 40,
"isPartOf": {
"name": "Annual Review of Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib24",
"authors": [
{
"givenNames": [
"MJ"
],
"familyNames": [
"Higley"
],
"type": "Person"
},
{
"givenNames": [
"D"
],
"familyNames": [
"Contreras"
],
"type": "Person"
}
],
"title": "Balanced excitation and inhibition determine spike timing during frequency adaptation",
"datePublished": {
"value": "2006",
"type": "Date"
},
"pageStart": 448,
"pageEnd": 457,
"isPartOf": {
"volumeNumber": 26,
"isPartOf": {
"name": "Journal of Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib25",
"authors": [
{
"givenNames": [
"DM"
],
"familyNames": [
"Iascone"
],
"type": "Person"
},
{
"givenNames": [
"Y"
],
"familyNames": [
"Li"
],
"type": "Person"
},
{
"givenNames": [
"U"
],
"familyNames": [
"Sumbul"
],
"type": "Person"
},
{
"givenNames": [
"M"
],
"familyNames": [
"Doron"
],
"type": "Person"
},
{
"givenNames": [
"H"
],
"familyNames": [
"Chen"
],
"type": "Person"
},
{
"givenNames": [
"V"
],
"familyNames": [
"Andreu"
],
"type": "Person"
},
{
"givenNames": [
"F"
],
"familyNames": [
"Polleux"
],
"type": "Person"
}
],
"title": "Whole-neuron synaptic mapping reveals local balance between excitatory and inhibitory synapse organization",
"datePublished": {
"value": "2018",
"type": "Date"
},
"isPartOf": {
"name": "bioRxiv",
"type": "Periodical"
},
"type": "Article"
},
{
"id": "bib26",
"authors": [
{
"givenNames": [
"N"
],
"familyNames": [
"Ishizuka"
],
"type": "Person"
},
{
"givenNames": [
"J"
],
"familyNames": [
"Weber"
],
"type": "Person"
},
{
"givenNames": [
"DG"
],
"familyNames": [
"Amaral"
],
"type": "Person"
}
],
"title": "Organization of intrahippocampal projections originating from CA3 pyramidal cells in the rat",
"datePublished": {
"value": "1990",
"type": "Date"
},
"pageStart": 580,
"pageEnd": 623,
"isPartOf": {
"volumeNumber": 295,
"isPartOf": {
"name": "The Journal of Comparative Neurology",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib27",
"authors": [
{
"givenNames": [
"N"
],
"familyNames": [
"Ishizuka"
],
"type": "Person"
},
{
"givenNames": [
"WM"
],
"familyNames": [
"Cowan"
],
"type": "Person"
},
{
"givenNames": [
"DG"
],
"familyNames": [
"Amaral"
],
"type": "Person"
}
],
"title": "A quantitative analysis of the dendritic organization of pyramidal cells in the rat hippocampus",
"datePublished": {
"value": "1995",
"type": "Date"
},
"pageStart": 17,
"pageEnd": 45,
"isPartOf": {
"volumeNumber": 362,
"isPartOf": {
"name": "The Journal of Comparative Neurology",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib28",
"authors": [
{
"givenNames": [
"O"
],
"familyNames": [
"Jensen"
],
"type": "Person"
},
{
"givenNames": [
"JE"
],
"familyNames": [
"Lisman"
],
"type": "Person"
}
],
"title": "Position reconstruction from an ensemble of hippocampal place cells: contribution of theta phase coding",
"datePublished": {
"value": "2000",
"type": "Date"
},
"pageStart": 2602,
"pageEnd": 2609,
"isPartOf": {
"volumeNumber": 83,
"isPartOf": {
"name": "Journal of Neurophysiology",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib29",
"authors": [
{
"givenNames": [
"H"
],
"familyNames": [
"Kim"
],
"type": "Person"
},
{
"givenNames": [
"S"
],
"familyNames": [
"Ährlund-Richter"
],
"type": "Person"
},
{
"givenNames": [
"X"
],
"familyNames": [
"Wang"
],
"type": "Person"
},
{
"givenNames": [
"K"
],
"familyNames": [
"Deisseroth"
],
"type": "Person"
},
{
"givenNames": [
"M"
],
"familyNames": [
"Carlén"
],
"type": "Person"
}
],
"title": "Prefrontal parvalbumin neurons in control of attention",
"datePublished": {
"value": "2016",
"type": "Date"
},
"pageStart": 208,
"pageEnd": 218,
"isPartOf": {
"volumeNumber": 164,
"isPartOf": {
"name": "Cell",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib30",
"authors": [
{
"givenNames": [
"VA"
],
"familyNames": [
"Klyachko"
],
"type": "Person"
},
{
"givenNames": [
"CF"
],
"familyNames": [
"Stevens"
],
"type": "Person"
}
],
"title": "Excitatory and feed-forward inhibitory hippocampal synapses work synergistically as an adaptive filter of natural spike trains",
"datePublished": {
"value": "2006",
"type": "Date"
},
"isPartOf": {
"volumeNumber": 4,
"isPartOf": {
"name": "PLOS Biology",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib31",
"authors": [
{
"givenNames": [
"J"
],
"familyNames": [
"Kremkow"
],
"type": "Person"
},
{
"givenNames": [
"A"
],
"familyNames": [
"Aertsen"
],
"type": "Person"
},
{
"givenNames": [
"A"
],
"familyNames": [
"Kumar"
],
"type": "Person"
}
],
"title": "Gating of signal propagation in spiking neural networks by balanced and correlated excitation and inhibition",
"datePublished": {
"value": "2010",
"type": "Date"
},
"pageStart": 15760,
"pageEnd": 15768,
"isPartOf": {
"volumeNumber": 30,
"isPartOf": {
"name": "Journal of Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib32",
"authors": [
{
"givenNames": [
"D"
],
"familyNames": [
"Lee"
],
"type": "Person"
},
{
"givenNames": [
"BJ"
],
"familyNames": [
"Lin"
],
"type": "Person"
},
{
"givenNames": [
"AK"
],
"familyNames": [
"Lee"
],
"type": "Person"
}
],
"title": "Hippocampal place fields emerge upon single-cell manipulation of excitability during behavior",
"datePublished": {
"value": "2012",
"type": "Date"
},
"pageStart": 849,
"pageEnd": 853,
"isPartOf": {
"volumeNumber": 337,
"isPartOf": {
"name": "Science",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib33",
"authors": [
{
"givenNames": [
"M"
],
"familyNames": [
"Lovett-Barron"
],
"type": "Person"
},
{
"givenNames": [
"GF"
],
"familyNames": [
"Turi"
],
"type": "Person"
},
{
"givenNames": [
"P"
],
"familyNames": [
"Kaifosh"
],
"type": "Person"
},
{
"givenNames": [
"PH"
],
"familyNames": [
"Lee"
],
"type": "Person"
},
{
"givenNames": [
"F"
],
"familyNames": [
"Bolze"
],
"type": "Person"
},
{
"givenNames": [
"XH"
],
"familyNames": [
"Sun"
],
"type": "Person"
},
{
"givenNames": [
"JF"
],
"familyNames": [
"Nicoud"
],
"type": "Person"
},
{
"givenNames": [
"BV"
],
"familyNames": [
"Zemelman"
],
"type": "Person"
},
{
"givenNames": [
"SM"
],
"familyNames": [
"Sternson"
],
"type": "Person"
},
{
"givenNames": [
"A"
],
"familyNames": [
"Losonczy"
],
"type": "Person"
}
],
"title": "Regulation of neuronal input transformations by tunable dendritic inhibition",
"datePublished": {
"value": "2012",
"type": "Date"
},
"pageStart": 423,
"pageEnd": 430,
"isPartOf": {
"volumeNumber": 15,
"isPartOf": {
"name": "Nature Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib34",
"authors": [
{
"givenNames": [
"Y"
],
"familyNames": [
"Luz"
],
"type": "Person"
},
{
"givenNames": [
"M"
],
"familyNames": [
"Shamir"
],
"type": "Person"
}
],
"title": "Balancing feed-forward excitation and inhibition via hebbian inhibitory synaptic plasticity",
"datePublished": {
"value": "2012",
"type": "Date"
},
"isPartOf": {
"volumeNumber": 8,
"isPartOf": {
"name": "PLOS Computational Biology",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib35",
"authors": [
{
"givenNames": [
"JC"
],
"familyNames": [
"Magee"
],
"type": "Person"
},
{
"givenNames": [
"EP"
],
"familyNames": [
"Cook"
],
"type": "Person"
}
],
"title": "Somatic EPSP amplitude is independent of synapse location in hippocampal pyramidal neurons",
"datePublished": {
"value": "2000",
"type": "Date"
},
"pageStart": 895,
"pageEnd": 903,
"isPartOf": {
"volumeNumber": 3,
"isPartOf": {
"name": "Nature Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib36",
"authors": [
{
"givenNames": [
"VN"
],
"familyNames": [
"Murthy"
],
"type": "Person"
},
{
"givenNames": [
"TJ"
],
"familyNames": [
"Sejnowski"
],
"type": "Person"
},
{
"givenNames": [
"CF"
],
"familyNames": [
"Stevens"
],
"type": "Person"
}
],
"title": "Heterogeneous release properties of visualized individual hippocampal synapses",
"datePublished": {
"value": "1997",
"type": "Date"
},
"pageStart": 599,
"pageEnd": 612,
"isPartOf": {
"volumeNumber": 18,
"isPartOf": {
"name": "Neuron",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib37",
"authors": [
{
"givenNames": [
"M"
],
"familyNames": [
"Okun"
],
"type": "Person"
},
{
"givenNames": [
"I"
],
"familyNames": [
"Lampl"
],
"type": "Person"
}
],
"title": "Instantaneous correlation of excitation and inhibition during ongoing and sensory-evoked activities",
"datePublished": {
"value": "2008",
"type": "Date"
},
"pageStart": 535,
"pageEnd": 537,
"isPartOf": {
"volumeNumber": 11,
"isPartOf": {
"name": "Nature Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib38",
"authors": [
{
"givenNames": [
"M"
],
"familyNames": [
"Okun"
],
"type": "Person"
},
{
"givenNames": [
"I"
],
"familyNames": [
"Lampl"
],
"type": "Person"
}
],
"title": "Balance of excitation and inhibition",
"datePublished": {
"value": "2009",
"type": "Date"
},
"isPartOf": {
"volumeNumber": 4,
"isPartOf": {
"name": "Scholarpedia",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib39",
"authors": [
{
"givenNames": [
"P"
],
"familyNames": [
"Poirazi"
],
"type": "Person"
},
{
"givenNames": [
"T"
],
"familyNames": [
"Brannon"
],
"type": "Person"
},
{
"givenNames": [
"BW"
],
"familyNames": [
"Mel"
],
"type": "Person"
}
],
"title": "Arithmetic of subthreshold synaptic summation in a model CA1 pyramidal cell",
"datePublished": {
"value": "2003",
"type": "Date"
},
"pageStart": 977,
"pageEnd": 987,
"isPartOf": {
"volumeNumber": 37,
"isPartOf": {
"name": "Neuron",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib40",
"authors": [
{
"givenNames": [
"F"
],
"familyNames": [
"Pouille"
],
"type": "Person"
},
{
"givenNames": [
"A"
],
"familyNames": [
"Marin-Burgin"
],
"type": "Person"
},
{
"givenNames": [
"H"
],
"familyNames": [
"Adesnik"
],
"type": "Person"
},
{
"givenNames": [
"BV"
],
"familyNames": [
"Atallah"
],
"type": "Person"
},
{
"givenNames": [
"M"
],
"familyNames": [
"Scanziani"
],
"type": "Person"
}
],
"title": "Input normalization by global feedforward inhibition expands cortical dynamic range",
"datePublished": {
"value": "2009",
"type": "Date"
},
"pageStart": 1577,
"pageEnd": 1585,
"isPartOf": {
"volumeNumber": 12,
"isPartOf": {
"name": "Nature Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib41",
"authors": [
{
"givenNames": [
"F"
],
"familyNames": [
"Pouille"
],
"type": "Person"
},
{
"givenNames": [
"M"
],
"familyNames": [
"Scanziani"
],
"type": "Person"
}
],
"title": "Enforcement of temporal fidelity in pyramidal cells by somatic feed-forward inhibition",
"datePublished": {
"value": "2001",
"type": "Date"
},
"pageStart": 1159,
"pageEnd": 1163,
"isPartOf": {
"volumeNumber": 293,
"isPartOf": {
"name": "Science",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib42",
"authors": [
{
"givenNames": [
"S"
],
"familyNames": [
"Ray"
],
"type": "Person"
},
{
"givenNames": [
"US"
],
"familyNames": [
"Bhalla"
],
"type": "Person"
}
],
"title": "PyMOOSE: interoperable scripting in Python for MOOSE",
"datePublished": {
"value": "2008",
"type": "Date"
},
"isPartOf": {
"volumeNumber": 2,
"isPartOf": {
"name": "Frontiers in Neuroinformatics",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib43",
"authors": [
{
"givenNames": [
"CE"
],
"familyNames": [
"Shannon"
],
"type": "Person"
}
],
"title": "A mathematical theory of communication",
"datePublished": {
"value": "1948",
"type": "Date"
},
"pageStart": 379,
"pageEnd": 423,
"isPartOf": {
"volumeNumber": 27,
"isPartOf": {
"name": "Bell System Technical Journal",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib44",
"authors": [
{
"givenNames": [
"Y"
],
"familyNames": [
"Shu"
],
"type": "Person"
},
{
"givenNames": [
"A"
],
"familyNames": [
"Hasenstaub"
],
"type": "Person"
},
{
"givenNames": [
"M"
],
"familyNames": [
"Badoual"
],
"type": "Person"
},
{
"givenNames": [
"T"
],
"familyNames": [
"Bal"
],
"type": "Person"
},
{
"givenNames": [
"DA"
],
"familyNames": [
"McCormick"
],
"type": "Person"
}
],
"title": "Barrages of synaptic activity control the gain and sensitivity of cortical neurons",
"datePublished": {
"value": "2003",
"type": "Date"
},
"pageStart": 10388,
"pageEnd": 10401,
"isPartOf": {
"volumeNumber": 23,
"isPartOf": {
"name": "The Journal of Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib45",
"authors": [
{
"givenNames": [
"MV"
],
"familyNames": [
"Tsodyks"
],
"type": "Person"
},
{
"givenNames": [
"H"
],
"familyNames": [
"Markram"
],
"type": "Person"
}
],
"title": "The neural code between neocortical pyramidal neurons depends on neurotransmitter release probability",
"datePublished": {
"value": "1997",
"type": "Date"
},
"pageStart": 719,
"pageEnd": 723,
"isPartOf": {
"volumeNumber": 94,
"isPartOf": {
"name": "PNAS",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib46",
"authors": [
{
"givenNames": [
"TP"
],
"familyNames": [
"Vogels"
],
"type": "Person"
},
{
"givenNames": [
"H"
],
"familyNames": [
"Sprekeler"
],
"type": "Person"
},
{
"givenNames": [
"F"
],
"familyNames": [
"Zenke"
],
"type": "Person"
},
{
"givenNames": [
"C"
],
"familyNames": [
"Clopath"
],
"type": "Person"
},
{
"givenNames": [
"W"
],
"familyNames": [
"Gerstner"
],
"type": "Person"
}
],
"title": "Inhibitory plasticity balances excitation and inhibition in sensory pathways and memory networks",
"datePublished": {
"value": "2011",
"type": "Date"
},
"pageStart": 1569,
"pageEnd": 1573,
"isPartOf": {
"volumeNumber": 334,
"isPartOf": {
"name": "Science",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib47",
"authors": [
{
"givenNames": [
"TP"
],
"familyNames": [
"Vogels"
],
"type": "Person"
},
{
"givenNames": [
"LF"
],
"familyNames": [
"Abbott"
],
"type": "Person"
}
],
"title": "Gating multiple signals through detailed balance of excitation and inhibition in spiking networks",
"datePublished": {
"value": "2009",
"type": "Date"
},
"pageStart": 483,
"pageEnd": 491,
"isPartOf": {
"volumeNumber": 12,
"isPartOf": {
"name": "Nature Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib48",
"authors": [
{
"givenNames": [
"M"
],
"familyNames": [
"Wehr"
],
"type": "Person"
},
{
"givenNames": [
"AM"
],
"familyNames": [
"Zador"
],
"type": "Person"
}
],
"title": "Balanced inhibition underlies tuning and sharpens spike timing in auditory cortex",
"datePublished": {
"value": "2003",
"type": "Date"
},
"pageStart": 442,
"pageEnd": 446,
"isPartOf": {
"volumeNumber": 426,
"isPartOf": {
"name": "Nature",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib49",
"authors": [
{
"givenNames": [
"WB"
],
"familyNames": [
"Wilent"
],
"type": "Person"
},
{
"givenNames": [
"D"
],
"familyNames": [
"Contreras"
],
"type": "Person"
}
],
"title": "Dynamics of excitation and inhibition underlying stimulus selectivity in rat somatosensory cortex",
"datePublished": {
"value": "2005",
"type": "Date"
},
"pageStart": 1364,
"pageEnd": 1370,
"isPartOf": {
"volumeNumber": 8,
"isPartOf": {
"name": "Nature Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib50",
"authors": [
{
"givenNames": [
"W"
],
"familyNames": [
"Woodson"
],
"type": "Person"
},
{
"givenNames": [
"L"
],
"familyNames": [
"Nitecka"
],
"type": "Person"
},
{
"givenNames": [
"Y"
],
"familyNames": [
"Ben-Ari"
],
"type": "Person"
}
],
"title": "Organization of the GABAergic system in the rat hippocampal formation: a quantitative immunocytochemical study",
"datePublished": {
"value": "1989",
"type": "Date"
},
"pageStart": 254,
"pageEnd": 271,
"isPartOf": {
"volumeNumber": 280,
"isPartOf": {
"name": "The Journal of Comparative Neurology",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib51",
"authors": [
{
"givenNames": [
"M"
],
"familyNames": [
"Xue"
],
"type": "Person"
},
{
"givenNames": [
"BV"
],
"familyNames": [
"Atallah"
],
"type": "Person"
},
{
"givenNames": [
"M"
],
"familyNames": [
"Scanziani"
],
"type": "Person"
}
],
"title": "Equalizing excitation-inhibition ratios across visual cortical neurons",
"datePublished": {
"value": "2014",
"type": "Date"
},
"pageStart": 596,
"pageEnd": 600,
"isPartOf": {
"volumeNumber": 511,
"isPartOf": {
"name": "Nature",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib52",
"authors": [
{
"givenNames": [
"O"
],
"familyNames": [
"Yizhar"
],
"type": "Person"
},
{
"givenNames": [
"LE"
],
"familyNames": [
"Fenno"
],
"type": "Person"
},
{
"givenNames": [
"M"
],
"familyNames": [
"Prigge"
],
"type": "Person"
},
{
"givenNames": [
"F"
],
"familyNames": [
"Schneider"
],
"type": "Person"
},
{
"givenNames": [
"TJ"
],
"familyNames": [
"Davidson"
],
"type": "Person"
},
{
"givenNames": [
"DJ"
],
"familyNames": [
"O'Shea"
],
"type": "Person"
},
{
"givenNames": [
"VS"
],
"familyNames": [
"Sohal"
],
"type": "Person"
},
{
"givenNames": [
"I"
],
"familyNames": [
"Goshen"
],
"type": "Person"
},
{
"givenNames": [
"J"
],
"familyNames": [
"Finkelstein"
],
"type": "Person"
},
{
"givenNames": [
"JT"
],
"familyNames": [
"Paz"
],
"type": "Person"
},
{
"givenNames": [
"K"
],
"familyNames": [
"Stehfest"
],
"type": "Person"
},
{
"givenNames": [
"R"
],
"familyNames": [
"Fudim"
],
"type": "Person"
},
{
"givenNames": [
"C"
],
"familyNames": [
"Ramakrishnan"
],
"type": "Person"
},
{
"givenNames": [
"JR"
],
"familyNames": [
"Huguenard"
],
"type": "Person"
},
{
"givenNames": [
"P"
],
"familyNames": [
"Hegemann"
],
"type": "Person"
},
{
"givenNames": [
"K"
],
"familyNames": [
"Deisseroth"
],
"type": "Person"
}
],
"title": "Neocortical excitation/inhibition balance in information processing and social dysfunction",
"datePublished": {
"value": "2011",
"type": "Date"
},
"pageStart": 171,
"pageEnd": 178,
"isPartOf": {
"volumeNumber": 477,
"isPartOf": {
"name": "Nature",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib53",
"authors": [
{
"givenNames": [
"LI"
],
"familyNames": [
"Zhang"
],
"type": "Person"
},
{
"givenNames": [
"AY"
],
"familyNames": [
"Tan"
],
"type": "Person"
},
{
"givenNames": [
"CE"
],
"familyNames": [
"Schreiner"
],
"type": "Person"
},
{
"givenNames": [
"MM"
],
"familyNames": [
"Merzenich"
],
"type": "Person"
}
],
"title": "Topography and synaptic shaping of direction selectivity in primary auditory cortex",
"datePublished": {
"value": "2003",
"type": "Date"
},
"pageStart": 201,
"pageEnd": 205,
"isPartOf": {
"volumeNumber": 424,
"isPartOf": {
"name": "Nature",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
},
{
"id": "bib54",
"authors": [
{
"givenNames": [
"M"
],
"familyNames": [
"Zhou"
],
"type": "Person"
},
{
"givenNames": [
"F"
],
"familyNames": [
"Liang"
],
"type": "Person"
},
{
"givenNames": [
"XR"
],
"familyNames": [
"Xiong"
],
"type": "Person"
},
{
"givenNames": [
"L"
],
"familyNames": [
"Li"
],
"type": "Person"
},
{
"givenNames": [
"H"
],
"familyNames": [
"Li"
],
"type": "Person"
},
{
"givenNames": [
"Z"
],
"familyNames": [
"Xiao"
],
"type": "Person"
},
{
"givenNames": [
"HW"
],
"familyNames": [
"Tao"
],
"type": "Person"
},
{
"givenNames": [
"LI"
],
"familyNames": [
"Zhang"
],
"type": "Person"
}
],
"title": "Scaling down of balanced excitation and inhibition by active behavioral states in auditory cortex",
"datePublished": {
"value": "2014",
"type": "Date"
},
"pageStart": 841,
"pageEnd": 850,
"isPartOf": {
"volumeNumber": 17,
"isPartOf": {
"name": "Nature Neuroscience",
"type": "Periodical"
},
"type": "PublicationVolume"
},
"type": "Article"
}
]
},
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Introduction\n",
"\n",
"Individual neurons in the brain can receive tens of thousands of excitatory (E) and inhibitory (I) synaptic inputs. Under normal conditions, the ratio of excitatory to inhibitory input remains invariant, a robust property of the nervous system, termed EI balance (@bib2; @bib3; @bib37; @bib38; @bib48). Disruption of balance is linked with several pathologies, including epilepsy, autism spectrum disorders and schizophrenia (@bib52).\n",
"\n",
"Theoretically, neurons in ‘detailed balanced’ EI networks receive balanced responses from all subsets of presynaptic inputs (@bib47), and neurons in ‘tightly balanced’ EI networks receive inputs balanced at fast (<10 ms) timescales (@bib14). Together, these properties constitute a ‘precisely balanced’ network (@bib23). This precise balance on all synaptic subsets can be exploited by the brain for ‘input gating’. In this process, neurons can be driven by selective shifts in EI ratios at specific inputs, while other inputs remain balanced in the background. This constitutes a flexible and instantaneous information channel local to the shifted synapses (@bib31; @bib47).\n",
"\n",
"Our current understanding of EI balance is based on measurements made at single neurons in response to various stimuli. Strong EI correlations have been seen in response to series of tones in auditory cortex (@bib48; @bib53; @bib54), whisker stimulation in somatosensory cortex (@bib49), during cortical up states in vitro (@bib44) and in vivo (@bib21), during gamma oscillations in vitro and in vivo (@bib3), and during spontaneous activity (@bib37). At the synaptic scale, the ratio of excitatory and inhibitory synapses on various dendrites of a neuron has been shown to be conserved (@bib25). However, the precision and presynaptic origin of balance is not well understood. It remains to be established if EI balance arises transiently from complex temporal dynamics of several presynaptic layers, if it requires summation of inputs from multiple presynaptic populations, or if it exists even at subsets of a single presynaptic population. This granularity of EI balance, of both presynaptic identity and number of inputs, can determine the precision with which synaptic inputs can be selected or ‘independently gated’ to affect postsynaptic activity.\n",
"\n",
"In this study, we address two key open questions in the field. First, can EI balance arise even in a single layer feedforward network, and if so, at what granularity of network subsets do postsynaptic cells experience balanced excitation and inhibition? Second, how do excitation and inhibition integrate to encode and communicate information at the postsynaptic neuron? We addressed these questions in vitro, to isolate the hippocampal network from background activity, and to deliver precisely controlled combinatorial stimuli. We stimulated channelrhodopsin-2 (ChR2) expressing CA3 neurons in several combinations using optical patterns, and measured responses in CA1.\n",
"\n",
"We report that hundreds of randomly chosen subsets of CA3 neurons provide excitatory and feedforward inhibitory inputs to CA1 cells with a close to identical ratio, demonstrating for the first time, precise balance (@bib23) in the brain. On examining the integration of excitation and feedforward inhibition, we found that inhibition arrives with a dynamically varying onset delay that decreases with increasing input amplitude. This leads to a characteristic initial linear portion in the neuronal input-output curve where the inhibition arrives too late to affect peak depolarization, and a progressively diminishing output as the EI delay decreases with increasing input. This novel gain control operation, termed Subthreshold Divisive Normalization (SDN) encodes input information in both amplitude and timing of the CA1 response.\n",
"\n",
"# Results\n",
"\n",
"In our study, we first utilize and characterize an optical stimulation protocol for CA3 pyramidal neurons, and measure intracellular responses at CA1 pyramidal neurons ([Figure 1](#fig1)). We then demonstrate precise EI balance for various combinations of CA3 inputs at CA1 using voltage clamp to separate the E and I components ([Figure 2](#fig2)). Next, we measure the depolarization at CA1 due to summation of E and I using different input combinations ([Figure 3](#fig3)), and show sublinearity of summation. Expansion of the range of inputs revealed divisive normalization and suggested that another factor such as inhibitory kinetics should be included to account for the sublinearity ([Figure 4](#fig4)). In [Figure 5](#fig5), we confirm that blocking inhibition leads to much reduced sublinearity of summation, and that inhibition scales linearly with stimulus amplitude. We then establish that inhibitory delay is crucial for explaining the sublinearity in SDN ([Figure 6](#fig6)). In [Figure 7](#fig7), we show that post-synaptic potential peak amplitude and timing both carry information about the summed stimulus amplitude, and show that this information carries over to spike timing. In [Figure 8](#fig8), we summarize the analysis and suggest how SDN could contribute to input gating in the hippocampus.\n",
"\n",
"figure: Figure 1.\n",
":::\n",
"\n",
"\n",
"## Stimulating CA3-CA1 network with hundreds of optical patterns.\n",
"\n",
"(**a**) Top, schematic of the CA3-CA1 circuit with direct excitation and feedforward inhibition. Bottom, image of a hippocampus slice expressing ChR2-tdTomato (red) in CA3 in a Cre-dependent manner. Optical stimulation grid (not drawn to scale) was centered at the CA3 cell body layer and CA1 neurons were patched. (**b**) Spike response map of CA3 neuron responses with one grid square active at a time. A CA3 neuron was patched and optically stimulated, in random spatio-temporal order, on the grid locations marked with grey border. This cell spiked (marked with number inside representing spike counts over four trials) in 5 out of 24 such one square stimuli delivered. (**c**) Heatmap of CA1 responses while CA3 neurons were stimulated with one square optical stimuli. Colormap represents peak V~m~ change averaged over three repeats. (**d**) Schematic of optical stimulus patterns. Examples of combinations of N-square stimuli where N could be 1, 2, 3, 5, 7 or 9 (in rows). (**e**) Spikes in response to four repeats for the circled square, in **b**. Spike times are marked with a black tick, showing variability in evoked peak times. Blue trace at the bottom represents photodiode measurement of the stimulus duration. Scale bar for time, same as **h**. (**f**) Distribution of jitter in spike timing (SD) for all stimuli for all CA3 cells (n = 8 cells). (**g**) CA3 spiking probability (fraction of times a neuron spiked across 24 stimuli, repeated thrice) is consistent over a single recording session. Randomization of the stimulus pattern prevented ChR2 desensitization. Circles, colored as in **d** depict spiking probability on each repeat of a stimulus set with connecting lines tracking three repeats of the set (n = 7 cells). (**h**) PSPs in response to three repeats of the circled square in **c**. Peak times are marked with an asterisk. Blue traces at the bottom represent corresponding photodiode traces for the stimulus duration. (**i**) Distribution of peak PSP amplitude variability (variance/mean) for all 1-square responses (n = 28 cells, stimuli = 695). (**j**) Histogram of peak amplitudes of all PSPs elicited by all 1-square stimuli, over all CA1 cells. Grey dotted line marks the mode (n = 38 cells, trials = 8845).\n",
":::\n",
"{#fig1}\n",
"\n",
"figure: Figure 1—figure supplement 1.\n",
":::\n",
"\n",
"\n",
"## Experiment design.\n",
"\n",
"(**a**) Schematic of the light path of patterned optical stimulation. Projector’s output was minified using a lens and introduced into the light path of the microscope by reflecting off of a dichroic. We could thus focus arbitrary minified patterns at the sample plane, and switch them at short intervals. (**b**) Heat map of CA3 neuron responses with one grid square active at a time, from the neuron in [Figure 1b](#fig1). Color in grid squares represents peak V~m~ change from baseline, averaged over trials when the neuron did not spike. Locations where the cell spiked all four times are marked with a cross. (**c–f**) Four different kinds of grids used for photostimulation. The grid was made sparser to avoid stimulation of the same region with light from nearby photostimulation squares. The four different grids were (from left to right): 24 square circular sparse grid ([Figure 1](#fig1)) (13 cells), 37 point sparse circular grid (6 cells), 13 × 13 dense grid (9 cells), 15 × 15 dense grid (10 cells). For one cell we used a fifth kind of grid: 112 point dense circular grid (not shown in figure).\n",
":::\n",
"{#fig1s1}\n",
"\n",
"## Optical stimuli at CA3 elicit subthreshold responses at CA1\n",
"\n",
"To provide a wide range of non-overlapping stimuli, we projected patterned optical stimuli onto channelrhodopsin-2 (ChR2) expressing CA3 neurons in acute hippocampal slices. We used CA3-cre mice to achieve CA3-specific localization of ChR2 upon injection of a Lox-ChR2 virus ([Figure 1a](#fig1), Materials and methods). We used a Digital Micromirror Device (DMD) projector (Materials and methods, [Figure 1—figure supplement 1](#fig1s1)) to generate spatiotemporal optical patterns in the form of a grid of several 16 um x 16 um squares, each square approximating the size of a CA3 soma (@bib27) ([Figure 1d](#fig1)). This grid was centered at the CA3 cell body layer, and extended to the dendritic layer ([Figure 1a](#fig1)). Each optical pattern consisted of 1 to 9 such randomly chosen grid squares, presented to CA3 cells as stimulus, at an inter-stimulus interval of 3 s ([Figure 1a,d](#fig1), Materials and methods). In a typical experiment, several randomly chosen stimulus patterns with different number of input squares were delivered to CA3, in three successive repeats. We first characterized how CA3 responded to the grid stimulation ([Figure 1b,e,f,g](#fig1)). CA3 neurons fired reliably with a < 2 ms jitter, calculated as the standard deviation of the time of first spike ([Figure 1f](#fig1)) (n = 8 CA3 cells, inputs = 52, median = 0.44 ms, N = 1 to 9 squares). No desensitization occurred during the timeframe of an experiment, and the probability of spiking remained constant between the three repeats ([Figure 1g](#fig1)) (n = 7 CA3 cells, N = 1 to 9 squares). Thus, we could stimulate CA3 with hundreds of distinct optical stimuli in each experiment.\n",
"\n",
"We then recorded postsynaptic potentials (PSPs) evoked at patched CA1 neurons while optically stimulating CA3 cells ([Figure 1c,h,i,j](#fig1)). A wide range of stimulus positions in CA3 excited CA1 neurons ([Figure 1c](#fig1)). Stimulation of CA3 elicited excitation and feedforward inhibition at CA1 ([Figure 1a](#fig1), [Figure 2](#fig2)). Most stimuli elicited subthreshold responses (N = 1 to 9 squares). Action potentials occurred in only 0.98% of trials (183 out of 18,668 trials, n = 38 cells, N = 1 to 9 squares). This helped rule out any significant feedback inhibition from CA1 interneurons for all our experiments. Restriction of ChR2 to CA3 pyramidal cells, coupled with the fact that ~99% of all recorded CA1 responses were subthreshold, ensured that the recorded inhibition was largely feedforward (disynaptic) ([Figure 1a](#fig1)). Responses to the same 1-square stimulus were consistent, 84.74% responses showed less than 0.5 variance by mean (695 stimuli, three repeats each, n = 28 cells, N = 1 square) ([Figure 1i](#fig1)). Notably, the distribution of all one square responses had a mode at 0.25 mV, which is close to previous reports of a 0.2 mV somatic response of single synapses in CA1 neurons (@bib35) (8845 trials, n = 38 cells, N = 1 square) ([Figure 1j](#fig1)).\n",
"\n",
"## Arbitrarily chosen CA3 inputs show precise EI balance at CA1\n",
"\n",
"To examine the relationship between excitation and inhibition, we voltage clamped CA1 neurons, first at the inhibitory (−70 mV) and then at the excitatory (0 mV) reversal potential to record Excitatory and Inhibitory Post Synaptic Currents (EPSCs and IPSCs) respectively. We first presented five different patterns of 5 squares each, at both of these potentials, and recorded EPSCs and IPSCs. We found strong proportionality between excitation and inhibition for every stimulus pattern ([Figures 1d](#fig1) and [2a](#fig2)). This suggested that inputs from even random groups of CA3 neurons may be balanced at CA1. Repeats with the same stimulus pattern gave consistent responses, but different patterns evoked different responses ([Figure 2a](#fig2), [Figure 2—figure supplement 1b](#fig2s1)). This indicated that the optically driven stimuli were able to reliably activate different subsets of synaptic inputs on the target neuron. Next, we asked, in what range of input strengths does random input yield balance? We presented five different patterns for each of 1, 2, 3, 5, 7 or 9 square combinations at both recording potentials. Surprisingly, all stimuli to a cell elicited excitatory and inhibitory responses in the same ratio, irrespective of response amplitude ([Figure 2b,c](#fig2)) (n = 13 CA1 cells, area under curve, mean R^2^ = 0.89 + /- 0.06 SD, [Figure 2—figure supplement 2](#fig2s2)). Notably, the mode of single-square responses was ~0.25 mV, close to single synapse PSP estimates (@bib35) ([Figure 1j](#fig1)). However, accounting for the low (~0.2) release probabilities (P~r~) at the CA3-CA1 synapse (@bib36), we should be able to see a single synapse response if approximately 1/P~r~ synapses were activated. Hence, we estimate that the granularity of the balance as resolved by our method is of the order of 5–10 synapses ([Figure 2—figure supplement 1d,e](#fig2s1)). The slope of the regression line through all stimulus-averaged responses for a CA1 cell was used to calculate the Inhibition/Excitation (I/E) ratio for the cell. IPSC/EPSC ratio was typically between 2 and 5 ([Figure 2d](#fig2)). The variability of I/E ratios over all stimuli for a cell was lower than the variability of all stimuli across cells (for 12 out of 13 cells, [Figure 2—figure supplement 1c](#fig2s1)). The high R^2^ values for all cells showed tight proportionality for all stimuli ([Figure 2e](#fig2)). The residual distribution remained symmetric for increasing numbers of spots, again showing that they were not affected by the number of stimulus squares presented ([Figure 2—figure supplement 1a](#fig2s1)). While feedforward inhibition is expected to increase with excitation, convergence of I/E ratios for randomly chosen inputs to a cell to a single number was unexpected, since shared interneurons consist of only about 10% of the total neuronal population (@bib50; @bib7).\n",
"\n",
"figure: Figure 2.\n",
":::\n",
"\n",
"\n",
"### Excitation and inhibition are tightly balanced for all stimuli to a CA1 cell.\n",
"\n",
"(**a**) Monosynaptic excitatory postsynaptic currents (EPSCs, at −70 mV) and disynaptic inhibitory postsynaptic currents (IPSCs, at 0 mV) in response to five different stimulus combinations of 5 squares each. All combinations show proportional excitatory and inhibitory currents over six repeats. Top, schematic of 5-square stimuli. (**b**) EPSCs and IPSCs are elicited with the same I/E ratio in response to six repeats of a combination, and across six different stimuli from 1 square to nine squares, for the same cell as in **a**. Top, schematic of the stimuli. (**c**) Area under the curve for EPSC and IPSC responses, obtained by averaging over six repeats, plotted against each other for all stimuli to the cell in **a, b**. Error bars are s.d. (**d**) Summary of I/E ratios for all cells (n = 13 cells). (**e**) Summary for all cells of R^2^ values of linear regression fits through all points. Note that 11 out of 13 cells had R^2^ greater than 0.9, implying strong proportionality. (**f**) Same as **e**, but with linear regression fits for 1 and 2 square responses, showing that even small number of synapses are balanced for excitation and inhibition (n = 9 cells). (**g**) Phase plot from the model showing how tuning of synapses (_ρ_) affects observation of EI balance (R^2^) for various values of variance/mean of the basal weight distribution. Changing the scale of the basal synaptic weight distributions against tuning parameter _ρ_ affects goodness of EI balance fits. Arrow indicates where our observed synaptic weight distribution lay. (**h**) Example of EI correlations (from data) for 1 and 2 square inputs for an example cell. Bottom, schematic of the stimuli. Excitation and inhibition are colored olive and purple, respectively. Error bars are s.d. (**i**) Examples of EI correlation (from model) for small number of synapses, from the row marked with arrow in **g**. The left and right curves show low and high correlations in mean amplitude when EI synapses are untuned (_ρ_ = 0) and tuned respectively (_ρ_ = 1) (A.U. = Arbitrary Units). Colors, same as **h**. Error bars are s.d.\n",
":::\n",
"{#fig2}\n",
"\n",
"figure: Figure 2—figure supplement 1.\n",
":::\n",
"\n",
"\n",
"### Detailed balance in CA3-CA1 feedforward network.\n",
"\n",
"(**a**) Plot of residuals for all inputs (colored by N-square as shown in [Figure 1d](#fig1)) of all cells, normalized by their standard deviation. The residuals were calculated by subtracting the actual values of inhibition from the values predicted by the regression line across the Excitation-Inhibition plot. These were then standardised by dividing individual residual values by the standard deviation across all stimuli of the same number of squares. We found that standardized residuals for different input squares were distributed symmetrically across 0, showing that the I/E ratios did not change substantially across the cell. Variability within stimulus repeats compared with that across all stimuli grouped by the number of squares, for excitation (left) and inhibition (right). For both, within stimulus variability is lower. (**b**) Variability over repeats of a stimulus (within) divided by variability across all stimuli of the same number of squares delivered to a cell (across), for excitation (left) and inhibition (right). Dotted line represents equal within and across variability. For both excitation and inhibition, within-stimulus variability is lower (n = 13 cells). (**c**) Standard deviation of I/E ratios across stimuli for a given cell was lower than that across all cells (dotted line) for 12 out of 13 cells. (**d**) Two example cells where EI balance can be observed for 1 and 2 square data from voltage clamped cells (Summary plot in [Figure 2f](#fig2)). (**e**) Two example cells where EI balance can be observed for one square data from current clamped cells, as shown in [Figure 5d](#fig5). (**f**) Summary of R^2^ for one square data from all current clamped cells (n = 9). (**g**) Histogram of correlation between standard deviations (s.d.) for excitation and inhibition for all stimuli over individual cells. Mean correlation = 0.61 (n = 13 cells). (**h**) A schematic for the detailed balance model. EI correlations increase with increase in _ρ_, as well as with decrease in variance of the distribution of basal excitatory synaptic weights. (**i**) Phase plot, similar to [Figure 2g](#fig2), but with correlations of standard deviations for excitation and inhibition. Arrow marks the row with our estimated synaptic weight distribution width. (**j**) Example cell with 0.72 correlation between s.d. of excitation and inhibition repeats over individual stimuli. Excitation and inhibition are colored olive and purple respectively. (**k**) Examples for s.d. correlation from the row marked with arrow, for untuned (left, _ρ_ = 0, corr = 0.22) and tuned (right, _ρ_ = 1, corr = 0.87) synapses (A.U. = Arbitrary Units). Colors, same as **i**. Error bars are s.d.\n",
":::\n",
"{#fig2s1}\n",
"\n",
"figure: Figure 2—figure supplement 2.\n",
":::\n",
"\n",
"\n",
"### Raw data from all cells showing precise balance between excitation and inhibition.\n",
"\n",
"Individual plots for area under the curve for excitation and inhibition recorded from all cells (except the display cell in [Figure 2c](#fig2)). These were measured by clamping the cells at inhibitory (−70 mV) and excitatory (0 mV) reversal potentials respectively. Cells exhibit close proportionality between excitation and inhibition. The I/E ratio (slope) and the R^2^ (goodness of fit) for each cell are mentioned in their individual plots. Error bars are s.d.\n",
":::\n",
"{#fig2s2}\n",
"\n",
"## Detailed balance requires co-tuning of EI weights\n",
"\n",
"We next tested the hypothesis that the observed correlation between excitatory and inhibitory inputs was due to an averaged sum over many untuned (globally balanced) synapses, as opposed to a much finer granularity of tuning between excitatory and inhibitory synaptic weights (detailed balance). To address this, we modelled excitatory and inhibitory synaptic weights to a neuron with different amounts of weight tuning, parameterized by rho (_ρ_), which takes values between 0 (no tuning or global balance) and 1 (detailed balance) (Materials and methods, [Figure 2—figure supplement 1h](#fig2s1)). For values between 0 and 1, _ρ_ determined the degree of correlation between the basal excitatory and inhibitory synaptic weights. To test if weight tuning was necessary to observe balance, we modeled the summation of synaptic inputs with the premise that excitatory and inhibitory afferents will be activated strictly proportionally (number balance). We then tested how mean and variance correlations between EI amplitudes changed with different degrees of weight tuning.\n",
"\n",
"We observed tight correlations between EI inputs without weight tuning, but only if the basal synaptic weight distribution was narrow. Further, for a narrow weight distribution, the change from global to detailed balance had little effect on mean EI amplitude correlations. In contrast, weight tuning was required to see EI balance for wider synaptic weight distributions, especially for stimuli which activated small numbers of synapses ([Figure 2g](#fig2)). We next calculated the width of the smallest responses (1-square GABAzine EPSP) as a proxy for the basal weight distribution ([Figure 3—figure supplement 1a](#fig3s1)). The observed responses were broadly distributed. With this basal weight distribution, the model exhibited EI balance only when the excitatory and inhibitory synaptic weights were co-tuned, that is, maintained at the same ratio (marked with arrow in [Figure 2g](#fig2), Materials and methods, [Figure 2—figure supplement 1h](#fig2s1)).\n",
"\n",
"With the reasoning developed above, we checked for EI balance in the smallest inputs in our datasets - 1 and 2 square data from voltage clamped cells (having five or more input patterns per cell) ([Figure 2—figure supplement 1d](#fig2s1)), and only one square from current clamped cells (24 inputs per cell) ([Figure 2—figure supplement 1e,f](#fig2s1)). We found that the responses corresponding to a few synapses per input were balanced ([Figures 1j](#fig1) and [2f,h](#fig2), [Figure 2—figure supplement 1d,e,f](#fig2s1)), suggesting tuning of excitatory and inhibitory weights.\n",
"\n",
"In addition, the model also predicted a tuning dependent change in the correlations of variability of excitation and inhibition amplitudes for repeats of the same stimulus. For a wide synaptic weight distribution, increase in tuning increased EI variability correlations ([Figure 2—figure supplement 1i,k](#fig2s1)). As with EI mean correlations ([Figure 2h,i](#fig2)), weight tuning had little effect in the case of narrow synaptic weight distributions. Again, our calculated synaptic weight distribution was in the range where strong variability correlations would be seen only if synaptic weights were tuned. We found strong correlations between excitatory and inhibitory standard deviation between six repeats of the same stimulus in our voltage-clamp dataset, suggesting that there is detailed balance in the network ([Figure 2—figure supplement 1g,j](#fig2s1)).\n",
"\n",
"Thus, we present three observations using small (one and two square) stimulus strengths: a wide basal weight distribution, correlated mean EI amplitude and correlated EI amplitude variability. Together, these are inconsistent with the hypothesis that EI balance can emerge with no other requirement than a proportional increase in number of EI afferents in a globally balanced network. This supports the existence of weight tuning and hence detailed balance in the CA3-CA1 network.\n",
"\n",
"Overall, we found stimulus-invariant proportionality of excitation and inhibition for any randomly selected input, over a large range of stimulus strengths from a single presynaptic network. In addition to detailed balance, we show below that there is tight balance, that is the timing of the balanced feedforward inhibition was within a few milliseconds of the excitation ([Figure 6g,h](#fig6)). Thus, we concluded that the CA3-CA1 circuit exhibits precise (both detailed and tight) balance (@bib23).\n",
"\n",
"## Combinatorial CA3 inputs sum sublinearly at CA1\n",
"\n",
"We next asked how CA3 inputs, that lead to balanced excitatory and feedforward inhibitory conductances, transform into membrane potential change at CA1 neurons. Based on anatomical studies, CA3 projections are likely to arrive in a distributed manner over a wide region of the dendritic tree of CA1 pyramidal neuron (@bib26) ([Figure 3a](#fig3)). While pairwise summation at CA1 has been shown to be largely linear in absence of inhibition (@bib12), the degree of heterogeneity of summation in response to spatially distributed excitatory and inhibitory synaptic inputs is not well understood (except, see @bib33) . To avoid biases that may arise from a single response measure during input integration (@bib39), we examined PSPs using four different measures ([Figure 3c](#fig3)). These were peak amplitude, area under curve (AUC), average membrane potential and area under curve till peak ([Figure 3c](#fig3)).\n",
"\n",
"figure: Figure 3.\n",
":::\n",
"\n",
"\n",
"### Excitatory and feed-forward inhibitory inputs from CA3 integrate sublinearly at CA1.\n",
"\n",
"(**a**) Schematic of a neuron receiving synaptic input distributed over its dendritic tree. (**b**) Schematic of input integration. Top, five 1-square stimuli presented individually, and a single 5-square stimulus comprising of the same squares. Bottom, PSPs elicited as a response to these stimuli. 5-square PSP can be larger (supralinear, orange), equal to (linear, black), or smaller (sublinear, grey) than the sum of the single square PSPs. (**c**) A PSP trace marked with the four measures used for further calculations. PSP peak, PSP area, area to peak and mean voltage are indicated. (**d**) Schematic of the input integration plot. Each circle represents response to one stimulus combination. ‘Observed’ (true response of 5 square stimulation) on Y-axis and ‘Expected’ (linear sum of 1 square responses) is on X-axis. (**e**) Most responses for a given cell show sublinear summation for a 5-square stimulus. The four panels show sublinear responses for four different measures (mentioned in **c**) for the same cell. The grey dotted line is the regression line and the slope of the line is the scaling factor for the responses for that cell. For peak (mV), area (mV.ms), average (mV), and area to peak (mV.ms); slope = 0.27, 0.23, 0.23, 0.18; R^2^ 0.57, 0.46, 0.46, 0.26, respectively. The responses to AUC and average are similar because of the similarity in the nature of the measure. (**f**) Distribution of Observed/Expected ratio of peaks of all responses for all 5-square stimuli (mean = 0.57, s.d. = 0.31), from all recorded cells pooled. 93.35% responses to 5-square stimuli were sublinear (2513 PSPs, n = 33 cells). (**g**) Distribution of slopes for peak amplitude of 5-square stimuli (mean = 0.38, s.d. = 0.22). Regression lines for all cells show that all cells display sublinear (<1 slope) summation (n = 33 cells).\n",
":::\n",
"{#fig3}\n",
"\n",
"figure: Figure 3—figure supplement 1.\n",
":::\n",
"\n",
"\n",
"### Summation at CA3-CA1 network is sublinear.\n",
"\n",
"(**a**) Responses to 1-square photostimulation at CA3 were similar in both Control and GABAzine conditions, except close to the tail of the distribution (n = 11, Control: 1092 trials, GABAzine: 1173 trials). This demonstrates that there is very little inhibitory effect on peak Vm with one square photostimulation. (**b**) Responses to 9-squares photostimulation lead to much larger responses in the presence of GABAzine than in its absence (n = 3, Control: 142 trials, GABAzine: 144 trials). Compare this with **a**. (**c**) Slope values for three other measures (area under curve, mean voltage and area to peak) of the observed PSP on five square stimulation in all cells. Sublinearity is seen in all four measures (n = 33 cells). The slope for first measure, peak Vm, is displayed in [Figure 3g](#fig3). Area under curve and Mean voltage panels look alike due to the similarity in the nature of the measure. (**d**) Summation plots analogous to those in [Figure 3f](#fig3) for the remaining 2, 3, 7 and 9 stimuli. Observed/Expected (O/E) ratio for most stimuli was less than 1, showing sublinear summation (grey).\n",
":::\n",
"{#fig3s1}\n",
"\n",
"We looked at input integration by presenting stimulus sets of 5 input squares to a given cell, with each stimulus set ranging from 24 to 225 combinations of inputs. We initially tested the center of our range of 1–9 squares (5-square inputs) before expanding the dataset to the full range ([Figure 4](#fig4)). We also recorded the responses to all squares of the grid individually (one square input). The one square PSP peak response amplitude with inhibition intact (control) was not distinguishable from that with inhibition blocked (GABAzine) (Materials and methods, [Figure 3—figure supplement 1a](#fig3s1)). As analyzed below ([Figure 6](#fig6)), we find that the apparent lack of effect of GABAzine for very small inputs is because inhibition arrives with a delay that does not affect the peak response of the neuron ([Video 1](#video1)). Since individual neurons may be targeted by more than one grid square ([Figure 1b](#fig1)), individual spots are not completely independent and may interact, especially given the spread in the CA3 pyramidal neuronal arbour. Our analyses show that this interaction does not have a strong or unidirectional effect on the responses of the combinations of squares ([Figure 4—figure supplement 1](#fig4s1), [Figure 5b,d](#fig5)). The ‘observed’ response for a given square combination was plotted against the ‘expected’ response, obtained by linearly summing 1-square responses constituting that combination ([Figure 3b,d](#fig3)). Perfectly linear summation would imply that a multi-square combination of inputs would elicit the same response as the sum of the responses to the individual squares (dotted line, [Figure 3d](#fig3)). [Figure 3e](#fig3) shows responses of a single cell stimulated with 126 distinct 5-square combinations. The ‘observed’ response was sublinear as compared to the ‘expected’ summed response, for most stimuli ([Figure 3e](#fig3)). For all the four measures in [Figure 3c](#fig3), CA3 inputs summed sublinearly at CA1 ([Figure 3e,g](#fig3), [Figure 3—figure supplement 1c](#fig3s1)). At this point, we hypothesised that the observed sublinearity might mostly be due to inhibition divisively scaling excitation, since excitatory and inhibitory conductances were proportional for all stimuli ([Figure 2](#fig2)). We later tested this hypothesis by blocking inhibition ([Figure 5](#fig5)). For all responses measured over all cells, 93.35% responses were individually sublinear, with distribution having mean 0.57 ± 0.31 (SD) ([Figure 3f](#fig3), [Figure 3—figure supplement 1d](#fig3s1)). The slope of the regression line, which indicated the extent of sublinearity, varied between cells, with mean 0.38 ± 0.22 (SD) (n = 33 cells) ([Figure 3g](#fig3)).\n",
"\n",
"\n",
"\n",
"Thus, we found that the CA3-CA1 network exhibits sublinear summation over a large number of inputs.\n",
"\n",
"## CA3-CA1 network performs Subthreshold Divisive Normalization\n",
"\n",
"We then tested how summation sublinearity scaled with a larger range of inputs. We noted that nonlinear functions can be observed better with a large range of inputs (@bib39), and therefore increased the stimulus range ([Figure 4—figure supplement 2](#fig4s2)). GABAergic inhibition has been shown to be responsible for sublinear summation when Schaffer collateral and perforant path inputs are delivered simultaneously to CA1 (@bib16). We hypothesized that the sublinearity within the CA3-CA1 network might also occur due to the effect of inhibition. In general, inhibition may interact with excitation to perform arithmetic operations like subtraction, division, and normalization (@bib11). In order to predict the operation performed by EI integration at the CA3-CA1 network, we created a composite phenomenological model to fit and test for the above three possibilities: subtractive inhibition, divisive inhibition, and divisive normalization ([Equation (1)](#equ1)). We later address the mechanism using a biophysical model ([Figure 6](#fig6)). [Equation (1)](#equ1) describes how inhibition controls the ‘observed’ response ($\\mathbf{\\theta }$) as a function of ‘expected’ response ($\\mathbf{\\epsilon }$), for the above three operations. Alpha (_α_) can be thought to be a subtractive inhibition parameter, beta (_β_) as a divisive inhibition parameter, and gamma (_γ_) a normalization parameter ([Figure 4a](#fig4)).\n",
"\n",
"$$\n",
"{\\displaystyle \\mathit{\\theta }=\\mathit{\\epsilon }-\\frac{\\beta \\mathit{\\epsilon }}{\\gamma +\\mathit{\\epsilon }}\\mathit{\\epsilon }-\\alpha }\n",
"$$\n",
"\n",
"Using the framework of [Equation (1)](#equ1), we asked what computation was performed at the CA3-CA1 network. We recorded from CA1 cells while stimulating CA3 with many combinations of 2, 3, 5, 7 or 9 squares ([Figure 4b](#fig4)). We selected cells with at least 50 input combinations, and pooled responses from all stimuli to a cell. Then, we fit [Equation (1)](#equ1) to the PSP amplitudes ([Figure 4b](#fig4)). From visual inspection, the subtractive inhibition model, $\\mathbf{\\theta }=\\mathbf{\\epsilon }-\\alpha$ (fixing _β_, _γ_ = 0) was a bad fit, since intercepts ($\\alpha$) were close to 0 ([Figure 4a](#fig4)).\n",
"\n",
"figure: Figure 4.\n",
":::\n",
"\n",
"\n",
"### Over a wide input range, integration of CA3 excitatory and feed-forward inhibitory input leads to SDN at CA1.\n",
"\n",
"(**a**) Three phenomenological models of how inhibition interacts with excitation and modulates membrane potential: (left to right) Subtractive Inhibition (SI), Divisive Inhibition (DI) and Divisive Normalization (DN). Note how parameters _α_, _β_ and _γ_ from [Equation (1)](#equ1) affect response output. (**b**) Divisive normalization seen in a cell stimulated with 2, 3, 5, 7 and 9 square combinations. DN and DI model fits are shown in purple and green, respectively. (**c**) Difference in Bayesian Information Criterion (BIC) values for the two models - DI and DN. Most differences between BIC for DI and DN were less than 0, which implied that DN model fit better, accounting for the number of variables used. Insets show raw BIC values. Raw BIC values were consistently lower for DN model, indicating better fit (Two-tailed paired t-test, p<0.00005, n = 32 cells). (**d**) Distribution of the parameter _γ_ of the DN fit for all cells (median = 7.9, n = 32 cells). Compare with **a**, **b** to observe the extent of normalization. (**e**) Distribution of the parameter beta of the DI fit for all cells (mean = 0.5, n = 32 cells). Values are less than 1, indicating sublinear behaviour.\n",
":::\n",
"{#fig4}\n",
"\n",
"figure: Figure 4—figure supplement 1.\n",
":::\n",
"\n",
"\n",
"### Interaction of squares does not affect summation unidirectionally.\n",
"\n",
"(**a**) Example cell showing PPF with electrical, but not with optical stimulation. Individual traces are in grey and black is the average trace. (**b**) Cross Pulse Ratio (Materials and methods) of 25 pairs of stimuli (from five photostimulation squares) presented to an example cell, different from that in **a**. Ratio less than one for self-self pairs, on the diagonal, implies lack of facilitation. (**c**) We restricted our analysis to non-bordering squares and fit the subthreshold divisive normalization model and checked for the value of the normalization parameter (_γ_). The degree of sublinearity and the input-output curve remained unchanged, as indicated by the similarity in the values of DN parameter _γ_, ruling out the hypothesis that interactions between neighboring squares account for the observed sublinearity in the SDN curve. (**d**) Median correlation (0.09) between distributedness of the photostimulation squares and the O/E ratio shows that distances between grid squares do not have a unidirectional relationship with the extent of sublinearity. (**e**) Median correlation (0.02) between how close the patterns were on the grid to the measured voltage response of CA1 again shows no unidirectional effect of physical proximity of patterns.\n",
":::\n",
"{#fig4s1}\n",
"\n",
"figure: Figure 4—figure supplement 2.\n",
":::\n",
"\n",
"\n",
"### Input range expansion for observing nonlinear summation and divisive normalization.\n",
"\n",
"(**a**) A large range of stimulus strengths is required to detect nonlinearity in summation and to characterize divisive normalization. Comparing a nonlinear function of x, f(x) = x^2^ (blue) with a linear function f(x)=0.2 x (orange). It is difficult to distinguish between the two possibilities (from x = 0. to x = 0.2). However, across a larger range of x, differences can be seen. (**b–l**) In order to sample uniformly and over a wider input range of a cell (x axis, labelled Expected sum), we did ‘range expansion’ for some cells (n = 6 cells) as follows. A 37-point grid ([Figure 1—figure supplement 1c](#fig1s1)) was used for these experiments. For a given cell, we first measured the 1-square responses to all 37 grid squares. Then, the 1-square responses were binned into 24 bins, and one member was picked from each bin randomly. In case there were fewer than than one member per bin, other bins were resampled without replacement until the total number was 24. Using this, we generated the distributions of expected sum of N-squares (N = 2, 3, 5, 7, 9). The distributions of the next N-square was truncated such that its minimum began at the maximum of the last N-square distribution (distributions in black, first column of all rows). For example, if the set of squares were 2, 3, 5, 7, 9, we picked the squares such that the distribution of 5 squares started from the end of 3 squares. We then sampled uniformly from this reduced distribution (distributions in grey, second column of all rows). This allowed us to uniformly sample from the entire theoretically possible range, and hence ensure that we observed summation over a wide range of inputs (as explained in **a**). (**m, n**) Responses from one cell without (**m**) and one with (**n**) range expansion are shown. Note that responses from different N-squares (shown in different colors, marked as in [Figure 1d](#fig1)), are overlapping in the Expected-axis for **m**, while the process of range expansion enforces exploration of a large input range in a non-overlapping manner in **n**. (**o**) Reduced chi-square fits confirm DN model fits better than DI. In addition to BIC ([Figure 4c](#fig4)), we used the reduced chi-square fit test to compare the two models of inhibitory computations (divisive inhibition (DI) and divisive normalization (DN)). Divisive normalization model ([Equation 3](#equ3)) almost always fit better than divisive inhibition model ([Equation 2](#equ2)) for reduced chi-square, as it did for BIC ([Figure 4b and c](#fig4)). This figure shows the reduced chi-square ratio plotted for DI and DN models (DI/DN). Values close to one imply good fits, while lesser and greater than one imply over and under-fitting the data respectively. Again, DN is on average closer to one (inset), and the reduced chi-square ratio is almost always larger than one, implying that DN model fit is better than a DI model fit (n = 32 cells).\n",
":::\n",
"{#fig4s2}\n",
"\n",
"By fixing _γ_ and _α_ to 0 in [Equation (1)](#equ1) we obtained the Divisive Inhibition (DI) model. In this form, _β_ can be thought of as I/E ratio. Increasing _β_ decreases the observed response ($\\mathbf{\\theta }$) ([Figure 4a](#fig4)).\n",
"\n",
"$$\n",
"{\\displaystyle \\mathit{\\theta }=\\mathit{\\epsilon }-\\text{}\\beta \\mathit{\\epsilon }}\n",
"$$\n",
"\n",
"Similarly, _β_ was fixed to 1 and _α_ to 0 to get the Divisive Normalization (DN) model. This form of the equation was inspired by the analogous canonical divisive normalization equation for firing rates (@bib11). Here, decrease in _γ_ implies increase in normalization ([Figure 4a](#fig4)).\n",
"\n",
"$$\n",
"{\\displaystyle \\mathit{\\theta }=\\mathit{\\epsilon }-\\frac{\\mathit{\\epsilon }}{\\gamma +\\mathit{\\epsilon }}\\mathit{\\epsilon }=\\text{}\\frac{\\gamma \\mathit{\\epsilon }}{\\gamma +\\mathit{\\epsilon }}}\n",
"$$\n",
"\n",
"We used least-squares polynomial regression to fit DI and DN models to our data. The goodness of fit for all cells was tested by comparing BIC (Bayesian Information Criterion) ([Figure 4c](#fig4)) and reduced chi-squares of the models ([Figure 4—figure supplement 2o](#fig4s2), Materials and methods). DN (_α_ = 0, _β_ = 1) was better than DI (_α_ = 0, _γ_ = 0) model in explaining the data (BIC: Two-tailed paired t-test, p<0.00005, reduced chi-square: Two-tailed paired t-test, p<0.00005, n = 32 cells).\n",
"\n",
"Subthreshold Divisive Normalization (SDN) can be clearly seen in [Figure 4b](#fig4), where observed responses to stimuli with 5 mV and 15 mV expected responses are very similar. This shows that SDN allows CA1 cells to integrate a large range of inputs before reaching spike threshold. Thus, testing with a larger range of inputs showed that the initial finding of constant I/E ratios from [Figure 2](#fig2) needed to be elaborated based on the observed response saturation with increasing input strength. Potential mechanisms for this could be nonlinear summation of excitation and inhibition at the soma (tested in [Figure 5](#fig5)) and inhibitory delays (examined in [Figure 6](#fig6)). In summary, we observed SDN as an outcome of integration of precisely balanced inputs in the CA3-CA1 network.\n",
"\n",
"figure: Figure 5.\n",
":::\n",
"\n",
"\n",
"### Blocking balanced inhibition at resting membrane potential attenuates SDN.\n",
"\n",
"(**a**) Top, schematic of experiment condition. Bottom, a cell showing divisive normalization in control condition. (**b**) Top, schematic of experiment condition with feedforward inhibition blocked (2 uM GABAzine). Bottom, responses of the same cell with inhibition blocked. The responses are much closer to the linear summation line (dashed). The blue lines in **a**, **b** are the fits of the DN model. The value of _γ_ of the fit increases when inhibition is blocked. (**c**) Parameter _γ_ was larger with GABAzine in bath (Wilcoxon rank sum test, p<0.05, n = 8 cells), implying reduction in normalization with inhibition blocked. (**d**) Excitation versus derived inhibition for all points for the cell shown in **a** (area under the curve) (Slope = 0.97, r-square = 0.93, x-intercept = 3.75e-5 mV.ms). Proportionality was seen for all responses at resting membrane potential. Top, ‘Derived inhibition’ was calculated by subtracting control PSP from the excitatory (GABAzine) PSP for each stimulus combination. (**e,f**) R^2^ (median = 0.8) and slope values (median = 0.7) for all cells (n = 8 cells), showing tight IPSP/EPSP proportionality, and slightly more excitation than inhibition at resting membrane potentials.\n",
":::\n",
"{#fig5}\n",
"\n",
"## CA3 feedforward inhibition is necessary for SDN\n",
"\n",
"We first verified our hypothesis that SDN results from feedforward inhibition in the CA3-CA1 network, and not from intrinsic properties of the CA1 neuron. We thus blocked inhibition and repeated the above experiment. We expected that SDN would be lost and linearity would be reinstated upon blocking inhibition.\n",
"\n",
"We recorded responses of CA1 cells to our array of optical stimuli ([Figures 1d](#fig1) and [5a](#fig5)), then applied GABAzine to the bath and repeated the stimulus array ([Figure 5b](#fig5)). We found that when inhibition was blocked, summation approached linearity ([Figure 5b,c](#fig5)). We compared the scaling parameter _γ_ of the divisive normalization model fit, for the above two conditions ([Equation (3)](#equ3)). The values of _γ_ were larger with inhibition blocked, indicative of approach to linearity (Wilcoxon rank-sum test, p<0.05, n = 8 cells) ([Figure 5c](#fig5)). While inhibition accounted in large part for the observed sublinear summation, the cells with inhibition blocked showed some residual sublinearity at high stimulus levels, which has been previously attributed to I~A~ conductance in CA1 neurons (@bib12). Based on the conductance equation ([Equation (5)](#equ5)), leak conductance also contributes in part to the residual sublinearity (Supplementary [Equations (6-8)](<#equ6 equ7 equ8>)). Thus, we confirmed that blocking inhibition reduced sublinearity, attenuating SDN.\n",
"\n",
"## Precise balance is also seen at resting membrane potential\n",
"\n",
"Then, we hypothesised that the membrane potential change evoked by inhibitory synaptic currents could be increasing non-linearly with increasing CA3 input, even though the I/E ratio of conductances would be consistent across the range of input strengths. To address this, we compared responses to identical patterns before and after GABAzine application. For a given cell, for each pattern, we subtracted the initial control response with inhibition intact from the corresponding response with inhibition blocked. This gave us the inhibitory component or ‘derived inhibition’ for each stimulus pattern ([Figure 5d](#fig5), inset). We found that all stimuli to a cell evoked proportional excitation and inhibition even when recorded at resting potential ([Figure 5d,e](#fig5)). Thus, we rejected our hypothesis of non-linear increase in inhibitory post-synaptic potential (IPSP) with excitatory post-synaptic potential at resting membrane potential (EPSP). Over the population, the median slope of the proportionality line was around 0.7, indicating that the EI balance was slightly tilted towards higher excitation than inhibition ([Figure 5f](#fig5)). IPSP/EPSP ratios ([Figure 5f](#fig5)) were smaller than IPSC/EPSC ratios ([Figure 2d](#fig2)) due to proximity of inhibition to its reversal (~−70 mV), than excitation to its reversal (~0 mV), at resting membrane potential (~−65 mV). Overall, we saw precise balance in evoked excitatory and inhibitory synaptic potentials for >100 combinations per neuron.\n",
"\n",
"## Advancing inhibitory onset with increasing input explains SDN\n",
"\n",
"We made a single compartment conductance model ([Figure 6—figure supplement 1a](#fig6s1), [Equation (5)](#equ5)) to analyze the mechanism of SDN. We first show a Hodgkin-Huxley (HH) type single compartment model (Materials and methods), where we have used data from our voltage clamp recordings ([Figure 2](#fig2)), as input to the model. Simulation with both excitation and inhibition produced curves resembling SDN, while only excitation gave a more linear response ([Figure 6a](#fig6), [Figure 6—figure supplement 2](#fig6s2)), hence reproducing the observations depicted in [Figure 5](#fig5). Again, fit parameter _γ_ was significantly higher for the cases without inhibition ([Figure 6—figure supplement 2](#fig6s2), Wilcoxon rank sum test, p<1e-4, n = 13). In order to dissect the mechanism, we wanted to have finer control over synaptic input parameters like kinetics and EI delay.\n",
"\n",
"figure: Figure 6.\n",
":::\n",
"\n",
"\n",
"### Conductance model predicts Excitatory-Inhibitory delay as an important parameter for divisive normalization.\n",
"\n",
"(**a**) Subthreshold responses from HH model, simulated with traces recorded from one voltage clamped cell ([Figure 2](#fig2)). Non-linearly saturating curve, similar to SDN, obtained by simulating with both excitation and inhibition synaptic conductances (black), while the response profile is much more linear with only excitation (red). Each black point is the median response of an excitation trace paired with six different repeats of inhibition for that combination. (**b**) PSP peak amplitude with both excitatory and balanced inhibitory inputs is plotted against the EPSP peak amplitude with only excitatory input. Model showed sublinear behaviour approximating divisive inhibition for I/E proportionality ranging from 1 to 6 when the inhibitory delay was static. Different colours show I/E ratios. (**c**) Same as in b, except the inhibitory delay was varied inversely with excitatory conductance (as shown in **e**). Initial linear zone and diminishing changes in PSP amplitude, indicative of SDN were observed, and the normalization gain was sensitive to the I/E ratio. δ~min~= 2 ms, k = 0.5 nS^−1^, and m = 8.15 ms. Note, the increased overlap in the initial zone (grey box) and the saturation of the PSP peaks in **c**, as compared to **b**. (**d**) Effect of changing EI delay, keeping I/E ratio constant (I/E ratio = 5). Divisive inhibition (green) seen while changing EI delay values from 2 to 10 ms. Divisive normalization (purple) emerges if delays are changed as shown in **e**. δ~min~= 2 ms, k = 0.5 nS^−1^, and m = 8.15 ms. (**e**) Inverse relationship of EI delays with excitation. Inhibitory delay was varied with excitatory conductance in [Equation (4)](#equ4) with δ~min~ = 2 ms, k = 2 nS^−1^, and m = 13 ms. (**f**) Data from an example cell showing the relationship of EI delays with excitation. The relationship is similar to the prediction in **e**. Points are binned averages. Error bars are s.d. (**g**) Data from all cells showing delay as a function of excitation. Different colors indicate different cells (n = 13 cells). Grey lines are linear regression lines through individual cells. (**h**) Traces (from a voltage clamped neuron) showing the decreasing EI delay with increasing amplitude of PSCs. Each trace is an average of 6 repeats.\n",
":::\n",
"{#fig6}\n",
"\n",
"figure: Figure 6—figure supplement 1.\n",
":::\n",
"\n",
"\n",
"### Sensitivity of SDN to EI balance and EI delay, and synaptic time courses used for model.\n",
"\n",
"(**a**) Equivalent circuit for the conductance model showing capacitive, excitatory, inhibitory, and leak components. (**b,c**) Subthreshold divisive normalization was sensitive to I/E ratio (**b**), and EI delay (**c**). SDN was lost when the relationship between I/E ratios for a given cell was permuted b. SDN was also lost when the relationship of delay to excitation ([Figure 6e](#fig6)) was permuted (**c**). (**d**) Histogram of slopes of linear regression lines through EI delays (${\\displaystyle {\\delta }_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}}$) vs excitatory conductance (${\\displaystyle {g}_{\\mathrm{e}\\mathrm{x}\\mathrm{c}}}$) for cells in [Figure 6g](#fig6). (**e**) Rise times (${\\displaystyle {\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}$) and decay times (${\\displaystyle {\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}}$) were extracted from optically evoked postsynaptic currents (EPSCs and IPSCs, [Figure 2](#fig2), [Figure 2—figure supplement 1](#fig2s1)) by fitting a difference of exponentials ([Figure 6](#fig6), [Equation (8)](#equ8)) (n = 13 cells). (**f**) Schematic showing that our observed delay function (right) can be obtained by thresholding of EPSPs (left) in the simple conductance model.\n",
":::\n",
"{#fig6s1}\n",
"\n",
"figure: Figure 6—figure supplement 2.\n",
":::\n",
"\n",
"\n",
"### HH model simulations with voltage clamped data show SDN.\n",
"\n",
"(**a**) Fits of parameter _γ_ ([Equation (3)](#equ3)) to peak Vm vs excitatory conductance for simulations of the HH model with synaptic conductances taken from voltage clamped cells. Black dots are medians for peak depolarization caused by excitatory conductance paired with inhibitory conductance for 6 repeats of the same stimulus. X-axis conductance values were scaled with a factor of 100 for fitting. (**b**) Difference in fit parameter _γ_ for the case with and without inhibition, black dots are means. The two distributions were significantly different (Wilcoxon rank sum test, p<1e-4, n = 13). The fit was susceptible to outliers in some cases.\n",
":::\n",
"{#fig6s2}\n",
"\n",
"With this in mind, we fit a function of difference of exponentials (Materials and methods) to our voltage clamp data to extract the peak amplitudes and kinetics of excitation and inhibition currents (Materials and methods). We used these and other parameters from literature ([Supplementary files 1](#supp1) and [2](#supp2)), and constrained the model to have EI balance, that is have maximum excitatory (${g}_{\\mathrm{e}\\mathrm{x}\\mathrm{c}}$) and inhibitory conductance (${g}_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}$) proportional to each other, with a given I/E ratio. To test for SDN, we simulated our model in the range of experimentally determined I/E ratios, ranging from 1 to 6.\n",
"\n",
"We observed that EI balance with static EI delay led to a slightly sublinear response which can be approximated with a divisive inhibition model ([Figure 6](#fig6)). In contrast, subthreshold divisive normalization (SDN) implies progressively smaller changes in peak PSP amplitude with increase in excitatory input. We surmised that without changing EI balance, SDN should result if the inhibitory onset delays were an inverse function of the excitation ([Figure 6e](#fig6), [Equation (4)](#equ4)). Hence, we simulated the model with dynamic delay, that is with values of inhibitory delay (${\\delta }_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}$) varying as a decreasing function of the excitation.\n",
"\n",
"$$\n",
"{\\delta }_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}={\\delta }_{\\mathrm{m}\\mathrm{i}\\mathrm{n}}+m{e}^{-k{g}_{\\mathrm{e}\\mathrm{x}\\mathrm{c}}}\n",
"$$\n",
"\n",
"Here, ${\\delta }_{\\mathrm{m}\\mathrm{i}\\mathrm{n}}$ is the minimum synaptic delay between excitation and inhibition, k sets the steepness of the delay change with excitation, and _m_ determines the maximum synaptic delay. In [Figure 6c](#fig6), we show that SDN emerged when we incorporated delays changing as a function of the total excitatory input to the model neuron.\n",
"\n",
"We then tested this model prediction. From the EPSC and IPSC curves, we extracted excitatory and inhibitory onsets (Materials and methods), and subtracted the average inhibitory onsets from average excitatory onsets to get inhibitory delay (${\\delta }_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}$) for each stimulus combination. We saw that ${\\delta }_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}$ indeed varied inversely with total excitation (${g}_{\\mathrm{e}\\mathrm{x}\\mathrm{c}}$) ([Figure 6f,g](#fig6)). Notably, the relationship of delay with conductance, with data from all cells pooled, seems to be a single inverse function, and might be a network property ([Figure 6g](#fig6), [Figure 6—figure supplement 1d](#fig6s1)). The input-dependent change in inhibitory delay could be attributed to delayed spiking of interneurons with small excitatory inputs, and quicker firing with larger excitatory inputs. We further illustrate that this delay function emerges naturally by simply applying a threshold to the rising curve of an EPSP at an interneuron ([Figure 6—figure supplement 1f](#fig6s1)). Thus, inhibition clamps down the rising EPSP at progressively earlier times, resulting in saturation of PSP amplitude when excitation is increased ([Figure 6c,d](#fig6), [Figure 8](#fig8)). In [Figure 8a and b](#fig8), we show using a schematic, how SDN emerges when inhibitory onset changes as an inverse function of input strength.\n",
"\n",
"We observed that we were also able to capture the initial linear regime observed in [Figure 4b](#fig4) by using the inverse relationship of delay with excitation in this conductance model. This can be understood as follows: at small excitatory input amplitudes, the EI delay is so large that inhibition arrives too late to affect the peak EPSP. At higher stimulus amplitudes the output response is now subjected to earlier, and hence increasingly effective inhibition, thus flattening the output curve (Appendix 1, [Video 1](#video1), [Figure 3—figure supplement 1a](#fig3s1), [Figure 6c](#fig6)).\n",
"\n",
"We then tested if SDN required both EI balance and dynamic EI delay. We obtained values for balanced ${g}_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}$ for each I/E ratio, and then shuffled the order of the balanced inhibitory vector with the excitation. This implied that the average I/E ratio was maintained over all stimuli, but not for individual stimuli. This shuffled set of inhibitory conductance with respect to excitation was used to calculate ${V}_{\\mathrm{m}\\mathrm{a}\\mathrm{x}}$ ([Figure 6—figure supplement 1b](#fig6s1)). Similarly, we obtained inhibitory delay (${\\delta }_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}$) corresponding to each value of excitation from the dynamic delay curve in [Equation (4)](#equ4) ([Figure 6e](#fig6)). We then shuffled the order of delays, keeping excitation in the same order ([Figure 6—figure supplement 1c](#fig6s1)). In both cases, SDN was strongly attenuated, implying that both EI balance and inverse scaling of inhibitory delay were necessary for SDN ([Figure 6—figure supplement 1b,c](#fig6s1), Supplementary [Equation (6) to (8)](#equ8)). Further, we transformed the membrane current equation ([Equation (5)](#equ5)) into the form which resembles divisive normalization equation (Appendix 1), and saw that in this form, _γ_ depends on the intrinsic properties of the neuron, and is modulated by delays and EI ratios.\n",
"\n",
"Thus, our analysis of a conductance model suggests that SDN could be a general property of balanced feedforward networks, due to two characteristic features: EI balance and inhibitory kinetics. Each of these variables may be subject to plasticity and modulation to attain different amounts of normalization ([Figure 8c,d](#fig8), [Figure 8—figure supplement 1](#fig8s1)).\n",
"\n",
"## Stimulus information is encoded both in amplitude and time\n",
"\n",
"We next asked if the temporally advancing inhibition ([Figure 6e–h](#fig6)) affected PSP peak time with increase in stimulus strength. We calculated the slope of the PSP peak times against the expected axis in the presence (Control) and absence of inhibition (GABAzine) for a given cell. If inhibition cut into excitation and resulted in advancing of peak times with increasing stimulus strength, the slope of peak times would be negative, as shown in [Figure 7a](#fig7). Conversely, when inhibition is blocked, slope of peak times is not expected to change much. We saw that for all cells, slope of the peak time with inhibition intact was lower than the slope in the case with inhibition blocked ([Figure 7b](#fig7)) (Wilcoxon Rank sum test (p=0.006), n = 8 cells).\n",
"\n",
"figure: Figure 7.\n",
":::\n",
"\n",
"\n",
"### Advancing inhibitory onset changes PSP peak time and spike time with increase in stimulus strength.\n",
"\n",
"(**a,b**) The PSP peak arrived earlier following larger input in the control case (black), but not with GABAzine in bath (red). Traces for an example cell, binned (20 bins for Expected sum axis) and averaged, for control (black) and with GABAzine in bath (red). (**c**) Slope of the peak time was more negative in presence of inhibition (control) than when inhibition was blocked (GABAzine) (n = 8 cells). (**d**) Three example traces from the cell in **g** showing the relationship of spikes in presence (black) and absence of inhibition (red). Spikes were produced by HH model, using synaptic conductances from voltage clamp data. The separation between spike times of the two conditions increased with decrease in input conductance (top to bottom). (**e**) Total mutual information of peak amplitude and peak timing with expected sum was not significantly different between Control and GABAzine case (Wilcoxon Rank sum test (<0.05), p=0.11, n = 7 CA1 cells). (**f**) Normalized mutual information between Expected Vm and peak time, Expected Vm and peak amplitude, and conditional mutual information between Expected Vm and peak time, given the knowledge of peak amplitude. Normalized information was calculated by dividing mutual information by total information for each cell (as shown in **d**). Peak times carried more information in the presence of inhibition, and peak amplitudes carried more information in the absence of inhibition. There was higher gain in information about the input with timing if the inhibition was kept intact (Wilcoxon Rank sum test (p<0.05), n = 7 (Pk time, Pk amp) and (p=0.05) n = 6 (Gain with time) CA1 cells). (**g**) Relationship of spike time with excitatory conductance, in the presence (black) and absence of inhibition (red), for HH model simulations. All black points are medians of spikes of each excitation trace paired with six different repeats of inhibition for that combination.\n",
":::\n",
"{#fig7}\n",
"\n",
"figure: Figure 7—figure supplement 1.\n",
":::\n",
"\n",
"\n",
"### Spike time changes with increasing input are steeper in presence of inhibition.\n",
"\n",
"(**a**) Plots from HH model, simulated using data from voltage clamp cells, show the relationship of spike time with conductance, with (black) and without inhibition (red). Each plot is produced using data from an individual cell. Each black point is the median of spike times evoked by excitation paired with six different repeats of inhibition for that stimulus combination. Red points are spike times produced for the same values of excitation, but without any inhibition. Separation of the curves was strongly dependent on the exact value of threshold. Using the natural threshold which emerged from the HH model (−56.3 mV), separation could be seen for about half the cells. (**b,c**) Value of fit parameters k,c when the relationship between spike times and excitatory were fit by c + k/x curve. There was significant difference (Wilcoxon rank sum test, p<0.05, n = 13) between the steepness of the curves for the cases with and without inhibition, where the absence of inhibition makes this curve steeper. There was no significant difference between the intercepts for these fits (Wilcoxon rank sum test, p>0.05, n = 13).\n",
":::\n",
"{#fig7s1}\n",
"\n",
"What does SDN mean for information transmission in balanced networks? While SDN allowed the cell to integrate a large range of inputs before reaching spiking threshold, it also resulted in diminishing changes in PSP peaks at larger inputs ([Figure 4b](#fig4)). This implied that information about the input was partially ‘lost’ from the PSP amplitude. However, PSP times to peak became shorter ([Figure 7a,b](#fig7)), hence potentially encoding some information about the input in this time variable ([Figure 7f](#fig7), [Figure 8b](#fig8)). In contrast, while the peak amplitudes seen with GABAzine predicted the input more reliably, peak times of EPSPs did not change much with input ([Figure 7b,f](#fig7)). Thus, PSP peak time may carry additional information about stimulus strength, when EI balance is maintained.\n",
"\n",
"We quantified this using an information theoretical framework (@bib43). We took linear sum of 1-square PSP peak amplitudes (Expected sum), as a proxy for input strength. We then calculated the mutual information between Expected sum and PSP peak amplitudes of the corresponding N-squares, and between Expected sum and PSP peak timing (Materials and methods). Using this, we asked, how is the information about the input divided between PSP peak amplitude and timing? The total mutual information of both peak amplitude and peak timing with expected sum was slightly lesser in the presence of inhibition, but this difference was statistically not significant ([Figure 7e](#fig7)) (Wilcoxon Rank sum test (<0.05), p=0.11, n = 7 cells). We found that peak timing had more information in presence of inhibition (control), and peak amplitude had more information in absence of inhibition (GABAzine) (Wilcoxon Rank sum test (<0.05), n = 7 cells) ([Figure 7f](#fig7)). Further, we asked, how better can we predict the input, with the knowledge of peak timing, when the peak amplitude is already known? We found that in the presence of inhibition, peak amplitude carried only a part of the total information about the input, and further knowledge of peak time substantially increased the total information. In contrast, in the absence of inhibition, peak amplitude carried most of the information about input, and there was very little gain in information with the knowledge of peak times ([Figure 7f](#fig7)) (Wilcoxon Rank sum test (=0.05), n = 6 cells).\n",
"\n",
"We then asked if the PSP peak time changes are also reflected in spike times. Since most of our stimuli elicited subthreshold responses, studying spiking required an artificial depolarization stimulus. From simulations we found that several parameters of the model (including resting membrane potential, membrane capacitance, synaptic conductances, EI ratio and delay, and spike threshold) could affect the mapping of subthreshold responses to spike timing, suggesting that this is a rich substrate for modulation. Keeping this caveat in mind, we tested the temporal profile of spikes with our model. We let the model cell spike in response to EI (similar to the Control condition) and only E (Gabazine condition). We observed that SDN translated to the spiking domain by encoding stronger stimulus amplitudes as shorter spike latencies, similar to the subthreshold responses. The presence of inhibition decreased the steepness of spike time with conductance ([Figure 7d,g](#fig7), [Figure 7—figure supplement 1](#fig7s1)). The separation between the two conditions was sensitive to the exact value of threshold. At threshold close to resting potential, the separation was low, because the cell spiked before the effect of inhibition set in. For a given threshold, a subset of the cells showed enough separation between conditions ([Figure 7d,g](#fig7), [Figure 7—figure supplement 1](#fig7s1)) and this value could be tuned to obtain maximum separation for each cell.\n",
"\n",
"Overall, these results suggest that with inhibition intact, input information is shared between amplitude and time, and knowledge of peak time and amplitude together contains more information about input than either of them alone.\n",
"\n",
"## Modulation of gating with SDN\n",
"\n",
"We next asked how the two basic parameters - I/E ratio and EI delay - modulated the degree of normalization and kinetics of the SDN curve ([Figure 8c,d](#fig8)). Using our conductance model, we measured the normalization parameter _γ_ (_α_ = 0, _β_ = 1, [Equation (1)](#equ1)) for a range of values of I/E ratio and delays, and found that normalization increased systematically with increase in I/E ratio as well as with increase in the steepness of the EI delay relationship ([Figure 8c](#fig8)). This implies that the degree of normalization of not only an entire neuron, but subsets of inputs to a neuron, could be dynamically altered by changing these parameters. In terms of gating, for a neuron with all inputs tightly balanced, any subset of inputs with reduction in I/E ratio will be gated ‘on’, corresponding to a condition of higher _γ_. Neurons can thus differentially gate and respond to specific inputs, while still retaining the capacity to respond to other input combinations.\n",
"\n",
"figure: Figure 8.\n",
":::\n",
"\n",
"\n",
"### Emergence of SDN from balanced excitation and inhibition, coupled with dynamic EI delays.\n",
"\n",
"(**a**) Schematic showing precisely balanced EPSPs (blue) and corresponding IPSPs (red) summing to produce PSPs (purple). The EPSPs and IPSPs increase in equal input steps. (**b**) Zooming into the portion in the rectangle in **a**. Excitation onset is constant, but inhibition onset changes as an inverse function of input or conductance (${\\displaystyle {g}_{\\mathrm{e}\\mathrm{x}\\mathrm{c}}}$), as shown in [Figure 6](#fig6). With increasing input, inhibition arrives earlier and cuts into excitation earlier for each input step. This results in smaller differences in excitatory peaks with each input step, resulting in SDN. The timing of PSP peaks (purple) becomes progressively advanced, whereas the timing of EPSP peaks (blue) does not, consistent with our results in [Figure 7](#fig7). (**c,d**) Normalization as a function of the two building blocks – EI balance (I/E ratio) and EI delays (interneuron recruitment kinetics, k, as predicted by the model. Larger values of both imply greater normalization and increased gating. Colors of the SDN curves depict the value of gamma (_γ_), as shown in the phase plot in **d**. White squares are values of _γ_ larger than 40, where almost no normalization occurs.\n",
":::\n",
"{#fig8}\n",
"\n",
"figure: Figure 8—figure supplement 1.\n",
":::\n",
"\n",
"\n",
"### PSP traces showing the effect of I/E ratio and inhibitory recruitment kinetics (k) on SDN.\n",
"\n",
"(**a**) Similar to [Figure 8c](#fig8), but with the simulated membrane potential traces replacing the Input-Output curves. Colorbar represents values of _γ_. (**b**) Input range expansion by SDN. Schematic shows comparison of various models of EI interaction for the range of inputs accommodated by a neuron, before reaching spike threshold.\n",
":::\n",
"{#fig8s1}\n",
"\n",
"# Discussion\n",
"\n",
"This study describes two fundamental properties of the CA3-CA1 feedforward circuit: balanced excitation and inhibition from arbitrary presynaptic CA3 subsets, and an inverse relationship of excitatory-inhibitory delays with CA3 input amplitude. We used optogenetic photostimulation of CA3 with hundreds of unique stimulus combinations and observed precise EI balance at individual CA1 neurons for every input combination. Stronger stimuli from CA3 led to proportional increase in excitatory and inhibitory amplitudes at CA1, and a decrease in the delay with which inhibition arrived. Consequently, larger CA3 inputs had shorter inhibitory delays, which led to progressively smaller changes in CA1 membrane potential. We term this gain control mechanism Subthreshold Divisive Normalization (SDN). This reduction in inhibitory delay with stronger inputs contributes to a division of input strength coding between PSP amplitude and PSP timing.\n",
"\n",
"## Precise balance in the hippocampus\n",
"\n",
"Our findings demonstrate that precise EI balance is maintained by arbitrary combinations of neurons in the presynaptic network, despite the reduced nature of the slice preparation, with no intrinsic network dynamics. This reveals exceptional structure in the connectivity of the network. Theoretical analyses suggest that networks can achieve detailed balance with inhibitory Spike Timing Dependent Plasticity (iSTDP) rules (@bib23; @bib34; @bib46). Such an iSTDP rule has been observed in the auditory cortex (@bib13). Given that balance needs to be actively maintained (@bib51), we suspect that similar plasticity rules (@bib23) may also exist in the hippocampus.\n",
"\n",
"Precisely balanced networks, with all input subsets balanced, are well suited for input gating (@bib4; @bib23). The finding that most silent CA1 cells can be converted to place cells for arbitrary locations predicts the existence of an input gating mechanism (@bib32), but the nature of this mechanism remains unknown. One prediction of precise balance is that inputs for multiple potential place fields may be balanced, and hence place field activity is gated ‘off’. Evoked depolarizations (@bib32) or dendritic plateau potentials (@bib8; @bib9), which potentiate the subset of active synapses, that is, change the I/E ratio (@bib20), can flip the gate ‘on’, thereby converting a silent cell to a place cell for that specific place field. This reasoning corroborates the observation of homogenous inhibition suppressing out-of-field heterogeneously tuned excitation (@bib20), while providing a finer, synaptic scale view of the gating mechanism.\n",
"\n",
"## EI delays and temporal coding\n",
"\n",
"In several EI networks in the brain, inhibition is known to suppress excitation after a short time delay, leaving a ‘window of opportunity’ for spiking to occur (@bib24; @bib41; @bib48). We have shown that balanced inhibitory input arrives with a delay modulated by the excitatory input in a feedforward circuit. This inverse relationship of EI delay with excitation has not been explicitly shown, although @bib22 report a decrease in EI delays with increase in whisker stimulation speed in layer 4 cells. We show that modulation of EI delay by excitation helps encode the input information in both amplitude and timing of the PSP ([Figure 7](#fig7)). Thus, large inputs could be represented with fewer spikes, while conserving input strength information in spike timing. In CA1, a classic example of such dual coding is theta phase precession (@bib28). In addition, spike times during sharp wave ripples, gamma oscillations and time cell representations are also precise up to ~10 ms, which is the range of the dynamic ‘window of opportunity’ we observe. This dynamic window also implies that the neuron can transition from temporal integration mode at small input amplitudes to coincidence detection at large input amplitudes (@bib19; @bib24; @bib48). Consistent with this range of spike-coding transformations, our simulations suggest that the precise mapping of subthreshold summation to spike timing can be effectively modulated by several cellular parameters as well as by details of input activity ([Figure 7d,g](#fig7)).\n",
"\n",
"## Subthreshold Divisive Normalization (SDN): a novel gain control mechanism\n",
"\n",
"We have introduced Subthreshold Divisive normalization (SDN) as a novel gain control mechanism arising from EI balance and dynamic EI delays. Our study was uniquely able to observe SDN because of the large range of inputs possible (@bib39) using patterned optical stimulation. While we observed no unidirectional correlation of the distance between input spots and their responses for most inputs ([Figure 4—figure supplement 1](#fig4s1)), a limitation of this stimulation design is that some of the inputs may not be fully independent due physical proximity of stimulus spots. SDN expands the dynamic range of inputs that a neuron can accommodate before reaching spike threshold ([Figure 8—figure supplement 1b](#fig8s1)). This is particularly useful for temporally coding, sparsely spiking neurons like CA1 (@bib1). So far, analogous gain control by divisive normalization has only been observed for firing rates of neurons (@bib11). This implies that the timescales of gain change in DN are averaged over periods of tens of milliseconds, over which rates change. As opposed to this, in SDN, gain of every input is normalized at synaptic (millisecond) timescales. Our results add a layer of subthreshold gain control in single neurons, to the known suprathreshold gain control at the population level in CA1 (@bib40). This two-step gain control implies that the dynamic range of the population may be wider than previously estimated. While most experimental observations of firing rate gain change have been explained by the phenomenological divisive normalization equation, the mechanistic basis for normalization has been unclear. SDN provides a biophysical explanation for phenomenological divisive normalization by connecting EI ratios and delays with gain control.\n",
"\n",
"I/E ratio can be changed by neuromodulation (@bib18; @bib17), by short term plasticity mechanisms (@bib5; @bib30; @bib45) and by disinhibition (@bib6). Although we show that EI delays are input amplitude dependent, they may also be modulated by external signals, or behavioural states such as attention (@bib29) ([Figure 8c,d](#fig8)). Such interneuron recruitment based changes have been shown to exist in thalamocortical neurons (@bib19). Dynamic regulation of EI delay has been theoretically explored in balanced networks (@bib10; @bib31) for temporal gating of transient inputs independently by amplitude and time. Thus, temporal gating by EI delays (@bib31), combined with the amplitude gating by detailed balance (@bib47) could be a powerful mechanism for gating signals (@bib31) in the hippocampal feedforward microcircuit.\n",
"\n",
"Several studies point toward the existence of precise EI balance in the cortex (@bib3; @bib37; @bib48; @bib49; @bib53; @bib54), and here we have shown it in the hippocampus. We propose that input strength dependent inhibitory delay change may be a general property of feedforward network motifs. Together, these suggest that precisely balanced feedforward networks are elegantly suited for controlling gain, timing and gating at individual neurons in neural circuits.\n",
"\n",
"# Materials and methods\n",
"\n",
"table: Key resources table\n",
":::\n",
"| Reagent type (species) or resource | Designation | Source or reference | Identifiers | Additional information |\n",
"| --------------------------------------------------- | --------------------------------------- | ------------------- | ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- | -------------------------------------- |\n",
"| Genetic reagent (_M. musculus_) | C57BL/6-Tg(Grik4-cre)G32-4Stl/J | Jackson Laboratory | Stock #:006474 | Dr. Susumu Tonegawa's laboratory, MIT |\n",
"| Strain, strain background (Adeno-associated virus) | AAV5.CAGGS.Flex.ChR2-tdTomato.WPRE.SV40 | Penn Vector Core | | |\n",
"| Software, algorithm | MOOSE simulator | @bib42 | [RRID:](https://www.google.com/url?q=https://github.com/BhallaLab/moose&sa=D&source=hangouts&ust=1552414403557000&usg=AFQjCNFghrTW2n5GnP6DqkvCmllzt7bwXQ)[SCR_008031](https://scicrunch.org/resolver/SCR_008031) | Dr. Upinder Bhalla's laboratory, NCBS |\n",
":::\n",
"{#keyresource}\n",
"\n",
"## Animals\n",
"\n",
"All experimental procedures were approved by the National Centre for Biological Sciences Institutional Animal Ethics Committee (Protocol number USB–19–1/2011), in accordance with the guidelines of the Government of India (animal facility CPCSEA registration number 109/1999/CPCSEA) and equivalent guidelines of the Society for Neuroscience. CA3-cre (_C57BL/6-Tg (Grik4-cre) G32-4Stl/J_ mice, Stock number 006474) were obtained from Jackson Laboratories. The animals were housed in a temperature controlled environment with a 14 hr light: 10 hr dark cycle, with _ad libitum_ food and water.\n",
"\n",
"## Virus injections\n",
"\n",
"21–30 days old male transgenic mice were injected with Lox-ChR2 (AAV5.CAGGS.Flex.ChR2-tdTomato.WPRE.SV40) virus obtained from University of Pennsylvania Vector Core. Injection coordinates used were −2.0 mm RC, ±1.9 mm ML, −1.5 mm DV. ~300–400 nl solution was injected into the CA3 region of left or right hemisphere with brief pressure pulses using Picospritzer-III (Parker-Hannifin, Cleveland, OH). Animals were allowed to recover for at least 4 weeks following surgery.\n",
"\n",
"## Slice preparation\n",
"\n",
"8–12 week (4–8 weeks post virus injection) old mice were anesthetized with halothane and decapitated post cervical dislocation. Hippocampus was dissected out and 350 um thick transverse slices were prepared. Slices (350 microns) were cut in ice-cold high sucrose ASCF containing (in mM): 87 NaCl, 2.5 KCl, 1.25 NaH~2~PO~4~, 25 NaHCO~3~, 75 sucrose, 10 glucose, 0.5 CaCl~2~, 7 MgCl~2~. Slices were stored in a holding chamber, in artificial cerebro-spinal fluid (aCSF) containing (in mM) - 124 NaCl, 2.7 KCl, 2 CaCl~2~, 1.3 MgCl~2~, 0.4 NaH~2~PO~4~, 26 NaHCO~3~, and 10 glucose, saturated with 95% O~2~/5% CO~2~. After at least an hour of incubation, the slices were transferred to a recording chamber and perfused with aCSF at room temperature.\n",
"\n",
"## Electrophysiology\n",
"\n",
"Whole cell recording pipettes of 2-5MO were pulled from thick-walled borosilicate glass on a P-97 Flaming/Brown micropipette puller (Sutter Instrument, Novato, CA). Pipettes were filled with internal solution containing (in mM): 130 K-gluconate, 5 NaCl, 10 HEPES, 1 EGTA, 2 MgCl~2~, 2 Mg-ATP, 0.5 Na-GTP and 10 Phosphocreatinine, pH adjusted to 7.3, osmolarity ~285 mOsm. The membrane potential of CA1 cells was maintained near −65 mV, with current injection, if necessary. GABA-A currents were blocked with GABAzine (SR-95531, Sigma) at 2 uM concentration for some experiments. Cells were excluded from analysis if the input resistance changed by more than 25% (measured for 15/39 cells) or if membrane voltage changed more than 2.5 mV (measured for 39/39 cells, maximum current injected to hold the cell at the same voltage was ±15 pA) of the initial value. For voltage clamp recordings, the K-gluconate was replaced by equal concentration Cs-gluconate. Cells were voltage clamped at 0 mV (close to calculated excitation reversal) and −70 mV (calculated inhibition reversal) for IPSC and EPSC recordings respectively. At 0 mV a small component of APV sensitive inward current was observed, and was not blocked during recordings. Cells were excluded if series resistance went above 25MO or if it changed more than 30% of the initial value, with mean series resistance being 15.7MO ± 4.5 MO s.d. (n = 13). For CA3 current clamp recordings, the cells were excluded if the V~m~ changed by 5 mV of the initial value. For whole-cell recordings, neurons were visualized using infrared microscopy and differential interference contrast (DIC) optics on an upright Olympus BX61WI microscope (Olympus, Japan) fitted with a 40X (Olympus LUMPLFLN, 40XW), 0.8NA water immersion objective. Recordings were acquired on a HEKA EPC10 double plus amplifier (HEKA Electronik, Germany) and filtered 2.9 kHz and digitized at 20 kHz.\n",
"\n",
"## Optical stimulation setup\n",
"\n",
"Optical stimulation was done using DMD (Digital Micromirror Device) based Optoma W316 projector (60 Hz refresh rate) with its color wheel removed. Image from the projector was miniaturized using a Nikon 50 mm f/1.4D lens and formed at the focal plane of the tube lens, confocal to the sample plane. The white light from the projector was filtered using a blue filter (Edmund Optics, 52532), reflected off of a dichroic mirror (Q495LP, Chroma), integrated into the light path of the Olympus microscope, and focused on sample through a 40X objective. This arrangement covered a circular field of around 200 micron diameter on sample. 2.5 pixels measured one micron at sample through the 40X objective. Light intensity, measured using a power meter, was about 15 mW/mm^2^ at sample surface. Background light from black screen usually elicited no or very little synaptic response at recorded CA1 cells. A shutter (NS15B, Uniblitz) was present in the optical path to prevent the slice from being stimulated by background light during the inter-trial interval. The shutter was used to deliver stimulus of 10–15 ms per trial. A photodiode was placed in the optical path after the shutter to record timestamps of the delivered stimuli.\n",
"\n",
"## Patterned optical stimulation\n",
"\n",
"Processing 2 was used for generating optical patterns. All stimuli were 16 micron squares subsampled from a grid. 16 micron was chosen since it is close to the size of a CA3 soma. The light intensity and square size were standardized to elicit typically one spike per cell per stimulus. The number of spikes varied to some extent based on the expression of ChR2, which varied from cell to cell. The switching of spots from one trial to next, at 3 s inter trial interval, prevented desensitization of ChR2 over successive trials ([Figure 1g](#fig1)).\n",
"\n",
"For a patched CA1 cell, the number of connected CA3 neurons stimulated per spot was estimated to be in the range of 0 to a maximum of 50 for responses ranging from 0 to 2 mV. These calculations assume a contribution of 0.2 mV per synapse (@bib35) and release probability of ~0.2 (@bib36). This number includes responses from passing axons, which could also get stimulated in our preparation.\n",
"\n",
"We did not observe any significant cross stimulation of CA1 cells. CA1 cells were patched and the objective was shifted to the CA3 region of the slice, where the optical patterns were then projected. CA1 cells showed no response to optical stimulation because (i) Use of CA3-cre line restricted ChR2 to CA3 cells, (ii) physical shifting of the objective away from CA1 also made sure that any leaky expression, if present, did not elicit responses. Using a cre-based targeted optogenetic stimulation combined with patterned optical stimulation, we designed an experiment which was both more specific and more effective at exploring a large stimulus space. Unlike electrical stimulation, optical stimulation specifically excited CA3 pyramidal neurons, and hence the recorded inhibition was largely feedforward. We believe this specificity was crucial to the finding that I/E ratios for all stimuli to a cell are conserved. Electrical stimulation does not distinguish between neuronal subclasses, and in particular fails to separate out the inhibitory interneurons. Since a key part of our findings emerged from being able to establish a temporal sequence of activation of interneurons, it was crucial to exclude monosynaptic stimulation of interneurons in the experimental design. Second, patterned optical stimulation allowed us to explore a grid of 225 stimulus points in CA3, thereby obtaining a wide array of stimulus combination with large dynamic range, without compromising on the specificity of stimulation ([Figure 1](#fig1), [Figure 1—figure supplement 1](#fig1s1)).\n",
"\n",
"We used four different stimulus grids ([Figure 1—figure supplement 1](#fig1s1)). All squares from a grid were presented individually (in random order) and in a stimulus set - randomly chosen combinations of 2, 3, 5, 7, or 9, with 2, 3 or 6 repeats of each combination. The order of presentation of a given N square combination was randomized from cell to cell.\n",
"\n",
"## Data analysis and code availability\n",
"\n",
"All analyses were done using custom written software in Python 2.7.12 (numpy, scipy, matplotlib and other free libraries) and MatlabR2012b. All error bars are standard deviations. All analysis codes are available as a free library at (<https://github.com/sahilm89/linearity>; copy archived at <https://github.com/elifesciences-publications/linearity>).\n",
"\n",
"## Pre-processing\n",
"\n",
"PSPs and PSCs were filtered using a low-pass Bessel filter at 2 kHz, and baseline normalized using 100 ms before the optical stimulation time as the baseline period. Period of interest was marked as 100 ms from the beginning of optical stimulation, as it was the typical timescales of PSPs. Timing of optical stimulation was determined using timestamp from a photodiode responding to the light from the projector. Trials were flagged if the PSP in the interest period were indistinguishable from baseline period due to high noise, using a two sample KS test (p-value<0.05). Similarly, action potentials in the interest period were flagged and not analyzed, unless specifically mentioned.\n",
"\n",
"## Feature extraction\n",
"\n",
"A total of four measures were used for analyzing PSPs and PSCs ([Figure 3c](#fig3)). These were mean, area under the curve, average and area to peak. This was done to be able to catch differences in integration at different timescales, as suggested by @bib39. Trials from CA1 were mapped back to the grid locations of CA3 stimulation for comparison of Expected and Observed responses. Grid coordinate-wise features were calculated by averaging all trials for a given grid coordinate.\n",
"\n",
"## Subthreshold divisive normalization model\n",
"\n",
"Different models of synaptic integration: Subtractive Inhibition, Divisive Inhibition, and Divisive Normalization models were obtained by constraining parameters in [Equation (1)](#equ1). The models were then fit to the current clamp dataset using lmfit. Reduced chi-squares ([Figure 4—figure supplement 2o](#fig4s2)) and Bayesian Information Criterion ([Figure 4c](#fig4)) were used to evaluate the goodness of fits of these models to experimental data.\n",
"\n",
"## Single-compartment model\n",
"\n",
"A single-compartment conductance-based model was created in Python using sympy and numpy. The model consisted of leak, excitatory and inhibitory synaptic conductances ([Equation (5)](#equ5), [Figure 6—figure supplement 1a](#fig6s1)) to model the subthreshold responses by the CA1 neurons.\n",
"\n",
"$$\n",
"{C}_{m}\\frac{d{V}_{m}}{dt}={g}_{\\mathrm{l}\\mathrm{e}\\mathrm{a}\\mathrm{k}}\\left({V}_{m}-{E}_{\\mathrm{l}\\mathrm{e}\\mathrm{a}\\mathrm{k}}\\right)+{g}_{\\mathrm{e}\\mathrm{x}\\mathrm{c}}\\left({V}_{m}-{E}_{\\mathrm{e}\\mathrm{x}\\mathrm{c}}\\right)+{g}_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}\\left({V}_{m}-{E}_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}\\right)\n",
"$$\n",
"\n",
"The parameters used for the model were taken directly from data, or literature ([Supplementary file 2](#supp2)). The synaptic conductances ${g}_{\\mathrm{e}\\mathrm{x}\\mathrm{c}}\\left(t\\right)$, and ${g}_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}\\left(t\\right)$ were modeled as difference of exponentials ([Equations (6) and (7)](<#equ6 equ7>)):\n",
"\n",
"$$\n",
"{\\displaystyle {\\displaystyle {g}_{\\mathrm{e}\\mathrm{x}\\mathrm{c}}\\left(t\\right)={\\displaystyle {\\stackrel{-}{g}}_{\\mathrm{e}\\mathrm{x}\\mathrm{c}}\\left(\\frac{{e}^{\\left(\\frac{-t}{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}}\\right)}-{e}^{\\left(\\frac{-t}{{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}\\right)}}{-{\\left(\\frac{{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}}\\right)}^{\\frac{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}}{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}-{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}}+{\\left(\\frac{{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}}\\right)}^{\\frac{{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}-{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}}}\\right)}}}\n",
"$$\n",
"\n",
"$$\n",
"{\\displaystyle {\\displaystyle {g}_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}\\left(t\\right)={\\displaystyle {\\stackrel{-}{g}}_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}\\left(\\frac{{e}^{\\left(\\frac{{\\delta }_{\\mathrm{i}\\mathrm{n}{\\mathrm{h}}^{-\\mathrm{t}}}}{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}}\\right)}-{e}^{\\left(\\frac{{\\delta }_{\\mathrm{i}\\mathrm{n}{\\mathrm{h}}^{-\\mathrm{t}}}}{{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}\\right)}}{-{\\left(\\frac{{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}}\\right)}^{\\frac{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}}{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}-{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}}+{\\left(\\frac{{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}}\\right)}^{\\frac{{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}-{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}}}\\right)}}}\n",
"$$\n",
"\n",
"For the divisive normalization case, the inhibitory delays (${\\delta }_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}$) were modeled to be an inverse function of ${g}_{\\mathrm{e}\\mathrm{x}\\mathrm{c}}\\left(t\\right)$ ([Equation (4)](#equ4)). In other cases, they were assumed to be constant and values were taken from [Supplementary file 2](#supp2).\n",
"\n",
"## HH-based single-compartment model\n",
"\n",
"A single-compartment Hodgkin Huxley model with parameters mentioned in [Supplementary file 3](#supp3) was simulated in MOOSE (@bib42) to analyze how measured synaptic conductances sum to cause CA1 somatic depolarization. To enable spiking, we included sodium and potassium delayed rectifier (KDR) channels in these neurons. Then, we drove this neuron with synaptic input as measured from voltage clamp data.\n",
"\n",
"## Measurement of synaptic conductances\n",
"\n",
"We calculated excitatory and inhibitory conductances using [Equation (5)](#equ5), while holding the neuron at inhibitory and excitatory reversal potentials respectively (@bib54, @bib3). To measure excitatory conductance (${g}_{\\mathrm{e}\\mathrm{x}\\mathrm{c}}$), we clamped the membrane to the inhibitory reversal potential (${E}_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}$). In the absence of a stimulus, the holding current gave us the value of leak current (${I}_{\\mathrm{l}\\mathrm{e}\\mathrm{a}\\mathrm{k}}$). Excitatory synaptic current (${I}_{\\mathrm{e}\\mathrm{x}\\mathrm{c}}$) was measured as the change in membrane current evoked by the input stimulus (${I}_{m}-{I}_{\\mathrm{l}\\mathrm{e}\\mathrm{a}\\mathrm{k}}$), from the baseline of holding current. We calculated the ${g}_{\\mathrm{e}\\mathrm{x}\\mathrm{c}}$ by dividing this stimulus evoked excitatory current by the excitatory driving force (${V}_{m}-{E}_{\\mathrm{e}\\mathrm{x}\\mathrm{c}}$). The same procedure was repeated at excitatory reversal to measure inhibitory conductance (${g}_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}$) for each stimulus.\n",
"\n",
"With this method, measurement of ${g}_{\\mathrm{e}\\mathrm{x}\\mathrm{c}}$ and ${g}_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}$ at corresponding clamped membrane voltages was independent of the absolute value of ${I}_{\\mathrm{l}\\mathrm{e}\\mathrm{a}\\mathrm{k}}$. However, we needed to obtain an estimate of leak conductance ${(g}_{\\mathrm{l}\\mathrm{e}\\mathrm{a}\\mathrm{k}})$ for the purposes of the model ([Equation (5)](#equ5)). We could not use the absolute value of ${I}_{\\mathrm{l}\\mathrm{e}\\mathrm{a}\\mathrm{k}}$ as measured in our voltage clamped neurons because of blockage of potassium channels with Cs internals. Hence, for use in our conductance model, ${g}_{\\mathrm{l}\\mathrm{e}\\mathrm{a}\\mathrm{k}}$ measurements were not taken from our voltage clamp data, and instead the value was taken from literature.\n",
"\n",
"## Fitting data\n",
"\n",
"Voltage clamp data was fit to a difference of exponential functions ([Equation (8)](#equ8), [Figure 6—figure supplement 1e](#fig6s1)) by a non-linear least squares minimization algorithm using lmfit, a freely available curve fitting library for Python. Using this, we obtained amplitudes ($\\overline{g}$), time course (${\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}$, ${\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}$) and onset delay from stimulus (${\\delta }_{\\mathrm{o}\\mathrm{n}\\mathrm{s}\\mathrm{e}\\mathrm{t}})$ for both excitatory and inhibitory currents. We then calculated inhibitory onset delay (${\\delta }_{\\mathrm{i}\\mathrm{n}\\mathrm{h}}$) by subtracting onset delay of excitatory from inhibitory traces.\n",
"\n",
"$$\n",
"{\\displaystyle {\\displaystyle g\\left(t\\right)=\\overline{g}\\left(\\frac{{e}^{\\left(\\frac{{\\delta }_{\\mathrm{o}\\mathrm{n}\\mathrm{s}\\mathrm{e}\\mathrm{t}}-t}{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}}\\right)}-{e}^{\\left(\\frac{{\\delta }_{\\mathrm{o}\\mathrm{n}\\mathrm{s}\\mathrm{e}\\mathrm{t}}-t}{{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}\\right)}}{-{\\left(\\frac{{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}}\\right)}^{\\frac{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}}{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}-{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}}+{\\left(\\frac{{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}}\\right)}^{\\frac{{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}{{\\tau }_{\\mathrm{d}\\mathrm{e}\\mathrm{c}\\mathrm{a}\\mathrm{y}}-{\\tau }_{\\mathrm{r}\\mathrm{i}\\mathrm{s}\\mathrm{e}}}}}\\right)}}\n",
"$$\n",
"\n",
"## Onset detection\n",
"\n",
"Onsets were also detected using three methods. Since we propose onset delays to be a function of the excitation peak, we avoided onset finding methods such as time to 10% of peak, which rely on peaks of the PSCs. We used threshold based (time at which the PSC crossed a threshold), slope based (time at which the slope of the PSC onset was the steepest) and a running window based method. In the running window method, we ran a short window of 0.5 ms, and found the time point at which distributions of two consecutive windows became dissimilar, using a two sample KS test. Ideally, with no input, the noise distribution across two consecutive windows should remain identical. All three methods gave qualitatively similar results.\n",
"\n",
"## Modeling detailed balanced synapses\n",
"\n",
"Synaptic inputs were modeled as sums of probabilistically activated basal synapses with synaptic strengths taken from a lognormal distribution with shape and scale parameters as given by our one square current clamp data (shape = −0.39, scale = 0.80). The width of the weight distribution was altered by changing the scale parameter. Probabilistic synaptic activation was modeled as a binomial process, with synaptic ‘release probability’ for excitatory and inhibitory inputs set at 0.2 and 0.8, respectively.\n",
"\n",
"Inhibitory inputs were generated with various degrees of correlation to the excitation, by shuffling the excitatory weights in differently sized bins, from one to the length of the excitatory weight vector, controlled by a parameter _ρ_. In this manner, as _ρ_ changed from 1 to 0, excitatory and inhibitory weight vectors changed from paired (detailed balance) to completely unpaired but with identical mean and variance of the weight distributions (global balance).\n",
"\n",
"These synapses could be engaged by delivering stimuli, with the number of synapses per stimulus sampled from a Poisson distribution with mean of 5 synapses per stimulus. The total number of excitatory and inhibitory synaptic inputs engaged by a stimulus were always identical. Each stimulus was repeated six times. The resultant means and standard deviations for excitatory and inhibitory inputs were plotted against each other to compare different degrees of correlation. The whole process was repeated 100 times, and correlations and r-squared values were averaged to generate the heatmaps.\n",
"\n",
"## Mutual information calculation\n",
"\n",
"Mutual information was calculated by non-parametric entropy estimation and histogram methods. NPEET (<https://github.com/gregversteeg/NPEET>) was used for non-parametric estimation of Mutual Information. The relationship between variables was shuffled 500 times to find the significance of the Mutual Information estimate. If the true value of MI was not larger than 90% of the distribution obtained by shuffling, mutual information was assumed to be 0. If the total information about the linear sum of one square responses using both peak amplitude and time could not be established with 90% confidence as described above, the cell was excluded from further analysis. We also used the histogram method to find the mutual information (data not shown), and saw a similar trend. Cells with fewer than 80 trials and less than 2 mV inter-quartile range in the linear sum from one square PSP were excluded from the analysis. The calculated linear sum from one square PSP peak amplitude responses, measured N-square peak amplitudes and time were binned with an equal number of bins. The number of bins were calculated using Sturges’ Rule, which selects the number of bins as 1 + 3.3 log n, where n is the total number of observations for a given neuron. Bin frequencies were divided by the total number of responses to get the probability of occurrence p(x) of each bin.\n",
"\n",
"Mutual Information was then calculated for all pairs of combinations between linear sum, peak amplitude and time using [Equation (9)](#equ9) and [(10)](#equ10).\n",
"\n",
"$$\n",
"MI\\left(X,Y\\right)=H\\left(X\\right)+H\\left(Y\\right)-H\\left(X,Y\\right)\n",
"$$\n",
"\n",
"where Shannon’s entropy $H\\left(X\\right)$ for a variable $X$, is given as:\n",
"\n",
"$$\n",
"{\\displaystyle H\\left(X\\right)=\\sum _{x\\phantom{\\rule{0.1667em}{0ex}}ϵX}-\\mathrm{p}(\\mathrm{x})\\mathrm{l}\\mathrm{o}\\mathrm{g}2\\mathrm{p}(\\mathrm{x})}\n",
"$$\n",
"\n",
"Further, conditional mutual Information was calculated to measure gain in information about input (linear sum) by knowledge of peak timing when peak amplitude is already known. It was calculated using [Equation 11](#equ11).\n",
"\n",
"$$\n",
"I\\left(X;Y|Z\\right)=H\\left(X,Z\\right)+H\\left(Y,Z\\right)-H\\left(X,Y,Z\\right)-H\\left(Z\\right)\n",
"$$\n",
"\n",
"Normalized mutual information was calculated by dividing mutual information between pairs of variables by the total information between all three variables ([Equation 12](#equ12)).\n",
"\n",
"$$\n",
"I\\left(X;Y,Z\\right)=H\\left(Z\\right)+H\\left(X,Y\\right)-H\\left(X,Y,Z\\right)\n",
"$$\n",
"\n",
"## Cross-pulse adaptation\n",
"\n",
"We individually presented five unique photostimulation spots in all possible pairwise combinations, with an inter-stimulus interval of 50 ms (@bib15), to test for the interaction using a Cross Pulse Adaptation protocol. We then compared the averages of ten repeats of the response for a given spot when it arrived second in the stimulus-pair, to when it came first. Hence, if there is facilitation caused due to the presence of the first spot, then we should observe that the response to the spot when it comes second is larger than when it comes first in the stimulus pair. To quantify this change, we calculated the ratio between the average response of the spot, when it arrives at the second place, to the response when it arrives at the first place. This gave us the Cross Pulse Ratio ([Figure 4—figure supplement 1b](#fig4s1)). A necessary internal control was that the self-self spot pairs should get facilitated. However, we observed lack of facilitation for self-self pairs, for all the cells we tested ([Figure 4—figure supplement 1](#fig4s1), n = 6 cells). To ensure that this effect was not due to a limitation of the preparation, we tested paired pulse facilitation with electrical stimulation on the same neuron which depressed with optical stimulation. We show that the neuron shows PPF with electrical, but not with optical stimulation ([Figure 4—figure supplement 1a](#fig4s1)). Unlike electrical stimulation, which strongly and briefly stimulates many axonal fibres, optical stimulation targets neurons with varying degrees of strengths, and incomplete recovery of ChR2 from desensitization at such short timescales may be the reason for the second pulse not being as effective as the first one. This interfered with our ability to measure paired pulse facilitation and introduced uncertainty in interpreting cross-pulse effects. This precluded further investigation using this approach.\n",
"\n",
"## Distributedness and physical distance between square patterns\n",
"\n",
"We calculated the effect of the interaction due to physical proximity of photostimulation squares on the responses. We defined a quantity distributedness, as the sum of the distance between all simultaneously stimulated spots from the combined centre of mass of these spots ([Figure 4—figure supplement 1d](#fig4s1)). We compared this to degree of sublinearity, that is the ratio between the Observed response (O) and the Expected sum (E) of individual squares. Thus, if the interaction between neighbouring squares caused sublinearity, we would see a positive correlation between the distributedness and O/E ratio (for the stimuli within an N-square set). Conversely, a negative correlation would imply supralinearity.\n",
"\n",
"We also checked for any interaction that may be taking place between two different optical stimulation patterns. To quantify this, we measured distances on the grid map between all spots in all pairs of patterns, and compared it against the Vm change they caused at CA1 ([Figure 4—figure supplement 1e](#fig4s1))."
]
}
]
}